SUMMARYThe extended finite element method (XFEM) is further improved for fracture analysis of composite laminates containing interlaminar delaminations. New set of bimaterial orthotropic enrichment functions are developed and utilized in XFEM analysis of linear-elastic fracture mechanics of layered composites. Interlaminar crack-tip enrichment functions are derived from analytical asymptotic displacement fields around a traction-free interfacial crack. Also, heaviside and weak discontinuity enrichment functions are utilized in modeling discontinuous fields across interface cracks and bimaterial weak discontinuities, respectively. In this procedure, elements containing a crack-tip or strong/weak discontinuities are not required to conform to those geometries. In addition, the same mesh can be used to analyze different interlaminar cracks or delamination propagation. The domain interaction integral approach is also adopted in order to numerically evaluate the mixed-mode stress intensity factors. A number of benchmark tests are simulated to assess the performance of the proposed approach and the results are compared with available reference results.