基于修正Mohr-Coulomb准则的弹塑性本构模型及其数值实施

针对Mohr-Coulomb准则高估岩土体抗拉性能的局限性,建立考虑最大拉应力准则的修正Mohr-Coulomb模型;系统地论述隐式本构积分算法的主要内容,推导相应的一致性刚度矩阵。以ABAQUS软件为平台,采用向后欧拉隐式应力积分算法编制了UMAT本构程序,对单轴拉伸试验和三轴压缩试验进行数值模拟,对比分析ABAQUS自带模型和自编模型的优劣,结果表明编写的修正Mohr-Coulomb模型能够有效地反映岩土介质的抗拉性能,弥补了ABAQUS自带模型的不足。     

岩石摩尔-库仑弹塑性损伤本构模型及其主应力隐式返回映射算法研究

在弹塑性损伤理论框架内考虑岩石的塑性变形机制和刚度退化,建立基于Mohr-Coulomb(M-C)屈服准则的弹塑性损伤模型,采用内变量即等效塑性应变表征岩石损伤变量的演化。由于M-C屈服准则在应力空间为一个六棱锥,在数值实施过程中六棱锥角点和棱线上的应力更新存在"奇异性"问题,角点光滑化方法可以处理该问题,但其不可避免的导致近似的计算结果。在M-C本构数值积分算法的基础上,推导弹塑性损伤本构方程的主应力空间隐式返回映射算法,包括弹性预测、塑性修正和损伤修正3个主要计算步骤。在塑性修正过程中,针对流动向量返回到主平面、左右棱线和尖点3种情况分别进行讨论,从主应力空间的角度出发解决"奇异性"问题。采用面向对象的编程方法,使用C++语言开发弹塑性损伤本构求解程序(RDM-C),并采用单轴压缩试验、地基和洞室算例对程序计算的结果进行分析和验证。研究结果表明,所建立的弹塑性损伤本构模型能够较好地描述岩石材料主要的力学和变形特性、塑性区和损伤区变化趋势。基于主应力空间的隐式积分算法所开发的程序可以进行岩土工程问题的数值分析,对现场施工提供指导和理论依据。     

单轴压缩非贯通节理岩石尖端塑性区及弹塑性断裂模型研究

基于Drucker-Prager（下简称D-P）准则,建立压缩载荷作用下的非贯通节理岩石的弹塑性断裂模型。针对节理岩石小范围屈服翼裂纹尖端塑性区,推导了D-P屈服准则的纯I、纯II及I、II复合型3种翼裂纹无量纲塑性区径长函数,并与Mises准则的塑性区进行对比;结果表明,D-P准则的I型和复合型塑性区较Mises屈服准则的塑性区大,且其II型及I、II复合型塑性区在翼裂纹上下表面不连续。进一步,引入断裂软化因子以表征节理岩石裂隙断裂扩展后的断裂软化规律,考虑非贯通节理岩石复合型断裂软化,是由于节理尖端翼裂纹应变能密度超过最小应变能密度导致其成核扩展引起的,提出用应变能密度的指数函数形式表征断裂软化变量的演化;塑性屈服函数采用Borja等的应力张量3个不变量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;建立节理岩石的弹塑性断裂本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性断裂模型的程序。以单轴压缩下非贯通节理岩石为例,分析岩石断裂韧度、节理摩擦系数和节理倾角等参数的影响,结果表明,所提出的弹塑性断裂模型与数值和试验结果比较吻合。     

岩石弹塑性损伤本构模型建立及在隧道工程中的应用

在实际隧道施工过程中,隧道开挖引起地下岩体应力重分布使得围岩的微裂纹扩展损伤,并伴随有塑性流动变形。在地下水环境中对于孔隙和微裂隙围岩介质受到应力作用时,在内部将产生高孔隙水压力影响岩石的力学性质,也改变了围岩的破坏模式。为了研究损伤引起的刚度退化和塑性导致的流动两种破坏机制的耦合作用,从弹塑性力学和损伤理论的角度出发,同时引入修正有效应力原理来考虑孔隙水压力的作用,建立基于Drucker-Prager屈服准则的弹塑性损伤本构模型;针对该本构模型推导了孔隙水压力作用下弹塑性损伤本构模型的数值积分算法-隐式返回映射算法,分别对预测应力返回到屈服面的光滑圆锥面或尖点奇异处两种可能的情况给出了详细的描述,隐式返回映射算法具有稳定性和准确性的特点;大多数弹塑性损伤模型中涉及参数多且不易确定的问题,采用反分析方法获得损伤参数,解决了损伤参数不易确定的难题;采用面向对象的编程方法,使用C++语言编制了弹塑性损伤本构求解程序,并对所建立的弹塑性损伤模型和所编程序进行了试验和数值两个方面的验证;最后将其在吉林抚松隧道工程中进行应用,模拟了塑性区和损伤区的发展变化。研究结果表明：所建立的弹塑性损伤本构模型能够较好地描述岩石的力学性能、塑性和损伤变化趋势,所编程序能够进行实际工程问题的模拟,对现场施工给予一定的指导。     

基于Coulomb准则的混凝土塑性损伤本构模型及其数值验证

以连续介质不可逆热力学为基础,采用了Mohr-Coulomb屈服准则,提出混凝土塑性损伤耦合的新的本构方程。在该模型中采用了塑性应变 、各向同性损伤标量D作为内变量。这个新的本构关系模型严格满足热力学的基本方程。以不同围压作用下混凝土试件的单轴压缩行为为例,采用开发的程序进行了局部水平上本构模型数值验证。结果表明,模型损伤演化数值结果符合试验趋势。     

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

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

Q1黄土的弹塑性损伤本构模型

在陕西洛惠渠五洞简易盾构法施工现场进行了取样，分析了Q1黄土的应力—应变曲线特性及物理性状，根据连续损伤力学和热力学定律，提出了Q1黄土的弹塑性损伤本构模型。并将模型计算结果与试验曲线进行了比较，模型验证结果表明，该模型能够较好模拟Q1黄土的软化和剪胀性。     

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

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

混凝土损伤塑性本构模型研究

提出了一种新的损伤塑性本构模型,强调了应力三轴比对塑性屈服的影响.采用本研究设计的专门用于本构校验的计算机软件,对上述模型的应力应变加载曲线进行了数值模拟,数值结果与实验结果吻合得比较好.     

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

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

广义Drucker-Prager强度准则

基于摩擦性材料的试验研究,在已有研究成果的基础上,提出了具有Drucker-Prager准则表达形式的广义非线性强度理论,即广义Drucker-Prager强度理论。该强度理论在偏平面上的强度函数是介于Drucker-Prager准则和Mohr-Coulomb准则外包线的光滑外凸曲线,在子午面上的强度函数为线性函数。该强度理论用一个表达式来统一描述材料的强度特性,包含或逼近了现有主要的非线性单一强度理论,并且在应力空间中偏应力和静水压力分离,容易与具体的应力应变模型相结合。     

泥岩弹塑性损伤本构模型及其参数辨识

以连续介质力学和不可逆热力学为基础,将损伤引入到修正的Mohr-Coulomb准则中,建立了泥岩弹塑性损伤本构模型反映泥岩软硬化行为,通过构建损伤势函数导出了泥岩的损伤演化方程,编制了泥岩弹塑性损伤本构模型及其参数反演程序。并根据非排水条件下泥岩三轴试验结果,采用多目标函数优化反分析法获得了泥岩本构模型参数。研究成果表明,所提出的弹塑性损伤本构模型能有效地描述泥岩在不同应力状态下的力学特性。     

内摩擦角对冻融岩石损伤本构模型的影响探讨

基于Weibull随机分布和损伤力学理论,从岩石变形全过程的特点入手,选用Drucker-Prager准则,建立能反映冻融岩石各变形阶段特征的损伤本构模型。通过冻融砂岩常规三轴压缩试验,运用Mohr应力圆,获得不同冻融循环次数下砂岩的黏聚力和内摩擦角;由模型理论曲线与冻融砂岩试验曲线对比分析,验证模型的合理性;通过变动内摩擦角,分析摩擦角对模型是否存在影响。结果表明：基于Drucker-Prager准则获得的理论曲线与试验曲线吻合较好,验证了模型的合理性;内摩擦角对分布变量的影响较大,且是线性关系,但对损伤本构模型没有影响或影响不大。研究成果对岩石本构关系的建立有较好的参考价值。     

基于摩尔-库仑准则的断续节理岩体复合损伤本构模型

针对工程实际中断续节理裂隙岩体的损伤本构模型,假设岩石微元强度服从Weibull随机分布,以摩尔-库仑破坏准则作为描述微元强度的表示方法,推导出细观损伤变量。利用能量和断裂力学理论,综合考虑节理几何特征及力学特性,推导宏观损伤变量计算公式。基于Lemaitre应变等效假设,考虑宏细观缺陷耦合作用,推导出复合损伤变量,建立了基于摩尔-库仑准则的宏细观缺陷耦合作用的断续裂隙岩体损伤本构模型。研究结果表明：（1）采用摩尔-库仑准则作为描述微元强度的统计分布变量建立的损伤模型能够较好地反映岩石内部缺陷分布和变形特征,该模型真实地反映岩石微元强度受应力状态的影响。（2）该模型建立的理论曲线与断续节理岩体的试验曲线吻合较好。（3）节理裂隙岩体宏观损伤变量及峰值强度随节理倾角的变化规律与综合考虑宏细观耦合作用下的损伤变量及裂隙岩体峰值强度随节理倾角的变化规律基本一致。（4）宏细观耦合作用下的等效弹性模量与节理贯通率呈非线性负相关;在节理倾角一定的情况下,损伤变量与节理长度呈非线性正相关;在贯通率较小时,岩体的宏观损伤变量与内摩擦角的关系呈线性负相关变化,贯通率达到一定程度,线性关系变成非线性关系。     

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

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

考虑损伤门槛的统计损伤本构模型研究

在统计损伤理论的基础上,从统计损伤本构模型的一般性公式出发,考虑到岩石的受力状态,提出统计损伤本构模型应考虑损伤门槛的影响,建立了考虑损伤门槛的统计损伤本构方程。以石膏角砾岩的常规三轴试验结果作为实例,进行了验证。理论计算结果与试验结果对比表明,考虑损伤门槛的统计损伤本构模型是更加合理的。     

土体非线弹性-塑性本构模型

把本质上属于亚弹性本构模型,在岩土工程界广为应用的Duncan-Chang模型与服从Drucker-PragerMohr-Column屈服准则的弹塑性本构模型相结合,推出了非线弹性-塑性的组合本构模型,以克服经典的弹塑性模型不能考虑岩土材料在塑性屈服前的非线性行为以及一般的Duncan-Chang模型不适用于应力水平接近于屈服或破坏状态等缺点。所建议的本构模型简单实用,能较好地综合利用工程地质勘察资料和常规的土工试验数据。针对实际边坡具体地层条件的分析表明,该模型及其非线性析算法的正确性和可行性,而且在岩土工程方面具有广泛的实际应用价值。     

基于Zienkiewicz-Pande屈服准则的弹塑性本构模型在FLAC3D中的二次开发及应用

Zienkiewicz-Pande屈服准则是Mohr-Coulomb准则的改进,在p-q子午面和π平面上都是光滑曲线,不存在尖点,不仅在数值迭代计算过程中易于处理,而且在一定程度上考虑了屈服曲线与静水压力的关系以及中主应力σ2对屈服的影响。该准则已应用于多个水电站地下洞室开挖支护计算中,但并未被大型数值计算软件所采用。根据弹塑性力学原理,详细推导了基于Zienkiewicz-Pande屈服准则的弹塑性本构模型在FLAC3D中的增量迭代计算格式;在FLAC3D提供的二次开发程序接口基础上,利用VC++程序语言实现了计算过程,并编译成动态链接库文件进行加载和调用。将自编的本构模型对一理想圆形隧洞进行模拟,计算结果与Mohr-Coulomb准则一致,验证了计算格式的正确性。将此本构模型应用于一大型地下洞室群的开挖围岩稳定分析中,得到的塑性区、位移场、应力场分布与Mohr-Coulomb准则计算结果分布规律相符,进一步验证了模型的适用性和可靠性。