In the present article, the idea of using the variable-order fractional-derivative thermoviscoelastic constitutive laws in dynamic stress and vibration analysis of the engineering structures, the required implementation backgrounds, and the relevant numerical solution procedures are investigated for the first time. In this regard, dynamic 3D stress and displacement fields and radial/transverse vibrations of transversely graded viscoelastic spinning thick plates/discs exposed to sudden thermoelastic loads are investigated. Instead of using the approximate plate theories, the exact thermoviscoelasticity theory is employed in the development of the governing equations. Since the variable fractional order is dependent on the localized deformation rates, the resulting thermoviscoelastic integro-differential equations are nonlinear. These equations are solved by utilizing a combination of the second-order backward/central/forward finite difference discretization of the spatial and time domains, numerical evaluation and updating of the Caputo-type fractional derivatives, updating the growing number of terms of the governing equations, and Picard's iterations. Various edge conditions are considered. Finally, comprehensive sensitivity analyses and various 3D plots are presented and discussed regarding the effects of the variable fractional order of the constitutive law, time variations of the nonuniformly distributed transverse loads, and edge conditions on the distributions and damping of the resulting displacement and stress components.