An efficient return mapping algorithm for elastoplasticity with exact closed form solution of the local constitutive problem