The scattering of elastic waves by microstructures in solids remains a fundamental issue both in theoretical researches and in engineering applications such as nondestructive detection and prediction of the overall dynamic properties of composites. Many works have been carried out on this subject based on the classical elastodynamics theory in literature  . The diffraction of elastic waves by a single inhomogeneity embedded in an elastic medium, as well as the dynamic stress concentration around the inhomogeneity, has been comprehensively reviewed by Pao and Mow  . These researches are mainly on the holes and inclusions of constant curvature in the whole space. Liu, et al.     used the complex variable function and the multi-polar coordinate methods to introduce the conformal mapping to analyze the scattering problems of the variable curvature holes. Shi, et al.    used the mirror method to extend the whole space problem to half space. Lee, et al.  discussed circular tunnels in half space
and obtained analytical solutions. Manoonan, et al.  used the weighted parameter method to solve the elastic half-space surface displacement of circular, elliptical and square inclusions by SH-wave. Yuan, et al.  studied the scattering and vibration of SH-wave by shallow cylindrical holes. Using a variety of different methods, many problems in life have been theoretically solved      .
With the rapid development of micro-nano-components, there is an increasing demand for understanding the mechanical behavior of small-sized materials and structures, which often differ distinctly from their macroscopic counterparts. As the volume of the object decreases, the ratio of surface area to volume increases, and the surface effect is significant, thus exhibiting mechanical behavior different from the macroscopic case. For example, geckos can walk freely on vertical walls and mosquitoes walk on the water. Based on Gurtin’s surface elasticity theory   , Sharma, et al.  studied the size dependence of the elastic field around the nano-cylindrical and nano-sphere inclusions in the whole space. Using the wave function expansion method, Wang, et al.   discussed the diffraction of plane compressional wave (P-wave) in nano-hole. Shen, et al.  discussed the influence of surface effects on the stress field around nano-inclusions. Ou, et al.    discussed the mechanical behaviors of inclusions and holes subjected to uniform loads at nano-scale.
In the present work, we discuss the scattering of plane SH-wave by an elliptical hole in elastic half-space at nano-scale. Based on the classical boundary conditions, the conformal mapping method is used to establish the boundary model and the stress concentration problem is analyzed by using the methods of wave function expansion and complex variable function.
2. Basic Equations
In the surface elasticity theory, a surface is regarded as negligibly thin membranes adhere to the bulk without slipping and it has material constants different from those of the bulk. The equilibrium and constitutive equations in the bulk of the solid read
where is the mass density of the material, t is the time. and are respectively shear modulus and Poisson’s ration, and are the stress tensor and strain tensor in the bulk material. The strain tensor is related to the displacement vector by
The equilibrium equations on the surface are expressed as 
, . (4)
where is the negative of the tangential component of the traction in the direction, and is the normal vector of the surface, denotes the curvature of the surface. The surface stresses are given for anisotropic surface by  
where is the second-rank tensor of surface strain, is the Kronecker delta, and and are surface elastic constants.
For a circular hole of the radius , according to the Equation (4) and Equation (5), we find
, . (6)
The boundary condition of stress around the circular hole can be obtained from Equation (6) as follows
with s being a dimensionless parameter reflecting the effect of the surface/interface at nano-scale. It is seen from Equation (8) that for a macroscopic hole with a big radius, , therefore the surface/interface effect can be neglected. However, when the radius of the hole shrinks to nano-scale, s becomes noticeable and the surface effect should be considered in the analysis   .
3. Theoretical Model and Analysis
Consider an elliptical nano-hole embedded in an elastic half-space, as shown in Figure 1. It is assumed that the medium is elastic, uniform, isotropic and linear. The long semi-axis of the elliptical hole is denoted by a, and the short semi-axis is denoted by b.
Figure 1. Scattering of SH-wave by an elliptical hole embedded in the half space.
If there is no hole and only incident plane wave exists in the half-space, the incident wave can be represented in the coordinate system as 
where represents the amplitude of the incident wave, i is the unit of imaginary number, is the number of the incident wave, represents the shear wave velocity of the medium, and are the mass density and shear modulus of the medium respectively.
For the steady-state problem, the time-factor is ignored; the wave displacement fulfills the following Helmholtz equation in the medium
The relationships between stress and displacement are
, . (11)
Introducing the complex plane , where the complex variable , . The Equation (10) can be expressed as follows
The incident wave passes through the half space to the boundary of the elliptical hole, which will generate scattering wave at the boundary of the hole. Using the virtual source method, the scattering wave function that satisfies the stress freeness at the horizontal boundary is expressed as 
where are the unknown coefficients and determined by the boundary conditions, is the first kind Hankel function of the n-th order.
Introducing conformal mapping function from elliptical hole to unit circular hole ( ), the boundary of the elliptical hole is mapped to the unit circle boundary as shown in Figure 2.
here , , , .
Substituting the conformal mapping function into the Equation (9), Equation (11), Equation (12) and Equation (13), we can obtain
Figure 2. Image plane after conformal transformation.
In order to improve the convergence speed and calculation accuracy, introducing Graf addition formula 
where is the first kind Bessel function of the n-th order. Substituting Equation (20) into Equation (18), we can obtain
4. Problem Solving
Substituting Equation (14) into Equation (16) and Equation (17), we can obtain stress expressions of the incident wave as
and substituting Equation (22) into Equation (16) and Equation (17), stress expressions of the scattering wave can be expressed as
where and are the radial stress, and are tangential stress, Re is the real part and Im is the imaginary part of a complex variable.
The total stress and total displacement are expressed as
The boundary conditions of the boundary of the half space and the unit circle hole respectively are
As shown in Equation (13), the scattering wave will be determined if the unknown coefficients are solved. Substituting Equation (23), Equation (24), Equation (25) and Equation (26) into Equation (32), we find
where the terms and are presented in Appendix A. It is evident that Equation (33) is a series of equations related to though there are still unknown coefficients. Multiplying both sides of Equation (33) with and integrating between the interval , we can obtain
where , .
A set of infinite algebraic equations with respect to the unknown constants will be obtained from Equation (34). It should be noted that it is mathematically convenient to truncate the infinite matrix in Equation (34) with during numerical calculation.
5. Problem Results and Analysis
Elastic-wave-induced dynamic stress concentration factor (DSCF) is an important parameter for engineering application. The SH-wave-induced DSCF is calculated as
where is the maximum amplitude of the incident stress.
To validate the present approach, a comparison with Liu’s conclusions is shown in Figure 3. It can be seen that the present results on DSCF around a circular hole in the half-space with a surface parameter are accordant with Liu’s conclusions  . Figure 4 shows the distributions of DSCF about the angle with different values of the surface parameter s. With the increase of s, the DSCF continuously decreases. It is seen that surface energy significantly affects DSCF around the hole.
We discuss the influences of surface effects, incident wave frequencies and the ration of the major and minor axis of the ellipse on the dynamic stress concentration factor near an elliptical hole now.
Figure 5 shows the distributions of DSCF about the angle , when , and . The results indicate that DSCF is symmetric about the y-axis. With the increasing of s, the DSCF continuously decreases around or . When , the DSCF in the ranges of is greater than the cases of and . When , the DSCF is zero.
Next, we discuss the distributions of DSCF around the elliptical hole at different frequencies, when . It is seen that when the surface/interface effect is taken into account, the dynamic stress depends not only on s but also on the frequencies. If and keep as a constant, then the incident SH-wave represents a static biaxial loading  .
For a high-frequency incident wave with . Figure 6 displays the distributions of DSCF around the elliptical hole for different interface parameter s. In this case, the DSCF is symmetric about the y-axis. With increasing s, the DSCF declines continuously in the angular range of . For a low-frequency incident wave with , the distributions of DSCF around the elliptical hole for different values s are plotted in Figure 7. It can be seen that DSCF decreases in the angular range of increasing s but slowly increases elsewhere. Comparing the high frequency and low frequency, we find that DSCF shows a decreasing trend as the frequency increases.
Figure 3. The distribution of DSCF near a circular hole for different wave number with .
Figure 4. Effect of interface parameter s on DSCF near a circular hole for , .
Figure 5. Effect of interface parameter s on DSCF near an elliptical hole for , .
Figure 6. Effect of interface parameter s on DSCF near an elliptical hole for , .
Figure 7. Effect of interface parameter s on DSCF near an elliptical hole for , .
In this paper, based on the theory of surface/interface elasticity, we focused on the observation that dynamic stress concentration factor of the elliptical hole. We used conformal mapping to transform an elliptical hole into a unit circular hole. Therefore, it simplifies the process of processing elliptical boundary conditions. The numerical solutions of the elastic fields induced by SH-wave near an elliptical nano-hole are obtained. The influence of surface/interface effect on the DSCF is discussed. We get the conclusion that DSCF around the hole is not only related to interface effects, but also related to frequencies.
The supports from the National Natural Science Foundation (Grant No. 11362009 and No. 11862014) are acknowledged.