基于Boyce-Parks-Argon Model，小弟我寫了一個含有damage criteria的VUMAT, 用來模擬polymer的斷裂失效。當我只用單個3D element進行debug，一切正常。當damage criteria觸發，相應控制單元刪除的狀態變量為0的時候，單元就刪除了。

但如果用含多個element的模型去模擬的時候（比如tensile bar），每次當model中有單元開始滿足damage criteria，模擬會顯示singular matrix issue并停止計算，但并不aborted。

請問各位大佬，這種情況是怎么回事。應該如何debug? 有了解的麻煩指導小弟一下。