This paper proposes a new strategy based on the proper generalized decomposition method to solve the forced vibration problem for bi-dimensional piezoelectric composite beams. The approach considers a classical harmonic space-frequency description of the dynamic problem and redefines it by approximating the displacement and electric potential fields as a sum of separated functions of x (beam axis coordinate), z (thickness coordinate) and (load frequency). The methodology consists in an iterative algorithm that solves three one-dimensional linear problems at each iteration. The result is a 2D solution in frequency domain with 1D computational complexity. Several numerical tests with wide range of slenderness ratios are considered in order to assess the validity of the method. Moreover, a study of different combinations of boundary conditions is carried out. The results in terms of modal parameters and frequency response functions are analysed and compared with exact elasticity solution and finite element simulations.