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

我写了一个分层土边坡的FLAC命令流，可是计算不出结果，运行了十几个小时还没结果，现在把它发出来，请各位大侠帮我看看毛病，谢谢了：） /xK5%cE>B  
##%R|P3  
.Yl*kG6r  
new x]({Po4  
;*****几何模型***** /eM_:H5  
impgrid Flacmodel_haitang.Flac3D %-zAV*>  
plot block group 9Y9pKTU  
plot add axes red #3:;&@#  
set gravity 0,0,-9.81 :+1S+w  
group soil range group 1 any   9!XW):  
group rock range group 2 any qO5.NIs  
;施加边界约束条件 <K8$00l
m  
fix x y z rang y -0.1 0.1 !PCw-&  
fix z range x 0 212 z -0.1 0.1 4o:
  
fix x range x -0.1 0.1 z 0 130 "@xF(fyg  
fix z range x 0 212 z 129.9 130.1 fk(l.A$  
fix x range x 211.9 212.1 z 0 130 S#M<d~rK  
;通过弹性求解生成初始应力场 3@k;"pFa<  
model elas @!92Ok  
pro bulk 1e10 she 1e10 range group soil R~-r8dWcw  
pro bulk 1e10 she 1e10 range group rock YF{KSGq  
solve fo 10 %cM2;a=2  
;位移和速度归零 # ;,b4O7@  
ini xdisp 0 ydisp 0 zdisp 0 )^&,[Q=i  
ini xvel 0 yvel 0 zvel 0 `r; .  
;定义材料特性 B_u+$Odo  
model mohr gnS0$kCJ:  
pro bulk 9.1e8 she 6.5e8 co 2.10e4 fric 28 ten 1e6  & Dh*>361y-  
dil 10 range group soil d2w;d&2S  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 ten 1e6  & N=) E$h  
dil 15 range group rock ~jd:3ip+!  
hist n=5 .Gq)@{o>  
hist unbal oa:GGW4Q  
plo hist 1 PVAs# ~  
set mech ratio 1.0e-6 %eoO3"//  
solve