Abstract
By combining the wave function expansion method with the auxiliary function, a closed series solution for the dynamic plane response is proposed considering horizontal circular heterogeneous topography under SH wave incidence. The displacement and stress residual is found to be small along the entire boundary, including the corner points; further, the solution can be used for the incidence of a highfrequency wave. In this study, the influence of the shear modulus ratio of a heterogeneous hill and a lower medium on the displacement amplitude spectrum and surface displacement is determined. The results show that the hardness of the medium in hill topography significantly influences outofplane surface motion. The surface displacement of the soft hill is significantly larger than that of the homogeneous hill, and the motion effects are enhanced considerably. In contrast, the hard hill weakens surface motion.
Highlights
 The number of series iteration is more with the increase of frequency of incident wave to ensure the calculation accuracy.
 The incident angle, dimensionless frequency, and shear modulus ratio of the heterogeneous hill and lower medium have a great influence on surface displacement.
 Because of the high calculation accuracy, it is possible to study the characteristics of surface.
1. Introduction
The site response problem is an important wave problem in seismic engineering. Several seismic damage data show that local terrains have a great impact on the scattering and diffraction of seismic waves. In recent years, largescale infrastructure development has been carried out in many cities globally, where raised terrains are common in engineering sites. Engineering structures near mountainous areas should consider the amplification effect of the terrain. Therefore, exploring the scattering law of seismic waves by heterogeneous protrusion has high theoretical significance and clear engineering reference value.
Various field measurements have shown that variations in topographic configuration influence the intensity of ground motion [15]. Numerous mathematical models provide insight into the behavior of ground shaking due to geologic and topographic variations. Most closedform analytical solutions are related to the irregularities of subsurface, such as horizontally stratified surface layers [6], semicircular canyons or valleys [7, 8], semielliptical canyons or valleys [9, 10], and circular underground cavities or tunnels [11, 12].
However, for the surface convex or hill topography, the scattering of elastic waves in hills is also a challenging problem due to mixed boundary value. According to earthquake investigators, the amplification effect of actual hill topography is much larger than the theoretical results. Some scholars speculate that it is caused by the difference in properties between hill topography and the underlying medium [13]. Based on the method of partition and auxiliary function, Yuan et al. [1416] proposed the auxiliary function and the external region Graf addition formula and obtained the closedlevel solution for the SHwave scattering by semicircular and arbitrary circular hill topography. Based on the above theory, Vincent [17, 18] used an improved wave function expansion method to resolve the analytical solution of the SH wave scattering of the semicircular convex terrain. Li and Yuan [19] used the auxiliary function method and wave function expansion to obtain an outofplane dynamic response of heterogeneous protrusions in halfspace. Tsaur et al. [20] solved the scattering of plane SH waves by arcarcconvex terrain and obtained the corresponding series solution. Tsaur [21] also studied the scattering of SH wave by semielliptical convex terrain and calculated the steadystate response. Amornwongpaibun et al. [22] proposed an analytical solution for the closedwave function of twodimensional scattering and diffraction of SH waves in a semielliptical shallow hill with a concentric elliptical tunnel in a semielastic halfspace. Liu et al. [23] investigated the dynamic interaction between a lined tunnel and a hill under plane SV waves using the indirect boundary element method.
In this study, an improved wave function method is used to solve the outofplane dynamic response of a halfspace heterogeneous hill. The wave function is expanded into a cosine function by Fourier series [11]. The solution is reduced to a set of infinite algebraic equations; a numerical solution can then be obtained from the truncation of the infinite equations. The calculation results show that the displacement and stress residual is minimal and can be ignored along the entire boundary, including the corner points. The calculation accuracy is quite high and the method can offer effective solutions for highfrequency waves. Based on the improved wave function method, the effect of key factors, such as the angle and frequency of the incident wave and material property of hill on the dynamic response, is investigated, and the conclusions will help with the seismic design of underground mountain structures.
2. Physical model
The crosssection of the model is shown in Fig. 1. It is identical to that considered by Yuan and Men [14]. It represents a homogeneous, isotropic, and fully elastic halfspace with a semicircular hill of radius $a$. ${\mu}_{0}$ and ${c}_{0}$ denote the shear modulus and wave velocity of the halfspace, respectively. ${\mu}_{1}$ and ${c}_{1}$ denote the shear modulus and wave velocity of the circular hill, respectively. The flat surface boundary is recorded as $\mathrm{\Gamma}$, and the height and halfwidth of the hill topography are $h$ and $b$, respectively. This can be written as $a=b=h$.
Fig. 1Physical model of a heterogeneous hill topography
Fig. 2Regional division of the model
As shown in Fig. 2, the physical model is divided into two regions. Region $\mathrm{\Theta}$ is the circular region including the circular hill, in which the upper and lower boundaries are denoted as $L$ and $\overline{L}$, respectively. Region $D$ forms the rest of the halfspace. The common boundary of regions $\mathrm{\Theta}$ and $D$ is $\overline{L}$, and the flat surface beyond the hill is denoted by $\mathrm{\Gamma}$. By dividing the original complex model into two completely independent small circular or flat boundary regions, the wave fields in these regions can be represented by the FourierBessel wave function. When combined with the continuity of displacement and stress between the two regions, the boundary conditions of the original model can be derived. The remaining boundary conditions for the waves are as follows:
(1) Stressfree conditions at hill boundary $L$ in region $\mathrm{\Theta}$.
(2) Stressfree conditions at flat surface $\mathrm{\Gamma}$ in the region $D$.
(3) Continuity of displacement and stress at $\overline{L}$.
3. Numerical model
3.1. Scattered waves: wave function series expansion
Assuming that the displacement of SH wave donated by ${u}^{\left(i\right)}$ is excited on the halfspace medium, it can be represented in Cartesian coordinate system $xoy$ by:
where $i=\sqrt{1}$, ${c}_{x}=c/\mathrm{c}\mathrm{o}\mathrm{s}\alpha $ and ${c}_{y}=c/\mathrm{s}\mathrm{i}\mathrm{n}\alpha $ are phase velocity of the incident wave in the $x$ and $y$ directions, respectively, $c$, $\omega $, $\alpha $, and ${u}_{0}$ denote the velocity, circular frequency, incidence angle, and amplitude of the incident SH wave.
When the incident wave ${u}^{\left(i\right)}$ goes through the hill, the total displacement field $u$ should satisfy the wave differential equation [24]:
The stress field should appropriate the stressfree boundary conditions, which is given by:
where:
where, $\mu $ is the shear modulus of the medium of half space soil.
The displacement field of the concave halfspace $D$ can be expressed as:
where ${u}^{\left(f\right)}$ denotes the displacement of free field excited by incident SH wave, ${u}^{\left(s\right)}$ denotes the scattering potential field due to the interaction between the fictitious hill boundary and half space.
The free field displacement ${u}^{\left(f\right)}$ consists of the incident wave displacement ${u}^{\left(i\right)}$ and the reflected wave displacement ${u}^{\left(r\right)}$, which is excited by ${u}^{\left(i\right)}$ at the flat surface of halfspace:
The equation of incident wave ${u}^{\left(i\right)}$ is as shown in Eq. (1). The equation of reflected wave ${u}^{\left(r\right)}$ can be represented as:
In Eq. (1) and Eq. (9), $\mathrm{e}\mathrm{x}\mathrm{p}\left(i\omega t\right)$ is the harmonic factor and can be omitted in subsequent wave expressions. $k=\omega /c$ denotes the wave number. Substituting both ${c}_{x}=c/\mathrm{c}\mathrm{o}\mathrm{s}\alpha $ and ${c}_{y}=c/\mathrm{s}\mathrm{i}\mathrm{n}\alpha $ into Eq. (1) and Eq. (9) gives:
According to the coordinate transformation relationship $x=r\u2022\mathrm{c}\mathrm{o}\mathrm{s}\theta $, $y=r\u2022\mathrm{s}\mathrm{i}\mathrm{n}\theta $. Eq. (10) will be thus be changed to Eq. (11):
Substituting Eq. (11) into Eq. (8) gives:
where ${J}_{n}(\u2022)$ is Bessel function of the first kind with order n, ${a}_{o,n}=2{\mathrm{\epsilon}}_{n}{i}^{n}\mathrm{c}\mathrm{o}\mathrm{s}n\alpha $, ${a}_{o,n}$ are the coefficients of freefield waves, ${\epsilon}_{0}=\text{1}$, ${\epsilon}_{n}=\text{2}$ for $n=$1, 2, 3….
Substituting Eq. (12) into Eq. (6):
where ${\sigma}_{0}={\mu}_{0}{k}_{0}{u}_{0}$, ${\sigma}_{0}$ is amplitude of the incident wave stress. It is derived from the expansion theorem of Bessel functions [1] and used extensively [11, 2529].
Solving the total wave Eq. (2) according to the separation variable method, the expression of the scattering displacement field generated by the semicircular boundary satisfy Eq. (2) and the boundary condition Eq. (3). It can be taken as:
where ${A}_{n}$ is an unknown coefficient of the new waves to be determined. ${H}_{n}^{\left(1\right)}(\u2022)$ is a Hankel function of the first kind with the order $n$.
Substituting Eq. (14) into Eq. (5):
The equation of the stressfree boundary condition Eq. (3) is always automatically satisfied on the boundary $\mathrm{\Gamma}$ ($\theta =0,\pi $). By substituting Eq. (15) into Eq. (6), we get:
The displacement field of the cohesive wave generated in region $\mathrm{\Theta}$ owing to the fictitious boundary $\overline{L}$ and the arcshaped hill surface $L$ can be expressed as:
where ${B}_{n}$ and ${C}_{n}$ are unknown coefficients of new waves to be determined.
The FourierBessel series expansions of all incident waves and scattered waves in each coordinate system are thus obtained.
3.2. A new analytical method after half range expansion of cosine function
To divide the boundary into two parts processed separately, the wave function ${u}^{\left(c\right)}$ must be expanded into an orthogonal trigonometric form in the upper and lower semiregions. The sine and the cosine functions are found to be orthogonal over the entire $2\pi $ range but not between the upper and lower $\pi $ ranges. In this case, problems can be solved by using the orthogonal cosine function in the upper and lower half ranges. This method has solved the analytical solution of the scattering of SH waves by semicircular hill topography [10].
3.2.1. Referencing half range expansion of cosine function
The cosine function is orthogonal in the half range $[0,\pi ]$ and $[\pi ,0]$. Any function $f\left(\theta \right)$ can be expanded by using the cosine series function in the range as:
where:
Using Eq. (19), the sine function can be expressed at half ranges of $[0,\pi ]$ and $[\pi ,0]$ as:
where, ${\epsilon}_{0}=$ 1, ${\epsilon}_{n}=$ 2, $n>$ 0, $m=$1, 2, 3,….
3.2.2. Reformulation of the relevant displacement field
The cohesive waves generated in the circular region $\mathrm{\Theta}$ owing to the fictitious boundary $\overline{L}$ and the arcshaped hill surface $L$ are represented as orthogonal cosine functions in the half ranges $[0,\pi ]$ and $[\pi ,0]$:
Substituting Eq. (21), Eq. (6) takes the form:
3.3. Transformation of boundary conditions
The FourierBessel expressions of various scattered waves in all regions are obtained. The boundary conditions are introduced to establish a system of equations for each group of undetermined coefficients.
(1) the continuity condition of displacement and stress at the boundary $\overline{L}$.
(2) stressfree boundary condition on circular hill surfaces:
Substituting the corresponding wave function into Eqs. (2325), we obtain:
3.4. Referencing equation elimination to determine the unknown coefficients
Eqs. (2628) contain three series of equations for solving three sets of unknowns ${A}_{n}$, ${B}_{n}$, and ${C}_{n}$. There are a variety of methods to solve the unknowns, and one of them is listed below.
From Eqs. (2728), we obtain:
Deformation of Eq. (29):
Substituting Eq. (29) into Eq. (26):
From Eq. (28) and Eq. (30), the following equations can be obtained:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\left[\frac{{J}_{n}\left({k}_{1}a\right)}{{J}_{n}^{\mathrm{\text{'}}}\left({k}_{1}a\right)}\frac{2\frac{{\mu}_{1}^{}}{{\mu}_{0}}{H}_{n}^{\left(1\right)}\left({k}_{0}a\right)}{{{H}_{n}^{\left(1\right)}}^{\mathrm{\text{'}}}\left({k}_{0}a\right)}\right]\frac{{\epsilon}_{n}}{\pi}\sum _{\begin{array}{l}m=1\\ n+m\in odd\end{array}}^{\mathrm{\infty}}{J}_{m}^{\mathrm{\text{'}}}\left({k}_{1}a\right)\frac{2m}{{m}^{2}{n}^{2}}{C}_{m}$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\frac{{\epsilon}_{n}}{\pi}\sum _{\begin{array}{l}m=1\\ n+m\in odd\end{array}}^{\mathrm{\infty}}{J}_{m}\left({k}_{1}a\right)\frac{2m}{{m}^{2}{n}^{2}}{C}_{m}.$
The solution is obtained by solving a set of infinite algebraic equations. The unknown coefficients are determined by equation truncation. The unknowns ${C}_{n}$, ${B}_{n}$,_{}and ${A}_{n}$ can be obtained from Eqs. (31), (28), and (29), respectively.
4. Numerical analysis
It is necessary to check convergence by calculating the displacement and stress residuals on every boundary. The residuals include:
(1) The continuity of displacement and stress at $\overline{L}$ between region $D$ and region $\mathrm{\Theta}$.
(2) The stressfree conditions at the hill boundary $L$ have the dimensionless displacement residual and the residual stress on the common boundary $\overline{L}$, respectively.
(3) Stress residual at the hill boundary $L$:
The stress residual of the free surface is strictly satisfied when the wave function is set and will not be discussed here.
The nature of the SH wave is determined by frequency and wave velocity. Wave velocity depends on the physical properties of halfspace, from which the wavelength can be determined. At the same time, the influence of waves on surface displacement has a great relationship with the geometrical characteristics of topography and the wavelength. To consider this influence, the dimensionless frequency $\eta $ is introduced and defined as the ratio of the hill diameter 2$a$ to the wavelength of the incident wave $\lambda $.
The formula is:
The characteristic length of topography takes the halfwidth of hill topography. The numerical results depend on the dimensionless frequency $\eta $, shear modulus ratio of the heterogeneous hill to adjacent medium ${\mu}_{1}/{\mu}_{0}$ and incidence angle of the wave $\alpha $.
The data analysis is divided into four sections from 4.14.4. They are convergence and accuracy analysis; a set of typical examples; the influence of shear modulus ratio of heterogeneous hill and lower medium on surface displacement; and the influence of shear modulus ratio of heterogeneous hill and lower medium on displacement amplitude spectrum.
4.1. Convergence and accuracy analysis
Fig. 3 shows the variation of residual of hill surface displacement with the iteration number at the dimensionless frequencies $\eta =$1.0, 3.0, 5.0, 10.0, for ${\mu}_{1}/{\mu}_{0}=$1/4, $\alpha =$90°, in which NUM is the iteration number of the series for each calculation, e is the root mean square of the iterative residuals. NUM reaches enough small value when it equals to 20, 35, 40, and 60, respectively. It can be seen that the number of series iteration is more with the increase of frequency of incident wave to ensure the calculation accuracy.
Fig. 3Series convergence analysis graph
a) Iterative displacement residual graph on surface of hill when $\eta =$1.0
b) Iterative displacement residual graph on surface of hill when $\eta =$3.0
c) Iterative displacement residual graph on surface of the hill when $\eta =$ 5.0
d) Iterative displacement residual graph on surface of the hill when $\eta =$10.0
Figs. 47 show the displacement and stress residual amplitudes at the boundaries under four different dimensionless frequencies when ${\mu}_{1}/{\mu}_{0}=$ 1/4. It can be seen that the residual on the entire boundary is so small that it is negligible.
Fig. 4Displacement and stress residual amplitude on various boundaries for μ1/μ0= 1/4, η= 1.0, α= 90°: a) displacement residual amplitude on L, b) stress residual amplitude on L¯, c) stress residual amplitude on L
a)
b)
c)
Fig. 5Displacement and stress residual amplitude on various boundaries μ1/μ0= 1/4, η= 3.0, α= 90°: a) displacement residual amplitude on L, b) stress residual amplitude on L¯, c) stress residual amplitude on L
a)
b)
c)
4.2. Typical examples
The amplitude of dimensionless displacement on the free surface is defined as:
where $Re\left(u\right)$ and $Im\left(u\right)$ are the real and imaginary part of surface displacement $u$, respectively.
Fig. 6Displacement and stress residual amplitude on various boundaries μ1/μ0= 1/4, η= 5.0, α= 90°: a) displacement residual amplitude on L, b) stress residual amplitude on L¯, c) stress residual amplitude on L
a)
b)
c)
Fig. 7Displacement and stress residual amplitude on various boundaries for μ1/μ0= 1/4, η= 10.0, α= 90°: a) displacement residual amplitude on L, b) stress residual amplitude on L¯, c) stress residual amplitude on L
a)
b)
c)
For convenience to compare the results of various numerical and analytical methods with the results given in this paper. Figs. 815 show eight typical examples. It shows that the incident angle, dimensionless frequency, and shear modulus ratio of the heterogeneous hill and lower medium have a great influence on surface displacement.
The calculation parameters include (1) shear modulus ratio of heterogeneous hill and lower medium ${\mu}_{1}/{\mu}_{0}=$ 1/4, 4); (2) dimensionless frequency: $\eta =$ 1.0, 3.0, 5.0, 10.0 (3) four angles of incidence: $\alpha =$ 0°, 30°, 60°, 90°. In this section, the abscissa of the surface displacement map is defined as $x/a$, where a is the radius of the semicircular hill.
It can be seen from Figs. 815 that as the frequency of incident wave increases, the surface displacement changes more severely, especially for the hill topography. When the incidence angle becomes small, the surface displacement on the left side of hill is more severe than that on the right side, which explains the protection effect of hill topography being similar with the actual vibration isolation trench. In addition, for the case of normal incidence of seismic waves, especially when the wavelength of the incident wave is close to the dimension of hill topography, the surface displacement is minimal near the corner of the hill, which is consistent with actual earthquake damage.
Fig. 8Surface displacement amplitude when the shear modulus ratio of the heterogeneous hill and the lower medium is (μ1/μ0= 1/4, η= 1.0)
a)$\alpha =$0°
b)$\alpha =$30°
c)$\alpha =$60°
d)$\alpha =$90°
Fig. 9Surface displacement amplitude when the shear modulus ratio of the heterogeneous hill and the lower medium is (μ1/μ0= 1/4, η= 3.0)
a)$\alpha =$0°
b)$\alpha =$30°
c)$\alpha =$60°
d)$\alpha =$90°
Fig. 10Surface displacement amplitude when the shear modulus ratio of the heterogeneous hill and the lower medium is (μ1/μ0= 1/4, η= 5.0)
a)$\alpha =$0°
b)$\alpha =$30°
c)$\alpha =$60°
d)$\alpha =$90°
Fig. 11Surface displacement amplitude when the shear modulus ratio of the heterogeneous hill and the lower medium is (μ1/μ0= 1/4, η= 10.0)
a)$\alpha =$0°
b)$\alpha =$30°
c)$\alpha =$60°
d)$\alpha =$90°
Fig. 12Surface displacement amplitude when the shear modulus ratio of the heterogeneous hill and the lower medium is (μ1/μ0= 4, η= 1.0)
a)$\alpha =$0°
b)$\alpha =$30°
c)$\alpha =$60°
d)$\alpha =$90°
Fig. 13Surface displacement amplitude when the shear modulus ratio of the heterogeneous hill and the lower medium is (μ1/μ0= 4, η= 3.0)
a)$\alpha =$0°
b)$\alpha =$30°
c)$\alpha =$60°
d)$\alpha =$90°
Fig. 14Surface displacement amplitude when the shear modulus ratio of the heterogeneous hill and the lower medium is (μ1/μ0= 4, η= 5.0)
a)$\alpha =$0°
b)$\alpha =$30°
c)$\alpha =$60°
d)$\alpha =$90°
Fig. 15Surface displacement amplitude when the shear modulus ratio of the heterogeneous hill and the lower medium is (μ1/μ0= 4, η= 10.0)
a)$\alpha =$0°
b)$\alpha =$30°
c)$\alpha =$60°
d)$\alpha =$90°
4.3. The influence of the shear modulus ratio of the heterogeneous hill to adjacent medium on surface displacement
According to the results from Sections 4.1 and 4.2, the shear modulus ratio of the heterogeneous hill and the lower medium has a significant effect on the surface displacement. Therefore, it is necessary to compare the shear modulus ratios of the heterogeneous hill with the lower medium at different time ranges. In this section, the abscissa of the surface displacement map is defined as $x/b$, where $b$ is the half width of the hill topography.
The calculation parameters then include (1) shear modulus ratio of heterogeneous hill and lower medium on semicircular hill ($h=a=b=$10), ${\mu}_{1}/{\mu}_{0}=$ 1/4, 1, 4; (2) dimensionless frequency $\eta =$1.0, 3.0, 5.0, 10.0; (3) two angles of incidence $\alpha =$ 0°, 90°.
Fig. 16Surface displacement under influence of different shear modulus ratios of heterogeneous hill and the lower medium (η= 1.0)
a)$\alpha =$0°
b)$\alpha =$90°
Fig. 17Surface displacement under influence of different shear modulus ratios of heterogeneous hill and the lower medium (η= 3.0)
a)$\alpha =$0°
b)$\alpha =$90°
Figs. 1619 show the amplitude of surface displacement near the hill in cases with the hill topography under three different hardness (${\mu}_{1}/{\mu}_{0}=$1/4, 1, 4)) with the different frequency ($\eta =$ 1.0, 3.0, 5.0, 10.0). In these figures, (a) and (b) indicate the effects of two different angles ($\alpha =$ 0°, 90°) of incident waves
It can be noted that the displacement amplitude of the free surface is constant at 2. Figs. 1619 indicate that the medium hardness of hill topography has a great influence on the surface amplification. Regardless of the incidence angle of the wave, compared with the homogeneous hill, the surface displacement of the soft hill increases significantly, and that of hard hill reduces. This difference between surface displacement is quite large, for vertical incidence of wave, when the medium hardness of hill changes from hard (${\mu}_{1}/{\mu}_{0}=\text{4}$) to soft (${\mu}_{1}/{\mu}_{0}=$ 1/4). For $\eta =$ 1.0 the maximum amplitude of hill top changes from 1 (0.5 times surface displacement of free field without hill) to 8 (4 times surface displacement of free field without hill); the maximum amplitude of hill top changes from 0.5 to 10 for $\eta =$ 3.0 and from 0.25 to 11 for $\eta =$ 10.0.
In addition, the location of the maximum displacement of the hill surface is mainly affected by the incidence angle of waves. The surface displacement on the hill midpoint is prone to amplification for the vertical incidence of the wave. For the grazing incidence of the wave, the amplification of hill surface displacements appears on the back top of the hill, and the location of the maximum displacement of hill surface will move toward the right hill surface with an increase in incidence angle. The motion of the horizontal surface at the hill back will reduce due to the protective effect of the hill, the surface motion fluctuation effect will enhance with increasing frequency of incident wave, and the displacement amplitude of the hill surface will change very sharply for the highfrequency wave.
Fig. 18Surface displacement under influence of different shear modulus ratios of heterogeneous hill and the lower medium (η= 5.0)
a)$\alpha =$0°
b)$\alpha =$90°
Fig. 19Surface displacement under influence of different shear modulus ratios of heterogeneous hill and the lower medium (η= 10.0)
a)$\alpha =$0°
b)$\alpha =$90°
4.4. Referencing the influences of shear modulus ratio of the heterogeneous hill and lower medium on the displacement amplitude spectrum
Figs. 2022 show the displacement amplitude spectrum at three positions of the hill surface (top, left corner, and right corner), in which ${\mu}_{1}/{\mu}_{0}$ is the ratio of the shear modulus of the heterogeneous hill to the adjacent medium, $\alpha $ is the incidence angle of the wave, and $\eta $ is the dimensionless frequency.
Fig. 20Displacement amplitude spectrum at the central point of the hill surface
a)$\alpha =$0°
b)$\alpha =$45°
c)$\alpha =$90°
Fig. 21Displacement amplitude spectrum at the left corner of the hill surface
a)$\alpha =$0°
b)$\alpha =$45°
c)$\alpha =$90°
Fig. 22Displacement amplitude spectrum at the right corner of the hill surface
a)$\alpha =$0°
b)$\alpha =$45°
c)$\alpha =$90°
The data from Figs. 2022 show the surface displacement of the soft hill surface differs from that of the homogeneous and hard hill as follows:
(1) On the ground fixed observation point, in the case of homogeneous and hard conditions, the displacement amplitude of the surface changes with the frequency as a smooth curve, while the amplitude of the surface displacement of the soft hill changes with frequency as a nonsmooth curve. As the hill topography gets softer, the fluctuation is more severe, implying that softer hill is more sensitive to the incident wave frequency.
(2) The soft hill area has a strong amplification effect on the upper surface motion, while the hard hill has a reducing effect on the ground motion. In addition, as the frequency increases, the amplitude of the surface displacement increases, with the maximum value tending to appear near $\eta =$ 4.
As is shown in Fig. 20, in the amplitude spectrum of the surface displacement of the hill top point under the normal wave ($\alpha =$90°), the maximum displacement amplitude of the top point of the homogeneous hill is 3.1, when the average displacement amplitude in the range of $0\le \eta \le \text{10}$ (the arithmetic mean of ${U}_{d}$ in the range of $0\le \eta \le \text{10}$) is 2.1. The maximum displacement amplitude of the soft hill top is 40 when ${\mu}_{1}/{\mu}_{0}=\text{1/4}$ and the average displacement amplitude in the range of $0\le \eta \le \text{10}$ is 7.8. Compared with the homogeneous conditions, these values are 13 and 3.7 times higher, respectively.
5. Conclusions
An improved wave function method is used innovatively for an incident SH wave in a halfspace heterogeneous hill. The results show that the displacement and stress residual is so small along the entire boundary, including the corner points, and can be ignored. The solution is effective for highfrequency waves because of the high calculation accuracy. The frequency is calculated to be $\eta =$10 in this paper, it is possible to study the characteristics of surface displacement at high frequencies. The medium hardness of hill topography has a great influence on the surface motion. Compared with the homogeneous hill, the surface displacement of the soft hill significantly increases, and the hard hill weakens the surface motion. It is suggested that seismic measures should be strengthened for engineering applications near the soft hill.
Acknowledgements
This work was supported by the Tianjin Research Program of Application Foundation Advanced Technology (17JCYBJC21700), the National Natural Science Foundation of China under grants (51878108), the key projects of Tianjin science and technology support program (17YFZCSF01140), and Tianjin Municipal Science and Technology Bureau (19PTZWHZ 00080).
References

Abramowitz M., Stegun I. A. Handbook of Mathematical Functions, with Formulas, Graphs, and Mathematical Tables. Dover, New York, 1972.

Gutenberg B. Effect of ground on earthquake motion. Bulletin of the Seismological Society of America, Vol. 47, Issue 3, 1957, p. 221250.

Kanai K., Suzuki T. Relation between the property of building vibration and the nature of ground (observation of earthquake motion at actual buildings). Earthquake Research Institute, Vol. 31, Issue 1, 1953, p. 305316.

Kanai K., Suzuki T., Yoshizawa S. Relation between the property of building vibration and the nature of ground (observation of earthquake motion at actual buildings) III. Earthquake Research Institute, Vol. 34, Issue 1, 1955, p. 6186.

Kanai K., Tanaka T., Yoshizawa S. Comparative studies of earthquake motions on the ground and underground (multiple reflection problem). Bulletin of the Earthquake Research Institute, Vol. 37, 1958, p. 5387.

Tsai N. C. Influence of Local Geology on Earthquake Ground Motion. Ph.D., Thesis, 1969.

Trifunac M. D. Surface motion of a semicylindrical alluvial valley for incident plane SH wave. Bulletin of Seismological and Society of America, Vol. 61, Issue 6, 1971, p. 17551770.

Trifunac M. D. Scattering of plane SH wave by a semicylindrical canyon. Earthquake Engineering and Structural Dynamics, Vol. 1, Issue 3, 1973, p. 267281.

Wong H. L., Trifunac M. D. Surface motion of a semielliptical alluvial valley for incident plane SH waves. Bulletin of the Seismological Society of America, Vol. 64, Issue 5, 1974, p. 13891408.

Wong H. L., Trifunac M. D. Scattering of plane SH waves by a semielliptical canyon. Earthquake Engineering and Structural Dynamics, Vol. 3, Issue 2, 1974, p. 157169.

Liang J. W., Luo H., Lee V. W. Scattering of plane SH waves by a circular arc hill with a circular tunnel. Acta Seismologica Sinica, Vol. 17, Issue 5, 2004, p. 549563.

Lee V. W. On deformations near circular underground cavity subjected to incident plane SH waves. Proceedings of the Symposium for Applications of Computer Methods in Engineering, Los Angeles, 1977, p. 951961.

Lee VW., Luo H., Liang J. W. Antiplane (SH) waves diffraction by a semicircular cylindrical hill revisited: an improved accurate wave series analytic solution. Journal of Engineering Mechanics, ASCE, Vol. 132, Issue 10, 2006, p. 11061114.

Yuan X., Men F. L. Scattering of plane SH waves by a semicylindrical hill. Earthquake Engineering and Structural Dynamics, Vol. 21, Issue 12, 1992, p. 10911098.

Yuan X., Liao Z. P. Surface motion of a cylindrical hill of circulararc crosssection for incident plane SH waves. Soil Dynamics and Earthquake Engineering, Vol. 15, Issue 3, 1996, p. 189199.

Yuan X., Liao Z. P. Scattering of plane SH waves by arbitrary circularly raised topography. Earthquake Engineering and Engineering Vibration, Vol. 2, 1996, p. 113.

Vincent W., Lee., Luo H., Liang J. W. Antiplane (SH) Waves diffraction by a semicircular cylindrical hill revisited: an improved analytic wave series solution. Journal of Engineering Mechanics, Vol. 132, Issue 10, 2006, p. 11061114.

Vincent W., et al. Diffraction of antiplane SH waves by a semicircular cylindrical hill with an inside concentric semicircular tunnel. Earthquake Engineering and Engineering Vibration, Vol. 3, 2004, p. 249262.

Li Y. L., Yuan X., Sun R. Out of plane dynamic response of a heterogeneous hill in half space: the closedform solution. Earthquake Engineering and Engineering Vibration, Vol. 25, Issue 3, 2005, p. 15.

Tsaur D. H., Chang K. H. Scattering and focusing of SH waves by a convex circulararc topography. Geophysical Journal of the Royal Astronomical Society, Vol. 177, Issue 1, 2009, p. 222234.

Tsaur D. H. Scattering and focusing of SH waves by a lower semielliptic convex topography. Bulletin of the Seismological Society of America, Vol. 101, Issue 5, 2011, p. 22122219.

Amornwongpaibun A., Luo H., Lee V. W. Scattering of antiplane (SH) waves by a shallow semielliptical hill with a concentric elliptical tunnel. Journal of Earthquake Engineering, Vol. 20, Issue 3, 2016, p. 20.

Liu Z., Zhang H., Cheng A., et al. Seismic interaction between a lined tunnel and a hill under plane SV waves by IBEM. International Journal of Structural Stability and Dynamics, Vol. 19, Issue 2, 2018, p. 247293.

Pao Y. H., Mow C. C., Achenbach J. D. Diffraction of elastic waves and dynamic stress concentrations. Journal of Applied Mechanics, Vol. 40, Issue 4, 1973, p. 872.

Lee V. W., Cao H. Diffraction of SV waves by circular cylindrical canyons of various depths. Journal of Engineering Mechanics, Vol. 11, Issue 8, 1989, p. 445456.

Lee V. W., Karl J. Diffraction of SV waves by underground, circular, cylindrical cavities. Soil Dynamics and Earthquake Engineering, Vol. 11, Issue 8, 1992, p. 445456.

Lee V. W., Karl J. Diffraction of elastic plane P waves by circular, underground unlined tunnels. European Earthquake Engineering, Vol. 6, Issue 1, 1993, p. 2936.

Lee V. W., Wu X. Application of the weighted residual method to diffraction by 2D canyons of arbitrary shape: I. Incident SH waves. Soil Dynamics and Earthquake Engineering, Vol. 13, Issue 5, 1994, p. 365373.

Lee V. W., Wu X. Application of the weighted residual method to diffraction by 2D canyons of arbitrary shape: Ⅱ. Incident P, SV, and Rayleigh waves. Soil Dynamics and Earthquake Engineering, Vol. 13, Issue 5, 1994, p. 355364.