在冲击过程的模拟中，出现了在冲击后半程中能量激增的问题，包括动能和总能量的异常增长，以及弹性应变能减为负值，使得计算出现错误。从结果上看，在冲击后半程会出现冲头将板带起的奇怪现象，并且冲头回弹的速度会大于碰撞之前的速度。下图为接触力曲线

下图为速度曲线

下图为动能曲线

下图为弹性应变能曲线

下图可看出动能和总能量的异常增长

我已经尝试过在不调用Vumat子程序时计算冲击过程，得到的结果是比较合理的，不会出现上述现象，所以模型以及接触设置应该不会有问题，此外，我还陆续改变了一些操作，包括：改变单元类型（C3D8R换成C3D8）、细化网格尺寸、缩小时间增量、将冲头从刚体换为可变形体、调整接触设置（由面面接触改为通用接触）等等，都没起到什么效果。

但是如果将网格划分的稀疏一些，结果会变得比较合理：动能不再异常增长，也不会出现回弹时板被带起现象。

请大家指导下可能是哪方面存在问题导致这种结果。