SinoTechIntel Academic Portal
Open AccessDOI: 10.1007/s40534-025-00405-4Original Research

Computational methods to predict RCF crack initiation hot spots in rails using critical plane SWT damage indicator parameter

KTH Royal Institute of Technology, Department of Engineering Mechanics, Stockholm, Sweden

Read Executive PreviewQuick FAQ
Computational methods to predict RCF crack initiation hot spots in rails using critical plane SWT damage indicator parameter
Graphical Abstract / Figure
Published In
Railway Engineering Science (铁道工程科学)
Published:January 15, 2026Edition:Vol 34, Issue 2 • pp. 100-112Citation:Jonathan Leung et al. (2026), Railway Engineering Science (铁道工程科学)
Strategic Intelligence Pillar
High-Entropy Alloys (HEAs): Microstructure, Tensile Ductility & Extreme Environment Performance
Explore Topic Pillar

Key Takeaways & Executive Findings

  • • • Maximum DIPSWT of 3.84×10−8 occurs at 38.59 mm from the lower gauge face corner, matching field-observed RCF crack locations between 35–45 mm, enabling targeted inspection and grinding schedules. • • Critical plane orientation of 35–37° relative to rolling direction and 33–44° to rail surface agrees with experimental measurements by Szablewski et al., validating the predictive capability for crack orientation. • • The loaded locomotive configuration yields the highest DIPSWT among all vehicle types, indicating that axle load and traction profile dominate RCF initiation under heavy haul operations. • • The multi-variable sampling technique reduces full loading spectra to representative traction profiles, cutting computational cost while preserving critical damage drivers for efficient RCF assessment.
Weekly Academic Intelligence

China Advanced Materials & Deep-Tech Radar

Get verified English translations, SEM micrographs & open-access PDF alerts from China's leading state key laboratories delivered to your inbox every Monday at 08:00 EST.

Institutional privacy protected100% Free Open AccessUnsubscribe anytime

Abstract

Rolling contact fatigue (RCF) crack initiation in rails remains a critical failure mode under increasing axle loads and line capacities. Existing predictive methods struggle to capture the combined influence of fluctuating creepage, contact positions, and load spectra on crack location and orientation. This study proposes a computational framework integrating multi-body simulations (MBS), finite element analysis (FEA), and critical plane approaches. A novel multi-variable sampling technique reduces complex loading spectra into representative traction profiles, which are then analyzed using the Smith–Watson–Topper damage indicator parameter (DIPSWT). The maximum DIPSWT identifies the critical plane and potential crack orientation. A case study on the Swedish heavy haul line Malmbanan, specifically a 384 m section of a R=450 m curve, incorporates measured traffic and loading conditions. Results show the highest DIPSWT for the locomotive with loaded payload configuration, reaching a cumulative maximum of 3.84×10−8 at 38.59 mm from the lower gauge face corner. The DIPSWT critical plane orientation (35–37° to rolling direction, 33–44° to rail surface) aligns with experimental measurements of RCF cracks near the gauge corner. This method enables efficient identification of RCF-prone conditions and crack orientations, complementing existing predictive tools.

1. Introduction

Rolling contact fatigue (RCF) cracking threatens rail network reliability as axle loads, line capacities, and speeds escalate. Predictive methods have historically been stymied by the multiscale nature of RCF: macroscale load variations from locomotive and wagon sets interact with nanoscale material inhomogeneities, producing crack initiation sites that are difficult to pinpoint. Existing commercial approaches often rely on simplified loading assumptions or empirical thresholds, failing to capture the combined effects of creepage, contact position, and load history on crack orientation and location.

This study addresses the bottleneck by integrating multi-body simulations, finite element analysis, and critical plane approaches with a novel multi-variable sampling technique. The method simplifies complex loading spectra into representative traction profiles, then applies the Smith–Watson–Topper damage indicator parameter (DIPSWT) to identify critical planes and potential crack orientations. A case study on the Swedish heavy haul line Malmbanan, using measured traffic and loading conditions, demonstrates that the predicted DIPSWT hot spots and crack orientations align with experimental observations, offering a computationally efficient tool for RCF risk assessment.

SinoTechIntel Interactive Document Reader
Page 1–5 of Preview
100%
Download Full PDF

Loading authentic research manuscript (Pages 1–5)...

Cite This Research Paper
Jonathan Leung, Saeed Hossein-Nia, Mårten Olsson, Carlos Casanueva (2026). Computational methods to predict RCF crack initiation hot spots in rails using critical plane SWT damage indicator parameter. Railway Engineering Science (铁道工程科学). https://doi.org/10.1007/s40534-025-00405-4
SinoTechIntel Academic & Legal Disclaimer

Research & Educational Purpose Only: The translations, structured abstracts, analytical annotations, and data reports provided by SinoTechIntelare intended exclusively for academic research, internal corporate R&D, and educational benchmarking. They do not constitute formal engineering, chemical safety, legal, or professional advice.

Copyright & Intellectual Property Notice: Original copyright of the underlying source articles and experimental data remains with the respective authors, institutions, and original publishing journals. SinoTechIntel claims intellectual property only over its proprietary translations, analytical syntheses, and AEO structured enhancements in accordance with international fair use and academic citation principles.

Frequently Asked Questions

What is the maximum DIPSWT value and its location, and how does it compare to field observations?

The maximum cumulative DIPSWT is 3.84×10−8, located at 38.59 mm from the lower gauge face corner. This falls within the experimentally observed RCF crack zone of 35–45 mm, confirming the method's accuracy for hot spot prediction.

How does the predicted critical plane orientation compare with experimental measurements?

The DIPSWT critical plane is oriented at 35–37° to the rolling direction and 33–44° to the rail surface, which agrees with experimental measurements of RCF cracks near the gauge corner by Szablewski et al., validating the orientation prediction.

Which vehicle configuration produces the highest RCF damage, and what are the implications for heavy haul operations?

The loaded locomotive configuration yields the highest DIPSWT, indicating that higher axle loads and traction forces drive RCF initiation. This suggests that operational limits should account for locomotive payload and traction profiles to mitigate rail damage.

What computational advantages does the multi-variable sampling technique offer over conventional full-spectrum analysis?

The technique reduces complex loading spectra into representative traction profiles, significantly lowering computational cost while preserving critical damage drivers. This enables efficient parametric studies and integration with existing predictive maintenance frameworks.

What are the limitations of the current study, and what future work is planned?

Limitations include the focus on a single curve radius and track geometry. Future work will explore shear stress-driven crack initiation, incorporate more comprehensive selector variable algorithms, and simulate varied track conditions to enhance generalizability.

Related Chinese Research & Cross-Citations

Research Citation2026
Punching test for mechanical characterization of asphalt railway sub-ballast

Punching test for mechanical characterization of asphalt railway sub-ballast

The adoption of asphalt sub-ballast layers in high-speed and high-capacity railway lines has gained traction across Europe and the United States, driven by structural, functional, and economic advantages. However, current mechanical characterization methods, inherited from road paving practice, fail to capture railway-specific interactions, particularly the ballast/sub-ballast interface and the granular behavior of the overlying unbound layer. This study introduces a novel punching test employing an adaptive indentation plate (AIP) to replicate the contact between ballast particles and the sub-ballast. Cylindrical asphalt specimens (150 mm diameter) were subjected to vertical point loads via the AIP at temperatures of 5, 20, and 35 °C and deformation rates of 5.08, 25.4, and 50.8 mm/min. Two asphalt mixes compliant with Italian sub-ballast standards—one with conventional B50/70 binder (SSB) and one with polymer-modified PmB 45/80-65 (HSB)—were tested to validate the procedure. Force–displacement curves revealed three interaction phases, and key parameters such as peak force (Fmax), displacement at peak force (uy,*), force at 1.5 mm deflection (Fu,1.5), and secant slope (m) between 1.25 and 1.75 mm were extracted. Results demonstrate that the punching test effectively evaluates resistance to plastic deformation and indentation behavior, offering insights beyond existing specifications. The lowest deformation rate (5.08 mm/min) best captures viscoelastic properties. This methodology provides a valuable tool for assessing bituminous sub-ballast performance under simulated railway interface conditions.

Examine Full Data & PDF
Research Citation2026
Dynamic Models and Analysis of Key Factors Influencing Stick–Slip Vibration in Disc Brake Systems

Dynamic Models and Analysis of Key Factors Influencing Stick–Slip Vibration in Disc Brake Systems

Stick–slip vibration in disc brake systems at low speeds arises from the interaction of multiple factors, posing persistent challenges for high-speed train safety and passenger comfort. This study establishes three- and four-degree-of-freedom (DOF) dynamic models that incorporate wheel–rail adhesion and nonlinear friction, validated through line testing. System stability, stick–slip bifurcation characteristics, and key influencing factors are analyzed via numerical simulation. Results demonstrate that the four-DOF model, by accounting for normal motion, avoids the over-evaluation of system stability inherent in the three-DOF model. In the three-DOF model, tangential stiffness is the primary factor inducing chaotic stick–slip vibration. In the four-DOF model, tangential stiffness predominantly affects vibration amplitude, while normal stiffness governs the onset of chaos. Damping exhibits minimal influence on chaotic stick–slip occurrence. Optimal parameter ranges are identified: brake disc rotational inertia of 5–9 kg·m² and 11–22 kg·m², and friction pad mass of 7–17 kg, which effectively mitigate chaotic stick–slip vibration. These findings provide quantitative guidelines for brake system design, enhancing operational reliability and reducing wear-related failures in high-speed rail applications.

Examine Full Data & PDF
Research Citation2026
Acoustic characteristic optimal design for railway steel–concrete composite bridge based on the RBFNN-NSGA-II algorithm

Acoustic characteristic optimal design for railway steel–concrete composite bridge based on the RBFNN-NSGA-II algorithm

Structure-borne noise from railway steel–concrete composite (SCC) bridges presents a persistent challenge due to high sound pressure levels across both low and high frequency ranges. This study establishes a hybrid finite element–statistical energy analysis (FE-SEA) numerical model to predict acoustic radiation from an SCC bridge. Field measurements validate the model with discrepancies of only 0.4 dB and 1.1 dB in overall sound pressure levels. Using uniform design sampling, a high-accuracy radial basis function neural network (RBFNN) surrogate is trained to map cross-sectional parameters to acoustic and cost objectives. The non-dominated sorting genetic algorithm (NSGA-II) then performs multi-objective constrained optimization, generating a Pareto frontier for sound power level (SWL) and material cost. The technique for order preference by similarity to an ideal solution (TOPSIS) selects the optimal parameter combination, achieving a 5 dB reduction in SWL and a 23.9% decrease in material cost. These results demonstrate that strategic cross-sectional adjustments can simultaneously mitigate noise and reduce expenditure, offering a practical framework for acoustic optimization in railway bridge design.

Examine Full Data & PDF
Research Citation2026
A modular simulation tool for fixed block and moving block railway signalling systems

A modular simulation tool for fixed block and moving block railway signalling systems

This paper presents a time-based modular simulation tool implemented in MATLAB® and Simulink® for evaluating fixed block (FB) and moving block (MB) railway signalling systems. The simulator integrates the Radio Block Centre (RBC), On-Board Unit (OBU), vehicle longitudinal dynamics, and a driver behavioural model to replicate realistic train operations. Two case studies are conducted on a main railway line with twelve stations, involving a leading and a trailing train with identical mechanical characteristics. An unexpected departure delay of the leading train at the fifth station reduces the headway, forcing the signalling system to intervene. Simulation A employs identical driver parameters for both trains, while simulation B assigns an aggressive and less responsive driver to the trailing train, resulting in a more pronounced signalling action. Results indicate that MB signalling reduces the trailing train's travel time by 70 s (2.5% of total time) compared to FB, albeit with a more jagged speed profile due to frequent activations. The authors caution that these findings are limited to the two case studies and that a comprehensive parametric investigation is necessary for generalised conclusions. The tool aims to reduce on-site testing of railway signalling systems by providing a flexible platform for scenario analysis.

Examine Full Data & PDF
Research Citation2026
Investigation on the Rolling Contact Fatigue Cracks Initiation of Subway Fixed Frogs Based on Transient Dynamics

Investigation on the Rolling Contact Fatigue Cracks Initiation of Subway Fixed Frogs Based on Transient Dynamics

The escalating traffic density and operational speeds of subway systems have intensified fatigue damage in turnout rails, particularly within the hazardous space of fixed frogs where wheel–rail dynamic interaction is exacerbated. This study addresses the rolling contact fatigue (RCF) crack initiation behavior of a No. 9 turnout fixed frog, a configuration widely deployed in subway networks. A three-dimensional explicit transient rolling contact finite element model was developed to simulate wheel–rail interaction under varying vehicle speeds and fastener vertical stiffness conditions. The analysis focused on crack initiation locations, angles, and fatigue life. Results demonstrate that the 30 mm top width cross-section of the nose rail is the most susceptible to fatigue cracking, with cracks initiating on the rail surface. The angle between the crack initiation surface and the lateral direction ranges from 70° to 95°, while the angle relative to the vertical direction remains difficult to predict. Higher vehicle speeds significantly reduce fatigue life, whereas fastener vertical stiffness exerts a minor influence. The calculated RCF crack initiation life is approximately 24,000 cycles across three stiffness conditions. Simulation outcomes align with field survey findings, validating the model's fidelity. The established methodology provides theoretical support for optimizing fixed frog structures and predicting fatigue life in subway turnouts.

Examine Full Data & PDF
Research Citation2026
Performance-based comprehensive functional damage probability assessment framework for high-speed railway bridge under earthquake

Performance-based comprehensive functional damage probability assessment framework for high-speed railway bridge under earthquake

Current seismic damage assessments for high-speed railway (HSR) bridges primarily focus on overall structural safety, lacking evaluations from multiple performance perspectives, which affects post-earthquake traffic decision-making. This study proposes a performance-based comprehensive functional damage probability assessment framework for high-speed railway simply supported bridges (HSRSSBs) under earthquakes. The framework categorizes bridge functions into three levels: post-earthquake traffic function (PTF), structural bearing function (SBF), and collapse resistance function (CRF), corresponding to operational, structural safety, and structural integrity requirements, respectively. By analyzing damage states of key bridge components during earthquakes, functional damage probability assessment indicators and classification thresholds are established according to various performance requirements. Damage probability calculations are conducted using the probability density evolution method and vulnerability method. Based on the relationship between damage probabilities at different functional levels, a comprehensive damage probability assessment framework considering the three-level performance requirements of HSRSSBs is developed, and the influence of varying pier heights on the functional damage probability relationship is examined. Results indicate that current HSRSSB designs meet all performance requirements under frequent earthquakes. Under design-level earthquake conditions, the SBF remains in a slight damage state, while the PTF exhibits varying degrees of damage, which worsens as pier height increases. The pier structure satisfies seismic demands even under rare earthquake conditions.

Examine Full Data & PDF