Next: FM-GBIEM Up: Crack problems for three-dimensional Previous: Crack problems for three-dimensional

## Variational Boundary Integral Equation

Multiplying (3.46) by a test function which satisfies the following condition:

and then integrating the resulting expressions on S, one obtains

 (3.100)

The non-regularised variational equation in (3.100) is used as the basis for the Galerkin BIEM.

Note that an arbitrary equilibrated stress field possesses an expression in terms of a stress function. Specifically, the stress'' computed from the fundamental solution can be expressed in terms of the stress function as (See Nishimura and Kobayashi [63])

 (3.101)

where eijk is the permutation symbol and the function is given by
 (3.102)

With (3.100) and (3.101), one derives the following variational equation given by (See Appendix L.3)

 (3.103)

The variational equation in (3.103) has been considered useful in the standard Galerkin BIEM since the kernel function is regularised to an integrable one. In the present investigation with FMM, however, the variational equation in (3.100) is more useful. It is because we are interested in applying FMM in evaluating contributions to the integrals in (3.100) or (3.103) from distant and , and because the kernel function in (3.100) decays more quickly than that in (3.103) as . Actually, FMM formulations with regularised integral equations have been tested in section 3.3 with collocation and were shown to introduce more moments than formulations without regularisation. Also, the numerical experiments in section 3.3 showed that the FMM with regularised integral equations is slower than that without regularisation. Notice, however, that the following identity holds for arbitrary smooth functions and having supports within S:

 = (3.104)

We shall use the regularised expression on the right-hand side of (3.104) in an auxiliary manner to compute the integral on the left-hand side of (3.104) with and , respectively, taken as one of the basis functions (sometimes referred to as global shape functions) for the crack opening displacement. This calculation will be needed in the 5th step of the FMM algorithm in 3.4.3 when the supports of and are not separated sufficiently.

Next: FM-GBIEM Up: Crack problems for three-dimensional Previous: Crack problems for three-dimensional
Ken-ichi Yoshida
2001-07-28