Time-varying loading is a frequently encountered loading type in geotechnical engineering. As the deformation of viscoelastic soil is related to its loading history, studying the viscoelastic problems under time-varying loads has important practical engineering significance. In this paper, the stress and displacement of a layered soil with fractional-order viscoelastic model under time-varying loads were solved using the complex variable method and the corresponding principle. Under the assumption of quasi-static and linear elasticity, this paper derives the quasi-static elastic solutions of a layered soil under time-varying strip loads. By introducing the fractional-order viscoelastic model and using the corresponding principle, the Laplace-domain analytical solutions are obtained. Finally, numerical methods are utilized to perform the Laplace inverse transform and obtain the solutions in the physical domain. The correctness of this paper is validated by comparing the numerical results with the ANSYS software, as well as by comparison with previous literature. Based on the experimental data, the material parameters of frozen soil at two temperatures are fitted. The settlement patterns of soil under two engineering loads were analyzed through case studies. By controlling variables, the influence of parameters of the fractional viscoelastic model on elastic aftereffect was investigated.