Functionally graded piezoelectric (FGP) material removes the brittleness issue of piezoceramics and small coupling coefficient of piezoelectric polymers. They also avoid the stress concentration and consequently crack propagation. This work addresses the finite element modeling of the functionally graded piezoelectric harvesters subjected to random vibrations. The finite element electro-mechanical governing equations of FGP beam are derived by using the generalized Hamilton’s principle under the assumptions of Euler–Bernoulli beam theory. The material properties of the beam are assumed to be varied in the thickness direction following a simple power law distribution. The random base excitation of FGP harvester is assumed to be the Gaussian white noise signal. Both time domain and frequency domain analysis are presented in this work.