A full 3D finite element analysis using adaptive refinement and PCG solver with back interpolation

In this paper, adaptive refinement finite element analyses were carried out for full 3D problems. In order to achieve an optimal computation cost and to eliminate the effects of singular points, the adaptive refinement procedure was applied in conjunction with a back interpolation for the construction of a good initial guess for the solution of the linear equations system. The combined use of the adaptive refinement procedure, the back interpolation scheme and the preconditioned conjugate gradient (PCG) solver lead to a significant reduction in the operations for the solution of the simultaneous linear equations in the adaptive refinement analysis. In many cases the number of iterations needed by the PCG solver to reach a converged solution is independent of the number of equations in the global system. The numerical results obtained indicated that for a series of carefully designed adaptive meshes, the computational cost required to solve the linear system could be made only proportional to the number of degrees of freedom in the mesh. To our knowledge, this is the first time that the global stiffness equations in 3D stress analysis have been solved with this efficiency. The same set of numerical results also showed that, in general, the total computational cost of the adaptive refinement procedure can usually be minimized by gradually reducing the target relative error of the solution during successive refinements rather than employing a constant target relative error throughout the whole adaptive analysis.

[1]  Kent L. Lawrence,et al.  A study of adaptively remeshed finite element problems using higher order tetrahedra , 1995 .

[2]  J. Bonet,et al.  An alternating digital tree (ADT) algorithm for 3D geometric searching and intersection problems , 1991 .

[3]  J. Hall,et al.  The multigrid method in solid mechanics: Part I—Algorithm description and behaviour , 1990 .

[4]  R. E. Hobbs,et al.  On using different finite elements with an automatic adaptive refinement procedure for the solution of 2‐D stress analysis problems , 1997 .

[5]  Rainald Löhner,et al.  Automatic unstructured grid generators , 1997 .

[6]  S. H. Lo,et al.  3D mesh refinement in compliance with a specified node spacing function , 1998 .

[7]  Kincho H. Law,et al.  A robust incomplete factorization based on value and space constraints , 1995 .

[8]  J. Z. Zhu,et al.  The superconvergent patch recovery and a posteriori error estimates. Part 1: The recovery technique , 1992 .

[9]  Nils-Erik Wiberg,et al.  Patch recovery based on superconvergent derivatives and equilibrium , 1993 .

[10]  Chi King Lee,et al.  On using different recovery procedures for the construction of smoothed stress in finite element method , 1998 .

[11]  S. Timoshenko,et al.  Theory of elasticity , 1975 .

[12]  Ted Belytschko,et al.  Superconvergent patch recovery with equilibrium and conjoint interpolant enhancements , 1994 .

[13]  I. D. Parsons,et al.  The multigrid method in solid mechanics: Part II—Practical applications , 1990 .

[14]  Chi King Lee,et al.  Automatic adaptive refinement finite element procedure for 3D stress analysis , 1997 .

[15]  Taeoh Lee,et al.  A SUPERCONVERGENT STRESS RECOVERY TECHNIQUE WITH EQUILIBRIUM CONSTRAINT , 1997 .

[16]  Robert Beauwens,et al.  HIGH-PERFORMANCE PCG SOLVERS FOR FEM STRUCTURAL ANALYSIS , 1996 .

[17]  P. Hood,et al.  Frontal solution program for unsymmetric matrices , 1976 .

[18]  Sh Lo,et al.  On robust implementation of the superconvergent patch recovery technique , 1993 .

[19]  Bruce M. Irons,et al.  A frontal solution program for finite element analysis , 1970 .

[20]  J. Peraire,et al.  Unstructured grid methods for compressible flows , 1992 .

[21]  O. C. Zienkiewicz,et al.  A simple error estimator and adaptive procedure for practical engineerng analysis , 1987 .

[22]  P. George Improvements on Delaunay-based three-dimensional automatic mesh generator , 1997 .

[23]  Chi King Lee,et al.  AUTOMATIC ADAPTIVE 3‐D FINITE ELEMENT REFINEMENT USING DIFFERENT‐ORDER TETRAHEDRAL ELEMENTS , 1997 .

[24]  Carlos A. Felippa,et al.  Solution of linear equations with skyline-stored symmetric matrix , 1975 .

[25]  Richard H. Crawford,et al.  Mesh rezoning of 2D isoparametric elements by inversion , 1989 .

[26]  Eliezer Mendelssohn,et al.  Solution of linear equations with a symmetrically skyline-stored nonsymmetric matrix , 1984 .

[27]  O. Axelsson Iterative solution methods , 1995 .

[28]  Graham H. Powell,et al.  Towards optimal in-core equation solving , 1974 .

[29]  Ivo Babuška,et al.  On the Rates of Convergence of the Finite Element Method , 1982 .

[30]  Bijan Boroomand,et al.  RECOVERY BY EQUILIBRIUM IN PATCHES (REP) , 1997 .

[31]  Chi King Lee,et al.  On constructing accurate recovered stress fields for the finite element solution of Reissner-Mindlin plate bending problems , 1998 .

[32]  J. Peraire,et al.  UNSTRUCTURED TETRAHEDRAL MESH GENERATION FOR THREE-DIMENSIONAL VISCOUS FLOWS , 1996 .

[33]  O. C. Zienkiewicz,et al.  The Finite Element Method: Basic Formulation and Linear Problems , 1987 .

[34]  I. Duff,et al.  The effect of ordering on preconditioned conjugate gradients , 1989 .

[35]  M. Papadrakakis,et al.  Accuracy and effectiveness of preconditioned conjugate gradient algorithms for large and ill-conditioned problems , 1993 .

[36]  Ivo Babuška,et al.  A MODEL STUDY OF THE QUALITY OF A POSTERIORI ERROR ESTIMATORS FOR FINITE ELEMENT SOLUTIONS OF LINEAR ELLIPTIC PROBLEMS, WITH PARTICULAR REFERENCE TO THE BEHAVIOR NEAR THE BOUNDARY , 1997 .

[37]  Ivo Babuška,et al.  A posteriori error analysis and adaptive processes in the finite element method: Part II—adaptive mesh refinement , 1983 .

[38]  R. E. Hobbs,et al.  Automatic adaptive refinement for shell analysis using nine-node assumed strain element , 1997 .

[39]  Thomas J. R. Hughes,et al.  LARGE-SCALE VECTORIZED IMPLICIT CALCULATIONS IN SOLID MECHANICS ON A CRAY X-MP/48 UTILIZING EBE PRECONDITIONED CONJUGATE GRADIENTS. , 1986 .