Second Order Differential Equation of Motion in Railways: The Variance of The Dynamic Component of Actions Due To The Sprung Masses of The Vehicles
Second Order Differential Equation of Motion in Railways: The Variance of The Dynamic Component of Actions Due To The Sprung Masses of The Vehicles
Second Order Differential Equation of Motion in Railways: The Variance of The Dynamic Component of Actions Due To The Sprung Masses of The Vehicles
net/publication/292988508
CITATIONS READS
2 126
1 author:
Konstantinos Giannakos
freelancer
57 PUBLICATIONS 165 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Minoan/Mycenaean Ashlar Technology, Mycenaeans’ Overseas Expansion, Attarissiya’s Raids and LCIIA-LCIIB Cyprus View project
All content following this page was uploaded by Konstantinos Giannakos on 05 February 2016.
http://www.iaras.org/iaras/journals/ijtam
Abstract: - In this paper the second order differential equation of motion is presented for the case of a railway
vehicle rolling on a railway track with defects/faults and its solution is presented for the Sprung (Suspended)
Masses of the vehicle that act indirectly on the track through the springs and dampers of the vehicle.
Key-Words: - Static/Dynamic Stiffness Coefficient, Sprung Masses, Unprung Masses, Fourier Transform,
Spectral Density, Variance, Standard Deviation, Dynamic Component of Actions.
30
simplified model in Fig. 2. An analysis of the Masses (axle, wheels and a fraction of the semi-
behavior of the Non Suspended Masses was suspended motive electromotor in the locomotives)
presented in [1], [2]. For the simulation of the which act without any damping directly on the track
Sprung (Suspended) Masses of the vehicles (that is panel. On the contrary the Sprung (Suspended)
of the masses located over the primary suspension Masses that are cited above the primary suspension
of the vehicle), the forces resulting from the forcing of the vehicle, act through a combination of springs
excitation imposed by the rail running table –in the and dampers on the track. A part of the track mass is
most of the cases– have small effect on the rolling also added to the unsprung masses, which
of the wheel, but for defects of very long participates in their motion ([3], [4]).
wavelength, the oscillations of the Suspended The rail running table has the shape of a wave
Masses become predominant since the oscillations that is not completely “rectilinear”, that is, it does
of the Non Suspended Masses decrease. The not form a perfectly straight line but contains faults/
Suspended/Sprung Masses are examined in the defects, varying from a few fractions of a millimeter
present paper. to a few millimeters, and imposes forced oscillation
on the railway vehicles that move on it; the
faults/defects are represented by the ordinate n in
Fig. 2. Moreover during the rolling of the wheel, a
2 Problem Formulation deflection y of the rail running table appears (see
[5]), since the support (track) is not undeflected.
The general form of the equation of a vehicle
moving on an infinitely long beam, on elastic
ground without damping, is described in the next
paragraph 2.2..
31
d2 y dy algebraically), we can approach the problem of the
( mNSM + mTRACK ) ⋅ 2
+ Γ ⋅ + hTRACK ⋅ y = random excitation, from cosine defect (V<<
dt dt
d 2n Vcritical=500 km/h).
= −mNSM ⋅ + ( mNSM + mSM ) ⋅ g (1)
dt 2
2.3 Fourier Transform and the Second
where: mNSM the Non-Suspended/Unsprung Order Differential Equation of Motion
Masses (NSM) of the vehicle in tonnes-mass, We begin from the hypothesis of a cosine form
mTRACK the mass of the track that participates in the defect on the rail running table of the form:
motion of the NSM (for its calculation see Ref. [3]
V ⋅t
and [4]), mSM the Suspended/Sprung Masses (SM) η = a ⋅ cos ω t = a ⋅ cos 2π ⋅ (4)
of the vehicle that are cited above the primary λ
suspension of the vehicle, Γ damping constant of the where: η the ordinate of the defect along the
track, hTRACK the total dynamic stiffness coefficient track (abscissa x), V the speed of the vehicle, t the
of the track (its calculation Eqn. 3 below), n the time and λ the wavelength of the defect, so:
fault ordinate of the rail running table, g the
2π 2π 2π Vt
acceleration of gravity and y the total deflection of T= ⇒ ωt = t ⇒ ωt = (5a)
the track. Furthermore: ω T λ
1 1 1 1 1 1 since the wheel overpasses the wavelength λ of
= + + + + (2) the defect, in:
ρtotal ρ rail ρ pad ρ sleeper ρballast ρ subgrade
λ
T= ⇒ λ = T ⋅V (5b)
where: ρi the static stiffness coefficients of the V
constitutive layers of the track, the quasi spring If we set:
constants of the layers and ρtotal the resultant total
static coefficient of the track and htrack the total mSM + mNSM dy dz
dynamic stiffness coefficient of the track given by: y=z+ ⋅g ⇒ =
hTRACK dt dt
1 ρtotal
hTRACK = ρ = ⋅ 4 E⋅ J ⋅ (3) its second derivative will be:
dyn−total
2⋅ 2
d2y d2z
with E, J the modulus of elasticity and the =
moment of inertia of the rail (steel) and ℓ the dt 2 dt 2
distance among the sleepers. mNSM + mSM
The phenomena of the wheel-rail contact and of where the quantity ⋅ g represents
hTRACK
the wheel hunting, particularly the equivalent the subsidence due to the static loads only, and z
conicity of the wheel and the forces of pseudo-glide,
are non-linear. In any case the use of the linear random (see [20]) due to the dynamic loads. Eqn (1)
system’s approach is valid for speeds lower than the
Vcritical≈500 km/h. The integration for the non-linear becomes:
model (wheel-rail contact, wheel-hunting and
pseudoglide forces) is performed through the Runge d2z dz d 2n
mNSM + Γ ⋅ + hTRACK ⋅ z = − m NSM ⋅ ⇒ (6a)
Kutta method ([7], p.94-95, 80, [8], p.98, see also dt 2 dt dt 2
[9], p.171, 351).
In Fig. 2 the rail running table depicts a
d 2 z d 2n dz
longitudinal fault/ defect of the rail surface. In the ⇒ mNSM 2 + 2 + Γ ⋅ + hTRACK ⋅ z = 0 (6b)
above equation, the oscillation of the axle is damped dt dt dt
after its passage over the defect. Viscous damping,
due to the ballast, enters the above equation under Since, in this case, we are examining the
the condition that it is proportional to the variation dynamic loads only (derived from the actions of the
of the deflection dy/dt. To simplify the Suspended and Non Suspended Masses) , in order to
investigation, if we ignore the track mass (for its approach their effect, we could narrow the study of
calculation Ref. [3] and [4]) in relation to the much equation (6b), by changing the variable:
larger Vehicle’s Non Suspended Mass and bearing d 2u d 2 n d 2 z
in mind that y+n is the total subsidence of the wheel u =n+ z⇒ = +
during its motion (since the y and n are added dt 2 dt 2 dt 2
32
Equation (6) becomes: where:
d 2u dz Γ 2 ⋅ ω 2 + hTRACK
2
mNSM + Γ ⋅ + hTRACK ⋅ z = 0 ⇒ (7a) B (ω ) =
2
(10b)
(m ⋅ ω 2 − h ) + Γ2 ⋅ ω 2
2
dt dt 2
NSM
d 2u d (u − n)
⇒ mNSM 2
+Γ⋅ + hTRACK ⋅ ( u − n ) = 0 (7b) B(ω) is a complex transfer function, the
dt dt
frequency response function, that makes it possible
to pass from the fault n to the u=n+Z. Practically it
where, u is the trajectory of the wheel over the
is verified also by the equation:
vertical fault (of ordinate n) in the longitudinal
B ( ω ) = H ( ω ) ⋅ G (ω ) =
2 2 2
profile of the rail.
If we apply the Fourier transform to the equation 2
hTRACK + Γ2 ⋅ ω 2
(6a) (see relevantly Ref. [10] for solving second = (10c)
(m ⋅ ω 2 − hTRACK ) + Γ 2 ⋅ ω 2
2
order differential equations with the Fourier NSM
transform):
passing from n to Z through H(ω) and afterwards
Γ ⋅ ( iω ) h from Z to n+Z through G(ω). This is a formula that
( iω ) ⋅ Z (ω ) + ⋅ Z (ω ) + TRACK ⋅ Z ( ω ) =
2
U (ω ) + Γ2 ⋅ ω 2 1 + 4ζ 2 ⋅ β 2
h2 B ( ω ) = Bn ( β ) =
2 2
G (ω ) = , G (ω )
2
= TRACK2 (10d)
(1 − β 2 ) + 4ζ 2 ⋅ β 2
(9)
Z (ω )
2
mNSM ⋅ ω 4
G(ω) is a complex transfer function, the The transfer function C(ω) of the second
frequency response function, that makes it possible
d 2 ( Z + n)
to pass from Z to Z+n. derivative of (Z+n) in relation to time: ,
If we name U the Fourier transform of u, N the dt 2
Fourier transform of n, p=2πiν=iω the variable of –that is the acceleration γ–, will be equal to ω∙Β(ω):
frequency and Δ̂Q the Fourier transform of ΔQ and
apply the Fourier transform at equation (7b):
p 2 ⋅ U ( ω ) = p 2 B (ω ) ⋅ p 2 ⋅ N (ω ) ⇒
d 2u du dn
Eq.(7) ⇒ mNSM +Γ⋅ + hTRACK ⋅ u = Γ ⋅ + hTRACK ⋅ n ⇒ (11a)
dt 2 dt dt
⇒ U (ω ) = p B ( ω ) ⋅ N ( ω )
2
(m NSM ⋅ p 2 + Γ ⋅ p + hTRACK ) ⋅ U = ( Γ ⋅ p + hTRACK ) ⋅ N ⇒ C (ω )
that is:
Γ ⋅ p + hTRACK
U (ω ) = ⋅ N (ω ) (10a) ω2
mNSM ⋅ p 2 + Γ ⋅ p + hTRACK C ( ω ) = ω 2 ⋅ B (ω ) = ωn2 ⋅ ⋅ B (ω ) ⇒
ωn2
B (ω )
⇒ C ( ω ) = ωn2 ⋅ β 2 ⋅ B (ω ) (11b)
33
The increase of the vertical load on the track due where λ is the wavelength of the defect. This
to the Non Suspended Masses, according to the means that circular frequency in space Ω is the wave
principle force = mass x acceleration, is given by: number k of the equation of oscillation, and [13]:
d 2u d 2 (n + Z ) ∞ ∞
∆Q = mNSM ⋅ 2 = mNSM ⋅
dt dt 2
(12) ∫ S ( Ω ) ⋅ d Ω = ∫ s (ω ) ⋅ dω ⇒
0 0
π 0
through the vehicle an output signal: the wheel
trajectory. The output spectral density and the input where n is the random variable of the defect
spectral density of the excitation are related through (input), z the subsidence of the wheel (output) and
equation [14], [5]: ΔQ the dynamic component of the Load that is
added to the Static Load of the wheel due to the Non
SOUTPUT (ω ) = H ( iω ) ⋅ S INPUT (ω )
2 Suspended Masses (output also).
(15a)
From these equations and the analytic form of
In order to relate the temporal spectrum with the the spectrum of the defects/faults, we can calculate
spectrum in space we use the following equation: the mean square value of the dynamic component of
the Load due to the Non Suspended Masses that is
2π Vt 2π
ω ⋅t = ⇒ω = ⋅V ⇒ ω = Ω ⋅V (15b) added to the relevant dynamic component of the
λ λ
34
+∞
A
Suspended Masses and the total dynamic component +∞
A 1 1
of the load is added to the Static Load of the wheel. ⇒ σ ( z) = ⋅
2
∫ ( x) ⋅ dx = − ⇒
π 0
3
2π ( B + Ω )2
0
+∞
A 1
σ ( z) = −
2
⋅ 2 2
⇒
2π B + 2 BΩ + Ω 0
A 1
Fig. 3 Power Spectral density S(Ω) –the black ⇒ σ 2 ( z) = − 0− 2 ⇒
2π B
curve– Variance (mean square value x̅2) –the shaded
area– of a function [5]. 1 A
σ 2 (z) = ⋅ (23)
2π B 2
From the power spectral density and the variance
functions and their definitions [5]:
If we examine only the much more severe case,
S ∆Q ( ω ) = S n ( ω ) ⋅ B ( ω )
2
(19a) for the case of the Non Suspended Masses, of the
defects of short wavelength, consequently large Ω
–like the undulatory wear– then we can omit the
∆Q = mNSM ⋅ ∆γ ⇒ σ 2 ( ∆Q ) = mNSM ⋅ σ 2 ( γ(19b)
)⇒ term B, and using Eqn. (15b):
σ 2 ( ∆Q ) A A A ⋅V 3
⇒ σ (γ ) =
2
(19c) Sn ( Ω ) = = = (24)
mNSM Ω3 1
⋅ ω 3 ω3
V3
and using the eqs. (19) and (11-12b): The term B characterizes the defects with large
wavelengths, for which the maintenance of track is
1 effective, and when we examine this kind of defects
σ 2 (γ ) = ⋅ σ 2 ( ∆Q ) ⇒
mNSM term B should be taken into account. Suspended
+∞ Masses should be examined for long wavelength
1
σ 2 (γ ) = ∫ B (ω ) ⋅ S n (ω ) ⋅ d ω (20)
2
⋅ defects.
m 2
NSM ⋅π 0 For the line “Les Aubrais – Vierzon”, the
2 +∞ parameters values are: B=0,36, A=2,1∙10-6 and S(Ω)
mNSM
σ 2 (γ ) = ∫ β 4 ⋅ ωn4 ⋅ B (ω ) ⋅ S n (ω ) ⋅ d ω =
2
⋅ is calculated in m3 and σ(z)=1,57 mm. The
mNSM ⋅ π 0
2
eigenfrequency of the Non Suspended Masses of the
σ 2 (γ ) = vehicles is approximately 30 – 40 Hz and even for
speeds of 300 km/h there are wavelengths less than
+∞
1 1 + 4ζ 2 ⋅ β 2
= ⋅ ∫β
4
⋅ ωn4 ⋅ ⋅ S n ( ω ) ⋅ d(21)
ω 3 m [17].
π (1 − β 2 ) + 4ζ 2 ⋅ β 2
2
0
From equations (16) and (24):
1 1 A ⋅V 3 A ⋅V 2
Sn (ω ) = ⋅ S (Ω) = ⋅ = 3 3 (25)
From the above equations and the analytical form V V ω 3
ωn ⋅ β
of the spectrum of the longitudinal defects/ faults of
the track we could effectively calculate the variance
(mean square value) of the dynamic component of
the Loads on the track panel due to the Non 3 The Case of Non Suspended Masses
Suspended Masses. All the results of measurements
on track in the French railways network show that This case has been analyzed in [1], [2], [5], [18].
the spectrum of defects in the longitudinal level has The interested reader should read the relevant texts.
the form [6], [17]:
A
Sn ( Ω ) = (22)
( B + Ω) 4 The Variance of the Suspended
3
35
then the conclusion that will be derived can be used AV 2 1 + 4ζ 2 ⋅ β 2
sγ (ω ) = ωn4 ⋅ ⋅ ⋅
(1 − β )
for more complicated cases of rolling of vehicles,
( BV + ω )
3 2 2
motion of car-bodies etc. Furthermore we consider + 4ζ 2 ⋅ β 2
the simplified model of Fig. 1 with one-floor mass- 1 + 4ζ ′2 ⋅ β ′2
spring-damper system rolling on a rail’s surface. ⋅ (27)
In order to calculate the power spectrum density (1 − β ′ )
2 2
+ 4ζ ′2 ⋅ β ′2
of the excitation sE(ω) from the excitation spectrum
of the wheel sυ(ω), we apply the Eqn. 15a with the The variance of the accelerations of the car-body
Eqns. 16, Eqn. 22 (the parameter B is not omitted) of the railway vehicles is given by (Eqns. 18):
and Eqns. 10 (in order to pass from the defect n to
n+Z): ∞
1
σ (γ ) = ⋅ sγ ( ω ) ⋅ d ω
π ∫0
2
(28)
sE (ω ) = B (ω ) ⋅ sυ ( ω ) ⇒
2
1 − 2
+ 4ζ ⋅
ωn ωn2 ∆QSM = mSM ⋅ ∆γ ⇒ σ 2 ( ∆QSM ) = mSM ⋅ σ 2 ( γ ) (29)
where: sE(ω) is the power spectrum density of the Finally an approximation could be used for the
excitation, ωn is always the eigenfrequency of the calculation of the variance of this part of the
Non Suspended Masses, ζ the damping coefficient dynamic component of the load (see [5], [6]):
of the track, sυ(ω) the spectrum of the excitation of
the wheel due to the track defects/faults and V − 40
│B(ω)│the modulus of the transfer function of the σ ( ∆QSM ) = ⋅ N L ⋅ Qwheel (30)
1000
motion of the wheel.
From the Eqns. 10 and the Eqns. 11, 12 and 13, where: Qwheel is the static wheel load, V is the
operational speed, and the coefficient NL is the
with the analysis cited above, we keep that, C(ω) ις
mean standard deviation of the longitudinal level
the transfer function of the second derivative of condition of the track, on a 300 m length
d 2 ( Z + n) approximately, for both rails is the mean standard
(Z+n) in relation to time: , that is the
dt 2 deviation of the longitudinal level condition of the
acceleration γ and it is equal to ω∙B(ω). track, on a 300 m length fluctuating between 0,7–
1,5 mm or more (see [6]; [19], p. 335–336); for the
In the case of the Suspended Masses │C(ω)│is
Greek network NL is estimated to fluctuate –
the modulus of the transfer function of the
mainly– between 1 and 1,5 [5].
accelerations of the car-body. Consequently for the
In more details NL, the average of the brutal
spectrum of the accelerations of the car-body we will
signal on a basis of approximately 300 m for the
use Eqn. 15a substituting the parameters ζ and ωn of
vertical and horizontal defects of the two rails, is the
the track with the relevant parameters ζ′ (damping
convolution:
coefficient) and ω′n (eigenfrequency) of the car-
body. +∞ x − x0
For the railway vehicles the eigenfrequencies ω′n N L ( x0 ) =
1
η l ( x )e
∫
300
dx (31)
of the car-body are in the area of 1 Hz, since with the 300 −∞
development of high-speeds it could arrive 10 Hz.
For the damping coefficient of the car-body of the
where ηl(x) is the value of the primary signal; in
railway vehicles two characteristic values of ζ′ could
practice a weighted average index which “crashes”
be used with reliability: 0,15 and 0,20 (see relevantly
less the isolated defects than one classic average and
[17] and [7]).
simulates roughly the “memory” of the vehicle (see
relevantly [19]).
sγ (ω ) = C (ω ) ⋅ s E (ω ) ⇒
2
36
5 Conclusions [5] Giannakos K. Actions on the Railway Track,
Papazissis publ., www.papazisi.gr, Athens,
In this paper the second order differential equation 2004.
of motion of the Sprung (Suspended) Masses of a [6] SNCF/Direction de l’ Equipement, Mecanique
railway vehicle rolling on the track with de la Voie, 1981.
defects/faults is investigated and its solution is [7] Alias J. La Voie Ferree – Techniques de
approached. The investigation is performed through Construction et Entretien, deuxieme edition,
the Fourier transform and the solution is verified Eyrolles, Paris, 1984.
from findings from a research program performed [8] Fortin J. La Deformee Dynamique de la Voie
by the Greek railways in collaboration with the Ferree, Revue Generale des Chemins de Fer
French state railways (SNCF) and further research (RGCF), 02/1982.
performed by the author. [9] Thompson D. Railway Noise and Vibration,
Elsevier, 2009.
[10] Wylie C.R. and Barrett L.C. Advanced
References: Engineering Mathematics, sixth edition,
[1] Giannakos K. Modeling the Influence of Short McGraw-Hill, Inc., USA, 1995.
Wavelength Defects in a Railway Track on the [11] Gent I. and Janin G. La Qualite de la Voie
Dynamic Behavior of the Non-Suspended Ferree, SNCF’s reprint.
Masses, Journal of Mechanical Systems and [12] Papoulis A. Signal Analysis, McGraw-Hill Inc.,
Signal Processing (Mech. Syst. Signal USA, 1977.
Process.), Elsevier, 2015, Vol. 68-69, 2016, pp. [13] Roddier F. Distributions et Transformation de
68-83. Fourier, Edisciense, Paris, 1971.
[2] Giannakos K. Track Defects and the Dynamic [14] Clough R.W. and Penzien J. Dynamics of
Loads due to Non-Suspended Masses of Structures, second edition, McGraw-Hill Intl.,
Railway Vehicles, NAUN International Journal Singapore, 1993.
of Mechanics, Vol.7, Issue 3, 2013, pp. 180- [15] Papoulis A. Probability, Random Variables
191. and Stochastic Processes, third edition,
[3] Giannakos K. Theoretical calculation of the McGraw-Hill Intl., Singapore, 1991.
track-mass in the motion of unsprung masses in [16] Spiegel M.R. Probability and Statistics,
relation to track dynamic stiffness and Schaum’s Outline Series, McGraw-Hill, New
damping, International Journal of Pavement York, greek edition, Athens, 1977.
Engineering (IJPE), Special Rail Issue “High- [17] Prud’Homme A. La Voie, Revue Generale des
Speed Railway Infrastructure: Recent Chemins de Fer (RGCF), Janvier, 1970, extrait
Developments and Performance”, Vol. 11, de RGCF.
Issue 4, August 2010, pp. 319-330. [18] Giannakos K. and Loizos A. Evaluation of
[4] Giannakos K. Influence of the track's damping actions on concrete sleepers as design loads –
on the track mass participating in the motion of Influence of fastenings, International Journal
the Non Suspended Masses of railway vehicles of Pavement Engineering (IJPE), Vol. 11, Issue
– theoretical calculation and comparison to 3, June 2010, pp. 197 – 213.
measurements, volume published in honour of [19] Janin G. La Maintenance de la Geometrie de la
(fs) professor George Giannopoulos, Aristotle Voie, Revue Generale des Chemins de Fer
University of Thessaloniki, 2012. (RGCF), Juin, 1982, pp. 331–346.
37