Study on mechanical properties and damage characteristics of booster explosives under static compression
-
摘要: 为研究聚黑-14C(JH-14C)传爆药静态压缩力学性能及损伤特性,开展准静态压缩实验,获得了不同应变率下的应力-应变曲线,建立了描述不同应变率下JH-14C力学行为的非线性本构模型;利用扫描电镜(SEM)对回收试样进行细微形貌观测,获得了准静态压缩JH-14C损伤特性的表征。结果表明:JH-14C压缩强度随应变率的升高而提高;实验与计算结果对照验证了本构模型的有效性;准静态压缩实验中,JH-14C主要损伤模式为脱湿和穿晶断裂。Abstract: In order to study the static compression mechanical properties and damage characteristics of the JH-14C booster explosive, quasi-static compressive experiments were performed on a testing machine equipped with an environmental chamber (INSTRON). According to the GJB 770B–2005 powder test method, dimensions of the cylindrical specimen were set as
∅ 12.5 mm×12.5 mm in the static compressive experiments. During compression, only one extensometer was used. All experiments were performed at a crosshead speed of 0.012 5, 0.062 5, 0.125, 0.625 and 1.25 mm/s at room temperature (25 °C), which led to a nominal strain rate of 0.001, 0.005, 0.01, 0.05 and 0.1 s−1, respectively. The average stress-strain values and standard deviations were calculated using five replicable experiments for each condition. The experimental results were compared with X0242 and PBX-9501, and the mechanical properties of JH-14C were analyzed. According to the mechanical properties of JH-14C at low strain rates, the original Ramberg-Osgood constitutive relationship was modified, and a nonlinear constitutive model including the strain rate term was established to describe the mechanical behavior of JH-14C at low strain rates. The micro morphologies of the recovered samples was observed by a scanning electron microscope (SEM) and compared with that of PBX-9501. The damage mode was analyzed to characterize the damage characteristics of JH-14C under quasi-static compression. The results show that the compressive strength of JH-14C increases with the increase of strain rate. The validity of the constitutive model was verified by comparing the experimental and calculated results. In the quasi-static compression experiments, the energetic particles and the binder were debonded. With the increase of the pressure, the original crack and the micro-crack formed by debonding on the energetic particles were converged and coalesced to form a macro crack, which led to the rupture and failure of the explosive.-
Key words:
- JH-14C /
- booster explosive /
- quasi-static /
- strain rate /
- constitutive model /
- damage mode
-
韧性金属柱壳在高能炸药驱动下的膨胀断裂(简称外爆)是一个快速而复杂的动态过程。从炸药起爆到柱壳破坏只需数十微妙, 在此过程中柱壳经受了强加载、高应变率、大变形和复杂的加卸载。
R.W.Gurney[1]给出了外爆下破片速度的经验公式并被广泛使用。G.I.Taylor[2]根据高速摄影的实验图像, 提出了著名的Taylor断裂判据。C.R.Hoggatt等[3]在回收的柱壳破片内壁面附近发现了绝热剪切带, 在Taylor思想的基础上, 提出外爆加载下在柱壳内壁面先生成剪切带并逐步向外扩展, 而后在拉应力控制下, 裂纹沿着剪切带失稳面由外向内传播的破坏机制。胡八一等[4]对外爆中发生的剪切带进行了系统研究。胡海波等[5]分析了柱壳破坏中出现的“单旋”剪切现象。汤铁钢等[6]开展了一系列实验, 研究柱壳破坏的应变率效应。李永池等[7]提出损伤演化方程, 模拟了中心线起爆方式下柱壳的破坏过程。
实验研究中, 炸药的起爆方式主要为单点起爆和平面起爆。但在理论分析和数值模拟中, 多数研究按照中心线起爆处理, 深入研究不同起爆方式对柱壳破裂影响的文章较少。本文中, 研究这3种起爆方式对柱壳破裂的影响。
1. 方法
外爆研究中, 炸药的起爆方式有3种:单点起爆、平面起爆和中心线起爆。单点起爆是指在炸药的一端插入雷管, 点燃雷管引爆炸药, 如图 1所示。平面起爆是指炸药一端各点同时起爆, 中心线起爆是指炸药柱中心线各点同时起爆。
对不同的起爆方式, 金属柱壳所经受的加载历程显然不同。但无论在理论研究还是在数值分析中, 为了便于理论推导和细化网格, 大多忽略结构的轴向影响, 把柱壳的破坏简化为二维平面应变问题, 并采用中心点起爆方式起爆炸药[6-8]。这种简化对中心线起爆的柱壳中部是合理的, 但大部分实验采用的是单点起爆和平面起爆[4-6]。此时这种简化是否也合理, 不同的起爆方式下柱壳所经受的加载历程又有何不同?本文中将对这些问题进行相应的数值模拟, 分析起爆方式对柱壳破裂的影响。
分别对3种不同的起爆方式, 使用Ls-Dyna软件模拟炸药爆轰驱动柱壳膨胀的过程。为了简化计算, 只计算了结构的四分之一, 计算模型如图 2。网格尺度为1 mm×1 mm×1 mm。
金属柱壳为4340钢, 强度采用Johnson-Cook模型:
Y=(A+Bˉεnp)(1+Cln˙ˉεp˙ˉε0)(1−(T−T0Tm−T0)m) (1) 式中:A、B、n、C和m是材料参数; T0为初始温度, Tm为熔点,
和
分别为等效塑性应变和等效塑性应变率。4340钢的参数分别为:密度ρ=7 830 kg/m3, 体积模量K=164 GPa, 剪切模量
炸药为固态TNT, 爆轰产物采用JWL状态方程:
p=A(1−ωR1V)e−R1V+B(1−ωR2V)e−R2V+ωEV (2) 式中:V为比体积, p为压力。TNT炸药的参数分别为:密度ρ=1 640 kg/m3, 爆速D=6 930 m/s, 爆压pCJ=27 GPa, A=374 GPa, B=3.23 GPa, R1=4.15, R2=0.95, ω=0.3。
2. 结果与分析
柱壳轴向长度为12 cm, 分别提取柱壳内表轴向2、4、6、8和10 cm处的径向位移曲线和爆轰产物对柱壳的压力曲线, 处理数据得到炸药对柱壳内表初始单位面积的累积做功曲线。为了便于比较分析, 对累积做功曲线进行平移, 使爆轰产生的压力到达柱壳内表相应位置的时刻为零时刻。
2.1 平面起爆
平面起爆得到的曲线如图 3所示。柱壳内表所受压力载荷峰值在2 cm处由于边界效应的影响略低, 其他各处的压力峰值接近。虽然各处压力峰值接近, 但卸载速度不同。接近起爆端处, 直接接触空气, 压力到达峰值后卸载速度快, 炸药对柱壳的做功也小于远离起爆点处。在20 μs时, 8 cm处的做功是2 cm处的约1.73倍。
2.2 单点起爆
单点起爆得到的曲线如图 4所示。距离起爆点越远, 柱壳内表所受压力载荷峰值越小。因此在初期(约2 μs以前), 做功也随距离起爆点的位置越远而越小。但是由于距离起爆点近的位置, 压力卸载速度快, 因此在后期, 做功随距离起爆点的位置越远而越大。在20 μs时, 8 cm处的做功是2 cm处的约1.56倍。
2.3 中心线起爆
中心线起爆得到的曲线如图 5所示。2和10 cm处的曲线几乎重合, 4、6和8 cm处的曲线几乎重合, 各处的内表压力峰值接近。在6.2 μs(压力曲线中为9.0 μs)时, 两侧稀疏波到达2和10 cm处, 此时该两处的压力, 特别是单位面积做功明显开始低于其他各处。炸药由中心线起爆后, 沿径向向外爆轰, 到达与柱壳的接触面, 冲击波在接触面反射, 向中心线会聚, 而后再向外反射传播, 在13.8 μs(压力曲线中为16.6 μs)时, 再次到达与柱壳的接触面, 二次加载柱壳。此时, 压力升高, 单位面积的做功也快速增加。这种重加载在平面起爆和单点起爆中没有出现。在20 μs时, 8 cm处的做功是2 cm处的约1.27倍。
2.4 不同起爆方式的差别
分别取平面起爆、单点起爆和中心线起爆下柱壳内表轴向8 cm处的数据, 对比3种不同的起爆方式下爆轰气体对柱壳作用的差异, 得到的曲线如图 6所示。
由柱壳内表的压力曲线可以看出, 不同起爆方式下柱壳内表的压力载荷历程差异显著。中心线起爆对应的压力峰值最大, 其次是单点起爆, 平面起爆下的压力峰值最小。达到峰值压力后, 不同起爆方式下压力的卸载速度也存在明显差异, 峰值压力最高的中心线起爆对应的压力卸载速度最快, 其次是单点起爆, 平面起爆下柱壳内表压力卸载速度最慢。根据应力波理论分析可知, 该压力峰值的差别是由爆轰波对金属柱壳的入射角不同所造成的。中心线起爆方式下, 爆轰波对金属柱壳正入射, 爆轰波阵面平行于柱壳内表面。单点起爆和平面起爆下, 爆轰波对金属柱壳均为斜入射, 而平面起爆下入射角最大, 因此压力峰值也最小。
从做功曲线可以看出, 对轴向8 cm处的数据, 在20 μs时中心线起爆比平面起爆下的做功低3.21%。即到过程末期, 不同起爆方式引起的总做功差异不大。但是做功的累积过程却存在显著差异。在10 μs时, 中心线起爆比平面起爆下的做功低19.5%。
由以上分析可以看出, 不同起爆方式下炸药对柱壳内表的压力、做功历史差别明显。即使同一起爆方式下, 轴向不同位置的压力、做功历史也存在显著不同。因此, 在数值模拟和理论分析中, 必须考虑起爆方式和轴向位置的影响。
3. 半解耦数值方法及分析
对外爆的数值模拟中, 炸药网格占总网格的比例很大, 计算效率低下[7-8]; 另外, 粗糙的网格划分不能描述柱壳典型破坏中出现的剪切带(实验中剪切带宽度约50 μm), 甚至不能区分两个相邻的剪切带(实验中剪切带间距约1 mm)[5-6]; 此外起爆方式和轴向位置对柱壳破裂带来的影响也不容忽视。基于此, 本文中提出一种炸药和柱壳计算半解耦的数值计算方法。
3.1 半解耦方法
首先, 对三维模型(包括炸药和柱壳)划分较粗网格, 不引入金属柱壳的破坏模型, 得到爆轰产物对柱壳内表的压力曲线。然后, 去除炸药, 作平面应变假定进行二维计算, 划分精细网格, 将粗网格得到的内表压力曲线施加于柱壳内表面, 同时引入金属柱壳的破坏模型。这样, 一方面充分细化柱壳网格使之能模拟剪切带等破坏现象, 另一方面可以准确研究不同起爆方式下柱壳不同位置的破坏。
3.2 断裂能量耗散
设柱壳破坏过程中形成的剪切带条数为N, 产生单位面积剪切断裂面耗能为λ, 柱壳向外运动的速度为v, 柱壳的厚度为δ, 密度为ρ, 柱壳内径为r, 剪切带间距为
则形成新断裂面的断裂耗能与柱壳的动能比为:
EsEk=√2λδN0.5ρδ×2πrv2=2λρLv2 (3) 汤铁钢等[6]对图 2所示尺寸的装置进行实验, 使用高速摄影技术拍摄了整个膨胀过程, 使用的金属材料是钢, 密度ρ=7 830 kg/m3, 测得剪切带间距L≈1 mm, 末速度v≈1.8 km/s, 单位表面积的表面能[8]约16 KJ/m2。计算得Es/Ek≈0.126%, 可见断裂面的表面能所引起的能量重分配对柱壳速度的影响微乎其微, 完全可以忽略不计。
3.3 柱壳膨胀断裂过程中的塑性应变能
分析柱壳在膨胀过程中的塑性应变能, 设初始时刻柱壳内径为r0, 外径为R0, 末状态柱壳内径为rf, 外径为Rf。设R0/r0=p, Rf/r0=q。忽略弹性应变, 做平面应变假定, 即εr=-εθ。r方向总应变
则总塑性功为:
Wt=L∫2π0dθ∫R0r02σεrrdr=Mσ(q2lnq2−p2lnp2−(q2−p2+1)ln(q2−p2+1))ρm(p2−1) (4) 式中:L为柱壳的轴向长度, ρm为金属的密度, σ为金属流动应力, 假设变形过程中流动应力保持恒值
为柱壳总质量。
在图 2构型中,
在柱壳贯穿破坏时刻, 取q=2[2-3], 则总塑性功Wt≈75 kJ。应该注意到, 在柱壳破坏开始之前, 已经消耗大量塑性功, 在破坏开始时刻, 取q=1.7, 则破坏开始前塑性功Wa=52 kJ。因此材料破坏过程中, 如不引入破坏模型, 带来的塑性功Wi=Wt-Wa=23 kJ。柱壳动能Ek=Mv2/2=842 kJ。总塑性功与总动能的比为8.9%, 但不引入破坏模型带来的塑性功与柱壳总动能的比仅为2.7%。所以, 破坏过程中的塑性功所引起的能量分配对柱壳速度的影响较小, 可忽略。
3.4 半解耦法计算结果
对断裂耗能和塑性应变能的分析表明, 破坏产生新断面的能量耗散和破坏过程中的塑性应变能对柱壳动能的影响较小, 可忽略。因此半解耦数值方法是合理的, 有理论依据的。采用此半解耦数值方法, 模拟图 2结构在单点起爆方式下轴向6 cm处壳体的破坏过程, 网格尺度为90 μm×90 μm。局部放大破坏图像如图 7所示, 它成功地再现了炸药驱动下柱壳膨胀断裂过程中出现的绝热剪切带。
4. 总结
采用数值模拟方法, 研究了3种不同起爆方式下爆轰气体产物对韧性金属柱壳内表的加载历程, 发现不同起爆方式下柱壳内表所经受的压力历程差别巨大, 中心线起爆峰值压力最大, 单点起爆其次, 平面起爆最小。理论分析认为, 造成这种差别的原因是爆轰波对柱壳的入射角不同。即使在同一起爆方式下, 柱壳内表轴向不同点的压力载荷历程也存在巨大差异。因此, 在理论分析和数值模拟中, 必须考虑起爆方式和轴向位置的影响。提出一种模拟炸药驱动下韧性金属柱壳膨胀破坏的半解耦数值方法, 该方法可以实现对柱壳精细网格划分条件下的破裂模拟, 并能研究不同起爆方式的影响。通过理论分析说明了这种半解耦方法的合理性。
-
表 1 模型拟合结果的相关系数
Table 1. Correlation coefficients of the model fitting results
˙ε/s−1 R2 0.1 0.990 5 0.05 0.997 5 0.01 0.998 7 0.005 0.992 5 0.001 0.984 2 -
[1] 邓琼, 叶婷, 苗应刚. 基于Hopkinson压杆实验技术研究火工品及含能材料的抗高过载能力 [J]. 火炸药学报, 2009, 32(6): 66–70. DOI: 10.14077/j.issn.1007-7812.2009.06.019.DENG Q, YE T, MIAO Y G. Study on overloading-resistibility of initiator and energetic materials based on the technique of Hopkinson pressure bar [J]. Chinese Journal of Explosive & Propellants, 2009, 32(6): 66–70. DOI: 10.14077/j.issn.1007-7812.2009.06.019. [2] RAE P J, PALMER S J P, GOLDREIN H T, et al. Quasi-static studies of the deformation and failure of PBX 9501 [J]. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2002, 458(2025): 2227–2242. DOI: 10.1098/rspa.2002.0967. [3] RAE P J. Quasi-static studies of the deformation, strength and failure of polymer bonded explosives [D]. Cambridge, UK: University of Cambridge, 2000. [4] RAE P J, GOLDREIN H T, PALMER S J P, et al. Quasi-static studies of the deformation and failure of β-HMX based polymer bonded explosives [J]. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2002, 458: 743–762. DOI: 10.1098/rspa.2001.0894. [5] HEIDER N, STEINBRENNER A, AURICH H, et al. A method for the determination of the viscoelastic relaxation function of reactive materials [J]. The European Physical Journal Special Topics, 2016, 225(2): 397–407. DOI: 10.1140/epjst/e2015-66666-x. [6] HEIDER N, STEINBRENNER A, WEIDEMAIER P, et al. Modelling the mechanical behaviour of PBX KS32 [C] // 44th Annual Conference ICT. Karlsruhe, Germany: Fraunhofer Institute for Chemical Technology, 2013. [7] DIENES J K, ZUO Q H, KERSHNER J D. Impact initiation of explosives and propellants via statistical crack mechanics [J]. Journal of the Mechanics and Physics of Solids, 2006, 54(6): 1237–1275. DOI: 10.1016/j.jmps.2005.12.001. [8] BENNETT J G, HABERMAN K S, JOHNSON J N, et al. A constitutive model for the non-shock ignition and mechanical response of high explosives [J]. Journal of the Mechanics and Physics of Solids, 1998, 46(12): 2303–2322. DOI: 10.1016/S0022-5096(98)00011-8. [9] 李硕, 袁俊明, 刘玉存, 等. 聚黑-14C的传爆装置冲击起爆实验及数值模拟 [J]. 火炸药学报, 2016, 39(6): 63–68,79. DOI: 10.14077/j.issn.1007-7812.2016.06.011.LI S, YUAN J M, LIU Y C, et al. Experiment and numerical simulation of shock initiation of JH-14C detonation device [J]. Chinese Journal of Explosives & Propellants, 2016, 39(6): 63–68,79. DOI: 10.14077/j.issn.1007-7812.2016.06.011. [10] 董理赢. 引信传爆序列殉爆反应特性研究 [D]. 太原: 中北大学, 2020. DOI: 10.27470/d.cnki.ghbgc.2020.000896.DONG L Y. Study on the characteristics of the sympathetic detonation reaction in fuze explosive trains [D]. Taiyuan, China: North University of China, 2020. DOI: 10.27470/d.cnki.ghbgc.2020.000896. [11] 张子敏, 许碧英, 仲凯, 等. 冲击载荷下JH-14C传爆药的动态响应实验研究 [J]. 火炸药学报, 2010, 33(1): 57–59,63. DOI: 10.14077/j.issn.1007-7812.2010.01.018.ZHANG Z M, XU B Y, ZHONG K, et al. Experimental study on the dynamic response of booster explosive JH-14C under impact load [J]. Chinese Journal of Explosives & Propellants, 2010, 33(1): 57–59,63. DOI: 10.14077/j.issn.1007-7812.2010.01.018. [12] 张子敏, 许碧英, 贾建新, 等. 基于Hopkinson杆技术分析典型传爆药的动态力学性能 [J]. 含能材料, 2012, 20(1): 62–66. DOI: 10.3969/j.issn.1006-9941.2012.01.015.ZHANG Z M, XU B Y, JIA J X, et al. Analysis on dynamic properties of typical boosters based on Hopkinson bars [J]. Chinese Journal of Energetic Materials, 2012, 20(1): 62–66. DOI: 10.3969/j.issn.1006-9941.2012.01.015. [13] GRAY Ⅲ G T, IDAR D J, BLUMENTHAL W R, et al. High- and low-strain rate compression properties of several energetic material composites as a function of strain rate and temperature [R]. Washington, DC, USA: Office of Scientific & Technical Information Technical Reports, 1998. [14] IDAR D J, THOMPSON D G, GRAY Ⅲ G T, et al. Influence of polymer molecular weight, temperature, and strain rate on the mechanical properties of PBX 9501 [J]. AIP Conference Proceedings, 2002, 620(1): 821–824. DOI: 10.1063/1.1483663. [15] RENGANATHAN K, RAO B N, JANA M K. Failure assessment on a strip biaxial tension specimen for a HTPB-based propellant material [J]. Propellants, Explosives, Pyrotechnics, 2015, 24(6): 349–352. DOI: 10.1002/(SICI)1521-4087(199912)24:6<349::AID-PREP349>3.0.CO;2-1. [16] OTTOSEN N S, RISTINMAA M. The mechanics of constitutive modeling [M]. Oxford UK: Elsevier Limited, 2005: 165-174. -