An efficient BE iterative-solver-based substructuring algorithm for 3D time-harmonic problems in elastodynamics

Abstract This work is concerned with the development of an efficient and general algorithm to solve frequency-domain problems modelled by the boundary element method based on a sub-region technique. A specific feature of the algorithm discussed here is that the global sparse matrix of the coupled system is implicitly considered, i.e. problem quantities are not condensed into interface variables. The proposed algorithm requires that only the block matrices with non-zero complex-valued coefficients be stored and manipulated during the analysis process. In addition, the efficiency of the technique presented is improved by using iterative solvers. The good performance of pre-conditioned iterative solvers for systems of equations having real-valued coefficients, well demonstrated in the literature, is confirmed for the present case where the system matrix coefficients are complex. The efficiency of the algorithm described here is verified by analysing a soil–machine foundation interaction problem. CPU time and accuracy are the parameters used for estimating the computational efficiency.

[1]  Webe João Mansur,et al.  Determination Of 3D Time Domain Responses In Layered Media By Using A Coupled BE/FE Process , 1998 .

[2]  M. H. Aliabadi,et al.  Out-of-core solver for large, multi-zone boundary element matrices , 1995 .

[3]  J. Telles,et al.  Multi-level hierarchical preconditioners for boundary element systems , 1993 .

[4]  R. A. Białecki,et al.  IN‐ AND OUT‐OF‐CORE BEM EQUATION SOLVER WITH PARALLEL AND NON‐LINEAR OPTIONS , 1996 .

[5]  Iterative methods for solving boundary element equations , 1987 .

[6]  J. H. Kane.  Boundary Element Analysis in Engineering Continuum Mechanics , 1994 .

[7]  Herbert A. Mang,et al.  A new method for evaluating singular integrals in stress analysis of solids by the direct boundary element method , 1985 .

[8]  J.C.F. Telles,et al.  Iterative solution of bem equations by GMRES algorithm , 1992 .

[9]  Kang C. Jea,et al.  Generalized conjugate-gradient acceleration of nonsymmetrizable iterative methods , 1980 .

[10]  J. M. Crotty.  A block equation solver for large unsymmetric matrices arising in the boundary integral equation method , 1982 .

[11]  David E. Keyes,et al.  Preconditioned Krylov solvers for BEA , 1994 .

[12]  Dimitri E. Beskos,et al.  Boundary Element Methods in Dynamic Analysis , 1987 .

[13]  F. Durbin,et al.  Numerical Inversion of Laplace Transforms: An Efficient Improvement to Dubner and Abate's Method , 1974, Comput. J..

[14]  Economical solution technique for boundary integral matrices , 1983 .

[15]  Γεωργιοσ Μανωλησ.  DYNAMIC RESPONSE OF UNDERGROUND STRUCTURES , 1980 .

[16]  O. Axelsson,et al.  Finite element solution of boundary value problemes - theory and computation , 2001, Classics in applied mathematics.

[17]  Webe João Mansur,et al.  Solution of BEM systems of equations via iterative techniques , 1992 .

[18]  D. Beskos,et al.  Boundary Element Methods in Elastodynamics , 1988 .

[19]  J. H. Kane,et al.  A fully symmetric multi‐zone Galerkin boundary element method , 1999 .

[20]  J. Wolf.  Dynamic soil-structure interaction , 1985 .

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

[22]  A. Papoulis.  A new method of inversion of the Laplace transform , 1957 .

[23]  George D. Manolis,et al.  A comparative study on three boundary element method approaches to problems in elastodynamics , 1983 .

[24]  W. Hackbusch.  Iterative Lösung großer schwachbesetzter Gleichungssysteme , 1991 .

[25]  Sunil Saigal,et al.  An arbitrary condensing, noncondensing solution strategy for large scale, multi-zone boundary element analysis , 1990 .

[26]  Dimitri E. Beskos,et al.  Boundary Element Methods in Dynamic Analysis: Part II (1986-1996) , 1997 .