A method, apparatus and equipment for generating parametric meshes for coupled stratigraphic models

Abstract

This invention provides a method, apparatus, and device for generating parametric meshes for a coupled geological structure model. The method includes: constructing a lining profile and determining the inner and outer edges of the coupled region; for the coupled region, calculating the centroid of the lining enclosure surface and the centroidal angles of each node on the inner and outer edges; unifying the density of the inner and outer edges of the coupled region, and segmenting relatively sparse edges and generating mesh points; generating radial mesh points of the coupled region based on a pre-programmed, progressively varying mesh size algorithm; dividing the coupled region into a planar mesh to obtain a two-dimensional planar mesh or a three-dimensional solid mesh, and combining this with the lining structure to generate a target mesh model. This invention effectively balances the computational precision and efficiency requirements of numerical analysis models, solving the problems of difficult mesh generation and poor mesh quality in existing related technologies.

Description

Technical Field

[0001] This invention relates to the field of geological structure modeling technology, and in particular to a method, apparatus and equipment for generating parametric meshes for coupled geological structure models. Background Technology

[0002] The existing engineering analogy method for tunnel structure design is increasingly unable to meet the increasingly complex and stringent requirements of tunnel structure design, especially in deep underground spaces where the synergistic bearing effect of surrounding rock and structure is significant, leading to a sharp increase in the demand for coupled modeling and analysis of strata and structure in tunnel structure design.

[0003] A prerequisite for coupled stratigraphy-structure calculations is the completion of coupled stratigraphy-structure modeling and mesh generation. However, existing techniques for mesh generation of coupled geological-structure models have the following problems:

[0004] 1. Complex geometric regions present significant challenges in mesh generation, sometimes even making it impossible. The mesh generation region of the stratigraphic-structure coupled model consists of the tunnel structure outline and the stratigraphic outline, such as... Figure 1 As shown, Figure 1 This is a schematic diagram of the mesh generation of the stratigraphic-structure coupled model. The inner contour lines of the region (i.e., the outer edge of the tunnel structure) are diverse, including straight-wall arched sections, single-centered circles, three-centered circles, five-centered circles, horseshoe shapes, and many others. The outer contour lines of the region (geological body edges) contain complex stratigraphic boundaries, making the region's shape extremely irregular. This presents difficulties for finite element mesh generation, and existing mesh generation methods may even fail, making it highly dependent on manual intervention.

[0005] 2. Lack of methods for controlling grid-scale distribution. Figure 2 This is a diagram illustrating a case of poor grid size distribution, such as... Figure 2 As shown, in the coupled strata-structure numerical model, design and construction technicians exhibit significant differences in their focus on different regions within the surface area. The tunnel structure and its surrounding strata are often the areas of primary concern for technicians, requiring a more precise mesh for accurate calculations and analysis. Conversely, the surrounding rock far from the tunnel structure does not require a small mesh, as this would severely consume computational resources and increase analysis time. Existing mesh replacement techniques mostly produce relatively uniform meshes, which fail to reflect the technicians' requirements for computational accuracy and efficiency.

[0006] 3. Mesh quality control is difficult, and mesh adaptability is lacking. Existing mesh generation / generation algorithms generally use triangles (surface elements) or tetrahedrons (volume elements) for complex geometries. Complex geometries are often mixed meshes of "triangles + quadrilaterals, or pure triangles (surface elements)" or "tetrahedrons + hexahedrons, or pure tetrahedrons (volume elements)". These meshes are generally suitable for finite element methods (e.g., ANSYS), but difficult to use for finite difference methods (e.g., FLAC). In numerical analysis of underground engineering, finite difference methods require relatively regular structured meshes, which necessitates that the mesh be entirely composed of quadrilaterals (surface elements) and hexahedrons (volume elements). In this case, the unstructured meshes generated by existing mesh generation algorithms are difficult to apply to numerical analysis using finite difference methods such as FLAC.

[0007] There is currently no effective solution to the problems of difficult mesh generation and poor mesh quality in existing related technologies. Summary of the Invention

[0008] This invention provides a parametric mesh generation method, apparatus, and equipment for stratigraphic structure coupling models, which solves the defects of existing related technologies such as high difficulty in mesh division and poor mesh quality, and realizes the generation of parametric, high-efficiency three-dimensional meshes for stratigraphic-structure coupling analysis models under complex geological bodies and complex lining contours.

[0009] In a first aspect, the present invention provides a method for generating parametric meshes for a coupled geological structure model, comprising:

[0010] Construct the lining profile and determine the inner and outer edges of the coupling region; both the inner and outer edges of the coupling region contain multiple nodes;

[0011] For the coupling region, calculate the centroid of the lining enclosure surface and the centroidal angles of each node on the inner and outer edges;

[0012] The density of the inner and outer edges of the coupling region is unified, and the relatively sparse edges are divided into segments and grid points are generated.

[0013] Based on a pre-programmed algorithm for progressively varying mesh sizes, radial mesh points are generated in the coupled region.

[0014] The coupling region is divided into planar parts to obtain a two-dimensional planar mesh or a three-dimensional solid mesh, and combined with the lining structure, a target mesh model is generated.

[0015] According to the present invention, a parametric mesh generation method for a coupled stratigraphic model is provided, which determines the inner boundary lines of the coupled region, including:

[0016] Obtain the required number of grid cells for the coupled region;

[0017] Based on the geometric features of the outer edge of the lining, the nodes of the inner edge line of the coupling region are directly drawn by approximating the arc with multiple straight line segments.

[0018] According to the parametric mesh generation method for a coupled geological structure model provided by the present invention, if the lining profile requires prior analysis, the inner boundary line of the coupled region is determined, including:

[0019] Obtain existing lining mesh information and determine the inner edge of the coupling region based on the lining mesh information.

[0020] According to the present invention, a parametric mesh generation method for a coupled stratigraphic model is provided, which determines the outer boundary of the coupled region, including:

[0021] Obtain stratigraphic information based on geological survey data or geological platform data;

[0022] The outer boundary line is generated based on the stratigraphic information of the coupling region.

[0023] According to the parametric mesh generation method for a coupled geological structure model provided by the present invention, the centroid of the lining enclosure surface is calculated, including:

[0024] Obtain the area of ​​the region enclosed by the inner edge of the coupling region, and determine the centroid coordinates of the lining enclosure surface.

[0025] According to the present invention, a parametric mesh generation method for a coupled stratigraphic model includes unifying the density of the inner and outer edges of the coupled region, comprising:

[0026] The relatively sparse edges among the inner and outer edges of the coupling region are defined as target edges; the nodes on the target edges are defined as target nodes.

[0027] A double traversal loop algorithm is set up to traverse the target edge and non-target edge, and mark the newly added nodes on the target edge.

[0028] For the target edge line, combined with the newly added marked nodes, a list of segmented segments after centroidal angle unification is formed.

[0029] According to the present invention, a parametric mesh generation method for a coupled stratigraphic model is provided, wherein the double traversal loop algorithm includes:

[0030] Traverse the target nodes along the target edge;

[0031] Traverse the non-target nodes of the non-target edge and compare the centroid angle of each non-target node with the current target node;

[0032] Based on the comparison results of the centroid angle, a new node needs to be added to mark the target edge.

[0033] According to the present invention, a parametric mesh generation method for a coupled stratigraphic model is provided, wherein the mesh size gradually varies stepwise, comprising:

[0034] Determine the gradient control parameters for the mesh size; the gradient control parameters include the number of mesh divisions and the mesh gradient ratio at each level;

[0035] Clearly define the start and end points of the grid segmentation;

[0036] Traverse the inner edge nodes of the coupling region and calculate the coordinates of the dividing points on each radial line to form radial grid points.

[0037] Secondly, the present invention also provides a parametric mesh generation device for a coupled geological structure model, comprising:

[0038] A construction module is used to construct the lining profile and determine the inner and outer edges of the coupling region; both the inner and outer edges of the coupling region contain multiple nodes.

[0039] The processing module is used to calculate the centroid of the lining enclosure surface and the centroidal angles of each node on the inner and outer edges for the coupling region.

[0040] A unified module is used to unify the density of the inner and outer edges of the coupling region, and to segment and evenly divide the relatively sparse edges and generate grid points.

[0041] The segmentation module is used to generate radial grid points in the coupled region based on a pre-programmed, progressively varying grid size segmentation algorithm.

[0042] The generation module is used to divide the coupling region into planes to obtain a two-dimensional planar mesh or a three-dimensional solid mesh, and combine it with the lining structure to generate a target mesh model.

[0043] Thirdly, the present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the parametric mesh generation method for the coupled stratigraphic model as described in the first aspect above.

[0044] Fourthly, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the parameterized mesh generation method for the coupled stratigraphic model as described in the first aspect above.

[0045] Fifthly, the present invention also provides a computer program product, including a computer program that, when executed by a processor, implements the parametric mesh generation method for the coupled geological structure model as described in the first aspect above.

[0046] Compared with the prior art, the present invention has the following beneficial effects:

[0047] The parametric mesh generation method for the stratigraphic structure coupling model provided by this invention achieves mesh regularization and high mesh size by unifying mesh density based on the centroid angle comparison of inner and outer geometric edges (boundaries). By setting a radial gradient ratio for the mesh size, a reasonable distribution of mesh precision is achieved, with fine mesh size and high mesh density for the tunnel lining structure and adjacent strata, while the mesh for distant strata far from the tunnel lining structure consumes less computational resources. This method highly balances the computational precision and efficiency requirements of the numerical analysis model, solving the problems of difficult mesh generation and poor mesh quality in existing related technologies. Furthermore, the stratigraphic-structure coupling model mesh generated by this method is a quadrilateral / hexahedral structured mesh, which is applicable not only to finite element method calculations represented by ANSYS, but also to finite difference method calculations represented by FLAC, which is widely used in underground engineering. The mesh has strong adaptability and a wide range of applications. Attached Figure Description

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

[0049] Figure 1 This is a schematic diagram of the mesh generation of the stratigraphic-structure coupling model.

[0050] Figure 2 This is a schematic diagram illustrating a case of poor grid size distribution;

[0051] Figure 3 This is a flowchart of the parametric mesh generation method for the coupled stratigraphic model provided by the present invention;

[0052] Figure 4 This is a schematic diagram of the outer edge grid nodes of the tunnel lining in an embodiment of the present invention;

[0053] Figure 5 This is a schematic diagram of the cross-sectional parameters of the five-centered circular tunnel in an embodiment of the present invention;

[0054] Figure 6 This is a schematic diagram of the central angle of the outer edge line within the coupling region in an embodiment of the present invention;

[0055] Figure 7 This is a schematic diagram illustrating the uniform density of outer edges in an embodiment of the present invention;

[0056] Figure 8 This is a schematic diagram of the mesh size gradually changing segmentation algorithm in an embodiment of the present invention;

[0057] Figure 9 This is a schematic diagram of the grid points in the coupling region in an embodiment of the present invention;

[0058] Figure 10 This is a two-dimensional grid diagram of the stratigraphic-structure coupling region in an embodiment of the present invention;

[0059] Figure 11 This is a schematic diagram of a three-dimensional mesh of the stratigraphic-structure coupling region in an embodiment of the present invention;

[0060] Figure 12 This is a schematic diagram of the three-dimensional mesh of the stratigraphic-structure coupling analysis model in an embodiment of the present invention;

[0061] Figure 13 This is a structural block diagram of the parametric mesh generation device for the coupled stratigraphic model provided by the present invention;

[0062] Figure 14 This is a schematic diagram of the structure of the electronic device provided by the present invention. Detailed Implementation

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

[0064] This invention provides a method for generating parametric meshes for coupled geological structures. Figure 3 This is a flowchart of the parametric mesh generation method for the coupled stratigraphic model provided by the present invention, such as... Figure 3 As shown, the method includes the following steps:

[0065] Step S301: Construct the lining profile and determine the inner and outer edges of the coupling region; both the inner and outer edges of the coupling region contain multiple nodes.

[0066] Step S302: For the coupling region, calculate the centroid of the lining enclosure surface and the centroid angle of each node on the inner and outer edges;

[0067] Step S303: Unify the density of the inner and outer edges of the coupling region, and divide the relatively sparse edges into segments and generate grid points.

[0068] Step S304: Based on a pre-programmed mesh size progressively varying segmentation algorithm, generate radial mesh points in the coupling region;

[0069] Step S305: Divide the coupled region into planar parts to obtain a two-dimensional planar mesh or a three-dimensional solid mesh, and combine it with the lining structure to generate a target mesh model.

[0070] In this method, firstly, the lining contour is constructed, and the inner and outer edges of the coupling region are generated, and mesh nodes are drawn. Then, the centroid of the lining enclosure surface and the centroidal angle of each node are calculated. In most cases, the number of mesh nodes on the inner and outer edges is inconsistent and enormous. To obtain a higher quality coupling mesh, the density of the inner and outer edges of the coupling region is unified, and the relatively sparse edges are then segmented and mesh points are generated. Next, based on a pre-programmed progressively varying mesh size subdivision algorithm, radial mesh points of the coupling region are generated. Finally, the coupling region is planarized to obtain a two-dimensional planar mesh or a three-dimensional solid mesh, which, combined with the lining structure, generates the target mesh model. Through the above process, a parameterized and highly efficient three-dimensional mesh generation of a stratigraphic-structural coupling analysis model under a complex geological body and complex lining contour is achieved. The generation process is highly automated, providing solid modeling technology support for intelligent design of underground spaces.

[0071] In the above process, by unifying the mesh density based on the centroidal angle comparison of the inner and outer geometric edges (boundaries), the mesh shape is regularized and the mesh size is relatively large. By setting the radial gradient ratio of the mesh size, a reasonable distribution of mesh precision is achieved, with fine mesh size and high mesh density for the tunnel lining structure and adjacent strata, while the mesh for distant strata far from the tunnel lining structure consumes less computational resources. This highly balances the computational precision and efficiency requirements of the numerical analysis model, solving the problems of difficult mesh generation and poor mesh quality in existing related technologies. In addition, the strata-structure coupling model mesh generated by this method is a quadrilateral / hexahedral structured mesh, which is not only applicable to finite element method calculations represented by ANSYS, but also to finite difference method calculations represented by FLAC, which is widely used in underground engineering. The mesh has strong adaptability and wide applicability.

[0072] In some embodiments, step S301, determining the inner edge of the coupling region, includes: obtaining the required number of grids for the coupling region; and, in combination with the geometric features of the outer edge of the lining, directly drawing the nodes of the inner edge of the coupling region by using multiple straight line segments to approximate the arc.

[0073] For example, based on the geometric characteristics of various tunnel cross-sections such as straight-wall arched sections, single-centered circles, three-centered circles, five-centered circles, and horseshoe shapes, grid nodes are constructed on the outer edge of the tunnel lining, such as... Figure 4 As shown, Figure 4 This is a schematic diagram of the outer edge grid nodes of the tunnel lining in an embodiment of the present invention. The total number of grids is [number missing]. n al .

[0074] Taking a five-centered circular cross-section as an example, the parameters and shape of the five-centered circular end face are as follows: Figure 5 As shown, Figure 5 This is a schematic diagram of the cross-sectional parameters of a five-centered circular tunnel in an embodiment of the present invention. Key parameters in the diagram include: O a ~ O d Let be the center of each arc segment; r a ~ r d Let be the radius of each arc segment; φ a ~ φ d Let be the radian angle of each arc segment; H The burial depth of deep-earth structures. W t The total width of the cross section H t This represents the total height of the cross-section.

[0075] The lining edge mesh nodes can be directly established using various programming languages ​​such as Python and C++, employing differential techniques to approximate the curve with multiple straight line segments. The number of straight line segments represents the total number of mesh nodes for the lining line. n al The specific steps are as follows:

[0076] 1. The first arc segment ( ):

[0077]

[0078] in, The angle is represented as The x-coordinate of the point. This represents the radius of the first arc segment. Indicates an angle in radians; The angle is represented as The z-coordinate of the point. H Indicates the burial depth of deep-earth structures.

[0079] 2. The second arc segment ( ):

[0080]

[0081] in, This represents the radian angle of the first arc segment. This indicates the radius of the second arc segment.

[0082] 3. The third arc segment ( ):

[0083]

[0084] in, This indicates the radian angle of the second arc segment. This indicates the radius of the third arc segment.

[0085] 4. The fourth arc segment ( ):

[0086]

[0087] in, This represents the radian angle of the third arc. This represents the radian angle of the fourth arc. This indicates the radius of the third arc segment.

[0088] The coordinates of the center of the fourth arc satisfy:

[0089]

[0090] in, This represents the coordinates of the center of the fourth arc. The coordinates of the fourth arc can be calculated using the following formula:

[0091]

[0092] In the above process, ~ , ~ , H , W t , H t These are the input parameters for parametric modeling techniques. When the actual project involves other types of cross-sections such as straight-walled arched sections, single-centered circles, three-centered circles, or horseshoe shapes, the corresponding geometric control parameters are the input parameters for generating the edge mesh within the coupled region.

[0093] In another implementation, if there is a need for preliminary analysis of the lining profile, the inner boundary line of the coupling region is determined, including: obtaining existing lining mesh information and determining the inner boundary line of the coupling region based on the lining mesh information.

[0094] In practical engineering, there may be situations where there is a need for preliminary analysis of tunnel lining, and the lining has already formed a mesh, requiring supplementary analysis of the strata-lining structure coupling. In this case, the user can import existing lining mesh information as the inner boundary line of the coupling region.

[0095] In some embodiments, step S301, determining the outer boundary of the coupling region, includes: obtaining stratigraphic information based on geological survey data or geological platform data; and generating the outer boundary based on the stratigraphic information of the coupling region.

[0096] In this embodiment, based on geological survey data or geological platform data, stratigraphic information is read in to form the outer contour of the coupling region, see... Figure 1 The outer edge of the coupling region.

[0097] Based on this, step S302 calculates the centroid of the lining enclosure surface, including: obtaining the area of ​​the area enclosed by the inner edge of the coupling region, and determining the centroid coordinates of the lining enclosure surface.

[0098] For example, the centroid coordinates x and z of the lining enclosure surface are calculated according to the following formula:

[0099]

[0100] in, x Indicates the centroid x coordinate, z Indicates the centroid z coordinate, Indicates coordinates as x The corresponding point z coordinate, A This indicates the area enclosed by the outer edge of the lining (the inner edge of the coupling area).

[0101] Furthermore, depending on whether a known point belongs to the inner or outer edge, calculate the centroid angle of each known point on the inner or outer edge respectively. Figure 6 This is a schematic diagram of the central angle of the outer and inner edges of the coupling region in an embodiment of the present invention, as shown below. Figure 6 As shown, the centroidal angle of each known node on the edge within the coupling region is... The centroidal angle of each known node on the outer edge of the coupling region is... .

[0102] In some embodiments, step S303, unifying the density of the inner and outer edges of the coupling region, includes: determining the relatively sparse edges among the inner and outer edges of the coupling region as target edges; the nodes on the target edges are target nodes; setting a double traversal loop algorithm to traverse the target edges and non-target edges and mark the newly added nodes on the target edges; for the target edges, combining the marked newly added nodes, forming a list of segmented numbers after centroidal angle unification.

[0103] Specifically, the double traversal loop algorithm includes: traversing the target nodes of the target edge; traversing the non-target nodes of the non-target edge and comparing the centroid angle of each non-target node with the current target node; and marking the target edge as needing to add a new node based on the centroid angle comparison results.

[0104] This embodiment can be applied to situations where the known points on the inner edge are dense, while those on the outer edge are sparse; conversely, it can also be applied. Considering that in practical engineering, the known points on the inner edge are mostly dense, while those on the outer edge are sparse, this case is used as an example to introduce the density unification algorithm of this embodiment.

[0105] To identify relatively sparse edges with known points, the following describes a double-traversal loop algorithm using outer edges as an example:

[0106] First step: Traverse the known points on the outer edge, corresponding to the... j Second traversal (known points on the outer edge) j Then, perform the second traversal.

[0107] The second step: Traverse the known points of the inner boundary line and compare each known point of the inner boundary line. i With the known points of the current outer boundary (first type) j The size of the centroidal angle. When When marking the outer edge, one new node needs to be added.

[0108] After completing the above two loops, the new node marking on the outer edge can be completed.

[0109] Meanwhile, the uniformity of edge density achieved through centroidal angles can realize good mesh size uniformity, fundamentally avoiding the possibility of distorted meshes and ensuring mesh quality.

[0110] For sparse edges, a list of segmentation numbers is generated after unifying the centroidal angles of the inner and outer edges. The following is a detailed explanation:

[0111] Figure 7 This is a schematic diagram illustrating the uniform density of outer edges in an embodiment of the present invention, as shown below. Figure 7For the outer edge shown, segment BC has only two known points, namely points B and C; while within the centroidal angle of segment BC, there are six known points on the inner edge. Therefore, using the aforementioned dual algorithm, four additional points can be marked on segment BC, achieving the following:

[0112] 1. The density of the inner and outer edges is uniform;

[0113] 2. The centroid angles of the grid points are evenly distributed.

[0114] Based on the above embodiments, the mesh size gradual change segmentation algorithm includes: determining the gradual change control parameters of the mesh size; the gradual change control parameters include the number of mesh segments and the gradual change ratio of the mesh; determining the starting point and ending point of the mesh segmentation; traversing the inner edge nodes of the coupling region and calculating the coordinates of the segmentation points on each radial line to form radial mesh points.

[0115] To achieve a reasonable distribution of mesh precision, a mesh size gradually varying subdivision algorithm was developed, such as... Figure 8 As shown, Figure 8 This is a schematic diagram of the mesh size gradually changing segmentation algorithm in an embodiment of the present invention, specifically:

[0116] 1. Input of gradual change control parameters.

[0117] This embodiment allows users to input grid size gradient control parameters: number of grid divisions. num Gradual scaling of the grid rat For grids with progressively varying scales rat The first one, far from the inner edge (tunnel lining structure) i The radial side length of each grid is l i , No. i +1 grid radial side length is l i+1 Then we have:

[0118]

[0119] in, rat This indicates the gradual change in grid scale, specifically the ratio of the lengths of two adjacent segments after division. l i+1 Indicates the first i +1 grid radial side length, l i Indicates the first i Each grid has a radial side length.

[0120] 2. Identify the known endpoints and denote the starting point as... M The endpoint is recorded as N .

[0121] 3. At this time, the first iThe coordinates of each dividing point can be determined recursively using the following formula:

[0122]

[0123] in, Indicates the first i Dividing points x coordinate, Indicates the first i Dividing points z coordinate, Indicates the starting point x coordinate, Indicates the end point x coordinate, Indicates the first i The grid scale gradually changes with each subdivision point starting at -1. The z-coordinate of the starting point. The z-coordinate of the endpoint. Indicates the first j The grid scale gradually changes with each dividing point -1.

[0124] 4. Traverse the inner edge nodes, and use the above formula to calculate the coordinates of the dividing points on each radial line to form grid points, such as... Figure 9 As shown, Figure 9 This is a schematic diagram of the grid points in the coupling region in an embodiment of the present invention.

[0125] In practical applications of this method, users can choose to generate either a two-dimensional planar mesh or a three-dimensional solid mesh, depending on their needs. When the user selects two-dimensional planar mesh generation, the coupling region is divided into two-dimensional planes, such as... Figure 10 As shown, Figure 10 This is a two-dimensional grid diagram of the stratigraphic-structure coupling region in an embodiment of the present invention. (Through...) Figure 10 It can be seen that the stratigraphic-structure coupled surface region partitioning algorithm of this method can generate a regular-shaped unit mesh with gradually changing radial dimensions, and all the meshes are quadrilateral structured meshes, which can be used for calculation of both finite element method and finite difference method numerical analysis models.

[0126] When the user selects 3D mesh generation, the coupled region is divided into 3D planes. Specifically: along the axial coordinates perpendicular to the tunnel cross-section... y For those that need to be modeled y By repeating the planar mesh generation process described above for cross-sections, mesh generation for 3D models of stratigraphic-structure coupling analysis under complex 3D terrain can be achieved, such as... Figure 11 As shown, Figure 11 This is a three-dimensional mesh diagram of the stratigraphic-structure coupling region in an embodiment of the present invention.

[0127] Will Figure 11 By adding the stratigraphic-structure coupled geological model mesh to the tunnel lining structure model, a high-quality mesh for generating a complete three-dimensional numerical analysis model of the tunnel structure and complex geology can be formed, such as... Figure 12 As shown, Figure 12 This is a schematic diagram of the three-dimensional mesh of the stratigraphic-structure coupling analysis model in an embodiment of the present invention.

[0128] In summary, this method can be widely applied to various scenarios involving coupled strata and tunnel structures, and it possesses certain technical advantages in scenarios requiring coupled strata and structure modeling and analysis, such as the construction of deep underground tunnels. This method can be extended to mesh generation and mesh fineness distribution control for various geometries containing complex inner and outer boundary lines and multi-connected geometric regions. Furthermore, the mesh generated by this method is applicable to both finite element method and finite difference method models.

[0129] The present invention also provides a parametric mesh generation device for a coupled stratigraphic model. The parametric mesh generation device for a coupled stratigraphic model provided by the present invention will be described below. The parametric mesh generation device for a coupled stratigraphic model described below can be referred to in correspondence with the parametric mesh generation method for a coupled stratigraphic model described above. Figure 13 This is a structural block diagram of the parametric mesh generation device for the coupled stratigraphic model provided by the present invention, as shown in the figure. Figure 13 As shown, the device includes:

[0130] Module 1301 is used to construct the lining profile and determine the inner and outer edges of the coupling region; both the inner and outer edges of the coupling region contain multiple nodes.

[0131] Processing module 1302 is used to calculate the centroid of the lining enclosure surface and the centroid angle of each node on the inner and outer edges for the coupled region.

[0132] The unified module 1303 is used to unify the density of the inner and outer edges of the coupling region, and to divide the relatively sparse edges into segments and generate grid points.

[0133] The segmentation module 1304 is used to generate radial grid points in the coupled region based on a pre-programmed progressively varying grid size segmentation algorithm.

[0134] The generation module 1305 is used to divide the coupled region into planes to obtain a two-dimensional planar mesh or a three-dimensional solid mesh, and to generate a target mesh model in combination with the lining structure.

[0135] In use, this device first constructs the lining outline and generates the inner and outer edges of the coupling region, drawing mesh nodes. Then, the processing module 1302 calculates the centroid of the lining enclosure surface and the centroidal angle of each node. In most cases, the number of mesh nodes on the inner and outer edges is inconsistent and enormous. To obtain a higher quality coupling mesh, the unification module 1303 unifies the density of the inner and outer edges of the coupling region, then segments the relatively sparse edges and generates mesh points. Next, the segmentation module 1304 generates radial mesh points in the coupling region based on a pre-programmed, progressively varying mesh size algorithm. Finally, the generation module 1305 divides the coupling region into planar parts, obtaining a two-dimensional planar mesh or a three-dimensional solid mesh, and combines this with the lining structure to generate the target mesh model. Through this process, a parameterized and efficient three-dimensional mesh generation of a stratigraphic-structure coupling analysis model under a complex geological body and complex lining outline is achieved. The generation process is highly automated, providing solid modeling technology support for intelligent design of underground spaces.

[0136] In the above process, the mesh density was unified by comparing the centroid angles of the inner and outer geometric edges (boundaries), achieving mesh regularization and a high mesh size. By setting a radial gradient ratio for the mesh size, a reasonable distribution of mesh precision was achieved, with fine mesh size and high mesh density for the tunnel lining structure and adjacent strata, while the mesh for distant strata far from the tunnel lining structure consumed fewer computational resources. This effectively balanced the computational precision and efficiency requirements of the numerical analysis model, solving the problems of difficult mesh generation and poor mesh quality in existing related technologies. Furthermore, the strata-structure coupling model mesh generated by this device is a quadrilateral / hexahedral structured mesh, which is applicable not only to finite element method calculations represented by ANSYS, but also to finite difference method calculations represented by FLAC, which is widely used in underground engineering. The mesh has strong adaptability and a wide range of applications.

[0137] Figure 14 An example is a schematic diagram of the physical structure of an electronic device, such as... Figure 14 As shown, the electronic device may include: a processor 1401, a communication interface 1402, a memory 1403, and a communication bus 1404, wherein the processor 1401, the communication interface 1402, and the memory 1403 communicate with each other via the communication bus 1404. The processor 1401 can call logical instructions in the memory 1403 to execute a parametric mesh generation method for a coupled geological structure model, the method including:

[0138] Construct the lining profile and determine the inner and outer edges of the coupling region; both the inner and outer edges of the coupling region contain multiple nodes.

[0139] For the coupled region, calculate the centroid of the lining enclosure surface and the centroidal angles of each node on the inner and outer edges;

[0140] The density of the inner and outer edges of the coupling region is unified, and the relatively sparse edges are segmented and grid points are generated.

[0141] Based on a pre-programmed algorithm for progressively varying mesh size, radial mesh points in the coupling region are generated.

[0142] The coupling region is divided into planar parts to obtain a two-dimensional planar mesh or a three-dimensional solid mesh, and combined with the lining structure, a target mesh model is generated.

[0143] Furthermore, the logical instructions in the aforementioned memory 1403 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, essentially, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0144] On the other hand, the present invention also provides a computer program product, which includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer is able to execute the parameterized mesh generation method for the coupled stratigraphic model provided by the above methods, the method comprising:

[0145] Construct the lining profile and determine the inner and outer edges of the coupling region; both the inner and outer edges of the coupling region contain multiple nodes.

[0146] For the coupled region, calculate the centroid of the lining enclosure surface and the centroidal angles of each node on the inner and outer edges;

[0147] The density of the inner and outer edges of the coupling region is unified, and the relatively sparse edges are segmented and grid points are generated.

[0148] Based on a pre-programmed algorithm for progressively varying mesh size, radial mesh points in the coupling region are generated.

[0149] The coupling region is divided into planar parts to obtain a two-dimensional planar mesh or a three-dimensional solid mesh, and combined with the lining structure, a target mesh model is generated.

[0150] In another aspect, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for generating parameterized meshes for coupled stratigraphic models provided by the methods described above, the method comprising:

[0151] Construct the lining profile and determine the inner and outer edges of the coupling region; both the inner and outer edges of the coupling region contain multiple nodes.

[0152] For the coupled region, calculate the centroid of the lining enclosure surface and the centroidal angles of each node on the inner and outer edges;

[0153] The density of the inner and outer edges of the coupling region is unified, and the relatively sparse edges are segmented and grid points are generated.

[0154] Based on a pre-programmed algorithm for progressively varying mesh size, radial mesh points in the coupling region are generated.

[0155] The coupling region is divided into planar parts to obtain a two-dimensional planar mesh or a three-dimensional solid mesh, and combined with the lining structure, a target mesh model is generated.

[0156] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0157] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0158] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for generating parametric meshes for a coupled geological structure model, characterized in that, The method comprises the following steps: constructing a lining profile and determining an inner edge line and an outer edge line of a coupling area; the inner edge line and the outer edge line of the coupling area each comprise a plurality of nodes; for the coupling area, calculating a centroid of a lining enclosing surface and a centroid angle of each node on the inner edge line and the outer edge line; uniformly adjusting the density of the inner edge line and the outer edge line of the coupling area, and segmentally and evenly dividing and generating grid points on the relatively sparse edge line; based on a pre-prepared grid size step-by-step gradual segmentation algorithm, generating radial grid points of the coupling area; plane-dividing the coupling area to obtain a two-dimensional plane grid or a three-dimensional solid grid, and combining the lining structure to generate a target grid model; uniformly adjusting the density of the inner edge line and the outer edge line of the coupling area comprises: determining a relatively sparse edge line of the inner edge line and the outer edge line of the coupling area as a target edge line; the nodes on the target edge line are target nodes; setting a double traversal loop algorithm to traverse and mark the target edge line and the non-target edge line, and add new nodes on the target edge line; for the target edge line, combining the marked new nodes to form a segmented segmentation number list after the centroid angle is adjusted; the grid size step-by-step gradual segmentation algorithm comprises: determining a gradual control parameter of the grid size; the gradual control parameter comprises a grid segmentation number and a grid step-by-step gradual ratio; defining a starting point and an ending point of the grid segmentation; traversing the inner edge line nodes of the coupling area, and calculating the coordinates of the segmentation points on each radial line one by one to form radial grid points; the double traversal loop algorithm comprises: traversing the target nodes of the target edge line; traversing the non-target nodes of the non-target edge line, and comparing the centroid angle of each non-target node with the current target node; according to the comparison result of the centroid angle, marking that a new node needs to be added on the target edge line.

2. The method of claim 1, wherein, determining the inner edge line of the coupling area comprises: obtaining the grid quantity requirement of the coupling area; combining the geometric characteristics of the lining outer edge, adopting a multi-segment straight line segment approximation arc line mode, and directly drawing the nodes of the inner edge line of the coupling area.

3. The method of claim 1, wherein, if the lining profile has a prior analysis requirement, determining the inner edge line of the coupling area comprises: obtaining existing lining grid information, and determining the inner edge line of the coupling area according to the lining grid information.

4. The method of claim 1, wherein, determining the outer edge line of the coupling area comprises: obtaining stratum information according to geological exploration data or geological platform data; based on the stratum information of the coupling area, generating the outer edge line.

5. The method for generating a parameterized grid for a formation structure coupling model according to claim 1, wherein, calculating the centroid of the lining enclosing surface comprises: obtaining the area of the inner edge line enclosing area of the coupling area, and determining the centroid coordinates of the lining enclosing surface.

6. An apparatus for generating a parameterized grid for a coupled model of a formation structure, for implementing the method for generating a parameterized grid for a coupled model of a formation structure according to any one of claims 1 to 5, characterized in that, comprise: a construction module, configured to construct a lining profile and determine an inner edge line and an outer edge line of a coupling area; the inner edge line and the outer edge line of the coupling area each comprise a plurality of nodes; a processing module, configured to calculate, for the coupling area, a centroid of a lining enclosing surface and a centroid angle of each node on the inner edge line and the outer edge line; a uniformity module, configured to uniformly adjust the density of the inner edge line and the outer edge line of the coupling area, and segmentally and evenly divide and generate grid points on the relatively sparse edge line; The segmentation module is configured to generate radial grid points of the coupling area based on a pre-programmed grid size step-by-step gradual segmentation algorithm. The generation module is configured to perform planar division on the coupling area to obtain a two-dimensional planar grid or a three-dimensional solid grid, and generate a target grid model in combination with a lining structure.

7. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor implements the method for generating a parameterized grid of a formation structure coupling model according to any one of claims 1 to 5 when executing the program.

Images