The paper addresses the numerical simulation of strain localisation in stiff clays that exhibit softening behaviour. An elastoplastic constitutive model developed to incorporate key features of stiff clay behaviour is described first. A non-local formulation is then introduced for the regularisation of the analysis of localisation. A series of analyses were conducted to explore relevant aspects of the numerical simulation of localisation. A 3D analysis was also performed to assess the suitability of the approach presented for 3D applications. Finally, application to the simulation of a laboratory test on Beaucaire marl results in an excellent reproduction of experimental observations.