6-1 V
级围岩施工过程模拟（1
）建立隧道1/4
圆周模型gen zonradcyl p0 0 0 0 p1 10.55 0 0 p2 0 1 0 p30 0 10.55 &
size 5 2 10 4 dim 5.55 5.55 5.555.55 rat 1 1 1 1.2 group outsiderocksy
plotblock group
plot addaxes red
gen zonecshell 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 11 1 group concretlinersy fill group insiderocksy
gen zonradcyl 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.55rat 1 1 1 1.2 group outsiderockxy
gen zonecshell 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 11 group concretlinerxy fill group insiderockxy
plotblock group
plot addaxes red
gen zonbrick p0 0 0 10.55 p1 10.55 0 10.55 p2 0 1 10.55 p3 0 0 25.55 &
size 5 28 group outsiderock1
gen zonbrick p0 0 0 -34.55 p1 10.55 0 -34.55 p2 0 1 -34.55 p3 0 0 -9 &
size 5 212 group outsiderock2
gen zonbrick p0 10.55 0 0 p1 50 0 0 p2 10.55 1 0 p3 10.55 0 10.55&
size 20 25 group outsiderock3
gen zonbrick p0 10.55 0 -9 p1 50 0 -9 p2 10.55 1 -9 p3 10.55 0 0 &
size 20 25 group outsiderock3
gen zonbrick 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 28 group outsiderock4
gen zonbrick p0 10.55 0 -34.55 p1 50 0 -34.55 p2 10.55 1 -34.55 p3 10.55 0 -9 &
size 20 212 group outsiderock5
gen zonereflect norm -1 0 0 orig 0 0 0
group
隧道右上 range gro insiderocksy z 0 5.05 x 0 5.05
group
隧道左上 range gro insiderocksy z 0 5.05 x 0 -5.05
group
隧道右下 range gro insiderockxy z 0 -5.05 x 0 5.05
group
隧道左下 range gro insiderockxy z 0 -5.05 x 0 -5.05
group
衬砌右上 range gro concretlinersy z 0 5.55 x 0 5.55
group
衬砌左上 range gro concretlinersy z 0 5.55 x 0 -5.55
group
衬砌右下 range gro concretlinerxy z 0 -5.55 x 0 5.55
group
衬砌左下 range gro concretlinerxy z 0 -5.55 x 0 -5.55
setgravity 0 0 -10 ;
设置重力加速度为z
方向-10
fix zrange z -34.56,-34.54 ;
设置底边界fix xrange x -50.01,-49.99 ;
设置左边界fix xrange x 49.99,50.01 ;
设置右边界fix yrange y -0.01 0.01 ;
设置前边界fix yrange y 0.99,1.01 ;
设置后边界model mohr ;
莫尔～库仑模型ini density 2000 ;
围岩的密度propbulk=7.14e8 shear=3.333e9 fric=25 coh=2e5 tension=1.0e5
Step6000 ;
求解6000
次Solve ;
求解计算 save V
级围岩.sav ;
计算结果保存在V
级围岩.sav
文件中plot cont zdisp ;
绘制竖向位移场，如图6-14
所示plot cont xdisp ;
绘制竖向位移场，如图6-15
所示plot cont szz ;
绘制竖向应力场，如图6-16
所示plot cont sxx ;
绘制竖向应力场，如图6-17
所示;1-
左上半断面开挖rest V
级围岩.sav
plot block group
plot addaxes red
ini xdis0.0 ydis 0.0 zdis 0.0
inidensity 2200 range group outsiderocksy any
prop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 &
tension=5e5, range group outsiderocksy
modelnull range group
隧道左上 any group
衬砌左上 any
step 2000
save stepp1.sav
plot contzdisp
rest stepp1.sav
sel shellid=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 shellid=1 range x -0.01 0.01 z -0.01 5.56
sel shellid=1 prop iso=(21.0e9,0.2) thick=0.25
step 2000
plot contzdisp
save stepp2.sav
;3-
左下半断面开挖rest stepp2.sav
plot block group
plot addaxes red
model null range group
隧道左下 any group
衬砌左下 any
step 2000
savestepp3.sav
plot contzdisp
rest stepp3.sav
sel shellid=1 range x -0.01 0.01 z -4.01 0.01
sel shellid=1 prop iso=(21.0e9,0.2) thick=0.25
step 2000
plot contzdisp
save stepp4.sav
;5-
右上半断面开挖rest stepp4.sav
plotblock group
plot addaxes red
model null range group
隧道右上 any group
衬砌右上 any
step 2000
savestepp5.sav
plot contzdisp
rest stepp5.sav
sel shellid=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 shellid=1 prop iso=(21.0e9,0.2) thick=0.25
step 2000
plot contzdisp
savestepp6.sav
;7-
右下半断面开挖rest stepp6.sav
plot block group
plot addaxes red
model null range group
隧道右下 any group
衬砌右下 any
step 2000
save stepp7.sav
plot contzdisp
rest stepp7.sav
seldelete shell id=1 range x -0.01 0.01
sel shellid=1 prop iso=(21.0e9,0.2) thick=0.25
step 2000
save stepp8.sav
plo blogro
plot contzdisp
；9-
二次衬砌reststepp8.sav
model elas range group
衬砌左上model elas range group
衬砌左下model elas range group
衬砌右上model elas range group
衬砌右下ini density 2600 range group
衬砌左上ini density 2600 range group
衬砌左下ini density 2600 range group
衬砌右上ini density 2600 range group
衬砌右下prop bulk=12.5e9, shear=16.667e9 range group
衬砌左上prop bulk=12.5e9, shear=16.667e9 range group
衬砌左下prop bulk=12.5e9, shear=16.667e9 range group
衬砌右上prop bulk=12.5e9, shear=16.667e9 range group
衬砌右下step 2000
plot contzdisp
save stepp9.sav
plot block group
plot addaxes red
