A finite element model based on the incremental endochronic theory for flexible pavement materials was developed in this study. Three grid systems with eight-node cubic isoparametric elements, and different loading steps were used to perform the calculations for a specimen of circular cylinder. The uniaxial stress experimental results on an asphalt mixture at $60^{circ}C$ in SHRP conducted by University of California at Berkeley were used to check the ability of the derived numerical model. Then, the numerical results showed isotropic response and deviatoric response on the specimen in a three dimensional manner, which provided a better understanding for a deformed flexible material under the specified loading conditions.