各位，小弟在做沖擊的時候用了層板以及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加了之后就會變得更大。
小弟真的不知道怎么弄了，請各位大神請教。前面應力更新以及破壞的部分沒有問題，就是能量這一塊不對勁。