摘要：Hoek-Brown（H．B）屈服准则存在奇异性，在岩石（体）应力状态的低（拉）围压区存在低估岩石抗拉强度的现象。为了在数值计算中处理H．B屈服准则的奇异性，解决H-B屈服准则在低（拉）围压条件下低估岩石抗拉强度的问题，采用主应力回映算法和不同应力状态选取不同的mi常数的方法，建立基于H．B屈服准则和考虑岩石屈服前的非线性特性的本构模型；在Ls-dyna的自定义材料接口程序中利用Fortran对模型进行二次开发，通过删除单元的方法模拟岩石试件的破裂过程，实现带破坏模式的低围压条件下的模型精确数值求解，获得岩石单轴压缩与巴西劈裂的三维破坏形态；岩石单轴压缩实验、巴西劈裂实验和声发射实验的三种力学特性实验结果与模型数值计算结果比较表明，实验结果与数值计算结果吻合较好，模型及其算法有较高的准确性。

关键词：岩石力学 非线性弹塑性 岩石本构模型 数值实现 

单位：重庆大学煤矿灾害动力学与控制国家重点实验室 重庆400044 重庆大学复杂煤气层瓦斯抽采国家地方联合工程实验室 重庆400044 中国矿业大学矿业工程学院 江苏徐州221116