Simulation of Train-Turnout Coupled Dynamics Using a Multibody Simulation Software

With the advancements of computing power, multibody simulation (MBS) tool is used to study not only train dynamics but also more realistic phenomena such as train-track coupled dynamics. However, train-turnout coupled dynamics within MBS is still hard to be found. In this paper, a train-turnout coupled model methodology using a MBS tool GENSYS is presented. Dynamic track properties of a railway track are identified through numerical receptance test on a simple straight track model. After that, the identified dynamic track properties are adopted in a switch and crossing (turnout) to simulate train-turnout coupled dynamic interaction including parameters such as rail bending stiffness and sleeper mass variation along the turnout. .e train-turnout coupled dynamic interaction is compared to the dynamic interaction simulated from a widely accepted moving mass trainturnout model. It is observed that the vertical and lateral normal forces for the new train-turnout coupled model and the conventional moving mass train-turnout model are in good agreement. In addition, the new train-turnout coupled model can provide additional track dynamics results. It is concluded that the train-turnout coupled model can provide a more realistic trainturnout dynamic interaction compared to the moving mass train-turnout model.


Introduction
A good understanding of the train and track dynamic interaction is important to optimise design and maintenance of the track components.Switches and crossings are of special interest because of the increased loading that causes extra wear and high maintenance costs.Train and track dynamic simulations are used to study how train or track reacts when they are subjected to different conditions such as different axle loads, travelling speed, and track geometry especially at the design and development stages.rough train and track dynamic analysis, it is possible to evaluate how modifications in train and/or track affect several factors such as riding comfort, vehicle or track accelerations, derailment risk, rolling contact fatigue, rail stresses, and rail deflection.On top of that, numerical simulations reduce reliance on experimental campaigns that in most cases can be very expensive and time consuming.
Two types of numerical methods are most commonly used in the quest of simulating train and/or track dynamics, that is, multibody simulation (MBS) method and finite element method (FEM).Since the train-track complicated when it comes to dynamics, geometry, and material behaviour, it is necessary to make simplifications to be able to build the model and run the simulations in a reasonable time.e MBS method is known to be efficient to study the train dynamics, but the tradeoff is the track model that has to be simplified.In the FEM, a very detailed train-track model can be modelled including bodies' structural flexibility; however, the computational effort is much higher compared to MBS. erefore, FEM in railway field is often used only for very detailed component modelling [1], detailed wheel-rail contact modelling [2,3], and so on.Sometimes, FEM is also used for train-track simulation as a whole, but simplifications are often made to reduce the computational effort [4,5].
MBS is preferable to study train-track dynamic interaction as a whole, especially for the complex railway track component such as switch and crossing.MBS tool dedicated for train-track analysis has a wheel-rail contact module that is able to model complex wheel-rail kinematic behaviour and the corresponding creep forces.Some authors, for example, Kassa and Nielsen [6], Lagos et al. [7], Wan et al. [8], Pålsson [9], and Lau and Kassa [10], used MBS to simulate trainturnout dynamic interaction to evaluate the maximum contact forces along a railway turnout, to study the influence of the turnout geometry, to optimise crossing nose profile, to optimise switch geometry, and to study the effect of yaw in a switch and crossing, respectively.As reported by Esveld [11], the track substructure has a direct influence on the dynamic wheel load, dynamic track stiffness, and track roughness.However, modelling railway track in MBS has always been a challenge [12].Moreover, the train-turnout model used by the abovementioned authors is moving mass train-turnout model of which the track is represented by a lumped mass or mass layers coupled using spring-damper elements, and the track follows the train's wheels throughout the simulation.Such track model is able to capture dynamics up to 200 Hz.Some frequencies such as sleeper passing frequency and frequencies induced by the rail bending stiffness are not able to be modelled.Moreover, constant sleeper mass is used even though reality sleeper mass is not the same along switches and crossings.
In the present paper, a train-turnout coupled model using a MBS software GENSYS is presented.Dynamic track properties are first determined using a simple straight track model through the receptance test by referring to existing receptance functions recorded in the literature.After that, the determined track properties are applied to a turnout model.e train-turnout coupled model is simulated, and the corresponding dynamic interaction is compared to the dynamic interaction from a moving mass train-turnout model.Some track dynamics of the train-turnout coupled model are also reported.Benefits and drawbacks of the model are discussed in the paper.

Train Model
Train model presented in this paper is a Norwegian passenger train represented by seven bodies, that is, a carbody, two bogies, and four wheelsets.Each of the bodies has six degrees of freedom (DOFs), that is, three translations and three rotations apart from the wheelsets which have only five DOFs with longitudinal translation being constrained.e bodies are coupled to each other with spring-damper elements.Unlike conventional moving mass model, the train model in this paper is modelled in a moving coordinate system relative to the track fixed coordinate system.Figure 1 illustrates a schematic drawing of the train and the track model.

Track Modelling.
ere have been some attempts to model railway track using MBS, for example, in [13][14][15].One of the most popular ways of modelling railway track is modelling rail of the track as a beam element using the Euler-Bernoulli beam theory.State of the general dynamic beam equation, rail in this case, can be expressed by the Euler-Lagrange equation: where ρ is the mass density, A is the cross-sectional area, E is the elastic modulus, and I is the second moment of the area of beam's cross section.Since Euler-Bernoulli beam can only be supplied with the load and provide deflection in vertical direction, q(x, t) is the external load and w(x, t) is the deflection of the beam in the vertical direction at time t along the longitudinal position x. e external load q(x, t) is the excitation and the corresponding reaction generated by the train (wheel) on the top of the beam and is expressed as where N s is the number of supports, which in this case is modelled as a body rigidly connected under the beam with vertical translation and rotation in pitch direction unconstrained, R i is the ith support reaction force, δ(•) is the Dirac delta function, x is the current position, x i is the support position, N w is the number of wheels travelling over the beam, F j is the wheel-rail contact force under the jth wheel, H(∼) is the unit step function, t is the current time step, t j is the inbound time of the jth wheel at the beam, and Δt j is the time travelled of the jth wheel through the beam.Since Euler-Lagrange equation can be expressed as where C is the viscous damping coefficient.rough modal superposition, the vertical deflection w(x, t) can be expressed by the following equation: where A h (t) is the hth mode time coefficient of vertical deflection and S z (h, x) is the hth mode shape function of vertical deflection.Substituting (5) into (4) yields the following equation: where is the model frequency, ξ h � C/2ρAω h is the damping ratio, and m h � ρA is the model mass of the hth mode for a single span beam.
Note that when defining the beam in GENSYS, Euler-Bernoulli beam function is a massless beam and the beam mass is only considered by the rail masses connected rigidly under the beam. is is because these masses are needed to act as bodies that allow more masses such as sleeper masses and ballast masses to be connected underneath them to 2 Modelling and Simulation in Engineering represent a more realistic track, Figure 1.
erefore, the Euler-Lagrange equation as stated in (1) is only completed when the rail masses rigidly connected under the beam are defined.e rail masses are constrained in every direction apart from vertical translation and pitch rotation.One should be careful when selecting the moment of inertia of the rail masses in pitch direction as it has direct consequences to the bending and deflection of the beam.For instance, if the moment of inertia of the rail masses is too high, the beam will be restricted from bending or the bending and the deflection of the beam will be unrealistically high if the moment of inertia is too low.
For multispan continuous beam, beam elements are represented by finite beam elements via matrix coupling.For more details of the nite beam elements, refer to Ottosen and Petersson [16].If needed, for example, in switch and crossing, the bending sti ness of the beam can be varied sectionally along the track.e corresponding model frequencies, damping ratio, and model masses can be expressed by complex trigonometric functions.Such approach has previously been successfully implemented by Sun et al. [17] to simulate vertical dynamic of train-rail bridge interaction.

Track Property Validation.
One of the most crucial information for the track modelling is the dynamic track properties.Although there is existing information of track properties such as the one recorded in [18], it applies only for the moving mass track model.erefore, a numerical track receptance test is performed on a straight track model (without vehicle) to identify the dynamic properties of the track.
e straight track model for the identi cation consists of 100 rail masses which connects rigidly under two beams (50 on each side).Each pair of rail masses (left and right) is connected to a sleeper mass (a total of 50 sleeper masses) with spring-damper elements representing railpads.Under each sleeper, three spring-damper elements, left, middle, and right, are used to represent the ballast exibility connecting the sleeper to the rigid ground.e arrangement of the sleeper to rigid ground coupling is such that the left and the right spring-damper elements are always 350 mm from the edge of the sleeper and another one is always in the middle of the sleeper.e sleeper masses are constrained in every direction apart from the vertical translation and roll rotation.
Two receptance tests, one on the rail above a sleeper and another in between two sleepers, are performed.e corresponding receptance functions are adjusted to match the receptance functions as described in the literature [19] by means of tuning the track parameters such as railpad and ballast sti ness and damping properties.
Since there is no rail mass in the middle of the supports, an approximated method has to be used to obtain the receptance function when the excitation is applied in the middle of two sleepers.e displacement of the rail in the middle of two sleepers, δ mid , can be approximated using the pitch angle and displacement of the adjacent rail masses: where δ i and δ i+1 , and χ i and χ i+1 are the vertical displacements and the pitch angles of the masses which rigidly connect under the beam adjacent to the point of observation, respectively, x is the support distance, and λ is the basis function which can be expressed as where where x r is the relative longitudinal position along the beam and x 0 is the distance to the position from the mass on the left of which in the case of observation in the middle of two masses, x r 0.5, and λ 1 , λ 2 , λ 3 , and λ 4 have a value of 0.5, 0.125, 0.5, and −0.125, respectively.Figure 2 shows the receptance functions of the track with the properties as listed in Table 1.It can be seen that the three principal resonance frequencies as mentioned in [19] are captured.e rst resonance, which is regarded as track resonance (rails and sleepers vibrate on the ballast) which usually occurs in the frequency range of 50 to 300 Hz, is captured at the frequency around 129 Hz. e second resonance frequency, which is regarded as the rails bouncing on the railpads and usually  Due to the approximation method ( 7)-( 9), the magnitude of the receptance function on the rail between two sleepers is slightly lower compared to the receptance function on the rail above the sleeper. is is because the displacement and rotation of the rail masses adjacent to the point of observation do not always have the same magnitude, and the compromise of the magnitudes leads to slight reduction of the overall magnitude.In addition, the magnitude of the receptance at the pin-pin frequency when the receptance function is observed on the rail between the two sleepers is not as drastic as stated in the literature. is is because the approximation is based on the de ection and rotation of the two adjacent rail masses at the point of observation, and the projection is restricted by the rotation magnitude of the rail masses.
A modal analysis of the track as shown in Figure 3 is performed to illustrate the corresponding mode shapes of the resonance frequencies.e resonance frequencies from the modal analysis coincide well with the principal resonance frequencies computed from the track receptance tests despite the linearization of the track model in the modal analysis.erefore, it is concluded that the track properties are valid to be applied to a more sophisticated turnout model.
In addition, a receptance test is performed on a moving mass track model as shown in Figure 4 to achieve a comparable track model with the full-track model.e moving mass track model consists of two rail masses coupled to a sleeper mass with a spring-damper element, and the sleeper mass is connected to rigid ground with two spring-damper elements.Besides, the sleeper mass is also attached to a rigid wall in the lateral direction using a spring-damper element with the value of 30 MN/m and 270 kNs/m, respectively.Apart from the vertical translation, the rail masses and sleeper masses are constrained in every direction.

Modelling and Simulation in Engineering
It can be seen that the track receptance function of the moving mass track model when similar values in Table 1 are adopted appears nearly the same as the receptance functions from the full-track model.However, the pin-pin resonance in the moving mass track model is missing because of the absence of the beam element.Moreover, one should also note that when the moving mass model is connected to a train model, the second resonance frequency will not appear because the inertia of the rail masses will not be considered.

Turnout Model.
e turnout model is a right turn 60E1-760-1 : 15 (nominal rail profile 60E1, curve radius 760 m, and turnout angle 1 : 15) turnout without rail inclination.e variation of the rail profile in the turnout is accounted for by 30 rail cross sections along the turnout.e total length of the turnout model is 120 m. e first 60 m of the track is a straight track modelled to damp the initial transient effect and the rest is the turnout model.is means that the turnout model consists of 400 rail masses, 200 sleeper masses, and 1000 spring-damper elements for the coupling of the rails to the sleepers and to the rigid ground.
Unlike normal track, the rail bending stiffness is varied according to the rail profile in contact with the wheels along the turnout.However, due to the constraints of the Euler-Bernoulli beam equation, the bending stiffness can only be varied sectionally.e bending stiffness of the beam sections is approximated by the mean of the bending stiffness at the closest adjacent measured rail profiles according to the construction drawing.In addition, the sleeper properties such as mass, dimension, and moment of inertias are also varied along the turnout according to the construction drawing.
Since the rails are positioned differently on the sleepers along the switch and crossing for facing and diverging route, the rails have to be positioned accordingly in the model for different moving directions (only diverging route is considered in this work).e rest of the track properties follow the validated track properties as listed in Table 1. Figure 5 is a schematic plan view of the turnout model without the 60 m straight section.

Train-Turnout Coupling
For the coupling of the train and turnout model, a stiff springdamper system 141 MN/m and 30 kNs/m, respectively, is adopted for each wheel of the train connected to a neglected inertia fictitious mass.e fictitious masses are needed in order to extract wheel-rail profile contact information and any track irregularity during the time simulation, that is, offline lookup table is used in the model.e fictitious mass is connected to the Euler-Bernoulli beams by a very stiff "railrail" contact stiffness k c , which has a value of 5000 MN/m assuming the wheel and the rail do not lost contact during the simulation.
For the normal wheel-rail contact problem, Hertz contact theory is adopted.e creepages and creep forces are solved by using the FASTSIM algorithm by Kalker based on the simplified theory of rolling contact.Unlike moving mass track model, the inertia of the rail masses does not need to be neglected because the wheel of the train in the full-track model rolls over the track without carrying the track under the wheel.

Numerical Results
Two simulations have been carried out, train-turnout coupled model and moving mass turnout model.Similar train model and values up to wheel-rail contact are used for the both simulations.e train travelled through the turnout at the speed of 80 km/h.e switch panel of the turnout starts at 0 m, and the crossing panel starts at around 47 m.
e sampling frequency for both simulations is 1000 Hz. Figure 6 shows the vertical normal force of the wheels on the leading wheelset.It can be seen that the vertical normal force is very similar in general for the both models.However, variation can be observed in switch and crossing panels where the dynamic impact is higher.is is probably because the track reacted to the dynamic impact differently and eventually affects the magnitude of the vertical normal force.Similar phenomena can be observed in the lateral normal force as illustrated in Figure 7. is can also be explained by the similarity of the vertical and lateral normal forces along the turnout when the dynamic interaction is not as significant.Figure 8 shows the zoomed-in vertical normal force.Sleeper passing effect can also be simulated in the trainturnout coupled model but not the moving mass turnout model.
Simulating the train-turnout coupled dynamic using this approach, additional information in connection to the track such as the time history of the sleeper, force exerted on the sleeper, sleeper acceleration, and force exerted on the ballast can also be observed.Figure 9 illustrates the forces from the left and the right rail exerted on the sleeper numbers 2, 40, and 80 near the switch panel, closure panel, and crossing panel, respectively.All the four wheel passages can be observed in the sleeper number 2 and number 40.However, only first two wheel passages can be observed at the sleeper number 80 due to the length of the track and the trailing bogie did not travel through it before the end of the simulation.It can be seen that the distribution of the forces is 6 Modelling and Simulation in Engineering even for the rst two wheel passage on the sleeper number 2, but the force gradually distributed more to the left wheel when the rear wheels passed the sleeper because of curvature experienced by the train.e curvature e ect becomes more obvious at the sleeper number 40.In addition to the curvature e ect, the dynamic impact due to the crossing pro le in the sleeper number 80 can be observed.
Figure 10 shows the maximum force exerted on the sleepers throughout the turnout.It can be seen that the maximum force exerted on the left rail is higher than that on the right rail because of the curvature of the track.Besides, highest force exerted on the sleeper is in the switch panel (sleeper number 1 to 20) and the lowest force exerted on the sleeper is in the crossing panel (sleeper number 70 to 90).
is is because the bending sti ness of the rail is higher in crossing panel compared to the switch panel.
e high bending sti ness of the rail in crossing panel prevents the force from transferring down to the sleepers even though the dynamic impact is high in this region.Although the dynamic impact on the rail in the switch panel is not as signi cant as in the crossing panel, the maximum force exerted on the sleepers is higher compared to the force exerted on the sleeper in the crossing panel due to the relatively lower bending sti ness of the rail.Figure 11 illustrates the forces transfer from the sleepers down to the ballast through the left, middle, and the right spring-damper elements under sleeper 2, 40, and 80.It can be seen that the force under the sleepers on the right is higher than the middle and the left in general. is is because of the formation of the turnout in which the rails are placed gradually rightward on the sleepers along the diverging route.
is e ect is more obvious at the sleeper number 80, which is Modelling and Simulation in Engineering near the crossing panel.e force exerted on the ballast on the right is greater than the force exerted on the ballast in the middle even though the magnitude of the impact force on the left rail (nearer to the middle of the sleeper) above the sleeper number 80 is higher.In addition, the magnitude of the initial force exerted on the ballast is also di erent along the turnout due to the variation of the sleeper masses along the turnout.e variation of the sleeper masses will also have e ect towards the dynamic of the track in general. is is especially crucial if the condition under the track is of interest such as track settlement prediction.It is because the mass of the sleeper might a ect the magnitude of the dynamic vertical force exerted on the ballast and eventually a ect the accuracy of the track settlement prediction especially in the switch and crossing panels where the dynamic e ect is signi cant.
Figure 12 shows the maximum force exerted on the ballast along the turnout.e maximum force exerted on the ballast on the left, middle, and right started at around the same magnitude 30 kN.Due to the formation of the turnout as mentioned earlier, the maximum force exerted on the ballast on the right increases gradually and vice versa for the left.However, the maximum force exerted on the ballast in the middle does not uctuate as much as the left and the right.Coincidently, at the sleeper number 13, the maximum force exerted on the ballast for all three points is very similar because of the balance of curvature and turnout formation e ect.In addition, due to the change of the sleeper geometry back to the ordinary geometry started from sleeper number 93, a reversal  e ect of the left and the right maximum force exerted on the ballast due to the carbody roll can also be observed.
As a comparison, the force exerted on the sleeper and the ballast of the leading wheelset for moving mass turnout model is also plotted, Figure 13.It can be seen that the time history for the vertical normal force, the force exerted on the sleeper, and the force exerted on the ballast are very similar.
is is because the vertical normal force is transmitted down to the sleeper and to the ground without a ecting by the rail bending sti ness, discretisation of the track, and the formation of the turnout.e forces exerted on the ballast have a slightly higher magnitude compared to the forces exerted on the sleeper due to the existence of the sleeper mass.However, this variation is constant throughout the turnout due to the consistency of the sleeper mass in the moving mass model.

Discussion
In most of the cases, railway track is not the major focus in MBS train-track dynamic analysis.is is because MBS tools dedicated for train-track dynamic analysis are mainly used for the train dynamic analysis and the dynamics contributed by the track is relatively low, assuming that the track geometry is more or less consistent such as tangent and curve tracks.However, when a more sophisticated track structure such as switch and crossing is involved, the track e ects such as rail bending sti ness and sleeper mass become crucial in order to simulate both train and track dynamics more accurately.e advantages of modelling railway track in MBS using this approach are listed as below: (1) Train and track models can be modelled as a whole.
(2) Rail bending sti ness can be modelled and varied.
(3) Sleeper mass and track sti ness can be varied.(4) Sleeper passing e ect can be modelled.
(5) Dynamic interaction within the track can be computed and observed.
In addition, modelling railway track in MBS can make other type of analysis in the track, that is, settlement analysis, inclusion of ballast characteristic, hanging sleeper e ect, track geometry irregularity, and track component optimisation, possible without oversimplifying the train or the track.
e downside of this approach is that generally in MBS, the initiation time increases with the increase of number of bodies in the simulation.For a switch and crossing model, many additional bodies are needed; for example, in this paper, additional 400 rail and 200 sleeper masses are needed alongside with the seven bodies for the train model.e initiation of the train-track model took 15 minutes, and the simulation took 11 minutes for a computer with processor Intel ® Core ™ i7-4600 CPU @ 2.10 GHz (4 CPUs) compared to the moving mass turnout model which took 3 minutes in total to simulate.
For that, a separate algorithm has been written to avoid the repetitive initiation process in which the initiation of the train-track model is saved during the rst simulation, so that for the subsequent simulation using the same train-track model, the initiation of the model can be avoided.

Conclusion
In this paper, a train-turnout coupled dynamics has been simulated using a multibody simulation software GENSYS.Track parameters such as rail bending sti ness, rail bending sti ness variation, and sleeper mass variation along the turnout are considered in the model.e dynamic track properties of the turnout model such as rail masses' pitch moment of inertia, railpads' sti ness and damping, and ballast's sti ness and damping are validated beforehand by Modelling and Simulation in Engineering comparing the receptance functions from a straight track model to the receptance functions recorded in the literature.
After that, the results of the train-turnout coupled model are compared to the results of a moving mass train-turnout model.e following conclusions are drawn: (1) e vertical and lateral normal forces for the trainturnout coupled model are in good agreement with moving mass train-turnout model.(2) e full-turnout model has influence towards the vertical and lateral normal forces especially in the switch panel and crossing panel where the dynamic impact is higher.(3) e force exerted on the sleepers of the train-turnout coupled model is lower compared to the moving mass track model because of the rail bending stiffness.Besides, the magnitude of the force exerted on the sleepers decreases with the increase of the rail bending stiffness.(4) e force exerted on the ballast is position of the trainload dependent.For example, in the present paper, the force on the right side of the sleepers is greater along the turnout as the train is travelling through the diverging direction.(5) e time history of the individual track component in the train-track coupled model can be analysed but not the moving mass track model.
It is also concluded that such approach can potentially provide a better train-track interaction result compared to moving mass track model because of the consideration of the track parameters.e future work is to validate the trainturnout model with a real turnout model by receptance test or comparing the components, that is, sleeper or rail accelerations from the real turnout with the numerical model.

Figure 1 :
Figure 1: Schematic drawing of the train-track model (not to scale): (a) side view and (b) rear view.v: wheelset velocity, EI r : rail bending stiffness, second subscript l: left position, k r : railpad stiffness, second subscript r: right position, c r : railpad damping, second subscript m: middle position, k b : ballast stiffness, subscript i: item number, c b : ballast damping, M c : carbody mass, z c : carbody vertical displacement, J c : carbody mass inertia, z b : bogie vertical displacement, M b : bogie mass, z w : wheelset vertical displacement, J b : bogie mass inertia, z s : sleeper mass vertical displacement, M w : wheelset mass, y c : carbody lateral displacement, J w : wheelset mass inertia, y b : bogie lateral displacement, M f : fictitious rail mass, y w : wheelset lateral displacement, M r : rail mass, z B : rail (beam) vertical displacement, J r : rail mass inertia, z r : rail mass vertical displacement, M s : sleeper mass, φ c : carbody roll rotation, J s : sleeper mass inertia, φ b : bogie roll rotation, k s : secondary suspension stiffness, φ w : wheelset roll rotation, c s : secondary suspension damping, φ s : sleeper mass roll rotation, k p : primary suspension stiffness, χ c : carbody mass pitch rotation, c p : primary suspension damping, χ b : bogie mass pitch rotation, k w : wheel-fictitious rail contact stiffness, χ r : rail mass pitch rotation, c w : wheel-fictitious rail contact damping, λ sl : sleeper distance, k c : fictitious rail-rail contact stiffness, and λ bl : ballast support distance.