隧道开挖过程的水-力耦合效应分析.docxVIP

隧道开挖过程的水-力耦合效应分析 在地下工程的建设和运营过程中，地下水通常受到影响。隧道工程通常都存在着不同程度的涌水或渗漏水情况, 水下隧道所处的环境使工程师不得不更加重视水的破坏作用。大量的工程实践证明, 隧道涌水的存在, 特别是施工期间, 不仅堵塞坑道、淹埋设备, 给隧道施工带来巨大的困难, 严重的还会造成人员伤亡。 隧道工程中由水引起的地质灾害是由水与岩土体相互作用产生的。隧道开挖后形成的二次应力场改变了原始地应力场分布及地下水渗流场分布, 从而影响岩体的结构, 引起岩体中地下水性态的改变。水的渗透力主要以静水压力和动水压力两种形式表现, 静水压力是一种表面力, 动水压力是一种体积力, 两种力叠加作用的结果可能使岩体引起劈裂扩展、剪切变形和位移, 增加岩体的孔隙度和连通性, 从而改变岩体的渗透性能, 最终影响岩体的稳定性。同时, 水的作用还要降低岩土体的强度。地下水是地下工程施工运营过程中诱发灾害的重要原因, 水下隧道所处的水环境更是不容忽视的, 应考虑渗流对隧道稳定性的影响。 应力对渗流的影响是应力通过改变岩体的体应变、孔隙率而影响岩体的渗透系数, 从而影响岩体的渗流场;渗流对岩体应力状态的影响是通过施加裂隙壁上的法向渗透压力 (面力) 和切向拖曳力 (面力) 而影响岩体的应力分布。 地下水渗流分析方法的发展基本经历了工程-试验研究阶段、解析数学方法研究阶段和数值模拟研究阶段。目前, 对岩体的力学性质、渗流规律及流-固耦合规律等开展了大量的研究工作[3,4,5,6,7,8,9,10,11,12], 但主要集中在渗流-应力耦合模型和地下工程涌突水预测方面, 而隧道工程施工过程的渗流-应力耦合模拟以及对模拟隧道开挖水力耦合计算渗流边界条件的合理确定方法研究较少, 一些问题还没有得到很好的解决。从现场测试及数值模拟结果来看, 隧道开挖引起的围岩渗流扰动带大于隧道开挖引起的应力扰动带。因而在建立隧道开挖问题的水力耦合模型以及在确定人为边界时, 应根据渗流场边界来确定计算模型的边界。本文把隧道开挖过程与渗流场边界的变化紧密联系起来, 根据渗流场边界建立合理的三维有限元分析模型, 研究水下隧道施工过程的渗流-应力耦合作用。 1 下隧道的渗流-应力耦合 应用MIDAS/GTS软件中的渗流-应力间接耦合来模拟水下隧道的渗流-应力耦合。先进行渗流分析, 由渗流分析中计算得到的孔隙水压力计算渗透力的效果, 再将得到的渗透力作用到模型上进行施工应力分析, 如此交替迭代得到最终结果。 1.1 孔隙水压力w 分析中使用的流动法则为达西定律, 建立的流动基本方程为: ??x(kx?Η?x)+??y(ky?Η?y)+??z(kz?Η?z)+Q=?Θ?t(1)??x(kx?H?x)+??y(ky?H?y)+??z(kz?H?z)+Q=?Θ?t(1) 式中:H——总水头; kx、ky、kz——x、y、z方向的渗透系数; Q——流量; Θ——体积含水率; t——时间。 该方程假定在任意位置、任意时刻微小体积的流入和流出的变化量与体积含水率的变化量相同。上述基本方程表现的是非稳定流的渗透方程, 而稳定流状态中流入和流出量随时间没有变化, 故公式右边为零: ??x(kx?Η?x)+??y(ky?Η?y)+??z(kz?Η?z)+Q=0(2)??x(kx?H?x)+??y(ky?H?y)+??z(kz?H?z)+Q=0(2) 体积含水率的变化依赖于应力状态的变化和岩土特性的变化。体积含水率的变化与孔隙水压力变化的关系如下: ?Θ=mv?uw(3)?Θ=mv?uw(3) 式中:mv——阻流系数。 总水头可视为压力水头和位置水头之和, 即: Η=uwγw+z(4) 式中:H——总水头; uw——孔隙水压; γw——水的容重; z——标高。 整理公式 (4) , 得到下面公式: uw=γw(Η-z)(5) 将式 (5) 代入式 (3) 得: ??x(kx?Η?x)+??y(ky?Η?y)+??z(kz?Η?z)+Q=mwγw?(Η-z)?t(6) 当标高一定时, 对时间的导函数为零, 则: ??x(kx?Η?x)+??y(ky?Η?y)+??z(kz?Η?z)+Q=mwγw?Η?t(7) 1.2 稳定流分析的基本原理 分析中使用加重残差的伽辽金 (Galerkin) 法表现的有限元方程如下: ∫V([B]Τ[C][B])dV{Η}+∫V(λ?Ν?Τ?Ν?)dV{Η}t=q∫A(?Ν?Τ)dA(8) 式中:[B]——动水坡度矩阵; [C]——单元渗透系数矩阵; {H}——节点水头向量; 〈N〉——形函数向量; q——单元边的单位重量; λ——非稳定流的阻流项,λ=mwγw; {H}t——随时间变化的水头,{Η}t=?h?t。 以上方程可简化为: [Κ]{Η}+[Μ