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

杆式钨合金弹超高速撞击薄靶的能量损耗

孙欢腾 李名锐 周刚 马坤 舒孝鸿

王志鹏, 韩志军, 王龙飞. 基于Galerkin法研究应力波作用下复合材料板的动力学失稳[J]. 高压物理学报, 2021, 35(5): 054204. doi: 10.11858/gywlxb.20210705
引用本文: 孙欢腾, 李名锐, 周刚, 马坤, 舒孝鸿. 杆式钨合金弹超高速撞击薄靶的能量损耗[J]. 高压物理学报, 2019, 33(6): 064106. doi: 10.11858/gywlxb.20190732
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: SUN Huanteng, LI Mingrui, ZHOU Gang, MA Kun, SHU Xiaohong. Energy Dissipation of Tungsten Alloys Cylindrical Rods Hypervelocity Impacting Thin Steel Target[J]. Chinese Journal of High Pressure Physics, 2019, 33(6): 064106. doi: 10.11858/gywlxb.20190732

杆式钨合金弹超高速撞击薄靶的能量损耗

doi: 10.11858/gywlxb.20190732
基金项目: 国家自然科学基金(11772269, 11402213)
详细信息
    作者简介:

    孙欢腾(1994-),男,博士研究生,主要从事超高速撞击能量研究. E-mail:780519755@qq.com

    通讯作者:

    周 刚(1964-),男,研究员,主要从事爆炸与冲击动力学研究. E-mail:gzhou@nint.ac.cn

  • 中图分类号: O385

Energy Dissipation of Tungsten Alloys Cylindrical Rods Hypervelocity Impacting Thin Steel Target

  • 摘要: 超高速撞击过程伴随着复杂的物理过程。为分析杆式圆柱形钨合金弹超高速撞击薄钢靶时的物理过程,采用AUTODYN/SPH数值仿真计算方法获得了撞击过程模型及每个光滑粒子流体动力学信息,并通过广度搜索破片识别程序识别每个破片所含粒子,利用MATLAB编程对破片粒子数据信息进行统计分析,获得弹靶撞击过程的变化特性、弹靶破片数量、相关能量随撞击时间的变化规律。通过分析发现:随着弹体撞击速度的增加,剩余弹体被严重侵蚀,且弹体能量损耗增加,弹体损失的能量主要转变为弹靶破片动能;计算得到了撞击20 μs时的能量损耗直方图,同时分析了发生撞击时靶板的能量变化过程,并简要描述了该过程。

     

  • 复合材料板具有较高的强度质量比、良好的耐腐蚀性和优异的可设计性,被广泛应用于航空航天和工业制造等领域[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.  Diagrammatic sketch of projectile impact target plate

    图  仿真结果与实验结果对比

    Figure  2.  Comparison between simulation and experimental results

    图  不同时刻弹体的侵蚀

    Figure  3.  Erosion of projectile at different time

    图  不同长度弹体以不同速度撞击时产生的弹体破片数量

    Figure  4.  Fragments’ number of projectile with different velocities and lengths

    图  不同长度弹体的撞击结果

    Figure  5.  Impact results of projectiles with different lengths

    图  不同长度的弹体撞击形成的破片云长度和宽度

    Figure  6.  Debris length and width of projectiles with different lengths

    图  能量耗散图

    Figure  7.  Image of energy dissipation

    图  耗散的能量损失与剩余弹体速度降低量

    Figure  8.  The dissipation of energy and the decline of projectile velocity

    图  耗散的总塑性功与内能

    Figure  9.  The dissipation of total plastic energy and internal energy

    图  10  靶板获能过程

    Figure  10.  The gain energy of target plate

    图  11  剩余弹体动能

    Figure  11.  Kinetic energy of residual projectile

    表  1  状态方程参数

    Table  1.   Parameters of equation of state

    MaterialSC0/(m·s–1)Γρ0/(kgm3)
    Tungsten alloy1.234 0401.6717.6
    Q345 Steel1.494 5692.177.83
    下载: 导出CSV

    表  2  钨合金的Steinberg-Guinan强度模型参数

    Table  2.   Steinberg-Guinan strength model parameters of tungsten alloy

    G0/GPaY0/GPaTm/KGpGT/(MPaK1)βnYp
    1321.44 5201.794–401.30.10.019 027
    下载: 导出CSV

    表  3  Q345钢的Johnson-Cook强度模型参数

    Table  3.   Johnson-Cook strength model parameters of Q345 steel

    A/GPaB/GPanCmTmelt/KT0/KG/GPa
    0.3740.795 70.454 50.015 860.885 61 75930080.47
    下载: 导出CSV

    表  4  实验与数值仿真结果

    Table  4.   Results of simulation and experiment

    MethodResidual projectile’s
    kinetic energy/(km·s–1)
    Residual projectile’s
    length/mm
    Experiment2.94512.510
    Simulation2.95211.823
    下载: 导出CSV
  • [1] 武强, 张庆明, 孙浩勇, 等. 超高速撞击下PTFE/Al含能材料薄板的载荷特性分析 [J]. 航天器环境工程, 2017, 34(1): 1–7. doi: 10.3969/j.issn.1673-1379.2017.01.001

    WU Q, ZHANG Q M, SUN H Y, et al. The loading characteristics of PTFE/Al energetic materialsunder hypervelocity impact [J]. Spacecraft Environment Engineering, 2017, 34(1): 1–7. doi: 10.3969/j.issn.1673-1379.2017.01.001
    [2] 迟润强, 庞宝君, 何茂坚. 球形弹丸超高速正撞击薄板破碎状态实验研究 [J]. 爆炸与冲击, 2009, 29(3): 231–236. doi: 10.3321/j.issn:1001-1455.2009.03.002

    CHI R Q, PANG B J, HE M J. Experimental investigation for deformation and fragmentation of spheres penetrating sheets at hypervelocity [J]. Explosion and Shock Waves, 2009, 29(3): 231–236. doi: 10.3321/j.issn:1001-1455.2009.03.002
    [3] 曲广吉, 韩增尧. 空间碎片超高速撞击动力学建模与数值仿真技术 [J]. 中国空间科学技术, 2002(5): 26–30. doi: 10.3321/j.issn:1000-758X.2002.05.005

    QU G J, HAN Z Y. Dynamical modeling and numerical simulation of hypervelocity space debris impact [J]. Chinese Space Science and Technology, 2002(5): 26–30. doi: 10.3321/j.issn:1000-758X.2002.05.005
    [4] 张伟, 庞宝君, 贾斌, 等. 弹丸超高速撞击防护屏碎片云数值模拟 [J]. 高压物理学报, 2004, 18(1): 47–52. doi: 10.3969/j.issn.1000-5773.2004.01.009

    ZHANG W, PANG B J, JIA B, et al. Numerical simulation of debris cloud produced by hypervelocity impact of projectile on bumper [J]. Chinese Journal of High Pressure Physics, 2004, 18(1): 47–52. doi: 10.3969/j.issn.1000-5773.2004.01.009
    [5] 童宗保, 王金相, 彭楚才. 2~3km/s初速下球形破片对金属薄靶的侵彻效应分析 [J]. 科学技术与工程, 2014, 14(33): 215–219. doi: 10.3969/j.issn.1671-1815.2014.33.040

    TONG Z B, WANG J X, PENG C C. Analysis of the effect of the circle fragment penetration the tine metal target with the initial velocity between 2km/s and 3km/s [J]. Science Technology and Engineering, 2014, 14(33): 215–219. doi: 10.3969/j.issn.1671-1815.2014.33.040
    [6] FA W K, JIE H. Test study on the performance of shielding configuration with stuffed layer under hypervelocity impact [J]. Acta Astronautica, 2016(127): 553–560.
    [7] SHI C L, YI L. Research on the technique of identifying debris and obtaining characteristic parameters of large-scale 3D point set [J]. International Journal of Impact Engineering, 2013(56): 27–31.
    [8] 唐恩凌, 徐名扬, 张庆明. 超高速撞击厚靶过程的能量分配研究 [J]. 固体力学学报, 2016, 37(2): 152–160.

    TANG E L, XU M Y, ZHANG Q M. Study on partitioning of energy in hypervelocity impact on thick target [J]. Chinese Journal o f Solid Mechanics, 2016, 37(2): 152–160.
    [9] Thomson W T. An approximate theory of armor penetration [J]. Journal of Applied Physics, 1955, 26: 80–82. doi: 10.1016/0734-743X(95)00023-4
    [10] WOODWARD R L. The penetration of metal targets by conical projectiles [J]. International Journal of Mechanical Sciences, 1978, 20: 349–359. doi: 10.1016/0020-7403(78)90038-3
    [11] WIERZBICKI T. Petalling of plates under explosive and impact loading [J]. International Journal of Impact Engineering, 1999, 22: 935–954. doi: 10.1016/S0734-743X(99)00028-7
    [12] WANG B, XIONG J. Energy absorption efficiency of carbon fiber reinforced polymer laminates under high velocity impact [J]. Materials and Design, 2013(50): 140–148.
    [13] 汪庆桃, 吴克刚, 陈志阳. 圆柱形长杆超高速正碰撞薄板结构破碎效应 [J]. 振动与冲击, 2017, 36(5): 53–60.

    WANG Q T, WU K G, CHEN Z Y. Fragmentation effect of a long cylindrical rod with a hypervelocity normally impacting a thin plate structure [J]. Journal of Vibration and Shock, 2017, 36(5): 53–60.
    [14] 汪庆桃, 吴克刚, 李必红, 等. 球形弹丸超高速碰撞破碎特性 [J]. 振动与冲击, 2014, 33(22): 46–50.

    WANG Q T, WU K G, LI B H, et al. Fragmentation properties of spherical projectiles impacting at hypervelocity [J]. Journal of Vibration and Shock, 2014, 33(22): 46–50.
    [15] 陈鸿, 周智炫, 黄洁. 带迎角钨杆斜侵彻铝板数值仿真研究 [J]. 实验流体力学, 2014, 28(2): 85–89. doi: 10.11729/syltlx2014pz29

    CHEN H, ZHOU Z X, HUANG J. Numerical simulation study on tungsten rod oblique penetrating aluminum plate with attack angle [J]. Journal of Experiments in Fluid Mechanics, 2014, 28(2): 85–89. doi: 10.11729/syltlx2014pz29
    [16] 王芳, 李磊, 金俊. 空间碎片超高速碰撞的数值方法研究 [J]. 中国安全科学学报, 2006, 16(5): 26–30. doi: 10.3969/j.issn.1003-3033.2006.05.005

    WANG F, LI L, JIN J. Study on numerical simulation of hypervelocity impact of debris in outer space [J]. China Safety Science Journal, 2006, 16(5): 26–30. doi: 10.3969/j.issn.1003-3033.2006.05.005
  • 加载中
图(11) / 表(4)
计量
  • 文章访问数:  7166
  • HTML全文浏览量:  3304
  • PDF下载量:  39
出版历程
  • 收稿日期:  2019-03-04
  • 修回日期:  2019-03-24

目录

/

返回文章
返回