CN111914330A - 基于图解试算法的土岩结合边坡稳定分析方法 - Google Patents

基于图解试算法的土岩结合边坡稳定分析方法 Download PDF

Info

Publication number
CN111914330A
CN111914330A CN202010787587.0A CN202010787587A CN111914330A CN 111914330 A CN111914330 A CN 111914330A CN 202010787587 A CN202010787587 A CN 202010787587A CN 111914330 A CN111914330 A CN 111914330A
Authority
CN
China
Prior art keywords
soil
rock
slope
sliding surface
safety factor
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
Application number
CN202010787587.0A
Other languages
English (en)
Other versions
CN111914330B (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.)
Hefei Municipal Design And Research Institute Co ltd
Original Assignee
Hefei Municipal Design And Research Institute Co 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 Hefei Municipal Design And Research Institute Co ltd filed Critical Hefei Municipal Design And Research Institute Co ltd
Priority to CN202010787587.0A priority Critical patent/CN111914330B/zh
Publication of CN111914330A publication Critical patent/CN111914330A/zh
Application granted granted Critical
Publication of CN111914330B publication Critical patent/CN111914330B/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/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
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

本发明公开了基于图解试算法的土岩结合边坡稳定分析方法；步骤如下：从坡脚起，采用极限平衡法得到第一最小安全系数Fs1；然后对上部土体假定圆弧滑动面采用条分法搜索最不利圆弧滑动面，并计算剩余下滑力；对下部岩体以第一最不利倾角α作图得到滑动结构面，将上滑动结构面所包含的上部土体的土条重力及其下滑力作为外荷载，得到第二最小安全系数Fs2；再从坡脚起，在岩体中对上部土体的最不利滑动面，以下部岩体内的切平面与土体的圆弧滑动面分离的岩土体为独立体，得到第三最小安全系数Fs3；取最小值作为土岩结合边坡的最终最小安全系数Fsmin。本发明符合传统土力学对土体圆弧滑动面的假定和传统岩石力学对岩体倾斜滑动面的假定。

Description

基于图解试算法的土岩结合边坡稳定分析方法
技术领域
本发明涉及岩土工程及边坡防护技术领域，特别涉及基于图解试算法的土岩结合边坡稳定分析方法。
背景技术
现有的边坡稳定分析方法中，用于土体时常假定圆弧滑动面采用条分法分析，用于岩体时常假定平面滑动面采用极限平衡法分析；当遇到土岩结合边坡，边坡实际滑动面形态显得更为复杂，在进行边坡稳定分析时，缺乏相应的理论支持。
在现有技术中,边坡稳定分析方法仅就单一的土质边坡或单一的岩质边坡,分别独立分析,对于土岩结合边坡常采将岩质边坡岩土地质参数等代为土质边坡进行圆弧滑动面分析，与岩质边坡的实际破坏形态不相符。
因此，如何使土岩结合边坡稳定分析方法更加符合传统边坡分析的力学假定，提高分析的精确性成为本领域技术人员急需解决的技术问题。
发明内容
有鉴于现有技术的上述缺陷，本发明提供基于图解试算法的土岩结合边坡稳定分析方法，实现的目的是分别将土质边坡与岩质边坡分别按圆弧滑动面、倾斜滑动面分开，并以剩余下滑力作为上述两者的力学联系，更符合传统边坡分析的力学假定。
为实现上述目的，本发明公开了基于图解试算法的土岩结合边坡稳定分析方法；步骤如下：
步骤1、从待分析的土岩结合边坡的下部岩体的坡脚起，对整个边坡采用极限平衡法进行分析，分别计算平面滑动面倾角为10°、20°、30°、35°、40°、45°、50°、60°、70°、75°、80°，以及45°+φ/2的情况下的安全系数Fs，其中φ为岩质边坡的内摩擦角，得到第一最小安全系数Fs1，以及与所述第一最小安全系数Fs1对应的第一最不利倾角α；所述第一最不利倾角α对应的平面滑动面即为第一最不利平面滑动面；
步骤2、对所述待分析的土岩结合边坡的上部土体假定圆弧滑动面采用条分法搜索最不利圆弧滑动面，并计算各土条剩余下滑力；
步骤3、对所述下部岩体以所述第一最不利倾角α作图得到滑动结构面，包含土体圆弧滑动面及岩体的倾斜滑动面，以所述滑动结构面分离的岩体为独立体，将上滑动结构面，即土体圆弧滑动面，所包含的上部土体的土条重力及其下滑力作为外荷载，再次采用极限平衡法分析，得到第二最小安全系数Fs2，以及与所述第二最小安全系数Fs2；Fs2对应倾角仍为α；
步骤4、从所述坡脚起，在岩体中对所述上部土体的最不利滑动面，即圆弧滑动面作切平面，以所述下部岩体内的切平面与土体的圆弧滑动面分离的岩土体为独立体，再次采用极限平衡法分析，得到第三最小安全系数Fs3，以及与所述第三最小安全系数Fs3；Fs3其对应倾角为β；所述第三最不利倾角β对应的平面滑动面即为第三最不利平面滑动面；
步骤5、取所述第一最小安全系数Fs1、所述第二最小安全系数Fs2、所述第三最小安全系数Fs3的最小值作为所述土岩结合边坡的最终最小安全系数Fsmin，与所述最终最小安全系数Fsmin对应的最终最不利倾角α或β对应的平面滑动面即为最终最不利平面滑动面。
优选的，所述上部土体的几何参数包括上部土质边坡的坡高H土和上部土质边坡的坡比i土；
所述上部土体的物理力学参数包括上部土质边坡的粘聚力c土和上部土质边坡的内摩擦角φ土。
优选的，所述下部岩体的几何参数包括下部岩质边坡的坡高H岩和下部岩质边坡的坡比i岩；
所述下部岩体的物理力学参数包括下部岩质边坡的粘聚力c岩和下部岩质边坡的内摩擦角φ岩。
本发明的有益效果：
本发明克服了土岩结合边坡分析时将岩质边坡岩土地质参数等代为土质边坡进行圆弧滑动面分析，与岩质边坡的传统力学假定不相符的缺点,采用岩体倾角试算与圆弧切平面作图法相结合的方法,搜索土岩结合边坡的最不利滑动面,既符合传统土力学对土体圆弧滑动面的假定,又符合传统岩石力学对岩体倾斜滑动面的假定。
以下将结合附图对本发明的构思、具体结构及产生的技术效果作进一步说明，以充分地了解本发明的目的、特征和效果。
附图说明
图1示出本发明一实施例的流程示意图。
图2示出本发明一实施例中采用试算法对下部岩体最不利平面滑动面搜索的示意图。
图3示出本发明一实施例中采用试算法对上层土体最不利圆弧滑动面搜索的示意图。
图4示出本发明一实施例中条分法土条静力分析图。
图5示出本发明一实施例中独立岩体静力分析图。
图6示出本发明一实施例中采用图解法圆弧切平面复合滑动面的示意图。
图7示出本发明一实施例中下滑力合力的示意图。
图8示出本发明一实施例中土岩独立体静力分析图。
具体实施方式
实施例
如图1所示，本发明的基于图解试算法的土岩结合边坡稳定分析方法一具体实施例的流程示意图。
具体地，如图2所示，步骤1、从下部岩体坡脚起，分别假定平面滑动面倾角为10°、20°、30°、35°、40°、45°、50°、60°、70°、75°、80°以及45°+φ/2(φ为岩体的内摩擦角)，对整个边坡采用极限平衡法进行分析，得到最不利平面滑动面(即最小安全系数Fs1对应的最不利倾角α)；
如图3所示，步骤2、对上部土体假定圆弧滑动面采用条分法搜索最不利圆弧滑动面；对下部岩体以最不利倾角α作出最不利平面滑动面；
如图4所示，步骤3、对上述最不利圆弧滑动面范围的各土条进行静力平衡分析，得到各土条的下滑力的合力在x、y方向的分量Exi、Eyi，图中第i土条重量G土i、相邻土条对第i土条的作用力Fi、Fi-1,圆弧滑动面对土条的阻滑力Ti；
如图5所示，然后，以上述平面滑动面，倾角α，分离的岩体为独立岩体，将该岩体正上方所包含的土条下滑力的总和(Ex＝∑Exi、Ey＝∑Eyi，i＝1,2…n)作为外荷载，将独立岩体重量G岩、Ex、Ey，对x’轴方向投影，求得合力即为滑动力，平面滑动面的阻滑力T与上述计算滑动力的比值即为该最不利平面滑动面的安全系数Fs2，岩体对应的平面倾角仍为α(如图5)；
如图6所示，步骤4、从坡脚对上部土体的最不利圆弧滑动面作切平面，切平面倾角β，与圆弧滑动面组成最不利复合滑动面，以最不利复合滑动面分离的岩土体为土岩独立体；
如图7所示，对复合滑动面圆弧段各土条进行静力平衡分析，得到各土条下滑力沿切平面切向(x”向)的合力Ex"和沿切平面法向(y”向)的合力Ey"；
如图8所示，对复合滑动面直线段，将上述下滑力合力Ex”、Ey”作为外荷载，将复合滑动面直线段土岩独立体重量G土、G岩、Ex”、Ey”，对x”轴方向投影，求得合力即为滑动力，平面滑动面的阻滑力T与上述计算滑动力的比值即为该最不利平面滑动面的安全系数Fs3，岩体对应的平面倾角为β。
步骤5、取所述第一最小安全系数Fs1、所述第二最小安全系数Fs2、所述第三最小安全系数Fs3的最小值作为所述土岩结合边坡的最终最小安全系数Fsmin，与所述最终最小安全系数Fsmin对应的最终最不利倾角αmin对应的平面滑动面即为最终最不利平面滑动面，其倾角为α或β。
本发明的原理如下：
首先单纯对岩体通过平面滑动倾角进行试算得到的最小安全系数Fs1,其对应倾角为α。
然后以α为基础，综合上部土体的圆弧滑动面，将圆弧滑动面作用在岩体的的滑动剩余力参与岩体的极限平衡分析，得到最小安全系数Fs2，也就是说Fs1、Fs2都是对应倾角α的。
随后从岩体坡脚开始对上述土体的圆弧滑动面作切平面,得到倾角β,再进行极限平衡分析,得到Fs3。
最后得到Fs1、Fs2、Fs3的最小值Fsmin，其对应的倾角为α或β之一。
通过上述技术手段将土质边坡圆弧滑动假定与岩质边坡极限平衡假定相结合，对上部土质边坡按圆弧滑动面进行剩余下滑力分析，对岩质边坡的破坏倾角采用作图法，并将上部土质边坡的剩余下滑力作为岩质边坡的外力，采用极限平衡法进行试算。
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思做出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。

Claims (3)

1.基于图解试算法的土岩结合边坡稳定分析方法；步骤如下：
步骤1、从待分析的土岩结合边坡的下部岩体的坡脚起，对整个边坡采用极限平衡法进行分析，分别计算平面滑动面倾角为10°、20°、30°、35°、40°、45°、50°、60°、70°、75°、80°，以及45°+φ/2的情况下的安全系数Fs，其中φ为岩质边坡的内摩擦角，得到第一最小安全系数Fs1，以及与所述第一最小安全系数Fs1对应的第一最不利倾角α；所述第一最不利倾角α对应的平面滑动面即为第一最不利平面滑动面；
步骤2、对所述待分析的土岩结合边坡的上部土体假定圆弧滑动面采用条分法搜索最不利圆弧滑动面，并计算各土条剩余下滑力；
步骤3、对所述下部岩体以所述第一最不利倾角α作图得到滑动结构面，包含土体圆弧滑动面及岩体的倾斜滑动面，以所述滑动结构面分离的岩体为独立体，将上滑动结构面，即土体圆弧滑动面，所包含的上部土体的土条重力及其下滑力作为外荷载，再次采用极限平衡法分析，得到第二最小安全系数Fs2，以及与所述第二最小安全系数Fs2；Fs2对应倾角仍为α；
步骤4、从所述坡脚起，在岩体中对所述上部土体的最不利滑动面，即圆弧滑动面作切平面，以所述下部岩体内的切平面与土体的圆弧滑动面分离的岩土体为独立体，再次采用极限平衡法分析，得到第三最小安全系数Fs3，以及与所述第三最小安全系数Fs3；Fs3其对应倾角为β；所述第三最不利倾角β对应的平面滑动面即为第三最不利平面滑动面；
步骤5、取所述第一最小安全系数Fs1、所述第二最小安全系数Fs2、所述第三最小安全系数Fs3的最小值作为所述土岩结合边坡的最终最小安全系数Fsmin，与所述最终最小安全系数Fsmin对应的最终最不利倾角α或β对应的平面滑动面即为最终最不利平面滑动面。
2.根据权利要求1所述的基于图解试算法的土岩结合边坡稳定分析方法，其特征在于，所述上部土体的几何参数包括上部土质边坡的坡高H土和上部土质边坡的坡比i土；
所述上部土体的物理力学参数包括上部土质边坡的粘聚力c土和上部土质边坡的内摩擦角φ土。
3.根据权利要求1所述的基于图解试算法的土岩结合边坡稳定分析方法，其特征在于，所述下部岩体的几何参数包括下部岩质边坡的坡高H岩和下部岩质边坡的坡比i岩；
所述下部岩体的物理力学参数包括下部岩质边坡的粘聚力c岩和下部岩质边坡的内摩擦角φ岩。
CN202010787587.0A 2020-08-07 2020-08-07 基于图解试算法的土岩结合边坡稳定分析方法 Active CN111914330B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010787587.0A CN111914330B (zh) 2020-08-07 2020-08-07 基于图解试算法的土岩结合边坡稳定分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010787587.0A CN111914330B (zh) 2020-08-07 2020-08-07 基于图解试算法的土岩结合边坡稳定分析方法

Publications (2)

Publication Number Publication Date
CN111914330A true CN111914330A (zh) 2020-11-10
CN111914330B CN111914330B (zh) 2022-06-17

Family

ID=73282973

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010787587.0A Active CN111914330B (zh) 2020-08-07 2020-08-07 基于图解试算法的土岩结合边坡稳定分析方法

Country Status (1)

Country Link
CN (1) CN111914330B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113235671A (zh) * 2021-05-07 2021-08-10 中国建筑西南勘察设计研究院有限公司 一种岩质边坡地基的承载力载荷试验系统及方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013009160A2 (en) * 2011-07-11 2013-01-17 Universiti Sains Malaysia A geometric method for predicting landslide disaste
CN105714738A (zh) * 2016-03-09 2016-06-29 华北水利水电大学 生态混凝土护坡复合结构稳定性研究方法
CN107330182A (zh) * 2017-06-28 2017-11-07 西北农林科技大学 基于增湿条件的强度折减计算边坡稳定安全系数的方法
CN109001021A (zh) * 2018-07-13 2018-12-14 中国水利水电科学研究院 采用离心模型试验评价岩土边坡安全储备的方法
CN109033681A (zh) * 2018-08-13 2018-12-18 福建省地质工程勘察院 一种边坡计算模型导入时的自动修复方法
CN109033642A (zh) * 2018-07-30 2018-12-18 青岛理工大学 一种边坡坡形改造稳定性优化测定方法
CN109359361A (zh) * 2018-09-30 2019-02-19 青岛理工大学 边坡失稳后果量化分析方法
CN109914379A (zh) * 2019-04-03 2019-06-21 昆明理工大学 一种基于条间法向力分布特征的边坡稳定性极限平衡计算方法
CN110532603A (zh) * 2019-07-22 2019-12-03 长安大学 一种黄土高边坡倾斜坡顶垂直裂隙计算方法
CN110765614A (zh) * 2019-10-24 2020-02-07 青岛理工大学 一种基于滑坡破坏形态的边坡风险综合评估方法
CN110990920A (zh) * 2019-11-25 2020-04-10 中国电建集团贵阳勘测设计研究院有限公司 抗滑桩设计中倾斜滑面滑坡剩余下滑力计算方法
CN111259473A (zh) * 2020-01-13 2020-06-09 安徽省交通勘察设计院有限公司 一种基于混合极限平衡法的涉水岸坡安全系数计算方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013009160A2 (en) * 2011-07-11 2013-01-17 Universiti Sains Malaysia A geometric method for predicting landslide disaste
CN105714738A (zh) * 2016-03-09 2016-06-29 华北水利水电大学 生态混凝土护坡复合结构稳定性研究方法
CN107330182A (zh) * 2017-06-28 2017-11-07 西北农林科技大学 基于增湿条件的强度折减计算边坡稳定安全系数的方法
CN109001021A (zh) * 2018-07-13 2018-12-14 中国水利水电科学研究院 采用离心模型试验评价岩土边坡安全储备的方法
CN109033642A (zh) * 2018-07-30 2018-12-18 青岛理工大学 一种边坡坡形改造稳定性优化测定方法
CN109033681A (zh) * 2018-08-13 2018-12-18 福建省地质工程勘察院 一种边坡计算模型导入时的自动修复方法
CN109359361A (zh) * 2018-09-30 2019-02-19 青岛理工大学 边坡失稳后果量化分析方法
CN109914379A (zh) * 2019-04-03 2019-06-21 昆明理工大学 一种基于条间法向力分布特征的边坡稳定性极限平衡计算方法
CN110532603A (zh) * 2019-07-22 2019-12-03 长安大学 一种黄土高边坡倾斜坡顶垂直裂隙计算方法
CN110765614A (zh) * 2019-10-24 2020-02-07 青岛理工大学 一种基于滑坡破坏形态的边坡风险综合评估方法
CN110990920A (zh) * 2019-11-25 2020-04-10 中国电建集团贵阳勘测设计研究院有限公司 抗滑桩设计中倾斜滑面滑坡剩余下滑力计算方法
CN111259473A (zh) * 2020-01-13 2020-06-09 安徽省交通勘察设计院有限公司 一种基于混合极限平衡法的涉水岸坡安全系数计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YUNPENG LIU: "Application of numerical simulation in the slope stability analysis of reservoir bank", 《IEEE》 *
吴柳东等: "不同极限平衡法确定边坡安全系数的比较与分析", 《水电能源科学》 *
蒋斌松等: "平面滑动边坡稳定性的解析计算", 《岩石力学与工程学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113235671A (zh) * 2021-05-07 2021-08-10 中国建筑西南勘察设计研究院有限公司 一种岩质边坡地基的承载力载荷试验系统及方法
CN113235671B (zh) * 2021-05-07 2022-12-02 中国建筑西南勘察设计研究院有限公司 一种岩质边坡地基的承载力载荷试验系统及方法

Also Published As

Publication number Publication date
CN111914330B (zh) 2022-06-17

Similar Documents

Publication Publication Date Title
Boon et al. New insights into the 1963 Vajont slide using 2D and 3D distinct-element method analyses
Ukritchon et al. Undrained limit analyses for combined loading of strip footings on clay
CN111814369B (zh) 一种可准确计算条间力倾角和边坡稳定安全系数的条分法
CN108920754B (zh) 基于应变软化剪切本构模型和位移变化的结构面控制边坡稳定性动态评价方法
CN111914330B (zh) 基于图解试算法的土岩结合边坡稳定分析方法
Little et al. Kinematics of oblique collision and ramping inferred from microstructures and strain in middle crustal rocks, central Southern Alps, New Zealand
CN106547942B (zh) 一种外倾顺层节理岩体被动侧向岩石压力的计算方法
Sarma et al. Determination of critical slip surface in slope analysis
Kumar et al. Stability determination for layered soil slopes using the upper bound limit analysis
CN106801422B (zh) 一种露天矿边坡形状优化设计方法
CN103485353A (zh) 基于全局最优化的边坡稳定性分析条分法
CN104200100B (zh) 一种基于滑动位移分析的三维边坡稳定性预测方法
Cowgill et al. The Akato Tagh bend along the Altyn Tagh fault, northwest Tibet 2: Active deformation and the importance of transpression and strain hardening within the Altyn Tagh system
Naylor The Casanova Complex of the Northern Apennines: A mélange formed on a distal passive continental margin
CN108733930A (zh) 一种反倾向边坡弯曲倾倒破坏面的搜索方法
CN108446251B (zh) 一种岩质边坡主动岩石压力的计算方法
CN112528490B (zh) 一种岩溶地基中独立基础的极限荷载的计算方法
Fredj et al. Influence of the failure surface choice on the safety factor value during slope stability studies
CN110110370B (zh) 基于极限分析下限定理的任意多面楔形体安全系数求解方法
Ismail et al. Comparison of factor of safety using different method of analysis for slope stability
CN107220470A (zh) 一种大倾角煤层工作面支架载荷的确定方法
Baker et al. Spatial distribution of safety factors: Cohesive vertical cut
Abduraimov et al. About the methods for calculation of the stability of slopes of ground massives
Holloway Cranial capacity of the hominine from Olduvai Bed I
CN109033662B (zh) 基于变形的顺层边坡开挖稳定性分析方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 230041 floors 14, 15, 16, 17 and 18 floors, Yingshang Road West and Yijing Road South, Hefei, Anhui Province

Applicant after: HEFEI MUNICIPAL DESIGN AND RESEARCH INSTITUTE Co.,Ltd.

Address before: 230041 floor 14, 15, 16, 17, 18, South riotous Nanguo office building, Yingshang Road West, Yijing Road South, Hefei City, Anhui Province

Applicant before: HEFEI MUNICIPAL DESIGN AND RESEARCH INSTITUTE Co.,Ltd.

SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant