Comptes Rendus
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.

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 Gˆ 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 R2 (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.

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 Gˆ. Pour calculer sa transformée de Fourier inverse, nous décomposons Gˆ 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.

Received:
Accepted:
Published online:
DOI: 10.1016/j.crme.2006.09.003
Keywords: Computational solid mechanics, Spectral Green's function, Seismology
Mot clés : Mécanique des solides numérique, Fonction de Green spectrale, Séismologie

Mario Durán 1; Eduardo Godoy 1; Jean-Claude Nédélec 2

1 Facultad de Ingeniería, Pontificia Universidad Católica de Chile, Casilla 306, Santiago 22, Chile
2 CMAP, École polytechnique, 91128 Palaiseau cedex, France
@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] J.-C. Nédélec Acoustic and Electromagnetic Equations. Integral Representations for Harmonic Problems, Applied Mathematical Sciences, vol. 144, Springer-Verlag, Berlin/New York, 2001

[2] P. Colton; R. Kress Integral Equations Methods in Scattering Theory, John Wiley, New York, 1983

[3] M. Bonnet Boundary Integral Equation Methods for Solids and Fluids, John Wiley & Sons Ltd., Chichester, 1995

[4] A.M. Linkov Boundary Integral Equations in Elasticity Theory, Kluwer Academic Publishers, Dordrecht, Boston, 2002

[5] M. Durán; I. Muga; J.-C. Nédélec The Helmholtz equation with impedance in a half-plane, C. R. Acad. Sci. Paris, Ser. I, Volume 340 (2005), pp. 483-488

[6] H. Bateman Tables of Integral Transformations, vol. I, McGraw–Hill Book Company, Inc., New York, 1954

[7] W.W. Bell 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

Cited by Sources:

Comments - Policy