Computation of wheel–rail contact force for non-mapping wheel–rail profile of Translohr tram

Translohr tram has steel wheels, in V-like arrangements, as guide wheels. These operate over the guide rails in inverted-V arrangements. However, the horizontal and vertical coordinates of the guide wheels and guide rails are not always mapped one-to-one. In this study, a simplified elastic method is proposed in order to calculate the contact points between the wheels and the rails. By transforming the coordinates, the non-mapping geometric relationship between wheel and rail is converted into a mapping relationship. Considering the Translohr tram’s multi-point contact between the guide wheel and the guide rail, the elastic-contact hypothesis take into account the existence of contact patches between the bodies, and the location of the contact points is calculated using a simplified elastic method. In order to speed up the calculation, a multi-dimensional contact table is generated, enabling the use of simulation for Translohr tram running on curvatures with different radii.     

A dynamic wheel–rail impact analysis of railway track under wheel flat by finite element analysis

Wheel–rail interaction is one of the most important research topics in railway engineering. It involves track impact response, track vibration and track safety. Track structure failures caused by wheel–rail impact forces can lead to significant economic loss for track owners through damage to rails and to the sleepers beneath. Wheel–rail impact forces occur because of imperfections in the wheels or rails such as wheel flats, irregular wheel profiles, rail corrugations and differences in the heights of rails connected at a welded joint. A wheel flat can cause a large dynamic impact force as well as a forced vibration with a high frequency, which can cause damage to the track structure. In the present work, a three-dimensional finite element (FE) model for the impact analysis induced by the wheel flat is developed by the use of the FE analysis (FEA) software package ANSYS and validated by another validated simulation. The effect of wheel flats on impact forces is thoroughly investigated. It is found that the presence of a wheel flat will significantly increase the dynamic impact force on both rail and sleeper. The impact force will monotonically increase with the size of wheel flats. The relationships between the impact force and the wheel flat size are explored from this FEA and they are important for track engineers to improve their understanding of the design and maintenance of the track system.     

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.     

A new method for wheel–rail contact force continuous measurement using instrumented wheelset

With the aim of improving the continuous measurement of wheel–rail contact force by instrumented wheelset, instead of solving the non-linear equations, we proposed a new method based on state space theory. With this new method, the wheel–rail contact force can be calculated by the recurrence relation and the signals from strain gauge bridges on wheel web. The implementation of continuous instrumented wheelset is quite general and simplified, due to the specific bridging scheme is not necessary. It means that continuous measurement of the contact force could be realised with a simple bridging scheme, even as simple as discrete instrumented scheme. In this work, we first demonstrated and discussed the effectiveness and accuracy of this new method by estimation results with the numerical simulations, and we also applied this new method to two field tests, including one was conducted in a loop test line using a high-speed train and the other one was conducted in an urban line with a light rail vehicle. In a word, this new method is proved to be an effective way to monitor the wheel–rail contact force of rail vehicle track system.     

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.     

Sensitivity of the wheel–rail contact interactions and Dang Van Fatigue Index in the rail with respect to irregularities of the track geometry

This paper investigates the effects of the track geometry irregularities on the wheel–rail dynamic interactions and the rail fatigue initiation through the application of the Dang Van criterion, that supposes an elastic shakedown of the structure. The irregularities are modelled, using experimental data, as a stochastic field which is representative of the considered railway network. The tracks thus generated are introduced as the input of a railway dynamics software to characterise the stochastic contact patch and the parameters on which it depends: contact forces and wheelset–rail relative position. A variance-based global sensitivity analysis is performed on quantities of interest representative of the dynamic behaviour of the system, with respect to the stochastic geometry irregularities and for different curve radius classes and operating conditions. The estimation of the internal stresses and the fatigue index being more time-consuming than the dynamical simulations, the sensitivity analysis is performed through a metamodel, whose input parameters are the wheel–rail relative position and velocity. The coefficient of variation of the number of fatigue cycles, when the simulations are performed with random geometry irregularities, varies between 0.13 and 0.28. In a large radius curve, the most influent irregularity is the horizontal curvature, while, in a tight curve, the gauge becomes more important.     

Non-steady state modelling of wheel–rail contact problem

Among all the algorithms to solve the wheel–rail contact problem, Kalker's FastSim has become the most useful computation tool since it combines a low computational cost and enough precision for most of the typical railway dynamics problems. However, some types of dynamic problems require the use of a non-steady state analysis. Alonso and Giménez developed a non-stationary method based on FastSim, which provides both, sufficiently accurate results and a low computational cost. However, it presents some limitations; the method is developed for one time-dependent creepage and its accuracy for varying normal forces has not been checked. This article presents the required changes in order to deal with both problems and compares its results with those given by Kalker's Variational Method for rolling contact.     

Estimation of wheel–rail friction for vehicle certification

In certification of new rail vehicles with respect to running characteristics, a wide variety of operating conditions needs to be considered. However, in associated test runs the wheel–rail friction condition is difficult to handle because the friction coefficient needs to be fairly high and the friction is also generally hard to assess. This is an issue that has been studied in the European project DynoTRAIN and part of the results is presented in this paper. More specifically, an algorithm for estimating the wheel–rail friction coefficient at vehicle certification tests is proposed. Owing to lack of some measurement results, the algorithm here is evaluated in a simulation environment which is also an important step towards practical implementation. A quality measure of the friction estimate is suggested in terms of estimated wheel–rail spin and total creep. It is concluded that, tentatively, the total creep should exceed 0.006 and the spin should be less than 1.0 m^?1 for the algorithm to give a good friction estimate. Sensitivity analysis is carried out to imitate measurement errors, but should be expanded in further work.     

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.     

Safety of a railway wheelset – derailment simulation with increasing lateral force

The motivation for this research is to make a comparison between dynamic results of a free railway wheelset derailment and safety limits. For this purpose, a numerical simulation of a wheelset derailment submitted to increasing lateral force is used to compare with the safety limit, using different criteria. A simplified wheelset model is used to simulate derailments with different adhesion conditions. The contact force components, including the longitudinal and spin effects, are identified in a steady-state condition on the verge of a derailment. The contact force ratios are used in a three-dimensional (3D) analytical formula to calculate the safety limits. Simulation results obtained with two contact methods were compared with the published results and the safety limit was identified with the two criteria. Results confirm Nadal’s conservative aspect and show that safety 3D analytical formula presents slightly higher safety limits for lower friction coefficients and smaller limits for high friction, in comparison with the simulation results with Fastsim.     

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.     

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

A thorough investigation of wheel–rail impact due to wheel flats is presented, together with a quantitative characterization of the main mechanisms and parameters. A criterion for the speed with respect to contact loss between wheel and rail is derived. In the subcritical speed regime, the magnitude of the impact force is shown to be directly related to the geometry of the flat, whereas in the transcritical speed regime a fictitious or effective flat depth exists, which decreases with the second order of the speed. In this domain, the position of impact shifts towards the end of the flat with increasing speed. The impact force increases with the second order of the speed in the subcritical speed regime and approximately the first order of the speed in the transcritical speed regime. The magnitude of the impact force is inversely proportional to the minimum circumferential curvature of the wheel tread defect in the subcritical speed regime, and to the effective minimum curvature in the transcritical case. The variation with the flat depth is less; the impact varies with the square root of the flat depth. The presented theory is in accordance with measurements reported in the literature and explains characteristic features in them.     

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.     

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.     

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.     

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.     

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.     

Dynamic model for the wheel–rail contact friction

Accurately estimating the coefficient of friction (CoF) is essential in modelling railroad dynamics, reducing maintenance costs, and increasing safety in rail operations. The typical assumption of a constant CoF is widely used in theoretical studies; however, it has been noticed that the CoF is not constant, but rather depends on various dynamic parameters and instantaneous conditions. In this paper, we present a newly developed three-dimensional nonlinear CoF model for the dry rail condition and test the CoF variation using this model with estimated dynamic parameters. The wheel–rail is modelled as a mass–spring–damper system to simulate the basic wheel–rail dynamics. Although relatively simple, this model is considered sufficient for the purpose of this study. Simulations are performed at a train speed of 20 m/s using rail roughness as an excitation source. The model captures the CoF extremes and illustrates its nonlinear behaviour and instantaneous dependence on several structural and dynamic parameters.     

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.