sel shell id=1 range

zsq

new
;1/4 上半部分
gen zon radcyl  p0 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 radcyl  p0 15.7 0 0 p1 15.7 0 -6.9  p2 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 brick  p0 0 0 -6.9  p1 6.9 0 -6.9  p2 0 1 -6.9  p3 0 0 0  size 6 2 5  group outrock1
gen zon brick  p0 24.5 0 -6.9  p1 65 0 -6.9  p2 24.5 1 -6.9  p3 24.5 0 0  size 20 2 5  group outrock1
gen zon brick  p0 0 0 0  p1 6.9 0 0  p2 0 1 0  p3 0 0 8.8  size 6 2 5  group outrock2
gen zon brick  p0 24.5 0 0  p1 65 0 0  p2 24.5 1 0  p3 24.5 0 8.8  size 20 2 5  group outrock2
;右隧道上侧土体网格
gen zon brick  p0 0 0 8.8  p1 6.9 0 8.8  p2 0 1 8.8  p3 0 0 36  size 6 2 10  group outrock3
gen zon brick  p0 24.5 0 8.8  p1 65 0 8.8  p2 24.5 1 8.8  p3 24.5 0 36  size 20 2 10  group outrock3
gen zon brick  p0 6.9 0 8.8  p1 24.5 0 8.8  p2 6.9 1 8.8  p3 6.9 0 36  size 10 2 10  group outrock3
;右隧道下侧土体网格
gen zon brick  p0 0 0 -24  p1 6.9 0 -24  p2 0 1 -24  p3 0 0 -6.9  size 6 2 8  group outrock4
gen zon brick  p0 24.5 0 -24  p1 65 0 -24  p2 24.5 1 -24  p3 24.5 0 -6.9  size 20 2 8  group outrock4
gen zon brick  p0 6.9 0 -24  p1 24.5 0 -24  p2 6.9 1 -24  p3 6.9 0 -6.9  size 10 2 8  group 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 radcyl  p0 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 radcyl  p0 15.7 0 0 p1 15.7 0 -6.9  p2 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 brick  p0 0 0 -6.9  p1 6.9 0 -6.9  p2 0 1 -6.9  p3 0 0 0  size 6 2 5  group outrock1
gen zon brick  p0 24.5 0 -6.9  p1 65 0 -6.9  p2 24.5 1 -6.9  p3 24.5 0 0  size 20 2 5  group outrock1
gen zon brick  p0 0 0 0  p1 6.9 0 0  p2 0 1 0  p3 0 0 8.8  size 6 2 5  group outrock2
gen zon brick  p0 24.5 0 0  p1 65 0 0  p2 24.5 1 0  p3 24.5 0 8.8  size 20 2 5  group outrock2
;右隧道上侧土体网格
gen zon brick  p0 0 0 8.8  p1 6.9 0 8.8  p2 0 1 8.8  p3 0 0 36  size 6 2 10  group outrock3
gen zon brick  p0 24.5 0 8.8  p1 65 0 8.8  p2 24.5 1 8.8  p3 24.5 0 36  size 20 2 10  group outrock3
gen zon brick  p0 6.9 0 8.8  p1 24.5 0 8.8  p2 6.9 1 8.8  p3 6.9 0 36  size 10 2 10  group outrock3
;右隧道下侧土体网格
gen zon brick  p0 0 0 -24  p1 6.9 0 -24  p2 0 1 -24  p3 0 0 -6.9  size 6 2 8  group outrock4
gen zon brick  p0 24.5 0 -24  p1 65 0 -24  p2 24.5 1 -24  p3 24.5 0 -6.9  size 20 2 8  group outrock4
gen zon brick  p0 6.9 0 -24  p1 24.5 0 -24  p2 6.9 1 -24  p3 6.9 0 -6.9  size 10 2 8  group 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 这个作用是什么，谢谢