• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

基于FLUENT-EDEM耦合的爆炸抛掷特性研究

顾强 张世豪 安晓红 张亚

潘亚豪, 宗周红, 钱海敏, 黄杰, 单玉麟. 钙质砂介质中爆炸波传播规律的试验研究[J]. 爆炸与冲击, 2023, 43(5): 053201. doi: 10.11883/bzycj-2022-0117
引用本文: 顾强, 张世豪, 安晓红, 张亚. 基于FLUENT-EDEM耦合的爆炸抛掷特性研究[J]. 爆炸与冲击, 2016, 36(5): 611-616. doi: 10.11883/1001-1455(2016)05-0611-06
PAN Yahao, ZONG Zhouhong, QIAN Haimin, HUANG Jie, SHAN Yulin. Experimental study on blast wave propagation in calcareous sand[J]. Explosion And Shock Waves, 2023, 43(5): 053201. doi: 10.11883/bzycj-2022-0117
Citation: Gu Qiang, Zhang Shihao, An Xiaohong, Zhang Ya. Research on throw characteristics of rock debris explosion based on coupled FLUENT-EDEM[J]. Explosion And Shock Waves, 2016, 36(5): 611-616. doi: 10.11883/1001-1455(2016)05-0611-06

基于FLUENT-EDEM耦合的爆炸抛掷特性研究

doi: 10.11883/1001-1455(2016)05-0611-06
详细信息
    作者简介:

    顾强(1961—),男,硕士,教授,hg14gu@163.com

  • 中图分类号: O383

Research on throw characteristics of rock debris explosion based on coupled FLUENT-EDEM

  • 摘要: 为了研究弹丸起爆后岩土碎屑的抛掷特性, 针对爆炸后侵彻通道内复杂的气固两相流动问题,提出了一种FLUENT软件与EDEM软件联合的数值模拟方案,分析了不同深径比条件下钻地通道内介质碎屑的抛掷过程。进一步研究了介质碎屑运移特性的影响因素。结果表明:介质碎屑的运移能力与爆炸荷载峰值成正比,与钻地通道深径比及碎屑粒径成反比。
  • 钙质砂通常是指由海洋生物形成的、CaCO3含量超过50%的颗粒状材料,具有孔隙多、形状不规则、易破碎、易胶结等特点,其工程力学性质与一般陆源砂有明显差异[1-2]。目前,对陆源砂抗爆性能已有较为系统的研究,如美军的《TM5-855-1》手册提供了不同密实程度砂土中爆炸波传播的经验公式[3]。已有砂土爆炸动力学方面的研究主要集中在两个方面[4-7]:(1)通过总结和分析爆炸作用下砂土内部的动力反应,得出应力波在砂土中的传播和衰减规律;(2)通过模型试验和数值分析,开展不同试验参数情况下,饱和砂土的爆炸密实特性和机制研究。钙质砂作为一种特殊的岩土介质,其工程力学性质与普通的石英砂有显著的差别,当前对钙质砂低应变率下的力学性质已经有了较为充分的研究[8-10] ,对钙质砂中高应变率下的研究多集中在颗粒破碎和耗能[11-14]上。Lv等[15]对钙质砂与石英砂开展的三轴试验显示,钙质砂的体积蠕变和剪切蠕变均明显大于石英砂;对钙质砂开展的霍普金森压杆实验显示,钙质砂的体积模量、压缩强度、平均波速均明显小于石英砂,对应力波的衰减效果更显著[16-17],因此不能将当前石英砂的研究结果直接用于钙质砂。

    伴随着我国海洋发展战略的实施,人们对强动载下钙质砂的动力特性研究愈加感兴趣,而钙质砂的爆炸试验研究刚刚起步。徐学勇[18]通过饱和钙质砂和石英砂的爆炸密实试验对爆炸波特性和表面沉降规律进行了研究。赵章泳等[19-20]对钙质砂进行了平面装药和球形装药试验研究,给出了非饱和钙质砂中爆炸波的弹塑性波速和应力衰减系数。石晗[21]进行了钙质砂地基野外爆炸试验,分析了钙质砂中土压力和孔隙水压力的变化规律,并比较了不同浸水状态下钙质砂的衰减作用。

    本文中通过开展钙质砂和石英砂的野外爆炸试验,对比分析两种砂在地面爆炸下的爆坑、爆炸波在介质中的弹塑性波速、升压时间、压力和冲量峰值,并给出钙质砂和石英砂的爆炸波衰减系数,研究成果可供相关工程地基设计和加固作参考。

    试验所用钙质砂(图1(a))取自南海某岛,主要成分为CaCO3,外观偏黄,且含有贝壳和珊瑚碎块。用作对比试验的石英砂(图1(b))主要成分为SiO2,外观为白色。钙质砂和石英砂试样的级配曲线如图2所示,由级配曲线可得钙质砂平均粒径d50 = 0.48 mm,不均匀系数Cu = 9.43,石英砂平均粒径d50 = 0.40 mm,不均匀系数Cu = 5.08。试验当天对回填砂取样测得钙质砂的密度为1.32 g/cm3,含水率为3.5 %,石英砂的密度为1.39 g/cm3,含水率为2.3 %。

    图  1  砂土试样
    Figure  1.  Sand specimen
    图  2  试验用砂土级配曲线
    Figure  2.  Particle size distribution

    坑底底部埋设有管廊结构模型,上方填土高度为0.83 m、宽度为4.0 m、总长度为6.0 m,其中长度方向每2.0 m为一节段,分别分层回填钙质砂、石英砂和黄土(见图3),每层厚度约为15 ~20 cm,并人工夯实,边缘处均用彩条布进行分隔。本文试验在钙质砂和石英砂回填段进行。

    图  3  试验现场
    Figure  3.  Test site

    起爆点设在各段长度方向中心点,炸药采用军用TNT药块,平放于地面。在爆心正下方距离地面0.35、0.50、0.65、0.83 m处放置土压力传感器(见图4)。

    图  4  土压力传感器布置
    Figure  4.  Locations of pressure sensors

    试验采用土压力传感器,量程为5~25 MPa。传感器输出的信号由泰斯特测试公司的TST5910N动态信号测试分析系统采集,采用瞬态采样,采样频率为1 MHz。

    表1汇总了本次爆炸试验的实际参数,钙质砂和石英砂因装药形状不同,爆心距有微小差距。

    表  1  爆炸试验参数设计表
    Table  1.  Measured parameters of explosion test
    试验编号土体类型装药质量/kg比例距离/(m·kg−1/3)
    H = 0.35 mH = 0.50 mH = 0.65 m
    CS-1钙质砂0.20.640.901.15
    CS-20.40.520.720.93
    CS-30.80.450.610.78
    CS-41.60.330.460.59
    SS-1石英砂0.20.620.881.13
    SS-20.40.510.710.92
    SS-30.80.400.570.73
    SS-41.60.340.470.60
    SS-53.20.270.370.48
    下载: 导出CSV 
    | 显示表格

    为验证试验重复性及数据的可靠性,对装药0.2、0.8 kg的爆炸试验分别进行两次,图5为两次试验在同一测点的压力和冲量时程曲线。表2中计算了两次重复试验的相对差值,可见第2次峰值压力略微大于第1次试验,且爆炸波波阵面到达时间更早,压力上升更快,原因可能是第2次试验的砂土密实度较第1次试验的更高。炸药爆炸产生的能量会破坏砂土颗粒的初始结构,使砂土孔隙水压力升高,颗粒结构在自重或其他外荷载作用下排列成更紧密的状态,这一原理已广泛应用于砂土地基处理[5]。另外,考虑到钙质砂颗粒易碎、在常应力下就有颗粒破碎发生的特性[22],在爆炸前后分别进行了采样筛分试验,绘制出原始砂样和多次爆炸后砂样的颗粒级配曲线(见图6),爆炸后大粒径颗粒含量明显减少,平均粒径d50 = 0.42 mm,小于初始试样平均粒径d50 = 0.48 mm,表明钙质砂颗粒在冲击荷载作用下发生了破碎,对峰值压力的衰减效果减弱。

    表  2  相同爆炸药量同测点峰值压力对比
    Table  2.  Comparison of peak pressure at the same location
    W/kg试验H = 0.35 mH = 0.50 mH = 0.65 m
    峰值压力/MPa相对差值/%峰值压力/MPa相对差值/%峰值压力/MPa相对差值/%
    0.210.8631.070.4622.930.3419.49
    21.130.570.40
    0.813.991.871.653.591.072.49
    24.061.711.10
    注:相对差值 = (第2次峰值压力-第1次峰值压力)/第1次峰值压力×100 %。
    下载: 导出CSV 
    | 显示表格
    图  5  相同爆炸药量同测点压力和冲量时程曲线对比图
    Figure  5.  Comparison of pressure time-history and impulse time-history curves at the same location
    图  6  爆炸试验前后钙质砂颗粒级配曲线
    Figure  6.  Particle size distribution before and after the explosion

    图7为试验后钙质砂和石英砂表面形成的爆坑,俯视呈两阶同心圆,爆坑断面示意图及尺寸定义见图8。一般来说,地面爆炸后,部分抛掷物会落回中心爆坑,应移除后测量深度h[23]表3列出了钙质砂和石英砂在W = 0.2,0.4,0.8 kg工况下的爆坑数据。爆坑尺寸通常与炸药质量、起爆深度和土体特性相关[24],根据本文试验结果,相同药量地面爆炸下,钙质砂爆坑总体直径和深度普遍小于石英砂,且钙质砂爆坑呈现出更加明显的分层现象。一般情况下,迅速膨胀的爆炸产物推动砂体运动形成周边隆起的抛掷型爆坑,而钙质砂由于颗粒易碎和形状不规则导致孔隙比较高,在爆炸近区的初始骨架结构被爆炸冲击压力破坏,形成深度较浅的塌陷区。

    图  7  爆炸后钙质砂和石英砂表面形成的爆坑
    Figure  7.  Craters on the surface of calcareous sand and silica sand after the explosions
    图  8  爆坑尺寸定义
    Figure  8.  Definition of crater sizes
    表  3  钙质砂和石英砂爆坑尺寸
    Table  3.  Crater sizes of calcareous sand and silica sand
    W/kg材料d/mmD/mmh/mmh/mm
    0.2钙质砂1645606818
    石英砂19559011190
    0.4钙质砂2366739129
    石英砂19567012590
    0.8钙质砂28870010530
    石英砂850150
    下载: 导出CSV 
    | 显示表格

    图9为本试验采集到的原始土压力时程曲线,钙质砂比例距离有效数据范围为0.33~1.46 m/kg1/3,石英砂比例距离范围为0.27~1.44 m/kg1/3。当深度达到0.65 m时,压力时程曲线出现第2个峰值,为下方管廊顶板反射所致。

    图  9  爆心正下方不同深度的压力时程曲线
    Figure  9.  Pressure time-history curves at different depths directly below the explosion center

    接近底面的测点反射波较为明显,计算冲量时会产生巨大误差,因此对每个工况取前两层测点的压力时程曲线进行积分,得到对应位置的冲量时程曲线。图10为试验所得钙质砂和石英砂各点的比例冲量稳定值及拟合曲线,可以看出,在对数坐标系下,比例冲量随比例距离呈线性衰减,拟合公式如下:

    图  10  冲量试验结果及其拟合曲线
    Figure  10.  Experimental results and fitting curves of the scaled impulse
    ICSW1/3=1.26(RW1/3)1.40
    (1)
    ISSW1/3=1.78(RW1/3)1.21
    (2)

    由于波阵面到时及压力峰值到时的特征与弹塑性波传播规律一致,根据弹塑性球面波的传播理论,可以将爆心距-波阵面到达时刻的斜率视为弹性波速,爆心距-压力到达峰值时刻的斜率视为塑性波速[19]图1112对钙质砂和石英砂爆炸工况的数据进行了线性拟合,波速拟合值汇总于表4,其中CS-3塑性波速明显异常,将其剔除。根据表中数据可知:钙质砂中弹性波速为236~300 m/s,塑性波速为188~294 m/s;石英砂中弹性波速为218~337 m/s,塑性波速为234~285 m/s。在相同爆炸工况下,弹性波速和塑性波速随深度变化不明显,但是伴随着多次爆炸试验,中远区砂土密实度会略微提高,实测弹性和塑性波速较前一工况略微提高。本文试验所用钙质砂含水率为3.5 %,高于石英砂含水率(2.3 %),因此相同药量下,钙质砂的弹性波速均略高于石英砂的。

    图  11  爆心距与波阵面到达时刻线性拟合曲线
    Figure  11.  Linear fitting curves of explosion center distance and arrival time of wave front
    图  12  爆心距与压力峰值到达时刻线性拟合曲线
    Figure  12.  Linear fitting curves of explosion center distance and arrival time of pressure peak
    表  4  波速拟合结果
    Table  4.  Fitting results of blast wave velocity
    试验弹性波速/(m∙s−1)拟合优度塑性波速/(m∙s−1)拟合优度
    CS-12360.9961880.998
    CS-22530.9872360.999
    CS-32940.982361*1.000
    CS-43000.9762941.000
    SS-12180.9952341.000
    SS-22470.9982341.000
    SS-32900.9992851.000
    SS-42790.9992851.000
    SS-53370.999
    注:*为异常值。
    下载: 导出CSV 
    | 显示表格

    图13为爆炸波升压时间与比例距离的关系。可以发现,在钙质砂中爆炸升压时间随比例距离的增加而增加,从0.17 ms 增加至0.81 ms;在石英砂中爆炸升压时间在0.07~0.09 ms之间,与比例距离无明显关系,且升压时间均小于钙质砂的,这是由于,随着传播距离的增加,钙质砂表现出较石英砂更强的衰减效应,使得爆炸波的应力上升过程变缓。

    图  13  爆炸波升压时间与比例距离关系
    Figure  13.  The relationship between rise time and scaled distance

    美国陆军技术手册《常规武器防护设计原理》(TM5-855-1)[3]中对地面常规爆炸给出了土中自由场峰值压力的半经验半理论计算公式:

    p0=48.77fρc(5.4RW1/3)n
    (3)

    式中:p0为峰值压力,f为耦合系数;ρc为土介质的波阻抗;R为目标点到爆心的距离;W为装药质量;n为衰减系数。

    其中,耦合系数f定义为:近地面爆炸与封闭爆炸在同一介质中所产生的地冲击大小的比值,可表示为:

    f=(σ,v,d,I,a)near(σ,v,d,I,a)contained
    (4)

    对式(1)两边同时取对数可得:

    lgp0=nlg(5.4RW1/3)+lg(48.77fρc)
    (5)

    由式(3)可知,在对数坐标系下,峰值压力和比例距离的直线斜率即为衰减系数n,截距与土介质的波阻抗相关。图14为对数坐标系下钙质砂和石英砂试验测点的峰值压力与比例距离的关系,线性拟合可得:

    图  14  峰值压力实测值和拟合曲线
    Figure  14.  Experimental results and fitting curves of peak pressure
    lg(pCS0)=2.859lg(RW1/3)0.386
    (6)
    lg(pSS0)=2.786lg(RW1/3)0.180
    (7)

    因此钙质砂和石英砂峰值压力衰减可分别按下式计算:

    pCS0=0.41(RW1/3)2.859
    (8)
    pSS0=0.66(RW1/3)2.786
    (9)

    由式(6)~(7)可知,试验所用的钙质砂衰减系数为2.86,石英砂为2.79,因此在钙质砂中的爆炸波衰减速度大于在石英砂中的。将试验得到的钙质砂和石英砂密度、平均波速和衰减系数代入式(3),试验为地面爆炸,耦合系数f均取0.14。具体参数如表5图15为钙质砂和石英砂试验值和公式计算值的线性拟合曲线。由图可知,计算值明显小于试验值,因此直接使用TM5公式(式(3))计算钙质砂和石英砂的峰值压力会存在较大误差。

    表  5  计算参数
    Table  5.  Parameters of calculation
    材料ρ/(kg∙m−3)c/(m∙s−1)nf
    钙质砂1 3202712.8590.14
    石英砂1 3902722.7860.14
    下载: 导出CSV 
    | 显示表格
    图  15  峰值压力试验值与计算值对比
    Figure  15.  Comparison of peak pressures measured value and calculated value

    (1)爆炸荷载作用下钙质砂颗粒发生了破碎,大粒径颗粒含量减少,试验后平均粒径d50 = 0.42 mm,小于初始试样平均粒径d50 = 0.48 mm。

    (2)地面爆炸作用下,钙质砂表面形成的爆坑在俯视下呈两阶同心圆,分别为小而深的抛掷型内坑及大而浅的塌陷型外坑。装药质量相同时,钙质砂爆坑直径和深度普遍小于石英砂。

    (3)钙质砂中爆炸波弹性波速范围为236~300 m/s,塑性波速为188~294 m/s;石英砂中弹性波速为218~337 m/s,塑性波速为234~285 m/s。弹性波速和塑性波速均随炸药质量的增加而提高,本文试验所用钙质砂含水率稍高于石英砂,钙质砂中的弹性波速稍大于石英砂中的。

    (4)爆炸波在钙质砂中升压时间随比例距离增加而增加,当爆炸比例距离从0.329 m/kg1/3增加至1.154 m/kg1/3时,升压时间从0.17 ms增加至0.81 ms;石英砂中升压时间基本保持不变,在0.07~0.09 ms之间,且较钙质砂中升压更快。

    (5)对于本文中使用的低含水率砂样,钙质砂衰减系数为2.859,峰值压力按照pCS0=0.41(R/W1/3)2.859衰减;石英砂衰减系数为2.786,峰值压力按照pSS0=0.66(R/W1/3)2.786衰减;TM5手册峰值压力公式计算值小于试验值。

  • 图  1  深径比γ=8/3时计算模型

    Figure  1.  Calculation model for a depth-diameter ratio γ as 8/3

    图  2  抛掷爆炸示意图

    Figure  2.  Schematic of throw explosion

    图  3  侵彻通道深度与抛掷指数关系曲线

    Figure  3.  Curves of channel depth vs. throw exponential

    图  4  不同深径比下对称轴上碎屑体积分布

    Figure  4.  Debris volume distribution on symmetrical axis with different γ

    图  5  钻地通道深径比对抛掷距离的影响

    Figure  5.  Variation of throw distance with γ

    图  6  不同载荷下碎屑体积分数云图

    Figure  6.  Debris volume fraction contour under different loading

    图  7  不同荷载峰值条件下碎屑体积分数分布曲线

    Figure  7.  Curves of volume fraction under different loading

    图  8  抛距x、抛高y随着荷载的变化

    Figure  8.  Variation of throw distance x and throw height y with loading

    图  9  不同碎屑直径下通道内碎屑体积分数云图

    Figure  9.  Debris volume fraction contours for different debris particle sizes

    图  10  不同碎屑粒径下碎屑体积分数分布曲线

    Figure  10.  Curves of volume fraction under different debris particle sizes

    图  11  抛距x、抛高y随着碎屑粒径的变化

    Figure  11.  Variation of throw distance x and throw height y with debris particle size

    表  1  材料参数

    Table  1.   Material parameters

    材料 ν G/Pa ρ/(g·cm-3)
    碎屑颗粒 0.2 5×107 2.6
    钻地通道 0.2 5×107 2.6
    弹体 0.3 7×1010 7.8
    下载: 导出CSV
  • [1] 王涛, 余文力, 王少龙, 等.国外钻地武器的现状与发展趋势[J].导弹与航天运载技术, 2005(5):51-56. doi: 10.3969/j.issn.1004-7182.2005.05.011

    Wang Tao, Yu Wenli, Wang Shaolong, et al. Present status and tendency of foreign earth-penetrting weapons[J]. Missiles and Space Vehicles, 2005(5):51-56. doi: 10.3969/j.issn.1004-7182.2005.05.011
    [2] 汪斌, 曹仁义, 谭多望.大质量高速动能弹侵彻钢筋混凝土的实验研究[J].爆炸与冲击, 2013, 33(1):98-102. doi: 10.3969/j.issn.1001-1455.2013.01.014

    Wang Bin, Cao renyi, Tan Duowang. Experimental study on penetration of reinforced concrete by a high-speed penetrator with large mass[J]. Explosion and Shock Waves, 2013, 33(1):98-102. doi: 10.3969/j.issn.1001-1455.2013.01.014
    [3] 袁竹林, 朱立平, 耿凡, 等.气固两相流动与数值模拟[M].南京:东南大学出版社, 2013.
    [4] 潘振海, 王昊, 王习东, 等.油砂干馏系统的DEM-CFD耦合模拟[J].天然气工业, 2008, 28(12):124-126. doi: 10.3787/j.issn.1000-0976.2008.12.037

    Pan Zhenhai, Wang Hao, Wang Xidong, et al. DEM-CFD coupled simulation of oil sands retorting in a rotary kiln[J]. Natural Gas Industry, 2008, 28(12):124-126. doi: 10.3787/j.issn.1000-0976.2008.12.037
    [5] 戚华彪, 周光正, 于福海.颗粒物质混合行为的离散单元法研究[J].化学进展, 2015, 27(1):113-124. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hxjz201501011

    Qi Huabiao, Zhou Guangzheng, Yu Fuhai. Researches on mixing of granular materials with discrete element method[J]. Progress in Chemistry, 2015, 27(1):113-124. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hxjz201501011
    [6] 心男.基于EDEM-FLUENT耦合的气吹式排种器工作过程仿真分析[D].长春: 吉林大学, 2013. http://cdmd.cnki.com.cn/Article/CDMD-10183-1013194817.htm
    [7] 王国强, 郝万军, 王继新.离散单元法及其在EDEM上的实践[M].西安:西北工业大学出版社, 2010.
    [8] 张奇.球形装药岩石抛掷爆破初速度的数值模拟[J].煤炭学报, 1993, 18(5):27-32. http://www.cnki.com.cn/Article/CJFDTOTAL-MTXB199305007.htm

    Zhang Qi. Numerical modelling of initial speed of throw blasting using spherical charge[J]. Journal of China Coal Society, 1993, 18(5):27-32. http://www.cnki.com.cn/Article/CJFDTOTAL-MTXB199305007.htm
  • 期刊类型引用(2)

    1. 陈加浩,李俊超,朱斌,卢强,汪玉冰,管龙华,赵凤奎. 干砂地基初始应力对爆炸波传播影响的模型试验. 浙江大学学报(工学版). 2024(03): 579-588+634 . 百度学术
    2. 孙善政,王炜,刘渊,卢浩,赖杰. 分层岩土介质中的爆炸地冲击横向传播规律. 火箭军工程大学学报. 2024(06): 12-22 . 百度学术

    其他类型引用(2)

  • 加载中
图(11) / 表(1)
计量
  • 文章访问数:  4503
  • HTML全文浏览量:  1347
  • PDF下载量:  532
  • 被引次数: 4
出版历程
  • 收稿日期:  2015-03-20
  • 修回日期:  2016-01-20
  • 刊出日期:  2016-09-25

目录

/

返回文章
返回