A Fast Multipole Method formulation for 3D elastodynamics in the frequency domain

The solution of the elastodynamic equations using boundary element methods (BEMs) gives rise to fully-populated matrix equations. Earlier investigations on the Helmholtz and Maxwell equations have established that the Fast Multipole (FM) method reduces the complexity of a BEM solution to N log₂ N per GMRES iteration. The present Note address the extension of the FM-BEM strategy to 3D elastodynamics in the frequency domain. Its efficiency and accuracy are demonstrated on numerical examples involving up to N = O (10⁶) nodal unknowns. To cite this article: S. Chaillat et al., C. R. Mecanique 335 (2007).