A transient FETI methodology for large‐scale parallel implicit computations in structural mechanics

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallellize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet and perhaps will never be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient than explicit codes when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

[1]  Charbel Farhat,et al.  An Unconventional Domain Decomposition Method for an Efficient Parallel Solution of Large-Scale Finite Element Systems , 1992, SIAM J. Sci. Comput..

[2]  H. V. D. Vorst,et al.  The rate of convergence of Conjugate Gradients , 1986 .

[3]  Charbel Farhat,et al.  Transient finite element computations on 65536 processors: The connection machine , 1990 .

[4]  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 .

[5]  T. Hughes,et al.  An element-by-element solution algorithm for problems of structural and solid mechanics , 1983 .

[6]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[7]  Charbel Farhat,et al.  A Lagrange multiplier based divide and conquer finite element algorithm , 1991 .

[8]  C. Farhat,et al.  A method of finite element tearing and interconnecting and its parallel solution algorithm , 1991 .

[9]  Charbel Farhat,et al.  Redesigning the Skyline Solver for Parallel/Vector Supercomputers , 1990, Int. J. High Speed Comput..

[10]  B. Parlett.  The Symmetric Eigenvalue Problem , 1981 .

[11]  Alan George,et al.  A linear time implementation of the reverse Cuthill-McKee algorithm , 1980, BIT.

[12]  Ted Belytschko,et al.  Finite element analysis on the connection machine , 1990 .

[13]  C. Felippa,et al.  A Hybrid Substructuring Method and an Adaptive Refinement Scheme for the Distributed Solution of Three-Dimensional Structural Problems , 1992 .

[14]  Åke Björck,et al.  Numerical Methods , 1995, Handbook of Marine Craft Hydrodynamics and Motion Control.

[15]  J. Z. Zhu,et al.  The finite element method , 1977 .

[16]  Charbel Farhat,et al.  A parallel active column equation solver , 1988 .

[17]  Jerome F. Hajjar,et al.  Parallel processing of central difference transient analysis for three-dimensional nonlinear framed structures , 1989 .

[18]  A. Booth.  Numerical Methods , 1957, Nature.

[19]  Charbel Farhat,et al.  A simple and efficient automatic fem domain decomposer , 1988 .

[20]  Charbel Farhat,et al.  Implicit time integration of a class of constrained hybrid formulations—Part I: Spectral stability theory , 1995 .

[21]  W. J. Stroud,et al.  An accurate nonlinear finite element analysis and test correlation of a stiffened composite wing panel , 1991 .

[22]  Charbel Farhat,et al.  Using a reduced number of Lagrange multipliers for assembling parallel incomplete field finite element approximations , 1992 .

[23]  O. Widlund,et al.  SOLVING ELLIPTIC PROBLEMS ON REGIONS PARTITIONED INTO SUBSTRUCTURES , 1984 .