In this chapter, a new type of Computational Grains is proposed to study the micro-electro-mechanical behavior of composite piezoelectric materials. This method is based on a new hybrid variational principle, and independently assumed displacements and electric potentials in each CG. Each CG can efficiently model a single physical grain of the composite material, thus saving a significant time of generating complex FEM mesh. Computational Grains can also model porous and composite piezoelectric materials even if the distribution of voids/inclusions is not symmetrical (which is assumption used with all unit cell models). Because the trial solutions are complete but do not satisfy the governing differential equations a priori, the formulation is very simple, and can accurately account for the local field concentrations efficiently and accurately. This is illustrated using different examples where the fields along the void/inclusion periphery are calculated, and the effective material properties of porous and composite materials are predicted, and compared with other analytical and computational models. The proposed CG in this chapter is expected to become a very powerful tool of direct numerical simulations of the micro/meso mechanics of composite piezoelectric materials, and can possibly lead to efficient multi-scale modeling of piezoelectric devices.