This paper presents a semi-analytical approach based on the shifted Legendre series for vibration modeling of complex axisymmetric shells. This method divides the structure along its generator into complex axisymmetric shells comprised solely of conical and spherical shell segments. A global cylindrical coordinate system is established to unify the coordinate system of complex axisymmetric shells. Coupling between segments and the simulation of arbitrary boundary conditions are achieved using artificial spring technology. Based on the first-order shear deformation theory (FSDT), the energy functional of each segment is derived. The circumferential displacement is expanded using Fourier series, while the meridional displacement is expanded using the shifted Legendre series at shifted Gauss-Lobatto nodes. Utilizing the inner product matrix and differential matrix of the shifted Legendre series, the integration and differentiation processes in the vibration characteristic equation are simplified into a product form, significantly enhancing computational efficiency. Through convergence analysis and comparison of the computational results of the present method with numerical simulation results, experimental results and literature results demonstrate the advantages of high convergence, high accuracy and fast solution speed of this solution method. On this basis, the vibration characteristics of ring-stiffened double-layer hemispherical shells are investigated. The results demonstrate the great applicability of this method for vibration modeling of complex axisymmetric shells in engineering.