In this paper, a methodology for modeling creep behavior of asphalt concretes (AC) under uniaxial static loading is presented. Three AC mixtures were produced using Coarse Matrix High Binder type C (CMHB?C) mixture type. Three different aggregate sources, namely a Hard Limestone (HL), a Soft Limestone (SL), and a Granite (G) and a PG 76?22 binder were used to prepare the mixes. The microstructure of AC was captured using X?Ray computed tomography (CT) and the Discrete Element Method (DEM) was utilized to describe the microstructure. An automated digital image processing technique called Volumetric?based Global Minima (VGM) thresholding algorithm was utilized to process the X?Ray CT images of AC cores. The viscoelastic properties of asphalt binders and their mastics were characterized by fitting the Burger models to Dynamic Shear Rheometer (DSR) mastic data. The DEM simulations results were compared favorably to experimental uniaxial creep data. It was found that the models underestimated the creep compliance of the mixtures in the primary creep stage and predicted quite accurately in the secondary creep stage.