- 积分
- 0
- 注册时间
- 2008-11-16
- 仿真币
-
- 最后登录
- 1970-1-1
|
;new
;
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 0.5 0 p3 0 0 20 size 50 1 15
gen zone brick p0 15 0 20 p1 23 0 20 p2 15 0.5 20 p3 23 0 30 p4 26 0.5 20 p5 23 0.5 30 p6 26 0 30 p7 26 0.5 30 size 8 1 10 ratio 1.03 1 1
gen zone brick p0 23 0 20 p1 60 0 20 p2 26 0.5 20 p3 26 0 30 p4 60 0.5 20 p5 26 0.5 30 p6 60 0 30 p7 60 0.5 30 size 30 1 10 ratio 1.03 1 1
;
fix x y z range z -0.1 0.1
fix x range x 59.9 60.1
fix x range x -0.1 0.1
fix y
;
model elas
prop density 1600 bulk 8.3e9 shear 3.8e9 range z 28 30
prop density 1700 bulk 16.6e9 shear 7.6e9 range z 18 28
prop density 1900 bulk 33.2e9 shear 15.2e9 range z 10 18
prop density 2200 bulk 16.6e10 shear 7.6e10 range z 0 10
set gravity 0 0 -10.0
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
;
sel cable id=1 begin=(25.7,0,29) end=(28.65,0,28.48) nseg=1
sel cable id=2 begin=(25.1,0,27) end=(28.05,0,26.48) nseg=1
sel cable id=3 begin=(24.5,0,25) end=(27.45,0,24.48) nseg=1
sel cable id=4 begin=(23.9,0,23) end=(26.85,0,22.48) nseg=1
sel cable id=5 begin=(23.3,0,21) end=(26.25,0,20.48) nseg=1
sel cable prop xcarea=7.06e-4 emod=200e9 ytens=1e10 gr_k=7e6 gr_coh=1e9
;
sel cable prop gr_per=0.314 gr_fric=25
;
plot add disp max=2e-2
plot add surface yellow
plot add sel cable force red
plot add sel cable grout slip
plot add sel geom node=off black black
Plot show
;
step 5000
;
save cable-sn.sav
;
model mohr
prop density 1600 bulk 8.3e9 shear 3.8e9 coh 20000 tens 1.0e6 fric 18 dila 18 range z 28 30
prop density 1700 bulk 16.6e9 shear 7.6e9 coh 22000 tens 1.0e6 fric 22 dila 22 range z 18 28
prop density 1900 bulk 33.2e9 shear 15.2e9 coh 32000 tens 1.0e6 fric 24 dila 24 range z 10 18
prop density 2200 bulk 16.6e10 shear 7.6e10 coh 200000 tens 1.0e6 fric 30 dila 30 range z 0 10
solve fos file slope3dfos1.sav associated
;
plot fos cont ssi outline on vel red
;
求助 为什么算到最后一步边坡系数是就不运行了 ? |
|