Polygonal Wear Mechanism of High-Speed Train Wheels Based on Lateral Friction Self-Excited Vibration

: This work details research on the formation mechanism of wheel polygonalization in high-speed trains and its effect factors by numerical modeling in order to prevent the increasingly prevalent problem of wheel polygonal wear. The lateral self-excited vibration model of a wheel was developed using the LuGre friction model and self-excited vibration theory. The properties of wheel self-excited vibration and the crucial condition of Hopf bifurcation were investigated; the process of wheel polygonal wear was simulated and the results were validated using ﬁeld tracking data. The results demonstrated that periodic self-excited vibration generated by Hopf bifurcation is a required condition for polygonal wheel attrition at a given speed. The wheel’s polygonal wear has the following characteristics: “Constant speed—Self-excited—Fixed frequency—Divisible.” The order of the polygon is determined by the ratio of the wheel lateral self-excited vibration frequency to its rotational frequency. Wheel polygonal wear was caused by the vertical dynamic force of the wheel rail. The ﬁndings of the study can serve as a theoretical foundation for the prediction and reduction of wheel polygonal wear.


Introduction
Wheel polygonalization, particularly high-order polygonalization, is becoming increasingly popular. This not only causes premature wheel scrapping, but also damages or causes failure of vehicle and track components, jeopardizing operational safety [1,2]. Wang, et al. [3] undertook field tracking experiments of high-speed EMUs from various railway lines, operating speeds, and train models, and found that wheelbase, fastener type, and ambient temperature are all strongly connected to polygonal wear. Tao, et al. [4] studied the polygonal wear characteristics of Chinese mainstream electric locomotive wheels. They discovered the process of polygonal wear generation and major influencing elements, and the influence of wheel re-profiling on the wheel polygon was also researched and addressed. Wu, et al. [5] conducted extensive field tests to examine the process of highorder polygonal wear on high-speed wheels. It was discovered that altering the running speed significantly prevented polygonal wear. Nielsen and Johansson [6] investigated the categorization, fundamental causes, and effects of abnormalities in long-wavelength wheel tread. Johansson, et al. [7] modelled the evolution of wheel polygonal wear and discovered that vertical track anti-resonance and P2resonance at 165 Hz were the sources of polygonal wear. Morys [8] discovered that the third-order initial polygonal wear arouses the first-order bending vibration of the wheelset, generating lateral creep and increasing the development of polygonal wear.
Chen, et al. [9] postulated that self-excited friction vibration of a wheel-rail system generates polygonal wear in 2011, and projected comparable consequences to that of a highspeed train. Zhao, et al. [10] investigated the features of wheel-rail self-excited vibration, as well as the mechanism and affecting variables of polygonal wheel wear. They claimed that polygonation was easily induced by friction self-excited vibration caused by saturation of wheel-rail creep force on a straight path. Cui, et al. [11] studied the dynamic cause of rail corrugation by the theories of self-excitation. Wu, et al. [12] investigated the effect of disc brake units on the polygonal wear generation of high-speed railway wheels by friction self-excited vibration. However, whether the condition of self-excited friction vibration can be accomplished in operation, and whether the wear would appear polygon after the self-excited friction vibration, has yet to be examined and validated.
The circumstances of polygonal wear of wheels were examined in this research using the wear work model, and the calculation equation of polygonal wear order of wheels was provided. The stability of wheel self-excited vibration was investigated using a dynamic model of wheel lateral self-excited vibration based on the LuGre friction model. The properties of wheel polygonal wear were discovered, and field data were validated. This establishes a theoretical foundation for comprehending the mechanics of polygonal wear.

Model of Wheel Circumference Wear Depth
Archard and wear work models of wheel material wear are examples. The study found that the wheel polygon development tendency of the two models is almost identical; however, the Archard model has a higher wear coefficient. The results of the wear work model are more accurate, and the computation efficiency is greater [13,14]. As a result, the wear work model is used for simulation analysis in this research. The wear mass ∆m is proportional to the wear work W w , according to the following model: where, K w and P w are the wear work coefficient and wear power, respectively, and τ, s, and A are the lateral stress in the contact spot element, the creep speed, the contact spot area, respectively. Assuming θ is a point of the circumference of wheel rolling, the radius is R(θ) and the radius of the last turn of the wheel is R(θ − 2π). The wear depth ∆r = R(θ − 2π) − R(θ). In addition, the wear mass is the product of the wear depth, contact spot area, and material density, that is, ∆m = ρA∆r.
The wear depth results of the global technique and the local method for the contact model concur well. However, the latter's calculation speed is double that of the former [15]. As a result, this research employs a global technique to investigate the contact parameters on the whole wheel-rail contact region, focusing on the average value of normal force, tangential force, and relative velocity. The tangential force and slip velocity of each element at the contact site are considered to have identical tangential stress τ e and slip velocity s e . Combining Equation (1) and r(θ) yields the equation for determining wear depth ∆r(θ).
where Fτ is tangential creep force, ρ is density of wheel material, and b is average width of the wear zone. The remaining radius of the wheel circle θ is

The Conditiom of Polygonal Wear
According to Equation (2), wheel wear is proportional to the tangential creep force. When the tangential creep force changes at a given frequency, the wear depth changes as well, resulting in wheel polygonation. The peak of wear happens just once in a tangential vibration cycle when the tangential creep force reaches its maximum.
The time interval between two consecutive wear peaks is 1/f 1 if the tangential vibration frequency is f 1 . The wheel rotation duration is 1/f 2 if the wheel rotational frequency is f 2 . The peak value of wear N within a rolling circle may then be calculated.
where V is the speed of the train and D is the nominal diameter of the wheel. There are two cases: (1) The wear depth at a certain position around the wheel is the same for each wheel rolling circle when the tangential vibration frequency is an approximate integer multiple of the wheel rotational frequency. The wear peak always emerges in certain regions, and the wheel circumference will become polygonal after a long period. (2) When the tangential vibration frequency is not an approximate integer multiple of the wheel rotational frequency, and the wear depth at a certain spot around the wheel varies within two rounds of wheel rolling. The last circle's highest point of wear may wear less. After a long period of operation, the wheel circumference will be consistent, and polygonal wear will not be visible.
According to the preceding study, tangential vibration is simply a required condition for polygonal wear. The ratio of wheel tangential vibration frequency to wheel rotational frequency determines the sequence of wheel polygonal wear. As a result, the tangential vibration condition of the wheel must be analyzed to establish the range of factors that might lead to polygonation.

Analysis of Lateral Self-Excited Vibration of Wheel-Rail Contact
Tangential force is also known as creep force, and it is further subdivided into lateral creep force and longitudinal creep force. We already investigated the influence of longitudinal creep force on polygonal wheel wear [16], thus this article focuses solely on the effect of lateral creep.

Model of Lateral Self-Excited Friction Vibration of Wheel
The wheel-rail system connects the sprung mass to the under rail foundation; the vibration of the wheelset and rail subsystem is the primary source of polygonal wheel degradation. As a result, as shown in Figure 1a, the model of wheel-rail lateral self-excited friction is created in this research, using the following assumptions: (1) The wheel-rail contact is simplified to cylinder-plane contact by ignoring the slope of the tread and the curvature of the rail. (2) Owing to the symmetry of the wheelset and track construction, only half of the wheel-rail subsystem is chosen. (3) The axle is considered a massless elastic entity, with the mass centered on the wheel.
The symmetrical constraint is placed on the symmetrical surface, and the connection mode of wheel and axle is equivalent to the lateral stiffness damping provided by the flexibility of 1/2 axle and 1 wheel. (4) The sprung mass of the vehicle body and the under foundation of the rail are simplified into wheel load P, which includes static load P0 and dynamic load ∆P, and the rail rigidity is assumed to be infinite.
Machines 2022, 10, 608 4 of 17 (4) The sprung mass of the vehicle body and the under foundation of the rail are simplified into wheel load P, which includes static load P0 and dynamic load ΔP, and the rail rigidity is assumed to be infinite.
(a) (b) In Figure 1b, m, cs, and ks represent the wheel equivalent mass, damping, and stiffness, respectively. Vb is the lateral linear velocity of the wheel relative to the rail surface, and Vb = Vsinψ, where V is the speed and ψ is the rolling angle of the wheelset.
While the lateral movement of the wheel on the rail is as simple as the turn of the belt pulley (as illustrated in Figure 1b), the creep force between the wheel and the rail provides the vibration energy input. Assume that the wheel's displacement to the original position (i.e., the position of the spring and damper when no force is applied) is x and the creep speed between the wheel and the belt (rail) is as follows: Then, the lateral vibration equation of the wheel is as follows: where Fτ is the wheel-rail creep force, and Fτ = μP. The wheel-rail creep force has a significant impact on the lateral vibration of the wheel-rail. There is a complex nonlinear relationship between the creep coefficient μ and the creep rate [17,18], as shown in Figure 2. When the creep speed is between 0 and 0.37 m/s, the creep coefficient rapidly increases, the creep force is static friction, and the slope of the creep coefficient curve is positive and steep. When the creep speed exceeds 0.37m/s, the wheel starts to slide, the creep force progressively changes from static to dynamic friction, and the creep coefficient drops as the creep speed increases. The slope of the creep coefficient curve is negative at this point and subsequently declines. In Figure 1b, m, c s , and k s represent the wheel equivalent mass, damping, and stiffness, respectively. V b is the lateral linear velocity of the wheel relative to the rail surface, and V b = Vsinψ, where V is the speed and ψ is the rolling angle of the wheelset.
While the lateral movement of the wheel on the rail is as simple as the turn of the belt pulley (as illustrated in Figure 1b), the creep force between the wheel and the rail provides the vibration energy input. Assume that the wheel's displacement to the original position (i.e., the position of the spring and damper when no force is applied) is x and the creep speed between the wheel and the belt (rail) is as follows: Then, the lateral vibration equation of the wheel is as follows: where F τ is the wheel-rail creep force, and F τ = µP. The wheel-rail creep force has a significant impact on the lateral vibration of the wheelrail. There is a complex nonlinear relationship between the creep coefficient µ and the creep rate [17,18], as shown in Figure 2. When the creep speed is between 0 and 0.37 m/s, the creep coefficient rapidly increases, the creep force is static friction, and the slope of the creep coefficient curve is positive and steep. When the creep speed exceeds 0.37 m/s, the wheel starts to slide, the creep force progressively changes from static to dynamic friction, and the creep coefficient drops as the creep speed increases. The slope of the creep coefficient curve is negative at this point and subsequently declines.
It can be seen that the creep force Fτ is a function of creep speed Vr: The creep force is expanded by first-order Taylor at the creep speed: Substituting Equation (9) into Equation (7), a vibration equation based on Taylor expansion of creep force is obtained.
That is m ..
x + (c s + c f e ) .
x + k s x = 0 Machines 2022, 10, 608 where C fe is the equivalent damping of creep force and is equal to the slope of the creep speed-creep coefficient curve. The equation has an approximate solution: where A is the constant, ω is circular frequency of the system, , and φ is the phase angle. According to Equation (12), α < 0, under the initial disturbance, the vibration of the wheel will gradually converge to zero. That is, in the positive slope of the creep coefficientcreep speed curve shown in Figure 2, c s + c fe > 0, and the wheel tends to be stable.
When α > 0, in the negative slope, and c s + c fe < 0. Wheel amplitude increases exponentially with time, and the wheel tends towards unstable dynamic sliding.
When α = 0, c s + c fe = 0, the wheel amplitude does not change with time, and the wheel tends towards unstable steady sliding, resulting in the periodic change in lateral creep force Fτ, which leads to polygonal wear of the wheel. The system obviously does not contain external excitation, but the wheel produces self-excited vibration. It can be seen that the creep force Fτ is a function of creep speed Vr: The creep force is expanded by first-order Taylor at the creep speed: Substituting Equation (9) into Equation (7), a vibration equation based on Taylor expa sion of creep force is obtained.
where Cfe is the equivalent damping of creep force and is equal to the slope of the creep speed-creep coefficient curve.
The equation has an approximate solution: )

LuGre Friction Model
The production of self-excited vibration in a wheel-rail system requires a negative slope of the creep coefficient curve. Self-excited vibration causes unstable continuous sliding of the wheel, resulting in polygonal wear. It is critical to correctly duplicate the force transmission between the wheel and the rail in order to examine the self-excited vibration of the wheel. The wheel lateral vibration is mainly dry friction motion, while the LuGre model presented by Canudas, et al. [19] accurately represents the "negative slope" in the friction process and fits the experimental results nicely. As a result, the LuGre model was utilized in this work to represent lateral creep force.
where σ 0 and σ 1 are the equivalent stiffness and damping coefficient of brush sliding displacement and sliding speed after differential treatment by friction force, respectively. σ 2 is the relative viscous damping coefficient, σ 3 is the wheel-rail contact stiffness, z is the elastic deformation of brush, V s is the Stribeck speed in steady-state friction character-istics, ϕ is Stribeck index of steady-state friction, F m is the maximum static friction force (F m = µ m P 0 sinψ), F s is the sliding friction (F s = µ s P 0 sinψ), µ m is the maximum static friction coefficient, and µ s is the dynamic friction coefficient.
The vibration differential equation of the system can be obtained by combining Equations (6), (7), and (13): x−V·sin ψ| The equation describes the process of wheel lateral vibration evolution, and its equilibrium solution corresponds to the system's stable state. The periodic solution indicates that the system generates self-excited vibration, which could result in wheel polygonation.

Wheel Rail Vertical Force Model
Sleeper spacing, pits, the initial polygonal of wheel tread, or wheel eccentricity can all cause a periodic change in vertical wheel-rail contact pressure. As a result, the wheel rail vertical force model is as follows: where P 0 is wheel load and k d and N d are the dynamic load coefficient and frequency ratio, respectively.

The Lateral Stiffness of the Wheel Is Determined
Before solving the above model, the modal mass, stiffness, and damping of the wheel must be determined. These parameters are related to the flexible wheel mode with the highest lateral vibration trend. The wheel is discretized into a vibration system with n concentrated masses, each with one degree of freedom in the U direction.
The motion equation of the system is follows: The structural displacement response can be transformed from physical coordinates to modal coordinates within the linear range.
The unit vector {I}, which represents the rigid body displacement of the system, can also be expressed as the linear superposition of each modal vector {ϕ i }, where γ i is the rigid body displacement mode participation factor, which indicates the contribution amount of each restricted mode when the system is not deformed and just follows the foundation rigid body displacement. The mass matrix [M] and stiffness matrix [K] satisfy the modal orthogonality, and the damping matrix [C] is assumed to also satisfy the modal orthogonality. Then, Equation (16) can be decoupled into n independent equations. where M i , C i , and K i are defined as follows: Solve Equation (16) to obtain q i , and then use Equation (17) to obtain {U}. The structural modal vector {ϕ i } only reflects the ratio of the displacements of each node, which can be multiplied by any constant factor, so the value of γ i is not fixed. The values of M i , C i , and K i maintain the same proportion relationship, as follows:

Analysis of System Stability
Letting x 1 = x and x 3 = z, Equation (14) will be Substituting .
x 3 = 0 into Equation (22), the equilibrium point is At the equilibrium point, Taylor expansion is done, and the linear component is inserted into Equation (22). The following is the related linear system state equation: The Jacobi matrix of the equation of state is The eigenroot equation of the matrix is According to the Routh-Hurwitz criterion, the determinant of the coefficients of Equation (27) is D. According to Equation (29) = 0, the critical conditions of system Hopf bifurcation are as follows: Substituting Equation (28) into Equation (30), the following can be obtained

Wheel Modal Parameters
The CAD model of CRH 3 wheels was created and loaded into CAE tools for the first 16 modes' analysis, the ninth of which is seen in Figure 3. The mode of this order can be observed to be wheel umbrella deformation, and the lateral deformation of the wheel hub is visible, which adds significantly to the lateral vibration. As a result of substituting the natural frequency (580.56 Hz) and modal mass (383 kg) of the mode of this order into Equation (20), the modal stiffness is 5.07 × 10 9 N/m. es 2022, 10, 608 natural frequency (580.56 Hz) and modal mass (383 kg) of Equation (20), the modal stiffness is 5.07 × 10 9 N/m.

System Bifurcation Point Is Determined
Draw crucial curves for the speed ratio β, damping ra

System Bifurcation Point Is Determined
Draw crucial curves for the speed ratio β, damping ratio ξ, and wheel load P 0 , as indicated in Figure 4. The parameter domain is divided into three sections by the critical border of parameters β and P: the area of β < β 1 and β > β 2 is stable, while the region of β 1 < β < β 2 is unstable. The instability zone grows as P 0 increases.

System Bifurcation Point Is Determined
Draw crucial curves for the speed ratio β, damping ratio ξ, and wheel load P0, as indicated in Figure 4. The parameter domain is divided into three sections by the critical border of parameters β and P: the area of β < β1 and β > β2 is stable, while the region of β1 < β < β2 is unstable. The instability zone grows as P0 increases.  The alterations in the ξ-β critical curve are comparable to those in the P0-β critical curve. The instability zone shrinks fast as ξ increases. The vibration equation is solved and the impact of each parameter is explored to validate this result and determine the changing characteristics at the critical point.

Speed
Taking the speed ratio β = 0.1~2.0 to investigate the vibration characteristics at various rotating speeds, the bifurcation diagram in Figure 5 shows that, when the speed is low, the wheel cannot create self-excited vibration. When β = 0.2, the stable equilibrium point loses stability, and Hopf bifurcation leads to periodic oscillation, i.e., the wheel generates self-excited vibration. The system becomes stable again at β = 1.8 and is increased further. The alterations in the ξ-β critical curve are comparable to those in the P 0 -β critical curve. The instability zone shrinks fast as ξ increases. The vibration equation is solved and the impact of each parameter is explored to validate this result and determine the changing characteristics at the critical point.

Speed
Taking the speed ratio β = 0.1~2.0 to investigate the vibration characteristics at various rotating speeds, the bifurcation diagram in Figure 5 shows that, when the speed is low, the wheel cannot create self-excited vibration. When β = 0.2, the stable equilibrium point loses stability, and Hopf bifurcation leads to periodic oscillation, i.e., the wheel generates self-excited vibration. The system becomes stable again at β = 1.8 and is increased further.
Machines 2022, 10, 608  The time domain diagrams and phase diagrams at different speeds in Figure 6 show that, when the speed is low (V = 36 km/h), the time domain diagrams increasingly converge. The stable focus in the phase diagram is the equilibrium point, and the wheel cannot create self-excited vibration. The time domain figure displays the constant amplitude and the phase diagram illustrates the elliptic limit cycle when V = 180 km/h. The system loses stability when the stable equilibrium point becomes unstable, and the wheel creates selfexcited vibration. When the speed is increased further, to V = 324 km/h, supercritical Hopf bifurcation returns to the system. The time domain diagram progressively converges, the limit cycle disappears in the phase diagram, and the equilibrium point reverts to being the stable point. The system is generally stable, and the wheel cannot produce self-excited vibration. The time domain diagrams and phase diagrams at different speeds in Figure 6 show that, when the speed is low (V = 36 km/h), the time domain diagrams increasingly converge. The stable focus in the phase diagram is the equilibrium point, and the wheel cannot create self-excited vibration. The time domain figure displays the constant amplitude and the phase diagram illustrates the elliptic limit cycle when V = 180 km/h. The system loses stability when the stable equilibrium point becomes unstable, and the wheel creates self-excited vibration. When the speed is increased further, to V = 324 km/h, supercritical Hopf bifurcation returns to the system. The time domain diagram progressively converges, the limit cycle disappears in the phase diagram, and the equilibrium point reverts to being the stable point. The system is generally stable, and the wheel cannot produce self-excited vibration.
As a result, while the vehicle is traveling at speeds ranging from 36 to 324 km/h, the wheel can create self-excited vibration, implying that polygonal wear may occur.

Wheel Load
The vibration characteristics of a wheel load at P0 = 65~80 kN were investigated using β = 1.6. The bifurcation diagram in Figure 7 demonstrates this. With the increase in wheel load, the stable equilibrium point loses stability when P0 = 67 kN, and Hopf bifurcation leads to periodic oscillation, that is, the wheel generates self-excited vibration.
When P0 < 67 kN, the time domain graph progressively converges, the phase diagram's equilibrium point is the stable focus, and the wheel cannot generate self-excited vibration. As the speed increases, the time domain graphic displays the constant amplitude waveform and the phase diagram reveals the elliptic limit cycle in P0 > 67 kN. The system loses stability when the stable equilibrium point becomes the unstable focus, and the wheel creates self-excited vibration. As a result, while the vehicle is traveling at speeds ranging from 36 to 324 km/h, the wheel can create self-excited vibration, implying that polygonal wear may occur.

Wheel Load
The vibration characteristics of a wheel load at P 0 = 65~80 kN were investigated using β = 1.6. The bifurcation diagram in Figure 7 demonstrates this. With the increase in wheel load, the stable equilibrium point loses stability when P 0 = 67 kN, and Hopf bifurcation leads to periodic oscillation, that is, the wheel generates self-excited vibration.
vibration. As the speed increases, the time domain graphic displays the constant amplitude waveform and the phase diagram reveals the elliptic limit cycle in P0 > 67 kN. The system loses stability when the stable equilibrium point becomes the unstable focus, and the wheel creates self-excited vibration.  When P 0 < 67 kN, the time domain graph progressively converges, the phase diagram's equilibrium point is the stable focus, and the wheel cannot generate self-excited vibration. As the speed increases, the time domain graphic displays the constant amplitude waveform and the phase diagram reveals the elliptic limit cycle in P 0 > 67 kN. The system loses stability when the stable equilibrium point becomes the unstable focus, and the wheel creates self-excited vibration.

Wheel Rail Vertical Dynamic Force
Using a speed ratio β = 1.6 and a wheel load P 0 = 80 kN, the vibration characteristics are investigated using dynamic load coefficients k d = 0.0025, 0.025, and 0.25, as well as frequency ratios n d = 1 and 3.
The spectrum of the system with the amplitude and frequency of the wheel load dynamic load shown in Figure 8a,b shows that the wheel dynamic load excites the same frequency response. The amplitude of the wheel dynamic load influences the low-frequency response amplitude, and the frequency of the wheel dynamic load influences the lowfrequency response frequency. Both have no effect on the size of higher-order frequency, implying that the vertical force of the wheel rail has no effect on the magnitude of fixed frequency. According to Figure 8c, the polygonal wear consists of a low order (3rd) with the same frequency as the wheel rail vertical force and a high order (20th) with the fixed frequency.
The bifurcation diagram of the system changing with the damping ratio shown in Figure 9 shows that, when the damping ratio is small, the wheels generate self-excited vibration; when the damping ratio is increased to 0.09, subcritical Hopf bifurcation occurs, the equilibrium point becomes the stable focus, and the wheel does not generate self-excited vibration.
The bifurcation diagram of the system changing with the damping ratio shown in Figure 9 shows that, when the damping ratio is small, the wheels generate self-excited vibration; when the damping ratio is increased to 0.09, subcritical Hopf bifurcation occurs, the equilibrium point becomes the stable focus, and the wheel does not generate self-excited vibration.
The bifurcation diagram of the system changing with the damping ratio sh Figure 9 shows that, when the damping ratio is small, the wheels generate self vibration; when the damping ratio is increased to 0.09, subcritical Hopf bifurcation the equilibrium point becomes the stable focus, and the wheel does not generate cited vibration.

Numerical Simulation Verification
The lateral creep force spectrum within the speed range of the self-excited vibration is drawn to verify the theoretical analysis conclusion that the wheel polygonal wear will evolve when the wheel rotational frequency is divided into the wheel lateral vibration frequency after the self-excited vibration of the wheel is generated, as shown in Figure 10a. The vibration frequency of the lateral creep force with periodic oscillation generated by Hopf bifurcation is independent of vehicle speed, and it is fixed around 588.47 Hz at various vehicle speeds. This finding is congruent with the conclusion in [20] that polygonal wheel wear is "frequency fixed". Given the wheel mass and lateral stiffness, the lateral vibration frequency is f 1 = 1 2π k s m = 580.4 Hz. It is clear that the frequency corresponds to the natural frequency of wheel lateral vibration.
Furthermore, within the speed range that creates self-excited vibration, the lateral creep force is substituted into the wheel circumferential wear formula, and the polygon wear characteristics of the wheel are shown. The wear characteristics in Figure 10b-i show that, when the vehicle speed is 36 or 324 km/h, that is, when the wheel does not have lateral self-excited vibration, the circumferential wear is uniform. The wheel wear will finally be uniform after long-term operation when the vehicle speed is 198, 234, and 270 km/h; however, there are evenly distributed wear peaks in the wheel circumferential direction only when the vehicle speed is 252, 288, and 308 km/h. After a long period of operation, obvious regular 24th, 21st, and 20th order polygons will be formed, and the product of wheel rotational frequency and wear order is approximately equal to 588 Hz, indicating that the wheel polygonal wear exhibits the "divisible" feature.
To summarize, when the wheel has polygon wear, it exhibits "self-excitation-constant speed-fixed frequency-divisible" properties. That is, the incidence and evolution of wheel polygonal wear must fulfill the criteria that, while the wheel is moving at a constant speed, it creates self-excited periodic vibration with a constant frequency, and the frequency of the self-excited vibration is divided by the wheel rotational frequency. The findings of the simulation validate the theoretical analysis.
by Hopf bifurcation is independent of vehicle speed, and it is fixed around 588.47 Hz at various vehicle speeds. This finding is congruent with the conclusion in [20] that polygonal wheel wear is "frequency fixed". Given the wheel mass and lateral stiffness, the lateral vibration frequency is Furthermore, within the speed range that creates self-excited vibration, the lateral creep force is substituted into the wheel circumferential wear formula, and the polygon

Actual Vehicle Tracking Verification
The CRH 3 EMU has a top speed of 300 km/h. According to field measurements, the wheel diameter life cycle has three stages of high-speed growth of polygonal wheel wear. The wheel diameters are 830, 875, and 915mm, respectively, corresponding to polygon wear of 18th, 19th, and 20th orders [2]. The rotational frequency f 2 and lateral vibration frequency f 1 of the matching wheelset may be estimated using the present vehicle recorded data, as given in Table 1.  Table 1 shows that polygonal wear occurs when the train is traveling at 300 km/h and the wheel diameters are 830, 875, and 915 mm, respectively. There is a constant frequency of 580 Hz in all three circumstances. The "fixed frequency" is equal to the multiple of the corresponding wheel rotational frequency, and the integral multiple causes 18th, 19th, and 20th order wheel polygon wear. It is clear from the measured data of current automobiles that the polygonal wear of wheels exhibits apparent "fixed frequency" features. This is consistent with the study result.

Conclusions
The properties of lateral self-excited vibration and polygonal wear of wheels are explored using the LuGre friction model and self-excited vibration theory, and the conditions for the generation and evolution of polygonal wear are disclosed. The following findings are reached after verifying the measured data of current vehicles: