A multigrid method for three-dimensional elasticity and algebraic convergence estimates

摘要

A two-level iterative method for 3D elasticity problems discretized by quadratic elements is proposed. The coarse level consists of linear elements. We develop a sharpened algebraic convergence theory for two-level methods with smoothing by Gauss-Seidel, which applies to more general situations as well. The convergence estimates do not deteriorate for strongly discontinuous coefficients. Observed convergence rates are compared with theoretical estimates. The method is competitive with a standard direct solver.