Numerical implementation of the Hoek-Brown material model with strain hardening

Each location requires its own return algorithm.The choice of return algorithm is guess upon based on the concept of boundary planes mentioned in Clausen & Damkilde (2008). Each return algorithm seeks to solve the relevant equations based on a Newton-Raphson iteration scheme.