一个应变软化的边坡例子

chengjia02

看了一篇论文的例子，想动手做一做，出错了，密度自己加的，想高手指导一下，谢谢
n
gen zone b p0 0 0 0 p1 9 0 0 p2 0 1 0 p3 0 0 6 size 45 5 30
pl grid
mo ss
prop bulk 1e8 she 3e7
prop coh 1.2e4 fric 20 ten 1e4
prop ctab 1 ftab 2 ttab 3
table 1 0 1.2e4 0.02 0
table 2 0 20 0.02 15
table 3 0 1e4 0.02 0
set grav 0 0 -9.81
water den 1000
group water_below range z(0,2)
group water_above range z(2,6)
gen separate water_below
interface 1 wrap water_below  water_above
plot add interface
water table face 0,0,2 0,1,2 6,1,2
ini dens 1800 range z(0,2)
ini dens 2000 range z(0,2)
fix x range x -.1 .1
fix x range x 8.9 9.1
fix y
fix z range z -.1 .1
model null range x 6 9 y 0 1 z 2 6
pl sur
set large
his unbal
pl con ssr
step 3000

[ 本帖最后由 chengjia02 于 2008-11-13 15:54 编辑 ]

chengjia02

截下来的图：

chengjia02

效果：剪应变云图

chengjia02

米人指教么，等待……

chengjia02

谁来帮帮我啊