复合材料冲击用了层板以及COHESIVE子程序能量不对

2901110298

各位，贤弟在做冲击的时候用了层板以及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加了之后就会变得更大。
小弟真的不知道怎么弄了，请各位大神请教。前面应力更新以及破坏的部分没有问题，就是能量这一块不对劲。

枯藤老树昏鸦

学习一下