Stress-testing method by PVDF gauge and its application in explosive test of porous material
-
摘要: 通过精确测量Hopkinson杆子弹速度,实现了对PVDF压力传感器的动态标定,测试数据线性度良好,误差不超过1.9%,得到动态灵敏度系数K=32.83 pC/N,采用500 g TNT对“钢板-泡沫铝-钢板”复合结构进行爆炸冲击加载,测量结构间应力波的传播情况。研究结果表明:电压测试信号可以较为准确地反映弹性波与塑性波的加载时间和传播速度,PVDF对弹性应力波段高频信号的动态响应灵敏准确,与理论数据的相对误差为3.5%。测得泡沫铝材料中塑性波的传播速度为590 m/s,A1-B1界面塑性波透射系数达到了0.53,远高于弹性波透射系数。从机理上对应力时程曲线中出现的特殊现象进行了阐述,为相关爆炸测试提供参考。Abstract: In this work, by accurately measuring the bullet speed of the Hopkinson bar, we succeeded in the dynamic calibration of the PVDF pressure sensor. The test data obtained have a good linearity and a deviation below 1.9%. The coefficient value K of dynamic sensitivity is 32.83 pC/N. By applying a blast loading of 500 g TNT to the composite structure of "steel-aluminum foam-steel", the propagation of the stress wave between its components was measured. The results show that the voltage signal can accurately represent the loading time and the propagation velocity of the elastic wave and plastic wave and the PVDF dynamic response of elastic stress wave with high-frequency signal is sensitive and accurate, with only a relative deviation of less than 3.5% from the theoretical data. The measured propagation velocity of the plastic wave in the aluminum foam is 590 m/s, the transmission coefficient of the plastic wave in A1-B1 interface is 0.53, which is much higher than that of the elastic wave. Special remarks were made concerning the unusual phenomena observed in the stress curve based upon our understanding of the mechanism, providing a reference for relevant explosion tests.
-
Key words:
- mechanics of explosion /
- explosion test /
- aluminum foam /
- PVDF sensor
-
近年来,炸药安全性受到了学者的重视。在多起意外事故中,炸药在运输途中跌落、钻孔操作时及与地面滑动等情况下发生爆炸[1],虽然尚未确定事故中准确机制,但摩擦很可能是最主要的点火机制。目前炸药摩擦一般是指狭义上的摩擦[2],在炸药摩擦过程中,伴随着众多的复杂现象,如塑性变形、相变、带电现象、化学反应、熔化和微动损伤等,而研究主要集中在热力学方面,即摩擦热带来的效应:摩擦将热量限制在移动的摩擦表面之间或在两表面之间的存在杂质的位置。当摩擦刺激达到足够的强度,温度将升到一定高度,这个温度足以使放热反应发生并超过放热过程,那么点火将可能发生,如存在约束条件也可能进一步转变成爆轰,对炸药安全性带来威胁。
为发展炸药摩擦安全性的评估方法,对炸药摩擦点火开展了许多相关研究。J.G.Glenn等[3]将炸药放置在圆柱管内,用活塞装置推动炸药使之在钢管内滑动,产生摩擦,称为IPFT(intense pressure and friction test)实验。在相同条件下,改变钢管内表面的光洁度,炸药发生了不同程度的反应。IPFT实验提供了摩擦点火的可重复实验。他们还建立了炸药与钢管摩擦放热的数学模型,并应用SPH方法模拟了动摩擦过程。A.Birk等[4]将COMP-B炸药约束在不同厚度的钢管中,活塞以不同的速率对它进行加载,使它发生反应。如果载荷速率较小,炸药点火经常在压力衰减至100 MPa开始,即在点火发生之前实际上只有很小的载荷,他们认为这表明是炸药与管壁的摩擦导致点火,实验显示摩擦效应是炸药点火的重要机制。在摩擦的数值模拟中,J.K.Dienes[5]假设摩擦面的温度达到熔点后就不再升高。D.Hoffman等[6]针对LX-04炸药-钢界面开展了摩擦系数测量,实验结果表明,由于摩擦热的生成,炸药表面会出现塑化和层间融化现象,形成摩擦弱化现象,即炸药-金属界面的摩擦系数随温度上升有减小现象,这个实验结果对安全性的数值模拟有效性有重要意义。Y.Q.Wu等[7]对HMX、RDX和PETN进行了摩擦实验和模拟,测量了各种晶体间摩擦的摩擦系数,给出升温公式,而且在实验图像中也发现了局部熔化区域。
我们考虑炸药的熔化效应,建立了一个炸药摩擦模型,对敏感炸药的摩擦感度进行了模拟研究[8-9]。本文中,对几种感度较弱炸药的摩擦感度进行数值模拟,以证明此模型的合理性和较广的适用性,并结合炸药活化能分析炸药在不同摩擦条件下感度的变化规律,以期能够为炸药摩擦安全性研究提供参考。
1. 摩擦点火模型
由于实验中观察到炸药的熔化现象,且许多炸药是低熔点(如TNT、DATB和NQ等炸药)[10],一般低于炸药的点火温度,所以建立的摩擦数学模型有必要加入熔化效应。参照文献[3]的炸药摩擦模型,建立了一维炸药摩擦点火模型[8-9]:摩擦功转换成热传入两个摩擦界面,保证摩擦面上温度相同,炸药点火应用经典的热爆炸理论[11],引入Arrhenius反应速率计算炸药化学反应热。忽略摩擦作用引起的炸药整体塑性变形导致炸药的升温。具体数学描述如下。
如图 1,炸药与某种材料介质之间有速度差,发生摩擦,炸药热反应扩散方程、材料热传导方程、在摩擦面上的条件分别为:
∂T1∂t=a1∂2T1∂x2+Qc1Ze−EaRT (1) ∂T2∂t=a2∂2T2∂x2 (2) x=L,q1+q2=q,T1=T2 (3) 式中:T为温度,ρ为密度,c为比热,Q为单位质量炸药的化学反应产生的热,a=k/ρc。Z为指前因子,Ea为活化能。q为在摩擦界面上产生的总热量,q1和q2分别为传入炸药和钢的热量。摩擦形成的总热量q=μpv,其中μ为动摩擦系数,p为炸药中的压力,v为钢与炸药之间的相对速度。
在摩擦模型中,加入熔化效应:在摩擦的开始阶段即未发生熔化阶段,摩擦力按照库伦摩擦定律计算。当摩擦持续造成炸药发生熔化时,根据热导理论,计算熔化界面移动,变为液体的炸药材质参数发生改变,但热导中不考虑液态炸药的对流效应。类似润滑效应,炸药摩擦系数会随着摩擦界面的温度升高而减小,可以认为摩擦系数随温度线性变化:
μ={(μ0−μm)Tm−TTm−T0+μmTs≤TmμmTs>Tm (4) 式中:μ0为初始温度下的摩擦系数,Tm为炸药的熔点。摩擦感度实验中,对炸药施加了一定的压力,金属与炸药为密实接触,因此μ0 =1。μm为炸药发生熔化后的摩擦系数。炸药熔化后,炸药中出现液体与固体的分界面。随着摩擦界面温度的升高,液体与固体的分界面向炸药内部移动。忽略由于相变引起的炸药密度的变化,相变的分界面方程为[12]:
kl∂T1∂x−ks∂Ts∂x=ρLdSdt (5) 式中:下标l表示液体,s表示固体。S是相变分界面,L是熔化潜热。需要补充说明的是,当炸药熔化后,一部分是液态,可认为两个摩擦面的温度不再相同。至此,建立了考虑熔化效应的一维炸药摩擦点火模型。
此模型还可以描述多种材料之间的炸药摩擦点火。如果摩擦时炸药外层加入其他材料(如黏合剂等),再与摩擦材料发生摩擦,如果中间的材料具有低熔点,同样也需要考虑测此材料的熔化效应。编制了一维炸药摩擦点火热传导程序,采用模块化设计,能够计算多种物质的热传导和熔化问题,通过了正确性验证,可以应用于多介质的炸药摩擦感度实验的数值模拟研究。
2. 数值模拟
2.1 数值模拟结果与实验对比
为了比较炸药感度强弱顺序,建立了炸药摩擦感度实验来测量炸药感度。炸药摩擦感度实验[10]的通常做法是,将炸药放在实验装置的两个滑柱中,用摆锤击打上滑柱,炸药与金属滑柱发生摩擦,观察是否发生爆炸(含燃烧,分解),以发生爆炸的概率表示摩擦感度。实验测量的结果,见表 1。
表 1 单质炸药的摩擦感度[10]Table 1. Friction sensitivity of simple explosive炸药 η/% 炸药 η/% HMX 92~100 TNT 4~6 PETN 92, 100 TATB 0~4 RDX 76±8 DATB 0~4 Tetryl 12 NQ 0 图 2为简化的炸药摩擦感度实验图。由于炸药摩擦点火过程很短且建立的是一维模型,可以忽略炸药变形和滑动过程炸药自身制备形式对摩擦的影响效应,认为炸药处在两个金属块中间。上边的金属块有相对速度3.8 m/s和位移,并施加了一定的压力390 MPa(由于实验中滑柱直径是炸药直径的10倍,所以炸药受到压力为表压的100倍)。由于实验中炸药与金属之间的位移远小于滑柱的直径,可以忽略侧向边界的热传导效应。由于炸药够厚,可以忽略下方金属块影响。至此,用一维摩擦点火数值模拟程序对4种摩擦感度较弱的炸药进行计算,炸药和金属摩擦参数见表 2。
材料 ρ/(kg·m-3) c/(J·kg-1·K-1) k/(J·m-1·s-1·K-1) km/(J·m-1·s-1·K-1) Q/(MJ·kg-1) Ea/(kJ·mol-1) Z/s-1 L/(kJ·kg-1) Tm/K TNT 1 654 1 062.76 0.260 3 0.266 20 1.26 144 2.51×1011 521 353.9 DATB 1 834 1 092.05 0.259 0 0.259 70 1.26 194 1.17×1015 344 559.0 TATB 1 930 899.58 0.543 9 0.182 34 2.51 251 3.18×1019 652 721.0 NQ 1 690 1 125.52 0.410 0 0.626 40 2.09 87.4 2.87×107 402 518.0 Fe 7 840 465 49.8 - - - - - - 由于部分炸药(如TATB和NQ)往往单个文献中的熔化后的参数(如熔化潜热)不全,因此结合两个文献资料的数据,所以可能存在部分矛盾,如熔点在不同的实验测量中结果不一致。本文中数值模拟尽量保证取同一资料的参数。
炸药熔化后,热传导系数发生改变,熔化后的热传导系数通过液体分子理论公式近似计算得到[14]。熔化后的摩擦系数μm难于确定,4种炸药TNT、DATB、TATB和NQ的摩擦系数分别取为0.15、0.15、0.13和0.15, 进行了计算。初始时,炸药和钢的温度均取为300 K。
图 3是4种炸药和钢摩擦界面上的温度随时间的变化曲线图。从图中可以看出,4种炸药温度随着时间增长而不断升高,到达某个时间后温度曲线会发生转折,急剧上升,接近垂直,温度达到极大值,可认为炸药发生了点火,此时就是点火时间(转折处的温度即为点火的最低温度)。4种炸药的点火温度不相同:TATB和DATB较低,NQ和TNT较高。从图 4可以看到,4种炸药由于不同的热学性能,温度上升曲线的斜率均不相同。点火最低温度均大于炸药的熔点,所以4种炸药均先于点火发生了熔化:熔化后物质吸热,炸药物态参数(仅考虑热传导系数)发生了变化,所以温度上升速率会发生变化,温度曲线在各自熔点附近有变化或转折。模型假设摩擦系数随时间线性变化,所以模拟结果中温度曲线表现为随时间较平顺的变化,是符合物理规律的。由计算结果可知,4种炸药点火时间分别为:0.65 ms(TNT)、0.77 ms(TATB)、0.83 ms(DATB)和1.37 ms(NQ),炸药的点火时间越短,则感度越强,所以摩擦感度的强弱顺序应为:η(TNT)>η(TATB)≈η(DATB)>η(NQ)。实验结果(见表 1)中,点火概率越高,越容易点火,感度越强。计算结果与实验结果的炸药强弱顺序相同,两者符合较好。这说明了此摩擦点火模型能够较好模拟感度较低的4种炸药的感度实验,而以往模拟得到4种炸药的点火时间分别为:0.25 ms(PETN)、0.40 ms(HMX)、0.42 ms(RDX)和0.46 ms(Tetryl),感度强弱顺序也符合实验结果,这8种炸药的计算结果证明了摩擦模型的正确性和较广的适应性。而从图 4可以看出,4种炸药的温度分布曲线在熔点处均出现转折,根据模型假设,转折幅度主要取决于固态和液态炸药的热传导系数之差以及熔化潜热。模型假设摩擦系数随时间线性变化,模拟结果中温度随时间较平顺变化,符合物理规律。从图 4还可以看到,发生点火时,熔化界面已经向炸药内部推进了一定距离,熔化区的厚度(摩擦面到熔化面的距离)也有区别,熔点高的薄(TATB的为1.07 μm),熔点低的厚(NQ的为17.76 μm),说明不同炸药熔化面移动速度也不同,影响因素可能包括了热传导系数和熔化潜热。
模型中可调的参数为熔化后的摩擦系数μm(无实验数据),4种炸药中仅TATB略小,其他均相同,说明模型对此参数依赖较小,模型具有较好的参数适应性。
2.2 数值模拟结果分析
由于炸药摩擦的带源项热传导方程无解析解,炸药的活化能反映了炸药反应的快慢,所以利用源项中包含活化能的炸药热分解反应速率k=Ze−EaRT,来分析炸药点火时间快慢因素。图 5为8种常用炸药的热分解速率,由图可见,热分解速率对数与时间倒数,呈线性关系。
8种炸药可以根据速率大小分为两组:第1组为HMX、RDX、PETN和Tetryl,温度越高其反应速度较快,第2组为TNT、TATB、DATB和NQ,温度越高其反应速度较慢。这种分组也符合摩擦实验的结果。第2组炸药速率对数直线在高温时有交点,速率大小顺序出现变化:NQ在温度低时速率快,温度高时速率慢;TATB温度低时速率慢,温度高时速率快;DATB处于二者之间。由于摩擦点火可能在不同的高温下发生,所以认为3种炸药在强度不同的摩擦条件下可能会出现点火顺序发生交换的情况,需要经过数值模拟验证。
利用炸药摩擦点火模型程序计算了3种炸药TATB、DATB和NQ在4种摩擦(压力为400、40、4和0.4 MPa,摩擦速度4 m/s)下的点火时间、点火温度和熔化区厚度,如图 6所示,具体数值见表 3。
表 3 3种炸药的点火时间、点火温度和熔化区厚度Table 3. Ignition time, temperature and molten zone thickness for three explosives炸药 p/MPa ti/s Ti/K hm/μm TATB 400 0.52×10-3 860 1.03 40 0.026 770 3.48 4 0.65 690 0 0.4 28 640 0 DATB 400 0.74×10-3 940 5.71 40 0.075 790 4. 11 4 4.4 690 203.5 0.4 206 620 712. 1 NQ 400 2.6 1 300 22.71 40 0.095 800 85. 87 4 2.7 640 242. 4 0.4 19 540 3 517.3 通过以上数据可看出:在高压下,TATB点火最快,其次是DATB,而NQ点火最慢;反之,在低压下,DATB最快,TATB居中,而NQ反而最快了。点火顺序随着摩擦强度变化,点火顺序发生了交换,证实了以上假设。同时,考察点火温度和熔化区厚度,发现TATB与DATB点火温度接近,随着压力减少而降低,而NQ的随着压力降低下降更快,主要是因为其反应速率随压力下降较快。有趣的是,3种炸药的熔化区厚度基本呈一个平行的状态,而不是交叉的关系,考察其原因,不仅受反应速率影响,还受到熔点和熔化吸能的影响。以上结果说明,炸药摩擦感度强弱顺序可能会受到摩擦条件的影响,目前的摩擦感度实验尚不能完全说明各种炸药的摩擦感度。
3. 结论
利用摩擦点火模型,对炸药的摩擦感度实验中摩擦感度较低的4种炸药进行了数值模拟。数值模拟结果中的点火时间快慢顺序代表摩擦感度强弱顺序,与实验结果符合较好,说明此理论模型考虑了炸药熔化及摩擦系数随温度变化是合理的,此模型能够适应广范围的炸药,具有较强的适应性,期望将来能对更多炸药适用,如混合炸药改进模型。
分析3种炸药热分解反应速率,发现在不同温度下,热分解反应速率顺序会发生交换,而数值模拟验证了此观点:在不同摩擦条件下,点火时间发生交换,即摩擦感度顺序发生交换,由此说明摩擦感度实验无法完全表示炸药感度强弱的顺序,即摩擦感度受到摩擦强度影响。
-
表 1 Hopkinson杆冲击测试数据
Table 1. Measurement data of Hopkinson bar impacting
传感器1 传感器2 v/(m·s-1) σ/MPa Ωmax/pC v/(m·s-1) σ/MPa Ωmax/pC 8.03 56 2 300 9.67 68 2 640 6.12 43 1 820 7.24 51 2 170 9.31 65 2 560 5.16 36 1 690 10.97 77 2 880 14.60 102 3 880 13.35 93 3 490 20.76 145 5 340 14.69 103 3 910 24.38 171 5 850 20.50 143 5 240 24.11 169 6 060 表 2 实验相关参数(炸高15 cm)
Table 2. Related parameters(blasting height: 15 cm)
泡沫铝 δ/% ρ/(g·cm-3) d/mm 传感器 布设位置 R/Ω B1 84.5 0.418 5 2.5 1 A1-B1 42.68 B2 80.8 0.518 4 2.0 2 B1-B2 43.02 B3 77.6 0.604 8 2.0 3 B2-B3 42.28 -
[1] Kawai H. The piezoelectricity of polyvinyldene fluoride[J]. Journal of Applied Physics, 1969, 8(7):975-976. doi: 10.1143/JJAP.8.975 [2] Graham R A, Anderson M U, Bauer F, et al. Piezoelectric polarization of the ferroelectric polymer PVDF from 10 MPa to 10 GPa: Studies of loading-path dependence[C]//Shock Compression of Condensed Matter-1991. Williamsburg, Virginia, 1992: 883-886. https://www.sciencedirect.com/science/article/pii/B978044489732950203X [3] Bauer F. Advances in piezoelectric PVDF shock compression sensors[C]//10th International Symposium on Electrets. Greece: IEEE Stands Office, 1999: 647-650. http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=832129 [4] Bauer F. PVDF shock sensors: Applications to polar materials and high explosives[J]. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 2000, 47(6):1448-1454. doi: 10.1109/58.883534 [5] Sokhanvar S, Dargahi J, Packirisamy M. Influence of friction on piezoelectric sensors[J]. Sensors & Actuators A Physical, 2008, 141(1):120-128. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=e7e4fd12958a1a28042606d90f80f6c4 [6] Shirinov A V, Schomburg W K. Pressure sensor from a PVDF film[J]. Sensors & Actuators A Physical, 2008, 142(1):48-55. doi: 10.1007-s10725-010-9536-z/ [7] 席道瑛, 郑永来.PVDF压电计在动态应力利测量中的应用[J].爆炸与冲击, 1995, 15(2):174-179. http://www.bzycj.cn/article/id/10573Xi Daoying, Zheng Yonglai. Application of PVDF gauges to dynamical stress measurements[J]. Explosion and Shock Waves, 1995, 15(2):174-179. http://www.bzycj.cn/article/id/10573 [8] 李焰, 张向荣, 谭红梅, 等.国产PVDF压电薄膜的冲击加载及卸载响应研究[J].高压物理学报, 2004, 18(3):261-266. doi: 10.3969/j.issn.1000-5773.2004.03.011Li Yan, Zhang Xiangrong, Tan Hongmei, et al. Response of homemade PVDF piezofilm under shock loading and unloading[J]. Chinese Journal of High Pressure Physics, 2004, 18(3):261-266. doi: 10.3969/j.issn.1000-5773.2004.03.011 [9] 赵继波, 谭多望, 张远平, 等.PVDF计在水中爆炸近场压力测试中的应用[J].火炸药学报, 2009, 32(3):1-4. doi: 10.3969/j.issn.1007-7812.2009.03.001Zhao Jibo, Tan Duowang, Zhang Yuanping, et al. Application of PVDF sensor in underwater explosive shock wave measurement at near-field[J]. Chinese Journal of Explosives & Propellants, 2009, 32(3):1-4. doi: 10.3969/j.issn.1007-7812.2009.03.001 [10] 王伟, 李小春.不耦合装药下爆炸应力波传播规律的试验研究[J].岩土力学, 2010, 31(6):1723-1729. doi: 10.3969/j.issn.1000-7598.2010.06.008Wang Wei, Li Xiaochun. Experimental study of propagation law of explosive stress wave under condition of decouple charge[J]. Rock and Soil Mechanics, 2010, 31(6):1723-1729. doi: 10.3969/j.issn.1000-7598.2010.06.008 [11] 任新见, 陈虎林.PVDF传感器在爆炸近区超压测量中的应用研究[J].振动与冲击, 2012, 31(12):146-149. http://d.old.wanfangdata.com.cn/Periodical/zdycj201212029Ren Xinjian, Chen Hulin. Application of PVDF sensor in shock wave measurement of near zone of explosion[J]. Journal of Vibration and Shock, 2012, 31(12):146-149. http://d.old.wanfangdata.com.cn/Periodical/zdycj201212029 [12] 王永刚, 胡时胜, 王礼立.爆炸荷载下泡沫铝材料中冲击波衰减特性的实验和数值模拟研究[J].爆炸与冲击, 2003, 32(6):516-523. doi: 10.3321/j.issn:1001-1455.2003.06.006Wang Yonggang, Hu Shisheng, Wang Lili. Shock attenuation in aluminum foams under explosion loading[J]. Explosion and Shock Waves, 2003, 32(6):516-523. doi: 10.3321/j.issn:1001-1455.2003.06.006 [13] 田杰.泡沫铝的冲击波衰减和抗爆震特性研究[D].合肥: 中国科学技术大学, 2006: 50-66. http://cdmd.cnki.com.cn/Article/CDMD-10358-2006102746.htm [14] 王礼立.应力波基础[M].北京:国防工业出版社, 1985. [15] 胡亚峰, 顾文彬, 刘建青, 等.应力波在泡沫铝夹层三明治板中传播规律的数值研究[J].爆破器材, 2014, 43(5):15-23. doi: 10.3969/j.issn.1001-8352.2014.05.004Hu Yafeng, Gu Wenbin, Liu Jianqing, et al. Numerical study of stress wave propagation rules in aluminum foam sandwich plates[J]. Explosive Materials, 2014, 43(5):15-23. doi: 10.3969/j.issn.1001-8352.2014.05.004 期刊类型引用(8)
1. 袁良柱,陈美多,谢雨珊,陆建华,王鹏飞,徐松林. 细观非连续介质的应力波传播研究. 爆炸与冲击. 2024(09): 50-61 . 本站查看
2. 刘巍,马宏昊,徐钦明,姚象洋,赵勇,杨科,杨辉,沈兆武. 浅埋炸药加载下含空穴泡沫铝夹芯板动态响应机制. 兵工学报. 2023(12): 3613-3621 . 百度学术
3. 杨建兴,张江涛,乔炎亮,张梅,翟鹏程. 硅橡胶动态压缩性能SHPB测试方法及其本构模型研究. 高分子通报. 2021(04): 27-34 . 百度学术
4. 孙权,于洋,陈宝成. 基于PVDF的压电传感器动态压力标定. 黑龙江科学. 2021(14): 24-28 . 百度学术
5. 舒波超,李卫东,黄遐,车志刚. 铝合金材料激光冲击喷丸力学响应有限元建模. 航空制造技术. 2020(03): 59-66 . 百度学术
6. 徐景林,顾文彬,刘建青,刘欣,胡云昊,韩阳明,符磊. PVDF的标定及在泡沫铝爆炸冲击实验中的应用. 中国测试. 2018(10): 101-106 . 百度学术
7. 王明洋,徐子章,岳松林,邱艳宇,刘家贵. PVDF压力计在相似材料动态力学参数测试中的应用. 北京理工大学学报. 2017(10): 991-997 . 百度学术
8. 曹州权,谢小柱,陈蔚芳,魏昕,胡伟,任庆磊. 液体辅助激光加工过程中压力检测及应用研究进展. 光电工程. 2017(04): 381-392+470 . 百度学术
其他类型引用(2)
-