压缩高地应力单侧软岩大变形隧道施工数值模拟分析.docVIP

文章编号： 高地应力单侧软岩大变形隧道施工数值模拟分析 作者 （单位） 摘要：以岩层倾角对高地应力单侧软岩大变形隧道位移场、塑性区和应力场的影响为研究目的。结合宜巴高速公路石门垭隧道施工过程，利用有限差分软件FLAC3D对隧道施工进行数值模拟分析，对比分析不同角度下全断面法和下导洞超前开挖法模拟计算结果。结果表明：软弱围岩区位移场、塑性区和应力场分布特征受岩层倾角影响较大，坚硬围岩区受其影响较小。两种开挖方式下隧道周边位移和塑性区分布差别不大；最终应力场分布特征与开挖方式无关，但开挖期间主应力、主应变、位移和破坏度并不相同。所得结论可为同类隧道的设计、施工和研究提供借鉴和参考。 关键词：岩层产状；软岩大变形；FLAC3D；数值模拟；下导洞超前开挖法 1 前言 随着我国交通事业的迅速发展，在深部岩体中修筑隧道工程已必不可少，随之而来的深部岩体所具有的特殊工程地质问题也更加突出。在深埋隧道勘察设计和施工过程中，高地应力的存在，是影响隧道稳定的重要因素，主要表现为硬岩岩爆和软岩大变形或塑性破坏 2 工程概况 石门垭隧道是宜巴高速公路控制性工程之一，为项目全线最长隧道。隧道采用分幅式，左幅起讫桩号ZK118+963~ZK126+487，总长7524.0m，右幅起讫桩号YK118+948~YK126+441总长7493.0m。隧道最大埋深约878m，属特长深埋隧道。 隧道施工至ZK120+000时，3 隧道施工数值模拟 3.1 计算模型 采用FLAC3D三维快速拉格朗日差分方法分析软件，根据地下结构的计算原理，隧道开挖影响范围为洞径的3~5倍x轴）长度取90m，竖直方向（y轴）取90m，纵向（z轴）沿隧道轴线方向取60m。围岩材料模型采用Mohr-Coulomb理想弹塑性模型开挖采用FLAC3D中的Null模型确定，为上覆岩层的平均加权容重，为上覆岩层总厚度[8-10]。计算模型如图3所示。 图 隧道计算模型 Fig.3 Computation model of tunnel 围岩 γ (kN/m3) E/GPa u φ/(°) C/MPa 1 Ⅳ级砂岩 23 4.5 0.32 32 0.6 2 Ⅲ级泥岩 24 9.9 0.3 40 0.8 因计算模型中岩性不唯一，将计算模型划分为三个区域，参数赋值分区、如下，其中角度θ为岩层倾斜线与x轴负方向所成的夹角： （1）当时，该区域内围岩为Ⅳ级砂岩，物理力学参数见表1； （2）当时，该区域内为岩性变化区，通过FLAC3D中的单元遍历程序对该区域进行赋值，物理力学参数由公式（1）~（5）计算得出； （1） （2） （3） （4） （5） 式中：为材料的粘聚力； 为材料的内摩擦角； 为材料的体积模量； 为材料的剪切模量； 为点到过点边界线的距离。 （3）当时，该区域内围岩为Ⅲ级泥岩，物理力学参数见表1。 图4所示为θ=60°和θ=120°时模型赋值块状示意图。 θ=60° θ=120° 图4 模型赋值块状示意图 Fig.4 Block diagram of model prop 3.3 开挖方案 针对无支护时全断面法、下导洞超前开挖法进行数值模拟，分析因隧道开挖引起的隧道围岩的应力、应变、位移情况及其随角度θ的变化趋势，判断围岩的稳定性。为模拟实际施工开挖顺序，采用如下开挖方案：全断面法采用15步开挖方案；下导洞超前开挖法采用下导洞先行开挖1515步的开挖方案。数值模拟中，隧道沿着z轴方向以3 m长度为一个开挖循环进行模拟[11-12]。 4 全断面法数值模拟结果 4.1 位移场特征 隧道水平收敛和拱顶下沉量测是监控量测的主要内容之一，也是隧道围岩应力状态变化的最直观反应[13]。隧道拱顶竖向位移、拱腰水平位移、拱肩竖向位移及水平位移计算值如表2所示。 表 2 隧道周边位移计算值 Table 2 Calculated displacement of surrounding rock (cm) θ /(°) 拱顶 拱腰 拱腰 拱肩 拱肩 Uy Ux Ux Ux Uy Ux Uy 0 -3.37 3.60 -3.60 2.57 -2.35 -2.57 -2.35 30 -3.55 4.46 -2.93 3.19 -2.68 -2.20 -2.22 45 -3.78 4.87 -2.71 3.53 -2.92 -2.12 -2.23 60 -4.10 5.19 -2.57 3.87 -3.22 -2.11 -2.29 90 -4.96 5.53 -2.47 4.49 -3.90 -2.25 -2.58 120