请各位大神不吝赐教，模拟中碰到了内聚力单元应力不退化问题。主要是多向载荷下失效问题。

在单向载荷下，内聚力本构我是用起来没啥问题，和网上教程一样。如下图所示。

但是在多向载荷下，碰到了达到损伤但是应力不退化的情况，如下图所示。我定义了二次应力准则作为损伤判据，BK能量准控制退化速度（BK能量准则大小应该只影响退化快慢）。S33,S13,S23单独最大损伤为130MPa，195MPa，195MPa。到了第四个点，某单元S33,S13,S23分别为--289MPa，-194MPa，-16MPa。由于在法向是压力，故而满足:

（0/130）^2+(-194/195)^2+(-16/195)^2>1 ，QUADSCRT=1发生失效。

但是在失效之后，我看到这个单元应力应变曲线我懵了。

可以看出，失效以后（黑色方块第四个点），虽然S33和S23没有变大，但是S13依然在增大。这导致我的模拟不正确，这该怎么处理呢