Next Article in Journal
Can Fulfillment of Social Responsibility Enable Enterprises to Innovate? The Role of Corporate Financialization and Agency Costs
Next Article in Special Issue
Life Cycle Assessment of Natural Zeolite-Based Warm Mix Asphalt and Reclaimed Asphalt Pavement
Previous Article in Journal
Status and Trends of Membrane Technology for Wastewater Treatment: A Patent Analysis
Previous Article in Special Issue
Enhancing Effect of Waste Engine Oil Bottom Incorporation on the Performance of CR+SBS Modified Bitumen: A Sustainable and Environmentally-Friendly Solution for Wastes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soft Computing Approach for Predicting the Effects of Waste Rubber–Bitumen Interaction Phenomena on the Viscosity of Rubberized Bitumen

Department of Civil Infrastructure and Environmental Engineering, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
Sustainability 2022, 14(21), 13798; https://doi.org/10.3390/su142113798
Submission received: 8 August 2022 / Revised: 15 September 2022 / Accepted: 12 October 2022 / Published: 24 October 2022
(This article belongs to the Special Issue Research and Application of Recycling Asphalt Materials)

Abstract

:
The ability to anticipate the effects of the interaction between waste rubber particles from end-of-life tires and bitumen can encourage the use of rubberized bitumen, a material with proven environmental benefits, in civil engineering applications. In this study, a predictive model of rubberized bitumen viscosity is presented for this purpose. A machine learning-based approach (Multi-Gene Genetic Programming—MGGP) and a more traditional multi-variable least square regression (MLSR) method are compared. The statistical analysis indicates that the robustness and the capability of the MGGP algorithm led to a better estimation of the rubberized bitumen’s viscosity. Additionally, the MGGP analysis returned an actual equation that could be easily implemented in any spreadsheet for an initial tuning of the production protocol based on the desired level of interaction between the rubber and bitumen.

1. Introduction

The use of waste tire rubber in bitumen for road paving applications has been evaluated extensively by several research groups around the world from both mechanical performance [1,2,3,4,5,6,7] and environmental perspectives [8,9]. Since the discovery of this technology in the early 1960s, the use of crumb rubber (CR) as a bitumen modifier has experienced a significant evolution. However, its field implementation is still limited compared to the claimed benefits. One of the reasons is that the interaction phenomena between rubber and bitumen are highly variable and the control of the modified bitumen production is troublesome.
In general terms, the interaction starts as soon as the CR particles come into contact with the asphalt bitumen at high temperatures (150–200 C). Due to the diffusion of the oily fractions of the bitumen into the rubber network, a layer of gel-like material forms in the outer part of the particles. The extension of this layer extends from the rubber–bitumen contact surface to the maximum point of penetration of the oily fraction into the rubber and leads to the swelling of the particles [10]. The presence of the rubber cross-linked chains prevents the degradation of the particles for a limited time only. Then, the high temperature and mechanical stirring action initiate the devulcanization and depolymerization process of the rubber. Some components (e.g., carbon black) are now released into the bitumen matrix and the fragmentation of the rubber particles starts [10,11,12,13]. The magnitude of the above-mentioned phenomena is extremely variable and the production protocol plays a fundamental role. Rubber particles, in fact, react in a time-temperature dependent manner [14]. Swelling and degradation can speed up or slow down by either increasing or decreasing the blending time and temperature [15,16]. The physical and chemical characteristics of the rubber and the chemical composition of the base binder also play a crucial role in the interaction. In fact, the higher the surface area, the higher the quantity of oil that can be absorbed by the rubber, if all other variables are kept constant. Then, the maximum degree of swelling can be reached in a shorter time [14].
According to the IBISWorld Industry Report on Tire and Rubber Recycling released in September 2018, in the US only, there are more than 100 crumb-rubber production facilities [17]. A study conducted in Italy by Ecopneus, reported 20 manufacturers selling crumb rubber in a territory that is a hundred times smaller than the US [18]. Each of these factories produces more than one type of crumb rubber, which leads to many different CR available on the market. When the purchase decision in merely economic, the same bitumen refinery may use different types of crumb rubber in the same paving season.
Regardless of the type of materials used, all interaction processes have something in common. In the first stage, the simultaneous reduction in the inter-particle distance due to the swelling action leads to an overall increase in viscosity [11], whereas in the second stage, as rubber starts to degrade, the particles become smaller with a consequent reduction in viscosity. Figure 1 shows the mixing phases and their effect on the viscosity of the modified bitumen. As such, viscosity, even though indirectly, can be an indicator of the interaction phenomena.
Where an asphalt plant is well equipped, material engineers can monitor the status of the modification by measuring the viscosity of the binder during the blending process and stop it when the desired degree of interaction has been met. However, the correct degree of interaction must be known in advance since the same viscosity can be obtained in the swelling and degradation phases (Figure 1) but the mechanical characteristics are substantially different [14]. The ability to anticipate the result of the blending process for a given combination of binder and crumb rubber and a given set of processing conditions can represent a useful tool for practitioners. This research contributes to the ongoing efforts on this subject [19,20,21] and presents a mathematical equation for predicting viscosity as a function of mixing time, mixing temperature, and crumb-rubber characteristics.
Machine learning methodologies are extensively used to develop predictive algorithms for a material’s characteristics [19,20,21]. Despite their significant contribution, the algorithms generated by some of these methodologies are not easily implementable (e.g., neural networks, fuzzy logic algorithms) since they necessitate a knowledge of programming. For this reason, the model presented in this paper is developed by a machine learning approach named Multi-Gene Genetic Programming (MGGP), which returns an actual equation that can be easily implemented in any worksheet [22,23,24].

2. Principle of Genetic Programming

In 1992, Koza introduced the genetic programming (GP) approach as a branch of genetic algorithms (GA) [25]. GP is a computation technique that creates and evolves computer programs by implementing the “Darwinian natural selection” principles to obtain the best solution for a given problem. The algorithm randomly generates a primary population of computer programs (e.g., mathematical equations) structured like trees, where function nodes (e.g., arithmetic operations, Boolean logic functions, mathematical functions) are used to connect terminal nodes (variables and numerical constants). The goodness-of-fit of all the individuals in the computer programs’ population is evaluated, and the ones with superior performance are selected and carried to the next generation. The next generation is then created starting from those individuals and evolves by genetic operations such as mutation, crossover, and reproduction. Since this is a randomly based process, some of the individuals may move to the next generation without any evolution during reproduction, whereas others can evolve extensively [26].
Multi-Gene Genetic programming (MGGP) is an extension of the GP approach. The above-mentioned trees, also known as genes, which usually contain nonlinear terms, are here combined linearly with each other using the least square method. This helps to incorporate the power of both linear and nonlinear regression approaches to enhance the performance of the algorithm. However, the number of genes and their depths must be controlled by the user to avoid a further increase in equation complexity and decrease in algorithm efficiency.
In this research, the MGGP algorithm implemented in the GPTIPS software was used [27] as it allows changing some parameters of the algorithm manually to find an acceptable balance between the running time, output complexity, and goodness-of-fit of the models. The software settings were initially chosen based on recommendations provided in the literature [23,24] and then adapted to the database of materials’ characteristics used in this study and described in the next paragraph. As indicated in Table 1, the initial population size is one of the main inputs of the analysis. This represents the number of equations initially developed by the software and carried over the analysis. The final output is the series of equations that reached the final evolutionary step, each of them characterized by its goodness-of-fit and expressional complexity. Therefore, the user can either decide to consider all the equations or analyze the equations in a certain range of these two values.

3. Materials and Methods

Four crumb rubbers (CRs) sampled from three end-of-life tire (ELT) processing plants were characterized in a laboratory in terms of particle size distribution, density, and surface area. Two recycling plants used the so-called ’ambient’ particle size reduction process, whereas the third one processed ELTs under cryogenic conditions. The nominal maximum size of the CR particles produced with the ambient process was 1 mm and 0.5 mm for the samples named CR1 and CR2, respectively. The cryogenic product, named CR4, had a nominal maximum size of 0.6 mm. CR1, CR2, and CR4 were used as collected in the plant. CR3, instead, was obtained in a laboratory as a sub-product of CR1 by eliminating particles retained at a 0.354 mm sieve. This additional step was necessary to have information about the interaction effects due to CR with the same chemical composition but different gradations.
CR gradation was determined by a sieve analysis performed in dry conditions adopting the ASTM series sieve (1, 0.841, 0.71, 0.589, 0.5, 0.354, 0.25, 0.177, 0.125, 0.088, 0.063 mm). The results were then fitted to the Weibull distribution:
P d = 1 e ( d / λ ) k
where d is the sieve opening (in mm), P d is the percentage passing to the sieve through the opening d, and λ and k are the Weibull distribution parameters. The fitting process allowed the use of λ and k as the synthetic indicators of the shape curve distribution and uniformity of the particles’ distribution, respectively.
The particles’ density was measured using the pycnometer method, with ethyl alcohol employed as the reference fluid since its density is lower than that of rubber and therefore prevents CR particles from floating to the surface.
The assessment of the surface area (SA) per unit mass was based on digital images of the particles retained at each sieve size generated using a stereomicroscope. Pictures captured using a digital camera were subsequently processed through the software ImageJ. The following equation, proposed by Santagata et al. [15], was used for the SA estimation:
S A m = ϕ ( 6 / ρ ) Σ i ( f 1 / d m , i )
where ϕ is a correction factor, which is a function of the particles’ shape and roughness, ρ is the density (in g/m 3 ), f i is the frequency (in decimal units) of the i-th single-size fraction, d m , i is the mean particle diameter (in mm) of the i-th fraction.
CRM binders were prepared to combine each CR with a standard Pen 50/70 bitumen at three different dosages (15%, 18.5%, and 22% by weight of the base bitumen). Rubber particles stored at room temperature were added to the pre-heated base bitumen in batches immersed in a thermostatic oil bath in order to maintain a constant temperature of 190 C. Blends were then stirred by an anchor-shaped spindle for a total time of 210 min at a speed of 600 RPM. During the mixing process, the CRM binders were sampled at predefined time intervals (5, 15, 25, 35, 45, 60, 80, 100, 120, 150, 180, 210 min) and immediately subjected to viscosity tests to assess the evolution of viscosity as a function of the mixing time. The same mixing, sampling, and testing procedures were used for the blends produced at 150 C and 170 C with a crumb-rubber dosage of 18.5% only.
Rotational viscosity was measured with a Brookfield viscometer DVIII-Ultra equipped with an SC4-27 spindle. The modified bitumen was sampled from the blending container, poured into the viscometer cup, and tested at a single temperature of 175 C (as specified in the ASTM D-6114). Three subsequent steps of rotational speeds characterized by different durations were considered. In the first step, the rotational speed was set at 10 RPM (corresponding to a shear rate of 3.4 s 1 ) for six minutes to allow the transition from mixing to testing temperature and reach the steady-state flow conditions. The second and third stages had a duration of one minute, and the rotational speed was increased to 20 RPM (6.8 s 1 ) and 50 RPM (17 s 1 ), respectively. Only data recorded at 20 RPM are considered in this study since they correspond to the shear rate indicated in the ASTM specifications. Other data were used for different purposes.

4. Development of the Predictive Model

The development of the predictive model has been divided into two phases. In Phase I, the results of a Multi-Variable Least Square Regression (MLSR) model were compared to a preliminary MGGP-based model to highlight the benefits of using the MGGP algorithm. In addition, the variables’ inter-dependency and possible anomalies were analyzed. The MGGP algorithm was then used for a refined analysis in Phase II. The quality of the prediction was assessed by calculating the coefficient of determination ( R 2 ), mean squared error (MSE), and mean absolute error (MAE) as follows:
R 2 = 1 i = 1 n ( h i t i ) 2 i = 1 n ( h i h ^ ) 2
RMSE = i = 1 n ( h i t i ) 2 n
MAE = i = 1 n h i t i 2 n
where h i is the experimental outputs, t i is the calculated outputs, h ^ is the average of the experimental outputs, t i is the average of the calculated outputs, and n is the number of samples.

4.1. Phase I: Preliminary Analyses with MGGP and MLSR Approaches

The inter-dependency of the variables was considered to avoid overstating the effects on the mathematical model. An R-value of 0.8 between two variables was specified in previous studies as a good threshold [28] to discriminate variables’ inter-dependency. The results indicated high correlations between three parameters: λ , k, and SA. Hence, a database should contain only one of them or none. Four distinct databases were created for the modeling process:
  • Set 1: Weibull parameter λ , density ( λ , g/cm 3 ), mixing temperature (T, C), mixing time (t, min), and rubber content (RC, %),
  • Set 2: Weibull parameter k, density, mixing temperature, mixing time, and rubber content,
  • Set 3: Surface area (SA, mm 2 ), density, mixing temperature, mixing time, and rubber content,
  • Set 4: Weibull parameter k, mixing temperature, mixing time, and rubber content.
The set of data used for the preliminary MLSR and MGGP analyses had a population of 222 viscosity test results that were randomly divided into two subsets: 80% of the population was utilized for training and the remaining 20% was considered for validating (testing) the models. Since the purpose of the study was not the description of these materials from a mechanical point of view, only the statistical description of the data collected and used for modeling is herein reported (Table 2).
The MLSR and MGGP models that resulted from the preliminary analysis are
η = 112 λ 1600 ρ + 37.6 T 3.9 t + 1001.5 R C 17588.6
η = 110 T 1688 λ + 215 t 1688 R C 1577 L n ( t 3 ) + 122 R C · L n ( t 3 ) 0.913 T · t 0.913 T · R C 0.913 t · R C 2 + 2.33 R C 3 + 2.33 ρ · t · R C 490
The best MLSR- and MGGP-based equations are both functions of the rubber content, mixing time, mixing temperature, rubber density, and Weibull parameter λ (Set 1).
Figure 2 shows the relationship between the predicted and measured viscosity values using the MLSR and MGGP approaches. It is evident that the MLSR approach, given its mathematical form, was unable to develop an adequate equation to forecast the viscosity of the rubberized binders. In Figure 2a, data points are presented with three different symbols to highlight the gaps in the predictions for values between 2500 cP and 3200 cP as well as in the range 5800–8200 cP. On the other hand, the MGGP model shows a good overall goodness of fit. However, a difference existed for viscosity values above 6000 cP, which were mainly recorded for the binders with the highest crumb-rubber content. Thus, some of the variables considered in the analysis have a specific significance and impact the different mixing conditions, and there was a certain level of particle swelling and degradation. For this reason, the threshold value of 6000 cP was selected to split the initial database and refine the analysis in Phase II.

4.2. Phase II: Refined MGGP Model

In this second round of the modeling process, 24 values of the initial database were removed to be used for the final validation of the predictive models. The remaining 198 individual values were divided among those below (153 values) and above (45 values) 6000 cP. The statistical descriptions of these two sub-databases are provided in Table 3 and Table 4. The percentages for the subdivisions of each database for training and testing purposes were kept as per the previous step of the analysis.
As reported in Table 1, 600 equations were developed for every generation of the analysis. The frequency of occurrence of each variable was calculated (frequency equal to one indicates that all equations contain that variable). The results of the sensitivity analysis are illustrated in Figure 3. The mixing time and rubber content had the most dominant effect on the CRM binder viscosity prediction, followed by the mixing temperature and rubber density. The gradation also played an important role in forecasting the viscosity values. None of the equations for viscosity above 6000 cP considered the surface area. As mentioned in the previous section, high-viscosity values are usually recorded for a high rubber content, which may be a condition in which the surface area could lose its significance. The effect of a single parameter on viscosity can also be analyzed by a parametric study. Figure 4 illustrates the impact on the forecast viscosity due to the variation in the rubber density, Weibull parameter λ , and k. The parametric study of the surface area is not shown for the reason mentioned above. The effect of the gradation changed considerably in the two scenarios, whereas the rubber density had the same effect. However, the inherent randomness of the MGGP algorithm did not allow for the control of or force the effect of a given parameter in order to have the desired trend, and the algorithm optimized the prediction models as a result of the superposition of all the effects of the variables analyzed.

4.3. Predictive Models and Performance Analysis

For each data set listed previously, the best equation was selected based on the goodness of fit and the complexity of the analytical equation. The statistics of the MGGP equations are reported in Table 5. Based on the sum of the MAE for the training and test data, it can also be concluded that the models obtained with the third and fourth sets of inputs had the highest accuracy for predicting the CRM binder viscosity below and over the viscosity threshold. The final predictive models can be then expressed as in Equation (8). The development of two distinct MGGP predictive models for the CRM binder viscosity increased the overall goodness of fit as shown in Figure 5a,b.
η = 469 T 11,900 ρ + 20.3 t + 186 R C + 1944 t + 2 R C 1.28 T t 9.06 ρ + L n ( k ) 1.28 T 2 + 0.524 ρ 1.5 · t 2 ρ + T 43,700 ( f o r η < 6000 cP ) 8.08 k + T + R C 2 + 270 t + 1900 R C + 62,000 L n ( t ) 25,800 t 7377 T R C + 7377 R C + 53,800 t + 1344 t · L n ( t ) R C 1,230,000 ) ( f o r η > 6000 cP )
Equation (8) was used for predicting two sets of data that were not included in the analysis described above. This was intended to be an independent validation of the model. Figure 6 shows a comparison of the measured and predicted viscosities using the final MGGP models.

5. Conclusions

The magnitude of the interaction mechanisms between CR and bitumen during the different stages of the blending process is highly variable. The ability to forecast the effect of the mixing process for a given combination of binder and crumb rubber and a given set of processing conditions may be of great interest to practitioners and promote the use of rubber from ELTs for paving applications. For this purpose, the Multi-Gene Genetic Programming algorithm was used to develop a predictive model of CR-modified binders’ viscosities using the physical properties of the crumb rubber and the mixing process characteristics as the inputs.
The results obtained by the application of the MGGP algorithm on the whole database, as well as those obtained using the Multi-Variable Least Square Regression (MLSR) method, showed that the overall goodness of fit of the MGGP model ( R 2 = 0.89) was excellent compared to the model obtained through the MLSR approach. In particular, the MLSR appeared to be inadequate due to its mathematical form. A detailed analysis of the results showed a significant difference in predicting viscosity values below ( R 2 = 0.72) and above ( R 2 = 0.81) 6000 cP, which were mainly collected on binders modified with a high CR content. This suggested that certain characteristics of both the crumb rubber particles and mixing process, may have a specific significance and impact during the different blending phases. Based on these findings, the MGGP model was refined in order to have separate models for viscosity above and below the threshold value of 6000 cP. A sensible improvement of the goodness of fit was obtained and the comparison of the measured and predicted viscosities of the dataset used for the validation process showed excellent performance. The application of the Multi-Gene Genetic Programming algorithm is very promising since it overcame the problems encountered when using the MLSR approach. Moreover, unlike other machine learning algorithms, the MGGP returned an actual equation, which can be used by practitioners to forecast the interaction effects in terms of viscosity and tune the production process accordingly.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data, models, or code that support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Xiang, L.; Cheng, J.; Que, G. Microstructure and performance of crumb rubber modified asphalt. Constr. Build. Mater. 2009, 23, 3586–3590. [Google Scholar] [CrossRef]
  2. Lee, S.J.; Akisetty, C.K.; Amirkhanian, S.N. The effect of crumb rubber modifier (CRM) on the performance properties of rubberized binders in HMA pavements. Constr. Build. Mater. 2008, 22, 1368–1376. [Google Scholar] [CrossRef]
  3. Lee, S.J.; Akisetty, C.K.; Amirkhanian, S.N. Recycling of laboratory-prepared long-term aged binders containing crumb rubber modifier. Constr. Build. Mater. 2008, 22, 1906–1913. [Google Scholar] [CrossRef]
  4. Chiu, C.T.; Lu, L.C. A laboratory study on stone matrix asphalt using ground tire rubber. Constr. Build. Mater. 2007, 21, 1027–1033. [Google Scholar] [CrossRef]
  5. Navarro, F.J.; Partal, P.; Martínez-Boza, F.; Gallegos, C. Thermo-rheological behaviour and storage stability of ground tire rubber-modified bitumens. Fuel 2004, 83, 2041–2049. [Google Scholar] [CrossRef]
  6. Shatnawi, S. Performance of Asphalt Rubber Mixes in California. Int. J. Pavement Eng. 2001, 2, 1–16. [Google Scholar] [CrossRef]
  7. Bahia, H.; Davies, R. Effect of Crumb Rubber Modifiers (CRM) in Performance-Related Properties of Asphalt Binders. Asph. Paving Technol. 1994, 63, 414–438. [Google Scholar]
  8. Yang, X.; Yiu, Z.; Perram, D.; Hand, D.; Zeyad, H.; Wei, W.; Luo, S. Emission analysis of recycled tire rubber modified asphalt in hot and warm mix conditions. J. Hazard. Mater. 2019, 365. [Google Scholar] [CrossRef]
  9. Farina, A.; Zanetti, M.C.; Santagata, E.; Blengini, G.A. Life cycle assessment applied to bituminous mixtures containing recycled materials: Crumb rubber and reclaimed asphalt pavement. Resour. Conserv. Recycl. 2017, 117, 204–212. [Google Scholar] [CrossRef]
  10. Abdelrahman, M.; Carpenter, S. Mechanism of interaction of asphalt cement with crumb rubber modifier. Transp. Res. Rec. J. Transp. Res. Board 1999, 1661, 106–113. [Google Scholar] [CrossRef]
  11. Peralta, J.; Silva, H.M.R.D.; Hilliou, L.; MacHado, A.V.; Pais, J.; Christopher Williams, R. Mutual changes in bitumen and rubber related to the production of asphalt rubber binders. Constr. Build. Mater. 2012, 36, 557–565. [Google Scholar] [CrossRef]
  12. Artamendi, I.; Khalid, H.A. Diffusion Kinetics of Bitumen into Waste Tyre Rubber (with Discussion). J. Assoc. Asph. Paving Technol. 2006, 75, 133–164. [Google Scholar]
  13. Lo Presti, D.; Airey, G.D.; Partal, P. Manufacturing Terminal and Field Bitumen-Tyre Rubber Blends: The Importance of Processing Conditions. Procedia—Soc. Behav. Sci. 2012, 53, 485–494. [Google Scholar] [CrossRef] [Green Version]
  14. Lanotte, M.A. Evaluation of the Interaction between Bitumen and Crumb Rubber from ELTs in Modified Binders Produced by “Wet” Technology. Doctoral Dissertation, Politecnico di Torino, Torino, Italy, 2014. (In Italian). [Google Scholar]
  15. Santagata, E.; Dalmazzo, D.; Lanotte, M.; Zanetti, M.C.; Ruffino, B. Relationship between Crumb Rubber Morphology and Asphalt Rubber Viscosity. In Proceedings of the 5th Asphalt Rubber International Conference, Munich, Germany, 23–26 October 2012; pp. 513–532. [Google Scholar]
  16. Shen, J.; Amirkhanian, S.; Xiao, F.; Tang, B. Surface area of crumb rubber modifier and its influence on high-temperature viscosity of CRM binders. Int. J. Pavement Eng. 2009, 10, 375–381. [Google Scholar] [CrossRef]
  17. Cook, D. IBISWorld Industry Report OD5022—Tire & Rubber Recycling in the US; Technical Report September 2018; IBISWorld Publisher: Melburne, Australia, 2018. [Google Scholar]
  18. Ecopneus Scpa. Sustainability Report 2017; Technical Report; Ecopneus Scpa: Milano, Italy, 2017. [Google Scholar]
  19. Xiao, F.; Putman, B.; Amirkhanian, S. Viscosity prediction of CRM binders using artificial neural network approach. Int. J. Pavement Eng. 2011, 12, 485–495. [Google Scholar] [CrossRef]
  20. Specht, L.P.; Khatchatourian, O.; Brito, L.; Ceratti, J. Modeling of Asphalt-rubber Rotational Viscosity by Statistical Analysis and Neural Networks. Mater. Res. 2007, 10, 69–74. [Google Scholar] [CrossRef]
  21. Thodesen, C.; Xiao, F.; Amirkhanian, S.N. Modeling viscosity behavior of crumb rubber modified binders. Constr. Build. Mater. 2009, 23, 3053–3062. [Google Scholar] [CrossRef]
  22. Iqbal, M.F.; Liu, Q.; Azim, I.; Zhu, X.; Yang, J.; Javed, M.F.; Rauf, M. Prediction of mechanical properties of green concrete incorporating waste foundry sand based on gene expression programming. J. Hazard. Mater. 2020, 384, 121322. [Google Scholar] [CrossRef]
  23. Gandomi, A.H.; Alavi, A.H. A new multi-gene genetic programming approach to non-linear system modeling. Part II: Geotechnical and earthquake engineering problems. Neural Comput. Appl. 2012, 21, 189–201. [Google Scholar] [CrossRef]
  24. Gandomi, A.H.; Alavi, A.H. A new multi-gene genetic programming approach to nonlinear system modeling. Part I: Materials and structural engineering problems. Neural Comput. Appl. 2012, 21, 171–187. [Google Scholar] [CrossRef]
  25. Koza, J.R. Genetic Programming: On the Programming of Computers by Means of Natural Selection; MIT Press: Cambridge, MA, USA, 1992; p. 819. [Google Scholar]
  26. Koza, J.R. Genetic Programming: On the Programming of Computers by Means of Natural Selection. Stat. Comput. 1994, 4, 87–112. [Google Scholar] [CrossRef]
  27. Searson, D.P.; Leahy, D.E.; Willis, M.J. GPTIPS: An Open Source Genetic Programming Toolbox for Multigene Symbolic Regression. In Proceedings of the International Multiconference of Engineers and Computer Scientists, Hong Kong, China, 17–19 March 2010; Volume I, pp. 77–80. [Google Scholar]
  28. Smith, G.N. Probability and Statistics in Civil Engineering. Collins Prof. Tech. Books 1986, 244, 10007808566. [Google Scholar]
Figure 1. Particle size distribution of the CRs employed for laboratory preparation of AR binders.
Figure 1. Particle size distribution of the CRs employed for laboratory preparation of AR binders.
Sustainability 14 13798 g001
Figure 2. Preliminary results of binder viscosity prediction using the (a) MLSR and (b) MGGP models.
Figure 2. Preliminary results of binder viscosity prediction using the (a) MLSR and (b) MGGP models.
Sustainability 14 13798 g002
Figure 3. Frequency of input parameters: (a) η < 6000 cP, and (b) η > 6000 cP.
Figure 3. Frequency of input parameters: (a) η < 6000 cP, and (b) η > 6000 cP.
Sustainability 14 13798 g003
Figure 4. Parametric study: (a) η < 6000 cP, and (b) η > 6000 cP.
Figure 4. Parametric study: (a) η < 6000 cP, and (b) η > 6000 cP.
Sustainability 14 13798 g004
Figure 5. Predictive models for (a) η < 6000 cP—Inputs: k, ρ , T, t, RC, and (b) η > 6000 cP—Inputs: k, T, t, RC.
Figure 5. Predictive models for (a) η < 6000 cP—Inputs: k, ρ , T, t, RC, and (b) η > 6000 cP—Inputs: k, T, t, RC.
Sustainability 14 13798 g005
Figure 6. Results of the validation process.
Figure 6. Results of the validation process.
Sustainability 14 13798 g006
Table 1. Parameter settings for MGGP algorithm.
Table 1. Parameter settings for MGGP algorithm.
ParameterSettings
Function set + , , × , ÷ , s q r , p o w e r , L n
Population size600
Number of generations100
Maximum number of genes allowed in an individual8
Maximum tree depth8
Tournament size25
Elitism0.01% of population
Crossover events0.85
High-level crossover0.2
Low-level crossover0.8
Mutation events0.1
Sub-tree mutation0.9
Replacing input terminal with another random terminal0.05
Gaussian perturbation of randomly selected constant0.05
Direct reproduction0.05
Ephemeral random constants[−10 10]
Table 2. Descriptive statistics of the model variables for the entire database.
Table 2. Descriptive statistics of the model variables for the entire database.
Parameter λ kSA ρ TtRC η
[mm 2 ][g/cm 3 ][ C][min][%][cP]
Mean0.4093.462251.201177.286.118.54883
Median0.4022.172211.196190.070.018.54295
Maximum0.5806.213861.223190.0210.022.015,715
Minimum0.2541.911131.180150.05.015.0945
Std.Dev.0.1141.691020.01616.2165.82.133057
Skewness0.1560.720.670.084−0.7440.5240.0041.619
Kurtosis2.0541.971.971.8081.9281.9842.7075.632
Sum90.776749,951266.739,34019,11041001.1 × 10 6
Σ Sq.Dev.2.8946272.300.05558,0689.5 × 10 6 10042.1 × 10 9
Observations222222222222222222222222
Table 3. Descriptive statistics of the model variables for η < 6000 cP.
Table 3. Descriptive statistics of the model variables for η < 6000 cP.
Parameter λ kSA ρ TtRC η
[mm 2 ][g/cm 3 ][ C][min][%][cP]
Mean0.4133.430219.01.202173.992.317.83544
Median0.4133.726171.01.205190.080.018.53300
Maximum0.5806.213386.01.223190.0210.022.05975
Minimum0.2541.908113.01.180150.05.015.0945
Std.Dev.0.1071.61395.50.01717.6370.41.881372
Skewness0.1560.7780.8390.032−0.3920.37−0.130.123
Kurtosis2.3592.1652.3871.6281.4081.682.841.825
Sum68.82524.233,502183.826,61014,12027155.4 × 10 5
Σ Sq.Dev.1.730395.31.4 × 10 6 0.0424.7 × 10 4 7.5 × 10 5 537.52.9 × 10 8
Observations153153153153153153153153
Table 4. Descriptive statistics of the model variables for η > 6000 cP.
Table 4. Descriptive statistics of the model variables for η > 6000 cP.
Parameter λ kSA ρ TtRC η
[mm 2 ][g/cm 3 ][ C][min][%][cP]
Mean0.3933.230232.51.200186.965.320.99549
Median0.4022.167221.01.196190.060.022.08245
Maximum0.5806.213386.01.223190.0150.022.015,715
Minimum0.2541.908113.01.180170.05.018.56025
Std.Dev.0.1081.580100.50.0167.33141.81.643219
Skewness0.2500.9880.6120.184−1.900.34−0.820.533
Kurtosis2.3572.5581.9321.7114.6132.191.671.792
Sum17.87145.510,46154.01841029409414.2 × 10 5
Σ Sq.Dev.1.730395.31.4 × 10 6 0.0424.7 × 10 4 7.5 × 10 5 537.52.9 × 10 8
Observations4545454545454545
Table 5. Summary of the statistical performance of the MGGP models.
Table 5. Summary of the statistical performance of the MGGP models.
Set 1: λ , ρ , T, t, RC
R 2 RMSEMAE
η > 6000 cPTrain0.89443359
Test0.88543475
η < 6000 cPTrain0.95641531
Test0.881071903
Set 2: SA, ρ , T, t, RC
R 2 RMSEMAE
η > 6000 cPTrain0.89437353
Test0.88534458
η < 6000 cPNot found
Set 3: k, ρ , T, t, RC
R 2 RMSEMAE
η > 6000 cPTrain0.89434328
Test0.88486391
η < 6000 cPTrain0.95678582
Test0.881118990
Set 4: k, T, t, RC
R 2 RMSEMAE
η > 6000 cPTrain0.86490386
Test0.85604518
η < 6000 cPTrain0.94749327
Test0.891040917
Gray highlights represent the selected models.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lanotte, M. Soft Computing Approach for Predicting the Effects of Waste Rubber–Bitumen Interaction Phenomena on the Viscosity of Rubberized Bitumen. Sustainability 2022, 14, 13798. https://doi.org/10.3390/su142113798

AMA Style

Lanotte M. Soft Computing Approach for Predicting the Effects of Waste Rubber–Bitumen Interaction Phenomena on the Viscosity of Rubberized Bitumen. Sustainability. 2022; 14(21):13798. https://doi.org/10.3390/su142113798

Chicago/Turabian Style

Lanotte, Michele. 2022. "Soft Computing Approach for Predicting the Effects of Waste Rubber–Bitumen Interaction Phenomena on the Viscosity of Rubberized Bitumen" Sustainability 14, no. 21: 13798. https://doi.org/10.3390/su142113798

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