CN115169752B - Coal mine roof initial pressure step pitch prediction method based on virtual work theory - Google Patents
Coal mine roof initial pressure step pitch prediction method based on virtual work theory Download PDFInfo
- Publication number
- CN115169752B CN115169752B CN202211076186.XA CN202211076186A CN115169752B CN 115169752 B CN115169752 B CN 115169752B CN 202211076186 A CN202211076186 A CN 202211076186A CN 115169752 B CN115169752 B CN 115169752B
- Authority
- CN
- China
- Prior art keywords
- top plate
- load
- virtual
- plastic
- roof
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Data Mining & Analysis (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Mathematical Optimization (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Business, Economics & Management (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Mining & Mineral Resources (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Evolutionary Biology (AREA)
- Marine Sciences & Fisheries (AREA)
- Game Theory and Decision Science (AREA)
- Probability & Statistics with Applications (AREA)
- Health & Medical Sciences (AREA)
- Algebra (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
The invention belongs to the field of mine pressure and rock stratum control, and discloses a method for predicting the initial pressure step pitch of a coal mine roof based on a virtual work theory. The method is simple and convenient, and can effectively predict the initial pressure step distance of the top plate, guide underground supporting and protecting work and reduce the damage of personnel and equipment caused by initial pressure.
Description
Technical Field
The invention belongs to the field of mine pressure and rock stratum control, and particularly relates to a method for predicting the primary pressure step pitch of a coal mine roof based on a virtual work theory.
Background
The basic roof is the roof located directly above or above the coal seam, also known as the old roof. When the working face is pushed forward from the open-off cut, the exposed area of the top plate is enlarged, the direct roof is collapsed to fill the goaf, and the old roof is still supported on the two coal walls completely to form a double-support beam component. When the span of the plate girder is increased to a certain range along with the advancing of the working face, the old roof is collapsed due to the self weight of the old roof and the action of the overlying rock stratum. At this time, the working face is no longer under the shield of the old roof, the roof rapidly sinks and is broken, and the first large-area caving of the old roof of the goaf is generally called primary caving. The pressure of the working face is increased due to the primary collapse of the old roof.
The roof accident is one of the main accidents in the underground coal mine, the damage is great, and the initial pressure is stronger and more harmful than the periodic pressure, so that the supporting and pressure reducing work can be reasonably carried out by effectively predicting the initial pressure step of the roof, and the accidents of personnel and equipment are reduced.
Disclosure of Invention
In view of the above, the invention aims to provide a method for predicting the initial pressure step pitch of a coal mine roof based on the virtual work theory, which is simple and convenient, can effectively predict the initial pressure step pitch of the roof, guides underground supporting and protecting work, and reduces the damage to personnel and equipment caused by the initial pressure.
In order to achieve the purpose of the invention, the technical scheme is as follows:
a coal mine roof initial pressure step pitch prediction method based on a virtual work theory is characterized by comprising the following steps of:
s1, determining the length of a working face, the width of a transportation gateway and a return air gateway, the average thickness of a top plate and the mechanical characteristics of top plate rocks according to the actual condition of the working face, wherein the top plate reaches a plastic limit state along with the propulsion of the working face in an underground coal mine, and a plastic hinge line in a four-slope top form appears;
s2, calculating virtual work done by external force:
virtual work done by external forceWVolume and load of "four-pitched roof" obtained by the cut-and-fill methodqThe product of (a) and (b) yields:
wherein the content of the first and second substances,equal to the volume of the destructive mechanism in the form of virtual displacement;
in the formula (I), the compound is shown in the specification,
W: the virtual work done for the external force is,
δ: the deflection of the center of the top plate,
a: the sum of the length of the working face and the width of the two gate grooves,
b: the step pitch is pressed for the first time,
q: covering a load on the top plate;
s3, calculating total virtual dissipation energy of internal force:
the total internal force virtual dissipation energy is the work of the limit bending moment at the top of the four-slope top and the plastic hinge line of the plate frame on the relative corner, namely:
in the formula (I), the compound is shown in the specification,
T: the internal force is always deficient and energy is dissipated,
k: the number of the plastic twisted wires is large,
l: the length of the plastic twisted wire is as long as,
φ: the relative angle of rotation is such that,
M p : plastic twisted wire unit length ultimate bending moment;
solving the relative rotation angle of the plastic twisted line by using space geometry:
substituting the above parameters to obtain total internal force deficiency energy dissipationTThe calculation formula of (c) is:
s4, obtaining the corresponding relation between the ultimate load and the plate size in the plastic limit state, and the physical and mechanical properties according to the condition that the virtual work done by the external force is equal to the total virtual dissipation of the internal force;
Has an average thickness ofhThe top plate of (a) is provided,,is obtained by substituting the yield limit of the roof strata;
The upper formula is the load under the plastic limit stateqRepresentative value of twisted wiretThe corresponding relationship of (a) to (b),tbyqTaking the minimum to obtain:,
In the formula (I), the compound is shown in the specification,the limit load is the limit load in the plastic limit state;
s5, calculating the step pitch of the first pressure:
the length of the plate is long along with the advance of the working face under the coal minebWhen the critical value is reached, the top plate enters a plastic limit state, and the load on the top plate is taken as the limit load, namelyq=q Pole The load on the top plate, the sum of the length of the working face and the width of the two crossroads, the physical and mechanical properties of the rock stratum of the top plate and the average thickness of the top plate are knownhThe length of the plate can be calculated by the following formulabValue of (2), plate lengthbNamely the initial pressure step distance, the calculation formula is as follows:
wherein, the first and the second end of the pipe are connected with each other,,hthe average thickness of the top plate is,qan overburden is applied to the roof.
The beneficial effects of the invention are: the invention provides a method for predicting the initial pressure step of a coal mine roof based on a virtual work theory, which comprises the steps of firstly determining the length of a working face, the width of a transportation gate way and a return air gate way, the average thickness of the roof and the mechanical characteristics of roof rocks according to the actual condition of the working face, then determining the load of an overlying rock layer supported by the roof of the working face according to a key layer theory, leading the roof to reach a plastic limit state along with the propulsion of the working face, generating a plastic hinge line in a four-slope roof form, and predicting the initial pressure step of the roof according to the corresponding relation between the load when the roof enters the plastic limit state and the size of the roof when the plastic hinge line appears. The method is simple and convenient, and can effectively predict the initial pressure step distance of the top plate, guide underground supporting and protecting work and reduce the damage of personnel and equipment caused by initial pressure.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this application, are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification. In the drawings:
FIG. 1 is a top view of a top plate in a plastic limit state, wherein the BA direction is a working surface advancing direction;
FIG. 2 is a schematic view of the volume of the four-slope top obtained by the cutting and filling method;
FIG. 3 is a perspective view of the top plate in a plastic limit state;
FIG. 4 shows the relative rotation angles of plastic twisted lines AE, BE, CF, DFθSolving an auxiliary graph;
FIG. 5 shows the relative rotation angle of the plastic hinge line EFαA schematic view;
FIG. 6 shows the opposite rotation angle of plastic hinge line AB and CDφA schematic diagram;
FIG. 7 shows the relative rotation angle of the plastic twisted lines AC and BDβSchematic representation.
Detailed Description
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present invention will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.
In order to make the technical solutions better understood by those skilled in the art, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only partial embodiments of the present application, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
As shown in fig. 1 to 7, a method for predicting the initial pressure step of a coal mine roof based on the virtual work theory, which assumes that the uniform thickness of the roof does not affect the stress analysis of the plate, comprises the following steps:
s1, determining the length of a working face, the width of a transportation gateway and a return air gateway, the average thickness of a top plate and the mechanical characteristics of top plate rocks according to the actual condition of the working face, wherein the top plate reaches a plastic limit state along with the propulsion of the working face in an underground coal mine, and a plastic hinge line in a four-slope top form appears;
s2, calculating virtual work done by external force:
virtual work done by external forceWVolume and load of "four-pitched roof" obtained by the cut-and-fill methodqThe product of (a) yields:
wherein, the first and the second end of the pipe are connected with each other,equal to the volume of the destructive mechanism in the form of virtual displacement;
the cut-and-fill method is characterized in that the volume of the triangular prism is subtracted by the volume of two tetrahedrons to obtain the volume of the' four-slope topVirtual work done by external force;
In the formula (I), the compound is shown in the specification,
δ: the flexibility of the center of the top plate,
a: the sum of the length of the working surface and the width of the two crossroads (the transportation crossroad and the return air crossroad),
b: the step pitch is pressed for the first time,
q: covering a load on the top plate;
s3, calculating total internal force virtual dissipation energy:
the total virtual dissipation energy of the internal force is the work of the limit bending moment at the tops of the four slope tops and the plastic hinge lines of the plate frames on the relative corner, namely:
in the formula (I), the compound is shown in the specification,
T: the general deficiency of internal force consumes the energy to dissipate,
k: the number of the plastic twisted wires is large,
l: the length of the plastic twisted wire is as long as,
φ: the relative angle of rotation is such that,
M p : plastic twisted wire unit length ultimate bending moment;
the relative rotation angles of the plastic twisted lines are solved by using the space geometry, and the actual relative rotation angles of the plastic twisted lines AE, BE, CF and DF are shown in figure 3θEstablishing a coordinate system as shown in FIG. 4,θthe included angle between the face ABE 'and the face BDF' E 'is obtained, the normal vector of the face ABE' and the normal vector of the face BDF 'E' are respectively obtained, and the included angle between the normal vectors and the dihedral angle is obtained through solutionThe lengths of the plastic twisted lines AE, BE, CF and DF being equal toTherefore, the virtual work done by the four plastic twisted wires is:
plastic hinge line EF relative angleαCalculated from a cross-section perpendicular to the direction of advancement of the working surface, according to the infinitesimal nature of the virtual displacementVirtual work done by plastic hinge EF;
Plastic twisted line AB, CD relative cornerφDerived from a section parallel to the direction of advancement of the working surface, and having the same principleDeficiency of the Qi;
In the same way, the relative rotation angle of the plastic twisted line AC and BDDeficiency of the Qi;
S4, obtaining the corresponding relation between the ultimate load and the plate size in the plastic limit state, and the physical and mechanical properties according to the condition that the virtual work done by the external force is equal to the total virtual dissipation of the internal force;
Because of infinitesimal virtual displacement, the deflection of the top plate center is infinitely close to zero, so that the following conditions are provided:
average thickness ofhThe top plate of (a) is provided,substituting for the yield limit of the roof strata to obtain:
a、bin the known manner, it is known that,q、tunknown, the above formula is the load under the plastic limit stateqRepresentative value of twisted wiretThe corresponding relationship of (2).
the upper formula is the corresponding relation between the lower limit load of the plastic limit state and the size, physical property and mechanical property of the top plate, the upper load of the top plate is known in the underground coal mine, and the plate length is increased along with the advance of the working surfacebWhen the critical value is reached, the top plate enters a plastic limit state, and the corresponding relation in the formula is met.
S5, calculating the initial pressure step:
the length of the plate is long along with the advance of the working face under the coal minebWhen the critical value is reached, the top plate enters a plastic limit state, and the load on the top plate is taken as the limit load, namelyq=q Pole(s) The physical and mechanical properties of the roof strata are known as the sum of the overburden load on the roof, the length of the working face and the width of the two crossroadsAverage thickness of the top platehThe length of the plate can be calculated by the above formulabValue of (2), plate lengthbNamely the initial pressure step.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement, component separation or combination and the like made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.
Claims (1)
1. A coal mine roof initial pressure step pitch prediction method based on a virtual work theory is characterized by comprising the following steps of:
s1, determining the length of a working face, the width of a transportation gateway and a return air gateway, the average thickness of a top plate and the mechanical characteristics of top plate rocks according to the actual condition of the working face, wherein the top plate reaches a plastic limit state along with the propulsion of the working face in an underground coal mine, and a plastic hinge line in a four-slope top form appears;
s2, calculating the virtual work done by the external force:
virtual work done by external forceWVolume and load of "four-pitched roof" obtained by the cut-and-fill methodqThe product of (a) yields:
wherein the content of the first and second substances,equal to the volume of the destructive mechanism in the virtual displacement mode;
in the formula (I), the compound is shown in the specification,
W: the virtual work done by the external force is provided,
δ: the flexibility of the center of the top plate,
a: the sum of the length of the working face and the width of the two gate grooves,
b: the step pitch is pressed for the first time,
q: covering a load on the top plate;
s3, calculating total internal force virtual dissipation energy:
the total internal force virtual dissipation energy is the work of the top of the four-slope top and the extreme bending moment at the plastic hinge line of the plate frame on the relative corner, and the total internal force virtual dissipation energyTThe calculation formula of (c) is:
s4, obtaining the corresponding relation between the ultimate load and the plate size in the plastic limit state, and the physical and mechanical properties according to the condition that the virtual work done by the external force is equal to the total virtual dissipation of the internal force;
Has an average thickness ofhThe top plate of (a) is provided,,is obtained by substituting the yield limit of the roof strata;
The upper formula is the load under the plastic limit stateqRepresentative value of twisted wiretThe corresponding relationship of (a) to (b),tbyqTaking the minimum to obtain:,
In the formula (I), the compound is shown in the specification,the limit load is the limit load in the plastic limit state;
s5, calculating the step pitch of the first pressure:
the length of the plate is long along with the advance of the working face under the coal minebWhen the critical value is reached, the top plate enters a plastic limit state, and the load on the top plate is taken as the limit load, namelyq=q Pole(s) The load on the top plate, the sum of the length of the working surface and the width of the two crossroads, the physical and mechanical properties of the rock stratum of the top plate and the average thickness of the top plate are knownhThe length of the plate can be calculated by the following formulabValue of (2), plate lengthbNamely the initial pressure step distance, the calculation formula is as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211076186.XA CN115169752B (en) | 2022-09-05 | 2022-09-05 | Coal mine roof initial pressure step pitch prediction method based on virtual work theory |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211076186.XA CN115169752B (en) | 2022-09-05 | 2022-09-05 | Coal mine roof initial pressure step pitch prediction method based on virtual work theory |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115169752A CN115169752A (en) | 2022-10-11 |
CN115169752B true CN115169752B (en) | 2022-12-23 |
Family
ID=83480599
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211076186.XA Active CN115169752B (en) | 2022-09-05 | 2022-09-05 | Coal mine roof initial pressure step pitch prediction method based on virtual work theory |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115169752B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116988831B (en) * | 2023-09-27 | 2023-12-12 | 太原理工大学 | Roof fracture form-based coal mine initial mining period I-shaped filling method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106285782A (en) * | 2016-08-30 | 2017-01-04 | 中国矿业大学(北京) | Bump method for early warning under a kind of complicated geological Environmental effect and system |
CN111305847A (en) * | 2019-12-06 | 2020-06-19 | 中国人民解放军军事科学院国防工程研究院工程防护研究所 | Dynamic ultimate strength control method for ultra-shallow-buried underground excavation construction of urban underground engineering |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102182453B (en) * | 2011-03-03 | 2014-01-29 | 中国石油集团钻井工程技术研究院 | Borehole wall sloughing analysis method |
CN114282362A (en) * | 2021-12-17 | 2022-04-05 | 晋能控股煤业集团有限公司 | Coal mine old goaf water inrush monitoring and judging method |
CN114266181B (en) * | 2021-12-24 | 2024-04-12 | 太原理工大学 | Vector limit balance safety coefficient calculation method based on projection direction extreme value |
-
2022
- 2022-09-05 CN CN202211076186.XA patent/CN115169752B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106285782A (en) * | 2016-08-30 | 2017-01-04 | 中国矿业大学(北京) | Bump method for early warning under a kind of complicated geological Environmental effect and system |
CN111305847A (en) * | 2019-12-06 | 2020-06-19 | 中国人民解放军军事科学院国防工程研究院工程防护研究所 | Dynamic ultimate strength control method for ultra-shallow-buried underground excavation construction of urban underground engineering |
Also Published As
Publication number | Publication date |
---|---|
CN115169752A (en) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115169752B (en) | Coal mine roof initial pressure step pitch prediction method based on virtual work theory | |
Guo et al. | A new method of predicting the height of the fractured water-conducting zone due to high-intensity longwall coal mining in China | |
Batugin et al. | Combined support mechanism of rock bolts and anchor cables for adjacent roadways in the external staggered split-level panel layout | |
Frith | A holistic examination of the load rating design of longwall shields after more than half a century of mechanised longwall mining | |
Aghababaei et al. | Risk analysis and prediction of floor failure mechanisms at longwall face in Parvadeh-I coal mine using rock engineering system (RES) | |
US3508368A (en) | Building panel having tapered countersunk end portion and method of erecting the same | |
Mark et al. | An empirical method for the design of chain pillars for longwall mining | |
Noroozı et al. | A model for determining the breaking characteristics of immediate roof in longwall mines | |
Ma et al. | Case study and design of steel set support for aged belt entry rehabilitation | |
Luo et al. | Cispm–ms: a Tool to predict surface subsidence and to study interactions associated with multi–seam mining Operations | |
CN109024988A (en) | A kind of double-layer plate shear wall structure put more energy into close rib being easily installed | |
Lind | Coal pillar extraction experiences in New South Wales | |
CN116988831B (en) | Roof fracture form-based coal mine initial mining period I-shaped filling method | |
Dahl et al. | Rock mechanics elements of coal mine design | |
Wang et al. | Influence of lateral structure and combined support system of split-level entries in the thick seam | |
Iannacchione | Behavior of a coal pillar prone to burst in the southern Appalachian Basin of the United States | |
CN218264001U (en) | Construction joint filling structure | |
Evanek et al. | Controlling crosscut damage in response to excessive levels of horizontal stress: case study at the Subtropolis Mine, Petersburg, OH | |
CN219535560U (en) | Protection mechanism for hot galvanizing bridge | |
CN104912590B (en) | Separated filling entry retaining method for sandstone roof medium-thickness coal seam | |
Kovalevska et al. | Geomechanics of overworked mine working support resistance in the laminal massif of soft rocks | |
CN214785376U (en) | Enhancement mode SBS waterproofing membrane | |
CN220318942U (en) | Assembled side fascia with big window entrance to a cave | |
Mandingaisa et al. | Stable span re-design and support optimization for a shallow hard-rock bord-and-pillar mine on the Great Dyke: the case of Unki Mine, Zimbabwe | |
Luo et al. | Effects of dewatering |
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 |