A method, system and device for calculating the entire process of slope instability based on particle method

Abstract

The present invention discloses a method, system and device for calculating the entire process of slope instability based on a particle method. The method of the present invention is based on the smoothed particle hydrodynamics (SPH) particle calculation method, combined with the Riemann large deformation theory, and realizes regional-scale landslide simulation by programming in C++ / CUDA. The method can truly reproduce the entire process of sliding surface generation, penetration, startup and accumulation, and can analyze the influence mechanism of rock and soil properties such as internal friction angle, cohesion, particle size, elastic modulus, and Poisson's ratio on the landslide process. At the same time, the method improves the computational efficiency by nearly 100 times compared with the traditional CPU programming efficiency through the parallel heterogeneous CPU / GPU programming architecture, and realizes the prediction of watershed-scale landslide flow velocity and accumulation range. It has strong adaptability and overcomes the defects of discrete element, finite element and other methods that are limited by simulation scale and computational efficiency, providing a new idea for the prediction of large-scale landslide disasters.

Description

Technical Field

[0001] The present invention belongs to the technical field of landslide prediction, and in particular relates to a method, system and equipment for calculating the entire process of slope instability based on a particle method. Background Art

[0002] Once a landslide occurs, it will have a massive impact and bury the landslide, posing a serious threat to residents' lives and property. Numerical simulation is a commonly used engineering safety analysis method. Traditional finite element methods are often used to simulate the startup phase of sliding surface formation. However, due to mesh distortion limitations, they struggle to simulate the large deformation phase and are unable to predict the landslide operation and accumulation phase, which is of greatest concern to risk management personnel. Discrete element methods can overcome the limitations of finite element methods and enable large deformation simulations, but their core contact calculation algorithm leads to low computational efficiency. Currently, their application is primarily focused on laboratory-scale mechanism research, making it difficult to conduct predictive analysis on large-scale slopes.

[0003] Therefore, realizing the simulation of rock and soil from small deformation to large deformation and solving the calculation problems caused by real terrain are the current problems that need to be solved urgently for on-site landslide prediction and analysis. Summary of the Invention

[0004] To address the problem that existing technologies have difficulty simulating rock and soil deformation from small to large deformation and are limited by the computational explosion caused by real terrain, the present invention provides a method, system and equipment for calculating the entire process of slope instability based on the particle method. Through multi-threaded parallel computing, based on the consideration of the constitutive model of the rock and soil, efficient numerical calculation of the entire process of landslide initiation, migration and accumulation is achieved.

[0005] According to one aspect of the present invention, a method for calculating the entire process of slope instability based on a particle physics method is provided, comprising:

[0006] Obtain the STL real terrain data of the slope site and import it into the smooth particle hydrodynamics calculation framework;

[0007] In the smooth particle fluid dynamics calculation framework, particles are generated and the calculation domain is divided and the particles are numbered. At the same time, the rock and soil parameter attributes are assigned to each particle.

[0008] With each particle corresponding to one thread, all particles are subjected to multi-threaded parallel calculations, where each particle performs the following cyclic calculations: the stress rate and adaptive time step of each rock and soil particle are calculated, the rock and soil particle stress is updated and stored; when the particle stress does not meet the ideal elastic-plastic criterion, the large deformation of the rock and soil is calculated, the particle velocity and density are updated and stored; time integration is performed according to the Verlet algorithm, the particle position is updated and stored; when the set single-step output time is reached, the stored particle stress, particle velocity, particle density and particle position are output;

[0009] Repeat the loop calculation process for each particle until the current physical simulation time reaches the preset simulation time.

[0010] As a further technical solution, the STL real terrain data of the slope site is obtained and imported into the smooth particle hydrodynamics calculation framework, including:

[0011] Obtain site point cloud data of the slope;

[0012] An initial slope terrain in STL format is generated according to the site point cloud data and is imported into a smoothed particle fluid dynamics calculation framework.

[0013] As a further technical solution, within the smooth particle hydrodynamics calculation framework, particles are generated, the computational domain is divided, and the particles are numbered. At the same time, geotechnical parameter attributes are assigned to each particle, including:

[0014] Generate SPH particles according to the set initial particle spacing and search radius, and obtain the terrain represented by SPH particles with fixed spacing;

[0015] Divide the entire computational domain into grids and number the particles and grids respectively;

[0016] Determine the site's geotechnical parameter properties and assign them to the SPH particles, and also determine the artificial sound velocity of the SPH particles.

[0017] As a further technical solution, the determined SPH particle artificial sound speed is greater than 10 times the maximum particle speed, or meets Where E represents the elastic modulus and ρ represents the particle density.

[0018] As a further technical solution, the stress rate and adaptive time step of each rock and soil particle are calculated, and the stress of the rock and soil particles is updated and stored, including:

[0019] Perform adjacent particle search, calculate the particles in the grid whose distance to the calculated particle is less than the set value, and determine the particle pairs in the search domain;

[0020] According to the ideal elastic-plastic constitutive model and combined with Riemann's large deformation theory, the stress rate of rock and soil particles is solved;

[0021] Calculate the time step;

[0022] According to the solved stress rate and time step, the stress information of the rock and soil particles is updated and stored.

[0023] As a further technical solution, when the particle stress does not meet the ideal elastic-plastic criterion, the large deformation of the rock and soil mass is calculated, and the particle velocity and density are updated and stored, including:

[0024] Solve the artificial viscosity and artificial stress of each rock and soil particle;

[0025] Solving the mass conservation and momentum conservation equations under SPH discretization based on the artificial viscosity, artificial stress, and updated stress information to obtain the density change rate and particle change rate;

[0026] Update particle velocity and density information and store them.

[0027] As a further technical solution, time integration is performed according to the Verlet algorithm to update and store the particle positions, including:

[0028] After all threads have completed calculations, the time integration is performed using the Verlet algorithm to update and store the particle position information.

[0029] According to one aspect of the present invention, a system for calculating the entire process of slope instability based on a particle method is provided, comprising a first processor and a second processor, wherein:

[0030] The first processor is used to obtain STL real terrain data of the slope site and import it into the smooth particle hydrodynamics calculation framework; under the smooth particle hydrodynamics calculation framework, generate particles and perform calculation domain division and particle numbering, and assign rock and soil parameter attributes to each particle;

[0031] The second processor is configured to receive information transmitted by the first processor and perform multi-threaded parallel calculations on all particles in a manner where each particle corresponds to one thread, wherein each particle performs the following cyclic calculations: calculating the stress rate and adaptive time step of each rock and soil particle, updating and storing the rock and soil particle stress; when the particle stress does not meet the ideal elastic-plastic criterion, calculating the large deformation of the rock and soil, updating and storing the particle velocity and density; performing time integration according to the Verlet algorithm, updating and storing the particle position; and outputting the stored particle stress, particle velocity, particle density, and particle position when the set single-step output time is reached;

[0032] The first processor is further configured to receive data output by the second processor and determine whether the current numerical simulation has actually reached a preset simulation time. If the preset simulation time has been reached, the sequential calculation process of the second processor is stopped.

[0033] According to one aspect of the present invention, a device for calculating the entire process of slope instability based on the particle method is provided, comprising a processor and a memory for storing a computer program that can be run on the processor, wherein the processor is configured to execute the steps of the method for calculating the entire process of slope instability based on the particle method when running the computer program.

[0034] According to one aspect of the present invention, a computer storage medium is provided, wherein a computer program is stored in the computer storage medium. When the computer program is executed by a processor, the steps of the method for calculating the entire process of slope instability based on the particle method are implemented.

[0035] Compared with the existing technology, the present invention provides an efficient numerical calculation method that can take into account the constitutive model of rock and soil and simulate the entire process of landslide initiation, migration and accumulation. Its beneficial effects are:

[0036] (1) Based on the Riemann large deformation theory, the present invention breaks through the limitation of the traditional method of landslide calculation that cannot consider small deformation and large deformation at the same time, and establishes a method that can predict the entire destruction process of slope rock and soil from small deformation of soil to large deformation of landslide.

[0037] (2) The method of the present invention is based on the ideal elastic-plastic rock and soil constitutive model, completes the simulation of the plastic zone generation and penetration process of the soil-rock mixture, and can accurately predict the location of the plastic zone generation and the penetration time of the rock and soil.

[0038] (3) The method of the present invention introduces the ideal elastic-plastic criterion into the smooth fluid particle dynamics, completes the sliding and accumulation process of the soil-rock mixture, and can predict the sliding speed and accumulation range of the landslide body.

[0039] (4) The method of the present invention takes into account the influence of the initial conditions of the rock and soil mass, and parameter sensitivity analysis can be carried out on the site rock and soil mass parameters (internal friction angle, cohesion, density, etc.) as needed.

[0040] (5) The method of the present invention adopts GPU parallel computing programming method, which avoids the long calculation time of traditional methods and can be used for rapid assessment of disaster sites.

[0041] (6) The method of the present invention adopts C++ / CUDA programming, which is easy to encapsulate and can be imported into large commercial software such as Multiphase ABAQUS after simple processing for more complex working condition calculations. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, a brief introduction will be given below to the drawings used in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0043] Figure 1 A schematic flow chart of a method for calculating the entire process of slope instability based on a particle method is provided in an embodiment of the present invention.

[0044] Figure 2Schematic diagram of the entire process of slope instability provided by an embodiment of the present invention.

[0045] Figure 3 A particle-based method-based method for calculating the entire process of slope instability provided by an embodiment of the present invention is used as a schematic diagram of an image of a landslide result in the southwest.

[0046] Figure 4 A schematic diagram of a method for calculating the entire process of slope instability based on a particle physics method implemented using CPU+GPU provided in an embodiment of the present invention.

[0047] Figure 5 A schematic diagram of a system for calculating the entire process of slope instability based on the particle method provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0048] It should be noted that:

[0049] Smoothed particle hydrodynamics (SPH) calculates particles with both positional and state information (stress, strain, velocity, etc.). This method facilitates parallel computation and can embed multiple constitutive models. This provides a promising approach for large-scale computation of rock and soil landslides in realistic terrain. Furthermore, the resulting landslide simulation method should be embeddable in complex water-soil and water-soil-solid numerical calculations, thus offering greater engineering practical value.

[0050] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention. In addition, the technical features in the various embodiments or single embodiments provided by the present invention are arbitrarily combined with each other to form a new technical solution. This combination is not restricted by the sequence of steps and / or structural composition mode, but must be based on the ability of ordinary technicians in this field to implement it. When the combination of technical solutions is contradictory or cannot be implemented, it should be deemed that this combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

[0051] The embodiment of the present invention provides a method for calculating the whole process of slope instability based on particle method, such as Figure 1As shown in the figure, first, the initial slope terrain in STL format is obtained and converted into terrain represented by SPH particles; then, the entire calculation domain is gridded and numbered, and soil parameters are assigned to the SPH particles; then, a multi-threaded parallel method is used to perform a cyclic calculation on each SPH particle, and each cycle process calculates the particle stress, velocity, density and position information, and updates and stores the information; finally, the actual physical simulation time is monitored, and the cyclic calculation process is stopped when the actual physical simulation time reaches the preset simulation time.

[0052] The embodiment of the present invention adopts a storage-in-CPU-computation approach, performing multi-threaded parallel computing in the GPU for each particle in the terrain, thus avoiding the problem of excessively long computation time in traditional methods and enabling rapid assessment of disaster sites. The calculation process for each particle is as follows: First, the stress of rock and soil particles is solved using the ideal elastic-plastic constitutive model combined with the Riemann large deformation theory, establishing a method that can predict the entire destruction process of the slope rock and soil from small deformation of the soil to large deformation of the landslide; then, the ideal elastic-plastic criterion (i.e., the Drucker-Parager criterion) is introduced into the smooth fluid particle dynamics to complete the sliding and accumulation process of the soil-rock mixture, and the sliding speed and accumulation range of the landslide can be predicted; then, based on the ideal elastic-plastic rock and soil constitutive model, the generation and penetration process of the plastic zone of the soil-rock mixture is simulated, and the location and penetration time of the plastic zone of the rock and soil can be accurately predicted.

[0053] The embodiment of the present invention provides an efficient numerical calculation method that can take into account the constitutive model of rock and soil and simulate the entire process of landslide initiation, migration, and accumulation. The method specifically includes the following steps:

[0054] Step 1. Obtain site point cloud data through remote sensing images, drone photogrammetry, laser scanners, etc., generate triangular STL file meshes, and establish the real slope terrain.

[0055] Step 2. Import the STL (Standard Tessellation Language) real terrain data into the particle method calculation framework, select the initial particle spacing Dp (SPH particle radius), determine the search radius (h = 1.4Dp), generate SPH particles in the real terrain, and use SPH particles with fixed spacing to represent the terrain. At the same time, divide the entire calculation domain into a grid with a side length of 2h, and number the particles and grids respectively.

[0056] Step 3. According to the site conditions, determine the internal friction angle (Φ), cohesion (c), elastic modulus (E), Poisson's ratio (e), density (ρ) and other geotechnical parameters and assign them to the SPH particles; determine the artificial sound velocity (Cs) of the SPH particles, write it into the .exe file compiled by Visual Studio by reading the XML file, and store all information in the CPU (central processing unit) memory in Visual Studio using the C++ language.

[0057] It should be noted that when assigning the parameters of rock and soil particles, the goal of fast assignment and change of assignment is achieved by writing XML reading files in Visual studio, and the input and calculation modules are separated to achieve the purpose of easy operation.

[0058] Furthermore, the artificial sound velocity has a regulating effect on the stability and efficiency of the calculation. Once the actual value of the rock and soil (600-700m / s) is taken, the calculation step size is too small, which affects the calculation efficiency. At the same time, the artificial sound velocity cannot be too small, which is easy to generate calculation noise. Therefore, it is generally believed that Cs>10Vmax, that is, the artificial sound velocity only needs to be greater than 10 times the maximum particle velocity, or meet That is, high calculation accuracy can be achieved while ensuring calculation efficiency.

[0059] Step 4. In Visual Studio, use CUDA to access empty GPU (graphics processing unit) memory and copy the CPU information to the GPU using pointers. Perform the computation on the GPU, placing each particle in a separate thread, with each thread performing the same computation steps.

[0060] In this embodiment of the present invention, information in the CPU is copied to the GPU via pointer transfer using the CUDA language in Visual Studio, preparing for rapid calculations on the GPU. This completes the core link of the CPU-GPU computational framework for simulating the entire process of slope instability, improving computational efficiency by hundreds of times.

[0061] Step 5. Calculate the stress (σ) of the soil particles. For a particular particle, search the nine grids (in 2D) or 27 grids (in 3D) adjacent to the particle's center. Calculate the particles within the grids whose distance from the particle is less than 2h. Use the kernel particle sorting method in CUDA to quickly identify particle pairs within the search domain.

[0062] Step 6. Solve the stress rate of rock and soil particles by combining the ideal elastic-plastic constitutive model with the Riemann large deformation theory

[0063]

[0064] In the formula is the stress rate, w is the moment of inertia, is the plastic strain rate operator, G is the shear modulus, K is the bulk modulus, is the strain rate, ψ is the dilatancy angle, J2 is the second stress tensor invariant, s is the shear stress, is the deviatoric stress tensor rate, δ is the Dirac factor, and s is the deviatoric stress tensor. The stress information of the particle is saved in the designated GPU.

[0065] Step 7. Calculate the time step. The stable variable time step is solved by the CFL condition, which can generally be determined by the following formula:

[0066]

[0067] Δt=C CFL min(Δt f ,h / C s )

[0068] Where f a is the maximum stress value, C CFL Generally, 0.05 is taken as the value, and h is the radius of the calculation core. After solving the stress rate and time step, the stress of each rock and soil particle can be updated.

[0069] Step 8. Substitute stress into the mass conservation and momentum conservation equations under SPH discretization, and in order to obtain stable and accurate stress under SPH method calculation, use the following formula:

[0070]

[0071]

[0072] Among them, v ij is the particle velocity difference, I is the unit matrix, g is the gravitational acceleration, ∏ ij For artificial viscosity, S ij is the artificial stress used to stabilize the calculation equation, is the divergence, ρ i is the density of the i-th particle, ρ j is the density of the jth particle, m j is the mass of the jth particle, W ij is the Winrider kernel function, σ i is the stress of the i-th particle, σ j is the stress of the jth particle. Thus, the ideal elastic-plastic constitutive relation of the soil is embedded in the SPH discrete framework, the calculation of the relative force is completed, and the velocity and density information of the particles are saved in the specified GPU.

[0073] When calculating high and steep slopes, we usually use stress correction. The correction method is to allow the particles to reach equilibrium using their own gravity in a short period of time, then set their displacement to zero, return their particle positions to their initial positions, and retain the stress distribution obtained at this time as the initial stress distribution, so as to quickly obtain the true stress distribution.

[0074] When solving the governing equations for conservation of mass and momentum and storing particle velocity (v) information, the embodiment of the present invention calls the syncthreads() function in CUDA to wait for all GPU threads to complete calculations. This ensures that out-of-bounds access does not occur during GPU calculations.

[0075] Specifically, solve the artificial viscosity (Π ij ), artificial stress (S ij ), when solving the artificial viscosity,

[0076]

[0077] Its coefficient α ΙΙ Usually the value is 0.5, β ΙΙ Generally, we do not consider taking 0, c ij is the average artificial sound speed of the particles, is the average density, h ij is the average nuclear radius, r ij is the interparticle distance, φ ij The calculation is as follows:

[0078]

[0079] S ij Artificial stress can be expressed as follows:

[0080]

[0081] in, R i and R j Calculated by the following formula:

[0082]

[0083] ∈ is a constant between [0,1], in this case 0.2, P i If the artificial viscosity coefficient is not 0.5 and the artificial stress coefficient is not 0.2 in the calculation, it will be difficult to achieve the goal of stable calculation speed.

[0084] Step 9. Perform time integration. Using the modified Verlet integration scheme, at each step, the interactions between the particle pairs are calculated, and then the variables of each particle are updated using the following explicit integral,

[0085]

[0086] Among them F a is the mass force, v a is the relative velocity of the particle pair, r a is the particle spacing, σ is stress, ε is strain, S is the stress rate, M is the strain rate, z is the current time, and Z+1 is the next time step. When solving for geotechnical structures, it's not necessary to update the particle density, but it is necessary to update the particle stress and strain. The new particle position information r is saved on the designated GPU, and the CUDA function -syncthreads() is called to wait for all GPU threads to complete the computation.

[0087] Step 10. Output the saved position, velocity, stress and other information as a binary file (VTK), and determine whether the current simulation time has reached the required simulation physical time. If not, repeat steps 5 to 10. If it has, exit the loop and end the calculation.

[0088] See also Figure 2 and Figure 3 ,in, Figure 2 A schematic diagram of the entire process of slope instability is given. Figure 3 A schematic diagram of the result image of a landslide in southwest China is given in which the particle-based method for calculating the entire process of slope instability described in an embodiment of the present invention is applied.

[0089] Table 1 GPU parallel computing

[0090]

[0091] It should be noted that the embodiment of the present invention calls GPU calculations to take advantage of the fact that the GPU has a large number of computing units. By assigning particle information to each minimum computing unit thread, the corresponding thread's position in the GPU is found when calling to perform calculations. Specifically, by taking advantage of the GPU's multi-core and multi-computing channel characteristics, the SPH particles are placed in separate minimum GPU computing units (threads). When calculating the interaction force and self-stress, independent calculations can be performed by finding the position of the corresponding thread in the GPU without the phenomenon of calculation waiting. The calling process is shown in Table 1.

[0092] As a preferred embodiment, Figure 4As shown, an embodiment of the present invention provides a method for calculating the entire process of slope instability based on a particle method, comprising the following steps:

[0093] (1) Get STL terrain data and set D in Visual Studio p = 10m (SPH particle radius), generate SPH particles, a total of 328439 particles, and number the particles sequentially;

[0094] (2) Determine the search radius h = 14 m, input the adaptive step size CFL = 0.05, gravity acceleration g = 9.8 N / kg, and density ρ = 2000 kg / m 3 , artificial sound velocity Cs = 120 m / s, soil internal friction angle Φ = 30°, cohesion 4 kPa, elastic modulus E = 20 MPa, Poisson's ratio e = 0.2;

[0095] (3) In Visual Studio, use the C++ language to select the free memory of the CPU (central processing unit) and store the particle position information (x, y, z) and attribute information (initial velocity, stress, soil elastic modulus, etc.).

[0096] (4) In Visual Studio, use CUDA to call the empty GPU (image processor) memory and copy the information in the CPU to the GPU through pointer passing.

[0097] (5) The calculation part is carried out in the GPU, each particle is placed in a separate thread, and each thread performs the same calculation steps.

[0098] (6) Calculate the stress rate of each rock particle in the thread As well as the adaptive time step (Δt), the stress (σ) of the rock and soil particles is updated and the stress (σ) information is stored.

[0099] (7) Solve the artificial viscosity (Π), artificial stress (S ij ), when solving for artificial viscosity,

[0100]

[0101] Its coefficient α ΙΙ Usually the value is 0.5, β ΙΙ Generally, we do not consider taking 0, c ij is the average artificial sound speed of the particles, is the average density, φ ij The calculation is as follows:

[0102]

[0103] S ij Artificial stress can be expressed as follows:

[0104]

[0105] in, R i and R j Calculated by the following formula:

[0106]

[0107] ∈ is a constant between [0,1], in this case 0.2. i is the hydrostatic pressure.

[0108] (8) Solve the main governing equations of conservation of mass and momentum, store the particle velocity (v) information, and call the -syncthreads() function in CUDA to wait for all GPU threads to complete the calculation.

[0109] (9) Update the particle position information through the Verlet algorithm and store it in the GPU.

[0110] (10) The particle stress information, velocity information, and position information in the GPU are stored in the CPU, and the memory in the GPU is deleted. The information in the CPU is imported and written into a binary file for storage.

[0111] (11) Check whether the current physical simulation time reaches the planned simulation time. If not, repeat steps 5 to 11. If satisfied, exit the calculation.

[0112] Based on the same inventive concept as the above method embodiment, the embodiment of the present invention further provides a slope instability whole process calculation system based on particle method, including a first processor and a second processor, such as Figure 5 As shown,

[0113] The first processor is used to obtain STL real terrain data of the slope site and import it into the smooth particle hydrodynamics calculation framework; under the smooth particle hydrodynamics calculation framework, generate particles and perform calculation domain division and particle numbering, and assign rock and soil parameter attributes to each particle;

[0114] The second processor is configured to receive information transmitted by the first processor and perform multi-threaded parallel calculations on all particles in a manner where each particle corresponds to one thread, wherein each particle performs the following cyclic calculations: calculating the stress rate and adaptive time step of each rock and soil particle, updating and storing the rock and soil particle stress; when the particle stress does not meet the ideal elastic-plastic criterion, calculating the large deformation of the rock and soil, updating and storing the particle velocity and density; performing time integration according to the Verlet algorithm, updating and storing the particle position; and outputting the stored particle stress, particle velocity, particle density, and particle position when the set single-step output time is reached;

[0115] The first processor is further configured to receive data output by the second processor and determine whether the current numerical simulation has actually reached a preset simulation time. If the preset simulation time has been reached, the sequential calculation process of the second processor is stopped.

[0116] In this embodiment of the present invention, the first processor can be implemented using a CPU, and the second processor can be implemented using a GPU. By adopting a CPU-based storage and GPU-based computation approach, this embodiment of the present invention performs multi-threaded parallel computation on the GPU for each particle in the terrain, avoiding the time-consuming computational overhead of traditional methods and enabling rapid disaster site assessments.

[0117] An embodiment of the present invention also provides a device for calculating the entire process of slope instability based on the particle method, comprising a processor and a memory for storing a computer program that can be run on the processor. When the processor is used to run the computer program, it executes the steps of the method for calculating the entire process of slope instability based on the particle method.

[0118] The memory in the embodiment of the present invention is used to store various types of data to support the operation of the device for calculating the entire process of slope instability based on the particle method. Examples of such data include any computer program for operating the device for calculating the entire process of slope instability based on the particle method.

[0119] An embodiment of the present invention further provides a computer storage medium storing a computer program. When the computer program is executed by a processor, the steps of the method for calculating the entire process of slope instability based on the particle method are implemented.

[0120] In summary, the present invention discloses a particle-based method for calculating the entire process of slope instability. This method is based on the smoothed particle hydrodynamics (SPH) particle calculation method, combined with the Riemann large deformation theory, and implements regional-scale landslide simulation through C++ / CUDA programming. This method can realistically reproduce the entire process of sliding surface generation, penetration, startup, and accumulation, and can analyze the influence of rock and soil properties such as internal friction angle, cohesion, particle size, elastic modulus, and Poisson's ratio on the landslide process. At the same time, this method uses a parallel heterogeneous CPU / GPU programming architecture, with computational efficiency nearly 100 times higher than traditional CPU programming efficiency, achieving basin-scale landslide velocity and accumulation range prediction. It has strong adaptability and overcomes the limitations of discrete element and finite element methods, which are limited by simulation scale and computational efficiency, providing new ideas for large-scale landslide disaster prediction.

[0121] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the technical solutions of the embodiments of the present invention.

Claims

1. A method for calculating the entire process of slope instability based on particle physics, characterized in that: include: Obtain the STL real terrain data of the slope site and import it into the smooth particle hydrodynamics calculation framework; In the smooth particle fluid dynamics calculation framework, particles are generated and the calculation domain is divided and the particles are numbered. At the same time, the rock and soil parameter attributes are assigned to each particle. With each particle corresponding to one thread, all particles are subjected to multi-threaded parallel calculations, where each particle performs the following cyclic calculations: the stress rate and adaptive time step of each rock and soil particle are calculated, the rock and soil particle stress is updated and stored; when the particle stress does not meet the ideal elastic-plastic criterion, the large deformation of the rock and soil is calculated, the particle velocity and density are updated and stored; time integration is performed according to the Verlet algorithm, the particle position is updated and stored; when the set single-step output time is reached, the stored particle stress, particle velocity, particle density and particle position are output; Repeat the loop calculation process for each particle until the current physical simulation time reaches the preset simulation time.

2. The method for calculating the entire process of slope instability based on the particle method according to claim 1 is characterized in that: Obtain the STL real terrain data of the slope site and import it into the smooth particle hydrodynamics calculation framework, including: Obtain site point cloud data of the slope; An initial slope terrain in STL format is generated according to the site point cloud data and is imported into a smoothed particle fluid dynamics calculation framework.

3. The method for calculating the entire process of slope instability based on the particle method according to claim 1 is characterized in that: In the smooth particle fluid dynamics calculation framework, particles are generated, the calculation domain is divided, and the particles are numbered. At the same time, the rock and soil parameter attributes are assigned to each particle, including: Generate SPH particles according to the set initial particle spacing and search radius, and obtain the terrain represented by SPH particles with fixed spacing; Divide the entire computational domain into grids and number the particles and grids respectively; Determine the site's geotechnical parameter properties and assign them to the SPH particles, and also determine the artificial sound velocity of the SPH particles.

4. The method for calculating the entire process of slope instability based on the particle method according to claim 3 is characterized in that: The artificial sound velocity of the SPH particles is determined to be greater than 10 times the maximum particle velocity, or Where E represents the elastic modulus and ρ represents the particle density.

5. The method for calculating the entire process of slope instability based on the particle method according to claim 1 is characterized in that: Calculate the stress rate and adaptive time step of each rock particle, update the rock particle stress and store it, including: Perform adjacent particle search, calculate the particles in the grid whose distance to the calculated particle is less than the set value, and determine the particle pairs in the search domain; According to the ideal elastic-plastic constitutive model and combined with Riemann's large deformation theory, the stress rate of rock and soil particles is solved; Calculate the time step; According to the solved stress rate and time step, the stress information of the rock and soil particles is updated and stored.

6. The method for calculating the entire process of slope instability based on the particle method according to claim 5 is characterized in that: When the particle stress does not meet the ideal elastic-plastic criterion, the large deformation of the rock mass is calculated, and the particle velocity and density are updated and stored, including: Solve the artificial viscosity and artificial stress of each rock and soil particle; Solving the mass conservation and momentum conservation equations under SPH discretization based on the artificial viscosity, artificial stress, and updated stress information to obtain the density change rate and particle change rate; Update particle velocity and density information and store them.

7. The method for calculating the entire process of slope instability based on the particle method according to claim 6 is characterized in that: Perform time integration according to the Verlet algorithm, update the particle position and store it, including: After all threads have completed calculations, the time integration is performed using the Verlet algorithm to update and store the particle position information.

8. A slope instability calculation system based on particle method, characterized by: comprising a first processor and a second processor, wherein: The first processor is used to obtain STL real terrain data of the slope site and import it into the smooth particle hydrodynamics calculation framework; under the smooth particle hydrodynamics calculation framework, generate particles and perform calculation domain division and particle numbering, and assign rock and soil parameter attributes to each particle; The second processor is configured to receive information transmitted by the first processor and perform multi-threaded parallel calculations on all particles in a manner where each particle corresponds to one thread, wherein each particle performs the following cyclic calculations: calculating the stress rate and adaptive time step of each rock and soil particle, updating and storing the rock and soil particle stress; when the particle stress does not meet the ideal elastic-plastic criterion, calculating the large deformation of the rock and soil, updating and storing the particle velocity and density; performing time integration according to the Verlet algorithm, updating and storing the particle position; and outputting the stored particle stress, particle velocity, particle density, and particle position when the set single-step output time is reached; The first processor is further configured to receive data output by the second processor and determine whether the current numerical simulation has actually reached a preset simulation time. If the preset simulation time has been reached, the sequential calculation process of the second processor is stopped.

9. A device for calculating the whole process of slope instability based on particle method, characterized in that: The invention comprises a processor and a memory for storing a computer program that can be run on the processor, wherein when the processor is used to run the computer program, the steps of the method for calculating the whole process of slope instability based on the particle method as described in any one of claims 1 to 7 are executed.

10. A computer storage medium storing a computer program, wherein: When the computer program is executed by a processor, the steps of the method for calculating the entire process of slope instability based on the particle method as described in any one of claims 1 to 7 are implemented.

Images