Applying System Dynamics of Discrete Supported Track to Analyze the Rail Corrugation Causation on Curved Urban Railway Tracks

Urban rail corrugation on curved tracks with small radii causes strong howling during operation, which has been bothering subway operating companies for many years. )erefore, revealing its causes and growth is important for the comfort and safety of subway operation. Current studies believe that the occurrence of rail corrugation is largely due to the resonant vibration of the wheel-rail system. However, little attention has been paid to the key causes of the track resonance and the practical prediction of the occurrence probability of rail corrugation on the certain track. )is paper intends to solve these above issues. Firstly, the practical model of predicting the rail corrugation growth is proposed based on the wheel-rail coupling interaction, the key causes of corrugation are investigated, and the sensitivity analysis is carried out, while the corrugation superposition model is introduced to the analyze the corrugation evolution as well as to validate the corrugation growth from the aspect of material friction and wear. Secondly, the impact of the key causes on the initiation and development of the rail corrugation is investigated based on the cosimulation. Finally, case studies validate the proposed theory model and method. )e results show that the practical prediction model for the rail corrugation growth proposed in this paper is able to estimate the occurrence possibility of rail corrugation on a specific track, and the superharmonic resonance of the track directly excited by passing vehicles eventually leads to rail corrugation. It is also found that shortwave corrugation develops more rapidly, and adjusting the support stiffness or sleeper spacing leads to fluctuations in the corrugation wavelength and its wear rate.


Introduction
Urban rail transit is an efficient and convenient mode of transportation that can transport a large number of passengers to and from remote urban areas [1,2]. e railway system is playing an increasingly important role in meeting urban travel demands [3]. e rapid development of urban transit has brought many positive consequences to the society and economy, as well as unavoidable environmental problems and potential safety hazards [4][5][6][7]. Rail corrugation is one of the most troublesome problems in urban rail transit [8]. e design of small radius curve on track will accelerate rail deterioration and increase maintenance cost [9]. In China, rail corrugation is very common in the railway lines, particularly observed in almost all the inner rails of the curved track where the radius is less than 350 m [10]. Rail corrugation exacerbates vibration and deterioration of vehicle components and causes annoying howling [11], while the failure of the vehicle components will cause great loss [12] and it has been puzzling subway operating companies for many years. Revealing its causes is useful to develop more effective maintenance to inhibit the growth of rail corrugation.
Current studies believe that the occurrence of rail corrugation is largely due to the resonant vibration of the wheel-rail system. en, how is the resonance of the wheel-rail system excited? How to evaluate the occurrence possibility of corrugation on a specific track simply and quickly? And what are the key factors affecting the formation and development of corrugation? ese are the main problems to be answered in this paper.
In this paper, the system dynamics of discrete supported track is applied to analyze the rail corrugation causation on curved urban railway tracks. e practical prediction model is proposed on the basis of the wheel-rail dynamics to estimate the rail corrugation growth as well as to investigate the key causes of the corrugation. In addition, the superimposed wear of the rail is analyzed to validate the calculated results of the corrugation growth from the aspect of material friction and wear.

Literature Review
e formation of the rail corrugation is the result of coaction of the wheel-rail interaction, rolling contact mechanics, and the wear properties of wheel-rail materials. e source of complexity often prevents finding the optimal solution [13][14][15][16][17][18]. Due to the complex and diverse causes of rail corrugation, there is currently no appropriate theoretical system that can give a reasonable explanation for all the rail corrugation observed on railway tracks [19]. In the last decades, the formation mechanism and the develop process of rail corrugation have been extensively studied. In [20,21], all the rail corrugations observed in railway transit were divided into six categories, namely, heavy-haul corrugation, light rail corrugation, track form-specific corrugation, P2 resonance corrugation, rutting corrugation, and roaring rail corrugation.
e formation of rail corrugation can be generally attributed to two mechanisms [22][23][24][25][26][27]: damage mechanism and constant frequency mechanism, while the damage mechanism is considered to be the friction and wear of the material, which is related to the wheel-rail interaction and the tread design of the rail itself [28]. In contrast, the constant frequency mechanism believed that all the categories of rail corrugation were associated with the resonant vibration of the wheel-rail system [20]. Most importantly, more and more researchers now believe that wheel-rail vibration is a critical factor in the formation and growth of rail corrugation.
In recent decades, most of the research studies aiming at the formation mechanism of rail corrugation are done from the perspective of the constant frequency mechanism.
ere are three important branches in the study of rail corrugation: the first one focuses on the rail corrugation related to the wheelset vibration, the second one concerns the rail corrugation due to the frictional self-excited vibration of wheelset or the track structure, and the third one studies the rail corrugation caused by the resonance of wheel-rail system. In the first branch, the wheelset vibration has a specific impact on rail corrugation formation. Tassilly in [24] revealed that rail corrugation was mainly related to the wheel vibration induced by rail damage, while the effect of the initial rail damage cannot be ignored. e authors in [29] held that the corrugation wavelength was affected by the first-order symmetric and antisymmetric bending vibration frequency of the front wheel. In [30], the cross excitation of wheel-axle torsional vibration and wheel-rail system vertical and lateral vibration causes corrugation on the rail surface. e influence of wheelset interactions on corrugation were discussed in [31], where the wheelset vibration played an essential role in the development of final corrugation. However, the above studies are carried out with comprehensive consideration of the line conditions of the rail corrugated section and are not applicable to all track lines, so the vibration response of the track structure cannot be ignored. In the second branch, the development of rail corrugation was eventually attributed to the self-excited vibration of wheelsets or track structure. For example, in [32], the corrugation growth was investigated and the results proved that the corrugation growth depends on the natural mode of the system. Kurzeck [33] also found that the friction self-excited vibration of the wheel-rail system was closely related to corrugation on the metro with small radius curve in Stuttgart, Germany. In [34,35], the theory that rail corrugation is caused by the friction self-excited vibration of wheelset was proposed. In [36], the generation of rail corrugation was attributed to the phenomenon of the self-excited vibration amplitude of the wheelset being greater than the length of the longitudinal axis of the wheel-rail contact spot. In [8], it considered that the main reason for corrugation appearance was a specific free vibration of the track structure. Nevertheless, selfexcited vibration-based mechanism was generally only applicable to the corrugation on the inner rail with small radius and cannot explain the corrugation occurring on the straight track. In the third branch, the resonance of the wheel-rail system contributed a lot to rail corrugation. In [37], the formation of corrugation was largely due to track resonance. In [38], wheel-rail resonance was the leading cause of corrugation. In [39], the flexibility difference of the rail and sleeper caused rail corrugation. Correa et al. used the corrugation growth as an index to analyze its development [40], and he noted that some resonance between the wheelset and track resulted in corrugation. e thesis in [41] showed that rail corrugation was related to the change in the difference of wheel-rail flexibility, discontinuous rail support, and rail support stiffness. However, the discovery of corrugation associated to wheel-rail system resonance in theory does not mean that corresponding corrugation will inevitably be formed in reality. Sometimes, the generation of corrugation is not only determined by the resonant vibration of the wheel-rail system but also limited by its growth rate.
In reality, some specific wheel-rail dynamics takes effect and leads to the development of the corresponding corrugation with the highest growth rate, while other corrugation is suppressed due to the low growth rate. In [42][43][44][45], it was found that the rail absorber could be applied to suppress the growth of rail corrugation due to the pinned-pinned resonance. In [46], the study showed that the growth of rail corrugation was slowed down with the help of the improvement of the material performance of rail. In [47], the corrugation growth was reduced due to the adjustment of the supporting stiffness of rail pad. In addition, the results in [48][49][50] showed that adjusting the friction coefficient between wheel and rail was also a technical method that could be widely used to reduce the corrugation growth. e above studies are all based on experiments or simulations. However, the research on the key causes of rail corrugation and its growth needs to provide strong evidence in theory and practice.
In summary, previous studies proposed the viewpoint that the inherent vibration of wheelsets or track structure results in corrugated wear. However, they did not analyze the key causes and how the inherent vibration was excited, while the investigation of restraining the corrugation growth lacked cogent evidences in theory. is paper concerns the causation and growth of the corrugation on curved tracks with small radii (CTSRs), and the purpose is to propose a practical rail corrugation growth model and discover its key causes in both theory and practice. Firstly, a practical growth model of corrugation based on the vibration theory of wheelrail system is proposed, from which the key causes are extracted and the sensitivity analysis is performed. Secondly, with the aid of the joint-simulation of vehicle-track dynamics and wheel/rail finite element model, the corrugation under different track conditions is calculated. On this basis, this paper proposes that the superharmonic resonance of the curved track directly excited by the passing vehicles eventually leads to corrugation. After that, the corrugation evolution, including the development of the wavelength and corrugation grade, is investigated by using the corrugation superposition model. At the same time, the calculated rail wear is used for the validation of the result of the rail corrugation growth. Finally, a case study is introduced. e main contribution of this paper lies in the following three aspects. e first one is the proposal of a practical rail corrugation growth model, which can be used to estimate the formation possibility of rail corrugation easily and quickly. e second one lies in that the sensitivity analysis of the track parameters is performed for the new corrugation growth model, where the curve radius and support stiffness are extracted as the two most important factors affecting the growth rate of corrugation. e last one is that a new view on the corrugation formation is introduced, which holds that the vehicle passing excitation caused by track discrete support structure can directly excite the track superharmonic resonance. And the superharmonic resonance eventually leads to rail corrugation. is result is completely different from the existing self-excitation theory and the other-excited theory (the natural vibration of the track structure or wheel-rail system, excited by the track irregularity, finally leads to corrugation). Here, the sleeper spacing and vehicle passing speed are the direct external causes of the resonance of track structure. It is useful for the prevention and maintenance of corrugated wear. is paper is outlined as follows: the methodology is presented in Section 3. In Section 4, the formation theory of rail corrugation is described. Section 5 focuses on the simulation validation of the corrugation evolution, and Section 6 presents the case study. Finally, the seventh section summarizes some conclusions.

Methodology
In this section, the practical prediction model of the corrugation growth based on wheel-rail dynamics and the superimposed wear model used for corrugation evolution analysis as well as the validation of corrugation growth are illustrated. After that, the simulation environment and the application of simulation data are stated. Finally, the key factors affecting the growth of rail corrugation are extracted from the derivation process of that, and the sensitivity analysis is carried out by using the cosimulation results.

Corrugation Growth Prediction Model.
e corrugation growth function was proposed by Valdivia [51] and was later improved in [24,52]. It can predict the trend of corrugation caused by the continuous development of irregularities at specific positions of tracks after thousands of trains pass by. At present, the wear growth rate has been used to predict the corrugated wear of the rail under various conditions [40,[53][54][55][56][57]. e model assumes that the wear rate is proportional to the friction force, and the variation of friction force under nominal steady-state condition is proportional to the change of creep. Assuming that the original irregularity wavelength of the line is L k , the frequency is f k , and the irregularity amplitude is Δz k . e change rate of irregularity after n trains passing by is expressed as follows [53]: where k is the wear coefficient of the material, ρ is the rail density, b is the length of the semiaxes of the wheel-rail contact patch, v is the train velocity, and ΔP frict (x, n) is the fluctuation part of the friction work. e fluctuation part of the friction work ΔP frict (x, n) can be expresses as follows [8]: where T y is the lateral creep force and C y represents the lateral creepage. en, equation (1) can be shown as follows: where λ k is the characteristic variable of equation (3), and it is a complex number that depends on the state variables of the system. e growth model of corrugation G c on position x is the real part of λ k and can be derived as follows. e detailed derivation process is explained in [8,31], and we will not explain it here. In this paper, the growth prediction function Discrete Dynamics in Nature and Society is simplified for corrugation on curved rail with small radii, and the result is validated by simulation: where N is the normal force of the wheel-rail contact, T y,0 are the actual part of the lateral creep force, C y,0 represent the real part of the lateral creepage, U y and U z are the lateral and vertical track receptance, respectively, and c is parameters related to the wheel-rail contact patch. See Table 1 for other symbols. e larger the amplitude of G c is, the more likely that corrugation is appearing on the line. It is known that positive values above 0.01 correspond to situations where rail corrugation wear certainly appears [52]. us, the frequency corresponding to the maximum amplitude is one of the main vibrations that generates the wave. e above model for the calculation of rail corrugation growth is a generalized calculation model, which is not suitable for direct estimation of rail corrugation growth in metro, and the parameters are complicated and unintuitive. erefore, the main influencing factors are highlighted to improve the above corrugated growth model. As a result, a more intuitive calculation model for the estimation of the corrugation growth is obtained. Firstly, the corrugation growth model at a specific position x can be described as follows: where K 0 is a constant as the material coefficient of rail, K 1 indicates the fluctuation of wheel-rail creep force, K 2 is the influencing of rail vertical dynamics, and K 3 is the influencing of rail lateral dynamics. And they are defined as follows: Considering Hertz's nonlinear elastic contact theory, the partial derivative from creep force T y to wheel track normal contact force N, contact spot size c and lateral creepage C y , and the vertical force between wheel and rail as well as the maximum lateral creepage can be determine as [8] zT y N � − μsgn C y C y C y,max where μ is the frictional coefficient, G is the wheel-rail contact constant, M s is the shear modulus of wheel-rail material, a, b is the length of the long and short semiaxis of the wheel-rail contact area, respectively, C 22 is the Kalker creep coefficient, and its value depends on the ratio of a, b.
As the fluctuation of creep force depends on the fluctuation of normal force, contact patch, and creepage [31], K 1 -K 3 can be further simplified as follows: where T is the ratio of creepage to maximum creepage: Dynamic calculation shows the magnitude of G − (3/2) d (1/2) U z is 10 4 , so K 2 can be further simplified as As the creep force on the small radius curve is saturated, that is, T y,0 � μN, then K 3 can be further simplified. And before that, the following problem should be solved: Note that the dynamic calculation shows the magnitude 4 Discrete Dynamics in Nature and Society en, K 3 can be reduced as follows: Finally, the corrugation growth model in equation (5) can be derived as follows: where Q Z indicates the factor of vertical vibration of track structure, Q N indicates the factor of wheel-rail normal contact, Q A indicates the factor of wheel-rail contact patch, Q Y indicates the factor of lateral vibration of track structure, and Q rz indicates the factor of wheel-rail contact deformation. e correlation coefficients in the corrugated growth model are calculated by referring to the dynamic parameters and relevant empirical parameters of wheel-rail materials (China metro rail CHN60 and wheelset UICS1002). erefore, rail corrugation growth model is finally simplified as follows: Remark 1. is paper introduces a practical prediction model of rail corrugation growth as shown in equation (15). It highlights the influence of the rail receptance (U z , U y ) and also considers the wheel-rail contact dynamics (N, a, b). In the following, the curved-track receptance (U z , U y ) is deduced and analyzed to extract the key track parameters associated to the corrugation growth on curve.

Rail Receptance Calculation Model.
e rail is simplified as a curved Euler-Bernoulli beam with equal distance and scattered point support as shown in Figure 1(a). All deformations are assumed to be so small that the linear theory is applied. Assuming that the curved beam has a constant cross section, the warpage resistance is negligible. When the wheel passes by the track at the speed v, the resonance wheel-rail force with angle frequency w f is generated due to the influence of discrete sleepers. Applying a unit resonance force e iw f t to the vertical and lateral directions of the rail, the lateral and vertical vibration equations of the curved beam are obtained as follows [58]:where u y and u z are the rail displacements in the y and z directions, respectively, φ is the torsional angle around the x-axis of a specific point at time t, I y and I z are the area moment of inertia of the y-axis and z-axis, respectively, I d is rail section torsional constant, and E is elastic modulus. And k y , k z , c y , and c z represent the stiffness and damping of the track structure in the y and z directions, respectively. See Table 1 for other symbols.
e deformation of equation (17) is derived with the help of the bending and torsion characteristics of the curved beam shown in Figure 1(b), where the torsion and vertical bending of curved beams are coupled with each other. When the beam microsegment OA has a slight bending angle θ around the y axis in the O section, the A section produces a torsion angle, and its meridian changes from DA to DA ′ , then where For small enough beam segments, Discrete Dynamics in Nature and Society z z In order to solve the differential equations in frequency domain, Fourier transform is applied with respect to time to equation (16) and equation (21) as follows [59]: where * indicates the expressions of the variables in the frequency domain, u * y and u * z are the displacement of point x in y, z direction in frequency domain, is the elastic modulus and shear modulus considering the damping of rail material, in which the η is the damping loss factor of the rail, k * y and k * z represent the lateral and vertical support complex stiffness of the fastener, and k * y � k y + ic y w, k * z � k z + ic z w, where w is the angular frequency of the track vibration.
In this paper, the generalized wave number method in [58] is introduced to solve the frequency response function of curved rail beam. According to the frequency domain mathematical mode superposition method, the dynamic response of rail beam can be obtained as follows: where where V * n (x, κ, w f ) is the mathematical mode of rail beam in frequency domain, and 2N + 1 modes are considered in this calculation.C n (κ, w f ) is the modal coordinate function vector, in which B y n (κ, w f ), B z n (κ, w f ) is the mathematical mode coordinate function of rail beam corresponding to the transverse and vertical displacement modes in frequency domain. κ is the generalized wave number [60], and κ � ((w − w f )/v).
Applying equation (24) into equation (22) and equation (23) and performing integral processing on it in the basic length (x ∈ [0, d f ]), the equations corresponding to the mth-order mode of the rail is where d f is the sleeper spacing. By sorting out equations (26) and (27), we can get where where P(κ, w f ) is the (4N + 2) × 1 order column vector with only two nonzero elements, and the j-th element P(j, 1) in the matrix can be expressed in MATLAB language as Discrete Dynamics in Nature and Society 7 is paper considers the steady-state fixed load, so v � 0 and x � x 0 , then substituting equation (28) into equation (24), and the receptance of the curved rail beams corresponding to the mth-order mode is obtained as follows: Remark 2. Equation (31) and equation (32) show that track receptance (U z , U y ) depends on the track material parameters (E, M s , I y , I z , I d ), support parameters under rail (k z , k y , c z , c y ), and the angular frequency of the external excitation (w f ). For a specified track, the track material properties are fixed, while the supporting stiffness and damping meets the design requirements of the damping ratio. e external excitation considered in this paper only refers to the vehicle passing excitation, so . erefore, the rail receptance (equation (31) and equation (32)) at a certain position x can finally be regarded as a function of the response frequency w, while its amplitude depends on the track structure variables R, d f , k z .

Sensitivity Analysis of Corrugation Growth
Model. e material parameters of the CHN60 rail of the metro track are introduced into equation (28) and equation (31) for model simplification, and the obtained rail receptance is substituted into equation (15) for the new corrugation growth model as follows: where σ � 6.28(m/d f ) and m is the modal order. It shows that the corrugation growth G c can be described as frequency f � (w/2π) dependent function, where the amplitude is limited by the contact parameters N, a, b, track parameters of R, d f , and k z . e fluctuation of variables a, b, N, U z , U y with track position x is calculated as shown in Figure 2. e discontinuity of railway track leads to significant fluctuation of the model variables, especially the rail receptance. It is found that the values of these variables fluctuate periodically along the running direction, and the period is about the length of sleeper spacing. Additionally, the track receptance fluctuates significantly with the position of the track as shown in Figure 2(b).
For the rail corrugation on curved-track with small radii, several factors affect the corrugation growth to some extent. In this paper, the dominant factors include the curve radius R, support stiffness k z , and the sleeper spacing d f . To investigate the relationship between the corrugation growth and each of these factors, the sensitivity analysis of the three track parameters is performed for the new corrugation growth model in equation (33). e contributions of these three factors to the corrugation growth are shown in Figure 3. Figure 3(a), corrugation growth G c can be described as an exciting frequency dependent function.

As presented in
is is because, according to the equation of the rail receptance U z and U y derived by us, the rail receptance is a simple harmonic wave, whose amplitude and period change with these three factors. In addition, the curve radius plays the dominant role in the growth of corrugated wear. With the increase of the curve radius, the value of G c decreases first and then decreases rapidly after the radius is greater than 600m. At the same time, the corrugation passing frequency decreases linearly. is can be explained as the large dynamic effect of the wheel-rail system on a curve with a small radius, which is easier to form rail corrugation and accelerate its development. And equation (33) also shows that G c is approximately inversely proportional to the power of the radius of curved-track. erefore, the corrugated wear rate decreases as the curvature increases. Figure 3(b) shows that the value of G c descends on the increase of the support stiffness. It indicates that the support stiffness k z has a great impact on the growth of corrugated wear as shown in equation (33). at is, G c is in direct proportional to k z but inversely proportional to the square of k y . Besides that, the corrugation frequency continues to decrease along with the increase of k z , while the rate of decrease slows down when k z 8 Discrete Dynamics in Nature and Society is larger than 40 MN/m. In comparison with the other two factors, the sleeper spacing d f has less influence on the corrugation growth but the dominant frequency of corrugated wear as shown in Figure 3(c). As the excitation frequency w f varies inversely with d f to some extent, it obtains through analysis that along with the increase of d f , the corrugation passing frequency and its growth rate decrease exponentially. However, the overall trend is gradually weakened when d f climbs to 0.65 m.

Remark 3.
In this section, the corrugation growth model is modified based on the wheel-rail dynamics. It proposes that the key factors influencing the rail corrugation growth are R, d f , and k z . With the help of sensitivity analysis, it can be concluded that R and k z contribute the most to the corrugated wear rate while k z and d f have relatively more contribution on the corrugation passing frequency (as well as the wavelength of rail corrugation). e theoretical deduction shows that the excitation of vehicles passing by is   indispensable for rail corrugation. In the following, the corrugation growth under different track conditions is calculated by using cosimulation, while the corrugation superposition model and simulation environment are introduced in the next subsection for corrugation evolution estimations and the validation of the calculated results of the corrugation growth from the aspect of material friction and wear.

Corrugation Superposition Model.
is section establishes the corrugation superposition model in view of material friction and wear, where the position of the troughs, peaks, and the main corrugated wear wavelength can be observed. According to the friction work theory in the Archard wear model [61], the rail wear mass Δm i produced is proportional to the friction work, and the wear rate in the same section is consistent, while in fact, compared with the wave crest, the rail trough wear is more serious. To solve this problem, the possible contact area near the contact point is discretized in this paper, and the wear calculation is carried out in a sufficiently small section. In the simulation, the sampling interval is 0.00007 s, the vehicle speed is calculated as 50 km/h, so the wear measurement point can be taken every 1 mm. And it is reasonable for the calculation of the medium and shortwave corrugation studied in this article. Assuming that the friction power P i (x, n) remains constant in a short time step Δt, the wear depth d i (x, n) can be expressed as where A is the wheel-rail contact area at position x. Add the wear depth and the original irregularity as the new rail irregularity to repeat the wear calculation, the accumulated wear W(x, n) after n trains passed is given as follows: where z r (n − 1) is the rail vertical irregularity after (n − 1) trains passed. e material wear coefficient k is the practical value of 2.103 × 10 13 m 3 /J, and the material density ρ of the 60 kg/m rail used in China is 7800 kg/m 3 . Alternatively, other unknown variables need to be determined by numerical simulation, including track receptance (U z , U y ) and wheelrail contact parameters (N, a, b, P i (x), A), etc. e simulation environment is stated in the next section.

Simulation Environment and Data Application.
In this section, the finite element model (FEM) [62] of the wheel-rail system is built for rail receptance calculation and vibration modal analysis. e vehicle-track coupling dynamic model is established to obtain the wheel-rail contact parameters for G c as well as the friction work variables for rail wear. e midpoint of circular curve is taken as the testing point.
is paper aims at the common short sleeper track structure, whose short sleepers are embedded in the track slab and consolidated together with no damping between them. e track bed has little effect on the high-frequency vibration ( > 50 Hz) induced by the shortwave grinding of the subway [63]. To facilitate modelling, the short sleeper is not considered separately and the lower structure is considered as a whole. e rail support stiffness and damping in the model are only the stiffness and damping of the fastener, so a relatively simple single-layer elastic point-support beam model can be used for analysis as shown in Figure 4(a). In this model, the load of the entire vehicle is evenly distributed among the four wheelsets. Furthermore, the FEM simplifies the loading method and boundary constraints to facilitate analysis.
e load added at both ends of the wheelset is applied to the axle journals at both ends of the axle through distributed coupling, and the lateral and vertical forces acting on the track are balanced by the track support force. e boundary constraints of wheelset are applied to the center of both ends of the axle, and the hinged-hinged constraints are added to the two ends of the rail. e track in FEM is long enough to ignore tail effects. e elastic fasteners and gaskets between the rail bottom and the foundation are simulated by nondimensional springs and damping elements. is model uses C3D8I hexahedral elements in ANSYS to divide the mesh, with a total of 578,752 nodes and 498,648 elements. e actual support structure under rail is modelled by changing the stiffness and damping of the spring as shown in Figure 4(b). e vehicle-track coupling dynamics model is built by using the dynamic software Universal Mechanism (UM) [64], where the B-type vehicle structure used by Guangzhou Metro is used as the modelling prototype. In this model, the Euler beam is adopted as the rail model, and the total length of the flexible track structure is 500 m. e track irregularity in this model is the measured irregularities of Guangzhou Metro Line 3, where no initial corrugation is found. In dynamic calculations, the KIK-Piotrowski contact model is used to simulate the normal parameters of wheel-rail contact [65]. e FAST-SIM algorithm is used to calculate the tangent contact parameters of the wheel-rail contact system. e friction power P i (x) and the contact area A can be directly obtained by vehicle dynamics simulation. e corrugated wear in the time domain can be calculated directly by substituting them into equation (35), from which one can observe the corrugated wear depth. By using the MATLAB toolbox, the power spectrum corresponding to the one-third octave band is calculated, where the evolution tendency of the corrugation grade (CG) with the cumulative number of vehicles passing by (CVP) and the passing frequency can be extracted.

Causation Analysis of Rail Corrugation
is section calculates the corrugation growth based on the cosimulation, and the relation between the corrugation passing frequency and the track natural vibration frequency is discussed.

e Corrugation Growth under Various Track Parameters.
With the help of ANSYS, the rail receptance (U y , U z ) is calculated for different track parameters. Due to space limitations, only the receptance of curved-track with different radii is shown in Figure 5.
ere is a peak response for both U y and U z for different track radii. e track with a small radius has a much smaller response, and the corresponding frequency of the peak is larger. Additionally, with the increase in the support stiffness, the vertical and lateral track receptance tend to decrease. Meanwhile, the corresponding frequencies move toward lower. Furthermore, the impact of the support spacing on the natural vibration cannot be ignored, especially for vibrations in the frequency range of 10 Hz-400 Hz. Dense sleepers lead to high-frequency short wave vibration with small amplitude, while large support spacing aggravates the low-frequency vibration of the track. e calculated results are in good agreement with those conclusions in [47,66]. e corrugation growth function is obtained as shown in Figure 6. It shows that the rail corrugation growth on a curved track with smaller radii achieves a larger value, which means that corrugation is more likely to occur on a small radius curve, and it is consistent with the frequent locations of the rail corrugation on-site. In fact, the maximum value of the corrugation growth G c for a track with a radius of 200 m is 3.83×10 − 3 , while the value for a track with a radius of 1000 m is 0.87×10 − 6 , which is much smaller. e primary frequency of the corrugation on the curved track with a radius of 200 m is 202 Hz. e primary frequency of the curved rail with radii of 300 m, 400 m, and 1000 m is 184 Hz, 168 Hz, and 78 Hz, respectively. As the running speed of the train is 50 km/h, the wavelengths of rail corrugation in the above four curves are obtained as follows: e corrugation growth rate G c variation to the support stiffness is shown in Figure 6(b). It shows that the track with a smaller support stiffness is more likely to form corrugation with a lower frequency and a longer wavelength, which verifies the conclusion in [47]. For the curved track with a support stiffness of k z1 , the primary frequency is 172 Hz, and the wavelength is 80 mm. When the support stiffness of the track increases to k z4 , the passing frequency is approximately 80 Hz, and the corrugation wavelength is only 173 mm. Figure 6(c) shows that reducing the sleeper spacing effectively suppresses the corrugated wear. More importantly, the sleeper spacing seriously affects the corrugation wavelength. e figure shows that for a track structure with a support spacing of 0.5 m, the corrugation frequency is 202 Hz. When the sleeper spacing increases to 0.8 m, the primary frequency of corrugation decreases to 150 Hz. e above results are consistent with the sensitivity analysis of the previous corrugation growth G c and also consistent with the experimental results in [66].

e Relation between Corrugation Passing Frequency and
Track Natural Vibration Frequency. To investigate the track vibration related to rail corrugation, the modal analysis of the curved track is carried out. Taking the modal analysis of curved tracks with different radius as an example due to space limitations, Figure 7 shows the typical vibration modal of the curved track. At 196 Hz, the fifth-order vertical bending vibration coupled with transverse torsional vibration occurs. As the natural frequency of 196 Hz is close to the resonance frequency (202 Hz) of the corrugation growth function G c in Figure 6(a), it reveals that the corrugation is directly related to this fifth-order vertical bending vibration.
Similarly, the 4th-order bending torsional vibration of 190 Hz in the curved track with a radius of 300 m results in corrugation with a wavelength of 75 mm. e vertical 4thorder bending vibration of 148 Hz in the curved track with a 400 m radius determines rail corrugation with the main wavelength of 94 mm. e long corrugation (wavelength 178 mm) on the inner rail of a 1000 m radius curve is caused by the first-order lateral bending vibration of the curved track at 82 Hz. erefore, how is the natural vibration of the track excited? It is certainly excited by the passing vehicles. For a curved track with sleeper spacing of 0.6 m, the excitation frequency caused by a vehicle passing at a speed of 50 km/h is 23 Hz f � (v/d f ). e natural frequency of the track related to corrugation is close to a certain frequency doubling of the vehicle passing excitation. For example, the natural frequencies of the above four curved tracks with different radii are 196 Hz, 190 Hz, 148 Hz, and 82 Hz, which are approximately nine times, eight times, six times, and four times the driving frequency, respectively.
For railway tracks with different sleeper spacings, the excitation frequency due to the vehicle passing by is calculated as shown in Table 2.
e results show that the primary frequency of rail corrugation in the curves with different sleeper spacings is quite different. However, it is very close to the natural vibration frequency of the track and a high-order harmonic frequency of the excitation frequency.

Remark 4.
e higher harmonic resonances of the track are excited when the vehicle passes by due to its discontinuous support characteristics, and it causes the initial corrugated wear on the rail surface. ese results are completely different from the existing self-excitation theory and the otherexcited theory. And a new corrugation formation theory is innovatively presented. It is the vehicle passing speed and the discontinuity of track structure that directly excite the superharmonic resonance of the track structure, which eventually leads to rail corrugation. e wavelength of the corrugation depends on the high-order harmonic resonance frequency and the vehicle passing speed.

Simulation of the Corrugation Evolution on a Curved Track Surface
In this section, the accumulated rail wear under different scenarios is calculated. e results are found to be in good agreement with the calculation results of corrugation growth. However, due to space limitations, this section only introduces the results of the superimposed wear on a curved inner rail with different radii for validation by simulation. Figure 8(a) shows the rail accumulated wear after 10 thousand vehicles pass by, where the curved rail with a small radius has a large wear on the surface. For example, the wear depth of the curved rail with a radius of 200 m is as deep as 0.117 mm, while that of the curve with a radius of 1000 m is only 0.00632 mm.
at is, the corrugation on the small radius curve evolves much faster than that on curves with larger radii, which is consistent with the theoretical results of the previous sensitivity analysis. e primary wavelength of 71 mm is close to the calculation results λ R200 � 69 mm of the corrugation growth model, as shown in Figure 6(a). Similar observations can be found for the other three rails.
As we calculated, with the decrease in the support stiffness of the track, the rail wear grows faster, and the shortwave wear becomes more prominent. Besides that, rail wear in different curves presents two characteristic frequencies. e larger frequency is the primary frequency of the corrugation, which decreases faster with the increase in the support stiffness. Furthermore, different sleeper spacing results in different wavelengths of corrugation, and shortwave wear expands faster on the tracks with dense supports. With the increase in sleeper spacing, the corrugation wavelength gradually increases. Similarly, the above wear frequencies are very close to the frequencies calculated in Figure 6(b) and 6(c).
To investigate the evolution of the wavelength and corrugation grade as the CVP steadily increases, the development trend of the CG at the primary passing frequency in different curves is calculated. And only the corrugation evolution with different radii is shown in Figure 9 due to space limitations.
It shows that the CG increases almost linearly with the increase of the CVP, where corrugation on a smaller curve develops faster. After a total of 20,000 trains pass by the track line, the corrugation level at the main frequency of 100 Hz is -42 dB for the curve with a radius of 200 m. It is only -62 dB in the curve with a radius of 1000 m. When the total traffic   corrugation, and its development is accelerated by the decrease in support stiffness. Moreover, different sleeper spacings have little impact on the development of that. erefore, the adjustment of the support stiffness or sleeper spacing affects the corrugation wavelength, and the development of corrugations is significantly affected by the support stiffness. In the next section, the corrugation measured from Guangzhou Metro is introduced to verify these studies.

Validation of the Rail Corrugation Causation and Evolution with Field Measurements
According to the quarterly monitoring data on-site, the internal track of the downward line from Shiqiao Station to Panyu Square Station of Guangzhou Metro Line 3 is a section with frequent corrugation. e measured track section selected in this paper is common short sleeper with monolithic track bed, and the curve radius is 300 m with the  Discrete Dynamics in Nature and Society superelevation of 100 mm. e medium and long waves are dominant, as shown in Figure 10. Figure 11(a) shows the measured rail corrugation obtained from four on-site tests with the help of the CAT corrugation inspection vehicle. During these four tests, no grinding treatment is carried out. e corrugation grows gradually with the increase of the CVP. e maximum corrugation depth evolves from 0.084 mm at the first inspection to 0.157 mm at the fourth inspection. e periodicity of the corrugation becomes increasingly significant from the first inspection to the fourth. Figure 11(b) shows the frequency spectrum of the rail corrugation within 600 Hz, where the primary and secondary frequency is 80 Hz and 160 Hz, separately. With the increase in CVP, these two frequencies of the corrugation hardly change. e running speed of vehicles in this section is approximately 60 km/h, and the corresponding wavelengths of the corrugation are calculated as 208 mm and 104 mm. e shortwave corrugation develops more rapidly with the increase in CVP. is is in good agreement with the previous simulation results.
e vibrating modal analysis of the measured track structure is shown in Figure 12.
e track structure is dominated by vertical second-order bending vibration with slight torsional vibration at 75 Hz, and the coupling of vertical high-order bending and severe torsion is observed at 160 Hz. e frequencies coincide with the corrugation dominant frequencies, which verifies the conclusion obtained from the above simulation that the generation of rail corrugation is related to the natural mode of the track.
How are these two natural modes excited? e sleeper spacing is 0.625 m, and the vehicle passing frequency is calculated to be 26 Hz. Additionally, the triple frequency (78 Hz) and sextuple-frequency (156 Hz) are close to the natural vibration frequencies of 75 Hz and 160 Hz, respectively. erefore, the vehicle passing incentive excited by track discontinuity support excites the high-order harmonic resonance of the track structure, which leads to rail corrugation. e theory proposed in this paper that corrugation is caused by the high-order harmonic resonance of track structures is proven.

Conclusions
is paper concerns the causation and growth of the rail corrugation on curved tracks with small radii, where a practical prediction model of rail corrugation growth is proposed. And the following conclusions are obtained. (1) e superharmonic resonance of the track is the key cause of rail corrugation, and this superharmonic resonance is directly excited by the passing excitation induced during the train negotiating the discrete support track. (2) e curve radius, supporting stiffness of fastener, and sleeper spacing are the three critical factors affecting the growth of corrugation. (3) e adjustment of the support stiffness or sleeper spacing leads to fluctuations in the corrugation wavelength and its growth rate, while reducing the support stiffness and the sleeper spacing can suppress the formation of rail corrugation. e above studies coincide well with the conclusions obtained in existing studies and are of great significance to the prevention and maintenance of corrugated wear, while the real-time detection of rail corrugation based on deep learning network [67,68] will be investigated in the future research.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest with respect to the research, authorship, and the publication of this article.