CN112685817A - 一种对锚固边坡风险进行量化分析的方法 - Google Patents
一种对锚固边坡风险进行量化分析的方法 Download PDFInfo
- Publication number
- CN112685817A CN112685817A CN202011572821.4A CN202011572821A CN112685817A CN 112685817 A CN112685817 A CN 112685817A CN 202011572821 A CN202011572821 A CN 202011572821A CN 112685817 A CN112685817 A CN 112685817A
- Authority
- CN
- China
- Prior art keywords
- slope
- rock
- anchoring
- sliding surface
- area
- 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
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
- E02D17/00—Excavations; Bordering of excavations; Making embankments
- E02D17/20—Securing of slopes or inclines
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
- E02D5/00—Bulkheads, piles, or other structural elements specially adapted to foundation engineering
- E02D5/74—Means for anchoring structural elements or bulkheads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
一种对锚固边坡风险进行量化分析的方法,涉及安全风险评估技术领域,包括如下步骤:构建施加锚杆后的边坡数学模型M0,计算M0的安全系数FS0,并找出滑动面S0;在M0的基础上,标记岩土体受锚固区影响的区域,撤去锚杆,增强受影响区域的力学参数,建立一个新模型Mx,利用Geo‑studio找出受锚固区影响的区域岩土体力学参数的范围值;将Mx运用光滑粒子流体动力学方法进行模拟,确定受锚杆影响的区域的岩土体力学参数的准确值;再次利用光滑粒子流体动力学方法进行模拟,结合蒙特卡罗方法对锚固边坡风险进行量化分析。本发明提出了一种对锚固边坡风险进行量化分析的方法,该方法可以对锚固边坡的风险进行量化分析。
Description
技术领域
本发明涉及安全风险评估技术领域,具体涉及一种对锚固边坡风险进行量化分析的方法。
背景技术
发明内容
为解决现有技术中存在的问题,本发明提出了一种对锚固边坡风险进行量化分析的方法,该方法可以对锚固边坡的风险进行量化分析。
为了达到上述目的,本发明的技术方案为:
一种对锚固边坡风险进行量化分析的方法,包括如下步骤:
步骤1、构建施加锚杆后的边坡数学模型M0,计算边坡数学模型M0的安全系数FS0,并找出此时安全系数FS0下的滑动面S0;
步骤2、在边坡数学模型M0的基础上,标记出边坡岩土体受锚固区影响的区域,撤去锚杆,增强受锚固区影响的区域的岩土体力学参数,建立一个新模型Mx,基于滑动面S0在边坡数学模型M0与新模型Mx下安全系数相等的原理,利用Geo-studio找出受锚固区影响的区域岩土体力学参数的范围值;
步骤3、将新模型Mx运用光滑粒子流体动力学方法进行模拟,根据其初始滑动面与滑动面S0一致,可以确定受锚杆影响的区域的岩土体力学参数的准确值;
步骤4、再次利用光滑粒子流体动力学方法进行模拟,结合蒙特卡罗方法对锚固边坡风险进行量化分析。
优选的,所述的步骤1中,所述的安全系数FS0通过Geo-studio中的极限平衡方法计算。
优选的,所述的步骤2中,将锚杆作用效果等效为岩土体自身强度增强,以此为依据增强受锚固区影响的区域的岩土体力学参数。
优选的,所述的步骤2中,确定受锚固区影响的区域岩土体力学参数的范围值的方法为:
优选的,所述的步骤3中,将新模型Mx运用光滑粒子流体动力学方法进行模拟时,用完全指定滑动面的方法计算与参数样本X1,X2,…,XN相应的安全系数Fs1,Fs2,…,FsN,并指定初始滑动面Si,且初始滑动面Si与滑动面S0相同,找出Fsi=FS0的样本记为Z1,Z2,…,Zm,并得到相应的受锚固区影响的区域的岩土体粘聚力和摩擦角ci、其中i=1,2,…,m,m≤N。
优选的,所述的步骤3中,确定受锚杆影响的区域的岩土体力学参数的准确值的方法为:
优选的,所述的步骤4中,再次利用光滑粒子流体动力学方法进行模拟时,将力学参数值cf,输入到所要计算的锚固边坡中,通过计算边坡施加锚杆后边坡土体的滑动方量,结合蒙特卡罗抽样方法计算边坡的失效概率,以此对锚固边坡风险进行量化分析。
本发明一种对锚固边坡风险进行量化分析的方法具有如下有益效果:
1.本发明对锚固边坡的风险分析提供了技术解决方案,通过考虑到锚固效应对边坡土体强度的影响,使锚固边坡的风险分析更具有实际参考意义。
2.本发明相比起现有技术的方法,通过考虑受锚固区影响的区域的岩土强度,不但能够确定锚固边坡的滑移体积,还可以确定锚固边坡失稳后的滑移状态。
3.通过对锚固边坡的风险进行量化评估,可以为锚固边坡的风险评估、安全防护提供有力的依据。
附图说明
图1、本发明流程图;
图2、锚固边坡安全系数和滑动面;
图3、等效替换锚杆边坡锚固区影响岩土体区域示意图;
图4、运用Geo-studio计算等效替换锚杆边坡安全系数;
图5、运用SPH计算等效替换锚固边坡滑动面和滑坡方量;
图6、算例边坡原始图;
图7、锚固边坡安全系数和滑动面S0;
图8、运用Geo-studio寻找等效替换锚杆区域粘聚力和摩擦角范围值;
图9、运用SPH寻找与S0最接近的初始滑动面;
图10、运用Geo-studio计算滑移土体体积;
图11、运用SPH计算等效替代锚杆锚固边坡滑坡方量;
具体实施方式
以下所述,是以阶梯递进的方式对本发明的实施方式详细说明,该说明仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
本发明的描述中,需要说明的是,术语“上”“下”“左”“右”“顶”“底”“内”“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以及特定的方位构造和操作,因此不能理解为对本发明的限制。
参考图1~11:
实施例1、
本发明一种对锚固边坡风险进行量化分析的方法,包括如下步骤:
步骤1、构建施加锚杆后的边坡数学模型M0,计算边坡数学模型M0的安全系数FS0,并找出此时安全系数FS0下的滑动面S0;
步骤2、在边坡数学模型M0的基础上,标记出边坡岩土体受锚固区影响的区域,撤去锚杆,增强受锚固区影响的区域的岩土体力学参数,建立一个新模型Mx,基于滑动面S0在边坡数学模型M0与新模型Mx下安全系数相等的原理,利用Geo-studio找出受锚固区影响的区域岩土体力学参数的范围值;
步骤3、将新模型Mx运用光滑粒子流体动力学方法(Smoothed ParticleHydrodynamics简写为SPH)进行模拟,根据其初始滑动面与滑动面S0一致,可以确定受锚杆影响的区域的岩土体力学参数的准确值;
步骤4、再次利用光滑粒子流体动力学方法进行模拟,结合蒙特卡罗方法对锚固边坡风险进行量化分析;
实施例2、
在实施例1的基础上,本实施例做出了进一步优化,具体为:
如图2所示,所述的步骤1中,所述的安全系数FS0通过Geo-studio中的极限平衡方法计算。
实施例3、
在实施例2的基础上,本实施例做出了进一步优化,具体为:
所述的步骤2中,将锚杆作用效果等效为岩土体自身强度增强,以此为依据增强受锚固区影响的区域的岩土体力学参数;
所述的步骤2中,确定受锚固区影响的区域岩土体力学参数的范围值的方法为:将受锚固区影响的区域的岩土体力学参数视为相互独立的随机变量,利用随机方法产生N个参数样本X1,X2,…,XN,样本相对应的粘聚力和摩擦角为c1、c2、cN、
如图3所示黑色区域,在分析时,将锚杆作用效果等效为岩土体自身强度增强,建立一个新模型Mx,黑色区域为等效替换锚杆后岩土体自身强度增强的区域。
实施例3、
在实施例2的基础上,本实施例做出了进一步优化,具体为:
如图4所示,所述的步骤3中,将新模型Mx运用光滑粒子流体动力学方法进行模拟时,用完全指定滑动面的方法计算与参数样本X1,X2,…,XN相应的安全系数Fs1,Fs2,…,FsN,并指定初始滑动面Si,且初始滑动面Si与滑动面S0相同,找出Fsi=FS0的样本记为Z1,Z2,…,Zm,并得到相应的受锚固区影响的区域的岩土体粘聚力和摩擦角ci、其中i=1,2,…,m,m≤N;
如图5所示,所述的步骤3中,确定受锚杆影响的区域的岩土体力学参数的准确值的方法为:对样本Z1,Z2,…,Zm进行滑动状态模拟,得到其初始滑动面P1,P2,…,Pm;在P1,P2,…,Pm中寻找与滑动面S0最接近的初始滑动面,以此确定锚固区影响区域岩土体的力学参数值cf,
实施例4、
在实施例3的基础上,本实施例做出了进一步优化,具体为:
所述的步骤4中,再次利用光滑粒子流体动力学方法进行模拟时,将力学参数值cf,输入到所要计算的锚固边坡中,通过计算边坡施加锚杆后边坡土体的滑动方量,结合蒙特卡罗抽样方法计算边坡的失效概率,以此对锚固边坡风险进行量化分析。
下面结合图6进行实例说明。
该边坡坡高10m,坡面水平投影长5m,边坡坡率为1:0.5。该边坡岩土体力学参数粘聚力c均值为5kPa、内摩擦角φ均值为25°、容重γ均值为19.5kN/m3;
运用本发明提出的方法对该案例进行分析,根据图6在Geo-Studio的Slope/W模块中建立锚固边坡几何模型M0,将岩土体力学参数的均值输入到模型中计算出此时锚固边坡的安全系数FS0=1.34,并标记出此时的滑动面S0,如图7所示;在原模型基础上,撤去锚杆,建立一个新模型Mx,如图3所示,黑色土条为锚固区影响的岩土体区域,接着将锚固区影响的区域岩土体的粘聚力c和摩擦角φ视为相互独立的随机变量,利用随机方法产生1000个参数样本X1,X2,…,X1000,样本相对应的粘聚力和摩擦角为c1、φ1,c2、φ2,…,c1000、φ1000;利用批处理计算在模型M0的滑动面S0下样本X1,X2,…,X1000相应的安全系数Fs1,Fs2,……,Fs1000,找出Fsi=FS0=1.34的样本(i=1,2,…,m,m≤1000),记为Z1,Z2,…,Z89,如图8所示,并能得到相应的锚固区影响的区域的岩土体的粘聚力c和摩擦角φ的范围值(40kPa≤c≤200kPa,25°≤φ≤55°);
运用SPH方法对样本Z1,Z2,…,Z89进行滑动状态模拟,得到其初始滑动面P1,P2,…,P89;在P1,P2,…,P89中寻找与S0最接近的初始滑动面,如图9所示,此时可以确定锚固区影响的区域岩土体的力学参数值cf=85kPa,φf=36°;
将上述得到的力学参数值cf=85kPa,φf=36°代替图3模型中黑色土体力学参数值,作为计算锚固边坡模型工具对锚固边坡风险量化分析。由于施工和岩土体参数的变异性,将图3所示灰色的土体的粘聚力c和摩擦角φ视为相互独立的随机变量,利用随机方法产生50个参数样本K1,K2,…,K50,再次运用SPH方法,计算边坡施加锚杆后每个边坡的滑移土体方量V1,V2,…,V50,如图11为该样本中某一边坡失稳滑移状态,统计滑移土体方量该50个样本中发生滑动的样本为17个,其中安全系数在0.7~0.8之间有3个,相应的平均滑移土体方量为172.59m3/m,安全系数在0.8~0.9之间有6个,相应的平均滑移土体方量为119.79m3/m,安全系数在0.9~1.0之间有8个,相应的平均滑移土体方量为76.21m3/m,最后统计滑移土体方量V=(172.59*3+119.79*6+76.21*8)/50=36.92m3/m;为了证明本发明方法的有效性,现进行对比分析,如图10所示为运用Geo-studio商用软件基于极限平衡方法计算得到的滑移土体体积为V0=31.07m3/m,本发明方法最终确定的边坡滑移土体体积为V=36.92m3/m,传统方法得到结果偏小。
因此,对比发现:传统的方法虽然能得出锚固边坡滑移体积,但锚杆对土体强度的增强没有考虑,而本发明考虑到锚杆对锚固区影响区域岩土体强度的增强,不但可以确定锚固边坡的滑移体积,而且还能得到锚固边坡失稳后的滑移状态。
Claims (8)
1.一种对锚固边坡风险进行量化分析的方法,其特征为:包括如下步骤:
步骤1、构建施加锚杆后的边坡数学模型M0,计算边坡数学模型M0的安全系数FS0,并找出此时安全系数FS0下的滑动面S0;
步骤2、在边坡数学模型M0的基础上,标记出边坡岩土体受锚固区影响的区域,撤去锚杆,增强受锚固区影响的区域的岩土体力学参数,建立一个新模型Mx,基于滑动面S0在边坡数学模型M0与新模型Mx下安全系数相等的原理,利用Geo-studio找出受锚固区影响的区域岩土体力学参数的范围值;
步骤3、将新模型Mx运用光滑粒子流体动力学方法进行模拟,根据其初始滑动面与滑动面S0一致,可以确定受锚杆影响的区域的岩土体力学参数的准确值;
步骤4、再次利用光滑粒子流体动力学方法进行模拟,结合蒙特卡罗方法对锚固边坡风险进行量化分析。
3.如权利要求2所述的一种对锚固边坡风险进行量化分析的方法,其特征为:所述的步骤1中,所述的安全系数FS0通过Geo-studio中的极限平衡方法计算。
4.如权利要求2或3所述的一种对锚固边坡风险进行量化分析的方法,其特征为:所述的步骤2中,将锚杆作用效果等效为岩土体自身强度增强,以此为依据增强受锚固区影响的区域的岩土体力学参数。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011572821.4A CN112685817B (zh) | 2020-12-24 | 2020-12-24 | 一种对锚固边坡风险进行量化分析的方法 |
PCT/CN2021/091443 WO2022134422A1 (zh) | 2020-12-24 | 2021-04-30 | 一种对锚固边坡风险进行量化分析的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011572821.4A CN112685817B (zh) | 2020-12-24 | 2020-12-24 | 一种对锚固边坡风险进行量化分析的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112685817A true CN112685817A (zh) | 2021-04-20 |
CN112685817B CN112685817B (zh) | 2022-06-14 |
Family
ID=75453445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011572821.4A Active CN112685817B (zh) | 2020-12-24 | 2020-12-24 | 一种对锚固边坡风险进行量化分析的方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN112685817B (zh) |
WO (1) | WO2022134422A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113239434A (zh) * | 2021-05-11 | 2021-08-10 | 青岛理工大学 | 一种确定锚杆衰减系数的新方法 |
WO2022134422A1 (zh) * | 2020-12-24 | 2022-06-30 | 青岛理工大学 | 一种对锚固边坡风险进行量化分析的方法 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115270266B (zh) * | 2022-08-03 | 2023-06-30 | 南昌大学 | 基于改进层次分析法的边坡监测状态风险判断方法 |
CN115048711B (zh) * | 2022-08-15 | 2023-01-10 | 中国矿业大学(北京) | 动静力锚固支护设计方法 |
CN115754225B (zh) * | 2022-10-21 | 2024-03-15 | 西南交通大学 | 基于加速度判断上下坡体相互作用的地震稳定性分析方法 |
CN117807690A (zh) * | 2024-03-01 | 2024-04-02 | 济南轨道交通集团有限公司 | 岩土边坡或基坑锚固体粘结强度的选取方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956317A (zh) * | 2016-05-18 | 2016-09-21 | 青岛理工大学 | 滑坡风险量化方法 |
CN109359361A (zh) * | 2018-09-30 | 2019-02-19 | 青岛理工大学 | 边坡失稳后果量化分析方法 |
CN110765614A (zh) * | 2019-10-24 | 2020-02-07 | 青岛理工大学 | 一种基于滑坡破坏形态的边坡风险综合评估方法 |
AU2020100405A4 (en) * | 2020-03-17 | 2020-04-30 | Qingdao university of technology | A slop risk comprehensive assessment method based on slope failures forms |
WO2020192118A1 (zh) * | 2019-03-28 | 2020-10-01 | 青岛理工大学 | 一种边坡滑动面积的评估方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030143038A1 (en) * | 2000-02-22 | 2003-07-31 | Babcock John W. | Multiple synthetic deformed bars and retaining walls |
CN111914444B (zh) * | 2020-06-30 | 2023-05-26 | 河海大学 | 基于粒子群优化和蒙特卡罗模拟的滑坡有限元区间分析法 |
CN112685817B (zh) * | 2020-12-24 | 2022-06-14 | 青岛理工大学 | 一种对锚固边坡风险进行量化分析的方法 |
-
2020
- 2020-12-24 CN CN202011572821.4A patent/CN112685817B/zh active Active
-
2021
- 2021-04-30 WO PCT/CN2021/091443 patent/WO2022134422A1/zh active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956317A (zh) * | 2016-05-18 | 2016-09-21 | 青岛理工大学 | 滑坡风险量化方法 |
CN109359361A (zh) * | 2018-09-30 | 2019-02-19 | 青岛理工大学 | 边坡失稳后果量化分析方法 |
WO2020192118A1 (zh) * | 2019-03-28 | 2020-10-01 | 青岛理工大学 | 一种边坡滑动面积的评估方法 |
CN110765614A (zh) * | 2019-10-24 | 2020-02-07 | 青岛理工大学 | 一种基于滑坡破坏形态的边坡风险综合评估方法 |
AU2020100405A4 (en) * | 2020-03-17 | 2020-04-30 | Qingdao university of technology | A slop risk comprehensive assessment method based on slope failures forms |
Non-Patent Citations (8)
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022134422A1 (zh) * | 2020-12-24 | 2022-06-30 | 青岛理工大学 | 一种对锚固边坡风险进行量化分析的方法 |
CN113239434A (zh) * | 2021-05-11 | 2021-08-10 | 青岛理工大学 | 一种确定锚杆衰减系数的新方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2022134422A1 (zh) | 2022-06-30 |
CN112685817B (zh) | 2022-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112685817B (zh) | 一种对锚固边坡风险进行量化分析的方法 | |
Yuan et al. | Development of an explicit smoothed particle finite element method for geotechnical applications | |
Cavalieri et al. | Dynamic soil-structure interaction models for fragility characterisation of buildings with shallow foundations | |
Cui et al. | Exploring the macro-and micro-scale response of an idealised granular material in the direct shear apparatus | |
Chandler et al. | Performance-based design in earthquake engineering: a multi-disciplinary review | |
Lam et al. | Time–history analysis of URM walls in out-of-plane flexure | |
Lee | Identification of multiple cracks in a beam using vibration amplitudes | |
Pan et al. | Model-independent strength-reduction factor for effect of spatial variability on tunnel with improved soil surrounds | |
Sadashiva et al. | Quantifying the seismic response of structures with flexible diaphragms | |
Engin et al. | Simplified numerical modelling of pile penetration–the press‐replace technique | |
Alibrandi | A response surface method for stochastic dynamic analysis | |
Wang et al. | Centrifuge testing on monotonic and cyclic lateral behavior of large-diameter slender piles in sand | |
Motra et al. | Quality assessment of soil bearing capacity factor models of shallow foundations | |
Beuth et al. | Large deformation analysis using a quasi-static material point method | |
Grimstad et al. | A non-associated creep model for structured anisotropic clay (n-SAC) | |
Mohapatra et al. | Collapse loads for rectangular foundations by three‐dimensional upper bound limit analysis using radial point interpolation method | |
Abate et al. | Numerical modelling of a shaking table test for soil-foundation-superstructure interaction by means of a soil constitutive model implemented in a FEM code | |
Staubach et al. | Mortar contact discretisation methods incorporating interface models based on Hypoplasticity and Sanisand: Application to vibratory pile driving | |
Wahrhaftig et al. | Buckling and free vibration analysis of non-prismatic columns using optimized shape functions and Rayleigh method | |
Soubra et al. | Bearing capacity of foundations resting on a spatially random soil | |
Bernhardt-Barry et al. | Analysis of the stress distribution in a laminar direct simple shear device and implications for test data interpretation | |
CN112528490A (zh) | 一种岩溶地基中独立基础的极限荷载的计算方法 | |
Rao et al. | Factorized high dimensional model representation for structural reliability analysis | |
Dabeet et al. | Simulation of cyclic direct simple shear loading response of soils using discrete element modeling | |
Hunt et al. | Numerical and physical modeling of the effect of the cone apex angle on the penetration resistance in coarse-grained soils |
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 |