This paper derives a semi-analytical solution to determine displacements and stresses in a thick cylindrical shell with variable thickness under uniform pressure based on disk form multilayers. The proposed study partitions the thick cylinder into disk-layer parts based on their thickness of the cylinder. According to the existence of shear stress in the thick cylindrical shell with variable thickness, the equations governing disk layers are acquired based on first shear deformation theory (FSDT), which are in the form of a set of general differential equations. In this study, the cylinder is partitioned into n different disks and n sets of differential equations are derived. The solution of these equations provides displacements and stresses based on the boundary conditions and continuity conditions between the layers. The results are compared with those obtained through the analytical solution and the numerical solution. For the purpose of the analytical solution, matched asymptotic method (MAM) and for the analytical solution, the finite element method (FEM are implemented.