Study on Anti-Fragment Impact Performance of Carbon Fiber Reinforced Plastics
-
摘要: 为探究碳纤维复合材料的抗破片侵彻性能,进行了8 g立方体钢破片侵彻厚度分别为5、10、15 mm碳纤维复合材料靶板试验,获取破片着靶速度,并回收破片、靶板,观测靶板破坏情况。根据试验情况进行数值模拟,探究破片侵彻碳纤维复合材料靶板的过程及碳纤维复合材料靶板的破坏机理,由此得到靶板破坏形式与破片速度的关系:破片速度大于弹道极限时,靶板以纤维的剪切破坏为主;破片速度逐渐降低时,纤维的拉伸破坏、基体破裂与纤维层分层多种破坏形式所占比重随着破片速度的变化而变化。Abstract: In order to explore the anti-penetration property of carbon fiber composite materials, a series of experiments need to be carried out. In experiments, 8 g cubic steel fragments penetrated carbon fiber reinforced plastics whose thickness are 5, 10 and 15 mm respectively. The velocity of fragments was obtained, the situation of target was observed, fragments and carbon fiber were recovered. The numerical simulations were carried out according to the situation of experiments. The process of fragments penetrating carbon fiber composite target was explored, and the failure mechanism of carbon fiber composite target was also explored. The relationship between the failure mode and velocity of fragments could be described as follows: when the fragment velocity exceeds the ballistic limit, the main form of target damage is fiber shear failure; when the fragment velocity is lower than the ballistic limit, the damage forms of the target include fiber shear failure, fiber tensile failure, matrix cracking and fiber layer delamination, and their proportions change with the fragment velocity.
-
弹药在储存、运输以及使用过程中,一旦遭受意外热刺激,则易发生燃烧、爆炸甚至爆轰,因此有必要对炸药的热安全性进行合理评估。在确保弹药安全的前提下提高炸药威力是目前亟待解决的问题。在混合炸药配方中添加氧化剂可改善其氧平衡,较大幅度地提升爆热。高氯酸铵(NH4ClO4,ammonium perchlorate,AP)具有密度较高、含氧量高、爆轰产物为气体、机械感度低等优点,是混合炸药配方的重要候选组分[1]。研究含有AP的浇注炸药的热响应特性对其推广应用以及新型钝感高聚物黏结炸药(plastic bonded explosives,PBX)研发具有重要意义。
国内外对含AP的含能材料热分解特性以及热安全性开展了较多研究。刘子如等[2]研究了AP-奥克托今(HMX)混合物的热分解特性,发现AP的存在使HMX的微商热重(derivative thermogravimetry,DTG)峰温大幅下降,证明AP与HMX在热分解中存在强烈的相互作用。Li等[3]研究了铝粉的加入对AP-HMX混合物热分解特性的影响,发现足量的铝粉可以消除AP与HMX之间的相互催化作用。陈中娥等[4]认为,AP热分解形成的孔隙是AP-HTPB推进剂慢速烤燃响应剧烈的原因。李苗苗等[5]对不同AP含量的HTPE推进剂进行了差示扫描量热法测试及慢速烤燃实验,发现随着AP含量的降低,推进剂的热稳定性和热安全性略有提升。Kim等[6]对含有AP的PBX炸药进行了快速和慢速烤燃实验,结果表明,在慢速烤燃中AP不参与点火,而在快速烤燃中点火由AP与硝胺类炸药共同参与。由此可见,AP在热刺激下的自热反应受到添加剂和温升速率等因素的影响,但是目前仍缺乏含AP的PBX炸药的热刺激响应特性相关研究[7]。
在炸药的热刺激响应特性实验中,温升速率对炸药的点火时间、点火位置以及装药的温度梯度等产生很大的影响。王沛等[8]对GHL炸药进行了慢速烤燃实验以及4种不同温升速率下的数值模拟,发现点火位置随着温升速率的增大从炸药内部移向炸药边缘。代晓淦等[9]进行了不同温升速率下PBX-2炸药的烤燃实验和数值模拟,发现温升速率越大,炸药中心温度在烤燃过程中的变化越小。刘静等[10]对模块装药进行了慢速烤燃实验和数值模拟,发现随着温升速率的提高,烤燃响应区域由1个中心环形区域变成关于中心对称的2个环形区域。由此可见,温升速率对炸药烤燃中的点火位置以及点火区域形状的影响较大,对获取热刺激下炸药的响应特性具有重要意义。
本研究将开展新型浇注炸药GOLA-1的慢速烤燃实验和数值模拟,分析GOLA-1炸药在不同条件下的热刺激响应特性及其在慢速烤燃中点火的主导因素,探讨温升速率对小型烤燃弹点火位置的影响。
1. 实验与数值模拟方法
1.1 小型烤燃弹实验
基于小型烤燃弹实验装置对GOLA-1炸药进行慢速烤燃实验。图1为小型烤燃弹实验测试系统示意图。该装置主要包括壳体、端盖、炸药、环形加热器、测温热电偶以及控温系统[11]。烤燃弹的外径为58 mm,长度为106 mm,壳体厚度为4 mm,药柱的直径为50 mm,长度为88 mm。药柱与端盖之间有高12 mm的空气域,用以缓解炸药受热膨胀对壳体的压力。端盖中心以及炸药上表面至炸药中心均开有直径为1 mm的小孔,用于放置K型热电偶,监测温升过程中炸药中心点的温度变化。壳体和端盖材质均为45钢,两者采用螺纹连接。装药壳体外侧包裹环形加热器,在慢速烤燃实验中温升速率为1.0 和1.5 K/min。加热器外部使用岩棉包裹进行保温,防止热量散失影响实验结果。在壳体与加热器之间放置两根K型热电偶,用来控制加热速率以及测量壳体外壁面的温度变化。
1.2 烤燃数值模拟
烤燃实验只能获得温度监测点的温度-时间关系,无法获取炸药内部的温度分布、点火位置、点火温度以及各组分的质量分数等信息。数值模拟可以较方便地获取以上信息,从而对炸药的热刺激响应特性进行较全面的评估。基于小型烤燃弹实验装置,建立如图2所示的数值模拟模型,其中端盖及壳体、空气域以及炸药分别用粉色、蓝色和红色表示。模型用六面体网格划分,网格尺寸为0.5 mm。
在直角坐标系下,热传导方程的形式为[12]
ρcp∂T∂t=λ∇2T+S (1) 式中:ρ为密度,kg·m−3;cp为比定压热容,J·kg−1·K−1;T为温度,K;λ为热导率,W·m−1·K−1;S为炸药的自热反应源项,W·m−3。
空气域与炸药之间的辐射热流qr采用P1辐射模型描述[13]
qr=13(a+σs)−Cσs∇G (2) 式中:a为吸收系数,m−1;σs为散射系数,m−1;C为线性各向异性相位函数系数;G 为入射辐射。
GOLA-1为HMX基浇注炸药,主要成分为HMX,其他成分包括Al、 AP、黏结剂(Kel-F)。组分HMX、Al、AP、Kel-F的质量比为44∶25∶20∶11,密度为1.837 g/cm3。对于炸药中的主要成分HMX,在慢速烤燃过程中,随着温度上升,HMX首先由β相转变为δ相,而后热分解生成最终的气态产物。基于过渡态理论,HMX在烤燃过程中的4步化学反应动力学模型[14]表示为:(1) β-HMX ↔ δ-HMX(反应1,单分子吸热反应);(2) β-HMX+δ-HMX ↔ δ-HMX(反应2,双分子吸热反应);(3) δ-HMX→产物(反应3,单分子吸热反应);(4) δ-HMX+产物→产物(反应4,双分子放热反应)。反应1和反应2描述β相HMX转变为δ相HMX的晶型转变过程,反应3和反应4描述最终产物生成的点火过程。各步反应速率方程为
r1=wβkThexp(TΔS1f−E1fRT)−wδkThexp(TΔS1r−E1rRT) (3) r2=wβwδkTh[exp(TΔS2f−E2fRT)−exp(TΔS2r−E2rRT)] (4) r3=wδZ3exp(−E3RT) (5) r4=wδwpZ4exp(−E4RT) (6) 式中:r1~r4为反应1~反应4的反应速率;wβ和wδ分别为β相和δ相HMX的质量分数,wp为产物的质量分数;k为玻尔兹曼常数;h为普朗克常数;R为普适气体常数;ΔS为活化熵,J∙mol−1∙K−1;E为活化能,J∙mol−1;下标f代表正反应,r代表逆反应;Z1~Z4为反应1~反应4的指前因子,s−1。
各组分质量分数的微分方程为
dwβdt=−r1−r2 (7) dwδdt=r1+r2−r3−r4 (8) dwpdt=r3+r4 (9) 单位时间内HMX热分解反应的反应源项SHMX为
SHMX=ρHMXQ1r1+ρHMXQ2r2+ρHMXQ3r3+ρHMXQ4r4 (10) 式中:ρHMX为HMX的密度,kg∙m−3;Q1~Q4为反应1~反应4的反应热,J∙kg−1。
Kel-F为含氟含氯黏结剂,分解时吸收热量,使用单步化学反应动力学模型描述其热分解过程。其热分解反应速率方程为[15]
r5=Z5exp(−Ea5RT) (11) 式中:r5、Z5、Ea5分别为Kel-F热分解反应的反应速率、指前因子和活化能。单位时间内Kel-F热分解反应的反应源项SKel-F为
SKel-F=ρKel-FQ5r5 (12) 式中:ρKel-F为Kel-F的密度,Q5为Kel-F热分解反应的反应热。
AP为氧化剂,其分解反应为放热反应,采用一步反应描述其热分解过程[16]
NH4ClO4→O2+H2O+HCl+N2 AP在240 ℃左右由斜方晶型(α-AP)转变为立方晶型(β-AP)[17]。在慢速烤燃的温升过程中,AP的温度几乎均在240 ℃以下,因此采用α-AP的热分解反应速率方程[18],即
r6=dαdt=Z6exp(−Ea6RT)αn(1−α)m (13) 式中:r6、Z6、Ea6分别为α-AP热分解反应速率、指前因子和活化能;α为AP的热分解转化率,初始阶段为零,完全分解为1;n和m为反应动力学参数,对于α-AP,n=1,m=1.38。根据文献[19],AP的热分解开始于150 ℃左右,因此在UDF用户子程序中设α在150 ℃以下为零,150 ℃以上为10−4,以启动对AP分解反应计算。设AP的密度为ρAP,则单位时间内AP热分解反应的反应源项SAP为
SAP=ρAPQ6r6 (14) 式中:Q6为AP热分解反应的反应热。
因为Al不参与点火前的热分解反应,所以不考虑铝粉的热分解反应源项,则炸药的总自热反应源项S为
S=0.44SHMX+0.11SKel-F+0.20SAP (15) 根据叠加原理,炸药的导热系数λ和比热容C分别表示为
λ=0.44λHMX+0.11λKel-F+0.20λAP (16) C=0.44CHMX+0.11CKel-F+0.20CAP (17) 表1和表2列出了炸药组分HMX、Kel-F和AP的反应动力学参数,表3列出了材料的物性参数。采用Fluent软件进行计算,设置温升边界条件为1.0和1.5 K/min,通过编写用户自定义函数UDF,将热导率、边界条件以及反应动力学参数加载到主程序中进行计算。
表 1 HMX的反应动力学参数Table 1. Reaction kinetic parameters of HMXReaction No. Z/s−1 ΔSf/
(J∙mol−1∙K−1)Ef/(kJ∙mol−1) ΔSr/ (J∙mol−1∙K−1) Er/(kJ∙mol−1) Q/(J∙g−1) 1 123.0 204 89.0 189 −25 2 −41.7 102 −75.2 865 −25 3 3.16×1016 200 −120 4 1.50×1014 173 3200 表 2 Kel-F和AP的反应动力学参数Table 2. Reaction kinetic parameters of Kel-F and APComponent Reaction No. Z/s−1 E/(kJ∙mol−1) Q/(J∙g−1) Kel-F 5 9.93×1017 272.0 −5870 AP 6 6.85×1012 146.3 297 表 3 材料物性参数Table 3. Parameters of material propertiesComponent Density/(kg∙m−3) C/(J·kg−1·K) λ/(W·m−1·K) HMX 1850 1004.26 0.5358 Kel-F 2020 1000.43 0.0527 AP 1950 1255 0.2760 Al 2719 871 1.39 Steel 8030 502.48 16.27 Air 1.225 1006.43 0.0242 2. 结果与讨论
2.1 烤燃实验
图3为不同温升速率下GOLA-1炸药点火后烤燃弹照片。图3(a)显示:点火后烤燃弹壳体膨胀变形,端盖被冲飞,烤燃弹底部断裂脱落,壳体螺纹完整;壳体表面有明显的烧灼痕迹,壳体内有少量未反应炸药。图3(b)显示了1.5 K/min温升速率下GOLA-1炸药点火前后烤燃弹的对比。结果显示,点火后烤燃弹壳体膨胀变形,端盖被冲飞,壳体底部与侧面有明显的烧灼痕迹,壳体内有少量未反应炸药。由实验现象判定,当温升速率为1.0和1.5 K/min时,GOLA-1炸药的热刺激响应等级均为燃烧,说明炸药在慢速烤燃条件下的热安全性较好。对比两种温升速率下炸药的响应程度,可以发现,温升速率为1.0 K/min时炸药的慢速烤燃响应更剧烈。
2.2 数值模拟
基于实验进行慢速烤燃数值计算,在1.0和1.5 K/min温升速率下实验与数值模拟得到的炸药中心测温点温度-时间曲线如图4和图5所示。可以看出,实验与数值模拟得到的温升曲线基本重合,点火温度和点火时间也基本吻合,说明数值模型及参数比较合理。图4显示,9900 s时实验测得的温升曲线的温升放缓,而后恢复正常温升速率。分析认为,此阶段的温升放缓是由于HMX在由β相转变为δ相的过程中吸收热量,温度约为451 K。在慢速烤燃实验中,炸药于11946 s时发生点火反应,点火时药柱中心点温度为488 K;计算得到的点火时间为11790 s,与实验结果的相对偏差为1.3%。图5显示:HMX相变造成的温升速率放缓出现在约6800 s,温度约为451 K;实验结果显示,炸药在约8207 s发生点火反应,点火时刻药柱中心温度约为480 K;计算得到的点火时间为8070 s,与实验结果的相对偏差为1.7%。由此可见,温升速率对点火时间的影响较大,随着温升速率的上升,点火时间缩短,然而温升速率对HMX的相变温度以及GOLA-1炸药的点火温度没有影响。
1.0和1.5 K/min温升速率下不同时刻炸药的温度分布如图6和图7所示。相比于炸药,壳体的比热容小而热导率大,因此壳体温度高于炸药温度,热量由壳体传递到装药内部。装药的部分热量散失至烤燃弹上部的空气域中,造成装药上部温度较低。图6(d)和图7(d)为点火时烤燃弹的温度分布,可见,点火区域位于装药底部的环形区域。这是由于药柱底部受到侧面壳体与底面的共同加热,使得装药底部温度相对较高。温度的升高使得炸药的自热分解反应速率上升,造成热量进一步累积,如此循环往复,致使点火最终发生在装药底部的环形区域。对比两者发现,在1.5 K/min的温升速率下装药的温度梯度较大,点火位置相对靠外。分析认为,这是由于炸药的导热系数远低于钢壳体,而温升速率相对较大,热量来不及传递到炸药内部造成的。
1.0 K/min温升速率下不同位置的温升曲线和质量分数曲线如图8、图9和图10所示。从图8可以看出,HMX在约9300 s发生相变,相变温度与中心点温度相同,均为451 K。对比中心测温点与点火点的温升曲线可以发现,点火点处HMX相变所造成的温度平台较小,且持续时间较短。点火点处的热量累积速率大于热量传递速率,温度上升较快,相同时刻点火点处的温度高于其他位置,从而使HMX相变时间早且速率快。临近点火时,由于HMX热分解反应释放大量的热量,造成点火点的温度急剧升高,点火温度约为517 K。图9和图10给出了中心测温点和点火点的HMX质量分数变化曲线。在点火时刻,β-HMX已经全部转变为δ-HMX,相变过程基本完成。点火时中心测温点仅有微量的HMX热分解产物,而点火位置则有质量分数为0.2的产物。说明在点火位置,第4步强放热反应已经开始进行,释放大量热量,导致温度急剧上升,主导了点火的发生。
1.5 K/min温升速率下不同位置的温升曲线和质量分数曲线如图11、图12和图13所示。图11中的点火点温度-时间曲线特征与图8基本相同,点火温度约为518 K。图12和图13显示:相比于温升速率为1.0 K/min的情况,点火时刻中心测温点处β-HMX的质量分数为0.04,δ-HMX的质量分数为0.96,表明部分β-HMX尚未转变为δ-HMX,第4步强放热反应尚未发生;而点火位置的HMX已完全相变,并且HMX分解的第4步强放热反应已经开始进行,有质量分数为0.19的产物生成,主导了点火的发生。
图14给出了温升速率为1.5 K/min时中心测温点处AP的转化率-时间曲线。转化率α与分解速率dα/dt均随着温度的升高而增大。在点火时刻,AP的转化率为1.06×10−4,表明此时中心点的AP基本不发生分解。7000 s左右,AP的分解速率放缓,而后快速上升,这是HMX在晶型转变过程中吸收热量而导致温升放缓造成的。图15为点火点处AP的转化率-时间曲线。可以看出:转化率在临近点火时刻迅速上升;点火时刻AP的转化率为2.26×10−4,相较于中心点处的转化率高约1倍。不论是中心测温点还是点火点,发生分解的AP的量均很少,表明在慢速烤燃过程中AP基本不参与反应,点火反应由HMX主导。这与Kim等[6]得到的含硝胺类及AP浇注炸药慢速烤燃实验结论一致。
2.3 温升速率对点火位置的影响
将实验与数值模拟相结合,分析温升速率对炸药热刺激响应特性的影响。在1.0和1.5 K/min两种温升速率的基础上,设定不同的温升速率,对GOLA-1炸药进行计算。当温升速率为0.5和0.4 K/min时,炸药点火位置的变化如图16所示。0.5 K/min时,点火位置仍为装药底部的环形区域,相比1.0和1.5 K/min的情况,环形区域缩小,并向装药的上部及中心移动。这是由于当温升速率下降时,炸药整体受热比较均匀,装药内部的温度梯度减小,相比于温升速率为1.0 K/min的情况,相同加热温度下炸药内部高温区域的体积较大,炸药的自热反应速率上升,炸药内部积累热量,导致点火发生在装药底部的环形区域。当温升速率降至0.4 K/min时,点火位置的形状发生改变,由环形变为扁的椭球形,并且位置移动至药柱下部中心。
为进一步分析温升速率对GOLA-1点火位置的影响,以0.05 K/min为间隔,进行0.1~0.4 K/min 7个不同温升速率的烤燃数值模拟。设温升速率为β,以壳体外侧底面的中心位置为坐标零点,沿烤燃弹中心轴线从壳体底面指向端盖为正方向,点火位置x为温度最高区域的中心,单位为m。温升速率对点火位置的影响如表4所示。随着温升速率的下降,点火位置向药柱上部移动。根据表4绘制点火位置随温升速率的散点图,如图17所示。观察图17发现,温升速率β与点火位置x之间存在线性关系,拟合关系式为
x=−0.05279β+0.04171 ,相关系数 R2=0.966。表 4 温升速率与点火位置的关系Table 4. Heating rate vs. ignition positionNo. β/(K·min−1) x/m 1 0.40 0.0212 2 0.35 0.0220 3 0.30 0.0273 4 0.25 0.0287 5 0.20 0.0298 6 0.15 0.0334 7 0.10 0.0374 3. 结 论
(1) 对新型HMX基含AP浇注炸药GOLA-1进行慢速烤燃实验,结果显示,当温升速率为1.0和1.5 K/min时炸药的烤燃响应等级均为燃烧,并且1.0 K/min下壳体的破坏更严重,响应更剧烈。从热刺激响应角度来看,温升速率越低,GOLA-1炸药的响应越剧烈。
(2) 开展了GOLA-1炸药烤燃数值模拟,得到的温升曲线与实验结果吻合较好。当温升速率为1.0 K/min时,GOLA-1炸药的点火时间与实验结果的相对偏差为1.3%,1.5 K/min时为1.7%,证明建立的烤燃模型及参数较合理。数值模拟结果显示,点火反应主要由HMX的第4步强放热反应主导,AP在慢速烤燃过程中基本不参与点火反应。
(3) 对不同温升速率下GOLA-1炸药的烤燃过程进行了数值模拟,结果表明:当温升速率为0.5 K/min时,点火仍发生在装药底部的环形区域,并向装药中心及上部移动;当温升速率下降至0.4 K/min时,点火区域的形状发生变化,由环形变为椭球形;随着温升速率继续下降,点火位置沿着药柱中心轴线向上部移动;当温升速率处在0.1~0.4 K/min之间时,温升速率与点火位置之间存在线性关系。
-
表 1 碳纤维复合材料的力学性能参数
Table 1. Mechanical properties of carbon fiber composites plate
σxt/
MPaE1t/
GPaσyt/
MPaE2t/
GPaσxc/
MPaE1c/
GPaσyc/
MPaE2c/
GPaτ12/
MPaG12/
GPaτJ/
MPaγ12 1755 138 48 8.36 1248 128 214 8.48 113 4.51 21.5 0.27 表 2 破片侵彻靶板的损伤区域统计
Table 2. Statistic of damage area of fragment penetrating target
Target thickness/mm Failure pattern Impact velocity/
(m·s–1)Damage area of the
projectile surface/
(mm×mm)Damage area of back
elastic surface/
(mm×mm)5 Local penetration 328 10.92×9.60 11.55×12.24 Complete penetration 340 13.76×10.50 23.10×25.10 10 Local penetration 398 10.63×10.21 29.80×31.42 Complete penetration 412 14.56×10.36 30.30×17.22 15 Local penetration 443 16.26×8.35 Complete penetration 453 15.49×11.22 33.23×24.05 表 3 钢破片材料参数
Table 3. Material parameters of steel fragment
ρ/(g·cm–3) G/GPa A/MPa B/MPa n c M TM/K TR/K 7.83 0.77 792 510 0.26 0.014 1.03 1793 294 表 4 试验和数值模拟得到的弹道极限误差
Table 4. Ballistic limit error of test and numerical simulation
Target thickness/mm Ballistic limit velocity/(m·s–1) Deviation/% Test Simulation 5 334 307.76 7.86 10 405 394.78 2.52 15 448 437.55 2.33 -
[1] ANUSE V S, SHANKAR K, VELMURUGAN R, et al. LVI and CAI analysis of woven carbon fiber reinforced composite laminates with different stacking sequence [J]. Key Engineering Materials, 2023, 969: 93–100. doi: 10.4028/p-n2McuG [2] 常新龙, 孙超凡, 齐重阳, 等. 集成于碳纤维复合材料的fsFBG响应特性研究 [J]. 兵器装备工程学报, 2024, 45(1): 208–216. doi: 10.11809/bqzbgcxb2024.01.027CHANG X L, SUN C F, QI C Y, et al. Response characteristics of fsFBG sensor based on carbon fiber composite laminate [J]. Journal of Ordnance Equipment Engineering, 2024, 45(1): 208–216. doi: 10.11809/bqzbgcxb2024.01.027 [3] HOSUR M V, ALEXANDER J, JEELANI S, et al. High strain compression response of affordable woven carbon/epoxy composites [J]. Journal of Reinforced Plastics and Composites, 2003, 22(3): 271–296. doi: 10.1177/0731684403022003844 [4] ABRATE S. Impact on laminated composites: recent advances [J]. Applied Mechanics Reviews, 1994, 47(11): 517–544. doi: 10.1115/1.3111065 [5] 金子明, 隋金玲, 张菡英, 等. 纤维增强复合防弹板研究进展及抗弹性能研究 [J]. 玻璃钢, 2001(1): 1–6.JIN Z M, SUI J L, ZHANG H Y, et al. Research progress and anti-elastic properties of fiber reinforced composite bulletproof panels [J]. Fiber Reinforced Plastics, 2001(1): 1–6. [6] MORYE S S, HINE P J, DUCKETT R A, et al. Modelling of the energy absorption by polymer composites upon ballistic impact [J]. Composites Science and Technology, 2000, 60(14): 2631–2642. doi: 10.1016/S0266-3538(00)00139-1 [7] LIU P F, ZHENG J Y. Recent developments on damage modeling and finite element analysis for composite laminates: a review [J]. Materials & Design, 2010, 31(8): 3825–3834. doi: 10.1016/j.matdes.2010.03.031 [8] ORIFICI A C, HERSZBERG I, THOMSON R S. Review of methodologies for composite material modelling incorporating failure [J]. Composite Structures, 2008, 86(1/2/3): 194–210. doi: 10.1016/j.compstruct.2008.03.007 [9] BØRVIK T, LANGSETH M, HOPPERSTAD O S, et al. Perforation of 12 mm thick steel plates by 20 mm diameter projectiles with flat, hemispherical and conical noses: part Ⅰ: experimental study [J]. International Journal of Impact Engineering, 2002, 27(1): 19–35. doi: 10.1016/S0734-743X(01)00034-3 [10] BØRVIK T, HOPPERSTAD O S, BERSTAD T, et al. Perforation of 12 mm thick steel plates by 20 mm diameter projectiles with flat, hemispherical and conical noses. part Ⅱ: numerical simulations [J]. International Journal of Impact Engineering, 2002, 27(1): 37–64. doi: 10.1016/S0734-743X(01)00035-5 [11] IQBAL M A, GUPTA N K. Ballistic limit of single and layered aluminium plates [J]. Strain, 2011, 47(Suppl 1): e205–e219. doi: 10.1111/j.1475-1305.2008.00601.x [12] GUPTA N K, IQBAL M A, SEKHON G S. Experimental and numerical studies on the behavior of thin aluminum plates subjected to impact by blunt- and hemispherical-nosed projectiles [J]. International Journal of Impact Engineering, 2006, 32(12): 1921–1944. doi: 10.1016/j.ijimpeng.2005.06.007 [13] 刘璐璐, 宣海军, 洪伟荣, 等. 锻纹机织复合材料高速冲击试验研究及理论分析 [J]. 中国科技论文, 2016, 11(10): 1169–1174. doi: 10.3969/j.issn.2095-2783.2016.10.020LIU L L, XUAN H J, HONG W R, et al. Experimental study and theoretical analysis on high velocity impact of satin woven carbon fiber reinforced composites [J]. China Sciencepaper, 2016, 11(10): 1169–1174. doi: 10.3969/j.issn.2095-2783.2016.10.020 [14] 邓云飞, 袁家俊, 徐建新. 半球形头弹不同角度冲击下编织复合材料板的侵彻特性 [J]. 复合材料学报, 2018, 35(4): 843–849. doi: 10.13801/j.cnki.fhclxb.20170704.002DENG Y F, YUAN J J, XU J X. Penetration characteristics of woven composite laminates impacted by hemispherical-nosed projectiles at different angles [J]. Acta Materiae Compositae Sinica, 2018, 35(4): 843–849. doi: 10.13801/j.cnki.fhclxb.20170704.002 [15] 邓云飞, 袁家俊, 徐建新. 刚性卵形头弹撞击角度对编织复合材料层合板侵彻特性的影响 [J]. 复合材料学报, 2018, 35(8): 2039–2045. doi: 10.13801/j.cnki.fhclxb.20170829.004DENG Y F, YUAN J J, XU J X. Influence of the impact angle of rigid ogival-nosed projectiles on the penetration characteristics of woven composite laminates [J]. Acta Materiae Compositae Sinica, 2018, 35(8): 2039–2045. doi: 10.13801/j.cnki.fhclxb.20170829.004 [16] 谢文波, 张伟, 姜雄文. 钢球斜侵彻碳纤维复合材料板的实验研究 [J]. 爆炸与冲击, 2018, 38(3): 647–653. doi: 10.11883/bzycj-2016-0289XIE W B, ZHANG W, JIANG X W. Oblique penetration on CFRPs by steel sphere [J]. Explosion and Shock Waves, 2018, 38(3): 647–653. doi: 10.11883/bzycj-2016-0289 [17] 彭捷, 张伟岐, 田锐, 等. 碳纤维层合板抗球形弹冲击动态响应特性试验研究 [J]. 复合材料科学与工程, 2020(6): 18–24, 56. doi: 10.3969/j.issn.1003-0999.2020.06.003PENG J, ZHANG W Q, TIAN R, et al. Experimental study on the dynamic response of carbon fiber laminates impacted by spherical projectile [J]. Composites Science and Engineering, 2020(6): 18–24, 56. doi: 10.3969/j.issn.1003-0999.2020.06.003 [18] 赵稼祥. 日本东丽公司碳纤维和石墨纤维的新进展 [J]. 宇航材料工艺, 1982(5): 37.ZHAO J X. New progress of carbon fiber and graphite fiber of toray corporation, japan [J]. Aerospace Materials & Technology, 1982(5): 37. [19] QI S B, HUANG G Y, ZHI X D, et al. External blast flow field evolution and response mechanism of single-layer reticulated dome structure [J]. Defence Technology, 2023, 24: 241–253. doi: 10.1016/j.dt.2022.04.003 -