The Stationary Motion of a One-Axle Vehicle Along a Circular Curve with Real Rail and Wheel Profiles

In this paper, we present a theory on the stationary motion of a one-axle railway vehicle along a circular curve in the presence of single- or double-point contact. The rail and the wheel profiles may be either stylized or real and as an example we take the profile combination UIC60 1:40 S1002. The mathematical model of the system is based on De Pater's first-order theory [1]. The geometrical contact problem between wheel and rail is solved by using a modified Newton-Raphson procedure. Both the cases with and without friction are considered. When friction is present, the non-linear Kalker creep law [6, 7] is used to describe the physical contact. For various values of the friction coefficient, the cant angle and the curvature of the track, the contact forces are presented as functions of the velocity parameter C v = V 2 / V 2 eq , where V is the velocity of the vehicle and V eq is the equilibrium velocity of the frictionless case. For the case of stylized profiles in which both the wheel treads and the wheel flanges are conical, and the rail cross sections are circular, we have determined the velocity range with single point contact in dependence on the friction coefficient, the conicity of the tread, the curvature of the track and the cant angle.     

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.     

Determination of the wheel rail contact patch in semi-Hertzian conditions

This paper deals with the subject of the semi-Hertzian contact, which is a way to represent the wheel rail contact in railways or roller bearing applications. The method is based on the interpenetration of the two underformed bodies' profiles. The first step deals with the problem of the shape ratio; it is proposed to compensate the two main curvatures to obtain the good ratio in Hertzian conditions. Then, Hertz and Kalker's equations are used to establish the stresses at the level of a strip. These stresses expressions are used directly in a contact model discretized in strips and tabulated as a function of the lateral displacement between the wheel and the rail. The validation is made by comparison to the previous multi-Hertzian model of the VOCOLIN software. A first test shows identical results in Hertzian conditions; a second one shows only a small difference in semi-Hertzian conditions like S1002/UIC60 1:40.     

Determination of the wheel rail contact patch in semi-Hertzian conditions

This paper deals with the subject of the semi-Hertzian contact, which is a way to represent the wheel rail contact in railways or roller bearing applications. The method is based on the interpenetration of the two underformed bodies' profiles. The first step deals with the problem of the shape ratio; it is proposed to compensate the two main curvatures to obtain the good ratio in Hertzian conditions. Then, Hertz and Kalker's equations are used to establish the stresses at the level of a strip. These stresses expressions are used directly in a contact model discretized in strips and tabulated as a function of the lateral displacement between the wheel and the rail. The validation is made by comparison to the previous multi-Hertzian model of the VOCOLIN software. A first test shows identical results in Hertzian conditions; a second one shows only a small difference in semi-Hertzian conditions like S1002/UIC60 1:40.     

Design of an optimised wheel profile for rail vehicles operating on two-track gauges

This article sets out an optimum synthesis methodology for wheel profiles of railway vehicles in order to secure good dynamic behaviour with different track configurations. Specifically, the optimisation process has been applied to the case of rail wheelsets mounted on double-gauge bogies that move over two different gauges, which also have different types of rail: the Iberian gauge (1668 mm) and the International Union of Railways (UIC) gauge (1435 mm). Optimisation is performed using Genetic Algorithms and traditional optimisation methods in a complementary way. The objective function used is based on an ideal equivalent conicity curve which ensures good stability on straight sections and also proper negotiation of curves. To this end, the curve is constructed in such a way that it is constant with a low value for small lateral wheelset displacements (with regard to stability), and increases as the displacements increase (to facilitate negotiation of curved sections). Using this kind of ideal conicity curve also enables a wheel profile to be secured where the contact points have a larger distribution over the active contact areas, making wear more homogeneous and reducing stresses. The result is a wheel profile with a conicity that is closer to the target conicity for both gauges studied, producing better curve negotiation while maintaining good stability on straight sections of track. The article shows the resultant wheel profile, the contact curves it produces, and a number of dynamic analyses demonstrating better dynamic behaviour of the synthesised wheel on curved sections with respect to the original wheel.     

Mathematical modelling of train–turnout interaction

The paper proposes a mathematical model of train–turnout interaction in the mid-frequency range (0–500 Hz). The model accounts for the effects of rail profile variation along the track and of local variation of track flexibility. The proposed approach is able to represent the condition of one wheel being simultaneously in contact with more than one rail, allowing the accurate prediction of the effect of wheels being transferred from one rail to another when passing over the switch toe and the crossing nose. Comprehensive results of train–turnout interaction during the negotiation of the main and the branch lines are presented, including the effect of wear of wheel/rail profiles and presence of track misalignment. In the final part of the paper, comparisons are performed between the results of numerical simulations and line measurements performed on two different turnouts for urban railway lines, showing a good agreement between experimental and numerical results.     

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.     

The effect of worn profile on wear progress of rail vehicle steel wheels over curved tracks

Dynamic performance, safety and maintenance costs of railway vehicles strongly depend on wheelset dynamics and particularly on the design of wheelset profile. This paper considers the effect of worn wheel profile on vehicle dynamics and the trend of wear in the wheels as a result of the vehicle movements. ADAMS/RAIL is used to build a multi-body system model of the vehicle. The track model is also configured as an elastic body. Measured new and worn wheel profiles are used to provide boundary conditions for the wheel/rail contacts. The fleet velocity profile taken during its normal braking is also used for the simulation. Wear numbers are calculated for different sets of wheels and the results compared with each other. Outcome of this research can be used for modifying dynamic performance of the vehicle, improving its suspension elements and increasing ride quality. It can also be further processed to reach to a modified wheel profile suitable for the fleet/track combination and for improved maintenance of the wheels. A major advantage of the computer models in this paper is the insertion of the wheel surface properties into the boundary conditions for dynamic modelling of the fleet. This is performed by regularly measuring the worn wheel profiles during their service life and by the calculation of the wear rate for individual wheels.     

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.     

Energy Dissipation Due to Vehicle/Track Interaction

SUMMARY

The effects of track irregularities and wheel profile on the amount of energy dissipated in railroad freight vehicles is examined. A nonlinear computational model is used to determine the average dissipation in the vehicle suspension and the wheel/rail contact patches. This dissipation is a component of the total resistance force acting on the vehicle. Parametric results are presented showing the effects of track geometry, wheel profile, suspension design, and hunting on train resistance. Track geometry studies consider the effects of track quality and curving. The AAR 1:20 wheel profile and the Heumann wheel profile are compared under various operating conditions. Compared with the Heumann profile, the AAR 1:20 profile is shown to have lower average resistance on good quality tangent track, but higher average resistance in steady curves. A trade-off exists between the two profiles when dynamic curve entry is considered.     

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.     

Energy Dissipation Due to Vehicle/Track Interaction

The effects of track irregularities and wheel profile on the amount of energy dissipated in railroad freight vehicles is examined. A nonlinear computational model is used to determine the average dissipation in the vehicle suspension and the wheel/rail contact patches. This dissipation is a component of the total resistance force acting on the vehicle. Parametric results are presented showing the effects of track geometry, wheel profile, suspension design, and hunting on train resistance. Track geometry studies consider the effects of track quality and curving. The AAR 1:20 wheel profile and the Heumann wheel profile are compared under various operating conditions. Compared with the Heumann profile, the AAR 1:20 profile is shown to have lower average resistance on good quality tangent track, but higher average resistance in steady curves. A trade-off exists between the two profiles when dynamic curve entry is considered.     

Dynamics of railway wagons subjected to braking torques on defective tracks

It is well known that track defects cause profound effects to the dynamics of railway wagons; normally such problems are examined for cases of wagons running at a constant speed. Brake/traction torques affect the speed profile due to the wheel–rail contact characteristics but most of the wagon–track interaction models do not explicitly consider them in simulation. The authors have recently published a model for the dynamics of wagons subject to braking/traction torques on a perfect track by explicitly considering the pitch degree of freedom for wheelsets. The model is extended for cases of lateral and vertical track geometry defects and worn railhead and wheel profiles. This paper presents the results of the analyses carried out using the model extended to the dynamics of wagons containing less ideal wheel profiles running on tracks with geometry defects and worn rails.     

Simulation of dynamic interaction between train and railway turnout

Dynamic train–track interaction is more complex in railway turnouts (switches and crossings) than that in ordinary tangent or curved tracks. Multiple contacts between wheel and rail are common, and severe impact loads with broad frequency contents are induced, when nominal wheel–rail contact conditions are disturbed because of the continuous variation in rail profiles and the discontinuities in the crossing panel. The absence of transition curves at the entry and exit of the turnout, and the cant deficiency, leads to large wheel–rail contact forces and passenger discomfort when the train is switching into the turnout track. Two alternative multibody system (MBS) models of dynamic interaction between train and a standard turnout design are developed. The first model is derived using a commercial MBS software. The second model is based on a multibody dynamics formulation, which may account for the structural flexibility of train and track components (based on finite element models and coordinate reduction methods). The variation in rail profile is accounted for by sampling the cross-section of each rail at several positions along the turnout. Contact between the back of the wheel flange and the check rail, when the wheelset is steered through the crossing, is considered. Good agreement in results from the two models is observed when the track model is taken as rigid.     

About Multi-Hertzian-Contact Hypothesis and Equivalent Conicity in the Case of S1002 and UIC60 Analytical Wheel/Rail Profiles

In a previous paper [1], a new method using 'multi-Hertzian' elasticity for identifying equivalent parameters of railway wheelsets rolling on defined tracks was presented. Another method also using Hertzian elasticity [2] has produced different results. Without knowing the origin of these discrepancies, the question was asked whether any of these methods could adequately describe actual non-Hertzian contacts. This paper, using as a basis Kalker's 'Contact' software which solves the non-Hertzian case for normal and tangential forces, answers this question, demonstrating that the assumption of summing Hertzian contacts is realistic provided that enough ellipses are considered. On this basis, the 'equivalent conicity' of the standard wheelset equipped with S1002 profiles rolling on UIC60 (1/40) rails is analysed in the range of small lateral displacements: around 2 mm it amounts to between 0.13 and 0.16. It is demonstrated that, due to actual wheel and rail elasticity, linearization of rigid profiles conicity, whatever method being used, is not realistic in practice.     

Parametric study of wheel transitions at railway crossings

Vehicle–track interaction at railway crossings is complex due to the discontinuity of the crossings. In this study, the effect of the local crossing geometry, the track alignment, and the wheel profiles on the wheel transition behaviour is investigated using the multi-body system software package VI-Rail. The transition behaviour is evaluated based on the location of the transition point along the crossing (and the location of impact), the contact pressure and the energy dissipation during the wheel–rail contact. A detailed parametric study of the crossing geometry has been performed, through which the most effective parameters for defining the crossing geometry are identified. These parameters are the cross-sectional shape of the nose rail, which can be tuned by one variable, and the vertical distance between the top of the wing rail and the nose rail. Additionally, a parametric study on the interaction influence of the crossing geometry, the track alignment and the wheel profile is performed using the design of experiments method with a two-level full factorial design. The longitudinal height profile of the crossing and the wheel profile are the most significant factors.     

Tangential force variation due to the bogie direction reversal procedure

A range of tangential forces is generated within the contact patch when a wheelset moves on the rail. These forces are intensified when incorporating curved tracks and motored axle rail vehicles [Arrus, P., de Pater, A.D. and Meyers, P., 2002, The stationary motion of a one-axle vehicle along a circular curve with real rail and wheel profiles. Vehicle System Dynamics, 37(1), 29–58]. The wheelset is subject to flange contact if an unbalanced force remains in a curve towards the high rail gauge face. The resultant force in the transverse direction includes the lateral force, the radial force, and the creep forces in addition to the effect of the frequent wheelset displacement due to the kinematic oscillation [Iwnicki, S., 2003, Simulation of wheel–rail contact forces. Fatigue Fracture Engineering Material Structure, 26, 887–900]. This article has focused on a potential variation in some of the forces cited when the wheelset is subject to backward and forward movements. A severe wear rate observed within the wheel flange region in Iranian Railways was investigated by operating a test bogie on a curvaceous track. An obvious improvement in the wear rate and wear pattern of the wheels was attained when the second test bogie encountered a bogie direction reversal procedure. This enhancement is considered in this article from the force analysis standpoint.     

Long-term wear and rolling contact fatigue behaviour of a conformal wheel profile designed for large radius curves

There are many reasons to optimise the wheel–rail interface through redesign or maintenance. Minimising wear and rolling contact fatigue (RCF) initiation on wheels and/or rails is often at the forefront of such considerations. This paper covers the design of a conformal wheel profile and its long-term wear and RCF performance to optimise the wheel–rail interface and subsequently reduce the occurrence of surface-initiated RCF on South Africa’s iron ore export line. A comparative study is performed using multibody dynamics simulation together with numerical wheel wear and RCF predictions. The advantages of a conformal wheel profile design are illustrated by evaluating the worn shape and resulting contact conditions of the conformal design. The conformal design has a steadier equivalent conicity progression and a smaller conicity range compared with the current wheel profile design over the wheel’s wear life. The combination of a conformal wheel profile design with 2?mm hollow wear and inadequate adherence to grinding tolerances often result in two-point contact, thereby increasing the probability of RCF initiation. The conformal wheel profile design proved to have wear and potential RCF benefits compared with the current wheel profile design. However, implementation of such a conformal wheel profile must be accompanied by improved rail grinding practices to ensure rail profile compliance.     

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.     

Modelling of wear and fatigue defect formation in wheel–rail contact

The model for analysing wear and fatigue defect formation is developed based on the approaches of contact and fracture mechanics. The model includes the solution of the contact problem for the wheel and rail to find the shape, size and position of the contact zones and the contact stresses and calculation of the surface wear and the function of damage accumulation in the rail and wheel. The wear rate and the worn-profile evolution of the wheel surface are calculated using both statistic and deterministic approaches to modelling of vehicle dynamics (tribo-dynamic modelling). The influence of the evolution of the wheel–rail profiles due to wear on the damage accumulation process is analysed. It is shown that for some values of the wear rate coefficient, the wear process can prevent the crack initiation under the wheel surface.