ABAQUS使用VUMAT子程序计算时出现的问题，为何单元会被无缘无故删除且有一些灰色单元出现？ 50

我现在做的是复合材料相关的低速冲击问题，使用VUMAT子程序，使用的本构关系是张超论文中的粘弹性本构模型，接触用的是面面接触，主面是冲头外表面，从面是材料的节点集合，但是在计算中出现冲头刚与材料接触就出错，且单元会被无缘无故删除且有一些灰色单元出现，我的VUMAT子程序如果使用简单的线弹性本构是可以进行计算的，我贴上我的图，还请相关的大佬们帮我看看这是什么情况，是什么导致这种情况的发生，在这里先谢谢各位了！如果还有什么需要的我在上传