Upper and lower bound limit analysis of plates using FEM and second-order cone programming

This paper presents two novel numerical procedures to determine upper and lower bounds on the actual collapse load multiplier for plates in bending. The conforming Hsieh–Clough–Tocher (HCT) and enhanced Morley (EM) elements are used to discrete the problem fields. A Morley element with enhanced moment fields is used. The constant moment fields is added a quadratic mode in which the pressure is equilibrated by corner loads only, ensuring that exact equilibrium relations associated with a uniform pressure can be obtained. Once the displacement or moment fields are approximated and the bound theorems applied, limit analysis becomes a problem of optimization. In this paper, the optimization problems are formulated in the form of a standard second-order cone programming which can be solved using highly efficient interior point solvers. The procedures are tested by applying it to several benchmark plate problems and are found good agreement between the present upper and lower bound solutions and results in the literature.