ABSTRACT
This paper deals with the numerical implementation of a hyperelastoplastic model for sands based on critical state theory (CST). Additional features introduced in the model capture sand behavior more realistically but also introduce new challenges in the numerical implementation of the model. We demonstrate that a fully implicit algorithm is achievable by showing how consistent linearization of the residual vectors results in a solution algorithm exhibiting asymptotic quadratic rate of convergence. We perform two representative numerical examples, at a stress point, to demonstrate the robustness of the proposed implicit scheme.