Subsurface rolling contact fatigue (RCF) failure is one of the primary failure modes in properly installed and lubricated rolling bearings. Its actual service life often exhibits significant scatter, posing a formidable challenge to the reliable life prediction and operational safety of bearings. This study establishes a macro-meso-coupled rolling contact fatigue model that accounts for crystalline anisotropy and grain topological structures. This model utilizes Voronoi tessellations and random Euler angles to construct a polycrystalline mesoscopic model, which is subsequently integrated with a macroscopic Hertzian contact finite element analysis to simulate the roller bearing loading cycles and determine the localized stress responses within the material. The results indicate that variations in macroscopic structural and operating parameters primarily affect the overall stress level of the subsurface RCF failure. The relative fatigue life of the bearing exhibits an exceptionally high sensitivity to changes in macroscopic and operating parameters. Specifically, an increase in radial load leads to an exponential decrease in relative life, with the Weibull slope ranging between 1.001 and 1.129, which is broadly consistent with the classical Lundberg–Palmgren experimental value of 1.125. Conversely, the heterogeneity of the mesoscopic crystalline structure strongly influences the statistical variance of localized extreme stresses. The scatter in bearing fatigue life demonstrates a much more pronounced sensitivity to mesostructural alterations; as the grain size increases from 10 μm to 40 μm, the Weibull slope drops from 1.041 to 0.784. This study provides an analytical basis for the reliable life prediction of rolling bearings.
Liu et al. (Sat,) studied this question.