[讨论]程序检查 [复制链接]

我写了一个分层土边坡的FLAC命令流，可是计算不出结果，运行了十几个小时还没结果，现在把它发出来，请各位大侠帮我看看毛病，谢谢了：） m7^a4  
[,|KVc=&H  
+6uun  
new B,|M  
;*****几何模型***** r/:s2oQ  
impgrid Flacmodel_haitang.Flac3D Yca9G?^\v  
plot block group u\&oiwSIP  
plot add axes red 2%1g%  
set gravity 0,0,-9.81 n4(w?,w}  
group soil range group 1 any   {HvR24
#  
group rock range group 2 any Vg6?a  
;施加边界约束条件 Af ^6  
fix x y z rang y -0.1 0.1 #=Q/<r.~G  
fix z range x 0 212 z -0.1 0.1 x-CYG?-x  
fix x range x -0.1 0.1 z 0 130 RYS]b[-xZz  
fix z range x 0 212 z 129.9 130.1 8p>%}LX/  
fix x range x 211.9 212.1 z 0 130 JB''Ujyi  
;通过弹性求解生成初始应力场 htlsU*x  
model elas -:cS}I  
pro bulk 1e10 she 1e10 range group soil !bT0kP$3}  
pro bulk 1e10 she 1e10 range group rock fC]+C(*d  
solve fo 10 v?n`kw  
;位移和速度归零 }B^s!y&b  
ini xdisp 0 ydisp 0 zdisp 0 cV)fe`Gg  
ini xvel 0 yvel 0 zvel 0 ZEUd?"gaR  
;定义材料特性 ,t61IU3"  
model mohr :a#]"z0  
pro bulk 9.1e8 she 6.5e8 co 2.10e4 fric 28 ten 1e6  & ]Fl+^aLS  
dil 10 range group soil %!p14c*J H  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 ten 1e6  & VH+^G)^)W  
dil 15 range group rock $:/y5zi  
hist n=5 [WX+/pm7>  
hist unbal 6SlE>b9tA  
plo hist 1 X1#D}  
set mech ratio 1.0e-6 : gv[X  
solve