Abstract
Aiming at the seismic response of plastic geogrid-reinforced embankments, with Zhounan Expressway as the research engineering background, a self-designed seismic-rainfall coupled slope model test system was designed and used to produce 1 : 20 scale plastic geogrid-reinforced embankments. Moreover, the physical model of the unreinforced embankment under Hanshin wave, Wenchuan wave, Tianjin wave, etc. was also studied to carry out comparative analysis on seismic response and dynamic response on test model. The dynamic characteristics and dynamic response of the embankment model were tested from low to high seismic intensity; the changes of the embankment’s natural frequency, damping ratio, acceleration at the measuring point, and dynamic earth pressure were analyzed; and the main influencing factors and damage to the embankment seismic response feature were discussed herein. The test results showed that the initial natural frequency of the reinforced embankment was 42.4% higher than that of the unreinforced embankment, and its initial damping ratio reduced by 19.4%. The attenuation effect of the natural frequency and damping ratio of the reinforced embankment with the loading history was significantly lower than that of the unreinforced embankment. Embankment reinforcement exhibited a very good inhibitory effect on the PGA amplification effect of the embankment, and the inhibitory effect on the interior of the slope was more significant than that on the slope. Moreover, the type of seismic wave, the amplitude of the seismic wave, and the frequency of the seismic wave significantly influenced the PGA amplification effect of the embankment. The peak dynamic soil pressure of the unreinforced embankment at the same location was significantly greater than that of the reinforced embankment. The two embankment models showed significantly different antivibration damage performance. After the peak acceleration of 2 m s-2 was loaded, no cracks were seen on the surface of the embankment model. When the peak acceleration of 3 m s-2 was loaded, on the slopes of the two embankment models, smaller cracks were observed in the middle and upper parts of the face. When the peak acceleration of 4 m s-2 was loaded, the failure of the unreinforced embankment model was obvious. Large cracks on the top of the slope could reach 16 mm in width, and 27 mm settlement appeared at the top, and the slope was convex. The reinforced embankment model was only on the slope shoulder. Moreover, there were fine cracks on the top, and the slope top settlement was less than 5 mm. The research results provide theoretical support for preventing and controlling the road embankment vibration diseases and improving highway durability design.
1. Introduction
China is located between the Pacific Rim seismic belt and the Himalayan-Mediterranean seismic belt. Since the 20th century, strong earthquakes have occurred in China and the frequency of earthquakes has been increasing year by year [1]. There are frequent reports on structural damage to highway subgrades during previous earthquakes. Moreover, subgrade subsidence, cracking, distortion, and slippage occur frequently. The dynamic response characteristics of subgrades under earthquakes have an important impact on highway safety and durability.
Reinforced embankments have gradually been widely used in engineering practice. Reinforcement can limit the lateral deformation of the embankment to a certain extent, disperse the additional load and settlement of the embankment, inhibit the transmission of seismic loads, and finally improve the seismic resistance of expressways to varying degrees [2–4]. Extensive research attention has been paid to the investigation of the dynamic response and seismic performance of reinforced soil slopes. Although there are few researches directly highlighting the dynamic problems of reinforced embankment slopes, numerous research reports are available on the study of reinforced slopes and retaining walls under earthquake action. For instance, Jianzhou et al. [5] used PLAXIS software numerical calculation method to study the dynamic characteristics of double-sided reinforced embankments under earthquake action and analyzed the maximum stress distribution of each layer of the embankment and the settlement form of the embankment. Lesniewska [6] used RES analysis software and the rigid-plastic theory of reinforced soil to analyze the differential boundary value of the bearing capacity of reinforced and unreinforced soil embankment slopes. Halder et al. [7] combined the lower bound finite element limit analysis method, anisotropic random field modeling method, and Monte Carlo simulation analysis to calculate the bearing capacity of the cohesive soil-reinforced embankment slope. Nouri et al. [8] used the limit equilibrium horizontal slice method to evaluate the pseudostatic acceleration and amplification effect of reinforced soil slopes and retaining walls. Furthermore, Lin et al. [9] compared the ground motion response of an unreinforced slope with that of a reinforced slope and found that the vertical acceleration magnification rate of the reinforced embankment slope was much smaller than that of the unreinforced embankment slope. Lihua et al. [10] conducted shaking table model experiments to study the dynamic response performance of unreinforced slope, waste tire strings and three-way geogrid composite reinforced slope, and tire strings and tire fragment composite-reinforced embankment slopes under earthquake action. The acceleration responses of the embankment slope model along the slope height distribution law under different seismic waves, acceleration peaks, and reinforcement methods were discussed. Hongwei et al. [11] carried out shaking table model tests for investigating the seismic performance of geobag-reinforced soil retaining walls under earthquake action. Lu [12] conducted a shaking table test on composite Gabion geogrid-reinforced soil retaining wall model with different similarity ratios, describing the macroscopic phenomenon of the test model under earthquake action, and they further studied the model retaining wall under earthquake action. Using the dynamic response, they tested the acceleration response and magnification of the backfill and tested the frequency spectrum characteristics at different heights of the model, the horizontal displacement of the wall, and the vertical settlement of the fill and finally compared and analyzed the parameters given in the design code. Jiang et al. [13] developed a combined panel reinforced soil retaining wall and conducted on-site filling tests and indoor shaking table model tests. Srilatha et al. [14] studied the influence of the fundamental vibration frequency on the dynamic response of unreinforced and reinforced soil slopes through shaking table tests and found that the acceleration and displacement responses of slopes with different reinforcement numbers and positions did not increase with frequency and linear growth. Bahadori et al. [15] conducted a shaking table test study on the liquefaction problem of the liquefiable soil layer reinforced with geogrid and geocomposite and found that the antiliquefaction settlement performance of geocomposite reinforcement was significantly better than that of geogrid reinforcement. Panah et al. [16] conducted a series of 1 g shaking table tests on the 80 cm high reinforced soil retaining wall model and studied the effect of the length, distribution, and shape (sawtooth and parallel) of the steel bar on the failure mode of the retaining wall. Furthermore, displacement and acceleration magnification factor has also been studied. Wartman et al. [17] analyzed the mechanism of permanent displacement of the slope under the action of seismic load and commented on the Newmark slider displacement calculation method. Koseki et al. [18] analyzed the seismic performance and potential failure mechanism of reinforced earth retaining walls. Edinçliler et al. [19] used the shaking table test method to study the seismic performance of slag as a filler for reinforced soil retaining walls.
Plastic reinforced soil embankments are widely used in engineering practice. Although scientific and technological researchers have carried out certain research work on the seismic effects of reinforced soil slopes, the seismic effects of plastic geogrid-reinforced embankments have rarely been investigated [20–25]. The current seismic design of highways nor does the code evaluate the seismic performance of plastic-reinforced embankments. In-depth study of the dynamic response of plastic-reinforced embankment under earthquake action and its failure characteristics was used to qualitatively and quantitatively analyze the seismic performance and reinforcement effect of plastic-reinforced embankment. This study provides important theories and engineering design/guidelines for preventing and treating earthquake-induced damage and destruction of embankment.
2. Model Test Design
2.1. Model Test System
In this study, a self-designed and customized seismic action slope model test system (see Figure 1), including hydraulic servo power loading system, shaking table, model box, control system, data acquisition, and analysis and processing system, was adopted to perform the tests. The dimension of the slope model box was , which can be used for various earthquake-induced slope scale model tests.

(a) Hydraulic power system

(b) Test model box

(c) Countertop
2.2. Engineering Background
The embankment model test considered the No. 8 section of Zhounan Expressway as the research engineering background, and the pile number section was K18+000–K42+600 and K42+600–K54+100. The upper part of the strata is dominated by Quaternary Holocene alluvial silt clay and silt soil, with some interbedded sand layers, and the upper part of the silty clay is soft plastic-plastic; and the lower part is an interbedded silty clay and silt soil. The sand layer is plastic-hard plastic. The earthquake fortification intensity is 7°. This is a section with a large embankment filling height, the maximum filling height of 10 m, the soil cohesion of the foundation bearing layer KPa, the internal friction angle , the natural density of 2074 kg m-3, and the embankment earth pressure. The actual control dry density is 1790 kg m-3, the cohesion force KPa, and the internal friction angle φ =24°.
2.3. Design and Production of Test Model
The embankment model design and the layout of the main measuring points are shown in Figure 2. Two types of embankment models were constructed: one was a filled embankment without reinforcement measures, and the other was a plastic-reinforced embankment. Owing to the symmetry of the cross section of the embankment, the two embankment models were designed in half width, and the geometric similarity ratio Cl was designed to be 20. Therefore, the total height of the prototype embankment after being scaled down was 800 mm, the cross-sectional direction of the foundation was 2000 mm, the thickness range was designed to be 300 mm, the height of the filled embankment was 500 mm, and the top cross-sectional dimension was 675 mm. The size of the foundation and the embankment along the road was 1500 mm according to the model box size. The slope angle of the embankment was the same as that of the prototype; using a broken line grading, the upper part was grading 1 : 1.5; and the lower part was grading 1 : 1.75. The ingredients and filling methods of the two embankment models were completely the same. The plastic reinforcement laying method of the reinforced embankment was horizontal through long bars, and the vertical spacing was 100 mm.

(a) Unreinforced embankment model design

(b) Reinforced embankment model design
The model similarity design takes the embankment geometry size and packing density as the basic control quantities and determines the similarity relationship of other physical parameters of the embankment according to the Buckingham theorem to calculate the similarity law of the embankment model. The physical model test similarity design is presented in Table 1.
After the embankment model was designed, the soil was sieved to remove oversized, layered, and compacted particles. The foundation soil was evenly filled and compacted in six layers, the filled embankment was evenly filled and compacted in five layers, and the plastic geogrid was spread between each layer of the reinforced embankment model. Each layer of soil was compacted to control the dry density of 1790 kg m-3. Corresponding sensors were embedded in the model making project based on the designed measuring points. The completed embankment model entities are shown in Figure 3.

(a) Unreinforced embankment model

(b) Reinforced embankment model
3. Seismic Wave Selection and Test Loading System
Three typical seismic wave types, namely, Hanshin wave (code-named KOB) Figure 4, Wenchuan wave (code-named WC) Figure 5, and Tianjin wave (code-named TJ) Figure 6, were used for test loading. The peak acceleration of various waves was controlled to be at 0.5, 1, 2, 3, and 4 m s-2 five-level loading, and various waves were processed with four types of time compression ratios. The peak acceleration of white noise (code named WTN) was controlled to be at 0.3 m s-2, as shown in Figure 7, to measure the initial dynamic characteristics of the model and its changes during the loading process. The loading of various seismic waves was carried out interspersed by magnitude, and the specific loading system is presented in Table 2.




4. Test Results and Analysis
Comparative analysis of the dynamic characteristics was carried out. Dynamic response and influencing factors of reinforced and unreinforced embankment models were studied, and the failure characteristics of the two embankment models were investigated.
4.1. Dynamic Characteristics of Embankment Model
The natural frequency and damping ratio are the main dynamic characteristic parameters of the embankment model. Before seismic waves were loaded in each working condition, the embankment model was tested with white noise with an amplitude of 0.3 m s-2, and the white noise transfer function of each measurement point was identified. Moreover, the average value of the identified natural frequency and damping ratio was taken as the characteristic natural frequency and damping ratio of the embankment model. The curves highlighting the variation in the natural frequency and damping ratio of the two embankment models with the loading history are shown in Figure 8.

(a)

(b)
Figure 8(a) demonstrates that the natural frequencies of the two embankment models gradually decrease with the loading process. The initial natural frequency of the unreinforced embankment is 12.63 Hz, until the seismic wave test conditions are loaded, and then its natural frequency is reduced to 9.23 Hz, corresponding to a drop of 26.9%. When the initial natural frequency of the reinforced embankment is 17.98 Hz, the load is completed. The frequency is reduced to 13.74 Hz, corresponding to a drop of 23.6%. Figure 8(b) illustrates that the damping ratio of the two embankment models gradually increases with the loading history. The initial damping ratio of the unreinforced embankment was 0.098, and the damping ratio increased to 0.146; an increase of 49.4% after the seismic wave test conditions was loaded. The initial damping ratio of the reinforced embankment was 0.079, and the damping ratio increased to 0.106 after the loading was completed. This corresponded to an increase of 37.3%. The test data showed that reinforcement could effectively improve the dynamic characteristics of the embankment. The initial natural frequency of the reinforced embankment was 42.4% higher than that of the unreinforced embankment, and its initial damping was reduced by 19.4% compared to the unreinforced embankment. After the multicondition loading process, the natural frequency decreased and the increase of damping ratio of the reinforced embankment model was smaller than those of the unreinforced embankment, Moreover, the dynamic damage resistance of the reinforced embankment was significantly better than that of the unreinforced embankment.
4.2. Dynamic Response of Embankment Model
The PGA amplification effect and dynamic earth pressure of the measuring points under different loading conditions were used as the main observation indexes of the dynamic response, and the analysis is presented in the subsequent sections.
4.2.1. Analysis of PGA Amplification Effect at Measuring Point
The PGA amplification factor of the measuring point was calculated according to the peak acceleration data of the measuring point in each working condition, that is, the ratio of the peak acceleration of each measuring point to the measured peak acceleration of the table. After comparing the data of multiple working conditions, the PGA amplification coefficients of the three seismic wave action measuring points exhibited similar characteristics. Herein, the KOB-18 and KOB-5 working conditions were used for analysis. The PGA amplification coefficients of the two embankment model measuring points were along the height. The change rule is shown in Figures 9 and 10.


The figure shows that from bottom to top along the height of the embankment, the PGA amplification coefficients of the slope measurement points of the two embankment models show an overall increasing trend, with the maximum value appearing at the top of the slope. The PGA amplification coefficients of the reinforced embankment at the same height are less than those of the unreinforced embankment. The analysis of two types of embankment model slope body measuring points PGA amplification factor indicates that the lower and upper quarter height range of the slope body shows a small increase. The upper part even shows a decrease in a certain height range, and the middle part exhibits the fastest increase. Moreover, no increase is observed at the same height of the reinforced embankment. The PGA amplification factor of the A11 measurement point in the slope was greater than that of the slope A4; however, the reinforced embankment did not exhibit this phenomenon, indicating that the reinforcement exerted an inhibitory effect on the PGA amplification effect of the embankment, in particular, the internal slope.
Based on the test data, the influences of seismic wave type, seismic wave amplitude, and seismic wave frequency on the PGA amplification factor are discussed below.
(1) Influence of Seismic Wave Type. In order to discuss the impact of seismic wave types on the acceleration response of the two embankment models, the KOB-10, WC-10, and TJ-10 working conditions were selected for analysis. The amplitude for all the three seismic waves was 2 m s-2, and the time compression ratio was 4.472. Figure 11, respectively, shows the same amplitude and time compression ratio, the PGA amplification coefficients of two embankment slopes A1, A3, A5, A7 and A8, A10, and A12 on the slope of the slope.

(a)

(b)
The figure illustrates that at the same amplitude and time compression ratio, the PGA amplification factor of the embankment measurement point under the action of the WC wave is the largest, followed by that of the KOB wave, and that of the TJ wave is the smallest. Moreover, PGA amplification effects of other points exhibit the same rule as the given measuring points. It shows that when the types of seismic waves are different, the PGA amplification coefficients of the measurement points at the same location of the embankment are also different. Furthermore, the PGA amplification factor of reinforced embankments is significantly lower than that of unreinforced embankments.
(2) Influence of Seismic Wave Amplitude. In order to discuss the impact of seismic wave amplitude on the acceleration response of the two embankment models, the Wenchuan waves WC-4, WC-8, WC-12, WC-16, and WC-20 with different amplitude operating conditions were selected for analysis. Moreover, the above-mentioned operating conditions show the same time compression ratio and seismic wave type; however, the acting amplitude is 0.5, 1, 2, 3, and 4 m s-2, respectively. Figure 12 shows the PGA amplification coefficients of the two embankment models and the variation curves of seismic wave amplitude.

(a)

(b)
The figure shows that the PGA amplification coefficients of the two embankment model measuring points generally decrease with the increase of the seismic wave amplitude, and the PGA amplification coefficient of the slope top A7 measuring point shows the largest change. The two embankment seismic wave amplitudes are 0.5 m s-2, and the PGA amplification factor of the A7 measuring point can reach 2.2 times that of the amplitude 4 m s-2. The influence of KOB seismic wave and TJ seismic wave amplitude on PGA amplification coefficient is similar to that of WC seismic wave.
(3) Seismic Wave Frequency. In order to discuss the influence of seismic wave frequency on the acceleration response of the two embankment models, the Hanshin waves KOB-5, KOB-6, KOB-7, and KOB-8 were selected for the analysis of different frequency conditions, and the seismic wave amplitude of this condition was 1 m s-2. The time compression ratios are 1, 201/4, 201/2, and 203/4, respectively. Figure 13, respectively, shows the PGA amplification coefficient of the two embankment measurement points versus frequency.

(a)

(b)
The figure shows that the PGA amplification factor of the unreinforced embankment at the same input seismic wave frequency is greater than that of the reinforced embankment. Moreover, the PGA amplification factor of the two embankment measurement points changes with the input seismic wave frequency. The maximum value of the PGA amplification factor is on the upper part of the slope. For the measurement points, the maximum value of the PGA amplification factor of the measurement point of the unreinforced embankment appears in the range of the smaller time compression ratio, while the maximum value of the reinforced embankment is in the medium time compression ratio range. The influence of WC seismic wave and TJ seismic wave frequency on PGA amplification coefficient is similar to that of KOB seismic wave. It indicates that the resonance frequency of reinforced embankments is significantly higher than that of unreinforced embankments, and the frequency range of natural seismic action is closer to the natural frequency of unreinforced embankments. Therefore, reinforced embankments can effectively improve the seismic behavior of embankments.
4.2.2. Analysis of Measuring Point Dynamic Earth Pressure
The dynamic earth pressure in the embankment under the action of earthquake is the main controlling factor to the damage and destruction of the embankment. In order to discuss the change law of the dynamic earth pressure along the height of the two embankment models, the KOB-5 and KOB-18 loading conditions were selected for analysis. The distribution curves of the peak dynamic earth pressure at the measuring point are shown in Figure 14.

The figure exhibits that the peak dynamic soil pressure at the measurement point of the unreinforced embankment at the same height is significantly greater than that of the reinforced embankment. For the KOB-5 and KOB-18 working conditions, the peak dynamic soil pressure at the Y5 measurement point of the unreinforced embankment is 1.7 times that of the reinforced embankment, respectively, and then 2.5 times. Clearly, the reinforcement can effectively reduce the peak dynamic earth pressure in the embankment. Moreover, the graph and curve reveal that reinforcement measures can effectively adjust the peak dynamic earth pressure distribution in the embankment, greatly reduce the dynamic earth pressure in the upper part of the embankment, and significantly improve the dynamic stability of the embankment.
4.3. Failure Characteristics of Embankment Model
During the test, the dynamic characteristics and dynamic response parameters of the embankment were collected and analyzed. At the same time, the deformation and crack failure characteristics of the embankment were observed and measured in detail. After the peak acceleration loading of 2 m s-2 was completed, no cracks were observed in the two embankment models. After the completion of the peak acceleration loading of 3 m s-2, fine cracks appeared on the upper and middle slopes of the two embankment models. When the peak acceleration of 4 m s-2 was loaded, the cracks in the unreinforced embankment model became wider. The cracks were mainly distributed in the middle, and the upper part of the embankment and the slope shoulder. The width of the larger cracks reached 16 mm, and some cracks could be seen deep into the third part of the slope. The top of the embankment model subsided to 27 mm, and the slope was convex. After the reinforced embankment model was loaded, only two small cracks appeared on the shoulder and top of the slope, and the slope top settlement was less than 5 mm. The damage and destruction patterns of the two embankment models are shown in Figure 15.

(a) Unreinforced embankment model top crack

(b) Unreinforced road embankment model shoulder crack

(c) Fine cracks on slope top of reinforced embankment model

(d) Fine cracks on the slope of the reinforced embankment model
5. Conclusions
Aiming at the seismic response of plastic geogrid-reinforced embankment, with Zhounan Expressway as the research engineering background, a comparative model test of plastic geogrid-reinforced embankment and unreinforced embankment was designed and further conducted, and the embankment model was tested in various earthquake scenarios. The dynamic characteristics and dynamic response under loading conditions were analyzed. The change laws of the embankment’s natural frequency, damping ratio, acceleration at the measuring point, and dynamic soil pressure were analyzed, and the main influencing factors and damage characteristics of the embankment seismic response were discussed. Through the systematic analysis of the test data, the following conclusions can be obtained: (1)The dynamic characteristics such as the initial natural frequency and damping ratio of the plastic-reinforced embankment were found to be significantly better than those of the unreinforced embankment. The reinforcement led to the increase in the initial natural frequency of the embankment by 42.4%, and the initial damping ratio decreased by 19.4%. The reinforcement reduced the embankment natural frequency by 3.3% with the loading history, and the damping ratio decreased by 12.1% with the loading history(2)Plastic reinforcement exhibited a good suppression effect on the PGA amplification effect of the embankment; in particular, the suppression effect on the interior of the slope was more significant than that on the slope. The type, amplitude, and frequency of seismic waves exhibited an important influence on the PGA amplification effect of the embankment(3)Plastic reinforcement significantly reduced the peak dynamic earth pressure of the embankment at the same location. The analysis of representative working conditions and measurement point data indicated that the peak dynamic earth pressure could be reduced to 40% of the unreinforced embankment. Reinforcement could not only effectively reduce the peak dynamic earth pressure but also effectively adjust the peak dynamic earth pressure distribution in the embankment, greatly reduce the dynamic earth pressure in the upper part of the embankment, and significantly improve the dynamic stability of the embankment(4)Reinforcement could effectively improve the antivibration damage performance of the embankment. After loading with a peak acceleration of 3 m s-2, relatively small cracks appeared in the upper part of the slope of the two embankment models. However, when the peak acceleration of 4 m s-2 was loaded, the cracks of the unreinforced embankment model evolved into cracks, and the cracks were mainly distributed in the embankment. In the middle and upper part and the slope shoulder, the width of the larger cracks reached 16 mm, and some cracks penetrated one-third of the height of the slope, and there was serious settlement at the top. The maximum settlement at the top of the embankment model was 27 mm, and the slope was convex. However, the reinforced embankment model only showed two small cracks on the slope shoulder and top, and the slope top settlement was less than 5 mm. Reinforcement can significantly improve the strong earthquake resistance of the embankment
Data Availability
The data are available and explained in this article; readers can access the data supporting the conclusions of this study.
Disclosure
I would like to declare on behalf of my coauthors that the work described is original research and has not been previously published.
Conflicts of Interest
The authors declare no conflict of interest. The manuscript is approved by all authors for publication.
Acknowledgments
This work was supported by the National Key Research & Development plan of China (2016YFC0600901).