000015657 001__ 15657
000015657 005__ 20161115135326.0
000015657 04107 $$aeng
000015657 046__ $$k2013-06-12
000015657 100__ $$aChaillat, S.
000015657 24500 $$aFast Multipole Accelerated Boundary Element Method for Problems in an Elastic Half-Space

000015657 24630 $$n34.$$pComputational Methods in Structural Dynamics and Earhquake Engineering
000015657 260__ $$bNational Technical University of Athens, 2013
000015657 506__ $$arestricted
000015657 520__ $$2eng$$aContext. This contribution is concerned with the development of the Fast Multipole accelerated Boundary Element Method (FM-BEM) to simulate the propagation of elastic waves in semi-infinite media. In [1], we previously developed a FM-BEM based on the elastic full-space fundamental solution. This method has been shown to be very efficient to simulate elastic wave in semi-infinite media, but nevertheless suffers from the drawback of requiring a discretization of the free surface. In practice, the free-surface was truncated at a chosen radius "large enough" for achieving a good accuracy in the target domain of study. A large number of BEM degrees of freedom (DOFs) was thus required on the free-surface for the sole purpose of enforcing the traction-free boundary condition. To avoid the truncation issue, one can use a fundamental solution that intrinsically satisfies a tractionfree boundary condition on the free-surface. This formulation thus reduces the overall size of the BE model since the free-surface no longer requires discretization. The derivation and implementation of the corresponding fundamental solution [3] are involved. In particular, unlike its full-space counterpart, the half-space fundamental solution cannot be expressed in terms of simpler kernels (e.g. Laplace or Helmholtz fundamental solutions) having already-known multipole expansions. Multipole expansions of the elastic half-space fundamental solution thus cannot be obtained in a simple way. Methodology. In [2], we proposed an expansion of the elastic half-space fundamental solution in a form which achieves the separation of variables required by the Fast Multipole Method. This expansion is based on a decomposition of the fundamental solution and of the single-layer potential into three terms: the full-space fundamental solution, the image full-space fundamental solution and a complementary term to satisfy the traction-free condition on the free-surface. In this communication, we present a new FM-BEM based on this treatment of the half-space fundamental solution. Numerical efficiency of the new FM-BEM. The accuracy of this new FM-BEM (which does not require meshing the free surface) is compared to that of the FM-BEM based on the elastic full-space fundamental solutions (which does require meshing the free surface). Moreover, the numerical efficiency of both approaches is compared on seismology-oriented examples such as the scattering of plane waves by a cavity embedded in an elastic half-space. Despite the additional computational effort required by the evaluation of the proposed FM-capable form of the half-space fundamental solution, the new approach is shown to reduce several-fold the overall analysis time, thus establishing the overall benefit brought by removing the free-surface mesh.

000015657 540__ $$aText je chráněný podle autorského zákona č. 121/2000 Sb.
000015657 653__ $$a

000015657 7112_ $$aCOMPDYN 2013 - 4th International Thematic Conference$$cIsland of Kos (GR)$$d2013-06-12 / 2013-06-14$$gCOMPDYN2013
000015657 720__ $$aChaillat, S.$$iBonnet, M.
000015657 8560_ $$ffischerc@itam.cas.cz
000015657 8564_ $$s204730$$uhttps://invenio.itam.cas.cz/record/15657/files/1210.pdf$$yOriginal version of the author's contribution as presented on CD, section: CD-MS 24 ADVANCES IN MODELING OF WAVE PROPAGATION AND APPLICATIONS
.
000015657 962__ $$r15525
000015657 980__ $$aPAPER