# Mechanism Analysis of Spalling Defect on Rail Surface under Rolling Contact Conditions

Authors: Rongshan Yang; Shihao Cao; Weixin Kang; Jiali Li; Xiaoyu Jiang
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

Under the wheel/rail contact loading conditions, the microcracks on the rail surface propagate, leading to spalling defect or rail fracture and threatening the travelling safety of high-speed railway directly. In order to analyze the mechanism of the crack propagation on the rail surface, the calculation model of the wheel/rail contact fatigue was established, and the variation of the stress intensity factor at the crack tip when the crack length was increased from 0.1 mm to 2 mm was obtained. Based on the mixed-mode fracture criterion and Paris growth theory, the mechanism of the crack propagation on the rail surface was analyzed. The results show that when the microcrack grows to macrocrack, the mode of the fatigue crack on the rail surface is mixed including sliding mode and open mode. With the increase of the crack length, the stress intensity factorKI increases first and then decreases gradually, and the relative dangerous location of the open-mode crack moves from the inner edge of the contact area to the outer edge, while the factor KII is increasing during the whole propagation process, and the relative dangerous location of the sliding-mode crack remains unchanged basically. The main failure mode of crack is open during the initial stage and then transforms into sliding mode with the crack length increasing. The crack tends to propagate upward and leads to spalling defect when the crack length is between 0.3 and 0.5 mm. This propagation path is basically identical with the spalling path of the service rail. The research results will provide a basis for improving the antifatigue performance of rail and establishing the grinding procedure.


---

## Body

## 1. Introduction

As the travelling base of the high-speed railway, rail plays a role of supporting and guiding the running of the train. However, due to the repeated extrusion caused by the contact stress between the wheel and the rail, microcrack initiation and propagation will appear on the surface of the rail. As a result, there will be rail head spalling or even fracture. With the freight volume and the train speed increasing, rolling contact fatigue (RCF) has become more and more serious, which affects the traffic safety directly [1]. For example, in 2000, a serious derailment accident happened to a British high-speed train because of the fractured rail [2]. Coincidentally, there were also rail fractures caused by the two oblique cracks in the Guangzhou-Shenzhen high-speed railway in 2002 [3].Many researches about the growth path of the wheel/rail contact fatigue have been reported. The review by Smith et al. [2, 4] indicated that there was an acute angle between the crack initiation and the direction of train running during the initial stage, and the growth angle increased with the increase of the depth, and, finally, the rail fractured in the direction perpendicular to the surface, as shown in Figure 1(a). Chen et al. [5] studied the spalling path of the rail head through slice analysis and concluded that when the cracks grew to a certain depth, they would change their propagation direction and then spread approximately along the right-about with train running direction towards the rail surface; consequently, spalling was formed in the shape of barb, as shown in Figure 1(b). Through analyzing the chemical composition and mechanical properties of rail with spalling defects, Zhou et al. [6] proposed that the generation of spalling defects had no direct relation with the rail quality, which was mainly caused by the excessive contact stress.Figure 1
Growth path of the rolling contact fatigue cracks: (a) fracture and (b) spalling on the rail surface.
(a) (b)In the aspect of experimental study, through the convex-concave contact configuration of a shaft running inside a ring, Christoph et al. [7] performed a laboratory test simulating the fatigue inclined cracks on the rail surface. Huang et al. [8, 9] characterized the wear of the rail surface via WR-1 wheel/rail rolling wear testing apparatus. Based on the field observation of rolling contact fatigue, Tao [10] investigated the cause of wheel tread spalling preliminarily by the use of shakedown theory.With respect to the mechanism of microcracks initiation, Langueh et al. [11, 12] proposed that stress concentration in the local region was the main cause of fatigue cracks initiation under the action of the contact stress. Dubourg et al. [13, 14] argued that shear stress played a significant role in the initiation of cracks. Jiang and Jin [15, 16] proposed that the large residual tensile stress on the rail surface could possibly be an important factor of microcracks initiation. As for the crack propagation mechanism, Ringsberg et al. [17–19] studied the microcracks propagation mechanism of the rail surface based on the Hertz theory and the wheel/rail tribology and found that the propagation of microcracks on surface was mainly open mode. With regard to the long cracks, Seo et al. [4, 20, 21] investigated that sliding propagation was their major mode, and their propagation rate increased as the friction coefficient between the wheel and the rail increased.Currently, most studies about the propagation mechanism of cracks are focused on macro long cracks, which ignore the stage when microcrack transforms into macrocrack. In fact, this stage is the main part of the entire propagation process, and most rail spalling occurs in this stage as well [10]. Therefore, based on the assumption that the wheel was rolling along the rail, this paper analyzed the propagation characteristics of the fatigue cracks on the rail surface from microcrack to macrocrack and made the fatigue spalling mechanism of rails more clear, which could provide the basis for the improvement of the rail fatigue performance and the establishment of grinding procedure. In this way the operation costs as well as risks will be cut down.


## 2. Force Analysis of Wheel/Rail Contact

With the development of production technology in rails, it is hard to find rolling contact fatigue (RCF) defects caused by nonmetallic inclusions in rail and hydrogen induced crack. However, in the high-speed passenger dedicated lines and heavy haul lines, the RCF defects that occurred on the rail head surface become more and more serious, these defects are mainly caused by the huge contact stress act on the rail surface [6].


### 2.1. Hertz Contact Theory

As the theoretical basis of calculating the wheel/rail contact stress, the Hertz contact theory assumes that the wheel/rail contact area is an ellipse, and the maximum contact pressureqmax is in the center of the ellipse, which is given by [18].(1)qmax=3P2πab,where P is the force act on the rail and a is the semimajor axis of the elliptical contact area, while b is the semiminor axis of the elliptical contact area. Then the values of a and b can be obtained by (2)a=m31-υ2P2EA+B3b=nma,where E is Young’s modulus and υ is Poisson’s ratio, while m and n denote the calculating coefficients of wheel/rail contact stress [22].As for the coefficientsA and B, they can be obtained through (3)A+B=121R1+1R2B-A=121R1-1R2,where R1 is the radius of the wheel, R2 is the crown radius of the rail head.The distribution of the contact pressure acting on the rail is(4)qx=qmax1-x2a2,The tangential friction force in the contact area can be obtained according to Coulomb’s law of friction [23]. When the wheel-rail contact is in an adhesion condition, the tangential friction force is described as(5)τx≤μqx.When the contact is in a sliding condition, the tangential friction force is(6)τx=μqx,where τ(x) is the tangential friction force in the contact area and μ is the friction coefficient of wheel and rail.See Figure2, where  X is the longitudinal direction of the rail, Y is the depth direction of the rail, a is the semimajor axis of the contact area, x is the distance from the center of the contact area to the crack, P is the force act on the rail, v is the rolling direction of the wheel, L is the crack length, θ is the angle between crack and rolling direction, σx is the horizontal normal stress of the rail surface, σy is the vertical normal stress of the rail surface, and τxy is the shear stress of the rail surface. And the numbers 1, 2, and 3 denote the rail microelements bearing the forces when the wheel is close to different positions, respectively.Figure 2
Distribution of contact forces and stress state of rail surface under the action of rolling wheel.


### 2.2. Stress State of the Rail Surface

When a wheel is rolling along the rail, the area of the contact area is approximately 100 mm2, and the load act on the contact area is about 50~200 kN, which could result in a huge contact stress [24]. Under the action of the normal contact pressure and the tangential friction force, the stress state of the rail surface can be seen from Figure 2. When the wheel rolls from position 1 to position 3, the horizontal tensile stress and shear stress will have periodic changes in direction. In general, this combined alternating stress is considered to be the main cause of the crack initiation and propagation on the rail surface.


### 2.3. Analysis of the Crack State

When the wheel is close to or away from the crack, the extrusion effect on the rail surface leads to the stretch of the surface material near the edges of the contact area. At the same time, the surface cracks of the rail will open under the action of tensile stress. While the wheel pressure acts on the crack, the crack will be closed because of the mutual extrusion of the crack surfaces, as shown in Figure3.Figure 3
State of crack on the rail surface.


## 2.1. Hertz Contact Theory

As the theoretical basis of calculating the wheel/rail contact stress, the Hertz contact theory assumes that the wheel/rail contact area is an ellipse, and the maximum contact pressureqmax is in the center of the ellipse, which is given by [18].(1)qmax=3P2πab,where P is the force act on the rail and a is the semimajor axis of the elliptical contact area, while b is the semiminor axis of the elliptical contact area. Then the values of a and b can be obtained by (2)a=m31-υ2P2EA+B3b=nma,where E is Young’s modulus and υ is Poisson’s ratio, while m and n denote the calculating coefficients of wheel/rail contact stress [22].As for the coefficientsA and B, they can be obtained through (3)A+B=121R1+1R2B-A=121R1-1R2,where R1 is the radius of the wheel, R2 is the crown radius of the rail head.The distribution of the contact pressure acting on the rail is(4)qx=qmax1-x2a2,The tangential friction force in the contact area can be obtained according to Coulomb’s law of friction [23]. When the wheel-rail contact is in an adhesion condition, the tangential friction force is described as(5)τx≤μqx.When the contact is in a sliding condition, the tangential friction force is(6)τx=μqx,where τ(x) is the tangential friction force in the contact area and μ is the friction coefficient of wheel and rail.See Figure2, where  X is the longitudinal direction of the rail, Y is the depth direction of the rail, a is the semimajor axis of the contact area, x is the distance from the center of the contact area to the crack, P is the force act on the rail, v is the rolling direction of the wheel, L is the crack length, θ is the angle between crack and rolling direction, σx is the horizontal normal stress of the rail surface, σy is the vertical normal stress of the rail surface, and τxy is the shear stress of the rail surface. And the numbers 1, 2, and 3 denote the rail microelements bearing the forces when the wheel is close to different positions, respectively.Figure 2
Distribution of contact forces and stress state of rail surface under the action of rolling wheel.


## 2.2. Stress State of the Rail Surface

When a wheel is rolling along the rail, the area of the contact area is approximately 100 mm2, and the load act on the contact area is about 50~200 kN, which could result in a huge contact stress [24]. Under the action of the normal contact pressure and the tangential friction force, the stress state of the rail surface can be seen from Figure 2. When the wheel rolls from position 1 to position 3, the horizontal tensile stress and shear stress will have periodic changes in direction. In general, this combined alternating stress is considered to be the main cause of the crack initiation and propagation on the rail surface.


## 2.3. Analysis of the Crack State

When the wheel is close to or away from the crack, the extrusion effect on the rail surface leads to the stretch of the surface material near the edges of the contact area. At the same time, the surface cracks of the rail will open under the action of tensile stress. While the wheel pressure acts on the crack, the crack will be closed because of the mutual extrusion of the crack surfaces, as shown in Figure3.Figure 3
State of crack on the rail surface.


## 3. Calculation Model and Methods of the Wheel/Rail Contact Fatigue

### 3.1. Calculation Model and Parameters

The type of the rail is 60 kg/m (U71Mn), whose mechanical properties are listed in Table1 [18]. The model is 176 mm in height and 600 mm in length. There is an angle of 60° between the crack and the rail surface, and the crack length varies from 0 to 2 mm. Without considering the influence of the foundation, the fixed constraint is applied at the bottom of the model. In order to describe the singularity of the crack tip, the singular element (PLANE 183) with eight nodes is adopted, and quarter-points singular degenerated elements are used to model the stress singularity around the crack tip. Then the finite element model for the whole rail and the local model at the crack tip are presented in Figure 4.Table 1
Mechanical properties of steel (U71Mn).
Type of rail Young’s modulus/GPa Poisson’s ratio Yield strength/MPa Tangent modulus/GPa U71Mn steel 210 0.3 550 21Figure 4
Finite element model of rail with crack on the surface.


### 3.2. Calculation Principle of the Stress Intensity FactorK

Based on the displacement extrapolation method, the stress intensity factor(K) at the crack tip can be calculated by the finite element software ANSYS [25–27]. The reference coordinate system is shown in Figure 5, and the relationship between the displacement of the nodes near the crack tip and the stress intensity factor is described as(7)u=KΙ4Gr2π2κ-1cos⁡θ2-cos⁡3θ2+KΙΙ4Gr2π2κ+3sin⁡θ2+sin⁡3θ2+orv=KΙ4Gr2π2κ+1sin⁡θ2-sin⁡3θ2-KΙΙ4Gr2π2κ-3cos⁡θ2+cos⁡3θ2+or,where u and v are the displacement components in the x and y directions, respectively. G is the shear modulus and K is the stress intensity factor, while κ is the coefficient related to Poisson’s ratio, and their relationship is shown in (8)κ=3-μ1+μforplanestress3-4μforplanestrain,where μ is Poisson’s ratio.Figure 5
Singular elements at the crack tip and their coordinates.When the angleθ equals ±180°, (7) can be simplified to(9)KΙ=Gvc-veκ+12πrKΙΙ=Guc-ueκ+12πr.Equation (9) can be used to evaluate the stress intensity factor K at the crack tip.


## 3.1. Calculation Model and Parameters

The type of the rail is 60 kg/m (U71Mn), whose mechanical properties are listed in Table1 [18]. The model is 176 mm in height and 600 mm in length. There is an angle of 60° between the crack and the rail surface, and the crack length varies from 0 to 2 mm. Without considering the influence of the foundation, the fixed constraint is applied at the bottom of the model. In order to describe the singularity of the crack tip, the singular element (PLANE 183) with eight nodes is adopted, and quarter-points singular degenerated elements are used to model the stress singularity around the crack tip. Then the finite element model for the whole rail and the local model at the crack tip are presented in Figure 4.Table 1
Mechanical properties of steel (U71Mn).
Type of rail Young’s modulus/GPa Poisson’s ratio Yield strength/MPa Tangent modulus/GPa U71Mn steel 210 0.3 550 21Figure 4
Finite element model of rail with crack on the surface.


## 3.2. Calculation Principle of the Stress Intensity FactorK

Based on the displacement extrapolation method, the stress intensity factor(K) at the crack tip can be calculated by the finite element software ANSYS [25–27]. The reference coordinate system is shown in Figure 5, and the relationship between the displacement of the nodes near the crack tip and the stress intensity factor is described as(7)u=KΙ4Gr2π2κ-1cos⁡θ2-cos⁡3θ2+KΙΙ4Gr2π2κ+3sin⁡θ2+sin⁡3θ2+orv=KΙ4Gr2π2κ+1sin⁡θ2-sin⁡3θ2-KΙΙ4Gr2π2κ-3cos⁡θ2+cos⁡3θ2+or,where u and v are the displacement components in the x and y directions, respectively. G is the shear modulus and K is the stress intensity factor, while κ is the coefficient related to Poisson’s ratio, and their relationship is shown in (8)κ=3-μ1+μforplanestress3-4μforplanestrain,where μ is Poisson’s ratio.Figure 5
Singular elements at the crack tip and their coordinates.When the angleθ equals ±180°, (7) can be simplified to(9)KΙ=Gvc-veκ+12πrKΙΙ=Guc-ueκ+12πr.Equation (9) can be used to evaluate the stress intensity factor K at the crack tip.


## 4. Calculation Results and Analysis

### 4.1. The Relative Dangerous Position of the Cracks

While the wheels of different axle load (10, 15, and 20 t) are rolling along the rail surface with cracks, the variation trends of the stress intensity factors (KI and KII) at the tip of cracks with different lengths (i.e., 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 1, 1.5, and 2 mm) are shown in Figure 6, where x denotes the distance from the center of the contact area to the crack.Figure 6
Position relationship between the stress intensity factors (KI and KII) and the contact area: (a) KI—10 t; (b) KII—10 t; (c) KI—15 t; (d) KII—15 t; (e) KI—20 t; (f) KII—20 t.
(a) (b) (c) (d) (e) (f)According to Figures6(a), 6(c), and 6(e), it is easy to find that, with the increase of the crack length, both the maximum value of KI and its position are changing. Under the condition that the axle load is 10 t, if the crack length is less than or equal to 0.2 mm, the inner edge of the contact area has the maximum value of KI. While the crack length is larger than 0.2 mm, there is the maximum value of KI at the outer edge of the contact area. And when the axle load is increased to 15 t, the maximum value of KI will occur at the inner edge of the contact area if the crack length is less than or equal to 0.5 mm. Similarly, when the crack length is larger than 0.5 mm, the outer edge of the contact area will have the maximum value of KI. When the axle load is 20 t, if the crack length is less than or equal to 0.7 mm, the inner edge of the contact area has the maximum value of KI. When the crack length is larger than 0.7 mm, there is the maximum value of KI at the outer edge of contact area, where the corresponding semilength values of the contact area for different axle load (10, 15, and 20 t) are 5.9, 7, and 8 mm, respectively. The results indicate that, with the increase of the crack length, the critical position where the crack has open failure is changed from the inner edge of the contact area to the outer edge.Figures6(b), 6(d), and 6(f) illustrate that, compared with KI, the maximum value of KII mainly occurs at the inner edge of the contact area. When the axle load is 10 t, the maximum value of KII varies between 4.45 and 4.75 mm. When the axle load is increased to 15 t, the maximum value of KII changes between 5.5 and 6.1 mm, while, under the axle load of 20 t, the maximum value of KII is changing between 6.4 and 6.6 mm. What is more, with the increase of the crack length, the position where the maximum value of KII occurs will be closer to the center of the contact area slightly.


### 4.2. Governing Factors of the Crack Propagation

Under the different axle load of 10 t, 15 t, and 20 t, the variation trends of the maximum stress intensity factors (KImax and KIImax) at the crack tip are shown in Figure 7.Figure 7
Variations ofKImax and KIImax with the increase of the crack length: (a) KImax and (b) KIImax.
(a) (b)As seen in Figure7, with the increase of the crack length, the value of KImax has a sharp increase at the beginning and then decreases gradually. When the crack length is 0.3 mm, the factor KImax has the maximum value. The possible explanation can be given by the distribution of the longitudinal stress (σx) on the rail surface, as shown in Figure 8(a). The maximum tensile stress occurs at the rail surface, and, with the increase of the depth, the tensile stress decreases significantly. While the depth exceeds a critical value, the rail will be in compression and the crack closure occurs. As for the factor KIImax, with the increase of the crack length, the value of KIImax also has a sharp increase first, and then its increase rate tends to be slow. The reason can be found from the shear stress distribution on the rail surface, with respect to the tensile stress caused by the wheel/rail contact, the maximum shear stress mainly occurs at a certain depth which is below the rail surface. When the axle load is 15 t, the maximum shear stress occurs at the edge of the contact area which is about 4 mm below the surface, as shown in Figure 8(b).Figure 8
Stress contour of the rail surface (the axle load is 15 t): (a)σx and (b) τxy.
(a) (b)To sum up, the rolling contact fatigue crack belongs to mixed-mode crack with open mode (I) as well as sliding mode (II). Furthermore, as for different length of the cracks, their governing factors of crack growth are not always the same. When the crack length is less than 0.3 mm, the fatigue crack is mainly open mode. While the crack length is larger than 0.3 mm, the sliding mode takes the leading role. And when the crack increases to a certain length, the crack tip tends to close, andKI is equal to zero. At the moment, the fatigue crack on the rail surface becomes pure sliding mode (II) without being open.


### 4.3. Propagation Rate of the Cracks

Since the wheel/rail contact fatigue crack belongs to mixed-mode crack (I and II), it is necessary to adopt the equivalent stress intensity factorKeff to exhibit the combined effect of the factors (KI and KII), as shown in [28](10)Keff=KΙ+KΙΙ,where KΙ and KΙΙ are stress intensity factors, which denote the open failure and sliding failure mode, respectively.The propagation rate of the crack can be obtained by the Paris-Erdogan model, as shown in [29, 30](11)dadN=CΔKm,where da/dN is the crack growth per cycle and C and m are material constants to be found experimentally, while ΔK is the stress intensity factor range, which can be determined by (12)ΔK=Kmax-Kmin,where Kmax is the maximum value of stress intensity factor K and Kmin is the minimum value of stress intensity factor K.When different axle load of the wheels (i.e. 10, 15, and 20 t) is rolling along the rail surface with cracks, the variation trends of the equivalent intensity factorsKeff and their ranges ΔKeff for different crack length are illustrated in Figure 9.Figure 9
Variations ofKeff and ΔKeff for different crack length: (a) Keff—10 t; (b) Keff—15 t; (c) Keff—20 t; and (d) ΔKeff.
(a) (b) (c) (d)According to Figures9(a), 9(b), and 9(c), it is clear that different length of the cracks all have two peak values of the equivalent stress intensity factor (Keff), and the distance between the two peaks increases with the increase of the crack length. The peak near the center of the contact area is mainly governed by the factor KΙΙ, while the peak far away from the center of the contact area is mainly controlled by the factor KΙ.From Figure9(d), it can be concluded that the equivalent stress intensity range ΔKeff has a sharp increase as the crack length increases. When the crack length is about 0.3 mm, there will be the first peak of ΔKeff occurring, and then it decreases. When the crack length is larger than 0.5 mm, the factor range increases again. That is to say, when the crack length is between 0.3 and 0.5 mm, because the rate the crack grows to a certain depth of the rail is less compared to growing to the rail surface, therefore, during this phase, the crack tends to grow towards the rail surface, and then the rail spalling occurs. This spalling path of the rail surface has a good agreement with the field observation, as shown in Figure 10 [31].Figure 10
Typical surface spalling path of the service rail.


### 4.4. Angles of the Crack Propagation

With respect to the mixed-mode crack (I and II), the propagation direction of the crack is not always the same. However, it tends to grow in the direction in which the maximum circle tensile stress occurs [32]. According to this criterion, the turning angle of the crack propagation can be obtained, as described in (13)θ=2tan-1⁡14KΙKΙΙ+14KΙKΙΙ2+8,where θ is the turning angle of the crack and its turning direction is illustrated in Figure 11.Figure 11
Turning direction of the crack based on the maximum circle tensile stress criterion.When the crack length is between 0.3 and 0.5 mm, the crack may turn to grow towards the rail surface in an opposite direction. Aimed at analyzing the reasons why the rail spalling is formed in the shape of barb, based on the maximum circle tensile stress criterion, the propagation directions of the cracks with different length (0.3, 0.4, and 0.5 mm) are predicted under different axle loading (10, 15, and 20 t). The calculation results are listed in Table2 and Figure 12.Table 2
Predicted propagation angles of the cracks based on the maximum circle tensile stress criterion.
Crack length/mm Axle load/t K e f f/MPa·m0.5 K I/MPa·m0.5 K I I/MPa·m0.5 Initial angleof the crack/° Turning angleof the crack/° Propagation angle of the crack/° 0.3 10 23.44 0.31 23.13 60 70.78 130.78 0.4 10 25.04 0.00 25.04 60 70.53 130.53 0.5 10 25.91 0.00 25.91 60 70.53 130.53 0.3 15 52.69 37.51 15.18 60 114.57 174.57 0.4 15 48.03 36.79 11.24 60 124.49 184.49 0.5 15 42.58 1.61 40.97 60 71.28 131.28 0.3 20 72.73 53.46 19.27 60 118.60 178.60 0.4 20 67.78 51.14 16.64 60 122.23 182.23 0.5 20 58.36 45.44 12.92 60 127.06 187.06Figure 12
Predicted propagation paths of the crack based on the maximum circle tensile stress criterion: (a) 10 t; (b) 15 t; and (c) 20 t.
(a) (b) (c)As is illustrated in Figure12, it can be found that when the axle load is 10 t, although the propagation direction of the crack has a certain turning, its general propagation tends downward, which is likely to fracture. While the axle load is 15 t or 20 t, the propagation direction will have a change and grow towards the rail surface. As a result, the spalling is formed in the shape of barb. In addition, according to Figures 12(b) and 12(c), compared with an axle load of 15 t, the crack is more likely to propagate towards the rail surface when the axle load is 20 t. This phenomenon may explain why the spalling of rail surface occurs more frequently in the line of heavy haul [5, 6, 31].


## 4.1. The Relative Dangerous Position of the Cracks

While the wheels of different axle load (10, 15, and 20 t) are rolling along the rail surface with cracks, the variation trends of the stress intensity factors (KI and KII) at the tip of cracks with different lengths (i.e., 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 1, 1.5, and 2 mm) are shown in Figure 6, where x denotes the distance from the center of the contact area to the crack.Figure 6
Position relationship between the stress intensity factors (KI and KII) and the contact area: (a) KI—10 t; (b) KII—10 t; (c) KI—15 t; (d) KII—15 t; (e) KI—20 t; (f) KII—20 t.
(a) (b) (c) (d) (e) (f)According to Figures6(a), 6(c), and 6(e), it is easy to find that, with the increase of the crack length, both the maximum value of KI and its position are changing. Under the condition that the axle load is 10 t, if the crack length is less than or equal to 0.2 mm, the inner edge of the contact area has the maximum value of KI. While the crack length is larger than 0.2 mm, there is the maximum value of KI at the outer edge of the contact area. And when the axle load is increased to 15 t, the maximum value of KI will occur at the inner edge of the contact area if the crack length is less than or equal to 0.5 mm. Similarly, when the crack length is larger than 0.5 mm, the outer edge of the contact area will have the maximum value of KI. When the axle load is 20 t, if the crack length is less than or equal to 0.7 mm, the inner edge of the contact area has the maximum value of KI. When the crack length is larger than 0.7 mm, there is the maximum value of KI at the outer edge of contact area, where the corresponding semilength values of the contact area for different axle load (10, 15, and 20 t) are 5.9, 7, and 8 mm, respectively. The results indicate that, with the increase of the crack length, the critical position where the crack has open failure is changed from the inner edge of the contact area to the outer edge.Figures6(b), 6(d), and 6(f) illustrate that, compared with KI, the maximum value of KII mainly occurs at the inner edge of the contact area. When the axle load is 10 t, the maximum value of KII varies between 4.45 and 4.75 mm. When the axle load is increased to 15 t, the maximum value of KII changes between 5.5 and 6.1 mm, while, under the axle load of 20 t, the maximum value of KII is changing between 6.4 and 6.6 mm. What is more, with the increase of the crack length, the position where the maximum value of KII occurs will be closer to the center of the contact area slightly.


## 4.2. Governing Factors of the Crack Propagation

Under the different axle load of 10 t, 15 t, and 20 t, the variation trends of the maximum stress intensity factors (KImax and KIImax) at the crack tip are shown in Figure 7.Figure 7
Variations ofKImax and KIImax with the increase of the crack length: (a) KImax and (b) KIImax.
(a) (b)As seen in Figure7, with the increase of the crack length, the value of KImax has a sharp increase at the beginning and then decreases gradually. When the crack length is 0.3 mm, the factor KImax has the maximum value. The possible explanation can be given by the distribution of the longitudinal stress (σx) on the rail surface, as shown in Figure 8(a). The maximum tensile stress occurs at the rail surface, and, with the increase of the depth, the tensile stress decreases significantly. While the depth exceeds a critical value, the rail will be in compression and the crack closure occurs. As for the factor KIImax, with the increase of the crack length, the value of KIImax also has a sharp increase first, and then its increase rate tends to be slow. The reason can be found from the shear stress distribution on the rail surface, with respect to the tensile stress caused by the wheel/rail contact, the maximum shear stress mainly occurs at a certain depth which is below the rail surface. When the axle load is 15 t, the maximum shear stress occurs at the edge of the contact area which is about 4 mm below the surface, as shown in Figure 8(b).Figure 8
Stress contour of the rail surface (the axle load is 15 t): (a)σx and (b) τxy.
(a) (b)To sum up, the rolling contact fatigue crack belongs to mixed-mode crack with open mode (I) as well as sliding mode (II). Furthermore, as for different length of the cracks, their governing factors of crack growth are not always the same. When the crack length is less than 0.3 mm, the fatigue crack is mainly open mode. While the crack length is larger than 0.3 mm, the sliding mode takes the leading role. And when the crack increases to a certain length, the crack tip tends to close, andKI is equal to zero. At the moment, the fatigue crack on the rail surface becomes pure sliding mode (II) without being open.


## 4.3. Propagation Rate of the Cracks

Since the wheel/rail contact fatigue crack belongs to mixed-mode crack (I and II), it is necessary to adopt the equivalent stress intensity factorKeff to exhibit the combined effect of the factors (KI and KII), as shown in [28](10)Keff=KΙ+KΙΙ,where KΙ and KΙΙ are stress intensity factors, which denote the open failure and sliding failure mode, respectively.The propagation rate of the crack can be obtained by the Paris-Erdogan model, as shown in [29, 30](11)dadN=CΔKm,where da/dN is the crack growth per cycle and C and m are material constants to be found experimentally, while ΔK is the stress intensity factor range, which can be determined by (12)ΔK=Kmax-Kmin,where Kmax is the maximum value of stress intensity factor K and Kmin is the minimum value of stress intensity factor K.When different axle load of the wheels (i.e. 10, 15, and 20 t) is rolling along the rail surface with cracks, the variation trends of the equivalent intensity factorsKeff and their ranges ΔKeff for different crack length are illustrated in Figure 9.Figure 9
Variations ofKeff and ΔKeff for different crack length: (a) Keff—10 t; (b) Keff—15 t; (c) Keff—20 t; and (d) ΔKeff.
(a) (b) (c) (d)According to Figures9(a), 9(b), and 9(c), it is clear that different length of the cracks all have two peak values of the equivalent stress intensity factor (Keff), and the distance between the two peaks increases with the increase of the crack length. The peak near the center of the contact area is mainly governed by the factor KΙΙ, while the peak far away from the center of the contact area is mainly controlled by the factor KΙ.From Figure9(d), it can be concluded that the equivalent stress intensity range ΔKeff has a sharp increase as the crack length increases. When the crack length is about 0.3 mm, there will be the first peak of ΔKeff occurring, and then it decreases. When the crack length is larger than 0.5 mm, the factor range increases again. That is to say, when the crack length is between 0.3 and 0.5 mm, because the rate the crack grows to a certain depth of the rail is less compared to growing to the rail surface, therefore, during this phase, the crack tends to grow towards the rail surface, and then the rail spalling occurs. This spalling path of the rail surface has a good agreement with the field observation, as shown in Figure 10 [31].Figure 10
Typical surface spalling path of the service rail.


## 4.4. Angles of the Crack Propagation

With respect to the mixed-mode crack (I and II), the propagation direction of the crack is not always the same. However, it tends to grow in the direction in which the maximum circle tensile stress occurs [32]. According to this criterion, the turning angle of the crack propagation can be obtained, as described in (13)θ=2tan-1⁡14KΙKΙΙ+14KΙKΙΙ2+8,where θ is the turning angle of the crack and its turning direction is illustrated in Figure 11.Figure 11
Turning direction of the crack based on the maximum circle tensile stress criterion.When the crack length is between 0.3 and 0.5 mm, the crack may turn to grow towards the rail surface in an opposite direction. Aimed at analyzing the reasons why the rail spalling is formed in the shape of barb, based on the maximum circle tensile stress criterion, the propagation directions of the cracks with different length (0.3, 0.4, and 0.5 mm) are predicted under different axle loading (10, 15, and 20 t). The calculation results are listed in Table2 and Figure 12.Table 2
Predicted propagation angles of the cracks based on the maximum circle tensile stress criterion.
Crack length/mm Axle load/t K e f f/MPa·m0.5 K I/MPa·m0.5 K I I/MPa·m0.5 Initial angleof the crack/° Turning angleof the crack/° Propagation angle of the crack/° 0.3 10 23.44 0.31 23.13 60 70.78 130.78 0.4 10 25.04 0.00 25.04 60 70.53 130.53 0.5 10 25.91 0.00 25.91 60 70.53 130.53 0.3 15 52.69 37.51 15.18 60 114.57 174.57 0.4 15 48.03 36.79 11.24 60 124.49 184.49 0.5 15 42.58 1.61 40.97 60 71.28 131.28 0.3 20 72.73 53.46 19.27 60 118.60 178.60 0.4 20 67.78 51.14 16.64 60 122.23 182.23 0.5 20 58.36 45.44 12.92 60 127.06 187.06Figure 12
Predicted propagation paths of the crack based on the maximum circle tensile stress criterion: (a) 10 t; (b) 15 t; and (c) 20 t.
(a) (b) (c)As is illustrated in Figure12, it can be found that when the axle load is 10 t, although the propagation direction of the crack has a certain turning, its general propagation tends downward, which is likely to fracture. While the axle load is 15 t or 20 t, the propagation direction will have a change and grow towards the rail surface. As a result, the spalling is formed in the shape of barb. In addition, according to Figures 12(b) and 12(c), compared with an axle load of 15 t, the crack is more likely to propagate towards the rail surface when the axle load is 20 t. This phenomenon may explain why the spalling of rail surface occurs more frequently in the line of heavy haul [5, 6, 31].


## 5. Conclusions

In this paper, in order to analyze the mechanism of the crack propagation on the rail surface, the calculation model of the wheel/rail contact fatigue was established, and the variation of the stress intensity factor at the crack tip when the crack length was increased from 0.1 mm to 2 mm was obtained. Based on the mixed-mode fracture criterion and Paris growth theory, the mechanism of the crack propagation on the rail surface was analyzed. The research results will provide a basis for improving the antifatigue performance of rail and establishing the grinding procedure. And the following conclusions can be drawn:(1)
In the process when the microcrack grows to macrocrack, the mode of the fatigue crack on the rail surface is mixed including open mode and sliding mode. And, with the increase of the crack length, the main propagation mode of the fatigue crack changes from open mode to sliding mode.(2)
With the increase of the crack length, the stress intensity factorKI increases rapidly and then decreases gradually, the relatively dangerous location of open-mode crack moves from the inner edge of the contact area to the outer edge. While the factor KII increases during the whole propagation process, and the relative dangerous location of sliding-mode crack remains the same roughly.(3)
The crack tends to propagate upward and leads to spalling defect when the crack length is between 0.3 and 0.5 mm. The propagation path is basically identical with the spalling path of the service rail in the field.



---

*Source: 7012710-2018-02-21.xml*