Subsurface crack mode II propagation parallel to the contact surface is a damage mechanism leading to dramatic failure in many components subjected to cyclic loading. A weight function (WF) was elaborated for calculating the applied mode II stress intensity factor (SIF) of a crack in a two‐dimensional half‐space in plane strain condition, for crack completely closed and frictionless contact between the crack faces. With respect to other methods, the WF allows faster SIF calculation, thus being suitable for simulation of many repeated load cycles and fatigue crack propagation. The WF was applied for simulating a case of rolling contact experiments found in the literature, and good agreement between experimental and numerical results was obtained, showing the effectiveness of the WF method in damage tolerant design.