一种基于动态强度折减DDA法的边坡稳定性分析系统的制作方法

本发明属于边坡稳定性分析领域，尤其涉及一种基于动态强度折减dda法的边坡稳定性分析系统。

背景技术：

目前，边坡的稳定性问题是我国水利水电、公路、铁路和矿产资源开发等建设工程中常见的岩土工程问题。自然滑坡、泥石流以及人类工程活动等引起的边坡灾害对我国经济建设和人民生命财产带来了巨大损失，因此边坡工程在各类工程建设中的地位是十分重要的，正确的评价边坡的稳定性，防患于未然，能确保生产建设与人民财产安全。

边坡稳定性分析中有两种常用的计算方法：极限平衡法和有限元法，但只能假设边坡为连续介质，而非连续性是岩体固有的属性，所以在模拟块体大位移、大变形方面独具优势的非连续变形分析(dda)，在边坡工程领域也得到迅速发展。过去研究中的强度折减法大都是对岩体统一折减，在实际情况下，主要是结构面的力学性质决定了岩体的力学性质，且岩体结构面破坏时并不是同时损伤同等程度，所以对边坡破坏过程中岩体结构面损伤不同程度进行动态表征具有十分重要的意义。

技术实现要素：

为解决上述技术问题，本发明提供了一种基于动态强度折减dda法的边坡稳定性分析系统，对不同破坏程度下岩体结构面动态折减不同的强度，以解决数值分析中岩体结构面损伤的程度不一性问题，使得边坡稳定性分析更加准确，以便工程人员及时采取措施治理危险边坡。

具体技术方案如下：

一种基于动态强度折减dda法的边坡稳定性分析系统，包括以下步骤：

(1)输入几何、物理、计算参数，确定边坡的几何尺寸、岩石容重、弹性模量、泊松比等；

(2)块体识别并进行接触检测，结构面切割岩体形成块体，对块体进行识别，并进行块体间的接触检测；

(3)形成子矩阵，不包括接触子矩阵，其中弹性子矩阵为：

式中，e、v为块体单元的弹性模量和泊松比；

初始应力子矩阵为：

式中，为块体初始常应力；

点荷载子矩阵为：

式中，(fx，fy)为点荷载力；

体积力荷载子矩阵为：

式中，(fx，fy)为体积力荷载，s为块体单元的面积；

(4)加减接触弹簧形成接触子矩阵，其中包括法向弹簧子矩阵与切向弹簧子矩阵；在每一时步内，都要重新确定弹簧的施加与否及弹簧的位置，需要反复生成求解总刚矩阵；刚性弹簧的施加与去除过程称之为开—合迭代；如果接触在法向弹簧上出现拉力，在刚性弹簧去除后两边将分离；如果顶点嵌入接触的另一边，要加入刚性弹簧，将其拉回接触面上；接触有三种状态：张开、滑动和锁定；根据相应的接触状态变化来对块体的位移变化进行更新。当接触状态变化分别为张开-张开、张开-滑动、张开-锁定、滑动-张开、滑动-滑动、滑动-锁定、锁定-张开、锁定-滑动、锁定-锁定时，对应位移变化为n＞0、n＜0，n＜0，n＞0、n＜0，n＜0，n＞0、n＜0，n＜0，其中：n是法向位移，n＞0为张开，是剪切位移矢量，t是剪切位移，||表示两个矢量方向相同；

(5)求解位移矢量平衡方程，每个块体的位移矢量包含六个变量：

d＝(u0,v0,r0,εx,εy,γxy)^t

其中(u0,v0)是块体内特定点(x0,y0)的刚性位移，r0是块体绕特定点(x0,y0)的旋转角度，εx、εy、γxy分别为该块体的正应变与切应变；块体内任意点(x,y)的位移(u,v)为：

基于最小势能原理构建的系统整体平衡方程为：

式中：系数矩阵中kij为6×6子矩阵；kii由块体单元的材料属性和几何参数决定；kij(i≠j)则由块体i和块体j间的接触条件而决定；[di]和[fi]是6×1子矩阵，di为块体i的变形变量(d1i,d2i,d3i,d4i,d5i,d6i)；fi是块体i上分配给6个变形变量的荷载；

(6)检查收敛性，系统根据计算结果是否收敛来决定如何运行下一步程序，若收敛，则根据接触判定点计算块体相对滑动位移并与位移阈值0.001m作对比，根据对比结果动态变化结构面剪切强度参数c和具体公式如下：

式中：c'和分别指折减后的结构面抗剪强度参数黏聚力和内摩擦角；fsr是折减系数，处于临界状态时的折减系数fsr即边坡滑落面安全系数fs，n指代块体间结构面相对位移满足阈值条件的次数；若不收敛，则返回至(4)；

(7)位移突变，则判定边坡失稳，并记录安全系数并检索出滑落面位置，反之，则返回(2)。

所述边坡滑落面安全系数fs可根据边坡滑落面安全系数大小判定边坡稳定性，安全系数越大，则边坡稳定性越好。

步骤(7)所述滑落面位置可根据系统检索出的滑落面位置，来对将来可能发生失稳的边坡区域进行提前开挖，以免边坡自行失稳，造成不必要的损失。

与现有技术相比，本发明具有如下有益技术效果：

现有的边坡稳定性分析系统一般采用静态均一性分析，而边坡失稳过程是一个动态变化过程，且各岩体结构面强度衰减程度是不一的，无法对而边坡失稳过程是一个动态变化过程，所以以往得到的边坡滑动面和相应的安全系数的准确性有待提高。本发明系统借鉴dda计算位移的优势，考虑块体间相对位移，提出动态强度折减dda法，对不同破坏程度下岩体结构面动态折减不同的强度，解决了边坡破坏过程中结构面损伤的程度不一性问题，并计算出与实际情况更相符安全系数与滑落面位置，可根据检索出的滑落面位置提前剥离边坡危险区域岩体，最大程度地降低损失。

附图说明

图1为基于动态强度折减dda法的边坡稳定性分析系统程序步骤流程图；

图2为实施例边坡原始状态示意图；

图3为经本发明系统计算后的实施例边坡滑落面位置示意图。

具体实施方式

下面结合附图对本发明进行详细说明，但本发明的保护范围不受实施例所限。

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。

本发明的一种基于动态强度折减dda法的边坡稳定性分析系统，包括以下内容：

(1)输入几何、物理、计算参数，确定边坡的几何尺寸、岩石容重、弹性模量、泊松比等，具体数值分别为：坡高为210m，坡底为197m，岩石容重为21000kn/m³，，弹性模量为1700mpa，泊松比为0.35，黏聚力为40kpa，内摩擦角为26°；

(2)块体识别并进行接触检测，结构面切割岩体形成块体，对块体进行识别，并进行块体间的接触检测；

(3)形成子矩阵，不包括接触子矩阵，其中弹性子矩阵为：

式中，e、v为块体单元的弹性模量和泊松比；

初始应力子矩阵为：

式中，为块体初始常应力；

点荷载子矩阵为：

式中，(fx，fy)为点荷载力；

体积力荷载子矩阵为：

式中，(fx，fy)为体积力荷载，s为块体单元的面积。

(4)加减接触弹簧形成接触子矩阵，其中包括法向弹簧子矩阵与切向弹簧子矩阵。在每一时步内，都要重新确定弹簧的施加与否及弹簧的位置，需要反复生成求解总刚矩阵。刚性弹簧的施加与去除过程称之为开—合迭代。如果接触在法向弹簧上出现拉力，在刚性弹簧去除后两边将分离；如果顶点嵌入接触的另一边，要加入刚性弹簧，将其拉回接触面上。接触有三种状态：张开、滑动和锁定。模式变化的判断准则为：

其中：n是法向位移，n＞0为张开，是剪切位移矢量，t是剪切位移，||表示两个矢量方向相同。

(5)求解位移矢量平衡方程，每个块体的位移矢量包含六个变量：

d＝(u0,v0,r0,εx,εy,γxy)^t

其中(u0,v0)是块体内特定点(x0,y0)的刚性位移，r0是块体绕特定点(x0,y0)的旋转角度，εx、εy、γxy分别为该块体的正应变与切应变。块体内任意点(x,y)的位移(u,v)为：

基于最小势能原理构建的系统整体平衡方程为：

式中：系数矩阵中kij为6×6子矩阵；kii由块体单元的材料属性和几何参数决定；kij(i≠j)则由块体i和块体j间的接触条件而决定；[di]和[fi]是6×1子矩阵，di为块体i的变形变量(d1i,d2i,d3i,d4i,d5i,d6i)；fi是块体i上分配给6个变形变量的荷载。

(6)检查收敛性，系统根据计算结果是否收敛来决定如何运行下一步程序，若收敛，则根据接触判定点计算块体相对滑动位移并与位移阈值0.001m作对比，根据对比结果动态变化结构面剪切强度参数c和具体公式如下：

式中：c'和分别指折减后的结构面抗剪强度参数黏聚力和内摩擦角；fsr是折减系数，处于临界状态时的折减系数fsr即边坡滑落面安全系数fs，n指代块体间结构面相对位移满足阈值条件的次数；若不收敛，则返回至(4)；

(7)位移突变，则判定边坡失稳，并记录安全系数并检索出滑落面位置，反之，则返回(2)。

根据相关设计规范规定，边坡在正常工况下的安全系数大于1.200～1.300时，可认为其满足稳定性要求。系统测得该边坡滑落面安全系数为1.189，不满足稳定性要求，需尽快对其采取措施进行治理。

根据图3所表征的滑落面位置，对预失稳的边坡区域进行提前开挖，避免边坡自行失稳，造成不必要的损失。

以上所述仅为本发明的较佳实施例，并不用以限制本发明的思想，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。