Mechanical transmission components are generally under the risk of contact fatigue, resulting in cracks initiated from the subsurface. When subject to compressive stress, these cracks might be forced to be closed, thus their faces are in partial slip contact. Under such conditions, the interactions of stick and slip zones cannot be directly determined, which is one of the greatest challenges of contact fatigue problems. In this study, the distributed dislocation technique is adopted to investigate the stick–slip friction along the alignment of shear cracks in arbitrary orientations in a half-plane under contact loading. The unknown dislocation density describing the crack face displacement is iteratively calculated based on the stress conditions and then introduced to update the surface gap between contacting bodies. The complex expressions of the influence coefficients used to determine the stress and displacement induced by the dislocations are derived and the FORTRAN codes are provided in the Supplementary Materials. A step-by-step loading method is proposed to progressively capture the slip history of crack faces. The model is also capable of simulating rolling contact problems. The results are comparable to those predicted using the finite element method with more robust convergence properties. The effects of friction coefficient, crack alignment, and locations of crack faces are discussed in detail. The model paves the foundation for further studies on competing mechanisms of crack propagation and crack face wear, and provides insights into the frictional fracture of materials under contact loading.