Next: Rotation of coefficients Up: Crack problems for three-dimensional Previous: Crack problems for three-dimensional

## New FM-BIEM

Suppose that we evaluate the contribution from a cell Cs to the coefficients of the local expansion in a cell Ct contained in the interaction list of Cs. The cells Cs and Ct are called a source cell and a target cell, respectively. We divide the interaction list of Cs into 6 lists called the uplist, downlist, northlist, southlist, eastlist and westlist of Cs. The uplist and downlist contain target cells located in +x3 and -x3directions of Cs, respectively. The northlist and southlist contain target cells located in +x2 and -x2 directions of Csexcept those in the uplist or downlist, respectively. The eastlist and westlist contain target cells located in +x1 and -x1directions of Cs except those in the uplist, downlist, northlist or southlist, respectively. In this section we discuss the case where target cells are included in the uplist of Cs and each of the cells is a cube having a volume of d3.

Suppose that a source point is located at (y1,y2,y3) in Cs and a target point at (x1,x2,x3) in Ct. It is easy to see that the following integral representation[61] holds:

 (4.1)

where we have assumed that the inequality x3 > y3 holds. The inner integral with respect to in (4.1) is computed with the trapezoidal rule and the outer integral with generalised Gaussian quadrature rules[83], yielding

 (4.2)

where is given by

is the error term and the numbers , M(k), Gaussian weights and nodes are given in Yarvin and Rokhlin[83]. These numbers are available from Netlib (http://www.netlib.org/pdes/multipole/). One may determine these parameters considering the required accuracy.

Noting (4.2) and the following formulae: (See (3.19) and Appendix B)

 = (4.3) = (4.4) = (4.5)

one can evaluate the integral in (3.22) in the following manner: (See Appendix J.1)

where X(k,j;O) is the coefficient of the exponential expansion at O, defined by

 (4.6)

This formula converts the multipole moments into the exponential expansion coefficients. We call the procedure given by (4.6) M2X (Multipole moment to(2) eXponential expansion) translation''. The coefficients of the exponential expansion is translated according to the following formula when the centre of the exponential expansion is shifted from O to : (See Appendix J.2)

 (4.7)

where stand for the components of the vector . We call the procedure given by (4.7) X2X (eXponential expansion to(2) eXponential expansion) translation''.

We next need to convert the coefficients of the exponential expansion into the coefficients of the local expansion (See (3.27)). Noting that the integral in (3.22) is expressed with as follows:

 = (4.8)

and using (3.26) and the following formula (See (3.18) and Hobson[43]):

 (4.9)

we obtain the following formula which converts the coefficients of the exponential expansion into the coefficient of the local expansion: (See Appendix J.3)

 (4.10)

We call the procedure given by (4.10) X2L (eXponential expansion to(2) Local expansion) translation''. In the new FM-BIEM we replace M2L translation given by (3.27) with M2X, X2X and X2L translation given by (4.6), (4.7) and (4.10).

Next: Rotation of coefficients Up: Crack problems for three-dimensional Previous: Crack problems for three-dimensional
Ken-ichi Yoshida
2001-07-28