A multi-stage return algorithm for solving the classical damage component of constitutive models for rocks, ceramics, and other rock-like media