Paper detail

Convergence in the incompressible limit of new discontinuous Galerkin methods with general quadrilateral and hexahedral elements

Standard low-order finite elements, which perform well for problems involving compressible elastic materials, are known to under-perform when nearly incompressible materials are involved, commonly exhibiting the locking phenomenon. Interior penalty (IP) discontinuous Galerkin methods have been shown to circumvent locking when simplicial elements are used. The same IP methods, however, result in locking on meshes of quadrilaterals. The authors have shown in earlier work that under-integration of specified terms in the IP formulation eliminates the locking problem for rectangular elements. Here it is demonstrated through an extensive numerical investigation that the effect of using under-integration carries over successfully to meshes of more general quadrilateral elements, as would likely be used in practical applications, and results in accurate displacement approximations. Uniform convergence with respect to the compressibility parameter is shown numerically. Additionally, a stress approximation obtained here by postprocessing shows good convergence in the incompressible limit.