The algorithm proposed in [9] for incremental elastoplasticity is extended and applied to shakedown analysis. Using the three field mixed finite element proposed in [22] a series of mathematical programming problems or steps, obtained from the application of the proximal point algorithm to the static shakedown theorem, are obtained. Each step is solved by an Equality Constrained Sequential Quadratic Programming (EC-SQP) tech- nique that allows a consistent linearization of the equations improving the computational efficiency
