Abstract. The available three-dimensional (3D) solutions reported for the analysis of Functionally Graded Piezoelectric (FGP) laminates are based on this simplifying assumption that the FGP layer consists of a number of homogeneous sub-layers. The accuracy of these formulations not only is dependent on the number of sub-layers, but also leads to inaccurate results in the prediction of higher natural frequencies of the FGP laminates. In the present paper, a 3D Peano series solution is developed for the cylindrical bending vibration of the FGP laminates. This novel formulation exactly satis es the equations of motion, the charge equation, and the boundary and interface conditions of the continuously nonhomogeneous piezoelectric layers. The obtained solution is exact because no a priori assumption for the displacement components and the electric potential along the thickness direction of FGP layers is introduced. The in uences of the di erent functionally gradient material properties and di erent electric boundary conditions on the natural frequencies and mode shapes of the FGP laminate have also been studied through examples. The present solution and its obtained numerical results can be employed to assess the accuracy of di erent FGP laminated beam/plate theories. It can also be used for FGP vibration behavior comprehension purposes.