The real working condition of the roller bearing is considered by quasi-dynamic method. Voronoi finite element method is adopted to simulate the microstructure of bearing steel. Stress distribution in micro-domain with an inclusion is obtained to determine the location of initial crack. The effect of inclusion elastic modulus, size and depth on crack propagation is analyzed.