CN113435087B - 一种分析洞穴围岩局部稳定性的方法 - Google Patents
一种分析洞穴围岩局部稳定性的方法 Download PDFInfo
- Publication number
- CN113435087B CN113435087B CN202110715849.7A CN202110715849A CN113435087B CN 113435087 B CN113435087 B CN 113435087B CN 202110715849 A CN202110715849 A CN 202110715849A CN 113435087 B CN113435087 B CN 113435087B
- Authority
- CN
- China
- Prior art keywords
- rock
- cave
- dimensional numerical
- numerical model
- model
- 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
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/04—Architectural design, interior design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及采矿技术领域,特别是分析洞穴围岩局部稳定性的方法,利用有限元模拟软件建立岩体的三维数值模型时,在三维数值模型中设置节理面,对三维数值模型进行模拟洞穴开挖后,计算洞穴围岩处节理面的抗剪安全系数,将抗剪安全系数与预设值进行比较,从而得出洞穴围岩的局部稳定性。本发明以对洞穴围岩的局部结构——节理面分析为切入点,建立包含节理面的岩体三维数值模型,计算出洞穴围岩各个节理面的抗剪安全系数,通过抗剪安全系数的分析,直观了解到洞穴围岩局部部位的稳定性情况,并可对局部失稳处采取相应的加固措施,避免对现场施工人员造成安全隐患,以及防止局部失稳状态的进一步扩展恶化。
Description
技术领域
本发明涉及采矿技术领域,特别是一种分析洞穴围岩局部稳定性的方法。
背景技术
现有分析洞穴围岩稳定性的方法包括如下步骤:利用有限元模拟软件建立岩体的三维数值模型;对三维数值模型赋予岩体本构模型,用以模拟岩体内的应力应变情况;对三维数值模型进行模拟洞穴开挖;将洞穴围岩作为整体进行受力分析。采用上述分析方法只能评估洞穴围岩是否会发生整体大面积塌陷、失稳,即只是对洞穴围岩的整体稳定性进行判断,即使整体稳定性满足要求,洞穴围岩局部处也可能出现失稳情况,如洞穴围岩局部发生岩土块脱落,如果不能对洞穴围岩局部稳定性进行有效分析并采取相应加固措施,局部的失稳情况会对现场施工人员造成安全隐患,甚至逐渐扩展成洞穴围岩的整体大面积塌陷、失稳,故亟待找到一种分析洞穴围岩局部稳定性的方法。
发明内容
本发明目的在于针对现有技术在分析洞穴围岩稳定性时只能评估洞穴围岩是否会发生整体大面积塌陷、失稳,却分析不了洞穴围岩局部岩体的失稳情况的问题,提供一种分析洞穴围岩局部稳定性的方法。
为了实现上述目的,本发明采用的技术方案为:
一种分析洞穴围岩局部稳定性的方法,包括如下步骤:
利用有限元模拟软件建立岩体三维数值模型,对建立完成的岩体三维数值模型进行网格划分;
对三维数值模型赋予岩体本构模型,用以模拟岩体内的应力应变情况;
对赋予岩体本构模型的岩体三维数值模型进行洞穴模拟开挖,利用有限元模拟软件建立岩体三维数值模型时,在岩体三维数值模型中设置节理面;
对岩体三维数值模型进行洞穴模拟开挖后,计算洞穴围岩处各个节理面的抗剪安全系数,将各个节理面抗剪安全系数与预设值进行比较,若某一节理面的抗剪安全系数大于或等于预设值,则该节理面处的洞穴围岩不发生局部失稳;若某一节理面的抗剪安全系数小于预设值,则该节理面的洞穴围岩发生局部失稳。
其中,选择节理面的参数如倾角和倾向对本领域技术人员而言是公知常识,根据现场实际情况勘探出节理面的参数这对本领域技术人员而言是常规技术手段。预设值可以根据实际情况中岩体所处的位置和工况选取,对本领域技术人员而言为公知常识。
本发明在分析洞穴围岩局部稳定性时引入岩体的局部结构——节理面,利用有限元模拟软件建立岩体三维数值模型时,在岩体三维数值模型中设置节理面,并计算出洞穴围岩各个节理面的抗剪安全系数,即可得出洞穴围岩局部处的安全系数。通过抗剪安全系数与预设值的分析,能够明确洞穴围岩不同部位的稳定性情况,并可对局部失稳处采取相应的加固措施,避免对现场施工人员造成安全隐患,以及防止局部失稳状态的进一步扩展恶化。
优选的,对赋予了岩体本构模型的三维数值模型进行洞穴模拟开挖包括步骤:
M1、对岩体三维数值模型设置边界条件,用以模拟岩体周围边界的限制情况;
M2、对岩体三维数值模型设置重力加速度,用以模拟岩体在重力作用下的自然沉降,以得到岩体三维数值模型的待开挖状态;
M3、在岩体三维数值模型的待开挖状态下进行洞穴模拟开挖。
优选的,在步骤M3后清除岩体三维数值模型所产生的位移和时间,以得到岩体三维数值模型的初始计算状态,在该初始计算状态下计算洞穴围岩处各个节理面的抗剪安全系数。
优选的,对岩体三维数值模型赋予的岩体本构模型为摩尔库伦本构模型。
优选的,划分岩体三维数值模型的网格长度为2-4m。
优选的,预设值为1.30。
综上,由于采用了上述技术方案,本发明的有益效果至少是:
1.本发明在分析洞穴围岩局部稳定性时引入岩体的局部结构特征——节理面,通过建立包含节理面的三维数值模型,计算洞穴围岩各个节理面的抗剪安全系数,通过抗剪安全系数与预设值的分析,能够明确洞穴围岩局部部位的稳定性情况,对局部失稳处采取相应的加固措施,避免对现场施工人员造成安全隐患。
附图说明
图1为块体模型示意图;
图2为节理面与块体模型间倾角示意图;
图3为节理面与块体模型间倾向示意图;
图4为实施例1中三组节理面位置关系图;
图5为实施例1中一种分析洞穴围岩局部稳定性的方法流程图;
图中标记:1-块体模型,2-洞穴模型,3-节理面,31-第一节理面,32-第二节理面,33-第三节理面。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
实施例1
如图1-5所示,为本实施例一种分析洞穴围岩局部稳定性的方法,包括如下步骤:
利用有限元模拟软件建立岩体三维数值模型,对建立完成的岩体三维数值模型进行网格划分;
对三维数值模型赋予岩体本构模型,用以模拟岩体内的应力应变情况;
对赋予岩体本构模型的岩体三维数值模型进行洞穴模拟开挖;
利用有限元模拟软件建立岩体三维数值模型时,在岩体三维数值模型中设置节理面3;
对岩体三维数值模型进行洞穴模拟开挖后,计算洞穴围岩处各个节理面3的抗剪安全系数,将各个节理面3抗剪安全系数与预设值进行比较,若某一节理面3的抗剪安全系数大于或等于预设值,则该节理面3处的洞穴围岩不发生局部失稳;若某一节理面3的抗剪安全系数小于预设值,则该节理面3的洞穴围岩发生局部失稳。
在本实施例中,建立包含节理面3的岩体三维数值模型包括步骤:
S1、可采用离散元数值模拟软件建立块体模型1,在其它可行的实施例中,也可用岩土工程中广泛采用的其它数值模拟方法如有限元法、有限差分法;
S2、如图1所示,可对块体模型1分为两组,一组为洞穴模型2,另一组为洞穴周围岩体模型,块体模型1和洞穴模型2的尺寸由现场勘探得出,在本实施例中,块体模型1长100m,宽70m,高100m。洞穴模型2的截面形状以及在块体模型1中的位置同样根据现场勘探结果设置,在本实施例中,洞穴模型2截面形状为二分之一椭圆,洞穴模型2位于块体模型1中部位置;
S3、对块体模型1设置节理面3,节理面3的参数包括倾角、倾向,以及节理面3之间的间距,节理面3设置的组数和上述参数由现场勘探结果设置,在本实施例中,设置3组节理面3,第一组倾角12度,倾向65度,间距10m,如图4中所示的第一节理面31;第二组倾角60度,倾向0度,间距10m,如图4中所示的第二节理面32;第三组倾角80度,倾向55度,间距10m,如图4中所示的第三节理面33。如图2所示,倾角α表示的是节理面3与xoy平面的夹角与x轴负方向所交的角,如图3所示,倾向β表示节理面3与xoz平面的夹角与x轴负方向所交的角,设置好的三组节理面3位置关系如图4所示。
S4、对块体模型1划分网格,以使模型能够进行受力分析,网格长度可设置为2-4m,如可选3m。
为了模拟土体真实的应力应变关系,需对建立完成的三维数值模型赋予岩体本构模型,如剑桥模型、修正剑桥模型或摩尔库伦本构模型。在本实施例中,对三维数值模型赋予摩尔库伦本构模型,并对岩体和节理面3赋予相应的强度参数,强度参数包括岩体的密度、体积模量、剪切模量、内摩擦角、粘聚力,节理面3的法向刚度、剪切刚度、内摩擦角、粘聚力。岩体和节理面3相应的强度参数由现场勘探得出,在本实施例中,模型岩体密度为2480kg/m3,体积模量为27.8Gpa,剪切模量为20.8Gpa,内摩擦角为35度、粘聚力为5.8Mpa、抗拉强度为4.6Mpa,节理面3的法向刚度为2.11Gpa,剪切刚度为0.68Gpa,内摩擦角为19.5度,粘聚力为3.7Mpa,抗拉强度为0.163Mpa。
在本实施例中,对赋予了岩体本构模型的三维数值模型进行模拟洞穴开挖包括步骤:
M1、对三维数值模型设置边界条件,以模拟岩体周围边界的限制情况。在本实施例中,三维数值模型边界条件为:模型底部固定,模型四周固定,模型顶部为自由边界,模拟地面状态;
M2、对三维数值模型设置重力加速度,用以模拟岩体内部的自然沉降,以得到岩体三维数值模型的待开挖状态;
M3、在岩体三维数值模型的待开挖状态下进行洞穴模拟开挖。进一步的,在步骤M3后清除岩体三维数值模型所产生的位移和产生所述位移经过的时间,以得到岩体三维数值模型的初始计算状态,在该初始计算状态下计算洞穴围岩处各个节理面3的抗剪安全系数,计算洞穴围岩处各个节理面3的抗剪安全系数的步骤如下:
利用软件遍历洞穴围岩处各个节理面3的面积和对应的切向力,将切向力除以对应的节理面3面积得到剪应力。提取各个节理面3相应的粘聚力和内摩擦角,并根据摩尔库伦本构模型计算出节理面3各个部分的抗剪强度,将抗剪强度除以对应节理面3所承受的剪应力,便得到各个节理面3相对应的抗剪安全系数。理论上抗剪安全系数等于1为极限稳定状态,小于1则会发生破坏,但是实际情况下为了保证安全,抗剪安全系数需大于1,可以根据《建筑边坡工程技术规范》GB50330-2013中边坡稳定性安全系数的规范标准作为参考,亦可根据设计人员根据实际情况自行设定。
表1:边坡稳定性安全系数规范标准
安全等级 | 一级边坡 | 二级边坡 | 三级边坡 |
一般工况安全系数 | 1.35 | 1.30 | 1.25 |
地震工况安全系数 | 1.15 | 1.10 | 1.05 |
临时工况安全系数 | 1.25 | 1.20 | 1.15 |
在本实施例中,在抗剪安全系数低于1.30的区域采取注浆加固和增加锚杆预紧力的支护措施。
实施例2
实施例1洞穴的形成过程是通过开挖形成的,实际中洞穴也可如溶洞的形成过程。本实施例与实施例1不同之处在于,洞穴形成的过程不是通过模拟开挖而成,而是通过模拟自然状态下溶洞的形成机制而形成的,其余分析洞穴围岩局部稳定性的方法与实施例1相同。
以上仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种分析洞穴围岩局部稳定性的方法,包括如下步骤:利用有限元模拟软件建立岩体三维数值模型,对建立完成的岩体三维数值模型进行网格划分;对三维数值模型赋予岩体本构模型,用以模拟岩体内的应力应变情况;对赋予岩体本构模型的岩体三维数值模型进行洞穴模拟开挖,其特征在于,利用有限元模拟软件建立岩体三维数值模型时,在岩体三维数值模型中设置节理面(3);对岩体三维数值模型进行洞穴模拟开挖后,计算洞穴围岩处各个节理面(3)的抗剪安全系数,将各个节理面(3)抗剪安全系数与预设值进行比较,若某一节理面(3)的抗剪安全系数大于或等于预设值,则该节理面(3)处的洞穴围岩不发生局部失稳;若某一节理面(3)的抗剪安全系数小于预设值,则该节理面(3)的洞穴围岩发生局部失稳;
对赋予了岩体本构模型的三维数值模型进行洞穴模拟开挖包括步骤:M1、对岩体三维数值模型设置边界条件,用以模拟岩体周围边界的限制情况;M2、对岩体三维数值模型设置重力加速度,用以模拟岩体在重力作用下的自然沉降,以得到岩体三维数值模型的待开挖状态;M3、在岩体三维数值模型的待开挖状态下进行洞穴模拟开挖;
在步骤M3后清除岩体三维数值模型所产生的位移和时间,以得到岩体三维数值模型的初始计算状态,在该初始计算状态下计算洞穴围岩处各个节理面(3)的抗剪安全系数;计算洞穴围岩处各个节理面(3)的抗剪安全系数的步骤如下:利用软件遍历洞穴围岩处各个节理面(3)的面积和对应的切向力,将切向力除以对应的节理面(3)面积得到剪应力,提取各个节理面(3)相应的粘聚力和内摩擦角,并根据摩尔库伦本构模型计算出节理面(3)各个部分的抗剪强度,将抗剪强度除以对应节理面(3)所承受的剪应力,便得到各个节理面(3)相对应的抗剪安全系数。
2.根据权利要求1所述的一种分析洞穴围岩局部稳定性的方法,其特征在于,对岩体三维数值模型赋予的岩体本构模型为摩尔库伦本构模型。
3.根据权利要求1所述的一种分析洞穴围岩局部稳定性的方法,其特征在于,划分岩体三维数值模型的网格长度为2-4m。
4.根据权利要求1所述的一种分析洞穴围岩局部稳定性的方法,其特征在于,预设值为1.30。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110715849.7A CN113435087B (zh) | 2021-06-24 | 2021-06-24 | 一种分析洞穴围岩局部稳定性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110715849.7A CN113435087B (zh) | 2021-06-24 | 2021-06-24 | 一种分析洞穴围岩局部稳定性的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113435087A CN113435087A (zh) | 2021-09-24 |
CN113435087B true CN113435087B (zh) | 2022-04-22 |
Family
ID=77754852
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110715849.7A Active CN113435087B (zh) | 2021-06-24 | 2021-06-24 | 一种分析洞穴围岩局部稳定性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113435087B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114428990B (zh) * | 2022-01-27 | 2023-04-28 | 中铁二院工程集团有限责任公司 | 一种基于AutoCAD自适应曲线趋势的自动填图方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105550441A (zh) * | 2015-12-12 | 2016-05-04 | 山东科技大学 | 一种基于连续介质的工程岩体破裂劣化数值模拟方法 |
WO2018149526A1 (en) * | 2017-02-17 | 2018-08-23 | Siemens Wind Power A/S | A wind turbine blade, a wind turbine and a method for producing a wind turbine blade |
CN109492262A (zh) * | 2018-10-16 | 2019-03-19 | 山东科技大学 | 一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 |
CN109493380A (zh) * | 2017-09-11 | 2019-03-19 | 重庆大学 | 一种岩石节理面抗剪试验中不规则剪切面面积的计算方法 |
CN111369156A (zh) * | 2020-03-11 | 2020-07-03 | 中铁(天津)隧道工程勘察设计有限公司 | 一种超大跨洞室岩石隧道开挖围岩局部稳定性评价方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107067333B (zh) * | 2017-01-16 | 2022-12-20 | 长沙矿山研究院有限责任公司 | 一种高寒高海拔高陡边坡稳定性监控方法 |
CN108959759A (zh) * | 2018-06-28 | 2018-12-07 | 河海大学 | 基于数值模拟确定柱状节理岩体各向异性比的方法 |
-
2021
- 2021-06-24 CN CN202110715849.7A patent/CN113435087B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105550441A (zh) * | 2015-12-12 | 2016-05-04 | 山东科技大学 | 一种基于连续介质的工程岩体破裂劣化数值模拟方法 |
WO2018149526A1 (en) * | 2017-02-17 | 2018-08-23 | Siemens Wind Power A/S | A wind turbine blade, a wind turbine and a method for producing a wind turbine blade |
CN109493380A (zh) * | 2017-09-11 | 2019-03-19 | 重庆大学 | 一种岩石节理面抗剪试验中不规则剪切面面积的计算方法 |
CN109492262A (zh) * | 2018-10-16 | 2019-03-19 | 山东科技大学 | 一种利用数值模拟分析非均匀分布裂隙巷道稳定性的方法 |
CN111369156A (zh) * | 2020-03-11 | 2020-07-03 | 中铁(天津)隧道工程勘察设计有限公司 | 一种超大跨洞室岩石隧道开挖围岩局部稳定性评价方法 |
Non-Patent Citations (4)
Title |
---|
3-Dimensional Modelling for the Rock Slope Excavation and Reinforcement by the Bolts in Civil Engineering;Wen-xiang Peng等;《2010 Second International Conference on Computer Engineering and Applications》;20100412;134-137 * |
甘肃红层工程地质特性与边坡稳定性研究;王骑虎;《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑》;20190215;C034-26 * |
节理特性对隧道围岩稳定性影响的研究;王贵君等;《河北工业大学学报》;20170215;103-107 * |
节理裂隙岩体隧道围岩变形及稳定性研究——以烧锅隧道为例;任杨茹;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20190215;C034-356 * |
Also Published As
Publication number | Publication date |
---|---|
CN113435087A (zh) | 2021-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Merrien-Soukatchoff et al. | Use of an integrated discrete fracture network code for stochastic stability analyses of fractured rock masses | |
CN105787206A (zh) | 一种公路隧道锚杆支护设计方法 | |
Jiang et al. | A new method for analyzing the stability of rock wedges | |
CN105550441A (zh) | 一种基于连续介质的工程岩体破裂劣化数值模拟方法 | |
Li et al. | Effects of in-situ stress on the stability of a roadway excavated through a coal seam | |
CN112488455A (zh) | 边坡稳定性分析与评价方法、系统、介质、设备、终端 | |
CN114297864B (zh) | 一种受陡缓倾角控制的碎裂松动岩体边坡稳定性分析方法 | |
CN111553101A (zh) | 一种隧道开挖上覆岩层开裂预报方法及围岩支护方法 | |
Ji et al. | Numerical simulation of hydraulic fracturing in earth and rockfill dam using extended finite element method | |
CN113435087B (zh) | 一种分析洞穴围岩局部稳定性的方法 | |
CN113536414B (zh) | 基于三维建模的岩质边坡稳定性分析方法、系统及介质 | |
Xia et al. | Study on instability mechanism of extraction structure under undercut space based on thin plate theory in block caving method | |
CN106372295B (zh) | 砂土岩溶地层中盾构与溶洞安全水平距离的确定方法 | |
Tan et al. | A new mining scheme for hanging-wall ore-body during the transition from open pit to underground mining: a numerical study | |
CN109543338A (zh) | 一种地震波作用下基坑稳定性安全系数的计算方法 | |
Yang et al. | Numerical analysis of the load transfer and deformation in a soil nailed slope | |
Zheng et al. | Smart and fast reinforcement design for anti-dip bedding rock slopes | |
Wei et al. | Numerical study of damage to rock surrounding an underground coal roadway excavation | |
CN115659741A (zh) | 一种sv波斜入射下山岭隧道洞口段结构地震反应分析方法 | |
Yan et al. | Structural Safety Evaluation of Tunnel Based on the Dynamic Monitoring Data during Construction | |
Wang et al. | Flowchart of DEM modeling stability analysis of large underground powerhouse caverns | |
Mahboubi et al. | Conventional and advanced numerical methods of rock slope stability analysis, a comparison study, gotvand dam right abutment (Iran) case study.” | |
Lemos | Computational topics in discrete element modeling | |
Qiu et al. | Stability Evaluation and Structural Parameters Optimization of Stope Based on Area Bearing Theory. Minerals 2022, 12, 808 | |
CN116502386B (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 |