代码： a`q">T% q
;上面部分
,Z tj
gen zone wedge size 3 5 2 p0 0,0,35.1 p1 -3,0,31.1 p2 0,10,35.1 p3 3,0,33.7 p4 -3,10,31.1 p5 3,10,33.7 group 2 -7' >Rw
;裂隙上方方的部分 JT[|l-\zo
gen zone brick size 2 5 13 p0 0,0,9.1 p1 3,0,7.3 p2 0,10,9.1 p3 -3,0,31.1 p4 & mZ7.#R*}
3,10,7.3 p5 -3,10,31.1 p6 3,0,33.7 p7 3,10,33.7 group 3 lmj73OB3
attach face range group 2 3 {\;CGoN|
;夹层下方的三棱柱 Go w_a'
gen zone wedge size 4 5 10 p0 0,0,9 p1 0,0,0 p2 0,10,9 p3 15,0,0 p4 0,10,0 p5 15,10,0 group 1 2P!Pbl<
;裂隙右方的 s7(mNpo
gen zone brick size 13 5 8 p0 15.1,0,0.1 p1 75.1,0,0.1 p2 15.1,10,0.1 p3 3.1,0,7.3 & f/*Xw {s#
p4 75.1,10,0.1 p5 3.1,10,7.3 p6 3.1,0,33.7 p7 3.1,10,33.7 group 4 _D$| lk-
r m+|xvZ4
;产生软弱夹层 9N5 &N3
interface 1 face range plane norm (0.6,0,1) origin (7.5,0,4.5) dist 0.1 `cy_@Z5A
;产生顶部裂隙 +7^%fX;3pW
interface 2 face range x 3 z 7.2000001 34 P9G c)$6{p
ini z add=-0.1 range group 2 3 a&.8*|w3
ini z add=-0.1 range group 4 3 5B0L.R
ini x add=-0.1 range group 4 fk#SD "iJ
2o6KVQ
;设置边界条件 Rzj5B\+Rk(
fix x y z range z -0.1 0.1 (gvaYKvr
fix y range y -0.1 0.1 'r(}7>~fC
fix y range y 9.9 10.1 nI|jUD +y
rVt6tx
;模型材料性质 ojT TYR{
model mohr \uQ(-ji
ini density=2685 K> lA6i7?
prop bulk=4.21e10 shear=1.95e10 coh=7e6 fric=39.2 tens=11.9e6 range group 1 #)GL%{Oa
prop bulk=4.21e10 shear=1.95e10 coh=7e6 fric=39.2 tens=11.9e6 range group 2 X*]uLgbl
prop bulk=4.21e10 shear=1.95e10 coh=7e6 fric=39.2 tens=11.9e6 range group 3 u~bk~ 3.I
prop bulk=4.21e10 shear=1.95e10 coh=7e6 fric=39.2 tens=11.9e6 range group 4 LUHj3H
int 1 prop ks=4.21e11 kn=4.21e11 coh=2.2e6 fric=33.6 &P&M6v+
int 2 prop ks=4.21e11 kn=4.21e11 fric=39.2 (] (:0H
; 设置重力 [,8@oM#
set large qg'RD]a> R
set gravity 0 0 -10 CbS- Rz:
\P^WUWY
;初始化 < :<E~anH
ini xdis=0 ydis=0 zdis=0 ]6GdB3?UVM
@4sv(HyDY
;monitor histories zy
hist unbal b5No>U) /
pLDs eEr<
;cycling ht!o_0{~
solve L'H'E,
pl int yel sket 1bYc^(z0
;计算得到x位移图 / /3iai
请问我这样建模对吗 U;`N:~|p#
我用力学计算左边的块体应该是不动的（或位移很小），与模拟的不符 ?`u Y*+u
