基于有限元法的浅埋无压隧洞稳定性分析.doc

基于有限元法的浅埋无压隧洞稳定性分析 戴跃华[1] ，李正良[2],梁蔚[3] （1.广东省水利电力规划勘测设计研究院，广东 广州；2.湖南省水利水电勘测设计研究总院，湖南 长沙；3. 水利部长江水利委员会长江科学院，湖北 武汉） 摘 要：浅埋隧洞由于上部覆盖层厚度小，顶拱自承能力不足，且地表处围岩地质条件往往较差，隧洞投入运行后，对外部荷载的扰动较为敏感，严重的将危及隧洞的安全运行。本文以某无压供水隧洞为例，采用有限元法，分析铁路修建后隧洞围岩的稳定性及衬砌结构的安全性。计算分析表明，在现有覆盖厚度下，修建铁路所新增的荷载对隧洞围岩及衬砌结构的影响是有限的，隧洞的安全稳定性是有保证的。 关 键 词：浅埋隧洞 围岩 衬砌 稳定 中图分类号：TU473.1 文献标识码：B 1前言 浅埋隧洞由于上部覆盖层厚度小，顶拱自承能力不足，且地表处围岩地质条件往往较差，隧洞投入运行后，对外部荷载的扰动较为敏感，严重的将危及隧洞的安全运行 某供水隧洞设计横断面为圆拱直墙式，净宽5.6m，净高7m，隧洞沿线沟谷发育，洞顶上的覆盖层厚度多在20m~30m。本次规划修建的铁路以平面20.7°的夹角经过隧洞正上方，隧洞与铁路在平面上形成交叉段（以下简称交叉段），隧洞顶距铁路路基约30m。 铁路路基以路堤形式通过，路堤依地形走势填高约3m～15m，采用分层填筑碾压施工，设计将基床范围内厚约3m的人工填土挖除并换填，对以下厚约10m的人工填土采用水泥粉煤灰碎石桩进行加固。铁路施工及投入运行后，交叉段隧洞将承受新增的路基填筑土及列车荷载，对隧洞稳定和结构安全构成一定的影响。 根据铁路勘探资料揭示，交叉段地层从上至下主要为：表层为人工素填土，成分主要为全风化～强风化的凝灰质砂岩与页岩，以及粉质粘土、耕植土的混合堆积，局部夹深灰色的凝灰质砂岩碎块，结构较松散，厚约14m；下部为粉质黏土，浅黄～褐红色，夹少量碎块石，硬塑，厚2m～3m；下伏基岩为白垩系凝灰质砂岩，局部夹泥质页岩，全风化～弱风化，其中全风化呈土状、砂土状，黄褐色、局部粉红色，厚约15m。 本文采用有限元法对交叉段隧洞围岩的稳定性及衬砌结构的安全性进行分析[1]，计算程序采用通用有限元计算软件ANSYS。 2 计算原理 根据工程实际情况，对隧洞开挖、衬砌施工及铁路荷载的施加按以下步骤进行模拟[2]：计算初始地应力场（自重应力场）→隧洞开挖→隧洞衬砌施工（考虑应力释放）→施加内、外水压力→铁路路基回填+施加铁路荷载。 (1) 初始地应力场的模拟 考虑到隧洞埋深较浅，地应力场的计算只考虑自重应力[3]，忽略构造应力。ANSYS程序可以通过输入初始地应力文件的方法将初始地应力作为一种荷载进行存储，将初始地应力场下的变形进行清零，在开挖条件下计算得到的变形是由开挖产生的净变形。 (2) 荷载释放系数 岩体开挖过程是开挖岩体部分的应力释放过程[4]，实际工程中是等开挖部分已经释放了一定的荷载，即开挖自由边界上的位移状态已经有了一定的改变，才进行临时支护，这个释放荷载用应力释放系数表示，参考其他类似工程的经验，并依据隧洞的实际开挖支护过程，应力释放系数取为0.5。 (3) 隧洞开挖、衬砌支护及路基回填的模拟 在ANSYS程序中，隧洞的开挖、喷砼支护及衬砌采用单元的“生”、“死”功能进行模拟，即通过将隧洞开挖单元“杀死”来模拟开挖，通过激活喷砼单元及衬砌单元模拟喷砼及衬砌的支护作用；铁路路基的单元在回填之前被“杀死”，当路基回填以后，通过激活该部分单元模拟路基的回填。 (4) 弹塑性有限元计算原理 由于交叉段隧洞主要以Ⅳ、Ⅴ类围岩为主，因此计算时须考虑材料的非线性，进行弹塑性有限元计算[5]，弹塑性有限元计算原理如下： 由虚功原理，单元平衡方程可表示为 （2-1） 式中{F}e为单元上的等效结点力；{(}e为单元结点位移列阵；[k]为单元刚度矩阵，由下式计算 （2-2） [D]ep为弹塑性矩阵，按照塑性势理论，有 （2-3） 其中[De]为弹性矩阵；[Dp]为塑性矩阵；A为与硬化有关的参数；F为塑性屈服函数；Q为塑性势函数，当F=Q时，称为关联流动法则。对于理想弹塑性材料，取A=0。 屈服函数采用Drucker-Prager准则。其表达式为： （2-4） 式中I1为应力张量第一不变量；J2为应力偏量第二不变量； （2-5） 其中c为岩土类摩擦型材料的粘聚力，φ为内摩擦角。 3 计算条件及模型 3.1 材料参数 交叉段地层材料参数详见表3-1；隧洞采用全断面混凝土衬砌（厚600mm），混凝土材料参数详见表3-2。 3.2 计算荷载 计算荷载主要包括：自重地应力场、围岩压力、内水压力、外水压力、路路基填筑土自重及列车荷载。 其中，内水压力，运行期洞内水深取5.47m；考虑到隧