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.　©　2016,　Beijing　University　of　Technology.　All　right　reserved.