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.     

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.     

A study on high-speed rolling contact between a wheel and a contaminated rail

A 3-D explicit finite element model is developed to investigate the transient wheel–rail rolling contact in the presence of rail contamination or short low adhesion zones (LAZs). A transient analysis is required because the wheel passes by a short LAZ very quickly, especially at high speeds. A surface-to-surface contact algorithm (by the penalty method) is employed to solve the frictional rolling contact between the wheel and the rail meshed by solid elements. The LAZ is simulated by a varying coefficient of friction along the rail. Different traction efforts and action of the traction control system triggered by the LAZ are simulated by applying a time-dependent driving torque to the wheel axle. Structural flexibilities of the vehicle–track system are considered properly. Analysis focuses on the contact forces, creepage, contact stresses and the derived frictional work and plastic deformation. It is found that the longitudinal contact force and the maximum surface shear stress in the contact patch become obviously lower in the LAZ and much higher as the wheel re-enters the dry rail section. Consequently, a higher wear rate and larger plastic flow are expected at the location where the dry contact starts to be rebuilt. In other words, contact surface damages such as wheel flats and rail burns may come into being because of the LAZ. Length of the LAZ, the traction level, etc. are varied. The results also show that local contact surface damages may still occur as the traction control system acts.     

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.     

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.     

Dynamic simulation of railway vehicles: wheel/rail contact analysis

The multibody simulation of railway vehicle dynamics needs a reliable and efficient method to evaluate the contact points between wheel and rail, because their positions have a considerable influence on the direction and intensity of the contact forces. In this work, an innovative semi-analytic procedure for the detection of the wheel/rail contact points (named the DIFF method) is presented. This method considers the wheel and the rail as two surfaces whose analytic expressions are known and is based on the idea that in the contact points the difference between the surfaces has local minima and is equivalent to solving an algebraic two-dimensional system. The original problem can be reduced analytically to a simple scalar equation that can be easily solved numerically (since the problem dimension is one, even elementary non-iterative algorithms can be efficient).     

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.     

A third-order approximation method for three-dimensional wheel–rail contact

Multibody train analysis is used increasingly by railway operators whenever a reliable and time-efficient method to evaluate the contact between wheel and rail is needed; particularly, the wheel–rail contact is one of the most important aspects that affects a reliable and time-efficient vehicle dynamics computation. The focus of the approach proposed here is to carry out such tasks by means of online wheel–rail elastic contact detection. In order to improve efficiency and save time, a main analytical approach is used for the definition of wheel and rail surfaces as well as for contact detection, then a final numerical evaluation is used to locate contact. The final numerical procedure consists in finding the zeros of a nonlinear function in a single variable. The overall method is based on the approximation of the wheel surface, which does not influence the contact location significantly, as shown in the paper.     

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.     

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.     

Transient Dynamic Characteristics of a Non-Pneumatic Mechanical Elastic Wheel Rolling Over a Ditch

The transient dynamic characteristic of a tire, which has a significant effect on vehicle handling stability and ride comfort, is difficult to study in detail because of its highly non-linear behavior. In this study, the transient dynamic characteristics of a non-pneumatic wheel, called the mechanical elastic wheel (MEW), which was rolling over a ditch were investigated by the explicit dynamic finite element (FE) method. A three-dimensional FE model of MEW considering geometric nonlinearity, material nonlinearity and large contact deformation between the wheel and the road, was established. For the validation of the accuracy and reliability of the FE model of MEW, the simulation and the experimental results of the radial stiffness and footprint of MEW were compared and analyzed. A dynamic simulation of the validated FE model of MEW rolling over a ditch was conducted using the ABAQUS/Explicit program. The equivalent stress and the contact stress generated during the process of the rolling MEW impacting the ditch were studied in detail. The effect of the rolling speed on the transient dynamic characteristics was also analyzed based on the simulation results. The simulation results could provide guidance for the optimization of the MEW structure and vehicle 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.     

On the interest of integrating vehicle dynamics for the ground propagation of vibrations: the case of urban railway traffic

This paper presents a complete numerical model for studying the vertical dynamics of the vehicle/track interaction and its impact on the surrounding soil, with the emphasis on vehicle modelling. A decoupling between the track and the soil is proposed, due to the difficulty of considering all the subsystem components. The train/track model is based on a multibody model (for the vehicle) and a finite element model (for the track). The soil is modelled using an infinite/finite element approach. Simulations of both models are carried out in the time domain, which is better able to simulate the propagation of the vibration waves and to take into account the possible nonlinearity of a component. The methodology is applied in the case of an urban tram track and validated with the available experimental data. Models for the tram, the track and the soil are described. Results from the complete model of the vehicle and a simple model, based on an axle load, are compared with experimental results and the benefits of a complete model in the simulation of the ground vibration propagation induced by railway vehicles are demonstrated. Moreover, a parametric study of the vehicle wheel type is conducted, which shows the advantage of a resilient wheel, for various rail defects.     

Dynamic vehicle–track interaction in switches and crossings and the influence of rail pad stiffness – field measurements and validation of a simulation model

A model for simulation of dynamic interaction between a railway vehicle and a turnout (switch and crossing, S&C) is validated versus field measurements. In particular, the implementation and accuracy of viscously damped track models with different complexities are assessed. The validation data come from full-scale field measurements of dynamic track stiffness and wheel–rail contact forces in a demonstrator turnout that was installed as part of the INNOTRACK project with funding from the European Union Sixth Framework Programme. Vertical track stiffness at nominal wheel loads, in the frequency range up to 20?Hz, was measured using a rolling stiffness measurement vehicle (RSMV). Vertical and lateral wheel–rail contact forces were measured by an instrumented wheel set mounted in a freight car featuring Y25 bogies. The measurements were performed for traffic in both the through and diverging routes, and in the facing and trailing moves. The full set of test runs was repeated with different types of rail pad to investigate the influence of rail pad stiffness on track stiffness and contact forces. It is concluded that impact loads on the crossing can be reduced by using more resilient rail pads. To allow for vehicle dynamics simulations at low computational cost, the track models are discretised space-variant mass–spring–damper models that are moving with each wheel set of the vehicle model. Acceptable agreement between simulated and measured vertical contact forces at the crossing can be obtained when the standard GENSYS track model is extended with one ballast/subgrade mass under each rail. This model can be tuned to capture the large phase delay in dynamic track stiffness at low frequencies, as measured by the RSMV, while remaining sufficiently resilient at higher frequencies.     

Extended applications of track irregularity probabilistic model and vehicle–slab track coupled model on dynamics of railway systems

Track irregularities are inevitably in a process of stochastic evolution due to the uncertainty and continuity of wheel–rail interactions. For depicting the dynamic behaviours of vehicle–track coupling system caused by track random irregularities thoroughly, it is a necessity to develop a track irregularity probabilistic model to simulate rail surface irregularities with ergodic properties on amplitudes, wavelengths and probabilities, and to build a three-dimensional vehicle–track coupled model by properly considering the wheel–rail nonlinear contact mechanisms. In the present study, the vehicle–track coupled model is programmed by combining finite element method with wheel–rail coupling model firstly. Then, in light of the capability of power spectral density (PSD) in characterising amplitudes and wavelengths of stationary random signals, a track irregularity probabilistic model is presented to reveal and simulate the whole characteristics of track irregularity PSD. Finally, extended applications from three aspects, that is, extreme analysis, reliability analysis and response relationships between dynamic indices, are conducted to the evaluation and application of the proposed models.     

The Elastic Cross-Influence between Two Quasi-Hertzian Contact Zones

In this paper a type of contact between two bodies is considered, which leads to the formation of two separate contact zones. The contact zones considered are Hertzian ellipses if the distance between them is large. When the distance between the zones is finite it is necessary to take into account the elastic cross-influence of the two zones. In this paper an approximate method is proposed which allows the determination of the actual contact zones as Hertz's ellipses but without disregarding the cross-influence. Numerical results are presented for two, two-humped bodies pressed against each other and for steady state rolling of a flanged wheel along a steel rail.

It is shown that neglecting the cross influence for a flanged wheel leads to a negligible error in the contact forces but the forces are transmitted between wheel and rail through contact patches which are generally more slender than with the cross-influence neglected.     

A novel method to model wheel–rail normal contact in vehicle dynamics simulation

An approximate analytical method is proposed for calculating the contact patch and pressure distribution in the wheel–rail interface. The deformation of the surfaces in contact is approximated using the separation between them. This makes it possible to estimate the contact patch analytically. The contact pressure distribution in the rolling direction is assumed to be elliptic with its maximum calculated by applying Hertz' solution locally. The results are identical to Hertz's for elliptic cases. In non-elliptic cases good agreement is achieved in comparison to the more accurate but computationally expensive Kalker's variational method (CONTACT code). Compared to simplified non-elliptic contact methods based on virtual penetration, the calculated contact patch and pressure distribution are markedly improved. The computational cost of the proposed method is significantly lower than the more detailed methods, making it worthwhile to be applied to rolling contact in rail vehicle dynamics simulation. Such fast and accurate estimation of contact patch and pressure paves the way for on-line modelling of damage phenomena in dynamics simulation packages.     

A simplified model of wheel/rail contact mechanics for non-Hertzian problems and its application in rail vehicle dynamic simulations

The presented model assumes semi-elliptical normal pressure distribution in the direction of rolling. The contact area is found by virtual penetration of wheel and rail. The normal pressure is calculated by satisfying contact conditions at the geometrical point of contact. The calculation is non-iterative, fast and completely reliable. It may be carried out on-line in MultiBody Systems (MBS) computer codes. The tests using the programme CONTACT by Kalker and experience from application in MBS codes show that the model is suitable for technical applications. The creep forces have been calculated with the FASTSIM algorithm, adapted for a non-elliptical contact area. Some applications in rail vehicle dynamics and wear simulation have been outlined.     

Wheel/Rail Non-linear Interactions With Coupling Between Vertical and Lateral Directions

Summary A theoretical model is developed to explore the high frequency wheel/rail interaction with coupling between the vertical and lateral directions. This coupling is introduced through the track dynamics due to the offset of the wheel/rail contact point from the rail centre line. Equivalent models of the railway track in the time domain are developed according to the rail vibration receptances in the frequency domain. The wheel is represented by a mass in each direction with no vertical-lateral coupling. The vertical wheel/rail interaction is generated through a non-linear Hertzian contact stiffness, allowing for the possibility of loss of contact between the wheel and rail. The lateral interaction is represented by a contact spring and a creep force damper in series and their values depend on the vertical contact force. The vibration source is the roughness on the wheel and rail contact surfaces which forms a relative displacement excitation in the vertical direction. Using the combined interaction model with this relative displacement excitation, the wheel/rail interactions with coupling between the vertical and lateral vibrations are simulated. It is found that the lateral interaction force caused by the offset is usually less than thirty percent of the vertical dynamic force. The lateral vibration of the rail is significantly reduced due to the presence of the lateral coupling, whereas the vertical interaction is almost unaffected by the lateral force.     

Influence of wheel–rail contact modelling on vehicle dynamic simulation

This paper presents a comparison of four models of rolling contact used for online contact force evaluation in rail vehicle dynamics. Until now only a few wheel–rail contact models have been used for online simulation in multibody software (MBS). Many more models exist and their behaviour has been studied offline, but a comparative study of the mutual influence between the calculation of the creep forces and the simulated vehicle dynamics seems to be missing. Such a comparison would help researchers with the assessment of accuracy and calculation time. The contact methods investigated in this paper are FASTSIM, Linder, Kik–Piotrowski and Stripes. They are compared through a coupling between an MBS for the vehicle simulation and Matlab for the contact models. This way the influence of the creep force calculation on the vehicle simulation is investigated. More specifically this study focuses on the influence of the contact model on the simulation of the hunting motion and on the curving behaviour.