A finite element method based on C0-continuous assumed gradients

This article presents an alternative approach to assumed gradient methods in FEM applied to three-dimensional elasticity. Starting from nodal integration (NI), a general C0-continuous assumed interpolation of the deformation gradient is formulated. The assumed gradient is incorporated using the principle of Hu-Washizu. By dual Lagrange multiplier spaces, the functional is reduced to the displacements as the only unknowns. An integration scheme is proposed where the integration points coincide with the support points of the interpolation. Requirements for regular finite element meshes are explained. Using this interpretation of NI, instabilities (appearance of spurious modes) can be explained. The article discusses and classifies available strategies to stabilize NI such as penalty methods, SCNI, alpha-FEM. Related approaches, such as the smoothed finite element method, are presented and discussed. New stabilization techniques for NI are presented being entirely based on the choice of the assumed gradient interpolation, i.e. nodal-bubble support, edge-based support and support using tensor-product interpolations. A strategy is presented on how the interpolation functions can be derived for various element types. Interpolation functions for the first-order hexahedral element, the first-order and the second-order tetrahedral elements are given. Numerous examples illustrate the strengths and limitations of the new schemes. Copyright (C) 2010 John Wiley & Sons, Ltd.