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

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

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

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

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

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

考虑应变率效应的广义Hoek-Brown动态强度准则

岩石强度特性具有明显的围压效应和应变率效应,研究复杂应力状态下动态强度准则对于岩石动力学分析至关重要。利用10~(-5)~10~(2.32) s~(-1)应变率范围内砂岩静动态间接拉伸和单、三轴压缩强度试验数据,进行了广义Hoek-Brown强度准则预测砂岩动态强度的适用性研究,分析了材料参数的应变率效应。结果表明,广义Hoek-Brown强度准则能够较好地描述同一应变率下较低拉应力区到较高压应力区砂岩强度规律,低应变率区单轴抗压强度c_s和参数a变化不大,高应变率区随应变率的增加非线性增大,而参数m几乎不受应变率影响。基于分析建立了单轴抗压强度c_s和参数a的应变率效应函数关系,提出了由低应变率到高应变率下的广义Hoek-Brown动态强度准则统一表达式,进一步探讨了岩石强度依赖应变率效应和围压效应的物理机制。     

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

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

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

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

材料缺陷对岩样变形局部化影响的数值模拟

在平面应变压缩条件下,采用拉格朗日元法研究了材料缺陷对岩样应变局部化及宏观力学行为的影响。在数值计算中,采用了莫尔-库仑与拉破坏复合的破坏准则,峰后岩石的本构关系为线性应变软化。对于理想岩样（不含任何缺陷）,从剪切应变率的等值线图及变形网格图发现,变形场相对试样垂直对称轴对称。在试样的边界上设置材料缺陷后,变形场的对称性被打破。材料缺陷附近是局部化现象启动的主要位置。与不含材料缺陷试样相比,含缺陷试样的局部化提前启动。当缺陷位于试样左边界中部附近时,试样内部出现多条剪切带,发生韧性剪切破坏,峰后应力-轴向应变曲线和应力-侧向应变曲线均倾向于韧性,岩样的稳定性增强;当缺陷位于试样上、下端面附近时,仅出现一条贯穿试样左右边界的剪切带,发生脆性剪切破坏,峰后两种曲线倾向于脆性,易发生失稳破坏。从数值结果中还观测到了剪切带的跳跃或迁移现象及剪切带的相互竞争、此消彼涨的现象。当岩样中不包含任何缺陷时,试样应力-变形曲线的峰值强度最高。随着缺陷的上移,峰值强度下降,直到峰值强度基本保持不变。缺陷越接近于岩样的下端,对岩样应力-变形曲线的峰值强度的影响越小。     

考虑围压影响的岩石峰后应变软化力学模型

试验表明,随着围压增加,岩石峰后残余强度和割线模量增加,如何模拟围压对岩石峰后应变软化力学行为的影响是岩石峰后力学行为研究的关键。利用Z. Fang[1]提出的峰后强度下降指数描述围压对岩石峰后残余强度和割线模量的影响,与摩尔-库仑模型结合,建立了考虑围压影响的岩石峰后应变软化力学模型。在FLAC环境下,利用Fish函数方法开发了相应的数值计算程序。数值算例研究了Tennessee大理岩在不同围压下的应变软化过程,模拟结果与试验数据基本一致,表明建立的模型合理,可以应用于岩石峰后软化过程的模拟。     

岩石材料应变软化模型及有限元分析

本文应用弹塑性理论,根据岩石的应力应变关系、峰值强度、残余强度和岩石材料的脆性特性,建立了理想的应变软化模型,并提出了相应的加载准则。文中提出的模型计算方便,所用参数简单,物理意义明确,便于编制有限元程序。     

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

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

椭圆-抛物线双屈服面模型参数灵敏度分析

利用两种土石料三轴CD试验数据确定两组椭圆-抛物双屈服面模型参数,研究了模型参数? 、KG、n、h、m、M1、M2和pr分别单独变化（而其他参数保持不变）对模型所反映的土的应力-应变关系的影响。结果表明,土体的体积应变对参数h、m比较敏感;? 、KG、n对剪胀性土的体变影响较大,而对剪缩性土体变影响不明显;? 、KG、n、h和m对模型反映的强度没有影响;pr变化对体变和强度影响均较小。M1和M2对强度和体变均有显著影响,且M1和M2中的小值决定模型的强度。定义了灵敏度,即一种用于表示模型参数对其所反映的应力-应变的影响程度的参数。分析表明,不同的参数,其灵敏度差异很大;不同的土类,参数灵敏度差异也较大。实际应用时,对灵敏度大的参数取值应特别慎重。     

高地应力下硬岩的本构模型研究

传统的本构模型在模拟高地应力下硬岩破坏的范围和深度方面并不理想。针对高应力下的拉西瓦花岗岩,通过真三轴压缩试验模拟了开挖时应力路径的演化。在摩尔-库仑强度准则、格里菲斯强度准则、德鲁克-普拉格准则、双剪统一强度理论和霍克-布朗强度准则的基础上,考虑十二面单元体主剪面上主剪应力、正应力和静水压力的共同作用,提出了一个有关高应力下硬岩的三剪强度准则。根据试验结果,利用遗传算法全局寻优的功能,搜索出了三剪强度准则的参数。预测样本和与摩尔-库仑强度准则、德鲁克-普拉格准则、双剪统一强度理论的比较表明,该三剪强度准则与试验结果吻合,适用于高地应力下的硬岩。根据试验得到的应力-应变关系曲线,建立了基于三剪强度准则和应变软化的弹-脆-塑性本构模型,并用FLAC3D提供的基于C++的用户自定义模型工具UDM,创建了用户自定义模型的动态链接库,嵌入了FLAC3D软件。采用遗传算法-FLAC方法,搜索得到了基于岩石试样和岩体的本构模型参数。计算结果与实测情况吻合,表明建立的本构模型适用于高地应力下的硬岩,为高地应力下硬岩地下工程的安全性和稳定性分析奠定了基础。     

基于一种脆性指标确定岩石残余强度

岩石的残余强度是岩石力学的重要指标,准确地评价岩石残余强度对于评价地下工程的稳定性以及优化岩体支护设计具有重要意义。基于岩石的三轴力学特性提出一种表征岩石峰后强度衰减行为的力学指标--岩石强度衰减系数,该指标可反映岩石的脆性程度,并提出岩石强度衰减系数与围压关系的幂函数模型。对22组不同成因的岩石常规三轴压缩试验数据进行幂函数模型参数拟合,发现不同岩石拟合所得参数离散性较大,分析其原因主要与岩石矿物组成和岩石结构特征等因素相关。在此基础上提出基于强度衰减方法确定岩石残余强度的方法,分析表明,该方法能够很好地拟合岩石残余强度试验数据,并能反映岩石结构性质对残余强度的影响。     

岩体应力应变曲线的围压及孔压效应和降雨触发滑坡的机制分析

从宏观分析与微观机制方面探讨了围压、孔隙水压力对岩体应力应变曲线峰后特性即脆－延性转变的影响。定义了一个峰后破坏割线模量 来表征孔隙水压力对岩体应力应变曲线转型的影响,并结合斜坡失稳的刚度判据,进一步分析了孔隙水压力对应力应变曲线峰后斜率的影响,随孔隙水压力的增加,峰后曲线斜率变陡,峰后刚度增大,即材料的均匀性、脆性增大与刚度比k减小。因此在围压等外界环境因素及系统内部条件不变时,刚度比 时就存在一个临界孔隙水压力 ,也就是说降雨等涨落因素引起滑面介质中的孔隙水压力大于这个临界值时斜坡就易于发生突变失稳。从而进一步加深认识了降雨等外部涨落因素对斜坡系统失稳的重要触发作用。