[Calcul de la fonction de Green de l'élasticité dans un demi-plan avec une condition aux limites d'impédance]
Dans cette Note, nous présentons une méthode numérique effective et précise de calcul de la fonction de Green G du système de l'élasticité en régime harmonique. Le problème est posé dans un demi-plan et les conditions aux limites sur la frontière de ce demi-plan sont du type impédance. Ce problème se rencontre dans l'exploitation des mines souterraines et en séismologie. La méthode de calcul est inspirée du travail de Duran et co-auteurs (2005) qui traite du calcul de la fonction de Green du problème de Helmholtz dans un demi-plan. Elle consiste à utiliser la transformée de Fourier partielle en la variable tangentielle, ce qui permet un calcul explicite de la fonction de Green spectrale . Pour calculer sa transformée de Fourier inverse, nous décomposons en une somme de deux termes. Le premier est associé à l'espace entier, tandis que le second prend en compte le demi-plan et les conditions aux limites. Le premier terme n'est autre que la fonction de Green du plan pour le système de l'élasticité en régime harmonique (cf. la thèse de J. Dompierre). Le second terme se sépare en trois termes, dont deux contiennent les singularités de la fonction spectrale (les pôles et les pseudo-pôles) tandis que le dernier est régulier et décroissant à l'infini. Nous calculons analytiquement la transformée de Fourier des deux termes singuliers et celle du dernier terme est obtenue par un algorithme type FFT. Nous présentons un résultat numérique. En outre, nous montrons que, sous certaines conditions, il est possible d'avoir une quatrièmme lenteur, qui pourrait propager une nouvelle onde de surface.
This Note presents an effective and accurate method for numerical calculation of the Green's function G associated with the time harmonic elasticity system in a half-plane, where an impedance boundary condition is considered. The need to compute this function arises when studying wave propagation in underground mining and seismological engineering. To theoretically obtain this Green's function, we have drawn our inspiration from the paper by Durán et al. (2005), where the Green's function for the Helmholtz equation has been computed. The method consists in applying a partial Fourier transform, which allows an explicit calculation of the so-called spectral Green's function. In order to compute its inverse Fourier transform, we separate as a sum of two terms. The first is associated with the whole plane, whereas the second takes into account the half-plane and the boundary conditions. The first term corresponds to the Green's function of the well known time-harmonic elasticity system in (cf. J. Dompierre, Thesis). The second term is separated as a sum of three terms, where two of them contain singularities in the spectral variable (pseudo-poles and poles) and the other is regular and decreasing at infinity. The inverse Fourier transform of the singular terms are analytically computed, whereas the regular one is numerically obtained via an FFT algorithm. We present a numerical result. Moreover, we show that, under some conditions, a fourth additional slowness appears and which could produce a new surface wave.
Accepté le :
Publié le :
Mot clés : Mécanique des solides numérique, Fonction de Green spectrale, Séismologie
@article{CRMECA_2006__334_12_725_0,
author = {Mario Dur\'an and Eduardo Godoy and Jean-Claude N\'ed\'elec},
title = {Computing {Green's} function of elasticity in a half-plane with impedance boundary condition},
journal = {Comptes Rendus. M\'ecanique},
pages = {725--731},
publisher = {Elsevier},
volume = {334},
number = {12},
year = {2006},
doi = {10.1016/j.crme.2006.09.003},
language = {en},
}
TY - JOUR AU - Mario Durán AU - Eduardo Godoy AU - Jean-Claude Nédélec TI - Computing Green's function of elasticity in a half-plane with impedance boundary condition JO - Comptes Rendus. Mécanique PY - 2006 SP - 725 EP - 731 VL - 334 IS - 12 PB - Elsevier DO - 10.1016/j.crme.2006.09.003 LA - en ID - CRMECA_2006__334_12_725_0 ER -
%0 Journal Article %A Mario Durán %A Eduardo Godoy %A Jean-Claude Nédélec %T Computing Green's function of elasticity in a half-plane with impedance boundary condition %J Comptes Rendus. Mécanique %D 2006 %P 725-731 %V 334 %N 12 %I Elsevier %R 10.1016/j.crme.2006.09.003 %G en %F CRMECA_2006__334_12_725_0
Mario Durán; Eduardo Godoy; Jean-Claude Nédélec. Computing Green's function of elasticity in a half-plane with impedance boundary condition. Comptes Rendus. Mécanique, Volume 334 (2006) no. 12, pp. 725-731. doi : 10.1016/j.crme.2006.09.003. https://comptes-rendus.academie-sciences.fr/mecanique/articles/10.1016/j.crme.2006.09.003/
[1] Acoustic and Electromagnetic Equations. Integral Representations for Harmonic Problems, Applied Mathematical Sciences, vol. 144, Springer-Verlag, Berlin/New York, 2001
[2] Integral Equations Methods in Scattering Theory, John Wiley, New York, 1983
[3] Boundary Integral Equation Methods for Solids and Fluids, John Wiley & Sons Ltd., Chichester, 1995
[4] Boundary Integral Equations in Elasticity Theory, Kluwer Academic Publishers, Dordrecht, Boston, 2002
[5] The Helmholtz equation with impedance in a half-plane, C. R. Acad. Sci. Paris, Ser. I, Volume 340 (2005), pp. 483-488
[6] Tables of Integral Transformations, vol. I, McGraw–Hill Book Company, Inc., New York, 1954
[7] Special Functions for Scientists and Engineers, Dover Publications, Inc., New York, 1968
[8] J. Dompierre, Équations intégrals en axisymétrie généralisée, application à la sismique entre puits, Thesis, École Centrale de Paris, 1993
Cité par Sources :
Commentaires - Politique
Radiation condition and uniqueness for the outgoing elastic wave in a half-plane with free boundary
Mario Durán; Ignacio Muga; Jean-Claude Nédélec
C. R. Math (2009)
The Helmholtz equation with impedance in a half-space
Mario Durán; Ignacio Muga; Jean-Claude Nédélec
C. R. Math (2005)
The Helmholtz equation with impedance in a half-plane
Mario Durán; Ignacio Muga; Jean-Claude Nédélec
C. R. Math (2005)