陈育民《FLAC/FLAC3D基础与工程实例》书上262页例子“边坡安全系数求解”，我完全照着他的命令流输进去的，为什么得出的结果和书上的不一样啊，算出来安全系数Fos = 0 ，这肯定是错的，但是不知道哪错了，我是完全按照他的命令流的，是不是书上写的命令流有错误啊？希望各位高手指点一下，感激不尽！ Dn&D!B
书中的命令流如下： @=}YTt q
new ),G= s Oo
;=================================== &e 6CJ
;建立网格模型 #wL
gen zone brick & z 9F fU
p0 0 0 0 p1 2 0 0 p2 0 0.5 0 p3 0 0 3 size 3 1 3 ])xx<5Jt4
gen zone brick & 1%4sHSN
p0 2 0 0 p1 20 0 0 p2 2 0.5 0 p3 2 0 3 & 1N7Kv4,
size 17 1 3 ratio 1.03 1 1 I!e} )Y
gen zone brick & ]QzGE8jp*
p0 2 0 3 p1 20 0 3 p2 2 0.5 3 p3 12 0 13 & I)A`)5="5
p4 20 0.5 3 p5 12 0.5 13 p6 20 0 13 & a}%#* J)!
p7 20 0.5 13 size 17 1 17 ratio 1.03 1 1 n2)q}_d
;==================================== MrZh09y
;设置边界条件 ]E\n9X-{
fix x y z range z -0.1 0.1 t2,A@2DU 2
fix x range x 19.9 20.1 ; ;L[e]Z
fix x range x -0.1 0.1 UfN&v >8f
fix y ]gYz
4OT
;==================================== KMI_zhyB
;初始地应力的生成 ~0beuK&p
model elas ^*P%=>zO
prop density 2000 bulk 3e9 shear 1e9 &?mD$Eo
set gravity 0 0 -10.0 Ecd;<$tk
solve N"nd*?
ini xdisp 0 ydisp 0 zdisp 0 DxUKUE
ini xvel 0 yvel 0 zvel 0 j@YU|-\qh
;=================================== \,u_7y2 c
;安全系数求解 u%w`:v7Yo(
model mohr GB$;n?
prop density 2000 bulk 1.0e8 shear 3.0e7 & X!e[GJ
coh 12380.0 tens 1.0e6 fric 20 dila 20 KXx32 b,~
solve fos file slope3dfos1.sav associated `DV.+>O-1
q@[Qj Gj@
