In this paper, Extended Finite Element Method is employed to solve problems, including weak and strong discontinuities. To this end, a level set framework is used to represent the discontinuities location and Heaviside and Branch function are included in standard finite element method. The case of two arbitrary curved cracks is solved numerically and SIF values at the crack tips are calculated based on the evaluation of the Crack Tip Opening Displacement (CTOD). Afterwards, J-integral methodology is adopted to evaluate the SIFs for isotropic and anisotropic bimaterial interface crack problems. Numerical results are verified with those presented in the literature.