level 1
Lm_心痛
楼主
如果是step太多,那么开挖后隧道位移已经平衡,那么加完锚杆以后应该位移不会再变了,结果怎么位移由70多厘米变成90厘米
程序如下:
NEW
gen zone brick p0 0 0 0 p1 10 0 0 p2 0 6 0 p3 0 0 10 size 10 6 10
plot reset
model elastic
prop shear =1.36e10 bulk=9.14e9
ini dens 2500
fix x ran x -0.1 0.1
fix x ran x 9.9 10.1
fix y ran y -0.1 0.1
fix y ran y 5.9 6.1
fix x y z ran z -0.1 0.1
app szz 1 ran z 9.9 10.1
ini szz 1
ini sxx 1
ini syy 1
set grav 0 0 -9.8
hist unbal
solve
save pingheng.sav
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model mo
prop bulk=1e8 shear=0.3e8 fric=35 coh=1e3 tens=1e3
model null range x 3 7 y 2 4 z 6 10
step 3000
save kaiwa.sav
sel cable id=1 beg 7 3 7 end 9 3 7 nseg 2
sel cable id=1 beg 7 3 8 end 9 3 8 nseg 2
sel cable id=1 beg 7 3 9 end 9 3 9 nseg 2
sel cable id =1 prop emod 2e9 ytens 1e8 xcarea 0.0314 gr_coh 1e6 gr_k 2e6 gr_per 1.0
hist unbal
step 3000
save jisuan.sav
2017年06月14日 04点06分
1
程序如下:
NEW
gen zone brick p0 0 0 0 p1 10 0 0 p2 0 6 0 p3 0 0 10 size 10 6 10
plot reset
model elastic
prop shear =1.36e10 bulk=9.14e9
ini dens 2500
fix x ran x -0.1 0.1
fix x ran x 9.9 10.1
fix y ran y -0.1 0.1
fix y ran y 5.9 6.1
fix x y z ran z -0.1 0.1
app szz 1 ran z 9.9 10.1
ini szz 1
ini sxx 1
ini syy 1
set grav 0 0 -9.8
hist unbal
solve
save pingheng.sav
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model mo
prop bulk=1e8 shear=0.3e8 fric=35 coh=1e3 tens=1e3
model null range x 3 7 y 2 4 z 6 10
step 3000
save kaiwa.sav
sel cable id=1 beg 7 3 7 end 9 3 7 nseg 2
sel cable id=1 beg 7 3 8 end 9 3 8 nseg 2
sel cable id=1 beg 7 3 9 end 9 3 9 nseg 2
sel cable id =1 prop emod 2e9 ytens 1e8 xcarea 0.0314 gr_coh 1e6 gr_k 2e6 gr_per 1.0
hist unbal
step 3000
save jisuan.sav