Stick–slip is considered the root cause of railway engineering phenomena such as squeal noise and corrugation. Little consensus regarding the actual physical description of stick–slip has been achieved because its manifestations cannot be explained by a single underlying mechanism. To investigate the generation mechanisms of stick–slip contact, this study reproduced wheel-rail stick–slip experimentally with an in-house test rig—V-Track and numerically with an explicit finite element method (FEM). The V-Track is capable of reproducing realistic wheel-rail dynamic interactions under well-controlled lab conditions, while the explicit FEM has been proven to be suitable for the simulation of dynamic contact and frictional instability. Crucial influential factors including wheel-rail lateral creepage, friction levels, and friction characteristics were varied in the experiments and simulations to examine their impacts on the occurrence of stick–slip. The study shows that the creepage level needs to be sufficiently high to generate stick–slip. Stick–slip can be eliminated by reducing friction to a very low level, whereas changing the friction characteristics from negative to positive may not work for stick–slip mitigation. Moreover, wheel-rail friction conditions cannot be sufficiently represented by a single parameter, i.e. coefficient of friction.

Frictional heat is generated at the wheel-rail interface during train operations, particularly under high slip ratios during acceleration and braking. Thermal effects can accelerate wear, induce plastic deformation, and contribute to thermal fatigue. Reliable modelling of wheel-rail contact that considers friction-induced thermal effects is desirable for the accurate prediction of wheel-rail interface deterioration. Several analytical and numerical models have been proposed to simulate thermal or thermomechanical wheel-rail loads but have rarely been validated, especially in high slip ratio scenarios where flash temperatures exceed 200 °C. This study develops and experimentally validates a three-dimensional thermomechanical finite element (FE) wheel-rail contact model for high slip ratio conditions, with contact temperatures reaching 360 °C. The model incorporates key mechanical parameters, including wheel loads, coefficients of friction, and slip ratios. Simulated rail surface temperatures across various slip ratios (5 %, 10 %, and 15 %) are compared with the flash temperatures measured with an onboard infrared thermal camera, showing good agreement with a maximum deviation of 9.9 %. This confirms the reliability of the model for simulating wheel-rail contact under thermal effects.

Short pitch corrugation is a typical rail defect that lacks a thorough understanding and adequate root-cause solutions. This paper aims to identify the damage mechanism of short pitch corrugation through a microstructural analysis of a field rail sample. This sample made of R260Mn pearlitic steel was taken from a straight section of the Dutch railway network, and its geometry and surface hardness variation along the corrugation were measured and analyzed. Eleven specimens, including both corrugated and non-corrugated zones, were sectioned from the rail sample and continuously examined using light optical microscopy, scanning electron microscopy and micro-hardness testing. The results indicate that the corrugation damage mechanism can be categorized into three stages: (1) pre-corrugation, characterized by uniform wear and plastic deformation; (2) corrugation initiation, dominated by differential wear; and (3) corrugation growth, involving both differential wear and plastic deformation accumulation. The initiation and growth of corrugation both contribute to an inhomogeneous distribution of plastic deformation layer (PDL) in the rail subsurface, which follows an approximately sinusoidal pattern, matching the corrugation geometry in both wavelength and phase. Consequently, the hardness also varies in phase with the corrugation geometry, with higher hardness values at corrugation peaks. In the non-corrugation zone, the PDL and hardness show relatively small and irregular fluctuations. This study also provides meaningful insights into rail grinding, suggesting that grinding should account for differential PDL thickness to prevent corrugation reoccurrence due to subsurface material inhomogeneity.

Short pitch corrugation is a major rail defect worldwide, and its development mechanism remains not fully understood. This work aims to better understand corrugation and validate some numerical predictions through extensive field monitoring. 105 corrugations on four sections of mainline tracks were continually monitored for five years. Comprehensive field data were collected, including photos, geometry, and vehicle-track dynamic responses. Numerical results of corrugation development agree with the field observations. It confirms that corrugation initiates with necessary initial excitation, consistently grows at fixed locations due to differential wear, and eventually reaches a limiting amplitude. Moreover, vehicle-track longitudinal vibrations are crucial to corrugation initiation, while vertical vibration plays an increasingly important role in corrugation growth.

This paper presents a methodology for detecting and monitoring short pitch corrugation (SPC) under varying measurement conditions using vertical and longitudinal axle box acceleration (ABA) measurements. The main objective of the detection algorithm is to determine the likelihood and approximate severity of SPC presence, providing insights for maintenance planning. The methodology combines a validated three-dimensional finite element (3D-FE) model of the ABA responses at SPC and signal processing techniques to extract meaningful data from the real-world on-board measurements. First, a 3D-FE vehicle-track model is validated and used to quantify the physical relationships in the time–frequency responses of ABA at SPC under different levels of corrugation severity and measurement speeds. Then, a measurement train is instrumented with multiple accelerometers to capture field data on ABA at SPC, which is validated with field inspections and Railprof measurements. Finally, the ABA responses are analyzed based on the number of signals detecting SPC and an assessment of severity based on impact energy due to SPC. The methodology is demonstrated by analyzing the track between Assen and Groningen on the Dutch rail network. Results show that the methodology accurately detects registered SPC locations. Further, a whole track analysis is conducted, from which the methodology proposes new locations and severities of SPC, providing crucial information for rail maintenance planning.

Polygonal wear is a type of damage commonly observed on the railway wheel tread. It induces wheel-rail impacts and consequent train/track components failure. This study presents a finite element (FE) thermomechanical wheel-rail contact model, which is able to cope with the three possible generation and development mechanisms of polygonal wear: initial defects, thermal effect, and structural dynamics. The polygonal wear-induced impact contact and further development of wear are simulated. The simulated elastic contact solutions are verified against the program CONTACT. Different material properties (elastic, elasto-plastic and elasto-plastic-thermo, i.e. with thermal softening) and initial polygonal profiles are then applied to the FE model to investigate the influence of wheel/rail material and wear amplitude on wheel-rail contact stress and wear development. The simulations indicate that the wheel-rail impact-induced temperature may reach up to 362 ℃ at the contact interface, and the high temperature at the contact area influences wheel-rail contact stress and wear depth.

Wheel-rail high-frequency interaction is closely related to the formation of railway short-wave defects. Finite element (FE) method has been widely used to simulate wheel-rail dynamic systems, but its validity in modelling high-frequency interaction has not been fully demonstrated in three dimensions (3D). This work aims at comprehensively validating the 3D FE modelling of wheel-rail high-frequency interaction using a downscale V-Track test rig. First, the FE model of the V-Track is developed that comprehensively includes the 3D track elasticity. The simulated track dynamic behaviours are validated against hammer tests, and the major vibration modes are analyzed employing modal analysis. Afterwards, the simulate wheel-rail dynamic responses are comprehensively compared with measurement results up to 10 kHz. Their characteristic frequencies are identified and correlated to the eigenmodes of the vehicle-track system. The results indicate that the proposed 3D FE model is capable of comprehensively and accurately simulating the 3D track dynamics and wheel-rail dynamic interaction of the V-Track up to 10 kHz. Rail vibrations dominate the wheel-rail dynamic contact within 10 kHz, while the wheel vibrations play an increasingly important role at higher frequencies and become decisive near the wheel eigenmode frequencies. The V-Track overall achieves dynamic similarity to the real vehicle-track system.

This paper aims at identifying the formation mechanism of short pitch rail corrugation using an improved configuration of a 1/5 scaled V-Track test rig. The loading conditions of the V-Track are designed to simulate the wheel-rail dynamic interactions on the tangent track. Corrugation initiation and consistent growth are successfully produced on the rail surface with a characteristic wavelength of about 6.0 mm, equivalent to the field one of 30 mm. By analysing the wheel-rail dynamic forces, the wheel and track dynamic behaviours across different rotations, it is found that the rail longitudinal vibration mode and its induced longitudinal dynamic force determine the corrugation formation with the necessary initial excitation.

Short pitch corrugation has been a problem for railways worldwide over one century. In this paper, a parametric investigation of fastenings is conducted to understand the corrugation formation mechanism and gain insights into corrugation mitigation. A three-dimensional finite element vehicle–track dynamic interaction model is employed, which considers the coupling between the structural dynamics and the contact mechanics, while the damage mechanism is assumed to be differential wear. Various fastening models with different configurations, boundary conditions, and parameters of stiffness and damping are built up and analysed. These models may represent different service stages of fastenings in the field. Besides, the effect of train speeds on corrugation features is studied. The results indicate: (1) Fastening parameters and modelling play an important role in corrugation formation. (2) The fastening longitudinal constraint to the rail is the major factor that determines the corrugation formation. The fastening vertical and lateral constraints influence corrugation features in terms of spatial distribution and wavelength components. (3) The strengthening of fastening constraints in the longitudinal dimension helps to mitigate corrugation. Meanwhile, the inner fastening constraint in the lateral direction is necessary for corrugation alleviation. (4) The increase in fastening longitudinal stiffness and damping can reduce the vibration amplitudes of longitudinal compression modes and thus reduce the track corrugation propensity. The simulation in this work can well explain the field corrugation in terms of the occurrence possibility and major wavelength components. It can also explain the field data with respect to the small variation between the corrugation wavelength and train speed, which is caused by frequency selection and jump between rail longitudinal compression modes.

Short pitch corrugation is a typical defect on rail surfaces that induces high level of noise and increases maintenance costs. Despite numerous research efforts, corrugation development mechanism has not yet been fully understood and root-cause solutions have not been developed. This work numerally simulates the rail corrugation in the V-Track test rig, aiming to better understand corrugation mechanism and also link the scaled laboratory tests to the full-scale reality. A three-dimensional finite element model of the V-Track is established to simulate the vehicle-track dynamic interaction. The fastening and ballast parameters are calibrated by fitting the simulated track receptances to hammer tests. Rail corrugation with a major wavelength of 5.7 mm is successfully reproduced using the FE model, which shares features similar to the experimentally produced corrugation in the V-Track. The numerical simulation demonstrates that a rail longitudinal compression mode at 790 Hz is the ‘wavelength-fixing’ mechanism of corrugation in the V-Track, agreeing with the experimental results. This work numerically verifies the dominance of the rail longitudinal vibration modes on corrugation formation.

While various train-borne techniques have been developed for measuring railway track stiffness, differentiating stiffness at different track layers remains a challenge. This study proposes a digital twin framework for the vehicle–track interaction system, which enables track stiffness evaluations based on axle box accelerations (ABA). The digital twin consists of a physics-based model, a model library and data-driven models. Compared to existing techniques, the proposed method simultaneously evaluates the stiffness of the railpad, sleeper and ballast layers at a sleeper spacing resolution, while being robust to varying track conditions, such as track irregularities and vehicle speeds. This is accomplished by employing a localized frequency-domain ABA feature capable of distinguishing between the characteristics of different track layers. Furthermore, track stiffness is evaluated in near real-time. This is achieved using a model library derived from physics-based simulations of a range of track conditions. Two data-driven models that can quickly select or interpolate model instances contained in the library are developed. During operation, the data-driven models use the measured ABA features as input and then infer the stiffness for the different track layers. The proposed method is applied to evaluate the track stiffness of a downscale test rig in a case study. The track stiffness evaluated by the proposed method is compared with that obtained through hammer tests and with the observations of the track component conditions. These comparisons show that the proposed method can capture the stiffness variations due to periodically fastened clamps and substructure misalignments at different speeds. In addition, the proposed method is demonstrated to be superior to the commonly used hammer test method for evaluating track stiffness under loaded conditions.

In practice, the assessment and treatment of rail corrugation are quantitatively based on the corrugation depth. Wheel–rail vertical forces (WRVF), as a direct reflection of wheel–rail interaction, can give expression to the corrugation depth and thus serve as a key parameter for assessing the corrugation. In this paper, we propose an evaluation method for rail corrugation based on the WRVF. First, a 3D wheel–rail dynamic finite element (FE) model was developed with typical parameters of CRTS II slab track and CRH3 vehicle for high-speed railways in China. The accuracy of the model was then validated with the measured WRVF data in the field. Second, using the validated model, the time–frequency domain distribution of WRVF (vehicle speed: 300 km/h) was obtained with consideration of the corrugation wavelength in the range of 40–180 mm. The non-linear least squares method and rational equation were used to fit the function between the large value of WRVF and the corrugation depth value under the conditions of different corrugation wavelengths. Next, effects of the Pinned–Pinned resonance frequency and vibration mode on the fitted parameters were analysed, by which an indicator for corrugation treatment (grinding) was proposed. Finally, the indicator was applied in the monitoring of rail corrugation for high-speed railway lines in the field. The results show that the misjudgement rate of rail grinding decisions (using the proposed indicator) is low with the accuracy at 92.5%. The proposed method can provide a basis for the rail corrugation evaluation and grinding decisions-making.

Short pitch corrugation is a quasi-periodic rail defect that induces a high level of noise and accelerates track degradation. This paper proposes a methodology to mitigate short pitch corrugation by rail constraint design, including four steps. In Step 1, corrugation is numerically reproduced by employing a three-dimensional (3D) finite element (FE) vehicle-track model with degraded fastenings. In Step 2, the corrugation initiation mechanism is identified by the operating deflection shapes (ODSs) approach. In Step 3, different types of rail constraints are designed and their effects on rail vibration modes are analysed. Then FE models of these rail constraints are built up and validated. In Step 4, rail constraint models from Step 3 are applied to the 3D FE vehicle-track interaction model and their effectiveness on corrugation mitigation is evaluated. The results indicate rail longitudinal compression modes and the induced longitudinal dynamic contact force dominate the initial differential wear and corrugation initiation. Based on this mechanism, a new rail constraint is designed in this work that can completely suppress longitudinal compression modes and significantly reduce the fluctuation amplitude of the longitudinal contact force so that corrugation can hardly initiate. This paper first points out a direction for field corrugation mitigation by strengthening the rail longitudinal constraint.

The coupled vehicle-track system (CVTS) dynamics have been extensively investigated for decades. However, the calculation accuracy of prevailing vehicle-track coupling models needs to be improved in the high frequency range due to the inappropriate model simplification and neglect of material nonlinearity. In this study, we propose a refined numerical model of the CVTS that considers the nonlinear properties of the railpads and primary suspension using the fraction derivative Zener model. Furthermore, we more realistically simulate the wheelset, rail and railpad configuration with the elastic axle, solid finite element and surface-support models, respectively, and improve the computation efficiency by employing the mode superposition method. The results demonstrate that the refined CVTS model is more accurate than the classical model in simulating vehicle-track coupling dynamics above 2 kHz. In particular, there are significant differences in the dynamic response of the elastic wheelset model compared to the rigid model over a broad frequency range, with an 11% difference in the bogie acceleration response at the first dominant frequency. When the railpads are modeled using the surface-support model, the rail acceleration differences exceed 41% near 1 kHz and 44% near 2650 Hz, compared to the point-support model. Additionally, the rail response at various locations across the rail cross section can be calculated using the finite element method in this refined model. Overall. the proposed CVTS model provides high accuracy and efficiency for random vibration analysis, especially in the high frequency domain.

This paper presents an experimental study on the development mechanism of short pitch corrugation using an innovative downscale V-Track test rig. Rail corrugation is reproduced in the V-Track with two wavelengths of 6.4 mm and 8.9 mm. Based on the measured wheel-rail contact forces, track dynamic behaviors, and observations, the consistent development mechanisms of short pitch corrugation are identified. It is found the longitudinal mode at 574 Hz and the induced longitudinal dynamic force are dominant for corrugation formation at the major wavelength of 6.4 mm. This paper provides the first experimental evidence indicating that rail longitudinal vibration mode and its induced longitudinal force contribute to the consistent development of short pitch corrugation.

Wheel-rail impact, which arises from structural discontinuities and short-wave defects, becomes more serious with higher train speeds and larger axle loads. The large impact load can accelerate the deterioration of vehicle-track components and induce a high level of impact noise. This paper aims to better understand the characteristics of wheel-rail impact and redevelop corresponding mitigation measures. First, a well-validated vehicle-track vertical interaction model with a singular rail defect is built up considering nonlinear Hertz contact. Then, the simulated wheel-rail impact force is characterized in both the time and frequency domains with different running speeds and defect geometries employing continuous wavelet transform. The identified characteristic frequencies are correlated to the track resonance modes. Afterward, a parameter sensitivity analysis of railpads, ballast, roadbed, and suspensions is performed to obtain the mitigation measures of wheel-rail impact. The results show that the wheel-rail impact force can be characterized by four stages in the time domain, the quasi-static stage before the impact, the forced vibration, the free vibration, and the quasi-static stage after the impact, respectively. Four characteristic frequencies are identified in the wheel-rail impact response: f₁ at 45 Hz, f₂ at 100 Hz, f₃ at 260 Hz, and f₄ at 810 Hz. Among them, f₄ has the dominant vibration energy and is determined together by the pinned-pinned resonance mode and the defect excitation frequency. Characteristic frequencies f₁, f₂, and f₃ correspond to the ballast, sleeper, and rail resonance modes, respectively, which are independent of the defect geometry. The increase of railpad stiffness can effectively reduce the maximum wheel-rail impact force and thus the impact factor. Larger railpad stiffness and damping can significantly reduce the dominant vibration energy at about 810 Hz of f₄. This work can contribute to the optimization of vehicle-track parameters for a new design of more impact-resistant railways.

This paper proposes a new hypothesis for the formation process of short pitch rail corrugation. An FE wheel-track dynamic model is utilized to verify the hypothesis by reproducing corrugation initiation and consistent growth. It is found longitudinal compression modes are responsible for corrugation initiation with necessary initial excitation that allows flexibility for longitudinal vibration. Consistency between longitudinal compression and vertical bending eigenfrequencies of the wheel-track system is required for consistent corrugation growth, which also determines maximum corrugation amplitude. Corrugation initiates by frequency selection instead of wavelength fixing. The proposed mechanism can explain field observations including the wavelength and periodicity of corrugation in the Netherlands, why corrugation forms on continuously-supported tracks where pinned-pinned resonance does not exist, and the small variation between the corrugation wavelength and train speed.

Short pitch corrugation is a (quasi-) periodic rail surface defect with shiny crests and dark valleys. It primarily occurs on tangent tracks or gentle curves with a typical wavelength in the range of 20-80 mm. Short pitch corrugation excites high-frequency wheel-rail dynamic contact forces and generates a high level of noise, which is a nuisance to both the passengers and the residents near the railway lines. The resulting large dynamic forces accelerate the degradation of the track components and may induce other rail defects (such as, squats), which increase the maintenance cost. The goal of this dissertation is to better understand the formation mechanism of short pitch corrugation and develop the root-cause solutions to mitigate it. Three steps are taken to achieve this goal: 1) identification and control of rail vibration modes which are crucial to short pitch corrugation formation; 2) design of a new rail constraint to mitigate short pitch corrugation; 3) experimental study of short pitch corrugation using an innovative V-Track test rig. Step 1 focuses on the identification and control of rail vibration modes. First, the vibration modes and dispersive waves of a free rail are simulated employing a finite element (FE) approach. The modal behaviors, wavenumber-frequency dispersion relations, and phase and group velocities of six types of propagative waves are derived and discussed in detail in 0-5 kHz. The operating deflection shape (ODS) approach distinguishes different types of rail vibration modes experimentally. A synchronized multiple-acceleration wavelet (SMAW) approach is proposed to experimentally study the propagation and dispersion characteristics of these waves. Both the laboratory and in-situ experimental results demonstrate the effectiveness of the ODS measurement for coupled rail mode identification and the SMAW approach for wave dispersion analysis. Afterward, the ODS and SMAW approaches are further applied to investigate rail vibration modes and wave propagation under fastening constraint. A three-dimensional (3D) FE rail-fastening model is also developed and validated against the ODS and SMAW measurement results. Subsequently, a sensitivity analysis of fastening parameters using this FE model is performed to gain insights into the control of rail vibrations. The results indicate that under fastening constraint, ODS measurement identifies vertical bending modes, longitudinal compression modes and lateral bending modes with shifted frequencies and significantly reduced vibration amplitude compared to free rail. Fastenings constrain the rail longitudinal vibrations less strongly compared to the vertical and lateral directions. The variation of fastening parameters can control rail mode frequencies and their vibration amplitudes, and influence the wave propagation velocities and attenuation along the rail. Step 2 proposes a methodology to design a new rail constraint to mitigate short pitch corrugation. First, a parametric investigation of fastenings is conducted to understand the corrugation development mechanism and gain insight for a new rail constraint design for corrugation mitigation. A 3D FE vehicle-track dynamic interaction model is employed, which considers the coupling between the structural dynamics and the contact mechanics, and the damage mechanism is assumed to be differential wear. Various fastening models with different configurations, boundary conditions, and dynamic parameters are built up and analyzed. The results indicate that the fastening longitudinal constraint to the rail is the major factor determining the corrugation development. The fastening vertical and lateral constraints influence corrugation features in terms of spatial distribution and wavelength components. The increase of fastening constraint in the longitudinal dimension helps to mitigate corrugation, and the inner fastening constraint in the lateral dimension is necessary for corrugation alleviation. Based on these insights, a methodology is proposed to mitigate short pitch corrugation by rail constraint design. First, short pitch corrugation is numerically reproduced employing a 3D FE vehicle-track interaction model. Then, the corrugation initiation mechanism is identified by examining the ODSs of rail longitudinal compression modes. Afterward, different rail constraints are designed, and their effects on longitudinal compression modes are analyzed. Models of these rail constraints are also built and validated. Finally, the rail constraint models are applied to the 3D FE vehicle-track interaction model, and their validity on short pitch corrugation mitigation is evaluated. It is found that a relative rigid constraint can completely suppress rail longitudinal compression modes and significantly reduce the fluctuation amplitude of the longitudinal contact force to mitigate corrugation. A direction is pointed out for corrugation mitigation in the field by strengthening the rail longitudinal constraint. Step 3 performs an experimental investigation of short pitch corrugation using the downscale V-Track test rig. First, a force measurement system named dynamometer is developed in the V-Track to measure the wheel-rail contact forces for short pitch corrugation experiments. The dynamometer consists of four 3-component piezo-electric force sensors and is mounted between the wheel assembly and the steel frame, enabling it to measure the forces transmitted from the wheel-rail interface to the frame. Static tests are first carried out to calibrate the dynamometer in three directions. Then, several tests are performed in the V-Track to examine the reliability and validity of the dynamometer for measuring the wheel-rail contact forces under running conditions. Experimental results show that the dynamometer is capable of reliably and accurately measuring these forces. Utilizing the measurement results from the dynamometer, the control of the wheel-rail contact forces in V-Track has also been achieved. Afterward, the V-Track test rig is used to investigate the development mechanism of short pitch corrugation experimentally. The loading conditions of the V-Track are designed to simulate the vehicle-track interaction on tangent tracks where short pitch corrugation mainly occurs in the field. Short pitch corrugation is successfully reproduced in the V-Track, and its spatial distribution, wavelength components, and hardness variation are captured by the 3D HandyScan and the hardness tests. Based on the measurement results of wheel-rail contact forces and track dynamic behaviors and observations, the development mechanism of short pitch corrugation is identified. It is found that rail longitudinal and lateral vibration modes contribute to the consistent development of short pitch corrugation. Overall, the major contribution of this dissertation is threefold: 1) a better understanding of vibration modes and wave propagation of the rail in free condition and under fastening constraint is obtained by ODS and SMAW measurement, which is essential to understand and mitigate short pitch corrugation; 2) a new rail constraint is designed which can effectively suppress rail longitudinal compression modes and mitigate short pitch corrugation; 3) experimental evidence is provided to demonstrate that initial excitation and longitudinal compression modes play a significant role in the consistent growth of short pitch corrugation.

By modifying friction to the desired level, the application of friction modifiers (FMs) has been considered as a promising emerging tool in the railway engineering for increasing braking/traction force in poor adhesion conditions and mitigating wheel/rail interface deterioration, energy consumption, vibration and noise. Understanding the effectiveness of FMs in wheel–rail dynamic interactions is crucial to their proper applications in practice, which has, however, not been well explained. This study experimentally investigates the effects of two types of top-of-rail FM, i.e. FM-A and FM-B, and their application dosages on wheel–rail dynamic interactions with a range of angles of attack (AoAs) using an innovative well-controlled V-track test rig. The tested FMs have been used to provide intermediate friction for wear and noise reduction. The effectiveness of the FMs is assessed in terms of the wheel–rail adhesion characteristics and friction rolling induced axle box acceleration (ABA). This study provides the following new insights into the study of FM: the applications of the tested FMs can both reduce the wheel–rail adhesion level and change the negative friction characteristic to positive; stick–slip can be generated in the V-Track and eliminated by FM-A but intensified by FM-B, depending on the dosage of the FMs applied; the negative friction characteristic is not a must for stick–slip; the increase in ABA with AoA is insignificant until stick–slip occurs and the ABA can thus be influenced by the applications of FM.

Polygonal wear is a common type of damage on the railway wheel tread, which could induce wheel-rail impacts and further components failure. This study presents a finite element (FE) thermomechanical model to investigate the causes of wheel polygonal wear. The FE model is able to cope with three possible causes of polygonal wear: thermal effect, initial defects, and structural dynamics. To analyse the influences of the three causes on wheel-rail contact stress and wear depth, different material properties (i.e., elastic, elasto-plastic, thermo-elasto-plastic with thermal softening), and wheel profiles (i.e., round and polygonal) were used in the FE model. The simulation indicates that a high temperature up to 264.20 ℃ could be induced by full-slip wheel-rail rolling contact when the polygonal profile is used. The thermal effect, similar to that induced by tread brake, may then have a significant influence on wheel-rail contact stress and wear depth. In addition, the involvement of initial defects, i.e., polygonal profile, causes wheel-rail impact contact and remarkably increases the contact stress and wear. By reliably considering all the three possible causes, the proposed FE model is believed promising for further explaining the generation mechanisms of wheel polygonal wear.