Multi-source interface coupling transition modeling method considering spatial incongruous arched caverns
By employing a multi-source interface coupling transition modeling method, combined with the OCTREE and TETRA algorithms, the problem of efficient and high-quality discretization of spatially anisotropic arched cavern structures in complex geotechnical engineering was solved, enabling rapid generation of fine mesh models and improving modeling efficiency and computational accuracy.
Patent Information
- Application Number
- CN202511736131.0
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-11-25
- Publication Date
- 2026-02-06
AI Technical Summary
Existing technologies are insufficient for efficiently and effectively discretely analyzing complex spatially oriented arched and circular cavern structures in geotechnical engineering, resulting in low simulation analysis efficiency and difficulty in meeting engineering analysis requirements.
A multi-source interface coupling transition modeling method is adopted, which combines the OCTREE and TETRA algorithms. Through geometric entity creation, region structured segmentation, cross-scale discretization and mesh quality diagnosis, a fine mesh model is generated to ensure the continuity of mesh quality and computational model.
It improves modeling efficiency, generates high-quality fine mesh models, solves the problems of low discretization efficiency and difficulty in meeting the needs of engineering analysis in traditional methods, simplifies the modeling process, and improves calculation speed and accuracy.
Smart Images
Figure CN121480189A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The application belongs to the technical field of geotechnical engineering, and specifically relates to a multi-source interface coupling transition modeling method considering spatial anisotropic arched caverns. BACKGROUND
[0002] With the accelerated construction of modern transportation and energy systems, numerous cross-sea tunnels, underground cavern structures, and water conservancy hubs and other major geotechnical projects have emerged. Among them, the cavern structure, as a key component of these large-scale engineering systems, undertakes important tasks such as local grouting reinforcement, maintenance, and functional health monitoring. Precise control of its safety state is crucial for the smooth operation of the main project.
[0003] With the continuous improvement of construction technology and scientific theory, geotechnical engineering construction scale is increasingly large, and structure design is becoming more complex. The unique nature of the cavern type continues to require higher numerical simulation analysis. For a long time, isoparametric elements have been widely used by engineering researchers due to their simplicity and versatility. However, their element shape is limited to hexahedron and degenerate type, making it difficult to discretize complex structures. For example, the modeling of large high earth-rock dams may take up to 80% of the total simulation time, severely restricting the efficiency of numerical analysis.
[0004] From the existing research progress, the classic tetrahedral discretization method (TETRA) has good adaptability to complex geometric boundaries and can achieve good automatic subdivision. However, due to the large volume of large geotechnical projects such as submarine tunnels and high earth-rock dams, the number of elements required for discretization is large, resulting in high memory consumption and long discretization time. In addition, the discretized model lacks key information such as horizontal layered filling boundaries and spatial joint boundaries, making it difficult to meet the simulation needs of structure staging and zoning construction. The recently developed OCTREE discretization method can achieve fast and fine cross-scale discretization of complex structures, providing an effective approach for efficient generation of complex structure grid models. However, when used for geotechnical engineering discretization considering spatial anisotropic arched and circular cavern structures, the local grid quality of the cavern component is poor, and the related performance needs to be further improved. To address the above problems, we propose a multi-source interface coupling transition modeling method considering spatial anisotropic arched caverns. SUMMARY
[0005] The purpose of this invention is to address the challenge of efficiently and effectively discretizing spatially oriented arched and circular caverns in complex structures such as high dams and marine engineering projects. It provides a multi-source interface coupling transition modeling method that considers spatially oriented arched caverns. This method effectively integrates the advantages of OCTREE's high-quality and efficient discretization and TETRA's strong adaptability to complex geometric boundaries, while avoiding the shortcomings of OCTREE in handling poor mesh quality for arched and circular structures and TETRA's difficulty in meeting the analytical requirements of simulation layering and joint boundaries. It can quickly generate fine mesh models, solving the problems of low discretization efficiency and difficulty in meeting the engineering analysis requirements of traditional methods.
[0006] This invention is implemented by considering a multi-source interface coupling transition modeling method for spatially anisotropic arched chambers, the method comprising: S10, Geometric Entity Creation and Region Structured Segmentation: Obtain engineering dimension parameters, create geometric entities based on the engineering dimension parameters, and perform region structured segmentation on the created geometric entities; S20, cross-scale OCTREE-TETRA coupled discretization: the main mesh and the cavern area mesh are created using OCTREE and hexahedral sweeping methods respectively. Through the conformal / non-conformal boundary fusion algorithm, the displacement coordination of the boundary nodes is forced to construct a stress-continuous overall calculation model. S30, Mesh Quality Diagnosis, Correction and Verification: Based on a pre-built multi-index mesh evaluation system, the fused mesh is subjected to quality detection and optimization to ensure that it meets the requirements of finite element analysis.
[0007] Preferably, the geometric entity creation and region structured segmentation method includes: S101, Obtain engineering dimension parameters, and create geometric entities based on engineering dimension parameters; When creating geometric entities based on engineering dimension parameters, an initial three-dimensional entity is generated by axial stretching based on the engineering dimension parameters; based on the upstream and downstream slope ratio of the dam body and the valley topography, the main structural entity model containing the clay core dam is constructed by positioning the working plane and Boolean cutting operations. S102, Perform region-structured segmentation on the created geometric entities; When performing regional structural segmentation on the created geometric entity, the geometric entity is divided into the main structural area and the spatially anisotropic arched cave area.
[0008] Preferably, the spatially oriented arched cavern area is divided into three structural regions according to its structural characteristics: a circular cross-section cavern body region A, a high-plasticity clay filling layer region B, and a geometric transition connection region C.
[0009] Preferably, the cross-scale OCTREE-TETRA coupled discretization method includes: S201, Load the created geometric entity, and perform cross-scale discretization of the main structure based on OCTREE; S202, perform hexahedral sweep discretization on the circular cross-section cavern body A area, the high plasticity clay filling layer B area, and the geometric transition connection area C area in the geometric entity; S203, obtains cross-scale discretization and hexahedral sweep discretization results, performs multi-source interface coupling transition discretization, and adaptively generates the transition region TETRA. S204 imports the Model-A, Tunnel-A, Tunnel-B, and Model-Transition-TETRA meshes into the mesh coupling system. Through a conformal / non-conformal boundary fusion algorithm, it forces the displacement coordination of the boundary nodes and constructs a stress-continuous overall calculation model.
[0010] Preferably, when performing cross-scale discretization of the main structure based on OCTREE, the created geometric entities are imported into the OCTREE mesh generator, and differentiated mesh sizes are set according to the mechanical response requirements of the upstream riprap zone, the core wall zone, and the downstream riprap zone to generate a cross-scale refined mesh Model-A, with the solid body of the central wall corridor marked as Model-B.
[0011] Preferably, when performing hexahedral sweep discretization on the circular cross-section cave body A area, the high-plasticity clay filling layer B area, and the geometric transition connection area C area in the geometric entity, a two-dimensional planar mesh is generated based on the intermediate cross-section mesh parameters, and the hexahedral sweep operation is performed along the cave axis direction according to the set size, respectively outputting the Tunnel-A three-dimensional mesh, Tunnel-B three-dimensional mesh, and Tunnel-C three-dimensional mesh corresponding to the circular cross-section cave body A area, the high-plasticity clay filling layer B area, and the geometric transition connection area C area.
[0012] Preferably, the method for discrete transition of multi-source interface coupling includes: S2031 automatically constructs conformal interface elements between Model-A and Model-B through a node topology mapping algorithm, extracts the core wall side node set and outputs it as a surface mesh file Interface-M-AB; S2032, similarly generate the Interface-M-BC file between Tunnel-B and Tunnel-C; S2033, a closed cavity surface is constructed based on Interface-M-AB and Interface-M-BC; S2034, execute the surface-to-solid transformation to generate the transition geometry Model-Transition; S2035 adopts the principle of one discrete part for the global line model, performs tetrahedral automatic meshing on the Model-Transition, and outputs the Model-Transition-TETRA mesh.
[0013] Preferably, the mesh quality diagnosis, correction, and verification method includes: S301 performs Jacobian determinant, twist, and warpage detection on the fused mesh; S302 performs topology optimization and iterative adjustment of node positions for low-quality cells; S303, the mesh file was verified to meet the convergence requirements of finite element analysis through trial calculations.
[0014] Compared with the prior art, the embodiments of this application have the following main advantages: The finite element modeling method for spatially oriented arched cavern structures proposed in this invention has the advantages of short processing time, high speed, and high mesh quality. It also ensures the consistency of the mesh around the spatially oriented arched cavern structure. The generated model can be solved directly using conventional methods without the need to introduce multi-point constraint equations or boundary stress-displacement coordination conditions to participate in the solution, thus improving the simplicity and efficiency of modeling and secondary modifications. It effectively solves the problems of long processing time and poor quality in finite element preprocessing for slender spatial arched structures, and can provide strong support for geotechnical engineering calculations of spatially oriented arched cavern structures.
[0015] This invention effectively combines the advantages of OCTREE's high-quality and efficient discretization and TETRA's strong adaptability to complex geometric boundaries, while avoiding the shortcomings of OCTREE in handling poor quality meshes of arched and circular structures and TETRA's difficulty in meeting the analysis requirements of simulation layering and seam boundaries. It can quickly generate fine mesh models and overcome the problems of low discretization efficiency and difficulty in meeting the engineering analysis requirements of traditional methods.
[0016] This invention can rapidly establish a high-quality, detailed model considering the entire system of foundation-cavity-dam, avoiding the shortcomings of traditional sub-model analysis or multi-point constraint analysis methods in considering the interaction between foundation and structure. Furthermore, the method is highly simple, operable, and universal, providing an effective approach for the efficient mesh generation of other similar complex spatial cavity systems. Attached Figure Description
[0017] Figure 1 This is a schematic diagram of the multi-source interface coupling transition modeling method for spatially anisotropic arched caverns provided by the present invention.
[0018] Figure 2 A schematic diagram of the initial three-dimensional solid modeling process in an embodiment of the present invention is shown.
[0019] Figure 3This diagram illustrates the region-structured segmentation and discretization of the created geometric entity.
[0020] Figure 4 A schematic diagram of cross-scale mesh generation is shown when the main structure is discretized across scales based on OCTREE.
[0021] Figure 5 A schematic diagram of the hexahedral mesh generation process for a hexahedral swept discrete cavern region is shown in an embodiment of the present invention.
[0022] Figure 6 A schematic diagram of the Model-Transition-TETRA mesh generation process for the connected transition region tetrahedron is shown.
[0023] Figure 7 A schematic diagram of the process of fusion transition connection zone and core wall / cavity interface in an embodiment of the present invention is shown.
[0024] Figure 8 A schematic diagram showing the dimensional information of a geometric model of an S-shaped tunnel in an embodiment of the present invention is provided.
[0025] Figure 9 This invention illustrates a multi-scale refined finite element mesh model of an S-shaped tunnel project coupled with OCTREE-TETRA in an embodiment of the present invention. Detailed Implementation
[0026] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs; the terminology used herein in the specification of the application is for the purpose of describing particular embodiments only and is not intended to be limiting of the application; the terms "comprising" and "having," and any variations thereof, in the specification, claims, and foregoing drawings of this application are intended to cover non-exclusive inclusion. The terms "first," "second," etc., in the specification, claims, or foregoing drawings of this application are used to distinguish different objects, not to describe a particular order.
[0027] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of this application. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.
[0028] Traditional discretization methods suffer from low efficiency and difficulty in meeting engineering analysis requirements. Furthermore, they struggle with efficient and high-quality discretization of spatially anisotropic arched and circular caverns in complex structures such as high dams and marine engineering projects. To address these issues, we propose a multi-source interface coupling transition modeling method considering spatially anisotropic arched caverns. The method first constructs a composite geometric entity based on engineering dimensional parameters and divides the structure into a main structural region and a cavern region according to analytical requirements. Then, the main structure mesh and cavern region mesh are created using OCTREE and hexahedral sweeping methods, respectively. Next, conformal interfaces are extracted between the dam body and the outer clay, and between the cavern mesh and the inner clay. Based on this, a closed cavity surface is constructed using the dual interfaces and transformed into a transition zone geometry. Then, a tetrahedral adaptive subdivision is performed on the transition zone geometry using the principle of one discretization fraction for the global line model. Finally, a non-conformal boundary fusion algorithm forces the transition zone mesh to be displaced and coordinated with the main structure mesh and the fine cavern mesh, forming a unified computational model. Finally, mesh quality diagnosis, correction, and verification are performed to obtain a computational mesh that meets engineering requirements. The finite element modeling method for spatially oriented arched cavern structures proposed in this invention has the advantages of short processing time, high speed, and high mesh quality; it effectively solves the problems of long processing time and poor quality of finite element preprocessing for slender spatial arched structures, and can provide strong support for geotechnical engineering calculations of spatially oriented arched cavern structures.
[0029] This invention provides a multi-source interface coupling transition modeling method considering spatially anisotropic arched chambers. Figure 1 This diagram illustrates the implementation flow of the multi-source interface coupling transition modeling method considering spatially anisotropic arched chambers. The method specifically includes: S10, Geometric Entity Creation and Region Structured Segmentation: Obtain engineering dimension parameters, create geometric entities based on the engineering dimension parameters, and perform region structured segmentation on the created geometric entities; The method for creating geometric entities and segmenting regions in a structured manner includes: S101, Obtain engineering dimension parameters, and create geometric entities based on engineering dimension parameters; When creating geometric entities based on engineering dimension parameters, an initial three-dimensional entity is generated by axial stretching based on the engineering dimension parameters; based on the upstream and downstream slope ratio of the dam body and the valley topography, the main structural entity model containing the clay core dam is constructed by positioning the working plane and Boolean cutting operations. S102, Perform region-structured segmentation on the created geometric entities; When performing regional structural segmentation on the created geometric entity, the geometric entity is divided into a main structural area and a spatially oriented arched cavern area. The central wall cavern area of the main structural area is filled with a solid body, Zone-M. The spatially oriented arched cavern area is divided into a circular cross-section cavern body area A (Zone-A), a high-plasticity clay filling layer area B (Zone-B), and a geometric transition connection area area C (Zone-C) according to its structural characteristics. Figure 2 A schematic diagram of the initial three-dimensional solid modeling process in an embodiment of the present invention is shown, while Figure 3 This diagram illustrates the region-structured segmentation and discretization of the created geometric entity.
[0030] S20, cross-scale OCTREE-TETRA coupled discretization: the main mesh and the cavern area mesh are created using OCTREE and hexahedral sweeping methods respectively. Through the conformal / non-conformal boundary fusion algorithm, the displacement coordination of the boundary nodes is forced to construct a stress-continuous overall calculation model. The cross-scale OCTREE-TETRA coupled discretization method includes: S201, Load the created geometric entities, and perform cross-scale discretization of the main structure based on OCTREE. During this cross-scale discretization, the created geometric entities are imported into the OCTREE mesh generator, and differentiated mesh sizes are set according to the mechanical response requirements of the upstream riprap zone, core wall zone, and downstream riprap zone, generating a cross-scale refined mesh Model-A. The solid body of the central wall corridor is labeled Model-B. Figure 4 This diagram illustrates the generation of cross-scale meshes when the main structure is discretized across scales based on OCTREE. S202, perform hexahedral sweep discretization on the circular cross-section cavern body A region, the high-plasticity clay filling layer B region, and the geometric transition connection region C region of the geometric entity; wherein, when performing hexahedral sweep discretization on the circular cross-section cavern body A region, the high-plasticity clay filling layer B region, and the geometric transition connection region C region of the geometric entity, a two-dimensional planar mesh is generated based on the intermediate cross-section mesh parameters, and the hexahedral sweep operation is performed along the cavern axis direction according to a set size, outputting the Tunnel-A three-dimensional mesh, Tunnel-B three-dimensional mesh, and Tunnel-C three-dimensional mesh corresponding to the circular cross-section cavern body A region, the high-plasticity clay filling layer B region, and the geometric transition connection region C region, respectively. Figure 5 A schematic diagram of the hexahedral mesh generation process for a hexahedral swept discrete cavern region is shown in an embodiment of the present invention.
[0031] S203 obtains cross-scale discretization and hexahedral sweep discretization results, performs transition discretization on multi-source interface coupling, and adaptively generates the transition region TETRA. Figure 6A schematic diagram of the Model-Transition-TETRA mesh generation process for the tetrahedral transition region is shown. S204 imports the Model-A, Tunnel-A, Tunnel-B, and Model-Transition-TETRA meshes into a mesh coupling system. Through a conformal / non-conformal boundary fusion algorithm, it forces displacement coordination at the boundary nodes, constructing a stress-continuous overall computational model. Figure 7 A schematic diagram of the process of fusion transition connection zone and core wall / cavity interface in an embodiment of the present invention is shown.
[0032] In this embodiment of the invention, the mesh coupling system can be a non-conformal mesh coupling system. Through the operations of combining bodies → merging boundaries → merging conformal edges and non-conformal edges, displacement coordination forced mapping is achieved, ensuring accurate matching of conformal boundary nodes, and thus generating a global-mesh model.
[0033] S30, Mesh Quality Diagnosis, Correction and Verification: Based on a pre-built multi-index mesh evaluation system, the fused mesh is subjected to quality detection and optimization to ensure that it meets the requirements of finite element analysis.
[0034] In a further preferred embodiment of the present invention, the method for discrete transition of multi-source interface coupling includes: S2031 automatically constructs conformal interface elements between Model-A and Model-B through a node topology mapping algorithm, extracts the core wall side node set and outputs it as a surface mesh file Interface-M-AB; When extracting the core wall side node set and outputting it as a surface mesh file Interface-M-AB, the shared node set between Model-A and Model-B can be automatically searched based on the spatial position tolerance algorithm. The automated operation of connecting nodes to form elements is achieved by copying node coordinates and modifying node coordinates along the element surface normal. Normal offset elements are generated directly in Model-A and Model-B. The node subset on the core wall side of the interface element is extracted through topology mapping technology, which is marked as Interface-M-AB.
[0035] S2032, similarly generate the interface file Interface-M-BC between Tunnel-B and Tunnel-C, thus obtaining the interface file between the high plasticity clay B region and the connecting transition C region, which is marked as Interface-M-BC here; S2033, a closed cavity surface is constructed based on Interface-M-AB and Interface-M-BC. The transition line and transition surface are constructed based on the boundary key points of Interface-M-AB and Interface-M-BC. Triangular surfaces are created in the end cavity area, and a closed cavity is formed by all the surface patches. S2034, perform surface-to-solid transformation to generate transition zone geometry Model-Transition, where all the triangular facets in the cavity formed by the above steps generate geometric solids, i.e. geometric solids connecting the transition zone, which are denoted as Model-Transition; S2035 adopts the principle of one discrete part for the global line model, performs tetrahedral automatic meshing on the Model-Transition, and outputs the Model-Transition-TETRA mesh.
[0036] In this embodiment, the principle of one discrete part of the global line model is adopted, and the tetrahedral discretization operation is performed to automatically generate the tetrahedral mesh of the Model-Transition region, which is denoted as Model-Transition-TETRA.
[0037] In a further preferred embodiment of the present invention, the mesh quality diagnosis, correction, and verification method includes: S301 performs Jacobian determinant, twist, and warpage detection on the fused mesh; S302 performs topology optimization and iterative adjustment of node positions for low-quality cells; S303, the mesh file was verified to meet the convergence requirements of finite element analysis through trial calculations.
[0038] The following engineering case is a preferred example of rapid and precise geotechnical engineering modeling of spatially oriented arched tunnel structures, which is a preferred example of this invention. Other examples, such as the principle framework, basic theory and ideas or implementation methods, that are the same as or similar to this example are all within the protection scope of this invention. Specifically, the engineering case is an S-shaped concrete-lined tunnel.
[0039] A circular cross-section is selected for an S-shaped concrete-lined tunnel. The main dimensions are as follows: the front and rear sections are arched upwards and downwards at certain angles, respectively. The overburden depth, the depth of the bedrock below, and the calculation range are all defined as A at both ends. The tunnel length is B, the tunnel outer diameter is C, the width of the swept hexahedral region is D, and the width of the transition zone is E. The overall structure is shown in [reference needed]. Figure 8 , Figure 8 This diagram illustrates the dimensional information of a geometric model of an S-shaped tunnel project according to an embodiment of the present invention. Applying the aforementioned modeling approach of the present invention, cross-scale fine modeling of this project was carried out, and the results are shown below.Figure 9 ,in, Figure 9 This invention illustrates a multi-scale refined finite element mesh model of an S-shaped tunnel project coupled with OCTREE-TETRA in an embodiment of the present invention.
[0040] In summary, this invention provides a multi-source interface coupling transition modeling method considering spatially anisotropic arched caverns. The finite element modeling method for spatially anisotropic arched cavern structures proposed in this invention has the advantages of short processing time, high speed, and high mesh quality, while ensuring the mesh coordination and consistency around the spatially anisotropic arched cavern structure. The generated model can be directly solved using conventional methods without the need to introduce multi-point constraint equations or boundary stress-displacement coordination conditions for the solution, thus improving the simplicity and efficiency of modeling and secondary modifications. It effectively solves the problems of long processing time and poor quality in finite element preprocessing of slender spatial arched structures, and can provide strong support for geotechnical engineering calculations of spatially anisotropic arched cavern structures.
[0041] This invention effectively combines the advantages of OCTREE's high-quality and efficient discretization and TETRA's strong adaptability to complex geometric boundaries, while avoiding the shortcomings of OCTREE in handling poor quality meshes of arched and circular structures and TETRA's difficulty in meeting the analysis requirements of simulation layering and seam boundaries. It can quickly generate fine mesh models and overcome the problems of low discretization efficiency and difficulty in meeting the engineering analysis requirements of traditional methods.
[0042] This invention can rapidly establish a high-quality, detailed model considering the entire system of foundation-cavity-dam, avoiding the shortcomings of traditional sub-model analysis or multi-point constraint analysis methods in considering the interaction between foundation and structure. Furthermore, the method is highly simple, operable, and universal, providing an effective approach for the efficient mesh generation of other similar complex spatial cavity systems.
[0043] It should be noted that, for the sake of simplicity, the foregoing embodiments are all described as a series of actions. However, those skilled in the art should understand that the present invention is not limited to the described order of actions, as some steps may be performed in other orders or simultaneously according to the present invention. Furthermore, those skilled in the art should also understand that the embodiments described in the specification are preferred embodiments, and the actions and modules involved are not necessarily essential to the present invention.
[0044] The above embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit the scope of protection of the invention. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on these embodiments, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art can still combine, add, delete, or otherwise adjust the features of the various embodiments of the present invention according to the circumstances without conflict or creative effort, thereby obtaining different technical solutions that do not fundamentally depart from the concept of the present invention. These technical solutions also fall within the scope of protection of the present invention.
Claims
1. A multi-source interface coupling transition modeling method considering spatially anisotropic arched chambers, characterized in that, The method includes: S10, Geometric Entity Creation and Region Structured Segmentation: Obtain engineering dimension parameters, create geometric entities based on the engineering dimension parameters, and perform region structured segmentation on the created geometric entities; S20, cross-scale OCTREE-TETRA coupled discretization: the main mesh and the cavern area mesh are created using OCTREE and hexahedral sweeping methods respectively. Through the conformal / non-conformal boundary fusion algorithm, the displacement coordination of the boundary nodes is forced to construct a stress-continuous overall calculation model. S30, Mesh Quality Diagnosis, Correction and Verification: Based on a pre-built multi-index mesh evaluation system, the fused mesh is subjected to quality detection and optimization to ensure that it meets the requirements of finite element analysis.
2. The multi-source interface coupling transition modeling method considering spatially anisotropic arched caverns as described in claim 1, characterized in that: The method for creating geometric entities and segmenting regions in a structured manner includes: S101, Obtain engineering dimension parameters, and create geometric entities based on engineering dimension parameters; When creating geometric entities based on engineering dimension parameters, an initial three-dimensional entity is generated by axial stretching based on the engineering dimension parameters; based on the upstream and downstream slope ratio of the dam body and the valley topography, the main structural entity model containing the clay core dam is constructed by positioning the working plane and Boolean cutting operations. S102, Perform region-structured segmentation on the created geometric entities; When performing regional structural segmentation on the created geometric entity, the geometric entity is divided into the main structural area and the spatially anisotropic arched cave area.
3. The multi-source interface coupling transition modeling method considering spatially anisotropic arched caverns as described in claim 2, characterized in that: The spatially oriented arched cavern area is divided into three structural regions based on its structural characteristics: a circular cross-section cavern body (region A), a high-plasticity clay filling layer (region B), and a geometric transition connection region (region C).
4. The multi-source interface coupling transition modeling method considering spatially anisotropic arched caverns as described in claim 3, characterized in that: The cross-scale OCTREE-TETRA coupled discretization method includes: S201, Load the created geometric entity, and perform cross-scale discretization of the main structure based on OCTREE; S202, the circular cross-section cavern body A area, the high plasticity clay filling layer B area, and the geometric transition connection area C area in the geometric entity are subjected to hexahedral sweep discretization; S203, obtains cross-scale discretization and hexahedral sweep discretization results, performs multi-source interface coupling transition discretization, and adaptively generates the transition region TETRA. S204 imports the Model-A, Tunnel-A, Tunnel-B, and Model-Transition-TETRA meshes into the mesh coupling system. Through a conformal / non-conformal boundary fusion algorithm, it forces the displacement coordination of the boundary nodes and constructs a stress-continuous overall calculation model.
5. The multi-source interface coupling transition modeling method considering spatially anisotropic arched caverns as described in claim 4, characterized in that: When performing cross-scale discretization of the main structure based on OCTREE, the created geometric entities are imported into the OCTREE mesh generator, and differentiated mesh sizes are set according to the mechanical response requirements of the upstream riprap zone, core wall zone, and downstream riprap zone to generate a cross-scale refined mesh Model-A, with the solid body of the central wall corridor marked as Model-B.
6. The multi-source interface coupling transition modeling method considering spatially anisotropic arched caverns as described in claim 5, characterized in that: When performing hexahedral sweep discretization on the circular cross-section cavern body A area, the high-plasticity clay filling layer B area, and the geometric transition connection area C area in the geometric entity, a two-dimensional planar mesh is generated based on the intermediate cross-section mesh parameters. The hexahedral sweep operation is performed along the cavern axis direction according to the set size, and the Tunnel-A three-dimensional mesh, Tunnel-B three-dimensional mesh, and Tunnel-C three-dimensional mesh corresponding to the circular cross-section cavern body A area, the high-plasticity clay filling layer B area, and the geometric transition connection area C area are output respectively.
7. The multi-source interface coupling transition modeling method considering spatially anisotropic arched caverns as described in claim 6, characterized in that: The method for discretizing the coupling transition of multi-source interfaces includes: S2031 automatically constructs conformal interface elements between Model-A and Model-B through a node topology mapping algorithm, extracts the core wall side node set and outputs it as a surface mesh file Interface-M-AB; S2032, similarly generate the Interface-M-BC file between Tunnel-B and Tunnel-C; S2033, a closed cavity surface is constructed based on Interface-M-AB and Interface-M-BC; S2034, execute the surface-to-solid transformation to generate the transition geometry Model-Transition; S2035 adopts the principle of one discrete part for the global line model, performs tetrahedral automatic meshing on the Model-Transition, and outputs the Model-Transition-TETRA mesh.
8. The multi-source interface coupling transition modeling method considering spatially anisotropic arched caverns as described in any one of claims 2-7, characterized in that: The mesh quality diagnosis, correction, and verification method includes: S301 performs Jacobian determinant, twist, and warpage detection on the fused mesh; S302 performs topology optimization and iterative adjustment of node positions for low-quality cells; S303, the mesh file was verified to meet the convergence requirements of finite element analysis through trial calculations.