有关超高速碰撞的能量问题

做了一个超高速碰撞的仿真模拟，采用自适应FEM-SPH算法，有两个疑问，（1）为什么总能量会增大？（有时候得出的结果是减小）（2）滑移能为什么会是负值？利用hypermesh软件和.messag文件都没有找到初始穿透，这是什么原因呢？希望各位老师和专家能解答一下，万分感谢！