陈育民《FLAC/FLAC3D基础与工程实例》书上262页例子“边坡安全系数求解”，我完全照着他的命令流输进去的，为什么得出的结果和书上的不一样啊，算出来安全系数Fos = 0 ，这肯定是错的，但是不知道哪错了，我是完全按照他的命令流的，是不是书上写的命令流有错误啊？希望各位高手指点一下，感激不尽！ f:S}h-AL&
书中的命令流如下： Rw/Ciw2@?
new m!;G/s*
;=================================== N_0pO<<cs
;建立网格模型 %
^e@`0L
gen zone brick & :fhB*SYK
p0 0 0 0 p1 2 0 0 p2 0 0.5 0 p3 0 0 3 size 3 1 3 3<+z46`?
gen zone brick & CNq[4T'~A
p0 2 0 0 p1 20 0 0 p2 2 0.5 0 p3 2 0 3 & SGQD ro=l
size 17 1 3 ratio 1.03 1 1 S3QaYq"v
gen zone brick & Jlz9E|*qV
p0 2 0 3 p1 20 0 3 p2 2 0.5 3 p3 12 0 13 & 1}`2\3,
p4 20 0.5 3 p5 12 0.5 13 p6 20 0 13 & !h?=Wv
==]
p7 20 0.5 13 size 17 1 17 ratio 1.03 1 1 ;e2Ij
;==================================== sLNNcj(Cy>
;设置边界条件 (,shiK[5f
fix x y z range z -0.1 0.1 ZxI]I1)
fix x range x 19.9 20.1 %Or2iuO%-,
fix x range x -0.1 0.1 s88y{o
fix y v3~FR,Kl
;==================================== Zct!/u9 Q
;初始地应力的生成 \PzN XQ$
model elas z1#oW f{*
prop density 2000 bulk 3e9 shear 1e9 ~C0 Pu.{o
set gravity 0 0 -10.0
K,6OGsh
solve L -YNz0A
ini xdisp 0 ydisp 0 zdisp 0 C]M7GHe1q
ini xvel 0 yvel 0 zvel 0 +eC3?B8rN
;=================================== &"xQ~05
;安全系数求解 E3bwy K!s
model mohr l`* ( f9Q
prop density 2000 bulk 1.0e8 shear 3.0e7 & EPv%LX_j
coh 12380.0 tens 1.0e6 fric 20 dila 20 4Q$!c{Y
r
solve fos file slope3dfos1.sav associated b 1 H7
<;aJ#qT
