2006 Volume 62 Issue 4 Pages 936-949
A method for the analysis of elastic wave field is presented by means of the volume integral equation. The usage of the unitarity of the Fourier transform as well as the piecewise constant basis functions lead to a large scale non-Hermite sparse matrix from the integralequation. It is found from the numerical analysis that the Bi-CGSTAB method which is one of the iterative methods for the linear algebraic equation can be applied successfully to the matrix obtained from the integral equation. The numerical solutions of the integral equation showed good agreement with those from the surface spherical harmonics expansion. The behaviour of the scattred elastic waves are found to be well explained by the numerical solutions as well as the spectral structure of the scattered waves.