This paper aims to incorporate the fractional derivative viscoelastic model into a finite element analysis. Firstly, based on the constitutive equation of the fractional derivative three-parameter solid model (FTS), the constitutive equation is discretized by using the Grünwald–Letnikov definition of the fractional derivative, and the stress increment and strain increment relationship and Jacobian matrix are obtained by using the difference method. Subsequently, we degrade the model to establish stress increment and strain increment relationships and Jacobian matrices for the fractional derivative Kelvin model (FK) and fractional derivative Maxwell model (FM). Finally, we further degrade the fractional derivative viscoelastic model to derive stress increment and strain increment relationships and Jacobian matrices for a three-component solid model and Kelvin and Maxwell models. Based on these developments, a UMAT subroutine is implemented in ABAQUS 6.14 finite element software. Three different loading modes, including static load, dynamic load, and mobile load, are analyzed and calculated. The calculations primarily involve a convergence analysis, verification of numerical solutions, and comparative analysis of responses among different viscoelastic models.

Currently, the research on the mechanical properties of rubber-modified asphalt mixtures primarily focuses on small-scale investigations, with insufficient exploration into the performance of rubber particles and their relationship with the mechanism and properties of modified asphalt mixtures. Limited studies have been conducted on large-scale rubber modification in asphalt mixtures. Due to frequent use and subsequent high damage to existing asphalt pavements, incorporating rubber-modified asphalt mixtures can partially alleviate premature deterioration. Dynamic modulus tests were conducted using MTS equipment under unconfined conditions to investigate the viscoelastic behavior of rubber-modified asphalt mixtures with high rubber content and elucidate the influence of rubber particle content on the elastic deformation and recovery capability. The dynamic mechanical properties of the mixtures were determined at different loading rates, temperatures, and types of rubber-modified asphalt mixtures. Based on the test data, variations in the dynamic modulus, phase angle, storage modulus, loss modulus, loss factor, and rut factor of the rubber-modified asphalt mixtures under different loading frequencies, temperatures, and types were analyzed. The results demonstrate the pronounced viscoelastic behavior of rubber-modified asphalt mixtures. The mixtures exhibit enhanced elasticity at low temperatures and high frequencies, while their viscosity becomes more prominent at high temperatures and low frequencies. Under constant test temperatures, an increase in load loading frequency leads to a higher dynamic modulus; conversely, a decrease in dynamic modulus is observed with increasing test temperatures. The dynamic modulus of ARHM-25 at a frequency of 10 Hz is found to be 12.99 times higher at 15 °C compared to that at 60 °C, while at 30 °C, the dynamic modulus at 25 Hz is observed to be 2.72 times greater than that at 0.1 Hz. Furthermore, the rutting resistance factors of the asphalt mixtures increase with loading frequency but decrease with temperature. The rutting factor for ARHM-13 at a frequency of 10 Hz is found to be 22.98 times higher at 15 °C compared to that at 60 °C, while at a temperature of 30 °C, the rutting factor for this material is observed to be 3.09 times greater at a frequency of 25 Hz than at 0.1 Hz. These findings suggest that rutting is most likely when vehicles drive at low speeds in hot weather conditions.

The bearing plates used in plate load test for highway engineering are typically rigid. However, due to limitations in obtaining the accurate distribution of compressive stress at the bottom of the bearing plate, there is often a significant deviation between the measured subgrade resilient modulus and the actual condition. To address this issue, a flexible bearing plate can be used to test the subgrade and obtain a more accurate resilient modulus. In this study, we use variance and degree of mean deviation to quantitatively evaluate the distribution uniformity of compressive stress. To create a rigid-flexible bearing plate that is similar to a flexible bearing plate, we explore the combinatorial design of steel plates and rubber mats. We examine factors such as the thickness (10, 20, and 30 mm) of the steel plate, elastic modulus (5, 10, and 20 MPa) and thickness (10, 20, and 30 mm) of the rubber mat, friction coefficient (μ:0, 0,2, 0.4, 0.6, 0.8, ∞) between the bearing plate and subgrade, and the combined shape characteristics of the rubber mat and steel plate. To reduce friction between the rubber mat and subgrade, we use lubricant, and through our design process, we develop a flexible bearing plate with relatively uniform compressive stress. Our computations show that when μ = 0.05, the variance is 0.0001, and the degree of mean deviation is 0.0780. These results indicate that the distribution uniformity of the compressive stress is very close to the uniform distribution load, which meets the necessary accuracy requirements for engineering applications.