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

我写了一个分层土边坡的FLAC命令流，可是计算不出结果，运行了十几个小时还没结果，现在把它发出来，请各位大侠帮我看看毛病，谢谢了：） =x9SvIm/tH  
6S#Y$2 P  
ZLsfF =/G  
new .\r=1HZ3  
;*****几何模型***** i$~2pr  
impgrid Flacmodel_haitang.Flac3D 6jv_j[[  
plot block group 4uQ\JD(*Eu  
plot add axes red q^:>sfd  
set gravity 0,0,-9.81 \8Mn[G9TL  
group soil range group 1 any   yc4?'k!  
group rock range group 2 any XCCh*qym  
;施加边界约束条件 %4F Q~  
fix x y z rang y -0.1 0.1 T&M*sydA  
fix z range x 0 212 z -0.1 0.1 6OBe^/ZRt  
fix x range x -0.1 0.1 z 0 130 _a;E>   
fix z range x 0 212 z 129.9 130.1 A1b</2  
fix x range x 211.9 212.1 z 0 130 Du
ESLMhz  
;通过弹性求解生成初始应力场 UKYQ @m  
model elas ~}uv4;0l]  
pro bulk 1e10 she 1e10 range group soil v"nN[_T  
pro bulk 1e10 she 1e10 range group rock do {E39  
solve fo 10 1&zvf4  
;位移和速度归零  ]/l"  
ini xdisp 0 ydisp 0 zdisp 0 H
GuU6@~hu  
ini xvel 0 yvel 0 zvel 0 j_` [Z  
;定义材料特性 )?5027^  
model mohr {FYWQ!L  
pro bulk 9.1e8 she 6.5e8 co 2.10e4 fric 28 ten 1e6  & AYhWeI+  
dil 10 range group soil
"I[a]T}/  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 ten 1e6  & nkI+"$Rz0  
dil 15 range group rock -s?f<f{  
hist n=5 8*nv+  
hist unbal ^!Y]l  
plo hist 1 `pMI@"m  
set mech ratio 1.0e-6 w %4SNR  
solve