一种基于下游构筑物受损程度的边坡安全度评价方法技术

本发明专利技术属于边坡安全评价与风险量化领域，尤其涉及一种基于下游构筑物受损程度的边坡安全度评价方法，包括如下步骤，步骤1：根据边坡材料力学参数的不确定性，假定生成一系列参数组合值

A Slope Safety Evaluation Method Based on Damage Degree of Downstream Structures

The present invention belongs to the field of slope safety evaluation and risk quantification, and in particular relates to a slope safety evaluation method based on the damage degree of downstream structures. It includes the following steps: step 1: according to the uncertainty of mechanical parameters of slope materials, a series of parameters combination values are assumed to be generated.

【技术实现步骤摘要】
一种基于下游构筑物受损程度的边坡安全度评价方法
本专利技术属于边坡安全评价与风险量化领域，尤其涉及一种基于下游构筑物受损程度的边坡安全度评价方法。
技术介绍
滑坡是一种非常常见的地质灾害之一，尤其在夏季，我国每年都会遭受降雨引发的滑坡灾害。发生滑坡之后，滑动的土体和岩体对边坡下游的建筑物或构筑物产生冲击作用，造成不同程度的损害。发生在我国比较大型的滑坡有，千将坪古滑坡，重庆武隆滑坡和深圳光明新区余泥渣土滑坡，给我国造成了不同程度的损失。因此，考虑边坡下游构筑物受损程度的边坡安全评价就显得尤为重要。在边坡安全评价领域，一般采用极限平衡方法和有限元强度折减法进行最小安全系数的搜索，利用最小安全系数指标来评价边坡安全程度。然而，对于下游存在构筑物或建筑物的边坡而言，在评价其安全程度时，因为没有结合滑动体对建筑物或构筑物的损害程度以及边坡构筑材料的不确定性，而仅依靠安全系数，故不能直观地进行评价，不利于滑坡风险的防治。所以目前在评价下游存在构筑物或者建筑物的边坡安全程度时，在评价结果的合理和有效性上亟待提高。
技术实现思路
根据以上现有技术的不足，本专利技术提供了一种基于下游构筑物受损程度的边坡安全度评价方法，其能有效解决现有技术中在评价下游存在构筑物或建筑物的边坡的安全度时不够有效、合理、直观，不利于滑坡风险防治的问题。本专利技术解决的技术问题采用的技术方案为：一种基于下游构筑物受损程度的边坡安全度评价方法，一种基于下游构筑物受损程度的边坡安全度评价方法，包括如下步骤：步骤1：针对所研究的边坡，因边坡材料力学参数的不确定性，基于材料力学参数的平均值和标准差σ，并根据材料力学参数3σ变化法则，确定力学参数的上下限值，并在此限值内按照均匀分布假定生成一系列参数组合值步骤2：基于第i组参数组合ci，i＝1,2...,N，计算其边坡的安全系数Fsi，并采用光滑粒子流体动力学方法计算边坡滑动体的运动轨迹，并记录每个粒子所代表土体单元的最终位置；步骤3：根据下游边坡构筑物的位置和步骤2中得到的边坡滑动体的运动轨迹中最远点粒子位置的相对关系，计算两者之间的距离di，i＝1,2...,N；步骤4：重复步骤2和步骤3，得到一系列距离值和一系列的安全系数步骤5：计算一系列距离值的平均值，记为dm，计算一系列安全系数的平均值，记为Fsm；步骤6：利用距离值的平均值dm和安全系数的平均值Fsm的乘积对下游存在构筑物边坡的安全度进行评价。通过光滑粒子流体动力学方法确定滑动体轨迹，能够确定构筑物受损程度，并综合考虑边坡构筑材料的不确定性以及安全系数，最终评价边坡安全程度，相对于现有技术而言，能更加直观的评价下游存在构筑物或建筑物的边坡的安全程度，有利于滑坡风险的防治。其中，优选方式为：在步骤2中，利用极限平衡方法中的简化毕晓普法计算边坡的安全系数Fsi，i＝1,2...,N。步骤3的具体实现为，下游边坡构筑物的位置坐标记为(xg,yg)，选取运动轨迹中最远点粒子，并确定其滑动后的位置坐标(xm,ym)；计算下游边坡构筑物和运动轨迹中最远点粒子两者之间的距离值di，其中，若运动轨迹中最远点粒子超过构筑物，则距离值di取负值，即否则，距离值di保持不变。在步骤5中，距离值的平均值dm具体为安全系数的平均值Fsm具体为本专利技术具有以下有益效果：本专利技术首先通过光滑粒子流体动力学方法计算滑坡滑动体的运动轨迹；其次，判断滑动体最远点至构筑物的距离，并根据超过和未超过构筑物来判断距离的正负，并用该正负距离值初步评价边坡的安全度；再次，考虑边坡材料参数的不确定性，变换不同的材料参数重复以上两步计算，并保存每次计算时滑动体运动最远点至构筑物的正负距离值；最后利用该正负距离值的均值和安全系数来综合评价边坡安全度，相比于现有技术，评价下游存在构筑物或者建筑物的边坡安全度时更加有效、合理、直观，有利于滑坡风险的防治。附图说明图1是本专利技术所述基于下游构筑物受损程度的边坡安全度评价方法的流程图；图2是本专利技术所提供实施例的粘性土边坡与构筑物位置示意图；图3是本专利技术所提供实施例的光滑粒子流体动力学方法计算所得运动轨迹图。具体实施方式下面结合附图对本专利技术做进一步描述。实施例一：如图1～图3所示，本专利技术所述的一种基于下游构筑物受损程度的边坡安全度评价方法，包括一种基于下游构筑物受损程度的边坡安全度评价方法，其特征在于：包括如下步骤，步骤S1：针对所研究的边坡，因边坡材料力学参数的不确定性，基于材料力学参数的平均值和标准差σ，并根据材料力学参数3σ变化法则，确定力学参数的上下限值，并在此限值内按照均匀分布假定生成一系列参数组合值步骤S2：基于第i组参数组合ci，i＝1,2...,N，计算其边坡的安全系数Fsi，并采用光滑粒子流体动力学方法计算边坡滑动体的运动轨迹，并记录每个粒子所代表土体单元的最终位置；步骤S3：根据下游边坡构筑物的位置和步骤2中得到的边坡滑动体的运动轨迹中最远点粒子位置的相对关系，计算两者之间的距离di，i＝1,2...,N；步骤S4：重复步骤2和步骤3，得到一系列距离值和一系列的安全系数步骤S5：计算一系列距离值的平均值，记为dm，计算一系列安全系数的平均值，记为Fsm；步骤S6：利用距离值的平均值dm和安全系数的平均值Fsm的乘积对下游存在构筑物边坡的安全度进行评价。通过光滑粒子流体动力学方法确定滑动体轨迹，能够确定构筑物受损程度，并综合考虑边坡构筑材料的不确定性以及安全系数，最终评价边坡安全程度，相对于现有技术而言，能更加直观的评价下游存在构筑物或建筑物的边坡的安全程度，有利于滑坡风险的防治。在步骤S2中，利用极限平衡方法中的简化毕晓普法计算边坡的安全系数Fsi，i＝1,2...,N。步骤S3的具体实现为，下游边坡构筑物的位置坐标记为(xg,yg)，选取边坡滑动体的运动轨迹中最远点粒子，并确定其滑动后的位置坐标(xm,ym)；计算下游边坡构筑物和运动轨迹中最远点粒子两者之间的距离值di，其中，若运动轨迹中最远点粒子超过构筑物，则距离值di取负值，即否则，距离值di保持不变。在步骤S5中，距离值的平均值dm具体为安全系数的平均值Fsm具体为以下是本专利技术所述方法在具体应用场景中的应用。步骤S1，如图2所示，某粘性土边坡，坡高4m，土的重度γ＝20kN/m3，内摩擦角土的粘聚力c均值为15.0kPa，考虑粘聚力c的不确定性，其标准差σ为1.67kPa，根据3σ法则，在10kPa到20kPa区间内，均匀生成N＝11个c值，即ci＝i+9.0,i＝1,2,....,11；步骤S2，在第一组参数c1下，即γ＝20kN/m3，c1＝10kPa，按照极限平衡方法中的简化毕晓普法计算其安全系数Fs1为0.73；将图2所示边坡离散为9020个直径为0.1m的圆形粒子，将以上三个参数输入光滑粒子流体动力学方法程序中，得到滑动体运动轨迹图，如图3所示；步骤S3，从滑动体运动轨迹图中可以读出，滑动体运动轨迹上最远点粒子的位置大致位于xm＝5，ym＝2.0；水塔坐标为xg＝10.0，yg＝2.0，因此滑动体运动轨迹上的最远点与水塔之间的距离值d1＝5.0，又因为该滑动体运动轨迹上的最远点粒子超过了构筑物(即水塔)，因此距离值取负值，即d1＝-5.0。本专利技术中，利用距离值d1的值来本文档来自技高网...

【技术保护点】
1.一种基于下游构筑物受损程度的边坡安全度评价方法，其特征在于：包括如下步骤，步骤1：针对所研究的边坡，因边坡材料力学参数的不确定性，基于材料力学参数的平均值和标准差σ，并根据材料力学参数3σ变化法则，确定力学参数的上下限值，并在此限值内按照均匀分布假定生成一系列参数组合值

【技术特征摘要】
1.一种基于下游构筑物受损程度的边坡安全度评价方法，其特征在于：包括如下步骤，步骤1：针对所研究的边坡，因边坡材料力学参数的不确定性，基于材料力学参数的平均值和标准差σ，并根据材料力学参数3σ变化法则，确定力学参数的上下限值，并在此限值内按照均匀分布假定生成一系列参数组合值步骤2：基于第i组参数组合ci，i＝1,2...,N，计算其边坡的安全系数Fsi，并采用光滑粒子流体动力学方法计算边坡滑动体的运动轨迹，并记录每个粒子所代表土体单元的最终位置；步骤3：根据下游边坡构筑物的位置和步骤2中得到的边坡滑动体的运动轨迹中最远点粒子位置的相对关系，计算两者之间的距离di，i＝1,2...,N；步骤4：重复步骤2和步骤3，得到一系列距离值和一系列的安全系数步骤5：计算一系列距离值的平均值，记为dm，计算一系列安全系数的平均值，记为Fsm；步骤6：利用距离值的平均...

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