FLAC3D计算隧道作业..docxVIP

在某Ⅳ级围岩中开挖一半圆拱直墙形隧道，隧道跨度10m，边墙高5m，隧道埋深500m，假设围岩为理想弹塑性材料，请采用有限元或有限差分方法分析以下问题： （1）自重应力场作用下隧道开挖后的拱顶下沉和边墙水平收敛大小以及围岩中的塑性区大小。 （2）若侧压系数为0.5—2.5，请分析构造应力对隧道拱顶下沉、边墙水平收敛大小以及塑性区的影响。（3）若开挖后采用锚喷支护，在隧道拱部和边墙布设系统锚杆，锚杆为全长锚固的金属锚杆，垂直于洞壁布设，间距1.5m，长度3.0m，直径25mm。喷射混凝土厚度100mm，标号为C20，请分析支护效果。本题采用FLAC3D软件建模计算分析 隧道未开挖时的立体模型 隧道开挖后的立体模型第一步，建模由于隧道的半径为5m，根据经验取6倍的隧道半径为围岩影响区，所以取30m为边界，划分网格的边长为0.5m ，本题只分析x-z平面上的受力及位移情况即可，建模命令流如下：new ;建立模型gen zone radcyl p0 0 0 0 p1 30 0 0 p2 0 1 0 p3 0 0 30 size 10 1 10 30 dim 5 5 5 5 ratio 1 1 1 1 group outsiderock ;右上圆形部分围岩gen zone cshell p0 0 0 0 p1 5 0 0 p2 0 1 0 p3 0 0 5 size 1 1 10 10 dim 4.9 4.9 4.9 4.9 rat 1 1 1 1 group concretliner fill group insiderock ;右上半圆衬砌gen zone radtun p0 0 0 0 p1 0 0 -30 p2 0 1 0 p3 30 0 0 size 10 1 10 30 dim 5 5 5 5 ratio 1 1 1 1 group outsiderock ;右下矩形部分围岩gen zone radtun p0 0 0 0 p1 0 0 -5 p2 0 1 0 p3 5 0 0 size 10 1 10 1 dim 4.9 4.9 4.9 4.9 ratio 1 1 1 1 group concretliner ;右下矩形部分衬砌gen zone brick p0 0 0 -4.9 p1 add 4.9 0 0 p2 add 0 1 0 p3 add 0 0 4.9 size 10 1 10 ratio 1 1 1 1 group insiderock ;隧道内部gen zon reflect dip 90 dd 90 orig 0 0 0 ;关于z轴对称plot block groupplot add axes blacksave jianmo.sav（1）自重应力场作用下隧道开挖后的拱顶下沉和边墙水平收敛大小以及围岩中的塑性区大小。围岩选用M-C模型，参数的选取：岩体密度（kg/m3）= 2200, 重力加速度（kg/N）=10, 体积模量(Pa)=5e10，剪切模量(Pa)= 1.1e10，摩擦角=30o，粘聚力（N）= 1.8e7，抗拉强度（N/m2）=1e7。在Z=30m平面处施加的应力为γh=22000×(500-30)=10.34 e6N/m2。命令流如下；;边界条件fix z range z -30.01 -29.99 fix x range x -30.01 -29.99 fix x range x 29.99 30.01 fix y range y -0.01 0.01 fix y range y 0.99 1.01 apply szz -10.34e6 range z 29.99 30.01 ;设置初始应力状态model mohr set gravity 0 0 -10 ;设置重力ini density 2200 ;围岩的密度prop bulk=5e10 shear=1.1e10 fric=30 coh=0.18e7 ten=1e7 ;体积模量、剪切模量、摩擦角、粘聚力、抗拉强度set mech ratio=1e-4 plot block statesolve ;求解自重应力场save Gravsol.sav plot block group plot add axes red ;毛洞开挖计算initial xdisp=0 ydisp=0 zdisp=0 model null range group insiderock any group concretliner anyset mech ratio=1e-4hist unbhist gp zdisp 0 0 5 hist gp xdisp 5 0 -2plot block state solve ;save maod