This paper investigates the implementation of the Extended-Matsuoka–Nakai yield criterion on a strict Limit Analysis finite element formulation. The current approach is based on a three-field mixed finite element model and the Alternating Direction Method of Multipliers optimization algorithm. With the support of duality principles two variants are derived, the lower bound and the upper bound element. The main novelty of this work is the development of an efficient iterative predictor–corrector scheme, customized for the Extended-Matsuoka–Nakai. This scheme is an indispensable requirement for this formulation. To conclude four numerical examples are presented to assess the effectiveness and efficiency of the numerical tool.
The paper describes the three-dimensional numerical implementation of the Lade-Duncan criterion in a finite element limit analysis (FELA) code. Validation is done using examples with a known solution. To conclude the proposed numerical tool is applied to the calculation of the ultimate bearing capacity of square footing.