The conventional frequency-domain boundary element method(BEM)based on discrete Fourier transform is difficult to complete transient analysis of low damped and undamped systems. To solve this problem
the exponential window method is combined with the frequency-domain BEM
and the precorrected fast Fourier transform(pFFT)is used to accelerate the frequency-domain BEM analysis. For solving the sequence of linear systems corresponding to the sampling frequencies with higher efficiency
an extrapolation method is proposed to obtain good initial guesses. In addition
a newly developed Krylov subspace recycling method is employed to solve this sequence of linear systems. Numerical examples demonstrate a dramatic reduction of iterations
which leads to higher efficiency and smaller memory consuming.
BANJAI L, SAUTER S. Rapid solution of the wave equation in unbounded domains[J]. SIAM Journal on Numerical Analysis, 2008, 47(1): 227-249.
CHAILLAT S, BONNET M, SEMBLAT J F. A multi-level fast multipole BEM for 3-D elastodynamics in the frequency domain[J]. Computer Methods in Applied Mechanics and Engineering, 2008, 197(49): 4233-4249.
XIAO Jinyou, YE Wenjing, CAI Yaxiong, et al. Precorrected FFT accelerated BEM for large-scale transient elastodynamic analysis using frequency-domain approach[J]. International Journal for Numerical Methods in Engineering, 2012, 90(1): 116-134.
PARKS M L, STURLER E, MACKEY G, et al. Recycling Krylov subspaces for sequences of linear systems[J]. SIAM Journal on Scientific Computing, 2006, 28(5): 1651-1674.
KAUSEL E, ROESSET J M. Frequency domain analysis of undamped systems[J]. Journal of Engineering Mechanics, 1992, 118(4): 721-734.
SAAD Y, MARTIN H S. GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems[J]. SIAM Journal on Scientific Computing, 1986, 7(3): 856-869.