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

我写了一个分层土边坡的FLAC命令流，可是计算不出结果，运行了十几个小时还没结果，现在把它发出来，请各位大侠帮我看看毛病，谢谢了：） Z_T~2t  
^vOEG;TR<-  
B5;94YIN  
new ^[2siG  
;*****几何模型***** eYv+tjIF  
impgrid Flacmodel_haitang.Flac3D ]Rmu+N|  
plot block group oL9ELtb]s  
plot add axes red F5w=tK  
set gravity 0,0,-9.81 b}hQU~,E  
group soil range group 1 any   =[gFa
B_H  
group rock range group 2 any 2D3mTpw  
;施加边界约束条件 H;FzWcm  
fix x y z rang y -0.1 0.1 Ka"1gbJ|  
fix z range x 0 212 z -0.1 0.1 P1`YbLER5  
fix x range x -0.1 0.1 z 0 130 9HlM0qE5b  
fix z range x 0 212 z 129.9 130.1 QX.U:p5C  
fix x range x 211.9 212.1 z 0 130 M IUB]  
;通过弹性求解生成初始应力场 hm1.UE  
model elas ;;EFiaA  
pro bulk 1e10 she 1e10 range group soil ;*20b@  
pro bulk 1e10 she 1e10 range group rock 8PG&/"K  
solve fo 10 ~AF' 6"A  
;位移和速度归零 FGpV ]p  
ini xdisp 0 ydisp 0 zdisp 0 grJ(z)c  
ini xvel 0 yvel 0 zvel 0 J]Q-#g'Z  
;定义材料特性 w&&)v~Y_  
model mohr =Q!V6+}nY^  
pro bulk 9.1e8 she 6.5e8 co 2.10e4 fric 28 ten 1e6  & .O{_^~w_q  
dil 10 range group soil Jp~[Dm  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 ten 1e6  & )LBbA  
dil 15 range group rock ,^@/I:  
hist n=5 L|A1bxt  
hist unbal XKT[8o<L  
plo hist 1 b1Ba}  
set mech ratio 1.0e-6 \@_?mL@=  
solve