网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

可折叠张拉整体‑薄膜系统的可控形态展开  PDF

  • 杨庶
西北工业大学自动化学院,西安710072

中图分类号: V414

最近更新:2021-10-28

DOI:10.16356/j.1005-2615.2021.S.020

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

针对可折叠薄膜航天器的特点和设计需求,设计一种由杆件、弦绳、薄膜构成的新型可折叠/展开的刚柔耦合多体系统:可折叠张拉整体‑薄膜系统。杆件可以提供充足的系统刚度,通过控制弦绳拉力可实现系统形态变化和形态控制,薄膜折叠技术可明显减小系统折叠形态的外形尺寸。采用非线性有限元方法建立系统动力学模型,用于描述系统运动中的几何非线性特性和薄膜折叠/展开过程中的材料非线性特性。基于准静态控制和H控制方法设计一种系统形态展开组合控制策略,仿真结果表明该控制策略可以在系统预应力构型发生切换的情况下实现系统形态展开,并且在展开形态下快速稳定系统并抑制系统振动。

可折叠薄膜航天器是未来大尺寸航天器的发展方向之一,通过系统折叠技术可以减小航天器在运载过程中的外形尺寸,降低对空天载具尺寸和运载能力的要求。薄膜是核心系统组件,可作为光学镜面、太阳帆帆面、雷达天线等任务组件。常规可折叠薄膜航天器一般采用薄膜和支架分离设计概念,以充气/柔性支架或刚性支架支撑薄膜,通过改变支架形态实现系统的折叠和展开。NASA于2017年发射的土壤水分主动被动探测(SMAP)卫星所安装的AstroMesh‑Lite薄膜反射[

1‑2]和DARPA主导的薄膜光学即时成像器(MOIRE)项目地面样[3‑4]是典型的基于刚性支架的可折叠薄膜航天器,其形态折叠和展开依赖于大量附加机械连接组件和作动器,不仅使得系统增重,系统展开失败的风险也将增加。AFRL研制的立方星衍射望远[5]、JPL和ILC Dover Inc.研制的太空充气雷达阵[6]、NASA充气天线实验中的充气薄膜天线[7]是典型的基于充气/柔性支架的可折叠薄膜航天器。充气/柔性支架的引入使得航天器折叠形态具有极小的外形尺寸,但在展开过程中缺乏足够的系统刚性,无法实现系统形态的精确控制,极易导致系统展开失[5,7]

可折叠张拉整体‑薄膜系统作为一种新型杆件‑弦绳‑薄膜预应力结[

8],杆件为系统提供足够的刚度,薄膜和弦绳产生的拉力通过预应力平衡使系统保持指定形态,弦绳作为内置作动器可以通过拉力/余长控制实现系统形态控制和折叠/展开。同时,薄膜折叠技术的引入,使得该系统具有更为紧凑的折叠形态,放宽了张拉整体‑薄膜系统中薄膜需要始终张紧的力学约[9]。显然,作为张拉整体‑薄膜系统(图1)的技术延伸和拓展,可折叠张拉整体‑薄膜系统这种新型刚柔耦合多体系统能够满足可折叠薄膜航天器的设计需求,克服常规可折叠薄膜航天器的诸多缺陷。

图1 典型的4杆张拉整体‑薄膜系统

Fig.1 A representative four‑bar tensegrity‑membrane system

本文针对一个可折叠张拉整体‑薄膜系统,采用非线性有限元方法建立系统动力学模型,用于描述系统形态展开中因系统组件大幅度平动和转动所导致的几何非线性和薄膜因折叠/展开所导致的材料非线性。基于准静态控制和H控制方法设计一种组合控制策略,用于实现从折叠形态到展开形态的系统形态可控展开,并采用数值仿真检验控制策略的性能。

1 可折叠张拉整体‑薄膜系统

本文研究一个2层、每层3杆、21根弦绳、顶层安装1张薄膜的可折叠张拉整体‑薄膜系统,如图2所示。第一层杆件(b1j)的下端点(A1j)通过球形铰链与系统底座相连,第二层杆件(b2j)的上端点(B2j)与薄膜的角点连接。惯性参考系oxyz的原点o与杆件b13的下端点A13重合,ox轴通过点A13和点A11oz轴与系统底座垂直,oy轴可通过右手定则确定。

图2 可折叠张拉整体‑薄膜系统

Fig.2 A foldable tensegrity‑membrane system

在可折叠张拉整体‑薄膜系统的形态设计和平衡点求解时,可以假设杆件为刚体,杆件bij的广义坐标如图3所示。对于本文研究的可折叠张拉整体‑薄膜系统,由于薄膜角点与第二层杆件上端点相连,并且弦绳两端分别与杆件端点相连,系统形态可由独立广义坐标向量qR21确定。

图3 杆件bij的广义坐标

Fig.3 Generalized coordinates of bar bij

q=[q11T  q12T  q13T  q21T  q22T  q23T]T (1)

式中:q1j=[θ1j  ψ1j]T为杆件b1j的广义坐标向量,q2j=[x2j  y2j  z2j  θ2j  ψ2j]T为杆件b2j的广义坐标向量。

由于悬垂线边缘能够减小薄膜在对称角点载荷下的褶皱现[

10],因此薄膜设计为具有悬垂线边缘的正三角形。采用文献[11‑12]中的薄膜螺旋折叠方法进行薄膜折叠,薄膜外形和折痕样式如图4所示。

图4 薄膜外形和折痕样式

Fig.4 Membrane shape and folding pattern

为了便于研究薄膜的刚体运动和弹性形变,引入浮动参考系ofxfyfzf。原点of与薄膜角点B23重合,ofxf轴始终通过薄膜角点B23B21ofzf轴与薄膜角点形成的平面垂直,ofyf轴由右手定则确定。薄膜的刚体运动可根据浮动参考系ofxfyfzf的位移和姿态变化确定,该浮动参考系也给定了薄膜的边界条件:在浮动参考系ofxfyfzf中,角点B23的位置不变,角点B21可沿ofxf轴产生位移Δ1u,角点B22可沿ofxf轴和ofyf轴分别产生位移Δ2uΔ2v,相应的薄膜角点载荷(F1uF2uF2v)如图4所示。

2 非线性有限元模型

为了描述可折叠张拉整体‑薄膜系统的运动特性,采用文献[

13]中的完整拉格朗日计算格式建立系统的非线性有限元模型。系统在t+Δt时刻的虚功可表达为

V          0t+ΔtS˜ijδε˜          0t+Δtε˜ijdV0=Rt+Δt (2)

式中:          0t+ΔtS˜ij为第二类Piola‑Kirchhoff应力张量的逆变分量,          0t+Δtε˜ij为Green‑Lagrange应变的协变分量,t+ΔtR为外力和外力矩作用所产生的虚功。

分别采用四节点薄膜元素和二节点杆元素对薄膜和杆件的动力学特性进行建模,考虑到t+ΔtR中包含惯性力和阻尼力产生的虚功,通过应用线性化和迭代计算格式,可将式(2)转化为如下的非线性有限元计算方程

MU¨t+ΔtU¨(k)+DU˙t+ΔtU˙(k)+KtK(k-1)ΔU(k)=t+ΔtR+Ft+ΔtF(k-1)t+ΔtU(k)=Ut+ΔtU(k-1)+ΔU(k)t+ΔtU(0)=Ut (3)

式中:M为质量矩阵,D=αM+βKtK(k-1)为基于Rayleigh阻尼模型的阻尼矩阵,刚度矩阵tK(k-1)=KtKL(k-1)+KtKNL(k-1)包含线性应变增量刚度矩阵tKL(k-1)和非线性应变增量刚度矩阵tKNL(k-1)t+ΔtU(k)t+Δt时刻第k次迭代的节点位移向量,t+ΔtR为外力和外力矩向量,t+ΔtF(k-1)为元素应力产生的节点载荷向量。

需要指出的是,式(3)中的非线性有限元模型在建立杆件动力学特性时,考虑了杆件的刚体运动和轴向弹性形变。相对于系统形态设计和平衡点求解所采用的杆件刚体假设,非线性有限元模型的精度和复杂程度明显提高,可以用于检验系统形态设计结果的准确性。

由于薄膜在系统折叠时为松弛状态,在系统展开时为张紧状态,本文采用Miller‑Hedgepeth非线性材料模[

14],建立薄膜在松弛、张紧、褶皱状态下的非线性本构关系。薄膜各点处的载荷状态可根据式(4)中的应力‑应变判别准则进行判断。

σ2>0ε10ε1>0 & σ20 (4)

式中:ε1为最大主应变,σ2为最小主应力。薄膜的本构关系可写为τ=Dmετ=τ11τ22τ12T为应力向量,ε=[ε11ε22γ12]T为应变向量。当薄膜处于松弛状态时,Dm=0。当薄膜处于张紧状态时,Dm=DT,其中DT可表示为

DT=E1-μ21μ0μ1000(1-μ)/2 (5)

式中:E为弹性模量,μ为泊松比。当薄膜处于褶皱状态时,Dm=DW,其中DW可表示为

DW=E421+P0Qμ21-PQQQ1 (6)

式中:P=(ε11-ε22)/(ε1-ε2)Q=γ12/(ε1-ε2)ε2为最小主应变。

该非线性有限元模型已用于多种张拉整体‑薄膜系统的形态设计和验[

9]、力学分[15]、形态控制仿[16],详细建模方法参见文献[8‑9, 13]。

3 系统展开的组合控制策略

为了实现可折叠张拉整体‑薄膜系统的可控形态展开,本文设计了一种基于准静态控制和H控制的混合控制策略。准静态控制基于系统预应力平衡条件,在系统平衡点区域内设计期望形态变化轨迹,进而给出对应的控制输入,即弦绳余长的变化规律。由于准静态控制是一种开环控制方法,在系统形态展开过程中将必然导致系统响应振荡现象。在系统形态展开接近终点时,引入一个H输出反馈控制器,不仅可以快速消除系统振动,并且可以提高系统在展开形态下的形态保持精度。

3.1 准静态控制

准静态控制采用开环控制概念,根据期望的系统平衡点轨迹参数生成控制输入,使系统沿该平衡点轨迹进行运[

17]。该控制方法是一种低成本控制方法,在航天器控[18‑20]、机器人控[21‑22]、结构控[17,23]等方面得到应用。由于该控制方法不引入反馈信号,系统在运动中将出现响应振荡现象,并且一般需要系统进行慢速运动以保证控制性能。

由于可折叠张拉整体‑薄膜系统是一种预应力系统,系统平衡点轨迹(即系统形态展开轨迹)需要根据系统平衡条件进行设计。考虑到系统展开过程中薄膜将出现从松弛到张紧的应力变化过程,可折叠张拉整体‑薄膜系统也将出现从张拉整体预应力构型到张拉整体‑薄膜预应力构型的切换现象。需要分别研究这两种预应力构型下的系统平衡条件,确定系统平衡点包线并设计期望的平衡点轨迹。

在张拉整体预应力构型下,系统仅依靠弦绳张力实现预应力平衡,系统平衡形态可由如下的系统平衡方程确[

8]

A(q)T=0 (7)

式中:T为弦绳拉力向量;q式(1)中定义的系统广义坐标向量;矩阵A(q)的元素可表示为Aij=lj/qilj为第j个弦绳的长度,qi为向量q的第i个元素。系统处于平衡形态时,弦绳拉力应大于0。对于本文研究的可折叠张拉整体‑薄膜系统,矩阵A(q)为方阵,根据式(7)可将系统在张拉整体预应力构型下的平衡条件写为如下形式

detA(q)=0TN(A(q))Tj>0 (8)

式中:N(A(q))代表矩阵A(q)的零空间,Tj为向量T的第j个元素。所有满足式(8)的系统平衡点构成了系统在张拉整体预应力构型下的平衡点区域。

在张拉整体‑薄膜构型下,薄膜和弦绳共同产生张力实现系统预应力平衡,可通过式(9)确定系统平衡形 [

8]

B(q)F=0 (9)

式中: F=TTF1uF2uF2vT;矩阵B(q)的元素可表示为Bij=dj/qidj为向量d的第j个元素,d=l1l21Δ1uΔ2uΔ2vTlk为第k个弦绳的长度,Δ1uΔ2uΔ2v为薄膜角点位移。根据式(9),系统在张拉整体‑薄膜预应力构型下的平衡条件可表示为

FN(B(q))Fj>0 (10)

式中:N(B(q))代表矩阵B(q)的零空间,Fj为向量F的第j个元素。所有满足式(10)的系统平衡点构成了系统在张拉整体‑薄膜预应力构型下的平衡点区域。系统平衡条件(式(7~10))的详细推导过程参见文献[

8‑9]。

在可折叠张拉整体‑薄膜系统的平衡点区域内,可选取并设计期望的形态展开轨迹,进而确定对应的弦绳拉力/余长变化规律。基于准静态控制概念,通过缓慢控制弦绳拉力变化,实现系统沿期望的形态展开轨迹进行运动,但系统响应将出现振荡现象。

3.2 H控制

为了消除准静态控制导致的系统响应振荡现象并提高系统在展开形态下的形态保持精度,采用H输出反馈控制器实现系统在展开形态的闭环控制。

在系统展开形态下,系统线性模型可写为

x˙ey=AB1B2C1D11D12C2D21D22xdu (11)

式中:x=[qTwTq˙Tw˙T]TwRm为薄膜上m个点处的横向挠度向量;控制输入向量u的各个元素为弦绳余长;向量dey分别为广义输入向量、广义误差向量和测量输出向量。

H输出反馈控制器K(s)=Ck(sI-Ak)Bk+Dk 的设计等价于求解如下的凸函数优化问[

24‑25]

find R, S, A¯k, B¯k, C¯k, Dks.t.RIIS0ΨA+ΨATΨBΨCTΨBT-γIΨDTΨCΨD-γI0 (12)

式中

ΨA=AR+B2C¯kA+B2DkC2A¯kSA+B¯kC2ΨB=B1+B2DkD21SB1+B¯kD21ΨC=[C1R+D12C¯kC1+D12DkC2]ΨD=D11+D12DkD21 (13)

控制器参数AkBkCkDk可通过式(14)确定

MNT=I-RSB¯k=NBk+SB2DkC¯k=CkMT+DkC2RA¯k=NAkMT+NBkC2R+SB2CkMT+       S(A+B2DkC2)R (14)

H控制器设计结构如图5所示,权重函数设计为如下形式

We=diag0.1s+10s+0.1I21×21, 0.05s+5s+0.25Im×mWu=1/0.05I21×21Wd1=Wd2=0.001I21×21Wd3=Wd4=0.001Im×m (15)

图5 H控制器设计结构图

Fig.5 System interconnections for H control design

4 仿真结果

本节根据非线性有限元模型,采用数值仿真方法研究可折叠张拉整体‑薄膜系统的形态展开,并检验所设计的组合控制策略的控制性能。

根据式(8)式(10)中两种系统预应力构型的平衡条件,基于系统对称条[

9],通过数值计算方法可确定系统平衡点区域,如图6所示。系统形态展开参考轨迹根据轨迹关键点E0E1E2E3生成,其中E0对应系统折叠形态,E3对应系统展开形态。当系统处于轨迹关键点E1时,薄膜展开但未张紧;当系统处于轨迹关键点E2时,薄膜张紧。在轨迹分段E2E3上,薄膜具有相同的应力状态。通过为每个轨迹关键点指定轨迹生成坐标rtraj的数值,可通过如下的二阶系统生成期望的系统形态展开轨迹。

r¨traj+2ξωr˙traj+ω2rtraj=ω2rtraj(0)=0r˙traj(0)=0 (16)

式中:ω=0.15ξ=1.0。该二阶系统可使期望轨迹具有平滑的变化规律,避免系统加速度的突变现象。

图6 系统平衡点区域和参考轨迹关键点

Fig.6 System equilibrium regions and key equilibriums

可折叠张拉整体‑薄膜系统的形态展开仿真结果如图7~10所示。

图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

图7图8中的杆件b11和杆件b21的响应曲线可以看出,准静态控制方法可以使系统跟随期望形态展开轨迹,但系统响应中存在明显的振荡。当系统从张拉整体预应力构型运动至张拉整体‑薄膜预应力构型时,由于薄膜出现应力突变,导致系统响应中的振荡幅值增大。当系统运动至展开形态时,控制器切换为H控制器,实现了系统的快速稳定和形态保持,可以看出该H控制器具有良好的控制性能。

图9所示的是系统进入张拉整体‑薄膜预应力构型后的薄膜横向挠度响应曲线,可以看出准静态控制方法使得薄膜出现持续的振动。当控制器切换为H控制器后,薄膜振动得到有效控制。从图10所示的弦绳拉力变化曲线可以看出,弦绳拉力在系统形态展开过程中均大于0,满足系统平衡条件。

图11所示的是系统展开过程的可视化结果,可以看出本文设计的组合控制策略能够控制可折叠张拉整体‑薄膜系统的形态变化,并使薄膜从折叠状态运动到展开状态。

图11 可折叠张拉整体‑薄膜系统展开过程

Fig.11 Deployment sequence of a foldable tensegrity‑membrane system

5 结  论

本文针对可折叠张拉整体‑薄膜系统,采用非线性有限元方法建立系统的动力学模型,结合准静态控制和H控制方法设计了系统形态展开的组合控制策略,获得了以下结论:

(1)非线性有限元方法能够描述可折叠张拉整体‑薄膜系统在形态展开过程中的几何非线性和材料非线性力学特性。

(2)准静态控制方法能够有效控制可折叠张拉整体‑薄膜系统的形态展开过程,使系统沿期望的形态展开轨迹进行运动,但系统响应中存在明显振荡。

(3)H控制方法能够在系统展开过程的终端实现系统快速稳定并消除系统振动,提高系统形态保持精度。

参考文献

1

ALVAREZ‑SALAZAR O SADAMS DMILMAN Met al. Pointing architecture of SMAP’s large spinning antenna[C]//Proceedings of Guidance, Navigation, and Control and Co‑located Conferences. [S.l.]AIAA2013. [百度学术

2

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.]AIAA2002. [百度学术

3

WALLER DCAMPBELL LBOMBER Jet al. MOIRE primary diffractive optical element structure deployment testing[C]//Proceedings of AIAA SciTech Forum. [S.l.]AIAA2015. [百度学术

4

TANDY W DATCHESON PDOMBER Jet al. MOIRE gossamer space telescope—Testing challenges and solutions[C]//Proceedings of AIAA SciTech Forum. [S.l.]AIAA2013. [百度学术

5

FOOTDALE J NMURPHEY T WPETERSON 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.]AIAA2013. [百度学术

6

CADOGAN D PLIN J KGRAHNE 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.]AIAA1999. [百度学术

7

FREELAND R EVEAL 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.]AIAA1998. [百度学术

8

YANG SSULTAN C. Deployment of foldable tensegrity‑membrane systems via transition between tensegrity configurations and tensegrity‑membrane configurations[J]. International Journal of Solids and Structures2019160103119. [百度学术

9

YANG SSULTAN C. Modeling of tensegrity⁃membrane systems[J]. International Journal of Solids and Structures201682125143. [百度学术

10

FANG HLOU MHSI Let 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.]AIAA2001. [百度学术

11

GUEST S DPELLEGRINO S. Inextensional wrapping of flat membranes[C]//Proceedings of the First International Seminar on Structural Morphology, Montpellier. France[s.n.]1992. [百度学术

12

REYNOLDS W DMURPHEY T W. Elastic spiral folding for flat membrane apertures[C]//Proceedings of AIAA SciTech Forum. [S.l.]AIAA2014. [百度学术

13

BATHE K J. Finite element procedures[M]. Upper Saddle River, New JerseyPrentice‑Hall1996. [百度学术

14

MILLER R KHEDGEPETH J M. An algorithm for finite element analysis of partly wrinkled membranes [J]. AIAA Journal19822012): 17611763. [百度学术

15

YANG SSULTAN C. A comparative study on the dynamics of tensegrity‑membrane systems based on multiple models[J]. International Journal of Solids and Structures2017113/1144769. [百度学术

16

YANG SSULTAN C. LPV control of a tensegrity‑membrane system[J]. Mechanical Systems and Signal Processing201795397424. [百度学术

17

SULTAN C. Close tracking of equilibrium paths[J]. International Journal of Robust and Nonlinear Control20182822092230. [百度学术

18

MANKALA K KAGRAWAL S K. Equilibrium‑to⁃equilibrium maneuvers of flexible electrodynamic tethers in equatorial orbits[J]. Journal of Spacecraft and Rockets2006433): 651658. [百度学术

19

STEPHANE RLUU TBERUBE A. Space shuttle testing of the TriDAR 3D rendezvous and docking sensor[J]. Journal of Field Robotics2012294): 535553. [百度学术

20

STEPHENS J PVOS G ABILIMORIA K Det al. Orion handling qualities during international space station proximity operations and docking[J]. Journal of Spacecraft and Rockets2013502): 449457. [百度学术

21

DU WMA SLI Bet al. Force analytic method for rolling gaits of tensegrity robots[J]. IEEE/ASME Transactions on Mechatronics2016215): 22492259. [百度学术

22

DENG JLIU YCHEN Wet al. A XY transporting and nanopositioning piezoelectric robot operated by leg rowing mechanism[J]. IEEE/ASME Transactions on Mechatronics2019241): 207217. [百度学术

23

PORTA J MHERNANDEZ‑JUAN S. Path planning for active tensegrity structures[J]. International Journal of Solids and Structures201678/794756. [百度学术

24

GAHINET PAPKARIAN P. A linear matrix inequality approach to H control[J]. International Journal of Robust and Nonlinear Control199444): 421448. [百度学术

25

IWASAKI TSKELTON R E. All controllers for the general H control problem: LMI existence conditions and state space formulas[J]. Automatica1994308): 13071317. [百度学术

您是第位访问者
网站版权 © 南京航空航天大学学报
技术支持:北京勤云科技发展有限公司
请使用 Firefox、Chrome、IE10、IE11、360极速模式、搜狗极速模式、QQ极速模式等浏览器,其他浏览器不建议使用!