Research Article  Open Access
Study on Structural Service Performance of HeavyHaul Railway Tunnel with Voided Base
Abstract
The structural design of heavyhaul railway tunnels still follows the design method of ordinary railway tunnels. Most of them do not take into account the influence of large axle load of 30 t or more, let alone such problems as void of surrounding rock under longterm dynamic loads. In order to analyze the dynamic response of heavyhaul railway tunnels under longterm reciprocating cyclic dynamic loads, considering the factors such as axle load of vehicle body, unsprung mass, and track irregularity, the vibration load timehistory curve of heavyhaul railway trains is determined, the threedimensional dynamics coupling model of dynamic loadtunnelsurrounding rock is established, and the fatigue life of the structure under different void conditions is analyzed based on the SN curve of concrete. According to the study, the loading, unloading, and vibration caused by train passing will lead to fluctuations in the vertical displacement response of the monitoring point. The peaks and valleys of the response timehistory curve correspond to the effect of the train wheels rolling through. When the void is 6 m wide and 10 cm thick, the vertical displacement of the inverted arch is increased by about 9 times, the peak velocity of the inverted arch is increased by about 3.8 times, and the maximum principal stress is increased by about 47.3%, compared with the condition without void. With the same void thickness, the vertical displacement and velocity curves of the inverted arch vary significantly with the increase of the void width. The width of the base void has a significant effect on the fatigue life of the structure of heavyhaul railway tunnels. Based on the operation requirement of 100year service life, the ultimate void width is 2 m.
1. Introduction
Heavyhaul railways feature great transport capacity and significant economic and social benefits. They represent the level of a country’s freight capacity. Since the 1960s, heavyhaul railways have gradually dominated the railway transportation in various countries. Unlike traditional railways, the undertrack structure and surrounding rock of heavyhaul railways will bear greater loads due to the significant increase in the axle load of trains. Combined with a series of factors such as stability of soft surrounding rocks and track irregularity, heavyhaul railway tunnels are more easily damaged under longterm vibration of trains, and new requirements for the structure of heavyload railway tunnels are therefore put forward. The impact of heavyhaul trains on the base structure is mainly manifested in the increase of the load on the track and the undertrack base foundation. Such load is different not only from the static load, but also from the seismic load. It is a reciprocating cyclic load that is applied for a long time. The study on the mathematical expression and timehistory curve of the train excitation load is mainly based on field measurements and theoretical analysis. Jenkins et al. [1] established the dynamic analysis model of depressed joint track and defined the two objective special forces in wheelrail impact vibration: highfrequency impact force P2 and lowfrequency response force P1. Lamaran et al. [2] analyzed forced vibration of rail and simulated the train load as the sum of the static load and additional dynamic load associated with wheel weight, rail support condition, sleeper spacing, and vehicle speed. One et al. [3] carried out actual measurements and proved that the wheelrail acting force is severer in a medium frequency range and mainly affects the dynamic response of the vehicle body in a highfrequency range.
The dynamic response of tunnel base structure under such cyclic dynamic load is also a key to the scholars’ study [4, 5]. In the early days when the computational capacity was limited, a reasonable analytical model was the main research approach. Metrikine et al. [6] established a typical plane analytical model and analyzed the response of the structure under constant point load, covariant point load, and random point load by simulating the tunnel in viscoelastic soil with the Euler beam. Balendra et al. [7] established a semianalytical semiplane strain model and adopted the substructure method to reduce the complexity of the problem. Forrest et al. [8, 9] established a pipeinpine model and analyzed the vibration effect during train passing by taking the outer elastomeric pipe with infinite radius as the stratum, the thin inner pipe as the lining, and the floating plate inside the thin shell as the filling (the plate and the shell are connected by a spring). The model was later improved by Hussein et al. [10]. With the improvement of computer performance and computational capacity, more and more scholars have proposed numerical models suitable for train vibration. The finite element method was widely used to study the tunnel vibration response [11–13]. Because both finite element and boundary element have their own advantages in dealing with lining structure and infinite soil, Peng et al. [14, 15] carried out numerical simulation for the underlying concrete structure of tunnel under the impact load of train vibration with the nonlinear dynamic finite element theory and predicted the service life of the underlying structure by using the Tepfers concrete singlelogarithmic fatigue equation. A finite elementboundary element coupling model was established to analyze the impact of train vibration on adjacent buildings [16–18].
Under longterm vibration loads of heavyhaul trains, the tunnel base will inevitably be subjected to cumulative fatigue damage or even failure and void of surrounding rock. Due to the large discreteness of concrete and many other influence factors, there is no effective method to predict the fatigue life of concrete, and existing studies usually predict the fatigue life of concrete with the cumulative damage theory. Bazant et al. [19] and Mazars et al. [20] established the early scalar fatigue damage model for concrete with theories related to fracture mechanics. Fardis et al. [21] further optimized the model and adapted it to concrete strain softening by modifying the flexibility matrix equation of bounding surface. AbuLebdeh et al. [22] proposed a constitutive model considering the plasticity of concrete and the damage boundary. The model simulates the state of stiffness degradation and softening of concrete under the cyclic action multiaxial stress. Wu et al. [23] introduced tensile and shear damage variables to describe concrete stiffness degradation at a macroscale and derived the mathematical expressions of plastic strain evolution law and Helmholtz free energy. Based on the theoretical system of Helmholtz free energy, Peng et al. [24] conducted an experimental study for concrete of inverted arch under different vibration amplitudes in a coupling environment of static and cyclic loads and figured out the static load level that can cause concrete damage and failure as well as the quantitative expression of damage evolution trends at different stages; Wang et al. [25] conducted the impact compression test for many times on 5 types of support concrete of different ages with SHPB test equipment. The peak damage at stable age indicates that the statistical damage model of the Weibull distribution can well reflect the damage characteristics of the material.
During the service of heavyhaul railway tunnels, the fatigue damage and void of surrounding rock of the base structure are not independent of each other but are concomitant, and there are incentive and superposition effects between the two. Structural damage [26] will lead to decrease of stiffness and greater deformation under stress, which will incur greater load and compression on the surrounding rock of the base as well as continuous accumulation of plastic deformation [27, 28] and aggravated void of the surrounding rock, thus changing the support conditions of the structure and causing greater damage to the structure subjected to upper dynamic loads. This paper studies the characteristics of the fullcycle temporalspatial response of the base structure of heavyhaul railway tunnels under repeated impact by establishing a threedimensional numerical dynamics coupling model of dynamic loadtunnelsurrounding rock and predicts the fatigue life of the tunnel structure under different void conditions by carrying out failure evaluation for the base structure under damagevoid excitation based on the SN curve of concrete. It is found that during the service of heavyhaul railway tunnels, the presence of voided base has a significant influence on the displacement and velocity response of the lining structure as well as the distribution and magnitude of the principal stress on the lining. In addition, the void width of the base has a great influence on various indicators of the lining structure, while the influence of void thickness is minor.
2. Numerical Model of HeavyHaul Railway Tunnels
2.1. Determination of Train Excitation Load
The vibration of a train is under the joint influence of the train and the track. It is mainly caused by the abrasion of wheels, the state of rail joints, and geometric irregularity. On the vertical longitudinal plane of the wheelrail system, there are a wide variety of vibration sources: local unevenness of the existing rail surface, such as rail head crushing or rail surface peeling, and periodic system irregularities, such as waveform track, corrugation, and eccentric wheels. Particularly, abrupt irregularities are more common at such places as depressed joints, joints with height difference, large rail joints, and welds. In addition, defects such as loose sleepers or hardened ballast bed may occur to the undertrack foundation, resulting in dynamic irregularity of uneven elasticity, and thereby creating an excitation load.
Considering the influence factors of train vibration and combined with the management value of track geometric irregularity in the UK, the vertical excitation force of the train can be simulated by an excitation function [1, 29]. This function takes into account the axle load of vehicle body, unsprung mass, track irregularity, driving speed, load periodicity, movement and combination of loads, and the effect of rails and sleepers on the dispersion and transmission of loads, as shown in the following formula [30]:where is the static load of wheels and , , and are the vibration loads corresponding to the typical values of ①–③ in Table 1. Taking as the unsprung mass of the train, then the corresponding vibration load should bewhere is the typical rise of arch and is the circular frequency of irregularity vibration wavelength at corresponding speed, corresponding to ①, ②, and ③ in Table 1, respectively. The formula iswhere is the running speed of train and is the typical wavelength corresponding to ①, ②, and ③ in Table 1.

The management value of track geometric irregularity in the UK is shown in Table 1 [31].
When calculating the train load, the unilateral static wheel weight should be taken: the vehicle load is 150_{ }kN and the unsprung mass is taken as M_{0 }=_{ }1200_{ }kg. Generally, the speed of heavyhaul railways does not exceed 100_{ }km/h. Therefore, the vibration load of the train at a speed _{ }=_{ }100_{ }km/h is taken into account for the simulation of the timehistory curve of train excitation load, as shown in Figure 1.
2.2. Material Mechanics Parameters and Boundary Condition
2.2.1. Material Mechanics Parameters
During dynamic analysis, the rock mass medium is taken as an ideal elasticplastic material which follows the Mohr–Coulomb strength criterion, and the largestrain deformation mode is adopted. The yield function iswhere , , and are the principal stresses of the material (, and the compressive stress is negative); is the cohesive force; is the tensile strength of the material; and is related to the internal friction angle, and its expression is
Shear damage occurs when the stress on a certain point within the rock mass is positive (); tensile damage occurs when .
The potential function in the Mohr–Coulomb model is described by and , defining the shear plastic flow and tensile plastic flow, respectively. Follows the nonassociated flow rule, while follows the associated flow rule:where is the dilatancy angle and .
The function h is defined in the stress plane:where and .
The positive and negative zones indicated by the function h are taken as the criteria for determining the tensile and shear plastic flow zones, as shown in Figure 2.
If the stress point falls into zone 1, then the shear yield occurs, and the stress should be adjusted to the curve corresponding to the yield function according to the flow rule determined by the shear potential function ; if the stress point falls into zone 2, then the tensile yield occurs, and the stress should be adjusted to the curve corresponding to the yield function according to the flow rule determined by the tensile potential function .
Composite lining is used for the tunnel support structure. C25 shotcrete is used for the primary support, C30 form concrete is used for the secondary lining, and the surrounding rocks are of redbed soft rock. See Table 2 for the specific material parameters of the model.

2.2.2. Dynamic Boundary Condition
Viscous boundary, which was first proposed by Lysmer and Kuhlemeyer [32], is adopted as the boundary condition of the model. This condition is equivalent to the situation where independent dampers are separately set at the normal and tangential positions of the artificial boundary to absorb the scattered wave by dissipating the energy of the vibration wave with the generated viscous damping force proportional to the moving speed. This method is easy to be implemented in dynamic timedomain analysis, and its validity has been widely proved in finite element and finite difference methods. The normal and tangential forces provided by the dampers arewhere is the density of rock mass media; and are the wave velocities of longitudinal and transverse waves in rock mass media, respectively; and and are the normal and tangential components of the particle vibration velocity due to fluctuations on the artificial boundary.
In the calculation model, a viscous boundary is applied to the bottom and the left and right boundaries, and normal constraint is imposed on the front and rear boundaries. The soil mass above the model is the free surface.
The depth of the tunnel is taken as 55.97 m, and the height and width of the model are 100 m. With reference to the wheelrail excitation force function provided in Formula (1), the duration of train excitation load applied on the model is 3.5 s. In the calculation of the train excitation load applied on the model, the concentrated force is converted into surface force loaded on the subgrade to simulate the train passing through the tunnel. See Figure 3 for the calculation model.
2.3. Analysis of HeavyHaul Railway Tunnel with Voided Base
The structure of the tunnel base is in an initial state at the beginning. With the increased times of train load, the surrounding rock of the base is gradually voided. This may occur symmetrically from the position near the track (void condition 1) or from the center of the inverted arch (void condition 2). With the occurrence of the void, the base structure of the tunnel gradually cracks from elastic deformation. With the occurrence of cracks, combined with the weakened constraint on the structure caused by void, the base structure will be subjected to greater deformation under train load, and such deformation will accelerate the formation of void in turn (void condition 3). Due to continuous development of void, the support condition of the base structure keeps changing and the number of cracks keeps increasing. Such structural damage and void of base develop with and promote each other, finally resulting in failure of base structure. See Figure 4 for the evolution concept.
In order to analyze the influence of the base void (Figure 5) on the dynamic response of the base structure, the calculation model compares and analyzes 8 void cases, as shown in Table 3. The dynamic response of the structure caused by the train load is mainly at the inverted arch. Therefore, monitoring points are arranged near the inverted arch, the arch springing, and the void area of the tunnel, as shown in Figure 6.

2.4. Prediction Methods for Fatigue Life of Concrete
Fatigue refers to the structural change process in which the local permanent and progressive development of the material occurs under the condition that alternating stress and strain are applied on a certain point or at some points of the material. Fatigue failure is the result that the fatigue damage is accumulated to a critical value. It may cause the accumulation of cracks or the bending fracture of the material after certain times of alternating loads are applied [33]. The fatigue life refers to the cycle times that the material is applied with stress or strain before the fatigue failure appears, or the time between the beginning of the loading and the occurrence of the fracture.
Under the action of the cyclic load, the material damage of the base structure of heavyhaul railway tunnels is more and more serious, the material properties deteriorate gradually with the increase of cycle times, and thus the bearing capacity of the components decreases gradually and finally damage appears. The cumulative fatigue damage model of concrete can be established to predict the service life. The prediction can provide a reference for the evaluation of the residual service life of the base structure of the existing heavyhaul railway tunnel and the maintenance and reinforcement of the damaged structure. The existing methods for fatigue life prediction of concrete mainly include nominal stress method, energy method, local stress and strain method, and stress field intensity method. The most widely used method is the nominal stress method, which is also known as SN curve method. This method is mainly applicable to the calculation of highcycle fatigue life. The specific analytical steps of the nominal stress method are as follows:(1)Determine the load spectrum and apply it to the structure to obtain the nominal stress spectrum of the structure(2)Determine the fatigue dangerous part of the structure and its nominal stress spectrum(3)Establish the SN curve equation of the structure or component(4)Determine the fatigue life of the structure
3. Results and Discussion
3.1. Effects on Displacement and Velocity of Base under Void Condition
The displacement and velocity response of tunnel structure under train load is mainly vertical, and the response of the inverted arch is the most intense. It is found that the displacement and velocity of the monitoring point A in the center of the inverted arch are the largest through monitoring the vertical displacement and velocity of each monitoring point. The point A is taken as the research object, the vertical displacement and velocity timehistory curves under different void dimensions are extracted, as shown in Figures 7 and 8.
(a)
(b)
(a)
(b)
Through the analysis on vertical displacement and velocity timehistory curves for the inverted arch of the tunnel, the findings are as follows.
The effects of loading, unloading, and vibration generated by train passing cause the fluctuation of vertical displacement response at lining monitoring points. The peaks and valleys of the timehistory curve for each response correspond to the rolling through effect of train wheels.
The vertical displacement and velocity curves of the inverted arch vary obviously with the increase of void width, which shows that the void width of the base under the rail has a great effect on the displacement and velocity of the inverted arch. However, such variation is minor with the increase of the void thickness, which shows that the void thickness of the base under the rail has little effect on the displacement of the inverted arch. The peak values of vertical displacement and velocity of the inverted arch in different cases are shown in Table 4.

The curves for relationship between the peak values of displacement and velocity of the inverted arch and the void area of the base are shown in Figures 9 and 10.
(a)
(b)
(a)
(b)
The peak value of vertical displacement of the inverted arch reaches 26 mm, and the vertical velocity reaches 8.2 cm/s when the void width of base reaches 6 m. The vertical displacement is increased by about 9 times, and the peak value of velocity is increased by about 3.8 times, compared with the condition without void. Thus, it can be seen that the vertical displacement and velocity of the inverted arch are obviously affected by the void of base.
The curve variation trends in Figures 9 and 10 also prove that the void depth has a greater effect on the displacement and velocity responses of the base structure than the void thickness. The peak value variations for vertical displacement and velocity of inverted arch are around 1 mm and 0.4 cm/s, respectively, with the increase of void thickness. In order to verify the rationality of this result, the peak value of vertical displacement and velocity of the inverted arch is calculated in the case of void width of 12 m since the void width is almost impossible to be greater than 12 m, and the calculation results are shown in Table 5.

It can be seen from Table 5 that, in the case of a void width of 12 m, displacement and velocity remain almost unchanged. Therefore, the void thickness has little effect on the peak value of displacement and velocity of the inverted arch when the void width is in the range between 5 m and 12 m.
3.2. Effects on Maximum Principal Stress of Base under Void Condition
The curves for maximum principal stress at the inverted arch in different cases are shown in Figure 11.
(a)
(b)
As shown in Figure 11, the response curve of maximum principal stress of the inverted arch varies greatly with the increase of void width, which shows that the void width of the base under the rail has a great effect on the principal stress of the inverted arch. However, the curve varies slightly with the increase of the void thickness, which shows that the void thickness of the base under the rail has little effect on the principal stress of the inverted arch.
The peak values of the maximum principal stress of the inverted arch in different cases are shown in Table 6. The variation of the peak value corresponding to the expansion of the void area of base is shown in Figure 12.

(a)
(b)
It can be seen from Table 6 and Figure 12 that the maximum principal stress of the inverted arch is affected obviously by the void of base under the rail and that the maximum principal stress of the inverted arch varies obviously with the increase of the void width. The peak value of maximum principal stress of the inverted arch reaches 1.65 MPa when the void width of the base reaches 6 m, increased by 0.53 MPa compared with the condition without void. Thus, it can be seen that the maximum principal stress of the inverted arch is affected obviously by the void width of the base. The peak value variation for maximum principal stress of the inverted arch is around 0.02 MPa with the gradual increase of void thickness. Like the displacement and velocity, the maximum principal stress of the inverted arch in the case of void width of 12 m is also obtained in Table 7.

Table 7 indicates that when the void width is 12 m, the peak value of maximum principal stress of the inverted arch is approximately 2.4 MPa regardless of the void thickness. Thus, it can be seen that void thickness has little effect on the peak value of principal stress of the inverted arch when the void width is in the range between 5 m and 12 m.
3.3. Analysis on Fatigue Life of HeavyHaul Railway Tunnel Structure under Train Vibration
It can be concluded through the analysis in the previous sections that the heavyhaul railway tunnel structure has the most obvious dynamic response at the inverted arch. Meanwhile, based on the analysis of the maximum principal stress, the structure is subjected to tensile stress at the inverted arch only, and it is subjected to compressive stress at other parts. The level of compressive stress is small for concrete, and it is generally considered that it will not cause fatigue failure. Therefore, the part that ultimately controls the fatigue life of heavyhaul railway tunnel structure under train load is the part where the inverted arch is subjected to tensile stress, and the maximum stress of this part has been determined.
The most important step in the nominal stress method is to determine the fatigue life SN curve of the structure or component. At present, the fatigue performance data of all kinds of components are obtained by a large number of tests so that the SN curve is determined. Many scholars have studied the fatigue performance of concrete material in tension. Tepfers et al. [34] put forward an equation similar to that of concrete compressive fatigue to describe concrete tensile fatigue on the basis of fatigue tests such as axial tension and splitting, as shown in the following formula:where , ; : upper limit value of fatigue tensile stress; : lower limit value of fatigue tensile stress; : static axial tensile strength of concrete; a_{t}, b_{t}: coefficients determined by the fatigue test, which are generally taken as a_{t} = 1.0 and b_{t} = 0.0685; and N: limit loading times.
The above equation is more accurate in describing the fatigue life of concrete under axial tension. But the base of heavyhaul railway tunnels is a structure subjected to bending tension, and its stress condition is more complicated than a structure subjected to axial tension. Sidoroff et al. [35] have proved that the above formula is also applicable to the bending tensile fatigue SN curve of concrete, but the taken value of the coefficient in the fatigue test is different. The bending tensile fatigue equation of concrete is obtained as follows [36, 37]:where ; : maximum tensile stress of concrete at lower edge under constantamplitude repeated load; : lower limit value of fatigue tensile stress; : static bending tensile strength of concrete; a_{r}, b_{r}: coefficients determined by the fatigue test, which are generally taken as a_{r} = 0.94–1.35 and b_{r} = 0.045–0.118.
The timehistory curves for the maximum principal stress of the inverted arch under different void conditions have been obtained through the calculation and analysis for the dynamic response of the structure under different void conditions. Table 8 shows the statistical results of the maximum and minimum values.

The maximum and minimum values of the maximum principal stress are, respectively, regarded as the maximum and the minimum values of the tensile stress of the structure. The bending tensile fatigue life of the concrete is calculated with the formula (10), where a_{r} = 0.94 and b_{r} = 0.049. Raphael et al. [38] has proved that the bending tensile strength of concrete is generally 35% higher than the straight tensile strength based on many test results. Therefore, the fatigue life of concrete is calculated by taking f_{r} = 2.2 MPa and setting Case 1 as the example:
The load action times corresponding to the fatigue failure:
Thus, it can be seen from the above calculation that the fatigue life of the structure under the condition without void is longer.
Similarly, the calculation results of the rest of the cases are statistically shown in Table 9. The variation of fatigue life in different cases is shown in Figure 13.

(a)
(b)
It can be seen from Table 9 and Figure 13 that the void width of the base has an obvious effect on the fatigue life of heavyhaul railway tunnel structure, and the fatigue life of the structure decreases sharply especially in the early stage of void development. The variation of fatigue life tends to be gentle when the void width reaches a certain value. It can be seen from Case 3 that the action times of fatigue load for the structure is only 7.08 × 10^{6} when the void reaches 3 m. If it is still estimated as per that the wheel set load of the Datong–Qinhuangdao Railway is 1,000 times per day, the service life of the tunnel is only 19.40 years. The allowable ultimate void width of base is 2 m based on the operation requirement of 100year service life.
The void thickness of the base will not have an obvious effect on the service life of the heavyhaul railway tunnel structure with certain void width, and the fatigue life of the structure is maintained at a longer level.
4. Conclusions
In this paper, the dynamic response characteristics of the heavyhaul railway tunnel structure under repeated impacts are studied based on the threedimensional numerical dynamics coupling model of dynamic loadtunnelsurrounding rock, and the fatigue life of the tunnel structure under different void conditions is predicted based on the SN curve of concrete. The main conclusions are as follows:(1)The effects of loading, unloading, and vibration generated by train passing cause the fluctuation of vertical displacement response at lining monitoring points. The peaks and valleys of the timehistory curve for each response correspond to the rolling through effect of train wheels.(2)The dynamic response of base with void is much larger than that without void: When the void width is 6 m and the void thickness is 10 cm, the vertical displacement of the inverted arch is about 10 times that under the condition without void, the peak velocity is about 4.9 times that under the condition without void, and the maximum principal stress is increased by about 47.3%.(3)The dynamic response of base is greatly affected by the void width under the condition of base with void. With the same void thickness, when the void width is 6 m, the peak values of displacement and velocity of the inverted arch are, respectively, about 4.26 and 3.90 times the values when the void width is 2 m.(4)The void width of base has an obvious effect on the fatigue life of the heavyhaul railway tunnel structure. The ultimate void width is 2 m according to the SN curve of concrete under bending tensile condition and based on the operation requirement of 100year service life.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interests regarding the publication of this paper.
Acknowledgments
The authors are very grateful for the funding from the National Key Research and Development Program of China (2016YFC0802202) and the General Program of National Natural Science Foundation of China (51678503).
References
 H. H. Jenkins, J. E. Stephenson, G. A. Clayton, G. W. Morland, and D. Lyon, “The effect of track and vehicle parameters on wheel/rail vertical dynamic forces,” Railway Engineering Journal, vol. 3, no. 1, pp. 2–16, 1974. View at: Google Scholar
 G. Lamaran and M. Derdas, “Evaluation of dynamics load on rail track sleepers based on vehicletrack modeling and analysis,” International Journal of Structural Stability, vol. 2, no. 3, 2002. View at: Google Scholar
 K. One and M. Yamada, “Analysis of railway track vibration,” Journal of Sound and Vibration, vol. 130, no. 2, pp. 269–297, 1989. View at: Google Scholar
 J. Lai, S. He, J. Qiu et al., “Characteristics of seismic disasters and aseismic measures of tunnels in Wenchuan earthquake,” Environmental Earth Sciences, vol. 76, no. 94, 2017. View at: Publisher Site  Google Scholar
 J. Qiu, X. Wang, and J. Lai, “Response characteristics and preventions for seismic subsidence of loess in Northwest China,” Natural Hazards, vol. 92, no. 3, pp. 1909–1935, 2018. View at: Publisher Site  Google Scholar
 A. V. Metrikine and A. C. W. M. Vrouwenvelder, “Surface ground vibration due to moving train in a tunnel: twodimensional model,” Journal of Sound and Vibration, vol. 234, no. 1, pp. 43–66, 2000. View at: Publisher Site  Google Scholar
 T. Balendra, C. G. Koh, and Y. C. Ho, “Dynamic response of buildings due to trains in underground tunnels,” Earthquake Engineering and Structural Dynamics, vol. 20, no. 3, pp. 275–291, 1991. View at: Publisher Site  Google Scholar
 J. A. Forrest and H. E. M. Hunt, “A threedimensional model for calculation of traininduced ground vibration,” Journal of Sound and Vibration, vol. 294, no. 1, pp. 678–705, 2006. View at: Publisher Site  Google Scholar
 J. A. Forrest and H. E. M. Hunt, “Ground vibration generated by trains in underground tunnels,” Journal of Sound and Vibration, vol. 294, no. 4, pp. 706–736, 2006. View at: Publisher Site  Google Scholar
 M. F. M. Hussein and H. E. M. Hunt, “A numerical model for calculating vibration from a railway tunnel embedded in a fullspace,” Journal of Sound and Vibration, vol. 305, no. 3, pp. 401–431, 2007. View at: Publisher Site  Google Scholar
 R. Thiede and H. G. Natke, “The influence of thickness variation of subway walls on the vibration emission generated by subway traffic,” in Proceedings of Soil Dynamics and Earthquake Engineering V: International Conference Soil Dynamics and Earthquake Engineering, Karlsruhe, Germany, September 1991. View at: Google Scholar
 F. Guan and I. D. Moore, “Threedimensional dynamic response of twin cavities due to traveling loads,” Journal of Engineering Mechanics, vol. 120, no. 3, pp. 637–651, 1994. View at: Publisher Site  Google Scholar
 K. H. Chua, K. W. Lo, and T. Balendra, “Building response due to subway train traffic,” Journal of Geotechnical Engineering, vol. 121, no. 11, pp. 747–754, 1995. View at: Publisher Site  Google Scholar
 L. M. Peng, C. H. Shi, J. Huang, and S. L. Liu, “Study on the fatigue life of the tunnel bed structure under train loads,” Journal of the China Railway Society, no. 1, pp. 82–85, 2007. View at: Google Scholar
 L. M. Peng, L. C. Huang, and S. L. Liu, “Study on dynamic response of railway tunnel bedding structure,” Journal of Vibration and Shock, vol. 25, no. 5, pp. 175–179, 2006. View at: Google Scholar
 D. Clouteau, M. Arnst, T. M. G. AlHussaini et al., “Freefield vibrations due to dynamic loading on a tunnel embedded in a stratified medium,” Journal of Sound and Vibration, vol. 283, no. 1, pp. 173–199, 2005. View at: Publisher Site  Google Scholar
 Z. Zhang, X. Shi, X. Wang, and H. Li, “Stability of NATM tunnel faces in soft surrounding rocks,” Computers and Geotechnics, vol. 96, pp. 90–102, 2018. View at: Publisher Site  Google Scholar
 Z. Zhang, H. Li, H. Yang, and B. Wang, “Failure modes and face instability of shallow tunnels under soft grounds,” International Journal of Damage Mechanics, Article ID 105678951877313, 2018. View at: Publisher Site  Google Scholar
 Z. P. Bazant and H. O. Byung, “Strainrate effect in rapid triaxial loading of concrete,” ACI Materials Journal, vol. 180, pp. 764–782, 1982. View at: Google Scholar
 J. Mazars, “Description of micro and macroscale damage of concrete structures,” Engineering Fracture Mechanics, vol. 25, no. 56, pp. 729–737, 1985. View at: Google Scholar
 M. N. Fardis and E. S. Chen, “A cyclic multiaxial model for concrete,” Computational Mechanics, vol. 1, no. 4, pp. 301–315, 1986. View at: Google Scholar
 T. M. AbuLebdeh, “Plasticitydamage model for concrete under cyclic multiaxial loading,” Journal of Engineering Mechanics, ASCE, vol. 119, no. 7, pp. 1465–1485, 1993. View at: Google Scholar
 J. Y. Wu, J. L. Li, and F. Rui, “An energy release ratebased plasticdamage model for concrete,” International Journal of Solids and Structures, vol. 43, no. 3, pp. 583–612, 2006. View at: Publisher Site  Google Scholar
 L. M. Peng, N. Liu, and C. H. Shi, “Experimental study on dynamic amplitudes for cumulate damage characteristics of tunnel invert concrete,” Journal of Railway Science and Engineering, vol. 13, no. 6, pp. 1091–1099, 2016. View at: Google Scholar
 S. M. Wang, X. B. Li, F. Q. Gong, and J. J. Zhu, “Experimental study on mechanical properties of different ages concrete under static and dynamic load,” Engineering Mechanics, vol. 30, no. 2, pp. 143–149, 2013. View at: Google Scholar
 Q. Yan, H. Chen, W. Chen et al., “Dynamic characteristic and fatigue accumulative damage of a cross shield tunnel structure under vibration load,” Shock and Vibration, vol. 2018, pp. 1–14, 2004. View at: Google Scholar
 J. Qiu, X. Wang, S. He et al., “The catastrophic landside in Maoxian county, Sichuan, SW China, on June 24, 2017,” Natural Hazards, vol. 89, no. 3, pp. 1485–1493, 2017. View at: Publisher Site  Google Scholar
 J. Qiu, Y. Xie, H. Fan, Z. Wang, and Y. Zhang, “Centrifuge modelling of twintunnelling induced ground movements in loess strata,” Arabian Journal of Geosciences, vol. 10, no. 493, pp. 1–14, 2017. View at: Publisher Site  Google Scholar
 T. M. Dawn and C. G. Stanworth, “Ground vibrations from passing trains,” Journal of Sound and Vibration, vol. 66, no. 3, pp. 355–362, 1979. View at: Publisher Site  Google Scholar
 C. Pan and G. N. Pande, “Preliminary deterministic finite element study on a tunnel driven in loess subjected to train loading,” China Civil Engineering Journal, vol. 17, no. 4, pp. 19–28, 1984. View at: Google Scholar
 B. Liang, The Dynamic Characteristics of High Speed Railway Subgrade and Geosynthetics’ Applications, Southwest Jiaotong University, Chengdu, China, 1998, PhD Dissertation.
 J. Lysmer and R. L. Kuhlemeyer, “Finite dynamic model for infinite media,” Journal of the Engineering Mechanics, ASCE, vol. 95, no. 4, pp. 859–877, 1969. View at: Google Scholar
 J. M. Hanson, C. A. Ballinger, and D. Linger, “Considerations for design of concrete structures, subjected to fatigue loading,” ACI Structural Journal, vol. 71, 1974. View at: Google Scholar
 R. Tepfers, “Tensile fatigue strength of plain concrete,” Journal of the American Concrete Institute, vol. 76, no. 8, pp. 919–933, 1979. View at: Google Scholar
 F. Sidoroff, Description of Anisotropic Damage Application to Elasticity, Springer, Heidelberg, Germany, 1981.
 J. W. Murdock and C. E. Kesler, “Effect of range of stress on fatigue strength of plain concrete beams,” Journal of American Concrete Institute, vol. 55, no. 2, pp. 221–232, 1959. View at: Google Scholar
 B. H. Oh, “Fatigue analysis of plain concrete in flexure,” Journal of Structural Engineering, vol. 112, no. 2, pp. 273–288, 1986. View at: Publisher Site  Google Scholar
 J. M. Raphael, “Tensile strength of concrete,” Journal of the American Concrete Institute, vol. 81, no. 2, pp. 158–165, 1984. View at: Google Scholar
Copyright
Copyright © 2018 Zhiqiang Zhang 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.