The purpose of the study is to implement a new model based on the Navier-Stokes equations for the characterization of landslide sediments interacting with a moving fluid. The model is implemented by combining Hypermesh, the LS-DYNA software and MATLAB. The results show that the main factors affecting the permeability of landslide sediments are the genetic mechanism, the structure and composition of materials, material lithology, and stress. The characteristics and mechanism of permeability changes are determined by adjusting the water levels of fluids. It is found that the permeability of landslide sediments increases at the front and decreases in the middle and at the rear of the considered domain. The magnitude of the internal flow for the whole system decreases first and then increases with a change in the water level during seepage.
A landslide is a common geological disaster in many countries. It usually makes rocks or soil slide down from slopes, which is caused by gravity. Landslides often occur in mountain areas because tunnels are not constructed and the preventive measures are not carried out there. For example, when the water in the reservoir drops, the seepage of the slopes on both sides of the reservoir changes [
However, the research on the characteristics of the seepage of landslide sediments is complex because its characteristics are affected by many factors, such as the geology, the fluid, and the mechanical properties of rocks. Therefore, an experiment is carried out to simulate the research on the pressure of landslide sediment seepage with a computer. Fluid simulation is conducted with a computer, and Fluid-Structure/solid Interaction (FSI) simulation is most complex. The solid of the slope is affected by the vibration and deformation of fluids during landslides, and the flow velocity is changed, which is called fluid-solid coupling. However, most of the studies on fluid-solid coupling are focused on the boundary between the fluid and solid, and their simulation algorithms are too simple [
Since the permeability has a great influence on the stability of slopes, the model of landslide sediments is implemented by using MATLAB to reconstruct its microstructure, and the change of the flow velocity of the landslide sediment seepage is simulated through built-in Navier-Stokes equations. After the change of the flow velocity caused by the concentration is analyzed, the post-processing software is used to process the calculated data so that the characteristics and influencing factors of the permeability of landslide sediments can be explored. Combined with the relationship between the cross-sectional velocity and the conditions of pressure, the seepage of landslide sediments is evaluated, and the rationality of the seepage law obtained by the simulation is verified.
Based on the previous studies and influencing factors of the permeability of landslide sediments, a microstructural model with pores of the landslide sediment is implemented by using three-dimensional modeling, the structure of the landslide sediment under pressure is analyzed, and the flow velocity is calculated according to the Navier-Stokes equation to study the characteristics of landslide sediments during seepage. After the simulation experiment, the fluid model implemented by the physical equation is discussed by analyzing its effectiveness and validity. The study provides a reference for the study of the pressure damage of landslide sediment seepage. The innovation of this study is to implement a model of the landslide sediment based on Navier Stokes equations to study the characteristics and mechanism of the permeability change of landslide sediments when the flow velocity changes.
With rapid development of computer technology, there appear many new research means and methods to study landslide sediments. Zheng et al. [
In short, the common methods used in the study of the existing landslide sediments are deep learning and a three-dimensional numerical simulation. Through the application of these new computer technologies, the hydrodynamic characteristics and the law of landslides are summarized, which is helpful to preventing and predicting landslides. The method used is the numerical simulation method.
The determination of the fluid motion direction is mainly achieved by solving the Navier-Stokes equation to calculate the flow velocity. Since the fluid mechanism is introduced in the 18th century, its theoretical analysis and mathematical calculations are widely used to study the motion of fluids [
Navier-Stokes equations include mass conservation equations and momentum equations. According to the characteristics of different fluids, different hypotheses are made and different Navier-Stokes equations are obtained. The mass conservation equation shows the continuity of the fluid. In rectangular coordinates, the general form of the continuity equation of fluids is:
In
For the fluid with a stable flow rate, the flow rate per unit time is regarded as a constant, which means that the density
The differential expression derived from the momentum conservation law is:
The finite difference method is used in this study.
The implementation of a three-dimensional model of landslide sediments is the first step to simulate the landslide sediment. And three-dimensional modeling methods mainly include the surface rendering method based on contour and the volume rendering method based on the equivalent surface [
Rock mass is composed of rock blocks with different shapes and structural planes. Therefore, the composition of sediments is very complex when landslides occur, which has a great influence on the study of physical and mechanical properties. When the crack size in the structural plane approaches a critical value, the structural parameters of the slope tend to be stable, which is called Representation Elementary Volume (REV). At present, most of the studies are mainly on the analysis of rock structure and mechanical properties to characterize the REV of rocks, and the studies on the the porosity is most. Porosity refers to the proportion of the pore volume to the total volume of porous materials, and it can reflect the structural characteristics of materials. Based on the REV, the structural characteristics of irregular shapes, like cylindrical and conical pores, are selected as the pore to simulate the model.
Joints, also known as fractures, are small fracture structures with no significant displacement on both sides of rock mass after its stress fracture, which have an important influence on the mechanical behavior of rock mass and cause the discontinuity, the heterogeneity, and the anisotropy of rock mass, as well as the change of rock mechanical parameters. Therefore, the equivalent calculation model of rock mass is implemented by combining the random joint three-dimensional network simulation technology with 3 Dimension Distinct Element Code (3DEC).
(1) First, a three-dimensional network model of random joints in the equivalent rock mass is implemented, and a block with a certain volume is established. The joint surface is determined according to the inclination of the joint, the point of the joint surface, and the penetration rate of the joint surface. Finally, the polygon joint surface is generated by cutting the block. The point passed by the joint surface is the coordinate of the center point of the disc joint. Joint surface penetration rate p can be calculated by
In the equation, a is the diameter of the joint disc, l is the length of the intact rock in the extension direction of the joint surface. A three-dimensional network model of random joints in the equivalent rock mass is implemented by using the built-in FISH language according to the random number of joint inclination, dip angle, joint center point coordinates, and penetration rate.
(2) Second, the equivalent calculation model of rock mass is implemented. In the equivalent model, the rock block is processed by using the strain-softening model, which is based on the Mohr-Coulomb criterion. The shear failure is conducted according to the non-associated flow rule, and the tensile failure follows the associated flow rule. In the strain-softening model, two softening parameters need to be defined: shear softening parameter k^{s} and tensile softening parameter k^{l}. k^{s} is the sum of plastic shear strain increment △k^{s}, and k^{l} is the sum of plastic tensile strain increment △k^{l}. △k^{s} and △k^{l} can be calculated by
In the above equations,
In the above equations,
The failure of joints is divided into shear failure, tensile failure, compression-shear failure, and tensile-shear failure, which can be simulated by the Coulomb slip joint model. The calculation principle is as follows: the normal force loss F^{n} and tangential force vector
In the above equations, k_{n} and k_{s} are the normal stiffness and tangential stiffness of joint surface,
When the joint does not slip in the elastic range, the maximum tensile force and shear force of the joint surface are calculated by the following equations:
In the above equations, T is the joint tensile strength, c and
When the joint is broken or sheared and the cohesion and tensile strength decreases to 0, the maximum tensile force of the joint is 0 and the maximum shear force is
Elastic deformation and plastic deformation have an important impact on the structural stability of the uniaxial compression simulation model. Because of this, it is necessary to select an appropriate elastic model and a plastic model. The plastic model usually constructs equations by Mohr Coulomb, which is simple and the parameter setting is clear:
The parameters in the model of landslide sediments are shown in
Particle size/mm | Percentage of the particles group/% | The mass of each particle/kg in model 1 | |
---|---|---|---|
1# | 2# | ||
40 or more | 0.02 | 0.04 | 11.236 |
20 | 0.0458 | 0.626 | 25.2635 |
10 | 0.1404 | 0.1784 | 78.2538 |
5 | 0.1755 | 0.2745 | 98.2876 |
2 | 0.1514 | 0.203 | 85.2652 |
1 | 0.0795 | 0.0579 | 44.9872 |
0.5 | 0.1241 | 0.0698 | 69.3652 |
0.25 | 0.1588 | 0.06 | 89.6332 |
0.075 | 0.0669 | 0.0342 | 37.1221 |
0.075 or less | 0.0376 | 0.0195 | 21.2522 |
The different types of the models of seepage is shown in
It has great significance to improve the quality and efficiency of finite element analysis by using a powerful, convenient and flexible finite element pre-and post-processing tool that can be easily exchanged with many CAD systems and finite element solvers. HyperMesh is a popular pre-processing software at present. It has powerful and high-performance in finite element mesh pre-processing and post-processing. The three-dimensional solid model can also be directly imported into it. HyperMesh can provide a set of advanced, perfect, and easy-to-use toolkits for implementing the models, which greatly improves work efficiency. Here, HyperMesh software is used to simulate the construction of the model of landslide sediments, and then the data are imported into LS-DYNA software for calculation. The simulation mainly studies the permeability distribution and flow velocity changes of landslide sediments by calculating the flow velocity of the permeation, as well as the pressure. The model is set to have three elements: fluid, air, and the structure of landslide sediments. Among them, fluid provides kinetic energy and makes air enter the inner surface of landslide sediments.
Lagrange is selected as the main algorithm of the structure of landslide sediments, ALE is used to overcome the complex deformation of the grid caused by the fluid passing through landslide sediments, and the Navier-Stokes equation is used to calculate the flow velocity of fluids:
Parameter C in the equation is the coefficient generated by the regional coordinate transformation.
Number | Parameters | Parameter settings |
---|---|---|
1 | ||
2 | ||
3 | ||
The fluid-solid coupling analysis needs to turn to ANSYS software. In the analysis process, the fluid and solid parts are separated. The first analysis is the load of the second analysis. If the analysis is completely coupled, the results of the second analysis will affect or become the load of the first analysis to couple the fluid and the solid. The analysis steps of fluid-solid coupling are as follows:
(1) Exporting solid CBD files and fluid CBD files in HyperMesh, and the fluid and solid grids are divided by HyperMesh.
(2) Importing the fluid grid.
(3) Setting analysis type, entering solid grid file, and setting transient analysis and time.
(4) Sorting fluid materials and setting properties.
(5) Setting the characteristics of the default domain in the flow field.
(6) Adding boundary conditions, which correspond to the boundary in the grid.
(7) Initialization.
(8) Solving control settings.
(9) Outputting control settings.
(10) Monitoring variable settings.
(11) Solving.
The algorithm is the “operator” in numerical simulation and model simulation. The setting of the algorithm affects the accuracy and the calculation time of the simulation. The commonly used algorithms based on LS-DYNA software include the ALE algorithm, the Lagrange algorithm, and the Euler algorithm [
According to the model of landslide sediments and the permeability model, the data during seepage are simulated by the software based on the relevant theory of the finite element, and the permeability distribution law and the stress mechanism of landslide sediments are explored.
The increased process of fluids in landslide sediments is simulated, and the flow velocity is recorded when the landslide sediment is saturated.
Number | Low water level (ml/min) | High water level (ml/min) | Flow velocity change at the low water level (%) | Flow velocity change at the high water level (%) |
---|---|---|---|---|
Initial value | 750 | – | – | – |
1 | 420 | 533 | −0.44 | – |
2 | 300 | 313 | −0.60 | −0.39 |
3 | 255 | 307 | −0.66 | −0.40 |
4 | 388 | 198 | −0.48 | −0.60 |
5 | 380 | 310 | −0.49 | −0.39 |
6 | 380 | 340 | −0.49 | −0.34 |
7 | 607 | 359 | −0.19 | −0.30 |
8 | 932 | 390 | 0.24 | −0.24 |
The table shows that the flow velocity at the low water level decreases by about 44% after the initial adjustment. And then the flow velocity continues to decrease after another two adjustments. It is relatively stable after the fourth to sixth adjustment and then continues to increase to 24%. However, the flow velocity at the high water level gradually decreases to 40%, and it begins to increase after the fifth adjustment but does not reach the initial value.
According to the simulation results of the model, the surface displacement and the deformation of the structure of landslide sediments are shown in
When the water level rises first, the structure of landslide sediments does not deform. However, the front part is gradually damaged, the surface becomes uneven, and the deformation of the middle and rear edges occurs as the flow velocity changes. The deformation of the front edge is large under the action of the water flow, and then the deformation becomes smaller and smaller gradually. Since the materials in the model are mainly soft rocks, the deformation speed and the degree of the slope are lower than hard rocks, and the deformation of the slope with little gravel is less than that of the slope with much gravel.
At the initial velocity, the fluid enters the internal structure of landslide sediments from the pores. At this time, there is no obvious difference between pores during seepage. As time passes by, the pores are gradually connected. When the time of the flow reaches 40 s, the flow will reach the rear part of the landslide sediment. The cross-section analysis of landslide sediments shows that the flow velocity in the each section of pores will experience a dynamic process of “rising first and then falling”. V_{Peak} and V_{Stability} represent the peak speed and the steady speed, respectively.
The deformation degrees of the model of landslide sediments have a great influence on the permeability of water injection. According to the calculation of the overall axial stress of the model and the surface displacement, the pressure dropping along the tube is linearly fitted by using the incompressible viscous fluid when it flows steadily in the rough tube. The results are shown in
The figure shows that the permeability of the model is lower than that of the original model. With the increase of the displacement speed, the shape variable of the model is also increasing. This is consistent with the change law of the permeability of porous media geological materials during deformation [
The axial stress
The figure shows that the values of the peak and the stable velocity increase first and then decrease with the increase of deformation degree in the structures of No. 1 and No. 2 pores. The gradual increase of deformation degree makes the peak and stable velocity values increase according to the stress distribution around the circular hole. In the structure of the No. 3 pore, the two parameters also gradually increase. According to the distribution and characteristics of the stress around the circular hole, the axial and radial stress affects the permeability of landslide sediments. In this case, the diameter of pores and the permeability decrease due to the axial deformation and increase due to the radial deformation. The results show that axial stress has a great influence on the change of the peak velocity in the structure of pores, and the stable change of the permeability is affected by the interaction of the axial stress and the radial stress. The deformation degree of the models does not change the flow velocity of pores. The increase of permeability in the front of the landslide sediment is affected by the overall deformation and leakage erosion, and the decline of the permeability in the middle and rear is due to landslides and collapses. After several tests on leakage, a more stable leakage passage of the slope is formed. Some pore leakage passages are compressed by the overall structure and continue to decrease. The other part of the pore structure is gradually penetrated by the fluid, and the leakage passage is continuously widened and the leakage ability is enhanced. And there are also some leakage passages and pore structures that decrease first and then increase. This shows that the structure of the model cannot change the trend of fluid leakage velocity in the pores over time, but it reduces the difference in the flow velocity.
According to the simulation experiment, the main influencing factors of the permeability of landslide sediments are the effect of seepage, the deformation of the internal structure, and the flow velocity. The effect of seepage on permeability is mainly manifested in latent erosion and blockage. Latent erosion is the particles of landslide sediments in the fluid, which increases the flow velocity caused by the increase of voids inside the landslide sediments. Usually, the front edge of landslide sediments is prone to subsurface erosion. Blockage refers to the deposition of fine particles in the structure of landslide sediments when landslide sediments are flowing with the fluid, which mainly occurs in the front edge of the landslide sediment, and the flow velocity is reduced.
According to the above data, it is found that the seepage flow of slopes changes with the change of water levels of fluids. Specifically, it decreases first and then increases, and the permeability changes similarly with the change of the seepage flow. If the seepage flow of landslide sediments is stable, the internal deformation makes the channel destroyed, and the local blockage of the landslide sediments makes the permeability of the whole slope decrease. The seepage channel will be re-formed in the slope after a seepage. In addition, the deformation and surface displacement occur in the interior of the slope, and the deformation cracks appear on the ground, which leads to the loose structure of landslide sediments and the low density and the enhanced permeability of the whole slope. This can be reflected by the increase of the seepage flow after the deformation of the slope. In terms of the influence of the water level on the permeability, the fluid force gradient in the slope is large and the permeability will increase when the water level increases, which accelerates the discharge of fine particles in the structure, resulting in seepage erosion.
The increase in the permeability of the front part of the landslide sediment is caused by the deformation and seepage erosion, while the decrease in the permeability of the middle and rear parts is due to the collapsibility of landslides. After the experiments, the slope is more stable. In the process of the simulation experiment, it is found that the change of the structure of pores is irregular, and the seepage channel of some pores receives the compression of the overall structure and decreases continuously; the other part of the structure of pores is gradually penetrated by the fluid, and the seepage channel in this part is continuously expanding and the seepage capacity is enhanced. However, some seepage channels decrease first and then increase. This indicates that the structure of the model cannot change the trend of the flow velocity of the fluid in the pore over time, but it reduces the difference in the flow velocity. Zhang et al. [
A simple structural model of landslide sediments is implemented based on the three-dimensional modeling technology to explore the influence of structural stress on the seepage of landslide sediments. On this basis, the simulation experiments of the stress and water injection during seepage are carried out. Through the analysis of the characteristics of the overall and local stress of different models, the following conclusions are drawn: (1) based on the basic principle of hydrodynamic equations, the fluid-solid coupling model is designed, and the Navier-Stokes equation is used to simulate the specific situation of fluids; (2) the results show that the internal deformation and expansion of landslide sediments are the main characteristics of the mechanical deformation of landslide sediments; (3) when the stress of internal structure of pores of landslide sediments increases, the shape of the internal structure of pores becomes the affecting factor of stress distribution. Through the analysis of various characteristics of the seepage and surface displacement of landslide sediments, it is found that the seepage of landslide sediments decreases first and then increases with the change of the water level of fluids, and the permeability also changes with the change of the permeability.
There are still some shortcomings of the study, which are summarized as follows: First, the design and solution of the Navier-Stokes equation should be more accurate, and more details should be satisfied to achieve a high-quality simulation effect; Second, only the influence of the porosity of the model during seepage is considered in the numerical simulation, which simplifies the design and shape of the structure of pores in the models, ignoring the influence of gravity. In the follow-up research, more accurate and reasonable algorithms should be introduced to optimize the research details of the characteristics of landslides during seepage, and the structure of pores in the model should be further studied to make the research results more universal and reliable.