CN103745131A - 边坡稳定性精细化评价方法 - Google Patents
边坡稳定性精细化评价方法 Download PDFInfo
- Publication number
- CN103745131A CN103745131A CN201410043002.9A CN201410043002A CN103745131A CN 103745131 A CN103745131 A CN 103745131A CN 201410043002 A CN201410043002 A CN 201410043002A CN 103745131 A CN103745131 A CN 103745131A
- Authority
- CN
- China
- Prior art keywords
- slope
- slope stability
- thrust
- broad sense
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000004458 analytical method Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000009466 transformation Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 2
- 230000001932 seasonal effect Effects 0.000 description 2
- 229910000571 Auricupride Inorganic materials 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000007616 round robin method Methods 0.000 description 1
Images
Landscapes
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
本发明提出了边坡稳定性精细化评价方法,属于边坡工程领域,其解决了边坡稳定性评价结果的可靠性问题,同时提高了边坡稳定性评价结果的精度。该方法主要应用广义将常规极限平衡评价方法,通过数学变换转化为任何状态下的广义剩余推力方法,在此基础上,基于线性规划理论、有限元学术方法和数值分析方法,建立微小单元体系,进行边坡稳定性的精细化评价,从而提高边坡稳定性评价结果的可靠性,真实反映边坡的安全状态。
Description
技术领域
本发明涉及边坡工程领域,具体涉及边坡稳定性评价技术与安全控制技术。
背景技术
边坡稳定性的评价结目前主要以刚体极限平衡法的计算结果为主要依据。对大型露天矿深部资源开发所形成的高边坡而言,常规计算方法不能将条块沿横竖两向都划分为微尺寸,由此,条块划分尺度大,导致计算结果精度低且结果非唯一,再加上岩体强度参数的不确定性,直接影响评价结果的可靠性,从而边影响到边坡设计的安全度与正常生产,因此,评价结果的可靠性亟待提高。
发明内容
本发明将稳定系数与下滑力之积视为给定安全状态条件下的下滑力,由此,任何边坡极限平衡评价方法都可以表征为下滑力和抗滑力,两者之差称为给定安全系数条件下的广义剩余推力。依掘有限元法,将边坡剖面划分为有限个微单元体,应用广义剩余推力法依次判别各个单元体滑动破坏的可能性及其安全性,然后应用线性规划理论和数值分析方法,构建边坡稳定性评价体系,从而可以同时求解边坡整体和局部的最小安全系数及其对应的最危险滑面,使计算结果唯一,由此,解决了高边坡稳定性评价结果因人而异的技术难题。当单元体划分的尺寸足够小时,评价结果无限接近现有评价理论的真值,提高了边坡稳定计算结果的可靠性和精准度。
首先,建立广义剩余推力法,△F=K*F1-F2,式中:△F为广义剩余推力;K为边坡设定的稳定系数;F1为抗滑力总和函数式,F2为下滑力总和函数式,且F1、F2对于不同的极限平衡计算方法,其计算函数式是不同的。
其次,基于有限元方法的学术思想,将边坡体划分为有限个微小单元体,线性规划方法和数值分析方法,解决精细化计算问题。具体某一过程计算过程如下:
依据临界滑移场原理,在某限定高度极限坡角时对应的Pmax=0,即P=f(θ*)=0时,θ*便为方程P=f(θ)=0的数值解,θ*值即为极限坡角。计算时先划分有解的区间(θa,θb),让其范围尽可能的大,一般θa,θb∈(0,90°),且要满足Pa=f(θa)<0,Pb=f(θb)>0,由数值解的存在定理知,极限坡角θ*必存在于(θa,θb)之间,取中点θ0=θa+θb/2将其分为均分,利用临界滑移场技术求出对应坡角θ0时的最大剩余推力为P0;若P0≠0,就是说θ0不是f(θ)的解,即不是极限坡角。继续搜索,检查P0是否与Pa同号,如果为同号,则数值解θ*(极限坡角)必大于θ0,这时令θa1=θ0,θb1=θb;否则,θ*必小于θ0;这时令θa1=θa,θb=θ0。不论出现何种情况,数值解的值域(θa1,θb1)仅为区间(θa,θb)前一半,而对缩小的数值解区间(θa1,θb1)可实行相同方法,进而得到系列有数值解区间 其中(θan,θbn)的中点为θn=θan+θbn/2,长度为θan-θbn=θa-θb/2n+1;也就是说n足够大时,最终区间必收缩于某一点θ*,从而求解出极限坡角。其误差按式(1)计算。
|θ*-θn|≤(θbn-θan)/2 (1)
对于通常的(0,90°)范围,迭代9次后误差大约为:
(90°-0)/29+1=90°/1024≈0.089°
计算精度比较精确,若初始的角度范围划分较小且包含极限边坡角,则更能快速的收敛;若给定坡底境界线,可同时求出极坡高。从而解决了精细化评价中某一计算步骤的计算结果,采用循环方法可以解决系统计算问题。
附图说明
图1为精细化评价单元体划分图示,b为条块跨度;B1B2为滑动面入口范围区;A2A1为滑动面出口范围区;(i,i)为状态点编号;
图2为48°坡角时边坡危险滑动场图示,坡面的数字为剩余推力;
图3为24°坡角时边坡危险滑动场图示;
图4为39.375°坡角时边坡临界滑动场图示。
具体实施方式
本方法的具体实施方法,可以结合附图和实例说明如下。
厦门某金铜矿露天采场边坡较大,东帮最高标高+1036m,露天采场底标高+148m,最大高差达880m。因此,保证矿山的安全生产既是前提也是不影响周边环境的关键部分。在设计给定安全系数Fs=1.15及坑底境界线不可变动的条件下,找到合理的极限边坡角非常困难,本专利提供的方法攻克了这一难点。
边坡精细化评价方法的详细步骤如下:
1)原设计48°边坡稳定性评价。将原设计边坡地质参数导入该方法的计算程序,利用广义剩余推力法理论(图1)求得边坡的危险滑动场如图2。由搜寻结果可知该坡角下最大剩余推力+393975kN,远大于零,边坡稳定状态不满足1.15要求,需要减缓坡角,以满足设计要求。
2)利用该方法,取坡角为24°时的危险滑移场如图3所示,无危险滑移面,说明最大剩余推力小于零,故坡角设计偏于保守。
3)根据该方法,依次对坡角进行优化,可以得到36°,42°,39°,40.5°,39.75°,39.375°时对应的危险滑移场。当坡角为39.375°时(图4),边坡最大生剩余推力接近零,当然可以取39.375°与39°之间的中点坡角继续迭代,直到最大剩余推力等于零为止,但由于剩余推力非常小,在这里取39°边坡作为安全系数Fs=1.15时的临界边坡,即极限坡角为39°,由此解决了给定安全系数和工程地质条件下的边坡优化设计问题。
Claims (3)
1.边坡稳定性精细化评价方法,其特征在于,该方法包括以下步骤:
S1,建立广义剩余推力计算方法;
S2,建立微小单元体系,构建精细化评价的计算过程。
2.如权利要求1所述的方法,其特征在于:
步骤S1中,建立广义剩余推力法,△F=K*F1-F2,式中:△F为广义剩余推力;K为边坡设定的稳定系数;F1为抗滑力总和函数式,F2为下滑力总和函数式,且F1、F2对于不同的极限平衡计算方法,其计算函数式是不同的。
3.如权利要求1所述的方法,其特征在于:
步骤S2中,基于有限元方法,将边坡体划分为有限个微小单元体,采用线性规划方法和数值分析方法,解决精细化计算问题。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410043002.9A CN103745131A (zh) | 2014-01-27 | 2014-01-27 | 边坡稳定性精细化评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410043002.9A CN103745131A (zh) | 2014-01-27 | 2014-01-27 | 边坡稳定性精细化评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103745131A true CN103745131A (zh) | 2014-04-23 |
Family
ID=50502148
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410043002.9A Pending CN103745131A (zh) | 2014-01-27 | 2014-01-27 | 边坡稳定性精细化评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103745131A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104715145A (zh) * | 2015-03-04 | 2015-06-17 | 中国科学院力学研究所 | 一种边坡的稳定安全性的分析方法 |
CN105178336A (zh) * | 2015-08-18 | 2015-12-23 | 辽宁工程技术大学 | 一种软弱基底排土场极限承载力的确定方法 |
CN105442620A (zh) * | 2015-11-17 | 2016-03-30 | 青岛理工大学 | 局部滑移边坡的抗滑桩桩位优化设计方法 |
CN105678012A (zh) * | 2016-01-27 | 2016-06-15 | 青岛理工大学 | 一种边坡变宽度条分的下滑推力评价方法 |
CN108665146A (zh) * | 2018-04-16 | 2018-10-16 | 浙江大学 | 大型露天矿山边坡稳定性评价精度确定方法 |
CN109115994A (zh) * | 2018-08-27 | 2019-01-01 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种斜坡震裂土体分块稳定性计算方法 |
CN110018292A (zh) * | 2019-04-19 | 2019-07-16 | 中国地质调查局西安地质调查中心 | 基于水-吸应力-局部安全场的黄土滑坡精细预警技术 |
CN110457855A (zh) * | 2019-08-21 | 2019-11-15 | 成都理工大学 | 部分淹没c′-φ′均质边坡临界库水位分析的解析方法 |
WO2020192118A1 (zh) * | 2019-03-28 | 2020-10-01 | 青岛理工大学 | 一种边坡滑动面积的评估方法 |
CN111985021A (zh) * | 2020-06-22 | 2020-11-24 | 大连海事大学 | 一种盾构开挖面的遍布节理流固耦合的安全度分析方法 |
CN117435891A (zh) * | 2023-12-20 | 2024-01-23 | 成都嘉新科技集团有限公司 | 一种土质边坡的稳定性模拟评估方法及系统 |
-
2014
- 2014-01-27 CN CN201410043002.9A patent/CN103745131A/zh active Pending
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104715145B (zh) * | 2015-03-04 | 2017-09-01 | 中国科学院力学研究所 | 一种边坡的稳定安全性的分析方法 |
CN104715145A (zh) * | 2015-03-04 | 2015-06-17 | 中国科学院力学研究所 | 一种边坡的稳定安全性的分析方法 |
CN105178336A (zh) * | 2015-08-18 | 2015-12-23 | 辽宁工程技术大学 | 一种软弱基底排土场极限承载力的确定方法 |
CN105178336B (zh) * | 2015-08-18 | 2017-02-01 | 辽宁工程技术大学 | 一种软弱基底排土场极限承载力的确定方法 |
CN105442620A (zh) * | 2015-11-17 | 2016-03-30 | 青岛理工大学 | 局部滑移边坡的抗滑桩桩位优化设计方法 |
CN105442620B (zh) * | 2015-11-17 | 2017-05-31 | 青岛理工大学 | 局部滑移边坡的抗滑桩桩位优化设计方法 |
CN105678012A (zh) * | 2016-01-27 | 2016-06-15 | 青岛理工大学 | 一种边坡变宽度条分的下滑推力评价方法 |
CN105678012B (zh) * | 2016-01-27 | 2018-07-27 | 青岛理工大学 | 一种边坡变宽度条分的下滑推力评价方法 |
CN108665146A (zh) * | 2018-04-16 | 2018-10-16 | 浙江大学 | 大型露天矿山边坡稳定性评价精度确定方法 |
CN108665146B (zh) * | 2018-04-16 | 2022-07-29 | 浙江大学 | 大型露天矿山边坡稳定性评价精度确定方法 |
CN109115994B (zh) * | 2018-08-27 | 2021-03-12 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种斜坡震裂土体分块稳定性计算方法 |
CN109115994A (zh) * | 2018-08-27 | 2019-01-01 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种斜坡震裂土体分块稳定性计算方法 |
AU2019438357B2 (en) * | 2019-03-28 | 2021-12-09 | Qingdao university of technology | Slope sliding area assessment method |
WO2020192118A1 (zh) * | 2019-03-28 | 2020-10-01 | 青岛理工大学 | 一种边坡滑动面积的评估方法 |
CN110018292A (zh) * | 2019-04-19 | 2019-07-16 | 中国地质调查局西安地质调查中心 | 基于水-吸应力-局部安全场的黄土滑坡精细预警技术 |
CN110457855A (zh) * | 2019-08-21 | 2019-11-15 | 成都理工大学 | 部分淹没c′-φ′均质边坡临界库水位分析的解析方法 |
CN111985021A (zh) * | 2020-06-22 | 2020-11-24 | 大连海事大学 | 一种盾构开挖面的遍布节理流固耦合的安全度分析方法 |
CN111985021B (zh) * | 2020-06-22 | 2024-04-19 | 大连海事大学 | 一种盾构开挖面的遍布节理流固耦合的安全度分析方法 |
CN117435891A (zh) * | 2023-12-20 | 2024-01-23 | 成都嘉新科技集团有限公司 | 一种土质边坡的稳定性模拟评估方法及系统 |
CN117435891B (zh) * | 2023-12-20 | 2024-02-27 | 成都嘉新科技集团有限公司 | 一种土质边坡的稳定性模拟评估方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103745131A (zh) | 边坡稳定性精细化评价方法 | |
Lazzaroni et al. | A model for crack propagation based on viscous approximation | |
Wang et al. | Boundary layers in incompressible Navier-Stokes equations with Navier boundary conditions for the vanishing viscosity limit | |
Obukhov et al. | Conservation laws in gravity: A unified framework | |
Jinzhou et al. | Applicable conditions and analytical corrections of plane strain assumption in the simulation of hydraulic fracturing | |
Wang et al. | An integral boundary value problem for nonlinear differential equations of fractional order on an unbounded domain | |
BR112015031904A2 (pt) | abordagem híbrida para casamento de histórico auxiliado em reservatórios grandes | |
Sakyi | A contrast between indirect boundary element and finite element methods. | |
Zhou et al. | Relationship between JZ25-1S transfer zone and hydrocarbon accumulation in Liaodongwan offshore of Bohai Bay Basin | |
Nishitani et al. | Drift displacement data based estimation of cumulative plastic deformation ratios for buildings | |
Xue-Song et al. | Theoretical analysis and numerical simulation of variable curvature friction pendulum isolation bearing | |
Wang et al. | Sensitivity analysis of slope stability influence factors based on BP neural network | |
CN102519457B (zh) | 基于角速度的欧拉角傅里埃近似输出方法 | |
CN103162688A (zh) | 基于角速度和勒让德近似的惯性测量单元欧拉角输出方法 | |
Cao et al. | Implicit expression solution of slope reliability on kinematical element method and response surface method | |
Lu et al. | Non-probabilistic reliability research on slope stability based on interval analysis | |
祝俊 et al. | A direction finding method for passive radar seeker using rolling phase interferometer | |
Tarasenko et al. | THE USE OF STATISTICAL CRITERIA FOR EVALUATION TEST OF DATA ON PROPERTIES OF INORGANIC SUBSTANCES | |
Li et al. | Another Solving Equation of Elastoplastic Problem Based on Stress Strain Relation | |
Seyedin et al. | Using Mutual Aggregate Uncertainty Measures in a Threat Assessment Problem Constructed by Dempster–Shafer Network | |
CN102506873A (zh) | 基于角速度的欧拉角拉盖尔近似输出方法 | |
Zamani et al. | Computation of Stress intensity factors of collinear interacting cracks in a simple lap joint structure | |
程帅 et al. | Impact analysis on finite element calculation error of the forced state of dam base surface | |
Zhong | A Symplectic Approach for Attitude Optimal Control of Spacecraft Solar Array Deployment | |
Irvash | The Metaphor of Vaghf (Dedication) about the Government |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140423 |
|
WD01 | Invention patent application deemed withdrawn after publication |