薄壁结构的抗冲击特性在航天航空、机械工程和交通运输等不同领域的设计中极具重要性^[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  不同层合方式的圆柱壳模型

Figure  1.  Cylindrical shell models withdifferent laminated ways

下载: 全尺寸图片 幻灯片

表  1  层合多孔圆柱壳的几何尺寸

Table  1.  Geometric dimensions of laminated cellular cylindrical shell

Number 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

下载: 导出CSV 

| 显示表格

图  2  不同孔单元数目的圆柱壳模型

Figure  2.  Cylindrical shell models withdifferent cells number

下载: 全尺寸图片 幻灯片

2.   实验与数值模拟验证

如图3所示，采用聚乳酸（PLA）材料制备实验所需的2层六边多孔圆柱壳与内凹多孔圆柱壳，使用万能材料试验机进行压缩实验，所用PLA材料的密度为1.25 g/cm³。实验采用位移加载的方式压缩85 mm，压缩速度为1 mm/min。使用有限元模拟软件ABAQUS完成准静态压缩数值模拟分析，同样施加85 mm位移荷载。有限元模型如图4所示，使用刚体板模拟上下压头，中间圆柱壳模型使用四面体单元，圆柱壳网格尺寸设置为2 mm。对下刚板完全固定，上刚板轴向压缩，考虑刚性板与圆柱壳的接触摩擦，摩擦因数取0.3。数值模拟分析使用的材料参数由压缩实验测定，如图5所示，圆柱体试样尺寸为$\varnothing$15 mm×45 mm。图5同时给出了PLA的应力-应变曲线，得到PLA的弹性模量E=1.57 GPa，泊松比$\,\mu$=0.35。

图  3  PLA层合多孔圆柱壳

Figure  3.  Laminated cellular cylindrical shell of PLA material

下载: 全尺寸图片 幻灯片

图  4  ABAQUS有限元模型

Figure  4.  Finite element model in ABAQUS

下载: 全尺寸图片 幻灯片

图  5  PLA材料的压缩实验与应力-应变曲线

Figure  5.  Compression experiment and stress-strain curve of PLA material

下载: 全尺寸图片 幻灯片

图6给出了实验与数值模拟的力-位移曲线对比。可以看到，实验与数值模拟结果在模型被完全压实之前吻合较好，峰值荷载的大小和出现的位置一致。其中六边圆柱壳的压缩力波动较大，而负泊松比内凹圆柱壳的压缩力波动平缓。在数值模拟中未考虑构件断裂，其误差影响主要出现在大应变阶段，这也是数值模拟的密实化压缩位移均小于实验值的原因，随着构件断裂，结构进一步被压缩密实。图7为实验与数值模型在相同应变时的变形模式对比。可以看到，数值模拟结果能够较好地模拟实验的变形模式。

图  6  实验与数值模拟得到的力-位移曲线对比

Figure  6.  Comparisons of compressive force-displacement curves between experiment and simulation

下载: 全尺寸图片 幻灯片

图  7  实验与数值模拟得到的变形模式对比

Figure  7.  Comparisons of deformation mode between experiment and simulation

下载: 全尺寸图片 幻灯片

综合图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为伪应变能（E_a）与内能（E_i）之比，其比值越小，说明伪应变能较小，即沙漏现象不明显。由图10可知，伪应变能与内能之比小于1%，说明数值模拟中没有出现明显的沙漏现象，可知有限元模拟分析中网格划分较合理。

图  8  3种尺寸网格的有限元仿真模型

Figure  8.  Finite element simulation models of three kinds of meshes

下载: 全尺寸图片 幻灯片

图  9  网格敏感性验证

Figure  9.  Validation of mesh sensitivity

下载: 全尺寸图片 幻灯片

图  10  伪应变能与内能之比

Figure  10.  Ratio of artificial strain energy to internal energy

下载: 全尺寸图片 幻灯片

3.2   冲击速度对变形模式的影响

数值模拟中采用2 mm网格，图11为六边多孔圆柱壳在不同冲击速度下的变形模式。从图11中可以看到，随着应变的增加，模型逐渐出现褶皱。$\varepsilon$=0.1时，不同速度的计算结果均出现中心对称的圆形褶皱，继续加载时，冲击速度为10 m/s的模型出现了不规则的变形模式，20、50、80 m/s冲击速度的模型则继续出现圆形褶皱。同样的现象可以在图12的内凹圆柱壳变形模式中观察到，冲击速度为10 m/s的内凹圆柱壳模型在第1个圆形褶皱之后也出现了不规则变形，速度提升之后则只出现圆形褶皱。研究^[16]表明，圆柱壳在轴向冲击载荷下的变形模式受冲击速度的影响，存在变形模式变化的临界速度。本研究中，当冲击速度超过某一临界速度时，层合多孔圆柱壳的变形模式由圆形褶皱混合不规则变形转变为仅有圆形褶皱。

图  11  六边多孔圆柱壳在不同冲击速度下的变形模式

Figure  11.  Deformation mode of hexagonal cellular cylindrical shell at different impact velocities

下载: 全尺寸图片 幻灯片

图  12  内凹多孔圆柱壳变形模式

Figure  12.  Deformation mode of re-entrant cellular cylindrical shell

下载: 全尺寸图片 幻灯片

3.3   层合方式的影响

由图11和图12中的变形模式可以看到，中高速冲击得到的变形模式较规则和稳定，选择冲击速度v=50 m/s分析层合方式与孔单元数目对层合多孔圆柱壳冲击吸能特性的影响。使用总吸能（E_t）、峰值压缩力（F_p）、平均压缩力（F_a） 3种评价指标分析层合多孔圆柱壳的冲击吸能特性。

总吸能E_t为结构受荷载变形吸收的能量，可由力-位移曲线中力对位移的积分得到

式中：F为冲击压缩力，$l$为冲击压缩的距离。

峰值压缩力F_p为结构在冲击压缩密实之前的压缩力峰值，通常出现在冲击荷载的初始阶段。平均压缩力F_a可以由总吸能得到

图13为不同层合方式多孔圆柱壳模型在冲击荷载作用下的力-位移曲线。从图13可以看出，不同层合方式圆柱壳模型的力-位移曲线总体比较相似，压缩力有两个较明显的峰值。对于六边多孔圆柱壳，4层圆柱壳初始峰值压缩力大于2层，3层压缩力最小，内凹多孔圆柱壳初始峰值压缩力均较接近。从六边与内凹圆柱壳的力-位移曲线可以看出，3层圆柱壳的第2个峰值压缩力远大于2、4层圆柱壳，甚至超过了初始峰值压缩力。在第2个峰值之后，4层圆柱壳的压缩力高于2、3层圆柱壳。

图  13  不同层合方式多孔圆柱壳的力-位移曲线

Figure  13.  Compressive force-displacement curves of cellular cylindrical shell with different laminated ways

下载: 全尺寸图片 幻灯片

计算位移达到80 mm时的峰值压缩力和平均压缩力，结果如表2所示。4层圆柱壳的总吸能与平均压缩力皆大于其他两种结构，3层圆柱壳的总吸能与平均压缩力最小，峰值压缩力却最大，分析是由内外多孔层导致的。对比正、负泊松比多孔圆柱壳的吸能特性，2、3、4层内凹圆柱壳的总吸能较六边圆柱壳分别提高17.4%、22.5%、18.3%，峰值压缩力相差不大，平均压缩力更高。

表  2  不同层合方式下的冲击吸能指标

Table  2.  Impact energy absorption parameters of different laminated way

Cellular 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

下载: 导出CSV 

| 显示表格

3.4   孔单元数目的影响

图14为不同孔单元数目圆柱壳模型在冲击载荷作用下的力-位移曲线。可以看到，相同质量下多孔圆柱壳的孔单元数目对力-位移曲线的影响较小，3种孔数的曲线都较接近。六边圆柱壳有2个明显的峰值压缩力，内凹圆柱壳较平缓，在初始峰值压缩力之后有3个较小的峰值压缩力。

图  14  不同孔数多孔圆柱壳的力-位移曲线

Figure  14.  Compressive force-displacement curves of cellular cylindrical shell with different cells number

下载: 全尺寸图片 幻灯片

同样计算总吸能E_t、峰值压缩力F_p和平均压缩力F_a 3种吸能指标，如表3所示。可以看到，孔数最多的六边圆柱壳吸能效果更好，内凹圆柱壳则是孔数最少的模型吸能效果更好。孔数对平均压缩力的影响较小。16×48孔圆柱壳的峰值压缩力较高，12×36孔圆柱壳的峰值压缩力相对较小。内凹圆柱壳的总吸能同样大于六边圆柱壳，其在力-位移曲线上也表现出更多的峰值。

表  3  不同孔数多孔圆柱壳的冲击吸能指标

Table  3.  Impact energy absorption parameters of of cellular cylindrical shell with different cells number

Cellular 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

下载: 导出CSV 

| 显示表格

4.   结　论

研究了不同层合方式、不同孔单元数目的层合多孔圆柱壳在冲击荷载下的吸能特性。在轴向压缩实验的基础上，利用 ABAQUS 有限元软件对多孔圆柱壳模型在准静态、不同冲击速度下的压缩行为进行了数值模拟。基于结构的变形模式，分析了不同层合方式与孔数对其承载能力、吸能特性的影响，得出以下结论。

(1) 随着冲击速度的提升，多孔圆柱壳从轴对称的三角变形模式逐渐变成中心对称的圆形变形模式，其中速度为10 m/s的模型同时具有三角变形模式和圆形变形模式。

(2) 4层圆柱壳有更好的吸能效果与更高的平均压缩力，3层圆柱壳的峰值压缩力远大于其他两种层合圆柱壳。

(3) 孔单元数目对吸能效果的影响与多孔构型有关，孔单元数目最多的六边圆柱壳的吸能效果最好，而在内凹圆柱壳中则相反。

(4) 同样参数下，具有负泊松比特性的内凹圆柱壳的吸能效果较六边圆柱壳更好，平均压缩力也更大。