提出了基于Euler向后积分的自适应子增量本构积分算法,推导了相应的一致切线模量矩阵;通过引入伪屈服函数(塑性势函数),提出了对屈服面角点应力区进行两个方向应力投射的本构积分算法,使超出屈服面的试应力收敛到角点;推导了两个投射方向的一致切线模量矩阵;采用赋小值方法解决0应力屈服的问题。用上述方法编制了基于D—P准则的理想弹塑性模型ABAQUS—UMAT子程序,并进行了算例验证。
Based on Euler backward integral method, an integration algorithm of constitutive equation is suggested for adaptive sub-stepping schemes. And the consistent tangent matrices are derived. According to the projection vector of trial stress beyond the yielding surface, the pseudo yielding function and plastic potential function is introduced for apex area of yielding surface. So, integration algorithms adopted of two projection vectors make trial stress in apex area to return to apex point of yielding surface. The consistent tangent matrices for two vectors are also derived. Zero stress point in the yielding surface is replaced by minimum value. An ABAQUS-UMAT subroutine using D-P criterion for linear elastic-perfectly plastic model is programmed; and two examples are analyzed to validate these methods.