我根据陈育民那本书上的例子改了一下，但是最后算出来的FOS值一直是0，跪求各位大神指导
pl bl gr
attach face
gen merge 10e-5
;对x y z三个方向的速度进行约束
fix z ran z -209.9 -210.1
fix x ran x -0.1 0.1
fix x ran x 629.9 630.1
fix y ran y -0.1 0.1
fix y ran y 379.9 380.1
;初始地应力生成
model elas
prop bulk 2.456e7 shear 1.069e7 dens 2000 ran group 5
prop bulk 8.33e8 shear 3.85e8 dens 2200 ran group 4
prop bulk 3.33e7 shear 1.11e7 dens 2050 ran group 3
prop bulk 2.78e7 shear 9.26e6 dens 2030 ran group 2
prop bulk 1.852e7 shear 7.58e6 dens 1950 ran group 1
set grav 0 0 -9.80
so
ini xdis 0 ydis 0 zdis 0 ;将节点位移清零
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零
;安全系数求解
model mohr
prop bulk 2.456e7 shear 1.069e7 coh 6.35e4 friction 27.5 dila 27.5 dens 2000 ran group 5
prop bulk 8.33e8 shear 3.85e8 coh 2e5 friction 35 dila 35 dens 2200 ran group 4
prop bulk 3.33e7 shear 1.11e7 coh 9e4 friction 20 dila 20 dens 2050 ran group 3
prop bulk 2.78e7 shear 9.26e6 coh 5.5e4 friction 26 dila 26 dens 2030 ran group 2
prop bulk 1.852e7 shear 7.58e6 coh 3.5e4 friction 23 dila 23 dens 1950 ran group 1
solve fos file slope3dfos1.sav associated ;采用关联流动法则求解安全系数(摩擦角等于剪胀角)
pause
pl fos cont ssi outline on vel red ;绘制安全系数、剪切应变增量及速度矢量图
pause
def _getfos ;定义获取高精度安全系数的函数
fos1=fos_f ;将当前安全系数值赋给fos_f
end ;函数定义结束
_getfos ;执行已定义的函数
print fos1 ;在屏幕上显示当前安全系数值
