边坡失稳后果量化分析方法技术

技术编号:20390053 阅读:86 留言:0更新日期:2019-02-20 02:55
本发明专利技术公开一种边坡失稳后果量化分析方法,首先通过极限平衡方法计算得到临界滑动面,其次利用光滑粒子流体动力学方法计算得到滑坡位移场,进而通过对比分析确定有效位移阈值,接着利用一系列边坡失稳工况拟合边坡失稳最终滑动面面积与其安全系数之间的关系式,最终利用拟合的关系式计算安全系数小于1情况下的滑动面积,并结合滑动面面积直方图来量化边坡失稳后果,克服传统方式边坡失稳后果量化效果差等缺陷,能够合理、有效量化边坡失稳后果,对于滑坡风险的防治具有重要意义。

【技术实现步骤摘要】
边坡失稳后果量化分析方法
本专利技术属于边坡灾害评价与防治
,具体涉及一种边坡失稳后果量化分析新方法,特别是一种基于光滑粒子流体动力学方法确定边坡失稳最终滑动面,进而与极限平衡方法计算所得最小安全系数建立关联,最终量化边坡失稳的后果。
技术介绍
在边坡灾害防治领域,一般采用边坡稳定分析方法,譬如极限平衡方法和有限元强度折减法进行最小安全系数的搜索,并利用其对应的临界滑动面来评估边坡失稳的后果。然而,对于濒临失稳或者已经失稳的边坡而言,极限平衡方法和有限元强度折减方法所得到的临界滑动面仅仅是滑坡发生初始时刻的起滑面,该起滑面与滑坡发生导致的最终滑动面存在一定差别,从而给边坡失稳后果量化带来误差,因此目前亟需一种能够合理、有效量化边坡失稳后果的方法。
技术实现思路
本专利技术所要解决的技术问题在于针对现有技术采用常规极限平衡方法或有限元强度折减方法所得到的临界滑动面仅仅是滑坡发生初始时刻的起滑面,与最终滑动面存在一定的差别,导致边坡失稳后果量化效果差等缺陷,提供一种用于边坡失稳后果量化的新方法。本专利技术是采用以下的技术方案实现的:一种边坡失稳后果量化分析方法,包括以下步骤:S1、对拟研究边坡的最小安全系数Fmin及其滑动面面积A进行分析计算;S2、将边坡原始力学参数除以最小安全系数Fmin得到折减参数,基于得到的折减参数,利用光滑粒子流体动力学方法计算与模拟得到边坡滑动的位移场,所述边坡原始参数包括边坡材料的力学参数粘聚力c和内摩擦角折减后的粘聚力参数cn=c/Fmin,内摩擦角S3、设定位移阈值δd,逐一判断每个粒子的位移值是否大于δd,若大于δd则该粒子视为滑动粒子,否则视为静止粒子,从位移场中选择位移值高于δd的粒子(即滑动粒子)并汇总其面积,将汇总得到的粒子面积记为B,若B不等于A,则变换不同的δd,重新计算滑动粒子的面积B,直至汇总的粒子面积B等于原始力学参数下边坡最小安全系数对应的滑动面面积A,此时的位移阈值δd即为有效的位移阈值δ0,即最终确定了有效位移阈值δ0;S4、根据粘聚力参数cn和内摩擦角设计m组致滑参数组合ci和i=1,2,……,m,首先用极限平衡方法计算其安全系数Fi,然后利用光滑粒子流体动力学方法计算相应的位移场,并根据有效位移阈值δ0计算每种工况下边坡失稳最终滑动面及其滑动面积Bi;S5、在EXCEL中以Fi为横轴,Bi为纵轴,绘制Fi~Bi之间的散点图,并采用拟合曲线拟合Fi~Bi之间的变化趋势,所得到的拟合函数为B=f(F);S6、生成符合假定概型的蒙特卡罗样本,针对每一样本,利用极限平衡方法计算其安全系数F,若安全系数F大于或等于1,则表明边坡未失稳;若安全系数F小于1,则利用函数曲线B=f(F)计算相应的滑动面积,最终通过滑动面积直方图量化边坡失稳后果。进一步的,所述步骤S4中,m的取值为6-12。进一步的,所述步骤S4中,ci=cn-(cn-cs)×i/m,其中,cs和为根据极限平衡法中Fmin=0.1确定的粘聚力值和内摩擦角值;基于ci,采用光滑粒子流体动力学方法,并利用有效位移阈值δ0确定滑动粒子面积,记为Bi;基于ci,采用极限平衡方法计算其安全系数,记为Fi。进一步的,所述步骤S6中,在EXCEL中利用NORMDIST函数生成正态分布变量,进而利用相关转换公式得到符合假定概型的蒙特卡罗样本值cj和j=1,2,……N,N为样本数目;基于cj和采用极限平衡方法计算其安全系数Fj,判断Fj是否小于1,若是,则利用B=f(Fj)计算边坡失稳后的滑动面积Bj;若否,则表明边坡未失稳,不计算滑动面积;当N个样本值全部计算完毕后,统计边坡失稳滑动面积的直方图,进行边坡失稳后果量化。与现有技术相比,本专利技术的优点和积极效果在于:本方案所提出的边坡失稳后果量化分析方法,首先通过极限平衡方法计算得到临界滑动面,其次利用光滑粒子流体动力学方法计算得到滑坡位移场,进而通过对比分析确定有效位移阈值,接着利用一系列边坡失稳工况拟合边坡失稳最终滑动面面积与其安全系数之间的关系式,最终利用拟合的关系式计算安全系数小于1情况下的滑动面积,并结合滑动面面积直方图来量化边坡失稳后果,不仅能够合理、有效量化边坡失稳后果,而且还能够为边坡支护与加固提供决策支持,根据滑动面积的平均值与标准差能够有效地进行灾害防治成本预算,对于滑坡风险的防治具有重要的意义。附图说明图1为本专利技术实施例所述边坡失稳后果量化分析方法流程框图;图2为本专利技术实施例中黏性土边坡示意图;图3为本专利技术实施例中通过光滑粒子动力学方法计算所得的位移云图;图4为本专利技术实施例有效位移阈值δ0确定示意图;图5为本专利技术实施例Fi~Bi拟合示意图;图6为传统方法边坡失稳后果量化示意图;图7为采用本专利技术实施例所述方法的边坡失稳后果量化示意图。具体实施方式为了能够更加清楚地理解本专利技术的上述目的、特征和优点,下面结合附图及实施例对本专利技术做进一步说明。需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。如图1所示,本实施例公开一种边坡失稳后果量化分析方法,包括以下步骤:S1、首先对待研究边坡进行稳定性验算,利用极限平衡方法计算边坡最小安全系数Fmin和相对应的滑动面S,并计算得到其对应的滑动面面积A;S2、将边坡原始力学参数除以最小安全系数Fmin得到折减参数,基于得到的折减参数,利用光滑粒子流体动力学方法计算与模拟得到边坡滑动的位移场,所述边坡原始参数包括边坡材料的力学参数粘聚力c和内摩擦角折减后的粘聚力参数cn=c/Fmin,内摩擦角S3、设定位移阈值δd(比如0.01m),逐一判断每个粒子的位移值是否大于δd,若大于δd则该粒子视为滑动粒子,否则视为静止粒子,从位移场中选择位移值高于δd的粒子(即滑动粒子)并汇总其面积,将汇总得到的粒子面积记为B,若B不等于A,则变换不同的δd,重新计算滑动粒子的面积B,直至汇总的粒子面积B等于原始力学参数下边坡最小安全系数对应的滑动面面积A,此时的位移阈值δd即为有效的位移阈值δ0,即最终确定了有效位移阈值δ0;S4、根据粘聚力参数cn和内摩擦角设计m组致滑参数组合ci和i=1,2,……,m,m的取值为6-12,ci=cn-(cn-cs)×i/m,cs和是根据极限平衡方法中Fmin=0.1所确定的粘聚力值和内摩擦角值;基于ci,采用光滑粒子流体动力学方法,并利用有效位移阈值δ0确定滑动粒子面积,记为Bi;基于ci,采用极限平衡方法计算其安全系数,记为Fi;S5、在EXCEL中以Fi为横轴,Bi为纵轴,绘制Fi~Bi之间的散点图,并采用拟合曲线拟合Fi~Bi之间的变化趋势,所得到的拟合函数为B=f(F);S6、在EXCEL中利用NORMDIST函数生成正态分布变量,进而利用相关转换公式生成符合假定概型的蒙特卡罗样本值cj和j=1,2,……N,N为样本数目;针对每一样本,基于cj和采用极限平衡方法计算其安全系数Fj,判断Fj是否小于1,若是,则利用B=f(Fj)计算边坡失稳后的滑动面积Bj;若否,则表明边坡未失稳,不计算滑动面积;当N个样本值全部计算完毕后,统计边坡失稳滑动面积的直方图,进行边坡失稳后果量化。如图2所示,某粘性土边坡,坡高2m,土的粘聚力c为8.0kpa,内摩擦角为0°,按照本文档来自技高网...

【技术保护点】
1.边坡失稳后果量化分析方法,其特征在于,包括以下步骤:S1、对拟研究边坡的最小安全系数Fmin及其滑动面面积A进行分析计算;S2、将边坡原始力学参数除以最小安全系数Fmin得到折减参数,基于得到的折减参数,利用光滑粒子流体动力学方法计算与模拟得到边坡滑动的位移场,所述边坡原始参数包括边坡材料的力学参数粘聚力c和内摩擦角

【技术特征摘要】
1.边坡失稳后果量化分析方法,其特征在于,包括以下步骤:S1、对拟研究边坡的最小安全系数Fmin及其滑动面面积A进行分析计算;S2、将边坡原始力学参数除以最小安全系数Fmin得到折减参数,基于得到的折减参数,利用光滑粒子流体动力学方法计算与模拟得到边坡滑动的位移场,所述边坡原始参数包括边坡材料的力学参数粘聚力c和内摩擦角折减后的粘聚力参数cn=c/Fmin,内摩擦角S3、设定位移阈值δd,从位移场中选择位移值高于δd的粒子并汇总其面积,将汇总得到的粒子面积记为B,若B不等于A,则变换不同的δd,直至汇总的粒子面积B等于原始力学参数下边坡最小安全系数对应的滑动面面积A,最终确定有效位移阈值δ0;S4、根据粘聚力参数cn和内摩擦角设计m组致滑参数组合ci和首先用极限平衡方法计算其安全系数Fi,然后利用光滑粒子流体动力学方法计算相应的位移场,并根据有效位移阈值δ0计算每种工况下边坡失稳最终滑动面及其滑动面积Bi;S5、绘制Fi~Bi之间的散点图,并采用拟合曲线拟合Fi~Bi之间的变化趋势,所得到的拟合函数为B=f(F);S6、生成符合假定概型的蒙特卡罗样本,针对每一样本,利用极限平衡方法计算其安全系数...

【专利技术属性】
技术研发人员:李亮翟明褚雪松袁长丰洪勇
申请(专利权)人:青岛理工大学
类型:发明
国别省市:山东,37

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1