Simulations of Rayleigh ’ s Wave on Curved Surface

Impulsive line load in a half-space (Lamb’s problem) can be solved with a closed form solution. This solution is helpful for understanding the phenomenon of Rayleigh’s waves. In this article, we use a boundary element method to simulate the solution of an elastic solid with a curved free surface under impact loading. This problem is considered difficult for numerical methods. Lamb’s problem is calculated first to verify the method. Then the method is applied on the problems with different surface curvatures. The method simulates the phenomenon of Rayleigh’s wave propagating on a curved surface very well. The results are shown in figures.


Introduction
The phenomenon of surface wave is interesting and important for many engineers and scientists.Impulsive line load in a half-space can be solved by analytic methods [1].Herewith scientists and engineers may get understanding with Rayleigh's wave.When the free surface is not flat the phenomenon of surface wave should be similar, but some detail can be different.On the case of curved surface, analytic solutions are no longer available.Therefore numerical methods become an alternative way to understand the phenomena.These problems are considered difficult for numerical methods.To the best of our knowledge, there are no related reports on this problem.
Simulating transient wave in elastodynamics needs mass computing time and huge memories.Boundary element methods are efficient for simulating elastodynamics [2,3].Recently personal computers (PC) become much more powerful and are equipped with more rams then ever.Practical problems can be simulated with a PC precisely.
In this article, we use a boundary element method to solve two dimensional elastodynamic problems with curved surfaces.The curved boundary is assumed to be an arc.The loading is an impulse.The numerical method is implemented with fortran programs and a PC.In Section 2, the mathematical problem is described.In the following section, we formulate the numerical method briefly.The results are shown in Section 4.

The Problem
Because our goal is to simulate the Rayleigh wave on a smooth surface near the loading point, the boundary of the 2D domain is modeled as a arc with curvature . The loading is a single line impact at the surface.
The schematic diagram is shown in Figure 1.
Boundary value problems with zero initial conditions and absent body force for 2D elastodynamics in plane strain condition are considered.The material is homogeneous, isotropic and linear elastic.
where  and  are the Lame constants an d  is the mass density of the elastic material.

,
The boundary conditions are where  is the stress tensor,   is the outward normal vector on ,  s is the coordinate on the boundary, is the Dirac delta function, and   s x presents the map from coordinate s to 2D spatial coordinates . The time-varied displacements on the curved surface are Rayleigh's waves.

The Method
We use a boundary element method to calculate the surface displacements.The boundary is approximated as a polygon.


There are two families of particular solutions of Navier's equation, Using Hooke's law, we have stress bases,

 
, , , , , Note that are vector fields and The bases have been derived in a close form [4].The the coefficients, and , will determined with the boundary condition (2).
x ik a y ik a Then, the approximated stress field is When a collocation method apply on Equation (7), the coefficients x n are as many as .It is difficult to calculate a precise elastodynamic solution on a personal computer with this formulation.Therefore we take the advantage of the symmetry of the boundary.
and cos sin sin cos sin cos Substituting Equations ( 8) and ( 9) into (7), we have Then, we decomposite the traction into normal and tangent directions. .Then apply the semicollocation method [4] to Equations (15) and (16).Rearranging the equations for and , we have the stepping equations. where and 1 if 0 and 0 0 e l s e 0.
The   where H is the Heaviside step function.
In formulae ( 17) and ( 18), the coefficients , ik j take the form of   ,00 . The usage of computer memory is enormously deduced to 2Nm.
Then the time-stepping technique is applied on Equations (17) and (18) to solve and on th step.
After solving the coefficients, and , the numerical displacement may be calculated by and the numerical interior stress by

The Results
In order to verify this method, we calculate the problem with very large first.For , the problem becomes Lamb's problem in which the exact surface displacements are available.The elastic half space The solutions for problems of a line impulse load on an elastic half plane were derived by Lamb.A modern treatment with integral transform technique was given by De Hoop, but the results were in complex function form.Nevertheless explicit form for surface displacements are available.The analytic solution for surface displacement can be found in page.614-626 of [1].
The schematic diagram of the geometry of the problem is shown in Figure 1.
In this example Poisson's ratio 0.25 . Therefore the shear and Rayleigh's wave speeds are and respectively, where is the dilatation wave speed. .Even though this problem is considered difficult for numerical methods, our results are precise.This example shows the method is applicable for impact problems.When 5 R  , the results are shown in Figure 3.In this case tangential displacement is almost the same with Lamb's problem, but the normal displacement is changed much.The normal displacement has a strong and short peak at the shear wave front and the head wave is enlarged.The head wave is the wave beyond the shear wave front [5].The dashed line is the displacement for R   .Figure 4 shows the displacements for 3 R  .The phenomenon is similar to the case of .5 R 

Conclusion
We use the boundary element method and write a fortran program running on a personal computer.In order to simulate the phenomenon of Rayleigh's wave propagating on a curved surface, the free surface is assumed to be a constant, 1 R .The method is verified by the problem of 1000 R  . Two examples, R = 5 and R = 3, calculated to demonstrate the phenomenon.On the shear wave front,   the wave is very different from that of Lamb's problem.

Figure 1 .
Figure 1.The schematic diagram of the geometry of the problem.

Figure 2
Figure 2 shows the displacements with vertical loading when 1000 R  .The dashed line shows exact displacement when R   .There is a Dirac delta function at the Rayleigh wave front in the exact solution.The vertical line indicate the arriving time of Rayleigh's wave.For this figure, 1 600 s  .Even though this pro-

Figure 2 .
Figure 2. The displacements with vertical loading when R = 1000.The dashed line shows exact displacement when R = ∞.The vertical red lines indicate the arrival time of Rayleigh's wave.There is a Dirac delta function at the Rayleigh wave front in the exact solution.In this figure, 1 600 s   is used.

Figure 3 .
Figure 3.The displacements when R = 5.The dashed line shows exact displacement when R = ∞.The vertical red lines indicate the arrival time of Rayleigh's wave.In this figure, 1 400 s   is used.

Figure 4 .
Figure 4.The displacements when R = 3.The red line shows exact displacement when R = ∞.The vertical dashed lines indicate the arrival time of Rayleigh's wave.In this figure, 1 400 s   is used.