Motivated by experimental observations, we carry out a numerical analysis of the two-stage crack growth under fretting fatigue by using an efficient and accurate boundary element method. To start with, the variation of stress field during a loading cycle is analyzed. Various values of friction coefficient in the contact zone are considered, which is shown to considerably affect the stress field. Then, by assuming crack initiation to occur in the shear mode, a surface-breaking crack is introduced to the specimen at the location of highest shear-stress amplitude. The crack-tip stress intensity factors (SIFs) are calculated for various crack lengths and at various crack angles ranging from 25Β°to 45Β°about the contact surface. It is shown that, for a loading ratio of 0.5, the cyclic mode-II SIF amplitude decreases with increasing crack length, whilst its mean value increases. It suggests that the (first-stage) shear crack would sooner or later become dormant, or switch to another mode that can provide continuous support of growth. Then, the first-stage shear crack is manually kinked into a second-stage opening crack, and the follow-on driving force is analyzed. It is shown that the kinking event is only favored after the first-stage crack has grown to a certain length. The present study thus provides insights in the mechanics of twostage crack growth that has been frequently observed in a typical dovetail joint under fretting fatigue. It also suggests an improved experimental setup to quantitatively investigate the fretting fatigue in dovetail joints.