陈育民《FLAC/FLAC3D基础与工程实例》书上262页例子“边坡安全系数求解”,我完全照着他的命令流输进去的,为什么得出的结果和书上的不一样啊,算出来安全系数Fos = 0 ,这肯定是错的,但是不知道哪错了,我是完全按照他的命令流的,是不是书上写的命令流有错误啊?希望各位高手指点一下,感激不尽! Xl@cHO=i
书中的命令流如下: 3\T2?w9u(
new ?g!V!VS2
;=================================== 52upoU>}2
;建立网格模型 iH^z:%dP
gen zone brick & [ sd;`xk
p0 0 0 0 p1 2 0 0 p2 0 0.5 0 p3 0 0 3 size 3 1 3 ~=R SKyzt
gen zone brick & qj cp65^
p0 2 0 0 p1 20 0 0 p2 2 0.5 0 p3 2 0 3 & >
iE!m
size 17 1 3 ratio 1.03 1 1 8kP3+
gen zone brick & }I`a`0/
p0 2 0 3 p1 20 0 3 p2 2 0.5 3 p3 12 0 13 & &rkEK4
p4 20 0.5 3 p5 12 0.5 13 p6 20 0 13 & bc
, p}
p7 20 0.5 13 size 17 1 17 ratio 1.03 1 1 rJ KZ)N{
;==================================== D&HV6#
;设置边界条件 5NJ4
fix x y z range z -0.1 0.1 hHqh{:q{v
fix x range x 19.9 20.1 hzk6rYg1
fix x range x -0.1 0.1 Kx_h1{
fix y wP"dZagpj
;==================================== v]B
L[/4
;初始地应力的生成 Qr
Wj>uR
model elas ;S xFp
prop density 2000 bulk 3e9 shear 1e9 "UwH\T4I
set gravity 0 0 -10.0 s,a}?W
solve czlFr|O;
ini xdisp 0 ydisp 0 zdisp 0 ^5r9 5
ini xvel 0 yvel 0 zvel 0 E^U0f/5
m
;=================================== sgE-`#
;安全系数求解 sB69R:U;
model mohr ?xK,mbFgl
prop density 2000 bulk 1.0e8 shear 3.0e7 & h*u`X>!!
coh 12380.0 tens 1.0e6 fric 20 dila 20 !mXxAo
solve fos file slope3dfos1.sav associated iAa;6mH
}w4QP+ x