岩石应变软化及渗透率演化模型和试验验证

为预测承载岩石的应变软化和渗透率演化,基于Gebdykes白云岩的三轴试验结果,分析了围压对岩石弹性模量、破坏应变、峰值强度、强度退化过程、残余强度和剪胀扩容的影响规律。将岩石变形全过程简化为3阶段,使用强度退化指数、脆性模量系数和扩容指数改进FLAC中的SS模型,建立了考虑围压影响的岩石应变软化模型。基于淮南潘一矿煤、凝灰岩、巴里坤砂岩、山西安家岭泥岩的渗透率与体积应变实验数据,建立了基于体积应变增透率的岩石渗透率演化模型,与改进SS模型结合,建立了考虑围压影响的岩石应变软化和渗透率演化模型。利用本文模型分别模拟了安家岭泥岩和Gebdykes白云岩的三轴压缩、渗透率演化和体积扩容过程,结果表明:1体积应变渗透率演化模型能较好地描述体积应变与渗透率的关系;2本文模型能较好地模拟围压对岩石残余强度、峰后强度退化过程和剪胀扩容的影响,能较准确预测承载岩石的渗透率演化。     

三峡花岗岩峰后力学特性及应变软化模型研究 张 帆，盛 谦，朱泽奇，张勇慧

 通过对三峡花岗岩进行常规三轴压缩试验,得到不同围压下的应力–应变全过程曲线。基于弹塑性理论,通过试验数据拟合屈服面,分析研究花岗岩强度参数与峰后应变软化参量的关系,得出在花岗岩的应变软化过程中,黏聚力c随应变软化参量的增大而快速减小,而内摩擦角j在应变软化的过程中几乎保持不变的结论。在此基础上建立花岗岩峰后应变软化模型,利用FLAC3D程序对花岗岩三轴压缩试验进行了数值模拟,其结果与试验数据比较吻合。     

围压影响的盐岩压缩力学特性及非线性本构模型

为了更准确地描述盐岩的力学特性,对盐岩开展多组单轴、三轴压缩试验。根据试验结果,多组试验数据得到的莫尔包络线具有非线性特征,并且盐岩力学特征对围压有很强的相关性。单轴及低于5 MPa围压时盐岩应力应变关系曲线表现弹性、应变硬化、软化3个过程,而当高于5MPa围压时盐岩延性特征显著,应力应变曲线硬化阶段比例加大,扩容点前移,而软化过程有一定程度的减弱。基于Hoek-Brown准则的软化模型将应力应变曲线峰前阶段改进为弹性非线性弹性模型,得到一种新的非线性本构模型。结果表明这种非线性本构模型可以很好地吻合不同围压作用的盐岩应力应变曲线。     

软岩三维弹黏塑性本构模型

为了有效地描述软岩的软化特征,首先推导了以压为正统一强度理论表达式,该理论可以考虑材料拉压强度不等及中间主应力效应,适合岩土工程界以压为正的使用习惯以及岩土材料拉压强度不等的特性;其次推导了以压为正统一强度理论π平面的极限线;第三,基于推导的π平面的极限线,结合Yin-Graham模型,推导了软岩的三维统一弹黏塑性本构模型表达式。根据针对取自日本石川县能登半岛端部的硅藻质泥岩所进行的不同围压、不同加载速率条件下的应变控制式和应力控制式固结不排水三轴试验资料确定了模型参数,进行了弹黏塑性数值模拟分析。并与三轴试验结果进行了对比分析,结果表明:本文建立的软岩三维统一弹黏塑性软化本构模型的数值模拟结果与试验结果非常接近,可以很好地模拟软岩的应变软化特征。     

基于3种破坏类型的岩石损伤软化统计模型

为研究岩石峰后力学特性和本构关系,进行岩石的单、三轴压缩全程试验。基于不同阶段和不同破坏类型分别建立"分段本构模型"和"基于破坏类型的损伤软化统计模型"。分段本构模型通过屈服点、峰值点和残余强度点将岩石的应力应变划分为3个阶段,建立分段本构关系函数;基于破坏类型的损伤软化统计模型将岩石的压缩破坏划分为低围压(含单轴)下的环向张拉破坏、高围压下的剪切破坏以及中低围压下的张剪组合破坏,分别选取环向应变、剪切应变和体积应变作为损伤控制变量建立3种破坏类型的本构模型。结果表明,基于破坏类型的损伤软化统计模型能更好地描述岩石峰前及峰后的本构关系,揭示岩石内部损伤的发展过程。     

砂岩力学特性及其改进Duncan-Chang模型

 为了研究砂岩的力学特性,对砂岩试件开展了不同围压下的常规三轴压缩试验。试验结果显示,随围压增加,砂岩峰值应力、峰值点应变及残余强度均逐渐增大;当围压低于15 MPa时,砂岩弹性模量随围压增加也逐渐增大,但增大幅度逐渐降低;当围压在15 MPa以上时,其弹性模量则与围压无关。为了描述砂岩破坏过程的应力–应变响应,提出一种改进的Duncan-Chang模型,并根据岩石应力–应变曲线峰值点处斜率为0的特点给出模型参数的确定方法。利用砂岩三轴压缩试验结果对模型合理性进行验证。预测曲线和试验结果对比显示,该模型能够准确描述砂岩应变软化特性和不同围压下砂岩破坏过程中除初始压密阶段以外的其余4个阶段,特别是能够反映砂岩破坏后的残余强度。对模型特性的进一步分析表明,除应变软化特性外,该模型还可模拟岩石在高围压下的应变硬化行为,具有较强的适应性。     

考虑围压效应和强度脆性跌落的岩石全应力–应变曲线统计损伤模拟方法

为了建立能够较准确地描述岩石全应力–应变曲线的统计损伤模拟方法,针对目前统计损伤本构模型在模拟围压效应和强度脆性跌落方面存在的严重缺陷,首先,基于三轴压缩条件下脆性岩石基本变形特征,将岩石视为由软物质和硬物质两种细观化材料串联组成,建立宏细观物质变形分析模型;然后,基于瞬时应变分析方法和统计损伤理论,分别提出软物质和硬物质的变形分析方法,从而建立能够考虑围压效应和强度脆性跌落的岩石全应力–应变曲线统计损伤本构模型,给出参数的确定方法;最后,通过模型验证与分析,模型不但能够描述岩石初始宏观变形特征与围压之间的关系,也能够消除岩石轴向峰后理论变形曲线存在的随机性,从而能够较准确地模拟出岩石全应力–应变曲线,表明方法具有一定的合理性和可行性。     

双剪统一弹塑性应变软化本构模型研究

 根据统一强度理论应力不变量的形式,确定所要建立模型的屈服函数以及塑性势函数,考虑土体强度随加载过程的逐渐发挥,确定土体硬化函数。采用非相关联流动法则建立复杂应力状态下的土体弹塑性本构方程,并且对模型中所产生的奇异性进行分析和处理,指出模型中的奇异性是屈服函数、塑性势函数的流动矢量、以及硬化模量所引起的。给出模型中所含参数的确定方法,采用膨胀性泥岩的常规三轴试验模型进行验证,验证结果表明,所建立弹塑性本构模型能够模拟泥岩的应力–应变关系,较好地反映土体的应变软化特性。     

软弱岩石峰后应变软化力学特性研究

 软弱岩石给采矿工程中巷道支护和维护带来一系列棘手的问题,深入研究软弱岩石受力变形、破坏的机制和规律,对于保证巷道围岩的安全和稳定具有十分重要的意义。通过对软弱泥岩进行常规三轴压缩试验,得到不同围压下的全应力–应变关系曲线,然后依据峰后岩石任意一点应力状态均满足Mohr-Coulomb极限破坏条件的假设,建立以广义黏聚力 和广义内摩擦角 两个状态参数来表征的软弱岩石后继屈服面模型。在此基础上,利用试验数据绘制岩石峰后不同软化状态时的几组莫尔应力圆,通过“切线法”得出莫尔强度包络线的拟合方程,进而确定出不同围压条件下的 和 值,并借助Matlab软件对广义黏聚力 、广义内摩擦角 与等效塑性剪切应变、围压之间的关系进行最小二乘曲面拟合,得出软弱岩石峰后力学参数的软化规律,结果表明：随着围压的增加,广义黏聚力 值呈快速增加的趋势,而广义内摩擦角 值则显著减小;广义黏聚力 受岩石软化程度的影响也十分明显,从岩石峰值状态到残余状态 值迅速降低,平均降低53.88%,而广义内摩擦角 值在该软化过程中则基本保持稳定。最后,将得到的广义黏聚力 和广义内摩擦角 的拟合方程嵌入到FLAC内置应变软化本构关系中,并利用FLAC3D软件对模型的正确性进行数值模拟验证,结果表明数值模拟曲线与试验曲线比较吻合。     

花岗岩峰后力学特性试验与模型研究

深部围岩在开挖卸载过程中表现出的峰后复杂力学特性一直是工程界十分关注的问题,深入研究岩石峰后力学行为对深部资源开采工程具有重要意义。以深部立井马头门工程为依托,通过室内试验方法研究花岗岩峰后力学特性,采用非线性拟合方法获得花岗岩峰后软化模量与围压的指数关系式,假定岩石的剪胀角为恒定值,基于塑性理论构建考虑围压及剪胀角影响的岩石峰后应变软化模型;以FLAC3D为平台开发数学模型并进行验证,通过构建马头门巷道数值模型,分析深部围岩在应变软化条件下的破坏特征规律。通过研究可知,花岗岩峰后破坏具有脆–延性转化趋势,在高围压条件下,岩石峰后表现出塑性软化破坏特征,岩石峰后软化模量随着围压的增大而减小;通过FLAC3D进行数值验证可知,构建的应变软化模型与试验数据基本吻合,所建立的应变软化模型具有较高的可靠性;通过数值模拟方法分析深部马头门巷道围岩破坏特征可知,巷道拱顶及拱脚等局部区域出现了塑性剪切应变,与现场巷道围岩破损位置及深度基本相同。