CN105117542A - 一种计算四边手性新型蜂窝轴向压缩应力的方法 - Google Patents

Abstract

本发明提供一种计算四边手性新型蜂窝轴向压缩应力的方法，包括以下步骤：步骤1、基于四边手性胞元蜂窝的几何构型，计算一个叠缩单元所吸收的总能量；步骤2、计算四个附加平面单元所吸收的总能量；步骤3、基于步骤1、2所计算得到的叠缩单元所吸收的总能量以及四个附加平面单元所吸收的总能量，根据能量守恒原理计算四边手性胞元蜂窝受轴向压缩下的轴向压缩应力。本发明的方法中，基于能量守恒定理，结合超折叠单元理论，给出了一种计算其轴向压缩应力的理论方法，并且运用有限元仿真技术对四边手性胞元蜂窝受轴向压缩进行模拟，从而验证理论计算方法的正确性。

Description

一种计算四边手性新型蜂窝轴向压缩应力的方法

技术领域

本发明专利涉及四边手性胞元蜂窝材料，具体讲的是一种计算四边手性新型蜂窝轴向压缩应力的方法。

背景技术

由于蜂窝材料质量轻并且具有良好的能量吸收特性，人们常用其制作各种能量吸收或缓冲装置，并且广泛应用于航天航空及车辆工程领域。蜂窝材料在受到轴向冲击作用时，能够吸收大量的冲击能量。随着蜂窝材料在防护结构中的应用日益广泛，冲击过程中蜂窝材料轴向压缩特性的研究已成为其性能设计的关键问题之一。

蜂窝结构的轴向冲击力学性能与其几何构型密切相关。为提高传统构型(正六边形、圆形、三角形、正方形等)蜂窝的吸能能力，国内外相关研究提出一种新型的双筋加强正六边形蜂窝。但当前的研究工作主要集中在预测其准静态力学性能及解释其变形机制上，对实际应用中其用作吸能装置材料时受到外部轴向冲击的动态力学特性尚缺乏了解。

发明内容

本发明目的在于提供一种计算四边手性新型蜂窝轴向压缩应力的方法.

本发明的上述目的通过独立权利要求的技术特征实现，从属权利要求以另选或有利的方式发展独立权利要求的技术特征。

为达成上述目的，本发明提出一种计算四边手性新型蜂窝轴向压缩应力的方法，包括：

步骤1、基于四边手性胞元蜂窝的几何构型，在四边手性胞元蜂窝的圆管部分受轴向压缩发生非轴对称屈曲时，计算一个叠缩单元所吸收的总能量；

步骤2、计算四个附加平面单元所吸收的总能量；

步骤3、基于步骤1、2所计算得到的叠缩单元所吸收的总能量以及四个附加平面单元所吸收的总能量，根据能量守恒原理计算四边手性胞元蜂窝受轴向压缩下的轴向压缩应力。

本发明与现有技术相比，其显著优点在于：本发明的计算方法，不局限于预测其准静态力学性能及解释其变形机制上，基于能量守恒原理，运用超折叠单元理论计算其轴向压缩应力。该计算方法具有解析表达式简单、精度高、可靠性高的优点，可以快速判断不同胞元尺寸的蜂窝吸能装置的吸能能力，为蜂窝吸能装置的设计提供有力的理论指导和技术支持，具有很好的工程价值及应用前景。

应当理解，前述构思以及在下面更加详细地描述的额外构思的所有组合只要在这样的构思不相互矛盾的情况下都可以被视为本公开的发明主题的一部分。另外，所要求保护的主题的所有组合都被视为本公开的发明主题的一部分。

结合附图从下面的描述中可以更加全面地理解本发明教导的前述和其他方面、实施例和特征。本发明的其他附加方面例如示例性实施方式的特征和/或有益效果将在下面的描述中显见，或通过根据本发明教导的具体实施方式的实践中得知。

附图说明

附图不意在按比例绘制。在附图中，在各个图中示出的每个相同或近似相同的组成部分可以用相同的标号表示。为了清晰起见，在每个图中，并非每个组成部分均被标记。现在，将通过例子并参考附图来描述本发明的各个方面的实施例，其中：

图1为四边手性胞元蜂窝几何构型的示意图。

图2为四边手性胞元蜂窝变形图。

图3为理想轴对称压缩变形模式示意图。

图4为环向屈曲变形简化示意图。

图5为简化的超折叠单元的整体变形图。

图6为简化的超折叠单元的延展单元与静态塑性铰线示意图。

具体实施方式

为了更了解本发明的技术内容，特举具体实施例并配合所附图式说明如下。

在本公开中参照附图来描述本发明的各方面，附图中示出了许多说明的实施例。本公开的实施例不必定意在包括本发明的所有方面。应当理解，上面介绍的多种构思和实施例，以及下面更加详细地描述的那些构思和实施方式可以以很多方式中任意一种来实施，这是因为本发明所公开的构思和实施例并不限于任何实施方式。另外，本发明公开的一些方面可以单独使用，或者与本发明公开的其他方面的任何适当组合来使用。

根据本发明的实施例，一种计算四边手性新型蜂窝轴向压缩应力的方法，总体上包括以下三个步骤：

步骤1、基于四边手性胞元蜂窝的几何构型，在四边手性胞元蜂窝的圆管部分受轴向压缩发生非轴对称屈曲时，计算一个叠缩单元所吸收的总能量；

步骤2、计算四个附加平面单元所吸收的总能量；

步骤3、基于步骤1、2所计算得到的叠缩单元所吸收的总能量以及四个附加平面单元所吸收的总能量，根据能量守恒原理计算四边手性胞元蜂窝受轴向压缩下的轴向压缩应力。

下面结合图1～图6所示，更加详细地描述前述各个步骤的具体实现。

如图1所示，对于四边手性胞元蜂窝结构，其几何构型如图1所示。图2给出了该蜂窝结构的变形模式图。四边手性胞元蜂窝的基本折叠单元可以看做是一个圆柱管和四个附加平面单元的组合。

[步骤1]

基于四边手性胞元蜂窝的几何构型，在四边手性胞元蜂窝的圆管部分受轴向压缩发生非轴对称屈曲时，计算一个叠缩单元所吸收的总能量

如图2所示，在四边手性胞元蜂窝变形过程中，圆柱管表现为轴向与环向同时出现波纹的非轴对称屈曲模式，而且屈曲并不是无规律的，它以一定的叠缩单元重复出现，只是每个叠缩单元之间有一定的扭转角。这种屈曲模式除发生轴向屈曲外，在环向截面也会由对称屈曲模式的圆形变为规则的三角形形状。一个叠缩单元变形所耗散的能量主要由三个部分组成，即轴向弯曲变形能E₁、环向弯曲变形能E₂以及伸张变形能E₃。

圆柱管压缩发生非轴对称屈曲大变形时管壁将沿三角形各边轴向折叠在一起，不考虑弹性变形影响，轴向弯曲变形能全部由各边上的塑性铰弯曲变形能组成。变形模式与轴对称屈曲相同，如图3，从而可得轴向弯曲变形能E₁表达为：

$E_1=4M_0\·\frac{\mathrm{\π}}{2}\·P---\left(1\right)$

其中t为胞元壁厚，M₀＝σ₀t²/4为管壁沿三角形各折叠边单位长度的弯曲塑性极限弯矩，其塑性流动应力σ₀取为P为环向塑性铰的总长度，即为三角形的周长，忽略变形时管壁的影响，假设变形是理想的，即变形后三角形的周长应该等于变形前圆管的周长，即：

P＝2πR(2)

其中R为圆管半径，将式(2)代入(1)则可得：

E₁＝σ₀π²Rt²(3)

图4为环向屈曲变形简化示意图，图中内圆表示为压缩变形前的圆管，外圆表示变形后三角形折叠单元的外接圆。三角形的周长和未变形的圆管周长相等，根据其几何特征关系，即：

通过上式可得半折叠波长：

$H=\frac{\sqrt{3}\mathrm{\π}R}{9}---\left(5\right)$

管壁发生环向弯曲屈曲时，圆管管壁平均分成三段圆弧，三段圆弧逐渐伸展成为三角形的三条边，并且每两段圆弧之间的夹角由原来的π逐渐弯曲折叠为π/3，此过程中消耗的能量即为弯曲变形能E₂，同时由于圆弧段长度大于直线长度，三角形的三个顶点将在展平过程中向外移动，如图4所示。环向屈曲变形能由三角形顶点处的塑性铰弯曲变形产生，则每个塑性铰的长度为忽略弹性变形影响，则每个折叠单元产生的环向屈曲变形能为：

$E_2=6\·M_0\·\frac{\mathrm{\π}}{3}\·\frac{\sqrt{3}\mathrm{\π}R}{9}=\frac{\sqrt{3}}{18}{\mathrm{\σ}}_0{\mathrm{\π}}^2{\mathrm{Rt}}^2---\left(6\right)$

由图2有限元模型模拟结果可以看出，塑性铰部分圆管的轴向屈曲变形并非按图3所示的形式变形，而是沿轴向被压缩，并向外凸起，压缩距离与图4所示棱边部分相同，凸出位移等于图4中顶点向外的径向位移。则轴向伸张变形能主要由圆管压缩变形所耗散的能量组成。

假设图5中圆管被压缩一个微小的位移，折叠胞壁与竖直方向夹角α也会产生一个微小增量dα，被压缩部分管壁平均应变为：

$\frac{2Hcos\mathrm{\α}-2Hcos(\mathrm{\α}+d\mathrm{\α})}{2H\cos \mathrm{\α}}=tan\mathrm{\α}\·\sin d\mathrm{\α}\≈tan\mathrm{\α}\·d\mathrm{\α}---\left(7\right)$

则轴向伸张变形能为：

$E_3=3\∫{\mathrm{dE}}_3=3{\∫}_0^{\mathrm{\π}/2}{\mathrm{\σ}}_{0}tan\mathrm{\α}\·\mathrm{\λ}t\·2Hcos\mathrm{\α}\·d\mathrm{\α}=\frac{8\sqrt{3}}{3}{\mathrm{\σ}}_0{\mathrm{\πRt}}^2---\left(8\right)$

式中：λ为塑性铰环向长度，取λ＝4t。

因此圆管部分受轴向压缩发生非轴对称屈曲时，一个叠缩单元所吸收的总能量为：

$E_{circular}=E_1+E_2+E_3={\mathrm{\σ}}_0{\mathrm{\π}}^2{\mathrm{Rt}}^2+\frac{\sqrt{3}}{18}{\mathrm{\σ}}_0{\mathrm{\π}}^2{\mathrm{Rt}}^2+\frac{8\sqrt{3}}{3}{\mathrm{\σ}}_0{\mathrm{\πRt}}^2=(\frac{\sqrt{3}+18}{18}\mathrm{\π}+\frac{8\sqrt{3}}{3}){\mathrm{\σ}}_0{\mathrm{\πRt}}^2---\left(9\right)$

[步骤2]

计算四个附加平面单元所吸收的总能量

单个附加平面单元的变形可以看作是由3个延展性三角形单元和3条静态塑铰线组成的一块翼缘板。如图5所示为简化的超折叠单元的整体变形图。图6为简化超折叠单元的延展单元及塑性铰线。

四边手性胞元蜂窝可以看做是多个超折叠单元的集合，而每个超折叠单元的整体变形包括附加平面单元的变形以及(蜂窝的)圆柱形的变形。

如图6所示，超折叠单元在变形后，采用阴影区表示在其角线附近形成的3个延展性单元，上下两个三角形为受压单元，中间的三角形为受拉单元。三条静态塑铰线分别位于板的上部、中部和下部，相应的旋转角分别为θ，2θ和θ。

其中弯曲变形能E_b可通过累计各条静态塑铰线处的能量耗散求得。对于单个翼缘板共有三条塑性铰，因此其弯曲变形能表示为：

$E_b=\underset{i=1}{\overset{3}{\Σ}}{M}_0{\mathrm{\θ}}_{i}C---\left(10\right)$

其中C为每条塑性铰的长度，θ_i为每条塑性铰的旋转角度。为了简化起见，假使轴向压缩距离为2H，则单个翼缘板被完全压平，此时三条塑性铰线的旋转角度分别为π和因此有

E_b＝2πM₀C(11)

单个波长压缩范围内所耗散的薄膜变形能E_m可以通过对拉伸和压缩区域的面积积分求得，即图6中阴影部分的面积。

$E_m={\∫}_s{\mathrm{\σ}}_{0}tds=\frac{1}{2}{\mathrm{\σ}}_0{\mathrm{tH}}^2---\left(12\right)$

根据公式(10)，四个附加平面基本折叠单元的塑性铰耗散能量可通过下式计算：

$E_{b-TC}=4\×\frac{\mathrm{\π}}{2}\×4\×\frac{L}{2}\·\frac{{\mathrm{\σ}}_0{t}^2}{4}={\mathrm{\π\σ}}_0{t}^{2}L---\left(13\right)$

L为胞元胞壁长度。假定每一个附加平面有着相同的变形模式，则四个附加平面基本折叠单元总薄膜变形能为：

E_m-TC＝4E_m＝2σ₀tH²(14)

将公式(5)代入上式，可得

$E_{m-TC}=\frac{2}{27}{\mathrm{\σ}}_0{\mathrm{\π}}^2{R}^{2}t---\left(15\right)$

[步骤3]

基于步骤1、2所计算得到的叠缩单元所吸收的总能量以及四个附加平面单元所吸收的总能量，根据能量守恒原理计算四边手性胞元蜂窝受轴向压缩下的轴向压缩应力

当四边手性胞元蜂窝受到轴向压缩时，根据能量守恒原理可得：

$q_{c-TC}{\mathrm{\δ}}_{e}2H=E_{b-TC}+E_{m-TC}+E_{circular}={\mathrm{\π\σ}}_0{t}^{2}L+\frac{2}{27}{\mathrm{\σ}}_0{\mathrm{\π}}^2{R}^{2}t+(\frac{\sqrt{3}+18}{18}\mathrm{\π}+\frac{8\sqrt{3}}{3}){\mathrm{\σ}}_0{\mathrm{\πRt}}^2---\left(16\right)$

将公式(5)以及R＝L/4代入上式，化简可得：

$q_{c-TC}=\frac{24}{\sqrt{3}\mathrm{\π}L}\[\frac{1}{216}{\mathrm{\σ}}_0{\mathrm{\π}}^2{L}^{2}t+(\frac{\sqrt{3}+18}{72}\mathrm{\π}+\frac{3+2\sqrt{3}}{3}){\mathrm{\σ}}_0{\mathrm{\πLt}}^2\]=\frac{\sqrt{3}}{27}{\mathrm{\σ}}_0\mathrm{\π}Lt+(\frac{6\sqrt{3}+1}{3}\mathrm{\π}+16+8\sqrt{3}){\mathrm{\σ}}_0{t}^2---\left(17\right)$

而四边手性胞元蜂窝的基本折叠单元所占的面积S(图1中阴影部分面积)为：

$S=\frac{5}{4}{L}^2---\left(18\right)$

所以四边手性胞元蜂窝的轴向准静态平均压缩应力为：

${\mathrm{\σ}}_{m-TC}=q_{c-TC}/S=\[\frac{4(6\sqrt{3}+1)}{15}\mathrm{\π}+\frac{64}{5}+\frac{32\sqrt{3}}{5}\]{\mathrm{\σ}}_0{(t/L)}^2+\frac{4\sqrt{3}}{135L}{\mathrm{\σ}}_0\mathrm{\π}t---\left(19\right)$

下面结合一些具体的实例，对前述方法得到的四边手性胞元蜂窝轴向压缩应力理论公式进行验证和说明。

采用非线性显式有限元软件ANSYS/LS-DYNA对轴向压缩下蜂窝结构的变形过程进行数值分析。在建立计算模型时，利用ANSYS参数化设计语言(APDL)生成蜂窝芯或蜂窝夹层板的有限元模型，并运用LS-DYNA970进行求解。

蜂窝夹层板置于两块刚性平面之间，作为压缩过程中的支撑平台，下端的刚性板(RIGID-WALL-GEOMETRIC-FLAT)完全固定，上端的刚性板(RIGID-WALL-GEOMETRIC-FLAT-MOTION)以0.5mm/s的速度向下压蜂窝结构。

计算中所有薄壁均用Belytschko-Tsay4节点四边形壳单元进行离散，沿厚度方向采用5个积分点，面内采用1个积分点，且假定蜂窝胞元之间的黏结不存在失效。仿真分析时，定义了蜂窝结构变形过程中自身结构的自动单面接触算法(Automaticsingle-surfacecontactalgorithm)以及刚性板与蜂窝材料之间的自动点面接触(Automaticnode-to-surfacecontact)，摩擦因子为0.2。

蜂窝材料采用铝合金AA6060T4，材料的力学性能参数为：杨氏模量E＝68.2GPa，屈服应力σ_y＝80MPa，极限应力σ_u＝173MPa，密度ρ＝2700kg/m3，泊松比υ＝0.3，幂指强化系数n＝0.23。采用LS-DYNA里的#123号材料模型(Modifiedpiece-wisedmaterialmodel)对蜂窝材料进行分析。

对于四边手性胞元蜂窝结构，其几何构型如图1所示。图2给出了该蜂窝结构的变形模式图。四边手性胞元蜂窝的基本折叠单元可以看做是一个圆柱管和四个附加平面单元的组合。

为验证理论公式的正确性，再运用Matlab软件随机生成20组结构参数，结构参数如表1所示。

表1为四边手性胞元蜂窝所有样本点平均压缩应力比较

运用公式(19)对四边手性胞元蜂窝轴向压缩进行了计算，并将理论值与相应的平台压缩应力有限元仿真值进行比较。其平均压缩应力的理论值和仿真值的比较如表1所示，误差在-5.44％～4.84％之间，属于可以接受的误差范围内，所以理论计算结果具有工程应用价值，可用于指导新型四边手性胞元蜂窝结构的轴向缓冲吸能装置设计。

虽然本发明已以较佳实施例揭露如上，然其并非用以限定本发明。本发明所属技术领域中具有通常知识者，在不脱离本发明的精神和范围内，当可作各种的更动与润饰。因此，本发明的保护范围当视权利要求书所界定者为准。

Claims (4)

1.一种计算四边手性新型蜂窝轴向压缩应力的方法，其特征在于，包括以下步骤：

步骤1、基于四边手性胞元蜂窝的几何构型，在四边手性胞元蜂窝的圆管部分受轴向压缩发生非轴对称屈曲时，计算一个叠缩单元所吸收的总能量；

步骤2、计算四个附加平面单元所吸收的总能量；

步骤3、基于步骤1、2所计算得到的叠缩单元所吸收的总能量以及四个附加平面单元所吸收的总能量，根据能量守恒原理计算四边手性胞元蜂窝受轴向压缩下的轴向压缩应力。

下面结合图1～图6所示，更加详细地描述前述各个步骤的具体实现。

2.根据权利要求1所述的计算四边手性新型蜂窝轴向压缩应力的方法，其特征在于，所述步骤1的具体实现包括：

将四边手性胞元蜂窝的基本折叠单元看做是一个圆柱管和四个附加平面单元的组合，在四边手性胞元蜂窝变形过程中，圆柱管表现为轴向与环向同时出现波纹的非轴对称屈曲模式，且屈曲以一定的叠缩单元重复出现，一个叠缩单元变形所耗散的能量由三个部分组成，即轴向弯曲变形能E₁、环向弯曲变形能E₂以及伸张变形能E₃，其中：

圆柱管压缩发生非轴对称屈曲大变形时管壁将沿三角形各边轴向折叠在一起，不考虑弹性变形影响，轴向弯曲变形能全部由各边上的塑性铰弯曲变形能组成，变形模式与轴对称屈曲相同，从而可得轴向弯曲变形能E₁表达为：

$E_1=4M_0\·\frac{\mathrm{\π}}{2}\·P---\left(1\right)$

其中，t为胞元壁厚，M₀＝σ₀t²/4为管壁沿三角形各折叠边单位长度的弯曲塑性极限弯矩，其塑性流动应力σ₀取为P为塑性铰的总长度，即为三角形的周长，忽略变形时管壁的影响，假设变形是理想的，即变形后三角形的周长应该等于变形前圆管的周长，即：

P＝2πR(2)

其中R为圆管半径，将式(2)代入(1)则可得：

E₁＝σ₀π²Rt²(3)

由于三角形的周长和未变形的圆管周长相等，根据其几何特征关系，即：

通过上式可得半折叠波长：

$H=\frac{\sqrt{3}\mathrm{\π}R}{9}---\left(5\right)$

管壁发生环向弯曲屈曲时，圆管管壁平均分成三段圆弧，三段圆弧逐渐伸展成为三角形的三条边，并且每两段圆弧之间的夹角由原来的π逐渐弯曲折叠为π/3，此过程中消耗的能量即为弯曲变形能E₂，同时由于圆弧段长度大于直线长度，三角形的三个顶点将在展平过程中向外移动，环向屈曲变形能由三角形顶点处的塑性铰弯曲变形产生，则每个塑性铰的长度为忽略弹性变形影响，则每个折叠单元产生的环向屈曲变形能为：

$E_2=6\·M_0\·\frac{\mathrm{\π}}{3}\·\frac{\sqrt{3}\mathrm{\π}R}{9}=\frac{\sqrt{3}}{18}{\mathrm{\σ}}_0{\mathrm{\π}}^2{\mathrm{Rt}}^2---\left(6\right)$

假设圆管被压缩一个微小的位移，折叠胞壁与竖直方向夹角α也会产生一个微小增量dα，被压缩部分管壁平均应变为：

$\frac{2Hcos\mathrm{\α}-2Hcos(\mathrm{\α}+d\mathrm{\α})}{2H\cos \mathrm{\α}}=tan\mathrm{\α}\·\sin d\mathrm{\α}\≈tan\mathrm{\α}\·d\mathrm{\α}---\left(7\right)$

则轴向伸张变形能为：

$E_3=3\∫{\mathrm{dE}}_3=3{\∫}_0^{\mathrm{\π}/2}{\mathrm{\σ}}_{0}tan\mathrm{\α}\·\mathrm{\λ}t\·2Hcos\mathrm{\α}\·d\mathrm{\α}=\frac{8\sqrt{3}}{3}{\mathrm{\σ}}_0{\mathrm{\πRt}}^2---\left(8\right)$

式中：λ为塑性铰环向长度，λ＝4t；

因此，圆管部分受轴向压缩发生非轴对称屈曲时，一个叠缩单元所吸收的总能量为：

$E_{circular}=E_1+E_2+E_3={\mathrm{\σ}}_0{\mathrm{\π}}^2{\mathrm{Rt}}^2+\frac{\sqrt{3}}{18}{\mathrm{\σ}}_0{\mathrm{\π}}^2{\mathrm{Rt}}^2+\frac{8\sqrt{3}}{3}{\mathrm{\σ}}_0{\mathrm{\πRt}}^2=(\frac{\sqrt{3}+18}{18}\mathrm{\π}+\frac{8\sqrt{3}}{3}){\mathrm{\σ}}_0{\mathrm{\πRt}}^2---\left(9\right)$

3.根据权利要求2所述的计算四边手性新型蜂窝轴向压缩应力的方法，其特征在于，所述步骤2的具体实现包括：

将单个附加平面单元的变形看作是由3个延展性三角形单元和3条静态塑铰线组成的一块翼缘板，四边手性新型蜂窝的超折叠单元在变形后，采用阴影区表示在其角线附近形成的3个延展性单元，上下两个三角形为受压单元，中间的三角形为受拉单元，三条静态塑铰线分别位于板的上部、中部和下部，相应的旋转角分别为θ，2θ和θ，

其中弯曲变形能E_b通过累计各条静态塑铰线处的能量耗散求得；

对于单个翼缘板共有三条塑性铰，因此总的弯曲变形能表示为：

$E_b=\underset{i=1}{\overset{3}{\Σ}}{M}_0{\mathrm{\θ}}_{i}C---\left(10\right)$

其中C为每条塑性铰的长度，θ_i为每条塑性铰的旋转角度，假使轴向压缩距离为2H，则单个翼缘板被完全压平，此时三条塑性铰线的旋转角度分别为π和因此有

E_b＝2πM₀C(11)

单个波长压缩范围内所耗散的薄膜变形能E_m通过对拉伸和压缩区域的面积积分求得：

$E_m={\∫}_s{\mathrm{\σ}}_{0}tds=\frac{1}{2}{\mathrm{\σ}}_0{\mathrm{tH}}^2---\left(12\right)$

根据公式(10)，四个附加平面基本折叠单元的塑性铰耗散能量可通过下式计算：

$E_{b-TC}=4\×\frac{\mathrm{\π}}{2}\×4\×\frac{L}{2}\·\frac{{\mathrm{\σ}}_0{t}^2}{4}={\mathrm{\π\σ}}_0{t}^{2}L---\left(13\right)$

L为胞元胞壁长度；假定每一个附加平面有着相同的变形模式，则四个附加平面基本折叠单元总薄膜变形能为：

E_m-TC＝4E_m＝2σ₀tH²(14)

将公式(5)代入上式，可得：

$E_{m-TC}=\frac{2}{27}{\mathrm{\σ}}_0{\mathrm{\π}}^2{R}^{2}t.---\left(15\right)$

4.根据权利要求3所述的计算四边手性新型蜂窝轴向压缩应力的方法，其特征在于，所述步骤3的具体实现包括：

当四边手性胞元蜂窝受到轴向压缩时，根据能量守恒原理可得：

$q_{c-TC}{\mathrm{\δ}}_{e}2H=E_{b-TC}+E_{m-TC}+E_{circular}={\mathrm{\π\σ}}_0{t}^{2}L+\frac{2}{27}{\mathrm{\σ}}_0{\mathrm{\π}}^2{R}^{2}t+(\frac{\sqrt{3}+18}{18}\mathrm{\π}+\frac{8\sqrt{3}}{3}){\mathrm{\σ}}_0{\mathrm{\πRt}}^2---\left(16\right)$

将公式(5)以及R＝L/4代入上式，化简可得：

$q_{c-TC}=\frac{24}{\sqrt{3}\mathrm{\π}L}\[\frac{1}{216}{\mathrm{\σ}}_0{\mathrm{\π}}^2{L}^{2}t+(\frac{\sqrt{3}+18}{72}\mathrm{\π}+\frac{3+2\sqrt{3}}{3}){\mathrm{\σ}}_0{\mathrm{\πLt}}^2\]=\frac{\sqrt{3}}{27}{\mathrm{\σ}}_0\mathrm{\π}Lt+(\frac{6\sqrt{3}+1}{3}\mathrm{\π}+16+8\sqrt{3}){\mathrm{\σ}}_0{t}^2---\left(17\right)$

而四边手性胞元蜂窝的基本折叠单元所占的面积S为：

$S=\frac{5}{4}{L}^2---\left(18\right)$

所以四边手性胞元蜂窝的轴向准静态平均压缩应力为：

${\mathrm{\σ}}_{m-TC}=q_{c-TC}/S=\[\frac{4(6\sqrt{3}+1)}{15}\mathrm{\π}+\frac{64}{5}+\frac{32\sqrt{3}}{5}\]{\mathrm{\σ}}_0{(t/L)}^2+\frac{4\sqrt{3}}{135L}{\mathrm{\σ}}_0\mathrm{\π}t---\left(19\right)$

该公式(19)的σ_m-TC即为四边手性胞元蜂窝受轴向压缩下的轴向压缩应力计算公式，根据此公式即可计算得出四边手性胞元蜂窝受轴向压缩下的轴向压缩应力。