Upper bound finite element limit analysis method with discontinuous quadratic displacement fields and remeshing in non-homogeneous clays

Abstract

As is known, limit analysis is a common way to study stability and failure mechanism of geotechnical engineering. This paper presents an upper bound finite element limit analysis (UB-FELA) formulation for non-homogeneous clays using discontinuous quadratic displacement fields and second-order cone programming (SOCP). The details are provided about the application of SOCP to the UB-FELA formulation for non-homogeneous clays. To reduce the calculation error, the analytic expressions are proposed concerning the rate of external work done by body forces and power dissipation caused by the plastic deformation following discontinuous quadratic displacement fields. With SOCP solver MOSEK, large-scale optimization problems of UB-FELA can be achieved in minutes via a desktop computer. The method is then applied to analyzing the stability of an undrained square tunnel and an undrained plane strain heading. The numerical results are compared with those reported in related studies. This method helps to obtain highly precise upper bound solution and clear failure mechanism of the computational domain.