求教：“边坡安全系数求解”陈育民《FLAC/FLAC3D基础与工程实例》书上262页例子

haobozhu

陈育民《FLAC/FLAC3D基础与工程实例》书上262页例子“边坡安全系数求解”，我完全照着他的命令流输进去的，为什么得出的结果和书上的不一样啊，算出来安全系数Fos = 0 ，这肯定是错的，但是不知道哪错了，我是完全按照他的命令流的，是不是书上写的命令流有错误啊？希望各位高手指点一下，感激不尽！
书中的命令流如下：
new
;===================================
;建立网格模型
gen zone brick
&
p0 0 0 0 p1 2 0 0 p2 0 0.5 0 p3 0 0 3 size 3 1 3
gen zone brick
&
p0 2 0 0 p1 20 0 0 p2 2 0.5 0 p3 2 0 3 &
size 17 1 3 ratio 1.03 1 1
gen zone brick
&
p0 2 0 3 p1 20 0 3 p2 2 0.5 3 p3 12 0 13 &
p4 20 0.5 3 p5 12 0.5 13 p6 20 0 13 &
p7 20 0.5 13 size 17 1 17 ratio 1.03 1 1
;====================================
;设置边界条件
fix x y z range z -0.1 0.1
fix x range x 19.9 20.1
fix x range x -0.1 0.1
fix y
;====================================
;初始地应力的生成
model elas
prop density 2000 bulk 3e9 shear 1e9
set gravity 0 0 -10.0
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
;===================================
;安全系数求解
model mohr
prop density 2000 bulk 1.0e8 shear 3.0e7 &
coh 12380.0 tens 1.0e6
fric 20 dila 20
solve fos file slope3dfos1.sav associated


希望有人指点一下小弟，感激不尽，谢谢！
QQ：1176890449

kekylmk

我也是计算出来安全系数为0.

xiaolij

用restore命令查看结果

492854523

果然，用restore查看就可以了，这是为什么呢？

dingdinggou110

同解啊 郁闷