Preconditioned Conjugate Gradient Methods for the Refined FEM Discretizations of Nearly Incompressible Elasticity Problems in Three Dimensions

Abstract

Nearly incompressible problems in three dimensions are the important problems in practical engineering computation. The volume-locking phenomenon will appear when the commonly used finite elements such as linear elements are applied to the solution of these problems. There are many efficient approaches to overcome this locking phenomenon, one of which is the higher-order conforming finite element method. However, we often use the lower-order nonconforming elements as Wilson elements by considering the computational complexity for three-dimensional (3D) problems considered. In general, the convergence of Wilson elements will heavily rely on the quality of the meshes. It will greatly deteriorate or no longer converge when the mesh distortion is very large. In this paper, the refined element method based on Wilson element is first applied to solve nearly incompressible elasticity problems, and the influence of mesh quality on the refined element is tested numerically. Its validity is verified by some numerical examples. By using the internal condensation method, the refined element discrete system of equations is deduced into the one which is spectrally equivalent to an 8-node hexahedral element discrete system of equations. And then, a type of efficient algebraic multigrid (AMG) preconditioner is presented by combining both the coarsening techniques based on the distance matrix and the effective smoothing operators. The resulting preconditioned conjugate gradient (PCG) method is efficient for 3D nearly incompressible problems. The numerical results verify the efficiency and robustness of the proposed method.