CN-121302955-B - Method and device for simulating flow characteristics of grouting material of coal mine based on LBM and DEM and computer equipment

CN121302955BCN 121302955 BCN121302955 BCN 121302955BCN-121302955-B

Abstract

A coal mine grouting material flow characteristic simulation method based on LBM and DEM comprises the steps of initializing a simulation space, determining a plurality of grid points in the simulation space, setting a plurality of simulation barriers in the simulation space, covering the grid points, traversing all the grid points, determining part of the grid points as solids, generating flow constraint limits based on the grid points determined as solids, initializing physical parameters of the LBM, model parameters of a D3Q19 model, macroscopic amounts of grouting fluid, grouting fluid distribution functions and particle parameters of the DEM, performing coupling iteration on the LBM and the DEM, generating flow dynamic data of grouting particles in grouting fluid in the simulation space, and performing visual display on the flow dynamic data in real time. The invention can realize the simulation test of the flow characteristic of the grouting material based on a computer, does not need to manufacture a complex testing device, can greatly reduce the cost and improves the testing efficiency.

Inventors

ZHOU AITAO
SONG FANGZHOU
WANG KAI
ZHU YUXUAN
Pan Bojie
Hao Yushuang

Assignees

中国矿业大学(北京)

Dates

Publication Date: 20260508
Application Date: 20250919

Claims (8)

1. The coal mine grouting material flow characteristic simulation method based on LBM and DEM is characterized by comprising the following steps: initializing an analog space, and determining a plurality of grid points in the analog space; Providing a plurality of simulation obstacles in the simulation space, wherein the simulation obstacles cover a plurality of grid points; After generating a flow constraint limit based on the grid points determined to be solid, ignoring the grid points located on a spatial boundary of the simulation space and the grid points determined to be solid and marking the remaining grid points as fluid points, thereby creating a plurality of fluid grids based on the fluid points; After the generating of the flow constraint limits based on the grid points determined to be solid, estimating, for each grid point, a possible critical type based on a possible discrete motion direction, the possible critical type being a spatial critical or a constraint critical, the spatial critical characterizing that the grid point is capable of impinging on the spatial boundary of the simulation space based on the possible discrete motion direction, the constraint critical characterizing that the grid point is capable of impinging on the flow constraint limits based on the possible discrete motion direction; Initializing physical parameters of the LBM (Lattice Boltzmann Method ), model parameters of the D3Q19 model, macroscopic quantity of grouting fluid, grouting fluid distribution function, and particle parameters of the DEM (DISCRETEELEMENT METHOD, discrete unit method); Performing coupling iteration on the LBM and the DEM to generate flow dynamic data of grouting particles in the grouting fluid in the simulation space; and carrying out visual display on the flow dynamic data in real time.
2. A method of modeling flow characteristics of a grouting material for a coal mine based on LBM and DEM as claimed in claim 1, wherein said method of initializing the modeling space comprises: Creating a three-dimensional coordinate system xyz, and determining a three-axis coordinate length of the three-dimensional coordinate system xyz; Determining a first point number in an x direction, a second point number in a y direction and a third point number in a z direction in the three-dimensional coordinate system xyz; Determining a number of the grid points based on the first points, the second points, and the third points, and generating the grid points based on the three-axis coordinate lengths; creating a spatial boundary with the grid points at edge locations, the spatial boundary enclosing the simulation space.
3. A method of modeling flow characteristics of a grouting material for a coal mine based on LBM and DEM as claimed in claim 1, wherein said method of providing a plurality of modeling obstacles in said modeling space comprises: Determining a number of the simulated obstacles and determining a geometry and a geometric dimension of the simulated obstacles; Generating a virtual entity based on the geometry and the geometry; Placing the virtual entity into the simulation space and causing the virtual entity to cover a plurality of the grid points; And when the number of the simulated obstacles is greater than or equal to two, adjusting the positions of the virtual entities in the simulation space until the distance between two adjacent virtual entities meets the preset obstacle distribution constraint condition.
4. A method of modeling flow characteristics of a coal mine grouting material based on an LBM and a DEM as claimed in claim 1, wherein the physical parameters of the LBM include a total time step, a first time step, a relaxation time, a relaxation factor, an interaction strength parameter and a potential function scaling factor; The model parameters of the D3Q19 model comprise a discrete speed vector array, a reverse speed index array and a weight array, wherein the reverse speed index array corresponds to elements in the weight array one by one; The grouting fluid macroscopic quantity comprises a fluid density value, a high-low density average value and a boundary adsorption force reference density; the particle parameters of the DEM include a second time step, a kinematic viscosity, a total number of particles, and particle physical characteristic parameters, wherein the particle physical characteristic parameters include a particle physical size and a particle mass.
5. A method of modeling flow characteristics of a grouting material for a coal mine based on LBM and DEM as claimed in claim 1, further comprising the steps of: A linear index is established for each of the grid points, the linear index including a spatial location characterizing the grid point location and a classification result of the grid point.
6. A coal mine grouting simulation system based on LBM and DEM, for implementing a coal mine grouting material flow characteristic simulation method based on LBM and DEM as claimed in any one of claims 1 to 5, the system comprising: a space creation module for initializing a simulation space, determining a plurality of grid points in the simulation space, setting a plurality of simulation obstacles in the simulation space, and generating a flow constraint limit based on the grid points determined to be solid; The coupling iteration module is used for carrying out coupling iteration on the LBM and the DEM to generate flow dynamic data of grouting particles in the grouting fluid in the simulation space; And the visual display module is used for visually displaying the flow dynamic data in real time.
7. A computer device, comprising: A memory for storing a computer program; a processor for reading and executing the computer program to implement a method for simulating flow characteristics of grouting material for coal mines based on LBM and DEM as claimed in any one of claims 1 to 5.
8. A storage medium storing a computer program which when executed implements a method of modeling flow characteristics of a grouting material for a coal mine based on LBM and DEM as claimed in any one of claims 1 to 5.

Description

Method and device for simulating flow characteristics of grouting material of coal mine based on LBM and DEM and computer equipment Technical Field The invention relates to the technical field of coal mine grouting material performance test, in particular to a coal mine grouting material flow characteristic simulation method and device based on LBM and DEM and computer equipment. Background Through researches for many years, coal mines in China form a roadway safe and efficient supporting technology based on geomechanical tests and taking anchoring and grouting as cores. In the roadway with simple shallow geological conditions and relatively complete surrounding rock, the anchor rod and anchor cable supporting mode is widely applied, and the technical and economic benefits are obvious. Along with the continuous increase of the coal mining depth, the deep mine of kilometers in China is more and more, and the control difficulty of surrounding rocks of a roadway is obviously increased. Under deep high stress and strong mining conditions, surrounding rock rheological property is strong, deformation is large, duration is long, deformation and damage of surrounding rock of a roadway are reduced to a certain extent due to high prestress and strong anchor bolt support, but the problems of anchor bolt breakage, support member failure, large deformation of surrounding rock of the roadway and the like still exist, and safe and efficient production of a coal mine is affected. The low-strength attribute characteristic of the soft surrounding rock under the action of high-stress strong mining is an inherent cause of strong deformation and damage. By means of high-pressure splitting grouting, active modification of the soft surrounding rock is achieved, strength characteristics, mechanical properties and bearing capacity of the soft surrounding rock are improved, the problem of low-strength strong deformation mechanical characteristic attribute of the soft surrounding rock is broken through, and therefore the problem of strong deformation and damage of deep high-stress soft surrounding rock is solved. However, the surrounding rock grouting under the deep mining condition of the coal mine needs to overcome high stress, surrounding rock strength, loss of grouting pressure caused by slurry viscosity and the like, so that a grouting system is required to have higher grouting pressure, the purpose of grouting reinforcement of the surrounding rock can be effectively achieved, and the surrounding rock deformation and damage aggravation caused by surrounding rock damage and damage are avoided. In the prior art, performance testing of grouting materials for grouting of coal mines is realized by means of single software simulation or a specially designed testing device, and parameters such as flow performance of the grouting materials are determined by carrying out grouting tests in a small range in a laboratory. The main problems of the method are that the structural design of the device is complex, and in the research and development process of grouting materials, when parameters such as components, proportions and the like are changed, the device needs to be retested, and is too complex. Disclosure of Invention In order to solve the defects in the prior art, the invention provides a coal mine grouting material flow characteristic simulation method, a coal mine grouting material flow characteristic simulation device and computer equipment based on LBM and DEM technologies, which can realize simulation test on grouting material flow characteristics based on a computer without manufacturing a complex test device, can greatly reduce cost, and can simulate test on different grouting materials by adjusting parameters, thereby greatly reducing operation complexity and greatly improving test efficiency. In order to achieve the above purpose, the invention adopts the following specific scheme: a coal mine grouting material flow characteristic simulation method based on LBM and DEM comprises the following steps: initializing an analog space, and determining a plurality of grid points in the analog space; Providing a plurality of simulation obstacles in the simulation space, wherein the simulation obstacles cover a plurality of grid points; traversing all of the grid points, determining some of the grid points therein as solid, and generating flow constraint limits based on the grid points determined as solid; Initializing physical parameters of the LBM (Lattice Boltzmann Method ), model parameters of the D3Q19 model, macroscopic quantity of grouting fluid, grouting fluid distribution function, and particle parameters of the DEM (DISCRETEELEMENT METHOD, discrete unit method); Performing coupling iteration on the LBM and the DEM to generate flow dynamic data of grouting particles in the grouting fluid in the simulation space; and carrying out visual display on the flow dynamic data in real time. Preferably, the method for initializing the s