Numerical modeling of sub-surface initiated spalling in rolling contacts

A 3D finite element model was developed to investigate the influence of microstructure topology on the stochastic nature of rolling contact fatigue. Grains of the material microstructure are modeled with random Voronoi tessellations. Continuum damage mechanics and mesh partitioning are implemented to capture the initiation and propagation phases of fatigue damage that lead to spalling. Simulated fatigue spalling is shown to progress similarly to experimental observations of rolling contact fatigue. The fatigue lives obtained with the model exhibit scatter on par with empirical measures and are fit well by 2 and 3-parameter Weibull distributions.