Abstract: Based on the general linear elastic theory of thin shells and the potential flow theory, considering the energy dissipation of passive constrained layer damping (PCLD) due to shear deformation of the viscoelastic core and the fluid-solid coupled interaction, a first order integrated differential matrix equation of the circular cylindrical shell with PCLD treatment was derived firstly. Each element in the state vector of the equation has clearly physical means, so it can be conveniently applied to solve dynamic problems of circular cylindrical shell with different boundary conditions. Secondly, by writing the fluid dynamic pressure as analytical form with undefined coefficients, and applying the fluid-solid interaction condition, the homogenous extended capacity precision integration method and the superposition principle, a semi-analytical and semi-numerical method with high efficient and precision is established. Comparing to the analytical solutions of the axissymmetrical free vibration of cylindrical shell (no water filled) fully treated with PCLD under the freely supported, the results obtained by this method is effective. Finally, the dynamic responses of the fluid-filled cylindrical shells with partially treated PCLD under harmonic ground motion are represented, and the damping effects of different thicknesses, lengths, and locations of PCLD and the complex shear modulus of the viscoelastic core are studied.