In case of rough contact fatigue, the accurate calculation of surface tractions is essential to the prediction of crack initiation life. Accurate Surface tractions influencing shear stress amplitude can be obtained by contact analysis based on the morphology of contact surfaces. In this study, to simulate rough contact under sliding condition, gaussian rough surface generated numerically in the previous study was used and to calculate clack initiation life in the substrate, dislocation pileup theory was used.