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.     

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.     

Tyre rolling kinematics and prediction of tyre forces and moments: part I – theory and method

A new method to describe tyre rolling kinematics and how to calculate tyre forces and moments is presented. The Lagrange–Euler method is used to calculate the velocity and contact deformation of a tyre structure under large deformation. The calculation of structure deformation is based on the Lagrange method, while the Euler method is used to analyse the deformation and forces in the contact area. The method to predict tyre forces and moments is built using kinematic theory and nonlinear finite element analysis. A detailed analysis of the tyre tangential contact velocity and the relationships between contact forces, contact areas, lateral forces, and yaw and camber angles has been performed for specific tyres. Research on the parametric sensitivity of tyre lateral forces and self-aligning torque on tread stiffness and friction coefficients is carried out in the second part of this paper.     

Shear forces in the contact patch of a braked-racing tyre

This article identifies tyre modelling features that are fundamental to the accurate simulation of the shear forces in the contact patch of a steady-rolling, slipping and cambered racing tyre. The features investigated include contact patch shape, contact pressure distribution, carcass flexibility, rolling radius (RR) variations and friction coefficient. Using a previously described physical tyre model of modular nature, validated for static conditions, the influence of each feature on the shear forces generated is examined under different running conditions, including normal loads of 1500, 3000 and 4500 N, camber angles of 0° and?3°, and longitudinal slip ratios from 0 to?20%. Special attention is paid to heavy braking, in which context the aligning moment is of great interest in terms of its connection with the limit-handling feel. The results of the simulations reveal that true representations of the contact patch shape, carcass flexibility and lateral RR variation are essential for an accurate prediction of the distribution and the magnitude of the shear forces generated at the tread–road interface of the cambered tyre. Independent of the camber angle, the contact pressure distribution primarily influences the shear force distribution and the slip characteristics around the peak longitudinal force. At low brake-slip ratios, the friction coefficient affects the shear forces in terms of their distribution, while, at medium to high-slip ratios, the force magnitude is significantly affected. On the one hand, these findings help in the creation of efficient yet accurate tyre models. On the other hand, the research results allow improved understanding of how individual tyre components affect the generation of shear forces in the contact patch of a rolling and slipping tyre.     

Investigation of Features of Tyre Rolling at Non-Small Velocities on the Basis of a Simple Tyre Model with Distributed Mass Periphery

SUMMARY

On the basis of the brush-type tyre model the paper considers the interaction between steady-state rolling deformable wheel and flat road surface as well as corresponding force and moment characteristics of the wheel.

At least two zones of sliding, anisotropic dry friction, sliding friction coefficient speed-dependent and instantaneous leap of the friction coefficient when transition from sliding to adhesion zone occurs, have been taken into account, as well as distributed peripheral mass of tyre, elasticity, pseudo-dry friction and damping properties in radial, tangential and lateral directions of the elements at the wheel periphery, including a visco-elastic belt. Vertical force distribution in the contact area is not supposed to be known in advance and follows from the calculation. As a result, sliding zone lengths, distributed forces in contact area, six components of generalized road reaction reduced to the wheel center, and rolling resistance moment are found as functions of vertical load, movement velocity, longitudinal and side slip, friction in contact area with road, stiffnesses, dry friction and damping in the tyre model elements and of distributed peripheral mass.

A computer program developed in Fortran and results of calculations are of particular interest for qualitative analysis including steady rolling of studded tyre and also racing car and aircraft tyres which peripheral mass shows itself in a special way because of great movement velocities.     

Smoothing dry friction damping by dither generated in rolling contact of wheel and rail and its influence on ride dynamics of freight wagons

Friction dampers of mechanical systems are frequently exposed to medium-frequency (M-F) dither generated in the surrounding environment. A dithered system of technical importance is the railway freight wagon with friction dampers in the primary suspension developing two-dimensional friction, where dither is generated by the rolling contact of wheel and rail. This paper presents some results of the investigation of the influence of dither on dry friction damping. This influence has been studied experimentally and theoretically, and the parameters of dither influencing dry friction damping have been indicated. An experimental set-up has been designed that allows investigating friction damping in the presence of dither. The experiments have shown that friction damping in the presence of the M-F dither behaves like viscous damping. This means that dither smoothes dry friction as far as damping is concerned. To investigate this phenomenon theoretically, a rheological model of dry friction has been proposed that is applicable to one- and two-dimensional friction. In the latter case, the model takes into account friction anisotropy. According to performed numerical simulations of freight wagon motion, with dither supplied to the model through measured vertical accelerations of axle boxes, smoothing dry friction by dither strongly influences ride dynamics of the wagon with friction dampers in the primary suspension. Smoothing dry friction by dither should be accounted for in numerical simulations of motion of vehicles with friction dampers in the primary suspension by employing a proper model of the two-dimensional friction and application of realistic dither generated by rolling contact.     

Comparative study of vehicle tyre–road friction coefficient estimation with a novel cost-effective method

This paper qualitatively and quantitatively reviews and compares three typical tyre–road friction coefficient estimation methods, which are the slip slope method, individual tyre force estimation method and extended Kalman filter method, and then presents a new cost-effective tyre–road friction coefficient estimation method. Based on the qualitative analysis and the numerical comparisons, it is found that all of the three typical methods can successfully estimate the tyre force and friction coefficient in most of the test conditions, but the estimation performance is compromised for some of the methods during different simulation scenarios. In addition, all of these three methods need global positioning system (GPS) to measure the absolute velocity of a vehicle. To overcome the above-mentioned problem, a novel cost-effective estimation method is proposed in this paper. This method requires only the inputs of wheel angular velocity, traction/brake torque and longitudinal acceleration, which are all easy to be measured using available sensors installed in passenger vehicles. By using this method, the vehicle absolute velocity and slip ratio can be estimated by an improved nonlinear observer without using GPS, and the friction force and tyre–road friction coefficient can be obtained from the estimated vehicle velocity and slip ratio. Simulations are used to validate the effectiveness of the proposed estimation method.     

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.     

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.     

An optimal pole-matching observer design for estimating tyre–road friction force

In this paper, considering the dynamical model of tyre–road contacts, we design a nonlinear observer for the on-line estimation of tyre–road friction force using the average lumped LuGre model without any simplification. The design is the extension of a previously offered observer to allow a muchmore realistic estimation by considering the effect of the rolling resistance and a term related to the relative velocity in the observer. Our aim is not to introduce a new friction model, but to present a more accurate nonlinear observer for the assumed model. We derive linear matrix equality conditions to obtain an observer gain with minimum pole mismatch for the desired observer error dynamic system. We prove the convergence of the observer for the non-simplified model. Finally, we compare the performance of the proposed observer with that of the previously mentioned nonlinear observer, which shows significant improvement in the accuracy of estimation.     

Analysis of lateral tyre friction dynamics1

The tyre friction model is a key part of the overall multi-body tyre dynamics model. The LuGre dynamic tyre friction model is analytically linearised for pure cornering conditions. The linearised model parameters are conveniently expressed as functions of static curve slope parameters. The linearised lateral force and self-aligning torque submodels are described by equivalent mechanical systems. The linearised model and equivalent system parameters are analysed for different slip angle and wheel centre speed operating points. An example of the application of linearised tyre friction model to tyre vibration analysis is presented as well.     

Alternative friction models for braking train dynamics

Braking dynamics plays an essential role in ensuring the comfort and running safety of trains. Due to the complexity of the phenomena involved and the presence of dry friction, studies have been carried out by numerical methods using smoothening. Based on the non-smooth dynamics framework, in the present paper, an efficient and comprehensive alternative model is defined. Set-valued friction of Coulomb's law type is accounted for and motion equations are formulated as a differential inclusion. Some of the fundamental issues of contact dynamics are briefly reviewed. Static friction forces which arise in buffers are computed in a very intuitive and efficient manner, using an original method involving generalised inverses of matrices. The corresponding algorithm is described. Numerical integration is done by an event-driven algorithm. Indeterminate system configurations can be appropriately handled. The number of vehicles on which tests have been carried out ranged from 5 to 50. Comparison with the smoothening method reveals significant differences. Computational efficiency is dramatically improved. Computation speed can rise by a few orders of magnitude. Specific phenomena such as stick–slip or persistent longitudinal forces can be evidenced. Computed values of the system states differ substantially. In addition to this, any wagon connection model can be embedded.     

Frequency-Dependent Creep Coefficients for Three-Dimensional Rolling Contact Problems

Kalker's creep coefficients for linear rolling contact problems are only valid in the steady state case. A method for the extension of linear contact mechanics into the high frequency range is presented. Frequency dependent creep coefficients are obtained, which can be interpreted as harmonic responses of the creep forces for harmonically varying creepages. The influence of different constant longitudinal creepages and of different shapes of the contact area on the harmonic responses are discussed. The analysis is based on Kalker [5]     

Derivation of Frequency Dependent Creep Coefficients Based on an Elastic Half-Space Model

Kalker's creep coefficients for linear rolling contact problems are only valid in the steady state case. A method for the extension of linear contact mechanics into the high frequency range is presented. Frequency dependent creep coefficients are obtained, which can be interpreted as harmonic responses of the creep forces for harmonically varying creepages. For linear investigations of contact phenomena in the high frequency range simple laws governing the Kalker coefficients C₁₁and C₂₂ are given.     

Introduction of a friction coefficient dependent on the slip in the FastSim algorithm

One of the main limitations of algorithms relating forces and creepages at the wheel/rail contact is the use of a friction coefficient independent of the slip. This paper overcomes this limitation through a modification of the FastSim algorithm (based on the Simplified Theory of Kalker). A friction law based on the local value of the slip is established and the required formulation of the local slip elsewhere in the contact area is presented. Some difficulties of the method and the solutions adopted by the authors are also presented. Finally, the achieved improvements are shown through comparison of the results obtained both with the original and the modified FastSim algorithms.     

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.     

A vehicle ABS adaptive sliding-mode control algorithm based on the vehicle velocity estimation and tyre/road friction coefficient estimations

A sliding-mode observer is designed to estimate the vehicle velocity with the measured vehicle acceleration, the wheel speeds and the braking torques. Based on the Burckhardt tyre model, the extended Kalman filter is designed to estimate the parameters of the Burckhardt model with the estimated vehicle velocity, the measured wheel speeds and the vehicle acceleration. According to the estimated parameters of the Burckhardt tyre model, the tyre/road friction coefficients and the optimal slip ratios are calculated. A vehicle adaptive sliding-mode control (SMC) algorithm is presented with the estimated vehicle velocity, the tyre/road friction coefficients and the optimal slip ratios. And the adjustment method of the sliding-mode gain factors is discussed. Based on the adaptive SMC algorithm, a vehicle's antilock braking system (ABS) control system model is built with the Simulink Toolbox. Under the single-road condition as well as the different road conditions, the performance of the vehicle ABS system is simulated with the vehicle velocity observer, the tyre/road friction coefficient estimator and the adaptive SMC algorithm. The results indicate that the estimated errors of the vehicle velocity and the tyre/road friction coefficients are acceptable and the vehicle ABS adaptive SMC algorithm is effective. So the proposed adaptive SMC algorithm can be used to control the vehicle ABS without the information of the vehicle velocity and the road conditions.     

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.     

Estimation of the tyre–road maximum friction coefficient and slip slope based on a novel tyre model

In this article, a new approach to estimate the vehicle tyre forces, tyre–road maximum friction coefficient, and slip slope is presented. Contrary to the majority of the previous work on this subject, a new tyre model for the estimation of the tyre–road interface characterisation is proposed. First, the tyre model is built and compared with those of Pacejka, Dugoff, and one other tyre model. Then, based on a vehicle model that uses four degrees of freedom, an extended Kalman filter (EKF) method is designed to estimate the vehicle motion and tyre forces. The shortcomings of force estimation are discussed in this article. Based on the proposed tyre model and the improved force measurements, another EKF is implemented to estimate the tyre model parameters, including the maximum friction coefficient, slip slope, etc. The tyre forces are accurately obtained simultaneously. Finally, very promising results have been achieved for pure acceleration/braking for varying road conditions, both in pure steering and combined manoeuvre simulations.     

Dynamic tyre friction models for combined longitudinal and lateral vehicle motion

An extension to the LuGre dynamic friction model from longitudinal to longitudinal/lateral motion is developed in this paper. Application of this model to a tyre yields a pair of partial differential equations that model the tyre-road contact forces and aligning moment. A comparison of the steady-state behaviour of the dynamic model with existing static tyre friction models is presented. This comparison allows one to determine realistic values of the parameters for the new dynamic model. Via the introduction of a set of mean states we reduce the partial differential equations to a lumped model governed by a set of three ordinary differential equations. Such a lumped form describes the aggregate effect of the friction forces and moments and it can be useful for control design and online estimation. A method to incorporate wheel rim rotation is also proposed. The proposed model is evaluated by comparing both its steady-state as well as its dynamic characteristics via numerical simulations. The results of the simulations corroborate steady-state and dynamic/transient tyre characteristics found in the literature.