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

我写了一个分层土边坡的FLAC命令流，可是计算不出结果，运行了十几个小时还没结果，现在把它发出来，请各位大侠帮我看看毛病，谢谢了：） PMfW;%I.  
Cz0FA]-g  
Ix-Mp   
new 8zp?WUb  
;*****几何模型***** $*ff]>#  
impgrid Flacmodel_haitang.Flac3D  >Uw:cq  
plot block group +<a\0FsD  
plot add axes red hzo> :U  
set gravity 0,0,-9.81 8K,X3a9  
group soil range group 1 any   aKDY_D  
group rock range group 2 any ubwM*P  
;施加边界约束条件 ev4[4T-(@  
fix x y z rang y -0.1 0.1 {
AD
d[V  
fix z range x 0 212 z -0.1 0.1 3`bQ0-D;  
fix x range x -0.1 0.1 z 0 130 +V4)><  
fix z range x 0 212 z 129.9 130.1 gJQ#j~'  
fix x range x 211.9 212.1 z 0 130 pF{jIXu  
;通过弹性求解生成初始应力场 tF:AnNp=  
model elas (BEe^]f  
pro bulk 1e10 she 1e10 range group soil O] @E8<?^  
pro bulk 1e10 she 1e10 range group rock fz(YP=@ZnP  
solve fo 10 XQo
\27Fo  
;位移和速度归零 vGkemJ^/  
ini xdisp 0 ydisp 0 zdisp 0 .PB!1C.}@  
ini xvel 0 yvel 0 zvel 0 2_q/
<8t  
;定义材料特性 V$?6%\M^*  
model mohr 'Aq^z%|  
pro bulk 9.1e8 she 6.5e8 co 2.10e4 fric 28 ten 1e6  & MgeC-XQM  
dil 10 range group soil MgXZN{  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 ten 1e6  & W_W!v&@E=  
dil 15 range group rock (/A.,8Ad  
hist n=5 .KrLvic  
hist unbal danPy2  
plo hist 1 D0Dz@25-  
set mech ratio 1.0e-6 BJjic%V  
solve