A simple computational model for modeling of subsurface crack growth under cyclic contact loading is presented. In this model, it is assumed that the initial fatigue crack will initiate in the region of the maximum equivalent stress at certain depth under the contacting surface. The position and magnitude of the maximum equivalent stress are determined by using the equivalent contact model, which is based on the Hertzian contact conditions with frictional forces. The virtual crack extension method is used for simulation of the fatigue crack growth from the initial crack up to the formation of the surface pit due to contact fatigue. The relationships between the stress intensity factor and crack length are then determined for various combinations of equivalent contact radii and loadings.