爆破

参数名称	颗粒	岩体壁面	平行键参数
密度/（kg·m^-3）	2650	2650	法向刚度/（GPa·m^-1）	30.52
杨氏模量/GPa	44.09	44.09	切向刚度/（GPa·m^-1）	12.208
泊松比	0.35	0.35	最大法向应力/MPa	15.76
恢复系数静摩擦系数	0.5	0.5
滚动摩擦系数	0.01	0.01

参数名称	颗粒	岩体壁面	平行键参数
密度/（kg·m^-3）	2650	2650	法向刚度/（GPa·m^-1）	30.52
杨氏模量/GPa	44.09	44.09	切向刚度/（GPa·m^-1）	12.208
泊松比	0.35	0.35	最大法向应力/MPa	15.76
恢复系数静摩擦系数	0.5	0.5
滚动摩擦系数	0.01	0.01

基于Fluent-EDEM流固耦合算法的水下台阶爆破爆堆形成过程研究

PDF下载

吴亮 ¹ , 余创 ¹ , 梁志坚 ¹ , 段卫东 ¹ , 陈明 ²

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

收起

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

基于Fluent-EDEM流固耦合算法的水下台阶爆破爆堆形成过程研究

全屏

吴亮¹, 余创¹, 梁志坚¹, 段卫东¹, 陈明²

作者信息

^1.武汉科技大学　理学院　湖北省智能爆破工程技术研究中心，武汉　430065

^2.武汉大学　水工岩石力学教育部重点实验室，武汉　430072

吴亮（1980-），男，博士、副教授，主要从事工程爆破与岩石动力学研究方面的工作，（E-mail）wuliangwust@sina.com。

WU Liang (1980-), male, Ph. D, associate professor, mainly engaged in research and teaching of engineering blasting and rock dynamics, (E-mail) wuliangwust@sina.com.

Study on Formation Process of Underwater Bench Blasting Pile based on Fluid-structure Coupling Algorithm of Fluent-EDEM

Liang WU¹, Chuang YU¹, Zhi-jian LIANG¹, Wei-dong DUAN¹, Ming CHEN²

Affiliations

^1.Engineering Technology Research Center in Intelligent Blasting of Hubei Province, College of Science, Wuhan University of Science and Technology, Wuhan 430065, China

^2.Key Laboratory of Rock Mechanics in Hydraulic Structural Engineering Ministry of Education, Wuhan University, Wuhan 430072, China

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

文章导航

摘要

收起

针对水下岩石台阶爆破过程中破碎块体复杂的受力运动及其堆积过程，采用FLUENT-EDEM流固耦合的数值计算方法，在颗粒间引入平行键与Hertz-Mindlin接触模型，通过巴西劈裂实验与数值计算的结果对比，标定了平行键的参数，并建立了水下台阶爆破模型，分析了水对岩体破碎和运动的力学作用机理，研究了水下台阶块体颗粒鼓包抛掷的整个运动过程，讨论了块体速度、块体分层堆积和流场特点。结果表明：块体颗粒堆积分层与其初始速度以及位置高度有关，台阶中部块体颗粒相对位置较高，在爆炸荷载作用下水平运动距离最远，上部堵塞部分的块体颗粒受爆炸荷载作用小，其主要在重力与浮力作用下翻滚覆盖在中部块体颗粒层之上。起爆初期爆炸能量主要是破碎岩石，水体无明显流速，随着块体颗粒的鼓包抛掷运动，台阶前表面附近的流体流线速度最大，当块体颗粒开始滚动滑落时，流线则主要受到颗粒的运动形式影响。另外，FLUENT-EDEM流固耦合的数值计算模型能精确模拟水下台阶爆破块体运动问题，可为水下岩体爆破工程提供新的研究手段。

关键词

水下爆破 / 流固耦合 / FLUENT / EDEM / 爆堆

Abstract

收起

To study the complex mechanical movement and accumulation process of broken blocks in the process of underwater rock bench blasting, a fluid-solid coupling numerical calculation method based on FLUENT-EDEM was adopted. Firstly, the parallel bond and Hertz-Mindlin contact model were introduced between the particles, and the parameters of the parallel bond were calibrated by comparing the results of Brazilian splitting experiment and numerical calculation. Furthermore, the underwater bench blasting model was established, and the mechanical action mechanism of water on rock mass breaking and movement was analyzed. Finally, the whole ejection process of particles from a submerged step block was studied, and characteristics of the block velocity, block stacking and flow field were discussed. The results show that the accumulation and stratification of block particles are related to their initial velocity and height, and the horizontal movement distance of the block particles in the middle step is the biggest under the explosion load action. Meanwhile, the block particles in the upper blocking part mainly roll and cover the middle block particle layer under the action of gravity and buoyancy, which are less affected by the explosion load. Additionally, the explosion energy is mainly used to break the rock at the initial detonation stage, and the water has no obvious flow velocity. With the bulging and throwing movement of block particles, the fluid streamline velocity near the front surface of the bench is largest. When the block particles start to roll and slide, the streamline is mainly affected by the movement form of the particles. Therefore, the fluid-solid coupling model by Fluent-EDEM can accurately simulate the block motion in the underwater bench blasting, which can provide a new research method for underwater rock blasting engineering.

Key words

underwater blasting / fluid solid coupling / FLUENT / EDEM / blasting pile

引用本文

吴亮, 余创, 梁志坚, 段卫东, 陈明. 基于Fluent-EDEM流固耦合算法的水下台阶爆破爆堆形成过程研究. 爆破, 2023 , 40 (4) : 37 -43 . DOI: 10.3963/j.issn.1001-487X.2023.04.005

Liang WU, Chuang YU, Zhi-jian LIANG, Wei-dong DUAN, Ming CHEN. Study on Formation Process of Underwater Bench Blasting Pile based on Fluid-structure Coupling Algorithm of Fluent-EDEM[J]. Blasting, 2023 , 40 (4) : 37 -43 . DOI: 10.3963/j.issn.1001-487X.2023.04.005

正文

收起

众所周知，岩石的破碎和运动是岩土爆破中的2个关键环节，其中破碎块体的运动常被称抛掷运动。当实际生产中矿岩与普通岩互为夹层时，矿体运动结果对分离矿石和废石的影响较大，有时对采矿效率的影响远大于机械设备^[1]。水下钻孔爆破作为工程爆破技术中的重要组成部分，广泛应用在扩展港口、疏浚航道、取水以及抢险救灾等工程中。水下与陆地岩石爆破的主要差异在于岩石表面接触介质的不同，直接影响了岩石的爆破破碎以及块体抛掷运动效果^[2,3]。由于影响水下爆破效果的因素错综复杂，很多工程只能凭有关经验公式与试验数据估算，往往没有达到破碎与堆积效果；另外，室内的模型实验成本高且周期长，而现场试验无法及时获取水下爆破破碎与堆积效果，且很多工程也不具备现场试验条件。因此，水下爆破破碎与堆积效果的预测工作亟待解决。

由于爆破作用过程的复杂性，且现有数学模型还不完善或存在局限性，目前工程中对爆破效果的预测还主要依靠经验判断，其准确性及可靠性难以满足要求^[4]。数十年来，数值模拟技术在爆破领域得到了广泛应用，但传统的连续介质模型很难对岩石爆炸冲击破碎后的块体运动进行模拟。YI采用颗粒模型（PBM）模拟了爆轰过程以及被爆的岩石等脆性介质^[5]，实现岩石爆破破碎过程，突破了连续介质模型难处理大变形的问题。PREECE结合并行计算技术采用三维离散元模型分析了数百个炮孔^[6]，离散单元数超过100万。安华明等采用有限—离散杂交元模拟了爆破漏斗产生的整个动态过程^[7]。余德运等采用流体与DDA块体之间网格覆盖的流固耦合技术模拟了爆破漏斗形成过程^[8]。苏都都等采用PFC 2D对爆破参数与爆堆形态的关系进行了研究^[4]，预测了台阶爆破的爆堆形态。冷振东等在岩体强度参数服从Weibull分布规律的基础上^[9]，采用3DEC离散元软件建立了台阶爆破的三维台阶爆破模型，模拟起爆位置对岩石的动态破碎与抛掷过程的影响规律。周旺潇等和YAN等在考虑不同区域块度尺寸分布的基础上^[10,11]，采用3DEC离散元模型讨论了的爆破漏斗与台阶爆破爆堆形态。可见，预测与控制爆破抛掷运动的数值软件的应用极大地提高了生产效率，但水下岩石爆破破碎及其抛掷运动目前鲜见报道。

近年来研究者开始采用CFD-DEM方法解决工程中的流固耦合问题，例如水下岩塞爆破和深海采矿管道固液混合流体^[12,13]。XU等基于Fluent-EDEM耦合方法模拟了管道堵塞过程^[14]，通过考虑连续流体相和离散颗粒相，初步研究了水合物的堵塞机理。HE等基于Fluent-EDEM引入滑移网格^[15]，该模型可以同时分析大型实体模型自由移动引起的流场和小颗粒的运动。FU等采用CFD-DEM耦合的计算方法模拟了渗流力与掌子面岩体颗粒的相互作用^[16]。WU等基于DEM中颗粒的接触模型^[17,18]，采用Fluent-EDEM耦合方法研究水下台阶爆破块体的运动过程及爆堆形态。随着计算软件的开发与应用，越来越多的研究者尝试采用各种数值计算的方法解决水下爆破的问题。因此，采用基于UDF开发的Fluent-EDEM耦合方法，考虑颗粒间的连接作用来研究水下岩石台阶爆破破碎及块体运动的问题，期望研究成果能促进水下台阶爆破技术发展。

1　计算方法

收起

1.1　流体基本原理

计算流体力学是利用数学的计算方法来实现对流场中的控制方程进行计算，并对有限元网格节点进行数值求解。当研究对象为不可压缩流体时，基本的定律有三个：质量守恒定律、动量守恒定律、能量守恒定律，如下式（1）和式（2）。在本文中没有能量的交换，故不考虑能量守恒方程。

(1)

(2)

式中：ρ为流体密度；t表示时间；

为流体速度矢量，右下标x、y、z分别表示速度矢量在X、Y、Z轴上的投影；∇为哈密顿算子；p表示流体介质压力；表示

加速度矢量；

为颗粒与流体的相互作用力；

为应力张量，表示流体的粘滞力。

1.2　DEM颗粒接触模型

DEM模型中颗粒受力运动遵循牛顿第二定律，包括了转动和平动，颗粒的运动方程（动量守恒方程）如下式（3）和（4）所示。

(3)

(4)

式中：v、ω、m分别为颗粒速度、角速度和质量，下标i，j表示两个相互接触的颗粒；F_ij为i和j颗粒间的接触力；F_g为颗粒间的引力；F_f为流体与颗粒间的相互作用力；I是颗粒转动惯量；T是颗粒间的力矩；T_f是流体引起的扭矩。

本文采用Hertz-Mindlin接触模型计算DEM颗粒之间的作用力，该模型具备滚动摩擦特性，如图1所示。Hertz-Mindlin接触模型的计算表达式见文献[19]。

1.3　DEM颗粒平行键模型

平行键模型（Parallel bond Model，PBM）可以用于模拟微观结构如何影响宏观行为，它将固体介质离散为具有固定尺寸的胶结颗粒材料，由于颗粒之间存在空隙，在固体介质模型形成平行键时，这些颗粒间的空隙可以视为混凝土材料中的自然孔隙或岩石中的自然裂缝，这些裂缝相互作用、合并形成宏观裂纹。平行键模型能表达颗粒间法向和切向的作用力，如图2所示，当法向或切向应力超过最大许可应力时平行键断裂，颗粒之间的相互作用回归到Hertz-Mindlin接触模型。

在平行键生成之后颗粒上的力和扭矩被重置为0，在模型开始计算时则根据公式（5）~（8）在每个时间步长中计算平行键受到的力和扭矩的增量。

(5)

(6)

(7)

(8)

式中：F_n和F_t为法向力和切向力；M_n和M_t则为法向扭矩和切向扭矩；v_n和v_t为颗粒的法向和切向速度；ω_n和ω_t为法向和切向角速度；S_n和S_t为法向刚度和切向刚度；δt为时间步长；A是平行键截面面积，计算公式如式（9）所示；J为平行键截面极惯性矩，计算公式见式（10）。

(9)

(10)

式中：R_B为平行键半径，是平行键合模型输入参数之一，取值一般大于颗粒半径。

当法向力或切向力超过其限定值时，颗粒间的平行键发生断裂

(11)

(12)

式中：σ_max、τ_max为最大法向应力和最大切向应力。

2　水下台阶爆破计算模型

收起

2.1　计算模型

参考实际台阶爆破工程，数值模型中上台阶面距离水面10 m、台阶面宽6 m，台阶高度为10 m，炮孔直径0.105 m，孔深10 m，装药长度为7 m，堵塞3 m，孔间距为4.5 m，最小抵抗线为3 m。本文数值模型以单段4个炮孔同时起爆为对象进行建模，如图3所示。图3（a）为Fluent模型，上表面为压力出口边界，压力值设置为一个大气压。图3（b）为EDEM模型，颗粒填充区域为长23.5 m、宽5.5 m、高10 m的长方体，颗粒采用直径0.21 m的球形颗粒，共填充146728颗。

2.2　模型参数

在EDEM模型中引入平行键的模型来模拟连续的岩体介质，平行键的模型除了包含摩擦系数在内的基本力学参数外还需要额外输入5个平行键参数。由于平行键的模型能体现微观上颗粒之间的力学行为，除了平行键半径取颗粒半径值外，至今的研究结果还不能明确地给出实验中获得的参数与平行键参数之间的关系式，因此为了获得准确的数值计算参数，需对平行键的参数进行标定。许多学者使用DEM进行巴西劈裂数值计算^[20-22]，通过与实验结果对比，标定平行键的参数。

在EDEM中建立直径50 mm、高100 mm的圆柱模型，投放的球形颗粒直径均为3 mm，共投放颗粒8473颗，共生成平行键29 258个。分别对圆柱颗粒模型进行巴西劈裂数值模拟，模型破坏结果如图4（a）所示，应力-应变曲线如图4（b）。

图4（a）中灰白色为颗粒间的平行键，红色代表已破坏的平行键，巴西劈裂模型的破坏模式为拉伸破坏。图4（b）中三组巴西劈裂实验曲线来自文献[23]的实验。三组巴西劈裂实验的平均峰值应力为2.46 MPa，而数值模型的峰值应力为2.44 MPa。说明基于平行键方法建立的材料模型能够在一定程度上反映实际工况中岩石、混凝土等脆性材料的力学行为。

综上，结合标定后的平行键模型参数，将堆积的颗粒建立平行键，台阶模型共生成平行键351 460个，具体参数如下表1所示。需要注意的是，同一组模型中含有平行键与无平行键在时间步长的选取上存在差异，通常在进行含有平行键的数值模型计算时选取两者中时间步长较小的一个^[24]，本文离散元的计算步长为3×10^-5 s，流体的计算步长是离散元的10倍，为3×10^-4 s。

3　模拟结果

收起

3.1　台阶爆破破碎与块体运动过程

数值模型中通过在炮孔位置设置高压区域实现爆炸荷载的冲击作用，块体颗粒运动速度与未破坏的平行键分布图见图5，该图显示了在有平行键作用下水下台阶块体颗粒鼓包抛掷的整个运动过程。从台阶爆破破碎来看，初始状态颗粒间的平行键完整密布，炮孔内爆炸荷载作用在炮孔壁后，炮孔壁介质颗粒向外传递爆炸荷载并向外运动，当颗粒间的平行键超过临界值时将断开，处于炮孔附近的平行键大量破裂，结果显示平行键的断裂是从炮孔壁逐步向外开始扩散的。由于台阶颗粒排列是随机的，所以平行键断裂扩展到临水面时并不是一个平面，同时考虑应力波的反射作用，最先破裂的部位在炮孔的中部，该部位的块体颗粒获得的动能较多，所以抛掷的速度也最大，而台阶底部和上部的颗粒平行键的断裂滞后，特别是台阶上部以及两侧残存有大块。

从块体运动来看，在0.0204 s时刻，炮孔附近的块体颗粒速度明显大于两侧靠近岩体壁面的颗粒，说明鼓包运动由最小抵抗线方向开始。在0.492 s时刻，表面块体颗粒抛掷的速度峰值为2.48 m/s，随着颗粒进一步剥离表面，峰值速度持续增加。在0.918 s时，中部爆破区域呈现塌陷趋势，而两侧未见崩坏，部分颗粒在重力与浮力的作用下进行抛掷运动，并在3 s时刻颗粒抛掷速度峰值达到4.5 m/s，随后速度峰值逐渐衰减。在2.004 s时刻，可以明显观察到块体颗粒群形成了具有一定坡度的堆积体，此时块体颗粒的运动形式滚落为主。随着块体颗粒的继续滚动，块体颗粒在台阶前表面形成了平缓的爆堆，并伴随有左侧与堵塞段大块堆积在上面。

3.2　堆积范围和分层结果分析

为了方便观察台阶前方块体颗粒的堆积层次分布规律，将台阶前堆积的块体颗粒群选取并进行上色，底部颗粒为蓝色、中部颗粒为红色、上部颗粒为绿色，见图6（a）。图6（b）为水下台阶爆堆中间截面最终堆积形状的侧视图。通过侧视图测出爆堆堆积形状的休止角约为20°。爆后块体颗粒的分层延续了其在原台阶模型中的分布形式，块体爆堆中蓝色块体颗粒位于最下层、红色块体颗粒位于中间、绿色块体颗粒位于上层，且中部红色块体颗粒运动距离最远，中部红色块体颗粒覆盖了大部分的底部蓝色块体颗粒，主要原因是块体颗粒堆积分层与其初始速度和位置高度有关，台阶中部块体颗粒相对位置较高，在爆炸荷载作用下水平运动距离最远，而上部堵塞部分的块体颗粒受爆炸荷载作用小，其主要在重力与浮力作用下翻滚覆盖在中部块体颗粒层之上。

3.3　块体运动速度和流体流场分析

为更具体地分析块体颗粒的堆积规律，在台阶前表面中轴线上选取6个测点，测点间距为2 m，每个测点均包含若干块体颗粒，见图3（b）所示，绘制颗粒的平均速度随时间变化的时程图，见图7，该速度方向垂直于台阶前表面。在炸药反向起爆之后，炸药起爆点位置处的2#和3#测点水平速度峰值最大，分别为3.4 m/s和3.0 m/s；位于底部的1#测点受到颗粒下沉碰撞影响，其峰值速度分别为2.1 m/s；台阶堵塞段的5#和6#测点的水平速度峰值曲线加速段斜率较低且到达最大速度峰值出现较晚，说明台阶堵塞段块体颗粒受爆炸荷载作用的强度较小，另外也说明该部分块体滚动滑落运动明显。随着块体颗粒的运动，颗粒的速度受水的阻力而衰减，在鼓包运动结束之后，块体颗粒以滚动滑落的运动形式为主。结合爆堆分层规律，并观察速度时程曲线后半段可知，块体颗粒的滚动滑落速度与其所处位置有关，位于底部的1#测点的块体颗粒由于受到底面的摩擦阻碍以及上层颗粒的覆盖作用，其块体颗粒的滑动速度较小，而4#~6#测点的滚动滑落速度随着高度的增加而衰减较慢。

图8为台阶模型中轴面的流线分布图。起爆初期爆炸能量主要是破碎岩石，水体无明显流速，随着块体颗粒的鼓包抛掷运动，台阶前表面附近的流体流线速度最大。在2.0 s时刻台阶上部块体颗粒溃散滑落，导致了其表面附近出现流线回流。随着块体颗粒的滚动滑落形成具有平缓坡度的爆堆，流线的回流形式进一步扩大，但是流线的峰值速度在逐渐衰减。当块体颗粒开始滚动滑落时，流体流线则主要受到颗粒的运动形式影响，流线峰值速度区域均出现在块体颗粒附近。

4　结论

收起

基于Fluent-EDEM耦合方法并采用平行键与Hertz-Mindlin接触模型模拟了水下岩体破碎和运动的力学作用机理，讨论了爆破块体抛掷速度、块体分层堆积和爆区流场特点，揭示了水下岩体爆破破碎块体堆积的运动特点。初步结论如下：

（1）水下台阶爆破块体颗粒堆积分层与其初始速度和位置高度有关，其本质是爆破能量密度分布与块体重力势能的外在表现。

（2）起爆初期爆炸能量主要是破碎岩石，水体无明显流速，随着块体颗粒的鼓包抛掷运动，最小抵抗线方向流体流线速度最大。当块体颗粒开始滚动滑落时，流线则主要受到颗粒的运动形式影响。

（3）采用平行键与Hertz-Mindlin接触模型能反映岩体爆破破碎过程，揭示大块产生的机理与范围。计算结果表明FLUENT-EDEM流固耦合的数值模型能较准确地模拟水下台阶爆破岩体破碎与块体运动问题，可为水下岩体爆破工程提供新的研究手段。

基金

收起

国家自然科学基金项目(51779193; 51979205)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

吴亮, 鲁帅, 许锋, 等. 矿岩爆破破碎机理、块度分布与测量技术研究动态[J]. 金属矿山, 2016(7): 47-53.

Liang

, LU

Shuai

, XU

Feng

, et al. Research trends about the mechanism of ore rock blasting fragmentation, the distribution of block size and measurement technique[J]. Metal Mine, 2016(7): 47-53. (in Chinese)

[2]

刘美山, 吴新霞, 张恒伟. 混凝土水下爆破炸药单耗试验分析[J]. 爆破, 2007(1): 10-13, 20.

LIU

Mei-shan

, WU

Xin-xia

, ZHANG

Heng-wei

. Experimental analysis on specific charge of underwater explosion of concrete[J]. Blasting, 2007(1): 10-13, 20. (in Chinese)

[3]

赵根. 深水条件下围堰拆除爆破技术研究[D]. 合肥: 中国科学技术大学, 2008.

ZHAO

Gen

. Study on technology of the cofferdam demolition blasting in deep water conditions[D]. Hefei: University of Science and Technology of China, 2008. (in Chinese)

[4]

苏都都, 严鹏, 卢文波, 等. 露天台阶爆破爆堆形态的PFC模拟[J]. 爆破, 2012, 29(3): 35-41.

Du-du

, YAN

Peng

, LU

Wen-bo

, et al. Prediction of muckpile profile for open bench blasting with PFC[J]. Blasting, 2012, 29(3): 35-41. (in Chinese)

[5]

, JOHANSSON

. Discrete element modelling of blast fragmentation of a mortar cylinder[C]//11th International Symposium on Rock Fragmentation by Blasting. Sydney: [sn], 2015: 793-800.

[6]

PREECE

D S

, TAWADROUS

, SILLING

S A

, et al. Modelling full-scale blast heave with three-dimensional distinct elements and parallel processing[C]//11th International Symposium on Rock Fragmentation by Blasting. Sydney: [sn], 2015: 127-136

[7]

安华明, LIU

H-Y

, 汪旭光, 等. 有限—离散杂交元岩石爆破漏斗数值模拟研究[J]. 爆破, 2016, 33(4): 1-8, 21.

Hua-ming

, LIU

H-Y

, WANG

Xu-guang

, et al. Hybrid finite-discreet element method for modelling of formation process of blasting crater[J]. Blasting, 2016, 33(4): 1-8, 21. (in Chinese)

[8]

余德运, 刘殿书, 王庆斌, 等. 基于流固耦合加载技术的DDA方法[J]. 爆破, 2016, 33(1): 6-11.

De-yun

, LIU

Dian-shu

, WANG

Qing-bin

, et al. An improved DDA method based on fluid-structure coupling loading technology[J]. Blasting, 2016, 33(1): 6-11. (in Chinese)

[9]

冷振东, 刘亮, 周旺潇, 等. 起爆位置对台阶爆破爆堆形态影响的离散元分析[J]. 爆破, 2018, 35(2): 50-55, 100.

LENG

Zhen-dong

, LIU

Liang

, ZHOU

Wang-xiao

, et al. Numerical investigation of initiation points on muck pile profile in bench blasting[J]. Blasting, 2018, 35(2): 50-55, 100. (in Chinese)

[10]

周旺潇, 严鹏, 郑炳旭, 等. 爆破漏斗形成过程数值模拟的几个关键问题[J]. 爆破, 2014, 31(3): 15-22.

ZHOU

Wang-xiao

, YAN

Peng

, ZHENG

Bing-xu

, et al. Key problems in simulation of formation process of blasting crater[J]. Blasting, 2014, 31(3): 15-22. (in Chinese)

[11]

YAN

Peng

, ZHOU

Wang-xiao

, LU

Wen-bo

, et al. Simulation of bench blasting considering fragmentation size distribution[J]. International Journal of Impact Engineering, 2016(90): 132-145.

[12]

Liang

, LIANG

Zhi-jian

, CHEN

Ming

, et al. Experiments and fluent-engineering discrete element method-based numerical analysis of block motion in underwater rock-plug blasting[J]. Applied Sciences-Basel, 2023(13): 348.

[13]

DAI

, ZHANG

Yan-yang

, LI

Xu-yang

. Numerical and experimental investigations on pipeline internal solid-liquid mixed fluid for deep ocean mining[J]. Ocean Engineering, 2020, 220(2): 108411.

[14]

DUAN

, SHI

Bo-hui

, WANG

Jun-ao

, et al. Simulation of the hydrate blockage process in a water-dominated system via the CFD-DEM method[J]. Journal of Natural Gas Science and Engineering, 2021(96): 104241.

[15]

, ANDREW

E B

, ALI

. Coupling CFD-DEM with dynamic meshing: A new approach for fluid-structure interaction in particle-fluid flows[J]. Powder Technology, 2018, 325: 620-631.

[16]

Yan-bin

, ZENG

De-qi

, XIONG

Hao

, et al. Seepage effect on failure mechanisms of the underwater tunnel face via CFD DEM coupling[J]. Computers and Geotechnics, 2022, 146: 104591.

[17]

Liang

, LIANG

Zhi-jian

, CHEN

Ming

. Experiments and simulation of block motion in underwater bench blasting[J]. Scientific Reports, 2023, 13(1): 4703.

[18]

Liang

, LIANG

Zhi-jian

, ZHANG

, CHEN

Ming

. The blasting pile formation process of rock mass bench blasting under water based on fluid solid coupling model of fluent and EDEM[C]//13# International Symposium on rock fragmentation by blasting, Hangzhou China, 2022: 539-547.

[19]

Xing-yue

, ZHAO

Ji-dong

. Dam-break of mixtures consisting of non-Newtonian liquids and granular particles[J]. Powder Technol, 2018, 338: 493-505.

[20]

WANG

Yu-han

, NGUYEN

. The effects of rock-infill interfacial properties on the compressive damage behaviour of flawed rocks: Results from a DEM study[J]. Theoretical and Applied Fracture Mechanics, 2022, 117: 103166.

[21]

Song-tao

, KARLOV

EK J

. Calibration and uniqueness analysis of microparameters for DEM cohesive granular material[J]. International Journal of Mining Science and Technology, 2022, 32(1): 121-136.

[22]

Zhi-jun

, ZHANG

Peng-lin

, FAN

Li-feng

, et al. Numerical study of the effect of confining pressure on the rock breakage efficiency and fragment size distribution of a TBM cutter using a coupled FEM-DEM method[J]. Tunnelling and Underground Space Technology, 2019, 88: 260-275.

[23]

WANG

, GAO

, JI

, et al. DEM analysis on the role of aggregates on concrete strength[J]. Computers and Geotechnics, 2020, 119: 103290.

[24]

WANG

Li-ge

, GE

Rui-huan

, CHEN

Xi-zhong

. Establishing an oblique impact breakage master curve using a DEM bonded contact model[J]. Computers and Geotechnics, 2020, 145: 104668.

2023年第40卷第4期

PDF下载

100

引用本文

BibTeX

文章信息

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

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

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2022-05-16

基金

National Natural Science Foundation of China(51779193; 51979205)

国家自然科学基金项目(51779193; 51979205)

作者信息

^1.武汉科技大学　理学院　湖北省智能爆破工程技术研究中心，武汉　430065

^2.武汉大学　水工岩石力学教育部重点实验室，武汉　430072

参考文献

分享链接

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

分享至

全文二维码

扫描看全文

引用本文

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