In order to describe the three-stage creep behavior of compressed asphalt mastic, a visco-elastoplastic damage constitutive model is proposed in this work. The model parameters are treated as quadratic polynomial functions with respect to stress and temperature. A series of uniaxial compressive creep experiments are performed at various stress and temperature conditions in order to determine these parameter functions, and then the proposed model is validated by comparison between the predictions and experiments at the other loading conditions. It is shown that very small permanent deformation at low stress and temperature increases rapidly with elevated stress or temperature and the damage may initiate in the stationary stage but mainly develops in the accelerated stage. Compared with the visco-elastoplastic models without damage, the predictions from the proposed model is in better agreement with the experiments, and can better capture the rate-dependency in creep responses of asphalt mastic especially below its softening point of 47 ℃