Theoretical Analysis and Experimental Evaluation of Road Dynamic Behavior on Different Subgrades

,e dynamic characteristics of a road on a subgrade filled with different materials under traffic loads are studied. An elastic double-layer plate model on a viscoelastic subgrade is established to describe the dynamic behavior of the road.,e load diffusion path in the form of a cone in the stabilised layer under single vibration source and three vibration sources is considered, and the equations for the displacement and effective soil pressure of the elastic double-layer plate on the viscoelastic subgrade are derived. ,e relevant empirical formulae are then properly revised to reflect the dynamic behavior of the road.,e dynamic response of the road is analysed by adjusting factors such as the excitation load, excitation frequency, and subgrade filling materials. ,e results show that gravel as a subgrade filling material has a better vibration reduction effect than sand bag. In addition, the test data not only have good consistency with the theoretical formulae but also validate the reliability of the empirical formulae.


Introduction
e dynamic responses of a road subjected to traffic loads have been widely studied in civil engineering.Most of the world's highways are paved using an asphalt layer, which is normally laid on a subgrade consisting of two or more layers of different materials.e movement of vehicles over the pavement may cause vibration and deformation in the subgrade, requiring high-cost maintenance and reducing the service life of highways.us, it is important to accurately predict the dynamic response of highway structures considering the vehicle-pavement interaction.Numerous numerical simulations and testing models have been proposed to predict the dynamic responses of highway pavement and subgrade.
Several simple discrete element models have been widely applied to investigate the effects of vehicle vibration on roads.Oscarsson and Dahlberg [1] developed a simplified numerical method to study the vertical dynamic characteristics of a railway track under moving loads.Kouroussis et al. [2] proposed a coupled lumped mass (CLM) model to study the dynamic interaction between a track and a soil.Zhai and Cai [3] developed a numerical model to simulate the dynamic response between a train and a railway track.Although these models do not accurately reflect the dynamic response of a highway structure, simple discrete models are acceptable for solving practical engineering problems.
Compared to the discrete element models, the viscoelastic models are more reasonable for the dynamic response analysis of highways or railway structures.Generally, an elastic structure is filled in layers such that the elastic model can more effectively reflect the characteristics of the structure, and analytical methods can be used to obtain the dynamic response of layered elastic structures.For example, Lu et al. [4] used the reflection matrix (TRM) method to study the dynamic characteristics of a layered road-subgrade system.e results showed that the TRM could effectively predict the stress and multilevel elastic deformation of the road structure.Lih and Mal [5] developed a multiple transformation technique to obtain double integral expressions for the displacement and stress in a composite laminate under dynamic loading.
ey found that this technique can more effectively determine the wave field in the laminate.Sheng et al. [6,7] established a series of parallel viscoelastic layers with an elastic half space or a rigid foundation to discuss the characteristics of vibration generated by the load.Karlström and Boström [8] used Fourier series to build a layered half-space model of a linear viscoelastic layer, and Karlström [9] compared the train speed with the corresponding constant speed.In addition, many methods have been used to conduct the dynamic analysis of layered viscoelastic media, such as finite element method (FEM) (Mulungye et al. [10]; Kim and Tutumluer [11]; Younesian and Sadri [12]; and Connolly et al. [13]), moving element method (MEM) (Koh et al. [14]), and boundary element method (BEM) (Rasmussen et al. [15]).
Because of the complexity of road structure layers, when analysing the dynamic response of the road structure, the vehicle load is generally simplified as a static load or as a harmonic load, and the surface layer is simplified as a thin plate or as a beam on a viscoelastic subgrade (Winkler or Kelvin).Kim and Roesset [16], Kim and McCullough [17], and Kim [18] studied the dynamic response of an infinitely large Kirchhoff plate on a viscoelastic foundation under harmonic loads.e results showed that, when considering viscosity of foundation, there were significant differences from those obtained with an elastic system.Sun [19][20][21] employed integral transformation and Green's function to derive the dynamic displacements of beams and plates under moving loads.
e maximum displacement response was found to occur at the centre of the moving load, and the propagation speed was the same as that of the moving load.Kargarnovin and Younesian [22] and Kargarnovin et al. [23] investigated the dynamic response of an infinitely long Timoshenko beam on a foundation (Pasternak) subjected to a simple harmonic load.It was found that the load distribution shape, load velocity, and frequency significantly affected the dynamic response of the beam.Huang and ambiratnam [24] studied the dynamic response of a moving acceleration point load and discussed the influence of the initial moving speed, acceleration, and initial load position on the dynamic response.In another study [25], they further discussed the dynamic response of plate structures on elastic foundations under different moving loads.e results showed that the foundation stiffness and the speed and frequency of moving loads have a significant effect on the dynamic response and resonance speed of the plate.Lu et al. [26] and Lu and Yao [27] investigated the dynamic displacement response of a continuousreinforced concrete pavement (CRCP) on a viscoelastic subgrade under harmonic moving loads.Zibdeh and Abu-Hilal [28] discussed the random dynamic response of composite coated beams under random moving loads.It was found that the speed, blade thickness, and coating direction significantly affect the vibration response of the beam.
Most concrete tests have been conducted to describe the dynamic response of rail traffic under dynamic loads.Nimbalkar and Indraratna [29] improved the performance of geosynthetics and rubber shockmat as the buffer layer on different subgrades under the action of rail traffic.Liu et al. [30] studied the vibration reduction effect of a soilbag under dynamic loads.
e results showed that the vertical and horizontal vibrations can be effectively reduced using the soilbag.Shaer et al. [31] simulated the settlement and dynamic behavior of a ballasted railway track under high-speed train load.It was found that the stiffness of the track can undergo significant variations during the experiment, and the settlement function for the acceleration of sleepers was also analysed.However, the dynamic response characteristics of road filling materials under vehicle loading conditions have not been clearly understood.
In this study, an elastic double-layer plate model on a viscoelastic subgrade is used to analyse the dynamic response of a road under the action of dynamic loads.An asphalt layer and a stabilised layer (cement stabilised gravel with 10% content) are considered as elastic support plates.By simplifying the vibration differential equation used for the elastic double-layer plate and using the Fourier transform and trigonometric series, the vertical displacement expression for the double plate on the viscoelastic subgrade is obtained.In addition, the load diffusion path in the form of a cone in the stabilised layer under the conditions of single vibration source and three vibration sources is considered, and its effect on the stiffness of the structural layer is analysed.Moreover, the test data were used to verify the rationality and stability of other empirical formulae.

Elastic Double-Layer Plate Model on
Viscoelastic Subgrade e Kirchhoff plate assumption is stable and reliable only if the thickness-to-width ratio of a plate is lower than 0.05.Based on the single-layer plate model [32], an elastic doublelayer plate model on a viscoelastic subgrade is assumed, as shown in Figure 1. e thicknesses of the asphalt and stabilised layers are low compared to their length and width.erefore, the behavior can be described using an elastic double-layer plate with a width 2b � 0.8 m.An orthogonal X-Y-Z coordinate system is set up at the pavement layer.e load acting on the surface of the asphalt layer is F(x, y, t) and the circular load with d 1 � 0.15 m.
As shown in Figure 1, the established model consists of an upper plate to account for the asphalt layer (with mass density M 1 per unit area and flexural rigidity D 1 ) and a lower plate to account for the stabilised layer (with mass density M 2 per unit area and flexural rigidity D 2 ).e flexural rigidity D is given by Eh 3 /12 ( where According to Liu et al. [33], when y � ±b, the boundary conditions of the plate can be simplified to Given the assumption that the plate is infinitely long in the road direction, when x ⟶ ±∞, Following the approach used by Zhang et al. [32], the displacement and load in equation ( 1) can be expanded using trigonometric series as follows: where α n � nπ/2b.
Substituting equations ( 4) and (5) into equation (1), we obtain the following linear differential equations: e solution to equation ( 6) is found to be where we can obtain c 1 ⟶ 0 and c 2 ⟶ 0. e solution to equation ( 6) can be rewritten as follows: For the considered load F n and stiffness K, equation ( 8) can be shown in the next two sections: Case A and Case B.

2.1.
Case A: Single Jack.According to the established ballast vibration model [34], this paper studies a load diffusion model in the stabilised layer, as shown in Figure 2. Single jack (vibration source) is considered for the circular load with diameter d 1 � 0.15 m, as shown in Figure 2(a).e thickness of the asphalt layer is lower.erefore, the load transfer largely affects the stabilised layer and the subgrade.
e load F n and stiffness of the asphalt layer K 1 can be expressed, respectively, as follows: where P is the excitation load, t is the load time, h 2 is the thickness of the stabilised layer, ζ n is a constant, α is the load distribution angle [34], and E 1 is the elastic modulus of the Asphalt layer Advances in Civil Engineering asphalt layer.In addition, i −1, ω 2πf, and f is the excitation frequency.
For the load transmission of single jack, as shown in Figure 2(a), the sti ness of the stabilised layer K 2 can be expressed as follows: where E 2 is the elastic modulus of the stabilised layer.Furthermore, the sti ness of the subgrade K 3 under single vibration source can be de ned as follows: where E 3 is the elastic modulus of the lling medium in the subgrade.

Case B: ree Jacks.
As shown in Figure 2(b), the load di usion model of the three jacks (vibration sources) is based on the assumption of single jack.In the case of the thickness of the stabilised layer, if the spacing between the adjacent vibration sources is small enough, an overlapping of the adjacent load distribution may occur.e height of the overlapping regions is de ned as h 3 .e sti ness of the stabilised layer K 2 consists of two parts in series: where where l is the vibration source spacing and β arccos(l/d 1 + 2h 2 tan α), ∈ (0 °, 90 °).e sti ness of the subgrade K 3 and load F n under three vibration sources can be rewritten as follows:

Test Devices and Materials
Figure 3 shows the model structures and layout of the scaled test, including an asphalt layer, a stabilised layer, four displacement meters, ve three-way capacitive accelerometers, and ten soil pressure cells.e displacement meters are arranged at the centre of the vibration source and at the boundary of the plate.e ve accelerometers and the ve central soil pressure cells are set in the vertical direction of the vibration source with a vertical distance of 0.15 m. e adjacent vertical spacing between the slope soil pressure cells is 0.15 m, and their horizontal distances to the load centre are 0.3, 0.4, 0.5, 0.6, and 0.7 m.In Figure 3, the dynamic load is simulated by varying the number of hydraulic jacks, and the adjacent spacing is 0.175 m.
e dynamic behavior of the subgrade is studied by measuring the acceleration and effective soil pressure.e displacement of the elastic doublelayer plate on the viscoelastic subgrade is analysed using the e gaps between the sand bags are filled with the same material as that in the bags; this filler material is used to stabilise the test components, as shown in Figure 3. e main factors that affect the dynamic response include the subgrade filling material, vibration load, vibration frequency, and number of jacks.In this experiment, the sand bags and gravel are selected as the subgrade filling medium.e excitation load and excitation frequency are, respectively, set to 2000, 3500, and 5000 N and 1, 5, 10, and 15 Hz. e conditions of single ack and three jacks are considered separately in this work.Table 1 lists the relevant parameters of the filling medium in the subgrade, mainly including the resilient modulus, Poisson's ratio, damping constant, internal friction angle, height, and particle size.Table 2 lists some parameters of the asphalt and stabilised layers.Tables 3-5 list some parameters of the displacement meter, the three-way capacitive accelerometer, and the soil pressure cell, respectively.

Validation and Comparison
To study the dynamic response of the road, the subgrade filled with sand bag and gravel is used for testing in the experiment.e purpose of changing the number of jacks is to study the superposition effect of regions under multiple vibration sources.By adjusting the excitation load and frequency through the control system, the influence of different loads on the dynamic response of the road can be obtained.
4.1.Displacement.Figure 4 shows the variation in the displacement with respect to time for excitation frequencies (f ) of 1, 5, 10, and 15 Hz (P � 5000 N).Figures 4(a) and 4(b) show the displacements of the elastic double-layer plate under single jack and three jacks conditions, respectively.Figure 4 shows that there is a relatively good agreement between the theoretical calculation results and the test data.e theoretical calculation curve reflects the correspondence between the frequency and the displacement, and the theoretical formulae under the single and three jacks conditions are respectively expressed as follows: where a 0.15 denotes the peak acceleration under different excitation frequencies, which are listed in Table 6.

Advances in Civil Engineering
Figure 5 shows the variation in the displacement with respect to time under different excitation loads (2000 N, 3500 N, and 5000 N) for an excitation frequency of 10 Hz.
For comparison, the theoretical results are also plotted in the figures, and the theoretical formulae for the two cases can be respectively written as follows: As shown in Figure 5, it is worth mentioning that the results are largely in good agreement.For a given number of jacks, the displacement of the elastic double-layer plate increases with the increase in the excitation load.Furthermore, the displacement under three jacks is higher than that under single jack.To further verify the practicability of the theoretical derivation, another filling medium in the subgrade (gravel) is selected, as shown in Figure 6.Tables 2 and 3 list the related parameters.Overall, the elastic double-layer plate model has a relatively good agreement with the experiments.However, because of the limited variables selected for the experiments, the numerical results cannot be concretely verified.In the following experiments, we discuss the theory of elastic double-layer plate on a viscoelastic subgrade in a more comprehensive manner.
Here, we now compare the test data under three jacks with three times that subjected to single jack (P � 5000 N, f � 10 Hz).As shown in Figure 7(a), the maximum displacements in the case of the sand bag are 0.43 mm and 0.89 mm for the former and latter, respectively.In addition, it is found that the maximum displacement in the case of gravel is, respectively, equal to 0.42 mm and 0.83 mm, as shown in Figure 7(b).It is noted that the displacement of the elastic double-layer plate under the three jacks condition cannot be estimated simply by multiplying the parameters associated with the single jack by three times.
An empirical formula for calculating the cumulative deformation Z of the elastic double-layer plate can be expressed as follows [35]: where Z 0 is defined as the rebound displacement when the deformation gradually becomes stable (P � 5000 N, f � 10 Hz).
Figure 8 shows the variation in the rebound displacement in the time interval 80-110 s.It is found that the rebound displacement Z 0 under the single jack and three jacks conditions are 0.2 and 0.31 m and 0.17 and 0.28 m when the sand bag and gravel are used as the filling medium in the subgrade, respectively.In addition, c and d are two parameters.e following equation is an alternative empirical approach to calculate Z [29]: where a and b are two parameters.8 Advances in Civil Engineering Figure 9 presents the comparisons of the displacements obtained from the empirical formulae ( 17) and (18), equation (8), and test data (sand bag).e test data have a certain nonlinear relationship with the theoretical formula.Moreover, there is a strong consistency between the empirical formulae and the test data under both single jack and three jacks conditions.To further verify the relationship between the theoretical formula, empirical formulae, and test data, an alternative filling medium in the subgrade (gravel) is selected, as shown in Figure 10.
e results have a strong consistency.In other words, the theoretical formula for the elastic doublelayer plate is thoroughly verified via this experiment.Moreover, empirical formulae ( 17) and ( 18) were effectively applied.

Effective Soil Pressure.
Under the above loading conditions, the effective soil pressure in the subgrade medium δ C can be expressed as follows:  Advances in Civil Engineering where H is the height of the vibration source, h is the embedded depth of the sensor, S is the effective loading area on the top of the subgrade, n is the number of jacks, and p and q are two parameters.It is also assumed that the initial effective soil pressure δ 0 � 0. e embedded depth of the soil pressure cells, h, is a parameter found to be important enough to affect the dynamic response of the subgrade.erefore, the depth h with values of 0.15, 0.3, 0.45, 0.6, and 0.75 m is set during the test process.
e simulations are conducted in load cycles of 200 s.Figures 11(a) and 11(b) show the comparisons of the effective soil pressure with respect to different depths under the single jack and three jacks conditions, respectively.As expected, the effective soil pressure decreases with the increase in the depth h.It is found that the theoretical formula most closely matches the test data when h � 0.45 m. e theoretical formulae, shown in Figure 11, for the single jack and three jacks conditions are expressed as follows, respectively: Figure 12 shows the effective soil pressure under different excitation loads obtained from the theoretical analysis compared to the test data for h � 0.45 m and f � 15 Hz. e results obtained from the theoretical formula are in good agreement with the test data under different excitation loads.In addition, the results show that the effective soil pressure increases with an increase in the excitation load under both single jack and three jacks conditions.
Generally, the effective soil pressure under single jack is lower than that under three jacks.Undoubtedly, this is a result of "effective superposition effect" induced by the joint action of multiple loads, thereby increasing the dynamic response of the subgrade.However, the effective soil pressure under the three jacks condition cannot be estimated simply by multiplying the parameters associated with the single jack by three times.
Figure 13(a) shows the effective soil pressure under different excitation frequencies obtained from the theoretical analysis compared with the test data.It is found that the effective soil pressure increases with an increase in the excitation frequency in the case of the single jack.However, the results obtained from the theoretical formula are different from the test data.e relative errors between the theoretical maximum effective soil pressure and the test data are 7.44%, 14.1%, 3%, and 6.7% for f � 1, 5, 10, and 15 Hz, respectively.e theoretical formula is written as follows: e parameter q is very important for determining the effective soil pressure.As shown in Figure 13(b), the variations in the effective soil pressure with respect to time under f � 1, 5, 10, and 15 Hz are, respectively, plotted when q is equal to −0.017, −0.017, −0.019, and −0.02.e test data are found to closely match the theoretical curves.
Figures 14(a with sand bag and gravel, respectively.e gaps between the sand bags are filled with the soil.e figures show that the effective soil pressure increases significantly with the increase in the excitation load, particularly when three jacks are used.In addition, the effective soil pressure increases with the increase in the excitation frequency when the Advances in Civil Engineering 13 effective soil pressure reaches the maximum value for an excitation frequency of 15 Hz.Further, when the filling materials are changed from sand bag to gravel, the effective soil pressure generally increases.In other words, filling the subgrade with gravel results is a greater vibration response compared to filling with the sand bag.An empirical formula to obtain the effective soil pressure of the subgrade medium δ B can be expressed as follows [36]: where g and j are two parameters.An alternative empirical approach to calculate the effective soil pressure of the subgrade medium δ D can be written as follows [37]: where H 0 is the height of the subgrade medium; the length of the plate b 0 � l 0 � 0.8 m; k p is the effective soil pressure constant (k p � 1 + sin θ/1 − sin θ, where θ is the internal friction angle of the subgrade medium); and a 1 and a 2 are two parameters.Figures 15 and 16 compare the measured effective soil pressures with those obtained using the empirical formulae for a subgrade with different filling media in the case of single jack.Most of the test data not only have a certain nonlinear relationship with equation ( 22) but also have a certain linear relationship with equation (23).Moreover, the relative error between most of the measured data and the two empirical formulae is low.
To further verify the relationship between the test data and empirical formulae, the analysis is performed using three jacks, as shown in Figures 17 and 18. e results have a strong consistency.In other words, the empirical formulae ( 22) and ( 23) can be used to effectively describe the dynamic response of the subgrade in this experiment.
e effective soil pressures of the central and slope areas are generally tested under four conditions: sand bag-single jack, sand bag-three jacks, gravel-single jack, and gravelthree jacks, as shown in Figure 19 (P � 5000 N, f � 10 Hz).
e effective soil pressure and attenuation rates are compared under the four conditions.It is found that the attenuation rate of the same filling medium in the horizontal direction under the single jack condition is generally higher than that under the three jacks condition.Moreover, the effective soil pressure of the central area is generally higher than that of the slope area. is result is observed mainly because the load effect in the slope area is smaller than that in the central area.In addition, the load effect gradually weakens as the distance from the vibration source increases.e figures show that the peak acceleration increases significantly with the increase in the excitation load, particularly for the excitation when three jacks are used.In addition, there is an obvious increase in the peak acceleration with the increase in the excitation frequency when the peak acceleration reaches the maximum value for an excitation frequency of 15 Hz.Taking three jacks as an example, when the filling medium changes from sand bag to gravel, the peak acceleration a 0.15 under P � 5000 N and f � 10 Hz changes from 2.31 m/s 2 to 2.15 m/s 2 .In other words, the subgrade filled with gravel has a greater vibration reduction effect than that filled with sand bag.

Acceleration.
Taking the single jack as an example, Figure 22(a) shows that the peak acceleration of the subgrade has a certain linear relationship with the distance from the vibration source (P � 3500 N, f � 5 Hz).
e fitting linear function can be expressed as follows:  Advances in Civil Engineering where y 0 is a coefficient related to the excitation frequency, excitation load, subgrade filling medium, and the number of jacks; λ represents the attenuation rate; and h is the distance from the vibration source.Table 7 lists the values of y 0 and λ.
Although the subgrade filled with gravel has a greater vibration reduction effect than that filled with sand bag, the attenuation rate of the latter is apparently higher than that of the former.In addition, it can be seen from the relative error that there is a strong agreement between the fitting formula and the test data.
To further verify the above conclusions, the analysis is performed under the three jacks condition, as shown in Figure 22(b).e results show that equation ( 24) can be used to fit the acceleration under the three jacks condition.

Conclusions
e dynamic behavior of a road on subgrades filled with sand bag and gravel is investigated through theoretical derivation and a series of laboratory experiments.e effects of factors, such as the excitation load, excitation frequency, and number of jacks, on the settlement and vibration are analysed.Based on the findings of this study, the following conclusions can be drawn: (1) Based on the double-layer plate model on the viscoelastic subgrade, an indoor traffic load model system was used to study the dynamic responses of subgrade filled with sand bag and gravel.Gravel as a subgrade filling material was found to have a better vibration reduction effect than sand bag.(2) e effect of the transfer path on the rigidity of the structural layer was analysed, and the theoretical calculation results were in good agreement with the test data.In other words, the proposed load diffusion model was more accurately verified via experiments.(3) e dynamic response of the road increased with an increase in the excitation load and frequency.It is also worth noting that the dynamic response obtained under the three jacks condition cannot be estimated simply by multiplying the parameters associated with the single jack by three times.

Figure 1 :
Figure 1: Elastic double-layer plate model on a viscoelastic subgrade.

Figure 3 :
Figure 3: Model structures and layout of scaled test.

Figure 4 :
Figure 4: Displacement versus time curves under different excitation frequencies (sand bag): (a) single jack and (b) three jacks.

Figure 5 :Figure 6 :
Figure 5: Displacement versus time curves under different excitation loads (sand bag): (a) single jack and (b) three jacks.

Figure 7 :
Figure 7: Variation in the displacement versus time curves: (a) sand bag and (b) gravel.

Figure 11 :Figure 12 :
Figure 11: Variation in the effective soil pressure versus time curve at different measurement points (sand bag): (a) single jack and (b) three jacks.

Figure 13 :Figure 14 :
Figure 13: Variation in the effective soil pressure versus time curve under different excitation frequencies (sand bag).
Figures 20 and 21 compare  the measured peak accelerations under excitations in the cases of single jack and three jacks, respectively, for subgrades filled with sand bag and gravel.

Figure 15 :Figure 16 :
Figure 15: Test results and empirical formulae at different measurement points in the case of a single jack (sand bag).

Figure 18 :Figure 17 :
Figure 18: Test results and empirical formulae at different measurement points in the case of three jacks (gravel).

Figure 20 :
Figure 20: Peak acceleration for subgrade filled with sand bag under (a) single jack and (b) three jacks.

Figure 19 :Figure 21 :Figure 22 :
Figure 19: Effective soil pressure of the central and slope areas under different conditions.

Table 3 :
Parameters of the displacement meter.

Table 1 :
Parameters of the subgrade.

Table 2 :
Parameters of the elastic double-layer plate.

Table 4 :
Parameters of the three-way capacitive accelerometer.

Table 5 :
Parameters of the soil pressure cell.

Table 7 :
Values of the fitting parameters under different conditions shown in Figure22.