Dynamics of the railway track and the underlying soil: the boundary-element solution, theoretical results and their experimental verification

A combined finite-element boundary-element method is presented in detail to calculate the dynamic interaction of the railway track and the underlying soil. A number of results are shown for ballasted and slab track, demonstrating the influence of the stiffness of the soil and the rail pads on the vertical compliance of the track. The compliance of the track is combined with a simple model of the vehicle giving the transfer function of vehicle-track interaction. An experimental verification of the theoretical results is achieved by harmonic and impulse excitation with and without static (train-) load and by combined measurements of train-track-soil interaction. A clear vehicle-track resonance is found for the slab track with elastic rail pads and for higher frequencies at highspeed traffic, the dynamic axle loads due to sleeper passage are reduced.     

Effects of rail thermal stress on the dynamic response of vehicle and track

A new method is proposed to obtain the dynamic responses of the vehicle–track coupling system under the conditions of rail thermal stress changes in high-speed railways. Exact models are established with different rail longitudinal forces, in which multibody dynamic models are used for vehicles and the direct stiffness method for structures. In order to provide a general, simple and flexible formulation to express longitudinal stress distribution, the accurate model of long slab track consists of many small units with parameters which can be initialised separately. The exact analytical equation of track frequency and modal function was obtained by the transition matrix method, which can be used in calculating the dynamic response of wheel–rail coupling model. The proposed model is verified through comparisons with other classical solutions. Under the influence of train velocities and track irregularities, the specific vibration performances that frequency shifted and amplitude peak enhanced with thermal force are demonstrated through examples. The results show that the response analyses of vehicle and track have great application potentiality for fast estimation of the rail longitudinal stress.     

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).     

Effect of temperature- and frequency-dependent dynamic properties of rail pads on high-speed vehicle–track coupled vibrations

The soft under baseplate pad of WJ-8 rail fastener frequently used in China’s high-speed railways was taken as the study subject, and a laboratory test was performed to measure its temperature and frequency-dependent dynamic performance at 0.3?Hz and at ?60°C to 20°C with intervals of 2.5°C. Its higher frequency-dependent results at different temperatures were then further predicted based on the time–temperature superposition (TTS) and Williams–Landel–Ferry (WLF) formula. The fractional derivative Kelvin–Voigt (FDKV) model was used to represent the temperature- and frequency-dependent dynamic properties of the tested rail pad. By means of the FDKV model for rail pads and vehicle–track coupled dynamic theory, high-speed vehicle–track coupled vibrations due to temperature- and frequency-dependent dynamic properties of rail pads was investigated. Finally, further combining with the measured frequency-dependent dynamic performance of vehicle’s rubber primary suspension, the high-speed vehicle–track coupled vibration responses were discussed. It is found that the storage stiffness and loss factor of the tested rail pad are sensitive to low temperatures or high frequencies. The proposed FDKV model for the frequency-dependent storage stiffness and loss factors of the tested rail pad can basically meet the fitting precision, especially at ordinary temperatures. The numerical simulation results indicate that the vertical vibration levels of high-speed vehicle–track coupled systems calculated with the FDKV model for rail pads in time domain are higher than those calculated with the ordinary Kelvin–Voigt (KV) model for rail pads. Additionally, the temperature- and frequency-dependent dynamic properties of the tested rail pads would alter the vertical vibration acceleration levels (VALs) of the car body and bogie in 1/3 octave frequencies above 31.5?Hz, especially enlarge the vertical VALs of the wheel set and rail in 1/3 octave frequencies of 31.5–100?Hz and above 315?Hz, which are the dominant frequencies of ground vibration acceleration and rolling noise (or bridge noise) caused by high-speed railways respectively. Since the fractional derivative value of the adopted rubber primary suspension, unlike the tested rail pad, is very close to 1, its frequency-dependent dynamic performance has little effect on high-speed vehicle–track coupled vibration responses.     

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.     

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.     

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.     

Probabilistic assessment of railway vehicle-curved track systems considering track random irregularities

The randomness of track irregularities directly leads to the random vibration of the vehicle–track systems. To assess the dynamic performance of a railway system in more comprehensive and practical ways, a framework for probabilistic assessment of vehicle-curved track systems is developed by effectively integrating a vehicle–track coupled model (VTCM), a track irregularity probabilistic model (TIPM) with a probability density evolution method (PDEM). In VTCM, the railway vehicle and the curved track are coupled by the nonlinear wheel–rail interaction forces, and through TIPM, the ergodic properties of random track irregularities on amplitudes, wavelengths and probabilities can be properly considered in the dynamic calculations. Lastly, PDEM, a newly developed method for solving probabilistic transmissions between stochastic excitations and deterministic dynamic responses, is introduced to this probabilistic assessment model. Numerical examples validate the correctness and practicability of the proposed models. In this paper, the results of probabilistic assessment are presented to illustrate the dynamic behaviours of a high-speed railway vehicle subject to curved tracks with various radii, and to demonstrate the importance of considering the actual status of wheel–rail contacts and curve negotiation effects in vehicle-curved track interactions.     

A model for vehicle–track random interactions on effects of crosswinds and track irregularities

A stochastic mathematical model is developed to evaluate the dynamic behaviours and statistical responses of vehicle–track systems when random system excitations including crosswinds and track irregularities are imposed. In this model, the railway vehicle is regarded as a multi-rigid-body system, the track system is modelled by finite element theory. These two systems are spatially coupled by the nonlinear wheel–rail contact forces and unsteady aerodynamic forces. The high efficiency and accuracy of this stochastic model are validated by comparing to the robust Monte-Carlo method. Numerical studies show that crosswinds have a great influence on the dynamic performance of vehicle–track systems, especially on transverse vibrations. When the railway vehicle initially runs into the wind field, it will experience a severe vibration stage, and then stepping into a relatively steady state where the fluctuating winds and track irregularities will play deterministic roles in the deviations of system responses. Moreover, it is found that track irregularities should be properly considered in the safety assessment of the vehicle even in strong crosswinds.     

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.     

Establishment and verification of three-dimensional dynamic model for heavy-haul train–track coupled system

For the long heavy-haul train, the basic principles of the inter-vehicle interaction and train–track dynamic interaction are analysed firstly. Based on the theories of train longitudinal dynamics and vehicle–track coupled dynamics, a three-dimensional (3-D) dynamic model of the heavy-haul train–track coupled system is established through a modularised method. Specifically, this model includes the subsystems such as the train control, the vehicle, the wheel–rail relation and the line geometries. And for the calculation of the wheel–rail interaction force under the driving or braking conditions, the large creep phenomenon that may occur within the wheel–rail contact patch is considered. For the coupler and draft gear system, the coupler forces in three directions and the coupler lateral tilt angles in curves are calculated. Then, according to the characteristics of the long heavy-haul train, an efficient solving method is developed to improve the computational efficiency for such a large system. Some basic principles which should be followed in order to meet the requirement of calculation accuracy are determined. Finally, the 3-D train–track coupled model is verified by comparing the calculated results with the running test results. It is indicated that the proposed dynamic model could simulate the dynamic performance of the heavy-haul train well.     

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.     

Theoretical and experimental excitation force spectra for railway-induced ground vibration: vehicle–track–soil interaction,irregularities and soil measurements

Excitation force spectra are necessary for a realistic prediction of railway-induced ground vibration. The excitation forces cause the ground vibration and they are themselves a result of irregularities passed by the train. The methods of the related analyses – the wavenumber integration for the wave propagation in homogeneous or layered soils, the combined finite-element boundary-element method for the vehicle–track–soil interaction – have already been presented and are the base for the advanced topic of this contribution. This contribution determines excitation force spectra of railway traffic by two completely different methods. The forward analysis starts with vehicle, track and soil irregularities, which are taken from literature and axle-box measurements, calculates the vehicle–track interaction and gets theoretical force spectra as the result. The second method is a backward analysis from the measured ground vibration of railway traffic. A calculated or measured transfer function of the soil is used to determine the excitation force spectrum of the train. A number of measurements of different soils and different trains with different speeds are analysed in that way. Forward and backward analysis yield the same approximate force spectra with values around 1 kN for each axle and third of octave.     

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.     

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.     

On vehicle/track impact at connection between a floating slab and ballasted track and floating slab track's effectiveness of force isolation

When a vehicle runs over the connection between a floating slab track (FST) and ballasted track, wheel/rail impact may occur because of the stiffness difference in the two kinds of track, and thus a transition sector is usually included at the connection to smoothen the stiffness change. This phenomenon is studied by numerical simulation using a time-domain model for an idealised case without such a transition to determine whether it is actually necessary. Calculation results show that the wheel/rail impact load is moderate for a light FST and increases with the vehicle speed or decreasing the natural frequency of the FST. From simulation the wheel/rail parametric excitation is observed, as a result of variation in the stiffness of the FST with the period of the single slab length. The wheel/rail load due to the parametric excitation also increases with the vehicle speed. In addition, good performance of vibration isolation can be seen for the FST in terms of the force transmitted to the infrastructure.     

Extended applications of track irregularity probabilistic model and vehicle–slab track coupled model on dynamics of railway systems

Track irregularities are inevitably in a process of stochastic evolution due to the uncertainty and continuity of wheel–rail interactions. For depicting the dynamic behaviours of vehicle–track coupling system caused by track random irregularities thoroughly, it is a necessity to develop a track irregularity probabilistic model to simulate rail surface irregularities with ergodic properties on amplitudes, wavelengths and probabilities, and to build a three-dimensional vehicle–track coupled model by properly considering the wheel–rail nonlinear contact mechanisms. In the present study, the vehicle–track coupled model is programmed by combining finite element method with wheel–rail coupling model firstly. Then, in light of the capability of power spectral density (PSD) in characterising amplitudes and wavelengths of stationary random signals, a track irregularity probabilistic model is presented to reveal and simulate the whole characteristics of track irregularity PSD. Finally, extended applications from three aspects, that is, extreme analysis, reliability analysis and response relationships between dynamic indices, are conducted to the evaluation and application of the proposed models.     

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.     

An efficient approach to the analysis of rail surface irregularities accounting for dynamic train–track interaction and inelastic deformations

A two-dimensional computational model for assessment of rolling contact fatigue induced by discrete rail surface irregularities, especially in the context of so-called squats, is presented. Dynamic excitation in a wide frequency range is considered in computationally efficient time-domain simulations of high-frequency dynamic vehicle–track interaction accounting for transient non-Hertzian wheel–rail contact. Results from dynamic simulations are mapped onto a finite element model to resolve the cyclic, elastoplastic stress response in the rail. Ratcheting under multiple wheel passages is quantified. In addition, low cycle fatigue impact is quantified using the Jiang–Sehitoglu fatigue parameter. The functionality of the model is demonstrated by numerical examples.     

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.