摘要

Hoek-Brown强度准则的屈服面与塑性势面在棱角处导数无定义,具有数值奇异性,采用圆角进行光滑过渡只能满足一阶导数连续性,而二阶导数仍然无定义,使得棱边上一致切线模量矩阵无法正确计算,导致有限元总体平衡方程组Newton-Raphson隐式迭代二阶收敛性丧失。提出基于C2阶连续函数的广义Hoek-Brown准则屈服面与塑性势面棱角圆化方法,使得棱角处函数曲面二阶连续可导,棱边上一致切线模量矩阵可精确计算。基于ABAQUS数值开发平台,采用FORTRAN语言编制Hoek-Brown准则理想弹塑性UMAT用户子程序,通过数值算例验证所提方法的正确性。

全文