摘要
针对可折叠薄膜航天器的特点和设计需求,设计一种由杆件、弦绳、薄膜构成的新型可折叠/展开的刚柔耦合多体系统:可折叠张拉整体‑薄膜系统。杆件可以提供充足的系统刚度,通过控制弦绳拉力可实现系统形态变化和形态控制,薄膜折叠技术可明显减小系统折叠形态的外形尺寸。采用非线性有限元方法建立系统动力学模型,用于描述系统运动中的几何非线性特性和薄膜折叠/展开过程中的材料非线性特性。基于准静态控制和H∞控制方法设计一种系统形态展开组合控制策略,仿真结果表明该控制策略可以在系统预应力构型发生切换的情况下实现系统形态展开,并且在展开形态下快速稳定系统并抑制系统振动。
可折叠薄膜航天器是未来大尺寸航天器的发展方向之一,通过系统折叠技术可以减小航天器在运载过程中的外形尺寸,降低对空天载具尺寸和运载能力的要求。薄膜是核心系统组件,可作为光学镜面、太阳帆帆面、雷达天线等任务组件。常规可折叠薄膜航天器一般采用薄膜和支架分离设计概念,以充气/柔性支架或刚性支架支撑薄膜,通过改变支架形态实现系统的折叠和展开。NASA于2017年发射的土壤水分主动被动探测(SMAP)卫星所安装的AstroMesh‑Lite薄膜反射
可折叠张拉整体‑薄膜系统作为一种新型杆件‑弦绳‑薄膜预应力结

图1 典型的4杆张拉整体‑薄膜系统
Fig.1 A representative four‑bar tensegrity‑membrane system
本文针对一个可折叠张拉整体‑薄膜系统,采用非线性有限元方法建立系统动力学模型,用于描述系统形态展开中因系统组件大幅度平动和转动所导致的几何非线性和薄膜因折叠/展开所导致的材料非线性。基于准静态控制和H∞控制方法设计一种组合控制策略,用于实现从折叠形态到展开形态的系统形态可控展开,并采用数值仿真检验控制策略的性能。
本文研究一个2层、每层3杆、21根弦绳、顶层安装1张薄膜的可折叠张拉整体‑薄膜系统,如

图2 可折叠张拉整体‑薄膜系统
Fig.2 A foldable tensegrity‑membrane system
在可折叠张拉整体‑薄膜系统的形态设计和平衡点求解时,可以假设杆件为刚体,杆件的广义坐标如

图3 杆件bij的广义坐标
Fig.3 Generalized coordinates of bar bij
(1) |
式中:为杆件的广义坐标向量,为杆件的广义坐标向量。
由于悬垂线边缘能够减小薄膜在对称角点载荷下的褶皱现

图4 薄膜外形和折痕样式
Fig.4 Membrane shape and folding pattern
为了便于研究薄膜的刚体运动和弹性形变,引入浮动参考系。原点与薄膜角点重合,轴始终通过薄膜角点和,轴与薄膜角点形成的平面垂直,轴由右手定则确定。薄膜的刚体运动可根据浮动参考系的位移和姿态变化确定,该浮动参考系也给定了薄膜的边界条件:在浮动参考系中,角点的位置不变,角点可沿轴产生位移,角点可沿轴和轴分别产生位移和,相应的薄膜角点载荷(、、)如
为了描述可折叠张拉整体‑薄膜系统的运动特性,采用文献[
(2) |
式中:为第二类Piola‑Kirchhoff应力张量的逆变分量,为Green‑Lagrange应变的协变分量,为外力和外力矩作用所产生的虚功。
分别采用四节点薄膜元素和二节点杆元素对薄膜和杆件的动力学特性进行建模,考虑到中包含惯性力和阻尼力产生的虚功,通过应用线性化和迭代计算格式,可将
(3) |
式中:为质量矩阵,为基于Rayleigh阻尼模型的阻尼矩阵,刚度矩阵包含线性应变增量刚度矩阵和非线性应变增量刚度矩阵, 为时刻第k次迭代的节点位移向量,为外力和外力矩向量,为元素应力产生的节点载荷向量。
需要指出的是,
由于薄膜在系统折叠时为松弛状态,在系统展开时为张紧状态,本文采用Miller‑Hedgepeth非线性材料模
(4) |
式中:为最大主应变,为最小主应力。薄膜的本构关系可写为,为应力向量,为应变向量。当薄膜处于松弛状态时,。当薄膜处于张紧状态时,,其中可表示为
(5) |
式中:为弹性模量,为泊松比。当薄膜处于褶皱状态时,,其中可表示为
(6) |
式中:,,为最小主应变。
该非线性有限元模型已用于多种张拉整体‑薄膜系统的形态设计和验
为了实现可折叠张拉整体‑薄膜系统的可控形态展开,本文设计了一种基于准静态控制和H∞控制的混合控制策略。准静态控制基于系统预应力平衡条件,在系统平衡点区域内设计期望形态变化轨迹,进而给出对应的控制输入,即弦绳余长的变化规律。由于准静态控制是一种开环控制方法,在系统形态展开过程中将必然导致系统响应振荡现象。在系统形态展开接近终点时,引入一个H∞输出反馈控制器,不仅可以快速消除系统振动,并且可以提高系统在展开形态下的形态保持精度。
准静态控制采用开环控制概念,根据期望的系统平衡点轨迹参数生成控制输入,使系统沿该平衡点轨迹进行运
由于可折叠张拉整体‑薄膜系统是一种预应力系统,系统平衡点轨迹(即系统形态展开轨迹)需要根据系统平衡条件进行设计。考虑到系统展开过程中薄膜将出现从松弛到张紧的应力变化过程,可折叠张拉整体‑薄膜系统也将出现从张拉整体预应力构型到张拉整体‑薄膜预应力构型的切换现象。需要分别研究这两种预应力构型下的系统平衡条件,确定系统平衡点包线并设计期望的平衡点轨迹。
在张拉整体预应力构型下,系统仅依靠弦绳张力实现预应力平衡,系统平衡形态可由如下的系统平衡方程确
(7) |
式中:为弦绳拉力向量;为
(8) |
式中:代表矩阵的零空间,为向量T的第j个元素。所有满足
在张拉整体‑薄膜构型下,薄膜和弦绳共同产生张力实现系统预应力平衡,可通过
(9) |
式中: ;矩阵的元素可表示为,为向量的第j个元素,,为第k个弦绳的长度,、和为薄膜角点位移。根据
(10) |
式中:代表矩阵的零空间,为向量F的第j个元素。所有满足
在可折叠张拉整体‑薄膜系统的平衡点区域内,可选取并设计期望的形态展开轨迹,进而确定对应的弦绳拉力/余长变化规律。基于准静态控制概念,通过缓慢控制弦绳拉力变化,实现系统沿期望的形态展开轨迹进行运动,但系统响应将出现振荡现象。
为了消除准静态控制导致的系统响应振荡现象并提高系统在展开形态下的形态保持精度,采用H∞输出反馈控制器实现系统在展开形态的闭环控制。
在系统展开形态下,系统线性模型可写为
(11) |
式中:,为薄膜上m个点处的横向挠度向量;控制输入向量的各个元素为弦绳余长;向量、和分别为广义输入向量、广义误差向量和测量输出向量。
H∞输出反馈控制器 的设计等价于求解如下的凸函数优化问
(12) |
式中
(13) |
控制器参数、、和可通过
(14) |
H∞控制器设计结构如
(15) |

图5 H∞控制器设计结构图
Fig.5 System interconnections for H∞ control design
本节根据非线性有限元模型,采用数值仿真方法研究可折叠张拉整体‑薄膜系统的形态展开,并检验所设计的组合控制策略的控制性能。
根据
(16) |
式中:,。该二阶系统可使期望轨迹具有平滑的变化规律,避免系统加速度的突变现象。

图6 系统平衡点区域和参考轨迹关键点
Fig.6 System equilibrium regions and key equilibriums
可折叠张拉整体‑薄膜系统的形态展开仿真结果如图

图7 杆件b11的响应曲线
Fig.7 Responses of bar b11

图8 杆件b21的响应曲线
Fig.8 Responses of bar b21

图9 薄膜横向挠度
Fig.9 Membrane transverse deflections

图10 弦绳拉力的变化曲线
Fig.10 Representative tendon force variations
从

图11 可折叠张拉整体‑薄膜系统展开过程
Fig.11 Deployment sequence of a foldable tensegrity‑membrane system
本文针对可折叠张拉整体‑薄膜系统,采用非线性有限元方法建立系统的动力学模型,结合准静态控制和H∞控制方法设计了系统形态展开的组合控制策略,获得了以下结论:
(1)非线性有限元方法能够描述可折叠张拉整体‑薄膜系统在形态展开过程中的几何非线性和材料非线性力学特性。
(2)准静态控制方法能够有效控制可折叠张拉整体‑薄膜系统的形态展开过程,使系统沿期望的形态展开轨迹进行运动,但系统响应中存在明显振荡。
(3)H∞控制方法能够在系统展开过程的终端实现系统快速稳定并消除系统振动,提高系统形态保持精度。
参考文献
ALVAREZ‑SALAZAR O S, ADAMS D, MILMAN M, et al. Pointing architecture of SMAP’s large spinning antenna[C]//Proceedings of Guidance, Navigation, and Control and Co‑located Conferences. [S.l.]: AIAA, 2013. [百度学术]
THOMSON M W. AstroMesh deployable reflectors for Ku‑ and Ka‑ band commercial satellites[C]//Proceedings of the 20th AIAA International Communication Satellite Systems Conference and Exhibit. [S.l.]: AIAA, 2002. [百度学术]
WALLER D, CAMPBELL L, BOMBER J, et al. MOIRE primary diffractive optical element structure deployment testing[C]//Proceedings of AIAA SciTech Forum. [S.l.]: AIAA, 2015. [百度学术]
TANDY W D, ATCHESON P, DOMBER J, et al. MOIRE gossamer space telescope—Testing challenges and solutions[C]//Proceedings of AIAA SciTech Forum. [S.l.]: AIAA, 2013. [百度学术]
FOOTDALE J N, MURPHEY T W, PETERSON M. Design and testing of self‑deploying membrane optic support structure using rollable composite tape springs[C]//Proceedings of the 54th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference. [S.l.]: AIAA, 2013. [百度学术]
CADOGAN D P, LIN J K, GRAHNE M S. The development of inflatable space radar reflectarrays[C]//Proceedings of the 40th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference. [S.l.]: AIAA, 1999. [百度学术]
FREELAND R E, VEAL G R. Significance of the inflatable antenna experiment technology[C]//Proceedings of the 39th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference and Exhibit. [S.l.]: AIAA, 1998. [百度学术]
YANG S, SULTAN C. Deployment of foldable tensegrity‑membrane systems via transition between tensegrity configurations and tensegrity‑membrane configurations[J]. International Journal of Solids and Structures, 2019, 160: 103‑119. [百度学术]
YANG S, SULTAN C. Modeling of tensegrity⁃membrane systems[J]. International Journal of Solids and Structures, 2016, 82: 125‑143. [百度学术]
FANG H, LOU M, HSI L, et al. Catenary systems for membrane structures[C]//Proceedings of the 42nd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference and Exhibit. [S.l.]: AIAA, 2001. [百度学术]
GUEST S D, PELLEGRINO S. Inextensional wrapping of flat membranes[C]//Proceedings of the First International Seminar on Structural Morphology, Montpellier. France: [s.n.], 1992. [百度学术]
REYNOLDS W D, MURPHEY T W. Elastic spiral folding for flat membrane apertures[C]//Proceedings of AIAA SciTech Forum. [S.l.]: AIAA, 2014. [百度学术]
BATHE K J. Finite element procedures[M]. Upper Saddle River, New Jersey: Prentice‑Hall, 1996. [百度学术]
MILLER R K, HEDGEPETH J M. An algorithm for finite element analysis of partly wrinkled membranes [J]. AIAA Journal, 1982, 20(12): 1761‑1763. [百度学术]
YANG S, SULTAN C. A comparative study on the dynamics of tensegrity‑membrane systems based on multiple models[J]. International Journal of Solids and Structures, 2017, 113/114: 47‑69. [百度学术]
YANG S, SULTAN C. LPV control of a tensegrity‑membrane system[J]. Mechanical Systems and Signal Processing, 2017, 95: 397‑424. [百度学术]
SULTAN C. Close tracking of equilibrium paths[J]. International Journal of Robust and Nonlinear Control, 2018, 28: 2209‑2230. [百度学术]
MANKALA K K, AGRAWAL S K. Equilibrium‑to⁃equilibrium maneuvers of flexible electrodynamic tethers in equatorial orbits[J]. Journal of Spacecraft and Rockets, 2006, 43(3): 651‑658. [百度学术]
STEPHANE R, LUU T, BERUBE A. Space shuttle testing of the TriDAR 3D rendezvous and docking sensor[J]. Journal of Field Robotics, 2012, 29(4): 535‑553. [百度学术]
STEPHENS J P, VOS G A, BILIMORIA K D, et al. Orion handling qualities during international space station proximity operations and docking[J]. Journal of Spacecraft and Rockets, 2013, 50(2): 449‑457. [百度学术]
DU W, MA S, LI B, et al. Force analytic method for rolling gaits of tensegrity robots[J]. IEEE/ASME Transactions on Mechatronics, 2016, 21(5): 2249‑2259. [百度学术]
DENG J, LIU Y, CHEN W, et al. A XY transporting and nanopositioning piezoelectric robot operated by leg rowing mechanism[J]. IEEE/ASME Transactions on Mechatronics, 2019, 24(1): 207‑217. [百度学术]
PORTA J M, HERNANDEZ‑JUAN S. Path planning for active tensegrity structures[J]. International Journal of Solids and Structures, 2016, 78/79: 47‑56. [百度学术]
GAHINET P, APKARIAN P. A linear matrix inequality approach to H∞ control[J]. International Journal of Robust and Nonlinear Control, 1994, 4(4): 421‑448. [百度学术]
IWASAKI T, SKELTON R E. All controllers for the general H∞ control problem: LMI existence conditions and state space formulas[J]. Automatica, 1994, 30(8): 1307‑1317. [百度学术]