Numerical integration of elastoplastic constitutive equations for geomaterials with extended hardening rules