ICCM Conferences, The 6th International Conference on Computational Methods (ICCM2015)

Font Size: 
Application of the quadrilateral area coordinates: A 4-node quadrilateral membrane element beyond MacNeal’s theorem
Pei-Lei Zhou, Song Cen

Last modified: 2015-06-28

Abstract


In this paper, the second form of quadrilateral area coordinate method (QACM-II) was employed to develop a novel 4-node quadrilateral membrane element insensitive to mesh distortion. First, two stress functions in terms of the QACM-II (S, T) are given to represent exact solutions of pure bending state in x and y directions. Their linear combination can reflect exact solutions for pure bending in all directions because the QACM-II is a local natural coordinate system. Thus, the basic analytical solutions of stresses, strains and displacements corresponding to the pure bending state for the isotropic and anisotropic plane problems can be obtained. Then, the unsymmetric element technique was adopted to formulate an unsymmetric 4-node, 8-DOF (2DOFs per node) quadrilateral membrane element. This approach is derived from the virtual work principle, and uses different test and trial functions. The test function employs the isoparametric coordinate interpolation to satisfy inter-element compatibility, while the trial function adopts a novel composite coordinate interpolation of the Cartesian coordinates and QACM-II to satisfy the completeness requirements in physical space. Due to the linear relationship between quadrilateral area and Cartesian coordinates, the order of trial displacement fields will not degrade when the meshes are distorted. Since the test displacement fields satisfy the compatibility requirements and the trial displacement fields contain analytical solutions of pure bending in all directions, the new element can produce exact solutions in both constant stress/strain patch test and pure bending problem in any directions, which means the well-known contradiction declared by MacNeal’s theorem: any 4-node, 8-DOF plane 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, has been perfectly broken through.

The authors would like to acknowledge the financial supports of the National Natural Science Foundation of China (Project No. 11272181), the Specialized Research Fund for the Doctoral Program of Higher Educationof China (Project No. 20120002110080), and Tsinghua University Initiative Scientific Research Program (Project No. 2014z09099).

An account with this site is required in order to view papers. Click here to create an account.