Formulation and Implicit Numerical Integration of a Kinematic Hardening Model for Unsaturated Soils

Currently, our understanding of material‐scale deterioration resulting from meteorologically induced variations in pore water pressure and its significant impact on infrastructure slopes is limited. To bridge this knowledge gap, we have developed an extended kinematic hardening constitutive model for unsaturated soils that refines our understanding of weather‐driven deterioration mechanisms in heterogeneous clay soils. This model has the capability of predicting the irrecoverable degradation of strength and stiffness that has been shown to occur when soils undergo wetting and drying cycles. The model is equipped with a fully coupled and hysteretic water retention curve and a hysteretic loading–collapse curve and has the capability to predict the irrecoverable degradation of strength and stiffness that occurs during cyclic loading of soils. Here, we employ a fully implicit stress integration technique and give particular emphasis to deriving a consistent tangent operator, which includes the linearisation of the retention curve. The proposed algorithm is evaluated for efficiency and performance by simulating various stress and strain‐driven triaxial paths, and the accuracy of the integration technique is evaluated through the use of convergence curves.