The Elastic Cross-Influence between Two Quasi-Hertzian Contact Zones

SUMMARY

In this paper a type of contact between two bodies is considered, which leads to the formation of two separate contact zones. The contact zones considered are Hertzian ellipses if the distance between them is large. When the distance between the zones is finite it is necessary to take into account the elastic cross-influence of the two zones. In this paper an approximate method is proposed which allows the determination of the actual contact zones as Hertz's ellipses but without disregarding the cross-influence. Numerical results are presented for two, two-humped bodies pressed against each other and for steady state rolling of a flanged wheel along a steel rail.

It is shown that neglecting the cross influence for a flanged wheel leads to a negligible error in the contact forces but the forces are transmitted between wheel and rail through contact patches which are generally more slender than with the cross-influence neglected.     

Tangential force variation due to the bogie direction reversal procedure

A range of tangential forces is generated within the contact patch when a wheelset moves on the rail. These forces are intensified when incorporating curved tracks and motored axle rail vehicles [Arrus, P., de Pater, A.D. and Meyers, P., 2002, The stationary motion of a one-axle vehicle along a circular curve with real rail and wheel profiles. Vehicle System Dynamics, 37(1), 29–58]. The wheelset is subject to flange contact if an unbalanced force remains in a curve towards the high rail gauge face. The resultant force in the transverse direction includes the lateral force, the radial force, and the creep forces in addition to the effect of the frequent wheelset displacement due to the kinematic oscillation [Iwnicki, S., 2003, Simulation of wheel–rail contact forces. Fatigue Fracture Engineering Material Structure, 26, 887–900]. This article has focused on a potential variation in some of the forces cited when the wheelset is subject to backward and forward movements. A severe wear rate observed within the wheel flange region in Iranian Railways was investigated by operating a test bogie on a curvaceous track. An obvious improvement in the wear rate and wear pattern of the wheels was attained when the second test bogie encountered a bogie direction reversal procedure. This enhancement is considered in this article from the force analysis standpoint.     

Solution of the Multiple Wheel and Rail Contact Dynamic Problem

SUMMARY

An unconventional method for calculating the forces developing in the wheel and rail contact patches of a railway vehicle has been implemented at the New Technology Laboratory of INRETS. It takes into account the elastic deformations of the materials in the Hertzian elliptical contact areas; the possibility of having simultaneously several contact patches on each wheel, is introduced in the simulation of the dynamic phenomena.

The theory is applied for a high speed bogie running on a perfectly straight track.     

Stochastic analysis of dynamic interaction between train and railway turnout

The performance of a railway turnout (switch and crossing) is influenced by a large number of input parameters of the complex train–turnout system. To reach a robust design that performs well for different traffic situations, random distributions (scatter) of these inputs need to be accounted for in the design process. Stochastic analysis methods are integrated with a simulation model of the dynamic interaction between train and turnout. For a given nominal layout of the turnout, using design of experiments methodology and a two-level fractional factorial screening design, four parameters (axle load, wheel–rail friction coefficient, and wheel and rail profiles) are identified to be the most significant. These parameters are further investigated using a three-level full factorial design and stochastic analysis. The random distributions of transverse wheel profile and set of transverse rail profiles along the switch panel are accounted for by the Karhunen–Loève expansion technique. The influence of the random distributions of the input parameters on the statistical outputs of wheel–rail contact forces, wear and rolling contact fatigue is assessed using Latin hypercube sampling to generate a number of stochastic load realizations.     

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.     

Critical Speed and Limit Speed in Phenomena of Lateral Instability on Railway Vehicles

SUMMARY

A comparison between theoretical calculations on dynamic lateral behaviour of railway vehicles and experimental results shows quite a sizeable difference between the calculated critical speed and the actual speed at which side impact phenomena will repeatedly occur between wheel flange and rail (running speed limit), such impact speed being remarkably lower than calculated.

Another typical experimental aspect is that the running speed limit will considerably vary for the same vehicle depending on the test track conditions. Such difference is usually attributed to alterations of the wheel-rail contact surfaces, only.

This paper will discuss some concurrent causes which may prove far from negligible, such as the effects of track defects, an amplification of the dynamic lateral displacement between wheel and rail on approaching the critical speed, the track mechanical properties, and in particular the track lateral rigidity.

The influence of some geometrical factors typical of the wheel-rail contact, such as side clearance and linearized conicity, will also be discussed. The approach is based on the application of statistical methods to dynamic linear systems.     

The vertical and the longitudinal dynamic responses of the vehicle–track system to squat-type short wavelength irregularity

The squat, a kind of rolling contact fatigue occurring on the rail top, can excite the high-frequency vehicle–track interaction effectively due to its geometric deviations with a typical wavelength of 20–40 mm, leading to the accelerated deterioration of a track. In this work, a validated 3D transient finite element model is employed to calculate in the time domain the vertical and the longitudinal dynamic contact forces between the wheel and the rail caused by squats. The vehicle–track structure and the wheel–rail continua are both considered in order to include all the important eigencharacteristics of the system related to squats. By introducing the rotational and translational movements of the wheel, the transient wheel–rail rolling contact is solved in detail by a 3D frictional contact model integrated. The contact filter effect is considered automatically in the simulations by the finite size of the contact patch. The present work focuses on the influences of the length, width and depth of a light squat on the resulted dynamic contact forces, for which idealised defect models are used. The growth of a squat is also modelled to a certain extent by a series of defects with different dimensions. The results show that the system is mainly excited at two frequencies separately in the vertical and the longitudinal dynamics. Their superposition explains the typical appearance of mature squats. As a squat grows up, the magnitude of the excited vibration at the lower frequency increases faster than the one at the higher frequency.     

A Theory of Wheelset Forces for Two Point Contact Between Wheel and Rail

This paper describes a quasistatic theory of wheelset forces for an important practical case of the wheelset rolling when one of the wheels touches the rail in two contact zones. One of these zones lies on the tread and the other on the wheel flange. For such contact the specific problem of finding the distribution of forces between the tread and flange arises. The simultaneous frictional rolling contact problems for both contact zones have been described with Kalker×apos;s non-linear theory and wheelset equilibrium equations.

The numerical results presented are for an individual wheelset on straight track, the distribution of forces being described for a wide range of loading conditions. The influence of steering on the distribution of forces has also been presented.

This theory can be easily extended for quasistatic curving of railway vehicles and may assist wear studies for vehicles with worn wheels.     

An inverse design method for rail grinding profiles

Inspired by the optimisation design method for restoration of worn wheel profiles, an inverse design method based on optimal rail grinding profiles is presented in this paper. To improve grinding quality, vehicle dynamic performance is chosen as the main criteria and rolling radii difference function is selected as the key factor (also main target function) determining dynamic performance. Grinding material to be removed is chosen as the auxiliary target aimed at extending rail service life. Besides that, wheel–rail contact distribution is also taken into consideration as an auxiliary target preventing stress concentration and fatigue growth. By introducing certain presuppositions, all the design targets will form an inverse design problem. This problem can be solved using hybrid discrete numerical methods. Considering different grinding requirements, two examples of grinding profile design for straight and curved track will be discussed. Results show that the presented method is efficient and effective. Practical implementation has been carried out at several grinding sites in China.     

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.     

Prediction of Wheel/Rail Profile Wear

SUMMARY

The alteration in wheel and rail profiles due to wear involves considerable vehicle and track-maintenance costs, and influences the loading capacity of the rails, as well as the operation safety and riding comfort of the vehicles. In the past twenty years a vehicle dynamics, contact mechanics and tribology based research work has emerged which is also recently continuous in an international scale, and this research is more and more intensive. Parallel to the growing possibilities of computer based analyses, several algorithms and numerical procedures have been elaborated, as well as measurement based experiments have been carried out to establish the reliable prediction of wear-caused wheel and rail profile alterations and to maximise the mileage performance by selecting the optimum vehicle system parameters for running gears operating on a selected railway line or a whole network under specified -in general inherently stochastic - traffic conditions. This paper takes an attempt to introduce the extended sphere of problems of wheel and rail wear prediction, as well as the latest results reflecting the present state of the art.     

The Available Methods to Calculate the Wheel/Rail Forces in Non Hertzian Contact Patches and Rail Damaging

Abstract

In the present paper, the methods used in railway dynamics to take into account the non Hertzian multi contact patches occurring on a wheel, are summarised. Some application cases are presented as well as an explanation of the so called “squat” rail deterioration.     

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.     

The design of a hold-off device to improve the lateral comfort of rail vehicles

ABSTRACT

Rail vehicles negotiating curves or in crosswinds are subjected to high lateral forces which provoke high displacements of the lateral suspension. As these displacements need to be limited due to gauging restrictions these forces cause the lateral suspension to reach the bumpstops and consequently the passenger comfort is significantly jeopardized. The paper presents the design of a pneumatic system that allows limiting the lateral displacement during curve negotiation (hold-off device). It describes the different phases of the design process starting from the definition of requirements to be fulfilled. The main components and the effect of their characteristics on the overall performance of the centring system are studied, and completed with an experimental analysis of the centring system. Finally, the described methodology is applied to a typical high speed rail vehicle. The results prove that the concept of a centring system which uses the same technology and components that are used in rail vehicles for the pneumatic height control system of secondary suspensions is possible. This fact is particularly interesting as the market offers this kind of components and has proven their reliability during many hours of service therefore the new hold-off system will be based on in-service validated components.     

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.     

Dynamic Response of a Six-axle Locomotive to Random Track Inputs

SUMMARY

Spectral analysis techniques are employed to analyze the dynamic response of a six-axle locomotive on tangent track to vertical and lateral random track irregularities. The locomotive is represented by a thirty-nine (39) degrees of freedom model. A linear model is employed by considering small displacements, linear suspension elements and a linear theory for the wheel-rail interaction. Power spectral densities of displacements, velocities and accelerations and the statistical average frequencies of the system are obtained for each degree of freedom. Comparison of the calculated dominating frequencies with existing experimental values shows good agreement. The technique of spectral analysis is an effective tool for model validation, and for the determination of rail vehicle response to track irregularities. The probability functions for the response can be used as a measure for the ride quality of rail vehicles and for the study of fatigue damage of components.

     

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.