This paper presents a three-dimensional constitutive model for natural clay that includes creep, anisotropy and structure, as well as a theoretical means to estimate the range for anisotropy- and structure-related parameters, as needed for parameter optimisation. Creep-SCLAY1S is an extension of the Creep-SCLAY1 model proposed by Sivasithamparam et al. (Comput Geotech 69:46–57, 2015 ) which includes the effects of bonding and destructuration. The model needs 14 model parameters, of which five are similar to those used in the modified Cam–Clay model. A method is developed to quantify the range for the three parameters related to structure and anisotropy that cannot be derived directly from experimental data. The theoretically derived range compares favourably with the values found in the literature. As a result, the model now can be used with more confidence, enabling sensitivity analysis and systematic parameter derivation with optimisation techniques.