Preprint
Technical Note

Response of a Triangular Hill Joined with a Semicircular Alluvial Valley under Incident SH Waves

Altmetrics

Downloads

73

Views

12

Comments

0

Submitted:

11 March 2024

Posted:

12 March 2024

You are already at the latest version

Alerts
Abstract
Under the same seismic load, different local topography will have different effects, which is an important research field in seismic engineering. Hills are connected with valleys is a very typical terrain. Taking this terrain as an example, its response to the most destructive SH waves in the earthquake is studied. By simplifying and dividing the model, the wave functions satisfying the wave equation and boundary conditions in each region are constructed by using the complex function and moving coordinate method. The equations satisfying the boundary continuity conditions are established. With the help of calculation software, the numerical results are obtained. The influence of different parameters on ground motion is discussed.
Keywords: 
Subject: Environmental and Earth Sciences  -   Geophysics and Geology

1. Introduction

Earthquake is a very serious natural disaster. It is a worldwide problem to accurately predict earthquakes. To reduce the damage caused by earthquakes, improving the seismic capacity of structures and evaluating the seismic damage is a key point in the field of earthquake engineering. Through previous studies, it is found that under the same seismic load, different local topography will have different effects.
Shear horizontal wave (SH wave) can make the surface shake in the horizontal direction, which is the main cause of earthquake disasters. Many scholars have done a lot of research on the response of different terrains to SH waves. Xu and Han [1] studied the scattering of plane SH waves by multiple semi-cylindrical canyons. Numerical results and discussions of displacements and phases on the surface of the canyons are given. By using the accurate region-matching technique to derive a rigorous series solution, Chang et al. [2] studied the scattering of SH waves by a circular sectorial canyon. Combining the region-matching technique and a periodic indirect boundary element method (PIBEM), Ba and Liang [3] discussed the scattering and diffraction of plane SH waves by periodic alluvial valleys embedded in a layered half-space. Yang et al. [4] investigated the scattering of SH waves around a circular canyon in radial inhomogeneous wedge space. By transforming the governing equation with variable coefficients into a standard Helmholtz equation, the corresponding analytical solution is derived. Qi et al. [5] discussed the scattering of steady-state shear horizontal guided waves from the elastic strip media with multiple semi-cylindrical depressions on the surface and obtained the analytical solution.
The real situation is usually very complex, which is the coupling of many different structures. In this paper, the scattering of SH waves on a hill joined with an alluvial valley is studied. The final solution of the surface response is obtained by using the analytical method.

2. Model Simplification and Partition

The hill is simplified as an isosceles triangle, and the alluvial valley is simplified as a semicircle. The simplified model and related parameters are shown in Figure 1. The horizontal surface, the hypotenuse of the triangle, and the boundary of the semicircle are marked as S, C, and S2, respectively. The mass density of the foundation and sediment is ρ1 and ρ2, and the shear modulus is µ1 and µ2. SH waves incident in the plane along the X1-O1-Y1 coordinate system at the angle of α.
The method of “partition and combination” is used to study this complex terrain. As shown in Figure 2, a semicircular virtual boundary S1 is created with O1 as the center and r1 as the radius. Taking S1 and S2 as the interface, the model is divided into 3 regions.
Wave functions satisfying the boundary conditions are constructed by using the wave function expansion method and coordinate moving method [6,7,8]. According to the continuous conditions of displacement and stress on the common boundary, the equations are established to solve the unknown coefficients in the wave function of each term. Finally, the dynamic response of each position on the surface can be obtained.

3. Construction of Wave Functions

For steady state, the governing equation of SH waves in isotropic homogeneous and continuous media can be expressed as [9,10]
Preprints 101101 i001
where W is the displacement function, k=ω/Cs, ω is frequency, Preprints 101101 i027 is the shear wave velocity.
In polar coordinates, the corresponding stress component is
Preprints 101101 i002

3.1. Wave Function in Region Ι

The zero-stress condition on the boundary C is
Preprints 101101 i003
The standing wave function satisfying the governing Equation (1) and boundary condition (3) can be expressed as [11]
Preprints 101101 i004
where Am is the constant to be solved, W0 is the displacement amplitude. p=π/(2θ0). Jmp(·) is a Bessel function with mp ranks.
The wave function needs to satisfy the continuity condition on S1 boundary [12,13]. According to the moving coordinate method, Equation (4) in Preprints 101101 i028 the coordinate system can be expressed as
Preprints 101101 i005
The corresponding stress can be expressed as
Preprints 101101 i006
where
P t j s = J t 1 k 1 s s s t 1 e i θ j J t + 1 k 1 s 1 t s s t + 1 e i θ j + J t 1 k 1 s 1 t s s 1 t e i θ j J t + 1 k 1 s s s t + 1 e i θ j

3.2. Wave Function in Region Ш

The boundary condition of the upper surface of sediment can be expressed as
Preprints 101101 i007
The standing wave function satisfying condition (8) and Equation (1) can be expressed as
Preprints 101101 i008
where Bm is the constant to be solved.
According to the moving coordinate method, Equation (9) in Preprints 101101 i029 the coordinate system can be expressed as
Preprints 101101 i009
The corresponding stress can be expressed as
Preprints 101101 i010

3.3. Scattering Waves in RegionⅡ

In region Ⅱ, there are two scattering waves on each boundary because the scattering waves generated by the two depressions interfere with each other. The functions that satisfy the governing equation and the boundary condition can be expressed in the form of multipolar coordinates as [14,15,16]
Preprints 101101 i011
Preprints 101101 i012
Preprints 101101 i013
Preprints 101101 i014
The corresponding stress can be expressed as
Preprints 101101 i015
Preprints 101101 i016
where Cm, Dm are the constant to be solved. H m 1 is Hankel function of the first kind, and
U t j s = H t 1 1 k 1 s s s t 1 e i θ j H t + 1 1 k 1 s 1 t s s t + 1 e i θ j + H t 1 1 k 1 s 1 t s s 1 t e i θ j H t + 1 1 k 1 s s s t + 1 e i θ j

3.4. Incident and Reflected Waves

In the complex plane Preprints 101101 i030, the incident and reflected waves can be expressed as
Preprints 101101 i017
The corresponding stress can be expressed as
Preprints 101101 i018
where Preprints 101101 i031, and
P t j s = J t 1 k 1 s s s t 1 e i θ j J t + 1 k 1 s 1 t s s t + 1 e i θ j + J t 1 k 1 s 1 t s s 1 t e i θ j J t + 1 k 1 s s s t + 1 e i θ j

4. Establishment of Equations

Wave functions in each region satisfy the continuous conditions of displacement and stress on the common boundary, thus the equations of unknown coefficients can be established as (22). Because there are many parameter variables, the solution is very complex. Due to the space limitation, this article only considers the special case of d=r2.
Preprints 101101 i019
Expanding (22) by Fourier series in the complex domain, the following equations are obtained
Preprints 101101 i020
The expression of each function is
Preprints 101101 i021
Preprints 101101 i022
Preprints 101101 i023
Preprints 101101 i024
(23) are four infinite series equations with unknown coefficients. Finally, the unknown coefficients are obtained by programming.

5. Results Analysis

There is only a standing wave in region Ⅰ and region Ⅲ, so the total wave field is
Preprints 101101 i025
The total wave field in region Ⅱ consists of three parts:
Preprints 101101 i026
The incident wave amplitude W0 is assumed to be 1.0. y1=0 represents the vertex of the hill, y1r1 represents the intersection of the hill and the horizontal plane, |y1|<r1 represents the points on the surface of the hill, y1/r1=d1 represents the center of the alluvial valley. Through numerical examples, the influence of different SH wave numbers (η=2r1/λ), incident angle, and other parameters on the amplitude of surface displacement is analyzed.
When the sediment is softer than the foundation, the variation of surface displacement amplitude with different wave numbers and incident angles is shown in Figure 3. When the wave number is small, the surface displacement changes little, which is similar to the quasi-static problem. With the increase of wave number, the change of surface displacement becomes more and more intense, and the maximum value of displacement amplitude gradually approaches the alluvial valley. When η≥0.5, the surface displacement amplitude in the valley is more than 10, which is 5 times that when the SH wave incident on the horizontal surface.
When the sediment is harder than the foundation, the surface displacement amplitude with different wave numbers and incident angles is shown in Figure 4, which is much smaller than in Figure 3. The surface displacement in the valley is generally smaller than that of the hill. The maximum value appears at the top of the hill (η=1, α=0°), which is 5.2. When α=0°, the ground motion of the hill increases with the decrease of the vertex angle, and the opposite is true for α>0°.
The average value of surface displacement on the hill is set as e1. When the SH waves only incident on the hill, the average value of surface displacement is set as e2. Let e=e1/e2. Taking the distance between hill and valley as abscissa and e as ordinate, the results are shown in Figure 5. With the increase of distance, e is closer to 1. When d1>30r1, the difference between e and 1 is less than 0.02. At this time, the interaction between the hill and the valley can be ignored.

6. Conclusions

Based on elastic wave dynamics and the knowledge of complex functions, the scattering of SH waves on a triangular hill joined with a semicircular alluvial valley is studied. The conclusions are as follows:
1) When the sediment is softer than the foundation, the surface displacement is relatively large as a whole. The maximum value appears near the center of the valley, which is more than 5 times when SH waves incident into the infinite half-space.
2) When the sediment is harder than the foundation, the influence of the valley on the hill is small. The maximum surface displacement appears at the top of the hill, which is more than 2 times when SH waves incident into the infinite half-space.
3) When the distance between the hill and the valley is more than 30 times the half-width of the bottom of the hill, the interaction can be ignored.

Author Contributions

Conceptualization, Jianze Liu and Guangzheng Wang; Funding acquisition, Fazhan Yang and Xuefeng Zhang; Methodology, Jianze Liu, Yang Li and Zhuang Wang; Project administration, Jianze Liu; Software, Jianze Liu and Guangzheng Wang; Visualization, Jianze Liu and Fazhan Yang; Writing – original draft, Jianze Liu and Zhuang Wang; Writing – review & editing, Jianze Liu, Yang Li and Zhuang Wang.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declared no potential conflicts of interest concerning the research, authorship, and/or publication of this article.

References

  1. Xu, Y.Y.; Han, F. Scattering of SH-waves by multiple semi-cylindrical canyons. Earthquake Engineering and Engineering Vibration 1992, 12, 12–18. [Google Scholar]
  2. Chang, K.H.; Tsaur, D.H.; Wang, J.H. Scattering of SH waves by a circular sectorial canyon. Geophysical Journal International 2013, 195, 532–543. [Google Scholar] [CrossRef]
  3. Ba, Z.N.; Liang, J.W. Dynamic response analysis of periodic alluvial valleys under incident plane SH-waves. Journal of Earthquake Engineering 2016, 21, 531–550. [Google Scholar] [CrossRef]
  4. Yang, Z.L.; Li, X.Z.; Song, Y.Q. Scattering of SH waves around circular canyon in inhomogeneous wedge space. Geophysical Journal International 2020, 223, 45–56. [Google Scholar] [CrossRef]
  5. Qi, H.; Yang, R.J.; Guo, J. Scattering of SH wave by multiple semi-cylindrical depressions in an elastic strip. Explosion and Shock Waves 2020, 40, 29–41. [Google Scholar]
  6. Shyu, W.S.; Teng, T.J.; Yeh, C.S. Surface motion of two canyons for incident SH waves by hybrid method. Procedia Engineering 2014, 79, 533–539. [Google Scholar] [CrossRef]
  7. Han, F.; Liu, D.K. Research of scattering of plane SH-wave by a semi-cylindrical canyon. Journal of Harbin University of Civil Engineering and Architecture 1990, 23, 24–31. [Google Scholar]
  8. Shashi, K.; Swapan, K.C. Influence of scattering of SH-waves in dynamic interaction of shear wall with soil layers. Earthquake Engineering and Engineering Vibration 2020, 19, 583–595. [Google Scholar]
  9. Lv, X.T.; An, J.B. Ground motion of a semi-cylindrical hill joined with a semi-cylindrical alluvial valley under incident SH-waves. Journal of Vibration and Shock 2014, 33, 127–131. [Google Scholar]
  10. Liu, Z.X.; Yu, Q.; He, Y. Effect of randomness of geotechnical medium on the seismic ground motion amplification effect of a sedimentary valley. Acta Seismologica Sinica 2017, 39, 764–777. [Google Scholar]
  11. Chen, S.H.; Zhang, Y.S. Dynamic responses of canyon with multiple arc-shaped layers under incidence of plane SV waves. Chinese Journal of Geotechnical Engineering 2017, 39, 1074–1081. [Google Scholar]
  12. Gai, Z.X.; Chen, C.X. Localized boundary element method for seismic wave surface scattering problems: 2D SH case. Chinese Journal of Geophysics 2020, 63, 2274–2280. [Google Scholar]
  13. Xu, H.N.; Zhang, Y.; Huang, Q.Y. Ground motion of underground complicated circular lining structure by out-plane wave. Earthquake Engineering and Engineering Dynamics 2019, 39, 154–159. [Google Scholar]
  14. Ma, R.; Li, Y.Q.; Jing, L.P. Influence of a semi-circular canyon site on ground surface displacement due to incidence of elastic waves. Earthquake Engineering Journal 2019, 41, 392–398. [Google Scholar]
  15. Liu, Q.J.; Wu, Z.Y.; Vincent, W.L. Scattering and reflection of SH waves around a slope on an elastic wedged space. Earthquake Engineering and Engineering Vibration 2019, 18, 255–266. [Google Scholar] [CrossRef]
  16. Zhou, F.X.; Pang, B.B.; Zhai, R.Z. Effect of circular tunnels on earthquake subground motions: Incidence of SH wave. Earthquake Engineering Journal 2018, 40, 873–878. [Google Scholar]
Figure 1. Numerical model of a triangular hill joined with a semi-circular alluvial valley.
Figure 1. Numerical model of a triangular hill joined with a semi-circular alluvial valley.
Preprints 101101 g001
Figure 2. Model of partition.
Figure 2. Model of partition.
Preprints 101101 g002
Figure 3. Surface displacement amplitude (μ1/μ2=6, ρ1/ρ2=3/2).
Figure 3. Surface displacement amplitude (μ1/μ2=6, ρ1/ρ2=3/2).
Preprints 101101 g003
Figure 4. Surface displacement amplitude (μ1/μ2=1/6, ρ1/ρ2=2/3).
Figure 4. Surface displacement amplitude (μ1/μ2=1/6, ρ1/ρ2=2/3).
Preprints 101101 g004
Figure 5. The influence of the alluvial valley on the hill with different distance.
Figure 5. The influence of the alluvial valley on the hill with different distance.
Preprints 101101 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated