MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2016/2717810 2717810 Research Article Discussion on the Time-Harmonic Elastodynamic Half-Space Green’s Function Obtained by Superposition http://orcid.org/0000-0002-3370-2952 Yuan Huina 1 http://orcid.org/0000-0002-9040-9701 Pan Ziyang 1 Baek Seungik State Key Laboratory of Hydroscience and Engineering Department of Hydraulic Engineering Tsinghua University Beijing 100084 China tsinghua.edu.cn 2016 722016 2016 03 12 2015 29 01 2016 01 02 2016 722016 2016 Copyright © 2016 Huina Yuan and Ziyang Pan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The time-harmonic elastodynamic half-space Green’s function derived by Banerjee and Mamoon by way of superposition is discussed and examined against another semianalytical solution and a numerical solution. It is shown that Banerjee and Mamoon’s solution gives infinite z-displacement response when the depth of the source goes to infinity, which is unreasonable and does not agree with other solutions. A possible problem in the derivation is that it is inappropriate to directly extend the results of Mindlin’s superposition method for the elastostatic half-space problem to the dynamic case. The superposition of the six full-space elastodynamic solutions does not satisfy the required boundary conditions of the half-space elastodynamic problem as in the static case and thus does not solve the dynamic half-space problem.

1. Introduction

The elastodynamic Green’s function for the half-space is fundamental to the application of boundary element method (BEM) to situations involving semi-infinite media. Various derivations of the elastic displacement due to a subsurface transient or time-harmonic point force can be found in the literature . Often expressed in Fourier-Bessel integral forms, numerical evaluation and application of these solutions are usually complex and time-consuming . The time-harmonic elastodynamic half-space Green’s function under discussion, proposed by Banerjee and Mamoon , is derived by extending the superposition technique devised by Mindlin  for the elastostatic half-space problem to the dynamic case of a periodic point force in a semi-infinite solid.

2. Banerjee and Mamoon’s Green’s Function

Consider the situation where the periodic force is normal to the free boundary of the half-space, as depicted in Figure 1. The semi-infinite solid is bounded by the plane z=0, with the positive z-axis pointing to the interior of the body. A periodic force Peiωt is applied at point (0,0,+c) in the positive z-direction, where ω is the circular frequency and i is the imaginary unit. The displacements at point ξ(x,y,z) are to be found. The distance between ξ and the real source point is given by R1=r2+z-c21/2, whereas the distance between  ξ and the image point (0,0,-c) is denoted by R2=r2+z-c21/2.

Periodic force normal to the boundary in the interior of a semi-infinite solid (after Banerjee and Mamoon ).

According to Banerjee and Mamoon (B&M for short) , the solution to this problem is composed of six individual components, corresponding to solutions to six problems in an infinite solid. The first two problems represent single periodic forces at (0,0,+c) and (0,0,-c), respectively; the third, fourth, and sixth problems represent a dynamic double force, a dynamic center of compression, and a dynamic doublet  at (0,0,-c), respectively; and the fifth problem represents a line of dynamic centers of compression extending from z=-c to z=-. To facilitate the ensuing discussion and save space, here only the z-displacement component of the sixth problem is given:(1a)uz6=K6A6R2z2+B6,(1b)K6=-c22πμ3-4ν,(1c)A6=-15c22s2R22+15c2sR2+sR2c2+6e-sR2/c2R23+c22c1215c12s2R22+15c1sR2+s2R22c12+4sR2c1+9e-sR2/c1R23,(1d)B6=3c22s2R22+3c2sR2+1e-sR2/c2R23-c22c123c12s2R22+3c1sR2+sR2c1+2e-sR2/c1R23,where μ is the shear modulus, ν is the Poisson’s ratio, c1 and c2 are the pressure and shear wave velocities, respectively, and s=iω is the Laplace transformed parameter (note that, in the digital copy of the original paper, the square in (1a) is missing).

3. A Possible Problem and Comparison with Other Solutions

A close inspection on (1a)–(1d) reveals that, after substituting A6 and K6 into (1a), the third term in the second component of A6 yields(2)-c222πμc143-4ν·c2z+c2R23·e-sR2/c1.When the depth of the periodic point force goes to infinity, that is, c, this term also goes to infinity. Further examination shows that this term cannot be cancelled out when summing up all the six solutions, which means that the z-displacement at ξ(x,y,z) goes to infinity as c. However, the displacements at a given point ξ(x,y,z) with finite z are expected to vanish when c, since the distance between source and receiver goes to infinity.

To demonstrate the discrepancy between B&M’s solution and other solutions, a test problem is used. The semi-infinite solid is characterized by mass density ρ=2000 kg/m3, shear modulus μ=180 MPa, and Poisson’s ratio ν=0.2. The excitation frequency f=10 Hz, and thus the shear wavelength is λ2=30 m. The amplitude of the periodic force is P=1 N. The free surface z-displacement responses at different radii as the depth of the source changes from c=0 to c=10λ2 are calculated using B&M’s solution and a semi-analytical solution presented by Maurel et al. in . Both solutions are coded in MATLAB. The results are compared in Figure 2, where the z-displacement calculated using the former solution oscillates and tends to increase when c increases, while that obtained using the latter solution decays quickly with increasing c.

Comparison of free surface responses at different radii for c=0 to c=10λ2: (a) B&M’s solution and (b) solution of Maurel et al. .

For further comparison, a numerical solution is computed using the commercial FEM software ABAQUS. As shown in Figure 3, a 2D axisymmetric model of size 2000 m × 2000 m is built. The top surface is free; the bottom and side surfaces are fixed. A periodic force P=1 N of frequency f=10 Hz is applied along the axis of symmetry in the positive z-direction. The P wave speed is approximately 490 m/s. To avoid reflections from the fixed boundaries, the simulation time is set to be 4 s. The time-dependent free surface z-displacement responses at radii r=λ2 and r=10λ2 are plotted in Figure 4, from which it can be seen that the surface points begin to oscillate periodically and stably in a short time after the initial arrival of the wave. Therefore, the discrete inverse Fourier transform is applied to the time series between 2 s and 4 s to obtain the displacement response at 10 Hz.

The 2D axisymmetric model for numerical simulation.

The time-dependent free surface z-displacement responses at different radii.

Figure 5 compares the frequency-domain z-displacement responses along the free surface when the source is buried at different depths computed using these three methods. In general, the numerical solution agrees with Maurel’s solution, while B&M’s solution is different. And the difference increases with source depth. In , Banerjee and Mamoon compared their solution with those obtained by Whittaker and Christiano  and Kobayashi and Nishimura  for the time-harmonic Boussinesq problem; that is, c=0. The results for a0=ωr/c2=2πr/λ2 between 0 and 2.0 are shown in Figures 7–10 in  and good agreement is displayed. The position of a0=2.0 is denoted by a vertical line in Figure 5(a), from which it can be seen that when c=0, B&M’s solution agrees reasonably well with other solutions as a0<2.0 and deviates from other solutions as a0>2.0.

Comparison of free surface responses computed using different methods and for different source depths: (a) c=0, (b) c=λ2/4, (c) c=λ2/2, and (d) c=λ2.

4. Discussion

The above comparison indicates that Banerjee and Mamoon’s solution for the elastodynamic half-space problem is incorrect. The reason might be that the results of Mindlin’s superposition method for the elastostatic problem cannot be simply extended to the dynamic case.

For the static problem, the coefficients of the six solution components are determined by imposing the boundary conditions and equilibrium condition . The boundary conditions for the free surface z=0 are (3)σzzz=0=σzrz=0=0,where σzz and σzr are the normal and shear stresses. The equilibrium condition is given by(4)P=-02πrσzzdr,z>c,where P is the static point force applied at (0,0,+c).

For the dynamic case, the boundary conditions still hold. However, whether these boundary conditions can be satisfied is not explicitly stated in . Furthermore, the equilibrium condition given by (4) will be inapplicable to the dynamic case due to the presence of the inertia force that is associated with acceleration. To consider the inertia force caused by time-harmonic excitation, the equilibrium condition can be rewritten as(5)P=-02πrσzzdr-0z02πrρω2uzdrdz,z>c,where uz is the displacement in z-direction.

The boundary conditions are examined first. The expressions of the stresses are derived. Since these expressions are lengthy, only the stress components σzz and σzr on the free surface, that is, z=0, for the first and fifth problems are presented here:(6a)σzz1=λ+2μK1A1-B1c2R12+λK1C1-B1r2R12,(6b)K1=14πμ,(6c)A1=9c22s2R12+9c2sR1+sR1c2+4ce-sR1/c2R13-c22c129c12s2R12+9c1sR1+3ce-sR1/c1R13,(6d)B1=15c22s2R12+15c2sR1+sR1c2+6ce-sR1/c2R13-c22c1215c12s2R12+15c1sR1+sR1c1+6ce-sR1/c1R13,(6e)C1=6c22s2R12+6c2sR1+2ce-sR1/c2R13-c22c126c12s2R12+6c1sR1+2ce-sR1/c1R13,where λ=2μν/(1-2ν) is the Lame constant:(7a)σzz5=λK5cf1dζ+λ+2μK5cf2dζ,(7b)K5=1-ν1-2νπμ3-4ν,(7c)f1=-A5r2R2+2B5;f2=-A5ζ2R2+B5,(7d)A5=15c22s2R2+15c2sR+sRc2+6e-sR/c2R3-c22c1215c12s2R2+15c1sR+s2R2c12+4sRc1+9e-sR/c1R3,(7e)B5=3c22s2R2+3c2sR+1e-sR/c2R3-c22c123c12s2R2+3c1sR+sRc1+2e-sR/c1R3,where R={r2+(z+ζ)2}1/2 is the distance from the virtue source point (0,0,-ζ) to the point ξ(x,y,z):(8a)σzr1=μK12B1c2R12-D1,(8b)D1=6c22s2R12+6c2sR1+sR1c2+3re-sR1/c2R13-c22c126c12s2R12+6c1sR1+2re-sR1/c1R13,(9a)σzr5=2μK5cf3dζ,(9b)f3=-A5ζrR2.As shown in (6a)–(9b), it is difficult to analytically determine where σzz or σzr equals zero on the free surface. To further simplify the expressions, let the depth of the source point be zero; that is, c=0. The expressions for the stress components can be rewritten as follows:(10)σzz1=σzz2=σzz3=σzz4=σzz6=0,(11)σzz5=λK50f1dζ+λ+2μK50f2dζ,(12a)σzr1=μK1-E1;σzr2=μK2-E1,(12b)K2=3-4ν4πμ,(12c)E1=6c22s2r2+6c2sr+src2+3e-sr/c2r2-c22c126c12s2r2+6c1sr+2e-sr/c1r2,(13)σzr3=σzr4=σzz6=0,(14)σzr5=2μK50f3dζ.From (10) and (11), it can be seen that although σzz of the other five problems are zero, the value of σzz5 cannot be evaluated analytically since it involves nonintegrable integrands. As for σzr given by (12a)–(14), σzr1, σzr2, and σzr5 are nonzero. Whether they can be cancelled out when summing up cannot be determined analytically either.

With (6a)–(14) in hand, σzz and σzr on the free surface are computed numerically using MATLAB. The distributions of σzz and σzr on the free surface for different source depths are displayed in Figures 6 and 7, respectively, from which it can be seen clearly that the normal and shear stresses are nonzero. Therefore, the boundary conditions for the free surface given by (3) are not satisfied, which means that the superposition of the six full-space elastodynamic solutions does not solve the dynamic half-space problem.

σ z z on the free surface for different source depths: (a) c=0 and (b) c=λ2.

σ z r on the free surface for different source depths: (a) c=0 and (b) c=λ2.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Lamb H. On the propagation of tremors over the surface of an elastic solid Philosophical Transactions of the Royal Society of London A 1904 203 1 42 Johnson L. R. Green's function for Lamb's problem Geophysical Journal International 1974 37 1 99 131 10.1111/j.1365-246x.1974.tb02446.x Kennett B. L. N. Seismic Wave Propagation in Stratified Media 1983 London, UK Cambridge University Press Cambridge Monographs on Mechanics and Applied Mathematics MR779287 Banerjee P. K. Mamoon S. M. A fundamental solution due to a periodic point force in the interior of an elastic half-space Earthquake Engineering and Structural Dynamics 1990 19 1 91 105 10.1002/eqe.4290190109 2-s2.0-0025199617 Guzina B. B. Fata S. N. Bonnet M. On the stress-wave imaging of cavities in a semi-infinite solid International Journal of Solids and Structures 2003 40 6 1505 1523 10.1016/S0020-7683(02)00650-9 ZBL1032.74576 2-s2.0-0037369306 Maurel A. Pagneux V. Barra F. Lund F. Interaction of a surface wave with a dislocation Physical Review B—Condensed Matter and Materials Physics 2007 75 22 224112 10.1103/physrevb.75.224112 2-s2.0-34347342721 Chaillat S. Bonnet M. A new fast multipole formulation for the elastodynamic half-space Green's tensor Journal of Computational Physics 2014 258 787 808 10.1016/j.jcp.2013.11.010 2-s2.0-84888401969 Mindlin R. D. Force at a point in the interior of a semi-infinite solid Journal of Applied Physics 1936 7 5 195 202 10.1063/1.1745385 2-s2.0-36849129839 Love A. E. H. A Treatise on the Mathematical Theory of Elasticity 1906 2nd London, UK Cambridge University Press Whittaker W. L. Christiano P. Dynamic response of plates on elastic half-space Journal of the Engineering Mechanics Division ASCE 1982 108 1 133 154 2-s2.0-0020090269 Kobayashi S. Nishimura N. Green's tensors for elastic half-spaces—an application of boundary integral equation method Memoirs of the Faculty of Engineering, Kyoto University 1980 42 2 228 241 2-s2.0-0019007366