CN112541257B - 一种基于溶蚀速率的灰岩边坡服役寿命的计算方法 - Google Patents

一种基于溶蚀速率的灰岩边坡服役寿命的计算方法 Download PDF

Info

Publication number
CN112541257B
CN112541257B CN202011422469.6A CN202011422469A CN112541257B CN 112541257 B CN112541257 B CN 112541257B CN 202011422469 A CN202011422469 A CN 202011422469A CN 112541257 B CN112541257 B CN 112541257B
Authority
CN
China
Prior art keywords
limestone slope
section
limestone
slope
structural plane
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
Application number
CN202011422469.6A
Other languages
English (en)
Other versions
CN112541257A (zh
Inventor
刘文连
李泽
眭素刚
王帮团
许汉华
唐果
李永红
樊亚红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Kunming Prospecting Design Institute of China Nonferrous Metals Industry Co Ltd
Original Assignee
Kunming University of Science and Technology
Kunming Prospecting Design Institute of China Nonferrous Metals Industry Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kunming University of Science and Technology, Kunming Prospecting Design Institute of China Nonferrous Metals Industry Co Ltd filed Critical Kunming University of Science and Technology
Priority to CN202011422469.6A priority Critical patent/CN112541257B/zh
Publication of CN112541257A publication Critical patent/CN112541257A/zh
Application granted granted Critical
Publication of CN112541257B publication Critical patent/CN112541257B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • 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

本发明公开了一种基于溶蚀速率的灰岩边坡服役寿命的计算方法，包括拟定灰岩边坡服役寿命计算的基本参数；进行灰岩边坡滑体的受力分析；根据溶蚀速率计算服役n年后灰岩边坡结构面贯通段和非贯通段的有效长度；根据结构面的贯通段的抗剪参数的溶蚀衰减速率计算服役n年后的结构面的贯通段的抗剪参数；建立灰岩边坡的长期赋役的稳定性非线性数学规划模型；求解灰岩边坡的服役寿命。本发明以灰岩边坡的岩体溶蚀速率计算服役期的贯通结构面的有效长度，以抗剪参数溶蚀衰减速率计算服役期的抗剪参数，基于目标函数、平衡方程和屈服条件建立了灰岩边坡的长期赋役的稳定性非线性数学规划模型，解决了预测灰岩边坡服役寿命数学建模的难题。

Description

一种基于溶蚀速率的灰岩边坡服役寿命的计算方法
技术领域
本发明涉及一种基于溶蚀速率的灰岩边坡服役寿命的计算方法，属于边坡稳定性分析技术领域。
背景技术
灰岩是石灰岩的简称，其是以方解石为主要成分的碳酸盐岩，按成因分类属于沉积岩。灰岩的主要化学成分是CaCO3，非常易溶蚀。灰岩边坡在我国西南地区广泛分布，灰岩边坡中的结构面在受水、大气的作用下会发生溶蚀，结构面的长度随着时间的推移其有效长度会增加，结构面的抗剪参数也会随着溶蚀而发生衰减。在灰岩边坡的服役过程中，由于灰岩岩体的溶蚀，结构面的长度、强度均发生变化，灰岩边坡的长期安全性会逐年下降。在工程中灰岩边坡在服役一定年限以后经常发生失稳破坏，给人民生命财产带来重大损失。
灰岩边坡服役寿命问题是一个复杂的力学和数学问题，要准确预测灰岩边坡的寿命，需要解决以下关键问题：(1)灰岩边坡结构面的长度和抗剪强度受岩体的溶蚀速度影响，需要在边坡稳定性模型中定量描述结构面长度和抗剪强度与岩体溶蚀速度的关系；(2)灰岩边坡的服役寿命与溶蚀速率直接相关，两者的关系呈现出高度非线性，力学建模较为复杂；(3)求解灰岩边坡的服役寿命是一个数学极大值问题，需要建立灰岩边坡服役寿命的数学规划模型。
在当前，工程界还没有完善的基于溶蚀速率的灰岩边坡服役寿命的计算方法，亟需建立求解灰岩边坡服役寿命的理论方法和数学模型，以能够定量、精确评价的灰岩边坡的长期安全性，为灰岩边坡的设计提供理论指导。
发明内容
本发明提供了一种基于溶蚀速率的灰岩边坡服役寿命的计算方法，以获得灰岩边坡的服役寿命，为评价灰岩边坡长期服役安全性的计算提供一种新方法。
本发明的技术方案是：一种基于溶蚀速率的灰岩边坡服役寿命的计算方法，所述方法步骤如下：
步骤1、拟定灰岩边坡服役寿命计算的基本参数；
步骤2、进行灰岩边坡滑体的受力分析；
步骤3、根据溶蚀速率计算服役n年后灰岩边坡结构面贯通段和非贯通段的有效长度；
步骤4、根据结构面的贯通段的抗剪参数的溶蚀衰减速率计算服役n年后的结构面的贯通段的抗剪参数；
步骤5、建立灰岩边坡的长期赋役的稳定性非线性数学规划模型；
步骤6、求解灰岩边坡的长期赋役的稳定性非线性数学规划模型，获得灰岩边坡的服役寿命。
所述拟定灰岩边坡服役寿命计算的基本参数，包括：
一、确定灰岩边坡的几何参数；
二、确定灰岩边坡的材料参数；
三、确定灰岩边坡岩体的溶蚀速率δD；
四、确定灰岩边坡中的结构面贯通段的抗剪参数溶蚀衰减速率。
所述灰岩边坡的几何参数，具体为：灰岩边坡的高度h，灰岩边坡滑体的顶部宽度b，灰岩边坡中结构面AB的倾角θ，灰岩边坡中结构面AB的贯通率k，灰岩边坡中结构面AB的长度L，灰岩边坡中结构面AB由贯通段BD、非贯通段AD组成；所述灰岩边坡的材料参数，具体为：灰岩边坡岩体的容重γr、灰岩边坡中结构面AB的贯通段BD的初始内摩擦角
Figure BDA0002823070540000021
灰岩边坡中结构面AB的贯通段BD的初始凝聚力cj；灰岩边坡岩体的初始内摩擦角
Figure BDA0002823070540000022
灰岩边坡岩体的初始凝聚力cr；所述灰岩边坡中的结构面贯通段的抗剪参数溶蚀衰减速率，具体为：灰岩边坡中结构面AB的贯通段 BD的内摩擦角溶蚀衰减速率
Figure BDA0002823070540000023
灰岩边坡岩体中结构面AB的贯通段BD的凝聚力溶蚀衰减速率
Figure BDA0002823070540000024
所述灰岩边坡滑体的受力分析，具体为：灰岩边坡滑体的重心处作用的自重力Gr，灰岩边坡结构面AB的贯通段BD上作用有法向力N1，灰岩边坡结构面AB的贯通段BD上作用有切向力S1，灰岩边坡结构面AB的非贯通段AD上作用有法向力 N2，灰岩边坡结构面AB的非贯通段AD上作用有切向力S2。
所述根据溶蚀速率计算服役n年后灰岩边坡结构面贯通段和非贯通段的有效长度，具体为：
一、根据溶蚀速率计算服役n年后灰岩边坡结构面AB的贯通段BD的有效长度
Figure BDA0002823070540000031
Figure BDA0002823070540000032
式中：L是灰岩边坡中结构面AB的长度，k是灰岩边坡中结构面AB的贯通率，δD是灰岩边坡岩体的溶蚀速率，n是灰岩边坡的服役寿命；
二、根据溶蚀速率计算服役n年后灰岩边坡结构面AB的非贯通段AD的有效长度
Figure BDA0002823070540000033
Figure BDA0002823070540000034
所述根据抗剪参数溶蚀衰减速率计算服役n年后灰岩边坡结构面的抗剪参数，具体为：
一、根据结构面的贯通段的内摩擦角溶蚀衰减速率计算服役n年后的结构面的贯通段的内摩擦角
Figure BDA0002823070540000035
具体为：
Figure BDA0002823070540000036
式中：
Figure BDA0002823070540000037
是灰岩边坡中的结构面AB的贯通段BD的初始内摩擦角，
Figure BDA0002823070540000038
是灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率，n是灰岩边坡的服役寿命；
二、根据结构面的贯通段的凝聚力溶蚀衰减速率计算服役n年后的结构面的贯通段的凝聚力
Figure BDA0002823070540000039
具体为：
Figure BDA00028230705400000310
式中：cj是灰岩边坡中结构面AB的贯通段BD的初始凝聚力，
Figure BDA00028230705400000311
是灰岩边坡中结构面AB的贯通段BD的凝聚力溶蚀衰减速率。
所述建立灰岩边坡的长期赋役的稳定性非线性数学规划模型，具体为：
一、建立目标函数
设灰岩边坡的服役寿命为目标函数，并求其最大值，目标函数如下：
Maximize:n
式中：n是灰岩边坡的服役寿命；Maximize表示“使最大”；
二、建立服役n年后灰岩边坡滑体的平衡方程，具体为：
①建立服役n年后灰岩边坡滑体的水平方向平衡方程，具体为：
N1 sinθ-S1 cosθ+N2 sinθ-S2 cosθ＝0
式中：N1是灰岩边坡结构面AB的贯通段BD上作用的法向力，S1是灰岩边坡结构面AB的贯通段BD上作用的切向力，N2是灰岩边坡结构面AB的非贯通段AD上作用的法向力，S2是灰岩边坡结构面AB的非贯通段AD上作用的切向力，θ是灰岩边坡结构面AB的倾角；
②建立服役n年后灰岩边坡滑体的竖直方向平衡方程，具体为：
N1 cosθ+S1 sinθ+N2 cosθ+S2 sinθ-Gr＝0
式中：Gr是灰岩边坡滑体的重心处作用的自重力，Gr＝γrhb/2，h是灰岩边坡的高度，b灰岩边坡滑体的顶部的宽度，γr是灰岩边坡岩体的容重；
③建立服役n年后灰岩边坡结构面AB的贯通段BD的屈服条件，具体为：
Figure BDA0002823070540000041
式中：
Figure BDA0002823070540000042
是服役n年后灰岩边坡结构面AB的贯通段BD的有效长度，
Figure BDA0002823070540000043
是服役n年后的结构面AB的贯通段BD的内摩擦角，
Figure BDA0002823070540000044
是服役n年后的结构面AB的贯通段BD 的凝聚力；
④建立服役n年后灰岩边坡结构面AB的非贯通段AD的屈服条件，具体为：
Figure BDA0002823070540000045
式中：
Figure BDA0002823070540000046
是服役n年后灰岩边坡结构面AB的非贯通段AD的有效长度，
Figure BDA0002823070540000047
是灰岩边坡岩体的初始内摩擦角，cr是灰岩边坡岩体的初始凝聚力；
⑤建立附加约束条件，具体为：
Figure BDA0002823070540000048
⑥建立灰岩边坡的长期赋役的稳定性非线性数学规划模型
将目标函数、灰岩边坡滑体的平衡方程、贯通段和非贯通段的屈服条件以及附加约束条件集成，得到灰岩边坡的长期赋役的稳定性模型如下：
Figure BDA0002823070540000051
式中：L是灰岩边坡中结构面AB的长度，k是灰岩边坡中结构面AB的贯通率，δD是灰岩边坡岩体的溶蚀速率，
Figure BDA0002823070540000052
是灰岩边坡中的结构面AB的贯通段BD的初始内摩擦角，
Figure BDA0002823070540000053
是灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率，cj是灰岩边坡中结构面AB的贯通段BD的初始凝聚力，
Figure BDA0002823070540000054
是灰岩边坡中结构面AB的贯通段 BD的凝聚力溶蚀衰减速率。
所述求解灰岩边坡的长期赋役的稳定性非线性数学规划模型，具体为：将已知参数L、k、δD、
Figure BDA0002823070540000055
cj、
Figure BDA0002823070540000056
θ代入灰岩边坡的长期赋役的稳定性非线性数学规划模型，以灰岩边坡的服役寿命n为目标函数，以
Figure BDA0002823070540000057
N1、S1、 N2、S2、Gr为决策变量，并使用“序列二次规划”求解该非线性数学规划模型，求解得到灰岩边坡的服役寿命n，以及决策变量
Figure BDA0002823070540000058
N1、S1、N2、S2、 Gr的计算结果；其中，L为灰岩边坡中结构面AB的长度，k为灰岩边坡中结构面AB 的贯通率，δD是灰岩边坡岩体的溶蚀速率，
Figure BDA0002823070540000059
是灰岩边坡中的结构面AB的贯通段 BD的初始内摩擦角，cj是灰岩边坡中结构面AB的贯通段BD的初始凝聚力，
Figure BDA00028230705400000510
是灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率，
Figure BDA00028230705400000511
是灰岩边坡中结构面AB的贯通段BD的凝聚力溶蚀衰减速率，θ是灰岩边坡结构面AB的倾角；
Figure BDA00028230705400000512
是服役n年后灰岩边坡结构面AB的贯通段BD的有效长度，
Figure BDA00028230705400000513
是服役n年后灰岩边坡结构面AB的非贯通段AD的有效长度，
Figure BDA00028230705400000514
是服役n年后的结构面AB的贯通段BD的内摩擦角，
Figure BDA0002823070540000061
是服役n年后的结构面AB的贯通段BD的凝聚力，N1是灰岩边坡结构面AB的贯通段BD上作用的法向力，S1是灰岩边坡结构面AB的贯通段BD上作用的切向力，N2是灰岩边坡结构面AB的非贯通段AD上作用的法向力，S2是灰岩边坡结构面AB的非贯通段AD上作用的切向力，Gr是灰岩边坡滑体的重心处作用的自重力。
本发明的有益效果是：
(1)本发明以灰岩边坡的岩体溶蚀速率计算服役期的贯通结构面的有效长度，以抗剪参数溶蚀衰减速率计算服役期的抗剪参数，基于目标函数、平衡方程和屈服条件建立了灰岩边坡的长期赋役的稳定性非线性数学规划模型，解决了预测灰岩边坡服役寿命数学建模的难题。
(2)本发明方法理论严谨、计算简单，获得的计算结果精度较高。
(3)本发明实现了基于岩体溶蚀速率和抗剪参数溶蚀衰减速率计算边坡的服役寿命，为灰岩边坡的长期安全性评价提供了理论依据。
附图说明
图1为本发明技术流程图；
图2为灰岩边坡计算模型示意图；
图3为灰岩边坡滑体的受力状态示意图。
具体实施方式
实施例1：如图1-图3所示，一种基于溶蚀速率的灰岩边坡服役寿命的计算方法，所述方法步骤如下：
步骤1、拟定灰岩边坡服役寿命计算的基本参数；
步骤2、进行灰岩边坡滑体的受力分析；
步骤3、根据溶蚀速率计算服役n年后灰岩边坡结构面贯通段和非贯通段的有效长度；
步骤4、根据结构面的贯通段的抗剪参数的溶蚀衰减速率计算服役n年后的结构面的贯通段的抗剪参数；
步骤5、建立灰岩边坡的长期赋役的稳定性非线性数学规划模型；
步骤6、求解灰岩边坡的长期赋役的稳定性非线性数学规划模型，获得灰岩边坡的服役寿命。
进一步地，可以设置所述拟定灰岩边坡服役寿命计算的基本参数，包括：一、确定灰岩边坡的几何参数；二、确定灰岩边坡的材料参数；三、确定灰岩边坡岩体的溶蚀速率δD；四、确定灰岩边坡中的结构面贯通段的抗剪参数溶蚀衰减速率。
进一步，可以设置所述灰岩边坡的几何参数，具体为：灰岩边坡的高度h，灰岩边坡滑体的顶部宽度b，灰岩边坡中结构面AB的倾角θ，灰岩边坡中结构面AB的贯通率k，灰岩边坡中结构面AB的长度L，灰岩边坡中结构面AB由贯通段BD、非贯通段AD组成；所述灰岩边坡的材料参数，具体为：灰岩边坡岩体的容重γr、灰岩边坡中结构面AB的贯通段BD的初始内摩擦角
Figure BDA0002823070540000071
灰岩边坡中结构面AB的贯通段 BD的初始凝聚力cj；灰岩边坡岩体的初始内摩擦角
Figure BDA0002823070540000072
灰岩边坡岩体的初始凝聚力cr；所述灰岩边坡中的结构面贯通段的抗剪参数溶蚀衰减速率，具体为：灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率
Figure BDA0002823070540000073
灰岩边坡岩体中结构面 AB的贯通段BD的凝聚力溶蚀衰减速率
Figure BDA0002823070540000074
本实施例采用如下数据：灰岩边坡的高度h取25.00m，灰岩边坡滑体的顶部宽度b取10.00m，结构面AB的倾角θ取53°，灰岩边坡中结构面AB的贯通率k取0.85，灰岩边坡中结构面AB的长度L取31.30m，灰岩边坡中结构面AB的由贯通段 BD、非贯通段AD组成；灰岩边坡岩体的容重γr取24.00kN/m3、灰岩边坡中的结构面AB的贯通段BD的初始内摩擦角
Figure BDA0002823070540000075
取32°、灰岩边坡中结构面AB的贯通段BD 的初始凝聚力cj取40.00kPa；灰岩边坡岩体的初始内摩擦角
Figure BDA0002823070540000076
取40.00°、灰岩边坡岩体的初始凝聚力cr取500.00kPa；δD取为取3.0％；灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率
Figure BDA0002823070540000077
取0.5％，灰岩边坡岩体中结构面AB的贯通段 BD的凝聚力溶蚀衰减速率
Figure BDA0002823070540000078
取0.6％。
进一步地，所述灰岩边坡滑体的受力分析，具体为：灰岩边坡滑体的受力情况如图2所示，灰岩边坡滑体的重心处作用的自重力Gr，灰岩边坡结构面AB的贯通段BD上作用有法向力N1，灰岩边坡结构面AB的贯通段BD上作用有切向力S1，灰岩边坡结构面AB的非贯通段AD上作用有法向力N2，灰岩边坡结构面AB的非贯通段AD上作用有切向力S2。
进一步地，所述根据溶蚀速率计算服役n年后灰岩边坡结构面贯通段和非贯通段的有效长度，具体为：
一、根据溶蚀速率计算服役n年后灰岩边坡结构面AB的贯通段BD的有效长度
Figure BDA0002823070540000081
Figure BDA0002823070540000082
式中：L是灰岩边坡中结构面AB的长度，k是灰岩边坡中结构面AB的贯通率，δD是灰岩边坡岩体的溶蚀速率，n是灰岩边坡的服役寿命；
二、根据溶蚀速率计算服役n年后灰岩边坡结构面AB的非贯通段AD的有效长度
Figure BDA0002823070540000083
Figure BDA0002823070540000084
式中：
Figure BDA0002823070540000085
是服役n年后灰岩边坡结构面AB的非贯通段AD 的有效长度，L是灰岩边坡中结构面AB的长度，k是灰岩边坡中结构面AB的贯通率，δD是灰岩边坡岩体的溶蚀速率，n是灰岩边坡的服役寿命。
进一步地，所述根据抗剪参数溶蚀衰减速率计算服役n年后灰岩边坡结构面的抗剪参数，具体为：
一、根据结构面的贯通段的内摩擦角溶蚀衰减速率计算服役n年后的结构面的贯通段的内摩擦角
Figure BDA0002823070540000086
具体为：
Figure BDA0002823070540000087
式中：
Figure BDA0002823070540000088
是服役n年后的结构面AB的贯通段BD的内摩擦角，
Figure BDA0002823070540000089
是灰岩边坡中的结构面AB的贯通段BD的初始内摩擦角，
Figure BDA00028230705400000810
是灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率，n是灰岩边坡的服役寿命；
二、根据结构面的贯通段的凝聚力溶蚀衰减速率计算服役n年后的结构面的贯通段的凝聚力
Figure BDA00028230705400000811
具体为：
Figure BDA00028230705400000812
式中：
Figure BDA00028230705400000813
服役n年后的结构面AB的贯通段 BD的凝聚力，cj是灰岩边坡中结构面AB的贯通段BD的初始凝聚力，
Figure BDA00028230705400000814
是灰岩边坡中结构面AB的贯通段BD的凝聚力溶蚀衰减速率，n是灰岩边坡的服役寿命。
进一步地，所述建立灰岩边坡的长期赋役的稳定性非线性数学规划模型，具体为：
一、建立目标函数
设灰岩边坡的服役寿命为目标函数，并求其最大值，目标函数如下：
Maximize:n
式中：n是灰岩边坡的服役寿命；Maximize表示“使最大”；
二、建立服役n年后灰岩边坡滑体的平衡方程，具体为：
①建立服役n年后灰岩边坡滑体的水平方向平衡方程，具体为：
N1 sinθ-S1 cosθ+N2 sinθ-S2 cosθ＝0
式中：N1是灰岩边坡结构面AB的贯通段BD上作用的法向力，S1是灰岩边坡结构面AB的贯通段BD上作用的切向力，N2是灰岩边坡结构面AB的非贯通段AD上作用的法向力，S2是灰岩边坡结构面AB的非贯通段AD上作用的切向力，θ是灰岩边坡结构面AB的倾角；
②建立服役n年后灰岩边坡滑体的竖直方向平衡方程，具体为：
N1 cosθ+S1 sinθ+N2 cosθ+S2 sinθ-Gr＝0
式中：N1是灰岩边坡结构面AB的贯通段BD上作用的法向力，S1是灰岩边坡结构面AB的贯通段BD上作用的切向力，N2是灰岩边坡结构面AB的非贯通段AD上作用的法向力，S2是灰岩边坡结构面AB的非贯通段AD上作用的切向力，θ是灰岩边坡结构面AB的倾角，Gr是灰岩边坡滑体的重心处作用的自重力，Gr＝γrhb/2，h是灰岩边坡的高度，b灰岩边坡滑体的顶部的宽度，γr是灰岩边坡岩体的容重；
③建立服役n年后灰岩边坡结构面AB的贯通段BD的屈服条件，具体为：
Figure BDA0002823070540000091
式中：N1是灰岩边坡结构面AB的贯通段BD上作用的法向力，S1是灰岩边坡结构面AB的贯通段BD上作用的切向力，
Figure BDA0002823070540000092
是服役n年后灰岩边坡结构面AB的贯通段BD 的有效长度，
Figure BDA0002823070540000093
是服役n年后的结构面AB的贯通段BD的内摩擦角，
Figure BDA0002823070540000094
是服役n年后的结构面AB的贯通段BD的凝聚力；
④建立服役n年后灰岩边坡结构面AB的非贯通段AD的屈服条件，具体为：
Figure BDA0002823070540000095
式中：N2是灰岩边坡结构面AB的非贯通段AD上作用的法向力，S2是灰岩边坡结构面AB的非贯通段AD上作用的切向力，
Figure BDA0002823070540000096
是服役n年后灰岩边坡结构面AB的非贯通段AD的有效长度，
Figure BDA0002823070540000101
是灰岩边坡岩体的初始内摩擦角，cr是灰岩边坡岩体的初始凝聚力；
⑤建立附加约束条件，具体为：
Figure BDA0002823070540000102
⑥建立灰岩边坡的长期赋役的稳定性非线性数学规划模型
将目标函数、灰岩边坡滑体的平衡方程、贯通段和非贯通段的屈服条件以及附加约束条件集成，得到灰岩边坡的长期赋役的稳定性模型如下：
Figure BDA0002823070540000103
式中：L是灰岩边坡中结构面AB的长度，k是灰岩边坡中结构面AB的贯通率，δD是灰岩边坡岩体的溶蚀速率，
Figure BDA0002823070540000104
是灰岩边坡中的结构面AB的贯通段BD的初始内摩擦角，
Figure BDA0002823070540000105
是灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率，cj是灰岩边坡中结构面AB的贯通段BD的初始凝聚力，
Figure BDA0002823070540000106
是灰岩边坡中结构面AB的贯通段 BD的凝聚力溶蚀衰减速率。
进一步地，所述求解灰岩边坡的长期赋役的稳定性非线性数学规划模型，具体为：将已知参数L、k、δD、
Figure BDA0002823070540000107
cj、
Figure BDA0002823070540000108
θ代入灰岩边坡的长期赋役的稳定性非线性数学规划模型，以灰岩边坡的服役寿命n为目标函数，以
Figure BDA0002823070540000109
N1、S1、N2、S2、Gr为决策变量，并使用“序列二次规划”求解该非线性数学规划模型，求解得到灰岩边坡的服役寿命n为45.32年，以及决策变量
Figure BDA00028230705400001010
N1、S1、N2、S2、Gr的计算结果如表1所示；其中，L为灰岩边坡中结构面AB的长度，k为灰岩边坡中结构面AB的贯通率，δD是灰岩边坡岩体的溶蚀速率，
Figure BDA0002823070540000111
是灰岩边坡中的结构面AB的贯通段BD的初始内摩擦角，cj是灰岩边坡中结构面AB的贯通段BD的初始凝聚力，
Figure BDA0002823070540000112
是灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率，
Figure BDA0002823070540000113
是灰岩边坡中结构面AB的贯通段BD的凝聚力溶蚀衰减速率，θ是灰岩边坡结构面AB的倾角；
Figure BDA0002823070540000114
是服役n年后灰岩边坡结构面AB的贯通段BD的有效长度，
Figure BDA0002823070540000115
是服役n年后灰岩边坡结构面AB的非贯通段AD的有效长度，
Figure BDA0002823070540000116
是服役 n年后的结构面AB的贯通段BD的内摩擦角，
Figure 1
是服役n年后的结构面AB的贯通段 BD的凝聚力，N1是灰岩边坡结构面AB的贯通段BD上作用的法向力，S1是灰岩边坡结构面AB的贯通段BD上作用的切向力，N2是灰岩边坡结构面AB的非贯通段 AD上作用的法向力，S2是灰岩边坡结构面AB的非贯通段AD上作用的切向力，Gr是灰岩边坡滑体的重心处作用的自重力。
表1实施例计算结果统计表
Figure BDA0002823070540000118
本发明的基本原理是：本发明以灰岩边坡为研究对象，进行灰岩边坡滑体的受力分析；根据溶蚀速率计算服役n年后灰岩边坡结构面贯通段和非贯通段的有效长度；根据抗剪参数溶蚀衰减速率计算服役n年后灰岩边坡结构面的抗剪参数；将灰岩边坡的服役寿命作为目标函数，结合灰岩边坡滑体的平衡方程、结构面贯通段和非贯通段的屈服条件，建立灰岩边坡的长期赋役的稳定性非线性数学规划模型，使用“序列二次规划法”求解数学规划模型，获得灰岩边坡的服役寿命。
上面结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (4)

1.一种基于溶蚀速率的灰岩边坡服役寿命的计算方法，其特征在于：所述方法步骤如下：
步骤1、拟定灰岩边坡服役寿命计算的基本参数；
步骤2、进行灰岩边坡滑体的受力分析；
步骤3、根据溶蚀速率计算服役n年后灰岩边坡结构面贯通段和非贯通段的有效长度；
步骤4、根据结构面的贯通段的抗剪参数的溶蚀衰减速率计算服役n年后的结构面的贯通段的抗剪参数；
步骤5、建立灰岩边坡的长期赋役的稳定性非线性数学规划模型；
步骤6、求解灰岩边坡的长期赋役的稳定性非线性数学规划模型，获得灰岩边坡的服役寿命；
所述根据溶蚀速率计算服役n年后灰岩边坡结构面贯通段和非贯通段的有效长度，具体为：
一、根据溶蚀速率计算服役n年后灰岩边坡结构面AB的贯通段BD的有效长度
Figure FDA0003986999900000011
Figure FDA0003986999900000012
式中：L是灰岩边坡中结构面AB的长度，k是灰岩边坡中结构面AB的贯通率，δD是灰岩边坡岩体的溶蚀速率，n是灰岩边坡的服役寿命；
二、根据溶蚀速率计算服役n年后灰岩边坡结构面AB的非贯通段AD的有效长度
Figure FDA0003986999900000013
Figure FDA0003986999900000014
所述根据抗剪参数溶蚀衰减速率计算服役n年后灰岩边坡结构面的抗剪参数，具体为：
一、根据结构面的贯通段的内摩擦角溶蚀衰减速率计算服役n年后的结构面的贯通段的内摩擦角
Figure FDA0003986999900000015
具体为：
Figure FDA0003986999900000016
式中：
Figure FDA0003986999900000021
是灰岩边坡中的结构面AB的贯通段BD的初始内摩擦角，
Figure FDA0003986999900000022
是灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率，n是灰岩边坡的服役寿命；
二、根据结构面的贯通段的凝聚力溶蚀衰减速率计算服役n年后的结构面的贯通段的凝聚力
Figure FDA0003986999900000023
具体为：
Figure FDA0003986999900000024
式中：cj是灰岩边坡中结构面AB的贯通段BD的初始凝聚力，
Figure FDA0003986999900000025
是灰岩边坡中结构面AB的贯通段BD的凝聚力溶蚀衰减速率；
所述建立灰岩边坡的长期赋役的稳定性非线性数学规划模型，具体为：
一、建立目标函数
设灰岩边坡的服役寿命为目标函数，并求其最大值，目标函数如下：
Maximize:n
式中：n是灰岩边坡的服役寿命；Maximize表示“使最大”；
二、建立服役n年后灰岩边坡滑体的平衡方程，具体为：
①建立服役n年后灰岩边坡滑体的水平方向平衡方程，具体为：
N1 sinθ-S1 cosθ+N2 sinθ-S2 cosθ＝0
式中：N1是灰岩边坡结构面AB的贯通段BD上作用的法向力，S1是灰岩边坡结构面AB的贯通段BD上作用的切向力，N2是灰岩边坡结构面AB的非贯通段AD上作用的法向力，S2是灰岩边坡结构面AB的非贯通段AD上作用的切向力，θ是灰岩边坡结构面AB的倾角；
②建立服役n年后灰岩边坡滑体的竖直方向平衡方程，具体为：
N1 cosθ+S1 sinθN2 cosθ+S2 sinθ-Gr＝0
式中：Gr是灰岩边坡滑体的重心处作用的自重力，Gr＝γrhb/2，h是灰岩边坡的高度，b灰岩边坡滑体的顶部的宽度，γr是灰岩边坡岩体的容重；
③建立服役n年后灰岩边坡结构面AB的贯通段BD的屈服条件，具体为：
Figure FDA0003986999900000026
式中：
Figure FDA0003986999900000027
是服役n年后灰岩边坡结构面AB的贯通段BD的有效长度，
Figure FDA0003986999900000028
是服役n年后的结构面AB的贯通段BD的内摩擦角，
Figure FDA0003986999900000031
是服役n年后的结构面AB的贯通段BD的凝聚力；
④建立服役n年后灰岩边坡结构面AB的非贯通段AD的屈服条件，具体为：
Figure FDA0003986999900000032
式中：
Figure FDA0003986999900000033
是服役n年后灰岩边坡结构面AB的非贯通段AD的有效长度，
Figure FDA0003986999900000034
是灰岩边坡岩体的初始内摩擦角，cr是灰岩边坡岩体的初始凝聚力；
⑤建立附加约束条件，具体为：
Figure FDA0003986999900000035
⑥建立灰岩边坡的长期赋役的稳定性非线性数学规划模型
将目标函数、灰岩边坡滑体的平衡方程、贯通段和非贯通段的屈服条件以及附加约束条件集成，得到灰岩边坡的长期赋役的稳定性模型如下：
Figure FDA0003986999900000036
式中：L是灰岩边坡中结构面AB的长度，k是灰岩边坡中结构面AB的贯通率，δD是灰岩边坡岩体的溶蚀速率，
Figure FDA0003986999900000037
是灰岩边坡中的结构面AB的贯通段BD的初始内摩擦角，
Figure FDA0003986999900000038
是灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率，cj是灰岩边坡中结构面AB的贯通段BD的初始凝聚力，
Figure FDA0003986999900000039
是灰岩边坡中结构面AB的贯通段BD的凝聚力溶蚀衰减速率；
所述求解灰岩边坡的长期赋役的稳定性非线性数学规划模型，具体为：将已知参数L、k、δD、
Figure FDA0003986999900000041
cj、
Figure FDA0003986999900000042
θ代入灰岩边坡的长期赋役的稳定性非线性数学规划模型，以灰岩边坡的服役寿命n为目标函数，以
Figure FDA0003986999900000043
N1、S1、N2、S2、Gr为决策变量，并使用“序列二次规划”求解该非线性数学规划模型，求解得到灰岩边坡的服役寿命n，以及决策变量
Figure FDA0003986999900000044
N1、S1、N2、S2、Gr的计算结果；其中，L为灰岩边坡中结构面AB的长度，k为灰岩边坡中结构面AB的贯通率，δD是灰岩边坡岩体的溶蚀速率，
Figure FDA0003986999900000045
是灰岩边坡中的结构面AB的贯通段BD的初始内摩擦角，cj是灰岩边坡中结构面AB的贯通段BD的初始凝聚力，
Figure FDA0003986999900000046
是灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率，
Figure FDA0003986999900000047
是灰岩边坡中结构面AB的贯通段BD的凝聚力溶蚀衰减速率，θ是灰岩边坡结构面AB的倾角；
Figure FDA0003986999900000048
是服役n年后灰岩边坡结构面AB的贯通段BD的有效长度，
Figure FDA0003986999900000049
是服役n年后灰岩边坡结构面AB的非贯通段AD的有效长度，
Figure FDA00039869999000000410
是服役n年后的结构面AB的贯通段BD的内摩擦角，
Figure FDA00039869999000000411
是服役n年后的结构面AB的贯通段BD的凝聚力，N1是灰岩边坡结构面AB的贯通段BD上作用的法向力，S1是灰岩边坡结构面AB的贯通段BD上作用的切向力，N2是灰岩边坡结构面AB的非贯通段AD上作用的法向力，S2是灰岩边坡结构面AB的非贯通段AD上作用的切向力，Gr是灰岩边坡滑体的重心处作用的自重力。
2.根据权利要求1所述的基于溶蚀速率的灰岩边坡服役寿命的计算方法，其特征在于：所述拟定灰岩边坡服役寿命计算的基本参数，包括：
一、确定灰岩边坡的几何参数；
二、确定灰岩边坡的材料参数；
三、确定灰岩边坡岩体的溶蚀速率δD；
四、确定灰岩边坡中的结构面贯通段的抗剪参数溶蚀衰减速率。
3.根据权利要求2所述的基于溶蚀速率的灰岩边坡服役寿命的计算方法，其特征在于：所述灰岩边坡的几何参数，具体为：灰岩边坡的高度h，灰岩边坡滑体的顶部宽度b，灰岩边坡中结构面AB的倾角θ，灰岩边坡中结构面AB的贯通率k，灰岩边坡中结构面AB的长度L，灰岩边坡中结构面AB由贯通段BD、非贯通段AD组成；所述灰岩边坡的材料参数，具体为：灰岩边坡岩体的容重γr、灰岩边坡中结构面AB的贯通段BD的初始内摩擦角
Figure FDA0003986999900000051
灰岩边坡中结构面AB的贯通段BD的初始凝聚力cj；灰岩边坡岩体的初始内摩擦角
Figure FDA0003986999900000052
灰岩边坡岩体的初始凝聚力cr；所述灰岩边坡中的结构面贯通段的抗剪参数溶蚀衰减速率，具体为：灰岩边坡中结构面AB的贯通段BD的内摩擦角溶蚀衰减速率
Figure FDA0003986999900000053
灰岩边坡岩体中结构面AB的贯通段BD的凝聚力溶蚀衰减速率
Figure FDA0003986999900000054
4.根据权利要求1所述的基于溶蚀速率的灰岩边坡服役寿命的计算方法，其特征在于：所述灰岩边坡滑体的受力分析，具体为：灰岩边坡滑体的重心处作用的自重力Gr，灰岩边坡结构面AB的贯通段BD上作用有法向力N1，灰岩边坡结构面AB的贯通段BD上作用有切向力S1，灰岩边坡结构面AB的非贯通段AD上作用有法向力N2，灰岩边坡结构面AB的非贯通段AD上作用有切向力S2。
CN202011422469.6A 2020-12-08 2020-12-08 一种基于溶蚀速率的灰岩边坡服役寿命的计算方法 Active CN112541257B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011422469.6A CN112541257B (zh) 2020-12-08 2020-12-08 一种基于溶蚀速率的灰岩边坡服役寿命的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011422469.6A CN112541257B (zh) 2020-12-08 2020-12-08 一种基于溶蚀速率的灰岩边坡服役寿命的计算方法

Publications (2)

Publication Number Publication Date
CN112541257A CN112541257A (zh) 2021-03-23
CN112541257B true CN112541257B (zh) 2023-02-10

Family

ID=75019196

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011422469.6A Active CN112541257B (zh) 2020-12-08 2020-12-08 一种基于溶蚀速率的灰岩边坡服役寿命的计算方法

Country Status (1)

Country Link
CN (1) CN112541257B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115791904B (zh) * 2023-01-09 2023-05-19 珠江水利委员会珠江水利科学研究院 一种砂浆溶蚀测试方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107330145A (zh) * 2017-05-26 2017-11-07 昆明理工大学 一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法
CN107330224A (zh) * 2017-07-24 2017-11-07 中国地质大学（武汉） 一种分条间作用力倾角非假定的边坡稳定性分析条分法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MX9200051A (es) * 1992-01-07 1993-07-01 Jose Luis Siller Franco Concector por friccion mejorado para anclar acero de refuerzo a tension en elementos de concreto preesforzado o reforzado.
CN107341285B (zh) * 2017-05-26 2020-08-25 昆明理工大学 一种地震力作用下节理岩质边坡倾倒破坏极限承载力分析下限法
CN108170898B (zh) * 2017-12-06 2021-12-17 昆明理工大学 一种节理岩质边坡可靠度分析下限法
CN108254782B (zh) * 2018-02-09 2019-11-05 中国地质大学（北京） 一种边坡地震破坏失稳概率的获取方法及系统
CN111750223A (zh) * 2020-06-22 2020-10-09 昆明理工大学 一种用于边坡稳定性分析的数据记录装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107330145A (zh) * 2017-05-26 2017-11-07 昆明理工大学 一种同时考虑岩体平动和转动效应的节理岩质边坡极限承载力分析下限法
CN107330224A (zh) * 2017-07-24 2017-11-07 中国地质大学（武汉） 一种分条间作用力倾角非假定的边坡稳定性分析条分法

Also Published As

Publication number Publication date
CN112541257A (zh) 2021-03-23

Similar Documents

Publication Publication Date Title
CN106855897B (zh) 适用于压力衰竭地层的井壁稳定的研究方法
CN103790568A (zh) 一种钻井参数与效率实时优化方法
CN112541257B (zh) 一种基于溶蚀速率的灰岩边坡服役寿命的计算方法
CN106295214B (zh) 一种矿井瓦斯浓度预警方法
CN105928813A (zh) 一种油井防砂筛管冲蚀寿命预测方法
O'Brien et al. A wind tunnel study of particle kinematics during crust rupture and erosion
CN108694254A (zh) 一种变产变压生产气井产量经验递减曲线分析方法
CN111475924B (zh) 一种降雨诱发变形的卸荷岩质边坡稳定性分析方法
US11781864B1 (en) Method and system for survey and stability evaluation of upstream tailing fill dam
CN110469299B (zh) 一种采油井注水开发见效效果评价方法
CN108867666B (zh) 基于开挖变形的结构面控制边坡稳定性评价方法
WO2019091134A1 (zh) 煤层顶板砂岩富水性等级预测方法
CN104405368A (zh) 高温高压气井井口测压适应性判别方法
CN116205064A (zh) 复杂地层隧道施工中盾构刀具服役期磨损量的预测方法
CN112182793A (zh) 一种气井防砂管冲蚀寿命预测的方法
CN106501147A (zh) 一种堤坝渗透稳定性与溃堤风险的测定方法
CN104790943B (zh) 油气储层含油性与孔隙性综合指数计算及储层判断方法
CN109470582B (zh) 露天矿山岩体结构面等效抗剪强度参数分级确定方法
CN109945822B (zh) 一种动态沉陷区房屋就地重建时间和区域测定方法
CN116821856A (zh) 一种基于岩碴比表面积的tbm破岩效率评价方法
CN105930660B (zh) 一种基于盲信息的堆积体滑坡安全评价方法
CN113536582B (zh) 一种获取溶洞顶板极限溶蚀耐久性的方法
CN104021280B (zh) 一种适用于尾粉砂发生管涌的临界水力梯度的计算方法
CN107060747B (zh) 钻井过程中钻遇裂缝发育带的预警方法及系统
CN113723026A (zh) 一种隧道掌子面前方煤层瓦斯涌出量预算方法

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