- 积分
- 0
- 注册时间
- 2004-12-1
- 仿真币
-
- 最后登录
- 1970-1-1
|
new
title 'SLOP SIMULATION'
;set large
set log on
set mecha ratio 1e-5
gen zon brick p0 1 0 0 p1 1 28 0 p2 0 0 0 p3 1 0 20 &
size 56 2 40 ratio 1.0 1.0 1.0
set grav 0.0 0.0 -10.0
mod mohr range z 18.5 20.0
pro bulk=2.14e6 she=0.39e6 dens 1.70e3 coh 10e3 range z 18.5 20.0
prop friction 10 dilation 0. tension 84.4e3 range z 18.5 20.0
mod mohr range z 17.5 18.5
pro bulk=4.14e6 she=0.74e6 dens 1.89e3 coh 25e3 range z 17.5 18.5
prop friction 16.5 dilation 0. tension 84.4e3 range z 17.5 18.5
mod mohr range z 16.5 17.5
pro bulk=8.18e6 she=1.75e6 dens 1.91e3 coh 16e3 range z 16.5 17.5
prop friction 12.5 dilation 0. tension 72.2e3 range z 16.5 17.5
mod mohr range z 15.0 16.5
pro bulk=2.87e6 she=0.61e6 dens 1.70e3 coh 16e3 range z 15.0 16.5
prop friction 15.5 dilation 0. tension 57.5e3 range z 15.0 16.5
mod mohr range z 0.0 15.0
pro bulk=2.04e6 she=0.44e6 dens 1.70e3 coh 15e3 range z 0.0 15.0
prop friction 11.5 dilation 0. tension 73.7e3 range z 0.0 15.0
fix z range z -0.01 0.01
fix y range y 27.99 28.01
fix y range y -0.01 0.01
fix x range x -0.01 0.01
fix x range x 0.99 1.01
solve
;solve for balance
ini state 0
ini xdisp 0.0 ydisp 0.0 zdisp 0.0
save initial.sav
rest initial.sav
;set pile
mod mohr range y 10.0 11.0 z 10.0 20.0
pro bulk=5.5e7 she=1.2e7 dens 2.0e3 coh 120e3 range y 10.0 11.0 z 10.0 20.0
prop friction 25.0 dilation 0. tension 300e3 range y 10.0 11.0 z 10.0 20.0
mod null range y 0.0 10.0 z 18.5 20.0
sel cable id=1 begin=(0.5, 10.0, 18.5) end=(0.5, 21.8, 16.4) nseg=12
sel cable prop xcarea=5.4e-3 emod=200e9 yTens=210e6 &
gr_k=1e7 gr_coh=15e3 ;
solve
save step1.sav |
|