Elastodynamic Green's functions for a piezoelectric structure represent the electro-mechanical response due to a steady-state point source as either a unit force or a unit charge. Herein, Green's functions for a laminated circular piezoelectric cylinder are constructed by means of the superposition of modal data from the spectral decomposition of the operator of the equations governing its dynamic behavior. These governing equations are based on a semi-analytical finite element formulation where the discretization occurs through the cylinder's thickness. Examples of a homogeneous PZT-4 cylinder and a two-layer cylinder composed of a PZT-4 material at crystal orientations of +/-30degrees with the longitudinal axis are presented. Numerical implementation details for these two circular cylinders show the convergence and accuracy of these Green's functions. (C) 2004 Elsevier Ltd. All rights reserved.