围堰流-固耦合分析的渗透边界条件研究

围堰流-固耦合分析的渗透边界条件研究

1流-固耦合计算方法在节水工程中，洪水广泛存在于水库中。堤坝渗漏会导致管涌，引起滑坡甚至溃坝。长期以来，堤坝的渗流研究只局限于浸润面以下的饱和区，而忽略了浸润面以上水体运动状态的研究。加之渗流过程的复杂性，其一般表现为渗透特性和边界条件的复杂性，如渗透的各向异性和非均匀性，排水孔边界及自由边界等。有关该问题的计算方法大致可以分为两类：一类是解析解；一类是数值解。数值解法通常有有限单元法或有限差分法。目前国内外采用有限元法进行分析研究。而有限元法在处理这类问题时往往会遇到这样或那样的困难，有关问题及其处理方法在相关的文献中均有所研究。目前，用差分法分析流-固耦合的研究并不多见，其中FLAC3D(FastLagrangianAnalysisofContinuain3Dimensions）是由美国ItascaConsultingGroupInc开发的三维显式有限差分法程序，它可以模拟岩土或其它材料的三维力学行为。FLAC3D采用了显示有限差分格式来求解控制微分方程，并应用混合单元离散模型，可以准确地模拟材料的屈服、流体的流动、软化直至大变形。尤其在材料的弹塑性分析、大变形分析以及流-固耦合分析等领域有其独到的优点。鉴于此，本文采用FLAC3D对水中填筑围堰进行流-固耦合分析，并针对有限差分中流体的特征时间和渗透边界的敏感性进行分析，以求得可靠的流-固耦合分析。2流体流动间的相互耦合FLAC3D模拟多孔介质（如土体）中流体流动时，流体的模拟独立于力学计算。其主要通过孔隙水压力的消散引起土体中位移的变化，这一过程包含两种力学效果。孔隙水压力的变化引起结构体中有效应力的变化；孔隙水压力的变化又引起流体区域的变化。流体在孔隙介质中的流动依据Darcy定律，流-固耦合过程满足Biot方程。流体在多孔介质中流动时，主要引起以下几个变量的改变，即孔隙压力、饱和状态和渗透流量。这些变量通过流体质点平衡方程，Darcy定律间的相互关系来描述流体的流动。本构方程主要表现孔隙压力、饱和状态、体积应变和温度的变化关系，进而实现温度-流体-固体三者间的耦合。其三者间的相互耦合有限差分方程见参文。在利用FLAC进行流-固耦合分析时，渗流对结构作用的好坏，不仅取决于基本方程的建立，而且在实现过程中还要考虑特征时间和边界条件两个重要的因素。这两个因素直接影响到结果的可信度，因此，分别加以介绍。(1)特征时间特征时间即为流体在孔隙介质达到稳态流动的最小时间，定义如下：式中cL为特征长度，。其中V是流体所流过区域的体积；S是流体流过区域的面积；c是流体的扩散率。式中k为渗透系数；M为Biot模量；α为Biot系数；α1=K+43G;K为排水介质体积模量；G为孔隙材料的剪切模量。当考虑孔隙材料的不可压缩性时，流体扩散率式中Kf是流体的体积模量；n是孔隙率。只有知道了流体的特征时间，才能在计算中确定流-固耦合的计算时间。否则，渗流的计算将得不到满足。(2)边界条件在任何计算中，边界条件是很重要的，它的好坏直接影响到计算的结果。在FLAC3D中，有4种类型的边界条件，它们分别是：(1)给定孔隙水压力。(2)给定边界外法线方向流量分量。(3)不透水边界。不透水边界程序中默认。(4)透水边界。透水边界采用如下形式给出：式中qn是边界外法线方向流量分量；k是渗透系数m/s；γ为水的容重；p是边界面处的孔隙水压力；ep是渗流出口处的孔隙水压力。其他边界条件，如部分边界上的常压力和常流量，则需通过FISH语言来实现。3基岩分布特征泰安抽水蓄能电站围堰填料采用上游水库开挖的砂粘土、砂石土的混合料，个别块石粒径达到50cm。围堰填筑采用自然堆坡，迎水坡坡比为1:1.5，背水坡坡比为1:1.7。堰顶顶部宽度为6m。堰体填料约为9.45×104m3，高压喷射灌浆防渗心墙约为13362延米。土工膜施工约2054m2，混凝土底座约为102m3，图1和图2分别给出了标准断面图和围堰平面图。围堰处的河床覆盖层由砂砾石层组成，一般厚度2～5m，透水性强，但已被防渗墙完全隔断。其下的全风化或强风化岩层厚度一般在7～15m之间，其中防渗心墙深入其中4.2m左右，没有完全隔断下部的强风化层。根据地质勘察得到的吕荣值可以换算出基岩的渗透系数，并按照渗透系数值划分为弱透水层和强透水层。堰下基岩上部一定范围内属弱透水层，是整个围堰中渗透性最强的连通区域，因此，渗流将主要在该部分基岩内发生。基岩分布具有明显的三维特征。主要材料性质见表1。由于围堰属于临时性建筑物，同时，考虑到工期要求，设计单位在初设时围堰下岩基不设置帷幕。本文渗流分析，按岩基不设帷幕条件计算，取防渗墙嵌入基岩为1.0m。4维流-固耦合分析方法根据围堰及其防渗结构的平面图和围堰的轴线图可以看出围堰堰高是随水深变化的，基岩的空间分布具有很明显的三维特性，因此，本文对初步选定的开挖和不开挖二个方案进行三维流-固耦合分析。4.1渗流分析的有限差分网格有限差分法计算深度取至基岩弱透水层一定深度，由于背水坡侧抽水后开挖约8m左右，故堰基取为12m，此时，计算深度已至弱透水层一定深度，可视为不透水层；另外，堰基两侧也各取一定深度（基岩深度取12m），以考虑在基岩内出现绕岩体的渗流情形，更远处的绕流将十分微弱，从而，可视为不透水层。渗流分析采用六面体单元，图3给出了含防渗心墙围堰的三维渗流分析的有限差分网格。弹塑性有限差分分析选用Mohr-Coulomb屈服准则，渗流分析选用Darcy定律。以围堰轴线为z轴，垂直地表方向为y轴，坐标原点选在S4处。考虑围堰主要用于挡水，为了节约计算，故在S2～S4方向水深小于2.95m处的围堰不进行计算。因此，S2～S4方向取为从245.07m到水深等于2.95m；同理，由于S5～S6段水位比较浅，最深处为2.95m，不会对围堰的稳定性造成危害，所以，S5～S6段没有被列入计算范围，只考虑S5～S4段，取为197.17m。围堰高度随地基基础的变化而变化。覆盖层也随地表高度的变化而变化，取其平均厚度为4m。以下均为风化层，取为8m。共18270个区域，16048个节点。浸润线的确定是通过透水地基上的均质土堤，由于受地基透水的影响，堤身浸润线降低。根据式(1)可分别求得开挖和不开挖时的特征时间分别为：不开挖时，tc=0.1×106s；开挖时，6ct×=108.1s。4.2渗流对设施的影响图4和图5分别给出了围堰及堰基在不开挖和开挖时的孔隙水压力和塑性区图，图中孔隙水压力等于0处是浸润线的位置。图4给出了堰体在不开挖时的塑性破坏区和孔隙水压力的变化关系，从图中可知，堰体由于受渗流的影响，致使堰体在背水坡坡脚处发生小范围塑性破坏，由于心墙的隔水作用，围堰和堰基处也发生了部分单元破坏，但没有形成贯通破坏。说明在背水坡处不会形成堰体和堰基的整体滑移破坏，但在坡脚处可能形成局部滑移。在迎水坡侧面由于渗流发生，加之防渗墙的作用，堰体发生大面积的塑性破坏，同时，与堰基的塑性破坏连在一起，导致堰体和堰基边坡共同形成反向滑弧破坏。如果将破坏区用一条线连起来，就形成了堰体发生滑移时的滑移线。从孔隙水压力看，以防渗墙为界，迎水坡侧的压力水头明显高于背水坡侧的水头压力，表明防渗墙对结构稳定的影响。围堰建成后，要进行抽水开挖，为了便于分析，也取最大水深的剖面如图5所示。首先堰体塑性区与不开挖相比没有太大的变化，在堰基处也没有形成连通的破坏区，说明堰体基础在开挖后，总体没有发生太大的变化，只是在开挖区附近产生了一定的影响。从塑性区图上看，开挖的影响只是局部的，没有超过预留的3m平台区。从孔隙水压力图上可见，水头线产生了明显的下降，尤其是在开挖区一侧，孔隙水压力变化比较明显，说明渗透流量明显增大，在背水坡侧，孔隙压力也有所下降，但不很明显。主要是有防渗墙的缘故，阻止了渗流的发生。从计算的渗流量来看，在不开挖条件下，围堰最大剖面单宽流量为6.27×10-4m2/s，在开挖条件下，围堰最大剖面单宽流量为1.52×10-3m2/s。二者相差比较大，说明开挖对堰体及基础的渗流产生很大的影响。另外，二者均是在不设帷幕条件下计算得到的流量，但流量都不是很大，表明围堰可不设置帷幕。5流-固耦合分析方法的应用(1)基于非均质各向异性的饱和-非饱和土的三维稳定渗流问题的基本控制方程，并给予正确的特征时间和边界条件，对泰安抽水蓄能电站进出水口围堰进行了三维流-固耦合分析，其计算结果具有良好的规律性，表明这一方法特别适用于分析形状复杂、结构分区多，且各分区渗透性差异较大的大型土石坝的流-固耦合分析。(2)对泰安抽水蓄能电站进