This paper presents an efficient meshless method for analyzing cracked piezoelectric structures subjected to mechanical and electrical loading. The method employs an element free Galerkin (EFG) formulation and an enriched basic function as well as special shape functions that contain discontinuous derivatives. Based on the moving least squares (MLS) interpolation approach, The EFG method is one of the promising methods for dealing with problems involving progressive crack growth. Since the method is meshless and no element connectivity data are needed, the burdensome remeshing procedure required in the conventional finite element method (FEM) is avoided. The numerical results show that the proposed method yields an accurate near-tip stress field in an infinite piezoelectric plate containing an interior hole. Another example is to study a ceramic multilayer actuator. The proposed model was found to be accurate in the simulation of stress and electric field concentrations due to the abrupt end of an internal electrode.