我現在做的是復合材料相關的低速沖擊問題，使用VUMAT子程序，使用的本構關系是張超論文中的粘彈性本構模型，接觸用的是面面接觸，主面是沖頭外表面，從面是材料的節(jié)點集合，但是在計算中出現沖頭剛與材料接觸就出錯，且單元會被無緣無故刪除且有一些灰色單元出現，我的VUMAT子程序如果使用簡單的線彈性本構是可以進行計算的，我貼上我的圖，還請相關的大佬們幫我看看這是什么情況，是什么導致這種情況的發(fā)生，在這里先謝謝各位了！如果還有什么需要的我在上傳