Summary: A practical two-surface plasticity model based on classical Dafalias/Popov and Krieg concepts was derived and implemented to incorporate yield anisotropy and three hardening effects for non-monotonous deformation paths: the Bauschinger effect, transient hardening and permanent softening. A simple-but-effective stress-update scheme avoiding overshooting was proposed and implemented. Constitutive parameters were fit to 5754-O aluminum alloy using uniaxial tension/compression data. Spring-back predictions using the resulting material model were compared with experiments and with single-surface material models which do not account for permanent softening. The two-surface model improved such predictions significantly as compared with single-surface models, while the differences between two-surface simulations and experiments were insignificant.