
A horizontal non-homogeneous field adversely affects the seismic resistance of both the utility tunnel and its internal pipes, with seismic waves obliquely incident on the underground structure causing more significant damages. To address these issues, this study, based on a viscous-spring artificial boundary, derives and validates the equivalent junction force formula for the horizontal non-homogeneous field. It then establishes a three-dimensional finite element model of the utility tunnel, pipes, and surrounding soil to obtain the acceleration and strain responses of the utility tunnel and its internal pipes under seismic loading. Finally, it investigates the impact of different incidence angles of shear waves (SV waves) on the response of the utility tunnel and its internal pipes. It was found that as the PGA increases from 0.1 to 0.4 g, both peak acceleration and strain of the utility tunnel and its internal pipes increase. The peak acceleration of the utility tunnel and pipes initially decreases and then increases with the angle of incidence, while the strain increases with the angle of incidence, reaching its peak value when the angle of incidence is 30°. The acceleration and strain responses of the utility tunnel and pipe are higher in sand than in clay, with the peak acceleration strongly correlating with the angle of incidence of ground shaking. The findings of this study provide valuable insights into the seismic design of horizontal non-homogeneous field utility tunnel systems.
Citation: Hui-yue Wang, Sha-sha Yu, De-long Huang, Chang-lu Xu, Hang Cen, Qiang Liu, Zhong-ling Zong, Zi-Yuan Huang. Seismic response of utility tunnel systems embedded in a horizontal heterogeneous domain subjected to oblique incident SV-wave[J]. AIMS Geosciences, 2025, 11(1): 47-67. doi: 10.3934/geosci.2025004
[1] | T. A. Shaposhnikova, M. N. Zubova . Homogenization problem for a parabolic variational inequality with constraints on subsets situated on the boundary of the domain. Networks and Heterogeneous Media, 2008, 3(3): 675-689. doi: 10.3934/nhm.2008.3.675 |
[2] | Ken-Ichi Nakamura, Toshiko Ogiwara . Periodically growing solutions in a class of strongly monotone semiflows. Networks and Heterogeneous Media, 2012, 7(4): 881-891. doi: 10.3934/nhm.2012.7.881 |
[3] | Junlong Chen, Yanbin Tang . Homogenization of nonlinear nonlocal diffusion equation with periodic and stationary structure. Networks and Heterogeneous Media, 2023, 18(3): 1118-1177. doi: 10.3934/nhm.2023049 |
[4] | Benjamin Contri . Fisher-KPP equations and applications to a model in medical sciences. Networks and Heterogeneous Media, 2018, 13(1): 119-153. doi: 10.3934/nhm.2018006 |
[5] | Thomas Geert de Jong, Georg Prokert, Alef Edou Sterk . Reaction–diffusion transport into core-shell geometry: Well-posedness and stability of stationary solutions. Networks and Heterogeneous Media, 2025, 20(1): 1-14. doi: 10.3934/nhm.2025001 |
[6] | Iryna Pankratova, Andrey Piatnitski . Homogenization of convection-diffusion equation in infinite cylinder. Networks and Heterogeneous Media, 2011, 6(1): 111-126. doi: 10.3934/nhm.2011.6.111 |
[7] | Feiyang Peng, Yanbin Tang . Inverse problem of determining diffusion matrix between different structures for time fractional diffusion equation. Networks and Heterogeneous Media, 2024, 19(1): 291-304. doi: 10.3934/nhm.2024013 |
[8] | Xavier Blanc, Claude Le Bris . Improving on computation of homogenized coefficients in the periodic and quasi-periodic settings. Networks and Heterogeneous Media, 2010, 5(1): 1-29. doi: 10.3934/nhm.2010.5.1 |
[9] | Avner Friedman . PDE problems arising in mathematical biology. Networks and Heterogeneous Media, 2012, 7(4): 691-703. doi: 10.3934/nhm.2012.7.691 |
[10] | Bendong Lou . Self-similar solutions in a sector for a quasilinear parabolic equation. Networks and Heterogeneous Media, 2012, 7(4): 857-879. doi: 10.3934/nhm.2012.7.857 |
A horizontal non-homogeneous field adversely affects the seismic resistance of both the utility tunnel and its internal pipes, with seismic waves obliquely incident on the underground structure causing more significant damages. To address these issues, this study, based on a viscous-spring artificial boundary, derives and validates the equivalent junction force formula for the horizontal non-homogeneous field. It then establishes a three-dimensional finite element model of the utility tunnel, pipes, and surrounding soil to obtain the acceleration and strain responses of the utility tunnel and its internal pipes under seismic loading. Finally, it investigates the impact of different incidence angles of shear waves (SV waves) on the response of the utility tunnel and its internal pipes. It was found that as the PGA increases from 0.1 to 0.4 g, both peak acceleration and strain of the utility tunnel and its internal pipes increase. The peak acceleration of the utility tunnel and pipes initially decreases and then increases with the angle of incidence, while the strain increases with the angle of incidence, reaching its peak value when the angle of incidence is 30°. The acceleration and strain responses of the utility tunnel and pipe are higher in sand than in clay, with the peak acceleration strongly correlating with the angle of incidence of ground shaking. The findings of this study provide valuable insights into the seismic design of horizontal non-homogeneous field utility tunnel systems.
Nonlinear partial differential equations (PDEs) impart multi-scale characteristics to the system, thereby allowing for a more accurate prediction of the transmission process of soliton solutions. In practical uses, nonlinear PDEs and soliton solutions are vital for characterizing various phenomena in science and engineering such as biology, physics, ocean engineering, and many more [1,2,3]. Various types of soliton solutions have been reported for integrable systems. For instance, horse-shoe like soliton and lump chain solitons have been studied for the elliptic cylindrical Kadomtsev–Petviashvili equation [4]. Yang et al. analyzed degenerating lump chains into anomalously scattered lumps for the Mel'nikov equation [5]. In literature [6], a series of ripple waves with decay modes for the (3+1)‑dimensional Kadomtsev–Petviashvili equation have been reported. Rogue wave solutions to the (3+1)-dimensional Korteweg-de Vries Benjamin-Bona-Mahony equation were studied via the Hirota bilinear approach [7]. The propagation features and interactions of Rossby waves soliton of the geophysical equation were studied [8]. Breather, lump, and its interaction solutions for the higher dimensional evolution equation were studied[9]. Multisoliton solutions for the variable coefficient Schrödinger equation has been explored in the literature [10]. Some other solitons solutions have been reported for the regularized long-wave equation [11], the Sharma-Tasso-Olver-Burgers equation [12], the modified Schrödinger's equation [13], the complex Ginzburg–Landau equation [14], the (2+1) dimensional Chaffee–Infante equation [15], and many more [16,17,18].
Stochastic differential equations (DEs) deal with phenomena having randomness or uncertainties. Stochastic DEs can be used in various field of science and engineering [19,20,21]. Solving stochastic nolinear PDEs is very challenging and hard due to randomness. Therefore, various methods have been introduced and implemented to derive solutions of stochatics PDEs such as the modified tanh method [22], the modified Kudrayshov technique [23], the Sardar subequation method [24], and many more [25,26].
Fractional operators (FOs) have been frequently used for modelling the physical phenomena in various fields due to its memory process [27,28,29]. In literature, several FOs have been constructed by researchers and scientists [30,31,32]. Most of them do not satisfy some properties such as the chain and quotient rules. A few years ago, Atangana [33] defined a local FO called beta derivative, which generalized the classical operator. The {beta derivative (BD)} is defined as follows:
DβxΨ(x)=dβΨdxβ=limh0→0Ψ(x+h0(x+1Γ(β))1−β−Ψ(x))h0,0<β≤1. |
Here, the BD has the following characteristics: For every real numbers, m and n:
(1)DβxΨ(x)=(x+1Γ(β))1−βdΨdx.(2)Dβx(mΨ+nΦ)=m(x+1Γ(β))1−βdΨdx+n(x+1Γ(β))1−βdΦdx.(3)Dβx(Ψ∘Φ(x))=(x+1Γ(β))1−βdΨdxΦ′(x)(Ψ′(x)).(4)DβxΨ(m)=0. |
The BD has been used for the analysis of soliton solutions with the fractional behavior of nonlinear PDEs [34,35,36]. This work modifies the Benjamin-Bona-Mahony equation (BBME) as follows:
Mt+6MDβxM+DβxxxM−ρDβxxMt=τ(M−ρDβxxM)dPdt, | (1.1) |
where ρ is real parameter, M=M(x,t) is a real valued wave profile, τ is the intensity of sound, and P=P(t) is a white noise having the following properties:
(i)P possesses constant trajectories.
(ii)P(0)=0.
(iii)P(tj+1)−P(tj) has a normal standard distribution.
When we consider τ=0 and β=1, we get the BBME as follows:
Mt+6MMx+Mxxx−ρMxxt=0. | (1.2) |
Benjamin, Bona, and Mahony examined equation (1.2) as an adjustment to the KdV equation. The BBME has been used to analyze the prorogation of long surface gravity pulses with small amplitudes. There are several studies on the BBME. For instance, BBME was studied by using the variational method [37], the deep learning method [38], the generalized exp-function method [39], and many more [40,41]. In [42], the authors have used the F-expansion method to study the solitary waves BBME under BD with white noise. In this paper, we use two advanced analytical methods to deduce more solitary waves solutions and to study the influence of the BD and the white noise.
This section provides the general procedure of the suggested approaches that one can use to find solitary and other waves solutions.
Here, we present the general procedure of the G′G′+G+A-expansion technique. Consider a PDE under space BD as follows
A1(M,∂βxM,∂tM,∂βx∂βxM,∂βx∂tM,∂t∂tM,⋯)=0, | (2.1) |
where A1 is a polynomial in M=M(x,t) and its partial derivatives. To use the proposed procedure, one should abide by the following:
Step 1. First using the wave transformation, one can obtain ODE as follows:
M(x,t)=M(ω1)eτP(t)−12τ2t, | (2.2) |
where ω1=ξ1β(x+1Γ(β))β+ξ2t. Additionally, ξ1 and ξ2 are referred to as the wave speed and the wave number, respectively. By inserting Eq (2.2) in Eq (2.1), the following will be obtained:
A1(M,M′,M′′,M′′′)=0, | (2.3) |
where the ordinary derivatives of different orders are indicated by primes.
Step 2. According to the proposed strategy, we examine the following form for the solution to Eq (2.3):
M(ω1)=ℵ∑i=0Fi(G′(ω1)G′(ω1)+G(ω1)+A)i, | (2.4) |
where Fi is the function of the polynomial's coefficients (G′G′+G+A)i,i=0,1,2,…,ℵ. Assume that G(ω1) is a function that fulfills the subsequent relation:
G′′+AG′+BG+BA=0. | (2.5) |
The value of ℵ can be determined using the homogeneous balance rule (HBR) between the highest nonlinear term and the highest order derivative in Eq (2.3).
Step 3. In this step, the result obtained from the substitution of Eq (2.4) into Eq (2.3) and the coefficients of various powers of (G′G′+G+A) should be compared in terms of A,B,ξ1,ξ2, and i=0,1,2,…,ℵ. Using Mathematica or any other mathematical package, one can determine the solution's values G in the term (G′G′+G+A), and ultimately for the principles of (G′G′+G+A), Fi and ω1. In doing so, the solution of Eq (2.2) can be obtained.
Here, we present the general procedure of applying the modified G′G2-expansion approach to obtain the wave solutions of a nonlinear PDE. This approach contains the following expansion:
M(ω1)=F0+ℵ∑i=1(Fi(G′(ω1)G(ω1)2)i+Si(G′(ω1)G(ω1)2)−i), | (2.6) |
where G(ω1) satisfies the following the equation:
G′′(ω1)=ΨG′(ω1)2G(ω1)2+ψG′(ω1)+2G′(ω1)2G(ω1)+ϖG(ω1)2, | (2.7) |
where Ψ,ψ, and ϖ are the arbitrary constants. Next, one should find the value of ℵ as previously mentioned. Then, substituting Eq (2.6) and using Eq (2.7) into Eq (2.3), one can obtain a differential equation in G(ω1). Then, collecting those terms which contain (G′G2)i, (i=0,1,2,…,n), and setting all the coefficients of (G′G2)i equal to zero, one can acquire a system of algebraic equations. Solving the obtained system can possibly result in the following families.
Family 1. If Ψϖ>0 and ψ=0, the we have the following:
G′G2=√Ψϖ(p1cos(ω1√Ψϖ)+p2sin(ω1√Ψϖ))ϖ(p2cos(ω1√Ψϖ)−p1sin(ω1√Ψϖ)), | (2.8) |
where p1,p2,Ψ, and ϖ are arbitrary constants.
Family 2. If Ψϖ<0 and ψ=0, then we have the following:
G′G2=−√Ψϖ(p1sinh(2ω1√Ψϖ)+p1cosh(2ω1√Ψϖ)+p2)ϖ(p1sinh(2ω1√Ψϖ)+p1cosh(2ω1√Ψϖ)+p2). | (2.9) |
Here, we explore the wave solutions for the proposed stochastic BBME under BD as given in Eq (1.1) with the following procedure:
M(x,t)=M(ω1)eτP(t)−12τ2t. | (3.1) |
Furthermore, we have the following:
Mt=(ξ2M′+τMPt+12τ2M−12τ2M)eτP(t)−12τ2t, | (3.2) |
and
DβxxMt=(ξ21ξ2M′′′+τPtξ21M′′)eτP(t)−12τ2t,DβxM=(ξ1M′)eτP(t)−12τ2t,DβxxxM=(ξ31M′′′)eτP(t)−12τ2t. | (3.3) |
Inserting Eq (3.1) into Eq (1.1) and using Eqs (3.2) and (3.3), we obtain the following:
ξ2M′+(ξ31−ρξ21ξ2)M′′′+6ξ1MM′e−12τ2tEeτP(t)=0. | (3.4) |
By considering P(t), the Gaussian process, and EeτP(t)=e12τ2t, then, Eq (3.4) becomes:
ξ2M′+(ξ31−ρξ21ξ2)M′′′+6ξ1MM′=0. | (3.5) |
Integrating Eq (3.5) one time while considering the integration constant to be zero, we obtain the following:
ζM+M′′+ηM2=0, | (3.6) |
where
ζ=ξ2ξ31−ρξ21ξ2,η=3ξ21−ρξ1ξ2. |
In Eq (3.6), by using the homogeneous balance principle, we obtain ℵ=2. Now, we have Eq (2.4) in the following form:
M1(ω1)=F0+F1(G′G′+G+A)+F2(G′G′+G+A)2. | (3.7) |
Inserting the solution of Eq (3.7) with Eq (2.5) into Eq (3.6), the polynomial of the left side will be in (G′G′+G+A)i,i=0,1,2⋯ℵ. By further equating the coefficients of various powers of (G′G′+G+A) to zero, we obtain a system of algebraic equations. Using Mathematica to solve the system of equations, we obtain the following sets:
{F0=ξ1ξ2(A2−12AB+4B(3B+2))−ξ2√ξ21(A2−4B)26ξ1√ξ21(A2−4B)2,F1=∓2ξ2(A−2B)(A−B−1)√ξ21(A2−4B)2,F2=2ξ2(−A+B+1)2√ξ21(A2−4B)2,ρ=ξ21√ξ21(A2−4B)2+ξ41ξ2ξ31. | (3.8) |
Now, inserting the parameter values presented in Eq (3.8) into Eq (3.7), we get the exact solutions of Eq (3.6) in the following two cases:
Set 1. For D=A2−4B>0, we have the following:
M(ω1)=(ξ1ξ2(A2−12AB+(12B2+8B))−ξ2√ξ21(A2−4B)26ξ1√ξ21(A2−4B)2−(2ξ2(A−2B)(A−B−1))(ν2e√Dω1(A−√D)+ν1(√D+A))√ξ21(A2−4B)2(ν2e√Dω1(−√D+A−2)+ν1(√D+A−2))(2ξ2(−A+B+1)2)(ν2e√Dω1(A−√D)+ν1(√D+A)ν2e√Dω1(−√D+A−2)+ν1(√D+A−2))2√ξ21(A2−4B)2)eτP(t)−12τ2t, | (3.9) |
where ν1 and ν2 remain constants.
Set 2. For D=A2−4B<0, we have the following:
M(ω1)=(ξ1ξ2(A2−12AB+(12B2+8B))−ξ2√ξ21(A2−4B)26ξ1√ξ21(A2−4B)2−(2ξ2(A−2B)(A−B−1))√ξ21(A2−4B)2(Aν2+ν1√−D)sin(√−D2)+(Aν1−ν2√−D)cos(√−D2)((A−2)ν2+ν1√−D)sin(√−D2)+((A−2)ν1−ν2√−D)cos(√−D2)(2ξ2(−A+B+1)2)√ξ21(A2−4B)2((Aν2+ν1√−D)sin(√−D2)+(Aν1−ν2√−D)cos(√−D2)((A−2)ν2+ν1√−D)sin(√−D2)+((A−2)ν1−ν2√−D)cos(√−D2))2)eτP(t)−12τ2t. | (3.10) |
Since the highest-order nonlinear term and the highest-order derivative term are balanced according to the homogenous balance principle in Eq (3.6), we know that the balance number is ℵ=2. Therefore, we have the following:
M(ω1)=F0+F1G′G2+F2(G′G2)2+S1G′G2+S2(G′G2)2. | (4.1) |
Inserting Eq (4.1) with aid of Eq (2.7) into Eq (3.6), and following the same procedure as earlier, we obtain the following:
F1=−2Ψξ21ψ4ρΨξ21ϖ+ρξ21(−ψ2)+1,F2=−2Ψ2ξ214ρΨξ21ϖ+ρξ21(−ψ2)+1,S1=0,S2=0,ξ2=ξ31(4Ψϖ−ψ2)4ρΨξ21ϖ+ρξ21(−ψ2)+1,F0=−2Ψξ21ϖρξ21(4ϖΨ−ψ2)+1. | (4.2) |
Putting the values of the parameters presented in Eq (4.1) into Eq (3.6) and making use of Eqs (2.8) and (2.9), we obtain the following exact solutions.
Family 1. If Ψϖ>0 and ψ=0, then we have the following:
M(ω1)=(−(2Ψ2ξ21)(√Ψϖ(p1cos(ω1√Ψϖ)+p2sin(ω1√Ψϖ))ϖ(p2cos(ω1√Ψϖ)−p1sin(ω1√Ψϖ))2)24Ψξ21ϖρ+1−2Ψξ21ϖ4Ψξ21ϖρ+1)eτP(t)−12τ2t. | (4.3) |
Family 2. If Ψϖ<0 and ψ=0, then we have the following:
M(ω1)=(−(2Ψ2ξ21)(−√Ψϖ(p1sinh(2ω1√Ψϖ)+p1cosh(2ω1√Ψϖ)+p2)(ϖ(p1sinh(2ω1√Ψϖ)+p1cosh(2ω1√Ψϖ)+p2))2)24Ψξ21ϖρ+1−2Ψξ21ϖ4Ψξ21ϖρ+1)eτP(t)−12τ2t. |
This portion of the present work graphically visualize the obtained solutions and presents some physical interpretations and discussions on the obtained results. In Figure 1, solution (3.9) with particular values (i.e, ν1=5,ν2=−.5,ξ1=−.2,ξ2=−1,A=3,B=2.6,τ=0,P=0) is visualized. In Figure 1, the value of β is varied while the noise intensity τ is considered as zero. The β is used as 1, 0.9, and 0.8 for subfigures (1a, 1d), (1b, 1e), and (1c, 1f), respectively. Here, we observed the dark soliton wave, where we see that as the fractional order decreases when the wave separation is increased.
Furthermore, Figure 2 shows the dynamics of the exact solution (2.2) by varying the noise intensity while keeping the β=0.95. Other parameters are used for the simulation of Figure 1. The τ is used as 0.1, 0.4, and 0.9 for subfigures (2a, 2d), which is (2b, 2e), and (2c, 2f), respectively. In Figure 2, one can observe the affects of noise on the dynamics of the solution, which is simulated here. Furthermore, the dynamics of the exact solution (3.10) are visualized in Figures 3 and 4 by varying β and τ, respectively. In the simulation of these figures, the parameters are selected in the form ν1=.5,ν2=1,ξ1=−.7,ξ2=.5,p1=2,p2=1,A=−4,B=0,τ=0,P=0; alternatively in Figure 3, the τ is considered as zero. and in Figure 4. the β is fixed as 0.95. The β is used as 1, 0.9, and 0.8 for subfigures (3a, 3d), (3b, 3e), and (3c, 3f), respectively. Similarly, τ is used as 0.2, 0.5, and 0.8 for subfigures (4a, 4d), (4b, 4e), and (4c, 4f), respectively. Here, we observed the interaction of the bright wave with a kink wave, where the amplitude of the bright wave decreases as the β decreases in the negative region of the spatial coordinate.
In Figure 5, the solution (3.9) with particular values (i.e, ν1=5,ν2=−.5,ξ1=−.2,ξ2=−1,A=3,B=2.6,τ=0, and P=0) is visualized. In Figure 5, the various values for β are considered, while the noise intensity τ is supposed to be zero. The β is considered as 1, 0.95, and 0.9 for subfigures (5a, 5d), (5b, 5e), and (5c, 5f), respectively. Here, we observed the hybrid bright-dark soliton wave, where we see that as the fractional order decreases when then amplitude of the dark solitons increases and the bright soliton is decreases.
Moreover, Figure 6 shows the dynamics of the exact solution (3.9) by varying the noise intensity while keeping the β=0.95. Other parameters are used for the simulation of Figure 5. The τ is used as 0.5, 0.6, and 0.9 for subfigures (6a, 6d), (6b, 6e), and (6c, 6f), respectively. In Figure 6, one can observe the affects of noise on the dynamics of the solution, which is simulated here; it can be seen that the highest and lowest amplitude areas become more random as τ increases.
Furthermore, the dynamics of the exact solution (3.10) are visualized in Figures 7 and 8 by varying β and τ, respectively. In the simulation of these figures, the parameters are selected in the form ρ=1,ϖ=−.1,ξ1=1,Ψ=1,p1=1,p2=1,P=0, and τ=0; alternatively, in Figure 7, the τ is considered as zero, and in Figure 8, the β is fixed as 0.95. The β is used as 1, 0.9, and 0.8 for subfigures (7a, 7d), (7b, 7e), and (7c, 7f), respectively. Similarly, τ is used as 0.05, 0.3, and 0.6 for subfigures (8a, 8d), (8b, 8e), and (8c, 4f), respectively. Here, we observed the periodic wave solution, where the amplitude of the periodic waves decreases as the β decreases in the negative region of the spatial coordinate. Furthermore, we see that the wave profile behaves more randomly in areas where the amplitude is either low or high. Thus, from these analyses, it can be noticed that the obtained results are more generalized than the solutions reported in previous papers. Indeed, when the BD operators equals one, the solution converges to the stochastic integer order solutions. If the intensity of the white noise is zero, then the solutions converge to a deterministic case. When β=1 and τ=0, the obtained solutions converge to the determinsitic case.
This study has explored the stochastic BBME with the BD, thereby incorporating multiplicative noise in the Itô sense. We have derived various analytical soliton solutions for these equations by utilizing two distinct expansion methods, both within the framework of beta derivatives. A fractional multistep transformation was employed to convert the equations into nonlinear forms with respect to an independent variable. After performing algebraic manipulations, the solutions were found to be trigonometric and hyperbolic trigonometric functions. Our analysis demonstrated that the wave behavior was influenced by the fractional-order derivative in the proposed equations, thus providing deeper insights into the wave composition as the fractional order increases or decreases. Additionally, we examined the effect of white noise on the propagation of wave solutions. This study has underscored the computational robustness and adaptability of the proposed approach to investigate various phenomena in the physical sciences and engineering.
Conceptualization: M.S.D.S. Methodology: K.A.A. Software: S.S. Validation: A.K. Formal analysis: A.K. Investigation: M.H. Writing-original draft preparation: K.A.A. Writing-review and editing: H.S., A.M.
The Researchers would like to thank the Deanship of Graduate Studies and Scientific Research at Qassim University for financial support (QU-APC-2024-9/1). The authors wish to extend their sincere gratitude to the Deanship of Scientific Research at the Islamic University of Madinah.
All authors declare no conflicts of interest in this paper.
[1] |
Canto-Perello J, Curiel-Esparza J (2013) Assessing governance issues of urban utility tunnels. Tunnelling Underground Space Technol 33: 82–87. https://doi.org/10.1016/j.tust.2012.08.007 doi: 10.1016/j.tust.2012.08.007
![]() |
[2] | American Public Works Association (1971) Feasibility of utility tunnels in urban areas. APWA Special Report No. 39. |
[3] |
Cano-Hurtado JJ, Canto-Perello J (1999) Sustainable development of urban underground space for utilities. Tunnelling Underground Space Technol 14: 335–340. https://doi.org/10.1016/S0886-7798(99)00048-6 doi: 10.1016/S0886-7798(99)00048-6
![]() |
[4] | Tang AP, Feng RC, Zhou XY, et al. (2008) Seismic response characteristics of shallow-buried utility tunnel systems. The 14th World Conference on Earthquake Engineering. |
[5] |
Jiang L, Chen J, Li J (2010) Seismic response of underground utility tunnels: shaking table testing and FEM analysis. Earthq Eng Eng Vib 9: 555–567. https://doi.org/10.1007/s11803-010-0037-x doi: 10.1007/s11803-010-0037-x
![]() |
[6] |
Chen J, Jiang L, Li J, et al. (2012) Numerical simulation of shaking table test on utility tunnel under non-uniform earthquake excitation. Tunnelling Underground Space Technol 30: 205–216. https://doi.org/10.1016/j.tust.2012.02.023 doi: 10.1016/j.tust.2012.02.023
![]() |
[7] | Sitharam TG, Anitha Kumari SD (2015) Numerical Simulations of Tunnels Using DEM and FEM. Paper presented at the 13th ISRM International Congress of Rock Mechanics, Montreal, Canada. |
[8] |
Konstandakopoulou FD, Beskou ND, Hatzigeorgiou GD (2021) Three-dimensional nonlinear response of utility tunnels under single and multiple earthquakes. Soil Dyn Earthquake Eng 143: 106607. https://doi.org/10.1016/j.soildyn.2021.106607 doi: 10.1016/j.soildyn.2021.106607
![]() |
[9] |
Yu H, Zhang Z, Chen J, et al. (2018) Analytical solution for longitudinal seismic response of tunnel liners with sharp stiffness transition. Tunnelling Underground Space Technol 77: 103–114. https://doi.org/10.1016/j.tust.2018.04.001 doi: 10.1016/j.tust.2018.04.001
![]() |
[10] |
Zhang J, He C, Geng P, et al. (2019) Improved longitudinal seismic deformation method of shield tunnels based on the iteration of the nonlinear stiffness of ring joints. Sustainable Cities Soc 45: 105–116. https://doi.org/10.1016/j.scs.2018.11.019 doi: 10.1016/j.scs.2018.11.019
![]() |
[11] |
Zhang W, Restrepo D, Crempien JG, et al. (2021) A computational workflow for rupture-to-structural-response simulation and its application to Istanbul. Earthquake Eng Struct Dyn 50: 177–196. https://doi.org/10.1002/eqe.3377 doi: 10.1002/eqe.3377
![]() |
[12] |
Zhang W, Taciroglu E (2021) 3D time‐domain nonlinear analysis of soil‐structure systems subjected to obliquely incident SV waves in layered soil media. Earthquake Eng Struct Dyn 50: 2156–2173. https://doi.org/10.1002/eqe.3443 doi: 10.1002/eqe.3443
![]() |
[13] |
Panji M (2023) A half-space TD-BEM model for a seismic corrugated orthotropic stratum. Eng Anal Boundary Elem 152: 655–677. https://doi.org/10.1016/j.enganabound.2023.04.032 doi: 10.1016/j.enganabound.2023.04.032
![]() |
[14] |
Kavandi P, Ganjian N, Panji M (2024) A DR-BEM approach for analyzing the transient SH-wave scattering problems: A comparative study. Eng Anal Boundary Elem 169: 105962. https://doi.org/10.1016/j.enganabound.2024.105962 doi: 10.1016/j.enganabound.2024.105962
![]() |
[15] |
Mojtabazadeh-Hasanlouei S, Panji M, Kamalian M (2022) Scattering attenuation of transient SH-wave by an orthotropic gaussian-shaped sedimentary basin. Eng Anal Boundary Elem 140: 186–219. https://doi.org/10.1016/j.enganabound.2022.04.023 doi: 10.1016/j.enganabound.2022.04.023
![]() |
[16] |
Mojtabazadeh-Hasanlouei S, Panji M, Kamalian M (2022) Attenuated orthotropic time-domain half-space BEM for SH-wave scattering problems. Geophys J Int 229: 1881–1913. https://doi.org/10.1093/gji/ggac032 doi: 10.1093/gji/ggac032
![]() |
[17] |
Li M, Zhang K, Meng K, et al. (2023) Seismic response study of shield tunnel with lateral karst cavity under SV waves. Structures 56: 105023. https://doi.org/10.1016/j.istruc.2023.105023 doi: 10.1016/j.istruc.2023.105023
![]() |
[18] |
Li M, Meng K, Zhou J, et al. (2024) Seismic response and damage analysis of shield tunnel with bottom Karst caves under oblique SV waves. Nat Hazards 120: 2731–2747. https://doi.org/10.1007/s11069-023-06302-5 doi: 10.1007/s11069-023-06302-5
![]() |
[19] |
Zarzalejos JM, Aznárez JJ, Padrón LA, et al. (2014) Influences of type of wave and angle of incidence on seismic bending moments in pile foundations. Earthquake Eng Struct Dyn 43: 41–59. https://doi.org/10.1002/eqe.2330 doi: 10.1002/eqe.2330
![]() |
[20] |
Huang JQ, Du XL, Zhao M, et al. (2017) Impact of incident angles of earthquake shear (S) waves on 3-D non-linear seismic responses of long lined tunnels. Eng Geol 222: 168–185. https://doi.org/10.1016/j.enggeo.2017.03.017 doi: 10.1016/j.enggeo.2017.03.017
![]() |
[21] |
Ye Z, Li Y, Li L, et al. (2025) Seismic response of a long shield tunnel crossing through multiple soil deposits. Comput Geotech 177: 106892. https://doi.org/10.1016/j.compgeo.2024.106892 doi: 10.1016/j.compgeo.2024.106892
![]() |
[22] |
Yan L, Haider A, Li P, et al. (2020) A numerical study on the transverse seismic response of lined circular tunnels under obliquely incident asynchronous P and SV waves. Tunnelling Underground Space Technol 97: 103235. https://doi.org/10.1016/j.tust.2019.103235 doi: 10.1016/j.tust.2019.103235
![]() |
[23] | Du XL, Li LY (2008) A Viscoelastic Artificial Boundary for Near-field Fluctuation Analysis of Saturated Porous Media. Chin J Geophys 51: 575–581. |
[24] | Qu J (2017) Dynamic response of comprehensive pipe gallery under oblique incident seismic wave. Harbin Institute of Technology. |
[25] | Zhao J, Lan WJ, Teng WQ, et al. (2021) Three-dimensional seismic response analysis of submarine immersed tube tunnel considering oblique incidence of P wave. Struct Eng 37: 73–81. |
[26] | Wang Z, Jing LP, Lu XY, et al. (2023) Comparative Study on Viscoelastic Artificial Boundary Elements and round Motion Input Methods. World Earthquake Eng 39(2). |
[27] | Wang F, Song Z, Liu Y, et al. (2023) Response Characteristics and Tensile Failure Evaluation of Asphalt Concrete Core Wall under Oblique Incidence of SV Wave Space. Chin J Geotech Eng 45: 1733–1742. |
[28] | Du XL, Xu HB, Zhao M (2015) Analysis of seismic response of high arch dam under SV wave oblique incidence. J Hydropower Gener 34: 139–145 |
[29] |
Huang ZY, Feng YZ, Tang AP, et al. (2023) Influence of oblique incidence of P-waves on seismic response of prefabricated utility tunnels considering joints. Soil Dyn Earthquake Eng 167: 107797. https://doi.org/10.1016/j.soildyn.2023.107797 doi: 10.1016/j.soildyn.2023.107797
![]() |
[30] |
Yu J, Wang ZZ (2021) The dynamic interaction of the soil-tunnel-building system under seismic waves. Soil Dyn Earthquake Eng 144: 106686. https://doi.org/10.1016/j.soildyn.2021.106686 doi: 10.1016/j.soildyn.2021.106686
![]() |
[31] | Zhou Tl, Dong CS, Li S, et al. (2024) Seismic response of tunnel with damping layer under oblique incident SV wave. World Earthquake Eng 40: 1–12. |
[32] |
Lysmer J, Kuhlemeyer RL (1969) Finite dynamic model for infinite media. J Eng Mech Div 95: 859–877. https://doi.org/10.1061/JMCEA3.0001144 doi: 10.1061/JMCEA3.0001144
![]() |
[33] | Xu CS, Hu ZY, Shi YB, et al. (2022) Seismic response law of prefabricated integrated utility tunnel interface under oblique incidence of SV wave. J Beijing Univ Technol 48: 1215–1225. |
[34] |
Wang GB, Yuan MZ, Miao Y, et al. (2018) Experimental study on seismic response of underground tunnel-soil-surface structure interaction system. Tunnelling Underground Space Technol 76: 145–159. https://doi.org/10.1016/j.tust.2018.03.015 doi: 10.1016/j.tust.2018.03.015
![]() |
1. | Peng E, Tingting Xu, Linhua Deng, Yulin Shan, Miao Wan, Weihong Zhou, Solutions of a class of higher order variable coefficient homogeneous differential equations, 2025, 20, 1556-1801, 213, 10.3934/nhm.2025011 |