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.     

Indirect method for wheel–rail force measurement and derailment evaluation

Wheel set flange derailment criteria for railway vehicles are derived and the influence of wheel–rail contact parameters is studied. An indirect method for wheel–rail force measurement based on these derailment evaluation criteria is proposed. Laboratory tests for the calibration of strain–force devices on the bearing box are carried out to determine the relationship between the applied force and the measured strain. The simulation package, SIMPACK, is used to develop a passenger car model to generate wheel–rail forces and vibration signals. Different cases are considered in this model to provide an accurate validation of the identified wheel–rail forces. A feasibility test is conducted in the Beijing Loop test line using a passenger car equipped with a set of strain gauges on the wheel set. The comparison of the force time history applied to the instrumented wheel set and that obtained using the indirect method is presented.     

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

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

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.     

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.     

An analytical mathematical method for calculation of the dynamic wheel–rail impact force caused by wheel flat

The simplified method to determine a vertical impact force of wheel with flat and rail interaction is presented in this article. The presented simplified method can be used to identify maximum contact force and its distribution in the contact length between the damaged wheel and the rail. The vertical impact force depends on geometrical parameters of the rail and wheel with flat, speed of vehicle and the angle of deviation of rail. This article demonstrates the influence of wheel with flat geometrical parameters, speed of vehicle to maximum contact force and its distribution in the contact zone. The obtained values of the simplified method for determination of a vertical contact force are compared with the results obtained from field measurements.     

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.     

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.     

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.     

Dynamics of an isolated railway wheelset with conformal wheel–rail interactions

ABSTRACT

In multi-body dynamics, we model a non-conformal wheel/rail contact at one centre point since this contact is flat and Hertzian. However, the quasi-conformal contact requires more points since the contact is curved and non-Hertzian. The methodology in solving these points during dynamic simulations is the basis of this development. In this online effort, first, we present a general contact joint and the gutter search method of Pascal and Jourdan [The rigid-multi-Hertzian method as applied to conformal contacts. USA: ASME; 2007] in the context of a multi-body approach. Next, by adopting the non-iterative approach, a subset of these points with positive profile interpenetrations is selected to idealise one curved contact by a set of multi-Hertzian patches for which the Hertz normal contact solution is available. Finally, the feasibility of this multibody approach together with its implementation in two different codes is evaluated by simulating the motion of an isolated wheelset with realistic inertia. On tangent tracks, the wheelset with non-conformal pairing displays self-excited unstable oscillations while it displays a stable behaviour below a critical speed with conformal profiles. A further study reveals the net friction losses at multiple patches within the curved contact being the reason for the stable behavior.     

Estimation of wheel–rail contact conditions and adhesion using the multiple model approach

This paper presents the development of a multiple model estimation approach for the identification of the adhesion limit to overcome the problem of the wheel slip/slide at the rail wheel-rail contact. The contact characteristics at the rail wheel-rail interface are both highly nonlinear and subject to changes due to exposure to external contaminations. The detection of adhesion and its changes is therefore scientifically challenging, but would provide a critical information in the control of trains to avoid undesirable wear of the wheels/track but also the safety compromise of rail operations. This study exploits the variations in the dynamic behaviour of the railway wheelset caused by the contact condition changes and applies a bank of Kalman filters designed at selected operation points for the adhesion estimation. A fuzzy logic system is then developed to identify the contact conditions by examining the residuals from the Kalman filters.     

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.     

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.     

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.