Circular Image

S.M.J.G. Erkens

info

Please Note

183 records found

Journal article (2026) - W. A.A.S. Premarathna, Kumar Anupam, M. Moenielal, Thijs Wensveen, Cor Kasbergen, Sandra M.J.G. Erkens
Accurate prediction of rolling resistance (RR) is essential for improving vehicle fuel efficiency and supporting policymakers in making sustainable environmental decisions. This study introduces a novel framework that integrates both data-driven and physics-based approaches to enhance RR prediction by incorporating tire-penetration level indicator, the Delta (δ) parameter. The research investigates the relationships between RR, the δ parameter, and texture properties to refine predictive modelling. A portable device was built to measure the in-field δ parameter using tire-pavement interaction. Machine learning (ML) techniques, including multiple linear regression (MLR), random forest regressor (RFR), artificial neural networks (ANN) and finite element method-based (FEM) tire-pavement interaction models were employed to develop and validate the framework. Findings from the FEM tire-pavement interaction model confirmed the reliability of the δ parameter. Exploratory data analysis (EDA) highlighted the strong correlation between texture metrices such as MPD, ETD, and RMS, reinforcing the δ parameter's role in tire-pavement interactions. Comparative analysis of different pavement surfaces revealed that worn surfaces contribute to higher δ parameter values and increased RR. The improvement resulting from the inclusion of the δ parameter is particularly evident in the ANN and RF models, confirming nonlinear interaction effects between tire penetration and surface texture. It was also observed that the obtained RR data follow a non-normal distribution, which most of the previous studies did not consider. A deeper statistical insight showed that the δ parameter has a significant impact on RRC prediction. The primary contribution of this study lies in demonstrating the feasibility of integrating a physics-based tire-pavement interaction parameter into ML models for rolling resistance prediction, thereby bridging mechanistic modelling and machine learning within pavement engineering. ...
Accurate replication of pavement textures is crucial for tire-pavement interaction studies in laboratory conditions, particularly at the macro-texture level. Existing replication methods often rely on costly equipment and lack standardized procedures for validating surface replicas. This study presents a novel methodology to replicate pavement surfaces and validation. The methodology integrates Mean Texture Depth (MTD) as a 3D texture indicator, silicone rubber casting technology, Finite Element (FE)-based 3D geometric calculations, and smartphone-video-based digital modelling using computer vision. Different asphalt mixtures with distinct surface textures were used for the development of the physical and digital replicas. Physical replicas were produced through silicone casting and epoxy molding, while digital replicas were generated using the Structure-from-Motion technique. Validation was conducted using the standard sand patch test for physical samples and FE-based MTD calculations for digital models. The results showed strong agreement between the replicas and the original pavement samples. Statistical analyses confirmed normal distribution, variance homogeneity, and no significant differences between digital and physical samples. The findings of this study are expected to benefit road authorities in the design and laboratory testing of pavement surfaces, with the aim of evaluating pavement performance. ...
Journal article (2025) - Mahmoud Khadijeh, Veronica Cerqueglini, Cor Kasbergen, Sandra Erkens, Aikaterini Varveri
Physics Informed Neural Networks (PINNs) have been rarely applied to solve multiphysics systems due to the inherent challenges in optimizing their complex loss functions, which typically incorporate multiple physics-based terms. This study presents a multistage PINN approach designed to efficiently solve coupled multiphysics systems with strong interdependencies. The multistage PINN progressively increases the complexity of the physical system being modeled, enabling more effective capture of coupling between different physics. The computational merits of this approach are demonstrated through two illustrative applications: prediction of asphalt aging and modeling of lid-driven cavity flow. Quantitative and qualitative comparisons with standard PINN and adaptive weight PINN approaches demonstrate the enhanced precision and computational efficiency of the proposed algorithm. The multistage PINN achieves a reduction in training time of more than 90% compared to standard PINNs while maintaining better alignment with the finite element method (FEM) solutions. The improvement in computational efficiency, coupled with enhanced accuracy, positions the multistage PINN as a powerful tool for addressing complex multiphysics problems across various engineering disciplines. The method’s ability to handle interactions between multiple physical processes, such as diffusion, chemical reactions, and fluid dynamics, makes it suitable for simulating long-term material behavior and complex fluid systems. ...
The chemo-mechanical properties of bitumen undergo significant alternations during aging and rejuvenation, posing challenges for accurately evaluating and enhancing rejuvenation efficiency in asphalt recycling. This study investigates how bitumen source, aging degree, rejuvenator type and dosage influence the chemical and rheological performance of rejuvenated bitumen. Comprehensive characterizations are performed using saturate, aromatic, resin, and asphaltene (SARA) fractionation, elemental analysis, gel permeation chromatography (GPC), and dynamic shear rheometer (DSR) tests. To elucidate chemo-rheological correlations, statistical techniques (Pearson correlation, analysis of variance (ANOVA), and Chi-square tests) are combined with artificial neural networks (ANN). Results indicate that the NB bitumen with more colloidal stability and less sulfur content exhibits the highest resistance to long-term aging. FB bitumen with 4.3 % sulfur achieves the best high-temperature deformation resistance with rutting failure temperature (RFT) higher than 80 °C, and TB bitumen exhibits the longest fatigue life. Rejuvenation using bio-oil is most effective on reducing relaxation time by up to 60 % and increasing creep compliance (Jnr3.2) by 1.7–2.5 times, depending on bitumen type. Rejuvenator dosage sensitivity for relaxation stress follows the trend: bio-oil < engine-oil < naphthenic-oil, while aromatic-oil shows variability depending on its source. Among the tested rejuvenators, bio-oil proves most effective, particularly for rejuvenating TB and FB bitumen. The ANN model demonstrates strong predictive performance for rheological properties, achieving R2 values between 0.90 and 0.98, with the highest accuracy observed for relaxation indices, followed by fatigue and rutting properties. ...

A python-based software package for initializing and running molecular dynamics simulations using LAMMPS

Journal article (2025) - Eli I. Assaf, Elsa Maalouf, Xueyan Liu, Peng Lin, Sandra Erkens
Scymol is a Python-based software package specifically designed to facilitate the setup and execution of molecular simulations in LAMMPS. It comes equipped with a user-friendly interface, which simplifies the process of initializing molecular systems and defining simulation parameters. Moreover, the software generates and executes LAMMPS simulation sequences, enabling researchers to establish comprehensive simulation schemes, such as heating or deformation cycles, in a single run. Through its successful application in diverse research projects and its modular design, Scymol demonstrates considerable promise as an indispensable tool for researchers aiming to carry out molecular dynamics simulations without sacrificing complexity or high-throughput capabilities in their methodologies. ...
Journal article (2025) - W. A.A.S. Premarathna, Kumar Anupam, M. Moenielal, Thijs Wensveen, Cor Kasbergen, Sandra M.J.G. Erkens
A major contributor to GHG emissions is the transportation sector, particularly pavement transport. The limited understanding of tire-pavement interactions leads to inaccurate predictions of these emissions, particularly from rolling resistance (RR). Traditional methods for predicting RR are constrained by their limited applicability and inability to account for the complex dynamics of tire-pavement interactions, resulting in poor prediction accuracy. These limitations make it challenging for policymakers to make proper decisions, as existing methods are manual and labour-intensive. This study aims to develop an automated system to capture tire-pavement interaction data using the Laser Crack Measurement System (LCMS). To the best of the authors' knowledge, no robust technique currently exists for automatically calculating tire penetration-related information from LCMS data to predict RR. Therefore, this research explores machine learning (ML)-based models to reduce uncertainties in existing approaches and enhance RR predictions using automated LCMS data. It examines the relationships between RR, tire penetration volume, and the characteristics of the Dutch pavement network, comparing the results with those of commonly used RR prediction models. The study introduces an automatic tire penetration calculation approach using LCMS data to assess the impact of tire penetration volume and depth on RR in relation to surface properties. The findings reveal that traditional empirical models show poor correlations between RR and texture indicators, whereas ML-based models significantly improve the accuracy of RR predictions. These results could inform the development of strategies to reduce GHG emissions from pavement transport, supporting global efforts to combat climate change and achieve the goals of the Paris Agreement. ...
Journal article (2025) - Kumar Anupam, Mohammadjavad Berangi, Juan Camilo Camargo, Cor Kasbergen, Sandra Erkens
Asphalt mixtures show complex mechanical behavior due to their heterogeneous structure. Traditionally, the mechanical characterization of asphalt mixture is done through laboratory testing or micromechanical modeling. While laboratory tests and micromechanical models provide reliable measurements and physical interpretability, they are often resource-intensive and demand extensive calibration. Recent advances in machine learning address some of the above issues by enabling accurate predictions, though often lacking physical interpretability and stability. Hence, this study aims to present a novel micromechanics-infused neural network (MINN) framework for predicting asphalt mixture stiffness. The framework embeds micromechanical principles derived from the modified Hirsch model into the neural network's loss function, allowing the model to learn from experimental data while adhering to micromechanics-based constraints. In this study, feature selection is performed using BorutaShap, and Bayesian optimization is applied for hyperparameter tuning. Results show that MINN improves prediction accuracy, interpretability, and robustness. ...
Journal article (2025) - Panos Apostolidis, Robbert Naus, Xueyan Liu, Remco Hermsen, Sandra Erkens, Tom Scarpas
Pavement materials that could enhance the mechanical properties of open-graded porous asphalt mixtures in long-term service periods could offer a solution to produce long-life pavements, causing a reduction of interventions' needs, as well as the associated disruptions to road users and user costs. One option to improve the longevity of open-graded porous mixtures is with the use of epoxy asphalt that, despite its high initial cost, offers enhanced longevity that might offset any future user and intervention costs. This study aimed to evaluate the durability of plant-produced epoxy-modified open-graded porous asphalt mixtures. A batch production plant was employed to produce loose mixtures, which were used to pave a test road in the Province of Gelderland, the Netherlands, and compact specimens in the laboratory. Control mixtures with a non-epoxy-modified asphalt binder were also produced in the same plant. The durability of laboratory- and field-compacted mixtures was evaluated by conducting indirect tensile tests before and after oven conditioning. Results illustrated that the epoxy-modified asphalt demonstrated the highest strength and stiffness values, while the strength was reduced after conditioning in a water bath with the retained strength within the allowable specification limits. This attribute was confirmed from drill cores obtained from the test road after one year in service. Also, the materials compacted in the field had slightly higher strength and stiffness values than the laboratory-produced mixtures. Although the results provided have illustrated the improvement of durability of open-graded porous asphalt with implementing epoxy modification, further evidence from the test road over the years is needed for validation. ...
Journal article (2025) - A. Jagadeesh, W. A.A.S. Premarathna, A. Kumar, C. Kasbergen, S. Erkens
This research addresses the critical issue of load transfer efficiency (LTE) in jointed plain concrete pavements (JPCP), with a specific focus on the role of dowel bars in ensuring optimal load transfer and providing a comfortable ride for vehicles. While experimental studies have investigated factors like joint width, slab thickness, concrete strength, and dowel bar size that influence LTE, they are limited in their ability to accurately replicate real-world conditions and can be time-consuming. To overcome these limitations, finite element modelling (FEM) is employed as a powerful tool for simulating complex loading conditions and analyzing stress and strain distributions in pavements. The primary objective of this research is to develop an advanced FE model that incorporates the forklift tire-pavement interaction, enabling precise analysis of complex loading conditions in industrial pavements and the impact of various rigid pavement parameters on load transfer. By explicitly considering the interaction between the tire and pavement, the proposed model will provide an extensive and robust numerical tool for designers and engineers. Additionally, this study represents a novel framework to integrate concrete pavement dowel bars and complex tire modelling using FEM. The developed methodology holds significant promise in optimizing the design of dowel bar systems and back-calculating the pavement parameters for rolling weight deflectometers. ...
Journal article (2025) - Eli I. Assaf, Xueyan Liu, Sandra Erkens
Previous work demonstrated that Random Forest Regressors (RFRs) could estimate the physical properties of bitumen using molecular descriptors derived from Molecular Dynamics (MD) simulations, thereby reducing the need for computationally intensive simulations. However, due to their decision-tree structure, RFRs lack true predictive capabilities, particularly for interpolation and extrapolation beyond the training data.

This study advances that foundation by employing Artificial Neural Networks (ANNs), which—when properly trained—can capture complex relationships with greater continuity and generalizability. Beyond simply replacing RFRs, we develop a fully automated framework for constructing Machine Learning Models (MLMs) to predict density and thermal expansion coefficients of bitumen. Using Optuna for hyperparameter optimization, we ensure that the information extracted from MD simulations is utilized effectively.

The resulting ANN models accurately reproduce MD-predicted densities, achieving R2>0.99, MSEs below 0.1 %, and maximum absolute errors below 5 % on test data. In addition to reducing computational cost, the models exhibit improved interpolation and extrapolation capabilities, enabling reliable predictions for properties, ranges, and compositions not explicitly simulated.

Key aspects of our approach include:
• Transitioning from RFRs to ANNs, improving generalization, interpolation, and predictive accuracy.
• Automated hyperparameter optimization, leveraging Optuna to maximize model efficiency.
• Expanding applicability, enabling property prediction for unseen compositions without additional MD simulations. ...
Asphalt binders are critical for asphalt pavement performance, and understanding their rheological behavior is essential for designing durable roadways. The complex shear modulus (G⁎) and phase angle (δ) are primary parameters characterizing binder rheology. This study introduces a novel hybrid machine learning model combining deep neural networks (DNN) and Gaussian process regression (GPR) to predict G⁎ and δ for bituminous binders and binder-filler systems (mastics). DNN excel at capturing complex, nonlinear relationships among eleven binder and thirteen mastic input parameters, including aging conditions, chemical and physical properties, and test parameters. However, standalone DNN struggle with small datasets, common at the binder scale, and lack inherent uncertainty quantification, limiting reliability in engineering applications. GPR improves DNN by refining predictions through probabilistic modeling, while providing uncertainty estimates, and enhancing accuracy with limited or noisy data. The hybrid model leverages DNN's feature extraction capabilities and GPR's ability to smooth predictions, significantly improving performance over standalone DNN. The hybrid model achieves high prediction accuracy, with R2 values of 0.997 for G⁎ and 0.947 for δ for binders, and 0.993 for G⁎ and 0.972 for δ for mastics, reducing G⁎ prediction error from 22.7% to 0.031% for fresh asphalt binder compared to standalone DNN. Feature importance analysis using random forest and SHAP techniques identifies test temperature, aging conditions, and penetration as key influencers of G⁎ and δ. This hybrid approach enhances the characterization of complex asphalt materials, offering pavement engineers a robust, reliable tool for predicting material behavior under diverse conditions. ...
Journal article (2024) - Eli I. Assaf, Xueyan Liu, Peng Lin, Shisong Ren, Sandra Erkens
This study explores the use of chemical descriptors derived from force field atom types to predict Fickian diffusion coefficients of rejuvenators in bitumen, utilizing machine learning models trained on data from 240 non-equilibrium molecular dynamics simulations. The simulations cover three bitumen types (NO, TO, FO), five aging degrees, and four temperatures (60 °C, 120 °C, 160 °C, 200 °C), capturing diffusion coefficients ranging from 0.0068e-10 m2/s in highly aged bitumens at 60 °C to 4.35e-10 m2/s in fresher samples at 200 °C. The MLM, built with 18 chemical descriptors for bitumen and rejuvenator sides, achieves an R2 of 0.97, accurately predicting diffusion across varied conditions. This approach abstracts away from the need for repeated MD simulations, enabling diffusion predictions even for systems outside the original dataset. The manuscript presents three case studies to illustrate how the model can be used for the iterative design of rejuvenators by optimizing molecular structures based on critical chemical features, such as rejuvenator oxygen content, bitumen sulfur content, and molecular weights. It also demonstrates how the model offers a practical framework for understanding the diffusion and performance of rejuvenators by linking time-dependent factors—such as concentration, depth, and rejuvenation time—with the bulk properties of bitumen-rejuvenator systems, facilitating industrial applications. ...
Journal article (2024) - Mohammadjavad Berangi, Bernardo Mota Lontra, Kumar Anupam, Sandra Erkens, Dave Van Vliet, Almar Snippe, Mahesh Moenielal
Inconsistencies between performance data from laboratory-prepared and field samples have been widely reported. These inconsistencies often result in inaccurate condition prediction, which leads to inefficient maintenance planning. Traditional pavement management systems (PMS) do not have the appropriate means (e.g., mechanistic solutions, extensive data handling facilities, etc.) to consider these data inconsistencies. With the growing demand for sustainable materials, there is a need for more self-learning systems that could quickly transfer laboratory-based information to field-based information inside the PMS. The article aims to present a future-ready machine learning-based framework for analyzing the differences between laboratory and field-prepared samples. Developed on the basis of data obtained from field and laboratory data, the gradient-boosting decision trees-based framework was able to establish a good relationship between laboratory performance and field performance (R2test > 80 for all models). At the same time, the framework could also show more complex relationships that are often not considered in practice. ...
Journal article (2024) - Eli I. Assaf, Xueyan Liu, Sandra Erkens
This study employs strain-controlled oscillatory deformations in Molecular Dynamics (MD) simulations to evaluate the dynamic properties of all-atom molecular systems, specifically targeting the SARA fractions of bitumen. Twelve molecular systems representing these fractions were modeled using the PCFF force field. The simulations effectively captured their viscoelastic properties across multiple frequency domains, including Elastic, Glassy, Rubbery, and Viscous responses. Reported storage and loss moduli range from thousands to tens of megapascals, with viscosities from tens to near-zero Pascal-seconds across various frequencies and temperatures, aligning well with experimental observations. Saturates and Aromatics were identified as the softest and most thermally susceptible fractions, while Resins and Asphaltenes were the stiffest and least susceptible. The study reveals that the relaxation time of all-atom molecular systems is significantly shorter than in experimental setups, necessitating careful comparison of stress-related phenomena across equivalent relaxation times. Although this allows for the exploration of response profiles in computationally tractable simulations, the nature of all-atom force fields and simulation algorithms introduces spatiotemporal scale discrepancies that must be addressed in future simulations involving the study of stress-related phenomena using MD. ...

A self-contained Python tool to generate atomistic systems of organic molecules using their SMILES notations

Journal article (2024) - Eli I. Assaf, Xueyan Liu, Peng Lin, Sandra Erkens
The advent of computational techniques, particularly atomistic simulations, has lessened the dependency on physical experiments in various scientific fields. Yet, the preparation complexity for simulations using platforms like LAMMPS and GROMACS persists. We introduce SMI2PDB, a Python tool that automates molecular systems assembly from SMILES to PDB format, easing molecular dynamics simulation setups. SMI2PDB manages molecule configuration and quantification effortlessly, establishes stable conformers, applies random rotations, and positions them in a simulation box with a Sobol sequence to reduce overlaps. This script facilitates the rapid preparation of complex organic mixtures for use in simulations, enhancing the exploration of novel materials. ...
Journal article (2024) - Chen Wang, Kumar Anupam, Cor Kasbergen, Sandra Erkens
The linear viscoelastic behavior of materials is represented using mechanical models of choice, which are further utilized in different numerical investigations, such as finite element simulations and discrete element simulations. Burger's model is one of the widely adopted mechanical models and remains highly favored in contemporary research due to its multiple advantages. Specifically, it excels in representing long-term creep and stress relaxation behavior in a relatively simplified manner. Accurate identification of the long-term behavior for the viscoelastic material, particularly asphalt concrete, is crucial, as it serves as a key indicator of asphalt pavement performance over its service life. However, past research studies show that the parameters of Burger's model should be back-calculated from experimental data only within a limited range of frequency, otherwise, the parameters fail to represent the true material behavior. To the best of the authors’ knowledge, there is no approach for researchers to obtain the critical frequency range in which the experiments should be performed. Therefore, this study proposes a novel framework to find the critical frequency range to obtain appropriate model parameters of Burger's model, to better characterize the viscoelastic behavior of the materials. To examine the framework, asphalt concrete mixtures are used as examples in this study. Necessary laboratory tests including complex modulus tests and stress relaxation tests, are performed on two distinctive types of asphalt concrete mixtures. The generalized Maxwell model with different number of Maxwell chains are used to evaluate the performance of Burger's model. Furthermore, since commercially available finite element packages generally do not have a direct built-in Burger's model, the article shows a way of implementing Burger's model in finite element simulation. The simulations corresponding to the laboratory tests are carried out in both frequency domain and time domain to thoroughly evaluate the performance of Burger's model. The optimal frequency range of 0.1–20 Hz for the examined mixtures is found to significantly improve the accuracy of the descriptive master curve. The results also suggest that the generalized Maxwell model requires a minimum of four Maxwell chains to maintain good performance in accurately characterizing the behavior of asphalt mixtures. However, adding more Maxwell chains beyond a critical limit may not provide significant benefits. Finite element simulations demonstrate that the stress relaxation behavior predicted by the obtained Burger's model parameters aligns more closely with experimental data over longer time intervals. This makes Burger's model a strong choice for aiding in the design of simulations for studies focused on the long-term behavior of materials. ...
Journal article (2024) - Shisong Ren, Xueyan Liu, Sandra Erkens
The low-temperature relaxation and fatigue cracking performance are two essential aspects in estimating the rejuvenation efficiency of recycling agents (RAs). This study aims to fundamentally investigate the effects of recycling agent type/dosage and aging degree of bitumen on thermodynamic and rheological properties of rejuvenated bitumen at low and intermediate temperatures. Molecular dynamics (MD) simulations are utilized to predict thermodynamic indices of rejuvenated bitumen, further linked to critical low-temperature and fatigue indicators from experiments. The results reveal that all RAs show a regeneration effect on fractional free volume (FFV), self-diffusion coefficient (DS), glass transition temperature (Tg), and surface free energy (γ). Bio-oil and engine-oil exhibit higher rejuvenation efficiency on these thermodynamic properties than naphthenic-oil and aromatic-oil. The aging degree of bitumen and temperature show significant effects on rejuvenation efficiency. It is recommended to use the FFV parameter to predict the relaxation properties of rejuvenated bitumen. However, these thermodynamic indicators inadequately differentiate between rejuvenators and softeners. Based on crossover parameter results, most recycling agents (bio-oil, engine oil, and naphthenic oil) in this study display softening characteristics. Only aromatic oil effectively rejuvenates the crossover modulus (Gc) of aged binder. Notably, engine oil demonstrates the least rejuvenation in crossover parameters for the recovery of aged bitumen. Further, γ demonstrates a strong association with both Glover-Rowe (G-R) and fatigue crack width C500 indices across all cases involving rejuvenated bitumen. This work will build a multi-scale evaluation framework on the rejuvenation effectiveness of recycling agents on low-temperature and fatigue performance of aged bitumen. ...
Journal article (2024) - Zili Wang, Panchamy Krishnakumari, Kumar Anupam, Hans van Lint, Sandra Erkens
The relationship between real-world traffic and pavement raveling is unclear and subject to ongoing debates. This research proposes a novel approach that extends beyond traditional correlation analyses to explore causal mechanisms between mixed traffic and raveling. This approach incorporates the causal discovery method, and is applied to five Dutch porous asphalt (PA) highway sites that have substantial data sets. Findings indicate a nonlinear relationship between traffic volume and raveling, with road age emerging as a shared contributor. The results also suggest that the degree to which different vehicle types contribute as a causal factor for raveling varies with carriageway configurations and lane characteristics. This underlines the need for targeted maintenance strategies. Challenges remain due to confounding correlations among traffic variables, necessitating further development of causal discovery models. This study may not conclusively resolve the debate on to what extent traffic contributes to raveling, but we argue we provide sufficient evidence against rejecting this hypothesis. ...

Automating LAMMPS data file generation from PDB molecular systems using Python, Rdkit, and Pysimm

Journal article (2024) - Eli I. Assaf, Xueyan Liu, Peng Lin, Sandra Erkens
Pdb2dat, developed in Python, is an open-source, self-contained utility that facilitates the conversion of PDB files into LAMMPS data files, catering to the need of initializing atomistic simulation from initial atomic configurations. It extracts molecular details from PDB files, uses Rdkit and Xyz2mol for bonding analysis and 3D conformer generation, and uses Pysimm for assigning force field types and charges. Designed to be lightweight and fully Pythonic, pdb2dat is suitable for use in privilege-limited high-throughput environments. The output details system topologies for use in MD simulations, significantly simplifying the preparatory steps needed by researchers to explore materials phenomena through LAMMPS. ...
Journal article (2024) - Sadaf Khalighi, Sandra Erkens, Aikaterini Varveri
Bituminous binders naturally age, affecting the properties and performance of asphalt pavements. The physical and chemical characteristics of binders are influenced by environmental factors, leading to a decline in their performance and durability. Therefore, it is essential to understand the mechanisms of binder environmental aging to design more resilient and long-lasting asphalt pavements. This study examines the effects of temperature, liquid, and vapor water on binder aging to develop more durable pavements. Aging was induced at three temperatures (60°C, 70°C, 85°C) under dry air, 90% relative humidity, and water immersion conditions. Field-aged samples were also analyzed to compare with laboratory-aged samples. Fourier-transform infrared spectroscopy (FTIR) and dynamic shear rheometer (DSR) were used to assess chemical and rheological changes. To assess the similarity between samples and identify the lab aging protocol closest to field aging, Hierarchical Cluster Analysis (HCA) was employed for data analysis. ...