四杆张拉整体单元的弯曲变形

罗阿妮，邹云涛，肖诗松，刘贺平

（哈尔滨工程大学 机电工程学院，黑龙江 哈尔滨 150001）

1 引言

张拉整体是由连续受拉索和离散的受压杆构成的空间结构［1］，所有受拉索连续、受压杆分散，且只承受轴向力，因而材料利用率高。在稳定构型中，杆、索构件的自由度均为零，并且破坏掉任意一根杆、索构件，整个结构都会发生坍塌。结构既不是冗余的，也不是过约束的，具有一定的连续性［2］，可以通过调整个别构件来改变整体刚度。不论是自然状态下还是受到冲击振动，不利因素都可以自动分散至整个结构，找到稳定状态。

张拉整体作为一种空间可展结构，具有独特的构型方式，被广泛应用于伸展臂［3-4］、载体运输工具［5］、桥梁、建筑［6-7］、以及多种艺术品［8］中。Schorr［9］等人利用曲杆代替直杆，通过控制模型的整体重心使张拉整体实现二维平面的滚动，张拉整体曲杆的建立方式为弯杆张拉整体弯角节点提供了帮助。Li［10］等人将张拉整体应用于交通吸振，提出一种混合张拉整体支撑柱，该支撑柱可有效减少因车辆碰撞或其他因素产生的不利振动，耗散因振动而产生的能量，将张拉整体优良的缓冲、吸振及自稳定性应用到实际工程中。Liu［11］等人将球形张拉整体部分索构件用弹簧代替，增强了张拉整体的缓冲性能，为张拉整体压缩提供一种方式，但整体刚度变小、稳定性变差。González 等人［12］通过驱动部分构件或全部构件实现模型从初始状态至完全压缩，找出模型在压缩过程中始终保持稳定的一种方式，为研究张拉整体在压缩过程中的节点坐标变化奠定基础。Kan［13］等人分析了多层张拉整体节点的展开速度和展开过程中的振幅和振动频率。他们的分析方法同样可应用于研究张拉整体的压缩过程，如压缩速度对结构的动态响应、压缩过程系统的稳定性等。

张拉整体具有质量轻、折展比大、缓冲抗振能力强、材料利用率高及使用条件不严格等特点，因此成为了空间结构领域的研究热点。但是张拉整体自身几何尺寸的限制，在某些情况下载物与压杆间会产生干涉，如果偏转角选择不当，自身杆构件也会发生干涉。这就对如何保持原有的结构特性，增大内部有效利用空间，提高缓冲性能，实现一定的自身折展［14-15］提出了更高的要求。本文对张拉整体轴向折叠方式及折展特性进行了理论分析及实验研究，提出了两种张拉整体轴向折叠的方法，不仅避免了杆构件间的干涉，而且提高了有效利用空间。

2 张拉整体基本单元的构建

2.1 四杆张拉整体基本单元

图1 为四杆张拉整体基本单元（彩图见期刊电子版）。构成单元的基本要素是：下底面圆上均匀分布4 个节点n1，n2，n3，n4，任意两个相邻节点与下底面圆圆心均成90°，4 个节点构成正方形；距离下水平面为h的上方以同样的形式均布4 个节点，并且上下水平面间形成一定单元内转角φ。蓝色表示受压杆构件，红色表示受拉索构件，构件均匀分布，并具有一定的对称性。每个节点都有4 个构件相连接，并且受力是平衡的。

图1 四杆张拉整体基本单元结构示意图Fig.1 Schemati diagram of 4-bar basic tensegrity unit

2.2 节点坐标

该系统的基本工作原理为：张拉整体节点坐标是构件张拉整体结构的基础，得到节点坐标便可知任意杆、索的空间位置和长度，节点坐标也是判断结构是否稳定的主要依据。四杆张拉整体共有8 个节点，下水平面节点坐标用Nd表示，上水平面节点用Nu表示，R为上下水平面外接球半径，节点矩阵表示如下：

式中i∈(0，1，2，3)。当四杆张拉整体单元稳定时，偏转角φ为：

2.3 连接矩阵

连接矩阵是将节点与节点间相互连接的介质，它能将节点坐标与矢量联系在一起。构件的两端点分别为起始点和终点，对于4 杆张拉整体，每节点有4 个构件与之相连，所组成的连接矩阵为连接杆构件与杆构件间的矩阵成为杆连接矩阵连接索构件与索构件间的矩阵成为索构件矩阵，它们的关系如下：

连接矩阵中元素仅由0，-1，1 构成，表1 为连接矩阵各元素的含义。

表1 连接矩阵中各元素的含义Tab.1 Element meanings of connection matrix

则上图中的杆连接矩阵为：

索连接矩阵为：

2.4 矢量矩阵

张拉整体构件矢量是将节点向向量进行转换的一种方式，通过终点坐标与起始点坐标做差得到构件矢量。构件矢量的向量长度即为杆构件长度，起始点向量已知，便可以表示出杆构件的空间位置。图2 为四杆张拉整体构件矢量示意图。

图2 张拉整体构件矢量示意图Fig.2 Schematic diagram of component vectors of tensegrity

杆构件的矢量矩阵为：

索构件的矢量矩阵分为上、下水平索以及分布在上下水平索中间的斜索3 部分，其表达式为：

3 四杆张拉整体轴向折叠

3.1 轴向折叠的意义

张拉整体是一种自平衡结构，任意构件受力时，该杆件的受力变形会发散至整体结构，使整体结构分散均匀受力，具有良好的稳定性。受拉构件分布在周围，受压构件分布在内部，有效地减小结构内部的损伤程度。由于只承受轴向力，并且构件材料利用率高，张拉整体的构件数目少、质量轻。

四杆张拉整体在构型不同时其力学性能是不同的，它共有3 种稳定构型。其中，第二种稳定构型杆构件间是相互干涉的，并且杆数越多相互干渉越严重。为了使张拉整体结构使用前占用空间小，需对结构进行轴向压缩。对张拉整体轴向加力时，索构件进一步被拉伸，杆构件进一步被压缩，但杆构件的轴向变化率较小，超过某根杆构件可承受强度，还可能造成整体模型的破坏。

3.2 轴向压缩方式

张拉整体轴向压缩过程中，杆构件进一步被压缩，索构件进一步被拉伸，因此轴向压缩主要有两种方案：

方案一：将受拉索构件更换为具有较大弹性的构件，例如在索构件中添加弹簧，将张拉整体斜索用弹簧代替［16］，如图3 所示，也可以将水平索用弹簧代替。

图3 张拉整体斜索等效替换Fig.3 Vertical string equivalent substitution of tensegrity

方案二：将杆构件复杂化，将单根直杆变成带有套管的组合构件，套管中间放置弹簧［16］结构，两端分别与杆构件相连，实现杆构件的压缩，如图4 所示。

图4 伸缩杆Fig.4 Telescopic bar

这两种方法均添加了额外的弹簧结构，当模型下降至全模型高度的60%时，杆构件间发生干渉，轴向无法继续压缩。

4 弯杆张拉整体设计

为了解决上述问题，本文将张拉整体直杆变成弯杆。张拉整体直杆变为弯杆不仅可以改变张拉整体承受轴向力时杆构件的受力方式，而且解决了直杆间干涉的问题。当施加轴向力时，整体构件会随着外力的增大发生明显的形变，这就使张拉整体在外力的作用下储存了一部分能量。很多领域将张拉整体作为缓冲装置，如NASA 航天局外太空探测器的缓冲投放装置、深海领域探测器的缓冲投放装置、陆地补给箱的缓冲投放装置等。张拉整体直杆变弯可实现轴向压缩，减小占用空间，并且提升缓冲性能，为张拉整体的伸展性及折叠性提供了一种新的思路。

4.1 模型节点的构建方式

弯杆张拉整体是在传统直杆平行面张拉整体结构的基础上，添加弯曲角度θ，将内部的直杆变弯。张拉整体弯杆节点在整体高度H的一半处，这样建立起来的弯杆张拉整体仍然为完美的对称结构，便于力学研究。

弯角处节点坐标的求解方式有二种。首先，为保证结构对称性，弯角节点选在α=H/2 平面处，且过杆构件端点nj，nj+p的中垂面β上。在这种条件约束下，弯角节点坐标被限制在一条直线上。图5 为弯角节点坐标与两平面间的关系，即直线Am为两平面相交所得。将直线njnj+p的中点设为A，在交线上除A点外任取一点m与两端点相连会形成∠njmnj+p。预设弯杆张拉整体弯曲角度，可以找到确定的两个点，令m'为m点关于点A的对称点。为了增大内部空间，m点与m'点中有一点是向外弯曲的，另一点则是向内部弯曲（舍弃）。这样每给出一个确定的弯曲角度就会的到一个确定的弯点，这种思路适用于任意杆数P的弯杆张拉整体。

图5 弯曲节点的坐标空间位置Fig.5 Coordinate space position of bending nodes

其次，为了方便实现弯杆张拉整体模型的绘制，通过向量合成方式得到图像。欲得到m点的准确坐标，需要得到nj A，∠mnj A。通过水平面的法向量，即向量［0，0，1］与njnj+p方向向量做外积得到Am方向向量，方向向量取模后便得到单位方向向量。该向量长度为：

A点坐标为(nj+nj+p)/2，从A点出发做向量Am，向量终点m即为弯杆张拉整体弯角处坐标。通过绘图命令将中点m分别与原直杆两端点相连，得到弯杆张拉整体结构示意图。该方法同样适用于任何构型，任何杆数的张拉整体。

弯杆角度为150°，H=270 mm，Ru=Rd=200 mm，P=4 的3 种结构示意图如图6 所示，从左到右单元内转角分别为135°，180°，225°。

图6 四杆弯曲张拉整体Fig.6 Four-bar bending tensegrity

4.2 基本参数之间的关系

张拉整体构型具有很强的规律性，已知一些参数，如上底面半径Ru，下底面半径Rd，高度H，偏转角φ，便可以得到p-1 种确定的构型。当构型确定后，张拉整体结构也唯一确定。然后，进一步给出弯曲方向及弯曲角度λ，那么弯角处的节点坐标也唯一。这里从节点坐标出发，推导出弯杆长度、弯曲角度和斜索长度之间的关系。

在得到弯曲前杆构件长度后，根据斜杆弯曲角度得到每一根弯杆的长度，如图7 所示。利用相同的方法对任意p杆，任意p-1 种构型进行推导，得到杆索长度的一般表达式。第j种构型节点1 斜索与节点p-（j-1）相连，节点1 杆构件与节点p+1 相连，单元内转角φ=π/2+jπ/p，弯杆角度为λ。节点1、节点p+1、节点p-（j-1）的坐标如下：

图7 弯杆长度Fig.7 Bending length

其中：i=p+1，p−(j−1)。

杆构件长度为：

斜索长度为：

弯杆张拉整体弯曲角度为λ，由三角函数关系得每一根弯杆长为：

4.3 结构分析

在张拉整体模型中节点所连接杆构件的数目称为张拉整体阶数，一阶、二阶张拉整体如图8所示。一阶张拉整体每个节点有且仅有一根杆构件相连，二阶张拉整体每个节点有且仅有两个杆构件相连。张拉整体的复杂程度随着阶数的增加而提高。

图8 不同阶数的张拉整体Fig.8 Tensegrity with different orders

弯杆张拉整体在直杆基础上，添加角度利用弯杆代替直杆，其复杂程度介于一阶与二阶之间。在中间平面弯曲节点处，每个节点连接两根杆构件，而对于未发生变化的上、下水平面处的节点仍连接一根杆构件。在同一张拉整体模型中，既有连接一根杆构件的节点，也有连接两根杆构件的节点，故可将弯杆张拉整体理解为一阶向二阶过渡的一种结构。相比于直杆张拉整体而言，弯杆张拉整体在中间平面增加了弯角节点，因而原有的节点矩阵、连接矩阵和构件矩阵维度等均会发生变化。弯曲节点表达式稍复杂，这里仅用mx，my表示。

对于P杆直杆张拉整体，节点数为2P，节点矩阵维数为3⋅2P；弯杆张拉整体将直杆一分为二，因而节点数变为3P，节点矩阵维数为3⋅3P。当P=4，j=1 时，节点矩阵为：

直杆连接矩阵2P个节点，P根杆，3P根索，直杆连接矩阵维数为：2P⋅4P；弯杆张拉整体索数量不变，杆数由P变为2P，因此弯杆连接矩阵维数为：3P⋅5P。当P=4，j=1 时，连接矩阵为：

构件矩阵是节点矩阵与连接矩阵的乘积。直杆张拉整体共P根杆构件，3P根索构件，共4P个向量，直杆构件矩阵维数为：3⋅4P；弯杆张拉整体共2P根杆，索构件数量不变，弯杆构件矩阵维数为：3⋅5P。当P=4，j=1 时，构件矩阵为：S=，仍然是弯杆节点矩阵与弯杆连接矩阵的乘积。

4.4 压缩过程中的弹性分析

图9 弯杆张拉整体压缩过程Fig.9 Compression process of bending bar tensegrity

对任意p杆，任意p−1 种构型进行推导，得到结构进行轴向压缩过程中杆、索长度及压缩量的表达式。任何构型节点1 均与节点p+1 以杆构件相连，但斜索连接方式与构型有关，第j种构型节点1 与2p−j+1 相连，同一根杆构件，下底面节点与上底面在下底面投影和下底面圆心间形成的角度称为单元内转角（下文均称单元内转角），其大小为：φ=π/2+jπ/p。节点坐标如下：

弯杆张拉整体轴向压缩过程中杆构件连接两点间距离相互靠近，而绳索伸长量远远小于杆构压缩量，在分析过程中将绳索受拉刚度看作无限大，即绳索不可伸长。在轴向压缩过程中，上下两端面有相对转动趋势。因此，在结构从高h压缩至上下水平面重合时，应综合分析高度和变化过程中偏转角的影响，如图10 所示。

图10 压缩过程参数变量分析Fig.10 Parameter variable analysis of compression pro⁃cess of bending bar tensegrity

图10 中，上底面节点C在下底面投影为C1，下底面节点A、投影点C1及下底面圆心O形成的夹角为φ1；当高度被完全压缩C点落在下底面C2，投影点C1，C2及底面圆心角O形成的 夹角为φ2。

压缩过程中，上底面相对于下底面做逆时针方向转动。上、下底面半径一定，高度越高，压缩时的相对转角越大，但高度过大时，转动至某一位置，继续压缩，斜索便会松弛。在压缩过程中，C点运动至线段MM'中间某一位置，假设运动至图中C3点，此时斜索AC为直角△AC3M斜边，两直角边分别为AM'和C3M'，此时为斜索最后绷紧时刻，继续压缩两直角边长度均会缩短，因而超过这一临界位置斜索会松弛，张拉整体结构松垮，失去轴向力而快速轴向压缩，结构不稳定。当发生上述不稳定情况时，杆长的压缩量超过整根原始杆长的50%，在实际模型搭建过程中是不会出现的，实际压缩过程中某根索发生松弛的主要原因是搭建模型精度不够高。只考虑稳定压缩过程，即C点不会转至C3便可完成结构完全压缩。

完全压缩后斜索AC变化为AC2，△AOC2的两边OA和OC2长度为端面半径R，由此得到：

则在压缩过程中上下端面的相对旋转角度为：

同理可以得出φ1。

家族兴旺，好事。讲个理儿什么的呼啦来一帮，气势上就能把对方压住。要不爷爷动不动就爱跟人显摆：我有四个儿子，四个孙子。说话的时候，那腰板挺得“杠杠”的。好嘛，“男人帮”的家族版原来在北京何老爷子家。瞧瞧，爷爷给孙子们起的名东南西北，他老人家统领四方！

C1点为C点在下底面的投影，C1点坐标如下：

则向量AC1的长度为：

在旋转过程中上下端的面相对转角为：

当张拉整体结构完全压缩，C点落在C2点，杆构件最终位置为线段BC2，C2点坐标，压缩后杆长BC2为：

将杆长压缩率ε定义为：

将P=4，高度h=270 mm，上下底面半径ru=rd=200 mm 的第一种构型代入上式，通过理论计算分别求得偏转角φ2=56.79°，压缩率ε=13.07%。

动力学软件柔性建模如图11 所示，测量杆构件从初始状态至完全压缩长度的变化量如图12所示。通过动力学软件柔性建模，得到偏转角和压缩率分别为56.6°和13.46%。

图11 弯杆张拉整体动力学建模Fig.11 Dynamic modeling of bending bar tensegrity

图12 压缩过程杆构件的长度变化Fig.12 Variation of bar length during compression

搭建相同参数的张拉整体模型，并对它进行轴向压缩，测量偏转角及杆长变化率，如图13 所示。压缩前杆长430 mm，压缩后杆长370 mm，得到杆长压缩率为13.95%，端面相对偏转角为50°。

图13 弯杆张拉整体模型测量数据Fig.13 Measuring data of tensegrity structure model

通过理论计算和动力学柔性建模对相同参数模型进行分析，误差控制在0.5%以内，以验证推导偏转角和杆长压缩率表达式的正确性。微小误差来源是仿真软件本身算法或刚体间的接触力等不确定因素。而实物模型上、下端面的偏转角误差在7°，偏大，可能是由于压缩过程中受拉索被拉长或微小误差的累计造成的。

4.5 压缩过程受力分析

仅得到弯杆张拉整体杆构件两节点间的相对位移长度，还不足得到弯杆的形变情况，形变情况还与弯曲角度，材料基本属性，如弹性模量、泊松比等有关。假设两点间距离为400 mm，形变后的长度为280 mm，分别对90°，120°，150°的弯杆张拉整体进行分析。

实际分析中，弯角通过三维建模后3D 打印得到，其打印材料为光敏树脂。通过施加力得到破坏时可施加的最大外力。图14 采用一端固定一段加力的方式。考虑到测量关节强度，将实际测量过程中的竹杆变为金属铝管。

图14 弯角强度校核Fig.14 Strength check of bending joint

对模型加力至28.5 N，外力关节并未发生变形，此时金属管发生明显弯曲，说明竹结构先发生断裂，弯曲关节强度大于杆构件强度。

在有限元软件中关节材料密度为3×102kg/m3，杨氏模量为2×105MPa，泊松比为0.3，屈服极限为50 MPa，初始应变为零；竹节类材料的参数则与多种因素有关，例如：竹节生长环境、竹节受光照时间、以及生长海拔等，另外对于同一根竹节不同部位的材料参数也是不同的，这里取竹节平均值近似研究。密度为1.1×102kg/m3，杨氏模量为9.6×103，泊松比为0.3，屈服极限为30 MPa，初始应变为零。建立长度相同弯曲角度不同的弯角模型，竹杆最大受力点产生的应力如图15 所示。

图15 不同弯角压缩应力云图Fig.15 Compression stress nephogram at different bend⁃ing angles

由图15 可知，最大应力值发生在弯角关节与杆构件连接处外侧，最大应力值分别为46，44，30 MPa。竹杆材料的屈服极限近似为30 MPa，由仿真结果可知，弯杆所成角度越小，材料越不易发生破坏，轴向折展性越好。其他尺寸的弯杆变形张拉整体可根据相同方法得到弯曲角度极限值，弯角小于极限值时变形四杆张拉整体均可完全压缩。对于上述材料，弯杆角度不大于90°时能够实现轴向完全压缩。

4.6 实物搭建

综合上述分析论述，结合张拉整体现有的理论基础，分别搭建单层和多层四杆弯曲张拉整体模型，实现了它们的轴向压缩与自动展开，如图16 所示。此方法可以用于搭建任意层数的弯杆张拉整体结构。

图16 弯杆张拉整体压缩前后的实物图Fig.16 Photos of bending bar before and after compres⁃sion

5 结论

本文针对张拉整体基本单元的要求，即张拉整体展开前占用空间小，展开后空间利用率高，提出了张拉整体沿轴向进行折叠压缩的必要性。给出了索伸长和杆压缩两种实现轴向折叠的方式，并将两种折叠方式进行对比，进而提出了弯杆张拉整体，得出了弯杆张拉整体节点坐标的构建方式和节点矩阵与各参数间的关系。对弯杆张拉整体压缩过程的杆长压缩率和上、下端面相对偏转角度进行理论推导和动力学柔性建模分析，将两种结果进行对比，验证了结果的正确性。对材料强度进行校核，利用有限元找出特定材料下弯曲角度的极限应力值和受力形变。结果表示：将直杆张拉整体变为弯杆张拉整体，在没有添加其他构件的基础上实现了轴向完全压缩，扩大内部空间的利用率，有效解决了杆构件间的干涉问题。高度为270 mm，半径为200 mm 的四杆张拉整体轴向完全压缩时，上、下底面的相对偏转角为57°，杆长压缩率为13%。利用竹节作为杆构件时，弯曲节点连接两杆构件所成角度不得大于90°。