Chirped fiber Bragg grating sensor for side detonation velocity measurement of the explosion
-
摘要: 提出了一种适用于啁啾光纤布拉格光栅CFBG(chirped fiber Bragg grating)传感器测量爆轰波速度的实验公式，该公式考虑了ASE(amplified spontaneous emission)光源的光谱、CFBG的平坦度对测量的影响；研究了CFBG的反射谱线宽与长度的理论关系，并设计了测量CFBG反射谱线宽的实验装置，从而可计算出CFBG长度。搭建了测量CFBG传感器系统常数的实验装置，并利用实验数据拟合了系统常数；用CFBG传感器测量了RDX/TNT(60/40)的旁侧爆轰波速度，与其爆速有3.4%的差异。
-
关键词:
- 爆炸力学 /
- 啁啾光纤布拉格光栅传感器 /
- 反射谱 /
- 爆轰波速度 /
- ASE光源
Abstract: An experimental equation of CFBG sensor for detonation velocity measurement which was based on the flatness of ASE light source and CFBG reflection spectrum, was quantitatively established for the first time. The CFBG length could be calculated by the theoretical relation between it and the CFBG reflection spectrum line-width, which was measured experimentally. The constant of the CFBG sensor was fitted by the experimental data. The side detonation-wave velocity of RDX/TNT (60/40) was measured with 3.4% error compared to the publication value. -
在常规武器系统中，聚能装药战斗部扮演着越来越重要的角色，提高炸药能量利用率是战斗部设计者追求的目标之一。带隔板装药不仅结构简单而且可以在不增加装药量以及不改变装药长径比的条件下，提高其对药型罩的压垮做功能力。典型的带隔板装药结构如图 1所示，根据惠更斯原理[1]，爆轰波经过隔板的绕射后将以O′点开始传播，以较小的入射角驱动药型罩形成聚能侵彻体。这类带隔板装药在驱动药型罩形成高速杆式射流、大长径比爆炸成型弹丸以及利用起爆控制技术形成多模毁伤元等方面应用非常广泛。
图 1 带隔板柱形装药结构及爆轰波示意图Figure 1. Schematic diagram for cylindrical charge with waveshaper and detonation wave早在20世纪70年代，汤明钧等[2]利用高速摄影测定了带隔板破甲弹中装药端面上爆轰波的爆速和爆压，得出马赫波的波速达到正常爆速的1.9倍，爆压达到正常CJ爆压的3倍。杨全中等[3]对带隔板装药爆轰波的相互作用进行了研究，并利用Whitham方法[4]计算了爆轰波传播过程中的波形和碰撞临界角等参数。B.B.Dunne等[5]指出了发生马赫反射的临界入射角取决于爆轰产物的多方指数，并给出了马赫杆压力等参数的计算方法。L.M.Hull等[6]计算了球面爆轰波的马赫反射参数，通过实验给出了PBX9501和PBXN-9等炸药的马赫反射参数。张俊秀等[7]进行了爆轰波马赫反射的实验研究，并根据修改的Whitham方法[4]计算了马赫反射的压力、爆速和增长角等参数。Zhang Xianfeng等[8]、张先锋等[9]分别利用理论计算和数值模拟分析了夹层装药爆轰波的传播过程，并基于质量守恒给出了爆轰波传播过程中马赫杆高度的计算模型。
本文中，基于三波理论[3]和Whitham方法[4]，对带隔板装药爆轰波相互作用的正规反射和马赫反射进行理论分析，并给出马赫杆高度的计算模型。在此基础上，利用AUTODYN有限元计算程序中的Lee-Tarver模型[10]对带隔板装药爆轰波传播过程进行数值模拟, 以期研究成果可为描述带隔板装药驱动典型介质提供理论基础和数据支持。
1. 带隔板装药爆轰波马赫反射理论模型
1.1 带隔板装药爆轰波相互作用
2个等强度爆轰波的正规斜碰撞，可视为其中一侧爆轰波对刚性壁面的正规斜碰撞，当入射角φⅠ较小时，流入界面I和界面R的粒子流的偏转角相等。入射角不断增大，达到马赫反射的临界入射角φⅠ, c时，偏转角θ2小于θ1。粒子流的不断积累迫使碰撞点P移动到Q点，便形成了图 2所示的马赫反射。基于F.Grasso等[11]关于爆轰波反射的描述，爆轰前反应区用带下标1的参数表示，反射波前反应区用带下标2的参数表示，超压爆轰反应区用带下标3的参数表示。
假设炸药完全反应且满足CJ条件，区域(1)中的爆轰参数可由炸药的CJ参数、初始密度ρ0、多方指数γ和入射角获得。粒子流通过区域(1)发生偏转流入区域(2)，偏转角和反射角满足下列关系:
tanθ1=[Ma21sin2(φⅡ+ε)−1]tan−1(φⅡ+ε)Ma21[(γ+1)/2−sin2(φⅡ+ε)]+1 (1) 式中Ma1=√1+(γ+1)2/(γtanφⅠ)2为马赫数; 发生正规反射时，ε=θ1；发生马赫反射时，ε=θ1-θ2。
发生正规反射和非正规反射时, 式(1)中ε的表达式是不同的，若用发生正规反射时ε=θ1求解发生非正规反射时式子(1)的表达式，在0＜φⅡ+ε＜π/2范围内必定无解。根据此条件可以获得发生马赫反射的临界入射角，不同炸药临界入射角的理论计算值和实验结果如表 1所示, 表中DCJ为炸药CJ爆轰速度。通过比较发现，两者比较接近，从而可以验证马赫反射临界入射角计算的正确性。通过计算还发现，发生马赫反射时的临界入射角与炸药多方指数成一一对应关系，随着多方指数的增大，发生马赫反射的临界入射角不断减小，如图 3所示，此结论与B.B.Dunne[5]给出的结果一致。
表 1 不同炸药装药中发生马赫反射时临界入射角的计算值与实验结果Table 1. Calculated and experimental results of critical incident angles corresponding to Mach reflection in different explosive charges
图 3 马赫反射临界入射角随多方指数的变化Figure 3. Critical incident angle for the onset of Mach reflection varying with polytropic exponent大量理论和实验研究发现[5, 11-12]，马赫杆是一个弯曲的曲面，交于对称轴P点，相切于爆轰波(I)Q点。在P点, 偏转角α=θ1-θ2=0°; 在Q点, α=θ1-θ2；马赫杆与对称轴切线的夹角从π/2变化到φⅠ。由于碰撞点P处的压力和密度等都远大于CJ值，难以准确获取，本文中借助炸药爆轰产物的状态方程p=p(E, v)来求得此处的压力, 其中p为压力，E为炸药能量，v为比容。由H.M.Sternberg等[12]的定义，产生马赫反射时爆轰产物中将从三波点产生一条滑动线。此直线通过的地方，压力是连续的，但速度和密度非连续，因此，在p3=p2处存在一个三波点。
1.2 马赫杆高度模型
基于Whitham方法[4, 13]获取马赫反射的增长角和高度等参数，并将马赫波简化为平面波。将爆轰波的每一元素沿每个微元射线管的传播当作固壁管道中的激波传播来处理，在二维情况下，爆轰波的位置和与其垂直的射线位置如图 4所示。爆轰波的传播过程可以看成是在变截面流管中的传播过程，当截面发生微小变化时，流动符合等熵流性质，得到马赫反射的A(z)关系式为：
图 4 马赫反射后爆轰波形位置与其垂直射线的关系简图Figure 4. Relationship of detonation wave and its vertical line after Mach reflectionAf(z)=A0f(1)=const (2) 式中:z=p/pCJ为爆轰波强度，且：
f(z)=z1/γ(2z−1)1/2B(z)C(z);C(z)=exp{γ−1/2sin−1[(3−γ)z−2(γ+1)z]} B(z)={[2(γ−1)(2z−1)]1/2+2[(γ−1)z+1]1/2}[2γ/(γ−1)]1/2 基于上述讨论，如图 4所示, 假设Δt时刻后，由爆轰波NS变成马赫波PQ，φⅠ为爆轰波入射角，χ为三波点增长角，N和Q点为不同时刻的三波点。设NS管道截面积为ACJ，PQ处管道截面积为AM，根据几何关系知:
tanφⅠ=(DCJDM+AMACJ){[1−(DCJDM)2][1−(AMACJ)2]}−1/2 (3) cosχsin(φⅠ−χ)=DMDCJ (4) sinχcos(φⅠ−χ)=AMACJ (5) 由式(5)知:
AMACJ=f(1)f(z) (6) 参考B.D.Lambourn等[14]的研究，将爆速DM和DCJ定义为:
DMDCJ=z(2z−1)1/2 (7) 上述计算中未考虑爆轰波后稀疏波的影响，计算结果误差较大，对Whitham方法[4, 13]中的A(z)关系式进行修正:
ACJAM=(2z−1)12k(z) (8) 式中:k(z)是慢变函数，一般取常数, 通常由实验确定。
根据张俊秀等[7]提供的入射角和爆速参数计算方法，通过方程(2)~(8)，得到马赫波参数DM、z和χ的计算结果如表 2所示。利用Whitham方法[4, 13]计算获得的马赫波速度和三波点增长角与实验结果[7]吻合较好，从而验证此方法是可行的。
φⅠ/(°) DCJ/(km·s-1) DM, exp/(km·s-1) DM, cal/(km·s-1) DM,exp−DM,calDM,exp/% z χexp/(°) χcal/(°) χexp−χcalχexp/% 46.5 4.14 5.89 5.86 0.5 3.42 1.5 1.55 -3.2 48.5 4.20 5.77 5.79 -0.3 3.22 2.27 2.13 6.6 49.5 4.48 6.07 6.11 -0.7 3.13 2.80 2.45 14.3 同一装药中炸药的CJ爆速为常数，三波点增长角只与入射角和k(z)有关。根据带隔板装药爆轰波传播的特点，将入射角离散化，网格点(n, n)上的值可由网格点(n-1, n-1)上的值求得，马赫杆的高度可以定义为:
{Hn=Rtanχn−1tanφn−1−Ln−1tanφn−1+Hn−11+tanχn−1tanφn−1Ln=Ln−1tanχn−1tanφn−1+(R−Hn−1)tanφn−11+tanχn−1tanφn−1 (9) 式中:R为隔板半径；Hn和Ln分别为马赫杆的高度和传播距离，n=1, 2, 3, …。
当n=1时，H0=1, L1=Rtanφ0。L.M.Hull[5]也认为，当入射角为临界入射角时，马赫杆高度为零或非常小。结合表 1的数据和方程(3)~(9)，可以确定马赫杆高度。
爆轰波绕射过隔板后，爆轰波不断向轴线方向汇聚，入射角不断增大。当入射角小于临界入射角时，发生正规反射，马赫杆高度为零。当入射角大于临界入射角时，发生马赫反射，马赫杆高度随着入射角的增大而增大，相同炸药中，随着k(z)的变化，马赫杆高度计算结果如图 5所示。马赫杆在不断增大的过程中，马赫杆上的压力和速度迅速降低，当马赫杆上的压力和速度降低到CJ值时，马赫杆消失，此时，马赫反射变成正常的CJ爆轰。
图 5 马赫杆高度随入射角变化的计算结果Figure 5. Calculated variation of Mach stem height with incident angle2. 带隔板装药爆轰波传播及马赫反射过程数值模拟
2.1 有限元计算方案
根据前述对带隔板装药爆轰波传播过程的分析，在保证隔板外侧爆轰波可以正常传播的条件下，隔板直径越大，调整爆轰波形状的能力越强，在装药结构中爆轰波马赫杆的形状就越明显。带隔板装药爆轰波的传播过程是一种多物质相互作用的大变形运动过程，因此本文中采用多物质耦合的Euler算法描述带隔板装药爆轰波的传播过程。图 6所示为装药的有限元模型。
图 6 带隔板装药数值模拟的有限元模型Figure 6. A finite element model for numerical simulation on charge with waveshaper2.2 材料模型及参数
凝聚炸药爆轰产物状态方程描述炸药爆轰达到CJ状态后爆轰产物中压力、温度和体积等物理量之间的关系，不仅体现炸药本身的做功能力，而且是研究爆炸力学各种问题的基础。本文中分别采用点火与增长模型的Lee-Tarver状态方程[10]和JWL状态方程[15]来描述爆轰波绕射过隔板后在炸药中的传播过程。炸药参数主要来自文献[10, 15]及AUTODYN中内置的参数。
为了获取环形爆轰波形并保证透射波不起爆主装药，隔板厚度应足够大，这里隔板厚度选20 mm。隔板直径选42 mm，药柱直径选50 mm，使得隔板外侧装药厚度大于爆轰波正常传播时装药的临界直径。隔板材料为环氧树脂，其状态方程为Puff状态方程[16]，主要参数分别为：ρ0=1.985 g/cm3，A1=17.89 MPa，A2=201.2 MPa，A3=-6.053 MPa，Grüneisen系数为1.72，膨胀系数为0.25。
2.3 数值模拟结果与分析
针对上述讨论，对主装药分别采用JWL[15]和Lee-Travel[10]2种不同的状态方程进行描述，爆轰波波形的模拟结果如图 7所示。基于JWL模型[15]所形成的爆轰波的波阵面存在多段线组合形式，如图 7(a)所示。产生该结果的主要原因是：由于JWL模型为程序起爆模式，而炸药微元的爆轰采用预先设定时间起爆，JWL模型在描述带隔板装药爆轰(特别是隔板直径较大)时存在局限性，较难精确描述炸药微元的超压爆轰行为。基于炸药冲击起爆点火与增长的Lee-Travel模型[10]形成的爆轰波形较光滑，如图 7(b)所示。通过对比发现，Lee-Travel模型[10]能更真实地描述带隔板装药爆轰波的传播过程。
图 7 采用不同的炸药状态方程模拟得到的爆轰波波形Figure 7. Detonation waveform of explosive simulated by different equations of state对带隔板装药中的主副装药分别采用JWL[15]和Lee-Travel[10]2种不同的状态方程进行描述，爆轰波传播过程中不同时刻的波形如图 8所示。装药起爆后，爆轰波开始以球面波向前传播，经过隔板绕射后变成汇聚爆轰波。汇聚爆轰波在装药轴线处发生相互碰撞反射，当t=86 μs时，爆轰波已经发生马赫反射，随着爆轰波的传播, 马赫杆的高度不断增加，如图 8(b)和(c)所示。随着爆轰波的传播，马赫杆的高度不断增加，但其压力和速度迅速降低，当马赫杆上的压力和爆轰速度等参数降低到CJ值时，马赫反射消失，马赫反射变成正常爆轰, 如图 8(d)所示。
为了进一步了解带隔板装药爆轰波传播过程中发生马赫发射后参数的变化，在装药轴线以及距轴线2.5、5.0、7.5和10.0 mm处分别设有一系列高斯点。以TNT炸药为例，通过设置高斯点可以获取爆轰产物的参数变化，当马赫反射发生后，随着爆轰波的传播，不同位置处压力和密度的变化见图 9~10。发生马赫反射时，装药轴线处(即马赫杆碰撞点处)的压力和密度均达到最大值，其中最大压力与CJ压力的比值达到4.3，最大密度与CJ密度的比为2.05，与Zhang Xianfeng等[8]对TNT装药发生马赫反射时最大压力和密度的理论计算值接近。随着爆轰波的传播，马赫杆碰撞点的压力和密度逐渐降低。通过对距轴线不同距离处设置高斯点发现，压力和密度先增大后减小，这说明马赫杆存在一定高度且是不断变化的。
图 9 马赫反射发生后，随着爆轰波的传播，在轴线不同位置处压力的变化Figure 9. Pressure at different positions of the axis after Mach reflection
图 10 马赫反射发生后，随着爆轰波的传播，在轴线不同位置处密度的变化Figure 10. Density at different positions of the axis after Mach reflection为了获取带隔板装药中马赫杆高度的变化，分别对带隔板装药中的Comp B、TNT和PBX9404炸药爆轰波的传播过程进行了数值模拟。随着爆轰波的传播，马赫杆高度的精确捕捉越来越难，其数值模拟结果和理论计算结果如图 11~13所示。
当入射角较小时，马赫杆的高度为零，随着入射角不断增大并达到临界入射角时，马赫反射发生。此时，马赫杆碰撞点的压力和密度均达到最大值，随爆轰波的传播，马赫杆高度不断增加，马赫杆上的压力和密度不断减小。当马赫杆上的压力和密度等降到CJ值时，马赫反射消失，变成正常爆轰。
3. 结论
通过对带隔板装药爆轰波马赫反射过程的理论分析以及数值模拟，得到的主要结论如下：
(1) 基于三波理论和修正的Whitham方法，得到了马赫反射临界入射角和三波点增长角等参数，最终获得了马赫杆高度的计算模型。计算结果表明：爆轰波相互作用发生马赫反射时，马赫杆的高度为零；随着爆轰波的传播，马赫杆的高度不断增加，马赫杆上的压力和速度降低，当其降低到CJ值时，马赫反射消失。
(2) 基于点火增长模型的Lee-Travel模型，带隔板装药爆轰波发生马赫反射后马赫杆高度的数值模拟结果与理论计算结果吻合较好，这说明利用Lee-Travel模型可以较准确地描述带隔板装药爆轰波的传播过程。
(3) 带隔板装药技术是一种利用爆轰波相互作用产生超压爆轰的装药技术。马赫杆高度模型的建立不仅丰富和完善了带隔板装药爆轰波马赫反射的传播规律，而且可为描述带隔板装药驱动药型罩的过程提供理论基础。
-
图 1 CFBG传感器测量爆轰波速度的结构示意图
Figure 1. Chirped fiber Bragg grating detonation velocity sensor system block diagram
图 4 标定CFBG传感器常数的实验布局图
Figure 4. Configuration for calibrating the constants of CFBG sensor
-
[1] Bdzil J B, Fickett W, Stewart D S. Detonation shock dynamics: A new approach to modeling multi-dimensional detonation waves[C]//Proceedings of the Ninth Symposium(International)on Detonation. Portland, OR, 1989: 730-742. [2] Aslam T D, Bdzil J B, Hill L G. Extensions to DSD theory: Analysis of PBX 9502 rate stick data[C]//Proceedings of the Eleventh International Detonation Symposium. Snowmass, CO, 1998: 21-29. [3] Cagliostro D J, Mandell D A, Schwalbe L A, et al. MESA 3-D calculations of armor penetration by projectiles with combined obliquity and yaw[J]. International Journal of Impact Engineering, 1990, 10(1): 81-92. https://www.sciencedirect.com/science/article/pii/0734743X90900506 [4] Walker F E. Initiation and detonation studies in sensitized nitromethane[J]. Acta Astronautica, 1979, 6(7): 807-813. https://www.sciencedirect.com/science/article/pii/0094576579900705 [5] Kasahara J, Hasegawa A, Nemoto T, et al. Performance validation of a single-tube pulse detonation rocket system[J]. Journal of Propulsion and Power, 2009, 25(1): 173-180. doi: 10.2514/1.37924 [6] Stewart D S, Yao J. The normal detonation shock velocity: Curvature relationship for materials with nonideal equation of state and multiple turning points[J]. Combustion and Flame, 1998, 113(1): 224-235. https://www.sciencedirect.com/science/article/pii/S0010218097001703 [7] Benterou J, Udd E, Wilkins P, et al. In-situ continuous detonation velocity measurements using fiber-optic Bragg grating sensors[C]//Proceedings of 34th International Pyrotechnics Seminar. Beaune, France, 2007: 309-322. [8] Udd E, Benterou J. High speed fiber optic grating sensor system: US 0028489[P]. 2009-01-29. [9] Benterou J, Udd E, Wilkins P, et al. In-situ continuous detonation velocity measurements using fiber-optic Bragg grating sensors[R]. LLNL-PROC-233137, 2007. [10] Udd E, Benterou J. Damage detection system with sub-microsecond resolution[C]//The 15th International Symposium on: Smart Structures and Materials & Nondestructive Evaluation and Health Monitoring. International Society for Optics and Photonics, 2008: 69330B-9. [11] Benterou J, Bennett C V, Cole G, et al. Internal detonation velocity measurements inside high explosives[C]//SPIE Defense, Security and Sensing. Orlando, FL, 2009. [12] Udd E, Benterou J, May C, et al. Review of high-speed fiber optic grating sensors systems[C]//Fiber Optic Sensors and Applications Ⅶ. Orlando, Florida, 2010: 76770B. [13] Benterou J, May C A, Udd E, et al. High Speed Measurements using fiber-optic Bragg grating sensors[C]//2011 SPIE Defense, Security, & Sensing Conference. Orlando, FL, 2011. [14] Rodriguez G, Sandberg G L, McCulloch Q, et al. Chirped fiber Bragg grating detonation velocity sensing[J]. Review of Scientific Instruments, 2013, 84(1): 015003. doi: 10.1063/1.4774112 [15] 王玉田, 曹丽芳.光纤Bragg光栅应力测量的研究[J].传感技术学报, 2007, 20(5): 1193-1195. http://www.cnki.com.cn/Article/CJFDTotal-CGJS200705053.htmWang Yu-fang, Cao Li-fang. Research of stress survey by fiber Bragg grating[J]. Chinese Journal of Sensors and Actuators, 2007, 20(5): 1193-1195. http://www.cnki.com.cn/Article/CJFDTotal-CGJS200705053.htm [16] 姜德生, 何伟.光纤光栅传感器的应用概况[J].光电子·激光, 2002, 13(4): 420-430. http://www.cnki.com.cn/Article/CJFDTotal-GDZJ200204025.htmJiang De-wei, He Wei. Review of applications for fiber Bragg grating sensors[J]. Journal of Optoelectronics·Laser, 2002, 13(4): 420-430. http://www.cnki.com.cn/Article/CJFDTotal-GDZJ200204025.htm [17] 张颖, 刘志国, 郭转运, 等.高灵敏度光纤光栅压力传感器及其压力传感特性的研究[J].光学学报, 2002, 22(1): 89-92. http://www.cqvip.com/Main/Detail.aspx?id=6382037Zhang Yin, Liu Zhi-guo, Guo Zhuan-yun. A high-sensitivity fiber grating pressure sensor and its pressure sensing characteristics[J]. Acta Optica Sinica, 2002, 22(1): 89-92. http://www.cqvip.com/Main/Detail.aspx?id=6382037 期刊类型引用(6)
1. 李昱锦，缪玉松，李斌，黄飞飞，葛帅，杨涛，吴佳怡. 对称双线性起爆下爆轰波的轴向传播特性. 爆破器材. 2025(01): 7-13 . 
百度学术2. 张朝平，张先锋，谈梦婷，侯先苇，熊玮，刘闯，顾晓辉. 聚能杆式射流侵彻混凝土和岩石靶体试验与数值模拟. 含能材料. 2023(08): 773-785 . 百度学术3. 李瑞，李伟兵，靳洪忠，王桂林，洪晓文，朱建军，李文彬，王晓鸣. 基于Jones-Wilkins-Lee状态方程的爆轰波相互作用参数理论分析. 兵工学报. 2019(03): 516-521 . 百度学术4. 潘建，张建军，侯云辉，冯雪磊. 带隔板装药EFP成型数值模拟与实验研究. 振动与冲击. 2019(18): 89-94+121 . 百度学术5. 潘建，张建军，侯云辉，刘燕宁. 带隔板装药爆轰波驱动药型罩数值模拟研究. 弹箭与制导学报. 2019(04): 149-154 . 百度学术6. 成乐乐，赵太勇，陈智刚，印立魁，王维占，郭光全，雷文星. 分层装药对杆式射流的影响分析. 爆破器材. 2018(04): 35-39 . 百度学术其他类型引用(5)
-
下载:
下载:
百度学术
计量
- 文章访问数: 3150
- HTML全文浏览量: 385
- PDF下载量: 303
- 被引次数: 11