Optimization design method, device and equipment for three-dimensional slope inclined slide-resistant pile and medium

Abstract

The invention relates to the field of computers, and provides an optimal design method, device and equipment for a three-dimensional slope inclined anti-slide pile and a medium. The method comprises the following steps: generating a dynamic library file of a soil constitutive model, and importing the dynamic library file into universal interactive modeling software; setting material attributes for the three-dimensional model, defining boundary conditions, and dividing grids; performing static analysis on the three-dimensional model to generate a visual toolkit file; extracting material point information in the visual toolkit file, and packaging the material point information into a material point file; determining pre-prepared material point method three-dimensional data analysis software; reading the material point file and carrying out dynamic analysis; generating an unstructured grid file of the visual toolkit file in the dynamic analysis process; extracting a slope failure state variable from the unstructured grid file of the visual toolkit file; and by taking the slope failure state variable as an optimization target, adjusting the three-dimensional model to obtain an optimal anti-slide pile design scheme.

Description

Technical Field

[0001] The present application relates to the field of computers, and more specifically, to an optimization design method, device, equipment and medium for three-dimensional slope-inclined anti-slide piles. Background Art

[0002] Landslide is the most widespread natural disaster. In order to improve the stability of the slope, reduce the possibility of landslide induction and the degree of damage after landslide, anti-slide piles are often used in engineering to reinforce the slope to improve the slope safety factor and control landslide disasters.

[0003] In recent years, domestic and foreign scholars have often used the finite element method and the finite difference method to analyze the stability of slopes reinforced with anti-slide piles. However, previous studies on the impact of anti-slide pile design on slopes have mainly focused on the slope stability stage. For large deformations of slopes, the disaster reduction effect of anti-slide pile design on the slope damage stage is rarely considered. Although improving the stability of slopes is of great significance to engineering protection, it is also particularly important to conduct risk assessment on landslide disasters managed by anti-slide piles after slope damage. However, the pile foundation reinforced slopes studied based on traditional methods are constrained to solve small deformation problems. The traditional method is used to numerically solve large deformation problems, which will cause the calculation to not converge due to mesh distortion, and it is impossible to capture the entire process of large deformation of the slope, thereby ignoring the effect of anti-slide pile design on slope damage management.

[0004] Therefore, the present application provides an optimization design method, device, equipment and medium for three-dimensional slope inclined anti-slide piles to solve one of the above-mentioned technical problems. Summary of the invention

[0005] The purpose of this application is to provide a method, device, equipment and medium for optimizing the design of three-dimensional slope anti-slide piles, which can solve at least one of the technical problems mentioned above. The specific solution is as follows: According to the specific implementation of the present application, in the first aspect, the present application provides an optimization design method for three-dimensional slope inclined anti-slide piles, comprising: Step S1: Based on the soil constitutive model generation software, a dynamic library file of the soil constitutive model is generated, and the dynamic library file is imported into the general interactive modeling software; Step S2: Based on the user interface of the general interactive modeling software, a three-dimensional model of the slope and the anti-slide pile is established, material properties are set for the three-dimensional model, boundary conditions are defined, and a grid is divided; Step S3: Based on the integrated environment of the general interactive modeling software and the soil constitutive model generation software, a static analysis is performed on the three-dimensional model to generate a visualization toolkit file; Step S4: Based on the data analysis visualization software, material point information in the visualization toolkit file is extracted, and the material point information is packaged into a material point file; Step S5: Determine the pre-prepared material point method three-dimensional data analysis software, and judge whether the size of the grid is consistent with the material point method three-dimensional data analysis software. The dimensions defined in the data analysis software are consistent. If they are inconsistent, return to step S2 and reset them, otherwise execute step S6; step S6: read the material point file based on the material point method three-dimensional data analysis software and perform dynamic analysis; step S7: based on the general interactive modeling software, soil constitutive model generation software and material point method three-dimensional data analysis software integrated environment, use the strength reduction method to generate an unstructured grid file of the visualization toolkit file during the dynamic analysis process; step S8: extract the slope failure state variable from the unstructured grid file of the visualization toolkit file through the data analysis visualization software, which is used to evaluate the disaster reduction effect of different anti-slide pile design schemes; step S9: take the slope failure state variable as the optimization target, adjust the three-dimensional model, and obtain the optimal anti-slide pile design scheme.

[0006] In one embodiment, adjusting the three-dimensional model includes: adjusting the length and inclination angle of the anti-slip pile.

[0007] In one embodiment, the method further comprises: simulating a large deformation process of the slope using a material point method; wherein a simulation result of the large deformation process simulation is used together with the failure state variable of the slope as an optimization target.

[0008] In one implementation, the slope failure state variables include energy, velocity and sliding distance.

[0009] According to a specific implementation of the present application, in a second aspect, the present application provides an optimization design device for three-dimensional slope inclined anti-slide piles, comprising: A generating unit is used to execute step S1, step S3, and step S7; the step S1: based on the soil constitutive model generating software, a dynamic library file of the soil constitutive model is generated, and the dynamic library file is imported into the general interactive modeling software; the step S3: based on the integrated environment of the general interactive modeling software and the soil constitutive model generating software, a static analysis is performed on the three-dimensional model to generate a visualization toolkit file; the step S7: based on the integrated environment of the general interactive modeling software, the soil constitutive model generating software and the material point method three-dimensional data analysis software, a strength reduction method is used to generate an unstructured grid file of the visualization toolkit file in the dynamic analysis process; the processing unit is used to execute step S2, step S4, step S6, step S8, and step S9; the step S2: based on the user interface of the general interactive modeling software, a three-dimensional model of the slope and the anti-slide pile is established, material properties are set for the three-dimensional model, boundary conditions are defined, and The method comprises the following steps: first, extracting the material point information in the visualization toolkit file based on the data analysis visualization software, and packaging the material point information into a material point file; second, reading the material point file based on the material point method three-dimensional data analysis software, and performing dynamic analysis; third, extracting the slope failure state variable from the unstructured grid file of the visualization toolkit file through the data analysis visualization software, for evaluating the disaster reduction effect of different anti-slide pile design schemes; fourth, adjusting the three-dimensional model with the slope failure state variable as the optimization target to obtain the optimal anti-slide pile design scheme; and a judgment unit for executing step S5; fourth, determining the pre-prepared material point method three-dimensional data analysis software, and judging whether the size of the grid is consistent with the size defined in the material point method three-dimensional data analysis software. If not, returning to step S2 to reset it, otherwise executing step S6.

[0010] In one implementation, the processing unit adjusts the three-dimensional model in the following manner: adjusting the length and inclination of the anti-slide pile.

[0011] In one embodiment, the processing unit is further used to: simulate the large deformation process of the slope using a material point method; wherein the simulation result of the large deformation process simulation is used together with the slope failure state variable as an optimization target.

[0012] In one implementation, the slope failure state variables include energy, velocity and sliding distance.

[0013] According to the specific implementation of the present application, in a third aspect, the present application provides an electronic device, comprising a memory, a processor, and a computer program stored on the memory, characterized in that the processor executes the computer program to implement any one of the methods described in the first aspect.

[0014] According to the specific implementation of the present application, in a fourth aspect, the present application provides a computer-readable storage medium having a computer program / instruction stored thereon, characterized in that when the computer program / instruction is executed by a processor, the method described in any one of the first aspects is implemented.

[0015] Compared with the prior art, the above solution of the embodiment of the present application has at least the following beneficial effects: The present application provides an optimization design method for three-dimensional slope inclined anti-slide piles. By integrating soil constitutive model generation software, material point method three-dimensional data analysis software, and general interactive modeling software, the material point method is combined to perform systematic optimization design of three-dimensional slope inclined anti-slide piles. The slope model can be generated and analyzed efficiently and accurately, and the stability of the slope can be comprehensively evaluated from static to dynamic, ensuring the reliability and effectiveness of the design, thereby significantly improving the stability and safety of the slope and reducing the risk and loss of landslide disasters. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] Figure 1 A flow chart showing an optimization design method for three-dimensional slope anti-slide piles; Figure 2 A schematic diagram of modeling from point to line to surface and then to body is shown; Figure 3 A schematic diagram of a modeling setting interface is exemplarily shown; Figure 4 A grid schematic diagram is exemplarily shown; Figure 5 A schematic diagram of a material point is shown as an example; Figure 6 A unit block diagram of a device for optimizing the design of three-dimensional slope anti-slide piles according to an embodiment of the present application is shown; Figure 7 It is a block diagram of an electronic device for optimizing the design of three-dimensional slope anti-slide piles according to an exemplary embodiment. DETAILED DESCRIPTION

[0017] In order to make the purpose, technical solutions and advantages of the present application clearer, the present application will be further described in detail below in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present application.

[0018] The terms used in the embodiments of the present application are only for the purpose of describing specific embodiments, and are not intended to limit the present application. The singular forms "a", "said" and "the" used in the embodiments of the present application and the appended claims are also intended to include plural forms, unless the context clearly indicates other meanings, and "multiple" generally includes at least two.

[0019] It should be understood that the term "and / or" used in this article is only a description of the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in this article generally indicates that the associated objects before and after are in an "or" relationship.

[0020] It should be understood that although the terms first, second, third, etc. may be used to describe in the embodiments of the present application, these descriptions should not be limited to these terms. These terms are only used to distinguish the descriptions. For example, without departing from the scope of the embodiments of the present application, the first may also be referred to as the second, and similarly, the second may also be referred to as the first.

[0021] As used herein, the words "if" and "if" may be interpreted as "at the time of" or "when" or "in response to determining" or "in response to detecting", depending on the context. Similarly, the phrases "if it is determined" or "if (stated condition or event) is detected" may be interpreted as "when it is determined" or "in response to determining" or "when detecting (stated condition or event)" or "in response to detecting (stated condition or event)", depending on the context.

[0022] It should also be noted that the term "includes", "comprising" or any other variation thereof is intended to cover non-exclusive inclusion, so that a commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprising a ..." do not exclude the existence of other identical elements in the commodity or device including the elements.

[0023] It should be particularly noted that any symbols and / or numbers in the specification that are not marked in the accompanying drawings are not drawing marks.

[0024] In order to make it easier to clarify the meaning of the Chinese interpretations in the claims, the English abbreviations or full names are shown below for each interpretation.

[0025] Soil constitutive model generation software (Anura3d); Material point method 3D data analysis software (MPM3D); General Interactive Modeling Software (GID); General interactive modeling software and soil constitutive model generation software integrated environment (GID-Anura3d); General interactive modeling software, soil constitutive model generation software and material point method 3D data analysis software integrated environment (GID-Anura3d-MPM3D); Dynamic library files (dll files); Visualization Toolkit files (vtk files); Data analysis and visualization software (ParaView); Material point files (mpm files).

[0026] In addition, the following embodiments of the present application may also involve some special English terms that are not shown in the figure. For example, MATLAB can be used to package material point information into a material point file mpm. Such contents whose corresponding explanations are not shown are all special terms in this field and can be directly referred to the explanations specified in this field.

[0027] The optional embodiments of the present application are described in detail below with reference to the accompanying drawings.

[0028] The embodiment provided in this application is an embodiment of an optimization design method for three-dimensional slope inclined anti-slide piles.

[0029] Combine the following Figure 1 The embodiments of the present application are described in detail.

[0030] Figure 1 A flowchart of an optimization design method for three-dimensional slope anti-slide piles is shown. Figure 1 As shown, the steps include the following steps S1 to S9.

[0031] Step S1: Based on the soil constitutive model generation software, a dynamic library file of the soil constitutive model is generated, and the dynamic library file is imported into the general interactive modeling software.

[0032] Step S2: Building a three-dimensional model of the slope and the anti-slide piles based on the user interface of the general interactive modeling software, setting material properties for the three-dimensional model, defining boundary conditions, and dividing the mesh.

[0033] In this application, the constructed three-dimensional model is a pile foundation reinforced three-dimensional slope model.

[0034] Step S3: Based on the integrated environment of the general interactive modeling software and the soil constitutive model generation software, a static analysis is performed on the three-dimensional model to generate a visualization toolkit file.

[0035] Step S4: extracting material point information in the visualization toolkit file based on the data analysis visualization software, and packaging the material point information into a material point file.

[0036] Step S5: Determine the pre-prepared material point method three-dimensional data analysis software, and judge whether the size of the grid is consistent with the size defined in the material point method three-dimensional data analysis software. If not, return to step S2 and reset it, otherwise execute step S6.

[0037] Step S6: The material point method-based three-dimensional data analysis software reads the material point file and performs dynamic analysis.

[0038] Step S7: Based on the integrated environment of the general interactive modeling software, the soil constitutive model generation software and the material point method three-dimensional data analysis software, the strength reduction method is used to generate an unstructured grid file of the visualization toolkit file during the dynamic analysis process.

[0039] Step S8: extracting slope failure state variables from the unstructured grid file of the visualization toolkit file through data analysis visualization software to evaluate the disaster reduction effects of different anti-slide pile design schemes.

[0040] Step S9: Taking the slope failure state variable as the optimization target, the three-dimensional model is adjusted to obtain the optimal anti-slide pile design solution.

[0041] In one embodiment of the present application, by generating a dynamic link library file of a soil constitutive model and importing it into a general interactive modeling software, the accuracy and flexibility of the soil constitutive model can be ensured, and the reliability and computational efficiency of the model can be improved. For example, Anura3d is used to generate a dll file of a soil constitutive model and import it into GID. The specific steps include compiling the Anura3d code, generating dll files of a Mohr-Coulomb model and a linear elastoplastic model, and placing these files in a specified path of the GID.

[0042] In one embodiment of the present application, by establishing a three-dimensional model of the slope and anti-slide piles in a general interactive modeling software, setting material properties, defining boundary conditions and dividing the grid, the accuracy and integrity of the model can be ensured, providing a reliable basis for subsequent analysis. For example, in the GID user interface, points, lines, surfaces and bodies are created in sequence to construct a three-dimensional model of the slope and anti-slide piles. Material parameters (such as density, elastic modulus, Poisson's ratio, etc.) are assigned to the model, boundary conditions (such as fixed boundaries, free boundaries, etc.) are defined, and the grid is divided.

[0043] In one embodiment of the present application, by statically analyzing the three-dimensional model and generating a vtk file, the rationality of the model can be preliminarily verified, providing a reference for subsequent dynamic analysis. For example, in the GID-Anura3d integrated environment, a static analysis is performed on the established three-dimensional model to generate a vtk file, which contains the geometric information of the model and the preliminary analysis results, which can be used for subsequent data processing and visualization.

[0044] In one embodiment of the present application, the material point information in the vtk file is extracted by data analysis visualization software and packaged into an mpm file, which can ensure the integrity and accuracy of the material point information and provide a reliable data basis for dynamic analysis. For example, ParaView is used to read the vtk file, extract the serial number, coordinates, mass, material number and other information of the material point, and package this information into an mpm file through MATLAB.

[0045] In one embodiment of the present application, the accuracy and reliability of subsequent dynamic analysis are ensured by judging the consistency of the grid size, and calculation errors caused by inconsistent grid sizes are avoided. For example, the pre-prepared MPM3D software is determined, and it is checked whether the grid size divided in the GID is consistent with the size defined in the MPM3D. If not, return to step S2 to reset the grid size to ensure that the two are consistent.

[0046] In one embodiment of the present application, the mpm file is read and dynamic analysis is performed by material point method three-dimensional data analysis software, which can accurately simulate the behavior of the slope under different conditions and provide a scientific basis for optimization design. For example, MPM3D is used to read the generated mpm file and perform dynamic analysis. MPM3D will output a vtu file containing the results of the dynamic analysis.

[0047] In one embodiment of the present application, by using the strength reduction method to generate a vtu file during the dynamic analysis process, the stability of the slope under different strength conditions can be evaluated, providing detailed analysis data for optimal design. For example, in the GID-Anura3d-MPM3D integrated environment, the strength reduction method is used to reduce the cohesion and internal friction angle of the soil parameters, induce slope landslide instability, and generate a vtu file.

[0048] In one embodiment of the present application, by extracting the slope failure state variables in the vtu file through data analysis visualization software, the disaster reduction effect of different anti-slide pile design schemes can be quantitatively evaluated, providing a scientific basis for optimizing the design. For example, ParaView is used to read the vtu file, extract the slope failure state variables (such as energy, speed, sliding distance, etc.), and analyze and compare them.

[0049] In one embodiment of the present application, the failure state variables of the slope are used as optimization targets, and the three-dimensional model is adjusted to find the optimal anti-slide pile design solution, thereby improving the stability and safety of the slope. For example, according to the extracted failure state variables of the slope, the length and inclination of the anti-slide pile are adjusted, and the above steps are repeated until the optimal anti-slide pile design solution is found.

[0050] The method provided in this application can efficiently and accurately optimize the design of slope tilt anti-slide piles by integrating GID, Anura3d and MPM3D, thereby improving the stability and safety of the slope. Among them, a three-dimensional model of the slope and anti-slide piles is established using GID, and a dll file of the soil constitutive model is generated by Anura3d and imported into GID. Static analysis is performed in GID to generate a vtk file, and ParaView is used to extract material point information and generate an mpm file. MPM3D reads the mpm file for dynamic analysis and generates a vtu file. Finally, ParaView is used to extract the slope failure state variables in the vtu file to evaluate the disaster reduction effects of different anti-slide pile design schemes. In addition, by adjusting the length and inclination of the anti-slide piles, the optimal anti-slide pile design scheme can be found to further improve the stability and safety of the slope. A three-dimensional model of the slope and anti-slide piles is established in GID, the length and inclination of the anti-slide piles are adjusted, and the above steps are repeated until the optimal design scheme is found. For example, different pile lengths (such as 8m, 9m, 10m) and inclination angles (such as 0°, 30°, 60°) can be set, and the optimal combination of pile length and inclination angle can be determined through dynamic analysis and extraction of failure state variables. Furthermore, by simulating the large deformation process of the slope through the material point method, the behavior of the slope under extreme conditions can be more accurately evaluated, providing more comprehensive data support for optimization design. The material point method is used in MPM3D to simulate the large deformation process of the slope and generate a vtu file. The large deformation process data in the vtu file is extracted and used together with the slope failure state variables as optimization targets to evaluate the effects of different anti-slide pile design schemes. Furthermore, by extracting the slope failure state variables (energy, velocity and sliding distance), the disaster reduction effects of different anti-slide pile design schemes can be quantitatively evaluated, providing a scientific basis for optimization design. Use ParaView to read the vtu file and extract the slope failure state variables (energy, velocity and sliding distance). By comparing these variables of different design schemes, the optimal anti-slide pile design scheme is determined. For example, the energy consumption, sliding distance and speed under different combinations of pile lengths and inclination angles can be compared to select the solution with the best disaster reduction effect.

[0051] In this application, the material point method can be used to analyze the large deformation of the slope.

[0052] In this application, the principle theoretical formula involved in step S7 is as follows: In the embodiment of the present application, the updated Lagrangian format is used to describe the motion of the object. During the entire motion process, heat exchange is not considered, and the control equations must satisfy the three major equations of conservation of mass, conservation of momentum, and conservation of energy: The mass conservation equation: ; Momentum conservation equation: ; Energy conservation equation: .

[0053] Among them, ρ, v, , e is the object's density, velocity, acceleration, and internal energy. bi is the force acting on the unit mass of the object. Dij and are the deformation rate tensor and the Cauchy stress tensor, respectively.

[0054] In this application, the following formula is used for the boundary conditions: ; ; Among them, n j is a unit vector representing the direction of the boundary normal. and Represent the boundary force and boundary velocity respectively. and They represent the given traction boundary and the given displacement boundary respectively.

[0055] In the embodiment of the present application, for the discrete material points, the continuum is in the solution domain. The inside is discretized into a group of mass points with mass. During the whole motion process, the mass carried by the mass points remains unchanged, so the mass conservation equation automatically satisfies its conditions. As for the momentum conservation equation, partial differential equations bring difficulties to direct numerical calculation in solving complex practical problems. Therefore, the weighted residual method is used to establish a weak form equivalent to the original differential equation and the boundary condition, and then the solution is approximated on this basis. Introducing virtual displacement , the displacement is 0 on the boundary condition, and the specific stress and specific boundary surface force are introduced , , combined with the given boundary conditions and the momentum conservation equation, the virtual work equation is obtained: ; Among them, the density of the continuum discretized into several particles can be approximated as: ; Where n is the total number of material points and mp is the mass of material point p. is the Dirac function. xp is the coordinate of the material point p. Transforming the set yields: ; Among them, h is the imaginary boundary layer thickness used to convert the surface integral on the left into a volume integral. By discretizing this relationship, the integral form is converted into the sum of the product of the value of the integrand at each particle and the volume represented by the particle.

[0056] In this application, during the entire process of solving the momentum, the mass point and the background grid are completely fixed without relative movement, and the mutual mapping of information between the mass point and the background grid is completed by finite element shape function Na(xi) interpolation: ; in, represents the coordinates of particle p at time t. Represents the variable of the background grid node. Na(xp) is the interpolation value of the shape function on node a at mass point p. Represents the coordinates of node a at time t. Through such interpolation relationship, the motion equation of the background grid node can be solved and converted into: ; Among them, m a is the node concentrated mass matrix, . is the acceleration of node a at time t. and are the internal force and external force of node a at time t respectively.

[0057] In this application, the optimization pile design ideas are as follows: Firstly, according to the actual slope project, different pile lengths and pile inclination angles are set independently.

[0058] Secondly, through strength reduction, slope landslide is achieved.

[0059] Finally, ParaView was used to extract the slope failure state variables (energy, velocity, and sliding distance) in the vtu file to compare and analyze different pile lengths, pile inclinations, and the combined effects of pile lengths and pile inclinations to achieve the optimal value for slope landslide disaster reduction losses, and the design parameters were used as the optimal pile design parameters.

[0060] In this application, the material point method theory of dual description of Lagrangian particles and Euler grids is adopted, and the GID-Anura3d-MPM3D combination method is adopted to simulate the influence of piles, different pile inclinations and pile lengths on the sliding distance of slope instability. The degree of damage to the slope is analyzed from the perspective of energy, and the difference in the degree of slope instability caused by inclined anti-sliding piles of different angles and lengths is compared and analyzed. For example, the sliding distance and deformation energy caused by large deformation of the slope can be effectively reduced by increasing the anti-sliding piles. When the inclination angle is 0°, as the pile length increases, the sliding distance, the horizontal velocity of the top and the foot of the slope all decrease. If the pile length is not properly selected, the anti-sliding piles may be unable to resist the energy accumulated during the deformation of the slope, thereby generating a greater impact. At the same time, the ability of the inclined anti-sliding piles to reinforce the slope to resist large deformation of the slope is greater than that of the vertical piles. In this application, under the condition of a pile length of 9m, the sliding distance, the top displacement, and the top and foot velocities of the slope generated by the inclined piles are all smaller than those generated by the vertical piles. Under the condition of 9m anti-slide piles, when the inclination angle is 50°, its reinforcement effect is better than that of 10m anti-slide piles for vertically reinforced slopes. Therefore, considering the relationship between pile length and inclination angle can not only effectively reduce disasters, but also reduce engineering costs.

[0061] In addition, in terms of slope disaster reduction, considering that the pile length is 9m, as the inclination angle increases from 0° to 70°, the sliding distance, deformation energy and velocity generated by the slope show a trend of first decreasing and then increasing. When the inclination angle is 60°, the slope disaster reduction effect is the best. In terms of slope stability maintenance: under the condition of a pile length of 10m, the inclination angle of the anti-sliding pile is in the range of 20° to 70°, and the sliding distance generated by the slope is reduced to within 1m. At 60°, the sliding distance, velocity and deformation energy generated by the slope are still very small. Therefore, without considering the pile length, when the inclination angle of the anti-sliding pile is 60°, the slope achieves the best disaster prevention and reduction effect.

[0062] In addition, the implementation details of the above-mentioned embodiment based on the use of software are described below in conjunction with the diagrams.

[0063] Step S1: (1) Anura3d can use foreign open source code, compile it, and generate the soil constitutive model Mohr-Coulomb model dll file and linear elastoplastic model dll file. (2) Put the generated dll file into the GID installation path problemtypes / Anura3d_2022.gid / dll file, so that the GID software can use the soil constitutive model. (3) GID is a general, adaptive and user-friendly application developed for geometric modeling, data input and visualization of result files of various numerical analysis programs. GID can be used as a pre-processing modeling and post-processing visualization analysis tool for calculation schemes such as finite element method, finite volume method, boundary element method, finite difference method or meshless algorithm. GID software has comprehensive geometric modeling, meshing, CAD data import, post-processing result display and other functions. GID adopts an operation mode similar to CAD, and it also has powerful batch processing capabilities. Users will feel unprecedented convenience and ease when using GID to create complex model problems. It is easy to operate, convenient and flexible, and intuitive and convenient. GID can handle modeling and simulation problems in various fields well, including solid and structural mechanics, fluid dynamics, electromagnetics, heat conduction and geomechanics.

[0064] Step S2: (1) Modeling is performed in the GID user interface. (2) The modeling sequence is from point to line to surface and then to volume, such as Figure 2 (3) Assign material parameters to the model, define boundary conditions, divide the grid, perform calculations, etc. Figure 3 shown.

[0065] Steps S3 and S4: (1) After calculation, a vtk file will be generated. The vtk model is a 3D model file represented in text format. It can represent point and surface information and can be stored in text format in a way that is easy for humans to read and understand. (2) ParaView is an open source, cross-platform data analysis and visualization program that analyzes and visualizes two-dimensional and three-dimensional data and supports file input in vtk format. (3) Use ParaView to read the vtk file and extract the serial number and coordinate information of the material point. Figure 4 A grid schematic diagram is shown exemplarily, based on Figure 4 Up to Grid Size, Figure 5 shows a schematic diagram of a material point, based on Figure 5 It can be seen that the model corresponds to each discrete material point. (4) The mpm file is automatically generated by packaging the data through MATLAB.

[0066] Step S5: (1) MPM3D is adapted from open source code, and the code can only read and run .mpm type files. (2) The content format of the mpm file can use the "dcell" item to define the network size, and define it to be consistent with the network size in step S4 to improve the calculation accuracy. In addition, different parameters can be labeled by the "num" item. For example, the "num" item is labeled with items "1" and "2", where item "1" is the soil parameter and item "2" is the anti-sliding pile parameter. (3) MPM3D reads the mpm file for calculation and generates a VTU file. VTU is the core file for data storage, which contains all the information of nodes, units, and result values. On this basis, visualization is performed using the Paraview tool.

[0067] Step S7: By using the strength reduction method, the soil parameters cohesion and internal friction angle are reduced to induce slope instability. Assuming the reduction factor k0, the soil parameters cf and for: ; ; Among them, c0 and is the initial cohesion and internal friction angle of the soil, k0 is the reduction factor k0, and also expresses the safety factor of the slope until it reaches its limit state.

[0068] Step S8: (1) Extract the VTU file generated by MPM3D calculation. (2) Use ParaView to read the VTU file, and you can select the state information of the point independently to extract information such as speed change, energy, real-time coordinates, and displacement change.

[0069] This scheme integrates GID, Anura3d and MPM3D software, and combines the material point method (MPM) to optimize the design of three-dimensional slope tilt anti-slide piles, which can efficiently and accurately simulate the behavior of slopes under different conditions, especially large deformation processes. By adjusting the length and inclination of the anti-slide piles and extracting slope failure state variables (such as energy, velocity and sliding distance), the disaster reduction effects of different design schemes can be quantitatively evaluated, thereby finding the optimal anti-slide pile design scheme, significantly improving the stability and safety of the slope, and reducing the risk and loss of landslide disasters.

[0070] The present application also provides an apparatus embodiment that is consistent with the above embodiment, which is used to implement the method steps of the above embodiment. The explanation based on the same name meaning is the same as the above embodiment, and has the same technical effect as the above embodiment, which will not be repeated here.

[0071] like Figure 6 As shown, the present application provides an optimization design device 600 for three-dimensional slope anti-slide piles, comprising: The generation unit 601 is used to execute step S1, step S3, and step S7. Step S1: Based on the soil constitutive model generation software, a dynamic library file of the soil constitutive model is generated, and the dynamic library file is imported into the general interactive modeling software. Step S3: Based on the integrated environment of the general interactive modeling software and the soil constitutive model generation software, a static analysis is performed on the three-dimensional model to generate a visualization toolkit file. Step S7: Based on the integrated environment of the general interactive modeling software, the soil constitutive model generation software, and the material point method three-dimensional data analysis software, an unstructured grid file of the visualization toolkit file is generated in the dynamic analysis process using the strength reduction method. The processing unit 602 is used to execute step S2, step S4, step S6, step S8, and step S9. Step S2: Based on the user interface of the general interactive modeling software, a three-dimensional model of the slope and the anti-slide pile is established, material properties are set for the three-dimensional model, boundary conditions are defined, and the mesh is divided. Step S4: Based on the data analysis visualization software, the material point information in the visualization toolkit file is extracted, and the material point information is packaged into a material point file. Step S6: Read the material point file based on the material point method three-dimensional data analysis software and perform dynamic analysis. Step S8: Extract the slope failure state variables from the unstructured grid file of the visualization toolkit file through the data analysis visualization software, which is used to evaluate the disaster reduction effect of different anti-slide pile design schemes. Step S9: Take the slope failure state variables as the optimization target and adjust the three-dimensional model to obtain the optimal anti-slide pile design scheme. Judgment unit 603 is used to execute step S5. Step S5: Determine the pre-prepared material point method three-dimensional data analysis software, and determine whether the size of the grid is consistent with the size defined in the material point method three-dimensional data analysis software. If not, return to step S2 to reset it, otherwise execute step S6.

[0072] In one implementation, the processing unit 602 adjusts the three-dimensional model in the following manner: adjusting the length and inclination of the anti-sliding piles.

[0073] In one embodiment, the processing unit 602 is further used to: simulate the large deformation process of the slope using a material point method, wherein the simulation result of the large deformation process simulation is used together with the slope failure state variable as an optimization target.

[0074] In one embodiment, the slope failure state variables include energy, velocity and sliding distance.

[0075] Regarding the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0076] Figure 7 It is a block diagram of an electronic device 700 for optimizing the design of three-dimensional slope anti-slide piles according to an exemplary embodiment.

[0077] like Figure 7 As shown, one embodiment of the present application provides an electronic device 700. The electronic device 700 includes a memory 701, a processor 702, and an input / output (I / O) interface 703. The memory 701 is used to store instructions. The processor 702 is used to call the instructions stored in the memory 701 to execute the method for optimizing the design of three-dimensional slope inclined anti-sliding piles in the embodiment of the present application. The processor 702 is connected to the memory 701 and the I / O interface 703, respectively, for example, through a bus system and / or other forms of connection mechanisms (not shown). The memory 701 can be used to store programs and data, including the program of the method for optimizing the design of three-dimensional slope inclined anti-sliding piles involved in the embodiment of the present application. The processor 702 executes various functional applications and data processing of the electronic device 700 by running the program stored in the memory 701.

[0078] In the embodiment of the present application, the processor 702 can be implemented in at least one hardware form of a digital signal processor (Digital Signal Processing, DSP), a field programmable gate array (Field Programmable Gate Array, FPGA), and a programmable logic array (Programmable Logic Array, PLA). The processor 702 can be a central processing unit (Central Processing Unit, CPU) or one or a combination of other forms of processing units with data processing capabilities and / or instruction execution capabilities.

[0079] The memory 701 in the embodiment of the present application may include one or more computer program products, and the computer program product may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory (cache), etc. The non-volatile memory may include, for example, read-only memory (ROM), flash memory, hard disk drive (HDD), or solid state drive (SSD), etc.

[0080] In the embodiment of the present application, the I / O interface 703 may be used to receive input instructions (such as digital or character information, and key signal input related to user settings and function control of the electronic device 700), and may also output various information to the outside (such as images or sounds). In the embodiment of the present application, the I / O interface 703 may include one or more of a physical keyboard, a function key (such as a volume control key, a switch key, etc.), a mouse, a joystick, a trackball, a microphone, a speaker, and a touch panel.

[0081] In some embodiments, the present application provides a computer-readable storage medium, which stores computer-executable instructions. When the computer-executable instructions are executed by a processor, any of the methods described above is performed.

[0082] In some embodiments, the present application provides a computer program product, which includes a computer program. When the computer program is executed by a processor, any of the methods described above is performed.

[0083] Although operations are described in a particular order in the drawings, this should not be understood as requiring that the operations be performed in the particular order shown or in serial order, or that all shown operations be performed to achieve desired results. In certain circumstances, multitasking and parallel processing may be advantageous.

[0084] The methods, devices, equipment and media of the present application can be implemented using standard programming techniques, using rule-based logic or other logic to implement various method steps. It should also be noted that the words "device" and "module" used herein and in the claims are intended to include implementations using one or more lines of software code and / or hardware implementations and / or devices for receiving input.

[0085] Any steps, operations or procedures described herein may be performed or implemented using one or more hardware or software modules, either alone or in combination with other devices. In one embodiment, the software module is implemented using a computer program product including a computer-readable medium containing computer program code, which can be executed by a computer processor to perform any or all of the described steps, operations or procedures.

[0086] The foregoing description of the implementation of the present application has been given for the purpose of illustration and description. The foregoing description is not exhaustive nor is it intended to limit the present application to the exact form disclosed, and various variations and modifications may exist according to the above teachings or may be obtained from the practice of the present application. These embodiments are selected and described in order to illustrate the principles of the present application and its practical application, so that those skilled in the art can utilize the present application in various embodiments and various modifications suitable for the specific purpose contemplated.

[0087] Regarding the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0088] It can be further understood that, unless otherwise specified, “connection” includes a direct connection without other components between the two, and also includes an indirect connection with other components between the two.

[0089] It is further understood that, although the operations are described in a specific order in the drawings in the embodiments of the present application, it should not be understood as requiring the operations to be performed in the specific order or serial order shown, or requiring the execution of all the operations shown to obtain the desired results. In certain environments, multitasking and parallel processing may be advantageous.

[0090] Those skilled in the art will readily appreciate other embodiments of the present application after considering the specification and practicing the invention disclosed herein. The present application is intended to encompass any modification, use or adaptation of the present application, which follows the general principles of the present application and includes common knowledge or customary techniques in the field of the present application that are not disclosed in the present application. The specification and embodiments are intended to be exemplary only, and the true scope and spirit of the present application are indicated by the following scope of rights.

[0091] It should be understood that the present application is not limited to the precise structures described above and shown in the drawings, and that various modifications and changes may be made without departing from the scope thereof. The scope of the present application is limited only by the scope of the appended claims.

[0092] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application 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 make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. An optimization design method for three-dimensional slope anti-slide piles, characterized in that: The following steps are involved: Step S1: generating a dynamic library file of the soil constitutive model based on the soil constitutive model generating software, and importing the dynamic library file into the general interactive modeling software; Step S2: establishing a three-dimensional model of the slope and the anti-slide pile based on the user interface of the general interactive modeling software, setting material properties for the three-dimensional model, defining boundary conditions, and dividing the grid; Step S3: Based on the integrated environment of general interactive modeling software and soil constitutive model generation software, static analysis is performed on the three-dimensional model to generate a visualization toolkit file; Step S4: extracting material point information in the visualization toolkit file based on data analysis visualization software, and packaging the material point information into a material point file; Step S5: determine the pre-prepared material point method three-dimensional data analysis software, and judge whether the size of the grid is consistent with the size defined in the material point method three-dimensional data analysis software, if not, return to step S2 to reset, otherwise execute step S6; Step S6: reading the material point file based on the material point method three-dimensional data analysis software and performing dynamic analysis; Step S7: Based on the integrated environment of the general interactive modeling software, the soil constitutive model generation software and the material point method three-dimensional data analysis software, the unstructured grid file of the visualization toolkit file is generated in the dynamic analysis process by using the strength reduction method; Step S8: extracting slope failure state variables from the unstructured grid file of the visualization toolkit file through the data analysis visualization software, for evaluating the disaster reduction effects of different anti-slide pile design schemes; Step S9: taking the slope failure state variable as the optimization target, adjusting the three-dimensional model to obtain the optimal anti-slide pile design solution.

2. The method according to claim 1, characterized in that The adjusting the three-dimensional model comprises: The length and inclination of the anti-slip pile are adjusted.

3. The method according to claim 1, characterized in that The method further comprises: The material point method is used to simulate the large deformation process of the slope; The simulation results of the large deformation process simulation are used together with the slope failure state variables as optimization targets.

4. The method according to claim 1, characterized in that The slope failure state variables include energy, velocity and sliding distance.

5. An optimization design device for three-dimensional slope anti-slide piles, characterized in that: include: A generation unit is used to execute step S1, step S3, and step S7; the step S1: based on the soil constitutive model generation software, a dynamic library file of the soil constitutive model is generated, and the dynamic library file is imported into the general interactive modeling software; the step S3: based on the general interactive modeling software and the soil constitutive model generation software integrated environment, a static analysis is performed on the three-dimensional model to generate a visualization toolkit file; the step S7: based on the general interactive modeling software, the soil constitutive model generation software and the material point method three-dimensional data analysis software integrated environment, a strength reduction method is used to generate an unstructured grid file of the visualization toolkit file during the dynamic analysis process; A processing unit, for executing step S2, step S4, step S6, step S8, and step S9; said step S2: establishing a three-dimensional model of a slope and an anti-slide pile based on the user interface of the general interactive modeling software, setting material properties for the three-dimensional model, defining boundary conditions, and dividing the grid; said step S4: extracting material point information in the visualization toolkit file based on data analysis visualization software, and packaging the material point information into a material point file; said step S6: reading the material point file based on the material point method three-dimensional data analysis software, and performing dynamic analysis; said step S8: extracting slope failure state variables from the unstructured grid file of the visualization toolkit file through the data analysis visualization software, for evaluating the disaster reduction effects of different anti-slide pile design schemes; said step S9: adjusting the three-dimensional model with the slope failure state variables as the optimization target to obtain the optimal anti-slide pile design scheme; The judgment unit is used to execute step S5; the step S5: determines the pre-prepared material point method three-dimensional data analysis software, and judges whether the size of the grid is consistent with the size defined in the material point method three-dimensional data analysis software, if not, returns to step S2 to reset, otherwise executes step S6.

6. The device according to claim 5, characterized in that The processing unit adjusts the three-dimensional model in the following manner: The length and inclination of the anti-slip pile are adjusted.

7. The device according to claim 5, characterized in that The processing unit is also used for: The material point method is used to simulate the large deformation process of the slope; The simulation results of the large deformation process simulation are used together with the slope failure state variables as optimization targets.

8. The device according to claim 5, characterized in that The slope failure state variables include energy, velocity and sliding distance.

9. An electronic device comprising a memory, a processor and a computer program stored in the memory, characterized in that: The processor executes the computer program to implement the method according to any one of claims 1 to 4.

10. A computer-readable storage medium having a computer program / instruction stored thereon, characterized in that: When the computer program / instructions are executed by a processor, the method according to any one of claims 1 to 4 is implemented.

Images