Linear smoothed polygonal and polyhedral finite elements

Amrita Francis, Alejandro Ortiz-Bernardin, Stéphane P A Bordas, Sundararajan Natarajan

Research output: Contribution to journalArticle

40 Citations (Scopus)

Abstract

The strain smoothing technique over higher order elements and arbitrary polytopes yields less accurate solutions than other techniques such as the conventional polygonal finite element method. In this work, we propose a linear strain smoothing scheme that improves the accuracy of linear and quadratic approximations over convex polytopes. The main idea is to subdivide the polytope into simplicial subcells and use a linear smoothing function in each subcell to compute the strain. This new strain is then used in the computation of the stiffness matrix. The convergence properties and accuracy of the proposed scheme are discussed by solving a few benchmark problems. Numerical results show that the proposed linear strain smoothing scheme makes the approximation based on polytopes able to deliver the same optimal convergence rate as traditional quadrilateral and hexahedral approximations. The accuracy is also improved, and all the methods tested pass the patch test to machine precision.

Original languageEnglish
Pages (from-to)1263-1288
Number of pages26
JournalInternational Journal for Numerical Methods in Engineering
Volume109
Issue number9
DOIs
Publication statusPublished - 2 Mar 2017

    Fingerprint

Cite this

Francis, A., Ortiz-Bernardin, A., Bordas, S. P. A., & Natarajan, S. (2017). Linear smoothed polygonal and polyhedral finite elements. International Journal for Numerical Methods in Engineering, 109(9), 1263-1288. https://doi.org/10.1002/nme.5324