基于Hoek-Brown准则的岩石应变软化模型研究

基于Hoek-Brown（H-B）准则,根据大量试验数据,研究提出一种具有广泛适用性的岩石应变软化模型。首先,利用40种硬岩和软岩三轴试验数据,验证了3参数幂函数型围压函数广泛适用于硬岩和软岩峰后临界塑性应变与围压的关系描述,并将其引入塑性内变量中;其次,比较分析当前基于H-B准则的岩石应变软化模型中软化参数的合理性,明确以GSI作为软化参数的方法最优,拟合结果正确且参数最少。然后,研究分析岩石峰后GSI随塑性内变量的演化规律,并建立一种普适的非线性演化模型。与其他演化模型相比,所建模型可更准确描述不同种类岩石峰后H-B强度参数的非线性演化规律,同时参数较少。接着,采用非关联流动法则,将所建GSI演化模型嵌入有限元软件ABAQUS中,实现软化模型的数值求解。最后,对多种岩石的室内试验模拟结果表明,所建模型可正确描述不同围压下硬岩和软岩的峰后应变软化变形及剪胀变形,具备良好的适用性。与定临界塑性应变或定剪胀角的应变软化模型模拟结果相比,所建模型模拟曲线可一致性地与试验曲线接近。     

基于Hoek-Brown准则的岩石应变软化模型研究

基于Hoek-Brown（H-B）准则，根据大量试验数据，研究提出一种具有广泛适用性的岩石应变软化模型。利用40种硬岩和软岩三轴试验数据，验证了3参数幂函数型围压函数广泛适用于硬岩和软岩峰后临界塑性应变与围压的关系描述，并将其引入塑性内变量中；其次，比较分析当前基于H-B准则的岩石应变软化模型中软化参数的合理性，明确以GSI作为软化参数的方法最优，拟合结果正确且参数最少。然后，分析岩石峰后GSI随塑性内变量的演化规律，并建立一种普适的非线性演化模型。与其他演化模型相比，所建模型可更准确地描述不同种类岩石峰后H-B强度参数的非线性演化规律，同时参数较少。接着，采用非关联流动法则，将所建GSI演化模型嵌入有限元软件ABAQUS中，实现软化模型的数值求解。最后，对多种岩石的室内试验模拟结果表明，所建模型可正确描述不同围压下硬岩和软岩的峰后应变软化变形及剪胀变形，具备良好的适用性。与定临界塑性应变或定剪胀角的应变软化模型模拟结果相比，所建模型模拟曲线可一致性地与试验曲线接近。     

软岩应变软化数值模型的建立与研究

根据三轴压缩试验的应力-应变曲线,运用Hoek提出的由主应力圆包络线确定粘聚力和内摩擦角等效值的方法和曲线拟合的方法,研究了围压对软岩峰后软化特性的影响和软岩的宏观物性参数峰后应变软化规律,建立了峰后应变软化数值模型。利用该数值模型对三轴试验进行数值模拟,其结果与试验数据吻合。     

基于改进Bingham模型的软岩参数非定常三维非线性黏弹塑性蠕变本构研究

传统线性元件组合模型难以描述岩石非线性加速蠕变阶段特性,给工程应用带来了巨大困难。提出一种新的应变触发非线性黏壶元件,并将其与Bingham模型串联,建立了可以描述岩石等速和加速蠕变特性的一维本构模型,进而推导出相应的三维蠕变本构模型,再将三轴压缩过程中岩石弹性模量的衰减方程引入该三维蠕变本构方程,得到一个能反映岩石蠕变全过程特性的三维非线性黏弹塑蠕变本构模型,根据泥岩蠕变试验数据对该模型参数进行了辨识。模型理论计算值与泥岩蠕变试验数据的对比表明改进的Bingham模型不仅可以充分反映软岩的等速蠕变特性,还可以很好地描述软岩的衰减和加速蠕变规律,且模型元件少,组合形式简单,为软岩非线性蠕变本构模型研究提供了新的借鉴。     

三轴压缩下粉砂质泥岩蠕变本构模型研究

依据粉砂质泥岩三轴压缩蠕变试验结果,分析了岩石的蠕变力学特性。选取线性黏弹性Burgers模型来描述岩石的衰减蠕变和稳定蠕变特性。引入非线性黏塑性元件,将其与Burgers模型串联起来,建立一个新的六元件非线性黏弹塑性Burgers蠕变本构模型。采用非线性Burgers模型来描述岩石的加速蠕变特性。应用Levenberg-Marquardt算法对蠕变试验曲线进行辨识,得出各级应力水平下岩石轴向与径向的蠕变模型参数。对比试验曲线和模型拟合曲线,二者吻合较好,表明线性Burgers模型与非线性Burgers模型可以准确描述岩石3个蠕变阶段的力学特性。模型参数辨识结果表明,在同一级应力水平下岩石轴向与径向的蠕变模型参数各自独立并不统一,模型参数取值时应考虑岩石蠕变的各向异性以及轴向与径向蠕变量、蠕变速率大小的差异。在前7级应力水平下,模型参数应取以轴向应变辨识得到的Burgers模型参数;在第8级应力水平下,模型参数应取以径向应变辨识得到的Burgers模型参数;在破裂应力水平下,模型参数应取以径向应变辨识得到的非线性Burgers模型参数。     

软弱破碎围岩隧道炭质页岩蠕变特性试验研究

采集贵广铁路天平山隧道典型软弱破碎围岩岩样,在室内开展三轴常规压缩和三轴压缩蠕变试验。基于岩石试验成果,分析炭质页岩的蠕变力学特性,低应力水平条件下岩石仅出现减速蠕变和等速蠕变阶段,在更高的应力水平条件下,岩石并未出现加速蠕变的现象,但围压在一定程度上提高岩石的屈服强度,弹性模量亦有随着围压增加而增大的趋势。据此,建立依托工程炭质页岩蠕变全过程的黏弹塑性应变软化蠕变力学模型,推导该蠕变力学模型在三维应力状态下的本构方程,并确定模型的参数,运用最小二乘法对模型的参数进行辨识。通过与试验曲线对比显示,所建蠕变模型能够描述依托工程炭质页岩的蠕变力学特征,可以用来研究该隧道软弱破碎围岩的蠕变性质和稳定性。     

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

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

模拟应变软化岩石三轴试验过程曲线

在经典弹塑性理论框架下,认为岩石材料在应变软化过程中遵守摩尔-库仑强度准则,塑性变形服从非关联流动法则。对两组岩石在三轴试验中呈现的应变软化行为进行了分析,获取了峰值前后强度参数的大小。采用应变软化模拟方法,将应力-应变曲线峰后应变软化阶段简化为一系列应力跌落和塑性流动的脆塑性过程,在给定的强度参数演化规律基础上,得到了岩石的应力-应变曲线、侧向应变-轴向应变曲线、体应变-轴向应变曲线和塑性体应变-轴向塑性应变曲线。模拟得到的应力-应变曲线与试验结果吻合较好,其他曲线与相关试验测试规律基本一致。研究成果对确定岩石强度参数与认识岩石峰后强度演化规律有一定的指导意义。     

膨胀土卸荷蠕变特性及其非线性蠕变模型

鉴于膨胀土滑坡往往表现为长期性、渐进性等与时间相关的特性,利用GDS应力路径三轴仪对膨胀土进行了三轴卸荷蠕变试验。试验结果表明：当偏应力较小时,膨胀土的蠕变曲线仅出现瞬时变形和衰减蠕变;当偏应力达到一定值时,其蠕变曲线也呈现衰减蠕变、稳态蠕变和加速蠕变3个阶段,但其加速蠕变阶段的特征与一般岩土体不同,其蠕变速度近乎常数。膨胀土的应力?应变等时曲线显示,膨胀土卸荷蠕变具有非线性特征,且其非线性程度与蠕变时间和应力水平相关,蠕变时间越长、应力水平越高,其非线性程度越高。基于非线性流变力学理论,提出了一种非线性四元件蠕变模型,将标准线性体与一个非线性黏壶串联,该模型可描述等围压三轴压缩应力状态下膨胀土轴向应变随时间的演变规律。根据膨胀土卸荷蠕变试验结果,采用曲线拟合法对三维非线性模型的参数进行反演识别。拟合曲线和试验曲线对比显示,两者吻合良好,说明该模型可以很好地描述膨胀土的蠕变特性。此外,基于该蠕变模型获取了膨胀土的临界破坏应力,其与常规剪切破坏应力的比值随着固结压力的减小而减小,表明越接近坡面的土层越容易发生蠕变破坏。     

岩石应变软化变形特性及损伤本构模型研究

对滇藏铁路三江段的花岗岩进行了单轴和常规三轴压缩试验,研究了岩石在不同围压条件下的变形特性及其损伤演化规律。试验结果表明:该岩石压缩变形过程中,内部裂隙的演化经历了无起裂、起裂并稳定扩展、加速扩展、扩展变缓以及扩展停止5个阶段,且围压对岩石变形特性有显著影响。在此基础上,提出了基于有损微元生成的损伤演化模型,建立了三轴压缩条件下岩石应变软化损伤本构模型。该模型能够统一描述岩石峰前及峰后全部变形特性及其损伤演化规律,且损伤方程形式简单、参数物理意义明确。理论结果和试验数据具有很好的一致性,且参数分析结果也表明了该模型的合理性。     

以屈服接近度分段函数表示的非线性流变模型的程序实现

基于三轴流变试验提出的非线性黏弹塑性流变模型存在蠕变门槛值、长期强度和非线性黏塑性体等组件。这种复杂的流变模型在进行程序实现时存在着诸多问题,比如一维的蠕变门槛值和长期强度推广至三维时难于表达、非线性黏塑性部分难于处理等问题。通过引入屈服接近度参数作为蠕变方程分段函数的判别准则,解决了三维情况下蠕变门槛值和长期强度的表达问题;详细推导了黏塑性部分的表达形式,并给出了三维中心差分形式,实现了流变模型在FLAC3D中的二次开发。通过算例验证了程序开发的正确性,对流变模型进行参数辨识,拟合曲线与试验曲线吻合较好,表明了所建模型的正确性与合理性。     

A frozen soil creep model with strength attenuation

A frozen soil creep model with shear strength attenuation was proposed based on a soft soil creep model. K ₀ compression and triaxial shear tests were conducted to obtain the model parameters. With triaxial creep tests, the performance of soft soil creep model was verified at first. It was shown that different creep stages can only be described separately, which was due to a constant failure line used in the model. After incorporating strength attenuation into the original model, relative position between stress point and failure line is time dependent, and the progressive development of creep strain from primary to tertiary stage can be captured reasonably.     

A viscoplastic model including anisotropic damage for the time dependent behaviour of rock

This paper presents a new constitutive model for the time dependent mechanical behaviour of rock which takes into account both viscoplastic behaviour and evolution of damage with respect to time. This model is built by associating a viscoplastic constitutive law to the damage theory. The main characteristics of this model are the account of a viscoplastic volumetric strain (i.e. contractancy and dilatancy) as well as the anisotropy of damage. The latter is described by a second rank tensor. Using this model, it is possible to predict delayed rupture by determining time to failure, in creep tests for example. The identification of the model parameters is based on experiments such as creep tests, relaxation tests and quasi‐static tests. The physical meaning of these parameters is discussed and comparisons with lab tests are presented. The ability of the model to reproduce the delayed failure observed in tertiary creep is demonstrated as well as the sensitivity of the mechanical response to the rate of loading. The model could be used to simulate the evolution of the excavated damage zone around underground openings. Copyright © 2005 John Wiley & Sons, Ltd.     

针对砂土应变软化强非线性问题的动态松弛有限元法研究

对于密实砂土峰值后应变软化特性引起的材料强非线性问题,利用常规有限元所采用的隐式算法很难求解,而动态松弛法在求解这类非线性问题方面却具有独特的优势。针对砂土应变软化强非线性问题,将动态松弛法与有限单元法相结合,提出了一种新的动态松弛有限元法。该方法根据动态松弛法的显式特性,由中心差分法导出了动态松弛有限元法的基本控制方程,并实现了对应力-应变空间中整个平衡路径的追踪。将动态松弛法嵌入到非线性有限元程序中,即可对砂土材料应变软化引起的强非线性问题进行有限元数值计算。该有限元程序中,应力更新采用回归映射算法。最后通过对砂土平面应变压缩试验进行有限元模拟,验证了动态松弛有限元法在求解材料强非线性问题方面的优越性。     

考虑含水率的红层泥岩蠕变特性及改进伯格斯模型

为揭示西南地区广泛分布的红层泥岩在不同含水条件下的蠕变特性,采用YSJ-01-00岩石三轴流变试验机进行三轴压缩蠕变试验。试验结果表明,（1）随着含水率的升高和轴向荷载的提升,蠕应变量与瞬时应变量的比值整体上都呈增大趋势,比值范围为0.117～0.985,平均比值维持在0.348～0.480之间;（2）红层泥岩的初始、稳态和极限加速蠕变速率都随含水率的升高而增大;（3）红层泥岩的长期强度随含水率的上升而降低,量化比较发现红层泥岩的非线性蠕变程度随含水率的上升而越得到体现。研究了红层泥岩蠕变的非线性特征,提出一个黏弹性模量 ,并引入一个非线性黏塑性体,建立一个新的改进后的伯格斯（Burgers）模型。基于BFGS算法和通用全局优化法对蠕变试验曲线进行辨识,确定模型参数,并进行模型可靠性分析,验证了所建模型的正确性和合理性。