Abstract. The paper improved S-FEMs formulations with an enriched displacement field, making use of modified Allman's shape functions. This mixed interpolation is the natural context in performing lower bound strategy for shakedown, limit analysis and elastoplastic analysis. The model takes advantages from the simplicity and few addressed requirements for good performances in nonlinear analysis. The simple assumption made for the stress field regards the convenience of using self-equilibrated stress interpolations in Cartesian coordinates. In the proposed composite elements the stress is discontinuous on the element and across their sides and the mesh of the elements is coincident with the discretization of the geometry. This stress interpolation is able to address the discontinuities in the plastic strain and, in such a way, to define in their description a finer mesh with respect the basic grid.