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

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

顾强 张世豪 安晓红 张亚

2024 年 1 期封面[J]. 爆炸与冲击, 2024, 44(1).
引用本文: 顾强, 张世豪, 安晓红, 张亚. 基于FLUENT-EDEM耦合的爆炸抛掷特性研究[J]. 爆炸与冲击, 2016, 36(5): 611-616. doi: 10.11883/1001-1455(2016)05-0611-06
Cover[J]. Explosion And Shock Waves, 2024, 44(1).
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软件联合的数值模拟方案,分析了不同深径比条件下钻地通道内介质碎屑的抛掷过程。进一步研究了介质碎屑运移特性的影响因素。结果表明:介质碎屑的运移能力与爆炸荷载峰值成正比,与钻地通道深径比及碎屑粒径成反比。
  • 金属及合金在强动载荷作用下将发生近表面破碎,产生大量高速运动碎片和粒子,对其他结构或器件造成损伤破坏。因此,金属及合金材料近表面的动态破碎是爆炸与冲击动力学领域高度关注的关键科学问题,对该问题的研究在装甲防护、航空航天、公共安全等方面具有重要应用价值[1-5]。Walsh等[6]和Asay等[7-9]早就开始了材料近表面的动态破碎研究。随着诊断技术的不断发展,研究人员通过实验研究发现,材料近表面动态破碎存在3种主要机制,即层裂、微层裂和微喷射,这3种机制强烈依赖于冲击压力、温度、相结构、材料微结构和加载波形等因素[10-12]。目前,相关研究基本都以简单一维应变加载为主,关于材料在复杂加载状态下动态破碎的研究还相对较少。因此,对复杂应变/应力条件下材料近表面动态破碎的认识严重不足。

    马赫透镜作为一种高压加载的特殊方法,通过冲击波的马赫反射和对碰可以在材料内部形成马赫杆加载,且在马赫透镜的中心位置产生三角波,它已经被用于高压状态方程[13-14]和疏松材料压实[15]方面的研究。该三角波与平板撞击、爆轰加载和激光冲击下产生的一维应变三角波存在显著不同。前者的压力幅值在一定的时间内随传播距离的变化几乎不变,而后者的压力幅值则会因为追赶稀疏波的影响随时间的变化急剧衰减。因此,马赫杆加载虽然是一维应变加载,但其应力状态更复杂,为研究复杂加载条件下材料的动态破碎提供了新的机遇。

    本文中,将数值模拟和实验相结合,利用大口径火炮在无氧铜样品中产生马赫杆加载,并利用激光粒子速度干涉仪测量样品自由面的粒子速度剖面,研究在马赫杆加载下无氧铜样品的近表面动态破碎行为,揭示在不同加载条件下无氧铜样品自由面的2种动态破碎机制。

    马赫透镜构型如图1所示,内外圆柱由不同波阻抗的金属组成。飞片高速撞击后,外层圆柱材料中的冲击波速度比内圆柱中的高,使外侧冲击波能透射入内层圆柱中,形成锥形收敛,会聚到中心轴线上,并与内圆柱中的冲击波叠加产生再加载,使内圆柱的局部压强和粒子速度显著提升。边界处的外侧冲击波入射可看作是斜冲击波的反射状态,通过斜冲击波的入射可以在内圆柱材料中心线上产生稳定的马赫杆[14]。虽然马赫杆加载宏观上表现为一维应变状态,但整个内圆柱中的冲击波阵面随时间和空间剧烈变化,因此其实质的加载状态更复杂。

    图  1  马赫透镜结构
    Figure  1.  Structure of Mach lens

    在内圆柱中,除了波阵面处的斜激波入射与反射,在界面处还有内外冲击波的相互作用,结果是向内侧圆柱中反射稀疏波。在入射波向中心汇聚时,反射的稀疏波也到达中心,跟随在马赫杆后。从而在马赫杆后形成低压卸载区,使冲击前沿压强到达峰值后又迅速下降,马赫杆波形就由方波转变成三角波。

    基于马赫透镜的原理分析,对马赫杆加载进行了模拟设计。利用有限元分析软件AUTODYN,采用Lagrange方法和二维旋转轴对称结构建立模型。材料动力学模型选用最常见的Mie-Grüneisen状态方程[16]和Johnson-Cook(JC)本构模型[17],其模型参数分别见表1[18-20]表2[17, 21-23]。模拟设计的两类马赫杆实验(Mach-1和Mach-2)参数如表3,分别对应在无氧铜中产生最高压力99.4和24.2 GPa。两类马赫杆实验的压力流场随时间的演化如图2所示。

    表  1  Mie-Grüneisen状态方程参数[18-20]
    Table  1.  Parameters of Mie-Grüneisen equations of state[18-20]
    材料ρ0/(g·cm−3)c0/(km·s−1)sγ来源
    304不锈钢7.904.571.49 1.93文献[18]
    LY12铝2.795.371.29 2.0 文献[19]
    无氧铜8.933.941.4892.02文献[18]
    TC4钛合金4.425.131.0281.23文献[20]
    下载: 导出CSV 
    | 显示表格
    表  2  Johnson-Cook本构模型参数[17, 21-23]
    Table  2.  Parameters of the Johnson-Cook constitutive model[17, 21-23]
    材料A/MPaB/MPanCmTm/K来源
    304不锈钢31010000.650.071.01673文献[21]
    LY12铝3696840.730.00831.7775文献[22]
    无氧铜902920.310.0251.091356文献[17]
    TC4钛合金8623310.340.0120.82110文献[23]
    下载: 导出CSV 
    | 显示表格
    表  3  设计参数
    Table  3.  Parameters of experimental design
    实验编号 材料 飞片速度/(km·s−1) 厚度/mm 直径/mm
    飞片 外圆 内圆 飞片 样靶 内圆 外圆
    Mach-1 304不锈钢 LY12 铝 无氧铜 1.40 3.0 16.0 4.8 38.0
    Mach-2 304不锈钢 TC4钛合金 无氧铜 0.50 12.0 26.0 14.0 45.0
    下载: 导出CSV 
    | 显示表格
    图  2  两类马赫杆加载下压力流场随时间演化的模拟结果
    Figure  2.  Numerically-simulated evolutions of pressure contours under two types of Mach stem loading

    图3为实验Mach-1中内圆柱中心线上等间距拉格朗日点的压力和粒子速度剖面模拟结果。0.7 μs前,冲击波为方波。0.7 μs后,峰值粒子速度开始提升,且波形开始变化。其变化来自两部分,分别是斜冲击的入射波和界面反射的稀疏波,导致在提升峰值粒子速度的同时,后续的卸载稀疏波也变得更陡峭,使冲击平台逐步缩短。2.2 μs时,粒子速度幅值基本保持不变,且能在较长时间内保持这种三角波的形式。

    图  3  实验Mach-1中圆柱试样靶中心线上等间距拉格朗日点的压强和速度模拟结果
    Figure  3.  Simulated velocity- and pressure-time evolutions of equidistant Lagrangian particles on the symmetrical centerline of the cylindrical specimen target in experiment Mach-1

    内外圆柱之间采用间隙配合,用环氧树脂粘结。为确保冲击实验时飞片与靶装置正面碰撞,将样品安装在有机玻璃支架上,经校准后用样品环将样品固定在炮口上,如图4(a)所示。火炮驱动304钢飞片以设计速度撞击马赫透镜装置。采用多普勒光探针系统(Doppler pin system, DPS)测量内外圆柱自由面的粒子速度,DPS探头分布见图4(b)~(c)。在实验Mach-1中,在内圆柱中心处设置一根DPS探针,用于测量马赫杆中心的运动历程。在铝套筒后表面距圆心3.0、5.0、7.0、9.0 mm处,布置系列DPS探针,用以判断碰撞的倾斜角度以及比较模拟结果与实验结果,并确定模拟需要的黏性等参数。在Mach-2实验中,碰撞速度较低,马赫杆的范围更大,因此在内圆柱的半径0、1.0、1.4、1.8、2.6和3.4 mm处布置了DPS探针,对应图4(c)中的点EFGHIJ,用于分析马赫杆的径向行为。飞片的碰撞速度也通过DPS测速确定。

    图  4  两类马赫杆加载实验设计
    Figure  4.  Experimental arrangements of two types of Mach stem loading

    DPS测量获得的飞片碰撞速度分别为1325和605 m/s。在实验Mach-1中,铝套筒后表面不同位置处的粒子速度剖面如图5所示。由同心圆上4个DPS信号的起跳时间间隔可知,碰撞倾斜角度为0.376°,因此该实验中飞片与靶装置的碰撞可以近似为完全正碰。在实测弹速下,铝套筒后自由面粒子速度的数值模拟结果和实验结果具有较好的一致性,由此证实实验中实现了预期的马赫杆加载。马赫杆稳定加载后的峰值压力分别达到95.75和32.38 GPa。结合模拟进行进一步分析发现,在自由面半径为3.0 mm的位置处,由于靠近内外圆柱的界面,受到内边界反射稀疏波的影响,冲击上升沿的峰值强度有所降低;在自由面半径为5.0 mm的位置处,冲击波基本没有受到来自外侧和内测边界处卸载稀疏波的影响,是一个完整的矩形冲击波;而在自由面半径为7.0、9.0 mm的位置处,冲击波则是受到了外侧稀疏波的影响,粒子速度剖面呈现三角波的特征,冲击波在到达自由面后开始卸载,粒子速度明显下降。图5(b)~(d)的粒子速度剖面在1.0 μs左右出现回弹上升,应该是发生了层裂,而且越接近外侧,层裂出现越早。

    图  5  实验Mach-1中铝套筒后自由面粒子速度剖面实验结果与模拟结果的比较
    Figure  5.  Comparison of rear free-surface particle velocity profiles of the aluminum sleeve in experiment Mach-1with the corresponding simulated ones

    实验Mach-1中内圆柱中心处的DPS测量结果见图6,可以看出无氧铜后自由面可能发生了微层裂,导致近表面物质碎化为细小碎片云并以接近自由面的速度脱离基体。这是由于稳定后的马赫杆中心处的压力由碰撞初始的28.92 GPa提升到了95.75 GPa,使无氧铜在卸载后进入近熔化区。

    图  6  实验 Mach-1中内圆柱中心处的DPS时谱图和自由面粒子速度剖面
    Figure  6.  DPS time-frequency spectrum and free-surface particle velocity profile at the center of the inner cylinder in experiment Mach-1

    在实验Mach-2中,由于马赫杆的范围更大,因此可以研究马赫杆不同位置的动力学过程。模拟设计中,点EFGH位于马赫杆平台范围内,而点IJ位于过渡区。如图7所示,粒子速度剖面都有明显的三角波层裂信号。为描述该层裂行为,数值模拟中采用最大拉应力断裂模型(设最大拉应力为2.6 GPa),可以获得实验Mach-2中不同位置处的粒子速度剖面。模拟结果与实验结果进行比较,可见两者的一致性较好,数值模拟较好复现了不同位置的层裂信号。该层裂行为是由于马赫杆产生的三角波加载到达自由面后经卸载产生拉伸应力所致。如图8所示,马赫透镜中内圆柱冲击波形成的层裂除了轴向应力的影响外还与其在径向分布的差异有关,不同位置的层裂规模和距离差异较大。内圆柱中的马赫杆冲击波阵面沿径向分布不均,冲击波阵面为锥状,点E~J的峰值速度和冲击波到达时间分别随半径增大而减小和提前(见图7)。这与平面加载的层裂行为存在显著的不同,平面加载的层裂在其径向上并无变化。中心线处的马赫杆虽然到达自由面的时间最迟,但强度最高,形成的三角波层裂位置也最靠前,层裂片速度最快,而两侧的冲击波强度相对较低,层裂位置也靠后。最终的层裂片沿径向各个位置的速度和厚度也不同,整体呈凸形(见图8)。

    图  7  实验Mach-2中不同位置处粒子速度剖面的实验结果与模拟结果的比较
    Figure  7.  Comparison of experimental and simulated results of particle velocity profiles at different positions in experiment Mach-2
    图  8  实验Mach-2中试样破坏形貌的模拟结果
    Figure  8.  Simulated failure morphologies of the specimen in experiment Mach-2

    以无氧铜为研究对象,模拟设计了2类不同冲击波压力的马赫杆加载实验。模拟结果与实验结果十分接近,马赫杆相关的主要特征也都比较符合,表明在实验上用马赫透镜实现了马赫杆加载。同时,观察到了无氧铜在2种不同马赫杆加载下产生的2种不同近表面破碎模式。在高压下,由于无氧铜卸载进入近熔化区,发生微层裂破坏。在低压下,无氧铜处于固态,三角波发生自由卸载形成层裂。与平面冲击产生的层裂不同,由于马赫杆冲击波状态径向分布不均匀,产生的层裂片呈凸形。

    该工作对进一步了解金属材料的动态破碎行为具有一定价值,并为研究复杂加载状态下的材料破坏提供了新的实验方法。

  • 图  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
  • 加载中
图(11) / 表(1)
计量
  • 文章访问数:  4480
  • HTML全文浏览量:  1343
  • PDF下载量:  524
  • 被引次数: 0
出版历程
  • 收稿日期:  2015-03-20
  • 修回日期:  2016-01-20
  • 刊出日期:  2016-09-25

目录

/

返回文章
返回