CN113405908A - Shallow shield tunnel roof stability judgment method - Google Patents
Shallow shield tunnel roof stability judgment method Download PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N9/00—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0001—Type of application of the stress
- G01N2203/0003—Steady
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0017—Tensile
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0019—Compressive
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
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
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:
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:
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:
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:
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:
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:
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:
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:
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:
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.
6. The method of claim 5, wherein when B is 0.5, the expression (1) and the expression (3) are squared simultaneously, resulting in:
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:
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.
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)
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 |
-
2021
- 2021-07-26 CN CN202110843766.6A patent/CN113405908B/en active Active
Patent Citations (3)
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)
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 |