In case of contact fatigue, the accurate calculation of surface tractions and subsurface stress is essential to the predication of crack initiation life. Surface tractions influencing shear stress amplitude have been obtained by contact analysis based on influence function. Subsurface stress has been obtained by using rectangular patch solutions. In this study, to simulate asperity contact under sliding condition, the tip of asperity was simulated by sphere and to calculate crack initiation life in the substrate, dislocation pileup theory was used.