Abstract:The theoretical research history of cylindrical thick shells’ stress and displacement is firstly reviewed, and a 3D mechanical model of length-limited cylindrical thick shells under non-axisymmetric moving loads is established. Then based on the assumption of radial shear strains with quadratic distribution and radial normal strains with linear distribution through the radial coordinate, displacement expressions satisfying boundary conditions are established. The non-axisymmetric loads are expressed by Heaviside function and Dirac function, and equilibrium differential equations containting unkowns are derived according to the minimum potential energy principle. After application of the Galerkin Method and the Modified Runge–Kutta–Fehlberg Mehtod, dynamic response of cylindrical thick shells under non-axisymmetric moving loads is obtained. This method is verified by a comparative analysis of the theoretical solutions with ANSYS numerical results from an example.