CN114969887B - 一种确定含弱层边坡稳定性随内排支挡变化规律的方法 - Google Patents
一种确定含弱层边坡稳定性随内排支挡变化规律的方法 Download PDFInfo
- Publication number
- CN114969887B CN114969887B CN202210393657.3A CN202210393657A CN114969887B CN 114969887 B CN114969887 B CN 114969887B CN 202210393657 A CN202210393657 A CN 202210393657A CN 114969887 B CN114969887 B CN 114969887B
- Authority
- CN
- China
- Prior art keywords
- inner row
- under
- supporting
- condition
- stability
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000011435 rock Substances 0.000 claims abstract description 31
- 230000000903 blocking effect Effects 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 18
- 238000005452 bending Methods 0.000 claims description 3
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 239000003245 coal Substances 0.000 abstract description 5
- 238000011084 recovery Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/23—Dune restoration or creation; Cliff stabilisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Geometry (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Structural Engineering (AREA)
- Operations Research (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Civil Engineering (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Probability & Statistics with Applications (AREA)
- Algebra (AREA)
- Architecture (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种确定含弱层边坡稳定性随内排支挡变化规律的方法，方法首先在未进行内排支挡的情况下和全部内排支挡的情况下分别采用二维极限平衡法计算含弱层边坡的稳定性系数；然后构建边坡最下岩体的简支梁模型，计算在未进行内排支挡的情况下和在进行部分内排支挡的情况下的简支梁理论定义下的含弱层边坡稳定性系数；接着计算两稳定性系数的比值，得到倍数m与最下岩体部分内排支挡长度L2之间的关系；最后确定含弱层边坡稳定性随内排支挡长度的变化规律，并绘制变化规律曲线。本发明的方法将空间力学问题转化为平面应变问题，简化了求解过程，并确定了含弱层边坡的稳定性随内排支挡长度的变化规律，为露天煤矿资源回采提供理论基础。
Description
技术领域
本发明涉及露天开采技术领域，尤其涉及一种确定含弱层边坡稳定性随内排支挡变化规律的方法。
背景技术
在某些地区露天煤矿赋存弱层，极易发生以圆弧为侧界面，以弱层为底界面的切层-顺层组合滑坡，治理该滑坡的最佳措施为内排压脚。当排弃物料无法实现全部压帮时，边坡具有显著的三维空间效应。而传统的二维极限平衡法，无法考虑边坡的三维空间效应；现有的三维极限平衡法也无法适用于该滑坡模式的分析；此外，由于数值模拟方法建模的复杂性较高，也难以在工程中推广。因此，迫切需要提出一种新的方法用来确定含弱层边坡稳定性随内排支挡的变化规律，为露天煤矿资源回采提供理论基础。
发明内容
针对上述现有技术的不足，本发明提供一种确定含弱层边坡稳定性随内排支挡变化规律的方法。
为解决上述技术问题，本发明所采取的技术方案是：一种确定含弱层边坡稳定性随内排支挡变化规律的方法，具体包括如下步骤：
步骤1：在未进行内排支挡的情况下，采用二维极限平衡法计算含弱层边坡的稳定性系数，记为Fs；
步骤2：在进行全部内排支挡的情况下，采用二维极限平衡法计算含弱层边坡的稳定性系数，记为Fs1；
步骤3：以含弱层边坡最下岩体为研究对象，构建该岩体的简支梁模型，并分别计算在未进行内排支挡的情况下和在进行部分内排支挡的情况下的简支梁理论定义下的含弱层边坡稳定性系数，分别记为Fs3和Fs4，具体过程如下：
步骤3.1：将边坡最下岩体剩余推力与最下岩体底界面抗剪力的合力记为P，边坡最下岩体走向长度记为L1，最下岩体部分内排支挡长度记为L2；
步骤3.2：在未进行内排支挡的情况下，计算简支梁理论定义下的含弱层边坡稳定性系数Fs3：
其中，σt为抗拉强度，wt为允许变形值；σ1max为未进行内排支挡的情况下简支梁的最大拉应力；w1max为未进行内排支挡的情况下简支梁的最大挠度；
步骤3.3：在进行部分内排支挡的情况下，计算简支梁理论定义下的含弱层边坡稳定性系数Fs4：
其中，σ2max为进行部分内排支挡的情况下简支梁的最大拉应力；w2max为进行部分内排支挡的情况下简支梁的最大挠度。
进一步的，所述未进行内排支挡的情况下简支梁的最大拉应力的计算公式为：
所述未进行内排支挡的情况下简支梁的最大挠度的计算公式为：
所述进行部分内排支挡的情况下简支梁的最大拉应力的计算公式为：
所述进行部分内排支挡的情况下简支梁的最大挠度的计算公式为：
其中，H为简支梁的高度，b为底宽，EI为简支梁的弯曲刚度。
步骤4：计算Fs4和Fs3的比值，得到倍数m与最下岩体部分内排支挡长度L2之间的关系，具体如下：
步骤5：比较m倍的Fs与Fs1的大小关系，确定含弱层边坡稳定性随内排支挡变化规律，具体如下；
其中，ΔF为内排支挡时边坡的稳定性系数与未进行支挡时边坡稳定性的差值。
步骤6：绘制含弱层边坡稳定性随内排支挡长度的变化规律曲线。
采用上述技术方案所产生的有益效果在于：本发明提供的方法充分利用二维极限平衡法操作简单、计算快捷的优势，从简支梁变形破坏与拉裂破坏角度出发，将空间力学问题转化为平面应变问题求解，克服了三维稳定性计算方法无法适用、数值模拟方法建模复杂等不足。本发明的方法还可用于露天煤矿回采端帮压覆资源时边坡形态参数优化。
附图说明
图1为本发明实施例中一种确定含弱层边坡稳定性随内排支挡变化规律的方法流程图；
图2为本发明实施例中未进行内排支挡时含弱层边坡稳定性系数计算结果示意图；
图3为本发明实施例中进行全部内排支挡时含弱层边坡稳定性系数计算结果示意图；
图4为本发明实施例中部分内排支挡的情况下含弱层边坡最下岩体的受力情况分析图；
图5为本发明实施例中未进行内排支挡时含弱层边坡最下岩体的简支梁模型示意图；
图6为本发明实施例中进行内排支挡时含弱层边坡最下岩体的简支梁模型示意图；
图7为本发明实施例中含弱层边坡最下岩体的截面示意图；
图8为本发明实施例中含弱层边坡稳定性随内排支挡变化规律曲线图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
本实施例中某露天煤矿东帮为含弱层边坡，边坡走向长度L1为1000m。
如图1所示，本实施例中一种确定含弱层边坡稳定性随内排支挡变化规律的方法如下所述。
步骤1：在未进行内排支挡的情况下，采用二维极限平衡法计算含弱层边坡的稳定性系数，记为Fs；
本实施例中，基于二维刚体极限平衡法计算未进行内排支挡时含弱层边坡稳定性系数Fs＝1.2，如图2所示。
步骤2：在进行全部内排支挡的情况下，采用二维极限平衡法计算含弱层边坡的稳定性系数，记为Fs1；
本实施例中，进行全部内排支挡时，采用二维极限平衡法计算含弱层边坡稳定性系数Fs1＝1.70，如图3所示。
步骤3：以含弱层边坡最下岩体为研究对象，构建该岩体的简支梁模型，并分别计算在未进行内排支挡的情况下和在进行部分内排支挡的情况下的简支梁理论定义下的含弱层边坡稳定性系数，分别记为Fs3和Fs4，具体过程如下：
步骤3.1：将边坡最下岩体剩余推力与最下岩体底界面抗剪力的合力记为P，边坡最下岩体走向长度记为L1，最下岩体部分内排支挡长度记为L2；
本实施例中，边坡最下岩体构建为简支梁模型后的受力情况如图4所示。
步骤3.2：在未进行内排支挡的情况下，计算简支梁理论定义下的含弱层边坡稳定性系数Fs3：
其中，σt为抗拉强度，wt为允许变形值；σ1max为未进行内排支挡的情况下简支梁的最大拉应力；w1max为未进行内排支挡的情况下简支梁的最大挠度；
本实施例中，在未进行内排支挡的情况下，含弱层边坡最下岩体构建的简支梁结构如图5所示，其中A、C分别为边坡最下岩体走向长度的起始点和结束点，FAy为A点的支撑反力，FCy为C点的支撑反力。
步骤3.3：在进行部分内排支挡的情况下，计算简支梁理论定义下的含弱层边坡稳定性系数Fs4：
其中，σ2max为进行部分内排支挡的情况下简支梁的最大拉应力；w2max为进行部分内排支挡的情况下简支梁的最大挠度。
本实施例中，在进行部分内排支挡的情况下，含弱层边坡最下岩体构建的简支梁结构如图6所示。
进一步的，所述未进行内排支挡的情况下简支梁的最大拉应力的计算公式为：
所述未进行内排支挡的情况下简支梁的最大挠度的计算公式为：
所述进行部分内排支挡的情况下简支梁的最大拉应力的计算公式为：
所述进行部分内排支挡的情况下简支梁的最大挠度的计算公式为：
其中，H为简支梁的高度，b为底宽，EI为简支梁的弯曲刚度，本实施例中简支梁的截面如图7所示。
步骤4：计算Fs4和Fs3的比值，得到倍数m与最下岩体部分内排支挡长度L2之间的关系，具体如下：
步骤5：比较m倍的Fs与Fs1的大小关系，确定含弱层边坡稳定性随内排支挡变化规律的方法，具体如下；
其中，ΔF为内排支挡时边坡的稳定性系数与未进行支挡时边坡稳定性的差值。
本实施例中，代入已知量计算得到：
步骤6：绘制含弱层边坡稳定性随内排支挡长度的变化规律曲线。
本实施例中，绘制含弱层边坡稳定性随内排支挡长度的变化规律曲线如图8所示，当L2＝400m时，mFs的值等于Fs1。
Claims (1)
1.一种确定含弱层边坡稳定性随内排支挡变化规律的方法，其特征在于，包括如下步骤：
步骤1：在未进行内排支挡的情况下，采用二维极限平衡法计算含弱层边坡的稳定性系数，记为Fs；
步骤2：在进行全部内排支挡的情况下，采用二维极限平衡法计算含弱层边坡的稳定性系数，记为Fs1；
步骤3：以含弱层边坡最下岩体为研究对象，构建该岩体的简支梁模型，并分别计算在未进行内排支挡的情况下和在进行部分内排支挡的情况下的简支梁理论定义下的含弱层边坡稳定性系数，分别记为Fs3和Fs4；
步骤4：计算Fs4和Fs3的比值，得到倍数m与最下岩体部分内排支挡长度L2之间的关系；
步骤5：比较m倍的Fs与Fs1的大小关系，确定含弱层边坡稳定性随内排支挡变化规律；
步骤6：绘制含弱层边坡稳定性随内排支挡长度的变化规律曲线；
所述步骤3的过程如下：
步骤3.1：将边坡最下岩体剩余推力与最下岩体底界面抗剪力的合力记为P，边坡最下岩体走向长度记为L1，最下岩体部分内排支挡长度记为L2；
步骤3.2：在未进行内排支挡的情况下，计算简支梁理论定义下的含弱层边坡稳定性系数Fs3：
其中，σt为抗拉强度，wt为允许变形值；σ1max为未进行内排支挡的情况下简支梁的最大拉应力；w1max为未进行内排支挡的情况下简支梁的最大挠度；
步骤3.3：在进行部分内排支挡的情况下，计算简支梁理论定义下的含弱层边坡稳定性系数Fs4：
其中，σ2max为进行部分内排支挡的情况下简支梁的最大拉应力；w2max为进行部分内排支挡的情况下简支梁的最大挠度；
所述未进行内排支挡的情况下简支梁的最大拉应力的计算公式为：
所述未进行内排支挡的情况下简支梁的最大挠度的计算公式为：所述进行部分内排支挡的情况下简支梁的最大拉应力的计算公式为：
所述进行部分内排支挡的情况下简支梁的最大挠度的计算公式为：其中，H为简支梁的高度，b为底宽，EI为简支梁的弯曲刚度；
所述步骤4的计算过程如下：
所述步骤5的确定过程如下：
其中，ΔF为内排支挡时边坡的稳定性系数与未进行支挡时边坡稳定性的差值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210393657.3A CN114969887B (zh) | 2022-04-15 | 2022-04-15 | 一种确定含弱层边坡稳定性随内排支挡变化规律的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210393657.3A CN114969887B (zh) | 2022-04-15 | 2022-04-15 | 一种确定含弱层边坡稳定性随内排支挡变化规律的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114969887A CN114969887A (zh) | 2022-08-30 |
| CN114969887B true CN114969887B (zh) | 2023-12-26 |
Family
ID=82977629
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210393657.3A Active CN114969887B (zh) | 2022-04-15 | 2022-04-15 | 一种确定含弱层边坡稳定性随内排支挡变化规律的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114969887B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115944276B (zh) * | 2023-02-22 | 2024-07-05 | 武汉大学人民医院（湖北省人民医院） | 胆总管纤维化等级确定方法、装置及其相关设备 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108061687A (zh) * | 2017-12-06 | 2018-05-22 | 中南大学 | 一种研究有潜在破坏面的岩土工程模拟试验平台 |
| CN108268978A (zh) * | 2018-01-26 | 2018-07-10 | 辽宁工程技术大学 | 一种露天矿端帮边坡形态的优化方法 |
| CN109145482A (zh) * | 2018-09-04 | 2019-01-04 | 辽宁工程技术大学 | 一种软弱基底内排土场支挡煤柱形态参数优化方法 |
| CN112836283A (zh) * | 2021-03-10 | 2021-05-25 | 辽宁工程技术大学 | 一种内排土场软基底破坏与回填大块范围的确定方法 |
-
2022
- 2022-04-15 CN CN202210393657.3A patent/CN114969887B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108061687A (zh) * | 2017-12-06 | 2018-05-22 | 中南大学 | 一种研究有潜在破坏面的岩土工程模拟试验平台 |
| CN108268978A (zh) * | 2018-01-26 | 2018-07-10 | 辽宁工程技术大学 | 一种露天矿端帮边坡形态的优化方法 |
| CN109145482A (zh) * | 2018-09-04 | 2019-01-04 | 辽宁工程技术大学 | 一种软弱基底内排土场支挡煤柱形态参数优化方法 |
| CN112836283A (zh) * | 2021-03-10 | 2021-05-25 | 辽宁工程技术大学 | 一种内排土场软基底破坏与回填大块范围的确定方法 |
Non-Patent Citations (3)
| Title |
|---|
| Limit equilibrium method (LEM) of slope stability and calculation of comprehensive factor of safety with double strength-reduction technique;DENG Dong-ping 等;《J. Mt. Sci.》;全文 * |
| 露天矿支撑―回填动态作用下有限边坡力学模型研究;周永利;周伟;陆翔;才庆祥;;煤炭工程(07);全文 * |
| 露天矿顺倾软岩边坡三维稳定性分析与形态优化;曹兰柱 等;《辽宁工程技术大学学报（自然科学版）》;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114969887A (zh) | 2022-08-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114969887B (zh) | 一种确定含弱层边坡稳定性随内排支挡变化规律的方法 | |
| WO2015161833A1 (zh) | 一种固体充填采煤液压支架设计流程及优化方法 | |
| CN113255174B (zh) | 考虑岩石动态强度和混合破碎模式的钻齿力学计算方法 | |
| CN107704675B (zh) | 超前管棚-钢拱架-锁脚锚杆一体化力学模型设计方法及其模型 | |
| CN113569442A (zh) | 一种基于rkpm-pd耦合算法的岩石裂纹扩展预测方法 | |
| CN111177901A (zh) | 基于Mohr-Coulomb准则的围岩破坏危险性数值模拟评价方法 | |
| CN111259569A (zh) | 区段煤柱临空侧向复合厚硬顶板最优破断位置确定方法 | |
| CN108681630B (zh) | 基于围岩压力连续的小净距隧道深浅埋判定方法 | |
| CN103195068A (zh) | 由岩土本构折线段确定基坑挖深的方法 | |
| CN108446251B (zh) | 一种岩质边坡主动岩石压力的计算方法 | |
| CN110044725B (zh) | 一种获取平面型筋-土界面拉拔力学行为的方法 | |
| CN111651912A (zh) | 一种基于断裂力学准则的水力裂缝扩展路径数值计算方法 | |
| CN109488331B (zh) | 成型盾构隧道管片及其受力测试方法 | |
| CN108664699B (zh) | 一种静扰动加载的刀柱残采区上行开采可行性判定方法 | |
| Su et al. | Mechanical behavior and failure characteristics of double-layer composite rock-like specimens with two coplanar joints under uniaxial loading | |
| CN112131732A (zh) | 基于可变荷载作用的刀柱法残采区稳定性分析方法 | |
| CN111680434A (zh) | 包含裂纹尖端非奇异应力项的裂缝扩展路径数值计算方法 | |
| CN108979619B (zh) | 综放开采超前支承压力峰值位置的确定方法 | |
| CN102637223B (zh) | 基于伪密度排序并考虑拔模制造约束的拓扑优化设计方法 | |
| CN107818200B (zh) | 超前小导管-钢拱架-锁脚锚杆一体化力学模型设计计算方法 | |
| CN115266348B (zh) | 一种基于动静载叠加试验的岩爆倾向性等级判断方法 | |
| Chatterjee et al. | Large deformation finite element analysis of vertical penetration of pipelines in seabed | |
| CN113533042B (zh) | 一种表征岩石应力与破裂的综合性指标计算方法及应用 | |
| Qi et al. | Development characteristics of the rock fracture field in strata overlying a mined coal seam group | |
| WO2019105177A1 (zh) | 一种充填采场矿压弱化表征方法 |
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 |