一种确定边坡临界滑动面的位移分析方法技术

本发明专利技术公开了一种确定边坡临界滑动面的位移分析方法，涉及边坡稳定性分析技术领域。包括步骤：选取待测坡体，并建模；使用强度折减法使边坡处于临界状态，确定离散点分布区域；在离散点分布区域设置若干条竖直线，在每条竖直线上均从头至尾分布着若干离散点，离散点分布区域的上部边界的直线上上也分布若干离散点；计算竖线上及边界的每条直线上的任意且相邻的两个离散点的连线的中点的位移变化率；分别找出竖线上及边界的直线上最大位移变化率的位置；平滑处理滑动面位置及滑动面进、出口的位置的连线。这种方法不仅保证了所判断滑动面的准确性、适用范围广且理论严谨、方法简单的要求，而且避免了人为性与经验性的弊端，更是有较强的工程实践意义。

【技术实现步骤摘要】

本专利技术涉及边坡稳定性分析

技术介绍
我国是世界上自然灾害最为严重的国家之一。有关研究表明，全球气候变暖对我国灾害风险分布和发生规律的影响将是全方位、多层次的，强台风更加活跃，暴雨洪涝灾害增多，发生流域性大洪水的可能性加大，强地震、局部强降雨引发的滑坡、山洪、泥石流等地质灾害也会加多，边坡治理的需求也将相应增大，要求也会更为精确严苛。目前针对边坡滑动面的确定，目前主要存在以下三种通过强度折减法来表征临界滑动面的方法:（1）从单元尺度近似表征滑动面，实质上此类方法仅表示剪切破坏带，并未准确刻画滑动面的位置；（2）选取计算结果中某一物理量的某条等值线作为滑动面。虽然该等值线可以具体表征滑动面的位置，但是等值线的选取具有较强的人为性和经验性；（3）根据数学力学模型进行严格的数值计算得到滑动面。该方法具有牢固的数学力学基础，但是仅仅只适用于摩尔库伦准则。上述方法均不是同时具备简单、快速且准确的有效的表征方法。
技术实现思路
本专利技术要解决的技术问题是针对上述现有技术的不足，提供一种确定边坡临界滑动面的位移分析方法，不仅保证了所判断滑动面的准确性、适用范围广且理论严谨、方法简单的要求，而且避免了人为性与经验性的弊端，更是有较强的工程实践意义。为解决上述技术问题，本专利技术所采取的技术方案是：一种确定边坡临界滑动面的位移分析方法，包括如下步骤：A、选取待测坡体，确定该坡体的尺寸，在计算机中建立该坡体的模型，并设置土体性质参数；B、使用强度折减法使边坡处于临界状态，随后确定在该坡体横断面上的离散点分布区域；C、在离散点分布区域设置若干条竖直线，在每条竖直线上均从头至尾分布着若干离散点，离散点分布区域的上部边界的直线上上也分布若干离散点；D、计算离散点分布区域内每条竖线上的任意且相邻的两个离散点的连线的中点的位移变化率，计算离散点分布区域上部边界的每条直线上的任意且相邻的两个离散点的连线的中点的位移变化率；E、分别找出离散点分布区域内每条竖线上及离散点分布区域上部边界的每条直线上最大位移变化率的位置，用于表征滑动面位置及滑动面进、出口的位置；F、平滑处理滑动面位置及滑动面进、出口的位置的连线，得到平滑的滑动面表征线。作为进一步的技术方案，所述步骤D中，采用中心差分法计算每条竖线上的任意且相邻的两个离散点及边界的每条直线上的任意且相邻的两个离散点的连线的中点计算该中点的位移变化率。作为进一步的技术方案，在同一条直线上的所述离散点均为等间距分布。作为进一步的技术方案，所述竖直线为等间距分布。采用上述技术方案所产生的有益效果在于：本专利技术的方法不仅保证了所判断滑动面的准确性、适用范围广且理论严谨方法简单的要求，而且避免了人为性与经验性的弊端，适用于任何破坏准则、编程简单、计算效率高，可以准确快速搜索边坡临界滑动面，更是有较强的工程实践意义，比如可根据本方法找到滑动面位置设计锚杆的长度、挡土墙入土深度等。附图说明图1是本专利技术的流程图；图2是本专利技术实施例的坡体模型图；图3是图2中坡体选取离散点分布区域及离散点的分布状态图；图4是图3中经过计算后得到的最大位移变化率的中点的分布图；图5是图3中计算相邻两点连线中点的位移变化率的原理图；图6是基于位移变化率确定坡体滑动面的原理图；图7是本专利技术测算结果与spencer法测算结果比较图。具体实施方式下面结合附图和具体实施方式对本专利技术作进一步详细的说明。如图1-7所示，为本专利技术一种确定边坡临界滑动面的位移分析方法，包括如下步骤：A、选取待测坡体，确定该坡体的尺寸，在计算机中建立该坡体的模型，并设置土体性质参数。B、使用强度折减法使边坡处于临界状态，强度折减法能够对复杂的地质、地貌进行计算，考虑了土体的本构关系，以及变形对应力的影响，是确定边坡临界状态的一种较准确的方法。随后确定在该坡体横断面上的离散点分布区域。该分布区域选取范围应当较大，包含较多的离散点可供测算。C、在离散点分布区域设置若干条竖直线，在每条竖直线上均从头至尾分布着若干离散点。离散点分布区域的上部边界的直线上上也分布若干离散点。为例方便计算，在同一条直线上的离散点均为等间距分布。D、离散点在临界状态处于时刻移动的状态，采用中心差分法计算离散点分布区域内每条竖线上的任意且相邻的两个离散点的连线的中点的位移变化率，采用中心差分法计算离散点分布区域上部边界的每条直线上的任意且相邻的两个离散点的连线的中点的位移变化率。滑动面附近的位移等值线最密集，滑动面上的点往往是深度方向上位移沿垂直滑面方向的位移变化率达到最大值的位置，因此，计算离散点连线中点的最大位移变化率的位置，就能够准确的表征处滑动面的位置。E、分别找出离散点分布区域内每条竖线上及离散点分布区域上部边界的每条直线上最大位移变化率的位置，用于表征滑动面位置及滑动面进、出口的位置。F、平滑处理滑动面位置及滑动面进、出口的位置的连线，得到平滑的滑动面表征线。对于中心差分法的建模及计算原理可以示例为：建立x-y总坐标系，使得滑动面处于x-y坐标系内，P点表示滑动面上任意一点，过P点建立ε-η局部坐标系，ε轴平行于P点切线方向。P点正上方及正下方分别设置离散点P1及P2，即P、P1及P2三点共线，且P是P、P1连线的中点。ε1、ε2分别表示P1、P2点沿平行于ε轴方向位移分量的大小，△h表示P1、P2点之间的距离在η轴方向的投影长度。因此，采用中心差分法计算P点沿垂直滑面方向的位移变化率（）如式（1）所示：（1）假设已知P1和P2的位移矢量在x-y坐标系中分别是和，坐标分别是和，则最终表达式如下：（2）如图5所示，滑动面把坡体分为滑动体和稳定体，稳定体的节点位移仅是单元的变形，而滑动体的节点位移由单元的变形和滑动位移组成，显然，滑动引起的节点位移要大于单元变形引起的节点位移。如图6所示，1-8号离散点位于相同的直线上，1-4号离散点位于滑动体中，5-8号离散点位于稳定体中，其中，4、5号离散点为靠近滑动面且分别位于滑动体和稳定体中的点。根据式（1）和式（2）的定义，可计算任意相邻两点连线中点的位移变化率，显然，4、5号连线中点的位移变化率要大于其它任意连续两点连线中点的位移变化率。因此，可认为4、5号连线中点ο是该直线与滑动面的交点。基于以上认识，若沿水平方向布置一系列竖直线，在每条竖直线上布置若干离散点，计算每条竖直线上任意相邻两个离散点本文档来自技高网...

【技术保护点】
一种确定边坡临界滑动面的位移分析方法，其特征在于包括如下步骤：A、选取待测坡体，确定该坡体的尺寸，在计算机中建立该坡体的模型，并设置土体性质参数；B、使用强度折减法使边坡处于临界状态，随后确定在该坡体横断面上的离散点分布区域；C、在离散点分布区域设置若干条竖直线，在每条竖直线上均从头至尾分布着若干离散点，离散点分布区域的上部边界的直线上上也分布若干离散点；D、计算离散点分布区域内每条竖线上的任意且相邻的两个离散点的连线的中点的位移变化率，计算离散点分布区域上部边界的每条直线上的任意且相邻的两个离散点的连线的中点的位移变化率；E、分别找出离散点分布区域内每条竖线上及离散点分布区域上部边界的每条直线上最大位移变化率的位置，用于表征滑动面位置及滑动面进、出口的位置；F、平滑处理滑动面位置及滑动面进、出口的位置的连线，得到平滑的滑动面表征线。

【技术特征摘要】
1.一种确定边坡临界滑动面的位移分析方法，其特征在于包括如下步骤：
A、选取待测坡体，确定该坡体的尺寸，在计算机中建立该坡体的模型，并设置土体性质
参数；
B、使用强度折减法使边坡处于临界状态，随后确定在该坡体横断面上的离散点分布区
域；
C、在离散点分布区域设置若干条竖直线，在每条竖直线上均从头至尾分布着若干离散
点，离散点分布区域的上部边界的直线上上也分布若干离散点；
D、计算离散点分布区域内每条竖线上的任意且相邻的两个离散点的连线的中点的位
移变化率，计算离散点分布区域上部边界的每条直线上的任意且相邻的两个离散点的连线
的中点的位移变化率；
E、分别找出离散点分布区域内每条竖线上及离散点分布区域上部边界的每条直线上
最大位移变化率的位置，用于表征滑动面位置及滑动面进、出口...

【专利技术属性】
技术研发人员：袁维，胡叶江，胡小莉，马斌，
申请(专利权)人：石家庄铁道大学，
类型：发明
国别省市：河北;13