Research Article  Open Access
Influence of the Hydraulic Boundary Condition between the Embankment and Saturated HalfSpace on the TrainInduced Ground Vibration
Abstract
A threedimensional analytical model was proposed to investigate the influence of the hydraulic boundary of the ground surface on the traininduced ground vibrations. The ground was simulated as a fully saturated poroelastic halfspace and the embankment as a rectangular elastic soil layer with finite width. The rails and the sleepers were modeled as Euler beams and a Kirchhoff plate, respectively. Two hydraulic boundaries of the interface between the embankment and the ground surface were simplified as two limiting cases, namely, the permeable and impermeable cases. The linearized dynamic equations of motion for the fully saturated poroelastic soil and the elastic embankment were solved by Fourier transform and Fourier series. The vertical velocity and pore pressure were firstly calculated in the frequencywavenumber domain and then transformed into the timespace domain by inverse Fourier transform. The influence of the hydraulic boundary at the ground surface on the traininduced ground vibration was specially investigated. It was found that the hydraulic boundary condition has a significant influence on groundborne vibrations for the highspeed moving train with high load excitation frequencies.
1. Introduction
The highspeed trains usually cause disturbance to residents and detrimental effects on the buildings nearby. In the past decades, researchers have made a lot of efforts to investigate the dynamic response of the ground generated by moving traffic loads.
Some researchers [1–3] have established several threedimensional models to study the dynamic response of a homogeneous elastic halfspace due to the moving constant load applied on the ground surface. In order to include the track, Metrikine and Dieterman [4] and Kim [5, 6] investigated the dynamic interaction between an Euler beam and the underlying elastic halfspace, and it was found that resonance phenomenon occurs when the moving load velocity approaches the critical velocity of the beamsoil model. Sun [7] presented a theoretical threedimensional model to simulate the vibrations of a rigid pavement resting on an elastic halfspace subjected to a moving point load for the subsonic, transonic, and supersonic speed, respectively. Using a simplified formulation of the railway track, namely, a layered beam structure resting on a layered elastic ground, Sheng et al. [8–10] proposed a series of analytical models to study the ground vibrations by solving the wave equations for the Euler beam and the viscoelastic halfspace.
All the works mentioned above treated the soil as a singlephase elastic and viscoelastic medium. However, underground water usually exists, which leads to the fact that the soil is saturated with water and becomes a twophase material. For the twophase medium, the coupling effects between the soil skeleton and the porefluid cannot be neglected when subjected to moving train loads [11]. If the halfspace is modeled as a saturated poroelastic medium by using Biot theory [12], the porefluid related parameters and hydraulic boundary conditions can be considered. By neglecting the coupling between the soil skeleton and the fluid, Siddharthan et al. [13] solved Biot equations of up formulation under the plane strain condition, and the displacement and pore pressure responses caused by the moving load with a low velocity were given. Based on Siddharthan work, Theodorakopoulo et al. [14] and Theodorakopoulo [15] solved the complete uwp equation and presented a twodimensional model to study the dynamic responses of the saturated soil under rectangular moving loads. Jin [16] established a threedimensional model to calculate the dynamic response of the saturated soil for the case of a point load moving on the ground surface. By adopting four scalar potential functions and Helmholtz decomposition theorem, Liu et al. [17] studied the dynamic response of saturated soil to a harmonic moving load under conditions of plane strain. Xu et al. [18] researched the vibrations of the infinite Euler beam on a multilayered saturated halfspace induced by moving loads on the beam, in which the effects of soil properties and the load velocity were studied. Gao et al. [19, 20] presented a trackground interaction model to evaluate the traininduced vibration in the saturated ground. The influence of permeability coefficient and shear wave velocity on the attenuation of ground vibration was analyzed in detail. Chahour et al. [21] carried out a spectral analysis of a railway track coupled to the multilayered poroelastic halfspace for the case of a moving harmonic load moving on the rails. However, in the above works the stress boundary condition between the beam and the underlying halfspace was simplified as uniform stress distribution in the transverse direction, and the displacement compatibility condition is only valid at the center line of the track, which could not model the real dynamic interaction between the railway track and the ground precisely when subjected to high frequency loads. Based on Biot’s theory, Cai et al. [22] and Cao et al. [23] established a trackembankmentsoil model with continuous displacement and stress compatibility conditions at the trackground interface to investigate the ground vibration induced by a moving train.
It should be noted that a fully permeable boundary condition is usually assumed at the ground surface or between the track structure and underlying halfspace in the above studies. However the hydraulic boundary condition between the trackembankment and ground may deteriorate during the service time due to the breakage of road bed filling particles; thus the contact surface between the track and halfspace is actually partially permeable. The partially permeable interface between the embankment and ground and the fully permeable free ground surface leads to a mixed boundary value problem, which is not easy to be solved analytically. To circumvent the difficulty of the mixed boundary value problem, two limiting cases, namely, fully permeable and fully impermeable ground surface, were used to bound the actual hydraulic boundary condition. A study by Zhou et al. [24] showed that the dynamic response of a saturated halfspace for the permeable and impermeable ground surface exhibits an obvious difference; however the dynamic interaction between the track system and the ground was not considered in this study. To date, few studies focus on the effects of the hydraulic boundary condition on the traininduced ground vibration. In order to make accurate prediction of traininduced vibrations, the effect of the hydraulic boundary condition on the traininduced ground vibration needs to be investigated.
In the present paper, the effects of the hydraulic boundary condition between the embankment and the saturated halfspace on the traininduced ground vibrations are investigated by a threedimensional semianalytical model. The model consists of rails, sleepers, an embankment, and a halfspace which are simulated as Euler beams, a Kirchhoff plate, an elastic layer, and a poroelastic soil medium, respectively. By using Fourier transform and Fourier series techniques, the linearized dynamic equations of motion for the embankment and fully saturated poroelastic halfspace were solved in the frequencywavenumber domain. The critical velocity of the embankmentground structure under different hydraulic boundary conditions was investigated. The effects of the hydraulic boundary condition on the dynamic responses of the ground generated by the moving trains with different selfexcitation frequencies were analyzed.
2. Governing Equations and Solutions
2.1. Governing Equations
A detailed trackembankmentground model is shown in Figure 1, consisting of a halfspace, an embankment, sleepers, and two rails. The embankment is modeled as a rectangular elastic layer with finite width 2 and finite thickness which is resting on a homogeneous poroelastic halfspace. The sleeper and rails are simulated as a Kirchhoff plate and EulerBernoulli beams. A highspeed train with 5 carriages is used to simulate the moving train loads and the detailed load amplitude and geometric distribution are presented in Figure 2. The drainage condition of interface between the embankment and the ground usually deteriorates during the service period due to the particle breakage and the invasion of fines. Thus, the interface is actually partially impermeable; then the hydraulic boundary of the ground surface becomes a mixed boundary with the partially permeable interface and the permeable free surface. In order to facilitate the analytical solutions of this model, two limiting cases, namely, the permeable and impermeable ground surface, are considered in the present study to bound the effects of the hydraulic boundary at the ground surface.
The Lamé equations of the embankment layer are given as [27]where the subscript e denotes the embankment and the dot over the variable represents the derivative with respect to time t. is the displacement vector and , are Lamé constants of the elastic layer.
The traction on a plane with the normal direction (k=x, y, z) is expressed as where is the gradient operator.
For the saturated poroelastic halfspace underneath the embankment, Biot’s theory is introduced:where the solid displacement vector is defined as and the fluid displacement vector is defined as ; λ and μ are Lamé constants of the solid skeleton; , , and are the mass densities of the saturated soil, fluid, and solid skeleton, respectively; , in which is the porosity; is a densitylike parameter that depends on and the geometry of the pores; α and are compressibility coefficients of the fluid and solid skeleton; , in which is the coefficient of fluid viscosity and is the coefficient of permeability.
The constitutive relations of the soil medium can be expressed aswhere is the total stress and is the pore water pressure; ; .
2.2. Boundary Conditions
The boundary conditions applicable to the problem depend on the contact conditions between the embankment and poroelastic halfspace, also the conditions between the embankment and the sleeper. Setting the boundary conditions on the sides of embankment layers () in the following special way is able to facilitate the series expansion of the displacement field in the embankment:where is half of the embankment width, and is the embankment thickness. In the above conditions, and are the natural ones, but is unphysical for enabling computation. The study by Karlström and Boström [25] showed that the unphysical condition can simulate the vertical displacement of the embankment well when only vertical loads are applied. The abovementioned conditions can be simplified asIn the present model, two rails are modeled as infinite EulerBernoulli beams in the vertical and transverse directions and as the rod equation in the longitudinal direction. is the sectional area, is the modulus of elasticity, is the density, c is width, and and are the area moment of inertia about the y and zaxis, respectively. The sleepers are simulated by a transversely isotropic Kirchhoff plate. is the density, is the modulus of elasticity in the direction, and in the direction. To simulate the longitudinal properties of the sleepers, is selected as zero. and are Poisson’s ratio in the yz plane and in the xk plane, respectively, in which denotes the directions normal to the direction. denotes the width of the sleepers, denotes the distance between the rails, h is the height of the sleepers, and is the axle load of the train. On the surface of the embankment (z=d), the stress conditions can be expressed aswhere , , , , , , , .
The longitudinal shear stress can be given by the rod equation of the rails:
The ground is assumed to be bonded with the embankment perfectly; hence the displacement and stress compatibility conditions can be expressed aswhere the displacement vector of the saturated halfspace is defined as and the stress vector is defined as .
For the case of permeable boundary at the ground surface, all stresses from the embankment are transmitted to the solid skeleton of poroelastic ground. Thus the pore pressure is zero at the interface. In contrast, for the case of impermeable boundary condition, stresses from the embankment are transmitted to both the fluid and soil skeleton of poroelastic ground, which implies that the pore water pressure is nonzero. The hydraulic boundary condition can be expressed as follows.
For the permeable boundary condition (z=0):
For the impermeable boundary condition (z=0):
2.3. Analytical Solution
In order to solve Lamé equations of the embankment layer, three potentials for longitudinal (P), horizontal transverse (SH), and vertical transverse (SV) waves, which are denoted by , , and , respectively, are employed to decompose the displacement field (for the detailed solution procedure, see Cao et al. [28]):
The Fourier transform pair with respect to is defined as
The Fourier transform pair with respect to is defined as
Meanwhile to meet the special boundary conditions of the sides of embankment in the direction, Fourier series expansions with respect to the three potentials mentioned above are introduced:where is wavenumber in the direction; and are unknown constants (n=1, 2, 3); and , in which denotes the compression wave speed of the embankment and is calculated as , and is the shear wave speed defined as ; , with and .
Adopting (21)(23) together with (16), the embankment displacements are expressed aswhereThe definition for the parameters , , , and is listed in Table 1.

Substituting (24)(26) into (2), the stresses in the embankment can be derived as follows:where, , and are given as follows:
Analytical solutions to the dynamic equations of motion for a fully saturated poroelastic medium (see (3)(4)) can be obtained by applying a triple Fourier transform.
The Fourier transform pair with respect to is defined as
The dynamic equations of saturated soil can be solved in the transformed domain and the solutions are expressed as
Substituting (42)(44) into (5), the stresses in the poroelastic halfspace are given aswhere A, B, C, and are the unknown constants. Other coefficients are given in the Appendix.
At the surface of the halfspace (z=0), the following equations are obtained by applying Fourier transform with respect to to the stress compatibility condition:where and are defined as
At the interface between the embankment and the poroelastic halfspace (z=0), the drainage conditions are assumed to be fully permeable or fully impermeable. The additional hydraulic boundary condition should be added.
For the permeable boundary condition (z=0):
For the impermeable boundary condition (z=0):
By applying (48)(50) and (53)(54), the unknown coefficients A, B, C, and can be expressed by and .
The application of the displacement continuity condition between the embankment and the poroelastic halfspace over the width of the embankment gives the remaining unknowns and . According to the orthogonality of trigonometric series, each mode of the series gives an independent equation, denoted as m’:where is the Neumann factor, defined as for m’ =0 and for m’ .
For the stress conditions at the top surface of the embankment (see (8)(10)), a similar inverse Fourier series with respect to over the width 2 is employed to obtain the following equations:where, , and can be calculated by the following integral formulas:where is the train axle load in the transformed domain. For the train load with a constant speed, it can be expressed aswhere is the train speed, is the load amplitude, and is the excitation frequency of the train load.
3. Numerical Results and Analysis
Referring to the work by Cao et al. [28], the parameters for the railway track and the embankment are listed in Tables 2 and 3. The parameters describing the properties of the saturated poroelastic halfspace are listed in Table 4. The parameters will be selected from Tables 2–4 unless otherwise specified. The effects of the hydraulic boundary condition between the embankment and saturated poroelastic halfspace on the dynamic response of the embankmentground system are investigated for two moving train speeds (V=70 km/h and V=250 km/h) below and above the critical speed, respectively.



3.1. Comparison with Existing Works
In order to verify the present work, results calculated by the present model are compared with those obtained by Karlström and Boström [25], in which the singlephase elastic material was used to model the halfspace instead of using a twophase saturated soil model. A viscoelastic halfspace can be modeled by choosing negligibly small values of the poroelastic parameters (, b, M, and α are set to be 10^{−4}). The vertical displacements at the central line of the top surface of embankment due to a moving train with speed of 70 km/h are presented in Figure 3(a). It can be seen that the two results are in good consistency. To further verify the correctness of the proposed analytical solution, the numerical results obtained by the present model are compared with the field measurements carried out by Madshus and Kaynia [26]. In the field test, the groundborne vibrations generated by X2000 trains running on a soft soil site at Lesgard in Sweden were monitored. In order to make the calculations obtained by the proposed method comparable with the field measurements, the soil and track parameters are selected as those used in the work by Madshus and Kaynia [26]. The measurement position is chosen at the central line of the top surface of embankment, and the train velocity is 185 km/h. As seen in Figure 3(b), in a general trend the results predicted by the proposed analytical model agree well with those measured in the field test.
(a)
(b)
3.2. Dynamic Response of the Ground Generated by Moving Constant Loads
In Figure 4, the maximum vertical velocity at the ground surface (y=0m, z=0m) versus the train speed is presented. Two hydraulic boundary cases, i.e., the permeable and impermeable ground surface, are considered. It is shown that the maximum vertical velocity at the ground surface varies with the train speed and reaches a sharp peak at a certain velocity, which is defined as the critical velocity of the embankmentground system. The critical velocity is near 1.05 V_{s} (V_{s} is the shear wave velocity of the ground); the reason why the critical velocity is higher than the shear wave velocity of the ground is that the embankment is more rigid than the saturated ground. No obvious difference in the vertical velocity for the permeable and impermeable boundary cases is observed at the low train speed range. However, when the train speed approaches the critical speed, the vertical velocity for the two cases becomes quite different.
In Figure 5, the time history of the traininduced vertical velocity at the ground surface (x=0, y=0, z=0) is presented at a train speed of 250 km/h. From Figure 5(a), it can be clearly observed that the vertical velocity reaches a positive peak value at the position of wheel pairs but a negative peak after the wheel pairs have passed by. The largest vertical velocity is 135 mm/s for the permeable hydraulic boundary and twice more than that for the impermeable hydraulic boundary. The frequency spectrum in Figure 5(b) presents a wavelike pattern with peaks which correspond to the geometric distribution of train axle loads under certain train speed. A significant distinction of the vertical velocity can be observed between the two cases, especially at the peaks. The results of Figures 4 and 5 can be explained as follows: for the high train speed, the frequency content spreads in the high frequency range, which implies that obvious relative motion between the solid skeleton and fluid occurs, and as a result the hydraulic boundary of the ground surface has a significant effect on the velocity response.
(a)
(b)
The maximum pore pressure versus depth under a constant moving load for V=70 km/h and V=250 km/h is shown in Figure 6. The effect of the hydraulic boundary condition on the pore water pressure is analyzed. The value of b, a parameter reflecting the permeability of the soil, is chosen to be 1.0 × 10^{6} N·s/m^{4}. It should be pointed out that is inversely proportional to the soil permeability. As shown in Figure 6, for the permeable case the maximum pore pressure increases with depth and reaches a peak value at z=2.5 m and then rapidly decreases to a constant value. For the impermeable case the maximum pore pressure decreases monotonously, and with the depths larger than z=20 m, the pore pressure response becomes very small. The comparison of Figures 6(a) and 6(b) shows that the maximum pore pressure increases with the train speed over the depth considered. For example, when the train speed increases from 70 km/h to 250 km/h, the peak pore pressure increases from 0.015 kPa to 0.22 kPa for the permeable case and from 0.04 kPa to 0.62kPa for the impermeable case. The ground with impermeable surface has a much larger pore water pressure than that with permeable surface. Furthermore, it can be seen that the pore pressure difference under two drainage conditions is mainly distributed in the region near the ground surface and decreases for deeper depths.
(a)
(b)
Figure 7 shows the maximum pore pressure against depth under a moving train for a lowpermeability soil (b=1.0 × 10^{9} N·s/m^{4}). By comparing Figure 6 with Figure 7, it can be found that the traininduced pore pressure in the lowpermeability soil is much larger than the pore pressure in the highpermeability soil.
(a)
(b)
3.3. Dynamic Response of the Ground Generated by Moving Harmonic Loads
In reality, the contact surface between the wheels and rails is usually uneven, and hence the harmonic loads with selffrequency are suitable for modeling the dynamic loads between the track and the rail [10]. In order to investigate the ground vibrations generated by the dynamic wheelrail interaction, the moving harmonic loads are used to simulate the train loads in this section.
The vertical velocity on ground surface generated by a moving train with selfexcitation frequency of 10 Hz is presented in Figures 8 and 9 for the train speeds 70 km/h and 250 km/h, respectively. For both train speeds, the vertical velocity for the permeable case is similar to the values for the impermeable case. Thus, the effect of the hydraulic boundary condition of the ground surface on the velocity response is quite limited for moving trains with a low selfexcitation frequency of 10 Hz.
(a)
(b)
(a)
(b)
The maximum pore pressure in the poroelastic ground with low and high soil permeability coefficients generated by a moving train with selfexcitation frequency =10 Hz is presented in Figures 10 and 11. In Figure 10(a), for the case of b=1.0 × 10^{6} N·s/m^{4} it can be seen that the maximum pore pressure for the impermeable case is larger than that for the permeable case near the ground surface. The values of pore pressure for two cases equal to each other at z=2.5 m, and for deeper depths the maximum pore pressure for the permeable case exceeds the pore pressure for the impermeable case. A similar trend can be observed in Figure 10(b). The traininduced pore pressure in a lowpermeability soil (b=1.0 × 10^{9} N·s/m^{4}) is presented in Figure 11. It is observed that the pore pressure response mainly occurs in the depths shallower than z=10 m, and the effect of the hydraulic boundary condition on the pore pressure is obvious only near the ground surface.
(a)
(b)
(a)
(b)
The vertical velocity on ground surface induced by a moving train with selfexcitation frequency of 30 Hz for two train speeds is presented in Figures 12 and 13. At the train speed V=70 km/h, as shown in Figure 12, the vertical velocity for the impermeable case is about 25% larger than that for the permeable case. The vertical velocity is mainly distributed in the frequency range of 2040 Hz, and most of the frequency components for the impermeable case are larger than those for the permeable case, especially at the peak values. In Figure 13(a), for the train speed of 250 km/h, a huge discrepancy of the vertical velocity between the two cases can be observed, and the waves propagate over a wider area and attenuate more slowly. Figures 12 and 13 show that the frequency components are mainly distributed around 30 Hz, and this is due to the wellknown Doppler effects. The frequency range can be calculated by =/2π(1±V/) where is the upper and lower limit of the frequency range.
(a)
(b)
(a)
(b)
The maximum pore pressure plotted against the depth to a moving train with selfexcitation frequency of 30 Hz is shown in Figure 14 for a high soil permeability (b =1.0 × 10^{6} N·s/m^{4}). In a general trend the pore pressure for the two cases decreases with depth and the fluctuation pattern is due to the wave propagations in the ground. When V=70 km/h, the local maximum value of the pore pressure for the permeable case is 100% larger than that for the impermeable case. When V=250 km/h, the difference of the local maximum value between the two hydraulic boundary conditions is 50%. The difference of the pore water pressure between the two cases disappears gradually as the depth increases. Figure 15 shows the pore pressure in the saturated ground with b=1.0 × 10^{9} N·s/m^{4}. It is observed that the magnitude of the pore pressure for the two cases is similar. This is due to the fact that the interface between the embankment and the halfspace presents an impermeable behavior in lowpermeability soil even assuming the hydraulic boundary condition to be permeable. As shown in Figure 15(b), the phase difference of the pore pressure is obvious at the high train speed of 250 km/h.
(a)
(b)
(a)
(b)
3.4. Influence of the Load Excitation Frequency on the Velocity Response
In order to investigate the effect of load excitation frequency on the velocity response specifically, the velocity response at the ground surface is plotted against the load selfexcitation frequency in Figures 16 and 17. In a general trend the vertical velocity to a moving train firstly increases with the increase of selfexcitation frequency and reaches its maximum value at a load frequency around 10 Hz and then decreases as the frequency increases further. The effect of the hydraulic boundary at the ground surface on the vertical velocity response becomes more significant at higher frequencies. The comparison between Figures 16 and 17 shows that the effect of the hydraulic boundary on the velocity response is reduced with the increase of b (a decrease of the soil permeability) due to the fact that relative motion between the soil skeleton and the fluid is restrained in the lowpermeability ground.
(a)
(b)
(a)
(b)
4. Conclusions
In this work, a threedimensional model is proposed to investigate the influence of hydraulic boundary condition between the embankment and poroelastic halfspace on the freefield response induced by moving trains. The traininduced vibration prediction model consists of a detailed railway track and a saturated ground. The saturated ground is governed by Biot’s theory, and the linearized dynamic equations of motion for fully saturated poroelastic halfspace and the embankment are solved by applying Fourier transform and Fourier series. Dynamic response of the ground with a permeable and impermeable surface induced by a moving train load is investigated and the effect of the hydraulic boundary condition on the vertical velocity and pore pressure in the ground is analyzed. In addition, a parametric analysis is performed for different train speeds, load frequencies, and soil permeability. The main conclusions of this paper are summarized as follows:(1)For moving trains with constant load amplitude, a critical velocity exists in the embankmentground structure, which is larger than the shear wave velocity of the ground. The vertical velocity of the ground with permeable surface is much larger than that with impermeable surface when the train speed approaches the critical speed.(2)Considering the train loads as moving harmonic loads, the vertical velocity of the ground reaches a peak value at an excitation frequency around 10 Hz. The effect of the hydraulic boundary of the ground surface on the vertical velocity response is limited at low load frequencies, while in the high frequency range the velocity response of the ground surface with impermeable boundary is significantly higher than that with permeable boundary.(3)The pore water pressure in the ground with impermeable surface is larger than that with permeable surface subjected to the moving train with a low selfexcitation frequency. As the load frequency increases, the permeable surface will result in higher pore water pressure in the ground than the impermeable surface.(4)The hydraulic boundary effects on the ground vibration are significant when the soil permeability is high, while its effects are reduced when the soil permeability is decreased.
Appendix
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work is supported by the National Key Research and Development Plan of China (Grant No. 2016YFC0800204), the National Natural Science Foundation of China (Grant Nos. 51578500 and 51708503), and the Natural Science Foundation of Zhejiang Province (Grant No. LY15E080009).
References
 G. Eason, “The stresses produced in a semiinfinite solid by a moving surface force,” International Journal of Engineering Science, vol. 2, no. 6, pp. 581–609, 1965. View at: Publisher Site  Google Scholar
 R. Kumar and P. Ailawalia, “Moving load response of micropolar elastic halfspace with voids,” Journal of Sound and Vibration, vol. 280, no. 35, pp. 837–848, 2005. View at: Publisher Site  Google Scholar
 Y. H. Huang, Pavement Analysis and Design, Prentice Hall, New Jersey, 1993.
 A. V. Metrikine and H. A. Dieterman, “Lateral vibrations of an axially compressed beam on an elastic halfspace due to a moving lateral load,” European Journal of Mechanics  A/Solids, vol. 18, no. 1, pp. 147–158, 1999. View at: Publisher Site  Google Scholar
 S. Kim, “Vibration and stability of axial loaded beams on elastic foundation under moving harmonic loads,” Engineering Structures, vol. 26, no. 1, pp. 95–105, 2004. View at: Publisher Site  Google Scholar
 S.M. Kim, “Influence of horizontal resistance at plate bottom on vibration of plates on elastic foundation under moving loads,” Engineering Structures, vol. 26, no. 4, pp. 519–529, 2004. View at: Publisher Site  Google Scholar
 S. Lu, “Analytical dynamic displacement response of rigid pavements to moving concentrated and line loads,” International Journal of Solids and Structures, vol. 43, no. 1415, pp. 4370–4383, 2006. View at: Publisher Site  Google Scholar
 X. Sheng, C. J. C. Jones, and D. J. Thompson, “A comparison of a theoretical model for quasistatically and dynamically induced environmental vibration from trains with measurements,” Journal of Sound and Vibration, vol. 267, no. 3, pp. 621–635, 2003. View at: Publisher Site  Google Scholar
 X. Sheng, C. J. C. Jones, and D. J. Thompson, “A theoretical study on the influence of the track on traininduced ground vibration,” Journal of Sound and Vibration, vol. 272, no. 3–5, pp. 909–936, 2004. View at: Publisher Site  Google Scholar
 X. Sheng, C. J. C. Jones, and D. J. Thompson, “A theoretical model for ground vibration from trains generated by vertical track irregularities,” Journal of Sound and Vibration, vol. 272, no. 3–5, pp. 937–965, 2004. View at: Publisher Site  Google Scholar
 B. Picoux, R. Rotinat, J. P. Regoin, and D. Le Houédec, “Prediction and measurements of vibrations from a railway track lying on a peaty ground,” Journal of Sound and Vibration, vol. 267, no. 3, pp. 575–589, 2003. View at: Publisher Site  Google Scholar
 M. A. Biot, “Theory of propagation of elastic waves in a fluidsaturated porous solid. I. Lowfrequency range,” The Journal of the Acoustical Society of America, vol. 28, pp. 168–178, 1956. View at: Publisher Site  Google Scholar  MathSciNet
 R. Siddharthan, Z. Zafir, and G. M. Norris, “Moving load response of layered soil. I: Formulation,” Journal of Engineering Mechanics, vol. 119, no. 10, pp. 2052–2071, 1993. View at: Publisher Site  Google Scholar
 D. D. Theodorakopoulos, A. P. Chassiakos, and D. E. Beskos, “Dynamic effects of moving load on a poroelastic soil medium by an approximate method,” International Journal of Solids and Structures, vol. 41, no. 7, pp. 1801–1822, 2004. View at: Publisher Site  Google Scholar
 D. D. Theodorakopoulos, “Dynamic analysis of a poroelastic halfplane soil medium under moving loads,” Soil Dynamics and Earthquake Engineering, vol. 23, no. 7, pp. 521–533, 2003. View at: Publisher Site  Google Scholar
 B. Jin, “Dynamic responses of a poroelastic half space generatedby high speed load,” Chinese Quarterly of Mechanics, vol. 25, no. 2, pp. 168–174, 2004. View at: Google Scholar
 G.B. Liu, P.C. Wang, Y.P. Chen, and H.L. Yao, “Ground vibration of a finite layer of soft foundation in the vicinity of a moving load,” Yantu Lixue/Rock and Soil Mechanics, vol. 27, no. 9, pp. 1607–1612, 2006. View at: Google Scholar
 B. Xu, J.F. Lu, and J.H. Wang, “Dynamic response of an infinite beam overlying a layered poroelastic halfspace to moving loads,” Journal of Sound and Vibration, vol. 306, no. 12, pp. 91–110, 2007. View at: Publisher Site  Google Scholar
 G. Y. Gao, J. F. He, C. B. Yang et al., “Ground vibration induced by trains moving on saturated ground using 2.5D FEM,” Chinese Journal of Geotechnical Engineering, vol. 33, no. 2, pp. 234–241, 2011. View at: Google Scholar
 F. Yu, Y.S. Diao, X.N. Tong, and Q.L. Zhang, “Damage identification of an offshore platform based on curvature of modal shape difference and BP neural network,” Journal of Vibration and Shock, vol. 30, no. 10, pp. 183–187, 2011. View at: Google Scholar
 K. Chahour, G. LefeuveMesgouez, and A. Mesgouez, “Spectral analysis of a railway track in contact with a multilayered poroviscoelastic soil subjected to a harmonic moving load,” Soil Dynamics and Earthquake Engineering, vol. 64, pp. 24–37, 2014. View at: Publisher Site  Google Scholar
 Y. Cai, Z. Cao, H. Sun, and C. Xu, “Effects of the dynamic wheelrail interaction on the ground vibration generated by a moving train,” International Journal of Solids and Structures, vol. 47, no. 17, pp. 2246–2259, 2010. View at: Publisher Site  Google Scholar
 Z. G. Cao, Y. Q. Cai, H. L. Sun, and C. J. Xu, “Dynamic responses of a poroelastic halfspace from moving trains caused by vertical track irregularities,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 35, no. 7, pp. 761–786, 2011. View at: Publisher Site  Google Scholar
 S. Zhou, C. He, and H. Di, “Dynamic 2.5D Green's function for a poroelastic halfspace,” Engineering Analysis with Boundary Elements, vol. 67, pp. 96–107, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 A. Karlström and A. Boström, “An analytical model for traininduced ground vibrations from railways,” Journal of Sound and Vibration, vol. 292, no. 12, pp. 221–241, 2006. View at: Publisher Site  Google Scholar
 C. Madshus and A. M. Kaynia, “Highspeed railway lines on soft ground: dynamic behaviour at critical train speed,” Journal of Sound and Vibration, vol. 231, no. 3, pp. 689–701, 2000. View at: Publisher Site  Google Scholar
 Z. Yuan, A. Boström, Y. Cai, and Z. Cao, “Analytical solution for calculating vibrations from twin circular tunnels,” Soil Dynamics and Earthquake Engineering, vol. 117, pp. 312–327, 2019. View at: Publisher Site  Google Scholar
 Z. Cao and A. Boström, “Dynamic response of a poroelastic halfspace to accelerating or decelerating trains,” Journal of Sound and Vibration, vol. 332, no. 11, pp. 2777–2794, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Zonghao Yuan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.