Reduced modified quadratures for quadratic membrane finite elements

Reduced integration is frequently used in evaluating the element stiffness matrix of quadratically interpolated finite elements. Typical examples are the serendipity (Q8) and Lagrangian (Q9) membrane finite elements, for which a reduced 2 × 2 Gauss–Legendre integration rule is frequently used, as opposed to full 3 × 3 Gauss–Legendre integration. This ‘softens’ these element, thereby increasing accuracy, albeit at the introduction of spurious zero energy modes on the element level. This is in general not considered problematic for the ‘hourglass’ mode common to Q8 and Q9 elements, since this spurious mode is non-communicable. The remaining two zero energy modes occurring in the Q9 element are indeed communicable. However, in topology optimization for instance, conditions may arise where the non-communicable spurious mode associated with the elements becomes activated. To effectively suppress these modes altogether in elements employing quadratic interpolation fields, two modified quadratures are employed herein. For the Q8 and Q9 membrane elements, the respective rules are a five and an eight point rule. As compared to fully integrated elements, the new rules enhance element accuracy due to the introduction of soft, higher-order deformation modes. A number of standard test problems reveal that element accuracy remains comparable to that of the under-integrated counterparts.