Abstract

Asphalt mixture is susceptible to moisture damage under the effect of freeze-thaw (F-T) cycles. In this paper, crumb rubber (CR) was used to modify stone mastic asphalt (SMA) and the effects of diatomite and styrene butadiene styrene (SBS) on antifreezing performances of crumb rubber modified SMA (CRSMA) were investigated. Regression analysis and modified grey model (MGM) were used to construct the prediction models for properties of modified mixtures. CRSMA, CR and diatomite modified SMA (CRDSMA), and CR and SBS modified SMA (CRSSMA) were prepared in laboratory, respectively. Process of F-T cycles was designed. Air void, indirect tensile strength (ITS), and indirect tensile stiffness modulus (ITSM) were measured to evaluate the antifreezing performances of CRSMA, CRDSMA, and CRSSMA. Results indicate that air voids increase with the increasing of F-T cycles. ITS and ITSM all decrease with the increasing of F-T cycles. The addition of diatomite and SBS can reduce the air void and improve the ITS and ITSM of CRSMA. CRSSMA presents the lowest air void, highest tensile strength, and largest stiffness modulus, which reveals that CRSSMA has the best F-T resistance among three different kinds of mixtures. Moreover, MGM (1, 2) models present more favorable accuracy in prediction of air void and ITS compared with regression ones.

1. Introduction

Reuse of waste materials can reduce the consumption of natural resources. Moreover, it can alleviate the waste accumulation, which can cause serious environmental problem [13]. Therefore, reuse of waste materials has got great attention. With the rapid development of automotive industry, the number of vehicles increases, and then amount of scrap tires also increases. According to the statistical data, over 1 billion tires are sold each year all around the world and millions of waste tires are generated [4, 5]. Waste tires are not easily disposed of and lead to great threat to public health and environment. Many solutions have been proposed for reuse of scrap tires [69].

Scrap tires can be recycled in the form of crumb rubber (CR), which is a profitable practice. La Rosa et al. [9] demonstrated that the addition of 50 phr of ground tire rubber (GTR) in styrene-isoprene-styrene (SIS) formulation presents favorable environmental benefits. Life cycle assessment (LCA) results indicate that the effects of virgin rubber (47.1%) and carbon black (34.3%) on environmental impacts of GTR are more obvious than tire grinding process (4%). Farina et al. [10] investigated the environmental performance of using crumb rubber in bituminous mixtures based on LCA. Results reveal that asphalt rubber obtained through wet process presents significant benefits for energy saving, environmental impact, human health, preservation of ecosystems, and minimization of resource depletion. The reductions of gross energy requirement and global warming potential range between 36% and 45% compared with standard paving solutions. Furthermore, CR can be used to modify asphalt, which is an effective method to consume them. Previous researches have demonstrated that CR is a successful modifier for asphalt mixture because of favorable compatibility and interaction between rubber and asphalt [5]. CR modified asphalt and mixture have many advantages. They can improve the rutting resistance at high temperature, crack resistance at low temperature, fatigue resistance, and elastic performance [4, 1113]. Moreover, Liang et al. [14] investigated the viscous properties, storage stability, and morphology of recycled tire scrap rubber modified asphalt (TSRMA). The results reveal that the addition of crumb rubber can increase viscosity of asphalt and improve rutting deformation resistance. However, it presents poor stability with the increasing of rubber particle size. Navarro and Gámez [4] presented the effect of crumb rubber on bearing capacity and cohesion of asphalt mixture. The results show that crumb rubber can increase the stiffness and stability of mixture, while it slightly reduces tensile strength. De Almeida Júnior et al. [15] evaluated the properties of asphalt modified by rubber and styrene butadiene styrene (SBS), respectively. The findings indicate that rubber has the potential performance to substitute SBS in pavement materials. However, it does not meet some of the requirements in standard specifications.

Diatomite is a widely used mineral with low cost and considerable storage, which has high absorptive capacity and stability [1622]. It has been widely used for asphalt modification in China [1621]. Mohamed Abdullah et al. [17, 18] studied the properties of diatomite modified asphalt. The results indicate that there is no chemical reaction between asphalt and diatomite. Diatomite can improve the rutting resistance of asphalt at high temperature. The aging resistances of modified asphalt are superior to pure one. Cheng et al. [19] demonstrated the antiaging properties of diatomite modified asphalt. The results reveal that diatomite is helpful for improving the high temperature stability and antiaging property of asphalt. Tan et al. [21] investigated the low temperature property of diatomite modified asphalt. The results show that diatomite particles distribute well in asphalt and diatomite modified mixture possesses better low temperature performance than neat one.

Currently, SBS copolymer is the most widely used modifier for asphalt, which can absorb the oil fractions and swell them [2325]. It can improve both high and low temperature performances of mixture. Khodaii and Mehrara [26] conducted dynamic creep tests for unmodified and SBS modified asphalt. SBS modified mixture possesses lower temperature susceptibility and higher compactibility. Sun et al. [27] found that reactive blending of SBS and asphalt can improve the hot storage stability, rheological properties, elasticity, and deformation resistance of modified asphalt. The limitations for SBS modified asphalt lie in its lower workability, higher production cost, and being easy to aging [28, 29]. Therefore, it is crucial to have some other waste materials to partially substitute for SBS in modified asphalt. CR can reduce the cost and possess favorable performances, which is a favorable substitute. In order to meet the requirements of high grade road pavement, modified asphalt with better performances can be prepared using CR and SBS composite modifier. Dong et al. [30] investigated the rheological behavior of SBS/CR composite modified asphalt. The results reveal that the addition of CR and SBS can effectively improve the viscoelastic property and temperature sensitivity of asphalt.

In addition, freeze-thaw (F-T) cycle is an important factor to influence the performances of road pavement in seasonal frozen regions, which cover 53.5% of China’s land area. Road pavement in this region is exposed to at least one F-T cycle each year [31]. CR has been widely used in concrete to provide F-T protection [32]. Meanwhile, Zhang et al. [33] found that rubber modified asphalt mixture has better antifreezing performance than the neat one. Current research results indicate that CR is an effective additive to improve the F-T resistance for road materials.

Grey theory was proposed by Deng in 1982 [34], which had been widely used in many fields for factor analysis and forecasting. Grey relational analysis (GRA) was widely used to evaluate the influential degrees of factors [35, 36]. Grey model is another important aspect in grey theory, which can achieve high prediction accuracy for uncertain systems. It can realize the prediction of system behavior using fewer data and less complicated mathematical calculation [37]. Ren et al. [38] predicted the long-term cumulative plastic deformation of subgrade under cyclic loads using grey model. Wang and Li [39] adopted grey model for pavement smoothness prediction. Grey model was also used for the predictions in environmental vulnerability assessment [40], natural disasters [41] and road traffic [42], and so forth. GM (1, 1) (Grey Model First-Order One Variable) is one of the most widely used grey models in practical application [43]. However, the limitation existing in traditional GM (1, 1) is that it needs equal interval for time series data. It cannot be used for time series data with nonequal interval [44].

As discussed above, reuse of CR and diatomite possesses favorable economic and environmental benefits. Furthermore, F-T action has significant effect on the properties of asphalt pavement in seasonal frozen regions. It will reduce the pavement life severely. The compound modified effects of CR and diatomite, CR, and SBS on asphalt mixture under F-T actions are still unknown. In order to improve the freezing resistance of asphalt pavement, crumb rubber modified stone mastic asphalt (CRSMA), crumb rubber and diatomite modified stone mastic asphalt (CRDSMA) and crumb rubber and SBS modified stone mastic asphalt (CRSSMA) were prepared, respectively. Effects of F-T cycles on three different kinds of asphalt mixtures were tested and evaluated using air voids, indirect tensile strength (ITS) test, and indirect tensile stiffness modulus (ITSM) test. Meanwhile, prediction analysis was used to construct the relationships between properties and F-T cycles in order to improve the convenience for performance analysis of mixture. Modified grey model MGM (1, 2) was applied for prediction of air void and ITS.

2. Materials

2.1. Raw Materials

Base asphalt AH-90 and SBS modified asphalt used for experiments were from Panjin Petrochemical Industry, Liaoning Province of China. Physical properties of AH-90 asphalt were measured and listed in Table 1. Properties of SBS modified asphalt were tested and given in Table 2. Aggregate and mineral filler were obtained from Changchun Municipal Mixing Plant in Jilin Province of China. Physical properties of andesite mineral aggregate were given in Table 3. Limestone powder was chosen as mineral filler. Corresponding properties of filler were listed in Table 4. Crumb rubber particle was obtained from Changchun Yuxing Rubber Materials Co., Ltd., Jilin Province, China. Its properties were shown in Table 5. Diatomite was produced by Changchun Diatomite Products Co., Ltd., Jilin Province, China. Its physical properties and particle distribution were listed in Tables 6 and 7.

2.2. Asphalt Mixture

Stone mastic asphalt (SMA) was developed in Germany during the 1960s. It has been used as the premium pavement surfacing course for heavy-duty pavements, high-speed motorways, and highways in America, Europe, China, and several other countries due to its excellent rutting resistance, fatigue resistance, durability, and lower noise characteristics compared with dense graded asphalt [45, 46]. In this research, SMA-13 gradation was designed based on the Course Aggregate Void Filling (CAVF) method [47]. The designed gradation was shown in Figure 1. According to previous research results, the diatomite content 15% by weight of asphalt and crumb rubber content 3% by weight of aggregate were adopted [15]. Optimum asphalt contents were determined by Marshall method. They are 6.3%, 6.5%, and 6.6% for CRSMA, CRDSMA, and CRSSMA, respectively.

In order to make crumb rubber and diatomite disperse homogeneous in the mixture, the preparation process was optimized. The designed preparation process for CRSMA was shown in Figure 2. Similar to the process in Figure 2, CRDSMA would be produced if filler was replaced with filler and diatomite; CRSSMA was obtained if neat asphalt was replaced with SBS modified asphalt. Detailed steps were explained as follows.

Step 1. Crumb rubber and aggregate were heated for 4 h at 180°C, respectively. They were blended together using mixing pot (BH-10, Beijing Zhonghangkegong Instrument Co., Ltd., China). The blending temperature is 165°C, speed is 80 r/min, and time is 90 s. This step is identical for CRSMA, CRDSMA, and CRSSMA.

Step 2. Asphalt was heated to 165°C and added to the mixture of crumb rubber and aggregate. They were blended together at 165°C for 90 s using mixing pot. The mixing speed is also set to be 80 r/min. In this step, neat asphalt was used for CRSMA and CRDSMA, while it was SBS modified asphalt for CRSSMA.

Step 3. Filler was heated for 4 h at 180°C and added to the mixture consisting of crumb rubber, aggregate, and asphalt (neat one for CRSMA and CRDSMA and SBS modified one for CRSSMA). They were blended together at the speed of 80 r/min and 165°C for 90 s using mixing pot. In this step, mineral powder was used for CRSMA and CRSSMA, while mineral powder and diatomite were used for CRDSMA.
Then, CRSMA, CRDSMA, and CRSSMA were prepared for testing. Marshall samples (101 mm diameter and 63.5 mm height) were prepared by Marshall hammer method. For each sample, its weight is 1200 g. According to the determined mixture proportion, aggregate 984.29 g, mineral powder 109.37 g, crumb rubber 32.81 g, and neat asphalt 73.53 g were used for preparation of each CRSMA sample, while aggregate 973.02 g, mineral powder 108.11 g, crumb rubber 32.43 g, neat asphalt 75.16 g, and diatomite 11.28 g were used for each CRDSMA sample and aggregate 981.23 g, mineral powder 109.02 g, crumb rubber 32.71 g, and SBS modified asphalt 77.04 g were used for each CRSSMA sample. The height of each sample was controlled to be () mm.

3. Experimental and Analytical Methods

3.1. Process of Freeze-Thaw

In this research, the detailed processes of freeze-thaw were designed and listed as follows:(1)The samples were divided into several groups. Each group consisted of three samples. One group of samples was placed at room temperature. Other groups of samples were immersed into water and placed in vacuum drying box for 15 min under vacuum degree 98.0 kPa. Then, vacuum drying box was adjusted to atmospheric condition. The samples were immersed into water for another 0.5 h.(2)Each sample was taken out of the tank and placed into plastic bag with 10 ml water. The bag was put in refrigerator for 16 h at −20°C until all the water in plastic bag was frozen.(3)The samples were removed from plastic bags and placed into water tank for 8 h with constant temperature 60°C until all the water was melted.

After (1), (2), and (3), one F-T cycle was completed. In this study, the number of F-T cycles was 15 for ITS test and 10 for air void and ITSM tests. After F-T cycles, the samples were removed from plastic bags and placed into water tank for 24 h with constant temperature 60°C. Then, they were put into water tank for 2 h with constant temperature 25°C and sample spacing was bigger than 10 mm. The treated samples were used for air void, ITS, and ITSM tests.

3.2. Air Void Test

In this paper, air voids for samples were tested by surface-dry condition method according to standard AASHTO T-166 [48]. Firstly, Marshall specimens (101 mm diameter and 63.5 mm height) were prepared through compaction method with 75 blows on each side of cylindrical samples. Secondly, masses for samples in dry, water, and surface-dry conditions were measured, respectively. Bulk specific gravity for samples can be calculated bywhere , , and are masses of sample in dry, water, and surface-dry conditions, respectively. Then, the theoretical maximum specific gravity () of samples was tested through vacuum sealing method. Finally, air voids of samples can be obtained by

Change ratio of air void was calculated bywhere is the air void of mixture after times of F-T cycles, .

3.3. Indirect Tensile Strength (ITS) Test

ITS test is an effective measure to evaluate the anticracking ability of asphalt mixture under low temperature [4, 21]. In this research, the test was carried out according to standard AASHTO T-283 at temperature of −10°C [48]. It is performed using a material testing machine (MQS-2) produced by Nanjing Ningxi Instrument Co., Ltd., Jiangsu Province, China. Loading rate was 1 mm/min. Marshall specimens (101 mm diameter and 63.5 mm height) were placed in the chamber at −10°C for 6 h before tests. The failure loads were recorded at the end of tests. ITS can be calculated bywhere ITS is the tensile strength of specimen; is the failure load; is the diameter of specimen; is the height of specimen.

In order to evaluate the F-T resistance ability of CRSMA, CRDSMA, and CRSSMA, change ratio of ITS before and after F-T cycles is given bywhere is the change ratio of ITS after times of F-T cycles. is ITS without F-T cycles, and is ITS after times of F-T cycles.

3.4. Indirect Tensile Stiffness Modulus (ITSM) Test

In order to investigate the bearing capacity of asphalt mixture, ITSM test was performed according to standard AASHTO TP-31 [48]. ITSM test was widely used to analyze the tensile properties of mixture, which can be further related to the cracking performance of pavement. In this research, a multifunctional testing machine (NU-14) was used, which was produced by Cooper Research Technology, Ltd., United Kingdom. In order to investigate the antifreezing performance of modified asphalt mixtures, the tests were conducted at 5°C and Marshall specimens (101 mm diameter and 63.5 mm height) were adopted. The specimens were placed in the chamber at the test temperature for 6 h before tests. The schematic diagram [4] for load pulse was shown in Figure 3. Rise time was 124 ms and duration period was 3.0 s. The target deformation in horizontal direction was 5 μm. The peak load was adjusted according to the test value of target deformation during the test. The stiffness modulus can be obtained bywhere is the stiffness modulus, MPa; is the peak load, N; is Poisson ratio, which is 0.25 at 5°C; is the height of specimen, mm; is the measured deformation in horizontal direction, mm.

Change ratio of before and after F-T cycles was calculated bywhere is the change ratio of stiffness modulus after times of F-T cycles, is stiffness modulus without F-T cycles, and is stiffness modulus after times of F-T cycles.

3.5. Modified Grey Method MGM

Traditional GM (1, 1) can realize prediction of time series data with equal interval. It cannot be used for time series data with nonequal interval. In this research, time series data is the number of F-T cycles, which is 0, 1, 3, 6, 10, and 15 with nonequal intervals. Therefore, the modified grey model (MGM) is used to construct the prediction model of air void and . The prediction procedure using MGM (Modified Grey Model First-Order m Variables) is briefly introduced in this paper and listed in Appendix.

4. Results and Discussion

4.1. Air Voids

Air void is an important characteristic of internal structure that affects the moisture damage of asphalt mixture. Air voids of CRSMA, CRDSMA, and CRSSMA under times of F-T cycles were tested, respectively. Corresponding () were calculated. Three specimens were measured for each case. For these measurements, mean value and standard deviation (SD) were calculated, respectively. Corresponding estimation of vertical error bar was determined for each case, which represents the standard deviation. Results of air voids (Mean ± SD) and for CRSMA, CRDSMA, and CRSSMA before and after F-T cycles are shown in Figures 4 and 5.

As can be seen from Figure 4, there is no significant increase or decrease in air void amplitude for each case. Air voids of CRSMA, CRDSMA, and CRSSMA are all below 7% after 10 times of F-T actions. Air voids increase with the increasing of F-T cycles, which agrees with the result obtained by Xu et al. [49]. Air void of CRSMA increases from 4.26% to 6.71% after 10 times of F-T cycles compared with that without F-T effect. It increases from 4.08% to 6.35% for CRDSMA and 4.09% to 6.01% for CRSSMA. The reasons lie in that F-T cycles can lead to not only expansion of existing voids but also formation of new voids [49]. Furthermore, air voids of CRSMA are generally larger than CRDSMA or CRSSMA, while air voids of CRSSMA are the lowest. According to the result by Song et al. [20], diatomite has larger void and absorption properties than mineral filler. It can effectively absorb the lower molecular group and lower polar aromatic molecules of asphalt to form thicker asphalt film and anchorage structure around aggregates and crumb rubber particles, which can increase the contact area and decrease the air voids between aggregates and rubber particles. According to the results of SBS modified asphalt by Larsen et al. [23] and Adedeji et al. [24], polymer phase and asphalt phase exist together. Polymer globules are dispersed homogenously in the continuous asphalt phase, which form a stable network structure. It can strengthen the adhesion between aggregate and rubber particle and decrease the air voids.

As shown in Figure 5, are the largest for CRSMA and CRDSMA. It reveals that rapid expansion of existing voids and formation of new voids happen after the first F-T cycle. However, change ratio of air void for CRSSMA increases before 5 times of F-T cycles. It is induced by the strong network structure formed by two twisted continuous phases [22], which delays the growth of air void. Change ratios of air voids for CRSMA, CRDSMA, and CRSSMA tend to be stable after several times of F-T actions. It is because the expansion of existing voids and formation of new voids come to steady state.

The analysis of variance (ANOVA) is the most commonly used statistical method to determine the significant effects of factors [50]. ANOVA and F-T tests were conducted to determine statistically significant process parameters of F-T cycles on air voids for CRSMA, CRDSMA, and CRSSMA. ANOVA results were listed in Table 8.

According to ANOVA results, values of CRSMA, CRDSMA, and CRSSMA are all greater than . The probability is 99% that F-T cycle possesses significant effects on air voids of CRSMA, CRDSMA, and CRSSMA. Besides, value of CRSMA is the largest, while value of CRSSMA is the smallest. The results indicate that F-T cycle presents the most statistically significant effect on air void of CRSMA and lowest one on CRSSMA. The additions of SBS and diatomite can reduce the influence of F-T cycles on air void.

4.2. ITS Results

and were calculated for CRSMA, CRDSMA, and CRSSMA after F-T cycles 0, 1, 3, 6, 10, and 15. Three specimens were measured for each case. For these measurements, mean value and standard deviation (SD) were calculated, respectively. Corresponding estimation of vertical error bar was determined for each case, which represents the standard deviation. Results of (Mean ± SD) and are shown in Figures 6 and 7.

As can be seen from Figures 6 and 7, ITS for three different kinds of asphalt mixtures all decrease with the increasing of F-T cycles and decrease rates become smaller. Meanwhile, of CRSMA, CRDSMA, and CRSSMA increase with the increasing of F-T cycles. The reasons lie in that repetitive F-T actions increase the air voids, which is conducive to water flow. Effect of water frost will cause microcrack formation and propagation in asphalt mixture. These kinds of serious damage lead to the reduction of ITS. With the increasing of F-T cycles, air voids tend to be stable. Therefore, the frost effect of water slows down. Moreover, ITS of CRSSMA is the highest, while ITS of CRSMA is the lowest. Slowing trends for CRDSMA and CRSSMA are more obvious than CRSMA. With the increasing of F-T cycles, of CRSSMA is the smallest, while it is the biggest for CRSMA. The reasons lie in that diatomite has good absorption ability of asphalt, which makes the adhesion capability of asphalt mortar increase and also the asphalt film thickness on aggregate surface [20, 35]. Therefore, the cohesion among aggregates is strengthened and the anti-F-T ability of mixture improves, which enhances the resistance to internal damage of mixture. For SBS, its elastomeric phase can absorb the oil fraction of asphalt and swell it. A continuous polymer phase forms that can improve the performances of asphalt [23, 24]. The network structure of SBS modified asphalt is more stable than diatomite modified one.

4.3. ITSM Results

of CRSMA, CRDSMA, and CRSSMA were tested and corresponding were calculated. Three specimens were measured for each case and the mean value was regarded as the result. Results for and are listed in Table 9.

As shown in Table 9, for CRSMA, CRDSMA, and CRSSMA decrease with the increasing of F-T cycles. increase with the increasing of F-T cycles. The reasons lie in that stiffness modulus is used to evaluate the bearing capacity of mixtures. Formation of microcrack under F-T action changes the internal structure and weakens its integrity of asphalt mixture [49], which reduces the bearing capacity. With the increasing of F-T cycles, the influence of internal damage is more obvious. Furthermore, CRSSMA possesses the highest stiffness modulus value, while CRSMA has the lowest one. It is because diatomite modified asphalt and SBS modified one have higher cohesion ability than the neat one [23], which can lower the degree of internal damage and also the bearing capacity for mixture. Similarly, network structure of SBS modified asphalt presents superior performance than diatomite modified one.

and ITS are important characteristics to represent the mechanical properties of asphalt mixture. Their correlations need to be demonstrated, which are shown in Figure 8. It can be observed that and ITS present favorable linear relationships. It can provide reference for the evaluation of mechanical properties of mixtures.

4.4. Performance Prediction

In practice, processes of F-T, air void, ITS, and ITSM tests are complicated and time-consuming. They make the evaluation of mixture properties inconvenient. Prediction analysis can be used to construct the relationships between mixture properties and F-T cycles, which is convenient for performance analysis of mixture. In this paper, prediction models of air void and ITS for CRSMA, CRDSMA, and CRSSMA were established based on regression analysis and MGM (1, 2), respectively. The prediction effects of two models were compared and analyzed.

4.4.1. Regression Analysis Model

The regression analysis has been one of the most popular methods for performance prediction in many fields from decades [5153]. Regression analysis consists of fitting a function to the data to discover how one or more variables (dependent variable) vary as a function of another (independent variable). In this study, independent variable is F-T cycles; dependent variables are air void and , respectively.

For CRSMA, is signed as variable and air void is . For CRDSMA, is signed as variable and air void is . For CRSSMA, is signed as variable and air void is . Regression models for CRSMA, CRDSMA, and CRSSMA are established based on the results of air voids and after 0, 1, 3, 6, and 10 times of F-T cycles. They are shown in Figures 9 and 10.

As can be seen from Figures 9 and 10, two order polynomials were used to construct the prediction models of for CRSMA, CRDSMA, and CRSSMA. values of the models for were higher than 0.97. Linear regression models were adopted for the prediction of air void for CRSMA and CRSSMA, and two order polynomials were used to for CRDSMA. values of the models for air void were higher than 0.95. All these models revealed close agreements between experimental results and predicted ones.

The prediction results of (, , ) and air void (, , ) are calculated by use of the established regression models. The relative errors are obtained for and air void after 15 times of F-T cycles, which can be calculated bywhere is relative error between test value and predicted one; and are test value and predicted one after th F-T cycles, respectively.

Corresponding root-mean-square error (RMSE) is also calculated, which is a frequently used measure of the differences between values predicted by a model and the values actually tested. The smaller RMSE is, the higher the prediction accuracy of model presents. RMSE can be calculated bywhere is the number of samples.

The prediction results of regression models are listed in Tables 1012.

As shown in Tables 10, 11, and 12, the relative errors for and after 15 times of F-T cycles are 1.15% and 3.09%, respectively. They are 2.98% and 19.94% for and and 2.00% and 5.08% for and . The prediction accuracy of is close to 20%, which is unsatisfactory. As for the maximum relative error, they are 2.75%, 7.75%, 22.55%, 19.94%, 11.08%, and 5.08% for , , , , , and , respectively. They are more than 5% except . In the aspect of RMSE, they are 0.78%, 0.18%, 1.96%, 0.58%, 1.23%, and 0.17% for , , , , , and , respectively.

4.4.2. MGM (1, 2) Model

For CRSMA, is signed as variable and air void is . MGM (1, 2) models for CRSMA are established based on the results of air voids and after 0, 1, 3, 6, and 10 times of F-T cycles.

For MGM (1, 2) models of CRSMA, the first-order differential equation is

Time response functions for and air void after 15 times of F-T cycles can be obtained as

The restored values for and air void after 15 times of F-T cycles can be calculated by

The prediction results of () and air void () are calculated by use of the established MGM (1, 2) model. The relative errors and RMSE are calculated for and air void after 15 times of F-T cycles. The prediction results of MGM (1, 2) model are listed in Table 13.

Prediction models of air voids and for CRDSMA and CRSSMA are also established based on MGM (1, 2) according to the calculation procedure of CRSMA. For CRDSMA, is signed as variable and air void is . The prediction results of () and air void () are calculated by use of the established MGM (1, 2) model. The relative errors and RMSE are calculated for and air void after 15 times of F-T cycles. The prediction results of MGM (1, 2) model are listed in Table 14.

For CRSSMA, is signed as variable and air void is . The prediction results of () and air void () are calculated by use of the established MGM (1, 2) model. The relative errors and RMSE are calculated for and air void after 15 times of F-T cycles. The prediction results of MGM (1, 2) model are listed in Table 15.

As shown in Tables 13, 14, and 15, the relative errors for and after 15 times of F-T cycles are 3.77% and 0.90%, respectively. They are 4.26% and 3.35% for and and 2.10% and 1.74% for and . It reveals that the established MGM (1, 2) models possess favorable accuracy in performance prediction of CRSMA, CRDSMA and CRSSMA. As for the maximum relative error, they are 3.77%, 0.90%, 5.64%, 3.35%, 3.82%, and 2.84% for , , , , , and , respectively. They are less than 5% except . In the aspect of RMSE, they are 1.01%, 0.03%, 0.91%, 0.11%, 0.61%, and 0.09% for , , , , , and , respectively.

Compared with regression models, the relative errors of MGM (1, 2) models in prediction of air void and after 15 times of F-T cycles are more stable and reliable. The maximum relative errors and RMSE of MGM (1, 2) models are much lower than regression ones except for of CRSMA. The results indicate that the differences between values predicted by MGM (1, 2) model and the values actually tested are smaller than regression model. MGM (1, 2) model presents more favorable accuracy in performance prediction of CRSMA, CRDSMA, and CRSSMA.

5. Conclusions

Reuse of waste materials is conducive to the protection of natural resources and environment. In this paper, CR and diatomite were recycled to modify SMA. CR, diatomite, and SBS were used to prepare CRSMA, CRDSMA, and CRSSMA. F-T effects on these mixtures were investigated. The following conclusions can be obtained:(1)Air voids of CRSMA, CRDSMA, and CRSSMA all increase with the increasing of F-T cycles. The addition of diatomite into CRSMA can reduce its air void because diatomite can effectively absorb the lower molecular group and lower polar aromatic molecules of asphalt to form thicker asphalt film and anchorage structure around aggregates and crumb rubber particles. As for the effect of SBS, it can also reduce the air void of CRSMA because of the strong network structure of SBS modified asphalt. CRSSMA presents the lowest air voids during the test. ANOVA results reveal that F-T cycle presents the lowest statistically significant effect on air void of CRSSMA.(2)Indirect tensile strength for three different kinds of asphalt mixtures all decrease with the increasing of F-T cycles. Diatomite and SBS present enhancement effects for indirect tensile strength of CRSMA. Anchorage structure of diatomite modified asphalt and continuous polymer phase of SBS modified asphalt can improve the cohesion among aggregates and enhances the resistance of mixture to internal damage. Results for change ratio of tensile strength before and after F-T cycles reveal that CRSSMA is the smallest and CRSMA is the biggest.(3)Indirect tensile stiffness modulus of CRSMA, CRDSMA, and CRSSMA decreases with the increasing of F-T cycles. Diatomite and SBS can effectively improve the stiffness modulus of CRSMA, which is because diatomite modified asphalt and SBS modified one have higher cohesion abilities than the neat one. CRSSMA possesses the highest stiffness modulus value. Stiffness modulus and indirect tensile strength present favorable linear relationship.(4)The differences between values predicted by MGM (1, 2) model and the values actually tested are smaller than regression model. It indicates that MGM model presents more favorable accuracy and can be used for performance prediction of asphalt mixtures.

In summary, the network structure of SBS modified asphalt is more stable than diatomite modified one. Pavement constructed using CRSSMA will have the better antifreezing performance than that using CRSMA or CRDSMA. However, the pavement constructed using CRDSMA will be a favorable choice considering its environmental and economic advantages. It not only presents good antifreezing performance but also realizes the reuse of waste materials including CR and diatomite.

Appendix

Modified Grey Method MGM

The prediction procedure using MGM (Modified Grey Model First-Order Variables) is briefly introduced as follows:

Step 1. Assume that the original series of data with variables are where is the nonnegative original time series data; is the number of variables. () is the th variable with nonequal interval; is the number of time series data for variables. is the data at time .

Step 2. One time accumulating generation operational sequences (1-AGO sequences) for are where is the AGO sequence of ; , is the interval between and .

Step 3. MGM () model can be obtained by the following first-order differential equation:Assume , .
Then, (A.3) can be expressed bywhere .

Step 4. The solution of (A.4) is which is called the time response function. In (A.5), matrix and vector are undetermined.

Step 5. Solution of and . The generating sequence for nearest-neighbor mean value with nonequal interval is where .

Difference equation for (A.5) can be obtained after discretization; it is

Least squares estimation sequence for MGM () is where

Matrix and vector can be calculated by

Step 6. Time response function of difference equation is

Step 7. The restored values can be obtained as follows:

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors express their appreciations for the financial supports of National Natural Science Foundation of China under Grants nos. 51408258 and 51578263; China Postdoctoral Science Foundation Funded Project (nos. 2014M560237 and 2015T80305); Fundamental Research Funds for the Central Universities (JCKYQKJC06), Transportation Science and Technology Project of Jilin Province (2015-1-11), and Science & Technology Development Program of Jilin Province.