岩石损伤软化统计本构模型及参数确定方法的新探讨

基于现有岩石损伤软化统计本构模型研究,通过探讨岩石损伤软化统计本构模型参数与岩石应力-应变全程曲线特征参数即峰值应力与应变的关系,建立起特定围压下模型参数与围压的解析表达式。引进岩石Mohr-Coulomb强度准则,建立不同围压下岩石应力-应变全程曲线峰值应力与围压之间的关系,再通过探讨不同岩石应力-应变全程曲线峰值应变与围压的关系,导出了具有普遍意义的不同围压下岩石峰值应变计算公式,从而建立岩石损伤软化统计本构模型参数确定的新方法,由此得到能够模拟不同围压下岩石应变软化全过程的统一损伤软化统计本构模型。该模型较同类模型具有参数少和易于确定等特点,理论计算和实测结果比较分析表明了该方法与模型的合理性。     

高应力区岩石统计损伤本构模型研究

为了反映高应力区岩石变形破坏全过程,通过引入双参数抛物线型Mohr强度准则,建立适用于高应力区岩体的微元强度度量方式,并假定由该度量方式得到的微元强度服从幂函数概率密度分布,结合连续损伤理论,考虑损伤变量修正,建立一个新的适用于高应力区岩体的统计损伤本构模型。以高应力区英安岩为验证对象,进行常规三轴压缩试验,依据岩石应力-应变曲线特征,给出模型参数确定方法。引用基于线性Mohr强度准则而建适用于浅部低应力水平岩体的相关统计损伤模型,利用该模型和所建模型对高应力区英安岩和相关文献中的高应力区砂岩、花岗岩试验数据进行验证,对比试验曲线和理论曲线,证明所建模型的可行性和适用性,并对英安岩累积损伤的扩展过程进行分析。     

基于改进Harris函数的岩石统计损伤软化模型

为了反映岩石变形全过程,通过改进Harris函数,假设岩石微元强度服从改进的Harris概率密度分布,结合连续损伤理论,考虑损伤变量修正,建立岩石统计损伤软化模型。采用YSJ-01-00岩石三轴流变试验机进行红层泥岩常规三轴压缩试验,基于红层泥岩试验数据对模型参数进行敏感性分析,阐述模型参数的物理意义,分析累积损伤的扩展过程。引用相关文献中二长花岗岩和砂岩的试验数据,利用本文模型对3种岩石材料的试验数据进行辨识验证,给出初始损伤系数的约束范围,对比3种岩石材料的试验曲线和理论曲线,证明本文模型的合理性和适用性。     

岩石的塑性、损伤及其本构表述

在工程条件下(常温、中低围压),岩石塑性变形和损伤变量的微观机制都是岩石构元中微裂纹的萌生和发育,因而塑性变形与损伤在岩石材料中应是同时出现的,它们的演化规律应是相互耦合的。在此基础上本文建立了岩石材料弹塑性损伤本构模型的理论框架。在这个框架中,由于损伤的出现,塑性变形率与屈服面是非正交的,而最后得到的本构矩阵具有对称性。非正交性与通常的实验资料一致。本构矩阵的对称性对数值分析和理论研究都是至关重要的。     

岩石统计渗流模型和统计损伤本构模型研究

采用Eshelby等效夹杂方法建立考虑渗流的岩石损伤本构模型是一种有效方法,但相关文献目前还极少。利用细观力学的Eshelby等效夹杂方法, 探索了考虑渗流和损伤的Helmholtz自由比能函数的确定,用连续介质损伤力学方法建立了相应的考虑渗流的岩石损伤统计本构模型,提出了考虑渗流和损伤过程的岩石破坏准则,建议了考虑损伤和应变的岩石渗透系数演化方程。与试验结果比较表明,所建立模型是合理的。     

基于新型损伤定义的岩石损伤统计本构模型探讨

针对几何损伤理论中损伤定义的局限性与不足,对岩石损伤定义进行了改进;在此基础上,将在应力作用下的岩石材料抽象为破坏与未破坏两部分,并根据这两部分不同的受力情况,建立出新型损伤模型;其次,利用统计损伤理论,建立出岩石统计损伤演化方程,进而建立起基于特定围压下的岩石损伤统计本构模型;然后,通过探讨模型参数与围压的经验关系,对模型进行合理修正,从而建立了反映不同围压条件下的岩石损伤统计本构模型,并通过探讨模型参数的物理意义,初步建立了反映岩石软化与硬化特性相互转化的临界物理力学参数的确定方法。该模型能够反映岩石软化与硬化特性相互转化的特征,与试验结果进行比较分析表明,该模型与实测结果吻合良好。     

一种岩石统计损伤本构模型的研究

基于岩石的非均质性,假设岩石微元粘聚力服从Weibull分布,微元破坏后其粘聚力降为0,结合统计理论和岩石剪切面上剪力平衡方程,建立了岩石发生剪切破坏时的损伤演化方程。该损伤演化方程考虑了破坏微元的抗压和抗剪能力,能够反映岩石峰后强度逐步丧失过程和残余强度。将该模型编入有限元程序,数值计算了不同围压下岩石应力-应变全过程曲线,与试验曲线比较取得了很好的一致性。     

基于残余强度修正的岩石损伤软化统计本构模型研究

岩石随着围压的增大,残余强度的增加幅度比峰值强度大,残余强度逐渐成为影响岩石全应力-应变曲线峰后段的主要因素,因此,建立岩石损伤软化统计本构模型时,对残余强度进行修正是非常必要的。基于岩石应变强度理论以及岩石微元强度服从Weibull随机分布的假设,考虑岩石峰后残余强度对损伤变量进行修正,在微元破坏符合Hoek-Brown屈服准则条件下,建立了能够反映岩石峰后软化特征的三维损伤统计本构模型;依据岩石试验全应力-应变曲线的几何特征,推导出本构模型参数的数学表达式,并基于花岗岩室内试验数据对模型参数进行探讨,研究了Weibull分布参数与本构模型的关系,探讨了损伤修正系数的取值和岩石累积损伤的扩展过程;将建立的本构模型理论曲线与4种岩石（斑状二长花岗岩、细晶大理岩、砂岩和粉砂质泥岩）不同围压下的常规三轴压缩试验曲线进行对比分析,结果表明,该模型能很好地描述岩石破裂过程的全应力-应变关系和表征岩石残余强度特征。这对于岩石损伤软化问题以及岩石加固处理措施的研究均具有重要的意义。     

考虑弹性模量变化的岩石统计损伤本构模型

为建立更为完善的岩石统计损伤本构模型,针对现有模型难以反映某些岩石在三轴试验条件下弹性模量随围压变化而改变的特性,从不同类型空隙对岩石变形性质的影响入手,深入分析了岩石弹性模量变化的机理,并在此基础上,依据弹性参数与岩石中未闭合裂隙数的关系建立了岩石弹性模量变化的分析方法。之后,引入统计损伤理论,考虑损伤阈值的影响和残余强度变形特征,建立了可模拟岩石变形全过程的统计损伤本构模型,并给出了基于三轴压缩试验曲线的模型参数确定方法。本文模型可反映岩石弹性模量随围压变化而改变的特性,且具备现有模型的优点,对岩石变形全过程模拟效果良好。最后,通过试验曲线与本文及同类模型理论曲线的对比分析,表明了本文模型的合理性与优越性。     

岩石损伤本构理论

本文在格森(Gerson)本构理论的基础上,引入了微观应变向宏观应变的转化因子,推出了一套损伤本构理论,建立了实际损伤模型,详细写出了叠代步骤,理论结果和试验结果吻合较好。     

基于统计强度理论的破裂岩体本构模型研究

借助损伤力学思想,基于统计强度理论,提出一种适用于深部工程破裂区破裂岩体的本构模型建立方法,并通过室内试验和数值试验进行了验证。将破裂岩体划分为无数微元立方体,微元立方体的强度与岩石破裂程度有关,且各立方体强度随机分布,故可用强度反映破裂岩体的破裂程度,据此提出一种破裂岩体本构模型建立方法。其中,根据岩石破裂面间的摩擦力做功等于材料破裂后释放的应变能,得到从力学角度定义的岩石破裂程度变量;另外,假定微元立方体强度分布服从Weibull分布,应力水平满足Hoek-Brown准则。利用泥质砂岩破裂岩体典型三轴试验结果,建立泥质砂岩破裂岩体本构模型,并进行了验证,结果表明模型计算曲线与试验曲线吻合度较好。利用离散元软件PFC进行了补充数值试验验证研究,证明了泥质砂岩破裂岩体理论模型的良好计算效果,进而证明提出的本构模型建立方法的可行性。     

考虑体积变化影响的岩石应变软硬化损伤本构模型及参数确定方法

针对现有岩石损伤模型的局限性与不足,将应力作用下的岩石抽象为空隙、损伤与未损伤材料三部分,以空隙率反映岩石体积或空隙的变化,以损伤变量或损伤因子反映岩石力学性态的改变程度,通过对各部分材料的微观受力及其与宏观应力和应变的关系进行分析,建立了可反映岩石变形过程中体积或空隙变化特性的新型岩石损伤模型;从岩石微元强度合理度量方法研究入手,引进统计损伤理论,建立了可考虑损伤阀值影响的岩石统计损伤演化模型,进而建立了能模拟应变软硬化全过程的岩石统计损伤本构模型,并提出了参数确定方法。该模型不仅能反映岩石应变软硬化特性,而且能反映岩石变形的初期压密和后期扩容的特性。通过实例分析与讨论,表明模型合理与可行。     

盐岩流变损伤特性及本构模型研究

借助盐岩单轴流变-声发射试验,分析了流变3个特征阶段与声发射事件数的对应关系,并基于Weibull分布构建了声发射事件数随加载时间的关系式,进而获得了损伤变量的演化关系式。对盐岩流变-声发射试验及流变-超声波试验数据拟合分析表明,该损伤变量表达式具有很好的适用性。进一步对损伤表达式中的参数进行了敏感性分析,揭示了形状参数、尺度参数对损伤值的影响规律。引入分数阶微积分理论及损伤理论,构建了基于损伤的Abel黏壶本构关系,通过替代西原模型中Newton黏壶的方法,求解获得了盐岩分数阶流变本构模型。由盐岩单轴流变试验曲线拟合分析可知,分数阶流变模型能很好地描述盐岩流变的3个特征阶段,尤其是加速流变阶段,且验证了西原模型为分数阶流变模型的一种特殊形式。     

基于统计损伤本构模型的改进桩-土接触面模型研究

建立准确合理的接触面本构模型是研究桩基承载力、沉降及变形等问题的前提。考虑桩-土界面初始剪切刚度的深度效应,提出基于统计损伤本构模型的无厚度接触面本构模型,根据剪切试验的数据对改进的接触面本构模型参数进行拟合,得到了不同法向应力下的模拟曲线。利用FLAC3D软件内嵌的FISH语言,对桩-土界面进行改进本构模型的二次开发,提出的无厚度接触面模型合理地模拟出桩-土界面刚度的非线性,成功实现了桩-土剪切试验过程的数值模拟。计算结果表明：提出的无厚度接触面本构模型能够合理地模拟桩-土界面的力学行为,可以描述接触面刚度的非线性变化,也适于程序化计算,拓展了接触面统计损伤本构模型的应用领域,为深入探讨桩基承载变形性状奠定了基础。     

含天然气水合物沉积物损伤统计本构模型

水合物含量、有效围压是影响含天然气水合物沉积物力学性质的主要因素,在忽略其他次要因素（包括水合物种类、试样颗粒大小、试验条件等）的情况下,水合物含量和有效围压是决定试样弹性模量的两个关键参数。在分析等效弹性模量与水合物含量相互关系的基础上,考虑有效围压的影响,建立了弹性模量与有效围压的幂函数关系;同时采用Drucker-Prager破坏准则来表示含天然气水合物沉积物微元强度,并假设其微元强度服从Weibull分布,从而建立了含天然气水合物沉积物的损伤统计本构模型,与不同有效围压下的试验结果及已有研究成果相比较,表明了所建模型能够很好地模拟三轴剪切条件下含水合物沉积物试样的应力-应变关系特性。此研究成果可对含天然气水合物沉积物工程性状的数值模拟提供参考。     

基于TCK损伤本构的岩石爆破效应数值模拟

脆性岩石在爆破载荷下的动态破碎行为是工程技术人员所关心的,采用数值法预估岩体破裂范围和损伤大小对施工安全具有重大意义,其中合理的岩石损伤模型是关键。假设岩体破坏近似服从Mises屈服准则,且考虑到岩体的应变硬化特性,把Taylor-Chen-Kuszmaul（TCK）模型中的拉裂损伤演化方程和材料双线性、弹塑性本构耦合到一起,并简洁地嵌入到有限元分析软件LS-DYNA中,通过半无限岩体中柱形药包爆破实例验证了其合理性与准确性。对球形药包爆破问题也进行了数值分析。该损伤模型及其数值模拟在工程中具有一定的参考价值。