Next Article in Journal
The Relationship between Innovative Technology and Driver’s Resistance and Acceptance Intention for Sustainable Use of Automobile Self-Driving System
Next Article in Special Issue
Seismic Performance of a Sliding Isolation Bridge System with a New Spring Re-Centering Device
Previous Article in Journal
Mining Height Effect and Application of Upper Protected Layer Mining Pressure Relief
Previous Article in Special Issue
Evaluation of the Surface Performance of Mortar Matrix Subjected to Sodium Chloride Solution Modified with Hybrid Nanosilica Cement Paste
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Propagation of Nonplanar SH Waves Emanating from a Fault Source around a Lined Tunnel

1
Key Laboratory of Geomechanics and Embankment Engineering of Ministry of Education, Hohai University, Nanjing 210024, China
2
College of Civil and Transportation Engineering, Hohai University, Nanjing 210098, China
3
School of Civil Engineering, Suzhou University of Science and Technology, Suzhou 215011, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(16), 10127; https://doi.org/10.3390/su141610127
Submission received: 20 June 2022 / Revised: 23 July 2022 / Accepted: 25 July 2022 / Published: 15 August 2022
(This article belongs to the Special Issue Soil Dynamics and Earthquake Engineering in Sustainability)

Abstract

:
An analytical solution is presented for scattering nonplanar SH waves emanating from a fault source using a lined tunnel. The lined tunnel is assumed to be an annular elastic solid in half-space. A simplified circular arc fault model is employed to represent the wave source. By means of the separation of variables method, all wave fields are given in terms of the wave function series with unknown coefficients. Taking advantage of the method of images, the zero-stress boundary condition on the horizontal ground surface is satisfied automatically. By applying Graf’s addition formula, a system of equations for seeking the unknowns is derived by taking advantage of the boundary conditions. The problem of wave scattering is finally solved after seeking solutions for the system of equations through standard matrix techniques. The effects of fault distance, fault curvature, and fault orientation are revealed with numerical results. It is found that the plane waves provide a good approximation to the fault-induced cylindrical waves when the source-receiver distance or fault radius of curvature is larger enough. Fault-induced topographic effects are strongly affected by source orientation.

1. Introduction

Nowadays, as so many underground structures, such as subways, water-conveyance tunnels, mineral prospecting, and military facilities, have been constructed all around the world, the dynamic response of underground tunnels subjected to seismic waves has attracted more and more attention from geophysicists, seismologists, and civil engineers. Long underground tunnels are more likely to cross faults. Investigations of the water conveyance tunnel of Shih-Gang Dam affected by the Chi-Chi Earthquake (1999) show that tunnels passing through the displaced fault zone suffered catastrophic damage, and the lining was sheared off [1]. Similar damages are found in numerous tunnels during the Wenchuan Earthquake (2008) due to fault displacements, e.g., Longxi tunnel, Shaohuoping tunnel, Longchi tunnel, Chediguan tunnel, and Taoguan tunnel [2]. Therefore, there is an urgent need to study the scattering of nonplanar SH waves emanating from a fault by underground tunnels. The seismic response of tunnels is affected by the incident seismic wave, the reflected wave caused by the ground surface, and the scattered wave caused by the tunnel. Therefore, the wave scattering effect should be considered in the study of the seismic response of tunnels.
There have been a series of analytical studies on the seismic response of tunnels. Pao and Mow [3] first proposed a wave function series solution for the dynamic stress concentration of circular cavities in the whole space under the incidence of elastic SH waves. Then, Lee and Trifunac [4] extended the solution to the half-space case. The wave scattering problem is much more complicated in the case of incident P or SV waves due to the presence of mode conversion when the incident wave is reflected. Lee and Karl [5] obtained a series of analytical solutions to the scattering of P and SV waves induced by an underground circular cavity using a large arc assumption. In addition, many other scholars have analytically studied the influence of tunnel shape [6,7,8]. However, numerous studies have been conducted on the effects excited by plane harmonic waves [3,4,5,6,7,8,9,10,11,12,13,14], while the incidence of near-source waves is still at a research stage. Most of the previous achievements of near-source effects consider incident waves to be line source cylindrical waves. In their pioneering work, Pao and Mow [3] investigated the dynamic stress concentration of a circular cavity in the entire space for incident cylindrical P waves. Smerzini et al. [15] studied the effect of underground cavities on earthquake ground motions under plane and cylindrical SH waves. Gao et al. [8] provided a rigorous solution for scattering plane and cylindrical SH waves using a horseshoe-shaped cavity. This assumption was then extended to studies of site effects on seismic waves. Gao and Zhang [16] provided a solution of wave functions for the cylindrical SH wave scattering induced by a symmetrical V-shaped canyon. Zhang et al. [17] investigated near-source site effects based on a partially filled semicircular alluvial valley. But most earthquakes occur on the faults, i.e., the Chi-Chi earthquake (1999) resulted from the reactivation of the Chelungpu Fault [1], and major shock of Wenchuan earthquake (2008) occurred on the Beichuan–Yingxiu Fault [18]. Although a line source assumption is effective in simulating a near-field earthquake, there is a notable scarcity of literature dealing with the effect of fault. Only Kara and Trifunac [19,20] investigated the scattering of semi-circular sedimentary valleys considering excitation as a circular arc fault.
The current work is motivated by a lack of studies regarding the fault-induced dynamic response of a lined tunnel. The method proposed by Kara and Trifunac [19,20] is extended herein to consider the configuration in which the fault has an arbitrary orientation. The main contributions of this paper are twofold. On the one hand, the solution to the scattering problem of a lined tunnel induced by a fault is first presented in this paper. Since the excitation by a fault is more real compared to the plane wave and the line source of the cylindrical wave, the results presented herein are very instructive for practical use. On the other hand, through the analysis of numerical results, one can gain a better understanding of the relationship between fault-induced nonplanar waves and plane waves.

2. Materials and Methods

The 2D model is depicted in Figure 1. It represents an infinitely long, circular-lined tunnel embedded in a homogeneous, isotropic, linearly elastic, semi-infinite medium. The lining is characterized by inner radius a, outer radius b, shear modulus μl and shear wave velocity cl, while the half-space is characterized by shear modulus μs and shear wave velocity cs. The center of the tunnel is located at point O1 with a distance d + b below the ground surface. The fault is considered an arch of circular O3 with radius a1. In coordinate system O, the fault is located at r = rf, θ = −αf. In coordinate system O3, the fault is located at r3 = a1, between the angles 2π − α2 and 2π − α1. To facilitate the solution to the problem, an imaginary tunnel and an imaginary fault with respect to the ground surface are introduced herein. The definitions of the five Cartesian and five polar coordinate systems are shown in Figure 1.
The out-of-plane wave fields of all regions in the half-space must satisfy the governing equation of motion, i.e., Helmholtz equation [21].
2 u r 2 + 1 r u r + 1 r 2 2 u θ 2 + k 2 u = 0
where k = ω/c is the shear wave number.
Using the method of separating variables [22,23], the general solution of Equation (1) can be written as:
u r , θ = n = I n C n k r e i n θ
where In is the complex coefficient to be determined, Cn is the Bessel function with nth order, which can be either Jn, Yn, Hn(1), Hn(2) depending on the physical conditions of the problem.
Considering the image method and the Sommerfeld radiation condition, the wave field for each region in the sub-space can be written as follows. The wave field in the lining, satisfying the Sommerfeld radiation condition, is the standing wave in the finite domain and is formed by the Bessel function of the first and second kinds, written as
u l = n = A n J n k l r 1 e i n θ 1 + n = B n Y n k l r 1 e i n θ 1
The scattering wave field generated by the tunnel and its image, characterizing the radial oscillation from the site to the infinite boundary, is formed by the Hankel function of the first and second kinds, written as:
u s = n = D n H n 1 k s r 1 e i n θ 1 + n = D n H n 1 k s r 2 e i n θ 2
Since the displacement is finite within O3, the wave field inside the circle O3 can be written as:
u i = n = E n J n k s r 3 e i n θ 3
Similar to the scattering wave field expressed in Equation (4), the wave field outside circle O3 can be written as:
u o = n = F n H n 1 k s r f e i n θ f + n = F n H n 1 k s r f e i n θ f
In addition, all the waves in the sub-region should satisfy the following boundary conditions.
(1)
Stress-free boundary conditions on the flat surface and the inner surface of the lining
μ s r u s + u o θ = 0 , θ = ± π 2
μ l u l r 1 = 0 , r 1 = a
(2)
The continuity of both displacement and stress fields on the outer surfaces of the lining
u l r 1 , θ 1 = u s r 1 , θ 1 + u o r 1 , θ 1 , r 1 = b
μ l u l r 1 = μ s u s + u o r 1 , r 1 = b
(3)
The continuity of stress field on the circle O3
u i r f , θ f = u s r f , θ f + u o r f , θ f , r f = a 1
μ s u i r f = μ s u s + u o r f , r f = b
(4)
Assuming there is a unit-amplitude dislocation with out of plane motion, the boundary condition on the fault can be written as:
u i r f , θ f u s r f , θ f + u o r f , θ f = f θ f , r f = a 1
where f(θf) is a function written as:
f θ f = H θ 2 π α 1 + H θ 2 π α 2
where H(θ) is the Heaviside function, such that
H θ = 0 , θ < 0 1 , θ 0
The Fourier series expansion of f(θf) can be written as:
f θ f = n = f n e i n θ f
where fn is the Fourier transform of f(θf), which can be written as:
f n = 1 2 π 0 2 π f θ f e i n θ f d θ f = i 2 π n e i n α 1 e i n α 2 , n 0 1 2 π α 1 α 2 , n = 0
To solve the problem, Graf’s addition formula is used to transform the Bessel function from one coordinate to another.
Because of the introduction of the imaginary tunnel, the tress-free boundary condition on the flat surface is satisfied automatically. Substituting the general solution into other boundary conditions by considering the orthogonal properties of the complex exponential function, the following five equations can be obtained:
A q J q k l a + B q Y q k l a = 0
μ l k l [ A q J q k l b + B q Y q k l b ] = μ s k s [ D q H q 1 k s b + n = D n H n + q 1 k s 2 d J q k s b e i n q π 2 + n = F n H n + q 1 k s R i J q k s b e i q + n φ + i 3 2 n q π + n = F n H n q 1 k s R r J q k s b e i q φ + π 2 + i n φ + 3 π 2 ]
A q J q k l b + B q Y q k l b = D q H q 1 k s b + n = D n H n + q 1 k s 2 d J q k s b e i n q π 2 + n = F n H n + q 1 k s R i J q k s b e i q + n φ + i 3 2 n q π + n = F n H n q 1 k s R r J q k s b e i q φ + π 2 + i n φ + 3 π 2
E q J q k s a 3 = F q H q 1 k s a 3 + n = F n H n + q 1 k s 2 y 0 J q k s a 3 e i 3 2 n q π + n = D n H n + q 1 k s R i J q k s a 3 e i n + q φ + i 3 2 n q π + n = D n H n q 1 k s R r J q k s a 3 e i q φ + 3 π 2 + i n φ + π 2
E q J q k s a 3 + f q = F q H q 1 k s a 3 + n = F n H n + q 1 k s 2 y 0 J q k s a 3 e i 3 2 n q π + n = D n H n + q 1 k s R i J q k s a 3 e i n + q φ + i 3 2 n q π + n = D n H n q 1 k s R r J q k s a 3 e i q φ + 3 π 2 + i n φ + π 2
The unknown coefficients An, Bn, Dn, En and Fn can be obtained by combining Equations (18)–(22) and truncating the infinite series in these equations to a finite number. Once the coefficients are known, the wave fields in the half-space are ready for numerical computation. A detailed flowchart for the summarization of the present method is shown in Figure 2.

3. Numerical Results and Discussion

As for earthquake engineering, researchers are more interested in the displacement amplitude |u|. It is defined as:
u = Re ( u ] 2 + I m ( u ] 2
where Re(·) and Im(·) are the real and imaginary parts of a complex expression, respectively.
Furthermore, it is convenient to define the dimensionless frequency η as:
η = ω b π c s = 2 b λ s = k s b π
where λs is the wavelength of the incident SH waves.
The parameters of the cavity studied below are chosen as a/b = 0.9, d/b = 3, and the material properties cs/cl = 2/5, μs/μl = 8/75 [15].
There has been strict mathematical proof that the line source cylindrical waves can be reduced to plane waves for a sufficiently large source–receiver distance [16,17]. If the excitation is changed as a finite fault, however, does the same conclusion stand? To facilitate the comparison, all transfer functions for excitation by the fault are normalized by letting the free-field displacement amplitude at point O equal to 2. The effects of fault distance on the displacement amplitudes at three representative points on the ground surface are shown in Figure 3. To illustrate the asymptotic relationship, the corresponding results for plane waves (gray dashed lines) are given as well. Figure 3a,b shows the displacement amplitudes for the vertical incidence and oblique incidence for η = 1. Figure 3c,d correspond to the results for η = 4. As expected, the displacements at three points induced by the fault approach the corresponding values of plane waves as the fault moves away. That is to say, similar to the line source cylindrical waves, fault-induced nonplanar waves are reduced to plane waves for a sufficiently large source-receiver distance as well.
To facilitate this further, the normalized displacement on the ground surface is compared with the previous solution induced by plane waves. Figure 4 shows the comparison of ground motion between the present assumption and that of plane waves presented by Lee & Trifunac (1979) [4]. The parameters are set as a/b = 0.9, d/b = 1.5, cs/cl = 1, μs/μl = 1/3, α = 0°, φ = 90°. One can see that the present results are in excellent agreement with those by Lee & Trifunac (1979) [4] for large source-receiver distance.
To investigate the effect of distance between fault and tunnel on the dynamic response of tunnel, Figure 5 and Figure 6 exhibit the dynamic stress concentration factor (DSCF) along the inner surface and outer surfaces of the lining for different distances between fault and tunnel. To keep the fault curvature and orientation consistent with each other, α = 90° is set for φ = 0°, α = 45° for φ = 45°, α = 0° for φ = 90°, and a1/b = 8 for each case. Figure 5 represents the case for η = 1, and Figure 6 represents the case for η = 2.
The definition of the DSCF in this study is written as [3]:
DSCF = σ θ z μ l k l 2
where σθz represents the shear stress in the lining.
One can observe bilateral symmetry in the distribution curves of the DSCF for vertical incident waves (α = 90° and φ = 0°). Most of the distribution curves are approximately in the shape of a clover with four local maxima. The distribution of the DSCF along the inner and outer surface of the lining is similar. However, their values are different. The DSCF along the lower part of the tunnel is larger than that along the upper part of the tunnel. The reason is that the existence of the tunnel obstructs the transmission of seismic waves. Similar trends can be seen in the case of horizontal incident waves (α = 0° and φ = 90°). The DSCF is larger along the left part of the tunnel than that along the right part of the tunnel. Similar to the pattern of incident plane waves, when the dimensionless frequency η increases, the ability of short waves to penetrate the tunnel becomes weaker and weaker. Hence, the shielding effect is more significant for η = 2. Comparing the DSCF of different source-receiver distances, one can find that the DSCF is larger for smaller source-receiver distances than for larger source-receiver distances overall. This means that the damage to the near-fault earthquake is larger for the far-field earthquake.
To investigate the effect of fault orientation α on the dynamic response of tunnels, Figure 7 and Figure 8 exhibit the DSCF along the inner surface and outer surfaces of the lining for different fault orientations. The parameters are set as a1/b = 8 and Ri/b = 30. Figure 7 represents the case for η = 1, and Figure 8 represents the case for η = 2. Fault orientation plays an important role in DSCF. Comparing Figure 7a to Figure 7d, one can find that the location of the maximum DSCF changed from the lower part of the tunnel to the upper part of the tunnel as α changed from 90° to 0°. Similar trends can be seen by comparing Figure 7c to Figure 7f. However, the location of the maximum DSCF is almost unchanged compared to Figure 7b to Figure 7e. That means, different from the plane wave assumption, the fault orientation should be considered for the seismic design of tunnels.
Figure 9 illustrates the ground surface displacement amplitudes as a function of the dimensionless distance x/b and dimensionless frequency η for the same location but different fault orientations. For the case of vertical incidence (αi = 0) in Figure 9a, the shaded area, where the displacement amplitudes reduce, appears on the ground surface right over the tunnel. What’s more, with the increase of dimensionless frequency η, the reduction of displacement amplitudes becomes more and more obvious. For cases of non-vertical incidence (αi ≠ 0, Figure 9d–f), one can also see the illuminated area on the left ground and the shaded area on the right ground. The motions on the ground surface decrease quickly from the illuminated area to the shaded area. Comparing the ground motion for different fault orientations, one can find that there is almost no difference in the ground motion for horizontal fault orientation ((α1 + α2)/2 = 90°) and vertical fault ((α1 + α2)/2 = 0°). However, the ground motion for the oblique fault ((α1 + α2)/2 = 60°) is obviously different from the vertical and horizontal cases. Moreover, the ground motion over the illuminated area for oblique faults is larger than that for vertical and horizontal faults. This means that the damage to the oblique fault is greater than the vertical and horizontal faults.

4. Conclusions

The scattering of nonplanar SH waves emanating from a finite source by a lined tunnel is studied in this paper. A simplified circular arc fault model is adopted. The influence of fault distance, fault curvature, and fault orientation are considered in the study. Similar to the line source cylindrical waves, the results of fault-induced nonplanar waves tend to be equivalent to those of plane waves for a sufficiently large source–receiver distance. Through the investigation of DSCF along the inner and outer surfaces of the lining, one can find that the distribution of the DSCF along the inner and outer surface of lining is similar. However, the values are slightly different. Fault orientation and source-receiver distance play an important role in the DSCF of lining. It was found that the tunnel can obstruct the transmission of seismic waves and make the DSCF smaller at the shaded side of the tunnel. The damage to the near-fault earthquake is larger for the far-field earthquake. Through the investigation of ground motion, one can find that fault orientation can significantly change ground motion for certain fault orientations. The damage to the oblique fault is greater than the vertical and horizontal faults.

Author Contributions

Conceptualization, N.Z. and D.D.; methodology, N.Z., D.D., Y.Z. (Yu Zhang) and X.C.; software, Y.Z. (Yunfei Zhang), Y.Z. (Yu Zhang) and X.C.; validation, N.Z., B.S. and D.D.; formal analysis, Y.Z. (Yunfei Zhang), Y.Z. (Yu Zhang) and X.C.; investigation, Y.Z. (Yunfei Zhang), Y.Z. (Yu Zhang) and X.C.; resources, N.Z. and D.D.; data curation, D.D.; writing—original draft preparation, Y.Z. (Yunfei Zhang); writing—review and editing, N.Z. and D.D.; visualization, Y.Z. (Yunfei Zhang); supervision, N.Z.; project administration, N.Z., D.D. and X.C.; funding acquisition, N.Z., D.D. and X.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant numbers 52078184 and 52108314; the Fundamental Research Funds for the Central Universities, grant number B210202042; the China Postdoctoral Science Foundation funded project, grant number 2020M681475; and the Natural Science Foundation of Jiangsu Province, grant number BK20220638.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

ainner radius of tunnel
bouter radius of tunnel
dburial depth of tunnel
Risource-receiver distance
μlshear modulus of lining
clshear wave velocity of lining
μsshear modulus of half-space
csshear wave velocity of half-space
klshear wave number of lining
ksshear wave number of half-space
λswavelength of incident SH waves
Jn, Yn, Hn(1), Hn(2)Bessel function
An, Bn, Dn, En, Fnunknown coefficients
ηdimensionless frequency

References

  1. Wang, W.L.; Wang, T.T.; Su, J.J.; Lin, C.H.; Seng, C.R.; Huang, T.H. Assessment of damage in mountain tunnels due to the Taiwan Chi-Chi earthquake. Tunn. Undergr. Space Technol. 2001, 16, 133–150. [Google Scholar] [CrossRef]
  2. Wang, Z.Z.; Gao, B.; Jiang, Y.; Yuan, S. Investigation and assessment on mountain tunnels and geotechnical damage after the Wenchuan earthquake. Sci. China Ser. E 2009, 52, 546–558. [Google Scholar] [CrossRef]
  3. Pao, Y.H.; Mow, C.C. Diffraction of Elastic Waves and Dynamic Stress Concentrations; Crane Russak: New York, NY, USA, 1973. [Google Scholar]
  4. Lee, V.W.; Trifunac, M.D. Response of tunnels to incident SH-waves. J. Eng. Mech. Div. 1979, 105, 643–659. [Google Scholar] [CrossRef]
  5. Lee, V.W.; Karl, J. Diffraction of SV waves by underground, circular, cylindrical cavities. Soil Dyn. Earthq. Eng. 1992, 11, 445–456. [Google Scholar] [CrossRef]
  6. Liang, J.W.; Luo, H.; Lee, V.W. Diffraction of plane SH waves by a semi-circular cavity in half-space. Earthq. Sci. 2010, 23, 5–12. [Google Scholar] [CrossRef]
  7. Tsaur, D.H.; Chang, K.H. Multiple scattering of SH waves by an embedded truncated circular cavity. J. Mar. Sci. Technol. 2012, 20, 9. [Google Scholar] [CrossRef]
  8. Gao, Y.F.; Dai, D.H.; Zhang, N.; Wu, Y.X.; Mahfouz, A.H. Scattering of plane and cylindrical SH waves by a horseshoe shaped cavity. J. Earthq. Tsunami 2017, 11, 1650011. [Google Scholar] [CrossRef]
  9. Gao, Y.F.; Chen, X.; Zhang, N.; Dai, D.H.; Yu, X. Scattering of plane SH waves induced by a semicylindrical canyon with a subsurface circular lined tunnel. Int. J. Geomech. 2018, 18, 06018012. [Google Scholar] [CrossRef]
  10. Chen, X.; Zhang, N.; Gao, Y.F.; Dai, D.H. Effects of a V-shaped canyon with a circular underground structure on surface ground motions under SH wave propagation. Soil Dyn. Earthq. Eng. 2019, 127, 105830. [Google Scholar] [CrossRef]
  11. Zhang, N.; Chen, X.; Gao, Y.F.; Dai, D.H. Analytical solution to scattering of SH waves by a circular lined tunnel embedded in a semi-circular alluvial valley in an elastic half–space. Tunn. Undergr. Space Technol. 2020, 106, 103615. [Google Scholar] [CrossRef]
  12. Zhang, N.; Zhang, Y.; Dai, D.H. Dynamic response of a functionally graded cylindrical tube with power-law varying properties due to SH-waves. Waves Random Complex Media 2021, 1–19. [Google Scholar] [CrossRef]
  13. Zhang, N.; Zhang, Y.; Gao, Y.F.; Dai, D.H.; Huang, C.X. Effect of imperfect interfaces on dynamic response of a composite lining tunnel with an isolation layer under plane P and SV waves. Soil Dyn. Earthq. Eng. 2021, 142, 106586. [Google Scholar] [CrossRef]
  14. Jin, L.; Zhu, J.; Zhou, W.; Liang, J.W.; Chen, G.X. 2D dynamic tunnel-soil-aboveground building interaction I: Analytical solution for incident plane SH-waves based on rigid tunnel and foundation model. Tunn. Undergr. Space Technol. 2022, 128, 104625. [Google Scholar] [CrossRef]
  15. Smerzini, C.; Aviles, J.; Paolucci, R.; Sánchez-Sesma, F.J. Effect of underground cavities on surface earthquake ground motion under SH wave propagation. Earthq. Eng. Struct. Dyn. 2009, 38, 1441–1460. [Google Scholar] [CrossRef]
  16. Gao, Y.F.; Zhang, N. Scattering of cylindrical SH waves induced by a symmetrical V-shaped canyon: Near-source topographic effects. Geophys. J. Int. 2013, 193, 874–885. [Google Scholar] [CrossRef]
  17. Zhang, N.; Gao, Y.F.; Yang, J.; Xu, C.J. An analytical solution to the scattering of cylindrical SH waves by a partially filled semi-circular alluvial valley: Near-source site effects. Earthq. Eng. Eng. Vib. 2015, 14, 189–201. [Google Scholar] [CrossRef]
  18. Cui, P.; Chen, X.Q.; Zhu, Y.Y.; Su, F.H.; Wei, F.Q.; Han, Y.S.; Liu, H.J.; Zhuang, J.Q. The Wenchuan earthquake (May 12, 2008), Sichuan province, China, and resulting geohazards. Nat. Hazards 2011, 56, 19–36. [Google Scholar] [CrossRef]
  19. Kara, H.F.; Trifunac, M.D. A note on plane-wave approximation. Soil Dyn. Earthq. Eng. 2013, 51, 9–13. [Google Scholar] [CrossRef]
  20. Kara, H.F.; Trifunac, M.D. Two-dimensional earthquake vibrations in sedimentary basins–SH waves. Soil Dyn. Earthq. Eng. 2014, 63, 69–82. [Google Scholar] [CrossRef]
  21. Zhang, N.; Pan, J.S.; Gao, Y.F.; Chen, X.; Dai, D.H.; Zhang, Y. Surface motion of an arbitrary number of semi-circular viscoelastic alluvial valleys for incident plane SH waves. Geophys. J. Int. 2022, 228, 1607–1620. [Google Scholar] [CrossRef]
  22. Dai, D.H.; El Naggar, M.H.; Zhang, N.; Wang, Z.B. Rigorous solution for kinematic response of floating piles to vertically propagating S-waves. Comput. Geotech. 2021, 137, 104270. [Google Scholar] [CrossRef]
  23. Dai, D.H.; El Naggar, M.H.; Zhang, N.; Wang, Z.B. Rigorous solution for kinematic response of floating piles subjected to vertical P-wave. Appl. Math. Model. 2022, 106, 114–125. [Google Scholar] [CrossRef]
Figure 1. 2D model of a lined tunnel subjected to a fault-induced nonplanar wave.
Figure 1. 2D model of a lined tunnel subjected to a fault-induced nonplanar wave.
Sustainability 14 10127 g001
Figure 2. Flowchart for the summarization of the present method.
Figure 2. Flowchart for the summarization of the present method.
Sustainability 14 10127 g002
Figure 3. Effects of source distance rf on surface displacement amplitudes at three positions.
Figure 3. Effects of source distance rf on surface displacement amplitudes at three positions.
Sustainability 14 10127 g003
Figure 4. Comparison of ground motion between the present assumption and that of plane waves presented by Lee & Trifunac (1979) [4].
Figure 4. Comparison of ground motion between the present assumption and that of plane waves presented by Lee & Trifunac (1979) [4].
Sustainability 14 10127 g004
Figure 5. Effects of distance between fault and tunnel on DSCF (η = 1).
Figure 5. Effects of distance between fault and tunnel on DSCF (η = 1).
Sustainability 14 10127 g005
Figure 6. Effects of distance between fault and tunnel on DSCF (η = 2).
Figure 6. Effects of distance between fault and tunnel on DSCF (η = 2).
Sustainability 14 10127 g006
Figure 7. Effects of fault orientation on DSCF (η = 1).
Figure 7. Effects of fault orientation on DSCF (η = 1).
Sustainability 14 10127 g007
Figure 8. Effects of fault orientation on DSCF (η = 2).
Figure 8. Effects of fault orientation on DSCF (η = 2).
Sustainability 14 10127 g008
Figure 9. 3D plots of ground surface displacement amplitudes for different fault orientations.
Figure 9. 3D plots of ground surface displacement amplitudes for different fault orientations.
Sustainability 14 10127 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, N.; Zhang, Y.; Dai, D.; Zhang, Y.; Sun, B.; Chen, X. Propagation of Nonplanar SH Waves Emanating from a Fault Source around a Lined Tunnel. Sustainability 2022, 14, 10127. https://doi.org/10.3390/su141610127

AMA Style

Zhang N, Zhang Y, Dai D, Zhang Y, Sun B, Chen X. Propagation of Nonplanar SH Waves Emanating from a Fault Source around a Lined Tunnel. Sustainability. 2022; 14(16):10127. https://doi.org/10.3390/su141610127

Chicago/Turabian Style

Zhang, Ning, Yunfei Zhang, Denghui Dai, Yu Zhang, Baoyin Sun, and Xin Chen. 2022. "Propagation of Nonplanar SH Waves Emanating from a Fault Source around a Lined Tunnel" Sustainability 14, no. 16: 10127. https://doi.org/10.3390/su141610127

APA Style

Zhang, N., Zhang, Y., Dai, D., Zhang, Y., Sun, B., & Chen, X. (2022). Propagation of Nonplanar SH Waves Emanating from a Fault Source around a Lined Tunnel. Sustainability, 14(16), 10127. https://doi.org/10.3390/su141610127

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop