1. Introduction
Fatigue cracking is one of the most common distresses of asphalt pavement. Crack initiation and propagation can reduce the bearing capacity of the pavement structure and allow the penetration of water into the structure, resulting in other distresses such as water damage [
1]. Fatigue cracking is mainly caused by repeated traffic loading. The repeated loading can lead to the micro-cracks forming and merging, eventually propagating into macro-cracks, i.e., alligator cracking, which can weaken the overall structural capacity of the asphalt pavement and deteriorate the bonding between asphalt and aggregates [
2]. Therefore, it is of significance to evaluate the fatigue behavior of the asphalt mixture, which could help prolong the fatigue life.
Most of the existing studies have primarily analyzed the fatigue behavior of the asphalt mixture using laboratory tests and summarized the mechanism of fatigue cracking [
3,
4,
5]. However, most of the previous studies focused on the properties of lab-fabricated asphalt mixture specimens and were limited when evaluating the properties of the core samples, which can lead to inefficient utilization of the core samples. Moreover, in practical pavement applications, asphalt pavement is constructed through three phases of compaction and ages under in-service conditions, which are subject to environmental fluctuations, such as temperature variation, precipitation, and ultraviolet rays. Thus, there is much difference in the internal structural distribution and aging between the lab-fabricated asphalt mixture specimens and those in the field. Consequently, many researchers have shifted focus to the asphalt field cores, studying the fracture properties of the cores drilled from the in-service pavement over various durations.
There are many lab tests to study the fatigue properties of the asphalt mixture, such as the indirect tensile (IDT) test, direct tensile (DT) test, overlay test (OT), four-point bending test, and semi-circular bending (SCB) test. Barman et al. conducted the IDT test to characterize the fatigue resistance of the asphalt mixture and proposed a simple data analysis approach [
6]. Luo et al. used the controlled-strain repeated direct tensile test to evaluate the fatigue cracking [
7]. Gu et al. conducted the overlay test to investigate the fracture properties of the field-aged asphalt concrete and found that the cracking resistance of the field reduced from 1st month to the 9th month [
8]. Kim et al. performed a four-point bending test to investigate the fatigue life of a total of ten asphalt mixtures, including hot-mix asphalt (HMA) and warm-mix asphalt (WMA) with different amounts of reclaimed asphalt pavement (RAP) and recycled asphalt shingles (RAS) [
9]. Du et al. performed the SCB test on the layer core samples drilled from five expressways to analyze the sensitivity of fracture energy to factors such as the equivalent single axle load, air void, service age, etc. [
10]. These studies demonstrate that the field core samples can reflect the asphalt pavement conditions, and the test results of core samples can be effectively used for the decision-making related to pavement maintenance actions.
Based on the results of the fatigue test of the asphalt mixture, the fatigue cracking models of the asphalt mixtures were developed. The fracture mechanics and the dissipated energy approach are most widely used to evaluate the fatigue resistance of the asphalt mixtures [
11]. However, the fatigue test results cannot describe the crack propagation of asphalt mixture at the mesoscopic level. Moreover, it is hard to validate the fatigue cracking models with only a limited number of core samples.
In recent years, researchers have attempted to use computer technology to simulate the fracture evolution of the asphalt mixture and investigate the various influencing factors on the fatigue behavior. The discrete element method (DEM) has been widely used in pavement engineering since discrete elements can reflect the discontinuous and non-uniform structural characteristics of asphalt mixtures. It can also help reveal the internal structural deformation, cracking, and other mechanical behaviors of asphalt mixtures. Ma et al. built a virtual specimen based on the DEM to estimate the fatigue life of an asphalt mixture and investigated the influence of air void on fatigue life [
12]. Xue et al. developed a new approach combining algorithmic techniques and DEM to perform a heterogeneous fracture simulation, and the study proved that the DEM could provide a valid understanding of the fracture behavior of materials so as to be used to diminish the need for numerous laboratory tests [
13]. Peng et al. adopted Python language and DEM to generate irregular particles and establish a three-dimensional (3D) discrete element model of asphalt surface to study the mechanical response under different working conditions [
14]. However, due to the limitation of obtaining the raw material parameters, most of the existing studies depended on lab-fabricated specimens and were limited to discrete element simulation of the fractures of the core samples. Moreover, it is unclear whether the calibration of the mesoscopic parameters for the simulation through the lab tests can be effectively applied to the limited number of core samples.
The objective of this study is to propose a reliable method to evaluate the fatigue behavior of the asphalt field cores based on discrete element modeling and to conduct mesoscopic contact parameter calibration through lab tests, including the uniaxial compression dynamic modulus test, SCB test, and SCB fatigue test, which can enhance the utilization efficiency of core samples and provide a reliable representation of the fatigue behavior of the core samples.
This paper is organized as follows. The following section presents the test samples and lab tests conducted to determine the mesoscopic contact parameters in the discrete element modeling. The next section describes the establishment of the virtual specimen using CT scanning of the asphalt field cores and image processing technologies, as well as mesoscopic contact parameter calibration and virtual fatigue tests using discrete element modeling. The fatigue life obtained from the lab test and virtual fatigue test is compared, and the virtual test results of force chains, crack evolution, and displacement are discussed in the following section. The final section summarizes the findings of this study.
5. Conclusions
In this study, the virtual SCB fatigue test was simulated by using discrete element modeling to evaluate the fatigue behavior of the asphalt field cores. The CT scan test was conducted to build the mesoscopic geometric model of the asphalt field cores. Additionally, the uniaxial compression dynamic modulus test and SCB test were performed to determine the parameters of the contact model in the virtual fatigue test. Based on the virtual SCB fatigue test, the displacement and contact forces, as well as crack growth, were analyzed. The main findings of this study can be drawn as follows.
- (1)
The evaluation methodology of fatigue behavior of the asphalt field cores based on the discrete element simulation was developed and can be used to enhance the effective usage of the field cores, which can help with the decision-making of pavement maintenance actions.
- (2)
The fatigue life simulated in the virtual fatigue test was consistent with that of the laboratory SCB fatigue test. The error between the simulated and test fatigue life was below 20%, which shows that the virtual fatigue test result is acceptable and reliable.
- (3)
It was found from the analysis of the force chain evolution process that concentrated tensile stress during loading can lead to crack initiation and propagation, ultimately resulting in material failure.
- (4)
The fatigue cracks in the asphalt mixture were observed as the three stages, i.e., crack initiation, crack propagation, and failure. It was found that the crack propagation stage consumes a significant portion of the fatigue life since tensile contact force mainly increases in this stage.
In this study, the discrete element modeling was restricted within 2D simulation due to the limited computational power. In future work, the 3D discrete element simulation will be performed to evaluate the fatigue behavior of the asphalt field cores, which could further improve the simulation accuracy. Additionally, future work will compare the asphalt field cores with the different aging times and investigate the difference in the fatigue behavior among them.