In this study, a new adhesive contact model is built upon a boundary element method (BEM) model developed by Pohrt and Popov (2015). The strain energy release rate (SERR) on the edge of the bonding interface is evaluated using Virtual Crack Closure Technique (VCCT) which is shown to have better accuracy and weaker mesh-size dependency than the closed-form SERR formula derived by Pohrt and Popov. A composite delamination criterion is proposed for crack nucleation and propagation. Numerical results predicted by the present model are in good agreement with the analytical solutions of two classic problems, namely, the axisymmetric parabolic contact and the sinusoidal waviness contact in the plane strain condition. The model of Pohrt and Popov can achieve a similar accuracy for the axisymmetric parabolic contact where the mesh grid is non-conforming to the crack front. Once the conforming mesh grid is used, the accuracy of their model is significantly deteriorated, especially at high work of adhesion and high mesh density. In both BEM models, however, the crack nucleation is found to be mesh-dependent which may be solved by introducing an upper limit for the tensile normal traction.