一个简单边坡的模拟，但感觉有点问题，还请大家指点一下，哪儿有问题？谢谢
以下是命令：new
;--------------------建立网格-----------------------
gen zone brick p0 274.989 0 56 p1 118.6168 0 56 p2 274.989 -1 56 p3 274.989 0 75 &
p4 118.6168 -1 56 p5 274.989 -1 75 p6 169.7607 0 75 p7 169.7607 -1 75 &
size 50,1,19 &
ratio 0.96,1,1.05
gen zone brick p0 274.989 0 50 p1 118.6168 0 50 p2 274.989 -1 50 p3 274.989 0 56 &
size 50,1,10 &
ratio 0.96,1,1
gen zone brick p0 274.989 0 0 p1 118.6168 0 0 p2 274.989 -1 0 p3 274.989 0 50 &
size 50,1,30 &
ratio 0.96,1,0.94
gen zone wedge p0 118.6168 0 50 p1 117.5333 0 50 p2 118.6168 -1 50 p3 118.6168 0 56 &
p4 117.5333 -1 50 p5 118.6168 -1 56 &
size 2,1,10 &
ratio 1,1,1
gen zone brick p0 118.6168 0 0 p1 117.5333 0 0 p2 118.6168 -1 0 p3 118.6168 0 50 &
size 1,1,30 &
ratio 1,1,0.94
gen zone wedge p0 117.5333 0 50 p1 116.4076 0 50 p2 117.5333 -1 50 p3 118.6168 0 56 &
p4 116.4076 -1 50 p5 118.6168 -1 56 &
size 2,1,10 &
ratio 1,1,1
gen zone brick p0 117.5333 0 0 p1 116.4076 0 0 p2 117.5333 -1 0 p3 117.5333 0 50 &
size 1,1,30 &
ratio 1,1,0.94
gen zone wedge p0 116.4076 0 50 p1 115.1069 0 50 p2 116.4076 -1 50 p3 118.6168 0 56 &
p4 115.1069 -1 50 p5 118.6168 -1 56 &
size 2,1,10 &
ratio 1,1,1
gen zone brick p0 116.4076 0 0 p1 115.1069 0 0 p2 116.4076 -1 0 p3 116.4076 0 50 &
size 1,1,30 &
ratio 1,1,0.94
gen zone wedge p0 115.1069 0 50 p1 113.5734 0 50 p2 115.1069 -1 50 p3 118.6168 0 56 &
p4 113.5734 -1 50 p5 118.6168 -1 56 &
size 2,1,10 &
ratio 1,1,1
gen zone brick p0 115.1069 0 0 p1 113.5734 0 0 p2 115.1069 -1 0 p3 115.1069 0 50 &
size 1,1,30 &
ratio 1,1,0.94
gen zone wedge p0 113.5734 0 50 p1 111.353 0 50 p2 113.5734 -1 50 p3 118.6168 0 56 &
p4 111.353 -1 50 p5 118.6168 -1 56 &
size 2,1,10 &
ratio 1,1,1
gen zone brick p0 113.5734 0 0 p1 111.353 0 0 p2 113.5734 -1 0 p3 113.5734 0 50 &
size 1,1,30 &
ratio 1,1,0.94
gen zone wedge p0 111.353 0 50 p1 108.1552 0 50 p2 111.353 -1 50 p3 118.6168 0 56 &
p4 108.1552 -1 50 p5 118.6168 -1 56 &
size 2,1,10 &
ratio 1,1,1
gen zone brick p0 111.353 0 0 p1 108.1552 0 0 p2 111.353 -1 0 p3 111.353 0 50 &
size 1,1,30 &
ratio 1,1,0.94
gen zone wedge p0 108.1552 0 50 p1 102.4661 0 50 p2 108.1552 -1 50 p3 118.6168 0 56 &
p4 102.4661 -1 50 p5 118.6168 -1 56 &
size 2,1,10 &
ratio 1,1,1
gen zone brick p0 108.1552 0 0 p1 102.4661 0 0 p2 108.1552 -1 0 p3 108.1552 0 50 &
size 2,1,30 &
ratio 1,1,0.94
gen zone brick p0 102.4661 0 0 p1 0 0 0 p2 102.4661 -1 0 p3 102.4661 0 50 &
size 20,1,30 &
ratio 1.05,1,0.94
attach face range z 49 51
;---------------------------------------------------
plot create slope
plot add axes
plot add sur
plot show
;--------------------给定模型参数-----------------------
model elas
prop bulk 1e7 shear 1e7
ini dens 1270 range z 68 75;PA-01 (上部也用pa-01)
ini dens 1560 range z 66 68;PA-02
ini dens 1620 range z 64 66;PA-03
ini dens 1790 range z 62 64;PA-04
ini dens 1720 range z 60 62;PA-05
ini dens 1810 range z 58 60;PA-06
ini dens 1860 range z 56 58;PA-07
ini dens 1640 range z 54 56;PA-08
ini dens 1480 range z 52 54;PA-09
ini dens 1470 range z 052;PA-10 (下部也用pa-10)
fix x range x -0.1 0.1
fix x range x 274.889 275.089
fix z range z -0.1 0.1
fix y
set grav 0,0,-9.81
set mech force=1e-3
solve
ini xdis 0 ydis 0
ini xvel 0 yvel 0
ini state 0
model moh
prop bulk 11.08e6shear 5.41e6fric 14.8coh 32.1e3tens 0range z 68 75;PA-01 (上部也用pa-01)
prop bulk 10.11e6shear 4.94e6fric 13.4coh 48.4e3tens 0range z 66 68;PA-02
prop bulk 5.81e6 shear 2.84e6fric 13.2coh 48.3e3tens 0range z 64 66;PA-03
prop bulk 9.27e6 shear 4.53e6fric 16.1coh 54.3e3tens 0range z 62 64;PA-04
prop bulk 10.54e6shear 5.15e6fric 15.6coh 55.0e3tens 0range z 60 62;PA-05
prop bulk 11.39e6shear 5.56e6fric 14.6coh 59.2e3tens 0range z 58 60;PA-06
prop bulk 10.05e6shear 4.91e6fric 16.1coh 55.5e3tens 0range z 56 58;PA-07
prop bulk 4.06e6 shear 1.98e6fric 14.1coh 45.6e3tens 0range z 54 56;PA-08
prop bulk 2.91e6 shear 1.42e6fric 13.7coh 40.6e3tens 0range z 52 54;PA-09
prop bulk 1.88e6 shear 0.92e6fric 14.6coh 39.7e3tens 0range z 052;PA-10 (下部也用pa-10)
;-------------------------------------------------
set mech force=1e-3
solve
save 20_20.sav
;---------------分步开挖并循环至平衡--------------------
model null range z 50 56 plane dip 30 dd 270 origin 118.6168,0,56 above
set mech force=1e-3
solve
save 20_30.sav
model null range z 50 56 plane dip 40 dd 270 origin 118.6168,0,56 above
set mech force=1e-3
solve
save 20_40.sav
model null range z 50 56 plane dip 50 dd 270 origin 118.6168,0,56 above
set mech force=1e-3
solve
save 20_50.sav
model null range z 50 56 plane dip 60 dd 270 origin 118.6168,0,56 above
set mech force=1e-3
solve
save 20_60.sav
model null range z 50 56 plane dip 70 dd 270 origin 118.6168,0,56 above
set mech force=1e-3
solve
save 20_70.sav
model null range z 50 56 plane dip 80 dd 270 origin 118.6168,0,56 above
set mech force=1e-3
solve
save 20_80.sav
;------------------使用强度折减法计算稳定系数--------------------------
restore 20_80.sav
solve fos file 20_80fos.sav associated
C:\Users\zh\Desktop\图片
这个是模型平衡后的图片，感觉显示的有点问题啊，求指导！
页:
[1]