“…More recently, several methods of modern numerical optimization have been successfully used to solve elasto-plastic problems: sequential quadratic programming method (Wieners, 2007), complementarity problem solvers (Tangaramvong et al, 2012;Zheng et al, 2020), accelerated gradient algorithms (Kanno, 2016;Shimizu and Kanno, 2020). Limit and shakedown analyses are problems for which an optimization approach of plasticity is natural, using second order cone programming technique and quadratic optimisation (Bisbos et al, 2005;Delbecq et al, 1977;Makrodimopoulos and Martin, 2005a,b;Mercier, 1976;Pastor, Thoré, et al, 2008;Pastor, Pastor, et al, 2015). In (Krabbenhøft, Lyamin, Sloan, and Wriggers, 2007;Krabbenhøft, Lyamin, Hjiaj, et al, 2005), developments are carried out for Mohr-Coulomb yield criterion in the context of cone variational inequality, extending the pioneering works of Berga and De Saxcé (1994) and Hjiaj et al (2003).…”