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.     

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.     

Creep force modelling for rail traction vehicles based on the Fastsim algorithm

The evaluation of creep forces is a complex task and their calculation is a time-consuming process for multibody simulation (MBS). A methodology of creep forces modelling at large traction creepages has been proposed by Polach [Creep forces in simulations of traction vehicles running on adhesion limit. Wear. 2005;258:992–1000; Influence of locomotive tractive effort on the forces between wheel and rail. Veh Syst Dyn. 2001(Suppl);35:7–22] adapting his previously published algorithm [Polach O. A fast wheel–rail forces calculation computer code. Veh Syst Dyn. 1999(Suppl);33:728–739]. The most common method for creep force modelling used by software packages for MBS of running dynamics is the Fastsim algorithm by Kalker [A fast algorithm for the simplified theory of rolling contact. Veh Syst Dyn. 1982;11:1–13]. However, the Fastsim code has some limitations which do not allow modelling the creep force – creep characteristic in agreement with measurements for locomotives and other high-power traction vehicles, mainly for large traction creep at low-adhesion conditions. This paper describes a newly developed methodology based on a variable contact flexibility increasing with the ratio of the slip area to the area of adhesion. This variable contact flexibility is introduced in a modification of Kalker's code Fastsim by replacing the constant Kalker's reduction factor, widely used in MBS, by a variable reduction factor together with a slip-velocity-dependent friction coefficient decreasing with increasing global creepage. The proposed methodology is presented in this work and compared with measurements for different locomotives. The modification allows use of the well recognised Fastsim code for simulation of creep forces at large creepages in agreement with measurements without modifying the proven modelling methodology at small creepages.     

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.     

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.     

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.     

Simulation of dynamic interaction between train and railway turnout

Dynamic train–track interaction is more complex in railway turnouts (switches and crossings) than that in ordinary tangent or curved tracks. Multiple contacts between wheel and rail are common, and severe impact loads with broad frequency contents are induced, when nominal wheel–rail contact conditions are disturbed because of the continuous variation in rail profiles and the discontinuities in the crossing panel. The absence of transition curves at the entry and exit of the turnout, and the cant deficiency, leads to large wheel–rail contact forces and passenger discomfort when the train is switching into the turnout track. Two alternative multibody system (MBS) models of dynamic interaction between train and a standard turnout design are developed. The first model is derived using a commercial MBS software. The second model is based on a multibody dynamics formulation, which may account for the structural flexibility of train and track components (based on finite element models and coordinate reduction methods). The variation in rail profile is accounted for by sampling the cross-section of each rail at several positions along the turnout. Contact between the back of the wheel flange and the check rail, when the wheelset is steered through the crossing, is considered. Good agreement in results from the two models is observed when the track model is taken as rigid.     

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.     

Simulation of dynamic interaction between train and railway turnout

Dynamic train-track interaction is more complex in railway turnouts (switches and crossings) than that in ordinary tangent or curved tracks. Multiple contacts between wheel and rail are common, and severe impact loads with broad frequency contents are induced, when nominal wheel-rail contact conditions are disturbed because of the continuous variation in rail profiles and the discontinuities in the crossing panel. The absence of transition curves at the entry and exit of the turnout, and the cant deficiency, leads to large wheel-rail contact forces and passenger discomfort when the train is switching into the turnout track. Two alternative multibody system (MBS) models of dynamic interaction between train and a standard turnout design are developed. The first model is derived using a commercial MBS software. The second model is based on a multibody dynamics formulation, which may account for the structural flexibility of train and track components (based on finite element models and coordinate reduction methods). The variation in rail profile is accounted for by sampling the cross-section of each rail at several positions along the turnout. Contact between the back of the wheel flange and the check rail, when the wheelset is steered through the crossing, is considered. Good agreement in results from the two models is observed when the track model is taken as rigid.     

Determination of the wheel rail contact patch in semi-Hertzian conditions

This paper deals with the subject of the semi-Hertzian contact, which is a way to represent the wheel rail contact in railways or roller bearing applications. The method is based on the interpenetration of the two underformed bodies' profiles. The first step deals with the problem of the shape ratio; it is proposed to compensate the two main curvatures to obtain the good ratio in Hertzian conditions. Then, Hertz and Kalker's equations are used to establish the stresses at the level of a strip. These stresses expressions are used directly in a contact model discretized in strips and tabulated as a function of the lateral displacement between the wheel and the rail. The validation is made by comparison to the previous multi-Hertzian model of the VOCOLIN software. A first test shows identical results in Hertzian conditions; a second one shows only a small difference in semi-Hertzian conditions like S1002/UIC60 1:40.     

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.     

Determination of the wheel rail contact patch in semi-Hertzian conditions

This paper deals with the subject of the semi-Hertzian contact, which is a way to represent the wheel rail contact in railways or roller bearing applications. The method is based on the interpenetration of the two underformed bodies' profiles. The first step deals with the problem of the shape ratio; it is proposed to compensate the two main curvatures to obtain the good ratio in Hertzian conditions. Then, Hertz and Kalker's equations are used to establish the stresses at the level of a strip. These stresses expressions are used directly in a contact model discretized in strips and tabulated as a function of the lateral displacement between the wheel and the rail. The validation is made by comparison to the previous multi-Hertzian model of the VOCOLIN software. A first test shows identical results in Hertzian conditions; a second one shows only a small difference in semi-Hertzian conditions like S1002/UIC60 1:40.     

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.     

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.     

Railway ground vibrations induced by wheel and rail singular defects

Railway local irregularities are a growing source of ground-borne vibration and can cause negative environmental impacts, particularly in urban areas. Therefore, this paper analyses the effect of railway track singular defects (discontinuities) on ground vibration generation and propagation. A vehicle/track/soil numerical railway model is presented, capable of accurately predicting vibration levels. The prediction model is composed of a multibody vehicle model, a flexible track model and a finite/infinite element soil model. Firstly, analysis is undertaken to assess the ability of wheel/rail contact models to accurately simulate the force generation at the wheel/rail contact, in the presence of a singular defect. It is found that, although linear contact models are sufficient for modelling ground vibration on smooth tracks, when singular defects are present higher accuracy wheel/rail models are required. Furthermore, it is found that the variation in wheel/rail force during the singular defect contact depends on the track flexibility, and thus requires a fully coupled vehicle/track/foundation model. Next, a parametric study of ground vibrations generated by singular rail and wheel defects is undertaken. Six shapes of discontinuity are modelled, representing various defect types such as transition zones, switches, crossings, rail joints and wheel flats. The vehicle is modelled as an AM96 train set and it is found that ground vibration levels are highly sensitive to defect height, length and shape.     

Simulations of roughness growth on rails - results from a 2D non-Hertzian, non-steady contact model

A new model for simulating rail roughness growth on tangent track is presented in this paper. The model consists of three relatively independent components: (1) a time-domain vehicle/track interaction model; (2) a 2D non-Hertzian and non-steady wheel/rail contact model; and (3) a wear model. Wheel/rail contact forces for a given initial roughness obtained from the vehicle/track interaction model are used by the contact model to calculate the contact patch size, normal pressure and tangential stresses with material removal assumed to be linearly proportional to the friction work in the contact patch. The roughness profile is updated and fed back into vehicle/track interaction model. The 2D contact model is initially compared with a 3D model for various wavelength of initial sinusoidal roughness. Long term roughness growth is then simulated with the 2D contact model. Simulation shows that all initial sinusoidal roughness of wavelengths between 20-100 mm are levelled out. The wavelength-fixing mechanism, that has previously been used to explain the cause of corrugation, is not found in the present investigations.     

History of wheel/rail contact mechanics: from Redtenbacher to Kalker

The history of wheel/rail contact mechanics is an integrated part of contact mechanics. Problems of wheel/rail contact (damage phenomena and influence of contact mechanics on vehicle dynamics, especially on vehicle stability) have been investigated since the middle of the 19th Century. In 1855, Redtenbacher was the first to consider head checking. Stability investigations started with Boedecker (1887) and were continued by Carter (1916); however, they were not recognized at that time. Stability analysis was only considered to be an integrated step in vehicle design in about 1955 when stability problems started to occur in practice. The scientific foundations of our present investigations are based on Heinrich Hertz, Frederick William Carter and Hans Fromm. The key contributors to research in the last fifty years in the field of wheel/rail contact mechanics are Joost Kalker and Ken Johnson.     

A study on high-speed rolling contact between a wheel and a contaminated rail

A 3-D explicit finite element model is developed to investigate the transient wheel–rail rolling contact in the presence of rail contamination or short low adhesion zones (LAZs). A transient analysis is required because the wheel passes by a short LAZ very quickly, especially at high speeds. A surface-to-surface contact algorithm (by the penalty method) is employed to solve the frictional rolling contact between the wheel and the rail meshed by solid elements. The LAZ is simulated by a varying coefficient of friction along the rail. Different traction efforts and action of the traction control system triggered by the LAZ are simulated by applying a time-dependent driving torque to the wheel axle. Structural flexibilities of the vehicle–track system are considered properly. Analysis focuses on the contact forces, creepage, contact stresses and the derived frictional work and plastic deformation. It is found that the longitudinal contact force and the maximum surface shear stress in the contact patch become obviously lower in the LAZ and much higher as the wheel re-enters the dry rail section. Consequently, a higher wear rate and larger plastic flow are expected at the location where the dry contact starts to be rebuilt. In other words, contact surface damages such as wheel flats and rail burns may come into being because of the LAZ. Length of the LAZ, the traction level, etc. are varied. The results also show that local contact surface damages may still occur as the traction control system acts.