Extensive experimental studies have demonstrated the time-dependent mechanical behaviors of frozen soil. Nonetheless, limited studies are focusing on the constitutive modeling of the time-dependent stress-strain behaviors of frozen clay soils at different subzero temperatures. The objective of this study is to numerically investigate the time-dependent behavior of frozen clay soils at a temperature range of 0 degrees C to - 15 degrees C. The Drucker-Prager model is adopted along with the Singh-Mitchell creep model to simulate time-dependent uniaxial compression and stress relaxation behaviors of frozen sandy clay soil. The numerical modeling is implemented through the finite element method based on the platform of Abaqus. The constitutive modeling is calibrated by a series of experimental results on laboratory-prepared frozen sandy clay soils, where the strain hardening, the post-peak softening, and stress relaxation behaviors are captured. Our results show that both the rate-dependent model and creep model should be adopted to characterize a comprehensive time-dependent behavior of frozen soils. The rate-dependent stress-strain behaviors heavily rely on the rate- and temperature-dependent hardening functions, where the creep strain provides a very limited contribution. Nevertheless, the creep strain should also be adopted when a long-term analysis or stress relaxation behavior is involved.