基于ABAQUS的弹脆塑性本构模型二次开发及其应用

在复杂应力状态下,大理岩、花岗岩等脆性岩石呈弹脆塑性破坏形式。利用ABAQUS提供的用户材料子程序UMAT的二次开发平台,编制完成了基于D-P准则的弹脆塑性模型的UMAT接口程序。论文简要介绍了ABAQUS有限元软件及其UMAT用户材料子程序,重点对基于D-P屈服准则的弹脆塑性本构模型的应力跌落公式进行了推导,并讨论了有关的应力更新算法。最后应用该程序对某引水隧洞工程进行了数值分析,并与理想弹塑性情况下的计算结果进行对比研究。结果表明:与理想弹塑性模型相比,弹脆塑性模型得到的洞周围岩变形和塑性区范围较大,与实际情况更为接近。

小湾水电站地下厂房洞室群弹脆塑性分析

针对小湾水电站枢纽区岩性特点与所处的地应力环境,采用脆塑性本构模型,分析了地下厂房洞室群开挖后围岩的变形与破损区;同时比较了使用脆塑性与弹塑性本构模型时计算结果的差异。由于地下厂房洞室群区域受到多条断层切割,洞室开挖完成后,在断层经过的部位,断层对洞室的变形与破损区有一定的控制作用。针对工程中所关心的断层对洞室稳定性的影响程度,通过断层材料的置换,定量地给出了在断层经过的区域、断层对洞室表面的变形的影响以及断层对洞室围岩的破损区的影响。

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

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

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

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