27.04.2014 Views

Simulation Modelling of the Railway Vehicle Guidance Mechanism

Simulation Modelling of the Railway Vehicle Guidance Mechanism

Simulation Modelling of the Railway Vehicle Guidance Mechanism

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

SYSTEMS DIVISION RESEARCH 2013 1<br />

<strong>Simulation</strong> <strong>Modelling</strong> <strong>of</strong> <strong>the</strong> <strong>Railway</strong> <strong>Vehicle</strong><br />

<strong>Guidance</strong> <strong>Mechanism</strong><br />

Priya Parthasarathy, Dr. Christopher Ward, Dr. Roger Dixon<br />

Abstract—This paper presents <strong>the</strong> initial stage <strong>of</strong> work<br />

carried out on <strong>the</strong> simulation <strong>of</strong> a simplified rail vehicle model.<br />

The stability analysis <strong>of</strong> <strong>the</strong> wheelset with a change in <strong>the</strong><br />

lateral track position has also been discussed. It puts forward<br />

<strong>the</strong> significance <strong>of</strong> yaw suspension stiffness on <strong>the</strong> dynamic<br />

stability <strong>of</strong> <strong>the</strong> wheelset. The results <strong>of</strong> <strong>the</strong> simplified rail<br />

vehicle modelling and a brief idea about <strong>the</strong> future work have<br />

been analysed.<br />

I. INTRODUCTION<br />

The initial phase <strong>of</strong> <strong>the</strong> work was to develop a simple rail<br />

vehicle model to understand <strong>the</strong> dynamics <strong>of</strong> a rail vehicle<br />

system. This involved simulating a linearised vehicle model<br />

with a single wheelset and a suspended mass. The wheelset<br />

is modelled in <strong>the</strong> lateral and yaw direction and <strong>the</strong><br />

suspended mass in <strong>the</strong> lateral direction. The lateral position<br />

<strong>of</strong> <strong>the</strong> wheelset and <strong>the</strong> mass suspension and <strong>the</strong> yaw angle<br />

<strong>of</strong> <strong>the</strong> wheelset have been determined for a 10mm step<br />

change in <strong>the</strong> lateral position <strong>of</strong> <strong>the</strong> track.<br />

II. BACKGROUND<br />

The low adhesion at <strong>the</strong> wheel rail contact causing huge<br />

problems to <strong>the</strong> railway network is one important aspect to<br />

be outlined. This project puts forward <strong>the</strong> real time<br />

estimation <strong>of</strong> low adhesion conditions with <strong>the</strong> use <strong>of</strong><br />

advanced monitoring techniques.<br />

The first step <strong>of</strong> this work is to develop simulation<br />

models <strong>of</strong> a simple rail vehicle and increasing <strong>the</strong><br />

complexity levels from simple linear models to nonlinear<br />

vehicle models as <strong>the</strong> project progresses.<br />

III. INITIAL STAGE OF WORK<br />

A linearised rail vehicle model comprising <strong>of</strong> a single<br />

wheelset and a suspended mass as shown in Fig 1[1] have<br />

been simulated to test <strong>the</strong> stability <strong>of</strong> <strong>the</strong> wheelset in<br />

response to <strong>the</strong> track displacement. The wheelset equations<br />

are a combination <strong>of</strong> suspension forces, creep forces<br />

generated at <strong>the</strong> wheel-rail contact and gravitational stiffness<br />

generated by <strong>the</strong> geometry <strong>of</strong> <strong>the</strong> wheel and rail.<br />

Fig 1 shows <strong>the</strong> simplified rail vehicle model with a single<br />

wheelset and a suspended mass. The suspension forces and<br />

<strong>the</strong> yaw moments, <strong>the</strong> lateral displacement <strong>of</strong> <strong>the</strong> wheelset<br />

and <strong>the</strong> mass are also shown in <strong>the</strong> figure below.<br />

It can be seen that <strong>the</strong>re is a 10mm step change in <strong>the</strong> track<br />

position for which <strong>the</strong> lateral positions and <strong>the</strong> yaw angle<br />

are observed.<br />

10mm<br />

Wheelset<br />

fΨ <br />

Ym<br />

kΨ <br />

y<br />

ky<br />

fy<br />

Suspended<br />

mass<br />

Track<br />

10mm<br />

Fig 1 Diagram <strong>of</strong> <strong>the</strong> simplified rail vehicle model with a single wheelset<br />

and suspended mass<br />

The differential equations <strong>of</strong> <strong>the</strong> lateral position and yaw<br />

angle <strong>of</strong> <strong>the</strong> wheelset and <strong>the</strong> mass lateral position are as<br />

shown below.<br />

y = 1 m − 2f !!<br />

v<br />

ψ = 1 I<br />

−2f !!<br />

y + r !λ<br />

l<br />

lλ<br />

r !<br />

− Wλ<br />

l<br />

y − r !λ<br />

l<br />

d − vψ − 2f !"<br />

v<br />

ψ<br />

y − d) + F !" (1)<br />

y − d + l! v ψ<br />

+ 2f !"<br />

v<br />

− 2f !"<br />

v<br />

y − r !λ<br />

l<br />

y + r !λ<br />

d − vψ<br />

l<br />

ψ + lWλψ − I !"vλ<br />

lr !<br />

y + I !"vλ<br />

d<br />

lr !<br />

+ M !" (2)<br />

y ! = !<br />

! !<br />

−F !" (3)<br />

1


SYSTEMS DIVISION RESEARCH 2013 2<br />

The suspension forces and <strong>the</strong> yaw moments are given as<br />

F !" = k ! y ! − y + f ! y ! − y (4)<br />

M !" = −k ! ψ − f ! ψ (5)<br />

A complete list <strong>of</strong> <strong>the</strong> parameters and states used in<br />

equations 1-5 are given in Table 1.<br />

Parameters Values Units<br />

f 11<br />

f 22<br />

f 23<br />

f 33<br />

f y<br />

fψ<br />

k y<br />

kψ<br />

m<br />

mm<br />

I<br />

Iwy<br />

l<br />

v<br />

W<br />

r0<br />

F sy<br />

M sψ<br />

y<br />

y m<br />

ψ<br />

λ<br />

d<br />

Longitudinal creep coefficient<br />

Lateral creep coefficient<br />

Lateral/spin creep coefficient<br />

Spin creep coefficient<br />

Lateral damper coefficient<br />

Yaw damper coefficient<br />

Lateral suspension stiffness<br />

Yaw suspension stiffness<br />

Wheelset mass<br />

Suspended mass<br />

Wheelset yaw inertia<br />

Wheelset roll inertia<br />

Half wheelset width<br />

Velocity<br />

Wheelset width<br />

Nominal rolling radius<br />

Suspension lateral force<br />

Yaw moment<br />

States<br />

Wheelset lateral position<br />

Suspended mass lateral position<br />

Wheelset yaw angle<br />

Conicity<br />

Lateral track displacement<br />

7.44e6<br />

6.79e6<br />

13.7e3<br />

0<br />

50e3<br />

0<br />

0.23e6<br />

2.5e6<br />

1250<br />

8000<br />

700<br />

200<br />

0.7452<br />

20<br />

12263<br />

0.45<br />

-<br />

-<br />

-<br />

-<br />

-<br />

-<br />

-<br />

N<br />

N<br />

N<br />

N<br />

Ns/m<br />

Ns/rad<br />

N/m<br />

N/rad<br />

kg<br />

kg<br />

kg/m 2<br />

kg/m 2<br />

m<br />

m/s<br />

N<br />

m<br />

N<br />

Nm/rad<br />

m<br />

m<br />

rad<br />

m<br />

Fig 3 Pole Zero plot to show <strong>the</strong> stability <strong>of</strong> <strong>the</strong> wheelset for a step change<br />

in <strong>the</strong> track position<br />

It is clear that <strong>the</strong> wheelset shows a stable response for a<br />

10mm lateral track displacement. This can be verified from<br />

fig 3 as <strong>the</strong> poles are located on <strong>the</strong> left half <strong>of</strong> <strong>the</strong> s-plane.<br />

The pole zero plot was obtained by state space analysis <strong>of</strong><br />

<strong>the</strong> rail vehicle model to determine <strong>the</strong> stability <strong>of</strong> <strong>the</strong><br />

system.<br />

Additionally, <strong>the</strong> effect <strong>of</strong> primary yaw stiffness (k ψ) on <strong>the</strong><br />

stability <strong>of</strong> <strong>the</strong> wheelset has been studied. This was<br />

performed by varying <strong>the</strong> value <strong>of</strong> <strong>the</strong> primary yaw stiffness<br />

and determining how <strong>the</strong> stability level varies with <strong>the</strong><br />

change in <strong>the</strong> stiffness. Initially, <strong>the</strong> value <strong>of</strong> yaw suspension<br />

stiffness was assigned as 2.5x10 6 Nm, which when <strong>the</strong><br />

system remained stable with a single oscillation damping<br />

out quickly as shown in fig 2.<br />

When <strong>the</strong> stiffness was reduced to a very low value <strong>of</strong> 2500<br />

Nm, <strong>the</strong> response became unsteady and <strong>the</strong> wheel sets were<br />

unstable. This can be made evident from fig 4 below.<br />

Table 1 Parameter values and states used in <strong>the</strong> simulation model<br />

Figure 2 shows <strong>the</strong> lateral position <strong>of</strong> <strong>the</strong> wheelset and <strong>the</strong><br />

mass and yaw angle for a step change <strong>of</strong> 10mm in <strong>the</strong> track<br />

at a step time <strong>of</strong> 1 second. This shows that <strong>the</strong> position <strong>of</strong><br />

<strong>the</strong> wheelset and <strong>the</strong> suspended mass follows <strong>the</strong> step<br />

change in <strong>the</strong> track position.<br />

Fig 4 Lateral position <strong>of</strong> <strong>the</strong> wheelset and mass, yaw angle for <strong>the</strong> change<br />

in primary yaw stiffness <strong>of</strong> 2.5x10 3 Nm<br />

The above result can be fur<strong>the</strong>r analysed using <strong>the</strong> pole zero<br />

analysis, which is as shown in fig 5.<br />

Fig 2 Lateral position <strong>of</strong> <strong>the</strong> wheelset and mass, yaw angle for <strong>the</strong><br />

corresponding track displacement<br />

2


SYSTEMS DIVISION RESEARCH 2013 3<br />

be verified using <strong>the</strong> pole zero analysis as in fig 7.<br />

It can be understood from <strong>the</strong> above results that <strong>the</strong> change<br />

in <strong>the</strong> primary yaw stiffness can have influence in <strong>the</strong><br />

stability <strong>of</strong> <strong>the</strong> wheelsets.<br />

Fig 5 Pole Zero plot shows <strong>the</strong> instability <strong>of</strong> <strong>the</strong> wheelset for <strong>the</strong> yaw<br />

stiffness <strong>of</strong> 2.5x10 3 Nm<br />

The response <strong>of</strong> <strong>the</strong> system was again observed by<br />

increasing <strong>the</strong> value <strong>of</strong> primary yaw stiffness to 2.5x10 9 .<br />

This has been carried out to examine <strong>the</strong> system’s stability<br />

for a massive range <strong>of</strong> stiffness.<br />

IV. FUTURE WORK<br />

The simplified rail vehicle modelling has been followed by<br />

<strong>the</strong> simulation <strong>of</strong> a linearised half vehicle model consisting<br />

<strong>of</strong> two wheelsets, a bogie and a half vehicle body. The next<br />

phase <strong>of</strong> <strong>the</strong> project will be increasing <strong>the</strong> complexity in <strong>the</strong><br />

vehicle suspension system and <strong>the</strong> contact mechanics in a<br />

non-linear vehicle model. One <strong>of</strong> <strong>the</strong> several advanced<br />

techniques will be applied for <strong>the</strong> condition monitoring <strong>of</strong><br />

<strong>the</strong> railway vehicle.<br />

V. CONCLUSION<br />

In order to understand <strong>the</strong> dynamic interaction between<br />

wheel and rail, <strong>the</strong> initial stage <strong>of</strong> <strong>the</strong> project was to design a<br />

simplified model <strong>of</strong> a rail vehicle. The results produced<br />

from <strong>the</strong> single wheelset and suspended mass modelling<br />

were as expected and <strong>the</strong> stability <strong>of</strong> <strong>the</strong> wheelset with <strong>the</strong><br />

step change in <strong>the</strong> track lateral position have also been<br />

examined. An estimator model will be developed for <strong>the</strong> real<br />

time estimation <strong>of</strong> low adhesion at <strong>the</strong> wheel rail interface.<br />

Advanced filtering techniques will be analysed for <strong>the</strong><br />

estimation <strong>of</strong> rail vehicle dynamics.<br />

ACKNOWLEDGEMENT<br />

This work is being funded by Engineering and Physical<br />

Sciences Research Council <strong>of</strong> UK (EPSRC).<br />

Fig 6 Wheelset and suspended mass response for <strong>the</strong> varied primary yaw<br />

stiffness <strong>of</strong> 2.5x10 9<br />

REFERENCES<br />

[1] Charles, G., Goodall, R., and Dixon, R. (2008b) “A Least Mean<br />

Squared Approach to Wheel Rail Pr<strong>of</strong>ile Estimation,” in Proceedings<br />

<strong>of</strong> <strong>the</strong> 4 th International Conference on <strong>Railway</strong> Condition Monitoring.<br />

[2] Charles, G., Goodall, R., Dixon, R., Model based condition<br />

monitoring at <strong>the</strong> wheel rail interface <strong>Vehicle</strong> System Dynamics,<br />

46(1), 415–430.<br />

[3] Garg, V. and Dukkipati, R., Dynamics <strong>of</strong> <strong>Railway</strong> <strong>Vehicle</strong> System,<br />

Academic Press, 1984.<br />

[4] Iwnicki S. (Ed.) Handbook <strong>of</strong> railway vehicles dynamics, 2006, p. 535<br />

(Tailor & Francis Group, Boca Raton, London, New York).<br />

[5] Kalman, R.E, A new approach to linear filtering and prediction,<br />

Tran.ASME.J.Basic Eng.82 (1960), pp.35-45.<br />

[6] Ward, C., Goodall, R., and Dixon, R., Wheel rail pr<strong>of</strong>ile condition<br />

monitoring, in Proceedings <strong>of</strong> <strong>the</strong> UKACC control conference,<br />

Coventry (2010a).<br />

[7] Ward, C., Goodall, R., Dixon, R., and Charles, G.A., Adhesion<br />

Estimation at <strong>the</strong> wheel rail interface using advanced model based<br />

filtering <strong>Vehicle</strong> System Dynamics Vol. 00, No. 00, January 2012, 1–<br />

23<br />

Fig 7 Pole Zero analysis for <strong>the</strong> corresponding change in yaw stiffness<br />

(2.5x10 9 )<br />

The wheelset takes a very long time (say few minutes) to<br />

follow <strong>the</strong> track’s lateral shift as shown in fig 6. The impact<br />

was due to <strong>the</strong> increasing in primary yaw stiffness to a huge<br />

range, making <strong>the</strong> system stiffer. However, <strong>the</strong> variation had<br />

not affected <strong>the</strong> stability <strong>of</strong> <strong>the</strong> rail vehicle system. This can<br />

3

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!