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.     

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.