代码： D&HV6#
;上面部分 PRNoqi3sY
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 1-RIN}CSd
;裂隙上方方的部分 A(]H{>PMy
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 & E( us'9c
3,10,7.3 p5 -3,10,31.1 p6 3,0,33.7 p7 3,10,33.7 group 3 `9k 0Gd
attach face range group 2 3 @
49nJi
;夹层下方的三棱柱 ie-vqLc
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 VLBE'3Qg 1
;裂隙右方的 bQ| V!mrN}
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 & yV) la@c
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 ,lCgQ0}<
sg E-`#
;产生软弱夹层 xkOpa,=FI
interface 1 face range plane norm (0.6,0,1) origin (7.5,0,4.5) dist 0.1 s+:= I
e
;产生顶部裂隙 8w({\=
interface 2 face range x 3 z 7.2000001 34 fO#vF.k%
ini z add=-0.1 range group 2 3 ;gC |
ini z add=-0.1 range group 4 =@F&o4) r
ini x add=-0.1 range group 4 fwzb!"!.@
r-,e;o>9
;设置边界条件 ( H+'X}1
fix x y z range z -0.1 0.1 gWY "w!f
fix y range y -0.1 0.1 Zo>]rKeV
fix y range y 9.9 10.1 L@JOGCYy
A. UUW
;模型材料性质 W2uOR{
'?
model mohr |W<+U
ini density=2685 p&VU0[LIC0
prop bulk=4.21e10 shear=1.95e10 coh=7e6 fric=39.2 tens=11.9e6 range group 1 :$MG*/Q
prop bulk=4.21e10 shear=1.95e10 coh=7e6 fric=39.2 tens=11.9e6 range group 2 qY 4#V k
prop bulk=4.21e10 shear=1.95e10 coh=7e6 fric=39.2 tens=11.9e6 range group 3 *,Bzc Z
prop bulk=4.21e10 shear=1.95e10 coh=7e6 fric=39.2 tens=11.9e6 range group 4 $=?@*p
int 1 prop ks=4.21e11 kn=4.21e11 coh=2.2e6 fric=33.6 QRLt9L
int 2 prop ks=4.21e11 kn=4.21e11 fric=39.2 1:Si,d,wh
; 设置重力 OT'[:|x ;
set large _G1gtu]
set gravity 0 0 -10 !x'/9^i~v
bI|2@H V2
;初始化 Z,iH y3`
ini xdis=0 ydis=0 zdis=0 dL|+d: v
u1xSp<59C
;monitor histories jY_T/233d
hist unbal #97h6m?
!%dN<%Ah
;cycling Fs[aa#v4B
solve | d6T/Uxo
pl int yel sket <3,<\ub
;计算得到x位移图 )E+'*e{cK
请问我这样建模对吗 b,8{ X<
我用力学计算左边的块体应该是不动的（或位移很小），与模拟的不符 %'0T Xr$
