各位，小弟在做冲击的时候用了层板以及cohesive单元的vumat子程序，只使用单层板的子程序的时候，总能量变化不大，但是加入了cohesive子程序之后，刚开始能量就变大了。程序中的能量公式为stressPower = half*(stressNew(km,1)*strainNew(km,1)+     1  stressNew(km,2)*strainNew(km,2)+
     2  stressNew(km,3)*strainNew(km,3))    

       enerInternNew(km) = enerInternOld(km) + stressPower

单层板的能量stresspower还除以了一个密度Density(k),cohesive加了之后就会变得更大。
小弟真的不知道怎么弄了，请各位大神请教。前面应力更新以及破坏的部分没有问题，就是能量这一块不对劲。