The problem of parametric oscillations of an isotropic viscoelastic shallow shell of variable thickness under periodic load is considered. It is believed that under the influence of specified load, the shallow shell allows displacements (in particular, deflections), commensurate with its thickness. In a geometrically nonlinear statement, taking into account the viscoelastic properties of material, a mathematical model of the problem has been developed using the classical Kirchhoff-Love hypothesis. Using the Bubnov-Galerkin method based on the polynomial approximation of the deflections, the problem is reduced to the study of the system of integro-differential equations, where time is the independent variable. The solution of the system of integrodifferential equations is determined by the proposed numerical method. Based on this method, a numerical solution algorithm is described. The Koltunov-Rzhanitsyn kernel with three different rheological parameters is chosen as a weakly singular kernel. At the same time, the effect of geometric nonlinearity, viscoelastic properties of material, as well as other physicomechanical and geometric parameters and factors (rheological parameters, thickness, initial shape imperfections, aspect ratios, boundary conditions, excitation coefficient) on the area of dynamic instability is taken into account. The results obtained in this study are in good agreement with the results and data obtained by other authors.