陈帅， 王兵， 朱绍伟

(1.哈尔滨工业大学 特种环境复合材料技术国家级重点实验室，黑龙江 哈尔滨 150001； 2.重庆大学 航空航天学院， 重庆 400044)

泊松比(Poisson′s ratio，PR)目前广泛用于表征材料的力学性能[1]。自然界中材料的泊松比大多介于0.2～0.5，泊松比为负值或为零的材料在自然界中较少，需要通过人为设计得到。

负泊松比(negative Poisson′s ratio, NPR)材料在拉伸时呈横向扩张，在压缩时呈横向收缩，故又被称为拉胀材料[2]。与具有正泊松比的传统材料相比，负泊松比材料具有优异的力学性能，如抗剪切性[3-4]、抗压痕性[5]、能量吸收和断裂韧性[6-7]等可以显著提高，这使其在各个领域具有广阔的应用前景[8-9]，特别是在航空航天、生物医疗、国防以及智能系统等领域。并已代表性地用于飞机涡轮发动机的叶片中。

自然界中的负泊松比材料确实存在，如奶牛乳头部分皮肤和猫的皮肤[10]。Ting等[11]的工作表明，对于各向异性材料，只要应变能密度函数是正定的，泊松比可以有任意正负值。因此，泊松比的负值不仅是假设的，而且可以在人造材料中实现。Lakes[3]首次人为造出具有负泊松比的聚氨酯泡沫。在后续研究中根据材料内部胞元的变形机制，负泊松比材料可分为内凹型[8]、手性型[12]、旋转型[13]和穿孔型[14]。

Gibson等[15]提出了内凹型蜂窝结构。Master等[16]发现了二维内凹结构的负泊松比行为。当在任一方向施加载荷时，对角肋的移动方式会导致另一个方向的负泊松比效应。另一种由于内凹性而表现出负泊松比效应的多胞结构是双箭头[17-18]和星形[19]。这些结构中的负泊松比机制类似于内凹结构。在拉伸过程中，双箭头和星形单元打开，导致拉胀效应。

Prall等[20]提出手性蜂窝结构的概念。手性是指镜像上不能重叠，手性结构由等半径的圆形元素通过与节点相切的直韧带(壁)连接而成。这种配置通过节点的旋转和相关韧带的弯曲提供了平面内的负泊松比效应，以及允许在组成材料的弹性范围内产生高应变的变形机制。这些结构的泊松比约为-1[12]。

Grima等[13]提出了带有旋转单元的负泊松比结构。在此类结构中，负泊松比行为是由通过铰链相互连接的刚性多边形的旋转获得的，如旋转正方形[13]、旋转矩形[21]、旋转平行四边形和菱形[22]、旋转三角形[23-25]和旋转四面体[26]。基于旋转刚性单元模型的灵感，负泊松比行为也在具有菱形或星形开孔[14]的片材中得到验证。

零泊松比(zero Poisson′s ratio, ZPR)材料在受单向拉压时，不会出现横向收缩和扩张。自然界中软木等材料所呈现的泊松比被观察到等于或接近于零[27-28]。零泊松比材料因其不同寻常的特性和在医学[29]、组织工程[30]和航空[31]等不同领域的潜在应用而备受关注。在航空领域，变形机翼中使用的柔性蒙皮需要提供足够的协同变形以承受空气载荷。此外，根据光滑度的要求，需要对柔性蒙皮的不变形方向进行限制。在这些方案中，泊松比为零的材料更为可取[32]。目前针对零泊松比材料的设计主要有2种思路：1)通过材料胞元的设计实现零泊松比效应；2)混合正、负泊松比胞元获得重组后的具有零泊松比效应的新型胞元。

程文杰等[32]提出一种新型十字型混合零泊松比蜂窝结构并进行了理论验证。何玉龙等[33]通过折纸结构设计出新型零泊松比结构，并对其进行了几何分析及其耗能性能研究。Huang等[34]提出了一种由六边形部分和薄板部分组成的新颖的平面内零泊松比蜂窝结构。六边形结构用于承受面外压缩并产生面内柔性，而薄板用于产生面外柔性。这2个部分对有效力学性能有不同的贡献，需要分别设计平面内和平面外的性能。

针对正、负泊松比胞元混合构型的设计也有广泛的研究。Olympio[35]提出的构型可以实现整体微元的零泊松比效应。当周期微元受到一个方向的载荷时，微元中的正、负泊松比胞元分别收缩(膨胀)和膨胀(收缩)相同的横向变形量，整个周期微元的等效横向变形为零，即表现出零泊松比效应。Gong等[36]设计了一种零泊松比蜂窝，可以实现沿2个正交方向的变形，并通过抑制非变形方向的泊松变形，避免了变形方向的有效刚度大幅增加。

本文基于曲杆/板组成的内凹式十边形多孔结构提出了一种正、负和零泊松比相互转换的新策略。即通过改变杆的曲率，结构的泊松比可实现为正、负和零。

1 结构的设计与制备

1.1 几何构型设计

本文提出的内凹式十边形蜂窝结构如图1(a)所示，由4×10个胞元组成。其内凹十边形胞元如图1(b)所示。胞元的每条边均由如图1(c)所示的正弦梁组成。当竖向力作用于蜂窝结构时，曲梁会偏转和弯曲，导致整体结构的大变形行为。通过改变曲梁的曲率参数可以在整个蜂窝结构中实现不同的泊松比效应。余弦梁的几何形状可以表示为:

(1)

式中：t为余弦梁的函数形式；L为余弦梁的长度；k为余弦梁的曲率。2个胞元之间均由横向直梁连接。横向直梁的存在可以联动整个结构的横向变形。

图1 内凹式十边形蜂窝结构的几何构型设计Fig.1 Geometric configuration design of re-entrant decagonal honeycomb structure

1.2 样品制备

制备样品的材料体系选用热塑性聚氨酯。热塑性聚氨酯是增材制造技术中最常用的聚合物之一，在世界上应用广泛。具有来源广、韧性好、强度高、成本低、耐高温等优点。且热塑性聚氨酯是一种成熟的环保材料，具有优良的耐磨性、耐老化性和安全性等优点。采用熔融沉积(FDM)3-D打印技术(HORI Z500，北京汇天威科技有限公司)对样品进行制备。图2(a)为样品的3-D打印示意图。制备的样品实物如图2(b)所示。这种增材制造技术有许多优点：可以制造不同材料和尺寸的复杂实体或腔体结构；结构的打印精度极高(可达0.05 mm)，制造成本较低。根据ASTM D638-14拉伸实验测试标准，获得了原材料热塑性聚氨酯的拉伸应力-应变曲线，如图3所示，其材料基本力学性能如表1所示。

图2 样品的制备过程示意Fig.2 Schematic diagram of the preparation process of the sample

图3 原材料热塑性聚氨酯的拉伸应力-应变曲线Fig.3 The stress-strain curve of the base material TPU

表1 热塑性聚氨酯的弹性属性Table1 The elastic properties of thermoplastic polyurethane

2 研究方法

2.1 有限元方法

为了研究和探索结构在竖向荷载作用下的变形特征及力学性能，采用商用有限元(FE)计算软件ABAQUS/Explicit求解器进行了数值模拟。首先通过建模软件SOILDWORKS建立了结构的三维CAD模型，并导入到ABAQUS中。在有限元模型中，由于结构的瞬变特性，考虑了几何非线性。采用适用于求解大应变问题的八节点线性减缩积分六面体单元(C3D8R)对蜂窝结构进行网格划分。采用线性减缩积分单元的优点是：当网格发生畸变时，分析精度不会受到很大影响[37-38]。在每个模型中，结构横截面上都有3个以上的单元，以满足网格收敛性要求。表1中列出的热塑性聚氨酯(TPU)材料的力学参数应用到有限元模型中。拉伸试验的有限元模型如图4所示。蜂窝上下表面用卡具固定，并采用体积中心有参考点的离散刚体对拉伸卡具进行模拟。在模拟过程中，卡具与试样之间进行绑定。同时，顶部参考点受竖向位移载荷作用，其他5个自由度均固定。同时，限制底部卡具基准点的所有自由度。为了消除边界效应，在结构中间设置4个参考点作为测量点，用于监测结构的横向变形和竖向变形。

图4 蜂窝结构有限元模型Fig.4 The FE model of the honeycomb structure

2.2 实验方法

为了验证有限元结果的有效性，采用INSTRON 5569试验机对结构沿垂直方向进行了准静态拉伸试验。准静态拉伸试样的实验装置如图5(a)所示，试样上下两端通过特制夹具进行固定，以便和试验机夹头进行连接。实验过程采用位移加载控制，恒定位移速率设置为1 mm/min，避免了动态影响。在试验过程中，监测试样的横向变形和纵向变形。为了准确测量试样的泊松比值，提出一种简易的DIC技术方法。拉伸试验前，在样品中心标记点，如图5(b)所示。在拉伸过程中，使用高分辨率数码相机对拉伸试样进行视频记录，如图5(a)所示。然后将视频裁剪成反映整个拉伸过程的一系列连续图片，利用图像处理技术得到每个标记点中心的运动轨迹，通过轨迹的计算即可得到试样的纵向应变-横向应变图，继而得到结构泊松比值，其整个方法过程如图6所示。实验测试结果也用于验证有限元结果。

图5 实验技术方法Fig.5 Experimental techniques and methods

3 结果分析与讨论

3.1 胞元构型的确定

图7对比了正常凸式六边形胞元、内凹式六边形胞元和内凹式十边形胞元的几何特征。可以看出结构的泊松比效应是由A、B和C3个节点的移动方向共同决定的。对于正常凸式六边形和内凹式六边形胞元，3个节点坐标在x方向的位置满足:xA=xC≠xB。在y向拉伸载荷下，节点B在x方向的坐标xB总是朝向xA和xc坐标方向移动。这种外凸(或内凹)的机制使正常六边形胞元和内凹式六边形胞元呈现相反的泊松比效应，分别表现出正泊松比效应和负泊松比效应。在内凹式十边形胞元中，3个节点的坐标在x方向的位置满足:xA=xC=xB。因为这种胞元同时存在内凹和外凸的机制，如果杆是直的，在y向拉伸载荷下节点B的坐标xB依然和节点A和C的坐标一致，即xA=xC=xB，此种情况下胞元的泊松比表现为零。但是，如果杆具有特殊的曲率，则在y方向拉伸情况下节点B的坐标xB将远离节点A和C的坐标(xA和xc)，即表现出正或负泊松比效应，其具体分析见下节。因此，本文选用具有不同曲率的内凹十边形胞元组成具有不同泊松比效应的蜂窝结构。

图6 简易的DIC技术方法Fig.6 Simple DIC technology method

3.2 曲率的影响

本文考虑了3种杆曲率内凹十边形胞元，如图8所示。在局部坐标系(s,t)中杆的曲率函数被表达为式(1)。为简化分析，把L视为一个常数。3种胞元杆的曲率k分别取为2、0和-2。限定t′(s)|s=0=t′(s)s=L=0，所以对于不同的曲率k，杆的末端总是具有相同的方向。

为了拓展工程应用，曲板也常作为内凹十边形多孔结构的胞元壁，且为说明胞元在拉伸载荷和压缩载荷下表现出同种性质的泊松比效应，本节中3种具有不同曲率k的胞元结构均由曲板组成，且受到压缩载荷作用，其变形情况如图8所示。可以看出，当k>0时,胞元结构呈现正的泊松比；当k=0时，胞元结构呈现零泊松比；当k<0时，胞元结构呈现负的泊松比。很明显，对于不同的k，胞元可以实现正、负和零泊松比的转换。

图7 胞元构型的对比Fig.7 Comparison of cell configurations

图8 曲率的影响Fig.8 The effect of curvature

3.3 实验结果与仿真结果的对比

本节在3.2节单胞研究的基础上进行了蜂窝结构的研究。蜂窝样品的结构参数为：曲梁长度L=10 mm，曲梁面内厚度T=1 mm，所有试样的面外厚度均为b=1.6 mm。图9展示了k=2情况下蜂窝结构在拉伸载荷下的响应分析。将实验得到的蜂窝结构纵向位移(dy)和横向位移(dx)响应与仿真结果进行了对比，如图9(d)所示。蜂窝结构在拉伸载荷下呈现出的变形形态如图9(a)、(b)所示。图9(e)为蜂窝结构在变形过程中表现出的泊松比值。从图9中可以看出，实验结果与有限元预测在泊松比量值和变形形态方面表现出定量一致性。其结果略区别于单胞结构，即当k>0时，蜂窝结构呈现正的泊松比，随着真实应变的增大，结构的泊松比呈现减小趋势，随后出现负的泊松比。结构的泊松比从正到零，再到负，是一条近乎线性的曲线。这是因为在小变形时，结构节点位置大致满足xA=xC=xB，非零泊松比是由杆的曲率效应造成的。变形较大时，杆件的曲率减小,xB偏离xA和xC，即出现位错。随着位错的增加和曲率的减小，内凹(或凸)机制将取代曲率效应成为影响泊松比的主要因素。在小变形条件下，有限元分析结果与实验结果吻合较好。在较大的变形量下，虽然泊松比的变化趋势相似，但有限元计算结果与实验结果有较大偏差。其偏差主要由以下2个原因造成的：1)熔融沉积(FDM)3-D打印技术在结构尺寸过大的情况下，打印质量存在一定缺陷；2)在有限元计算中把材料假设成了理想弹性，而结构在大变形情况下会有部分塑性变形发生，这也是造成结果偏差的原因之一。

图9 k=2情况下蜂窝结构在拉伸载荷下的响应分析Fig.9 Response analysis of honeycomb structure under tensile load whenk=2

图10展示了k=0情况下蜂窝结构在拉伸载荷下的响应分析。图10(d)为实验得到的蜂窝结构纵向位移(dy)和横向位移(dx)响应与仿真的对比结果。蜂窝结构在拉伸载荷下呈现出的变形形态如图10(a)、(b)所示。图10(e)为蜂窝结构在变形过程中表现出的泊松比值。从图10中可以看出，实验结果与有限元预测在泊松比量值和变形形态方面表现出定量一致性。即当k=0时，蜂窝结构随着真实应变的增大始终呈现近乎零的泊松比。这是因为在结构变形时，节点位置坐标大致始终满足xA=xC=xB。

图10 k=0情况下蜂窝结构在拉伸载荷下的响应分析Fig.10 Response analysis of honeycomb structure under tensile load whenk=0

图11展示了k=-2情况下蜂窝结构在拉伸载荷下的响应分析。图11(d)为实验得到的蜂窝结构纵向位移(dy)和横向位移(dx)响应与仿真的对比结果。蜂窝结构在拉伸载荷下呈现出的变形形态如图11(a)(b)所示。图11(e)为蜂窝结构在变形过程中表现出的泊松比值。从图11中可以看出，实验结果与有限元预测在泊松比量值和变形形态方面表现出定量一致性。即当k<0时，结构的泊松比从负到零，再到正。其原因同样是：在小变形时，结构节点位置大致满足xA=xC=xB，非零泊松比是由杆的曲率效应造成的；大变形时，内凹(或凸)机制将取代曲率效应成为影响泊松比的主要因素。

图11 k=-2情况下蜂窝结构在拉伸载荷下的响应分析Fig.11 Response analysis of honeycomb structure under tensile load whenk=-2

4 结论

1)当曲率k>0时，结构呈现正的泊松比，随着真实应变的增大，结构的泊松比呈现减小趋势，随后出现负的泊松比。结构的泊松比从正到零，再到负，是一条近乎线性的曲线。

2)当曲率k=0时，结构在整个变形中始终表现为近乎零的泊松比。

3)当曲率k<0时，结构呈现负的泊松比，随着真实应变的增大，结构的泊松比转变为正。即结构的泊松比从负到零，再到正。

不同曲率下内凹十边形蜂窝结构的泊松比绝对值随着有效真实应变的增大，逐渐减小，而后出现了相反的泊松比。对这种现象进行了分析，并提出了一种曲率效应。结果表明，在小变形条件下，杆件的非零泊松比是由杆件的曲率效应引起的。而在大变形条件下，泊松比的转换是由结构内凹(或凸)的机制取代曲率效应成为主导机制引起的。