A robust preconditioner for higher order finite element discretizations in linear elasticity.

Based on the auxiliary space method, a preconditioner is studied in this paper for linear systems of equations arising from higher order finite element (FEM) discretizations of linear elasticity equations. The main idea, which is proposed by Xu ( Computing 1996; 56:215-235) for the scalar PDE, is to construct the preconditioner as a combination of a smoother and a coarse level solver, where the systems of equations arising from lower order FEM discretizations are used in the coarse level solver. It is theoretically shown that the condition number of the preconditioned systems is uniformly bounded with respect to both the problem size and moderate Poisson's ratio. When the Poisson's ratio is near the limit of 0.5, we have presented some numerical tests for the case of fourth-order FEM discretization in a combination with quadratic conforming FEM as a coarse space. The results are almost robust when Poisson's ratio is near the limit of 0.5. Copyright © 2011 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]