第6章 命令流 隧道及地下工程解析方法 李围

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

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