In this study, the simulation of contact fatigue based on stress analysis is conducted under Elastohydrodynamic Lubrication (EHL) state. To predict a crack initiation life accurately, it is necessary to calculate contact stress and subsurface stresses accurately. Contact stresses are obtained by contact analysis of a semi-infinite solid based on the use of influence functions and the subsurface stress is obtained using rectangular patch solutions. The numerical algorithm using newton-rapson method was constructed to calculate the EHL pressure. Based on these stress values, three multiaxial high-cycle fatigue criteria are used. As a result, the effects of EHL on contact fatigue life are calculated.