The Multi-Objective Optimization Design and Hydrothermal Performance Evaluation of Anhydrous Calcium Sulfate Whisker and Polyester Fiber Compound Modified Asphalt Mixture in Hot-Humid Areas

In hot and humid climates, asphalt pavements frequently encounter environmental factors such as elevated temperatures and rainfall, leading to rutting deformations and potholes, which can affect pavement performance. The primary objective of this study was to enhance the hydrothermal characteristics of asphalt mixtures through an investigation into the impact of anhydrous calcium sulfate whisker (ACSW) and polyester fiber (PF) on the hydrothermal properties of asphalt mixtures. In this paper, a central composite concatenation design (CCC) was employed to determine the optimal combination of ACSW and PF contents, as well as the asphalt aggregate ratio (AAR). Each influencing factor was assigned three levels for analysis. The evaluation indexes included dynamic stability, retained Marshall stability, and tensile strength ratio. Using the analysis methods of variance and gray correlation degree analysis, the hydrothermal properties of the asphalt mixture were examined in relation to the contents of ACSW, PF, and AAR based on the CCC results. Consequently, the optimal mix design parameters for composite modified asphalt mixture incorporating ACSW and PF were determined. The results indicated that the asphalt mixtures with hydrothermal qualities exhibited optimal performance in terms of 4.1% ARR, 11.84% ACSW, and 0.4% PF. The interaction between AAR and ACSW content had a greater effect on the dynamic stability and tensile strength ratio of the asphalt mixture, whereas the incorporation of ACSW and PF had a greater effect on the retained Marshall stability of the asphalt mixture. Among the three contributing factors, AAR exhibited the strongest relationship with the hydrothermal characteristics of the asphalt mixture, followed by the ACSW content; the correlation of PF content was the lowest. Therefore, to enhance the hydrothermal characteristics of the asphalt mixture, it is important to conduct a full evaluation of the constituents of ACSW and PF, along with the AAR in hot–humid regions.


Introduction
Asphalt pavements in different climatic regions have various performance requirements for asphalt mixtures [1][2][3][4].In order to comprehensively analyze the performance of asphalt mixtures, it is important to have a thorough comprehension of the meteorological attributes pertaining to the specific geographical region in which the road is situated.Asphalt pavements in hot-humid areas are often influenced by high temperatures, rain, and other natural conditions and are prone to rutting, and potholes in the process of use, which result in the overall functionality and durability of the asphalt pavements being significantly compromised.Therefore, researchers have occasion to conduct research on methods that would enhance the properties of asphalt pavements to be suitable for hot-humid areas.
Current research shows that there is a certain correlation between the roadway effectiveness of modified asphalt mixtures and the proportion parameters (externally modified amount and asphalt aggregate ratio).The results presented by Hong R.B. showed that when the polyester fiber (PF) content was 0.4% and the gangue powder replacement rate was 50%, the interaction between asphalt, PF, and gangue powder caused asphalt mixtures to exhibit excellent low-temperature anti-cracking properties [5].Zhao X.W. studied the effect of Sasobit and Deurex on the characteristics of crystalline tetra-alpha-ethylene-modifying asphalt mixes and found that 4% Sasobit +2% Deurex was the optimal amount for asphalt mixtures, and they had better mechanical properties and economic benefits [6].P. Ahmedzade found that when the asphalt content was 5.3%, asphalt mixtures made from 8% tall oil bitumen and 6% styrene-butadiene-styrene had the best physical and mechanical properties [7].Yang B. found that the best performance of large-particle-size-modified asphalt mixtures was achieved at elevated temperatures when the asphalt aggregate ratio was 3.23% and the rock asphalt content was 20% [8].Zhai R.X. showed that 5% nano-CaCO 3 and 4% SBR exhibited exceptional resistance to rutting deformation under elevated temperatures [9].In the study conducted by Gong Y.F., it was observed that the composite-modified asphalt mixtures retained their pore properties in optimal water stability even after undergoing freeze-thaw cycles with a 3.9% basalt fiber, 5.1% nano-TiO 2 /CaCO 3 content, and 5.67% asphalt aggregate ratio [10].Yan K.ZH.found that 15% waste tire rubber and 4% amorphous polyolefin was the optimal content for achieving the best performance for asphalt mixtures [11].Zhang H.T. concluded that asphalt mixtures modified using 4% SBS and 15% rubber were appropriate for asphalt pavements located in high-temperature, heavy traffic regions, and asphalt mixtures containing 4% SBS and 3% SBR were deemed suitable for cold-humid regions [12].Cheng Y.CH.found that diatomite and basalt-fiber-compositemodified asphalt mixtures were more resistant to seasonal permafrost than the base asphalt mixtures in areas with freeze-thaw cycle damage [13].Zhe H. demonstrated that a 2% thermoplastic polyurethane and 6% amorphous poly-alpha-olefin-modified asphalt mixtures had better rutting and moisture resistance than base asphalt mixtures [14].Therefore, optimizing the mix design of modified asphalt mixtures to adapt to the climatic conditions in hot-humid regions has a beneficial impact on the enhancement of asphalt pavements' resistance to high temperature and water damage.
A common method for creating a mix ratio design for externally mixed and modified asphalt is to first determine the optimum asphalt aggregate ratio of asphalt mixtures under different externally mixed mixtures using the Marshall test, and then to carry out mixture performance tests such as the wheel tracking test, low bending test, immersed Marshall test, and freeze-thaw cycle test.Asphalt mixture proportion parameters with better overall performance are then obtained according to the results of the various performance tests for asphalt mixtures under different external mixing amounts.Nevertheless, this method usually requires a large amount of experimentation, resulting in the squandering of both resources and time.Thus, researchers have conducted relevant research on convenient methods applicable to the optimization of the design parameters of asphalt mixture ratios, such as response surface design and multifactorial orthogonal design.Fan T.T. used the central composite circumscribed design (CCC) method, a type of response surface design, to optimize the design parameters of polyester fiber and anhydrous calcium sulfate whisker composite modified asphalt mixtures with better low-temperature properties [15].Ahmed I. Nassar utilized response surface design to optimize the mixing ratio of an emulsified asphalt mixture [16].Taher B.M. used response surface design to optimize the content of asphalt and polyethylene terephthalate in asphalt mixtures, and the results showed that the Marshall properties of the modified asphalt mixture achieved the best performance after adding 5.88% asphalt and 0.18% polyethylene terephthalate [17].Wei Y. revealed the relationship between the bone glue/polyurethane content and the properties of an asphalt binder using response surface design, and the optimal contents were 6.848% bone glue and 2.759% polyurethane [18].Cheng Y.CH conducted a comprehensive study on the optimization design and testing of a diatomite-and basalt-fiber-compound-modified asphalt mixture.Through the application of grey correlation grade analysis, the ideal mixing scheme was determined.The experimental findings indicated that the incorporation of 14% diatomite, 0.32% basalt fibers, and 5.45% asphalt aggregate resulted in the improved high-and low-temperature performance of the asphalt mixture [19].Du J.CH. studied an approach for conducting grey relational regression analysis on hot-mix asphalt [20].Lyu ZH.H. conducted an analysis of the performance of cold recycled mixtures.Additionally, an objective theory was optimized for the material content based on a multi-index-weighted grey methodology [21].Slebi Acevedo, C.J. used multi-standard decision-making technology to study the functional and mechanical properties of polyolefin-aramid fiber and hydrated-lime-modified porous asphalt mixtures [22,23].
Therefore, this article conducts a study on the optimization of mix proportion design and hydrothermal-performance-influencing factors of the asphalt mixture with the addition of ACSW and PF.Three mix ratio design parameters (namely, ACSW content, PF content, and asphalt aggregate ratio) and three hydrothermal performance evaluation indexes (namely, dynamic stability, retained Marshall stability, and tensile strength ratio) were selected.The central composite circumscribed design (CCC) and grey correlation grade analysis (GCGA) were used to perfect the mix ratio design characteristics of modified asphalt mixtures in hot-humid areas.The effects of the design parameters on the hydrothermal performance of the asphalt mixtures were explored; meanwhile, the correlation degree between the design parameters and the hydrothermal performance indexes was evaluated using GCGA.

Materials
SK-90# asphalt (penetration grade 80/100) was used, and its properties are shown in Table 1.The aggregate was generated in Xianyang City, Shaanxi Province.According to Chinese standard JTG E20-2011 [24], 20 sets of asphalt mixtures were prepared with an aggregate particle size characterized by a nominal maximum of 19 mm, as shown in Tables 2 and 3.The characteristics of ACSW and PF are listed in Tables 4 and 5, respectively.The rutting test is a commonly employed technique for assessing the durability of asphalt mixtures in terms of their ability to withstand permanent deformation when subjected to repetitive loading at elevated temperatures.According to Chinese standard T0719, the evaluation of the anti-deformation of asphalt mixtures under high temperatures was conducted using the rutting test.Specimens of 300 × 300 × 50 mm were made using the wheel-rolling method according to Chinese standard T0703.Before testing, the samples were placed in an oven at 60 • C for 5 h.Then, the rutting test was conducted under 0.7 MPa rubber wheel pressure, a 42 times/min loading rate, and a test temperature of 60 • C. The value of dynamic stability (DS) was calculated based on Equation (1).
where N is the loading speed, 42 times/min; t 1 and t 2 are the test times, usually 45 min and 60 min, respectively; d 1 and d 2 are the surface rut depth of the asphalt mixture specimens in correspondence with the test times t 1 and t 2 , mm; and c 1 and c 2 are the test parameters, and their values are 1.0, respectively.

Immersed Marshall Test
The immersed Marshall test is the main method employed for ascertaining the moisture sensitivity of asphalt mixtures under high-temperature immersion.The test was carried out according to Chinese standard T0709.Eight samples were generated for each mixture using the Marshall method, and these samples were evenly divided into two groups.The first group of specimens should be soaked in water at 60 • C for 30 min and tested for the stability value MS 0 .The second group of specimens should be soaked in water at 60 • C for 48 h and tested for the stability value MS 1 .The assessment of moisture stability in the asphalt mixtures was conducted by obtaining the retained Marshall stability (RMS), as shown in Equation (2).

Freeze-Thaw Splitting Test
The moisture stability of the titanium combination was assessed using a frozen splitting test, conducted in line with Chinese standard T0729.A total of eight samples were generated using the Marshall approach and subsequently separated into two groups, resulting in an equal distribution.A single group of samples was subjected to ambient temperatures ranging from around 20 to 25 • C.An additional set of samples underwent a freezing process for 16 h at −18 • C after water saturation.These samples were subsequently immersed in a solution at 60 • C for 24 h.Following this, the two sets of specimens were immersed in water at a temperature of 25 • C for 2 h.Subsequently, an indirect tensile test was conducted, employing an indirect tensile load with a loading rate of 50 mm/min.The tensile strength ratio (TSR) is determined by calculating the average split stretch strength before (S before ) and after the freeze-thaw cycle (S after ), as demonstrated in Equation ( 3).

TSR =
S a f ter S be f ore × 100%

Grey Correlation Grade Analysis
To ascertain the most suitable composition of ACSW and PF in the asphalt mixture, the CCC results were examined using GCGA [25].The grey correlation coefficient pertaining to the reference sequence g 0 = (g 0 (1), g 0 (2), g 0 (3)) and the comparison sequence g i = (g i (1), g i (2), g i (3)), ξ i (m) are defined as Equation ( 4).
where ξ i (m) is the correlation coefficient between the m-th evaluation index of the i-th comparison group and the associated reference group.max Then, the grey correlation grade λ i is derived by calculating the weighted sum of the grey correlation coefficients, as outlined in Equation (5).
where Z D , Z R , and Z T represent the weighting of DS, RMS, and TSR, respectively.

Experimental Design
CCC, a type of response surface design, is used to assess the interaction between independent and response variables within the context of an experiment.In the CCC technique, at least two numerical inputs are required and are varied in the range of alpha (α) to five (−α, −1, 0, +1, +α) (α > 1) phases.The CCC method offers the benefit of incorporating both sequential and rotational aspects within the framework of experimental design.In this paper, CCC was chosen to evaluate the association between the independent and response variables with the aim of identifying the most suitable experimental formulation.The number of test times can be expressed as (2 m + 2m + n), where m represents the count of independent variables and n is the number of center points.The value of α can be determined by 2 m/4 .
The experimental design employed in this study necessitates the execution of 20 experimental trials.It encompassed eight factorial points, six star points, and six center points, all of which were conducted at five levels.The design parameters were as follows: m = 3, n = 6, and α = 1.682.The decision to establish six center points was made in order to account for potential experimental errors by the inclusion of additional replications.Three independent variables were the asphalt aggregate ratio (AAR), ACSW content (related to asphalt-binding mass ratio), and PF content (relative to asphalt mix mass ratio), as presented in Table 6.Three response variables, DS, RMS, and TSR, were chosen to optimize the content of ACSW and PF, and the AAR of the asphalt mixture.Table 7 presents the experimental design and corresponding outcomes.Design-Expert 10.0 software was used to perform a response model (as shown in Equation ( 7)) and statistical regression analysis of the experimental program to determine the best combination of independent variable levels.
where the response variable is denoted by R. The coded independent parameters are represented by F i and F j .The coefficients for the constant, primary, interaction, and quadratic terms are denoted by µ 0 , µ i , µ ij , and µ ii , respectively.The random error is denoted by η.

Statistical Modeling
To study the interaction between different independent variables and their impact on the responsive variable, analysis of variance (ANOVA) was used.When the p-value was less than 0.05, there was a significant difference between the model and the factor variable.The results of the ANOVA are presented in Tables 8 and 9, respectively.The ANOVA results for the quadratic model of the DS are presented in the first row of Table 8.The R-squared value of 0.9389 and the adjusted R-squared value of 0.8839 indicated that the fit of the DS model was strong.Furthermore, Adeq.precision was adopted to evaluate the signal-to-noise ratio of the response model.It is worth noting that the regression model was considered optimal when the Adeq.precisionvalue exceeded 4. Based on the findings shown in Table 8, it can be observed that the Adeq.precision of the DS model was 13.292, suggesting the favorable match of the model.The statistical significance of the DS model is supported by the p-value, which is less than 0.0001.The ANOVA results in Table 8 demonstrate that the models for RMS and TSR exhibited a good fit.
The significance of Lack of Fit is used to evaluate the model fit.The p-values of the Lack of Fit for the DS, RMS, and TSR models, as presented in Table 9, were 0.4540, 0.0552, and 0.2001, respectively, which suggests that all models were well fitted.

Optimization Results Based on CCC
For the hot-humid area, the target values listed in Table 10 were chosen based on Chinese standard JTG F40-2004 [26].The larger the values of response variables, the better the anti-deformation and anti-water damage of the asphalt mixture.The most suitable values of the layout parameters were determined by calculating Equations ( 8)- (10).The ramps of CCC analysis are depicted in Figure 1.Each ramp is associated with a specific point that represents the intended aim for both the independent and dependent variables.The meaning of the red dot in Figure 1 is the value of the ideal design parameters of asphalt mixture, and the blue dot is the value of the predicted values of hydrothermal performance of asphalt mixture.Based on the data presented in Figure 1, it can be observed that the ideal design parameters for the ACSW and PF compound-modified asphalt mixture were 4.1% AAR, 11.84% ACC, and 0.42% PFC.Table 11 shows the optimal predicted values and laboratory values of response variables of the asphalt mixtures.The laboratory tests were conducted three times.The assessment of the deviation rate between the predicted and laboratory results was conducted using Equation (11).Table 11 shows the optimal predicted values and laboratory values of response variables of the asphalt mixtures.The laboratory tests were conducted three times.The assess-ment of the deviation rate between the predicted and laboratory results was conducted using Equation (11).In Table 11, the deviation rates of all three response variables were lower than 2%, which showed that there was no discernible distinction between the predicted values and laboratory values.This demonstrated that it was feasible to enhance the hydrothermal characteristics of the modified asphalt mixture by optimizing its mix proportion using the CCC-RSM approach.

Optimization Results Based on GCGA
The GCGA is a scientific methodology utilized to evaluate the correlation between components by analyzing the similarity in their developmental tendencies.This approach holds significant potential as a viable answer for research programs focused on multiobjective optimal design.The challenges in drawing accurate findings when making comparisons arise from the varying physical interpretations of the measured indicators and the lack of standardization in data dimensions.
Typically, to obtain better hydrothermal properties of the asphalt mixture, the index values of high-temperature performance and water stability are needed to reach the maximum.In GCGA, the reference series and the comparative series need to be established first.The maximum values of DS, RMS, and TSR are defined as reference series; the values of DS, RMS, and TSR are set to 2500 times/mm, 100%, and 100%, respectively.The results of 20 schemes are defined as comparative series.Secondly, the numerical dimensionless test is carried out using Equation (12), and the dimensionless results are shown in Table 12.
where x n represents the dimensionless result of the n-th group, n = 1, 2, . . ., 20; and f n represents the test result of group n. ξ i (m) of the dimensionless test results is computed based on Equation ( 4), as presented in Table 12.In the proportion design of the asphalt mixture in hot and humid regions, the high-temperature performance index and water stability index of the asphalt mixture are equally important.Consequently, the performance indicators are assigned equal weights in order to derive the most ideal contents of ACSW and PF, as well as the optimal asphalt aggregate ratio.The weights assigned to the performance indicators are determined as follows: Z D = 0.5, Z R = 0.25, and Z T = 0.25.The λ is calculated using Equation ( 5), as presented in Table 12.Table 13 displays the average correlation grade associated with each evaluation and level.In the GCGA algorithm, as per its specifications, the greater the correlation rank, the closer the corresponding factor level is to the optimal value.As can be seen from Table 13, the factors correspond to the optimal ratios of ACC = 11%, PFC = 0.4%, and AAR = 4.0%.Hence, in accordance with the multi-objective results of CCC, modified asphalt mixture performs better at high temperature and has greater moisture susceptibility when the ACSW content is 11%, the PF content is 0.4%, and the asphalt aggregate ratio is 4.0%.
In addition, compared with the results of CCC and GCGA, as shown in Table 14, the deviation rates of the optimal design parameters of the modified asphalt mixtures were smaller than 7%, which revealed that it was feasible to carry out the multi-objective optimized design of asphalt mixtures using the GCGA method.Dynamic stability (DS) was adopted to evaluate the high-temperature performance of the asphalt mixture, and the larger the DS value, the better its high-temperature perfor-mance.The DS results of the ACSW-and PF-compound-modified asphalt mixture (APCRA) at different AAR, ACC, and PFC are shown in Figure 2. Figure 2 shows that the coupling effects of AAR and ACC, as well as AAR and PFC, exhibited more significant effects on the DS values of APCRA than the interaction between ACC and PFC, which is in accordance with the findings in Table 9.In Figure 2a,c, the DS Figure 2 shows that the coupling effects of AAR and ACC, as well as AAR and PFC, exhibited more significant effects on the DS values of APCRA than the interaction between ACC and PFC, which is in accordance with the findings in Table 9.In Figure 2a,c, the DS values of APCRA first increased as the ACSW content increased from a lower level (9.0%) to a higher level (13.0%), and then stabilized with a further increase in the AAR, reaching a maximum value when the AAR was 4.1%.When the AAR was constant, the DS value of APCRA was less affected by the ACSW content.In Figure 2b, the DS value of APCRA first increased at constant PFC and then increased with AAR, and this pattern of change became more pronounced when the PFC was at a low level (0.3%) compared to a high level (0.5%).

Moisture Susceptibility
The evaluation of moisture sensitivity in APCRA was conducted using two parameters: residual Marshall stability (RMS) and tensile strength ratio (TSR).The larger the RMS or TSR value, the better the moisture damage resistance of APCRA.

RMS
Figure 3 shows the RMS results of APCRA at different AAR, ACC, and PFC.The coupling effect of ACC and PFC on the RMS values of APCRA was more significant than the effect of AAR and ACC or AAR and PFC on the RMS values of APC.As shown in Figure 3a,c, when the AAR or PFC was kept constant, the RMS values of APCRA increased with the increase in ACC.When the value of ACC was kept constant, the RMS values of APCRA increased and then decreased with the increase in PFC.In addition, an increase in AAR did not result in a significant change in the RMS values of APCRA.
the effect of AAR and ACC or AAR and PFC on the RMS values of APC.As shown in Figure 3a,c, when the AAR or PFC was kept constant, the RMS values of APCRA increased with the increase in ACC.When the value of ACC was kept constant, the RMS values of APCRA increased and then decreased with the increase in PFC.In addition, an increase in AAR did not result in a significant change in the RMS values of APCRA.
From Figure 3a,b, the RMS values of APCRA did not change significantly with the increase in AAR at the same PFC.Meanwhile, at the same AAR, the RMS values of APCRA increased and then decreased with the increase in PFC.The RMS values of APRA reached a relative maximum when the PFC was at a medium level (0.4%) and the AAR was at a high level (4.5%).From Figure 3a,b, the RMS values of APCRA did not change significantly with the increase in AAR at the same PFC.Meanwhile, at the same AAR, the RMS values of APCRA increased and then decreased with the increase in PFC.The RMS values of APRA reached a relative maximum when the PFC was at a medium level (0.4%) and the AAR was at a high level (4.5%).

TSR
The TSR results of APCRA at different AAR, ACC, and PFC content levels are shown in Figure 4. Figure 4a shows that the coupling relationship between ACC and AAR had a significant effect on the TSR values of APCRA.When the ACC value was at a low level (9.88%), the TSR values gradually increased with the increase in AAR, and when the ACC value was at a high level (13.88%), the TSR values gradually decreased with the increase in AAR.In addition, when the coupling values of ACC and AAR changed from 9.88% (ACC) and 3.7% (AAR) to 13.88% (ACC) and 4.7% (AAR), the TSR values increased first and then decreased.
Figure 4b shows that the TSR values of APCRA gradually increased when the coupling values of AAR and PFC changed from 0.3% (PFC) and 3.7% (AAR) to 0.5% (PFC) and 4.7% (AAR).Figure 4c illustrates that when the value of PFC was lower than 0.4%, the TSR values of APCRA increased with the increase in ACC, whereas when the value of PFC was higher than 0.4%, the TSR values of APC gradually decreased with the increase in ACC. and then decreased.
Figure 4b shows that the TSR values of APCRA gradually increased when the coupling values of AAR and PFC changed from 0.3% (PFC) and 3.7% (AAR) to 0.5% (PFC) and 4.7% (AAR).Figure 4c illustrates that when the value of PFC was lower than 0.4%, the TSR values of APCRA increased with the increase in ACC, whereas when the value of PFC was higher than 0.4%, the TSR values of APC gradually decreased with the increase in ACC.

Effect Degree of Factors on Performance Indexes
In this paper, GCGA was used to evaluate the correlation rank of the influencing factors with the performance indexes.The AAR, ACC, and PFC values were selected as comparison sequences, while the DS, RMS, and TSR results were adopted as reference sequences.Using Equations ( 4) and ( 12), the dimensionless results and ξi(m) were calculated as shown in Tables 15 and 16.

Effect Degree of Factors on Performance Indexes
In this paper, GCGA was used to evaluate the correlation rank of the influencing factors with the performance indexes.The AAR, ACC, and PFC values were selected as comparison sequences, while the DS, RMS, and TSR results were adopted as reference sequences.Using Equations ( 4) and ( 12), the dimensionless results and ξ i (m) were calculated as shown in Tables 15 and 16.
From Table 16, it can be observed that the grey correlation degrees of the different factors on various performance indicators were different.The grey correlation degrees of DS were as follows: AAR (0.660) > PFC (0.614) > ACC (0.593); on RMS: AAR (0.788) > ACC (0.748) > PFC (0.696); and on TSR: AAR (0.810) > ACC (0.710) > PFC (0.683).This indicates that the effect of AAR on the high-temperature and moisture sensitivity of the asphalt mixtures was more significant than that of ACC and PFC.In addition, ACC mainly affected the moisture sensitivity of the asphalt mixture, while PFC mainly affected the high-temperature performance of the asphalt mixture.anhydrous calcium sulfate whisker and polyester fiber compound modified asphalt mixtures in hot-humid areas.(2) The proportions of 11.8% anhydrous calcium sulfate whisker, 0.42% polyester fiber, and 4.1% asphalt aggregate can produce an asphalt mixture with better hydrothermal performance.
(3) The asphalt aggregate ratio has a greater influence on the hydrothermal performance, polyester fibers mainly affect the moisture susceptibility, and anhydrous calcium sulfate whiskers mainly affect the high-temperature properties of the asphalt mixture.
Recommendations: Feasibility research on the application of calcium sulfate whiskers as a renewable solid-waste resource in pavement materials has just begun.This paper studied the influence of calcium sulfate whiskers on the high-temperature deformation resistance and water damage resistance of an asphalt mixture and puts forward an optimization method for the mix ratio design of a modified asphalt mixture based on hygrothermal performance.However, the influence and mechanism of calcium sulfate whiskers on the low-temperature crack resistance and durability of asphalt mixtures, as well as the optimization method of the mix ratio design of modified asphalt mixtures considering physical and mechanical properties and road performance, need to be further studied.

i max m |x 0
(m) − x i (m)| is the maximum difference between the reference and comparison series.max i max m |x 0 (m) − x i (m)| is the minimum difference value between the reference sequence and comparative sequences.

1 .Figure 2 .
Figure 2. Influence of three factors on DS: (a) AAR and ACC; (b) ARR and PFC; and (c) ACC and PFC.

Figure 2 .
Figure 2. Influence of three factors on DS: (a) AAR and ACC; (b) ARR and PFC; and (c) ACC and PFC.

Figure 3 .
Figure 3. Influence of three factors on RMS: (a) AAR and ACC; (b) ARR and PFC; and (c) ACC and PFC.

Figure 3 .
Figure 3. Influence of three factors on RMS: (a) AAR and ACC; (b) ARR and PFC; and (c) ACC and PFC.

Figure 4 .
Figure 4. Influence of three factors on TSR: (a) AAR and ACC; (b) ARR and PFC; and (c) ACC and PFC.

Figure 4 .
Figure 4. Influence of three factors on TSR: (a) AAR and ACC; (b) ARR and PFC; and (c) ACC and PFC.

Table 2 .
Passing rate and physical characteristics of aggregates.

Table 3 .
Physical characteristics of mineral powder.

Table 5 .
Properties of PF.

Table 6 .
Experimental factors and levels of CCC.

Table 7 .
Results of CCC.

Table 8 .
ANOVA results of CCC.

Table 9 .
ANOVA results for experimental variables.

Table 10 .
Target values of response variables.

Table 11 .
Validation laboratory results of CCC.

Table 11 .
Validation laboratory results of CCC.

Table 12 .
Analysis results of grey correlation degree.

Table 13 .
Average relational between levels of each factor and response.

Table 14 .
Deviation rate of CCC and GCGA results.