EEDHD: Compute the positive and negative part of the elastic energy density associated with a strain tensor EEDMinus = SmoothPositiveSquare(-trace(e)) AI.I/2/N^2 EEDPlus = Ae.e/2 - EEDMinus with N = 3 for 3D and N = 2 for 2D, where e is the strain tensor, A is the Hookes law tensor
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(MEF90DefMechSplitHD), | intent(in) | :: | self | |||
| class(MEF90HookesLaw), | intent(in) | :: | HookesLaw | |||
| class(mef90Mat), | intent(in) | :: | phi | |||
| real(kind=C_DOUBLE), | intent(out) | :: | EEDPlus | |||
| real(kind=C_DOUBLE), | intent(out) | :: | EEDMinus | |||
| integer(kind=C_INT32_T), | intent(inout) | :: | ierr |