内容发布更新时间 : 2024/11/6 8:16:59星期一 下面是文章的全部内容请认真阅读。

第七章 命令流按照顺序进行

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 groupoutrock4

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

;3-左隧道下半断面开挖 rest stepp2.sav

model null range group inrockzx any group liningzx any step 2000

save stepp3.sav plot add axes red plot block group plot cont zdisp plot cont szz

;4-左隧道二次衬砌 rest stepp3.sav

model elas range group liningzs model elas range group liningzx

ini density 2600 range group liningzs ini density 2600 range group liningzx

prop bulk=12.5e9, shear=16.667e9 range group liningzs prop bulk=12.5e9, shear=16.667e9 range group liningzx

step 2000

save stepp4.sav plot add axes red plot block group plot cont zdisp plot cont szz

;5-右隧道上半断面开挖 rest stepp4.sav

ini density 2400 range group outrockys any

prop bulk=6.0e9, shear=10.0e9 fric=45, coh=1.2e6 tension=2e6 & range group outrockys any

model null range group inrockys any group liningys any step 2000

save stepp5.sav plot cont zdisp ;6-初期支护 rest stepp5.sav

sel shell id=1 range x 9.89 21.51 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

plot cont zdisp save stepp6.sav plot add axes red plot block group plot cont zdisp plot cont szz

;7-左隧道下半断面开挖 rest stepp6.sav

model null range group inrockyx any group liningyx any step 2000

save stepp7.sav plot cont zdisp

;8-左隧道二次衬砌 rest stepp7.sav

model elas range group liningys model elas range group liningyx

ini density 2600 range group liningys ini density 2600 range group liningyx

prop bulk=12.5e9, shear=16.667e9 range group liningys prop bulk=12.5e9, shear=16.667e9 range group liningyx step 2000

save stepp8.sav plot add axes red plot block group