# Theoretical 3D Model for Quasistatic Critical Derailment Coefficient of Railway Vehicles and a Simplified Formula

Authors: Ping Wang; Jian Wang; Xiaochuan Ma; Daolin Ma; Jingmang Xu; Yao Qian
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

The formula for the critical derailment coefficient concerning wheelset yaw angles and wheel-rail creep forces is deduced based on the three-dimensional (3D) force equilibrium relationship in the critical wheel derailment state under quasistatic assumption. The change of critical derailment coefficient and wheel-rail contact patch normal force/creep force as wheelset yaw angles change under the influence of the friction coefficient, maximum flange angle, and net wheel weight is analyzed according to the Kalker linear creep theory and Shen-Hedrick-Elkins creep theory. Analysis shows that the wheel-rail friction coefficient and maximum wheel flange contact angle can significantly influence the critical wheel derailment coefficient, further proving the conservative results when the critical Nadal derailment coefficient is adopted in analyzing wheel derailment under small wheelset yaw angles. To realize easy calculation and application of critical 3D derailment coefficients, the ratio of lateral creep force to longitudinal creep force of wheel-rail contact patches under critical quasistatic wheel derailment conditions is deduced. A simplified calculation method of critical derailment coefficients is presented based on this. The calculation accuracy is verified, proving that it can satisfy engineering application.


---

## Body

## 1. Introduction

Safety is the most important issue of rail transportation, and preventing derailment is a fundamental requirement for ensuring safe rail transportation. The openness of constraints between wheels and rails determines the objective existence of vehicle derailment. Derailment has been a problem throughout the development of rail transportation. However, considering the complexity of the wheel-rail system, no full explanation has been presented for the mechanism of vehicle derailment. As rail transportation becomes faster and heavier loads are carried, derailments can produce disastrous results; further understanding of derailment has become more and more important.As far back as 1896, Nadal, a French scientist [1] derived the force equilibrium conditions of the contact points when wheels show a tendency to climb on the rails under quasistatic assumption according to Coulomb’s friction law and introduced the famous Nadal derailment evaluation criterion. However, in 1976, Gilchrist and Brickle [2] reviewed the Nadal derailment evaluation criterion according to the Hertz contact theory and Duvorol creep force table, pointing out that it, to a certain extent, was conservative under small wheel-rail attack angles and agreed that the Nadal derailment evaluation criterion was appropriate in evaluating derailment under large wheel-rail attack angles. Sweet and Karmel [3–5] simulated wheel derailment under quasistatic assumption according to the two-degree-of-freedom (2-DOF) and three-degree-of-freedom (3-DOF) dynamic wheelset derailment models and obtained consistent results compared with those of the 1 : 5 single wheelset derailment model test. Weinstock [6] suggested that rail derailment was related to not only the derailment coefficient at the flange contact side but also the derailment coefficient at the nonflange contact. Therefore, he pointed out that the sum of the derailment coefficients at both sides of the wheelset could be adopted as a criterion in evaluating wheel derailment, and, to a certain extent, corrected the conservative Nadal derailment evaluation criterion under small attack angles. Elkins and Shust [7, 8] investigated the influence of friction coefficient and wheel-rail attack angle on rail climbing of wheels. They argued that rail climbing of wheels depended on the vehicle running distance when derailment coefficient is out of safe region rather than the duration of the unsafe derailment coefficient. Barbosa [9, 10] comprehensively considered the wheel-rail creep action based on a 3D wheel derailment model established under quasistatic assumption and accordingly improved the Nadal derailment evaluation criterion that was conservative under small wheel-rail attack angles. Braghin et al. [11] studied the derailment of single wheelsets through a rolling rig, and, based on theoretical research, pointed out that the derailment evaluation criterion requiring longitudinal creep forces and attack angles was a new concept for measuring derailment coefficients. Zeng and Wu [12] analyzed the 3D force equilibrium conditions of wheelsets during derailment and presented the derailment evaluation criterion concerning both wheelset derailment coefficient and wheel load reduction rate. Kuo and Lin [13], based on the Weinstock derailment criterion as well as the force equilibrium conditions of wheelset surfaces, introduced a new derailment evaluation criterion and corrected those deficiencies of the Weinstock derailment criterion.Although many studies have been presented, the problem of vehicle derailment has still not been solved and different derailment evaluation criteria still do not satisfactorily coincide with each other. To further understand mechanism of derailment, we carried out this research. In this paper, we establish the force analyzing framework of 3-dimensional critical wheel quasistatic derailment in Section2. Then the formula of the critical derailment coefficient is deduced, and the change of critical derailment coefficient and normal force/creep force of wheel-rail contact patches with wheelset yaw angles under the influence of the friction coefficient, maximum flange angle, and net wheel weight are analyzed according to the Kalker linear creep theory and Shen-Hedrick-Elkins creep theory, as shown in Sections 3 and 4. In addition, to realize fast calculation and application of critical 3D derailment coefficients, we present a simplified calculation method of critical derailment coefficients in Section 5. Finally, conclusions are addressed in Section 6.


## 2. Quasistatic Wheel Derailment Analysis Model

### 2.1. Coordinates of Wheel-Rail System

The wheel-rail contact relation is a link for coupling wheels and rails, while the coordinates of a wheel-rail system are crucial for establishing the wheel-rail relation. Figure1 shows the rail coordinates (OT-XT-YT-ZT), wheelset coordinates (Ow-Xw-Yw-Zw), and wheel-rail contact patch coordinates (OL,R-XL,R-YL,R-ZL,R). The rail coordinates are fixed on the central line of the rail and move along the central line of the rail at a certain speedV. The wheelset coordinates are fixed on the centroid of the wheelset. They move with the wheelset and have rotational and translational DOFs (degrees of freedom) in relation to the rail coordinates; the contact patch coordinates are fixed on the center of the wheel-rail contact patch and move with the wheelset.Figure 1
Coordinates of wheel-rail system.Transform relation between wheelset coordinates and rail coordinates is(1)iTjTkT=cos⁡ψ-sin⁡ψ0sin⁡ψcos⁡ϕcos⁡ψcos⁡ϕ-sin⁡ϕsin⁡ψsin⁡ϕcos⁡ψsin⁡ϕcos⁡ϕiwjwkw,where ψ is the wheelset yaw angle; ϕ is the wheelset rolling angle; [iT,jT,kT]T is the direction vector of the rail coordinates; [iw,jw,kw]T is the direction vector of the wheelset coordinates.Transform relation between the left/right wheel-rail contact patch coordinates and wheelset coordinates is(2)iwjwkw=1000cos⁡δL-sin⁡δL0sin⁡δLcos⁡δLiLjLkLiwjwkw=1000cos⁡δRsin⁡δR0-sin⁡δRcos⁡δRiRjRkR,where δL and δR are the wheel contact angles of the left and right wheel-rail contact points, respectively; [iL,jL,kL]T, [iR,jR,kR]T are the direction vectors of the left/right wheel-rail contact patch coordinates.The transform relation between the contact point coordinates and rail coordinates is the following after substituting (2) into (1):(3)iTjTkT=TL,RiL,RjL,RkL,R,where TL,R is the rotation matrix between the contact patch coordinates and rail coordinates (4)TL,R=cos⁡ψ-sin⁡ψcos⁡δL,R±sin⁡ψsin⁡δL,Rsin⁡ψcos⁡ϕcos⁡ψcos⁡ϕcos⁡δL,R∓sin⁡ϕsin⁡δL,R∓cos⁡ψcos⁡ϕsin⁡δL,R-sin⁡ϕcos⁡δL,Rsin⁡ψsin⁡ϕcos⁡ψsin⁡ϕcos⁡δL,R±cos⁡ϕsin⁡δL,R∓cos⁡ψsin⁡ϕsin⁡δL,R+cos⁡ϕcos⁡δL,R.Sinceψ and ϕ are small, TL,R can be simplified as follows: (5)TL,R=1-ψcos⁡δL,R±ψsin⁡δL,Rψcos⁡δL,R∓ϕsin⁡δL,R∓sin⁡δL,R-ϕcos⁡δL,R0ϕcos⁡δL,R±sin⁡δL,R∓ϕsin⁡δL,R+cos⁡δL,R.


### 2.2. Forces on Contact Patch and Critical Derailment Coefficient

Suppose the right wheel of the wheelset is about to derail and its flange contacts with the side surface of the rail; see Figure2. In this figure, F, Q, and P are the longitudinal force, lateral force, and vertical force on the contact point in the rail coordinate system, respectively. Tx, Ty, and N are the longitudinal creep force, lateral creep force, and normal force on the contact patch in the contact patch coordinate system, respectively.Figure 2
Forces on contact patch.According to the force equilibrium conditions on the wheel-rail contact patch as well as the transform relation between the rail coordinates and wheel-rail contact patch coordinates, we obtain(6)FQP+TRTxTy-N=0.The ratio of the wheel-rail lateral force to the vertical force at the contact patch can be obtained according to (5):(7)QP=-Txsin⁡ψcos⁡ϕ-Tycos⁡ψcos⁡ϕcos⁡δR+sin⁡θsin⁡δR+Ncos⁡ψcos⁡ϕsin⁡δR-sin⁡ϕcos⁡δR-Txsin⁡ψsin⁡ϕ-Tycos⁡ψsin⁡ϕcos⁡δR-cos⁡ϕsin⁡δR+Ncos⁡ψsin⁡ϕsin⁡δR+cos⁡ϕcos⁡δR.Since the yaw angle and rolling angle are small, then (7) can be further simplified as (8)QP=-Tycos⁡δR+Nsin⁡δRTysin⁡δR+Ncos⁡δR,where the ratio of lateral wheel-rail force to vertical force is called the derailment coefficient (DC) of the wheels. When the flange contact angle δR reaches the maximum, the gradient at this position during rail climbing will be the highest, the resistance for the wheels during rail climbing will be the highest and the wheels will be liable to derail when passing this position. Therefore, when δR is the maximum flange contact angle, the derailment coefficient is considered as the critical value of wheel derailment.Nadal thought that when the wheels were subject to a critical derailment state, the wheel-rail tangential force would satisfy Coulomb’s friction law. He used a two-dimensional (2D) model for analyzing derailment, so the contribution of the longitudinal creep force to the tangential force was ignored, then(9)Ty=μN,where μ is the wheel-rail friction coefficient.The formula of critical Nadal derailment coefficient can be obtained by substituting (9) in (8):(10)QPNadal=-μ+tan⁡δRμtan⁡δR+1.


## 2.1. Coordinates of Wheel-Rail System

The wheel-rail contact relation is a link for coupling wheels and rails, while the coordinates of a wheel-rail system are crucial for establishing the wheel-rail relation. Figure1 shows the rail coordinates (OT-XT-YT-ZT), wheelset coordinates (Ow-Xw-Yw-Zw), and wheel-rail contact patch coordinates (OL,R-XL,R-YL,R-ZL,R). The rail coordinates are fixed on the central line of the rail and move along the central line of the rail at a certain speedV. The wheelset coordinates are fixed on the centroid of the wheelset. They move with the wheelset and have rotational and translational DOFs (degrees of freedom) in relation to the rail coordinates; the contact patch coordinates are fixed on the center of the wheel-rail contact patch and move with the wheelset.Figure 1
Coordinates of wheel-rail system.Transform relation between wheelset coordinates and rail coordinates is(1)iTjTkT=cos⁡ψ-sin⁡ψ0sin⁡ψcos⁡ϕcos⁡ψcos⁡ϕ-sin⁡ϕsin⁡ψsin⁡ϕcos⁡ψsin⁡ϕcos⁡ϕiwjwkw,where ψ is the wheelset yaw angle; ϕ is the wheelset rolling angle; [iT,jT,kT]T is the direction vector of the rail coordinates; [iw,jw,kw]T is the direction vector of the wheelset coordinates.Transform relation between the left/right wheel-rail contact patch coordinates and wheelset coordinates is(2)iwjwkw=1000cos⁡δL-sin⁡δL0sin⁡δLcos⁡δLiLjLkLiwjwkw=1000cos⁡δRsin⁡δR0-sin⁡δRcos⁡δRiRjRkR,where δL and δR are the wheel contact angles of the left and right wheel-rail contact points, respectively; [iL,jL,kL]T, [iR,jR,kR]T are the direction vectors of the left/right wheel-rail contact patch coordinates.The transform relation between the contact point coordinates and rail coordinates is the following after substituting (2) into (1):(3)iTjTkT=TL,RiL,RjL,RkL,R,where TL,R is the rotation matrix between the contact patch coordinates and rail coordinates (4)TL,R=cos⁡ψ-sin⁡ψcos⁡δL,R±sin⁡ψsin⁡δL,Rsin⁡ψcos⁡ϕcos⁡ψcos⁡ϕcos⁡δL,R∓sin⁡ϕsin⁡δL,R∓cos⁡ψcos⁡ϕsin⁡δL,R-sin⁡ϕcos⁡δL,Rsin⁡ψsin⁡ϕcos⁡ψsin⁡ϕcos⁡δL,R±cos⁡ϕsin⁡δL,R∓cos⁡ψsin⁡ϕsin⁡δL,R+cos⁡ϕcos⁡δL,R.Sinceψ and ϕ are small, TL,R can be simplified as follows: (5)TL,R=1-ψcos⁡δL,R±ψsin⁡δL,Rψcos⁡δL,R∓ϕsin⁡δL,R∓sin⁡δL,R-ϕcos⁡δL,R0ϕcos⁡δL,R±sin⁡δL,R∓ϕsin⁡δL,R+cos⁡δL,R.


## 2.2. Forces on Contact Patch and Critical Derailment Coefficient

Suppose the right wheel of the wheelset is about to derail and its flange contacts with the side surface of the rail; see Figure2. In this figure, F, Q, and P are the longitudinal force, lateral force, and vertical force on the contact point in the rail coordinate system, respectively. Tx, Ty, and N are the longitudinal creep force, lateral creep force, and normal force on the contact patch in the contact patch coordinate system, respectively.Figure 2
Forces on contact patch.According to the force equilibrium conditions on the wheel-rail contact patch as well as the transform relation between the rail coordinates and wheel-rail contact patch coordinates, we obtain(6)FQP+TRTxTy-N=0.The ratio of the wheel-rail lateral force to the vertical force at the contact patch can be obtained according to (5):(7)QP=-Txsin⁡ψcos⁡ϕ-Tycos⁡ψcos⁡ϕcos⁡δR+sin⁡θsin⁡δR+Ncos⁡ψcos⁡ϕsin⁡δR-sin⁡ϕcos⁡δR-Txsin⁡ψsin⁡ϕ-Tycos⁡ψsin⁡ϕcos⁡δR-cos⁡ϕsin⁡δR+Ncos⁡ψsin⁡ϕsin⁡δR+cos⁡ϕcos⁡δR.Since the yaw angle and rolling angle are small, then (7) can be further simplified as (8)QP=-Tycos⁡δR+Nsin⁡δRTysin⁡δR+Ncos⁡δR,where the ratio of lateral wheel-rail force to vertical force is called the derailment coefficient (DC) of the wheels. When the flange contact angle δR reaches the maximum, the gradient at this position during rail climbing will be the highest, the resistance for the wheels during rail climbing will be the highest and the wheels will be liable to derail when passing this position. Therefore, when δR is the maximum flange contact angle, the derailment coefficient is considered as the critical value of wheel derailment.Nadal thought that when the wheels were subject to a critical derailment state, the wheel-rail tangential force would satisfy Coulomb’s friction law. He used a two-dimensional (2D) model for analyzing derailment, so the contribution of the longitudinal creep force to the tangential force was ignored, then(9)Ty=μN,where μ is the wheel-rail friction coefficient.The formula of critical Nadal derailment coefficient can be obtained by substituting (9) in (8):(10)QPNadal=-μ+tan⁡δRμtan⁡δR+1.


## 3. Calculation of Creep Force

Under quasistatic assumption, when the wheelset speed and acceleration are ignored, the wheel-rail rigid creepage can be calculated according to the following [14]:(11)ξx=1-rwRr0ξy=-ψsec⁡δRξs=-sin⁡δRr0,where ξx, ξy, and ξs are the longitudinal creepage, lateral creepage, and spin creepage, respectively. rwR is the real rolling radius of the wheels at the right wheel-rail contact point, while r0 is the nominal rolling radius of the wheels. According to (11), the creepage obtained under quasistatic assumption relates only to the wheel-rail contact geometric parameter.According to the Kalker linear creep theory [15], the relation between the creep force and creepage can be expressed as follows: (12)Tx′Ty′M′=EabC11000C22abC230-abC23abC33ξxξyξs,where Tx′, Ty′, and M′ are the Kalker linear longitudinal/lateral creep force and spin creep force, respectively; E is Young’s modulus of the wheel and rail materials; a and b are the longitudinal and lateral semiaxle lengths of the Hertz elliptical contact patch, respectively; C11, C22, C23, and C33 are the Kalker creep coefficients relating to the ratio of a and b, which are shown in Table 1.Table 1
Relation between creep coefficient anda/b.
a
/
b
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
—
C11
3.31
3.37
3.44
3.53
3.62
3.72
3.81
3.91
4.01
—
C22
2.52
2.63
2.75
2.88
3.01
3.14
3.28
3.41
3.54
—
C23
0.473
0.603
0.715
0.823
0.929
1.03
1.14
1.25
1.36
—
C33
8.82
4.27
2.96
2.32
1.93
1.68
1.50
1.37
1.27
—
b
/
a
1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
C11
4.12
4.22
4.36
4.54
4.78
5.10
5.57
6.34
7.78
11.7
C22
3.67
3.81
3.99
4.21
4.50
4.90
5.48
6.40
8.14
12.8
C23
1.47
1.59
1.75
1.95
2.23
2.62
3.24
4.23
6.63
14.6
C33
1.19
1.11
1.04
0.965
0.892
0.819
0.747
0.674
0.601
0.526According to the Hertz contact theory, the longitudinal/lateral semiaxlea and b of the elliptical wheel-rail contact patch are(13)a=m31-σ22EA+BN1/3b=n31-σ22EA+BN1/3,where σ is Poisson’s ratio, m and n are the semiaxle correction coefficients of the ellipses, which can obtained according to θ of (14) and Table 2:(14)θ=arccos⁡B-AA+B,where (15)A+B=121RTR+1RwR+1rwRB-A=121RTR-1RwR-1rwR,where RTR and RwR are the radii of the rail and wheel profile at the wheel-rail contact.Table 2
Relation betweenm/n and θ.
θ
0°
10°
20°
30°
35°
40°
45°
50°
m
∞
6.612
3.778
2.731
2.397
2.130
1.926
1.754
n
0
0.319
0.408
0.493
0.530
0.567
0.604
0.641
θ
55°
60°
65°
70°
75°
80°
85°
90°
m
1.611
1.486
1.378
1.284
1.202
1.128
1.061
1.000
n
0.678
0.717
0.759
0.802
0.846
0.893
0.944
1.000The Kalker linear creep force model is only suitable when the creep is low since it is supposed that no relative slip occurs in the contact patch. Under flange contact, the creepage increases, and even a full slip state can occur in the contact patch, so the linear creep force must be corrected. In this paper, the Shen-Hedrick-Elkins nonlinear creep model [16] has been used for correcting the wheel-rail creep force; see the following:(16)TxTyM=TsTs′Tx′Ty′M′,where Ts′ is the comprehensive creep force calculated according to the Kalker linear creep theory; see (17); Ts is the corrected comprehensive creep force; see (18).(17)Ts′=Tx′+Ty′(18)Ts=μNTs′μN-13Ts′μN2+127Ts′μN3,Ts′≤3μNμN,Ts′>3μN.According to the wheel-rail creep force calculation model, when the wheel-rail contact parameter and wheelset yaw angle are given, the creep force only relates to the normal force. Since contact point vertical force equilibrium can be obtained based on (6), then (19)TyNϕcos⁡δR+sin⁡δR+Ncos⁡δR-ϕsin⁡δR-P=0.The following can be obtained when the rolling angle is ignored:(20)TyNsin⁡δR+Ncos⁡δR-P=0.When the wheel-rail vertical forceP is given, the wheel-rail normal force N can be calculated by way of iteration according to (20) and Shen-Hedrick-Elkins nonlinear creep model and then the creep force and critical derailment coefficient can be obtained.


## 4. Results and Discussion

Table3 shows the wheel-rail contact parameters obtained under quasistatic critical wheel derailment conditions. The wheel-rail force and critical derailment coefficient as the wheelset yaw angle changes are shown in Figures 3 and 4.Table 3
Wheel-rail contact parameters obtained under critical derailment conditions.
P
E
δ
R
σ
r
0
r
R
R
W
R
R
T
R
μ
70 kN
2.06 × 105 Mpa
70°
0.25
0.4575 m
0.475 m
0.014 m
0.013 m
0.25Figure 3
Wheel-rail force of contact patch under critical derailment.Figure 4
Critical derailment coefficient.As shown in Figure3, when the wheelset yaw angle is a negative value, the wheel-rail normal force will significantly increase with the increase of the wheelset yaw angle. When the yaw angle increases from 0° to −1.5°, the wheel-rail normal force will linearly increase. With the increase of the wheelset yaw angle, the longitudinal creep force will increase and then decrease and reach the maximum when the yaw angle is about −0.5°. The lateral creep force is a negative value and increases significantly with the increase of the wheelset yaw angle, indicating that it is opposite to the wheel derailment direction shown in Figure 2 and hinders the wheel derailment. This mode is called wheel-slide derailment. When the wheelset yaw angle is a positive value, the wheel-rail normal force and longitudinal creep force, to a certain extent, decrease with the increase of the wheelset yaw angle; however when the wheelset yaw angle is over 1.5°, their amplitudes basically do not change. Moreover, when the positive wheelset yaw angle becomes high, the longitudinal creep force basically tends toward 0. The lateral creep force, to a certain extent, increases with the increase of the wheelset yaw angle, and its amplitude is always positive, coinciding with the wheel derailment direction, indicating that it can prompt wheel derailment. This mode is called wheel-climb derailment.As shown in Figure4, the critical derailment coefficient under a negative wheelset yaw angle is significantly higher than that under a positive yaw angle. This means that wheels are liable to derail under positive yaw angles, because the influence of the lateral creep forces on wheel derailment under negative wheelset yaw angles is different from that under positive yaw angles. In addition, Figure 4 also shows the ratio of the critical derailment coefficient to the limited Nadal derailment coefficient. It can be seen that when the positive wheelset yaw angle is over 1.5°, their amplitudes are almost the same, because at this time the longitudinal creep force amplitude is close to 0, and the lateral creep force amplitude is close to the Coulomb friction force. However, when the positive wheelset yaw angle is low or the angle is negative, the limited Nadal derailment coefficient is conservative to a certain extent, because the longitudinal creep force cannot be ignored, and the direction of the lateral creep force must to be considered.For given profiles of wheels and rails, parameters such asσ, r0, RwR, and RTR do not change much or influence wheel derailment much during operation. However, the wheel-rail friction coefficient can be significantly influenced by changes of season and weather, and the wheel flange can wear to some extent, causing change of the maximum flange contact angle. In addition, the axle load also significantly changes due to different vehicle loads. Hence, the study should focus on the influence of the above three parameters on the critical derailment coefficient.Equation (8) can be simplified to be the same as that for calculating the Nadal derailment coefficient for the sake of easy analysis: (21)QP=-μe+tan⁡δRμetan⁡δR+1,where μe=Ty/N, called the equivalent friction coefficient.Figures5 and 6 show the influence of changing wheel-rail friction coefficient on the critical derailment coefficient. The equivalent friction coefficient can be decreased by decreasing the wheel-rail friction coefficient; that is, the ratio of lateral creep force to normal force can be decreased. However, the direction of the lateral creep force is the same as that of the wheelset yaw angle, so, under a positive wheelset yaw angle, the positive effect of the lateral creep force on wheel derailment can be reduced by decreasing the wheel-rail friction coefficient, which is beneficial for preventing wheel derailment. However, under a negative wheelset yaw angle, the decrease of wheel-rail friction coefficient can obstruct the positive effect of the lateral creep force on wheel derailment, which is detrimental for preventing wheel derailment. Therefore, care should be taken to prevent wheel derailment through decreasing the wheel-rail friction coefficient. Though an over low friction coefficient is good for reducing derailment risks under a positive yaw angle, the difference between the limited wheel derailment coefficient under a negative yaw angle and the critical derailment coefficient under a positive yaw angle can be reduced, causing possible wheel-slide derailment.Figure 5
Equivalent friction coefficient under different friction coefficients.Figure 6
Critical derailment coefficient under different friction coefficients.Figures7 and 8 show the influence of changing maximum flange contact angle on the equivalent friction coefficient and critical derailment coefficient. When the yaw angle is kept within -1.5°~1.5°, the equivalent friction coefficient will increase to a certain extent with the increase of the maximum flange contact angle. However, when the absolute value of the wheelset yaw angle is over 1.5°, the equivalent friction coefficient will basically not change, because when the wheelset yaw angle becomes high, the longitudinal creep force will be close to 0, and the lateral creep force will be close to the Coulomb friction force. As shown in the figures, no matter whether the wheelset yaw angle is positive or negative, the critical derailment coefficient can increase with the increase of the maximum flange contact angle, reducing the wheel derailment risk.Figure 7
Equivalent friction coefficient under different maximum flange contact angles.Figure 8
Critical derailment coefficient under different maximum flange contact angles.Figures9 and 10 show the influence of the wheel load on the equivalent friction coefficient and critical derailment coefficient. It can be seen that the change of wheel load basically has no influence on the equivalent friction coefficient and critical derailment coefficient.Figure 9
Equivalent friction coefficient under different wheel loads.Figure 10
Critical derailment coefficient under different wheel loads.


## 5. Simplified Calculation Method

According to (18), the ratio of the Kalker linear comprehensive creep force to 3 times the Coulomb friction force is a key factor in determining whether the creep of wheel-rail contact patches reaches saturation. The formula presented above depends on many parameters and cannot facilitate further understanding of rail-wheel derailment. In this chapter, we present a simplified calculation model which is easy to use while maintaining accuracy.


### 5.1. Presentation of Simplified Calculation Method

Figure11 shows how the ratio of the Kalker linear comprehensive creep force to 3 times the Coulomb friction force changes with the yaw angle under different friction coefficients, maximum flange contact angles, and net wheel loads. According to the analysis, when the wheelset is subject to critical derailment, all ratios of the comprehensive Kalker linear creep force to 3 times the Coulomb friction force are higher than 1, indicating that the creep force within the wheel-rail contact patch at the flange contact side reaches saturation at that moment.Figure 11
Ratios of the Kalker linear creep force to 3 times the Coulomb friction force. (a) Different friction coefficients, (b) different maximum flange contact angles, and (c) different net wheel loads.
(a)
(b)
(c)According to the Shen-Hedrick-Elkins creep theory, when the creep force within the contact patch reaches saturation, the relation between the wheel-rail longitudinal/lateral creep force and normal force is(22)Tx2+Ty2=μN.Suppose that the included angle betweenTx and Ty is α; then (23)tan⁡α=TyTxTx=μNcos⁡αTy=μNsin⁡α.Substitute (23) into (21), then the critical derailment coefficient can be expressed as follows: (24)QP=-μe+tan⁡δRμetan⁡δR+1μe=μsin⁡α.The critical derailment coefficient can be calculated according to (24) if the included angle between the longitudinal and lateral creep forces can be obtained. The ratio of the lateral creep force to the longitudinal creep force under critical derailment can be obtained according to (12):(25)TyTx=EabC22ξy+Eab3/2C23ξsEabC11ξx=C22ξy+ab1/2C23ξsC11ξx.The following can be obtained by substituting (11) and (13) into (25): (26)TyTx=K1ψ+K2N1/3,where(27)K1=C22sec⁡δRC11rwR/r0-1,K2=C23sin⁡δRC11rwR-r0mn1/231-σ22EA+B1/3.According to (26), when the contact parameters under critical wheel derailment are known, K1 and K2 can be determined, and the ratio of the lateral creep force to the longitudinal creep force only relates to the wheelset yaw angle and normal force. Figure 12 shows the amplitudes of K1 and K2 under different maximum flange contact angles. It can be seen that the amplitude of K1 is far higher than that of K2, and the maximum flange contact angle can significantly influence K1 but does not influence K2 much.Figure 12
Amplitudes ofK1 and K2 under different maximum flange contact angles. (a) K1; (b) K2.
(a)
(b)In addition to the influence of wheel-rail contact parameters, the ratio of the lateral creep force to the longitudinal creep force can also be influenced by the normal force. According to Section4, the normal force is mainly influenced by the friction coefficient and maximum flange contact angle. Figure 13 shows how the amplitudes of K1ψ and K2N1/3 of (26) change with the wheelset yaw angle when the maximum flange contact angles are 60° and 70°, respectively. According to Figure 13, K1 is far higher than K2, causing the amplitude of K1ψ to be significantly higher than that of K2N1/3, so it can be deduced that the ratio of the longitudinal creep force and lateral creep force has a strong linear relation to the wheelset yaw angle rather than the friction coefficient.Figure 13
Amplitudes ofK1ψ and K2N1/3 under different friction coefficients. (a) δR = 60°; (b) δR = 70°.
(a)
(b)To verify the above deduction, Figure14 shows how the ratio of the lateral creep force to the longitudinal creep force changes with the wheelset yaw angle under different friction coefficients when the maximum flange contact angles are 60° and 70°, respectively. Under a positive wheelset yaw angle, the friction coefficient basically has no influence on their ratio. However, under a negative wheelset yaw angle, the friction coefficient can influence their ratio to some extent, because under the positive attack angle the wheel-rail normal force is low and does not change much with the friction coefficient, but, under the negative wheelset yaw angle, the wheel-rail normal force is significantly higher than that under the positive wheelset yaw angle. As a whole, K1 is far higher than K2 under critical derailment, so the ratio of the lateral creep force to the longitudinal creep force changes with the wheelset yaw angle linearly and intensely and is not influenced much by the friction coefficient, indicating that the deduction of this paper is reasonable and correct.Figure 14
Ratio of lateral creep force to longitudinal creep force. (a)δR = 60°; (b) δR = 70°.
(a)
(b)In order to simplify the calculation of critical derailment coefficient, based on the analysis of the strong linear relation between the lateral creep force and longitudinal creep force under critical derailment, this paper proposed that the linear relationship, as shown in (27), can be adopted to fit the ratio of the lateral creep force to the longitudinal creep force. Figure 15 shows the comparison of this fitting to raw simulation data. (28)TyTx=Cψ+D.Figure 15
Comparison between linear fitting and theoretical result for ratio of lateral creep force to longitudinal creep force. (a)δR = 60°; (b) δR = 70°.
(a)
(b)The simplified evaluation formula of 3D derailment coefficients limit under quasistatic assumption can be obtained by substituting (28) into (24), as shown in(29)QP=-μe+tan⁡δRμetan⁡δR+1μe=μCψ+D1+Cψ+D2.


### 5.2. Verification of Accuracy

To verify the difference between the critical derailment coefficients calculated according to the simplified derailment coefficient calculation method and the accurate calculation method, Figure16 shows the comparison between them under different wheelset yaw angle and friction coefficients. The errors between them is calculated according to (30), as shown in Figure 17. (30)Error  rate=formula29-formula8formula8100%.Figure 16
Comparison between critical derailment coefficients calculated according to the simplified formula and accurate calculation formula. (a)δR = 60°; (b) δR = 70°.
(a)
(b)Figure 17
Analysis of error between critical derailment coefficients calculated according to the simplified formula and accurate calculation formula. (a)δR = 60°; (b) δR = 70°.
(a)
(b)As shown in Figure16, the critical derailment coefficients calculated according to the simplified formula (29) mentioned in this paper and the quasistatic formula (8) under critical derailment conditions are close to each other, especially when the wheelset yaw angle is high, and the two results approach each other. As shown in Figure 17, with increase of the friction coefficient, the error rate between the critical derailment coefficients calculated according to the simplified formula (29) and the quasistatic formula (8) under critical derailment conditions will increase to a certain extent; however, both of the errors are controllable within -5%~5%, meeting the requirement of engineering application. Thus the simplified formula for wheel derailment evaluation criterion under quasistatic assumption is reasonable, feasible, and reliable.


### 5.3. Simplified Formula Calculation Table

As is well known, wheel-rail wear is unavoidable with the long-term operation of vehicles. The abrasion of wheel flanges and rail gauge corners can cause decrease of the maximum flange contact angle, influencing the operational safety. In order to make the simplified calculation method of critical derailment coefficients more practical, the linear fitting parameters for maximum flange contact angles within60°~70° are calculated in this paper; see Table 4.Table 4
Simplified calculation parameters for critical derailment coefficients.
Maximum flange contact angle (°)
C
D
60
1.00178
0.19069
62
1.06522
0.20025
64
1.13845
0.21123
66
1.22337
0.22473
68
1.3206
0.24474
70
1.4535
0.24839During practical application, the proper parameters can be obtained according to the measured maximum wheel flange contact angle and the data in Table4, and the corresponding critical derailment coefficient can be calculated according to the simplified calculation formula, reducing the working load and promoting the application of 3D derailment coefficient limit.


## 5.1. Presentation of Simplified Calculation Method

Figure11 shows how the ratio of the Kalker linear comprehensive creep force to 3 times the Coulomb friction force changes with the yaw angle under different friction coefficients, maximum flange contact angles, and net wheel loads. According to the analysis, when the wheelset is subject to critical derailment, all ratios of the comprehensive Kalker linear creep force to 3 times the Coulomb friction force are higher than 1, indicating that the creep force within the wheel-rail contact patch at the flange contact side reaches saturation at that moment.Figure 11
Ratios of the Kalker linear creep force to 3 times the Coulomb friction force. (a) Different friction coefficients, (b) different maximum flange contact angles, and (c) different net wheel loads.
(a)
(b)
(c)According to the Shen-Hedrick-Elkins creep theory, when the creep force within the contact patch reaches saturation, the relation between the wheel-rail longitudinal/lateral creep force and normal force is(22)Tx2+Ty2=μN.Suppose that the included angle betweenTx and Ty is α; then (23)tan⁡α=TyTxTx=μNcos⁡αTy=μNsin⁡α.Substitute (23) into (21), then the critical derailment coefficient can be expressed as follows: (24)QP=-μe+tan⁡δRμetan⁡δR+1μe=μsin⁡α.The critical derailment coefficient can be calculated according to (24) if the included angle between the longitudinal and lateral creep forces can be obtained. The ratio of the lateral creep force to the longitudinal creep force under critical derailment can be obtained according to (12):(25)TyTx=EabC22ξy+Eab3/2C23ξsEabC11ξx=C22ξy+ab1/2C23ξsC11ξx.The following can be obtained by substituting (11) and (13) into (25): (26)TyTx=K1ψ+K2N1/3,where(27)K1=C22sec⁡δRC11rwR/r0-1,K2=C23sin⁡δRC11rwR-r0mn1/231-σ22EA+B1/3.According to (26), when the contact parameters under critical wheel derailment are known, K1 and K2 can be determined, and the ratio of the lateral creep force to the longitudinal creep force only relates to the wheelset yaw angle and normal force. Figure 12 shows the amplitudes of K1 and K2 under different maximum flange contact angles. It can be seen that the amplitude of K1 is far higher than that of K2, and the maximum flange contact angle can significantly influence K1 but does not influence K2 much.Figure 12
Amplitudes ofK1 and K2 under different maximum flange contact angles. (a) K1; (b) K2.
(a)
(b)In addition to the influence of wheel-rail contact parameters, the ratio of the lateral creep force to the longitudinal creep force can also be influenced by the normal force. According to Section4, the normal force is mainly influenced by the friction coefficient and maximum flange contact angle. Figure 13 shows how the amplitudes of K1ψ and K2N1/3 of (26) change with the wheelset yaw angle when the maximum flange contact angles are 60° and 70°, respectively. According to Figure 13, K1 is far higher than K2, causing the amplitude of K1ψ to be significantly higher than that of K2N1/3, so it can be deduced that the ratio of the longitudinal creep force and lateral creep force has a strong linear relation to the wheelset yaw angle rather than the friction coefficient.Figure 13
Amplitudes ofK1ψ and K2N1/3 under different friction coefficients. (a) δR = 60°; (b) δR = 70°.
(a)
(b)To verify the above deduction, Figure14 shows how the ratio of the lateral creep force to the longitudinal creep force changes with the wheelset yaw angle under different friction coefficients when the maximum flange contact angles are 60° and 70°, respectively. Under a positive wheelset yaw angle, the friction coefficient basically has no influence on their ratio. However, under a negative wheelset yaw angle, the friction coefficient can influence their ratio to some extent, because under the positive attack angle the wheel-rail normal force is low and does not change much with the friction coefficient, but, under the negative wheelset yaw angle, the wheel-rail normal force is significantly higher than that under the positive wheelset yaw angle. As a whole, K1 is far higher than K2 under critical derailment, so the ratio of the lateral creep force to the longitudinal creep force changes with the wheelset yaw angle linearly and intensely and is not influenced much by the friction coefficient, indicating that the deduction of this paper is reasonable and correct.Figure 14
Ratio of lateral creep force to longitudinal creep force. (a)δR = 60°; (b) δR = 70°.
(a)
(b)In order to simplify the calculation of critical derailment coefficient, based on the analysis of the strong linear relation between the lateral creep force and longitudinal creep force under critical derailment, this paper proposed that the linear relationship, as shown in (27), can be adopted to fit the ratio of the lateral creep force to the longitudinal creep force. Figure 15 shows the comparison of this fitting to raw simulation data. (28)TyTx=Cψ+D.Figure 15
Comparison between linear fitting and theoretical result for ratio of lateral creep force to longitudinal creep force. (a)δR = 60°; (b) δR = 70°.
(a)
(b)The simplified evaluation formula of 3D derailment coefficients limit under quasistatic assumption can be obtained by substituting (28) into (24), as shown in(29)QP=-μe+tan⁡δRμetan⁡δR+1μe=μCψ+D1+Cψ+D2.


## 5.2. Verification of Accuracy

To verify the difference between the critical derailment coefficients calculated according to the simplified derailment coefficient calculation method and the accurate calculation method, Figure16 shows the comparison between them under different wheelset yaw angle and friction coefficients. The errors between them is calculated according to (30), as shown in Figure 17. (30)Error  rate=formula29-formula8formula8100%.Figure 16
Comparison between critical derailment coefficients calculated according to the simplified formula and accurate calculation formula. (a)δR = 60°; (b) δR = 70°.
(a)
(b)Figure 17
Analysis of error between critical derailment coefficients calculated according to the simplified formula and accurate calculation formula. (a)δR = 60°; (b) δR = 70°.
(a)
(b)As shown in Figure16, the critical derailment coefficients calculated according to the simplified formula (29) mentioned in this paper and the quasistatic formula (8) under critical derailment conditions are close to each other, especially when the wheelset yaw angle is high, and the two results approach each other. As shown in Figure 17, with increase of the friction coefficient, the error rate between the critical derailment coefficients calculated according to the simplified formula (29) and the quasistatic formula (8) under critical derailment conditions will increase to a certain extent; however, both of the errors are controllable within -5%~5%, meeting the requirement of engineering application. Thus the simplified formula for wheel derailment evaluation criterion under quasistatic assumption is reasonable, feasible, and reliable.


## 5.3. Simplified Formula Calculation Table

As is well known, wheel-rail wear is unavoidable with the long-term operation of vehicles. The abrasion of wheel flanges and rail gauge corners can cause decrease of the maximum flange contact angle, influencing the operational safety. In order to make the simplified calculation method of critical derailment coefficients more practical, the linear fitting parameters for maximum flange contact angles within60°~70° are calculated in this paper; see Table 4.Table 4
Simplified calculation parameters for critical derailment coefficients.
Maximum flange contact angle (°)
C
D
60
1.00178
0.19069
62
1.06522
0.20025
64
1.13845
0.21123
66
1.22337
0.22473
68
1.3206
0.24474
70
1.4535
0.24839During practical application, the proper parameters can be obtained according to the measured maximum wheel flange contact angle and the data in Table4, and the corresponding critical derailment coefficient can be calculated according to the simplified calculation formula, reducing the working load and promoting the application of 3D derailment coefficient limit.


## 6. Conclusion

In this paper, force analysis for 3D critical derailment is carried out, and a simplified method for calculating the critical derailment coefficient is deduced based on Kalker linear creep theory and Shen-Hedrick-Elkins creep theory. Based on the theoretical analysis, the following conclusions are reached.(1) Under critical wheel derailment conditions, derailments modes of wheel-climb and wheel-slide can be distinguished according to the direction of the lateral creep force. Since the direction of lateral creep force is the same as that of wheel derailment in case of wheel-climb derailment, which can enhance the probability of derailment. Thus, derailment of wheel-climb is more likely to happen and particular attentions should be paid to it.(2) When the wheelset yaw angle is small, the longitudinal creep force will account for a large portion of the tangential force of wheel-rail contact patches. That means if its influence on wheel derailment is neglected, the calculated critical derailment coefficient will be relatively low, further proving that critical Nadal derailment coefficients are conservative in evaluating wheel derailment under small wheelset yaw angles.(3) Wheel-rail friction coefficients and maximum wheel flange contact angles can greatly influence the wheel derailment coefficient limit. Decrease of wheel-rail friction coefficient is helpful in avoiding wheel-climb derailment, but it can cause a decrease of the limit of wheel-slide derailment at the same time. Thus, thorough calculation should be carried out when trying to prevent wheel-slide derailment by decreasing the wheel-rail friction coefficient in order to prevent wheel-climb derailment.(4) Under critical conditions of quasistatic wheel derailment, the tangential force within a wheel-rail contact patch calculated under different friction coefficients, maximum flange contact angles, and net wheel loads is analyzed according to the Kalker linear creep theory and Shen-Hedrick-Elkins creep theory. Results show that the creep force within the wheel-rail contact patch has saturated and matched the Coulomb friction law under these conditions.(5) Hence a simplified method for calculating critical derailment coefficients is presented through linearly fitting the ratio of the lateral creep force to the longitudinal creep force, as inspired by numerous calculation. The error rate for the results calculated according to the simplified formula can be kept within ±5% when compared to that calculated according to the accurate formula, indicating that it is reliable and meets requirements of engineering application.



---

*Source: 7910753-2018-01-17.xml*