sel shell id=1 range
new;1/4 上半部分
gen zon radcylp0 15.7 0 0 p1 24.5 0 0 p2 15.7 1 0 p3 15.7 0 8.8 &
size 5 2 10 4 dim 5.8 5.8 5.8 5.8 rat 1 1 1 1.2 group outrockys
gen zone cshell p0 15.7 0 0 p1 21.5 0 0 p2 15.7 1 0 p3 15.7 0 5.8 &
size 1 2 10 4 dim 5.4 5.4 5.4 5.4 rat 1 1 1 1 group liningys fill group inrockys
gen zon reflect dip 90 dd 270 orig 15.7 0 0
;1/4 下半部分
gen zon radcylp0 15.7 0 0 p1 15.7 0 -6.9p2 15.7 1 0 p3 24.5 0 0 &
size 5 2 10 4 dim 3.9 5.8 3.9 5.8 rat 1 1 1 1.2 group outrockyx
gen zone cshell p0 15.7 0 0 p1 15.7 0 -3.9 p2 15.7 1 0 p3 21.5 0 0 &
size 1 2 10 4 dim 3.5 5.4 3.5 5.4 rat 1 1 1 1 group liningyx fill group inrockyx
gen zon reflect dip 90 dd 270 orig 15.7 0 0 range group liningyx any &
group inrockyx any group outrockyx any
;右隧道中地层网格
gen zon brickp0 0 0 -6.9p1 6.9 0 -6.9p2 0 1 -6.9p3 0 0 0size 6 2 5group outrock1
gen zon brickp0 24.5 0 -6.9p1 65 0 -6.9p2 24.5 1 -6.9p3 24.5 0 0size 20 2 5group outrock1
gen zon brickp0 0 0 0p1 6.9 0 0p2 0 1 0p3 0 0 8.8size 6 2 5group outrock2
gen zon brickp0 24.5 0 0p1 65 0 0p2 24.5 1 0p3 24.5 0 8.8size 20 2 5group outrock2
;右隧道上侧土体网格
gen zon brickp0 0 0 8.8p1 6.9 0 8.8p2 0 1 8.8p3 0 0 36size 6 2 10group outrock3
gen zon brickp0 24.5 0 8.8p1 65 0 8.8p2 24.5 1 8.8p3 24.5 0 36size 20 2 10group outrock3
gen zon brickp0 6.9 0 8.8p1 24.5 0 8.8p2 6.9 1 8.8p3 6.9 0 36size 10 2 10group outrock3
;右隧道下侧土体网格
gen zon brickp0 0 0 -24p1 6.9 0 -24p2 0 1 -24p3 0 0 -6.9size 6 2 8group outrock4
gen zon brickp0 24.5 0 -24p1 65 0 -24p2 24.5 1 -24p3 24.5 0 -6.9size 20 2 8group outrock4
gen zon brickp0 6.9 0 -24p1 24.5 0 -24p2 6.9 1 -24p3 6.9 0 -6.9size 10 2 8group outrock4
gen zon reflect dip 90 dd 270 orig 0 0 0
plot add axes red
plot block group
;设置重力加速度
set gravity 0 0 -10
;给定边界条件
fix z range z -24.01,-23.99
fix x range x -65.01,-64.99
fix x range x 64.99,65.01
fix y range y -0.01 0.01
fix y range y 0.99,1.01
;0-求解自重应力场
model mohr
ini density 2200
;围岩的密度
prop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 tension=5e5
;体积、剪切、摩擦角、凝聚力、抗拉强度
step 6000
save Gsol.sav
plot cont zdisp
plot cont szz
;1-左隧道上半断面开挖
rest Gsol.sav
ini xdis 0.0 ydis 0.0 zdis 0.0
ini density 2400 range group outrockzs any
prop bulk=6.0e9, shear=10.0e9 fric=45, coh=1.2e6 tension=2e6 &
range group outrockzs any
model null range group inrockzs any group liningzs any
step 2000
save stepp1.sav
;2-初期支护
rest stepp1.sav
new
;1/4 上半部分
gen zon radcylp0 15.7 0 0 p1 24.5 0 0 p2 15.7 1 0 p3 15.7 0 8.8 &
size 5 2 10 4 dim 5.8 5.8 5.8 5.8 rat 1 1 1 1.2 group outrockys
gen zone cshell p0 15.7 0 0 p1 21.5 0 0 p2 15.7 1 0 p3 15.7 0 5.8 &
size 1 2 10 4 dim 5.4 5.4 5.4 5.4 rat 1 1 1 1 group liningys fill group inrockys
gen zon reflect dip 90 dd 270 orig 15.7 0 0
;1/4 下半部分
gen zon radcylp0 15.7 0 0 p1 15.7 0 -6.9p2 15.7 1 0 p3 24.5 0 0 &
size 5 2 10 4 dim 3.9 5.8 3.9 5.8 rat 1 1 1 1.2 group outrockyx
gen zone cshell p0 15.7 0 0 p1 15.7 0 -3.9 p2 15.7 1 0 p3 21.5 0 0 &
size 1 2 10 4 dim 3.5 5.4 3.5 5.4 rat 1 1 1 1 group liningyx fill group inrockyx
gen zon reflect dip 90 dd 270 orig 15.7 0 0 range group liningyx any &
group inrockyx any group outrockyx any
;右隧道中地层网格
gen zon brickp0 0 0 -6.9p1 6.9 0 -6.9p2 0 1 -6.9p3 0 0 0size 6 2 5group outrock1
gen zon brickp0 24.5 0 -6.9p1 65 0 -6.9p2 24.5 1 -6.9p3 24.5 0 0size 20 2 5group outrock1
gen zon brickp0 0 0 0p1 6.9 0 0p2 0 1 0p3 0 0 8.8size 6 2 5group outrock2
gen zon brickp0 24.5 0 0p1 65 0 0p2 24.5 1 0p3 24.5 0 8.8size 20 2 5group outrock2
;右隧道上侧土体网格
gen zon brickp0 0 0 8.8p1 6.9 0 8.8p2 0 1 8.8p3 0 0 36size 6 2 10group outrock3
gen zon brickp0 24.5 0 8.8p1 65 0 8.8p2 24.5 1 8.8p3 24.5 0 36size 20 2 10group outrock3
gen zon brickp0 6.9 0 8.8p1 24.5 0 8.8p2 6.9 1 8.8p3 6.9 0 36size 10 2 10group outrock3
;右隧道下侧土体网格
gen zon brickp0 0 0 -24p1 6.9 0 -24p2 0 1 -24p3 0 0 -6.9size 6 2 8group outrock4
gen zon brickp0 24.5 0 -24p1 65 0 -24p2 24.5 1 -24p3 24.5 0 -6.9size 20 2 8group outrock4
gen zon brickp0 6.9 0 -24p1 24.5 0 -24p2 6.9 1 -24p3 6.9 0 -6.9size 10 2 8group outrock4
gen zon reflect dip 90 dd 270 orig 0 0 0
plot add axes red
plot block group
;设置重力加速度
set gravity 0 0 -10
;给定边界条件
fix z range z -24.01,-23.99
fix x range x -65.01,-64.99
fix x range x 64.99,65.01
fix y range y -0.01 0.01
fix y range y 0.99,1.01
;0-求解自重应力场
model mohr
ini density 2200
;围岩的密度
prop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 tension=5e5
;体积、剪切、摩擦角、凝聚力、抗拉强度
step 6000
save Gsol.sav
plot cont zdisp
plot cont szz
;1-左隧道上半断面开挖
rest Gsol.sav
ini xdis 0.0 ydis 0.0 zdis 0.0
ini density 2400 range group outrockzs any
prop bulk=6.0e9, shear=10.0e9 fric=45, coh=1.2e6 tension=2e6 &
range group outrockzs any
model null range group inrockzs any group liningzs any
step 2000
save stepp1.sav
;2-初期支护
rest stepp1.sav
sel shell id=1 range x -21.51 -9.89 cyl end1 -15.7 0 0 end2 -15.7 1 0 rad 5.8
sel shell id=1 prop iso=(21.0e9,0.2) thick=0.20
step 2000
save stepp2.sav
sel shell id=1 prop iso=(21.0e9,0.2) thick=0.20
step 2000
save stepp2.sav
sel shell id=1 range x -21.51 -9.89 cyl end1 -15.7 0 0 end2 -15.7 1 0 rad 5.8 range x -21.51 -9.89 这个作用是什么，谢谢
页:
[1]