A survey of wheel–rail contact models for rail vehicles

Accurate and efficient contact models for wheel–rail interaction are essential for the study of the dynamic behaviour of a railway vehicle. Assessment of the contact forces and moments, as well as contact geometry provide a fundamental foundation for such tasks as design of braking and traction control systems, prediction of wheel and rail wear, and evaluation of ride safety and comfort. This paper discusses the evolution and the current state of the theories for solving the wheel–rail contact problem for rolling stock. The well-known theories for modelling both normal contact (Hertzian and non-Hertzian) and tangential contact (Kalker's linear theory, FASTSIM, CONTACT, Polach's theory, etc.) are reviewed. The paper discusses the simplifying assumptions for developing these models and compares their functionality. The experimental studies for evaluation of contact models are also reviewed. This paper concludes with discussing open areas in contact mechanics that require further research for developing better models to represent the wheel–rail interaction.     

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.     

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.     

Comments on ‘the Kalker book of tables for non-Hertzian contact of wheel and rail’

The ‘simple double-elliptical contact’ (SDEC) approach by Piotrowski et al. [The Kalker book of tables for non-Hertzian contact of wheel and rail. Vehicle Syst Dyn. 2017;55:875–901] generates a-symmetrical contact patches in an elegant way. This allows to extend the table-based approach for the wheel–rail creep force calculation towards non-Hertzian contact geometry. This is an important line of research, because FASTSIM is intricate for non-Hertzian contacts, whereas CONTACT requires long calculation times.

Here, we comment on the further motivation that's provided for the approach. According to the authors, ‘the spin creepage generates longitudinal creep force in non-symmetric, non-elliptical contacts’, which is ‘completely lost’ when using elliptical regularisation. We demonstrate that this mainly depends on the choice of contact origin, and that the interaction is much reduced if different choices are made. This suggests that elliptical regularisation may be viable still, if the details are properly worked out. Furthermore, we introduce the spin center and the free-rolling position as means to extend the table-based approach towards more general non-Hertzian circumstances.     

Assessment of a semi-Hertzian method for determination of wheel–rail contact patch

Wheel–rail contact calculations are essential for simulating railway vehicle dynamic behavior. Currently, these simulations usually use the Hertz contact theory to calculate normal forces and Kalker's ‘FASTSIM’ program to evaluate tangential stresses. Since 1996, new methods called semi-Hertzian have appeared: 5 Kik, W. and Piotrowski, J. A fast approximate method to calculate normal load at contact between wheel and rail and creep forces during rolling. Paper presented at the 2nd Mini-conference on Contact Mechanics and Wear of Rail/Wheel Systems. July29–31, Budapest.  [Google Scholar] 7 Ayasse, J. B., Chollet, H. and Maupu, J. L. 2000. Paramètres caractéristiques du contact roue-rail. Rapport de Recherche INRETS n225, ISSN 0768–9756 (in French) [Google Scholar] (STRIPES). These methods attempt to estimate the non-elliptical contact patches with a discrete extension of the Hertz theory. As a continuation of 2 Ayasse, J. B and Chollet, H. 2005. Determination of the wheel–rail contact patch in semi-Hertzian conditions. Vehicle System Dynamics, 43(3) [Google Scholar], a validation of the STRIPES method for normal problem computing on three test cases is proposed in this article. The test cases do not fulfill the hypothesis required for the Hertz theory. Then, the Kalker's FASTSIM algorithm is adapted to STRIPES patch calculus to perform tangential forces computation. This adaptation is assessed using Kalker's CONTACT algorithm.     

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.     

Systems-on-chip approach for real-time simulation of wheel–rail contact laws

This paper presents the development of a systems-on-chip approach to speed up the simulation of wheel–rail contact laws, which can be used to reduce the requirement for high-performance computers and enable simulation in real time for the use of hardware-in-loop for experimental studies of the latest vehicle dynamic and control technologies. The wheel–rail contact laws are implemented using a field programmable gate array (FPGA) device with a design that substantially outperforms modern general-purpose PC platforms or fixed architecture digital signal processor devices in terms of processing time, configuration flexibility and cost. In order to utilise the FPGA's parallel-processing capability, the operations in the contact laws algorithms are arranged in a parallel manner and multi-contact patches are tackled simultaneously in the design. The interface between the FPGA device and the host PC is achieved by using a high-throughput and low-latency Ethernet link. The development is based on FASTSIM algorithms, although the design can be adapted and expanded for even more computationally demanding tasks.     

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.     

Design and optimisation of wheel–rail profiles for adhesion improvement

This paper describes a study for the optimisation of the wheel profile in the wheel–rail system to increase the overall level of adhesion available at the contact interface, in particular to investigate how the wheel and rail profile combination may be designed to ensure the improved delivery of tractive/braking forces even in poor contact conditions. The research focuses on the geometric combination of both wheel and rail profiles to establish how the contact interface may be optimised to increase the adhesion level, but also to investigate how the change in the property of the contact mechanics at the wheel–rail interface may also lead to changes in the vehicle dynamic behaviour.     

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.     

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.     

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.     

Effect of temperature- and frequency-dependent dynamic properties of rail pads on high-speed vehicle–track coupled vibrations

The soft under baseplate pad of WJ-8 rail fastener frequently used in China’s high-speed railways was taken as the study subject, and a laboratory test was performed to measure its temperature and frequency-dependent dynamic performance at 0.3?Hz and at ?60°C to 20°C with intervals of 2.5°C. Its higher frequency-dependent results at different temperatures were then further predicted based on the time–temperature superposition (TTS) and Williams–Landel–Ferry (WLF) formula. The fractional derivative Kelvin–Voigt (FDKV) model was used to represent the temperature- and frequency-dependent dynamic properties of the tested rail pad. By means of the FDKV model for rail pads and vehicle–track coupled dynamic theory, high-speed vehicle–track coupled vibrations due to temperature- and frequency-dependent dynamic properties of rail pads was investigated. Finally, further combining with the measured frequency-dependent dynamic performance of vehicle’s rubber primary suspension, the high-speed vehicle–track coupled vibration responses were discussed. It is found that the storage stiffness and loss factor of the tested rail pad are sensitive to low temperatures or high frequencies. The proposed FDKV model for the frequency-dependent storage stiffness and loss factors of the tested rail pad can basically meet the fitting precision, especially at ordinary temperatures. The numerical simulation results indicate that the vertical vibration levels of high-speed vehicle–track coupled systems calculated with the FDKV model for rail pads in time domain are higher than those calculated with the ordinary Kelvin–Voigt (KV) model for rail pads. Additionally, the temperature- and frequency-dependent dynamic properties of the tested rail pads would alter the vertical vibration acceleration levels (VALs) of the car body and bogie in 1/3 octave frequencies above 31.5?Hz, especially enlarge the vertical VALs of the wheel set and rail in 1/3 octave frequencies of 31.5–100?Hz and above 315?Hz, which are the dominant frequencies of ground vibration acceleration and rolling noise (or bridge noise) caused by high-speed railways respectively. Since the fractional derivative value of the adopted rubber primary suspension, unlike the tested rail pad, is very close to 1, its frequency-dependent dynamic performance has little effect on high-speed vehicle–track coupled vibration responses.     

High-frequency random vibration analysis of a high-speed vehicle–track system with the frequency-dependent dynamic properties of rail pads using a hybrid SEM–SM method

A hybrid Spectral Element Method (SEM)–Symplectic Method(SM) method for high-efficiency computation of the high-frequency random vibrations of a high-speed vehicle–track system with the frequency-dependent dynamic properties of rail pads is presented. First, the Williams-Landel-Ferry (WLF) formula and Fractional Derivative Zener (FDZ) model were, respectively, applied for prediction and representation of the frequency-dependent dynamic properties of Vossloh 300 rail pads frequently used in China's high-speed railway. Then, the proposed hybrid SEM–SM method was used to investigate the influence of the frequency-dependent dynamic performance of Vossloh 300 rail pads on the high-frequency random vibrations of high-speed vehicle–track systems at various train speeds or different levels of rail surface roughness. The experimental results indicate that the storage stiffness and loss factors of Vossloh 300 rail pad increase with the decrease in dynamic loads or the increase in preloads within 0.1–10,000?Hz at 20°C, and basically linearly increase with frequency in a logarithmic coordinate system. The results computed by the hybrid SEM–SM method demonstrate that the frequency-dependent viscous damping of Vossloh 300 rail pads, compared with its constant viscous damping and frequency-dependent stiffness, has a much more conspicuous influence on the medium-frequency (i.e. 20–63?Hz) random vibrations of car bodies and rail fasteners, and on the mid- (i.e. 20–63?Hz) and high-frequency (i.e. 630–1250?Hz) random vibrations of bogies, wheels and rails, especially with the increase in train speeds or the deterioration of rail surface roughness. The two sensitive frequency bands can also be validated by frequency response function (FRF) analysis of the proposed infinite rail–fastener model. The mid and high frequencies influenced by the frequency-dependent viscous damping of rail pads are exactly the dominant frequencies of ground vibration acceleration and wheel rolling noise caused by high-speed railways, respectively. Even though the existing time-domain (or frequency-domain) finite track models associated with the time-domain (or frequency-domain) fractional derivative viscoelastic (FDV) models of rail pads can also be used to reach the same conclusions, the hybrid SEM–SM method in which only one element is required to compute the high-order vibration modes of infinite rail is more appropriate for high-efficiency analysis of the high-frequency random vibrations of high-speed vehicle–track systems.     

Images of the Year

前几期，本刊陆续转载了由ENR杂志评选出的2012年“施工照片竞赛”的前20名作品。近期整理了后面部分，竟发现不少丝毫不逊色于前者的惊艳之作，现摭拾若干，使其免遭遗珠之憾。     

An Analysis of Vehicular Engines From the Viewpoint of the Tenth-Five-Year Plan

It is stimulated in the Tenth-Five-Year-Plan period that by the year of 2010, the Chinese auto industry should become the pillar industry by means of integration of competition and self-development. According to incomplete statistics, by 2005, the demand will be: 3.6 million units of motor vehicle engines, of which petrol engines will be 2.288 million units, diesel engines 1.232 million units and the substitute engines 80,000 units. And the     

Guidelines for the Development of the Auto Industry

Ⅰ. Cars The emphasis of this sector must be on developing economy cars which meet the demands of safety, energy conservation, emission regulations and individual use, and further raising the share in the total auto production. Environmentally-friendly cars for use as taxis should also be developed. By the end of the Tenth Five-Year-Plan period, gasoline     

Survey of the development of Chinese heavy trucks

Heavy trucks are tending to get bigger. The market share for those with a payload of more than 10 tons is increasing rapidly while that for below 10 tons is contracting. The scramble for markets of the most-favored vehicle types is keenest among......