[求助]高手帮帮忙 啊 ！！谢谢  啊！！ [复制链接]

代码： ykH?;Xu  
;上面部分 |G~LJsXW!v  
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   k, &*d4  
;裂隙上方方的部分 p [4/Nq,c  
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 & 3*"
$E_%  
    3,10,7.3 p5 -3,10,31.1 p6 3,0,33.7 p7 3,10,33.7  group 3 hW
$B;  
attach face range group 2  3 ^\Nsx)Y;  
;夹层下方的三棱柱 V~tq _  
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 Cb7f-Eag  
;裂隙右方的 w?d~c*4+  
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 & tI|?k(D  
        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 QM=M<~<Voh  
5e~\o}
]  
;产生软弱夹层 8Gzc3  
interface 1 face range plane norm (0.6,0,1)  origin (7.5,0,4.5) dist 0.1 _(&^M[O  
;产生顶部裂隙 hn#i,XnY  
interface 2 face range x 3 z 7.2000001  34  $J>GCY  
ini z add=-0.1 range group 2 3 \?{nP6=  
ini z add=-0.1 range group 4 dkt'~  
ini x add=-0.1 range group 4 vcy}ZqWBO  
O6NgI2[O  
;设置边界条件 '3O@Nxof4  
fix x y z  range z -0.1 0.1 8rAOs\ys  
fix y range y -0.1 0.1 Mp^%.m  
fix y range y 9.9 10.1 0
'wB':v  
xAw$bJj~s  
;模型材料性质 qvy
~b  
model mohr XZLo*C!MG  
ini density=2685 Ci0:-IS  
prop   bulk=4.21e10  shear=1.95e10  coh=7e6  fric=39.2 tens=11.9e6 range group 1 @tWyc%t  
prop   bulk=4.21e10  shear=1.95e10  coh=7e6  fric=39.2 tens=11.9e6 range group 2 ~0,Utqy  
prop   bulk=4.21e10  shear=1.95e10  coh=7e6  fric=39.2 tens=11.9e6 range group 3 (|I:
d!>:U  
prop   bulk=4.21e10  shear=1.95e10  coh=7e6  fric=39.2 tens=11.9e6  range group 4 s9>f5u?d
K  
int 1 prop  ks=4.21e11   kn=4.21e11  coh=2.2e6  fric=33.6 1-r#v  
int 2 prop  ks=4.21e11   kn=4.21e11   fric=39.2 1T a48  
; 设置重力 L!Iu\_{q  
set large `9n%Dy<  
set gravity 0 0 -10 SH2|xn  
Q3vC^}Dmr  
;初始化 Ol1[o  
ini xdis=0 ydis=0 zdis=0 4d#w}  
U8KB@E  
;monitor histories 0tv"tA;  
hist unbal +!6aB|-  
ce{(5IC  
;cycling "rOe J~4 X  
solve M>VT$!Lx  
pl int yel sket `=g9Rg/<  
;计算得到x位移图 o7)<pfif  
请问我这样建模对吗 Lg0Vn&k  
我用力学计算左边的块体应该是不动的（或位移很小），与模拟的不符 w$WN` =  

看不明白你写的是是什么东东