Simulations of roughness growth on rails - results from a 2D non-Hertzian, non-steady contact model

A new model for simulating rail roughness growth on tangent track is presented in this paper. The model consists of three relatively independent components: (1) a time-domain vehicle/track interaction model; (2) a 2D non-Hertzian and non-steady wheel/rail contact model; and (3) a wear model. Wheel/rail contact forces for a given initial roughness obtained from the vehicle/track interaction model are used by the contact model to calculate the contact patch size, normal pressure and tangential stresses with material removal assumed to be linearly proportional to the friction work in the contact patch. The roughness profile is updated and fed back into vehicle/track interaction model. The 2D contact model is initially compared with a 3D model for various wavelength of initial sinusoidal roughness. Long term roughness growth is then simulated with the 2D contact model. Simulation shows that all initial sinusoidal roughness of wavelengths between 20-100 mm are levelled out. The wavelength-fixing mechanism, that has previously been used to explain the cause of corrugation, is not found in the present investigations.     

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.     

A locomotive–track coupled vertical dynamics model with gear transmissions

A gear transmission system is a key element in a locomotive for the transmission of traction or braking forces between the motor and the wheel–rail interface. Its dynamic performance has a direct effect on the operational reliability of the locomotive and its components. This paper proposes a comprehensive locomotive–track coupled vertical dynamics model, in which the locomotive is driven by axle-hung motors. In this coupled dynamics model, the dynamic interactions between the gear transmission system and the other components, e.g. motor and wheelset, are considered based on the detailed analysis of its structural properties and working mechanism. Thus, the mechanical transmission system for power delivery from the motor to the wheelset via gear transmission is coupled with a traditional locomotive–track dynamics system via the wheel–rail contact interface and the gear mesh interface. This developed dynamics model enables investigations of the dynamic performance of the entire dynamics system under the excitations from the wheel–rail contact interface and/or the gear mesh interface. Dynamic interactions are demonstrated by numerical simulations using this dynamics model. The results indicate that both of the excitations from the wheel–rail contact interface and the gear mesh interface have a significant effect on the dynamic responses of the components in this coupled dynamics system.     

A new method for wheel–rail contact force continuous measurement using instrumented wheelset

With the aim of improving the continuous measurement of wheel–rail contact force by instrumented wheelset, instead of solving the non-linear equations, we proposed a new method based on state space theory. With this new method, the wheel–rail contact force can be calculated by the recurrence relation and the signals from strain gauge bridges on wheel web. The implementation of continuous instrumented wheelset is quite general and simplified, due to the specific bridging scheme is not necessary. It means that continuous measurement of the contact force could be realised with a simple bridging scheme, even as simple as discrete instrumented scheme. In this work, we first demonstrated and discussed the effectiveness and accuracy of this new method by estimation results with the numerical simulations, and we also applied this new method to two field tests, including one was conducted in a loop test line using a high-speed train and the other one was conducted in an urban line with a light rail vehicle. In a word, this new method is proved to be an effective way to monitor the wheel–rail contact force of rail vehicle track system.     

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.     

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.     

Sensitivity of the wheel–rail contact interactions and Dang Van Fatigue Index in the rail with respect to irregularities of the track geometry

This paper investigates the effects of the track geometry irregularities on the wheel–rail dynamic interactions and the rail fatigue initiation through the application of the Dang Van criterion, that supposes an elastic shakedown of the structure. The irregularities are modelled, using experimental data, as a stochastic field which is representative of the considered railway network. The tracks thus generated are introduced as the input of a railway dynamics software to characterise the stochastic contact patch and the parameters on which it depends: contact forces and wheelset–rail relative position. A variance-based global sensitivity analysis is performed on quantities of interest representative of the dynamic behaviour of the system, with respect to the stochastic geometry irregularities and for different curve radius classes and operating conditions. The estimation of the internal stresses and the fatigue index being more time-consuming than the dynamical simulations, the sensitivity analysis is performed through a metamodel, whose input parameters are the wheel–rail relative position and velocity. The coefficient of variation of the number of fatigue cycles, when the simulations are performed with random geometry irregularities, varies between 0.13 and 0.28. In a large radius curve, the most influent irregularity is the horizontal curvature, while, in a tight curve, the gauge becomes more important.     

High-speed wheel/rail contact determining method with rotating flexible wheelset and validation under wheel polygon excitation

A rotating flexible wheelset model is developed and integrated into a vehicle/track dynamic model. Flexible wheelset modes with natural frequencies less than 1000?Hz are considered in the wheelset modelling. An innovation of the paper is that wheel/rail rolling contact calculation considers the effect of the wheelset flexibility and the rotating effect. By introducing two half dummy rigid wheelsets the rolling contact between the flexible wheelset and the two rails can be transformed to that between a rigid wheelset and the rails. As an extension application, the wheel OOR (Out-Of-Round) wears with the 11th, 15th, and 17th orders are used to the vehicle system dynamic model with rigid, flexible and rotating-flexible wheelset model. The results of the three models are compared to study the influence of wheelset flexibility and rotation. The ‘online searching contact method’ developed in this paper is compared with the traditional contact method with considering the rotating flexible wheelset. And then a measured OOR is used to excite the rotating flexible wheelset, the response of which is analysed and verified.     

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.     

The vertical and the longitudinal dynamic responses of the vehicle–track system to squat-type short wavelength irregularity

The squat, a kind of rolling contact fatigue occurring on the rail top, can excite the high-frequency vehicle–track interaction effectively due to its geometric deviations with a typical wavelength of 20–40 mm, leading to the accelerated deterioration of a track. In this work, a validated 3D transient finite element model is employed to calculate in the time domain the vertical and the longitudinal dynamic contact forces between the wheel and the rail caused by squats. The vehicle–track structure and the wheel–rail continua are both considered in order to include all the important eigencharacteristics of the system related to squats. By introducing the rotational and translational movements of the wheel, the transient wheel–rail rolling contact is solved in detail by a 3D frictional contact model integrated. The contact filter effect is considered automatically in the simulations by the finite size of the contact patch. The present work focuses on the influences of the length, width and depth of a light squat on the resulted dynamic contact forces, for which idealised defect models are used. The growth of a squat is also modelled to a certain extent by a series of defects with different dimensions. The results show that the system is mainly excited at two frequencies separately in the vertical and the longitudinal dynamics. Their superposition explains the typical appearance of mature squats. As a squat grows up, the magnitude of the excited vibration at the lower frequency increases faster than the one at the higher frequency.     

An inverse railway wagon model and its applications

An inverse wagon model was developed to estimate wheel–rail contact forces using only measurements of wagon body responses as inputs. The purpose of this work was to provide mathematical modelling to embed in low-cost devices that can be mounted on each freight wagon in a large wagon fleet. To minimize cost, complication, and the maintenance inconvenience of these devices, the constraint is imposed that transducers and connections are limited to locations on the wagon body. Inputs to the inverse model developed include only vertical and lateral translational accelerations and angular accelerations of roll, pitch, and yaw of the wagon body. The model combines the integration and partial modal matrix (PMM) techniques together to form an IPMM method. Besides wheel–rail contact forces some motion quantities such as the lateral and yaw displacements of wheelset are also predicted. Results from the inverse model were compared with data from full scale laboratory suspension tests for vertical suspension excitations. The inverse model was also compared with results from simulations completed in VAMPIRE^® for more complicated track input profiles. The model results and the applications of the model are discussed.     

Short wavelength rail corrugation and non-steady-state contact mechanics

J. J. Kalker has been the first to consider non-steady-state or transient contact mechanics. Based on Kalker the second author developed a linear contact model for the non-steady-state rolling contact of a wheel running over slightly corrugated rails. The theoretical investigations are concentrated on linear, non-steady-state contact mechanics superimposed to a nonlinear reference state. The reference state is given by the running behaviour of a wheelset due to traction, curving or hunting. For the linear, non-steady-state analysis Kalker's theory has to be modified to predict wear rates in dependency of the corrugation wavelengths. As a result corrugations are only amplified in the range between 2 and 10 cm. Therefore, non-steady-state contact mechanics and wear are responsible for a wavelength fixing mechanism. Structural mechanics of the rail indicate that wavelength in this range is predominantly amplified.     

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.     

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.     

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.     

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.     

Dynamic Train-Track Interaction: Variability Attributable to Scatter in the Track Properties

A numerical method to simulate vertical dynamic interaction between a rolling train and a railway track has been used to investigate the influence of stochastic properties of the track structure. A perturbation technique has been used to investigate the influence of the scatter in selected track properties. The train-track interaction problem has been numerically solved by use of an extended state-space vector approach in conjunction with a complex modal superposition for the whole track structure. All numerical simulations have been carried out in the time-domain with a moving mass model. Properties such as rail pad stiffness, ballast stiffness, dynamic ballast-subgrade mass and sleeper spacing have been studied. To obtain sufficient statistical information from track structures, full-scale measurements in the field and laboratory measurements have been carried out. The influence of scatter in the track properties on the maximum contact force between the rail and the wheel, the maximum magnitude of the vertical wheelset acceleration, and the maximum sleeper displacement have been studied. Mean values and standard deviations of these quantities have been calculated. The effects of the variation of the investigated track properties are discussed.     

A non-Hertzian method for solving wheel–rail normal contact problem taking into account the effect of yaw

A novel approach is proposed in this paper to deal with non-Hertzian normal contact in wheel–rail interface, extending the widely used Kik–Piotrowski method. The new approach is able to consider the effect of the yaw angle of the wheelset against the rail on the shape of the contact patch and on pressure distribution. Furthermore, the method considers the variation of profile curvature across the contact patch, enhancing the correspondence to CONTACT for highly non-Hertzian contact conditions. The simulation results show that the proposed method can provide more accurate estimation than the original algorithm compared to Kalker’s CONTACT, and that the influence of yaw on the contact results is significant under certain circumstances.     

An efficient approach to the analysis of rail surface irregularities accounting for dynamic train–track interaction and inelastic deformations

A two-dimensional computational model for assessment of rolling contact fatigue induced by discrete rail surface irregularities, especially in the context of so-called squats, is presented. Dynamic excitation in a wide frequency range is considered in computationally efficient time-domain simulations of high-frequency dynamic vehicle–track interaction accounting for transient non-Hertzian wheel–rail contact. Results from dynamic simulations are mapped onto a finite element model to resolve the cyclic, elastoplastic stress response in the rail. Ratcheting under multiple wheel passages is quantified. In addition, low cycle fatigue impact is quantified using the Jiang–Sehitoglu fatigue parameter. The functionality of the model is demonstrated by numerical examples.