Loading [MathJax]/jax/output/SVG/jax.js

可变形定向破片战斗部在不同展开模式下的性能数值模拟

杜伊杨 黄炫宁 杨鸿铭 周晴 李伟兵

王志鹏, 韩志军, 王龙飞. 基于Galerkin法研究应力波作用下复合材料板的动力学失稳[J]. 高压物理学报, 2021, 35(5): 054204. doi: 10.11858/gywlxb.20210705
引用本文: 杜伊杨, 黄炫宁, 杨鸿铭, 周晴, 李伟兵. 可变形定向破片战斗部在不同展开模式下的性能数值模拟[J]. 高压物理学报, 2022, 36(3): 035103. doi: 10.11858/gywlxb.20210887
WANG Zhipeng, HAN Zhijun, WANG Longfei. Dynamic Instability of Composite Plate under Stress Wave Based on Galerkin Method[J]. Chinese Journal of High Pressure Physics, 2021, 35(5): 054204. doi: 10.11858/gywlxb.20210705
Citation: DU Yiyang, HUANG Xuanning, YANG Hongming, ZHOU Qing, LI Weibing. Numerical Simulation of Deformable Directional Fragment Warhead Performance in Different Expansion Modes[J]. Chinese Journal of High Pressure Physics, 2022, 36(3): 035103. doi: 10.11858/gywlxb.20210887

可变形定向破片战斗部在不同展开模式下的性能数值模拟

doi: 10.11858/gywlxb.20210887
基金项目: 国家级大学生创新创业训练计划(202010288007Z)
详细信息
    作者简介:

    杜伊杨(1999-),男,本科,主要从事战斗部设计研究. E-mail:duyiyang@njust.edu.cn

    通讯作者:

    李伟兵(1982-),男,博士,教授,博士生导师,主要从事高效毁伤技术研究.E-mail:njustlwb@163.com

  • 中图分类号: O383; TJ410.33

Numerical Simulation of Deformable Directional Fragment Warhead Performance in Different Expansion Modes

  • 摘要: 为使战斗部具有多种定向毁伤模式并实现一定程度上的可控毁伤,提出了一种扇形装药的可变形定向破片战斗部,该战斗部可实现轴向展开和侧向展开2种模式。采用AUTODYN软件进行破片场的数值模拟。首先,基于战斗部单元体分析获得了距离轴心25 mm处的最佳起爆点位置;其次,对整个战斗部进行分析,在轴向展开模式下分析了轴向展开角度对破片飞散速度、破片数目和破片空间分布的影响,发现轴向展开角在60°~75°范围内毁伤效果较佳;最后,在侧向展开模式下分析了整个战斗部的破片速度和破片空间分布情况,结果表明破片具有明显的定向飞散特性。

     

  • 复合材料板具有较高的强度质量比、良好的耐腐蚀性和优异的可设计性,被广泛应用于航空航天和工业制造等领域[1]。在实际使用中复合材料板经常受到不同形式的冲击荷载,从而产生振动和屈曲问题,因此冲击载荷作用下复合材料板的动力稳定性问题备受关注。

    近年来关于复合材料板的研究越来越多,尤其是冲击荷载作用下复合材料板的动力稳定性问题研究[2],对工程部件结构设计和使用具有重要的意义。Sun等[3]研究了在加热环境中应力波对功能梯度圆柱壳轴向冲击屈曲的影响;毛柳伟等[4]对弹性直杆在应力波作用下的动力分叉屈曲进行了分析与探讨,提出了求解应力波作用下直杆动力屈曲的数值方法;Lepik[5]讨论了在应力波影响下轴向压缩的弹塑性梁的屈曲;Zhang等[6]分析了不确定初始几何缺陷对薄板屈曲的影响;Abdelaziz等[7]利用双曲线剪切变形理论,分析了在各种边界条件下复合材料板的弯曲变形和屈曲;Kouchakzadeh等[8]采用线性和旋转弹簧的均匀分布来模拟边界条件,对矩形层压复合板的屈曲进行了分析;Czapski等[9]通过数值和实验方法,研究了残余应力对压缩至破坏期间薄壁层压板屈曲性能的影响。

    在实际工程中,复合材料板多应用于振动环境,其在动力响应下的动态特性和振动分析是必不可少的,因而对该类材料的振动屈曲研究至关重要。Kuo[10]研究了两种非均匀分布纤维复合材料板的振动屈曲问题,Villarreal等[11]对典型正交异性板的本征频率和振动屈曲进行了理论分析,Eftekhari等[12]提出了通过组合应用有限元方法和微分正交方法求解矩形板的振动屈曲问题,Rehman等[13]探讨了壳体结构的缺陷和损坏对结构振动屈曲的影响,Sayyad等[14]将三角剪切变形理论应用于复合板的变形和振动屈曲研究。

    关于复合材料板的振动屈曲问题已开展了较多的研究,但大多未考虑应力波效应对振动屈曲的影响,而动力屈曲一般与应力波相联系且具有局部发生的特点,研究含初始缺陷的复合材料板能更好地揭示实际工程中复合材料板在不同工况下发生动力屈曲的机理。基于此,本研究利用Kirchhoff薄板理论和Hamilton原理,考虑应力波效应,建立含初始几何缺陷的四边简支复合材料板的振动控制方程,得到板的屈曲临界荷载表达式,在此基础上通过数值计算讨论初始几何缺陷、振型函数初相位、铺层角度、屈曲模态阶数和铺层层数对复合材料板振动屈曲临界荷载的影响,为工程实际提供理论依据。

    复合材料板在x = La处为固定边界条件,其余3边为简支边界条件,在z = 0的中性面上受x方向的面内阶跃荷载N作用,如图1所示。板在z方向上含初始几何缺陷w1,且xyz方向的位移分别为uvw。复合材料板的长度、宽度和厚度分别为LaLbh,由n层单层板组成,θ为单层板的铺层角度,即纤维材料铺设方向与x方向的夹角(见图1)。

    图  1  复合材料板结构示意图
    Figure  1.  Schematic diagram of composite plate structure

    根据Kirchhoff薄板理论及经典弹性理论,复合材料板的位移与应变、弹性曲面的曲率和扭率的表达式为

    {w=w1+w0u=u0z(ww1)xv=v0z(ww1)y (1)
    {ε0x=u0xε0y=v0yγ0xy=u0y+v0x (2)
    {κx=2(ww1)x2κy=2(ww1)y2κxy=22(ww1)xy (3)
    {εx=ε0x+zκxεy=ε0y+zκyγxy=γ0xy+zκxy (4)

    式中:u0v0w0分别为复合材料板在xyz方向上的中面位移,ε0xε0yγ0xy分别为复合材料板中面应变分量,κxκyκxy为中面的曲率和扭率,εxεyγxy为复合材料板任意一点的应变。

    板的内力(NxNyNxy)与内力矩(MxMyMxy)之间的关系为

    [NxNyNxyMxMyMxy]=[A11A12A16B11B12B16A12A22A26B12B22B26A16A26A66B16B26B66B11B12B16D11D12D16B12B22B26D12D22D26B16B26B66D16D26D66][ε0xε0yγ0xyκxκyκxy] (5)

    式中:AijBijDij(i,j=1,2,6)分别表示板的拉伸刚度、耦合刚度和弯曲刚度系数[15],表达式如下

    {Aij=nk=1¯Qkij(hkhk1)Bij=12nk=1¯Qkij(h2kh2k1)Dij=13nk=1¯Qkij(h3kh3k1) (6)

    式中:¯Qkij为复合材料板第k层的偏轴刚度系数

    ¯Q=P1Q(P1)T (7)

    式中:P为坐标转换矩阵,Q为刚度矩阵。

    P=[cos2θsin2θ2sinθcosθsin2θcos2θ2sinθcosθsinθcosθsinθcosθcos2θsin2θ] (8)
    Q=[Q11Q120Q12Q22000Q66] (9)

    考虑材料为正交各向异性材料,设E1E2G12μ12μ21分别为板材料xy方向的拉压弹性模量、剪切弹性模量、主泊松比和副泊松比,则有

    Q11=E11μ12μ21,Q22=E21μ12μ21,Q12=μ12E21μ12μ21=μ21E11μ12μ21,Q66=G12 (10)

    复合材料板在左端受面内冲击荷载N作用(见图1)时,应力波沿x方向在板内传播,其应力变化如图2所示。

    图  2  应力波传播示意图
    Figure  2.  Schematic diagram of stress wave propagation

    当应力波传播至波阵面位置Lcr(临界长度)时,板发生振动屈曲,板的内力Nt和应力波波速c分别表示为

    Nt={N0xLcr0x>Lcr,c=A11ρh (11)

    板发生振动屈曲时的变形能可以表示为

    U=12Lcr0Lb0(Nxε0x+Nyε0y+Nxyγ0xy+Mxκx+Myκy+Mxyκxy)dxdy=12Lcr0Lb0[Nxu0x+Nyv0y+Nxy(u0x+v0y)+Mx(2w0x2)+My(2w0y2)+Mxy(22w0xy)]dxdy (12)

    发生振动屈曲时的动能(考虑转动惯量)可以表示为

    T=12h/2h/2Lcr0Lb0ρ(k)[(ut)2+(vt)2+(w0t)2]dxdydz=12Lcr0Lb0[I0(u0t)2+I0(v0t)2+I0(w0t)22I1u0t2w0xt2I1v0t2w0yt+I2(2w0xt)2+I2(2w0yt)2]dxdy (13)

    式中:(I0,I1,I2)=Nkk=1hkhk1ρ(k)(1,z,z2)dzρ(k)为第k层材料的密度。

    板发生振动屈曲时的外力功可以表示为

    W=12Lcr0Lb0Nt(wx)2dxdy (14)

    考虑Hamilton变分原理,即

    δt2t1(TU+W)dt=0 (15)

    将式(12)~式(14)代入式(15)并进行变分计算,可得

    Nxx+NxyxI02u0t2+I13w0xt2=0 (16)
    Nyy+NxyyI02v0t2+I13w0yt2=0 (17)
    2Mxx2+2Myy222MxyxyNt(2w0x2+2w1x2)=I02w0t2+I1(3u0xt2+3v0xt2)I2(4w0x2t2+4w0y2t2) (18)

    对于正交各向异性正规对称正交铺设的复合材料板,其刚度矩阵满足[16]

    D16=D26=0,A16=A26=0,Bij=0 (19)

    根据 Kirchhoff 薄板理论及经典弹性理论,薄板中面在变形过程中没有伸长变形,将板的本构关系代入式(16)~式(18)中,略去含u0v0的项,得到复合材料板在面向阶跃荷载激励下的控制方程

    D114w0x4+(2D12+4D66)4w0x2y2+D224w0y4+Nt(2w0x2+2w1x2)I02w0t2+I2(4w0x2t2+4w0y2t2)=0 (20)

    设三边简支和应力波传播到Lcrw0(Lcr,y,t)=w0(Lcr,y,t)=0)时的振型函数[17]

    w0(x,y,t)=Rij[siniπxLcr+ii+1sin(i+1)πxLcr]sinjπyLbsin(ωt+φ) (21)

    由屈曲模态确定的缺陷分布形式是板结构最有可能发生的屈曲形式,能够很好地确定结构的缺陷敏感性[18]。对于复合材料板在制造过程中出现的初始几何缺陷,引入屈曲模态的ε倍变形作为初始几何缺陷[19],可以表示为

    w1(x,y)=εRij[siniπxLcr+ii+1sin(i+1)πxLcr]sinjπyLb (22)

    式中:ij为屈曲模态阶数,i,j=1,2,3,Rij为板的第(i, j)阶模态幅值;ε表示初始几何缺陷系数。

    根据式(21)和式(22),利用棣莫弗公式对控制方程式(20)中的各项进行求导并化简,得到

    D11(πLcr)4[i4+(i+1)4]sin(ωt+φ)+(2D12+4D66)(πLcr)2[i2+(i+1)2](jπLb)2sin(ωt+φ)+2D22(jπLb)4sin(ωt+φ)Nt(πLcr)2[i2+(i+1)2][sin(ωt+φ)+ε]+2I0ω2sin(ωt+φ)+I2ω2{(πLcr)2[i2+(i+1)2]+2(jπLb)2}sin(ωt+φ)=0 (23)

    根据式(23)可以得到N的表达式

    N={D11(πLcr)4[i4+(i+1)4]+(2D12+4D66)(πLcr)2[i2+(i+1)2](jπLb)2+2D22(jπLb)4+2I0ω2+I2ω2{(πLcr)2[i2+(i+1)2]+2(jπLb)2}}sin(ωt+φ)(πLcr)2[i2+(i+1)2][sin(ωt+φ)+ε] (24)

    板发生屈曲时,临界条件为ω=0[20],代入式(24)可得振动屈曲临界荷载为

    Ncr=D11(πLcr)2[i4+(i+1)4]sinφ[i2+(i+1)2](sinφ+ε)+(2D12+4D66)(jπLb)2sinφ(sinφ+ε)+2D22(j2πL2b)2L2crsinφ[i2+(i+1)2](sinφ+ε) (25)

    利用MATLAB数值分析应力波未反射时初始几何缺陷、初相位、铺层角度、屈曲模态阶数、铺设厚度以及铺层层数对复合材料板振动屈曲临界荷载的影响,使用的材料参数见表1[21]

    表  1  复合材料板参数[21]
    Table  1.  Material parameters of composite plate[21]
    E1/GPaE2/GPaG12/GPaμ12La/mLb/m
    140.08.65.00.350.600.50
    下载: 导出CSV 
    | 显示表格

    图3显示了复合材料板的模态阶数i为1、2、3,j为1时的屈曲模态。当板的x方向模态增大时,x方向的屈曲模态第一峰值增大且波数增加,而y方向的屈曲模态呈正对称分布。模态阶数的增加使板振动的屈曲模态变得更复杂。

    图  3  x方向模态取值增大时板的屈曲模态
    Figure  3.  Buckling mode of composite plate with increasing mode value in x direction

    设置7组算例,分别以初始几何缺陷、初相位、铺层角度、xy两个方向屈曲模态阶数、铺层层数及铺设厚度为变量进行算例分析,研究以上因素对板振动屈曲临界荷载的影响,算例参数见表2

    表  2  算例分析参数表
    Table  2.  Example analysis parameter table
    GroupInitial defect
    coefficient
    Order of modeLaying angle/(°)Initial phaseNumber of
    layers laid
    Thickness
    of the plate/m
    x directiony direction
    AVariablei = 1j = 1[0, 0, 0, 0, 0]π/250.01
    B0.1Variablej = 1[0, 0, 0, 0, 0]π/250.01
    C0.1i =1Variable[0, 0, 0, 0, 0]π/250.01
    D0.1i =1j = 1Variableπ/250.01
    E0.1i =1j = 1[0, 0, 0, 0, 0]Variable50.01
    F0.1i =1j = 1[0, 0, 0, 0, 0]π/2Variable0.01
    G0.1i =1j = 1[0, 0, 0, 0, 0]π/25Variable
    下载: 导出CSV 
    | 显示表格

    将A组数据代入式(25)中,可以得到不同初始缺陷系数对复合材料板振动屈曲的影响,如图4所示。由Lcr-Ncr曲线可知:在应力波传播过程中,Ncr呈指数下降,分为敏感区和非敏感区。应力波在Lcr < 0.4 m区域传播时,Ncr的变化较陡峭,该区域为敏感区;应力波在Lcr > 0.4 m区域传播时,Ncr的变化趋于平缓,该区域为非敏感区,因此敏感分界点为0.4 m。当选取的初始缺陷系数增大时,临界荷载Ncr也随之增大。在敏感区,初始缺陷系数对临界荷载Ncr的影响较大,且随应力波传播呈减小趋势。此外,初始缺陷系数对非敏感区的影响较小。图4表明,初始几何缺陷系数越大,板越容易发生屈曲。

    图  4  不同初始缺陷系数条件下NcrLcr的关系曲线
    Figure  4.  Relationship between Ncr and Lcr under different initial defect coefficients

    将B组数据代入式(25)中,可以得到不同x方向模态阶数对复合材料板振动屈曲的影响,如图5所示。由Lcr-Ncr曲线可知:当选取的x方向模态阶数增大时,临界荷载Ncr随之明显增大。在敏感区,x方向模态阶数对临界荷载Ncr的影响很大,并随应力波的传播不断减小,到达非敏感区之后影响较小并趋于稳定。图5表明,x方向模态阶数的增加会显著增大板的屈曲临界荷载。

    图  5  x方向模态阶数不同时NcrLcr的关系曲线
    Figure  5.  Relationship between Ncr and Lcr with different order of modes in x direction

    将C组数据代入式(25)中,可以得到不同的y方向模态阶数对复合材料板振动屈曲的影响,如图6所示。由Lcr-Ncr曲线可知:当选取的y方向模态阶数增大时,临界荷载Ncr也随之变大。在应力波传播过程中,在敏感区y方向模态阶数对临界载荷基本没有影响,而在非敏感区有极小的影响。图6表明,y方向模态阶数的变化对板屈曲临界荷载基本没有影响。

    图  6  y方向模态阶数不同时NcrLcr的关系曲线
    Figure  6.  Relationship between Ncr and Lcr with different order of modes in y direction

    将D组数据代入式(25)中,得到不同铺层角度对复合材料板振动屈曲的影响,如图7所示。由Lcr-Ncr曲线可知:在敏感区,不同的铺设角度对临界荷载Ncr的影响较大,且随应力波的传播不断减小,到达非敏感区后趋于稳定。图7表明,铺层角度小的单层板的层数越多,板的临界屈曲荷载越大,说明复合材料板的铺设角度直接影响板的屈曲临界荷载。

    图  7  不同铺层角度条件下NcrLcr的关系曲线
    Figure  7.  Relationship between Ncr and Lcr under different laying angles

    将E组数据代入式(25)中,得到不同初相位对复合材料板振动屈曲的影响,如图8所示。由Lcr-Ncr曲线可知:振型函数的初相位越大,对应的临界荷载越大。在敏感区,振型函数的初相位对临界荷载Ncr的影响较小,且随应力波的传播不断减小;到达非敏感区之后,影响趋于平缓。图8表明,振型函数的初相位越大,板的屈曲临界荷载越大。

    图  8  不同初相位条件下NcrLcr的关系曲线
    Figure  8.  Relationship between Ncr and Lcr under the condition of initial phase of different mode functions

    将F组数据代入式(25)中,得到不同铺层层数对复合材料板振动屈曲的影响,如图9所示。由Lcr-Ncr曲线可知:当按照不同层数铺设时,敏感区的临界荷载Ncr的变化较大,且随应力波的传播不断减小;到达非敏感区之后变化较小并趋于平缓。图9表明,对于厚度固定且对称铺设的板,当铺设层数达到7时,其屈曲荷载随层数增加趋于稳定。

    图  9  不同铺层层数下NcrLcr的关系曲线
    Figure  9.  Relationship between Ncr and Lcr under different laying modes

    将G组数据代入式(25)中,得到不同铺设厚度对复合材料板振动屈曲的影响,如图10所示。由Lcr-Ncr曲线可知:板的铺设厚度越大,对应的临界荷载越大。在敏感区,不同的板厚对临界荷载Ncr的影响很大,且随应力波的传播不断减小,到达非敏感区后趋于稳定。图10表明,复合材料板的铺设厚度将直接决定板的屈曲临界荷载。

    图  10  不同铺设厚度下NcrLcr的关系曲线
    Figure  10.  Relationship between Ncr and Lcr under different thicknesses

    基于Kirchhoff薄板理论和Hamilton变分原理,建立了具有初始几何缺陷的四边简支复合材料板的振动控制方程。采用伽辽金法,选取符合边界条件的振型函数求解控制方程,得到屈曲临界载荷表达式。数值计算结果表明:应力波在未发生反射前的传播过程中,复合材料板的振动屈曲临界载荷随着临界长度的增大、铺设厚度的减小、初始几何缺陷系数的增大、振型函数初相位的减小而减小;复合材料板的各层铺层角度与荷载作用方向的夹角越小,屈曲临界载荷越大,当对称铺设层数达7层时,临界荷载趋于稳定。研究结果可为工程中复合材料板的结构设计与应用提供一定的参考。

  • 图  战斗部结构示意图

    Figure  1.  Schematic diagram of warhead structure

    图  展开模式示意图

    Figure  2.  Schematic diagram of expansion modes

    图  战斗部的简化模型示意图(单位:mm)

    Figure  3.  Schematic diagram of warhead model (Unit: mm)

    图  选取的起爆点位置示意图

    Figure  4.  Schematic diagram of initiation point locations

    图  轴向展开角示意图

    Figure  5.  Schematic diagram of axial-expansion angle

    图  侧向展开角示意图

    Figure  6.  Schematic diagram of lateral-expansion angle

    图  不同起爆点位置下破片的飞散速度

    Figure  7.  Fragments’ mean velocity for different initiation points

    图  10 μs时不同起爆点位置下壳体所受压强分布云图

    Figure  8.  Pressure distribution on warhead case for different initiaion points at a time delay of 10 μs

    图  不同起爆点位置下的破片数统计

    Figure  9.  Fragment quantity for different initiation points

    图  10  100 μs时单元体破片的飞散情况

    Figure  10.  Fragment dispersion of unit constituting warhead at different initiation points at a time delay of 100 μs

    图  11  不同展开角下破片数统计

    Figure  11.  Fragment quantity at different axial-expansion angles

    图  12  不同轴向展开角下破片位置的散点图

    Figure  12.  Scatter plots of fragments’ locations at different axial-expansion angles

    图  13  不同轴向展开角下破片的平均速度及其Z轴的分量

    Figure  13.  Mean velocities and velocities in Z axis of fragments at different expansion angles

    图  14  破片密集区

    Figure  14.  Dense area of fragments

    图  15  不同轴向展开角下破片的分布情况

    Figure  15.  Fragment dispersions at different axial-expansion angles

    图  16  破片盲区锥顶角示意图

    Figure  16.  Schematic diagram of the cone angle for the zone without fragment

    图  17  不同轴向展开角下破片盲区的锥顶角

    Figure  17.  Cone angles of zone without fragment at different axial-expansion angles

    图  18  侧平面破片飞散角示意图

    Figure  18.  Schematic diagram of fragments’ dispersion angle in lateral plane

    图  19  不同轴向展开角度下破片的侧平面飞散角

    Figure  19.  Fragments’ dispersion angles in lateral plane for different axial-expansion angles

    图  20  不同轴向展开角下破片速度随周向角的分布情况

    Figure  20.  Fragments’ velocity distribution in versus of ϕ for different axial-expension angles

    图  21  不同轴向展开角下破片速度随侧平面飞散角的分布情况

    Figure  21.  Fragments’ velocity distributions in versus of δ for different axial-expansion angles

    图  22  侧向展开模式下破片飞散情况

    Figure  22.  Fragment dispersion in lateral-expanding mode

    图  23  侧向展开模式下破片位置散点图

    Figure  23.  Scatter plot of fragments’ locations in lateral-expanding mode

    图  24  侧向展开模式下破片速度分布

    Figure  24.  Velocity distribution of fragments in lateral-expanding mode

    表  1  材料模型

    Table  1.   Material models

    ComponentMaterialEquation of stateStrength modelFailure model
    ExplosiveComp. BJWLHydro
    CasingSteel 4340LinearJohnson-CookPrincipal strain
    下载: 导出CSV

    表  2  Comp. B炸药材料参数

    Table  2.   Material parameters of Comp. B

    ρ/(g·cm−3)D/(m·s−1)pCJ/GPa A/GPaB/GPaR1R2ω
    1.717798029.5524.27.6784.21.10.34
    下载: 导出CSV

    表  3  4340钢材料参数

    Table  3.   Material parameters of steel 4340

    ρ/(g·cm−3)IYS/GPaHC/GPaHESRCRSRBM/GPaSM/GPa
    7.830.7920.510.260.014115981.8
    下载: 导出CSV

    表  4  密集区破片的平均速度

    Table  4.   Mean velocity of fragments in dense area

    β/(°)Mean velocity/(m∙s−1)
    d=10 mmd=20 mmd=30 mmd=40 mmd=50 mmd=60 mm
    60118011821225128713111323
    75130213321343136613761401
    90141314911503151515241527
    下载: 导出CSV

    表  5  侧平面飞散角在不同取值区间的破片占比

    Table  5.   Distribution of fragments at different intervals of dispersion angle in the lateral plane

    β/(°)δmin/(°)δmax/(°)xδ/%
    δ<10°10°
    −20°
    20°
    −30°
    30°
    −40°
    40°
    −50°
    50°
    −60°
    60°
    −70°
    70°
    −80°
    80°
    −90°
    δ>90°
    15−3.5433.4733.2641.3422.253.15000000
    304.2455.05 6.927.1445.7631.03 4.02 5.130000
    4516.7669.0704.345.8428.3543.2913.42 4.76000
    6030.9884.5700013.1213.6242.1625.45 3.342.310
    7549.9991.84000010.6021.6627.1929.959.45 1.15
    9047.75101.000000 0.17 8.9525.9923.7529.95 11.19
    下载: 导出CSV

    表  6  飞散角γη在不同取值区间内的碎片数占比

    Table  6.   Fragment distributions at different intervals of dispersion angles γ and η

    xγ/%γmax/(°)
    0°−5°5°−10°10°−15°15°−20°20°−25°25°−30°30°−35°35°−40°40°−45°
    57.6823.419.883.292.320.730.850.611.2344.24
    xη/%ηmax/(°)
    0°−5°5°−10°10°−15°15°−20°20°−25°25°−30°30°−35°35°−40°40°−45°
    54.6330.005.242.325.981.590.240031.34
    下载: 导出CSV
  • [1] 崔瀚, 张国新. 定向战斗部研究现状及展望 [J]. 飞航导弹, 2019(3): 84–89.
    [2] 李记刚, 余文力, 王涛. 定向战斗部的研究现状及发展趋势 [J]. 飞航导弹, 2005(5): 25–29. doi: 10.3969/j.issn.1009-1319.2005.05.008
    [3] 张天光. 美英定向战斗部的研究与应用 [J]. 航空兵器, 2002, 9(3): 38–41. doi: 10.3969/j.issn.1673-5048.2002.03.013
    [4] SAVAGE J C, O'NEAL J K, BROWN R A, et al. Powered low cost autonomous attack system: cooperative, autonomous, wide-area-search munitions with capability to serve as non-traditional ISR assets in a network-centric environment [C]//Proceedings of SPIE 5791, Laser Radar Technology and Applications X. Orlando: SPIE, 2005.
    [5] 谭多望, 张振宇, 王志兵, 等. 爆炸变形战斗部模型试验研究 [J]. 爆炸与冲击, 2002, 22(3): 247–251. doi: 10.3321/j.issn:1001-1455.2002.03.010

    TAN D W, ZHANG Z Y, WANG Z B, et al. Investigation of a detonatively deformable warhead [J]. Explosion and Shock Waves, 2002, 22(3): 247–251. doi: 10.3321/j.issn:1001-1455.2002.03.010
    [6] MILLER C W. Warhead having selectable axial effects: United States Patent Application, 20150033971 [P]. 2015-02-05.
    [7] 凌琦, 何勇, 何源, 等. 扇形复合装药驱动破片定向飞散的数值模拟 [J]. 高压物理学报, 2017, 31(05): 557–565.

    LING Q, HE Y, HE Y, et al. Numerical simulation of directed scattering of fragments driven by sector-shaped double-layer charge [J]. Chinese Journal of High Pressure Physics, 2017, 31(05): 557–565.
    [8] 洪晓文, 李伟兵, 李文彬, 等. 展开角度及起爆位置对轴向展开式定向战斗部性能的影响 [J]. 含能材料, 2018, 26(5): 383–389.

    HONG X W, LI W B, LI W B, et al. Effect of expanding angle and initiation position on the performance of axial-expanding directional warhead [J]. Chinese Journal of Energetic Materials, 2018, 26(5): 383–389.
    [9] 陈闯, 杨丽, 阮凌飞, 等. 展开式EFP战斗部复合破片场成型的数值模拟 [J]. 火炸药学报, 2019, 42(3): 295–302.

    CHEN C, YANG L, RUAN L F, et al. Numerical simulation on composite fragment field formation of expansion EFP warhead [J]. Chinese Journal of Explosives & Propellants, 2019, 42(3): 295–302.
    [10] 赵宇哲, 李健, 马天宝. 展开式定向战斗部展开过程实验研究 [J]. 高压物理学报, 2016, 30(2): 116–122.

    ZHAO Y Z, LI J, MA T B. Experiment on spread processes of the spreadable aimed warhead [J]. Chinese Journal of High Pressure Physics, 2016, 30(2): 116–122.
    [11] 耿荻, 马天宝, 宁建国. 定向战斗部爆炸驱动规律研究 [J]. 高压物理学报, 2013, 27(5): 685–690. doi: 10.11858/gywlxb.2013.05.005

    GENG D, MA T B, NING J G. Study on laws of explosive driven behaviors of aimed warhead [J]. Chinese Journal of High Pressure Physics, 2013, 27(5): 685–690. doi: 10.11858/gywlxb.2013.05.005
    [12] 马征, 宁建国, 马天宝. 展开型定向战斗部展开过程的理论计算与数值模拟 [J]. 弹箭与制导学报, 2007, 27(1): 111–114. doi: 10.3969/j.issn.1673-9728.2007.01.033

    MA Z, NING J G, MA T B. Dynamic analysis and numerical simulation for spread process of a evolvable aimed warhead [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2007, 27(1): 111–114. doi: 10.3969/j.issn.1673-9728.2007.01.033
    [13] AUTODYN, Autodyn user’s manual [M]. Release 19.0, 2019.
  • 加载中
图(25) / 表(6)
计量
  • 文章访问数:  1762
  • HTML全文浏览量:  483
  • PDF下载量:  43
出版历程
  • 收稿日期:  2021-09-30
  • 修回日期:  2021-10-15
  • 录用日期:  2021-10-20
  • 刊出日期:  2022-05-30

目录

/

返回文章
返回