[发明专利]一种受陡缓倾角控制的碎裂松动岩体边坡稳定性分析方法

权利要求书

1.一种受陡缓倾角控制的碎裂松动岩体边坡稳定性分析方法，其特征在于：包括以下步骤：

步骤一：收集研究区岩体结构发育特征，通过统计分析寻找出陡缓倾角组合成长大裂隙控制结构面，及各类结构面的统计分布规律；结合现场工程地质勘察结果，建立坝肩边坡的二维数值计算模型；

步骤二：随机结构面特征参数获取及模拟过程，通过调查获取确定性结构面的倾向、倾角的信息，再对确定性结构面通过内聚理论进行分组，确定每组的统计参数，代入分析计算公式，解析出概率密度分布函数统计参数特征值；再根据结构面分布统计体的密度，模拟出结构面数量N，运用概率分布密度函数模拟出随机结构面的统计参数特征值，输入UDEC程序中，生成随机分布结构面；

步骤三：在UDEC程序中将调查获得的陡缓倾角组合长大裂隙数据信息导入离散元模型，输入岩体力学参数和结构面强度参数，计算碎裂边坡稳定性；

步骤四：在UDEC程序中模拟分析天然条件下边坡变形特征，包括位移场和塑性区变化情况，同时分析边坡失稳破坏模式。

2.根据权利要求1所述的受陡缓倾角控制的碎裂松动岩体边坡稳定性分析方法，其特征在于：所述步骤二具体如下：

(1)岩体结构面模拟过程中，采用Beacher圆盘模型，通过圆盘中心的空间坐标(x₀，y₀，z₀)来确定结构面的位置，及圆盘与直角坐标系的夹角α，β来确定结构面倾向倾角，圆盘半径r用于确定结构面的大小，岩体结构面模型按照如下表达式计算：

A(x-x₀)+B(y-y₀)+c(z-z₀)＝0 (式1)

(x-x₀)²+(y-y₀)²+(z-z₀)²r (式2)

式中：A＝sinαsinβ；B＝sinαcosβ；C＝cosα；

(2)随机结构面模拟过程中，首先选取确定性结构面统计参数特征值中的最大和最小值，生成平均密度分布函数，然后加入Monte Carlo随机数，生成随机结构面空间产状，具体如下：

令F(θ)＝μ，取反函数，解方程得到：

θ′＝(θ_max-θ_min)μ+θ_min，0≤θ≤π (式5)

式中：θ′、为随机结构面空间产状特征值，μ、v为(0，1)范围内Monte Carlo随机数。

3.根据权利要求1所述的受陡缓倾角控制的碎裂松动岩体边坡稳定性分析方法，其特征在于：所述步骤三具体如下：

(1)导入陡-缓倾角组合长大裂隙坐标信息(x_i，y_i)及倾向、倾角，将陡-缓倾角组合长大裂隙结构面设置弱接触；

(2)Crack添加节理，将随机分布模型所得随机结构面空间产状特征值θ′、等分布特征参数输入；

(3)输入岩体力学参数γ，C，及长大裂隙、结构面强度参数C，Kn，Ks。

(4)设定块体位移监测，使用“cell-space detection”模式跟踪位移及下落块体的潜在接触情况。

4.根据权利要求1所述的受陡缓倾角控制的碎裂松动岩体边坡稳定性分析方法，其特征在于：所述步骤四具体如下：

(1)模型计算采用弹塑性模型，屈服条件遵循莫尔—库伦准则；力学边界通过固定边界的形式进行约束，即采用约束速度的方式，分别将模型底部竖直方向(Y方向，yvel＝0)、两侧水平方向(X方向，xvel＝0)的速度设置为零；

(2)utility访问不同变量，监测即将运动的块体垂直和水平位移；设置重力9.81；run-solve自动检测平衡；Slove工具实现岩体稳定性评价计算；

(3)计算达到稳定状态后，得到岩体周边的位移及塑性区变化情况，通过不同时间步位移云图及塑性区云图分析边坡失稳模式。