An unsymmetric 4‐node, 8‐DOF plane membrane element perfectly breaking through MacNeal's theorem

Publisher: John Wiley & Sons Inc

E-ISSN: 1097-0207|103|7|469-500

ISSN: 0029-5981

Source: INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Vol.103, Iss.7, 2015-08, pp. : 469-500

Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.

Previous Menu Next

Abstract

SummaryAmong numerous finite element techniques, few models can perfectly (without any numerical problems) break through MacNeal's theorem: any 4‐node, 8‐DOF membrane element will either lock in in‐plane bending or fail to pass a C0 patch test when the element's shape is an isosceles trapezoid. In this paper, a 4‐node plane quadrilateral membrane element is developed following the unsymmetric formulation concept, which means two different sets of interpolation functions for displacement fields are simultaneously used. The first set employs the shape functions of the traditional 4‐node bilinear isoparametric element, while the second set adopts a novel composite coordinate interpolation scheme with analytical trail function method, in which the Cartesian coordinates (x,y) and the second form of quadrilateral area coordinates (QACM‐II) (S,T) are applied together. The resulting element US‐ATFQ4 exhibits amazing performance in rigorous numerical tests. It is insensitive to various serious mesh distortions, free of trapezoidal locking, and can satisfy both the classical first‐order patch test and the second‐order patch test for pure bending. Furthermore, because of usage of the second form of quadrilateral area coordinates (QACM‐II), the new element provides the invariance for the coordinate rotation. It seems that the behaviors of the present model are beyond the well‐known contradiction defined by MacNeal's theorem. Copyright © 2015 John Wiley & Sons, Ltd.