CN110245429B - 基于简化Bishop法的环形凸坡稳定性评价方法 - Google Patents

基于简化Bishop法的环形凸坡稳定性评价方法 Download PDF

Info

Publication number
CN110245429B
CN110245429B CN201910524263.5A CN201910524263A CN110245429B CN 110245429 B CN110245429 B CN 110245429B CN 201910524263 A CN201910524263 A CN 201910524263A CN 110245429 B CN110245429 B CN 110245429B
Authority
CN
China
Prior art keywords
annular
section
sliding surface
convex slope
annular convex
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
Application number
CN201910524263.5A
Other languages
English (en)
Other versions
CN110245429A (zh
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 GEOLOGICAL ENVIRONMENTAL MONITORING INSTITUTE
Zhengye Engineering & Investment Inc Ltd
Original Assignee
CHINA GEOLOGICAL ENVIRONMENTAL MONITORING INSTITUTE
Zhengye Engineering & Investment Inc 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 GEOLOGICAL ENVIRONMENTAL MONITORING INSTITUTE, Zhengye Engineering & Investment Inc Ltd filed Critical CHINA GEOLOGICAL ENVIRONMENTAL MONITORING INSTITUTE
Priority to CN201910524263.5A priority Critical patent/CN110245429B/zh
Publication of CN110245429A publication Critical patent/CN110245429A/zh
Application granted granted Critical
Publication of CN110245429B publication Critical patent/CN110245429B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

本发明公开了一种基于简化Bishop法的环形凸坡稳定性评价方法，其实施过程如下：将环形凸坡分为若干环形条块，获取每个环形条块的重力W i 、滑面面积A 1i ；将环形凸坡对称轴向任意一个方向拉伸形成一个剖面，每个环形条块在剖面上形成一个断面，获取剖面上每个环形条块断面面积A 2i 、滑面倾角θ i ；通过下列公式迭代计算边坡安全系数F s ；考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，对简化Bishop法引入新的假定进行改进，经过改进的简化Bishop法可以用于评价环形凸坡的稳定性，且计算过程简单，为环形凸坡的稳定性评价提供了一种计算结果更为合理的方法。

Description

基于简化Bishop法的环形凸坡稳定性评价方法
技术领域
本发明涉及边坡稳定性评价方法，尤其涉及基于简化Bishop法的环形凸坡稳定性评价方法。
背景技术
在山区工程建设及滑坡灾害预测分析中，会遇到各种形状的边坡，如考察边坡在水平面内的形状，可将其分为凸形、凹形和直线形，边坡的空间形状对其稳定性无疑会有影响。严格来说，边坡稳定性分析属于空间问题，采用三维分析方法更符合实际情况，而工程上对边坡稳定性的评价，一般采用二维极限平衡法，这类方法对于直线形边坡计算精度尚可，但对于环形凸坡这种空间效应显著的边坡计算结果误差较大。如何分析环形凸坡的稳定性是边坡稳定性评价中需要解决的问题。
本发明以目前工程中常用的二维极限平衡分析方法简化Bishop法为基础，通过考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，将简化Bishop法进行修正以适用环形凸坡，使计算结果更符合实际情况。
发明内容
针对上述问题，本发明要解决的技术问题是提供一种基于简化Bishop法的环形凸坡稳定性评价方法，以解决现有二维极限平衡法在评价环形凸坡稳定性方面的不足。
本发明采用的技术方案为：基于简化Bishop法的环形凸坡稳定性评价方法，其实施过程如下：
步骤一：将环形凸坡分为若干环形条块，获取每个环形条块的重力Wi、滑面面积A1i；
步骤二：将环形凸坡对称轴向任意一个方向拉伸形成一个剖面，每个环形条块在剖面上形成一个断面，获取剖面上每个环形条块断面面积A2i、滑面倾角θi；
步骤三：通过下列公式迭代计算边坡安全系数Fs；
Figure BDA0002097678510000021
Figure BDA0002097678510000022
Figure BDA0002097678510000023
公式中，ci为第i个环形条块的滑面粘聚力；
Figure BDA0002097678510000024
为第i个环形条块的滑面内摩擦角；σ拉i为第i个环形条块土体抗拉强度；Ti为第i个环形条块轴向拉力产生的抗滑力；Ri为第i个环形条块轴向拉力产生的抗滑力Ti到圆弧滑面圆心的力矩；R为圆弧滑面半径。
其中，步骤三中的公式基于简化Bishop法并考虑环形条块的轴向拉力对抗滑力的贡献，除了简化Bishop法的基本假定外，新引入1条假定：第i个环形条块轴向拉力产生的抗滑力Ti，其作用点位于该环形条块重心上。
本发明的有益效果：考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，对简化Bishop法引入新的假定进行改进，经过改进的简化Bishop法可以用于评价环形凸坡的稳定性，且计算过程简单，为环形凸坡的稳定性评价提供了一种计算结果更为合理的方法。
附图说明
为了易于说明，本发明由下述的具体实施及附图作以详细描述。
图1为本发明实施方式中环形条块i的结构示意图；
图2为本发明实施方式中环形条块和等效环形条块的结构示意图；
图3为本发明实施方式中等效环形条块简化Bishop法受力分析图；
图4为本发明实施例中环形凸坡三维模型计算参数图；
图5为本发明实施例中环形凸坡断面计算参数图。
具体实施方式
下面将结合本发明实施例中对本发明技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本具体实施方式基于简化Bishop法的环形凸坡稳定性评价方法，其实施过程如下：取图4中任意一个环形条块i，如图1所示。
由土体抗拉强度σ拉i得到环形条块土体轴向抗拉力FNi，计算公式如公式(1)所示。
FNi＝σ拉iA2i (1)
环形条块土体轴向抗拉力FNi会产生一个环形均布线荷载qi，qi方向水平且指向对称轴，qi计算公式如公式(2)所示。
Figure BDA0002097678510000031
公式中，ri为第i个环形条块的重力Wi到对称轴的距离。
把公式(1)带入(2)得到公式(3)。
Figure BDA0002097678510000032
现建立一个截面与环形条块相同且长度与环形条块周长相等的长条形条块，将其命名为等效环形条块i，并将环形均布线荷载qi施加在等效环形条块i上如图2所示，之后将等效环形条块i作为研究对象。
引入安全系数Fs，等效环形条块i的抗滑力Ti计算公式如公式(4)所示，Ti方向水平且指向坡内。
Figure BDA0002097678510000041
由于等效环形条块i的抗滑力Ti的存在，使得环形凸坡的安全系数较长直边坡会更高一点。将Ti引入到简化Bishop法中，取整个等效环形条块i作为研究对象，将其所受的外力投影到重心所在的剖面上如图3所示，O为圆弧滑面的圆心。
针对整个等效环形条块i，由竖向合力∑Fz＝0得到公式(5)。
Figure BDA0002097678510000042
公式中，Ni为等效环形条块i的滑面A1i上的法向力，Ti0为等效环形条块i在滑面A1i上的抗滑剪力。
针对整个等效环形条块i，由力矩求和∑MO＝0得到公式(6)。
∑WiRsinθi＝∑Ti0R+∑TiRi (6)
由摩尔库伦强度准则并引入安全系数Fs可得Ti0，如公式(7)所示。
Figure BDA0002097678510000043
将公式(5)代入公式(7)整理得到公式(8)。
Figure BDA0002097678510000044
将公式(8)代入公式(6)整理得到公式(9)。
Figure BDA0002097678510000051
Figure BDA0002097678510000052
实施例：基于简化Bishop法的环形凸坡稳定性评价方法包含以下步骤：步骤一：环形凸坡三维模型计算参数图如图4所示，共分为10个环形条块，所有环形条块的滑体容重为25kN/m3。所有环形条块的滑面粘聚力ci均为33kPa，内摩擦角
Figure BDA0002097678510000053
均为35°，所有土体抗拉强度σ拉i均为32kPa。10个环形条块的重力Wi和滑面面积A1i，如表1所示。
表1环形条块的重力Wi和滑面面积A1i
条块编号 1 2 3 4 5
W<sub>i</sub>(kN) 46125 100475 141125 102275 204125
条块编号 6 7 8 9 10
W<sub>i</sub>(kN) 194400 192525 166150 161325 117925
条块编号 1 2 3 4 5
A<sub>1i</sub>(m<sup>2</sup>) 1084 884 848 512 912
条块编号 6 7 8 9 10
A<sub>1i</sub>(m<sup>2</sup>) 816 828 800 1008 2572
步骤二：环形凸坡断面计算参数图如图5所示，每个环形条块断面面积A2i、滑面倾角θi，如表2所示。
表2环形条块断面面积A2i和滑面倾角θi
条块编号 1 2 3 4 5
A<sub>2i</sub>(m<sup>2</sup>) 5.9 13.7 20.4 15.6 33
条块编号 6 7 8 9 10
A<sub>2i</sub>(m<sup>2</sup>) 34 36.4 34 39 29.2
条块编号 1 2 3 4 5
θ<sub>i</sub>(°) 10 16 21 27 30
条块编号 6 7 8 9 10
θ<sub>i</sub>(°) 36 42 49 55 65
步骤三：通过公式进行迭代计算安全系数Fs，
Figure BDA0002097678510000061
Figure BDA0002097678510000062
Figure BDA0002097678510000063
共迭代7次，每次迭代结果分别为1.290，1.399，1.432，1.441，1.444，1.445，1.445。通过迭代计算，得到最终安全系数Fs＝1.445。
本发明考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，对简化Bishop法引入新的假定进行改进，经过改进的简化Bishop法可以用于评价环形凸坡的稳定性，且计算过程简单，为环形凸坡的稳定性评价提供了一种计算结果更为合理的方法。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。

Claims (2)

1.基于简化Bishop法的环形凸坡稳定性评价方法，其特征在于：其实施过程如下：
步骤一：将环形凸坡分为若干环形条块，获取每个环形条块的重力Wi、滑面面积A1i；
步骤二：将环形凸坡对称轴向任意一个方向拉伸形成一个剖面，每个环形条块在剖面上形成一个断面，获取剖面上每个环形条块断面面积A2i、滑面倾角θi；
步骤三：通过下列公式迭代计算边坡安全系数Fs；
Figure FDA0002559622710000011
Figure FDA0002559622710000012
Figure FDA0002559622710000013
公式中，ci为第i个环形条块的滑面粘聚力；
Figure FDA0002559622710000014
为第i个环形条块的滑面内摩擦角；σ拉i为第i个环形条块土体抗拉强度；Ti为第i个环形条块轴向拉力产生的抗滑力；Ri为第i个环形条块轴向拉力产生的抗滑力Ti到圆弧滑面圆心的力矩；R为圆弧滑面半径。
2.根据权利要求1所述的基于简化Bishop法的环形凸坡稳定性评价方法，其特征在于：所述的步骤三中的公式基于简化Bishop法并考虑环形条块的轴向拉力对抗滑力的贡献，除了简化Bishop法的基本假定外，新引入1条假定：第i个环形条块轴向拉力产生的抗滑力Ti，其作用点位于该环形条块重心上。
CN201910524263.5A 2019-06-18 2019-06-18 基于简化Bishop法的环形凸坡稳定性评价方法 Active CN110245429B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910524263.5A CN110245429B (zh) 2019-06-18 2019-06-18 基于简化Bishop法的环形凸坡稳定性评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910524263.5A CN110245429B (zh) 2019-06-18 2019-06-18 基于简化Bishop法的环形凸坡稳定性评价方法

Publications (2)

Publication Number Publication Date
CN110245429A CN110245429A (zh) 2019-09-17
CN110245429B true CN110245429B (zh) 2020-09-04

Family

ID=67887704

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910524263.5A Active CN110245429B (zh) 2019-06-18 2019-06-18 基于简化Bishop法的环形凸坡稳定性评价方法

Country Status (1)

Country Link
CN (1) CN110245429B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110263422A (zh) * 2019-06-18 2019-09-20 贵州正业工程技术投资有限公司 基于简单平面滑动法的环形凸坡稳定性评价方法
CN110489935B (zh) * 2019-06-18 2021-01-15 贵州正业工程技术投资有限公司 基于简化Bishop法的群拉力作用球冠型边坡稳定性评价方法
CN110598354B (zh) * 2019-09-25 2021-04-02 贵州正业工程技术投资有限公司 基于Janbu法的圆弧型凹坡稳定性评价方法
CN110674550B (zh) * 2019-09-25 2021-03-26 贵州正业工程技术投资有限公司 基于传递系数法的圆弧型凹坡稳定性评价方法
CN110689969B (zh) * 2019-09-25 2021-03-26 贵州正业工程技术投资有限公司 基于简单平面滑动法的圆弧型凹坡稳定性评价方法
CN111368458B (zh) * 2020-03-23 2021-04-02 青岛理工大学 一种基坑开挖边坡安全系数的计算方法
CN112597569B (zh) * 2020-12-17 2023-05-16 贵州正业工程技术投资有限公司 基于简化毕肖普法的复合地基填方边坡稳定性评价方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1598857A (zh) * 2004-09-21 2005-03-23 北京科技大学 一种优化露天矿山边坡的方法
CN101514553A (zh) * 2009-04-03 2009-08-26 重庆交通大学 基于极限平衡理论和应力分析的土坡稳定分析方法
CN106021853A (zh) * 2016-05-09 2016-10-12 青岛理工大学 一种空间变异特性下边坡稳定可靠度分析软件的开发方法
US9880081B1 (en) * 2017-03-07 2018-01-30 Ramesh Chandra Gupta Expandable jacket for triaxial, unconfined and uniaxial compression tests and test device for three-dimensional consolidation and settlement tests
CN109190593A (zh) * 2018-09-21 2019-01-11 重庆大学 基于凹凸类别划分的山区道路沿线边坡稳定性初步判别方法
CN109598013A (zh) * 2018-09-30 2019-04-09 青岛理工大学 推移式滑坡最危险滑移面及抗滑桩最优桩位的确定方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104318103B (zh) * 2014-10-23 2017-12-15 中国科学院、水利部成都山地灾害与环境研究所 一种滑坡灾害监测预警降雨阈值判定方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1598857A (zh) * 2004-09-21 2005-03-23 北京科技大学 一种优化露天矿山边坡的方法
CN101514553A (zh) * 2009-04-03 2009-08-26 重庆交通大学 基于极限平衡理论和应力分析的土坡稳定分析方法
CN106021853A (zh) * 2016-05-09 2016-10-12 青岛理工大学 一种空间变异特性下边坡稳定可靠度分析软件的开发方法
US9880081B1 (en) * 2017-03-07 2018-01-30 Ramesh Chandra Gupta Expandable jacket for triaxial, unconfined and uniaxial compression tests and test device for three-dimensional consolidation and settlement tests
CN109190593A (zh) * 2018-09-21 2019-01-11 重庆大学 基于凹凸类别划分的山区道路沿线边坡稳定性初步判别方法
CN109598013A (zh) * 2018-09-30 2019-04-09 青岛理工大学 推移式滑坡最危险滑移面及抗滑桩最优桩位的确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种考虑土条侧向强度的边坡稳定性分析新方法;孙加平等;《人民长江》;20170628;全文 *
张石高速公路某边坡稳定性分析与支护设计;刘亚峰等;《山西建筑》;20090120;全文 *

Also Published As

Publication number Publication date
CN110245429A (zh) 2019-09-17

Similar Documents

Publication Publication Date Title
CN110245429B (zh) 基于简化Bishop法的环形凸坡稳定性评价方法
CN110197047B (zh) 基于Janbu法的环形凸坡稳定性评价方法
CN110232248B (zh) 基于传递系数法的环形凸坡稳定性评价方法
CN110598353B (zh) 基于简单平面滑动法的环形凸坡稳定性评价方法
CN110633541B (zh) 基于Janbu法的群拉力作用球冠型边坡稳定性评价方法
CN110750873B (zh) 基于传递系数法的群拉力作用球冠型边坡稳定性评价方法
CN103294898A (zh) 一种计算反应堆全堆芯单棒功率的方法
CN102944872A (zh) 雷达散射截面近场到近场的变换方法
CN110887737B (zh) 一种复合材料加筋壁板压损强度试验确定方法
CN111651806B (zh) 基于惯性力法的常规输电塔设计风载荷计算方法
CN110689970B (zh) 基于简化Bishop法的圆弧型凹坡稳定性评价方法
CN110489935B (zh) 基于简化Bishop法的群拉力作用球冠型边坡稳定性评价方法
CN110598354B (zh) 基于Janbu法的圆弧型凹坡稳定性评价方法
CN107270825A (zh) 一种大圆盾构隧道直径解算方法
CN110674550B (zh) 基于传递系数法的圆弧型凹坡稳定性评价方法
CN110633542B (zh) 基于平面滑动法的群拉力作用球冠型边坡稳定性评价方法
CN101833597A (zh) 用快速边界元法得到大型复杂飞行器电场分布的方法
Imran et al. On counting polynomials of certain polyomino chains
CN104408238B (zh) 一种桁式半硬壳结构轴压稳定性的数值分析方法
CN204296786U (zh) 一种工装手推车
Rhaly Jr A comment on coposinormal operators
Hung et al. On the rate of convergence in limit theorems for random sums via Trotter-distance
CN110784888B (zh) 多层网络的识别方法及装置
CN110689969B (zh) 基于简单平面滑动法的圆弧型凹坡稳定性评价方法
Huang et al. A new spatial interpolation method based on cross-sections sampling

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