星形蜂窝结构本构关系的建立方法、介质和设备

本发明涉及星形蜂窝结构技术领域，特别涉及一种星形蜂窝结构本构关系的建立方法、介质和设备。

背景技术：

负泊松比(negativepoisson’sratio,npr)多胞元材料在在拉伸时横向膨胀，压缩时横向收缩，由于其独特的性能，使负泊松比多胞材料表现出优异的力学性能，比如轻质、比强度高、隔音隔热性能好，吸热效率高的特点，在航空航天，车辆舰船，包装工程等领域已有广泛应用。近年来，一种负泊松比星形节点周期性蜂窝结构，因其受到荷载作用而表现出独特的变形特征引起了相关学者的关注。但要在工程上应用该类新型周期性结构，首先要对其宏微观力学性能进行充分的研究和认识。

在工程设计和分析中，星形蜂窝结构复杂的几何形式和各项异性的力学性能给学者带来很大的苦难，通常的试验方法耗时费力又不经济。随着计算科学的发展，采用通用有限元程序对大型复杂结构进行结构建模及力学仿真已成为现实，然后，在结构系统建模过程中，详细的建模又大大增加模型规模和难度，使得计算效率显著降低。同时，复杂多样的结构及边界条件也使得模型具有一定的局限性。因此，工程技术人员往往希望将非连续的星形结构等效替代为连续的正交各向异性体，从而用材料等效模型和板桥单元来模拟星形结构，达到降低建模难度并提高计算效率。

技术实现要素：

本发明的第一目的在于克服现有技术的缺点与不足，提供一种星形蜂窝结构本构关系的建立方法，该方法可以针对结构的实际变形、几何特征及材料杨氏模量，准确方便地建立出反映同一类负泊松比星形节点周期性蜂窝结构的本构关系。

本发明的第二目的在于提供一种计算机可读存储介质。

本发明的第三目的在于提供一种计算设备。

本发明的第四目的在于提供一种星形蜂窝结构本构关系的建立装置。

本发明的第一目的通过下述技术方案实现：一种星形蜂窝结构本构关系的建立方法，包括如下步骤：

s1、从由一系列星形胞元组成的星形蜂窝结构中切出一个星形胞元；

s2、当沿着水平方向或竖直方向单向拉伸星形胞元时，考虑到星形胞元具有对称性，对于星形胞元的研究，直接取附加边界约束的四分之一胞元进行力学分析；

考虑四分之一胞元受力特性以及边界约束，将四分之一胞元的力学分析过程进一步简化为斜壁板的受力变形；然后根据斜壁板的受力变形、星形胞元结构的几何特征及材料杨氏模量，运用能量法求解出胞元结构的等效杨氏模量以及等效泊松比；

s3、当沿着与水平方向成角方向单向拉伸星形胞元时，将星形胞元的力学分析过程简化为星形胞元在水平方向受力和竖直方向受力的叠加，通过求解在x、y方向上的受力变形，即可进一步计算出胞元结构的等效杨氏模量以及等效泊松比。

优选的，四分之一胞元具有相连的两个斜壁板，一个自其中一个斜壁板水平延伸以连接其他星形胞元的水平壁板，以及一个自另一个斜壁板竖直延伸以连接其他星形胞元的竖直壁板；

斜壁板的受力变形分析及求解等效杨氏模量、等效泊松比的过程具体如下：

先求出斜壁板的弯矩方程、轴力方程，以及弯矩、轴力对斜壁板受力的偏导数；

然后基于弯矩方程、轴力方程及其对拉力的偏导数，运用能量法求得斜壁板在水平方向和竖直方向的位移，进而求得四分之一胞元在水平方向和竖直方向的位移，根据四分之一胞元在水平方向的位移计算水平方向的应变，根据四分之一胞元在竖直方向的位移计算竖直方向的应变；

将水平和竖直方向的应变代入杨氏模型计算公式和泊松比计算公式，即可得出四分之一胞元的等效杨氏模量以及等效泊松比，星形胞元的等效杨氏模量以及等效泊松比与四分之一胞元的等效杨氏模量以及等效泊松比相等。

更进一步的，当星形胞元在水平方向上单向受拉时，定义x方向代表水平方向，y方向代表竖直方向，z方向代表与x方向和y方向相垂直的方向，连接水平壁板的斜壁板为第一斜壁板，连接竖直壁板的斜壁板为第二斜壁板，星形胞元的等效泊松比vxy和等效杨氏模量ex求解过程如下：

s211、四分之一胞元在x方向上的有效横截面积ax、总拉伸力fx表示为：

式中，l2是第一斜壁板的长度；l3是第二斜壁板的长度；l4是竖直壁板的长度；θ1是水平壁板与第一斜壁板的夹角；θ2是第二斜壁板与竖直壁板的夹角；b是胞元壁板的厚度；σx为四分之一胞元在水平方向上的应力；t为胞元结构在z方向上的深度；

s212、忽略由于剪切胞元引起的挠曲变形，第一斜壁板的轴力方程和弯矩方程及其对各力的偏导数表示为：

其中，fbc为第一斜壁板的轴力；mbc(x1)为第一斜壁板上任一点的弯矩；fb为第一斜壁板与水平壁板的共同端点上的虚加力，fb＝0；x1为第一斜壁板上任意一点至第一斜壁板与水平壁板共同端点的长度；fx表示第一斜壁板在x方向上受到的力；

第二斜壁板的轴力方程和弯矩方程及其对各力的偏导数表示为：

其中，fcd为第二斜壁板的轴力；mcd(x2)为第二斜壁板上任一点的弯矩；x2为第二斜壁板上任意一点至第二斜壁板与第一斜壁板共同端点的长度；

s213、水平壁板中用于连接其他星形胞元的端点在y方向上的位移等于水平壁板与第一斜壁板的共同端点在y方向上的位移，水平壁板中用于连接其他星形胞元的端点在x方向上的位移等于水平壁板拉伸长度加上水平壁板与第一斜壁板的共同端点在x方向上的位移，由能量法计算出水平壁板中用于连接其他星形胞元的端点在水平方向上的位移δax和在竖直方向上的位移δay，表示为：

其中，es为胞元结构的材料杨氏模量；i为胞元结构的极惯性矩；a为胞元壁板的横截面积；

s214、计算出四分之一胞元在x方向上的应变εx′、在y方向上的应变εy′，表示为：

其中，l1为水平壁板的长度；

s215、四分之一胞元的等效泊松比vxy和等效杨氏模量ex表示为：

更进一步的，随着壁板夹角θ1、θ2增大，星形蜂窝结构在x方向的杨氏模量ex减小，在x方向的泊松比νxy增大；

随着水平壁板长度l1增大，杨氏模量ex增大，泊松比νxy减小；

随着第一斜壁板长度l2增大，杨氏模量ex减小，泊松比νxy先减小后增大；

随着第二斜壁板长度l3增大，杨氏模量ex增大，泊松比νxy减小；

随着竖直壁板长度l4增大，杨氏模量ex减小，泊松比νxy增大。

更进一步的，对于星形胞元胞元，当其在竖直方向上单向受拉时，等效泊松比vyx和等效杨氏模量ey的求解过程如下：

s221、四分之一胞元在y方向上的有效横截面积ay、总拉伸力fy表示为：

式中，σy为四分之一胞元在竖直方向上的应力；

s222、第一斜壁板、第二斜壁板的轴力方程和弯矩方程及其对各力的偏导数表示为：

其中，f′bc为第一斜壁板的轴力；f′cd为第二斜壁板的轴力；m′cd(x4)为第二斜壁板上任一点的弯矩；m′bc(x3)为第一斜壁板上任一点的弯矩；x3为第一斜壁板上任意一点至第一斜壁板与第二斜壁板共同端点的长度；x4为第二斜壁板上任意一点至第二斜壁板与垂直壁板共同端点的长度；fd为第二斜壁板与竖直壁板的共同端点上的虚加力，fd＝0；fy表示第二斜壁板在y方向上受到的力；

s223、竖直壁板中用于连接其他星形胞元的端点在x方向上的位移等于竖直壁板与第二斜壁板的共同端点在x方向上的位移，竖直壁板中用于连接其他星形胞元的端点在y方向上的位移等于竖直壁板拉伸长度加上竖直壁板与第二斜壁板的共同端点在y方向上的位移，由能量法计算出竖直壁板中用于连接其他星形胞元的端点在水平方向上的位移δex和在竖直方向上的位移δey，表示为：

s224、四分之一胞元在x方向上的应变εx″、在y方向上的εy″表示为：

s225、四分之一胞元等效泊松比vyx和等效杨氏模量ey表示为：

更进一步的，随着壁板夹角θ1、θ2增大，星形蜂窝结构在y方向单轴受拉力情况下的杨氏模量ey减小，在y方向单轴受拉力情况下的泊松比νyx增大；

随着水平壁板长度l1增大，杨氏模量ey减小，泊松比νyx增大；

随着第一斜壁板长度l2增大，杨氏模量ey增大，泊松比νyx增大；

随着第二斜壁板长度l3增大，杨氏模量ey减小，泊松比νyx先减小后增大；

随着竖直壁板长度l4增大，杨氏模量ey增大，泊松比νyx减小。

更进一步的，对于星形胞元，当其在与x方向成角的方向上单向受拉时，胞元等效泊松比和等效杨氏模量的求解过程如下：

s311、星形胞元在与水平方向成角方向上的有效横截面积总拉伸力表示为：

式中，为星形胞元在与水平方向成角方向上的应力；

s312、星形胞元总拉伸力在x方向上的投影在y方向上的投影表示为：

s313、星形胞元在x方向上的应变在y方向上的应变在与水平方向成角方向上的应变在与水平方向成角方向上的应变表示为：

s314、星形胞元等效泊松比和等效杨氏模量表示为：

本发明的第二目的通过下述技术方案实现：一种计算机可读存储介质，存储有程序，所述程序被处理器执行时，实现本发明第一目的所述的星形蜂窝结构本构关系的建立方法。

本发明的第三目的通过下述技术方案实现：一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，所述处理器执行存储器存储的程序时，实现本发明第一目的所述的星形蜂窝结构本构关系的建立方法。

本发明的第四目的通过下述技术方案实现：一种星形蜂窝结构本构关系的建立装置，包括：胞元获取模块、单轴受力分析模块和倾斜受力分析模块，其中，

胞元获取模块，用于从由一系列星形胞元组成的星形蜂窝结构中切出一个星形胞元；

单轴受力分析模块，用于在沿着水平方向或竖直方向单向拉伸星形胞元情况下，考虑到星形胞元具有对称性，对于星形胞元的研究，直接取附加边界约束的四分之一胞元进行力学分析；考虑四分之一胞元受力特性以及边界约束，将四分之一胞元的力学分析过程进一步简化为斜壁板的受力变形；然后根据斜壁板的受力变形、星形胞元结构的几何特征及材料杨氏模量，运用能量法求解出胞元结构的等效杨氏模量以及等效泊松比；

倾斜受力分析模块，用于在沿着与水平方向成角方向单向拉伸星形胞元情况下，将星形胞元的力学分析过程简化为星形胞元在水平方向受力和竖直方向受力的叠加，通过单轴受力分析模块求解胞元在x、y方向上的受力变形，即可进一步计算出胞元结构的等效杨氏模量以及等效泊松比。

本发明相对于现有技术具有如下的优点及效果：

(1)本发明设计了星形蜂窝结构本构关系的建立方法，将星形蜂窝结构视为一种具有周期特性的各向异性材料，在小变形情况下，对于沿x方向或y方向单向拉伸胞元结构的情况，可根据胞元受力特性以及边界约束，将星形蜂窝结构的分析转化为分析斜壁板的受力变形；对于沿与x方向成角的方向单向拉伸胞元结构的情况，可将星形蜂窝结构的分析简化为x方向单向拉伸和y方向单向拉伸的叠加，然后针对结构的实际变形、几何特征及材料杨氏模量，运用能量法可得出一系列方程组并对其进行解耦，即可求出结构的位移、杨氏模量和泊松比的等效公式。可见，本发明可以针对结构的实际变形、几何特征及材料杨氏模量，准确方便地建立出反映同一类负泊松比星形节点周期性蜂窝结构的本构关系。

(2)本发明所构建的星形蜂窝结构本构关系考虑到了多种多样的结构参数，因此所等效的代表胞元具有良好的代表性，能够用于模拟星形蜂窝结构进行研究分析。相比于数值方法和试验方法，本发明方法具有更好的计算精度和适用范围，且更便于使用。

附图说明

图1为负泊松比星形节点周期性蜂窝结构的示意图。

图2为胞元基本计算参数的示意图。

图3为星形蜂窝结构在x方向上承受拉荷载时的整体受力示意图。

图4为胞元在x方向上承受拉荷载时的整体受力示意图。

图5(a)为四分之一胞元在x方向上承受拉荷载时的示意图。

图5(b)为只包含斜壁板的结构在x方向上承受拉荷载时的示意图。

图6为星形蜂窝结构在y方向上承受拉荷载时的整体受力示意图。

图7为胞元在y方向上承受拉荷载时的整体受力示意图。

图8(a)为四分之一胞元在y方向上承受拉荷载时的示意图。

图8(b)为只包含斜壁板的结构在y方向上承受拉荷载时的示意图。

图9为星形蜂窝结构在与x方向成角方向单向拉伸的整体受力示意图。

图10为理论模型与有限元模型的分析结果对比示意图。

图11为星形蜂窝结构参数θ1和θ2对本构关系的影响示意图。

图12为星形蜂窝结构参数l1对本构关系的影响示意图。

图13为星形蜂窝结构参数l2对本构关系的影响示意图。

图14为星形蜂窝结构参数l3对本构关系的影响示意图。

图15为星形蜂窝结构参数l4对本构关系的影响示意图。

具体实施方式

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

实施例1

本实施例公开了一种星形蜂窝结构本构关系的建立方法，包括如下步骤：

s1、从星形蜂窝结构中切出一个星形胞元，当沿着水平方向或竖直方向单向拉伸星形胞元，考虑到星形胞元是轴对称的结构，其具有对称性，因此，对于星形胞元的研究，可直接取附加边界约束的四分之一胞元进行力学分析。

星形蜂窝结构可参见图1、图3、图6和图9，其由一系列八边的星形胞元组成，切出来的星形胞元可参见图2、图4和图7，四分之一胞元可参见图5(a)和图8(a)。

如图5(a)和图8(a)所示，四分之一胞元具有相连的两个斜壁板bc、cd，一个自斜壁板bc水平延伸以连接其他星形胞元的水平壁板ab，以及一个自斜壁板cd竖直延伸以连接其他星形胞元的竖直壁板de。

在图5(a)中，当四分之一胞元水平受拉时，水平壁板ab约束竖直位移，竖直壁板de约束水平方向位移，e点约束竖直方向位移。在图8(a)中，当四分之一胞元竖直受拉时，竖直壁板de约束水平方向位移，水平壁板ab约束竖直方向位移，a点约束水平方向位移。

这里，定义x方向代表水平方向(参见图1中的x轴)，y方向代表竖直方向(参见图1中的y轴)，z方向代表与x方向和y方向相垂直的方向，并将连接水平壁板的斜壁板bc定义为第一斜壁板，将连接竖直壁板的斜壁板cd定义为第二斜壁板。

如图2所示，星形胞元结构的几何特征具体为：水平壁板ab的长度为l1，第一斜壁板bc的长度为l2，第二斜壁板cd的长度为l3，竖直壁板de的长度为l4，水平壁板ab与第一斜壁板bc的夹角为θ1，第二斜壁板cd与竖直壁板de的夹角为θ2，每个胞元壁板的横截面积均为a，每个胞元壁板的厚度均为b，胞元结构在z方向上的深度为t。另外，设胞元结构的材料杨氏模量为es，胞元结构的极惯性矩为i。

s2、当沿着水平方向或竖直方向单向拉伸星形胞元时，可简化为对四分之一胞元进行力学分析，考虑四分之一胞元受力特性以及边界约束，对四分之一胞元的力学分析过程可进一步简化为斜壁板的受力变形，可参见图3～图5(b)和图6～图8(b)。其中，为方便描述，令整个胞元水平受力2fx或竖直受力2fy，以便在取四分之一胞元进行分析时，受力可以是fx或fy。

根据斜壁板的受力变形、星形胞元结构的几何特征及材料杨氏模量，运用能量法求解出胞元结构的等效杨氏模量以及等效泊松比：

先求出斜壁板的弯矩方程、轴力方程，以及弯矩、轴力对斜壁板受力的偏导数；

然后基于弯矩方程、轴力方程及其对拉力的偏导数，运用能量法求得斜壁板在水平方向和竖直方向的位移，进而求得四分之一胞元在水平方向和竖直方向的位移，根据四分之一胞元在水平方向的位移计算水平方向的应变，根据四分之一胞元在竖直方向的位移计算竖直方向的应变；

将水平和竖直方向的应变代入杨氏模型计算公式和泊松比计算公式，即可得出四分之一胞元的等效杨氏模量以及等效泊松比，星形胞元的等效杨氏模量以及等效泊松比与四分之一胞元的等效杨氏模量以及等效泊松比相等。

具体来说，1)当星形胞元在水平方向上单向受拉时，星形胞元的等效泊松比vxy和等效杨氏模量ex求解过程如下：

s211、四分之一胞元在x方向上的有效横截面积ax、总拉伸力fx表示为：

s212、和弯曲引起的挠曲变形相比，可以忽略由于剪切胞元引起的挠曲变形，因此，第一斜壁板bc的轴力方程和弯矩方程及其对各力的偏导数表示为：

其中，fbc为第一斜壁板的轴力；mbc(x1)为第一斜壁板上任一点的弯矩；如图5(b)，fb为第一斜壁板与水平壁板的共同端点即b点上的虚加力(即人为定义的虚加的力)，fb＝0；x1为第一斜壁板bc上任意一点至第一斜壁板与水平壁板共同端点即b点的长度，例如b点对应长度x1＝0，c点对应长度x1＝l2；fx表示第一斜壁板在x方向上受到的力；

第二斜壁板cd的轴力方程和弯矩方程及其对各力的偏导数表示为：

其中，fcd为第二斜壁板的轴力；mcd(x2)为第二斜壁板上任一点的弯矩；x2为第二斜壁板cd上任意一点至第二斜壁板与第一斜壁板共同端点即c点的长度，例如c点对应长度x2＝0，d点对应长度x2＝l3；

s213、图5(a)中水平壁板中用于连接其他星形胞元的端点即a点在y方向上的位移等于图5(b)中b点在y方向上的位移，a点在x方向上的位移等于水平壁板ab拉伸长度加上图5(b)中b点在x方向上的位移，由能量法计算出a点在水平方向上的位移δax和在竖直方向上的位移δay，即计算出四分之一胞元在水平方向和竖直方向的位移，表示为：

其中，es为胞元结构的材料杨氏模量；

s214、由于x方向上的应变为x方向上的位移和x方向上的投影长度的比值，y方向上的应变为y方向上的位移和y方向上的投影长度的比值，因此计算出四分之一胞元在x方向上的应变εx′、在y方向上的应变εy′，表示为：

s215、四分之一胞元的等效泊松比vxy和等效杨氏模量ex表示为：

同样，2)对于星形胞元胞元，当其在竖直方向上单向受拉时，等效泊松比vyx和等效杨氏模量ey的求解过程如下：

s221、四分之一胞元在y方向上的有效横截面积ay、总拉伸力fy表示为：

式中，σy为四分之一胞元在竖直方向上的应力；

s222、第一斜壁板bc，第二斜壁板cd的轴力方程和弯矩方程及其对各力的偏导数表示为：

其中，f′bc为第一斜壁板的轴力；f′cd为第二斜壁板的轴力；m′cd(x4)为第二斜壁板上任一点的弯矩；m′bc(x3)为第一斜壁板上任一点的弯矩；x3为第一斜壁板bc上任意一点至第一斜壁板与第二斜壁板共同端点即c点的长度，例如b点对应长度x3＝l2，c点对应长度x3＝0；x4为第二斜壁板cd上任意一点至第二斜壁板与垂直壁板共同端点即d点的长度，例如c点对应长度x4＝l3，d点对应长度x4＝0；如图8(b)，fd为第二斜壁板与竖直壁板的共同端点即d点上的虚加力，fd＝0；fy表示第二斜壁板在y方向上受到的力；

s223、图8(a)中竖直壁板中用于连接其他星形胞元的端点即e点在x方向上的位移等于图8(b)中d点在x方向上的位移，e点在y方向上的位移等于竖直壁板de拉伸长度加上图8(b)中d点在y方向上的位移，由能量法计算出e点在水平方向上的位移δex和在竖直方向上的位移δey，即计算出四分之一胞元在水平方向和竖直方向的位移，表示为：

s224、四分之一胞元在x方向上的应变εx″、在y方向上的εy″表示为：

s225、四分之一胞元等效泊松比vyx和等效杨氏模量ey表示为：

s3、如图9，当沿着与水平方向成角方向单向拉伸星形胞元时，可将星形胞元的力学分析过程简化为星形胞元在水平方向受力和竖直方向受力的叠加，通过求解在x、y方向上的受力变形，即可进一步计算出胞元结构的等效杨氏模量以及等效泊松比

等效杨氏模量以及等效泊松比的求解过程如下：

s311、星形胞元在与水平方向成角方向上的有效横截面积总拉伸力表示为：

式中，为星形胞元在与水平方向成角方向上的应力；

s312、星形胞元总拉伸力在x方向上的投影在y方向上的投影表示为：

s313、星形胞元在x方向上的应变在y方向上的应变在与水平方向成角方向上的应变在与水平方向成角方向上的应变表示为：

s314、星形胞元等效泊松比和等效杨氏模量表示为：

为验证上述方法构建的本构关系理论模型，本实施还构建了星形蜂窝结构的有限元模型来进行分析比较：

图10为本实施例理论模型与有限元模型的对比分析结果，图10中的圆点和三角形分别代表有限元模型分析得到的a点在x方向受到拉力时所发生的x方向和y方向上的位移，实线和虚线分别代表本实施例方法构建的本构关系理论模型分析得到的a点在x方向受到拉力时所发生的x方向和y方向上的位移。从图10可以看出，圆点刚好在实线上，三角形刚好在虚线上，可见，本实施例理论模型与相对精细的有限元模型结果吻合良好。

图11为本实施例中结构参数θ1和θ1对星形蜂窝结构的影响示意图。其中，结构及材料参数设为：θ1分别为45°，50°，55°，60°，65°，70°，75°，θ2分别为45°，50°，55°，60°，65°，70°，75°，l1＝10mm，l2＝10mm，l3＝10mm，l4＝10mm，t＝1mm，b＝1mm，es＝200000mpa。从结果可以看出：1、随着壁板夹角θ1、θ2增大，星形蜂窝结构在x方向的杨氏模量ex和在y方向的杨氏模量ey减小；2、随着壁板夹角θ1、θ2增大，x方向泊松比νxy和y方向泊松比νyx增大。

图12为本实施例中结构参数l1对星形蜂窝结构的影响示意图。其中，结构及材料参数为：θ1＝60°，θ2＝60°，l1分别为7，8，9，10，11，12，13mm，l2＝10mm，l3＝10mm，l4＝10mm，t＝1mm，b＝1mm，es＝200000mpa。从结果可以看出：1、随着壁板长度l1增大，x方向杨氏模量ex增大，y方向杨氏模量ey减小；2、随着壁板长度l1增大，x方向泊松比νxy减小，y方向泊松比νyx增大。

图13为本实施例中结构参数l2对星形蜂窝结构的影响示意图。其中，结构及材料参数为：θ1＝60°，θ1＝60°，θ2＝60°，l1＝10mm，l2分别为7，8，9，10，11，12，13mm，l3＝10mm，l4＝10mm，t＝1mm，b＝1mm，es＝200000mpa。从结果可以看出：1、随着壁板长度l2增大，x方向杨氏模量ex减小，y方向杨氏模量ey增大；2、随着壁板长度l2增大，x方向泊松比νxy先减小后增大(l2为9mm左右时泊松比最小)，y方向泊松比νyx增大。

图14为本实施例中结构参数l3对星形蜂窝结构的影响示意图。其中，结构及材料参数为：θ1＝60°，θ2＝60°，l1＝10mm，l2＝10mm，l3分别为7，8，9，10，11，12，13mm，l4＝10mm，t＝1mm，b＝1mm，es＝200000mpa。从结果可以看出：1、随着壁板长度l3增大，x方向杨氏模量ex增大，y方向杨氏模量ey减小；2、随着壁板长度l3增大，x方向泊松比νxy减小，y方向泊松比νyx先减小后增大(l3为9mm左右时泊松比最小)。

图15为本实施例中结构参数l4对星形蜂窝结构的影响示意图。其中，结构及材料参数为：θ1＝60°θ1＝60°，θ2＝60°，l1＝10mm，l2＝10mm，l3＝10mm，l4分别为7，8，9，10，11，12，13mm，t＝1mm，b＝1mm，es＝200000mpa。从结果可以看出：1、随着壁板长度l4增大，x方向杨氏模量ex减小，y方向杨氏模量ey增大；2、随着壁板长度l4增大，x方向泊松比νxy增大，y方向泊松比νyx减小。

从图11至图15可得，本实施例方法所构建的星形蜂窝结构本构关系考虑到了多种多样的结构参数，所等效的代表胞元具有良好的代表性，能够用于模拟星形蜂窝结构进行研究分析。

实施例2

本实施例公开了一种计算机可读存储介质，存储有程序，所述程序被处理器执行时，实现实施例1所述的星形蜂窝结构本构关系的建立方法，具体如下：

s1、从由一系列星形胞元组成的星形蜂窝结构中切出一个星形胞元；

s2、当沿着水平方向或竖直方向单向拉伸星形胞元时，考虑到星形胞元具有对称性，对于星形胞元的研究，直接取附加边界约束的四分之一胞元进行力学分析；

考虑四分之一胞元受力特性以及边界约束，将四分之一胞元的力学分析过程进一步简化为斜壁板的受力变形；然后根据斜壁板的受力变形、星形胞元结构的几何特征及材料杨氏模量，运用能量法求解出胞元结构的等效杨氏模量以及等效泊松比；

s3、当沿着与水平方向成角方向单向拉伸星形胞元时，将星形胞元的力学分析过程简化为星形胞元在水平方向受力和竖直方向受力的叠加，通过求解在x、y方向上的受力变形，即可进一步计算出胞元结构的等效杨氏模量以及等效泊松比

本实施例中的计算机可读存储介质可以是磁盘、光盘、计算机存储器、只读存储器(rom，read-onlymemory)、随机存取存储器(ram，randomaccessmemory)、u盘、移动硬盘等介质。

实施例3

本实施例公开了一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，所述处理器执行存储器存储的程序时，实现实施例1所述的星形蜂窝结构本构关系的建立方法，具体如下：

s1、从由一系列星形胞元组成的星形蜂窝结构中切出一个星形胞元；

s2、当沿着水平方向或竖直方向单向拉伸星形胞元时，考虑到星形胞元具有对称性，对于星形胞元的研究，直接取附加边界约束的四分之一胞元进行力学分析；

考虑四分之一胞元受力特性以及边界约束，将四分之一胞元的力学分析过程进一步简化为斜壁板的受力变形；然后根据斜壁板的受力变形、星形胞元结构的几何特征及材料杨氏模量，运用能量法求解出胞元结构的等效杨氏模量以及等效泊松比；

s3、当沿着与水平方向成角方向单向拉伸星形胞元时，将星形胞元的力学分析过程简化为星形胞元在水平方向受力和竖直方向受力的叠加，通过求解在x、y方向上的受力变形，即可进一步计算出胞元结构的等效杨氏模量以及等效泊松比

本实施例中所述的计算设备可以是台式电脑、笔记本电脑、智能手机、pda手持终端、平板电脑或其他具有处理器功能的终端设备。

实施例4

本实施例公开了一种星形蜂窝结构本构关系的建立装置，可实现实施例1所述的星形蜂窝结构本构关系的建立方法，所述建立装置包括：胞元获取模块、单轴受力分析模块和倾斜受力分析模块。

具体来说，胞元获取模块，用于从由一系列星形胞元组成的星形蜂窝结构中切出一个星形胞元；

单轴受力分析模块，用于在沿着水平方向或竖直方向单向拉伸星形胞元情况下，考虑到星形胞元具有对称性，对于星形胞元的研究，直接取附加边界约束的四分之一胞元进行力学分析；考虑四分之一胞元受力特性以及边界约束，将四分之一胞元的力学分析过程进一步简化为斜壁板的受力变形；然后根据斜壁板的受力变形、星形胞元结构的几何特征及材料杨氏模量，运用能量法求解出胞元结构的等效杨氏模量以及等效泊松比；

倾斜受力分析模块，用于在沿着与水平方向成角方向单向拉伸星形胞元情况下，将星形胞元的力学分析过程简化为星形胞元在水平方向受力和竖直方向受力的叠加，通过单轴受力分析模块求解胞元在x、y方向上的受力变形，即可进一步计算出胞元结构的等效杨氏模量以及等效泊松比。

在此需要说明的是，本实施例的装置仅以上述各功能模块的划分进行举例说明，在实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。