A variety of hot mix asphalt mixtures are used in highway and runway pavement construction. Each mixture caters to specific needs. Mixtures differ from one another in the type and percentage of aggregates and asphalt used, and their response can be markedly different, and thus there is a need to develop constitutive models that can differentiate between the different kinds of mixtures. In this paper, we outline a general procedure for the constitutive modeling of bituminous mixtures. We illustrate the efficacy of this approach by means of an application to sand asphalt. The governing equations for this special problem reduce to a stiff nonlinear ordinary differential equation and this is solved numerically using Gear's method. We compare the results of the predictions of the model that we have developed with the compressive creep experiments carried out by Wood and Goetz on a typical sand asphalt mixture and find them to be in good agreement. © ASCE.