CN113392561B - Method and device for realizing energy field analysis of surrounding rock around roadway - Google Patents
Method and device for realizing energy field analysis of surrounding rock around roadway Download PDFInfo
- Publication number
- CN113392561B CN113392561B CN202110728033.8A CN202110728033A CN113392561B CN 113392561 B CN113392561 B CN 113392561B CN 202110728033 A CN202110728033 A CN 202110728033A CN 113392561 B CN113392561 B CN 113392561B
- Authority
- CN
- China
- Prior art keywords
- strain
- energy
- roadway
- model
- representing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims abstract description 104
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000004458 analytical method Methods 0.000 title claims abstract description 35
- 238000005381 potential energy Methods 0.000 claims abstract description 69
- 230000010354 integration Effects 0.000 claims abstract description 9
- 230000008859 change Effects 0.000 claims description 31
- 238000009825 accumulation Methods 0.000 claims description 11
- 239000000463 material Substances 0.000 claims description 10
- 230000003014 reinforcing effect Effects 0.000 claims description 9
- 239000003245 coal Substances 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 abstract description 7
- 238000013461 design Methods 0.000 abstract description 7
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 14
- 229910000831 Steel Inorganic materials 0.000 description 5
- 239000010959 steel Substances 0.000 description 5
- 230000002776 aggregation Effects 0.000 description 4
- 238000004220 aggregation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical group C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000009529 body temperature measurement Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000002076 thermal analysis method Methods 0.000 description 1
- 230000005641 tunneling Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Devices Affording Protection Of Roads Or Walls For Sound Insulation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method and a device for realizing energy field analysis of surrounding rocks around a roadway, wherein the method comprises the following steps: FLAC (FLAC) 3D Constructing model units corresponding to surrounding rocks around a roadway on software, performing triple integration on the total strain potential energy density of the model units to obtain strain energy of the system, and assigning the strain energy of the system to the FLAC 3D The internal cohesion display module of the software displays the current energy field to form a variable cloud picture, and the energy field is analyzed according to the variable cloud picture; the invention has the advantages that: the method for analyzing the energy field of the surrounding rock around the roadway is beneficial to analyzing the surrounding rock of the roadway, guiding site construction and optimizing roadway support design.
Description
Technical Field
The invention relates to the field of energy field analysis of surrounding rocks of coal mine roadways, in particular to a method and a device for realizing energy field analysis of surrounding rocks around roadways.
Background
In elastoplastics, when an elastic body is subjected to an external force, deformation is inevitably generated, and meanwhile, the potential energy of the external force is also changed. When an external force is slowly applied to the object (without causing the object to generate acceleration motion), the external force is regarded as static force, the kinetic energy of the system can be ignored, and meanwhile, the consumption of other energy (such as heat energy and the like) is omitted, so that the change of the potential energy of the external force is completely converted into strain energy (potential energy) to be stored in the object.
When the roadway is excavated or the working face is stoped, surrounding rock is in a natural balance state, namely an original rock stress state, and meanwhile, the energy stored in the surrounding rock is in a balance state, namely an original energy field, which is regarded as a static energy field, after the roadway is excavated or the working face is stoped, the original stress balance is destroyed, the internal stress of the surrounding rock is redistributed according to the constitutive relation of rock media to form a new balance state, namely a secondary stress field, the surrounding rock inevitably generates strain in the process of stress rebalancing of the surrounding rock, and the energy field in the surrounding rock is redistributed to form a new energy field. The method realizes visualization of the energy field of surrounding rocks around the coal mine tunnel, analyzes the energy aggregation state around the excavated body, can be used for stability analysis of surrounding rocks of the tunnel, prediction of surrounding rock impact tendency and impact occurrence positions, prediction of damage of anchor cable supporting structures of the tunnel anchor rods and the like, guides site construction, optimizes tunnel supporting design, and pre-judges the possibility of tunnel impact and impact occurrence positions by combining surrounding rock impact tendency index values, so that measures can be taken in time, and occurrence of disasters such as supporting instability, surrounding rock impact damage and the like can be prevented.
FLAC 3D Finite difference program FLAC being two-dimensional 2D Can simulate the stress characteristics of three-dimensional structures of soil, rock and other materials and analyze plastic flow. The actual structure is fitted by adjusting the polyhedral cells in the three-dimensional grid. The unit material can adopt linear or nonlinear constitutive model, under the action of external force, when the material is subjected to yielding flow, the grid can be correspondingly deformed and moved (largeDeformation mode). FLAC (FLAC) 3D By adopting the explicit Lagrangian algorithm and the mixed-discrete partition technology, the plastic damage and flow of the material can be simulated very accurately. Because the stiffness matrix is not required to be formed, a large-scale three-dimensional problem can be solved based on a smaller memory space.
Traditional FLAC 3D Numerical simulation analysis mainly analyzes stress cloud patterns, surface displacement, plastic region distribution and the like, for example, chinese patent application No. 202010229878.8 discloses a method based on FLAC 3D The roadway surrounding rock thermophysical parameter inversion method of the numerical software is characterized by comprising the following steps of: 1. underground on-site temperature measurement provides a foundation for parameter inversion; 2. giving out a simulation working condition by using an orthogonal test method; 3. establishing and using FLAC 3D The thermal analysis module of the software establishes a thermal conduction numerical calculation model; 4. obtaining a calculation result of each working condition; 5. giving out an optimal value of the thermal physical parameters of the surrounding rock by using a neural network prediction model; 6. verifying the calculation result of the inverse analysis, substituting the parameters into FLAC 3D The model is calculated and compared with the actual measurement results. According to the method, the consumption of manpower and material resources is reduced, the measured temperature value of the roadway is provided, the surrounding rock thermal parameters are reversely calculated according to the measured value of a laboratory, errors in actual operation can be effectively reduced, possibility is provided for accurately calculating the surrounding rock temperature field of the roadway, but the energy field analysis of surrounding rocks around the roadway is not carried out, so that the stability analysis of the surrounding rocks of the roadway, the prediction of the impact tendency and the impact occurrence part of the surrounding rocks and the prediction of the damage of the anchor cable supporting structure of the roadway anchor rod cannot be carried out, the on-site construction is difficult to guide, and the roadway supporting design is difficult to optimize.
Disclosure of Invention
The technical problem to be solved by the invention is that the prior art lacks a method for analyzing the energy field of surrounding rocks around a roadway, so that the stability analysis of the surrounding rocks of the roadway, the prediction of the impact tendency and the impact occurrence position of the surrounding rocks and the prediction of the damage of the anchor cable supporting structure of the anchor rod of the roadway cannot be carried out, the on-site construction is difficult to guide, and the optimization of the supporting design of the roadway is difficult.
The invention solves the above problems by the following technical meansThe technical problems are as follows: a method of enabling energy field analysis of surrounding rock around a roadway, the method comprising: FLAC (FLAC) 3D Constructing model units corresponding to surrounding rocks around a roadway on software, performing triple integration on the total strain potential energy density of the model units to obtain strain energy of the system, and assigning the strain energy of the system to the FLAC 3D And a cohesion display module built in the software displays the current energy field to form a variable cloud picture, and the energy field is analyzed according to the variable cloud picture.
The invention utilizes FLAC 3D The software performs triple integration on the total strain potential energy density of the model unit to obtain the strain energy of the system, so as to obtain the current energy field, and the current energy field passes through FLAC 3D The built-in graphic display module of the software displays the current energy field to form a variable cloud picture, the energy field is analyzed according to the variable cloud picture, the method for analyzing the energy field of surrounding rocks around the roadway is provided, the analysis of the surrounding rocks of the roadway is facilitated, the site construction is guided, and the roadway support design is optimized.
Further, the strain energy of the system is obtained by the following steps:
by the formulaObtaining work done by external force to deform the model unit in the x direction, wherein epsilon x Representing the strain of the model element in the x-direction, sigma x Representing the positive stress of the model element in the x-direction;
by the formulaObtaining work done by external force to deform the model unit in y direction, wherein epsilon y Representing the strain of the model element in the y-direction, σ y Representing the positive stress of the model element in the y-direction;
by the formulaObtaining work done by external force to deform the model unit in the z direction, wherein epsilon z Representing the strain of the model element in the z-direction, sigma z Representing the positive stress of the model element in the z direction;
by the formula w=w x +W y +W z Obtaining the total work of the external force on the model unit in the three-dimensional direction;
the total work of the model unit in three-dimensional deformation is converted into the strain energy of the system, so that the strain energy formula of the system is constructed
Wherein U is 0x Represents strain energy per unit volume formed by deformation of the model unit in the x direction, U 0y Represents strain energy per unit volume formed by deformation of the model element in the y direction, U 0z Represents strain energy per unit volume formed by deformation of the model unit in the z direction, U 0 Representing the total strain potential energy density of the model element.
Further, the deformation of the object obtained by the plastic theory can be decomposed into two parts, one part is the change of volume and the other part is the change of shape, so that the total strain potential energy density of the model unit is the sum of the strain energy stored in the unit volume due to the change of volume and the strain energy stored in the unit volume due to the change of shape.
Further, by the formulaAcquiring strain energy stored in a unit volume due to a volume change, wherein σ 0 Representing the stress before deformation of the model element, i.e. the original stress, K is the bulk modulus and K is a constant.
Further, by the formulaAcquiring strain energy stored in a unit volume due to a shape change, wherein σ 1 Representing the preset strain force, sigma, of the model unit in the x direction 2 Representing the preset strain force, sigma, of the model unit in the y direction 3 Representing model elements in the z-directionUpward preset strain force, G represents shear modulus and G takes a constant.
Further, by the formula
And obtaining the total strain potential energy density of the model unit.
Further, the analyzing the energy field according to the variable cloud image includes: and (5) finding a potential energy gathering area according to the variable cloud chart, and reinforcing and supporting surrounding rocks around the roadway where the potential energy gathering area is located.
Further, the larger the energy concentration of the potential energy accumulation area is, the higher the stress level of the corresponding position is, the larger the displacement of the surrounding rock is, and the larger the damage to the surrounding rock and the supporting structure of the surrounding rock is.
Further, the method for reinforcing and supporting surrounding rocks around the roadway where the potential energy gathering area is located comprises the following steps: the material strength of the supporting structure of surrounding rocks around the roadway is improved, the surrounding rocks around the roadway are injected along the empty coal pillar, or the supporting roof plate of the surrounding rocks around the roadway is designed into an arch structure.
The invention also provides a device for realizing the energy field analysis of surrounding rocks around a roadway, which comprises: strain energy acquisition module of system for FLAC 3D Constructing a model unit corresponding to surrounding rock around a roadway on software, and carrying out triple integration on the total strain potential energy density of the model unit to obtain the strain energy of the system;
a variable cloud image acquisition module for assigning strain energy of the system to FLAC 3D A cohesion display module built in the software displays a variable cloud picture formed by the current energy field;
and the analysis module is used for analyzing the energy field according to the variable cloud image.
Further, the strain energy of the system is obtained by the following steps:
by the formulaThe work done by the external force to deform the model unit in the x direction is obtained,wherein ε x Representing the strain of the model element in the x-direction, sigma x Representing the positive stress of the model element in the x-direction;
by the formulaObtaining work done by external force to deform the model unit in y direction, wherein epsilon y Representing the strain of the model element in the y-direction, σ y Representing the positive stress of the model element in the y-direction;
by the formulaObtaining work done by external force to deform the model unit in the z direction, wherein epsilon z Representing the strain of the model element in the z-direction, sigma z Representing the positive stress of the model element in the z direction;
by the formula w=w x +W y +W z Obtaining the total work of the external force on the model unit in the three-dimensional direction;
the total work of the model unit in three-dimensional deformation is converted into the strain energy of the system, so that the strain energy formula of the system is constructed
Wherein U is 0x Represents strain energy per unit volume formed by deformation of the model unit in the x direction, U 0y Represents strain energy per unit volume formed by deformation of the model element in the y direction, U 0z Represents strain energy per unit volume formed by deformation of the model unit in the z direction, U 0 Representing the total strain potential energy density of the model element.
Further, the deformation of the object obtained by the plastic theory can be decomposed into two parts, one part is the change of volume and the other part is the change of shape, so that the total strain potential energy density of the model unit is the sum of the strain energy stored in the unit volume due to the change of volume and the strain energy stored in the unit volume due to the change of shape.
Further, by the formulaAcquiring strain energy stored in a unit volume due to a volume change, wherein σ 0 Representing the stress before deformation of the model element, i.e. the original stress, K is the bulk modulus and K is a constant.
Further, by the formulaAcquiring strain energy stored in a unit volume due to a shape change, wherein σ 1 Representing the preset strain force, sigma, of the model unit in the x direction 2 Representing the preset strain force, sigma, of the model unit in the y direction 3 Representing the preset strain force of the model element in the z direction, G represents the shear modulus and G takes a constant.
Further, by the formula
And obtaining the total strain potential energy density of the model unit.
Further, the analysis module is further configured to: and (5) finding a potential energy gathering area according to the variable cloud chart, and reinforcing and supporting surrounding rocks around the roadway where the potential energy gathering area is located.
Further, the larger the energy concentration of the potential energy accumulation area is, the higher the stress level of the corresponding position is, the larger the displacement of the surrounding rock is, and the larger the damage to the surrounding rock and the supporting structure of the surrounding rock is.
Further, the method for reinforcing and supporting surrounding rocks around the roadway where the potential energy gathering area is located comprises the following steps: the material strength of the supporting structure of surrounding rocks around the roadway is improved, the surrounding rocks around the roadway are injected along the empty coal pillar, or the supporting roof plate of the surrounding rocks around the roadway is designed into an arch structure.
The invention has the advantages that: the invention utilizes FLAC 3D Software for triple integral acquisition system of total strain potential energy density of model unitEnergy conversion to obtain the current energy field through FLAC 3D The built-in graphic display module of the software displays the current energy field to form a variable cloud picture, the energy field is analyzed according to the variable cloud picture, the method for analyzing the energy field of surrounding rocks around the roadway is provided, the analysis of the surrounding rocks of the roadway is facilitated, the site construction is guided, and the roadway support design is optimized.
Drawings
FIG. 1 is a flow chart of a method for implementing energy field analysis of surrounding rock around a roadway according to an embodiment of the present invention;
FIG. 2 is a schematic diagram showing deformation of an ABCD unit in an x direction in a method for realizing energy field analysis of surrounding rocks around a roadway according to an embodiment of the present invention;
FIG. 3 is a graph showing a potential energy aggregation distribution diagram of a method for implementing energy field analysis of surrounding rock around a roadway in a specific application scenario according to an embodiment of the present invention;
FIG. 4 is an energy distribution diagram after arch roadway excavation in a method for realizing energy field analysis of surrounding rock around a roadway according to an embodiment of the present invention;
FIG. 5 is an energy distribution diagram after rectangular roadway excavation in a method for realizing energy field analysis of surrounding rock around a roadway according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a situation of concentrated energy when there is no support in a roadway in a method for analyzing an energy field of surrounding rock around the roadway according to an embodiment of the present invention;
fig. 7 is a schematic diagram of a situation of concentrated potential energy when 2 anchor bolts are used for supporting a roadway in a method for analyzing an energy field of surrounding rock around the roadway according to an embodiment of the present invention;
fig. 8 is a schematic diagram of a potential energy concentration situation when 4 anchor bolts are used for supporting a roadway in the method for analyzing the energy field of surrounding rock around the roadway according to the embodiment of the invention;
fig. 9 is a schematic diagram of a potential energy concentration situation when 8 anchor bolts are used for supporting a roadway in the method for analyzing the energy field of surrounding rock around the roadway according to the embodiment of the invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions in the embodiments of the present invention will be clearly and completely described in the following in conjunction with the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Example 1
As shown in fig. 1, a method for implementing energy field analysis of surrounding rock around a roadway includes: FLAC (FLAC) 3D Constructing model units corresponding to surrounding rocks around a roadway on software, performing triple integration on the total strain potential energy density of the model units to obtain strain energy of the system, and assigning the strain energy of the system to the FLAC 3D And a cohesion display module built in the software displays the current energy field to form a variable cloud picture, and the energy field is analyzed according to the variable cloud picture. The following describes in detail the implementation of the invention:
the strain energy of the system is obtained by the following steps: as shown in FIG. 2, the external force acting on the ABCD units is AD and sigma on the CB side x Thus passing through the formulaObtaining work done by external force to deform the model unit in the x direction, wherein epsilon x Representing the strain of the model element in the x-direction, sigma x Representing the positive stress of the model element in the x-direction; the deformation in the y direction has no external force, so that no work is done, and therefore, all work done in the x direction is converted into strain energy, and similarly, all work done by the external force for the deformation of the model unit in the y direction and all work done by the external force for the deformation of the model unit in the z direction are converted into strain energy.
By the formulaObtaining work done by external force to deform the model unit in y direction, wherein epsilon y Representation modelStrain, sigma, of the cell in the y-direction y Representing the positive stress of the model element in the y-direction;
by the formulaObtaining work done by external force to deform the model unit in the z direction, wherein epsilon z Representing the strain of the model element in the z-direction, sigma z Representing the positive stress of the model element in the z direction;
by the formula w=w x +W y +W z Obtaining the total work of the external force on the model unit in the three-dimensional direction;
the total work of the model unit in three-dimensional deformation is converted into the strain energy of the system, so that the strain energy formula of the system is constructed
Wherein U is 0x Represents strain energy per unit volume formed by deformation of the model unit in the x direction, U 0y Represents strain energy per unit volume formed by deformation of the model element in the y direction, U 0z Represents strain energy per unit volume formed by deformation of the model unit in the z direction, U 0 Representing the total strain potential energy density of the model element.
The deformation of the object can be decomposed into two parts according to the plastic theory, wherein one part is the change of volume and the other part is the change of shape, so that the total strain potential energy density of the model unit is the sum of the strain energy stored in the unit volume due to the change of volume and the strain energy stored in the unit volume due to the change of shape.
By the formulaAcquiring strain energy stored in a unit volume due to a volume change, wherein σ 0 Representing the stress before deformation of the model element, i.e. the original stress, K is the bulk modulus and K is a constant.
By the formulaAcquiring strain energy stored in a unit volume due to a shape change, wherein σ 1 Representing the preset strain force, sigma, of the model unit in the x direction 2 Representing the preset strain force, sigma, of the model unit in the y direction 3 Representing the preset strain force of the model element in the z direction, G represents the shear modulus and G takes a constant. Sigma is added to 0 、σ 1 、σ 2 Sigma (sigma) 3 Is through FLAC 3D Results of software simulations, i.e., at FLAC 3D When the model units corresponding to surrounding rocks around the roadway are constructed on the software, the stress of the model units in all directions is simulated, so sigma 0 、σ 1 、σ 2 Sigma (sigma) 3 Are known constants, which can be passed through FLAC 3D The software simulation results may also be directly preset.
By the formula
And obtaining the total strain potential energy density of the model unit.
The analyzing the energy field according to the variable cloud image comprises the following steps: and (5) finding a potential energy gathering area according to the variable cloud chart, and reinforcing and supporting surrounding rocks around the roadway where the potential energy gathering area is located. The following are examples of 3 specific optimized supports.
Example one: and (3) optimizing the high-drainage roadway support of a mining group 17236, wherein 17236 is the number of the high-drainage roadway of the mining group. Fig. 3 is a 17236 high-suction roadway potential energy aggregation distribution chart, and table 1 is a roadway surrounding rock potential energy aggregation distribution chart.
Table 1 roadway surrounding rock potential energy gathering and distributing table
| Position of | Potential energy density mean (J/m) 3 ) | Potential energy distribution range (m) 2 ) |
| Left upper | 351310 | 4.6 |
| Right upper | 340731 | 14 |
| Left shoulder nest | 264100 | 1.7 |
As shown in the graph, strain potential energy is mainly accumulated at two sides and a left shoulder nest of a high-suction roadway and is asymmetrically distributed, the potential energy range of a right side is larger than that of the left side, but the potential energy value of the left side is slightly larger than that of the right side, particularly, compared with the top plate of the right shoulder nest, the top plate of the left side shoulder nest has obvious potential energy accumulation, and the potential energy accumulation area possibly has a destructive effect on an anchor rod at the left side shoulder nest, so that the anchor rod breakage phenomenon at the left side shoulder nest of a roadway behind is to be reinforced and examined in the roadway tunneling process, broken anchors are found to be timely repaired, and meanwhile, the support of the right side is reinforced.
Example two: analysis and countermeasure for broken anchor of shoulder pit of rectangular roadway side part
And judging the energy distribution rule around the rectangular roadway by using an energy method. After the tunnel is excavated, the energy concentration is generated. The level of concentration is related to the stress of the surrounding rock and the displacement generated, and the higher the energy concentration is, the higher the stress level at the position is, the larger the displacement of the surrounding rock is, and the greater the damage to the supporting structure and the surrounding rock is.
Fig. 4 is an energy distribution diagram after the arch-shaped roadway is excavated, and fig. 5 is an energy distribution diagram after the rectangular roadway is excavated, it can be clearly seen that after the arch-shaped roadway is excavated, the energy concentration is mainly distributed at two bottom corners of the roadway, which is 8 times of a roadway bottom plate and about 2 times of a shoulder pit. The position with the greatest energy concentration of the rectangular roadway is formed in the shoulder pit and the two sides of the roadway, particularly the side of the rectangular roadway along the air, and the shoulder pit has obviously high energy, and the energy value is 1.5 multiplied by 105, which is the main reason for causing the anchor rod of the shoulder pit of the rectangular roadway to break. Countermeasures that can be taken:
1) The material strength is improved. Because the shearing strength of the deformed steel bar anchor rod is low, the stress damage of surrounding rock cannot be effectively resisted, and the deformed steel bar anchor rod can be replaced by a steel strand short anchor rope. The shearing strength of the steel strand anchor cable is 6.2 times of that of the anchor rod, and the steel strand anchor cable can meet the strength requirement under the current ground stress and engineering conditions.
2) Grouting along the empty coal pillar, increasing the self strength of the coal pillar, improving the deformation resistance of surrounding rock and enhancing the stability of the supporting structure.
3) The stress distribution of the arch roadway has obvious advantages compared with that of a rectangular roadway, the arch roadway can be used as a reference along the gob-side entry, the top plate forms an approximate arch structure, and the stress concentration level is reduced.
Example three: 2.5m shield tunnel anchor bolt support energy method analysis
Fig. 6 is a schematic diagram of the situation of potential energy concentration without support, fig. 7 is a schematic diagram of the situation of potential energy concentration with 2 anchor bolts, fig. 8 is a schematic diagram of the situation of potential energy concentration with 4 anchor bolts, and fig. 9 is a schematic diagram of the situation of potential energy concentration with 8 anchor bolts. The middle wing-like area in fig. 6 to 9 is a low potential energy accumulation area, the fan-shaped areas on two sides of the wing-like area are high potential energy accumulation areas, the potential energy accumulation distribution of four working conditions is basically consistent, a closed low potential energy accumulation area is formed around a roadway, strain potential energy is mainly accumulated at 1.3m positions in two sides of the roadway, and when surrounding rock of the side of the roadway is broken, side wall shoring should be timely reinforced to avoid side wall instability.
Through the technical scheme, the invention utilizes FLAC 3D The software performs triple integration on the total strain potential energy density of the model unit to obtain the strain energy of the system, so as to obtain the current energy field, and the current energy field passes through FLAC 3D Graphic display with built-in softwareThe module displays the current energy field to form a variable cloud chart, analyzes the energy field according to the variable cloud chart, provides an energy field analysis method for surrounding rocks around the roadway, is beneficial to analyzing the surrounding rocks of the roadway, guides site construction and optimizes roadway support design.
Example 2
The invention also provides a device for realizing the energy field analysis of surrounding rocks around a roadway, which comprises: strain energy acquisition module of system for FLAC 3D Constructing a model unit corresponding to surrounding rock around a roadway on software, and carrying out triple integration on the total strain potential energy density of the model unit to obtain the strain energy of the system;
a variable cloud image acquisition module for assigning strain energy of the system to FLAC 3D A cohesion display module built in the software displays a variable cloud picture formed by the current energy field;
and the analysis module is used for analyzing the energy field according to the variable cloud image.
Specifically, the strain energy of the system is obtained by the following steps:
by the formulaObtaining work done by external force to deform the model unit in the x direction, wherein epsilon x Representing the strain of the model element in the x-direction, sigma x Representing the positive stress of the model element in the x-direction;
by the formulaObtaining work done by external force to deform the model unit in y direction, wherein epsilon y Representing the strain of the model element in the y-direction, σ y Representing the positive stress of the model element in the y-direction;
by the formulaObtaining work done by external force to deform the model unit in the z direction, wherein epsilon z Representing the strain of the model element in the z-direction, sigma z Representing the positive stress of the model element in the z direction;
by the formula w=w x +W y +W z Obtaining the total work of the external force on the model unit in the three-dimensional direction;
the total work of the model unit in three-dimensional deformation is converted into the strain energy of the system, so that the strain energy formula of the system is constructed
Wherein U is 0x Represents strain energy per unit volume formed by deformation of the model unit in the x direction, U 0y Represents strain energy per unit volume formed by deformation of the model element in the y direction, U 0z Represents strain energy per unit volume formed by deformation of the model unit in the z direction, U 0 Representing the total strain potential energy density of the model element.
More specifically, the deformation of the object obtained by the plastic theory can be decomposed into two parts, one part is the change of volume and the other part is the change of shape, so that the total strain potential energy density of the model unit is the sum of the strain energy stored in the unit volume due to the change of volume and the strain energy stored in the unit volume due to the change of shape.
More specifically, by the formulaAcquiring strain energy stored in a unit volume due to a volume change, wherein σ 0 Representing the stress before deformation of the model element, i.e. the original stress, K is the bulk modulus and K is a constant.
More specifically, by the formulaAcquiring strain energy stored in a unit volume due to a shape change, wherein σ 1 Representing the preset strain force, sigma, of the model unit in the x direction 2 Representing the preset strain force, sigma, of the model unit in the y direction 3 Representing model elements inThe preset strain force in the z direction, G represents the shear modulus and G takes a constant.
More specifically, by the formula
And obtaining the total strain potential energy density of the model unit.
Specifically, the analysis module is further configured to: and (5) finding a potential energy gathering area according to the variable cloud chart, and reinforcing and supporting surrounding rocks around the roadway where the potential energy gathering area is located.
More specifically, the larger the energy concentration of the potential energy accumulation area is, the higher the stress level of the corresponding position is, the larger the displacement of the surrounding rock is, and the larger the damage to the surrounding rock and the supporting structure of the surrounding rock is.
More specifically, the method for reinforcing and supporting surrounding rocks around the roadway where the potential energy gathering area is located comprises the following steps: the material strength of the supporting structure of surrounding rocks around the roadway is improved, the surrounding rocks around the roadway are injected along the empty coal pillar, or the supporting roof plate of the surrounding rocks around the roadway is designed into an arch structure.
The above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (5)
1. A method of performing an energy field analysis of surrounding rock around a roadway, the method comprising: FLAC (FLAC) 3D Constructing model units corresponding to surrounding rocks around a roadway on software, performing triple integration on the total strain potential energy density of the model units to obtain strain energy of the system, and assigning the strain energy of the system to the FLAC 3D The internal cohesion display module of the software displays the current energy field to form a variable cloud picture, and the energy field is divided according to the variable cloud pictureSeparating out;
the strain energy of the system is obtained by the following steps:
by the formulaThe external force is acquired at +.>Work done by directional deformation, wherein +.>Representing model element in->Strain in the direction, ++>Representing model element in->Positive stress in the direction;
by the formulaThe external force is acquired at +.>Work done by directional deformation, wherein +.>Representing model element in->Strain in the direction, ++>Representing model element in->Positive stress in the direction;
by the formulaThe external force is acquired at +.>Work done by directional deformation, wherein +.>Representing model element in->Strain in the direction, ++>Representing model element in->Positive stress in the direction;
by the formulaObtaining the total work of the external force on the model unit in the three-dimensional direction;
the total work of the model unit in three-dimensional deformation is converted into the strain energy of the system, so that the strain energy formula of the system is constructed
Wherein->Indicating that the model unit is->Strain energy per unit volume formed by directional deformation, +.>Indicating that the model unit is->Strain energy per unit volume formed by directional deformation, +.>Indicating that the model unit is->Strain energy per unit volume formed by directional deformation, +.>Representing the total strain potential energy density of the model unit;
the deformation of the object can be decomposed into two parts according to the plastic theory, wherein one part is the change of volume and the other part is the change of shape, so that the total strain potential energy density of the model unit is the sum of the strain energy stored in the unit volume due to the change of volume and the strain energy stored in the unit volume due to the change of shape;
by the formulaAcquiring strain energy stored in a unit volume due to volume change, wherein +.>Representing the stress before deformation of the model element, i.e. the original stress,/->Is of bulk modulus and>taking a constant;
by the formulaAcquiring strain energy stored in a unit volume due to shape change, wherein +.>Representing model element in->Preset strain force in direction, +.>Representing model element in->Preset strain force in direction, +.>Representing model element in->Preset strain force in direction, +.>Represents the shear modulus and>taking a constant;
by the formulaAnd obtaining the total strain potential energy density of the model unit.
2. A method of performing energy field analysis of surrounding rock around a roadway according to claim 1, wherein the analyzing the energy field according to the variable cloud pattern comprises: and (5) finding a potential energy gathering area according to the variable cloud chart, and reinforcing and supporting surrounding rocks around the roadway where the potential energy gathering area is located.
3. The method for realizing the energy field analysis of surrounding rocks around a roadway according to claim 2, wherein the larger the energy concentration of the potential energy accumulation area is, the higher the stress level of the corresponding position is, the larger the displacement of the surrounding rocks is, and the larger the damage to the surrounding rocks and the supporting structure of the surrounding rocks is.
4. The method for realizing the energy field analysis of surrounding rocks around a roadway according to claim 2, wherein the method for reinforcing and supporting the surrounding rocks around the roadway at the position of the potential energy accumulation zone comprises the following steps: the material strength of the supporting structure of surrounding rocks around the roadway is improved, the surrounding rocks around the roadway are injected along the empty coal pillar, or the supporting roof plate of the surrounding rocks around the roadway is designed into an arch structure.
5. An apparatus for performing energy field analysis of surrounding rock around a roadway using the method of any one of claims 1-4, the apparatus comprising: strain energy acquisition module of system for FLAC 3D Constructing a model unit corresponding to surrounding rock around a roadway on software, and carrying out triple integration on the total strain potential energy density of the model unit to obtain the strain energy of the system;
a variable cloud image acquisition module for assigning strain energy of the system to FLAC 3D A cohesion display module built in the software displays a variable cloud picture formed by the current energy field;
and the analysis module is used for analyzing the energy field according to the variable cloud image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110728033.8A CN113392561B (en) | 2021-06-29 | 2021-06-29 | Method and device for realizing energy field analysis of surrounding rock around roadway |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110728033.8A CN113392561B (en) | 2021-06-29 | 2021-06-29 | Method and device for realizing energy field analysis of surrounding rock around roadway |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113392561A CN113392561A (en) | 2021-09-14 |
| CN113392561B true CN113392561B (en) | 2024-02-23 |
Family
ID=77624427
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110728033.8A Active CN113392561B (en) | 2021-06-29 | 2021-06-29 | Method and device for realizing energy field analysis of surrounding rock around roadway |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113392561B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106053331A (en) * | 2016-05-27 | 2016-10-26 | 华北理工大学 | Method for migrating and quantitating roadway surrounding rock strain energy under action of excavation unloading |
-
2021
- 2021-06-29 CN CN202110728033.8A patent/CN113392561B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106053331A (en) * | 2016-05-27 | 2016-10-26 | 华北理工大学 | Method for migrating and quantitating roadway surrounding rock strain energy under action of excavation unloading |
Non-Patent Citations (9)
| Title |
|---|
| 在围压作用下硬包体模型若干力学性质的研究;陈宇坤, 聂永安, 崔晓峰, 王文, 王进忠;华北地震科学;19990630(02);全文 * |
| 大断面采动剧烈影响煤巷变形破坏机制与控制技术;何富连;张广超;;采矿与安全工程学报;20160515(03);全文 * |
| 岩石变形破坏过程中的能量传递和耗散研究;赵忠虎;谢和平;;四川大学学报(工程科学版);20080315(02);全文 * |
| 巷道开挖围岩能量释放与偏应力应变能生成的分析计算;潘岳;王志强;吴敏应;;岩土力学;20070420(04);全文 * |
| 李家壕矿软岩回采巷道围岩稳定性分析与支护技术研究;王烨;;煤炭工程;20170828(S2);全文 * |
| 窄煤柱综放回采巷道围岩稳定性分析及控制技术;王虎胜;黄肖;席朝东;杨亚峰;贺冲冲;刘斌慧;;中国煤炭;20150922(09);全文 * |
| 采动巷道稳定性数值模拟分析;贾照远;尚剑飞;刘金辉;;建井技术;20191225(06);全文 * |
| 采用边坡稳定性强度折减法分析弹性应变能突变判据;柴红保;曹平;林杭;赵延林;;中南大学学报(自然科学版);20090826(04);全文 * |
| 高强度开采工作面煤岩灾变的推进速度效应分析;杨胜利;王兆会;蒋威;杨敬虎;;煤炭学报;20160315(03);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113392561A (en) | 2021-09-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yazdani et al. | Displacement-based numerical back analysis for estimation of rock mass parameters in Siah Bisheh powerhouse cavern using continuum and discontinuum approach | |
| Liu et al. | Full 3D modelling for effects of tunnelling on existing support systems in the Sydney region | |
| Liu et al. | Experimental study and numerical simulation on dynamic soil‐structure interaction under earthquake excitations | |
| Fraldi et al. | Evaluation of impending collapse in circular tunnels by analytical and numerical approaches | |
| Hsiung | A case study on the behaviour of a deep excavation in sand | |
| Li et al. | Investigation on passive pile responses subject to adjacent tunnelling in anisotropic clay | |
| Kulatilake et al. | Investigation of stability of a tunnel in a deep coal mine in China | |
| Li et al. | Study on excavation sequence of pilot tunnels for a rectangular tunnel using numerical simulation and field monitoring method | |
| Ding et al. | Stability of large underground caverns excavated in layered rock masses with steep dip angles: a case study | |
| CN109271738A (en) | It is a kind of for obtaining the Method for Numerical Inversion of roadway surrounding rock Weibull distribution parameter | |
| Do et al. | Verifying discontinuous deformation analysis simulations of the jointed rock mass behavior of shallow twin mountain tunnels | |
| CN109492262A (en) | A method of utilizing numerical simulation analysis non-uniform Distribution crack Drift stability | |
| Yang et al. | Numerical study of group effects for pile groups in sands | |
| Yang et al. | Time and technique of rehabilitation for large deformation of tunnels in jointed rock masses based on FDM and DEM numerical modeling | |
| CN105550441A (en) | Engineering rock mass fracture degradation numerical simulation method based on continuous medium | |
| CN111814234A (en) | Construction treatment method for broken rock mass flood discharge tunnel | |
| Liu et al. | Reliability assessment on stability of tunnel-type anchorages | |
| Wang et al. | Seismic responses and damage mechanisms of the structure in the portal section of a hydraulic tunnel in rock | |
| CN116703244B (en) | Mining subsidence area treatment effect and comprehensive evaluation method | |
| Li et al. | Layout and length optimization of anchor cables for reinforcing rock wedges | |
| Zhang et al. | A simplified analytical solution for deformation behavior of existing tunnels subjected to influences of landslides | |
| Zhang et al. | Effects of soil arching on behavior of composite pile supporting foundation pit | |
| CN113392561B (en) | Method and device for realizing energy field analysis of surrounding rock around roadway | |
| Yang et al. | The influence of different excavation methods on deep foundation pit and surrounding environment | |
| Lei et al. | Study on train vibration response and cumulative deformation of double arch tunnel in Kast foundation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |