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

  • Authors:
  • G. R. Liu;H. Nguyen-Xuan;T. Nguyen-Thoi;X. Xu

  • Affiliations:
  • Singapore-MIT Alliance (SMA), E4-04-10, 4 Engineering Drive 3, Singapore 117576, Singapore and Center for Advanced Computations in Engineering Science (ACES), Department of Mechanical Engineering, ...;Singapore-MIT Alliance (SMA), E4-04-10, 4 Engineering Drive 3, Singapore 117576, Singapore and Department of Mechanics, Faculty of Mathematics and Computer Science, University of Science HCM, Viet ...;Center for Advanced Computations in Engineering Science (ACES), Department of Mechanical Engineering, National University of Singapore, 9 Engineering Drive 1, Singapore 117576, Singapore and Depar ...;Center for Advanced Computations in Engineering Science (ACES), Department of Mechanical Engineering, National University of Singapore, 9 Engineering Drive 1, Singapore 117576, Singapore

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2009

Quantified Score

Hi-index 31.45

Visualization

Abstract

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. 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@aFEM) is then formulated by using the constructed strain field and the Galerkin-like weakform for solid mechanics problems. The implementation of the S@aFEM is straightforward and no additional parameters are used. We prove theoretically and show numerically that the S@aFEM 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@aFEM can always produce both lower and upper bounds to the exact solution in the energy norm for all elasticity problems by properly choosing an @a. In addition, a preferable-@a 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@aFEM/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@aFEM.