Abstract

Dynamic loads such as the train vibration load usually act on the shield tunnel lining in the long term, which could make the initial flaws in shield segment propagate and gradually weaken the robustness of the tunnel structure. In this paper, a three-dimensional numerical model of shield tunnel lining structure with the initial defect is built to study its dynamic reaction and fatigue crack propagation under the train vibration load. Furthermore, the damage to intact shield segment caused by train vibration load is studied by employing the rain-flow counting method and the Miner damage theory, and a rational fatigue life estimation for the concrete shield tunnel lining is finally made. Results show that crack propagation is influenced by both the train speed and train axle, the higher the train speed, the longer the final crack, and train axle has a larger influence than train speed on the crack propagation in shield tunnel segment. The shield tunnel lining structure of Nanjing Metro Line 5 can meet the demand of working for a hundred years under the current working conditions.

1. Introduction

As one kind of the public transportation, subway has played a more and more important role in the contemporary society nowadays, especially in the modern metropolises. The main structure of the subway system is the so-called shield tunnel lining structure, which usually consists of plenty of concrete shield segments. The shield segments are manufactured in the factory and transported to the construction site in a long distance, and they are usually assembled one by one in situ; thus, some initial flaws such as cracks may come into being in these shield segments. Furthermore, the lining structure of the shield tunnel is under train vibration load day and night, making these flaws propagate ceaselessly. All of these produce a large threat to the shield tunnel system. Moreover, even if the shield tunnel lining is intact without any defect, it has a significant and realistic value to make a rational life estimation of the shield tunnel lining.

During the past several decades, numerous numerical and experimental research studies have been conducted on shield tunnel structures. The Railway Engineering Society of the United Kingdom [1] studied the vertical wheel-rail force of the train by both theoretical analysis and experiment in 1970s. Newton and Clark [2] compared the Timoshenko beam model with the Euler beam model and concluded that the predicted wheel-rail forces of the train using these two models are close to each other, while the Euler beam model is much simpler than the other one. Yang and Hung [35] proposed a simplified 2.5D finite element-infinite element coupling model, using finite element to simulate the near field and infinite element to simulate the far field, and obtained quite good results. Sheng et al. [6, 7] proposed a semispace numerical method for a circular tunnel based on the wavenumber finite and boundary element method. Degrande et al. [8] developed a finite element-boundary element coupling method to calculate subway vibration and proposed a highly efficient numerical model. Forrest and Hunt [9, 10] proposed a pipe-in-pipe model to study the train-induced ground vibration in underground tunnel, and the computational efficiency of this model is rather good. Seleznev et al. [11, 12] developed a method to study an elastic half-space containing a cylindrical tunnel. Gardien and Stuit [13] proposed a model consisting of three submodels to analyse the tunnel vibration induced by train and to reduce the numerical burden using finite element analysis. Other research studies on tunnel vibration included the application of Green’s functions [14] and three-dimensional finite element simulations of long tunnels in half-space [1517]. For the studies of practical projects in engineering, Lai et al. [18] conducted a study concerned with the assessment of the vibration induced by the passage of commuter trains running in a tunnel placed underground of the city of Rome, using both the numerical simulation and field monitoring. Huang et al. [19] studied the dynamic characteristics of the tunnel invert for high-speed railways. Connolly et al. [20] analysed 1500 records of ground vibration from 17 high-speed railway stations in Europe and found that the soil property was the main factor which had an influence on the ground vibration. Yi et al. [21] and Wei et al. [22] studied the dynamic reactions of Guangzhou Metro Line by both the indoor model tests and numerical simulations. Besides, a series of numerical simulations on damage and fatigue analysis of some representative structures in civil engineering have been conducted [2327], and lots of useful results have been obtained for practical engineering.

The shield tunnel lining segment is usually precast with concrete in factories. Concrete is a kind of man-made material which is brittle and sometimes fragile. Lots of studies have been conducted on the fracture and fatigue properties of concrete in the past century [2832]. Besides, Shah [33] acquired the fracture parameters for the concrete fracture model using the three-point bending test. Kazemi [34] proposed the process zone model and applied it to the fracture analysis of concrete structures. Karihaloo and Nallathambi [35] found that the fracture toughness of concrete was not determined by the shape or dimension of the specimen but by the size effect model with mixed variables. Xu et al. [36, 37] proposed the double-K fracture criterion for concrete, which has been widely used for the fracture analysis of concrete structures since then. By employing the double-K fracture criterion of concrete, Tian et al. [38, 39] conducted prototype load tests and three-dimensional numerical simulations to investigate the fracture property and failure model of concrete shield segments.

In this paper, a three-dimensional numerical model of the shield tunnel lining structure based on Nanjing Metro Line 5 is built in ABAQUS to study the structure reaction and fatigue crack growth under train vibration load. Besides, damage to the shield segment induced by the train vibration load is also calculated by employing the Miner damage theory and S-N fatigue law of concrete, and the fatigue life of the concrete shield lining is rationally estimated.

2. Numerical Model and Train Vibration Load

2.1. Numerical Model of the Shield Lining

The numerical model of the shield tunnel lining is built according to the real shield tunnel of Nanjing Metro Line 5 in China. Thickness of the segment is 350 mm, and the inner diameter of the tunnel lining is 5500 mm. The intact lining ring consists of six components, which include one sealing-roof segment, two adjacent segments, and three standard segments. The joints of the adjacent segments are connected by M30 bolts. The cross section of the tunnel lining ring is illustrated in Figure 1, in which A denotes the standard segment, B denotes the adjacent segment, and K denotes the sealing-roof segment. Due to the complexity of rebars in the concrete shield tunnel segments, they are not considered in the current numerical simulation.

The numerical models of the shield tunnel are built in FEM software ABAQUS, which are shown in Figures 24. Note that there are circumferential bolts and longitudinal bolts in the numerical model of the shield tunnel lining, as shown in Figure 4, and they are “tied” with the concrete in the numerical simulation.

Concrete grade of this tunnel segment is C50, and its material properties are listed in Table 1.

2.2. The Train Vibration Load

The train vibration load of the subway is a very complex variable which is influenced by many factors, such as the train axle, speed of the train, rail irregularity, and so on. The best way to obtain train vibration load is to conduct in situ monitoring, despite the fact that it usually costs plenty of human labour as well as economic burden. In this numerical simulation, the empirical formula recommended by the Railway Engineering Society of the United Kingdom is employed to calculate the train vibration [1], which is expressed as follows:where denotes the static load of the wheel axel and denote three vibration loads with different frequencies.

The vibration loads can be expressed as follows:where denotes the mass of the train, denotes the arch rise, and denotes the circular frequency of the vibration wave. The circular frequency is expressed aswhere denotes the speed of the train and denotes the wavelength.

The related parameters for the aforementioned formulas are listed in Table 2.

The axle for the train of Nanjing Metro Line 5 is 16 t, and the speed of the train for the numerical simulation is 80 km/h. The wavelength of the carriage smoothness is 10 m, and the versine is 5 mm. The wavelength for waveform wearning is 0.5 m, and its versine is 0.005 mm. The time-history curve of the train vibration load is depicted in Figure 5. The aforementioned train vibration loads are applied at the position of wheel track at the bottom segment of the tunnel lining, as described by equation (1), and their time-history curve is illustrated in Figure 5.

3. Crack Propagation

3.1. The Fatigue Criterion

Lots of fatigue laws have been proposed and applied for the fatigue analysis of various kinds of materials in engineering during the past several decades; among them, the Paris law [40] is the most widely used one due to its simplicity. The general form of Paris law is usually expressed as follows:where a denotes the crack length (mm); N denotes the fatigue load cycles; denotes the crack increment during each load cycle; ΔK denotes the stress intensity factor range, with the unit of ; and C and n are constants related to the material properties.

The Paris fatigue law is often used for metals, but research studies on concrete [41] show that it could be applied to concrete material as well, with the following expression:where denotes the initial crack depth and D denotes thickness of the concrete specimen.

As the stress intensity factor can depict the stress field near the crack tip, it is the most important variable for fracture and fatigue analysis of concrete structures. The double-K fracture criterion [37] is employed in this study for crack propagation simulation of the concrete shield tunnel ling. For concrete C50, its measured initial fracture toughness is , and the ultimate fracture toughness is .

3.2. The Crack Propagation Path

In this section, the numerical simulation is conducted on the standard segment. The dimension for the initial crack is 300 mm × 25 mm, and the number of load cycles is set to 108. The train axle weighs 16 t, and the standard speed of the train is 80 km/h.

The numerical model with an initial crack is shown in Figure 6, in which the area with red colour represents the initial crack. The final crack is shown in Figure 7 under the cyclic train vibration load. Since the initial crack is so small, only the segment which contains the crack is shown in Figures 69.

3.3. The Influence of Train Speed

The influence of train speed on the crack propagation is studied in this section. The dimension for the initial crack is 300 mm × 25 mm, and the number of load cycles is set to 108. The train axle weighs 16 t, and the train speeds are set to four levels, which include 100 km/h, 80 km/h, 60 km/h, and 40 km/h. The final crack propagation path under four different train speeds is shown in Figure 8.

As shown in Figure 8, crack propagation shows the same path overall under different train speeds, while the final crack lengths are a little different. The higher the train speed, the longer the final crack.

3.4. The Influence of  Train Axle

The influence of train axle on the crack propagation is studied in this section. Dimension of the initial crack and number of load cycles are the same as those in Section 3.3. The train speed is fixed at 80 km/h, and the train axles are set to four levels, which include 18 t, 16 t, 14 t, and 12 t. The final crack propagation path under four different train axles is shown in Figure 9.

As shown in Figure 9, crack propagation shows the same path overall under different train axles, while the final crack lengths are very different. The final crack under train axle of 18 t is so long that it almost cuts through the shield segment, while the final crack under train axle of 12 t is quite short and only develops to a local crack.

4. Fatigue Life Estimation

4.1. The Rain-Flow Counting Method

The rain-flow counting method [42] is widely used to deal with the time-history of load in engineering, which could account for the memory characteristics of the material. In the rain-flow counting method, the time-history of load is decomposed into several loading cycles which include the loading amplitude and loading mean value, in order to make a rational fatigue life estimation of the structure based on the fatigue law of the material.

Figure 10 shows a schematic diagram for the rain-flow counting method. The rain-flow persistently flows from top to bottom along the slope, until it encounters a larger peak value or the upper rain-flow. The fundamental principle of the rain-flow counting method is illustrated by the following figure.

4.2. Miner Damage Theory and S-N Equation

The Miner damage theory assumes that if the structure is subjected to different cyclic stress , the damage caused by cyclic loading is . Structure fatigue failure occurs when the cumulative damage of the structure under various stress levels reaches 1. That is to say, to prevent the fatigue failure of the structure, the following equation must be satisfied:where DTotal is the total damage of the structure, is the applied times for the cyclic loading , and is the life expectancy of the structure only subjected to cyclic loading .

For the concrete shield tunnel lining with no initial flaws, it is a challenging subject to predict its crack initiation and make a rational life estimation. In this section, the steel reinforcements are assumed to be reliable in the long term, so the fatigue life of the shield tunnel structure is dominated by its concrete material. For its concrete with strength grade C50, the fatigue equation proposed by Guangyi Zhao is employed in this study, which is expressed as follows:where Smax denotes the maximum tension stress, , and N denotes the loading cycles.

4.3. Prediction of Fatigue Life

The basic procedure for fatigue life prediction of the shield lining is as follows:(a)Use ABAQUS to obtain the time-history of stress of the shield tunnel lining, at the speed of 80 km/h of the train.(b)Analyse the amplitude, average value, and the corresponding cyclic times for one period of the stress by the rain-flow counting method and then calculate the damage of the tunnel lining using the Miner damage theory and S-N equation of the concrete.(c)Make an estimation of the fatigue life for the shield lining based on the calculated damage above.

Numerical simulations show that the standard segment at the bottom of the shield tunnel lining ring has the highest stress level under the train vibration load. So, the standard segment at the bottom of the shield tunnel mainly determines the fatigue life of the structure. The time-history of stress of the standard segment at the bottom is illustrated in Figure 11.

The time-history curve of the stress in the standard segment at the bottom is analysed using the rain-flow counting method, and the amplitude, average value, and loading cycles of the stress are finally obtained, which is illustrated in Figure 12.

The damage to the shield tunnel segment caused by train vibration load for one time can be obtained by the above stress analysis, as in Figure 12, using the Miner damage theory and S-N equation of the concrete, which is estimated as 2.42 × 10−8, and its corresponding fatigue life lgN is about 7.62.

The train of Nanjing Metro Line 5 departs 129 times a day and consists of 6 carriages. So, the cyclic time of the train vibration load for 100 years is about N100 = 6 × 129 × 365 × 100 = 28251000, and its corresponding lgN100 = 7.451 < 7.62. Therefore, the shield tunnel lining structure of Nanjing Metro Line 5 can meet the demand of working for a hundred years under such working conditions.

5. Conclusions

In this paper, a three-dimensional numerical model of the shield tunnel lining structure is built to investigate the structure reaction and fatigue crack growth under train vibration load. Furthermore, damage of the shield segment caused by train vibration load is studied by employing the Miner damage theory and S-N fatigue law of the concrete, and thus a rational fatigue life estimation for the concrete shield tunnel lining can be finally made. The following conclusions can be drawn:(1)Crack propagation shows the overall same path under different train speeds, while the final crack lengths are a little different. The higher the train speed, the longer the final crack.(2)Train axle has a larger influence than train speed on the crack propagation. The final crack under a heavier train axle is so long that it almost cuts through the shield segment, while the final crack under a lighter train axle is so short and it only develops to a local crack.(3)The damage of the shield segment caused by train vibration load can be determined by the stress analysis using the rain-flow counting method, and the results show that the shield tunnel lining structure of Nanjing Metro Line 5 can meet the demand of working for a hundred years under the current working conditions.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was financially supported by the National Natural Science Foundation of China (grant no. 51808114) and Natural Science Foundation of Jiangsu Province (grant no. BK20170670).