Next Article in Journal
Evaluation of English–Slovak Neural and Statistical Machine Translation
Previous Article in Journal
Bio-Inspired Morphing Tail for Flapping-Wings Aerial Robots Using Macro Fiber Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Case Report

Prediction of Theoretical Derailments Caused by Cross-Winds with Frequency

1
Department of Rolling Stock System Engineering, Seoul National University of Science and Technology, Seoul 01811, Korea
2
Department of Railway Safety Engineering, Seoul National University of Science and Technology, Seoul 01811, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(7), 2947; https://doi.org/10.3390/app11072947
Submission received: 9 February 2021 / Revised: 11 March 2021 / Accepted: 22 March 2021 / Published: 25 March 2021
(This article belongs to the Section Mechanical Engineering)

Abstract

:
In this paper, theoretical derailment equations for cross-wind with frequency were derived to assess running safety. For a KTX (Korean high-speed train) unit, the wheel unloading ratios, which are the criteria for evaluating derailments in UIC (International union of railways) and TSI (Technical Specification for Interoperability) regulations, were calculated through the formula under the driving regulations according to cross-wind speeds, and the theoretical results were compared and evaluated through a multibody dynamics (MBD) simulation. In addition, the wheel unloading ratios were calculated for various frequencies of cross-winds. As a result of the formula and MBD, the wheel unloading ratios were shown to increase rapidly regardless of the dampers in suspension when the cross-wind frequency and the natural frequency of a vehicle were in agreement. Finally, we calculated the changes of wheel unloading ratio for different track gauges and found that these theoretical equations could calculate more accurate results than the existing Kunieda’s formula. The formula derived in this study has the advantage of considering various variables, such as fluctuant cross-winds, rail irregularities, and derailment behaviors, which were not considered in previous studies or Kunieda’s formula. It could be used for setting suspensions or railway vehicle specifications in the initial design stage.

1. Introduction

There are a lot of derailment causes for railway vehicles, such as improper wheel–rail interactions, track irregularities, and cross-winds. Recently, the speeds of railway vehicles including freight trains have increased but the weights of railway vehicles have decreased. For these reasons, cross-winds have increasingly become the cause of derailment accidents [1]. In addition, the strength of typhoons or gusts has increased due to climate change. So the risk of derailments by strong winds has increased [2,3]. Derailment accidents by cross-wind frequently occur all around the world, including Japan [4], China [5], and so forth, as shown in Figure 1.
As derailment accidents by strong winds have become serious, the risks of derailment by winds are evaluated by criteria given in standards like EN (European Standard), JIS (Japanese Industrial Standards) at the design stage of railway vehicles. In Europe, the safety of cross-winds by CWC (characteristic wind curves) in EN standard has usually been evaluated [6,7], and the EN and UIC regulations provide a model for wind forces acting on a vehicle under various conditions such as a wind force of Chinese hat type and wind pressure acting on a vehicle when entering a tunnel, as shown in Figure 2 [8,9].
Since the 1970s in Japan, Masharu Kunieda [10] proposed Kunieda’s formula to evaluate the risk of derailment by cross-wind because derailment accidents had been frequently caused by strong cross-winds.
Since the shape of railway vehicles is a box type which can be severely affected by cross-wind, much research on derailment by cross-wind has been carried out.
The studies of cross-winds acting on railway vehicles can be divided into two fields: (1) aerodynamics, which considers the wind flow and force effects on vehicles, and (2) vehicle dynamics, which evaluates the risk of derailments of trains using the force effects derived from aerodynamics.
  • Aerodynamics
Zhuang Sun measured the wind force acting on a train when it passed by a windbreak breach, and studied the effects according to the gap size of the windbreak breach through simulations [11]. Mulugeta Biadgo Asress carried out simulations on the ICE 2 high-speed train for a range of various yaw angles to compare with wind tunnel test data, and evaluated the changes of wind flows according to the yaw angles through simulation [12]. S. Giappino studied the aerodynamic coefficients through wind tunnel tests on down-scale models, which were on a lightweight train and a high-speed train, and evaluated rollover risks using the CWC based on static equilibrium [13].
Research on aerodynamics has been carried out through tests, using real or down scale models, and CFD (computational fluid dynamics) simulations. The tests using a real vehicle were carried out only in a special case because of cost and test bed problems. In many studies, down-scale tests and CFD simulations were performed in order to obtain wind force acting on a vehicle or aerodynamic characteristics of a vehicle. Whereas down-scale tests yield highly reliable values of the measured forces, CFD simulation reduces the cost of studies on phenomena, such as the pressure field and the wind flow [14].
  • Vehicle dynamics
Takahiro Hosoi [15] used MBD (multibody dynamics) simulation to obtain the derailment coefficient for all wheels when cross-wind was applied in a curved track. In addition, the derailment coefficients according to vehicle speeds and wind speeds were analyzed. C. Baker [16] performed down-scale tests and CFD simulations to obtain the aerodynamic coefficients of a passenger car and a freight car. The wind forces were calculated based on the results and applied to the MBD simulation. The risk of derailment was assessed under various conditions reflecting wind forces and actual track conditions. Masaharu Kunieda studied the rollover derailment that could be caused by cross-wind and/or centrifugal force using physical parameters such as the center of gravity, aerodynamic coefficient, and the center of wind force on a vehicle body [10]. It was confirmed that the theoretical results obtained through the formula and the statical experimental ones obtained through the full-scale real tests were in good agreement [17]. The formula has been used in Japan to theoretically predict a rollover derailment by cross-winds until now [18]. S. W. Nam researched the rollover derailment safety of the KTX based on Kunieda’s formula, considering the wind forces generated by a train running from the opposite direction and the cross-wind on a curved track [19].
Research on derailments is carried out using MBD simulations and theoretical calculations. By the way, it is difficult to obtain accurate results through theoretical calculation because the behavior of a railway vehicle could be obtained through various variables and nonlinear equations. So, in many studies, the running safety has been evaluated using MBD simulation to obtain results similar to real situations. However, a theoretical derailment analysis is necessary to derive good estimations on suspensions at the initial design stage, because every detail of information on a vehicle must be prepared for MBD simulations, which give a better result.
Although a lot of studies on derailments by cross-wind were conducted, the derailment study on cross-winds with frequencies generated by geographical features or wayside structures is insufficient.
A cross-wind with exciting frequency can act on vehicles due to railside structures that are repeatedly placed, as in Figure 3, when running on a bridge. However, the theoretical derailment equations due to cross-wind in the previous studies could only statically evaluate the derailment risk under a constant wind speed condition [10,19], and the fluctuant wind or gust conditions have been evaluated using a multibody dynamics simulation [15,16,20].
In this paper, a theoretical formula on the risk of derailment was derived considering cross-winds with varying wind forces generated by structures displaced at regular intervals on the trackside. Then, the relationship of the cross-wind frequency to the natural frequency of a vehicle was evaluated. The main content in this study is to predict the behavior of a railway vehicle using theoretical calculations, and the aerodynamic coefficients and cross-wind condition were referred to in previous studies [21,22,23] and the regulation [24] on driving conditions of high-speed trains in Korea.
To consider the dynamic effects of cyclic cross-winds and wheelsets running on tracks with irregularities, we combined this formulation with the single wheelset derailment formula derived by Koo and others in previous studies [25,26,27]. So, this derailment formula was improved for running train derailments, which is different from Kunieda’s formula or previous derailment theories for statically standing train derailments.

2. Vehicle Behaviors under a Constant Cross-Wind

2.1. Formula for a Constant Cross-Wind Acting on a Car Body

Since railway vehicles run on rails, the lateral and vertical forces applied to the rails can be used to assess derailment. When a cross-wind acts on a car body, the forces acting on suspensions of the vehicle must be known to evaluate its derailment. These forces are closely related to the data of the vertical and horizontal displacements and the roll angle of the body.
The free body diagram of a railway vehicle is shown in Figure 4 [28].
The free body diagram is divided into a car body, a bogie, and a wheelset when a cross-wind force with an upper angle of φ is applied at a distance of “d” from the center of gravity of the body.
To take into account the equivalent amount to one wheelset, it is assumed that the mass and rotational inertias of the car body are 1/4 and those of the bogie are 1/2 [29].
For the primary suspension, the vertical and lateral stiffness values of rubber-metal springs installed at an axle box were denoted as k1l and k1v, respectively, because each axle box is directly connected with the bogie. For the secondary suspension, the vertical and lateral stiffness values of air springs connecting the bogie and the car body were denoted as k2l and k2v considering their half values to each wheelset.
The static equilibrium equation of a vehicle unit to predict the spring forces transmitted by the cross-wind can be summarized as follows.
  • Static equilibrium equation for the car body:
  F x 2 = P t cos + m 2 ( v 2 r cos φ g sin φ ) F R H 2 F L H 2 = 0
  F y 2 = m 2 ( v 2 r sin φ + g cos φ ) P t sin F R V 2 F L V 2 = 0
  M G 2 = P t ( d cos + l sin ) + ( F R H 2 + F L H 2 ) h 2 + ( F L V 2 F R V 2 ) e 2 = 0
  • Static equilibrium equation for the bogie:
  F x 1 = m 1 ( v 2 r cos φ g sin φ ) + F R H 2 + F L H 2 F R H 1 F L H 1 = 0
  F y 1 = m 1 ( v 2 r sin φ + g cos φ ) + F R V 2 + F L V 2 F R V 1 F L V 1 = 0
  M G 1 = ( F R H 2 + F L H 2 ) h 1 u + ( F R V 2 F L V 2 ) e 2 + ( F R H 1 + F L H 1 ) h 1 l + ( F L V 1 F R V 1 ) e 1 = 0
  • Static equilibrium equation for a wheelset:
  F x 0 = F L H 1 + F R H 1 + m 0 ( v 2 r cos φ g sin φ ) R R H = 0
  F y 0 = m 0 ( v 2 r sin φ + g cos φ ) + F L V 1 + F R V 1 R R V R L V = 0
  M G 0 = R R H h 0 + ( F R V 1 F L V 1 ) e 1 + ( R L V R R V ) e 0 = 0
The forces of the secondary suspension springs generated by the relative displacements between the body and the bogie are as follows.
F R H 2 = F L H 2 = k 2 l ( x 2 x 1 h 2 θ 2 + h 1 u θ 1 )
F R V 2 = k 2 v ( y 2 y 1 + e 2 θ 2 e 2 θ 1 )
F L V 2 = k 2 v ( y 2 y 1 e 2 θ 2 + e 2 θ 1 )
The forces of the primary suspension springs generated by the relative displacements between the bogie and the wheelset are as follows.
F R H 1 = F L H 1 = k 1 l ( x 1 x 0 h 1 l θ 1 )
F R V 1 = k 1 v ( y 1 y 0 + e 1 θ 1 e 1 θ 0 )
F L V 1 = k 1 v ( y 1 y 0 e 1 θ 1 + e 1 θ 0 )
The contact forces between the wheel and the rail can be calculated from the primary suspensions deformations as follows.
R R H = k 0 l ( x 0 h 0 θ 0 )
R R V = k 0 v ( y 0 + e 0 θ 0 )
R L V = k 0 v ( y 0 e 0 θ 0 )

2.2. Dynamic Equations for a Cyclical Cross-Wind Acting on the Body

Winds in nature can generally blow at a constant speed but sometimes in the form of momentary gusts in the EN regulation. However, a cross-wind with various frequency components acts on railway vehicles by means of wayside structures repeatedly positioned or a cross-running train from the opposite direction.
In this study, we assumed that the cyclical cross-wind consists of one sine component. We evaluated how it affects derailment of a railway vehicle.
We can consider a wind force with frequency as the summation of a constant wind force of P t and an alternating wind force of F w sin ω t as shown in Figure 5.
In order to derive the dynamic equations for the cross-wind with frequency, the free body diagram in term of a deformed shape is shown in Figure 6. In this step, the centrifugal force and gravity were not considered because their effects were considered in the previous step in Figure 4.
And then, the final forces of suspensions and the behavior of a train can be calculated by superimposing the two results obtained in the previous step (Figure 4) and this step (Figure 6). Here we assumed the cyclic effects are small compared to those of the constant force.
The equation of motion for the car body:
  F x 2 = F t ( sin ω t ) cos F R H 2 F L H 2 = m 2 x ¨ 2
  F y 2 = F t ( sin ω t ) sin F R V 2 F L V 2 = m 2 y ¨ 2
  M G 2 = F t ( d cos + l sin ω t ) + h 2 ( F R H 2 + F L H 2 ) + e 2 ( F L V 2 F R V 2 ) = I 2 θ ¨ 2
The forces of the secondary suspension springs generated by relative displacements between the body and the bogie can be obtained by the Equations (10)~(12), but the following assumptions are made to simplify these equations.
k 1 l k 2 l ,   k 1 v k 2 v     x 2 x 1 ,   y 2 y 1 ,   θ 2 θ 1
From the Equations (19) and (21), the natural frequencies (ω2a, ω2b) for the lateral and rolling directions of the body are derived as the following Equations (22) and (23).
In general, since most of the cross-wind acting on a vehicle has low frequencies, the first natural frequency of the car body on suspensions becomes an important factor to assess derailment safety.
w 2 a 2 = 1 2 { ( 2 k 2 l m 2 + 2 k 2 l h 2 2 I 2 + 2 k 2 v e 2 2 I 2 ) ( 2 k 2 l m 2 + 2 k 2 l h 2 2 I 2 + 2 k 2 v e 2 2 I 2 ) 4 2 k 2 l m 2 2 k 2 v e 2 2 I 2 }
w 2 b 2 = 1 2 { ( 2 k 2 l m 2 + 2 k 2 l h 2 2 I 2 + 2 k 2 v e 2 2 I 2 ) + ( 2 k 2 l m 2 + 2 k 2 l h 2 2 I 2 + 2 k 2 v e 2 2 I 2 ) 4 2 k 2 l m 2 2 k 2 v e 2 2 I 2 }
From the free body diagram of Figure 6, the equations of motion of the bogie and the wheelset can be derived as follows.
  • The equation of motion for the bogie frame:
  F x 1 = F R H 2 + F L H 2 F R H 1 F L H 1 = m 1 x ¨ 1
  F y 1 = F R V 2 + F L V 2 F R V 1 F L V 1 = m 1 y ¨ 1
  M G 1 = h 1 u ( F R H 2 + F L H 2 ) + h 1 l ( F R H 1 + F L H 1 ) + e 2 ( F R V 2 F L V 2 ) + e 1 ( F L V 1 F R V 1 ) = I 1 θ ¨ 1
  • The equation of motion for the wheelset:
  F x 0 = F R H 1 + F L H 1 R R H = m 0 x ¨ 0
  F y 0 = F R V 1 + F L V 1 R R V R L V = m 0 y ¨ 0
  M G 1 = R R H h 0 + e 1 ( F R V 1 F L V 1 ) + e 0 ( F L V F R V ) = I 0 θ ¨ 0
It is possible to predict the forces of suspensions and wheel/rails, and the behavior of a vehicle by solving the simultaneous differential Equations (19)–(21), and (24)–(29).

3. A Dynamic Simulation Model to Evaluate the Formulation

3.1. A Model for Multibody Dynamics Simulation

The commercial software “Simpack” was used to evaluate the derived equations [30]. Here, Simpack is a general multibody dynamics (MBD) simulation software enabling users to simulate the nonlinear motion of any mechanical or mechatronic system.
The simulation model was based on the power car of the KTX. The compositions of the primary and secondary suspensions are shown in Figure 7. The characteristics of the suspensions are shown in Table 1, Table 2 and Table 3 [21,22]. In addition, the simulation model to evaluate the theoretical equations is shown in Figure 8. Table 4 shows the height of the center of vehicle gravity, and the masses of the body, the bogie, and the wheelset [23].
As mentioned above, the theoretical equations were derived as a 2D formula simplified for a real vehicle. So, the vehicle parameters must be converted for a 2D quarter model on one wheelset. Table 5 shows the modified input data for 2D dynamic simulations matching the 2D quarter model.

3.2. Evaluation of Natural Frequency of the Vehicle

The natural frequency of the power car is a unique characteristic determined by its weight and suspensions, and it is a factor closely related to the derailment safety when considering cross-winds with low frequency.
The natural frequencies for lateral and rolling directions are derived as shown in Equations (22) and (23) using the motion equations of the vehicle in Section 2.2.
The natural frequencies of the derived formula and the MBD model were compared in Table 6. Figure 9 show the behaviors of the vehicle in which the exciting frequency component coincides with its natural frequency.

4. Evaluation of Derailment Safety under Cross-Wind

4.1. Operating Regulations of the KTX According to Wind Speed

The operating regulations of the KTX for cross-wind speeds are enforced, as shown in Table 7, to prevent derailments [24].
Table 8 shows the aerodynamic coefficient and side areas of the KTX, and air density. The values related to the aerodynamics of the KTX are based on the previous research and data from vehicle manufacturers [19]. Equation (30), which is based on Bernoulli’s Equation, is used to calculate aerodynamic forces due to cross-winds [31].
The regulations of Table 7 for wind speeds are converted into the equivalent lateral forces by using Equation (30) and the aerodynamic characteristics of the KTX in Table 8.
F w i n d = C y × A × ρ v 2 2

4.2. Criteria of Derailment

In advanced railway countries, such as Europe, North America, and Japan, there are various technical criteria to prevent derailments, as in Table 9 [32,33,34,35].
Although there are many ways to judge derailments, the ratio of the lateral and the vertical forces (P/Q) is generally used for a wheel climb-up derailment. The wheel unloading ratio ( Q/Q0) is used to assess a rollover derailment. As cross-winds usually cause a rollover derailment, the rollover risk is assessed using the wheel unloading ratio.
The definitions of P and Q are as follows: P is for the lateral force, Q for the vertical force.

4.3. Derailment Safety According to the Cross-Wind Regulations

The ratios of wheel unloading between simulations and the theoretical formula were compared under the driving conditions (Table 7), in the case that wind forces are applied in a horizontal direction ( = 0 ° in Figure 4 and Figure 6) on the ground.
Table 10 shows the ratio of wheel unloading when constant cross-winds act on the KTX power car on the straight track and curved tracks (lateral acc. 0.5 [m/s2]; this is the case of 1.0 [m/s2] referred to in the EN regulation [EN 14067]). Here it is assumed that there are no lateral forces transmitted to the wheelset due to rail irregularities.
Figure 10 shows the contact forces of wheel/rail and the ratio of wheel unloading. When the wind speed of about 65 [m/s] acts on the vehicle, the wheel unloading ratio exceeds the TSI criteria for the rollover derailment ( Q / Q < 0.9 ).
On the other hand, basic lateral forces acting on running wheelsets take place due to wheel conicity and rail irregularities, even though there is no cross-wind. In the previous study, the KTX has a measured coefficient of derailment (P/Q = 0.315) during normal running on a high-speed rail [36,37]. Using the measured coefficient of derailment and Equation (31), the P/Q can be converted to an additional lateral force acting on the wheelset, as given in Equation (32).
The ratios of wheel unloading according to cross-wind speeds are compared in Table 11. These results were obtained considering both the wheel–rail interactions and the lateral forces due to cross-winds on the running vehicle.
P Q = ( F L H 1 + F R H 1 ) ( 1 + μ β ) μ ( 1 + γ ) F L V 1 + μ γ F R V 1 ( 1 + γ ) F R V 1 γ F L V 1 + β ( F L H 1 + F R H 1 )
F H = [ P Q ( 1 + γ ) μ γ ] F L V 1 + [ μ ( 1 + γ ) P Q γ ] F R V 1 [ ( 1 + μ β ) P Q β ]
Here, β = h 0 e 1 , γ = 2 e 0 e 1 (see Figure 4).
h0: Radius of wheel
e0: Track width (distance between two contact points of both treads)
e1: Distance between the contact point of a wheel tread and the primary suspension
μ: Friction coefficient

4.4. Derailment Safety Due to the Cross-Wind Frequencies

When a constant wind force acts on a vehicle, the ratio of wheel unloading must be within 0.9 not to be derailed on the straight track and the curved tracks (lateral acc. 0.5 [m/s2], 1.0 [m/s2] in the EN). However, wind speeds in nature are not constant. In addition, the winds generated by wayside structures on the track or other opposite trains make exciting frequencies [38].
In this section, the ratios of wheel unloading according to the winds with frequency were compared using the theoretical formula derived in Section 2 and the dynamic simulations.
The conditions to compare the ratios of wheel unloading according to different wind frequencies are shown in Table 12.

4.4.1. The 2D Dynamic Vehicle Model without Dampers

In the 2D dynamic vehicle model of Section 2, the dampers in suspension are not considered to formulate theoretical equations. Therefore, if the lateral natural frequency of the vehicle and the frequency of the cross-wind are matched, resonance phenomena will occur and the amplitude will be continuously increased as time goes by.
The resonance phenomena can easily be described using the dynamic equation of a simple spring-mass system. When a periodic external force is applied to the system without dampers, the dynamic equation is shown in Equation (33), and its general solution is shown in Equation (34) [39].
m x ¨ + k x = F 0 cos ω t
x ( t ) = A cos ω 0 t + B sin ω 0 t + F 0 m ( ω 0 2 ω 2 ) cos ω t
A and B are constants determined by the initial conditions and ω 0 is a natural frequency.
ω 0 = k / m
If the frequency of external forces ( ω ) is equal to or very near to that of the vehicle ( ω 0 ), the solution of the dynamic equation can be converted as shown in Equation (36).
x ( t ) = A cos ω 0 t + B sin ω 0 t + F 0 2 m ω 0 2 t sin ω 0 t
Because of the term t sin ω 0 t , the solution continuously increases with the passage of time. This phenomenon is called a resonance as shown in Figure 11.
However, the actual springs in a vehicle suspension must not be extended longer than the predesigned distance not to be permanently deformed or damaged over the design limit. So, the spring constant (k) must be assumed as a nonlinear over the linear elastic range of a normal spring, which means Equation (33) is invalid in an excessive amplitude condition.
Figure 12 and Table 13 show the results that the ratios of the wheel unloading of the simulation and the current formula without dampers are rapidly increased at a resonance condition.
However, the results using Kunieda’s formula are constant regardless of the cross-wind frequency because the formula cannot consider it. Also, the formula cannot take into account the change of ratio of wheel unloading by rail irregularities. Thus, when rail irregularity (P/Q = 0.315) is applied, as shown in Figure 12, the value of the ratio of wheel unloading is low compared to other results, even in low frequency.
Even though the ratio of the wheel unloading can have 1 as the maximum value, the TSI criteria regards over 0.9 as a derailment. Here, the simulation contact forces were filtered using 1.5 [Hz] low-pass filter to remove high-frequency components caused by numerical factors.
Although oscillatory lateral forces between wheel–rails are generated by the wheel conicity in the simulation, our formulations could not consider it. These oscillation effects gave higher values in simulation compared to the results of our formula in Table 13.

4.4.2. The 2D Dynamic Vehicle with Dampers

The amplitude (X) and the phase angle ( δ ) are obtained using the dynamic equation with a spring damper and an external force with a constant frequency. The dimensionless amplitude (X) and the phase angle ( δ ) are expressed as shown in Equations (37) and (38).
X k F 0 = 1 [ 1 ( ω ω 0 ) 2 ] 2 + [ 2 ζ ( ω ω 0 ) ] 2
tan δ = 2 ζ ( ω ω 0 ) 1 ( ω ω 0 ) 2
Here, the used symbols are defined in Table 14.
The above equations show that the amplitude and phases are the functions of frequency ratio ( ω / ω 0 ) and damping ratio ( ζ ), and they can be shown in Figure 13.
In general, an actual train has dampers, as shown in Figure 14, which is expected to reduce the amplitude at resonant frequencies as described above. The characteristics of dampers installed in the KTX were applied for accurate results because the dampers can affect the simulation results. Figure 15 shows the simulation results between the models with dampers and without dampers.
Although the ratios of the wheel unloading showed a little difference in frequencies other than resonance, the results showed a good agreement and exceeded the TSI criteria of 0.9 in the frequency near to the resonance regardless of dampers.
Even though the characteristics of dampers can affirmatively affect the running safety, the phenomenon of a resonance is critical to derailments regardless of dampers. In a design stage, it is important to evaluate how much a cross-wind with frequency near to the natural frequency of a vehicle can affect the running safety regardless of dampers.

4.5. Derailment Safety According to Track Gauges

Track gauges are generally classified as narrow gauges, standard gauges, and wide gauges. Various track gauges are used around the world.
High-speed trains operating above 300 [km/h] generally use the standard gauge of 1435 [mm]. However, in order to compare the effects of the track gauges, simulations and theoretical formulas were calculated according to the track gauges. The risks of derailment were assessed considering the different gauges under the driving speed and the cross-wind condition in Table 15.
Figure 16a shows the change of the ratios of wheel unloading when a train runs on a straight track according to the gauges. The results of our formula, Kunieda’s formula, and MBD simulations showed that the larger the gauges, the less the risks.
The deviations between our formula and the MBD simulation were only about 0.04 on the straight track. It is shown that the theoretical equations of this paper were more accurate to predict derailment risks than Kunieda’s formula when considering variables such as rail irregularity which are not considered in existing theory.
Figure 16b shows the change of the ratios of wheel unloading when a train runs on the curved track according to the gauges. The results show the decreasing tendency like a form of a quadratic curve. Although the overall trend of the results using Kunieda’s formula and the theory of this paper was similar, the wheel unloading ratios using Kunieda’s formula were relatively large. In addition, the results of simulations and the theoretical ones of this research showed very similar values, but our theoretical solutions could give more accurate results than those of Kunieda’s formula.

5. Conclusions

Even though there was a lot of theoretical research to predict the risk of derailment due to cross-winds, most of them could not consider the variables like wheel–rail interactions and frequencies of cross-winds. So, there occurred some differences compared with actual results.
In this study, a new theoretical formula was proposed to predict the risk of derailment due to cross-winds with frequency considering wheel–rail interactions. The proposed formula was verified using MBD simulations.
The following conclusions were obtained:
(1)
The static equilibrium equations for the body, the bogie, and the wheelset were derived under a constant force of winds applied to the vehicle using the 2D free body diagram. The wheel–rail and suspension forces, and the displacements and rotation angle of the vehicle were obtained by solving the equations for each static variable.
(2)
The dynamic equations for the body, the bogie, and the wheelset were derived under a wind force with frequency components applied to the vehicle using the 2D free body diagram. The wheel–rail and suspension forces, and the displacements and rotation angle of the vehicle were obtained by solving the equations for each dynamic variable. In addition, the lateral and rolling natural frequencies of the vehicle were derived by solving the motion equations for the lateral and rolling directions.
(3)
The wheel–rail forces obtained from summing the dynamic and static ones were improved using a known derailment coefficient to consider wheel–rail interactions during a normal train running without winds.
(4)
The multibody dynamics software (Simpack) model was generated using the specifications of a Korean high-speed train (KTX) to evaluate our theoretical formula. For the operational regulations of the KTX under cross-wind conditions, the ratios of wheel unloading on the dynamic simulation and our theoretical formula were compared and well matched.
(5)
The ratio changes of wheel unloading were investigated under the presence or absence of suspension dampers considering the natural frequency of a vehicle and the cross-wind frequency. As a result, the ratios were found to exceed the TSI criteria regardless of the presence or absence of dampers if the natural frequency of the vehicle and the cross-wind frequency were matched.
(6)
The ratio changes of wheel unloading were investigated according to track gauges using the simulations, Kunieda’s formula, and our theoretical one. Although all the results were similar on the straight and curved tracks, the wheel unloading ratios using Kunieda’s formula were larger than the simulations and ours which are in a good agreement. So, our theoretical formula gives more accurate results than Kunieda’s formula, because it considers various operational conditions like a cross-wind frequency and a normal running derailment coefficient due to wheel–rail interaction.
Therefore, this formula has strong advantages over the previous studies or the commonly used Nadal or Kunieda’s formula because it considers various variables, such as fluctuant cross-winds, rail irregularities, and dynamic derailment behaviors. It would be possible to predict well the derailment risk under actual operating conditions. It could be used for setting suspensions or railway vehicle specifications in the initial design stage. However, this formula has the disadvantage of being more complicated when compared to the Nadal or Kunieda’s formula.

Author Contributions

Conceptualization, J.K. (JunWoo Kim), I.S. and J.K. (JeongSeo Koo); methodology J.K. (JunWoo Kim), I.S. and J.K. (JeongSeo Koo); software J.K. (JunWoo Kim) and I.S.; validation J.K. (JunWoo Kim); writing—review and editing, J.K. (JunWoo Kim) and I.S.; supervision, J.K. (JunWoo Kim); project administration, J.K. (JeongSeo Koo); funding acquisition, J.K. (JeongSeo Koo). All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Research Program funded by SeoulTech (Seoul National University of Science and Technology).

Data Availability Statement

The data of the rolling stock used to support findings of this study have been stated in the references [19,21,22,23].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Deng, Y.; Xiao, X. Effect of Cross-Wind on High-Speed Vehicle Dynamic Derailment. Logistics 2009, 5, 2287–2293. [Google Scholar]
  2. Gerd, M.; Peter, D.; Burkhard, S.-W. Effects of Strong Cross Winds on High-Speed Trains: A Methodology for Risk Assessment and Development of Countermeasures; Deutsche Bahn AG Research & Technology Center: Berlin, Germany, 2001; pp. 1–15. [Google Scholar]
  3. Bae, H.U.; Min, K.H.; Go, Y.S.; Cheon, H.S.; Lim, N.H. Collision Behavior Analysis between Derailed Train and Structure According to Collision Analysis Condition. J. Korean Soc. Hazard Mitig. 2017, 17, 45–54. [Google Scholar] [CrossRef]
  4. Tamura, Y. Wind-induced damage to buildings and disaster risk reduction. In Proceedings of the APCWE-VII, Taipei, Taiwan, 8–12 November 2009. [Google Scholar]
  5. Zhai, W.; Yang, J.; Li, Z.; Han, H. Dynamics of high-speed train in crosswinds based on an air-train- track interaction model. Wind Struct. 2015, 20, 143–168. [Google Scholar] [CrossRef]
  6. DB Netz AG. Ausgewählte Maßnahmen und Anforderungen an das Gesamtsystem Fahrweg/Fahrzeug-Aerodynamik/Seitenwind; DB Netz AG: Frankfurt, Germany, 2006. [Google Scholar]
  7. Baker, C.; Cheli, F.; Orellano, A. Nicolas Paradot, Carsten Proppe and Daniele Rocchi, “Cross-wind effects on road and rail vehicles”. Veh. Syst. Dyn. 2009, 47, 983–1022. [Google Scholar] [CrossRef]
  8. European Committee for Standardization. EN 14067-5, Railway Applications—Aerodynamics—Part 5: Requirements and Test Procedures for Aerodynamics in Tunnels; European Committee for Standardization: Brussels, Belgium, 2010. [Google Scholar]
  9. European Committee for Standardization. EN 14067-6, Railway Applications—Aerodynamics—Part 6: Requirements and Test Procedures of Cross Wind Assessment; European Committee for Standardization: Brussels, Belgium, 2016. [Google Scholar]
  10. Kunieda, M. Theoretical Study on the Mechanics of Overturn of Rail Rolling Stock. Rtri Railw. Tech. Res. Rep. 1972, 793, 173–186. [Google Scholar]
  11. Sun, Z.; Dai, H.; Hemida, H.; Li, T.; Huang, C. Safety of High-Speed Train Passing by Windbreak Breach with Different Sizes. Veh. Syst. Dyn. 2020, 58, 1935–1952. [Google Scholar] [CrossRef]
  12. Biadgo Asress, M.; Svorcan, J. Numerical investigation on the aerodynamic characteristics of high-speed train under turbulent crosswind. J. Mod. Transp. 2014, 22, 225–234. [Google Scholar] [CrossRef] [Green Version]
  13. Giappino, S.; Rocchi, D.; Schito, P.; Tomasini, G. Cross wind and rollover risk on lightweight railway vehicles. J. Wind Eng. Ind. Aerodyn. 2016, 153, 106–112. [Google Scholar] [CrossRef]
  14. Yu, M.; Jiang, R.; Zhang, Q.; Zhang, J. Crosswind stability evaluation of high-speed train using different wind models. Chin. J. Mech. Eng. 2019, 32, 1–13. [Google Scholar] [CrossRef] [Green Version]
  15. Hosoi, T.; Tanifuji, K. Effect of Crosswind on Derailment of Railway Vehicles Running on Curved Track at Low Speed. Int. J. Railw. 2012, 5, 93–101. [Google Scholar] [CrossRef] [Green Version]
  16. Baker, C.; Hemida, H.; Iwnicki, S.; Xie, G.; Ongaro, D. Integration of Crosswind Forces into Train Dynamic Modelling. Proc. Inst. Mech. Eng. Part F J. Rail Rapid Transit 2011, 225, 154–164. [Google Scholar] [CrossRef] [Green Version]
  17. Yu, H.; Shimomura, T.; Tanifuji, K. Full-scale experiment on the behavior of a railway vehicle being subjected to lateral force. J. Mech. Syst. Transp. Logist. 2010, 3, 35–43. [Google Scholar]
  18. Yu, H.; Misu, Y.; Kurihara, T. Atsushi Moriyama and Makoto Shimamura, Study of New Methods for Train Operation Control in Strong. Wind. Jr East Tech. Rev. 2011, 19, 31–36. [Google Scholar]
  19. Nam, S.W. Theoretical Analysis on Overturn Safety of Train affected by Wind Pressure. J. Korean Soc. Railw. 2012, 15, 537–542. [Google Scholar] [CrossRef]
  20. Xiao, X.B.; Ling, L.; Xiong, J.Y.; Zhou, L.; Jin, X.S. Study on the safety of operating high-speed railway vehicles subjected to crosswinds. J. Zhejiang Univ. Sci. A 2014, 15, 694–710. [Google Scholar] [CrossRef] [Green Version]
  21. Hong, T.D. A Study on the Improvement of KTX Bogie Instability. Master’s Thesis, University of Seoul, Seoul, Korea, 2006; pp. 1–88. [Google Scholar]
  22. Eom, B.G.; Han, B.Y.; Kang, B.B.; Lee, H.S. A Study on 1/5 Scaled Bogie Design of KTX-Sancheon using Similarity Laws. In Proceedings of the Korean Society for Railway Spring Conference, Korea, Mokpo, 17–19 May 2012; pp. 1772–1777. [Google Scholar]
  23. Kim, M.S. A Study on Derailment Prediction of Railway Vehicle Due to Cross-Wind Using the Derailment Theory of Single Wheelset. Master’s Thesis, Seoul National University of Science & Technology, Seoul, Korea, 2014; pp. 1–59. [Google Scholar]
  24. Choi, J.H.; Park, T.W.; Sim, K.S.; Kwak, M.H.; Lee, D.H. Stability Evaluation on Aerodynamics of High Speed Railway Train. Trans. Korean Soc. Noise Vib. Eng. 2012, 22, 244–252. [Google Scholar] [CrossRef]
  25. Koo, J.S.; Choi, S.Y. Theoretical development of a simplified wheelset model to evaluate collision-induced derailments of rolling stock. J. Sound Vib. 2012, 331, 3172–3198. [Google Scholar] [CrossRef]
  26. Oh, H.S. A Study on Mechanical Parameters Influencing Derailment of a Wheelset. Master’s Thesis, Seoul National University of Science & Technology, Seoul, Korea, 2013; pp. 1–51. [Google Scholar]
  27. Kim, M.S.; Kim, G.Y.; Kim, H.T.; Koo, J.S. Theoretical cross-wind speed against rail vehicle derailment considering the cross-running wind of trains and the dynamic wheel-rail effects. J. Mech. Sci. Technol. 2016, 30, 3487–3498. [Google Scholar] [CrossRef]
  28. Kim, G.Y.; Koo, J.S.; Rho, J.H. Theoretical Evaluation of Rollingstock Derailment Safety under Cyclical Crosswind Condition. J. Mech. Sci. Technol. 2018, 42, 95–101. [Google Scholar]
  29. Choi, S.Y. Development of a Wheelset Derailment Model for Prediction of Derailment Behavior after Train Collision. Master’s Thesis, Seoul National University of Science & Technology, Seoul, Korea, 2011; pp. 1–54. [Google Scholar]
  30. Dassault System Simulia Corp. Simpack Documentation; Dassault System Simulia Corp.: Velizy-Villacoublay, France, 2017. [Google Scholar]
  31. Kwon, H.B.; You, W.H.; Cho, T.H. Measurement of Aerodynamic Loads on Railway Vehicles under Crosswind. Trans. Korean Soc. Mech. Eng. A 2011, 35, 91–98. [Google Scholar] [CrossRef]
  32. Yi, S. The Effect Law of the Curve Parameters of a High-Speed Railway on Vehicle-Line Dynamic Performance. Dyn. Anal. High-Speed Railw. Alignment 2018, 73–108. [Google Scholar] [CrossRef]
  33. Wang, P. Wheel-Rail Relation Design. Des. High-Speed Railw. Turnouts 2015, 97–162. [Google Scholar] [CrossRef]
  34. Soper, D. The Aerodynamics of a Container Freight Train; Springer: Berlin/Heidelberg, Germany, 2016; pp. 10–32. [Google Scholar]
  35. Official Journal of the European Union. COMMISSION DECISION of 21 February 2008 Concerning a Technical Specification for Interoperability Relating to the ‘Rolling Stock’ Sub-System of the Trans-European High-Speed Rail System. Available online: https://eur-lex.europa.eu (accessed on 26 March 2008).
  36. Kim, J.C.; Ham, Y.S. Analysis on Running Safety for KTX Vehicle. Korean Soc. Railw. 2007, 10, 473–479. [Google Scholar]
  37. Ham, Y.S. The Evaluation of Running Safety in the Double Connection Operation for KTX-Sancheon Train. In Proceedings of the the Korean Society of Mechanical Engineers Autumn Conference, Jeju, Korea, 28–30 October 2010; pp. 155–156. [Google Scholar]
  38. Sun, Z.; Dai, H.; Gao, H.; Li, T.; Song, C. Dynamic performance of high-speed train passing windbreak breach under unsteady crosswind. Int. J. Veh. Mech. Mobil. 2018, 408–424. [Google Scholar] [CrossRef]
  39. Inman, D.J. Engineering Vibrations; PEARSON: London, UK, 2012; pp. 109–121. [Google Scholar]
Figure 1. Derailment accidents in Japan by strong cross-winds. [4].
Figure 1. Derailment accidents in Japan by strong cross-winds. [4].
Applsci 11 02947 g001
Figure 2. Tunnel pressure profile at a position just behind the nose of a train.
Figure 2. Tunnel pressure profile at a position just behind the nose of a train.
Applsci 11 02947 g002
Figure 3. Structures installed with regular intervals.
Figure 3. Structures installed with regular intervals.
Applsci 11 02947 g003
Figure 4. The free body diagram of a vehicle under a cross-wind.
Figure 4. The free body diagram of a vehicle under a cross-wind.
Applsci 11 02947 g004
Figure 5. The wind force with a frequency component.
Figure 5. The wind force with a frequency component.
Applsci 11 02947 g005
Figure 6. The free body diagram under the cyclic wind force with a specific frequency.
Figure 6. The free body diagram under the cyclic wind force with a specific frequency.
Applsci 11 02947 g006
Figure 7. The compositions of suspensions.
Figure 7. The compositions of suspensions.
Applsci 11 02947 g007
Figure 8. The 3D simulation model (Simpack).
Figure 8. The 3D simulation model (Simpack).
Applsci 11 02947 g008
Figure 9. The vehicle behavior under the resonant frequencies.
Figure 9. The vehicle behavior under the resonant frequencies.
Applsci 11 02947 g009
Figure 10. The ratios of wheel unloading to cross-wind speeds under the overturning derailment condition.
Figure 10. The ratios of wheel unloading to cross-wind speeds under the overturning derailment condition.
Applsci 11 02947 g010
Figure 11. The amplitude under a resonance condition.
Figure 11. The amplitude under a resonance condition.
Applsci 11 02947 g011
Figure 12. The comparison of the ratios of wheel unloading according to the frequencies of a cross-wind.
Figure 12. The comparison of the ratios of wheel unloading according to the frequencies of a cross-wind.
Applsci 11 02947 g012
Figure 13. Nondimensional values of amplitude depending on damping factors.
Figure 13. Nondimensional values of amplitude depending on damping factors.
Applsci 11 02947 g013
Figure 14. The motor bogie of the KTX.
Figure 14. The motor bogie of the KTX.
Applsci 11 02947 g014
Figure 15. The comparison of the ratios of wheel unloading between the models with dampers and without dampers (simulation results).
Figure 15. The comparison of the ratios of wheel unloading between the models with dampers and without dampers (simulation results).
Applsci 11 02947 g015
Figure 16. The ratios of wheel unloading according to track gauges.
Figure 16. The ratios of wheel unloading according to track gauges.
Applsci 11 02947 g016
Table 1. Specifications of the primary suspensions.
Table 1. Specifications of the primary suspensions.
SuspensionDouble Coil SpringGuide SpringVertical Oil Damper
specificationsX direction0.581320.0-
Y direction0.58134.5-
Z direction0.750.250.01
T direction0.10.044-
P direction0.10.03-
W direction0.10.473-
UnitsMN/m
MN·m/rad
MN/m
MN·m/rad
MN·s/m
Table 2. Specifications of an anti-yaw damper.
Table 2. Specifications of an anti-yaw damper.
Specification
Velocity [m/s]0.00.00260.010.020.3
Force [kN]0.011.011.612.019.0
Table 3. Specifications of the secondary suspensions.
Table 3. Specifications of the secondary suspensions.
SuspensionCoil SpringVertical Oil
Damper
Lateral Oil DamperPivot Joint
SpecificationsX direction0.15--10.29
Y direction0.15-0.1-
Z direction0.6340.02--
T direction0.1---
P direction0.1---
W direction0.1---
UnitsMN/m
MN·m/rad
MN/m
MN·m/rad
MN·s/mMN/m
Table 4. Specifications of the power car.
Table 4. Specifications of the power car.
Values
WheelsetWheel weight1745[kg]
Moment of inertia740[kg·m2]
Height of weight center0.46[m]
Bogie
frame
Bogie weight2420[kg]
Moment of inertia972.5[kg·m2]
Height of weight center0.56[m]
Car bodyCar body weight54,960[kg]
Moment of inertia13,246.36[kg·m2]
Height of weight center1.72[m]
Height of the wind force center on the car body2.015[m]
Table 5. Input data for 2D simulations.
Table 5. Input data for 2D simulations.
ParametersValuesParametersValues
m 0 1745[kg] h 0 0.46[m]
I 0 740[kg·m2] e 0 0.7175[m]
m 1 1210[kg] h 1 u 0.067[m]
I 1 972.5[kg·m2] h 1 l 0.435[m]
e 1 1[m]
m 2 13,740[kg] d 0.295[m]
I 2 13,246[kg·m2] l 1.335[m]
e 2 1.025[m] h 2 1.093[m]
k 0 l 100[kg/mm] k 0 v 100[kg/mm]
k 1 l 9.58[kg/mm] k 1 v 1.25[kg/mm]
k 2 l 0.51[kg/mm] k 2 v 0.74[kg/mm]
Table 6. The comparison of the resonant frequencies between the theoretical formula and Simpack.
Table 6. The comparison of the resonant frequencies between the theoretical formula and Simpack.
Resonant Frequencies [Hz]
Lateral DirectionRolling Direction
Theory0.942.51
Simpack0.902.77
Table 7. The regulations on driving conditions according to wind speeds.
Table 7. The regulations on driving conditions according to wind speeds.
Wind Speed [m/s]Driving Conditions
45~Operation stop
40~45Under 90 [km/h]
30~40Under 170 [km/h]
~30Decelerated operation
Table 8. The specifications of the high-speed train.
Table 8. The specifications of the high-speed train.
ParametersValues
Aerodynamic side force coefficient( C y )0.6
Side area( A )70.5[m2]
Air density( ρ )1.275[kg/m3]
Table 9. The various evaluation criteria of derailment.
Table 9. The various evaluation criteria of derailment.
Country/
Regulation
Method/CriteriaNote
EU ( P / Q ) s m u < 0.8 2 m moving average of derailment coefficient
UIC P / Q < 0.8 Derailment coefficient
TSI Q / Q 0 < 0.9 Ratio of wheel unloading
North
America
P / Q < 1.0 Derailment coefficient
  | P / Q | < 1.0 ~ 1.5 Weinstock coefficient
Q / Q 0 < 0.9 Ratio of wheel unloading
RTRI
(Japan)
Q / Q 0 < 1.0 Ratio of wheel unloading
Table 10. The ratio of wheel unloading ( Q / Q ) between the theoretical formula and the simulation for cross-winds.
Table 10. The ratio of wheel unloading ( Q / Q ) between the theoretical formula and the simulation for cross-winds.
Track ConditionWind Speeds [m/s]Ratios of Wheel UnloadingDifferences
TheorySimulation
Straight track0000
300.1200.1190.001
400.2330.2200.013
450.2920.2760.016
Curved track
(Lateral Acc. = 0.5 m/s2)
00.1280.1200.008
300.2480.2480
400.3430.3530.010
450.3980.3950.003
Curved track
(Lateral Acc. = 1.0 m/s2)
00.2560.2420.014
300.3760.3640.012
400.4710.4630.008
450.5260.5260
Table 11. The ratios of wheel unloading ( Q / Q ) between the theoretical formula and the simulation for cross-wind along with dynamic wheel–rail interactions (P/Q = 0.315).
Table 11. The ratios of wheel unloading ( Q / Q ) between the theoretical formula and the simulation for cross-wind along with dynamic wheel–rail interactions (P/Q = 0.315).
Track ConditionWind Speeds
[m/s]
Ratios of
Wheel Unloading
Differences
TheorySimulation
Straight track00.1980.1730.025
300.3180.3040.014
400.4130.4030.010
450.4680.4620.006
Curved track
(Lateral Acc. = 0.5 m/s2)
00.3230.3050.018
300.4430.4290.014
400.5380.5280.010
450.5930.5870.006
Curved track
(Lateral Acc. = 1.0 m/s2)
00.4480.4540.006
300.5680.5790.011
400.6630.6810.018
450.7180.7390.021
Table 12. The conditions to compare the derailment coefficients according to cross-wind frequencies.
Table 12. The conditions to compare the derailment coefficients according to cross-wind frequencies.
ParametersValues
TrackStraight
Rail irregularity (P/Q)0.315
Vehicle velocity300[km/h]
WindAverage speed45[m/s]
Fluctuation speed ± 4.5[m/s]
Table 13. The comparison between the theoretical formulas (current study and Kunieda’s formula) and simulations.
Table 13. The comparison between the theoretical formulas (current study and Kunieda’s formula) and simulations.
Frequency
[Hz]
TheorySimulation
(P/Q = 0.315)
Current Study (P/Q = 0.315)Kunieda’s Formula
(P/Q = 0)
0.10.4880.3550.528
0.30.4960.3550.530
0.50.5100.3550.532
0.70.5480.3550.579
0.90.9080.3551.000
1.10.5850.3550.636
1.30.5250.3550.580
1.50.5070.3550.569
Table 14. The symbols used in the damped vibration equation.
Table 14. The symbols used in the damped vibration equation.
Natural frequency ω 0 = k m
Critical damping c c = 2 m k = 2 m ω 0
Damping factor ζ = c c c
c ω k = c c c   c c ω k = 2 ζ ω ω 0
Table 15. The symbols used in the damped vibration equation.
Table 15. The symbols used in the damped vibration equation.
ParametersValues
WindSpeed30[m/s]
Angle0[°]
Vehicle speed300[km/h]
Rail irregularity (P/Q)0.315
TrackStraightCurved
Lateral acc.-1[m/s2]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, J.; Song, I.; Koo, J. Prediction of Theoretical Derailments Caused by Cross-Winds with Frequency. Appl. Sci. 2021, 11, 2947. https://doi.org/10.3390/app11072947

AMA Style

Kim J, Song I, Koo J. Prediction of Theoretical Derailments Caused by Cross-Winds with Frequency. Applied Sciences. 2021; 11(7):2947. https://doi.org/10.3390/app11072947

Chicago/Turabian Style

Kim, JunWoo, InHo Song, and JeongSeo Koo. 2021. "Prediction of Theoretical Derailments Caused by Cross-Winds with Frequency" Applied Sciences 11, no. 7: 2947. https://doi.org/10.3390/app11072947

APA Style

Kim, J., Song, I., & Koo, J. (2021). Prediction of Theoretical Derailments Caused by Cross-Winds with Frequency. Applied Sciences, 11(7), 2947. https://doi.org/10.3390/app11072947

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop