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, 7 (STRIPES). These methods attempt to estimate the non-elliptical contact patches with a discrete extension of the Hertz theory. As a continuation of 2, 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.     

The Kalker book of tables for non-Hertzian contact of wheel and rail

A new regularisation of non-elliptical contact patches has been introduced, which enables building the look-up table called by us the Kalker book of tables for non-Hertzian contact (KBTNH), which is a fast creep force generator that can be used by multibody dynamics system simulation programs. The non-elliptical contact patch is regularised by a simple double-elliptical contact region (SDEC). The SDEC region is especially suitable for regularisation of contact patches obtained with approximate non-Hertzian methods for solving the normal contact problem of wheel and rail. The new regularisation is suitable for wheels and rails with any profiles, including worn profiles.

The paper describes the new procedure of regularisation of the non-elliptical contact patch, the structure of the Kalker book of tables, and parameterisation of the independent variables of the tables and creep forces.

A moderate volume Kalker book of tables for SDEC region suitable for simulation of modern running gears has been computed in co-simulation of Matlab and program CONTACT.

To access the creep forces of the Kalker book of tables, the linear interpolation has been applied.

The creep forces obtained from KBTNH have been compared to those obtained by program CONTACT and FASTSIM algorithm. FASTSIM has been applied on both the contact ellipse and the SDEC contact patch. The comparison shows that KBTNH is in good agreement with CONTACT for a wide range of creepage condition and shapes of the contact patch, whereas the use of FASTSIM on the elliptical patch and SDEC may lead to significant deviations from the reference CONTACT solutions.

The computational cost of calling creep forces from KBTNH has been estimated by comparing CPU time of FASTSIM and KBTNH. The KBTNH is 7.8–51 times faster than FASTSIM working on 36–256 discretisation elements, respectively.

In the example of application, the KBTNH has been applied for curving simulations and results compared with those obtained with the creep force generator employing the elliptical regularisation. The results significantly differ, especially in predicted creepages, because the elliptical regularisation neglects generation of the longitudinal creep force by spin creepage.     

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.     

Some New Results in Rolling Contact

In the present paper three problems in the simplified theory of rolling contact are investigated. As to the first problem, three benchmark loadings, derived by Kalker in 1973 for Hertzian rollingcontact, are in existence. Each of these loadings gives rise to a value of the flexibility parameter of the simplified theory. These values are combined to a single, creepage dependent value of the flexibility, which appears to have an error of at most 10 to 15%. Secondly, the law of Coulomb is generalised by introducing two values of the coefficient of friction. The FASTSIM algorithm is adapted to that, and it is found that the traction, and hence the displacement, show a discontinuity inside the contact area. The discontinuity in the displacement is removed by introducing damping in the constitutive relations of the simplified theory. The damping constant is determined experimentally. When the damping coefficient decreases, the damped solution tends to the solution obtained directly without damping. This establishes the correctness of the latter, but it does not compare well with the complete theory as implemented by the program CONTACT. Thirdly and finally, it is shown that inertial effects may be neglected at speeds of around 100 km/h and also for much higher rolling velocities.     

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.     

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.     

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.     

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.     

Wheel life prediction model – an alternative to the FASTSIM algorithm for RCF

In this article, a wheel life prediction model considering wear and rolling contact fatigue (RCF) is developed and applied to a heavy-haul locomotive. For wear calculations, a methodology based on Archard's wear calculation theory is used. The simulated wear depth is compared with profile measurements within 100,000?km. For RCF, a shakedown-based theory is applied locally, using the FaStrip algorithm to estimate the tangential stresses instead of FASTSIM. The differences between the two algorithms on damage prediction models are studied. The running distance between the two reprofiling due to RCF is estimated based on a Wöhler-like relationship developed from laboratory test results from the literature and the Palmgren-Miner rule. The simulated crack locations and their angles are compared with a five-year field study. Calculations to study the effects of electro-dynamic braking, track gauge, harder wheel material and the increase of axle load on the wheel life are also carried out.     

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.     

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.     

Modification of the semi-Hertz wheel–rail contact method based on recalculating the virtual penetration value

ABSTRACT

Wheel–rail contact calculation is of vital importance in vehicle system dynamics. In the existing methods of wheel–rail contact calculation, the finite element method and Kalker’s CONTACT program, which are based on the complementary energy principle, are the two methods with accuracy recognised. However, because of its very slow calculation speed, it cannot meet the requirement of online calculation, so a variety of fast non-elliptic algorithms have been proposed. The semi-Hertz method, which is recognised for its great contributions to the fast wheel–rail contact calculation, is based on the concept of virtual penetration. The calculation of virtual penetration is crucial to evaluate the shape and normal pressure distribution of the contact patch. In practice, the virtual penetration is related to the curvature of the whole contact patch; however, the range of the contact patch is determined by the value of penetration. Such an interaction leads the calculation into a dead loop. In the semi-Hertz method, the penetration is calculated by the Hertz parameters of the initial contact point. Thus, the practical range of the method is limited. In this paper, a fast-iterative method for solving virtual penetration is proposed, and a reliable value of virtual penetration can be obtained under any lateral wheel–rail relative curvature variation with good stability and speed. The normal and tangential solutions are analysed with different methods in this paper.     

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 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.     

Development and validation of a wear model for the analysis of the wheel profile evolution in railway vehicles

The numerical wheel wear prediction in railway applications is of great importance for different aspects, such as the safety against vehicle instability and derailment, the planning of wheelset maintenance interventions and the design of an optimal wheel profile from the wear point of view. For these reasons, this paper presents a complete model aimed at the evaluation of the wheel wear and the wheel profile evolution by means of dynamic simulations, organised in two parts which interact with each other mutually: a vehicle's dynamic model and a model for the wear estimation. The first is a 3D multibody model of a railway vehicle implemented in SIMPACK?, a commercial software for the analysis of mechanical systems, where the wheel–rail interaction is entrusted to a C/C++user routine external to SIMPACK, in which the global contact model is implemented. In this regard, the research on the contact points between the wheel and the rail is based on an innovative algorithm developed by the authors in previous works, while normal and tangential forces in the contact patches are calculated according to Hertz's theory and Kalker's global theory, respectively. Due to the numerical efficiency of the global contact model, the multibody vehicle and the contact model interact directly online during the dynamic simulations.

The second is the wear model, written in the MATLAB® environment, mainly based on an experimental relationship between the frictional power developed at the wheel–rail interface and the amount of material removed by wear. Starting from a few outputs of the multibody simulations (position of contact points, contact forces and rigid creepages), it evaluates the local variables, such as the contact pressures and local creepages, using a local contact model (Kalker's FASTSIM algorithm). These data are then passed to another subsystem which evaluates, by means of the considered experimental relationship, both the material to be removed and its distribution along the wheel profile, obtaining the correspondent worn wheel geometry.

The wheel wear evolution is reproduced by dividing the overall chosen mileage to be simulated in discrete spatial steps: at each step, the dynamic simulations are performed by means of the 3D multibody model keeping the wheel profile constant, while the wheel geometry is updated through the wear model only at the end of the discrete step. Thus, the two parts of the whole model work alternately until the completion of the whole established mileage. Clearly, the choice of an appropriate step length is one of the most important aspects of the procedure and it directly affects the result accuracy and the required computational time to complete the analysis.

The whole model has been validated using experimental data relative to tests performed with the ALn 501 ‘Minuetto’ vehicle in service on the Aosta–Pre Saint Didier track; this work has been carried out thanks to a collaboration with Trenitalia S.p.A and Rete Ferroviaria Italiana, which have provided the necessary technical data and experimental results.     

Out-of-round railway wheels—a study of wheel polygonalization through simulation of three-dimensional wheel–rail interaction and wear

The objective of this study is to develop a tool for investigation of wheel tread polygonalization with radial irregularities including 1 to 20 wavelengths around the circumference of the wheel. Therefore, an existing multibody system model for simulation of general three-dimensional train–track interaction (accounting for frequencies up to several kHz) is extended with rolling contact mechanics according to FASTSIM. Furthermore, the model is also modified to allow for general wheel–rail profiles. The numerical model uses the concept of an iteration scheme including simulation of dynamic train–track interaction in the time domain coupled with a long-term wear model. A demonstration example including a bogie of a subway train travelling on a straight track is presented. In the example, an initial wheel out-of-roundness (OOR) is applied to the wheels. This irregularity is based on an amplitude spectrum derived from measurements on new wheels. Simulation results show that the most important wavelength-fixing mechanisms of the wheel OOR are (i) the vertical resonance of the coupled train–track system at approximately 40 Hz (the P2 resonance) and (ii) the frequency region including the lowest vertical track antiresonance at 165 Hz, where the dynamic track stiffness is high. Only a straight track is studied, but the model allows for asymmetric train motion on such a track.     

Out-of-round railway wheels—a study of wheel polygonalization through simulation of three-dimensional wheel-rail interaction and wear

The objective of this study is to develop a tool for investigation of wheel tread polygonalization with radial irregularities including 1 to 20 wavelengths around the circumference of the wheel. Therefore, an existing multibody system model for simulation of general three-dimensional train-track interaction (accounting for frequencies up to several kHz) is extended with rolling contact mechanics according to FASTSIM. Furthermore, the model is also modified to allow for general wheel-rail profiles. The numerical model uses the concept of an iteration scheme including simulation of dynamic train-track interaction in the time domain coupled with a long-term wear model. A demonstration example including a bogie of a subway train travelling on a straight track is presented. In the example, an initial wheel out-of-roundness (OOR) is applied to the wheels. This irregularity is based on an amplitude spectrum derived from measurements on new wheels. Simulation results show that the most important wavelength-fixing mechanisms of the wheel OOR are (i) the vertical resonance of the coupled train-track system at approximately 40 Hz (the P2 resonance) and (ii) the frequency region including the lowest vertical track antiresonance at 165 Hz, where the dynamic track stiffness is high. Only a straight track is studied, but the model allows for asymmetric train motion on such a track.