This paper presents domain form of the interaction integrals based on three independent formulations for computation of the stress intensity factors and electric displacement intensity factor for cracks in functionally graded piezoelectric materials. Conservation integrals of J-type are derived based on the governing equations for piezoelectric media and the crack tip asymptotic fields of homogeneous piezoelectric medium as auxiliary fields. Each of the formulation differs in the way auxiliary fields are imposed in the evaluation of interaction integral and each of them results in a consistent form of the interaction integral in the sense that extra terms naturally appears in their derivation to compensate for the difference in the chosen crack tip asymptotic fields of homogeneous and functionally graded piezoelectric medium. The additional terms play an important role of ensuring domain independence of the presented interaction integrals. Comparison of the numerically evaluated intensity factors through the three consistent formulations with those obtained using displacement extrapolation method is presented by means of two examples. © 2008 Elsevier Ltd. All rights reserved.