An analysis is presented for the free vibration of angle-ply laminated circular cylindrical thick shells having any combinations of boundary conditions. For this purpose, the maximum strain and kinetic energies of the shell are evaluated analytically according to the first-order shear deformation theory (FSDT), and the frequency equation is derived by minimizing the Lagrange functional with respect to unknown coefficients of the displacement functions. These displacement functions are written in the power series which satisfy any combination of the geometrical boundary conditions automatically. In numerical examples, the natural frequencies and the mode shapes of the angle-ply laminated thick circular cylindrical shells are calculated, and the effects of the lamination parameters (fiber orientation, number of layers) and the boundary conditions on the Vibration are studied.