隧底岩溶对围岩稳定性影响的数值模拟研讨.pdfVIP

隧底岩溶对围岩稳定性影响的数值模拟研究 吴梦军 (重庆交通科研设计院重庆400067) 摘要隧道底部存在溶洞时，将对整个隧道围岩的应力、应变重分布产生较大的影 响。本文利用有限元分析程序对隧底存在溶洞的公路隧道围岩稳定性进行 了一系列的数值模拟研究，探讨了隧底溶洞大小及位置对围岩稳定性的影 响，总结了公路隧道围岩稳定性与隧底岩溶关系的一些规律。 关键词公路隧道岩溶围岩稳定性数值模拟 1引言 岩溶是指碳酸盐岩等可溶性岩层在水的化学和物理作用下产生沟槽、裂隙和空洞，以及由 于空洞顶板塌落使地表产生陷穴、洼地等现象和作用的总称。当公路隧道底部存在有岩溶时， 因受岩溶的发育及其与隧道间的距离等因素的影响，隧道开挖引起的应力、应变重分布将与无 岩溶公路隧道有着极大的区别。对于岩溶地区公路隧道技术的研究，国内有部分学者对岩溶 区公路隧道的施工工艺和方法进行过探讨【1j【2儿3]【4|。而在岩溶对围岩稳定性影响方面，文 献bj对拱顶和拱腰上岩溶对围岩稳定性的影响进行了一定的研究，而对于隧道底部岩溶对围 岩稳定性的影响研究，国内进行比较少。本文正是针对这一现状，就公路隧道底部岩溶对隧道 围岩稳定性影响，利用大型通用有限元软件ANSYS进行了一系列的数值模拟研究，分析了不 同溶洞大小及其与隧道间不同距离情况下围岩的应力、应变、塑性区等的分布特征，得到了围 岩位移与溶洞大小、溶洞与隧道间距离等的关系，为今后岩溶地区公路隧道的设计与施工提供 了一定的依据与指导。 2计算模型与计算参数 2．1计算模型 隧道工程问题为一无限域问题，而有限单元法是在有限区域内对结构进行离散，为使这种 离散不产生较大的误差，必须取足够大的计算范围。但计算范围太大，计算工作量及计算时间 又较大。根据文献【6J可得，对于地下工程开挖后的应力应变，仅在洞室周围距洞室中心点3。 5倍隧道开挖宽度(或高度)范围内存在实际影响。在3倍宽度处的应力变化一般在10％以 下，在5倍宽度处一般在3％以下。本次模拟计算范围水平方向两边各取3倍洞跨宽；垂直方 向，上边界取隧道净高的2倍，下边界为隧道净宽的3倍。模型边界条件采用位移边界，其下 36 边界垂直、水平位移均固定，左右边界固定水平位移。岩体单元采用程序中的PLANE42四边 形等参单元。划分网格后的模型如图1。 图1有限元计算模型及网格图 2．2计算参数 物理力学参数按《公路隧道设计规范》中IⅡ类围岩的参数选取，具体如表1。 表1模型物体力学参数 容重r(kN／m3) 泊松比肚 弹性模量E。(GPa) 内聚力c(kPa) 内摩擦角+(。) 19 O．35 2．1 10．O 45 岩体的初始应力主要包括岩体的自重应力和地质构造应力，地质构造应力只有在构造运 动比较强烈的情况下才予考虑。本次模拟计算中，围岩岩体初始应力场仅考虑了自重应力，按 隧道埋深90m计算，则由dI：埘(式中盯l为垂直均布应力，7为围岩容重，日为隧道埋深)计 算出垂直均布压力仃，=1．71MPa。 本数值模拟研究仅为研究隧道底部岩溶对围岩整体稳定性的影响，未考虑施工过程中 围岩的动态受力情况，开挖方法采用全断面开挖，模拟方法采用“生”、“死”单元来进行。当 杀死单元时，其质量、弹性模量、应力、应变等均被设置为0，从而可以实现对隧道开挖的模 拟。 3模拟结果与分析 本次数值模拟共进行了10种工况计算，分别模拟了不同溶洞大小及溶洞与隧道间不同间 距条件下围岩的应力、应变情况。通过模拟计算，可得到围岩的位移场、应力场、应变场及塑性 区分布等的等值线图及矢量图等，图3为部分竖向应变等值线图(图中R表示溶洞半径，三表 示隧道底与溶洞顶间的距离)。另外，对隧道拱顶、拱腰及拱脚处位移进行分析后得到竖向位 移曲线如图4。通过对各工况下围岩应力、应变、塑性区等的分析比较，可得隧道底部岩溶对 隧道围岩稳定的一些影响规律。 37