Abstract

In order to evaluate the seismic vulnerability of a railway bridge, a nonlinear finite element model of typical three-span continuous beam bridge on the Sichuan-Tibet railway in China was built. It further aimed at performing a probabilistic seismic demand analysis based on the seismic performance of the above-mentioned bridge. Firstly, the uncertainties of bridge parameters were analyzed while a set of finite element model samples were formulated with Latin hypercube sampling method. Secondly, under Wenchuan earthquake ground motions, an incremental dynamic method (IDA) analysis was performed, and the seismic peak responses of bridge components were recorded. Thirdly, the probabilistic seismic demand model for the bridge principal components under the prerequisite of two different kinds of bearing, with and without seismic isolation, was generated. Finally, comparison was drawn to further ascertain the effect of two different kinds of bearings on the fragility components. Based on the reliability theory, results were presented concerning the seismic fragility curves.

1. Introduction

In recent years, with the rapid development of China’s western railway network, a large number of railway lines have emerged in the western regions. With a full length of 1,629 km, Sichuan-Tibet railway is currently under construction. To span across topographically complex regions with high mountains and deep valleys, three-span continuous beam railway bridge structure has been widely employed. Due to topographical reasons, these bridges are often provided with long span and large height. Besides, there is also very large interpier height difference, with some reaching up to 20–30 meters. For such irregular continuous beam bridges, ever-increasing pier height and span have brought about many technological problems to seismic design. In addition, there have been frequent seismic activities in the western regions in recent years, such as 2008 Wenchuan earthquake (M8.0), 2010 Yushu earthquake (M7.1), and 2013 Lushan earthquake (M7.0). These factors would seriously challenge the seismic safety of typical three-span continuous beam bridges on the Sichuan-Tibet railway.

With the development of seismic theories for bridge structures, seismic vulnerability analysis has become an important approach to evaluate the seismic behavior of bridge structures [1]. Seismic vulnerability refers to the conditional probability that seismic damage of a structure or component exceeds respective seismic capacity. At present, three types of fragility curves are mainly analyzed, namely, analytical fragility curves, empirical fragility curves, and expert opinion fragility curves [2]. In reality, there has been lack of seismic disaster data on bridge structures in many regions, due to which scholars have made more extensive and in-depth researches on the creation of analytical seismic fragility curves. In recent years, scholars have carried out a series of studies on typical bridges with the seismic vulnerability theory. Moreover, based on the reliability theory, Hwang et al. [3] analyzed the seismic vulnerability of typical three-span continuous beam highway bridge in east-central United States while Nielson and DesRoches [4] carried out vulnerability analysis on the components and systems of typical simply supported beam bridges in central and southeastern United States. Gianmarco and Giannini [5] analyzed the seismic fragility of multispan simply supported beam bridge based on the response surface methodology taking into account the uncertainties of structural parameters. Besides, Choi [6] et al. carried out seismic vulnerability analysis on four typical types of beam bridges in south-central United States, in which the influence of different components on the bridge system vulnerability was taken into consideration. Wang [7] made vertical and horizontal seismic fragility analyses on typical 32 m simply supported beam railway bridges in China, so as to further study the effects of ordinary bearings and seismic isolated ones on vulnerability. Zhou and Wang [8] analyzed the seismic vulnerability of common double-T small-span continuous bridges in China for evaluating the seismic performance of such bridges, while a seismic vulnerability analysis was also conducted by Feng [9] on piers and bearings of a continuous beam highway bridge. In addition, through the use of numerical simulation, Zhang [10] analyzed the seismic vulnerability of three-span continuous beam highway bridge. However, the above studies mainly focused on highway or railway simply supported beam bridges or small-span continuous beam highway bridges with low piers, while little attention has been paid to the seismic vulnerability of common high-pier, large-span continuous beam railway bridges. When it comes to the effects of different bearing types on the seismic vulnerability of continuous beam bridges, relevant studies have been seldom presented. Therefore, there is an urgent need for the research on the seismic vulnerability of typical three-span continuous beam railway bridges.

This paper studies a three-span continuous beam bridge (40 m + 64 m + 40 m) on the Sichuan-Tibet railway that is under construction. Firstly, a 3D finite element model of this typical bridge was built, followed by the creation of ground motion-bridge analysis sample database considering the uncertainties of bridge model parameters and ground motions. And then a dynamic nonlinear time-history analysis was performed on the sample database, based on which a probabilistic seismic demand model for the bridge principal components was also constructed. Finally, the seismic fragility curves for dangerous components of the bridge were generated, and the seismic behavior was evaluated, along with the comparison made about the effects of two different bearing types on bridge structure vulnerability.

2. Seismic Vulnerability Analysis Theories

2.1. Probabilistic Seismic Demand Model

The main purpose of probabilistic seismic demand models was to establish the relationship between ground motion intensity and structural seismic demand for a certain type of structure. Based on the seismic demand followed by lognormal distribution, Cornell et al. [11] recommended that the logarithmic medians of demand and ground motion intensity shall obey a linear relation of bridge structure and the ground motion intensity as follows:where (0, 1) is standard normal random variable; refers to the conditional logarithmic standard deviation of bridge structure seismic demand; and parameters A and B could be solved by the least square method. Linear regression was widely used in vulnerability analysis [12, 13]. However, the premise for the adoption of the method can be described as follows: ground motion intensity IM and median seismic displacement demand follow a linear relation in logarithmic space. Under the action of the th ground motion, the residual between actual seismic demand of piers and predicted value by demand model (1) wasThen, the parameter in (1) was

2.2. Seismic Vulnerability of Bridge Components

Seismic vulnerability refers to the conditional probability that a structural component reaches or exceeds a particular damage state under a given ground motion. In terms of the conditional probability that seismic demand of bridge component is greater than the seismic capacity, it was shown in the following equation:where refers to the conditional probability to exceed a particular damage state. When the seismic demand D and seismic capacity C of bridge component obeyed lognormal distribution, the probability for the bridge component to reach or exceed a particular damage state shall also obey the lognormal distribution. Thus, based on the lognormal cumulative probability density function, bridge component vulnerability function (Kennedy and Ravindra [14]) could be obtained as follows:where indicates the median seismic capacity of bridge component; and is the logarithmic standard deviation of seismic capacity of bridge component. According to literature [1], , where COV means the coefficient of variation. represents the standard normal distribution function. And other parameters have been defined in (1), (2), and (3). Demand model (1) of bridge structure can be substituted into (5), thus gettingwhere and mean the median and logarithmic standard deviation of the bridge component vulnerability function, respectively. , .

Based on the above vulnerability analysis theories, seismic fragility analysis on three-span continuous beam railway bridges could thus be implemented.

2.3. Flow for Vulnerability Analysis

To perform seismic vulnerability analysis on typical three-span continuous beam railway bridge, a flowchart for seismic vulnerability evaluation was presented based on the above analysis theories (see Figure 1).

3. Case Study of a Bridge

3.1. Project Overview

With the Brahmaputra River Bridge under construction in Lang County on the Sichuan-Tibet railway taken as an example, seismic vulnerability of typical continuous beam railway bridge was discussed based on the analytical flow in Figure 1. The superstructure of bridge was three-span variable cross-section continuous box beams with span arrangement of 40 m + 64 m + 40 m and material of C55 concrete. Piers were round-ended hollow piers made of C35 concrete. Heights of piers 1~4# were 35, 37, 40, and 7 m, respectively (Figure 2), of which pier 2# was braking pier. Moreover, those piers had a longitudinal reinforcement ratio of 1.36% and a stirrup reinforcement ratio of 0.85%, while both the longitudinal bars and stirrups used grade HRB335 steel. This paper studied the vulnerability of bridge structure under two different bearing types, which are detailed in Table 1. For spherical steel bearings, the design horizontal force of longitudinal expansion bearings accounted for 5% of the design vertical bearing capacity, while the design horizontal force of fixed bearings occupied a percentage of 30% of the design vertical bearing capacity. In terms of the hyperboloid spherical seismic isolated bearings (hereinafter referred to as the seismic isolated bearings), the design horizontal forces of both longitudinal expansion and fixed bearings accounted for 20% of the design vertical bearing capacity.

3.2. Finite Element Model

In recent years, numerical analysis approaches have played a crucial role in the seismic analysis of structures. Among the analysis approaches mentioned above, the Open System for Earthquake Engineering Simulation (OpenSEES) developed by the US Pacific Earthquake Engineering Research Center has gradually gained widespread application. In this paper, OpenSEES was also used to perform nonlinear time-history analysis about the three-span continuous beam bridge.

An existing study on earthquake damage data has shown that, for most bridge structures, main beams were in substantially linear elastic state under seismic action [15]. In the present study, the main beams of finite element model were simulated with elastic beam elements, with self-weight and secondary dead loads taken into consideration. As piers were prone to damage and failure during earthquakes, force beam column elements based on flexibility method were thus employed for the simulation [16]. Moreover, elemental nonlinearity in the model was mainly reflected in the nonlinear stress-strain relation of fiber materials. Unconstrained and constrained concrete materials were defined with Kent-Scott-Park model [17] (Concrete02 material) in OpenSEES, while reinforcement materials were defined with Giuffre-Menegotto-Pinto model [18] (Steel02 material). For simulation of bearings, bilinear ideal elastoplastic spring elements were used based on the method in literature [19]. This paper intended to investigate the effects of spherical steel bearings and seismic isolated bearings on the seismic vulnerability of components, for which the definition of preyield and postyield stiffness of two types of bearings was required according to their specific design parameters [20] (Table 1). Specific parameter defining method is detailed in the literature [1]. Pile foundation was simulated with equivalent boundary elements. That is, pile group effect was equivalently simulated using six springs at the bottom of bearing platform. In terms of the stiffness of the six springs, it consists of vertical stiffness, axial and transverse direction lateral stiffness, torsional stiffness around vertical axis, and two types of horizontal direction flexural stiffness. Equivalent stiffness of each pile foundation could be solved according to Appendix D, “computation theory of elastic resistance of pier and abutment foundations considering soil,” under Specifications for Design of Railway Bridge and Culvert Foundations (TB 10002.5-2005), in [21].

3.3. Uncertainties of Model Parameters
3.3.1. Selection of Ground Motions

Given that the bridge studied herein was a typical three-span continuous beam bridge on the Sichuan-Tibet railway, actually measured data of ground motions in the region should be selected as far as possible. In this way, the seismic vulnerability analysis results of bridge were closer to the engineering practice. Thus, the actual ground motion data of Wenchuan earthquake was selected to perform the seismic vulnerability analysis. When Wenchuan earthquake occurred in 2008, over 400 stations in the China Digital Strong Motion Network received the data about ground motion, including 46 observation stations within Sichuan with recorded PGA of ground motions greater than 0.1 g. After the comprehensive consideration of the intensity of ground motions, integrity of seismic records, and the spectral characteristics of near-field ground motions, 90 ground motion records from 30 observation stations were selected as the original inputs for the nonlinear dynamic analysis of the bridge studied. Table 2 lists the names of the 30 stations and peak ground acceleration, while Figure 3 illustrates the acceleration response spectra of various ground motions and their mean spectrum.

3.3.2. Uncertainties of Bridge Structure Parameters

For railway bridge structures, due to the effects of various factors during construction process (such as raw materials, construction methods, environments, and construction workers), the material properties, component behavior, and load effect of actual bridge structure may face certain uncertainties, which thus will lead to uncertain seismic responses and capacity of bridge structures. Thus, uncertainties of these factors need to be considered in the seismic vulnerability analysis herein. According to the structural characteristics of a typical three-span continuous beam bridge on the Sichuan-Tibet railway, and based on Pan [22] and Tavares [23] et al.’s findings and construction features in China, uncertainty parameters in the bridge model and their distribution were identified as shown in Table 3.

3.4. Determination of Bridge Structure Damage Indices

Before the generation of seismic fragility curves for the typical three-span continuous beam railway bridge, damage indices of each bridge component need to be determined. For piers, after comprehensively considering the characteristics of three-span continuous beam bridge and consulting literature findings [3, 6, 24], curvature ductility ratio was finalized as the pier damage index. Specific method for damage index calculation is detailed in the literature [22]. The curvature ductility ratio, being less than 1.00, implies elastic state and no damage for the piers, while other larger values mean more serious damage states. And the curvature ductility ratio, being larger than about 4.6, implies complete destruction for the common piers. As for bearings, since there is yet no descriptive research on quantifying the damage degree of spherical steel railway bearings or seismic isolated railway bearings, we thus identified the relations between damage degree of two bearing types and bearing displacements during study after consulting relevant literatures [25, 26]. And those bearing displacements depend on the specific bearing devices with different displacement capacities. The displacement, being less than the design value, implies no damage for the bearings, while other larger values mean more serious damage states. The displacement of 5 times the design value is assumed as the complete destruction state for the bearings. Tables 4 and 5 show the median damage indices corresponding to various damage states of three-span continuous beam railway bridge piers and bearings. Meanwhile, seismic capacities of components were assumed to follow lognormal distribution, while the definition of uncertainties for each damage state was also needed. Nielson [1] used the coefficient of variation (COV) to describe such uncertainties, assuming that COV value was 0.246 for the former two damage states and was 0.472 for the latter two damage states. In this paper, the above assumption was employed to obtain the logarithmic standard deviations of bridge component ultimate states, as shown in Tables 4 and 5.

3.5. Probabilistic Seismic Demand Analysis

To generate fragility curves for various components of the typical three-span continuous beam bridge on the Sichuan-Tibet railway, study shall be conducted on the correlation between seismic demand of each component and ground motion intensity index, so as to further clarify their variation relations. Hence, probabilistic seismic demand analysis on each component was required, with the analytical procedure shown as follows.

Firstly, based on the 90 actually measured ground motion records of Wenchuan earthquake selected in Section 3.3.1, and with the employment of IDA analysis [27], 90 ground motions were amplitude modulated to 10 different PGA intensities from 0.1 to 1.0 g at a step size of 0.1 g to form sample database consisting of a total of 900 ground motions.

To reduce sampling error, Latin hypercube sampling [28] was used to sample uncertainty parameters in the model described in Section 3.3.2, so as to generate 900 typical bridge analysis samples.

Ground motion samples and bridge samples were paired one by one to generate ground motion-bridge analysis sample database. Then, OpenSEES was used to perform time-history analysis on the sample database and output the maximum dynamic response of dangerous bridge parts, including the curvature ductility ratio of pier bottom section and the maximum relative displacement of bearing .

Probabilistic seismic demand model of each component was calculated by regression analysis on the maximum responses of components and ground motion intensity index.

Through probabilistic seismic demand analysis on the typical three-span continuous beam railway bridge according to the above procedure, demand model could be obtained for each component. Due to space limitations, only the calculation results of bridge components with spherical steel bearings could be given herein (Table 6).

3.6. Seismic Vulnerability Analysis

Based on the above probabilistic seismic demand analysis results of bridge components, as well as median damage indices of various components, fragility curves were constructed for each component using formula (6). Figures 4 and 5 illustrate the seismic vulnerability curves of piers and bearings corresponding to four damage states under two bearing types at PGA = 0–1.0 g. As can be seen, for the two different bearing types, the vulnerability curves for different components, however, shared the similar shapes. However, damage probabilities corresponding to different damage states had a slight difference. Vulnerability function parameters of components corresponding to slight and moderate damage states are shown in Table 7.

Analysis of Figure 4 showed that, for bridge piers, when PGA was 0.3 g (seismic intensity of 8 degrees), destruction probability of 2# pier (braking pier) was around 30% under slight damage state. For other non-braking piers, destruction probabilities were all lower than 10% in slight damage state. In comparison, the damage probabilities were almost zero in case of moderate, severe, and complete damage states. It can be indicated that the continuous beam bridge piers achieved a good seismic performance, and that braking pier was most susceptible to damage among all piers, with a destruction probability far higher than other piers under the same ground motion intensity. Thus, seismic isolation and reduction measured should be employed to reduce the seismic damage probability of 2# pier (braking pier).

For bridge bearings, when PGA was 0.3 g (seismic intensity of 8 degrees), probability of slight damage was lower than 40%; probability of moderate damage was also between 10 and 20%, while damage probability was almost zero in severe and complete destruction states. It can be indicated that the spherical steel bearings of the continuous beam bridge were susceptible to slight and moderate damage. In addition, through the comparison of vulnerability between piers and corresponding bearings, it was shown that the bearings were more prone to damage.

Figure 5 showed that, for bridge piers, when PGA was 0.3 g (seismic intensity of 8 degrees), probability of slight damage was all lower than 21% and probability of moderate damage was lower than 11%, while damage probability was almost zero in severe and complete damage states, which indicated that the continuous beam bridge piers had a good seismic performance. Moreover, through the comparison of seismic vulnerability between piers, it was shown that braking pier was most susceptible to damage among all piers. However, compared with the spherical steel bearings, seismic vulnerability difference between piers became smaller after the use of seismic isolated bearings. It was mainly because, after the use of the seismic isolated bearings, horizontal bearing capacities of bearings were all 20% of vertical bearing capacities, which varied greatly from the horizontal bearing capacity of spherical steel bearings [29]. As a result, the piers were distributed more uniformly, along with fewer differences concerning the damage probabilities of piers.

As for bridge bearings, there was small difference in damage probability between them. When PGA was 0.3 g (seismic intensity of 8 degrees), probability of slight damage was lower than 32%; while probability of moderate damage was lower than 15%, thus indicating that the seismic isolated bearings of the continuous beam bridge were susceptible to slight damage. In addition, it can be found through the comparison of vulnerability between piers and corresponding bearings that the bearings were more prone to damage. This finding was similar to the calculation results under spherical steel bearing conditions.

To further understand the seismic vulnerability of each component, the median exceedance probability ground motion intensity index was used to describe the component vulnerability. That is, in case of certain state of damage, components with smaller median PGA were more prone to seismic damage and failure.

Figures 6 and 7 show the histograms for comparison of median exceedance probability PGA between various components in slight and moderate damage states with two bearing types.

Figure 6 shows that under spherical steel bearing conditions for slight damage state, median PGA values were the smallest for 1# and 3# bearings, namely, 0.36 and 0.39 g, respectively. It indicated that 1# and 3# bearings were most susceptible to slight damage, followed by 2# braking pier prone to slight damage. For moderate damage state, 2# pier (braking pier) had the smallest median PGA, with a value of 0.51 g, thus indicating that 2# pier was most susceptible to moderate damage, followed by 1# pier bearing, 3# pier bearing, 3# pier, 1# pier, and 4# pier. Besides, vulnerability ranking results somewhat differed between slight and moderate damage states.

Figure 7 shows that, under seismic isolated bearing conditions for slight damage state, median PGA ranged between 0.41 and 0.83 g for various components. With median PGA values being 0.41 and 0.45 g, respectively, 3# and 1# pier bearings were most susceptible to slight damage, followed by 2# pier, 3# pier, 1# pier, and 4# pier bearing. For moderate damage state, 2# and 3# piers were most susceptible to damage and failure, followed by 3# pier bearing, 1# pier bearing, 1# pier, 4# pier bearing, and 4# pier, which somewhat differed from the results of slight damage state.

Based on the comparison between Figures 6 and 7, it was shown that, by the action of different bearing types, there was an obvious difference with the median PGA of components in various damage states. Therefore, it was necessary to conduct an in-depth study on the seismic vulnerability of the same components under two bearing types; besides, it was also necessary to investigate the differences between the two, so as to provide the basis for seismic design of typical three-span continuous beam railway bridges.

3.7. Comparative Analysis under Different Bearing Types

For studying the effects of two types of bearings on the seismic behavior of different components, variation trends relating to damage probability deviations of various vulnerable components with PGA under two bearing type conditions were presented (Figure 8).

Figures 8(a)8(c) show that, for various bearings, damage probability would experience a decrease under seismic isolated bearing condition. In slight, moderate, and severe damage states, damage probability deviations corresponding to 1# and 3# bearings all presented the first increase at first, followed by the decrease trend with an increase of PGA. For 4# pier bearing, its damage probability deviations in four corresponding damage states all increased with an increase of PGA.

For 3# bearing, its maximum damage deviations in various damage states would not exceed 5%, which meant that the damage probabilities of the 3# pier bearing in various damage states with the presence of seismic isolated bearing were not significantly lower than that with spherical steel bearing. That is, there was no qualitative difference in the vulnerability results of the component under two different bearing types of conditions. This was mainly due to the large difference in the displacement control parameter between the two bearing types (Table 5). With possible improvement of displacement control parameter of the seismic isolated bearing (e.g., adjusting to displacement parameter 150 mm to be the same as the spherical steel bearing), there would be qualitative difference in the reduction rate of damage probability of the seismic isolated bearing.

For 1# and 4# side pier bearings, their maximum damage probability deviations in slight, moderate, and severe damage states all exceeded 10%. In terms of 1# side bearing, the maximum damage probability deviation in severe damage state would also exceed 10%. It indicated that, for side bearings (1# and 4# bearings), the seismic isolated bearing had excellent seismic reduction and isolation effects despite the different displacement control parameters between the spherical steel bearing and the seismic isolated one.

Figures 8(d)8(g) show that, for 1#, 3#, and 4# piers, their damage probability deviations were negative, while, for 2# pier (braking pier), its damage deviations in various damage states were all positive. It indicated that, after the use of seismic isolated bearing, damage probabilities of 2# pier (braking pier) for corresponding damage states all experienced a decrease by seismic action, while, for the remaining piers (1#, 3#, and 4# piers), damage probabilities somewhat experienced an increase in various damage states. This was mainly because the horizontal ultimate bearing capacity increased by 15% for 1#, 3#, and 4# movable pier bearings after the use of the seismic isolated bearing compared with the use of spherical steel bearing. When it comes to 2# pier fixed bearing, on the other hand, its horizontal ultimate bearing capacity decreased by 10% after the use of the seismic isolated bearing compared with the application of spherical steel bearing. Thus, under ground motion, damage energy borne by 1#, 3#, and 4# piers (non-braking piers) would increase, and, correspondingly, their damage probability would also increase [30, 31]. At that time, due to the lowered damage energy borne by 2# pier, the damage probability under ground motion, therefore, would decrease for 2# pier.

For 2# pier (braking pier), its seismic damage probability could be reduced by around 13% at maximum in slight to moderate damage states after the employment of the seismic isolated bearing, which occurred at a PGA of around 0.4 g. In severe and complete damage states, the seismic isolated bearing could reduce the damage probability by 7% and 2% at maximum, respectively, which occurred at 1.0 g. For 1#, 3#, and 4# piers (non-braking piers), moreover, their damage probabilities in slight to moderate damage states all increased after the use of seismic isolated bearing, with a maximum increase of around 20%. However, as shown in Figure 5, after the application of seismic isolated bearing, 2# pier (braking pier) was still the most susceptible to damage [3234]. Nevertheless, seismic vulnerabilities of the other three piers were closer to the braking pier, thus proving that, after the use of the seismic isolated bearing, four piers could more uniformly withstand the earthquake damage energy.

4. Conclusions

Based on IDA analysis and taking into account uncertainties in the ground motion and bridge structure parameters, seismic vulnerability analysis is performed on the structure of a typical three-span continuous beam bridge on the Sichuan-Tibet railway. Moreover, comparative studies have been conducted on the effects of two types of bearings (spherical steel bearing and seismic isolated bearing) on the vulnerability of dangerous bridge components. The main conclusions are drawn as follows:

Piers of the typical three-span continuous beam bridge on the Sichuan-Tibet railway have a good seismic performance. In case of the earthquake with PGA less than 0.3 g (seismic intensity of 8 degrees), slight damage probability of 2# pier (braking pier) was lower than 30%. When it comes to other non-braking piers, however, the probability was lower than 10%. With good seismic performance, severe damage or complete destruction is almost unlikely for piers.

Through the comparative study on the vulnerability of piers and their corresponding bearings, it was found that bearings are more susceptible to damage than piers, thus suggesting that bearings are seismic damage prone components in the continuous beam railway bridges.

With the use of spherical steel bearing, vulnerability of 2# pier (braking pier) is significantly higher than the other three piers. After the application of seismic isolated bearing, vulnerability gained a reduction for the 2# pier (braking pier), while there would be an increase for the other three piers. However, the vulnerability of braking pier is still the greatest. Besides, at this time, fragility curves of the other three piers are relatively close to the braking pier, which indicates that, with the use of seismic isolated bearing, four piers are able to withstand the earthquake damage energy more uniformly.

Based on the analysis of deviation plots of damage probability of bridge pier bearings, the damage probability deviations are more distinct for side pier bearings. However, 3# central pier bearing has very small damage probability deviations with the seismic isolated bearing. Therefore, the adjustment of displacement control parameter for the seismic isolated bearings is recommended, in order to further reduce the seismic vulnerability of each bearing.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research is jointly supported by the National Natural Science Foundation of China (51338004, 51778635, 51308549, 51378504, and 51478475), the Natural Science Foundations of Hunan Province (2015JJ3159), the Special Fund of Strategic Leader in Central South University (2016CSU001), the Innovation-Driven Plan in Central South University (2015CX006), the Experimental Foundations of the New Railway from Chengdu to Lanzhou (CLRQT-2015-010), and the Graduate Innovative Project in Central South University (2016zzts399). The support is greatly appreciated.