Numerical Simulation of Anti-Penetration of Al/CFRP/Hybrid Honeycomb Aluminum Composite Sandwich Multilayer Structure
-
摘要: 结合金属/复合材料层合结构的抗侵彻能力,基于混合蜂窝结构低成本、高韧性以及在低速冲击下吸能的特点,设计了一种Al/CFRP(carbon fiber reinforced plastics)/混合蜂窝铝复合夹芯多层结构,旨在利用各层结构特点,逐步降低弹体速度,高效吸收弹体动能,以达到防护效果。为探究Al/CFRP/混合蜂窝铝复合夹芯多层结构在弹体侵彻下的损伤演化规律及吸能特性,开展了Al/CFRP/混合蜂窝铝复合夹芯多层结构在弹体侵彻下的数值分析,探讨了冲击能量对多层结构抗侵彻性能的影响。结果表明:与Al/CFRP复合结构相比,引入混合蜂窝铝后,结构给予弹体的反作用力增大,在能量不变的情况下,弹体作用板的时间变短。在Al/CFRP/混合蜂窝铝复合夹芯多层结构抗侵彻过程中,Al板和CFRP芯层主要抵抗侵彻以降低弹体速度,混合蜂窝铝主要是吸能。在 40 J的冲击能量下,结构总吸能为36.79 J,比吸能为0.217 J/g,蜂窝铝芯层吸能占主要部分,吸能比率为30.3%;随着冲击能量的增大,蜂窝铝芯层的吸能比率增至56.2%,即冲击能量较大时蜂窝铝芯层的吸能效果更好。Abstract: Due to the low cost, high toughness and energy absorption characteristics of hybrid honeycomb structure under low velocity impact, an Al/carbon fiber reinforced plastics (CFRP)/hybrid honeycomb aluminum composite sandwich multilayer structure was designed. The kinetic energy of the projectile was supposed to be effectively absorbed and the protection was supposed to be achieved through gradually reducing the velocity of the projectile layer by layer. In order to investigate the damage evolution law and energy absorption characteristics, numerical analysis was carried out, and the impact energy effect on the penetration resistance of multilayer structure was discussed. It is found that, compared with the Al/CFRP composite structure, the reaction force given by the structure becomes larger for hybrid honeycomb aluminum. Hence, with an identical energy, the time of the projectile acting on the plate becomes shorter. In the process of anti-penetration of Al/CFRP/hybrid aluminum honeycomb composite sandwich multilayer, the Al plate and CFRP core mainly resist the penetration, and the honeycomb aluminum mainly absorbs the energy of the projectile. When the impact energy is 40 J, the total absorbed energy is 36.79 J, and the specific energy absorption is 0.217 J/g, the honeycomb aluminum core layer absorbs the main part of the energy with the proportion of 30.3%; as the impact energy increases, the proportion increases to 56.2%. This indicates that the energy absorption of the honeycomb aluminum core layer is better when the impact energy increases.
-
薄壁结构的抗冲击特性在航天航空、机械工程和交通运输等不同领域的设计中极具重要性[1-3]。随着工程领域对材料轻量化的要求,多孔结构因具有质量轻、高比能、高容积效率等优势[4-5],常作为缓冲、吸能元件,被广泛用作飞机、导弹及汽车等行业的结构件。
王璠等[6]通过对纤维增强复合材料圆柱壳的吸能特性进行试验,研究了铺层方式和引发方式对该结构吸能效果的影响和压溃过程的破坏机理。罗国强等[7]采用离散元方法模拟了多种孔洞排布方式的PMMA多孔材料在冲击加载过程中早期孔洞塌缩破坏、应力分布与粒子速度等冲击响应行为,得到了孔洞破坏形式、应力集中区域等结果。Bich等[8]提出了一种有缺陷偏心加筋功能梯度圆柱薄壳在轴压作用下的非线性静力和动力屈曲分析方法。Zhang等[9]从理论、实验与数值模拟3个方面论证了多胞金属管的吸能效果高于单胞金属管。
随着多孔材料在工程领域应用的逐渐深入,将多孔材料壳体作为吸能结构成为国内外学者关注的焦点[10-12],诸多学者对轻质多孔结构的力学特性开展了研究[13-15]。因此,研究多孔圆柱壳结构在冲击载荷下的吸能以及破坏特性具有重要的工程背景和一定的理论价值。
本研究以正、负泊松比两种多孔构型的层合圆柱壳模型为基础,以层合方式和孔单元数目为变量设计模型。在轴向压缩实验的基础上,使用有限元模拟软件ABAQUS对设计的模型进行冲击模拟分析。通过分析数值模拟结果,得到参数变化对层合多孔圆柱壳吸能特性的影响。
1. 模型介绍
研究层合多孔圆柱壳的层合方式以及孔单元数目对其冲击吸能特性的影响,在保证质量不变的前提下,设计不同层合方式和孔数的六边与内凹多孔圆柱壳模型。图1为3种不同层合方式的圆柱壳模型,分别为外多孔层、内实心层的2层层合圆柱壳,内外为多孔层、中间为实心层的3层圆柱壳,以及多孔层与实心层相间排列的4层圆柱壳(最外为多孔层)。图2为改变孔单元数目的圆柱壳模型,轴向孔数分别为8、12、16,其对应的周向孔数为24、36、48,表1给出了各个模型的几何尺寸。
表 1 层合多孔圆柱壳的几何尺寸Table 1. Geometric dimensions of laminated cellular cylindrical shellNumber of axial cell Number of circumferential cell Thickness of each layer/mm Height/mm Diameter/mm 8 24 2/2 106 78.4 12 36 2/2 106 78.4 16 48 2/2 106 78.4 8 24 1/2/1 106 78.4 8 24 1/1/1/1 106 78.4 2. 实验与数值模拟验证
如图3所示,采用聚乳酸(PLA)材料制备实验所需的2层六边多孔圆柱壳与内凹多孔圆柱壳,使用万能材料试验机进行压缩实验,所用PLA材料的密度为1.25 g/cm3。实验采用位移加载的方式压缩85 mm,压缩速度为1 mm/min。使用有限元模拟软件ABAQUS完成准静态压缩数值模拟分析,同样施加85 mm位移荷载。有限元模型如图4所示,使用刚体板模拟上下压头,中间圆柱壳模型使用四面体单元,圆柱壳网格尺寸设置为2 mm。对下刚板完全固定,上刚板轴向压缩,考虑刚性板与圆柱壳的接触摩擦,摩擦因数取0.3。数值模拟分析使用的材料参数由压缩实验测定,如图5所示,圆柱体试样尺寸为
∅ 15 mm×45 mm。图5同时给出了PLA的应力-应变曲线,得到PLA的弹性模量E=1.57 GPa,泊松比μ =0.35。图6给出了实验与数值模拟的力-位移曲线对比。可以看到,实验与数值模拟结果在模型被完全压实之前吻合较好,峰值荷载的大小和出现的位置一致。其中六边圆柱壳的压缩力波动较大,而负泊松比内凹圆柱壳的压缩力波动平缓。在数值模拟中未考虑构件断裂,其误差影响主要出现在大应变阶段,这也是数值模拟的密实化压缩位移均小于实验值的原因,随着构件断裂,结构进一步被压缩密实。图7为实验与数值模型在相同应变时的变形模式对比。可以看到,数值模拟结果能够较好地模拟实验的变形模式。
综合图6和图7可以看出,结构的变形可以划分为5个阶段:第1阶段,结构受压缩荷载一直达到初始峰值荷载,此阶段可以认为是弹性变形阶段,结构未发生明显形变;第2阶段,结构发生屈服,承受的压缩力快速下降,塑性变形明显,出现第1个褶皱;结构出现褶皱之后,结构承载力在波动中开始上升,达到第2个较低的峰值荷载,为第3阶段;第4阶段,结构承载力继续下降,出现第2个褶皱,与第1个褶皱十字相叠;第5阶段,结构基本压实,承载力持续上升。
3. 冲击仿真模拟
3.1 有限元网格及算法验证
使用经过实验测试的PLA材料参数,对六边与内凹圆柱壳有限元模型进行网格敏感性验证,分别使用尺寸为2、3、4 mm的四面体单元划分网格,图8为不同尺寸网格的六边圆柱壳有限元模型。图9为两种圆柱壳在50 m/s冲击速度下的力-位移曲线。可以看到,采用2和3 mm网格计算的力-位移曲线基本接近,4 mm网格则有一些差别,考虑计算精度与效率之后选择了精度更高的2 mm网格进行计算。沙漏现象是存在于有限元分析中的一种零能变形模式,产生零应变和应力。伪应变能为控制沙漏变形所耗散的主要能量,图10为伪应变能(Ea)与内能(Ei)之比,其比值越小,说明伪应变能较小,即沙漏现象不明显。由图10可知,伪应变能与内能之比小于1%,说明数值模拟中没有出现明显的沙漏现象,可知有限元模拟分析中网格划分较合理。
3.2 冲击速度对变形模式的影响
数值模拟中采用2 mm网格,图11为六边多孔圆柱壳在不同冲击速度下的变形模式。从图11中可以看到,随着应变的增加,模型逐渐出现褶皱。
ε =0.1时,不同速度的计算结果均出现中心对称的圆形褶皱,继续加载时,冲击速度为10 m/s的模型出现了不规则的变形模式,20、50、80 m/s冲击速度的模型则继续出现圆形褶皱。同样的现象可以在图12的内凹圆柱壳变形模式中观察到,冲击速度为10 m/s的内凹圆柱壳模型在第1个圆形褶皱之后也出现了不规则变形,速度提升之后则只出现圆形褶皱。研究[16]表明,圆柱壳在轴向冲击载荷下的变形模式受冲击速度的影响,存在变形模式变化的临界速度。本研究中,当冲击速度超过某一临界速度时,层合多孔圆柱壳的变形模式由圆形褶皱混合不规则变形转变为仅有圆形褶皱。3.3 层合方式的影响
由图11和图12中的变形模式可以看到,中高速冲击得到的变形模式较规则和稳定,选择冲击速度v=50 m/s分析层合方式与孔单元数目对层合多孔圆柱壳冲击吸能特性的影响。使用总吸能(Et)、峰值压缩力(Fp)、平均压缩力(Fa) 3种评价指标分析层合多孔圆柱壳的冲击吸能特性。
总吸能Et为结构受荷载变形吸收的能量,可由力-位移曲线中力对位移的积分得到
Et=∫l0F(x)dx (1) 式中:F为冲击压缩力,
l 为冲击压缩的距离。峰值压缩力Fp为结构在冲击压缩密实之前的压缩力峰值,通常出现在冲击荷载的初始阶段。平均压缩力Fa可以由总吸能得到
Fa=Et/l (2) 图13为不同层合方式多孔圆柱壳模型在冲击荷载作用下的力-位移曲线。从图13可以看出,不同层合方式圆柱壳模型的力-位移曲线总体比较相似,压缩力有两个较明显的峰值。对于六边多孔圆柱壳,4层圆柱壳初始峰值压缩力大于2层,3层压缩力最小,内凹多孔圆柱壳初始峰值压缩力均较接近。从六边与内凹圆柱壳的力-位移曲线可以看出,3层圆柱壳的第2个峰值压缩力远大于2、4层圆柱壳,甚至超过了初始峰值压缩力。在第2个峰值之后,4层圆柱壳的压缩力高于2、3层圆柱壳。
计算位移达到80 mm时的峰值压缩力和平均压缩力,结果如表2所示。4层圆柱壳的总吸能与平均压缩力皆大于其他两种结构,3层圆柱壳的总吸能与平均压缩力最小,峰值压缩力却最大,分析是由内外多孔层导致的。对比正、负泊松比多孔圆柱壳的吸能特性,2、3、4层内凹圆柱壳的总吸能较六边圆柱壳分别提高17.4%、22.5%、18.3%,峰值压缩力相差不大,平均压缩力更高。
表 2 不同层合方式下的冲击吸能指标Table 2. Impact energy absorption parameters of different laminated wayCellular units Layer number Et/J Fp/kN Fa/kN Hexagon 2 1 206.0 29.9 12.2 3 1 032.1 43.9 11.0 4 1 414.3 32.6 14.5 Re-entrant 2 1 416.7 30.8 13.7 3 1 264.6 48.6 11.8 4 1 673.5 32.6 17.0 3.4 孔单元数目的影响
图14为不同孔单元数目圆柱壳模型在冲击载荷作用下的力-位移曲线。可以看到,相同质量下多孔圆柱壳的孔单元数目对力-位移曲线的影响较小,3种孔数的曲线都较接近。六边圆柱壳有2个明显的峰值压缩力,内凹圆柱壳较平缓,在初始峰值压缩力之后有3个较小的峰值压缩力。
同样计算总吸能Et、峰值压缩力Fp和平均压缩力Fa 3种吸能指标,如表3所示。可以看到,孔数最多的六边圆柱壳吸能效果更好,内凹圆柱壳则是孔数最少的模型吸能效果更好。孔数对平均压缩力的影响较小。16×48孔圆柱壳的峰值压缩力较高,12×36孔圆柱壳的峰值压缩力相对较小。内凹圆柱壳的总吸能同样大于六边圆柱壳,其在力-位移曲线上也表现出更多的峰值。
表 3 不同孔数多孔圆柱壳的冲击吸能指标Table 3. Impact energy absorption parameters of of cellular cylindrical shell with different cells numberCellular units Cell number Et/J Fp/kN Fa/kN Hexagon 8×24 1 206.0 29.9 12.2 12×36 1 169.7 25.8 11.2 16×48 1 266.4 29.4 11.9 Re-entrant 8×24 1 416.7 30.8 13.7 12×36 1 383.4 28.9 13.4 16×48 1 394.0 39.6 13.3 4. 结 论
研究了不同层合方式、不同孔单元数目的层合多孔圆柱壳在冲击荷载下的吸能特性。在轴向压缩实验的基础上,利用 ABAQUS 有限元软件对多孔圆柱壳模型在准静态、不同冲击速度下的压缩行为进行了数值模拟。基于结构的变形模式,分析了不同层合方式与孔数对其承载能力、吸能特性的影响,得出以下结论。
(1) 随着冲击速度的提升,多孔圆柱壳从轴对称的三角变形模式逐渐变成中心对称的圆形变形模式,其中速度为10 m/s的模型同时具有三角变形模式和圆形变形模式。
(2) 4层圆柱壳有更好的吸能效果与更高的平均压缩力,3层圆柱壳的峰值压缩力远大于其他两种层合圆柱壳。
(3) 孔单元数目对吸能效果的影响与多孔构型有关,孔单元数目最多的六边圆柱壳的吸能效果最好,而在内凹圆柱壳中则相反。
(4) 同样参数下,具有负泊松比特性的内凹圆柱壳的吸能效果较六边圆柱壳更好,平均压缩力也更大。
-
表 1 7075-T651铝合金的材料参数以及Johnson-Cook本构模型和失效模型参数
Table 1. Material parameters, Johnson-Cook constitutive model and failure model parameters of 7075-T651 aluminum alloy
ρ/(kg·m−3) E/GPa v A/MPa B/MPa n m Tm/K 2810 71.7 0.33 520 477 0.52 1.61 893 Tt/K D1 D2 D3 D4 D5 ˙ε0/s−1 293 0.096 0.049 3.465 0.016 1.099 0.0005 表 2 40 J冲击能量下峰值应力与吸能对比
Table 2. Comparison of peak stress and energy absorption when the impact energy is 40 J
Categories Peak force/kN Energy absorption/J This paper 3.48 21.42 Ref. [14] 3.58 22.08 Error/% 2.8 2.9 ρ/(kg·m−3) E11/GPa E12/GPa E13/GPa G12/GPa G13/GPa G23/GPa ν12 ν13 1560 141 9.7 9.7 5.2 5.2 3.4 0.34 0.34 ν23 Xt/MPa Xc/MPa Yt/MPa Yc/MPa S12/MPa S13/MPa S23/MPa 0.44 2703 1737 81 312 57 57 57 表 4 Al/CFRP/混合蜂窝铝复合夹芯多层结构的吸能情况
Table 4. Energy absorption of Al/CFRP/hybrid honeycomb aluminum composite sandwich multilayer structure
Core layer Energy absorption/J Mass/g Energy absorption ratio/% Specific energy absorption/(J·g−1) Upper aluminum plate 9.38 28.1 27.9 0.33 CFRP 9.59 23.4 28.5 0.41 Honeycomb aluminum 10.22 89.3 30.4 0.11 Lower aluminum plate 4.47 28.1 13.2 0.16 表 5 冲击能量的影响
Table 5. Effect of impact energy
Impact energy/J Peak force/kN Maximum displacement/mm Energy absorption/J Specific energy absorption/(J·g−1) 40 12.7 7.0 36.8 0.218 80 16.4 9.0 71.8 0.425 120 17.1 11.6 107.0 0.633 200 17.2 16.1 162.3 0.959 -
[1] 孙卫兵. 纤维增强复合材料层合板抗高速破片侵彻性能研究 [D]. 武汉: 武汉理工大学, 2020.SUN W B. Research on penetration resistance of fiber reinforced composite laminates under high-speed fragments [D]. Wuhan: Wuhan University of Technology, 2020. [2] LI L J, SUN L Y, WANG T K, et al. Repeated low-velocity impact response and damage mechanism of glass fiber aluminium laminates [J]. Aerospace Science and Technology, 2019, 84: 995–1010. doi: 10.1016/j.ast.2018.11.038 [3] 马小敏, 李世强, 李鑫, 等. 编织Kevlar/Epoxy复合材料层合板在冲击荷载下的动态响应 [J]. 爆炸与冲击, 2016, 36(2): 170–176. doi: 10.11883/1001-1455(2016)02-0170-07MA X M, LI S Q, LI X, et al. Dynamic response of woven Kevlar/Epoxy composite laminates under impact loading [J]. Explosion and Shock Waves, 2016, 36(2): 170–176. doi: 10.11883/1001-1455(2016)02-0170-07 [4] 金子明, 沈峰, 曲志敏, 等. 纤维增强复合材料抗弹性能研究 [J]. 纤维复合材料, 1999, 16(3): 5–9.JIN Z M, SHEN F, QU Z M, et al. A study of anti-ballistic properties of FRP [J]. Fiber Composites, 1999, 16(3): 5–9. [5] ZHOU J J, WEN P H, WANG S N. Numerical investigation on the repeated low-velocity impact behavior of composite laminates [J]. Composites Part B: Engineering, 2020, 185: 107771. doi: 10.1016/j.compositesb.2020.107771 [6] JAROSLAW B, BARBARA S, PATRYK J. The comparison of low-velocity impact resistance of aluminum/carbon and glass fiber metal laminates [J]. Polymer Composites, 2016, 37(4): 1056–1063. doi: 10.1002/pc.23266 [7] RYAN S, SCHAEFER F, DESTEFANIS R, et al. A ballistic limit equation for hypervelocity impacts on composite honeycomb sandwich panel satellite structures [J]. Advances in Space Research, 2008, 41(7): 1152–1166. doi: 10.1016/j.asr.2007.02.032 [8] CHRISTIANSEN E L. Design and performance equations for advanced meteoroid and debris shields [J]. International Journal of Impact Engineering, 1993, 14(1): 145–156. doi: 10.1016/0734-743X(93)90016-Z [9] ZHANG D H, FEI Q G, ZHANG P W. Drop-weight impact behavior of honeycomb sandwich panels under a spherical impactor [J]. Composite Structures, 2017, 168: 633–645. doi: 10.1016/j.compstruct.2017.02.053 [10] MORADA G, OUADDAY R, VADEAN A, et al. Low-velocity impact resistance of ATH/epoxy core sandwich composite panels: experimental and numerical analyses [J]. Composites Part B: Engineering, 2017, 114: 418–431. doi: 10.1016/j.compositesb.2017.01.070 [11] GUO K L, ZHU L, LI Y G, et al. Experimental investigation on the dynamic behaviour of aluminum foam sandwich plate under repeated impacts [J]. Composite Structures, 2018, 200: 298–305. doi: 10.1016/j.compstruct.2018.05.148 [12] XU M M, HUANG G Y, DONG Y X, et al. An experimental investigation into the high velocity penetration resistance of CFRP and CFRP/aluminium laminates [J]. Composite Structures, 2018, 188: 450–460. doi: 10.1016/j.compstruct.2018.01.020 [13] 朱倩. 纤维金属层板抗高速冲击性能及损伤机理研究 [D]. 镇江: 江苏大学, 2020.ZHU Q. Study on impact resistance and damage mechanism of fiber metal laminates under high velocity impact [D]. Zhenjiang: Jiangsu University, 2020. [14] 李毅翔. Al/CFRP混杂层合板抗低速冲击性能研究 [D]. 长沙: 湖南大学, 2020.LI Y X. Study on low velocity impact resistance of Al/CFRP hybrid structures [D]. Changsha: Hunan University, 2020. [15] HASHIN Z. Failure criteria for unidirectional fiber composites [J]. Journal of Applied Mechanics, 1980, 47(2): 329–334. doi: 10.1115/1.3153664 [16] 刘礼平, 段科好, 徐卓, 等. 碳纤维增强树脂基复合材料层合板胶螺混合连接失效机制 [J]. 复合材料学报, 2023, 40(1): 592–602. doi: 10.13801/j.cnki.fhclxb.20220215.001LIU L P, DUAN K H, XU Z, et al. Failure mechanism of carbon fiber reinforced polymer bonded-bolted hybrid connection [J]. Acta Materiae Compositae Sinica, 2023, 40(1): 592–602. doi: 10.13801/j.cnki.fhclxb.20220215.001 [17] XU M C, LIU D B, WANG P D, et al. In-plane compression behavior of hybrid honeycomb metastructures: theoretical and experimental studies [J]. Aerospace Science and Technology, 2020, 106: 106081. doi: 10.1016/j.ast.2020.106081 -