For accurate prediction of mechanical behaviour of cemented clay, a realistic description of its constitutive relation is important. In this study, within the framework of critical state soil mechanics, the bounding surface concept is employed to model the behaviour of cement stabilized soil. The plastic modulus proposed is load path-dependent where radial mapping rule is used for the imaginary points on the bounding surface. This bounding surface model, which includes tensile strength due to cementation effect, is extended to general stress space for implementing it in ABAQUS. The parameters for the model can be easily obtained from conventional triaxial drained and undrained tests and 1D consolidation test. With different boundary conditions, simulations are carried out for drained and undrained triaxial tests for a wide variety clay with different cement contents. In simulation, key attention is given to stress–strain behaviour and ultimate strength of cement stabilized clay. © 2020, Springer Nature Singapore Pte Ltd.