Next: Numerical procedure Up: Three-dimensional scattering of elastic Previous: Integral equation

## FM-BIEM

In FM-BIEM our first step is to expand the fundamental solution (3.139) into a series of products of functions of and those of . To this end we notice that the fundamental solution (3.139) can be rewritten as

 (3.141)

where we have used the following identity:

Using (3.120) and (3.142), we obtain a series expansion of the fundamental solution as follows:

where the function is given by (3.121) and the functions and are defined as

 = (3.142) = (3.143)

In (3.144) and (3.145) the function is given by (3.122).

Now we compute the integral on the right-hand side of (3.140) over a subset of S denoted by Sy for an which is away from Sy. Using (3.143) we obtain

 = (3.144)

where MTr,n,m(O) and MLn,m(O) are the multipole moments centred at O defined as

 MTr,n,m(O) = (3.145) MLn,m(O) = (3.146)

In (3.146) we have assumed that the inequality ( ) holds. Notice that MTr,n,m(O) has three components and MLn,m(O) has one component for a fixed pair of nand m and therefore the number of the multipole moments is 4. The number of the multipole moments is thus seen to be the same as that in three-dimensional elastostatics. Fujiwara [22] uses FMM to compute the double-layer potential in three-dimensional elastodynamics and in Fujiwara's formulation the number of the multipole moments is 6. If one adopts our formulation, the number of the multipole moments for the double-layer potential is 4 (details of the difference between our formulation and Fujiwara's one are given later). Also, Fukui and Inoue [24] use 4 multipole moments in FM-BIEM for ordinary problems in two-dimensional elastodynamics. However, replacing erspwith e3sp in (3.147), one can find that the number of the multipole moments is 2 in our formulation for two-dimensional elastodynamics.

The multipole moments are translated according to the following formulae as the centre of the multipole expansion is shifted from O to O':

 (3.147)

 (3.148)

where we have used (3.123).

One can evaluate the integral on the right-hand side of (3.140) with the local expansion in the following manner:

 = (3.149)

where and are the coefficients of the local expansion centred at defined as

 (3.150)

 (3.151)

and the functions and are obtained by replacing with in (3.144) and (3.145). In the derivation of (3.152) and (3.153) we have used (3.124) and have assumed that the inequality holds.

The coefficients of the local expansion are translated according to the following formulae as the centre of the local expansion is shifted from to :

 (3.152)

 (3.153)

where we have used (3.123).

Notice that M2M translations in (3.149) and (3.150), M2L translations in (3.152) and (3.153) and L2L translations in (3.154) and (3.155) have the same form as the M2M translation in (3.129), M2L translation in (3.131) and L2L translation in (3.132) in Helmholtz's equation. Also, the derivation of translation formulae in elastodynamics is similar to that in Helmholtz's equation.

We have thus prepared all the formulae needed for FM-BIEM. Using formulae presented above and the algorithm described in chapter 2, one can implement FM-BIEM.

Next: Numerical procedure Up: Three-dimensional scattering of elastic Previous: Integral equation
Ken-ichi Yoshida
2001-07-28