Armando Nunes Antão

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.

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 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.