CN-121997646-A - Finite element-smooth point interpolation coupling simulation method suitable for deformation movement of coal mining rock stratum

CN121997646ACN 121997646 ACN121997646 ACN 121997646ACN-121997646-A

Abstract

The application relates to the technical field of data simulation processing, in particular to a finite element-smooth point interpolation coupling simulation method suitable for deformation movement of a coal mining rock stratum, which comprises the steps of inputting rock stratum parameter information of a research area into a three-dimensional calculation model, firstly simulating continuous deformation of the rock stratum by using a finite element method, marking a unit with deformation parameters greater than or equal to a deformation threshold value in the rock stratum as a unit to be converted, and then continuing smooth point interpolation simulation on a region to be converted corresponding to the unit to be converted by adopting a smooth point interpolation method; and carrying out finite element simulation on the area corresponding to the unit with the deformation parameter smaller than the deformation threshold by using a finite element method, and coupling a smooth point interpolation simulation process and a finite element simulation process to output a dynamic simulation result of the rock stratum deformation of the research area. The method and the device can analyze the deformation movement rule of the overburden under the influence of mining based on the finite element-smooth point interpolation coupling simulation mode, so that the deformation movement process of the overburden can be accurately simulated.

Inventors

Qin jiayu
XU NENGXIONG
QIN YAN
YANG TIANXIAO
LIU LINAN
MIAO CHENGYU

Assignees

中国地质大学（北京）

Dates

Publication Date: 20260508
Application Date: 20260108

Claims (13)

1. The finite element-smooth point interpolation coupling simulation method suitable for deformation movement of the coal mining rock stratum is characterized by comprising the following steps of: Inputting rock stratum parameter information of a target research area into a pre-constructed three-dimensional calculation model of the target research area, simulating continuous deformation of the rock stratum by using a finite element method, and determining deformation parameters corresponding to target units in the rock stratum; Marking a unit with the deformation parameter larger than or equal to a preset deformation threshold value in the target unit as a unit to be converted, and continuing to perform smooth point interpolation simulation on a region to be converted corresponding to the unit to be converted by adopting a smooth point interpolation method; And carrying out finite element simulation on the region corresponding to the unit with the deformation parameter smaller than the preset deformation threshold by using the finite element method, coupling the smooth point interpolation simulation process and the finite element simulation process, and dynamically simulating the deformation of the rock stratum of the target research area in a mode of finite element-smooth point interpolation coupling simulation so as to output a dynamic simulation result of the rock stratum deformation of the target research area.
2. The method of claim 1, further comprising, prior to simulating the continuous deformation of the formation of the target zone of interest using a finite element method: based on engineering geological data, establishing a three-dimensional geological generalization model of the target research area, and splitting the three-dimensional geological generalization model to obtain an initial three-dimensional calculation model; And applying displacement boundary conditions to the initial three-dimensional calculation model, and solving an initial stress field and a displacement field of the initial three-dimensional calculation model to construct the three-dimensional calculation model meeting preset conditions.
3. The method according to claim 1, further comprising, before continuing the smooth point interpolation simulation for the to-be-converted area corresponding to the to-be-converted unit by using a smooth point interpolation method: acquiring mass center coordinates of all units to be converted, and establishing a minimum bounding box containing the mass center coordinates of all units to be converted; determining a plurality of tetrahedrons corresponding to the centroid coordinates of all the units to be converted in the minimum bounding box; and identifying the tetrahedrons as final units to be converted, and determining the areas to be converted by utilizing the final units to be converted.
4. The method according to claim 3, wherein the performing the smooth point interpolation simulation on the to-be-converted area corresponding to the to-be-converted unit by using a smooth point interpolation method includes: Acquiring the mass center of a tetrahedron where each triangular surface in the target tetrahedron grid model is located, and sequentially connecting the mass center with three vertexes of the corresponding triangular surface to construct a smooth domain of each triangular surface; and continuing to perform the smooth point interpolation simulation on the to-be-converted area corresponding to the to-be-converted unit based on the smooth area of each triangle surface.
5. The method of claim 1, wherein dynamically simulating deformation of the formation of the target zone of interest based on the finite element-smooth point interpolation coupled simulation comprises: Identifying a coupling interface between a finite element simulation region and a smooth point interpolation simulation region in the target research region; The coupling information transfer between the finite element simulation area and the smooth point interpolation simulation area is realized; And dynamically simulating the deformation of the rock stratum of the target research area based on the coupling interface and the coupling information transmission.
6. A finite element-smooth point interpolation coupling simulation device suitable for deformation movement of a coal mining rock stratum, which is characterized by comprising: The finite element simulation module is used for inputting the rock stratum parameter information of the target research area into a pre-constructed three-dimensional calculation model of the target research area, simulating continuous deformation of the rock stratum by using a finite element method, and determining deformation parameters corresponding to target units in the rock stratum; The smooth point interpolation simulation module is used for marking a unit, in the target unit, with the deformation parameter being greater than or equal to a preset deformation threshold value as a unit to be converted, and continuing to perform smooth point interpolation simulation on a region to be converted corresponding to the unit to be converted by adopting a smooth point interpolation method; And the coupling module is used for carrying out finite element simulation on the area corresponding to the unit with the deformation parameter smaller than the preset deformation threshold by utilizing the finite element method, coupling the smooth point interpolation simulation process and the finite element simulation process, and carrying out dynamic simulation on the deformation of the rock stratum of the target research area in a mode of finite element-smooth point interpolation coupling simulation so as to output a dynamic simulation result of the rock stratum deformation of the target research area.
7. The apparatus as recited in claim 6, further comprising: The model construction module is used for establishing a three-dimensional geological generalization model of the target research area based on engineering geological data before simulating the continuous deformation of the rock stratum of the target research area by utilizing a finite element method, and splitting the three-dimensional geological generalization model to obtain an initial three-dimensional calculation model; And the calculation initialization module is used for applying displacement boundary conditions to the initial three-dimensional calculation model before simulating the continuous deformation of the rock stratum of the target research area by using a finite element method, and solving an initial stress field and a displacement field of the initial three-dimensional calculation model so as to construct the three-dimensional calculation model meeting preset conditions.
8. The apparatus as recited in claim 6, further comprising: The searching bounding box module is used for acquiring the barycenter coordinates of all the units to be converted before the smooth point interpolation simulation is continuously carried out on the areas to be converted corresponding to the units to be converted by adopting a smooth point interpolation method, and establishing a minimum bounding box containing the barycenter coordinates of all the units to be converted; The tetrahedron searching module is used for determining a plurality of tetrahedrons corresponding to the centroid coordinates of all the units to be converted in the minimum bounding box before the smooth point interpolation simulation is carried out on the areas to be converted corresponding to the units to be converted by adopting a smooth point interpolation method; and the conversion module is used for identifying the tetrahedrons as final units to be converted before the smooth point interpolation simulation is continuously carried out on the areas to be converted corresponding to the units to be converted by adopting a smooth point interpolation method, and determining the areas to be converted by utilizing the final units to be converted.
9. The apparatus of claim 8, wherein the smooth point interpolation modeling module comprises: The smooth domain construction unit is used for obtaining the mass center of the tetrahedron where each triangular surface in the target tetrahedron grid model is located, and connecting the mass center with three vertexes of the corresponding triangular surface in sequence to construct the smooth domain of each triangular surface; and the smooth point interpolation simulation unit is used for continuing to perform the smooth point interpolation simulation on the to-be-converted area corresponding to the to-be-converted unit based on the smooth domain of each triangle surface.
10. The apparatus of claim 6, wherein the coupling module comprises: the coupling interface identification unit is used for identifying a coupling interface between the finite element simulation area and the smooth point interpolation simulation area in the target research area; The coupling information transfer unit is used for realizing coupling information transfer between the finite element simulation area and the smooth point interpolation simulation area; and the dynamic simulation unit is used for dynamically simulating the deformation of the rock stratum of the target research area based on the coupling interface and the coupling information transmission.
11. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor executing the program to implement the finite element-smooth point interpolation coupled simulation method adapted for deformation movement of a coal formation as claimed in any one of claims 1 to 5.
12. A computer readable storage medium having stored thereon a computer program, the program being executable by a processor for implementing a finite element-smooth point interpolation coupling simulation method adapted to deformation movements of a coal formation according to 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 for implementing a finite element-smooth point interpolation coupled simulation method adapted for deformation movement of a coal-bearing formation according to any one of claims 1-5.

Description

Finite element-smooth point interpolation coupling simulation method suitable for deformation movement of coal mining rock stratum Technical Field The application relates to the technical field of data simulation processing, in particular to a finite element-smooth point interpolation coupling simulation method suitable for deformation movement of a coal mining rock stratum. Background China is the largest coal production and consumption world, and coal mining plays a vital role in national energy safety, economic development and industrial structures. However, coal mining can cause severe disturbance to rock stratum, induce rock stratum movement and earth surface subsidence, and can cause great threat to building construction, traffic engineering, hydraulic engineering and the like in mining areas, and bring outstanding social, economic and environmental problems. Therefore, the accurate analysis of rock deformation and surface subsidence caused by mining is an important precondition for mining planning and design, and has very important significance for safe production and construction of mines. For a long time, students at home and abroad pay attention to the problem of rock stratum movement deformation induced by underground mining, and a theoretical analysis method, a physical model test method and a numerical simulation method are mainly formed. Because of strong adaptability, large calculation scale and high calculation efficiency, the numerical simulation methods such as finite element and finite difference become an important means for analyzing the mechanical behaviors of the induced rock stratum in coal mining. Meanwhile, the finite element and finite difference methods have grid dependence, and grid distortion can occur when the problems are analyzed, so that a simulation result is far from reality and even cannot be calculated. Therefore, in the related art, the exploitation is simulated by adopting an equivalent exploitation mode, taking the simulation of the caving process as an example, a common processing method is to determine the range of the caving zone in advance according to an empirical formula, then assume that the rock mass of the caving zone is an elastoplastic material or a discrete material, and determine the mechanical parameters of the rock mass according to experimental results or the empirical formula. By the method for adjusting the rock mechanical parameters of the damaged area, equivalent simulation of discontinuous behavior can be realized. However, in the related art, the existing numerical simulation methods such as finite element and finite difference have a grid dependence problem when the complex deformation problem of the overburden is handled. Due to grid distortion, continuous deformation and discontinuous deformation are difficult to consider, deviation exists between a simulation result and an actual situation, deformation movement rules of overburden rock under the influence of mining are difficult to analyze, and stability and accuracy of rock stratum movement deformation simulation are reduced. Disclosure of Invention The present application is based on the inventors' knowledge and knowledge of the following problems: For a long time, students at home and abroad pay attention to the problem of rock stratum movement deformation induced by underground mining, and a theoretical analysis method, a physical model test method and a numerical simulation method are mainly formed. Because of strong adaptability, large calculation scale and high calculation efficiency, the numerical simulation methods such as finite element and finite difference become an important means for analyzing the mechanical behaviors of the induced rock stratum in coal mining. However, the methods such as finite element and finite difference still have defects when simulating the problem of rock stratum movement deformation induced by coal mining, and mainly include deformation, cracking, even caving and compaction of the rock stratum due to severe disturbance in the coal mining process. The above process comprises continuous-discontinuous process, and the methods of finite element, finite difference and the like still have the following problems when simulating the above process: 1) The finite element method, the finite difference method and the like are suitable for analyzing continuous problems, and discontinuous behaviors such as cracking, caving and the like cannot be simulated. 2) The finite element and finite difference method has grid dependence, as shown in fig. 1, grid distortion occurs when the problems are analyzed, so that a simulation result is far from reality and cannot be calculated. To solve the above problems, researchers typically simulate mining in an equivalent mining manner. Taking the simulation of the caving process as an example, a common processing method is to determine the range of the caving zone in advance according to an empirical formula, then assume the rock m