基于有限元法的导流洞内塌落机制研究.docxVIP

基于有限元法的导流洞内塌落机制研究

1有限元分析方法

在分析平面滑动围岩环境下的洞室稳定性时，传统的有限元朝法以连续介质力学为基础，难以反映滑动介质在开挖和扰动后的不均匀坍塌趋势。同时，根据均匀变形压力负荷原理分析了洞室围岩的应力场，但该方程不反映在开挖和扰动后的每一个破裂面上的随机滑动和崩溃。

目前,规范求解围岩作用在支护结构上的压力常采取经验简化方法,如普氏分析法、泰沙基分析方法等,这些经验方法简单实用,但求得的围岩压力值没有考虑隧道断面的形状影响和局部的压力变化,更无法考虑隧道围岩复杂的地质条件和地形地表条件,在隧道结构设计中有一定的局限性。孙钧等[1~10]将有限元法与规范法对比,说明了通常情况下有限元法分析结果的合理性。

规范中的求解支护结构内力采用荷载结构法或边值法,也是常用的分析方法,其缺憾在于没有(结构力学法)或简单地考虑(边值法)支护与围岩的相互作用,即仅考虑支护结构,而将围岩对支护的变形约束之反作用用事先假定分布模式的方法予以确定,再作用上很可能并不存在的山岩压力等荷载,按结构力学的拟塌落拱荷载方法分析出支护结构的各个截面的内力,忽略了周围岩土体与支护结构的相互作用;岩土体非线性、非均匀性及各向异性等对支护结构的特殊影响在结构荷载法或边值法中更是无法考虑的。

有限元法可以方便地考虑以上问题,而且在通常情况下(如非浅埋、非离散、破碎情况),有限元的分析结果正在被越来越多的隧洞设计专家所采用,然而笔者认为:常规有限元法在分析中将围岩看作连续体,对于相对埋深较大且围岩较好的地下洞室,其上覆岩体较完整,支护结构主要承受围岩的较均匀变形压力作用,此时支护结构易发生压碎性破坏,对此种情况下的围岩稳定性进行常规有限元分析较符合工程实际情形,但对于低地应力区浅埋的松散围岩体却有明显的局限性;最具说服力的是鲁地拉导流洞进口段,洞径为15~17m、埋深为15~35m,节理裂隙发育,有经验的设计人员与本文作者均估计此处围岩稳定性与喷层安全性均可能出问题。但应用常规有限元分析之后的结果却令人不敢相信:喷层最大压应力不超过4.5MPa,无拉应力,最大位移为5mm,稳定性无任何问题!甚至连续降低围岩变形模量与强度参数c,?,喷层的应力并无明显增大,然而现场施工与监测却发现喷层出现了不小的拉应力,变形也早已超过10mm,可见常规的有限元分析结果显示出喷层支护结构的内力与围岩变形均很小,远小于工程经验与实测所认可的程度。针对这一反常的工程现象,本文结合近期诸多工程实例及其相应的数值分析方法与理念,认为基于连续介质理论上(即仅考虑围岩的变形荷载)的传统有限元法对这种松散浅埋围岩不再适用,从而将结构荷载法与有限元法结合,针对松散浅埋围岩提出了一种新的荷载模拟方法,即拟塌落拱有限元法和拟塌落拱离散化有限元法,试图将塌落荷载的概念、半连续松散块体变形的弹性或弹塑性性状、围岩的不均匀性、不连续性及不同初始地应力作用下对支护结构的变形影响均纳入统一的分析模型中予以考虑,为工程实际提供更可靠且更合理的参考。

2基于有限元法的分析方法

浅埋洞室围岩常常属于半连续、半离散的裂隙化岩体介质,完全离散的块体分析方法或离散元法较适用于这种半裂隙化的洞室围岩稳定性分析,但也存在局限性。笔者曾提出过一种分析此种半连续、可变形的裂隙化围岩的分析方法——离散化有限元法,但有限元法通常是以洞室开挖过程中相应部分的初始地应力释放和转移引起的变形压力为主要荷载的,一般不单独考虑塌落岩体荷载,这一点也被广大设计与分析人员所认同。

但针对如前所述这一反常的工程实际问题,本文提出了以下2种数值模拟方法,即拟塌落拱有限元模型和拟塌落拱离散化有限元模型,为了保证有限元计算的围岩压力荷载与规范中结构荷载法提出的围岩压力的一致性,根据规范规定的深埋和浅埋隧洞围岩压力确定方法首先确定隧洞的深浅埋,然后再确定塌落拱变形的破坏范围,如图1所示,最后根据规范中所确定的塌落范围建立考虑边界条件的离散化有限元模型。

2.1结构-支护相互作用模拟

在塌落拱的侧滑面边界上建立可张开或可滑动的摩擦界面单元,模拟可能的塌落荷载沿着侧滑面的下滑,从而对支护结构产生压力的机制,模型简图如图2(a)所示。

该方法主要是针对浅埋但围岩较为完整的洞室,模拟方法是:先根据泰沙基塌落拱理论确定塌落围岩的范围,将这一塌落范围内岩体的自重荷载作用在预先布置的支护结构单元上,应用有限元法计算支护结构在此荷载下发生的变形与内力,其与传统的结构荷载法不同的是它可以较全面准确地模拟围岩与支护结构的相互作用;与传统的有限元法不同的是它以塌落拱荷载代替了“自重场+二次开挖应力场”的有限元模式,特别适合不具备连续性的浅埋节理化围岩的情况。

2.2半连续缺失变形性状的弹塑性分析

对于松散的围岩体,不仅要考虑到塌落拱边界上的