Finite element-smooth point interpolation coupling simulation method suitable for coal mining rock stratum deformation movement
By using a coupled finite element method and smooth point interpolation, elements with deformation parameters greater than a threshold are marked and smooth point interpolation simulations are performed. This solves the mesh dependency problem of finite element and finite difference methods when simulating rock strata deformation in coal mining, and achieves a more stable and accurate simulation of rock strata deformation.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- CHINA UNIV OF GEOSCIENCES (BEIJING)
- Filing Date
- 2026-01-08
- Publication Date
- 2026-05-08
AI Technical Summary
Existing numerical simulation methods such as finite element method and finite difference method have mesh dependency problems when simulating rock strata movement and deformation induced by coal mining. They are difficult to take into account both continuous and discontinuous deformation, resulting in deviations between simulation results and actual conditions, and reducing the stability and accuracy of rock strata movement and deformation simulation.
A coupled simulation method of finite element method and smooth point interpolation is adopted. In the finite element method, elements with deformation parameters greater than the threshold are marked as elements to be converted, and smooth point interpolation is used to simulate them. At the same time, the finite element method and smooth point interpolation process are coupled to dynamically simulate the deformation of rock strata.
It improves the stability and accuracy of simulation results, effectively analyzes the deformation and movement patterns of overlying rock under the influence of mining, solves the grid dependency problem, and improves computational efficiency and simulation accuracy.
Smart Images
Figure CN121997646A_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of data simulation and processing technology, and in particular to a finite element-smooth point interpolation coupled simulation method suitable for the deformation and movement of coal mining strata. Background Technology
[0002] As the world's largest producer and consumer of coal, China relies heavily on coal mining for its national energy security, economic development, and industrial structure. However, coal mining causes severe disturbances to rock strata, inducing rock movement and surface subsidence. This inevitably poses a significant threat to buildings, transportation infrastructure, and water conservancy projects in mining areas, leading to prominent social, economic, and environmental problems. Therefore, accurately analyzing rock deformation and surface subsidence caused by mining is a crucial prerequisite for mining planning and design, and is of paramount importance for the safe production and construction of mines.
[0003] For a long time, scholars at home and abroad have paid close attention to the problem of rock strata movement and deformation induced by underground mining, and have mainly developed theoretical analysis methods, physical model test methods, and numerical simulation methods. Due to their strong adaptability, large computational scale, and high computational efficiency, numerical simulation methods such as finite element method (FEM) and finite difference method (FD) have become important tools for analyzing the mechanical behavior of rock strata induced by coal mining. However, FEM and FD methods are suitable for analyzing continuous problems, but cannot simulate discontinuous behaviors such as cracking and caving; at the same time, FEM and FD methods have mesh dependence, which can lead to mesh distortion when analyzing the above problems, causing the simulation results to deviate from reality, or even become uncalculated.
[0004] Therefore, related technologies employ equivalent mining methods to simulate mining. Taking the simulation of caving processes as an example, a common approach is to predetermine the extent of the caving zone based on empirical formulas, then assume the rock mass in the caving zone is an elasto-plastic or granular material, and determine the mechanical parameters of the rock mass according to experimental results or empirical formulas. This method, by adjusting the mechanical parameters of the rock mass in the failure zone, can achieve equivalent simulation of discontinuous behavior.
[0005] However, existing numerical simulation methods such as the finite element method and finite difference method suffer from mesh dependence when dealing with complex overburden deformation problems. Limited by mesh distortion, it is difficult to simultaneously account for continuous and discontinuous deformation, leading to discrepancies between simulation results and actual conditions. This makes it difficult to analyze the deformation and movement patterns of overburden under mining influence, reducing the stability and accuracy of rock strata movement and deformation simulations, which urgently needs to be addressed. Summary of the Invention
[0006] This application is based on the inventor's understanding and insights into the following issues: For a long time, scholars at home and abroad have paid close attention to the problem of rock strata movement and deformation induced by underground mining, and have mainly developed theoretical analysis methods, physical model test methods, and numerical simulation methods. Due to their strong adaptability, large computational scale, and high computational efficiency, numerical simulation methods such as finite element and finite difference have become important means to analyze the mechanical behavior of rock strata induced by coal mining.
[0007] However, methods such as the finite element method (FEM) and finite difference method (FDDM) still have limitations in simulating rock strata movement and deformation induced by coal mining. These limitations mainly stem from the fact that during coal mining, due to severe disturbances, rock strata can deform, crack, and even collapse and compact. These processes simultaneously involve continuous and discontinuous processes, and methods like the FEM and FDDM still face the following challenges when simulating them: 1) Finite element method, finite difference method and other methods are suitable for analyzing continuous problems, but cannot simulate discontinuous behaviors such as cracking and collapse.
[0008] 2) Finite element and finite difference methods exhibit mesh dependence, such as... Figure 1 As shown, mesh distortion occurs when analyzing the above problems, causing the simulation results to deviate from reality, or even become uncalculateable.
[0009] To address the aforementioned issues, researchers typically employ equivalent mining simulation methods. Taking the simulation of caving as an example, a common approach is to predetermine the extent of the caving zone using empirical formulas, then assume the rock mass within the caving zone is an elasto-plastic or granular material, and determine the rock mass's mechanical parameters based on experimental results or empirical formulas. This method, by adjusting the rock mass's mechanical parameters in the failure zone, can achieve an equivalent simulation of discontinuous behavior. However, the equivalent simulation method still suffers from the following problems: 1) The extent of the failure zone (fracture zone, caving zone) is predetermined, which does not reflect the actual situation; 2) When encountering large deformations, this type of method still exhibits grid dependence, making calculations difficult and requiring urgent improvement.
[0010] This application provides a finite element-smooth point interpolation coupled simulation method suitable for coal mining strata deformation and movement. This addresses the mesh dependency problem inherent in existing numerical simulation methods such as the finite element method and finite difference method when dealing with complex overburden deformation. Limited by mesh distortion, it is difficult to simultaneously account for continuous and discontinuous deformation, leading to deviations between simulation results and actual conditions. This makes it difficult to analyze the deformation and movement patterns of overburden under mining influence, reducing the stability and accuracy of strata movement and deformation simulations.
[0011] The first aspect of this application provides a finite element-smooth point interpolation coupled simulation method suitable for coal mining strata deformation and movement, comprising the following steps: First, based on exploration data, a three-dimensional calculation model of the target study area is established, and the calculation parameters of the strata are obtained. Then, the continuous deformation of the strata is simulated using the finite element method, and the deformation parameters corresponding to the target elements in the strata are determined. Second, elements in the target elements whose deformation parameters are greater than or equal to a preset deformation threshold are marked as elements to be converted, and the smooth point interpolation method is used to continue the smooth point interpolation simulation of the regions to be converted corresponding to the elements to be converted. Simultaneously, the finite element method is used to perform the finite element simulation on the regions corresponding to elements whose deformation parameters are less than the preset deformation threshold. The smooth point interpolation simulation process and the finite element simulation process are coupled to dynamically simulate the deformation of the strata in the target study area based on the finite element-smooth point interpolation coupled simulation method, so as to output the dynamic simulation results of the strata deformation in the target study area.
[0012] Optionally, in one embodiment of this application, before simulating the continuous deformation of the rock strata in the target study area using the finite element method, the method further includes: establishing a three-dimensional geological generalization model of the target study area based on engineering geological data, and subdividing the three-dimensional geological generalization model to obtain an initial three-dimensional calculation model; applying displacement boundary conditions to the initial three-dimensional calculation model, and solving the initial stress field and displacement field of the initial three-dimensional calculation model to construct the three-dimensional calculation model that meets preset conditions.
[0013] Optionally, in one embodiment of this application, before continuing to perform smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted using the smooth point interpolation method, the method further includes: obtaining the centroid coordinates of all units to be converted, and establishing a minimum bounding box containing the centroid coordinates of all units to be converted; determining multiple tetrahedrons corresponding to the centroid coordinates of all units to be converted in the minimum bounding box; identifying the multiple tetrahedrons as the final units to be converted, and using the final units to be converted to determine the region to be converted.
[0014] Optionally, in one embodiment of this application, the step of using a smooth point interpolation method to continue the smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted includes: obtaining the centroid of the tetrahedron containing each triangular face in the target tetrahedral mesh model, connecting the centroid sequentially to the three vertices of the corresponding triangular face to construct a smooth domain for each triangular face; and based on the smooth domain of each triangular face, continuing the smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted.
[0015] Optionally, in one embodiment of this application, the dynamic simulation of the deformation of the rock strata in the target study area based on the finite element-smooth point interpolation coupled simulation method includes: identifying the coupling interface between the finite element simulation region and the smooth point interpolation simulation region in the target study area; realizing the transmission of coupling information between the finite element simulation region and the smooth point interpolation simulation region; and dynamically simulating the deformation of the rock strata in the target study area based on the coupling interface and the transmission of coupling information.
[0016] A second aspect of this application provides a finite element-smooth point interpolation coupled simulation device suitable for coal mining strata deformation and movement, comprising: a finite element simulation module, used to input strata parameter information of a target study area into a pre-constructed three-dimensional calculation model of the target study area, simulate the continuous deformation of the strata using the finite element method, and determine the deformation parameters corresponding to target units in the strata; a smooth point interpolation simulation module, used to mark units in the target units whose deformation parameters are greater than or equal to a preset deformation threshold as units to be converted, and use the smooth point interpolation method to continue smooth point interpolation simulation on the regions to be converted corresponding to the units to be converted; and a coupling module, used to use the finite element method to perform the finite element simulation on the regions corresponding to units whose deformation parameters are less than the preset deformation threshold, and to couple the smooth point interpolation simulation process with the finite element simulation process, so as to dynamically simulate the deformation of the strata in the target study area based on the finite element-smooth point interpolation coupled simulation, and output the dynamic simulation results of the strata deformation in the target study area.
[0017] Optionally, in one embodiment of this application, the apparatus further includes: a model building module, used to establish a three-dimensional geological generalization model of the target study area based on engineering geological data before simulating the continuous deformation of the rock strata in the target study area using the finite element method, and to divide the three-dimensional geological generalization model to obtain an initial three-dimensional calculation model; and a calculation initialization module, used to apply displacement boundary conditions to the initial three-dimensional calculation model before simulating the continuous deformation of the rock strata in the target study area using the finite element method, and to solve the initial stress field and displacement field of the initial three-dimensional calculation model to construct the three-dimensional calculation model that meets preset conditions.
[0018] Optionally, in one embodiment of this application, the apparatus further includes: a bounding box search module, configured to obtain the centroid coordinates of all units to be converted and establish a minimum bounding box containing the centroid coordinates of all units to be converted before continuing to perform smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted using the smooth point interpolation method; a tetrahedron search module, configured to determine a plurality of tetrahedrons corresponding to the centroid coordinates of all units to be converted in the minimum bounding box before continuing to perform smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted using the smooth point interpolation method; and a conversion module, configured to identify the plurality of tetrahedrons as the final units to be converted and use the final units to be converted to determine the region to be converted before continuing to perform smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted using the smooth point interpolation method.
[0019] Optionally, in one embodiment of this application, the smooth point interpolation simulation module includes: a smooth domain construction unit, used to obtain the centroid of the tetrahedron in which each triangular face is located in the target tetrahedral mesh model, and connect the centroid to the three vertices of the corresponding triangular face in sequence to construct a smooth domain for each triangular face; and a smooth point interpolation simulation unit, used to continue performing the smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted based on the smooth domain of each triangular face.
[0020] Optionally, in one embodiment of this application, the coupling module includes: a coupling interface identification unit for identifying the coupling interface between the finite element simulation region and the smooth point interpolation simulation region in the target study area; a coupling information transmission unit for transmitting coupling information between the finite element simulation region and the smooth point interpolation simulation region; and a dynamic simulation unit for dynamically simulating the deformation of the rock strata in the target study area based on the coupling interface and the coupling information transmission.
[0021] A third aspect of this application 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 finite element-smooth point interpolation coupled simulation method for deformation and movement of coal strata as described in the above embodiments.
[0022] A fourth aspect of this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described finite element-smooth point interpolation coupled simulation method suitable for coal strata deformation and movement.
[0023] A fifth aspect of this application provides a computer program product, including a computer program that, when executed, is used to implement the above-described finite element-smooth point interpolation coupled simulation method suitable for coal mining strata deformation and movement.
[0024] This application embodiment can input the rock strata parameter information of the study area into a three-dimensional calculation model. First, the continuous deformation of the rock strata is simulated using the finite element method, and the elements in the rock strata whose deformation parameters are greater than or equal to the deformation threshold are marked as elements to be converted. Then, the smooth point interpolation method is used to continue the smooth point interpolation simulation of the regions to be converted corresponding to the elements to be converted. Next, the finite element method is used to perform finite element simulation on the regions corresponding to the elements whose deformation parameters are less than the deformation threshold, and the smooth point interpolation simulation process and the finite element simulation process are coupled to output the dynamic simulation results of the rock strata deformation in the study area. Thus, based on the finite element-smooth point interpolation coupled simulation method, the deformation and movement law of the overburden under the influence of mining can be analyzed, and the deformation and movement process of the overburden can be accurately simulated. This solves the problem in related technologies where existing numerical simulation methods such as finite element and finite difference have mesh dependency problems when dealing with complex deformation problems of overburden. Due to mesh distortion, it is difficult to take into account both continuous and discontinuous deformation, resulting in deviations between the simulation results and the actual situation, making it difficult to analyze the deformation and movement law of the overburden under the influence of mining, and reducing the stability and accuracy of the rock strata movement and deformation simulation.
[0025] Additional aspects and advantages of this application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this application. Attached Figure Description
[0026] The above and / or additional aspects and advantages of this application will become apparent and readily understood from the following description of the embodiments taken in conjunction with the accompanying drawings, wherein: Figure 1 This is a schematic diagram of grid distortion in coal mining simulation using related technologies; Figure 2 This is a flowchart of a finite element-smooth point interpolation coupled simulation method suitable for the deformation and movement of coal mining strata, according to an embodiment of this application. Figure 3 A schematic diagram of a three-dimensional calculation model of a goaf area for a specific implementation of this application; Figure 4 This is a schematic diagram of the vertical displacement of the rock strata after mining 50m in a specific implementation of this application; Figure 5 This is a schematic diagram of the equivalent plastic strain of the rock strata after mining 50m in a specific implementation of this application; Figure 6 This is a schematic diagram of the simulated region division for a specific implementation of this application; Figure 7This is a schematic diagram of the internal and boundary smooth regions of a specific embodiment of this application; Figure 8 This is a schematic diagram of the simulation region coupling in a specific embodiment of this application; Figure 9 This is a schematic diagram of a finite element-smooth point interpolation coupled simulation device suitable for coal mining strata deformation and movement, according to an embodiment of this application. Figure 10 This is a schematic diagram of the structure of an electronic device provided according to an embodiment of this application. Detailed Implementation
[0027] The embodiments of this application are described in detail below. Examples of these embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and intended to explain this application, and should not be construed as limiting this application.
[0028] The following describes, with reference to the accompanying drawings, a finite element-smooth point interpolation coupled simulation method suitable for coal mining strata deformation and movement, according to embodiments of this application. Regarding the related technologies mentioned in the background section, existing numerical simulation methods such as the finite element method and finite difference method suffer from mesh dependency issues when dealing with complex overburden deformation problems. Limited by mesh distortion, it is difficult to simultaneously account for continuous and discontinuous deformation, leading to deviations between simulation results and actual conditions. This makes it difficult to analyze the deformation and movement patterns of overburden under mining influence, reducing the stability and accuracy of rock strata movement and deformation simulation. This application provides a finite element-smooth point interpolation coupled simulation method suitable for coal mining rock strata deformation and movement. In this method, the rock strata parameter information of the study area can be input into a three-dimensional calculation model. First, the continuous deformation of the rock strata is simulated using the finite element method, and the elements in the rock strata with deformation parameters greater than or equal to the deformation threshold are marked as elements to be converted. Then, the smooth point interpolation method is used to continue the smooth point interpolation simulation of the regions to be converted corresponding to the elements to be converted. Next, the finite element method is used to perform finite element simulation on the regions corresponding to the elements with deformation parameters less than the deformation threshold, and the smooth point interpolation simulation process and the finite element simulation process are coupled to output the dynamic simulation results of rock strata deformation in the study area. Thus, based on the finite element-smooth point interpolation coupled simulation method, the deformation and movement patterns of overburden under mining influence can be analyzed, and the deformation and movement process of overburden can be accurately simulated. This solves the problem of mesh dependency in existing numerical simulation methods such as finite element and finite difference methods when dealing with complex overburden deformation. Limited by mesh distortion, it is difficult to simultaneously account for continuous and discontinuous deformation, leading to deviations between simulation results and actual conditions. This makes it difficult to analyze the deformation and movement patterns of overburden under mining influence, reducing the stability and accuracy of rock strata movement and deformation simulations.
[0029] Specifically, Figure 2 This is a schematic flowchart illustrating a finite element-smooth point interpolation coupled simulation method suitable for the deformation and movement of coal mining strata, provided in an embodiment of this application.
[0030] like Figure 2 As shown, the finite element-smooth point interpolation coupled simulation method suitable for coal mining strata deformation and movement includes the following steps: In step S201, the rock strata parameter information of the target study area is input into the pre-constructed three-dimensional calculation model of the target study area, the continuous deformation of the rock strata is simulated using the finite element method, and the deformation parameters corresponding to the target element in the rock strata are determined.
[0031] In the embodiments of this application, the target study area is the area where overburden deformation simulation analysis needs to be performed, such as a coal mine; the target unit is the deformation unit in the rock strata.
[0032] It is understood that, in the embodiments of this application, rock strata parameter information of the study area, such as mechanical parameters, can be input into the three-dimensional calculation model pre-constructed in the following steps, and the continuous deformation of the rock strata can be simulated using the finite element method to determine the deformation parameters corresponding to the target element in the rock strata.
[0033] In the initial stage of mining, the finite element method is used to simulate the continuous deformation of the rock strata, and the equivalent mining method is used to simulate mining. Specifically, for the area to be mined, the rock mass in the area is assumed to be an elastoplastic material, and the mechanical parameters of the rock mass are determined according to experimental results or empirical formulas. In the equivalent mining simulation process, in order to better reflect the crushing process of the broken rock mass in the goaf and caving zone, the double yield model is used for calculation, while for the remaining parts, the Mohr-Coulomb constitutive model is used for calculation.
[0034] In addition, in order to simulate the failure process of rock strata such as cracking and collapse, it is necessary to further reduce the rock mass mechanical parameters to better simulate discontinuous mechanical behavior using continuous media. According to relevant research, when the calculation parameters are reduced to the mechanical parameters after the rock mass is subjected to compressive failure, the simulation results can be more consistent with the actual situation. Among them, it is usually appropriate to reduce the rock mass elastic modulus to 16.7% of the original parameter.
[0035] Optionally, in one embodiment of this application, before simulating the continuous deformation of the rock strata in the target study area using the finite element method, the method further includes: establishing a three-dimensional geological generalization model of the target study area based on engineering geological data, and subdividing the three-dimensional geological generalization model to obtain an initial three-dimensional calculation model; applying displacement boundary conditions to the initial three-dimensional calculation model, and solving the initial stress field and displacement field of the initial three-dimensional calculation model to construct a three-dimensional calculation model that meets preset conditions.
[0036] In actual implementation, the embodiments of this application can establish a three-dimensional geological generalization model of the study area and perform meshing on the three-dimensional geological generalization model. Specifically, the embodiments of this application can establish a three-dimensional geological generalization model of the study area based on engineering geological data and use tetrahedral meshing to obtain an initial three-dimensional calculation model. The initial three-dimensional calculation model has a length of 1200m, a height of 255m, and a thickness of 10m. After tetrahedral meshing, the number of nodes is 6889 and the number of elements is 20771.
[0037] Next, in this embodiment, displacement boundary conditions can be applied to the four sides of the initial three-dimensional calculation model. Specifically, the displacement at the bottom of the initial three-dimensional calculation model is fixed at 0, the normal displacement around the model is fixed at 0, and the top is free. Then, the mechanical parameters of the rock strata are obtained, and the initial stress field and displacement field of the initial three-dimensional calculation model are solved using an elastic constitutive method. The displacement field is then zeroed out to reduce errors, thereby determining... Figure 3 The three-dimensional computational model shown is designed to improve computational efficiency.
[0038] In step S202, the elements in the target element whose deformation parameters are greater than or equal to the preset deformation threshold are marked as elements to be converted, and the smooth point interpolation method is used to continue to perform smooth point interpolation simulation on the region to be converted corresponding to the elements to be converted.
[0039] It is understood that, in the initial stage of mining, the finite element method described above can be used for simulation, such as... Figure 4 As shown, after a certain mining step is completed, the degree of deformation of each unit in the rock strata is calculated, in order to... Figure 5 Taking equivalent plastic strain as an example, when the equivalent plastic strain of a certain element exceeds a set threshold, such as when the equivalent plastic strain of a certain element is not 0, that is, when the element undergoes plastic deformation, and when the equivalent plastic strain exceeds the set threshold (e.g., greater than or equal to 5%), the element is marked as an element to be converted. Then, as... Figure 6 As shown, the smooth point interpolation method can be used to continue the smooth point interpolation simulation of the elements in the region to be converted in the following steps corresponding to the element to be converted. The calculation information can be directly inherited from the tetrahedral background mesh in the region to be converted, which effectively improves the calculation efficiency of the simulation.
[0040] Optionally, in one embodiment of this application, before continuing to perform smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted using the smooth point interpolation method, the method further includes: obtaining the centroid coordinates of all units to be converted and establishing a minimum bounding box containing the centroid coordinates of all units to be converted; determining multiple tetrahedrons corresponding to the centroid coordinates of all units to be converted in the minimum bounding box; identifying the multiple tetrahedrons as the final units to be converted and using the final units to be converted to determine the region to be converted.
[0041] As one possible implementation, this embodiment of the application can determine the region to be converted after the excavation calculations in the above steps are completed. Specifically, based on the centroid coordinates of all units to be converted, a hexahedron with the smallest volume that can contain all centroids is established, i.e., the minimum bounding box. Then, all tetrahedrons whose centroids are within the bounding box are identified as units to be converted, thereby determining the region to be converted. It should be noted that the determination of the region to be converted is dynamic, not predetermined. For other regions, i.e., regions not to be converted, such as... Figure 6 In the finite element simulation region, the finite element method is still used for simulation, which effectively improves the simulation accuracy and computational efficiency.
[0042] Optionally, in one embodiment of this application, a smooth point interpolation method is used to continue smooth point interpolation simulation of the region to be converted corresponding to the unit to be converted, including: obtaining the centroid of the tetrahedron in which each triangular face is located in the target tetrahedral mesh model, connecting the centroid to the three vertices of the corresponding triangular face in sequence to construct a smooth domain for each triangular face; and based on the smooth domain of each triangular face, continuing smooth point interpolation simulation of the region to be converted corresponding to the unit to be converted.
[0043] In some embodiments, this application employs a smooth point interpolation method to simulate the region to be transformed. Unlike the finite element method, the smooth point interpolation method requires the construction of a smooth domain, which can be established based on a tetrahedral background mesh within the region to be transformed. Specifically, based on each triangular face in the tetrahedral mesh model, the centroid of the tetrahedron containing the triangular face is taken, and the centroid is sequentially connected to the three vertices of the triangular face, thus forming the smooth domain to which that triangular face belongs. For example... Figure 7 As shown in (a), if the triangular face is not on the boundary, then the smooth region to which this triangular face belongs is a hexahedron; Figure 7 As shown in (b), when the triangular face is on the boundary, the smooth domain of this triangular face is a tetrahedron. Therefore, based on the smooth domain of each triangular face, the region to be transformed corresponding to the unit to be transformed can continue to be simulated by smooth point interpolation, which effectively improves the simulation calculation efficiency.
[0044] In step S203, the finite element method is used to perform finite element simulation on the region corresponding to the element whose deformation parameter is less than the preset deformation threshold. The smooth point interpolation simulation process is coupled with the finite element simulation process. The deformation of the rock strata in the target study area is dynamically simulated based on the finite element-smooth point interpolation coupled simulation method, so as to output the dynamic simulation results of the rock strata deformation in the target study area.
[0045] It is understood that the embodiments of this application can utilize the finite element method to continue finite element simulation of the region corresponding to the element whose deformation parameter is less than a certain deformation threshold (e.g., less than 5%), and couple the smooth point interpolation simulation process with the finite element simulation process. Thus, the embodiments of this application can simulate the continuous deformation, cracking, caving and compaction process of the overburden in the target study area based on the finite element-smooth point interpolation coupled simulation method. This helps to analyze the deformation and movement law of the overburden under the influence of mining and provides technical support for ensuring the engineering construction of the mining area.
[0046] Optionally, in one embodiment of this application, the deformation of the rock strata in the target study area is dynamically simulated based on the finite element-smooth point interpolation coupled simulation method, including: identifying the coupling interface between the finite element simulation region and the smooth point interpolation simulation region in the target study area; realizing the coupling information transmission between the finite element simulation region and the smooth point interpolation simulation region; and dynamically simulating the deformation of the rock strata in the target study area based on the coupling interface and the coupling information transmission.
[0047] As one possible implementation, the computational model in this embodiment is divided into a smooth point interpolation simulation region and a finite element simulation region. In this case, since the shape functions of both the finite element method and the smooth point interpolation method possess the Kronecker delta property, they can be directly coupled without requiring additional techniques. It should be noted that, as... Figure 8 As shown, for the corresponding coupled node, its relevant attribute values need to be inherited from both regions.
[0048] In summary, the embodiments of this application can first use the finite element method to simulate the initial stage of mining in the target study area, and constantly monitor the deformation of the elements in the rock strata. When the deformation of the elements reaches a certain threshold, the elements exceeding the threshold are automatically marked, and the marked elements are simulated using the smooth point interpolation method. Other areas are still simulated using the finite element method. Then, the smooth point interpolation simulation and the finite element simulation process are coupled to achieve dynamic simulation. Compared with related technologies, this application has good stability, no mesh dependency, and high computational efficiency. It can effectively simulate the deformation and movement process of the overburden, which helps to analyze the deformation and movement law of the overburden under the influence of mining, and provides technical support for ensuring the engineering construction of the mining area.
[0049] According to the finite element-smooth point interpolation coupled simulation method for deformation and movement of coal mining strata proposed in this application, the strata parameter information of the study area can be input into a three-dimensional calculation model. First, the continuous deformation of the strata is simulated using the finite element method, and elements with deformation parameters greater than or equal to the deformation threshold are marked as elements to be converted. Then, the smooth point interpolation method is used to continue the smooth point interpolation simulation of the regions corresponding to the elements to be converted. Next, the finite element method is used to simulate the regions corresponding to elements with deformation parameters less than the deformation threshold, and the smooth point interpolation simulation process and the finite element simulation process are coupled to output the dynamic simulation results of the strata deformation in the study area. Therefore, based on the finite element-smooth point interpolation coupled simulation method, the deformation and movement law of overburden under the influence of mining can be analyzed, and the deformation and movement process of overburden can be accurately simulated. This solves the problem in related technologies where existing numerical simulation methods such as finite element and finite difference methods have mesh dependency issues when dealing with complex overburden deformation problems. Limited by mesh distortion, it is difficult to simultaneously consider continuous and discontinuous deformation, leading to deviations between simulation results and actual conditions, and making it difficult to analyze the deformation and movement law of overburden under the influence of mining.
[0050] Secondly, referring to the accompanying drawings, a finite element-smooth point interpolation coupled simulation device suitable for coal mining strata deformation and movement, based on an embodiment of this application, is described.
[0051] Figure 9 This is a block diagram of a finite element-smooth point interpolation coupled simulation device suitable for coal mining strata deformation and movement, according to an embodiment of this application.
[0052] like Figure 9 As shown, the finite element-smooth point interpolation coupled simulation device 10 suitable for coal mining strata deformation and movement includes: a finite element simulation module 100, a smooth point interpolation simulation module 200, and a coupling module 300.
[0053] Specifically, the finite element simulation module 100 is used to input the rock strata parameter information of the target study area into the pre-constructed three-dimensional calculation model of the target study area, simulate the continuous deformation of the rock strata using the finite element method, and determine the deformation parameters corresponding to the target elements in the rock strata.
[0054] The smooth point interpolation simulation module 200 is used to mark the elements in the target element whose deformation parameters are greater than or equal to the preset deformation threshold as elements to be converted, and to continue to perform smooth point interpolation simulation on the region to be converted corresponding to the elements to be converted using the smooth point interpolation method.
[0055] The coupling module 300 is used to perform finite element simulation on the region corresponding to the element whose deformation parameter is less than the preset deformation threshold using the finite element method, and to couple the smooth point interpolation simulation process with the finite element simulation process. In order to perform dynamic simulation of the deformation of the rock strata in the target study area based on the finite element-smooth point interpolation coupled simulation, the dynamic simulation results of the rock strata deformation in the target study area are output.
[0056] Optionally, in one embodiment of this application, the apparatus 10 of this application embodiment further includes: a model building module and a calculation initialization module.
[0057] The model building module is used to establish a three-dimensional geological generalization model of the target study area based on engineering geological data before simulating the continuous deformation of the rock strata in the target study area using the finite element method, and to divide the three-dimensional geological generalization model to obtain an initial three-dimensional calculation model. The calculation initialization module is used to apply displacement boundary conditions to the initial three-dimensional calculation model before simulating the continuous deformation of the rock strata in the target study area using the finite element method, and to solve the initial stress field and displacement field of the initial three-dimensional calculation model in order to construct a three-dimensional calculation model that meets the preset conditions.
[0058] Optionally, in one embodiment of this application, the apparatus 10 of this application embodiment further includes: a search bounding box module, a search tetrahedron module, and a conversion module.
[0059] The bounding box search module is used to obtain the centroid coordinates of all units to be transformed and establish the minimum bounding box containing the centroid coordinates of all units to be transformed before continuing to perform smooth point interpolation simulation on the region to be transformed corresponding to the unit to be transformed using the smooth point interpolation method.
[0060] The tetrahedron search module is used to determine multiple tetrahedrons corresponding to the centroid coordinates of all units to be transformed within the minimum bounding box before continuing the smooth point interpolation simulation of the region to be transformed corresponding to the unit to be transformed using the smooth point interpolation method.
[0061] The conversion module is used to identify multiple tetrahedrons as the final units to be converted before continuing to perform smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted using the smooth point interpolation method, and to determine the region to be converted using the final units to be converted.
[0062] Optionally, in one embodiment of this application, the smooth point interpolation simulation module 200 includes: a smooth domain construction unit and a smooth point interpolation simulation unit.
[0063] The smooth domain construction unit is used to obtain the centroid of the tetrahedron containing each triangular face in the target tetrahedral mesh model, and connects the centroid to the three vertices of the corresponding triangular face in sequence to construct the smooth domain of each triangular face.
[0064] The smooth point interpolation simulation unit is used to continue performing smooth point interpolation simulation on the region to be transformed corresponding to the unit to be transformed, based on the smooth domain of each triangular face.
[0065] Optionally, in one embodiment of this application, the coupling module 300 includes: a coupling interface recognition unit, a coupling information transmission unit, and a dynamic simulation unit.
[0066] Among them, the coupling interface recognition unit is used to identify the coupling interface between the finite element simulation region and the smooth point interpolation simulation region in the target study area.
[0067] The coupling information transmission unit is used to realize the coupling information transmission between the finite element simulation region and the smooth point interpolation simulation region.
[0068] The dynamic simulation unit is used to dynamically simulate the deformation of rock strata in the target study area based on the coupling interface and coupling information transmission.
[0069] It should be noted that the foregoing explanation of the embodiment of the finite element-smooth point interpolation coupling simulation method suitable for the deformation and movement of coal mining strata also applies to the finite element-smooth point interpolation coupling simulation device suitable for the deformation and movement of coal mining strata in this embodiment, and will not be repeated here.
[0070] According to the embodiment of this application, the finite element method-smooth point interpolation coupled simulation device suitable for coal mining strata deformation and movement can input the strata parameter information of the study area into a three-dimensional calculation model. First, the finite element method is used to simulate the continuous deformation of the strata, and elements with deformation parameters greater than or equal to the deformation threshold are marked as elements to be converted. Then, the smooth point interpolation method is used to continue the smooth point interpolation simulation of the regions corresponding to the elements to be converted. Next, the finite element method is used to simulate the regions corresponding to the elements with deformation parameters less than the deformation threshold, and the smooth point interpolation simulation process and the finite element simulation process are coupled to output the dynamic simulation results of the strata deformation in the study area. Therefore, based on the finite element-smooth point interpolation coupled simulation, the deformation and movement law of the overburden under the influence of mining can be analyzed, and the deformation and movement process of the overburden can be accurately simulated. This solves the problem in related technologies where existing numerical simulation methods such as finite element and finite difference methods have mesh dependency issues when dealing with complex overburden deformation problems. Limited by mesh distortion, it is difficult to simultaneously consider continuous and discontinuous deformation, leading to deviations between simulation results and actual conditions, and making it difficult to analyze the deformation and movement law of the overburden under the influence of mining.
[0071] Figure 10 A schematic diagram of the structure of an electronic device provided in an embodiment of this application. The electronic device may include: The memory 1001, the processor 1002, and the computer program stored on the memory 1001 and capable of running on the processor 1002.
[0072] When the processor 1002 executes the program, it implements the finite element-smooth point interpolation coupled simulation method for deformation and movement of coal mining strata provided in the above embodiments.
[0073] Furthermore, electronic devices also include: Communication interface 1003 is used for communication between memory 1001 and processor 1002.
[0074] The memory 1001 is used to store computer programs that can run on the processor 1002.
[0075] The memory 1001 may include high-speed RAM memory, and may also include non-volatile memory, such as at least one disk storage device.
[0076] If the memory 1001, processor 1002, and communication interface 1003 are implemented independently, then the communication interface 1003, memory 1001, and processor 1002 can be interconnected via a bus to complete communication between them. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. Buses can be categorized into address buses, data buses, control buses, etc. For ease of representation, Figure 10 The bus is represented by a single thick line, but this does not mean that there is only one bus or one type of bus.
[0077] Optionally, in a specific implementation, if the memory 1001, processor 1002, and communication interface 1003 are integrated on a single chip, then the memory 1001, processor 1002, and communication interface 1003 can communicate with each other through an internal interface.
[0078] The processor 1002 may be a central processing unit (CPU), an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of this application.
[0079] This embodiment also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described finite element-smooth point interpolation coupled simulation method suitable for the deformation and movement of coal mining strata.
[0080] This embodiment also provides a computer program product, including a computer program that, when executed, is used to implement the above-described finite element-smooth point interpolation coupled simulation method suitable for the deformation and movement of coal mining strata.
[0081] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.
[0082] Furthermore, 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 at least one of that feature. In the description of this application, "N" means at least two, such as two, three, etc., unless otherwise explicitly specified.
[0083] Any process or method described in the flowchart or otherwise herein can be understood as representing a module, segment, or portion of code comprising one or N executable instructions for implementing custom logic functions or processes, and the scope of the preferred embodiments of this application includes additional implementations in which functions may be performed not in the order shown or discussed, including substantially simultaneously or in reverse order depending on the functions involved, as should be understood by those skilled in the art to which embodiments of this application pertain.
[0084] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (such as a computer-based system, a processor-included system, or other system that can fetch and execute instructions from, an instruction execution system, apparatus, or device). For the purposes of this specification, "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transmit programs for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include: an electrical connection having one or more wires (electronic device), a portable computer disk drive (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and editable read-only memory (EPROM or flash memory), fiber optic devices, and portable optical disc read-only memory (CDROM). Alternatively, the computer-readable medium may be paper or other suitable media on which the program can be printed, since the program can be obtained electronically by optically scanning the paper or other medium, followed by editing, interpreting, or otherwise processing as necessary, and then stored in a computer memory.
[0085] It should be understood that the various parts of this application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiments, the N steps or methods can be implemented using software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, it can be implemented using any one or a combination of the following techniques known in the art: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.
[0086] Those skilled in the art will understand that all or part of the steps of the methods in the above embodiments can be implemented by a program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, the program includes one or a combination of the steps of the method embodiments.
[0087] Furthermore, the functional units in the various embodiments of this application can be integrated into a processing module, or each unit can exist physically separately, or two or more units can be integrated into a module. The integrated module can be implemented in hardware or as a software functional module. If the integrated module is implemented as a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium.
[0088] The storage medium mentioned above can be a read-only memory, a disk, or an optical disk, etc. Although embodiments of this application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting this application. Those skilled in the art can make changes, modifications, substitutions, and variations to the above embodiments within the scope of this application.
Claims
1. A finite element-smooth point interpolation coupled simulation method suitable for coal mining strata deformation and movement, characterized in that, Includes the following steps: The rock strata parameter information of the target study area is input into the pre-constructed three-dimensional calculation model of the target study area. The continuous deformation of the rock strata is simulated using the finite element method, and the deformation parameters corresponding to the target element in the rock strata are determined. The units in the target unit whose deformation parameters are greater than or equal to a preset deformation threshold are marked as units to be converted, and the smooth point interpolation method is used to continue to perform smooth point interpolation simulation on the region to be converted corresponding to the units to be converted. The finite element method is used to perform finite element simulation on the region corresponding to the unit whose deformation parameter is less than the preset deformation threshold. The smooth point interpolation simulation process is coupled with the finite element simulation process to dynamically simulate the deformation of the rock strata in the target study area based on the finite element-smooth point interpolation coupled simulation method, so as to output the dynamic simulation results of the rock strata deformation in the target study area.
2. The method according to claim 1, characterized in that, Before simulating the continuous deformation of the rock strata in the target study area using the finite element method, the following steps are also included: Based on engineering geological data, a three-dimensional geological generalization model of the target study area is established, and the three-dimensional geological generalization model is divided to obtain an initial three-dimensional calculation model. Displacement boundary conditions are applied to the initial three-dimensional calculation model, and the initial stress field and displacement field of the initial three-dimensional calculation model are solved to construct the three-dimensional calculation model that satisfies the preset conditions.
3. The method according to claim 1, characterized in that, Before performing smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted using the smooth point interpolation method, the following steps are also included: Obtain the centroid coordinates of all units to be transformed, and construct the minimum bounding box containing the centroid coordinates of all units to be transformed; Determine the multiple tetrahedrons corresponding to the centroid coordinates of all the units to be transformed in the minimum bounding box; The plurality of tetrahedrons are identified as the final units to be converted, and the regions to be converted are determined using the final units to be converted.
4. The method according to claim 3, characterized in that, The step of using a smooth point interpolation method to further simulate the smooth point interpolation of the region to be converted corresponding to the unit to be converted includes: Obtain the centroid of the tetrahedron containing each triangular face in the target tetrahedral mesh model, and connect the centroid to the three vertices of the corresponding triangular face in sequence to construct a smooth domain for each triangular face; Based on the smooth region of each triangular face, the smooth point interpolation simulation is continued for the region to be converted corresponding to the unit to be converted.
5. The method according to claim 1, characterized in that, The method based on finite element-smooth point interpolation coupling simulation is used to dynamically simulate the deformation of the rock strata in the target study area, including: Identify the coupling interface between the finite element simulation region and the smooth point interpolation simulation region in the target study area; To achieve coupled information transfer between the finite element simulation region and the smooth point interpolation simulation region; Based on the coupling interface and the coupling information transmission, the deformation of the rock strata in the target study area is dynamically simulated.
6. A finite element-smooth point interpolation coupled simulation device suitable for coal mining strata deformation and movement, characterized in that, include: The finite element simulation module is used to input the rock strata parameter information of the target study area into the pre-constructed three-dimensional calculation model of the target study area, simulate the continuous deformation of the rock strata using the finite element method, and determine the deformation parameters corresponding to the target element in the rock strata. The smooth point interpolation simulation module is used to mark the units in the target unit whose deformation parameters are greater than or equal to a preset deformation threshold as units to be converted, and to continue to perform smooth point interpolation simulation on the regions to be converted corresponding to the units to be converted using the smooth point interpolation method. The coupling module is used to perform finite element simulation on the region corresponding to the unit whose deformation parameter is less than the preset deformation threshold using the finite element method, and to couple the smooth point interpolation simulation process with the finite element simulation process, so as to dynamically simulate the deformation of the rock strata in the target study area based on the finite element-smooth point interpolation coupled simulation, and output the dynamic simulation results of the rock strata deformation in the target study area.
7. The apparatus according to claim 6, characterized in that, Also includes: The model building module is used to establish a three-dimensional geological generalization model of the target study area based on engineering geological data before simulating the continuous deformation of the rock strata in the target study area using the finite element method, and to divide the three-dimensional geological generalization model to obtain an initial three-dimensional calculation model. The calculation initialization module is used to apply displacement boundary conditions to the initial three-dimensional calculation model before simulating the continuous deformation of the rock strata in the target study area using the finite element method, and to solve the initial stress field and displacement field of the initial three-dimensional calculation model in order to construct the three-dimensional calculation model that meets the preset conditions.
8. The apparatus according to claim 6, characterized in that, Also includes: The bounding box search module is used to obtain the centroid coordinates of all units to be converted and establish the minimum bounding box containing the centroid coordinates of all units to be converted before continuing to perform smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted using the smooth point interpolation method. The tetrahedron search module is used to determine multiple tetrahedrons corresponding to the centroid coordinates of all units to be converted in the minimum bounding box before continuing to perform smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted using the smooth point interpolation method. The conversion module is used to identify the plurality of tetrahedrons as the final conversion units before continuing to perform smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted using the smooth point interpolation method, and to determine the region to be converted using the final conversion units.
9. The apparatus according to claim 8, characterized in that, The smooth point interpolation simulation module includes: A smooth domain construction unit is used to obtain the centroid of the tetrahedron containing each triangular face in the target tetrahedral mesh model, and connect the centroid to the three vertices of the corresponding triangular face in sequence to construct the smooth domain of each triangular face; The smooth point interpolation simulation unit is used to continue performing the smooth point interpolation simulation on the region to be converted corresponding to the unit to be converted, based on the smooth domain of each triangular face.
10. The apparatus according to claim 6, characterized in that, The coupling module includes: A coupling interface identification unit is used to identify the coupling interface between the finite element simulation region and the smooth point interpolation simulation region in the target study area. A coupling information transmission unit is used to realize the coupling information transmission between the finite element simulation region and the smooth point interpolation simulation region; The dynamic simulation unit is used to dynamically simulate the deformation of the rock strata in the target study area based on the coupling interface and the coupling information transmission.
11. An electronic device, characterized in that, include: The memory, the processor, and the computer program stored in the memory and capable of running on the processor, the processor executing the program to implement the finite element-smooth point interpolation coupled simulation method for deformation and movement of coal mining strata as described in any one of claims 1-5.
12. A computer-readable storage medium having a computer program stored thereon, characterized in that, The program is executed by the processor to implement the finite element-smooth point interpolation coupled simulation method for deformation and movement of coal mining strata as described in any one of claims 1-5.
13. A computer program product, comprising a computer program, characterized in that, The computer program is executed by a processor to implement the finite element-smooth point interpolation coupled simulation method for deformation and movement of coal mining strata as described in any one of claims 1-5.