Discrete element method for geological disaster simulation based on contact between spherical particles and triangular meshes
By determining the contact type between spherical particles and triangular mesh in the discrete element method and calculating the contact force in detail, the problem of discontinuity of contact force calculation in the prior art is solved, and the accuracy and stability of geological disaster simulation are improved.
Patent Information
- Application Number
- CN202510037929.X
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2025-01-09
- Publication Date
- 2025-05-06
- Estimated Expiration
- 2045-01-09
AI Technical Summary
The existing discrete element method ignores the topological relationship between triangle meshes in disaster simulation, resulting in discontinuity of contact force calculations, affecting the accuracy of simulation results.
By determining the type of contact between spherical particles and the triangular mesh (surface contact, edge contact and point contact), the contact force is calculated refinely to ensure a smooth transition between adjacent triangular mesh.
It improves the accuracy of geological disaster simulation, ensures the continuity of contact force calculation and energy conservation, and avoids errors and "splash" in calculations.
Smart Images
Figure CN119494244B_ABST
Abstract
Description
Technical Field
[0001] The present disclosure relates to the technical field of disaster simulation, and in particular to a discrete element geological disaster simulation method based on contact between spherical particles and triangular meshes. Background Art
[0002] In the process of disaster simulation using discrete element method, the disaster body is usually represented as a series of closely contacting spherical particles, and the terrain surface is represented as a series of interconnected triangular meshes. On this basis, in order to finely simulate the movement evolution process of the disaster, the contact force calculation between the spherical particles and the triangular meshes is crucial.
[0003] In the process of calculating contact force, related technologies simplify or even ignore the topological relationship between triangular meshes, resulting in the inability to smoothly transition the contact force between adjacent triangles, making the disaster simulation results less accurate or even erroneous. Summary of the invention
[0004] In order to overcome the problems existing in the related art, the present disclosure provides a discrete element geological disaster simulation method based on the contact between spherical particles and triangular meshes.
[0005] According to a first aspect of an embodiment of the present disclosure, a discrete element geological disaster simulation method based on contact between spherical particles and triangular meshes is provided, comprising:
[0006] Constructing a geological disaster model, wherein the geological disaster model includes a terrain surface and a disaster body, wherein at least a portion of the terrain surface and the disaster body are in contact with each other;
[0007] Dividing the terrain surface into a plurality of interconnected triangular grids, and discretizing the disaster body into a plurality of mutually contacting spherical particles;
[0008] For a spherical particle and a triangular mesh that are in contact with each other, determining the contact type between the spherical particle and the triangular mesh, and calculating the contact force between the spherical particle and the triangular mesh according to the contact type; wherein the contact type includes: surface contact, edge contact and point contact;
[0009] According to the contact force between the mutually contacting spherical particles and the triangular grid, the interaction between the disaster body and the terrain surface and the movement evolution process of the disaster body are simulated.
[0010] In some embodiments, determining the contact type between the spherical particle and the triangular mesh comprises:
[0011] For spherical particles and triangular meshes that are in contact with each other, the surface projection point of the center of the spherical particle on the triangular mesh is calculated;
[0012] If the surface projection point is inside the triangular mesh, the contact type between the spherical particle and the triangular mesh is determined to be surface contact.
[0013] In some embodiments, the method further comprises:
[0014] For the spherical particles and triangular meshes that are in contact with each other, if the surface projection point is not inside the triangular mesh, the edge projection point of the sphere center on the straight line where each edge of the triangular mesh is located is calculated respectively;
[0015] If any edge projection point of the spherical center falls on the edge of the triangular mesh, the contact type between the spherical particle and the triangular mesh is determined to be edge contact.
[0016] In some embodiments, the method further comprises:
[0017] For the spherical particles and triangular meshes that are in contact with each other, if the contact type between the spherical particles and the triangular meshes is neither surface contact nor edge contact, the contact type between the spherical particles and the triangular meshes is determined to be point contact.
[0018] In some embodiments, the step of calculating the contact force between the spherical particle and the triangular mesh according to the contact type includes:
[0019] In the case where the contact type between the spherical particle and the triangular mesh is edge contact, respectively obtaining the normal vector of each adjacent triangular mesh of the target edge, wherein the target edge is an edge on which any edge projection point of the sphere center can fall;
[0020] Based on the angle between the normal vector of each adjacent triangular mesh and the contact normal vector of the contact point, the edge distribution coefficient corresponding to each adjacent triangular mesh is calculated respectively;
[0021] Based on the edge partition coefficient corresponding to each adjacent triangular mesh, the contact force between the spherical particle and each adjacent triangular mesh is calculated respectively.
[0022] In some embodiments, the step of calculating the contact force between the spherical particle and the triangular mesh according to the contact type includes:
[0023] When the contact type between the spherical particle and the triangular mesh is point contact, respectively obtain the normal vector of each adjacent triangular mesh of the target point, wherein the target point is the vertex in the triangular mesh that contacts the spherical particle;
[0024] Based on the angle between the normal vector of each adjacent triangular mesh and the contact normal vector of the contact point, the point allocation coefficient corresponding to each adjacent triangular mesh is calculated respectively;
[0025] Based on the point distribution coefficient corresponding to each adjacent triangular mesh, the contact force between the spherical particle and each adjacent triangular mesh is calculated respectively.
[0026] In some embodiments, the method further comprises:
[0027] Based on the hierarchical bounding box method, the contact detection between the triangular mesh and the spherical particles is performed to obtain multiple groups of candidate triangular meshes and spherical particles;
[0028] For each group of candidate triangular meshes and spherical particles, determine whether the distance from the projection point of the center of the spherical particle on the plane where the triangular mesh is located to the center of the spherical particle is less than or equal to the radius of the spherical particle; if so, determine that the spherical particle is in contact with the triangular mesh.
[0029] According to a second aspect of an embodiment of the present disclosure, a discrete element geological disaster simulation device based on contact between spherical particles and triangular meshes is provided, comprising:
[0030] A modeling module, used for constructing a geological disaster model, wherein the geological disaster model includes a terrain surface and a disaster body, wherein at least a portion of the terrain surface and the disaster body are in contact with each other;
[0031] A subdivision module, used for dividing the terrain surface into a plurality of interconnected triangular grids, and discretizing the disaster body into a plurality of mutually contacting spherical particles;
[0032] A contact calculation module is used to determine the contact type between the spherical particles and the triangular mesh that are in contact with each other, and calculate the contact force between the spherical particles and the triangular mesh according to the contact type; wherein the contact type includes: surface contact, edge contact and point contact;
[0033] The simulation module is used to simulate the interaction between the disaster body and the terrain surface and the movement evolution process of the disaster body according to the contact force between the mutually contacting spherical particles and the triangular grid.
[0034] According to a third aspect of an embodiment of the present disclosure, there is provided an electronic device, including:
[0035] Processor; and
[0036] A memory, configured to store executable instructions of the processor;
[0037] The processor is configured to execute the method described in the first aspect by executing the executable instructions.
[0038] According to a fourth aspect of an embodiment of the present disclosure, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the method described in the first aspect is implemented.
[0039] The solution provided by the embodiment of the present disclosure can divide the terrain surface into multiple triangular meshes and discretize the disaster body into multiple spherical particles after constructing a geological disaster model including a terrain surface and a disaster body. Subsequently, according to the contact type between the mutually contacting spherical particles and the triangular meshes, the contact force between the spherical particles and the triangular meshes is refinedly calculated, thereby ensuring the continuity of the contact force calculation between adjacent triangular meshes, thereby improving the accuracy of geological disaster simulation. BRIEF DESCRIPTION OF THE DRAWINGS
[0040] Figure 1 A schematic flow chart of a geological disaster simulation method in an embodiment of the present disclosure is shown.
[0041] Figure 2 A schematic diagram of a geological disaster model in an embodiment of the present disclosure is shown.
[0042] Figure 3 A schematic diagram of a triangular mesh in an embodiment of the present disclosure is shown.
[0043] Figure 4 A schematic diagram showing the projection relationship between the center of a spherical particle and a triangular mesh in an embodiment of the present disclosure.
[0044] Figure 5 A schematic diagram showing the contact relationship between spherical particles and triangular meshes in an embodiment of the present disclosure is shown.
[0045] Figure 6 A schematic diagram showing the application of the geological disaster simulation method in an embodiment of the present disclosure.
[0046] Figure 7 A schematic structural diagram of a geological disaster simulation device in an embodiment of the present disclosure is shown.
[0047] Figure 8 A schematic structural diagram of an electronic device in an embodiment of the present disclosure is shown. DETAILED DESCRIPTION
[0048] Exemplary embodiments will be described in detail herein, examples of which are shown in the accompanying drawings. When the following description refers to the drawings, the same numbers in different drawings represent the same or similar elements unless otherwise indicated. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present disclosure. Instead, they are merely examples of devices and methods consistent with some aspects of the present disclosure as detailed in the appended claims.
[0049] Discrete Element Method (DEM) is a commonly used numerical calculation and analysis method in the fields of granular materials, industrial engineering, agricultural machinery, heavy machinery, rock and soil mechanics and geological disasters. Especially in the analysis of geological disasters such as collapse, landslide and debris flow, DEM numerical calculation and analysis method is often used in the analysis of its disaster dynamic process, disaster mechanism and disaster prevention and mitigation.
[0050] In DEM, spherical particles are the most commonly used particle form in discrete elements due to their simple contact algorithm and high computational efficiency. When using DEM for geological disaster analysis, the disaster body (such as landslide body, debris flow material, etc.) can be discretized into a series of spherical particles.
[0051] Complex terrain can be geometrically viewed as a complex three-dimensional space surface. In order to characterize the terrain, in the field of three-dimensional geometry, it is usually discretized into a series of interconnected (shared edges) triangular meshes. Therefore, the contact calculation between spherical particles and triangular meshes is often encountered in discrete element calculations. The spherical particle-triangular mesh contact geometry algorithm is of great significance to the accuracy and stability of numerical calculations. In the algorithm, it is necessary to consider the continuity of the contact force information when the spherical particle crosses from one triangular mesh to another adjacent triangular mesh. However, the contact algorithm between spherical particles and triangular meshes in the existing algorithms is usually greatly simplified, ignoring the topological relationship between triangular meshes, resulting in large fluctuations and discontinuities in force when the spherical particle enters from one triangle to another adjacent triangle with shared edges, and leads to non-conservation of energy during the calculation process, and even "splashing" of spherical particles due to calculation errors.
[0052] In view of this, the solution provided in the embodiment of the present disclosure divides the contact types between spherical particles and triangular meshes into three types, and calculates the contact force between the spherical particles and the triangular meshes according to the specific contact types, so that the calculation of the contact force can be adjusted as the contact relationship between the spherical particles and the triangular meshes changes, thereby ensuring a smooth transition of the contact force and ensuring energy conservation in the calculation.
[0053] Next, specific embodiments of the present disclosure will be described in detail.
[0054] First, a geological disaster simulation method is provided in an embodiment of the present disclosure, and the method can be executed by any computer.
[0055] Figure 1 A schematic diagram of a geological disaster simulation method in an embodiment of the present disclosure is shown. Figure 1 As shown, the geological disaster simulation method provided in the embodiment of the present disclosure includes the following steps S101 to S104.
[0056] S101, constructing a geological disaster model, wherein the geological disaster model includes a terrain surface and a disaster body, and the terrain surface and the disaster body are in contact with each other in at least a portion of the area.
[0057] It should be noted that the geological disaster model can be a three-dimensional geometric model. When a geological disaster occurs, the rock (soil) body that has not moved is the terrain surface, and the rock (soil) body that has moved is the disaster body.
[0058] Figure 2 The schematic diagram of the geological disaster model in the embodiment of the present disclosure is shown as an example. For the convenience of observation, Figure 2 The geological disaster model is presented in a two-dimensional cross-section. In the geological disaster model, the disaster body is located on the upper part of the terrain surface and contacts the terrain surface, and the disaster body can slide along the terrain surface.
[0059] For example, for the terrain surface, the construction of the terrain surface in the geological hazard model can be completed by obtaining high-precision Digital Elevation Model (DEM) data. These data provide detailed information on the terrain surface and can help analyze the terrain undulation, slope, slope direction, etc.
[0060] For disaster bodies, historical disaster event data can be obtained, such as the spatial distribution, impact range, and occurrence frequency of landslides and mudslides, or disaster bodies can be identified and calibrated through remote sensing images and field surveys.
[0061] S102, dividing the terrain surface into a plurality of interconnected triangular grids, and discretizing the disaster body into a plurality of mutually contacting spherical particles.
[0062] For example, please refer to Figure 2 , among the multiple spherical particles used to represent the disaster body, each spherical particle is in close contact with each other and is stacked on the terrain surface. Figure 3 As shown, among the multiple triangular meshes used to represent the terrain surface, each triangular mesh is interconnected with (shares a common edge with) at least one other triangular mesh other than the triangular mesh.
[0063] S103, for the spherical particles and the triangular mesh that are in contact with each other, determining the contact type between the spherical particles and the triangular mesh, and calculating the contact force between the spherical particles and the triangular mesh according to the contact type.
[0064] Among them, the contact types include: surface contact, edge contact and point contact.
[0065] In the disclosed embodiments, the contact force between the spherical particle and the triangular mesh can be calculated in different ways according to the contact type. When the contact type is edge contact or point contact, the contact force can be dispersed to the triangular meshes adjacent to the contact point during the calculation process, thereby avoiding the calculated contact force from jumping when the spherical particle moves between different triangular meshes.
[0066] S104, simulating the interaction between the disaster body and the terrain surface and the movement evolution process of the disaster body according to the contact force between the mutually contacting spherical particles and the triangular mesh.
[0067] For example, the discrete element method can be used to simulate the interaction between the disaster body and the terrain surface and the movement evolution process of the disaster body. Specifically, after obtaining the contact force between the spherical particles and the triangular grid, the force between the unit and other units in all directions and the external forces caused by other physical fields acting on the unit can be calculated, and the acceleration of the unit can be obtained according to Newton's second law of motion. Then, the time integration is performed to obtain the velocity and displacement of the unit, and then the physical quantities such as velocity, acceleration, angular velocity, linear displacement and angle of rotation of all units at any time are obtained. These physical quantities can be used to simulate the interaction between the disaster body and the terrain surface and the movement evolution process of the disaster body.
[0068] In some embodiments, before executing S103 and S104, contact determination may be performed on the spherical particles and the triangular meshes in advance, so as to screen out the spherical particles and triangular meshes that are in contact with each other.
[0069] For example, firstly, contact detection can be performed on the triangular meshes and the spherical particles based on a bounding volume hierarchy (BVH) method to obtain multiple groups of candidate triangular meshes and spherical particles.
[0070] Among them, the candidate triangular meshes and spherical particles can be understood as triangular meshes and spherical particles with a high possibility of contact. It can be understood that the hierarchical bounding box method... Based on the hierarchical bounding box method, the contact detection process can be accelerated. By checking the overlap of the bounding boxes layer by layer, those triangular meshes that are unlikely to be in contact with the spherical particles can be quickly eliminated, and only some adjacent triangular faces that may be in contact with the spherical particles are retained, and for these triangular faces, fine contact judgments are carried out one by one, thereby improving the efficiency of contact detection.
[0071] Then, for each group of candidate triangular meshes and spherical particles, it can be determined whether the distance from the projection point of the center of the spherical particle on the plane where the triangular mesh is located to the center of the spherical particle is less than or equal to the radius of the spherical particle.
[0072] If so, it is determined that the spherical particle and the triangular mesh are in contact with each other, and the contact type can be further determined.
[0073] If not, it is determined that there is no contact between the spherical particle and the triangular mesh.
[0074] In some embodiments, please refer to Figure 4 , Figure 4 Schematic diagram showing the projection relationship between the center of a spherical particle and a triangular grid in an embodiment of the present disclosure. Figure 4 In the figure, A, B, and C are the three corner points of the triangular mesh. is the normal vector of the triangular mesh △ABC, S is the center of the spherical particle, and P is the projection point of the center S in the triangular mesh.
[0075] The following will be combined Figure 4 , details how to determine the contact type between spherical particles and triangular meshes, and further explains how to calculate the contact force.
[0076] ①First, determine whether the contact type is surface contact.
[0077] For example, for a spherical particle and a triangular mesh that are in contact with each other, the surface projection point of the center of the spherical particle on the triangular mesh may be calculated.
[0078] If the surface projection point is inside the triangular mesh, the contact type between the spherical particle and the triangular mesh is determined to be surface contact.
[0079] For more details, please refer to Figure 4 According to the geometric relationship in the triangle mesh △ABC, the normal vector And any corner point (for example, point A), the position of the projection point P can be obtained by calculation.
[0080] (1)
[0081] in, Reason Figure 4 The vector from the origin of the coordinate system of the triangular mesh shown points to point P, is the vector pointing to point S from the origin of the same coordinate system.
[0082] If point P satisfies the conditions shown in equations 2-a, 2-b and 2-c, point P is determined to be inside the triangular mesh △ABC, and the spherical particle is in surface contact with the triangular mesh (e.g. Figure 5 Otherwise, point P is outside the triangular mesh △ABC, and the contact judgment of edge contact type continues.
[0083] (2-a)
[0084] (2-b)
[0085] (2-c)
[0086] For example, when the contact type is surface contact, the contact force between the spherical particle and the triangular mesh is as follows:
[0087] (3)
[0088] (4)
[0089] in, is the contact normal force, is the contact tangential force, is the contact modulus, is the material Poisson’s ratio, is the radius of the sphere, is the embedding depth of the spherical particle in the normal direction of the triangle, is the tangential displacement increment of the contact point between the spherical particle and the triangular mesh, It is the historical tangential force.
[0090] At the same time, when the spherical particle makes surface contact with the triangular mesh △ABC, the contact direction is the normal direction of the triangular mesh and the contact point is the above projection point .
[0091] ② If the contact type is not surface contact, continue to determine whether the contact type is edge contact.
[0092] Exemplarily, for a spherical particle and a triangular mesh that are in contact with each other, if the surface projection point is not inside the triangular mesh (ie, the contact type is non-surface contact), the edge projection point of the sphere center on the straight line where each edge of the triangular mesh is located is calculated respectively.
[0093] If any edge projection point of the spherical center falls on the edge of the triangular mesh, the contact type between the spherical particle and the triangular mesh is determined to be edge contact.
[0094] Specifically, the following calculations may be performed on the three edges of the triangular mesh respectively.
[0095] Please continue to refer to Figure 4 The geometric relationship in For example, the following formula can be used to calculate the center S of the sphere at the edge Projection point on the straight line :
[0096] (5)
[0097] If point Located in Ling , and there is no surface contact with any triangular mesh, then the contact type between the spherical particle and the triangular mesh is determined to be edge contact (such as Figure 5 shown in green area).
[0098] In this case, the contact point is the projected point , the contact normal vector By point Pointing to the center of the sphere S:
[0099] (6)
[0100] For example, when the contact type is edge contact, the contact force between the spherical particle and the triangular mesh can be calculated as follows:
[0101] The normal vector of each adjacent triangular mesh of the target edge is obtained respectively, and the target edge is the edge on which any edge projection point of the sphere center can fall.
[0102] Based on the angle between the normal vector of each adjacent triangular mesh and the contact normal vector of the contact point, the edge allocation coefficient corresponding to each adjacent triangular mesh is calculated respectively.
[0103] Based on the edge partition coefficient corresponding to each adjacent triangular mesh, the contact force between the spherical particle and each adjacent triangular mesh is calculated respectively.
[0104] Specifically, spherical particles and prismatic particles are still Take edge contact as an example, assuming and Respectively with the edge The surface normal vectors of two adjacent triangle meshes are , , then the spherical particles and the triangular mesh The contact normal force and contact tangential force are as follows:
[0105] (7)
[0106] (8)
[0107] (9)
[0108] (10)
[0109] In the above formulas, is the contact normal force, is the contact tangential force, is the contact modulus, is the material Poisson’s ratio, is the radius of the sphere, is the embedding depth of the spherical particle in the normal direction of the triangle, is the tangential displacement increment of the contact point between the spherical particle and the triangular mesh, is the historical tangential force, is the contact normal vector, is the edge partition coefficient, Used to represent the angle between two vectors.
[0110] It can be seen that Equation (7) and Equation (8) respectively add the edge distribution coefficient based on the above Equation (1) and Equation (2). , so that the contact force will be distributed according to the angle between the contact direction and the normals of the two triangles.
[0111] By setting the distribution coefficient , so that the center of the ball is Contact transition to surface In the process of contact, The proportion of contact force on the ball gradually decreases, and The proportion of contact force with the ball gradually increases to ensure a smooth transition of contact force.
[0112] For example, based on the same principle, we can also get the sphere and triangle mesh The contact force between them is not elaborated in detail in the embodiment of the present disclosure.
[0113] ③ If the contact type is not surface contact and the contact type is not edge contact, then the contact type can be determined to be point contact (such as Figure 5 shown in the middle blue area).
[0114] For example, please refer to Figure 4 , taking the contact between a spherical particle and the corner point A of the triangular mesh △ABC as an example, when this contact occurs, the following equation (11) can be obtained.
[0115] (11)
[0116] The contact point is , the contact normal vector By point Point to the center of the sphere:
[0117] (12)
[0118] For example, similar to edge allocation, when the contact type is point contact, the contact force between the spherical particle and the triangular mesh can be calculated as follows:
[0119] Respectively obtain the normal vector of each adjacent triangular mesh of the target point, wherein the target point is a vertex in the triangular mesh that contacts the spherical particle;
[0120] Based on the angle between the normal vector of each adjacent triangular mesh and the contact normal vector of the contact point, the point allocation coefficient corresponding to each adjacent triangular mesh is calculated respectively;
[0121] Based on the point distribution coefficient corresponding to each adjacent triangular mesh, the contact force between the spherical particle and each adjacent triangular mesh is calculated respectively.
[0122] Specifically, the spherical particles and corner points For example, when point contact occurs, the contact force can be calculated based on the contact force with the corner point. The normal angles of the connected triangle meshes are distributed, and the normal angles of the triangle meshes connected to the point are The triangular meshes are , , …, , and their surface normal vectors are , , …, , then the spherical particles and the triangular mesh The contact normal force and contact tangential force of (k=1, 2, ..., m) are as follows:
[0123] (13)
[0124] (14)
[0125] (15)
[0126] (16)
[0127] In the above formulas, is the contact normal force, is the contact tangential force, is the contact modulus, is the material Poisson’s ratio, is the radius of the sphere, is the embedding depth of the spherical particle in the normal direction of the triangle, is the tangential displacement increment of the contact point between the spherical particle and the triangular mesh, is the historical tangential force, is the contact normal vector, For the The point distribution coefficient of each face, Used to represent the angle between two vectors.
[0128] It can be understood that by calculating in the above manner, it can be ensured that when the spherical particles transition between the triangular meshes adjacent to the contact corner points, the contact force transitions smoothly without abrupt changes.
[0129] For ease of understanding, the following will combine Figure 6 , taking landslide as an example, the application process of geological disaster simulation using the method provided by the embodiment of the present disclosure is described in detail.
[0130] 1. Establish landslide hazard model
[0131] For example, a three-dimensional geological disaster model including the landslide body and the terrain surface can be established based on the geological survey data of the landslide body using three-dimensional modeling software (such as Catia, Gmsh, etc.).
[0132] 2. Import the established 3D geological disaster model into discrete element calculation and analysis software (such as CoSim), and discretize the terrain surface into triangular meshes, and the landslide body into a series of closely contacted spheres. At the same time, check each triangular face of the triangular mesh that constitutes the terrain to ensure that the normals of all triangular faces point to the side of the landslide body. Figure 6 As shown in (a), the landslide body represented by spherical particles and the terrain surface represented by triangulated network generated by a landslide.
[0133] 3. Apply initial loads such as gravity to each unit in the geological hazard model and configure material parameters.
[0134] 4. Based on the method provided in the present disclosure, the contact force between the landslide body and the terrain surface is calculated, and finally the following is obtained: Figure 6 (b) shows the shape of the landslide body after it slides down and accumulates along the terrain surface, thereby providing a quantitative assessment and analysis of the scope of the landslide disaster.
[0135] Based on the same inventive concept, the present disclosure also provides a geological disaster simulation device, such as the following embodiment. Figure 1 The method embodiment shown is similar to that shown in FIG. 1 , so the implementation of the geological disaster simulation device embodiment can refer to the above Figure 1 The implementation of the illustrated method embodiment will not be described in detail in detail.
[0136] Figure 7 A schematic diagram of the structure of a geological disaster simulation device in an embodiment of the present disclosure is shown. The geological disaster simulation device can be applied to an electronic device, and the electronic device includes a plurality of computing units.
[0137] Specifically, Figure 7As shown, the geological disaster simulation device 700 includes: a modeling module 701, a segmentation module 702, a contact calculation module 703, and a simulation module 704.
[0138] Modeling module 701, used to construct a geological disaster model, the geological disaster model includes a terrain surface and a disaster body, and the terrain surface and the disaster body are in contact with each other at least partially;
[0139] A partitioning module 702 is used to divide the terrain surface into a plurality of interconnected triangular grids, and to discretize the disaster body into a plurality of mutually contacting spherical particles;
[0140] The contact calculation module 703 is used to determine the contact type between the spherical particles and the triangular mesh that are in contact with each other, and calculate the contact force between the spherical particles and the triangular mesh according to the contact type; wherein the contact type includes: surface contact, edge contact and point contact;
[0141] The simulation module 704 is used to simulate the interaction between the disaster body and the terrain surface and the movement evolution process of the disaster body according to the contact force between the mutually contacting spherical particles and the triangular mesh.
[0142] In some embodiments, the contact calculation module 703 is specifically used to calculate the surface projection point of the center of the spherical particle on the triangular mesh for the spherical particles and the triangular mesh that are in contact with each other; if the surface projection point is inside the triangular mesh, the contact type between the spherical particle and the triangular mesh is determined to be surface contact.
[0143] In some embodiments, the contact calculation module 703 is also used to, for spherical particles and triangular meshes that are in contact with each other, calculate the edge projection points of the sphere center on the straight line where each edge of the triangular mesh is located if the surface projection point is not inside the triangular mesh; if any edge projection point of the sphere center falls on the edge of the triangular mesh, the contact type between the spherical particle and the triangular mesh is determined to be edge contact.
[0144] In some embodiments, the contact calculation module 703 is further used to determine that the contact type between the spherical particle and the triangular mesh that are in contact with each other is point contact if the contact type between the spherical particle and the triangular mesh is neither surface contact nor edge contact.
[0145] In some embodiments, the contact calculation module 703 is also used to, when the contact type between the spherical particle and the triangular mesh is edge contact, respectively obtain the normal vector of each adjacent triangular mesh of the target edge, the target edge being the edge on which any edge projection point of the sphere center can fall; based on the angle between the normal vector of each adjacent triangular mesh and the contact normal vector of the contact point, respectively calculate the edge distribution coefficient corresponding to each adjacent triangular mesh; based on the edge distribution coefficient corresponding to each adjacent triangular mesh, respectively calculate the contact force between the spherical particle and each adjacent triangular mesh.
[0146] In some embodiments, the contact calculation module 703 is also used to, when the contact type between the spherical particle and the triangular mesh is point contact, respectively obtain the normal vector of each adjacent triangular mesh of the target point, where the target point is the vertex in the triangular mesh that is in contact with the spherical particle; based on the angle between the normal vector of each adjacent triangular mesh and the contact normal vector of the contact point, respectively calculate the point distribution coefficient corresponding to each adjacent triangular mesh; based on the point distribution coefficient corresponding to each adjacent triangular mesh, respectively calculate the contact force between the spherical particle and each adjacent triangular mesh.
[0147] In some embodiments, the contact calculation module 703 is also used to perform contact detection on triangular meshes and spherical particles based on a hierarchical bounding box method to obtain multiple groups of candidate triangular meshes and spherical particles; for each group of candidate triangular meshes and spherical particles, determine whether the distance from the projection point of the center of the spherical particle on the plane where the triangular mesh is located to the center of the spherical particle is less than or equal to the radius of the spherical particle; if so, determine that the spherical particle and the triangular mesh are in contact with each other.
[0148] Refer to the following Figure 8 800 capable of implementing the embodiments of the present disclosure is described. Figure 8 The electronic device 800 shown is merely an example and should not bring any limitation to the functions and scope of application of the embodiments of the present disclosure.
[0149] like Figure 8 As shown, the electronic device 800 is in the form of a general-purpose computing device. The components of the electronic device 800 may include, but are not limited to: at least one processor 810, at least one memory 820, and a bus 830 connecting different system components (including the memory 820 and the processor 810).
[0150] The memory stores program codes, which can be executed by the processor 810, so that the processor 810 executes the steps according to various exemplary embodiments of the present disclosure described in the above “Exemplary Method” section of the present disclosure.
[0151] In some embodiments, the processor 810 may further perform the following steps of the above method embodiment:
[0152] Constructing a geological disaster model, the geological disaster model includes a terrain surface and a disaster body, and the terrain surface and the disaster body are in contact with each other at least in part;
[0153] The terrain surface is divided into multiple interconnected triangular grids, and the hazard body is discretized into multiple spherical particles that are in contact with each other;
[0154] For the spherical particles and triangular meshes that are in contact with each other, the contact type between the spherical particles and the triangular meshes is determined, and the contact force between the spherical particles and the triangular meshes is calculated according to the contact type; wherein the contact type includes: surface contact, edge contact and point contact;
[0155] According to the contact force between the mutually contacting spherical particles and the triangular meshes, the interaction between the disaster body and the terrain surface as well as the movement evolution process of the disaster body are simulated.
[0156] The memory 820 may include a readable medium in the form of a volatile memory, such as a random access memory (RAM) 8201 and / or a cache memory 8202 , and may further include a read-only memory (ROM) 8203 .
[0157] The memory 820 may also include a program / utility 8204 having a set (at least one) of program modules 8205, such program modules 8205 including but not limited to: an operating system, one or more application programs, other program modules, and program data, each of which or some combination may include an implementation of a network environment.
[0158] Bus 830 may be a bus representing one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, a processor, or a local bus using any of a variety of bus architectures.
[0159] The electronic device 800 may also communicate with one or more external devices 840 (e.g., keyboards, pointing devices, Bluetooth devices, etc.), one or more devices that enable a user to interact with the electronic device 800, and / or any device that enables the electronic device 800 to communicate with one or more other computing devices (e.g., routers, modems, etc.). Such communication may be performed via an input / output (I / O) interface 850. Furthermore, the electronic device 800 may also communicate with one or more networks (e.g., a local area network (LAN), a wide area network (WAN), and / or a public network, such as the Internet) via a network adapter 860. Figure 8As shown, the network adapter 860 communicates with other modules of the electronic device 800 via the bus 830. It should be understood that, although not shown in the figure, other hardware and / or software modules can be used in conjunction with the electronic device 800, including but not limited to: microcode, device drivers, redundant processors, external disk drive arrays, RAID systems, tape drives, and data backup storage systems.
[0160] Through the description of the above implementation, it is easy for those skilled in the art to understand that the example implementation described here can be implemented by software, or by software combined with necessary hardware. Therefore, the technical solution according to the implementation of the present disclosure can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (which can be a CD-ROM, a USB flash drive, a mobile hard disk, etc.) or on a network, and includes a number of instructions to enable a computing device (which can be a personal computer, a server, a terminal device, or a network device, etc.) to execute the method according to the implementation of the present disclosure.
[0161] In an exemplary embodiment of the present disclosure, a computer-readable storage medium is also provided, which may be a readable signal medium or a readable storage medium. A program product capable of implementing the above-mentioned method of the present disclosure is stored thereon. In some possible implementations, various aspects of the present disclosure may also be implemented in the form of a program product, which includes a program code, and when the program product is run on a terminal device, the program code is used to cause the terminal device to execute the steps according to various exemplary implementations of the present disclosure described in the above-mentioned "Exemplary Method" section of the present disclosure.
[0162] Those skilled in the art will readily appreciate other embodiments of the present disclosure after considering the specification and practicing the invention disclosed herein. The present disclosure is intended to cover any variations, uses or adaptations of the present disclosure, which follow the general principles of the present disclosure and include common knowledge or customary techniques in the art that are not disclosed in the present disclosure. The description and examples are intended to be exemplary only, and the true scope and spirit of the present disclosure are indicated by the appended claims.
Claims
1. A discrete element geological disaster simulation method based on contact between spherical particles and triangular grids, characterized in that: include: Constructing a geological disaster model, wherein the geological disaster model includes a terrain surface and a disaster body, wherein at least a portion of the terrain surface and the disaster body are in contact with each other; Dividing the terrain surface into a plurality of interconnected triangular grids, and discretizing the disaster body into a plurality of mutually contacting spherical particles; For a spherical particle and a triangular mesh that are in contact with each other, determining the contact type between the spherical particle and the triangular mesh, and calculating the contact force between the spherical particle and the triangular mesh according to the contact type; wherein the contact type includes: surface contact, edge contact and point contact; According to the contact force between the mutually contacting spherical particles and the triangular mesh, the interaction between the disaster body and the terrain surface and the movement evolution process of the disaster body are simulated; The step of calculating the contact force between the spherical particle and the triangular mesh according to the contact type includes: When the contact type between the spherical particle and the triangular mesh is edge contact, respectively obtain the normal vector of each adjacent triangular mesh of the target edge, wherein the target edge is an edge on which any edge projection point of the sphere center can fall; Based on the angle between the normal vector of each adjacent triangular mesh and the contact normal vector of the contact point, the edge distribution coefficient corresponding to each adjacent triangular mesh is calculated respectively; Based on the edge distribution coefficient corresponding to each adjacent triangular mesh, the contact force between the spherical particle and each adjacent triangular mesh is calculated respectively; The calculation method of the edge distribution coefficient corresponding to any adjacent triangular mesh is as follows: , In the formula, is the edge distribution coefficient corresponding to the triangular mesh, is the contact normal vector of the contact point, is the normal vector of the triangle mesh, is the normal vector of another triangle mesh adjacent to the target edge, Used to represent the angle between two vectors.
2. The method according to claim 1, characterized in that The determining of the contact type between the spherical particle and the triangular mesh comprises: For spherical particles and triangular meshes that are in contact with each other, the surface projection point of the center of the spherical particle on the triangular mesh is calculated; If the surface projection point is inside the triangular mesh, the contact type between the spherical particle and the triangular mesh is determined to be surface contact.
3. The method according to claim 2, characterized in that The method further comprises: For the spherical particles and triangular meshes that are in contact with each other, if the surface projection point is not inside the triangular mesh, the edge projection point of the sphere center on the straight line where each edge of the triangular mesh is located is calculated respectively; If any edge projection point of the spherical center falls on the edge of the triangular mesh, the contact type between the spherical particle and the triangular mesh is determined to be edge contact.
4. The method according to claim 3, characterized in that The method further comprises: For the spherical particles and triangular meshes that are in contact with each other, if the contact type between the spherical particles and the triangular meshes is neither surface contact nor edge contact, the contact type between the spherical particles and the triangular meshes is determined to be point contact.
5. The method according to claim 4, characterized in that The step of calculating the contact force between the spherical particle and the triangular mesh according to the contact type further includes: When the contact type between the spherical particle and the triangular mesh is point contact, respectively obtain the normal vector of each adjacent triangular mesh of the target point, wherein the target point is the vertex in the triangular mesh that contacts the spherical particle; Based on the angle between the normal vector of each adjacent triangular mesh and the contact normal vector of the contact point, the point allocation coefficient corresponding to each adjacent triangular mesh is calculated respectively; Based on the point distribution coefficient corresponding to each adjacent triangular mesh, the contact force between the spherical particle and each adjacent triangular mesh is calculated respectively; Among them, the point allocation coefficient corresponding to each adjacent triangle mesh is calculated as follows: , In the formula, For the The point allocation coefficients corresponding to adjacent triangle meshes are: is the contact normal vector of the contact point, For the The normal vectors of adjacent triangle meshes, For the The normal vectors of adjacent triangle meshes, Used to represent the angle between two vectors. is the total number of adjacent triangle meshes.
6. The method according to any one of claims 1 to 5, characterized in that: The method further comprises: Based on the hierarchical bounding box method, the contact detection between the triangular mesh and the spherical particles is performed to obtain multiple groups of candidate triangular meshes and spherical particles; For each group of candidate triangular meshes and spherical particles, determine whether the distance from the projection point of the center of the spherical particle on the plane where the triangular mesh is located to the center of the spherical particle is less than or equal to the radius of the spherical particle; if so, determine that the spherical particle is in contact with the triangular mesh.
7. A discrete element geological disaster simulation device based on contact between spherical particles and triangular grids, characterized in that: include: A modeling module, used for constructing a geological disaster model, wherein the geological disaster model includes a terrain surface and a disaster body, wherein at least a portion of the terrain surface and the disaster body are in contact with each other; A subdivision module, used for dividing the terrain surface into a plurality of interconnected triangular grids, and discretizing the disaster body into a plurality of mutually contacting spherical particles; A contact calculation module is used to determine the contact type between the spherical particles and the triangular mesh that are in contact with each other, and calculate the contact force between the spherical particles and the triangular mesh according to the contact type; wherein the contact type includes: surface contact, edge contact and point contact; A simulation module, used for simulating the interaction between the disaster body and the terrain surface and the movement evolution process of the disaster body according to the contact force between the mutually contacting spherical particles and the triangular mesh; The step of calculating the contact force between the spherical particle and the triangular mesh according to the contact type includes: When the contact type between the spherical particle and the triangular mesh is edge contact, respectively obtain the normal vector of each adjacent triangular mesh of the target edge, wherein the target edge is an edge on which any edge projection point of the sphere center can fall; Based on the angle between the normal vector of each adjacent triangular mesh and the contact normal vector of the contact point, the edge distribution coefficient corresponding to each adjacent triangular mesh is calculated respectively; Based on the edge distribution coefficient corresponding to each adjacent triangular mesh, the contact force between the spherical particle and each adjacent triangular mesh is calculated respectively; The calculation method of the edge distribution coefficient corresponding to any adjacent triangular mesh is as follows: , In the formula, is the edge distribution coefficient corresponding to the triangular mesh, is the contact normal vector of the contact point, is the normal vector of the triangle mesh, is the normal vector of another triangle mesh adjacent to the target edge, Used to represent the angle between two vectors.
8. An electronic device, characterized in that: include: processor; as well as A memory, configured to store executable instructions of the processor; The processor is configured to perform the method of any one of claims 1 to 6 by executing the executable instructions.
9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 6 is implemented.