Safety Evaluation of the Wear Life of High-Speed Railway Bridge Bearings by Monitoring Train-Induced Dynamic Displacements

Bridge bearings experience numerous small-amplitude displacements under environmental loads. ,e continuous cyclic accumulations of these small-amplitude displacements will result in severe wear on the poly-tetra-fluoro-ethylene (PTFE) plates in the bridge bearings, which seriously endangers the service life of bearings. Traditional method directly uses the linear wear rate of cumulative displacements in a short period to evaluate the wearing life, but the linear wear rate only in a short period such as several days may not represent the characteristics in the whole bridge service life. Hence, this research takes the spherical steel bearings of the Nanjing Dashengguan Yangtze River Bridge as a study object. ,e cumulative dynamic displacement (CDD) under the action of a single train and the cumulative bearing travel (CBT) under the continual actions of many trains are studied using the monitored longitudinal displacement data from spherical steel bearings. Furthermore, the probability statistics and the Monte Carlo sampling simulation for CDD are studied, and the safety evaluation method for bearing wear life in the real environment is proposed using a reliability index regarding the failure probability of monitored CBTover the wear limit during service lifetime. In addition, safety evaluation on the bearing wear life was performed to assess the condition of spherical steel bearings in the real service environment. ,e results can provide an important reference for analysis on the bearing wear life of long-span railway bridge structures.


Introduction
As the service time of railway bridges increases, some bridge components may suffer from the aging and damage problems due to the influence of external environment loads.Among these bridge components, the bridge bearing is an important load-bearing component in large-span railway bridge structures, which can suffer from a series of diseases, such as the disengagement of bearings, the uneven compression, the twisting and fracturing of the bearing plate, and the weld cracking, which directly threaten the safety of railway bridges [1][2][3][4][5].For example, Wei et al. revealed that the fixed bearings of continuous high-speed railway bridges were vulnerable under earthquakes because they transferred most of the seismic forces between superstructure and piers, and the sliding friction action of fixed bearings can protect other components from severe damage under earthquakes [2]; the research results of Oh et al. showed that the rubber material characteristics and spherical shape of spherical elastomeric bearing were important parameters in reducing vibrations of railway bridge [3]; Huang et al. emphasized that bearing is the important part that connects the superstructure and substructure of a railway bridge, in which the occurrence of defects affects running safety [4].erefore, it is necessary to further study the bearing damage of high-speed railway bridges subjected to environmental loads to ensure safe operation.
Current research on bearing damage in high-speed railway bridges mainly concentrated on the following three aspects: (1) the causes of bearing damage and prevention methods [4,6,7].For example, Yan conducted a detailed study on the cause analysis and remediation measures for the roller bearings of the Sunkou Yellow River Bridge [7], and Huang et al. researched a rapid treatment technique for the damaged rocking axle bearing of a railway simply supported beam bridge [4].(2) e damage recognition methods of bridge bearings [1,8,9].For instance, Zhan et al. utilized the dynamic response of a bridge under an impact load to diagnose the bearing disengagement of existing rubber bearings [1]; Qiao conducted an in-depth study on a damage identification technique for bridge bearings based on the change rate of natural frequency [9].(3) e influence of bearing damage on the mechanical properties of bridge structure [3,[10][11][12].For instance, one study revealed that the increased vibration level of a railway bridge could cause bridge bearing damage, and the rubber material characteristics and spherical shape were important parameters in reducing bridge vibration [3]; another study showed that the load-carrying fillet weld around the gusset plate of diagonal bracing at the bridge bearing was predicted to be the most fatigued critical detail among the various structural details, and a fatigue assessment was carried out for a composite railway bridge on which high-speed trains traveled [12].
e latest studies revealed that bridge bearings experienced numerous small-amplitude dynamic displacements under dynamic loads [13,14].Continuous accumulations of small-amplitude dynamic displacements can cause severe wear to poly-tetra-fluoro-ethylene (PTFE) plates of bridge bearings, which seriously decreases the service lifetime of bearings [13,14].For instance, according to the field results, the small but rapid girder movements contributed the most to the large cumulative movements, resulting in accelerated fatigue damage and bearing wear [13,14].However, current research mainly focuses on highway bridge bearings rather than railway bridge bearings.Large-span high-speed railway bridges have the special characteristics of being designed for high-speed trains.For example, the Dashengguan Yangtze River Bridge is one large-span steel truss arch bridge with a 336-meter-long main span and six tracks, and the designed high-speed train passed speed can reach 300 km/h [15].
ese conditions can easily induce a large number of smallamplitude dynamic displacements, which can result in serious wear damage.However, current research has not clearly revealed the cumulative wear characteristics and the safety evaluation method for railway bridge bearings in the service lifetime.
erefore, this research takes the spherical steel bearings of the Nanjing Dashengguan Yangtze River Bridge as a study object as shown in Figure 1 to investigate the cumulative dynamic displacement (CDD) under the action of a single train and the cumulative bearing travel (CBT) under the continual actions of many trains using monitored longitudinal displacement data of spherical steel bearings.Furthermore, the probability statistics and the Monte Carlo sampling simulation for CDD are studied, and the safety evaluation method for the bearing wear life is proposed using a reliability index regarding the failure probability of monitored CBT over the wear limit in the service lifetime.In addition, the safety evaluation for the bearing wear life was performed to assess the condition of spherical steel bearings in real service environment.Traditional method directly uses the linear wear rate of cumulative displacements in a short period to evaluate the wearing life, and such linear wear rate only in a short period such as several days may not represent the characteristics in the whole bridge service life, so the novelty of this paper is the safety evaluation method for the bearing wear life using a reliability index integrated with probability statistics analysis and Monte Carlo sampling simulation in the whole bridge life.

Description of Bridge Health Monitoring.
e Nanjing Dashengguan Yangtze River Bridge is a large-span highspeed railway steel truss arch bridge.It is designed for river crossing of the Beijing-Shanghai high-speed railway and the Shanghai-Wuhan-Chengdu high-speed railway, as shown in Figure 2. e upper structure of this bridge is a six-span continuous steel truss arch girder, which is composed of steel truss arch and steel bridge deck [16].In detail, the steel truss arch is composed of chord members, diagonal web members, vertical web members, and horizontal and vertical bracings; the steel deck is composed of top plate and transverse stiffening girder.In addition, a total of seven groups of spherical steel bearings (i.e., bearings 1∼7) are used to support the steel truss arch girder, where the bearing 4 is totally fixed in all the three translational directions while the other six bearings allow longitudinal thermal expansion in the steel truss arch girder.erefore, 12 LVDT displacement sensors (namely linear variable differential transformers) are installed at the upstream and downstream sides of the six groups of bearings (i.e., B 1,u ∼ B 3,u , B 5,u ∼ B 7,u , B 1,d ∼ B 3,d , and B 5,d ∼ B 7,d ) for long-term monitoring of the longitudinal displacement data, where "u" indicates the upstream side, and "d" indicates the downstream side.In detail, the voltages in the displacement sensors which are termed as analog signals are amplified, filtered, and then transformed to the digital signals in the data acquisition station at bridge site, which can be remotely monitored, displayed, and downloaded through Internet in the health monitoring center far from  the bridge site.e sampling frequency of the 12 displacement sensors is 1 Hz.Because the Nanjing Dashengguan Yangtze River Bridge is one symmetrical structure, only the bearings on the left half (i.e., B 1,u ∼ B 3,u and B 1,d ∼ B 3,d ) are selected for analysis.

Cumulative Dynamic Displacement Caused by a Passed
Train.
e monitored longitudinal displacement data at the ith spherical steel bearing B i,u are denoted by D i,u ; the monitored longitudinal displacement data at the ith spherical steel bearing B i,d are denoted by D i,d ; i � 1, 2, 3.Because the time histories of longitudinal displacements have similar change trends [16], only D 1,u is selected for analysis.e time history of D 1,u from March to October in 2017 is shown in Figure 3(a), showing typical seasonal characteristics over eight months.Furthermore, the time history of D 1,u on May 20th, 2017 is plotted in Figure 3(b), showing typical sinusoidal characteristics.ese variation characteristics are very similar to those of temperature field [17], indicating that the displacement is mainly induced by temperature field.In addition, Figure 3(b) shows the time history of D 1,u from 16:00 to 18:00 on May 20th, 2017.It is clear that the longitudinal displacements contain numerous small-amplitude dynamic displacements with random variability that are mainly induced by dynamic loads.
Because the monitored longitudinal displacement data contain static longitudinal displacements induced by the temperature field and dynamic longitudinal displacements induced by dynamic loads as well, the dynamic longitudinal displacements need to be extracted from the monitored data first.Note that the static and dynamic longitudinal displacements are in different frequency bands.
e static longitudinal displacement is greatly influenced by temperature changes, so the cycle period is about one day (i.e., the corresponding frequency is 1/86400 Hz).To ensure that the frequency of static longitudinal displacements is completely within the band range, the upper limit of the frequency band is designed to be 5/86400 Hz, namely (0, 5/86400] [13,14].
e frequency of dynamic longitudinal displacements is far higher than the frequency of static longitudinal displacements, and relevant research has indicated that the frequency bandwidth of dynamic longitudinal displacements is within the range of (5/86400, 0.5] [13,14].e reason for choosing 0.5 Hz as the upper limit is that the dynamic displacements with frequencies over 0.5 Hz are mainly induced by electromagnetic interference [13,14], which are not in the range of concern.
e wavelet packet decomposition method can use a pair of filters to decompose the monitored longitudinal displacement into different frequency bands scale by scale [18][19][20].Figure 4 shows the tree structure of the wavelet packet coefficients within 13 scales after layer-by-layer decomposition, where x j,m denotes the wavelet packet coefficient in the mth frequency band of the jth scale (m � 0, 1, 2, . . ., 2 j − 1), and f s denotes the sampling frequency with f s � 1 Hz in this research.e wavelet packet coefficient in the first frequency band (0, 1/2 14 ], which is close to (0, 5/86400], is selected and reconstructed to obtain static displacements; the wavelet packet coefficients in the other frequency bands (1/2 14 , 1/2], which are close to (5/86400, 0.5], are selected and reconstructed to obtain dynamic displacements. Taking Figure 3(b) as an example, the time history of D 1,u on May 20th, 2017 is decomposed by the wavelet packet method, and the decomposition result is shown in Figure 5.It is clear that the wavelet packet decomposition method can effectively decompose the monitored longitudinal displacement data into two components: the first component caused by temperature field and the second component (namely small-amplitude dynamic displacements) caused by dynamic loads.Furthermore, the cumulative displacements of two components are calculated to investigate the contribution of two components to cumulative displacements, which are 47.1 mm and 773.4 mm, respectively.Obviously, the second component caused by dynamic loads contributes the most to cumulative displacements.In detail, the dynamic loads mainly include train and wind loads.What should be mentioned is that there is no train passing the Dashengguan Yangtze River Bridge from 0 am to 6 am for bridge inspection, so the first component during this period is mainly induced by wind load; in addition, the second component during the rest of time is mainly induced by both wind and train loads, as shown in Figure 5

Shock and Vibration
because the influence of train load is obviously higher than wind load.Hence, the small-amplitude dynamic displacement caused by train load is specifically studied in this research.
e continuous accumulation of small-amplitude dynamic displacements will induce serious fatigue wear, which can severely decrease the service lifetime of bearings.us, it is necessary to carry out study on the cumulative characteristics of small-amplitude dynamic displacements.e CDD caused by a single passed train was calculated to evaluate the contribution of this train to the cumulative wear.e CDD caused by the kth passed train, which is denoted by N(k), is calculated by referring to the equation in Table 1, where  D(k) denotes the dynamic displacement caused by the kth passed train;  D j (k) and  D j+1 (k) denote the jth and (j + 1)th values in  D(k), respectively, where j � 1, 2, . . ., U − 1; U denotes the number of monitored values in  D(k).N(k) indicates the cumulative travel of bearing caused by the kth passed train.For example, one time history of dynamic displacement of B 3,u caused by one train is plotted in Figure 6, and the number of U in the train-passed segment of the plot is 48, which can be clearly captured, so the value of N(k) is 20.3 mm by calculating the CDD in the train-passed segment.e monitored data from August 10th to August 31st in 2017 is selected for analysis.A total of 3495 trains passed over the Dashengguan Yangtze River Bridge over this period, and each train corresponds to one value of N(k).
erefore, 3495 values of N(k) for each bearing at upstream and downstream sides (i.e., B 1,u ∼ B 3,u and B 1,d ∼ B 3,d ) are calculated as shown in Figure 7, where N i,u and N i,d denote a set of N(1), N(2), . . ., N(3495) { } for the bearings B i,u and B i,d , respectively, as shown in Table 1.It can be seen that all the CDDs show uniform and random characteristics.Furthermore, the average value of 3495 CDDs for each N i,u or N i,d , which is denoted by mean(N i,u ) or mean(N i,d ), respectively, is calculated by referring to the equations in Table 1, with the results shown in Figure 7.It can be seen that mean(N 3,d ) has the largest value about 7.020 mm while mean(N 1,u ) has the smallest value about 2.904 mm.

Cumulative Bearing Travel Caused by Many Passed
Trains.Each value of N i,u or N i,d corresponds to one passed train, and furthermore, the 3495 values for each B i,u or B i,d are accumulated to obtain CBT, which indicates the influence of the continuous actions of passed trains on bearing wear.e calculation formula for CBT is expressed as follows:

Signal
x 0,0    Shock and Vibration where M(p) denotes the CBT caused by p passed trains and p � 1, 2, . . ., 3495.e M(p) values at the upstream and downstream are shown in Figure 8, where M i,u denotes the CBT of B i,u and M i,d denotes the CBT of B i,d (i � 1, 2, 3).It can be seen that all the CBTs have obvious linear correlations with the number of passed trains, but the linear growth rates are different for each M i,u or M i,d , where the linear growth rate of M 3,d is the fastest, and its maximum CBT caused by 3495 passed trains is 24360 mm.

Probability Statistics Analysis
e number of monitored CDDs are not enough to satisfy the requirement of large sample for safety evaluation of bearing wear life.erefore, it is necessary to perform sampling simulation of CDDs using probabilistic statistics of the monitored CDDs to obtain adequate data.
3.1.Cumulative Probability Characteristics.Figure 7 shows that N i,u and N i,d contain obvious uniform and random characteristics, so N i,u and N i,d can be treated as statistical variables for probabilistic statistics analysis.In detail, the cumulative probability (CP) is used to describe the probabilistic statistics characteristics of CDDs.First, the cumulative probabilities of N i,u and N i,d are calculated based on probability statistics theory, and then the cumulative distribution function is fitted in a least-squares manner using the cumulative probabilities of N i,u and N i,d , which is illustrated in two steps as follows: (1) e cumulative probabilities of N i,u and N i,d are calculated as follows: where P(N i,u (j)) denotes the CP of the jth value N i,u (j) in N i,u ; P(N i,d (j)) denotes the CP of the jth

Variable
Calculation equation

Shock and Vibration
value N i,d (j) in N i,d ; NUM(N i,u ) denotes the total number of N i,u ; NUM(N i,d ) denotes the total number of N i,d ; NUM(N i,u < N i,u (j)) denotes the number of CDDs in N i,u that are smaller than N i,u (j); and NUM(N i,d < N i,d (j)) denotes the number of CDDs in N i,d that are smaller than N i,d (j).For instance, the CP of N 1,u is shown in Figure 9.
(2) e cumulative distribution function is fitted in a leastsquares manner using the cumulative probabilities of N i,u and N i,d .Considering that the true cumulative distribution function of N i,u and N i,d is unknown, three possible cumulative distribution functions (i.e., normal distribution (ND), Weibull distribution (WD), and general extreme value distribution (GEVD)) are chosen for least-square fitting, and then the cumulative distribution function with the best fitting effect is selected as the true cumulative distribution function of N i,u and N i,d .Taking N 1,u as an example, the fitting results of the three cumulative  6 Shock and Vibration distribution functions are shown in Figure 10.It can be seen that GEVD has the best fitting results, so GEVD is chosen to describe the cumulative probability characteristics of N 1,u .e formula is expressed as where G(N 1,u ) denotes the cumulative distribution function of N 1,u , which follows the generalized extreme value distribution; r denotes the shape parameter; a denotes the scale parameter; and b denotes the location parameter.e fitted values of these three parameters for G(N 1,u ) are r � 0.3058, a � 0.8163, and b � 2.0294.
Based on the method above, the best cumulative distribution functions of N 1,u − N 3,u and N 1,d − N 3,d are all GEVD, as shown in Figures 10(a)-10(f ).It can be seen that GEVD can well describe the cumulative probability characteristics of N 1,u − N 3,u and N 1,d − N 3,d .

Sampling Simulation Method.
Since N i,u and N i,d follow GEVDs, their GEVDs can be utilized for Monte Carlo sampling to obtain a demanded number of CDDs.Specifically, if the demanded number of CDDs is K, K values are 8 Shock and Vibration uniformly and randomly sampled from the interval (0, 1), and then the demanded number of CDDs can be simulated by substituting the K values into the cumulative probability R in the following equations: where N i,us and N i,ds denote the simulated CDDs of N i,u and N i,d , respectively.e values of r, a, and b have been calculated in the Section 3.2.However, N i,u and N i,d are embedded in G(N i,u ) and G(N i,ds ), so they cannot be directly solved from Equations (4a) and (4b).Considering that G(N i,us ) and G(N i,ds ) are monotonically increasing functions, N i,us and N i,ds can be determined by the Newton iteration method as follows: where N n+1 i,us and N n i,us denote the (n + 1)th and nth iteration results of N i,us , respectively; N n+1 i,ds andN n i,ds denote the (n + 1) th and nth iteration results of N i,ds , respectively; G ′ (N n i,us ) and G ′ (N n i,ds ) denote the first-order derivatives of GEVD (i.e., the probability density functions of GEVD).e iteration will terminate when the absolute difference between two adjacent iterations is less than 10 −2 .e last iteration result is treated as the values of N i,us and N i,ds .

e Sampling Simulation Result.
Based on the sampling simulation method above, the simulation results of 3495 N 1,us are shown in Figure 11(a).
e mean value of the simulated result is 2.820 mm, which is very close to the mean value of N i,u , i.e., 2.904 mm.Furthermore, the CP of the simulated result is calculated and then compared with the monitored CP, as shown in Figure 11(b).It can be seen that the simulated CP coincides with the monitored CP, which verifies the accuracy of the sampling simulation method.
Furthermore, the simulated 3495 values of N i,us and N i,ds are accumulated using Equation ( 1) to obtain the simulated CBT (i.e., M i,us and M i,ds ) as shown in Figure 12.It can be seen that the simulated M i,us and M i,ds present good linear growth characteristics with the number of passed trains.To verify the accuracy of the simulated linear growth rates, a comparison between the monitored and simulated linear growth rates was conducted as shown in Figure 12.It can be seen that the simulated linear growth rates have agreement with the monitored ones, which further verifies the accuracy of the sampling simulation method.

Safety Evaluation Analysis
e CBT relates to the bearing wear condition.If the CBTs exceed the wear limit in the service lifetime, the spherical steel bearing is totally damaged and cannot be used for longitudinal thermal expansion any longer.erefore, it is necessary to analyze whether the CBTs at bridge site can exceed the wear limit in the service lifetime.

Evaluation Method.
e failure probabilities of spherical steel bearings in the service lifetime (i.e., the probabilities of CBT over the wear limit in the service lifetime) are calculated as follows: or where P f denotes the failure probability; [M and N i,ds in the service lifetime, respectively, which can be obtained through sampling simulation as introduced in the Section 3.2; and J denotes the total number of passed trains in the service lifetime.F(M i,us ) and F(M i,us ) are unknown, which should be first determined by the following two steps: (1) By referring to the sampling simulation method in the Section 3.2, the values of N i,us or N i,ds in the service life are simulated through Monte Carlo sampling.
en, the values of N i,us or N i,ds are summed using Equation (6c) or (6d), respectively, to obtain one value of M i,us or M i,ds . is step is repeated W times to obtain W values of M i,us or M i,ds ; (2) e cumulative probability characteristics for the W values of M i,us or M i,ds are fitted using cumulative distribution function.e true cumulative distribution function of M i,us and M i,ds is unknown, so three possible cumulative distribution functions (namely, ND, WD, and GEVD) are considered, and the cumulative distribution function with the best fitting effect is selected as the true cumulative  After F(M i,us ) and F(M i,us ) are determined, [M i,us ] and [M i,ds ] are substituted in Equations (6a) and (6b) to obtain the value of P f .Furthermore, the reliability level β can be calculated using P f as follows: If the calculated reliability level β is larger than the target reliability level [β] (i.e., β > [β]), it can be inferred that the bearing is still in a good service state; otherwise (i.e., β ≤ [β]), it can be inferred that the bearing is in a poor service state and should be replaced in time.According to the Eurocode [21], the target reliability level [β] for newly built bridges with a 100-year design lifetime is 3.8, which corresponds to a failure probability P f � 7.2 × 10 −5 .However, the Dashengguan Railway Bridge has been in operation for 8 years, and it is not a newly built bridge.From the viewpoint of bridge reliability, a reduction in the uncertainties of the load and resistance parameters decreases the failure probability of existing bridge structures, indicating the possibility of adopting a lower reliability level for the evaluation of existing bridges rather than the reliability level that should be used for newly built bridges [22].According to Kotes' research [22], the target reliability level [β] with a 90-year remaining lifetime is 3.773, so the target reliability level [β] of the Dashengguan Railway Bridge with a 92-year remaining lifetime is 3.778 using the linear interpolation calculation, namely, [β] � 3.778.

Evaluation Results
. Figure 8 reveals that M 3,d has the fastest growth, so M 3,d is taken as an example to calculate its failure probability and evaluate the wear condition of B 3,d in the service lifetime.Before evaluation, three parameters should be determined, namely, J, [M 3,ds ], and W. For the parameter J, the designed service life of the Nanjing Dashengguan Yangtze River Bridge is 100 years, and the number of passed trains in 2017 was 84343, so the value of J is 8434300; for the parameter [M 3,ds ], the China code specifies that the linear wear rate of PTFE in a spherical steel bearing should be less than 15 μm/km and that the thickness of PTFE should be more than 7 mm and less than 8 mm (i.e., only a 1 mm wear thickness is allowed for an 8 mm thickness of PTFE) [23], so the value of [M 3,ds ] is 66.7 km for the spherical steel bearing; for the parameter W, a larger W can provide better evaluation results, and the value of W is 1000 for evaluation in this research.
After the three parameters are determined, the bearing wear life is evaluated by the following seven steps:  6) e reliability level β is 5.612 calculated by Equation (7), which is larger than the target reliability level 3.778, indicating that bearing B 3,d does not reach the wear limit in the service life.(7) Because the linear growth rate of CBT for B 3,d is faster than those of the other spherical steel bearings, it can be inferred that the other spherical steel bearings do not reach the wear limit as well in the service life.

Conclusions
Based on the monitored dynamic displacement data of steel truss arch bridge bearings, the CDD under the action of a single train and the CBT under the continual actions of many passed trains are investigated.Furthermore, the probability statistics and the Monte Carlo sampling simulation method for CDD are studied, and the safety evaluation method for the bearing wear life is proposed using a reliability index regarding the failure probability of monitored CBT over the wear limit in the service lifetime.e main conclusions are drawn as follows: (1) rough monitored data analysis, the monitored longitudinal displacement data mainly consist of temperature-induced static displacements and traininduced dynamic displacements.e wavelet packet decomposition method can effectively extract the small-amplitude dynamic displacements caused by train loads from the monitored data.All the CDDs contain uniform and random characteristics, and all the CBTs have an obvious linear correlation with the number of passed trains, but their linear growth rates are different.(2) rough probability statistics analysis, GEVD can well describe the cumulative probability characteristics of CDD. e Monte Carlo sampling simulation was performed using GEVD to simulate CDD and CBT, and the results show that the simulated cumulative probabilities of CDD have agreement with the monitored ones of CDD, and the simulated linear growth rates of CBT have agreement with the monitored ones of CBT.(3) rough safety evaluation analysis, the safety evaluation method for the bearing wear life is proposed using a reliability index, and it is finally judged whether the CBTs can exceed the wear limit in the service lifetime.e evaluation results show that all the spherical steel bearings do not exceed the wear limit in the service life.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Figure 1 :
Figure 1: e components of spherical steel bearing.
(b).It can be inferred from Figure5(b) that train load is the major influence factor for cumulative displacements

Figure 2 :
Figure 2: Locations of the longitudinal displacement sensors m.

Figure 4 :
Figure 4: Tree structure of the wavelet packet coefficients.

Figure 3 :
Figure 3: e monitored result of D 1,u .(a) From March to October in 2017.(b) On May 20th, 2017.

Figure 5 :
Figure 5: Decomposition result of monitored longitudinal displacement.(a) Component caused by temperature field.(b) Component caused by dynamical loads.

Figure 6 :
Figure 6: e dynamic displacement of B 3,u caused by one train.

10
Shock and Vibration distribution function of M i,us and M i,ds after leastsquare fitting.

Figure 13 :
Figure 13: e 1000 values of M 3,ds and their cumulative probabilities.(a)1000M 3,ds .(b) Cumulative probability of M 3,ds and its fitted results.

Table 1 :
Calculation of variables N(k), N i,u , N i,d , mean(N i,u ), and mean(N i,d ).
NUM(N i.d ): total number of N i,d NUM(N i,u < N i,u (j)): number of CDDs in N i,u that are smaller than N i,u (j) NUM(N i,d < N i,d (j)): number of CDDs in N i,d that are smaller than N i,d (j) CP of the jth value N i,u (j) in N i,u P(N i,d (j)):CP of the jth value N i,d (j) in N i,d P