CN105160165A - 一种应变软化型边坡动态稳定性评价方法 - Google Patents
一种应变软化型边坡动态稳定性评价方法 Download PDFInfo
- Publication number
- CN105160165A CN105160165A CN201510528282.7A CN201510528282A CN105160165A CN 105160165 A CN105160165 A CN 105160165A CN 201510528282 A CN201510528282 A CN 201510528282A CN 105160165 A CN105160165 A CN 105160165A
- Authority
- CN
- China
- Prior art keywords
- soil block
- side slope
- safety coefficient
- strain
- soil
- 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
Links
Abstract
本发明公开了一种基于瑞典条分法以及简化应变软化理论的边坡动态稳定性评价方法,属于地质工程技术领域。本发明将边坡竖直条分为n个土块,土块的应变作为参数,根据应变与应力的简化关系,确定土块的应力大小,再基于瑞典条分法,计算边坡的稳定性。本发明考虑了边坡的动态渐进破坏过程,根据破坏过程中土块应变发生变化,从而确定边坡的破坏过程中安全系数的变化情况。本发明原理简单,计算方便,适用于应变软化型边坡的动态稳定性分析。
Description
技术领域
本发明涉及以一种边坡稳定性分析的评价方法,特别涉及一种应变软化型边坡动态稳定性评价方法,属于地质工程技术领域。
背景技术
我国是一个山区众多的国家,边坡的不稳定而导致的滑坡会对人民安全以及经济建设造成巨大的损失。正确评价边坡的稳定性,并采取相应的加固措施,能够在很大程度上减少这种损失。对边坡的稳定性进行分析时,主要有条分法和有限元这两种方法,条分法由于概念简单、计算效率高等优点,仍然被广泛使用。
目前,在采用条分法对边坡的稳定性进行分析时,一般没有考虑到土体的应变软化特性。在采用土体应力指标的峰值强度来计算边坡稳定性时,计算结果偏小;在采用残余强度计算边坡稳定性时,计算得到的安全系数又偏大。同时,在采用条分法对边坡进行稳定性分析时,一般很少考虑边坡的动态渐进破坏过程,只是得到一个单一的安全系数。因此,在实际应用过程中,条分法还是具有一定的局限性。
发明内容
本发明的目的在于提供一种应变软化型边坡动态稳定性评价方法,本发明在使用条分法对边坡的稳定性进行分析时,充分考虑了边坡的动态渐进破坏过程,能够得到在破坏过程中一系列安全系数的变化情况;本发明考虑不同土质的软化模量对边坡稳定性的影响。
本发明采用的技术方案,一种应变软化型边坡动态稳定性评价方法,包括以下步骤:
步骤一:建立数值计算模型,确定边坡的计算参数:边坡几何大小,土体的重度γ,土体应力—应变软化关系中,土体的摩擦角黏聚力的峰值强度cp,残余强度cγ,以及峰值强度和残余强度所对应的应变γp,γγ。
步骤二:将边坡滑体从下到上分为n块,边坡下部发生剪切破坏,边坡上部同时发生拉破坏。
步骤三:边坡下部第1个土块由于应变软化开始发生剪切破坏,此时土块1的应变γi为γp<γi<γγ,取此时土块1的安全系数可表示为:
土块n由于发生拉破坏,黏聚力c和摩擦角都取0,此时土块n的安全系数为0;
其余土块2~n-1均未发生破坏,应变γi为γi≤γp,此时土块2~n-1的安全系数可表示为:
可得此时边坡整体的安全系数。
步骤四:边坡下部第1个土块的应变γi为γi≥γp,土块1的安全系数表示为
土块2的应变γi为γp<γi<γγ,取此时土块2的安全系数可表示为:
土块n和土块n-1此时均发生拉破坏,土块n和土块n-1的安全系数均为0;
其余土块3~n-2均未发生破坏,应变γi为γi≤γp,此时土块3~n-2的安全系数可表示为:
可得此时破坏后边坡的安全系数。
步骤五:根据边坡动态渐进破坏过程,边坡上部拉裂缝不再发展,边坡剪切破坏继续逐渐发展。
重复步骤四,边坡下部土块1和土块2的应变γi为γi≥γp,土块1和土块2的安全系数表示为
土块3的应变γi为γp<γi<γγ,取此时土块3的安全系数可表示为:
由于边坡拉裂缝不再发展,土块n和土块n-1此时均发生拉破坏,土块n和土块n-1的安全系数均为0;
其余土块4~n-2均未发生破坏,应变γi为γi≤γp,此时土块4~n-2的安全系数可表示为:
可得此时破坏后边坡的安全系数。
边坡上部拉裂缝不再发展,边坡剪切破坏继续逐渐发展,依次计算直至整个滑动面贯通,可得到边坡渐进破坏过程中,安全系数的变化。根据边坡的安全系数变化,判断边坡稳定的情况,从而指导边坡的加固。
本发明提供的一种应变软化型边坡动态稳定性评价方法,在条分法的基础上考虑了应变软化型边坡的渐进破坏过程,能计算出边坡破坏过程中一系列安全系数的变化,克服了条分法只有单一的安全系数的缺陷,改善了条分法单纯只使用土体应力指标的峰值强度或者残余强度来计算边坡的安全系数时,造成计算结果与实际不符的情况。该评价方法概念清晰,原理简单,计算方便,在实际工程中具有良好的应用发展前景。
附图说明
为了对本发明做进一步详细说明,下面结合附图和实施例进行介绍。
图1为本发明的流程示意图;
图2为本发明实施例简化的应力应变关系图;
图3为本发明实施例边坡分条示意图;
图4为本发明实施例边坡渐进破坏过程示意图。
图5为本发明实施例安全系数变化过程示意图。
具体实施方式
本发明旨在提供一种应变软化型边坡动态稳定性评价方法。下面结合工程实例及附图对本发明具体实施例进行详细描述。
如图1所示,本发明实施例提供的应变软化型边坡动态稳定性评价方法,包括以下步骤:
步骤一:如图2及图3所示,边坡高20m,坡角为45°。土体参数:重度为γ=18kN/m3,考虑土体的应变软化特性,摩擦角黏聚力峰值强度为cp=40kPa,对应的应变γp=1×10-3,黏聚力残余强度为cγ=10kPa,对应的应变γγ=3×10-3。
步骤二:将边坡滑体从下到上分为32条,每2个土条为1个土块,共16个土块。
步骤三:边坡下部第1个土块由于应变软化开始发生剪切破坏,此时土块1的应变γi为γp<γi<γγ,取土块1安全系数的计算可带入:
土块16由于发生拉破坏,黏聚力c和摩擦角都取0,此时土块16的安全系数为0。
其余土块2~15均未发生破坏,应变γi为γi≤γp。此时土块2~15安全系数的计算可带入:
此时边坡整体的安全系数为:Fs=1.229。
步骤四:边坡下部第1个土块的应变γi为γi≥γp,土块1安全系数的计算带入:
土块2的应变γi为γp<γi<γγ,取此时土块2安全系数的计算可带入:
土块16和土块15此时均发生拉破坏,土块16和土块15的安全系数均为0。
其余土块3~14均未发生破坏,应变γi为γi≤γp,此时土块3~14安全系数的计算可带入:
此时破坏后边坡的安全系数为:Fs=1.127。
步骤五:根据边坡动态渐进破坏过程,边坡上部拉裂缝不再发展,边坡剪切破坏继续逐渐发展。
重复步骤四,边坡下部土块1和土块2的应变γi为γi≥γp,土块1和土块2的安全系数表示为
土块3的应变γi为γp<γi<γγ,取此时土块3的安全系数可表示为:
由于边坡拉裂缝不再发展,土块16和土块15此时均发生拉破坏,土块16和土块15的安全系数均为0。
其余土块4~14均未发生破坏,应变γi为γi≤γp,此时土块4~14的安全系数可表示为:
此时破坏后边坡的安全系数Fs=1.102。
边坡上部拉裂缝不再发展,边坡剪切破坏继续逐渐发展,依次计算直至整个滑动面贯通,可得到边坡渐进破坏过程中,安全系数的变化。
计算结果列于表1:
表1安全系数的变化
工况 | 安全系数Fs |
1 | 1.229 |
2 | 1.127 |
3 | 1.102 |
4 | 1.078 |
5 | 1.053 |
6 | 1.028 |
7 | 1.002 |
8 | 1.097 |
9 | 0.947 |
10 | 0.917 |
11 | 0.888 |
12 | 0.856 |
13 | 0.821 |
14 | 0.784 |
本发明实施例中安全系数的变化过程如图5所示,渐进破坏过程中滑动面的发展如图4所示。
根据边坡的安全系数变化的示意图,判断边坡稳定的情况,从而指导边坡的加固。
Claims (1)
1.一种应变软化型边坡动态稳定性评价方法,其特征在于,包括以下步骤:
步骤一:建立数值计算模型,确定边坡的计算参数:边坡几何大小,土体的重度γ,土体应力—应变软化关系中,土体的摩擦角黏聚力的峰值强度cp,残余强度cγ,以及峰值强度和残余强度所对应的应变γp,γγ;
步骤二:将边坡滑体从下到上分为n块,边坡下部发生剪切破坏,边坡上部同时发生拉破坏;
步骤三:边坡下部第1个土块由于应变软化开始发生剪切破坏,此时土块1的应变γi为γp<γi<γγ,取此时土块1的安全系数表示为:
土块n由于发生拉破坏,黏聚力c和摩擦角都取0,此时土块n的安全系数为0;
其余土块2~n-1均未发生破坏,应变γi为γi≤γp,此时土块2~n-1的安全系数表示为:
可得此时边坡整体的安全系数;
步骤四:边坡下部第1个土块的应变γi为γi≥γp,土块1的安全系数表示为:
土块2的应变γi为γp<γi<γγ,取此时土块2的安全系数可表示为:
土块n和土块n-1此时均发生拉破坏,土块n和土块n-1的安全系数均为0;
其余土块3~n-2均未发生破坏,应变γi为γi≤γp,此时土块3~n-2的安全系数可表示为:
可得此时破坏后边坡的安全系数;
步骤五:根据边坡动态渐进破坏过程,边坡上部拉裂缝不再发展,边坡剪切破坏继续逐渐发展;
重复步骤四,边坡下部土块1和土块2的应变γi为γi≥γp,土块1和土块2的安全系数表示为
土块3的应变γi为γp<γi<γγ,取此时土块3的安全系数表示为:
由于边坡拉裂缝不再发展,土块n和土块nγ1此时均发生拉破坏,土块n和土块n-1的安全系数均为0;
其余土块4~n-2均未发生破坏,应变γi为γi≤γp,此时土块4~n-2的安全系数可表示为:
可得此时破坏后边坡的安全系数;
边坡上部拉裂缝不再发展,边坡剪切破坏继续逐渐发展,依次计算直至整个滑动面贯通,得到边坡渐进破坏过程中,安全系数的变化;根据边坡的安全系数变化,判断边坡稳定的情况,从而指导边坡的加固。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510528282.7A CN105160165B (zh) | 2015-08-25 | 2015-08-25 | 一种应变软化型边坡动态稳定性评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510528282.7A CN105160165B (zh) | 2015-08-25 | 2015-08-25 | 一种应变软化型边坡动态稳定性评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105160165A true CN105160165A (zh) | 2015-12-16 |
CN105160165B CN105160165B (zh) | 2017-11-21 |
Family
ID=54801021
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510528282.7A Active CN105160165B (zh) | 2015-08-25 | 2015-08-25 | 一种应变软化型边坡动态稳定性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105160165B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105804094A (zh) * | 2016-03-10 | 2016-07-27 | 成都理工大学 | 一种岩质滑坡的危险性判断方法及其应用 |
CN106295017A (zh) * | 2016-08-15 | 2017-01-04 | 河海大学 | 一种以变形量为失稳判据的开挖土体稳定性分析方法 |
CN107169271A (zh) * | 2017-04-26 | 2017-09-15 | 中国地质大学(武汉) | 一种应变软化边坡的稳定性动态评价方法 |
CN117688780A (zh) * | 2024-01-30 | 2024-03-12 | 中交一航局第三工程有限公司 | 一种计算二维多级堆载边坡局部和全局安全系数的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103485353A (zh) * | 2013-09-24 | 2014-01-01 | 昆明理工大学 | 基于全局最优化的边坡稳定性分析条分法 |
CN104331607A (zh) * | 2014-10-22 | 2015-02-04 | 西安长庆科技工程有限责任公司 | 一种边坡稳定分析方法 |
CN104631453A (zh) * | 2014-12-25 | 2015-05-20 | 青岛理工大学 | 边坡锚杆加固最优入射倾角的测定方法 |
-
2015
- 2015-08-25 CN CN201510528282.7A patent/CN105160165B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103485353A (zh) * | 2013-09-24 | 2014-01-01 | 昆明理工大学 | 基于全局最优化的边坡稳定性分析条分法 |
CN104331607A (zh) * | 2014-10-22 | 2015-02-04 | 西安长庆科技工程有限责任公司 | 一种边坡稳定分析方法 |
CN104631453A (zh) * | 2014-12-25 | 2015-05-20 | 青岛理工大学 | 边坡锚杆加固最优入射倾角的测定方法 |
Non-Patent Citations (3)
Title |
---|
CHEN GUO-QING ET AL.: "Progressive Modelling of the Gravity-induced Landslide Using the Local Dynamic Strength Reduction Method", 《J.MT.SCI.》 * |
陈国庆 等: "基于动态和整体强度折减法的边坡稳定性分析", 《岩石力学与工程学报》 * |
陈国庆 等: "边坡渐进破坏的动态强度折减法研究", 《岩土力学》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105804094A (zh) * | 2016-03-10 | 2016-07-27 | 成都理工大学 | 一种岩质滑坡的危险性判断方法及其应用 |
CN106295017A (zh) * | 2016-08-15 | 2017-01-04 | 河海大学 | 一种以变形量为失稳判据的开挖土体稳定性分析方法 |
CN106295017B (zh) * | 2016-08-15 | 2019-10-25 | 河海大学 | 一种以变形量为失稳判据的开挖土体稳定性分析方法 |
CN107169271A (zh) * | 2017-04-26 | 2017-09-15 | 中国地质大学(武汉) | 一种应变软化边坡的稳定性动态评价方法 |
CN107169271B (zh) * | 2017-04-26 | 2019-10-18 | 中国地质大学(武汉) | 一种应变软化边坡的稳定性动态评价方法 |
CN117688780A (zh) * | 2024-01-30 | 2024-03-12 | 中交一航局第三工程有限公司 | 一种计算二维多级堆载边坡局部和全局安全系数的方法 |
CN117688780B (zh) * | 2024-01-30 | 2024-04-12 | 中交一航局第三工程有限公司 | 一种计算二维多级堆载边坡局部和全局安全系数的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105160165B (zh) | 2017-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105160165A (zh) | 一种应变软化型边坡动态稳定性评价方法 | |
Bai et al. | A new double reduction method for slope stability analysis | |
CN105260536B (zh) | 焊趾处热点应力的计算方法 | |
CN106777520A (zh) | 一种埋入式抗滑桩的设计方法 | |
CN108627137A (zh) | 一种滑坡变形预测计算方法 | |
CN110929384A (zh) | 一种基于综采工作面的矿压大数据实时分析系统及方法 | |
CN105844423A (zh) | 一种沉管隧道运营阶段健康状态的评估方法 | |
CN103790628A (zh) | 一种综采工作面顶板灾害预警评价方法 | |
WO2021031461A1 (zh) | 一种含结构面的拱坝坝肩岩体失稳判断方法 | |
CN107657124A (zh) | 一种考虑高边坡强卸荷的预应力锚索锚固力损失计算方法 | |
CN110887749B (zh) | 结构面抗剪强度尺寸效应试验试样尺寸确定方法 | |
CN103940626A (zh) | 在役正交异性钢桥面板疲劳开裂后剩余使用寿命评估方法 | |
CN204001776U (zh) | 现浇防撞护栏钢筋定位胎架 | |
CN106845014A (zh) | 一种轨道板温度裂缝宽度的计算方法 | |
CN103218494B (zh) | 一种建立材料或土体弹塑性本构模型的方法 | |
CN110044725A (zh) | 一种获取平面型筋-土界面拉拔力学行为的方法 | |
GB2486388A (en) | Apparatus and method for predicting vertical stress fields | |
CN106529004B (zh) | 一种隧道钢筋混凝土衬砌截面钢筋面积的设计方法 | |
CN106150454B (zh) | 海上油田矢量化配产配注方法 | |
CN106677079A (zh) | 基于拱桥拱脚水平推力减小的连续拱桥加固方法 | |
Yin | A constitutive model with two yield surfaces for soils | |
CN104453807B (zh) | 一种油田注采井组水驱窜流判别方法 | |
CN103512812A (zh) | 一种粗粒土直剪试验抗剪强度取值方法 | |
Grimstad et al. | Stability analyses of quick clay using FEM and an anisotropic strain softening model with internal length scale | |
CN105095585B (zh) | 锚杆抗浮地下室整体协同设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Chen Guoqing Inventor after: Tang Peng Inventor after: Zhao Cong Inventor after: Li Le Inventor after: Huang Runqiu Inventor before: Chen Guoqing Inventor before: Zhao Cong Inventor before: Li Le Inventor before: Huang Runqiu |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |