CN104239637A - 一种离散元爆堆形态模拟方法 - Google Patents
一种离散元爆堆形态模拟方法 Download PDFInfo
- Publication number
- CN104239637A CN104239637A CN201410472310.3A CN201410472310A CN104239637A CN 104239637 A CN104239637 A CN 104239637A CN 201410472310 A CN201410472310 A CN 201410472310A CN 104239637 A CN104239637 A CN 104239637A
- Authority
- CN
- China
- Prior art keywords
- discrete
- lumpiness
- rockmass
- model
- rock mass
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000004088 simulation Methods 0.000 claims abstract description 27
- 239000011435 rock Substances 0.000 claims abstract description 22
- 238000005422 blasting Methods 0.000 claims abstract description 19
- 238000011439 discrete element method Methods 0.000 claims abstract description 8
- 230000000877 morphologic effect Effects 0.000 claims description 18
- 238000013467 fragmentation Methods 0.000 claims description 15
- 238000006062 fragmentation reaction Methods 0.000 claims description 15
- 238000004880 explosion Methods 0.000 claims description 10
- 238000005192 partition Methods 0.000 claims description 9
- 239000002360 explosive Substances 0.000 claims description 4
- 238000010008 shearing Methods 0.000 claims description 3
- 238000013316 zoning Methods 0.000 abstract 1
- 238000006073 displacement reaction Methods 0.000 description 4
- 238000005474 detonation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F42—AMMUNITION; BLASTING
- F42D—BLASTING
- F42D3/00—Particular applications of blasting techniques
- F42D3/04—Particular applications of blasting techniques for rock blasting
-
- 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/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种离散元爆堆形态模拟方法,包括步骤:步骤1,将完整岩体的体积模量和剪切模量分别赋给人工节理的法向刚度和和切向刚度;步骤2,针对完整岩体模型,采用Harries爆破模型预测距炮孔中心不同距离处各点的预测爆破块度;步骤3,将完整岩体模型划分区域,以各区域中所有点的预测爆破块度的平均值为该区域离散块的预设块度;步骤4,根据各区域离散块的预设块度设置人工节理,并将人工节理设为倾斜正交,采用离散元法对完整岩体模型进行离散;步骤5,采用离散元法对离散后的完整岩体模型进行爆堆形态模拟。本发明方法操作简单,可显著提高爆堆形态模拟的效率和精准度,更符合工程实际,应用前景广阔。
Description
技术领域
本发明属于岩土工程中爆破数值模拟技术领域,尤其涉及一种离散元爆堆形态模拟方法。
背景技术
目前,工程爆破中爆堆形态模拟是爆破数值模拟中最重要的研究课题之一。随着计算机仿真技术发展,适用于离散元的数值模拟方法逐渐成为研究爆破的有力工具。离散元将不连续面看作离散块的边界条件,可以模拟大变形,允许离散块间的相对运动,并且不需要满足位移连续和变形协调条件,非常适用于爆堆形态模拟。
模拟爆堆形成过程时,需要将完整岩体预先离散为数目众多的离散块。所以,确定离散块的块度是爆堆形态模拟的基础。传统的离散元爆堆形态模拟方法中是将岩体模型离散为大小均匀的离散块,而离散块的预设块度的取值依靠经验或重复试验确定。因此,传统离散元爆堆形态模拟方法无法体现爆破本质,不仅工作量大,而且爆堆形态模拟的准确性也不高。
发明内容
针对现有技术存在的不足,本发明提供了一种更客观、更高效的离散元爆堆形态模拟方法。
为解决上述技术问题,本发明采用如下的技术方案:
一种离散元爆堆形态模拟方法,包括步骤:
步骤1,将完整岩体的体积模量和剪切模量分别赋给人工节理的法向刚度和和切向刚度;
步骤2,针对完整岩体模型,采用Harries爆破模型预测距炮孔中心不同距离处各点的预测爆破块度;
步骤3,将完整岩体模型划分区域,以各区域中所有点的预测爆破块度的平均值为该区域离散块的预设块度;
步骤4,根据各区域离散块的预设块度设置人工节理,并将人工节理设为倾斜正交,采用离散元法对完整岩体模型进行离散;
步骤5,采用离散元法对离散后的完整岩体模型进行爆堆形态模拟。
上述步骤2进一步包括子步骤:
2.1根据岩体属性和炸药属性,采用Harries爆破模型获得炮孔壁上应变K;
2.2根据炮孔壁上应变K,采用Harries爆破模型获得距炮孔中心不同距离处任意点的切应变ε;
2.3根据任意点切应变ε和岩体的动态极限抗拉应变T,采用Harries爆破模型获得爆破后距炮孔中心不同距离处的径向裂缝条数N;
2.4以距炮孔中心R处的相邻径向裂缝间距离为距炮孔中心距离R处各点的预测爆破块度L。
作为优选方案,将步骤4获得的模拟结果与实际获得的爆堆形状和爆破块度进行对比,并根据对比结果调整完整岩体模型划分的区域数量。
与现有技术相比,本发明具有如下特点和有益效果:
(1)基于岩石爆破的内在机理预设离散块的块度,相比传统方法更客观、更符合工程实际,能显著提高爆堆形态模拟的精准度。
(2)采用本发明方法获得的离散块数量更少,仅有传统方法获得的离散块数量的一半,因此,可大幅减少模拟的运算时间和资源占有,显著提高爆堆形态模拟效率,应用前景广阔。
(3)操作简单,采用Excel即可实现本发明,可有效减少工作量。
附图说明
图1为本发明方法的具体流程图;
图2为具体实施方式中完整岩体模型示意图;
图3为具体实施方式中完整岩体模型划分区域示意图;
图4为具体实施方式中完整岩体模型的离散结果示意图。
具体实施方式
本发明首先将完整岩体的力学参数值赋给人工节理的力学参数,保证人工离散化不影响岩体的力学性质;然后,根据哈里斯爆破理论计算岩体模型各处的预测爆破块度;依照数值计算需要将模型分区,以各区内的平均爆破块度作为离散块体的块度;最后,可根据爆堆模拟的结果来修正预设块度或模型分区。
本发明具体流程图见图1,下面将结合图1进一步说明本发明方法步骤。
步骤1,在离散元程序中为人工节理的力学参数赋值,将完整岩体的力学参数值赋给人工节理的力学参数,例如,将完整岩体的体积模量和剪切模量分别赋给人工节理的法向刚度和和切向刚度,以保证人工离散化不影响岩体的力学性质。完整岩体的力学参数值指爆破模拟参数,根据工程资料获得。
步骤2,基于Harries爆破模型(哈里斯爆破模型)预测离散块的块度。
本发明将预测爆破块度作为离散完整岩体模型的根据,针对完整岩体模型采用Harries爆破模型预测距离炮孔中心R处的爆破块度,具体如下:
(1)计算炮孔壁上应变K为:
式(1)中,μ是岩体泊松比,属于岩体力学参数,根据工程资料获得;P0表示炮孔壁上爆压,可通过计算获得,本具体实施中采用C-J爆轰理论计算炮孔壁上爆压P0;ρ为岩体密度;vp为岩体纵波声速,纵波声速为岩体自身属性,根据工程资料获得;γ为炸药的等熵指数,可根据炸药出厂参数计算获得。
(2)计算爆破时距炮孔中心R处任意点上的切应变ε。
爆破时,爆破应力波的传播使得与炮孔轴线垂直的平面内岩体质点产生径向位移,在半径为r的球形装药上,距炮孔中心R处任意点上的切应变ε为:
式(2)中,K为炮孔壁上应变;β为应变波吸收系数,根据实验资料或工程资料获得;r为炮孔装药半径。
(3)计算爆破后距炮孔中心R处的径向裂缝条数N。
爆破应力波作用下岩体作径向位移,当由径向位移衍生的切应变ε超过岩体的动态极限抗拉应变T时,岩石形成径向裂缝,动态极限抗拉应变T为岩体自身属性,根据工程资料获得。距炮孔中心R处的径向裂缝条数N为:
(4)以距炮孔中R处的相邻径向裂缝间距离为距炮孔中R处各点的预测爆破块度L:
步骤3,将完整岩体模型划分为n个区域,各区域的离散块预设块度采用预测爆破块度的平均值。
完整岩体模型区域的划分根据使用者的具体需求进行,一般划分为矩形区域,见图3。区域数量越多则爆破模拟准确度越高,但同时会增加运算时间及运算量。本具体实施中将完整岩体模型划分为6~10个区域。
对于区域i,其面积为Ai,则区域i内的预设离散块度:
式(5)中,∫LdAi表示区域i中各点的预测爆破块度之和。
步骤4,根据离散块的预设块度设置人工节理,并对完整岩体模型进行预离散,,预离散后的完整岩体模型见图4。
本发明中人工节理设为倾斜正交,即节理倾角设置为±45°。
步骤5,采用离散元法对预离散后的完整岩体模型进行爆堆形态模拟,将模拟结果与实际获得的爆堆形状和爆破块度进行对比,并根据对比结果调整完整岩体模型划分区域的数量。若模拟结果和实际情况相差较大,则增加完整岩体模型划分区域数量。
以上所述的具体实施例,对本发明的目的、技术方案和有益结果进行了进一步的详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之类,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种离散元爆堆形态模拟方法,其特征在于,包括步骤:
步骤1,将完整岩体的体积模量和剪切模量分别赋给人工节理的法向刚度和和切向刚度;
步骤2,针对完整岩体模型,采用Harries爆破模型预测距炮孔中心不同距离处各点的预测爆破块度;
步骤3,将完整岩体模型划分区域,以各区域中所有点的预测爆破块度的平均值为该区域离散块的预设块度;
步骤4,根据各区域离散块的预设块度设置人工节理,并将人工节理设为倾斜正交,采用离散元法对完整岩体模型进行离散;
步骤5,采用离散元法对离散后的完整岩体模型进行爆堆形态模拟。
2.如权利要求1所述的离散元爆堆形态模拟方法,其特征在于:
步骤2进一步包括子步骤:
2.1根据岩体属性和炸药属性,采用Harries爆破模型获得炮孔壁上应变K;
2.2根据炮孔壁上应变K,采用Harries爆破模型获得距炮孔中心不同距离处任意点的切应变ε;
2.3根据任意点切应变ε和岩体的动态极限抗拉应变T,采用Harries爆破模型获得爆破后距炮孔中心不同距离处的径向裂缝条数N;
2.4以距炮孔中心R处的相邻径向裂缝间距离为距炮孔中心距离R处各点的预测爆破块度L。
3.如权利要求1所述的离散元爆堆形态模拟方法,其特征在于,还包括:
将步骤4获得的模拟结果与实际获得的爆堆形状和爆破块度进行对比,并根据对比结果调整完整岩体模型划分的区域数量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410472310.3A CN104239637A (zh) | 2014-09-16 | 2014-09-16 | 一种离散元爆堆形态模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410472310.3A CN104239637A (zh) | 2014-09-16 | 2014-09-16 | 一种离散元爆堆形态模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104239637A true CN104239637A (zh) | 2014-12-24 |
Family
ID=52227691
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410472310.3A Pending CN104239637A (zh) | 2014-09-16 | 2014-09-16 | 一种离散元爆堆形态模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104239637A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339547A (zh) * | 2016-08-26 | 2017-01-18 | 中国水利水电科学研究院 | 一种爆破的数值模拟方法 |
CN109992841A (zh) * | 2019-03-11 | 2019-07-09 | 长江水利委员会长江科学院 | 一种爆破块度空间全方位分区耦合数值仿真方法 |
CN112991539A (zh) * | 2021-03-11 | 2021-06-18 | 广东锡源爆破科技股份有限公司 | 一种基于离散元的爆堆及块度分布的三维模拟方法 |
CN113284241A (zh) * | 2021-03-11 | 2021-08-20 | 广东锡源爆破科技股份有限公司 | 一种堆形为Weibull分布的爆堆表面岩块三维模拟方法 |
CN113343441A (zh) * | 2021-05-21 | 2021-09-03 | 武汉理工大学 | 基于筛分法仿真的岩石爆破块度测量方法 |
CN114218831A (zh) * | 2021-12-15 | 2022-03-22 | 武汉市市政工程机械化施工有限公司 | 一种通用爆破数值模拟方法 |
EP4153937A4 (en) * | 2020-05-21 | 2024-07-10 | Dyno Nobel Inc | METHODS, DEVICES AND SYSTEMS FOR MOVING BLASTS WITH DIFFERENT ELEMENTS |
-
2014
- 2014-09-16 CN CN201410472310.3A patent/CN104239637A/zh active Pending
Non-Patent Citations (1)
Title |
---|
周旺潇等: "爆破漏斗形成过程数值模拟的几个关键问题", 《爆破》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339547A (zh) * | 2016-08-26 | 2017-01-18 | 中国水利水电科学研究院 | 一种爆破的数值模拟方法 |
CN106339547B (zh) * | 2016-08-26 | 2019-04-23 | 中国水利水电科学研究院 | 一种爆破的数值模拟方法 |
CN109992841A (zh) * | 2019-03-11 | 2019-07-09 | 长江水利委员会长江科学院 | 一种爆破块度空间全方位分区耦合数值仿真方法 |
CN109992841B (zh) * | 2019-03-11 | 2022-12-06 | 长江水利委员会长江科学院 | 一种爆破块度空间全方位分区耦合数值仿真方法 |
EP4153937A4 (en) * | 2020-05-21 | 2024-07-10 | Dyno Nobel Inc | METHODS, DEVICES AND SYSTEMS FOR MOVING BLASTS WITH DIFFERENT ELEMENTS |
CN112991539A (zh) * | 2021-03-11 | 2021-06-18 | 广东锡源爆破科技股份有限公司 | 一种基于离散元的爆堆及块度分布的三维模拟方法 |
CN113284241A (zh) * | 2021-03-11 | 2021-08-20 | 广东锡源爆破科技股份有限公司 | 一种堆形为Weibull分布的爆堆表面岩块三维模拟方法 |
CN113284241B (zh) * | 2021-03-11 | 2022-03-15 | 广东锡源爆破科技股份有限公司 | 一种堆形为Weibull分布的爆堆表面岩块三维模拟方法 |
CN113343441A (zh) * | 2021-05-21 | 2021-09-03 | 武汉理工大学 | 基于筛分法仿真的岩石爆破块度测量方法 |
CN114218831A (zh) * | 2021-12-15 | 2022-03-22 | 武汉市市政工程机械化施工有限公司 | 一种通用爆破数值模拟方法 |
CN114218831B (zh) * | 2021-12-15 | 2024-04-09 | 武汉市市政工程机械化施工有限公司 | 一种通用爆破数值模拟方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104239637A (zh) | 一种离散元爆堆形态模拟方法 | |
He et al. | An innovative approach for gob-side entry retaining in thick coal seam longwall mining | |
CN104361211B (zh) | 岩石统计损伤本构模型的构建和应用方法 | |
Zhu et al. | Numerical investigation of blasting-induced crack initiation and propagation in rocks | |
CN105787206A (zh) | 一种公路隧道锚杆支护设计方法 | |
CN114722614B (zh) | 一种基于二维-三维协同的爆破动力响应模拟方法 | |
CN113255175A (zh) | 一种层状软岩爆破数值模拟方法 | |
CN101787898A (zh) | 一种煤与瓦斯突出的预报方法 | |
CN114297864B (zh) | 一种受陡缓倾角控制的碎裂松动岩体边坡稳定性分析方法 | |
Mahabadi et al. | Development of a new fully-parallel finite-discrete element code: Irazu | |
CN115859714A (zh) | 一种基于fem-dem联合仿真的岩石爆破全过程模拟方法 | |
CN105443102A (zh) | 一种水平井层内爆炸裂缝体建模方法 | |
CN103745040A (zh) | 一种露天矿边坡爆破过程稳定性的分析方法 | |
Li et al. | Numerical Simulation of Blast Vibration and Crack Forming Effect of Rock‐Anchored Beam Excavation in Deep Underground Caverns | |
Gharehdash et al. | Numerical investigation on fracturing of rock under blast using coupled finite element method and smoothed particle hydrodynamics | |
CN115711091A (zh) | 一种煤矿井下水力造穴工艺模拟优化方法与系统 | |
CN115510778A (zh) | 一种陆相页岩储层无限级压裂工艺优化方法及系统 | |
Hebert et al. | Three-dimensional simulation of cave initiation, propagation and surface subsidence using a coupled finite difference–cellular automata solution | |
Yu et al. | Numerical investigation of blast-induced rock movement characteristics in open-pit bench blasting using bonded-particle method | |
CN107545078A (zh) | 流固耦合模型的构建方法及储藏缝洞变形的分析方法 | |
WO2023124938A1 (zh) | 测试完井管柱的射孔爆轰分析方法及装置 | |
Cao et al. | Study on the guiding mechanism of multiple and empty holes under explosion load | |
Yu et al. | Numerical experimental study on ore dilution in sublevel caving mining | |
CN115705454A (zh) | 一种基于相场法的裂缝扩展模拟压裂设计优化方法 | |
CN103077556B (zh) | 油井出砂的三维数值模型设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20141224 |