摘 要 ：如今直升机已广泛应用于军事和民用领域，但目前对其抗高速冲击性能的研究较少。通过对球型破片高速冲击蜂窝夹芯结构的数值仿真，研究了蜂窝夹芯结构在高速冲击下的动态响应过程。结果表明，在破片的高速冲击作用下，蜂窝夹芯结构的破坏形式不再像低速冲击一样，而是以局部破坏为主。同时，芯层结构也由于参与程度较少，其对抗冲击作用的影响远小于面板，但芯层的存在仍优于实心结构。基于夹芯结构的重量和防护效果，分析几何结构同抗高速冲击的防护能力之间的关系。

关键词 ：蜂窝夹芯结构 ；抗冲击性能 ；数值仿真 ；耐撞性

中图分类号 ：V214.6 文献标识码 ：B 文章编号 ：1671-2064(2023)11-0071-06

0引言

直升机由于其具有平移飞行、悬停飞行和垂直起降的能力，使其成为所有飞行器中用途最广泛的类型之一 ^[1]。迄今为止，直升机已经广泛应用于军事和民用领域，然而，近几十年来，直升机安全问题越来越受到关注。对直升机事故调查发现，飞机耐撞性被视为关键问题之一。但直到越南战争时期才引起足够的重视，第一批优先考虑防撞设计的是 UH-60 黑鹰和 AH-64 阿帕奇直升机 ^[2]。在之后的开发中，耐撞性设计变得越来越重要。而在各种典型的能量吸收结构中，蜂窝结构由于具有高比强度、比刚度和显著的能量吸收性能，在大多数工程领域得到了广泛应用 ^[3]。迄今为止，在耐撞性研究中，人们开发了许多蜂窝夹层结构来提高能量吸收能力。

目前Sun 等人 ^[4] 研究了基于一阶和二阶顶点的分层蜂窝的平面外耐撞性行为。结果表明，一阶和二阶蜂窝的比能量吸收分别提高了 81.3% 和 185.7%，而峰值力没有增加太多。Ma 等人 ^[5] 研究了仿生自相似规则分层蜂窝在面外冲击载荷下的耐撞性。分层单元组织可以加强材料强度，从而提高抗压强度和能量吸收能力。湖南大学^[6] 和东南大学团队 ^[7] 关于蜂窝夹层结构几何属性对其抗冲击的影响有了系统性的研究。

但是，直升机在飞行过程中往往遭受高速冲击，而在这方面的研究目前还相对较少。所以本文通过 ABAQUS 建立了蜂窝夹芯结构的有限元模型，并模拟了破片侵彻夹芯结构的过程，分析了在破片高速贯穿夹芯板时，结构各部分的损伤情况。并进一步对比了各部分参数的变化对其吸能的影响，旨在进一步了解夹芯结构的面板、芯层的参数和夹芯板耐撞性之间的关系。

1数值仿真模型及参数设置

如图 1 所示，前后面板为边长 150mm 的正方形，厚度1mm ；芯层采用边长为 4mm，厚度为 0.07mm 的正六边形蜂窝，高度为 15mm。破片采用球型破片，其直径为 12mm，质量为 7.05g。由于主要研究夹芯板的动态响应过程，所以不考虑破片自身的变形情况，将破片设置为钢体结构。

在ABAQUS 中，网格的数量决定了计算结果的准确度以及计算的速度。综合考虑下，仅在夹芯板与破片接触区域进行网格的局部细化，细化后的尺寸为 0.5mm，其余地方采用渐变的方式，如图 2 所示。模型的边界条件为刚性固定，面板和芯层之间的接触关系设定为Tie 接触，破片和夹芯结构之间的碰撞关系设定为通用接触，其中法线方向为硬接触，切线方向为摩擦接触，系数为 0.3。

在本次数值仿真中，面板的材料为 5052 铝，芯层的材料为 3003 铝，在 ABAQUS 中的塑性模型和失效模型选用 Johnson-Cook 本构方程 ^[8-9]。公式如下。

式中，材料参数的值取自《国际冲击工程》中公布的数据^[9]，具体的参数如表 1 所示。

2数值计算方法的验证

为了验证数值仿真方法的合理性，本文将仿真结果同Sun 等人 ^[10] 的试验结果进行对比。如图 3 所示，夹芯板的尺寸与模型中选用的保持一致。在试验过程中，夹芯板分别在 76.2m/s、138m/s、172.4m/s 3 种速度下出现了前面板凹陷、破片反弹、前面板被贯穿，破片陷入芯层中和前后面板均被贯穿 3 种毁伤形式。所以在数值仿真中也将破片模拟了这 3 种速度，能更为准确地验证ABAQUS 仿真结果的合理性。

破片在数值仿真和试验过程中的剩余速度结果如表 2 所示，对比分析可得，ABAQUS 基本上能够模拟出破片冲击夹芯板后的剩余速度。由于在制造夹芯板时，面板和芯层之间采用胶水粘连的工艺，从而会多一层黏合层，并且试验中螺栓固定的程度并非完全刚性固定。而在有限元模型中为理想状态。因此在入射速度较小时，仿真结果与试验结果有一定的误差，但在入射速度增大时，其影响均会减小。

夹芯板的面板毁伤情况如图 4 所示，通过对比试验和仿真的面板变形结果来看，ABAQUS 能够较好地模拟夹芯板面板的破坏情况。在侵彻速度为 76.2m/s 时，面板仅出现凹陷变形，其变形的挠度大小基本上与试验结果一致 ；在侵彻速度为 138m/s，部分穿透 ；在侵彻速度为172.4m/s 时，完全穿透 ；从贯穿破孔孔径来看，其数据也十分接近。所以，可以使用 ABAQUS 来对夹芯板的毁伤情况进行有限元分析。

本文主要研究夹芯结构的抗高速冲击性能，所以选用完全贯穿的夹芯板毁伤情况进行对比，如图 5 所示。在破片完全侵彻夹芯板时，前面板收到侵彻而破碎，冲击区域的芯层被完全破坏，远离冲击区域的芯层受到挤压屈曲变形，并与后面板出现脱粘，而这些破坏情况在数值仿真中也得到了体现。综合 3 种情况，说明可以使用 ABAQUS 来仿真夹芯结构在高速破片侵彻下的动态响应过程。

3计算结果与分析

3.1物理过程分析

数值仿真比起试验能够更为直观地观察夹芯板内部的毁伤过程。如图 6 所示，当破片首先撞击到前面板时，使面板冲塞变形。同时面板挤压芯层使芯层出现屈曲破坏， 并与前面板发生挤压破坏。随后后面板遭到贯穿，形成充塞孔，附带周围相连的芯层受到撕裂而出现“脱粘”。但同时，远离冲击位置的部分仍保持完整的结构，破片对夹芯板的毁伤影响能力仅存在于冲击区域内。

3.2能量吸收特性

从图 7 中可以看出，整个破坏过程均遵循能量守恒定律，破片的动能转换为了夹芯板的塑性应变能。在 3 种速度的破片侵彻下，由于夹芯板的毁伤情况不同，夹芯板各部分吸收的能量不同，主要是当速度较低时，后面板并未参与整体抗冲击过程，所以其几乎不吸收能量，而前面板由于发生变形较大，所以其吸收的能量也始终最多，而随速度的增大，后面板参与抵抗冲击后，其吸收的能量就超过了芯层。

结合图 7 与图 4 可知，夹芯结构的吸能情况与其自身的毁伤情况有关，说明夹芯板是通过塑性应变吸收破片的动能，而芯层的存在连接了前后两个面板，使其整体参与抵抗破片的冲击，因此，可以推断夹芯板的抗冲击性能优于两块面板单独抵抗。但由于在高速冲击下，夹芯板的破坏是局部性破坏，所以面板仍然是抗冲击的主要部分。

数值仿真可以改变破片侵彻的速度，通过数值仿真设置了一系列不同侵彻速度下的夹芯板吸能情况对比，如图8 所示。从中可以看出，当破片贯穿结构后，结构的吸收能量就达到了极值，此时再增加入射速度，并不会对结构的吸能情况造成影响，相反入射速度越快，其破坏范围越集中，芯层所贡献的抗冲击性能效果反而有所降低。而前后面板均产生冲塞孔的毁伤情况，所以其吸收能量基本保持一致。

4蜂窝夹芯板抗冲击性能参数研究

4.1几何参数对弹道极限速度的影响

弹道极限速度表示破片完全贯穿目标所需要的最小冲击速度 ^[11]，因此，在本次研究中将其作为评价夹芯板抗高速冲击能力的标准之一。通过在 ABAQUS 中修改各部位的参数，得到了夹芯板在面板、芯层在不同几何参数下的弹道极限速度，如图 9 所示。从中可以看出，作为主要抗冲击的面板，随着厚度的增加，弹道极限速度也跟着增加。而改变芯层的几何参数，则对抗冲击性能影响较小。增加胞元壁厚和减小胞元边长等效于提高了破片粉碎芯层所需要的动能，使得更多的芯层参与抗冲击作用，所以在一定程度上提高了弹道极限速度，但是影响很小，而增加芯层高度则几乎没有影响。

4.2几何参数对能量吸收的影响

为了研究高速冲击下夹芯板的动态响应，在 ABAQUS 中将破片的速度设置为 250m/s，实现所有工况下夹芯板均完全贯穿。得到的结果如图 10 所示，当面板厚度增加时，面板自身对破片动能的吸收也会增加，同时破片侵彻面板后的速度会降低，从而使结果更趋近于低速冲击，带动更多的芯层参与抵抗破片的侵彻，从而提高了夹芯板整体结构以及各个部分的能量吸收效率。而仅改变芯层参数时，因为芯层自身结构强度很弱，在完全贯穿的情况下，芯层相对于面板影响很小，所以即使改变芯层的几何参数，也因为自身比重太低而对整个结果影响很小。

4.3几何参数对比吸能的影响

在直升机的结构设计过程中，重量是不可忽略的元素，这就要求应用于直升机上的结构尽可能满足轻量化设计。所以在研究夹芯板的抗高速冲击性能时，也需要考虑质量对吸能效果的影响，比吸能则是一个重要的指标。如图 11 所示，增加了质量因素之后，得到了不同几何参数下蜂窝夹芯板的比吸能结果。从前文研究可得，增加面板厚度能够提高夹芯板的吸能效率，但同时也会增加夹芯板的重量，而从图 11（a）可以看出，当面板厚度增加超过1.6mm 之后，再增加厚度，反而比吸能开始降低，这说明此时吸能效率的增加开始小于质量的增加。而改变芯层的几何参数，对结构的质量和吸能都影响较小，所以整体的结果变化不大，但从图 11（b）、11（c）、11（d）中也能看出，在波动的状态下也能有一个比吸能较高的结果。

综合研究弹道极限速度、吸能情况以及重量的因素可以得出，蜂窝夹芯板的抗高速冲击能力主要与面板有关， 最直接提高防撞能力的方法就是增加面板的厚度，但同时增加面板厚度会使夹芯结构质量也随之增加，增加到一定程度之后反而会降低吸能效率。而芯层对夹芯板抗高速冲击的性能影响较弱，所以改变芯层的几何参数对整体抗冲击作用也影响较小。

5结语

本文通过 ABAQUS 建立了蜂窝夹芯结构的有限元模型，对其抗高速冲击性能进行了数值仿真研究，得出以下几点结论。

（1）与低速冲击不同，高速冲击下的蜂窝夹芯板主要的毁伤情况是前后面板形成冲塞孔，芯层结构出现变形，并且冲击区域的芯层会与面板完全脱粘。（2）破片在高速冲击的过程中，蜂窝夹芯板通过自身的塑性应变来吸收破片的动能，所以在贯穿整体后，前后面板的破坏形式一致，变形一致，吸能量也比较接近。（3）在抗高速冲击时，起主要作用的是夹芯板的面板，所以直接提高夹芯结构的抗高速冲击性能的方式就是增加面板厚度，但结合轻量化设计时，则需要考虑质量对其的影响。（4）芯层虽然对夹芯板在低速冲击下的防护性能有影响，但夹芯结构在高速冲击过程中，由于损伤更为局部化，而芯层自身在夹层结构中占比较低，所以改变其几何参数对整体效果影响较少。

参考文献

[1]Yang X,Ma J,Wen D,et al.Crashworthy Design and Energy Absorption Mechanisms for Helicopter Structures: A Systematic Literature Review[J].Progress in Aerospace Sciences,2020,(114):1-31.

[2]YANG PEI,BIFENG SONG.Several Key Techniques for Civil Helicopter Crashworthiness Safety Design and Simulation[C].//Asia Simulation Conference 2008/the 7th International Conference on System Simulation and Scientific Computing(2008亚洲仿真会议/第七届系统仿真和科学计算国际会议)(ICSC'2008)论文集,2008:1308-1311.

[3]Zhang Y,Lu M,Wang C H,et al.Out-of-plane Crashworthiness of Bioinspired Self-similar Regular Hierarchical Honeycombs, Compos.Struct,2016,(144):1-13.

[4]Taylor C M,Smith C W,Miller W,et al.Functional Grading in Hierarchical Honeycombs:Density Specific Elastic Performance,Compos[J].Struct,2012,(94):2296–2305.

[5]Yahaya M,Ruan D,Lu G,et al.Response of Aluminium Honeycomb Sandwich Panels Subjected to Foam Projectile Impact-an Experimental Study[J].International Journal of Impact Engineering,2015,(75):100–109.

[6]Sun G,Chen D,Huo X,et al.Experimental and Numerical Studies on Indentation and Perforation Characteristics of Honeycomb Sandwich Panels[J].Composite Structures, 2018,184(JAN.1):110-124.

[7]张大海.冲击载荷下蜂窝及其夹芯结构力学行为研究[D].南京: 东南大学,2018.

[8]Johnson G R.A Constitutive Model and Data for Metals Subjected to Large Strains,High Strain Rates and High Temperature[J].1983.

[9]Børvik T,Forrestal M J,Hopperstad O S,et al.Perforation of AA5083-H116 aluminium plates Plates with conicalConical- nose steel Steel projectiles Projectiles-calculations Calculations [J].International Journal of Impact Engineering.2009,(36): 426–37.

[10]Sun G Y,Chen D D,Wang H X,et al.High-velocity Impact Behaviour of Aluminium Honeycomb Sandwich Panels with Different Structural Configurations[J].International Journal of Impact Engineering,2018,(122):119-136.

[11]Hongxu Wang,Karthik Ram Ramakrishnan,Krishna Shankar.Experimental Study of the Medium Velocity Impact Response of Sandwich Panels with Different Cores[J].Materials & Design,2016,99(5):68-82.

Dynamic Response of Honeycomb Sandwich Panels Under High Velocity Impact

HU Fangliang1,YU Sheng2

(1.China Helicopter Research and Development Institute, Jingdezhen Jiangxi 333001;

2.Wuhan Bureau of Haizhuang Military Representative Office in Nanchang, Nanchang Jiangxi 333000)

Abstract:Helicopters have been widely used in military and civil fields, but there is little research on their high-speed impact performance.The dynamic response process of honeycomb sandwich structure under high velocity impact is studied by numerical simulation of honeycomb sandwich structure under high velocity impact of spherical fragments. The results show that under the high-speed impact of fragments, the failure mode of honeycomb sandwich structure is no longer the same as that of low-speed impact, but mainly local failure.At the same time, the influence of the core structure on the impact resistance is far less than that of the panel because of its less participation, but the existence of the core layer is still better than that of the solid structure.Finally, the relationship between the geometric structure and the protection ability against high-speed impact is obtained by synthesizing the weight and protection effect of the sandwich structure.

Key words: honeycomb sandwich panels;impact resistance;numerical simulation;crashworthiness