New findings on limit state analysis with unstructured triangular finite elements

• Published in: International journal for numerical methods in engineering 2013-04, Vol.94 (4), p.400-417
• Main Authors: Preisig, M., Prévost, J. H.
• Format: Article
• Language: English
• Subjects: Computation; Displacement; Exact sciences and technology; Finite element method; finite element methods; Fracture mechanics (crack, fatigue, damage...); Fundamental areas of phenomenology (including applications); Inelasticity (thermoplasticity, viscoplasticity...); Limit load; Locking; low-order elements; Mathematical analysis; Mathematical models; Mathematics; mesh generation; Meshing; Methods of scientific computing (including symbolic computation, algebraic computation); Numerical analysis. Scientific computation; Physics; plasticity; Sciences and techniques of general use; Solid mechanics; stability; Static elasticity (thermoelasticity...); Structural and continuum mechanics
• ISSN: 0029-5981; 1097-0207
• DOI: 10.1002/nme.4455

SUMMARY The emergence of high performance unstructured mesh generators enables the use of low‐order triangles and tetrahedras for a wide variety of problems, where no structured quadrangular mesh can be easily constructed. Some of the results acknowledged by the community of computational mechanics therefore need to be carefully looked at under the light of these recent advances in unstructured and adaptive meshing. In particular, we question evidence of locking during plastic flow, as it has been presented by several authors in the past. We claim that locking can also manifest through stress oscillations. Overshoot of limit loads or displacement locking can be explained by other causes: Insufficient mesh refinement or comparison with analytical solutions that are based on assumptions that are inconsistent with the finite element model. We present a stabilized Galerkin mixed method that allows accurate and converging limit loads to be computed for any elastic‐plastic problem using standard equal‐order interpolation for both displacements and mean stress for low‐order linear triangular elements. The method is attractive for simulations where the cost of meshing is more important than the added computational cost of an additional nodal unknown required by the mixed formulation. Copyright © 2013 John Wiley & Sons, Ltd.