Consistent tangent stiffness matrix

The consistent (algorithmic) stiffness matrix

$$J=\frac{\partial\Delta \boldsymbol{\sigma}}{\partial\Delta \boldsymbol{\epsilon}}$$

is in the umat calculated numerically by means of a direct perturbation of the strain increment $\Delta \boldsymbol{\epsilon}$.