找回密码
 注册
Simdroid-非首页
查看: 126|回复: 1

[工程实例] 一个简单边坡的模拟,但感觉有点问题,还请大家指点一下,哪儿有问题?谢谢

[复制链接]
发表于 2015-8-31 15:03:11 | 显示全部楼层 |阅读模式 来自 青海西宁
以下是命令:

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  A-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 0  52  ;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.08e6  shear 5.41e6  fric 14.8  coh 32.1e3  tens 0  range z 68 75  ;PA-01 (上部也用pa-01)
prop bulk 10.11e6  shear 4.94e6  fric 13.4  coh 48.4e3  tens 0  range z 66 68  ;PA-02
prop bulk 5.81e6   shear 2.84e6  fric 13.2  coh 48.3e3  tens 0  range z 64 66  ;PA-03
prop bulk 9.27e6   shear 4.53e6  fric 16.1  coh 54.3e3  tens 0  range z 62 64  ;PA-04
prop bulk 10.54e6  shear 5.15e6  fric 15.6  coh 55.0e3  tens 0  range z 60 62  ;PA-05
prop bulk 11.39e6  shear 5.56e6  fric 14.6  coh 59.2e3  tens 0  range z 58 60  ;PA-06
prop bulk 10.05e6  shear 4.91e6  fric 16.1  coh 55.5e3  tens 0  range z 56 58  ;PA-07
prop bulk 4.06e6   shear 1.98e6  fric 14.1  coh 45.6e3  tens 0  range z 54 56  ;PA-08
prop bulk 2.91e6   shear 1.42e6  fric 13.7  coh 40.6e3  tens 0  range z 52 54  ;PA-09
prop bulk 1.88e6   shear 0.92e6  fric 14.6  coh 39.7e3  tens 0  range z 0  52  ;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\图片

 楼主| 发表于 2015-8-31 15:05:27 | 显示全部楼层 来自 青海西宁
Simdroid开发平台

这个是模型平衡后的图片,感觉显示的有点问题啊,求指导!

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-9-23 20:13 , Processed in 0.027716 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表