Mesh division method, device and equipment for earth and rockfill dam model, medium and product

Abstract

The invention relates to the technical field of water conservancy projects, and discloses an earth and rockfill dam model mesh generation method, device and equipment, a medium and a product, and the earth and rockfill dam model mesh generation method comprises the steps of obtaining an initial three-dimensional earth and rockfill dam model, according to the stress structure of the target earth and rockfill dam, performing preliminary grid division on the initial three-dimensional earth and rockfill dam model to obtain a first earth and rockfill dam model; performing mesh refinement processing on the initial cutting belt of the first earth and rockfill dam model to obtain a second earth and rockfill dam model; correcting the overlapping area of the second earth and rockfill dam model to obtain a third earth and rockfill dam model; performing topological geometric adaptation on grids in the third earth and rockfill dam model to obtain a fourth earth and rockfill dam model; the grid physical attributes of the first earth and rockfill dam model are mapped to the fourth earth and rockfill dam model to obtain a fifth earth and rockfill dam model.According to the method, the initial three-dimensional earth and rockfill dam model is subjected to preliminary division, cutting belt refinement, overlapping area correction and other treatment, and the accuracy of grid division is improved.

Description

Mesh generation methods, devices, equipment, media, and products for earth-rock dam models. Technical Field

[0001] This invention relates to the field of water conservancy engineering technology, specifically to a method, apparatus, equipment, medium, and product for mesh generation of earth-rock dam models. Background Technology

[0002] When analyzing earth-rock dams in pumped storage power stations, it is usually necessary to simulate the earth-rock dam. With the development of computational mechanics, finite element analysis (FEA) has been widely used in the design of earth-rock dams and the simulation analysis of other engineering structures.

[0003] However, traditional mesh generation techniques suffer from insufficient computational accuracy and low efficiency when dealing with complex geometries and large-scale 3D models. In particular, in complex structures such as earth-rock dams, traditional mesh generation methods often face the challenge of balancing accuracy and computational efficiency, resulting in inaccurate mesh generation. Summary of the Invention

[0004] This invention provides a method, apparatus, equipment, medium, and product for mesh generation of earth-rock dam models, in order to solve the problem of inaccurate mesh generation caused by traditional mesh generation techniques in related technologies.

[0005] In a first aspect, the present invention provides a mesh generation method for an earth-rock dam model, comprising: obtaining an initial three-dimensional earth-rock dam model corresponding to a target earth-rock dam; performing preliminary mesh generation on the initial three-dimensional earth-rock dam model according to the stress structure of the target earth-rock dam to obtain a first target three-dimensional earth-rock dam model; the initial three-dimensional earth-rock dam model is an earth-rock dam simulation model obtained by three-dimensional modeling of the target earth-rock dam; refining the mesh of the initial cutting zone of the first target three-dimensional earth-rock dam model to obtain a second target three-dimensional earth-rock dam model; the initial cutting zone is a transition mesh zone around the key area in the first target three-dimensional earth-rock dam model; refining the mesh of the second target three-dimensional earth-rock dam model... The overlapping areas of the earth-rock dam model are corrected to obtain the third target three-dimensional earth-rock dam model; the overlapping area is the area where the mesh in the second target three-dimensional earth-rock dam model overlaps with the preset reference mesh; wherein, the preset reference mesh is at least one or a combination of the tool mesh, template mesh, and previous version reference mesh that intersects with the second target three-dimensional earth-rock dam model; the mesh in the third target three-dimensional earth-rock dam model is topologically adapted to obtain the fourth target three-dimensional earth-rock dam model; the physical properties of the mesh of the first target three-dimensional earth-rock dam model are mapped onto the fourth target three-dimensional earth-rock dam model to obtain the fifth target three-dimensional earth-rock dam model.

[0006] The mesh generation method for earth-rock dam models of this invention involves obtaining an initial three-dimensional earth-rock dam model corresponding to the target earth-rock dam, performing preliminary mesh generation on the initial three-dimensional earth-rock dam model based on the stress structure of the target earth-rock dam to obtain a first target three-dimensional earth-rock dam model, and carrying out preliminary mesh layout based on the actual stress structure characteristics of the target earth-rock dam to ensure that the density distribution of the initial mesh is accurately adapted to the stress gradient and stress concentration areas of the dam body. This invention further refines the mesh of the initial cutting zone of the first target three-dimensional earth-rock dam model to obtain a second target three-dimensional earth-rock dam model. For the transition mesh zone around the key areas in the first target three-dimensional earth-rock dam model, refinement is performed to ensure a smooth transition between the meshes of key and non-key areas in the first target three-dimensional earth-rock dam model, making the mesh gradient of the first target three-dimensional earth-rock dam model more reasonable. Finally, this invention corrects the overlapping areas of the second target three-dimensional earth-rock dam model to obtain a third target three-dimensional earth-rock dam model, avoiding calculation conflicts or errors caused by mesh overlap and improving the standardization and compatibility of the mesh. This invention performs topological geometric adaptation on the mesh in the third target 3D earth-rock dam model to obtain the fourth target 3D earth-rock dam model. This adaptation solves the problem of irregular mesh shape, making the mesh geometry and element connection relationships more consistent with computational requirements, thus improving computational convergence and result accuracy. This invention maps the physical properties of the mesh from the first target 3D earth-rock dam model to the fourth target 3D earth-rock dam model to obtain the fifth target 3D earth-rock dam model. Mapping physical properties ensures the correspondence between the mesh and actual engineering properties, making the final mesh not only rationally shaped but also accurately carrying the physical information of the target earth-rock dam, providing a realistic physical basis for subsequent stress and deformation simulations. Compared with related technologies, this invention, through innovative cutting zone concepts, overlap operations, and topological completion steps, achieves complete restoration of the dam material partitioning, filling partitioning, and boundary geometric partitioning of the target earth-rock dam through cutting operations rather than densification. This ensures the accuracy and reliability of numerical calculations, improves the accuracy of mesh generation, and enhances the efficiency of mesh generation.

[0007] In one optional implementation, a preliminary meshing of the three-dimensional earth-rock dam model is performed based on the stress structure of the target earth-rock dam to obtain a first target three-dimensional earth-rock dam model. This includes: dividing the initial three-dimensional earth-rock dam model into multiple engineering regions based on the stress structure of the target earth-rock dam; performing preliminary meshing on each engineering region according to a corresponding preset meshing rule to obtain the first target three-dimensional earth-rock dam model; the preset meshing rule is the basic division standard for different engineering regions.

[0008] In one optional implementation, the initial cutting zone of the first target three-dimensional earth-rock dam model is refined into a mesh to obtain the second target three-dimensional earth-rock dam model. This includes: judging multiple regions in the first target three-dimensional earth-rock dam model according to preset key indicators to obtain multiple key regions; obtaining an initial cutting zone based on the transition mesh zone around each key region; obtaining multiple initial mesh elements that intersect with the initial cutting zone; determining a transition element that matches the cutting surface based on each initial mesh element; converting the transition element into a standard tetrahedral element; and refining the standard tetrahedral element to obtain the second target three-dimensional earth-rock dam model.

[0009] In one optional implementation, the overlapping region of the second target three-dimensional earth-rock dam model is corrected to obtain the third target three-dimensional earth-rock dam model. This includes: extracting the tool mesh intersecting with the second target three-dimensional earth-rock dam model as a preset reference mesh; determining the overlapping region based on the overlapping region between the mesh in the second target three-dimensional earth-rock dam model and the preset reference mesh; removing ineffective cells in the overlapping region using mesh basis functions; and locally refining and merging the mesh in the overlapping region after removing ineffective cells to obtain the third target three-dimensional earth-rock dam model. The local refining and merging are limited to the overlapping region and / or the initial cutting zone, and the non-overlapping regions are not subjected to overall re-division.

[0010] In one optional implementation, topological geometric adaptation is performed on the mesh in the third target three-dimensional earth-rock dam model to obtain the fourth target three-dimensional earth-rock dam model. This includes: completing the topological information of the mesh in the third target three-dimensional earth-rock dam model based on the original node coordinates of each mesh in the third target three-dimensional earth-rock dam model, the geometric displacement correction vector of adjacent meshes, and the topological adjustment matrix; and smoothing the area corresponding to the panel dam and the area corresponding to the rockfill material in the third target three-dimensional earth-rock dam model based on the material stiffness matrix and the interpolation weight coefficient to obtain the fourth target three-dimensional earth-rock dam model.

[0011] In one optional implementation, mapping the mesh physical properties of the first target three-dimensional earth-rock dam model to the fourth target three-dimensional earth-rock dam model to obtain the fifth target three-dimensional earth-rock dam model includes: mapping the mesh physical properties of the first target three-dimensional earth-rock dam model to the fourth target three-dimensional earth-rock dam model based on the mesh physical properties and volume weighting coefficient to obtain the fifth target three-dimensional earth-rock dam model.

[0012] Secondly, the present invention provides a mesh generation device for an earth-rock dam model, comprising: a preliminary mesh generation module, used to acquire an initial three-dimensional earth-rock dam model corresponding to a target earth-rock dam, and to perform preliminary mesh generation on the initial three-dimensional earth-rock dam model according to the stress structure of the target earth-rock dam to obtain a first target three-dimensional earth-rock dam model; the initial three-dimensional earth-rock dam model is an earth-rock dam simulation model obtained by three-dimensional modeling of the target earth-rock dam; a mesh refinement module, used to refine the mesh of the initial cutting zone of the first target three-dimensional earth-rock dam model to obtain a second target three-dimensional earth-rock dam model; the initial cutting zone is a transition mesh zone around the key area in the first target three-dimensional earth-rock dam model; and an overlapping area correction module, used to... The overlapping region of the second target 3D earth-rock dam model is corrected to obtain the third target 3D earth-rock dam model; the overlapping region is the area where the mesh in the second target 3D earth-rock dam model overlaps with the preset reference mesh; wherein, the preset reference mesh is at least one or a combination of the tool mesh, template mesh, and previous version reference mesh that intersects with the second target 3D earth-rock dam model; the topology geometry adaptation module is used to perform topology geometry adaptation on the mesh in the third target 3D earth-rock dam model to obtain the fourth target 3D earth-rock dam model; the physical property mapping module is used to map the physical properties of the mesh of the first target 3D earth-rock dam model to the fourth target 3D earth-rock dam model to obtain the fifth target 3D earth-rock dam model.

[0013] Thirdly, the present invention provides an electronic device, comprising: a memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the computer instructions to perform the mesh generation method for the earth-rock dam model described in the first aspect or any corresponding embodiment thereof.

[0014] Fourthly, the present invention provides a computer-readable storage medium storing computer instructions for causing a computer to execute the mesh generation method for an earth-rock dam model according to the first aspect or any corresponding embodiment thereof.

[0015] Fifthly, the present invention provides a computer program product, including computer instructions for causing a computer to execute the mesh generation method for an earth-rock dam model according to the first aspect or any corresponding embodiment described above. Attached Figure Description

[0016] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0017] Figure 1 is a schematic diagram of an application scenario according to an embodiment of the present invention; Figure 2 is a schematic diagram of the first process of a mesh generation method for an earth-rock dam model according to an embodiment of the present invention; Figure 3 is a schematic diagram of the second process of a mesh generation method for an earth-rock dam model according to an embodiment of the present invention; Figure 4 is a structural block diagram of a mesh generation device for an earth-rock dam model according to an embodiment of the present invention; Figure 5 is a schematic diagram of the hardware structure of an electronic device according to an embodiment of the present invention. Detailed Implementation

[0018] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0019] It is understood that before using the technical solutions disclosed in the various embodiments of the present invention, users should be informed of the types, scope of use, and usage scenarios of the personal information involved in the present invention and their authorization should be obtained in accordance with relevant laws and regulations through appropriate means.

[0020] The terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of this invention, "a plurality of" means two or more, unless otherwise explicitly specified.

[0021] As an optional application scenario of this invention, as shown in Figure 1, the grid division system of the earth-rock dam model may include at least one terminal device and at least one server. Figure 1 exemplarily shows that the system includes a computer 101, a mobile terminal 102 and a server 103, and the terminal devices such as the computer 101 and the mobile terminal 102 are connected to the server 103 through the network 110.

[0022] Specifically, the terminal device can be a smartphone, tablet, laptop, PDA, desktop computer, game console, smart TV, smart wearable device, in-vehicle terminal, VR (Virtual Reality) device, AR (Augmented Reality) device, etc. Server 103 can be a standalone physical server, a server cluster, a distributed system, or a cloud server providing cloud services. Network 110 can be a wired or wireless network, examples of which include, but are not limited to, the Internet, corporate intranet, local area network, wide area network, mobile communication network, and combinations thereof.

[0023] This invention provides a mesh generation method for an earth-rock dam model. By performing preliminary mesh generation, refining cutting zones, and correcting overlapping areas on the initial three-dimensional earth-rock dam model, the accuracy of mesh generation can be improved.

[0024] According to an embodiment of the present invention, a method for mesh generation of an earth-rock dam model is provided. It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions. Furthermore, although a logical order is shown in the flowchart, in some cases, the steps shown or described may be executed in a different order than that shown here.

[0025] This embodiment provides a mesh generation method for an earth-rock dam model, which can be used in computer equipment. Figure 2 is a flowchart of the first type of mesh generation method for an earth-rock dam model according to an embodiment of the present invention. As shown in Figure 2, the process includes the following steps: Step S201, obtain the initial three-dimensional earth-rock dam model corresponding to the target earth-rock dam, and perform preliminary mesh generation on the initial three-dimensional earth-rock dam model according to the stress structure of the target earth-rock dam to obtain the first target three-dimensional earth-rock dam model; the initial three-dimensional earth-rock dam model is an earth-rock dam simulation model obtained by three-dimensional modeling of the target earth-rock dam.

[0026] Among them, the target earth-rock dam is an earth-rock dam that needs to be meshed and analyzed in actual engineering. Earth-rock dam is a common dam structure in water conservancy projects. Its core is a water-retaining structure made of local materials (or mixtures) such as soil, sand, gravel, and pebbles, which are filled by layering and compaction. For example, the target earth-rock dam is an earth-rock dam of a pumped storage power station with a height of 180m. The dam body is composed of core wall, rockfill, face panel and foundation rock mass. The original finite element model contains about 250,000 three-dimensional mesh elements and is divided into 8 material zones and 12 filling zones.

[0027] In some optional implementations, the initial three-dimensional earth-rock dam model is a virtual model obtained by digitally reconstructing the geometric information such as the shape and size of the target earth-rock dam using three-dimensional modeling software, which can be finite element modeling software. The initial three-dimensional earth-rock dam model is preprocessed to partition and simplify the geometry to reduce unnecessary details, thereby improving the efficiency of subsequent mesh generation. The initial three-dimensional earth-rock dam model that is subsequently pre-meshed is the preprocessed initial three-dimensional earth-rock dam model.

[0028] In some alternative implementations, stress structures are used to characterize the stress distribution of the target earth-rock dam under stress. For example, stress structures include dam body structures, core wall or panel structures, rockfill structures, contact interface structures, etc.

[0029] In some alternative implementations, the first target three-dimensional earth-rock dam model is an earth-rock dam model with a basic mesh obtained after preliminary mesh generation.

[0030] In some optional implementations, obtaining an initial three-dimensional earth-rock dam model corresponding to the target earth-rock dam includes: modeling the target earth-rock dam using finite element modeling software to obtain an initial three-dimensional earth-rock dam model.

[0031] In some optional implementations, the initial three-dimensional earth-rock dam model is preliminarily meshed based on the stress structure of the target earth-rock dam, including: analyzing the stress characteristics of the target earth-rock dam, and performing preliminary meshing of the initial three-dimensional earth-rock dam model based on the stress structure obtained from the stress characteristics analysis.

[0032] Step S202: Refine the mesh of the initial cutting zone of the first target three-dimensional earth-rock dam model to obtain the second target three-dimensional earth-rock dam model; the initial cutting zone is the transition mesh zone around the key area in the first target three-dimensional earth-rock dam model.

[0033] The initial cutting zone is the transition zone between the critical area and the non-critical area. For example, the initial cutting zone can be the connection zone between the dam seepage prevention zone (critical area) and the ordinary dam body (non-critical area). The second target three-dimensional earth-rock dam model is a mesh model refined from the initial cutting zone.

[0034] In some optional implementations, the initial cutting zone of the first target three-dimensional earth-rock dam model is refined into a mesh to obtain a second target three-dimensional earth-rock dam model. This includes: using a mesh refinement tool to further subdivide the mesh of the initial cutting zone, reducing the size of the mesh cells, and obtaining the second target three-dimensional earth-rock dam model.

[0035] For example, after refining the mesh of the initial cutting zone of the first target three-dimensional earth-rock dam model, the number of newly added mesh elements in the second target three-dimensional earth-rock dam model is about 30,000, accounting for 12% of the total number of mesh elements in the model.

[0036] Step S203: Correct the overlapping area of ​​the second target three-dimensional earth-rock dam model to obtain the third target three-dimensional earth-rock dam model; the overlapping area is the area where the mesh in the second target three-dimensional earth-rock dam model overlaps with the preset reference mesh.

[0037] Among them, the preset reference grid is a grid template that is set in advance and meets the engineering requirements or calculation standards. The preset reference grid is at least one or a combination of the tool grid, template grid, and previous version reference grid that intersects with the second target three-dimensional earth-rock dam model; the third target three-dimensional earth-rock dam model is a grid model after the overlapping area has been corrected.

[0038] In some optional implementations, the overlapping areas of the second target three-dimensional earth-rock dam model are corrected to obtain the third target three-dimensional earth-rock dam model, including: adjusting the mesh cells of the overlapping areas to obtain the third target three-dimensional earth-rock dam model.

[0039] Step S204: Perform topological geometric adaptation on the mesh in the third target three-dimensional earth-rock dam model to obtain the fourth target three-dimensional earth-rock dam model.

[0040] Among them, topological geometry adaptation is to adjust the topological relationships and geometric shape of the mesh to meet the calculation requirements; the fourth objective, the three-dimensional earth-rock dam model, is the mesh model after topological geometry adaptation.

[0041] In some optional implementations, the mesh in the third target three-dimensional earth-rock dam model is topologically adapted to obtain the fourth target three-dimensional earth-rock dam model, including: checking whether the connections between the meshes in the third target three-dimensional earth-rock dam model are reasonable and whether the shapes are distorted, and adjusting the unreasonable or distorted meshes to obtain the fourth target three-dimensional earth-rock dam model.

[0042] Step S205: Map the mesh physical properties of the first target three-dimensional earth-rock dam model onto the fourth target three-dimensional earth-rock dam model to obtain the fifth target three-dimensional earth-rock dam model.

[0043] Among them, the mesh physical properties are the physical parameters corresponding to the target earth-rock dam material. For example, the mesh physical properties include: the density, elastic modulus, Poisson's ratio, etc. of the dam body soil. The fifth target three-dimensional earth-rock dam model is the final three-dimensional earth-rock dam model after meshing.

[0044] In some alternative implementations, mapping the mesh physical properties of the first target three-dimensional earth-rock dam model to the fourth target three-dimensional earth-rock dam model includes: associating the material physical parameters corresponding to the first target three-dimensional earth-rock dam model with each mesh cell of the fourth target model.

[0045] The mesh generation method for earth-rock dam models provided in this embodiment obtains an initial 3D earth-rock dam model corresponding to the target earth-rock dam. Based on the stress structure of the target earth-rock dam, the initial 3D earth-rock dam model is preliminarily meshed to obtain a first target 3D earth-rock dam model. The mesh is then preliminarily laid out according to the actual stress structure characteristics of the target earth-rock dam, ensuring that the density distribution of the initial mesh accurately matches the stress gradient and stress concentration areas of the dam body. This embodiment further refines the mesh of the initial cutting zone of the first target 3D earth-rock dam model to obtain a second target 3D earth-rock dam model. The transition mesh zone around the key areas in the first target 3D earth-rock dam model is also refined, ensuring a smooth transition between the key and non-key areas, making the mesh gradient of the first target 3D earth-rock dam model more reasonable. Finally, this invention corrects the overlapping areas of the second target 3D earth-rock dam model to obtain a third target 3D earth-rock dam model, avoiding calculation conflicts or errors caused by mesh overlap and improving the standardization and compatibility of the mesh. This invention performs topological geometric adaptation on the mesh in the third target 3D earth-rock dam model to obtain the fourth target 3D earth-rock dam model. This adaptation solves the problem of irregular mesh shape, making the mesh geometry and element connection relationships more consistent with computational requirements, thus improving computational convergence and result accuracy. This invention also maps the physical properties of the mesh from the first target 3D earth-rock dam model to the fourth target 3D earth-rock dam model to obtain the fifth target 3D earth-rock dam model. Mapping physical properties ensures the correspondence between the mesh and actual engineering properties, resulting in a final mesh that is not only morphologically sound but also accurately reflects the physical information of the target earth-rock dam, providing a realistic physical basis for subsequent stress and deformation simulations. Compared with related technologies, this invention, through innovative cutting zone concepts, overlap operations, and topological completion steps, achieves complete restoration of the dam material partitioning, filling partitioning, and boundary geometric partitioning of the target earth-rock dam through cutting operations rather than densification. This ensures the accuracy and reliability of numerical calculations, improves the accuracy of mesh generation, and enhances the efficiency of mesh generation.

[0046] This embodiment provides a mesh generation method for an earth-rock dam model, which can be used in computer equipment. Figure 3 is a second flowchart of the mesh generation method for an earth-rock dam model according to an embodiment of the present invention. As shown in Figure 3, the process includes the following steps: Step S301, obtain the initial three-dimensional earth-rock dam model corresponding to the target earth-rock dam, perform preliminary mesh generation on the initial three-dimensional earth-rock dam model according to the stress structure of the target earth-rock dam, and obtain the first target three-dimensional earth-rock dam model; the initial three-dimensional earth-rock dam model is an earth-rock dam simulation model obtained by three-dimensional modeling of the target earth-rock dam.

[0047] Specifically, step S301 includes: step S3011, dividing the initial three-dimensional earth-rock dam model into multiple engineering areas based on the stress structure of the target earth-rock dam.

[0048] Based on the stress characteristics and material zoning of the target earth-rock dam structure, the initial three-dimensional earth-rock dam model is divided into multiple engineering regions. For example, these multiple engineering regions include: the dam foundation region, the core wall or panel region, the rockfill region, the contact interface region, and other regions. Specifically, the dam foundation region is the area in direct contact with the foundation rock mass, bearing its own weight and overlying load; the core wall or panel region is a key structure for seepage prevention and deformation control, with the core wall being the seepage barrier in the middle of the dam body and the panel being the seepage barrier on the upstream face of the dam; the rockfill region is the area filled with granular materials with larger particle sizes such as boulders and gravel; the contact interface region is the interface region between the core wall and the rockfill, and between the panel and the rockfill; and the other regions are the upper part of the dam body far from the dam toe and interface.

[0049] Step S3012: Perform preliminary meshing for each engineering area according to the corresponding preset meshing rules to obtain the first target three-dimensional earth-rock dam model; the preset meshing rules are the basic division standards for different engineering areas.

[0050] Each engineering area corresponds to a preset grid division rule. The preset grid division rule can be set according to the actual situation. For example, for the dam foundation area, the preset grid division rule can be: the grid size is controlled within 1 / 200 to 1 / 300 of the target earth-rock dam height; for the core wall or panel area, the preset grid division rule can be: if the thickness is less than 1 / 50 of the target earth-rock dam height, then more than 5 layers of grids need to be arranged along the thickness direction; for the rockfill area, the preset grid division rule can be: use a coarser grid, but it needs to be locally refined near the core wall or panel and dam toe area; for the contact interface area, the preset grid division rule can be: if the interface thickness is less than 0.5m, at least 3 layers of grids should be set in the interface normal direction; for other areas, the preset grid division rule can be: a relatively coarse grid of more than 1 / 100 of the dam height can be used.

[0051] Step S302: Refine the mesh of the initial cutting zone of the first target three-dimensional earth-rock dam model to obtain the second target three-dimensional earth-rock dam model; the initial cutting zone is the transition mesh zone around the key area in the first target three-dimensional earth-rock dam model.

[0052] Specifically, step S302 includes: step S3021, judging multiple regions in the first target three-dimensional earth-rock dam model according to preset key indicators to obtain multiple key regions, and obtaining an initial cutting zone based on the transition mesh zone around each key region.

[0053] Among them, the preset key indicators are pre-set indicators used to divide key areas. For example, the preset key indicators may include: geometric feature indicators, stress / strain distribution indicators, seepage characteristic indicators, and engineering experience sensitive indicators.

[0054] In some optional implementations, each preset key indicator corresponds to a key area judgment criterion. For example, for geometric feature indicators, if the local radius of curvature R < 5h (where h is the average grid size), or the polygonal line angle... If the angle is >30°, it is considered a critical region. Regarding stress / strain distribution indices, preliminary static analysis is used; if the quotient of the stress change rate and average stress between adjacent grids is greater than 20%, or the maximum principal strain gradient of the grid is greater than [missing value], then [missing value]. If the seepage gradient is greater than 80% of the critical seepage gradient, it is considered a critical area. For engineering experience-sensitive indicators, the dam toe area (when the dam height is >100m), the contact surface between the core wall and the rockfill, and the upstream and downstream slope toes are all automatically identified as critical areas.

[0055] In some optional implementations, each region in the first target three-dimensional earth-rock dam model is judged based on multiple preset key indicators to obtain multiple key regions.

[0056] In some optional implementations, an initial cutting band is obtained based on the transition mesh band around each critical region, including: generating an initial cutting band with a bandwidth of a preset size along the interface normal around the critical region; wherein the preset size is 1 or 2 mesh sizes.

[0057] Specifically, the size of the mesh cells within the initial cutting zone should be controlled to be 1 / 2 to 1 / 3 of that in the main zone; the aspect ratio of the cells should be less than or equal to 3, and the Jacobian ratio should be greater than or equal to 0.3; the cells outside the initial cutting zone should gradually transition to the main grid, with the cell size growth rate not exceeding 1.5 times to avoid abrupt changes; the total proportion of newly added mesh cells should be less than or equal to 15%, ensuring that the critical response error of the dam body is less than or equal to 5%, while improving the computational efficiency by more than or equal to 30%.

[0058] Step S3022: Obtain multiple initial mesh elements that intersect with the initial cutting strip, and determine the transition element that matches the cutting surface based on each initial mesh element.

[0059] In this process, a transitional cell, either a pyramid or prism, matching the cutting surface is generated in each initial mesh cell to ensure a seamless connection between the cutting interface and the original mesh boundary.

[0060] Step S3023: Convert the transition element into a standard tetrahedral element, refine the standard tetrahedral element into a mesh, and obtain the second target three-dimensional earth-rock dam model.

[0061] Among them, the standard tetrahedral element is the most commonly used basic element in finite element method (a triangular pyramid composed of 4 vertices), which has strong computational stability.

[0062] In some optional implementations, the transition unit is converted into a standard tetrahedral unit based on its cross-sectional shape. For example, if the cross-sectional shape of the transition unit is triangular, the pyramid unit is directly split into a tetrahedron; if the cross-sectional shape of the transition unit is quadrilateral, it is first divided into two triangles and then a tetrahedron is generated; if the cross-sectional shape of the transition unit is a polygon with more than four sides, the polygon is divided into multiple triangles using the Ear-Clipping algorithm or the Delaunay algorithm, and then the corresponding tetrahedron is generated.

[0063] In some alternative implementations, the following constraints must be met during the conversion of transition elements to standard tetrahedral elements: mesh volume conservation error less than 0.5%, element aspect ratio less than or equal to 3, and Jacobian greater than or equal to 0.3.

[0064] In some optional implementations, mesh refinement of the standard tetrahedral elements includes: determining a target size based on the original size of the standard tetrahedral elements, the error estimate of the standard tetrahedral elements, and a target error threshold; and refining the standard tetrahedral elements according to the target size. For example, the formula for determining the target size is:

[0065] in, The target size for a standard tetrahedral element. The original dimensions of a standard tetrahedral unit. This is the error estimate for a standard tetrahedral element. The target error threshold.

[0066] In some alternative implementations, for critical areas (the contact surface between the core wall and the riprap, the dam toe, and the panel interface): the mesh is refined to 1 / 3 to 1 / 2 of the main mesh unit size; for ordinary areas: the original unit mesh size is maintained or only refined by 1 level; the unit mesh growth rate is controlled within the range of less than or equal to 1.5 times to avoid abrupt size changes.

[0067] In some alternative implementations, during the mesh refinement process of standard tetrahedral elements, it is necessary to ensure that the proportion of newly added elements in the initial cutting zone is less than or equal to 15%, and to ensure that the error of the response result in the key area is less than or equal to 5%, thereby improving the overall computational efficiency by more than or equal to 30%.

[0068] Step S303: Correct the overlapping area of ​​the second target three-dimensional earth-rock dam model to obtain the third target three-dimensional earth-rock dam model; the overlapping area is the area where the mesh in the second target three-dimensional earth-rock dam model overlaps with the preset reference mesh.

[0069] Specifically, step S303 includes: step S3031, extracting the tool mesh that intersects with the second target three-dimensional earth-rock dam model as a preset reference mesh, and determining the overlapping area based on the overlapping area between the mesh in the second target three-dimensional earth-rock dam model and the preset reference mesh.

[0070] The formula for determining the degree of overlap in the overlapping regions is:

[0071] in, This represents the mesh in the 3D earth-rock dam model representing the second objective. With preset reference grid The degree of overlap of the overlapping regions is defined by the integration region, which is called the Intersection. A is an area element, representing an infinitesimally small area within the overlapping region.

[0072] Step S3032: Use mesh basis functions to remove ineffective cells in the overlapping region, and then locally refine and merge the mesh in the overlapping region after removing ineffective cells to obtain the third target three-dimensional earth-rock dam model.

[0073] The mesh basis functions include mathematical functions that describe the mesh shape in the second objective 3D earth-rock dam model. Mathematical functions describing the shape of a preset reference grid. The precision of the overlapping region is controlled by the grid basis function; local refinement and merging are limited to the overlapping region and / or the initial cutting zone, and the non-overlapping region is not re-divided as a whole.

[0074] In some alternative implementations, the conservation of geometric and physical quantities is ensured by overlapping integrals, while invalid cells are eliminated. If the overlapping area is less than 5% of the area / volume of the original mesh cell, the mesh cell is marked as invalid and removed to avoid generating malformed meshes.

[0075] In some optional implementations, the overlapping mesh cells must satisfy Jacobian greater than or equal to 0.3, aspect ratio less than or equal to 3, and minimum angle greater than or equal to 10°. For mesh cells that do not meet the conditions, local subdivision or node smoothing is automatically triggered.

[0076] In some optional implementations, the mesh in the overlapping area after removing invalid cells is locally refined and merged, including: refining the mesh cells with an error estimate greater than 0.6 (splitting them into 2 to 4 tetrahedrons); merging the cells with an error estimate less than 0.2 to reduce redundant cells. At the same time, while balancing accuracy and efficiency, the proportion of newly added mesh cells is less than or equal to 15% of the total number of mesh cells.

[0077] In some optional implementations, before locally refining and merging the mesh in the overlapping region after removing invalid cells, the mesh generation method for the earth-rock dam model further includes: geometrically projecting the boundary nodes of the overlapping region to ensure boundary continuity; and using Laplacian smoothing or weighted volume smoothing to make the rate of change of adjacent mesh cell size less than or equal to 1.5.

[0078] In some optional implementations, after locally refining and merging the mesh in the overlapping area after removing non-valid cells, the mesh generation method for the earth-rock dam model further includes: determining the physical quantities of the new mesh cells after local refining and merging based on the physical quantities of the original mesh cells before local refining and merging, the volume of the original mesh cells, and the total volume of the overlapping area, and configuring the physical quantities of the new mesh cells after local refining and merging onto the corresponding new mesh cells.

[0079] For example, the formula for determining the physical quantities of the new mesh cells after local refinement and merging is:

[0080]

[0081] in, The physical quantities of the new mesh cells after local encryption and merging. To determine the physical quantities of the original mesh cells before local densification and merging, For the first The weights of the original grid cells, For the first The volume of each original mesh cell This represents the total volume of the overlapping region.

[0082] In some optional implementations, the meshing method for the earth-rock dam model also includes: automatically detecting abnormal mesh elements with a slenderness ratio less than 0.2 and a Jacobian ratio less than 0.3, and sequentially performing node smoothing, local re-meshing, and template replacement to ensure computational stability.

[0083] Step S304: Perform topological geometric adaptation on the mesh in the third target three-dimensional earth-rock dam model to obtain the fourth target three-dimensional earth-rock dam model.

[0084] Specifically, step S304 includes: step S3041, completing the topological information of the grids in the third target three-dimensional earth-rock dam model based on the original node coordinates of each grid in the third target three-dimensional earth-rock dam model, the geometric displacement correction vector of adjacent grids, and the topological adjustment matrix.

[0085] Specifically, based on the original node coordinates of each mesh in the third-target 3D earth-rock dam model, the geometric displacement correction vectors of adjacent meshes, and the topology adjustment matrix, the corrected coordinates are determined. The topological information of the meshes in the third-target 3D earth-rock dam model is then completed based on these corrected coordinates. For example, the formula for determining the corrected coordinates is:

[0086] in, To correct the coordinates, The original node coordinates of each mesh in the 3D earth-rock dam model for the third objective are given. In the third objective of the three-dimensional earth-rock dam model, the first The geometric displacement correction vector of adjacent grids of each grid. In the third objective of the three-dimensional earth-rock dam model, the first A topology adjustment matrix for each grid, used to record the adjacency relationships of grids, faces, and edges.

[0087] In some optional implementations, the topological information of the mesh in the three-dimensional earth-rock dam model of the third target is completed, which also includes: updating the adjacency list (the adjacency list is used to record the adjacent edges, faces and other meshes of each mesh) through the topological adjustment matrix, and filling in the connection relationships of the cells, edges and faces that were broken after the overlapping area was processed (for example, a face that was originally only connected to 1 cell is connected to 2 cells after completion, ensuring that the face is a common face of two cells).

[0088] In some alternative implementations, the compliance of the topology is verified using Euler's formula for a 3D mesh. Euler's formula is:

[0089] in, The number of vertices. Let be the number of sides. For the number of faces, The number of body units.

[0090] In some alternative implementations, incremental repair is triggered for meshes that do not satisfy Euler's formula: point patching, edge addition, face creation, and volume element closure.

[0091] In some optional implementations, edge completion: when new intersections are generated in overlapping regions, if the original topology lacks corresponding edges, the algorithm automatically generates and adds them to the adjacency list; face completion: if new polygonal faces are generated after cell cutting, ear-clipping triangulation is performed first, and then the cell topology structure is added; volume cell repair: when non-closed volumes appear, they are reconstructed through local Boolean operations to make the cell boundaries form closed polyhedra; redundant information cleanup: isolated nodes and dangling edges are deleted to avoid disrupting the overall coherence.

[0092] Step S3042: Based on the material stiffness matrix and interpolation weight coefficients, smooth the area corresponding to the panel dam and the area corresponding to the rockfill in the third target three-dimensional earth-rock dam model to obtain the fourth target three-dimensional earth-rock dam model.

[0093] In this design, the panel of the concrete dam is a rigid, thin structure, while the rockfill area is a flexible, granular fill area. Using a weighted interpolation formula, the equivalent material parameter matrix of the corresponding areas of the concrete dam and the rockfill area is calculated, allowing for a smooth transition from the rigidity of the panel to the flexibility of the rockfill. For example, the weighted interpolation formula is as follows:

[0094] in, This represents the equivalent material parameter matrix for the region corresponding to the panel dam and the region corresponding to the rockfill. This refers to the number of layers (usually 3 to 5 layers). For the first interpolation weight coefficients of the layer, For the first Material stiffness matrix of the layer.

[0095] In some alternative implementations, to ensure a smooth and continuous transition of stiffness from the panel to the riprap, The value can be:

[0096] in, For the first interpolation weight coefficients of the layer, It is an exponential function. This is the attenuation coefficient, used to control the rate of weight attenuation, and its value is between 1 and 2. This refers to the layer where the current transition unit resides.

[0097] In some optional implementations, the mesh elements between the panel dam area and the rockfill area in the three-dimensional earth-rock dam model of the third objective are transition elements. Through the transition elements, the stiffness and displacement transfer between the panel dam and the rockfill or bedrock are made continuous, avoiding numerical singularities caused by abrupt stiffness changes. Normal stiffness, tangential stiffness and friction coefficient are introduced at the interface, and these contact parameters are inherited through the transition elements to ensure consistent contact surfaces. The thickness of the transition elements is generally 0.5 meters to 1.0 meters, which can effectively reduce the risk of interface cracking or numerical oscillation.

[0098] In this embodiment of the invention, the transition unit retains the high stiffness of the panel while avoiding local instability caused by excessive differences in properties between it and the riprap. The stress distribution error of the mesh element within the transition layer is less than or equal to 5%, the displacement continuity error is less than or equal to 3%, and the overall computational efficiency is improved by approximately 25%.

[0099] Step S305: Map the mesh physical properties of the first target three-dimensional earth-rock dam model onto the fourth target three-dimensional earth-rock dam model to obtain the fifth target three-dimensional earth-rock dam model.

[0100] Specifically, step S305 includes: step S3051, mapping the mesh physical properties onto the fourth target three-dimensional earth-rock dam model based on the mesh physical properties and volume weighting coefficient of the first target three-dimensional earth-rock dam model to obtain the fifth target three-dimensional earth-rock dam model.

[0101] In the process of mesh cutting, in order to ensure the physical consistency of the computational model, the physical properties of the original elements (such as material parameters, boundary conditions, load information, etc.) must be mapped and inherited to the new mesh elements, so as to ensure that the cut mesh and the original mesh are consistent in physical properties.

[0102] In some alternative implementations, the mesh physical properties are mapped onto the fourth objective 3D earth-rock dam model as follows: For material parameters, a volume-weighted approach is used, allocating parameters based on the proportion of overlapping volumes between the new and original elements. For boundary conditions, if the cutting involves boundary elements, area weighting is used for interpolation to ensure the continuity of the boundary conditions. For construction zone / load labels, the maximum coverage principle is adopted, i.e., the original element zone information with the largest proportion in the new element is selected.

[0103] In some optional implementations, the physical properties of the mesh in the first target three-dimensional earth-rock dam model are mapped to the physical properties of the mesh in the fourth target three-dimensional earth-rock dam model based on the physical properties of the mesh in the first target three-dimensional earth-rock dam model and the volume weighting coefficient. This includes: determining the physical properties inherited by the new mesh units in the fourth target three-dimensional earth-rock dam model based on the physical properties of the mesh in the first target three-dimensional earth-rock dam model and the volume weighting coefficient; configuring the physical properties inherited by the new mesh units in the fourth target three-dimensional earth-rock dam model into the corresponding new mesh units; and obtaining the fifth target three-dimensional earth-rock dam model.

[0104] For example, the formula for determining the physical properties inherited by new mesh elements in the fourth objective 3D earth-rock dam model is as follows:

[0105] in, The physical properties inherited by new mesh elements in the fourth objective 3D earth-rock dam model. The first objective is the three-dimensional earth-rock dam model. Volume weighting coefficients for each grid cell The first objective is the three-dimensional earth-rock dam model. Physical properties of each grid cell.

[0106] In some alternative implementations, the middle three-dimensional earth-rock dam model of the first target is determined. The formula for the volume weighting coefficient of each grid cell is:

[0107]

[0108] in, The first objective is the three-dimensional earth-rock dam model. The volume weighting coefficient of each grid cell is used to represent the volume overlap ratio. For the volume of the new mesh element in the fourth objective 3D earth-rock dam model, The first objective is the three-dimensional earth-rock dam model. The volume of each grid cell, when boundary conditions are involved. Replace with area, the rest of the format remains the same.

[0109] In some optional implementations, the meshing method for the earth-rock dam model further includes automatically detecting abnormal or malformed mesh elements after the initial cutting zone is generated and the mesh is generated. The judgment rules are as follows: if the aspect ratio is greater than 5, it is judged as an abnormal mesh element; if the minimum Jacobian determinant is less than 0.3, it is judged as an abnormal mesh element; if the minimum interior angle is less than 5°, it is judged as an abnormal mesh element; if the volume of the new element mesh deviates from the theoretical value by more than 5%, it is judged as an abnormal mesh element; a mesh element that meets any of the conditions is considered an abnormal mesh element.

[0110] In some optional implementations, abnormal mesh elements are corrected using the following formula:

[0111] in, Corrected coordinates for abnormal mesh cells. These are the original node coordinates of the constant mesh element. This is a correction factor, with a value range of constants, and can be 0.5. This is a correction vector calculated based on geometric boundaries or adjacent elements.

[0112] If the quality indicators are still not met after correction, the following processing sequence is executed: node smoothing (Laplacian / volume weighting); local re-meshment (decompose the abnormal elements into multiple tetrahedrons / hexahedrons); template replacement (replace the abnormal mesh elements with standard elements from the standard element library).

[0113] In some alternative implementations, during the mesh cell cutting process, in order to maintain the geometric and physical properties of the boundary, it is necessary to correct the boundary nodes after cutting. The corrected boundary node coordinates are:

[0114] in, These are the corrected boundary node coordinates. These are the original boundary node coordinates before cutting. This is the repair coefficient, used to control the intensity of the repair. Its value ranges from 0 to 1, and is typically between 0.5 and 0.8. The boundary correction vector represents the geometric offset between the cut node and the original boundary, determined by the shortest distance vector from the cut boundary node to the original boundary surface / curve.

[0115] In some optional implementations, after correcting the cut boundary nodes, the geometric offset of the nodes should satisfy:

[0116] in, These are the corrected boundary node coordinates. These are the original boundary node coordinates before cutting. For local feature dimensions, the boundary condition propagation error is kept within the range of less than 1%.

[0117] In this embodiment of the invention, the cut boundary nodes are corrected to ensure that the cut boundary points can be restored to their original geometric boundaries (such as the dam surface and dam foundation outline). This ensures that boundary conditions (such as water pressure, contact conditions, and support constraints) can still be accurately transmitted after repair. This is achieved by introducing a repair coefficient. This avoids numerical oscillations that may result from direct forced projection.

[0118] In some alternative implementations, for large-scale earth-rock dam models with tens of thousands of elements, a recursive densification algorithm is used for mesh discretization, and the degree of mesh densification is controlled by the following formula:

[0119] in, For the process After each iteration, the overall mesh refinement level is determined; the higher the value, the finer the mesh. The number of iterations for encryption. For the first The grid encryption factor at the next iteration To decay the weights and ensure that the encryption level decreases with each iteration, we avoid over-refinement.

[0120] In this embodiment of the invention, the mesh element size is gradually reduced through a progressively decreasing densification process, while keeping the total number of elements within a computable range. The key response results (displacement, stress, seepage field gradient) after iteration deviate from the reference solution by less than or equal to 5%, and the total number of elements increases by no more than 20% of the original model, ensuring engineering feasibility.

[0121] In some alternative implementations, the meshing method for the earth-rock dam model further includes: automatically adjusting the mesh density based on changes in the displacement field, using the following formula:

[0122] in, The adjusted grid density distribution, The initial density distribution of the base grid, For displacement field, The Laplace operator for the displacement field reflects local curvature and strain changes. This is an adaptive factor used to control the adjustment range, typically ranging from 0.3 to 0.7.

[0123] In some alternative implementations, by Local error indices are calculated; regions with large errors are automatically refined, while regions with small errors are appropriately merged. To avoid abrupt changes, the rate of change of mesh size is controlled to be less than or equal to 1.5 times. After each round of calculation, the mesh is re-evaluated and updated based on the new displacement field until the error converges.

[0124] In this embodiment of the invention, the grid cells of the key area (dam toe, core wall interface) are refined to 1 / 2 to 1 / 3 of the basic grid; the overall number of cells increases by no more than 20%; the error of the calculation result relative to the whole-domain refined grid is less than or equal to 5%, and the efficiency is improved by more than or equal to 30%.

[0125] The mesh generation method for earth-rock dam models in this invention is highly specific to dam construction: unlike general algorithms, it is designed for the complex partitioning and contact characteristics of earth-rock dams. It is highly efficient: no overall re-meshing is required; local cutting is sufficient for model updates, significantly improving computational efficiency. It ensures accurate partitioning: after cutting, material partitions, filling partitions, and geometric boundary partitions are automatically restored, avoiding manual operation. It provides numerical reliability: integrated topological and physical completion ensures the geometric and physical consistency of the model after cutting. It exhibits good boundary adaptability: the boundary restoration method guarantees accuracy and stability under high water levels and seismic conditions.

[0126] The mesh generation method for earth-rock dam models in this invention, through innovative concepts of cutting zones, overlap operations, topology completion, and tetrahedral mesh filling, overcomes the bottlenecks of traditional mesh generation techniques in complex earth-rock dam models, improving computational efficiency and ensuring computational accuracy. This technology not only meets the computational needs of complex projects such as earth-rock dams but also has broad application prospects. By using cutting operations rather than mesh refinement, it achieves complete restoration of dam material partitions, filling partitions, and boundary geometric partitions. Furthermore, through topology, physical completion, and dam-specific boundary correction mechanisms, it ensures the accuracy and reliability of numerical calculations. This method is particularly suitable for pumped storage power stations and large earth-rock dam projects, providing an efficient and reliable solution for complex dam modeling and operational safety evaluation.

[0127] This embodiment also provides a mesh generation device for an earth-rock dam model, which is used to implement the above embodiments and preferred embodiments; details already described will not be repeated. As used below, the term "module" can refer to a combination of software and / or hardware that performs a predetermined function. Although the device described in the following embodiments is preferably implemented in software, hardware implementation, or a combination of software and hardware, is also possible and contemplated.

[0128] This embodiment provides a mesh generation device for an earth-rock dam model, as shown in Figure 4, including: a preliminary mesh generation module 401, used to obtain an initial three-dimensional earth-rock dam model corresponding to the target earth-rock dam, and to perform preliminary mesh generation on the initial three-dimensional earth-rock dam model according to the stress structure of the target earth-rock dam to obtain a first target three-dimensional earth-rock dam model; the initial three-dimensional earth-rock dam model is an earth-rock dam simulation model obtained by three-dimensional modeling of the target earth-rock dam.

[0129] The mesh refinement module 402 is used to refine the mesh of the initial cutting zone of the first target three-dimensional earth-rock dam model to obtain the second target three-dimensional earth-rock dam model; the initial cutting zone is the transition mesh zone around the key area in the first target three-dimensional earth-rock dam model.

[0130] The overlapping region correction module 403 is used to correct the overlapping region of the second target three-dimensional earth-rock dam model to obtain the third target three-dimensional earth-rock dam model; the overlapping region is the area where the mesh in the second target three-dimensional earth-rock dam model overlaps with the preset reference mesh; wherein, the preset reference mesh is at least one or a combination of the tool mesh, template mesh, and previous version reference mesh that intersects with the second target three-dimensional earth-rock dam model.

[0131] The topology geometry adaptation module 404 is used to perform topology geometry adaptation on the mesh in the third target three-dimensional earth-rock dam model to obtain the fourth target three-dimensional earth-rock dam model.

[0132] The physical property mapping module 405 is used to map the mesh physical properties of the first target three-dimensional earth-rock dam model to the fourth target three-dimensional earth-rock dam model to obtain the fifth target three-dimensional earth-rock dam model.

[0133] In some optional implementations, the preliminary mesh generation module 401 includes: an engineering area division unit, used to divide the initial three-dimensional earth-rock dam model according to the stress structure of the target earth-rock dam to obtain multiple engineering areas.

[0134] The preliminary mesh division unit is used to perform preliminary mesh division on each engineering area according to the corresponding preset mesh division rules to obtain the first target three-dimensional earth-rock dam model; the preset mesh division rules are the basic division standards for different engineering areas.

[0135] In some optional implementations, the mesh refinement processing module 402 includes: a region judgment unit, used to judge multiple regions in the first target three-dimensional earth-rock dam model according to preset key indicators, to obtain multiple key regions, and to obtain an initial cutting zone based on the transition mesh zone around each key region.

[0136] The transition element determination element is used to obtain multiple initial mesh elements that intersect with the initial cutting strip, and to determine the transition element that matches the cutting surface based on each initial mesh element.

[0137] The mesh refinement unit is used to convert the transition unit into a standard tetrahedral unit, and then refine the standard tetrahedral unit to obtain the second target three-dimensional earth-rock dam model.

[0138] In some optional implementations, the overlapping area correction module 403 includes: an overlapping area determination unit, used to extract the tool mesh that intersects with the second target three-dimensional earth-rock dam model as a preset reference mesh, and to determine the overlapping area based on the overlapping area between the mesh in the second target three-dimensional earth-rock dam model and the preset reference mesh.

[0139] Local processing units are used to remove ineffective units in overlapping regions using mesh basis functions, and to locally refine and merge the mesh in the overlapping regions after removing ineffective units, to obtain the third target three-dimensional earth-rock dam model; wherein, local refinement and merging are limited to the overlapping regions and / or the initial cutting zone, and the non-overlapping regions are not re-divided as a whole.

[0140] In some optional implementations, the topology geometry adaptation module 404 includes: an information completion unit, used to complete the topology information of the grids in the third target three-dimensional earth-rock dam model based on the original node coordinates of each grid in the third target three-dimensional earth-rock dam model, the geometric displacement correction vector of adjacent grids, and the topology adjustment matrix.

[0141] The smoothing unit is used to smooth the area corresponding to the panel dam and the area corresponding to the rockfill in the third objective three-dimensional earth-rock dam model according to the material stiffness matrix and interpolation weight coefficient, so as to obtain the fourth objective three-dimensional earth-rock dam model.

[0142] In some optional implementations, the physical property mapping module 405 includes: a property mapping unit, used to map the mesh physical properties to the fourth target three-dimensional earth-rock dam model according to the mesh physical properties and volume weighting coefficient of the first target three-dimensional earth-rock dam model, so as to obtain the fifth target three-dimensional earth-rock dam model.

[0143] The mesh generation device for earth-rock dam models provided in this embodiment of the invention can execute the mesh generation method for earth-rock dam models provided in any embodiment of the invention, and has the corresponding functional modules and beneficial effects of the method. Further functional descriptions of the above modules and units are the same as in the corresponding embodiments described above, and will not be repeated here.

[0144] Figure 5 is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention.

[0145] Referring specifically to Figure 5, a schematic diagram of a suitable electronic device for implementing embodiments of the present invention is shown below. The electronic device may include a processor (e.g., a central processing unit, a graphics processing unit, etc.) 501, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 502 or a program loaded from memory 508 into random access memory (RAM) 503. The RAM 503 also stores various programs and data required for the operation of the electronic device. The processor 501, ROM 502, and RAM 503 are interconnected via a bus 504. An input / output (I / O) interface 505 is also connected to the bus 504.

[0146] Typically, the following devices can be connected to I / O interface 505: input devices 506 including, for example, touchscreens, touchpads, keyboards, mice, cameras, microphones, accelerometers, gyroscopes, etc.; output devices 507 including, for example, liquid crystal displays (LCDs), speakers, vibrators, etc.; memory devices 508 including, for example, magnetic tapes, hard disks, etc.; and communication devices 509. Communication device 509 allows electronic devices to communicate wirelessly or wiredly with other devices to exchange data. Although Figure 5 shows an electronic device with various devices, it should be understood that it is not required to implement or have all the devices shown, and more or fewer devices may be implemented or have alternatively.

[0147] In particular, according to embodiments of the present invention, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of the present invention include a computer program product comprising a computer program carried on a non-transitory computer-readable medium, the computer program containing program code for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via a communication device 509, or installed from a memory 508, or installed from a ROM 502. When the computer program is executed by the processor 501, it performs the functions defined in the mesh generation method for the earth-rock dam model according to embodiments of the present invention.

[0148] The electronic device shown in Figure 5 is merely an example and should not impose any limitations on the functionality and scope of use of the embodiments of the present invention.

[0149] This invention also provides a computer-readable storage medium. The methods described above according to embodiments of the invention can be implemented in hardware or firmware, or implemented as computer code that can be recorded on a storage medium, or implemented as computer code downloaded via a network and originally stored on a remote storage medium or a non-transitory machine-readable storage medium and then stored on a local storage medium. Thus, the methods described herein can be processed by software stored on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. The storage medium can be a magnetic disk, optical disk, read-only memory, random access memory, flash memory, hard disk, or solid-state drive, etc.; further, the storage medium can also include combinations of the above types of memory. It is understood that computers, processors, microprocessor controllers, or programmable hardware include storage components capable of storing or receiving software or computer code. When the software or computer code is accessed and executed by the computer, processor, or hardware, the mesh generation method for the earth-rock dam model shown in the above embodiments is implemented.

[0150] A portion of this invention can be applied as a computer program product, such as computer program instructions, which, when executed by a computer, can invoke or provide the methods and / or technical solutions according to the invention through the operation of the computer. Those skilled in the art will understand that the forms in which computer program instructions exist in a computer-readable medium include, but are not limited to, source files, executable files, installation package files, etc. Correspondingly, the ways in which computer program instructions are executed by a computer include, but are not limited to: the computer directly executing the instructions, or the computer compiling the instructions and then executing the corresponding compiled program, or the computer reading and executing the instructions, or the computer reading and installing the instructions and then executing the corresponding installed program. Here, the computer-readable medium can be any available computer-readable storage medium or communication medium accessible to a computer.

[0151] Although embodiments of the invention have been described in conjunction with the accompanying drawings, those skilled in the art can make various modifications and variations without departing from the spirit and scope of the invention, and such modifications and variations all fall within the scope defined by the appended claims.

Claims

1. A mesh generation method for an earth-rock dam model, characterized in that, The method includes: obtaining an initial three-dimensional earth-rock dam model corresponding to the target earth-rock dam; performing preliminary meshing on the initial three-dimensional earth-rock dam model according to the stress structure of the target earth-rock dam to obtain a first target three-dimensional earth-rock dam model; the initial three-dimensional earth-rock dam model is an earth-rock dam simulation model obtained by three-dimensional modeling of the target earth-rock dam; refining the mesh of the initial cutting zone of the first target three-dimensional earth-rock dam model to obtain a second target three-dimensional earth-rock dam model; the initial cutting zone is a transition mesh zone around the key area in the first target three-dimensional earth-rock dam model; and refining the mesh of the second target three-dimensional earth-rock dam model. The overlapping region is corrected to obtain a third target 3D earth-rock dam model; the overlapping region is the area where the mesh in the second target 3D earth-rock dam model overlaps with a preset reference mesh; wherein, the preset reference mesh is at least one or a combination of a tool mesh, a template mesh, and a previous version reference mesh that intersects with the second target 3D earth-rock dam model; the mesh in the third target 3D earth-rock dam model is topologically adapted to obtain a fourth target 3D earth-rock dam model; the physical properties of the mesh in the first target 3D earth-rock dam model are mapped onto the fourth target 3D earth-rock dam model to obtain a fifth target 3D earth-rock dam model.

2. The method according to claim 1, characterized in that, The step of performing preliminary meshing of the three-dimensional earth-rock dam model based on the stress structure of the target earth-rock dam to obtain the first target three-dimensional earth-rock dam model includes: dividing the initial three-dimensional earth-rock dam model according to the stress structure of the target earth-rock dam to obtain multiple engineering regions; performing preliminary meshing of each engineering region according to a corresponding preset meshing rule to obtain the first target three-dimensional earth-rock dam model; the preset meshing rule is the basic division standard for different engineering regions.

3. The method according to claim 1 or 2, characterized in that, The step of refining the initial cutting zone of the first target three-dimensional earth-rock dam model to obtain the second target three-dimensional earth-rock dam model includes: judging multiple regions in the first target three-dimensional earth-rock dam model according to preset key indicators to obtain multiple key regions; obtaining an initial cutting zone based on the transition mesh zone around each key region; obtaining multiple initial mesh elements that intersect with the initial cutting zone; determining a transition element that matches the cutting surface based on each initial mesh element; converting the transition element into a standard tetrahedral element; and refining the standard tetrahedral element to obtain the second target three-dimensional earth-rock dam model.

4. The method according to claim 1 or 2, characterized in that, The step of correcting the overlapping region of the second target 3D earth-rock dam model to obtain the third target 3D earth-rock dam model includes: extracting the tool mesh intersecting with the second target 3D earth-rock dam model as a preset reference mesh; determining the overlapping region based on the overlapping region between the mesh in the second target 3D earth-rock dam model and the preset reference mesh; removing ineffective cells in the overlapping region using mesh basis functions; and locally densifying and merging the mesh in the overlapping region after removing the ineffective cells to obtain the third target 3D earth-rock dam model; wherein the local densification and merging are limited to the overlapping region and / or the initial cutting zone, and the non-overlapping regions are not re-divided as a whole.

5. The method according to claim 1 or 2, characterized in that, The step of performing topological geometric adaptation on the mesh in the third target three-dimensional earth-rock dam model to obtain the fourth target three-dimensional earth-rock dam model includes: completing the topological information of the mesh in the third target three-dimensional earth-rock dam model based on the original node coordinates of each mesh, the geometric displacement correction vector of adjacent meshes, and the topological adjustment matrix; and smoothing the area corresponding to the panel dam and the area corresponding to the rockfill material in the third target three-dimensional earth-rock dam model based on the material stiffness matrix and the interpolation weight coefficient to obtain the fourth target three-dimensional earth-rock dam model.

6. The method according to claim 1 or 2, characterized in that, Mapping the mesh physical properties of the first target three-dimensional earth-rock dam model onto the fourth target three-dimensional earth-rock dam model to obtain the fifth target three-dimensional earth-rock dam model includes: mapping the mesh physical properties and volume weighting coefficient of the first target three-dimensional earth-rock dam model onto the fourth target three-dimensional earth-rock dam model to obtain the fifth target three-dimensional earth-rock dam model.

7. A mesh generation device for an earth-rock dam model, characterized in that, The device includes: a preliminary mesh generation module, used to acquire an initial three-dimensional earth-rock dam model corresponding to the target earth-rock dam, and to perform preliminary mesh generation on the initial three-dimensional earth-rock dam model according to the stress structure of the target earth-rock dam to obtain a first target three-dimensional earth-rock dam model; the initial three-dimensional earth-rock dam model is an earth-rock dam simulation model obtained by three-dimensional modeling of the target earth-rock dam; a mesh refinement module, used to refine the mesh of the initial cutting zone of the first target three-dimensional earth-rock dam model to obtain a second target three-dimensional earth-rock dam model; the initial cutting zone is a transition mesh zone around the key area in the first target three-dimensional earth-rock dam model; and an overlapping region correction module, used to refine the mesh of the second target three-dimensional earth-rock dam model. The overlapping areas of the dam model are corrected to obtain a third target 3D earth-rock dam model; the overlapping area is the area where the mesh in the second target 3D earth-rock dam model overlaps with a preset reference mesh; wherein, the preset reference mesh is at least one or a combination of a tool mesh, a template mesh, and a previous version reference mesh that intersects with the second target 3D earth-rock dam model; a topology geometry adaptation module is used to perform topology geometry adaptation on the mesh in the third target 3D earth-rock dam model to obtain a fourth target 3D earth-rock dam model; a physical property mapping module is used to map the physical properties of the mesh of the first target 3D earth-rock dam model to the fourth target 3D earth-rock dam model to obtain a fifth target 3D earth-rock dam model.

8. An electronic device, characterized in that, include: A memory and a processor are interconnected, the memory stores computer instructions, and the processor executes the computer instructions to perform the mesh generation method for the earth-rock dam model according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions for causing the computer to execute the mesh generation method for the earth-rock dam model according to any one of claims 1 to 6.

10. A computer program product, characterized in that, Includes computer instructions for causing a computer to execute the mesh generation method for the earth-rock dam model according to any one of claims 1 to 6.

Images