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.     

Deterioration mechanisms in the wheel–rail interface with focus on wear prediction: a literature review

Wheel–rail interface management is imperative to railway operation and its maintenance represents a major share of the total maintenance cost. In general, the course of events usually called wear is a complicated process involving several modes of material deterioration and contact surface alteration. Thus material removal or relocation, plastic flow and phase transformation may take place at, just below, or in-between the contacting surfaces. A higher degree of predictability of deterioration mechanisms and a firm basis for optimisation of the wheel–rail system are anticipated to reveal a great potential for cost savings. Wear in the sense of material loss and related wheel–rail profile evolution represents one of several modes of damage. The purpose of this survey is to explore research on wear simulation, to some degree extended to neighbouring disciplines. It is believed that a cross-disciplinary approach involving, for instance, adhesive and abrasive wear, surface plasticity, and rolling contact fatigue opens new perspectives to improved damage prediction procedures.     

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.     

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.     

An alternative to FASTSIM for tangential solution of the wheel–rail contact

In most rail vehicle dynamics simulation packages, tangential solution of the wheel–rail contact is gained by means of Kalker's FASTSIM algorithm. While 5–25% error is expected for creep force estimation, the errors of shear stress distribution, needed for wheel–rail damage analysis, may rise above 30% due to the parabolic traction bound. Therefore, a novel algorithm named FaStrip is proposed as an alternative to FASTSIM. It is based on the strip theory which extends the two-dimensional rolling contact solution to three-dimensional contacts. To form FaStrip, the original strip theory is amended to obtain accurate estimations for any contact ellipse size and it is combined by a numerical algorithm to handle spin. The comparison between the two algorithms shows that using FaStrip improves the accuracy of the estimated shear stress distribution and the creep force estimation in all studied cases. In combined lateral creepage and spin cases, for instance, the error in force estimation reduces from 18% to less than 2%. The estimation of the slip velocities in the slip zone, needed for wear analysis, is also studied. Since FaStrip is as fast as FASTSIM, it can be an alternative for tangential solution of the wheel–rail contact in simulation packages.     

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.     

Further study on the wheel–rail impact response induced by a single wheel flat: the coupling effect of strain rate and thermal stress

A comprehensive dynamic finite-element simulation method was proposed to study the wheel–rail impact response induced by a single wheel flat based on a 3-D rolling contact model, where the influences of the structural inertia, strain rate effect of wheel–rail materials and thermal stress due to the wheel–rail sliding friction were considered. Four different initial conditions (i.e. pure mechanical loading plus rate-independent, pure mechanical loading plus rate-dependent, thermo-mechanical loading plus rate-independent, and thermo-mechanical loading plus rate-dependent) were involved into explore the corresponding impact responses in term of the vertical impact force, von-Mises equivalent stress, equivalent plastic strain and shear stress. Influences of train speed, flat length and axle load on the flat-induced wheel–rail impact response were discussed, respectively. The results indicate that the maximum thermal stresses are occurred on the tread of the wheel and on the top surface of the middle rail; the strain rate hardening effect contributes to elevate the von-Mises equivalent stress and restrain the plastic deformation; and the initial thermal stress due to the sliding friction will aggravate the plastic deformation of wheel and rail. Besides, the wheel–rail impact responses (i.e. impact force, von-Mises equivalent stress, equivalent plastic strain, and XY shear stress) induced by a flat are sensitive to the train speed, flat length and axle load.     

Modelling of wheels and rail discontinuities in dynamic wheel–rail contact analysis

A classification of wheel–rail contact is given. Difference is made between modelling of a running wheel with continuous single-point-contact, as is common practice in wheel–rail contact analysis, and a wheel with transient double- or multi-point-contact, which may occur for rail irregularities with curvatures larger than that of the wheel circumference. It is shown that application of the first model for these irregularities will strongly underestimate the contact forces as it does not describe occurring mechanisms correctly. Further, it is shown that in principle it is not possible to describe the second type of contact fully correct with a lumped wheel model. Both wheel models are formulated mathematically for some basic contact cases. Afterwards, results are applied to a linear track model. Analytical closed-form solutions are found in the frequency domain for arbitrary type of contact and numerically transformed to the time domain. Finally, the necessity is shown to avoid situations where transient multiple-point-contact may occur (like rail joints) in practice.     

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.     

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.     

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.     

The role of the contact geometry in wheel–rail impact due to wheel flats: Part II

A classification of wheel flats according to the different stages of their growth is given, along with the characteristic features of the dynamic wheel–rail interaction for each category. Mathematical expressions and frequency spectra of the corresponding wheel mass trajectories are derived. Difference is made between the subcritical and the transcritical speed regime. A criterion is derived for contact loss for worn flats. Simulations show that the dynamic wheel–rail interaction is governed by the track stiffness for low train speeds or long flat lengths; for high speeds and/or short flat lengths the interaction is governed by the inertial properties of the wheel and the rail. For a given flat geometry, nonlinearities in the relationship between the impact magnitude and the train speed occur in the stiffness-dominated speed domain, whereas this relationship is approximately linear in the inertia-governed domain. In the latter domain, the impact magnitude is found to be linearly dependent upon the maximum trajectorial curvature or inversely linearly dependent on the minimum circumferential wheel tread curvature. The above relationships are valid for the subcritical speed regime, in which no contact loss occurs. Different contributions from the literature are compared with respect to the established relationship between impact magnitude and speed. Significant differences are found, due to insufficiently defined parameters and conditions. Conditions are derived for a consistent application of the so-called equivalent rail indentation in experiments with wheel flats, and the indirect strain registration method for measuring dynamic wheel–rail contact forces is reviewed.     

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.     

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.     

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.     

A fast-simplified wheel–rail contact model consistent with perfect plastic materials

A method is described which is an extension of rolling contact models with respect to plasticity. This new method, which is an extension of the STRIPES semi-Hertzian (SH) model, has been implemented in a multi-body-system (MBS) package and does not result in a longer execution time than the STRIPES SH model [J.B. Ayasse and H. Chollet, Determination of the wheel–rail contact patch in semi-Hertzian conditions, Veh. Syst. Dyn. 43(3) (2005), pp. 161–172]. High speed of computation is obtained by some hypotheses about the plastic law, the shape of stresses, the locus of the maximum stress and the slip. Plasticity does not change the vehicle behaviour but there is a need for an extension of rolling contact models with respect to plasticity as far as fatigue analysis of rail is concerned: rolling contact fatigue may be addressed via the finite element method (FEM) including material non-linearities, where loads are the contact stresses provided by the post-processing of MBS results [K. Dang Van, M.H. Maitournam, Z. Moumni, and F. Roger, A comprehensive approach for modeling fatigue and fracture of rails, Eng. Fract. Mech. 76 (2009), pp. 2626–2636]. In STRIPES, like in other MBS models, contact stresses may exceed the plastic yield criterion, leading to wrong results in the subsequent FEM analysis. With the proposed method, contact stresses are kept consistent with a perfect plastic law, avoiding these problems. The method is benchmarked versus non-linear FEM in Hertzian geometries. As a consequence of taking plasticity into account, contact patch area is bigger than the elastic one. In accordance with FEM results, a different ellipse aspect ratio than the one predicted by Hertz theory was also found and finally pressure does not exceed the threshold prescribed by the plastic law. The method also provides more exact results with non-Hertzian geometries. The new approach is finally compared with non-linear FEM in a tangent case with a unidirectional load and a complete slip: when plasticity is taken into account, and for large adhesion values, friction forces have an influence on the size of the contact patch. The proposed approach enables also to assess extensively the level of plasticity along a track through an indicator associated with a given yield stress.     

A linear complementarity method for the solution of vertical vehicle–track interaction

A new method is proposed for the solution of the vertical vehicle–track interaction including a separation between wheel and rail. The vehicle is modelled as a multi-body system using rigid bodies, and the track is treated as a three-layer beam model in which the rail is considered as an Euler-Bernoulli beam and both the sleepers and the ballast are represented by lumped masses. A linear complementarity formulation is directly established using a combination of the wheel–rail normal contact condition and the generalised-α method. This linear complementarity problem is solved using the Lemke algorithm, and the wheel–rail contact force can be obtained. Then the dynamic responses of the vehicle and the track are solved without iteration based on the generalised-α method. The same equations of motion for the vehicle and track are adopted at the different wheel–rail contact situations. This method can remove some restrictions, that is, time-dependent mass, damping and stiffness matrices of the coupled system, multiple equations of motion for the different contact situations and the effect of the contact stiffness. Numerical results demonstrate that the proposed method is effective for simulating the vehicle–track interaction including a separation between wheel and rail.     

A robust non-Hertzian contact method for wheel–rail normal contact analysis

A modified Kik–Piotrowski (MKP) model is proposed in this paper for an accurate and robust calculation of wheel–rail normal contact problem. The presented method is able to consider the relationship between the elastic deformation of a line and the normal pressure distribution within the contact patch. A novel shape correction method is put forward to correctly describe the elastic deformation of the contact patch. Taking the results estimated by Kalker’s variational method and Kik–Piotrowski method as references, the proposed method is validated by three contact cases, including the assumed standardised non-Hertzian contact and the two-point contact, as well as the contact behaviours based on three actual wheel–rail profiles. The simulation results indicate that, compared with Kik–Piotrowski method, the proposed MKP method achieves better agreement with Kalker’s variational method. Moreover, the MKP method can avoid the abrupt change of wheel–rail normal force due to the sudden transfer of the contact point, which contributes to a better computational stability.