In this study, the simulation of rolling contact fatigue based on stress analysis is conducted under Elastohydrodynamic Lubrication 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-infinile solid based on the use of influence functions and the subsurface stress field is obtained using rectangular patch solutions. And a numerical algorithm using Newton-Rapson method was constructed to calculate the Elastohydrodynamic lubrication pressure. Based on these stress values, several multiaxial high-cycle fatigue criteria are used and the critical loads corresponding to fatigue limits are calculated.