Chaos in a railway bogie

  • Published on
    13-Aug-2016

  • View
    214

  • Download
    0

Transcript

Acta Mechanica 61, 89--107 (1986) ACTA MECHANICA 9 by Springer-Verlag 1986 Chaos in a Railway Bogie By Ch. Kaas-Petersen, 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 wheel-rail contact forces and nonlinear wheel flange-rail 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. Kaas-Petersen: 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 kg-m 2 and of the frame it is Is~ ~ 6780kg-m ~. The roll motion inertia of the frame is I Ir ~ 6780 kg-m 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 wheel-rail contact forces, the wheel flange-rail 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. Kaas-Petersen: 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 non-ideal 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 8Chaos 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 kN-s/m, k~ = 3.646 M_N/m, k~ = 0.3333 MN/m, D~ = 29.2 kN-s/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. Kaas-Petersen: 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) , OChaos 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,X2--x3-state 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. Kaas-Petersen: 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 . . . . . e-in} 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, @, l98 Ch. Kaas-Pe~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 bi-periodie solution. When we have bifurcation to a torus we can determine the stability of the bifurcating hi-periodic 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]. Bi-periodic 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. Kaas-Petersen: (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 bi-periodic 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. Kaas-Petersen: 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 bi-periodic solution. But we know, that the bifurcating bi-periodic 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 bi-periodic 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 E-01 - -0 .378693 E-01 - -0 .128728 E-01 - -0 .797676 E 00 - -0 .879224 E 00 - -0 .882330 E 0O - -0 .247277 E-01 - -0 .526221 E-01 - -0 .970951 E-02 - -0 .377638 E 00 0.291463 E 00 0.379779 E 00 0.641233 E-01 0.635954 E-01 - -0 .769662 E-02 0.552782 E 00 - -0 .765668 E 00 - -0 .482514 E 00 - -0 .257410 E-01 - -0 .468173 E-01 - -0 .109587 E-01 - -0 .486861 E 00 0.655995 E-01 0.480364 E 00 0.503040 E-01 0.210775 E-01 - -0 .732148 E-02 - -0 .238742 E-03 - -0 .413832 E 00 0.440516 E-01 - -0 .201531 E-01 - -0 .437 275 E-01 - -0 .868 681 E-02 - -0 .505552 E 00 0.852630 E-01 0.450742 E 00 0.453666 E-02 - -0 .470156 E-02 - -0 .810282 E-02 - -0 .607333 E-01 - -0 .874878 E-01 0.519966 E-01 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 E-02 v 1 --~ 0.7787 E-02 x~ ---- 0.61018 69662 E 00 ve ~ 0.2672 E 00 x 3 = 0.86043 53557 E-02 v a ~ 0.5340 E-02 xa ~ 0.45971 55997 E 00 v 4 ~ 0.8825 E-01 x 5 ~ - -0 .34572 21142 E-02 v 5 ~ - -0 .8159 E-02 x 6 ~ 0.33816 00148 E-01 v 6 = - -0 .8993 E 00 x~ ~ 0.96269 91552 E-02 v 7 ~ 0.3927 E-02 x s ~ 0.50304 33566 E 00 v s = 0.1175 E 00 x 9 ~ - -0 .43702 73831 E-02 v~ = 0.1683 E-02 xl0 ~ 0.37582 34619 E 00 vl0 = 0.2527 E 00 xl l = 0.97162 16562 E-02 vll = 0.3511 E-02 x~2 ~ 0.34642 77889 E 00 vl~ ~ - -0 .1668 E 00 x13 ~ - -0 .36333 62197 E-02 via ~ - -0 .1685 E-03 x~4 ---- 0.42574 63258 E 00 vl~ ~ 0.8018 E-01 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. Kaas-Petersen: 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, 752--762 (1972). [2] True, H., Kaas-Petersen, Chr.: A bifurcation analysis of nonlinear oscillations in railway vehicles. Vehicle System Dynamics 12, 5--6 (1983), and the full proceedings ed. by Paeejka, J. B., pp. 655--665. Lisse: Swets and Zeitlinger 1984. [3] Kaas-Petersen, Chr.: Dynamical behaviour of railway vehicles, Proe. EUROMECH Colloquium 190, pp. 41--43, (1985). [4] Kaas-Petersen, Chr., True, H. : Periodic, bi-periodie and chaotic dyneomical behavior of railway vehicles. Vehicle System Dynamics 14, 194--196 (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, 338--340 (1964). [6] Nagurka, M. L.: Curving performance of rail passenger vehicles, PhD-thesis, MIT, 1984. Chaos in a Railway Bogie 107 [7] Haken, H. : Advanced synergetics. Berlin--Heidelberg--New York--Tokyo: Springer 1983. [8] Liehtenberg, A. J., Lieberman, M. A.: Regular and stochastic motion. New York-- Heidelberg--Berlin: Springer 1983. [9] Wolf, A., Swift, J. B., Swinney, H. L., Vastano, J. A. : Determining Lyapunov expo- nents from time series. Physica D 16, 285--317 (1985). [10] Kaas-Petersen, 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 York--Berlin--Heidelberg--Tokyo: 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, 136--148 (1983). [13] Gear, C. W., Osterby, 0.: Solving ordinary differential equations with discontinuities. ACM Trans. Math. Software 10, 23--44 (1984). Dr. Ch. Kaas-Petersen Laboratory of Applied Mathematical Physics The Technical University o/Denmark DK-2800 Lyngby Denmar~