Dynamic Failure Mode and Dynamic Response of High Slope Using Shaking Table Test

A shaking table test was performed to study the dynamic response and failure modes of high slope. Test results show that PGA amplification coefficients increased with increasing elevation and the PGA amplification coefficient of the concave slope was slightly larger than that of the convex slope.-e slope type affected the dynamic response of the slope.-e elevation amplification effect of the concave slope under seismic load was more significant than that of the convex slope; thus, the concave slope was more unstable than the convex slope. Additionally, the PGA amplification coefficient measured on the slope surface was always larger than that inside the slope, and the data show an increasing trend with the broken line.-e dynamic amplification effect of the high slope was closely related to the natural frequency of the slope. Within a certain range, the higher the frequency, the more significant the amplification effect. -e dynamic failure process of concave and convex slopes was studied through tests. Findings indicate that the dynamic failure modes of the concave slope are characterized by shoulder collapse, formation of the sliding surface, and integral sliding above the slope line. Dynamic failure modes of the convex slope are mainly slips in the soil layer and collapse of the slope near the slope line.


Introduction
Southwest China is an earthquake-prone area, especially in Yunnan Province, a province with high earthquake occurrence.Red beds, a rock mass with poor engineering characteristics such as low strength, poor permeability, and easy weathering, are widely distributed in Yunnan.e area hosts typical landslide-prone strata.Earthquake landslides have become the most common geological hazard in the construction of these areas [1][2][3][4].Earthquake landslides are extremely dangerous and have caused severe casualties and property losses worldwide.In 2008, an 8.0-magnitude earthquake struck Wenchuan County, Sichuan Province, China, triggering massive slope damage, including more than 56,000 landslides resulting in 256 barrier lakes and nearly 20,000 deaths [5][6][7][8].e Yushu earthquake in China in 2010 triggered more than 2,000 landslides with a total area of 1.194 square kilometers [9].On July 22, 2013, a magnitude 5.9 earthquake occurred in Minxian, China, which also triggered more than 2,000 slope damages [10].erefore, slope stability has become a prominent technical problem in the construction of southwest Yunnan.
e dynamic stability of slopes is an important topic in geotechnical engineering [11].e study of slope dynamic response law is the foundation of slope dynamic stability analysis, including responses of acceleration, velocity, displacement, and stress-strain equivalents [12,13].Accurate understanding of the dynamic response characteristics of slopes under seismic loads can provide scientific and theoretical guidance for seismic designs in slope engineering.
Methods for studying seismic responses of slopes mainly include numerical simulation and model testing.e numerical calculation method can simulate the actual failure mechanism of landslides [14,15].Toki et al. used the finite element method to conduct slope dynamic stability analysis [16].Jiang et al. employed numerical methods to study the dynamic failure process of the model [17].However, some simple numerical values or correction coefficients cannot be used to accurately simulate actual situations due to the complexity of projects, leading to deviations in experimental conclusions.In recent years, researchers have begun to use the shaking table model test to study dynamic slope response.is test system can reproduce earthquake excitation and collect test data automatically and accurately.erefore, the shaking table test is a direct method used to reveal structural dynamic response and failure modes; the approach has been widely used in seismic response research on various structures [18][19][20].
Clough and Pirtz carried out a preliminary shaking table test study on the stability of slopes under earthquake conditions [21].Wang and Zhang used a vibration simulation test to study the dynamic stability of a rock slope [22].Wang and Lin used a shaking table test to analyze the initiation and displacement of landslides under earthquakes [23].Yang et al. investigated the disaster phenomenon of the Wenchuan earthquake landslide based on a shaking table test [24].Zhao et al. examined the mechanisms of deformation and instability of seismic slopes via a shaking table test [25].Matsuo used a large shaking table to study the dynamic response and failure modes of a geosynthetics-reinforced earth retaining wall [26].rough a large-scale shaking table test, Xu et al. studied the dynamic response law of a slope [27].Yang et al. investigated the dynamic failure mechanism and dynamic response characteristics of a counter-bedding slope under earthquake excitation [28].Che et al. used a shaking table test to analyze the influence of discontinuous joints on a high and steep slope in an earthquake [29].Lee et al. used shaking table model tests to study seismic dynamic response characteristics of dry sand and saturated sedimentary sand; they also discussed the boundary effect of the layered model box [30].Based on energy dissipation theory, Kokusho and Ishizawa studied the failure mechanism of slopes under seismic loads using a small-scale shaking table test, and the Newmark method was employed to calculate and analyze the failure mechanism [31].Ling et al. introduced a high-rise concrete model with a 2.7-meter-thick transfer slab into seismic action representing small, medium, and large earthquakes [32].Wartman et al. carried out a physical model experiment with a 1 g shaking table [33].Lin and Wang studied the dynamic response of slopes excited by sinusoidal waves (loading frequency of 8.9 Hz) using homogeneous sandy soil (water content of 5.1%) model slopes.Results showed that the slopes exhibited elastic characteristics when the acceleration amplitude was less than 0.4 g, and the amplification effect of slopes was significant as the nonlinear characteristics of soil became stronger [34].
Based on a field investigation of typical high slopes in southwestern Yunnan, we summarized the failure mechanism of model slopes via a shaking table test.On the basis of measured data, the effects of landform, slope type, elevation, and ground motion parameters on the dynamic response of the slope were analyzed in detail.

Study Area
Dabao and Baolong expressways are located in the red bed area of Yunnan province (Figure 1).e lithology of the red bed soft rock area is mainly composed of mudstone and sandstone; some rocks are complex lithologic assemblages formed by the interbedding of hard and soft rocks.Basic characteristics of slope height, slope morphology, and lithology combination were obtained through field investigation.
An investigation of 97 slopes along the Dabao and Baolong expressways revealed many high rock slopes, 80% of which exceeded 20 m and many of which were irregularly shaped.ere are a few types of convex and concave slopes: the upper part of the concave slope is steep, and the lower part is gentle, while the slope is protected by anchor frame beams; the upper part of the convex slope is gentle and the lower part is almost vertical, and the rock surface is weathered.
Figure 2 presents a common lithology combination diagram of soft and hard interbedded slopes.e upper part and surface of the rock mass are weathered and exhibit a vegetation layer, and two muddy interlayers exist inside the slope.ese muddy interlayers formed under various geological tectonic action and long-term groundwater action.
e weak interlayer plays a vital role in controlling slope stability.3) was adopted in the experiment.e platform size measured 2.2 × 2 m with a maximum load capacity of 40 tons.

Experimental Method
e maximum acceleration at full load was ±27.7 m•s −2 , and the working frequency range was 0.01-30 Hz; specific parameters are listed in Table 1.
A rigid model box with an internal size of 1.84 m × 1.3 m × 1.1 m was adopted in this test, welded from a steel plate and angle steel.e rigid model box was constructed of plexiglass on both sides to allow for observation of the slope failure process during vibration.e model box was lined with a 5 cm thick foam cushion on the back wall to reduce the boundary effect of the box. is method to reduce the model box effect has been adopted by many scholars [35,36].A wooden baffle was fixed in the center of the model box with angle steel and bolts and used to separate the two models.
e baffle was made of waterproof and smooth wooden materials.

Test Model.
e experimental model did not correspond fully to the characteristics of the high rock slope in southwestern Yunnan; rather, the model was simplified and manufactured according to common characteristics of slopes in this area.Two slope models (concave and convex) were included in this experiment.e scale of each model was 1 : 30; the height of each model was 1 m, the width of the top of the slope was 0.4 m, and the width of the bottom of the slope was 1.2 m.Both models adopted 1 : 0.5 and 1 : 1 slope rates.Taking barite powder, quartz sand, and gypsum as Figure 4 shows the completed models and the key steps of model construction.
Model construction is inherently complicated, and every step should be precise.e model-making process in this paper involved the following steps: (1) Proportionally mixed model materials were packed into the model box in batches, with 10 cm as a layer and then filled and compacted in layers.e boundary of each compacted layer was marked with vermilion-colored sand so model deformation and settlement could be observed through plexiglass on the side of the model box.
(2) When compacted to the interface of the designed weak interlayer, a layer of clay with a thickness of 2 cm was laid on the compacted surface, and the side was marked with blue sand.According to the design, four layers of clay were spread during compaction.(3) In the model-making process, the sensor was buried in the design position and small red plastic particles were sprayed on the upper surface; thus, the sensor could be easily found without damage when removing the model.(4) e slope surface was cut to make the shape coincide with the designed slope shape.

Similarity Relation and Similar Materials.
Based on π theorem [37][38][39], the similarity criterion was derived through dimensional analysis.Because of the complexity of the model, it would be difficult to meet similarity of all parameters.It was thus necessary to select the main control parameters as fundamental dimensions according to the test purpose and model characteristics.is test took geometric length, density, and acceleration as fundamental dimensions; similar constants of the test appear in Table 2.

Shock and Vibration
Although many materials could be used for modeling, we could not identify an ideally similar material.erefore, it is essential to understand the properties of test materials.Similar materials should meet the requirements of low price, stable mechanical properties, easy formation, simple measurement, uniformity, and isotropy.Similar materials generally consist of aggregates, cementing agents, and auxiliary materials in certain proportions.Barite powder    Shock and Vibration and quartz sand were selected in this study as coarse aggregates along with gypsum as the cementing agent and water as an additive.Clay was chosen to simulate the interface of the weak interlayer.e physical and chemical properties of aggregates were stable, and the barite powder was of great weight, thus meeting material requirements discussed above.Properties of similar materials are shown in Table 3.
To obtain an appropriate material ratio, several indoor geotechnical tests were conducted to acquire the physical and mechanical parameters of mixtures with different ratios.
e density of the model material was measured using a ring knife test: the model material was cut with a certain volume ring knife and weighed.e density of the material was the ratio of mass to volume of the material in the ring knife.e water content of the material was determined by a drying method: model material of a certain quality was weighed and placed into a drying oven and maintained at 105-110 degrees Celsius until the material reached a constant weight.
e ratio of the lost water quality to material quality after drying represented the water content.e cohesion and internal friction angle of materials were measured using a consolidated undrained shear test; the cylindrical specimen was coated with latex film and placed in a sealed pressure chamber.e material was equalized by injection of hydraulic pressure.At the same time, the drain valve was opened to consolidate the specimen fully, after which the drain valve was closed to apply axial pressure to destroy the specimen.Experimental results were plotted as σ − τ curves, and values of c and ϕ were obtained.e shear modulus of materials with different ratios was measured by a dynamic triaxial test.By comparing the physical and mechanical indices of materials in different proportions, the proportioning scheme closest to the parameters of similar models was found.Finally, a 32 : 56 : 9 : 8 ratio of barite powder: quartz sand: gypsum powder: water completed the slope model.Table 4 shows the main physicomechanical parameters of mudstone and similar materials; the indoor geotechnical test is depicted in Figure 5.

Measurement System and Input Motion.
A total of 18 acceleration sensors and 4 displacement sensors were used in the experiment.Accelerometers A1 to A16 were symmetrically distributed in two models as shown in Figure 6.Displacement sensors S1 and S2 were attached on the top of the slope to compare with the displacement obtained by integrating the acceleration time-history curve.Sensors A0 and A0′ and sensors S0 and S0′ were symmetrically arranged in the model; accelerometer A0′ and displacement sensor S0′ were installed on the top of the model box to monitor the dynamic response of the box.Acceleration sensor A0 and displacement sensor S0 were arranged on the shaking table surface.Acceleration sensors A15 and A16 were embedded in concave and convex slopes, respectively, 10 cm away from the box wall and 30 cm away from the bottom of the box; these sensors were used to detect the effect of the flexible boundary of the model box.Sensors A1, A7, and A13 were arranged on the same vertical line of the concave slope, 20 cm away from the box wall.Sensors A3, A9, and A14 were installed on corresponding points of the convex slope.Sensors A2, A5, A8, and A11 were arranged at elevations of 20 cm, 40 cm, 60 cm, and 80 cm, respectively, in the convex slope.Sensors A4, A6, A10, and A12 were installed at symmetrical measuring points on the convex slope.All sensors were waterproof.Figure 6 shows the full dimensions of the model slope and layout of the acceleration monitoring points.
For the loading scheme, elucidated in Table 5, the test adopted two wave patterns of an El Centro wave and sine wave with a time compression ratio of 5.48.e sine wave adopted four frequency grades: 5 Hz, 7.5 Hz, 10 Hz, and 12.5 Hz.
e peak acceleration value in the test loading scheme was determined according to the corresponding relationship between seismic fortification intensity and basic acceleration of ground motion in the Code for Seismic Design of Buildings (GB50011-2010).
e seismic fortification intensity in southwestern Yunnan is VII-IX.Peak acceleration levels of 0.1 g, 0.2 g, 0.3 g, and 0.4 g were designed corresponding to the fortification intensity, and the loading conditions of 0.6 g, 0.8 g, and 1.0 g were added.
Considering that the sine wave was more destructive, only

48
Note.L, ρ, and a are controlling parameters.* Fundamental dimensions.Shock and Vibration 0.1 g and 0.3 g levels with different frequencies were adopted in the test.White noise is an important means of measuring dynamic characteristics of a structure.To understand the dynamic characteristics of the model before and after excitation, white noise sweeping was applied at different stages.Model deformation was observed and recorded 15 minutes after every condition was implemented.

Dynamic Characteristics of Model Slope
Dynamic characteristics of the slope model included natural frequency and damping ratio.Dynamic characteristics were determined through the shaking table test, and variations in dynamic characteristics with the input seismic parameters of the platform were studied.Before the experiment was conducted, dynamic characteristics of the model box and treated boundary conditions were checked to prevent the "model box effect" from influencing the test results.

Natural Frequency Analysis of Model Box.
e rigid model box used in the experiment had its own natural frequency.When the natural frequency of the model box was close to that of the input seismic wave, a resonance phenomenon occurred.e model box produced large vibration, which destroyed the stability of the model slope and affected the test accuracy.6

Shock and Vibration
To avoid resonance between the rigid model box and the loaded seismic wave, a natural frequency experiment was carried out on the model box before the test.e natural frequency of the model box was 43.467 Hz, determined by a hammer test (Figure 7), which was much larger than the frequencies of seismic waves designed in the loading scheme and the natural frequencies of the model slope; thus, further testing avoided the resonance phenomenon.

Boundary Effect Analysis of Model Box.
As a semi-infinite body, a slope has no boundary in theory.In a shaking table test, however, the model can only be placed in a container of limited size.Because seismic waves propagate in soil, a series of reflections and refractions occur on the wall of the box along with changes in the vibration form of the system, which will introduce certain errors to the test, namely, the "model box effect."Effective simulation of soil boundary conditions should cause the model in the container to deform similarly to the original model in the free field under an earthquake to reduce the influence of boundary conditions.For this reason, the boundary effect of the model box wall was treated before the test; the specific method was described above.
e boundary effect was tested using special measuring points.Accelerometers A2 and A4 were compared with A0 to analyze the friction boundary effect; accelerometers A15 and A16 were compared with A1 and A3 to analyze the flexible boundary effect.As shown in Figure 8, peak acceleration at the A2 and A4 measurement points was significantly higher than that at A0. Waveforms at the A2 and A4 measured points did not lag obviously behind that of A0, indicating that the bottom of the model and the boundary of the model box were in good contact and exhibited no relative sliding.e frictional boundary was well treated and effectively simulated boundary conditions of real slope.
Figure 9 reveals that the time-history curves were in good agreement; therefore, the flexible boundary of the model box wall was well treated in the test.e difference between peak acceleration values of two measuring points near the box wall at the same elevation point was not substantial, confirming that the flexible boundary essentially met the test requirements to reduce the influence of boundary effect of model box on test results.

Modal Parameter Analysis.
Modal parameter analysis is an important part of dynamic signal processing in shaking table tests.e purpose is to determine the modal damping ratio, natural frequency, and other modal parameters from obtained signals.According to the loading scheme, white noise with 0.05 g amplitude was used to excite the slope.e frequency response function was obtained by the acceleration response time history of the measuring points.en the dynamic characteristics of the model, such as natural frequency and damping ratio, were obtained by modal parameter identification technology.
Taking the convex slope as an example, the modal parameters of measuring points excited by white noise under various working conditions are listed in Table 6 with the variation trend illustrated in Figure 10.
e natural frequencies and damping ratios of different points were identical under the same conditions, indicating that slope modal parameters were not related to elevation.As the white noise gradually increased, the first-order natural frequencies of the model exhibited a gradual decline while the damping ratios gradually increased.is phenomenon occurred because the dynamic shear strength and dynamic shear modulus of rock and soil decreased, whereas the nonlinear characteristics strengthened under vibration excitation.

Dynamic Response of Model Slope
e dynamic response of a slope under seismic load can be used to judge slope stability, representing an important means of predicting high-slope deformation and failure.
e dynamic response of the model slope induced acceleration, velocity, and displacement.Examinations of seismic damage have revealed that seismic inertia force related to acceleration is the main cause of slope deformation and instability; many scholars thus take acceleration as the object of dynamic analysis [40][41][42][43][44]. is  e PGA amplification coefficient P AX was selected as the index to describe the acceleration response law, defined as the ratio of the peak acceleration value of each measuring point A X to the peak acceleration value measured on the shaking table surface A 0 .e acceleration amplitudes of El Centro waves were 0.1 g, 0.2 g, 0.3 g, 0.4 g, 0.6 g, 0.8 g, and 1.0 g, respectively.

Dynamic Response of Acceleration. According to
e acceleration amplification effect near the slope line and slope shoulder of the two models was most obvious.Taking ELC-3 as an example, the peak acceleration at the shoulder and slope line of the two slopes increased by 102% and 59% and 114% and 109% compared 8 Shock and Vibration with those measured on the shaking table, and the peak acceleration at other measuring points increased to varying degrees.

Dynamic Response of Velocity and Displacement.
Like the PGA amplification coefficient, the PGV amplification coefficient and PDG amplification coefficient are defined as the

Shock and Vibration
ratio of a certain measured point value in the model to the measured point value on the platform.Taking ELC-3 as an example, the PGV and PGD amplification coefficients of the slope shoulder, slope line, and slope toe were obtained by analyzing the velocity and displacement time-history curves.According to Figure 11, the PGV amplification coefficient measured on the convex slope shoulder was 1.72 whereas those measured on the slope line and slope toe were between 1.2 and 1.6.e PGV amplification coefficient measured on the concave slope shoulder was 1.81, and those measured on the slope line and slope toe were approximately 1.5.e velocity dynamic response on the two slopes was remarkable.e PGD amplification coefficient of the convex and concave slopes from toe to shoulder showed an increasing trend.Coefficients of the convex and concave slope were 1.25 and 1.29, respectively.ese findings suggest that the PGV and PGD amplification effects of the two slopes were each significant.
e velocity and displacement magnification effect of the concave slope was clearly stronger than that of the convex slope.
Due to the error of the test equipment itself, the accuracy of displacement values obtained by integration needs to be verified by displacement sensor S1 and S2.Taking ELC-1 as an example, by analyzing the displacement time-history curve obtained by displacement sensor, the displacement peak value of sensor S1 on convex slope is 1.35 mm and that of sensor S2 on concave slope is 1.40 mm.
e displacement peak values obtained by integration at the shoulder of concave and convex slopes are 1.36 mm and 1.41 mm, respectively (since there is no acceleration sensors at the top of the slope, the data of acceleration sensor at the shoulder nearest to the slope top were used).By comparison, it can be found that the displacement peak values obtained by displacement sensors are little different from those obtained by integration.erefore, the instantaneous displacement peak values integrated by acceleration timehistory curve can be used as the maximum displacement value of slope [45][46][47].

Influencing Factors of Dynamic Response of High Slope.
e dynamic response of a high slope under seismic load is closely related to the geometric parameters of the slope, mechanical parameters of the soil, and seismic load parameters.Studying the influences of various parameters on the dynamic response of the slope can further reveal the failure mechanism of slope instability.In this experiment, a set of rock slope models with different slope shapes was designed, and four layers of weak interlayer structural planes were designed inside the slope to ensure the model structure approximated the high slope in southwestern Yunnan.Acceleration data and the effects of slope shape, landforms, elevation, and seismic parameters on the dynamic response of acceleration were examined.
Taking ELC-3 and SIN-1 as examples (Figure 12), similarities were identified between the two slope models: PGA amplification coefficients of the convex and concave slopes exhibited an overall upward trend along the elevation with a broken line.
e PGA amplification coefficients of the two slopes declined sharply near the slope line, especially at the measuring points (40 cm high) below the slope line.e PGA amplification coefficients of the slopes rose sharply at the slope shoulder, which was substantially higher than the other points.In terms of differences, the PGA amplification coefficients measured at each point on the concave slope were slightly larger than those on the convex slope, especially measured above the slope line.e amplification effect on the concave slope shoulder (80 cm high) was more obvious than that on the concave slope; therefore, the concave slope shoulder sustained damage before the convex slope.
e influences of topography and elevation on the dynamic response of slopes under seismic loads have been studied extensively.Regarding landforms, the slope surface is controlled by factors such as the free surface; the magnification effect differs from that of the rock mass in the slope in the horizontal direction.e effect of elevation is due to the elevation amplification effect.Taking ELC-1 and SIN-1 as examples, the curve of PGA amplification coefficients varied with elevation.Figure 13 shows that the PGA amplification coefficients measured on both slope surfaces were significantly larger than those inside the slopes, and the difference was larger near the slope line at the same elevation.e gap narrowed gradually with an increase in elevation, and the two curves essentially retained a parallel upward trend.Under the ELC-1 condition, PGA amplification coefficients measured inside the slope below the slope line of the convex slope largely remained within 1.5.
e amplification effect was more obvious near the slope shoulder, with coefficients close to 1.8 exhibiting a sharp upward trend.
e PGA amplification coefficients measured on the slope surface exceeded 1.8 near the slope line, which demonstrated an initial downward trend followed by a dramatic increase near the line.All PGA amplification coefficients measured on the concave slope surface were greater than those inside the concave slope.Acceleration amplification peaked near the concave slope shoulder, with PGA amplification coefficients reaching approximately 2.0.
Figure 14 depicts the curve of PGA amplification coefficients along elevation under the SIN-1.Results show that amplification effect of convex slope surface was remarkable in the horizontal direction, and PGA amplification coefficients at each elevation point were greater than those inside the slope.Two curves grew in broken lines, exhibiting nonlinear characteristics, and the PGA amplification coefficients changed markedly at the slope line and slope shoulder.
Measuring points on the slope surface were not in the same vertical direction, hence three measuring points inside   Shock and Vibration the slope being selected as objects for analysis.Taking ELC-1 and SIN-1 as examples (Figure 15), the PGA amplification coefficients increased in line with elevation; they declined near the slope line under certain conditions.e influences of slope type, a weak interlayer, and ground motion parameters caused the magnification effect along the elevation to show an upward trend and minor reduction at certain positions.
ese findings indicate that the magnification effect of the concave slope was more significant.
Slope dynamic response is closely related to seismic characteristics of the slope, including the seismic wave type, seismic amplitude, frequency, and other characteristics.e experiment included two kinds of seismic wave input, an El Centro wave and sine wave, and took ELC-1 and SIN-1 as examples (Figure 16).
e PGA amplification coefficients corresponding to A4, A6, and A12 under SIN-1 were obviously larger than those under ELC-1.e acceleration amplification effect therefore differed under varying seismic waves, potentially due to the distinct spectral characteristics of different seismic waves.
Two loading schemes were designed for seismic waves.Under an El Centro wave, the acceleration amplitudes were divided into 6 grades: 0.1 g, 0.2 g, 0.3 g, 0.4 g, 0.6 g, and 0.8 g.Under a sine wave, the acceleration amplitudes were 0.1 g and 0.3 g.As indicated in Figures 17 and 18, PGA amplification coefficients showed an increasing trend at small amplitudes (0.1 g-0.3 g), remained stable between 0.4 g and 0.6 g, and then declined once the amplitude exceeded 0.6 g.
e PGA amplification coefficients of A3 exhibited an initial increase from 0.1 g to 0.3 g followed by a decline, and the range was wider.A9 and A14 increased initially followed by a stable stage and finally fell.
e PGA amplification 12 Shock and Vibration coefficients of the convex slope surface showed a small fluctuation in the early stage and continued increasing before gradually decreasing.Amplification coefficients measured inside the concave slope increased gradually from 0.1 g to 0.3 g, and A13 increased the most.In the later period, PGA amplification coefficients decreased gradually, and A1 declined most remarkably.All measuring points on the slope showed the above characteristics; that is, with an increase in ground motion amplitudes, PGA amplification coefficients exhibited an initial increase followed by a stable stage and ultimate decline.e greater the acceleration amplitudes, the higher the dynamic stress level; therefore, the shear modulus declined while the damping ratio increased, and the vibration-insulating effect on earthquake waves of the slope strengthened.
Figure 19 presents the variation law of PGA amplification coefficients along elevation under SIN-1 and SIN-5.Findings show that PGA amplification coefficients measured inside the slope and on the slope surface generally increased,   Shock and Vibration whereas coefficients measured on the slope surface declined near the slope line.PGA amplification coefficients under 0.3 g acceleration amplitude were larger than those under 0.1 g, indicating that the vibration intensity had a significant effect on acceleration amplification.Yet the curves of the two conditions were mostly consistent, implying that vibration intensity exerted little effect on the dynamic response distribution of the slope along the elevation.ese results reveal that the distribution law of PGA amplification coefficients of the slope was similar under different seismic loads.Vibration intensity had little effect on the dynamic response law of the slope along elevation, but it had a clear influence on the acceleration magnification effect.
According to the principle of seismic dynamics, when the load frequency is the same as the natural frequency of the slope, the slope will exhibit resonance.Natural frequencies of our model were determined according to Fourier spectrum analysis.A self-vibration test of the model box was conducted before the test, providing a basis for the reasonable selection of test frequency.e loading scheme was designed for the sine wave frequency at 5, 7.5, 10, and 12.5 Hz to analyze the effect of frequency on the dynamic response of the slope.
Figure 20(a) presents the curve of the PGA amplification coefficient of the convex slope at various frequencies under different sine waves at 0.1 g acceleration amplitude.With an increase in frequency, the PGA amplification coefficients measured on A9 and A14 showed an upward trend, whereas that of A3 remained stable first and then declined sharply.
e PGA amplification coefficients of A4 and A12 on the slope surface retained an initial slight upward tendency followed by a sharp decline.
e PGA amplification coefficients of A6 and A10 exhibited an upward trend along the elevation; the growth was rapid and then slowed.Figure 20(b) shows the curve of PGA amplification coefficients of the concave slope with different frequencies under sine waves with 0.1 g acceleration amplitude.e PGA amplification coefficients measured in the concave slope rose with an increase in frequency.e coefficients measured on A2 and A11 exhibited a slow growth trend, and those on A5 and A8 demonstrated a generally upward trend with a slight decline around 10 Hz.

Phenomenon of Dynamic Damage.
is section describes the dynamic failure process of the model slope under seismic loading.e concave slope showed no obvious deformation phenomena under the first few working conditions.Figure 21(a) shows a slight slip of the upper slope under condition SIN-5.
Figure 21(b) reveals that from working conditions SIN-6 to SIN-7, cracks in the model slope developed further.Cracks appeared on both sides of the slope shoulder under SIN-6.Under ELC-6 and ELC-7 conditions, a 5 cm long crack appeared on the top of the slope, extending from the back of the model to the slope surface at a 45-degree angle, and slight dislocation occurred at the interface of the muddy interlayer.
From conditions SIN-8 to SIN-9 (Figure 21(c)), the slope crevices ran through the lower part of the slope line to form a slip surface, and the upper slope sheared along the muddy interface to the free surface with bulges in the slope toe before the whole slope was destroyed.
For the convex slope, under the first few working conditions, no obvious damage was observed.When loading SIN-6, cracks appeared in the upper part of the shoulder, extending from the edge of the slope to the inner part of the slope; cracks also appeared at the interface of the weak interlayer.Under condition ELC-7, cracks in the layered  e dynamic failure modes of the concave slope included collapse of the whole slope near the slope line under strong earthquake conditions and shear sliding along the weak interlayer.At the same time, longitudinal cracks in the slope connected with horizontal cracks, leading to model destruction.e concave slope demonstrated an obvious collapse in the upper part of the slope line, and large cracks appeared in the slope body; then, the slope toe bulged.e convex slope collapsed near the slope line, and many interlaced cracks appeared inside the slope.

Influence Factors of Slope
e overall stability of the convex slope was thus stronger than that of the concave slope under seismic load, and the concave slope was more sensitive to the acceleration dynamic response than the convex slope.Moreover, slippage occurred in varying degrees at the interface of the four mud layers, with clearer slippage at the upper part of the slope than at the lower part.On one hand, the effect of elevation amplification was remarkable; on the other, muddy interlayers comprised the weak area of the slope.ese areas were subject to seismic loads, which could easily induce overall slope failure.
e influences of seismic load conditions on slope failure were primarily reflected in amplitude and duration.Under the first few working conditions with relatively small amplitudes, the amplitude and duration exerted little effect on the model with nearly no deformation.With an increase in amplitude and duration, the dynamic response of the slope model increased, and cracks and slips appeared.e cracks extended and expanded until the model was destroyed.
e stronger the vibration intensity of the seismic waves on the slope, the greater the seismic inertia force and the more damage to the slope.e total energy of ground motion input into the slope increased with an increase in duration, leading to accumulating slope damage.Under continuous loading, slope deformation became increasingly large, and slope displacement rose continuously, eventually resulting in slope failure.When the ground motion parameters were identical, the destructive force of the sine wave on the slope model was notably stronger than that of the El Centro wave [48,49].

Conclusions
rough a shaking table test, the failure modes and dynamic response law of a high slope were studied.Based on test data, the effects of landform, slope type, elevation, and ground motion parameters on dynamic response of the slope were analyzed in depth.e following conclusions were drawn: (1) e differences in PGA amplification coefficients between symmetrically measured points on the Shock and Vibration convex and concave slopes were not overly large.PGA amplification coefficients of the concave slope were slightly larger than those of the convex slope.Slope type influenced the acceleration response; the convex slope was more stable than the concave slope.Coefficients measured on the slope surface were always greater than those inside the slope.e PGA amplification coefficients measured inside the slope and on the slope surface increased by line, exhibiting nonlinear characteristics.
(2) Seismic wave type had a significant influence on dynamic response of the slope.e spectrum characteristics of different seismic waves varied greatly, leading to different acceleration amplification effects.Under the same conditions, the acceleration amplification effect of the sine wave was more significant than that of the El Centro wave.With an increase in seismic amplitude, PGA amplification coefficients increased initially followed by a stable stage and later declined.With an increase in frequency, the dynamic response of acceleration presented a nonlinear change.e dynamic amplification effect of the slope was closely related to the natural frequency.Within a certain range, the higher the frequency, the more significant the amplification effect.
(3) Different slope types introduced diverse failure modes, destruction positions, and degrees of failure.e concave slope was dominated by shoulder collapse, whereas the convex slope primarily collapsed near the slope line.Structural characteristics of the slope controlled slope deformation and failure.e position of the muddy interlayer was weak and prone to collapsing and slipping, which reduced slope stability.(4) Higher vibration intensity induced greater slope deformation and failure.e duration was closely related to the total energy input by ground motion.e duration increased and ground motion energy accumulated continuously, resulting in continuous accumulation of slope damage and ultimate instability.

2
Shock and Vibration aggregates, the model was built by stratified compaction according to the ratio of 32 : 56 : 9 by controlling water content.Four layers of a muddy interlayer interface (made of clay) were included with a thickness of 2 cm in the model slope, placed 0.2 m, 0.5 m, 0.7 m, and 0.9 m away from the bottom of the box, respectively.e height of the bottom step was 15 cm, which was paved and compacted on the bottom of the box to increase model stability.Concave and convex slopes were placed in the same model box, separated by partitions and subjected to the same input motion.

Figure 1 :
Figure 1: Location map of the study area.

Figure 4 :
Figure 4: Photos of completed models and the key steps of model construction.(a) View 1 of finished models.(b) View 2 of finished models.(c) Laying of muddy interlayer interface.(d) Compaction of model.(e) Cutting ring method.

Remarks: ( 1 )
Microseismic tests were carried out for the models by inputting white noise at the beginning of the test to measure the dynamic characteristics of the model.(2) White noise was input into the models after loading the seismic wave to measure changes in the dynamic characteristics of the model.(3) Cross-input of two seismic waves was used to observe and compare the effects of different types of seismic waves on the model slope.Shock and Vibrationexperiment mainly studies the response laws of acceleration, velocity, and displacement of the model slope.In the test, we monitored acceleration responses of the slope from different measuring points by inputting seismic waves with different ground motion parameters.e mass and volume of displacement sensor are too large to be installed inside the slope.It could affect the failure modes of the slope and the accuracy of the test.erefore, the sensors embedded in the model are all acceleration sensors.Velocity and displacement time-history curves were obtained by integrating acceleration time-history curves.

Figure 9 :
Figure 9: (a) Comparison of acceleration time-history curves between A15 and A1.(b) Comparison of acceleration time-history curves between A16 and A3.

Figure 10 :
Figure 10: Variation rules of modal parameters of model.

Figure 18 :
Figure 18: PGA amplification coefficients of concave slope under different El Centro wave.

Figure 19 :
Figure 19: (a) PGA amplification coefficients of convex slope under different sine waves.(b) PGA amplification coefficients of concave slope under different sine waves.

Figure 20 :Figure 21 :
Figure 20: (a) PGA amplification coefficients of convex slope under sine waves with different frequencies.(b) PGA amplification coefficients of concave slope under sine waves with different frequencies.

Figure 21 :
Figure 21: Illustration of the breakage process of concave slope under seismic loading.(a) First stage of model failure.(b) Second stage of model failure.(c) ird stage of model failure.

Figure 22 :
Figure 22: Illustration of the breakage process of convex slope under seismic loading.

Table 1 :
Figure 3: Shaking table system.(a) Shaking table surface.(b) Dynamic system of shaking table.Technical parameters of shaking table.

Table 2 :
e similarity ratios of the test model.

Table 3 :
Types and properties of model materials.

Table 4 :
Main physicomechanical parameters of prototype and model materials.

Table 5 :
e test loading scheme.

Table 7 ,
the peak acceleration of each measuring point on the slope enlarged to varying degrees relative to the measured point on the shaking table.

Table 6 :
First-order modal parameters for convex slope internal measurement points.

Table 7 :
Peak acceleration of each measuring point under El Centro waves (m/s 2 ).