A novel Galerkin-like weakform and a superconvergent alpha finite element method (S αFEM) for mechanics problems using triangular meshes

A carefully designed procedure is presented to modify the piecewise constant strain field of linear triangular FEM models, and to reconstruct a strain field with an adjustable parameter α. A novel Galerkin-like weakform derived from the Hellinger–Reissner variational principle is proposed for establ...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics Vol. 228; no. 11; pp. 4055 - 4087
Main Authors: Liu, G.R., Nguyen-Xuan, H., Nguyen-Thoi, T., Xu, X.
Format: Journal Article
Language:English
Published: Kidlington Elsevier Inc 20.06.2009
Elsevier
Subjects:
ISSN:0021-9991, 1090-2716
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:A carefully designed procedure is presented to modify the piecewise constant strain field of linear triangular FEM models, and to reconstruct a strain field with an adjustable parameter α. A novel Galerkin-like weakform derived from the Hellinger–Reissner variational principle is proposed for establishing the discretized system equations. The new weak form is very simple, possesses the same good properties of the standard Galerkin weakform, and works particularly well for strain construction methods. A superconvergent alpha finite element method (S αFEM) is then formulated by using the constructed strain field and the Galerkin-like weakform for solid mechanics problems. The implementation of the S αFEM is straightforward and no additional parameters are used. We prove theoretically and show numerically that the S αFEM always achieves more accurate and higher convergence rate than the standard FEM of triangular elements (T3) and even more accurate than the four-node quadrilateral elements (Q4) when the same sets of nodes are used. The S αFEM can always produce both lower and upper bounds to the exact solution in the energy norm for all elasticity problems by properly choosing an α. In addition, a preferable- α approach has also been devised to produce very accurate solutions for both displacement and energy norms and a superconvergent rate in the energy error norm. Furthermore, a model-based selective scheme is proposed to formulate a combined S αFEM/NS-FEM model that handily overcomes the volumetric locking problems. Intensive numerical studies including singularity problems have been conducted to confirm the theory and properties of the S αFEM.
Bibliography:ObjectType-Article-2
SourceType-Scholarly Journals-1
ObjectType-Feature-1
content type line 23
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2009.02.017