The proposed Kirchhoff–Love shell stress resultant plasticity model extends a previously reported model for plates by complementing the constitutive law of elastoplasticity with membrane effects. This enhanced model is designed for bending dominant settings with small to moderate membrane forces. It is thus implemented in a purpose-built nonlinear mixed Eulerian–Lagrangian finite element scheme for the simulation of sheet metal roll forming. Numerical experiments by imposing artificial strain histories on a through-the thickness element are conducted to test the model against previously reported stress resultant plasticity models and to validate it against the traditional continuum plasticity approach that features an integration of relations of elastoplasticity in a set of grid points distributed over the thickness. Results of actual roll forming simulations demonstrate the practicality in comparison to the computationally more expensive continuum plasticity approach.