Dynamic responses of a high-speed railway car due to wheel polygonalisation

The polygonal wear around the wheel circumference could pose highly adverse influences on the wheel/rail interactions and thereby the performance of the vehicle system. In this study, the effects of wheel polygonalisation on the dynamic responses of a high-speed rail vehicle are investigated through development and simulations of a comprehensive coupled vehicle/track dynamic model. The model integrates flexible slab track, wheelsets and axle boxes subsystem models so as to account for elastic deformations caused by impact loads induced by the wheel polygonalisation. A field-test programme was undertaken to acquire the polygonal wear profile and axle box acceleration response of a high-speed train, and the data are used to demonstrate the validity of the coupled vehicle/track system model. Subsequently, the effects of wheel polygonalisation are evaluated in terms of wheel/rail impact forces, axle box vertical acceleration and dynamic stress developed in the axle considering different amplitudes and harmonic orders of the polygonal wear. The results suggest that the high-order wheel polygonalisation can give rise to high-frequency impact loads at the wheel/rail interface, and excite some of the vibration modes of the wheelset and the axle box leading to high-magnitude axle box acceleration and dynamic stress in the wheelset axle.     

High frequency railway vehicle-track dynamics through flexible rotating wheelsets

Some railway problems, such as the corrugation of rails or the impact caused by a wheelflat, are associated with a vehicle–track coupled dynamic phenomenon. Models for the analysis of these problems must account for the structural vibrations of the track components (rails and sleepers), but the most adequate approach for the wheelset has not been sufficiently investigated until present. The wheelset can be considered as an undeformable solid, as an elastic structure where the rotation effects are neglected, or as a rotating flexible solid. In order to fill this gap, this article presents a methodology to use the structural vibrations of a rotating wheelset in high-frequency railway dynamics analyses. The model makes use of Eulerian modal coordinates, a formulation that provides very low computational cost. The method is applied in this article to a wheelflat impact calculation and a vehicle running on a corrugated track. The results show the importance of the more realistic model in the simulations, mainly in certain frequencies.     

Wheel/Rail Contact: Geometrical Study*

Before trying to ascertain the precise nature of the Wheel-Rail contact (internal stresses and strains, wear and friction, security against a derailment, dynamic behaviour of the vehicle, etc.) the geometrical problem must necessarily be solved. That is, for each position of the wheelset (this is defined by 6 parameters, of which only 4 are independent) the two dependent parameters and the coordinates of the points of contact of each wheel, and rail must be obtained. A new method is proposed of obtaining the spatial position of a wheelset with reference to the rails, from the most general point of view.     

Method for obtaining the wheel–rail contact location and its application to the normal problem calculation through ‘CONTACT’

This work presents a robust methodology for calculating inter-penetration areas between railway wheel and rail surfaces, the profiles of which are defined by a series of points. The method allows general three-dimensional displacements of the wheelset to be considered, and its characteristics make it especially suitable for dynamic simulations where the wheel–rail contact is assumed to be flexible. The technique is based on the discretisation of the geometries of the surfaces in contact, considering the wheel as a set of truncated cones and the rail as points. By means of this approach, it is possible to reduce the problem to the calculation of the intersections between cones and lines, the solution for which has a closed-form expression. The method has been used in conjunction with the CONTACT algorithm in order to solve the static normal contact problem when the lateral displacement of the wheelset, its yaw angle and the vertical force applied in the wheelset centroid are prescribed. The results consist of smooth functions when the dependent coordinates are represented as a function of the independent ones, lacking the jump discontinuities that are present when a rigid contact model is adopted. Example results are shown and assessed for the normal contact problem for different lateral and yaw positions of the wheelset on the track.     

Out-of-round railway wheels—a study of wheel polygonalization through simulation of three-dimensional wheel–rail interaction and wear

The objective of this study is to develop a tool for investigation of wheel tread polygonalization with radial irregularities including 1 to 20 wavelengths around the circumference of the wheel. Therefore, an existing multibody system model for simulation of general three-dimensional train–track interaction (accounting for frequencies up to several kHz) is extended with rolling contact mechanics according to FASTSIM. Furthermore, the model is also modified to allow for general wheel–rail profiles. The numerical model uses the concept of an iteration scheme including simulation of dynamic train–track interaction in the time domain coupled with a long-term wear model. A demonstration example including a bogie of a subway train travelling on a straight track is presented. In the example, an initial wheel out-of-roundness (OOR) is applied to the wheels. This irregularity is based on an amplitude spectrum derived from measurements on new wheels. Simulation results show that the most important wavelength-fixing mechanisms of the wheel OOR are (i) the vertical resonance of the coupled train–track system at approximately 40 Hz (the P2 resonance) and (ii) the frequency region including the lowest vertical track antiresonance at 165 Hz, where the dynamic track stiffness is high. Only a straight track is studied, but the model allows for asymmetric train motion on such a track.     

Influence of polygonal wear of railway wheels on the wheel set axle stress

The coupled vehicle/track dynamic model with the flexible wheel set was developed to investigate the effects of polygonal wear on the dynamic stresses of the wheel set axle. In the model, the railway vehicle was modelled by the rigid multibody dynamics. The wheel set was established by the finite element method to analyse the high-frequency oscillation and dynamic stress of wheel set axle induced by the polygonal wear based on the modal stress recovery method. The slab track model was taken into account in which the rail was described by the Timoshenko beam and the three-dimensional solid finite element was employed to establish the concrete slab. Furthermore, the modal superposition method was adopted to calculate the dynamic response of the track. The wheel/rail normal forces and the tangent forces were, respectively, determined by the Hertz nonlinear contact theory and the Shen–Hedrick–Elkins model. Using the coupled vehicle/track dynamic model, the dynamic stresses of wheel set axle with consideration of the ideal polygonal wear and measured polygonal wear were investigated. The results show that the amplitude of wheel/rail normal forces and the dynamic stress of wheel set axle increase as the vehicle speeds rise. Moreover, the impact loads induced by the polygonal wear could excite the resonance of wheel set axle. In the resonance region, the amplitude of the dynamic stress for the wheel set axle would increase considerably comparing with the normal conditions.     

Development and validation of a wear model for the analysis of the wheel profile evolution in railway vehicles

The numerical wheel wear prediction in railway applications is of great importance for different aspects, such as the safety against vehicle instability and derailment, the planning of wheelset maintenance interventions and the design of an optimal wheel profile from the wear point of view. For these reasons, this paper presents a complete model aimed at the evaluation of the wheel wear and the wheel profile evolution by means of dynamic simulations, organised in two parts which interact with each other mutually: a vehicle's dynamic model and a model for the wear estimation. The first is a 3D multibody model of a railway vehicle implemented in SIMPACK?, a commercial software for the analysis of mechanical systems, where the wheel–rail interaction is entrusted to a C/C++user routine external to SIMPACK, in which the global contact model is implemented. In this regard, the research on the contact points between the wheel and the rail is based on an innovative algorithm developed by the authors in previous works, while normal and tangential forces in the contact patches are calculated according to Hertz's theory and Kalker's global theory, respectively. Due to the numerical efficiency of the global contact model, the multibody vehicle and the contact model interact directly online during the dynamic simulations.

The second is the wear model, written in the MATLAB® environment, mainly based on an experimental relationship between the frictional power developed at the wheel–rail interface and the amount of material removed by wear. Starting from a few outputs of the multibody simulations (position of contact points, contact forces and rigid creepages), it evaluates the local variables, such as the contact pressures and local creepages, using a local contact model (Kalker's FASTSIM algorithm). These data are then passed to another subsystem which evaluates, by means of the considered experimental relationship, both the material to be removed and its distribution along the wheel profile, obtaining the correspondent worn wheel geometry.

The wheel wear evolution is reproduced by dividing the overall chosen mileage to be simulated in discrete spatial steps: at each step, the dynamic simulations are performed by means of the 3D multibody model keeping the wheel profile constant, while the wheel geometry is updated through the wear model only at the end of the discrete step. Thus, the two parts of the whole model work alternately until the completion of the whole established mileage. Clearly, the choice of an appropriate step length is one of the most important aspects of the procedure and it directly affects the result accuracy and the required computational time to complete the analysis.

The whole model has been validated using experimental data relative to tests performed with the ALn 501 ‘Minuetto’ vehicle in service on the Aosta–Pre Saint Didier track; this work has been carried out thanks to a collaboration with Trenitalia S.p.A and Rete Ferroviaria Italiana, which have provided the necessary technical data and experimental results.     

Investigation of the effects of sleeper-passing impacts on the high-speed train

The sleeper-passing impact has always been considered negligible in normal conditions, while the experimental data obtained from a High-speed train in a cold weather expressed significant sleeper-passing impacts on the axle box, bogie frame and car body. Therefore, in this study, a vertical coupled vehicle/track dynamic model was developed to investigate the sleeper-passing impacts and its effects on the dynamic performance of the high-speed train. In the model, the dynamic model of vehicle is established with 10 degrees of freedom. The track model is formulated with two rails supported on the discrete supports through the finite element method. The contact forces between the wheel and rail are estimated using the non-linear Hertz contact theory. The parametric studies are conducted to analyse effects of both the vehicle speeds and the discrete support stiffness on the sleeper-passing impacts. The results show that the sleeper-passing impacts become extremely significant with the increased support stiffness of track, especially when the frequencies of sleeper-passing impacts approach to the resonance frequencies of wheel/track system. The damping of primary suspension can effectively lower the magnitude of impacts in the resonance speed ranges, but has little effect on other speed ranges. Finally, a more comprehensively coupled vehicle/track dynamic model integrating with a flexible wheel set is developed to discuss the sleeper-passing-induced flexible vibration of wheel set.     

Carbody elastic vibrations of high-speed vehicles caused by bogie hunting instability

In particular locations of the high-speed track, the worn wheel profile matched up with the worn rail profile will lead to an extremely high-conicity wheel–rail contact. Consequently, the bogie hunting instability arises, which further results in the so-called carbody shaking phenomenon. In this paper, the carbody elastic vibrations of a high-speed vehicle in service are firstly introduced. Modal tests are conducted to identity the elastic modes of the carbody. The ride comfort and running safety indices for the tested vehicle are evaluated. The rigid–flexible coupling dynamic model for the high-speed passenger car is then developed by using the FE and MBS coupling approach. The rail profiles in those particular locations are measured and further integrated into the simulation model to reproduce the bogie hunting and carbody elastic vibrations. The effects of wheel and rail wear on the vehicle system response, e.g. wheelset bifurcation graph and carbody vibrations, are studied. Two improvement measures, including the wheel profile modification and rail grinding, are proposed to provide possible solutions. It is found that the wheel–rail contact conicity can be lowered by decreasing wheel flange thickness or grinding rail corner, which is expected to improve the bogie hunting stability under worn rail and worn wheel conditions. The carbody elastic vibrations caused by bogie hunting instability can be further restrained.     

A new method for determining wheel–rail multi-point contact

A new method for wheel–rail multi-point contact is presented in this paper. In this method, the first- and the second-order derivatives of the wheel–rail interpolation distance function and the elastic wheel–rail virtual penetration are used to determine multiple contact points. The method takes account of the yaw angle of the wheelset and allows the identification of all possible points of contact between wheel and rail surfaces with an arbitrary geometry. Static contact geometry calculations are first carried out using the developed method for both new and worn wheel profiles and with a new rail profile. The validity of the method is then verified by simulations of a coupled vehicle and track system dynamics over a small radius curve. The simulation results show that the developed method for multi-point contact is efficient and reliable enough to be implemented online for simulations of vehicle–track system dynamics.     

An investigation on wheel/rail impact dynamics with a three-dimensional flat model

Wheel flat is one kind of railway train wheelset defects. It has great influence on wheel/rail dynamics and damages. In most of the presented studies, wheel/rail impact velocity or rolling radius variation of the wheel because of flat spot was taken into account to study the wheel/rail impact dynamics. In this paper, a three-dimensional wheel flat model considering the length, width and depth of the flat spot is established. Including the wheel rotation and wheel/rail contact geometry, a high-speed vehicle–rail coupling system dynamics model is developed to investigate the effect of the wheel flat on the wheel/rail dynamics. With time integration method of the models, the impact dynamics of the wheel/rail system with three types of flat width and five kinds of flat length are obtained. The results show that the width, the length of the wheel flat and the width/length ratio have a great influence on the wheel/rail impact dynamics. The wheel/rail impact dynamics of the flat with large width is more severe than with small width as the flat length is fixed. When a flat spot occurs, the permissible length of the wheel defect, needless to action, is 25?mm in maximum. The speed safety domain with three kinds of flat width/length ratio of a vehicle is obtained according to the wheel/rail vertical force limitation.     

Out-of-round railway wheels—a study of wheel polygonalization through simulation of three-dimensional wheel-rail interaction and wear

The objective of this study is to develop a tool for investigation of wheel tread polygonalization with radial irregularities including 1 to 20 wavelengths around the circumference of the wheel. Therefore, an existing multibody system model for simulation of general three-dimensional train-track interaction (accounting for frequencies up to several kHz) is extended with rolling contact mechanics according to FASTSIM. Furthermore, the model is also modified to allow for general wheel-rail profiles. The numerical model uses the concept of an iteration scheme including simulation of dynamic train-track interaction in the time domain coupled with a long-term wear model. A demonstration example including a bogie of a subway train travelling on a straight track is presented. In the example, an initial wheel out-of-roundness (OOR) is applied to the wheels. This irregularity is based on an amplitude spectrum derived from measurements on new wheels. Simulation results show that the most important wavelength-fixing mechanisms of the wheel OOR are (i) the vertical resonance of the coupled train-track system at approximately 40 Hz (the P2 resonance) and (ii) the frequency region including the lowest vertical track antiresonance at 165 Hz, where the dynamic track stiffness is high. Only a straight track is studied, but the model allows for asymmetric train motion on such a track.     

Impact Loads due to Wheel Flats and Shells

A Finite Element (FE) model of vehicle-track system is employed to duplicate the experiments carried out by British Rail and CP Rail System. The theoretical results of the wheel/rail contact forces, rail-pad forces and strains in the rail showed very good correlation to the experimental data. Extensive results are compared with experimental data in the time domain for through validation of the developed model. The characteristics of the impact loads due to wheel flats and shells are investigated based on the validated FE model. The study shows that the shape and size of flat or shell, axle load, vehicle speed and rail-pad stiffness mainly affect the impact loads. Adding elastomeric shear pads on the wheelset bearing does not reduce the wheel/rail dynamic contact force but it may reduce the dynamic force on the bearing. Reducing rail-pad stiffness to a certain level on a concrete-tie track may significantly reduce the dynamic load and the force transmitted to the concrete tie.     

Near surface deformation of rails and wheels in a curve and during acceleration and braking

A hierarchical concept for describing the rolling/sliding contact of a wheel on a rail is proposed, which is developed to model the near surface deformation of the contacting partners depending on the vehicle system dynamics. In principal the vehicle dynamics simulation delivers the input for a static finite element model for the elastic-plastic rolling/sliding contact investigation. To include also more realistic rough surfaces in the calculations, which is absolutely necessary to predict the highly deformed surface layer as observed in wheels and rails, a further refinement of the mesh is necessary and performed. The output values of the static finite element model are now boundary conditions for the finite element micro-model. The results of this micro-model correspond well with the micro-structural investigations of wheel and rail surfaces.     

Research into the problem of wheel tread spalling caused by wheelset longitudinal vibration

This study mainly focuses on the mechanism of wheel tread spalling through wheelset longitudinal vibration that has been often neglected. Analysis of two actual cases of the wheel tread spalling problem leads to the conclusion that the wheel tread spalling is closely related to the wheelset longitudinal vibration in some locomotives, and many of these problems can be reasonably explained if the wheelset longitudinal vibration is considered. For better understanding of some abnormal wheel spalling problems, the formations of the wheelset longitudinal vibration and the wheel/rail contact parameters were analysed in the initial wheel tread spalling. With the preliminary analytical results, the wheelset longitudinal dynamic behaviour, the characteristics of wheel/rail contact and the mechanics in the condition of the wheelset longitudinal vibration were further studied quantitatively. The results showed that the wheelset longitudinal vibration changed not only the limit of these parameters and the position of principal stress, but also the direction of the principal stress on the surface of wheel/rail contact patch. It is likely that the significant stress changes provoke too much stress on the surface of wheel/rail contact patch, cause fatigue in wheel/rail contact patch and eventually lead to wheel tread spalling. The results of these studies suggest that the suppression of the wheelset longitudinal vibration extends wheel/rail life and the addition of a vertical damper with an ahead angle provides a possible solution to the wheel spalling problem.     

A method to determine the two-point contact zone and transfer of wheel–rail forces in a turnout

A practical method to determine the zone of two contact points and the transfer of wheel–rail forces between two rails in a turnout is presented in this paper. The method is based on a wheel–rail elastic penetration assumption and used to study a turnout system for a 200 km/h high-speed railway in China. Rail profiles in a number of key sections in the turnout are identified first, and profiles in other sections are then obtained by interpolation between key sections. The track is modelled as flexible with rails and sleepers represented by beams and the interaction between the vehicle and turnout is simulated for cases of the vehicle passing the turnout. Results are mainly presented for two-point contact positions and the characteristics of the wheel–rail forces transference. It is found that the heights of the switch and crossing rail top have significant effects on the wheel–rail contact forces. Finally, the optimised top height for the crossing rails is proposed to reduce the system dynamic force in the turnout system.     

Numerical investigation on effect of the relative motion of stock/switch rails on the load transfer distribution along the switch panel in high-speed railway turnout

In railway turnout, the stock rail and switch rail are separated to enable the vehicle changing among the tracks, and they are provided with different rail resilience level on the baseplate. Therefore, there will be vertical relative motion between stock/switch rails under the wheel loads, and the relative motion will affect consequentially the wheel–rail contact conditions. A method is developed to investigate the effect of the relative motion of stock/switch rails on the load transfer distribution along the switch panel in high-speed railway turnout. First, the rigid wheel–rail contact points of stock/switch rails are calculated based on the trace line method, and then the contact status is determined by the presented equations, finally, the distribution of wheel–rail contact forces of stock/switch rails is obtained based on the continuity of interface displacements and forces which using an approximate surface deformation method. Some parametric studies have been performed, such as the lateral displacement of wheel set, the vertical contact forces, the wheel profiles and the vertical stiffness of rail pad. The results of the parametric study are presented and discussed.     

Railway ground vibrations induced by wheel and rail singular defects

Railway local irregularities are a growing source of ground-borne vibration and can cause negative environmental impacts, particularly in urban areas. Therefore, this paper analyses the effect of railway track singular defects (discontinuities) on ground vibration generation and propagation. A vehicle/track/soil numerical railway model is presented, capable of accurately predicting vibration levels. The prediction model is composed of a multibody vehicle model, a flexible track model and a finite/infinite element soil model. Firstly, analysis is undertaken to assess the ability of wheel/rail contact models to accurately simulate the force generation at the wheel/rail contact, in the presence of a singular defect. It is found that, although linear contact models are sufficient for modelling ground vibration on smooth tracks, when singular defects are present higher accuracy wheel/rail models are required. Furthermore, it is found that the variation in wheel/rail force during the singular defect contact depends on the track flexibility, and thus requires a fully coupled vehicle/track/foundation model. Next, a parametric study of ground vibrations generated by singular rail and wheel defects is undertaken. Six shapes of discontinuity are modelled, representing various defect types such as transition zones, switches, crossings, rail joints and wheel flats. The vehicle is modelled as an AM96 train set and it is found that ground vibration levels are highly sensitive to defect height, length and shape.     

Inverse identification of wheel–rail contact forces based on observation of wheel disc strains: an evaluation of three numerical algorithms

In this paper, three numerical algorithms for the identification of wheel–rail contact forces based on measured wheel disc strains on an instrumented railway wheelset are discussed and compared. The three algorithms include one approach resting on static calibration, one that is applying a Kalman filter and the third is exploiting an inverse identification scheme. To demonstrate and evaluate the alternative methods, two load cases including periodic excitation by sinusoidal wheel–rail irregularities and transient excitation by an insulated rail joint are considered. Based on a previously presented vehicle–track interaction model in the time domain, load scenarios are defined by taking the calculated vertical wheel–rail contact forces as the reference force to be re-identified by the proposed algorithms. The reference contact forces are applied on a finite element model of the wheel to generate synthetic observation data, that is, radial strains at the positions of the strain gauges, serving as input to the identification procedures. It is concluded that the inverse identification scheme leads to superior accuracy at higher computational cost. If on-line implementation and evaluation is required, the Kalman filter generates better accuracy than the static calibration approach.