数值模拟方法在隧道围岩稳定性分析中应用.docVIP

数值模拟方法在隧道围岩稳定性分析中应用 　　摘要： 随着计算机技术的发展，数值模拟方法在工程建设领域得到了广泛的应用。本文以火石岗隧道为依托，在工程施工全过程监控量测的基础上，采用数值模拟方法，对浅埋隧道洞口段围岩稳定性进行研究，结果表明隧道开挖后，围岩向临空面的位移以竖向沉降为主，挤压两侧围岩向着负收敛方向变形发展，支护结构施作后，支护反力使得围岩应力和位移变形速率减小，逐渐趋于平稳，同时隧道开挖对目标断面应力影响的范围在11米左右，对围岩的变形影响的范围在6米左右。 　　关键词：数值模拟，浅埋隧道，围岩稳定，拱顶沉降，周边收敛 　　中图分类号：TU457 文献标志码：A 　　1引言 　　 近年来，数值分析在隧道工程领域的应用越来越广泛，成为隧道工程研究设计的重要手段。数值模拟分析具有很多优点，主要有：①可以模拟复杂的地质条件、复杂的工程结构以及复杂的荷载、边界条件;②在隧道工程开挖过程中，如若用数值软件进行模拟的话，就能从应力应变云图、变形矢量图、位移变化曲线图等图中直接明了地观察岩土体变形过程中的应力场、位移场的变化。与现场模型试验相比，数字模拟不用实际的工程材料、工程仪器以及具体的试验方案，而且数值模拟及时方便的调整相关的模拟参变量的大小，也能适时的停止模拟，观察某一阶段的应力应变，总的来说，数值模拟的效果有时甚至要远远好于现场模型试验 [1]。余存鹏以明垭子软岩隧道为依托，通过FLAC3D数值模拟分析了现场施工引起的隧道围岩变形值，根据位移评判依据来评判隧道的稳定性[2];尚岳全等基于流固耦合理论，利用有限差分软件对含有破碎带的隧道围岩在饱水条件下的开挖稳定性进行分析，得到不同倾角的破碎带在开挖前后的渗流场特性、主应力特性和塑性区特性等结果，并在此基础上分析了地下水的存在对隧道围岩稳定性的影响[3];姚军，王国才等基于新奥法的基本原理要求，采用数值分析开展在不同地应力释放条件下围岩稳定性影响的研究，结果表明，地应力释放越大，锚杆承担的荷载越小，围岩的塑性区发展范围越大[4];廖军，杨万霞等采用有限元模拟分析某一公路隧道的施工过程，研究在不同的工况条件下，隧道围岩的稳定性，根据分析结果为隧道施工选择了合理的开挖施工方法[5]。因此，在隧道工程中，通过采用数值模拟方法研究施工过程中围岩的应力、应变和位移变化，进而分析研究隧道施工过程中的稳定性，具有重要的现实意义。 　　2工程概况 　　 火石岗隧道为贵州省境内仁怀至赤水高速公路第6合同段中的一条中长距离分离式隧道。隧道建筑界限净宽10.25m，净高5m。左幅起讫桩号ZK38+273～ZK38+800段，全长527米，最大埋深100米;右幅起讫桩号YK38+300～YK38+827段，全长527米，最大埋深115米。火石岗隧道左幅进口ZK38+267～ZK38+400段，长133米，隧道埋深较浅，约0～50米，覆盖层较薄。如图1所示。 　　 　　 　　 　　 　　 　　 　　 　　 　　 图1 火石岗隧道洞口地形图 　　 本文依托火石岗隧道左线进口段，在工程施工全过程监控量测的基础上，采用数值模拟方法，对浅埋隧道洞口段围岩稳定性进行研究。 　　3围岩稳定性数值分析 　　 隧道工程中的围岩稳定性问题属于三维空间问题，而三维数值模拟能够将时间和空间效应较真实地再现出来。根据隧道实际工程情况，基于FLAC3D数值分析软件，建立三维模型，模拟分析隧道施工过程中围岩的应力应变，并将分析结果与监控量测结果相对比，相互验证，为以后的支护设计提供相应的参考。 　　 3.1计算模型建立 　　 本模型边界尺寸取值：隧道长度取80米，开挖隧道洞径12米，开挖隧道洞高6米，模型左边界距洞边距离为5倍洞高，模型右边界取至5倍洞高，模型上边界取至地表，下边界距洞底为5倍洞高。隧道横坡坡度为47°左右，纵坡坡度约为26°左右。 　　 模型位移边界条件：左边界、右边界、后边界、地面边界法向位移约束，前边界隧道开挖面以下采用法向位移约束。模型共划分30586个节点，81231个单元。三维计算模型如下图2所示： 　　图2 火石岗隧道三维数值计算模型图 　　 3.2参数选取 　　 计算采用三维弹塑性分析，围岩材料采用摩尔-库仑模型，初衬及二次衬砌采用弹性模型，锚杆用提高加固区的围岩参数进行模拟。围岩参数根据火石岗隧道设计说明并参照公路隧道设计规范综合考虑确定。对于衬砌的参数，由于钢拱架与喷射混凝土实际上是紧裹在一起，共同变形、共同受力，所以钢拱架根据钢筋混凝土计算原理采用等效截面计算，即将钢拱架弹性模量折算给混凝土面层，计算公式如下式1。 　　 （1） 　　 式中：E―折算后混凝土弹性模量;E0―原混凝土弹性模量;Eg―钢拱架弹性模量;S