当有单元满足VUMAT中的damage criteria, 出现singular matrix问题

基于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? 有了解的麻烦指导小弟一下。