Next Article in Journal
Personalized Video Summarization: A Comprehensive Survey of Methods and Datasets
Previous Article in Journal
Modification of Insulating Properties of Surfaces of Dielectric High-Voltage Devices Using Plasma
Previous Article in Special Issue
Field Reconnaissance and Earthquake Vulnerability of the RC Buildings in Adıyaman during 2023 Türkiye Earthquakes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Seismic Response of a Water-Conveyance Tunnel Considering Non-Uniform Longitudinal Subsurface Geometry and Obliquely Incident SV-Waves

1
Key Laboratory of Geotechnical Mechanics and Engineering of the Ministry of Water, Changjiang River Scientific Research Institute, Wuhan 430010, China
2
School of Civil Engineering and Architecture, Wuhan University of Technology, Wuhan 430070, China
3
China Communications Construction Company Second Harbor Engineering Company Ltd., Wuhan 430040, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(11), 4398; https://doi.org/10.3390/app14114398
Submission received: 15 March 2024 / Revised: 4 May 2024 / Accepted: 15 May 2024 / Published: 22 May 2024
(This article belongs to the Special Issue Earthquake Engineering and Seismic Risk)

Abstract

:
The longitudinal seismic response characteristics of a shallow-buried water-conveyance tunnel under non-uniform longitudinal subsurface geometry and obliquely incident SV-waves was studied using the numerical method, where the effect of the non-uniform longitudinal subsurface geometry due to the existence of a local one-sided rock mountain on the seismic response of the tunnel was focused on. Correspondingly, a large-scale three-dimensional (3D) finite-element model was established, where different incidence angles and incidence directions of the SV-wave were taken into consideration. Also, the non-linearity of soil and rock and the damage plastic of the concrete lining were incorporated. In addition, the wave field of the site and the acceleration response as well as damage of the tunnel were observed. The results revealed the following: (i) a local inclined subsurface geometry may focus an obliquely incident wave due to refraction or total reflection; (ii) a tunnel in a site adjacent to a rock mountain may exhibit a higher acceleration response than a tunnel in a homogeneous plain site; and (iii) damage in the tunnel in the site adjacent to a rock mountain may appear in multiple positions, and the effect of the incidence angle on the mode of compressive deformation and damage of the lining is of significance.

1. Introduction

Several earthquake events have caused heavy human tolls and substantial property damages [1,2,3] in recent years, and, thus, the seismic response characteristics of tunnels have attracted extensive attention. Many researchers have studied the dynamic response of underground tunnels using different methods [4], and, due to the complexity of wave propagation and scattering in soil media, the dynamic analysis of underground structures is always carried out economically and under realistic conditions with the aid of the numerical method [5,6,7].
The angle of incidence is a critical factor that determines the seismic tunnel response [8,9,10,11,12,13], and the influence of the incidence direction and the angles of the seismic wave on the response of a tunnel was first noticed in the anti-seismic studies of underground tunnels. Through a series of studies and analyses, Huang et al. [7] stated that the incidence direction of earthquake waves has significance on the seismic responses of long lined tunnels, and it is important to include the incidence angles of earthquake waves when mapping the seismic risk of tunnels. Li and Song [14] developed a novel numerical method for analyzing the longitudinal seismic response of tunnels under an asynchronous earthquake wave input, based on which the influence of the incidence angle and incidence direction on the longitudinal seismic response of tunnels was discussed. Considering obliquely incident seismic waves, Zhu et al. [15,16] carried out a study on the dynamic response of tunnels using a 2.5D coupled finite-element—boundary-element model. Zhou et al. [17] conducted a 3D seismic response analysis of an immersed tunnel in a horizontally layered site subjected to obliquely incident waves based on the precise dynamic stiffness matrix of the soil layer and half-space via combined viscous–spring boundary and equivalent node stress methods. Zhu et al. [18] investigated the effects of pulse-like ground motions on tunnels in saturated poroelastic soil for obliquely incident seismic waves using a 2.5D modeling technique. Considering 3D obliquely incident P- and SV-waves, Yang et al. [19] calculated the 3D seismic response of a lined tunnel embedded in a half-space using the 2.5D finite/infinite-element method. Huang et al. [20] investigated the dynamic performance of prefabricated utility tunnels considering the oblique incidence of P-waves. Huang et al. [21] developed a viscous–slip interface model to reflect the contact between the surrounding rock and tunnels, and the dynamic stress concentrations and surface displacements of two shallowly buried lining tunnels were examined using IBIEM, during which the effects of the angle on the medium’s responses were considered. Zhang et al. [22] evaluated the seismic response of an immersed tunnel in an ocean space under obliquely incident earthquake waves.
In addition, the responses of a tunnel in the longitudinal direction should not be neglected [23], especially when considering the tunnel lining damage which has been caused by recent earthquakes, and the longitudinal seismic response is another important perspective in the field of anti-seismic analyses of tunnels. Considering the longitudinal variation in ground properties, Park et al. [24] developed a longitudinal displacement profile-based procedure for simulating the dynamic response of tunnels and performed a series of pseudo-static 3D finite-element analyses of tunnels under spatially varying ground motion. Based on the generalized response displacement method, Miao et al. [25,26] studied the longitudinal seismic response of a metro tunnel in a site with non-uniform soil properties and found that the inter-segment opening of a tunnel in soft soil is larger than that in hard rock. Liu et al. [27] proposed a new multi-scale numerical model and calculated the longitudinal seismic response of a shield tunnel through a typical site with a soft–hard strata junction, and it was found that the dynamic response of the shield tunnel rose dramatically near the soft–hard strata interface because of the non-uniform deformation of the soil in that area. Yuan et al. [28], Yu et al. [29], and Yan et al. [30] carried out a series of shaking-table tests to investigate the longitudinal seismic response of a tunnel under non-uniform ground motion.
As a matter of fact, many circular underground structures such as cavities and tunnels have been constructed near various canyons in seismic zones [31], and the local subsurface geometry may possess a remarkable effect on the seismic response of underground structures, and the longitudinal non-uniform subsurface geometry (i.e., subsurface geometry variation) may affect the longitudinal seismic response characteristics of tunnels. However, although studies on the seismic response analyses of tunnels are extensive, the longitudinal dynamic response and damage characteristics of tunnels under non-uniform longitudinal subsurface geometry are not well known and understood. Therefore, in the present paper, taking the incidence angle and direction of SV-waves into account, the authors carried out a numerical study of this problem. The rest of this paper is organized as follows. The definition of the problem discussed in this paper is presented in Section 2. In Section 3, the setup of the numerical models is presented in detail. Subsequently, Section 4 offers the results of the longitudinal dynamic responses of the tunnels. Finally, concluding remarks are given in Section 5.

2. Problem Definitions

The Dianzhong Water Diversion Project is the largest water diversion project under construction in China. It is located in the Yunnan province, China, and passes through Lijiang city, the Dali Prefecture, the Chuxiong Prefecture, Kunming city, Yuxi city, and the Honghe Prefecture. The total length of the project is about 664 km, more than 90% of which is tunnel engineering. What is noticeable is that multiple intervals of tunnel engineering in Kunming city are in soft soils, significantly near a rock mountain, as shown in Figure 1, and the burial depth of the tunnel is generally small (tens of meters). In other words, the tunnel might pass through a plain soft-soil site first, then a soft-soil site near a rock mountain, and, subsequently, a plain soft-soil site again, which leads to the fact that the longitudinal non-uniform subsurface geometry along the tunnel, induced by a local rock mountain, may possess a considerable influence on the seismic response of said tunnel. The existence of a local one-sided rock mountain near the tunnel may affect the propagation of incident seismic waves and, thus, brings about differences in the seismic responses between tunnels with and without adjacent rock mountains. The mechanism and regularity of this influence are studied in this paper, where different incidence directions and angles of SV-waves have been considered, as shown in Figure 2.
In summary, the influence of subsurface geometry variations along a tunnel’s path on the seismic response and damage characteristics of said tunnel are the main focus of this paper. It is hereby noted that the original and practical subsurface geometry was not completely reproduced in the numerical model in this paper, because the local elevation and site information was classified, and, thus, the authors could not attain or publish the corresponding detailed information at present. Therefore, simplified models were established instead.

3. Establishing the Numerical Model

3.1. General Model Information

Figure 3 and Figure 4 demonstrate the dimensions of the one-sided mountain–plain site tunnel model and the corresponding mesh, respectively. The one-sided mountain–plain site tunnel model was 210 m, 100 m, and 50 m in length, width, and depth, respectively, while the mountain was 90 m, 40 m, and 50 m in length, width, and height. A circular underground tunnel was set in the model, and the length, outer radius, lining thickness, and burial depth (the distance between the site surface and the tunnel center) of the tunnel were 210 m, 2.5 m, 0.3 m, and 20 m, respectively. In terms of mesh, the C3D8R element was used to discretize the model, and the maximum size of the mesh for the whole model should be set to 1/6~1/12 of the wavelength in soft soil [32]. Herein, eight elements for the wavelength were assumed and employed (i.e., lmaxVs/8fmax, where lmax is the maximum element length; Vs denotes the shear wave velocity; and fmax is the maximum accurate frequency, set to 25 Hz in this study). The element size was set to 1.5 m for the whole site model, which can assure the adequate transmission of a wave in soil and rock. Although this size is much too small for rock material and lowers the calculation efficiency, it can make sure that seismic waves can propagate in a rock medium more adequately and brings about no influence on the results. The element number of the one-sided mountain–plain site tunnel model was 371340. Moreover, the maximum time increment was set to 0.00005 s. As for the tunnel lining, the C3D8R element was also used, and the element size was set to 1.5 m as well so that the nodes of the soil overlapped those of the lining. The coupled Lagrangian–Eulerian method was employed to describe the water–tunnel lining interaction and the sloshing effect of water under seismic waves [33]. The Eulerian element was utilized to discrete the water domain. In this model, the normal direction contact passes via the contact compressive stress transfer interaction, and the element nodes on the surfaces satisfy Hook’s Law and the Harmonized Condition of Displacement. Tangential contact obeys Coulomb’s law of friction. To obtain the ultimate stress to satisfy Coulomb’s law, the following formula was used:
τ crit = μ P ,
where μ and P denote the coefficient of friction and the normal contact pressure, respectively. Herein, the soil–tunnel lining interface friction coefficient was set to 0.4 [33,34].
In this study, the SV-wave was considered, and two sets of numerical simulations were designed, each one consisting of four cases with different incidence angles, as shown in Table 1. The incidence angles that were taken into consideration were 11.31°, 18.43°, 26.57°, and 33.69°, which corresponded to following vector directions: (x, y) = (1, 5), (1, 3), (1, 2), and (1, 1.5), respectively. Two different incidence directions were set, i.e., left-side incidence and right-side incidence. In the cases with left-side incidence, the wave passed through the rock mountain first and then the tunnel (see Figure 2a), while the wave passed through the tunnel first and then the rock mountain in the cases with right-side incidence (see Figure 2b).

3.2. Constitutive Model

The modified Davidenkov viscoelastic dynamic constitutive model was used to describe the dynamic characteristics of soil and rock [35], through the following expressions:
τ τ c = G max γ γ c 1 H γ γ c 2 n a ,
2 n a γ 0 2 B = γ ex ± γ c 2 B 1 R R ,
R = 1 τ ex ± τ c G max ( γ ex ± γ c ) 1 A ,
where τ and γ are the shear stress and the strain, respectively; Gmax is the initial shear modulus; A, B, and γ0 are the dimensionless fitting parameters; na is the coefficient controlling the scale of hysteresis loop; τc and γc denote the shear stress and the strain at the last stress reversal point, respectively; and the τex and γex are the shear stress and strain at the last extreme value point, respectively. The parameters of this model for soil and rock are shown in Table 2.
The tunnel lining’s concrete was assumed to be C30, and, to describe the compressive and tensile damage and degrading of the concrete, the concrete damage–plasticity constitutive model proposed by Jeeho and Fenves [36] was adopted. The corresponding material property parameters for the tunnel’s concrete lining are shown in Table 3. The physico-mechanical parameters of the fluid in the present analysis are described in Table 4.

3.3. Input Ground Motion and Boundary Conditions

Since most ground motion records have a long duration, the artificial simulation of ground motion was used in this paper to generate ground motion with a short duration and reduce the time cost of the calculation. The classical spectral representation method [37] was adopted to simulate the input ground motion:
f j ( t ) = A ( t ) l = 1 N 2 S 0 ( ω l ) Δ ω cos ( ω l t + φ l )
with S 0 ( ω l ) being the auto-power spectral density (PSD) and Δ ω denoting the bandwidth and N the number of frequency intervals. ω l is the lth value of frequency ω and can be calculated using ω l = l Δ ω , and Δω can be calculated using Δ ω = ω u / N , where ω u represents an upper cut-off frequency beyond which the elements of auto-PSD may be assumed to be zero for any time instant. Herein, ω u = 25 Hz, and N = 2048. φ l is a random phase angle uniformly distributed in [0, 2π], and the relevant probability density function can be expressed as f ( x ) = 1 ( 2 π 0 ) . Moreover, the following intensity envelope function was employed to proceed with the simulation of temporal non-stationary ground motion:
A ( t ) = ( t / t 0 ) 2 1 exp ( 0.155 ( t t n ) ) 0 t t 0 t 0 < t t n t n < t T
with t 0 = 2 s , and t n = 4 s .
The Clough–Penzien model [38] was utilized to establish the PSD mentioned above, with the following expression:
S 0 ( ω ) = S 0 1 + 4 ζ g 2 ω ω g 2 1 ω ω g 2 2 + 4 ζ g 2 ω ω g 2 × ω ω f 4 1 ω ω f 2 2 + 4 ζ f 2 ω ω f 2 .
with ω g and ζ g denoting the central frequency and damping ratio of the Tajimi–Kanai PSD function, and ω f and ζ f being the central frequency and damping ratio of the high-pass filter. And, according to Ref. [39],
ω g = 10 π ; ζ g = 0.6 ; ω f = 0.5 π ; ζ f = ζ g ;
S 0 = σ 2 π ω g 2 ζ g + 1 2 ζ g ; σ = 1.5   m / s 2
The duration of the simulated ground motions was set to 20 s and the time interval to 0.01 s. The simulated ground motion is shown in Figure 5. It should be noted that the peak value of the simulated ground motion was about 2.58 m/s2, and, clearly, this ground motion was a strong earthquake motion.
The viscous–elastic artificial boundary [7,33,40] composed of linear springs and dampers was used on the lateral and bottom sides of the site model. The velocity and displacement time-histories shown in Figure 5 were employed to derive the equivalent nodal force on the model boundary nodes.

4. Results and Analysis

4.1. Analysis of Wave Field

Figure 6 presents the diagrams of the wave field in the mountain–plain site for the SV-wave. Based on Snell’s law, the incidence angle, and the wave velocities in soil and rock, the angles of the transmitted SV- and P-waves at the soil–rock interface for the left-side-incident SV-wave are shown in Table 4, and the angles of the reflected waves for right-side incidence are in Table 5. Under the considered incidence angles, total reflection would occur on the left or right interface between soil and rock for left- or right-side-incident SV-waves. Since the right side of the mountain was inclined, total reflection was more likely to show up on the right side of the mountain under the right-side incidence of the wave. In the subsequent text, if lacking a special statement, the interface stands for the soil–rock interface on the right side of the mountain. Moreover, based on Snell’s law, the incidence angles considered herein, and the wave velocities in soil and rock, the critical angle that reflected P-wave can show up is 35.26° on the interface; however, the angle between the incident SV-wave and the normal direction of the interface was beyond the above-mentioned critical angle, and, thus, no reflected P-wave could be produced on the interface. Additionally, for the left-side-incident SV-wave, the tunnel was subjected to the incident wave, the transmitted converted waves coming from the rock mountain side, and the reflected converted waves from the site’s surface. For the right-side-incident SV-wave, the tunnel was mainly subjected to the incident wave, the reflected SV-wave on the soil–rock interface after total reflection, the reflected converted waves coming from the site’s surface, and the subsequent reflected waves. Due to the obstacle induced by total reflection, part of the waves could not go into the rock mountain from the soil, and, thus, multiple reflections would occur in the soil for right-side incidence. Taking all the above-mentioned information into consideration and based on Figure 6, it can be found that, for both left-side incidence and right-side incidence, the tunnel’s lining may be subjected to multiple waves pointing to the tunnel, and an inclined subsurface geometry focuses both left-side and right-side obliquely incident waves.
Figure 7 displays the displacement contour of the site for the SV-wave with an incidence angle = 33.69°, where the left, bottom region of Figure 7b is enlarged to show the wave’s front. Note that the observation time instant during which the incident wave just passes through the rock mountain has been chosen to observe the propagation direction of the transmitted wave. Firstly, it can be seen that the displacement contour is obliquely layered, and this demonstrates that the waves are, indeed, obliquely incident to the site. Secondly, as shown in the enlarged zone in Figure 7b, as the wave propagates, the inclined angle of the wave’s front, near the soil–rock interface on the left side of the mountain, increases gradually, and this might suggest that, when arriving at the interface between soil and rock, waves encounter some obstacles. The wave’s front tends to up-warp at the soil–rock interface in comparison to the initial incidence, which suggests that total reflection indeed occurs. Thirdly, the propagation direction in the rock mountain and that in the soil site behind the mountain are noticeably different due to refraction. Based on Snell’s law, the shear wave velocity in soil and rock, and the incidence angle, the transmission angle is equal to 13° for the SV-wave transmitting from the rock to the soil. Considering that the slope angle of the right side of the mountain is 80°, the angle between the transmitted wave coming from the rock mountain to the soil site and the horizon line should be 23°. Then, the propagation direction of the transmitted wave (solid line with a single arrow) derived from the theoretical calculation is given in Figure 7b and compared with the wave’s front in the displacement contour obtained by numerical simulation. It can be seen that the wave’s front of the transmitted wave in the numerical results is just perpendicular to the theoretical propagation direction. In the meantime, the predefined direction of the incident wave (solid line with double arrows) is also presented in Figure 7, and it can be found that the numerical results are in good agreement with the theoretical propagation direction. Taking all the above-mentioned information into account, the analysis of the wave field shown in Figure 6 can be verified, and the validity of this study can also be verified to some extent.

4.2. Acceleration Response of Tunnel Lining

The acceleration responses of the tunnel under SV-waves with different incidence directions and angles were calculated and analyzed. Note that, since the longitudinal seismic response of the tunnel is the object of concern, the longitudinal envelope of the acceleration response is presented and analyzed herein. As a matter of fact, the seismic responses at all the outer element nodes at the arch vault, arch bottom, left arch, and right arch along the tunnel model were extracted, and then the envelopes of the seismic response for the arch vault, arch bottom, left arch, and right arch were obtained and are displayed in Figure 8, Figure 9, Figure 10 and Figure 11, respectively. According to Table 5 and Table 6, the propagation and vibration directions of the various waves in the site for incidence angles = 11.31° and 33.69° are drawn in Figure 12 and Figure 13, respectively, to facilitate a relevant analysis.
Figure 8, Figure 9, Figure 10 and Figure 11 show that the acceleration response at the arch vault, arch bottom, left arch, and right arch of the tunnel’s lining in the vertical direction are ultimately lower than those in the horizontal direction. This can be attributed to the fact that the vibration direction of the SV-wave is perpendicular to the propagation direction, and the incidence angles taken into consideration are lower than 45°, and, thus, the vertical component of motion is smaller. As the incidence angle increases, the horizontal component of the acceleration response of the lining decreases, while the vertical component increases gradually, because, as the incidence angle increases, the horizontal component of the incident SV-wave decreases, whereas the vertical component increases. The above-mentioned results are in agreement with theoretical knowledge, which verifies the validity of this study to some extent.
Under either left incidence or right incidence, the horizontal and vertical acceleration responses of the 60~150 m tunnel interval (i.e., the tunnel interval near the rock mountain) are noticeably higher than those of other intervals. In other words, the acceleration response of the tunnel in the plain–mountain site is greater than that in the plain site. As stated in Section 4.1, the wave field of the soil site located on the right side of the mountain for right incidence is remarkably complex due to the total reflection at the soil–rock interface and surface reflection, and this results in the fact that the tunnel’s lining may be subjected to multiple waves, which aggregate the vibration of the lining. Similarly, for left incidence, the transmitted SV- and P-waves produced by the left-side-incident SV-wave at the soil–rock interface possess different propagation and vibration directions from the incident wave that comes from the boundary beneath the tunnel (see Figure 12 and Figure 13), and, thus, the superposition of these waves may aggregate the vibration of the tunnel’s lining. However, in the plain-soil site, the site conditions are uniform, and no total reflection or refraction occur. Therefore, compared to the tunnel in the plain site, the tunnel interval adjacent to the rock mountain may be subjected to multiple waves possessing different directions of propagation and vibration due to obliquely incident wave and non-uniform subsurface geometry conditions, and, thus, the acceleration response of the tunnel in the 60~150 m range is higher.
Moreover, the horizontal and vertical acceleration responses of the arch vault, arch bottom, left arch, and right arch of the 60~150 m tunnel interval under left-side incidence are lower than those under right-side incidence. This is because multiple reflections occur at the site’s surface and soil–rock interface in the cases with right-side incidence, and, thus, multiple conversions of waves may occur to a great extent. Therefore, the complex superposition of waves induces the aggravation of the lining’s acceleration response. Regarding left-side incidence, the propagation direction of the transmitted and reflected waves generally points to the lateral model boundary, and, thus, it is easier for the waves to run outside of the site, meaning that the degree of multiple reflections is lower.

4.3. Damage Assessment of Tunnel Lining

The tensile damage contours in the longitudinal view of the tunnel lining under right-side-incident SV-waves for the different incidence angles considered are shown in Figure 14, where the selected time intervals are (i) the moment when the damage first appeared and (ii) the end time of the total duration of the input wave. As it can be seen in the longitudinal dimension, the position in which the damage first showed up is located in the tunnel interval of 60~150 m. This is because the existence of a local rock mountain brings about the reflection of waves (see Figure 12 and Figure 13), and the superposition of the reflected wave and the incident wave aggravates the deformation of the lining and then damages it, whereas, in a plain site, the above-mentioned reflection and superposition of waves are absent. In the circumferential dimension, the positions are different among the cases: as the incidence angle decreases, the position in which the damage first showed up tends to move from the vault towards the hance. Since the damage of the tunnel’s lining is highly related to deformation, the deformation of the lining under right-side-incident SV-waves with different incidence angles at the moment when the damage first appeared was observed, as shown in Figure 15, where the scaling factor is set to 300 so that the deformation of the lining can be investigated more clearly and two profiles are selected for each case therein: (i) the middle of the plain site with a longitudinal distance = 30 m; and (ii) the middle of the plain–mountain site with a longitudinal distance = 105 m. As it can be seen, thanks to the 300 scaling factor used, the ovaling deformation of the lining seems to be significant, and the compressive deformation is approximately perpendicular to the incidence direction of the SV-wave. To better observe the relation between the ovaling deformation of the tunnel lining and the incidence angles, the major axis (dash line) of the ovaling lining in the profile for the middle of the plain–mountain site with a longitudinal distance = 105 m is given for each case. Clearly, the inclined angle of the major axis tends to be greater with a decreasing incidence angle, meaning that the part of lining which is initially subjected to severe tension, i.e., the position where tensile damage first appears, gradually changes from the vault to the hance with a decreasing incidence angle.
In addition, for each case with right-side incidence, the final damage area reduces as the incidence angle decreases. As shown in Figure 13b, when the incidence angle = 33.69°, the tunnel lining is subjected to a directly incident SV-wave and reflected SV-waves from the soil–rock interface and the site’s surface, and the difference in the angle between the incident and reflected waves essentially induces the tunnel lining to experience compressive deformation in two directions. The fact that tensile damage shows up at the arch vault, arch bottom, left arch, and right arch of the 60~150 m tunnel interval can well verify the above-mentioned statement. By contrast, when the incidence angle = 11.31° (see Figure 12b), the difference in the angle between the incident and reflected waves is much smaller, and, thus, the tunnel lining is approximately subjected to a unidirectional deformation, and, thus, tensile damage only appears in a single, diagonal direction. In summary, as the incidence angle decreases, the mode of compressive deformation changes from multi-directional compression to unidirectional compression. To verify this statement, the deformation of tunnel lining in the profile for longitudinal distance = 105 m at time interval = 6.0 s under right-side-incident SV-wave was observed, as shown in Figure 16. It can be obviously found that at the same time interval the deformation modes of lining under right-side-incident SV-wave with predefined various incidence angles are significantly different, and when the incidence angle = 33.69°, the shape of deformed lining tends to be rectangular. This demonstrates that the lining may be subjected to multi-directional deformation.
The tensile damage contours in the longitudinal view of the tunnel’s lining under left-side-incident SV-waves for the different incidence angles considered and different time intervals are shown in Figure 17. The deformation of the lining under left-side-incident SV-waves with different incidence angles at the moment when damage first appears is shown in Figure 18. Similar to the cases with right-side incidence, the part of lining that is initially subjected to severe tension gradually changes from the vault to the hance with varying incidence angles in the cases with left-side incidence. However, what is different from the right-side cases is that, in all four cases of left-side incidence, the tunnel lining displays damage in multiple directions. As shown in Table 5, the transmitted angle of the SV- and P-waves at the soil–rock interface varies slightly with varying incidence angles, meaning that the propagation direction of the transmitted SV- and P-waves is always significantly different from that of the incident waves (see Figure 12 and Figure 13). Therefore, for all the incidence angles considered, the tunnel lining always experiences compression in multiple directions. Thus, the damage distribution of the lining for the cases with left incidence is always multi-directional.

5. Conclusions

Considering longitudinal non-uniform subsurface geometry conditions and different incidence directions and arbitrary angles of the SV-wave, the longitudinal seismic response and damage characteristics of a water-conveyance tunnel were investigated in this paper using the numerical method. The main conclusions can be summarized as follows: (i) an inclined subsurface geometry may focus both left-side and right-side obliquely incident waves due to refraction or total reflection; (ii) a tunnel in a site adjacent to a rock mountain may exhibit higher horizontal and vertical acceleration responses than a tunnel in a homogeneous plain site due to the refraction or reflection induced by the rock mountain; and (iii) a tunnel in a site adjacent to a rock mountain may be subjected to compression in multiple directions, thus damage may appear in multiple circumferential directions, and, as the incidence angle decreases, the mode of tensile damage for right-side incidence may change from multi-directional damage to unidirectional damage, while that of tensile damage for left-side incidence may always be multi-directional.
The chosen incidence angles in this study were all smaller than the critical incidence angle of the SV-wave, and, thus, the horizontal component of the input wave and the seismic response were higher than those in the vertical direction. The calculation and analysis results of cases in which the vertical component is higher may be needed in the future.

Author Contributions

Conceptualization, A.C.; methodology, E.Y.; software, Y.R.; validation, Y.R.; investigation, M.L.; resources, A.C.; data curation, E.Y.; writing—original draft preparation, E.Y.; writing—review and editing, Z.L.; visualization, M.L.; supervision, E.Y.; project administration, A.C.; funding acquisition, E.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the Knowledge Innovation Program of Wuhan-Shuguang (No. 2022010801020242), the National Natural Science Foundation of China (Nos. 51908236, 51978304, 51778260, and 52079009), and the China Postdoctoral Science Foundation (2019M652649).

Data Availability Statement

Data is contained within the article.

Acknowledgments

The authors also thank Yu Miao from the Huazhong University of Science and Technology and Gen Zhao from the Changjiang River Scientific Research Institute for their help. The first author is grateful for his wife, Yan Yan, for the support and care.

Conflicts of Interest

Author Ang Cao was employed by the company China Communications Construction Company Second Harbor Engineering Company Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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. Shen, Y.; Gao, B.; Yang, X.; Tao, S. Seismic damage mechanism and dynamic deformation characteristic analysis of mountain tunnel after Wenchuan earthquake. Eng. Geol. 2014, 180, 85–98. [Google Scholar] [CrossRef]
  3. Li, T. Damage to mountain tunnels related to the Wenchuan earthquake and some suggestions for a seismic tunnel construction. Bull. Eng. Geol. Environ. 2012, 71, 297–308. [Google Scholar] [CrossRef]
  4. Hashash, Y.M.; Hook, J.J.; Schmidt, B.; John, I.; Yao, C. Seismic design and analysis of underground structures. Tunn. Undergr. Space Technol. 2001, 16, 247–293. [Google Scholar] [CrossRef]
  5. Stamos, A.A.; Beskos, D.E. Dynamic analysis of large 3-D underground structures by the BEM. Earthq. Eng. Struct. Dyn. 1995, 24, 917–934. [Google Scholar] [CrossRef]
  6. Manolis, G.D.; Dineva, P.S. Elastic waves in continuous and discontinuous geological media by boundary integral equation methods: A review. Soil. Dyn. Earthq. Eng. 2015, 70, 11–29. [Google Scholar] [CrossRef]
  7. Huang, J.-Q.; Du, X.-L.; Zhao, M.; Zhao, X. Impact of incident angles of earthquake shear (S) waves on 3-D non-linear seismic responses of longlined tunnels. Eng. Geol. 2016, 222, 168–185. [Google Scholar] [CrossRef]
  8. de Barros, F.C.P.; Luco, J.E. Diffraction of obliquely incident waves by a cylindrical cavity embedded in a layered viscoelastic half-space. Soil. Dyn. Earthq. Eng. 1993, 12, 159–171. [Google Scholar] [CrossRef]
  9. Huang, J.; Du, X.; Jin, L.; Zhao, M. Impact of incident angles of P waves on the dynamic responses of long lined tunnels. Earthq. Eng. Struct. Dyn. 2016, 45, 2435–2454. [Google Scholar] [CrossRef]
  10. Stamos, A.A.; Beskos, D.E. 3-D seismic response analysis of long lined tunnels in halfspace. Soil. Dyn. Earthq. Eng. 1996, 15, 111–118. [Google Scholar] [CrossRef]
  11. Pitilakis, K.; Tsinidis, G. Performance and seismic design of underground structures. In Earthquake Geotechnical Engineering Design; Maugeri, M., Soccodato, C., Eds.; Springer: Berlin/Heidelberg, Germany, 2014; pp. 279–340. [Google Scholar] [CrossRef]
  12. Kiureghian, A.D. A coherency model for spatially varying ground motions. Earthq. Eng. Struct. Dyn. 1996, 19, 533–547. [Google Scholar] [CrossRef]
  13. Yan, L.Q.; Abbas, H.; Li, P.; Song, E.X. A numerical study on the transverse seismic response of lined circular tunnels under obliquely incident asynchronous P and SV waves. Tunn. Undergr. Space Technol. 2020, 97, 103235. [Google Scholar] [CrossRef]
  14. Li, P.; Song, E.X. Three-dimensional numerical analysis for the longitudinal seismic response of tunnels under an asynchronous wave input. Comput. Geotech. 2015, 63, 229–243. [Google Scholar] [CrossRef]
  15. Zhu, J.; Li, X.J.; Liang, J.W. 2.5D FE-BE modelling of dynamic responses of segmented tunnels subjected to obliquely incident seismic waves. Soil. Dyn. Earthq. Eng. 2022, 163, 107564. [Google Scholar] [CrossRef]
  16. Zhu, J.; Liang, J.; Ba, Z. A 2.5D equivalent linear model for longitudinal seismic analysis of tunnels in water-saturated poroelastic half-space. Comput. Geotech. 2019, 109, 166–188. [Google Scholar] [CrossRef]
  17. Zhou, X.; Liang, Q.; Zhang, Y.; Liu, Z.; He, Y. Three-Dimensional Nonlinear Seismic Response of Immersed Tunnel in Horizontally Layered Site under Obliquely Incident SV Waves. Shock Vib. 2019, 2019, 1–17. [Google Scholar] [CrossRef]
  18. Zhu, J.; Li, X.J.; Liang, J.W.; Kontoe, S.; He, Q.; Chen, S. Effects of pulse-like ground motions on tunnels in saturated poroelastic soil for obliquely incident seismic waves. Soil Dyn. Earthq. Eng. 2023, 173, 107564. [Google Scholar] [CrossRef]
  19. Yang, Y.-B.; Wang, X.; Zhou, Z.; Zhang, X.; Li, Z. Seismic analysis of underground tunnels subjected to 3D oblique incident P and SV waves by 2.5D approach. Undergr. Space 2023, 12, 271–286. [Google Scholar] [CrossRef]
  20. Huang, Z.Y.; Feng, Y.Z.; Tang, A.P. Influence of oblique incidence of P-waves on seismic response of prefabricated utility tunnels considering joints. Soil. Dyn. Earthq. Eng. 2023, 167, 107797. [Google Scholar] [CrossRef]
  21. Huang, L.; Liu, Z.X.; Wu, C.; Liang, J. The scattering of plane P, SV waves by twin lining tunnels with imperfect interfaces embedded in an elastic half-space. Tunn. Undergr. Space Technol. 2019, 85, 319–330. [Google Scholar] [CrossRef]
  22. Zhang, G.L.; Wang, P.; Zhao, M.; Du, X.; Zhao, X. Seismic structure-water-sediment-rock interaction model and its application to immersed tunnel analysis under obliquely incident earthquake. Tunn. Undergr. Space Technol. 2021, 109, 103758. [Google Scholar] [CrossRef]
  23. Yasuda, N.; Asakura, T. Seismic response of a cylindrical tunnel with construction joints subjected to longitudinal ground displacement. Tunn. Undergr. Space Technol. 2020, 102, 103408. [Google Scholar] [CrossRef]
  24. Park, D.; Sagong, M.; Kwak, D.Y.; Jeong, C.G. Simulation of tunnel response underspatially varying ground motion. Soil Dyn. Earthq. Eng. 2009, 29, 1417–1424. [Google Scholar] [CrossRef]
  25. Miao, Y.; Yao, E.; Ruan, B.; Zhuang, H. Seismic response of shield tunnel subjected to spatially varying earthquake ground motions. Tunn. Undergr. Space Technol. 2018, 77, 216–226. [Google Scholar] [CrossRef]
  26. Miao, Y.; Yao, E.; Ruan, B.; Zhuang, H.; Chen, G.; Long, X. Improved hilbert spectral representation method and its application to seismic analysis of shield tunnel subjected to spatially correlated ground motions. Soil Dyn. Earthq. Eng. 2018, 111, 119–130. [Google Scholar] [CrossRef]
  27. Liu, L.; Xu, C.; Du, X.L.; Iqbal, K. Longitudinal seismic response of shield tunnel: A multi-scale numerical analysis. Tunn. Undergr. Space Technol. 2023, 138, 105163. [Google Scholar] [CrossRef]
  28. Yuan, Y.; Yu, H.T.; Li, C.; Yan, X.; Yuan, J.T. Multi-point shaking table test for long tunnel ssubjected to non-uniform seismic loadings–Part I: Theory and validation. Soil Dynam Earthq. Eng. 2018, 108, 177–186. [Google Scholar] [CrossRef]
  29. Yu, H.T.; Yuan, Y.; Xu, G.; Su, Q.; Yan, X.; Li, C. Multi-point shaking table test for long tunnels subjected to non-uniform seismic loadings—Part II: Application to the HZM immersed tunnel. Soil Dyn. Earthq. Eng. 2018, 108, 187–195. [Google Scholar] [CrossRef]
  30. Yan, X.; Yu, H.; Yuan, Y.; Yuan, J. Multi-point shaking table test of the free field under non-uniform earthquake excitation. Soils Found. 2015, 55, 985–1000. [Google Scholar] [CrossRef]
  31. Zhang, N. 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]
  32. Karabalis, D.L.; Beskos, D.E. Numerical methods in earthquake engineering. In Computer Analysis and Design of Earthquake Resistant Structures: A Handbook; Beskos, D.E., Anagnostopoulos, S.A., Eds.; Computational Mechanics Publications: Southampton, UK, 1997; pp. 1–102. [Google Scholar]
  33. Sun, B.B.; Deng, M.; Zhang, S.; Cui, W.; Wang, C.; Yu, L.; Cao, K. Inelastic dynamic analysis and damage assessment of a hydraulic arched tunnel under near-fault SV waves with arbitrary incoming angles. Tunn. Undergr. Space Technol. 2020, 104, 103523. [Google Scholar] [CrossRef]
  34. Zhuang, H.Y.; Ren, J.; Miao, Y.; Jing, L.; Yao, E.; Xu, C. Seismic Performance Levels of a Large Underground Subway Station in Different Soil Foundations. J. Earthq. Eng. 2021, 25, 2808–2833. [Google Scholar] [CrossRef]
  35. Zhao, D.F.; Ruan, B.; Chen, G.X. Validation of the modified irregular loading-reloading rules based on Davidenkov skeleton curve and its equivalent strain algorithm implemented in ABAQUS. Chin. J. Geotech. Eng. 2017, 39, 888–895. (In Chinese) [Google Scholar]
  36. Jeeho, L.; Fenves, G.L. Plastic-damage model for cyclic loading of concrete structures. J. Eng. Mech. 1998, 124, 892–900. [Google Scholar]
  37. Shinozuka, M.; Jan, C.M. Digital simulation of random processes and its applications. J. Sound Vib. 1972, 25, 111–128. [Google Scholar] [CrossRef]
  38. Clough, R.W.; Penzien, J. Dynamics of Structures; McGraw Hill: New York, NY, USA, 1993. [Google Scholar]
  39. Bi, K.M.; Hao, H. Modelling and simulation of spatially varying earthquake ground motions at sites with varying conditions. Probabilistic Eng. Mech. 2012, 29, 92–104. [Google Scholar] [CrossRef]
  40. Du, X.; Zhao, M.; Wang, J. A stress artificial boundary in FEA for near-field wave Problem. Chin. J. Theory Appl. Mech. 2006, 38, 49–56. [Google Scholar]
Figure 1. The spatial location relationship between the tunnel and the mountain for the Kunming interval of the Dianzhong Water Diversion Project.
Figure 1. The spatial location relationship between the tunnel and the mountain for the Kunming interval of the Dianzhong Water Diversion Project.
Applsci 14 04398 g001
Figure 2. Schematic diagram for the local one-sided mountain–plain site tunnel model: (a) left side incidence; and (b) right side incidence.
Figure 2. Schematic diagram for the local one-sided mountain–plain site tunnel model: (a) left side incidence; and (b) right side incidence.
Applsci 14 04398 g002
Figure 3. Dimensions of the one-sided mountain tunnel model.
Figure 3. Dimensions of the one-sided mountain tunnel model.
Applsci 14 04398 g003
Figure 4. Mesh of the one-sided mountain tunnel model.
Figure 4. Mesh of the one-sided mountain tunnel model.
Applsci 14 04398 g004
Figure 5. Simulated ground motions: (a) acceleration time–history; (b) velocity time–history; (c) displacement time–history; and (d) the PSD of the simulated ground motion.
Figure 5. Simulated ground motions: (a) acceleration time–history; (b) velocity time–history; (c) displacement time–history; and (d) the PSD of the simulated ground motion.
Applsci 14 04398 g005
Figure 6. Wave propagation diagram for SV-waves: (a) left-side incidence; and (b) right-side incidence.
Figure 6. Wave propagation diagram for SV-waves: (a) left-side incidence; and (b) right-side incidence.
Applsci 14 04398 g006
Figure 7. Displacement contour under incident SV–wave for (a) plain site and left-side incidence at t = 0.27 s; (b) plain–mountain site and left–side incidence at t = 0.27 s; (c) plain site and right–side incidence at t = 0.36 s; and (d) plain–mountain site and right-side incidence at t = 0.36 s. (Unit: m).
Figure 7. Displacement contour under incident SV–wave for (a) plain site and left-side incidence at t = 0.27 s; (b) plain–mountain site and left–side incidence at t = 0.27 s; (c) plain site and right–side incidence at t = 0.36 s; and (d) plain–mountain site and right-side incidence at t = 0.36 s. (Unit: m).
Applsci 14 04398 g007aApplsci 14 04398 g007b
Figure 8. Seismic response at the arch vault of the tunnel: (a) acceleration response in the horizontal direction; and (b) acceleration response in the vertical direction.
Figure 8. Seismic response at the arch vault of the tunnel: (a) acceleration response in the horizontal direction; and (b) acceleration response in the vertical direction.
Applsci 14 04398 g008
Figure 9. Seismic response at the arch bottom of the tunnel: (a) acceleration response in the horizontal direction; and (b) acceleration response in the vertical direction.
Figure 9. Seismic response at the arch bottom of the tunnel: (a) acceleration response in the horizontal direction; and (b) acceleration response in the vertical direction.
Applsci 14 04398 g009
Figure 10. Seismic response at the left arch of the tunnel: (a) acceleration response in the horizontal direction; and (b) acceleration response in the vertical direction.
Figure 10. Seismic response at the left arch of the tunnel: (a) acceleration response in the horizontal direction; and (b) acceleration response in the vertical direction.
Applsci 14 04398 g010
Figure 11. Seismic response at the right arch of the tunnel: (a) acceleration response in the horizontal direction; and (b) acceleration response in the vertical direction.
Figure 11. Seismic response at the right arch of the tunnel: (a) acceleration response in the horizontal direction; and (b) acceleration response in the vertical direction.
Applsci 14 04398 g011
Figure 12. Diagram of the wave field for (a) left-side-incident SV-wave and (b) right-side-incident SV-wave with an incidence angle = 11.31° in the plain–mountain site.
Figure 12. Diagram of the wave field for (a) left-side-incident SV-wave and (b) right-side-incident SV-wave with an incidence angle = 11.31° in the plain–mountain site.
Applsci 14 04398 g012
Figure 13. Diagram of the wave field for (a) left-side-incident SV-wave and (b) right-side-incident SV-wave with an incidence angle = 33.69° in the plain–mountain site.
Figure 13. Diagram of the wave field for (a) left-side-incident SV-wave and (b) right-side-incident SV-wave with an incidence angle = 33.69° in the plain–mountain site.
Applsci 14 04398 g013
Figure 14. Tensile damage of tunnel lining under right-side-incident SV-wave: (a) incidence angle = 33.69° and t = 2.25 s; (b) incidence angle = 33.69° and t = 20 s; (c) incidence angle = 26.57° and t = 2.31 s; (d) incidence angle = 26.57° and t = 20 s; (e) incidence angle = 18.43° and t = 2.26 s; (f) incidence angle = 18.43° and t = 20 s; (g) incidence angle = 11.31° and t = 2.21 s; and (h) incidence angle = 11.31° and t = 20 s.
Figure 14. Tensile damage of tunnel lining under right-side-incident SV-wave: (a) incidence angle = 33.69° and t = 2.25 s; (b) incidence angle = 33.69° and t = 20 s; (c) incidence angle = 26.57° and t = 2.31 s; (d) incidence angle = 26.57° and t = 20 s; (e) incidence angle = 18.43° and t = 2.26 s; (f) incidence angle = 18.43° and t = 20 s; (g) incidence angle = 11.31° and t = 2.21 s; and (h) incidence angle = 11.31° and t = 20 s.
Applsci 14 04398 g014
Figure 15. Tunnel lining deformation at the time when damage initially shows up for right-side-incident SV-wave with the following incidence angles: (a) 33.69°; (b) 26.57°; (c) 18.43°; and (d) 11.31° (scaling factor = 300; unit: m).
Figure 15. Tunnel lining deformation at the time when damage initially shows up for right-side-incident SV-wave with the following incidence angles: (a) 33.69°; (b) 26.57°; (c) 18.43°; and (d) 11.31° (scaling factor = 300; unit: m).
Applsci 14 04398 g015
Figure 16. Tunnel lining deformation at t = 6.0 s for right-side-incident SV-wave with the following incidence angles: (a) 33.69°; (b) 26.57°; (c) 18.43°; and (d) 11.31° (scaling factor = 300; unit: m).
Figure 16. Tunnel lining deformation at t = 6.0 s for right-side-incident SV-wave with the following incidence angles: (a) 33.69°; (b) 26.57°; (c) 18.43°; and (d) 11.31° (scaling factor = 300; unit: m).
Applsci 14 04398 g016
Figure 17. Tensile damage of tunnel lining under left-side-incident SV-wave: (a) incidence angle = 33.69° and t = 2.25 s; (b) incidence angle = 33.69° and t = 20 s; (c) incidence angle = 26.57° and t = 2.31 s; (d) incidence angle = 26.57° and t = 20 s; (e) incidence angle = 18.43° and t = 2.26 s; (f) incidence angle = 18.43° and t = 20 s; (g) incidence angle = 11.31° and t = 2.21 s; and (h) incidence angle = 11.31° and t = 20 s.
Figure 17. Tensile damage of tunnel lining under left-side-incident SV-wave: (a) incidence angle = 33.69° and t = 2.25 s; (b) incidence angle = 33.69° and t = 20 s; (c) incidence angle = 26.57° and t = 2.31 s; (d) incidence angle = 26.57° and t = 20 s; (e) incidence angle = 18.43° and t = 2.26 s; (f) incidence angle = 18.43° and t = 20 s; (g) incidence angle = 11.31° and t = 2.21 s; and (h) incidence angle = 11.31° and t = 20 s.
Applsci 14 04398 g017
Figure 18. Tunnel lining deformation at the time when damage initially shows up for left-side-incident SV-wave with the following incidence angles: (a) 33.69°; (b) 26.57°; (c) 18.43°; and (d) 11.31° (scaling factor = 300; unit: m).
Figure 18. Tunnel lining deformation at the time when damage initially shows up for left-side-incident SV-wave with the following incidence angles: (a) 33.69°; (b) 26.57°; (c) 18.43°; and (d) 11.31° (scaling factor = 300; unit: m).
Applsci 14 04398 g018
Table 1. Sets of numerical simulation.
Table 1. Sets of numerical simulation.
No.ModelIncidence DirectionWave TypeIncidence Angles
Set-11One-sided mountain tunnel modelLeft-side incidenceSV-wave33.69°
226.57°
318.43°
411.31°
Set-25Right-side incidence33.69°
626.57°
718.43°
811.31°
Table 2. Values of the fitting parameters of the modified Davidenkov model.
Table 2. Values of the fitting parameters of the modified Davidenkov model.
Soil TypesDensity
(kg/m3)
Poisson’s RatioShear Wave Velocity
(m/s)
Compressive Wave Velocity
(m/s)
Initial Shear Modulus
(MPa)
Fitting Parameters
ABγ0 (×10−4)
Medium-coarse sand20000.25280.18485.281571.200.377.4
Bed rock25220.2899.821469.4020421.300.4010.0
Table 3. Material property of the tunnel’s concrete lining.
Table 3. Material property of the tunnel’s concrete lining.
ParameterDescriptionValue
ρ Density2450 kg/m3
E 0 Elastic modulus30,000 MPa
μ Poisson’s ratio0.18
Ψ Dilation angle36.31°
ε 0 Flow potential eccentricity0.1
σ b 0 / σ c 0 The ratio of initial equibiaxial compressive yield stress
to initial uniaxial compressive yield stress
1.16
KThe ratio of the second stress invariant on the tensile meridian0.6667
V c Viscosity coefficient0.0005
ω c Compression stiffness recovery parameter1.0
ω t Tensile stiffness recovery parameter0.0
Table 4. Physico-mechanical parameters of the water inner tunnel.
Table 4. Physico-mechanical parameters of the water inner tunnel.
MaterialDensity
(kg/m3)
Dynamic Viscosity Coefficient
(Ns/m3)
Equation of State Parameter
(m/s)
Water10000.00131400
Table 5. Angles of the transmitted SV- and P-waves at the soil–rock interface for a left-side-incident SV-wave.
Table 5. Angles of the transmitted SV- and P-waves at the soil–rock interface for a left-side-incident SV-wave.
Incidence Angle/°Angle between Incident Wave and Normal Direction of Soil–Rock Interface/°Transmitted Angle of SV-Wave/°Transmitted Angle of Converted P-Wave/°
33.6946.3113.0122.95
25.6754.3314.6525.98
18.4361.5715.8928.31
11.3168.6916.8630.16
Table 6. Angles of the reflected SV- and P-waves at the soil–rock interface for a right-side-incident SV-wave.
Table 6. Angles of the reflected SV- and P-waves at the soil–rock interface for a right-side-incident SV-wave.
Incidence Angle/°Angle between Incident Wave and Normal Direction of Soil–Rock Interface/°Reflected Angle of SV-Wave/°Reflected Angle of Converted P-Wave/°Transmitted Angle of SV- or P-Wave/°
33.6966.3166.31No reflected converted P-waveNo transmitted wave
25.6774.3374.33
18.4381.5781.57
11.3188.6988.69
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.

Share and Cite

MDPI and ACS Style

Yao, E.; Rao, Y.; Liu, M.; Liu, Z.; Cao, A. Study on the Seismic Response of a Water-Conveyance Tunnel Considering Non-Uniform Longitudinal Subsurface Geometry and Obliquely Incident SV-Waves. Appl. Sci. 2024, 14, 4398. https://doi.org/10.3390/app14114398

AMA Style

Yao E, Rao Y, Liu M, Liu Z, Cao A. Study on the Seismic Response of a Water-Conveyance Tunnel Considering Non-Uniform Longitudinal Subsurface Geometry and Obliquely Incident SV-Waves. Applied Sciences. 2024; 14(11):4398. https://doi.org/10.3390/app14114398

Chicago/Turabian Style

Yao, Erlei, Yu Rao, Meishan Liu, Zhifang Liu, and Ang Cao. 2024. "Study on the Seismic Response of a Water-Conveyance Tunnel Considering Non-Uniform Longitudinal Subsurface Geometry and Obliquely Incident SV-Waves" Applied Sciences 14, no. 11: 4398. https://doi.org/10.3390/app14114398

APA Style

Yao, E., Rao, Y., Liu, M., Liu, Z., & Cao, A. (2024). Study on the Seismic Response of a Water-Conveyance Tunnel Considering Non-Uniform Longitudinal Subsurface Geometry and Obliquely Incident SV-Waves. Applied Sciences, 14(11), 4398. https://doi.org/10.3390/app14114398

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