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

我写了一个分层土边坡的FLAC命令流，可是计算不出结果，运行了十几个小时还没结果，现在把它发出来，请各位大侠帮我看看毛病，谢谢了：） J[}j8x?r  
]HWeVhG  
7oC8ID  
new odL*_<Z  
;*****几何模型***** &am<_Tn*3  
impgrid Flacmodel_haitang.Flac3D gR/?MJ(v  
plot block group SOPair <r  
plot add axes red )W_akUL  
set gravity 0,0,-9.81 y=Eb->a){  
group soil range group 1 any   #y&
5pP:@  
group rock range group 2 any ByE@4+9  
;施加边界约束条件 ShQ!'[J  
fix x y z rang y -0.1 0.1 TQ![  
fix z range x 0 212 z -0.1 0.1 e6*,MnqBh  
fix x range x -0.1 0.1 z 0 130 JV,h1/a("  
fix z range x 0 212 z 129.9 130.1 papMC"<g$  
fix x range x 211.9 212.1 z 0 130 r7R39#  
;通过弹性求解生成初始应力场 n+ H2cl }  
model elas \hBzQ%0  
pro bulk 1e10 she 1e10 range group soil BJ_"FG  
pro bulk 1e10 she 1e10 range group rock X+{brvM<  
solve fo 10 y ~-v0/  
;位移和速度归零 N~b0b;e  
ini xdisp 0 ydisp 0 zdisp 0 Ca5LLG  
ini xvel 0 yvel 0 zvel 0 X6}W]  
;定义材料特性 j0LZ )V  
model mohr XlXt,  
pro bulk 9.1e8 she 6.5e8 co 2.10e4 fric 28 ten 1e6  & F7<u1Rx]  
dil 10 range group soil AP`1hz4].-  
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 ten 1e6  & %INkuNa8\  
dil 15 range group rock du0o4~-  
hist n=5 cD9U^SOS  
hist unbal 60n>FQ<  
plo hist 1 X!{K`~DRX  
set mech ratio 1.0e-6 Y9-F\t=~  
solve