Currently, few studies are focused on the stationary random vibration for composite laminated shell structures of revolution (CLSSR), including composite laminated cylindrical shell (CLCY), composite laminated conical shells (CLCO), and composite laminated annular plates (CLAP). To fill this void corresponding to the above research in the literatures, a combination of the spectro-geometric method (SGM) and pseudo-excitation method (PEM) was developed to construct the theoretical model within the first-order shear deformation theory (FSDT). The different boundary restraints and coupling conditions were achieved by taking the appropriate stiffness values of artificial springs, and the thermal effect induced by thermal load was considered. Moreover, the Rayleigh–Ritz method was employed to deduce the governing differential equation. Further, the solution accuracy of the established model was assessed by comparing the obtained results with those from the literatures and the finite element method (FEM). Finally, the effect of specific parameters (i.e. fiber angle, temperature value and ply number) on the stationary random response of CLSSR was explored. According to the results, the proposed method proved effective for predicting the stationary random response characteristics of CLCY, CLCO, and CLAP in a thermal environment.