A new method for implicit integration of the Mohr-Coulomb non-smooth multisurface plasticity models is presented, and Koiter’s requirements are incorporated exactly within proposed algorithm. Algorithmic numerical complexities identified introduced by nonsmooth intersections surfaces; then, a projection contraction algorithm applied to solve classical Kuhn–Tucker complementary equations which p...