CN108763697A - 一种岩质边坡中岩块失效概率的计算方法 - Google Patents
一种岩质边坡中岩块失效概率的计算方法 Download PDFInfo
- Publication number
- CN108763697A CN108763697A CN201810478606.4A CN201810478606A CN108763697A CN 108763697 A CN108763697 A CN 108763697A CN 201810478606 A CN201810478606 A CN 201810478606A CN 108763697 A CN108763697 A CN 108763697A
- Authority
- CN
- China
- Prior art keywords
- side slope
- rock side
- block movement
- rigid block
- sillar
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 126
- 238000000205 computational method Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 230000014509 gene expression Effects 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000009977 dual effect Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000010008 shearing Methods 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/23—Dune restoration or creation; Cliff stabilisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
本发明公开了一种岩质边坡中岩块失效概率的计算方法,本发明方法为:步骤1、拟定岩质边坡的计算参数;步骤2、建立岩质边坡的极限状态函数;步骤3、使用刚性块体单元离散岩质边坡;步骤4、建立求解岩块失效概率的上限法线性规划模型;步骤5、使用蒙特卡洛法求解上限法线性规划模型,并计算岩块的失效概率。本发明以二维岩质边坡为研究对象,将塑性极限分析上限定理、刚性块体单元离散技术、线性规划方法和蒙特卡洛方法结合起来建立岩块失效概率计算的上限法数学规划模型;通过该模型可以有效地用于统计岩块的失效概率,并获得岩质边坡中所有岩块的失效概率和对应的失效模式。
Description
技术领域
本发明涉及一种岩质边坡中岩块失效概率的计算方法,属于岩质边坡工程可靠度分析技术领域。
背景技术
在岩土工程领域进行岩质边坡的可靠度分析已经非常流行,但在此领域还存在一些不足,比如:(1)当前的岩质边坡可靠度分析中使用的边坡稳定性计算方法主要以刚体极限平衡法、有限单元法为主,然而刚体极限平衡法需要事先假设岩质边坡破坏的滑裂面、有限单元法需要模拟复杂的节理网络和本构关系,这些都影响了岩质边坡可靠度分析的准确性;(2)在岩质边坡的失效概率方面,主要是基于岩质边坡整体稳定的安全系数以计算边坡整体的失效概率为主;(3)由于参数的变异性,岩质边坡一般会存在多种的失效模式,当前研究成果对岩质边坡多种失效模式的研究存在不足;
鉴于此,本发明基于国家自然科学基金项目(51564026)进行研究工作。
发明内容
本发明提供了一种岩质边坡中岩块失效概率的计算方法,以获得岩质边坡中岩块失效概率,为岩块失效概率的计算提供一种新的方法。
本发明的技术方案是:一种岩质边坡中岩块失效概率的计算方法,所述方法具体步骤如下:
步骤1、拟定岩质边坡的计算参数;
步骤2、建立岩质边坡的极限状态函数;
步骤3、使用刚性块体单元离散岩质边坡;
步骤4、建立求解岩块失效概率的上限法线性规划模型;
步骤5、使用蒙特卡洛法求解上限法线性规划模型,并计算岩块的失效概率。
所述岩质边坡的计算参数包括岩质边坡的几何参数、岩质边坡的地质参数、岩体材料的物理力学参数、结构面抗剪参数统计值、边坡的边界荷载信息。
所述岩质边坡的极限状态函数为:
式中:Z是岩质边坡的极限状态函数,λγ(t)是与第t个抗剪参数的蒙特卡洛随机量相对应的第t个容重超载系数的随机量;是与第t个抗剪参数随机变量相关的极限容重随机量;cr(t)是结构面凝聚力的第t个蒙特卡洛随机量;是结构面摩擦角的第t个蒙特卡洛随机量;t=(1,…,N),N是抗剪参数的蒙特卡洛随机量的数量。
所述步骤3具体为:在总体坐标系(x,y)下,任意一个刚性块体单元i形心Ci处作用有速度向量ui=[ui,vi]T,其中ui是刚性块体单元i形心Ci处沿x方向的速度、vi是刚性块体单元i形心Ci处沿y方向的速度,相邻两个刚性块体单元i和刚性块体单元j之间的结构面k的局部坐标系为(nk,sk),使用刚性块体单元离散岩质边坡以后,岩质边坡变成刚性块体单元+结构面的几何系统;其中nk沿结构面法线方向并与x轴正向的夹角为αk,sk沿结构面切线方向,nk轴与sk轴符合右手系法则。
所述步骤4包括:①建立结构面的塑性流动约束条件;②建立刚性块体单元的速度边界条件;③建立目标函数;④建立岩块的失效概率计算公式;⑤根据极限状态函数、①、②、③和④建立的函数建立求解岩块失效概率的上限法线性规划模型。
所述步骤4具体为:
①建立结构面的塑性流动约束条件为:
Dkuij-Nkλk=0(1)
式中: K1是岩质边坡中所有结构面的数量;uij=[ui,vi,uj,vj]T,ui是刚性块体单元i形心Ci处沿x方向的速度、vi是刚性块体单元i形心Ci处沿y方向的速度,uj是刚性块体单元j形心Cj处沿x方向的速度、vj是刚性块体单元j形心Cj处沿y方向的速度;是结构面k的摩擦角的第t个蒙特卡洛随机量;分别是结构面k的两个非负塑性乘子,αk是结构面k的局部坐标系的nk轴与总体坐标系的x轴之间的夹角且满足逆时针为正;
②建立刚性块体单元的速度边界条件:
Tbub=0(2)
式中:Tb是边界b的坐标转换矩阵;ub是边界b的速度向量;b=(1,…,K3),K3是岩质边坡中边界速度等于0的边界单元的总数量;
③目标函数:
在考虑容重超载时,整个岩质边坡的外功功率等于内功功率的约束条件如下:
式中:K1是岩质边坡中所有结构面的数量,K2是岩质边坡中所有刚性块体单元的数量,是结构面k的凝聚力的第t个蒙特卡洛随机量,lk是结构面k的长度,Ai是刚性块体单元i的面积,γa是岩质边坡当前实际的容重,vi是刚性块体单元i形心Ci处沿y方向的速度,分别是结构面k的两个非负塑性乘子;
假设则式(3)等价于以下两式:
将容重超载系数作为目标函数,根据上限定理需要求解容重超载系数的最小值;结合式(4),定义目标函数如下:
上式中:Minimize表示求“最小”;
④岩块的失效概率计算公式:
刚性块体单元i的失效功能函数如下:
式中:i=(1,…,K2),K2是岩质边坡中所有刚性块体单元的数量;t=(1,…,N),N是抗剪参数的蒙特卡洛随机量的数量,ui(t)是第t次蒙特卡洛计算的刚性块体单元i形心处的速度向量;λγ(t)是第t次蒙特卡洛计算的岩质边坡的容重超载系数的随机量,“or”表示逻辑运算“或者”,“and”表示逻辑运算“与”;
根据式(7),岩质边坡中每个岩块的失效概率按下式计算:
式中:Pf,i是刚性块体单元i的失效概率;i=(1,…,K2),K2是岩质边坡中所有刚性块体单元的数量;
⑤建立求解岩块失效概率的上限法线性规划模型:
所述步骤5具体为:
S1、使用蒙特卡洛方法生成结构面的N个凝聚力的蒙特卡洛随机量[cr(t),(t=1,…,N)]和内摩擦角的蒙特卡洛随机量然后将cr(t)、从t=1到t=N循环,逐次代入求解岩块失效概率的上限法线性规划模型,使用“对偶单纯形法”求解线性规划问题,求解得到N个容重超载系数的随机变量[λγ(t),(t=1,…,N)]和与之对应的刚性块体单元的速度场;
S2、根据容重超载系数的随机量[λγ(t),(t=1,…,N)]和刚性块体单元的速度场,计算所有刚性块体单元的失效功能函数[Ii(t),i=(1,…,K2),t=(1,…,N)];
S3、根据所有刚性块体单元的失效功能函数[Ii(t),i=(1,…,K2),t=(1,…,N)]计算每个刚性块体单元的失效概率[Pf,i,i=(1,…,K2)],并统计分析岩质边坡的失效模式。
本发明的有益效果是:本发明以二维岩质边坡为研究对象,将塑性极限分析上限定理、刚性块体单元离散技术、线性规划方法和蒙特卡洛方法结合起来建立岩块失效概率计算的上限法数学规划模型;通过该模型可以有效地用于统计岩块的失效概率,并获得岩质边坡中所有岩块的失效概率和对应的失效模式。
附图说明
图1本发明的技术路线图;
图2刚性块体单元示意图;
图3相邻刚性块体单元之间的结构面示意图;
图4实施例1土岩质边坡的几何形状示意图(单位:m);
图5实施例1岩质边坡刚性块体单元离散示意图;
图6实施例1岩质边坡失效模式1速度矢量图;
图7实施例1岩质边坡失效模式2速度矢量图。
具体实施方式
实施例1:如图1-7所示,一种岩质边坡中岩块失效概率的计算方法,所述方法具体步骤如下:
步骤1、拟定岩质边坡的计算参数;
步骤2、建立岩质边坡的极限状态函数;
步骤3、使用刚性块体单元离散岩质边坡;
步骤4、建立求解岩块失效概率的上限法线性规划模型;
步骤5、使用蒙特卡洛法求解上限法线性规划模型,并计算岩块的失效概率。
流程图如图1所示。
进一步地,可以设置所述方法步骤如下:
步骤一、拟定岩质边坡的基本信息:根据二维岩质边坡的工程资料,拟定进行计算分析的主要计算参数,包括:岩质边坡的几何参数(边坡的高度、坡顶宽度、坡面的坡角等)、岩质边坡的地质参数(结构面的倾角)、岩体材料的物理力学参数(岩体的容重等)、结构面抗剪参数统计值(结构面的凝聚力和摩擦角的均值、标准差、变异系数)、边坡的边界荷载信息。
二维岩质边坡如图4所示,岩质边坡的几何参数:边坡HD的高度30m、坡顶FD宽度17.7m、坡面AF的坡角50°;岩质边坡的地质参数:结构面AD的倾角35°、结构面BE的倾角90°、结构面CG的倾角75°;岩体材料的物理力学参数:岩体的容重27kN/m3;三条结构面的凝聚力的均值取65kPa、三条结构面的凝聚力变异系数取0.3,三条结构面的摩擦角的均值取34°、三条结构面摩擦角的变异系数取0.3;边坡无其他外荷载作用。
步骤二、建立岩质边坡的极限状态函数。
本发明使用容重超载的方式使岩质边坡达到极限状态,定义岩质边坡的容重超载系数如下:
上式中:λγ是边坡的容重超载系数;γa是岩质边坡当前实际的容重,γc是边坡达到极限状态时的极限容重。
为了求解岩质边坡容重超载系数的分布规律,本发明设岩质边坡的极限状态函数为:
上式中:Z是岩质边坡的极限状态函数,是与抗剪参数随机变量相关的容重超载系数的随机变量;是与抗剪参数随机变量相关的极限容重随机变量;cr是结构面凝聚力的随机变量,是结构面摩擦角的随机变量。
本发明将岩质边坡中结构面的凝聚力cr和摩擦角设为相互独立的随机变量,并假设两者均符合对数正态分布,并假设边坡其它的参数为确定值;然后根据结构面凝聚力和摩擦角的均值、变异系数、标准差,使用蒙特卡洛方法生成结构面的N个凝聚力的蒙特卡洛随机量[cr(t),(t=1,…,N)]和内摩擦角的蒙特卡洛随机量则式(2)可写为:
上式中:λγ(t)是与第t个抗剪参数的蒙特卡洛随机量相对应的第t个容重超载系数的随机量;是与第t个抗剪参数随机变量相关的极限容重随机量;cr(t)是结构面凝聚力的第t个蒙特卡洛随机量;是结构面摩擦角的第t个蒙特卡洛随机量;t=(1,…,N),N是抗剪参数的蒙特卡洛随机量的数量。
具体的,将岩质边坡中结构面的凝聚力和摩擦角设为相互独立的随机变量,并假设两者均符合对数正态分布,并设岩体容重为确定值;据结构面凝聚力和摩擦角的均值、变异系数、标准差,使用蒙特卡洛方法生成结构面的N=10000个凝聚力的蒙特卡洛随机量[cr(t),(t=1,…,10000)]和内摩擦角的蒙特卡洛随机量并根据式(3)建立岩质边坡的极限状态函数Z。
步骤三、使用刚性块体单元离散岩质边坡:岩质边坡由于受到结构面(节理、裂隙等)的切割而被分成若干岩块,岩质边坡的承载能力主要受控于结构面的几何分布和强度。本发明使用刚性块体单元来离散岩质边坡,任意一个刚性块体单元i如图2所示,在总体坐标系(x,y)下,刚性块体单元i形心Ci处作用有速度向量ui=[ui,vi]T,其中ui是刚性块体单元i形心Ci处沿x方向的速度、vi是刚性块体单元i形心Ci处沿y方向的速度。相邻两个刚性块体单元i和刚性块体单元j之间的结构面k如图3所示,定义结构面k的局部坐标系为(nk,sk),其中nk沿结构面法线方向并与x轴正向的夹角为αk、sk沿结构面切线方向,nk轴与sk轴符合右手系法则。使用刚性块体单元离散岩质边坡以后,岩质边坡变成刚性块体单元+结构面的几何系统。
使用刚性块体单元离散岩质边坡(如图5所示),共计离散为4个刚性块体单元、5条结构面,刚性块体单元4为边界上速度等于0(固定不动)的单元,即K1=5、K2=4、K3=4。
步骤四、建立求解岩块失效概率的上限法线性规划模型。
1.结构面的塑性流动约束条件
相邻刚性块体单元i和刚性块体单元j之间的结构面k如图3所示,结构面k的变形必须满足变形协调条件,同时结构面k上的不连续的法向和切向速度间断值必须符合关联流动准则。经过推导,结构面k的塑性流动约束条件可为:
Dkuij-Nkλk=0(4)
式中:
K1是岩质边坡中所有结构面的数量;uij=[ui,vi,uj,vj]T,ui是刚性块体单元i形心Ci处沿x方向的速度、vi是刚性块体单元i形心Ci处沿y方向的速度,uj是刚性块体单元j形心Cj处沿x方向的速度、vj是刚性块体单元j形心Cj处沿y方向的速度;是结构面k的摩擦角的第t个蒙特卡洛随机量;分别是结构面k的两个非负塑性乘子,αk是结构面k的局部坐标系的nk轴与总体坐标系的x轴之间的夹角(逆时针为正)。
2.刚性块体单元的速度边界条件
根据上限定理,岩质边坡边界上的刚性块体单元应满足已知的速度边界条件,速度等于0的边界上的块体单元的速度边界条件为:
Tbub=0(5)
式中:Tb是边界b的坐标转换矩阵;ub是边界b的速度向量;b=(1,…,K3),K3是岩质边坡中边界速度等于0的边界单元的总数量。
3、目标函数
在考虑容重超载时,整个岩质边坡的外功功率等于内功功率的约束条件如下:
上式中:K1是岩质边坡中所有结构面的数量,K2是岩质边坡中所有刚性块体单元的数量,是结构面k的凝聚力的第t个蒙特卡洛随机量,lk是结构面k的长度,Ai是刚性块体单元i的面积,γa是岩质边坡当前实际的容重,vi是刚性块体单元i形心Ci处沿y方向的速度,分别是结构面k的两个非负塑性乘子。
式(6)是一个非线性表达式,为了避免求解非线性数学规划问题,本发明假设则式(6)等价于以下两式:
在建立上限法线性规划模型时,本发明将容重超载系数作为目标函数,根据上限定理需要求解容重超载系数的最小值。结合式(7),定义目标函数如下:
上式中:Minimize表示求“最小”。
4.岩块的失效概率计算公式
根据上限法理论,当刚性块体单元形心处的速度大于0时,其相对于边界上速度为0的刚性块体单元会产生相对速度,此时刚性块体单元发生失效。因此,当刚性块体单元同时满足以下两个条件时岩块发生失效:(1)岩质边坡容重超载系数小于1.0;(2)刚性块体单元形心的速度大于0。本发明定义刚性块体单元i的失效功能函数如下:
上式中:i=(1,…,K2),K2是岩质边坡中所有刚性块体单元的数量;t=(1,…,N),ui(t)是第t次蒙特卡洛计算的刚性块体单元i形心处的速度向量;λγ(t)是第t次蒙特卡洛计算的岩质边坡的容重超载系数的随机量,“or”表示逻辑运算“或者”,“and”表示逻辑运算“与”。
根据式(10),岩质边坡中每个岩块(刚性块体单元)的失效概率可按下式计算:
上式中:Pf,i是刚性块体单元i的失效概率(单位:%);i=(1,…,K2),K2是岩质边坡中所有刚性块体单元的数量。
5.求解岩块失效概率的上限法线性规划模型
根据极限状态函数式(3),目标函数式(9),约束条件式(4)、(5)、(8),失效概率计算公式(10)、(11),可得到求解岩块失效概率的上限法线性规划模型如下:
根据式(12)建立实施例1的求解岩块失效概率的上限法线性规划模型。
步骤五、使用蒙特卡洛法求解实施例上限法线性规划模型,并计算岩块的失效概率:
S1、使用蒙特卡洛方法生成结构面的N个凝聚力的蒙特卡洛随机量[cr(t),(t=1,…,N)]和内摩擦角的蒙特卡洛随机量然后将cr(t)、从t=1到t=N循环,逐次代入求解岩块失效概率的上限法线性规划模型,使用“对偶单纯形法”求解线性规划问题,求解得到N个容重超载系数的随机变量[λγ(t),(t=1,…,N)]和与之对应的刚性块体单元的速度场;
S2、根据容重超载系数的随机量[λγ(t),(t=1,…,N)]和刚性块体单元的速度场,计算所有刚性块体单元的失效功能函数[Ii(t),i=(1,…,K2),t=(1,…,N)];
S3、根据所有刚性块体单元的失效功能函数[Ii(t),i=(1,…,K2),t=(1,…,N)]计算每个刚性块体单元的失效概率[Pf,i,i=(1,…,K2)],并统计分析岩质边坡的失效模式。
计算结果是:岩块1的失效概率是9.07%、岩块2和岩块3的失效概率均是0.45%、岩块4的失效概率0.00%;在10000次蒙特卡洛模拟计算的结果中,经过统计分析此岩质边坡存在两种失效模式(分别如图6、图7所示),失效模式1中只有岩块1发生失效,失效次数为862次;失效模式2中岩块1、岩块2、岩块3同时发生失效,失效次数为45次;
表1实施例1岩质边坡的岩块失效概率统计表
上面结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (7)
1.一种岩质边坡中岩块失效概率的计算方法,其特征在于:所述方法具体步骤如下:
步骤1、拟定岩质边坡的计算参数;
步骤2、建立岩质边坡的极限状态函数;
步骤3、使用刚性块体单元离散岩质边坡;
步骤4、建立求解岩块失效概率的上限法线性规划模型;
步骤5、使用蒙特卡洛法求解上限法线性规划模型,并计算岩块的失效概率。
2.根据权利要求1所述的岩质边坡中岩块失效概率的计算方法,其特征在于:所述岩质边坡的计算参数包括岩质边坡的几何参数、岩质边坡的地质参数、岩体材料的物理力学参数、结构面抗剪参数统计值、边坡的边界荷载信息。
3.根据权利要求1所述的岩质边坡中岩块失效概率的计算方法,其特征在于:所述岩质边坡的极限状态函数为:
式中:Z是岩质边坡的极限状态函数,λγ(t)是与第t个抗剪参数的蒙特卡洛随机量相对应的第t个容重超载系数的随机量;是与第t个抗剪参数随机变量相关的极限容重随机量;cr(t)是结构面凝聚力的第t个蒙特卡洛随机量;是结构面摩擦角的第t个蒙特卡洛随机量;t=(1,…,N),N是抗剪参数的蒙特卡洛随机量的数量。
4.根据权利要求1所述的岩质边坡中岩块失效概率的计算方法,其特征在于:所述步骤3具体为:在总体坐标系(x,y)下,任意一个刚性块体单元i形心Ci处作用有速度向量ui=[ui,vi]T,其中ui是刚性块体单元i形心Ci处沿x方向的速度、vi是刚性块体单元i形心Ci处沿y方向的速度,相邻两个刚性块体单元i和刚性块体单元j之间的结构面k的局部坐标系为(nk,sk),使用刚性块体单元离散岩质边坡以后,岩质边坡变成刚性块体单元+结构面的几何系统;其中nk沿结构面法线方向并与x轴正向的夹角为αk,sk沿结构面切线方向,nk轴与sk轴符合右手系法则。
5.根据权利要求1所述的岩质边坡中岩块失效概率的计算方法,其特征在于:所述步骤4包括:①建立结构面的塑性流动约束条件;②建立刚性块体单元的速度边界条件;③建立目标函数;④建立岩块的失效概率计算公式;⑤根据极限状态函数、①、②、③和④建立的函数建立求解岩块失效概率的上限法线性规划模型。
6.根据权利要求5所述的岩质边坡中岩块失效概率的计算方法,其特征在于:所述步骤4具体为:
①建立结构面的塑性流动约束条件为:
Dkuij-Nkλk=0 (1)
式中:
K1是岩质边坡中所有结构面的数量;uij=[ui,vi,uj,vj]T,ui是刚性块体单元i形心Ci处沿x方向的速度、vi是刚性块体单元i形心Ci处沿y方向的速度,uj是刚性块体单元j形心Cj处沿x方向的速度、vj是刚性块体单元j形心Cj处沿y方向的速度;是结构面k的摩擦角的第t个蒙特卡洛随机量;分别是结构面k的两个非负塑性乘子,αk是结构面k的局部坐标系的nk轴与总体坐标系的x轴之间的夹角且满足逆时针为正;
②建立刚性块体单元的速度边界条件:
Tbub=0 (2)
式中:Tb是边界b的坐标转换矩阵;ub是边界b的速度向量;b=(1,…,K3),K3是岩质边坡中边界速度等于0的边界单元的总数量;
③目标函数:
在考虑容重超载时,整个岩质边坡的外功功率等于内功功率的约束条件如下:
式中:K1是岩质边坡中所有结构面的数量,K2是岩质边坡中所有刚性块体单元的数量,是结构面k的凝聚力的第t个蒙特卡洛随机量,lk是结构面k的长度,Ai是刚性块体单元i的面积,γa是岩质边坡当前实际的容重,vi是刚性块体单元i形心Ci处沿y方向的速度,分别是结构面k的两个非负塑性乘子;
假设则式(3)等价于以下两式:
将容重超载系数作为目标函数,根据上限定理需要求解容重超载系数的最小值;结合式(4),定义目标函数如下:
上式中:Minimize表示求“最小”;
④岩块的失效概率计算公式:
刚性块体单元i的失效功能函数如下:
式中:i=(1,…,K2),K2是岩质边坡中所有刚性块体单元的数量;t=(1,…,N),N是抗剪参数的蒙特卡洛随机量的数量,ui(t)是第t次蒙特卡洛计算的刚性块体单元i形心处的速度向量;λγ(t)是第t次蒙特卡洛计算的岩质边坡的容重超载系数的随机量,“or”表示逻辑运算“或者”,“and”表示逻辑运算“与”;
根据式(7),岩质边坡中每个岩块的失效概率按下式计算:
式中:Pf,i是刚性块体单元i的失效概率;i=(1,…,K2),K2是岩质边坡中所有刚性块体单元的数量;
⑤建立求解岩块失效概率的上限法线性规划模型:
7.根据权利要求1所述的岩质边坡中岩块失效概率的计算方法,其特征在于:所述步骤5具体为:
S1、使用蒙特卡洛方法生成结构面的N个凝聚力的蒙特卡洛随机量[cr(t),(t=1,…,N)]和内摩擦角的蒙特卡洛随机量然后将cr(t)、从t=1到t=N循环,逐次代入求解岩块失效概率的上限法线性规划模型,使用“对偶单纯形法”求解线性规划问题,求解得到N个容重超载系数的随机变量[λγ(t),(t=1,…,N)]和与之对应的刚性块体单元的速度场;
S2、根据容重超载系数的随机量[λγ(t),(t=1,…,N)]和刚性块体单元的速度场,计算所有刚性块体单元的失效功能函数[Ii(t),i=(1,…,K2),t=(1,…,N)];
S3、根据所有刚性块体单元的失效功能函数[Ii(t),i=(1,…,K2),t=(1,…,N)]计算每个刚性块体单元的失效概率[Pf,i,i=(1,…,K2)],并统计分析岩质边坡的失效模式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810478606.4A CN108763697B (zh) | 2018-05-18 | 2018-05-18 | 一种岩质边坡中岩块失效概率的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810478606.4A CN108763697B (zh) | 2018-05-18 | 2018-05-18 | 一种岩质边坡中岩块失效概率的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108763697A true CN108763697A (zh) | 2018-11-06 |
CN108763697B CN108763697B (zh) | 2022-05-20 |
Family
ID=64007127
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810478606.4A Active CN108763697B (zh) | 2018-05-18 | 2018-05-18 | 一种岩质边坡中岩块失效概率的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108763697B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109753700A (zh) * | 2018-12-21 | 2019-05-14 | 昆明理工大学 | 一种随机库水位作用下的土石坝可靠度分析上限法 |
CN109902345A (zh) * | 2019-01-23 | 2019-06-18 | 中国水利水电科学研究院 | 基于极限分析下限定理的二维边坡极限荷载求解方法 |
CN111611728A (zh) * | 2020-06-29 | 2020-09-01 | 中国有色金属工业昆明勘察设计研究院有限公司 | 昔格达地层岩质边坡后缘滑裂面位置的稳定性计算方法 |
CN111651900A (zh) * | 2020-06-29 | 2020-09-11 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种昔格达地层岩质边坡稳定性计算的拟动力上限法 |
CN112115530A (zh) * | 2020-08-05 | 2020-12-22 | 昆明理工大学 | 一种土质边坡的点失效概率的计算方法 |
CN113012399A (zh) * | 2021-02-25 | 2021-06-22 | 中原工学院 | 一种降雨型滑坡预警方法及系统 |
CN113832992A (zh) * | 2021-09-09 | 2021-12-24 | 武汉大学 | 基于边坡系统失稳风险的抗滑桩设计优化方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090248312A1 (en) * | 2008-03-31 | 2009-10-01 | Sheng-Yuan Hsu | Integration of geomechanics and seismic analysis for passive seismic feasibility analysis |
CN104318121A (zh) * | 2014-11-07 | 2015-01-28 | 中铁第四勘察设计院集团有限公司 | 一种路基稳定可靠度的计算方法 |
CN106126796A (zh) * | 2016-06-20 | 2016-11-16 | 昆明理工大学 | 一种带岩桥的岩质边坡极限承载力的塑性极限分析下限法 |
CN106372314A (zh) * | 2016-08-31 | 2017-02-01 | 广东省交通规划设计研究院股份有限公司 | 一种线性工程岩质边坡稳定性半定量自动评价系统 |
CN107220401A (zh) * | 2017-04-12 | 2017-09-29 | 中国地质大学(武汉) | 基于并行蒙特卡洛法的边坡可靠性参数获取方法及装置 |
CN107330146A (zh) * | 2017-05-26 | 2017-11-07 | 昆明理工大学 | 一种同时考虑平动和转动效应的节理岩质边坡极限承载力分析上限法 |
CN107330145A (zh) * | 2017-05-26 | 2017-11-07 | 昆明理工大学 | 一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法 |
CN107992687A (zh) * | 2017-12-06 | 2018-05-04 | 昆明理工大学 | 基于蒙特卡洛方法的岩质边坡可靠度塑性极限分析上限法 |
-
2018
- 2018-05-18 CN CN201810478606.4A patent/CN108763697B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090248312A1 (en) * | 2008-03-31 | 2009-10-01 | Sheng-Yuan Hsu | Integration of geomechanics and seismic analysis for passive seismic feasibility analysis |
CN104318121A (zh) * | 2014-11-07 | 2015-01-28 | 中铁第四勘察设计院集团有限公司 | 一种路基稳定可靠度的计算方法 |
CN106126796A (zh) * | 2016-06-20 | 2016-11-16 | 昆明理工大学 | 一种带岩桥的岩质边坡极限承载力的塑性极限分析下限法 |
CN106372314A (zh) * | 2016-08-31 | 2017-02-01 | 广东省交通规划设计研究院股份有限公司 | 一种线性工程岩质边坡稳定性半定量自动评价系统 |
CN107220401A (zh) * | 2017-04-12 | 2017-09-29 | 中国地质大学(武汉) | 基于并行蒙特卡洛法的边坡可靠性参数获取方法及装置 |
CN107330146A (zh) * | 2017-05-26 | 2017-11-07 | 昆明理工大学 | 一种同时考虑平动和转动效应的节理岩质边坡极限承载力分析上限法 |
CN107330145A (zh) * | 2017-05-26 | 2017-11-07 | 昆明理工大学 | 一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法 |
CN107992687A (zh) * | 2017-12-06 | 2018-05-04 | 昆明理工大学 | 基于蒙特卡洛方法的岩质边坡可靠度塑性极限分析上限法 |
Non-Patent Citations (3)
Title |
---|
A.J.LI等: "Rock slope stability analyses using extreme learning neural network and terminal steepest descent algorithm", 《AUTOMATION IN CONSTRUCTION》 * |
张小艳等: "基于塑性极限分析上限法理论的土质边坡可靠度分析", 《岩土力学》 * |
苏树明: "基于蒙特卡洛法与极限分析上限法的边坡失效概率分析", 《路基工程》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109753700A (zh) * | 2018-12-21 | 2019-05-14 | 昆明理工大学 | 一种随机库水位作用下的土石坝可靠度分析上限法 |
CN109753700B (zh) * | 2018-12-21 | 2022-09-02 | 昆明理工大学 | 一种土石坝可靠度分析上限法 |
CN109902345A (zh) * | 2019-01-23 | 2019-06-18 | 中国水利水电科学研究院 | 基于极限分析下限定理的二维边坡极限荷载求解方法 |
CN109902345B (zh) * | 2019-01-23 | 2021-02-12 | 中国水利水电科学研究院 | 基于极限分析下限定理的二维边坡极限荷载求解方法 |
CN111611728A (zh) * | 2020-06-29 | 2020-09-01 | 中国有色金属工业昆明勘察设计研究院有限公司 | 昔格达地层岩质边坡后缘滑裂面位置的稳定性计算方法 |
CN111651900A (zh) * | 2020-06-29 | 2020-09-11 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种昔格达地层岩质边坡稳定性计算的拟动力上限法 |
CN111651900B (zh) * | 2020-06-29 | 2023-07-04 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种昔格达地层岩质边坡稳定性计算的拟动力上限法 |
CN112115530A (zh) * | 2020-08-05 | 2020-12-22 | 昆明理工大学 | 一种土质边坡的点失效概率的计算方法 |
CN112115530B (zh) * | 2020-08-05 | 2022-12-30 | 昆明理工大学 | 一种土质边坡的点失效概率的计算方法 |
CN113012399A (zh) * | 2021-02-25 | 2021-06-22 | 中原工学院 | 一种降雨型滑坡预警方法及系统 |
CN113832992A (zh) * | 2021-09-09 | 2021-12-24 | 武汉大学 | 基于边坡系统失稳风险的抗滑桩设计优化方法 |
CN113832992B (zh) * | 2021-09-09 | 2023-02-03 | 武汉大学 | 基于边坡系统失稳风险的抗滑桩设计优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108763697B (zh) | 2022-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108763697A (zh) | 一种岩质边坡中岩块失效概率的计算方法 | |
CN106126796B (zh) | 一种带岩桥的岩质边坡极限承载力的塑性极限分析下限法 | |
Cascini et al. | LiABlock_3D: a software tool for collapse mechanism analysis of historic masonry structures | |
CN107330145B (zh) | 一种考虑转动效应的岩质边坡极限承载力分析下限法 | |
Zhang et al. | BIM-and 4D-based integrated solution of analysis and management for conflicts and structural safety problems during construction: 1. Principles and methodologies | |
Malena et al. | Collapse mechanism analysis of historic masonry structures subjected to lateral loads: A comparison between continuous and discrete models | |
CN100573593C (zh) | 滑坡灾变过程时空预测的智能方法 | |
Zheng et al. | An improved Monte Carlo simulation method for discontinuity orientations based on Fisher distribution and its program implementation | |
CN104239687B (zh) | 一种基于航天产品信号传递路径的可靠性建模与评估方法 | |
CN107992687A (zh) | 基于蒙特卡洛方法的岩质边坡可靠度塑性极限分析上限法 | |
CN108170899A (zh) | 一种土质边坡可靠度分析上限法 | |
CN104281730B (zh) | 一种大转动变形的板壳结构动响应的有限元分析方法 | |
CN105893774A (zh) | 生物多样性保护优先区域绿色发展生态安全评价指数的计算机建模方法 | |
CN107194024B (zh) | 基于支持向量机的边坡可靠性参数获取方法及装置 | |
CN108170898B (zh) | 一种节理岩质边坡可靠度分析下限法 | |
Xie et al. | Multi-city sustainable regional urban growth simulation—MSRUGS: a case study along the mid-section of Silk Road of China | |
Wang et al. | Selection of an appropriate mechanized mining technical process for thin coal seam mining | |
CN107391882A (zh) | 一种拱坝的拱座抗滑稳定安全度计算方法 | |
CN113761732A (zh) | 一种基于强化学习的一类多扰动车间柔性调度建模与优化的方法 | |
Qu et al. | An explicit time integration scheme of numerical manifold method | |
Sun et al. | A three-dimensional procedure for evaluating the stability of gravity dams against deep slide in the foundation | |
Wang et al. | A novel evaluation method for the stability of construction sites on an abandoned goaf: a case study | |
Korolev et al. | Basic stages of creating a BIM model for transport infrastructure objects | |
CN113221399A (zh) | 一种铁路隧道装配式衬砌构件划分方案决策方法 | |
Yang et al. | Two-dimensional stability analysis of a soil slope using the finite element method and the limit equilibrium principle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |