The purpose of this study is to present a new method for estimating the fatigue life of the screw blade in the screw sand washing machine. To ensure the accuracy of numerical simulation, the loading area and the value of load are determined by means of the theoretical analysis. To ascertain the location of the stress peak and stress range, the static analysis of the screw blade is executed via the finite element method. To reduce the research cost and ensure the feasibility of the research method, Markov chain Monte Carlo (MCMC) is employed to simulate the random load with the Gauss distribution on the screw blade. In addition, the Nondominated Sorting Genetic Algorithm (NSGA-II) is utilized to find out an optimum variation coefficient of the stress, aiming at guaranteeing the precision of the random load. The rainflow cycle extrapolation is adopted to generate the fatigue load spectrum closer to the real condition, taking account of the possibility of the extreme loads caused by overload occurrence. Subsequently, the rainflow matrix after extrapolation, the estimated P-S-N curve, Goodman stress correction method and Miner's rules are made use of assessing the service life of the screw blade. In particular, the effects of the surface roughness, residual stresses and fatigue notch factors on the fatigue life are taken into consideration. Ultimately, the non-linear surface fitting technique is used to obtain the equation concerning the fatigue life of the screw blade versus residual stresses and fatigue notch factors. The numerical results show that the stress peak is in the root of the screw blade and the service life of the screw blade declines exponentially with growing residual stresses and fatigue notch factors.