This work investigates the numerical simulation of the sub-interface crack in the limited permeable piezoelectric bi-material. The mechanical and dielectric equilibrium is satisfied by incorporating electrostatic tractions on the crack surface. A crack opening model is used to incorporate the crack face tractions in the framework of the extended finite element method (XFEM). The level-set method is used in association with XFEM to track the discontinuity in the discretized domain. The sixfold electro-mechanical enrichment functions are adopted to get the effect of arbitrary poling direction in the context of XFEM. The electric displacement intensity factor (EDIF) is obtained using the electro-mechanical interaction integral in conjunction with asymptotic crack tip variables. The effect of poling direction and permittivity of the medium are studied to understand the behavior of EDIF.