Next Article in Journal
Metal Nanoparticles as Green Catalysts
Next Article in Special Issue
Rheological Properties, Compatibility, and Storage Stability of SBS Latex-Modified Asphalt
Previous Article in Journal
Ultrasound-Mediated Atom Transfer Radical Polymerization (ATRP)
Previous Article in Special Issue
Study of the Self-Healing Performance of Semi-Flexible Pavement Materials Grouted with Engineered Cementitious Composites Mortar based on a Non-Standard Test
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Aggregate Mesostructure on Permanent Deformation of Asphalt Mixture Using Three-Dimensional Discrete Element Modeling

1
School of Transportation, Southeast University, Nanjing 210096, China
2
School of Architecture Engineering, Nanjing Institute of Technology, Nanjing 211167, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(21), 3601; https://doi.org/10.3390/ma12213601
Submission received: 26 September 2019 / Revised: 24 October 2019 / Accepted: 31 October 2019 / Published: 2 November 2019
(This article belongs to the Special Issue Sustainable Designed Pavement Materials)

Abstract

:
This paper investigated the effects of aggregate mesostructures on permanent deformation behavior of an asphalt mixture using the three-dimensional (3D) discrete element method (DEM). A 3D discrete element (DE) model of an asphalt mixture composed of coarse aggregates, asphalt mastic, and air voids was developed. Mesomechanical models representing the interactions among the components of asphalt mixture were assigned. Based on the mesomechanical modeling, the uniaxial static load creep tests were simulated using the prepared models, and effects of aggregate angularity, orientation, surface texture, and distribution on the permanent deformation behavior of the asphalt mixtures were analyzed. It was proven that good aggregate angularity had a positive effect on the permanent deformation performance of the asphalt mixtures, especially when approximate cubic aggregates were used. Aggregate packing was more stable when the aggregate orientations tended to be horizontal, which improved the permanent deformation performance of the asphalt mixture. The influence of orientations of 4.75 mm size aggregates on the permanent deformation behavior of the asphalt mixture was significant. Use of aggregates with good surface texture benefitted the permanent deformation performance of the asphalt mixture. Additionally, the non-uniform distribution of aggregates had a negative impact on the permanent deformation performance of the asphalt mixtures, especially when aggregates were distributed non-uniformly in the vertical direction.

1. Introduction

With the increase of traffic and heavy vehicles, rutting has become one of the main distresses of asphalt pavement, which seriously affects driving safety and comfort. The permanent deformation behavior of asphalt mixtures has become an important issue for road researchers. Within asphalt mixtures, aggregates account for more than 90% by weight and 80% by volume. Aggregates, as the main component of asphalt mixtures, have a significant impact on the permanent deformation behavior of asphalt mixtures. Previous studies on the effects of aggregates on the permanent deformation behavior of asphalt mixtures have been mainly related to aggregate gradations. It is difficult to investigate the effects of the mesostructural characteristics of aggregates on permanent deformation behavior of asphalt mixtures, due to the limitations of available laboratory tests.
In previous studies, X-ray computed tomography (CT) techniques and image processing technology have been used for aggregate mesostructural analysis. Masad proposed computer-automated image analysis procedures to quantify the internal structure of asphalt concrete in terms of aggregate orientation, aggregate contacts, and air void distribution, and developed a finite-element model of asphalt concrete mesostructure to study the influence of localized strain distribution on the mechanical response of asphalt concrete [1,2,3]. Kose investigated strain distribution within a binder using digitized images analyzed using finite-element procedures [4]. Wang developed a method to quantify the local volume fractions of voids and their spatial gradients using X-ray tomography imaging and image analysis [5,6]. You predicted the dynamic modulus of asphalt mixtures using both the two-dimensional and three-dimensional discrete element method, generated using X-ray computed tomography [7]. However, pre-compacted specimens should have been prepared in the laboratory for control purposes in the studies mentioned above. Variability of the components within specimens have inevitably occurred, and the distribution, morphological characteristics of aggregates, and air voids are not the same even within two mixture specimens of the same kind. Therefore, there are many factors that might disrupt investigation of the effect of a single mesostructure on the behavior of asphalt mixtures. Moreover, the mesostructures of the components cannot be desirably controlled in the laboratory, and the results to date have been consequently unconvincing. Therefore, it is necessary to control the internal structures of asphalt mixtures in order to accurately analyze the effects of a single mesostructure.
At present, commonly used numerical simulation methods at the meso-scale include the finite element method (FEM) [8], the discrete element method (DEM) [9], and the combined finite–discrete element method (FDEM) [10]. In recent years, DEM, as a numerical simulation method used to solve the problem of discontinuity media, has been used to analyze the mechanical properties of asphalt mixtures. Chang developed a model called ASBAL (TRUBAL for ASphalt) based on the discrete element method by modifying the TRUBAL program to simulate hot-mix asphalt mixtures [9]. Buttlar presented a microfabric discrete-element modeling approach (MDEM) for modeling asphalt concrete mesostructures using image analysis techniques [11]. You predicted the asphalt mixture complex modulus in extension/compression across a range of test temperatures and load frequencies using the MDEM approach, and simulated and analyzed the creep responses of an asphalt mixture with a 3D-mesostructure-based DE viscoelastic model [12,13]. Abbas presented a methodology for analyzing the viscoelastic response of asphalt mixtures using the DEM [14]. Liu developed a viscoelastic model of asphalt mixtures using the discrete element method, where the viscoelastic behaviors of asphalt mastics are represented by a Burger’s model [15]. Chen built a mesomechanical model to investigate the stiffness anisotropy of asphalt concrete using the DEM [16,17]. Ma conducted simulated wheel-tracking tests on asphalt mixtures, predicted the rutting deformation of asphalt mixtures, and analyzed the mesomechanical response of aggregate skeletons within the asphalt mixtures during the tests [18,19]. Ma also investigated the effects of different parameters related to air voids on the creep behavior of asphalt mixtures based on a 3D DEM [20]. Ding proposed a new modeling method to reconstruct hollow shapes of aggregate particles using the DEM to accurately characterize the mesostructures of aggregates and efficiently predict the mechanical properties of aggregate skeletons [21]. It has been proven that the DEM can accurately capture the internal structure of asphalt mixtures, and desirably control the mesostructural characteristics of aggregates, asphalt binder, and air voids within the mixtures [22,23,24,25]. However, few researchers have focused on the permanent deformation behavior of asphalt mixtures, especially the effects of aggregate mesostructure on the permanent deformation behavior of asphalt mixtures using the DEM.
The objective of this study was to investigate the effects of aggregate mesostructures on the permanent deformation behavior of asphalt mixtures using the DEM. A 3D DE model of asphalt mixture including aggregates, asphalt mastics, and air voids was developed using the DEM software Particle Flow Code in three dimensions (PFC3D). Three-dimensional discrete element simulations of uniaxial static creep tests of asphalt mixtures were carried out. By accurately controlling the mesostructural characteristics of the aggregates, such as angularity, orientation, surface texture, and distribution, the effects of aggregate mesostructures on the permanent deformation behavior of the asphalt mixtures were carefully analyzed.

2. Materials and Methods

2.1. Materials

According to the Marshall mix design method based on the Chinese specifications [26], a dense-graded asphalt mixture, AC-20, was prepared in the laboratory. Aggregate gradation with a nominal maximum aggregate size of 19 mm is shown in Figure 1. Based on the Marshall mix design method, the volumetric properties of the AC-20 asphalt mixture were determined as shown in Table 1. Aggregates smaller than 2.36 mm, mineral filler, and asphalt binder were mixed as asphalt mastic with an asphalt content of 11.5%. The air void content of the asphalt mastic was supposed to be zero to maximize its flowability.

2.2. Laboratory Tests

The uniaxial static creep test is one of the simplest and most practical test methods, and it is an effective way to investigate the permanent deformation behaviors of asphalt mixtures and mastics at high temperatures. In this test, an instantaneous load is applied to a cylinder specimen in the axial direction and the load is kept constant, and the creep curve of the asphalt mixtures or mastics can be obtained. Uniaxial static creep tests were conducted in this paper to evaluate the permanent deformation behaviors of asphalt mixtures and mastics at the temperature of 60 °C using a universal test machine (UTM, IPC Global, Melbourne, Australia) [27]. Cylindrical specimens of asphalt mixtures with diameter of 100 mm and height of 150 mm were prepared using a gyratory compactor for the uniaxial static creep test. Cylindrical asphalt mastic specimens with diameter and height of 100 mm were prepared by vibration. The applied axial stresses for the asphalt mixture and asphalt mastic were set at 0.7 MPa and 0.07 MPa, respectively, in this study. The asphalt mastic specimens and the testing process are shown in Figure 2. Lab test results were used for parameters of DE modeling and to be compared with simulation test results, and are presented in the following sections. The experimental program for the asphalt mixture and mastics is summarized in Table 2.

3. Discrete Element Modeling of Asphalt Mixtures

3.1. Discrete Element Modeling

An asphalt mixture, as a multiphase composite, is composed of aggregates, asphalt binder, and air voids. If the particle size of the fine aggregates and mineral filler is small, it will lead to a large increase of discrete elements in the DE model, which is bound to significantly reduce the computational efficiency if the fine aggregate and mineral powder are considered fully in DE models of asphalt mixtures. Therefore, the asphalt mixture was simplified into coarse aggregates (bigger than 2.36 mm), asphalt mastics (asphalt binder, fine aggregates smaller than 2.36 mm and mineral filler), and air voids to improve the calculation efficiency. The asphalt mastic was considered as a homogeneous material within the DE model.
During modeling, the spatial range of the mixture model was first constructed by “wall”. The quantity of coarse aggregates in each grade was calculated according to aggregate gradation, asphalt content and air void content. The coarse aggregate balls were then delivered into the spatial range of the mixture model constructed by “wall”. Force was generated between overlapped balls, as there were overlaps between the coarse aggregate balls during delivery. This force was eliminated by the “cycle” command. The coarse aggregate balls within the spatial range of the mixture model are shown in Figure 3a.
A regular array of discrete elements was then filled into the mixture space as the base of coarse aggregates and asphalt mastic, as shown in Figure 3b. The radius of the discrete element was set to 1 mm, considering the calculation efficiency.
It has been proven that the geometry of coarse aggregates has a significant effect on the permanent deformation behavior of asphalt mixtures [28,29]. Therefore, the coarse aggregate geometry will directly determine the accuracy of the simulation results. In this paper, coarse aggregate was simplified as an irregular polyhedron. Random planes were used to cut a cube or a sphere to generate irregular polyhedral aggregates using a user-defined program. By traversing the regularly packed discrete elements, the positional relationship between the regular packing discrete elements and the irregular polyhedral aggregates could be evaluated. Discrete elements belonging to irregular polyhedral aggregate were regarded as aggregate elements, and were set as a clump, as seen in Figure 4. The original coarse aggregate balls were then deleted. Discrete elements outside the irregular polyhedron aggregate were considered as asphalt mastic. The initially developed 3D DE model of asphalt mixture is shown in Figure 5a.
Considering the complexity of air void distribution within an asphalt mixture specimen, the air voids were assumed to be randomly distributed [30]. Asphalt mastic elements in the model were traversed and randomly deleted and regarded as air voids. The distribution of air voids is shown in Figure 5b.

3.2. Mesomechanical Models and Parameters

There are four types of contact within asphalt mixtures, including contacts between aggregate elements, contacts between adjacent aggregates, contacts between asphalt mastic elements, and contact between asphalt mastic and aggregates, as seen in Figure 6. In PFC3D, mesomechanical models are used to describe the contact behaviors between different components within asphalt mixtures. In this study, the mesomechanical models used in the model of asphalt mixtures included the stiffness model, slipping model, bonding model, and Burger’s model.
Due to the high stiffness of an aggregate, it can be approximately regarded as an elastic material. In this study, the stiffness model and slipping model were used to characterize the mesomechanical behavior between adjacent aggregates. As coarse aggregates within the model of asphalt mixtures were set as clumps, it was unnecessary to assign the mesomechanical model within aggregates. The mesomechanical parameters of the stiffness model could be obtained from the macro properties of the aggregates, as shown in Equations (1) and (2) [15,31,32]. The macro parameters for the aggregates are shown in Table 3 [9,12,33].
E = k n 4 R ,   k s = k n 2 ( 1 + υ )
μ c = μ a
where E is the apparent Young’s modulus of the aggregates, kn and ks are the stiffness in the normal and shear direction, respectively, R is the discrete element radius, υ′ is the aggregate Poisson’s ratio, μc is the friction coefficient between aggregates, and μa is the friction coefficient of aggregates.
Asphalt mixtures show a macro viscoelastic behavior due to the viscoelastic characteristic of the asphalt mastic. Therefore, the mesomechanical model of the asphalt mastic directly affects the macro properties of asphalt mixtures. The meso Burger’s model in PFC3D was well able to describe the mechanical properties of viscoelastic materials and was used to characterize the viscoelastic properties of the asphalt mastic in this study, as shown in Figure 7. It has been proven that there is a conversion between the parameters of the meso Burger’s model and the macro Burger’s model [15,31], as shown in Equations (3) and (4). Parameters of meso Burger’s model could then be obtained from the macro Burger’s model parameters. The macro Burger’s model parameters of asphalt mastic were obtained by uniaxial static creep test, as shown in Table 3.
K m n = E 1 L ,   C m n = η 1 L ,   K k n = E 2 L ,   C k n = η 2 L
K m s = E 1 L 2 ( 1 + υ ) ,   C m s = η 1 L 2 ( 1 + υ ) ,   K k s = E 2 L 2 ( 1 + υ ) ,   C k s = η 2 L 2 ( 1 + υ )
where E1, η1, E2, and η2, are parameters of the macro Burger’s model; Kmn, Cmn, Kkn, and Ckn are parameters of the meso Burger’s model in the normal direction; Kms, Cms, Kks, and Cks are parameters of the meso Burger’s model in the shear direction; L is the distance between adjacent discrete elements; and υ is Poisson’s ratio of the asphalt mastic.
The contact behavior between aggregates and the asphalt mastic can be characterized by the equivalent meso Burger’s model, as shown in Figure 8. The equivalent meso model parameters can be obtained from the macro parameters of the aggregates and asphalt mastic, as expressed in Equations (5) and (6) [15,31]:
K m n = 2 E E 1 E + E 1 L ,   C m n = 2 η 1 L ,   K k n = 2 E 2 L ,   C k n = 2 η 2 L
K m s = 2 E E 1 E ( 1 + υ ) + 2 E 1 ( 1 + υ ) L ,   C m s = η 1 L ( 1 + υ ) ,   K k s = E 2 L ( 1 + υ ) ,   C k s = η 2 L ( 1 + υ )
where Kmm, Cmm, Kkn, and Ckn are the parameters of the equivalent meso Burger’s model between the aggregates and mastic in the normal direction and Kms, Cms, Kks, and Cks are the parameters of the equivalent meso Burger’s model between the aggregates and the mastic in the shear direction.

3.3. Simulation of Uniaxial Creep Test

Based on the above-mentioned DE modeling procedures for asphalt mixtures, the simulated uniaxial creep test was carried out using PFC3D. The results were compared with laboratory test results under the same conditions. Figure 9 shows the axial strain curve obtained from the simulation and lab test, and it can be seen that the simulation curve was very close to the curve of the laboratory test. This indicates that the simulation of uniaxial creep test conducted using DEM could precisely estimate the permanent deformation behavior of asphalt mixtures.

4. Results and Discussion

4.1. Effect of Aggregate Angularity

The aggregate angularity could represent the morphological characteristics of aggregates, and significantly affect the skeleton stability, interlocking force, and shear resistance of aggregates. It has been proven that the aggregate angularity is closely related to the strength of asphalt mixtures, especially the permanent deformation resistance and shear resistance [28,29].
Some previous studies about aggregate angularity have been carried out, and some indexes to describe the aggregate angularity have been proposed. These angularity indexes are mainly divided into two categories, one denoting the roundness of aggregate corners and the other one denoting the roundness of the overall outline of the aggregate [34]. However, most of these angularity indexes are in two dimensions. Moreover, the measurement and calculation of the first type of angularity index is very complex, while the second type could be affected by the overall outline of aggregates. Therefore, the angularity indexes mentioned above cannot fully represent the angularity characteristics of aggregates. In this study, the ratio of surface area of irregular aggregate to equivalent ellipsoid was proposed to characterize the aggregate angularity. The equivalent ellipsoid of an aggregate has the same volume as the aggregate. It was considered that the equivalent ellipsoid of aggregates could accurately reflect the overall outline, and this angularity index could significantly reduce the influence of aggregate outline on angularity index, as expressed in Equation (7). The larger AI is, the better the aggregate angularity is.
A I = ( S S e l l i p s e ) 3
where AI is the angularity index of aggregates, S is surface area of aggregates, and Sellipse is surface area of equivalent ellipsoid.
As mentioned above, irregular polyhedral aggregates were generated by cutting a cube or a sphere with random planes. Thus, the equivalent ellipsoid of an irregular polyhedral aggregate can be approximated to an equivalent sphere, and then Equation (8) can be further expressed as follows:
A I = ( S S s p h e r e ) 3 = ( S 4 π ( 3 V / 4 π 3 ) 2 ) 3 = S 3 36 π V 2
where Ssphere is surface area of equivalent sphere and V is the volume of aggregates.
Coarse aggregates with different angularities were generated by cutting a cube or a sphere with random planes. Some coarse aggregates with different angularities are shown in Figure 10, in which their angularity become better with the increasing number. DE models of asphalt mixtures with different aggregate angularities were then obtained.
In order to calculate the surface area of aggregates, the number of discrete elements on aggregate surface and the number of discrete elements that made up the aggregates were counted. The angularity index of the aggregates was then calculated by Equation (9). The method used to count the number of discrete elements on the aggregate surface and the number of discrete elements that make up aggregates is shown in Figure 11.
A I = ( 4 R 2 n S ) 3 36 π ( 8 R 3 n V ) 2 = n S 3 36 π n V 2
where ns is the number of discrete elements on the aggregate surface, nv is the number of discrete elements that make up the aggregate, and R is radius of the discrete element.
The angularity of a single aggregate can be obtained by Equation (9). However, there are many aggregates within asphalt mixtures, and the angularities of each aggregate are different. It is necessary to evaluate the overall angularity of the aggregates within asphalt mixtures. In this study, the volumetrically weighted average of the angularities of each aggregate was used to describe the overall angularity of the aggregates within the asphalt mixtures, as expressed in Equation (10).
A I ¯ = i = 1 n V i × A I i i = 1 n V i
where A I ¯ is the overall angularity of aggregates, AIi is the angularity of aggregate i, Vi is the volume of aggregate i, and n is the number of aggregates.
Simulations of uniaxial static creep tests were then conducted on the DE models of asphalt mixtures with different overall aggregate angularities. The results are shown in Figure 12. The permanent deformation performance of asphalt mixtures improved with the increasing aggregate angularity, but the improvement weakened gradually. The axial deformation of the asphalt mixture was the smallest when the aggregates were closest to a cube. This is consistent with the requirement of selecting aggregates with shapes close to a cube. Therefore, aggregates with good angularities should be selected to improve the permanent deformation resistance of asphalt pavement.

4.2. Effect of Aggregate Orientation

Previous studies about aggregate orientation have generally defined the long axis of the equivalent ellipse of aggregate in two dimensions as the long axis of aggregate, and used the angle between the long axis and the horizontal direction to represent the aggregate orientation [18,19,20,21]. In this study, the aggregate orientation was defined in three dimensions. The long axis of the equivalent ellipsoid of the aggregate was taken as the long axis, and the angles α, β, and γ between the long axis and X axis, Y axis, and Z axis were used to represent the aggregate orientation in three dimensions, as shown in Figure 13.
In order to investigate the effect of aggregate orientation on the permanent deformation performance of asphalt mixtures, a subroutine was written in PFC3D with the built-in “Fish” language. The long axis of the aggregate was set to be parallel to the YOZ surface, and the aggregate orientation was described by the angle between the long axis of the aggregate and the XOY surface. Some coarse aggregates with different orientations are shown in Figure 14.
Graded coarse aggregates with different orientations were generated and randomly put into the cylinder space of the asphalt mixture model. The unbalanced forces between the aggregates were eliminated by the “cycle” command, so that the coarse aggregates within the asphalt mixtures could reach equilibrium. The angular velocities of the coarse aggregates were fixed to keep the aggregate orientations constant during modeling.
The effects of the orientations of aggregates of different sizes on the mechanical properties of the asphalt mixtures are different. To address the size effect, 30 cylinder DE models of asphalt mixtures with different orientations of aggregates with different sizes were developed. The aggregate orientations were varied from 0°, 30°, 45°, 60°, to 90°. The rest of the aggregates remained spherical. Figure 15 shows the coarse aggregates with different orientations.
Simulations of uniaxial static creep tests were then conducted on the prepared DE models. The results are shown in Figure 16.
It can be seen from Figure 16 that the axial strains of the mixtures increased gradually with the increasing aggregate orientations. The orientations of 19–26.5mm and 16–19 mm aggregates had insignificant effects on the permanent deformation behavior of asphalt mixtures. The orientations of 13.2–16 mm, 9.5–13.2 mm, and 2.36–4.75 mm aggregates had more significant effects than those of 19–26.5 mm and 16–19 mm. The orientation of 4.75–9.5 mm aggregates had the greatest effect on the permanent deformation behavior of the asphalt mixtures. The reason is that the aggregates tended to rotate horizontally under vertical loads when the aggregate orientation was large. Thus, the aggregate skeleton was unstable and led to large axial deformation of the asphalt mixtures. When the aggregate orientation was closer to horizontal, the aggregate skeleton was more stable, and the axial deformation of the mixture specimens was smaller. It has been shown that the percentage pass of 4.75 mm plays an important role in controlling the aggregate skeleton, and the influence is next to that of air void content on the permanent deformation behavior of asphalt mixtures. This may be the reason why the orientation of 4.75–9.5 mm aggregates had such a great influence. Therefore, different kinds of rollers should be used during rolling compaction of asphalt pavement, which could force the aggregate orientation to be more horizontal. Thus, the asphalt pavement would reach a more stable state and its permanent deformation performance could be significantly improved.

4.3. Effect of Aggregate Surface Texture

The surface texture of aggregates reflects the roughness of aggregate surface. The rougher the aggregate surface, the better the surface texture. Internal friction between aggregates and adhesion between aggregate and asphalt binder are closely related to the surface texture of aggregates. It has been proven that the surface texture of aggregates has a significant impact on the performance of asphalt mixtures, especially the permanent deformation performance. The randomly created irregular polyhedral aggregates in the DE model could only reflect geometry and angularity, without characterizing the surface texture of aggregates. In this study, the friction coefficient between aggregates was used to indirectly represent the surface texture of aggregates. The higher the friction coefficient, the rougher the surface texture. A series of DE models of asphalt mixture specimens with different aggregate friction coefficients assigned were developed. The internal structures of the developed DE models were the same, and only the aggregate friction coefficients of polyhedron aggregates in the DE models were different. Simulations of uniaxial static creep tests were conducted on the prepared DE models. The results are shown in Figure 17.
It can be seen from Figure 17 that the permanent deformation performance of the asphalt mixtures was improved with the increasing friction coefficient, but the increment was gradually lowered. The permanent deformation resistance of the asphalt mixtures increased with rougher aggregate surface texture. Therefore, aggregates with good surface texture should be used to improve the permanent deformation performance of asphalt pavement.

4.4. Effect of Aggregate Distribution

Aggregates are ideally uniformly distributed within asphalt mixture for optimum performance. However, insufficient mixing force or mixing time may cause non-uniform distribution of aggregates. Aggregate segregation may also occur during paving. Early distresses of asphalt pavement such as cracking, pothole, blooding, and rutting are closely related to the material inhomogeneity caused by variable construction quality. The non-uniform distribution of aggregates greatly influences the internal structure of asphalt mixtures and consequently lowers the permanent deformation performance of asphalt mixtures.
To evaluate the uniformity of aggregate distribution, a cubic space containing asphalt mixtures was constructed by “wall” and the specimen was divided into three parts evenly in both the horizontal and vertical directions, as shown in Figure 18. Coarse aggregates with different sizes were added to each of the three parts, as shown in Figure 19 to create non-uniform distribution. The total volume of aggregates within the three parts remained constant. The distribution percentages of the aggregate volume within each part are shown in Table 4.
A quantitative indicator of aggregate distribution was proposed to measure the uniformity of distribution of aggregate. The volume of aggregates within each part was recorded and the coefficient of variation (var) of aggregate distribution was calculated as follows:
V j = i = 1 m V i
V ¯ = j = 1 n V j / n
v a r j = ( V j V ¯ ) 2 / ( n 1 ) / V ¯
v a r = j = 1 n v a r j
where Vi is the volume of aggregate i within part j, Vj is the total volume of aggregates within part j, V ¯ is the average aggregate volume of all divided parts, varj is the coefficient of variation of the total volume of aggregates within part j, and var is the coefficient of variation of the total volume of aggregates within specimen.
Figure 20 and Table 5 below present the results, in which the axial strain was calculated by measuring the axial deformation and dividing it by the initial axial length.
As shown in Figure 20 and Table 5, axial deformation of the specimens increased as the coefficient of variation of aggregate distribution increased using both partition methods. Aggregate distribution had a significant impact on the deformation performance of the asphalt mixtures. Non-uniformity in the vertical direction had much more impact on the deformation of the specimen than that in the horizontal direction, as shown in Figure 20.

5. Conclusions

In this study, three-dimensional discrete element (DE) models of asphalt mixtures composed of coarse aggregates, asphalt mastic, and air voids was developed using PFC3D, and uniaxial static load creep tests were conducted on the prepared models to investigate the effects of aggregate mesostructure on the permanent deformation behavior of the asphalt mixtures. The effects of aggregate angularity, orientation, surface texture, and distribution on the permanent deformation behavior of asphalt mixtures were carefully analyzed. The following conclusions were drawn:
(1) The constructed DE models of asphalt mixtures could accurately capture the mesostructures of aggregates, such as angularity, orientation, surface texture, and distribution. The simulation of the uniaxial creep test conducted using the DEM could precisely estimate the permanent deformation behavior of asphalt mixtures.
(2) The permanent deformation performance of the asphalt mixtures was improved with increasing aggregate angularity. The axial deformation of the asphalt mixtures was smallest when the aggregates were close to cubical. Good aggregate angularity had a positive effect on the permanent deformation performances of asphalt mixtures, especially when the aggregates were nearly cubical.
(3) Aggregate packing was more stable when the orientations tended to be horizontal, which improved the permanent deformation performances of asphalt mixtures. The percentage pass of 4.75 mm played an important role in controlling the aggregate skeleton and had a great impact on the permanent deformation behaviors of the asphalt mixtures.
(4) The permanent deformation performance of the asphalt mixtures was improved with increasing friction coefficient, indicating that aggregates with good surface texture benefitted the permanent deformation performance of the asphalt mixtures.
(5) Axial deformation of the asphalt mixture specimen increased with increasing coefficient of variation of aggregate distribution. Non-uniform distribution of the aggregates had a negative impact on the permanent deformation performance of the asphalt mixtures, especially when the aggregates were distributed non-uniformly in the vertical direction.
In the future research steps, a viscoplastic mesomechanical model will be developed to further describe the viscoplastic behavior of the permanent deformation behavior of asphalt mixtures.

Author Contributions

Conceptualization, D.Z. and J.Z.; Data curation, D.Z., L.G.; Formal analysis, D.Z., L.G. and J.Z.; Funding acquisition, D.Z.; Methodology, D.Z. and J.Z.; Writing—original draft, D.Z.; Writing—review & editing, J.Z.

Funding

The study is financially supported by Postdoctoral Research funding Project of Jiangsu (2019K138), the Natural Science Foundation of Jiangsu (BK20140109) and the Research Foundation of Nanjing Institute of Technology (YKJ201512).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Masad, E.; Muhunthan, B.; Shashidhar, N.; Harman, T. Quantifying laboratory compaction effects on the internal structure of asphalt concrete. Transp. Res. Rec. 1999, 1681, 179–185. [Google Scholar] [CrossRef]
  2. Masad, E.; Muhunthan, B.; Shashidhar, N.; Harman, T. Internal structure characterization of asphalt concrete using image analysis. J. Comput. Civ. Eng. 1999, 13, 88–95. [Google Scholar] [CrossRef]
  3. Masad, E.; Somadevan, N. Microstructural finite-element analysis of influence of localized strain distribution on asphalt mix properties. J. Eng. Mech. 2002, 128, 1105–1114. [Google Scholar] [CrossRef]
  4. Kose, S.; Guler, M.; Bahia, H.; Masad, E. Distribution of strains within hot-mix asphalt binders: Applying imaging and finite-element techniques. Transp. Res. Rec. 2000, 1728, 21–27. [Google Scholar] [CrossRef]
  5. Wang, L.B.; Frost, J.D.; Shashidhar, N. Microstructure study of westrack mixes from x-ray tomography images. Transp. Res. Rec. J. Transp. Res. Board 2001, 1767, 85–94. [Google Scholar] [CrossRef]
  6. Wang, L.; Wang, Y.; Mohammad, L.; Harman, T. Voids distribution and performance of asphalt concrete. Int. J. Pavements 2002, 1, 22–33. [Google Scholar]
  7. You, Z.; Adhikari, S.; Kutay, M.E. Dynamic modulus simulation of the asphalt concrete using the X-ray computed tomography images. Mater. Struct. 2009, 42, 617–630. [Google Scholar] [CrossRef]
  8. Ding, X.; Ma, T.; Gu, L.; Zhang, Y. Investigation of surface micro-crack growth behavior of asphalt mortar based on the designed innovative mesoscopic test. Mater. Des. 2020, 185, 108238. [Google Scholar] [CrossRef]
  9. Chang, K.-N.G.; Meegoda, J.N. Micromechanical Ssimulation of hot mix asphalt. J. Eng. Mech. 1997, 123, 495–503. [Google Scholar] [CrossRef]
  10. Rougier, E.; Knight, E.E.; Broome, S.T.; Sussman, A.J.; Munjiza, A. Validation of a three-dimensional finite-discrete element method using experimentalresults of the Split-Hopkinson pressure bar test. Int. J. Rock Mech. Min. Sci. 2014, 70, 101–108. [Google Scholar] [CrossRef]
  11. Buttlar, W.G.; You, Z. Discrete element modeling of asphalt concrete: Microfabric approach. Transp. Res. Rec. J. Transp. Res. Board 2001, 1757, 111–118. [Google Scholar] [CrossRef]
  12. You, Z.; Buttlar, W.G. Discrete element modeling to predict the modulus of asphalt concrete mixtures. J. Mater. Civ. Eng. 2004, 16, 140–146. [Google Scholar] [CrossRef]
  13. You, Z.; Liu, Y.; Dai, Q. Three-dimensional microstructural-based discrete element viscoelastic modeling of creep compliance tests for asphalt mixtures. J. Mater. Civ. Eng. 2011, 23, 79–87. [Google Scholar] [CrossRef]
  14. Abbas, A.; Masad, E.; Papagiannakis, T.; Harman, T. Micromechanical modeling of the viscoelastic behavior of asphalt mixtures using the discrete-element method. Int. J. Géoméch. 2007, 7, 131–139. [Google Scholar] [CrossRef]
  15. Liu, Y.; Dai, Q.; You, Z. Viscoelastic model for discrete element simulation of asphalt mixtures. J. Eng. Mech. 2009, 135, 324–333. [Google Scholar] [CrossRef]
  16. Chen, J.; Zhang, Q.; Wang, H.; Wang, L.; Huang, X. Numerical investigation into the effect of air voids on the anisotropy of asphalt mixtures. J. Wuhan Univ. Technol.-Mater. Sci. Ed. 2017, 32, 473–481. [Google Scholar] [CrossRef]
  17. Chen, J.; Pan, T.; Huang, X. Numerical investigation into the stiffness anisotropy of asphalt concrete from a microstructural perspective. Constr. Build. Mater. 2011, 25, 3059–3065. [Google Scholar] [CrossRef]
  18. Ma, T.; Zhang, D.; Zhang, Y.; Hong, J. Micromechanical response of aggregate skeleton within asphalt mixture based on virtual simulation of wheel tracking test. Constr. Build. Mater. 2016, 111, 153–163. [Google Scholar] [CrossRef]
  19. Ma, T.; Zhang, D.; Zhang, Y.; Wang, S.; Huang, X. Simulation of wheel tracking test for asphalt mixture using discrete element modelling. Road Mater. Pavement Des. 2016, 19, 367–384. [Google Scholar] [CrossRef]
  20. Ma, T.; Zhang, D.; Zhang, Y.; Zhao, Y.; Huang, X. Effect of air voids on the high-temperature creep behavior of asphalt mixture based on three-dimensional discrete element modeling. Mater. Des. 2016, 89, 304–313. [Google Scholar] [CrossRef]
  21. Ding, X.; Ma, T.; Huang, X. Discrete-element contour-filling modeling method for micromechanical and macromechanical analysis of aggregate skeleton of asphalt mixture. J. Transp. Eng. Part B: Pavements 2019, 145, 04018056. [Google Scholar] [CrossRef]
  22. Zhang, Y.; Ma, T.; Ling, M.; Huang, X. Mechanistic sieve-size classification of aggregate gradation by characterizing load-carrying capacity of inner structures. J. Eng. Mech. 2019, 145, 04019069. [Google Scholar] [CrossRef]
  23. Zhang, Y.; Ma, T.; Ling, M.; Zhang, D.; Huang, X. Predicting dynamic shear modulus of asphalt mastics using discretized-element simulation and reinforcement mechanisms. J. Mater. Civ. Eng. 2019, 31, 04019163. [Google Scholar] [CrossRef]
  24. Zhang, Y.; Ma, T.; Huang, X.; Zhao, Y.; Hu, P. Algorithms for generating air-void structures of idealized asphalt mixture based on three-dimensional discrete-element method. J. Transp. Eng. Part B: Pavements 2018, 144, 04018023. [Google Scholar] [CrossRef]
  25. Zhang, Y.; Ma, T.; Ding, X.; Chen, T.; Huang, X.; Xu, G. Impacts of air-void structures on the rutting tests of asphalt concrete based on discretized emulation. Constr. Build. Mater. 2018, 166, 334–344. [Google Scholar] [CrossRef]
  26. JTG F40-Technical Specification for Construction of Highway Asphalt Pavement; Ministry of Transport of the People’s Republic of China: Beijing, China, 2011.
  27. JTG E20-Standard Test Methods of Bitumen and Bituminous Mixtures for Highway Engineering; Ministry of Transport of the People’s Republic of China: Beijing, China, 2011.
  28. Pan, T.; Tutumluer, E.; Carpenter, S.H. Effect of coarse aggregate morphology on permanent deformation behavior of hot mix asphalt. J. Transp. Eng. 2015, 132, 580–589. [Google Scholar] [CrossRef]
  29. Stakston, A.D.; Bahia, H.U.; Bushek, J.J. Effect of fine aggregate angularity on compaction and shearing resistance of asphalt mixtures. Transp. Res. Rec. J. Transp. Res. Board 2002, 1789, 14–24. [Google Scholar] [CrossRef]
  30. Masad, E.; Jandhyala, V.K.; Dasgupta, N.; Somadevan, N.; Shashidhar, N. Characterization of air void distribution in asphalt mixes using x-ray computed tomography. J. Mater. Civ. Eng. 2002, 14, 122–129. [Google Scholar] [CrossRef]
  31. Itasca Consulting Group. PFC3D Version 4.0; Itasca Consulting Group: Minneapolis, MN, USA, 2008. [Google Scholar]
  32. Thornton, C. The conditions of failure of a face-centered cubic array of uniform rigid spheres. Géotechnique 2015, 29, 441–459. [Google Scholar] [CrossRef]
  33. Dai, Q.; You, Z. Prediction of creep stiffness of asphalt mixture with micromechanical finite-element and discrete-element models. J. Eng. Mech. 2007, 133, 163–173. [Google Scholar] [CrossRef]
  34. Erdoğan; Turhan, S.; Fowler, D.; Garboczi, E.J. Determination of Aggregate Shape Properties Using X-ray Tomographic Methods and the Effect of Shape on Concrete Rheology; Research Report ICAR 106-1; International Center for Aggregate Research: Austin, TX, USA, 2005. [Google Scholar]
Figure 1. Gradation of the AC-20 asphalt mixture.
Figure 1. Gradation of the AC-20 asphalt mixture.
Materials 12 03601 g001
Figure 2. (a) The asphalt mastic specimens; (b) the uniaxial static creep test testing process.
Figure 2. (a) The asphalt mastic specimens; (b) the uniaxial static creep test testing process.
Materials 12 03601 g002
Figure 3. (a) Graded coarse aggregates; (b) uniform-sized packed discrete elements.
Figure 3. (a) Graded coarse aggregates; (b) uniform-sized packed discrete elements.
Materials 12 03601 g003
Figure 4. Irregularly shaped polyhedron aggregates.
Figure 4. Irregularly shaped polyhedron aggregates.
Materials 12 03601 g004
Figure 5. (a) Discrete element (DE) model of asphalt mixture; (b) air voids within the DE model.
Figure 5. (a) Discrete element (DE) model of asphalt mixture; (b) air voids within the DE model.
Materials 12 03601 g005
Figure 6. Interactions among the components within the asphalt mixture.
Figure 6. Interactions among the components within the asphalt mixture.
Materials 12 03601 g006
Figure 7. Burger’s models. (a) Macro behavior; (b) micro behavior in the normal direction; (c) micro behavior in the shear direction.
Figure 7. Burger’s models. (a) Macro behavior; (b) micro behavior in the normal direction; (c) micro behavior in the shear direction.
Materials 12 03601 g007
Figure 8. Equivalent meso Burger’s model for contacts between aggregates and asphalt mastic.
Figure 8. Equivalent meso Burger’s model for contacts between aggregates and asphalt mastic.
Materials 12 03601 g008
Figure 9. Results of discrete element model (DEM) prediction and laboratory measurement.
Figure 9. Results of discrete element model (DEM) prediction and laboratory measurement.
Materials 12 03601 g009
Figure 10. Irregular polyhedral aggregates with different angularities. AI: angularity index.
Figure 10. Irregular polyhedral aggregates with different angularities. AI: angularity index.
Materials 12 03601 g010
Figure 11. Procedure for quantifying the discrete element number of aggregates: (a) quantifying the number of aggregate surfaces; (b) quantifying the number that makes up aggregates.
Figure 11. Procedure for quantifying the discrete element number of aggregates: (a) quantifying the number of aggregate surfaces; (b) quantifying the number that makes up aggregates.
Materials 12 03601 g011
Figure 12. Axial strains with different overall aggregate angularities.
Figure 12. Axial strains with different overall aggregate angularities.
Materials 12 03601 g012
Figure 13. Orientation of ellipsoidal aggregate.
Figure 13. Orientation of ellipsoidal aggregate.
Materials 12 03601 g013
Figure 14. Ellipsoidal aggregates with different orientations: (a) 0°; (b) 30°; (c) 45°.
Figure 14. Ellipsoidal aggregates with different orientations: (a) 0°; (b) 30°; (c) 45°.
Materials 12 03601 g014
Figure 15. Coarse aggregates with different aggregate orientations within asphalt mixtures: (a) 0°; (b) 30°; (c) 45°; (d) 60°; (e) 90°.
Figure 15. Coarse aggregates with different aggregate orientations within asphalt mixtures: (a) 0°; (b) 30°; (c) 45°; (d) 60°; (e) 90°.
Materials 12 03601 g015
Figure 16. Axial strains with different aggregate orientations.
Figure 16. Axial strains with different aggregate orientations.
Materials 12 03601 g016
Figure 17. Axial strains with different aggregate surface textures.
Figure 17. Axial strains with different aggregate surface textures.
Materials 12 03601 g017
Figure 18. Partitioning within asphalt mixture specimen: (a) partitioning in the horizontal direction; (b) partitioning in the vertical direction.
Figure 18. Partitioning within asphalt mixture specimen: (a) partitioning in the horizontal direction; (b) partitioning in the vertical direction.
Materials 12 03601 g018
Figure 19. Aggregates distributed non-uniformly within DE models of asphalt mixtures: (a) aggregates distributed non-uniformly in the horizontal direction; (b) aggregates distributed non-uniformly in the vertical direction.
Figure 19. Aggregates distributed non-uniformly within DE models of asphalt mixtures: (a) aggregates distributed non-uniformly in the horizontal direction; (b) aggregates distributed non-uniformly in the vertical direction.
Materials 12 03601 g019
Figure 20. Axial strains with different coefficient of variation of aggregate distribution.
Figure 20. Axial strains with different coefficient of variation of aggregate distribution.
Materials 12 03601 g020
Table 1. Volumetric properties of the AC-20 asphalt mixture.
Table 1. Volumetric properties of the AC-20 asphalt mixture.
Asphalt Content (%)Bulk Density (g/cm3)VV (%)VMA (%)VFA (%)
4.32.4324.013.872.2
Table 2. Experimental program for the asphalt mixture and mastics.
Table 2. Experimental program for the asphalt mixture and mastics.
MaterialsSpecimen SizeTestLoad/MPa
Asphalt mixtureΦ 100 mm × H 150 mm (cylinder)
L 100 mm × L 150 mm (cubic)
Uniaxial static creep test0.7
Asphalt masticΦ 100 mm × H 100 mm (cylinder)Uniaxial static creep test0.07
Table 3. Macro parameters for asphalt mastic and aggregates.
Table 3. Macro parameters for asphalt mastic and aggregates.
E1 (MPa)η1 (MPa·s)E2 (MPa)η2 (MPa·s)υE (GPa)μaυ
0.568973.1630.39627.8950.555.50.50.35
Table 4. Volume distribution percentages within different parts.
Table 4. Volume distribution percentages within different parts.
Model No.Volume Distribution Percentage within Each Part (%)
Part IPart IIPart III
Partitioned in the horizontal direction
#1303040
#2253045
#3203050
#4153055
Partitioned in the vertical direction
#5303040
#6253045
#7203050
#8153055
Table 5. Axial strains when aggregates were distributed non-uniformly.
Table 5. Axial strains when aggregates were distributed non-uniformly.
Model No.Coefficient of VariationAxial Strain ε (%)ε/ε0
Partitioned in the horizontal direction
#10.2830.151.12
#20.4950.161.24
#30.7070.201.56
#40.9190.282.15
Partitioned in the vertical direction
#50.2830.171.31
#60.4950.221.70
#70.7070.382.90
#80.9190.644.90

Share and Cite

MDPI and ACS Style

Zhang, D.; Gu, L.; Zhu, J. Effects of Aggregate Mesostructure on Permanent Deformation of Asphalt Mixture Using Three-Dimensional Discrete Element Modeling. Materials 2019, 12, 3601. https://doi.org/10.3390/ma12213601

AMA Style

Zhang D, Gu L, Zhu J. Effects of Aggregate Mesostructure on Permanent Deformation of Asphalt Mixture Using Three-Dimensional Discrete Element Modeling. Materials. 2019; 12(21):3601. https://doi.org/10.3390/ma12213601

Chicago/Turabian Style

Zhang, Deyu, Linhao Gu, and Junqing Zhu. 2019. "Effects of Aggregate Mesostructure on Permanent Deformation of Asphalt Mixture Using Three-Dimensional Discrete Element Modeling" Materials 12, no. 21: 3601. https://doi.org/10.3390/ma12213601

APA Style

Zhang, D., Gu, L., & Zhu, J. (2019). Effects of Aggregate Mesostructure on Permanent Deformation of Asphalt Mixture Using Three-Dimensional Discrete Element Modeling. Materials, 12(21), 3601. https://doi.org/10.3390/ma12213601

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop