Abstract: Cohesive zone model(CZM) was applied to numerical analysis of wavelet finite element method(WFEM) by using the excellent localization properties of interval B-spline wavelet, interface elements of wavelet cohesive zone were constructed by regarding interval B-spline wavelet scaling functions as the interpolation functions, stiffness matrix of interface elements of wavelet cohesive zone was derived, and strain energy release rate(SERR) of interface crack was calculated by using virtual crack closure technique(VCCT), and the quasi-static analysis of interface crack growth was achieved by taking β-Κ fracture criterion. The numerical analyzed results of SERR by WFEM and conventional finite element method(CFEM) were compared with theoretical results, and the results show that the SERR calculated by WFEM and CFEM are 96.60 J/m² and 101.43 J/m², respectively, and the relative errors of SERR between numerical and theoretical results of two method are 1.85% and 3.06%, respectively. It clearly indicates that WFEM can obtain higher accuracy and efficiency of calculating interface crack growth with less numbers of elements and nodes. On this basis, we studied the influence of initial interface crack length and bi-material elastic modulus ratio on interface crack growth. It is found that the equivalent stress at interface crack tip tends to increase with the increase of initial interface crack length. The difference of bi-material elastic modulus ratio is greater, interface crack is easier to expand, and the length of interface crack growth is also greater, therefore, interface crack growth can be slowed by adjusting bi-material elastic modulus ratio.