Abstract:Natural soft clay tends to be anisotropy developed during deposition, one dimensional consolidation, and any subsequent straining. Anisotropy of clays affects its shear strength, stress strain response, yield surface declination, and so on. An anisotropic bounding surface model is developed based on Wheeler’s elastoplastic model (S CLAY1). The developed model is programmed in ABAQUS software by writing UMAT subroutine, which is integrated in return mapping algorithm. The developed bounding surface model is validated with undrained triaxial shear test of Kaolin clay. The study shows that the developed model can well simulate the stress strain behaviors and the change of excess pore pressure of normally consolidated and over consolidated clays. Meanwhile, the developed bounding surface model can also efficiently depict the initial anisotropy and stress induced anisotropy of consolidated natural soft clay.