Loading [MathJax]/jax/element/mml/optable/BasicLatin.js

不同类型装药侵彻安全性数值模拟

白晨 杨昆 吴艳青 高洪泉 薛海蛟

白晨, 杨昆, 吴艳青, 高洪泉, 薛海蛟. 不同类型装药侵彻安全性数值模拟[J]. 高压物理学报, 2021, 35(6): 065101. doi: 10.11858/gywlxb.20210754
引用本文: 白晨, 杨昆, 吴艳青, 高洪泉, 薛海蛟. 不同类型装药侵彻安全性数值模拟[J]. 高压物理学报, 2021, 35(6): 065101. doi: 10.11858/gywlxb.20210754
BAI Chen, YANG Kun, WU Yanqing, GAO Hongquan, XUE Haijiao. Numerical Simulation for PBX Charges Safety of Different Types During Penetration[J]. Chinese Journal of High Pressure Physics, 2021, 35(6): 065101. doi: 10.11858/gywlxb.20210754
Citation: BAI Chen, YANG Kun, WU Yanqing, GAO Hongquan, XUE Haijiao. Numerical Simulation for PBX Charges Safety of Different Types During Penetration[J]. Chinese Journal of High Pressure Physics, 2021, 35(6): 065101. doi: 10.11858/gywlxb.20210754

不同类型装药侵彻安全性数值模拟

doi: 10.11858/gywlxb.20210754
基金项目: 国家自然科学基金(11872119);博士后科学基金(BX20200046,2020M680394)
详细信息
    作者简介:

    白 晨(1997-),男,硕士研究生,主要从事弹药安全性研究. E-mail:906587639@qq.com

    通讯作者:

    吴艳青(1974-),女,博士,教授,主要从事高能炸药细观力学研究. E-mail:wuyqing@bit.edu.cn

  • 中图分类号: O341; TJ55

Numerical Simulation for PBX Charges Safety of Different Types During Penetration

  • 摘要: 为研究侵彻过程典型装药的力学响应和损伤点火过程,采用高聚物黏结炸药微裂纹-微孔洞力热化学耦合细观模型,对侵彻过程中PBX装药的力学响应、应力波传播情况和损伤-温升机理进行了研究,标定了两类炸药的本构模型参数,对比分析了压装PBX04和浇注GOFL-5两类装药的力学-损伤-温升响应的差异性。结果表明:GOFL-5炸药的屈服强度、硬化模量、初始微裂纹密度和微裂纹尺寸均小于PBX04炸药;加载初期,压装药头部微裂纹损伤高于浇注药;整个侵彻过程中,两类炸药的微裂纹损伤较严重区域均为装药头部和尾部;剪切裂纹热点为压装药主导的温升机制,且浇注药GOFL-5在侵彻过程中的温升较压装药PBX04更低。

     

  • 侵彻战斗部主装药广泛采用PBX炸药,由于服役环境可能同时存在多种外界刺激,炸药自身的力学响应会影响侵彻过程中装药的安全性和战斗部的毁伤能力[1],因此PBX装药对弹药的整体性能有着决定性意义。而在目标侵彻过程中,侵彻战斗部的内部装药面临复杂的力学环境[2],同时复杂的装药结构也会对内部装药产生影响。在这些外部刺激下,炸药中的局部损伤处可能产生热点,导致意外点火的发生,进而可能转化成为更剧烈的爆燃或者爆轰,引发多种事故。因此研究各种复杂侵彻环境下装药的力学性能与损伤响应,对保障弹药在整个服役过程中的安全性与可靠性具有重要意义。

    对于侵彻过程中内部装药的安全性而言,建立能够准确描述炸药力学响应和损伤行为的本构模型,是通过数值模拟方法研究侵彻装药安全性问题的关键。PBX炸药动态损伤本构模型可以划分为两类:第一类是基于黏弹性、黏塑性等连续介质损伤力学基础理论,考虑应变率、温度对PBX炸药力学行为的影响而建立的宏观唯象经验本构模型[3-9];第二类则考虑微缺陷是PBX炸药中一类重要的细观结构,是基于材料力学损伤-点火行为和细观缺陷演化与热点形成之间的联系而建立的PBX炸药细观损伤力学本构模型[10-13]。近年来,国内外学者基于以上模型对侵彻装药安全性进行了广泛的研究。张馨予等[14]将孔隙压塌损伤、炸药晶体破碎损伤、黏结剂脱黏等多种细观损伤形式耦合到炸药宏观本构模型中,研究了侵彻环境下弹体装药的损伤分布情况。石啸海等[15-16]基于内聚裂纹模型,模拟了战斗部侵彻半无限大混凝土过程中PBX装药的动态力学响应及损伤演化,并且对某种缩比弹侵彻混凝土靶板进行了数值模拟,从过载、裂纹宽度、裂纹含量等角度比较了弹头形状对装药损伤的影响,结果表明控制弹头曲径比有利于减小装药损伤。成丽蓉等[17]基于裂纹摩擦、孔洞塌缩两种热点生成机制细观模型,开展了侵彻单层和多层典型靶板时战斗部装药的动态响应、损伤演化及热点生成对比研究。Li等[18]基于微裂纹动态损伤模型,通过数值模拟方法预测了实弹侵彻实验中PBX1314的损伤分布、点火位置与点火区域的反应进程,结果表明,侵彻过程中经多次撞击载荷作用,PBX1314尾端面局部区域达到点火临界条件而发生点火,反复撞击过程中PBX1314内部微裂纹摩擦引起的能量局部化是引发点火的重要原因。

    上述模型往往仅考虑了剪切裂纹热点等单一缺陷机制的细观力热响应过程,在研究复杂载荷条件下不同类型PBX炸药损伤-热点主导机制的自适应能力方面,模型还需要进一步改进。此外,上述研究工作多围绕压装型炸药PBX展开,针对浇注型PBX炸药的研究较少,由于浇注类PBX炸药和压装类PBX炸药的材料组分、物理状态和成型工艺不同,力学性能与点火特性也存在差异,因此对比其在侵彻环境下的响应特性,对于战斗部装药材料的选取和装药结构设计具有重要意义。基于此,本研究应用前期发展的PBX炸药微裂纹-微孔洞力热化学耦合细观模型[18],考虑微裂纹-微孔洞两种细观缺陷演化对炸药损伤-热点形成的影响,通过动态分离式霍普金森压杆(SHPB)实验校核模型参数,并利用数值模拟方法分析两类典型装药(压装PBX04和浇注GOFL-5)在弹体侵彻混凝土薄板过程中的应力波传播、损伤演化和温升响应情况,为深入理解侵彻过程中装药的力学-损伤-点火响应提供参考。

    PBX炸药微裂纹-微孔洞力热化学耦合细观模型(CMM)中考虑了拉伸张开、剪切张开、纯剪切、剪切摩擦、摩擦自锁5种微裂纹演化模式,以及微孔洞坍塌与扭曲变形两种演化模式,如图1所示。同时,CMM模型还囊括了剪切裂纹热点与孔洞坍塌热点子模型,具有复杂应力状态-微缺陷演化模式自判断能力,以及两种微缺陷热点自启动能力。CMM模型中总体应力、应变分别分解为偏量部分与体量部分,依次建立微裂纹相关偏量本构关系与微孔洞相关体量本构关系,二者通过Gurson屈服准则进行耦合,通过微裂纹、微孔洞演化方程更新相关变量,建立剪切裂纹摩擦热点与孔洞坍塌热点子模型,模型细节与算法实现详见文献[1920]。

    图  1  CMM模型所考虑的微缺陷演化机制与变形机制示意图
    Figure  1.  Conceptual diagram of all kinds of considered mechanisms in the current model

    将总体偏应变(ε)分解为黏弹性应变(εve)、塑性应变(εp)以及由于微裂纹扩展所引起的微裂纹应变(εcr) 3部分,分别描述PBX炸药内黏结剂基体材料所表现出的黏弹性变形、塑性变形以及由于微裂纹扩展对材料力学行为的影响。

    PBX炸药黏弹性变形由广义Maxwell模型描述,微裂纹张开/剪切扩展引起的裂纹应变由SCRAM模型描述,材料的偏量本构关系可表示为

    ˙S=2GA0(˙ε˙εp)B0(S+C0) (1)

    式中:S为应力偏量,G为剪切模量,A0B0C0表示为

    A0=11+αe3(ˉc/a)3,B0=αe(ˉc/a)2˙ˉc/a1+αe3(ˉc/a)3,C0=Nn=1Sn/τnαe(ˉc/a)2˙ˉc/a (2)

    式中:ˉc为微裂纹平均尺寸;a为微裂纹初始特征尺寸,a3=6GβSnτn为广义黏弹性体元中第n个Maxwell黏弹性单元中的偏应力和松弛时间。

    αe={3p05νp<0,β=64π(1ν)N015(2ν)G (3)

    式中:p为压力,N0为初始裂纹密度,ν为泊松比。

    为描述材料内微孔洞演化所发生的不可逆损伤对PBX炸药体积变形的影响,采用孔隙率相关状态方程

    p(ρ,e,f)=(1f)[ρsc20ηs(1sηs)2(1Γsηs2)+Γsρses] (4)

    式中:f为孔隙率;ρ、e分别为孔隙材料的密度与比内能,密实材料对应的ρs = ρ/(1 − f),es = eηs = 1 − ρs0/ρsρs0为密实材料的初始密度;Γs为Grüneisen系数;c0s为材料参数。

    含孔隙PBX材料的塑性变形采用经典Gurson模型描述,模型中材料屈服面与von-Mises等效应力σe和压力p相关

    F(σe,p,f)=(σeYM)2+2fcosh(3p2YM)f21=0 (5)

    考虑动态加载下材料硬化效应与应变率效应,密实材料(f = 0)的屈服强度可表示为

    YM=[σ0+h(ˉεpM)n][1+Cln(1+˙ε)] (6)

    式中:σ0˙ε0=10 s−1加载下所对应的初始屈服应力,ˉεpM为等效塑性应变,˙ε=˙ε/˙ε0为等效应变率,h为硬化模量,C为应变率相关系数。

    基于Griffith能量释放率裂纹扩展准则,微裂纹扩展方程可表示为

    ˙ˉc=˙cmax (7)

    式中:{\dot c_{\max }}为最大扩展速度,\bar \gamma 为材料比表面能,{g_{\rm{dom}}}\left( {{{\sigma}},\bar c} \right)为主裂纹对应的能量释放率。

    由于材料内微裂纹方向分布具有随机性,存在临界微裂纹方向,其对应的能量释放率最大,即该方向的微裂纹在最小施加应力下最先发生失稳扩展,因此定义为主裂纹。主裂纹方向的确定与当前应力状态相关(\sigma 1\sigma 3平面),微裂纹扩展对应的材料整体损伤度定义为{d_{\rm{cr}}} = {{{{\bar c}^3}} / {\left( {{a^3} + {{\bar c}^3}} \right)}}

    考虑孔洞坍塌引起的孔隙率减小,以及孔洞扭曲引起的孔隙率增加两种变形机制,孔隙率演化方程可表示为

    \dot f = \left( {1 - f} \right)\dot \varepsilon _V^{\rm{p}} + f{k_{\rm{w}}}\omega \left( \sigma \right)\frac{{{s_{ij}}\dot e_{ij}^{\rm{p}}}}{{{\sigma _{\rm{e}}}}} (8)

    式中:\dot \varepsilon _V^{\rm{p}} = {{\dot \varepsilon _{\rm{kk}}^{\rm{p}}} / 3},为塑性体积应变率;kw为剪切相关材料参数;\omega 与应力状态相关,取值范围为0 ≤ \omega ≤ 1。

    \omega \left( {{\sigma}} \right) = 1 - {\left( {\frac{{27{{{J}}_3}}}{{2{{\sigma}} _{\rm{e}}^3}}} \right)^2} (9)

    式中:J3为表示应力张量的第三不变量。

    采用一维热传导方程来描述剪切裂纹表面及其周围区域摩擦生热、熔化、点火与传热等热力学过程

    \rho {c_V}{\dot T_{\rm{hs}}} = \frac{\partial }{{\partial x}}\left( {\kappa \frac{{\partial {T_{\rm{hs}}}}}{{\partial x}}} \right) + \rho {Q_{\rm{r}}}Z{{\rm{e}}^{ - E/(R{T_{\rm{hs}}})}} + \varphi {\mu _v}{\dot \varepsilon_{\rm{m}} ^2} (10)

    式中:x为沿微裂纹法向的坐标轴;等式右端3项分别表示热传导项、化学反应释放热量以及熔化区域(ThsTm)液相黏性流动生热;ThscVκQr分别为微裂纹热点温度、比定容热容、热传导系数以及单位质量化学反应放热;ZE为Arrhenius反应速率方程参数;\varphi (0 ≤ \varphi ≤ 1)为熔化百分数;{\;\mu _v}为熔化区域黏性系数;\dot \varepsilon_{\rm{m}}= \varphi {v_{\rm{c}}}/l为熔化区域的剪切应变率,vc为微裂纹滑动速度,熔化区域宽度l由熔化区域的传播速度确定,\dot l(t) = \Delta \left( {k{T_x}} \right)/\rho L,其中\Delta \left( {k{T_x}} \right)表示固液分界面热流量梯度。

    微孔洞周围材料温度(Tvo)分布情况由一维球形热传导方程进行描述

    \rho {c_V}{\dot T_{\rm{vo}}} = k\left[ {\frac{{{\partial ^2}{T_{\rm{vo}}}}}{{\partial {r^2}}} + \frac{2}{r}\frac{{\partial {T_{\rm{vo}}}}}{{\partial r}}} \right] + \dot w_{\rm{vp}}^* + \rho {Q_{\rm{r}}}Z{{\rm{e}}^{ - E/(R{T_{\rm{vo}}})}} (11)

    其中,孔洞周围单位体黏塑性功生成速率可表示为

    \dot w_{\rm{vp}}^*\left( r \right) = s_{ij}^*\dot \varepsilon _{ij}^* = \frac{{ - 2{Y_{\rm{M}}}\dot \varepsilon _V^{\rm{p}}}}{{{r^3}}}\frac{{1 - {f_0}}}{{1 - {f_{\rm{vc}}}}}b_0^3 + \frac{{4\eta {{\left( {\dot \varepsilon _V^{\rm{p}}} \right)}^2}}}{{{r^6}}}{\left( {\frac{{1 - {f_0}}}{{1 - {f_{\rm{vc}}}}}} \right)^2}b_0^6 (12)

    式中:等号右边两项分别表示基体材料的塑性效应与黏性效应所产生的功率,b0为微孔洞初始外径,\eta 为微孔洞周围基体材料相关的黏性系数。

    应用SHPB方法对中国工程物理研究院化工材料研究所提供的压装和浇注两类PBX炸药进行动态力学性能测试,浇注类复合炸药选取GOFL-5炸药(HMX、FOX-7、黏结剂的质量分数分别为50%、35%、15%),样品尺寸为20 mm × 20 mm,密度为1.75 g/cm3,初始孔隙率为1.01%。压装类复合炸药选取PBX04炸药(HMX、黏结剂的质量分数分别为95%、5%),密度为1.82 g/cm3,初始孔隙率为1.08%。在中国工程物理研究院化工材料研究所物理与力学性能实验室开展SHPB实验,实验中采用直径为20 mm的铝杆,子弹长度为300 mm,入射杆、透射杆长度分别为2 000、1 500 mm。Yang等[19]在前期研究了CMM模型中微缺陷参数({\bar c_0}, N0, \bar \gamma , f0, kw)对材料力学性能的影响规律。基于此,依次对CMM模型中的广义Maxwell黏弹性参数(G1G5, τ1τ5)、塑性参数({\sigma _0}, h, C, n)、微裂纹模型参数({\bar c_0}, N0, \bar \gamma, m)、微孔洞模型参数(f0, kw)以及Mie-Grüneisen状态方程参数(c0, s, Г)进行优化调整,使得标定的参数能够最大程度地符合计算与实验应力-应变曲线。压装炸药PBX04与浇注炸药GOFL-5标定曲线如图2所示,相关材料参数如表1所示,其中 {\mu _{\rm{s}}} 为静摩擦系数。由表1可知:两种炸药材料的黏弹性模量G1G5存在数量级差异,GOFL-5炸药的塑性参数({\sigma _0}, h, C)均远小于PBX-5炸药,GOFL-5炸药的初始微裂纹密度与微裂纹尺寸均小于PBX-5炸药。

    图  2  计算与实验得到的两类炸药在不同应变率下的应力-应变曲线
    Figure  2.  Calculated and experimental stress-strain curves of two kinds of explosives at different strain rates
    表  1  GOFL-5与PBX04材料参数
    Table  1.  Material parameters for GOFL-5 and PBX04
    Material\;\rho 0/(kg·m−3)G/GPa\nuG1/MPaG2/MPaG3/MPaG4/MPaG5/MPa{\tau{_1^{-1} } } /{\rm{s} }{^{-1} }
    GOFL-51 7500.550.3 167 30.45 90.03 185.6 120.00
    PBX041 8208.250.31 9401 175.001 521.001 909.01 688.00
    Material{\tau{_2^{-1} } } /{\rm{s} }{^{-1} }{\tau{_3^{-1} } } /{\rm{s} }{^{-1} }{\tau{_4^{-1} } } /{\rm{s} }{^{-1} }{\tau{_5^{-1} } } /{\rm{s} }{^{-1} }\sigma{{_0}}/MPaCh/MPan{\bar c}0/μm
    GOFL-57.32 × 1037.32 × 1047.32 × 1057.32 × 106 2.20.76 4.50.4530
    PBX049.00 × 1039.00 × 1049.00 × 1052.00 × 10640.00.101500.01.0030
    MaterialN0/cm−3\bar \gamma /(J·m−2){\dot c_{\max }}/(m·s–1)\;\mu_{\rm{s}}mf0/(m·s–1)kwc0/(m·s–1)s \varGamma
    GOFL-5 30.53000.35.00.012.01 0000.46 0.89
    PBX043001.43000.55.00.012.02 5002.26 1.50
    下载: 导出CSV 
    | 显示表格

    对比压装炸药与浇注炸药的应力-应变曲线,可以看到PBX04炸药峰值应力远高于GOFL-5炸药,但其破坏应变远小于GOFL-5炸药。两种材料的力学行为差异主要与其内部微缺陷数量和黏结剂含量相关。压装类PBX04炸药内微裂纹含量较多,且黏结剂含量较少,因此材料表现为准脆性材料破坏特征。浇注类GOFL-5炸药内部微缺陷数量较少,且对材料起到增韧效应的黏结剂含量较多,表现为韧性材料破坏特征。

    建立二维侵彻混凝土靶板的计算模型,设置轴对称边界条件和加载条件,建立1/2模型进行计算,选取侵彻速度为800 m/s,弹体尺寸如图3所示,侵彻的混凝土靶板厚度为0.5 m。计算过程中,由于侵彻过程中靶体内温度升高的区域是有限的,靶体外边界的温度不会受到侵彻过程中靶体温度变化区域的影响,因此在靶体边界处设置温度边界条件,边界温度设置成294 K。

    图  3  有限元计算模型
    Figure  3.  Finite element calculation model

    弹体材料本构模型使用Johnson-Cook模型。该模型是一个经验型的黏塑性本构模型,由Johnson和Cook在1983年首先提出,多用于描述金属材料在高载荷、高应变率和高温下的应力-应变关系,模型形式简单,本构参数少,能较好地描述材料的加工硬化效应、应变率效应和温度软化效应。在侵彻速度不高时弹体变形较小,因此本研究运用Johnson-Cook模型来描述壳体材料4340钢的力学响应,其主要材料参数见表2,其中:\rho为密度,c为比热容,α为热膨胀系数,Tm为熔化温度,\kappa 为导热系数,M为温度软化系数,N为应变强化参数,C1C4D1D5为材料常数。混凝土靶板的材料模型使用HJC本构模型[21]。该模型是针对混凝土材料提出的一种率相关损伤型本构模型,适用于大应变、高应变率和高压情况,其主要材料参数见表3,其中:ABN为材料硬化参数,FC为材料准静态单轴抗压强度,K1K3为材料体积模量参数,pCpLUCUL为压力参数,D1D2为损伤参数,\varepsilon_{\rm{min}}为最小断裂应变。

    表  2  弹体材料参数
    Table  2.  Parameters of projectile material
    Physical properties Johnson-Cook model
    \;\rho /(g·cm−3)c/(J·kg−1·K−1)\kappa /(kW·m−1·K−1)\alpha /(m·K−1)Tm/℃G/GPa NM
    7.8247838.113.24 × 10−51 793.15 774.970.26 1.03
    下载: 导出CSV 
    | 显示表格
    Johnson-Cook model Damage model
    C1/MPaC2/MPa C3 C4 D1 D2 D3 D4 D5
    792.21509.521.4 0 –0.8 2.1 –0.5 20.0 0.61
    下载: 导出CSV 
    | 显示表格
    表  3  混凝土靶板材料参数
    Table  3.  Material parameters of concrete plate target
    Physical properties HJC model
    \;\rho /(g·cm−3)c/(J·kg−1·K−1)κ/(kW·m−1·K−1)α/(m·K−1)G/GPaFC/MPaABNC
    2.286541.764.32 × 10−5 16.4040.680.751.650.767.0 × 10−3
    下载: 导出CSV 
    | 显示表格
    HJC model Damage model
    pC/MPaUC pL/GPa UL K1/GPa K2/GPa K3/GPa D1 D2 \varepsilon min
    13.565.80 × 10−4 1.05 0.10 17.40 38.80 29.80 0.03 1.0 0.01
    下载: 导出CSV 
    | 显示表格

    为分析药柱内压缩波传播情况,初始时刻沿药柱中心线由装药头部至尾部等间隔依次选取1#~5#位置,如图4(b)所示,图4(a)图5(a)分别为压装(PBX04)与浇注炸药(GOFL-5)内1#~5#位置处的压力变化曲线。由图4(a)可知,弹体与混凝土靶板初始接触时,弹体受到强烈的压缩波作用,装药头部压力最大,入射压缩波S1分别于104和179 μs到达装药中部(3#位置)和装药尾部(5#位置),当应力波传播至尾部自由端面时反射拉伸波,两种应力波相互抵消,装药尾部压力值最低。930 μs时,装药尾部在弹性恢复与拉伸波的作用下与弹体内壁面发生撞击,尾部压力迅速升高。由图5(a)可知,入射压缩波S1到达装药中部和装药尾部的时刻分别为123和213 μs,对比图4(a)图5(a)可知,相同侵彻速度下,浇注药内压缩波S1的波速(约2.11 km/s)小于压装药(约2.57 km/s)。

    图  4  压装药PBX04的压力演化情况
    Figure  4.  Pressure evolution of PBX04
    图  5  浇注药GOFL-5的压力演化情况
    Figure  5.  Pressure evolution of GOFL-5

    不同时刻压装、浇注药内压力云图演化如图4(b)图5(b)所示。随着压缩波的传播,药柱整体受压区域逐渐扩大,装药尾部与弹体内壁产生相对间隙,浇注GOFL-5炸药由于材料强度较低,抵抗变形的能力较弱,侵彻过程中装药尾部的最大间隙(dmax)较压装药更大(GOFL-5,dmax = 1.61 cm;PBX04,dmax = 0.68 cm)。结合图4(a)图5(a)可知,浇注GOFL-5炸药在侵彻过程中产生了较大的变形,装药尾部与弹体内壁发生撞击时尾部压力峰值(pmax)也更高(GOFL-5,pmax = 234 MPa;PBX04,pmax = 128 MPa)。

    图6为压装、浇注药柱内微裂纹损伤(dcr)演化情况。由图6(a)可知,加载初期压装药微裂纹损伤主要集中在头部,随着压缩波的传播,微裂纹损伤逐渐扩展至装药中部,装药尾部与弹体内壁发生撞击后尾部微裂纹损伤较严重。由图6(b)可知,浇注药的流动性较好,在侵彻过程中产生了较大的变形,尾部装药受到向药柱中心的拉应力作用产生裂纹扩展损伤。整个侵彻过程中,两类炸药微裂纹损伤较严重区域均为装药头部和尾部,由于开始侵彻时压装药内入射压缩波强度高于浇注药,同时压装药初始微裂纹密度较高,加载初期压装药头部微裂纹损伤高于浇注药。

    图  6  不同时刻两类炸药微裂纹损伤演化情况
    Figure  6.  Evolution of microcrack damage of two kinds of explosives

    图7为不同时刻压装、浇注药柱内不同位置处微孔洞损伤(dvo = f0ft)随时间演化曲线。对比图6图7可知,微孔洞坍塌时间尺度相对于微裂纹扩展时间尺度明显更小,微孔洞随着入射压缩波传播而发生坍塌,孔洞坍塌损伤发生的时刻与入射压缩波到达的时刻基本一致,而由于裂纹扩展损伤主要受拉伸和剪切状态影响,入射压缩波刚到达装药对应位置时并未发生明显的微裂纹扩展损伤。

    图  7  不同时刻两类炸药的微孔洞损伤演化
    Figure  7.  Evolution of microvoid damage of two kinds of explosives

    PBX04和GOFL-5内不同位置微裂纹与微孔洞热点温度随时间演化曲线分别如图8图9所示。由图8(a)可知,压装药微裂纹热点温升主要集中在装药头部和装药前端,1#和2#处的微裂纹热点峰值分别为362和487 K。根据CMM模型中由应力二轴度可确定5种微裂纹状态(拉伸张开、剪切张开、纯剪切、剪切摩擦、摩擦自锁),5种微裂纹演化模式如图10所示,由此可进一步分析温升较高位置处的微裂纹状态。整个侵彻过程中,1#和2#处5种微裂纹状态的频率分布如图11所示。侵彻过程中弹体内装药处于压力主导的应力状态,微裂纹大多处于摩擦自锁状态而非剪切裂纹扩展状态,同时热点密度较小可能引起热点湮灭现象,因此微裂纹热点温度达到峰值后逐渐降低。对比图8(a)图9(a)可知,由于GOFL-5材料的初始微裂纹密度、微裂纹尺寸以及微裂纹扩展速率较低,因此GOFL-5材料中剪切裂纹引起的热点温升较PBX04更低。结合图8(b)图9(b),随着入射压缩波的传播孔洞发生坍塌,与GOFL-5不同的是,PBX04装药尾部微孔洞温升峰值主要来自弹体撞击造成的孔洞坍塌,而坍塌孔洞周围黏塑性功引起的温升较小(约20 K),不足以引起点火。对比图8图9可知:对于压装药而言,裂纹摩擦相比孔洞坍塌引起的温升更高,剪切裂纹热点为压装药主导的温升机制;对于浇注药而言,两种热点机制引起的温升差别不大。

    图  8  不同时刻PBX04中微裂纹与微孔洞相关热点温度曲线
    Figure  8.  Microcrack and microvoid related hotspot temperature evolution curves for PBX04
    图  9  不同时刻GOFL-5中微裂纹与微孔洞相关热点温度曲线
    Figure  9.  Microcrack and microvoid related hotspot temperature evolution curves for GOFL-5
    图  10  微裂纹不同扩展形态示意图
    Figure  10.  Schematic diagram of different propagation patterns of microcracks
    图  11  PBX04装药中不同位置的微裂纹状态频率分布
    Figure  11.  Frequency distribution of microcracks state in different positions of PBX04 charge

    应用PBX炸药微裂纹-微孔洞力热化学耦合细观模型,研究了压装PBX04和浇注GOFL-5两类典型装药在弹体侵彻混凝土薄板过程中应力波的传播、损伤演化和温升响应情况,对比分析了相同侵彻条件下两类炸药力学-损伤-温升响应的差异性,得到以下主要结论。

    (1) 根据实验曲线标定了PBX04和GOFL-5微裂纹-微孔洞本构模型参数,两种炸药材料的弹性模量存在数量级差异,GOFL-5炸药的屈服强度、硬化模量、初始微裂纹密度和微裂纹尺寸均小于PBX04炸药。

    (2) 加载初期压装药头部微裂纹损伤高于浇注药,而浇注药的流动性较好,侵彻过程中产生了较大的变形,当装药尾部和壳体内表面发生撞击时形成高压区,整个侵彻过程中两类炸药微裂纹损伤较严重的区域均为装药头部和尾部,装药损伤较严重的区域往往容易引起能量聚集,进而在这些局部高温区形成热点,在侵彻弹体设计时应作为重点防护区域。

    (3) 通过计算可知,800 m/s速度侵彻混凝土薄板条件下,两种装药材料均未发生点火,裂纹摩擦热点为压装药主导的温升机制,而对于浇注药而言,两种热点机制引起的温升差别不大,且浇注药GOFL-5在侵彻过程中的温升较压装药PBX04更低。

  • 图  CMM模型所考虑的微缺陷演化机制与变形机制示意图

    Figure  1.  Conceptual diagram of all kinds of considered mechanisms in the current model

    图  计算与实验得到的两类炸药在不同应变率下的应力-应变曲线

    Figure  2.  Calculated and experimental stress-strain curves of two kinds of explosives at different strain rates

    图  有限元计算模型

    Figure  3.  Finite element calculation model

    图  压装药PBX04的压力演化情况

    Figure  4.  Pressure evolution of PBX04

    图  浇注药GOFL-5的压力演化情况

    Figure  5.  Pressure evolution of GOFL-5

    图  不同时刻两类炸药微裂纹损伤演化情况

    Figure  6.  Evolution of microcrack damage of two kinds of explosives

    图  不同时刻两类炸药的微孔洞损伤演化

    Figure  7.  Evolution of microvoid damage of two kinds of explosives

    图  不同时刻PBX04中微裂纹与微孔洞相关热点温度曲线

    Figure  8.  Microcrack and microvoid related hotspot temperature evolution curves for PBX04

    图  不同时刻GOFL-5中微裂纹与微孔洞相关热点温度曲线

    Figure  9.  Microcrack and microvoid related hotspot temperature evolution curves for GOFL-5

    图  10  微裂纹不同扩展形态示意图

    Figure  10.  Schematic diagram of different propagation patterns of microcracks

    图  11  PBX04装药中不同位置的微裂纹状态频率分布

    Figure  11.  Frequency distribution of microcracks state in different positions of PBX04 charge

    表  1  GOFL-5与PBX04材料参数

    Table  1.   Material parameters for GOFL-5 and PBX04

    Material\;\rho 0/(kg·m−3)G/GPa\nuG1/MPaG2/MPaG3/MPaG4/MPaG5/MPa{\tau{_1^{-1} } } /{\rm{s} }{^{-1} }
    GOFL-51 7500.550.3 167 30.45 90.03 185.6 120.00
    PBX041 8208.250.31 9401 175.001 521.001 909.01 688.00
    Material{\tau{_2^{-1} } } /{\rm{s} }{^{-1} }{\tau{_3^{-1} } } /{\rm{s} }{^{-1} }{\tau{_4^{-1} } } /{\rm{s} }{^{-1} }{\tau{_5^{-1} } } /{\rm{s} }{^{-1} }\sigma{{_0}}/MPaCh/MPan{\bar c}0/μm
    GOFL-57.32 × 1037.32 × 1047.32 × 1057.32 × 106 2.20.76 4.50.4530
    PBX049.00 × 1039.00 × 1049.00 × 1052.00 × 10640.00.101500.01.0030
    MaterialN0/cm−3\bar \gamma /(J·m−2){\dot c_{\max }}/(m·s–1)\;\mu_{\rm{s}}mf0/(m·s–1)kwc0/(m·s–1)s \varGamma
    GOFL-5 30.53000.35.00.012.01 0000.46 0.89
    PBX043001.43000.55.00.012.02 5002.26 1.50
    下载: 导出CSV

    表  2  弹体材料参数

    Table  2.   Parameters of projectile material

    Physical properties Johnson-Cook model
    \;\rho /(g·cm−3)c/(J·kg−1·K−1)\kappa /(kW·m−1·K−1)\alpha /(m·K−1)Tm/℃G/GPa NM
    7.8247838.113.24 × 10−51 793.15 774.970.26 1.03
    下载: 导出CSV
    Johnson-Cook model Damage model
    C1/MPaC2/MPa C3 C4 D1 D2 D3 D4 D5
    792.21509.521.4 0 –0.8 2.1 –0.5 20.0 0.61
    下载: 导出CSV

    表  3  混凝土靶板材料参数

    Table  3.   Material parameters of concrete plate target

    Physical properties HJC model
    \;\rho /(g·cm−3)c/(J·kg−1·K−1)κ/(kW·m−1·K−1)α/(m·K−1)G/GPaFC/MPaABNC
    2.286541.764.32 × 10−5 16.4040.680.751.650.767.0 × 10−3
    下载: 导出CSV
    HJC model Damage model
    pC/MPaUC pL/GPa UL K1/GPa K2/GPa K3/GPa D1 D2 \varepsilon min
    13.565.80 × 10−4 1.05 0.10 17.40 38.80 29.80 0.03 1.0 0.01
    下载: 导出CSV
  • [1] 唐明峰, 李明, 蓝林钢. 浇注PBX力学性能的研究进展 [J]. 含能材料, 2013, 21(6): 812–817. doi: 10.3969/j.issn.1006-9941.2013.06.024

    TANG M F, LI M, LAN L G. Review on the mechanical properties of cast PBXs [J]. Chinese Journal of Energetic Materials, 2013, 21(6): 812–817. doi: 10.3969/j.issn.1006-9941.2013.06.024
    [2] 李媛媛, 南海. 国外浇注PBX炸药在硬目标侵彻武器中的应用 [J]. 飞航导弹, 2012(11): 88–91.

    LI Y Y, NAN H. Application of foreign casting PBX explosive in hard target penetration weapons [J]. Cruise Missile, 2012(11): 88–91.
    [3] 吴会民, 卢芳云, 卢力, 等. 三种含能材料力学行为应变率效应的实验研究 [J]. 含能材料, 2004, 12(4): 40–43.

    WU H M, LU F Y, LU L, et al. Experimental studies on strain rate effects of mechanical behaviors of energetic materials [J]. Chinese Journal of Energetic Materials, 2004, 12(4): 40–43.
    [4] 韩小平, 张元冲, 沈亚鹏, 等. Comp. B复合炸药动态压缩力学性能和本构关系的研究 [J]. 实验力学, 1996, 11(3): 303–310.

    HAN X P, ZHANG Y C, SHEN Y P, et al. Dynamic behavior and constitutive model of Comp. B explosive [J]. Journal of Experimental Mechanics, 1996, 11(3): 303–310.
    [5] 罗景润. PBX的损伤、断裂及本构关系研究[D]. 绵阳: 中国工程物理研究院, 2001.

    LUO J R. Study on damage, fracture and constitutive relation of PBX [D]. Mianyang: China Academy of Engineering Physics, 2001.
    [6] MA X, LI X G, ZHENG X X, et al. Weak shock loadings induce potential hot spots formation around an intergranular pore [J]. Journal of Applied Physics, 2017, 121(11): 115102. doi: 10.1063/1.4978355
    [7] MA X, LI X G, ZHENG X X, et al. Crack initiation and potential hot-spot formation around a cylindrical defect under dynamic compression [J]. Journal of Applied Physics, 2017, 122(18): 185104. doi: 10.1063/1.4998679
    [8] 李英雷, 李大红, 胡时胜, 等. TATB钝感炸药本构关系的实验研究 [J]. 爆炸与冲击, 1999, 19(4): 353–359. doi: 10.3321/j.issn:1001-1455.1999.04.011

    LI Y L, LI D H, HU S S, et al. Experimental study on constitutive relation of TATB insensitive explosive [J]. Explosion and Shock Waves, 1999, 19(4): 353–359. doi: 10.3321/j.issn:1001-1455.1999.04.011
    [9] 李俊玲, 卢芳云, 傅华, 等. 某PBX炸药的动态力学性能研究 [J]. 高压物理学报, 2011, 25(2): 159–164. doi: 10.11858/gywlxb.2011.02.012

    LI J L, LU F Y, FU H, et al. Research on the dynamic behavior of a PBX explosive [J]. Chinese Journal of High Pressure Physics, 2011, 25(2): 159–164. doi: 10.11858/gywlxb.2011.02.012
    [10] TAN H. The cohesive law of particle/binder interfaces in solid propellants [J]. Progress in Propulsion Physics, 2011, 2: 59–66.
    [11] WANG J C, LUO J R. Predicting the effective elastic properties of polymer bonded explosives based on micromechanical methods [J]. Journal of Energetic Materials, 2018, 36(2): 211–222. doi: 10.1080/07370652.2017.1354098
    [12] BENELFEKKAH A, FRACHON A, GRATTON M, et al. Analytical and numerical comparison of discrete damage models with induced anisotropy [J]. Engineering Fracture Mechanics, 2014, 121/122: 28–39. doi: 10.1016/j.engfracmech.2014.03.022
    [13] PICART D, BENELFELLAH A, BRIGOLLE J L, et al. Characterization and modeling of the anisotropic damage of a high-explosive composition [J]. Engineering Fracture Mechanics, 2014, 131: 525–537. doi: 10.1016/j.engfracmech.2014.09.009
    [14] 张馨予, 吴艳青, 黄风雷. PBX装药弹体侵彻混凝土薄板的数值模拟 [J]. 含能材料, 2018, 26(1): 101–108. doi: 10.11943/j.issn.1006-9941.2018.01.013

    ZHANG X Y, WU Y Q, HUANG F L. Numerical simulation on the dynamic damage of PBX charges filled in projectiles during penetrating thin concrete targets [J]. Chinese Journal of Energetic Materials, 2018, 26(1): 101–108. doi: 10.11943/j.issn.1006-9941.2018.01.013
    [15] 石啸海, 余春祥, 戴开达, 等. 侵彻过程中弹头形状对PBX炸药损伤的影响 [J]. 弹箭与制导学报, 2019, 39(3): 81–85.

    SHI X H, YU C X, DAI K D, et al. The influence of nose shape to dynamic damage of PBX charge during the penetration process [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2019, 39(3): 81–85.
    [16] 石啸海, 戴开达, 陈鹏万, 等. 战斗部侵彻过程中PBX装药动态损伤数值模拟 [J]. 中国测试, 2016, 42(10): 138–142. doi: 10.11857/j.issn.1674-5124.2016.10.026

    SHI X H, DAI K D, CHEN P W, et al. Numerical simulation of dynamic damage of PBX charge during the warhead penetration process [J]. China Measurement and Test, 2016, 42(10): 138–142. doi: 10.11857/j.issn.1674-5124.2016.10.026
    [17] 成丽蓉, 汪德武, 贺元吉. 侵彻单层和多层靶时战斗部装药损伤及热点生成机理研究 [J]. 兵工学报, 2020, 41(1): 32–39. doi: 10.3969/j.issn.1000-1093.2020.01.004

    CHENG L R, WANG D W, HE Y J. Research on the damage and hot-spot generation in explosive charges during penetration into single- or multi-layer target [J]. Acta Armamentarii, 2020, 41(1): 32–39. doi: 10.3969/j.issn.1000-1093.2020.01.004
    [18] LI X, LIU Y Z, SUN Y. Dynamic mechanical damage and non-shock initiation of a new polymer bonded explosive during penetration [J]. Polymers, 2020, 12(6): 1342. doi: 10.3390/polym12061342
    [19] YANG K, WU Y Q, HUANG F L. Damage and hotspot formation simulation for impact-shear loaded PBXs using combined microcrack and microvoid model [J]. European Journal of Mechanics-A/Solids, 2020, 80: 103924. doi: 10.1016/j.euromechsol.2019.103924
    [20] YANG K, WU Y Q, HUANG F L. Microcrack and microvoid dominated damage behaviors for polymer bonded explosives under different dynamic loading conditions [J]. Mechanics of Materials, 2019, 137: 103130. doi: 10.1016/j.mechmat.2019.103130
    [21] DANCYGIER A N. Concrete strength effect on penetration and perforation resistance to impact of non-deforming projectiles [C]//8th International Symposium on Plasticity and Impact Mechanics. Delhi, India, 2003: 826–832.
  • 加载中
图(11) / 表(5)
计量
  • 文章访问数:  1010
  • HTML全文浏览量:  639
  • PDF下载量:  79
出版历程
  • 收稿日期:  2021-03-27
  • 修回日期:  2021-04-16

目录

/

返回文章
返回