ABSTRACT
An implicit integration algorithm has been developed to simulate stress-strain response of unsaturated soil under controlled-suction multiaxial stress states. The algorithm supports numerical analyses in the π-plane using a mixed-control constitutive driver, along with a generalized cam-clay model, incorporating a third stress invariant, Lode angle θ, within a constant-suction scheme. True triaxial data (σ₁ ≥ σ₂ ≥ σ₃ ) from constant- suction triaxial extension (TE), triaxial compression (TC), and simple shear (SS) tests, conducted on cubical specimens of silty sand, were used for tuning of the developed algorithm. The Willam-Warnke surface was adapted to the unsaturated case for simulation of soil response in three-invariant stress space (p: q: θ). Agreement between the observed and predicted soil responses highlight the potential of the developed algorithm for analyses of boundary-value problems involving unsaturated soil deposits.