Armando Nunes Antão

The Hoek–Brown failure criterion has been widely applied to predict the strength of rock masses, demonstrating its relevance in diverse geotechnical contexts. This paper presents a novel 3D numerical implementation of the Hoek–Brown criterion in a finite-element limit analysis code. The proposed formulation is unique in its ability to produce strict upper and lower bounds for 3D problems, providing more accurate and reliable predictions of failure mechanisms compared to previous formulations. The validity of the formulation is demonstrated through comparisons with known analytical solutions or other authors’ numerical solutions. Furthermore, the proposed numerical tool is used to determine the stability of shallow circular tunnels in rock masses, highlighting its practical applicability in engineering design.

In the current paper, we propose a new three-dimensional strength criterion for rocks expressed in terms of the first principal stress invariant, I1, and the second and third invariants of the deviatoric stress tensor, J2 and J3. The design of this constitutive model conjugates the characteristics of two of the most well-known and widely used criteria in geotechnical engineering: Hoek-Brown and Matsuoka-Nakai. Its material parameters can be calibrated based on conventional axisymmetric compression and extension tests. Experimental polyaxial test data from a dozen different rock types were used to validate the current criterion.

This paper aims at contributing towards a better understanding of the non-uniform elastoplastic torsion mechanism of I-section beams. The particular case of cantilevers subjected to an end torque is analysed, which constitutes a simple yet interesting problem, since the maximum torque is very close to the so-called Merchant upper bound (MUB), with added independent maximum bishear and Saint-Venant torques. Consequently, it turns out that the maximum torque can be significantly higher than that for uniform plastic torsion. Besides the MUB, several solutions are presented and compared, namely (i) a stress resultant-based solution stemming from the warping beam theory differential equilibrium equation and (ii) solutions obtained with several beam finite elements that allow for a coarse/refined description of warping. It is found that all models are in very close agreement in terms of maximum torque (including the MUB) and stress resultants. However, the beam finite elements that allow for bishear, even with a simplified warping function, are further capable of reproducing quite accurately the stress field, as a comparison with a 3D solid finite element solution shows. Although the paper is primarily concerned with the small displacement case, the influence of considering finite rotations is also addressed.

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.

A numerical implementation of the upper-bound theorem of limit analysis is applied to determine two-dimensional (2D) and three-dimensional (3D) active horizontal earth pressure coefficients considering seismic actions through a horizontal seismic coefficient. Results are obtained for vertical wall, horizontal soil, different friction angles of the soil, soil-to-wall friction ratios, horizontal seismic coefficients and wall width-to-height ratios. The few cases for which 3D active earth pressure coefficients are available in the literature using upper-bound methods were used for comparison with the corresponding earth pressure coefficients obtained in this study. This showed a general improvement of these results, which allows expecting a good accuracy for the set of cases studied. The ratios between the 3D and 2D horizontal active earth pressure coefficients are found to be practically independent of the soil-to-wall friction ratio. An equation is proposed for calculating these ratios. This equation can be easily used in the design of geotechnical structures requiring the determination of 3D active earth pressure coefficients.

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.

This paper deals with the determination of upper and lower bounds for the three-dimensional undrained stability of shallow tunnels. The tunnel is circular and a distance between its face and its lining is considered. The soil shear strength is modeled using the Tresca criterion. Results of the upper and lower bounds of the stability number are presented, for several geometrical and resistance configurations and their comparison with previous results is made, showing the clear improvement obtained. Finally, equations approaching the stability number are proposed.

Computational limit analysis methods invariably lead to the need to solve a mathematical programming problem. The alternating direction method of multipliers (ADMM) is one versatile and robust technique to solve non-linear convex optimization problems that has recently found applications in a wide range of fields. Its solution scheme, based on an operator splitting algorithm, is not only easy to implement but also suitable to efficiently solve large-scale variational problems. Starting from the ADMM framework, we derive a strict upper bound finite element formulation using a two-(primal)-field approximation, one for the velocity field and the other for the plastic strain rate field. Next, following a similar approach, we develop a novel strict lower bound formulation. Here, the two-(primal)-field model is based on a redundant approximation of the stress field. Duality principles are then explored in order to unify these two formulations.The effectiveness of this approach is demonstrated on test problems and, to conclude, some considerations are made about the performance results.

This paper describes a novel upper-bound formulation of limit analysis. This formulation is an innovative variant of an existing two-field mixed formulation based on the augmented Lagrangian method also developed by the authors. A natural approach is used to describe the deformation of each finite element. Furthermore, and in contrast to the previous formulation, two independent field approximations are now both used to define the velocity field, defined globally and at element level. It is shown that this feature allows a governing system of uncoupled linear equations to be obtained. Some numerical examples in plane strain conditions are presented in order to illustrate the current model performance. In conclusion, the potential and advantages of this new approach are discussed.

A strict upper-bound limit analysis finite element formulation is used to estimate shape factors s_gamma and s_q for determining the bearing capacity of shallow foundations using the classic bearing capacity formula. The finite element formulation uses a quadratic approximation for the velocity field, an extension of a previously published Augmented Lagrangian formulation with a linear velocity field, and was implemented for a parallel processing environment. Results from determining the limit loads under three-dimensional conditions are presented and compared with previously published data. The results obtained allow a strict upper-bound determination of the shape factors. Furthermore, a practical proposal for these factors is made and compared with other proposals made by other authors.