×

Large-scale vectorized implicit calculations in solid mechanics on a Cray X-MP/48 utilizing EBE preconitioned conjugate gradients. (English) Zbl 0606.73096

The major costs of large implicit finite element calculations, particularly in three dimensions, arise from computing solutions to systems of linear equations. Direct methods, i.e., those based upon Gaussian elimination, can easily require prohibitively large amounts of both CPU time and storage, even on current supercomputers. Iterative procedures avoiding the formation and factorization of a global system of equations can circumvent these difficulties. The element-by-element (EBE) preconditioned conjugate gradients (PCG) algorithm is presented in the context of a vectorized implementation within the production nonlinear stress analysis code Nike3D. Due to continued confusion as to the ease of vectorizing finite element procedures, we include examples of the main EBE subroutines in their entirety.
The concept of a fractal dimension of a finite element mesh is introduced, and proves useful in characterizing the efficiency of this iterative algorithm with respect to a variable band, active column direct method. Sample calculations on a Cray X-MP/48 with solid-state storage device (SSD) illustrate the economy and range of applicability of EBE/PCG. Asymptotic cost formulae derived for two linear problems underscore differences between the direct and iterative algorithms for large problems and lead to predictions of problem size limitations imposed by the computing environment.

MSC:

74S05 Finite element methods applied to problems in solid mechanics
74-04 Software, source code, etc. for problems pertaining to mechanics of deformable solids
65Yxx Computer aspects of numerical algorithms
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] (CRAY X-MP series model 48 mainframe reference manual. CRAY X-MP series model 48 mainframe reference manual, Publication No. HR-0097 (1984), Cray Research, Inc: Cray Research, Inc Bloomington, MN)
[2] (CFT: The CRAY-1 FORTRAN compiler. CFT: The CRAY-1 FORTRAN compiler, Publication No. SR-0009 (1984), Cray Research, Inc: Cray Research, Inc Bloomington, MN)
[3] Golub, G. H.; Van Loan, C. F., Matrix Computations (1983), The Johns Hopkins University Press: The Johns Hopkins University Press Baltimore, MD · Zbl 0559.65011
[4] Goudreau, G. L.; Hallquist, J. O., Recent developments in large-scale Lagrangian hydrocode technology, Comput. Meths. Appl. Mech. Engrg., 33, 725-757 (1982) · Zbl 0493.73072
[5] Goudreau, G. L.; Bailey, R. A.; Hallquist, J. O.; Murray, R. C.; Sackett, S. J., Efficient large-scale finite element calculations in a CRAY environment, (UCRL-89385 (1984), University of California)
[6] Goudreau, G. L., Large Scale Computations, (Nemat-Nasser, S.; Asoro, R. J.; Hegemier, G. A., Theoretical Foundations for Large-Scale Computations for Nonlinear Material Behavior (1984), Martinus Nijhoff: Martinus Nijhoff Dordrecht, The Netherlands), 65-91
[7] Hageman, L. A.; Young, D. M., Applied Iterative Methods (1981), Academic Press: Academic Press New York · Zbl 0459.65014
[8] Hallquist, J. O., User’s manual for DYNA2D—An explicit two-dimensional hydrodynamic finite element code with interactive rezoning, (UCID-18756, Rev. 1 (1982), University of California)
[9] Hallquist, J. O., NIKE2D—A vectorized, implicit, finite deformation, finite element code for analyzing the static and dynamic response of 2-D solids, (UCID-19677 (1983), University of California) · Zbl 0319.70015
[10] Hallquist, J. O., NIKE3D: An implicit, finite deformation, finite element code for analyzing the static and dynamic response of three-dimensional solids, (UCID-18822, Rev. 1 (1984), University of California) · Zbl 0319.70015
[11] Hallquist, J. O.; Goudreau, G. L.; Benson, D. J., Sliding interfaces with contact-impact in large-scale Lagrangian computations, Comput. Meths. Appl. Mech. Engrg., 51, 107-135 (1985) · Zbl 0567.73120
[12] Hallquist, J. O.; Benson, D. J., DYNA3D user’s manual (nonlinear dynamic analysis of structures in three dimensions), (UCID-19592, Rev. 2 (1986), University of California)
[13] Hestenes, M. R.; Stiefel, E., Method of conjugate gradients for solving linear systems, J. Res. Nat. Bur. Standards, 49, 409-436 (1952) · Zbl 0048.09901
[14] Hughes, T. J.R.; Pister, K. S.; Taylor, R. L., Implicit-explicit finite elements in nonlinear transient analysis, Comput. Meths. Appl. Mech. Engrg., 17/18, 159-182 (1979) · Zbl 0413.73074
[15] Hughes, T. J.R.; Levit, I.; Winget, J. M., Implicit, unconditionally stable algorithms for heat conduction analysis, ASCE J. Engrg. Mech. Div., 109, 576-585 (1983)
[16] Hughes, T. J.R.; Levit, I.; Winget, J. M., An element-by-element solution algorithm for problems of structural and solid mechanics, Comput. Meths. Appl. Mech. Engrg., 36, 241-254 (1983) · Zbl 0487.73083
[17] Hughes, T. J.R.; Winget, J.; Levit, I.; Tezduyar, T., New alternating direction procedures in finite element analysis based upon EBE approximate factorizations, (Atluri, S. N.; Perrone, N., Computer Methods for Nonlinear Solids and Structural Mechanics, AMD-4 (1983), ASME: ASME New York), 75-109 · Zbl 0563.73052
[18] Hughes, T. J.R.; Raefsky, A.; Muller, A.; Winget, J. M.; Levit, I., A progress report on EBE solution procedures in solid mechanics, (Taylor, C.; Hinton, E.; Owen, D. R.J.; Oñate, E., Numerical Methods for Nonlinear Problems, Vol. 2 (1984), Pineridge Press: Pineridge Press Swansea, U.K), 18-26
[19] Hughes, T. J.R., Numerical implementation of constitutive models: Rate-independent deviatoric plasticity, (Nemat-Nasser, S.; Asaro, R. J.; Hegemier, G. A., Theoretical Foundations for Large-Scale Computations for Nonlinear Material Behavior (1984), Martinus Nijhoff: Martinus Nijhoff Dordrecht, The Netherlands), 29-57
[20] Lyzenga, G. A.; Raefsky, A., Finite elements and the method of preconditioned conjugate gradients on a concurrent processor, (Fox, J.; Lyzenga, G. A., The Cal Tech Concurrent Processor (1986)), appear.
[21] Newmark, N. M., A method of computation for structural dynamics, ASCE J. Engrg. Mech. Div., 85, 67-94 (1959)
[22] Nour-Omid, B., A preconditioned conjugate gradient method for solution of finite element equations, (Liu, W. K.; Belytschko, T.; Park, K. C., Innovative Methods for Nonlinear Problems (1984), Pineridge Press: Pineridge Press Swansea, U.K), 17-40
[23] Ortiz, M.; Pinsky, P. M.; Taylor, R. L., Unconditionally stable element-by-element algorithms for dynamics problems, Comput. Meths. Appl. Mech. Engrg., 36, 223-239 (1983) · Zbl 0501.73068
[24] Taylor, R. L.; Wilson, E. L.; Sackett, S. J., Direct solution of equations by frontal and variable band, active column methods, (Wunderlich, W.; Stein, E.; Bathe, K. J., Nonlinear Finite Element Analysis in Structural Mechanics (1981), Springer: Springer Berlin), 521-552 · Zbl 0457.73064
[25] Varga, R. S., Matrix Iterative Analysis (1962), Prentice-Hall: Prentice-Hall Englewood Cliffs, NJ · Zbl 0133.08602
[26] Winget, J. M.; Hughes, T. J.R., Solution algorithms for nonlinear transient heat conduction analysis employing element-by-element iterative strategies, Comput. Meths. Appl. Mech. Engrg., 52, 711-815 (1985) · Zbl 0579.73119
[27] Young, D. M., Iterative Solution of Large Linear Systems (1971), Academic Press: Academic Press New York · Zbl 0204.48102
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.