泥岩弹塑性损伤本构模型及其参数反演.docxVIP

泥岩弹塑性损伤本构模型及其参数反演 1 泥岩弹塑性损伤本构模型的研究 软岩的设计性质明显不同于硬岩和软土壤。特别是在水环境的作用下，软岩的物理力学特性会发生很大变化。对复杂多变的软岩工程,实施合理有效的工程控制,必须突破现有的岩体力学和土力学的学科框架,研究软岩的工程力学特征,构建能够考虑水化作用的本构关系。 Robinet针对饱和黏土岩的非线性膨胀问题建立了热弹塑性模型,并采用两类塑性准则描述黏土岩的力学行为;Chiarelli对泥岩进行了三轴压缩试验,指出塑性变形和各向异性损伤是相互耦合的,对矿物结构和含水率对泥岩力学性质的影响也进行了研究,并在此基础上建立了用于描述泥岩的弹塑性损伤模型,采用数值分析方法对试验结果进行了验证。Philippe Kolmayer建立了描述泥岩峰后效应、脆性断裂、软化以及膨胀效应的损伤本构模型,并给出了模型的实现方法及数值算例。Nathalie Conil在研究泥岩物理性质和力学性质的基础上,应用应变等效原理研究了塑性与损伤的耦合效应,引入了热动力势能描述泥岩的非线性行为。Conil根据泥岩的试验结果,针对材料力学性质的降低建立了各向异性孔隙塑性模型,采用Drucker-Prager屈服准则,损伤以二阶张量的形式表述,用于描述损伤对水-力耦合的影响,并给出了数值算例进行验证。Chun-liang Zhang对用于储存核废料围岩的硬黏土进行了热-水-力三场耦合试验,指出THM三场耦合试验非常复杂,伴随着遇水膨胀、热膨胀、热效应影响岩石的变形和强度以及裂隙的自愈合现象,认为孔隙水的吸附效应在THM耦合中起主导作用。Jia基于试验和细观力学机理,以孔隙塑性和连续介质损伤力学理论建立了在饱和和非饱和条件下泥岩的弹塑性损伤模型,并通过数值模拟对试验结果进行了验证。Abou-Chakra Guery基于试验结果对泥岩的弹塑性损伤本构问题进行探讨,并与室内三轴试验进行了对比,证实了模型的有效性。 泥岩由于其低渗透性、良好的蠕变性和遇水损伤自修复的特性,被认为是储存核废料的备选场地之一。本文以Boom Clay泥岩为研究对象,根据非排水条件下泥岩三轴试验结果,在弹塑性力学和不可逆热力学的基础上,建立了泥岩弹塑性损伤本构模型,导出了损伤演化方程,编制了弹塑性损伤本构模型及其参数反演程序,采用优化反分析法获得了泥岩本构模型参数。 2 泥岩应力-应变本构模型 Boom clay泥岩主要是由黏土矿物组成(30%~70%,平均值为55%),非黏土部分主要由石英(25%)和长石组成,其中黏土矿物由伊利石(50%)、蒙脱石(30%)、高岭石(10%)以及伊利石-蒙脱石胶合物组成,基本力学参数如表1所示。 对Boom clay泥岩进行了非排水各向同性固结试验(简称为CIU),在0.89~5.42 MPa压力下进行固结,固结程序为:以0.5 MPa的载荷增量施加到预定的围压,然后打开阀门进行反压饱和,一旦固结和孔隙压力耗散完成,对试样增加轴向应力进行剪切试验(如图1)。 根据试验结果将泥岩应力-应变过程分为4个阶段(如图2),即OA、AB、BC、CD 4段: (1)OA段,该段的应力-应变曲线基本成直线,为弹性变形阶段,A的应力称为初始屈服应力σc0。 (2)AB段,是岩石微裂隙开始产生、扩展的阶段,称为非线性弹性阶段,主应力差与轴向应变之间基本满足双曲线关系,B点的应力称为峰值强度σcu,也就是通常意义上的强度。 (3)BC段,称为应变软化段,在峰值强度之后,随着应变的增加,应力下降,岩石发生应变软化。轴向压力使试件形成破裂面,强度降低,应变增长。C点的应力称为残余强度σcr。 (4)CD段,称为塑性流动阶段,随着塑性变形的持续发展,最终强度不再降低,达到破碎、松动的残余强度,这个阶段可以认为是理想的塑性阶段。 3 mohsomb/点带能量损伤模型有关公式 损伤对材料弹塑性的影响体现在对弹性刚度的影响和对加载函数、塑性势函数及软-硬化参数的影响。本文构建弹塑性损伤模型: (1)在峰前区域,泥岩残余变形较小,曲线基本上呈非线性弹性形式,采用弹性损伤模型描述其力学行为; (2)软化现象开始时,泥岩峰值强度满足Mohr-Coulomb强度准则,在损伤的过程中的强度也满足Mohr-Coulomb强度准则; (3)泥岩的残余强度准则也满足Mohr-Coulomb强度准则。 根据热力学第二定律可知,损伤的过程实质上是能量耗散过程或不可逆热力学过程。损伤演化方程可以借助损伤耗散势函数得到,在应变空间内定义一个能量指标为 取损伤势函数为 假定损伤流动与损伤势函数的梯度方向相同,增量损伤变量可以写为 式中表示损伤乘子。 损伤势函数的相容方程为 于是可得到 最后得到损伤演化方程为 3.1 弹性损伤演化方程 由泥岩三轴试验曲线可以看出,在变形初