CN113405908A - Shallow shield tunnel roof stability judgment method - Google Patents

Shallow shield tunnel roof stability judgment method Download PDF

Info

Publication number
CN113405908A
CN113405908A CN202110843766.6A CN202110843766A CN113405908A CN 113405908 A CN113405908 A CN 113405908A CN 202110843766 A CN202110843766 A CN 202110843766A CN 113405908 A CN113405908 A CN 113405908A
Authority
CN
China
Prior art keywords
expression
tunnel
excavation width
limit value
value range
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
CN202110843766.6A
Other languages
Chinese (zh)
Other versions
CN113405908B (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.)
China Construction Fifth Engineering Bureau Co Ltd
Original Assignee
China Construction Fifth Engineering Bureau Co Ltd
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 China Construction Fifth Engineering Bureau Co Ltd filed Critical China Construction Fifth Engineering Bureau Co Ltd
Priority to CN202110843766.6A priority Critical patent/CN113405908B/en
Publication of CN113405908A publication Critical patent/CN113405908A/en
Application granted granted Critical
Publication of CN113405908B publication Critical patent/CN113405908B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N9/00Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0001Type of application of the stress
    • G01N2203/0003Steady
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0017Tensile
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0019Compressive
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Biochemistry (AREA)
  • Evolutionary Computation (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Civil Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Structural Engineering (AREA)
  • Architecture (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)

Abstract

The invention provides a method for judging stability of a top plate of a shallow shield tunnel, which comprises the steps of determining parameters according to preset conditions and setting a preset value of half-excavation width of the tunnel; establishing an expression of a molar strength envelope curve based on parameters, establishing a shearing strength expression of a generalized HockBrownian rule according to a limit analysis principle, further establishing an expression of the tunnel semi-excavation width, calculating a limit value range of the tunnel semi-excavation width, and determining the limit value range of the tunnel semi-excavation width by combining each expression and each parameter; and finally, comparing the predicted value with the limit value range, and if the predicted value is within the limit value range, determining that the tunnel top plate is stable. The potential collapse criterion formula of the shallow tunnel shield is solved according to rigorous mathematical formula derivation, accurate reference is provided for preventing tunnel collapse, a simple, convenient and practical method for determining the limit width of the top plate of the shallow tunnel is provided, and related parameters are easy to determine.

Description

Shallow shield tunnel roof stability judgment method
Technical Field
The invention relates to the technical field of disaster prevention and control of geotechnical and urban underground engineering, in particular to a method for judging stability of a shallow shield tunnel roof.
Background
In the past decades, China has a large number of subways and tunnel constructions, wherein the construction of the subway projects can not avoid the shallow address section. However, accidents such as tunnel collapse are caused, and the tunnel collapse is mostly generated before the tunnel is built, so that stability analysis is necessary before excavation is performed to reduce the risk of safety accidents, prevent the accidents in the future, and have important significance for ensuring the sustainable development of economy and society in China.
The most critical issue for tunnel roof stability analysis is determining its excavation width. The existing patents or other achievements mostly adopt the HockBrownian criterion to analyze the tunnel instability mechanism, but the HockBrownian parameter is difficult to determine in the actual engineering, and the problem still cannot be solved well at present.
Disclosure of Invention
In view of the above, it is necessary to provide a method for determining stability of a top plate of a shallow shield tunnel.
A shallow shield tunnel roof stability judgment method comprises the following steps: according to preset conditions, specific parameters of shield surrounding rocks are determined, and the value of the half-excavation width of the tunnel is preset as a preset value, wherein the parameters comprise uniaxial compressive strength of a rock mass, uniaxial tensile strength of the rock mass and volume weight of the rock mass; establishing an expression of a molar strength envelope curve according to the parameters; establishing a shearing strength expression of a generalized HockBrown criterion according to a limit analysis principle, and further establishing an expression of tunnel half-excavation width; determining a limit value range of the tunnel half-excavation width according to the expression of the tunnel half-excavation width based on the expression of the molar strength envelope curve and the shear strength expression of the generalized HockBrown criterion; and comparing the preset value of the tunnel half-excavation width with the limit value range of the tunnel half-excavation width, and determining that the tunnel top plate is stable when the preset value belongs to the limit value range.
In one embodiment, the expression of the molar intensity envelope is specifically:
τ2=n(σ+σt); (1)
wherein, tau represents the shear strength shear stress of rock-soil body, sigma is positive stress, sigmatThe uniaxial tensile strength of the rock mass is shown, n is a coefficient, and the expression of n is as follows:
Figure BDA0003179779420000021
wherein σcThe uniaxial compressive strength of the rock mass.
In one embodiment, the shear strength expression of the generalized hockbron criterion is specifically:
τ=±Aσc[(σt+σ)σc -1]B; (3)
wherein A and B are HockBrown parameters respectively.
In one embodiment, the expression of the tunnel half-excavation width specifically includes:
Lu=AB-B(1+B)B(γ)-1σc 1-Bσt B; (4)
wherein L isuThe critical width of the tunnel half excavation, namely the limit value range, and gamma is the volume weight of the rock mass.
In one embodiment, according to the expression (1) and the expression (3), B may be equal to 0.5, so the expression (4) is:
Figure BDA0003179779420000022
where γ is ρ g, ρ is density, g is gravitational acceleration.
In one embodiment, when B is 0.5, the expression (1) and the expression (3) are squared simultaneously, and then:
Figure BDA0003179779420000023
on the basis of the expression (6), and in combination with the expression (5), the range of the limit value of the tunnel half excavation is obtained as follows:
Figure BDA0003179779420000024
where γ is ρ g, ρ is density, g is gravitational acceleration.
In one embodiment, the establishing of the shear strength expression of the generalized hockbron criterion according to the limit analysis principle further establishes an expression of the tunnel half-excavation width, specifically: establishing a shearing strength expression of a generalized HockBrown criterion according to a limit analysis principle; and establishing an expression of the tunnel half-excavation width according to the shearing strength expression of the generalized HockBrownian criterion based on the limit analysis principle.
In one embodiment, the comparing the preset value of the tunnel half-cut width with the limit value range of the tunnel half-cut width, and when the preset value falls within the limit value range, after the tunnel roof is determined to be stable, the method further includes: and when the preset value exceeds the limit value range, resetting the preset value of the tunnel half-excavation width, so that the preset value belongs to the limit value range.
The invention has the beneficial effects that: according to the method for judging the stability of the top plate of the shallow shield tunnel, parameters are determined according to preset conditions, and the preset value of the half-excavation width of the tunnel is set; establishing an expression of a molar strength envelope curve based on parameters, establishing a shearing strength expression of a generalized HockBrownian rule according to a limit analysis principle, further establishing an expression of the tunnel semi-excavation width, calculating a limit value range of the tunnel semi-excavation width, and determining the limit value range of the tunnel semi-excavation width by combining each expression and each parameter; and finally, comparing the predicted value with the limit value range, and if the predicted value is within the limit value range, determining that the tunnel top plate is stable. The potential collapse criterion formula of the shallow tunnel shield is solved according to rigorous mathematical formula derivation, accurate reference is provided for preventing tunnel collapse, a simple, convenient and practical method for determining the limit width of the top plate of the shallow tunnel is provided, and related parameters are easy to determine.
Drawings
Fig. 1 is a scene model diagram of a shallow tunnel shield roof stability determination method in an embodiment;
fig. 2 is a schematic flow chart of a method for determining stability of a top plate of a shallow shield tunnel according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings by way of specific embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The method for judging the stability of the top plate of the shallow shield tunnel can be applied to a scene model diagram shown in figure 1. Wherein, 1-a collapse zone; 2-a rectangular tunnel; 3-circular tunnel; l isu-the critical width of the tunnel half excavation, i.e. the limit value range of the tunnel half excavation; the actual half-excavation width of the L-tunnel is a preset value in the following; h is the collapse zone height; xoy is as shown to establish a coordinate system(ii) a (x) is the tunnel roof collapse curve; h is the tunnel burial depth, (here, the rectangular tunnel burial depth).
In one embodiment, as shown in fig. 2, a method for determining stability of a shallow shield tunnel roof is provided, which includes the following steps:
s110, specific parameters of shield surrounding rocks are determined according to preset conditions, and the value of the half-excavation width of the tunnel is preset as a preset value, wherein the parameters comprise uniaxial compressive strength of a rock mass, uniaxial tensile strength of the rock mass and volume weight of the rock mass.
Specifically, the preset conditions are that specific parameters of the shield surrounding rock, including uniaxial compressive strength sigma of the rock mass, are determined according to construction drawings and other datacUniaxial tensile strength of rock masstThe volume weight γ of the rock mass.
S120, establishing an expression of the molar strength envelope curve according to the parameters.
In one embodiment, the expression of the molar intensity envelope in step S120 is specifically:
τ2=n(σ+σt); (1)
wherein, tau represents the shear strength shear stress of rock-soil body, sigma is positive stress, sigmatThe uniaxial tensile strength of the rock mass is shown, n is a coefficient, and the expression of n is as follows:
Figure BDA0003179779420000041
wherein σcThe uniaxial compressive strength of the rock mass.
S130, according to the limit analysis principle, a shearing strength expression of the generalized HockBrownian rule is established, and an expression of the tunnel half-excavation width is further established.
In one embodiment, step S130 specifically includes: establishing a shearing strength expression of a generalized HockBrown criterion according to a limit analysis principle; based on the principle of limit analysis, an expression of the tunnel half-excavation width is established according to a shearing strength expression of the generalized HockBrown criterion.
In one embodiment, the shearing strength expression of the generalized hockbron criterion in step S130 is specifically:
τ=±Aσc[(σt+σ)σc -1]B; (3)
wherein A and B are HockBrown parameters respectively.
In an embodiment, the expression of the tunnel half-excavation width in step S130 is specifically:
Lu=AB-B(1+B)B(γ)-1σc 1-Bσt B; (4)
wherein L isuThe critical width of the tunnel half excavation, namely the limit value range, and gamma is the volume weight of the rock mass.
Specifically, a shearing strength expression of the generalized HockBrown criterion is established according to the limit analysis principle, and then checking calculation is carried out according to the shearing strength expression of the generalized HockBrown criterion, so that an expression of the tunnel half-excavation width can be obtained.
S140, determining the limit value range of the tunnel half-excavation width according to the expression of the tunnel half-excavation width based on the expression of the molar strength envelope curve and the shearing strength expression of the generalized HockBrown criterion.
In one embodiment, according to expression (1) and expression (3), B may be equal to 0.5, so expression (4) is:
Figure BDA0003179779420000051
where γ is ρ g, ρ is density, g is gravitational acceleration.
In one embodiment, when B is 0.5, expression (1) and expression (3) are squared simultaneously, which may result in:
Figure BDA0003179779420000052
on the basis of the expression (6) and in combination with the expression (5), the range of the limit value of the tunnel half excavation can be obtained as follows:
Figure BDA0003179779420000053
where γ is ρ g, ρ is density, g is gravitational acceleration. Specifically, based on the expression of the molar strength envelope and the shear strength expression of the generalized hockbron criterion, the molar strength envelope is obtained first, and the actual B is 0.5, so that substituting B into 0.5 can simplify the expression of the tunnel half-excavation width. And then, squaring two sides according to an expression of a molar strength envelope curve and a shearing strength expression of a generalized HockBrown criterion, combining the squaring two sides, and finally obtaining an expression of the tunnel half-excavation width without calculating HockBrown parameters by combining the simplified expression of the tunnel half-excavation width, and finally calculating the limit value range of the tunnel half-excavation width according to the parameters.
S150, comparing the preset value of the tunnel half-excavation width with the limit value range of the tunnel half-excavation width, and determining that the tunnel roof is stable when the preset value belongs to the limit value range.
Specifically, if the preset value falls within the limit value range, the tunnel corresponding to the preset value is stable and will not collapse.
In one embodiment, after step S150, the method further includes: and when the preset value exceeds the limit value range, resetting the preset value of the tunnel half excavation width, so that the preset value is within the limit value range. Specifically, if the preset value exceeds the limit value range, the preset value needs to be reset, and the limit value range is stable and cannot collapse for the tunnel corresponding to the half excavation width of the tunnel in the range, and if the preset value exceeds the limit value range, the corresponding tunnel is unstable and is likely to collapse.
In one embodiment, as shown in FIG. 1, the uniaxial compressive strength and tensile strength of the rock mass are respectively 0.5MPa, 0.06MPa and 26kN/m3. Then, according to the expression (7), the tunnel half-excavation limit width can be determined to be 6.71m, and the limit value range is<6.71 m. Based on the same parameter setting, calculating the obtained half-excavation width limit value by adopting an ANSYS numerical model, wherein the maximum value area of the plastic zone is 6.59 m; using FLAC numerical modulusThe maximum plastic zone area of the half-excavation width limit value obtained by the model calculation is 6.63 m. It can be seen that the calculation result of the method provided by the invention is close to the result obtained by the reliable method of the prior mature technology, which shows that the method of the invention is accurate and feasible.
In the embodiment, parameters are determined according to preset conditions, and the preset value of the tunnel half-excavation width is set; establishing an expression of a molar strength envelope curve based on parameters, establishing a shearing strength expression of a generalized HockBrownian rule according to a limit analysis principle, further establishing an expression of the tunnel semi-excavation width, calculating a limit value range of the tunnel semi-excavation width, and determining the limit value range of the tunnel semi-excavation width by combining each expression and each parameter; and finally, comparing the predicted value with the limit value range, and if the predicted value is within the limit value range, determining that the tunnel top plate is stable. The potential collapse criterion formula of the shallow tunnel shield is solved according to rigorous mathematical formula derivation, accurate reference is provided for preventing tunnel collapse, a simple, convenient and practical method for determining the limit width of the top plate of the shallow tunnel is provided, and related parameters are easy to determine.
It will be apparent to those skilled in the art that the modules or steps of the invention described above may be implemented in a general purpose computing device, they may be centralized on a single computing device or distributed across a network of computing devices, and optionally they may be implemented in program code executable by a computing device, such that they may be stored on a computer storage medium (ROM/RAM, magnetic disks, optical disks) and executed by a computing device, and in some cases, the steps shown or described may be performed in an order different than that described herein, or they may be separately fabricated into individual integrated circuit modules, or multiple ones of them may be fabricated into a single integrated circuit module. Thus, the present invention is not limited to any specific combination of hardware and software.
The foregoing is a more detailed description of the present invention that is presented in conjunction with specific embodiments, and the practice of the invention is not to be considered limited to those descriptions. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.

Claims (8)

1. A shallow shield tunnel roof stability judgment method is characterized by comprising the following steps:
according to preset conditions, specific parameters of shield surrounding rocks are determined, and the value of the half-excavation width of the tunnel is preset as a preset value, wherein the parameters comprise uniaxial compressive strength of a rock mass, uniaxial tensile strength of the rock mass and volume weight of the rock mass;
establishing an expression of a molar strength envelope curve according to the parameters;
establishing a shearing strength expression of a generalized HockBrown criterion according to a limit analysis principle, and further establishing an expression of tunnel half-excavation width;
determining a limit value range of the tunnel half-excavation width according to the expression of the tunnel half-excavation width based on the expression of the molar strength envelope curve and the shear strength expression of the generalized HockBrown criterion;
and comparing the preset value of the tunnel half-excavation width with the limit value range of the tunnel half-excavation width, and determining that the tunnel top plate is stable when the preset value belongs to the limit value range.
2. The method of claim 1, wherein the molar intensity envelope is expressed by:
τ2=n(σ+σt); (1)
wherein, tau represents the shear strength shear stress of rock-soil body, sigma is positive stress, sigmatThe uniaxial tensile strength of the rock mass is shown, n is a coefficient, and the expression of n is as follows:
Figure FDA0003179779410000011
wherein σcThe uniaxial compressive strength of the rock mass.
3. The method according to claim 2, characterized in that the shear strength expression of the generalized HockBrownian criterion is in particular:
τ=±Aσc[(σt+σ)σc -1]B; (3)
wherein A and B are HockBrown parameters respectively.
4. The method according to claim 3, wherein the expression of the tunnel half-excavation width is:
Lu=AB-B(1+B)B(γ)-1σc 1-Bσt B; (4)
wherein L isuThe critical width of the tunnel half excavation, namely the limit value range, and gamma is the volume weight of the rock mass.
5. The method according to claim 4, wherein B is 0.5 according to the expression (1) and the expression (3), so that the expression (4) is:
Figure FDA0003179779410000021
where γ is ρ g, ρ is density, g is gravitational acceleration.
6. The method of claim 5, wherein when B is 0.5, the expression (1) and the expression (3) are squared simultaneously, resulting in:
Figure FDA0003179779410000022
on the basis of the expression (6), and in combination with the expression (5), the range of the limit value of the tunnel half excavation is obtained as follows:
Figure FDA0003179779410000023
where γ is ρ g, ρ is density, g is gravitational acceleration.
7. The method according to claim 1, wherein the shear strength expression of the generalized HockBrownian criterion is established according to the principle of limit analysis, and the expression of the tunnel half-excavation width is further established, specifically:
establishing a shearing strength expression of a generalized HockBrown criterion according to a limit analysis principle;
and establishing an expression of the tunnel half-excavation width according to the shearing strength expression of the generalized HockBrownian criterion based on the limit analysis principle.
8. The method of claim 1, wherein the comparing the preset value of the tunnel half-excavation width with a limit value range of the tunnel half-excavation width, and when the preset value falls within the limit value range, after the tunnel roof is determined to be stable, further comprises:
and when the preset value exceeds the limit value range, resetting the preset value of the tunnel half-excavation width, so that the preset value belongs to the limit value range.
CN202110843766.6A 2021-07-26 2021-07-26 Shallow shield tunnel roof stability judgment method Active CN113405908B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110843766.6A CN113405908B (en) 2021-07-26 2021-07-26 Shallow shield tunnel roof stability judgment method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110843766.6A CN113405908B (en) 2021-07-26 2021-07-26 Shallow shield tunnel roof stability judgment method

Publications (2)

Publication Number Publication Date
CN113405908A true CN113405908A (en) 2021-09-17
CN113405908B CN113405908B (en) 2022-10-21

Family

ID=77687576

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110843766.6A Active CN113405908B (en) 2021-07-26 2021-07-26 Shallow shield tunnel roof stability judgment method

Country Status (1)

Country Link
CN (1) CN113405908B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108681523A (en) * 2018-05-06 2018-10-19 北京工业大学 Small-clear-distance tunnel depth, shallow embedding critical depth computational methods
CN109026056A (en) * 2018-08-13 2018-12-18 长安大学 The horizontal layered rock tunnel critical span calculation method of anchoring beam structure
CN110222405A (en) * 2019-05-30 2019-09-10 长安大学 A kind of Stability of Slope In Rock test method based on reduction strength envelope reduction method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108681523A (en) * 2018-05-06 2018-10-19 北京工业大学 Small-clear-distance tunnel depth, shallow embedding critical depth computational methods
CN109026056A (en) * 2018-08-13 2018-12-18 长安大学 The horizontal layered rock tunnel critical span calculation method of anchoring beam structure
CN110222405A (en) * 2019-05-30 2019-09-10 长安大学 A kind of Stability of Slope In Rock test method based on reduction strength envelope reduction method

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
ALI GHORBANI ET AL.: "A comprehensive solution for the calculation of ground reaction curve in the crown and sidewalls of circular tunnels in the elastic-plastic-EDZ rock mass considering strain softening", 《TUNNELLING AND UNDERGROUND SPACE TECHNOLOGY》 *
G.D.MANEV ET AL.: "ESTIMATION OF THE ENVELOPE OF MAXIMUM MOHRS CIRCLES PREPARING THE STRENGTH TEST CERTIFICATE OF ROCKS", 《JOURNAL OF MINING SCIENCE》 *
LIAN-HENG ZHAO ET AL.: "Variational analysis of the ultimate pullout capacity of shallow circular anchor plates in rock foundations based on the Hoek-Brown failure criterion", 《INTERNATIONAL JOURNAL OF ROCK MECHANICS AND MINING SCIENCES》 *
LINJIAN MA ET AL.: "Applicability of a new modified explicit three-dimensional Hoek-Brown failure criterion to eight rocks", 《INTERNATIONAL JOURNAL OF ROCK MECHANICS&MINING SCIENCES》 *
M. FRALDI ET AL.: "Limit analysis of collapse mechanisms in cavities and tunnels according to the Hoek–Brown failure criterion", 《INTERNATIONAL JOURNAL OF ROCK MECHANICS & MINING SCIENCES》 *
X.L. YANG ET AL.: "Collapse mechanism of shallow tunnel based on nonlinear Hoek–Brown failure criterion", 《TUNNELLING AND UNDERGROUND SPACE TECHNOLOGY》 *
汪斌 等: "高应力下岩石非线性强度特性的试验验证", 《岩石力学与工程学报》 *

Also Published As

Publication number Publication date
CN113405908B (en) 2022-10-21

Similar Documents

Publication Publication Date Title
CN106897510A (en) A kind of bridge structure 3-D seismics vulnerability analysis method
Ma et al. Mechanical behavior of a coal seam penetrated by a karst collapse pillar: mining-induced groundwater inrush risk
Krejsa et al. The use of the direct optimized probabilistic calculation method in design of bolt reinforcement for underground and mining workings
CN106593524A (en) Solid filling exploitation earth surface subsidence dynamic prediction method
CN103577642A (en) Method for determining safe distance between concealed karst cave and newly built mountain tunnel
CN107305700A (en) Judge coal and gas whether the method that can be protruded
CN109446659A (en) A kind of method that integrated forecasting heading emits heights of roofs and hidden danger classification
CN112488455A (en) Slope stability analysis and evaluation method, system, medium, equipment and terminal
Stark et al. Selecting minimum factors of safety for 3D slope stability analyses
CN113405908B (en) Shallow shield tunnel roof stability judgment method
Villaescusa et al. Probabilistic estimate of rock mass static and dynamic demands for underground excavation stabilisation
Deng et al. Numerical simulation on the evolution of mining-induced fracture network in a coal seam and its overburden under the top coal caving method
CN107315868A (en) A kind of design method and system of the base pit engineering stability based on Revit secondary development
Lazar et al. A closed form solution for seismic risk assessment incorporating intensity bounds
CN105970975B (en) A kind of base pit stand construction with basement building
CN115618614B (en) Stope old roof triangle area movement judging method, system, equipment and storage medium
CN113435087B (en) Method for analyzing local stability of cave surrounding rock
CN113094905B (en) Calculation method and system suitable for multi-middle-section continuous empty area support key points
Wei et al. Numerical study of damage to rock surrounding an underground coal roadway excavation
Zhu et al. A Risk Assessment Model for Dam Combining the Probabilistic and the Nonprobabilistic Methods
Wang et al. Analysis of surface deformation and settlement characteristics caused by tunnel excavation and unloading
CN110130984B (en) Mine water inrush flooding estimation method and computing device
CN112699559A (en) Safety evaluation method and system based on settling rate information entropy safety value
Nan et al. A predictive model of mining collapse extent and its application
CA2361460A1 (en) Method for performing wedge analysis for assessing wedge instabilities in underground openings

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