隧道围岩弹-脆-塑性应力的三剪统一解

基于三剪统一强度准则,采用更符合脆性岩石峰后强度特性的弹-脆-塑性模型,推导了考虑第二主应力和脆性软化共同影响的隧道围岩应力解,并与广义Hoek-Brown强度准则解进行了比较。研究结果表明:第二主应力、围岩材料模型、脆性软化均对围岩应力大小及分布具有显著影响。本文准则参数b=0时的应力解析解与广义Hoek-Brown强度准则应力解的峰值和分布规律均有较好的一致性,说明了本文应力解的合理性;当岩体强度具有明显的第二主应力效应时,广义Hoek-Brown强度准则已不再适用,但本文准则参数b>0时的应力解则具有很好的适用性。同时,理想弹-塑性模型低估了围岩塑性区范围,在准则参数b分别为0、0.05、0.1的情况下,弹-脆-塑性模型的围岩塑性区半径R相比理想弹-塑性模型平均可以增大45.5%;考虑第二主应力效应可以更加充分发挥岩石材料的强度潜能,随着准则参数b的增大,即第二主应力σ2效应的增大,围岩塑性区半径R和临界支护力py均不断减小,相比准则参数b=0时最大值分别减小了18.9%和10.8%。峰后粘聚力cr对围岩应力的影响也很显著,随着峰后粘聚力cr的增加,塑性区半径R不断减小,cr=0.11 MPa时的塑性区厚度比cr=0.055 MPa时减小43.3%。

基于统一强度理论的扰损围岩巷道脆塑性解答

钻爆法施工常造成岩质巷道扰损区围岩性能劣化。针对扰损围岩、扰损外围岩的弹塑性状态下不同的扰损围岩巷道力学模型,采用统一强度理论、弹-脆-塑性模型和非关联流动法则,建立扰损围岩巷道有关应力、位移和塑性区半径的脆塑性解答,探讨巷道力学模型的转换路径及判定方法,并将结果进行对比、验证,获得各因素的影响规律。研究表明:提出的巷道脆塑性解答合理反映了中间主应力、围岩扰损范围和峰后强度、剪胀特性的综合影响,通过参数组合变换可退化为多种解答并得到正确性验证,具有重要的理论意义和应用前景;巷道力学模型判定与中间主应力、围岩扰损范围密切相关,显著影响巷道应力-位移分布和围岩特征曲线;支护压力和巷道稳定变形随着围岩峰后强度下降、剪胀特性参数增加而明显增大。

考虑轴向力和渗透力时圆形隧道广义Hoek-Brown解

为了研究渗透水压力和轴向应力共同作用时隧道围岩的应力和位移变化趋势,将圆形隧道简化为轴对称模型,假定渗透场以渗透体积力作用在原应力场,以围岩开挖断面为假定平面,引入垂直于该平面的轴向应力。基于广义Hoek-Brown强度准则和非关联流动法则,推导出考虑轴向应力和渗透场共同作用时弹-脆-塑性围岩的应力和位移非线性解,采用数值算例分析了轴向应力和渗透力共同作用时隧道围岩塑性区应力场和位移场的分布规律。计算结果表明:与无渗透水压力作用下的模型相比,渗透力作用使得围岩塑性区各点位移增大,并且内外水头差越大,位移增大越明显。当轴向应力为中主应力时,围岩塑性区半径和塑性区各点应力增大,轴向应力为大主应力和小主应力时,围岩塑性区半径和塑性区应力变化较小。因此,渗透力的存在不利于隧道围岩的稳定性,并且轴向应力的大小对于富水地区隧道围岩的应力和位移分布具有较大影响。在施工设计时考虑渗透力以及轴向应力的共同影响对于保证隧道围岩稳定性具有重要意义。

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

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

基于环向应力岩爆判据的不同级别岩爆区分布的数值模拟研究

为了避免莫尔-库仑和虎克-布朗屈服准则各自的缺陷,采用C++语言对FLAC进行了二次开发,将它们复合在一起,提出了一种弹-脆-塑性的本构模型。采用徐林生及Russense判据(环向应力与单轴抗压强度之比)作为岩爆判据,研究了在不同围压条件下,圆形洞室开挖之后围岩再次平衡时已发生的各种级别岩爆区的分布规律,对严重、中等及轻微岩爆区的关系进行了讨论。将数值计算结果与现场观测结果之间的差异及共同点进行了定性分析。研究发现,随着围压的增加,围岩中发生各种级别岩爆的单元数不仅增多,而且向围岩深部迁移。