CN106643610B - A kind of calculation method of coal body plastic zone width - Google Patents
A kind of calculation method of coal body plastic zone width Download PDFInfo
- Publication number
- CN106643610B CN106643610B CN201611247319.XA CN201611247319A CN106643610B CN 106643610 B CN106643610 B CN 106643610B CN 201611247319 A CN201611247319 A CN 201611247319A CN 106643610 B CN106643610 B CN 106643610B
- Authority
- CN
- China
- Prior art keywords
- coal body
- coal
- plastic zone
- zone width
- stress
- 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
- 239000003245 coal Substances 0.000 title claims abstract description 111
- 238000004364 calculation method Methods 0.000 title claims abstract description 17
- 238000005065 mining Methods 0.000 claims abstract description 29
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 239000000463 material Substances 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 description 5
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/22—Fuels; Explosives
- G01N33/222—Solid fuels, e.g. coal
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Food Science & Technology (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of calculation method of coal body plastic zone width, comprising: establishes coal body plastic zone width computation model;Coal body mining height under current mining condition, buried depth, face guard are obtained to the Support Resistance of coal wall, coal seam and roof and floor coefficient of friction, the factor of stress concentration, internal friction angle, cohesive force, bulk density;Using history mining condition as reference, different Lode's parameters are taken to calculate coal body plastic zone width, and the calculated result of different Lode's parameters is compared analysis with practical coal body plastic zone width, take Lode's parameters corresponding to minimal difference as the Lode's parameters in coal body plastic zone width computation model;The coal body plastic zone width that working face is newly exploited in model prediction is calculated using coal body plastic zone width.The method of the present invention introduces Lode's parameters μ to react the true stress of coal body, to improve the accuracy of prediction on the basis of D-P criterion.
Description
Technical field
The invention belongs to coal production technical field, in particular to a kind of calculation method of coal body plastic zone width.
Background technique
In coal production, coal body plastic zone width (i.e. abutment pressure distribution feature) is to roof collapse, roadway maintenance, coal
Column width etc. has great influence, and production safe and efficient to safeguard work face is of great significance.Therefore coal body plastic zone
It is always the content extremely paid close attention to both at home and abroad that width, which calculates, successively proposes many calculation formula, A.H Wilson theory, big plate are split
Gap theory, mole-coulomb criterion, D-P criterion etc..
However, the calculation formula that above-mentioned theory and criterion obtain respectively has advantage and disadvantage.A.H Wilson theory is according to Britain
The big probable value that matter condition and exploitation example take, therefore the versatility of its formula is restricted;The analysis model mistake of big plate crack theory
In simple, many influence factors are not taken into account, so that formula is distorted;Mole-coulomb criterion has clear physical concept and is applicable in
The advantages that facilitating is used widely, but has ignored the influence of intermediate principal stress, cannot explain well coal body surrender or
It destroys;D-P criterion is the further research to mole-coulomb criterion, is calculating side though considering the influence of intermediate principal stress
Just influence of the intermediate principal stress to yield function is but had ignored, therefore its calculated value always big probable value.
It can be seen from the above, the equal existing defects of the theoretical calculation of the above method, it is therefore necessary to propose a kind of more accurate property
Calculation method, improve existing drawback, lay the foundation for the accurate coal body plastic zone width that calculates, to guarantee working face safety
Efficiently production.
Summary of the invention
In view of the problems of the existing technology, the present invention provides a kind of calculation method of coal body plastic zone width.
Technical scheme is as follows:
A kind of calculation method of coal body plastic zone width, comprising:
Step 1 establishes coal body plastic zone width computation model;
Wherein, coal body mining height M, buried depth H, face guard to the Support Resistance p of coal wall, coal seam and roof and floor coefficient of friction f, answer
Power coefficient of concentration K is obtained from scene;Internal friction angleCohesive force c, bulk density γ experiment measure;Intermediate parametersμ is
Lode's parameters react the stress of coal body, | μ |≤1, material parameter α, k is according to internal friction angleCohesive force c is calculated;
Step 2 obtains coal body mining height M under current mining condition, buried depth H, face guard to Support Resistance p, the coal of coal wall
Layer and roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c, bulk density γ;
Step 3, according to internal friction angleMaterial parameter α, k is calculated with cohesive force c;
Step 4 takes different Lode's parameters μ, in conjunction with the coal body mining height M under known history mining condition, buried depth H, shield
Side board is to the Support Resistance p of coal wall, coal seam and roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c,
Bulk density γ calculates the coal body plastic zone width under history mining condition using coal body plastic zone width computation model;
Step 5, the coal body plastic zone width under the calculated history mining condition of more different Lode's parameters μ and history
The difference of practical coal body plastic zone width under mining condition, using Lode's parameters μ corresponding to minimal difference as coal body plasticity
Lode's parameters μ in sector width computation model;
Step 6 calculates the coal body plastic zone width that working face is newly exploited in model prediction using coal body plastic zone width.
The step 1, comprising:
Step 1-1, force analysis is carried out to cell cube in coal body, establishes the stress balance equation of cell cube;
In formula: what document generally believed coal wall is laterally working face, can regard the scope of freedom, therefore level suffered by coal body as
Stress σxTo be much smaller than vertical stress σy, σx、σyFirst principal stress suffered by cell cube and third principal stress are respectively indicated, M- is adopted
It is high;C- cohesive force;The coal seam f- and roof and floor coefficient of friction;
Step 1-2, the consideration intermediate principal stress σ parallel with working face inclined direction2Influence to yield function introduces anti-
The Lode's parameters μ of stress is reflected, | μ |≤1, intermediate principal stress σ is obtained according to plane strain problems2:
Enable σ1=-σx, σ3=-σy, then
Step 1-3, when coal body reaches yield situation, using D-P criterion;
In formula: the first stress invariant I1=σ1+σ2+σ3;
Second deviator of stress invariant
Step 1-4, cell cube first principal stress and third principal stress relational expression are obtained by D-P criterion;
Step 1-5, by relational expressionBring the stress balance equation of cell cube intoAnd combine boundary condition x=0, σx=p, p are Support Resistance of the face guard to coal wall, are obtained vertical
Stress σyExpression formula;
Step 1-6, maximal support pressure expression formula σ=K γ H according to coal body on interface, obtains coal body plastic zone
Width computation model:
The utility model has the advantages that the calculation method of coal body plastic zone width provided by the invention introduces Lip river on the basis of D-P criterion
Moral parameter μ is to react the true stress of coal body, to improve the accuracy of prediction.
Detailed description of the invention
Fig. 1 is the stress analysis schematic diagram of cell cube in the specific embodiment of the invention;
Fig. 2 is the calculation method flow chart of coal body plastic zone width in the specific embodiment of the invention.
Specific embodiment
The present invention is described in further detail With reference to embodiment.
Present embodiment provides a kind of calculation method of coal body plastic zone width as shown in Figure 2, comprising:
Step 1 establishes coal body plastic zone width computation model;
Wherein, coal body mining height M, buried depth H, face guard to the Support Resistance p of coal wall, coal seam and roof and floor coefficient of friction f, answer
Power coefficient of concentration K is obtained from scene;Internal friction angleCohesive force c, bulk density γ experiment measure;Intermediate parametersμ is
Lode's parameters react the stress of coal body, | μ |≤1, material parameter α, k is according to internal friction angleCohesive force c is calculated;
The step 1, comprising:
Step 1-1, assume coal body homogeneous, force analysis is carried out to cell cube in coal body, as shown in Figure 1, establishing cell cube
Stress balance equation;
In formula: what document generally believed coal wall is laterally working face, can regard the scope of freedom, therefore level suffered by coal body as
Stress σxTo be much smaller than vertical stress σy, σx、σyFirst principal stress suffered by cell cube and third principal stress are respectively indicated, M- is adopted
It is high;C- cohesive force;The coal seam f- and roof and floor coefficient of friction;
Step 1-2, the consideration intermediate principal stress σ parallel with working face inclined direction2Influence to yield function introduces anti-
The Lode's parameters μ of stress is reflected, | μ |≤1, intermediate principal stress σ is obtained according to plane strain problems2:
Enable σ1=-σx, σ3=-σy, then
Step 1-3, when coal body reaches yield situation, using D-P criterion;
In formula: the first stress invariant I1=σ1+σ2+σ3;
Second deviator of stress invariant
Step 1-4, cell cube first principal stress and third principal stress relational expression are obtained by D-P criterion;
Step 1-5, by relational expressionBring the stress balance equation of cell cube intoAnd combine boundary condition x=0, σx=p, p are Support Resistance of the face guard to coal wall, are obtained vertical
Stress σyExpression formula;
Step 1-6, maximal support pressure expression formula σ=K γ H according to coal body on interface, obtains coal body plastic zone
Width computation model:
By coal body plastic zone width computation model it is found that influencing the physics that coal body plastic zone width factor includes mining height, coal
Mechanical property, Support Resistance etc. are related, also related with Lode's parameters μ.
Traditional D-P criterion calculation formula (i.e. Lode's parameters μ=0):
After coal body plastic zone width computation model and traditional D-P criterion calculation formula are compared it is found that introducing Lode's parameters μ
Calculated value is not less than traditional calculating formula, in engineering practice can be relatively safe.
Step 2 obtains coal body mining height M under current mining condition, buried depth H, face guard to Support Resistance p, the coal of coal wall
Layer and roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c, bulk density γ;(see by taking infrared ray light shine as an example
Table 1)
1. working face mining condition analysis table of table
Step 3, according to internal friction angleMaterial parameter α, k is calculated with cohesive force c, is shown in Table 2;
Table 2 α, k calculate Data-Statistics
Step 4 takes different Lode's parameters μ, in conjunction with the coal body mining height M under known history mining condition, buried depth H, shield
Side board is to the Support Resistance p of coal wall, coal seam and roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c,
Bulk density γ calculates the coal body plastic zone width (being shown in Table 3) under history mining condition using coal body plastic zone width computation model;
The different Lode's parameters calculated values of table 3 and measured value statistical form
Step 5, the coal body plastic zone width under the calculated history mining condition of more different Lode's parameters μ and history
The difference of practical coal body plastic zone width under mining condition, using Lode's parameters μ corresponding to minimal difference as coal body plasticity
Lode's parameters μ in sector width computation model, (being shown in Table 4);
4 calculated value of table and measured value statistical form
As shown in Table 4, infrared ray light shine should choose Lode's parameters | μ |=0.8, and calculated value is seen slightly larger than mind east scene
Measured value not only improves the accuracy for calculating data, but also relatively safe in engineering practice, therefore compared to traditional calculating formula, exist
Biggish advantage, it follows that Lode's parameters should be selected in analysis infrared ray light shine coal body plastic zone width | μ |=0.8 is counted
It calculates.
Step 6 calculates the coal body plastic zone width that working face is newly exploited in model prediction using coal body plastic zone width.
The present invention mainly on the basis of D-P criterion, considers influence of the intermediate principal stress to yield function, makes up traditional coal body
The defect of Width of Plastic Zone calculation formula, and its engineering practice proves there is good reliability, is the safe and efficient life of working face
Production lays the foundation.
Claims (1)
1. a kind of calculation method of coal body plastic zone width characterized by comprising
Step 1 establishes coal body plastic zone width computation model;
Wherein, coal body mining height M, buried depth H, face guard are to the Support Resistance p of coal wall, coal seam and roof and floor coefficient of friction f, stress collection
Middle COEFFICIENT K is obtained from scene;Internal friction angleCohesive force c, bulk density γ experiment measure;Intermediate parametersμ is Lip river
Moral parameter reacts the stress of coal body, | μ |≤1, material parameter α, k is according to internal friction angleCohesive force c is calculated;
Step 2, obtain coal body mining height M under current mining condition, buried depth H, face guard to the Support Resistance p of coal wall, coal seam with
Roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c, bulk density γ;
Step 3, according to internal friction angleMaterial parameter α, k is calculated with cohesive force c;
Step 4 takes different Lode's parameters μ, in conjunction with coal body mining height M, the buried depth H, face guard under known history mining condition
To the Support Resistance p of coal wall, coal seam and roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c, bulk density
γ calculates the coal body plastic zone width under history mining condition using coal body plastic zone width computation model;
Step 5, the coal body plastic zone width under the calculated history mining condition of more different Lode's parameters μ and history are dug up mine
Under the conditions of practical coal body plastic zone width difference, Lode's parameters μ corresponding to minimal difference is wide as coal body plastic zone
Spend the Lode's parameters μ in computation model;
Step 6 calculates the coal body plastic zone width that working face is newly exploited in model prediction using coal body plastic zone width;
The step 1, comprising:
Step 1-1, force analysis is carried out to cell cube in coal body, establishes the stress balance equation of cell cube;
In formula: σx、σyRespectively indicate first principal stress suffered by cell cube and third principal stress, M-mining height;C-cohesive force;f—
Coal seam and roof and floor coefficient of friction;
Step 1-2, the consideration intermediate principal stress σ parallel with working face inclined direction2Influence to yield function, introduce reflection by
The Lode's parameters μ of power state, | μ |≤1, intermediate principal stress σ is obtained according to plane strain problems2:
Enable σ1=-σx, σ3=-σy, then
Step 1-3, when coal body reaches yield situation, using D-P criterion;
In formula: the first stress invariant I1=σ1+σ2+σ3;
Second deviator of stress invariant
Step 1-4, cell cube first principal stress and third principal stress relational expression are obtained by D-P criterion;
Step 1-5, by relational expressionBring the stress balance equation of cell cube intoAnd combine boundary condition x=0, σx=p, p are Support Resistance of the face guard to coal wall, are obtained vertical
Stress σyExpression formula;
Step 1-6, maximal support pressure expression formula σ=K γ H according to coal body on interface, obtains coal body plastic zone width
Computation model:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611247319.XA CN106643610B (en) | 2016-12-29 | 2016-12-29 | A kind of calculation method of coal body plastic zone width |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611247319.XA CN106643610B (en) | 2016-12-29 | 2016-12-29 | A kind of calculation method of coal body plastic zone width |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106643610A CN106643610A (en) | 2017-05-10 |
CN106643610B true CN106643610B (en) | 2019-05-14 |
Family
ID=58835929
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611247319.XA Active CN106643610B (en) | 2016-12-29 | 2016-12-29 | A kind of calculation method of coal body plastic zone width |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106643610B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107665285B (en) * | 2017-10-31 | 2020-09-08 | 河南理工大学 | Analytical method for determining reasonable position of coal roadway under left coal pillar based on sensitive factors |
CN109441537B (en) * | 2018-10-19 | 2020-12-04 | 中国矿业大学 | Design method for coal mine underground reservoir coal pillar dam body |
CN110096809B (en) * | 2019-04-30 | 2023-03-14 | 中煤科工集团重庆研究院有限公司 | Modeling method for material unstable roadway rock burst based on double-yield contour model |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104594899A (en) * | 2014-11-24 | 2015-05-06 | 山西潞安环保能源开发股份有限公司 | Method for determining reasonable width of gob-side entrydriving narrow coal pillar |
CN105825018A (en) * | 2016-03-18 | 2016-08-03 | 山东科技大学 | Design method for close-distance coal roadway supporting structure and supporting parameters below coal mine goaf |
-
2016
- 2016-12-29 CN CN201611247319.XA patent/CN106643610B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104594899A (en) * | 2014-11-24 | 2015-05-06 | 山西潞安环保能源开发股份有限公司 | Method for determining reasonable width of gob-side entrydriving narrow coal pillar |
CN105825018A (en) * | 2016-03-18 | 2016-08-03 | 山东科技大学 | Design method for close-distance coal roadway supporting structure and supporting parameters below coal mine goaf |
Non-Patent Citations (3)
Title |
---|
《基于Drucker-Prager屈服准则的煤壁内塑性区宽度的确定》;李占魁 等;《煤矿安全》;20130531;第44卷(第5期);第1页第1栏最后一段至第3页第1栏第1段,表1,图1、2 * |
《基于Lode参数的应力状态描述及D-P系列准则》;楚锡华 等;《武汉理工大学学报》;20090831;第31卷(第16期);全文 * |
《考虑峰后应变软化与扩容的圆形巷道围岩弹塑性D-P准则解》;张小波 等;《采矿与安全工程学报》;20131130;第30卷(第6期);第2页第2栏最后一段至第3页第1栏倒数第3段 * |
Also Published As
Publication number | Publication date |
---|---|
CN106643610A (en) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106643610B (en) | A kind of calculation method of coal body plastic zone width | |
WO2020119177A1 (en) | Wall continuous mining and continuous filling water-preserved coal mining method, and water resource migration monitoring and water disaster early warning method | |
WO2016095462A1 (en) | TS-q METHOD FOR EVALUATING FLOOR WATER INRUSH RISK IN DEEP MINE COAL SEAM MINING | |
WO2020048137A1 (en) | Method for determining internal stress of solid-cut-and-fill material | |
CN103590846B (en) | The ratification method of the reasonable working resistance of a kind of working surface hydraulic support | |
RU2671502C2 (en) | Structural discrimination indexes of ordovician limestone top filling zones and determination method | |
CN109117586B (en) | Bedding rock slope three-dimensional geological model building and stability evaluation method | |
CN104111483A (en) | Exploration method of faulted basin slope oil reservoir distribution | |
CN103225509B (en) | Large mining height fully-mechanized mining stope roof plate classification and supporting resistance determination method | |
CN107060761B (en) | A kind of nearly horizontal mining in combined mining thick coal seam stope base object model first breaking space computing method | |
CN106593524A (en) | Solid filling exploitation earth surface subsidence dynamic prediction method | |
CN106815445A (en) | A kind of shallow buried coal seam section coal pillar LOAD FOR model and its computational methods | |
CN106503284B (en) | Shale gas horizontal well horizontal segment gas-bearing formation produces gas evaluation method | |
CN101845815A (en) | Method for calculating and processing uprush plastic failures of confined water foundation pits of soft soil strata | |
CN104750940B (en) | Dynamic strength design method for cemented filling body of underground stope | |
CN105447319A (en) | Determination method for different geological interface calculation parameters of sliding blocks of arch dam abutment | |
CN106121721B (en) | A kind of tight roof Face Ground Pressure Behavior stage division | |
CN109059833B (en) | Artesian water leads a liter band method for determining height under temperature-action of seepage-stress coupling | |
CN111159643B (en) | Method for influencing ground subsidence by Weibull time function based on measured data | |
CN111912953B (en) | Deep-well mining slope stability determination method based on excavation amount monitoring | |
CN107387039A (en) | Utilize the method for the continuous pressure measurement data inversion carbon dioxide flooding leading edge of multiple spot | |
CN106777655B (en) | Method and device for calculating equivalent mining thickness ratio of coal mine goaf of overhead transmission line | |
CN106503854A (en) | Longwall top coal caving coal seam top covering rockmass is across splitting high predicted method | |
Pavlenkova et al. | The upper mantle structure of the Northern Eurasia from the seismic profiling with nuclear explosions | |
Ilyukhin et al. | The use of the finite element method for ensuring efficient and safe extraction of minerals |
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 |