Slope simulation method and related equipment

Abstract

The invention provides a side slope simulation method, which comprises the following steps: acquiring side slope simulation parameters, and assigning a preset side slope simulation model according to the side slope simulation parameters to obtain an initial model; and solving the initial model according to a smoothed particle fluid dynamic method to obtain an updated model. And calculating the particle state according to the updating model to obtain slope stress characteristics. And evaluating the slope stress characteristics according to a preset simulation rule to obtain a slope simulation result. The initial model is solved according to a smoothed particle fluid dynamics method to obtain an updated model, so that the use of traditional grid division is avoided, and the method has obvious advantages in processing the problem of large deformation; the particle state is calculated according to the updating model to obtain the slope stress characteristics, the slope stress characteristics are evaluated according to the failure rule and the contact criterion, the slope breaking condition is calculated to obtain the slope simulation result, and therefore simulation of slope rock collapse large deformation is achieved.

Description

Technical Field

[0001] This application relates to the technical field of hydraulic engineering, and particularly relates to a slope simulation method and related equipment. Background Art

[0002] Since the construction of the Yarlung Zangbo River Hydropower Station is imminent, under the change of the reservoir water level, the rock slope near the hydropower station is weathered on the surface and fissures are developed due to surface weathering, resulting in the enrichment of water pressure in the rock fissures, leading to the collapse of the rock slope and the formation of rockfalls, which pose a huge safety hazard to the hydropower project.

[0003] Currently, the finite element method is mostly used to simulate slopes in such scenarios. However, in the actual application process, it is found that the traditional finite element simulation scheme is prone to mesh distortion when solving scenarios with large deformations, resulting in a large difference between the simulation results and the reality and being difficult to directly apply. Summary of the Invention

[0004] In view of this, the purpose of this application is to propose a slope simulation method and related equipment.

[0005] Based on the above purpose, this application provides a slope simulation method, including: obtaining slope simulation parameters, assigning values to a preset slope simulation model according to the slope simulation parameters to obtain an initial model; solving the initial model according to the smoothed particle hydrodynamics method to obtain an updated model; calculating the particle state according to the updated model to obtain slope stress characteristics; evaluating the slope stress characteristics according to a preset simulation rule to obtain a slope simulation result; wherein, the simulation rule includes a failure criterion and a contact criterion.

[0006] In some embodiments, the step of solving the initial model according to the smoothed particle hydrodynamics method to obtain an updated model specifically includes: calculating the relationship between adjacent particles in the initial model to obtain initial particle information; optimizing the initial particle information according to the basis spline function to obtain the updated model.

[0007] In some embodiments, the step of calculating the relationship between adjacent particles in the initial model to obtain initial particle information specifically includes: setting a smoothing kernel radius, pairing adjacent particles in the initial model according to the smoothing kernel radius; calculating the initial particle information according to the linked list search method; the initial particle information is calculated by the following formula: where i represents the serial number of the base particle; j represents the serial number of the particle adjacent to the base particle; N represents the total number of particles in the initial model; f(x i ) represents the field function of the base particle i; Denote the derivative of the field function of the base particle i; ρ represents the density of a single particle; m represents the mass of a single particle; h represents the smooth kernel radius of the base particle i; W(x i -x j ,h) represents the smooth kernel function; Denote the gradient of the smooth kernel function.

[0008] In some embodiments, calculating the slope stress characteristics by calculating the particle state according to the update model specifically includes: calculating the force condition of the particles in the update model to obtain slope mechanical information; solving the slope mechanical information according to the preset water volume information to obtain slope seepage information; calculating the slope stress characteristics according to the preset reservoir water level and the slope seepage information.

[0009] In some embodiments, solving the slope mechanical information according to the preset water volume information to obtain slope seepage information specifically includes: solving the slope mechanical information according to the cubic law to obtain the seepage coefficient; calculating the water level boundary, and calculating the steady-state seepage field according to the seepage coefficient and the water level boundary; calculating the precipitation boundary, and calculating the transient seepage field according to the steady-state seepage field and the precipitation boundary.

[0010] In some embodiments, evaluating the slope stress characteristics according to the preset simulation rules to obtain the slope simulation result specifically includes: evaluating the slope stress characteristics according to the failure criterion to obtain particle failure information; calculating the particle position according to the particle failure information to obtain particle position information; calculating the particle contact condition according to the contact criterion and the particle position information to obtain the slope simulation result.

[0011] In some embodiments, evaluating the slope stress characteristics according to the failure criterion to obtain particle failure information specifically includes: evaluating the damage condition of the particles in the update model according to the Drucker-Prager criterion and the slope stress characteristics to obtain an evaluation result; solving the particle control equation according to the evaluation result to obtain the particle failure information; the particle control equation is calculated by the following formula: Wherein, Denote the particle to be evaluated; W ij represents the smooth kernel function of the particle to be evaluated; ψ is a constant, and ψ = 0 in response to the evaluation result that the particle to be evaluated is damaged.

[0012] The present application also provides a slope simulation device, including: a parameter adjustment module, configured to obtain slope simulation parameters, assign values to a preset slope simulation model according to the slope simulation parameters to obtain an initial model; a model update module, configured to solve the initial model according to the smoothed particle hydrodynamics method to obtain an updated model; a stress calculation module, configured to calculate particle states according to the updated model to obtain slope stress characteristics; and a result output module, configured to evaluate the slope stress characteristics according to preset simulation rules to obtain a slope simulation result.

[0013] The present application also provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, where when the processor executes the program, the method described in any one of the above is implemented.

[0014] The present application also provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the method described in any one of the above.

[0015] As can be seen from the above, the slope simulation method provided by the present application inputs slope simulation parameters into a slope simulation model, initializes the slope simulation model to obtain an initial model, and then solves the initial model according to the smoothed particle hydrodynamics method to obtain an updated model, thereby avoiding the use of traditional mesh generation and having obvious advantages in dealing with large deformation problems; then calculates particle states according to the updated model to obtain slope stress characteristics, and finally evaluates the slope stress characteristics according to failure rules and contact criteria, calculates the situation of slope fragmentation to obtain a slope simulation result, so as to realize the simulation of large deformation of slope rockfall. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solutions in the present application or related technologies, the following will briefly introduce the drawings required for use in the description of the embodiments or related technologies. Obviously, the drawings in the following description are only embodiments of the present application. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0017] Figure 1 is a flowchart of the slope simulation method provided by an embodiment of the present application;

[0018] Figure 2 is a flowchart of the updated model calculation method provided by an embodiment of the present application;

[0019] Figure 3 is a flowchart of the initial particle information calculation method provided by an embodiment of the present application;

[0020] Figure 4 Schematic flow chart of the slope stress characteristic calculation method provided by the embodiment of the present application;

[0021] Figure 5 Schematic flow chart of the seepage information calculation method provided by the embodiment of the present application;

[0022] Figure 6 Schematic flow chart of the slope simulation result generation method provided by the embodiment of the present application;

[0023] Figure 7 Schematic flow chart of the particle failure information calculation method provided by the embodiment of the present application;

[0024] Figure 8 Schematic structural diagram of the slope simulation device provided by the embodiment of the present application;

[0025] Figure 9 Schematic diagram of boundary condition setting provided by a specific embodiment of the present application;

[0026] Figure 10 Schematic diagram of a more specific hardware structure of an electronic device provided by the embodiment of the present application. Detailed implementation manners

[0027] To make the objectives, technical solutions and advantages of the present application clearer and more understandable, the present application will be further described in detail below with reference to specific embodiments and the accompanying drawings.

[0028] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the embodiments of the present application should have the ordinary meanings understood by those of ordinary skill in the art to which the present application belongs. The "first", "second" and similar terms used in the embodiments of the present application do not indicate any order, quantity or importance, but are only used to distinguish different components. The terms such as "including" or "comprising" mean that the elements or objects appearing before this word cover the elements or objects listed after this word and their equivalents, without excluding other elements or objects. The terms such as "connected" or "coupled" are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. The terms such as "upper", "lower", "left" and "right" are only used to represent relative positional relationships, and when the absolute position of the object being described changes, the relative positional relationship may also change accordingly.

[0029] It can be understood that, before using the technical solutions of the various embodiments of the present disclosure, the types, usage scopes, usage scenarios, etc. of the personal information involved will be informed to the user in an appropriate manner and the user's authorization will be obtained.

[0030] For example, when responding to an active request from a user, a prompt message is sent to the user to clearly prompt the user that the operation requested by the user will require obtaining and using the user's personal information. Thus, the user can autonomously choose whether to provide personal information to software or hardware such as an electronic device, an application program, a server, or a storage medium that performs the operations of the present disclosure technical solution based on the prompt message.

[0031] As an optional but non-limiting implementation manner, the manner of sending a prompt message to the user in response to receiving an active request from the user may be, for example, in the form of a pop-up window, and the prompt message may be presented in text in the pop-up window. In addition, the pop-up window may also carry selection controls for the user to select "agree" or "disagree" to provide personal information to the electronic device.

[0032] It can be understood that the above process of notifying and obtaining user authorization is only illustrative and does not limit the implementation manner of the present disclosure. Other manners that meet relevant laws and regulations can also be applied to the implementation manner of the present disclosure.

[0033] This application provides a slope simulation method, as Figure 1 shown, including:

[0034] Step S1, obtain slope simulation parameters, and assign values to a preset slope simulation model according to the slope simulation parameters to obtain an initial model.

[0035] In this embodiment, the slope simulation parameters are the basic mechanical parameters of the slope that need to be simulated and calculated, and may include material parameters (hydraulic parameters and strength) and boundary conditions (rainfall boundary, water level boundary, stress boundary), etc.; the position information of particles, the initial velocity of particles, and the initial stress of particles; the strength parameters mainly include elastic modulus, Poisson's ratio, cohesion, tensile strength, and internal friction angle. Permeability coefficients such as the normal stiffness of cracks, the water yield of rock masses, the bulk modulus of water bodies, the permeability coefficient of the matrix, the initial crack width, and the normal deformation coefficient of cracks, the dynamic viscosity of water bodies, etc.) and boundary conditions (reservoir water level boundary, rainfall boundary, and stress boundary), etc.

[0036] Step S2, solve the initial model according to the smoothed particle hydrodynamics method to obtain an updated model.

[0037] Among them, since the displacement of the fragmented rock mass generated when the slope fails is relatively large, it is easy to generate mesh distortion when using the finite element method for solution. Therefore, in this application, the smoothed particle hydrodynamics method is used to solve the initial model, thereby avoiding the calculation error caused by the finite element mesh deformation.

[0038] Step S3, calculate the particle state according to the updated model to obtain the slope stress characteristics.

[0039] Among them, by calculating the stress characteristics of the slope, the force conditions of each particle in the updated model are determined, and whether the particle is damaged is judged according to the force conditions, so as to obtain the simulation results after the slope failure.

[0040] Step S4, evaluate the slope stress characteristics according to the preset simulation rules to obtain the slope simulation results. Among them, the simulation rules include failure criteria and contact criteria.

[0041] By setting the failure criteria and contact criteria, it is judged whether the slope fails, and the collapse of the rock mass fragments after failure is judged according to the contact criteria.

[0042] In some embodiments, as Figure 2 shown, step S2 specifically includes:

[0043] Step S21, calculate the relationship between adjacent particles in the initial model to obtain the initial particle information.

[0044] Step S22, optimize the initial particle information according to the basis spline function to obtain the updated model.

[0045] In this embodiment, by calculating the relationship between the base particles and adjacent particles in the initial model, the mechanical interaction relationship between each particle and other particles is obtained, and the action relationship of the base particles on adjacent particles is calculated according to the basis spline function, so as to obtain the updated model.

[0046] In some embodiments, as Figure 3 shown, step S21 specifically includes:

[0047] Step S201, set the smoothing kernel radius, and pair the adjacent particles in the initial model according to the smoothing kernel radius.

[0048] Step S202, calculate the initial particle information according to the linked list search method. The initial particle information is calculated by the following formula:

[0049]

[0050] Among them, i represents the serial number of the base particle. j represents the serial number of the particle adjacent to the base particle. N represents the total number of particles in the initial model. f(x i ) represents the field function of the base particle i. represents the derivative of the field function of the base particle i. ρ represents the density of a single particle. m represents the mass of a single particle. h represents the smoothing kernel radius of the base particle i. W(x i -x j ,h) represents the smoothing kernel function. represents the gradient of the smoothing kernel function.

[0051] In this embodiment, the calculation process of the initial particle information can also be optimized by the B-spline function (basis spline function) to eliminate the influence of particle instability during compression. The B-spline function is expressed by the following formula:

[0052]

[0053] where λ is a constant, and its value is different for one-dimensional, two-dimensional, and three-dimensional space problems [18, 19]. For example, for the B-spline function, 1 / h, 15 / (7πh^2), and 3 / (2πh^2) are taken respectively; p = r_ij / h represents the normalized distance, and r_ij is the distance between particles i and j.

[0054] In some embodiments, as Figure 4 shown, step S3 specifically includes:

[0055] Step S31, calculating the force on the particles in the updated model to obtain the slope mechanical information.

[0056] In this embodiment, the slope mechanical information may include: density change information, artificial viscosity information, internal stress change information, and external stress change information, which are expressed by the following formula:

[0057]

[0058] where m represents the mass of the particle; ρ represents the density of the particle; t represents the calculation time; v represents the velocity of the particle; σ represents the total stress tensor; α, β represent the coordinate directions; δ represents the Dirac function, П ij represents the artificial viscosity, which is used to reduce non-physical oscillations; T ij represents the artificial damping, which is used to reduce tensile instability; f i α represents the additional force on the base particle, such as gravity.

[0059] In the above formula, the total stress tensor σ αβ is expressed by the following formula:

[0060] σ αβ =-Pδ αβ +τ αβ ;

[0061] where P represents the isotropic pressure; τ represents the anisotropic shear stress; δ αβ is the Kronecker symbol.

[0062] In the above formula, the pressure-volume-energy state equation of the isotropic pressure P is:

[0063]

[0064] where η = ρ / ρ0 - 1, ρ0 represents the initial density of the particles, P H represents the Hugoniot curve function, and Γ is the Gruneisen parameter. In solid mechanics, for anisotropic shear stress τ, the stress rate and the strain rate are proportional to each other, and the proportionality coefficient is the shear modulus G:

[0065]

[0066] where ε γγ and ε αβ are the strain rate tensors, defined as

[0067]

[0068] To make the information of the material consistent with the strain, the Jaumann rate is introduced, and the shear stress rate can be expressed as:

[0069]

[0070] where R αγ and R βγ are the torsion rate tensors, defined as

[0071]

[0072] Step S32, solve the slope mechanical information according to the preset water volume information to obtain the slope seepage information.

[0073] Among them, since the slope needs to be in direct contact with the water body, when calculating the stress characteristics of the slope, the influence of water seepage on the slope also needs to be considered.

[0074] Step S33, calculate the slope stress characteristics according to the preset reservoir water level and the slope seepage information.

[0075] Among them, since the slope directly bears the pressure of the water body in the reservoir, when calculating the stress characteristics of the slope, the pressure of the reservoir water body on the slope also needs to be calculated according to the reservoir water level. By superimposing the seepage forces generated by the reservoir water level and rainfall onto the Cauchy stress of the total stress tensor, the full coupling characteristics of rainfall-seepage-stress are formed, and the formula is as follows:

[0076]

[0077] In some embodiments, as Figure 5 shown, step S32 specifically includes:

[0078] Step S301, solve the slope mechanical information according to the cubic law to obtain the seepage coefficient.

[0079] In fractured rock masses, deformation fractures are the main channels for groundwater flow. Therefore, when the stress conditions of the fracture system change, the permeability of the rock will change significantly. There is a strong correlation between fracture permeability and fracture aperture, and it is controlled by the effective normal stress. The water flow in rock fractures is simulated using the cubic law, considering the influence of shear stress on the fracture surface, and thus updating the entire seepage field as follows:

[0080]

[0081] The fluid flow in the fracture is an implicit equation. By solving the Laplace operator, the seepage velocity of the fluid and the change in water head can be obtained, and thus the change in the seepage field can be obtained, based on the seepage velocity considering SPH discretization and the seepage water head h i The updated equation is as follows:

[0082]

[0083] k si refers to the permeability coefficient, and μ is the dynamic viscosity. Where σ1 is the maximum principal stress, σ3 is the minimum principal stress, and b0 is the initial crack width. K n is the stiffness coefficient of the fracture, where v b is the Poisson's ratio, β w is the normal deformation coefficient, S s is the specific yield of the rock. P w = i gh i represents the water pressure.

[0084] Based on this, we can calculate the permeability coefficient here according to the maximum principal stress and minimum principal stress values at this point. The maximum principal stress σ1 and the minimum principal stress σ3 can be obtained through the total stress tensor as follows:

[0085]

[0086] Among them, σ xx and σ yy are obtained by taking different x and y superscripts of σ in the total stress tensor αβ as shown.

[0087] Step S302, calculate the water level boundary, and calculate the steady-state seepage field according to the seepage coefficient and the water level boundary.

[0088] On the premise of considering the gravity effect, combine the seepage velocity calculation formula with the pressure head difference Δh ij =h i -h j -y i -y j and rewrite it as:

[0089]

[0090] Meanwhile, set the initial water level boundary according to the formula. According to Dupuit's formula (Harr 1962), the free surface is given by the following formula:

[0091]

[0092] Based on this surface, the initial water level line and the steady-state seepage field can be quickly generated.

[0093] Step S303: Calculate the precipitation boundary, and calculate the transient seepage field based on the steady-state seepage field and the precipitation boundary.

[0094] Among them, in the rainfall-type slope, since the influence of the rainfall boundary on the seepage field needs to be considered, the seepage velocity calculation formula can be further rewritten as:

[0095]

[0096] Among them, Q rainfull is the rainfall infiltration boundary, which can be transformed into Q rainfull = q rainfull where S is the infiltration area and q rainfull is the infiltration intensity.

[0097] In some embodiments, as Figure 6 shown, step S4 specifically includes:

[0098] Step S41: Evaluate the slope stress characteristics according to the failure criterion to obtain particle failure information.

[0099] Step S42: Calculate the particle positions according to the particle failure information to obtain particle position information.

[0100] Step S43: Calculate the particle contact situation according to the contact criterion and the particle position information to obtain the slope simulation result.

[0101] In this embodiment, judge whether the slope stress characteristics meet the failure conditions of the slope material according to the failure criterion, and then calculate the movement trajectories of the damaged particles according to the contact criterion, so as to obtain the final simulation result.

[0102] In some embodiments, as Figure 7 shown, step S41 specifically includes:

[0103] Step S401: Evaluate the damage condition of the particles in the updated model according to the Drucker-Prager criterion and the slope stress characteristics to obtain an evaluation result.

[0104] Step S402: Solve the particle control equation according to the evaluation result to obtain particle damage information.

[0105] The particle control equation is calculated by the following formula:

[0106]

[0107] where, represents the particle to be evaluated. W ij represents the smoothing kernel function of the particle to be evaluated. ψ is a constant, and in response to the evaluation result that the particle to be evaluated is damaged, ψ = 0.

[0108] In this embodiment, usually, rocks fail in a brittle manner. All particles in these materials have the same properties as the base material; therefore, the initiation and propagation of particle damage are determined by the failure criterion of the rock material. The Drucker-Prager yield criterion with tension cutoff is widely used in rock fracture mechanics and has achieved good results. By using the Drucker-Prager yield criterion with tension cutoff, it is possible to determine whether a particle fails in tension or shear. The Drucker-Prager yield criterion is selected to determine the initiation and propagation of shear cracks, and its form is:

[0109]

[0110] I1 = σ xx + σ yy + σ zz ;

[0111] where, σ xx , σ yy and σ zz mean that the superscripts α and β of σ are taken as x, y, z simultaneously.

[0112]

[0113] where, I1 represents the first invariant of the stress tensor; J2 represents the second invariant of the stress tensor; α φ , k c are Drucker-Prager constants.

[0114] Considering the judgment of tensile cracks, the Drucker-Prager criterion is improved as follows:

[0115]

[0116] where R tis the tensile strength. Therefore, the Drucker-Prager criterion with tensile truncation can determine the cracking behavior caused by tension or shear. σ1 is the maximum principal stress, and its calculation method is Once the rock particles reach the failure state, the particles are no longer controlled by the kernel function and can be regarded as falling off from the matrix. And the material properties are consistent with the matrix, thus providing a path for crack propagation. According to the improved Drucker-Prager criterion, the fracture state of the solid and the attributes of the crack, whether it is tensile or shear failure, can be obtained.

[0117] After tensile or shear failure, the damage degree is calculated based on the number of damaged particles occupying the total number of particles within the particle search domain to obtain the damage degree at that point. The damage degree is calculated as follows:

[0118] The damage degree of this particle is counted as the ratio of all damaged particles within twice the kernel radius to the total number of particles within twice the kernel radius of this particle as follows:

[0119]

[0120] where ∑N damge represents the number of all damaged particles within twice the kernel radius, and ∑N total represents the total number of particles within twice the kernel radius.

[0121] After the SPH particles are damaged, a particle channel appears here. Therefore, the SPH governing equations considering damaged and undamaged particles can be rewritten as follows. Here we adopt two sets of smoothing kernel functions. For the stress part, the one used is By introducing the damage factor and then deriving the governing equations for damaged particles and intact particles. By introducing a parameter ψ to improve the derivative of the kernel function, when the local particles satisfy the above failure criterion, ψ = 0. When the local particles do not satisfy the above failure criterion, ψ = 1.

[0122]

[0123] Then the improved governing equations for intact particles and damaged particles are as follows:

[0124]

[0125] In the formula: U represents undamaged particles, represents damaged particles. In the formula, W ij (i.e., W(x i -x j ,h)) and (i.e., ) are the smoothing kernel function of IP and its gradient respectively.

[0126] The damaged particles are converted into discrete particles (DPs), with each damaged particle corresponding to a discrete particle (DPs). Contact pairs are formed between DP particles and DP particles, and between IP particles and DP particles. Therefore, when the distance between the damaged particles and the intact particles reaches a certain condition, new contacts are generated between them. The condition is as follows:

[0127]

[0128] where r i and r d are the smoothing kernel radii of the intact particle and the damaged particle, and l id is the distance between the two particles. When I n > 0 and reaches the set value, an interaction contact force is generated between the two types of particles.

[0129] Taking the intact particle as an example, the resultant force F i acting on it is composed of the contact force and other external forces . Among them, the contact force can be calculated according to the geometric shape and material properties of the particles, and can be decomposed into the normal contact force in the direction of the particle connection line and the perpendicular shear contact force . The shear force is related to the shear increment ΔI s between the particles, and its expression is:

[0130]

[0131] where n is the unit outward normal vector, Δv i is the velocity change of the intact particle after contact, and s n and s s are the normal and shear stiffnesses respectively, which can be calculated by the following formulas:

[0132]

[0133] where E is the shear modulus and κ is the ratio of shear and normal stiffnesses. According to the Coulomb friction hypothesis, the shear forces before and after sliding can be unified into the following expression:

[0134]

[0135] Derive the control equations of SPH considering contact, and finally obtain the fracture seepage and damage process of the rock slope under rainfall and reservoir water level conditions. Since it is an instantaneous transformation, the fracture surface is transformed into a contact surface.

[0136]

[0137] Particle contact pair search calculates the contact force between discontinuous particles (DPs). The calculation method follows Newton's second law, and the final control equation is:

[0138]

[0139] In some optional embodiments, cracks are formed after particle destruction. Based on the cubic law, the permeability coefficient in the cracks changes, thereby affecting the seepage field between the cracks and the matrix again. Since the smooth kernel function of the particles is improved in the rock stress layer, but remains unchanged for the seepage layer, the final seepage-stress-fracture-contact equation based on different kernel functions can be obtained as shown below:

[0140]

[0141] In some optional implementations, the final velocity and position information of the particle is calculated according to the above formula to form a single loop. The leapfrog algorithm is used to perform time step integration, and the position, velocity, temperature, head and density changes of each elementary particle are obtained recursively from the previous step. Specifically:

[0142]

[0143] Where Δt is a single time step, t0 is the initial time, and Δv i is the acceleration, is the rate of change of water head. If the cycle is not completed, return to step S2 and continue the cycle until the predetermined total time step is reached.

[0144] It can be seen from the above embodiments of the present application that the present application avoids the use of traditional grid division by inputting slope simulation parameters into a slope simulation model, initializing the slope simulation model to obtain an initial model, and then solving the initial model according to the smooth particle fluid dynamics method to obtain an updated model, which has obvious advantages in dealing with large deformation problems; then the particle state is calculated according to the updated model to obtain the slope stress characteristics, and finally the slope stress characteristics are evaluated according to the failure rule and contact criterion, and the slope crushing situation is calculated to obtain the slope simulation result, thereby realizing the simulation of large deformation of slope rock collapse.

[0145] As a specific implementation method, based on the plane strain assumption, a seepage stress damage rock slope model without support structure and various boundary conditions such as Figure 9Shown. To fully understand the influence of seepage and stress coupling failure. In the HMD-SPH model, the intermittent joints evenly distributed in the upper part of the slope are about 2.8 m long, 0.2 m thick, with an average joint spacing of about 1.2 m and an inclination angle of 45°. It is assumed that the influence of joints and rainfall intensity on the osmotic pressure distribution is considered. The total number of SPH solid particles is 19037, and the number of virtual particles is about 1100 (the particles of the three-layer displacement boundary are used to prevent particle boundary truncation). The initial particle spacing of the particles is set to 0.2 m. To prevent the false viscosity caused by too large a smoothing kernel radius, for stress calculation, we set the smoothing kernel radius to about 0.9 times the initial particle spacing (hsml = 0.18 m), but for seepage field calculation, we use 1 times the initial particle spacing (hsml = 0.2 m). Considering the stability of the calculation, the time step is set to 1×10-5 step / s. The physical parameters such as strength parameters and hydraulic parameters are shown in Table 1. All physical parameters can be obtained by on-site collection.

[0146] Regarding the boundary setting, it can be seen that there are a total of five boundaries in this simulation (two displacement boundaries and three seepage boundaries): The application of the seepage boundary is the key to ensuring a reasonable seepage simulation, which is divided into the steady seepage field of the reservoir water level and the transient seepage field caused by the rainfall boundary. First, apply the steady seepage field. The left free surface is the reservoir, set the left water head boundary height to 13 m, and the right water head boundary height to 8 m. According to the Dupuit formula (Harr 1962) The initial stable field can be directly applied, where the bottom and the upper part of the low water head on the right side of the slope are non-seepage boundaries. The final stable osmotic pressure distribution diagram is as follows according to Darcy's seepage law and the cubic law. On this basis, apply natural gravity, and at the same time set two displacement boundaries, where the left and right sides are lateral constraint boundaries, and the bottom is a vertical and lateral full constraint boundary. According to the stress-seepage characteristics, the gravity redistribution under the steady reservoir water level seepage field conditions can be obtained. Finally, wait until the above stress field and seepage field reach stability (about 10000 steps). Apply the heavy rainfall boundary, and at the same time consider the rainfall gravity effect and the fracture distribution. According to Darcy's law and the cubic law, the instantaneous rainfall seepage field can be obtained. According to the stress-seepage characteristics, the stress distribution under the instantaneous seepage field conditions can be obtained. At the same time, the gravitational acceleration in the slope is amplified by 20 times the original to accelerate the crack initiation and failure process of the slope. After continuing the simulation for 10000 steps. Considering the possible boundary truncation, all stress and seepage boundaries are three-layer, and both seepage and displacement boundaries can be superimposed and reused because they can be solved based on two different smoothing kernel functions.

[0147] It should be noted that the method of the embodiment of the present application can be executed by a single device, such as a computer or a server. The method of this embodiment can also be applied to a distributed scenario and completed by multiple devices cooperating with each other. In this case of a distributed scenario, one of the multiple devices can only execute one or more steps of the method of the embodiment of the present application, and these multiple devices will interact with each other to complete the described method.

[0148] It should be noted that some embodiments of the present application have been described above. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims can be executed in a different order from that in the above embodiments and still achieve the desired result. Additionally, the processes depicted in the drawings do not necessarily require the specific order or sequential order shown to achieve the desired result. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0149] Based on the same inventive concept, corresponding to the method of any of the above embodiments, the present application further provides a slope simulation device.

[0150] Refer to Figure 8 , the slope simulation device includes:

[0151] A parameter adjustment module 100, configured to obtain slope simulation parameters and assign values to a preset slope simulation model according to the slope simulation parameters to obtain an initial model.

[0152] A model update module 200, configured to solve the initial model according to the smoothed particle hydrodynamics method to obtain an updated model.

[0153] A stress calculation module 300, configured to calculate the particle state according to the updated model to obtain slope stress characteristics.

[0154] A result output module 400, configured to evaluate the slope stress characteristics according to a preset simulation rule to obtain a slope simulation result.

[0155] For the convenience of description, when describing the above device, it is divided into various modules according to functions and described separately. Of course, when implementing the present application, the functions of each module can be implemented in one or more software and / or hardware.

[0156] The device of the above embodiment is used to implement the corresponding slope simulation method in any of the foregoing embodiments and has the beneficial effects of the corresponding method embodiments, which will not be elaborated here.

[0157] Based on the same inventive concept, corresponding to the method of any of the above embodiments, the present application further provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, it implements the slope simulation method described in any one of the above embodiments.

[0158] Figure 10 FIG. shows a more specific schematic diagram of the hardware structure of the electronic device provided in this embodiment. The device may include: a processor 1010, a memory 1020, an input / output interface 1030, a communication interface 1040, and a bus 1050. Among them, the processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040 are communicatively connected to each other inside the device through the bus 1050.

[0159] The processor 1010 may be implemented in a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, etc., and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this specification.

[0160] The memory 1020 may be implemented in the form of a ROM (Read Only Memory), a RAM (Random Access Memory), a static storage device, a dynamic storage device, etc. The memory 1020 may store an operating system and other application programs. When implementing the technical solutions provided in the embodiments of this specification through software or firmware, the relevant program codes are stored in the memory 1020 and are called and executed by the processor 1010.

[0161] The input / output interface 1030 is used to connect to an input / output module to implement information input and output. The input / output module may be configured as a component in the device (not shown in the figure) or externally connected to the device to provide corresponding functions. Among them, the input device may include a keyboard, a mouse, a touch screen, a microphone, various sensors, etc., and the output device may include a display, a speaker, a vibrator, an indicator light, etc.

[0162] The communication interface 1040 is used to connect to a communication module (not shown in the figure) to implement communication interaction between this device and other devices. Among them, the communication module may communicate through a wired method (such as USB, network cable, etc.) or through a wireless method (such as a mobile network, WIFI, Bluetooth, etc.).

[0163] The bus 1050 includes a path for transmitting information among various components of the device, such as the processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040.

[0164] It should be noted that although the above device only shows the processor 1010, the memory 1020, the input / output interface 1030, the communication interface 1040, and the bus 1050, in the specific implementation process, the device may also include other components necessary for normal operation. In addition, those skilled in the art can understand that the above device may also only include the components necessary to implement the solution of the embodiments of this specification, and does not necessarily include all the components shown in the figure.

[0165] The electronic device of the above embodiment is used to implement the corresponding slope simulation method in any of the foregoing embodiments, and has the beneficial effects of the corresponding method embodiments, which will not be elaborated here.

[0166] Based on the same inventive concept, corresponding to the method of any of the above embodiments, the present application also provides a non-transitory computer-readable storage medium storing computer instructions for causing the computer to execute the slope simulation method as described in any of the above embodiments.

[0167] The computer-readable medium of this embodiment includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. The information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette tapes, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible by a computing device.

[0168] The computer instructions stored in the storage medium of the above embodiment are used to cause the computer to execute the slope simulation method as described in any of the above embodiments, and have the beneficial effects of the corresponding method embodiments, which will not be elaborated here.

[0169] It should be noted that the embodiments of the present application can also be further described in the following manner:

[0170] A slope simulation method, including:

[0171] Obtain slope simulation parameters, and assign values to a preset slope simulation model according to the slope simulation parameters to obtain an initial model.

[0172] Solve the initial model according to the smoothed particle hydrodynamics method to obtain an updated model.

[0173] Calculate the particle state according to the updated model to obtain the slope stress characteristics.

[0174] Evaluate the slope stress characteristics according to the preset simulation rules to obtain the slope simulation results. Among them, the simulation rules include failure criteria and contact criteria.

[0175] In some embodiments, solving the initial model according to the smoothed particle hydrodynamics method to obtain an updated model specifically includes:

[0176] Calculate the relationship between adjacent particles in the initial model to obtain initial particle information.

[0177] Optimize the initial particle information according to the basis spline function to obtain an updated model.

[0178] In some embodiments, calculating the relationship between adjacent particles in the initial model to obtain initial particle information specifically includes:

[0179] Set the smoothing kernel radius, and pair adjacent particles in the initial model according to the smoothing kernel radius.

[0180] Calculate the initial particle information according to the linked list search method. The initial particle information is calculated by the following formula:

[0181]

[0182] Among them, i represents the serial number of the base particle. j represents the serial number of the particle adjacent to the base particle. N represents the total number of particles in the initial model. f(x i ) represents the field function of the base particle i. represents the derivative of the field function of the base particle i. ρ represents the density of a single particle. m represents the mass of a single particle. h represents the smoothing kernel radius of the base particle i. W(x i -x j ,h) represents the smoothing kernel function. represents the gradient of the smoothing kernel function.

[0183] In some embodiments, calculating the particle state according to the updated model to obtain the slope stress characteristics specifically includes:

[0184] Calculate the force condition of the particles in the updated model to obtain the slope mechanical information.

[0185] Solve the slope mechanical information according to the preset water volume information to obtain the slope seepage information.

[0186] Calculate the slope stress characteristics based on the preset reservoir water level and slope seepage information.

[0187] In some embodiments, solve the slope mechanical information according to the preset water volume information to obtain the slope seepage information, specifically including:

[0188] Solve the slope mechanical information according to the cubic law to obtain the seepage coefficient.

[0189] Calculate the water level boundary, and calculate the steady-state seepage field according to the seepage coefficient and the water level boundary.

[0190] Calculate the precipitation boundary, and calculate the transient seepage field according to the steady-state seepage field and the precipitation boundary.

[0191] In some embodiments, evaluate the slope stress characteristics according to the preset simulation rules to obtain the slope simulation results, specifically including:

[0192] Evaluate the slope stress characteristics according to the failure criterion to obtain the particle failure information.

[0193] Calculate the particle position according to the particle failure information to obtain the particle position information.

[0194] Calculate the particle contact situation according to the contact criterion and the particle position information to obtain the slope simulation result.

[0195] In some embodiments, evaluate the slope stress characteristics according to the failure criterion to obtain the particle failure information, specifically including:

[0196] Evaluate the damage condition of the particles in the updated model according to the Drucker-Prager criterion and the slope stress characteristics to obtain the evaluation result.

[0197] Solve the particle control equation according to the evaluation result to obtain the particle failure information.

[0198] The particle control equation is calculated by the following formula:

[0199]

[0200] Wherein, represents the particle to be evaluated. W ij represents the smooth kernel function of the particle to be evaluated. ψ is a constant, and in response to the evaluation result that the particle to be evaluated is damaged, ψ = 0.

[0201] A slope simulation device, comprising:

[0202] A parameter adjustment module, configured to obtain slope simulation parameters, and assign values to a preset slope simulation model according to the slope simulation parameters to obtain an initial model.

[0203] A model update module, configured to solve the initial model according to the smoothed particle hydrodynamics method to obtain an updated model.

[0204] A stress calculation module, configured to calculate the particle state according to the updated model to obtain the slope stress characteristics.

[0205] A result output module, configured to evaluate the slope stress characteristics according to a preset simulation rule to obtain a slope simulation result.

[0206] An electronic device includes a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, the method described in any one of the above is implemented.

[0207] A non-transitory computer-readable storage medium stores computer instructions for causing a computer to execute the method described in any one of the above.

[0208] Those of ordinary skill in the art should understand that the discussion of any embodiment above is only exemplary and is not intended to imply that the scope of the present application (including the claims) is limited to these examples; under the idea of the present application, the technical features between the above embodiments or different embodiments can also be combined, the steps can be implemented in any order, and there are many other variations in different aspects of the embodiments of the present application as described above. For the sake of brevity, they are not provided in detail.

[0209] In addition, for the sake of simplicity of explanation and discussion, and in order not to make the embodiments of the present application difficult to understand, the known power / ground connections to integrated circuit (IC) chips and other components may or may not be shown in the provided drawings. In addition, the device may be shown in block diagram form to avoid making the embodiments of the present application difficult to understand, and this also takes into account the fact that the details of the implementation of these block diagram devices are highly dependent on the platform on which the embodiments of the present application are to be implemented (i.e., these details should be completely within the understanding of those skilled in the art). In the case where specific details (such as circuits) are set forth to describe the exemplary embodiments of the present application, it will be apparent to those skilled in the art that the embodiments of the present application can be implemented without these specific details or with variations of these specific details. Therefore, these descriptions should be considered illustrative rather than restrictive.

[0210] Although the present application has been described in connection with specific embodiments of the present application, many alternatives, modifications, and variations of these embodiments will be apparent to those of ordinary skill in the art based on the foregoing description. For example, other memory architectures (such as dynamic RAM (DRAM)) may be used with the embodiments discussed.

[0211] Embodiments of the present application are intended to cover all such substitutions, modifications, and variations that fall within the broad scope of the appended claims. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiments of the present application shall be included within the protection scope of the present application.

Claims

1. A slope simulation method, comprising: Obtaining slope simulation parameters, and assigning values to a preset slope simulation model according to the slope simulation parameters to obtain an initial model; Solving the initial model according to the smoothed particle hydrodynamics method to obtain an updated model; Calculating the particle state according to the updated model to obtain slope stress characteristics; Evaluating the slope stress characteristics according to preset simulation rules to obtain a slope simulation result; wherein, the simulation rules include a failure criterion and a contact criterion.

2. The slope simulation method according to claim 1, wherein, The step of solving the initial model according to the smoothed particle hydrodynamics method to obtain an updated model specifically includes: Calculating the relationship between adjacent particles in the initial model to obtain initial particle information; Optimizing the initial particle information according to the basis spline function to obtain the updated model.

3. The slope simulation method according to claim 2, wherein, The step of calculating the relationship between adjacent particles in the initial model to obtain initial particle information specifically includes: Setting a smoothing kernel radius, and pairing adjacent particles in the initial model according to the smoothing kernel radius; Calculating the initial particle information according to the linked list search method; the initial particle information is calculated by the following formula: Among them, i represents the serial number of the base particle; j represents the serial number of the particle adjacent to the base particle; N represents the total number of particles in the initial model; f(x i ) represents the field function of the base particle i; represents the derivative of the field function of the base particle i; ρ represents the density of a single particle; m represents the mass of a single particle; h represents the smoothing kernel radius of the base particle i; W(x i -x j , h) represents the smoothing kernel function; represents the gradient of the smoothing kernel function.

4. The slope simulation method according to claim 1, wherein, The step of calculating the particle state according to the updated model to obtain slope stress characteristics specifically includes: Calculating the force condition of particles in the updated model to obtain slope mechanical information; Solving the slope mechanical information according to preset water volume information to obtain slope seepage information; Calculating the slope stress characteristics according to the preset reservoir water level and the slope seepage information.

5. The slope simulation method according to claim 4, wherein, The step of solving the slope mechanical information according to the preset water volume information to obtain slope seepage information specifically includes: Solving the slope mechanical information according to the cubic law to obtain a seepage coefficient; Calculating a water level boundary, and calculating a steady-state seepage field according to the seepage coefficient and the water level boundary; Calculating a precipitation boundary, and calculating a transient seepage field according to the steady-state seepage field and the precipitation boundary.

6. According to the slope simulation method described in claim 1, the step of evaluating the slope stress characteristics according to preset simulation rules to obtain a slope simulation result specifically includes: Evaluating the slope stress characteristics according to the failure criterion to obtain particle failure information; Calculating the particle position according to the particle failure information to obtain particle position information; Calculating the particle contact condition according to the contact criterion and the particle position information to obtain the slope simulation result.

7. The slope simulation method according to claim 6, wherein, The step of evaluating the slope stress characteristics according to the failure criterion to obtain particle failure information specifically includes: Evaluating the damage condition of particles in the updated model according to the Drucker-Prager criterion and the slope stress characteristics to obtain an evaluation result; Solving the particle control equation according to the evaluation result to obtain the particle failure information; The particle control equation is calculated by the following formula: Among them, represents the particle to be evaluated; W ij represents the smoothing kernel function of the particle to be evaluated; ψ is a constant, and in response to the evaluation result that the particle to be evaluated is damaged, ψ = 0.

8. A slope simulation device, comprising: A parameter adjustment module, configured to obtain slope simulation parameters, and assign values to a preset slope simulation model according to the slope simulation parameters to obtain an initial model; A model update module, configured to solve the initial model according to the smoothed particle hydrodynamics method to obtain an updated model; A stress calculation module, configured to calculate the particle state according to the updated model to obtain the slope stress characteristics; A result output module, configured to evaluate the slope stress characteristics according to a preset simulation rule to obtain a slope simulation result.

9. An electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein when the processor executes the program, the method according to any one of claims 1 to 7 is implemented.

10. A non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the method according to any one of claims 1 to 7.

Images