Abstract

The advancement of technology and the increasing growth of environmental threats have led experts in the field of transportation to find alternatives for hot mix asphalts (HMAs). One of these alternatives is low-energy asphalt (LEA) mixtures. LEA mixture consists of hot asphalt binder and coarse aggregates and wet sand, which, having a lower mixing temperature, has caused a significant reduction in energy consumption and pollutants emitted from asphalt binder. According to previous studies, the influence of crumb rubber (CR) on the rutting and moisture performances of LEAs has not been investigated so far. In this regard, due to the softness of this type of asphalt mixture and its weakness against rutting and moisture, in this study, as an innovation, CR was used to modify LEA mixtures in three weight percentages of 10, 15, and 20. Considering the high moisture sensitivity of LEA mixtures due to the presence of water, the ability of CR to overcome this weakness and its potential to be replaced with antistriping additives were also investigated using the indirect tensile strength (ITS) test. In addition, in order to evaluate the rutting resistance, LEA mixtures were subjected to the dynamic creep test at a temperature of 54.4°C and at two stress levels of 270 and 310 kPa. According to ITS test results, it was found that although the modified LEA mixtures with 10 and 15% CR met the minimum requirement, the tensile strength ratio (TSR) value recorded for mentioned mixtures was still lower than that for HMA. Therefore, it is recommended to use antistripping additives in the use of LEA-containing CR, especially in areas with a high probability of moisture damage. The results of creep curves and flow numbers (FNs) extracted from four different models of three-stage, Fnest, Francken, and stepwise showed that although the LEA sample without additive had lower rutting resistance than HMA, the modification of this mixture to 15% CR significantly increased the rutting resistance. Moreover, considering the variety of variables in the process of making LEA compared to HMA, the comparison of FN determination models and the introduction of the most stable model were conducted in this research. According to the laboratory and statistical analysis results, it was found that the Francken model had less sensitivity to the variables, so it is suggested to determine FN for LEA mixtures.

1. Introduction

Today, with the advancement of technology, negative effects, such as environmental pollution, are the most important threat to human societies. The emission of greenhouse gases has resulted in global warming and climate change [13]. One of these polluting factors is the production of hot mix asphalt (HMA), which has a high consumption of fossil fuels and releases many pollutants due to the high production temperature (150–180°C). Increasing concerns about the preservation and protection of the environment led many researchers to find alternatives for HMAs [4, 5]. One of these efforts is to achieve warm mix asphalt (WMA), which has led to the production of asphalt mixture at a temperature of 110–140°C by using cleaner technologies [6]. Recently, with more attention to sustainability issues in the road pavement industry, there have been innovations to produce mixtures at a temperature of 60–100°C. This asphalt mixture is called half-warm mix asphalt, which has presented many environmental benefits with acceptable performance [7].

LEA mixtures are part of half-warm mix asphalts, which were produced at a range of 80 to 100°C by LEA-Co in France in 2006 [8]. This asphalt mixture is produced in such a way that after separating the aggregates into two sections, coarse and fine, first the coarse aggregates are heated to 170°C, and then, they are mixed with hot asphalt binder. After ensuring proper asphalt coating, fine aggregates are added at ambient temperature with desired humidity. The contact of moisture with a hot mixture produces foam, increases the asphalt binder volume, and results in a slow decrease in viscosity. Among the benefits of production, this asphalt mixture is a significant saving in energy consumption. Since fine aggregates are added wet to the mixture, there is no need to dry them, and therefore, fuel consumption is lower in the production process. Also, low production temperature results in less emission of greenhouse gases from asphalt binder. In addition, the reduction of asphalt binder viscosity due to the foam formation and its less aging due to the low production temperature causes the mixture to soften, and as a result, the mixture has better compactability. Greater safety during execution and less consumption of solvent to clean the equipment at the end of daily operations are other benefits of using LEA [9].

Despite the above advantages, since the moisture sensitivity of asphalt mixtures affects their performance, the investigation of moisture damage is a fundamental challenge in their use. Moisture damage occurs due to weak adhesion and lack of cohesion in asphalt mixtures. Half-warm asphalt mixtures and especially technologies based on asphalt binder foam have less adhesion between asphalt binder and aggregates [10]. Also, less asphalt binder coating on aggregates due to the low viscosity of asphalt binder causes lower integrity of these mixtures. Studies conducted on LEA mixtures show that moisture damage is one of the weaknesses of this asphalt mixture. In this regard, in the evaluation of moisture resistance by Zelelew et al., by determining indicators, such as Humburg stripping slope and Humburg inflection point, it was indicated that LEA mixtures have a lower moisture resistance than HMA mixtures [11]. Also, the results of the moisture resistance test in the research conducted by Olard et al. showed that the LEA sample has a lower moisture resistance than the HMA sample. In this evaluation, the r/R ratio (in which r is the compressive strength of samples stored in water and R is the compressive strength of samples stored in the air) in the Duriez test (NF P 98 251-4) was reported to be almost 10% smaller than the ratio obtained for HMA sample [12]. In addition, in the research carried out by Some et al., the evaluation of the moisture sensitivity of asphalt mixtures according to EN 12697-12 indicated a lower moisture resistance for LEA than HMA [13]. Therefore, due to the low moisture resistance of LEA and due to the presence of water in its production process, the use of an additive capable of increasing adhesion and reducing moisture sensitivity is emphasized to compensate for this deficiency.

Also, the other important weakness of half-warm asphalt mixtures, including LEA, is the high sensitivity of these asphalt mixtures against rutting. Due to the low mixing and compaction temperatures, the asphalt binder faces less aging, so this mixture is softer compared to HMA and has less rutting resistance [14]. Rutting or permanent deformation is the accumulation of permanent strains at the place of vehicle wheel loading, which leads to a decrease in pavement performance and service quality [15]. Despite extensive research on the rutting resistance of half-warm mix asphalt, little research has been performed on the susceptibility of LEA mixtures against this failure. Among the studies conducted in this regard, we can note a study presented in the NCHRP report 691, in which the dynamic creep test was used in order to evaluate the rutting resistance, and four asphalt samples including LEA, HMA, a sample modified with Advera, and a sample modified with Sasobit were tested. According to the findings, the highest value of flow number (FN) was assigned to the sample modified with Sasobit, and the lowest value of this parameter was recorded for the LEA sample. In this evaluation, the HMA sample showed much higher resistance to rutting by recording FN about twice the LEA sample. Also, the sample containing Advera was in a better condition than LEA with a slight difference in FN value. In the other part of this research, the rutting resistance of LEA was evaluated by conducting the dynamic modulus test and comparing it with HMA. The master curve drawn from the test results showed that LEA is significantly softer than HMA. In this research, by increasing one level in the high-temperature grade LEA sample, it showed better rutting resistance, so that the rutting depth for this mixture reached 0.05 inches with a significant reduction [16]. Zelelew et al. studied the rutting resistance of LEA by conducting a wheel track test using the Hamburg method. For this purpose, the number of wheel passes regarding 20 mm rut depth was evaluated. According to the results, HMA has the highest rutting depth followed by warm mix asphalt (WMA) samples containing Sasobit and Advera. In this evaluation, the sample containing Gencor and LEA sample, respectively, showed the weakest rutting resistance [11].

In the process of urbanization and construction [17, 18], a large number of waste materials are produced, one of which is waste tires in the field of automobile manufacturing [19]. Urban infrastructure is important to promote sustainable development [2022]. These waste rubber materials are difficult to degrade naturally, and incineration and landfill will do some harm to the environment and also bring financial losses to societies. Crumb rubber (CR) is one of the most significant asphalt modifiers utilized in several years in the pavement industry, which is reclaimed from the tires of vehicles. CR application in asphalt mixtures has been examined for several years as a suitable additive for enhancing the properties of unmodified mixtures and crumb tire removals. Previous research showed that the use of CR as a thermoplastic elastomer in virgin binders enhances the characteristics of asphalt mixtures to rutting, moisture, and fatigue, decreases the maintenance and rehabilitation costs of mixtures, improves pavement life, and reduces traffic noises [23]. By mixing CR and asphalt binder, CR is immersed in the asphalt binder and quickly undergoes a physicochemical reaction. This two-stage reaction includes swelling and degradation, which is a function of mixing time and temperature. At first, with increasing temperature, aromatic oils are absorbed by polymer chains of rubber and physically cause swelling of rubber particles. By increasing the size of the rubber particles and decreasing the oil part of the asphalt binder, the interparticle distance decreases. This causes the formation of a gel-like structure in the asphalt binder and increases its viscosity greatly. By increasing the temperature and mixing time, the swelling reaches the point where the rubber undergoes degradation. This part of the reaction has a chemical nature and includes devulcanization and depolymerization, which occurs by increasing the temperature up to 260°C after a few hours and a shear of 8000 rpm [24]. On the other hand, the studies conducted regarding the use of CR to modify asphalt binder show that the complete reaction of CR in asphalt binder has a negative effect on its characteristics. Also, raising the temperature leads to the destruction of asphalt binder networking. Therefore, the greatest effect of CR in asphalt binder modification occurs when devulcanization and depolymerization do not take place completely, but the asphalt binder is sufficiently swollen and the reaction is carried out at a partial level of chemical reaction [25].

Moreover, regarding other characteristics of LEA mixtures, Carter et al. examined the compatibility of LEA, HMA, and WMA samples made with Sasobit and Evotherm. In this research, LEA and samples modified with Sasobit were the most compatible samples, followed by HMA and samples made with Evotherm [26]. Olard and Gaudefroy compared the performance and environmental characteristics of LEA and HMA by determining the complex modulus at 15°C and 10 Hz according to European Standard NF EN 12697-26, and it was found that LEA has a lower complex modulus than HMA. This result was presented due to the low mixing and compaction temperatures, as a result of less aging of the asphalt binder and the softness of the mixture. In addition, performing a fatigue test showed lower fatigue resistance of LEA compared to HMA. Another part of this research was the environmental evaluation of the samples. In this evaluation, the published total organic compounds value for LEA was reported to be half of the value recorded for HMA [12]. In another study, Gregory and Harder achieved a significant reduction of pollutants in the LEA production process. According to the results, it was found that LEA production is associated with a 51% reduction in suspended particles, 21% nitrogen oxide, and 82% carbon monoxide. Also, the amount of sulfur oxide and formaldehyde emissions was reported to decrease by 46% and 95%, respectively. Another result of this research was the reduction of fuel consumption so that the energy consumption for the production of LEA showed a 33% reduction compared to HMA [27].

2. Research Objectives

The purpose of this research is to evaluate the rutting and moisture resistance of LEA mixtures. Considering that the effect of additives in LEA modification is the research gap of previous studies, the use of CR to increase the performance characteristics of LEA is considered an innovation of this research. In this regard, the ability of CR to overcome moisture sensitivity and its potential to be replaced with adhesion-enhancing additives is investigated as the first objective of this study by performing the indirect tensile strength test. Also, in order to examine the rutting resistance, the dynamic creep test is used, and by determining FN extracted from four conventional methods, the rutting resistance of asphalt mixtures is compared with each other, and the effect of CR in increasing the rutting resistance of LEA is investigated as the second objective of this research. Moreover, due to the variety of variables in the LEA construction process, the comparison of FN determination models and the introduction of the most stable model are other objectives of this study. In this regard, the trend of FN changes obtained from different models is examined, and finally, the dependence of the models on the mentioned variables is evaluated with statistical analysis.

3. Methodology

3.1. Materials

In this study, limestone aggregate was used to produce LEA mixtures. Table 1 shows the characteristics of used aggregates. The gradation of aggregates, as illustrated in Figure 1, was according to ASTM standards the nominal and maximum sizes of the aggregates are 12.5 mm and 19 mm, respectively. Also, a base asphalt binder with a 60/70 penetration grade prepared from the Jey Oil Refinery was selected. Table 2 shows the physical characteristics of the asphalt binder.

3.2. Additive

According to the studies conducted on CR and the advantages of the wet process in the production of rubberized asphalt mixtures, in this research, this method was used. Also, because of the considerable influence of crumb rubber on the properties of asphalt mixtures in the range of 10–20% by the weight of the binder, CR percentage was limited to this range [28]. Table 3 presents the physical and chemical properties of CR.

On the other hand, research that has been performed regarding the performance of different types of CR shows that CR-modified asphalt binders produced by ambient grinding method with a smaller size have better performance, such as stability, viscosity, stiffness, and moisture and rutting resistance [29]. Therefore, in this research, CR produced by the mentioned method was used with 40 mesh size, and Table 4 represents the gradation of used CR based on ASTM D 5644.

In addition, the asphalt binder temperature during mixing was selected at 180°C, and a high shear mixer with 2000 rpm was used for 2 hours to achieve a homogeneous mixture. Table 5 presents the characteristics of asphalt binders modified with CR in the mentioned percentages.

In this research, two methods of a softening point difference in aluminum tubes and morphology evaluation of the mixture were used to ensure the compatibility of the mixture. In the first method, according to ASTM D7173, modified binders were heated to a temperature of 163°C and poured into aluminum tubes and kept in an oven for 48 hours. In the next step, they were placed at a temperature of −18°C for 4 hours. Finally, each sample was cut into three equal parts, and the softening point was determined from their upper and lower parts. As can be seen from Table 6, although the softening point difference is increasing with the increase of CR percentage, the values obtained are lower than 2.5°C. Therefore, it can be assured that the obtained mixtures have good homogeneity and acceptable storage stability [30]. Another common method to evaluate the compatibility of polymeric asphalt binders is to examine their morphology. For this purpose, preparing scanning electron microscopy (SEM) images is an effective tool, so showing the size, increasing the volume of polymer particles, and how they are distributed inside the asphalt binder are considered advantages of their application. Figure 2 represents an SEM image of the asphalt binder sample modified with 15% CR. As can be seen, CR particles are distributed evenly inside the binder, so CR particles have a homogeneous distribution inside the binder, and no agglomeration is observed in them. In addition, the close distance of CR particles to each other indicates the fact that these particles are sufficiently swollen so that their interparticle distance is decreased. Therefore, it can be expected that CR and asphalt binder have a homogeneous mixture, so they will have good compatibility and storage stability.

3.3. Mix Design

In this research, the Superpave volumetric mix design based on the AASHTO R35 standard was used to design the samples. Also, according to the studies conducted on LEA and limiting the maximum size of fine aggregates between 2 and 3 mm, sieve No. 8 was selected for the separation of fine aggregates and coarse aggregates [11, 12, 27].

Considering the necessity of mixing temperature in the range of 80 to 100°C, the choice of factors, such as the coarse aggregate temperature, asphalt binder temperature, and moisture content of fine aggregates, is decisive in LEA construction. According to NCHRP report 691, the moisture content of fine aggregates is recommended between 3 and 4%. Also, the coarse aggregates should be placed at 170°C for 4 hours, and fine aggregates should be used at room temperature, and the asphalt binder is utilized at 135°C [16]. Therefore, in this study, according to NCHRP report 691, firstly, the aggregates were divided into two categories, fine and coarse, by sieve No. 8. Then, enough water was added to the fine fraction, which included 44% of the weight of aggregates, and the final moisture content reached 3.5%. Simultaneously, coarse aggregates and asphalt binder were heated to 170 and 135°C, respectively. Afterward, the coarse aggregates were mixed with the asphalt binder in a wire mixer. After ensuring the complete coating of the obtained mixture, fine aggregates were added to the mixture. Sudden contact of moisture with the hot mixture caused the formation of foam according to Figure 3. At this stage, the mixing of the mixture was continued until fine aggregates were completely coated. After the mixing operation, the desired samples were compacted according to the AASHTO T312 standard, and their characteristics were determined, as indicated in Table 7, in which is the maximum theoretical specific gravity, is the bulk specific gravity, is the air voids, VMA is the voids in mineral aggregates, VFA is the voids filled with asphalt, and OBC is the optimum binder content.

3.4. Experimental Program
3.4.1. Moisture Sensitivity Test

The most common method for the evaluation of moisture sensitivity is the modified Lottman test, which is performed according to AASHTO T283. In this regard, six compacted cylindrical samples are made with air voids of 6–8%. Then, the samples are divided into two categories. The first category is considered unconditioned or control samples, and the second one, as conditioned samples, after obtaining 70–80% saturation and being placed at −18°C for 16 hours, is kept in a water bath at 60°C temperature for 24 hours. Finally, by performing the indirect tensile strength (ITS) test with a loading rate of 50 mm/min at a temperature of 25°C and recording ITS values, the tensile strength ratio (TSR) is determined through equation (1), where and (kPa) are the indirect tensile strength for the conditioned and dry samples, respectively [31].

3.4.2. Dynamic Creep Test

In this research, the uniaxial dynamic creep test was performed according to NCHRP report 465 by applying a repetitive shear haversine loading by a universal testing machine (UTM). First, the test samples were made in 7% air void and a cylindrical mold according to AASHTO T312, and after coring, samples with a diameter of 100 mm and a height of 150 mm were extracted. Since permanent deformation occurs at high temperatures, in this research, the test temperature was selected as 54.4°C according to the recommendation of the aforementioned report. The samples were then transferred to the loading frame, and LVDTs were installed on it to record the deformations created during the test. First, a constant load equal to 5% of the final load was applied. Then, haversine loading was introduced to the samples with a period of 0.1 seconds of loading and 0.9 seconds of rest time, and the loading continued until the failure of the samples or reaching 10000 cycles. In this research, according to the recommendation of the aforementioned report, the uniaxial dynamic creep test was performed at the stress level of 207 kPa. Also, in order to evaluate the rutting susceptibility of LEA to the increase in stress level, all samples were also tested at the stress level of 310 kPa. Then, ascending diagrams consisting of accumulated strain against loading cycles were drawn for each of the samples. These diagrams, which consist of three different areas, are known as creep curves. In the first area, as the loading cycle increases, the strain increases at a decreasing rate. In the second area, the strain increases at a constant rate, and finally, in the third area, the strain increases at an increasing rate. The number of cycles related to the starting point of this area is known as flow number (FN), and it is used to evaluate the rutting resistance of asphalt mixtures [32]. Many models have been presented to determine the location of FN, the definition of four of which has been discussed in the following.

(1) Three-Stage Model. This model was presented to interpret the creep diagram by Zhou et al. According to the proposed model, the power law model according to equation (2) was used to interpret the first area of creep curves. Also, a simple linear model according to equation (3) was used to interpret the second area. On the other hand, due to the rapid increase in the permanent strain after the second area, a model according to equation (4) was presented to describe the third area [33].in which a, b, c, d, and f are materials constants, is the accumulated plastic strain at the tertiary stage beginning, is the accumulated plastic strain at the secondary area beginning, is the number of cycles at the tertiary stage beginning, is the number of cycles at the secondary stage beginning, is plastic strain, and N is the number of cycles.

(2) FNest Method. This method was presented by Archilla et al. [34]. First, the strain rate is measured using equation (5), where N is the cycle number, is the permanent strain at the N + 1 cycle and is the permanent strain at the N − 1 cycle.

Then, the five-point moving average period technique is performed to smooth the strain rate values obtained from equation (5). Also, the range of load cycles applied is defined by equation (6), where FNMAP5 is the minimum point of this strain rate, Nlb is the lower number of cycles, Nub is the highest number of cycles, C is 0.5, and Nmax is the maximum cycle number.

In this method, it is assumed that NUB ≤ Nmax. Hence, when NUB > Nmax, NUB is considered equal to Nmax. By using the obtained data and the above explanations, equation (7) is used for fitting the data, where α, β, and γ denote probability distribution parameters.

It is worth mentioning that these parameters are obtained from the solver tool in Excel. By determining the above parameters, FN is obtained.

(3) Francken Model. Billgiri et al. used a combined mathematical model for determining FN [35], which is derived from the Francken model and is defined as equation (8), where A, B, C, and D are regression coefficients, N is the loading cycle, and εp (N) is the permanent strain.

These parameters are determined by using a statistical package and regression technique. Then, the strain slope is obtained by the following equation:

Finally, FN is obtained with the second differentiation of εp (N). In this case, FN is at the point that the values obtained from equation (10) change their sign from negative to positive.

(4) Stepwise Method. Goh and You presented a simple and step-by-step method to identify the beginning of the third area by modifying the traditional method [36]. According to the mentioned method, in the first step, the permanent strain obtained from the test is smoothed by the reallocating technique. This operation is performed by assuming that the permanent strain either remains constant or increases during the loading cycles. In the second step, the strain rate is measured according to equation (11), and eventually, FN is obtained by identifying the location of the min amount of strain rate.where N is the loading cycle, and ε is the permanent strain.

4. Results and Discussion

4.1. Moisture Sensitivity

As can be seen from ITS values represented in Table 8, the LEA sample has a lower ITS value than the HMA sample in both dry and conditioned states. The reason for this is the presence of many pores due to the confinement of part of the water vapor created in the process of making LEA in the asphalt binder and, as a result, easier penetration of water to the surface of aggregates [37]. The modification of asphalt binder by CR increases the viscosity of the asphalt binder, so it results in a thicker coating of asphalt binder on aggregates and its greater adhesion, and consequently, an increase in the resistance of the mixture and its integrity. Therefore, with the addition of CR up to 15% of binder weight, this asphalt mixture experiences a significant increase in ITS value, so that the ITS value recorded for the LEA15 sample in dry and conditioned conditions was 1.52 and 1.77 times this parameter, respectively, compared to LEA sample. With the increase in the percentage of CR from 15 to 20%, due to the excessive increase in asphalt binder viscosity, and consequently, the decrease in the efficiency of the mixture, the LEA20 sample experienced a decrease in coverability, adhesion, and cohesion, and therefore, ITS value for it in both dry and wet conditions was significantly reduced compared to LEA15 sample.

Another part of the information presented in Table 8 is related to TSR values recorded for control and LEA samples. In this evaluation, as expected, the LEA sample with a TSR value of 0.62 showed high moisture sensitivity. The low moisture resistance of this sample is due to the low viscosity of the asphalt binder, low quality of aggregate coating, and as a result, easier penetration of water due to the presence of many pores left in the bituminous coating because of the creation of water vapor during mixing [38]. By adding CR to LEA samples, the moisture resistance increased significantly, so that the application of 10% additive resulted in a 27% increase in TSR value in the LEA10 sample compared to the LEA sample without additive. With the increase in the percentage of CR, a downward trend was observed in the TSR values of LEA samples, so that LEA15 and LEA20 samples with a TSR of 0.72 and 0.63, respectively, faced a 10% and 25% decrease in the value of this parameter compared to LEA10 sample. The decrease in moisture resistance in the mentioned samples was due to the excessive stiffness of the asphalt binder and its inability to properly cover the aggregates and more opportunity for moisture to penetrate their surface. The studies performed on the performance of CR-modified samples against moisture damage indicate that these types of asphalt mixtures, despite their high resistance to rutting, have low moisture resistance, and only a small content of CR can lead to a decrease in moisture sensitivity [15, 3941]. The point to consider in this evaluation is the inability of CR to compensate for the drop in the moisture resistance of LEA samples compared to HMA samples [39]. As can be seen from the results, although LEA10 and LEA15 samples met the lowest recommended value for this test, the highest TSR recorded in this evaluation was assigned to the HMA sample. Therefore, CR, despite its proven benefits, seems not to be a suitable alternative to antistripping additives for LEA, especially in areas with high rainfall and situations with a high probability of moisture damage. Therefore, antistripping additives should be applied in the production of CR-modified LEA mixtures in the mentioned situations.

4.2. Creep Curves

Figures 4(a) and 4(b) present the permanent strain diagrams according to the loading cycle obtained from the dynamic creep test at the stress level of 207 and 310 kPa, respectively. It is evident that LEA15 and LEA20 samples have a better performance against rutting by recording lower accumulated strain during loading cycles. LEA and LEA20 samples showed the weakest performance, respectively, and the HMA sample is placed between the two mentioned groups. By adding CR to the asphalt binder, the rubber polymer chains absorb the aromatic part of the asphalt binder, which creates a gel-like state in the asphalt binder and increases its elasticity and viscosity. These interactions strengthen the internal network of asphalt binder and cause its thicker coating on aggregates, increasing cohesion and, ultimately, higher rutting resistance of mixtures. As can be seen from the diagrams, by increasing the CR percentage, the LEA20 sample represented a drop in resistance to rutting. This is considered due to the excessive absorption of asphalt binder’s light oils by rubber particles and, as a result, excessive increase of its viscosity, the results of which were also observed in previous studies with increasing CR from 15 to 20% [23, 42]. Hence, as the viscosity increases, the workability of the mixture is reduced during the mixing time, and consequently, this results in a weaker coverage of aggregates and low cohesion and integrity of the mixture. Another noteworthy point in this evaluation is the low resistance of LEA samples against permanent deformation compared to other samples. The excessive weakness of this sample due to the low mixing temperature compared to the HMA sample is the lower viscosity and aging of the asphalt binder, which has led to the softening of the mixture and the reduction of stiffness.

4.3. Flow Number

Figures 5(a) and 5(b) represent FN values obtained from the three-stage, FNest, Francken, and stepwise models at two levels of stress 207 and 310 kPa, respectively. It is obvious from the results that, apart from the stepwise model, other methods were successful in determining FN in all test samples, and the mentioned method was only able to determine FN in the HMA sample, LEA20 sample at 207 kPa level of stress, and LEA, LEA10, and LEA15 samples at 310 kPa level of stress. According to the results obtained from the four mentioned models, the LEA sample has the lowest FN value in both stress levels as expected. The low rutting resistance of the LEA sample is due to the less aging of the asphalt binder and, as a result, its excessive softness [16]. By using CR in LEA samples and its physicochemical reaction with the asphalt binder, the viscosity of the asphalt binder is increased, and as a result, the coating on the aggregates is improved. This has caused an increase in the adhesion between the aggregates and has finally led to an increase in the rutting resistance of LEA samples. Also, with the start of the chemical reaction of the asphalt binder and CR, new cross-linking bonds have been formed between the asphalt binder and CR, which has created a three-dimensional reinforcing network in the structure of the mixture. Therefore, the use of CR in the asphalt binder leads to an increase in the integrity and elasticity of the mixture, so that by using 10% CR, the LEA10 sample has a significant increase in FN and has a higher FN value than the LEA sample and even HMA sample. By increasing the amount of CR to 15%, the LEA sample still experienced an upward trend in FN, the result of which can be clearly observed in the study of other researchers [23, 43]. In this research, the modification of asphalt binder to the mentioned amount has resulted in an increase in FN value by 2 times for the LEA15 sample compared to the HMA sample in all the mentioned models. By increasing the amount of CR to 20%, the LEA sample faced a significant decrease in the amount of FN, similar to previous research [44]. The drop in rutting resistance of the LEA20 sample is due to the high viscosity of asphalt binder, the low mixing temperature of the LEA production process, and the low ability of asphalt binder in coating aggregates, which has resulted in a decrease in FN value compared to HMA sample.

On the other hand, comparing FN results obtained from the mentioned methods shows that the stepwise model provides the highest FN value in all test samples. Also, in the situation where FN obtained from the three-stage model has higher values, FNest and Francken models provided the lowest value of this index with relatively similar FNs compared to other methods.

Figure 6 shows the changes in FN with respect to the increase in the stress level. It is necessary to mention that due to the inability of the stepwise model to provide FN for all samples, it was not possible to calculate FN changes for it. As can be seen from the results, with the increase in the stress level from 207 to 310 kPa, FN ratio changes have presented a different trend, which is also true in the other three models. Although LEA15 provides the highest FN and is the most resistant sample against rutting, LEA10 has the lowest FN ratio. Therefore, this sample has the lowest decrease in FN with the increase in the stress level and provides the lowest sensitivity to the increase in the stress level. By increasing the amount of CR to 15%, the effect of the additive can still be observed in reducing the sensitivity to the stress level increase, so that LEA15 has a lower value of this ratio than the remaining three samples. The important point in this evaluation is the almost identical process of the three-stage and Francken models in presenting the order of the results of the rutting susceptibility. Also, determining the amount of correlation between different models according to Figure 7 shows that the three-stage and Francken models have the highest correlation coefficient () value, so it can be concluded that these models have more correlation to each other. But the important point about the FN ratio extracted from the three-stage and Francken models is that the LEA20 sample, despite having higher FN values than the LEA sample, has shown more sensitivity to the stress level. This can be due to the negative effect of increasing the additive percentage and, as a result, the excessive absorption of the binder oil part by the rubber particles and the brittleness of this sample in the face of high stress levels.

Figures 812 show the location of FN in the creep curves obtained from the test samples. As can be seen from the results, in LEA, LEA10, and LEA15 samples at the stress level of 207 kPa and HMA and LEA20 samples at the stress level of 310 kPa, even though the creep curves have entered the third area after passing through the second area, the stepwise model was not able to determine FN. In addition, in HMA and LEA20 samples at 207 kPa level of stress and LEA10, LEA, and LEA15 samples at the stress level of 310 kPa, FN obtained by the mentioned method is located in the third area, and its location is contradictory to the definition of FN based on the junction of the second and third areas. Therefore, it seems that the stepwise model is not a reliable method to determine FN in this research. Also, by observing the location of FN obtained from other methods, this fact can be achieved although it seems that FN obtained from the three-stage model is closer to the third area in all test samples, and the results of the other two methods are generally located on the straight line; however, the judgment in this regard is difficult, and choosing the best model requires more investigation. Therefore, in the next section, statistical analysis has been conducted to determine the best model and compare FN values extracted from the mentioned models.

4.4. Statistical Analysis
4.4.1. Investigating the Methods of Determining FN

For determining the appropriateness and correlation of results obtained from various models, the ANOVA test was used according to Tables 9 and 10 in two levels of stress 207 and 310 kPa, respectively, and the mean difference of FN obtained from different models was compared.

According to the comparison, the lowest mean difference in both stress levels was assigned to Fnest and Francken models. Therefore, the mentioned models showed the highest correlation in FN values, the result that was also observed in the research conducted by Ameri et al. by achieving similar FN values compared to other models [45]. Another noteworthy point in this comparison is the high correlation between the results of the stepwise model and the three-stage model. According to the obtained results, the mentioned models provide the lowest mean difference, and therefore, they have the best fit compared to others. Studies conducted by other researchers such as Goh et al. also represent the high correlation of the mentioned models by calculating value [36, 46].

4.4.2. Between-Sample Variability

A between-sample variability is an approach that is used to account for the variability of results for similar replicates. This solution is used when only one method is used to calculate the results [45]. Table 11 and Figure 13 present the coefficient of variation (CV) of FN for asphalt samples in each of the three methods of three-stage, FNest, and Francken. It should be mentioned that in this analysis, the stepwise model is removed due to the inability to determine FN for most of the samples and providing insufficient data. In this research, due to the existence of many variables in the construction of asphalt mixture samples, such as mixing and compaction temperatures, humidity, additive percentage, air voids, and asphalt binder percentage, as well as changes in the stress level, between-sample variability can be a suitable method to determine a model with the lowest sensitivity to the mentioned variables [45]. As can be seen from the results, at the stress level of 207 kPa, the three-stage model provides the highest CV in determining FN. Therefore, the mentioned method does not have enough stability and has the most changes in the results. At this stress level, the Francken model shows fewer changes than other models; therefore, the mentioned model has less dependence on the variety of asphalt samples and the changes in stress levels, and it can be considered the most stable model compared to others. Also, at the stress level of 310 kPa, the highest CV belongs to the Fnest model. At this stress level, although the three-stage model shows less dependence on the variables, the Francken model still has the lowest CV, the result of which was also obtained in the research conducted by Ameri et al. [45]. Therefore, as it is suggested in AASHTO T378, the use of the Francken model is recommended to determine the FN of LEA mixtures since it is less sensitive to variables.

5. Conclusion

In this research, by modifying LEA with CR in three weight percentages of 10, 15, and 20, evaluating the moisture and rutting resistance of LEA samples was conducted. Also, the comparison of FN determination models and the introduction of the most stable model were performed. According to the obtained results, it was found that:(i)In the moisture resistance results, although LEA10 and LEA15 samples met the minimum requirement in the ITS test, the TSR value recorded for these samples was still lower than that for the HMA sample. Therefore, it is recommended to utilize antistripping additives when using LEA-containing CR, especially in areas with a high probability of moisture damage.(ii)Investigating creep curves obtained from the dynamic creep test at 54.4°C and two stress levels of 207 and 310 kPa showed that LEA without additive had the highest amount of accumulated strain during the loading cycle compared to other samples. Also, by adding CR, a significant decrease in the accumulated strain of LEA samples was observed, so that the creep curves recorded for LEA10 and LEA15 samples were in a much better condition than HMA. By increasing the amount of CR to 20%, the accumulated strain experienced a significant increase.(iii)FN values obtained by different methods at two levels of 207 and 310 kPa showed that LEA15 and LEA10 samples had the highest FN value, respectively. In this evaluation, HMA was in a better condition than LEA20, and LEA was recognized as the weakest sample with the lowest FN value.(iv)The evaluation of the rutting susceptibility of laboratory samples to an increase in the stress level showed that LEA10 and LEA20 samples, respectively, provide the lowest FN changes with respect to the stress level increase. On the other hand, considering the variety of variables in the construction process of LEA compared to HMA, the comparison of FN determination models showed that the Francken model had less dependence on the variety of LEA samples and stress levels than other models, and it is recommended as the most stable model to determine FN for LEA samples.(v)For future studies, various nanomaterials and fibers can be applied in this regard [47, 48]. In addition, other failures of LEA mixtures containing CR and nanomaterials can be investigated [49, 50]. The use of electric vehicles and clean fuels can also help reduce the carbon footprint in the pavement industry [5153]. In addition, other statistical analyses can be incorporated with the proposed approaches [5456]. Various modeling methods and algorithms are recommended for further investigation in the continuation of this research [5759]. Different validation methods are recommended in this regard [60, 61]. Optimization algorithms can also be used to obtain the optimal additive content [6264]. Moreover, deep learning techniques can be applied to present more accurate results [65, 66]. Other types of asphalt binders with different penetration grades are also recommended to be used in a study to compare their results.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

In this study, Iranian governmental organizations have not been partners and sponsors, and this study is purely studious.

Conflicts of Interest

The authors declare that they have no conflicts of interest.