Analytical model to calculate the contact fatigue life of rough surface is presented in this paper. The effect of surface roughness can be calculated by this model. Computational method and the theoretical basis are also discussed. Contact stresses are obtained by contact analysis of a semi-infinite solid based on the use of influence functions; the subsurface stress field is obtained using rectangular patch solutions. Mesoscopic multiaxial fatigue criterion which can yield satisfactory results for non-proportional loading is then applied to predict fatigue damage. Suitable counting method and damage rule were used to calculate the fatigue life of random loading caused by rough surface. As a result of analysis the relationship between the life and the roughness as well as the most probable depth of the crack initiation is calculated.