陈育民《FLAC/FLAC3D基础与工程实例》书上262页例子“边坡安全系数求解”，我完全照着他的命令流输进去的，为什么得出的结果和书上的不一样啊，算出来安全系数Fos = 0 ，这肯定是错的，但是不知道哪错了，我是完全按照他的命令流的，是不是书上写的命令流有错误啊？希望各位高手指点一下，感激不尽！ 1c$c e+n~
书中的命令流如下： 8`q"] BQN
new / ? Bu^KX
;=================================== t>"|~T$9
;建立网格模型 h9 + 7 6
gen zone brick & 8ZDWaq8^2N
p0 0 0 0 p1 2 0 0 p2 0 0.5 0 p3 0 0 3 size 3 1 3 EL9JM}%0v
gen zone brick & ,*;g+[Bhpl
p0 2 0 0 p1 20 0 0 p2 2 0.5 0 p3 2 0 3 & tSaD =# v
size 17 1 3 ratio 1.03 1 1 l q)[
gen zone brick & l, j0n0h.
p0 2 0 3 p1 20 0 3 p2 2 0.5 3 p3 12 0 13 & \}7xgQ>oV
p4 20 0.5 3 p5 12 0.5 13 p6 20 0 13 & p1,.f&(f
p7 20 0.5 13 size 17 1 17 ratio 1.03 1 1 g"D:zK)
;==================================== DXfQy6k'
;设置边界条件 v0+$d\mP4<
fix x y z range z -0.1 0.1 05;J7T<
fix x range x 19.9 20.1 f7 zGz
fix x range x -0.1 0.1 2O0< /^Z%E
fix y *BQy$dfE
;==================================== 4I;$a;R!
;初始地应力的生成 "%+9p6/
model elas ?;p45y~n%
prop density 2000 bulk 3e9 shear 1e9 k
i{8f
set gravity 0 0 -10.0 n si&r
solve S!7|vb*ko
ini xdisp 0 ydisp 0 zdisp 0 p 4l B#
ini xvel 0 yvel 0 zvel 0 C0'_bTfB
;=================================== us^J!
s7
;安全系数求解 =OooTZb:x-
model mohr R,W
w/D
prop density 2000 bulk 1.0e8 shear 3.0e7 & E)( Rhvij
coh 12380.0 tens 1.0e6 fric 20 dila 20 I yt.`z
solve fos file slope3dfos1.sav associated x]|-2t
);V.le}% (
