VUMAT的收敛问题

yuanqiuping

本人编了一个chaboche本构的VUMAT（本人在计算ΔP的时候用的是牛顿迭代），对一个单元进行拉压循环加载计算，第一个半循环进入苏醒没有问题，但是后面的循环中只要一进入塑性，应力马上就跳回到弹性阶段然后经过一段弹性计算后才能顺利的进入塑性。我不知道是哪里出了问题，请大侠指点，附件有计算得到的应力应变曲线