温度场耦合.DOCVIP

大坝及其周围地质体中渗流-温度-应力耦合的研究综述 周慧敏 梅明荣 (工程力学系 河海大学 南京 210098) 摘要:分析了大坝周围地质岩体中渗流场、温度场与应力场的相互作用，对三场耦合重要参数及控制方程进行了系统的分类与总结,并对工程使用进行归纳展望。 关键词：大坝地质岩体、渗流场、温度场、应力场、耦合分析 1 引 言 据统计，30 % ~40%的水电工程大坝失事是由地下水流渗透压力、温度应力引起岩体变形，而岩体变形又使其渗透性能发生变化，即渗流场、温度场和应力场（THM）相互作用是溃坝的主要原因。考虑单一力的作用解决问题效果往往并不理想，各种作用力相互交叉和渗透十分突出。因此，在设计时为了保证大坝长期安全运行有必要进行大坝-岩体系统THM 全耦合分析研究。近年来，THM研究方兴未艾，在核废料储存库防污染问题上已经扩展到THMC。经过大量阅读和深入思考，作者对工程中THM耦合分析的研究成果进行系统的评述与小结。 2 渗流场、温度场和应力场之间的相互作用 渗流场－应力场耦合：坝体-岩体系统中渗流场的变化引起渗流作用力包括渗流体积力和渗透压的改变,其中的应力状态改变使渗透系数改变进而改变渗流场。有两种可能：（1）贯通裂隙之中水压力的传递过程和扩容劈裂作用,诱使裂纹不断地扩展,对岩石强度和破裂模式的影响十分显著，裂隙的扩展导致孔隙水压力和梯度分布的非均匀性,重新改变应力分布形式,相互贯通的裂隙作为新的渗流通道，将使渗透系数增大；（2）裂隙岩体受法向、侧向应力作用时，岩体变得密实，导致裂隙岩体渗透系数降低。 温度场－渗透场耦合：一方面，坝体-岩体系统中渗流场改变，将使渗透水流参与坝体-岩体系统中的热量传递与交换，从而影响温度场的分布。另一方面，坝体-岩体系统中温度场的改变，形成的温度梯度引起水的运动改变坝体-岩体系统渗透系数。而岩石的弹性模量及强度都是温度的函数,因此温度的改变也将导致坝体-岩体系统应力场的变化。 温度场－应力场耦合：当温度场发生变化时, 将在坝体-岩体系统内产生热应力，与温度有关的固体特性也要发生变化。研究表明，岩石的导热系数随温度的升高而下降,比热、导温系数则随温度的升高而增加，裂隙岩体温度升高，裂隙封闭，开通裂隙数密度减少，导水性降低，温度场对应力场的影响还体现在温度的变化将导致孔隙流体压力较大的变化；当应力场发生变化时,存在岩体内部力学能量转换等。 渗透场－温度场－应力场耦合：裂隙岩体三场之间的耦合作用通过地下水的渗流运动、裂隙岩体的变形及热量的转移得以实现。当裂隙岩体内发生地下水的渗流、裂隙结构面的变形及热量的转移(包括热传导和热对流)时，在合适的地质环境下渗流场、温度场与应力场之间会产生耦合作用。 3 耦合场的支配方程 以渗流力学、传热学理论、连续介质力学为基础，耦合应用了如质量守恒、能量守恒等理论，假定岩体中水的流动符合Darcy定律： （1） 其中：V为渗流场中任一点的渗流速度，K为渗透系数，J为水力梯度。 假定岩体中热传导符合Fourier定律： （2） 其中T—温度，J—热流密度，K—为导热率。 能量守恒方程：根据热力学第一定律，单位体积中单位时间内，由外界传入系统的能量与内部热源产生的能量之和等于物质内能的增量与对外界做功之和。 3.1 渗流场-应力场耦合 双场耦合方程多采用机理分析法。此两场耦合方程已经较为成熟，控制方程格式基本如下[1]： （3） （1）渗流连续性方程，根据质量守恒定理，假设介质变形-应力为线弹性关系，将连续方程、状态方程、Darcy定律结合起来：水压梯度引起的流动和渗流源之和恒等于流体在介质内的累积率。所谓的流体在介质内的累积率分别由密度变化与孔隙度变化引起。孔隙变化有三个直接原因：体积应变；液体孔隙压力引起的颗粒大小改变；有效应力引起的颗粒大小改变。其中后两个影响相对较小可以忽略。 （2）弹性力学应力与应变协调方程，渐变性材料应力引起的应变用位移表示，包括弹、塑性应变所引起的应力。其建立的基础需要形成合适的能够描述受力表面变化的硬化规律及能够描述裂隙材料本构关系的流动规律。 裂隙岩体渗透系数是最重要的水力参数，裂隙岩体渗流渗透系数的非均匀性非常突出，导致其各向异性非常明显，岩体裂隙隙宽、单一裂隙体表比等影响因素难以确定。叶源新[2]提出渗透系数k是孔隙变化量的函数，参照己有的试验成果，设定渗透系数与孔隙率成指数关系： （4） （5） 式中：a为耦合参数，表征应力、应变对渗透系数的影响程度，由试验确定；p为孔隙水压力；为孔隙变化量；为体积应变；Q