Abstract: To get a more effective and accurate method to study the fracture behavior of the piezoelectric materials,a computational formula based on the interaction integral method to calculate the crack-tip field intensity factors was obtained in this paper.This method was successfully applied to the boundary element analysis software.With this software,the fracture behavior of PZT-4 made piezoelectric board's central crack and edge crack under different loading conditions was analyzed in different polarization directions.The path-independent characteristics of the interaction integral method was verified.The crack-tip field intensity factors obtained were compared well with the results of both the J integral method and the displacement interpolating method.This proves that this method is effective and accurate.