A method is presented to determine the forced responses of piezoelectric cylinders using weighted sums of only certain exact solutions to the equations of motion and the Gauss electrostatic conditions. One infinite set of solutions is chosen such that each field variable is expressed in terms of Bessel functions that form a complete set in the radial direction. Another infinite set of solutions is chosen such that each field variable is expressed in terms of trigonometric functions that form a complete set in the axial direction. Another solution is used to account for the electric field that can exist even when there is no vibration. The weights are determined by using the orthogonal properties of the functions and are used to satisfy specified, arbitrary, axisymmetric boundary conditions on all the surfaces. Special cases including simultaneous mechanical and electrical excitation of cylinders are presented. All numerical results are in excellent agreement with those obtained using the finite element software ATILA. For example, the five lowest frequencies at which the conductance and susceptance of a stress-free cylinder, of length 10 mm and radius 5 mm, reach a local maximum or minimum differ by less than 0.01% from those computed using ATILA. © 2005 Acoustical Society of America.