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.     

Study on the derailment behaviour of a railway wheelset with solid axles in a railway turnout

ABSTRACT

Dynamic wheel–rail interaction in railway turnouts is more complicated than on ordinary track. In order to evaluate the derailment behaviour of railway wheelsets in railway turnouts, this paper presents a study of dynamic wheel–rail interaction during a wheel flange climbs on the turnout rails, by applying the elasticity positioning wheelset model. A numerical model is established based on a coupled finite element method and multi-body dynamics, and applied to study the derailment behaviour of a railway wheelset in both the facing and trailing directions in a railway turnout, as well as dynamic wheel–turnout rail interaction during the wheel flange climbing on the turnout rails. The influence of the wheel–rail attack angle and the friction coefficient on the dynamic derailment behaviour is investigated through the proposed model. The results show that the derailment safety for a wheelset passing the railway turnout in facing direction is significantly lower than that for the trailing direction and the ordinary track. The possibility of derailment for the wheelset passing the railway turnout in facing and trailing directions at positive wheel–rail attack angles will increase with an increase in the attack angles, and the possibility of derailment can be reduced by decreasing the friction coefficient.     

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.     

Numerical investigation on effect of the relative motion of stock/switch rails on the load transfer distribution along the switch panel in high-speed railway turnout

In railway turnout, the stock rail and switch rail are separated to enable the vehicle changing among the tracks, and they are provided with different rail resilience level on the baseplate. Therefore, there will be vertical relative motion between stock/switch rails under the wheel loads, and the relative motion will affect consequentially the wheel–rail contact conditions. A method is developed to investigate the effect of the relative motion of stock/switch rails on the load transfer distribution along the switch panel in high-speed railway turnout. First, the rigid wheel–rail contact points of stock/switch rails are calculated based on the trace line method, and then the contact status is determined by the presented equations, finally, the distribution of wheel–rail contact forces of stock/switch rails is obtained based on the continuity of interface displacements and forces which using an approximate surface deformation method. Some parametric studies have been performed, such as the lateral displacement of wheel set, the vertical contact forces, the wheel profiles and the vertical stiffness of rail pad. The results of the parametric study are presented and discussed.     

Influence of wheel and rail profile shape on the initiation of rolling contact fatigue cracks at high axle loads

The influence of wheel and rail profile shape features on the initiation of rolling contact fatigue (RCF) cracks is evaluated based on the results of multi-body vehicle dynamics simulations. The damage index and surface fatigue index are used as two damage parameters to assess the influence of the different features. The damage parameters showed good agreement to one another and to in-field observations. The wheel and rail profile shape features showed a correlation to the predicted RCF damage. The RCF damage proved to be most sensitive to the position of hollow wear and thus bogie tracking. RCF initiation and crack growth can be reduced by eliminating unwanted shape features through maintenance and design and by improving bogie tracking.     

Development and validation of a wear model for the analysis of the wheel profile evolution in railway vehicles

The numerical wheel wear prediction in railway applications is of great importance for different aspects, such as the safety against vehicle instability and derailment, the planning of wheelset maintenance interventions and the design of an optimal wheel profile from the wear point of view. For these reasons, this paper presents a complete model aimed at the evaluation of the wheel wear and the wheel profile evolution by means of dynamic simulations, organised in two parts which interact with each other mutually: a vehicle's dynamic model and a model for the wear estimation. The first is a 3D multibody model of a railway vehicle implemented in SIMPACK?, a commercial software for the analysis of mechanical systems, where the wheel–rail interaction is entrusted to a C/C++user routine external to SIMPACK, in which the global contact model is implemented. In this regard, the research on the contact points between the wheel and the rail is based on an innovative algorithm developed by the authors in previous works, while normal and tangential forces in the contact patches are calculated according to Hertz's theory and Kalker's global theory, respectively. Due to the numerical efficiency of the global contact model, the multibody vehicle and the contact model interact directly online during the dynamic simulations.

The second is the wear model, written in the MATLAB® environment, mainly based on an experimental relationship between the frictional power developed at the wheel–rail interface and the amount of material removed by wear. Starting from a few outputs of the multibody simulations (position of contact points, contact forces and rigid creepages), it evaluates the local variables, such as the contact pressures and local creepages, using a local contact model (Kalker's FASTSIM algorithm). These data are then passed to another subsystem which evaluates, by means of the considered experimental relationship, both the material to be removed and its distribution along the wheel profile, obtaining the correspondent worn wheel geometry.

The wheel wear evolution is reproduced by dividing the overall chosen mileage to be simulated in discrete spatial steps: at each step, the dynamic simulations are performed by means of the 3D multibody model keeping the wheel profile constant, while the wheel geometry is updated through the wear model only at the end of the discrete step. Thus, the two parts of the whole model work alternately until the completion of the whole established mileage. Clearly, the choice of an appropriate step length is one of the most important aspects of the procedure and it directly affects the result accuracy and the required computational time to complete the analysis.

The whole model has been validated using experimental data relative to tests performed with the ALn 501 ‘Minuetto’ vehicle in service on the Aosta–Pre Saint Didier track; this work has been carried out thanks to a collaboration with Trenitalia S.p.A and Rete Ferroviaria Italiana, which have provided the necessary technical data and experimental results.     

Precision analysis and dynamic stability in the numerical solution of the two-dimensional wheel/rail tangential contact problem

ABSTRACT

In this paper the two-dimensional contact problem is analysed through different mesh topologies and strategies for approaching equations, namely; the collocation method, Galerkin, and the polynomial approach. The two-dimensional asymptotic problem (linear theory) associated with very small creepage (or infinite friction coefficient) is taken as a reference in order to analyse the numerical methods, and its solution is tackled in three different ways, namely steady-state problem, dynamic stability problem, and non-steady state problem in the frequency domain. In addition, two elastic displacements derivatives calculation methods are explored: analytic and finite differences. The results of this work establish the calculation conditions that are necessary to guarantee dynamic stability and the absence of numerical singularities, as well as the parameters for using the method that allows for maximum precision at the minimum computational cost to be reached.     

A study of the relationship between Magic Formula coefficients and tyre design attributes through finite element analysis

Pacejka's Magic Formula Tyre Model is widely used to represent force and moment characteristics in vehicle simulation studies meant to improve handling behaviour during steady-state cornering. The experimental technique required to determine this tyre model parameters is fairly involved and highly sophisticated. Also, total test facilities are not available in most countries. As force and moment characteristics are affected by tyre design attributes and tread patterns, manufacturing of separate tyres for each design alternative affects tyre development cycle time and economics significantly. The objective of this work is to identify the interactions among various tyre design attributes-cum-operating conditions and the Magic Formula coefficients. This objective is achieved by eliminating actual prototyping of tyres for various design alternatives as well as total experimentation on each tyre through simulation using finite element analysis. Mixed Lagrangian–Eulerian finite element technique, a specialized technique in ABAQUS, is used to simulate the steady-state cornering behaviour; it is also efficient and cost-effective. Predicted force and moment characteristics are represented as Magic Formula Tyre Model parameters through non-linear least-squares fit using MATLAB. Issues involved in the Magic Formula Tyre Model representation are also discussed. A detailed analysis is made to understand the influence of various design attributes and operating conditions on the Magic Formula parameters. Tread pattern, tread material properties, belt angle, inflation pressure, frictional behaviour at the tyre–road contact interface and their interactions are found to significantly influence vehicle-handling characteristics.