RFPA渗流原理..docVIP

渗流－损伤耦合方程 当单元的应力状态或者应变状态将满足某个给定的损伤阀值时，单元开始损伤，损伤单元的弹性模量有下式表达： (3.16) 式中：D 损伤变量，E 和 E0 分别是损伤单元和无损单元的弹性模，.这些参数假定都是标量。下面以单轴压缩和拉伸本构关系为例，介绍单元的渗透率－损伤耦合方程。 当单元的剪应力达到Mohr-Coulomb 损伤阈值时： (3.17) 式中：ф —摩擦角； fc —单轴抗压强度，损伤变量D按下式表达： (3.18) 式中：fcr —残余强度，其余的参数见图3.6，由实验可知[21-25,28]，损伤引起的试件渗透突跳系数增大，ξ=5，ξ的大小由实验给出，单元渗透系数的描述按下式： (3.19) 当单元达到抗拉强度ft损伤阈值时： (3.20) 损伤变量D按下式表达： (3.21) 式中：ftr —残余强度， 其余的参数见图3.6， 单元渗透系数的描述按下式： (3.22) 由于实验结果是在压缩状态下得到的，由于拉伸状态下渗透率的实验结果很少，所以本文的模型中假设应力和渗透率的关系满足负指数方程，并延伸到拉伸坐标轴。另外，当单元处于多轴应力状态并且满足摩尔库仑准则时，我们可以用最大压缩主应变来代替单轴压应变， 用平均主应力代替，这样就可以将以上表述的一维压应力作用下的本构关系推广到三维，所以单元的强度和其残余强度也按照广义虎克定律按相同的比例提高。 (1) 弹性状态： 当单元处于弹性状态，式3.12中=1，描述单元的渗透性只和应力状态成负指数函数关系。是一个跟孔隙率和孔隙结构相关的参数，在弹性应力状态下，材料的孔隙形状、大小发生变化，不存在数量、排列和位置上的变化，所以可以认为，=0表明单元为隔水单元，孔隙率n=0，没有孔隙水压的作用；=1表明单元为孔隙单元，孔隙率n=1，有孔隙水压的作用；01表明单元的孔隙率0n1，由于岩石的和n的定量函数关系存在较大的争论，本模型认为参数通过实验确定。同理，β反映了一定孔隙结构下的应力敏感因子，也只能通过实验确定。 (2) 损伤状态： 当单元应力诱发单元损伤时，其弹性模量按具有残余强度的弹脆性损伤本构方程进行弱化，式3.12中参数=5，描述单元损伤后渗透性有一定的提到，提高程度由实验确定；由于单元损伤软化，β增大，反映应力对渗透性敏感程度增加，也只能通过实验确定；保持初始值不变。由于拉伸损伤的渗透性演化规律较难通过实验得到，所以在模型中，单元参数演化的拉伸损伤和剪切损伤同样处理。 (3) 断裂状态： 当拉应变达到极限拉应变()，时，单元的完全丧失承载能力和刚度，其弹性模量赋予小值，和文献[29]对裂纹单元的处理方法一致，画黑形成裂纹单元，实验研究表明[30]，岩石试件达到峰值后，形成贯通裂纹的试件比形成非贯通剪切带试件的渗透系数还要高上百倍。所以在模型中认为裂纹单元的渗透系数和孔压系数和损伤单元相比，显著增加，式3.12中=1000，变成裂纹单元后，=1。岩石峰值后的渗透特性，取决于岩石性质，ξ值应该由实验确定。 (4) 裂纹闭合状态： 当裂纹单元的压应变达到极限压应变（）时，裂纹处于闭合状态，单元弹性模量随着压应力的增加而增高，按式3.23计算。文献表明[31]，闭合裂纹的渗透率明显减小，式3.12中=0.01，渗透率减小，=0，没有孔隙水压的作用。 E = (3.23) 在本文的数值模型中，考虑可单元性质的非均匀分布，就可以得到材料应力分布的非均匀性，同时确定了单元的以上四种基本状态，通过弹性模量和强度等力学性质的弱化和渗透系数、水压力的提高描述单元的损伤和破坏后的渗流—应力—损伤耦合作用机制，进而可以研究断裂的局部化以及裂纹萌生、扩展乃至渐进破坏的全过程。本文的模型中的裂纹模拟与分布裂纹模型和钝裂纹带模型中的方法类似，假定裂纹分散在整个单元，不用在模型中引入专门的奇异型单元，单元损伤后不用调节单元节点的排列顺序以及数目，程序上很容易实