爆破

粒径/mm	0~20	20~40	40~60	60~80	80~100	100~200	200~400	400~600	600以上
占比/%	7.51	7.22	6.25	6.17	10.53	18.01	11.77	14.87	17.67

粒径/mm	0~20	20~40	40~60	60~80	80~100	100~200	200~400	400~600	600以上
占比/%	7.51	7.22	6.25	6.17	10.53	18.01	11.77	14.87	17.67

裂隙岩体精细化数值模型构建与爆破模拟

PDF下载

叶海旺 ¹ , 韦文蓬 ¹ , 周汉红 ² , 余梦豪 ¹ , 李兴旺 ¹ , 雷涛 ¹ , 温颖 ¹ , 王其洲 ¹ , 石斌宏 ³ , 于燕 ³ , 张生 ³ , 钟传山 ³ , 高玉文 ⁴ , Abdou Mohamed Abdelkader HASSAN ¹

爆破 | 理论与技术探索 2023,40(4): 44-51

收起

爆破 | 理论与技术探索 2023, 40(4): 44-51

裂隙岩体精细化数值模型构建与爆破模拟

全屏

叶海旺¹, 韦文蓬¹, 周汉红², 余梦豪¹, 李兴旺¹, 雷涛¹, 温颖¹, 王其洲¹, 石斌宏³, 于燕³, 张生³, 钟传山³, 高玉文⁴, Abdou Mohamed Abdelkader HASSAN¹

作者信息

^1.武汉理工大学　资源与环境工程学院，武汉　430070

^2.武汉市公安局　治安管理局，武汉　430077

^3.华新水泥股份有限公司，武汉　430200

^4.五矿勘查开发有限公司，北京　100044

叶海旺（1971-），男，教授、博士、博导，主要从事采矿、爆破、安全等方面的教学和科研工作，（E-mail）yehaiwang369@hotmail.com。

YE Hai-wang (1971-), male, professor, Ph. D, doctoral supervisor, mainly engaged in teaching and scientific research in the fields of mining, blasting and safety, (E-mail) yehaiwang369@hotmail.com.

通讯作者:

雷涛（1983-），男，博士、讲师，主要从事采矿、安全、数字化矿山等方面的教学和科研工作，（E-mail）leitao539@163.com。

Refined Numerical Model Construction and Blasting Simulation of Fractured Rock Mass

Hai-wang YE¹, Wen-peng WEI¹, Han-hong ZHOU², Meng-hao YU¹, Xing-wang LI¹, Tao LEI¹, Ying WEN¹, Qi-zhou WANG¹, Bin-hong SHI³, Yan YU³, Sheng ZHANG³, Chuan-shan ZHONG³, Yu-wen GAO⁴, Abdou Mohamed Abdelkader HASSAN¹

Affiliations

^1.School of Resources and Environment Engineering, Wuhan University of Technology, Wuhan 430070, China

^2.Security Administration of Wuhan Public Security Bureau, Wuhan 430077, China

^3.Huaxin Cement Co., Ltd., Wuhan 430200, China

^4.Minmetals Exploration and Development Co., Ltd., Beijing 100044, China

出版时间: 2023-12-01 doi: 10.3963/j.issn.1001-487X.2023.04.006

文章导航

摘要

收起

裂隙岩体爆破数值模拟过程中，裂隙分布模型的真实性，是影响数值模拟结果的关键因素之一，现有节理裂隙构建方法难以实现三维分布上的复杂表征或是实现操作过于繁杂。为探索在LS-DYNA软件中构建三维分布上真实、复杂表征的节理裂隙模型简单可行的操作方法，通过MATLAB软件对爆破数值模型K文件进行解析重组，按照节理裂隙三维分布规律进行勾勒并赋予节理裂隙材料本构，构建出符合节理裂隙分布规律的节理岩体三维精细化数值模型。以某露天石灰石矿山为实例，开展节理裂隙三维分布规律统计分析并在台阶爆破数值模型中进行还原，基于裂隙岩体精细化数值模型，开展露天台阶爆破数值模拟试验和现场爆破试验对比研究。结果表明：采用解析重组方法构建的数值节理模型与实际节理裂隙的误差小于13%；节理裂隙面使岩体内的损伤分布规律发生了变化，与完整岩体相比，损伤范围增加了12.04%，粒径处于0~100 mm的占比降低了8.11%；爆破模拟得到的损伤结果与现场破碎效果接近，统计粒径在0~100 mm之间的碎块的占比，误差为4.16%。解析重组法构建三维节理裂隙，可行、操作简便，数值模拟结果接近实际。

关键词

裂隙岩体 / 裂隙分布规律 / 地质模型 / LS-DYNA / 爆破数值模拟 / 爆破试验

Abstract

收起

The authenticity of fracture distribution model is one of the key factors during numerical simulation of blasting in jointed rock mass, which would obviously affect the numerical simulation results. It is hard to represent the complex three-dimensional joint distribution in the existing joint construction method. To explore a simple and feasible operation method for constructing the complex 3D joint model in LS-DYNA software, a K file of blasting numerical model was analyzed and reorganized by MATLAB software. Furthermore, a 3D refined numerical model for jointed rock mass was constructed by the 3D joint distribution law and the constitutive joint model parameters. Finally, a statistical analysis of the three-dimensional joint distribution law was carried out in an open-pit limestone mine, and the joints were reconstructed in the numerical model of a bench blasting. Consequently, a comparative study of the numerical simulation and the field blasting test for open-pit bench blasting was carried out. The results show that the error between the joints built in the numerical model and the actual joints is less than 13%. The joint surface changes the damage distribution of the rock mass. Compared with the intact rock mass, the damage rock mass range increases by 12.04%, and the proportion of fragments with the size of 0~100 mm decreases by 8.11%. The damage results obtained by blasting simulation are close to the field rock breaking effect, and the percentage error of fragments with the size of 0~100 mm is 4.16%. The analytical reconstruction method is feasible and easy to represent the complex three-dimensional joint distribution, and the numerical results are close to experimental results.

Key words

jointed rock mass / joint distribution law / geological model / LS-DYNA / blasting simulation / blasting test

引用本文

叶海旺, 韦文蓬, 周汉红, 余梦豪, 李兴旺, 雷涛, 温颖, 王其洲, 石斌宏, 于燕, 张生, 钟传山, 高玉文, Abdou Mohamed Abdelkader HASSAN. 裂隙岩体精细化数值模型构建与爆破模拟. 爆破, 2023 , 40 (4) : 44 -51 . DOI: 10.3963/j.issn.1001-487X.2023.04.006

Hai-wang YE, Wen-peng WEI, Han-hong ZHOU, Meng-hao YU, Xing-wang LI, Tao LEI, Ying WEN, Qi-zhou WANG, Bin-hong SHI, Yan YU, Sheng ZHANG, Chuan-shan ZHONG, Yu-wen GAO, Abdou Mohamed Abdelkader HASSAN. Refined Numerical Model Construction and Blasting Simulation of Fractured Rock Mass[J]. Blasting, 2023 , 40 (4) : 44 -51 . DOI: 10.3963/j.issn.1001-487X.2023.04.006

正文

收起

裂隙岩体爆破时，节理裂隙分布特性是影响岩体爆破破碎效果的关键因素。在爆炸应力作用下，节理裂隙面易产生应力集中、破坏放大及引导破坏等效应，对岩体动力学响应特性及破坏规律具有重要影响^[1-3]。节理裂隙的非均质性、不连续性、各向异性和不可全窥性^[4]，使得节理裂隙岩体的爆破动力响应特性研究工作难度增大。目前，数值模拟是研究裂隙岩体爆破响应特性的主流手段之一，而不规则节理裂隙构建的准确性和合理性将直接影响到数值模拟结果的精度。

在ANSYS中直接建立节理裂隙是开展节理岩体爆破数值模拟最简便的方法，金李通过该种方式构建了长宽高为10 m×5 m×10 m的节理岩体^[5]，马小满^[6]、周文海等则构建规则^[7]、条带状的节理裂隙，规则分布的节理对应力传播、损伤破坏规律的影响较为清晰、直观。叶海旺^[8]、吴立辉等通过在ANSYS中建立模型之后利用LS/PREPOST对模型进行修改的方式^[9]，构建出含X型节理裂隙的爆破数值模型，完整表述出了节理裂隙倾角、间距等二维平面内的分布特征。基于3DEC、UDEC、PFC等软件构建节理裂隙地质模型基础上，ZL Wang^[10]、唐海^[11]、汤云坤提出将爆破应力曲线加载到上述软件中开展数值模拟的方法^[12]，实现了节理裂隙复杂、精细化表征下的爆破数值模拟。直接法和修改法可以实现二维节理地快速构建，节理岩体爆破数值试验中简便可行，但在倾向、倾角、间距和延伸长度等节理裂隙信息的复杂表征上较为困难。且岩体中因炸药爆轰而产生的应力随着传爆方向、传爆距离的变化而变化^[13]，需加载大量爆破荷载曲线才能实现节理裂隙岩体的爆破数值模拟，其操作过程繁琐。现有裂隙岩体爆破数值模拟，在节理裂隙的构建上，难以实现其在三维分布上的复杂表征或是实现操作过于繁杂。

基于已有研究成果，提出爆破数值模型解析重组建模方法，构建符合实际节理裂隙分布规律的岩体爆破数值模型，开展节理裂隙岩体的露天台阶爆破数值模拟研究，将模拟结果与现场试验结果进行对比分析。

1　精细化裂隙岩体数值模型构建方法

收起

1.1　节理裂隙解析重组建模方法简述

采用MATLAB软件进行编程，实现爆破数值模型K文件的解析，得到数据化的爆破数值模型库，按照节理裂隙产状（倾角、倾向、间距、延伸长度）分布函数对数据库进行筛分、重组，修改其材料特征，将修改后的数据库编译成新的爆破数值模型K文件。详细流程如图1所示。

1.2　节理裂隙参数定义

节理裂隙的三维空间分布特性参量倾角α、倾向β、间距γ、节理裂隙迹长l通过采样、统计分析可得，节理裂隙延伸长度w通过节理裂隙迹长l计算得到。

以绝对地理坐标角度34°为X轴正方向，以垂直坡面向外（304°）为Y轴正方向，即绝对地理坐标轴顺时针旋转34°得到ANSYS数值模拟软件中建模时的相对坐标轴，r（x₀，y₀，z₀）为节理裂隙面中心点，节理裂隙面与单位块体的交点分别为A、B、C、D，EF为水平面上的走向线。倾向向量

位于水平面上，垂直于走向线EF，

与x轴负方向的顺时针夹角为β，节理裂隙面与水平面的夹角为α。数值模型中节理裂隙倾向、倾角定义如图2所示。节理裂隙间距γ为岩层出露断面上相邻节理的距离，节理裂隙迹长l为岩层出露断面上节理裂隙的长度，延伸长度w为节理裂隙三维空间内的延展长度。

1.3　解析数据重组算法

在数据库所表征的三维模型内虚构出测线截面及测线，根据边坡测线法得到节理裂隙间距分布函数，在测线上筛选节理裂隙中心点，以节理裂隙中心点为起点绘制出倾角单位向量

，走向单位向量

（如图2所示），在裂隙中心点、

、

及延伸长度w等控制参数下，计算出节理裂隙面A、B、C、D四个点的坐标，确定节理裂隙平面。

(1)

(2)

走向单位向量

水平垂直于倾向单位向量

，通过向量计算得到

为

(3)

通过向量运算得到节理裂隙中心点至A、B、C、D四个点的向量为

(4)

A点坐标计算公式为

(5)

式中：x_a、y_a、z_a为向量

的x、y、z的值，x_A、y_A、z_A为点A的坐标。

同理求出B、C、D的坐标。并通过A、B、C、D及裂隙中心点r五个点确定节理裂隙面公式

(6)

式中：a、b、c、d为节理裂隙面方程参数。

采用三角投影法对单元节点进行判别是否在节理裂隙面上。点i、j、k为节理裂隙面上的任意点，p点为单元节点P在节理面上的投影点。三维空间内任意一点P（x₁+y₁，z₁）在平面S上的投影点

的计算公式为

(7)

若Δijp、Δikp、Δjkp面积之和等于Δijk，则表示P点位于节理裂隙面范围内。设置容差μ和R，对位于节理裂隙面范围内节点进一步筛选。P点到平面垂直距离d大于容差μ或p点到节理裂隙中心点x距离R大于

时，该节点不在节理裂隙平面内。

将筛选出的节点进行块体匹配，形成节理裂隙面块体组，修改其材料特征，编译并输出K文件，构建精细化节理岩体数值模型。

2　精细化裂隙岩体数值模型构建实例

收起

以云南省某石灰岩露天矿山为工程背景，对节理裂隙分布规律进行统计分析，按照节理裂隙分布函数对数值模型解析重组，构建精细化节理岩体地质模型。

2.1　节理裂隙分布规律求解方法

测线法、岩心取样法和钻孔成像法等方法采集到的节理裂隙分布信息为平面内的节理裂隙分布特征，不能完整描述节理裂隙三维空间分布特征。对于倾向、倾角相互独立的天然节理裂隙，可采用三维分形函数求解得到其三维分布累积概率^[14,15]，通过函数逼近求解方法得到最优分布函数及参数。岩层出露断面的节理裂隙迹长与节理裂隙延伸长度存在紧密联系^[16]，三维空间内，节理裂隙可近似看成圆形^[17]，圆盘直径即为延伸长度，可从测线法得到的节理裂隙迹长分布函数进行推算。

2.2　节理裂隙分布规律

采用测线法及钻孔电视成像法进行节理裂隙信息采集，测线法采集总长度72 m，钻孔电视成像法采集总长度378 m。台阶边坡坡面及节理裂隙分布如图3（a）所示，钻孔成像采集到的节理裂隙如图3（b）所示；将倾向和倾角这两个数组进行Pearson卡方检验（χ²），显著性概率为0.819，大于显著水平（0.05），即节理裂隙倾向和倾角的分布相互独立。通过三维分形函数求解得到各自三维分布累积概率^[14]，采用Kolmogorov-Smirnov逼近方法得到最优分布函数及参数。通过延伸长度与迹长之间的关系式求解得到节理裂隙延伸长度分布函数^[17]；三维空间内节理裂隙分布规律为：倾向β~N（211.1，17²），倾角α~N（75.5，13²），间距γ~N（1.1，5²），迹长l~F（2，1.52），延伸长度

。

2.3　精细化节理裂隙岩体地质模型

按照节理裂隙倾向、倾角、间距、延伸长度分布函数将地质模型进行解析重组，得到节理裂隙岩体精细化地质模型（图3c），地质模型尺寸为24 m×13 m×13.5 m；图3（a）为实际台阶边坡坡面，测线长度24 m，台阶坡面出露节理裂隙24条，节理裂隙密度为1~2条/m，图3（d）为地质模型台阶坡面，地质模型含68条节理裂隙，坡面上出露节理裂隙27条，节理裂隙密度为1~2条/m。台阶坡面与地质模型坡面均存在多条相互平行、交叉的节理裂隙，两者在数量、节理裂隙密度、分布形式等特征上接近；将测线法及解析重组法得到的节理裂隙产状信息绘制成等值线云图（图4a、b），两者节理裂隙的优势倾角、倾向以及分布的离散性接近。绘制重影对比图（如图4c），投影面积重叠率为87.2%。

3　基于精细化裂隙岩体模型的爆破数值模拟

收起

3.1　裂隙岩体爆破数值模型

3.1.1　爆破数值模型

结合矿山现场，在地质模型外围增加厚20 m的岩体，模拟半无限岩体中的爆破条件。在地质模型中构建单排三个垂直孔，钻孔直径110 mm，孔距为4.3 m，台阶高度12 m，钻孔超深1.5 m，装药长度10.5 m，堵塞长度3 m，底盘抵抗线5.3 m，逐孔起爆，数值模型单元网格尺寸为100 mm，台阶爆破模型如图5所示。设立不加载节理裂隙的完整岩体和优势倾角为75°的节理岩体两个工况，节理岩体的台阶坡面节理裂隙分布如图3（d）所示，台阶横截面节理裂隙分布如图6所示。通过不同截面的损伤分布对比，研究节理裂隙分布对爆破损伤扩展规律的影响，a~h为截面。截面的布置如图5所示。

3.1.2　材料本构模型

岩体选用Riedel-Hiermaier-Thoma（RHT）本构模型^[18]，其基本物理力学参数由试验所得。炸药为膨化硝铵炸药，采用Jones-Wilkens-Lee（JWL）状态方程描述。节理裂隙采用*MAT_PLASTIC_KINEMATIC材料本构模型^[7,19]，本构模型及材料参数如表1所示。

3.2　数值模拟结果及分析

从不同截面的损伤分布图（图7a、b）对比可知，钻孔轴线损伤范围平均长度为15.53 m，损伤轮廓线较规整，大致与台阶坡面平行；节理岩体的损伤范围平均长度为16.36 m，但是损伤轮廓线不规则。在台阶表面，节理岩体的损伤分布更为密集，损伤带分割成的块体大小更均匀。在岩体内部，节理裂隙岩体爆破过程中钻孔周边损伤破碎圈范围减小，完整岩体中损伤破碎圈直径为2.73 m，约为钻孔直径的25倍；节理岩体损伤破碎圈直径为1.85 m，约为钻孔直径的17倍，缩小了32.2%。从截面b、d、f的对比可知，在节理裂隙表面产生明显的损伤增大，局部破碎效果增强；同时，在节理裂隙面两侧的损伤分布密度存在较大的落差，一侧为片状分布的损伤区域，另一侧则是带状分布的损伤区域。完整岩体的损伤分布以钻孔为中心，损伤分布密度随着距钻孔的距离增大而逐渐减小，由片状的损伤分布逐渐变成带状损伤分布，即损伤后的岩体由破碎的小块逐渐向大块变化，粒径逐渐增大；而节理岩体的损伤分布较为离散，钻孔周边压缩损伤范围减小，岩体损伤分布范围增大，整体的损伤呈现片状与带状交织分布的规律，即损伤后的岩体大块与小块交错分布。

Chang-ping Yi进行爆破数值模拟块度统计分析时^[20]，将岩体划分为破碎岩体和未破碎岩体，当岩体损伤值D≥0.6时，岩体完全破碎，破碎岩体产生的块体其粒径小于数值模拟中的单元网格尺寸。数值模型的单元网格尺寸为100 mm，即可用损伤（D≥0.6）体积占比来表述爆破后粒径处于0~100 mm的体积占比。以损伤轮廓线、台阶轮廓线为界进行爆破体积统计，完整岩体和节理岩体的爆破体积分别为1193.01 m³、1336.61 m³，损伤（D≥0.6）体积分别为526 m³和481 m³，得到整个爆破区域碎块粒径处于0~100 mm的占比分别为44.09%、35.98%。节理岩体爆破时，爆破体积增加了12.04%，而0~100 mm的碎块占比较完整岩体降低了8.11%。

4　裂隙岩体爆破试验

收起

在矿山1775平台进行爆破试验，台阶坡面倾角75°，台阶高度12 m，节理裂隙分布规律如2-（2）节所述。采用垂直孔爆破方案，钻孔直径110 mm，孔距、排距为4.8 m×4.1 m，逐孔起爆方案，孔间延期50 ms，排间延期110 ms。爆破效果如图8所示，爆堆表面块度分布较为均匀，块度适中；而爆堆内部出现局部破碎效果增强，形成较多细小块体，同时在细小块体周围出现明显大块。表明，节理岩体爆破过程中，节理裂隙面两侧的破碎效果存在较大的落差。采用两个直径18.5 cm的皮球进行尺寸标定，通过Split-Desktop软件获得爆堆粒径分布，如表2所示，粒径在0~100 mm之间的碎块占比为37.68%。

对比现场爆破试验和数值模拟结果可知，损伤分布规律与爆破块度分布规律非常接近。根据数值模拟损伤分析，统计得到的0~100 mm的块体占比为35.98%，与现场爆破试验得到0~100 mm的块体占比为37.68%相比，误差为4.16%。

5　结论

收起

（1）通过采样得到的节理裂隙分布特征与解析重组法构建地质模型时得到的节理裂隙分布特征进行对比，二者倾向、倾角、间距、坡面出露节理裂隙条数、分布形式等特征接近，节理裂隙分布等值线云图面积重叠率达到87.2%。解析重组法生成的节理裂隙三维模型接近实际，操作方法可行、简便。

（2）节理裂隙面对岩体损伤破坏存在引导、控制作用，与完整岩体相比，爆破体积提高了12.04%，爆堆粒径处于0~100 mm的占比降低了8.11%。与现场爆破破碎效果相比，爆破模拟损伤分析得到粒径0~100 mm的占比误差为4.16%。解析重组法构建的三维节理裂隙数值模拟结果接近实际。

基金

收起

国家重点研发计划(2020YFC1909602; 2021YFC2902901)
湖北省重点研发计划(2021BCA152)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

杨超, 蒋淏南, 岳六一, 等. 法向卸荷条件下非贯通节理岩体剪切力学特性试验研究[J]. 岩土工程学报, 2022, 44(9): 1741-1750, 11.

YANG

Chao

, JIANG

Hao-nan

, YUE

Liu-yi

, et al. Shear mechanical properties of rock mass with discontinuous joints under unloading normal stress[J], Chinese Journal of Geotechnical Engineering, 2022, 44(9): 1741-1750, 11. (in Chinese)

[2]

李地元, 韩震宇, 孙小磊, 等. 含预制裂隙大理岩SHPB动态力学破坏特性试验研究[J]. 岩石力学与工程学报, 2017, 36(12): 2872-2883.

Di-yuan

, HAN

Zhen-yu

, SUN

Xiao-lei

, et al. Characteristics of dynamic failure of marble with artificial flaws under split hopkinson pressure bar tests[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(12): 2872-2883. (in Chinese)

[3]

王卫华, 李坤, 王小金, 等. SHPB加载下含不同倾角裂隙的类岩石试样力学特性[J]. 科技导报, 2016, 34(18): 246-250.

WANG

Wei-hua

, LI

Kun

, WANG

Xiao-jin

, et al. Mechanical properties of rock like specimens with different Inclined fractures under SHPB loading[J]. Science & Technology Review, 2016, 34(18): 246-250. (in Chinese)

[4]

宋杰, 孟凡震, 岳祝凤, 等. 节理裂隙岩体力学特性的研究进展[J]. 青岛理工大学学报, 2021, 42(5): 140-150, 154.

SONG

Jie

, MENG

Fan-zhen

, YUE

Zhu-feng

, et al. Research progress in mechanical properties of jointed and fractured rock mass[J]. Journal of Qingdao University of Technology, 2021, 42(5): 140-150, 154. (in Chinese)

[5]

金李, 卢文波, 陈明, 等. 节理岩体的爆破松动机理[J]. 爆炸与冲击, 2009, 29(5): 474-480.

JIN

, LU

Wen-Po

, CHEN

Ming

, et al. Blasting loosening mechanism of jointed rock mass[J]. Explosion and Shock Waves, 2009, 29(5): 474-480. (in Chinese)

[6]

马小满, 陈明, 严鹏, 等. 爆破作用下节理岩体应变能变化规律及影响研究[J]. 爆破, 2018, 35(3): 1-6, 67.

Xiao-man

, CHEN

Ming

, YAN

Peng

, et al. Study on strain energy processes and influence on jointed rock mass under blasting excavation[J]. Blasting, 2018, 35(3): 1-6, 67. (in Chinese)

[7]

周文海, 胡才智, 包娟, 等. 含节理岩体爆破过程中应力波传播与裂纹扩展的数值研究[J]. 力学学报, 2022, 54(9): 2501-2512.

ZHOU

Wen-hai

, HU

Cai-zhi

, BAO

Juan

, et al. Numerical study on crack propagation and stress wave propagation during blasting of rock mass[J], Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(9): 2501-2512. (in Chinese)

[8]

叶海旺, 王皓永, 雷涛, 等. 基于孔间延时优化的骨料用石灰岩爆破粉矿率控制[J]. 爆破, 2019, 36(4): 43-48, 68.

Hai-wang

, WANG

Hao-yong

, LEI

Tao

, et al. Rock fragmentation distribution control of limestone blasting for aggregate based on optimization of delay time between blasting hole[J]. Blasting, 2019, 36(4): 43-48, 68. (in Chinese)

[9]

王建超, 梁卫龙. 老龙山石灰岩矿深孔爆破参数优化试验[J]. 爆破, 2022, 39(4): 80-84, 99.

Li-hui

, LI

Hong-Wei

, YANG

Sai-Qun

, et al. Experimental study on influence of weak interlayer on rock blasting crack propagation[J]. Blasting, 2022, 39(4): 80-84, 99. (in Chinese)

[10]

WANG

Z L

, KONIETZKY

. Modelling of blast-induced fractures in jointed rock masses[J]. Engineering Fracture Mechanics, 2009, 76(12): 1945-1955.

[11]

唐海, 马谕杰, 夏祥, 等. 不同位置爆源作用下硐室围岩损伤规律的数值模拟研究[J/OL]. 爆破: 1-13.http://kns.cnki.net/kcms/detail/42.1164.TJ.20221115.1531.020.html.

TANG

Hai

, MA

Yu-jie

, XIA

Xiang

, et al. Numerical simulation study on damage law of chamber surrounding rock under different blasting sources[J/OL]. Blasting: 1-13.http://kns.cnki.net/kcms/detail/42.1164.TJ.20221115.1531.020.html. (in Chinese)

[12]

汤云坤. 充填节理岩体中应力波传播特性的颗粒流离散元数值模拟研究[D]. 西安: 长安大学, 2021.

TANG

Yun-kun

. Numerical simulation of stress wave propagation characteristics in filled jointed rock masses by particle flow discrete element method[D]. Xi'an: Chang'an University, 2021. (in Chinese)

[13]

高启栋, 卢文波, 冷振东, 等. 岩石爆破中孔内起爆位置对爆炸能量传输的调控作用研究[J]. 岩土工程学报, 2020, 42(11): 2050-2058.

GAO

Qi-dong

, LU

Wen-bo

, LENG

Zhen-dong

, et al. Regulating effect of detonator location in blast-holes on transmission of explosion energy in rock blasting[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(11): 2050-2058. (in Chinese)

[14]

赵萌, 唐辉明, 詹红兵, 等. 求解岩体裂隙产状三维概率分布的数值方法[J]. 地球科学, 2022, 47(4): 1470-1482.

ZHAO

Meng

, TANG

Hui-ming

, ZHAN

Hong-bing

, et al. A Numerical method for solving the three-dimensional probability distribution of rock mass fracture orientations[J]. Earth Science, 2020, 47(4): 1470-1482. (in Chinese)

[15]

OLIVIER

Fouché

, JEAN

Diebol

. Describing the geometry of 3d fracture systems by correcting for linear sampling bias[J]. Mathematical Geology, 2004, Vol.36(1): 33-63.

[16]

梅涛, 杨春和, 王贵宾, 等. 岩体节理大小估算新方法研究[J]. 岩石力学与工程学报, 2008(S2): 3503-3508.

MEI

Tao

, YANG

Chun-he

, WANG

Gui-bin

. et, al. A new method for estimating rock joint size[J]. Journal of Rock Mechanics and Engineering.2008(S2): 3503-3508. (in Chinese)

[17]

汪进超, 王川婴, 唐新建, 等. 王益腾. 基于钻孔摄像技术的岩体节理大小估算方法[J]. 岩土力学, 2017, 38(9): 2701-2707.

WANG

Jin-chao

, WANG

Chuan-ying

, TANG

Xin-jian

, et, al. A method for estimating rock mass joint size using borehole camera technique[J]. Rock and soil mechanics, 2017, 38(9): 2701-2707. (in Chinese)

[18]

HUO,

Xiao-feng

, SHI,

Xiu-zhi

, QIU,

Xian-yang

, et al. A study on raise blasting and blast-induced vibrations in highly stressed rock masses[J]. Tunneling & Underground Space Technology, 2022, Vol.123: 1-15.

[19]

袁尔君. 灰岩矿台阶爆破孔间延时对粉矿率影响及应用研究[D]. 武汉: 武汉理工大学, 2020.

YUAN

Er-jun

. Influence of delay time between holes interval of limestone mine bench blasting on ore powder rate and its application[D]. Wuhan: Wuhan University of Technology, 2020. (in Chinese)

[20]

Chang-ping

, SJÖBERG

Jonny

, JOHANSSON

Daniel

. Numerical modelling for blast-induced fragmentation in sublevel caving mines[J]. Tunneling & Underground Space Technology, 2017, Vol.68(1): 167-173.

2023年第40卷第4期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.3963/j.issn.1001-487X.2023.04.006

接收时间：2022-11-24
首发时间：2026-03-19
出版时间：2023-12-01

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2022-11-24

基金

National Key Research and Development Program(2020YFC1909602; 2021YFC2902901)

国家重点研发计划(2020YFC1909602; 2021YFC2902901)

Key R&D Program of Hubei Province(2021BCA152)

湖北省重点研发计划(2021BCA152)

作者信息

^1.武汉理工大学　资源与环境工程学院，武汉　430070

^2.武汉市公安局　治安管理局，武汉　430077

^3.华新水泥股份有限公司，武汉　430200

^4.五矿勘查开发有限公司，北京　100044

通讯作者:

雷涛（1983-），男，博士、讲师，主要从事采矿、安全、数字化矿山等方面的教学和科研工作，（E-mail）leitao539@163.com。

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/bp/CN/10.3963/j.issn.1001-487X.2023.04.006

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

RHT本构模型	ρ₀/（kg·m^-3） 2730	f_c/MPa 88.81	G/GPa 27.77	0.0518	0.083	β_C 0.032	β_t 0.036
膨化硝铵炸药本构模型	931	D/（m·s^-1） 3200	P_CJ/GPa 5.9	A 49.4	B 1.89	R₁ 3.9	R₂ 1.11
节理裂隙本构模型	1800	E₀/GPa 28	μ 0.241	σ/GPa 0.25	E_tan/MPa 0.025	β 0.5

RHT本构模型	ρ₀/（kg·m^-3） 2730	f_c/MPa 88.81	G/GPa 27.77	0.0518	0.083	β_C 0.032	β_t 0.036
膨化硝铵炸药本构模型	931	D/（m·s^-1） 3200	P_CJ/GPa 5.9	A 49.4	B 1.89	R₁ 3.9	R₂ 1.11
节理裂隙本构模型	1800	E₀/GPa 28	μ 0.241	σ/GPa 0.25	E_tan/MPa 0.025	β 0.5

RHT本构模型	ρ₀/（kg·m^-3） 2730	f_c/MPa 88.81	G/GPa 27.77	0.0518	0.083	β_C 0.032	β_t 0.036
膨化硝铵炸药本构模型	931	D/（m·s^-1） 3200	P_CJ/GPa 5.9	A 49.4	B 1.89	R₁ 3.9	R₂ 1.11
节理裂隙本构模型	1800	E₀/GPa 28	μ 0.241	σ/GPa 0.25	E_tan/MPa 0.025	β 0.5