1 //- Increment of Green finite strain tensor
2 DEpsilon = symm(gradDU) + 0.5*symm((gradDU & gradU.T()) + (gradU & gradDU.T()) + (gradDU & gradDU.T()));
4 //- Increment of second Piola-Kirchhoff stress tensor
5 DSigma = 2*mu*DEpsilon + lambda*(I*tr(DEpsilon));