基于Burgers模型考虑损伤的非定常岩石蠕变模型

Burgers模型只能描述岩石蠕变的前两个阶段,为了描述蠕变全过程,考虑蠕变参数的时间效应及损伤带来的影响,采用非定常黏性元件取代Burgers模型中串联的定常黏性元件,使其能描述加速蠕变阶段。根据损伤变量的变化特征,假定了一个函数,经过Lemaitre应变等效原理,代入Kelvin模型可得到能描述蠕变衰减阶段的模型。试验数据拟合结果显示,改进的Burgers模型能很好地描述加速蠕变阶段,模型拟合曲线与试验曲线基本吻合,相关系数高,参数取值也在合理范围内;非定常Burgers模型更适合描述不同应力下的岩石蠕变特征曲线。     

考虑硬化和损伤效应的盐岩蠕变本构模型研究

盐岩是储存化石能源和高放射性核废料的理想介质,研究盐岩蠕变力学特性对地下盐岩储库的安全营运具有重要意义。为合理反映盐岩蠕变过程存在的损伤和硬化两种机制,以元件组合模型为基础,结合分数阶微积分理论建立了一种考虑损伤和硬化效应的盐岩蠕变本构模型。该模型采用考虑时效损伤的弹性元件描述盐岩加载初期的损伤变形,采用分数阶村山体描述盐岩衰减蠕变阶段的黏弹塑性蠕变力学行为;通过引入描述盐岩屈服强度随时间强化的硬化函数来反映盐岩硬化机制;采用瞬时塑性元件描述不可恢复的瞬时变形。基于Kachanov蠕变损伤定律与Lemaitre应变等价原理构造了一种带有应变触发的非线性黏壶元件,该元件能够较好描述盐岩加速蠕变阶段的非线性变形特征。基于组合模型理论推导了考虑硬化和损伤效应的盐岩一维、三维蠕变方程;通过分析已有盐岩单、三轴蠕变试验等时应力-应变曲线特征,确定了村山体启动应力阈值;结合等时应力-应变曲线和蠕变试验数据对模型中的参数进行了辨识。结果表明：建立的蠕变本构模型仅用一组参数即能描述不同应力状态下盐岩的蠕变力学特性,可为预测盐岩蠕变变形特征提供一定的理论依据。     

基于硬化-损伤机制的贫煤蠕变本构模型研究

岩石蠕变是其内部硬化效应和损伤效应共同作用演化的结果。通过分级加载下贫煤单轴压缩蠕变试验,基于试验中瞬时弹性模量增大和黏滞系数减小的现象,分析贫煤蠕变中的硬化-损伤机制。结果表明,较低应力下煤样的变形仅为瞬时应变,煤样仅有硬化效应;当应力达到蠕变起始应力阈值时煤样逐渐出现损伤,次生裂纹出现并扩展,硬化效应与损伤效应共同作用下煤样表现出衰减蠕变和稳定蠕变特征。长期蠕变过程中煤样强度经历先硬化再弱化的过程,最终因累计损伤过大而失稳破坏,依据各级应力下煤样的变形特征,构建了符合其变形规律的黏弹塑性组合模型,并引入了硬化函数和损伤函数,推导得到其蠕变本构方程。采用该模型描述分级加载下贫煤单轴压缩蠕变变形曲线,试验曲线和模型曲线吻合良好。该模型元件简单,物理意义明确,能很好地反映贫煤单轴压缩下的硬化–损伤蠕变特性。     

考虑温度损伤的盐岩蠕变本构关系研究

根据统计力学原理,以分形岩石力学为桥梁,对盐岩在温度与应力耦合作用下蠕变特性进行了研究,导出了考虑围压效应的损伤变量表达式。综合前人的成果,推导出温度–应力耦合下的盐岩损伤方程。对盐岩在温度与三轴应力共同作用下的蠕变试验数据进行研究,分析了盐岩蠕变过程中衰减蠕变阶段、稳态蠕变阶段和加速蠕变阶段的各自特征。对于广义Bingham 蠕变模型在衰减和稳态蠕变阶段引入一非线性函数,在加速蠕变阶段引入损伤,建立了盐岩考虑温度损伤的蠕变本构关系。利用盐岩的三轴蠕变试验数据对新的蠕变损伤模型参数进行辩识与比较分析。结果表明,该模型能够很好地反映在不同温度作用下盐岩的衰减蠕变阶段、稳态蠕变阶段和加速蠕变阶段的规律,验证了该模型的合理性和正确性。     

一种基于修正西原模型的冻土蠕变本构关系

冻土在高应力水平下通常会表现出加速蠕变特性。经典西原（Nishihara）模型可以反映衰减蠕变和稳态蠕变,但不适用于描述加速蠕变。通过对Nishihara模型中黏性元件进行修正,将黏弹性元件中的定常的黏滞系数修正为时间的函数,并在黏塑性元件中采用非定常、非牛顿黏性元件,得到的蠕变柔量可以考虑应力水平和时间的影响,使其能够反映其加速蠕变的过程。用提出的模型分别对一维和三维应力、不同应力水平下的冻土蠕变试验数据进行了拟合,拟合结果能够反映冻土从衰减蠕变和稳定蠕变状态进入加速蠕变的过程,说明了该模型对于描述冻土在不同应力下不同蠕变状态的适用性。     

深部岩体多孔介质流变模型的研究

考虑到传统流变模型在反映载荷水平影响和描述加速蠕变阶段方面的缺陷,在改进传统流变元件模型的基础上,引入非理想黏性元件,建立了在数学上能体现载荷水平影响并能描述加速蠕变的新的改进流变模型;同时从定义出发,考虑到孔隙度的动态性,对孔隙度和通用有效应力公式进行修正;鉴于岩体的孔隙性和不均匀性,把岩体看作多孔介质,应用多孔介质理论中的有效应力代替改进岩土流变模型中的恒应力,计入温度和孔隙压的影响,最终建立了深部（高温度、高地压和高孔隙压）岩体多孔介质流变模型。     

深部人工冻土黏弹塑损伤耦合本构模型研究及其验证

采用改进的Mohr-Coulomb屈服准则(Zienkiewicz-Pande准则), 建立了能描述人工冻土3个阶段蠕变特性的黏弹塑损伤耦合本构模型. 基于该本构模型建立了模型的柔度矩阵, 将本构模型嵌入到ADINA有限元程序中. 计算表明: 数值模拟得到的人工冻土蠕变特性曲线与实测值吻合良好, 最大误差为4.8%, 说明基于改进的Mohr-Coulomb准则建立的黏弹塑损伤耦合本构模型能够准确的描述人工冻土剪切特性和3个阶段的蠕变特性.获得的本构模型具有参数少、容易确定等优点. 深部人工冻土黏弹塑损伤耦合本构模型的建立, 对于人工冻土结构物长期稳定分析和预测具有重大的理论指导意义.     

基于断裂及高温损伤的岩石蠕变模型研究

为了反映热-力耦合作用下岩石蠕变变形的全过程,依据断裂力学原理提出了岩石裂纹扩展的临界损伤应力和一个新的可描述岩石在稳态蠕变阶段与临界损伤应力相关的非线性黏性分量,在传统西原模型和Burgers模型的基础上,将指数形式的损伤变量、临界损伤应力以及与其有关的非线性黏性分量引入到流变微分方程,通过叠加原理推导了考虑温度效应的单轴和三轴压缩条件下岩石的流变本构关系,建立了岩石的热-力耦合损伤蠕变本构模型。利用不同温度、不同应力条件下花岗岩的三轴蠕变试验曲线和本文蠕变模型的计算曲线进行比较,结果表明本文蠕变模型能够较好地模拟岩石在初始瞬态、稳态和加速蠕变阶段全过程的变形规律,验证了所建模型的有效性和合理性。该模型为分析高温、高应力环境下岩石工程的长期变形和稳定情况提供了理论依据。     

人工冻结黏土广义开尔文蠕变本构模型模糊随机优化

冻土的蠕变本构模型是冻结法建井的重要参考资料。在-5℃、-10℃、-15℃和-20℃下,针对人工冻结黏土进行不同应力等级的单轴蠕变试验,获得不同温度、不同应力等级条件下蠕变曲线的变化规律。通过模糊随机迭代搜索改进传统的Guess-Newton算法,得到模糊随机Guess-Newton算法步骤,利用改进的算法优化广义开尔文本构模型参数,获得各温度和应力下的优化蠕变模型。工程试验表明:模糊随机Guess-Newton算法能有效地优化广义开尔文本构模型参数,使模型更好地拟合冻土各阶段的蠕变值,准确地表征了冻土的蠕变特性,同时改进后的算法较传统的算法效率更高,收敛更快。      

基于硬化参量的盐岩蠕变疲劳本构模型

盐岩具有较好的蠕变特性和损伤自恢复性,被公认为储能或油气储备的理想介质。对盐岩复杂力学行为进行准确表征和预测是保障盐穴地下空间利用工程安全性的基础。通过定义硬化参量等特征因子,建立了新的能够考虑复杂加卸载路径的盐岩蠕变疲劳本构模型。基于盐岩变形位错机制,引入双曲线阻尼元件,作为表征岩石硬化程度的状态变量;通过硬化参量的演化,考虑加卸载历史对盐岩变形行为的影响。借鉴经典Norton模型的应力-应变关系建立蠕变疲劳本构的基本框架关系。通过引入裂隙扩展因子,假设初始形核长度,考虑材料的断裂韧度,修正邻近破坏阶段(加速变形阶段)的应力-应变关系。该模型能够很好地预测常规蠕变、循环加卸载、下限间隔循环加卸载、梯形波蠕变循环加卸载等常见复杂加卸载路径下的塑性变形特征,且能够较好地表征恒荷载蠕变与循环加卸载之间的相互影响。新的蠕变疲劳本构模型中大部分参数具有较为明确物理意义,参数a表征盐岩稳态变形阶段应力与变形率的关系因子,参数b为确定盐岩第一阶段减速变形阶段的关系因子,参数d₀和μ_d分别表示初始裂纹形核量和裂隙扩展速率因子,共同修正模型临界破坏阶段的应力...     

Nonlinear Variation Parameters Creep Model of Rock and Parametric Inversion

Traditional Nishihara model is a linear constant parameters model, and it is difficult to describe nonlinear accelerated creep characteristics of rock with this model. Therefore, an improved Nishihara model is proposed in the way of replacing, cascading elements and substituting parameters in the original model. Firstly, fractional differential elements are used to replace viscous elements in the traditional Nishihara model; then a viscous plastic element with assumption exponential relationship between the rock nonlinear rheological deformation and time power function is cascaded; finally, the viscosity coefficient of the traditional model is rewritten to a variable parameter related to stress and time, thus obtaining the nonlinear variation parameters creep model of rock. The improved model is verified with data from creep test to frozen soft rock during Cretaceous period, which shows that this model can describe the rock nonlinear accelerated creep characteristics more accurately and more comprehensively.     

基于分数阶导数理论的Q3黄土流变本构模型研究

荷载的时间效应诱发的黄土蠕变和渐进破坏会显著影响其工程构筑物的长期服役性能。特别是在我国西北地区黄土广泛分布,在荷载的长期作用下一些典型工程建构筑物(如西安地铁开挖的黄土隧道工程)可能因黄土蠕变和渐进破坏而失稳,进而对工程建设的顺利实施以及工程建设完成后其安全运营维护产生极为不利影响。可见,对黄土流变特性进行分析并揭示其流变机理进而建立其本构模型具有重要的理论和实际意义。为此,以Q₃黄土为研究对象,基于分数阶导数理论,通过对典型蠕变变形过程进行分析,提出了一种可模拟加速蠕变阶段非线性变化特征的分数阶元件模型,对其进行理论分析进而建立了分数阶改进西原模型并推导出本构方程。在此基础上对Q₃黄土的原状与重塑样开展了三轴分级循环加-卸载流变试验,以验证分数阶改进西原模型的有效性。相应的理论计算结果表明：分数阶改进西原模型既能实现对Q₃黄土减速蠕变、等速蠕变和加速蠕变3个流变阶段的模拟,又可弥补整数阶改进西原模型无法描述加速蠕变阶段的不足,而且相比于整数阶改进西原模型,其在减速蠕变及卸载阶段的预测效果更佳。研究结果为考虑黄土流变问...     

煤岩蠕变模型与破坏特征试验研究

在MTS岩石伺服试验机上对山东某矿3煤进行了短时流变试验,并对试验结果进行了分析。结果表明：该煤岩蠕变门槛系数为0.915、流变系数为0.383;煤岩流变破坏呈现空间不均匀性,整体以塑性破坏为主,但存在剪切主导破坏面;煤岩流变过程中有微破裂破坏即损伤存在,并且这种损伤在较低的应力水平就有所发生,大量微破坏的积累与贯通导致了煤岩试件破坏;煤岩蠕变过程中不仅存在损伤软化,而且存在蠕变硬化,蠕变的过程取决于这两种机制竞争的综合结果;可以用3次多项式经验蠕变模型较好地描述该煤岩的蠕变特性,西原模型可以较好地描述该煤岩的初始蠕变和等速蠕变阶段,但不能描述加速蠕变阶段。     

粘弹塑本构模型及用于冻土数值计算的柔度矩阵

随着寒区经济发展及冻结法施工技术在工程中的推广应用,对冻土蠕变的深入研究显得尤为重要.对冻土单轴流变试验结果分析,用非线性牛顿体替代线性牛顿体,从而对西原模型进行了改进,得到了冻土三向应力状态下的本构方程;并利用D-P屈服准则推导出应用与数值计算的柔度矩阵.通过有限元程序的二次开发,将本构模型添加到大型非线性程序ADINA中.用三轴蠕变试验检验了粘弹塑本构模型,获得相应的试验拟合曲线.结果表明:三轴蠕变试验与理论计算结果吻合良好,说明粘弹塑本构模型可以描述冻土在高应力下反映的非衰减性蠕变变形特征.     

基于西原加速模型的煤体蠕变特性试验

为克服传统西原模型无法描述煤体加速蠕变特征的缺陷,针对煤体加速蠕变阶段中黏-弹-塑性流动特征,利用一个非线性的黏-弹-塑性模型来代替传统西原体模型中的理想黏-塑体,得到改进的西原加速模型,并推导该模型在等围压三轴条件下的蠕变方程。通过开展煤样三轴蠕变试验,采用模型辨识的方法,利用改进的西原加速蠕变模型对煤体蠕变试验数据进行了拟合分析。结果表明,拟合结果与试验数据基本吻合,较好地反映出了煤体非衰减蠕变特征;并对改进模型的稳定性进行讨论,加速蠕变模型中的弹性模量E₃具有明确的物理意义,可将E₃<0作为煤体发生蠕变失稳的判断依据。      

抛物线型屈服面人工冻土蠕变本构模型研究

利用自行研制的W3Z-200型冻土三轴压缩蠕变仪,对淮南地区深表土冻土进行了三轴蠕变试验。根据试验结果,提出以标准黏弹塑性模型为基础,并耦合温度自由度,应用抛物线型屈服准则对模型中的黏塑性项进行改进,建立了新的本构模型来描述高围压复杂应力状态下冻土的蠕变变形特征。基于ADINA有限元软件平台,开发了冻土本构子程序。对淮南某冻结井筒开挖过程进行了数值模拟,计算结果验证了冻土抛物线型黏弹塑本构模型的正确性和合理性。     

人工冻土黏弹塑性本构参数反分析研究

通过大量的人工冻结黏土蠕变试验研究,获得在高偏应力下人工冻结黏土具有显著的蠕变特性：当应力水平较低时,人工冻土发生黏弹性蠕变变性;当应力水平较高时,人工冻土发生黏塑性蠕变变性。因此,以西原模型为基础,增加一个黏弹性元件项和用抛物线屈服函数代替塑性屈服项,构建了高应力下抛物线型黏弹塑性蠕变本构模型。改进的西原模型可以较好描述人工冻土剪切蠕变特性。基于小生境技术的遗传算法和人工冻土黏弹塑性模型有限元理论,在FEPG（有限元程序生成器）自动生成平台上研发了位移反分析程序。利用现场实测数据对深井冻结井筒开挖过程进行了位移反分析,获得了与试验相一致的冻土本构力学模型参数,数值分析结果与实测数据之间的误差在6%以内,反演结果验证了模型的可行性。冻结井筒开挖过程的位移反分析获得的人工冻土本构模型参数,能反映冻结壁施工过程和有效冻结范围内冻土本构参数平均值,因而具有广泛的代表性,对冻结壁应力场、位移场和稳定性预测具有重要工程意义。     

冻粘土蠕变损伤耦合本构关系研究

以冻土试验为基础,提出了冻粘土服从DP屈服准则的粘塑损伤变量,通过推导得到了相关联流动法则下的冻土粘弹塑性损伤耦合本构方程.为了研究粘弹塑性和损伤的耦合作用,将建立的冻土材料粘弹塑性损伤耦合本构模型,通过用户子程序嵌入到ADINA有限元中.结果表明:粘弹塑损伤耦合本构模型能较好的模拟冻粘土加速蠕变过程,与实测蠕变吻合较好.