Based on the theory of continuous damage mechanics, this paper proposes a fatigue damage evolution model and numerical calculation method considering the influence of vacuum chemical heat treatment process, and studies the rotational bending fatigue damage analysis and life prediction method of M50NiL bearing steel. Firstly, the constitutive model, fatigue damage evolution equation and material parameter calibration method in the theoretical model are given. Secondly, based on the ABAQUS platform, the damage mechanics-finite element numerical calculation method of fatigue damage analysis considering the influence of hardened layer is realized by writing the UMAT subroutine; then, the rotational bending fatigue test of M50NiL bearing steel treated by vacuum chemical heat treatment of quenching and tempering, carburizing and carbonitriding is carried out, and the influence of heat treatment process on fatigue performance is analyzed; finally, based on the proposed fatigue damage model and numerical calculation method, the rotational bending fatigue life of M50NiL bearing steel is predicted, and the results are compared with the experimental results to verify the correctness of the proposed method.