PurposeConsidering the size-dependent influence ignored by classical continuum mechanics, a new non-classical Euler–Bernoulli beam model is proposed in this paper. The new fractional viscoelastic nanobeam model is set up using the fractional Kelvin–Voigt viscoelastic model and Hamilton’s principle. And the new model studies the total effects of nonlocal elasticity, modified couple stress, and surface energy.MethodsThe model represented the fractional integral-partial differential governing equation is solved by Galerkin’s and predictor–corrector methods. First, the effects of nonlocal elasticity, modified couple stress, surface energy, and their coupling impact on the nonlinear time response of free vibration of fractional viscoelastic nanobeam are analyzed. Then, in the frame of nonlocal couple-stress elasticity and surface energy theory, the effects of different parameters on the nonlinear time responses of free and forced vibration of fractional viscoelastic nanobeam are discussed.Results and ConclusionThe numerical results show that the fractional order must be considered in the modeling of viscoelastic nanobeam, and the system’s damping enhances by the increase of fractional order. Moreover, owing to the correlation between fractional order and excitation frequency, the nonlinear time responses of fractional order to free and forced vibration are different.