Piezoelectric materials are used in a wide range of engineering applications as sensors, actuators, and transducers due to their intrinsic properties. However, analyzing these devices is not straightforward because of the material anisotropy and coupled electromechanical nature of the problems. In this study, a stabilized non-ordinary state-based peridynamic model for linear piezoelectricity is proposed. Fully coupled electromechanical simulation is conducted using an implicit method, and the tangent stiffness matrix is obtained via perturbation technique. Three two-dimensional electromechanical problems with combined static loading and various boundary conditions are used to validate the proposed technique. Then, the crack-opening displacements, and induced electric field on the crack-face are presented for a center-cracked PZT-5H plate. Hence, this study demonstrates that the proposed technique is capable of modeling linear piezoelectric materials and provides a promising method to study fracture and damage progression in piezoelectric materials.