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

阶跃信号电探针诊断爆轰加载下金属微喷现象

文雪峰 王晓燕 王健 任国武 洪仁楷 胡杨

李英雷, 刘明涛, 陈艳, 张世文, 汤铁钢. 线起爆膨胀柱壳实验加载及诊断技术[J]. 爆炸与冲击, 2022, 42(12): 124101. doi: 10.11883/bzycj-2021-0484
引用本文: 文雪峰, 王晓燕, 王健, 任国武, 洪仁楷, 胡杨. 阶跃信号电探针诊断爆轰加载下金属微喷现象[J]. 爆炸与冲击, 2019, 39(7): 074102. doi: 10.11883/bzycj-2018-0104
LI Yinglei, LIU Mingtao, CHEN Yan, ZHANG Shiwen, TANG Tiegang. Technologies for loading and diagnosis of expanding cylinder experiments with linearly-initiated explosives[J]. Explosion And Shock Waves, 2022, 42(12): 124101. doi: 10.11883/bzycj-2021-0484
Citation: WEN Xuefeng, WANG Xiaoyan, WANG Jian, REN Guowu, HONG Renkai, HU Yang. Characterizing detonation-induced micro-jetting by using step signal electric probe[J]. Explosion And Shock Waves, 2019, 39(7): 074102. doi: 10.11883/bzycj-2018-0104

阶跃信号电探针诊断爆轰加载下金属微喷现象

doi: 10.11883/bzycj-2018-0104
基金项目: 国家自然科学基金(11502253);中国工程物理研究院科学技术发展基金(2015B0201002)
详细信息
    作者简介:

    文雪峰(1988- ),男,硕士,助理研究员,scu_wxf@163.com

    通讯作者:

    王晓燕(1976- ),女,硕士,副研究员,wangxy101s@163.cm

  • 中图分类号: O383

Characterizing detonation-induced micro-jetting by using step signal electric probe

  • 摘要: 针对目前缺乏用于狭小测试空间复杂结构金属样品的微喷现象诊断技术的现状,探索使用阶跃信号电探针诊断爆轰加载下金属样品微喷现象的方法。设计阶跃信号电探针测试技术,开展仿真验证电探针在微喷物质导通下K+RX模式的放电机理。在爆轰微喷实验中通过电探针信号观察到准连续状态微喷区及其存在的两种动作演化过程:阶梯放电曲线表征微喷物质随靠近后界面密度逐步增加,多次放电现象表征微喷物质从连续状态拉升变为离散状态。采用微射流模型描述准连续状态微喷区物质状态,通过电平曲线计算出被测物质等效电阻,再通过等效电阻计算微射流的等效尺寸,从而可描述准连续状态微喷区物质的密度。
  • 对爆炸加载壳体膨胀断裂的研究主要集中在膨胀断裂机理和破片尺寸及其速度分布2个方面。在破片尺寸及速度分布研究方面:Gurney[1]提出了破片速度预估模型;Mott[2]、Grady等[3]、Hopson等[4]、Zhou等[5]和郑宇轩等[6]基于一维理论提出了破片尺寸分布模型[2-6];对于二维和三维情况,理论预测的破片尺寸与实际差异较大[7-8]。在壳体膨胀断裂机理方面,研究主要涉及柱壳的环向拉伸和剪切2种断裂模式,给出了断裂判据以及直观的裂纹扩展演化描述[9-11]。其基本依据是壳体内的应力分布和材料断裂判据。从研究情况看,壳体膨胀断裂的萌生位置分别有内表面、外表面、内表面附近等不同认识,对应的扩展路径和方式也有不同认识[9-14]。这些认识差异的起源,既有壳体内应力分布演化复杂,也有诊断信息不足造成分析难以收敛的原因。

    对于壳体内的应力分布,早期研究一般采用二维轴对称模型的解析解来简化分析滑移爆轰加载柱壳,未考虑滑移爆轰带来的沿柱壳轴向的应力梯度[9-11];随着数值模拟技术的应用,研究者发现滑移爆轰加载柱壳内的空间应力梯度分布具有显著、不可忽略的影响[15-17]。因此,将二维轴对称简化模型用于滑移爆轰加载的膨胀柱壳断裂分析,存在较大的认识偏差。对于数值模拟技术,虽然可以提供裂纹萌生前的应力分布描述,但是受材料模型参数获取难度的限制,其应用相对比较困难。因此,膨胀壳体断裂研究目前主要还是依赖实验诊断和简化应力分析开展。

    近年来发展起来的爆炸丝起爆技术(简称线起爆技术)[18-19]利用大电流通过金属丝并使其等离子体化,然后沿丝的长度方向同步起爆柱壳内部装填的炸药,从而能够在膨胀柱壳中实现一维轴对称应力的加载。在此状态下,就可以使用一维轴对称应力分析模型对线起爆膨胀柱壳的断裂问题进行有效分析。

    在壳体膨胀断裂诊断方面,常用手段包括高速分幅照相[11, 16]、高速狭缝扫描照相[20]、X射线透射照相[21]、DPS (Doppler detection system)测速[16]。其中,高速分幅照相是常用诊断手段之一,可观测壳体表面发生的皱褶、冒烟(爆轰产物泄漏)以及外径尺寸等特征信息。以往研究一般采用冒烟特征作为壳体断裂诊断标准[11],也有一些研究者将皱褶特征作为壳体断裂的诊断标准[22]。按照胡八一等[11]的认识,从壳体表面皱褶发展到冒烟需要几微秒到十几微秒。高速狭缝扫描照相可观测空间固定位置上的壳体轮廓投影的不连续变化(如冒烟、宏观断裂等)和外径尺寸信息。X射线透射照相可获取壳体沿照相投影方向的密度分布,进而通过密度梯度变化识别壳体的裂纹宽度及分布、外径尺寸等特征信息。X射线透射照相的优点是不受撞击或者爆炸发光影响,缺点是图像的对比度和清晰度较差,定量分析精度低[23]。DPS测速是近年来发展起来的常用诊断手段,用于测量壳体表面的垂直运动速度。DPS测速的优点是可以长时间精确测量。受柱壳弧度影响,上述照相诊断手段一般只能准确监测柱壳圆周部分角度范围内的断裂状态,因此诊断结果可能存在一定程度(或严重)的偏差。

    对于均匀承载柱壳,断裂将使局部壳体承载失效,进而导致壳体的应力、应变和速度出现显著的非均匀分布。中低应变率下的传统断裂诊断方法是采用粘贴在裂纹附近的应变计监测应变扰动来判读断裂时刻[24-25]。由于爆轰加载强度较高,不满足应变计的使用要求,因此可以采用一定数量、沿空间分布的DPS探头来监测壳体表面速度分布的演化,实现柱壳圆周范围内的断裂诊断。

    采用线起爆技术对装填粉末PETN (pentaerythritol tetranitrate)炸药的金属柱壳做一维柱面加载,装置结构见图1。线起爆金属丝安装在柱壳中轴位置。炸药与柱壳之间填充延展性较好的尼龙来约束粉末装药,同时抑制尼龙层碎裂而破坏加载均匀性。装填炸药的密度控制在(1.0~1.1)×103 kg/m3范围内,直径为15 mm。柱壳内径为40 mm,外径为48 mm。装药、尼龙层和柱壳高度相同,均为160 mm。柱壳材料为304钢和45钢。每种柱壳材料各开展1次实验。

    图  1  线起爆实验装置
    Figure  1.  The experimental device with linear initiation

    304钢柱壳实验仅采用DPS探头测量柱壳外壁的径向速度,以监测柱壳膨胀和断裂状态。在柱壳80 mm高度处,沿环向0°、90°、180°、270°方向各布置1个DPS探头,监测径向速度的对称性;基于柱壳沿高度方向的中心对称原则,在沿环向90°方向、80~140 mm高度范围内,按照20 mm等高度间隔布置了4个DPS探头,监测径向速度的一致性。

    45钢柱壳实验除采用DPS探头测量柱壳外壁速度外,还增加了高速分幅照相检测内容,以直观观测柱壳膨胀和断裂的宏观形态。其中,照相光路占据沿环向0°方向,因此该方向上的DPS探头取消。另外,考虑柱壳端部边侧稀疏对径向速度的影响,取消了沿环向90°方向、140 mm高度处的DPS探头。其余DPS探头布局与304钢柱壳实验的相同。此外,在45钢柱壳外壁20~140 mm高度范围内,按照20 mm等高度间隔分别画了蓝色环形高度标线。

    在时序控制方面,以金属丝通电起爆时刻为零时刻,并通过系统设置保障DPS测速和高速分幅照相与金属丝起爆的时基相同。两相邻高速分幅照片的时间间隔约为2 μs。

    304钢柱壳的径向速度监测结果见图2。受柱壳端部边侧稀疏卸载影响,140 mm高度处的速度明显低于柱壳中部其他高度处的速度。45钢柱壳的速度监测结果和高速分幅照相结果分别见图3图4图4(b)中5条裂纹带的最早出现时刻依次为39.2 μs (45°)、43.3 μs (0°)、47.4 μs (270°)、49.4 μs (315°)和51.5 μs (90°)。各相邻裂纹带的间距估算约为2 cm。受柱壳弧度影响,其中270°和90°裂纹带的最早出现时刻判读明显滞后。

    图  2  304钢柱壳外壁的径向速度曲线
    Figure  2.  The radial velocity curves of the outer surface of the 304 steel cylinder
    图  3  45钢柱壳外壁的径向速度曲线
    Figure  3.  The radial velocity curves of the outer surface of the 45 steel cylinder
    图  4  45钢柱壳的高速分幅照相结果
    Figure  4.  The high-speed framing photography results of the 45 steel cylinder

    图4(a)的有效段内,柱壳轮廓基本为直圆柱状。而在有效段外,圆柱轮廓有轻微的直径收缩。这一形态与图2中速度曲线积分获得的膨胀柱壳轮廓一致,即如图5中黑色曲线所示的中部平直凸出、两侧略滞后的鼓形轮廓。图6显示了图2中各速度曲线的对应位移状态。表1列出了柱壳有效段内各测点处的速度起跳时刻, t11t13分别为沿环向90°方向,在80、100和120 mm高度处柱壳速度的起跳时刻; t21t24分别为沿环向0°、90°、180°和270°方向,在80 mm高度处柱壳各测点的速度起跳时刻。

    图  5  304钢柱壳的初始轮廓与运动35 μs后的轮廓
    Figure  5.  The contours of the outer-surface of the 304 steel cylinder at the initial state and the deformed state after 35 μs
    图  6  304钢柱壳沿环向90°方向、不同高度的外壁位移状态
    Figure  6.  The displacement of the outer surface of the 45 steel cylinder at various heights and the circle angle of 90°
    表  1  304钢和45钢柱壳外壁各测点的速度曲线起跳时刻
    Table  1.  Jump-up times in velocity curves of the outside surfaces of the 304 steel and 45 steel cylinders
    柱壳材料t11/μst12/μst13/μst21/μst22/μst23/μst24/μs
    304钢10.2710.1710.0910.2710.2710.1410.25
    45钢12.1311.9711.8712.1312.1112.06
    下载: 导出CSV 
    | 显示表格

    加载状态分析用于确认加载满足一维柱面状态,从而为承载结构的应力分析提供简单的状态基础。从图2(a)和3(a)来看,80~120 mm高度范围的各速度曲线在柱壳断裂前基本保持重合,说明加载响应处于良好的柱面状态。而从图2(b)和图3(b)来看,沿环向各角度的速度曲线也基本保持重合,说明加载响应的轴对称性良好。因此,可以定性判断2种钢柱壳均处于良好的一维柱面状态。

    在定量判断方面,考虑加载状态偏离一维柱面将造成波形前沿出现一定程度的空间离散。定义空间离散度为单位长度波形前沿的最大距离差。为了便于实验监测,假定加载波在柱壳内稳定传播,将波形前沿的最大距离差转化为最大时刻差,得到空间离散度:

    ψ=ΔsΔL=cΔtΔL (1)

    式中:ψ为空间离散度,Δs为波形前沿最大距离差,ΔL为Δs的统计长度,Δt为波形前沿最大时刻差(或速度起跳时刻晃动),c为应力波在柱壳中沿径向传播的速度。

    按照一般小量考虑,可接受的空间离散度应不超过5%。因此,由式(1)可得速度起跳时刻晃动的要求如下:

    ΔtΔLc×5% (2)

    考虑柱壳在自由膨胀阶段的静水压相对较低,加载波的传播速度可简单近似为材料声速(计为5 km/s)。对于外壁半径为24 mm、有效段长80 mm的柱壳段,分别按照沿外壁环绕一周、沿轴向80 mm的长度统计柱壳有效段的波形前沿最大距离差,则可以通过式(2)计算出满足一维柱面状态的径向速度起跳时刻晃动上限:沿环向的时刻晃动Δt≤0.24 μs,沿轴向的时刻晃动Δt≤0.80 μs。从表1的数据来看,2种钢柱壳沿环向不同角度的速度起跳时刻晃动(0.13、0.07 μs)和沿轴向不同高度处的速度起跳时刻晃动(0.18、0.26 μs)均小于上述2个限值,表明柱壳近似处于一维柱面加载状态。

    图2(b)、3(b)显示,柱壳外壁的速度曲线呈现多峰值加载特征。根据应力波分析,在各峰值前的加速阶段,尼龙层对柱壳做碰撞加载;在峰值后的减速阶段,柱壳与尼龙层脱离、自由膨胀并在环向拉伸应力作用下近似线性减速。减速阶段的柱壳自由膨胀特征表现为速度峰值前后的起跳时刻间隔大于应力波在尼龙层内/外壁来回反射一次的时间(约11 μs)以及各减速阶段相近的加速度值(−(4.0~5.0) ×106 m/s2)。

    在图3(b)的第2个峰值后的减速阶段,45钢柱壳在初始约5 μs时段内保持均匀速度分布,并且环向承载水平与第1个减速阶段的近似相等(a′a1)。从37 μs时刻开始,沿环向90°方向的速度曲线呈现出较长时间的迅速下降趋势。假设此时段的柱壳仍为完整承载体,则该方向上的环向承载水平明显高于其余2个方向的(a″≈2a′),而其余2个方向的环向承载水平未发生变化。考虑到柱壳处于自由膨胀状态,其环向承载水平仅由柱壳承载能力控制。因此,由a″≈2a′关系可进一步推断该方向的柱壳承载力出现大幅度的阶跃强化。该推断结果与柱壳材料流动应力连续演化和柱壳质量沿环向均匀分布或连续演化(对应颈缩状态)的认识相悖,因此说明假设错误。柱壳在沿环向90°方向附近应发生了初始断裂并丧失环向承载能力,导致该方向的柱壳径向速度不受环向拉伸应力控制,而其余2个方向的柱壳则保持完整、环向承载水平未发生变化。

    图2(b)的第3个峰值时刻开始,沿环向0°方向的304钢速度曲线呈现出较长时间的水平演化趋势(加速度值为零),而其他3个方向的速度曲线则呈现自由膨胀的减速趋势(a′a1a2)。同样假设304钢柱壳为完整承载体,则沿环向0°方向的环向承载水平衰减至零,而其余3个方向的环向承载水平与之前各减速阶段的近似相等。由此推断,柱壳在沿环向0°方向附近发生断裂、丧失环向承载能力,而其余各方向上的柱壳则保持完整;如假设不成立,同样可以得出上述结论,因此说明该结论是唯一、真实的。

    为了直观展示分布式表面速度诊断方法的效果,将高速分幅照相诊断的45钢柱壳断裂情况与速度诊断结果进行了对比。在空间关系上,速度诊断的初始断裂位置与图4(b)中沿环向90°方向的裂纹带位置重叠,说明速度诊断的断裂位置是正确的;在时间关系上,速度诊断的初始断裂时刻(37.0 μs)明显早于同方向上照相诊断的裂纹最早出现时刻(约51.5 μs),并稍早于照相诊断的初始断裂时刻(39.2 μs),说明速度诊断方法可消除柱壳弧度影响,速度诊断的初始断裂时刻为柱壳的初始断裂时刻。

    表2是基于图2(b)3(b)中的速度曲线诊断得到的304钢和45钢柱壳的初始断裂参数。其中,断裂时刻根据速度曲线簇的分叉(或演化趋势变化)时刻判读;断裂应变通过断裂时刻前的速度曲线积分获得的位移量计算;平均应变率为柱壳外表速度起跳至断裂时刻的应变率平均值,用来近似表征柱壳开始承载变形至断裂时刻的应变率水平。从表2的断裂参数看,45钢和304钢柱壳的一维柱面动态拉伸断裂性能存在明显差异。45钢柱壳的断裂应变相对较小,其延展性弱于304钢柱壳。

    表  2  钢柱壳的初始断裂参数
    Table  2.  Initial fracture parameters of steel cylinders
    材料断裂时刻/μs断裂应变/%平均应变率/s−1
    304钢51370.9×104
    45钢37241.0×104
    下载: 导出CSV 
    | 显示表格

    采用电爆炸丝起爆PETN粉末炸药并驱动尼龙对金属柱壳加载的方式,结合分布式表面速度监测和高速分幅照相监测,开展了2种钢柱壳的膨胀断裂实验。通过实验研究,实现了金属柱壳的一维柱面均匀加载,建立了可监测柱壳圆周范围内初始断裂信息的分布式表面速度诊断方法,并获得了45钢和304钢柱壳的初始断裂参数(含断裂应变、平均应变率),具体结论如下。

    (1)基于线起爆粉末炸药并驱动尼龙加载的方式,可在金属柱壳中部、半柱高范围内实现一维柱面均匀加载。柱壳外表速度起跳时刻的晃动满足加载状态判据要求。

    (2)基于均匀承载壳体断裂引起的局部承载失效将导致均匀分布的速度曲线簇出现分叉(或演化趋势变化)的原理,分布式表面速度诊断方法可准确获取柱壳圆周范围内的初始断裂信息。标记沿柱壳环向分布的径向速度曲线簇出现分叉(或演化趋势变化)的时刻为断裂时刻,标记与速度曲线簇分叉的曲线对应区域为断裂位置。

    (3)在相同的一维柱面动态膨胀加载条件下,45钢柱壳的断裂应变(或延展性)低于304钢柱壳的。

    感谢金山、但加坤、郭昭亮、陈浩玉、罗振雄、张振涛、莫俊杰、李军、赵延安、谢明强、王远在实验设计、实施及结果分析方面给予的帮助。

  • 图  1  实验原理示意图

    Figure  1.  Schematic of principle of experiment

    图  2  阶跃信号电探针原理

    Figure  2.  Principle of an electric probe for detecting step signals

    图  3  电探针信号仿真

    Figure  3.  Simulated electric probe signal

    图  4  电探针信号实验结果

    Figure  4.  Experimental signals by electric probes

    图  5  物质状态界面提取

    Figure  5.  Extracting the mass state boundary

    图  6  Sn样品体密度

    Figure  6.  Volumetric density of Sn sample

    图  7  微喷物质密度计算

    Figure  7.  Calculation of micro-jetting mass density

  • [1] DURAND O, SOULARD L. Power law and exponential jetting size distributions from the dynamic fragmentation of shock-loaded Cu and Cn metals under melt conditions [J]. Journal of Applied Physics, 2013, 114: 194902. doi: 10.1063/1.4832758
    [2] MONFAREED S K, ORO D M, GROVER M, et al. Experimental observation on the links between surface perturbation parameter and shock-induced mass ejection [J]. Journal of Applied Physics, 2016, 116: 0636504.
    [3] ASAY J R. Material ejection from shock-loaded free surfaces of aluminums and lead: SAND-76-0542 [R]. 1976.
    [4] VOGAN W S, ANDERSON W W, GROVER M, et al. Piezoelectric characterization of jetting from shocked tin surfaces [J]. Journal of Applied Physics, 2005, 98: 113508. doi: 10.1063/1.2132521
    [5] SIGNOR L, DRAGON A, et al. Investigation of fragments size resulting from dynamic fragmentation in melted state of laser shock-loaded tin [J]. International Journal of Impact Engineering, 2010, 37: 887–900. doi: 10.1016/j.ijimpeng.2010.03.001
    [6] BUTTLER W T, ORO D M, OLSON R T, et al. Second shock jetting measurement with an explosive driven two-shockwave drive [J]. Journal of Applied Physics, 2014, 116: 103519. doi: 10.1063/1.4895053
    [7] MCMILLAN C, WHISPKY R. Holographic measurement of jetting from shocked metal surface [J]. SPIE High Speed Photography and Photonics, 1988, 1032: 553–557.
    [8] 汪伟, 李作友, 李欣竹, 等. 用超高速阴影摄影技术研究微喷现象 [J]. 应用光学, 2008, 29(4): 526–529. doi: 10.3969/j.issn.1002-2082.2008.04.010

    WANG Wei, LI Zuoyou, LI Xinzhu, et al. Study on micro-jet on ultra-high speed shadow photography [J]. Journal of Applied Optics, 2008, 29(4): 526–529. doi: 10.3969/j.issn.1002-2082.2008.04.010
    [9] ZELLER M B, MCNEIL W V, GRAY G T, et al. Surface preparation methods to enhance dynamic surface property measurements of shocked metal surfaces [J]. Journal of Applied Physics, 2008, 103: 083521. doi: 10.1063/1.2906107
    [10] 黄正平. 爆炸与冲击电测技术 [M]. 北京: 国防工业出版社. 2006: 617−623.
    [11] 王翔, 贾路峰, 傅秋卫, 等. 宽脉冲网络信号源及应用 [J]. 高压物理学报, 2005, 19(3): 279–283. doi: 10.3969/j.issn.1000-5773.2005.03.015

    WANG Xiang, JIA Lufeng, FU Qiuwei, et al. Broad pulse forming circuit and its appl ication [J]. Chinese Journal of High Pressure Physics, 2005, 19(3): 279–283. doi: 10.3969/j.issn.1000-5773.2005.03.015
    [12] 文雪峰, 王晓燕, 王健, 等. 可识别微层裂前界面的阶跃信号电探针测试技术 [J]. 爆炸与冲击, 2018, 38(2): 309–315. DOI: 10.11883/bcycj-2016-0271.

    WEN Xuefeng, WANG Xiaoyan, WANG Jian, et al. A step-signal electric probe technology for recognasing the front surface of micro-spall [J]. Explosion and Shock Waves, 2018, 38(2): 309–315. DOI: 10.11883/bcycj-2016-0271.
    [13] SEIFTER A, FURLANETTO M R, GROVER M, et al. Use of IR pyrometry to measure free-surface temperature of partially melted tin as a function of shock pressure [J]. Journal of Applied Physics, 2009, 105: 123526. doi: 10.1063/1.3153973
  • 期刊类型引用(1)

    1. 张伟,孙毅萌,杨大峰,许宁宁. 池火灾环境下喷涂涂料液化石油气储罐热响应机理试验研究. 河南科技大学学报(自然科学版). 2024(04): 49-56+6-7 . 百度学术

    其他类型引用(0)

  • 加载中
图(7)
计量
  • 文章访问数:  5093
  • HTML全文浏览量:  1913
  • PDF下载量:  35
  • 被引次数: 1
出版历程
  • 收稿日期:  2018-03-30
  • 修回日期:  2018-08-24
  • 网络出版日期:  2019-06-25
  • 刊出日期:  2019-07-01

目录

/

返回文章
返回