Finite element implementation of an isotach elastoplastic constitutive model for soft soils

An isotach elastoplastic constitutive model devised by Yang et al. (2016) and referred to as the Hunter Clay (HC) model attempts to capture a number of key behaviours of soft soils within a critical state framework, namely destructuration, fabric anisotropy and rate dependency, the latter often manifesting in creep settlement. Finite element implementation of the HC model is a useful means by which its application to practical problems can be facilitated. However, there are a number of significant challenges associated with the translation of isotach elastoplastic models into a finite element setting. In this paper, a detailed discussion of these challenges is undertaken and a new finite element implementation of the HC model is subsequently developed. This includes sophisticated numerical integration algorithms which employ automatic time substepping for solution of the governing finite element equations. The ability of the implemented HC model to predict the mechanical behaviour of soft soils under 1D compression is investigated via simulation of laboratory tests carried out on Ballina clay by Pineda et al. (2016) and Parkinson (2018).