何 强，王勇辉，史肖娜，顾 航，陈 宇

（1.江苏科技大学机械工程学院，江苏镇江212003；2.江苏理工学院机械工程学院，江苏常州213001）

薄壁结构质量轻并且具有较好的能量吸收特性，常被用来生产各种能量吸收或缓冲装置，已被广泛应用于航天航空、车辆工程及国防装备等领域。薄壁结构在轴向冲击载荷作用下通过胞壁的塑性变形吸收大量的冲击能量，所以研究薄壁结构的轴向压缩特性具有重要的工程意义。

McFarland[1]最先提出正六边形多胞结构轴向准静态平均压缩应力的计算方法，Wierzbicki[2]修正了McFarland 提出的基本折叠模式，并根据超折叠单元理论[3]对正六边形多胞结构的轴向平均压缩应力和折叠波长进行了理论求解；Tran 等[4]通过理论和数值手段获得了多种薄壁管在冲击载荷作用下的轴向压缩应力理论模型；尹汉锋等[5]基于超折叠单元理论对几种胞元构型蜂窝的平均压缩应力进行求解，并开展耐撞性优化设计。

薄壁结构的几何构型对其轴向冲击载荷作用下的吸能特性有着较大的影响。为提高普通构型（三角形、正六边形、圆形、正方形等）薄壁管的吸能能力，许多学者对一些新型层级薄壁结构产生了兴趣。Sun 等[6]将规则六边形多胞结构的每个顶点替换为一个较小的六边形结构，并分析了其轴向冲击载荷作用下的耐撞性能。Mousanezhad 等[7]提出了具有高刚度和高韧性的蜘蛛网层级多胞结构，并研究了结构参数对其面内及面外力学性能的影响。Sun 等[8]用同性子结构替换了规则六边形多胞结构的实心胞壁，研究发现该层级特性的引入能大幅度提高其面内刚度。张越等[9]运用数值模拟方法研究了二阶自相似四边形蜂窝结构参数对其面外动态压缩性能的影响。于国际等[10]研究了二阶层级六边形蜂窝的面内动态压缩性能。赖燕辉等[11]基于多级蜂窝构型法分析了层级蜂窝的弹性模量等力学参数。

受Sierpinski三角形的启发，本文中将Sierpinski分形特性引入薄壁吸能管的层级设计，从而提出一种具有Sierpinski层级特性的新型薄壁多胞管(Sierpinski hierarchical tube,SHT)，并对SHTs在轴向冲击作用下的变形模式和吸能特性进行数值模拟分析。进一步基于能量守恒原理对SHTs在轴向压缩下的平均压缩应力理论模型进行求解，以期为新型薄壁构型轴向缓冲吸能装置的设计提供指导。

1 SHTs的层级结构设计

分形结构可以通过多种方法产生，Sierpinski[12]于1916年提出了一个典型的自相似集，并命名为Sierpinski 三角形。如图1所示，Sierpinski 垫是一个由若干个三角形构成的自相似几何结构。

基于Sierpinski[12]的研究工作，本文研究的具有Sierpinski 层级特性的多胞管是通过胞壁连接规则三角形薄壁管胞壁中点形成的。在每个SHT中有3i个边长为li的单位三角形，通过重复这一迭代过程能够得到更高阶的SHTs，i表示层级数，同时普通三角形薄壁管可以定义为零阶SHTs。第零～三阶SHTs的结构示意图如图1所示。

图1 SHT几何结构Fig.1 Geometrical structuresof SHTs

构成SHTs的单位三角形薄壁长度随着层级参数i的增大而变短，这就会给数值模拟带来很大的麻烦，因此本文中只考虑到三阶SHTs。本文中l0为90 mm，多胞管轴向长度h为250 mm。考虑3组壁厚，即t0=1.0，1.3，1.6 mm，薄壁管的相对密度分别为0.077、0.100和0.123。

由于Sierpinski层级结构特征的引入，本文中所研究的SHTs结构变得更复杂。然而，SHT 作为一种典型的薄壁结构，其仍可以被划分为若干个典型单元。如图2所示，所有的SHTs均可看作由2种基本单元组成，即V 形角单元和K 形角单元。

图2 SHTs结构角单元示意图Fig.2 Schematic angle elements of SHTs

2 SHTs的轴向压缩模拟分析

图3 SHT 轴向压缩有限元模型Fig.3 A finite element model for dynamic axial compression of an SHT

2.1 SHTs轴向压缩有限元模型

利用通用显式非线性有限元分析软件LSDYNA 模拟轴向冲击载荷作用下SHTs的动态压缩特性，计算模型如图3所示。构建上下端面刚性板，当顶端刚性板沿轴向以10 m/s的恒定速度冲击薄壁试件时，底端刚性板固定。为了准确地模拟薄壁管的大变形，胞壁采用Beltschko-Tsay 四边形壳单元，单元厚度方向采用五点积分，面内采用单点积分。基体材料为铝合金AA6060T4，材料的杨氏模量EY=68.2 GPa，屈服应力σy=80 MPa，极限应力σu=173 MPa，密度ρs=2 700 kg/m3，泊松比μ =0.3，幂指强化因数n=0.23[13]。由于铝合金对应变率不敏感，本文中未考虑应变率的影响。薄壁管胞壁采用自动单面接触算法来考虑自身变形产生的接触；薄壁管与刚性板之间采用自动点-面接触算法。数值模型中的静摩擦因数及动摩擦因数均取0.2。

2.2 有限元模型验证

为验证该模型的可靠性，首先对薄壁方形管[14]的压缩行为进行模拟。通过有限元计算后处理可以提取冲击端应力应变曲线，名义应力σ可表示为：

表1 薄壁方形管动态平均压缩力模拟结果与实验数据[14]的比较Table 1 Comparison between simulated and experimental mean dynamic compressiveforces[14]for thin-walled square tubes

模拟结果与实验数据[14]的比较见表1，误差均小于10%，表明模拟结果与实验结果吻合较好。

2.3 变形模式和塑性机理

图4给出了 ρ ¯=0.077的不同SHTs结构在轴向冲击载荷作用下不同时刻的典型变形图。由于多胞管的底部被固定，胞壁均从冲击端开始发生塑性变形，随着冲击端的不断下压，胞元逐层发生折叠变形。研究结果还表明，具有Sierpinski 层级特性的多胞管在轴向压缩过程中均呈现出轴对称渐进屈曲模式。轴对称渐进屈曲是理想的屈曲模式，已有研究表明薄壁结构在受到轴向冲击载荷作用时也多发生轴对称的向外屈曲变形[2-5]，并且塑性变形主要集中在管壁的连接处。

图4 不同SHTs 的轴向压缩变形Fig.4 Deformation of different SHTs under axial dynamic crushing

与此同时，Sierpinski 层级特征的引入导致胞元节点段数的增加。为保持结构的相对密度不变，胞壁厚度会变小，从而在某种程度上削弱了薄壁结构的弯曲变形能。与规则三角形管相比，引入Sierpinski 分级特性的多胞管变形吸能最显著的区别就是折叠单元的数量。SHTs完全折叠单元的数量会随着层级数的增加而增加。Sierpinski 分级特性的引入大大缩短了胞壁变形的半折叠波长，这也就意味着在压缩过程中更多的塑性变形能被耗散，增强了薄壁结构的抗压缩能力。

Abramowicz 等[15]指出，薄壁管在塑性变形过程中有2类基本的变形单元，即非延展性基本单元和延展性基本单元，如图5所示。图5中2H为折叠单元高度，图5(a)中延展吸能区域阴影面积为S*，图5(b)中延展吸能区域阴影面积为S#。

图5 两类基本单元Fig.5 Two kinds of basic elements

前文介绍SHTs结构时提及到作为一种典型的薄壁结构，其仍可以被划分为若干个典型单元，即如图2所示，所有的SHTs均可看作由2种基本单元组成，即V 形角单元和K 形角单元。

图6给出了2种角单元的典型变形轮廓图，通过对变形模式的仔细观察能得到SHT 结构在整个变形过程中的2种典型变形模式。这些变形模式主要取决于相邻胞壁的相对运动。对于V 形角单元，在每一个折叠单元变形过程中，2个胞壁都向同一方向移动，其变形轮廓如图6（a）所示，这属于传统的非延展性变形模式（Ⅰ型）。

对于K 形单元，4个胞壁在折叠变形过程中均向外移动，如图6（b）所示，此时相邻胞壁间的变形机理属于延展性变形模式（Ⅱ型）。

图6 SHTs变形机理Fig.6 Folding mechanisms of SHTs

3 SHTs的轴向压缩应力计算

图7为薄壁结构的胞壁受轴向冲击载荷作用时的变形示意图，其中H被称为胞元折叠单元半波长塑性铰长度，δ为轴向压缩长度，q为冲击载荷。

轴对称渐进屈曲是理想的屈曲模式，已有研究[2-5]表明薄壁结构在受到轴向冲击载荷作用时也多发生轴对称渐进屈曲变形。学者们基于简化超折叠单元（simplified super folded element,SSFE）理论求解了大量该变形模式下薄壁管的塑性坍塌问题[3-5]。

图7 薄壁结构在轴向冲击载荷下的变形示意图Fig.7 Schematic of thebuckling of a thin-walled structure under axial impact load

SHTs是一种典型的薄壁结构，通过数值模拟和分析可知其在压缩载荷作用下发生轴向渐进屈曲模式变形，这就启发我们运用SSFE理论来推导SHTs的轴向压缩力。为确定SHTs的轴向压缩应力，还需作如下假设：

(1)薄壁管材料塑性好，可视为理想刚塑性材料；

(2)如图7所示，薄壁管在轴向压缩载荷作用下发生向外的轴对称渐进屈曲，变形过程中各个折叠的塑性铰长度相等，为2H；

(3)在变形过程中，同一胞元上下表面保持平行，即胞元的各边变形量相同；

(4)压缩过程中胞壁间粘接强度足够大，不发生破裂，可忽略粘接对薄壁结构力学性能的影响。

3.1 SHTs的弯曲耗能计算

图8 基本折叠单元凸缘充分压缩示意图Fig.8 Schematic full compression of a basic folding element flange

3.2 角单元的拉伸耗散能

3.2.1 V 形角单元

图6给出了简化超折叠单元理论中的基本

折叠单元延展吸能示意图，本文中考虑了SSFE理论中的2种拉伸变形模式。图5(a)所示为变形模式，延展吸能区由3个三角形阴影部分组成，通过对阴影面积积分可求得该变形模式下单个胞壁的拉伸耗散能：

Chen 等[16]指出2个互相连接的胞壁对角单元的拉伸形变能具有相似的贡献。因此，在非延展变形的情况下，直角单元的拉伸耗散吸能可表示为单个角折叠单元拉伸耗散吸能的2倍：

Tran 等[4]在角单元的拉伸耗散吸能方面做了大量的研究，指出与直角单元相比，V 形角单元的拉伸耗散吸能较小。如图9所示，V 形角单元在完全塑性坍塌过程中的拉伸耗散吸能可表示为：

图9 直角角单元与V 形角单元之间薄膜耗散能关系Fig.9 Membrane energy relationship between the rectangular angleelement and the V-angle element

3.2.2 K 形角单元

通过对变形轮廓图的分析，不难发现K 形角单元的变形机理比直角角单元的变形机理复杂。将K 形角单元的薄膜耗散吸能简化为由一个V 形角单元和两个附加平面的薄膜耗散吸能组成。K 形角单元中的每个胞壁都具有相似的变形模式。K 形角单元中包含的V 形角单元与单个V 形角单元具有相同的抗压强度，需注意的是根据仿真结果（图7）可知此时的V 形角单元处于延展性变形模式，如图6（b）所示。阴影部分为延展吸能区，通过积分阴影区域来评估完全塌陷期间每个胞壁的薄膜耗散吸能：

图10给出了在K 形角单元的拐角处形成的延展吸能三角形单元。K 形角单元中包含的V 形角单元的薄膜耗散能可表示为：

图10 变形模式Fig.10 Deformation modes

3.3 SHTs的动态平均压缩应力计算

普通三角形薄壁管，即零阶SHTs是由3个V 形角单元组合而成。将式（12）和（15）代入式（9），平均压缩力Pm,0th的理论计算方程可写为：

3.4 理论模型模拟验证

本文中t0分别取1.0、1.3和1.6 mm，胞壁长l0取90 mm，分别对这些新型SHT多胞结构进行轴向压缩的数值模拟。

用式(29)～(32)计算出动态压缩应力的理论值，并将其与模拟结果进行比较，如表2所示，当t0=1.0mm（ ρ ¯=0.077）时，一阶、二阶及三阶SHTs的动态压缩应力较普通三角形薄壁管分别升高了85.8%、138.2%和183.8%。这表明：将Sierpinski 层级特性应用到薄壁管的设计中，可以有效提高薄壁管的吸能特性。4种薄壁管的动态压缩应力模拟值和理论值的差异分别为-2.64%～3.42%、-5.23%～1.54%、-2.77%～3.14%和-1.19%～3.67%，结果基本吻合，进一步验证了本文动态压缩应力理论推导方法的可行性，说明理论计算结果具有工程应用价值，可用于指导新型薄壁结构轴向缓冲吸能装置的设计。

4 结论

将Sierpinski分形结构引入到薄壁吸能管的层级设计中，提出了一种具有Sierpinski 层级特性的新型薄壁多胞管（SHT），并通过数值模拟和理论方法研究其受轴向冲击载荷作用下的变形模式和能量吸收特性，研究结果表明：

（1）具有Sierpinski层级特性的SHTs在轴向压缩过程中均呈现出理想的轴对称渐进屈曲模式。SHTs完全折叠单元的数量会随着层级数的增加而增加。Sierpinski 分级特性的引入大大缩短了胞壁变形的半折叠波长，压缩过程中更多的塑性变形能量被耗散，增强了薄壁结构的抗压缩能力。

（2）SHTs均可看作由两种基本单元组成，即V 形角单元和K 形角单元。V 形角单元在每一个折叠单元变形过程中，两个胞壁都向同一方向移动，属于传统的非延展性变形模式（Ⅰ型）。K 形角单元4个胞壁在折叠变形过程中均向外移动，相邻胞壁间的变形机理属于延展性变形模式（Ⅱ型）。

（3）基于能量守恒理论和塑性铰理论建立了SHTs的塑性屈曲理论模型，获得了压溃过程中的塑性弯曲耗散能、薄膜耗散能、轴向压缩应力的理论表达式，模拟结果与理论计算结果基本吻合，验证了本文中动态压缩应力理论推导方法是可行的，理论计算结果具有工程应用价值。

（4）在相同的相对密度下，一阶、二阶及三阶SHTs的动态压缩应力较普通三角形薄壁管分别增高了85.8%、138.2%和183.8%。将Sierpinski 层级特性引入到薄壁管的设计中，能够有效提高薄壁管的耐撞性能，这可为新型吸能元件的研究和设计提供参考。