Abstract
Based on long-term inspection data about track geometric irregularity of 32 m-span simply supported girder bridge in China Railway High-Speed for nearly 10 years, it is found that there is continuous periodic track longitudinal level (TLL) with different wavelengths, thereinto, simply supported girder bridges creep and ballastless track slab upwarp are the primary reasons. Moreover, there are two types of variation trends for bridge creep of 32 m simply-supported bridges with ballastless track, that is, stable type and increasing type, and the maximum
Keywords
Introduction
It has been 10 years for the operation of China high-speed railway. Various deficiencies of infrastructure and track structure occur, which cause obvious track geometric irregularity, thereby affecting ride comfort and running safety. For simply supported girder bridge with continuous equal span, bridge creep will cause continuous periodic TLL 1 (track longitudinal level, also referred to as longitudinal irregularity or profile). Under various actions of temperature2,3 creep or train vibration load, 4 track slab deformations would also occur, such as arching or upwarp. Wang 5 studied the influence of bridge pier settlement and girder bridge creep, which are two typical types of long-term bridge deformation, on the vibration of high-speed trains. Based on the vehicle-bridge coupling dynamics theory, Yang et al., 6 Gou et al. 7 investigated the dynamic response of vehicle-track-bridge structure and running performance of vehicle when there is bridge creep. Xia et al. 8 studied the influence of PC beam deformation induced by creep and temperature effect on the dynamic responses and running safety of HSR train-bridge system. Zhao et al. 9 developed two finite element models of CRTS II ballastless track superstructure were to simulate temperature deformations of the precast concrete track slab. Based on dynamic inspection data of TLL, Tian et al. 1 analyzed the causes of periodic TLL of simply supported girder bridge, and proposed the control standard for bridge creep via the vehicle-track-bridge vertical coupling vibration analysis model. Ren et al. 10 studied influence of cement asphalt mortar debonding on the damage distribution and mechanical responses of CRTS I prefabricated slab, Li et al. 11 studied the mechanical behavior characteristics of CRTS II track slab on bridge and the mean value, standard deviation, and variation coefficient of periodic track irregularity. Zhu 12 revealed the effect of train–track interaction and environment loads on the mechanical characteristic variation of ballastless tracks in high-speed railways. To capture the mechanical characteristic variation of ballastless tracks under environment loads such as the temperature load, Zhong 13 established a three-dimensional finite element model of slab track to study the deformation and interface stress of a slab track under daily changing temperature. Wang et al. 14 evaluated the mechanical performance of the CA mortar of the slab track interlayer under extreme temperatures, by conducting the flexural strength and fatigue tests and a transient thermal simulation. Cai et al. 15 investigated the arching mechanism of slab joints under high temperature conditions considering the buckling instability of track structure, and the damage process of joint concrete and its effects on the arching instability were analyzed with the help of the field investigation.
TLL caused by bridge creep and slab upwarp deformation is continuous and periodic, which has a great influence on ride comfort and running safety. However, there are few researches on the track state influenced by bridge creep and slab upwarp deformation simultaneously. The variation trend of TLL caused by bridge creep and slab upwarp deformation for 32 m simply supported bridge obtained using 10 years’ high-speed comprehensive inspection train data. To make the track management more targeted and ensure the safety and comfort for high-speed railway operation, the threshold values are proposed based on the simulation results of 3D vehicle-track-bridge coupling dynamic model considering current operation and maintenance situation. 16
Analysis of TLL on the bridge
For 32 m-span simply supported girder bridge with ballastless track slab in high-speed railway, the bridge creep and slab upwarp deformation are two main factors affecting the TLL. Figure 1 presents the schematic of ballastless track slab on 32 m-span simply supported girder bridge and common slab deficiencies. Due to significant difference between structures, TLL with different wavelength would be caused by both deficiencies. 17 For example, the TLL in C line (represents a high-speed railway line that has been in operation for several years, A line and B line are the same) of high-speed railway with CRTS I caused by bridge creep and track slab deformation is shown by the black line in Figure 2(a). As shown in Figure 2(b) of the analysis results by power spectral density, 18 the wavelength of TLL caused by bridge creep is 32.7 m, and the values caused by track slab deformation are mainly 4.67 and 5.45 m, which is basically the same as the length of bridge and CRTS I. Thus the individual influence of bridge creep and track slab deformation on TLL can be identified by low-pass, high-pass, or band-pass filtering methods because of the different wavelengths, as shown in Figure 2(a).

Ballastless track slab and deficiencies of 32 m-span simply supported girder bridge in high-speed railway: (a) layout drawing of ballastless track slab on 32 m-span simply supported girder bridge, (b) mud pumping and rising soil of CRTS I track slab, and (c) the broken wide and narrow joints of CRTS II track slab.

The TLL caused by bridge creep of 32 m-span simply supported girder bridge and track slab upwarp deformation in C line of high-speed railway and power spectral density: (a) the TLL of in C line of high-speed railway and (b) the PSD of TLL of left and right rail.
Analysis of TLL caused by bridge creep
Thirty-two metre simply supported PC girders with box sections are widely used in China Railway High-Speed. Under the interaction of prestressed steel and dead load, each section is in an eccentric compression state. In a section, the compressive stress varies along the height linearly, with larger compressive stress in the upper flange, and less stress in the lower flange, even tensile stress. Due to the creep characteristics of concrete, under the reverse bending moment caused by prestressed steel, the creep occurs, 19 as shown in Figure 3. The main factors affecting the creep include the arrangement of prestressed tendons, elastic modulus, strength of concrete, the age at prestressing and time between prestressing and application, and the age at of non-structural loads. Moreover, the actual train load is lower than the design value, which also intensifies the amplitude of bridge arching.

Diagram for creep of simply supported girder bridge.
Slab ballastless track slab and rail are laid on 32 m-span simply supported girder bridges. Slab ballastless track is consolidated with bridge deck. The vertical stiffness of bridge is much larger than the bending stiffness of slab ballastless track and rail. The deformation caused by bridge creep can be transmitted to the rail surface through ballastless track. As a continuous elastic body, there is limited mitigation effect for rails on the beam end angle of simply supported girder bridge. 6 However, due to the poor bending stiffness of rail, the difference between the vertical deformation of rail and bridge creep deformation is small. 1 Therefore, the creep characteristics of simply supported girder bridge can be obtained through the TLL. Bridge creep will cause continuous periodic TLL, as shown in Figure 4, when the TLL is large, the standard deviation of TLL data will exceed the threshold value 11.

Influence on TLL of continuous 32 m simply supported girder bridge creep.
Ballastless track slab is laid on the simply supported bridges of high-speed railway, with the increase of service time, slab upwarp deformation might occur. Because two structural deformations cause different wavelength components, band-pass filtering of track dynamic longitudinal level was carried out, and the components with wavelength at 10–42 m were refined. The wavelength of angle irregularity caused by creep at beam end is relatively short, and band-pass filtering will change the real situation of TLL of beam joint. Meanwhile, considering the influence of different settlement of piers, the median of the filtered peak to the two valleys of the original waveform is taken as the
Where

Schematic diagram of TLL and
High-speed comprehensive inspection trains are used to regular detect track geometric irregularity in China.20,21 By nearly 10 years, massive track geometric irregularity data have been obtained. There is certain temperature effect for simply supported girder bridge, that is, the TLL is larger in summer and smaller in winter. To better ensure running safety, the maximum value envelope curve was built to capture the development trend of bridge creep. As seen in Figures 6 and 7, the
Where z(t) represents the variation of

Development of a 32 m-span simply supported girder bridge creep in A line of high speed railway.

Development of a 32 m-span simply supported bridge creep in B line of high-speed railway.
Figures 6 and 7 show the trend of
Tian et al.
1
pointed out that the theoretical calculation of the maximum bridge creep deformation is about 6.2–7.3 mm for 32 m-span simply supported girder bridges. Based on the statistics result of
Analysis of upwarp deformation characteristics of track slab on the bridge
There are various forms of slab ballastless track in China Railway High-Speed, 25 such as CRTS I,26–28 CRTS II,29,30 CRTS III, 31 and other forms. Under the action of train traffic load and external environment effect, ballastless track slab is prone to various deficiencies, such as mud-pumping, mortar seam opening, joint damage, and upper arch of track slab end. 32 According to track dynamic detection data, the slab arching has become one of the most prominent deficiencies, and the superposition of slab arching and bridge creep would intensify the dynamic response of vehicle.
CRTS II is laid on the simply supported girder bridge in B line of high-speed railway. The track geometric irregularity from October 2011 to August 2020 are shown in Figure 8(a). The TLL of three-span simply supported girder bridge near K532+500 (equal to 532.500 km, same bellow) not only includes the 32.7 m wavelength component caused by bridge creep, but also includes the 6.5 m wavelength component, which is equal to the length of CRTS II. This point can be proved by the power spectral density of TLL around K532+400–K532+600, as shown in Figure 8(b). Meanwhile, it can be seen from Figure 8(c) that the slab upper arching at K532+496 occurs in July 2017, and the

(a) The superposition curve of TLL of B line of high-speed railway. (b) Power spectral density of TLL in K532+400–K532+600. (c) The variation of the
Model of vehicle-track-substructure interaction
The accurate three-dimensional finite element model of vehicle-track-slab-bridge is established, as shown in Figure 9. And the dynamic wheel-rail coupling relationship, the reasonable equations and solving methods are established. The influence of slab and bridge deficiencies on TLL is taken as the input excitation of the whole model, and the dynamic simulation calculation is carried out to provide the premise for evaluating the driving performance and structural service performance.

Vehicle-track-bridge dynamic interaction model.
Vehicle Model
There are multiple vehicles for a train, each of which contains one body, two bogies, and four wheelsets, and each of the seven parts is considered as a rigid body. The analysis of two suspension train is carried out, and 35 independent degrees of freedom are considered, including lateral motion, vertical motion, roll motion, pitch motion, and yaw motion, respectively, as shown in Table 1. The vehicle dynamic analysis model is presented in Figure 10. The physical meaning of each variable in Figure 10 is shown in reference. 33 The CRH380A is used as the dynamic simulation model. To improve calculation efficiency, one motor car and one trail car are selected, and they are different in quality, size and the suspension parameter.
Degrees of freedom of passenger vehicle dynamics model.

Vehicle dynamic analysis model: (a) side view of vehicle model and (b) front view of vehicle model.
Dynamic model of Substructure
Compared with other slab types, the length of CRTS I track slab is the shortest, and the wavelength of TLL caused by slab deformation is also the shortest, which is more unfavorable to running safety. Therefore, this paper establishes a three-dimensional finite element model of 10-span 32 m simply supported girder bridge-CRTS I-seamless rail based on a commonly used finite element software, and the model is called substructure including rail, fastening system, track slab, base slab, and bridge structure components, as shown in Figure 11(a). The detailed model is shown in Figure 11(b) and (c). Some material characteristic parameters of the bridge can refer to this paper. 34 Furthermore, the length of every bridge span is 32.6 m, the length between the cross sections of the two bearings is 32 m, and the gap width between the center line of every bridge span is 0.01 m, finally the overall length of rail is 900 m considering the influence of train entering and leaving the bridge, one about 300 m track is added on both sides of bridges, then the train will go from one end of the track to the other. And the height of every pier is 20 m (a average value of high-speed railway pier, and pier height has little influence on TLL and driving performance, so the average value is used for calculation.). Moreover, one end of the track fastener is connected to the rail and the other is consolidated on the ground. In the model, the beam element is used to simulate the rail of CHN60, and spring-damping element simulates the constraint of fastener and rubber on rail, the constraint between sleeper and track slab, and the connection between track slab and base slab instead of the CA mortar mixture, with the fastener spacing being 0.629 m, in addition, the rigid connection is used to restrict motion between the base slab and bridge. According to structure and stress characteristics, the track slab, base slab, and simply supported girder bridge are modeled by beam element with 12 freedoms for two nodes. 35 The section parameters of simply supported girder bridge are introduced in literature. 8 The main-slave relation simulates the bridge bearing without considering the friction of the active bearing, and the piers bottom are fixed on the ground, leaving out the friction of the active bearing. Consequently, the whole rail-track-bridge model has 9042 nodes and 11,375 elements. In addition, the structural parameters of track slab on bridge are listed in Table 2.

Three-dimensional finite element model of rail-track-bridge: (a) Overall model drawing, (b) detail drawing of the bridge and (c) detail drawing of the track.
The parameters of CRTS I track slab.
The values in brackets express the track slab size of the bridge end.
Rayleigh damping is adopted for the damping matrix of track system and substructure. Meanwhile, the influence of fasteners on the stiffness matrix and damping matrix is considered according to the principle “seat by number.” 36 The Rayleigh damping calculation can be expressed as:
where
where
The vertical static stiffness of rail pad is 25 MN/m, and the vertical damping is 75 kNs/m. The lateral static stiffness of rail pad is 20 MN/m, and the lateral damping is 60 kNs/m. 33
Wheel-rail contact principle
The Hertz nonlinear elastic contact theory is adopted in the normal direction of wheel-rail, the vertical force between wheel and rail can be expressed as:
Where G is the wheel rail contact constant, the units is m/N2/3.
First, the longitudinal creep force
where
Then, the Shen-Hedrick-Elkins theory is used to make nonlinear corrections for the longitudinal and lateral creep slip forces between wheel and rail, and the calculation details and parameters description are presented in the literature. 33
Selection of TLL
Creep of prestressed concrete simply supported girder bridge and slab upwarp deformation will cause additional rail deformation (also referred to as additional TLL). To obtain the mapping relationship between two deformations and TLL, based on the finite element model of vehicle-track-bridge shown in Section “Dynamic model of substructure,” the specified forced displacement is arranged, accordingly, the rail deformation curve is the additional TLL. The bridge creep deformations of 5, 6, 7, 8, 9, and 10 mm are applied to the mid-span of each bridge, and the rail deformation curve caused by the bridge creep is obtained, as shown in Figure 12. Similarly, the slab upwarp deformations of 1, 2, 3, 4, and 5 mm are applied to the middle of each track slab, and the rail deformation curve caused by the slab upwarp deformation is obtained, as seen in Figure 13. Two kinds of rail deformations are input as the additional track irregularity for the vehicle-track-substructure dynamic interaction analysis.

Rail vertical deformation caused by creep of simply supported girder bridge.

Rail vertical deformation caused by slab upwarp deformation.
According to the ballastless track spectrum of China Railway High-Speed, 24 the calculation formula is shown as:
where
The frequency-time domain conversion method is used to generate random track irregularity, and the generated random track irregularities of longitudinal level and alignment are shown in Figure 14. Then the additional TLL caused by the bridge creep and track slab upwarp are superimposed with the random irregularity of longitudinal level as wheel-rail excitation. Results of bridge creep of 7 mm superimposed random irregularity of longitudinal level are shown in Figure 15. Results of track slab upwarp of 4 mm superimposed random irregularity of longitudinal level are shown in Figure 16. Results of bridge creep of 7 mm superimposed track slab upwarp of 4 mm and random irregularity of longitudinal level are shown in Figure 17.

Random track irregularity.

Superposition of bridge creep of 7 mm and random track irregularity.

Superposition of slab upwarp deformation of 4 mm and random track irregularity.

Superposition of bridge creep of 7 mm, slab upwarp deformation of 4 mm, and random track irregularity.
Coupled equation and numerical solution method for vehicle-track-substructure
Vehicle, track, and bridge constitute the dynamic simulation analysis system, which is coupled by wheel-rail force. The system vibration should obey the force equilibrium and deformation compatibility at any time. The motion equation of the system can be expressed as:
Where
where
Where
The Newmark-
Determination threshold of bridge creep and slab upwarp deformation based on running performance and structural service performance
Based on vehicle-track-substructure dynamic model, track longitudinal level caused by bridge creep, slab upwarp deformation and superposition of both are determined to be various calculation cases, with the train speeds at 200, 230, 250, 280, 300, 330, and 350 km/h for each case. The car body vertical vibration acceleration, Sperling index, rate of wheel load reduction (ΔP/P), derailment coefficient, and wheel-rail vertical force are used to evaluate the running performance, and the detailed reference standard mentioned above is presented in the literature. 37 Finally, the vertical displacement at midspan and vertical accelerations at midspan and beam end for tracks, slabs and bridges are employed to investigate the structure service performance. Based on running safety and ride comfort, the deformation thresholds are proposed.
Case study of bridge creep
According to current situation of bridge creep for 32 m-span simply supported girder bridge for high-speed railway and track geometric irregularity presented in Section “Selection of TLL,” the

Vertical acceleration time curve of car body and its FFT: (a) time curve and (b) FFT.

Rate of wheel load reduction time curve and its FFT when the bridge creep is 7 mm: (a) time curve and (b) FFT.
From Figure 20(a) to (e), it can be seen that the

Vehicle dynamic response: (a) Vertical acceleration of car body, (b) Sperling/vertical, (c) rate of wheel load reduction, (d) derailment coefficient, and (e) wheel-rail vertical force.
Maximum dynamic response results of structure.
Analysis of track slab deformation
According to current situation for slab ballastless track deficiencies for high-speed railway and track geometric irregularity presented in Section “Selection of TLL,” the

Vertical acceleration time curve of car body and its FFT: (a) time curve and (b) FFT.

Rate of wheel load reduction time curve and its FFT: (a) time curve and (b) FFT.
It can be obtained from Figure 23, When the amplitude of track slab deformation reaches 5 mm and train speed is 350 km/h, the ΔP/P is 0.82, beyond the safety threshold of 0.8. With the same track slab deformation, the derailment coefficient is 0.81 and 0.82 at train speed of 330 and 350 km/h, respectively, exceeding the safety threshold of 0.8. The wheel-rail vertical force increases rapidly with the increase of

Vehicle dynamic response: (a) vertical acceleration of car body, (b) Sperling/vertical, (c) rate of wheel load reduction, (d) derailment coefficient, and (e) wheel-rail vertical force.
Maximum dynamic response results of structure.
Analysis of superposition of bridge creep and slab upwarp deformation
Based on the above bridge creep and slab upwarp deformation, the

Vertical acceleration time curve of car body and its FFT: (a) time curve and (b) FFT.

Rate of wheel load reduction time curve and its FFT: (a) time curve and (b) FFT.
From Sections “Case study of bridge creep” to “Analysis of track slab deformation” and Figure 26, it can be seen that the creep and slab deformation have great influence on ride comfort. The bridge creep has little influence running safety, however, the slab deformation has great influence on running safety. When the

Vehicle dynamic response: (a) vertical acceleration of car body, (b) Sperling index/vertical, (c) rate of wheel load reduction, (d) derailment coefficient, and (e) wheel-rail vertical force.
The vertical vibration acceleration of car body and Sperling index increase with the increase of both deformations and train speed. When the
Maximum dynamic response results of structural.
Conclusions
(1) To make the track management more targeted and ensure the safety and comfort for high-speed railway operation, 10 years’ high-speed comprehensive inspection train data are analyzed, and the variation trends and level for bridge creep and ballastless track slab upwarp are given, which can give more clear understanding for bridge creep of 32 m simply-supported bridges with ballastless track and the ballastless track slab, and can guide the design, construction, and management. The main conclusions are as follows. Simply supported girder bridge creep and ballastless track slab upwarp are two important reasons for periodic track longitudinal level in high-speed railway. And these deformation mainly affect the vertical vibration response of vehicle.
(2) There are two types of variation trends for bridge creep of 32 m simply-supported bridges with ballastless track China Railway High-Speed, that is, stable type and increasing type. The maximum
(3) The maximum
(4) The car bodies are more sensitive to bridge creep and produce mainly low-frequency vibration, but the rate of wheel load reduction is more sensitive to the track deformation and shows the high-frequency vibration.
(5) Bridge creep affects ride comfort significantly. For 32 m-span simply supported girder bridges of 350 km/h high-speed railway without slab upwarp deformation, it is proposed that the
(6) Due to the short wavelength of TLL, the slab upwarp deformation has influence on both the ride comfort and running safety heavily. For 350 km/h high-speed railway with small bridge creep (less than 5 mm), it is proposed that the
(7) When the bridge creep and slab upwarp deformation exist simultaneously, the
Footnotes
Handling Editor: Chenhui Liang
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work is financial supported by the National Natural Science Foundation (No. 52278465), China State Railway Group Corporation Limited (No. SY2021G002), and China Academy of Railway Sciences Corporation Limited science foundation (No. 2021YJ250). These supports are gratefully acknowledged.
