A novel integral equation method is developed in this paper for the analysis of two-dimensional general piezoelectric cracked bodies. In contrast to the conventional boundary integral methods based on reciprocal work theorem, the present method is derived from Stroh’s formalism for anisotropic elasticity in conjunction with Cauchy’s integral formula. The proposed boundary integral equations contain generalized boundary displacement (displacements and electric potential) gradients and generalized tractions (tractions and electric displacement) on the non-crack boundary, and the generalized dislocations on the crack lines. The boundary integral equations can be solved using Gaussian-type integration formulas without dividing the boundary into discrete elements. The crack-tip singularity is explicitly incorporated and the generalized intensity factors can be computed directly. Numerical examples of generalized stress intensity factors are given to illustrate the effectiveness and accuracy of the present method.