flac3d边坡稳定性程序，如何分析 [复制链接]

用FLAC3D程序写了个边坡稳定性的程序，得到的结果，不会分析啊，是哪写错了吗？请高手指点下………… tJ^p}yxO  
new O52/fGt  
gen zone brick p0=(0,0,0) p1=(3,0,0) p2=(0,1,0) p3=(0,0,3) size 2 1 3 x"b'Pmw  
gen zone brick p0=(3,0,0) p1=(35,0,0) p2=(3,1,0) p3=(3,0,3) size 17 1 3 :AzT=^S  
gen zone brick p0=(3,0,3) p1=(35,0,3) p2=(3,1,3) p3=(18,0,13) p4=(35,1,3) p5=(18,1,13) p6=(35,0,13) p7=(35,1,13) size 17 1 10 P 2WAnm  
fix x y z range z -0.1 0.1 l!tR
<$|  
fix x range x 34.9 35.1 IbI0".o  
fix x range x -0.1 0.1 sy
cAAmH<  
fix y &&
&-P\3  
model elastic 4,)9@-|0R  
property density 2000 bulk 3e9 shear 1e9 u9!  ?  
set gravity 0 0 -10.0 L_Ok?9$  
solve 4'O,xC  
initial xdisplacement 0 ydisplacement 0 zdisplacement 0 ?9~^QRLT  
initial xvelocity 0 yvelocity 0 zvelocity 0 ?\o~P  
model mohr pkn^K+<n,  
property density 2050 bulk 1e8 shear 3e7 cohesion 12380.0 tension 0.98e6 friction 15 dila 15 HA,o2jZ?In  
solve iXMJ1\!q\|  
plot contour ssi outline on vel red L I<S  
9+@h2"|N4*  

你的命令流没什么错误，不知道你遇到什么问题了

P(shbi@  
你的命令流没什么错误，不知道你遇到什么问题了 OijuOLt  
这是你的结果图 NxHUOPAJc  

描述:计算结果

图片:无标题.jpg

aC,adNub  
Z}S[
fN8  

第二个求解命令solve应改成solve fos 进行安全系数求解

这个solve fos命令，运算时间超级长，模型复杂了就算的很艰难了。用书上的强度折减法，在设置体积模量和剪切模量的时候，土的数量级取10^6~7次方，很难收敛，取10^8~9计算容易收敛，但不知会不会出现问题