Shock and Vibration
Volume 2016 (2016), Article ID 1450895, 19 pages
http://dx.doi.org/10.1155/2016/1450895

Research Article

Nonstationary Random Vibration Performance of Train-Bridge Coupling System with Vertical Track Irregularity

Xiaozhen Li,¹ Yan Zhu,1,2 and Zhibin Jin¹

¹School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China
²School of Civil and Environmental Engineering, The University of New South Wales, Sydney, NSW 2052, Australia

Received 6 August 2015; Revised 2 October 2015; Accepted 8 October 2015

Academic Editor: Vadim V. Silberschmidt

Copyright © 2016 Xiaozhen Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In order to study the random vibration performance of trains running on continuous beam bridge with vertical track irregularity, a time-domain framework of random analysis on train-bridge coupling system is established. The vertical rail irregularity is regarded as a random process. A multibody mass-spring-damper model is employed to represent a moving railway car and the bridge system is simulated by finite elements. By introducing the pseudo excitation algorithm into the train-bridge interaction dynamic system, expressions of the mean value, standard deviation, and power spectral density of the nonstationary random dynamic responses of bridge and vehicles are derived. Monte-Carlo simulations are implemented to validate the presented method. A comprehensive analysis of the train-bridge coupling system with vertical track irregularity is conducted focusing on the effect of the randomness of the vertical rail irregularity on the dynamic behavior of the running train and the three-span continuous concrete bridge. Moreover, stochastic characteristics of the indicator for assessing the safety and the riding quality of the railway cars running on continuous beam bridge are carried out, which may be a useful reference in the dynamic design of the bridge.

1. Introduction

The passage of moving trains leads to the dynamic amplification of a bridge; meanwhile, the motion of a bridge deck affects the kinetic characteristic of trains. The study of train-bridge coupling system commenced from the mid-19th century [1]. After that, a considerable number of investigations have been conducted on the dynamic behaviour of bridges and trains [2–15]. The majority of the literature has been devoted to the analysis of the bridge and trains vibrations when the system is subjected to deterministic loads. In conventional dynamic analysis of train-bridge coupling system, a few samples of the rail irregularity are adopted which are obtained from some transformation by use of the power spectral density (PSD) of the track irregularity. Strictly speaking, most obtained results may be sample values. Therefore, the rail irregularity, being regarded as one kind of excitation of a train-bridge coupling system, is essentially a random process and its effects on the dynamic performance of the train-bridge coupling system should be extensively investigated.

The key issue for determining the statistical characteristics of random dynamic responses of the railway bridge system, when high-speed trains are passing through, is the complexity of the coupling system. For instance, nonlinearity exists at the contacts between the wheels and the rail, and the implementation of the theories of nonstationary random vibrations to the coupling system is difficult. A few researchers [16–18] have attempted to investigate random vibration responses of the bridge using the so-called moving load model. However, without the train model, random vibration responses of the train and the interactions between the bridge and the train cannot be developed.

Monte-Carlo method has been adopted in random dynamic analysis of the train-bridge coupling system. Using Monte-Carlo simulations, many samples of the track irregularity are obtained and then analyses according to different parameters of bridges and trains are performed [19, 20]. Xia and Zhang selected 20 samples of the rail irregularity, and then the kinetic properties of a single span box beam and the vehicles were obtained by using the samples as the excitation of the system [21]. The calculations from different samples show great differences. Majka and Hartnett [22] investigated the dynamic responses of an existing railway bridge subjected to moving trains and discussed the effect of different service trains and the random track irregularity. Monte-Carlo simulations are easily carried out but require enough number of samples to ensure the convergence of calculations, which are very time consuming and sometimes are not implementable in engineering practice especially for complex systems or problems.

In addition to performing Monte-Carlo simulations, other approaches have also been developed to analyse stochastic dynamic problems. Jin [23] employed the covariance method to calculate both displacement and accelerations of the coupling system. Pseudo excitation method (PEM) proposed by Lin et al. [24] has been widely used in linear analysis of random vibration owing to its accuracy and high efficiency. In PEM, the random excitation is transformed into pseudo excitation; then, stochastic equations of motion can be solved using numerical integration. Zhang et al. [25] investigated the influences of train speed and the three kinds of track irregularities on the three-dimensional train-bridge system based on PEM.

Many of the investigations on the stochastic characteristic of the bridge are on the simply supported beam. The integrity of the deck of multispan continuous beam is better than that of simple supported girder, so the continuous beam is widely used in Chinese high-speed railway bridge. However, most of the research on the dynamic response of the multispan continuous beam bridge is deterministic [6, 26, 27]. In this paper, the random vibration model of the time varying train-bridge coupling system is developed. The vehicle is treated as a combination of several rigid bodies which are connected by springs and dampers. The bridge is modelled using the finite element method. The vertical rail irregularity which is treated as internal excitation is constructed as pseudo excitation. Then, the random vibration equations of the bridge system and the vehicle system are integrated into a coupling system based on the equilibrium of the forces and the deformation compatibility. The stochastic characteristics of the bridge and the vehicles are obtained in the time domain by an iterative procedure at each time step. We focus on the effect of the randomness of the vertical rail irregularity on the dynamic responses of the continuous beam bridge and the running rail cars. Furthermore, the dispersal of the safety and the riding comfort of the vehicles are investigated.

The remaining parts of this paper are organised as follows. Section 2 presents the train-bridge coupling system used in this study. Section 3 summarises the pseudo excitation method for time-invariant system. Pseudo load and solution method for time-variant train-bridge system are derived in Sections 4 and 5. Section 6 investigates the safety and riding quality of the railway cars. The developed model is validated by using Monte-Carlo method in Section 7. A comprehensive analysis of the train-bridge coupling system with vertical track irregularity is conducted in Section 8. Conclusions are drawn in Section 9.

2. Train-Bridge Coupling System

In the train-bridge coupling system, a multibody model is used to represent a moving railway car and a spatial beam element is adopted to simulate the bridge system as shown in Figure 1. Here, the subscripts , , and represent the car body, front bogie, and back bogie, respectively; and represent the vertical displacement and rotation about -axis; , , and denote the mass of car body, front bogie, and back bogie, respectively; the primary suspension system is demonstrated by four linear springs with stiffness and four dampers with damping coefficient ; the secondary suspension system is represented by two linear springs with stiffness and two dampers with damping coefficient ; , , and denote the moment of inertia of the car body, front bogie, and back bogie, respectively.

Figure 1: Model of train-bridge coupling system.

In this study, the railway car is modelled by a multibody connected by the primary and secondary springs with 6 degrees of freedom as shown in Figure 1. The dynamic equations of the vehicle system can be written in the matrix form as

where subscript represents the vehicle. , , and are the mass, damping, and stiffness matrices of the vehicle [28], , , and are the acceleration, velocity, and displacement of the vehicle, respectively, and

is the interaction vector applied to the vehicle due to the vibration of the track; that is,

where and denote force and moment, respectively; subscripts and denote vertical and horizontal direction, respectively. Since the car body and track have no direct contact, and equal zero.

A spatial beam element model is applied to simulate the bridge system. Equations of motion of the bridge system can then be expressed as

where subscript denotes the bridge; , , and are the mass, damping, and stiffness matrices of the bridge [21, 28]; , , and are the acceleration, velocity, and displacement of the bridge system with

is the interaction vector applied to the bridge deck due to the vibration of the moving train.

The interaction between the train and bridge can be calculated through the compatibility condition of the displacement and equilibrium equations of the wheel/rail forces (see Figures 2 and 3).

Figure 2: Compatibility condition of the displacement.

Figure 3: Wheel/rail forces.

In Figure 2, - is the coordinate system of the bridge model, - is the coordinate system of the wheelset model, and - is the coordinate system of the track model. Assuming the wheelset and the track are constantly connected, the displacement of the track centre equals displacement of the wheelset , and

where subscript denotes the bridge; is the vertical displacement of the bridge; denotes the eccentricity of the track centre; denotes the rotation along the axis of the cross section; is the vertical rail irregularity.

In Figure 3, subscript denotes the wheelset; and are the suspension force owing to the primary spring and the vertical wheel/rail force exerted on every wheelset, respectively. The mechanical equilibrium equation of the wheelset in vertical direction can be expressed as

(7) can be rewritten as

The vertical wheel/rail force vector exerted on the rail, that is, (), and the vertical wheel/rail force vector exerted on wheelset () are action and reaction, respectively. Hence, they have the same magnitudes but are in the opposite direction. Meanwhile, assuming that is fully transferred to the bridge element, the single wheel load vector of the bridge element can be written as

where is the lateral wheel/rail force exerted on the rail; is the moment around -axis exerted on the rail; and is the vertical distance between the centroid of the bridge section and the centre of the rail.

After converting the single wheel load vector () to the nodal load vector, in (10) can be written as

where is the total number of the wheelsets.

3. Pseudo Excitation Method (PEM) for Time-Invariant System

Equation (11) shows the principle of PEM [29]. Here, denotes the frequency response function. As shown in (a), if the random force is equal to , then the corresponding harmonic response can be expressed as. Hence,

Pseudo excitation shown in (b) is constructed [24] as follows:

The corresponding pseudo response is given by

where denotes the pseudo expression of . The autopower spectral density of the dynamic response can be taken as

where superscript “” denotes the complex conjugation.

If the two pseudo variables are obtained, the cross-power spectral density can be taken as

Accordingly, the PSD matrix of the dynamic responses can be expressed as

If the dynamic responses are linear with the excitation, the pseudo excitation method can be applied. For linear structures, the stationary random responses can be determined by means of harmonic vibration analyses, while the nonstationary random responses can be computed by means of an arbitrary direct integration scheme. Therefore, all these random response analyses become simple by using the PEM [30].

4. PEM Method for Time-Variant Train-Bridge Coupling System

In this paper, vertical track irregularity is considered as random progress with zero mean. The kinetic equation of the coupled system can be expressed as follows:

where , , and are the mass, damping, and stiffness matrix, respectively; , , and are the acceleration, velocity, and displacement vector, respectively; is the deterministic load vector mainly generated by the weight of the vehicles; and is the random load vector caused by the vertical track irregularity.

Equation (17) can be solved by using Duhamel integration with respect to over the period when the train moves on the bridge; it obtains

The mean value of the responses can be written as

As, (19) can be simplified as

The autocovariance of the responses can be calculated by

Substituting (18) and (19) into (21) yields

can be written as

where is the transition matrix; is the vertical profile irregularity vector; that is, ( is the sum of the wheelsets and ,  ).

Substituting (23) into (22) yields

Here, the Wiener-Khintchine relation is adopted:

Then, the autocovariance of the responses becomes

where superscript denotes the complex conjugation; is the PSD of the vertical rail irregularity; and

Assuming that

(26) can be rewritten as

Let ; we have

From (25) and (30), it can be obtained that

where is the autopower spectral density of the dynamic responses.

Furthermore,

From (28), can be regarded as the response vector caused by the excitation vector ; thus, the pseudo excitation of the coupling system can be instituted [28, 29]:

where superscript denotes “pseudo.”

Therefore, the responses vector induced by the pseudo load vector as indicated in (33) can be obtained by

So the auto-PSD matrix of the response vector can be expressed as [28]

From the principle of the pseudo excitation method and (3), the auto-PSD matrix of the response vector can also be written as

Therefore,

The expression of (32) is the same as (37). Consequently, the assumption of the pseudo load vector expressed in (33) is logical.

5. Solutions to the Random Characteristics of Train-Bridge Coupling System

The problems presented in the above section can be solved by applying the numerical integration and programming to the following analysis process:

(1)Average value can be obtained by using (20).(2)PSD and SD of the dynamic responses of the train-bridge system are obtained as shown in the flowchart (Figure 4).

Figure 4: Flowchart of the dynamic responses of the train-bridge system.

6. Safety and Riding Quality of the Railway Cars

When normal contact conditions between the wheels and the rail are satisfied, the safety of the running vehicles can be guaranteed. However, under the worst load combination, the derailment or overturn accident may happen. Several parameters can be used to evaluate the safety of the trains, such as derailment coefficient, wheel unloading rate, and lateral swing force. In this paper, we choose the wheel unloading rate as the safety index of the railway cars. Riding quality can show the influence of the vibration of the train on the comfort of the passengers. In ISO2631, Janeway and Jacking comfort curves and Sperling index are usually employed to evaluate the riding quality of the running railway cars. In this study, Sperling index is used to assess the passengers’ riding comfort.

6.1. Wheel Unloading Rate

When a train is moving with high speed, the dynamic wheel weight will decrease during the upward movement. Meanwhile, the wheelset may derail owing to the lateral relative displacement between the bridge deck and the wheelset in spite of small lateral wheel-rail interaction force [31]. Wheel unloading rate can be defined as the dynamic wheel load reduction amount divided by the average static wheel load, that is, [21].

6.2. Sperling Index

Sperling index of the comfort can be expressed as [21]

where is acceleration (cm/s²); is vibration frequency (Hz); is the correction coefficient of frequency; that is,

7. Verification of the Presented Computational Model

In this section, a single span simple supported bridge is adopted to confirm the rationality of the method presented in Section 4, and only one ICE high-speed motor train [31] passes through the bridge with a speed of 220 km per hour. The span is 32 m. The cross-section properties of the bridge are taken as follows: elasticity modulus is 3.45 × 10¹⁰ N/m²; Poisson’s ratio is 0.2; ratio of damping is 0.02; area is 7.5 m²; inertia moment in , , and direction is 35.4 m⁴, 9 m⁴, and 70.158 m⁴, respectively.

The German low-interference track irregularity spectrum is applied; that is,

where is the spatial frequency; is roughness constant of track; that is, = 4.032 × 10⁻⁷; and are cut-off frequency; that is, and .

Wavelength coverage of the rail irregularity is from 1 m to 80 m. Eccentricity of the line equals 2.5 m; Rayleigh damping is adopted in this study. The railway cars travel about 100 m on the same line as on the bridge before they pass through the bridge and then travel about 100 m again after getting off the bridge. The time step of the integration procedures is 0.0005 s.

The trigonometric series method is adopted [21] in the modelling of the rail irregularity; that is,

where is the sequence of the vertical track irregularity; , where and are, respectively, the upper and the lower limits of the spatial frequency:

is the th phase and is in accordance with the principle of uniform distribution.

To verify the results obtained by the presented method, 500, 5000, 10000, and 20000 Monte-Carlo simulations are implemented. Figures 5 and 6 show that the results calculated by the presented method, denoted by SRCTB, agree well with those obtained by using Monte-Carlo simulation method; however, the latter requires more time in the tentative calculation efforts. When the number of samples reaches 500, the statistical characteristics of the bridge using Monte-Carlo method are satisfactory as shown in Figure 5. When sample number increases to 10000, the mean values of the car body and the first bogie obtained by using Monte-Carlo method agree very well with the calculations by applying (20). As to the standard deviation (SD) of the dynamic responses about the vehicle, 20000 samples should be required as shown in Figure 6. In Figures 5 and 6, the -axis title “Distance (m)” means the distance between the first wheelset and the starting point of the bridge.

(a)

(b)

(c)

(d)

Figure 5: Statistical characteristics of the bridge response in the midpoint. (a) Mean value of the vertical displacement. (b) Mean value of the acceleration. (c) SD of the vertical displacement. (d) SD of the acceleration.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 6: Statistical characteristics of the railway car. (a) Mean value of the vertical displacement of the car body. (b) Mean value of vertical acceleration of the car body. (c) Mean value of the vertical displacement of the front bogie. (d) Mean value of the vertical acceleration of the front bogie. (e) SD of the vertical displacement of the car body. (f) SD of the vertical acceleration of the car body. (g) SD of the vertical displacement of the front bogie. (h) SD of the vertical acceleration of the front bogie.

The comparison of the robustness of the dynamic responses and the computation efficiency between Monte-Carlo method and the presented method (denoted by SRCTB) is listed in Table 1. We find that SRCTB method has excellent robustness and computational efficiency.

Table 1: Comparison of robustness and computational efficiency.

Method
Computational time (h)
Robustness
Efficiency

Bridge responses
Vehicle response

Monte-Carlo
(500 samples)
0.9
Fair
Poor
Fair

Monte-Carlo
(5000 samples)
9.8
Fair
Poor
Poor

Monte-Carlo
(10000 samples)
19.4
Excellent
Fair
Poor

Monte-Carlo
(20000 samples)
27.8
Excellent
Excellent
Poor

SRCTB
0.5
Excellent
Excellent
Excellent

8. Numerical Examples

In this section, a three-span nonuniform prestressed continuous beam bridge along the Chinese high-speed railway line is chosen to demonstrate the random properties of the train-bridge coupling system. Vehicles adopted are ICE high-speed train [30]. The bridge model is shown in Figure 7.

Figure 7: Three-span continuous beam bridge.

The cross section of the beam is a single cell box girder as shown in Figure 8. The elastic modulus is 3.55 × 10¹⁰ N/m² and Poisson’s ratio is 0.2.

Figure 8: Cross section of the beam.

The wavelength coverage of the vertical rail irregularity is from 1 m to 80 m; the vehicle velocity is 275 km/h; eccentricity of the rail line equals 2.5 m; Rayleigh damping is adopted; the trains travel about 100 m on the same line as on the bridge before they pass through the bridge, and they travel about 100 m again after getting off the bridge.

8.1. Stochastic Characteristic Analysis of the Continuous Beam Bridge

As shown in Figure 9(a), vehicles loaded in the first span induce downward displacement in the first span and upward displacement in the third span.

(a)

(b)

Figure 9: Displacement of the midpoint in the side span. (a) Mean value. (b) SD.

From Figures 9(b), 10(b), and 11(a), it can be observed that the peak values of the standard value of the displacement in the midpoint of the middle span are larger than those in the midpoint of the two side spans. However, as to the peak values of the standard value of the acceleration, the opposite calculations are obtained.

(a)

(b)

Figure 10: Acceleration of the midpoint in the side span. (a) Mean value. (b) SD.

(a)

(b)

Figure 11: Dynamic responses of the midpoint in the middle span. (a) Mean value and SD of the displacement. (b) (a) Mean value and SD of the acceleration.

As shown in Figures 9, 10, and 11, the absolute values of the mean about the bridge vertical dynamic responses are much larger than the SD of them. From (20), we believe that the vertical responses of the midpoint in the three spans are mainly affected by the certain load, that is, the wheel load.

As shown in Figure 12, the majority of the PSDs about the displacement in the midpoint of the first span and the same place in the third span cover the same frequency range. As to the corresponding place in the middle span, its PSD is larger than that of the side spans.

(a)

(b)

(c)

Figure 12: PSD of the bridge displacement of (a) midpoint of the first span, (b) midpoint of the middle span, and (c) midpoint of the third span.

From Figure 13, it can be seen that the acceleration PSDs of the midpoint in the three spans do not have similar frequency range.

(a)

(b)

(c)

Figure 13: PSD of the bridge acceleration in the midpoint of (a) first span, (b) middle span, and (c) third span.

The PSDs have different varying rules with each other as shown in Figures 12 and 13. The nonstationary features do exist in the dynamic responses of bridge.

8.2. Stochastic Characteristic Analysis of Vehicles

From Figures 14–17, it can be found that the mean values of dynamic responses are smaller than the standard deviations of them of railway cars. In Figure 14, the mean and the SD of the displacement of the two bogies vary with time in the same pattern but the values of the front bogie are smaller.

(a) 

(b)

Figure 14: Displacement of the first railway car. (a) Mean value. (b) Standard deviation.

(a)

(b)

Figure 15: Displacement of the car bodies. (a) Mean value. (b) Standard deviation.

(a)

(b)

Figure 16: Displacement of the bogies. (a) Mean value. (b) Standard deviation.

(a)

(b)

Figure 17: Acceleration of the vehicles. (a) Mean value. (b) Standard deviation.

The mean and SD of the displacement as to the four car bodies vary with time in almost the same pattern. When all the four railway cars are on the bridge deck and the third vehicle passes around 0.4 times of the middle span length, the SD of the displacement about the car bodies reach their maximum value (Figure 15). Table 2 gives the position information where each vehicle is the SD of the displacement about the car body reaching its maximum.

Table 2: Position corresponding to the maximum SD values.

Vehicle number
Position of the vehicle

1
0.6 times of the middle span length

2
0.6 times of the middle span length

3
0.4 times of the middle span length

4
0.1 times of the middle span length

Note: the distance is counted from the starting point of the middle span.

Figures 16 and 17 show the mean and the standard deviation of the dynamic responses of the bogies varying with time in a similar pattern, but with different peak value and corresponding different position on the deck. The mean and the SD of the acceleration about the car bodies are smaller than the bogies.

As shown in Figure 18, at every frequency point, the PSDs of number 1 railway car vary with time when the vehicle passes through the bridge in the beginning; then the values are close to stationary. The PSD values decrease with the frequency when the time is fixed, shown in Figure 18(a). We can see the same varying pattern in the values of PSD about the other three railway cars as shown in Figures 18(b), 18(c), and 18(d). The PSD of the displacement of the bogies and the PSD of the acceleration of the vehicles have similar varying pattern to the PSD of the displacement of the car bodies.

(a)

(b)

(c)

(d)

Figure 18: PSD of the displacement of the car bodies in the centre of the (a) first car body, (b) second car body, (c) third car body, and (d) fourth car body.

8.3. Discussion of the Safety and Riding Quality of the Railway Cars

As for the vertical wheel/rail force of the first railway car, the mean values of four wheelsets are around their static wheel load, that is, 156.8 kN (motor car) or 143.08 kN (trail car). As shown in Figures 19(a) and 19(b), we find that the vertical wheel/rail force is affected by the static wheel load. The mean and the SD of the wheel load reduction rate are very small and slightly vary with time as shown in Figures 19(c) and 19(d). We also find that the wheel load reduction rate is affected by the static wheel load from Figures 19(c) and 19(d).

(a)

(b)

(c)

(d)

Figure 19: Vertical force and wheel load reduction rate as for number 1 railway car. (a) Mean value of the vertical wheel/rail force. (b) SD of the vertical wheel/rail force. (c) Mean value of the wheel load reduction rate. (d) SD of the wheel load reduction rate.

As shown in Table 3, the mean value and the standard deviation of the vertical Sperling index of the car body are smaller than those of the bogie, owing to the dynamic loads not being exerted directly on the car body. The Sperling index has high-level dispersal, because its mean value is smaller than its SD.

Table 3: Sperling index of the vehicles.

Vehicle number
Car body
Bogie

Mean
SD
Mean
SD

1
0.751
0.858
1.081
1.763

2
0.655
0.823
1.014
1.948

3
0.659
0.823
1.033
1.938

4
0.682
0.869
1.014
1.795

9. Conclusions

Random vibration analysis of train-bridge coupling system is implemented by extending the pseudo excitation method (PEM). Expressions of the mean value, standard deviation, and power spectral density of the nonstationary random dynamic responses of bridge and vehicles are derived. Evaluation indicators for assessing the safety and riding quality of the railway cars are developed. Results produced by the redeveloped PEM agree with those generated by the Monte-Carlo simulation method, but the presented method requires less computational efforts.

By taking a high-railway bridge with three-span continuous beam as a case study, the whole histories of an ICE3 train running on the bridge subjected to vertical random rail irregularity are investigated. The following observations and conclusions may be drawn from the present numerical study:

(1)The location of vehicle affects the bridge responses.(2)The randomness of the dynamic responses of the midpoint in the middle span is different from that in the side span. The vertical displacement of the midpoint is mainly affected by the certain load, that is, the wheel load.(3)Nonstationary features can be found from the PSD of the responses of the bridge.(4)Randomness exists in the dynamic responses of the vehicles.(5)The value of the vertical wheel/rail force and wheel load reduction is mainly affected by the static wheel load.(6)The mean value and the standard deviation of the vertical Sperling index of the car body are smaller than those of the bogie due to the vibration suppression caused by the bogie. The Sperling index has high-level dispersal, which should be considered in the analysis and design of the coupling system.

Random analysis of train-bridge coupling system with rail irregularities is a very complex problem, which is related to the simulation of the random process, the implementation of the theory of random vibration on the time-varying coupling system, efficiency of the solution in the actual engineering application, and so on. In this paper, only a preliminary study on train-bridge coupling system with vertical rail irregularity is performed and illustrated by a case study. The algorithm presented in this paper and conclusions from the numerical example may be used as some reference for the dynamic design of the continuous beam bridge.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 51308470), Fundamental Research Funds for the Central Universities (Grant no. SWJTU12CX062), and Australian Research Council through Discovery Projects.

References

R. Willis, Appendix to the Report of the Commissioners Appointed to Inquire into the Application of Iron to Railway Structures, H.M. Stationary Office, London, UK, 1849.

L. Frýba, “Effect of a force moving at variable speed along a beam,” Acta Technica CSAV, vol. 18, no. 1, pp. 54–68, 1973.

L. Frýba, Dynamics of Railway Bridges, Telfor, London, UK, 2nd edition, 1996.

L. Frýba, Vibration of Solids and Structures under Moving Loads, Thomas Telford, 3rd edition, 1999.

C. H. Lee, M. Kawatani, C. W. Kim, N. Nishimura, and Y. Kobayashi, “Dynamic response of a monorail steel bridge under a moving train,” Journal of Sound and Vibration, vol. 294, no. 3, pp. 562–579, 2006.

J. W. Kwark, E. S. Choi, Y. J. Kim, B. S. Kim, and S. I. Kim, “Dynamic behavior of two-span continuous concrete bridges under moving high-speed train,” Computers and Structures, vol. 82, no. 4-5, pp. 463–474, 2004.

W. Lacarbonara and V. Colone, “Dynamic response of arch bridges traversed by high-speed trains,” Journal of Sound and Vibration, vol. 304, no. 1-2, pp. 72–90, 2007.

C. H. Lee, C. W. Kim, M. Kawatani, N. Nishimura, and T. Kamizono, “Dynamic response analysis of monorail bridges under moving trains and riding comfort of trains,” Engineering Structures, vol. 27, no. 14, pp. 1999–2013, 2005.

Y. L. Xu, N. Zhang, and H. Xia, “Vibration of coupled train and cable-stayed bridge systems in cross winds,” Engineering Structures, vol. 26, no. 10, pp. 1389–1406, 2004.

Y. B. Yang, C. L. Lin, J. D. Yau, and D. W. Chang, “Mechanism of resonance and cancellation for train-induced vibrations on bridges with elastic bearings,” Journal of Sound and Vibration, vol. 269, no. 1-2, pp. 345–360, 2004.

Y. B. Yang, K. C. Chang, and Y. C. Li, “Filtering techniques for extracting bridge frequencies from a test vehicle moving over the bridge,” Engineering Structures, vol. 48, pp. 353–362, 2013.

Y.-B. Yang and Y.-S. Wu, “A versatile element for analyzing vehicle-bridge interaction response,” Engineering Structures, vol. 23, no. 5, pp. 452–469, 2001.

H. Zeng and C. W. Bert, “Dynamic amplification of bridge/vehicle interaction: a parametric study for a skewed bridge,” International Journal of Structural Stability and Dynamics, vol. 3, no. 1, pp. 71–90, 2003.

C. Xia, H. Xia, N. Zhang, and W. Guo, “Effect of truck collision on dynamic response of train-bridge systems and running safety of high-speed trains,” International Journal of Structural Stability and Dynamics, vol. 13, no. 3, Article ID 1250064, 2013.

S.-I. Kim and N.-S. Kim, “Dynamic performances of a railway bridge under moving train load using experimental modal parameters,” International Journal of Structural Stability and Dynamics, vol. 10, no. 1, pp. 91–109, 2010.

L. Frýba, “Non-stationary response of a beam to a moving random force,” Journal of Sound and Vibration, vol. 46, no. 3, pp. 323–338, 1976.

A. Garinei and G. Risitano, “Vibrations of railway bridges for high speed trains under moving loads varying in time,” Engineering Structures, vol. 30, no. 3, pp. 724–732, 2008.

G. Ricciardi, “Random vibration of beam under moving loads,” Journal of Engineering Mechanics, vol. 120, no. 11, pp. 2361–2378, 1994.

F. T. K. Au, J. J. Wang, and Y. K. Cheung, “Impact study of cable-stayed railway bridges with random rail irregularities,” Engineering Structures, vol. 24, no. 5, pp. 529–541, 2002.

H. Xia, G. De Roeck, H. R. Zhang, and N. Zhang, “Dynamic analysis of train-bridge system and its application in steel girder reinforcement,” Computers and Structures, vol. 79, no. 20-21, pp. 1851–1860, 2001.

H. Xia and N. Zhang, Dynamic Interaction of Vehicles and structures, Science Press, Beijing, China, 2005.

M. Majka and M. Hartnett, “Dynamic response of bridges to moving trains: a study on effects of random track irregularities and bridge skewness,” Computers & Structures, vol. 87, no. 19-20, pp. 1233–1252, 2009.

Z. B. Jin, Train-track-bridge system and train-bridge stochastic vibration [Ph.D. thesis], Southwest Jiaotong University, 2007.

J. H. Lin, W. Shen, and F. W. Williams, “Accurate high-speed computation of non-stationary random structural response,” Engineering Structures, vol. 19, no. 7, pp. 586–593, 1997.

Z. C. Zhang, J. H. Lin, Y. H. Zhang, W. P. Howson, and F. W. Williams, “Non-stationary random vibration analysis of three-dimensional train-bridge systems,” Vehicle System Dynamics, vol. 48, no. 4, pp. 457–480, 2010.

V. De Salvo, G. Muscolino, and A. Palmeri, “A substructure approach tailored to the dynamic analysis of multi-span continuous beams under moving loads,” Journal of Sound and Vibration, vol. 329, no. 15, pp. 3101–3120, 2010.

D. Mu and D.-H. Choi, “Dynamic responses of a continuous beam railway bridge under moving high speed train with random track irregularity,” International Journal of Steel Structures, vol. 14, no. 4, pp. 797–810, 2014.

Y. Zhu, Train-bridge coupling system stochastic theory and application [Ph.D. thesis], Southwest Jiaotong University, 2011.

J. H. Lin and Y. H. Zhang, Pseudo Excitation Method in Random Vibration, Science Press, Beijing, China, 2006.

J. Lin, Y. Zhao, and Y. Zhang, “Accurate and highly efficient algorithms for structural stationary/non-stationary random responses,” Computer Methods in Applied Mechanics and Engineering, vol. 191, no. 1-2, pp. 103–111, 2001.

W. M. Zhai, Vehicle-Track Coupling Dynamic, Science Press, Beijing, China, 2007.