An implicit integration algorithm has been refined to predict the stress–strain–strength response of unsaturated soil under suction-controlled, multiaxial stress paths that are not achievable in a conventional cylindrical cell. The algorithm supports numerical analyses in a deviatoric plane by using a mixed control constitutive driver, in conjunction with a generalized Cam-Clay model that also incorporates the influence of a third stress invariant, or Lode-angle θ , within a constant-suction scheme. True triaxial data from a previously accomplished series of suction-controlled triaxial compression, triaxial extension, and simple shear tests on 10-cm cubical specimens of silty sand, were used for the tuning and validation of the refined algorithm. The elliptical Willam–Warnke surface was adopted for simulation of unsaturated soil response in three-invariant stress space. Reasonably satisfactory agreement was observed between experimental and predicted deviatoric stress versus principal strain response for ...