岩体水力学基础(五)——岩体渗流场与应力场耦合的裂隙网络模型

岩体是由裂隙和其间的岩块组成。当岩块处于近乎隔水状态时，岩体中的渗流可看作为裂隙网络流。本文通过分析岩体裂隙网络渗流场、应力场以及它们之间的相互力学关系，提出了岩体渗流场与应力场耦合的裂隙网络模型及数值计算方法。     

岩体裂隙系统渗流场与应力场耦合模型

岩体系统具有复杂的结构。一般认为，岩体系统是非均质各向异性不连续的多相介质体系。当岩体以裂隙为主，且其分布较密集时，可将岩体系统看作等效连续多相介质体系。本文运用等效连续介质理论，提出了两种岩体裂隙系统渗流场与应力场耦合模型：一是以渗透水压力与隙变形关系、应力与渗透系统数关系为基础，建立渗透系数张量计算公式，进而建立等效效连续介质渗流为数学模型。以裂隙岩体应变张量分析为基础，建立裂隙岩体效应力张量     

岩体水力学基础（四）：岩体渗流场与应力场耦合的等效连续介质模型

本文把裂隙岩体看作等效连续介质，用渗透系数张量描述岩体的渗透性能；用应力张量描述岩体的力学性能。介绍岩体渗流场与应力场耦合的等效连续介质模型。     

岩体水力学基础（五）

岩体是由裂隙和其间的岩块组成。当岩块御于近乎隔水状态时，岩体中的渗流可看作裂隙网络流。本文通过分析岩体裂隙网络渗流场、应力场以及它们之间的相互力学关系，提出了岩体渗流与应力场耦合的裂隙网络模型及数值计算方法。     

岩体渗流场与温度场耦合的连续介质模型

孔隙型岩体、密集裂隙型岩体以及孔隙—密集裂隙型岩体的渗流问题均可视为连续介质或等效连续介质渗流问题。本文从理论上分析了这种类型岩体渗流场与温度场相互影响、相互作用的机理，提出了岩体渗流场与温度场耦合分析的连续介质数学模型，并讨论了该数学模型的有限元数值求解方法。     

泄洪雾雨区裂隙岩质边坡饱和-非饱和渗流场与应力场耦合分析

岩质滑坡除与地质条件、降雨等因素有关外,泄洪雾雨也成为其不可忽视的催化剂。通过裂隙岩体渗透性与应力关系,裂隙岩体的弹塑性本构关系以及渗流场的控制方程,建立起饱和非饱和渗流场与应力场耦合的等效连续介质模型。编制了有泄洪雾雨影响的饱和-非饱和渗流与应力耦合三维有限元程序,通过迭代求解达到两场耦合目的。以实际工程为例,根据边坡的地表信息、开挖信息、岩层分界信息以及排水洞和帷幕信息等,建立起边坡地质模型。通过耦合计算,详细分析了耦合后边坡岩体的变形、应力以及塑性区等。计算结果表明,雾雨区的边坡稳定与渗流场变化有着重要关系,在评价边坡稳定与否时,要考虑渗流场对其影响。计算成果为实际工程的安全评价提供了一定的科学依据。     

裂隙网络非连续介质渗流场与温度场耦合分析研究

在岩体水力学研究中,一个重要问题就是岩体渗流场与温度场的耦合分析。笔者关注了国内外相关的研究,对其发展的趋势进行了展望。系统综述了裂隙岩体渗流-温度耦合的研究情况,阐明了裂隙岩体三种渗流模型的特点、适用范围及裂隙网络非连续介质模型的研究概况,分析了裂隙岩体渗流-温度耦合作用研究及工程应用。最后,提出基于裂隙网络非连续介质渗流温度耦合分析是目前亟待解决的课题,并指出今后需要深入研究的问题。     

考虑动力压力裂隙网络岩体渗流应力耦合分析

渗流对裂隙网络岩体的作用力包括渗透静水压力和渗透动力压力（裂隙壁切向拖曳力）两部分。较全面地引入渗流与应力之间的相互作用关系，同时考虑渗透静水压力和切向拖曳力的作用，以裂隙（等效）隙宽与岩体应变的关系为桥梁，建立隙裂网络岩体渗流场与应力场耦合分析的数学模型，并采用有限元数值方法进行某工程实例的双场耦合分析，验证所建立模型的可靠性。     

核废料贮存围岩介质THM耦合过程的力学分析

针对核废料贮库围岩介质中的THM耦合作用下的岩体力学响应的特点,分析了THM耦合过程,利用O'C onnell的裂隙岩体等效弹性模量关系式中的裂隙密度概念意义,建立了温度作用下的裂隙密度与O da提出的裂隙张量之间的关系式;以裂隙结构面的开度、岩体裂隙数(包括受温度影响开通裂隙数)、裂隙连通率、附加应力、剪切膨胀等为研究对象,给出了具有THM耦合的渗流系数张量;在理论分析的基础上,建立了应力、应变、温度和热力学压力的关系方程式,给出了核废料贮库裂隙岩体介质热-液-力耦合的定解方程。为探讨核废料贮存裂隙岩体水热耦合迁移以及应力响应特征,提供了力学模型。     

裂隙岩体渗流场的无单元法模拟及渗流特性研究

裂隙岩体因含有发育程度不同的裂隙、节理和断层等不连续面,致其渗透性具有各向异性、不连续性等特点,因此传统的有限元法对分布密集的裂隙岩体渗流场求解有一定的难度。本文提出了采用无单元Glaerkin法求解有自由面裂隙渗流问题,并推导了无单元法求解渗流场的基本方程和积分格式,给出了应用罚函数法处理渗流边界条件和自由面处理方法。采用IDL语言编制了二维无单元法计算软件LIDAREFM。文中以北京怀柔桥梓镇某裂隙岩体边坡渗流场计算为例,研究了复杂裂隙共同作用下渗流场特性和自由面分布,讨论了不同开度、不同连通程度的裂隙对渗流场的影响。研究结果表明:无单元法可以较好地解决有密集裂隙的岩体渗流场的求解问题,实现了裂隙处结点任意加密以及积分网格的独立布置,避免了对有自由面和裂隙穿越的子域的重新处理,简化了渗流问题的求解过程。     

裂隙网络非连续介质滲流场与温度场耦合分析研究

在岩体水力学研究中，一个重要问题就是岩体渗流场与温度场的耦合分析。笔者关注了国内外相关的研究，对其发展的趋势进行了展望。系统综述了裂隙岩体渗流一温度耦合的研究情况，阐明了裂隙岩体三种渗流模型的特点、适用范围及裂隙网络非连续介质模型的研究概况，分析了裂隙岩体渗流一温度耦合作用研究及工程应用。最后，提出基于裂隙网络非连续介质渗流温度耦合分析是目前亟待解决的课题，并指出今后需要深入研究的问题。     

高渗压条件下压剪岩石裂纹断裂损伤演化机制研究

探讨了高渗透压作用下压剪岩石裂纹的起裂规律及分支裂纹尖端应力强度因子的演变规律,建立了高渗压下裂隙岩体发生拉剪破坏的临界水压力值和初裂强度判据,分析了不同渗压作用下裂纹的扩展情况表明,渗透压的存在加剧了分支裂纹的扩展,高渗透压作用下分支裂纹扩展由稳定扩展变成不稳定扩展,并导致分支裂纹尖端岩桥剪切破坏,同时考虑分支裂纹的相互作用,建立了高渗透压作用下压剪岩石裂纹体岩桥剪切贯通的断裂破坏力学模型,最后依据裂隙岩体的损伤力学效应研究了岩体的初始损伤及损伤演化柔度张量,提出了高渗压下压剪岩石裂纹渐进破坏的损伤演化方程。该理论为定量研究高渗压下裂隙岩体的失稳破坏提供了理论依据。     

渗透压力对裂隙岩体损伤破坏的研究

基于Betti能量互易定理推导了含水裂隙岩体初始损伤柔度张量和损伤演化方程,建立了考虑渗透压力的裂隙岩体脆弹性断裂损伤本构模型。并提出采用渗透压力附加柔度张量的概念来定义渗透压力对岩体力学特性的影响。工程算例结果表明：当岩体中有地下水活动时,考虑渗透压力对损伤区的影响是十分必要的。     

渗透压力作用下加锚裂隙岩体围岩稳定性研究

岩体中裂隙水的存在加剧了岩体结构围岩损伤;锚杆作为岩体工程的重要支护手段,被广泛应用于岩体加固工程中。结合无损伤材料的柔度张量的概念,用裂隙附加柔度张量来表示裂隙对岩体的损伤影响;引入渗透压力附加柔度张量的概念来定义渗透压力对裂隙岩体强度的损伤影响;利用附加刚度来反映锚杆对裂隙岩体的加固作用,附加刚度求逆即得柔度张量。对渗透压力作用下锚固裂隙岩体渗流场与损伤场相互作用及耦合机制进行了深入研究;基于能量互易定理,自洽理论等,分别在压剪和拉剪应力状态下推导了渗透压力作用下加锚裂隙岩体等效损伤模型,建立了相应的计算方法,应用半解耦方法对理论模型进行了有限元程序化。在有限元程序中,将渗透体积力化为等效节点力,参与各单元间应力调整且表现为岩体变形过程中刚度的降低;锚杆锚固力等效为单元节点力,并体现为岩体整体变形刚度的提高。结合工程实际,着重讨论渗透压力作用对隧道围岩稳定性的影响。分析表明,渗透压力对岩体围岩变形破坏及稳定性的影响是不可忽略的。     

低温环境下含表面裂隙硬岩温度场及冻胀演化过程分析

寒区含表面裂隙硬岩与大气连通,易受外界低温环境引起裂隙内部水分凝结,造成裂隙带附近形成温度场突变,进而随水冰相变产生冻胀力,二者随着冻结活动处于动态演化之中,涉及到温度场、渗流场及应力场等多场问题。目前的研究大多依据数值模拟结果,描述裂隙岩体温度场、应力场分布特征及规律,考虑从基础理论解析计算角度,探究低温环境下裂隙岩体温度场、冻胀力演化规律鲜有报道。针对含表面裂隙硬岩在低温环境下发生的冻结现象,从传热学理论、相变理论、弹塑性力学及断裂力学理论等多角度出发,探讨随着冻结深入含表面裂隙硬岩内部温度场及冻胀荷载演化过程,提出对应解析表达式,引入多场耦合数值程序Comsol-Multiphysics进行数值模拟验证,结果表明,文中提出的温度场及冻胀荷载演化解析表达式与数值模拟结果具有较好吻合度,可满足对含表面裂隙半无限板模型温度场、应力场演化过程分析,对于裂隙岩体在低温环境下温度场、应力场演化及热-力耦合分析具有较好借鉴价值。     

岩体系统渗流场与应力场耦合的广义双重介质模型的应用研究

本文通过实验资料分析, 提出了裂隙岩体应力与渗流之间的分形几何关系式。以研究区岩体系统结构为基础, 提出了岩体系统应力场与渗流场耦合的广义双重介质模型, 并将该模型应用于某水利水电工程的渗流与控制中, 对坝区岩体进行了应力场和渗流场的耦合分析, 取得了比较满意的结果。