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

Li, Bailong; Wang, Changming; Li, Yanying; Zhang, Shuhua

doi:https://doi.org/10.1155/2021/9098250

Advances in Civil Engineering

On this page

Abstract Introduction Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Assessment and Prevention Methods of Collapse Disasters

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 9098250 | https://doi.org/10.1155/2021/9098250

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

Bailong Li,¹Changming Wang,¹Yanying Li,¹and Shuhua Zhang²

Academic Editor: Xiaohu Zhang

Received17 Jun 2021

Accepted11 Oct 2021

Published23 Oct 2021

Abstract

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. The 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. The 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. The dynamic time-history curve of the coupling numerical analysis method for the flexible retaining structure is consistent with the results of the existing literature. The 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. The 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.

1. 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. The Tianshan Mountains, Kunlun Mountains, Himalayas, Qinling Mountains, Hengduan Mountains, and Changbai Mountains in China are all high-risk areas of debris flow. The total area of debris flow distribution is 4.3 million km², of which 1.3 million km² 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]. The 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. The results proved that the performance of solid-liquid mixture was significantly different from that of single-phase flow. The 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–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–13]. The 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. Therefore, 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]. The continuous environment calculation method relies on the continuous distribution of mass and deformation of the object under study. The equations describing the dynamic characteristics are established. These equations are discretized on the computational domain based on grids or particles and solved by combining the initial and boundary conditions [18, 22]. The 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. This method is more mature and has been verified in theory and experiment. The 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. The 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–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. The 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. The 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. The results showed that the debris flow impulsive force was significantly influenced by the slope and the solid-liquid ratio. The 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. The 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. The 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. There 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. The 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.

2. 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. The 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. The realization principle and control equations of numerical simulation are introduced in this chapter.

2.1. 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–36].

The expression of force on particles is shown as follows:

Among which, and represent the overlapping distances of normal direction and tangential direction, respectively, and and denote normal stiffness and tangential stiffness, respectively.

The resulting external force and the resulting external moment on the particles are expressed as follows:

Among which, 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, represents the particle elastic model and is the particle elastic strain.

2.2. SPH Control Equation

In this paper, SPH is used to simulate the flow of debris in the liquid phase. The 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–40].

The approximate function of particles is shown in the following formula:where is the kernel function.

Kernel function 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.

The SPH method bases on the integral formula of moving particles , , among which is the spatial position vector of particle , which moves along the velocity field .

At this time, the approximation function of particle can be expressed by the following formula:where, and the gravity of the particle changes in proportion to the divergence of the fluid.

2.3. SPH-DEM-FEM Coupling Control Equation

The 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–43]:where and represent the normal and the tangential contact force.

The DEM-FEM coupling equation is shown as follows:

The 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 and are the force on each other and , , and denote the torque.

3. Calculation Model

This paper selects Nanjiaogou as the prototype. The 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. The 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°.

3.1. Geometric Model

The 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. The specific parameters of the geometric model are displayed in Table 1. The 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].

The research object of this paper is mainly the impact force and dynamic response of the flexible retaining structure. Then, 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.

The flexible retaining structure consists of flexible retaining nets and support cables. The four support cables are numbered as NO.1, NO.2, NO.3, and NO.4 from the bottom to top. The distance between the support cables is 0.5 m. The 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. The 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. The support cable is formed by a steel wire winding with a diameter of 20 mm, having an equivalent diameter of 0.1 m. The flexible retaining net and the supports cable are fixedly connected to bear the force in a better way. The specific mechanical parameters of the flexible retaining net and support cable are shown in the literature [24, 28, 29].

3.2. 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. The 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³. Based on the field survey and laboratory test, we measured that the particle size range of provenance is 0.2∼0.4 m. The size range of DEM particles set by us is consistent with the field test results. The 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. This paper uses discrete element particle simulation. The 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.

The object of this paper is the low-speed moving dilute debris flow, without considering the effect of viscous force for the time being. The density of water is 1000 kg/m³, and the friction force between fluids is 0.12. The value of friction coefficient will be described in detail in Section 3.3. The SPH model is used for analysis, and the constitutive model employs the NULL material model, with the viscosity coefficient of 0.001. The 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₀ of 150,000. To sum up, there are totally 4351 particles in the model. The SPH_DE_COUPLING contact model is adopted for the interaction between SPH particles and DEM particles, and the AUTOMATIC_NODES_TO_SURFACE model is adopted for the contact between SPH particles and channels. The parameters of fluids adopt the numerical simulation parameters in the literature [19].

The support cable of flexible retaining structure is modelled by cable element, with the density of 7850 kg/m³ and yield strength of 1770 MPa. Because the flexible retaining structure bears the dynamic impact, the material constitutive model adopts the PLASTIC_KINEMATIC model, with the ultimate strain of 0.05. The flexible retaining net is simulated by BEAM element, with the density of 7850 kg/m³ and yield strength of 1380 MPa. Because the flexible retaining structure bears the dynamic impact, the constitutive model of the material adopts the PLASTIC_KINEMATIC model. The specific mechanical parameters of flexible retaining nets and support cables can be found in the literature [24, 28, 29].

3.3. Boundary Conditions

The channel and the flexible retaining structure are firmly connected with the ground. The bottom and both sides of the debris channel are completely fixed and capped. The flexible retaining net and the supporting cable are connected with the ground through the anchor cable. Therefore, the flexible retaining net and support cables are fully fixed in this paper.

The coefficient of friction between debris flow solids is typically 0.1–0.5 [16–18] by analyzing and comparing the experimental and simulation results in a large amount of relevant literature. According to similar studies, the trial-and-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-and-error 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. The analysis step time is 10 s. The debris flow begins to slide with an initial velocity of 0 m/s². The gravitational acceleration of both the debris flow fluid and the debris flow solid particles is 9.8 m/s². Totally, the calculation gets 1974743 incremental steps. The CPU time is 22 h 3 min 43 s.

4. Establishment of Calculation Model

4.1. Impacting Process

The 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. The 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. The 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. The numerical simulation well reflects the complex dynamic interaction of particle-fluid-structure. 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. The fluid is slowly discharged.

(a)

(b)

(c)

(d)

4.2. Climbing and Silting Process

The 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. The 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 particle-fluid-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. The 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.

(a)

(b)

(c)

(d)

4.3. Analysis of Displacement Change of Flexible Retaining Structure over Time

The displacement change over time in the middle position of four support cables (NO.1, NO.2, NO.3, and NO.4) of flexible retaining structure under different slopes is monitored in Figure 6. The results reveal that when the slope is 50 degrees and t = 3 s, the displacement of the middle position of NO.4 support cable of flexible retaining structure reaches the maximum, and the maximum value is 0.265 m. Then, it decreases rapidly to 0.160 m after reaching the peak value. When the slope is 15 degrees, 20 degrees, 30 degrees, and 40 degrees, the peak displacement value of the middle position of the support cable of the flexible retaining structure is 0.0283 m, 0.071 m, 0.164 m, and 0.221 m, respectively. Subsequently, all decrease to dead displacement, which are 0.0146 m, 0.055 m, 0.122 m, and 0.170 m, respectively. When the slope is 15°, 20°, 30°, 40°, and 50°, the difference between the maximum peak displacement and static displacement of the flexible barrier structure is 0.0137 m, 0.016 m, 0.042 m, 0.051 m, and 0.105 m, respectively. Moreover, the simulation results are identical to the numerical simulation and test results performed by Leonardi et al. [24–27], showing that the above method is suitable for simulating the flexible retaining structure impacted by debris flow.

(a)

(b)

(c)

(d)

(e)

4.4. Analysis of Impulsive Force Change over Time

The 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 is the dynamic pressure coefficient, represents the density, is the impulsive velocity, is the width, and h₀ is the height of debris flow.

The calculation results are obtained from the numerical calculation method in this paper, and the existing empirical formulas are shown in Table 2. The 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. Therefore, the numerical analysis method is more reasonable than the empirical methods.

4.5. Maximum Deformation of Flexible Retaining Structure

The maximum deformation of flexible retaining structure with different slopes (15°, 20°, 30°, 40°, and 50°) can be found in Figure 8. The 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, 0.071 m, 0.164 m, 0.221 m, and 0.265 m, respectively. Clearly, it can be seen that the support cable can restrain the flexible retaining net from deformation. Because the flexible retaining net and the supporting cables are fixed to each other, the deformation of the two is consistent with each other. In addition, the deformation of the flexible retaining net between the supporting cables is bigger than the deformation of the supporting cable, showing a convex shape.

5. 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. The 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. The 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. The displacement change of flexible retaining structure over time is analyzed here. Moreover, the simulation results and laws are consistent with the existing research results.

The 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%. Therefore, the numerical analysis method is more reasonable than the empirical method. The peak impulsive force of debris flow is positively correlated with the slope of the channel, and the larger the slope is, the greater the peak impulsive force of debris flow is. Then, 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. The 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. The results of the coupled numerical analysis method are more convincing than the empirical formula method.

The maximum deformation of flexible retaining structure is analyzed, and the deformation is positively associated with the slope of the channel. This 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. This 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. The 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. The displacement change over time between the flexible retaining structure and rigid retaining structure is different by two orders of magnitude. The 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. The maximum deformation position is located in the middle of the flexible retaining structure, which gradually decreases to both sides. The support cable plays an obvious role in restraining the flexible retaining net and acts as a major carrier in bearing and transmitting force. The 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.

6. 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.(2)We compare the dynamic response analysis of flexible retaining structure impacted by debris flow under different slopes, monitoring the displacement change of the middle position of four support cables (NO.1, NO.2, NO.3, and NO.4) of flexible retaining structure over time. It can be found that the simulation results and laws are consistent with the existing research results.(3)We compare the change law of the impulsive force of debris flow under different slopes. 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. The 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

The data for this study are available on request to the corresponding author.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was funded by the Natural Science Foundation of China (grant no. 41972267).

References

J. L. Miao, W. Z. Zhang, and J. Y. Zhou, “Numerical simulation of the accumulation state of viscous debris flow by smooth particle hydrodynamics method,” Journal of Natural Disasters, vol. 22, no. 6, pp. 125–130, 2013.
View at: Google Scholar
M. Jakob, O. Hungr, and D. M. Jakob, “Debris-flow hazards and related phenomena,” Environmental and Engineering Geoscience, vol. 8, no. 1, pp. 165–171, 2007.
View at: Google Scholar
C. A. Dowling and P. M. Santi, “Debris flows and their toll on human life: a global analysis of debris-flow fatalities from 1950 to 2011,” Natural Hazards, vol. 71, no. 1, pp. 203–227, 2014.
View at: Publisher Site | Google Scholar
C. Tang, N. Rengers, T. W. J. van Asch, Y. H. Yang, and G. F. Wang, “Triggering conditions and depositional characteristics of a disastrous debris flow event in Zhouqu city, Gansu province, northwestern China,” Natural Hazards and Earth System Sciences, vol. 11, no. 11, pp. 2903–2912, 2011.
View at: Publisher Site | Google Scholar
W. Ashwood and O. Hungr, “Estimating total resisting force in flexible barrier impacted by a granular avalanche using physical and numerical modeling,” Canadian Geotechnical Journal, vol. 53, no. 10, pp. 1700–1717, 2016.
View at: Publisher Site | Google Scholar
C. W. W. Ng, D. Song, C. E. Choi, L. H. D. Liu, and J. S. H. Kwan, “Impact mechanisms of granular and viscous flows on rigid and flexible barriers,” Canadian Geotechnical Journal, vol. 54, no. 2, pp. 188–206, 2016.
View at: Publisher Site | Google Scholar
T. Haas, L. Braat, J. R. Leuven, I. R. Lokhorst, and M. G. Kleinhans, “Effects of debris flow composition on runout, depositional mechanisms, and deposit morphology in laboratory experiments,” Journal of Geophysical Research: Earth Surface, vol. 120, no. 9, pp. 19493–1972, 2015.
View at: Google Scholar
X. B. Zhou, Debris Flow’s Harm to Bridge Engineering and its Prevention, Tongji University, Shanghai, China, 2006.
Z. Tao, M. Li, C. Zhu, M. He, X. Zheng, and S. Yu, “Analysis of the critical safety thickness for pretreatment of mined-out areas underlying the final slopes of open-pit mines and the effects of treatment,” Shock and Vibration, vol. 2018, no. 2, pp. 1–8, 2018.
View at: Publisher Site | Google Scholar
X. J. Yang, J. M. Wang, C. Zhu, M. C. He, and Y. Gao, “Effect of wetting and drying cycles on microstructure of rock based on SEM,” Environmental Earth Sciences, vol. 78, no. 6, pp. 1866–6280, 2019.
View at: Publisher Site | Google Scholar
Z. Liang, C. Wang, D. Ma, and K. U. J. Khan, “Exploring the potential relationship between the occurrence of debris flow and landslides,” Natural Hazards and Earth System Sciences, vol. 21, no. 4, pp. 1247–1262, 2021.
View at: Publisher Site | Google Scholar
L. Ban, C. Zhu, C. Qi, and Z. Tao, “New roughness parameters for 3D roughness of rock joints,” Bulletin of Engineering Geology and the Environment, vol. 78, no. 6, pp. 4505–4517, 2019.
View at: Publisher Site | Google Scholar
C. Zhu, K. Zhang, H. Cai et al., “Combined application of optical fibers and CRLD bolts to monitor deformation of a pit-in-pit foundation,” Advances in Civil Engineering, vol. 2019, no. 1, Article ID 2572034, 16 pages, 2019.
View at: Publisher Site | Google Scholar
A. Leonardi, F. K. Wittel, M. Mendoza, and J. Herrmann, “Multiphase debris flow simulations with the discrete element method coupled with a lattice-Boltzmann fluid,” in Proceedings of the III International Conference on Particle-based Methods Fundamentals and Applications, pp. 276–287, Stuttgart (Germany), 2013.
View at: Google Scholar
A. Armanini, “On the dynamic of debris flows,” Recent Development on Debris Flows, Springer, Heidelberg, Germany, 1997.
View at: Google Scholar
C. Lichtenan, Die berechnung von sperren in beton und eisenbeton, kolloquium on torrent dams odc 384.3, Mitteilungen der Forstlichen Bundes-Versuchsanstalt, Wien, 1973.
A. Albaba, S. Lambert, F. Kneib, B. Chareyre, and F. Nicot, “DEM modeling of a flexible barrier impacted by a dry granular flow,” Rock Mechanics and Rock Engineering, vol. 50, no. 11, pp. 3029–3048, 2017.
View at: Publisher Site | Google Scholar
W. Anthony, G. Laurence, and V. Guillaume, “Grains3D, a flexible DEM approach for particles of arbitrary convex shape part I: numerical model and validations,” Powder Technology, vol. 224, pp. 374–389, 2012.
View at: Google Scholar
C. Scheidl, M. Chiari, R. Kaitna et al., “Analysing debris-flow impact models, based on a small scale modelling approach,” Surveys in Geophysics, vol. 34, no. 1, pp. 121–140, 2013.
View at: Publisher Site | Google Scholar
Z. Tao, H. Zhang, C. Zhu, Z. Hao, X. Zhang, and X. Hu, “Design and operation of App-based intelligent landslide monitoring system: the case of Three Gorges Reservoir Region,” Geomatics, Natural Hazards and Risk, vol. 10, no. 1, pp. 1209–1226, 2019.
View at: Publisher Site | Google Scholar
Z. G. Tao, C. Zhu, M. C. He, and K. M. Liu, “Research on the safe mining depth of anti-dip bedding slope in Changshanhao Mine,” Geomechanics and Geophysics for Geo-Energy and Geo-Resources, vol. 36, no. 6, pp. 1–20, 2020.
View at: Google Scholar
Z. Liang, C.-M. Wang, Z.-M. Zhang, and K.-U.-J. Khan, “A comparison of statistical and machine learning methods for debris flow susceptibility mapping,” Stochastic Environmental Research and Risk Assessment, vol. 34, no. 11, pp. 1887–1907, 2020.
View at: Publisher Site | Google Scholar
Z. Tao, Q. Geng, C. Zhu et al., “The mechanical mechanisms of large-scale toppling failure for counter-inclined rock slopes,” Journal of Geophysics and Engineering, vol. 16, no. 3, pp. 541–558, 2019.
View at: Publisher Site | Google Scholar
A. Leonardi, F. K. Wittel, M. Mendoza, R. Vetter, and H. J. Herrmann, “Particle-fluid-structure interaction for debris flow impact on flexible barriers,” Computer-Aided Civil and Infrastructure Engineering, vol. 31, no. 5, pp. 323–333, 2016.
View at: Publisher Site | Google Scholar
A. Leonardi, F. K. Wittel, M. Mendoza, and J. Hans, “Multiphase debris flow simulations with the discrete element method coupled with a lattice-Boltzmann fluid,” in Proceedings of the III International Conference on Particle-based Methods Fundamentals and Applications, pp. 276–287, Stuttgart (Germany), 2013.
View at: Google Scholar
A. Leonardi, F. K. Wittel, M. Mendoza, and H. J. Herrmann, “Lattice-Boltzmann method for geophysical plastic flows,” Recent Advances In Modeling Landslides And Debris Flows, Springer, Cham, Switzerland, 2015.
View at: Publisher Site | Google Scholar
A. Leonardi, F. K. Wittel, M. Mendoza, and H. J. Herrmann, “Coupled DEM-LBM method for the free-surface simulation of heterogeneous suspensions,” Computational Particle Mechanics, vol. 1, no. 1, pp. 3–13, 2014.
View at: Publisher Site | Google Scholar
X. Li and J. Zhao, “A unified CFD-DEM approach for modeling of debris flow impacts on flexible barriers,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 42, no. 14, pp. 1643–1670, 2018.
View at: Publisher Site | Google Scholar
X. Li and J. Zhao, “Dam-break of mixtures consisting of non-Newtonian liquids and granular particles,” Powder Technology, vol. 338, pp. 493–505, 2018.
View at: Publisher Site | Google Scholar
C. Liu, Z. Yu, and S. Zhao, “Quantifying the impact of a debris avalanche against a flexible barrier by coupled DEM-FEM analyses,” Landslides, vol. 17, no. 1, pp. 33–47, 2020.
View at: Publisher Site | Google Scholar
C. Liu, Z. X. Yu, and L. R. Luo, “Dynamic behavior of a concrete dam impacted by debris flows with rock,” Journal of Vibration and Shock, vol. 38, no. 14, pp. 161–168, 2019.
View at: Google Scholar
B. Li, C. Wang, Y. Li et al., “Dynamic response analysis of retaining dam under the impact of solid-liquid two-phase debris flow based on the coupled SPH-DEM-FEM method,” Geofluids, vol. 2020, Article ID 6635378, 12 pages, 2020.
View at: Publisher Site | Google Scholar
H.-k. Chen, H.-m. Tang, and Y.-i. Chen, “Research on method to calculate velocities of solid phase and liquid phase in debris flow,” Applied Mathematics and Mechanics, vol. 27, no. 3, pp. 399–408, 2006.
View at: Publisher Site | Google Scholar
Q.-X. Meng, W.-Y. Xu, H.-L. Wang, X.-Y. Zhuang, W.-C. Xie, and T. Rabczuk, “DigiSim - an open source software package for heterogeneous material modeling based on digital image processing,” Advances in Engineering Software, vol. 148, Article ID 102836, 2020.
View at: Publisher Site | Google Scholar
C. Zhu, M.-c. He, M. Karakus, X.-h. Zhang, and Z. Guo, “The collision experiment between rolling stones of different shapes and protective cushion in open-pit mines,” Journal of Mountain Science, vol. 18, no. 5, pp. 1391–1403, 2021.
View at: Publisher Site | Google Scholar
Q. Yin, J. Y. Wu, C. Zhu et al., “Shear mechanical responses of sandstone exposed to high temperature under constant normal stiffness boundary conditions,” Geomechanics and Geophysics for Geo-Energy and Geo-Resources, vol. 7, no. 2, 2021.
View at: Publisher Site | Google Scholar
W. Wang, G. Chen, Z. Han, S. Zhou, H. Zhang, and P. Jing, “3D numerical simulation of debris-flow motion using SPH method incorporating non-Newtonian fluid behavior,” Natural Hazards, vol. 81, no. 3, pp. 1981–1998, 2016.
View at: Publisher Site | Google Scholar
Q. Wang, Z. Jiang, B. Jiang, H. Gao, Y. Huang, and P. Zhang, “Research on an automatic roadway formation method in deep mining areas by roof cutting with high-strength bolt-grouting,” International Journal of Rock Mechanics and Mining Sciences, vol. 128, Article ID 104264, 2020.
View at: Publisher Site | Google Scholar
Y. Wang, W. K. Feng, R. L. Hu, and C. H. Li, “Fracture evolution and energy characteristics during marble failure under triaxial fatigue cyclic and confining pressure unloading (FC-CPU) conditions,” Rock Mechanics and Rock Engineering, vol. 54, no. 2, pp. 799–818, 2021.
View at: Publisher Site | Google Scholar
Q. Wang, Y. Wang, M. He et al., “Experimental research and application of automatically formed roadway without advance tunneling,” Tunnelling and Underground Space Technology, vol. 114, no. 3, Article ID 103999, 2021.
View at: Publisher Site | Google Scholar
F. Wu, R. Gao, J. Liu, and C. Li, “New fractional variable-order creep model with short memory,” Applied Mathematics and Computation, vol. 380, Article ID 125278, 2020.
View at: Publisher Site | Google Scholar
F. Wu, H. Zhang, Q. Zou, C. Li, J. Chen, and R. Gao, “Viscoelastic-plastic damage creep model for salt rock based on fractional derivative theory,” Mechanics of Materials, vol. 150, Article ID 103600, 2020.
View at: Publisher Site | Google Scholar
A. Li, F. Dai, Y. Liu, H. Du, and R. Jiang, “Dynamic stability evaluation of underground cavern sidewalls against flexural toppling considering excavation-induced damage,” Tunnelling and Underground Space Technology, vol. 112, Article ID 103903, 2021.
View at: Publisher Site | Google Scholar
J. S. H. Kwan, “Supplementary technical guidance on design of rigid debris-resisting barriers,” Geotechnical Engineering Office, Civil Engineering Dept, Chennai, India, 2012.
View at: Google Scholar

Copyright

Copyright © 2021 Bailong Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies