Numerical Simulation of Anti-Penetration Performance of Body Armor against Small Tungsten Sphere
-
摘要: 为研究防弹衣抗小钨球侵彻的性能,结合试验,利用有限元分析软件LS-DYNA建立了小钨球侵彻防弹衣的数值模型。在此基础上,对侵彻过程进行了数值模拟,分析了防弹衣的破坏机理,并探讨了凯夫拉(Kevlar)与超高分子量聚乙烯(Ultra-high molecular weight polyethylene, UHMWPE)混杂配比对防弹衣抗侵彻性能的影响。研究结果表明:在小钨球侵彻作用下,防弹衣迎弹面主要发生纤维剪切破坏,背弹面主要发生纤维拉伸断裂破坏,并伴随着一定的层间分层破坏。随着着靶速度的提高,纤维的拉伸及分层破坏程度降低;与单一Kevlar制作的防弹衣相比,采用面板Kevlar、背板UHMWPE混杂结构的防弹衣抗侵彻性能更好。当Kevlar/UHMWPE的体积配比分别为1∶1、1∶2和1∶4时,防弹衣的抗侵彻性能分别提高3.7%、5.3%和4.4%,质量分别减少14.1%、18.8%和22.5%。综合考虑防弹衣的抗侵彻性能和重量,采用Kevlar/UHMWPE混杂配比为1∶2的防弹衣结构最佳;在弹道极限附近,采用Kevlar/UHMWPE混杂结构的防弹衣的吸能效果优于单一Kevlar结构,且随着着靶速度的提高,两者的吸能差异逐渐减小。研究结果对防护装备的优化设计具有一定的参考价值。Abstract: In order to study the anti-penetration performance of body armor against small tungsten sphere, combined with the experiment, the numerical model of small tungsten sphere penetrating into body armor was established by using the finite element analysis software LS-DYNA. On this basis, the penetration process was numerically simulated, the damage mechanism of body armor was analyzed, and the influence of Kevlar and Ultra-high molecular weight polyethylene (UHMWPE) hybrid ratio on the anti-penetration performance of body armor was discussed. The investigation results show that the shear failure of fiber mainly occurs on the impact face of body armor, and the tensile fracture of fiber mainly occurs on the back of body armor, accompanied by a certain delamination failure under the penetration effect of small tungsten sphere. The tensile failure and delamination failure degree of fiber decrease with the increase of the impact velocity of small tungsten sphere. Compared with the body armor made of Kevlar, the body armor with hybrid structure of Kevlar placed on the impact face and UHMWPE placed on the back has better anti-penetration performance. When the hybrid ratio of Kevlar and UHMWPE is 1∶1, 1∶2 and 1∶4, the anti-penetration performance of body armor is improved by 3.7%, 5.3% and 4.4% respectively, and the mass of body armor is reduced by 14.1%, 18.8% and 22.5% respectively. In the comprehensive consideration of anti-penetration performance and weight of body armor, the Kevlar/UHMWPE structure whose fiber hybrid ratio is 1∶2 is the best. When the impact velocity is near the ballistic limit, Kevlar/UHMWPE hybrid structure has better energy absorption effect than single Kevlar structure. However, with the increase of the impact velocity, the difference of energy absorption between them decreases gradually. The research results have a certain reference value for the optimization design of protective equipment.
-
Key words:
- body armor /
- small tungsten sphere /
- anti-penetration performance /
- hybrid composites
-
负泊松比蜂窝结构又称拉胀结构,因具有许多常规结构不具备的力学特性[1],而成为研究热点。蜂窝材料具有较高的相对刚度、强度和高效的能量吸收能力,在抗剪切、抗屈曲、提高硬度以及抗疲劳等方面拥有独特的优越性[2-3],在一些应用领域中发挥着关键作用,如汽车、航空、军事、医学领域[4]。多孔结构的力学性能主要取决于细观上的拓扑结构。近年来,通过改变细观结构,人们发现负泊松比结构具有很多特殊优势,因而被广泛应用[5]。
马芳武等[6]研究了一种内凹三角形负泊松比结构,通过改变内凹角度,分析了冲击端和固定端的平台应力和能量吸收能力,并与内凹六边形进行了对比。Zhang等[7]分析了内凹六边形蜂窝在两个正交方向上的后继屈服拉伸行为,同时考虑结构的塑性影响以及孔壁的非线性行为分析模型,提出了单胞结构的塑性铰变形机制,得到了单胞结构的应力-应变曲线。Li等[8-9]对内凹蜂窝结构进行分级、强化,并将正弦曲线引入内凹蜂窝结构,得到了新的改进模型,进而分析了结构的泊松比和能量吸收变化。邓小林等[10]研究了全参数化的正弦曲线蜂窝结构,以不同振幅、不同厚度建立模型,研究了蜂窝结构在不同冲击速度下的动力响应,发现正弦曲线蜂窝较常规六边形蜂窝有更好的能量吸收效果。崔世堂等[11]利用有限元模拟方法研究了负泊松比蜂窝结构面内冲击动力学特性,发现平台应力和结构的比吸能随冲击速度的增大而增高,随胞元扩展角的增大而降低。陈鹏等[12]研究了具有零泊松比特征的半凹角蜂窝结构,并将其与正六边形蜂窝和内凹负泊松比蜂窝在面内冲击荷载作用下的抗冲击性能进行对比分析,数值结果表明,半凹角蜂窝的抗冲击性能介于正六边形蜂窝和内凹蜂窝之间。Hu等[13]通过理论分析和数值模拟,研究了内凹角度和壁长对内凹负泊松比蜂窝在大变形下的单轴动态冲击性能的影响,推导出冲击过程中平均冲击应力的经验公式。Zhang等[14]通过有限元模拟,研究了内凹蜂窝x方向的平面内动态冲击行为,发现内凹蜂窝的面内动态性能不仅与冲击速度和边缘厚度有关,还受蜂窝壁角的影响。Li等[15]通过单轴和双轴压缩模拟以及理论分析,研究了正六边形蜂窝结构的面内压缩动态力学性能,分析了双轴压缩的变形模式,结果表明:相比单轴冲击,双轴冲击下在x和y方向的真实应力增强,能量吸收能力也得到了提高,且完全致密化阶段比单轴压缩阶段更平滑。此外,Li等[16]研究了六边形、内凹、混合3种蜂窝模型在单、双轴冲击下的面内动态力学性能,结果表明:正交双轴冲击下,六边形蜂窝表现出3种变形模式,内凹和混合型蜂窝没有明显的过渡模式,由于负泊松比效应的影响,内凹蜂窝具有较差的耗能能力。
值得注意的是,自然界中的蜂窝结构和人造蜂窝结构在细观上总存在一定的缺陷,从而引起结构的不规则性,力学性能也会发生一定的变化。Ajdari等[17]通过数值模拟研究了正六边形和不规则二维蜂窝的平面内动态冲击问题,分析了孔壁缺失和空间扰动形成的结构微观不规则性对力学性能的影响。Alkhader等[18]用函数定义六边形蜂窝、随机Voronoi泡沫以及正方形和三角形拓扑结构等多种二维拓扑结构的不规则程度,以研究其单轴压缩响应,结果表明,相对于以弯曲为主的结构,以拉伸为主的结构有表现出灾难性屈服后软化反应的趋势,而不规则性则会导致更多的弯曲现象。Liu等[19]对内凹蜂窝材料的面内动态冲击过程进行了数值模拟,并在此基础上定义了内凹蜂窝结构的不规则性,分析发现,在准静态下不规则的内凹蜂窝比规则的正六边形蜂窝能吸收更多的能量,但这种情况在高速撞击下逆转。Zheng等[20]通过数值模拟研究了坐标扰动和Voronoi随机模型两种不规则模型与正六边形蜂窝在不同冲击速度下的变形模式和平台冲击力,得到不规则性结构更具复杂性的结论。Zhu等[21]研究了孔的不规则性对二维随机泡沫弹性性能的影响,构造了不规则度不同的周期性随机结构,并通过数值模拟确定了其有效弹性性能,结果表明,二维随机泡沫体形状越不规则,有效弹性模量和剪切模量越大,在一定的压比相对密度下,体积模量越小。
综上所述,实际中蜂窝结构往往是不规则的,且易受双轴冲击载荷作用。而关于不规则结构在双轴冲击下的研究较少,为此本工作将针对不规则内凹负泊松比蜂窝结构在双轴冲击下的面内冲击响应,分析规则度和冲击速度对结构变形影响的规律。
1. 计算模型
1.1 不规则内凹蜂窝结构有限元模型的建立
采用如图1所示的节点扰动方法来建立不规则内凹蜂窝的有限元模型。
如图1所示,将规则的内凹六边形蜂窝结构的每个节点按照式(1)中的方法进行随机扰动
{yi=y0+Δysinθmxi=x0+Δxcosθm (1) 式中:
θm 为角度随机值(0∘⩽θm⩽360∘ ),x0 、y0 为节点的原始坐标,xi 、yi 为扰动后的节点坐标。为了保证随机扰动之后模型的棱壁不会重叠,需要对Δx 、Δy 进行如式(2)的限制{μ=Δx2+Δy20⩽μ⩽μm (2) 式中:
μ 为节点扰动的随机长度;μm 为扰动的最大长度,0⩽μm⩽l1/2 。内凹负泊松比蜂窝结构的不规则度可以定义为K=2μml1 (3) 式中:
l1 为规则蜂窝结构的最短棱壁长度。假设蜂窝结构所有棱壁的厚度均相同,则可通过改变棱壁的厚度来调节蜂窝结构的相对密度。本研究采用15%的相对密度进行分析,图2显示了部分模型。
图2中内凹蜂窝结构的相对密度
Δρ 可以表示为Δρ=ρ∗ρs=N∑i=1li×tL1×L2 (4) 此外,对于规则的内凹负泊松比蜂窝,其相对密度
Δρr 也可以表示为Δρr=12tl1(l2/l1+2)cosα(l2/l1+sinα) (5) 式中:
ρ∗ 为模型的密度,ρs 为基体材料的密度,li 为各个孔壁的长度,t 为孔壁的厚度,N 为孔壁的数量,L1 、L2 为整个蜂窝结构的长度和宽度,l2为规则蜂窝结构的最长棱壁长度。1.2 有限元模型
采用ABAQUS/EXPLICIT软件进行分析。模型的边界条件设置:在两个正交方向上,将模型置于两块刚性板之间、底部刚性板之上,底部和左端的刚性板作为固定端, 顶部和右端作为冲击端,冲击速度为3~100 m/s,同时约束内凹蜂窝结构的面内自由度,如图3所示。建立的内凹蜂窝结构的主要参数为L1 = 129.9 mm, L2 = 120.0 mm, l1 = 5 mm, l2 = 10 mm,
θ = 60°。由于蜂窝铝具有高强度和高刚度的良好力学性能,本研究采用铝合金作为基体材料,主要参数为:密度ρ = 2700 kg/m3,弹性模量E = 72 GPa,泊松比为0.33,屈服强度σy = 103 MPa,并采用线性强化模型,图4为结构基体材料的本构关系,其中Et为切线模量,σs 为线性强化模型的屈服强度。蜂窝细胞数量为15 × 15,可保证材料不受尺寸效应的影响。所有单元均采用4节点壳单元进行网格划分,网格单元尺寸为0.5 mm,节点数为28660,网格数为 19 540,建立无摩擦和通用接触。2. 结果与讨论
2.1 双轴冲击下的变形模式
为了对双轴冲击条件进行分类,采用与双轴冲击有关的参数
λ ,表达式为λ=vx/vy ,其中vx 、vy 分别为x 、y 方向的冲击速度。这里只讨论λ=1 的情况,显然λ=1 时为等双轴冲击。首先研究了不同规则度(K = 0, 0.6, 1.0)的内凹负泊松比在不同冲击速度(6、50和100 m/s)下的变形模式。图5、图6和图7给出了内凹蜂窝结构的结构变形情况。需要说明的是,为更好地展示变形结果,每隔约10%的应变截取一张变形模态图,同时为了清晰、规律地显示图像,所有图形都设置了相同的大小。从图5可以看出,对于规则的内凹蜂窝,在等双轴低速冲击过程中,结构首先在交叉处棱壁堆积,从而使内部先形成四边形,结构整体的变形在近端和远端都较均匀。这与文献[15]中内凹蜂窝的变形是一致的,也验证了本模型的有效性。持续的压缩使孔壁进一步堆积形成局部致密化,结构的致密过程主要是局部致密。由于负泊松比效应的影响,材料在一个方向受压时,其另一个正交方向会出现颈缩。因此,在双向冲击受压的情况下,结构会更早进入完全密实阶段。随着冲击速度的增大,结构从冲击端(上部和右端)开始密实,而固定端几乎没有变形。随着应变增加,致密向固定端传递,直至完全进入密实化。从图5中第2行和第3行图像还可以看出,随着冲击速度的增大,蜂窝结构的下端会产生部分“翘起”现象,这是由于负泊松比效应的影响会导致结构颈缩,且结构与固定端端部是无绑定约束,从而造成这类现象。
与规则蜂窝不同的是,不规则蜂窝结构在低速冲击下,其内部不会形成较为规则的四边形。这是由于不规则度的存在使结构棱壁处的堆叠也变得不规则。此外,从图5~图7中
ε=0.5 列可以看出,由于不规则度的引入,结构的变形模式由局部密实转变为整体密实,从而使内凹蜂窝结构在相同压缩程度下,密实化程度明显降低。在高速冲击下(v = 100 m/s),从图5~图7中可以看出,不规则程度越高,冲击端的致密程度越大。这是因为高速冲击下,结构在冲击端的密实主要是棱壁的弯曲折叠过程,随着不规则度的增加,棱壁的弯曲折叠受到的约束增大,向固定端传递的速度也会降低,所以不规则蜂窝结构的密实过程会更长,而在冲击端密实程度也会更高。此外,从图6、图7中ε=0.6 、v = 100 m/s对应的变形情况可以看出,固定端还有尚未变形进入密实的孔,说明不规则蜂窝结构具有较长的平台阶段,能够承受更大的压缩变形。2.2 不同冲击速度下的应力-应变曲线
图8和图9给出了蜂窝结构在双轴冲击下两个正交方向冲击端的名义应力-应变曲线,其中名义应力
σ 通过冲击端的反力除以对应截面的原始面积获得,名义应变ε 通过冲击位移除以对应的原长获得。从图中可以看出,内凹蜂窝结构在不同方向上的σ-ε 曲线均表现出典型多孔材料在受压时所具有的弹性阶段、平台阶段和密实阶段3部分。从图8和图9中v = 6 m/s时的曲线可以看出,对于K = 0时的应力-应变曲线,在应变接近0.4处,结构变形的平台阶段均出现一个上升的阶梯,并且x方向最明显。结合2.1节关于变形模态的分析,认为这主要是由于在等低速双轴冲击下内凹蜂窝结构变形主要经历两种棱壁堆叠过程,即堆叠形成四边形以及四边形的进一步弯曲堆叠。由于第1步的堆叠,棱壁基本不会屈曲,主要是旋转折叠,因此这一平台阶段的应力水平较低;第2步的堆叠主要是棱壁的屈曲折叠,所以此阶段的应力水平较高。从图8和图9中也可以看出,K = 0时,结构会更早进入密实化阶段,而不规则度的引入使结构拥有较长的平台阶段,密实化阶段出现滞后现象,此现象与2.1节中变形模态的分析结果是一致的。随着冲击速度的增大,平台阶段的应力升高,说明结构的能量吸收能力随着冲击速度的增大而增强。
蜂窝结构的平台应力一般表示为
σp=1εd−ε0∫εdε0σ(ε)dε (6) 式中:
σp 为平台应力;ε0 为对应初始应力峰值的名义应变;εd 为锁定应变,为蜂窝结构密实化阶段所对应的应变;σ(ε) 为名义应力-应变曲线。图10给出了不规则度不同的内凹蜂窝结构在两个正交方向上不同冲击速度下的平台应力变化趋势。从图10可以看到:随着冲击速度的增大,两个方向上的平台应力值都会上升;对于K = 0的规则蜂窝结构,其在两个方向上的平台应力相差较大,这是结构的各向异性所导致的。引入不规则度时,在高速冲击下两个方向上的平台应力变化大小及趋势都较接近,说明结构的各向异性降低,这一点从2.1节内凹蜂窝结构的变形模态中也可以看出。
2.3 内凹蜂窝结构的能量吸收能力
在动态冲击过程中,能量主要由材料的塑性变形消耗。采用比塑性耗散能表征单位质量的能量吸收能力,表达式为
W=EPEDM (7) 式中:
EPED 为塑性耗散能,可以从有限元分析软件中直接获得;M为结构的质量。图11给出了内凹蜂窝结构在6、50和100 m/s 3种不同冲击速度下的比塑性能量耗散与
y 方向冲击应变的关系。从图11中可以看出,当应变较低时,比塑性耗散能W 上升较缓慢,且所有曲线基本重合。这表明在早期,不规则度对内凹蜂窝结构的影响较小。随着压缩程度的增加,W 增加的速率变大,且K = 0时,W 增加得最快,表明结构开始进入密实阶段,这是由结构的负泊松比效应引起的。对于不规则蜂窝结构,曲线上升得较缓慢,表明不规则度的引入使结构的平台阶段延长,结构具有更强的能量吸收能力。3. 结 论
采用有限元方法研究了具有不同不规则度内凹负泊松比结构的面内双轴冲击响应,得到了以下结论。
(1)内凹蜂窝结构的变形受冲击速度的影响。随着冲击速度的提高,蜂窝结构的变形逐渐转向逐层致密,受结构负泊松比效应的影响,在等高速双轴压缩时,结构的固定端会有局部“翘起”现象。此外,由于不规则度的引入,在低速冲击下,结构的密实化过程从局部致密转变为整体致密,从而导致在相同的压缩程度下,结构的密实化程度降低。
(2)随着冲击速度的增大,平台阶段的应力上升,能量吸收能力更强,比塑性耗散能也上升。不规则度的引入延长了平台阶段,降低了结构的各向异性程度,从而提高了结构的能量吸收能力。
-
表 1 试验数据
Table 1. Experimental data
No. Impact velocity/(m·s−1) Residual velocity/(m·s−1) Perforation diameter/mm Perforation 1 694.8 0 2.90 No 2 714.8 0 2.92 No 3 725.3 0 2.92 No 4 732.7 15.2 2.96 Yes 5 748.4 118.3 2.98 Yes 6 811.5 189.6 2.98 Yes 表 2 钨球材料模型参数
Table 2. Material model parameters of tungsten sphere
ρ/(g·cm−3) E/GPa μ σ0/MPa Et/MPa β C P εf 13.7 367 0.303 1 506 792 1 3.9 6 1.2 表 3 Kevlar材料模型参数
Table 3. Material model parameters of Kevlar
ρ/(g·cm−3) E1/GPa E2/GPa E3/GPa μ12 μ13 μ23 G12/GPa G13/GPa G23/GPa 1.35 21 21 4.6 0.31 0.14 0.14 1.2 1.2 1.2 Kf/GPa SC/GPa XT/GPa YT/GPa YC/GPa α SN/GPa S13/GPa S23/GPa 2 0.25 1.8 1.8 1.4 0.5 0.55 0.55 0.55 表 4 仿真值与试验值的比较
Table 4. Comparison between simulation results and experimental results
No. Impact velocity/(m·s−1) Residual velocity/(m·s−1) Perforation diameter/mm Ballistic limit/(m·s−1) Error/% Experiment Simulation Experiment Simulation Experiment Simulation 1 694.8 0 0 2.90 2.83 729 732 0.41 2 714.8 0 0 2.92 2.86 3 725.3 0 0 2.92 2.81 4 732.7 15.2 19.2 2.96 2.85 5 748.4 118.3 108.0 2.98 2.87 6 811.5 189.6 208.3 2.98 2.86 ρ/(g·cm−3) E1/GPa E2/GPa E3/GPa μ12 μ13 μ23 G12/GPa G13/GPa G23/GPa 0.97 30.7 30.7 1.97 0.008 0.044 0.044 0.73 0.67 0.67 Kf/GPa SC/GPa XT/GPa YT/GPa YC/GPa α SN/GPa S13/GPa S23/GPa 2.2 0.36 3.0 3.0 2.5 0.5 0.95 0.95 0.95 表 6 小钨球侵彻不同结构防弹衣弹道极限的仿真结果
Table 6. Simulation results of ballistic limit of small tungsten sphere penetrating into body armor with different structures
Type of tungsten sphere Thickness of body armor/mm Structure of body armor Ballistic limit/(m·s−1) 0.16 g, ∅2.8 mm 9 K1U1 759.0 K1U2 771.0 K1U4 764.0 表 7 小钨球侵彻不同结构防弹衣的仿真结果
Table 7. Simulation results of small tungsten sphere penetrating into body armor with different structures
Type of tungsten sphere Thickness of body armor/mm Impact velocity/(m·s−1) Residual velocity/(m·s−1) K1U0 K1U1 K1U2 K1U4 0.16 g, ∅2.8 mm 9 750.0 108.9 0 0 0 775.0 162.1 114.6 60.2 80.1 800.0 199.5 187.6 164.6 177.0 825.0 241.6 229.9 206.5 219.6 850.0 282.9 274.3 254.5 265.4 -
[1] FLANAGAN M P, ZIKRY M A, WALL J W, et al. An experimental investigation of high velocity impact and penetration failure modes in textile composites [J]. Journal of Composite Materials, 1999, 33(12): 1080–1103. doi: 10.1177/002199839903301202 [2] SIKARWAR R S, VELMURUGAN R, GUPTA N K. Influence of fiber orientation and thickness on the response of glass/epoxy composites subjected to impact loading [J]. Composites Part B: Engineering, 2014, 60: 627–636. doi: 10.1016/j.compositesb.2013.12.023 [3] NAIK N K, SHRIRAO P, REDDY B C K. Ballistic impact behaviour of woven fabric composites: parametric studies [J]. Materials Science and Engineering: A, 2005, 412(1/2): 104–116. [4] MAJZOOBI G H, ZAHERIB F M. Numerical and experimental study of ballistic response of kevlar fabric and Kevlar/Epoxy composite [J]. International Journal of Engineering, Transactions B: Applications, 2017, 30(5): 791–799. [5] JORDAN J B, NAITO C J. An experimental investigation of the effect of nose shape on fragments penetrating GFRP [J]. International Journal of Impact Engineering, 2014, 63: 63–71. doi: 10.1016/j.ijimpeng.2013.08.002 [6] 邓云飞, 蔡雄峰, 曾宪智, 等. 弹体头部形状对碳纤维层合板抗冲击性能影响实验研究 [J]. 振动与冲击, 2020, 39(7): 86–92.DENG Y F, CAI X F, ZENG X Z, et al. Tests for effects of projectile nose shape on anti-impact performance of carbon fiber reinforced plates [J]. Journal of Vibration and Shock, 2020, 39(7): 86–92. [7] 谢文波, 张伟, 姜雄文. 钢球斜侵彻碳纤维复合材料板的实验研究 [J]. 爆炸与冲击, 2018, 38(3): 647–653.XIE W B, ZHANG W, JIANG X W. Oblique penetration on CFRPs by steel sphere [J]. Explosion and Shock Waves, 2018, 38(3): 647–653. [8] 秦建兵, 韩志军, 刘云雁, 等. 复合材料层合板侵彻行为的研究 [J]. 振动与冲击, 2013, 32(24): 122–126.QIN J B, HAN Z J, LIU Y Y, et al. Penetration behavior of composite laminated plates [J]. Journal of Vibration and Shock, 2013, 32(24): 122–126. [9] 肖露, 程建芳, 柴晓明, 等. 层间混杂复合材料的弹道侵彻性能研究 [J]. 浙江理工大学学报, 2013, 30(4): 471–476.XIAO L, CHENG J F, CHAI X M, et al. Study on ballistic penetration property of interply hybrid composite material [J]. Journal of Zhejiang Sci-Tech University, 2013, 30(4): 471–476. [10] 肖文莹, 李想, 郭万涛, 等. Kevlar/UHMWPE混杂纤维复合材料抗弹性能 [J]. 玻璃钢/复合材料, 2019(9): 79–84.XIAO W Y, LI X, GUO W T, et al. Resistance performance of Kevlar/UHMWPE fiber hybrid composites [J]. Fiber Reinforced Plastics/Composites, 2019(9): 79–84. [11] 秦溶蔓, 朱波, 乔琨, 等. 陶瓷/纤维复合材料层间混杂结构对装甲板抗侵彻性能的影响 [J]. 材料导报, 2020, 34(9): 18183–18187.QIN R M, ZHU B, QIAO K, et al. Effect of hybrid structure of ceramic/fiber composite material on penetration resistance of armor plate [J]. Materials Reports, 2020, 34(9): 18183–18187. [12] 黄献聪. 军用防弹衣的性能需求与发展方向[C]//2005现代服装纺织高科技发展研讨会. 北京: 北京纺织工程学会, 2005: 396−402.HUANG X C. Performance requirements and development direction of military body armor [C]//2005 Symposium on High Tech Development of Modern Clothing and Textile. Beijing: Beijing Textile Engineering Society, 2005: 396−402. [13] 纤维增强复合材料抗破片模拟弹性能试验方法 V50法: GB/T 32497–2016 [S]. 北京: 中国标准出版社, 2016.Test method of ballistic resistance against fragment simulating projectiles for fiber-reinforced composites—V50 method: GB/T 32497–2016 [S]. Beijing: China Standard Press, 2016. [14] 刘国繁. 层合结构复合材料抗弹机理研究及模拟仿真[D]. 南京: 南京航空航天大学, 2015: 15−16.LIU G F. Study on numerical simulation of ballistic damage mechanisms of laminated composites [D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2015: 15−16. [15] 唐昌州, 智小琦, 徐锦波, 等. 小钨球对防弹衣加松木靶的侵彻研究 [J]. 高压物理学报, 2020, 34(5): 055101.TANG C Z, ZHI X Q, XU J B, et al. Research on small tungsten spheres penetrating into pine target with body armor [J]. Chinese Journal of High Pressure Physics, 2020, 34(5): 055101. [16] CHANG F K, CHANG K Y. Post-failure analysis of bolted composite joints in tension or shear-out mode failure [J]. Journal of Composite Materials, 1987, 21(9): 809–833. doi: 10.1177/002199838702100903 [17] CHANG F K, CHANG K Y. A progressive damage model for laminated composites containing stress concentrations [J]. Journal of Composite Materials, 1987, 21(9): 834–855. doi: 10.1177/002199838702100904 [18] 胡年明, 陈长海, 侯海量, 等. 高速弹丸冲击下复合材料层合板损伤特性仿真研究 [J]. 兵器材料科学与工程, 2017, 40(3): 66–70.HU N M, CHEN C H, HOU H L, et al. Simulation on damage characteristic of composite laminates under high-velocity projectile impact [J]. Ordnance Material Science and Engineering, 2017, 40(3): 66–70. [19] 乔咏梅, 余铜辉. 软质防弹衣结构分析与性能研究 [J]. 警察技术, 2017(3): 81–84.QIAO Y M, YU T H. Structure analysis and performance research of soft body armor [J]. Police Technology, 2017(3): 81–84. -