CN110232248B - 基于传递系数法的环形凸坡稳定性评价方法 - Google Patents
基于传递系数法的环形凸坡稳定性评价方法 Download PDFInfo
- Publication number
- CN110232248B CN110232248B CN201910524272.4A CN201910524272A CN110232248B CN 110232248 B CN110232248 B CN 110232248B CN 201910524272 A CN201910524272 A CN 201910524272A CN 110232248 B CN110232248 B CN 110232248B
- Authority
- CN
- China
- Prior art keywords
- annular
- ith
- bar
- convex slope
- section
- 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
Links
Images
Classifications
-
- 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
-
- 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
-
- 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/06—Power analysis or power optimisation
Abstract
本发明公开了一种基于传递系数法的环形凸坡稳定性评价方法，其实施过程如下：将环形凸坡分为若干环形条块，获取每个环形条块的重力W i 、滑面面积A 1i ；将环形凸坡对称轴向任意一个方向拉伸形成一个剖面，每个环形条块在剖面上形成一个断面，获取剖面上每个环形条块断面面积A 2i 、滑面倾角θ i ；通过下列公式迭代计算边坡安全系数F s 。考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，对传递系数法引入新的假定进行改进，经过改进的传递系数法可以用于评价环形凸坡的稳定性，且计算过程简单，为环形凸坡的稳定性评价提供了一种计算结果更为合理的方法。
Description
技术领域
本发明涉及边坡稳定性评价方法，尤其涉及基于传递系数法的环形凸坡稳定性评价方法。
背景技术
在山区工程建设及滑坡灾害预测分析中，会遇到各种形状的边坡，如考察边坡在水平面内的形状，可将其分为凸形、凹形和直线形，边坡的空间形状对其稳定性无疑会有影响。严格来说，边坡稳定性分析属于空间问题，采用三维分析方法更符合实际情况，而工程上对边坡稳定性的评价，一般采用二维极限平衡法，这类方法对于直线形边坡计算精度尚可，但对于环形凸坡这种空间效应显著的边坡计算结果误差较大。如何分析环形凸坡的稳定性是边坡稳定性评价中需要解决的问题。
本发明以目前工程中常用的二维极限平衡分析方法传递系数法为基础，通过考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，将传递系数法进行修正以适用环形凸坡，使计算结果更符合实际情况。
发明内容
针对上述问题，本发明要解决的技术问题是提供基于传递系数法的环形凸坡稳定性评价方法，以解决现有传递系数法在评价环形凸坡稳性方面的不足。
本发明采用的技术方案为：基于传递系数法的环形凸坡稳定性评价方法，其实施过程如下：
步骤一：将环形凸坡分为若干环形条块，获取每个环形条块的重力Wi、滑面面积A1i；
步骤二：将环形凸坡对称轴向任意一个方向拉伸形成一个剖面，每个环形条块在剖面上形成一个断面，获取剖面上每个环形条块断面面积A2i、滑面倾角θi；
步骤三：通过下列公式迭代计算边坡安全系数Fs；
Pn＝0
Pi＝Pi-1ψi-1+Ti-R1i/Fs-R2i cosθi/Fs
Ti＝Wi sinθi
R2i＝2πσ拉iA2i
式中，ci为第i个环形条块的滑面粘聚力；
为第i个环形条块的滑面内摩擦角；Pn为第n个环形条块剩余下滑力；Pi为第i个环形条块与第i+1个环形条块剩余下滑力；Ti为第i个环形条块重力产生的下滑力；ψi-1为第i-1个环形条块对第i个环形条块的传递系数；R1i为第i个环形条块在滑面A1i上的抗滑剪力；R2i为第i个环形条块轴向拉力产生的抗滑力；σ拉i为第i个环形条块土体抗拉强度。
其中，步骤三中的公式基于传递系数法并考虑环形条块的轴向拉力对抗滑力的贡献，除了传递系数法的基本假定外，新引入1条假定：第i个环形条块轴向拉力产生的抗滑力R2i，作用点位于该环形条块重心上。
本发明的有益效果：考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，对传递系数法引入新的假定进行改进，经过改进的传递系数法可以用于评价环形凸坡的稳定性，且计算过程简单，为环形凸坡的稳定性评价提供了一种计算结果更为合理的方法。
附图说明
为了易于说明，本发明由下述的具体实施及附图作以详细描述。
图1为本发明实施方式中环形条块i的结构示意图；
图2为本发明实施方式中环形条块和等效环形条块的结构示意图；
图3为本发明实施例中环形凸坡三维模型计算参数图；
图4为本发明实施例中环形凸坡断面计算参数图。
具体实施方式
下面将结合本发明实施例中对本发明技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本具体实施方式基于传递系数法的环形凸坡稳定性评价方法，其实施过程如下：取图3中任意一个环形条块i，如图1所示。
由土体抗拉强度σ拉i得到环形条块土体轴向抗拉力FNi，计算公式如公式(1)所示。
FNi＝σ拉iA2i (1)
环形条块土体轴向抗拉力FNi会产生一个环形均布线荷载qi，qi方向水平且指向对称轴，qi计算公式如公式(2)所示。
公式中，ri为第i个环形条块的重力Wi到对称轴的距离。
把公式(1)带入(2)得到公式(3)。
现建立一个截面与环形条块相同且长度与环形条块周长相等的长条形条块，将其命名为等效环形条块i，并将环形均布线荷载qi施加在等效环形条块i上如图2所示，之后将等效环形条块i作为研究对象。
等效环形条块i的抗滑力R2i计算公式如公式(4)所示，R2i方向水平且指向坡内。
R2i＝2πσ拉iA2i (4)
由于等效环形条块i的抗滑力R2i的存在，使得环形凸坡的安全系数较长直边坡会更高一点。将R2i引入到传递系数法中，取整个等效环形条块i作为研究对象。
原始传递系数法计算公式如公式(5)所示
Pn＝0 (5-1)
Pi＝Pi-1ψi-1+Ti-R1i/Fs (5-2)
Ti＝Wi sinθi (5-4)
将抗滑力R2i加入到公式(5-2)中，R2i的法向分力和切向分力分别为R2i sinθi和R2icosθi得到公式(6)。
Pn＝0 (6-1)
Pi＝Pi-1ψi-1+Ti-R1i/Fs-R2i cosθi/Fs (6-2)
Ti＝Wi sinθi (6-4)
R2i＝2πσ拉iA2i (6-6)
实施例：步骤一：环形凸坡三维模型计算参数图如图3所示，共分为10个环形条块，所有环形条块的滑体容重为25kN/m3。所有环形条块的滑面粘聚力ci均为33kPa，内摩擦角
均为35°；环形条块土体抗拉强度σ拉i均为32kPa。10个环形条块的重力Wi和滑面面积A1i，如表1所示。
表1环形条块的重力Wi和滑面面积A1i
| 条块编号 | 1 | 2 | 3 | 4 | 5 |
| W<sub>i</sub>(kN) | 13975 | 44575 | 72600 | 101325 | 108000 |
| 条块编号 | 6 | 7 | 8 | 9 | 10 |
| W<sub>i</sub>(kN) | 115375 | 119550 | 100350 | 80200 | 35650 |
| 条块编号 | 1 | 2 | 3 | 4 | 5 |
| A<sub>1i</sub>(m<sup>2</sup>) | 532 | 476 | 656 | 632 | 640 |
| 条块编号 | 6 | 7 | 8 | 9 | 10 |
| A<sub>1i</sub>(m<sup>2</sup>) | 848 | 835 | 935 | 1036 | 1240 |
步骤二：环形凸坡断面计算参数图如图4所示，每个环形条块断面面积A2i、滑面倾角θi，如表2所示。
表2环形条块断面面积A2i和滑面倾角θi
步骤三：通过公式进行迭代计算安全系数Fs，
Pn＝0
Pi＝Pi-1ψi-1+Ti-R1i/Fs-R2i cosθi/Fs
Ti＝Wi sinθi
R2i＝2πσ拉iA2i
共迭代7次，每次迭代结果分别为4.80，4.90，5.0，5.05，5.052，5.054，5.056。通过迭代计算，得出最终安全系数Fs＝5.056。
本发明考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，对传递系数法引入新的假定进行改进，经过改进的传递系数法可以用于评价环形凸坡的稳定性，且计算过程简单，为环形凸坡的稳定性评价提供了一种计算结果更为合理的方法。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。
Claims (2)
1.基于传递系数法的环形凸坡稳定性评价方法，其特征在于：其实施过程如下：
步骤一：将环形凸坡分为若干环形条块，获取每个环形条块的重力Wi、滑面面积A1i；
步骤二：将环形凸坡对称轴向任意一个方向拉伸形成一个剖面，每个环形条块在剖面上形成一个断面，获取剖面上每个环形条块断面面积A2i、滑面倾角θi；
步骤三：通过下列公式迭代计算边坡安全系数Fs；
Pn＝0
Pi＝Pi-1ψi-1+Ti-R1i/Fs-R2icosθi/Fs
Ti＝Wisinθi
R2i＝2πσ拉iA2i
公式中，ci为第i个环形条块的滑面粘聚力；
为第i个环形条块的滑面内摩擦角；Pn为第n个环形条块剩余下滑力；Pi为第i个环形条块与第i+1个环形条块剩余下滑力；Ti为第i个环形条块重力产生的下滑力；ψi-1为第i-1个环形条块对第i个环形条块的传递系数；R1i为第i个环形条块在滑面A1i上的抗滑剪力；R2i为第i个环形条块轴向拉力产生的抗滑力；σ拉i为第i个环形条块土体抗拉强度。
2.根据权利要求1所述的基于传递系数法的环形凸坡稳定性评价方法，其特征在于：步骤三中的公式基于传递系数法并考虑环形条块的轴向拉力对抗滑力的贡献，除了传递系数法的基本假定外，新引入1条假定：第i个环形条块轴向拉力产生的抗滑力R2i，作用点位于该环形条块重心上。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910524272.4A CN110232248B (zh) | 2019-06-18 | 2019-06-18 | 基于传递系数法的环形凸坡稳定性评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910524272.4A CN110232248B (zh) | 2019-06-18 | 2019-06-18 | 基于传递系数法的环形凸坡稳定性评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110232248A CN110232248A (zh) | 2019-09-13 |
| CN110232248B true CN110232248B (zh) | 2020-08-07 |
Family
ID=67860055
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910524272.4A Active CN110232248B (zh) | 2019-06-18 | 2019-06-18 | 基于传递系数法的环形凸坡稳定性评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110232248B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110263422A (zh) * | 2019-06-18 | 2019-09-20 | 贵州正业工程技术投资有限公司 | 基于简单平面滑动法的环形凸坡稳定性评价方法 |
| CN110263421A (zh) * | 2019-06-18 | 2019-09-20 | 贵州正业工程技术投资有限公司 | 基于传递系数法的群拉力作用球冠型边坡稳定性评价方法 |
| CN110674550B (zh) * | 2019-09-25 | 2021-03-26 | 贵州正业工程技术投资有限公司 | 基于传递系数法的圆弧型凹坡稳定性评价方法 |
| CN110689969B (zh) * | 2019-09-25 | 2021-03-26 | 贵州正业工程技术投资有限公司 | 基于简单平面滑动法的圆弧型凹坡稳定性评价方法 |
| CN113216223A (zh) * | 2021-05-26 | 2021-08-06 | 江苏富邦环境建设集团有限公司 | 一种矿山边坡生态修复方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100577918C (zh) * | 2007-12-20 | 2010-01-06 | 马向前 | 一种治理水土流失的流域堰系及其构筑方法 |
| CN103942446B (zh) * | 2014-04-30 | 2017-02-22 | 湖北工业大学 | 基于牵引式斜坡变形破坏机理的稳定性分析和预测预警方法 |
| CN108491575A (zh) * | 2018-02-11 | 2018-09-04 | 华北水利水电大学 | 一种基于数字地形进行水库涉水边坡安全系数的计算方法 |
| CN109598013A (zh) * | 2018-09-30 | 2019-04-09 | 青岛理工大学 | 推移式滑坡最危险滑移面及抗滑桩最优桩位的确定方法 |
| CN109657291A (zh) * | 2018-11-28 | 2019-04-19 | 青岛理工大学 | 滑坡排水孔优化设计参数的测定方法 |
| CN109583135A (zh) * | 2018-12-24 | 2019-04-05 | 重庆大学 | 针对s型沟槽填方场地边坡稳定性的极限平衡分析方法 |
| CN109858128B (zh) * | 2019-01-23 | 2020-12-22 | 中国水利水电科学研究院 | 基于块分系统的边坡稳定性二维极限分析下限解法 |
-
2019
- 2019-06-18 CN CN201910524272.4A patent/CN110232248B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110232248A (zh) | 2019-09-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110232248B (zh) | 基于传递系数法的环形凸坡稳定性评价方法 | |
| CN110245429B (zh) | 基于简化Bishop法的环形凸坡稳定性评价方法 | |
| CN110197047B (zh) | 基于Janbu法的环形凸坡稳定性评价方法 | |
| CN110598353B (zh) | 基于简单平面滑动法的环形凸坡稳定性评价方法 | |
| CN110633541B (zh) | 基于Janbu法的群拉力作用球冠型边坡稳定性评价方法 | |
| CN104063903A (zh) | 三维实体模型的四面体网格生成方法和装置 | |
| CN110750873B (zh) | 基于传递系数法的群拉力作用球冠型边坡稳定性评价方法 | |
| CN110887737B (zh) | 一种复合材料加筋壁板压损强度试验确定方法 | |
| CN112067261B (zh) | 一种载荷谱的获取方法及系统 | |
| CN110598354B (zh) | 基于Janbu法的圆弧型凹坡稳定性评价方法 | |
| CN103065015A (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
| Nachmias et al. | Non-amenable Cayley graphs of high girth have p_c<p_u and mean-field exponents | |
| CN110674550B (zh) | 基于传递系数法的圆弧型凹坡稳定性评价方法 | |
| CN105608719A (zh) | 一种基于两阶段投影调整的快速ct图像重建方法 | |
| CN103366095B (zh) | 一种基于坐标变换的最小二乘拟合信号处理方法 | |
| CN106570282A (zh) | 一种办公建筑空调能耗拆分方法 | |
| CN101833597A (zh) | 用快速边界元法得到大型复杂飞行器电场分布的方法 | |
| CN104392128A (zh) | 一种三合金纳米粒子催化剂稳定结构的实现方法 | |
| CN110633542B (zh) | 基于平面滑动法的群拉力作用球冠型边坡稳定性评价方法 | |
| Turnage | Performance of soils under tire loads-application of test results to tire selection for off-road vehicles | |
| Wang | The Spatial Evolution Analysis of Manufacturing Industrial Gravity Centre of Guangdong Province | |
| CN105982644A (zh) | 一种心脏三维标测方法 | |
| CN104408238B (zh) | 一种桁式半硬壳结构轴压稳定性的数值分析方法 | |
| CN110689969B (zh) | 基于简单平面滑动法的圆弧型凹坡稳定性评价方法 | |
| Hung et al. | On the rate of convergence in limit theorems for random sums via Trotter-distance |
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 |