基于Janbu法的环形凸坡稳定性评价方法与流程

本发明涉及边坡稳定性评价方法，尤其涉及基于janbu法的环形凸坡稳定性评价方法。
背景技术：
：在山区工程建设及滑坡灾害预测分析中，会遇到各种形状的边坡，如考察边坡在水平面内的形状，可将其分为凸形、凹形和直线形，边坡的空间形状对其稳定性无疑会有影响。严格来说，边坡稳定性分析属于空间问题，采用三维分析方法更符合实际情况，而工程上对边坡稳定性的评价，一般采用二维极限平衡法，这类方法对于直线形边坡计算精度尚可，但对于环形凸坡这种空间效应显著的边坡计算结果误差较大。如何分析环形凸坡的稳定性是边坡稳定性评价中需要解决的问题。本发明以目前工程中常用的二维极限平衡分析方法janbu法为基础，通过考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，将janbu法进行修正以适用于环形凸坡，使计算结果更符合实际情况。技术实现要素：针对上述问题，本发明要解决的技术问题是提供基于janbu法的环形凸坡稳定性评价方法，以解决现有二维极限平衡法在评价环形凸坡稳定性方面的不足。本发明采用的技术方案为：基于janbu法的环形凸坡稳定性评价方法，其实施过程如下：步骤一：将环形凸坡分为若干环形条块，获取每个环形条块的重力wi、滑面面积a1i。步骤二：将环形凸坡对称轴向任意一个方向拉伸形成一个剖面，每个环形条块在剖面上形成一个断面，获取剖面上每个环形条块断面面积a2i、滑面倾角θi。步骤三：通过下列公式迭代计算边坡安全系数fs。式中，ci为第i个环形条块的滑面粘聚力；为第i个环形条块的滑面内摩擦角；σ拉i为第i个环形条块土体抗拉强度；ri为第i个环形条块轴向拉力产生的抗滑力。其中，步骤三中的公式基于janbu法并考虑环形条块的轴向拉力对抗滑力的贡献，除了janbu法的基本假定外，新引入3条假定：(1)忽略环形条块间的切向力；(2)第i个环形条块轴向拉力产生的抗滑力ri，作用点位于该环形条块重心上；(3)忽略环形条块的力矩平衡。本发明的有益效果：考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，对janbu法引入新的假定进行改进，经过改进的janbu法可以用于评价环形凸坡的稳定性，且计算过程简单，为环形凸坡的稳定性评价提供了一种计算结果更为合理的方法。附图说明为了易于说明，本发明由下述的具体实施及附图作以详细描述。图1为本发明实施方式中环形条块i的结构示意图；图2为本发明实施方式中环形条块和等效环形条块的结构示意图；图3为本发明实施方式中等效环形条块janbu法受力分析图；图4为本发明实施例中环形凸坡三维模型计算参数图；图5为本发明实施例中环形凸坡断面计算参数图。具体实施方式下面将结合本发明实施例中对本发明技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。本具体实施方式基于janbu法的环形凸坡稳定性评价方法，其实施过程如下：取图4中任意一个环形条块i，如图1所示。由土体抗拉强度σ拉i得到环形条块土体轴向抗拉力fni，计算公式如公式(1)所示。fni＝σ拉ia2i(1)环形条块土体轴向抗拉力fni会产生一个环形均布线荷载qi，qi方向水平且指向对称轴，qi计算公式如公式(2)所示。公式中，ri为第i个环形条块的重力wi到对称轴的距离。把公式(1)带入(2)得到公式(3)。现建立一个截面与环形条块相同且长度与环形条块周长相等的长条形条块，将其命名为等效环形条块i，并将环形均布线荷载qi施加在等效环形条块i上如图2所示，之后将等效环形条块i作为研究对象。引入安全系数fs，等效环形条块i的抗滑力ri计算公式如公式(4)所示，ri方向水平且指向坡内。由于等效环形条块i的抗滑力ri的存在，使得环形凸坡的安全系数较长直边坡会更高一点。将ri引入到janbu法中，取整个等效环形条块i作为研究对象，将其所受的外力投影到重心所在的剖面上如图3所示。针对整个等效环形条块i，由竖向合力∑fz＝0得到公式(5)。式中，ni为等效环形条块i在滑面a1i上的法向力，ti为等效环形条块i在滑面a1i上的抗滑剪力。针对整个等效环形条块i，由水平合力∑fx＝0得到公式(6)。δpi＝pi+1-pi＝ticosθi+ri-nisinθi(6)公式中，pi为土条间法向力。公式(5)带入公式(6)得到公式(7)。δpi＝ti(cosθi+sinθitanθi)-witanθi+ri(7)由摩尔库伦强度准则并引入安全系数fs可得ti，如公式(8)所示。公式(5)带入公式(8)整理得到公式(9)。将公式(9)带入公式(7)得到公式(10)。由∑δpi＝0整理得到公式(11)～(13)。实施例：步骤一：环形凸坡三维模型计算参数图如图4所示，共分为10个环形条块，所有环形条块的滑体容重为25kn/m3。所有环形条块的滑面粘聚力ci均为33kpa，内摩擦角均为35°，所有土体抗拉强度σ拉i均为32kpa。10个环形条块的重力wi和滑面面积a1i，如表1所示。表1环形条块的重力wi和滑面面积a1i条块编号12345wi(kn)4612510047512725082475154675条块编号678910wi(kn)143325144350131925144875117925条块编号12345a1i(m2)108488410646201052条块编号678910a1i(m2)8728006726842572步骤二：环形凸坡断面计算参数图如图5所示，每个环形条块断面面积a2i、滑面倾角θi，如表2所示。表2环形条块断面面积a2i和滑面倾角θi条块编号12345a2i(m2)5.913.718.412.625条块编号678910a2i(m2)2527.32732.429.2条块编号12345θi(°)1016434343条块编号678910θi(°)4343434365步骤三：通过公式进行迭代计算安全系数fs，共迭代8次，每次迭代结果分别为1.132，1.189，1.211，1.219，1.222，1.223，1.224，1.224。通过迭代计算，得出最终安全系数fs＝1.224。本发明考虑环形凸坡各环形条块的轴向拉力对抗滑力的贡献，对janbu法引入新的假定进行改进，经过改进的janbu法可以用于评价环形凸坡的稳定性，且计算过程简单，为环形凸坡的稳定性评价提供了一种计算结果更为合理的方法。尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。当前第1页12