Published on
14Aug2016
View
214
Download
0
Transcript
Acta Mechanica 61, 89107 (1986) ACTA MECHANICA 9 by SpringerVerlag 1986
Chaos in a Railway Bogie
By
Ch. KaasPetersen, Lyngby, Denmark
With 10 Figures
(Received September 13, 1985; revised October ld, 1985)
Summary
Railway bogies may perform sustained lateral oscillations  the hunting motion  when the speed reaches a certain value. We examine the hunting motion in the complex Cooperrider bogie and find that the wheel flanges have a strong influence on the behaviour. When the bogie has flangeless wheels we have found a symmetry breaking bifurcation, by which we mean, that there is a transition from a symmetric periodic oscillation to an asymmetric periodic oscillation. The periodic motions are examined using the residual map. When the bogie has flanged wheels (the flange is represented as a dead band spring), we do instead find chaotic behaviour confirmed by a positive Liapunov number. The dynamical equations for the bogie model are strongly nonlinear, and we use computer methods to examine the dynamical behaviour.
1. Introduction
We examine the dynamical motion of a railway bogie running with constant velocity V along a straight, horizontal and perfect track. The resulting dynamical motion is a coupling between seven basic motions: The lateral motion (i.e. dis placement transversal to the track) and yaw motion (i.e. rotation around a vertical axis) of the front wheel axle, the lateral and yaw motions of the rear wheel axle, and the lateral, yaw, and roll motions (i.e. rotation around a horizontal axis parallel to the track) of the bogie frame. The equations of motion are nonlinear because of nonl inear wheelrail contact forces and nonlinear wheel flangerail forces. I f the motion of the bogie is a periodic oscillation we have the socalled hunting behaviour. We are interested in the asymptot ic motion, i.e. the motion when the transient has died away.
We state the equations decribing the dynamical motion of the bogie. The equations are due to Cooperrider [1]. The equations are formulated as a system of 14 first order autonomous nonlinear ordinary differential equations. We deter mine the steady state periodic motion by solving these equations numerically
90 Ch. KaasPetersen:
and use the residual map to eliminate the transient. We describe how we define and use the residual map. We next examine the steady state periodic motion in dependence of the velocity V of the bogie.
When the wheels are flangeless we have found only periodic motions of which some are symmetric and some are asymmetric. These two types of motions have been found to be connected to each other in a symmetry breaking bifurcation. We show this in a solution diagram and we show both types of periodic motions.
With flanged wheels we have found periodic motions and chaotic motions. A motion is said to be chaotic when the maximal Liapunov number is positive. We describe what the Liapunov number means geometrically and how to determine it.
Some preliminary results have already been presented. In [2] we have dealt with the Hopf bifurcation where the periodic motion is bifurcating from the stationary solution. In [3], [4] some further results on the periodic and chaotic motions have been given. But because of some errors in the equations we used, we have redone and extended the analysis, the results of which are presented in this paper.
2. The Bogie Model
The bogie we shall describe is the complex bogie due to Cooperrider [1]. The bogie consists of three stiff elements: Two wheel axles and a bogie frame. The wheels are rolling on the rails. The wheel axles are connected to the frame with springs and dampers. The frame is connected to the car body with springs and dampers. The car body is moving with constant velocity V along the track.
We allow 7 degrees of freedom in this model, namely lateral and yaw motion of the front axle, the rear axle, and the frame, and roll motion of the frame. By lateral motion we mean displacement in. a horizontal plane orthogona] to the track direction; by yaw motion we mean rotation in a horizontal plane around a vertical axis; and by roll motion we mean rotation in a vertical plane around a horizontal axis parallel to the track. In a coordinate system moving with the bogie these 7 state variables are assumed to be small and will be denoted ql . . . . . q7; see Fig. 1. We want to examine ql, ..., q7 as functions of time t. In order to do this we for mulate the dynamical equations using Newton's second law. The mass of the wheel axle is mw ~ 1022 kg and of the frame it is ms ~ 2918 kg. The moment of inertia for the yaw motion of the wheel axle is I~y ~ 678 kgm 2 and of the frame it is Is~ ~ 6780kgm ~. The roll motion inertia of the frame is I Ir ~ 6780 kgm 2.
The geometry of the bogie is described by the parameters a = 0.716 m being half the track gauge, b  1.074 m being half the axle distance and the distances dl ~ 0.620 m, d~ ~ 0.680 m, h~ ~ 0.0762 m and h2 ~ 0.6584 m from the springs or dampers to u centre of gravity; see Fig. 1.
There are three types of forces acting on this bogie: The wheelrail contact forces, the wheel flangerail forces and the spring and damper forces,
Chaos in a Railway Bogie 91
f
2D~ k 6 2k5 ~I/11/ / / / / / / t l / / / / / I / / / / / / / / /
V >
~V

2d 2
~f~111111III//IIf1111I///I/1111///////////////////////////////////////~
DI~ k5 :~k6 k5*"~Dt ~ ~k4
kz, k , D2 h2
ql
Fig. 1. The bogie model with the notation used. The wheels are shown flangeless
92 Ch. KaasPetersen:
The wheels and the rails are in contact and create contact forces   here we describe how these forces are determined. The wheels and the rails are made of steel with shear modulus G = 8.08. l0 s N/m 2. The wheels are conical with eonicity ~ = 0.05 and the radius r0 = 0.4572 m (when the wheel is centred on the track). The railhead is circular with radius 0.21 m. Hertz theory then tells, that the contact area is elliptical, and we have found the major and minor semi axes of that ellipse to be 6.578 mm and 3.934 mm respectively, when the vertical load (equal to one eights of the total weight of one railway vehicle i.e. two bogies and a car body) in one contact area is 66 670 N. Then we can use Vermeulen and Johnson's theory [5] to the nonideal rolling motion. They show, that in the contact area the wheel slides relative to the rail. This sliding is denoted creep and results in creep forces acting on the wheel. I f $, is the creep lateral to the track and ~v is the creep longitudinal to the track, then Vermeulen and Johnson derive, that
the creep force lateral to the track: Fx     
the creep force longitudinal to the track : Fy =    
Here ~:~ is the resulting creep
~ FR
~ = V(~/~,) ~ +($~/~)~,
and we have found T1 : 0.54219 and ~ : 0.60252, cf. [1], [5]. The resulting force F~ is given by
1 1 u ~4 u 3 for u 3, /~N
and here ~ = 0.15 is the friction coefficient, 2t is the vertical toad in the contact: area, and a, b are the semiaxes of the ellipse. (Thus G:~ab= 6.563 MN and ~N = 10 kN.)
The flanges of the wheels will contact the rail when the lateral excursion exceeds the clearance 6 = 0.0091 m, and usually we have a two point contact problem [6]. However, Cooperrider [1] reduces this to a one point contact problem by using a stiff spring with a dead band. The flange force FT is then
[ k0 . (u  8 ) for 8
Chaos in a l~ailway Bogie 93
The spring (damper) forces depend linearly on the relative displacements
(velocities). We have used
kl = 1.823 MN/m,
k4 = 0.1823 MN/m,
D1 = 20.0 kNs/m,
k~ = 3.646 M_N/m,
k~ = 0.3333 MN/m,
D~ = 29.2 kNs/m.
k3 = 3.646 MN/m,
ks = 2.710 MN.m,
The governing equations are listed below. (Here a = 0.716 m and b = 1.074 m.)
m,,(~ + Ax + 2F~f + Fr(ql) = 0
I,~,,~ + A~ + 2aNy! = 0
m~3 + A~ + 2F~,+ Fr(q3) = 0
Iwv~, + A4 + 2aF~ = 0
mill5   AI   A2 + As = 0
Itv~s   bA~ + bA~   Aa   A , +As =0
where
A1 = 2kt(ql   q5  bqs   htqT)
A~ = 2kl(q~   qs + bq6   hlq,)
A3 = 2k2dle(q2   q6)
A4 = 2k2d12(q~  qs)
A~ = 2D~(~5   h.i~,) + 2~,(q5   h~qT)
A, = ksqs
A~ = 2Dld2~ + 2ksd2~q~ + 4kadl2q~
Fxt, Fvl are creep forces for front axle, and the creep terms are ~e = ~I/V   q2,
~v = a~/V + ~q~/ro
F~,, Fu, are creep forces for rear axle, and the creep terms are ~:x = ~3/V   q4,
~, = ~. lV + ~q~/~o
z~T is the flange force.
We note in passing, that Cooperrider's simple bogie model [1]   a two degrees of freedom bogie with lateral and yaw motion of the full bogie as one stiff system   can be derived from the equations above by the substitution
(q, q2, q~, q,, qs, qs, qT) > (q5 + bqs, qs, q5   bqs, q~, qs, qs, 0).
94: Ch. KaasPetersen:
With the substitution x l = ql, x2 ~ ~[1, xa ~ q~, x4 = ~'2 . . . . we obtain a system of 14 ordinary differential equations (ODEs) of first order. We shall use the velocity V as parameter and take all other geometrical and physical quantities as constants. The dynamical system above can therefore be rewritten on the form
dx  = 2=/ (x ;V) , O
Chaos in a Railway Bogie 95
d v
t=o
v d 
d 6
Fig. 2. A piece o~ a chaotic solution starting in x~ at time t = 0. The chaotic solution remains bounded in the Xl,X2x3state space. Nearby solutions (initially the distance away but after a time interval ~ a distance d away from the chaotic solution) are shown.
The unit vectors v are shown; the vectors w are not shown
time, where we compute the Vector
w = [ r x0 + ~v)  r Xo)]l~.
Since 6 is Small, then w is an approx imat ion to the image of v for the flow l inearized
around the chaotic solution, w is now renormal ized to a unit vector
d= llwll,
v = w/d ,
where d is the factor with which v has been magnif ied. A f irst approx imat ion to
the growth rate is therefore
~u] = _1 In (d).
We now take the point ~b(~; x0) as init ial point, and solve the ODEs for a new time
nterva l of length v, then take a point close to it, ~b(v; x0) k 8v where v is the
vector determined after the f irst renormal izat ion, f ind the new vector w, the growth factor d, the vector v and the second approx imat ion to the exponential
divergence: ~[2] = 1 In (d). And so on and so forth. T
The largest L iapunov number can then be computed as the averaged value of
al l the growth rates
1 ~r
.Y*co ~' j= l
The number of renormal izat ion steps needed depends on the convergence rate of the series, of how close x0 is to the chaotic solution, and the vector v we choose to point in the direct ion of exponent ia l divergence. The renormal izat ion t ime
96 Ch. KaasPetersen:
should be short enough to catch the exponent ia l growth, and long enough to avoid cancel lat ion of digits in the computat ions of w.
A periodic solution will have a  0 within computat ional accuracy, and the
vector v will converge to the velocity vector of the solution. A quasi periodic
solution will also have a ~ 0 and the vector v will converge to a vector, tangent
to the corresponding toms surface in the state space.
4. Per iod ic Bchav iour
For a periodic solution we can do much better than comput ing the L iapunov
number. The characterist ic feature of a periodic solution is that the solution
repeats itself after a t ime interval T, the period. And since the dynamica l system
is autonomous, if we know one point of the solution, we can take that point as
init ial po int a t t ime t  0 and solve the ODEs to generate all of the periodic
solution.
..N
"/qJffR:u) "~/
~ig. 3. Illustration of construction of the residual map Q in R a. In the hypersurface H we take u as an initial point at time t = 0. The solution ~b(t; u) starting in u generates the curve shown. At time t ~ T that curve hits H again. The difference vector q is different from the zero vector so u is not a point on the periodic solution. If H (i.e. x~ and np) and u are chosen suitably, then Newton's method can be used to find a point x/, such that ~b(TR; xi,) ~ xp and xp is then a point on the periodic solution, xiD is not shown on the
sketch
We will now introduce (Fig. 3) the residual map Q, which makes use of the feature above. Assume for a moment, that we know xg to be a point close to the periodic solution, and let r~p be a unit vector paral lel to the veloci ty vector / (x~) . Then we
can define the hypersufface H through xH with np as normal. Let u be a point in H. (We could take u ~ xH.) The solut ion of the ODEs start ing in u will after a certain t ime interval Ti~ (this return t ime will be close to the per iod T) again be
in H ; the return point is r u). Let q be the difference vector
q = ~(T~; u)   u .
Chaos in a Railway Bogie 97
q is a vector entirely in H, so we can define the residual map Q mapping points in
H into vectors in H:
Q: H ~~ H, Q(u) = q.
We see, that if xp is a point in H on the periodic solution, then Te '= T, and Q(xp) = ~b(T; xp)  xp = 0 and thus a periodic solution is a zero point of Q. Note that Q is defined on H, and we have
H = x~r ~ span {el . . . . . ein}
where e~ .. . . . e,_~ is an orthonormal basis orthogonal to np. Another way to represent a point u in H is by the equation
(~  x . , ~p) = 0 ,
where (, .) is the usual euclidean scalar product. Usually we do not know XP, but if q~ is a sufficiently good approximation then
we can use Newton's method to find xp. We than need the derivative DQ of Q. Now DQ is an (n  1) (n  1) matrix, where the i 'th element of the ]'th column is given by the expression
([Q(u ~ ~ei)  Q(u)]/d, @,
l
98 Ch. KaasPe~ersen:
Continuation methods and stability analysis can be combined to determine critical points or bifurcation points, where the periodic solution changes its stability. Changes in the stability can come about in several ways. If one real eigenvalue lies on the unit circle in the point + 1, we may either have a symmetry breaking bifurcation or a saddle node bifurcation or no periodic solution any more (we are then in a ttopf bifurcation point). If one real eigenvalue lies in 1 we have a period doubling bifurcation. If one pair of complex conjugated eigenvalues lies on the unit circle we have bifurcation to a torus, also denoted a biperiodie solution. When we have bifurcation to a torus we can determine the stability of the bifurcating hiperiodic solution by the method described in [10]. All these types of bifurcations appear in the bogie model. Further details on the different types of bifurcations can be found in [11]. Biperiodic solutions can also be reformulated as zero point problems of a certain residual map [10]. However, we have not yet any results for the torus solutions in the bogie model treated in this paper.
5. Numerical Details
We have solved the ODEs with the routine LSODA [12] whose most important feature is that it automatically switches between stiff and nonstiff solution methods whenever needed. (It also has variable steplength and variable order.) When flange forces are included, the function expressions in the ODEs are continuous but not differentiable. We have used LSODA in this case as well, but we are aware, that a better method exists [13]. We have used absolute error control with tolerance 10 s, working in double precision on an IBM 3033.
In order to determine periodic solutions we determine the point r in the return plane H as a zero point of the function
(qb(Te; u)  xH, ni~), with T~ close to T.
We accept a point q~(TR ; u) as a zero point, when
ll < *.
with s~ = 10 11.
The residual map Q was said to have a zero point xp when llQ(xp)ll
Chaos in a Railway Bogie 99
6. The Resu l ts
We first present the results for the bogie with f langeless wheels. The solution d iagrams for the periodic solutions are shown in Fig. 4. We show the ampl i tudes
(i.e. the maximal posit ive value of the lateral d isplacements) and the periods
of the periodic solutions in dependence of the veloc i ty V. A solid line means, that the periodic solut ion is stable, and a dashed line, that it is unstable. In A
a Hopf b i furcat ion takes place, where the zero solut ion x~ ~ x~ . . . . . x~4 = 0
max xl max x5 0,06  0.10 
0.05 
0 ,04
0.03 
0 .02
0.01 
0.00
I t ! t I t I t
_ IA ! !
0 50 100 I I I
150 200 250
0.08 
0.06 
0.04  B
i 0.02 t
t
0.00 tlA I J I 1 J
0 50 100 150 200 250
V V
max x9 0.06
0.04 B [
I I t
0.02 i t t t l
0.0o !~A ! 0 50 100
C
K__
T 0,6 
0.5
0~   B ~ C. .
t I t t \A
I I I 0 .2  I i i i i 150 200 250 0 50 100 150 200 250
V V Fig. 4. Solution diagrams for bogie with flangeless wheels. We show the maximal latera l displacement of the front axle xl, of the rear axle xs, and of the bogie frame x 9 in meters and period of the motion T in seconds versus the velocity V in meters per second. A, B and C are bifurcation points, and V A ~ 66.7 m/s, V B ~ 49.0 m/s, V c = 176.7 m/s. The stable symmetric solution at V = 175 m/s is shown in Fig. 5, and the stable asymmetric
solution at V = 190 m/s is shown in Fig. 6
7*
100 Ch. KaasPetersen:
(corresponding to the bogie running centered on the track) becomes unstable
and generates a subcritical periodic solution, which in a saddle node bifurcation
at B becomes stable. So here we have a hysteresis phenomenon: We will see in
the physical system a jump from the zero solution up to a periodic solution
when the velocity VA is reached, but we have to decrease the velocity to V~
in order to jump back to the zero solution. The periodic solution remains stable
unt i l G where we have a symmetry breaking bifurcation. We see on :Fig. 5 a
stable symmetric periodic solution for a velocity below Vc and we see in :Fig. 6
a stable asymmetric periodic solution for a velocity above Ft.
x l 0.05 
0.04 
0.05 
0.02 
0.01 
0.00 
0 .01 
 0.02 
0 .03 
0 .04  0 .05 
0
J A i
X5 0.06 
0.04 
0.02 
0 .00 
0 .02 
0 .04 
, I 0 .06 
2 0
t
  i
1
i !
2
t
0.10
0.05
0.00
0 .05
0 .10
X5
t! 0
i ' i
1 2
t
0.06
0.04
0,02
0.00
0 .02
0 .04
0 ,06
X9
J d v " i ' ~ i ' ~ i = i
0 I
Fig. 5. The stable symmetric solution at V = 175 m/s. The initial conditions at t ~ 0 are given in Table 1. x a is the yaw of the front axle in radians; xl, xa, and x~ as explained in Fig. 4. The maximal values of x z, xs, and x 9 and the period T of this solution have been
used to make one point on the solution diagram in Fig. 4
Chaos in a Railway Bogie 101
Next we present the results for the bogie with flanged wheels. The solution diagrams are shown in Fig. 7. From these drawings we see, that the flanges stabilize the unstable periodic solution generated in the t Iopf bifurcation pointA. The hysteresis phenomenon is therefore taking place in a narrow parameter interval. The solution curve has no unique tangent in point B, this is because the wheel flange force has no unique tangent for x = ~. In C we have a saddle node bifurcation where the path of periodic solutions bends back. We have not yet followed that path further. We did find a jump to another stable periodic solution. When we followed that path for decreasing values of the velocity we found in D a bifurcation point, where a biperiodic solution bifurcates. We have
0.05
0.04
0.03
0.02
0.01
0.00
0 .01
0.02
0 .03
0.04 0 .05
x l
0
X5 0.10 
i I 1
i ! 2
t
X3 0.06 
0 .04 
0.02 
0 .00 
0 .02 
0 .04 
0 .06 t ' .... 1 0 1 2
0.05 
0.00 
0 .05
0 .10
0 1 2
t
X9 0.06 
0.04 
0.02 
0 .00 
0 .02 
0.04 
0 .08 ,
0
i  J 1 2
t Fig. 6. The stable asymmetric solution at V = 190 m/s, else same text as to Fig. 5
102
max x l 0.02 
0.01
0.00
Ch. KaasPetersen:
E
! i I ! I I IA
i i I i t 50 100 150 200 250
V
max x5 0.02 
0.01 
0,00
C
I I I t I
I" t i t I 'i
50 100 150 200 250
V
max x9 0.020 
0 .015
0.010
0,005
0.000
T 0.4 /
; [ I i ! I
0,3
0.2
B
C F
',A = i i i .... / 0.1 I t I I !
0 50 100 150 200 250 0 50 100 150 200 250
V V
Fig. 7. Solution diagrams for bogie with flanged wheels. A, B, C, D, E, and F are bi furcation points, and V A = 66.7 m/s, VB = 65.4 m/s, V c = 113.9 m/s, VD = 112.2 m/s, V~ = 203.7 m/s, Vr  218.3 m/s. The stable periodic solution at F  175 m/s is shown
in Fig. 8
neither followed the path of periodic solutions for velocities below VD nor the biperiodic solution. But we know, that the bifurcating biperiodic solution is unstable, so its path bends to the right (not indicated on Fig. 7). The path of periodic solutions starting in D consists of stable solutions until E. On Fig. 8 we have shown the stable periodic solution for V = 175 m/s. Point E is a bifur cation point where unstable biperiodic solutions bifurcate (not indicated on Fig. 7). We have followed the path of unstable periodic solutions up to V ~ 250 m/s. In point F we met a bifurcation point where a period doubling bifurcation
Chaos in a Railway Bogie 103
0,02
0,01
0,00
0 ,01
0 .02 a
0,03
0,02
0.01
0.00
0 .01
0 .02
t I I l 0 ,03 2 3 4 5
x3
7
0 1 2 3 4 5
t t
x5 0.02 
0,01 
0 .00 
0 .01 
0.02 0
x9
t t t t 2 3 4 5
t
0.03 
0.02
0.01
0.00 
0 .01
0 .02 
0 ,03 '''~ 0
 I i I I i 1 1 2 3 4 5
Fig. 8. The stable solution at V = 175 m/s for the bogie with flanged wheels. From this solution we get one point on the solution diagrams in Fig. 7
took place. We have not yet determined the unstable periodic solut ion with the double period.
At the velocit ies V = 225 m/s and V = 250 m/s, both above VE, we have
found chaotic behaviour. The chaotic solut ion at V ~ 250 m/s is shown in Fig. 9. I t is there seen, that it is the rear axle which performs a seemingly random motion. Wi th the init ial condit ion and the vector v as given in Table 2 we show
in Fig. 10 the convergence of the L iapunov number. We cannot tell def init ively why we have chaotic motion, but the point F will p resumably p lay a key role.
This is because of the Fe igenbaum scenario [7], [8] where a periodic solution
104 Ch. Kaas Petersen :
Table 1. In i t ia l condition x ~ (xl, . . . , xla ) at t = 0 and period T o/ the periodic solutions shown
Fig. 5, V  175 m/s Fig. 6, V = 190 m/s :Fig. 8, V = 175 m/s
Xl ~ g l
X2 ~ ~1
X 3 ~ q~
X4 ~ q2
X5  q3
x~ = Cs X7 ~ ~4
X8 ~ ~4
X 9 ~ ~t 5
xn = qs
X12 ~ ~6
X13 ~ q7
X14 = ~7
T
0.165547 E01  0 .378693 E01  0 .128728 E01  0 .797676 E 00  0 .879224 E 00  0 .882330 E 0O  0 .247277 E01  0 .526221 E01  0 .970951 E02  0 .377638 E 00 0.291463 E 00 0.379779 E 00
0.641233 E01 0.635954 E01  0 .769662 E02 0.552782 E 00  0 .765668 E 00  0 .482514 E 00
 0 .257410 E01  0 .468173 E01  0 .109587 E01  0 .486861 E 00 0.655995 E01 0.480364 E 00
0.503040 E01 0.210775 E01  0 .732148 E02  0 .238742 E03  0 .413832 E 00 0.440516 E01  0 .201531 E01  0 .437 275 E01  0 .868 681 E02  0 .505552 E 00 0.852630 E01 0.450742 E 00
0.453666 E02  0 .470156 E02  0 .810282 E02  0 .607333 E01  0 .874878 E01 0.519966 E01
0.402 510 0.400730 0.224451
Table 2. In i t ia l condition x at t = 0 and vector v ]or chaotic solution at V = 250 m/s shown in .Fig. 9
xl ~ 0.48702 38147 E02 v 1 ~ 0.7787 E02 x~  0.61018 69662 E 00 ve ~ 0.2672 E 00 x 3 = 0.86043 53557 E02 v a ~ 0.5340 E02 xa ~ 0.45971 55997 E 00 v 4 ~ 0.8825 E01 x 5 ~  0 .34572 21142 E02 v 5 ~  0 .8159 E02 x 6 ~ 0.33816 00148 E01 v 6 =  0 .8993 E 00 x~ ~ 0.96269 91552 E02 v 7 ~ 0.3927 E02 x s ~ 0.50304 33566 E 00 v s = 0.1175 E 00 x 9 ~  0 .43702 73831 E02 v~ = 0.1683 E02
xl0 ~ 0.37582 34619 E 00 vl0 = 0.2527 E 00 xl l = 0.97162 16562 E02 vll = 0.3511 E02 x~2 ~ 0.34642 77889 E 00 vl~ ~  0 .1668 E 00 x13 ~  0 .36333 62197 E02 via ~  0 .1685 E03 x~4  0.42574 63258 E 00 vl~ ~ 0.8018 E01
makes a cascade of per iod doub l ing b i fu rcat ions ; the d i s t inc t b i fu rcat ion po in ts
accumulates , ar id beyond that accumulat ion po in t chaot i c behav iour is observed .
I f th i s scenar io ac tua l ly took p lace here , then the per iod ic so lu t ions wou ld a l l
be unstab le , but they may never the less have a s t rong in f luence on the chaot i c
mot ion .
Chaos in a Railway Bogie 105
0.,02
0.01
0.00
0 .01
0 .02
xl
i I 2 3 4
I 5
t
x3 0.03 
0.02 
0.01 
0.00 
0 .01 
0 .02
0.03 0
I i i i i 1 2 3 4 5
t
0.02
0.01
0.00
0 .01
0.02
x5
i
I I I i I 0 1 2 3 4 5
t
x9 0.03 
o.01 
o.0o 
0.01 
O.O2 
0.03 J i l I I 0 1 2 3 4 5
Fig. 9. The chaotic solution at V = 250 m/s. The initial conditions st t = 0 are listed in Table 2
7. Conc lus ions
We have examined the complex Cooperr ider bogie with f langed and with flangeless wheels running with constant velocity V along a straight track. The flanges, when present, were represented as dead band springs. The dynamica l mot ion of this bogie is described by a system of 14 first order autonomous non
l inear ODEs, and we have used computers to obta in the dynamica l motion. Apar t from the t topf bifurcation, where the s tat ionary zero solut ion bifurcates to an unstable periodic solution, the solution d iagrams for these two bogies are
106 Ch. KaasPetersen:
3
0 50 100 150 200
N F~. 10. Conve~ence of Liapunov number. The initial conditions at t = 0 and the unit
vector v are listed in Table 2
very different indicating the strong influence of the flange forces. Especially
we have found chaotic behaviour when the wheels are flanged, and the chaotic
behaviour is seen on the rear wheel axle.
Acknowledgments
The first version of the bogie chaos was presented on a poster session at the XVIth ICTAM meeting, August 1984. I thank Dr. O. Mahrenholtz for allowing me to present a second version [3] at the EUROMECH Colloquium 190, October 1984. Later, Dr. Sauvage, SNCF, kindly pointed out some errors in [2] so the calculations were redone and extended.
References
[1] Cooperrider, N. K.: The hunting behavior of conventional railway trucks. ASME J. Eng. Industry 94, 752762 (1972).
[2] True, H., KaasPetersen, Chr.: A bifurcation analysis of nonlinear oscillations in railway vehicles. Vehicle System Dynamics 12, 56 (1983), and the full proceedings ed. by Paeejka, J. B., pp. 655665. Lisse: Swets and Zeitlinger 1984.
[3] KaasPetersen, Chr.: Dynamical behaviour of railway vehicles, Proe. EUROMECH Colloquium 190, pp. 4143, (1985).
[4] KaasPetersen, Chr., True, H. : Periodic, biperiodie and chaotic dyneomical behavior of railway vehicles. Vehicle System Dynamics 14, 194196 (1985), and the full proceedings (to be published).
[5] Vermeulen, P. J., Johnson, K. L. : Contact of nonspherieal elastic bodies transmitting tangential forces. J. Appl. Mech. 81, 338340 (1964).
[6] Nagurka, M. L.: Curving performance of rail passenger vehicles, PhDthesis, MIT, 1984.
Chaos in a Railway Bogie 107
[7] Haken, H. : Advanced synergetics. BerlinHeidelbergNew YorkTokyo: Springer 1983.
[8] Liehtenberg, A. J., Lieberman, M. A.: Regular and stochastic motion. New York HeidelbergBerlin: Springer 1983.
[9] Wolf, A., Swift, J. B., Swinney, H. L., Vastano, J. A. : Determining Lyapunov expo nents from time series. Physica D 16, 285317 (1985).
[10] KaasPetersen, Chr. : Computation, continuation, and bifurcation of torus solutions for dissipative maps and ordinary differential equations. To be published in Physica D.
[11] Guckenheimer, J., Holmes, P.: Nonlinear oscillations, dynamical systems, and bi furcations of vector fields. New YorkBerlinHeidelbergTokyo: Springer 1983.
[12] Petzold, L.: Automatic selection of methods for solving stiff and nonstiff systems of ordinary differential equations. SIAS[ J. Sei. Star. Comput. 4, 136148 (1983).
[13] Gear, C. W., Osterby, 0.: Solving ordinary differential equations with discontinuities. ACM Trans. Math. Software 10, 2344 (1984).
Dr. Ch. KaasPetersen Laboratory of Applied Mathematical Physics
The Technical University o/Denmark DK2800 Lyngby
Denmar~