A generalized backward Euler scheme for the integration of a mixed isotropic-kinematic hardening model for clays