Next Article in Journal
Fabrication of Graphitized Carbon Fibers from Fusible Lignin and Their Application in Supercapacitors
Next Article in Special Issue
Effects of Bitumen Thickness on the Aging Behavior of High-Content Polymer-Modified Asphalt Mixture
Previous Article in Journal
Establishing a Framework for Fused Filament Fabrication Process Optimization: A Case Study with PLA Filaments
Previous Article in Special Issue
Review of the Application of Microwave Heating Technology in Asphalt Pavement Self-Healing and De-icing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Objective Optimization of Epoxy Resin Adhesive for Pavement Toughened by Self-Made Toughening Agent

1
CRCC Yunnan Investment Co., Ltd., Kunming 650299, China
2
Faculty of Civil Engineering and Mechanics, Kunming University of Science and Technology, Kunming 650500, China
*
Author to whom correspondence should be addressed.
Polymers 2023, 15(8), 1946; https://doi.org/10.3390/polym15081946
Submission received: 13 March 2023 / Revised: 15 April 2023 / Accepted: 18 April 2023 / Published: 19 April 2023
(This article belongs to the Special Issue Application of Polymer Materials in Pavement Design)

Abstract

:
Epoxy resin adhesive for pavement is often insufficient in flexibility and toughness. Therefore, a new type of toughening agent was prepared to overcome this shortcoming. To achieve the best toughening effect of a self-made toughening agent on an epoxy resin adhesive, its ratio to the epoxy resin needs to be optimally selected. A curing agent, a toughening agent, and an accelerator dosage were chosen as independent variables. The epoxy resin’s adhesive tensile strength, elongation at break, flexural strength, and flexural deflection were used as response values to establish a single-objective prediction model of epoxy resin mechanical property indexes. Response surface methodology (RSM) was used to determine the single-objective optimal ratio and analyze the effect of factor interaction on epoxy resin adhesive’s performance indexes. Based on principal component analysis (PCA), multi-objective optimization was performed using gray relational analysis (GRA) to construct a second-order regression prediction model between the ratio and gray relational grade (GRG) to determine the optimal ratio and to validate it. The results showed that the multi-objective optimization using response surface methodology and gray relational analysis (RSM-GRA) was more effective than the single-objective optimization model. The optimal ratio of epoxy resin adhesive was 100 parts of epoxy resin, 160.7 parts curing agent, 16.1 parts toughening agent, and 3.0 parts accelerator. The measured tensile strength was 10.75 MPa, elongation at break was 23.54%, the bending strength was 6.16 MPa, and the bending deflection was 7.15 mm. RSM-GRA has excellent accuracy for epoxy resin adhesive ratio optimization and can provide a reference for the epoxy resin system ratio optimization design of complex components.

Graphical Abstract

1. Introduction

The steel bridge deck pavement system is an important part of the steel box girder bridge, whose waterproof adhesive layer is critical. Its failure is one of the main factors causing the problem of steel deck pavement [1]. The failure of the waterproof adhesive layer leads to the separation or slippage between the pavement layer and the bridge deck. This allows the steel structure of the steel box girder bridge to be eroded by water, thus causing serious damage to the main structure of the bridge. Using materials with excellent mechanical properties, stability, and deformability in relation to the waterproof bonding layer is a critical measure to avoid adhesive failure and to prolong the life of the steel bridge deck pavement system [2]. At present, two kinds of waterproof adhesives for steel bridge deck pavement are commonly used, namely, epoxy resin adhesive and methacrylate adhesive [3,4]. Epoxy resin adhesive has excellent mechanical properties, high bonding strength, good impermeability, strong corrosion resistance, and other advantages [5,6], and it has been extensively studied.
Zeng et al. [7] studied the components of epoxy resin binder by single factor control variable method. With the increase in the amount of toughening agent, the elongation at break increased, but the tensile strength and elastic modulus decreased. Chen et al. [8] studied the effects of different raw-material compositions on the performance of epoxy resins for road surfaces. They determined reasonable material compositions, as well as ratios, and found that the epoxy resin raw-material composition had a significant impact on the performance of epoxy resin, and its bonding performance index was proportional to the mechanical performance index. Xi et al. [9] studied three types of epoxy resins for waterproof binding and found that the tensile strength and elongation at break data for all three were quite different. Liu et al. [10] studied the basic properties of epoxy resin adhesive for steel bridge deck pavement under different environmental conditions. The tensile strength of the self-made epoxy resin adhesive was 3.7 MPa, the elongation at break could reach 180%, and it had good high- and low-temperature performance, bonding performance, and water permeability resistance. Zhou et al. [11] summarized the results of relevant research on epoxy resin adhesives. They found that the elongation at break of epoxy resin adhesives with high tensile strength is low, and their strength and deformability are often negatively correlated. How to balance the interaction of each component and how to achieve the balance of strength and deformability, two important mechanical indexes, as well as how to prepare epoxy resin with high strength, strong deformability, and good toughness as a waterproof binder, have become the foci of this research.
Response surface methodology (RSM) is an optimization method, integrating experimental design and mathematical modeling proposed by mathematicians Box and Wilson [12], which can be used to solve multifactor and multilevel continuous response problems. Compared with orthogonal experimental and uniform designs, it has the advantage of high accuracy and can also analyze the interaction between influencing factors [13,14,15]. Epoxy resin binder has several performance indicators with different dimensions to be evaluated that are difficult to balance, which leads to the challenge of determining their optimal ratio. Gray relational analysis (GRA) is introduced to solve multi-objective response problems; it is suitable for solving those involving complex relationships among multiple objectives and factors and can optimize multi-objective responses. Compared with most scholars’ single-factor control variable method [16,17], the RSM–GRA method has the advantages of analyzing the interactions of multiple factors and obtaining the best material parameters quickly and accurately [18].
In this study, the dosage of the curing agent, self-made toughening agent, and accelerator were taken as independent variables, and the tensile strength, elongation at break, bending strength, and bending deflection of epoxy resin adhesive were taken as response values. The response surface optimization test was designed by the response surface methodology Box–Behnken design (RSM-BBD) method. RSM was used to establish the single-objective prediction model of each response. The influence of factor interactions on the response values was analyzed, and the single objective optimal ratio was determined. On this basis, the gray relational analysis (GRA) method was introduced to solve the problem of multi-objective response optimization, and the gray relational degree (gray relational grade (GRG)) prediction model was established to optimize the above indicators to obtain the optimal proportion scheme of epoxy resin binder with comprehensive properties. Our results provide a reference for the design of an epoxy resin system with complex components and may guide the further development and engineering application of self-made acrylate copolymer toughening agents.

2. Materials and Methods

2.1. Raw Materials

The epoxy resin used in the study was E-51 bisphenol A epoxy resin, whose epoxy equivalent is 184~195 g/eq, and its viscosity is 10,000~16,000 mPa∙s (25 °C), and this was produced by Baling Petrochemical Company, Sinopec., Yueyang, China. The curing agent was tung oil anhydride (TOA) with a viscosity of 5000~15,000 mPa∙s (25 °C), and the anhydride equivalent is 120, and it was produced by Shandong Jiaying Chemical Technology Co., Ltd., Qingdao, China. The accelerator was 2,4,6-tris (dimethylaminomethyl) phenol (DMP-30) with an amine value of 600~630 mg/g and viscosity of 100~300 mPa∙s (25 °C), which was produced by Changzhou Shanfeng Chemical Co., Ltd., Changzhou, China.

2.2. Synthetic Toughening Agent

The synthesis reaction was carried out in two steps. Step 1: polyether glycol (PTMG) was dehydrated and vacuumed at 110 °C for 2 h. Then, when the temperature was reduced to 60 °C, a certain amount of toluene diisocyanate (TDI) was added, and the molar ratio of TDI to PTMG was 2:1. The reaction was carried out at 80 °C for 2 h to obtain the NCO-terminated polyurethane prepolymer. Step 2: after the prepolymer was prepared, 1,4-butanediol diglycidyl ether (BDDGE) and a small amount of 2-ethyl-4-methylimidazole (EMI) were added successively, and the reaction temperature was kept at about 160 °C. Samples were taken every 1 h for infrared testing, and the reaction was continued until the NCO infrared absorption peak disappeared. The epoxy value of the final product was 0.23, and it was marked as the self-made epoxy resin toughening agent for pavements (SM-EPT).

2.3. Sample Preparation

The epoxy resin and the toughening agent were preheated in a vacuum drying oven (DZF-1, Beijing Yongguangming Medical Instrument Co., Ltd. Beijing, China) at 80 °C and dried in a vacuum for 30 min, then mixed, and stirred at 80 °C for 5 min at 1000 R/min to obtain component A. The curing agent and the accelerator were mixed and stirred for 5 min at 80 °C and 1000 R/min to obtain component B. The components A and B were mixed and stirred for 10 min at 60 °C and 500 R/min to obtain the epoxy resin binder. Then, the binder was placed in a vacuum drying oven, defoamed at 60 °C for 20 min, and poured into a PTFE mold. The geometric size of the specimen is shown in Figure 1. It was cured at 100 °C/12 h + 120 °C/12 h in a constant-temperature drying oven (101-0ES, Beijing Yongguangming Medical Instrument Co., Ltd. Beijing, China) and then naturally cooled to room temperature for demolding and performance testing. The prepared specimen is shown in Figure 2.

2.4. Test Method

The mechanical properties of the specimens were tested according to the test method for properties of resin castings (GB/T 2567-2021). An electronic universal testing machine (DDL100, Changchun Research Institute of Mechanical Sciences Co., Ltd., Changchun, China) was used to continuously record the tensile and bending loads. Extensometer (YYU-25/50, Steel Yanke Testing Technology Co., Ltd., Beijing, China) was used to measure the elongation within the gauge length of the tensile test piece. The deflection at the mid-span of the bending specimen was measured by a displacement meter (JC-LVDT type, Liyang Jincheng Testing Instrument Factory, Liyang, China). The loading speed and the bending tests were 10 mm/min. Five samples were tested in each group to obtain the average value.

2.5. Response Surface Experimental Design

Consulting with the relevant research results [19,20,21], the dosage range of the curing agent was 150~170%, the dosage range of toughening agent was 10~20%, and the dosage of the accelerator was 2~4%. Following the RSM-BBD method, the amount of curing agent, toughening agent, and accelerator (calculated by the mass percentage of epoxy resin, represented by A, B, and C, respectively) were taken as the independent variables, and the tensile strength ( σ t ), elongation at break ( ε t ), bending strength ( σ f ), and bending deflection ( S ) of epoxy resin adhesive were taken as the response values. A three-factor, three-level response surface test was designed. The factor code and level design are shown in Table 1. The data were analyzed by Design-Expert8.0 statistical software, and the response surface design and test results are shown in Table 2.

3. Results and Discussion

3.1. Construction of the Single-Objective Optimization Model

3.1.1. Model Construction and Model Verification

Multivariate quadratic regression fitting was carried out on the test data to obtain the regression models of σ t , ε t , σ f , and S, as shown in Equations (1)–(4).
σ t = 468.105 + 5.302 × A + 2.248 × B + 21.595 × C 0.002 × A × B 0.094 × A × C + 0.170 × B × C 0.015 × A 2 0.069 × B 2 1.159 × C 2
ε t = 1024.470 + 12.993 × A 5.158 × B + 32.626 × C + 0.032 × A × B 0.037 × A × C 0.315 × B × C 0.042 × A 2 0.053 × B 2 3.465 × C 2
σ f = 294.172 + 3.270 × A + 0.770 × B + 18.830 × C 0.003 × A × B 0.075 × A × C 0.031 × B × C 0.0095 × A 2 0.036 × B 2 1.004 × C 2
S = 198.850 + 2.251 × A 0.075 × B + 4.994 × C + 0.005 × A × B + 0.013 × A × C 0.122 × B × C 0.084 × A 2 + 0.0018 × B 2 0.824 × C 2
where σ t is tensile strength in MPa; ε t is elongation at break in %; σ f is bending strength in MPa; S is bending deflection in mm; A is the amount of curing agent,%; B is the amount of toughening agent in %; and C is the dosage of accelerator in %;
To explore the significance of the influence between the factors (independent variables) of the response surface regression model and the response values, the regression model was subjected to analysis of variance, as shown in Table 3. The model was evaluated using the F test. The larger the F value is, the smaller the p value is, which means that the probability of the invalid hypothesis of the model is also smaller, and the model is more significant [22]. The F values of the regression models σ t , ε t , σ f , and S were 66.50, 49.57, 46.32, and 40.49, respectively, and the corresponding p values were less than 0.0001, indicating that the statistical significance of the four models was very high. The p values, corresponding to the F values of the four regression models, were greater than 0.05, indicating that the lack of fit caused by error was not significant. The four models could well describe the relationship between response and factors.
Table 3 shows the relationship between three independent variables (factors) and materials’ response values. It can be seen that the order of influence of the three factors A (curing agent dosage), B (toughening agent dosage), and C (accelerator dosage) on the tensile strength model was A > B > C, the order of influence of the three factors on the bending strength model was A > C > B, and factor A was the most significant in both models. The results show that the amount of curing agent was the main factor affecting the strength of epoxy resin adhesive. A, B, and C influenced the elongation at break model in the order of B > C > A and influenced the bending deflection model in the order of B > A > C. Factor B was the most significant factor in both models, which indicates that the content of toughening agent was the main factor affecting the deformation ability of epoxy resin adhesive. In the two models of tensile strength and flexural strength, AC was significant (p < 0.05), while AB and BC were not (p > 0.05), indicating that the interaction of curing agent dosage A and accelerator dosage C had significant effects on tensile strength and flexural strength. The order of factor interaction significance in the elongation at break model was AB > BC > AC, and the order of factor interaction significance in the bending deflection model was BC > AB > AC.
Table 4 shows the statistical analysis of the fitting accuracy of the four regression models σ t , ε t , σ f , and S . The R2 of the four models is close to 1, indicating that the correlation between the predicted value and the actual value of the four regression models was good. The difference between the calibration coefficient of determination (R2Adj) and the prediction coefficient of determination (R2Pred) of the four models was less than 0.2, and the coefficient of variation (CV) for all of them was less than 10%. The signal-to-noise ratio was far greater than 4, which further shows that the fitting error of the four regression models σ t , ε t , σ f , and S was small, and the model fitting effect was good [23,24]. Figure 3 shows the comparison results of actual values and predicted values of models σ t , ε t , σ f , and S . The predicted values of the four regression models are close to the actual values, and the average deviations between the actual values and the predicted values are 1.68%, 2.23%, 2.57%, and 1.60%, respectively, indicating that the reliability of the model fitting is high. The above analysis shows that the four regression models can accurately describe the functional relationship between the response values and the factors and analyze and predict the test results.

3.1.2. Analysis of Three-Dimensional Response Surface Interaction Effect

According to the test results and the variance analysis results in Table 3, the tensile strength σ t , elongation at break ε t , bending strength σ f , and bending deflection S of the epoxy resin adhesive were all affected by the interaction between the factors.
In the regression model of tensile strength σ t , the F value of the interaction term AC was the largest, and the p value was the smallest, which indicates that the interaction of curing agent dosage A and accelerator dosage C had the most significant effect on the tensile strength of epoxy resin adhesive. Figure 4 is a three-dimensional response surface showing the effect of the interaction between the curing agent dosage A and the accelerator dosage C on the tensile strength ( σ t ), with the toughening agent dosage B of 15%. When the amount of accelerator was constant, the tensile strength increased first and then decreased with the increase in the amount of curing agent. The main reason is that when the amount of curing agent is too small, the proportion of epoxy is too high, and a large number of ether bonds are formed in the reaction process, resulting in many irreversible cross-linking points. Thus, the transesterification reaction does not occur easily in the curing system. When the cross-linking density and uniformity of the curing system were reduced, strength decreased [25,26]. When the amount of curing agent was constant, the tensile strength of epoxy resin binder increased first and then decreased with the increase in the amount of accelerator, which is consistent with the research results of Liang Ming et al. [19]. The reason is that a proper amount of DMP-30 accelerator can catalyze the anhydride group, promote the curing and cross-linking of epoxy resin, and relatively improve its mechanical properties. However, the excessive amount of accelerator will make the curing system release too much heat in unit time, resulting in phase separation of the curing system, reducing the cross-linking density, and, thus, reducing the mechanical properties of the resin [27,28,29]. In summary, simultaneously increasing the amount of curing agent and accelerator within a certain range can effectively improve the tensile strength of epoxy resin binder.
Figure 5a is a three-dimensional response surface, showing the effect of the interaction between the curing agent dosage A and the toughening agent dosage B on the elongation at break ( ε t ), with the accelerator dosage C of 3%. Figure 5b is a three-dimensional response surface of the effect of the interaction between the toughening agent content B and the accelerator content C on the elongation at break ( ε t ), with a curing agent content of 160%. When the dosage of toughening agent B was constant, the elongation at break increased first and then decreased with the increase in the dosage of curing agent A or the dosage of accelerator C. When the dosage of curing agent A or accelerator C was constant, the elongation at break increased with the increase in the dosage of the toughening agent. The variance analysis of the model ε t results from Table 3 obtained the p values of 0.0036 and 0.0037 for AB and BC, respectively, indicating that the interaction of toughening agent B and curing agent A, as well as the interaction between toughening agent B and accelerator C, had significant effects on the elongation at break. The interaction between toughening agent dosage B and curing agent dosage A had the most significant effect on the elongation at break. A and B, as well as B and C interactions, had significant effects on elongation at break, which may be due to the decrease in cross-link density of the curing system when the toughening agent is used, but the increase in cross-link density can be caused by the addition of proper amounts of curing agent and accelerator. The interactions between the amount of toughening agent and the amount of curing agent, and between the amount of toughening agent and the accelerator, significantly affected the elongation at break [30,31]. Thus, in summary, to improve the elongation at break of epoxy resin adhesive, it is necessary to comprehensively consider the effects of the interaction of the amounts of toughening agent and curing agent, as well as the quantities of toughening agent and accelerator on the elongation at break of epoxy resin adhesive in the process of ratio optimization.
Figure 6 is a three-dimensional response surface with the effect of the interaction of the curing agent dosage A and the accelerator dosage C on the flexural strength ( σ f ), with the toughening agent dosage B of 15%. It can be seen from Figure 6 that, with the increase in the amount of curing agent or accelerator, the bending strength of the epoxy resin binder increased first and then decreased. According to the variance analysis results of the σ f model in Table 3, the AC term was significant (p = 0.0002 < 0.05), and the response surface was steep when the curing agent dosage A and the accelerator dosage C changed at the same time, indicating that the interaction of the curing agent dosage A and the accelerator dosage C had a significant effect on the bending strength. The reason is that the tertiary amine produced by the decomposition of the DMP-30 accelerator can react with the anhydride group in the anhydride curing agent to produce carboxylate anion, which can catalyze the anhydride curing agent and make the system more easily cross-linked and cured. However, the excessive accelerator will quickly cross-link and cure the system, resulting in the unreacted chain segment that cannot continue participating in the reaction. Macroscopically, the interaction between the amount of DMP-30 accelerator and the amount of anhydride curing agent affects the mechanical properties of epoxy resin adhesive [19,27,28].
Figure 7 is a three-dimensional response surface with the effect of the interaction of the toughening agent dosage B and the accelerator dosage C on the bending deflection ( S ) of the epoxy resin adhesive with the curing agent dosage A of 160%. With the increase in the content of toughening agent, the range of change in bending deflection was greater when the content of accelerator was 2% than when it was 4%, indicating that the sensitivity of bending deflection of epoxy resin adhesive to the amount of toughening agent decreased with the increase in the amount of accelerator. According to the variance analysis results of the bending deflection S model in Table 3, the BC term was significant (p = 0.001 < 0.05), indicating that the interaction of the toughening agent B and the accelerator C had a significant effect on the bending deflection.

3.1.3. Single-Objective Optimization Results

The single-objective response was optimized by Design-Expert 8.0 software. The maximum tensile strength of epoxy resin adhesive was 10.61 MPa, the corresponding amount of curing agent was 163.1%, the amount of toughening agent was 14.3%, and the amount of accelerator was 2.8%. The maximum elongation at break of the epoxy resin adhesive was 24.60%, and the corresponding amounts of curing agent, toughening agent, and accelerator were 156.9%, 18.8%, and 3.1%, respectively. The maximum bending strength of the epoxy resin adhesive was 6.17 MPa, and the amounts of curing agent, toughening agent, and accelerator were 162.5%, 14.9%, and 3.1%, respectively. The maximum bending deflection of epoxy resin adhesive was 7.42 mm, the corresponding amount of curing agent was 159.0%, the amount of toughening agent was 20.0%, and the amount of accelerator was 2.8%.

3.2. Construction of Multi-Objective Optimization Model

3.2.1. Calculation Process of Gray Correlation Degree

Gray relational analysis (GRA) can convert multiple optimization objectives into gray relational values by reducing dimensions and then optimizing the gray relational values [32,33]. The larger the gray correlation value is, the better the corresponding response is [34]. The calculation process of the gray correlation degree is as follows:
(1) Normalization. The tensile strength ( σ t ), elongation at break ( ε t ), bending strength ( σ f ), and bending deflection ( S ) of the epoxy resin adhesive were normalized to eliminate the effect of dimensions on the analysis. The bigger all four are, the better, and the normalization formula is shown in Equation (5).
N = y min ( y ) max ( y ) min ( y )
where N is the normalized value of each response, max ( y ) is the maximum value of the actual response, min ( y ) is the minimum value of the actual response, and y is the actual value of each group of tests.
(2) Calculation of gray correlation coefficient. The gray relational coefficient (GRC) represents the relationship between the test result and the optimal solution [35], and the calculation formulas are shown in Equations (6) and (7)
GRC = Δ min + ξ Δ max Δ + ξ Δ max
Δ = 1 N
where Δ represents the deviation sequence, ξ is the judgment coefficient, ξ [ 0 , 1 ] , and, in this study, ξ is 0.5.
(3) Response weight calculation. To obtain the gray correlation coefficient of each response, it is necessary to calculate the influence weight. Principal component analysis (PCA) quantitatively analyzes the weight of the contribution rate of each target to the response by reducing the dimensionality [36]. With the help of the PCA analysis module of Minitab software, the influence weights of tensile strength σ t , elongation at break ε t , bending strength σ f , and bending deflection S on gray correlation degree are calculated.
(4) Calculation of gray correlation degree. Gray relational degree (GRG) is the weighted sum of gray relational coefficients. The higher the gray relational degree is, the better the corresponding response is. The calculation formula is presented in Equation (8)
GRG = i = 1 n β i GRC
where i = 1 n β i = 1 , and β i is the weight of the i th response, calculated by PCA.

3.2.2. Calculation Result and Analysis of Gray Correlation Degree

The test results are normalized by Equation (5), and GRC is calculated by Equations (6) and (7). The influence weights of tensile strength σ t , elongation at break ε t , bending strength σ f , and bending deflection S on the gray correlation degree are further calculated by PCA. The results are shown in Table 5. Table 6 shows the results of the gray correlation degree calculated according to Equation (8).

3.2.3. Construction and Optimization of GRG Response Model

In order to optimize the proportion of epoxy resin binder, the mapping relationship between curing agent dosage A, toughening agent dosage B, accelerator dosage C, and GRG must be established. In this paper, Design-Expert8.0 software was used to establish the second-order mathematical prediction model of GRG, as shown in Equation (9). The comparison between the actual and predicted values of GRG is shown in Figure 8. It can be seen that all data are evenly distributed on a straight line and on both sides, and the average deviation between the fitting value and the actual value was 3.42%, indicating that the fitting effect of the GRG prediction model is good. According to the variance analysis results of the GRG prediction model in Table 7, the model had p < 0.0001, indicating that the model was highly significant. The R2 of the GRG prediction model was 98.39%, the difference between R2Adj and R2Pred was less than 0.2, and the coefficient of variation CV was less than 5%, indicating that the model had high reliability and good fitting degree and could be used for subsequent prediction and optimization.
GRG = 68.414 + 0.811 × A + 0.100 × B + 2.204 × C + 5.211 × 10 4 × A × B 5.963 × 10 3 × A × C 3.427 × 10 3 × B × C 2.495 × 10 3 × A 2 5.388 × 10 3 × B 2 0.199 × C 2
To reflect the influence of the interaction of each factor on the GRG response value, the three-dimensional response surface of GRG under the interaction of different factors was established using the mathematical prediction model, as shown in Figure 9. The shape of the response surface of the GRG model under the interaction of different factors is a quadratic paraboloid with a downward opening, indicating a maximum value of GRG in the test range. According to the variance analysis results of the GRG response model in Table 7, the influence of the three factors on the GRG model was in the order of B > A > C. The interaction of curing agent dosage A and accelerator dosage C had the most significant influence on GRG. The optimized GRG response surface model was analyzed using Design-Expert 8.0 software, and the optimal gray correlation degree was 0.885, corresponding to a curing agent dosage of 160.7%, toughening agent dosage of 16.1%, and accelerator dosage of 3.0%.

3.3. Experimental Verification

Through the analysis of the results, the proportion of the optimal gray correlation degree ( GRG max ) was determined, and the proportion corresponding to the optimal gray correlation degree ( GRG max ) was selected for comparison with the proportion corresponding to the maximum tensile strength σ t - max , the maximum elongation at break ε t - max , the maximum bending strength σ f - max , and the maximum bending S max deflection obtained by the single-objective prediction model.
The ratio obtained from the single-objective prediction model in Section 3.1.3 and the ratio obtained from the multi-objective optimization model in Section 3.2.3 were experimentally verified, and the experimental verification results are shown in Table 8. Figure 10 shows the stress–strain curve for the tensile test and the load-deflection curve for the bending test. The analysis revealed that the tensile strength obtained by the ratio σ t - max was the largest, the elongation at break obtained by the ratio ε t - max was the largest, the bending strength obtained by the ratio σ f - max was the maximum, and the bending deflection obtained by the ratio S max was the maximum. It can be found that the actual values of tensile strength, elongation at break, bending strength, and bending deflection deviate from the predicted values by 2.7%, 2.0%, 2.8%, and 1.1%, respectively, with actual values higher than the predicted values. This proved the feasibility of the single objective prediction model.
GRG max was compared to σ t - max , elongation at break increased by 19.31%, bending strength increased by 1.81%, and bending deflection increased by 31.43%. GRG max was compared to ε t - max , tensile strength increased by 24.70%, bending strength increased by 16.22%, and bending deflection increased by 4.38%. GRG max was compared to σ f - max , the tensile strength was increased by 5.29%, the elongation at break was increased by 10.41%, and the bending deflection was increased by 27.00%. GRG max was compared to S max , the tensile strength was improved by 30.46%, the elongation at break was improved by 6.04%, and the bending strength was improved by 22.95%.
The gray relational degrees of GRG max , σ t - max , ε t - max , σ f - max , and S max were 0.830, 0.794, 0.503, 0.626, and 0.408, respectively, and the gray relational degree of GRG max was the highest. The results showed that the optimum ratio of epoxy resin binder was obtained by multi-objective optimization of response surface methodology and gray relational grade. The optimum ratio was 160.7% for the curing agent, 16.1% for the toughening agent, and 3.0% for the accelerator.

4. Conclusions

(1)
The amount of curing agent is the main factor affecting the strength of epoxy resin adhesive, and the amount of toughening agent is the main factor affecting its deformability. The interaction of curing agent dosage and accelerator dosage has the most significant effect on the tensile strength and bending strength of epoxy resin binder. The interaction of toughening agent dosage and curing agent dosage has the most significant effect on the elongation at break of the epoxy resin binder, and the interaction of toughening agent dosage and accelerator dosage has a significant effect on bending deflection.
(2)
The optimal ratio of epoxy resin binder obtained by RSM-GRA multi-objective optimization was a curing agent dosage of 160.7%, toughening agent dosage of 16.1%, and an accelerator dosage of 3.0%. An epoxy resin binder with a tensile strength of 10.75 MPa, elongation at break of 23.54%, bending strength of 6.16 MPa, and a bending deflection of 7.15 mm can be prepared with this ratio. The test results show that the RSM-GRA multi-objective optimization model is accurate, effective, and has future application significance for optimizing the epoxy resin binder ratios.
(3)
Compared with the single-objective optimization model, in this study, the RSM-GRA multi-objective optimization model was used to obtain the largest gray relational grade (GRG) of the performance indicators of the epoxy resin binder, and the corresponding epoxy resin binder had the best comprehensive performance. The RSM-GRA multi-objective optimization method used in this paper can not only optimize the proportion of epoxy resin binder, but it also provide a reference for the proportion optimization of other complex epoxy resin components.

Author Contributions

H.B.: conceptualization, methodology, writing—original draft, and investigation. L.G.: data curation and writing—review and editing. H.H.: conceptualization, methodology, data curation. S.Z.: conceptualization, methodology. Z.L.: funding acquisition and project administration, supervision. All authors have read and agreed to the published version of the manuscript.

Funding

Science and Technology Project of Transportation Department of Yunnan Province: Science and Education Division [2021] No.39.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank all project research team members for their help with technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Z.Q.; Tao, J.; Zhang, S.T. Experiment and evaluation on performance of epoxy asphalt waterproof cohesive layer on bridge deck pavement. J. Chang. Univ. (Nat. Sci. Ed.) 2011, 31, 1–6. (In Chinese) [Google Scholar]
  2. Xu, Y.; Lv, X.; Ma, C.; Liang, F.; Qi, J.; Chou, Z.; Xu, S. Shear Fatigue Performance of Epoxy Resin Waterproof Adhesive Layer on Steel Bridge Deck Pavement. Front. Mater. 2021, 7, 469. [Google Scholar] [CrossRef]
  3. Kumar, P.; Patnaik, A.; Chaudhary, S. A review on application of structural adhesives in concrete and steel-concrete composite and factors influencing the performance of composite connections. Int. J. Adhes. Adhes. 2017, 77, 1–14. [Google Scholar] [CrossRef]
  4. Zhang, M.; Hao, P.; Men, G.; Liu, N.; Yuan, G. Research on the compatibility of waterproof layer materials and asphalt mixture for steel bridge deck. Constr. Build. Mater. 2021, 269, 121346. [Google Scholar] [CrossRef]
  5. Xiang, Q.; Xiao, F. Applications of epoxy materials in pavement engineering. Constr. Build. Mater. 2020, 235, 117529. [Google Scholar] [CrossRef]
  6. Paluvai, N.R.; Mohanty, S.; Nayak, S.K. Synthesis and Modifications of Epoxy Resins and Their Composites: A Review. Polym. Plast. Technol. Eng. 2014, 53, 1723–1758. [Google Scholar] [CrossRef]
  7. Zeng, G.D.; Wen, G.X.; Li, J.C.; Ouyang, T.Z.; Yuan, M.; Zhou, Y.; Li, Q.; Wu, C.F. Development of waterproof binder for road and bridge curing flexible epoxy resin at room temperature. Highway 2021, 66, 310–315. (In Chinese) [Google Scholar]
  8. Chen, Q.; Lu, Y.; Wang, C.; Han, B.; Fu, H. Effect of raw material composition on the working performance of waterborne epoxy resin for road. Int. J. Pavement Eng. 2022, 23, 2380–2391. [Google Scholar] [CrossRef]
  9. Xi, L.; Huang, W.R. Study on the application of epoxy resin in waterproofing layer of bridge deck pavement. J. China Foreign Highw. 2016, 36, 292–295. (In Chinese) [Google Scholar]
  10. Liu, P.; Liu, Y.G.; Hao, Z.H.; Sheng, X.Y.; Li, L. Effect of environmental conditions on bonding performance of self-made epoxy resin adhesive. Appl. Chem. Ind. 2019, 48, 2032–2035. (In Chinese) [Google Scholar]
  11. Zhou, L.; Zhang, D.; Li, X.; Gao, Z.; Chen, Q.; Wang, C. Overview: Application of Resin Waterproof Adhesive Materials in Bridge Deck Pavement in China. Adv. Civ. Eng. 2022, 2022, 2320374. [Google Scholar] [CrossRef]
  12. Box, G.E.; Wilson, K.B. On the experimental attainment of optimum conditions. In Breakthroughs in Statistics: Methodology and Distribution; Springer: Berlin/Heidelberg, Germany, 1992; pp. 270–310. [Google Scholar]
  13. Li, L.; Zhang, Q.; He, Q.; Hu, X.B. Application of response surface methodology in experiment design and optimization. Res. Explor. Lab. 2015, 34, 41–45. (In Chinese) [Google Scholar]
  14. Cheng, J.L.; Zheng, M.; Lou, J.Q. Comparison of several common optimal experimental design methods. Res. Explor. Lab. 2012, 31, 7–11. (In Chinese) [Google Scholar]
  15. Li, Z.; Lu, D.; Gao, X. Optimization of mixture proportions by statistical experimental design using response surface method—A review. J. Build. Eng. 2021, 36, 102101. [Google Scholar] [CrossRef]
  16. Chen, X.Y.; Huang, W.D.; Zhang, W.J.; Lai, Z.P.; Lian, G.F. Multiple targets technology optimization based grey relative analysis of 18Ni300 die steel formed by selective laser melting. Chin. J. Lasers 2020, 47, 341–351. (In Chinese) [Google Scholar]
  17. Singh, O.P.; Kumar, G.; Kumar, M. Role of Taguchi and grey relational method in optimization of machining parameters of different materials: A review. Acta Electron. Malays. (AEM) 2019, 3, 19–22. [Google Scholar] [CrossRef]
  18. Deshmukh, S.S.; Jadhav, V.S.; Shrivastava, R. Review on Single and Multi-objective Optimization Process Parameters of EDM Using Taguchi Method and Grey Relational Analysis. In Proceedings of the 9th International Conference of Materials Processing and Characterization (ICMPC), Hyderbad, India, 8–10 March 2019; pp. 3856–3866. [Google Scholar]
  19. Liang, M.; Su, L.P.; Qiu, Z.M.; Xin, X.; Yao, Z.Y.; Ma, C.Y.; Ding, X.M. Effects of DMP-30 on curing kinetics and mechanical properties of epoxy resin/anhydride system. J. China Univ. Pet. (Ed. Nat. Sci.) 2021, 45, 175–181. (In Chinese) [Google Scholar]
  20. Fu, P.; Tan, X.; Xiao, L.H.; Nie, X.A.; Huang, J.R.; Zhang, L. Study on high-performance tung oil-based anhydride epoxy curing agent. Resin 2021, 36, 22–30. (In Chinese) [Google Scholar]
  21. Zou, S.L.; Jiang, S.F.; Shen, R.D. Properties of phenolic/epoxy resin curing system modified with acrylic block copolymer. Adhesion 2016, 37, 51–56. (In Chinese) [Google Scholar]
  22. Wang, J.W.; Wang, W. Response surface based multi-objective optimization of basalt fiber reinforced foamed concrete. Mater. Rep. 2019, 33, 4092–4097. (In Chinese) [Google Scholar]
  23. He, D.P.; Pan, Z.Q.; Wang, H.G.; Du, W.X. Optimizing the content of activated crumb rubber and SBR composite modified asphalt by response surface methodology-grey relation analysis. N. Chem. Mater. 2022, 50, 187–191. (In Chinese) [Google Scholar]
  24. Van Den Bergh, D.; Van Doorn, J.; Marsman, M.; Draws, T.; van Kesteren, E.-J.; Derks, K.; Dablander, F.; Gronau, Q.F.; Kucharský, Š.; Gupta, A.R.K.N.; et al. A Tutorial on Conducting and Interpreting a Bayesian ANOVA in JASP. Annee Psychol. 2020, 120, 73–96. [Google Scholar] [CrossRef]
  25. Zhou, D.W.; Zhai, J.X.; Li, S.; Quan, X.J.; Huo, F.; Li, G. Study on curing behavior and solubility of epoxy/anhydride/DMPA system. Chem. Res. Appl. 2017, 29, 1760–1764. (In Chinese) [Google Scholar]
  26. Smallenburg, F.; Leibler, L.; Sciortino, F. Patchy particle model for vitrimers. Phys. Rev. Lett. 2013, 111, 188002. [Google Scholar] [CrossRef] [PubMed]
  27. Han, Y.; Wang, Z.; Zhao, S.; Wang, J. AC impedance function of electrochemical working station as novel curing degree monitor method: A model curing system of epoxy/anhydride/DMP-30. Measurement 2019, 145, 600–610. [Google Scholar] [CrossRef]
  28. Huang, C.; Ge, Z.; Zhao, B.B.; Wang, Z.; Luo, Y.J. Effects of DMP-30 on curing behavior of epoxy resin/maleicanhydride systems. J. Chem. Eng. Chin. Univ. 2017, 31, 197–204. (In Chinese) [Google Scholar]
  29. Fan, M.; Liu, J.; Li, X.; Jue, C.; Zhang, J. Curing behaviors and properties of an extrinsic toughened epoxy/anhydride system and an intrinsic toughened epoxy/anhydride system. Thermochim. Acta 2013, 554, 39–47. [Google Scholar] [CrossRef]
  30. Liu, Z.L.; Li, H.F.; Gu, J.Y.; Wang, D.Z.; Qu, C.Y.; Yang, H.D. Properties of epoxy resin modified with acrylic block copolymer. Polym. Mater. Sci. Eng. 2018, 34, 24–31+36. (In Chinese) [Google Scholar]
  31. Heng, Z.; Chen, Y.; Zou, H.; Liang, M. Simultaneously enhanced tensile strength and fracture toughness of epoxy resins by a poly(ethylene oxide)-block-carboxyl terminated butadiene-acrylonitrile rubber dilock copolymer. RSC Adv. 2015, 5, 42362–42368. [Google Scholar] [CrossRef]
  32. Huang, W.; Mao, X.; Wu, Q.; Zhang, J. Experimental Investigation on the Shear Characteristics of Frozen Silty Clay and Grey Relational Analysis. Sustainability 2023, 15, 180. [Google Scholar] [CrossRef]
  33. Wen, P.; Wang, C.; Chen, M.; Chai, Z. Engineering Property Evaluation and Multiobjective Parameter Optimization of Argillaceous Gangue–Filled Subgrade Based on Grey Relational Analysis. J. Mater. Civ. Eng. 2023, 35, 04023007. [Google Scholar] [CrossRef]
  34. Hong, Q.; Shi, Y.Y.; Lu, D.N.; Guo, Y.M. Multi-response parameter optimization for the composite tape winding process based on grey relational analysis and response surface methodology. Acta Mater. Compos. Sin. 2019, 36, 2822–2832. (In Chinese) [Google Scholar]
  35. Xie, Z.; Liu, X.; Liu, Z.; Chen, P.; Lai, B.; Zhan, B.; Lao, J. Human Reliability Analysis of Virus Detection Equipment Based on Entropy Weighting Method and GRA. In Proceedings of the 7th International Conference on Condition Monitoring of Machinery in Non-Stationary Operations (CMMNO), Guangzhou, China, 11–13 June 2021; pp. 263–266. [Google Scholar]
  36. Jiang, R.C.; Ci, S.K.; Liu, D.W.; Sun, H.X.; Wang, D.F. Ply optimization of carbon fiber reinforced plastic control arm based on grey relational analysis. Acta Mater. Compos. Sin. 2022, 39, 390–398. (In Chinese) [Google Scholar]
Figure 1. The geometric size of the specimen: (a) tensile test specimen, (b) bending test specimen.
Figure 1. The geometric size of the specimen: (a) tensile test specimen, (b) bending test specimen.
Polymers 15 01946 g001
Figure 2. The prepared specimen: (a) tensile test specimen, (b) bending test specimen.
Figure 2. The prepared specimen: (a) tensile test specimen, (b) bending test specimen.
Polymers 15 01946 g002
Figure 3. Comparison of actual and prediction values: (a) model σ t , (b) model ε t , (c) model σ f , and (d) model S .
Figure 3. Comparison of actual and prediction values: (a) model σ t , (b) model ε t , (c) model σ f , and (d) model S .
Polymers 15 01946 g003
Figure 4. Response surface of the tensile strength model ( σ t ) under the interaction of factors A and C.
Figure 4. Response surface of the tensile strength model ( σ t ) under the interaction of factors A and C.
Polymers 15 01946 g004
Figure 5. Response surface of elongation at break ( ε t ) under the (a) interaction of factors A and B, as well as the (b) interaction of factors A and C.
Figure 5. Response surface of elongation at break ( ε t ) under the (a) interaction of factors A and B, as well as the (b) interaction of factors A and C.
Polymers 15 01946 g005
Figure 6. Response surface of bending strength ( σ f ) under the interaction of factors A and factor C.
Figure 6. Response surface of bending strength ( σ f ) under the interaction of factors A and factor C.
Polymers 15 01946 g006
Figure 7. Response surface of bending deflection ( S ) under the interaction of factor B and factor C.
Figure 7. Response surface of bending deflection ( S ) under the interaction of factor B and factor C.
Polymers 15 01946 g007
Figure 8. Comparison of actual and prediction values of GRG.
Figure 8. Comparison of actual and prediction values of GRG.
Polymers 15 01946 g008
Figure 9. GRG three-dimensional response surface under the interaction of factors (a) A and B, (b) A and C, and (c) C and B.
Figure 9. GRG three-dimensional response surface under the interaction of factors (a) A and B, (b) A and C, and (c) C and B.
Polymers 15 01946 g009
Figure 10. (a) Stress–strain curve of tensile test, (b) load–deflection curve of bending test.
Figure 10. (a) Stress–strain curve of tensile test, (b) load–deflection curve of bending test.
Polymers 15 01946 g010
Table 1. Test factors and levels.
Table 1. Test factors and levels.
LevelABC
Curing Agent/%Toughening Agent/%Accelerator/%
−1150102
0160153
1170204
Table 2. Response surface test design and test results.
Table 2. Response surface test design and test results.
NO.LevelTest Results
ABCTensile Strength/MPaElongation at Break/%Flexural Strength/MPaFlexural Deflection/mm
11107.7022.585.036.30
2−1106.1621.353.576.60
30−117.7720.344.935.28
400010.7021.686.265.86
5−10−16.2714.702.704.50
6−1017.9616.994.654.70
70−1−18.4713.793.853.62
800010.2221.666.205.93
91−108.4014.164.673.94
1000010.4023.175.966.32
1100010.3922.426.095.80
1210−19.4212.825.143.50
1300010.5022.135.955.55
14−1−106.4719.263.705.20
1501−17.1822.973.757.00
161017.3413.644.084.20
170116.8223.234.216.22
Table 3. Results of ANOVA of response surface model.
Table 3. Results of ANOVA of response surface model.
Model σ t Model ε t Model σ f Model S
FpFpFpFp
Model66.50<0.000149.57<0.000146.32<0.000140.49<0.0001
A65.71<0.000119.060.003352.650.000223.360.0019
B19.280.0032117.34<0.00010.990.3526162.85<0.0001
C3.840.090922.650.002116.810.00467.900.0261
AB0.560.480418.440.00361.370.28054.600.0692
AC51.890.00020.990.351851.590.00021.250.3009
BC0.420.536718.210.00372.190.182529.700.0010
A2145.06<0.0001137.38<0.000185.83<0.000159.780.0001
B2182.64<0.000113.820.007578.29<0.000117.900.0039
C282.52<0.000193.06<0.000196.58<0.000156.980.0001
Lack of Fit3.830.11351.930.26703.950.10860.470.9113
Table 4. Fit statistics for response surfaces of the model.
Table 4. Fit statistics for response surfaces of the model.
ModelR2R2AdjR2PredCV/%Adeq. Precision
σ t 0.98840.97360.85813.1321.516
ε t 0.98460.96470.84413.8318.947
σ f 0.98350.96230.79594.4122.049
S 0.98120.95690.93984.2019.992
Table 5. Weightiness of response.
Table 5. Weightiness of response.
Principal ComponentEigenvalueWeightiness/%
σ t 2.39359.813
ε t 1.48237.038
σ f 0.0892.221
S 0.0370.928
Sum4100
Table 6. Gray correlation degree calculation results.
Table 6. Gray correlation degree calculation results.
NO.GRCGRG
σ t ε t σ f S
10.4310.8890.5910.7140.607
20.3330.7350.3980.8140.488
30.4370.6430.5720.5040.517
41.0000.7711.0000.6060.911
50.3390.3790.3330.4120.354
60.4530.4550.5250.4320.455
70.5040.3550.4250.3410.446
80.8250.7680.9670.6210.806
90.4970.3650.5280.3640.447
100.8830.9890.8560.7200.920
110.8800.8650.9130.5930.873
120.6390.3330.6140.3330.523
130.9190.8260.8520.5470.879
140.3490.5670.4100.4930.433
150.3920.9520.4151.0000.606
160.4030.3520.4490.3850.385
170.3691.0000.4650.6920.608
Table 7. ANOVA results for response surface model of GRG.
Table 7. ANOVA results for response surface model of GRG.
SourceSum of SquaresMean SquareFp
Model0.6100.06847.64<0.0001
A6.709 × 10−36.709 × 10−34.700.0668
B0.0270.02718.980.0033
C1.628 × 10−41.628 × 10−40.110.7455
AB2.715 × 10−32.715 × 10−31.900.2103
AC0.0140.0149.960.0160
BC1.175 × 10−31.175 × 10−30.820.3945
A20.260.26183.57<0.0001
B20.0760.07653.500.0002
C20.170.17116.85<0.0001
Error9.991 × 10−31.428 × 10−3--
Lack of Fit1.820 × 10−36.068 × 10−40.300.8269
R2 = 0.9839R2Adj = 0.9633R2Pred = 0.9327CV/% = 6.26
Table 8. Response optimal ratio experimental comparison table.
Table 8. Response optimal ratio experimental comparison table.
NO.A/%B/%C/%Tensile Strength/MPaElongation at Break/%Flexural Strength/MPaFlexural Deflection/mmGRG
GRG max 160.716.13.010.7523.546.167.150.830
σ t - max 163.114.32.810.9019.736.055.440.794
ε t - max 156.918.83.18.6225.105.306.850.503
σ f - max 162.514.93.110.2121.326.345.630.626
S max 159.020.02.88.2422.205.017.500.408
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ba, H.; Guo, L.; Huan, H.; Zhang, S.; Lin, Z. Multi-Objective Optimization of Epoxy Resin Adhesive for Pavement Toughened by Self-Made Toughening Agent. Polymers 2023, 15, 1946. https://doi.org/10.3390/polym15081946

AMA Style

Ba H, Guo L, Huan H, Zhang S, Lin Z. Multi-Objective Optimization of Epoxy Resin Adhesive for Pavement Toughened by Self-Made Toughening Agent. Polymers. 2023; 15(8):1946. https://doi.org/10.3390/polym15081946

Chicago/Turabian Style

Ba, Huaiqiang, Luxin Guo, Haiyang Huan, Shibo Zhang, and Zhiwei Lin. 2023. "Multi-Objective Optimization of Epoxy Resin Adhesive for Pavement Toughened by Self-Made Toughening Agent" Polymers 15, no. 8: 1946. https://doi.org/10.3390/polym15081946

APA Style

Ba, H., Guo, L., Huan, H., Zhang, S., & Lin, Z. (2023). Multi-Objective Optimization of Epoxy Resin Adhesive for Pavement Toughened by Self-Made Toughening Agent. Polymers, 15(8), 1946. https://doi.org/10.3390/polym15081946

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

Article Metrics

Back to TopTop