abaqus vumat子程序损伤演化阶段（应力更新）问题

我想编写一个关于Johnson-Cook损伤的vumat子程序，然后做了一个单轴拉伸的案例，但是进入损伤阶段后力和位移曲线震荡严重

我想知道损伤演化阶段应力如何更新，只需要改变弹塑性阶段的刚度矩阵吗