gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 & size 4 10 6 4 dim 4 4 4 4 rat 1 1 1 1 fill gen zon reflect dip 90 dd 90 orig 15 0 0 gen zon reflect dip 0 ori 0 0 0 ;辅助隧道建模
gen point id 1 (2.969848,0.0,-0.575736) gen point id 2 (2.969848,50.0,-0.575736)
gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 & p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 point 1 p10 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fill
gen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 & p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 point 1 p11 point 2 & size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fill ;隧道边界网格
gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 & p5 15 50 50 p6 50 0 50 p7 50 50 50 &
p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 & size 6 10 3 10 rat 1 1 1 1.1
gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 & p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 & size 3 10 10 rat 1 1 1.1
gen zon reflect dip 0 ori 0 0 0 range x 0 23 y 0 50 z 8 50 gen zon reflect dip 0 ori 0 0 0 range x 23 50 y 0 50 z 0 50 group main1 range cyl end1 15 0 0 end2 15 25 0 rad 4 group main2 range cyl end1 15 25 0 end2 15 50 0 rad 4 group service range cyl end1 0.0,0.0,-0.575736 & end2 0.0,50.0,-0.575736 rad 3.0 save tun0.sav rest tun0.sav
; 初始应力场模量 model mohr
prop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5 apply szz -1.4e6 range z 49.9 50.1 fix z range z -50.1 -49.1 fix x range x -.1 .1 fix x range x 49.9 50.1 fix y range y -.1 .1 fix y range y 49.9 50.1
ini sxx -1.4e6 syy -1.4e6 szz -1.4e6 hist unbal
hist gp xdis 3,0,-1 hist gp zdis 0,0,2 hist gp xdis 3,25,-1 hist gp zdis 0,25,2
step 1000 save tun1.sav
; 辅助隧道开挖25 m计算 rest tun1.sav
ini xdis 0.0 ydis 0.0 zdis 0.0
model null range group service y 0,25 step 1000 save tun2.sav
; 在辅助隧道上施加衬砌 restore tun2.sav
sel shell id=1 range cyl end1 0 0 -1 end2 0 25 -1 rad 3 sel shell prop iso=(25.3e9, 0.266) thick = 0.5
sel node fix y xr zr range y -0.1 0.1 ; symmetry cond. sel node fix x yr zr range x -0.1 0.1 ; symmetry cond. model mohr range group service y 25 50
prop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5 ; 主隧道开挖25m
model null range group main1 step 1000 save tun3.sav
plot set rotation 0 0 60 plot cont zdisp outl on
6-1 V级围岩施工过程模拟
（1）建立隧道1/4圆周模型
gen zon radcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p3 0 0 10.55 &
size 5 2 10 4 dim 5.55 5.55 5.55 5.55 rat 1 1 1 1.2 group outsiderocksy plot block group plot add axes red
gen zone cshell p0 0 0 0 p1 5.55 0 0 p2 0 1 0 p3 0 0 5.55 &
size 1 2 10 4 dim 5.05 5.05 5.05 5.05 rat 1 1 1 1 group concretlinersy fill group insiderocksy gen zon radcyl p0 0 0 0 p1 0 0 -9 p2 0 1 0 p3 10.55 0 0 &
size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxy gen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 &
size 1 2 10 4 dim 3.5 5.05 3.5 5.05 rat 1 1 1 1 group concretlinerxy fill group insiderockxy plot block group plot add axes red
gen zon brick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 & size 5 2 8 group outsiderock1
gen zon brick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 & size 5 2 12 group outsiderock2
gen zon brick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55 & size 20 2 5 group outsiderock3
gen zon brick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 &
size 20 2 5 group outsiderock3
gen zon brick p0 10.55 0 10.55 p1 50 0 10.55 p2 10.55 1 10.55 p3 10.55 0 25.55 & size 20 2 8 group outsiderock4
gen zon brick p0 10.55 0 -34.55 p1 50 0 -34.55 p2 10.55 1 -34.55 p3 10.55 0 -9 & size 20 2 12 group outsiderock5
gen zon reflect dip 90 dd 90 orig 0 0 0
plot set rotation 30 0 30 ;显示三维图，绕x轴负方向转30度 set gravity 0 0 -10 ;设置重力加速度为z方向-10 fix z range z -34.56,-34.54 ;设置底边界 fix x range x -50.01,-49.99 ;设置左边界 fix x range x 49.99,50.01 ;设置右边界 fix y range y -0.01 0.01 ;设置前边界 fix y range y 0.99,1.01 ;设置后边界 model mohr ;莫尔～库仑模型 ini density 2000 ;围岩的密度
prop bulk=7.14e8 shear=3.333e9 fric=25 coh=2e5 tension=1.0e5 Setp 6000 ;求解6000次 Solve ;求解计算
save Gsol.sav ;计算结果保存在Gsol.sav文件中 plot cont zdisp ;绘制竖向位移场，如图6-14所示 plot cont xdisp ;绘制竖向位移场，如图6-15所示 plot cont szz ;绘制竖向应力场，如图6-16所示 plot cont sxx ;绘制竖向应力场，如图6-17所示 ;1-左上半断面开挖 rest Gsol.sav plot block group plot add axes red
ini xdis 0.0 ydis 0.0 zdis 0.0
ini density 2200 range group outsiderocksy any group outsiderocksz any prop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 &
tension=5e5, range group outsiderocksy any group outsiderocksz any model null range group insiderocksz any group concretlinersz any step 2000
save stepp1.sav plot cont zdisp rest stepp1.sav
sel shell id=1 range x -0.01 0.01 z -0.01 5.56 cyl end1 0 0 0 end2 0 1 0 rad 5.55 sel shell id=1 range x -0.01 0.01 z -0.01 5.56 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25 step 2000
plot cont zdisp save stepp2.sav ;3-左下半断面开挖 rest stepp2.sav
plot block group plot add axes red
model null range group insiderockxz any group concretlinerxz any step 2000
save stepp3.sav plot cont zdisp rest stepp3.sav
sel shell id=1 range x -0.01 0.01 z -4.01 0.01 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25 step 2000
plot cont zdisp save stepp4.sav ;5-右上半断面开挖 rest stepp4.sav plot block group plot add axes red
model null range group insiderocksy any group concretlinersy any step 2000
save stepp5.sav plot cont zdisp rest stepp5.sav
sel shell id=1 range x -0.01 5.56 z -0.01 5.56 cyl end1 0 0 0 end2 0 1 0 rad 5.55 sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25 step 2000
plot cont zdisp save stepp6.sav ;7-右下半断面开挖 rest stepp6.sav plot block group plot add axes red
model null range group insiderockxy any group concretlinerxy any step 2000
save stepp7.sav plot cont zdisp rest stepp7.sav
sel delete shell id=1 range x -0.01 0.01
sel shell id=1 prop iso=(21.0e9,0.2) thick=0.25 step 2000
save stepp8.sav rest stepp8.sav
model elas range group concretlinerxy model elas range group concretlinerxz model elas range group concretlinersy model elas range group concretlinersz
ini density 2600 range group concretlinerxy ini density 2600 range group concretlinerxz ini density 2600 range group concretlinersy ini density 2600 range group concretlinersz
prop bulk=12.5e9, shear=16.667e9 range group concretlinerxy prop bulk=12.5e9, shear=16.667e9 range group concretlinerxz prop bulk=12.5e9, shear=16.667e9 range group concretlinersy prop bulk=12.5e9, shear=16.667e9 range group concretlinersz step 2000
plot cont zdisp save stepp9.sav plot block group plot add axes red
6.2 IV级围岩施工过程模拟
new
;1/4 上半部分
gen zon radcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p3 0 0 10.55 size 5 2 10 4& dim 5.55 5.55 5.55 5.55 rat 1 1 1 1.2 group outsiderocksy
gen zone cshell p0 0 0 0 p1 5.55 0 0 p2 0 1 0 p3 0 0 5.55 size 1 2 10 4&
dim 5.15 5.15 5.15 5.15 rat 1 1 1 1 group concretlinersy fill group insiderocksy ;1/4 下半部分
gen zon radcyl p0 0 0 0 p1 0 0 -9 p2 0 1 0 p3 10.55 0 0 &
size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxy gen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 size 1 2 10 4&
dim 3.6 5.15 3.6 5.15 rat 1 1 1 1 group concretlinerxy fill group insiderockxy ;1/2绘制上下地层网格
gen zon brick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 size 5 2 8 group outsiderock1
gen zon brick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 size 5 2 12 group outsiderock2
;1/2右侧中间部分土体网格
gen zon brick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55 size 20 2 5 group outsiderock3
gen zon brick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 size 20 2 5 group outsiderock3
;1/2右上侧和下侧土体网格
gen zon brick p0 10.55 0 10.55 p1 50 0 10.55 p2 10.55 1 10.55 p3 10.55 0 25.55 size 20 2 8 group outsiderock4
gen zon brick p0 10.55 0 -34.55 p1 50 0 -34.55 p2 10.55 1 -34.55 p3 10.55 0 -9 size 20 2 12 group outsiderock5 ;对称得出完整模型
gen zon reflect dip 90 dd 90 orig 0 0 0