2. 清华大学航天航空学院工程力学系应用力学教育部重点实验室, 北京 100084
2. School of Aerospace, AML, Department of Engineering Mechanics. Tsinghua University, Beijing 100084, China
由于炸药宏观力学变形导致能量耗散所产生的温升很低, 无法使得其发生化学反应, 因此目前普遍认为热点的产生来自炸药内局部区域的热量积聚, 最终导致发生化学反应。炸药内部“热点”生成机理常见的有[1]:内孔的塌陷、摩擦生热、绝热剪切等。
围绕热点生成机理已开展了大量的实验和理论研究, 并建立了多种模型, 如粘塑性孔洞塌缩模型[2-3]、绝热剪切生热模型[4-6]、炸药内部颗粒及粘结剂摩擦生热模型[7]等。在“热点”生成机理中, Massoni[8]等认为粘塑性孔洞塌缩生热可产生较多的能量, 更容易产生热点。众多学者针对一维粘塑性孔洞塌缩模型进行了深入的研究, Kim[9]将粘塑性孔洞塌缩热点生成模型与低压下炸药点火模型相连, 研究了PBX-9404炸药冲击起爆过程。Whitworth[10]在Carroll-Holt[11]塑性孔洞模型基础上, 加入了孔洞粘-塑性运动过程产生的温升。梁增友[12]基于Kim粘塑性孔洞塌缩热点模型原理, 建立了三项式整体化学反应速率方程, 描述了PBX炸药在孔隙率、颗粒尺寸等损伤变化下的冲击起爆过程。张振宇[13]在粘塑性孔洞模型中考虑了熔化效应, 研究了多孔低熔点TNT炸药热点生成过程。温丽晶[14]在一维粘塑性孔洞塌缩模型基础上发展了双球壳模型, 用于描述炸药固相颗粒与粘结剂混合的细观结构, 给出了不同组分对热点产生的影响。
以上孔洞塌缩模型使用时, 均考虑孔洞受一直持续的外压作用, 而实际中孔洞受力非常复杂。侵彻战斗部装药抗高过载试验结果表明[15], 战斗部在撞击目标过程中内部装药将经受拉压交变等复杂的力学环境, 会产生裂纹、疏松等预损伤, 容易产生热点, 引起装药提前爆燃。为探索炸药在复杂冲击环境下热点生成机理, 本研究在一维粘塑性孔洞塌缩模型基础上, 考虑了孔洞在拉伸及压缩等复杂受力下的运动模式, 引入炸药力学性能随温度及损伤的影响, 建立了炸药孔洞塌缩热点生成模型, 分析了炸药在持续压缩及拉伸与压缩交变两种不同外力作用下热点的生成模式以及规律。
2 炸药孔洞塌缩热点生成模型 2.1 一维粘塑性孔洞塌缩模型参考Whitworth[10]一维粘塑性孔洞塌缩模型理论, 孔洞内径随时间变化的运动方程为:
$ -\rho \left( {a\ddot a + 1.5{{\dot a}^2}} \right) = {p_{\rm{s}}} + {p_{\rm{g}}} + {p_\nu }-{p_{{\rm{y}}}} $ | (1) |
式中, ρ为材料密度, g·cm3; a为孔洞内径, μm; ps为施加在孔洞的外力, pg为孔洞内气体压力, pv为粘性力, py为弹塑性应力, 其中的+ps表示孔洞承受压缩, -ps表示孔洞承受拉伸, MPa;
在一维粘塑性孔洞塌缩模型中, 孔洞运动过程可分为三个阶段:
(1) 弹性阶段(α3≥α≥α1)
$ {p_\nu } = 0, {p_{\rm{y}}} = \frac{{4\mu \left( {{\alpha _0}-\alpha } \right)}}{{3\alpha \left( {\alpha-1} \right)}}, {\alpha _1} = \frac{{2\mu {\alpha _0} + Y}}{{2\mu + Y}}, {\alpha _3} = \frac{{2\mu {\alpha _0}-Y}}{{2\mu - Y}} $ | (2) |
(2) 弹-粘塑性过渡阶段, 压缩状态时(α1≥α≥α2), 拉伸状态时(α4≥α≥α3)
$ {p_\nu } = 12\eta {a^2}\dot a\int_a^d {\frac{1}{{{r^4}}}{\rm{d}}r}, {\alpha _2} = \frac{{2\mu {\alpha _0}}}{{2\mu + Y}}, {\alpha _4} = \frac{{2\mu {\alpha _0}}}{{2\mu-Y}} $ | (3) |
$ {p_{\rm{y}}} = \frac{2}{3}Y\left( {1-\frac{{2\mu |\left( {{\alpha _0}-\alpha } \right)|}}{{Y\alpha }} + {\rm{ln}}\left\{ {\frac{{2\mu \left| {\left( {{\alpha _0}-\alpha } \right)} \right|}}{{Y\left( {\alpha - 1} \right)}}} \right\}} \right) $ | (4) |
$ d = \sqrt[3]{{\frac{{2\mu B}}{Y}}}, B = \frac{{\alpha _0^3\left| {\left( {{\alpha _0} -\alpha } \right)} \right|}}{{{\alpha _0} -1}} $ | (5) |
(3) 完全粘塑性阶段, 压缩状态时(α2≥α≥1), 拉伸状态时(α≥α4)
$ {p_\nu } = 12\eta {a^2}\dot a\int_a^b {\frac{{{\rm{d}}r}}{{{r^4}}}}, {p_{\rm{y}}} = \frac{2}{3}Y{\rm{ln}}\left( {\frac{\alpha }{{\alpha-1}}} \right) $ | (6) |
式中,μ为剪切模量, MPa; Y为屈服强度, MPa; η为粘性系数, MPa·μs; d为孔洞内弹性区域与粘塑性区域边界, μm。图 1为孔洞模型示意图, a0为孔洞初始内径, μm; b0为孔洞初始外径, μm; b为对应内径a时的孔洞外径, μm; α为孔隙率, 初始孔隙率α0= b03/ b03-a03, α1、α2、α3、α4分别为不同时刻孔隙率, %。
一般情况下, 炸药材料力学性能受到应变强化效应、应变率强化效应、温度软化效应、材料损伤等的影响[16-17], 具体见公式(7)。
$ \sigma = Y{\left( {1 + \frac{\varepsilon }{{{\varepsilon _y}}}} \right)^n}{\left( {1 + b\dot \varepsilon } \right)^m}h\left( {{T^*}} \right)\left( {1-D} \right) $ | (7) |
式中, σ为应力, MPa; ε为应变, εy为屈服应变(εy=Y/E); n为应变强化因子;
$ \sigma = \left( {{p_\nu }-{p_{\rm{y}}}} \right)h\left( {{T^*}} \right)\left( {1-D} \right) $ | (8) |
由于炸药材料在发生点火过程中温升幅度较大(300~800 K), 故孔洞周围炸药力学性能随温升的变化不可忽略, 假设温度升高对炸药热软化效应服从线性关系[16-17], 如公式(9), 其中c为热软化系数, 一旦温度上升到1/c, 材料完全丧失了承载能力。在孔洞运动过程中, 局部温升的增长主要源于内部塑性应变, 得到了h(T*) 随应变ε的变化关系(公式(10))。
$ h\left( {{T^*}} \right) = 1-c{T^*} $ | (9) |
$ h{(T^*)} = {\rm{exp}}\left\{ {-\frac{{\beta c}}{{\rho {C_v}}}\left( {{p_\nu }-{p_{\rm{y}}}} \right)\varepsilon } \right\} $ | (10) |
式中, β为塑性功转化率(一般取为0.95); Cv为热容系数, MPa·cm3·g-1·K-1; ε为应变。
假设外力作用下孔洞材料损伤服从weibull统计分布[16-17], 具体表达式为:
$ D\left( \varepsilon \right) = 1-{e^{-\varepsilon }} $ | (11) |
$ \varepsilon = 2{\rm{ln}}\left( {r/{r_0}} \right) $ | (12) |
式中, r为孔洞内某点当前的坐标,μm; r0为该点的初始坐标, μm。
利用公式(10)和(11)分别获得了炸药孔洞内不同位置r处温度热软化效应h(T*)以及损伤因子D随应变ε的变化关系, 通过对h(T*)和D(ε) 在整个孔洞内积分, 可获得材料的应力σ:
$ \sigma = \int_{{\varepsilon _y}}^\infty {\frac{{\left( {{p_\nu }-{p_{\rm{y}}}} \right)h\left( {{T^*}} \right)\left( {1-D} \right)}}{{{e^{\frac{3}{2}\varepsilon }}-1}}d\varepsilon } $ | (13) |
因此在考虑材料温度热软化和损伤效应情况下, 公式(1)可表示为:
$ -\rho \left( {a\ddot a + 1.5{{\dot a}^2}} \right) = {p_{\rm{s}}} + {p_{\rm{g}}} + \sigma $ | (14) |
孔洞内径塌缩运动而产生的温升[10]可表示为:
$ \rho {{C}_{v}}\dot{T}=12\eta \left( \frac{\text{ }{{u}^{2}}}{{{r}^{2}}} \right)+2Y\left( \frac{\left| \text{ }u \right|}{r} \right),u=\frac{\dot{a}{{a}^{2}}}{{{r}^{2}}} $ | (15) |
式中,
前面给出了复杂受力环境下, 考虑炸药材料受温度和损伤影响的孔洞塌缩温升过程。通常认为炸药温度超过点火温度Tc后将发生点火燃烧过程, 即孔洞内产生炸药反应气体压力pg。文中孔洞点火过程采用Ward、Son和Brewster[18]提出的WSB公式描述, 其中反应过程中的质量变化率(
$ \dot m({T_s}) = \left[{\frac{{{A_{\rm{c}}}T_s^2{e^{-{E_{\rm{c}}}/{T_s}}}}}{{{E_{\rm{c}}}\left( {{T_s}-{T_0}-{Q_{\rm{c}}}/2} \right)}}} \right]{^{\frac{1}{2}}} $ | (16) |
当炸药点火产生气体后, 孔洞内的气体压力变化服从绝热气体膨胀规律:
$ {p_{\rm{g}}} = {p_{{\rm{g}}0}}{\left( {{\rho _g}/{\rho _{{\rm{g0}}}}} \right)^\gamma } $ | (17) |
$ \frac{{{\rm{d}}{\rho _{{\rm{g}}}}}}{{{\rm{d}}t}} =-3\frac{{{\rho _{\rm{g}}}}}{a}\left( {\dot {a-\dot m}{\rho _{\rm{g}}}} \right) $ | (18) |
式中, 下标c代表炸药压缩阶段, 下标g代表气体阶段, Ac为反应率因子, s-1; Ec为激活能, kJ·mol-1; Qc为化学能, kJ·kg-1; T0为炸药固体温度, K; Ts为炸药表面温升, K; pg为气体压力, Pa; ρg为气体密度, kg·m-3。
3 复杂受力环境下热点生成分析选用PBX9404为研究对象, 具体材料参数见参考文献[10]。当炸药孔洞内径最高温升超过点火温度Tc=800 K时, 开始孔洞内点火过程计算, 采用WSB方程和国际单位制, 求得孔洞内气体密度和压强。WSB方程计算过程中相关参数见参考文献[18]。
3.1 压缩外力作用图 2给出ps随时间变化情况。图 3为压缩外力ps=500, 400, 300 MPa时孔洞内径a变化情况, 对比发现500 MPa和400 MPa均可导致热点产生, 300 MPa下孔洞内径塌缩较为缓慢, 未产生热点, 并且外界压力越高塌缩越快。当ps=500, 400 MPa时, 孔洞在外界压力ps以及内部气体压力pg的共同作用下运动, 当内部气体压力超过外界压力时, 孔洞内径开始向外快速扩展, 热点产生。图 4为不同压力下孔洞由于粘塑性运动产生的温升, 图 5为热软化及损伤因子(即h(T*)(1-D))变化情况。
图 3~图 5给出了在不同压缩外力下孔洞热点产生过程, 在外力作用下孔洞内径开始缩小, 应变不断增大, 应变导致了软化和损伤因子的增大, 降低了材料性能, 加剧了孔洞内径塌缩, 使得温升快速增加, 孔洞内部产生热点。当温升超过点火温度, 孔洞内不断产生化学反应气体, 随着气体压力的升高, 孔洞内径开始向外膨胀, 孔洞点火过程开始。
3.2 压缩与拉伸交变外力作用在战斗部侵彻目标过程中, 战斗部内装药经受复杂的力学环境, 最典型承受拉伸与压缩交变的外力作用, 图 6假设两种典型的交变外力作用, 图 6a首先是拉伸外力ps=-300 MPa作用接着压缩外力ps=+300 MPa作用(记为ps=-300 MPa+300 MPa), 图 6b首先是压缩外力+300 MPa作用接着拉伸外力-300 MPa作用(记为ps=+300 MPa-300 MPa), 图 6中每次拉伸或压缩外力作用持续时间假设为0.5 μs。
下面说明在压缩和拉伸交变作用下如何利用一维粘塑性孔洞塌缩模型进行计算。孔洞在第一个0.5 μs内先受到拉伸外力ps=-300 MPa的作用(图 6a), 内径由10 μm增大到11 μm(图 7), 温度升高到320 K(图 8), 孔洞塑性应变达到0.25(图 9), 图 10中孔洞热软化及损伤因子增大为0.1。在第二个0.5 μs内, 孔洞受到压缩外力的作用, 此时初始内径a、热软化及损伤综合因子取第一个0.5 μs结束时的11 μm和0.1, 孔洞内径在压缩外力作用下减小至10 μm, 温度升高到350 K, 塑性应变为0.5, 热软化及损伤综合因子增大为0.18。在第三个0.5 μs内, 孔洞再次受到拉伸外力的作用, 孔洞初始内径取第二个0.5 μs结束时的10 μm, 由于损伤的累积作用, 孔洞内径在第三循环结束时增大至13 μm, 温度升高到400 K。以此类推, 在第七个循环即3.2 μs时孔洞内径温升超过点火温度800 K, 孔洞内产生气体, 使得孔洞内径急剧增大, 产生热点。对比图 7~图 10还可以看出, 图 6a和图 6b两种不同压缩与拉伸交变作用模式, 首先受拉伸作用比压缩作用更容易产生热点。
图 11分别为拉伸与压缩交变外力ps=300, 250, 200 MPa条件下孔洞的运动过程, 可以看出300 MPa与250 MPa条件下均可产生热点, 但200 MPa条件下未能产生热点。图 12为在拉伸与压缩交变外力ps=+300 MPa-300 MPa作用下以及持续的压缩ps=+300 MPa作用下孔洞运动过程, 可见在拉伸和压缩交变作用比持续压缩作用更容易产生热点。
对比3.1和3.2节计算结果, 可见在拉伸与压缩交互作用下, 拉伸作用使得孔洞内发生塑性变形, 增加了损伤因子, 同时也增大孔洞孔隙率α, 最终导致孔洞塌缩加速并产生高温。因此在战斗部装药安全性研究中, 拉伸与压缩交变外力作用使得孔洞更容易产生热点。该模型可为装药安全性研究提供技术手段, 后续将模型拓展应用到战斗部侵彻目标等复杂的力学过程中。
4 结论建立了复杂冲击环境下孔洞塌缩热点生成模型, 用该模型分析了非均质炸药PBX-9404孔洞分别在压缩、拉伸与压缩交变外力作用下, 热点产生过程及规律, 得出如下结论:
(1) 新模型考虑了孔洞在拉伸、压缩下的不同运动规律, 以及材料性能随温度软化及损伤影响, 故孔洞塌缩模型适用于模拟复杂外载受力环境下热点产生。
(2) 单一持续压缩作用, 400 MPa条件下孔洞可产生热点, 而300 MPa不能产生热点; 拉伸与压缩交变作用, 250 MPa条件下孔洞可产生热点。
(3) 对比单一持续压缩、拉伸与压缩交变作用不同受力状态, 拉伸与压缩交变作用比单一压缩更容易产生热点。
[1] |
Field J E, Bourne N K, Palmer S J P, et al. Hot-spot ignition mechanisms for explosives and propellants[J].
Philos Trans R Soc, 1992, 339(1654): 269-283. DOI:10.1098/rsta.1992.0034 |
[2] |
Kang J, Butler P B, Baer M R. A thermomechanical analysis of hot spot formation in condensed-phase energetic materials[J].
Combust and Flame, 1992, 89(2): 117-139. DOI:10.1016/0010-2180(92)90023-I |
[3] |
Khasainov B A, Borisov A A, Ermolaev B S, et al. Two-phase visco-plastic model of shock initiation of detonation in high density pressed explosives[C]//In Proceedings of the 7th Symposium International on Detonation, Annapolis, MD, 1981: 435-447.
|
[4] |
Frey R B. The initiation of explosive charges by rapid shear[C]// In Proceedings of the 7th Symposium International on Detonation, Annapolis, MD, 1981: 36-42.
|
[5] |
Grady D E, Kipp M E. The growth of inhomogeneous thermoplastic shear[C]//In Proceedings of the International Conference on Mechanical and Physical Behavior of Materials Under Dynamics Loading, Paris, 1985.
|
[6] |
Kipp M E. Modeling granular explosives detonation with shear band concepts[C]//In Proceedings of the 8th Symposium International on Detonation, Albuquerque, NM, 1985: 35-41.
|
[7] |
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: 2303-2322. DOI:10.1016/S0022-5096(98)00011-8 |
[8] |
Massoni J, Saurel R. A mechanistic model for shock initiation of solid explosives[J].
Physics of Fluids, 1999, 11(3): 710-735. DOI:10.1063/1.869941 |
[9] |
Kim K. Development of a model of reactive rates in shocked composite explosives[C]//In Proceedings of the 9th Symposium International on Detonation, 1989: 593-603.
|
[10] |
Whitworth N J. Mathematical and numerical modeling of shock initiation in heterogeneous solid explosives[D]. Cranfield University, 2008.
|
[11] |
Carroll M M, Holt A C. Static and dynamic pore-collapse relations for ductile porous materials[J].
Journal of Applied Physics, 1972, 43(4): 1626-1635. DOI:10.1063/1.1661372 |
[12] |
梁增友, 黄风雷, 张振宇. 损伤炸药的冲击起爆数值模拟[J].
北京理工大学学报, 2006, 26(12): 1047-1051. LIANG Zeng-you, HUANG Feng-lei, ZHANG Zhen-yu. Numerical simulation of damaged explosive in shock detonation[J]. Transactions of Beijing Institute of Technology, 2006, 26(12): 1047-1051. DOI:10.3969/j.issn.1001-0645.2006.12.004 |
[13] |
张振宇, 浣石, 卢芳云, 等. 多孔TNT炸药中热点形成的粘塑性塌缩机理[J].
含能材料, 1994, 2(2): 36-42. ZHANG Zhen-yu, HUAN Shi, LU Fang-yun, et al. Viscoplastic collapse mechanism of hot-spot formation in porous TNT explosives[J]. Journal of Energetic Materials(Hanneng Cailiao), 1994, 2(2): 36-42. |
[14] |
温丽晶, 段卓平, 张振宇, 等. 弹粘塑性双球壳塌缩热点反应模型[J].
高压物理学报, 2011, 25(6): 493-500. WEN Li-jing, DUAN Zhuo-ping, ZHANG Zhen-yu. An elastic/viscoplastic pore collapse model of double layered hollow sphere for hot-spot ignition in shocked explosives[J]. Chinese Journal of High Pressure Physics, 2011, 25(6): 493-500. DOI:10.11858/gywlxb.2011.06.003 |
[15] |
高金霞, 赵卫刚, 郑腾. 侵彻战斗部装药抗过载技术研究[J].
火工品, 2008, 5(4): 4-7. GAO Jin-xia, ZHAO Wei-gang, ZHEN Teng. Study on the antioverloading technique for penetrating warhead charge[J]. Initiators and Pyrotechnic, 2008, 5(4): 4-7. |
[16] |
Wei Z G, Batra R C. Damage model for anisotropic materials, and its application to analysis of stability and spallation[J].
International Journal of Impact Engineering, 2007, 34(11): 1780-1796. DOI:10.1016/j.ijimpeng.2006.10.002 |
[17] |
Wu X Y, Ramesh K T, Wright T W. The effects of thermal softening and heat conduction on the dynamic growth of voids[J].
International Journal of Solids and Structure, 2003, 40(17): 4461-4478. DOI:10.1016/S0020-7683(03)00214-2 |
[18] |
Ward M J, Son S F, Brewster M Q. Steady deflagration of HMX with simple kinetics: a gas phase chain reaction model[J].
Combustion and Flame, 1998, 114(3): 556-568. |
The new “hot-spot”model under the complex impact environment was established based on the one-dimensional elastic-viscoplastic collapse model. The effects of temperature and damage on the mechanical properties of the explosives were taken into account in the formation of hot-spot.