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

用FLAC3D程序写了个边坡稳定性的程序，得到的结果，不会分析啊，是哪写错了吗？请高手指点下………… C(qqGK{  
new Rs_bM@  
gen zone brick p0=(0,0,0) p1=(3,0,0) p2=(0,1,0) p3=(0,0,3) size 2 1 3 
Gqvj  
gen zone brick p0=(3,0,0) p1=(35,0,0) p2=(3,1,0) p3=(3,0,3) size 17 1 3 ()i8 Qepo}  
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 !)FM/Xj,o  
fix x y z range z -0.1 0.1 BuII|j  
fix x range x 34.9 35.1 8pp^ w  
fix x range x -0.1 0.1
Nz %{T  
fix y Z<
T%:F  
model elastic ~ x- R78'  
property density 2000 bulk 3e9 shear 1e9 Ke@zS9  
set gravity 0 0 -10.0 O RAKg.49  
solve #Y6'Q8gf  
initial xdisplacement 0 ydisplacement 0 zdisplacement 0 of!Bz  
initial xvelocity 0 yvelocity 0 zvelocity 0 WB<_AIt+  
model mohr SO^:6GuJ  
property density 2050 bulk 1e8 shear 3e7 cohesion 12380.0 tension 0.98e6 friction 15 dila 15 wyvrNru<l4  
solve F$T@OT6  
plot contour ssi outline on vel red -*t4(wT|j  
:yO.Te F  

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

". #=_/op  
你的命令流没什么错误，不知道你遇到什么问题了 *<Yn  
这是你的结果图 X;25G  

描述:计算结果

图片:无标题.jpg

'i#m%D`dt  
oVk*G  

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

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