隧道及地下工程FLAC解析方法 - 命令流

2-1定义一个FISH函数 .................................................................................................................. 2 2-2使用一个变量 ............................................................................................................................ 2 2-3对变量和函数的理解 ................................................................................................................ 2 2-4获取变量的历史记录 ................................................................................................................ 2 2-5用FISH函数计算体积模量和剪砌模量 .................................................................................. 3 2-6 在FLAC输入中使用符号变量 ............................................................................................... 3 2-7 控制循环 ................................................................................................................................... 4 2-8 拆分命令行 ............................................................................................................................... 4 2-9 变量类型 ................................................................................................................................... 5 2-10 IF条件语句 .............................................................................................................................. 5 2-11 索单元自动生成 ..................................................................................................................... 5 2-12圆形隧道开挖模拟计算 .......................................................................................................... 6 4-1数组 ............................................................................................................................................ 8 4-2函数操作 .................................................................................................................................... 9 4-3函数删除与重定义 .................................................................................................................. 10 4-4字符串 ...................................................................................................................................... 10 4-5马蹄形隧道网格 ...................................................................................................................... 11 4-6复杂形状网格生成 .................................................................................................................. 11 4-7网格连接 .................................................................................................................................. 12 4-8立方体洞穴网格生成 .............................................................................................................. 13 4-9球体洞穴网格生成 .................................................................................................................. 14 4-10应力边界 ................................................................................................................................ 15 4-11改变应力边界条件 ................................................................................................................ 15 4-12位移边界 ................................................................................................................................ 17 4-13不考虑重力影响的均匀应力 ................................................................................................ 19 4-14考虑应力梯度的均匀材料 .................................................................................................... 19 4-15考虑应力梯度的非均匀材料 ................................................................................................ 20 4-16非均匀网格应力初始化 ........................................................................................................ 20 4-17不规则自由面应力初始化 .................................................................................................... 21 4-18非均网格内部压实 ................................................................................................................ 21 4-19模型改变后初始应力变化 .................................................................................................... 22 4-20应力与孔隙压力的初始化 .................................................................................................... 22 4-21加载顺序 ................................................................................................................................ 23 6-1 V级围岩施工过程模拟........................................................................................................... 25 6.2 IV级围岩施工过程模拟 ......................................................................................................... 28 6.3 III级围岩施工过程模拟 .......................................................................................................... 30 第七章 命令流按照顺序进行 ....................................................................................................... 32 ;右隧道中地层网格........................................................................................................................ 32 ;右隧道上侧土体网格.................................................................................................................... 32 ;右隧道下侧土体网格.................................................................................................................... 32 第七章 命令流按照顺序进行 ....................................................................................................... 35

2-1定义一个FISH函数

new def abc

abc = 25 * 3 + 5 End print abc

2-2使用一个变量

new def abc hh = 25

abc = hh * 3 + 5 End Print hh Print abc

2-3对变量和函数的理解

new def abc hh = 25

abc = hh * 3 + 5 End

set abc=0 hh=0 print hh print abc print hh new def abc

abc = hh * 3 + 5 end

set hh=25 print abc

set abc=0 hh=0 print hh print abc print hh

2-4获取变量的历史记录

new

gen zone brick size 1 2 1

model mohr

prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10 fix x y z range y -0.1 0.1

apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45 plot block group def get_ad

ad1 = gp_near(0,2,0) ad2 = gp_near(1,2,0) ad3 = gp_near(0,2,1) ad4 = gp_near(1,2,1) end get_ad def load

load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end hist load

hist gp ydis 0,2,0 step 1000

plot his 1 vs -2

2-5用FISH函数计算体积模量和剪砌模量

new

def derive

s_mod = y_mod / (2.0 * (1.0 + p_ratio))

b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end

set y_mod = 5e8 p_ratio = 0.25 derive

print b_mod print s_mod

2-6 在FLAC输入中使用符号变量

New

def derive

s_mod = y_mod / (2.0 * (1.0 + p_ratio))

b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end

set y_mod = 5e8 p_ratio = 0.25 derive

gen zone brick size 2,2,2 model elastic

prop bulk=b_mod shear=s_mod print zone prop bulk print zone prop shear

2-7 控制循环

New def xxx

sum = 0 prod = 1

loop n (1,10) sum = sum + n prod = prod * n end_loop end xxx

print sum, prod new

gen zone brick p0 (0,0,0) p1 (-10,0,0) p2 (0,10,0) p3 (0,0,-10) model elas

plot set rotation 0 0 45 plot block group def install

pnt = zone_head

loop while pnt #null z_depth = -z_zcen(pnt)

y_mod = y_zero + cc * sqrt(z_depth)

z_prop(pnt, ’shear’) = y_mod / (2.0*(1.0+p_ratio)) z_prop(pnt, ’bulk’) = y_mod / (3.0*(1.0-2.0*p_ratio)) pnt = zone_next(pnt) end_loop end

set p_ratio=0.25 y_zero=1e7 cc=1e8 install

2-8 拆分命令行

new ;example of a sum of many things def long_sum

temp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10 long_sum = temp + v11 + v12 + v13 + v14 + v15 end

2-9 变量类型

new

def haveone

aa = 2 bb = 3.4

cc = ’Have a nice day’ dd = aa * bb

ee = cc + ’, old chap’ end haveone print fish

2-10 IF条件语句

new def abc

if xx > 0 then abc = 1000 else

abc = -1000 end_if end

set xx = 10 print abc set xx = 0 print abc

2-11 索单元自动生成

new

gen zone brick size 10 3 5 plot set rotation 0 0 45 plot block group def place_cables

loop n (1,5)

z_d = float(n) - 0.5 command

sel cable beg 0.0,1.5,z_d end end_command end_loop end

place_cables

plot grid sel geom red

7.0,1.5,z_d nseg 7

end ;

moduli_data

gen zone brick size 2,2,2 model elastic

prop bulk = b_mod shear = s_mod print p_ratio y_mod b_mod s_mod pause

print zone prop bulk pause

print zone prop shear

4-5马蹄形隧道网格

new

gen zone radcyl size 5 10 6 12 rat 1 1 1 1.2 p0 0,0,0 p1 100,0,0 p2 0,200,0 p3 0,0,100 gen zone radtun size 5 10 5 12 rat 1 1 1 1.2 p0 0,0,0 p1 0,0,-100 p2 0,200,0 p3 100,0,0 gen zone reflect dip 90 dd 270 origin 0,0,0 plot set rotation 30 0 30 plot block group

plo add axes red

4-6复杂形状网格生成

new

; 主隧道网格

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

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 ;衬砌网格

sel shell range cyl end1 0 0 -1 end2 0 50 -1 rad 3 ; 隧道外围边界网格

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 ;对称得出1/2模型

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 service range cyl end1 0 0 -1 end2 0 50 -1 rad 3 group main range cyl end1 15 0 0 end2 15 50 0 rad 4 ;对称得出完整模型

gen zon reflect dip 90 dd 270 ori 0 0 0 ;网格显示

plot set rotation 30 0 30 plot block group

4-7网格连接

new

gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,2 gen zone brick size 8 8 4 p0 0,0,2 p1 4,0,2 p2 0,4,2 p3 0,0,4 attach face range z 1.9 2.1 model elas

prop bulk 8e9 shear 5e9 fix z range z -.1 .1 fix x range x -.1 .1 fix x range x 3.9 4.1 fix y range y -.1 .1 fix y range y 3.9 4.1

apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2 hist unbal ;网格显示

plot set rotation 30 0 30 plot block group solve

save att.sav

plot cont zdisp outl on

单一网格（小）程序 new

gen zone brick size 8 8 8 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elas

prop bulk 8e9 shear 5e9

fix z range z -.1 .1 fix x range x -.1 .1 fix x range x 3.9 4.1 fix y range y -.1 .1 fix y range y 3.9 4.1

apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2 hist unbal ;网格显示

plot set rotation 30 0 30 plot block group solve

save noattx.sav

plot cont zdisp outl on

单一网格（大）程序 new

gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,4 model elas

prop bulk 8e9 shear 5e9 fix z range z -.1 .1 fix x range x -.1 .1 fix x range x 3.9 4.1 fix y range y -.1 .1 fix y range y 3.9 4.1

apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2 hist unbal ;网格显示

plot set rotation 30 0 30 plot block group solve

save noattx.sav

plot cont zdisp outl on

4-8立方体洞穴网格生成

new def parm

rad=4.0 len=10.0 in_size=6 rad_size=10 end parm

gen zone radbrick edge len size in_size in_size in_size rad_size &

rat 1.0 1.0 1.0 1.2 dim rad rad rad

;对称得出完整模型

gen zon reflect dip 0 dd 0 ori 0 0 0 gen zon reflect dip 90 dd 270 ori 0 0 0 ;网格显示

plot set rotation 30 0 30 plot block group

4-9球体洞穴网格生成

new def parm

rad=4.0 len=10.0 in_size=6 rad_size=10 end parm

def make_sphere

p_gp=gp_head

loop while p_gp#null

px=gp_xpos(p_gp) py=gp_ypos(p_gp) pz=gp_zpos(p_gp)

dist=sqrt(px*px+py*py+pz*pz) if dist>0 then

k=rad/dist ax=px*k ay=py*k az=pz*k

maxp=max(px,max(py,pz)) k=len/maxp bx=px*k by=py*k bz=pz*k

u=(maxp-rad)/(len-rad)

gp_xpos(p_gp)=ax+u*(bx-ax) gp_ypos(p_gp)=ay+u*(by-ay) gp_zpos(p_gp)=az+u*(bz-az) end_if

p_gp=gp_next(p_gp) end_loop end

gen zone radbrick edge len size in_size in_size in_size rad_size &

rat 1.0 1.0 1.0 1.2 dim rad rad rad make_sphere

;对称得出完整模型

gen zon reflect dip 0 dd 0 ori 0 0 0 gen zon reflect dip 90 dd 270 ori 0 0 0 ;网格显示

plot set rotation 15 0 30 plot block group

4-10应力边界

new

gen zone brick size (4,4,4) p0 (0,0,0) p1 (4,0,0) & p2 (0,4,0) p3 (2,0,3.464) model elastic

prop bulk 1e8 shear .3e8

apply nstress -1e6 range plane dip 60 dd 270 origin 0.1,0,0 above step 1 ;网格显示

plot set rotation 30 0 30 plot block group plot grid red fap green

4-11改变应力边界条件

new

gen zone brick size 6 6 6 model elas

prop bulk 1e8 shear 7e7 fix x range x -0.1 0.1 def superstep

loop ns (1,n_steps) x_stress = stress_inc command

apply sxx add x_stress range x 5.9,6.1 y 0,6 z 0,2 step 100

end_command end_loop end

set n_steps=100 stress_inc=-1e3 hist zone sxx 6,0,0 plot create sxx_hist plot add hist 1 plot show Superstep ;网格显示

plot set rotation 30 0 30 plot block group plot grid red fap green new

gen zone brick size 6 6 6 model elas

prop bulk 1e8 shear 7e7 fix x range x -0.1 0.1 def x_stress

x_stress = stress_inc * step end

set stress_inc = -1e3

apply sxx 1.0 hist x_stress range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0 hist x_stress step 100 ;网格显示

plot set rotation 30 0 30 plot grid red fap green new

gen zone brick size 6 6 6 model elas

prop bulk 1e8 shear 7e7 fix x range x -0.1 0.1 table 1 0,0 100,-1e5

apply sxx 1.0 hist table 1 range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0 step 100;网格显示 plot set rotation 30 0 30 plot grid red fap green new

gen zone brick size 6,6,6 p1 6,0,-1 model elas

prop bulk 8e9 shear 5e9

apply sxx -2e6 range x -0.1 0.1 apply sxx -2e6 range x 5.9 6.1 step 500

plot bound fap plot grid disp new

gen zone brick size 5,5,5 model elas

prop bulk 8e9 shear 5e9 set grav 0 0 -10 ini dens 1000

fix x range x -.1 .1 fix x range x 4.9 5.1 fix y range y -.1 .1 fix y range y 4.9 5.1

ini szz -5e4 grad 0 0 -1e4 app szz -5e4 range z -.1 .1 solve

model null range x 1,4 y 1 4 z 3 5 step 100

plot set plane dip 90 dd 180 origin 0,2.5,0 plot add surf plane behind yell plot add vel plane behind black plot set rotation 30 0 30 plot block group plot show

4-12位移边界

new

gen zone brick size 4,4,4 p3 2,0,3.464 model elas

prop bulk 1e8 shear .3e8

macro left_boun ’plane dip 60 dd 270 origin 0.1,0,0 above’ macro right_boun ’plane dip 60 dd 270 origin 3.9,0,0 below’ apply nvel 0.1 plane dip 60 dd 270 range left_boun apply nvel 0.1 plane dip 120 dd 90 range right_boun step 1 ;网格显示

plot set rotation 30 0 30 plot block group plot grid vel new new

gen zone brick size 4 4 4 model elastic

prop bulk 1e8 shear .3e8

apply nstress -1e6 plane dip 0 dd 0 range y 3.9 4.1 apply nvel 0.0 plane dip 0 dd 0 range y -.1 .1

apply nvel 0.0 plane dip 0 dd 0 range x -.1 .1 step 10 ;网格显示

plot set rotation 30 0 30 plot block group plot grid vel new

gen zone brick size 10 5 5 mod el

prop shear 1e8 bulk 2e8

fix x y z range x -.1 .1 y 0 5 z 0 5 fix x y z range x 0 10 y 0 5 z -.1 .1 fix x y z range x 9.9 10.1 y 0 5 z 0 5 table 1

def find_add head = null p_gp = gp_head

loop while p_gp # null x_pos = gp_xpos(p_gp) if x_pos = width then new = get_mem(2) mem(new) = head mem(new+1) = p_gp head = new endif

p_gp = gp_next(p_gp) endloop end

set width=10.0 find_add def apply_vel while_stepping ad = head

loop while ad # null p_gp = mem(ad+1)

gp_xvel(p_gp) = vel_max * gp_zpos(p_gp) / height

gp_zvel(p_gp) = -vel_max * (gp_xpos(p_gp) - width) / height ad = mem(ad) endloop end

set large vel_max=1e-2 height=5.0 step 100 ;网格显示

plot set rotation 30 0 30 plot block group plot grid vel

4-13不考虑重力影响的均匀应力

new ;set grav 0 0 0

gen zone brick size 6 6 6 model elas

ini sxx=-5e6 syy=-1e7 szz=-2e7 apply sxx=-5e6 range x -0.1 0.1 apply sxx=-5e6 range x 5.9 6.1 apply syy=-1e7 range y -0.1 0.1 apply syy=-1e7 range y 5.9 6.1 apply szz=-2e7 range z -0.1 0.1 apply szz=-2e7 range z 5.9 6.1 prop bulk 8e9 shear 5e9 ini dens 2000 fix x range z -.1 .1 step 10

plot set rotation 30 0 30 ;网格三维显示 plot cont zdisp outl on

4-14考虑应力梯度的均匀材料

new

gen zone brick size 10 10 10 p1 20,0,0 p2 0,20,0 p3 0,0,20 model mohr

prop bulk 5e9 shear 3e9 fric 35 ini density 2500 set gravity 0,0,-10

fix x y z range z -0.1 0.1 ini szz -5.0e6 grad 0,0,2.5e4 ini syy -2.5e6 grad 0,0,1.25e4 ini sxx -2.5e6 grad 0,0,1.25e4 apply szz -4.5e6 range z 19.9 20.1 apply szz -5.0e6 range z -0.1 0.1

apply sxx -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1 apply sxx -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1 apply syy -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1 apply syy -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1 step 10

;网格三维显示

plot set rotation 30 0 30

plot cont zdisp outl on

4-15考虑应力梯度的非均匀材料

new

gen zone brick size 10 10 10 &

p0 0,-25,0 p1 20,-25,0 p2 0,0,0 p3 0,-25,20 model elas

prop bulk 5e9 shear 3e9

ini density 1600 range y -10,0 ini density 2000 range y -15,-10 ini density 2200 range y -25,-15 set gravity 0,-10,0 fix x range x -.1 .1 fix x range x 19.9 20.1 fix z range z -.1 .1 fix z range z 19.9 20.1 fix y range y -25.1 -24.9

ini syy 0.0 grad 0,1.6e4,0 range y -10,0 ini syy 4e4 grad 0,2.0e4,0 range y -15,-10 ini syy 7e4 grad 0,2.2e4,0 range y -15,-25 step 10

;网格三维显示

plot set rotation 30 0 30 plot cont zdisp outl on

4-16非均匀网格应力初始化

new

gen zone radcyl size 3 8 4 5 fill p1 10,0,0 p2 0,10,0 p3 0,0,10 mode elastic

prop shear 3e8 bulk 5e8 fix x range x -.1 .1 fix x range x 9.9 10.1 fix y range y -.1 .1 fix y range y 9.9 10.1 fix z range z -.1 .1

ini szz = -2.5e5 grad 0,0,2.5e4 ini density 2500 set grav 0,0,-10 step 10

;网格三维显示

plot set rotation 30 0 30 plot cont zdisp outl on

;围岩的密度

ini density 1800

;体积、剪切、摩擦角、凝聚力、抗拉强度

prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 step 3000 save Gsol.sav ;第1步 rest Gsol.sav

initial xdisp=0 ydisp=0 zdisp=0 model null range group inrz y 0 6 model elas range group linz y 0 6

ini density 2500 range group linz y 0 6

prop bulk=16.67e9,shear=12.5e9 range group linz y 0 6 step 3000

save step1.sav ;第2步

rest step1.sav

model null range group inrz y 6 12 model elas range group linz y 6 12

ini density 2500 range group linz y 6 12

prop bulk=16.67e9,shear=12.5e9 range group linz y 6 12 step 3000

save step2.sav ;第3步

rest step2.sav

model null range group inrz y 12 18 model elas range group linz y 12 18

ini density 2500 range group linz y 12 18

prop bulk=16.67e9,shear=12.5e9 range group linz y 12 18 step 3000

save step3.sav ;第4步

rest step3.sav

model null range group inrz y 18 24 model elas range group linz y 18 24

ini density 2500 range group linz y 18 24

prop bulk=16.67e9,shear=12.5e9 range group linz y 18 24 step 3000

save step4.sav ;第5步

rest step4.sav

model null range group inrz y 24 30 model elas range group linz y 24 30

ini density 2500 range group linz y 24 30

prop bulk=16.67e9,shear=12.5e9 range group linz y 24 30 step 3000

save step5.sav ;第6步

rest step5.sav

model null range group inrz y 30 36 model elas range group linz y 30 36

ini density 2500 range group linz y 30 36

prop bulk=16.67e9,shear=12.5e9 range group linz y 30 36 step 3000

save step6.sav ;第1步

rest step6.sav

model null range group inry y 0 6 model elas range group liny y 0 6

ini density 2500 range group liny y 0 6

prop bulk=16.67e9,shear=12.5e9 range group liny y 0 6 step 3000

save step7.sav ;第2步

rest step7.sav

model null range group inry y 6 12 model elas range group liny y 6 12

ini density 2500 range group liny y 6 12

prop bulk=16.67e9,shear=12.5e9 range group liny y 6 12 step 3000

save step8.sav ;第3步

rest step8.sav

model null range group inry y 12 18 model elas range group liny y 12 18

ini density 2500 range group liny y 12 18

prop bulk=16.67e9,shear=12.5e9 range group liny y 12 18 step 3000

save step9.sav ;第4步

rest step9.sav

model null range group inry y 18 24 model elas range group liny y 18 24

ini density 2500 range group liny y 18 24

prop bulk=16.67e9,shear=12.5e9 range group liny y 18 24 step 3000

save step10.sav ;第5步

rest step10.sav

model null range group inry y 24 30 model elas range group liny y 24 30

ini density 2500 range group liny y 24 30

prop bulk=16.67e9,shear=12.5e9 range group liny y 24 30 step 3000

save step11.sav ;第6步

rest step11.sav

model null range group inry y 30 36 model elas range group liny y 30 36

ini density 2500 range group liny y 30 36

prop bulk=16.67e9,shear=12.5e9 range group liny y 30 36 step 3000

save step12.sav new

4-17不规则自由面应力初始化

gen zone brick size 15 15 10 p0 0,0,0 edge=100.0 model elastic

prop shear 3e8 bulk 5e8 def mountain gp = gp_head

loop while gp # null

zz = sqrt(gp_xpos(gp)?2 + gp_ypos(gp)?2)

dz = 0.06 * sin(0.2 * zz + 100.0) ; Sum Fourier terms for dz = dz + 0.06 * sin(0.22 * zz - 20.3) ; quasi-random surface dz = dz - 0.04 * sin(0.33 * zz + 33.3) ; topology. gp_zpos(gp) = 0.5 * gp_zpos(gp) * (1.0 + dz) gp = gp_next(gp) end_loop end

mountain

fix x range x -.1 .1

fix x range x 99.9 100.1 fix y range y -.1 .1

fix y range y 99.9 100.1 fix z range z -.1 .1 set grav 0,0,-10 ini density=2000

ini szz=-2.0e6 (grad 0,0,2.0e4) sxx=-4.0e6 (grad 0,0,4.0e4) syy=-4.0e6 (grad 0,0,4.0e4) step 100

plot set rotation 30 0 30

plot cont zdisp outl on;网格三维显示

4-18非均网格内部压实

new

gen zone brick size 8 8 10 ratio 1.2 1 1 model mohr ini dens 2000

prop bulk 2e8 shear 1e8 prop fric 30

fix x range x -.1 .1 fix x range x 7.9 8.1 fix y range y -.1 .1 fix y range y 7.9 8.1 fix z range z -.1 .1

;ini szz -2.0e5 grad 0,0,2e4 ;ini sxx -1.5e5 grad 0,0,1.5e4 ;ini syy -1.5e5 grad 0,0,1.5e4

set grav 10 step 1000 ;pause

;prop tens 1e10 coh 1e10 ;step 750

;prop tens 0 coh 0 ;step 250

;网格三维显示

plot set rotation 30 0 30 plot cont zdisp outl on

4-19模型改变后初始应力变化

new

gen zone brick size 5 5 5 model elas

prop sh 2e8 bu 3e8 fix x y z range z -.1 .1 set grav 0 0 -10 ini dens 2000 Solve

;网格三维显示

plot set rotation 30 0 30 plot cont zdisp outl on ;模型改变

model mohr range x 0 2 y 0 5 z 0 2

prop sh 2e8 bu 3e8 fric 35 range x 0 2 y 0 5 z 0 2 Step 2000

;网格三维显示

plot set rotation 30 0 30 plot cont zdisp outl on

4-20应力与孔隙压力的初始化

new

config fluid

gen zone brick size 8 5 10 model elas model fl_iso ini dens 2000

prop bulk 1e9 shear 5e8 prop poros 0.5 perm 1e-10 ini fmod 2e9 ini fdensity 1e3 ini sat 0

ini sat 1 range z -.1 5.1

set grav 0 0 -10 fix x range x -.1 .1 fix x range x 7.9 8.1 fix y range y -.1 .1 fix y range y 4.9 5.1 fix z range z -.1 .1

ini pp 5.e4 grad 0,0,-1.e4 range z 0.0 5. ini szz -20e4 grad 0,0,20e3

ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5 solve

;网格三维显示

plot set rotation 30 0 30 plot cont zdisp outl on

new

gen zone brick size 1 1 10 model elas

ini dens 2500 range z 0 5 ini dens 2250 range z 5 6 ini dens 2000 range z 6 10 prop bulk 1e9 shear 5e8 set grav 10

water dens 1000

water table ori 0 0 5 normal 0 0 1 fix x range x -.1 .1 fix x range x 7.9 8.1 fix y range y -.1 .1 fix y range y 4.9 5.1 fix z range z -.1 .1

ini szz -20e4 grad 0,0,20e3

ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5 solve

;网格三维显示

plot set rotation 30 0 30 plot cont zdisp outl on

4-21加载顺序

new

; 主隧道建模

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

;设置重力加速度 set gravity 0 0 -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 ;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 xdisp plot cont szz plot cont sxx ;new

;1-上半断面开挖 rest Gsol.sav

ini xdis 0.0 ydis 0.0 zdis 0.0

model null range group insiderocksz any group concretlinersz any & group insiderocksy any group concretlinersy any step 2000

save stepp1.sav plot cont zdisp ;2-初期支护 rest stepp1.sav

sel shell id=1 range x -5.56 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.20 step 20

save stepp2.sav ;3-下半断面开挖 rest stepp2.sav

model null range group insiderockxz any group concretlinerxz any & group insiderockxy any group concretlinerxy any step 2000

save stepp3.sav plot cont zdisp ;4-二次衬砌

rest stepp3.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

save stepp4.sav

6.3 III级围岩施工过程模拟

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.25 5.25 5.25 5.25 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.7 5.25 3.7 5.25 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

;设置重力加速度 set gravity 0 0 -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 ;0-求解自重应力场 model mohr ini density 2400 ;围岩的密度

prop bulk=6.0e9, shear=10.0e9 fric=45, coh=1.2e6 tension=2e6 ;体积、剪切、摩擦角、凝聚力、抗拉强度 step 6000

save Gsol.sav ;1-全断面开挖 rest Gsol.sav

ini xdis 0.0 ydis 0.0 zdis 0.0

model null range group insiderocksz any group concretlinersz any & group insiderocksy any group concretlinersy any model null range group insiderockxz any group concretlinerxz any & group insiderockxy any group concretlinerxy any step 2000

save stepp1.sav ;2-初期支护 rest stepp1.sav

sel shell id=1 range group x -5.56 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.15 step 2000

save stepp2.sav ;3-二次衬砌 rest stepp2.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

save stepp3.sav plot block group new

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

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 plot cont zdisp plot cont szz new

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

new

gen zon radcyl p0 9 0 0 p1 18 0 0 p2 9 36 0 p3 9 0 9 &

size 4 24 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outr gen zone cshell p0 9 0 0 p1 12 0 0 p2 9 36 0 p3 9 0 3 &

size 1 24 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group liny fill group inr gen zon reflect dip 90 dd 90 orig 9 0 0 gen zon reflect dip 0 dd 0 ori 9 0 0 ;1/2隧道地层部分

gen zon brick p0 18 0 -9 p1 36 0 -9 p2 18 36 -9 p3 18 0 9 size 10 12 8 group outr1

gen zon brick p0 18 0 9 p1 36 0 9 p2 18 36 9 p3 18 0 13 size 10 12 2 group outr2

gen zon brick p0 0 0 9 p1 18 0 9 p2 0 36 9 p3 0 0 13 size 8 12 2 group outr2 gen zon brick p0 18 0 -23 p1 36 0 -23 p2 18 36 -23 p3 18 0 -9 size 10 12 8 group outr3

gen zon brick p0 0 0 -23 p1 18 0 -23 p2 0 36 -23 p3 0 0 -9 size 8 12 8 group outr3

gen zon reflect dip 90 dd 270 orig 0 0 0 group linz range group liny x -12.1 -5.99 ;绘制模型图 plot add axes red

plot set rotation 20 0 45 plot block group ;0-自重应力场求解 set gravity 0 0 -10 ;给定边界条件

fix z range z -23.01,-22.99 fix x range x -36.01,-35.99 fix x range x 35.99,36.01 fix y range y -0.01 0.01 fix y range y 35.99,36.01 model mohr

本文来源：https://www.bwwdw.com/article/ws9v.html