Implicit integration scheme applied to a multi-surface constitutive model