Modelling,validation and analysis of a three-dimensional railway vehicle–track system model with linear and nonlinear track properties in the presence of wheel flats

This paper presents dynamic contact loads at wheel–rail contact point in a three-dimensional railway vehicle–track model as well as dynamic response at vehicle–track component levels in the presence of wheel flats. The 17-degrees of freedom lumped mass vehicle is modelled as a full car body, two bogies and four wheelsets, whereas the railway track is modelled as two parallel Timoshenko beams periodically supported by lumped masses representing the sleepers. The rail beam is also supported by nonlinear spring and damper elements representing the railpad and ballast. In order to ensure the interactions between the railpads, a shear parameter beneath the rail beams has also been considered into the model. The wheel–rail contact is modelled using nonlinear Hertzian contact theory. In order to solve the coupled partial and ordinary differential equations of the vehicle–track system, modal analysis method is employed. Idealised Haversine wheel flats with the rounded corner are included in the wheel–rail contact model. The developed model is validated with the existing measured and analytical data available in the literature. The nonlinear model is then employed to investigate the wheel–rail impact forces that arise in the wheel–rail interface due to the presence of wheel flats. The validated model is further employed to investigate the dynamic responses of vehicle and track components in terms of displacement, velocity, and acceleration in the presence of single wheel flat.     

Influence of some vehicle and track parameters on the environmental vibrations induced by railway traffic

A study is performed on the influence of some typical railway vehicle and track parameters on the level of ground vibrations induced in the neighbourhood. The results are obtained from a previously validated simulation framework considering in a first step the vehicle/track subsystem and, in a second step, the response of the soil to the forces resulting from the first analysis. The vehicle is reduced to a simple vertical 3-dof model, corresponding to the superposition of the wheelset, the bogie and the car body. The rail is modelled as a succession of beam elements elastically supported by the sleepers, lying themselves on a flexible foundation representing the ballast and the subgrade. The connection between the wheels and the rails is realised through a non-linear Hertzian contact. The soil motion is obtained from a finite/infinite element model. The investigated vehicle parameters are its type (urban, high speed, freight, etc.) and its speed. For the track, the rail flexural stiffness, the railpad stiffness, the spacing between sleepers and the rail and sleeper masses are considered. In all cases, the parameter value range is defined from a bibliographic browsing. At the end, the paper proposes a table summarising the influence of each studied parameter on three indicators: the vehicle acceleration, the rail velocity and the soil velocity. It namely turns out that the vehicle has a serious influence on the vibration level and should be considered in prediction models.     

Influence of uneven rail irregularities on the dynamic response of the railway track using a three-dimensional model of the vehicle–track system

A mathematical model of the vehicle–track interaction is developed to investigate the coupled behaviour of vehicle–track system, in the presence of uneven irregularities at left/right rails. The railway vehicle is simplified as a 3D multi-rigid-body model, and the track is treated as the two parallel beams on a layered discrete support system. Besides the car-body, the bogies and the wheel sets, the sleepers are assumed to have roll degree of freedom, in order to simulate the in-plane rotation of the components. The wheel–rail interface is treated using a nonlinear Hertzian contact model, coupling the mathematical equations of the vehicle–track systems. The dynamic interaction of the entire system is numerically studied in time domain, employing Newmark's integration method. The track irregularity spectra of both the left/right rails are taken into account, as the inputs of dynamic excitations. The dynamic responses of the track system induced by such irregularities are obtained, particularly in terms of the vertical (bounce) and roll displacements. The numerical model of the present research is validated using several benchmark models reported in the literature, for both the smooth and unsmooth track conditions. Four sample profiles of the measured rail irregularities are considered as the case studies of excitation sources, examining their influences on the dynamic behaviour of the coupled system. The results of numerical simulations demonstrate that the motion of track system is significantly influenced by the presence of uneven irregularities in left/right rails. Dynamic response of the sleepers in the roll direction becomes more sensitive to the rail irregularities, as the unevenness severity of the parallel profiles (quantitative difference between left and right rail spectra) is increased. The severe geometric deformation of the track in the bounce–pitch–roll directions is mainly related to such profile unevenness (cross-level) in left/right rails.     

A linear complementarity method for the solution of vertical vehicle–track interaction

A new method is proposed for the solution of the vertical vehicle–track interaction including a separation between wheel and rail. The vehicle is modelled as a multi-body system using rigid bodies, and the track is treated as a three-layer beam model in which the rail is considered as an Euler-Bernoulli beam and both the sleepers and the ballast are represented by lumped masses. A linear complementarity formulation is directly established using a combination of the wheel–rail normal contact condition and the generalised-α method. This linear complementarity problem is solved using the Lemke algorithm, and the wheel–rail contact force can be obtained. Then the dynamic responses of the vehicle and the track are solved without iteration based on the generalised-α method. The same equations of motion for the vehicle and track are adopted at the different wheel–rail contact situations. This method can remove some restrictions, that is, time-dependent mass, damping and stiffness matrices of the coupled system, multiple equations of motion for the different contact situations and the effect of the contact stiffness. Numerical results demonstrate that the proposed method is effective for simulating the vehicle–track interaction including a separation between wheel and rail.     

Effect of curved track support failure on vehicle derailment

In order to investigate the effect of curved track support failure on railway vehicle derailment, a coupled vehicle–track dynamic model is put forward. In the model, the vehicle and the structure under rails are, respectively, modelled as a multi-body system, and the rail is modelled with a Timoshenko beam rested on the discrete sleepers. The lateral, vertical, and torsional deformations of the beam are taken into account. The model also considers the effect of the discrete support by sleepers on the coupling dynamics of the vehicle and track. The sleepers are assumed to move backward at a constant speed to simulate the vehicle running along the track at the same speed. In the calculation of the coupled vehicle and track dynamics, the normal forces of the wheels/rails are calculated using the Hertzian contact theory and their creep forces are determined with the nonlinear creep theory by Shen et al [Z.Y. Shen, J.K. Hedrick, and J.A. Elkins, A comparison of alternative creep-force models for rail vehicle dynamic analysis, Proceedings of the 8th IAVSD Symposium, Cambridge, MA, 1984, pp. 591–605]. The motion equations of the vehicle/track are solved by means of an explicit integration method. The failure of the components of the curved track is simulated by changing the track stiffness and damping along the track. The cases where zero to six supports of the curved rails fail are considered. The transient derailment coefficients are calculated. They are, respectively, the ratio of the wheel/rail lateral force to the vertical force and the wheel load reduction. The contact points of the wheels/rails are in detail analysed and used to evaluate the risk of the vehicle derailment. Also, the present work investigates the effect of friction coefficient, axle load and vehicle speed on the derailments under the condition of track failure. The numerical results obtained indicate that the failure of track supports has a great influence on the whole vehicle running safety.     

Fundamentals of vehicle–track coupled dynamics

This paper presents a framework to investigate the dynamics of overall vehicle–track systems with emphasis on theoretical modelling, numerical simulation and experimental validation. A three-dimensional vehicle–track coupled dynamics model is developed in which a typical railway passenger vehicle is modelled as a 35-degree-of-freedom multi-body system. A traditional ballasted track is modelled as two parallel continuous beams supported by a discrete-elastic foundation of three layers with sleepers and ballasts included. The non-ballasted slab track is modelled as two parallel continuous beams supported by a series of elastic rectangle plates on a viscoelastic foundation. The vehicle subsystem and the track subsystem are coupled through a wheel–rail spatial coupling model that considers rail vibrations in vertical, lateral and torsional directions. Random track irregularities expressed by track spectra are considered as system excitations by means of a time–frequency transformation technique. A fast explicit integration method is applied to solve the large nonlinear equations of motion of the system in the time domain. A computer program named TTISIM is developed to predict the vertical and lateral dynamic responses of the vehicle–track coupled system. The theoretical model is validated by full-scale field experiments, including the speed-up test on the Beijing–Qinhuangdao line and the high-speed running test on the Qinhuangdao–Shenyang line. Differences in the dynamic responses analysed by the vehicle–track coupled dynamics and by the classical vehicle dynamics are ascertained in the case of vehicles passing through curved tracks.     

On the effect of unsupported sleepers on the dynamic behaviour of a railway track

The effect of unsupported sleepers on the dynamic behaviour of a railway track is studied based on vehicle–track dynamic interaction theory, using a model of the track as a Timoshenko beam supported on a periodic elastic foundation. Considering the vehicle's running speed and the number of unsupported sleepers, the track dynamic characteristics are investigated and verified in the time and frequency domains by experiments on a 1:5 scale model wheel–rail test rig. The results show that when hanging sleepers are present, leading to a discontinuous and irregular track support, additional wheel–rail interaction forces are generated. These forces increase as further sleepers become unsupported and as the vehicle's running speed increases. The adjacent supports experience increased dynamic forces which will lead to further deterioration of track quality and the formation of long wavelength track irregularities, which worsen the vehicles’ running stability and riding comfort. Stationary transfer functions measurements of the dynamic behaviour of the track are also presented to support the findings.     

Prediction of the interaction between a simple moving vehicle and an infinite periodically supported rail – Green's functions approach

This paper herein describes the interaction between a simple moving vehicle and an infinite periodically supported rail, in order to signalise the basic features of the vehicle/track vibration behaviour in general, and wheel/rail vibration, in particular. The rail is modelled as an infinite Timoshenko beam resting on semi-sleepers via three-directional rail pads and ballast. The time-domain analysis was performed applying Green's matrix of the track method. This method allows taking into account the nonlinearities of the wheel/rail contact and the Doppler effect. The numerical analysis is dedicated to the wheel/rail response due to two types of excitation: the steady-state interaction and rail irregularities. The study points out to certain aspects regarding the parametric resonance, the amplitude-modulated vibration due to corrugation and the Doppler effect.     

Vertical random vibration analysis of vehicle–track coupled system using Green's function method

A vertical vehicle–track coupled dynamic model, consisting of a high-speed train on a continuously supported rail, is established in the frequency-domain. The solution is obtained efficiently by use of the Green's function method, which can determine the vibration response over a wide range of frequency without any limitations due to modal truncation. Moreover, real track irregularity spectra can be used conveniently as input. The effect of the flexibility of both track and car body on the entire vehicle–track coupled dynamic response is investigated. A multi-body model of a vehicle with either rigid or flexible car body is defined running on three kinds of track: a rigid rail, a track stiffness model and a Timoshenko beam model. The results show that neglecting the track flexibility leads to an overestimation of both the contact force and the whole vehicle vibration response. The car body flexibility affects the ride quality of the vehicle and the coupling through the track and can be significant in certain frequency ranges. Finally, the effect of railpad and ballast stiffness on the vehicle–track coupled vibration is analysed, indicating that the stiffness of the railpad has an influence on the system in a higher frequency range than the ballast.     

Simulation of vertical dynamic vehicle–track interaction using a two-dimensional slab track model

The vertical dynamic interaction between a railway vehicle and a slab track is simulated in the time domain using an extended state-space vector approach in combination with a complex-valued modal superposition technique for the linear, time-invariant and two-dimensional track model. Wheel–rail contact forces, bending moments in the concrete panel and load distributions on the supporting foundation are evaluated. Two generic slab track models including one or two layers of concrete slabs are presented. The upper layer containing the discrete slab panels is described by decoupled beams of finite length, while the lower layer is a continuous beam. Both the rail and concrete layers are modelled using Rayleigh–Timoshenko beam theory. Rail receptances for the two slab track models are compared with the receptance of a traditional ballasted track. The described procedure is demonstrated by two application examples involving: (i) the periodic response due to the rail seat passing frequency as influenced by the vehicle speed and a foundation stiffness gradient and (ii) the transient response due to a local rail irregularity (dipped welded joint).     

High-frequency random vibration analysis of a high-speed vehicle–track system with the frequency-dependent dynamic properties of rail pads using a hybrid SEM–SM method

A hybrid Spectral Element Method (SEM)–Symplectic Method(SM) method for high-efficiency computation of the high-frequency random vibrations of a high-speed vehicle–track system with the frequency-dependent dynamic properties of rail pads is presented. First, the Williams-Landel-Ferry (WLF) formula and Fractional Derivative Zener (FDZ) model were, respectively, applied for prediction and representation of the frequency-dependent dynamic properties of Vossloh 300 rail pads frequently used in China's high-speed railway. Then, the proposed hybrid SEM–SM method was used to investigate the influence of the frequency-dependent dynamic performance of Vossloh 300 rail pads on the high-frequency random vibrations of high-speed vehicle–track systems at various train speeds or different levels of rail surface roughness. The experimental results indicate that the storage stiffness and loss factors of Vossloh 300 rail pad increase with the decrease in dynamic loads or the increase in preloads within 0.1–10,000?Hz at 20°C, and basically linearly increase with frequency in a logarithmic coordinate system. The results computed by the hybrid SEM–SM method demonstrate that the frequency-dependent viscous damping of Vossloh 300 rail pads, compared with its constant viscous damping and frequency-dependent stiffness, has a much more conspicuous influence on the medium-frequency (i.e. 20–63?Hz) random vibrations of car bodies and rail fasteners, and on the mid- (i.e. 20–63?Hz) and high-frequency (i.e. 630–1250?Hz) random vibrations of bogies, wheels and rails, especially with the increase in train speeds or the deterioration of rail surface roughness. The two sensitive frequency bands can also be validated by frequency response function (FRF) analysis of the proposed infinite rail–fastener model. The mid and high frequencies influenced by the frequency-dependent viscous damping of rail pads are exactly the dominant frequencies of ground vibration acceleration and wheel rolling noise caused by high-speed railways, respectively. Even though the existing time-domain (or frequency-domain) finite track models associated with the time-domain (or frequency-domain) fractional derivative viscoelastic (FDV) models of rail pads can also be used to reach the same conclusions, the hybrid SEM–SM method in which only one element is required to compute the high-order vibration modes of infinite rail is more appropriate for high-efficiency analysis of the high-frequency random vibrations of high-speed vehicle–track systems.     

Distributed support modelling for vertical track dynamic analysis

The finite length nature of rail-pad supports is characterised by a Timoshenko beam element formulation over an elastic foundation, giving rise to the distributed support element. The new element is integrated into a vertical track model, which is solved in frequency and time domain. The developed formulation is obtained by solving the governing equations of a Timoshenko beam for this particular case. The interaction between sleeper and rail via the elastic connection is considered in an analytical, compact and efficient way. The modelling technique results in realistic amplitudes of the ‘pinned–pinned’ vibration mode and, additionally, it leads to a smooth evolution of the contact force temporal response and to reduced amplitudes of the rail vertical oscillation, as compared to the results from concentrated support models. Simulations are performed for both parametric and sinusoidal roughness excitation. The model of support proposed here is compared with a previous finite length model developed by other authors, coming to the conclusion that the proposed model gives accurate results at a reduced computational cost.     

Influence of long-wavelength track irregularities on the motion of a high-speed train

Vertical track irregularities over viaducts in high-speed rail systems could be possibly caused by concrete creep if pre-stressed concrete bridges are used. For bridge spans that are almost uniformly distributed, track irregularity exhibits a near-regular wave profile that excites car bodies as a high-speed train moves over the bridge system. A long-wavelength irregularity induces low-frequency excitation that may be close to the natural frequencies of the train suspension system, thereby causing significant vibration of the car body. This paper investigates the relationship between the levels of car vibration, bridge vibration, track irregularity, and the train speed. First, this study investigates the vibration levels of a high-speed train and bridge system using 3D finite-element (FE) transient dynamic analysis, before and after adjustment of vertical track irregularities by means of installing shimming plates under rail pads. The analysis models are validated by in situ measurements and on-board measurement. Parametric studies of car body vibration and bridge vibration under three different levels of track irregularity at five train speeds and over two bridge span lengths are conducted using the FE model. Finally, a discontinuous shimming pattern is proposed to avoid vehicle suspension resonance.     

Dynamic Interaction between Wheel and Track - A Parametric Search towards an Optimal Design of Rail Structures

The dynamic vertical interaction between a moving rigid wheel and a flexible railway track is investigated. A round and smooth wheel tread and an initially straight and noncorrugated rail surface are assumed in the present optimization study. A symmetric linear three-dimensional beam structure model of a finite portion of the track is suggested including rail, pads, sleepers and ballast with spatially nonproportional damping. The full interaction problem is numerically solved by use of an extended state-space vector approach in conjunction with a complex modal superposition for the track. Transient bending stresses in sleepers and rail are calculated. The influence of seven selected track parameters on the dynamic behaviour of the track is investigated. A two-level fractional factorial design method is used in the search for a combination of numerical levels of these parameters making the maximum bending stresses a minimum.     

Semi-active suspension systems for railway vehicles using magnetorheological dampers. Part I: system integration and modelling

In this paper, it is aimed to investigate semi-active suspension systems using magnetorheological (MR) fluid dampers for improving the ride quality of railway vehicles. A 17-degree-of-freedom (DOF) model of a full-scale railway vehicle integrated with the semi-active controlled MR fluid dampers in its secondary suspension system is proposed to cope with the lateral, yaw, and roll motions of the car body, trucks, and wheelsets. The governing equations combining the dynamics of the railway vehicle integrated with MR dampers in the suspension system and the dynamics of the rail track irregularities are developed and a linear quadratic Gaussian (LQG) control law using the acceleration feedback is adopted, in which the state variables are estimated from the measurable accelerations with a Kalman estimator. In order to evaluate the performances of the semi-active suspension systems based on MR dampers for railway vehicles, the random and periodical track irregularities are modelled with a uniform state-space formulation according to the testing data and incorporated into the governing equation of the railway vehicle integrated with the semi-active suspension system. Utilising the governing equations and the semi-active controller developed in this paper, the simulation and analysis are presented in Part II of this paper.     

Dynamic Interaction between Wheel and Track - A Parametric Search towards an Optimal Design of Rail Structures

SUMMARY

The dynamic vertical interaction between a moving rigid wheel and a flexible railway track is investigated. A round and smooth wheel tread and an initially straight and noncorrugated rail surface are assumed in the present optimization study. A symmetric linear three-dimensional beam structure model of a finite portion of the track is suggested including rail, pads, sleepers and ballast with spatially nonproportional damping. The full interaction problem is numerically solved by use of an extended state-space vector approach in conjunction with a complex modal superposition for the track. Transient bending stresses in sleepers and rail are calculated. The influence of seven selected track parameters on the dynamic behaviour of the track is investigated. A two-level fractional factorial design method is used in the search for a combination of numerical levels of these parameters making the maximum bending stresses a minimum.     

Experimental and numerical investigation of the effect of rail corrugation on the behaviour of rail fastenings

This paper presents the results of a detailed investigation of the effects of rail corrugation on the dynamic behaviour of metro rail fastenings, obtained from extensive experiments conducted on site and from simulations of train–track dynamics. The results of tests conducted with a metro train operating on corrugated tracks are presented and discussed first. A three-dimensional (3D) model of the metro train and a slab track was developed using multi-body dynamics modelling and the finite element method to simulate the effect of rail corrugation on the dynamic behaviour of rail fastenings. In the model, the metro train is modelled as a multi-rigid body system, and the slab track is modelled as a discrete elastic support system consisting of two Timoshenko beams for the rails, a 3D solid finite element (FE) model for the slabs, periodic discrete viscoelastic elements for the rail fastenings that connect the rails to the slabs, and uniformly viscoelastic elements for the subgrade beneath the slabs. The proposed train–track model was used to investigate the effects of rail corrugation on the dynamic behaviour of the metro track system and fastenings. An FE model for the rail fastenings was also developed and was used to calculate the stresses in the clips, some of which rupture under the excitation of rail corrugation. The results of the field experiments and dynamics simulations provide an insight into the root causes of the fracture of the clips, and several remedies are suggested for mitigating strong vibrations and failure of metro rail fastening systems.     

A Model Reduction Procedure for the Dynamic Analysis of Rail Vehicles Subjected to Linear Creep Forces

The set of differential equations governing the motion of an unrestrained coned wheelset travelling on a tangent section of track and acted upon by creep forces arising from the contact between wheel and rail are, in the terminology of numerical analysis, extremely "stiff". This stiffness can be attributed to the existence of two negative real eigenvalues in the solution of the eigenproblem associated with the linearized equations of motion. Compared with the two complex conjugate eigenvalues that complete this solution, the real eigenvalues have large magnitudes and necessitate that relatively. small timesteps be used in order to obtain an accurate numerical integration of the full set of equations of motion. However, by truncating the set of left and right eigenvectors to eliminate these real eigenvalues in a modal analysis of the wheelset, it was found that their contribution to the overall dynamic response is negligible. This same modal truncation approach was then applied to the sub-structured equations of motion for a simple rail vehicle system consisting of two wheelsets connected to a main body by linear springs and dampers. Essentially, the physical degrees of freedom for each wheelset substructure were replaced by a single complex coordinate obtained from the previous normal modes analysis. Using this model reduction procedure, accurate numerical results for the motion of the rail vehicle were generated several times faster than the results obtained by numerically integrating the full set of differential equations directly.     

Vibrations Generated by Rail Vehicles: A Mathematical Model in the Frequency Domain

Movement of railway vehicles generates mechanical vibrations of a wide range of frequency. Depending on track materials, dissipation in form of viscous and hysteretic damping is present, and stiffness depends on strain-rate. In a previous paper (Castellani et al., 1998), a mathematical model to describe track materials has been developed in the frequency domain. The present paper applies this model, and attempts an analytical formulation of vehicle-track and soil interaction in the frequency domain. Rail vibrations during the passage of a vehicle are generated by three families of forces: a) the weight of the moving vehicle, b) the inertial reaction of the vehicle under the effect of corrugations over an undeformable rail, and, c) the vehicle inertial forces due to displacements of the rail. The first two groups of forces do not depend on the rail displacement, and the related mathematical formulation is a simple problem of forces at a mobile point of application. Formulation of the vehicle inertial forces, related to the rail vibration, requires reference to the acceleration of the rail, as seen by an observer in motion with the vehicle itself. Moreover, it is necessary to express the equilibrium equation of two dynamic systems, the vehicle and the track, at a the movable point of contact. There is no straight numerical procedure to solve this equation in the frequency domain. In the paper two theoretical propositions (Fryba, 1988; Grassie et al., 1982) are revisited with reference to the effect of the transit of a single wheel. Fryba infers that, in the absence of corrugations, the forces c) are null. Grassie et al. (1982) present a mathematical formulation of the interaction between wheel and rail, at mobile point of contact. At each position, the interaction force is of impulsive type. They presume that for a corrugation of harmonic type, of wavelength ?, the wheel is subject to a harmonic motion, of the frequency f = V/?, where V is the wheel velocity. All other frequency components, due to the impulse, are disregarded. Both these assumptions are shown to be inconsistent from a theoretical point of view, however they suggest suitable approaches to the solution.