A new method for the incremental analysis of elastoplastic associated materials is presented. The method fully retains all the equations and variables of the problems at the same level and uses a sequential quadratic programming with equality constraints to solve in an efficient and robust fashion the elastoplastic step equations derived by means of a suitable mathematical programming formulation of the problem. The new proposal is compared with standard strain driven formulations which use a return mapping by closest point projection schemes. The numerical tests performed show a good performance and a great robustness of the proposed formulation also in the case of multi–surface elastoplasticity