A numerical approach to the yield strength of shell structures

• Lower and upper bound finite elements for the limit analysis of shells are developed.
• A method for computing generalized strength criteria in membrane-bending interaction is proposed.
• Optimization problems are solved using second-order cone programming tools.
• The method is illustrated on different examples.

This work investigates the formulation of lower and upper bound finite elements for the yield design (or limit analysis) of shell structures. The shell geometry is first discretized into triangular planar facets so that previously developed lower bound equilibrium and upper bound kinematic plate finite elements can be coupled to membrane elements. The other main novelty of this paper relies on the formulation of generalized strength criteria for shells in membrane-bending interaction via an implicit upscaling procedure. This formulation provides a natural strategy for constructing lower and upper bound approximations of the exact shell strength criterion and are particularly well suited for a numerical implementation using second-order cone programming tools. By combining these approximate strength criteria to the previously mentioned finite elements, rigorous lower and upper bound ultimate load estimates for shell structures can be computed very efficiently. Different numerical examples illustrate the accuracy as well as the generality and versatility of the proposed approach.