Optimization Design of the Mix Ratio of a Nano-TiO2/CaCO3-Basalt Fiber Composite Modified Asphalt Mixture Based on Response Surface Methodology

This research optimizes the mix ratio of nano-TiO2/CaCO3 (NTC)-basalt fiber (BF) composite modified asphalt mixture. Based on the Box–Behnken method and the response surface method, a three-factor and three-level test was designed. The input indicators were the asphalt–aggregate ratio, NTC content, and BF content. The output indicators were the density, air voids, Marshall stability, flow value, voids in mineral aggregate (VMA), and voids filled with asphalt (VFA) values of the asphalt mixture. The response surface model was established according to the test response index value. Then, the function was fitted through multiple regression equations and multivariate analysis of variance was performed. Finally, according to the specification requirements and engineering needs, the selected conditions of each response value were determined to optimize the asphalt–aggregate ratio and the contents of NTC and BF, and the predicted values were verified through the measured data. The test results show that the optimal contents of NTC and BF and the optimal asphalt–aggregate ratio were 5.1%, 3.9%, and 5.67%, respectively. By comparing the measured Marshall test index value with the predicted value, the minimum relative error was 0.096% and the maximum error was 6.960%. The results show that response surface methodology can be used to optimize the mix ratio of composite modified asphalt mixtures.


Introduction
In recent years, with the increase of traffic flow and the deterioration of the environment, it is becoming increasingly difficult for existing asphalt pavement to meet requirements, meaning modified asphalt technology is very important. Improving the durability and stability of pavement through modification of asphalt is a new challenge [1]. With the emergence of nano-modified and fiber-modified asphalts, the performance of asphalt pavement has been effectively improved, and the service life of asphalt pavement has also been improved.
Due to the large specific surface area of the nanomaterials, after being added to the asphalt, the nanomaterials can improve the cohesion and also increase the viscosity and high temperature stability of the asphalt, thereby improving the high temperature performance and water stability of the asphalt mixture. Chang et al. [2] added different proportions of nanometer calcium carbonate to SBS (styrene-butatriene-styrene) modified asphalt and analyzed the high-and low-temperature stability of the composite modified asphalt via changes of the nanometer calcium carbonate content. Manias et al. [3] used silicate nanomaterials and polymers to prepare a composite modified asphalt. low-temperature crack resistance of BF-modified asphalt mixture was significantly improved. Based on the three-dimensional model and numerical analysis, Zhang et al. [26,27] studied the rheological and viscoelastic properties of basalt-fiber-reinforced asphalt mastic and mortar. The numerical analysis showed that the BF caused the stress redistribution of the asphalt mortar and reduced the stress of the asphalt mortar.
Compared with conventional regression model approach, the response ratio optimized by the response surface method is more accurate and is becoming increasingly popular in the construction industry [28]. The response surface method is a collection of statistical methods used for designing experiments, developing models, assessing the effects of the experimental factors, and optimizing the process [29,30]. Wang et al. [31] used the response surface methodology to study the effects of basalt fiber content, length, asphalt-aggregate ratio, and other factors on the volume and strength properties of SBS-modified asphalt mixture. Chavez et al. [32] applied the response surface method to evaluate the aging responses of asphalt binder. Moghaddam et al. [33] introduced the response surface method into the mix design for an asphalt mixture. According to the laboratory test results, it was finally determined that the optimal mix ratio was 5.88% whetstone and 0.18% PET (Polyethylene terephthalate) content. Based on the extensive literature, the response surface method can be successfully applied to study asphalt binders and mixtures.
In this study, an experimental scheme for a nano-TiO 2 /CaCO 3 (NTC)-BF composite modified asphalt mixture was designed based on the response surface methodology. The NTC content, BF content, and asphalt-aggregate ratio were independent variables, while the density, Marshall stability, flow value, air voids, voids filled with asphalt (VFA), and voids in mineral aggregate (VMA) were responses or dependent variables. The relationships between preparation parameters and Marshall test indices were analyzed to evaluate the effects of the NTC content, BF content, and asphalt-aggregate ratio. A design optimization of the NTC-BF composite modified asphalt mixture was proposed and validated with experimental results.

Asphalt
The matrix asphalt used in this paper was AH #90 from Panjin of Liaoning Province, China. The penetration, softening point, ductility, Brinell viscosity, and density of the asphalt were measured according to the "Testing Regulations for Asphalt and Asphalt Mixtures of Highway Engineering" (JTG E20-2011) [34]. The specific test results are listed in Table 1. In Table 1, the Brinell viscosity is the apparent viscosity of the asphalt measured with a Brookfield viscometer according to the specification JTG E20-2011 [34]. The BF selected in the experiment comprised chopped fibers produced by Jilin Jiuxin Basalt Industry Co., Ltd. The test results of basalt fiber index parameters and the limits in the "Classification, Gradation, and Designation of Basalt Fiber" (GB/T 38111-2019) regulations [35] are shown in Table 2. The basic technical indicators are listed in Table 2. In addition to the basic technical indicators listed in Table 2, the BF also has low hydrophilicity, high heat resistance, and high oil absorption (oil absorption rate of 6.154%), meaning the BF can be applied in composite modified asphalt.

Nanomaterials
The selected NTC material (TC-816) was produced by Tiancheng High-Tech Nano-Composite Co., Ltd. It is a white material made of titanium dioxide (TiO 2 ) and inorganic synthetic materials via chemical methods. Its molecular formula is TiO 2 /CaCO 3 . The ratio of TiO 2 /CaCO 3 is approximately 1:3.5. Its appearance is shown in Figure 1. The selected NTC material is chemically stable and hardly reacts with other elements or compounds at room temperature, and is insoluble in water. The support particles are rod-shaped, and the outer titanium dioxide particles are in the form of sheets (diameter 50-60 nm), which have a regular structure and excellent hiding power. According to the determination, the 45 micron sieve residue of the NTC material is 0.02% and the stacking density is 0.4 g/ml.

Nanomaterials
The selected NTC material (TC-816) was produced by Tiancheng High-Tech Nano-Composite Co., Ltd. It is a white material made of titanium dioxide (TiO2) and inorganic synthetic materials via chemical methods. Its molecular formula is TiO2/CaCO3. The ratio of TiO2/CaCO3 is approximately 1:3.5. Its appearance is shown in Figure 1. The selected NTC material is chemically stable and hardly reacts with other elements or compounds at room temperature, and is insoluble in water. The support particles are rod-shaped, and the outer titanium dioxide particles are in the form of sheets (diameter 50-60 nm), which have a regular structure and excellent hiding power. According to the determination, the 45 micron sieve residue of the NTC material is 0.02% and the stacking density is 0.4 g/ml.

Aggregates and Fillers
The coarse and fine aggregates selected in this paper were all from the Shuangyang quarry in Changchun City, Jilin Province, and the ore powder came from the limestone ore powder produced in Shuangyang Shixi, Changchun City, Jilin Province. The coarse and fine aggregates used in the test were screened step-by-step, which ensured the accuracy of the test and reduced the impact of the discreteness of the aggregates. According to the relevant requirements of the "Highway Engineering Aggregate Test Regulations" (JTG E42-2005) [36], the physical properties of coarse aggregates, fine aggregates, and mineral powder were tested. The results are shown in Tables 3-5.

Aggregates and Fillers
The coarse and fine aggregates selected in this paper were all from the Shuangyang quarry in Changchun City, Jilin Province, and the ore powder came from the limestone ore powder produced in Shuangyang Shixi, Changchun City, Jilin Province. The coarse and fine aggregates used in the test were screened step-by-step, which ensured the accuracy of the test and reduced the impact of the discreteness of the aggregates. According to the relevant requirements of the "Highway Engineering Aggregate Test Regulations" (JTG E42-2005) [36], the physical properties of coarse aggregates, fine aggregates, and mineral powder were tested. The results are shown in Tables 3-5.

Response Surface Methodology
The response surface methodology was originally used to optimize biological processes. When there are multiple influencing factors, the relationship between the influencing factors can be obtained by regression design method, so as to optimize the test process. Under the premise of considering random errors, this method selects fewer test points and reduces the amount of tests. After the test results are obtained, the results are fitted and a continuous variable regression surface model is established to finally obtain the optimal range of influence factors. Compared with the orthogonal test, the obtained response surface is a continuous surface model with greater visibility, and the model can be used to predict points outside the test range. This article only briefly introduces the response surface method. For details, see the literature [37].
The response surface methodology is constructed on the basis of the least-squares method. The optimization process is as follows. In the first step, assume the relationship between a corresponding value and multiple impact factors: In the formula, y is the response value; x 1 , x 2 , . . . . . . , x k are the influence factors; ε is the random variable.
If there is an approximately linear relationship between the response value and the impact factor, then make a first-order Taylor expansion of Formula (1): In the Equation, β i is the slope or linear effect of the influence factor x i . Fit the influence factor with the above Equation.
If the fitting is successful, an appropriate influence factor can be selected to obtain the best response value. If it is found that the relationship between y and xi cannot be expressed by the linear Taylor expansion, then the second-order transformation of Equation (2) needs to be further performed to obtain the following second-order model Equation: In the Equation, β i is the linear effect of the independent variable x i ; β ij is the interaction between the independent variables; B ii is the secondary effect of x i . This can be expressed as a matrix: In the Equation, X = (x 1 , x 2 , · · · , x k ) T ; b = (β 1 , β 2 , · · · , β k ) T is the regression coefficient matrix; B is the k-th order symmetric matrix.
If the fitting is successful and the response reaches the best, then the following Equation is satisfied: The solution X in Equation (5) is called the stable point, and the transformation can be obtained as: Substituting Equation (6) into Equation (4), the predicted response value at the stable point is obtained as: There are many types of design models for response surface analysis, such as the central composite, Box-Behnken design, D-optimal design, and uniform design models. The specific design schemes of the various models are different, meaning the numbers of tests required are also different. It is necessary to select the appropriate model for the design according to the actual needs.

Test Design
The AC-13 gradation of asphalt mixtures is a fine-grained, continuous gradation, and is mostly used for pavement surface layer. Therefore, the asphalt mixture in this paper adopted the AC-13 gradation. The gradation of asphalt mixtures, including the upper limit, lower limit, median, and synthetic gradation curve values of the aggregates, are shown in Figure 2. The synthetic gradation was used in this study to mix the aggregates of multiple grades in proportion and close to the median value of the gradation. The gradation curve of the asphalt mixture is the same as in reference [38]. For the experimental design, this paper adopted the Box-Behnken method. NTC and BF are asphalt admixtures. The upper and lower limits of admixtures and the ranges of asphalt ratios were determined according to the literature [10,[39][40][41]. The selected three impact factor levels are shown in Table 6.
was used in this study to mix the aggregates of multiple grades in proportion and close to the median value of the gradation. The gradation curve of the asphalt mixture is the same as in reference [38]. For the experimental design, this paper adopted the Box-Behnken method. NTC and BF are asphalt admixtures. The upper and lower limits of admixtures and the ranges of asphalt ratios were determined according to the literature [10,[39][40][41]. The selected three impact factor levels are shown in Table 6.    In this study, the asphalt-aggregate ratio is measured by weight.
The Box-Behnken test design has an impact factor range of 3-7, which can better evaluate the non-linear correlation between the response index and the impact factor. Under the same conditions, the number of tests required is less than other methods. In the Box-Behnken test design, the number of tests required for the three impact factors was 17. There were six response output indicators that needed to be determined in the experiment, namely the density, Marshall stability, flow value, air voids, VMA, and VFA of the asphalt mixture. The impact factors of the 17 groups of experiments are set out in Table 7 below.

Production of Asphalt Mixture Test Piece
The asphalt mixtures were prepared according to JTG E20-2011 [34]. When preparing modified asphalt, the control asphalt and NTC were mixed using blending equipment at a speed of 4000 r/min for 45 min. The temperature for the mixing process was 150 • C. When preparing the asphalt mixture, the aggregates and the BF were mixed at a temperature of 160 • C. Then, the modified asphalt was poured in and mixed until the aggregates were coated. Next, the weighted fillers were added and mixed well. The Marshall specimens in this study were prepared at a height of 63.5 ± 1.3 mm and diameter of 101.6 ± 0.2 mm. The double-sided compaction occurred 75 times.

Measurement and Calculation of the Volume Index
According to the requirements of JTG E20-2011 [34], the dry surface method was used to determine the relative gross weight density, gross bulk density, air voids, mineral air voids, and VFA of the asphalt mixture. The specific calculation Equation was as follows: In the Equation, γ f is the relative density of the wool volume of the specimen, which is dimensionless; m a is the air weight of the dried specimen; m f is the surface dry weight of the specimen; m w is the weight of the test piece in water (g); ρ f is the test piece of the gross volume density; ρ w is the density of the water at 25 • C, the value of which is 0.9971 g/cm 3 .
In the Equation, γ f is the maximum theoretical density of the asphalt mixture, which is dimensionless; γ t is the relative density of the asphalt mixture gross volume, which is dimensionless.
In the Equation, VMA is the aggregate gap ratio of asphalt mixture; γ sb is the relative density of the synthetic wool volume of the mixture, which is dimensionless; p s is the sum of the percentage of various mineral materials in the total mass of the asphalt mixture.
In the Equation, VFA represents the voids filled with asphalt in the mixture.

Determination of Marshall Stability Test Index
According to the requirements of JTG E20-2011 [34], after the test specimen was kept in a constant temperature water bath at 60 ± 1 • C for 30-40 minutes, the Marshall tester was used to test the Marshall stability and flow value of the test piece.

Results and Discussion
The test results are shown in Table 8. After obtaining the statistics for the test data, the maximum value, minimum value, average value, and standard deviation of each output index could be obtained, which could be used as the basis for further analysis. The results are shown in Table 9. For the convenience of calculation, test groups with the same impact factor level are calculated using the same test data.

Density
By analyzing the linear model (linear), two-factor model (2FI), and quadratic equation model (quadratic), the best fitting equation model was found and the coefficient of determination R 2 was used to evaluate the fit.
The impact factor was fitted to the three model equations, and SST (sum of squares total) of each model was calculated. Then, the SST was decomposed: Therefore, the total dispersion was finally decomposed into the following Equation: In the Equation, SSR is the sum of squares within the group (regression square sum); SSE is the sum of squares between groups (sum of squared residuals).
The solvability coefficient R 2 is as follows: The closer SSR is to SST, the smaller the residual is and the closer R 2 is to 1, indicating a higher degree of fit. In order to exclude the influence of too many variables on the degree of fit, R 2 was adjusted by the following Equation to obtain the predicted R 2 : In the Equation, n − k − 1 is the degree of freedom of the sum of squares between groups; n − 1 is the degree of freedom of the sum of squares of the total dispersion. Table 10 lists the comparison results after fitting with various equations. For the test model: If the parameter is significantly non-zero, the following null hypothesis and alternative hypothesis are made: It can be found that when R 2 approaches 0, the F value also approaches 0; when R 2 approaches 1, the F value approaches infinity. Therefore, the test hypothesis H 0 is equivalent to the test R 2 = 0. Tables 11-13 list the results for the analysis of variance for each fitted model.
According to the analysis of the above tables, it is found that the quadratic equation had the highest degree of fitting, which was recommended. The variance of each variable fitted by the quadratic equation is analyzed below, and the results are shown in Table 14.
In the Equation, A is the amount of BF; B is the amount of NTC; C is the asphalt-aggregate ratio. The distribution of predicted and measured values is shown in Figure 3.
According to the results of analysis of variance, the fitting equation of the response index can be obtained: Density=1.78911+5.94747E−004×A+0.020×B+0.205×C−8.289E−005×A×B+0.0029×A×C+2.53744E−00 4×B×C−0.0024862×A 2 −1.62926E−003×B 2 −0.019378×C 2 In the Equation, A is the amount of BF; B is the amount of NTC; C is the asphalt-aggregate ratio. The distribution of predicted and measured values is shown in Figure 3. Figure 3 shows the relationship between the response index predicted by the final equation and the measured values. The smaller the distance between the point and the diagonal, the higher the accuracy of the prediction; the greater the degree of deviation from the top to the bottom, the greater the error. It can be seen that there are seven sets of measured values that coincide with the predicted values, and there are six sets of measured values that are distributed on the diagonal line but not far away. There are no outliers in the whole distribution, which indicates that the method of obtaining the predicted value via the fitting equation is reliable. Three-dimensional coordinates can only represent the effects of two impact factors at the same time, so when analyzing the interaction of the two factors, the third impact factor is taken as the  Three-dimensional coordinates can only represent the effects of two impact factors at the same time, so when analyzing the interaction of the two factors, the third impact factor is taken as the median value within the range and is fixed. The contour plot and response surface plot between the density and each influencing factor are shown in Figure 4.
It can be seen from Figure 4 that the location with the highest density is in the middle dark parts of the contour plot and the response surface plot. There are five measured points in each surface graph, indicating how well the predicted surface fits the actual values. It can be found from the surface graph that the predicted values fit the surface better. Increasing the amount of NTC and BF makes the mixture density increase first and then decrease, however the two factors influence each other. Under the premise that the asphalt-aggregate ratio is 5.5%, the effect of the fiber content on density gradually increases during the increase of the NTC content from 1% to 7%. In the process of increasing the fiber content from 3.5% to 9%, the effect of NTC on the density decreases slightly. In addition, as the asphalt-aggregate ratio increases, the density of the mixture increases first and then decreases.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 12 of 21 median value within the range and is fixed. The contour plot and response surface plot between the density and each influencing factor are shown in Figure 4. It can be seen from Figure 4 that the location with the highest density is in the middle dark parts of the contour plot and the response surface plot. There are five measured points in each surface

Air Voids
The comparison of the fitting degrees of various fitting equations in the air voids analysis is consistent with the method selected in the density analysis process in the previous section, and will not be repeated in this section. Table 15 is the variance test situation of each variable after fitting through the quadratic equation.
The distribution of the predicted and measured values is shown in Figure 5. It can be seen from the above figure that the fit between the predicted value and the measured value is good, indicating that the fitting equation was selected accurately. There are 5 test groups with air void contents greater than 6% in the figure, exceeding the "Technical Specifications for Highway Asphalt Pavement Construction" (JTG F40-2004) regulations [42]; therefore, it is necessary to limit the scope when optimizing the impact factor.
Using the same method as Section 3.1.1, the contour plot and the response surface plot between the air voids and each influencing factor were drawn. The results are shown in Figure 6. It can be seen from the above figure that the fit between the predicted value and the measured value is good, indicating that the fitting equation was selected accurately. There are 5 test groups with air void contents greater than 6% in the figure, exceeding the "Technical Specifications for Highway Asphalt Pavement Construction" (JTG F40-2004) regulations [42]; therefore, it is necessary to limit the scope when optimizing the impact factor.
Using the same method as Section 3.1.1, the contour plot and the response surface plot between the air voids and each influencing factor were drawn. The results are shown in Figure 6. It can be seen from the above figure that the fit between the predicted value and the measured value is good, indicating that the fitting equation was selected accurately. There are 5 test groups with air void contents greater than 6% in the figure, exceeding the "Technical Specifications for Highway Asphalt Pavement Construction" (JTG F40-2004) regulations [42]; therefore, it is necessary to limit the scope when optimizing the impact factor.
Using the same method as Section 3.1.1, the contour plot and the response surface plot between the air voids and each influencing factor were drawn. The results are shown in Figure 6. According to Figure 6, the addition of BF leads to an increase of the air voids in the asphalt mixture, while the addition of NTC reduces the air voids in the asphalt mixture; the interaction between NTC and the BF modifier is obvious. On the premise that the asphalt-aggregate ratio is 5.5%, when the NTC content is 1.0% the fiber content increases from 1.0% to 7.0%, while the air voids only increase by about 1.0%. When the NTC content is 9.0%, the fiber content increases from 1.0% to 7.0% and the air voids increase by about 3.0%. When the fiber content is 1.0%, during the increase of the NTC content the air voids decreases by about 4.0% and the effect is more obvious. When the fiber content is 7.0%, the increase of NTC reduces the air voids by about 1.0% and the effect is greatly reduced. The increase of the asphalt-aggregate ratio makes the air voids of the mixture decrease first and then increase. When the asphalt-aggregate ratio is relatively high or low the influences of fiber and NTC are reduced, and when the asphalt-aggregate ratio is moderate the effects of the two are amplified.

Marshall Stability
The quadratic equation was selected as the design model for the fitting calculation. Then, the variance test was performed. The results are shown in Table 16.  According to Figure 6, the addition of BF leads to an increase of the air voids in the asphalt mixture, while the addition of NTC reduces the air voids in the asphalt mixture; the interaction between NTC and the BF modifier is obvious. On the premise that the asphalt-aggregate ratio is 5.5%, when the NTC content is 1.0% the fiber content increases from 1.0% to 7.0%, while the air voids only increase by about 1.0%. When the NTC content is 9.0%, the fiber content increases from 1.0% to 7.0% and the air voids increase by about 3.0%. When the fiber content is 1.0%, during the increase of the NTC content the air voids decreases by about 4.0% and the effect is more obvious. When the fiber content is 7.0%, the increase of NTC reduces the air voids by about 1.0% and the effect is greatly reduced. The increase of the asphalt-aggregate ratio makes the air voids of the mixture decrease first and then increase. When the asphalt-aggregate ratio is relatively high or low the influences of fiber and NTC are reduced, and when the asphalt-aggregate ratio is moderate the effects of the two are amplified.

Marshall Stability
The quadratic equation was selected as the design model for the fitting calculation. Then, the variance test was performed. The results are shown in Table 16.
The distribution of the predicted and measured values is shown in Figure 7. The distribution of the predicted and measured values is shown in Figure 7. It can be seen from the above figure that the fit between the predicted value and the measured value obtained by the equation is good, indicating that the fitting equation was selected accurately. Using the same method as Section 3.1.1, the contour plots and the response surface plots between the Marshall stability and each influence factor were drawn. The results are shown in Figure 8. It can be seen from the above figure that the fit between the predicted value and the measured value obtained by the equation is good, indicating that the fitting equation was selected accurately. Using the same method as Section 3.1.1, the contour plots and the response surface plots between the Marshall stability and each influence factor were drawn. The results are shown in Figure 8. It can be seen from the surface graph that the incorporation of fiber and NTC significantly improves the Marshall stability of the asphalt mixture, but as the amount of both increases the Marshall stability first increases and then decreases. When the content of the fiber or NTC is too high or too low, the influence of other modifiers on the Marshall stability will become greater. The asphaltaggregate ratio has a greater influence on the Marshall stability. The Marshall stability of the mixture will increase first and then decrease with the increase of the asphalt-aggregate ratio. It can be seen from the surface graph that the incorporation of fiber and NTC significantly improves the Marshall stability of the asphalt mixture, but as the amount of both increases the Marshall stability first increases and then decreases. When the content of the fiber or NTC is too high or too low, the influence of other modifiers on the Marshall stability will become greater. The asphalt-aggregate ratio has a greater influence on the Marshall stability. The Marshall stability of the mixture will increase first and then decrease with the increase of the asphalt-aggregate ratio.

Three Other Response Indicators
The three response indicators are analyzed above. The analysis process of the other three response indicators is similar, and their respective fitting equations are given here.
The fitting equation of the flow value is: Flow value = 0.90820 − 0.021562 × A − 0.029375 × B + 0.36354 × C The fitting equation of the VMA is: The fitting equation of the VMA is:

Input Index Optimization Based on Response Surface Fitting Model
It can be seen from the analysis of the variance of each response index that the influence of each input index on the response value can be fitted by the equation. It is difficult to obtain the optimal solutions of the three impact factors through the Marshall test. Therefore, it was necessary to set the range and expected value of each response index in the program. Then, the continuous surface model was fitted and a series of input index combinations were predicted.
The response output index range was selected with reference to the relevant index range in the Marshall test and the selection of each expected value was adjusted according to the situation. According to the regulations in JTG F40−2004 [42], the response output index range was selected. There is no requirement for the density range in JTG F40−2004 [42] and the range target value was set as the maximum density. The required range of air voids was 3.0-5.5%, and the target value was set at a median value of 4%. The target value of Marshall stability was set as the maximum value. The specified range of flow value was 2.0 mm-4.0 mm and the target value was set as the median value of 3.0 mm. The provisions of VMA were not less than 14%. The requirements of VFA were 65-75%. The expected value of the mixture with the highest Marshall stability and the best Marshall stability is the desired result. The predicted values after optimization under the optimal conditions are shown in Table 17.

Model Validation
In order to verify the accuracy of the prediction results of the response surface methodology, the optimum asphalt-aggregate ratio (5.67%) and the mixture ratio of BF (3.9%) and NTC (5.1%) were selected as the formed specimens. The sample preparation process was the same as Section 2.4.1. The Marshall results were compared and verified. There were 6 test pieces in each group in the verification tests. The results are shown in Table 18. According to the data obtained from the verification test and the prediction data of the response surface method, the errors of the predicted response value and the experimental value are shown in Table 18. The maximum error was 6.960% and the minimum error was 0.096%. This shows that in the research on the mix ratio of the asphalt concrete, it was feasible to use the response surface method to optimize the selection of the best impact factor.

Conclusions
Based on the response surface method, this research adopted the Box-Behnken design method to carry out a three−factor, three−level design. Density, air voids, Marshall stability, flow value, VMA, and VFA values were used as response indicators, while the asphalt-aggregate ratio, NTC, and fiber contents were taken as influencing factors. This method is more logical and predictable when determining the optimal content of various factors. This paper draws the following conclusions through experimental analysis: (1) Through the Box-Behnken design method, a three−factor function was fitted to each response indicator. It is concluded that the relationship between density, air voids, Marshall stability, VMA, VFA, and the three influencing factors is a quadratic function, and the relationship between the flow value and the three influencing factors is linear; (2) Through analysis of variance and testing of the fitted model, it is concluded that the fitted model has a higher significance level. Comparing the predicted value with the measured value, it is found that the difference between them is small. This shows that the model has high fitting accuracy and prediction ability, and can accurately reflect the relationship between each input index and output index; (3) Through comprehensive consideration of the response index values, the three influencing factors are optimized and predicted, and the optimal content of NTC, BF, and the optimal asphalt-aggregate ratio are 5.1%, 3.9%, and 5.67%, respectively. The response index value is obtained through model test prediction, while the measured data is compared with the predicted value to obtain a minimum relative error of 0.096% and a maximum error of 6.960%. This shows that the response surface method can be used to design, analyze, optimize, and predict the mix ratio of the asphalt mixture.