Dynamic Response Study of Impulsive Force of Debris Flow Evaluation and Flexible Retaining Structure Based on SPH-DEM-FEM Coupling

Flexible retaining structure is demonstrated to be an effective measure for debris flow prevention in mountainous areas, which can effectively separate water and stone, reduce particle mass, and dissipate kinetic energy. In order to explore the impulsive force and dynamic response of flexible retaining structure impacted by solid-liquid two-phase debris flow, a complex dynamic interaction model of particle-fluid-structure has been established by employing the SPH-DEM-FEM coupling numerical analysis method. In the present study, the process, impulse force, and dynamic response of a flexible retaining structure subjected to debris flow under different slopes were investigated, respectively. ,e numerical results are compared with those of the calculation formula of the peak impulsive force of the semiempirical debris flow. Meanwhile, the effectiveness of the coupling numerical simulation is verified. ,e simulation results show that the coupled SPH-DEM-FEM numerical analysis method can visualise the impact of the solid-liquid two-phase debris flow on the flexible retaining structure, reproducing the impact, retaining, water-stone separation, climbing, silting, and deposition process of debris flow. ,e dynamic time-history curve of the coupling numerical analysis method for the flexible retaining structure is consistent with the results of the existing literature.,e debris flow evaluation results of flexible retaining structure impacted by solid-liquid two-phase debris flow are in an order of magnitude with the empirical formula results. Moreover, the results are reliable. ,e obtained results have a certain reference value for the study on the impulsive force and dynamic response of the flexible retaining structure impacted by solid-liquid two-phase debris flow and the engineering design of the flexible retaining structure.


Introduction
Debris flow is a solid-liquid two-phase mixed fluid containing a large number of rubbles and sediments, showing the movement characteristics of viscous laminar flow or dilute turbulent flow. It is the result of the combined action of topography, meteorology, hydrology, soil, and vegetation in the basin [1]. Debris flows are characterized by rapid flow velocity, large flow velocity, sudden eruption, and amazing destructive power. It is a special fluid-solid coupling material formed by rainfall confluence. Limited by a certain slope groove, it pours down to the accumulation area at a significantly fast speed, which can cause great damage to a variety of buildings along the way [2].
Debris flow disasters are extensively distributed around the world, while China is a mountainous country and the mountainous area accounts for approximately 69% of the total land area. e Tianshan Mountains, Kunlun Mountains, Himalayas, Qinling Mountains, Hengduan Mountains, and Changbai Mountains in China are all high-risk areas of debris flow. e total area of debris flow distribution is 4.3 million km 2 , of which 1.3 million km 2 is an intense active area [3]. In China, there have been many major debris flow disasters, among which, the most notable one is the August 8 Zhouqu debris flow disaster in Gansu Province, with the severe consequences of 1,481 deaths, 1,824 injuries, 284 missing persons, and nearly 20,000 people affected [3]. e data show that debris flow composed of particles and fluids has great kinetic energy and impulsive force and can cause devastating damage to the structure [4]. Debris flow study is relatively complex, and previous studies are generally simplified into pure particle matter and fluid. For example, Ashwood and Hungr [5] used quartz sand and rubble for performing the physical simulation test to study the impact of debris flow. Ng [6] et al. studied the impact performance of dry sand and viscous fluid by the conducting centrifugal test. However, some scholars have studied the mixture of solid-liquid debris flow. Haas [7] et al. studied the solid-liquid mixture composed of rubble, sand, clay, and water. e results proved that the performance of solid-liquid mixture was significantly different from that of single-phase flow. e interaction between solid particles and fluid of debris flow would significantly change the movement characteristics and impulsive force of debris flow. It is challenging to describe and quantify the formation, movement, impact, and deposition process of solid-liquid two-phase debris flow due to solid-liquid interaction, which has always been the difficulty of study.
In addition to the formation, movement, impact, and deposition process of debris flow, the prevention and control measures of debris flow have also dominated a vital role in the protection of life and property [8][9][10]. As for the design of debris flow retaining structure, the impact of debris flow should be considered; otherwise, structural damage will be caused and the safety of life and property downstream will be seriously threatened [11][12][13]. e retaining structure can generally be divided into rigid retaining structure and flexible retaining structure [14], and the peak impulsive force of debris flow plays a key role in the design of retaining dam. At present, it is generally based on the fluid statics model, fluid dynamics model, or hybrid model [15]. In particular, the fluid dynamics model has been widely applied in engineering design due to the reliability of data. However, because of the different conditions, the value of empirical coefficient varies greatly among different scholars, seriously affecting the accurate prediction of the peak impulsive force of debris flow [16].
In recent years, the flexible retaining dam becomes increasingly popular because of its simple installation, high construction efficiency, good economy, and small impact on the environment [17]. Compared with the rigid retaining structure, the flexible retaining structure can generate large deformation under the impact, resulting in a better energy dissipation mechanism, thereby lowering the peak impulsive force [18]. However, the dynamic coupling among the particle-fluid flexible retaining structure is significantly complex, which depends on the hydrodynamic mechanism, the particle movement characteristics, and the deformation and mechanical characteristics of the flexible retaining structure, as well as the energy conversion and dissipation between them. It remains a highly challenging topic. At present, there is no reasonable theoretical basis for the design of flexible retaining structure. erefore, a limited test and numerical simulation are carried out and the engineering is only designed and constructed by experience. Although the failure analysis of flexible retaining dam under debris flow impact is rarely discussed, this is the key issue of engineering design [19].
One of the main aspects of the study of the coupling effects of debris flow and flexible support structure is the numerical simulation method, which is mainly divided into continuous environment calculation method, discrete environment method, and mixed environment method [17,20,21]. e continuous environment calculation method relies on the continuous distribution of mass and deformation of the object under study. e equations describing the dynamic characteristics are established. ese equations are discretized on the computational domain based on grids or particles and solved by combining the initial and boundary conditions [18,22]. e main continuous medium calculation methods are the finite element method and the smooth particle kinematics method, or in other words, both debris flow and flexible retaining structures are simulated by the finite element method or the SPH method. is method is more mature and has been verified in theory and experiment. e discrete medium method in the analysis of debris flow and flexible retaining structure is mainly the discrete element method [19,23]. Leonardi et al. [24] employed the DEM method to study the coupling effect of debris flow and flexible retaining structure and compared with the experimental results. e results showed that the DEM method could better reflect the interaction between debris flow and flexible retaining structure and was close to the experimental results, thus verifying the effectiveness of numerical analysis method.
In order to make full use of the advantages of different methods, the coupling numerical analysis method of debris flow and retaining structure has become a research hotspot in recent years. Leonardi et al. [25][26][27] used the LBM-DEM method to address the effect of the debris flow on the structure and adopted the coupling numerical analysis method to conduct a series of analysis and numerical verification, proving the robustness of the method in various cases. e debris flow impact of different flexible retaining structures was carried out, and the failure mode of the flexible retaining structure was explored. Li and Zhao [28,29] adopted the coupled CFD-FEM method to simulate the debris flow, as well as the impact on rigid and flexible retaining dams. e model revealed that the actual simulation of the debris flow needs to consider the solid-liquid interaction. In addition, the effect of erosion on the impulsive force of debris flow was taken into account. e results showed that the debris flow impulsive force was significantly influenced by the slope and the solid-liquid ratio. e impact law and failure mode of flexible retaining dam are studied by employing the coupling numerical analysis method. Liu et al. [30,31], respectively, used SPH-FEM and DEM-FEM to investigate the impact characteristics of debris flow, large stone, and retaining dam as well as the coupling effect of debris flow and flexible retaining structure. e results showed that the coupling numerical method could well reproduce the impact-climbing-silting process of debris flow and the impact issue of flexible retaining structure. Li et al. [32] used the SPH-DEM method to study the dynamic performance and impact force of debris flow on a rigid retaining structure. Although the coupling numerical analysis method provides the basis and reference for the theory and method of this paper, there is no research and discussion on the mechanical properties and impulsive force of flexible retaining structure. e flow process of debris flow is shown in Figure 1 [3].
In summary, there is no reasonable theoretical basis for the design of flexible retaining structure. Instead, the design and construction of the project only depend on experience, and no corresponding specification for the prevention and control of debris flow by flexible retaining dams has been formulated in various countries. ere are very few research studies on dynamic response analysis of flexible retaining structure impacted by debris flow. Currently, the research only concentrates on the retaining structure affected by large rock, debris flow fluid, or debris flow fluid along with large rock. Only few numerical analysis methods are used to analyze the coupling effect of debris flow and flexible retaining structure. e SPH-DEM-FEM numerical analysis method is applied to investigate the dynamic response of a retaining structure impacted by debris flow here, providing a certain reference for debris flow prevention engineering practice.

Numerical Calculation Theory
LS-DYNA code is adopted for SPH-DEM-FEM numerical simulation, in which SPH is used to simulate liquid debris flow, DEM is used to simulate solid particles, and FEM is used to simulate flexible retaining structure and debris flow channel. e SPH-DEM-FEM coupling analysis method is adopted to establish a complex dynamic interaction model of particle-fluid-structure, which can truly simulate the dynamic interaction between actual debris flow and flexible retaining structure. e realization principle and control equations of numerical simulation are introduced in this chapter.

DEM Control Equation.
Debris flow particles in the solid phase are modelled by DEM. An elastic model is adopted for the contact between solid-phase debris flow particles [33][34][35][36].
e expression of force on particles is shown as follows: Among which, λ and δ represent the overlapping distances of normal direction and tangential direction, respectively, and k n and k t denote normal stiffness and tangential stiffness, respectively. e resulting external force and the resulting external moment on the particles are expressed as follows: Among which, f s represents the friction between particles.
An elastic material model is adopted for the interaction force between particles, and the constitutive model of the material is shown as follows: Among which, E e represents the particle elastic model and ε e is the particle elastic strain.

SPH Control Equation.
In this paper, SPH is used to simulate the flow of debris in the liquid phase. e basic concept of this method is to regard the fluid (or solid) as continuous one and describe it by interacting particle groups. Each particle point carries various physical parameters, including mass and velocity [37][38][39][40].
e approximate function of particles is shown in the following formula: where W is the kernel function. Kernel function W can be obtained by function θ, and their relation is shown in the following formula: where d is the spatial dimension and h is the smooth length that varies with time and space. In this paper, the smooth kernel function is shown as follows: where C is a generalized constant depending on the spatial dimension. e SPH method bases on the integral formula of moving particles ((x i (t))), i ⊂ 1, 2, . . . , N { }, among which x i (t) is the spatial position vector of particle i, which moves along the velocity field v.

Advances in Civil Engineering
At this time, the approximation function of particle can be expressed by the following formula: wherew j � (m j /ρ j ), and the gravity of the particle w j changes in proportion to the divergence of the fluid.

SPH-DEM-FEM Coupling Control Equation.
e DEM-FEM coupling equation is used to describe the complex dynamic interaction between solid-phase particles and a flexible containment structure, which can be found in the following equation [18,[41][42][43]: where f n,ij and f t,ij represent the normal and the tangential contact force. e DEM-FEM coupling equation is shown as follows: e SPH-DEM coupling control equation is used to describe the complex dynamic interaction between the particles in the solid phase of the debris flow and the fluid in the liquid phase of the debris flow [18], as presented in the following equation: where

Calculation Model
is paper selects Nanjiaogou as the prototype. e rugged topography in the ditch is sharp, the gullies can be found anywhere, and the slope of the main ditch formation area is 19°∼47°, which belong to a vast debris flow ditch. e detailed geological information of Nanjiaogou can be found in [32].
In order to study the impulsive force and dynamic response of flexible retaining structure impacted by debris flow, a complex dynamic interaction model of particle-fluid-structure is established in the present study. Considering that the debris flow impulse force and the dynamic response of the flexible retaining structure are mainly studied in this paper and the actual gully is basically in trapezoidal distribution, to facilitate the study, the debris flow gully is simplified into trapezoidal one, and retaining and governance measures are performed on the main gully formation area with the steep slope. Considering that the slope of the main ditch formation area is 19°∼47°, the average slope of the main gully formation area is confirmed at 30°. In this paper, the impulsive force and dynamic response of flexible retaining structure are analyzed, when the slope of gully is 15°, 20°, 30°, 40°, and 50°.

Geometric Model.
e geometric model of the debris flow is shown in Figure 2, and the geometric model of the flexible retaining structure is presented in Figure 3. e specific parameters of the geometric model are displayed in Table 1. e parameters such as debris flow channel and DEM particle density, elastic modulus, and Poisson's ratio are shown in the literature [32], and the specific mechanical parameters of the flexible retaining net and support cable are shown in the literature [24,28,29].
e research object of this paper is mainly the impact force and dynamic response of the flexible retaining structure. en, according to the field survey, the channel is roughly trapezoidal distribution and the bottom and the side of the gully are then simplified into the trapezoidal section according to the on-site situation. Furthermore, in view of the fact that the slope of the main ditch formation area is 19°∼47°, the average value 30°is selected; debris flow is composed of solid particles and fluids, which impacts the flexible retaining structure under the action of gravity, reproducing scenery of the impact, blocking by the flexible retaining structure, climbing along the flexible retaining structure, and water-rock separating in the process of impacting by debris flow. And this paper studies the impulsive force and dynamic response of flexible retaining structure impacted by solid-liquid two-phase debris flow. e flexible retaining structure consists of flexible retaining nets and support cables. e four support cables are numbered as NO.1, NO.2, NO.3, and NO.4 from the bottom to top. e distance between the support cables is 0.5 m. e geometric model of the flexible retaining structure is shown in Figure 3. With a diameter of 20 mm, the flexible retaining net is formed by steel wires winding. e hexagonal side length of the retaining net is 0.1 m, which can smoothly discharge debris flow and effectively prevent debris flow particles from flowing out. e support cable is formed by a steel wire winding with a diameter of 20 mm, having an equivalent diameter of 0.1 m. e flexible retaining net and the supports cable are fixedly connected to bear the force in a better way. e specific mechanical parameters of the flexible retaining net and support cable are shown in the literature [24,28,29].

Constitutive Model of Materials.
Because the channel rock mass in this paper relatively remains stable, the debris flow channel is simulated by rigid shell material. In this paper, the RIGID material model is adopted, with 6056 shell elements in total. e parameters of debris flow channel adopt the numerical simulation parameters in the literature [19].
According to the field sampling test, the measured average density of solid particles in debris flow source is 2650 kg/m 3 . Based on the field survey and laboratory test, we measured that the particle size range of provenance is 0.2∼0.4 m. e size range of DEM particles set by us is consistent with the field test results. e solid-liquid ratio of debris flow exerts a great influence on the characteristics of debris flow. In this paper, the ratio of solid to liquid is 1:1.
is paper uses discrete element particle simulation. e LS-DYNA ELASTIC model is adopted for the material constitutive, the LS-DYNA DE_BY_PART contact model is adopted for the interaction between the particles, and the LS-DYNA DE_SURFACE_COUPLING model is used for the contact between the particles and the channel. In addition, the LS-DYNA DE_TO_BEAM_COUPLING model is used for the contact between particles and flexible retaining structure, with a total of 2236 particles in the abovementioned models. e object of this paper is the low-speed moving dilute debris flow, without considering the effect of viscous force for the time being. e density of water is 1000 kg/m 3 , and the friction force between fluids is 0.12. e value of friction coefficient will be described in detail in Section 3.3. e SPH model is used for analysis, and the constitutive model employs the NULL material model, with the viscosity coefficient of 0.001. e constitutive model of water is involved in this paper, while the constitutive model uses the EOS_-MURNAGHAN model, and the constant of state equation GAMMA is 7.0, with K 0 of 150,000. To sum up, there are totally 4351 particles in the model. e SPH_DE_COU-PLING contact model is adopted for the interaction between SPH particles and DEM particles, and the AUTOMA-TIC_NODES_TO_SURFACE model is adopted for the contact between SPH particles and channels. e parameters of fluids adopt the numerical simulation parameters in the literature [19].

Boundary Conditions.
e channel and the flexible retaining structure are firmly connected with the ground. e bottom and both sides of the debris channel are completely fixed and capped. e flexible retaining net and the supporting cable are connected with the ground through the anchor cable. erefore, the flexible retaining net and support cables are fully fixed in this paper. e coefficient of friction between debris flow solids is typically 0.1-0.5 [16][17][18] by analyzing and comparing the experimental and simulation results in a large amount of relevant literature. According to similar studies, the trialand-error method is generally used to explore and compare the numerical simulation results with actual cases. When the peak velocity data of debris flow numerical simulation are close to the experimental or actual data, the parameters are regarded to be more accurate. After adopting the trial-anderror method and reviewing similar literature, it can be found that the simulated results are close to the actual results, when the coefficient of friction between the debris flow solids and the channel and the retaining wall is 0.3 and the coefficient of friction between the debris flow fluid and the channel and the flexible retaining structure is 0.12.
We use the dynamic explicit method of LS-DYNA code for performing the analysis. e analysis step time is 10 s. e debris flow begins to slide with an initial velocity of 0 m/ s 2 . e gravitational acceleration of both the debris flow fluid and the debris flow solid particles is 9.8 m/s 2 . Totally, the calculation gets 1974743 incremental steps. e CPU time is 22 h 3 min 43 s.

Impacting Process.
e impact, retaining, water-rock separation, and deposition process of flexible retaining structure impacted by debris flow are shown in Figure 4. When t � 0.9 s, the debris flow starts to move, and at this time, "dragon head, dragon body, and dragon tail" have basically taken shape. Debris flow fluid and solid particles yield complex dynamic coupling interaction. When t � 2.3 s, the solid-liquid two-phase debris flow reaches the bottom of the flexible retaining structure. At this moment, the maximum velocity of debris flow reaches 10.52 m/s, and the particle-fluid-structure joint is forced to produce complex dynamic interaction and impact the flexible retaining structure. Meanwhile, the "dragon head" of debris flow impacts the flexible retaining structure, climbing high and separating water-rock, and the solid particles remain at the bottom of the retaining structure under the action of retaining structure. At the same time, most of the fluid is discharged from the retaining structure. When t � 3.5 s, the "dragon body" impacts the flexible retaining structure. Most of the debris flow solid particles are retained at the bottom, and the fluid is discharged. e debris flow will continue to impact the flexible retaining structure. At this time, because of the climbing effect, the bottom of the flexible retaining structure produces a dead load due to particle deposition, while the top part is still impacted by the dynamic impact of the debris flow. e SPH method can greatly deal with the splash and diffusion of fluid particles and participate in complex dynamic interaction with flexible retaining structure and debris flow particles. e numerical simulation well reflects the complex dynamic interaction of particle-fluidstructure. When t � 5.2 s, the debris flow silts back and the particles finally deposit at the bottom of the retaining dam to produce a dead load. e fluid is slowly discharged.

Climbing and Silting Process.
e climbing and silting process of the flexible retaining structure impacted by debris flow is illustrated in Figure 5. When t � 2.5 s, debris flow particles reach the bottom of the flexible retaining structure and start to impact the debris flow retaining dam and climb high. At this time, the water-rock separation has already taken shape and the fluid is discharged through the retaining structure. e flexible retaining structure is deformed. Besides, the most deformed parts are the middle and bottom of the flexible retaining structure. When t � 3.0 s, debris flow particles continue to climb high and some even silt back in certain part. At this time, it can be observed that particlefluid-structure shapes complex dynamic interaction, and the current climbing height is 1.3 m. When t � 3.5 s, the "dragon body" impacts the retaining dam, and the debris flow particles climb to the highest level, with obvious back silting and particle splashing. e climbing height is as high as 2.1 m. When t � 5.2 s, the debris flow silts back, finally producing a dead load at the bottom of the retaining dam. At this moment, the height of silting is the height of the retaining structure.  [24][25][26][27], showing that the above method is suitable for simulating the flexible retaining structure impacted by debris flow.

Analysis of Impulsive Force Change over Time.
e impulsive force change of flexible retaining structures with different slopes (15°, 20°, 30°, 40°, and 50°) over time is displayed in Figure 7. It can be observed from the figure that when the slope is 50 degrees and t � 2.9 s, the impulsive force of the debris flow on the retaining dam reaches its peak value, which is 6089 kN, and then decreases rapidly to 3832 kN. When the slope is 15°, 20°, 30°, and 40°, the peak values of impulsive force of the retaining dam are 987 kN, 1983 kN, 3690 kN, and 4625 kN, respectively, and subsequently all decrease to dead loads. With the increase of the slope, the start time of the debris flow impact force is 3.4 s, 3.2 s, 3 s, 2.5 s, and 2 s, respectively. In addition, the peak time of the debris flow impact force is 4.9 s, 4.7 s, 3.9 s, 3.2 s, and 2.9 s, respectively.
We compare the numerical calculation results according to the semiempirical formula of impulsive force peak value of the debris flow based on hydrodynamic theory proposed by Kwan [44]: where κ � aFr b is the dynamic pressure coefficient, ρ represents the density, v is the impulsive velocity, w is the width, and h 0 is the height of debris flow.  Resultant Velocity 1.051e+01 9.459e+00 8.408e+00 7.357e+00 6.306e+00 5.255e+00 4.204e+00 3.153e+00 2.102e+00 1.051e+00 0.000e+00

(c)
Resultant Velocity 6.842e+00 6.158e+00 5.474e+00 4.789e+00 4.105e+00 3.421e+00 2.737e+00 2.053e+00 1.368e+00 6.842e-01 0.000e+00 (d) Advances in Civil Engineering e calculation results are obtained from the numerical calculation method in this paper, and the existing empirical formulas are shown in Table 2. e calculation error ranges from 12.1% to 27.4%. Comparing the numerical calculation method adopted in the present study with the empirical values, the impulsive force results of the two methods are relatively close. However, the empirical methods do not consider the water-rock separation and the complex dynamic interaction of particle-fluid-structure. erefore, the numerical analysis method is more reasonable than the empirical methods.

Maximum Deformation of Flexible Retaining Structure.
e maximum deformation of flexible retaining structure with different slopes (15°, 20°, 30°, 40°, and 50°) can be found in Figure 8. e maximum deformation of flexible retaining structure is positively associated with the channel slope, and the maximum deformation position is located in the middle of flexible retaining structure, which is consistent with the observed results in practice. When the slope is 15°, 20°, 30°, 40°, and 50°, the maximum peak displacement values in the middle of the flexible retaining structure are 0.0283 m,

Discussion
According to the results, the SPH-DEM-FEM coupling numerical analysis method vividly simulates the complex dynamic interaction of particle-fluid-structure, reproducing the impact, retaining, water-stone separation, climbing, silting, and deposition process of debris flow. e SPH method can greatly handle the splashing and diffusion of fluid particles and participate in complex dynamic interaction with flexible retaining structure and debris flow particles. e water-stone separation phenomenon can well reproduce the real debris flow blocking process, which is difficult to be reflected by other numerical analysis methods, which can also build up a good foundation for better evaluating the accuracy of debris flow impact assessment of flexible retaining structure. e displacement change of flexible retaining structure over time is analyzed here. Moreover, the simulation results and laws are consistent with the existing research results. e analysis of impulsive force change over time is carried out. Although the analysis results are similar to the existing empirical formula, the empirical formula does not consider the water-rock separation and the complex dynamic interaction between particle-fluid-structure, and the error between the numerical simulation results and the empirical formula ranges from 12.1% to 27.4%. erefore, the numerical analysis method is more reasonable than the empirical method. e peak impulsive force of debris flow is positively correlated with the slope of the channel, and the   Advances in Civil Engineering larger the slope is, the greater the peak impulsive force of debris flow is. en, the impulsive force of debris flow is decreased to a dead load, with enormous decrease in amplitude. When the slope is small, the peak impulsive forces of debris flow are close to a dead load. e impulsive force of debris flow will decrease slowly in the later stage, which is caused by the fluid removal from the flexible retaining structure due to the water-rock separation. It cannot be presented by employing the empirical formula method and traditional numerical analysis method. e results of the coupled numerical analysis method are more convincing than the empirical formula method. e maximum deformation of flexible retaining structure is analyzed, and the deformation is positively associated with the slope of the channel. is is because when the slope of the channel is sharp, the gravitational potential energy of debris flow is relatively large, and once the debris flow takes shape, the kinetic energy and impulsive force produced by debris flow will be far greater than when the slope is small. Li et al. [32] explore the impulsive force and dynamic performance of rigid retaining structure impacted by debris flow by adopting the SPH-DEM-FEM method. By comparing impulsive force and dynamic performance of flexible retaining structure impacted by solid-liquid two-phase debris flow discussed in this paper, the dynamic performance and impact characteristics of flexible retaining structure and rigid retaining structure are quite different. Because of the water-rock separation and the remarkable ability to absorb energy, the impulsive force of the flexible retaining structure is much smaller than that of the rigid retaining structure.
is also indicates that the traditional impulsive force formula of debris flow is not applicable to the flexible retaining structure, which is also the research significance of the present study. e results show that the mechanical performance of flexible retaining structure is more reasonable than that of rigid retaining structure.
Since the rigidity of the flexible retaining structure is much smaller than that of the rigid retaining structure, the dynamic characteristics between them are significantly different. e displacement change over time between the flexible retaining structure and rigid retaining structure is different by two orders of magnitude. e recovery ability after deformation of the flexible retaining structure is far less than that of the rigid retaining structure, requiring us to accordingly monitor the deformation of the flexible retaining structure in order to prevent dam collapse from bringing devastation to the downstream. e maximum deformation position is located in the middle of the flexible retaining structure, which gradually decreases to both sides. e support cable plays an obvious role in restraining the flexible retaining net and acts as a major carrier in bearing and transmitting force. e results of numerical simulation are compared with the deformation of real flexible retaining structure, showing good consistence.
In this article, the terrain conditions are simplified, considering neither the impact of the unevenness of the movement path on the motion characteristics of the debris flow and the dynamic characteristics of the flexible retaining structure nor the impact of the loose accumulation layer on the movement path on the impulsive force of the debris flow and the dynamic characteristics of the flexible retaining structure. All the abovementioned problems are the research difficulties and can be further studied in future research.

Conclusion
(1) Based on the coupling numerical method, this article considers the complex dynamic interaction of particle-fluid-structure, reproducing the whole process of flexible retaining structure impacted by debris flow, such as impact, retaining, water-rock separation, climbing, siltation, and deposition. However, the empirical formula does not consider the water-rock separation and the complex dynamic interaction of particle-fluid-structure, and thus, the error of empirical formula results and the numerical simulation results ranges from 12.1% to 27.4%. As a result, the numerical analysis method is more reasonable than the empirical method. (4) We analyze the maximum deformation of the flexible retaining structure, concluding that the deformation is positively associated with the slope of the channel. e maximum deformation is located in the middle of the flexible retaining structure, and the support cable can obviously restrain the flexible retaining net.
Data Availability e data for this study are available on request to the corresponding author.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.