CN111414576B - 一种边坡安全系数不迭代求解方法 - Google Patents
一种边坡安全系数不迭代求解方法 Download PDFInfo
- Publication number
- CN111414576B CN111414576B CN202010251996.9A CN202010251996A CN111414576B CN 111414576 B CN111414576 B CN 111414576B CN 202010251996 A CN202010251996 A CN 202010251996A CN 111414576 B CN111414576 B CN 111414576B
- Authority
- CN
- China
- Prior art keywords
- slope
- sliding
- block
- force
- blocks
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000002689 soil Substances 0.000 claims abstract description 30
- 238000012546 transfer Methods 0.000 claims abstract description 9
- 230000003068 static effect Effects 0.000 claims abstract description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims 2
- 231100000817 safety factor Toxicity 0.000 abstract description 18
- 238000011156 evaluation Methods 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012271 agricultural production Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
- E02D17/00—Excavations; Bordering of excavations; Making embankments
- E02D17/20—Securing of slopes or inclines
Abstract
Description
技术领域
本发明属于边坡安全评价领域,尤其涉及一种边坡安全系数不迭代求解方法。
背景技术
开展边坡安全程度评价并及时进行支护是岩土工程领域最为重要的工作之一,究其原因在于:一旦发生滑坡之后,滑动的岩体和土体会对边坡下游的公路、建筑物、工农业生产设施造成不同程度的破坏。每年我国都会发生很多滑坡灾害,尤其是在雨季或者生产活动频繁的时期,不稳定因素诱发滑坡,给我国安全生产和经济活动造成了不可挽回的损失。因此,如何预先评价边坡的安全程度就显得尤为重要。
在边坡安全评价领域,通常采用极限平衡方法进行最小安全系数的搜索,利用搜索得到的最小安全系数来评价边坡安全程度。在给定滑动面之后,传统上采用迭代法来求解该滑动面的安全系数,对于某些滑动面,迭代法会无法收敛从而不能求解安全系数,这种不收敛现象会误导最小安全系数搜索的方向,从而给出不符合实际的安全系数结果。因此目前亟需一种能够合理、有效确定边坡安全系数的不迭代求解方法。
发明内容
根据以上现有技术的不足,本发明提供了一种边坡安全系数不迭代求解方法,本方法能够合理、有效地确定边坡安全系数。
本发明解决的技术问题采用的技术方案为:
本发明提供了一种边坡安全系数不迭代求解方法,包括如下步骤:
步骤2、假定相邻块体间的作用力方向,记作用力Ti-1为块体Bi-1对块体Bi的作用力,作用力Ti为块体Bi+1对块体Bi的作用力,通过静力平衡条件构建相邻块体间作用力的传递方程Ti=f(Ti-1);
步骤5、从集合F'中选择最大值的根作为可行滑动面S的安全系数。
其中,优选方式为:
步骤1的具体实现过程包括如下步骤:
1.1、根据勘察报告建立边坡几何模型,确定边坡高度H,边坡坡角β,边坡坡面线p(x),并查询获取边坡土层的物理力学参数;
步骤2的具体实现过程为:
块体Bi的底部与水平线夹角记为αi,块体Bi上作用有自重Wi,块体Bi底部中心点处的滑床反力记为Ni、块体抗滑力记为Ri,将自重Wi、滑床反力Ni、块体抗滑力Ri分别向滑床反力Ni、块体抗滑力Ri方向上投影,构建块体间作用力的传递方程:Ti=f(Ti-1)。
本发明具有以下有益效果:本方法通过将滑动土体分成若干垂直块体,并逐一分析块体间力的平衡条件,进而通过边界条件建立关于安全系数的多项式方程,通过逐一判断方程的根是否满足块体力矩平衡条件最终确定安全系数,避免了现有技术中求取安全系数时不能收敛导致求取不准确的问题,求得的安全系数结果有效、合理。
附图说明
图1是本发明的流程框图;
图2是本发明所提供实施例中滑动块体作用力示意图;
图3是本发明所提供实施例中某典型边坡示意图;
图4是本发明所提供实施例中边坡划分块体示意图;
图5是块体B1作用力详细图;
图6是块体B2作用力详细图。
具体实施方式
下面结合附图对本发明做进一步描述。
实施例一:
如图1~图6所示,本发明所提供的一种边坡安全系数不迭代求解方法,包括如下步骤:
步骤1、将边坡的滑动土体划分成n个的垂直块体,具体的步骤为:
1.2针对一可行的滑动面S,滑动面S与边坡坡面线p(x)构成的滑动土体,自坡面线p(x)转折点向下做垂直线,止于滑床线,分别确定垂直线段的上下端点,将滑动土体划分成n个的垂直块体,记为i的值沿着滑动方向增加,n为正整数,1≤i≤n。
步骤2、假定相邻块体间的作用力方向,记Ti-1为块体Bi-1对块体Bi的作用力,Ti为块体Bi+1对块体Bi的作用力,通过静力平衡条件构建相邻块体间作用力的传递方程,具体步骤为:
块体Bi的底部与水平线夹角记为αi,块体Bi上作用有自重Wi,块体Bi底部中心点处的滑床反力记为Ni、块体Bi的抗滑力记为Ri,将自重Wi、滑床反力Ni、块体抗滑力Ri分别向滑床反力Ni、块体抗滑力Ri方向上投影,构建块体间作用力的传递方程:
步骤5、从集合F'中选择最大值的根作为可行滑动面S的安全系数。
本方法通过将滑动土体分成若干垂直块体,并逐一分析块体间力的平衡条件,进而通过边界条件建立关于安全系数的多项式方程,通过逐一判断方程的根是否满足块体力矩平衡条件最终确定安全系数,避免了现有技术中求取安全系数时不能收敛导致求取不准确的问题,求得的安全系数结果有效,过程合理。
以下是本发明所述方法在具体场景中的应用,以证明本方法的有效性:
根据业界常用的滑动面生成方法,生成可行滑动面S如图4所示,滑动面S与坡面线p(x)构成了滑动体。为简单起见,自坡面线转折点(a)向下做垂直线,止于滑床线,分别确定垂直线段的上、下端点为(a)、(c),将滑动体划分成B1和B2两个块体,如图4中所示,B1和B2两个块体之间用竖直虚线分隔。B1块体的自重W1=194.5kN,底部与水平线夹角α1=58.8°,B1块体底部粘聚力总和C1=198.8kN,B2块体的自重W2=284.6kN,底部与水平线夹角α2=40.8°,B2块体底部粘聚力总和C2=198.8kN。根据图2,将Bi块体上的作用力向滑床反力Ni方向上投影得到,Ni=Wi cosαi+Ti-1 sin(αi-1-αi),假设块体底部抗滑反力Fs为滑动面S的安全系数。将Bi块体上的作用力向Ti方向上投影,得到作用力Ti的传递方程:
根据滑动土体的边界条件T0=0,得到:
依次得到:
根据滑动土体的边界条件T2=0,得到关于Fs的多项式方程如下:
344.2×Fs2-519.8Fs+26.5=0,求解该方程,得到方程的两个根Fs1=1.46,Fs2=0.05,分别将这两个根代入块体力矩平衡条件,验证是否满足。
验证过程具体如下:如图5、6所示,B1块体上作用有W1、T1、N1、R1,块体三个顶点(a)、(b)、(c)以及W1作用点(d)、T1作用点(e)以及N1、R1作用点(f)的坐标见表1。
表1 B1块体控制点及作用点详细信息
(a) | (b) | (c) | (d) | (e) | (f) | |
x | 10.0 | 13.4 | 10.0 | 11.1 | 10.0 | 可调 |
y | 15.0 | 15.0 | 9.3 | 13.1 | 可调 | 可调 |
B2块体上作用有W2、T1、N2、R2,块体三个顶点(a)、(g)、(c)以及W2作用点(h)、T1作用点(e)以及N2、R2作用点(j)的坐标见表2。
表2 B2块体控制点及作用点详细信息
(a) | (c) | (g) | (h) | (j) | (e) | |
x | 10.0 | 10.0 | 5.0 | 8.3 | 可调 | 10.0 |
y | 15.0 | 9.3 | 5.0 | 9.8 | 可调 | 可调 |
以Fs1=1.46为例,对于B1、B2块体,分别对(f)、(j)点取矩,通过调整点(f)的位置,(在线段(b)(c)上移动)、点(j)的位置(在线段(g)(c)上移动)以及点(e)的位置(在线段(a)(c)上移动)能够满足力矩平衡条件,因此Fs1=1.46保留。以B1为例,先假设点(f)位于线段(b)(c)中点,根据B1力矩平衡确定(e)点的位置,然后再通过B2的力矩平衡确定点(j)的位置,若确定的(e)点位置不在点(a)、(c)之间或者点(j)的位置不在(g)、(c)之间则称之为不合理的(e)点或(j)点,则改变点(f)的位置,直至找到合理的(e)或者(j),若无法找到合理的(e)点或者(j)点,则视为不满足力矩平衡条件。进行相同操作,发现Fs2=0.05不满足力矩平衡条件,因此舍弃。最终选择1.46作为滑动面安全系数,如果划分的块体数量更多时,计算方法与以上两个块体时的方法相同,此处不再赘述。
为了证明本发明方法的有效性,现进行对比分析,传统的边坡安全评价方法中,推荐使用摩根斯坦法和斯宾塞法计算可行滑动面的安全系数。利用摩根斯坦法和斯宾塞法计算本例中滑动面S的安全系数,结果发现:迭代均不收敛,无法给出有意义的安全系数,在最小安全系数搜索时,会摒弃滑动面S,因此会误导搜索的方向。本发明方法利用块体作用力传递方程,并结合滑动土体的边界条件建立关于安全系数的多项式方程。通过方程根与块体力矩平衡的综合比较,可以得出合理的安全系数。因此,对比发现:采用迭代策略求解安全系数的传统方法存在不收敛的可能性,这种可能性会误导最小安全系数的搜索方向,不利于边坡稳定安全系数的计算,从而影响滑坡风险的防治。通过实例对比分析验证了本发明的有效性。
以上所述为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书以及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (3)
1.一种边坡安全系数不迭代求解方法,其特征在于,包括如下步骤:
步骤2、假定相邻块体间的作用力方向,记作用力Ti-1为块体Bi-1对块体Bi的作用力,作用力Ti为块体Bi+1对块体Bi的作用力,通过考虑块体自重、滑床反力、块体抗滑力在滑床反力方向和块体抗滑力方向上的静力平衡条件构建相邻块体间作用力的传递方程Ti=f(Ti-1);
步骤5、从集合F'中选择最大值的根作为可行滑动面S的安全系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010251996.9A CN111414576B (zh) | 2020-04-01 | 2020-04-01 | 一种边坡安全系数不迭代求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010251996.9A CN111414576B (zh) | 2020-04-01 | 2020-04-01 | 一种边坡安全系数不迭代求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111414576A CN111414576A (zh) | 2020-07-14 |
CN111414576B true CN111414576B (zh) | 2021-03-30 |
Family
ID=71491794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010251996.9A Expired - Fee Related CN111414576B (zh) | 2020-04-01 | 2020-04-01 | 一种边坡安全系数不迭代求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111414576B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101514553A (zh) * | 2009-04-03 | 2009-08-26 | 重庆交通大学 | 基于极限平衡理论和应力分析的土坡稳定分析方法 |
CN103485353A (zh) * | 2013-09-24 | 2014-01-01 | 昆明理工大学 | 基于全局最优化的边坡稳定性分析条分法 |
CN107330224A (zh) * | 2017-07-24 | 2017-11-07 | 中国地质大学(武汉) | 一种分条间作用力倾角非假定的边坡稳定性分析条分法 |
CN109446616A (zh) * | 2018-10-18 | 2019-03-08 | 东北电力大学 | 一种均质边坡可靠性评价方法 |
CN110083883A (zh) * | 2019-04-04 | 2019-08-02 | 中国水利水电科学研究院 | 基于极限分析下限定理的三维边坡安全系数求解方法 |
CN110222405A (zh) * | 2019-05-30 | 2019-09-10 | 长安大学 | 一种基于降低强度包线折减法的岩质边坡稳定性测试方法 |
WO2019176835A1 (ja) * | 2018-03-13 | 2019-09-19 | 日本電気株式会社 | 斜面監視システム、斜面監視方法及び記録媒体 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110378056A (zh) * | 2019-07-25 | 2019-10-25 | 四川大学 | 一种用于边坡地质力学模型的边坡稳定性测定方法及应用 |
-
2020
- 2020-04-01 CN CN202010251996.9A patent/CN111414576B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101514553A (zh) * | 2009-04-03 | 2009-08-26 | 重庆交通大学 | 基于极限平衡理论和应力分析的土坡稳定分析方法 |
CN103485353A (zh) * | 2013-09-24 | 2014-01-01 | 昆明理工大学 | 基于全局最优化的边坡稳定性分析条分法 |
CN107330224A (zh) * | 2017-07-24 | 2017-11-07 | 中国地质大学(武汉) | 一种分条间作用力倾角非假定的边坡稳定性分析条分法 |
WO2019176835A1 (ja) * | 2018-03-13 | 2019-09-19 | 日本電気株式会社 | 斜面監視システム、斜面監視方法及び記録媒体 |
CN109446616A (zh) * | 2018-10-18 | 2019-03-08 | 东北电力大学 | 一种均质边坡可靠性评价方法 |
CN110083883A (zh) * | 2019-04-04 | 2019-08-02 | 中国水利水电科学研究院 | 基于极限分析下限定理的三维边坡安全系数求解方法 |
CN110222405A (zh) * | 2019-05-30 | 2019-09-10 | 长安大学 | 一种基于降低强度包线折减法的岩质边坡稳定性测试方法 |
Non-Patent Citations (2)
Title |
---|
Stability analysis of statically indeterminate blocks in key block theory and application to rock slope in Jinping-I Hydropower Station;Guanhua Sun 等;《Engineering Geology》;20150224;第186卷;第57-67页 * |
基于整体稳定性分析法的边坡临界滑动面搜索;孙聪 等;《岩土力学》;20130930;第34卷(第9期);第2853-2857页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111414576A (zh) | 2020-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111368458B (zh) | 一种基坑开挖边坡安全系数的计算方法 | |
CN111445156B (zh) | 基于变权重模糊综合评价的偏压隧道施工安全评价方法 | |
CN110335352B (zh) | 一种机载激光雷达点云的双基元多分辨率层次滤波方法 | |
WO2020173253A1 (zh) | 页岩气压裂水平井簇间距优化方法 | |
CN111414576B (zh) | 一种边坡安全系数不迭代求解方法 | |
CN105913332A (zh) | 一种确定致密气藏开发指标的方法及装置 | |
CN108548730A (zh) | 基于传递系数法和表面位移的滑坡稳定性评价方法 | |
CN110633541B (zh) | 基于Janbu法的群拉力作用球冠型边坡稳定性评价方法 | |
CN106530397A (zh) | 一种基于稀疏剖面地质轮廓线的地质面三维重建方法 | |
CN106530396A (zh) | 一种矿山地质编录数据三维处理方法 | |
CN112489209B (zh) | 一种崩塌滚石灾害情景再现方法 | |
CN109033642B (zh) | 一种边坡坡形改造稳定性优化测定方法 | |
CN112069579B (zh) | 基于dem数字地形分析的土石坝变形震害定量评估方法 | |
CN107886573A (zh) | 一种复杂地质条件下边坡三维有限元网格生成方法 | |
CN108446251A (zh) | 一种岩质边坡沿陡倾的外倾软弱结构面破坏时的主动岩石压力的计算方法 | |
CN110599595B (zh) | 一种地质构造的三维表征方法 | |
CN112116709A (zh) | 一种提高地形表达精度的地形特征线处理方法 | |
CN109783991B (zh) | 一种无已知底滑面的滑坡滑动过程模拟方法 | |
CN104007249A (zh) | 一种岩体应力状态分析方法 | |
CN110457772A (zh) | 一种结合平面曲率和最陡下坡方向的dem流向估计方法 | |
TWI534340B (zh) | Tilt Angle Transformation and Energy Dissipation Structure of Semi - automatic Tunnel | |
CN106203681A (zh) | 一种数据驱动的公共交通导向型城市区域划分与挖掘方法 | |
CN115374625A (zh) | 地震滑坡危险性检测方法、装置及介质 | |
CN112016210B (zh) | 形状不规则坚硬物体冲击接触点搜寻与冲击力测量方法、山洪泥石流大块石冲击力测量方法 | |
CN104376601B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210330 |
|
CF01 | Termination of patent right due to non-payment of annual fee |