The symplectic method is introduced for boundary-condition problems of finite viscoelastic cylinders. On the basis of the state space formalism and the use of the Laplace integral transform, the general solution of the governing equations, zero- and nonzero-eigenvalue eigenvectors, are obtained. Since the eigenvectors are expressed in concise analytical forms, the adjoint symplectic relation of the Laplace domain is generalized to the time domain. Therefore, the particular solution and the eigenvector expansion method can be discussed directly in the eigenvector space of the time domain, without employing the iterative application of the inverse Laplace transformation. Using this method, various boundary conditions, the particular solution of nonhomogeneous equations, especially the interfacial continuity conditions of composite materials, can be conveniently described by combinations of the eigenvectors.