CN-122021052-A - Calculation method and system for tunnel roof anchoring supporting impact resistance

CN122021052ACN 122021052 ACN122021052 ACN 122021052ACN-122021052-A

Abstract

The invention provides a calculation method and a calculation system for the impact resistance of an anchoring support of a roadway roof, relates to the technical field of calculation for the impact resistance of the anchoring support, acquires geological and support parameters, combines a three-dimensional model with construction of a stress coupling matrix, realizes high-precision simulation of dynamic response of each anchor rod of the roadway roof under various vibration waves, further obtains more fitting actual stress distribution, can better simulate different situations possibly occurring in a region where the roadway is located by adopting simulation setting of random vibration sources and multiple types of vibration waves, combines probability distribution fitting, can scientifically evaluate failure risks of each anchor rod, and performs impact grade division by counting proportion of risk grids, thereby providing quantized risk judgment basis for overall anchoring support safety. The scheme breaks through the limitation of traditional static analysis, enhances the evaluation capability of the impact resistance of the roadway support and improves the reliability of the evaluation result.

Inventors

Nie jigang
WANG CHAO
ZHAO YONGZHI
GONG ZHILI
CAO XU
DUAN WENPENG
DING KEKE
ZHANG DAOWEN
HOU HUAYING
LIU ZHENGFANG
NIE SHILONG
ZHANG XUEJING
WANG NANNAN
HUANG SHUAI
WANG XIAO
Sha Nanxing
LIU HAIRUI
SONG QINGDE
LIU ERCENG
YAO XIN
Zheng Yushuai
TIAN YAKUI
LU LIN

Assignees

河南龙宇能源股份有限公司

Dates

Publication Date: 20260512
Application Date: 20260226

Claims (9)

1. The method for calculating the impact resistance of the tunnel roof anchoring support is characterized by comprising the following specific steps of: S1, collecting geological parameters of a roadway and supporting parameters of an anchoring supporting structure, building a three-dimensional model according to the geological parameters and the supporting parameters, equally dividing a roadway top plate into a plurality of groups of grids by taking each anchor rod in the anchoring supporting structure as a center, and sequentially applying impact loads to different grids to carry out simulation analysis; S2, sequentially collecting the maximum stress suffered by each grid in simulation analysis, calculating stress coupling coefficients among different grids according to the maximum stress of each grid so as to quantify the mechanical influence among different grids, and generating a stress coupling matrix according to the stress coupling coefficients so as to reflect the stress distribution at the top plate of the roadway; S3, generating a geological active region in a three-dimensional model based on a historical seismic record and a geological theory of the roadway, randomly setting a plurality of groups of vibration sources in the geological active region to simulate the impact process of vibration waves on a roadway top plate, sequentially collecting time sequence parameters of each group of vibration waves transmitted to the roadway top plate, wherein each group of vibration sources corresponds to a random vibration wave type; and S4, based on the time sequence parameters and the stress coupling matrix, calculating the maximum stress born by different grids in each impact, fitting probability distribution of the maximum stress of each grid, generating failure probability of each grid corresponding to the anchor rod according to the probability distribution, and finally calculating the impact resistance grade of the whole roadway roof anchoring support to evaluate the impact resistance capability of the roadway roof anchoring support.
2. The method for calculating the impact resistance of the tunnel roof anchoring support according to claim 1, wherein the geological parameters comprise rock mass type, rock mass strength and rock stratum thickness of the area where the tunnel is located, and the support parameters comprise the number of anchor rods, the support positions, the anchor rod size and the anchor rod material.
3. The method for calculating the impact resistance of the tunnel roof anchoring support according to claim 2, wherein when impact loads are applied to different grids for simulation analysis, the impact loads are applied to the whole grid in a surface load mode, and the maximum stress at the joint of the anchor rod and the tunnel roof is taken as the maximum stress of the corresponding grid.
4. The method for calculating the impact resistance of the tunnel roof anchoring support according to claim 3, wherein the specific flow of the step S2 is as follows: s201, sequentially collecting the maximum stress of each grid, and comparing the maximum stress of each grid with the reference stress before impact load is applied to obtain the stress increment of each grid; s202, for each grid applying impact load, calculating the ratio between the stress increment of other grids and the stress increment of the grid, and taking the ratio as the stress coupling coefficient of the grid to the other grids respectively; S203, constructing a stress increment matrix based on the stress coupling coefficient, wherein the stress increment matrix is used for reflecting the stress increment of the roadway top plate when a single grid is subjected to impact load, the elements of the stress increment matrix correspond to the grids one by one, and the number of the elements is the same as that of the grids; S204, superposing a plurality of groups of stress increment matrixes to obtain a composite increment matrix so as to reflect the stress increment of the roadway roof when a plurality of grids are subjected to impact load, constructing a reference matrix based on the reference stress of each grid, and superposing the composite increment matrix and the reference matrix to obtain a stress coupling matrix.
5. The method of claim 4, wherein for any element in the stress delta matrix, the value of the element is equal to the stress delta of the grid multiplied by the stress coupling coefficient of the grid relative to other grids.
6. The method for calculating the impact resistance of the tunnel roof anchor support according to claim 1, wherein the historical seismic records comprise historical vibration source positions, historical vibration wave types and historical vibration source energy, and the vibration wave types comprise P waves, S waves and surface waves; when the impact process of the vibration wave to the roadway roof is simulated, the vibration source energy is set randomly, and the random range is between the maximum and minimum values of the historical vibration source energy.
7. The method for calculating the impact resistance of the tunnel roof anchoring support according to claim 4, wherein the time interval of the time sequence parameter is from the time when the vibration wave just propagates to the tunnel roof to the time when the vibration wave completely leaves the tunnel roof, and the time sequence parameter is the grid number of the grid impacted by the vibration wave at different time points and the real-time stress of the grid impacted by the vibration wave.
8. The method for calculating the impact resistance of the tunnel roof anchoring support according to claim 7, wherein the specific flow of the step S4 is as follows: s401, for single impact, calculating direct stress increment of different grids according to corresponding time sequence parameters by combining vibration source positions and vibration source energy, and calculating real-time stress received by each grid by combining the direct stress increment and a stress coupling matrix, and obtaining the maximum stress received by each grid in the impact according to the real-time stress increment; S402, carrying out probability distribution estimation on the maximum stress of each grid in multiple impacts, fitting a cumulative distribution function of the maximum stress of each grid, and calculating the failure probability of the corresponding anchor rod of the grid based on the cumulative distribution function, wherein the failure probability and the maximum stress are positively correlated; And S403, calibrating grids with failure probability larger than a preset probability threshold as risk grids, counting the proportion of the risk grids in all grids, comparing the proportion with a preset proportion threshold, and outputting the impact resistance level of the whole tunnel roof anchoring support according to the comparison result, wherein the lower the proportion of the risk grids is, the higher the impact resistance level is, and the stronger the impact resistance is.
9. The system for calculating the impact resistance of the tunnel roof anchoring support is characterized by being suitable for the method for calculating the impact resistance of the tunnel roof anchoring support according to any one of claims 1-8, and specifically comprising the following steps: the three-dimensional simulation module is used for collecting geological parameters of a roadway and supporting parameters of an anchoring support, building a three-dimensional model according to the geological parameters and the supporting parameters, equally dividing the roadway top plate into a plurality of groups of grids by taking each anchor rod in the anchoring support as a center, and sequentially applying impact loads to different grids to carry out simulation analysis; the data analysis module is used for sequentially collecting the maximum stress suffered by each grid in simulation analysis, calculating stress coupling coefficients among different grids according to the maximum stress of each grid so as to quantify the mechanical influence among different grids, and generating a stress coupling matrix according to the stress coupling coefficients so as to reflect the stress distribution at the roadway roof; The system comprises a shock simulation module, a vibration source and a vibration source, wherein the shock simulation module is used for generating a geological active area in a three-dimensional model based on historical seismic records and geological theory of a roadway, randomly setting a plurality of groups of vibration sources in the geological active area to simulate the shock process of vibration waves to a roadway top plate, sequentially collecting time sequence parameters of each group of vibration waves transmitted to the roadway top plate, and each group of vibration sources corresponds to a random vibration wave type; the capacity evaluation module is used for calculating the maximum stress received by different grids in each impact based on the time sequence parameters and the stress coupling matrix, fitting probability distribution of the maximum stress of each grid, generating failure probability of the corresponding anchor rod of each grid according to the probability distribution, and finally calculating the impact resistance level of the whole roadway roof anchoring support to evaluate the impact resistance capacity of the roadway roof anchoring support.

Description

Calculation method and system for tunnel roof anchoring supporting impact resistance Technical Field The invention relates to the technical field of calculation of impact resistance of anchoring supports, in particular to a calculation method and a system of impact resistance of a tunnel roof anchoring support. Background Under complex geological conditions, the factors such as heterogeneity, crack distribution and faults of the rock mass cause significant differences in vibration wave propagation paths and intensities, the traditional design and analysis methods depend on static loads or simplified dynamic models, and the mutual influence among different areas of a roadway roof in the vibration wave propagation process is difficult to accurately reflect, so that the stress distribution of the whole model is inaccurate. In addition, in the prior art, the statistical analysis of the stress extremum of the anchor rod under the impact of multiple vibration waves is insufficient, a reliable failure probability assessment and risk classification system is lacked, and scientific basis is difficult to provide for engineering decision. In the prior art, publication number CN111259542A discloses a calculation method of the impact resistance of a tunnel roof anchoring support, which comprises the steps of firstly calculating the energy before the tunnel roof anchor rod and the anchor cable are subjected to limit deformation and damage respectively, then calculating the energy limit value absorbed by a roof support system before the damage according to the interval distance, and judging the impact resistance of the current roof anchoring net cable support according to the comparison result, wherein the energy applied by the roof to the support system is the sum of kinetic energy and potential energy when rock burst occurs, and in a critical state, the energy is equal to the energy limit value of the support system, so that the minimum speed for causing the failure of the roof support system is obtained, and the rock mass maximum vibration speed corresponding to the maximum mine vibration energy level obtained by monitoring is obtained. According to the scheme, although the judgment of the impact resistance of the support system can be realized, the roadway top plate is regarded as a whole, calculation is carried out based on the static energy limit value, and when the impact occurs in a part of the roadway top plate, the mutual influence among all areas is ignored, so that the local response of the anchor rod is difficult to be carefully reflected, the whole stress distribution of the model is inaccurate, and the final evaluation result accuracy is low. The above information disclosed in the background section is only for enhancement of understanding of the background of the disclosure and therefore it may include information that does not form the prior art that is already known to a person of ordinary skill in the art. Disclosure of Invention The invention aims to provide a method and a system for calculating the impact resistance of a tunnel roof anchoring support, so as to solve the problems in the background technology. In order to achieve the above purpose, the present invention provides the following technical solutions: The method for calculating the impact resistance of the tunnel roof anchoring support comprises the following specific steps: S1, collecting geological parameters of a roadway and supporting parameters of an anchoring supporting structure, building a three-dimensional model according to the geological parameters and the supporting parameters, equally dividing a roadway top plate into a plurality of groups of grids by taking each anchor rod in the anchoring supporting structure as a center, and sequentially applying impact loads to different grids to carry out simulation analysis; S2, sequentially collecting the maximum stress suffered by each grid in simulation analysis, calculating stress coupling coefficients among different grids according to the maximum stress of each grid so as to quantify the mechanical influence among different grids, and generating a stress coupling matrix according to the stress coupling coefficients so as to reflect the stress distribution at the top plate of the roadway; S3, generating a geological active region in a three-dimensional model based on a historical seismic record and a geological theory of the roadway, randomly setting a plurality of groups of vibration sources in the geological active region to simulate the impact process of vibration waves on a roadway top plate, sequentially collecting time sequence parameters of each group of vibration waves transmitted to the roadway top plate, wherein each group of vibration sources corresponds to a random vibration wave type; and S4, based on the time sequence parameters and the stress coupling matrix, calculating the maximum stress born by different grids in each impact, fitting probability distribution of the m