FLAC3D3.0在某隧道工程开挖支护中的应用 U7g,@/Qx @>r3=s.Q
隧道建模命令流入下： ieWXr4@: Oy,7>vWQI
set log on (|<+yQ,@> @O!BQ^'hk#
set logfile yang.log +eKLwM *{t]fds
gen zon radcyl p0 0 0 0 p1 9.0 0 0 p2 0 50 0 p3 0 0 8 & /F|VYl^_ eS qKXmH[m
size 4 20 6 4 dim 6 5 6 5 rat 1 1 1 1 group 围岩 L\bc R Bb,l.w
gen zon cshell p0 0 0 0 p1 6.0 0 0 p2 0 50 0 p3 0 0 5.0 & 3KkJQ5a W E-cq1)
size 4 20 6 4 dim 5.6 4.6 5.6 4.6 rat 1 1 1 1 group 初期支护 =1)9>=} N)kZ2|oD
gen zon cshell p0 0 0 0 p1 5.6 0 0 p2 0 50 0 p3 0 0 4.6 & <1x u&Z7 <}\!FuC
size 4 20 6 4 dim 5.0 4.0 5.0 4.0 rat 1 1 1 1 group 二次衬砌 fill group 原岩 vJ 28A t",=]k
gen zon radcyl p0 0 0 0 p1 0 0 -8.0 p2 0 50 0 p3 9.0 0 0 & **1=|aa: :>:F6Db"U
size 4 20 6 4 dim 3 6 3 6 rat 1 1 1 1 group 围岩2 Y(T$k9%}+ I%^Bl:M
gen zon cshell p0 0 0 0 p1 0 0 -3.0 p2 0 50 0 p3 6.0 0 0 & M A%g-} `b11,lg
size 4 20 6 4 dim 2.6 5.6 2.6 5.6 rat 1 1 1 1 group 仰拱初期支护 !g|[A7<| Qj1q x;S
gen zon cshell p0 0 0 0 p1 0 0 -2.6 p2 0 50 0 p3 5.6 0 0 & |GPY bxzc )p$a1\ ~m
size 4 20 6 4 dim 2 5 2 5 rat 1 1 1 1 group 仰拱二次衬砌 fill group 仰拱原岩 ,#0#1k<Dm HX; JO[0
gen zone reflect normal -1 0 0 FS!vnl8` :i?7RouO
gen zone radtun p0 0 0 0 p1 45 0 0 p2 0 50 0 p3 0 0 20 & N>sH T =_ ,P%a0\
size 3 20 3 12 dim 9 8 9 8 rat 1 1 1 1.1 group 围岩3 D&],.N GOKca%DT=
gen zon reflect dip 0 ori 0 0 0 range x 0 9 y 0 50 z 8 20 3YUF\L]yyw X
j>?P/=Z
gen zon reflect dip 0 ori 0 0 0 range x 9 45 y 0 50 z 0 20 _tQM<~Y]u\ }{R?i,j(
gen zon reflect dip 90 dd 270 ori 0 0 0 range x 0 9 y 0 50 z 8 20 aHitPPlq c Y C@@?
gen zon reflect dip 90 dd 270 ori 0 0 0 range x 0 9 y 0 50 z -8 -20 .i#'IS0c U<YP@?w
gen zon reflect dip 90 dd 270 ori 0 0 0 range x 9 45 y 0 50 z -20 20 <Y;w I#C ^^9O9]
gen zon brick p0 -45 0 -20 p1 -45 0 -40 p2 -45 50 -20 p3 45 0 -20 & WK{F wM-I*<L>
size 5 20 6 rat 1.1 1 1 group 围岩4 f&Juq8s_0 vKcZgIR
save tun_model.sav PT4iy< c:J;Q){Xz
假设围岩岩体符合mohr-coulomb本构模型，给围岩赋参数命令流如下， lpbcpB ]4uY<9VL
; mohr-coulomb model sv6 m)pwh Y<]A 5cm
model mohr @)ls+}=Y -4%{Jb- 1
def derive @qB>qD~WsD 0^0Q0A
s_mod1=E_mod1/(2.0*(1.0+p_ratio1)) I*n]8c {Rm N1'%
b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1)) tY#Zl 54~{ ^$N}[1
s_mod2=E_mod2/(2.0*(1.0+p_ratio2)) E,Xl8rC hT=6XO od4
b_mod2=E_mod2/(3.0*(1.0-2.0*p_ratio2)) XQ8q)B= "^Ax}Jr
end 3jHE,5m Smt&/~7D%
set E_mod1=0.6e9 p_ratio1=0.27 E_mod2=0.8e9 p_ratio2=0.26 qI<c47d;q d)R7#HLZ7
derive #8y"1I=i& L@|xpq
prop bulk b_mod1 shear s_mod1 cohe 1.8e6 tens 0.8e6 fric 30 range z 4.5 20 tqjjn5! 8\N`2mPt
prop bulk b_mod2 shear s_mod2 cohe 2.8e6 tens 1.0e6 fric 35 range z -40 4.5 e7plL^^` lpy ( un
ini dens=2300 D$E#:[ Vr<eU> W
set grav 0 0 -10 ]@}o"Td u l-A'
; boundary and initial conditions _nSEp >]L (GeOD V?U
apply szz -1.4e6 range z 19.9 20.1 ubc k{\. z]twh&^1L
fix z range z -40.1 -39.1 ?yu@eo JO3x#1~;_
fix x range x -45.1 -44.9 b`f6(6 33M10
1X{6
fix x range x 44.9 45.1 :z6? fn~Jc~[G|
fix y range y 49.9 50.1 (!Xb8rV0_ !_No\ O
hist unbal r$Ik* R aqw;T\GI+~
hist gp xdis 6.0,0,0 5r d t xJ);P.
hist gp zdis 0,0,5 muKu@nshL @@ 1Sxv_
hist gp xdis 6.0,50,0 2 #cw_Ua 5x2m ]u
hist gp zdis 0,50,5 K1- RJj\L 6EX_IDb
plot hist 3 ZVit] 3hd Knjg`f
solve ?9?o8! RT`jWWh*Lo
save tun_nature.sav X7gB.=\X [z2jR(+`U
对后面计算而言，模型建立时岩体在开挖前认为位移已经终了，因此需要对位移进行“清零”，而应力可以保留。假设隧道先开挖上断面，中间不设支护，直到进尺50m，那么此时位移和应力的分布情况可用如下命令流 $"#M:V @ -p-<mC@<&S
ini xdis=0 ydis=0 zdis=0 !~ j9Oc^ h1(GzL%i_
plot re y{XNB}E RBA{!
model null range group 原岩 M3dUGM fz[-pJ5[
model null range group 二次衬砌 5r;)P po \# hp,XV>
model null range group 初期支护 =w:H9uj6F f\|R<3 L
set large ZY;g)`E1 fF|m~#y
hist unbal h,]lN'JG{ G-DvM6T
plot hist 3 R80|q#h,] 7 { n>0@_
solve ~u?rjkSFoh Rxf.@E
save tun_ext1.sav 472 'P 6^+T_{gl
如果在开挖后适时对隧道进行锚喷混凝土初期支护，该命令流如下： 37V$Qb_ (6Y.|u]bq
restore tun_nature.sav Exat_ L'? !eC]=PoY
ini xdis=0 ydis=0 zdis=0 g](&H$g 3d>3f3D8;
plot re U/7jK40 -~|E (ys
model null range group 初期支护 >PGsY[N Ec/-f `8
model null range group 二次衬砌 V0gu0+u~R Ih;I&D+e;
model null range group 原岩 / ]nrxT smJ#.I6/L
step 100 (k%r_O 6 <5xlP:Cx
hist unbal 4&~*;an7 ~ Pm[ Ud
plot hist 3 rE\&FVx @hG]Gs[,o
mo el range group 初期支护 `oE.$~' ;bMmJ>[l-
prop bulk 1.33e9 shear 0.8e9 cohe 2.2e6 tens 1.2e6 fric 32 range z 4.5 7.0 x -7.0 7.0 Ge76/T%{Q UNC%<=
prop bulk 10.9e9 shear 8.9e9 dens 2500 range group 初期支护 M/}i7oS] $q%l)]+
set large QX9['B< `J26Y"]P
solve _%aT3C}k ssx #\
save tun_ext2.sav
