Abstract

By analyzing the measured data of rail corrugation in a small radius curve track with Cologne fasteners, the typical passing frequencies of corrugation were obtained. Then, according to actual line conditions, the vehicle-track coupled model and rail material friction and wear model were established and simulated, and rail corrugation characteristics of the curved section were studied. The results show that the variation amplitude of creep force on the inside of a guiding wheelset is large, and it coincides with that of the saturated creep force partially and the coincident part appears periodically. The variation amplitude of creep force on the outside of the guiding wheelset is small, and it is approximately equal to that of the saturated creep force. In addition, the variation range of the wear depth curve on inner rail is mainly changing from 0 μm to 0.0014 μm, showing periodic wavy wear, and when the creep force is equal to the saturated creep force, the wear depth reaches the peak. The variation range of the wear depth curve on outer rail is mainly changing from 0.0005 μm to 0.0008 μm, showing uniform wear. Combined with the predicted profile of rail wear, it can be seen that the inner rail mainly suffers from corrugation on the top of rail and the outer rail mainly suffers from uniform wear on the side, which results in serious corrugation on inner rail and slight corrugation on outer rail. The frequency characteristics analysis of wear shows that the inner rail wear has characteristic frequencies similar to the passing frequencies. Modal analysis results show that the vibration of inner rail at the characteristic frequencies is greater than that of outer rail, which lead to the corrugation of the corresponding frequencies more easily. The wear growth rates at the characteristic frequencies are relatively large, which indicates that the wear at the corresponding frequencies will continue to develop and eventually form corrugation.

1. Introduction

Corrugation is a kind of defect on rail heads, which is characterized by periodic irregularity of rail surface, and its shape is similar to that of wave. Corrugation occurs to varying degrees in metro systems all over the world, especially in small radius curve sections [14]. The wavelength of corrugation is generally between 20 mm and 300 mm, and the wave depth is not more than 1 mm [5]. Generally speaking, if the recurrence period of rail corrugation in a small radius curve section is not less than half a year, it can be regarded as normal corrugation. If rail corrugation occurs on a large radius curve section or straight line and the recurrence period is less than four months, it is considered as abnormal corrugation [6]. Rail corrugation has a serious negative impact on vehicles, tracks, and surrounding environment.

Since rail corrugation was first discovered in 1895, people have observed and investigated corrugation for more than 100 years. During this period, a lot of results on theoretical models, indoor tests, and field tests of corrugation had been presented. Table 1 lists a part of researchers and their corresponding research results.

The existing results on corrugation are mainly based on time domain (nonlinear) and frequency domain (linear). Based on the nonlinear vehicle-track interaction, the nonlinear factors were introduced into the calculation model of rail corrugation for the first time by Igeland and Ilias [12]. Nielsen [23] used the nonlinear rail wear model to study the generation of corrugation on straight lines, and validated the model with the field test on Dutch railways. Sun et al. [24, 25] proposed a simplified feedback model of corrugation including the delay of wheel transfer time. Stability analysis was used to determine the growth of corrugation when multiple wheelsets passed. Sheng et al. [26] used an unsteady two-dimensional method and a three-dimensional FASTSIM model to calculate the wheel-rail tangential creep force and predicted the development trend of corrugation through repeated loading on wheelsets. Chen et al. [27] founded the finite element time domain model of the wheelset passing sharp radius curve. The simulation results showed that the wheel-rail self-induced vibration is an important cause of curve corrugation. Tassilly and Vincent [28] studied the transfer function among the wheel, rail initial roughness, and wear rate by using the linear model of corrugation. It was illustrated that the increase of rail corrugation in some frequency bands was mainly due to the combined effect of creepage between wheels and rails, wheel-rail parameters, and initial irregularity of rail surface. Based on the previous linear theoretical model, Müller [29] assumed that the wheel is in a static state, and the rail surface irregularity was loaded between wheels and rails to simulate the dynamic characteristics of wheels passing track. Li [30] used the finite element model to analyze the relationship between passing frequencies of corrugation and natural frequencies of track structure. However, subject to the complexity of interaction between wheels and rails and the influence of various factors, it is still an open problem on building convincing explanations for the formation and evolution mechanisms of corrugation.

Based on the measured data of corrugation obtained from our experiments, a serious corrugation phenomenon appears on the curve section of the Cologne Egg fastener track, as an unusually severe one in comparison with the ones on the straight line. More generally, the inner rail of the curve track with a radius less than 350 m produces corrugation of almost 100%, while outer rail seldom produces corrugation or has slight corrugation. In this paper, the corrugation characteristics of inner and outer rails in the small radius curve section of the Cologne fastener track are investigated by both computational simulations and analysis of the measured data. Then, the vehicle-track coupled model and the rail material friction and wear model are established, and the wheel-rail creep characteristics are analyzed in time domain. Combined with the prediction profile of rail wear, the phenomena of serious corrugation on inner rail and slight corrugation on outer rail in the small radius curve section are explained. Finally, by using the frequency domain method, the development characteristics of wear are analyzed by using wear growth rate and the wear growth rates at the characteristic frequencies to develop and eventually form corrugation are discussed.

2. Characteristics of Measured Corrugation

The operating vehicle in the corrugation section is a metro type A vehicle. The curve radius is 350 m and foundation under rails is the integral track bed structure. The irregularity of measured corrugation (inner rail) and one-third octave wavelength spectrum of irregularity grades are shown in Figures 1(a) and 1(b). The characteristic wavelengths shown in Figure 1(b) are 20 mm, 31.5 mm, 63∼80 mm. The running speed of the vehicle is about 40 km/h, and the passing frequency corresponding to the specific wavelength can be obtained according to the following equation:where is the speed. The passing frequencies are

Further, the field measurement is carried out as shown in Figure 2; the wavelength ranges of corrugation are about 20∼35 mm and 55∼85 mm, which are close to that listed in Figure 1(b), indicating that the numerical calculation result is correct.

3. Vehicle-Track Space Coupled Model

3.1. Vehicle Model

The main parameters of a metro type A vehicle are shown in Table 2. The software UM (Universal Mechanism) is used to build the vehicle model of a multibody system, referring to the literature [31] for details.

3.2. Track Model

The track model is a flexible track model and parameter values of the track structure are shown in Table 3. The detailed establishment process of the track model is shown in the literature [31].

3.3. Wheel-Rail Contact Model

The wheel-rail contact model adopted the Kik–Piotrowski contact model, which has the characteristics of fast calculation speed and high reliability [32]. The vehicle-track space coupled model diagram is shown in Figure 3.

4. Wear Calculation Model

The friction and wear calculation model based on friction work theory is used to calculate rail surface wear. The model can consider the effect of spin creepage on wear. Assuming that the friction power remains unchanged within a shorter time step , the equation for calculating friction work is as follows:where is the friction power, which is the sum of the product of creep forces and creepages in all directions, as shown in the following equation:

In the UM software, the simulation can output the friction power in the contact area at any time, calculate the friction work in the corresponding time step, and then calculate the total wear amount at that time.

Considering Archard’s wear model, the mathematic model of the rail material wear can be expressed as follows:where is the total wear amount produced at position and stands for the wear ratio coefficient, which is 1 × 10−9 kg/Nm [30].

After the total wear amount is obtained, the wear depth of rail surface corresponding to the position at that time can be calculated according to the following equation:where represents the area of contact zone and is the density of rail material, taking 7800 kg/m3.

In the abovementioned vehicle-track model, the four wheels are running on each rail at the same time, and the wear depth corresponding to each wheel is accumulated at the fixed position of the rail so that the wear depth under single operation can be expressed as follows:where is the wear depth of each wheel (see Figure 3) at the fixed rail position. By adding the wear depth calculated by the single operation as irregularity to rail surface, a new rail profile can be obtained after the single operation:where is the initial irregularity of rail surface. In what follows, the next wear amount under the influence of previous wear can be obtained by calculating again. Repeating the process, the superimposed wear on rail surface can be obtained for many times. The flow chart of wear superposition calculation is shown in Figure 4.

By using the Fourier transform of rail profile after times operation, the corresponding rail profile amplitude in frequency domain can be obtained. Compared with the initial irregularity, the wear growth rate at frequency is defined as follows [33, 34]:where is the profile amplitude of initial irregularity of rail surface after the Fourier transform, is the sensitivity of wear, is the sensitivity of transverse frictional power induced wear variations to vertical wheel-rail contact deflection variations, is the modal sensitivity of the wheel-rail displacement to a change in longitudinal rail profile, is the number of wheel passing, is the modal damping ratio, is the modal natural frequency, is the angular frequency, for the frequency domain, is the ratio of the current speed to the mean speed, and is the probability distribution of . According to equation (9), the wear growth rate is a function of frequency . If wear growth rate is positive, it means that rail surface irregularity will continue to develop at the corresponding frequencies. The larger the numerical values, the faster the development, and finally the corrugation will be formed.

5. Rail Corrugation Characteristics in Curve Section

5.1. Creep Characteristics

The radius of research curve line is 350 m and the vehicle speed is 40 km/h. In order to investigate the wheel-rail creep characteristics, the software UM is used to transform the data of measured corrugation into the irregularity spectrum, and then the generated irregularity spectrum is applied to left and right rail surfaces. The creep force and saturated creep force between guiding wheelsets and rails and corresponding rail wear are simulated and calculated. The data curves are shown in Figures 5 and 6.

The results in Figures 5 and 6 can be summarized as follows:(1)The variation amplitude of creep force on the inside of the guiding wheelset is large, which is partially coincident with that of the saturated creep force and the coincident part appears periodically. The variation amplitude of creep force on the outside of the guiding wheelset is small, which is approximately equal to that of saturated creep force.(2)The variation period of the wear depth curve on inner rail is consistent with that of the creep force curve (there is a certain phase difference), and the variation range changes from 0 μm to 0.0014 μm and shows a periodic wavy wear. When the creep force is equal to the saturated creep force, the wear depth reaches the peak. The variation range of the wear depth curve on outer rail is small, which mainly changes from 0.0005 μm to 0.0008 μm, and implies a kind of uniform wear. Therefore, the inner rail tends to have corrugation while the outer rail tends to have uniform corrugation, which is consistent with the fact that corrugation mostly occurs on inner rail of the small radius curve section.

5.2. Prediction Profile of Rail Wear

In order to further explain the serious corrugation on the inner rail and slight corrugation on the outer rail, the wear position of the rail profile is predicted and analyzed. According to the literature [35], the maximum rail wear depth 0.1 mm is selected as the basis for profile renewal, and the rail profile is smoothed.

Taking the update times as 10, the profile contrast diagram and cumulative wear contrast diagram of the inner and outer rails are obtained as shown in Figures 7 and 8.

The analysis results show that the wear of outer rail distributes in the middle of the rail head and shoulder area (−10 mm, 35 mm), and the wear of the shoulder area (20 mm, 30 mm) is large, with the maximum value reaching 0.40 mm (rail transverse position 23.1 mm). The wear distribution of the inner rail head tread is symmetrical, and the wear of the middle part of the rail head (−10 mm, 10 mm) is relatively large, with the maximum value reaching 0.35 mm (rail transverse position −1.1 mm).

In order to analyze the reasons for the difference of wears between inner and outer rails, the wheel-rail geometric contact relationship in the simulation process is studied. Since the geometric parameters involved in the wheel-rail space geometric contact relationship are in the principle functions of the wheelset transverse displacement and contact angle (angle between wheel-rail contact surface and horizontal surface), several geometric parameters can be simplified to functions of wheelset transverse displacement and contact angle [36] when studying the difference of wear distribution between inner and outer rails. Figure 9 shows the change of wheelset transverse displacement during the operation of vehicle, and Figure 10 shows the change of the corresponding contact angle. When the vehicle is passing on the small radius curve section, the wheelset traverse displacement is large, and it is in the range of 9∼10 mm, and the contact angle changes slightly, approaching zero.

Based on the matching relationship between the LMA wheel tread and CN60 rail profile, the transient contact solutions of the wheelset passing curve with different transverse displacements and contact angles can be obtained, as shown in Figure 11. The analysis shows that when the wheelset transverse displacement reaches 9∼10 mm and the contact angle approaches zero, the wheel flange attaches to the rail. Increasing the wheelset transverse displacement will cause the wheel flange to move toward outer rail, resulting in a large outer wheel-rail contact angle, which is prone to rail side wear. The inner wheel-rail contact angle is smaller, and wear mainly occurs on the middle part of the rail head. This can also be derived from the wheel-rail contact diagram shown in Figure 12.

Furthermore, combined with the wheel-rail contact patch diagram shown in Figure 13, it can be obtained that the outer wheel-rail contact patch is mainly located in the shoulder area, resulting in large wear in the corresponding area, which is consistent with the wear distribution of outer rail shown in Figure 7. The inner wheel-rail contact patch mainly distributes in the middle part of the rail head, making the wear of the corresponding area large, which is also consistent with the wear distribution of inner rail shown in Figure 8. Therefore, from the angle of wheel-rail geometric contact relationship, the difference of wear distribution between inner and outer rails in a small radius curve can be well explained.

According to the predicted profile of rail wear, the outer rail wear mainly occurs on the rail shoulder, while the inner rail wear mainly occurs on the middle of the rail head. Combining the conclusion in Section 5.1 (the inner rail tends to have corrugation while the outer rail tends to have uniform wear), the inner rail mainly has corrugation on the top of rail, and the outer rail mainly has uniform wear on the side, which results in serious corrugation on inner rail and slight corrugation on outer rail.

5.3. Frequency-Domain Characteristics of Wear Growth

According to the data of measured corrugation in Section 2, the maximum passing frequency is about 556 Hz, so the frequency range in frequency domain analysis is set to 0∼600 Hz. The circular curve (R = 350 m) with a length of 200 m is selected for simulation, the measured rail surface irregularity is taken as the initial irregularity to calculate the wear amount (the additional wear generated by train running, considering 4 wheel runs on a rail) under a single operation, and the characteristic frequencies of rail corrugation (inner rail) can be obtained by the Fourier transform of wear results, as shown in Figure 14.

As can be seen from Figure 14, the basic frequency of characteristic frequencies of rail wear is 35.1 Hz, which is caused by the wheel-rail low frequency impact force P2. When the wheel passes the low joint, weld, corrugation section, or wheel scratch, there will be two impact forces: high-frequency impact force P1 and low-frequency impact force P2. Force P1 is the impact force caused by the prominent surface defects such as joints and wheel flat. The P1 resonance can be regarded as the vibration caused by the wheel-rail contact deformation between unsprung mass and rail, and its frequency range is 500 ∼ 1000 Hz. It has a very short acting time and usually causes damage to the surfaces of wheels and rails. Force P2 is the impact force caused by the inevitable track irregularity and wheel circumferential irregularity. The P2 resonance can be regarded as the vibration of unsprung mass and track as the whole mass on the resilient foundation under rails, and its frequency range is 20∼100 Hz. It has a long acting time, which can be transmitted to the vehicle parts upward and to the slab downward [37, 38]. Besides, it can be seen from Figure 14 that all other arising frequencies are harmonics of the frequency of 35 Hz, which is considered to be caused by the main frequency separation phenomenon of wheel-rail system responses under the wheel rotation effect [39, 40]. Meanwhile, the inner rail wear contains characteristic frequencies (142 Hz, 177 Hz, 355 Hz, and 560 Hz) similar to the passing frequencies of measured corrugation (139∼176 Hz, 353 Hz, and 556 Hz).

Due to fixed frequency characteristics of corrugation, the vibration characteristics of the wheelset-track structure are generally considered as the main cause of corrugation occurrence at the specific frequencies. Using the finite element software ABAQUS, the track structure model with fastener spacing 0.6 m, track centerline length 36 m, curve radius 350 m, and wheelset model with a diameter 840 mm is established, as shown in Figure 15. The wheel tread is the LMA tread and cross-section size of the slab is 2700 mm × 300 mm.

Through the modal analysis of the abovementioned model, the natural frequencies of the wheelset-track system and corresponding modes can be obtained, and the modes close to the passing frequencies are shown in Figure 16.

It can be seen from the analysis that the vibration of the inner rail is more violent than that of the outer rail at characteristic frequencies of the wheel-track system, which makes the inner rail more likely to produce the corrugation of corresponding frequencies. This is consistent with the conclusion in Section 5.2 that the inner rail corrugation is serious and the outer rail corrugation is slight.

In order to measure the development level of wear at different frequencies, the wear growth rates at different frequencies of the wear amount (inner rail) relative to the initial irregularity are calculated according to equation (9), as shown in Figure 17.

It can be seen from Figure 17 that the wear growth rates at the characteristic frequencies of 142 Hz, 177 Hz, 355 Hz, and 560 Hz are relatively high, and these frequencies are similar to the passing frequencies. This indicates that the wear of corresponding frequencies will continue to develop and finally form corrugation.

6. Conclusions

Based on the measured data in the practical work conditions, the vehicle-track model and the rail material friction and the wear model are built. Further, corrugation characteristics of the Cologne Egg fastener track in the small radius curve are investigated. The following conclusions can be drawn:(1)The variation amplitude of creep force on the inside of the guiding wheelset is large, which is partially coincident with that of the saturated creep force and the coincident part appears periodically. The variation amplitude of creep force on the outside of the guiding wheelset is small, which is approximately equal to that of the saturated creep force. In addition, the variation period of the wear depth curve on inner rail is consistent with that of the creep force curve (there is a certain phase difference), the variation range is mainly changing from 0 μm to 0.0014 μm, showing periodic wavy wear, and when the creep force is equal to the saturated creep force, the wear depth reaches the peak. The variation range of the wear depth curve on outer rail is small, mainly changing from 0.0005 μm to 0.0008 μm, showing uniform wear.(2)Combining with the prediction profile of rail wear, it can be seen that side wear mainly occurs on the outer rail, while the corrugation tends to occur on inner rail, which results in serious corrugation on inner rail and slight corrugation on outer rail.(3)The inner rail wear in the small radius curve contains the characteristic frequencies (142 Hz, 177 Hz, 355 Hz, and 560 Hz) similar to the passing frequencies (139∼176 Hz, 353 Hz, and 556 Hz). The simulated results of modal analysis show that the vibration of inner rail is more violent than that of outer rail at the characteristic frequencies of the wheel-track system, which makes inner rail more likely to produce corrugation of corresponding frequencies. The wear growth rates at characteristic frequencies are relatively high, which indicates that the wear at corresponding frequencies will continue to develop and eventually form corrugation.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 11772230).