Analysis of Rock β-Dynamic Parameters and the Stability of Earthquake Dangerous Rocks Based on PFC

Mesoparameters of rock materials are the main factors affecting the macromechanical properties of dangerous rock slopes. Based on the principle of particle flow and synthetic rockmass technology (SRM), the influence of mesoparameters onmacromechanical properties is investigated by calibrating mesoparameters of rock materials at depth for a rock sequence in Beichuan Qiang Autonomous County, Sichuan Province, China. By combining these parameters with conventional and dynamic cycle triaxial tests, sensitivity analysis of rock β-parameters was completed. As a result, the reliability of mesoparameters in the simulation of dangerous rocks is strengthened, providing a basis to examine the failure mechanism of earthquake dangerous rocks in this region. Results indicate that, in the triaxial test, sandstone failed in tension, and brittleness gradually weakened as confining pressure increased. Mudstone recorded shear failure, and the characteristic value of brittle attenuation showed a V-shaped change with increasing confining pressure. Under cyclic loading, cracks had a degrading effect on the damping ration (β) and the damping coefficient (C) of sandstone. Mudstone recorded relatively low β and low brittleness whilst sandstone had high β and high brittleness. In rock materials, βn is more sensitive than βs in mechanical properties. When the value of the βn-parameter was between 0.2 and 0.3 and the value of the βs-parameter was between 0.2 and 0.6, rock brittleness was more stable, and the reflected macroscopic mechanical properties were the most authentic. By using a deepened mesoparameter trial adjustment method, the failure mode of the Particle Flow Code (PFC) dangerous rock model near provincial highway 205, simulated under conditions for theWenchuan earthquake, indicated a tensile fracture-horizontal slip failure.-e simulated failure mode was consistent with that of real dangerous rocks, with the failure trend being concentrated between the first and the third layer of the rock mass.


Introduction
Sichuan Province, China, located in the Helan-Liupan-Longmen-Hengduan Mountain seismic zone, experiences frequent earthquakes and rockslides. In this area, earthquake-induced rock slope instability often occurs, resulting in significant economic and property loss. Dangerous rock slopes have been identified as the main formation resulting in earthquake-type slope collapse disasters. Rock masses are generally cut by different structural planes. When the main control structural plane is located at the rear of the dangerous rock block, it gradually penetrates the rock mass until the mass breaks under the action of multiple factors. Under this condition, the dangerous rock block becomes unstable and is destroyed. Although previous investigations have predominantly examined instability and failure of dangerous rocks, studies examining the stability of dangerous rocks in earthquake areas are typically based on the principle of particle flow; reliable rock material mesoparameter calibration methods are urgently needed. is area of research has become one of the main research areas and technical challenges in this field [1,2]. From the perspective of rock materials, the discrete element method has become important to solve critical rock simulation technology, and the key condition to determine the reliability of the method is the selection of mesoparameters in the model.
For seismic slopes, Chen [3] deduced the chain development process of dangerous rocks in the ree Gorges Reservoir area. From the perspective of geomorphology, they revealed the influence of rock cavern weathering on the stability of dangerous rock slopes, proposing the dangerous rock chain pattern development. Based on fracture mechanics and damage mechanics, Tang et al. [4,5] discussed the sequence of dangerous rock caving on a cliff with a weak base, establishing a calculation method for the time of dangerous rock caving. Yan et al. [6] used the QUIVER code to realize the response analysis of earthquake gentle slopes, identifying slope response and sliding displacement law under the action of an earthquake. e movement characteristics of dangerous rock slopes under different seismic conditions were simulated by Huang et al. [7] using the discontinuous deformation analysis (DDA) method. e influence of seismic loads on the movement characteristics of dangerous rock mass collapse blocks was shown using the DDA method. Zhao et al. [8] used particle flow PFC2D software to dynamically simulate the instability process of rock slopes, revealing the influence of the distribution of weak structural planes on the instability and failure of rock slopes. Tang et al. [9] used a two-dimensional discrete element simulation to study the formation mechanism of the Caoling landslide in 1941.
For studies of rock materials, Cong et al. [10] proposed that the macroscopic mechanical characteristics of materials are closely related to microscopic parameters. ese parameters can be combined with marble indoor loading and unloading tests to determine the mesoparameters suitable for rock materials (such as marble). He et al. [11] comparatively analyzed the damping characteristics of sandstone, conglomerate, and glutenite under different stress paths, highlighting the evolution of the damping ratio (β) and the damping coefficient (C) with the number of cycles. e influence of fractures on rock mechanical properties proposed by Huang et al. [12] has important practical significance for the engineering stability of fractured rock masses. In addition, strength, deformation parameters, and failure modes of rock materials with discontinuous, nonparallel double fissures were analyzed in detail. e mesoparameters of PFC2D numerical sandstone samples from the mesolevel aspects of displacement, crack evolution, and distribution were studied by Zhou et al. [13]. Tang and Zhang [14][15][16] studied the shear behavior of rock fractures from the perspective of temperature and provided new ideas for the stability evaluation of surrounding rock in temperature-dependent underground engineering. Tang and Jiao [17] also proposed that, in the contact behavior of rock joints, the three-point peak (3PP) standard seems to be the most appropriate method to identify the contour peaks of rock joints. Zou et al. [18] used the extended finite element method (XFEM) to simulate the propagation of hydraulic fractures in heterogeneous reservoirs. ey proposed that the hard blocks in the reservoir promote the growth of hydraulic fractures, but limit the growth of fracture widths.
As the majority of previous investigations considered the two-dimensional particle flow method, dangerous rock slope simulation technology lacks reliability due to insufficient factors being considered, resulting in low applicability of earthquake dangerous rock simulation technology in earthquake areas. Furthermore, the mechanical characteristics and mesoscopic parameters of rock materials need further evaluation. In particle flow rock material simulation, the parallel bond contact model is generally used to assign values to the contact between rock particles. Since Poisson's ratio μ, elastic modulus E * , stiffness ratio k * , and other mesoparameters will directly affect the macromechanical properties of rock materials under numerical tests, they have become the main objects considered in the trial adjustment processes [19][20][21][22][23][24][25]. However, under the dynamic action of a rock slope, the damping ratio (β) is an important parameter that affects the dynamic response of a rock [11]. is influence on the macromechanical properties of rock materials has generally been overlooked in previous investigations. erefore, our study includes crack tracking, strength attenuation index (brittleness characteristics), and β-dynamic parameter analysis on the basis of traditional parameter adjustment methods. By selecting different mesoparameters from previous studies, we investigated normal critical damping ratio β n and shear critical damping ratio β s for sensitivity analysis. As a result, rock mesoparameters that can better reflect true materials are obtained, providing a strong basis for simulating the seismic damage of dangerous rocks and analyzing the stability of earthquake dangerous rocks.

Sample Equipment and Sample Preparation
2.1.1. Sample Equipment. Mudstone and sandstone mechanical tests in this study were performed on a RMT-150C rock mechanics test system ( Figure 1) which can realize automatic data collection and processing. e maximum axial load of the testing machine is 1000 kN, the vibration frequency is 0.2-2 Hz, the specimen size is Φ (50 × 100)-(100 × 300) mm, the maximum three-axis confining pressure is 100 MPa, the confining pressure rate is 0.001-1 MPa/s, the displacement measurement range is ±50 mm, and the strain measurement range is ±0.01. e source vibration wave can be sine wave, triangle wave, square wave, oblique wave, etc.

Sample Preparation.
One set of mudstone and three sets of sandstone ( Figure 2) were used in this investigation. e static mechanics test was divided into one set of mudstone specimens (No. A) and one set of sandstone specimens (No. B). e dynamic test was divided into one group of conventional sandstone specimens (No. C) and one group of fractured sandstone specimens (No. D) (Figure 3). Rock samples were collected from dangerous rock on Provincial Highway 205 in Beichuan Qiang Autonomous County. e test sample was a cylinder with a diameter of 50 mm and a height of 100 mm. e error range was controlled within ±2 mm, meeting the requirement of a 2 : 1 height-to-diameter ratio specified by the International Society of Rock Mechanics. e cutting position of the cracked sandstone was set as the top central axis; the cutting depth was 50 mm, and the gap width was 2 mm.  [26,27]. e failure mechanism was that the material experiences shear failure under the action of normal stress. e maximum shear stress τ m it bears was determined by the cohesion and the angle of internal friction, expressed as where c is the cohesive force, σ is the normal stress on the shear failure surface, and φ is the internal friction angle. In order to further reflect the brittleness of the rock, we examined the mechanical index that characterizes the strength attenuation behavior of rock samples after the peak [28]. e specific algorithm was as follows: where D s is the intensity attenuation coefficient and the value range is [0, 1]; Δσ is the intensity attenuation value; σ P 1 is the peak intensity; and σ r is the residual intensity.

Dynamic Cycle Triaxial
Test. e dynamic cycle triaxial test was used to study the dynamic mechanical characteristics of the rock samples under three-dimensional stress. e test simulates not only the static stress state of the rock on-site, but also the actual seismic effect.
is test provides a reference for studying the dynamic parameters of sandstone and the dynamic stability of dangerous rocks. In the cyclic loading process of an ideal elastic material, the stress and strain time history of the elastic body is synchronized. Rock is a typical inhomogeneous material, with  e cement interface between mineral particles is also a relatively weak surface. e area of the hysteresis loop can reflect the amount of energy consumed in the cyclic loading process. Under ideal elasticity, the vibration damping ratio β can be calculated according to the area of the hysteresis loop; the average slope of the hysteresis loop reflects the magnitude of the dynamic elastic modulus E d (Figure 4) [29]. e values of damping ratio β, damping coefficient C, and E d can be calculated as follows [30]: where A is the area of the hysteresis loop ABCDA; A s is the area of triangle AOE; X is the amplitude of vibration; ω is the angular frequency of vibration; σ d max is the maximum axial dynamic stress of the hysteresis loop; and ε d max is the maximum axial dynamic strain of the hysteresis loop. As dangerous rocks in this area are mainly composed of sandstone, this cyclic test was only conducted on sandstone samples. Confining pressures were 2 MPa, 4 MPa, and 6 MPa, and the vibration frequency was 1 Hz; the loading rate was 60 kN/min. e loading process was divided into six stages of cyclic loading: (1) when the hoop load reached the predetermined confining pressure, the first stage was loaded from 0 to 20 kN and then loaded cyclically between 0 and 20 kN at a frequency of 1 Hz 30 times and then unloaded to 0. (2) e second stage was loaded from 0 to 40 kN, loaded cyclically from 0 to 40 kN at a frequency of 1 Hz 30 times, and finally unloaded to 0. Stages 3 to 6 continued in this format until 120 kN was used in stage 6. During the step-bystep loading process, if the test piece was damaged, the test ended.

Numerical ree-Axis Servo Test.
In this study, the PFC3D simulated triaxial test was used to conduct microscopic failure analysis and mesoparameter calibration on two kinds of rocks. e entire three-axis numerical simulation loads the sample by specifying the top and bottom wall boundary constraint speeds, determining the stress and strain state of the sample through the built-in FISH function of PFC3D. e established loading model is shown in Figure 5. According to the real triaxial test, the model generates particles in a cylindrical wall with a height of 10 m and a diameter of 5 m. In preliminary model calculations, it was found that if the particles are generated according to the actual gradation, the number of particles will reach 1 million. However, the existing calculation conditions cannot complete the numerical analysis of such a large number of particles. In this case, after considering the impact of the particle size effect on the test results [31], the radius of the particles is appropriately expanded, and the radius is 0.2-0.25 m. e purpose of setting the sample size in this way is to keep the particle size of the numerical triaxial sample and the numerical slope model at the same size (that is, the unit is m). At the same time, it not only improves the calculation efficiency, but also does not affect the accuracy of the simulation results. Under the condition of constant confining pressure, the upper and lower loading plates moved towards each other at a speed of 10 mm/min until the specimen was broken. e stress-strain curve was automatically recorded throughout the process. In addition, the confining pressure was loaded according to the amount of 2 MPa, 4 MPa, 6 MPa, and 8 MPa in the indoor test. Finally, the overall project view and flowchart of the analysis steps is shown in Figure 6.   Advances in Civil Engineering

Conventional Triaxial Test
Analysis. e macromechanical properties of sandstone and mudstone under the two triaxial test methods are shown in Table 1. Comparison of results under the two methods indicates that cohesion c and the friction angle φ were essentially consistent. Among these, the cohesive force error was controlled within 0-0.02 MPa, and the friction angle error was controlled within 0.3°-0.4°, thereby preliminarily determining the rationality of the microscopic parameters.
In the process of the numerical triaxial test, the damaged contact links were expressed in the form of disc joints (DNF_disk), thereby forming damage cracks. e change law of strength attenuation value Δσ, strength attenuation coefficient D s , and failure characteristics of the two rock samples is shown in Figure 7. With an increase of confining pressure, Δσ for sandstone recorded a gradual trend of decrease, consistent with the findings of Zuo et al. [26]. e D s curve for sandstone recorded an obvious sharp decline (0.67 to 0.37) after a confining pressure of 2 MPa, after which the decline continued more slowly. is result indicates that the brittleness of sandstone gradually weakens with an increase of confining pressure, and the weakening ability of brittleness increases after a confining pressure of 2 MPa is reached. Δσ and D s mudstone curves recorded a V-shaped change with increasing pressure (Figure 7(b)). When the curve dropped to a 4 MPa confining pressure, Δσ and D s increased at the same time, by 7.05 MPa and 0.15, respectively; after 6 MPa, the rate of increase declined. is result indicates that, under a low confining pressure, due to large initial porosity, the brittleness of mudstone initially weakens and ductility is enhanced. Mudstone brittleness then increases, and ductility decreases at a high confining pressure.
Finally, the difference in Δσ and D s between the two rocks is due to the difference in initial porosity and initial pore water content between sandstone and mudstone. It can be said that, under the same external conditions, the strength of mudstone is generally smaller than that of sandstone, and the initial porosity and initial pore water ratio are also greater than those of mudstone. erefore, when the confining pressure of mudstone continues to increase, the pores will be compacted in the middle, the pore water will be discharged, and a new arrangement and combination of rock particles will be formed inside the rock.
is makes the strength decay curve different from sandstone.
Analysis of the partial failure characteristics of the two rock samples shows that the failure mode of sandstone is a tensile failure; mudstone has a shear failure.
ere were numerous tensile cracks in the microscopic joints of the sandstone, with tensile fractures causing the main damage among individual grains. When the confining pressure was 2 MPa, mesocracks were generated from the bottom of the sandstone which gradually developed towards the top. Due to the low confining pressure, macroscopic shear cracks did not penetrate the top, and the top mesocracks developed slowly. When the confining pressure was 8 MPa, the bottom crack gradually expanded to the top, forming a penetrating macroscopic shear crack. e sandstone samples failed under this confining pressure; the number of cracks gradually increased as the confining pressure increased, finally forming rock failure characterized by a multiple crack staggered development. Mudstone specimens mostly recorded shear cracks in the microscopic joints, with the majority of fractures occurring between particles being shear failures. Under a confining pressure of 2 MPa, a single shear failure was recorded for the mudstone,  Advances in Civil Engineering having fewer shear failure cracks and being characterized by obvious brittleness. e mudstone sample also recorded a few tension cracks from the base upwards. When the confining pressure increased to 8 MPa, microscopic cracks at the base and top gradually increased, propagating towards the middle. At this time, shear failure of mudstone cracks increased and became more obvious.
e failure characteristics and strength attenuation law of the two kinds of rocks in the indoor test and the PFC test were basically the same, adhering to the law proposed by Chen [32]. erefore, rationality of the material mesoscopic parameters was further verified.

Triaxial Cyclic Loading Test and β-Parameter.
e sandstone cyclic triaxial test dynamic stress-strain curve and its failure characteristics ( Figure 8) indicate that, under the action of dynamic circulation, the peak stress intensity of the sandstone was obviously lower than that of the conventional triaxial. is result indicates that the compressive strength of the sandstone was weakened under dynamic load. Sandstone without cracks has instantaneous elastic deformation at the moment of loading, and its strain value is generally between 0 and 0.2%; the strain value of sandstone with cracks is between 0 and 0.1%. Under a certain confining pressure, sandstone gradually exhibits obvious creep properties as load applied by stepwise loading increases. Eventually, under the action of the last level of load, the sandstone appeared to be destroyed by an accelerated creep stage. Our results indicate that as the confining pressure increased, the number of loading levels also increased. As the confining pressure increased from 0 to 6 MPa, the axial load increased when creep failure occurred. As the increase of confining pressure  improved the compressive strength of the sandstone, the occurrence of cracks resulted in the axial strain of the sandstone to increase. Under a confining pressure of 2 MPa, the number of cycles of sandstone with fractures was 36 times lower than that of sandstone without fractures, and the maximum peak stress declined by 10 MPa (25%). At a confining pressure of 4 MPa, the number of cycles of sandstone with fractures was reduced by four times compared with sandstone without fractures, and the maximum peak stress was reduced by 10 MPa (20%). At a confining pressure of 6 MPa, the number of cycles of sandstone with fractures was 31 times lower than that of sandstone without fractures, and the maximum peak stress declined by 10 MPa (16.7%). is result indicates that while the confining pressure increased, the yield strength of sandstone with fractures was one lower than that of sandstone without fractures, and it was more prone to failure under dynamic loads.
Analysis of failure characteristics indicates that sample C-02 has low confining pressure, low compressive strength, and low creep, and the sample is accompanied by partial tensile failure and primarily shear failure. As the confining pressure increased, the number of cycles increased stepwise and creep increased. Samples C-04 and C-08 recorded obvious simultaneous effects of tension and shear. In addition, tensile crack penetrated the entire sample, intersecting with the shear crack to form a Y-shaped crack in C-08. Results for cracked sandstone (Figure 8(b)) indicate that, when the low confining pressure is cyclically loaded, the crack expands outwards in the form of a tensile crack along with the reserved fracture, finally intersecting with the shear crack and failing.
In the cyclic loading stage, the damping ratio β and the damping coefficient C of the rock sample changed, and the condition of the cracks will affect β and C of the sandstone. By using equations (4) and (5), β and C were calculated. In order to facilitate the analysis of the influence of cracks on damping parameters, two sets of hysteresis loops with different cracks under a high confining pressure of 6 MPa were selected to calculate β and C. Finally, the relationship between the damping parameters of the two sets of sandstone samples and the number of cycles was calculated ( Figure 9). Results in Figure 9 indicate that, under the same confining pressure, the damping coefficient C increases linearly with an increase in the number of cycles, and it is not affected by the increase of the load amplitude. In addition, the damping coefficient C of sandstone with fractures is five times lower than that of sandstone without fractures. e damping ratio β noticeably fluctuates under the action of cyclic loading with a confining pressure of 6 MPa, gradually decreasing as the number of cycles increases. e β of fractured sandstone is generally 0.1% lower than that of sandstone without fractures. Under a higher confining pressure, the damping coefficient C varies with cycle times in the opposite direction to the damping ratio β. is law is consistent with the results of Zhou [29], indicating that cracks have a degrading effect on the damping coefficient C of sandstone. is is because during the trial adjustment of the mesoparameters, those that affect rock brittleness mainly depend on the critical damping ratio β. e critical damping of mudstone is relatively low, and its brittleness is low in numerical tests; as the critical damping of sandstone is relatively high, it is therefore more brittle. Based on these conclusions, the following sensitivity analysis of β-mesoparameters of sandstone and mudstone was carried out.
Results indicate that β n has an important influence on the peak strength of deviatoric stress of mudstone and sandstone. When β n had a value of 0 to 1.0, the stress-strain curves of the two materials recorded an obvious "threesegment" distribution, and the curve trend was different. For mudstone, the increase of β n resulted in peak deviatoric stress and residual stress to increase; results for strain recorded no significant influence. When β n was 0, 0.1, 0.5, 0.6, 0.7, 0.8, and 1.0, the peak deviatoric stress of mudstone was the largest, being 10-12 MPa larger than the peak deviator stress under actual mudstone parameters. is result demonstrates that the value of β n has an amplifying effect on the peak deviatoric stress of mudstone. When the β n value was 0.2 and 0.3, the peak deviatoric stress and initial deviatoric stress of mudstone were moderate and consistent with the deviatoric stress under actual parameters, and the difference range was within 2 MPa. is indicates that the mesoscopic parameters of mudstone under β n are the closest to the actual effect. When the β n value was 0.4 and 0.9, the strength of mudstone material was the lowest, and the peak deviatoric stress and initial deviatoric stress were 10-12 MPa lower than the actual value. e β n value has a degrading effect on the strength of mudstone, and the initial deviatoric stress direction is abnormal, making the initial deviatoric stress value negative.
For sandstone, an increase in β n reduced overall peak deviatoric stress and residual stress. e peak strain fluctuated slightly under different β n values, and peak strain was the largest when β n was 0.3. Considering the strength range, when β n was 0.1, 0.4, and 0.9, the peak deviatoric stress of sandstone was the largest, being 5-6 MPa larger than the peak deviatoric stress under actual sandstone parameters, and the initial deviatoric stress was also larger than the actual one. is result indicates that the β n value has an amplification effect on the peak deviatoric stress of sandstone, and the amplification value was within 10-12 MPa. When β n was 0.2 and 0.3, the initial deviatoric stress and peak deviatoric stress of sandstone were closest to the actual value. Although the error was within 2 MPa, the strain was larger than the actual value. When adjusting these parameters, it is recommended that after these values are tested and adjusted, other mesoparameters need to be considered to continue to improve the strain calibration. When β n was 0.5, 0.6, 0.7, and 1.0, the peak deviator stress was too small. In the same way, the β n parameter under this condition had a degrading effect on sandstone and, with abnormal initial conditions, it is not recommended to be used in parameter sensitivity test adjustments.
For the process of judging whether the particle flow rock material is close to the real material, only using the peak stress and the values of c and φ as the judgment basis is not recommended. Rock ductility and brittleness should also be considered. Based on this, we also considered the sensitivity of β-parameters under the intensity attenuation value and intensity attenuation coefficient. e influence of β n parameter on material brittleness ( Figure 11) indicates that, within the range 0-0.3, the strength attenuation value of sandstone is 12-15 MPa larger than that of mudstone, and a slight upward "U"-shaped fluctuation occurs with the attenuation coefficient. Here, the β n -parameter has little effect on the strength attenuation value and attenuation coefficient of the two materials, that is, the brittleness and ductility of the material remain unchanged. When β n increased to 1.0, the intensity attenuation value and attenuation coefficient significantly fluctuated as the damping ratio increased. e attenuation value of sandstone increased to 42.24 MPa when β n had a value of 0.4, and the attenuation coefficient increased to 0.68. At this time, the mudstone attenuation value did not fluctuate.
is result indicates that the β n value greatly improved the brittleness of the two kinds of rocks, and sandstone was more sensitive. When β n reached 0.5, the attenuation value of sandstone significantly dropped to 22.19 MPa (the attenuation coefficient dropped to 0.52); the attenuation value of mudstone began to increase to 25.63 MPa, and the attenuation coefficient decreased. Although this result shows that β n at this moment had greatly improved the ductility of sandstone and decreased its brittleness, its sensitivity to mudstone was not high. By taking all of our results into consideration, brittleness of the two rock materials was relatively stable when the β n value was 0-0.3 and the parameter value ranged between 0 and 1.0. At the same time, the trial adjustment of the β n -parameter between 0.2 and 0.3 was closer to the actual value, and the macroscopic mechanical properties were more stable.
As seen in Figure 12, the peak deviatoric stress of mudstone and sandstone materials did not significantly change with an increase in β s , and it was consistent with the deviatoric stress under the actual macroscopic test. is result indicates that the shear damping ratio had no great influence on the macroscopic yield strength of the two materials, and there was no need to consider its influence on the peak deviatoric stress during the mesoscale trial adjustment. e influence of β s -parameter on material brittleness is shown in Figure 13. β s -parameter had a certain influence on the brittleness of the two rock materials. For mudstone, the strength attenuation value and attenuation coefficient showed a "wave-shaped" change as β s increased. When β s was 0.1, the maximum attenuation value and attenuation coefficient reached 19.33 MPa and 0.54, respectively, and the minimum attenuation value and attenuation coefficient were 15.1 MPa and 0.43 when β s was 0.7, respectively. In the process of parameter trial adjustment, the range can be set within 0.1∼0.7, and the selection was based on the strength attenuation coefficient obtained from the indoor test. For sandstone, as β s increased, the strength attenuation curve recorded an "arch" change. When the β s value was 0.2, the brittleness of sandstone increased from 0.53 to 0.65 before remaining stable; when β s was 0.6, brittleness decreased to 0.53. According to the attenuation coefficient obtained from the indoor test, it was observed that the β s -parameter was within the range of 0.2 to 0.6 as a reasonable and sensitive test adjustment range. At this time, brittleness of the sandstone material was closest to the real material.  in Table 2. With reference to previous studies [33,34], the structure of surface mesoparameters (Table 3) was confirmed.

Model Establishment.
In this study, we used a typical steep slope of dangerous rock on Provincial Highway 205, Beichuan Qiang Autonomous County, Sichuan Province, China (Figure 14), as a case study. e 2008 Wenchuan Earthquake in Sichuan caused a series of instabilities, resulting in the collapse of dangerous rocks posing a significant threat to the safe operation of the provincial highway and nearby residents. After a detailed geological survey, the steep slope of dangerous rock (about 23 m high) was restored. e upper rock mass is composed of feldspar sandstone with good integrity, about 16 m high, with a natural bulk density of 24.7 kg/m 3 and a fracture toughness of 26.0 MPa m 1/2 ; the lower rock mass is mudstone with poor integrity, being about 7 m high with a natural bulk density of 25.3 kg/m 3 . e fracture penetration rate is 25%, and four sandstone layers were identified in the rock mass, including two main control structural planes. According to the chain law of rock slope [35], the external macrochain was called the No. 1 chain, and the internal macrochain was the No. 2 chain. Each macrochain contains four microchains, and the dangerous rock masses are numbered from bottom to top (the first chain was 1-1, 1-2, 1-3, and 1-4, and the second chain was 2-1, 2-2, 2-3, and 2-4).
According to survey results of on-site dangerous rocks, reasonable generalizations were made for this investigation. PFC3D was used to establish a particle flow model for dangerous rock slopes in three dimensions ( Figure 15). Here, the formation of the rock joint surface adopted artificial synthetic rock mass technology [36], and the linear parallel bond model was selected as the internal contact model of the rock; the smooth-joint contact model was the contact model of the joint surface between the rocks. e seismic model under the coupling action of transverse waves and longitudinal waves has an amplifying effect on stability, and the probability of both acting at the same point at the same time is small. It was therefore decided to use the first 20 s acceleration curve (Figure 16) of the real EW transverse wave monitored in Wolong during the Wenchuan earthquake as the seismic input. When applying seismic action, in order to prevent seismic waves from reflecting at the boundary of the model, a viscous boundary condition is required [37]. In PFC3D, this was achieved by setting dampers in the normal and tangential directions of the two sides and the bottom boundary of the model, as Table 2: Mesoscopic parameters of sandstone and mudstone materials.  The fracture penetration rate shown in Figure 17. During the loading process, seismic acceleration was transformed into boundary acceleration, resulting in the model particles to vibrate.

Failure Mode Analysis.
e destruction process of the dangerous rock model under a 20 s earthquake is shown in Figure 18. Model results indicate that, between a seismic shear wave time of 0 and 6 s, the rock mass did not undergo significant deformation or damage. At the same time, crack development in the lower mudstone was not recorded and slope stability was generally good. Comparison with the particle velocity cloud map indicated that the maximum particle movement velocity was concentrated on the bottom rock mass 1-1 and rock mass 1-2 at 3 s. At 6 s, the maximum particle velocity was transferred to the second and third layers of rock, indicating that, before the rock mass was destroyed, the failure trend was concentrated between the first and third layers. When calculation time reached 9 s, the seismic acceleration curve gradually peaked (see Figure 8). Here, the internal structural surface contact model of Nos. 1-1, 2-2, and 3-2 dangerous rocks was locally fractured, and cracks in rock mass 1-2 gradually developed. Approximately 10 s later, the 1-2 rock mass fractured. In addition, the top rock layer began to move, accompanied by the fracture of the smooth joint model between the layers. e particle velocity at this moment was the largest relative to the velocity at other moments, and the failure trend was concentrated between the first and third rock mass layers. When the earthquake rupture time was 11.5 s, the internal vertical rock mass fell outward as a whole, colliding with the external vertical rock mass, being pulled apart. Because of this movement, fissures between the rock masses were further enlarged. In addition, due to gravity and seismic action of the upper rock mass, particles in the lower mudstone base began to disintegrate. At 20 s, due to the interaction of adjacent rock blocks at the bottom under the earthquake, the top rock block (No. 4-1) broke locally. As this rock block was located on the lower rock mass, it therefore fluctuated with the earthquake. At the same time, the rest of the rock masses, except for No. 4-2, had a tendency to bulge outwards as a whole. However, based on the support of the upper and lower rock masses, no rock mass collapse occurred. In summary, the dangerous rock of Provincial Highway 205 in the high-intensity area will undergo slight instability failure under the action of a traverse EW-wave 20 s strong earthquake, and the failure mode will be a tensile fracture-horizontal slip failure.
In order to further verify the rationality of the rock material mesoparameters and the seismic damage of the PFC dangerous rock model, we introduced two other dangerous rock failure maps in the area (Figure 19). It can be seen that after the Wenchuan earthquake, the failure mode of the real dangerous rock was basically consistent with the failure mode of the earthquake simulated in our study.

Displacement Analysis.
Monitoring the vertical and horizontal displacement of the rock block in this study provides further information on the failure mode and stability of the rock block. Model displacement in the vertical z-direction and horizontal y-direction ( Figure 20) indicate that there is no sign of rock collapse or falling during the earthquake action period; the displacement value was recorded to be very small. Negative vertical displacement (distributed between 0 and 0.4 m) and   Figure 18 indicate that the rock block was damaged and fractured during this period. Although the mudstone base cracked and stability gradually decreased, caving and rolling failure did not occur. Based on the damage trend, it can be judged that rock blocks 4-1 and 1-1 will be the most likely to fall. It should be noted that the vertical displacements of No. 1-2 and No. 3-2 rock masses are large, possibly due to adjacent rock masses colliding and squeezing each other, resulting in surface rock mass particles to break and fall off.

Analysis of Unbalanced Force of the Rock Block.
By analyzing the failure modes of various slope models under earthquake conditions, it is accepted that the dynamic interaction of rock blocks is an important factor causing rock   mass collapse. In PFC, the unbalanced force of the particle not only represents the size of the contact force chain around the particle but also includes the tension and compression of the interlayer load and adjacent rock blocks. erefore, the unbalanced force in each rock block can represent the magnitude and direction of the load on each rock block. e unbalanced force fluctuation curve of each rock block, obtained using in situ monitoring, is shown in Figure 21. In the vertical z-direction, the rock mass unbalanced force predominantly fluctuates between 8 and 13 seconds. is wave phenomenon is caused by rocks caving and colliding with each other. e vertical upward peak unbalanced force generated by the fracture was 230 kN, generated by rock block 4-1; the vertical downward peak unbalanced force was 350 kN, generated by rock blocks 1-2 and 3-2. e maximum horizontal y-direction unbalanced force was distributed in rock blocks 1-2, and the peaks were concentrated in the maximum seismic acceleration interval. It can be seen that

Conclusions
Based on the principle of discrete element particle flow, we conducted an in-depth analysis of the mechanical characteristics and mesoparameters (β-parameters) of rock materials through dynamic and static triaxial tests of sandstone and mudstone. Findings from this investigation enhance the reliability of rock material mesoparameter calibration technology. Furthermore, a more realistic and reliable dangerous rock model compared to previous studies was presented, with seismic stability analysis undertaken. Conclusions from this study are as follows: (1) e brittleness of sandstone gradually weakened with an increase in confining pressure, with brittleness weakening more severely after a confining pressure of 2 MPa. In addition, the characteristic value of mudstone brittleness attenuation recorded a V-shaped change with increasing confining pressure. Under low confining pressure, the brittleness of mudstone was preferentially weakened and ductility was enhanced. Under high confining pressure, the brittleness of mudstone was enhanced and ductility decreased. During the PFC numerical three-axis test, mudstone particles were recorded to undergo shear failure and the sandstone particles underwent tensile failure. e failure characteristics of the two types of rocks in the PFC numerical test are in accordance with results from an indoor test. (2) Fissures have a degrading effect on β and C of sandstone. When the rock fails, the dynamic compressive strength of sandstone with cracks is one lower than that of sandstone without cracks. e peak stress of the compressive strength of the former is reduced by 16.7%∼25% compared with the latter.
is finding indicates that sandstone, especially sandstone with many cracks, is more prone to damage under dynamic load. Under the same confining pressure, the damping coefficient increases linearly with the increase of the number of cycles, and it is not affected by the increase of the loading amplitude.
e damping coefficient of sandstone with fractures is generally lower than that of sandstone without fractures by a factor of five.
(3) During the trial adjustment of mesoparameters, the peak deviatoric stress and c − φ value are not the only criteria for judging the authenticity of materials. e effect of material brittleness is equally important to mesoparameters and cannot be ignored. In addition, the mesoparameters of rock brittleness mainly depend on the damping ratio β. Mudstone has a relatively low β and low brittleness whereas sandstone has a high β and high brittleness. In rock materials, β n is more sensitive than β s . e value of the β n -parameter was between 0.2 and 0.3, and the value of β s -parameter was between 0.2 and 0.6. Rock brittleness was more stable and the reflected macroscopic mechanical properties were more realistic.
With an increase of β s , Δ s and D s of mudstone indicated "wave-shaped" changes. In contrast, sandstone recorded "arch" changes. Advances in Civil Engineering 15 slip failure. e simulated failure mode was basically consistent with results recorded for the real dangerous rock. Combined with displacement analysis, the failure trend was concentrated between the first layer to the third layer of the rock mass. e failure characteristics of the model in the vertical and horizontal directions were small, and the horizontal influence was relatively large. In addition, the dangerous rock did not eventually collapse and fall in a large area. It is judged that rock blocks 4-1 and 1-1 are likely to fall in the later stage of an earthquake. Combined with analysis of the unbalanced force, collisions of rock blocks 4-1, 1-2, and 3-2 are the most violent under the action of the earthquake, especially rock block 1-2.
Data Availability e table and figure data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.