In this study, using high cycle fatigue (HCF) criteria, the simulation of rolling contact fatigue is conducted under elliptical contact. The HCF criteria fall into three categories: the critical plane approach, the stress invariant approach and the approach based on the mesoscopic scale. The accurate calculation of contact stresses and subsurface stresses is essential to the prediction of crack initiation life. Contact stresses are obtained by contact analysis of a semi-infinite solid based on the use of influence functions and the subsurface stress field is obtained using rectangular patch solutions. The simulation results show that the critical load is decreasing rapidly and the site of crack initiation also moves rapidly to the surface from the subsurface when the friction coefficient exceeds a specific value for all of three fatigue criteria.