Fast multipole method applied to elastostatic BEM-FEM coupling

Authors:
Massimiliano Margonari;Marc Bonnet
Affiliations:
Department of Mechanical and Structural Engineering, University of Trento, Via Mesiano 77, 38050 Trento, Italy;Laboratoire de Mécanique des Solides (UMR 7649 CNRS), ícole Polytechnique, F-91128 Palaiseau cedex, France
Venue:
Computers and Structures
Year:
2005

Citing 8
Cited 3

GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems

SIAM Journal on Scientific and Statistical Computing
A new method for the coupling of finite element and boundary element discretized subdomains of elastic bodies

Computer Methods in Applied Mechanics and Engineering
A fast algorithm for particle simulations

Journal of Computational Physics
Parallel Preconditioning with Sparse Approximate Inverses

SIAM Journal on Scientific Computing
Iterative methods for solving linear systems

Iterative methods for solving linear systems
A Sparse Approximate Inverse Preconditioner for Nonsymmetric Linear Systems

SIAM Journal on Scientific Computing
Programming the Boundary Element Method

Programming the Boundary Element Method
Iterative Methods for Sparse Linear Systems

Iterative Methods for Sparse Linear Systems

Speeding-up Finite Element analyses by replacing the linear equation solver with a Boundary Element code. Part 1: 2D linear elasticity

Computers and Structures
The fast multipole boundary element methods (FMBEM) and its applications in rolling engineering analysis

Computational Mechanics
A fast elasto-plastic formulation with hierarchical matrices and the boundary element method

Computational Mechanics

Quantified Score

Hi-index	0.00

Visualization

Abstract

BEM-FEM coupling is desirable for three-dimensional problems involving specific features such as (i) large or unbounded media with linear constitutive properties, (ii) cracks, (iii) critical parts of complex geometry requiring accurate stress analyses. However, for cases with a BEM discretization involving a large number N"B"E"M of degrees of freedom, setting up the BEM contribution to the coupled problem using conventional techniques is an expensive O(N"B"E"M^2) task. Moreover, the fully-populated BEM block entails a O(N"B"E"M^2) storage requirement and a O(N"B"E"M^3) contribution to the solution time via usual direct solvers. To overcome these pitfalls, the BEM contribution is formulated using the fast multipole method (FMM) and the coupled equations are solved by means of an iterative GMRES solver. Both the storage requirements and the solution times are found to be close to O(N"B"E"M). A preconditioner based on the sparse approximate inverse of the BEM block is shown to improve the convergence of the GMRES solver. Numerical examples involving N"B"E"M=O(10^5-10^6) unknowns, run on a PC computer, are presented; they include the Eshelby inclusion (as a validation example), a many-inclusion configuration, and a dam structure.