AbstractIn this study, the dislocation‐based model is developed to study the interaction between time‐harmonic elastic waves and multiple interface cracks in 2 bonded dissimilar piezoelectric layers. In this model, cracks are represented by a distribution of so‐called electro‐elastic dislocations whose density is to be determined by satisfying the boundary conditions. Using the Fourier transform, this formulation leads to 2 singular integral equations, which can be solved numerically for the densities of electro‐elastic dislocations on a crack surface. The formulation is used to determine dynamic field intensity factors for multiple interface cracks without limitation of number of cracks. The dynamic field intensity factors are then calculated for both permeable and impermeable crack, and finally, numerical results are presented to illustrate the variation of these quantities with the electromechanical coupling, crack spacing, and the frequency of loading.