To improve computational accuracy and efficiency of the energy release rate for piezoelectric body containing a crack and to reduce the workload of program compiling and debugging, a formulation based on the four node piezoelectric element was presented for electromechanical virtual crack closure technique. Based on the proposed formulation, a fracture of piezoelectric element was implemented with ABAQUS user defined subroutine UEL with the idea of dummy nodes. The element can be solved independently with energy release rate component, which is validated by comparing the analytical solution with the theory of a center cracked piezoelectric with different materials, crack length and mesh for energy release rates. Numerical examples show that this new fracture piezoelectric element is easy to use with high precision and is insensitive to the mesh size.