研究高压引水隧洞衬砌透水设计.docVIP

研究高压引水隧洞衬砌透水设计 　　摘要：目前国内在高压引水隧洞衬砌的透水设计中没有较为详细和完善的计算方法，而且现阶段大多采用的是限裂设计和荷载结构法的方法对水工高压隧洞的衬砌进行不透水设计，但是不透水设计并不适合于高压引水隧洞，它将导致过高的配筋量。因此，基于体力理论，一种安全合理的有限元计算透水设计方法得到应用。此种方法不仅可以模拟运营检修工况和隧洞的施工过程，还可以分析应力场和渗流场的直接耦合作用，并计算出衬砌裂缝和配筋量。 　　关键词：高压引水; 隧洞衬砌; 透水设计 　　中图分类号：S611 文献标识码：A 文章编号： 　　1引言 　　现阶段大多采用的是限裂设计和荷载结构法的方法对水工高压隧洞的衬砌进行不透水设计，但是不透水设计并不适合于高压引水隧洞，它将导致过高的配筋量。现在在设计时出现体力理论，本文以此为基础，提出有限元计算透水设计方法，此种方法不仅可以模拟运营检修工况和隧洞的施工过程，还可以分析应力场和渗流场的直接耦合作用，并计算出衬砌裂缝和配筋量，可以较为科学有效的帮助进行高压引水隧洞的透水设计。 　　2耦合分析方法和设计原理 　　当分析裂隙岩石的渗流问题时，将岩石中的裂隙水流按流量等效原则抽象为等效连续介质考虑时，实际流速比达西流速要大若干个量级；对于恒定渗流问题，与时间无关。对单一介质或由多种介质组成的非均匀介质，只要各介质渗透系数相对值正确，实际流速与达西流速无论差别有多大，都不会对用水头表示的渗流场计算结果有任何影响，而只对渗流量有影响。因此，可以用等效连续介质模型分析恒定渗流问题。 　　采用等效连续介质模型进行渗流场和应力场的耦合分析时，分为两种情况：①介质的渗透系数为常量，与应力环境无关；②介质的渗透系数是应力环境的函数。对于工程问题，孔隙介质在荷载作用下发生变形，当荷载不是太大、孔隙的变形相对较小、不会对其渗透系数有明显影响时，可将孔隙介质作为常量处理。随着时间的增长，介质内的水逐渐排走，孔隙压力逐渐消散，有效应力随之变化，这表明了渗流场与应力场的相互作用。对于平面问题，基本方程如下： 　　以位移和孔隙压力表示的平衡微分方程为 　　 　　式中：G 和μ 分别为剪切模量和泊松比；u 为孔隙压力；γ 为重度； w 为位移；Vε 为体应变，Vε= 　　 　　根据变形协调条件，体积变形的速率为单位时间流入单位体积的净流量，并引入达西定律，得到连续性方程为 　　 　　式中：K 为 x和 z 方向的渗透系数；wγ 为水的重度；t 为时间。 　　联立式（1）、（2），组成应力场和渗流场耦合分析的基本方程，即著名的比奥（Biot）方程[5－6]，其有限元方程表示为 　　 　　式中： K 为刚度矩阵； K ′为与节点孔隙压力对应的节点力；K 为节点孔隙压力差产生的水力坡降引起边界排出的水量；δ 和β 分别为节点位移和超静孔隙压力。 　　式（3）右端 R 为各节点的外荷载，在固结过程中若荷载变化，则只要在相应时刻增加 R 即可；S 为前一时刻各节点的体积压缩量，由下式求得 　　 　　在求解时，从 t = 0开始，每次增加一个 Δt ，用式（3）解得δ 和β ，并求出 S ，用于求解下一时刻对应的值。如此反复至任意时刻，同时由δ 又可解出应力和应变的分布。这种方法将渗流场与应力场直接耦合求解，渗透系数为常量，无需进行渗流场和应力场反复迭代，只要按时间过程连续求解，即可得到全部结果，且每一时间步内不存在收敛问题。本文研究的重点是在稳定渗流情况下的隧洞衬砌。根据上述内容，将围岩等效为连续介质来模拟，不考虑介质渗透系数与应力环境的关系，采用渗流场和应力场的直接耦合方式进行隧洞在不同工况下的耦合计算分析。 　　3 隧洞衬砌的透水设计 　　3.1基本假设 　　哈德布特水电站工程高压引水隧洞全长13.55km, 沿途Ⅰ, Ⅱ类围岩约占 68%, Ⅲ, Ⅳ类围岩约占 32% , 与各类围岩相应的衬砌总长都不短, 分析中将其均归类为平面应变问题。 　　鉴于高压引水隧洞衬砌结构的主要荷载为水荷载, 计算中对Ⅰ, Ⅱ类围岩均假设与地层的初始接触应力为零, 即认为地层开挖后围岩变形均瞬时发生; 对Ⅲ, Ⅳ类围岩, 则取为在假设洞室开挖后立即施作衬砌( 即释放荷载取为 100% ) 的条件下, 按连续介质力学模型计算所得的地层与衬砌间接触应力的 20% ( Ⅲ类围岩) ～30%( Ⅳ类围岩) , 以包容由围岩地层变形的弹塑性及( 或) 时效特征引起的衬砌结构的内力。这类荷载习称形变压力, 研究表明实际发生的形变压力有可能比假设数据略高, 然而由于这类荷载的作用方向与渗流体积力相反, 计算中取用低值将可使结构偏于安全。 　　完建工况通常不起控制作用, 故分析中仅对高压引水隧洞的充水运营开始