Novel numerical strategy for solving strongly coupled elastoplastic damage models with explicit return algorithms: Application to geomaterials