Methods and devices for simulating landslide surge waves

Abstract

This invention provides a method and apparatus for simulating surging waves in deformable landslides. The method includes: Step 1, constructing geological, landslide, and water body models; Step 2, setting model boundary conditions; Step 3, performing deformable landslide calculations, using the DEM method to simulate the landslide, where the entire landslide area is composed of DEM particles, and the forces between particles are calculated based on different contact models, with acceleration and velocity updated; Step 4, performing surging wave calculations, using the SPH method to simulate the water body, where the entire water body area is discretized into several fluid SPH particles, and parameters such as density gradient and acceleration are calculated; Step 5, directly adding the calculated acceleration of the coupling force of the landslide DEM particles and the acceleration of the reaction force acting on the SPH particles to Steps 3 and 4 for overall iterative calculation and displacement update; Step 6, performing displacement correction, directly adding the displacement correction term to the displacement calculation in Step 5 to update the particle positions within one time step.

Description

Technical Field

[0001] This invention belongs to the field of landslide surge disaster simulation technology, specifically relating to a method and apparatus for simulating surges in deformable landslides. Background Technology

[0002] Riverbank landslides caused by earthquakes or water level changes can induce surges during their entry into the water. Landslide-induced surges are often considered a significant secondary disaster, with their impact area potentially exceeding that of the landslide itself. Numerical simulation methods are widely used as a common engineering analysis approach for landslide surge prediction and disaster inversion. Traditional finite element method (FEM) and Eulerian-Lagrange methods are currently the most commonly used mesh-based modeling and calculation methods. However, because both the landslide body and the water body undergo significant deformation during surge evolution, existing numerical methods suffer from accuracy loss, mesh distortion, or jagged "sawtooth" distortion of the free surface of the water body due to the influence of mesh size. Therefore, how to accurately simulate landslide bodies and water bodies is a pressing problem that needs to be solved.

[0003] Current research techniques are mainly based on numerical models that assume traditional inviscid fluids, selectively neglecting the influence of viscous forces within the fluid. However, the viscous properties of fluids have a certain impact on the velocity field during water movement and wave propagation, causing their influence characteristics in the far field and their effect on landslide motion to still need improvement. Furthermore, how to address the non-physical diffusion of particles on the free surface under the meshless method framework in the face of high-speed landslide wave conditions also requires further optimization. Summary of the Invention

[0004] This invention is made to solve the above-mentioned problems, and aims to provide a method and apparatus for simulating deformation landslide surge waves, which can accurately simulate large deformations of soil and water bodies and the propagation of near-field and far-field surge waves without relying on grids and fully considering fluid viscosity parameters and accurately simulating free liquid surfaces.

[0005] To achieve the above objectives, the present invention adopts the following solution:

[0006] <Method>

[0007] like Figure 1 As shown, the present invention provides a method for simulating swell waves in deformable landslides, characterized by comprising the following steps:

[0008] Step 1: Based on the actual landslide and water conditions to be simulated, construct geological, landslide, and water models, and input the mechanical and material parameters of the landslide and water bodies;

[0009] Step 2: Based on the actual landslide and water conditions to be simulated, set the model boundary conditions. Apply gravity to the landslide and water models without additional constraints. Set fixed boundaries for the slope and water body boundaries, i.e., acceleration is 0.

[0010] Step 3: Perform deformation landslide calculation. The deformation landslide simulation adopts the discrete element method (DEM). The entire landslide area is composed of DEM particles. The force between the particles is calculated based on different contact models. The acceleration and velocity of the particles are updated according to Newton's second law.

[0011] Step 4: Perform water surge calculations. The water body simulation uses the Smooth Particle Flow (SPH) method. The entire water body region is discretized into several fluid SPH particles, and parameters such as density gradient and acceleration are calculated. The particle acceleration is calculated using the following formula:

[0012]

[0013]

[0014] In the formula, subscript i represents the target SPH particle; subscript j represents other SPH particles within the influence domain of the target particle; N represents the total number of particles within the influence domain; p represents the pressure; m represents the mass of the SPH particle; and v represents the particle velocity. The gradient operator is the smooth kernel function; γ is a constant, typically taken as 7; ρ represents density; ρ0 represents standard density, taken as the density of water body 1000 kg / m³. 3 ; α and β represent different directions, μ i μ j Describe the viscosity coefficients of particle i and particle j respectively, δ αβ The symbol represents the Kronecker tensor; g represents gravitational acceleration; the arrow on the symbol indicates a vector.

[0015] Step 5: Perform coupling calculation. First, search for DEM particles within the influence range of fluid SPH particles and assign them the same parameters as standard SPH particles (including density, mass, smooth length, etc.) for combined calculation. The acceleration of the coupling force of landslide DEM particles calculated by buoyancy, drag force, and lubrication force and the acceleration of their reaction force acting on SPH particles are directly added to steps 3 and 4 for overall iterative calculation and displacement update.

[0016] Step 6: Perform displacement correction for fluid SPH particles. Due to the lack of a free surface kernel function and uneven particle distribution within the internal influence domain, this invention, based on the aforementioned steps, further introduces a free surface correction method in landslide surge simulation to reduce particle distribution distortion and void error near the free surface when high-speed landslides impact the solid-liquid coupling interface. This includes the following sub-steps:

[0017] Step 6-1: Determine whether the particle to be judged is located on the free liquid surface based on whether there are DEM particles within the smooth length of each direction of the particle to be judged. If DEM particles are present, the particle is determined to be located on the free liquid surface; otherwise, the particle is determined to be located below the free liquid surface.

[0018] Step 6-2: Calculate different correction terms based on whether the particles are located on the free surface. The formulas are shown below:

[0019]

[0020] In the formula, dt represents the time step; v max h represents the maximum flow velocity within the flow field. i W represents the smooth length of particle i; ij Represents the kernel function calculation values ​​for particles i and j;

[0021] Step 6-3: Add the displacement correction term directly to the displacement calculation in step 5 to update the particle position within one time step.

[0022] Preferably, the deformation landslide surge simulation method provided by the present invention may also have the following feature: the pressure p is calculated in step 4 as follows:

[0023]

[0024] Preferably, the deformation landslide surge simulation method provided by the present invention may also have the following feature: in step 4, the density time gradient dρ / dt (density gradient) is calculated using the following formula:

[0025]

[0026] In the formula, ξ is the dissipation coefficient; h represents the smooth length of the SPH particle; and c0 represents the sound velocity of the entire flow field.

[0027] Preferably, the deformation landslide surge simulation method provided by the present invention may also have the following feature: in step 5, the standard SPH particles are SPH particles with a density set to ρ0 and a volume the same as the initially arranged particles.

[0028] <device>

[0029] Furthermore, the present invention also provides a deformation landslide surge simulation device, capable of automatically implementing the above-mentioned <method>, characterized in that it includes:

[0030] The modeling module constructs geological, landslide, and water models based on the actual landslide and water conditions to be simulated, and inputs the mechanical and material parameters of the landslide and water bodies.

[0031] The parameter setting module allows you to set the model boundary conditions based on the actual landslide and water conditions to be simulated. Gravity is applied to the landslide and water models without additional constraints, and fixed boundaries are set between the slope and the water body.

[0032] The deformation landslide simulation module performs deformation landslide calculations using the discrete element method. The entire landslide area is formed by DEM particles, and the forces between the particles are calculated based on different contact models, updating the acceleration and velocity of the particles.

[0033] The water surge simulation module performs water surge calculations. The water simulation employs the smoothed particle flow (SPH) method, discretizing the entire water region into several fluid particles (SPH particles) and calculating parameters such as density gradient and acceleration. Particle acceleration is calculated using the following formula:

[0034]

[0035]

[0036] In the formula, subscript i represents the target SPH particle; subscript j represents other SPH particles within the influence domain of the target particle; N represents the total number of particles within the influence domain; p represents the pressure; m represents the mass of the SPH particle; and v represents the particle velocity. The gradient operator is the smooth kernel function; γ is a constant; ρ represents the density; ρ0 represents the standard density; α and β represent different directions, and μ i μ j Describe the viscosity coefficients of particle i and particle j respectively, δ αβ The symbol represents the Kronecker tensor; g represents gravitational acceleration; the arrow on the symbol indicates a vector.

[0037] The coupled calculation module first searches for DEM particles within the influence domain of fluid SPH particles and assigns them the same parameters as standard SPH particles for merging calculation. The acceleration of the coupled force of the landslide DEM particles and its reaction force, calculated by buoyancy, drag, and lubrication, are added to the deformation landslide simulation module and the water surge simulation module for overall iterative calculation and displacement update.

[0038] The displacement correction module performs displacement correction of fluid SPH particles using the following steps 6-1 to 6-3.

[0039] Step 6-1: Determine whether the particle to be judged is located on the free liquid surface based on whether there are DEM particles within the smooth length of each direction of the particle to be judged. If DEM particles are present, the particle is determined to be located on the free liquid surface; otherwise, the particle is determined to be located below the free liquid surface.

[0040] Step 6-2: Calculate different correction terms based on whether the particles are located on the free surface. The formulas are shown below:

[0041]

[0042] In the formula, dt represents the time step; v max h represents the maximum flow velocity within the flow field. i W represents the smooth length of particle i; ij Represents the kernel function calculation values ​​for particles i and j;

[0043] Step 6-3: Use the displacement correction term as the current displacement in the displacement calculation in the coupled calculation module to update the particle position within one time step;

[0044] The control module communicates with the modeling module, parameter setting module, deformation landslide simulation module, water surge simulation module, coupled calculation module, and displacement correction module to control their operation.

[0045] Preferably, the deformation landslide surge simulation device provided by the present invention may also have the following features: an input display module, which is communicatively connected to the control unit, for allowing the user to input operation commands and display them accordingly.

[0046] Preferably, the deformation landslide surge simulation device provided by the present invention may also have the following features: the input display module can display prompts to allow users to input or import actual landslide and water body information, and based on this information, the modeling module and parameter setting module can construct geological, landslide and water body models and set model boundary conditions.

[0047] Preferably, the deformation landslide surge simulation device provided by the present invention may also have the following features: the input display unit can display the constructed model as a three-dimensional model diagram according to the corresponding operation instructions, and display the set boundary conditions at the corresponding positions of the three-dimensional model diagram. It can also display the model update and simulation status as a static change trend diagram or demonstrate it as a dynamically changing three-dimensional model diagram.

[0048] Preferably, the deformation landslide surge simulation device provided by the present invention may also have the following feature: the pressure p is calculated using the following formula in the water surge simulation module:

[0049]

[0050] Preferably, the deformation landslide surge simulation device provided by the present invention may also have the following feature: in the water surge simulation module, the formula for calculating the density time gradient dρ / dt is as follows:

[0051]

[0052] In the formula, ξ is the dissipation coefficient; h represents the smooth length of the SPH particle; and c0 represents the sound velocity of the entire flow field.

[0053] The role and effect of invention

[0054] The landslide surge simulation method and apparatus provided by this invention, without the need for mesh calculation, considers the true viscosity of the water body and accurately considers the viscous dissipation during surge propagation. It accurately simulates the deformation of the landslide body and the true shape of the water body to simulate the landslide surge, solving the problem of inaccurate prediction of surge propagation and opposite bank rise in landslide surge calculation. It also corrects the free liquid surface, which can effectively reduce the particle distribution distortion and void error near the free liquid surface when high-speed landslides impact the solid-liquid coupling interface. It accurately simulates the soil entering the water and the surge propagation pattern (accurately simulates large deformation of soil and water body and near-field and far-field surge propagation). The constructed landslide surge simulation model can be used in numerical calculations for the prediction and prevention of complex geological secondary disasters and has good engineering application value. Attached Figure Description

[0055] Figure 1 The flowchart is a process for simulating deformable landslide surge waves based on a meshless algorithm, as per the present invention.

[0056] Figure 2 This is a schematic diagram of particle arrangement and water pressure cloud map of the initial conditions of slow (a) and fast (b) landslides according to embodiments of the present invention;

[0057] Figure 3 The diagram shows the flow regime and velocity field results of a slow landslide surge in an embodiment of the present invention, where (a) corresponds to 0.23s, (b) corresponds to 0.41s, and (c) corresponds to 0.52s.

[0058] Figure 4 This is a graph showing the variation of surge height along the path of a slow landslide, as described in an embodiment of the present invention.

[0059] Figure 5 The figures show the comparison results of applying the free surface correction method of the present invention to a case study of rapid landslides, wherein (a) is the result without using the free surface correction method, and (b) is the result using the free surface correction method of the present invention.

[0060] Figure 6 This is a comparison chart of the surge height changes at specific measuring points in a rapid landslide according to an embodiment of the present invention. Detailed Implementation

[0061] The following describes in detail, with reference to the accompanying drawings, the specific implementation scheme of the deformation landslide surge simulation method and apparatus of the present invention.

[0062] <Example>

[0063] like Figure 1 As shown, the deformation landslide surge simulation method used in this embodiment includes the following steps (step I corresponds to steps 1 and 2 in the "Summary of the Invention" section, steps II to VI correspond to steps 3 to 5 in the "Summary of the Invention" section, step VII corresponds to step 6 in the "Summary of the Invention" section, and the content already described in detail in the "Summary of the Invention" section will not be repeated):

[0064] Step I, as Figure 2 As shown, a model is constructed based on the actual conditions of slow and fast landslides to be simulated, the model parameters are matched, and the water body range and landslide body range are determined.

[0065] Step II, Pre-processing: Using Discrete Element Method (DEM) software, fill the isosceles trapezoidal region with DEM particles and export the coordinates of the particle centers to generate a DEM.dat file. Using Hypermesh mesh generation software, after determining the target mesh size, divide the water body into regions and export the sub-mesh node coordinates to generate an SPH.dat file. Use the program's input file to read two coordinate data files, number each point sequentially, and store them. Add other mechanical and geometric parameters. Add retrieval attribute parameter 'a' to the DEM particles and retrieval attribute parameter 'b' to the SPH particles, i.e., type[1~n]. DEM ] = a, type[(n DEM +1)~(n DEM +n SPH )]=b, and the input parameters for the particles are shown in Tables 1 and 2 below.

[0066] Table 1. Particle property parameters for slow landslide examples

[0067]

[0068] Table 2 Particle property parameters for rapid landslide examples

[0069]

[0070]

[0071] In step III, the number of the two types of boundary particles is less than that of the computational domain particles. They are arranged using loop code, and the addition of mechanical and geometric parameters is the same as in step II.

[0072] Step IV: Based on the particle attribute pointer type(), the particles are distinguished. DEM particles enter the DEM calculation module, and SPH particles undergo SPH iterative calculation.

[0073] Step V: The optional contact models for DEM particles are the linear contact model and the simplified Hertz-Mindlin model. The particle updates its position and velocity according to Newton's second law, and updates the contact force according to the force and displacement laws. The contact force is then converted into a force and additional torque acting on the particle's center of mass. The formulas for calculating the particle's linear acceleration and angular acceleration are as follows:

[0074]

[0075]

[0076] In the formula, m is the mass of the DEM particle, v is the velocity, and F is the velocity. c The resultant force of the contact forces, g is the acceleration due to gravity, and F is the resultant force of the contact forces. df→s F bf→s M represents the coupling force acting on the particle. c With M r Let represent the sliding friction torque and rolling friction torque of the particles, respectively, and α be the global damping coefficient of the particle system, which is taken as 0 in this simulation.

[0077] Step VI, the governing equations for SPH particle motion are the SPH discretized form of the Navier-Stokes equations:

[0078]

[0079] In the formula, j represents a particle in the influence domain of particle i, N represents the total number of particles in the influence domain, ρ is the particle density, v is the particle velocity, p is the particle pressure (due to the non-negativity of pressure, if the pressure value calculated according to the equation of state is less than 0, it should be taken as 0), c0 is taken as 10 times the highest flow velocity in the flow field (8.5 in the slow landslide example and 1300 in the fast landslide example), and ρ0 is the default fluid density (1000 kg / m³ in water). 3 ), where γ is usually set to a value of 7.

[0080] Step VII: Perform displacement correction for fluid SPH particles.

[0081] In the above calculation process, the iteration of the motion parameters of DEM particles and the iteration of the motion parameters and density of SPH are carried out in the same iteration program. The iteration process follows the leapfrog method. The force, acceleration and other modules are calculated at the time t and time t+Δt nodes, and the velocity, displacement and other parameters after iteration are calculated at the time nodes t-0.5Δt and t+0.5Δt.

[0082] After the calculation is complete, export the required DEM, SPH coordinates, velocity, pressure, and other parameters to generate a .dat file, and import it into post-processing software (using Tecplot as an example). Figure 3 , Figure 4 It can be seen that the correctness of the model can be verified by analyzing the flow regime, landslide trajectory, and surge height. For example... Figure 5 It is evident that free surface correction has a very significant effect on closing non-physical voids in free surfaces, such as... Figure 6 It can be seen that by comparing the changes in swell height at different locations, the method proposed in this invention has higher accuracy in simulating swell propagation and shore climb.

[0083] In summary, this invention proposes a deformation landslide surge simulation method that does not rely on grid computing and fully considers fluid viscosity parameters and accurately simulates the free surface. It can accurately simulate large deformations of soil and water bodies and the propagation of near-field and far-field surges without relying on grids and fully considering fluid viscosity parameters and accurately simulating the free surface. Furthermore, the parameters selected in this invention can all be obtained through on-site material measurements, empirical formulas, and other methods, making it easy to implement. It can be used as a laboratory-scale auxiliary numerical simulation method and can also be applied to numerical calculations for the prediction and prevention of complex geological secondary disasters.

[0084] Furthermore, in this embodiment, a deformation landslide surge simulation device capable of automatically implementing the above-mentioned method of the present invention is also provided. The device includes a modeling module, a parameter setting module, a deformation landslide simulation module, a water surge simulation module, a coupling calculation module, a displacement correction module, and an input display module.

[0085] The modeling module performs the steps described in step 1 above, constructing geological, landslide, and water models based on the actual landslide and water conditions to be simulated, and inputting the mechanical and material parameters of the landslide and water bodies.

[0086] The parameter setting module performs the steps described in step 2 above, setting the model boundary conditions according to the actual landslide and water conditions to be simulated.

[0087] The deformation landslide simulation module performs the steps described in step 3 above, calculates the deformation landslide, and updates the acceleration and velocity of the particles.

[0088] The water surge simulation module performs the steps described in step 4 above to calculate water surges. The water simulation adopts the smoothed particle flow method, and the entire water body region is discretized into several fluid SPH particles. The density gradient and acceleration of the fluid SPH particles are calculated.

[0089] The coupled calculation module performs the steps described in step 5 above. First, it searches for DEM particles within the influence range of fluid SPH particles and assigns them the same parameters as standard SPH particles for merging calculation. The acceleration of the coupled force of the landslide DEM particles and its reaction force acting on the SPH particles are calculated through buoyancy, drag, and lubrication and added to the deformation landslide simulation module and the water surge simulation module for overall iterative calculation and displacement update.

[0090] The displacement correction module performs the steps described in step 6 above, corrects the displacement of the fluid SPH particles, and uses the correction term as the current displacement in the displacement calculation of the coupled calculation module to update the particle position within one time step.

[0091] The input display module allows users to input operation commands and displays the corresponding information. For example, it can display prompts allowing users to input or import actual landslide and water body information. Based on this information, the modeling and parameter setting modules can construct geological, landslide, and water body models and set model boundary conditions. The input display module can also display the constructed model as a 3D model diagram based on the corresponding operation commands, displaying the set boundary conditions at the appropriate locations on the 3D model diagram. Furthermore, it can display model updates and simulation progress as a static trend chart or as a dynamically changing 3D model diagram.

[0092] The control module is communicatively connected to the modeling module, parameter setting module, deformation landslide simulation module, water surge simulation module, coupled calculation module, displacement correction module, and input display module, controlling their operation.

[0093] The above embodiments are merely illustrative examples of the technical solutions of the present invention. The deformation landslide surge simulation method and apparatus involved in the present invention are not limited to the contents described in the above embodiments, but are defined by the scope of the claims. Any modifications, additions, or equivalent substitutions made by those skilled in the art based on these embodiments are within the scope of protection claimed by the claims of the present invention.

Claims

1. A method of simulating a surge of a transformed landslide, characterized by, The method comprises the following steps: Step 1, according to the actual landslide body and water body to be simulated, a geological, landslide model and water body model are constructed, and mechanical and material parameters of the landslide body and water body are inputted; Step 2, according to the actual landslide body and water body to be simulated, a model boundary condition is set, a gravity action is applied to the landslide and water body model, there is no additional constraint, and a fixed boundary is set for the slope and water body boundary; Step 3, deformation landslide calculation is performed, the deformation landslide simulation adopts a discrete element method, a whole landslide area is formed by DEM particles, forces between the particles are calculated based on different contact models, and accelerations and velocities of the particles are updated; Step 4, water surge calculation is performed, the water body simulation adopts a smooth particle hydrodynamics method, a whole water body area is discretized into a plurality of fluid SPH particles, and particle accelerations are calculated based on the following formula: In the formula, subscript i represents the target SPH particle; subscript j represents other SPH particles in the influence domain of the target particle; N represents the total number of particles in the influence domain; p represents pressure; m represents the mass of the SPH particle; and v represents the particle velocity. is a gradient operator of a smoothing kernel function; γ is a constant; and ρ represents density. p0 denotes the standard density; a, b denote different directions, m i , m j denote the viscosity coefficients of particle i and particle j, respectively, s αβ is the Kronecker tensor; g denotes the acceleration of gravity; the arrow on the symbol denotes a vector; Step 5, coupling calculation is performed, first, DEM particles within an influence domain of the fluid SPH particles are searched, and the same parameters of the standard SPH particles are given to the DEM particles for combined calculation, coupling forces of the landslide DEM particles obtained through calculation of buoyancy, drag force and lubrication force are acted on accelerations of the SPH particles and added to the accelerations in steps 3 and 4, and whole iterative calculation and displacement updating are performed; Step 6, displacement correction of the fluid SPH particles is performed, including the following sub-steps: Step 6-1, whether the target SPH particle is located on a free surface is judged according to whether there is a DEM particle in each direction smooth length of the target SPH particle, if there is a DEM particle, it is judged that the target SPH particle is located on the free surface, otherwise, it is judged that the target SPH particle is located below the free surface; Step 6-2, different correction term calculations are performed according to whether the target SPH particle is located on the free surface, and the formula is as follows: where dt represents the time step; v max represents the maximum flow velocity within the flow field; h i represents the smooth length of particle i; W ij represents the kernel function calculation value of particles i, j; Step 6-3, the displacement correction term is used as the current displacement for displacement calculation in step 5, and the particle position in one time step is updated.

2. The deformation landslide surge simulation method according to claim 1, characterized in that: wherein In step 4, the pressure p calculation formula is as follows:

3. The deformation landslide surge simulation method according to claim 1, characterized in that: wherein, In step 4, the density gradient dρ / dt calculation formula is as follows: In the formula, ξ is a dissipation coefficient; h represents the smooth length of the SPH particle; c0 represents the sound speed of the whole flow field.

4. The deformation landslide surge simulation method according to claim 1, characterized in that: wherein, In step 5, the standard SPH particle is an SPH particle with a density of ρ0 and a volume same as that of the initial arranged particle; the same parameters of the standard SPH particle include density, mass and smooth length.

5. A device for simulating a surge of a landslide, characterized in that, It comprises: A modeling module, according to the actual landslide body and water body to be simulated, a geological, landslide model and water body model are constructed, and mechanical and material parameters of the landslide body and water body are inputted; A parameter setting module, according to the actual landslide body and water body to be simulated, a model boundary condition is set, a gravity action is applied to the landslide and water body model, there is no additional constraint, and a fixed boundary is set for the slope and water body boundary; The deformed landslide simulation module performs deformed landslide calculation, adopts the discrete element method, and forms a whole landslide area by DEM particles. The forces between the particles are calculated based on different contact models to update the acceleration and velocity of the particles. The water surge simulation module performs water surge calculation. The simulation of the water body adopts the smoothed particle hydrodynamics method. The whole water body area is discretized into a plurality of fluid SPH particles. The particle acceleration is calculated based on the following formula: where subscript i represents the target SPH particle; subscript j represents other SPH particles within the influence domain of the target particle; N represents the total number of particles within the influence domain; p represents pressure; m represents the mass of the SPH particle; and v represents the particle velocity; is the gradient operator of the smoothing kernel function; γ is a constant; and ρ represents density. p0 denotes the standard density; a, b denote different directions, m i , m j denote the viscosity coefficients of particle i and particle j, respectively, s αβ is the Kronecker tensor; g denotes the acceleration of gravity; the arrow on the symbol denotes a vector; The coupling calculation module first searches for DEM particles within the influence domain of the fluid SPH particles, and merges and calculates them by assigning the same parameters as the standard SPH particles. The coupling force acceleration of the landslide DEM particles obtained by calculating the buoyancy, drag force, and lubrication force acts on the acceleration of the SPH particles and is added to the deformed landslide simulation module and the water surge simulation module for overall iterative calculation and displacement updating. The displacement correction module performs displacement correction of the fluid SPH particles by adopting the following steps 6-1 to 6-3. Step 6-1: Determine whether the target SPH particle is located on the free surface according to whether there are DEM particles within the smooth length of each direction of the target SPH particle. If there are DEM particles, it is determined that the target SPH particle is located on the free surface. Otherwise, it is determined that the target SPH particle is located below the free surface. Step 6-2: Perform different correction term calculations according to whether the target SPH particle is located on the free surface. The formula is as follows: where dt represents the time step; v max represents the maximum flow velocity within the flow field; h i represents the smooth length of particle i; W ij represents the kernel function calculation value of particles i, j; Step 6-3: Use the displacement correction term as the current displacement for displacement calculation in the coupling calculation module to update the particle position within one time step. The control module is in communication with the modeling module, the parameter setting module, the deformed landslide simulation module, the water surge simulation module, the coupling calculation module, and the displacement correction module, and controls their operation.

6. The apparatus of claim 5, wherein, Further comprising: The input display module is in communication with the control unit and is used for inputting operation instructions and displaying corresponding information.

7. The apparatus of claim 6, wherein, Further comprising: The input display module can display prompt information to let the user input or import actual landslide body and water body information, and use the information to make the modeling module and the parameter setting module construct a geological model, a landslide model, a water body model, and set boundary conditions.

8. The deformed landslide surge simulation device according to claim 6, wherein: wherein, The input display module can display the constructed model in the form of a three-dimensional model diagram and display the set boundary conditions at the corresponding positions of the three-dimensional model diagram according to the corresponding operation instructions. The model updating and simulation conditions can be displayed in the form of a static change trend diagram or demonstrated in the form of a dynamic change three-dimensional model diagram.

9. The deformed landslide surge simulation device according to claim 5, wherein: wherein, In the water surge simulation module, the pressure p calculation formula is as follows:

10. The deformed landslide surge simulation device according to claim 5, wherein: wherein, In the water surge simulation module, the time gradient dρ / dt of the density calculation formula is as follows: In the formula, ξ is a dissipation coefficient; h represents the smoothing length of the SPH particle; c0represents the sound speed of the entire flow field.

Images