在挤压复材（含单元删除和cohesive层）时，怎么避免沙漏？ 100

大佬们，问一下我的芯棒挤压衬套，衬套挤压复材板孔的模型，怎么减少沙漏。我用的单元是C3D8R，之前换了一下C3D8I，不知道为什么伪应变能与内能出现了负值，然后稳定在40%以内，反而还增加了。我的模型有渐进损伤和cohesive层，两者有单元删除功能，这个单元删除是不是会影响这个比值？我的子程序中的刚度退化对沙漏是不是会有影响？我的子程序中没有涉及能量的编写，这个会不会也是其中的一个原因。我应该怎么操作才能降低沙漏？