基于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中,实现软化模型的数值求解。最后,对多种岩石的室内试验模拟结果表明,所建模型可正确描述不同围压下硬岩和软岩的峰后应变软化变形及剪胀变形,具备良好的适用性。与定临界塑性应变或定剪胀角的应变软化模型模拟结果相比,所建模型模拟曲线可一致性地与试验曲线接近。     

基于裂隙岩样的多组贯穿裂隙岩体峰后应力-应变曲线研究

在岩土工程中,裂隙岩体经常处于峰后变形状态,研究裂隙岩体的峰后应力-应变关系对预测结构的稳定性有重要意义。基于峰后软化阶段强度参数的逐渐演化行为,首先提出一个求岩石峰后应力-应变关系和裂隙峰后应力-切向位移关系的一般方法。然后采用摩尔-库仑强度准则,以岩石的最大主应变和裂隙的切向位移作为软化参数,假设强度参数为软化参数的分段线性函数,分岩体沿裂隙滑移破坏和沿岩石剪切破坏两种情况,提出多组贯穿裂隙岩体峰后应力-应变关系式的求法。最后,在算例中,分岩体沿裂隙滑移破坏和沿岩石剪切破坏两种情况给出了裂隙岩体的峰后应力-应变曲线,讨论了裂隙的平均间距、法向刚度和剪切刚度对峰后应变的影响。结果表明,裂隙的平均间距、法向刚度和剪切刚度越小,裂隙岩体的轴向应变越大。     

硬岩弹塑性变形破坏过程中强度参数及剪胀角演化模型研究

根据大量试验资料,研究硬岩弹塑性变形破坏过程中强度参数及剪胀角演化规律,并建立相应演化模型。分析26例硬岩临界塑性应变与围压之间的关系,提出一种具有广泛适用性的3参数幂函数型围压函数,并将其引入塑性内变量中。基于Mohr-Coulomb屈服准则,分析塑性强化-软化变形过程中强度参数及剪胀角随塑性内变量的演化规律,并建立高斯型强度参数演化模型和考虑围压及塑性内变量的剪胀角演化模型。考虑到工程应用中往往将硬岩峰前视为线弹性变形,进一步研究提出峰后强度参数及剪胀角演化模型。对比多组硬岩试验结果与模拟结果,发现应用所建塑性强化-软化变形过程中强度参数及剪胀角演化模型,可正确模拟不同围压下硬岩塑性强化-软化变形过程中轴向变形及环向变形规律;应用所建峰后强度参数及剪胀角演化模型,也可合理模拟不同围压下硬岩峰后塑性软化变形过程中轴向变形及环向变形规律。上述力学模型对不同种类硬岩均适用,具有良好的应用前景。     

泊松比对岩样破坏模式及全部变形特征的影响

利用编写的计算岩样全部变形特征的FISH函数, 采用FLAC模拟了泊松比不同时单缺陷岩石试样的破坏及全部变形特征。在峰前及峰后, 本构模型分别取为线弹性模型及莫尔库仑剪破坏与拉破坏复合的应变软化模型。高泊松比使岩样发生由单一剪切破坏向复杂破坏转变、破坏区域的面积增加、剪切带倾角降低, Coulomb、Roscoe及Arthur理论对此无法解释。不同泊松比时计算得到的峰前应力-轴向应变曲线、应力-侧向应变曲线、侧向应变-轴向应变曲线、体积应变-轴向应变曲线的线性阶段与平面应变压缩条件下的线弹性解吻合。若泊松比超过1/3, 通过计算得到的平面应变压缩泊松比可大于0.5, 这被数值模拟确认。泊松比的增加使峰后的侧向应变-轴向应变曲线、体积应变-轴向应变曲线、计算得到的泊松比-轴向应变曲线变得不陡峭, 使峰后的应力-侧向应变曲线变得陡峭, 使破坏的前兆变得不明显。      

多孔隙岩石加卸载力学特性及本构模型研究

为了研究岩石在加载-卸载过程中的应力-应变关系,以砂岩为例,对其进行常规三轴加卸载试验。分析了峰后卸载阶段岩石的非线性特性,对岩石的损伤变量进行定义,给出了峰后卸载过程中用于描述应力-应变关系的弹性模量模型。通过分析加载-卸载过程中的轴向应变与径向应变的关系,得到了卸载过程中泊松比模型。引入D-P塑性模型,针对砂岩的塑性硬化特性,对硬化函数进行修正,建立了与等效塑性应变相关联的损伤模型。将计算模型矩阵化后进行数值计算。在此过程中得到如下结论:多孔隙岩石在加载过程中表现出明显的非线性特征,随着体应力的增大,岩石的弹性模量逐渐增大。岩石峰后卸载过程中,当轴向应力大于围压时,应力-应变可以利用峰前弹性阶段的弹性模量模型乘以连续性因子进行描述。随着等效塑性应变的增大,泊松比先增大后减小,最终趋于稳定。峰后卸载过程中,等效塑性应变不发生变化,此时泊松比保持不变。利用提出的本构模型进行了数值计算,数值计算结果与试验结果进行对比,结果表明,提出的模型能够反映出岩石在峰后卸载过程中的应力-应变规律。     

弹性模量对岩样破坏前兆及全部变形特征的影响

利用FLAC内嵌语言编制的计算平面应变压缩岩样轴向、侧向、体积应变及泊松比的FISH函数,计算了弹性模量不同时单缺陷岩样的全部变形特征,研究了弹性模量对岩样的破坏过程及前兆的影响。在峰前及峰后,岩石的本构模型分别取为线弹性模型及莫尔-库仑剪破坏与拉破坏复合的应变软化模型。当弹性模量不是较高时,岩样自始至终仅出现一条倾斜的剪切带。当弹模较高时,最终剪切应变仅强烈集中于贯通岩样的剪切带内部。剪切带倾角在Arthur与Coulomb倾角之间,且随弹模的增加而降低,经典理论对此不能解释。随着弹模的增加,峰值应力增加,峰前的应力-轴向应变曲线变得陡峭,峰后的应力-轴向应变曲线的斜率几乎不变。随着弹模的增加,峰值强度所对应的轴向应变及侧向应变的大小均降低;应力-侧向应变曲线软化段变得陡峭。随着弹模的降低,侧向应变-轴向应变曲线、泊松比-轴向应变曲线及体积应变-轴向应变曲线在峰前发生转折（偏离线性状态）的程度越来越大;非弹性轴向应变增加;材料缺陷附近的最大剪切应变增量增加。岩样破坏的前兆随着弹性模量的降低而逐渐增强。     

强度参数对初始随机缺陷岩样全部变形的影响

在单轴平面应变压缩条件下, 采用FLAC模拟了初始内聚力及内摩擦角对具有随机材料缺陷岩样轴向、侧向、体积变形及由侧向应变及轴向应变计算得到的泊松比的演变的影响。采用编写的FISH函数于试样内部规定随机缺陷并计算其全部变形特征。密实的岩石服从莫尔库仑剪破坏与拉破坏复合的破坏准则, 破坏之后呈现应变软化—理想塑性行为。随着强度参数(初始内聚力及内摩擦角)的增加, 应力峰值及对应的侧向应变的值提高, 体积应变的峰值增加, 应力峰值所对应的计算得到的泊松比稍有增加。由于被剪切带所分割的毗邻块体之间的相对滑动, 纵然扩容角为零, 试样在峰后的变形阶段, 仍然可以观察到体积膨胀现象, 这与作者的理论分析结果一致。当强度参数降低时, 侧向应变—轴向应变曲线、体积应变—轴向应变曲线及计算得到的泊松比—轴向应变曲线的峰前非线性部分变得短暂。剪切带倾角的数值解的上限低于Coulomb理论, 下限在Roscoe理论附近波动, Arthur理论的预测结果与本文剪切带倾角的数值解更接近.      

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

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

加载速度对岩样全部变形特征的影响

利用编制的计算平面应变压缩岩样轴向、侧向、体积应变及泊松比的FISH函数,采用FLAC模拟了加载速度对剪切带图案及岩样全部变形特征的影响。在峰前及峰后,本构模型分别取为线弹性及莫尔-库仑剪破坏与拉破坏复合的应变软化模型。加载速度较低及适中时,岩样发生单剪切破坏,剪切带倾角及宽度不受加载速度影响,应力-轴向应变曲线及应力与侧向应变曲线软化段的斜率不依赖于加载速度;高加载速度使岩样发生X型剪切破坏,两种曲线软化段较平缓;在相同的轴向应变时,高加载速度使剪切带长度降低。随着加载速度的增加,岩样失稳破坏的前兆越来越明显,当加载速度较高时,前兆反而不明显,这是由于应力存在较大的波动,导致不正确地估计了应力峰值所对应的轴向应变。在应变软化阶段,高加载速度使侧向应变与轴向应变曲线、泊松比与轴向应变曲线及体积应变与轴向应变曲线变平缓,也使体积应变与轴向应变曲线的峰值及对应的轴向应变增加。     

基于裂纹体应变的深埋大理岩加、卸荷破坏机制研究

为了弄清深埋硬岩的加、卸荷破坏机制,系统开展了深埋大理岩的常规三轴试验、保持 不变的卸围压试验和变 的卸围压试验等。以裂纹体应变为主要分析变量,结合体应变、等效塑性应变等参量深入分析上述不同应力路径下硬岩的破坏过程。运用裂纹体应变-轴向应变曲线、等效塑性应变-轴向应变曲线和轴向应力-应变曲线来解释岩石破坏过程所产生的现象与规律。结果表明：在裂纹闭合阶段岩样裂纹闭合的变化程度（裂纹体应变改变的大小）可以反映岩石的初始损伤程度;裂纹闭合阶段对岩石弹性模量的计算有重要的影响,需要根据合适的阶段划分,消除裂纹闭合阶段对弹性模量计算的影响,进而得到较为准确的弹性模量;在卸围压点处岩石的可闭合程度突然增加,裂纹体应变-轴向应变曲线发生突变;岩石卸荷破坏过程中裂纹扩展存在滞后性。研究成果有助于进一步理解深埋硬岩的加、卸荷破坏过程和机制,为深埋隧洞的灾害防治提供理论依据。     

应变软化模拟与圆形隧道衬砌分析

在岩土工程中,有许多类岩体呈现应变软化特性。在应变软化过程中,如果多个材料强度参数分析发生变化,就很难给出软化问题的解析解。已有的文献研究表明,在经典的弹塑性理论框架下,模拟岩土介质的应变软化过程依然值得深入研究。利用应变软化材料强度参数随塑性应变增加而弱化的特点,提出了将应变软化过程简化为一系列的脆性应力跌落与塑性流动过程的思想,应变软化行为的模拟转化为一系列脆塑性过程的分析。将这种模拟应变软化过程的思想方法用于求解均匀初始应力场中应变软化介质内的圆形隧道开挖问题,研究了衬砌对减少岩土开挖扰动的作用以及考虑应变软化对衬砌内力的影响。数值结果表明,对于应变软化岩土材料,忽略其强度软化的行为是偏于危险的。     

煤岩两体模型变形破坏数值模拟

采用拉格朗日元法,在弹性岩石与弹性-应变软化煤体所构成的平面应变两体模型的上、下端面上不存在水平方向摩擦力条件下,模拟了模型的破坏过程、岩石高度对模型及煤体全程应力-应变曲线、煤体变形速率、煤体破坏模式及剪切应变增量分布的影响。结果表明,当模型的全程应力-应变曲线达到峰值时煤体内部的剪切带图案已经十分明显,在模型的应变硬化阶段,煤体中的应变局部化可视为模型失稳破坏的前兆,随岩石高度的增加,模型应力-应变曲线的软化段变得陡峭,这与单轴压缩条件下的解析解在定性上是一致的;煤体应力-应变曲线的软化段变得平缓,煤体消耗能量的能力增强;弹性阶段煤体的变形速率降低;煤体内部的剪切应变增量增加。煤体应力-应变曲线的软化段的斜率、弹性阶段煤体的变形速率、煤体内部的剪切应变增量及塑性耗散能都受岩石高度的影响,说明了岩石几何尺寸对煤体的影响(煤岩相互作用)是不容忽视的。     

冻结粉土三轴压缩变形破坏与能量特征分析

对不同温度和不同围压下的青藏冻结粉土进行了三轴压缩与加卸载试验,得到冻结粉土应力-应变关系曲线、抗压强度等力学参数随温度与围压变化的关系。结果表明,冻结粉土典型应力-应变曲线在低围压下大致可以分为线弹性、峰前塑性变形与峰后软化3个阶段。当? 3 < 3.0 MPa时,应力-应变曲线具有明显的峰后软化现象,随围压的增大,软化现象逐渐减弱,当? 3 达到14 MPa,应变软化现象重新变得明显;冻结粉土的强度与变形模量均随围压的增加先升后降;低围压作用下冻结粉土体积随轴向应变的增加先缩后胀,而高围压下体积变形只有体缩;低围压下冻结粉土体积塑性变形耗散能先是随着体积塑性变形增大而增大,之后由于剪胀而减少,高围压下体积塑性变形耗散能始终增加;剪切塑性变形耗散能与塑性剪应变之间近似成抛物线的关系。