Discussion on the Time-Harmonic Elastodynamic Half-Space Green’s Function Obtained by Superposition
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.
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 [1–4]. Often expressed in Fourier-Bessel integral forms, numerical evaluation and application of these solutions are usually complex and time-consuming [5–7]. 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 , with the positive -axis pointing to the interior of the body. A periodic force is applied at point in the positive -direction, where is the circular frequency and is the imaginary unit. The displacements at point are to be found. The distance between and the real source point is given by , whereas the distance between and the image point is denoted by .
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 and , respectively; the third, fourth, and sixth problems represent a dynamic double force, a dynamic center of compression, and a dynamic doublet  at , respectively; and the fifth problem represents a line of dynamic centers of compression extending from to . To facilitate the ensuing discussion and save space, here only the -displacement component of the sixth problem is given:where is the shear modulus, is the Poisson’s ratio, and are the pressure and shear wave velocities, respectively, and 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 and into (1a), the third term in the second component of yieldsWhen the depth of the periodic point force goes to infinity, that is, , 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 -displacement at goes to infinity as . However, the displacements at a given point with finite are expected to vanish when , 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 kg/m3, shear modulus MPa, and Poisson’s ratio . The excitation frequency Hz, and thus the shear wavelength is m. The amplitude of the periodic force is N. The free surface -displacement responses at different radii as the depth of the source changes from to 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 -displacement calculated using the former solution oscillates and tends to increase when increases, while that obtained using the latter solution decays quickly with increasing .
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 N of frequency Hz is applied along the axis of symmetry in the positive -direction. The 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 -displacement responses at radii and 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.
Figure 5 compares the frequency-domain -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, . The results for between 0 and 2.0 are shown in Figures 7–10 in  and good agreement is displayed. The position of is denoted by a vertical line in Figure 5(a), from which it can be seen that when , B&M’s solution agrees reasonably well with other solutions as and deviates from other solutions as .
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 are where and are the normal and shear stresses. The equilibrium condition is given bywhere is the static point force applied at .
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 aswhere is the displacement in -direction.
The boundary conditions are examined first. The expressions of the stresses are derived. Since these expressions are lengthy, only the stress components and on the free surface, that is, , for the first and fifth problems are presented here:where is the Lame constant:where is the distance from the virtue source point to the point :As shown in (6a)–(9b), it is difficult to analytically determine where or equals zero on the free surface. To further simplify the expressions, let the depth of the source point be zero; that is, . The expressions for the stress components can be rewritten as follows:From (10) and (11), it can be seen that although of the other five problems are zero, the value of cannot be evaluated analytically since it involves nonintegrable integrands. As for given by (12a)–(14), , , and are nonzero. Whether they can be cancelled out when summing up cannot be determined analytically either.
With (6a)–(14) in hand, and on the free surface are computed numerically using MATLAB. The distributions of and 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.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
H. Lamb, “On the propagation of tremors over the surface of an elastic solid,” Philosophical Transactions of the Royal Society of London A, vol. 203, pp. 1–42, 1904.View at: Google Scholar
B. L. N. Kennett, Seismic Wave Propagation in Stratified Media, Cambridge Monographs on Mechanics and Applied Mathematics, Cambridge University Press, London, UK, 1983.View at: MathSciNet
A. E. H. Love, A Treatise on the Mathematical Theory of Elasticity, Cambridge University Press, London, UK, 2nd edition, 1906.
W. L. Whittaker and P. Christiano, “Dynamic response of plates on elastic half-space,” Journal of the Engineering Mechanics Division ASCE, vol. 108, no. 1, pp. 133–154, 1982.View at: Google Scholar
S. Kobayashi and N. Nishimura, “Green's tensors for elastic half-spaces—an application of boundary integral equation method,” Memoirs of the Faculty of Engineering, Kyoto University, vol. 42, no. 2, pp. 228–241, 1980.View at: Google Scholar