Abstract

Limited supplies of natural aggregates for highway construction, in addition to increasing processing costs, time, and environmental concerns, have led to the use of various reclaimed/recycled materials. Reclaimed asphalt pavement (RAP) and recycled concrete aggregate (RCA) have prospective uses in substantial amounts in base and subbase layers of flexible pavement in order to overcome the increasing issue of a shortage of natural aggregates. This research presents the development of an empirical model for the estimation of resilient modulus value (MR) on the basis of CBR values using experimental results obtained for 52 remoulded granular samples containing natural aggregates, RCA, and RAP samples. Statistical analysis of the suggested model shows promising results in terms of its strength and significance when t-test was applied. Additionally, experimental results also show that MR value increases in conjunction with an increase in RAP contents, while the trend for the CBR value is the opposite. Statistical analysis of simulation results using PerRoad and KenPave demonstrates that addition of RAP contents in the subbase layer of flexible pavements significantly improves its performance when considering resistance against rutting and fatigue. However, results of repeated load triaxial tests show that residual accumulative strain under a certain range of loading conditions increases substantially due to the addition of RAP materials, which may be disadvantageous to the serviceable life of the whole pavement structure.

1. Introduction

The highway construction industry is responsible for almost 30% of global air pollution and greenhouse gas emissions and contributes roughly a quarter of the total fossil fuel consumption across the world [1]. Replacing natural or virgin aggregates with high-quality recycled materials has considerable potential to reduce the carbon footprint of the road/pavement construction industry. The overall financial and environmental savings due to the replacement of natural aggregates with recycled materials can rationalise the stabilisation cost in pavement applications. Hence, low-carbon and low-cost substitutes for conventional aggregates are actively being sought by researchers worldwide [2]. Reclaimed asphalt pavement (RAP) [3, 4] and recycled concrete aggregate (RCA) [57] have been reported as the most commonly recycled materials used in different layers of flexible pavements, while use of recycled bricks [8], recycled glass [911], and fly ash [12] has also been documented by many researchers and institutions.

In general, RAP materials are a blend of coarse and fine aggregates and bitumen obtained from aged or expired asphalt pavements. Throughout the world, major use of RAP material in a surfacing layer is commonly termed as hot mix asphalt (HMA). RAP has been used in the surface layer of flexible pavements in combination with natural aggregates in different percentages, extending up to 80% in some cases [13]. Most of the researchers suggest a typical range of 20–50% [12, 14, 15]. This shows that RAP materials should also be used in base/subbase layers of pavement in addition to using them to make blends with natural granular material in HMA applications. RCA can be obtained from the revamping or demolition of different types of structures, including commercial and residential buildings or other civil engineering structures [1618]. Many researchers have found that the engineering properties of RCA are either comparable to, or even better than, typical natural aggregates used in road construction [19, 20]. This implies that RCA products have potentially suitable applications in subbase or base course. It should also be noted that RCA can be used as premium base course materials [21].

Conventional pavement design usually depends on the “California bearing ratio” (CBR) of the soil/aggregate used in pavement structure, while the resilient modulus value (MR) of unbound aggregates is the fundamental input parameter required in the mechanistic empirical design/analysis of pavement structures. The most reliable, and hence most desirable, way to determine the resilient moduli is through repeated triaxial load testing. However, because of the difficulties encountered with the test procedure, including time consumption and the economy of the project, other laboratory tests, such as CBR, would also be considered if a reliable correlation could be established for the estimation of the MR value. In the existing literature, there are many studies to correlate resilient modulus values to those with CBR values determined for the natural granular materials. However, studies incorporating RAP/RCA instead of natural materials are very limited. Furthermore, a rational or coherent mechanism of the correlation development could probably not be identified, owing to the fact that the mechanics of both of these tests are starkly different from each other.

The specific objectives of this paper are as follows:(1)To develop a rational model for the prediction of a resilient modulus value of unbound granular materials containing RAP/RCA as a major component on the basis of CBR values;(2)To evaluate the performance of blended samples used in subbase layers through the computer software packages PerRoad and KenPave;(3)To evaluate the long-term performance of blended samples under a range of cyclic loading conditions.

Furthermore, a brief literature review on correlations for the estimation of MR values, with major emphasis on the development of a correlation between CBR value and resilient modulus, is presented in the next section of this research paper.

2. Existing Regression Models for the Prediction of MR Value

From the existing literature, regression models for the prediction of MR value for granular material can be categorically divided into four types:

Class I. In this category, the models are based on single strength or stiffness parameter of soil such as(i)CBR values [2226];(ii)R-value [2731];(iii)Unconfined compressive strength [32];(iv)Undrained compressive strength [33].

Class II. In this category, the models are based on soil properties and stress state, for instance,(i)Bulk stress and index properties of the soil [34];(ii)Unconfined compressive strength and index properties of the soil [35];(iii)R-value and index properties of the soil [36].

Regression models based on this methodology have markedly varying intricacy and acceptability in the research community [32, 3739].

Class III. In this approach, resilient modulus value for a certain soil is obtained by considering a certain type of stress invariant or set of stress invariants, for instance,(i)Bulk stress [27];(ii)Confining pressure and deviator stress [40]; and bulk stress and atmospheric pressure [41];(iii)Bulk stress, atmospheric pressure, and deviator stress [42];(iv)Bulk stress, atmospheric pressure, and octahedral shear stress [43];(v)Atmospheric pressure, octahedral normal stress, and octahedral shear stress [44].

In this category, the model parameters are given simple numerical values.

Class IV. There are also certain constitutive equations for the estimation of resilient modulus values derived from considering soil’s physical properties incorporated in model parameters in addition to stress invariant [4548].

Since this research is focused on the establishment of a correlation between CBR value and resilient modulus value, it is logical to further explain such attempts in a historical perspective.

A number of researchers, including Porter [49, 50], Hight and Stevens [51], and Fleming and Rogers [52] pointed out that the CBR tends to be a bearing value (more of a parameter in terms of strength) rather than a support value (in terms of recoverable behaviour) of materials. Thompson and Robnett [53] could not find a suitable correlation between CBR and MR; Hight and Stevens [51] stated that the CBR does not correlate consistently with either strength or stiffness; and Sukumaran et al. [54] opined that there is an apparent wide variation in the MR value that can be obtained using the CBR, which depends on many factors. On the other hand, Lister and Powell [55] felt positive that the CBR can be related, within reasonable limits, to subgrade stiffness. Hossain [56] believed that the CBR test is still one of the most widely used tests for evaluating the competency of pavement subgrade; however, there are variations in the procedure followed by different agencies (e.g., in terms of size of mould, compaction techniques, and efforts), and it was found that [56] correlations between resilient modulus values and all test results (including the CBR) were not statistically significant. Garg et al. [57] opined that the CBR value can be converted to a resilient modulus; a strong trend is apparent in the correlation but there is a lot of scatter. A few of the existing correlations based on CBR values are given in Table 1.

3. Material Characterisation

For this research, four different types of RAP samples (designated as “RAP(1),” “RAP(2),” “RAP(3),” and “RAP(4)”), three different types of natural aggregates (designated as “A,” “F,” and “W”), and one RCA sample were used. Based on the gradation curves, natural aggregate A is finer than other natural aggregate (F and W) while natural aggregate W is the most coarser among them. Each of the natural aggregate samples, as well as the RAP and RCA materials, contained crushed limestone of subangular to angular shape. Elongated and flat particles in the samples/materials were not more than 6% as per ASTM D 4791.

A summary of material characterisations in terms of detailed gradation properties (D10/D30/D50 grain size diameter corresponding to 10/30/50 percent finer by mass; Cc coefficient of curvature; Cu coefficient of uniformity), compaction characteristics, effective shear strength parameters, and physical properties of recovered bitumen binders has been presented in Table 2. Further detail on material characterisation can be found in Arshad and Ahmed [60] and Arshad [61]. Since focus of the testing campaign consists of the CBR test and the resilient modulus test, the same has been briefly explained in the following subsections.

Table 3 shows the matrix of the testing program including the resilient modulus tests, the CBR tests, and the repeated load tests.

3.1. CBR Test

In conventional pavement design, the CBR value is an important parameter used to determine the thickness of various layers of the pavement structure. Usually, the higher the CBR value, the better the performance of the pavement is, with regard to both stiffness and strength. This implies that the CBR value can be used as a parameter to evaluate the suitability of a soil for use as pavement construction material. For this study, standard 3-point CBR tests were performed on the natural/RAPs and blended samples under unsoaked conditions, to simulate the moisture content at which the resilient modulus tests were performed. A schematic diagram of the CBR test is shown in Figure 1.

The test equipment primarily consisted of a:(1)cylindrical mould having an inner diameter of 150 mm and a height of 175 mm;(2)spacer disc of 148 mm in diameter and 47.7 mm in height;(3)special surcharge weights;(4)metallic penetration piston of 50 mm diameter and minimum of 100 mm in length;(5)loading machine with a capacity of at least 50 kN and equipped with a movable head or base that travels at a uniform rate of 1.25 mm/min.

The CBR value is defined as the ratio of stress required for the circular piston to penetrate, at the rate of 1.25 mm/min, the soil mass in the cylinder to the standard stress that is required for the corresponding penetration of a standard material, i.e., like limestone found in California. Further procedural details on the CBR value calculations can be found in AASHTO T 193.

3.2. Resilient Modulus Test

Resilient modulus (MR) is defined as the ratio of cyclic axial stress to recoverable axial strain (). The cylindrical test specimen is compacted at a desired density and is subjected to cyclic axial stress at a given confining pressure within a conventional triaxial cell. Resilient modulus tests for this research were conducted as per the guidelines specified in AASHTO T 307-99(2004), for which a haversine-shaped loading waveform is mandatory to simulate traffic loading. Each load cycle of this waveform essentially consists of 0.1 seconds load duration and a 0.9 second rest period. Figure 2 shows a typical repetitive load/stress pulse along with the generated residual deformation curve in the time domain. The regular loading sequence in AASHTO T 307-99 consists of 15 different stages, each one having 100 load cycles, following the execution of the conditioning stage of 750 load cycles. Each of the loading stages has a particular combination of confining stress, maximum axial stress, and cyclic deviator stress as shown in Table 4. Further procedural details of the test including sample preparation and installation technique; electromechanics of the loading system (servo-controlled electrohydraulic MTS testing machine); specification of the used load cells and LVDTs; and particulars of the data acquisition system can be found in several works [6062].

4. Test Results and Discussion

This section presents a discussion on trends obtained for the CBR tests and the resilient modulus tests performed on a number of reconstituted granular samples as identified in Table 2.

4.1. Sample Results of the CBR Test

Sample results of the CBR test in terms of compaction effort and the percentage of RAP content on the CBR values are presented in Figure 3. From this figure, it can be inferred that the CBR value decreases noticeably in correspondence to the increasing content in the blended samples incorporating granular (A) and RAP(1). For instance, the blend containing 25% RAP(1) achieves only 74% of the CBR value achieved by the aggregate (A) (natural material) at the same level of compaction effort, i.e., 65 blows. Similarly, the corresponding value for blends containing 50% RAP(1), 75% RAP(1), and 100% RAP(1) is limited to 56.25%, 32.5%, and 22.5%, respectively. Likewise, the trend can also be observed for the other two compaction efforts consisting of 10 and 30 blows.

4.2. Effect of RAP Contents on Resilient Modulus (MR)

Forty-eight blended samples were prepared by mixing the four types of RAP materials with natural granular samples A, F, W, and RCA in proportion with RAP contents of 25%, 50%, and 75% for each. Figures 4(a)4(d) show the effects of the addition of these RAPs on the measured MR values over a range of bulk stresses as specified in AASHTO T 307-99 (68). Variation of MR values with bulk stress can be approximated through trend lines based on a power function having the coefficient of determination (R2) in the range of 0.95–0.99. From these trend lines, it is evident that MR values increase not only due to the increase in bulk stress but also in combination with the addition of RAP contents, i.e., blended samples give higher MR values when compared to those obtained from the natural samples of granular A, F, W, and RCA. More specifically, for instance, at a bulk stress of ∼673 kPa, MR value for the 100% granular W is 320 MPa while the corresponding values for blended samples incorporating 25%, 50%, and 75% of RAP(1) are 340 MPa, 360 MPa, and 390 MPa, respectively.

A similar trend in MR values was observed at lower levels of bulk stress, or more precisely, over the entire range of bulk stresses considered during the testing campaign. In some cases, ∼50% increase in MR values was observed for the blends containing 75% RAP contents. In general, for most of the blended samples containing 25% and 50% RAP contents, the corresponding increase in MR values was in the range of 5–15% and 10–20%, respectively. An important point is that the addition of RAPs to RCA induced the most significant increase in the resilient moduli when compared with the increase in MR when RAP was added to granular samples. Similar observations have also been documented by many researchers, including Kim and Labuz [12], MacGregor et al. [63], Alam et al. [64], and Bennert and Maher [65].

5. Development of Correlation between Resilient Modulus and CBR Values

5.1. Proposed Correlation and Theoretical Background

For this study, an attempt has been made to correlate the MR values with the CBR values on the basis of a common value of bulk stress identified during both types of tests. This should be emphasize that CBR values tend to decrease with addition of RAP contents while reverse is the situation in case of resilient modulus values. Such trends are primarily due to the fact that loading for the resilient modulus test is dynamic in nature while in the case of the CBR test, it is virtually static.

Fundamentally, axial stress is applied during the CBR tests through a plunger (σpa) in addition to an axial surcharge weight (σsa) as shown in Figure 5. However, lateral stress on the walls of the CBR mould can be estimated on the basis of the lateral earth pressure coefficient (Ka) as described in classical soil mechanics, such thatwhere is the effective angle of internal friction of the soil sample (blend) under consideration.

It is interesting to note that the 3-point CBR test encompasses a wide range of dry and bulk density in the sample, while on the other hand, each resilient modulus test uses a unique value of the sample’s dry and bulk density. To estimate a unique value of the bulk stress during the CBR test which is analogous to the bulk stress value identified from resilient modulus test, the following procedure was adopted:(1)Determine the dry density of the resilient modulus sample(2)Estimate the CBR value corresponding to that dry density value(3)Determine the axial stress () value corresponding to the CBR value(4)Calculate the bulk stress using the following relation:

Using a unique value of bulk stress on the basis of equation (2), the corresponding resilient modulus value from the actual experimental data (MR test) was matched and then a correlation between CBR value (point 2) and the resilient modulus value, in terms of the power function having a coefficient of determination ∼0.81, can be given as

Table 5 shows the experimental data and the estimated values of MR based on the four steps explained above and equation (3).

Figure 6 shows the capabilities of the proposed regression model on the basis of the distribution of the data points, which in fact involves the experimentally determined and estimated MR values along the unity line, giving a 95% prediction and confidence interval. This figure illustrates that there is an almost equal distribution of the data points on both sides of the unity line. Additionally, the entire set of the data points is confined within the 95% prediction interval, which is defined as the interval around the linear regression line such that 95% of the predicted values will fall in this interval. Further details on the mathematical framework of the prediction and the confidence interval can be found in standard textbooks dealing with statistics and probability analysis, such as those authored by Wonnacott and Ronald [66], Penny and Roberts [67], and Dybowski and Roberts [68].

5.2. Statistical Analysis of the Proposed Model

A Pearson’s correlation coefficient (r), as given by equation (4), is a measure of the strength and direction of linear association that exists between two continuous variables. More specifically, for the drawn line of best fit for sample data points of two variables, Pearson’s correlation indicates how well the data points fit this new model/line of best fit, and its numerical value indicates how far away all these data points are to the line of best fit (i.e., how well the data points fit this new model/line of best fit) [69, 70]:where and  = experimental and predicted values, respectively, for the ith output,  = average of experimental outputs, and n = size of sample.

The statistical value of r may range from −1 for a perfect negative linear relationship (inversely related) to +1 for a perfect positive linear relationship (directly related). In general, from theoretical point of view, the strength of the linear relationship can be categorised as “very strong,” “moderately strong,” and “fairly strong” for the corresponding numerical values of r in the range 0.8–1.0, 0.6–0.8, and 0.3–0.5 [71]. It should be emphasized that if r = 0.0 it does not necessarily mean that the two variables have no relationship. In order to look into the real strength of the correlation, usually, a statistical test of significance is performed, as discussed in [72]. For this study, three different types of statistical tests were applied to assess the significance of the developed correlation presented in equation (3) and Table 5.

Case 1. t-Test for the assessment of the implication of coefficient of correlation (r) at a specific degree of freedom and level of significance [73].
This provides the researcher with some idea of how large a correlation coefficient must be before it can be considered as demonstrating that there really is a relationship between two variables (in our case the (MR)measured and CBR values). It may be the situation that two variables are related by chance, and a hypothesis test for r allows the researcher to decide whether the observed r could have emerged by chance or not. The null hypothesis is that there is no relationship between the two variables. That is, if ρ is the true correlation coefficient for the two variables and when all population values have been observed, then the null hypothesis can be given asThe alternative hypothesis could be written aswhereas the standardised t-test for the null hypothesis that r is equal to zero can be written aswhere is the number of paired observations in the given sample.
The null hypothesis is evaluated by comparing the t-statistic of equation (7) with t-critical (2.01) obtained from t distribution having the degree of freedom.

Case 2. Direct comparison of the coefficient of correlation value (0.9) with the r-critical value (0.27) obtained from a statistical table corresponding to a specific degree of freedom and level of significance [61, 74].

Case 3. It often becomes mandatory to statistically evaluate the difference between two datasets obtained by two different sources. As in our case, one set of MR values was obtained experimentally using the AASHTO T 307-99(2004) protocol and the other set of MR values was obtained using the proposed correlation. A standard t-test was performed to evaluate the difference between the paired values of (MR)measured and (MR)estimated at a specific degree of freedom and level of significance [74, 75].
Table 6 summarises the statistical analysis of the above three cases.

6. Assessment of Pavement Performance Containing RAP Content in Its Subbase Layer

To examine the performance of the pavement containing natural aggregates mixed with RAP and used as subbase materials, the following analyses were performed using computer software simulation:(1)To access the likelihood that critical pavement responses exceed predefined thresholds using computer program PerRoad [76];(2)To determine the stresses and strains at critical locations in the pavement structure using computer program KenPave developed by the University of Kentucky [77].

PerRoad is a Monte Carlo-based simulation software used to develop probability-based analysis for flexible pavement. This software can easily demonstrate the influence of the MR values of the pavement material on the cumulative damage factor (CDF). It also demonstrates the likelihood that critical pavement responses could exceed predefined thresholds, which are the horizontal tensile strain of 70 microns at the bottom of the asphalt concrete (which is linked with fatigue cracking) and the vertical compressive strain of 200 microns at the top of the subgrade (which is associated with the structural rutting) [78].

6.1. Pavement Structure and Material Characteristics

Figure 7 shows a typical cross section of flexible pavements, which consists of a hot mix asphalt layer supported by the unbound base, unbound subbase, and compacted subgrade. The thickness of each layer generally depends on the traffic load or more specifically, the equivalent single axle load (ESAL) during the proposed life cycle of the road.

In this parametric study, the focus is placed on how the properties of granular subbase materials are changed by the addition of a certain amount of RAP and their effect on the pavement’s performance. As such, the resilient modulus of the granular base and the asphalt concrete is fixed. The structure of the pavement to be simulated and the properties of materials in different layers are summarised in Table 7. In this particular study, the subbase material will be one of the aggregates or aggregate/RAP blends tested during the experimental studies of this research.

The properties of the different pavement layers used in the analysis are shown in Table 7. The bulk stresses at the top and bottom of the granular subbase layer were determined using the computer program KenPave. For the selected HMA resilient moduli (MR = 5000 MPa), the bulk stresses were found in the range of 27 kPa to 30 kPa and 14 kPa to 18 kPa at the top and bottom of the subbase layers, respectively, for the wheel load of 40 kN. In the simulations, values of MR are used corresponding to the bulk stress of 100 kPa.

6.2. Traffic Load

For this parametric study, the traffic data for urban interstate highways, as recommended by AASHTO, were used. In the case of the PerRoad software, traffic is separated by axle type: single axle, tandem axle, tridem axle, and steer axle. After determining the percentage of each axle type in the total traffic, traffic is then subdivided into weight classes in 2 kip intervals. The following is the summary of the traffic data: the average annual daily traffic (AADT) is 1000 vehicles with 10% being trucks; annual growth rate of traffic is 4%; the directional distribution is 50%; and the percentages of single, tandem, and tridem axles are 55.73%, 42.66%, and 1.61%, respectively. The rest of the traffic loading characteristics, in terms of the distribution of vehicle types and axle weights, are described in Figure 8. These values were used as input in PerRoad 3.5. Axle weights were used to evaluate the response of the pavement layers in terms of stresses, strains, and deformations at critical locations of certain layers of the flexible pavement.

6.3. Pavement Performance Criteria

In the design of conventional flexible pavements, pavement sections are designed corresponding to a cumulative damage factor (CDF) of 1.0, which corresponds to a terminal level of pavement damage. For perpetual pavements, however, it is recommended that the CDF is equal to 0.1 at the end of its design life [76]. The fatigue and rutting algorithms developed at Mn/ROAD [79] for the calibration of flexible pavement performance equations were used to predict pavement damage in cases where pavement responses exceeded these thresholds. These correlations are as follows:where Nf = number of load repetitions when fatigue failure occurs, Nr = number of load repetition when rutting failure occurs,  = the horizontal tensile strain at the bottom of the HMA, and  = the vertical compressive strain at the top of the subgrade.

It should be noted that the damage of rutting estimated by PerRoad only takes into account the permanent deformation in the subgrade. Any rutting associated with the deformation of granular subbase materials is not considered. According to the test results in this study, use of RAP in granular subbase layers may induce substantial residual deformation. Further investigation should be performed to get a better understanding of its influence on the rutting of flexible pavements.

6.4. Simulation Results and Discussion

Simulation results are obtained and discussed in terms of(1)The likelihood of the tensile strain at the bottom of the HMA exceeding 70 με;(2)The likelihood of the vertical strain at the top of the subgrade soil exceeding 200 με;(3)The number of years it could take to reach a CDF of 0.1 (the threshold level) for fatigue damage;(4)The number of years it could take to reach a CDF of 0.1 for damage induced by rutting.

Figure 9 illustrates the concept frequency distribution for strain both within and outside the threshold limits.

We first examine the likelihood of critical pavement responses exceeding the predefined thresholds when the resilient modulus of the asphalt concrete is 5000 MPa, and different aggregates or aggregate-RAP blends are used in the granular subbase layer. Figure 10 summarises (1) the likelihood of the tensile strain at the bottom of the HMA remaining within the critical value of 70 με and (2) the likelihood of the vertical strain at the top of the subgrade soil exceeding 200 με. From this figure, it can be interpreted that all of the blends result in better pavement performance than the natural aggregate both in terms of the tensile strains at the bottom of the HMA and the vertical compressive strains at the top of the subgrade soil. For instance, the natural aggregates and RCA on average have an 81.3% likelihood that horizontal strain at the bottom of the HMA will remain within the critical limit of 70 με.

On the other hand, the likelihood for blended materials containing 25%, 50%, and 75% of RAP materials in the blended samples may reach (on average) 85.75%, 89.37%, and 93.77% (respectively) chance of remaining within the limit. Similarly, the likelihood that the compressive strain at the top of the subgrade will not exceed the critical value of 200 microns is 88.28% (natural samples), 93.46% (25% RAP material blends), 95.57% (50% RAP material blends), and 98.2% (75% RAP material blends).

Figure 11 shows the number of years required to reach a CDF of 0.1 when fatigue is controlled by the horizontal tensile strain at the bottom of the subgrade, and the HMA and vertical structural rutting are controlled by the vertical compressive strain at the top of the subgrade. The number of years required to reach a CDF of 0.1, in general, increases in line with the increase in the quantity of RAP in the blend. For example, the number of years required to reach a CDF of 0.1 in terms of fatigue damage at the bottom of HMA is 40.33 years for natural samples, while the corresponding figure for the blends containing 75% RAP is 51.68 years on average. A similar trend regarding an increase in the number of years required to reach a CDF of 0.1 in terms of structural damage at the top of the subgrade was also observed. Furthermore, for all the cases, the number of years to reach a CDF of 0.1 was more than 40 years, but it is clear that fatigue will be the predominant pavement failure mode.

7. Long-Term Effect of Cyclic Loading on Blended Samples

In order to explore the long-term effects of cyclic loading, 9 repeated load triaxial tests were performed under a range of cyclic loading conditions and varying percentages of RAP contents. Figure 12 presents the variation of the residual strain of the tested samples subjected to 20000 load cycles under two different values of confining pressures of 34.5 kPa and 137.9 kPa each, while the corresponding cyclic deviator stress was maintained at 31.05 kPa and 372 kPa. From this figure, it can be inferred that the presence of RAP contents increases the residual strain considerably for both of the confining pressure scenarios. For instance, for the repeated load test conducted at the confining pressure of 34.5 kPa, the residual strain for 100% natural aggregate F is 0.12% after 20000 load cycles, while for the blend containing 50% RAP(2) and 75% RAP(2), the corresponding figures reach 0.22% and 0.60%, i.e., an addition of 50% RAP increases the residual strain by a margin of almost 83% while the addition of 75% RAP increases the residual strain by almost 400%. Similarly, for the repeated load tests conducted at a confining pressure of 137.9 kPa on samples containing natural aggregate W and the blends with RAP(1) at 50% and 75%, the residual strain after 20000 load cycles is 55% and 300% higher when compared with the corresponding value for the 100% natural aggregate W. Furthermore, it is evident from the figures that almost 60% of the total residual strain occurred during the first 2000 load cycles out of the 20000 total load cycles applied across all the cases.

However, it should be emphasized that there is a tendency for the elastic shakedown to be less pronounced for blended samples when compared to pure natural aggregates. For instance, for the blended sample 50% W and 50% RAP(1), the rate of increase of strain becomes 0.000154% for the load cycles in the range of 2000–20000. On the other hand, this figure remained limited to 0.000056% for the sample containing 100% natural aggregate W under the same loading condition.

Figure 13 shows the effect of the ratio of “cyclic deviator stress to the confining stress ” on the residual strain for the blended sample containing 50% A and 50% RAP(3). From this figure, it can be inferred that the ratio has a substantial effect on the residual strain generated under cyclic loading. For instance, at , the residual strain after 6000 load cycles is limited to 0.06%; however, this value reaches 0.11% and 0.21% at and , respectively. The samples tend to stabilise more quickly under a lower value of when compared to the stabilising tendency at the higher values of and , which were also investigated in this research. A 50–60% residual strain occurs during the first 1000 load cycles out of the total 6000 load cycles applied during the experimental campaign, irrespective of the value.

8. Summary and Conclusions

(1)The blended samples (i.e., RAP combined with natural granular materials) result in higher MR values than those obtained for natural granular samples under the same loading conditions, while variations between MR values and the applied bulk stresses during the resilient modulus tests can be approximated through power law having coefficient of correlation (r) value in the range of 0.97–0.99.(2)The CBR values of blended samples decrease with the addition of RAP contents; however, a clear decreasing trend in conjunction with an increasing percentage of RAP contents could not be found.(3)The new model for the estimation of MR values includes the stress state through the experimentally determined CBR having a coefficient of correlation “r” equal to approximately 0.9, with fair distribution of the data points (MR measured and MR estimated) about the unity line indicating that the proposed regression model has a moderate to strong correlation.(4)Statistical analysis of simulation results based on PerRoad and KenPave software demonstrates that the addition of RAP contents in the subbase layer of the flexible pavement significantly improves its performance against rutting and fatigue.(5)Residual strain during the long-term repeated load triaxial test was found to increase in line with an increase in the ratio of “cyclic deviator stress to the confining stress ” in correlation to an increase in the percentage of RAP contents in the blended sample;(6)From the author’s perspective, the addition of RAP contents beyond a certain limit (∼50%) may prove to be detrimental for the overall performance of flexible pavement structure.

Data Availability

The experimental data used to support the findings of this study are included within the article in the form of graphs and tables.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.