Internal energy is computed in LS-DYNA based on the six components of stress and strain (tensorial values).
The calculation is done incrementally for each element as follows:
(IE)new = (IE)old + sum over all six directions of (stress * incremental strain * volume)
The internal energies of all the elements is summed to give the total internal energy.
In addition, energy dissipated due to stiffness damping (*DAMPING_PART_STIFFNESS) is added
to the internal energy if RYLEN=2 in *CONTROL_ENERGY.
The calculation of energies is done for user-defined materials outside the material subroutine.
To combat this spurious effect,