摘要
Hoek-Brown强度准则的屈服面与塑性势面在棱角处导数无定义,具有数值奇异性,采用圆角进行光滑过渡只能满足一阶导数连续性,而二阶导数仍然无定义,使得棱边上一致切线模量矩阵无法正确计算,导致有限元总体平衡方程组Newton-Raphson隐式迭代二阶收敛性丧失。提出基于C2阶连续函数的广义Hoek-Brown准则屈服面与塑性势面棱角圆化方法,使得棱角处函数曲面二阶连续可导,棱边上一致切线模量矩阵可精确计算。基于ABAQUS数值开发平台,采用FORTRAN语言编制Hoek-Brown准则理想弹塑性UMAT用户子程序,通过数值算例验证所提方法的正确性。
Derivative is not defined at the edge corners of yield and plastic potential surfaces based on the generalized Hoek-Brown strength criterion and numerical singularity has been caused there. Smooth transitions with circular surface only satisfy the first derivative continuity; the second derivative is still undefined on the edges; so that the consistent tangent modulus can not be calculated correctly there. The implicit solution of nonlinear equilibrium equations in finite element method(FEM) with Newton-Raphson method is unable to achieve full quadratic convergence. C2 continuous rounding of yield and plastic potential surfaces in the π-plane has been introduced; so that the surfaces can be second-order continuous derivative on the edges; and the consistent tangent modulus matrix can be calculated accurately. Subroutine UMAT for the generalized Hoek-Brown linear elastic-perfectly plastic constitutive model has been coded with FORTRAN in ABAQUS software. Numerical examples are used to verify the correctness of the proposed method.
出处
《岩土力学》
EI
CAS
CSCD
北大核心
2018年第S1期477-487,共11页
Rock and Soil Mechanics
基金
国家重点基础研究发展计划(973)项目(No.2015CB057905)
国家自然科学基金资助项目(No.51779253
No.41672319)
湖北省自然科学基金资助项目(No.2017CFB725)~~