A numerical method based on the finite fracture mechanics model was proposed in this study to accurately predict the failure of irregular arranged multi-bolt composite repair. The stress distribution and the stress intensity factors required in the model were obtained using a two-stage analysis strategy, which can effectively lower the FE model scale and computational cost while keeping high computational accuracy. In the two-stage analysis strategy, the bolt load distribution was calculated using the global analysis method. After obtaining the fastener load from the first-level model, detailed stress and stress intensity factor analysis were conducted around individual fasteners and cutouts, which were then used to evaluate the failure behaviour of the composite panels. To validate the numerical method, tensile tests were carried out on three kinds of bolted composite repairs with circular, elliptical and oblong cutouts. Failure mode, failure location and failure strength predicted by the proposed numerical method were consistent with the experimental results, illustrating the effectiveness and applicability of the numerical method on predicting the failure behaviour of bolted composite repairs.