CN110222405B - 一种基于降低强度包线折减法的岩质边坡稳定性测试方法 - Google Patents
一种基于降低强度包线折减法的岩质边坡稳定性测试方法 Download PDFInfo
- Publication number
- CN110222405B CN110222405B CN201910460844.7A CN201910460844A CN110222405B CN 110222405 B CN110222405 B CN 110222405B CN 201910460844 A CN201910460844 A CN 201910460844A CN 110222405 B CN110222405 B CN 110222405B
- Authority
- CN
- China
- Prior art keywords
- formula
- reduction
- intensity
- strength
- fos
- 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.)
- Active
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/32—Circuit design at the digital level
- G06F30/333—Design for testability [DFT], e.g. scan chain or built-in self-test [BIST]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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/06—Power analysis or power optimisation
-
- 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)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Mining & Mineral Resources (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Structural Engineering (AREA)
- Operations Research (AREA)
- Paleontology (AREA)
- Computing Systems (AREA)
- Civil Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种基于降低强度包线折减法的岩质边坡稳定性测试方法,包括以下步骤:步骤一、根据Hoek–Brown准则,建立折减后强度与折减前强度的含参量数学关系,将折减后强度的确定问题转换为参量的确定问题,获得折减后强度的表达式;步骤二、推导参量和折减后强度的数学表达式,确定表达式当中的折减后强度;步骤三、分析参量可解范围,优化折减后强度的表达式来进行岩质边坡稳定性测试。本发明通过数学推导得到折减后强度的直接表达式,实际应用时能够在数值软件内完成运算,无需再导出至统计软件执行拟合运算,从而使得计算难度以及工作量得到大幅降低。
Description
技术领域
本发明属于岩质边坡稳定性测试领域,具体涉及一种基于降低强度包线折减法的岩质边坡稳定性测试方法,通过折减强度的直接表达式完成降低强度包络线运算。
背景技术
岩体强度通常遵循Hoek–Brown准则,强度折减法是计算岩质边坡稳定性的重要工具。到目前为止,对适用于Hoek–Brown准则的强度折减法,国内外学者已经开展了一定程度的研究,大致可以分为等效折减法、局部线性化折减法、降低强度包线折减法和参数折减法四类。其中,降低强度包线折减法物理含义最清晰,但计算繁杂不利应用。
Hoek–Brown强度准则可表示为:
式中,σci为完整岩石的单轴抗压强度;σ1表示破坏时的最大主应力;σ3表示破坏时的最小主应力。s、a由地质强度指标GSI决定;mb由地质强度指标GSI和材料常数mi共同决定:
式中,mi值反映岩石的软硬程度,且仅与岩石类型有关。地质强度指标GSI值与风化情况、岩石结构有关。D值反映扰动程度。降低强度包线折减法由Hammah提出,见文献《Theshear strength reduction method for the generalized Hoek-Brown criterion》。
降低强度包线折减法建立在Balmer研究成果的基础上。根据Balmer的研究,主应力包络线上某一点相对应的法向应力σn和剪应力τ可由以下方程确定:
对于Hoek–Brown准则,σn、τ和σ1、σ3之间的关系可以写为:
然后,按照以下步骤进行强度折减:
(1)绘制τ-σn剪应力强度包络线(图1中的线A),根据折减系数FOS对剪应力τ折减,得到折减后的强度包络线(图1中的线B);(2)根据Hoek–Brown准则对折减后的强度包络线进行拟合,得到折减后强度;(3)采用折减后强度进行稳定性分析。
得到折减后剪应力τred的表达式:
式中,σci red、mb red、sred和ared是折减后强度包络线的强度参数,即,折减后强度。
为确定折减后强度,Hammah等人提出一种近似估计的方法,其步骤如下:(1)确定围压上限σ3max、下限σ3min。将σ3max、σ3min代入式(25),计算σnmax、σnmin;(2)确定折减后强度包络线(图1中的线B)下方区域的面积(σnmax、σnmin之间,横坐标轴以上,强度包络线以下的区域);(3)将σci red、mb red、sred和ared作为自变量,根据式(25)计算σn,根据式(27)计算τred。采用最小二乘法等拟合方法,得到折减后强度包络线的拟合曲线(图1中的线C),使拟合曲线与折减后强度包络线下方区域的面积相等。
Hammah提出的降低强度包线折减法虽理论上可行,但实际应用过于繁琐。每一个迭代步都需要导出数据到统计软件(如OriginPro、R、SPSS等软件),利用最小二乘等拟合方法计算折减后强度,再将折减后强度导入数值分析软件(如FLAC3D、ANSYS等软件)试算,如此反复试算直至模型稳定性在某一强度状态下恰好达到极限平衡状态。
发明内容
本发明的目的在于针对上述现有技术中的问题,提供一种基于降低强度包线折减法的岩质边坡稳定性测试方法,能够在完成降低强度包络线运算的基础上大幅降低工作量。
为了实现上述目的,本发明采用的技术方案包括以下步骤:
步骤一、根据Hoek–Brown准则,建立折减后强度与折减前强度的含参量数学关系,将折减后强度的确定问题转换为参量的确定问题,获得折减后强度的表达式;
步骤二、推导参量和折减后强度的数学表达式,确定表达式当中的折减后强度;
步骤三、分析参量可解范围,优化折减后强度的表达式来进行岩质边坡稳定性测试。
所述的步骤一根据Hoek–Brown准则,对于GSI>30的非破碎岩体,a取0.5得到:
式中,GSI为地质强度指标;GSIr为折减后地质强度指标;mi为材料常数,物理含义为岩石劈裂强度与单轴无侧限抗压强度之比;s、mb表示折减前强度;sr、mb r表示折减后强度。
根据式(1)与式(2)之间的关系得到:
式中,参量k=(GSI-100)-1·(GSIr-100);
根据式(3)得到折减后强度的表达式为:
所述的步骤二推导参量和折减后强度的数学表达式的具体步骤为:
定义中间函数ψ1、ψ2,其中,ψ1为折减前强度的函数,ψ2为折减后强度的函数;
根据式(5)、式(6)的定义,ψ1、ψ2均为正值;ψ1已知,ψ2待解。
由式(5)得:
将式(4)代入式(6)得:
将式(7)代入式(8)得:
将式(9)整理为:
exp2(k)+ξ1exp(k)+ξ2=0 (10)
由式(10)得:
由式(11)知,k由ξ1、ξ2决定,由式(10)知,ξ1、ξ2由ψ1、ψ2决定;若ξ2=0,则k无解。
所述的ψ2通过如下方式进行推导求解:
根据Hammah方法,折减后剪应力τred的表达式为:
将式(5)、式(6)代入式(12)得到:
将式(13)整理为:
(ψ2+2)2-κ(ψ2+2)+κ=0 (14)
式中,κ=(ψ1+1)-1(ψ1+2)2FOS2;
由式(14)得到:
ψ2为正值,因此,将式(15)写为:
根据式(5)、(16)能够确定ψ1、ψ2,根据ψ1、ψ2能够确定式(10)中的ξ1、ξ2;根据ξ1、ξ2能够确定式(11)中的k,将k代入式(4)即能够得到折减后强度。
将折减后强度的表达式进行如下整理:
将ψ1写为κ和FOS的函数:
对比式(16)、(18)得到:
分析式(19)发现,当FOS>1时,ψ1<ψ2;当FOS<1时,ψ1>ψ2;当FOS=1时,ψ1=ψ2;
根据式(10)中ξ2的定义,发现:
结合式(11)、(20),k无解实际上对应于ξ2=0,亦即FOS=1的情况;若出现这种情况,可能导致程序运行中断甚至崩溃,因此采用将FOS求解范围限制为0<FOS<1和FOS>1。
与现有技术相比,本发明具有如下的有益效果:通过根据Hoek–Brown准则建立折减后强度与折减前强度的含参量数学关系,将折减后强度的确定问题转换为参量的确定问题,然后推导参量与折减后强度的数学表达式,最后分析参量取值范围,优化计算方法。本发明通过数学推导得到折减后强度的直接表达式,实际应用时能够在数值软件内完成运算,无需再导出至统计软件执行拟合运算,从而使得计算难度以及工作量得到大幅降低。
附图说明
图1降低强度包线折减法的原理示意图;
图2边坡计算图示;
图3安全系数计算结果统计图;
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
本发明基于降低强度包线折减法的岩质边坡稳定性测试方法,包括以下步骤:
步骤一、建立折减后强度与折减前强度的含参量数学关系,将折减后强度的确定问题转换为参量的确定问题,根据Hoek–Brown准则,对于非破碎岩体(GSI>30),a的取值范围介于0.5和0.52之间。假设a=0.5,得到:
发现式(1)与式(2)之间存在联系:
式中,k=(GSI-100)-1·(GSIr-100)。
根据式(3)得到折减后强度的表达式:
Hammah提出的降低强度包线折减法需要经过繁杂的拟合运算确定折减后强度,应用本发明式(4)可以通过基本数学运算大幅简化这一过程。式(4)的关键是,如何确定参量k。
步骤二、推导参量和折减后强度的数学表达式。
定义中间函数ψ1、ψ2:
根据式(5)、式(6)的定义,ψ1、ψ2均为正值。
由式(5)得:
将式(4)代入式(6)得:
将式(7)代入式(8)得:
式(9)可整理为:
exp2(k)+ξ1exp(k)+ξ2=0 (10)
由式(10)得:
由式(11)知,k由ξ1、ξ2决定。由式(10)知,ξ1、ξ2由ψ1、ψ2决定。根据式(5)、式(6)的定义,ψ1为折减前强度的函数,ψ2为折减后强度的函数。因此,ψ1已知,ψ2待解。
为确定ψ2,作如下推导:
根据Hammah方法,折减后剪应力τred的表达式为:
将式(5)、式(6)代入式(12)得:
式(13)可整理为:
(ψ2+2)2-κ(ψ2+2)+κ=0 (14)
式中,κ=(ψ1+1)-1(ψ1+2)2FOS2。
由式(14)得:
ψ2为正值,因此,式(15)可以写为:
根据式(5)、(16)可确定ψ1、ψ2。根据ψ1、ψ2可确定式(10)中的ξ1、ξ2。根据ξ1、ξ2可确定式(11)中的k。将k代入式(4)可得到折减后强度。
为方便应用,将折减后强度的表达式整理如下:
步骤三、分析参量取值范围,优化计算方法。
在式(10)中,ξ2由ψ1、ψ2决定,若ξ2=0,则k无解。为分析这种情况,需根据式(14)中κ的定义和式(5)中ψ1为正值的结论,将ψ1写为κ和FOS的函数:
对比式(16)、(18)得到:
分析式(19)发现,当FOS>1时,ψ1<ψ2;当FOS<1时,ψ1>ψ2;当FOS=1时,ψ1=ψ2。
根据式(10)中ξ2的定义,发现:
结合式(11)、(20)可知,k无解实际上对应于ξ2=0,亦即FOS=1的情况。若出现这种情况,可能导致程序运行中断甚至崩溃。一种可行的解决方法是,将FOS求解范围限制为0<FOS<1和FOS>1。
以上所述仅仅是本发明较佳的实施例,并不用以对本发明做任何形式上的限定,本领域技术人员应当理解的是,在不脱离本发明精神和原则的前提下,本发明的技术方案还可以进行若干简单的修改和替换,这些修改和替换也均会落入由所提交权利要求划定的范围之内。
Claims (2)
1.一种基于降低强度包线折减法的岩质边坡稳定性测试方法,其特征在于,包括:
步骤一、根据Hoek–Brown准则,建立折减后强度与折减前强度的含参量数学关系,将折减后强度的确定问题转换为参量的确定问题,获得折减后强度的表达式;
步骤二、推导参量和折减后强度的数学表达式,确定表达式当中的折减后强度;
步骤三、分析参量可解范围,优化折减后强度的表达式来进行岩质边坡稳定性测试;
所述的步骤一根据Hoek–Brown准则,对于GSI>30的非破碎岩体,a取0.5得到:
式中,GSI为地质强度指标;GSIr为折减后地质强度指标;mi为材料常数,物理含义为岩石劈裂强度与单轴无侧限抗压强度之比;s、mb表示折减前强度;sr、mb r表示折减后强度;
根据式(1)与式(2)之间的关系得到:
式中,参量k=(GSI-100)-1·(GSIr-100);
根据式(3)得到折减后强度的表达式为:
所述的步骤二推导参量和折减后强度的数学表达式的具体步骤为:
定义中间函数ψ1、ψ2,其中,ψ1为折减前强度的函数,ψ2为折减后强度的函数;
根据式(5)、式(6)的定义,ψ1、ψ2均为正值;ψ1已知,ψ2待解;
由式(5)得:
将式(4)代入式(6)得:
将式(7)代入式(8)得:
将式(9)整理为:
exp2(k)+ξ1exp(k)+ξ2=0 (10)
由式(10)得:
由式(11)知,k由ξ1、ξ2决定,由式(10)知,ξ1、ξ2由ψ1、ψ2决定;若ξ2=0,则k无解;
所述的ψ2通过如下方式进行推导求解:
根据Hammah方法,折减后剪应力τred的表达式为:
将式(5)、式(6)代入式(12)得到:
将式(13)整理为:
(ψ2+2)2-κ(ψ2+2)+κ=0 (14)
式中,κ=(ψ1+1)-1(ψ1+2)2FOS2;
由式(14)得到:
ψ2为正值,因此,将式(15)写为:
根据式(5)、(16)能够确定ψ1、ψ2,根据ψ1、ψ2能够确定式(10)中的ξ1、ξ2;根据ξ1、ξ2能够确定式(11)中的k,将k代入式(4)即能够得到折减后强度;
将折减后强度的表达式进行如下整理:
将ψ1写为κ和FOS的函数:
对比式(16)、(18)得到:
分析式(19)发现,当FOS>1时,ψ1<ψ2;当FOS<1时,ψ1>ψ2;当FOS=1时,ψ1=ψ2;
根据式(10)中ξ2的定义,发现:
2.根据权利要求1所述基于降低强度包线折减法的岩质边坡稳定性测试方法,其特征在于,结合式(11)、(20),k无解实际上对应于ξ2=0,亦即FOS=1的情况;若出现这种情况,将FOS求解范围限制为0<FOS<1和FOS>1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910460844.7A CN110222405B (zh) | 2019-05-30 | 2019-05-30 | 一种基于降低强度包线折减法的岩质边坡稳定性测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910460844.7A CN110222405B (zh) | 2019-05-30 | 2019-05-30 | 一种基于降低强度包线折减法的岩质边坡稳定性测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110222405A CN110222405A (zh) | 2019-09-10 |
CN110222405B true CN110222405B (zh) | 2022-10-04 |
Family
ID=67818587
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910460844.7A Active CN110222405B (zh) | 2019-05-30 | 2019-05-30 | 一种基于降低强度包线折减法的岩质边坡稳定性测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110222405B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111414576B (zh) * | 2020-04-01 | 2021-03-30 | 青岛农业大学 | 一种边坡安全系数不迭代求解方法 |
CN113405908B (zh) * | 2021-07-26 | 2022-10-21 | 中国建筑第五工程局有限公司 | 一种浅埋盾构隧道顶板稳定性判断方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103235889A (zh) * | 2013-05-02 | 2013-08-07 | 林杭 | 一种Hoek-Brown准则的直接强度折减方法 |
CN103310097A (zh) * | 2013-05-21 | 2013-09-18 | 河海大学 | 边坡稳定评价的改进双强度折减法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105335607B (zh) * | 2015-10-12 | 2017-06-16 | 湖北工业大学 | 一种边坡渐进破坏潜在滑动面的计算方法 |
-
2019
- 2019-05-30 CN CN201910460844.7A patent/CN110222405B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103235889A (zh) * | 2013-05-02 | 2013-08-07 | 林杭 | 一种Hoek-Brown准则的直接强度折减方法 |
CN103310097A (zh) * | 2013-05-21 | 2013-09-18 | 河海大学 | 边坡稳定评价的改进双强度折减法 |
Non-Patent Citations (2)
Title |
---|
基于Hoek-Brown准则的非等比强度折减方法;韩龙强等;《岩土力学》;20161010;全文 * |
基于Hoek-Brown强度折减法的边坡稳定性图表法研究;孙超伟等;《岩石力学与工程学报》;20180110(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110222405A (zh) | 2019-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110222405B (zh) | 一种基于降低强度包线折减法的岩质边坡稳定性测试方法 | |
Nagtegaal et al. | Some computational aspects of elastic‐plastic large strain analysis | |
CN111008502A (zh) | 一种数字孪生驱动的复杂装备故障预测方法 | |
CN107392786B (zh) | 基于支持向量机的矿井光纤光栅监测系统缺失数据补偿方法 | |
EP2390800A1 (en) | Sheet metal forming failure prediction using numerical simulations | |
WO2020125682A1 (zh) | 一种应用随钻数据计算岩石强度的方法和系统 | |
CN103616179A (zh) | 一种基于缺陷建模的传动齿轮疲劳寿命评估方法 | |
EP2402876A2 (en) | Method and apparatus for correlating simulation models with physical devices based on correlation metrics | |
CN114279842B (zh) | 一种确定岩石裂纹起裂应力和损伤应力的方法及系统 | |
CN105975589A (zh) | 一种高维数据的特征选择方法及装置 | |
CN113159162A (zh) | 一种基于信息融合与灰色关联度的故障诊断方法及系统 | |
CN111090957B (zh) | 一种高温结构危险点应力-应变计算方法 | |
CN109902326B (zh) | 一种有限元仿真实验效果测评方法 | |
CN106067032A (zh) | 用于数据异常监测的传感器选择方法的评价方法 | |
CN109885875B (zh) | 一种涡轮叶片概率损伤容限分析方法 | |
CN101701935A (zh) | 耐火砖内部质量在线检测装置 | |
Zhao et al. | Weight functions and stress intensity factors for pin‐loaded single‐edge notch bend specimen | |
CN106647274A (zh) | 一种连续生产过程中运行工况稳态判别方法 | |
CN116522096A (zh) | 基于动作捕捉的三维数字孪生内容智能制作方法 | |
CN102914970B (zh) | 工业过程控制系统中不能直接测量的性能参数的在线数据驱动估计方法 | |
CN112466416B (zh) | 一种结合镍基合金先验知识的材料数据清洗方法 | |
CN115310705A (zh) | 瓦斯涌出量的确定方法及其装置、计算机可读存储介质 | |
CN114418314A (zh) | 一种拉链式压裂施工时效性分析细节多模型求解方法 | |
CN112906130A (zh) | 一种基于小样本数据的结构可靠性评估方法 | |
CN111241753A (zh) | 一种基于bp神经网络的硬岩tbm掘进速率预测方法 |
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 |