- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Table of Contents

Advances in Mechanical Engineering

Volume 2013 (2013), Article ID 296106, 19 pages

http://dx.doi.org/10.1155/2013/296106

## The Effect of First-Order Bending Resonance of Wheelset at High Speed on Wheel-Rail Contact Behavior

State Key Laboratory of Traction Power, Southwest Jiaotong University, Chengdu 610031, China

Received 29 June 2013; Accepted 3 September 2013

Academic Editor: Xiaoting Rui

Copyright © 2013 Shuoqiao Zhong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The first-order bending deformation of wheelset is considered in the modeling vehicle/track coupling dynamic system to investigate its effect on wheel/rail contact behavior. In considering the effect of the first-order bending resonance on the rolling contact of wheel/rail, a new wheel/rail contact model is derived in detail in the modeling vehicle/track coupling dynamic system, in which the many intermediate coordinate systems and complex coordinate system transformations are used. The bending mode shape and its corresponding frequency of the wheelset are obtained through the modal analysis by using commercial software ANSYS. The modal superposition method is used to solve the differential equations of wheelset motion considering its flexible deformation due to the first-order bending resonance. In order to verify the present model and clarify the influence of the first-order bending deformation of wheelset on wheel/track contact behavior, a harmonic track irregularity with a fixed wavelength and a white-noise roughness are, respectively used as the excitations in the two models of vehicle-rail coupling dynamic system, one considers the effect of wheelset bending deformation, and the other does not. The numerical results indicate that the wheelset first-order bending deformation has an influence on wheel/rail rolling contact behavior and is easily excited under wheel/rail roughness excitation.

#### 1. Introduction

The development of high-speed railway is associated with dynamic problems, for example, the wear of their wheels/rails. When a train with polygonal wheels due to wear operates on a track at high speed, the train and the track vibrate fiercely in an extensive frequency range. The wheel polygonal wear excites the strong dynamic behavior of the vehicle/track structures in mid- and high-frequency ranges. In the extensive experiments conducted on China metro lines by Jin et al. [1], the test results showed that all of the wheels of one type of metro trains in service present the 9th-order polygonal wear, and the excitation frequency of polygonal wheel was close to the first-order bending resonance frequency of the wheelset. The root of the 9th order polygonal wear was attributed to the first-order bending resonance occurring of wheelset in operation, as indicated in Figure 1. Now, the problem of wheel polygonal wear is very common on the existing high-speed trains in China. The mechanism of this phenomenon has not been explained perfectly until now, and it also has become one of the most difficult problems of the wheel-rail system to be solved. The wheel polygon due to wear usually takes a long time and so far cannot be reproduced through the numerical simulation due to the existing insufficiently advanced theoretical models for railway vehicle coupling with track.

The dynamics analysis of railway vehicle by using the existing rigid multibody dynamics models for vehicle/track coupling system is limited in the low-frequency domain of 0–20 Hz [2]. An effective way to broaden the analysis frequency range is to take the structural flexibility of vehicle/track into consideration in modeling the vehicle/track system [3–6]. And the key part in modeling structural flexibility of the vehicle/track is to consider the effect of the flexibility of wheelset on the rolling contact of wheelset/rail.

Many scholars have carried out a lot of researches on the vehicle/track dynamic system modeling that considers the effect of flexible wheelset. In [5], the wheelset was treated as a lumped and discrete structure by Popp et al. In [4, 6–8], the existing models of the flexible wheelset were continuous models. A vehicle/track dynamic model developed by Meinders considering wheelset flexibility was based on the cosimulation of finite element analysis software ANSYS and multibody dynamics simulation software NEWEUL [9, 10]. Baeza et al. [11–13] proposed a numerical dynamic system model originating from a rotating elastic cylinder model, taking inertial effect and effect of Coriolis force into account.

In terms of wheel/rail rolling contact modeling, the effect of wheelset deformation in mid- and high-frequency ranges should not be neglected [14]. Though many literatures presented different methods to model flexible wheelsets, the detailed derivations of rolling contact geometry of a flexible wheelset and a pair of flexible rails and the detailed calculations and discussions on the effect of the high-frequency bending of wheelset on the rolling contact behavior were not given. This is a very difficult research. The influence exerted by the wheelset bending formation is much concerned with wheelset polygonal wear. In [15], the effect of wheelset flexible deformation on the creepages and the creep forces in static state was indicated to be significant. The flexible deformation of a wheelset in service always contains mid- and high-frequency components. The rolling contact behavior of the wheel/rail caused by wheelset deformation would be in mid- and high-frequency ranges, which brings a great difficulty into the on-line calculation of wheel/rail contact geometry. And the existing models for the wheel/rail contact geometry are not suitable for the model with flexible wheelsets and need to be further improved.

In this paper, a new wheel/rail contact model is developed to consider the effect of the wheelset first-order bending deformation on the wheel/rail rolling contact behavior. The main object is to analyze the effect of wheelset flexibility on the wheel/rail contact behavior mainly in mid-frequency range. The mid-frequency components of the excitation from the wheel/rail interaction system are filtered by the suspension system, and the wheel/rail mid-frequency excitation cannot easily excite the flexibility of the car-body and bogies. Hence, the influence of the structural flexibility of car-body and bogies is ignorable in this present paper. In Section 2, the model of the flexible wheelset and the wheel/rail contact model are introduced. In Section 3, a harmonic track irregularity with a fixed wavelength and white-noise roughness of the rail are adopted as the excitations to analyze the effect of the wheelset flexibility.

#### 2. Modeling

##### 2.1. Vehicle/Track Model

Conventional vehicle system models consist of one carriage, two bogies, and four wheelsets, which are considered as rigid bodies or lumped masses connected by spring-damping systems. A ballasted track is usually modeled as a triple-layer model of discrete elastic support. The rail is modeled as a Timoshenko beam on elastic point supporting foundation. The two ends of the calculation rail are simply supported. The vertical and lateral bending deformations and twisting of the rails are taken into account. The sleepers are modeled as rigid bodies, and the ballasts are regarded as discrete equivalent bodies. The equivalent spring-damping systems are used as the connections between the rails and the sleepers, the sleepers and the equivalent ballast bodies, and the ballast bodies and the subgrade. Figure 2 illustrates such a model of the vehicle and the track. Its detailed description and the derivation of motion equations refer to the published papers, for example, [16–18]. The important parameters of the vehicle and the track used in this paper calculation are given in Table 1.

##### 2.2. Flexible Wheelset Model

The flexible wheelset modeling mainly considers the contribution of the first-order bending deformations of the wheelset axle in the two planes, perpendicular and parallel to the track level plane, respectively, and neglects their crossing effect. The wheelset modes are calculated by using the finite element software ANSYS. The corresponding resonance frequency of the first-bending, shown in Figure 3, is about 117 Hz. This mode function will be used in the analytical solution to the partial differential equations of the wheelset deformation based on the modal superposition method [16].

Figure 3 indicates the force system on the flexible wheelset. and are, respectively, the left and right points of the application of the primary suspension forces on the axle. and are, respectively, the left and right contact points of wheel/rail. indicates the origin of the coordinate system - that is a coordinate system with a translational motion along with a tangent track central line at the vehicle speed. If the vehicle speed is a constant, this coordinate system is an inertial coordinate system, therefore, regarded as an absolute coordinate system (geodetic coordinate system).

Based on the above assumptions that the wheels are regarded as the rigid bodies and perpendicular to the axle central lines at their centers, their spatial poses are determined by the positions of nominal rolling circle centers and the directions of the normal lines through nominal rolling circle centers on the deformed axle centerline. The wheelset axle is modeled as an Euler-Bernoulli beam. So according to the force analysis diagram of the wheelset given in Figure 3, the left and right wheel/rail forces are, respectively, translated to the nominal circle centers, and . Such a force translation introduces the additional moments applied to the axle. The additional moments are the products of the wheel/rail forces and the distance vector between the wheel/rail contact points and the rolling circle centers. Then, the force analysis diagrams in the two planes of the wheelset axle could be obtained, as shown in Figures 4(a) and 4(b), based on which the differential equations of the axle can be established.

The partial differential equations of the wheelset in bending deformation in the two planes, (1) and (2), give the equations for the perpendicular and the parallel planes to the track level plane, respectively, as follows

In (1) and (2), and are, respectively, given as

and are, respectively, given as

The function in (3)~(5) reads as

The boundary conditions of (1) and (2) are given as where “” means one-order partial derivative with respect to coordinate.

Based on the modal superposition principle, the solutions to (1) and (2) are, respectively, given as

Using the normalized modal functions which are mutually orthogonal, we have that and are, respectively, given as

In addition, the contribution of the rigid modes to the wheelset motion is also taken into account. The total displacements of the wheelset in the two planes, which are, respectively, parallel and perpendicular to the track level plane, read as

The elastic mode number is given as , and the rigid model number is given as . The rigid modal functions read as

The wheelset bending deformation in the two planes has an influence on the primary suspension forces in the directions of - and -axes, which read, respectively, as

The lateral primary suspension forces on the four wheelsets read as

The effect of the wheelset bending deformation on the lateral primary suspension forces is neglected. They are mainly influenced by the wheelset rigid motion in the -direction. Substituting (11) into (13), the components of the primary suspension forces in the directions of - and -axes are written as follows:

The flexible wheelset model is successfully coupled with the bogie through (14) and (15). The calculation formulas of the wheel/rail contact forces are to be discussed below.

##### 2.3. Wheel/Rail Rolling Contact Model

###### 2.3.1. Description of Coordinate Systems for Wheel/Rail Rolling Contact Calculation

It is complex and tedious to carry out the rolling contact behavior calculation of a flexible wheelset running on a pair of rails in service. The calculation is carried out with the aid of the transformations of multicoordinate systems. A wheel/rail rolling contact model considering flexible wheelset effect is much different from that without considering flexible wheelset effect.

Figure 5 indicates a detailed contact geometry description of a flexible wheelset in rolling contact with a pair of rails, the global coordinate system, and the local coordinate systems used in the rolling contact calculation. In Figure 4, the rolling angles, , are, actually, the projections of the rotation angles of the normal line of the axle cross-section at the wheel nominal rolling circle centers in the - plane, but the yawing angles of the wheels, which are the projections in the - plane, are not shown in the figure and can be described in the top view of Figure 5 that is omitted here.

The left and right contact points are regarded as, respectively, the origins and of the local coordinates -. and are the centers of the left and right instant rolling circles, respectively. The two auxiliary lines are created, and they are, respectively, perpendicular to the left and right instant circles and pass through the centers and . These two auxiliary lines, respectively, intersect with the common planes of the wheel/rail contact surfaces, and points and are, respectively, their intersection points, as indicated in Figure 5. Making a line connects point to the left contact point and the other line connects point to the left contact point . The included angle between line and line and the one between line , and line are exactly the contact angles and , respectively. Each coordinate system in Figures 2 and 4 has the following detailed description and requirements. (1): this coordinate system is illustrated in Section 2.2. Its base vector is that . Its unit vectors, , , and are, respectively, in the longitudinal, vertical and lateral directions of the tangent track. (2)-: it is a translating coordinate system that is always parallel to the coordinate system -, as indicated in Figure 5. It translates with the wheelset mass center. Its base vector is that .(3)-: they are the local coordinate systems attached to the two rigid wheels, as indicated in Figure 5. Their origins and are at the centers of the nominal rolling circles of the left and right wheels, respectively, which means that their poses are determined by the wheel poses and their axes are always perpendicular to the corresponding nominal rolling circle planes. The nominal rolling circles of the left and right wheels coincide with the planes - and -, respectively. But the rolling motions of the left and right wheels are around the axes and , respectively, and their yawing motions are around the axes and , respectively. So the axes and are always parallel to the horizontal plane, and the axes and are always parallel to the vertical plane. Their base vectors are that . The unit vectors , , and indicate the directions of the axes , , and , respectively.(4)-: their axes and , respectively, coincide with the lines and , as indicated in Figure 5. The planes - and - coincide with the common tangent planes of the contact surfaces of the left and right wheels/rails, respectively. The axes and coincide with the common normal lines of the contact surfaces of the left and right wheels/rails, respectively. Their base vectors are that . (5)-: these two coordinate systems are used to describe the configurations and motions of the left and right rail cross-sections which are assumed to be rigid. Their origins are fixed at the left and right rail cross-session centers, and , and their rotation and translation are determined by those of the rail cross-sections. Their base vectors are, that, .

The transformations from the base vector of coordinate system - to those of the coordinate systems - read as

The transformation from the wheel coordinate systems - to the coordinate system - is expressed in (17), in which are the coordinates of any point described in the coordinate system - or -, and is the coordinates of the same point described in the absolute coordinate system -. indicates the vector form the origin of - to the point . are the coordinates of the nominal rolling circle centers in the coordinate system -. It is noted that in (17) the bases are the same as , as indicated in Figure 4. The transformation in (17) includes two times of translations and one rotation of the coordinate systems -. Consider the following:

The transformations from the base vectors of the local coordinate systems - at the contact points to the absolute coordinate system read as where

The transformations from the bases of the local coordinate systems - and - at the rail cross-section centers to those of the coordinate system - (or -) are given as

The wheel/rail rolling contact calculation consists of three parts, wheel/rail contact points, wheel/rail normal forces, and wheel/rail tangential forces. They are elaborated as follows.

###### 2.3.2. Wheel/Rail Contact Point Calculation

Before the calculation of the wheel/rail contact points, the given discrete data points of the profiles of the wheel treads are usually described in the local coordinate system associated with the wheel tread, and the rail top profile described in the local coordinate system associated with the rail [19]. In the calculation, they have to be transformed into the expression in the absolute coordinate system - through a series of the coordinate system transformations which include rotating and translating the coordinate system.

The prescribed profile of the wheel is usually described in the local coordinate system associated with the wheel tread. The local coordinate system associated with the wheel tread is, actually, parallel to the rigid wheel coordinate system, -, or -, and their origins are at the lowest points of the nominal rolling circles. In order to obtain the coordinates of the wheel tread profiles described in the absolute coordinate system -, firstly, the coordinates of the profiles of the wheels described in the local coordinate systems associated with the wheel treads have a translation that equals the radius of the nominal rolling circles, and can be described in -. Then, they are transformed into the description in the absolute coordinate system - by using formula (17) and with the help of the wheelset coordinate system -.

The transformation from - to - is through translating and rotating the coordinates of the rail top profiles described in -. The translation is just done along the vector or , and (20) is the transformation from the base vector of the coordinate systems - to that of the coordinate system -.

The existing models for the calculation of rigid wheel/rail contact geometrical relationship need to be further improved to be applicable to the contact geometry calculation of a flexible wheelset in contact with a pair of rails. Figure 6 illustrates the spatial contact geometrical relationship of a wheel of a rigid wheelset in contact with a rail [19], in which the contact state on the right side is shown. In Figure 6, is the contact point of the right wheel/rail, is the center of the rigid wheelset, the left and right wheels have the same rolling angle and yawing angle, the wheelset axle centerline is perpendicular to the instant rolling circle plane, is the intersection point of the wheelset axle centerline and the common tangent line of the contact surfaces, is the distance between the point and the rolling circle, and are, respectively, the center and radius of the rolling circle, is the intersection point of the common normal line of the contact surfaces and the wheelset axle centerline, - is the coordinate system of the rigid wheelset, and are, respectively, the projections of on the - plane, so the line is the projection of the axle centerline on the - plane.

The wheel contact point trace line, on which the contact point is, in the absolute coordinate system - can be calculated with (21) according to [18] as follows: where are the coordinates of the instant rolling circle center in the absolute coordinate system, and they are given as follows: where are described in the coordinate system - defined as a rigid wheelset coordinate system. So the series of the wheel rolling circles can be determined with .

But, for the contact between a flexible wheelset and a pair of rails, due to the axle deformation, as indicated in Figure 5, strictly speaking, is not constant, and the directioncosine components, , of the central lines of the left and right wheels should be modified. in (22) and (23) should be replaced by the direction cosine components, , of the normal lines at the instant rolling circle centers in -; namely, in the calculation of , the rolling and yawing angles of the wheelset should be replaced by of the left and right wheels of the flexible wheelset. It should be noted that the normal lines through the instant rolling circle centers are, actually, those through the centers of the nominal rolling circles since the wheels are assumed to be rigid. For the flexible wheelset, the coordinates of the centers of the left and right instant rolling circles should be written as, based on formula (23) and Figure 6, where are the distances between the center of the flexible wheelset and the centers of the left and right instant rolling circles, respectively. And , , and are the increments of the wheelset center translation in the -, -, and -directions, respectively. They are caused by the bending deformation of the flexible wheelset. Comparing (23) with (24), for the flexible wheelset, the calculation of the coordinates of the wheelset center and the instant rolling circle center considers the combined effect of the rigid motion and flexible deformation of the wheelset. But its stretching and compressing deformations caused by the axle bending deformation in the axle direction are very small, compared with the axle bending deformation amplitudes. So is neglected. Since the analyzed wheelset is not a powered wheelset, the component along the track of the wheel/rail force is quite small, compared with the vertical component, the bending deformation of the wheelset is very small in a horizontal plane of the track. So caused by the wheelset bending deformation can also be neglected. It should be noted that are obtained through the dynamic analysis using the traditional multibody dynamic model of the vehicle/track, which is introduced in the present paper.

Once the contact circles and the rail top profiles are described in the absolute coordinate system -, it will be easy to calculate the minimum distances and between the wheel treads and rail top surfaces, as well as the coordinates of the wheel/rail contact points. The wheel/rail contact angles can be subsequently calculated with the known coordinates of the wheel/rail contact points and the given wheel/rail profiles.

###### 2.3.3. Wheel/Rail Normal Forces

The vertical instant compressions of wheel/rail at the contact points, and , are written as

The instant normal compressions, and , are calculated according to (26) and the contact geometrical relationship shown in Figure 6, and read they as The left and right wheel/rail normal forces are calculated by Then, the normal force should be decomposed into three components in the -, -, and -directions in establishing the motion equations of the wheelset and the rails.

###### 2.3.4. Wheel/Rail Creep Forces and Spin Moments

The wheel/rail creep forces are calculated by using the model developed by Shen-Hedrick-Elkins, as described in [18]. In the calculation of each time step, the wheel/rail creepages , , and should be calculated in advance. The existing creepage calculation formulas need to be improved and can be used to calculate the creepages of a flexible wheelset in rolling contact with a pair of rails. The creepages represent the ratios of the relative velocities of the contact surfaces at the contact points to the forward velocities, , of the instant rolling circles of the wheels, as shown in (28). Consider the following: where , , and are, respectively, the relative longitudinal velocity, the relative lateral velocity, and angular velocities. The relative velocities are firstly calculated in the absolute coordinate system and then transformed into the expressions in the local coordinate systems associated with the contact surfaces, which will be discussed in the following.

can be calculated in two ways. The first way is to make a summation of the rigid motion and the flexible deformation of the wheelset, and the other way is to take the derivatives of the coordinates of the centers of the instant rolling circles with respect to time directly. But, actually, the difference between and is quite small due to the very small effect of the rigid yaw, motion and the bending deformation of the wheelset on the forward speeds, , of the instant rolling circles. So, by using formulas (28), for simplicity, is replaced with . is the speed of the vehicle or the nominal forward speed of the wheelset center.

Based on the composite motion theorem of rigid body motion, the absolute velocity of a particle is the resultant vector of the relative velocity and the transport velocity , given in (29). In other words, is the velocity of the moving point, described in the moving coordinate system, and is the velocity of the moving coordinate system relative to the absolute coordinate system. Consider the following:

The wheel coordinate systems -, as shown in Figure 5, can be regarded as the moving coordinate systems. So the absolute velocities of the left and right contact points , respectively, on the left and right wheels, have the following expressions: where and are, respectively, the transport velocities and the relative velocities. The relative velocities are the velocities of the wheel contact points relative to the origins of the moving coordinate systems. The transport velocities are the velocities of two points in -, and these two points coincide with the contact points , respectively. So read as where indicate the vectors . indicate the derivatives of the vectors with respect to time. are the vectors . represent the translational velocities of the centers of the nominal rolling circles of the left and right wheels, described in the absolute coordinate system. are determined by the angular velocity vectors of the left and right wheels, and , as indicated in (32). It is noted that is just the rotation of the transport coordinate systems, -, relative to the absolute coordinate system -, as shown in Figure 5. Consider the following:

It should also be noted that do not include the rolling () of the rigid wheelset over the track and neglect the rolling angular velocities of the wheels due to the twisting deformation of the wheelset, which produce the longitudinal creepages. The relative velocities are the velocities of the left and right contact points described in the local coordinate systems -, due to the rolling of the wheelset over the two rails. They are written as where Equation (34) is expressed in the local coordinate systems -.

Substituting (31) and (33) into (30), the equation could be deduced as follows: It should be noted that should be transformed into the expression in absolute coordinate system - through using transformation matrices , as indicated in (16), before being added to the calculation.

The absolute velocities of the contact points on the left and right rails are written as where and are, respectively, the lateral and vertical velocities of the left and right rail cross-sections and and are, respectively, the rates of changes of the lateral and vertical irregularities on the left and right railheads.

Using (35) and (36), the velocity differences between the contact surfaces of the wheels and the rails can be written as The velocity differences , are expressed in the absolute coordinate system, shown in (37), and then need to be transformed into the expressions in the local coordinate systems - through (18).

Then, plugging (18) into (37), the velocity differences are

Equation (38) is plugged into (28) to calculate the lateral and longitudinal creepages. The calculation of the spin creepages by using the third formula in (28) neglects the effect of the wheelset axle bending since its effect on the spin creepages is quite small.

#### 3. Results and Discussions

In order to verify the present model considering the effect of the first-order bending modes of the wheelset in the two planes mentioned above on the rolling contact behavior, a harmonic track irregularity with fixed wavelength of 1.18 m is used. The excitation frequency of the irregularity is 117 Hz at the speed of 500 km/h, which is exactly equal to the resonance frequency of the wheelset first-order bending mode. The fluctuating amplitude of the irregularity is 0.05 mm. The calculation cases include the speeds of 200 km/h, 300 km/h, 400 km/h, and 500 km/h. The parameters used in the calculation are listed in Table 1 [20].

In Figure 7, the right vertical axis denotes the difference (DIFF) between the normal forces calculated by the models of the rigid (RW) and flexible (EW) wheelsets. The lines of DIFF in Figure 7 are obtained by subtracting the lines of RW from the corresponding lines of EW. The numerical computation benefits from fast double precision. Figure 7 shows that the amplitudes of the total vertical wheel/rail forces increase with the speed increase, as indicated by the red and blue broken lines. The amplitude of the normal force difference denoted by DIFF, that is, the normal force variation caused by the first-order wheelset bending deformation, is about 0.2 N when the speed is less than 500 km/h, as shown in Figures 7(a), 7(b), and 7(c). In the case of 500 km/h, as shown in Figure 7(d), the maximum amplitude of the difference reaches about 1.2 N, which is 6 times larger than those of the cases of Figures 7(a), 7(b), and 7(c). This is caused by the excited resonance of first-order bending mode. An increase in the wheel/rail normal force caused by the first-order bending resonance is quite small, compared with the static wheel/rail forces. But such a small difference of the force applied on wheels at a fixed frequency could cause a polygonal wear on high-speed wheels in a long-time service. For instance, now the high-speed vehicles in service widely use the wheelsets with diameters in the range of 840 mm–920 mm, and the operating speed of the high-speed trains is 300 km/h. Large amounts of wheel tread roughness measured from the sites show that the trilateral wear on the wheels is very common due to the wheelset first-order bending resonance excited probably.

Figure 8 indicates the PSDs (power spectral densities) of the vertical forces, corresponding to those analyzed in time domain in Figure 7. The frequencies of the dominant vibration components (at peaks) in Figure 8 are the excitation frequencies of the harmonic track irregularity, which are, respectively, 46 Hz, 69 Hz, 92 Hz and 117 Hz at the speeds of 200 km/h, 300 km/h, 400 km/h, and 500 km/h. This indicates that track irregularity plays a dominant role in the influence on wheel/rail vertical force. In other words, the influence of the wheelset flexibility is not significant if the resonance frequency of the first-order bending mode is far away from the passing irregularity frequency. The lines of DIFF illustrate the PSDs of the differences between the vertical forces applied on the rigid and flexible wheelsets. In Figures 7(a), 7(b), and 7(c), three frequency components dominate the PSDs of the force differences. The first is that Hz, which is the natural frequency of the wheel-rail contact system independent of the operating speed. The second is the passing harmonic track irregularity frequency, that is equal to, respectively, 46 Hz, 69 Hz, 92 Hz and 117 Hz, corresponding to the speeds of 200 km/h, 300 km/h, 400 km/h and 500 km/h. The third is the first-order bending resonance frequency, . At the speed of 500 km/h, the peak value, as seen in Figure 8(d), is about two orders of magnitude larger than the values at the frequency of 117 Hz at the other three speeds. So the difference in the dynamic behavior of the rigid and flexible wheelsets has been clearly understood.

The differences of the lateral wheel/rail forces and their PSDs have almost the same features. So the discussion on them is omitted in this paper.

The excited wheelset bending mode makes a contribution to the wheel/rail interaction to some extent, as known from the above analysis on the wheel/rail forces. The first-order bending resonance causes the oscillations of the wheel/rail contact points and the relative velocity of the wheel/rail contact surfaces. Figure 9 shows the changes of contact points on the rigid and flexible wheelsets with the travelling distance. In Figure 8, the vertical axis represents the coordinates of the wheel contact points in the direction of axis of the local coordinate systems -. The black heavy lines indicate the traces of the contact points on the rigid wheelset and the fine red lines denote those on the flexile wheelset, which are, respectively, indicated by RW and EW. The contact points on the rigid wheelset are closer to the nominal rolling circles, and the contact points on the flexible wheelset move toward their corresponding field sides due to the bending deformation. The first-order bending resonance of the wheelset in the plane perpendicular to the track level causes the opposite rolling angles of the left and right wheels, which are nearly symmetrical about the wheelset centre. This mode shape is similar to the downward bending deformation of the wheelset axis, as shown in Figure 3. With the increase of the speed the passing irregularity frequency gradually approaches the wheelset resonance frequency of the first-order bending mode, so that the trace of the contact points oscillates more and more seriously. When the vehicle speeds up to 500 km/h, as displayed in Figure 9(d), the first-order bending resonance of the wheelset occurs and the left and right contact points on the elastic wheelset fiercely fluctuate around the traces of the contact points on the rigid wheelset. The fluctuating amplitude reaches 0.1 *μ*m, and is five times larger than those at the speeds of 200 km/h 300 km/h and 400 km/h.

It is very interesting that the contact point traces of the rigid wheelset do not present the obvious oscillation at the passing frequency of the harmonic regularity excitation, as indicated by the black lines of RW. However, the RW lines show the position oscillations of the contact points on the rigid wheelset at a low-frequency about 22 Hz. The frequency is probably the natural frequency of the bogie in the lateral direction. The resonance of the bogie is caused by the harmonic track irregularity. When the speed is 400 km/h, the low frequency oscillation reaches the largest, as indicated in Figure 9(c), and when the speed is 500 km/h, the oscillation disappears, as indicated by the coarse black lines in Figure 9(d). The further discussion on dynamic behavior of the bogie does not belong to the task of this paper.

The above analysis shows the evident impact of the wheelset bending deformation on the contact points. The wheel/rail contact forces rely on not only the contact point positions but also the relative speeds of the contact surfaces represented by the creepages. There is a small difference between the creepages calculated with the models of the elastic and rigid wheelsets, as indicated in Figure 10.

As well known, rail corrugation and wheel polygonization are mainly related to the sensitive resonances of the wheel/rail system [1, 21, 22]. In order to verify that the first-order bending resonance of the wheelset is active and find the sensitive resonance frequencies of the wheel/rail systems, the present model of the high-speed vehicle coupled with the track is used to analyze the dynamic behaviour of the wheel/rail systems under a white-noise excitation and at different speeds. The analyzed frequency range is from 0 to 500 Hz and the operation speed is from 200 km/h to 400 km/h. Figure 11 shows the wheel/rail force responses analyzed in frequency-domain. Figures 11(a) and 11(b) illustrate the linear spectra of the lateral and vertical wheel/rail forces of the rigid wheelset, respectively, and Figures 11(c) and 11(d) indicate the corresponding results of the flexible wheelset. In Figure 11(b), indicates the resonance frequency of the wheel/rail system, as mentioned above. It is about 31 Hz, and independent of the operation speed. In Figures 11(a) and 11(b), is the passing frequency of the sleepers and is the twice passing frequency of the sleepers or called harmonic frequency due to the excitation of the multisleepers with an equal pitch between the neighboring sleepers. In Figures 11(c) and 11(d), denotes the bending resonance frequency of the wheelset. equals 117 Hz, and it is evident that the resonance frequency is independent of the vehicle speed. In all the figures of Figure 11, always appears at about 380 Hz. It is one of the resonance frequencies of the system consisting of the track and the bogie of the vehicle. It will be further analyzed in the following.

Carefully making a comparison of Figures 11(a), 11(b), 11(c), and 11(d) shows the interesting results, in which one should carefully pay attention to the vertical coordinate scales of the figures. Under the present excitation of the white-noise, the wheel/rail force response has high peaks at the passing frequencies of the sleepers, or the wheel/rail force response to the excitation of the periodically discrete sleepers is clearly characterized by using the existing model of the coupling vehicle/track without considering the effect of the first-order bending resonance of wheelset, as indicated in Figures 11(a) and 11(b). But the wheel/rail force response has very low peaks at the passing frequencies of the sleepers in Figures 11(c) and 11(d) in which the results are calculated by using the model considering the effect of the flexible wheelset. Actually, the wheel/rail force response to the excitation of the periodically discrete sleepers is overwhelmed by the responses at the other resonance frequencies. So the first-order bending resonance of the wheelset plays a dominant role in the wheel/rail force response under the present excitation of the white-noise and in the analyzed speed range.

The peaks at Hz are always remarkable in the results obtained by the two models, or the vibration at about 380 Hz always makes a great contribution to the wheel/rail force response. This is because is the resonance frequency of the track including the effect of the sprung mass and the bogie frame of the vehicle, and in the employed two models the track modelling is the same. In order to find the mode shape of the track resonance, a finite element model for the track with a bogie is developed and used to calculate the resonance frequencies and the corresponding mode shapes of them. In Figure 12, the two figures indicate the mode shapes of the track with a bogie, the corresponding frequencies of which are, respectively, 374 Hz and 390 Hz. They are close to 380 Hz.

#### 4. Conclusions

This paper presents a new vehicle/track coupling dynamics model considering the effect of the first-order bending resonance of the flexible wheelset. A new wheel/rail contact geometry model is derived in detail and considers the interaction of a flexible wheelset and a pair of rails. So the existing coupling dynamics model of railway vehicle/track is further improved. In order to verify the new model and clarify the influence of the first-order bending resonance of the flexible wheelset on the wheel/rail contact behavior, the new and existing models are used to analyze the difference of the wheel/rail contact behavior under the excitation of a harmonic track irregularity with a fixed wavelength and a white-noise at different operating speeds. The results calculated by the two models are compared to understand the effect of the wheelset first-order bending on the rolling contact behavior of the wheel/rail system. The conclusions are given as follows.(1)In the analysis cases that indicate the different operating speeds and a harmonic irregularity and white-noise excitations, the first-order bending mode of the wheelset in service can be excited easily. When the passing irregularity frequency equals the natural frequency of the first-order bending mode of wheelset, the resonance of the wheel/rail systems occurs. In such a situation, the wheelset bending has a great influence on the wheel/rail rolling contact behavior.(2)When the first-order bending mode in the vertical plane is excited at the speed of 500 km/h, the wheel/rail forces and contact points oscillate more seriously than those at the other speeds at the first-order bending resonance frequency. Besides, the resonance easily occurs between the flexibility wheelset and the rails when the passing frequency of the track irregularity excitation is equal or close to the wheelset bending resonance frequency at any speed. Though the numeric differences between the dynamic responses of the rigid and flexible wheelsets are small, they could contribute to the development of the polygonal wear on the wheel treads after a long-time service of the wheelset.(3)Developing the coupling dynamics model of railway vehicle and track considering the effect of structural flexibility of the key parts is a complex subject. The present work is just at a preliminary stage. The next step will be to consider modeling the high-order bending resonances of wheelset, wheel flexibility and powered wheelset flexibility in the wheel/rail rolling contact model.

#### Nomenclature

(=1–4): | The number of wheelsets |

(=1-2): | The number of bogies |

: | The components in - and -directions of the spin moment on the left wheel |

: | The components in - and -directions of the spin moment on the right wheel |

: | Components in -, -, and -directions of wheel/rail normal force on left wheel |

: | The components in -, -, and -directions of wheel/rail normal force on the right wheel |

: | The components in -, -, and -directions of wheel/rail tangent force on the left wheel |

: | The components in -, -, and -directions of wheel/rail tangent force on the right wheel |

: | The components in -, -, and -directions of primary suspension forces on the left sides of a wheelset |

: | The components in -, -, and -directions of primary suspension forces on the right sides of a wheelset |

: | Time |

: | The displacement of the bending deformation in the directions of -axis |

: | The displacement of the bending deformation in the directions of -axis |

: | Young’s modulus of the wheelset material |

: | The wheelset density |

: | The equivalent cross-sectional area of the axle |

: | The gravity of the wheelset |

: | The inertia moments of the cross-section with respect to the -axis and the -axis, respectively |

: | The concentrated forces on the wheelset axle in the - and -directions |

: | The concentrated moments on the wheelset axle, respectively, in the negative -axis direction and the -axis direction |

: | The coordinates describing the distances from the wheelset axle to the contact points of the left and right wheels/rails, respectively |

: | The coordinates describing the distances from the axle to the left and right primary suspensions, respectively |

: | The th normalized modal functions describing the bending deformations of the axle in the planes parallel and perpendicular to the track level plane, respectively |

: | The th generalized coordinates of the bending deformations in the directions of the -axis and the -axis, respectively |

: | The second-order derivatives of and with respect to time |

: | The generalized coordinates corresponding to the flexible modes |

: | The generalized coordinates corresponding to the rigid modes |

: | The th circular frequencies |

: | The flexible modal functions normalized in the directions of the -axis and the -axis, respectively |

: | The rigid modal functions normalized in the directions of the -axis and the -axis, respectively |

: | The total number of the modes in the calculation |

: | The total number of the rigid modal functions |

: | The total number of the flexible modal functions |

: | The half distance between the left and right action points of the primary suspension forces on the wheelset |

: | The half distance between the left and right action points of the primary suspension forces on the bogie |

: | The half of the length of the wheelset axle |

: | The yawing angle of the bogie |

: | The pitching angles of the bogie |

: | The rolling angles of the bogie |

: | The primary suspension stiffness in -, -, and -directions, respectively |

: | The primary suspension damping in -, -, and -directions, respectively |

: | The vertical distance between the mass centers of the bogie and the wheelset |

: | The left and right wheel/rail contact angles |

: | The rolling angles of the left and right wheels |

: | The yawing angles of the left and right wheels |

: | The rail cant |

: | The left and right rail rolling angles |

: | The abscissa ordinate of the instant rolling circle center |

: | The cosine components of the axle centerline in the -, -, and -directions, respectively, of a rigid wheelset |

: | The cosine components of the axle centerline in the -, -, and -directions, respectively, of left and right wheels of a flexible wheelset |

: | The rolling angle of a rigid wheelset |

: | The yawing angle of a rigid wheelset |

: | The coordinates of a rigid wheelset center |

: | The relative vertical displacements of the left and right wheels of the flexible wheelset, respectively |

: | The initial vertical compressions of the left and right wheels at |

: | The vertical displacements of the left and right wheels caused by the wheelset bending deformation occurred in the transverse section of the track |

: | The vertical irregularities of the left and right rails |

: | The left and right wheel/rail normal compressions |

: | The left and right wheel/rail normal forces |

: | The longitudinal, lateral, and spin creepages |

: | Forward velocities for left and right wheels |

: | The relative longitudinal and lateral velocities |

: | The relative angular velocities |

: | The absolute velocity |

: | A relative velocity |

: | A transport velocity |

: | The transport velocity used in creepage calculation of the flexible wheelset |

: | The relative velocities used in creepage calculation of the flexible wheelset |

: | The angular velocity vectors of the left and right wheels |

: | The velocity differences of the contact points on left and right rail and wheel, respectively. |

#### Acknowledgments

The present work is supported by the National Natural Science Foundation of China (U1134202), the National Basic Research Program of China (2011CB711103), and the Program for Changjiang Scholars and Innovative Research Team in university (IRT1178 and SWJTU12ZT01).

#### References

- X. S. Jin, L. Wu, J. Y. Fang, S. Q. Zhong, and L. Ling, “An investigation into the mechanism of the polygonal wear of metro train wheels and its effect on the dynamic behavior of a wheel/ rail system,”
*Vehicle System Dynamics*, vol. 50, no. 12, pp. 1817–1834, 2012. - K. L. Knothe and S. L. Grassie, “Modelling of railway track and vehicle/track interaction at high frequencies,”
*Vehicle System Dynamics*, vol. 22, no. 3-4, pp. 209–262, 1993. View at Scopus - T. B. Meinders, P. Meinke, and W. O. Schiehlen, “Wear estimation in flexible multibody systems with application to railroads,” in
*Proceedings of the ECCOMAS Thematic Conference on Multi-Body Dynamics*, pp. 21–24, Madrid, Spain, June 2005. - T. Szolc, “Medium frequency dynamic investigation of the railway wheelset-track system using a discrete-continuous model,”
*Archive of Applied Mechanics*, vol. 68, no. 1, pp. 30–45, 1998. View at Scopus - K. Popp, H. Kruse, and I. Kaiser, “Vehicle-track dynamics in the mid-frequency range,”
*Vehicle System Dynamics*, vol. 31, no. 5-6, pp. 423–464, 1999. View at Scopus - N. Chaar,
*Wheelset structural flexibility and track flexibility in vehicle/track dynamic interaction [Ph.D. thesis]*, Royal Institute of Technology, Stockholm, Sweden, 2007. - M. Meywerk, “Polygonalization of railway wheels,”
*Archive of Applied Mechanics*, vol. 69, no. 2, pp. 105–120, 1999. View at Publisher · View at Google Scholar · View at Scopus - T. Szolc, “Medium frequency dynamic analysis of various designs of the high-speed railway bogies,” in
*Proceedings of the 7th Mini Conference on Vehicle System Dynamics, Identification and Anomalies*, pp. 255–262, November 2000. View at Scopus - T. Meinders, “Modeling of a railway wheelset as a rotating elastic multibody system,”
*Machine Dynamics Problems*, vol. 20, pp. 209–219, 1998. - T. Meinders and P. Meinker, “Rotor dynamics and irregular wear of elastic wheelsets,” in
*System Dynamics and Long-Term Behaviour of Railway Vehicles, Track and Subgrade*, vol. 6 of*Lecture Notes in Applied Mechanics*, pp. 133–152, 2002. - J. Fayos, L. Baeza, F. D. Denia, and J. E. Tarancón, “An Eulerian coordinate-based method for analysing the structural vibrations of a solid of revolution rotating about its main axis,”
*Journal of Sound and Vibration*, vol. 306, no. 3–5, pp. 618–635, 2007. View at Publisher · View at Google Scholar · View at Scopus - L. Baeza, P. Vila, A. Roda, and J. Fayos, “Prediction of corrugation in rails using a non-stationary wheel-rail contact model,”
*Wear*, vol. 265, no. 9-10, pp. 1156–1162, 2008. View at Publisher · View at Google Scholar · View at Scopus - L. Baeza, P. Vila, G. Xie, and S. D. Iwnicki, “Prediction of rail corrugation using a rotating flexible wheelset coupled with a flexible track model and a non-Hertzian/non-steady contact model,”
*Journal of Sound and Vibration*, vol. 330, no. 18-19, pp. 4493–4507, 2011. View at Publisher · View at Google Scholar · View at Scopus - A. A. Shabana, “Flexible multibody dynamics: review of past and recent developments,”
*Multibody System Dynamics*, vol. 1, no. 2, pp. 189–222, 1997. View at Scopus - X. Jin, P. Wu, and Z. Wen, “Effects of structure elastic deformations of wheelset and track on creep forces of wheel/rail in rolling contact,”
*Wear*, vol. 253, no. 1-2, pp. 247–256, 2002. View at Publisher · View at Google Scholar · View at Scopus - X. Xiao, X. Jin, Z. Wen, M. Zhu, and W. Zhang, “Effect of tangent track buckle on vehicle derailment,”
*Multibody System Dynamics*, vol. 25, no. 1, pp. 1–41, 2011. View at Publisher · View at Google Scholar · View at Scopus - X. Xiao, X. Jin, and Z. Wen, “Effect of disabled fastening systems and ballast on vehicle derailment,”
*Journal of Vibration and Acoustics*, vol. 129, no. 2, pp. 217–229, 2007. View at Publisher · View at Google Scholar · View at Scopus - X. Xiao, X. Jin, Y. Deng, and Z. Zhou, “Effect of curved track support failure on vehicle derailment,”
*Vehicle System Dynamics*, vol. 46, no. 11, pp. 1029–1059, 2008. View at Publisher · View at Google Scholar · View at Scopus - W. M. Zhai,
*Vehicle/Track Coupling Dynamics*, Chinese Science Press, 3rd edition, 2007. - S. G. Zhang,
*Gaosu Lieche Sheji Fangfa Yanjiu*, Study on Design Method of High-Speed Trains, China Railway Publishing House, Beijing, China, 2009. - W. Li, X. Li, L. Wu et al., “Investigation into mechanism of metro rail corrugation using experimental and theoretical methods,” in
*Proceedings of the 22th Symposium of the International Association of Vehicle System Dynamics (IAVSD '11)*, Manchester, UK, 2011. - X. S. Jin, B. Q. Wu, L. Wu, J. Y. Fang, and S. Q. Zhong, “An investigation into the mechanism of the polygonal wear of metro train wheels and its effect on the dynamical behavior of wheel/rail system,” in
*Proceedings of the 22th Symposium of the International Association of Vehicle System Dynamics (IAVSD '11)*, Manchester, UK, 2011.