Published: 30 June 2014

Analytical solution for dynamic response of curved rail subjected to moving train

Hougui Zhang1
Weining Liu2
Kefei Li3
Dan Su4
1, 2, 3School of Civil Engineering, Beijing Jiaotong University, Beijing, 100044, China
3Beijing MTR Construction Administration Corporation, Beijing, 100037, China
4, 1Departament of Civil and Environment Engineering, Rutgers, The State University of New Jersey, New Jersey, 08854, USA
Corresponding Author:
Hougui Zhang
Views 123
Reads 56
Downloads 1213

Abstract

The objective of this paper is to present an analytical solution for the dynamic response of curved rail. The detailed solution was derived for the out-of-plane vibration response of periodically supported curved Timoshenko beam subjected to moving loads. The accuracy of the solution was validated based on the results from previous study. Furthermore, the train/track interaction model was introduced into this solution to calculate the rail dynamic response in Beijing metro. The results presented herein indicated that the solution provided accurate results in comparison with both the result from previous study and the measurement data collected from Beijing metro and it can be used to derive the dynamic response for similar situations.

1. Introduction

Rail corrugation, which is found to be very common for curved rails, is becoming an international issue for the railway industry worldwide. In order to investigate the rail corrugation due to the train/track dynamic interaction that was found in Beijing metro, the solution for the dynamic responses of curved rail under moving train loads is needed. In order to focus on the dynamic behavior of track structure, the rail was simplified as a periodically supported Timoshenko beam along horizontally curve, while other external factors, such as the super elevation and the wheel/rail contact condition were not considered in this paper.

Previously, some researchers have been conducted for the out-of-plane vibration of horizontally curved beams. Rao (1971) developed the governing differential equations of motion for the vibration of circular rings based on Hamilton’s principle [1]. Kirkhope (1976), Silva and Urgueira (1988) calculated the nature frequency of out-of-plane curved beams according to Dynamic Reciprocity Theorem [2, 3]. Wang T. M., (1980) also calculated the nature frequency for out-of-plane vibrations of continuous curves beams [4]. Kawakami M. (1995) has investigated both the in-plane and out-of-plane free vibrations of curved beams with variable sections [5]. Yang Y. B. (2001) derived an analytic solution for a horizontal curved simple supported beam subjected to vertical and horizontal moving loads [6]. However, the dynamic response of periodically supported curved Timoshenko beam subjected to moving loads has rarely been studied.

In comparison with the approaches mentioned above, the objective of this paper is to present a solution for the out-of-plane dynamic response of curved rail, which can be simplified as periodically supported curved Timoshenko beam under moving trains load.

In this paper, the general dynamic response of rail induced by moving load along curved path on an elastic semi-infinite space is obtained based on Duhamel Integral and Dynamic Reciprocity Theorem. Regarding periodic track structure, the general dynamic response equation in the frequency domain is simplified in a form of summation within the track sleeper period instead of integral. The transfer function of the curved beam is solved using the transfer matrix approach.

2. Dynamic response of curved beam subjected to moving loads (CBSML)

2.1. Moving load in the semi-infinite space

In the elastic semi-infinite space, as shown in Fig. 1, the vertical load g(t) moves along a curved path with the radius of R, the initial position of θ0, and the angular speed of c. The vertical dynamic displacement of receiver ξ at time-domain can be obtained based on Duhamel integral (Jia Y. X., 2009) [7]:

1
uξ,t=-+gτhzξ,θτ,t-τdτ,

where, u(ξ,t) is the vertical vibration displacement of receiver ξ, and right hand side of Eq. (1) represents the convolution integral of the time history of the moving load g(t) and the vertical transfer function hz(ξ,θ(τ),t-τ) between the time-dependent load position θ(τ) and receiver ξ. Note that, t-τ<0, hz(ξ,θ(τ),t-τ)0.

After transforming the time t to the circular frequency ω, using the Dynamic Reciprocal Theorem and the Forward Fourier Transform, the response displacement in the frequency domain can be expressed as:

2
u^ξ,ω=-+gτh^zθ0+cτ,ξ,ωexp-iωτdτ,

where, h^z(θ0+cτ,ξ,ω)=-+hz(θ0+cτ,ξ,t-τ)exp(-iω(t-τ))dt is the transfer function in the frequency domain. And the superscript “^” is used to indicate the expression in the frequency domain, similarly hereinafter.

Fig. 1Semi-infinite space subjected to moving load

Semi-infinite space subjected to moving load

Fig. 2Curved track subjected to moving load

Curved track subjected to moving load

2.2. Moving load on track structure

Since the curved track is symmetric, as shown in Fig. 2, only half of the curved track was considered during analysis. The curved track was subdivided into a number of track cells with the length of θcell, which is the space between sleepers. A vertical load g(t) traversed on the track structure, with the angle speed of c.

According to the relativity of motion, the load moved forward passing over one cell, equivalents to the load itself does not move, while the observation point moves in the opposite direction passing over a cell. Thus, the dynamic response in frequency domain can be simplified in a form of summation with the track sleeper spacing instead of integral, by converting the moving of the load on the rail to the moving of pick-up point by a specific sleeper spacing, which has been proved by the Floquet Transformation (Jia Y. X., 2009) [7].

At time stamp t, the load position in the global coordinate system can be expressed as: θ=θ0+ct, where, θ0 is the initial position of the load.

Fig. 3Local coordinate system

Local coordinate system

The local coordinate system θ~ is set up in track basic cell, as shown in Fig. 3. The relationship between the global coordinate system and the local coordinate system can be expressed as follows:

θ~=θ-nθθcell, θ~ξ=θξ-nξθcell, θ~0=θ0-n0θcell,

where, “~” is used to indicate the expression of the local coordinate system here and throughout this paper, nθ, nξ, n0 are the numbers of basic track cells θcell between the origin and the load position θ, between the origin and the pick-up point ξ, between the origin and the initial load position θ0 in the global coordinates, respectively. Thus, the following equation can be developed:

3
t=θ-θ0c=τ~+nθ-n0θcellc,
4
τ~=θ~-θ~0c.

According to Eq. (4), when the load moved by one sleeper spacing on the track, the vibration response at receiver ξ can be expressed as:

5
u^ξ,ωθcell=-+gτh^zθ0+cτ,ξ,ωexp-iωτdτ=0θcellcgτ~+nθ-n0θcellch^zθ0+cτ~+nθ-n0θcell,ξ,ω
*exp-iωτ~+nθ-n0θcellcdτ~.

In fact, nθ will change when a sleeper spacing was passed over. However, since the curve is not infinitely long and the angle of the curve is θt, nθ would change from n0 to n0+θt/θcell.

Appling Eq. (4), the expression of the time can be transformed to the expression of space, which can be expressed as:

6
u^ξ,ω=nθ=n0n0+θt/θcell1cθ~0θ~0+θcellgθ~-θ~0+nθ-n0θcellch^zθ~+nθθcell,ξ,ω
*xp-iωθ~-θ~0+nθ-n0θcellcdθ~.

Eq. (6) is the dynamic response of the track structure under vertical moving load in the frequency domain.

2.3. Transfer function of curved track

As proved in previous research (Jia Y. X., 2009) [7], the transfer function h^z(θ~,θ~ξ+(nθξ-nθ)θcell,ω) can be solved as the product of the state variables S(θ~,ω) at the load excitation point and the transfer function of the periodically supported beam, that can be divided into several basic track cells θcell. Besides that, the transfer function of basic track cell θcell can also be solved as the product of the transfer function of the curved beam and the support under the curved beam, using transfer matrix approach (Sun J. P., 2009) [8] as described below.

2.3.1. The transfer matrix of the curved beam

The curved track is simulated as periodically supported horizontal curved Timoshenko beam. The support under rail is modeled as mass-spring-damper element. For an infinitesimal element of curved beam as shown in Fig. 4, with the length measured along the neutral axis of the curved beam denoted as s, the x-, y- and z- axes are defined as tangential directions, radial and transverse direction, respectively and the origin of the coordinates moves along the neutral axis of the beam. In addition, u is transverse deflection, α is the slope due to pure bending, φ is the angle of torsion, R is the radius and θ is the central angle corresponds to the curve element. It is assumed the cross-section properties and material properties are constant along the beam. The shearing force Qz, bending moment My and the torsion moment Mx are also shown in the Fig. 4.

Fig. 4The coordinates of the curved beam element

The coordinates of the curved beam element

The shear deformation was considered for the analysis of an infinitesimal element ds in the curved beam, α can be expressed:

7
α=v+ux-φxR,

where, ν is the transverse shear angle.

And twisting angle γ can be calculated by:

8
γ=φx+1Rux.

Furthermore, the force-displacement relationship of the curve beam can be obtained from the following equations:

9
Qx=KGAv,
10
My=-EIyαx-φR=EIyφR-αx,
11
Mx=-EIs2γ2x+GIdφx+αR,
12
Bi=EIsγx,

where: E is the Young’s modulus, G is the shear modulus; K is the shear correction factor, Iy is the vertical bending moment of inertia, Id is free torsion moment of inertia, Is is polar moment of the cross-section and A is the sectional area.

Considering the homogeneous beam with infinite degrees of freedom, the dynamic equilibrium equations of the infinitesimal element of curved beam can be obtained, from its equilibrium condition listed below:

13
Qzx=ρA2ut2,
14
Myx=ρIy2αt2+Qx-MxR,
15
Mxx=ρIs2φt2+MyR,
16
Bix=-Mx+GIdγ,

where: Qz is the shearing force, My is bending moment, Mx is torsion moment, Bi is double warping moment, γ is warping angle and ρ is the mass per unit volume.

The state vector of any point in the curve beam can be expressed as:

17
S={Qz,My,Mx,Bi,u,α,φ,γ}T.

Thus, Eqs. (7)-(17) can be expressed using matrixes:

18
Sx=AS,

where:

A=0000-ρAω200010-1R00-ρIyω20001R0000-ρIsω2000-10000GId-1KGA00001xR00-1EIy00001R01RKGA0000-1R-xR210001EIs0000.

The general solution of Eq. (18) can be settled as:

19
Sx=eAxS0,

where, S0 is a constant matrix in the solution.

The curved beam can be divided into many infinitesimal elements, with the length of Δx, and then xk can be expressed as:

20
xk=kΔx, k=1, 2, 3,,
21
xk+1=xk+Δx.

Then:

22
Sxk+1=TiΔxSxk,

where, Ti(Δx)=eAΔx.

Based on the precise integration method of the exponential matrix (Sun J. P. 2009) [8]:

23
TiΔx=eAΔx=eAΔx2N2N=eAτ2N,

where, τ=Δx/2N, N= 20.

2.3.2. The transfer matrix of the support

For the periodically supported track structure, the periodic support is simulated as double-layer mass-spring-damper system, in which rail pad and sleeper pad are both modeled as spring-damper element, the sleeper is modeled as concentrate mass between the rail pad and sleeper pad. The double-layer support is calculated as a spring-damper element, as shown in Fig. 5, of which the composite stiffness kv can be expressed as:

24
kv=ckr(cksb-Msω2)ckr+(cksb-Msω2),

where, kr, ks, kb are the stiffness of rail pad, sleeper pad and subgrade, respectively. cr, cs, cb are the damping of rail pad, sleeper pad, and subgrade , respectively Ms is the sleeper mass.

Fig. 5The spring-damper element under the curved beam

The spring-damper element under the curved beam

Considering an infinitesimal element of the curved beam on the support as shown in Fig. 5, the state vectors of the two sides of the infinitesimal element are defined as follows.

The left side:

25a
SiL={QizL,MiyL,MixL,BiL,uiL,αiL,φiL,γiL}T.

While the right side:

25b
SiR={QizR,MiyR,MixR,BiR,uiR,αiR,φiR,γiR}T.

According to on the continuity requirements:

25c
uiR=uiL, αiR=αiL, φiR=φiL, γiR=γiL, QizR=QizL-ksleeperωui,
MiyR=MiyL, MizR=MizL, BiR=BiL,

where, ksleeper is the composite stiffness of the sleeper that was simplified as spring-damper element. Eq. (25) can be expressed as:

26
SiR=TsuppSiL,

where:

Tsupp=1000-ksleeperω00001000000001000000001000000001000000001000000001000000001.

2.3.3. Initial state vector of the curved beam under unit load

1) Unit load between two sleepers.

As shown in Fig. 6, the state vectors of the double sides of the curve beam element are defined as SL and SR for the left side and for the right side, respectively.

Based on transfer matrix, SR can be calculated as:

27
SR=TiΔxSL,
28
SR-SL=P.

Fig. 6Mechanical analysis of the beam element

Mechanical analysis of the beam element

Fig. 7Mechanical analysis of the beam element with support

Mechanical analysis of the beam element  with support

2) Unit load on the sleeper:

29
-QzL+QzR+1=Psleeper=Xksleeperω,
30
SR=TiΔxSL,

where, Ti(Δx)=Ti(Δx/2)LTsuppTi(Δx/2)R.

Thus the state vector can be settled, as shown in Fig. 7. With the initial state variables and the transfer function of the curved beam settled, the dynamic response of the periodically supported curved track structure under the moving load can be solved.

3. Model validation

Based on the mathematical model described above, the calculation program is developed. In order to validate this model, the rail dynamic responses were calculated and compared with a special case from the literature review. Furthermore, assembled with a well proved vehicle/track interaction model, the rail dynamic response was calculated and compared with the measurement data collected from Beijing metro.

3.1. Model validation based on previous research

Yang Y. B. (2001) [6] proposed an analytical model to calculate the vibration of simple supported curved beam subjected to moving load, as shown in Fig. 8. In order to validate the model in this paper, the example was recalculated with the same input data as below.

Fig. 8Simple supported curved beam

Simple supported curved beam

Where, the properties of beam cross section, a=5 m,b=1.8 m,A=ab=9 m2,Ix=ab3/12=2.43 m4,Iy=ba3/12=18.75 m4,Jθ=Ix+Iy=21.18 m4, the length of beam L=24 m, the subtended angle α-=30°=π/6, the radius R=45.84 m, Young’s modulus E=32.3×109 N/m2, the Poisson ratio v=0.2, shear modulus G=E/21+v, load P=9.8×29.9×103 N, moving speed Vp=40 m/s, and damping ζd=0.

As shown in Fig. 9, the mid-span vibration displacement of the simple supported curved beam under moving load is obtained based on the solution developed above.

Fig. 9Displacement response of the simple supported beam

Displacement response of the simple supported beam

Fig. 10Calculated results by Yang Y. B. (2001)

Calculated results by Yang Y. B. (2001)

As shown in Fig. 9 and Fig. 10, it is proved that the calculated displacement response has good agreement with the example given by Yang Y. B. (2001) [6], which confirms the reliability of the presented solution.

3.2. Model validation based on measurement data

3.2.1. Assembling the vehicle and CBSML model

In order to further verify the application of the CBSML model, a calibrated vehicle model [9-11] (Fig. 11) was employed into the assembled model to simulate the vertical moving loads on the curved rail. To assemble the vehicle model with CBSML model, Hertz wheel/rail contact model (Fig. 12) is used and the array of the moving load on the rail beam is shown in Fig. 13.

The calculation input data of vehicle model are listed in Table 1, while in the track structure model, kr= 4.2 MN/m, cr= 0.05 MN·s/m, train speed V= 50 km/h, and radii of the curve R= 400 m.

Fig. 11Vehicle model

Vehicle model

Fig. 12Wheel/Rail Hertz contact model

Wheel/Rail Hertz contact model

Fig. 13Array of train loads on the rail beam

Array of train loads on the rail beam

Table 1Parameters of vehicle model

Items
Symbol
Unit
Value
Mass of carriage
Mc
t
35.0
Carriage mass moment of inertia
Jc
tm2
1700
Length of carriage
L
m
19
Mass of bogie
Mt
t
4.60
Bogie mass moment of inertia
Jt
tm2
9.62
Bogie stiffness
Kt
kN/m
2080
Bogie damping
Ct
kN∙s/m
240
Space between wheels
l1
M
2.20
Space between bogies
l2
M
12.6
Wheel mass
Mw
t
1.42
Wheel stiffness
Kw
kN/m
2450
Wheel damping
Cw
kN∙s/m
24

3.2.2. Measurement setup

A campaign of measurement was conducted to estimate the vibration isolation performance of different track forms at the beginning of opening in certain lines in Beijing metro. Therefore, the pass-by rail dynamic responses were collected when both the vehicle and track structure were in good condition. A sharp curve (R= 400 m, V= 50 km/h) was selected as the control measurement section. The fastening systems were directly fixed on the concrete slab without sleeper. The design static stiffness of the fastening system is 4-8 MN/m, but the average static stiffness of 3 systems selected in random from this section is 4.2 MN/mm according to the laboratory test. The detailed measurement condition is listed in Table 2.

Table 2Measurement condition of pass-by rail dynamic response

Item
Value
Radii of curve
400 m
Static stiffness of fastener
4-8 MN/m
Space of fastener
625 mm
Track gauge
1435 mm
Rail
CHN 60 (GB 2585-2007): 60 kg/m
Car marshalling
6 cars (T-M-T-M-M-T)
Axle load
140 kN
Wheel load
70 kN
Bogies space
2.3 m
Train speed at test site
50 km/h

The accelerometer was mounted on the rail with a magnetic support in both vertical and lateral direction, as shown in Fig. 14. However, in this paper, only the vertical rail dynamic response was considered to demonstrate the calculation model.

Fig. 14Accelerometer mounting

Accelerometer mounting

The experimental configuration for the vibration measurements consists of Lance-LC0123T piezoelectric accelerometers and INV3018C dynamic data acquisition system. The accelerometer has the sampling frequency between 0.2 and 11000 Hz, the sensitivity of 26.0 mV/g and the maximum measuring acceleration of 200 g. During the field measurement, the sampling frequency was setup to 5120 Hz, while the sampling distance is 0.625 Hz. The post-processing of the measured data in this paper was carried out in the software DASP, and Hanning window was taken into consideration.

3.2.3. Result comparison

Then the calculated dynamic response results are compared with the measurement data collected from Beijing metro (see Fig. 15) [12].

As shown in Fig. 15, it is observed that there are six wave vibration signals from the calculated time history, each of which corresponds to one car that passed-by. But in raw measurement data, it is hard to clearly distinguish the dynamic response of every car from the irregular signal. It could result from the different of actual axle loads, since the passengers distribute in different cars. However, certain degree of agreement was still observed between raw testing data and analytical results in terms of occurring time and magnitude of the vibration.

Thus, in order to better interpret the testing data and further compare the testing and analytical results, the data was transformed to frequency domain since the frequency spectrum is more interpretable in the analysis of train/track interaction and the estimation of the level of vibration energy. In this study, the frequency spectrum is presented in the form of 1/3 octave band.

Fig. 15Vibration acceleration of curved track under moving train (V= 50 km/h, R= 400 m)

Vibration acceleration of curved track under moving train (V= 50 km/h, R= 400 m)

a) Time history (measurements)

Vibration acceleration of curved track under moving train (V= 50 km/h, R= 400 m)

b) Time history (calculation)

Vibration acceleration of curved track under moving train (V= 50 km/h, R= 400 m)

c) Frequency spectrum in 1/3 octave band

From Fig. 15(c), it is clear that the calculated result shows the same tendency with measurement result, especially in the range of 0-25 Hz. As is hardly decayed along the transit path, the low frequency (0-25 Hz) train-induced vibration would cause ground vibration, which generates perceivable vibration and re-radiated noise that would discomfort the occupants in the nearby building. It is indicated that the analytical model could be used to predict the ground vibration problem induced by rail transit systems.

Since the irregularities of track structure and roughness of wheel/rail contact surface were not considered in calculation model, differences were found in the range of 31.5-200 Hz. However, the tendency of both curves is similar and the accuracy of the model would be improved if the roughness of the wheel/rail contact surface is added to the model.

4. Conclusions

An analytical solution was derived for the out-of-plane vibration response of periodically supported curved Timoshenko beam subjected to moving loads. This study is aiming at providing a solution for rail dynamic response of curved rail thereby establish the relationship between the dynamic behaviour and rail corrugation.

The analytical model for vibration of simple supported curved beam under moving load was validated based on results from previous research. Furthermore, in order to further validate the presented solution, a train/track dynamic model was introduced into the assembled model and the calculated results were compared with the measurement data collected from Beijing Metro. The comparison proved that the presented solution could provide accurate results and can be used in future studies.

The advantage of the present approach is that it provides clear mechanism insights into the various vibration phenomena induced by vehicles, in particular, the phenomena of resonance and cancellation, and allows us to identify the key parameters involved.

References

  • Rao S. S. Effects of transverse shear and rotatory inertia on the coupled twist-bending vibrations of circular rings. Journal of Sound and Vibration, Vol. 16, Issue 4, 1971, p. 551-566.
  • Kirkhope J. Out-of-plane vibration of thick circular ring. Journal of the Engineering Mechanics Division, Vol. 102, Issue EM2, 1976, p. 239-247.
  • Silva J. M. M., Urgueira A. P. V. Out-of-plane dynamic response of curved beams – an analytical model. International Journal of Solids and Structures, Vol. 14, Issue 3, 1988, p. 271-284.
  • Wang T. M., Nettleton R. H., Keita B. Natural frequencies for out-of-plane vibrations of continuous curved beams. Journal of Sound and Vibration, Vol. 68, Issue 3, 1980, p. 427-436.
  • Kawakami M., Sakiyama T., Matsuda H., Morita C. In-plane and out-of-plane free vibrations of curved beams with variable sections. Journal of Sound and Vibration, Vol. 187, Issue 3, 1995, p. 381-401.
  • Yang Y. B., Wu C. M. Dynamic response of a horizontally curved beam subjected to vertical and horizontal moving loads. Journal of Sound and Vibration, Vol. 242, Issue 3, 2001, p. 519-537.
  • Jia Y. X. Study on analytical model of coupled vehicle & track and effect to environment by metro train-induced vibrations. Department of Urban Rail Transit, Beijing Jiao tong University, 2009, (in Chinese).
  • Sun J. P., Li Q. N. Precise transfer matrix method for solving earthquake response of curved Box Bridge. Journal of Earthquake Engineering and Engineering Vibration, Vol. 29, Issue 4, 2009, p. 139-146.
  • Yinxuan Jia, Weining Liu, et al. Influence of vibration on the existing metro structure induced by trains operated on Beijing underground connecting line. China Railway Science, Vol. 29, Issue 5, 2008, p. 72-77, (in Chinese).
  • Yinxuan Jia, Weining Liu, et al. Vibration effect on surroundings induced by passing trains in spatial overlapping tunnels. China Railway Science, Vol. 32, Issue 2, 2009, p. 104-109, (in Chinese).
  • Hougui Zhang Study on effect on nearby metro structures due to train induced vibrations on underground diameter Line in Beijing. The Thesis in Candidacy for the Master Degree of Bridge and Tunneling Engineering, 2007, (in Chinese).
  • Kefei Li, Weining Liu, et al. In-situ test and analysis on the vibration mitigation measures of the elevated line in Beijing metro line 5. China Railway Science, Vol. 30, Issue 4, 2009, p. 25-29, (in Chinese).

About this article

Received
11 March 2014
Accepted
20 April 2014
Published
30 June 2014
Keywords
analytical solution
dynamic response
out-of-plane curved beam
moving loads
Acknowledgements

The presented research in this paper was sponsored by National Science Foundation of China (No. 51378001), the Fundamental Research Funds for the Central Universities (No. 20110009120023) and the Fundamental Research Funds for the Central Universities (No. 2013JBM061).