CN106643610A - Method for computing widths of plastic regions of coal bodies - Google Patents

Method for computing widths of plastic regions of coal bodies Download PDF

Info

Publication number
CN106643610A
CN106643610A CN201611247319.XA CN201611247319A CN106643610A CN 106643610 A CN106643610 A CN 106643610A CN 201611247319 A CN201611247319 A CN 201611247319A CN 106643610 A CN106643610 A CN 106643610A
Authority
CN
China
Prior art keywords
coal body
coal
stress
lode
parameters
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.)
Granted
Application number
CN201611247319.XA
Other languages
Chinese (zh)
Other versions
CN106643610B (en
Inventor
邸帅
王继仁
张春华
郝朝瑜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liaoning Technical University
Original Assignee
Liaoning Technical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN201611247319.XA priority Critical patent/CN106643610B/en
Publication of CN106643610A publication Critical patent/CN106643610A/en
Application granted granted Critical
Publication of CN106643610B publication Critical patent/CN106643610B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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 invention provides a method for computing the widths of plastic regions of coal bodies. The method includes building coal body plastic region width computing models; acquiring coal body mining heights, buried depths, coal wall support resistance of side wall protection plates, friction coefficients of coal seams, roofs and floors, stress concentration coefficients, internal friction angles, cohesive force and volume weights under current mining conditions; computing the widths of the plastic regions of the coal bodies by the aid of different Lode parameters according to historical mining conditions, comparing computation results of the different Lode parameters to the actual widths of the plastic regions of the coal bodies, analyzing the computation results and the actual widths of the plastic regions of the coal bodies and utilizing the Lode parameters corresponding to the minimum difference values as Lode parameters for the coal body plastic region width computing models; predicting the widths of plastic regions of coal bodies of new mining working faces by the aid of the coal body plastic region width computing models. The method has the advantages that the Lode parameters mu are introduced on the basis of D-P (Drucker-Prager) criteria, so that actual stress states of the coal bodies can be reflected, and the prediction accuracy can be improved.

Description

A kind of computational methods of coal body plastic zone width
Technical field
The invention belongs to coal production technical field, more particularly to a kind of computational methods of coal body plastic zone width.
Background technology
In coal production, coal body plastic zone width (i.e. abutment pressure distribution feature) is to roof collapse, roadway maintenance, coal The aspects such as post width have material impact, and production safe and efficient to safeguard work face is significant.Therefore coal body plastic zone The content that width calculation is always extremely paid close attention to both at home and abroad, successively proposes many computing formula, and A.H Wilson theories, big plate split Gap theory, mole-coulomb criterion, D-P criterions etc..
However, the computing formula that above-mentioned theory and criterion are obtained respectively has pluses and minuses.A.H Wilson theories are according to Britain ground The big probable value that matter condition and exploitation example take, therefore the versatility of its formula is restricted；The theoretical analysis model mistake in big plate crack In simple, many influence factors are not taken into account so that formula distortion；Mole-coulomb criterion has clear physical concept and is suitable for The advantages of facilitating, is used widely, but have ignored the impact of intermediate principal stress, it is impossible to explain well coal body surrender or Destruction；D-P criterions are the further research to mole-coulomb criterion, though in view of the impact of intermediate principal stress, for calculating side Just but have ignored impact of the intermediate principal stress to yield function, therefore its calculated value always big probable value.
From the foregoing, the equal existing defects of the theoretical calculation of said method, it is therefore necessary to propose a kind of more accuracy Computational methods, improve it is existing the drawbacks of, be accurately calculate coal body plastic zone width lay the foundation, so as to ensure working face safety Efficiently production.
The content of the invention
For the problem that prior art is present, the present invention provides a kind of computational methods of coal body plastic zone width.
Technical scheme is as follows：
A kind of computational methods of coal body plastic zone width, including：
Step 1, set up coal body plastic zone width calculation model；
Wherein, coal body mining height M, buried depth H, face guard to the Support Resistance p of rib, coal seam and roof and floor coefficient of friction f, should Power coefficient of concentration K is from scene acquisition；Internal friction angleCohesive force c, unit weight γ experiment is measured；Intermediate parametersμ is Lode's parameters, the stress of reaction coal body, | μ |≤1, material parameter α, k is according to internal friction angleCohesive force c is calculated；
Support Resistance p, coals of the coal body mining height M, buried depth H, face guard under step 2, the current mining condition of acquisition to rib Layer and roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c, unit weight γ；
Step 3, according to internal friction angleMaterial parameter α, k is calculated with cohesive force c；
Step 4, different Lode's parameters μ are taken, with reference to the coal body mining height M under known history mining condition, buried depth H, shield Side board is to the Support Resistance p of rib, coal seam and roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c, Unit weight γ, using coal body plastic zone width calculation model the coal body plastic zone width under history mining condition is calculated；
Coal body plastic zone width and history under the history mining condition that step 5, relatively more different Lode's parameters μ are calculated The difference of the actual coal body plastic zone width under mining condition, using the Lode's parameters μ corresponding to minimal difference as coal body plasticity Lode's parameters μ in sector width computation model；
Step 6, the coal body plastic zone width that working face is newly exploited using the model prediction of coal body plastic zone width calculation.
Described step 1, including：
Step 1-1, force analysis is carried out to cell cube in coal body, set up the stress balance equation of cell cube；
In formula：What document generally believed rib is laterally working face, can regard the scope of freedom, therefore the level suffered by coal body as Stress σxTo be much smaller than vertical stress σy, σx、σyThe first principal stress and third principal stress suffered by cell cube is represented respectively, and M- is adopted It is high；C- cohesive force；F- coal seams and roof and floor coefficient of friction；
Step 1-2, consideration intermediate principal stress σ parallel with working face incline direction2Impact to yield function, introduces anti- The Lode's parameters μ of stress is reflected, | μ |≤1 obtains intermediate principal stress σ according to plane strain problems2：
Make σ1=-σx, σ3=-σy, then
Step 1-3, when coal body reaches yield situation, using D-P criterions；
In formula：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 criterions；
Step 1-5, by relational expressionBring the stress balance equation of cell cube intoAnd with reference to boundary condition x=0, σx=p, p are Support Resistance of the face guard to rib, obtain vertical Stress σyExpression formula；
Step 1-6, the maximal support pressure formula σ according to coal body on interface=K γ H, obtain coal body plastic zone Width calculation model：
Beneficial effect：The computational methods of the coal body plastic zone width that the present invention is provided, on the basis of D-P criterions, introduce Lip river Moral parameter μ to react the real stress of coal body, so as to improve the accuracy of prediction.
Description of the drawings
Fig. 1 is the force analysis schematic diagram of cell cube in the specific embodiment of the invention；
Fig. 2 is the computational methods flow chart of coal body plastic zone width in the specific embodiment of the invention.
Specific embodiment
With reference to specific embodiment, the present invention is described in further detail.
Present embodiment provides a kind of computational methods of coal body plastic zone width as shown in Figure 2, including：
Step 1, set up coal body plastic zone width calculation model；
Wherein, coal body mining height M, buried depth H, face guard to the Support Resistance p of rib, coal seam and roof and floor coefficient of friction f, should Power coefficient of concentration K is from scene acquisition；Internal friction angleCohesive force c, unit weight γ experiment is measured；Intermediate parametersμ is Lode's parameters, the stress of reaction coal body, | μ |≤1, material parameter α, k is according to internal friction angleCohesive force c is calculated；
Described step 1, including：
Step 1-1, hypothesis coal body homogeneous, carry out force analysis, as shown in figure 1, setting up cell cube to cell cube in coal body Stress balance equation；
In formula：What document generally believed rib is laterally working face, can regard the scope of freedom, therefore the level suffered by coal body as Stress σxTo be much smaller than vertical stress σy, σx、σyThe first principal stress and third principal stress suffered by cell cube is represented respectively, and M- is adopted It is high；C- cohesive force；F- coal seams and roof and floor coefficient of friction；
Step 1-2, consideration intermediate principal stress σ parallel with working face incline direction2Impact to yield function, introduces anti- The Lode's parameters μ of stress is reflected, | μ |≤1 obtains intermediate principal stress σ according to plane strain problems2：
Make σ1=-σx, σ3=-σy, then
Step 1-3, when coal body reaches yield situation, using D-P criterions；
In formula：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 criterions；
Step 1-5, by relational expressionBring the stress balance equation of cell cube intoAnd with reference to boundary condition x=0, σx=p, p are Support Resistance of the face guard to rib, obtain vertical Stress σyExpression formula；
Step 1-6, the maximal support pressure formula σ according to coal body on interface=K γ H, obtain coal body plastic zone Width calculation model：
From coal body plastic zone width calculation model, coal body plastic zone width factor is affected to include mining height, the physics of coal Mechanical property, Support Resistance etc. are relevant, also relevant with Lode's parameters μ.
Traditional D-P criterion computing formula (i.e. Lode's parameters μ=0)：
Contrast coal body plastic zone width calculation model and traditional D-P criterions computing formula understand, after introducing Lode's parameters μ Calculated value is not less than traditional calculating formula, can be relatively safe in engineering practice.
Support Resistance p, coals of the coal body mining height M, buried depth H, face guard under step 2, the current mining condition of acquisition to rib Layer and roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c, unit weight γ；(see by taking infrared ray light shine as an example Table 1)
The working face mining condition analysis table of table 1.
Step 3, according to internal friction angleMaterial parameter α, k is calculated with cohesive force c, 2 are shown in Table；
Table 2 α, k calculate Data-Statistics
Step 4, different Lode's parameters μ are taken, with reference to the coal body mining height M under known history mining condition, buried depth H, shield Side board is to the Support Resistance p of rib, coal seam and roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c, Unit weight γ, using coal body plastic zone width calculation model the coal body plastic zone width (being shown in Table 3) under history mining condition is calculated；
The different Lode's parameters calculated values of table 3 and measured value statistical form
Coal body plastic zone width and history under the history mining condition that step 5, relatively more different Lode's parameters μ are calculated The difference of the actual coal body plastic zone width under mining condition, using the Lode's parameters μ corresponding to minimal difference as coal body plasticity Lode's parameters μ in sector width computation model, (being shown in Table 4)；
The calculated value of table 4 and measured value statistical form
As shown in Table 4, infrared ray light shine should choose Lode's parameters | μ |=0.8, and its calculated value is seen slightly larger than god east scene Measured value, both improves the accuracy for calculating data, and the relatively safety in engineering practice, therefore compared to traditional calculating formula, exists Larger advantage, it follows that should be counted from Lode's parameters | μ |=0.8 in analysis infrared ray light shine coal body plastic zone width Calculate.
Step 6, the coal body plastic zone width that working face is newly exploited using the model prediction of coal body plastic zone width calculation.
The present invention is mainly on the basis of D-P criterions, it is considered to impact of the intermediate principal stress to yield function, makes up traditional coal body The defect of Width of Plastic Zone computing formula, and its engineering practice proves there is good reliability, is the safe and efficient life of working face Product lays the foundation.

Claims (2)

1. a kind of computational methods of coal body plastic zone width, it is characterised in that include：
Step 1, set up coal body plastic zone width calculation model；
Wherein, coal body mining height M, buried depth H, face guard are to the Support Resistance p of rib, coal seam and roof and floor coefficient of friction f, stress collection Middle COEFFICIENT K is from scene acquisition；Internal friction angleCohesive force c, unit weight γ experiment is measured；Intermediate parametersμ is Lip river moral Parameter, the stress of reaction coal body, | μ |≤1, material parameter α, k is according to internal friction angleCohesive force c is calculated；
Step 2, the coal body mining height M obtained under current mining condition, buried depth H, face guard to the Support Resistance p of rib, coal seam with Roof and floor coefficient of friction f, stress concentration factor K, internal friction angleCohesive force c, unit weight γ；
Step 3, according to internal friction angleMaterial parameter α, k is calculated with cohesive force c；
Step 4, different Lode's parameters μ are taken, with reference to coal body mining height M, buried depth H, face guard under known history mining condition Support Resistance p, coal seam and roof and floor coefficient of friction f, stress concentration factor K, internal friction angle to ribCohesive force c, unit weight γ, using coal body plastic zone width calculation model the coal body plastic zone width under history mining condition is calculated；
Coal body plastic zone width and history mining under the history mining condition that step 5, relatively more different Lode's parameters μ are calculated Under the conditions of actual coal body plastic zone width difference, using the Lode's parameters μ corresponding to minimal difference as coal body plastic zone width Lode's parameters μ in degree computation model；
Step 6, the coal body plastic zone width that working face is newly exploited using the model prediction of coal body plastic zone width calculation.
2. computational methods of coal body plastic zone width according to claim 1, it is characterised in that described step 1, bag Include：
Step 1-1, force analysis is carried out to cell cube in coal body, set up the stress balance equation of cell cube；
dσ x d x - 2 M ( c + fσ y ) = 0
In formula：σx、σyThe first principal stress and third principal stress suffered by cell cube, M-mining height are represented respectively；C-cohesive force；f— Coal seam and roof and floor coefficient of friction；
Step 1-2, consideration intermediate principal stress σ parallel with working face incline direction2Impact to yield function, introduces reflection and receives The Lode's parameters μ of power state, | μ |≤1 obtains intermediate principal stress σ according to plane strain problems2：
Make σ1=-σx, σ3=-σy, then
σ 2 = 1 2 ( σ 1 + σ 3 ) + μ 2 ( σ 1 - σ 3 )
Step 1-3, when coal body reaches yield situation, using D-P criterions；
αI 1 + J 2 1 / 2 - k = 0
In formula：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 criterions；
σ y = 2 k l - 3 α + l + 3 α l - 3 α σ x , l = 3 + μ 2 / 3
Step 1-5, by relational expressionBring the stress balance equation of cell cube into And with reference to boundary condition x=0, σx=p, p are Support Resistance of the face guard to rib, obtain vertical stress σyExpression formula；
σ y = - c f + ( l + 3 α l - 3 α p + c f + 2 k l - 3 α ) e 2 f ( l + 3 α ) M ( l - 3 α ) x
Step 1-6, the maximal support pressure formula σ according to coal body on interface=K γ H, obtain coal body plastic zone width Computation model：
CN201611247319.XA 2016-12-29 2016-12-29 A kind of calculation method of coal body plastic zone width Active CN106643610B (en)

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 true CN106643610A (en) 2017-05-10
CN106643610B 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)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107665285A (en) * 2017-10-31 2018-02-06 河南理工大学 One kind determines to leave coal road rational position analysis method under coal column based on sensitive factor
CN109441537A (en) * 2018-10-19 2019-03-08 中国矿业大学 A kind of coal mine underground reservoir coal column design of dam body method
CN110096809A (en) * 2019-04-30 2019-08-06 中煤科工集团重庆研究院有限公司 The modeling method pressed based on double construction instability type roadway bumps for surrendering equivalent line model

Citations (2)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
张小波 等: "《考虑峰后应变软化与扩容的圆形巷道围岩弹塑性D-P准则解》", 《采矿与安全工程学报》 *
李占魁 等: "《基于Drucker-Prager屈服准则的煤壁内塑性区宽度的确定》", 《煤矿安全》 *
楚锡华 等: "《基于Lode参数的应力状态描述及D-P系列准则》", 《武汉理工大学学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107665285A (en) * 2017-10-31 2018-02-06 河南理工大学 One kind determines to leave coal road rational position analysis method under coal column based on sensitive factor
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
CN109441537A (en) * 2018-10-19 2019-03-08 中国矿业大学 A kind of coal mine underground reservoir coal column design of dam body method
CN110096809A (en) * 2019-04-30 2019-08-06 中煤科工集团重庆研究院有限公司 The modeling method pressed based on double construction instability type roadway bumps for surrendering equivalent line model
CN110096809B (en) * 2019-04-30 2023-03-14 中煤科工集团重庆研究院有限公司 Modeling method for material unstable roadway rock burst based on double-yield contour model

Also Published As

Publication number Publication date
CN106643610B (en) 2019-05-14

Similar Documents

Publication Publication Date Title
WO2020119177A1 (en) Wall continuous mining and continuous filling water-preserved coal mining method, and water resource migration monitoring and water disaster early warning method
WO2019233114A1 (en) Method for evaluating risk of water damage from strata water inrush of overlying strata in mining area
Wang et al. Evolution mechanism of water-flowing zones and control technology for longwall mining in shallow coal seams beneath gully topography
AU2015345707B2 (en) Method for determining earth surface interpenetrated crack distribution and air leakage characteristics in shallow burial coal mining
CN106643610A (en) Method for computing widths of plastic regions of coal bodies
WO2021109922A1 (en) Method for designing width of isolation coal pillar and gasifier in underground gasification under aquifer facing overlying strata
CN105926569A (en) Method for quantitatively evaluating site stability of old goaf in coalmine based on settlement monitoring data
Kovalevs’Ka et al. Research of stress-strain state of cracked coal-containing massif near-the-working area using finite elements technique
CN103728655A (en) Stope face impact risk pre-mining pre-evaluation method
CN105426619A (en) Analysis method for determining longitudinal bending rigidity of shield tunnel
CN105787185B (en) A kind of best coal column mistake of shallow buried coal seam group working face away from design method
CN115062391B (en) Main-passenger water separation method for measured data of water measuring weir behind dam
CN111475924A (en) Unloading rock slope stability analysis method for rainfall induced deformation
CN107886185A (en) A kind of quiet water yield computational methods of mine based on spatial analysis and system
CN106251053A (en) The suitability assessment method of high voltage power transmission transmission tower is built on earth's surface, a kind of sinking land in coalmining areas
Maity et al. Contributions of burst-sweep cycles to Reynolds shear stress over fluvial obstacle marks generated in a laboratory flume
CN104866664A (en) Method for predicting construction difficulty of water seepage prevention curtain-grouting
Zhang et al. A new formula based on computational fluid dynamics for estimating maximum depth of scour by jets from overflow dams
CN104596735B (en) Optimal arrangement method of bridge group
CN108229050B (en) Method for simply and conveniently calculating soil body reinforcing effect between adjacent tamping points under dynamic compaction action
CN106407569A (en) A calculation method for a surface subsidence value under the condition of thick unconsolidated formation thin bedrock
CN108062447A (en) A kind of method predicted diving precipitation and influenced on building settlement
CN108133108A (en) A kind of method predicted artesian water precipitation and influenced on building settlement
Jin et al. Machine learning-based identification of segment joint failure in underground tunnels
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