CN108920789A - 一种柔性航天器多级复合控制动力学建模方法 - Google Patents
一种柔性航天器多级复合控制动力学建模方法 Download PDFInfo
- Publication number
- CN108920789A CN108920789A CN201810634463.1A CN201810634463A CN108920789A CN 108920789 A CN108920789 A CN 108920789A CN 201810634463 A CN201810634463 A CN 201810634463A CN 108920789 A CN108920789 A CN 108920789A
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- node
- input
- finite element
- vector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000002131 composite material Substances 0.000 title claims abstract description 23
- 239000013598 vector Substances 0.000 claims description 87
- 239000011159 matrix material Substances 0.000 claims description 54
- 230000009466 transformation Effects 0.000 claims description 40
- 238000006073 displacement reaction Methods 0.000 claims description 32
- 238000004458 analytical method Methods 0.000 claims description 14
- 238000013016 damping Methods 0.000 claims description 11
- 230000007246 mechanism Effects 0.000 claims description 7
- 239000004576 sand Substances 0.000 claims description 7
- 238000002156 mixing Methods 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 6
- 238000003491 array Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims 4
- 230000010354 integration Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 206010020751 Hypersensitivity Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种柔性航天器多级复合控制动力学建模方法,包括如下步骤:(1)、建立柔性航天器一体化有限元模型;(2)、以主动指向超静平台与载荷和星体连接节点、星体执行机构节点为输入节点,载荷敏感器节点、星体敏感器节点、主动指向超静平台与载荷和星体连接节点为输出节点,建立航天器基于有限元方法的结构动力学方程;(3)、将航天器基于有限元方法的结构动力学方程变换为航天器模态坐标下的结构动力学方程,并改写成航天器的状态空间方程;(4)、对状态空间方程进行输入输出变换,得到以输入节点作用力为输入,以输出节点运动学状态为输出的柔性航天器多级复合控制动力学模型,用于控制系统仿真分析。
Description
技术领域
本发明属于航天器控制领域,涉及一种柔性航天器多级复合控制动力学建模方法。
背景技术
当前航天器对指向控制提出了三超(超高精度超高稳定度超敏捷)的要求,仅基于卫星姿态控制的单级控制技术已经难以满足要求。通过在航天器星体和载荷之间安装具有振动隔离、扰振抑制和精确指向调节的柔性环节,即主动指向超静平台,有效解决了航天器控制中的“稳、快、准”的突出矛盾问题,易于实现未来航天器的超高精度超高稳定度超敏捷控制目标,所述主动指向超静平台可以是由多个带作动器的支杆组成。
为验证添加主动指向超静平台后的控制效果,需要尽可能地利用现有模型,加入主动指向超静平台,进行一体化控制系统分析。然而,在进行建立一体化模型和进行控制系统分析时,面临以下问题。一是,在现有卫星控制系统中,星体平台和载荷平台都为刚性连接,在进行三超控制系统分析时,往往采用刚体模型,不能准确反映出各阶模态的响应。二是,已有的有限元模型为整体模型,无法直接添加主动指向超静平台模型进行控制系统设计与分析。三是,现有控制系统分析方法都只针对整星模型,无法验证添加主动指向超静平台后的控制效果。
基于此需要研究柔性航天器多级复合控制动力学建模方法。
发明内容
本发明解决的技术问题是:针对现有整星有限元模型无法进行三超控制系统分析的问题,提出了一种柔性航天器多级复合控制动力学建模方法。
本发明的技术解决方案是:一种柔性航天器多级复合控制动力学建模方法,所述柔性航天器包括载荷、星体和主动指向超静平台三部分,所述主动指向超静平台为六自由度并联平台,用于连接星体和载荷,该方法包括如下步骤:
(1)、建立包含主动指向超静平台的航天器有限元模型,即航天器一体化有限元模型;
(2)、以主动指向超静平台与载荷和星体连接节点、星体执行机构节点为输入节点,载荷敏感器节点、星体敏感器节点、主动指向超静平台与载荷和星体连接节点为输出节点,建立航天器基于有限元方法的结构动力学方程;
(3)、引入模态坐标变换公式,将航天器基于有限元方法的结构动力学方程变换为航天器模态坐标下的结构动力学方程,并进一步改写成航天器模态坐标下的状态空间方程;
(4)、对模态坐标下的状态空间方程进行输入输出变换,得到以输入节点作用力为输入,以输出节点运动学状态为输出的柔性航天器多级复合控制动力学模型,用于控制系统仿真分析。
所述航天器基于有限元方法的结构动力学方程为:
其中,x为航天器一体化有限元模型所有节点位移向量,为n×1维向量,n为自由度数,n=R×N,N为航天器一体化有限元模型节点数,R为单个节点的自由度,M为航天器一体化有限元模型质量阵;D为航天器一体化有限元模型阻尼阵,G为航天器一体化模型陀螺阵;K为航天器一体化模型刚度阵;v为航天器一体化有限元模型中输入节点的控制力向量,为nv×1维向量,nv为R×Nv,Nv为输入节点的个数,w为输入扰振力向量,为nv×1维向量;Bv与Bw分别为控制力向量和输入扰振力向量输入变换矩阵,y为输出节点的状态向量,C为输出变换矩阵。
所述航天器模态坐标下的结构动力学方程为:
式中,q为模态坐标,为航天器一体化有限元模型的特征向量振型矩阵,由求解的航天器一体化有限元模型的特征方程得到的m(m≤n)个特征值λi对应的m个特征向量构成,i=1~m,m为模态数。
所述航天器模态坐标下的状态空间方程为:
其中:
其中,xs为模态坐标q表示的状态向量,为xs的导数,q为模态坐标,为q的导数,Λ为求解航天器一体化有限元模型的特征方程得到的m(m<<n)个特征值λi,i=1~m构成的对角矩阵,v为航天器一体化有限元模型中输入节点的控制力向量,为nv×1维向量,nv为R×Nv,Nv为输入节点的个数,w为输入扰振力向量,为nv×1维向量,σ、g、bv、bw、c分别为航天器一体化有限元模型模态坐标下的阻尼矩阵、陀螺矩阵、控制力向量输入变换矩阵、扰振力向量输入变换矩阵和输出变换矩阵。
所述步骤(5)的具体实现为:
(5.1)、忽略航天器一体化有限元模型陀螺阵G,设Bv和C均为单位阵,忽略Bw,设定航天器一体化有限元模型阻尼阵σ=ξωn,并将这些量代入模态坐标下的状态空间方程,得到航天器物理坐标下的状态空间方程:
其中,us为每个输入节点的六维作用力和力矩,Xs为模态坐标q表示的状态向量,ys为每个输出节点的位移和速度;
(5.2)、取有限元模型输入节点的三维力和三维力矩作为输入变量u,输出节点位移、角位移、速度、角速度、主动指向超静平台与载荷和星体连接节点的距离作为输出变量y;
(5.3)、将us变换成u两次左乘形式,即
us=B1B2u
(5.4)、将ys变换成y两次左乘形式,即:
y=C1C2ys
(5.5)、将步骤(5.3)和步骤(5.4)变换后的us和ys代入步骤(5.1)得到的物理坐标下的状态空间方程,得到输入节点作用力为输入,以输出节点运动学状态为输出的柔性航天器多级复合控制动力学模型。
所述步骤(5.2)中输入变量u为:
其中,ub为作用于星体执行机构节点的三维力和三维力矩向量,Fa为作用于主动指向超静平台与载荷和星体连接节点的三维力和三维力矩向量,即主动指向超静平台作动杆输出力。
所述步骤(5.2)中输出变量y为:
其中,Xp为载荷敏感器节点的三维位移和角位移向量,Xb为星体敏感器节点的三维位移和角位移向量,为载荷敏感器节点的三维速度和角速度向量,为星体敏感器节点的三维速度和角速度向量,δl为主动指向超静平台与载荷和星体连接节点的距离。
步骤(5.3)中B2变换矩阵表达式为:
式中,I1为单位阵,为R×R的单位阵,R为单个节点的自由度,n1=R+R′,主动指向超静平台与载荷和星体连接节点数相等,R'为主动指向超静平台与载荷和星体连接节点数之和,I2为单位阵。
B1变换矩阵表达式为:
其中,n3=R′×R+2×R,n4=R+R′,假设第Ni个节点上的控制力向量为其中中的元素分别为作用于该节点,方向沿节点坐标系X、Y、Z轴的力,分别为作用于该节点,方向绕节点坐标系X、Y、Z轴的力矩,I1,R×1=[1,0,0,0,0,...,0]T表示将输入变量u中主动指向超静平台作动杆输出力Fa直接作为主动指向超静平台与载荷和星体连接节点X轴方向上的输入力
C2变换矩阵表达式为
其中,I1,1×R=[1,0,0,0,0,...,0],n4=R′+2×R,n5=4×R+R′/2,n6=4×R+2×R′。
本发明与现有技术相比的优点在于:
(1)、与刚体模型相比,本发明中所有部件均采用有限元模型建模,能够准确反映出各部件各阶模态的响应;
(2)、本发明直接以执行机构输出力的作用点(如:主动指向超静平台与载荷和星体连接点)和执行机构输出力矩的作用点(星体执行机构节点)为输入节点,以传感器的测量点(如:载荷敏感器节点、星体敏感器节点、主动指向超静平台与载荷和星体连接点)为输出节点,建立航天器模态坐标下的结构动力学方程,物理意义明确直观;
(3)、本发明对模态坐标下的状态空间方程进行输入输出变换,得到以输入节点作用力为输入,以输出节点运动学状态为输出的柔性航天器多级复合控制动力学模型,可直接用于控制系统仿真分析;
(4)、本发明各级模型均采用物理坐标,可直接进行多级复合控制器设计,提供了一种基于载荷-主动指向超静平台-星体多级复合控制系统分析的有效手段。
附图说明
图1为本发明方法的流程图;
图2为本发明对已有模型改进示意图;
图3为本发明方法一般构型的主动指向超静平台坐标系图;
具体实施方式
以下结合附图与具体实施例对本发明进行详细说明。
如图1所示,本发明提供了一种柔性航天器多级复合控制动力学建模方法,包括如下步骤:
(1)、建立包含主动指向超静平台的航天器有限元模型,即航天器一体化有限元模型;
(1.1)、采用有限元分析方法,对载荷和星体两部分建模,分别得到载荷模型和星体模型;
如图2所示,所述柔性航天器包括载荷(上平台)、星体(下平台)和主动指向超静平台三部分,星体包括了安装于星体的挠性附件所述挠性附件包括太阳帆板、天线。
所述主动指向超静平台为六自由度并联平台,也可以是三自由度并联平台,用于连接星体和载荷。本发明实施例采用6根作动杆组成桁架结构的主动指向超静平台连接星体和载荷。
(1.2)、将满足设计构型和刚度系数要求的主动指向超静平台有限元模型,添加到载荷模型和星体模型之间,共同构成包含主动指向超静平台的航天器有限元模型,即航天器一体化有限元模型。
(2)、将满足设计构型和刚度系数要求的主动指向超静平台有限元模型,添加到载荷模型和星体模型之间,共同构成包含主动指向超静平台的航天器有限元模型,即航天器一体化有限元模型;
所述航天器基于有限元方法的结构动力学方程为:
其中,x为航天器一体化有限元模型所有节点位移向量,为n×1维向量,n为自由度数,n=R×N,N为航天器一体化有限元模型节点数,R为单个节点的自由度,M为航天器一体化有限元模型质量阵,为对称正定矩阵,反映了结构的质量特性;D为航天器一体化有限元模型阻尼阵,G为航天器一体化模型陀螺阵,为反对称矩阵,反映了结构上安装的高速转动部件,如飞轮、控制力矩陀螺等产生的陀螺效应;K为航天器一体化模型刚度阵;K、D分别反映了结构的刚度、阻尼特性,当航天器处于在轨自由状态时,二者均为半正定矩阵;v为航天器一体化有限元模型中输入节点的控制力向量,为nv×1维向量,nv为R×Nv,Nv为输入节点的个数,w为输入扰振力向量,为nv×1维向量;Bv与Bw分别为控制力向量和输入扰振力向量输入变换矩阵,y为输出节点的状态向量,C为输出变换矩阵。
(3)、以主动指向超静平台与载荷和星体连接节点、星体执行机构节点为输入节点,载荷敏感器节点、星体敏感器节点、主动指向超静平台与载荷和星体连接节点为输出节点,建立航天器基于有限元方法的结构动力学方程;
本发明实施例中,航天器一体化有限元模型输出节点包括:载荷敏感器节点、星体敏感器节点、主动指向超静平台与载荷和星体连接节点(即作动杆两端节点)。航天器一体化有限元模型输入节点包括:主动指向超静平台与载荷和星体连接节点、星体执行机构节点。
所述航天器模态坐标下的结构动力学方程为:
式中,q为模态坐标,为航天器一体化有限元模型的特征向量振型矩阵,由求解的航天器一体化有限元模型的特征方程得到的m(m≤n)个特征值λi对应的m个特征向量构成,i=1~m,m为模态数,等于航天器控制所关心的模态阶数。
所述航天器一体化有限元模型的特征方程为:
(4)、引入模态坐标变换公式,将航天器基于有限元方法的结构动力学方程变换为航天器模态坐标下的结构动力学方程,并进一步改写成航天器模态坐标下的状态空间方程;
模态坐标变换公式为:
x=Φq
定义模态坐标下的阻尼矩阵σ、陀螺矩阵g、控制力向量输入变换矩阵bv、扰振力向量输入变换矩阵bw、输出变换矩阵c:
2σ=ΦTDΦ,g=ΦTGΦ
bv=ΦTBv,bw=ΦTBw,c=CΦ
将航天器模态坐标下的结构动力学方程用状态空间形式进行改写,得到航天器模态坐标下的状态空间方程为:
其中:
其中,xs为模态坐标q表示的状态向量,为2m×1向量,为xs的导数,q为模态坐标,为q的导数,Λ为求解航天器一体化有限元模型的特征方程得到的m(m<<n)个特征值λi,i=1~m构成的对角矩阵,v为航天器一体化有限元模型中输入节点的控制力向量,为nv×1维向量,nv为6×Nv,Nv为输入节点的个数,w为输入扰振力向量,为nv×1维向量,σ、g、bv、bw、c分别为航天器一体化有限元模型模态坐标下的阻尼矩阵、陀螺矩阵、控制力向量输入变换矩阵、扰振力向量输入变换矩阵和输出变换矩阵。us为输入控制力和扰振力向量的组合。As为模态坐标下的状态变换矩阵,Bs为模态坐标下的输入变换矩阵,Cs为模态坐标下的输出变换矩阵;ys为模态坐标下的输出向量;
(5)、对模态坐标下的状态空间方程进行输入输出变换,得到以输入节点作用力为输入,以输出节点运动学状态为输出的柔性航天器多级复合控制动力学模型,用于控制系统仿真分析。
(5.1)、忽略航天器一体化有限元模型陀螺阵G,设Bv和C均为单位阵,忽略Bw,设定航天器一体化有限元模型阻尼阵σ=ξωn,并将这些量代入模态坐标下的状态空间方程,得到航天器物理坐标下的状态空间方程:
其中,us为每个输入节点的六维作用力和力矩,xs为模态坐标q表示的状态向量,ys为每个输出节点的位移和速度;ξ为给定的经验值0.005。
根据航天器一体化有限元模型的特征向量振型矩阵Φ,计算得到模态坐标下的控制力向量输入变换矩阵bv、输出变换矩阵c为:
bv=ΦT,c=Φ
us为每个输入节点的六维作用力和力矩,为nv×1维向量,ys为每个输出节点的位移和速度,为2ny×1维向量,ny=输出节点数×6。因此,状态空间表示的是以模态坐标为内部状态,以输入节点作用力为输入,以输出节点运动学状态为输出的系统。
本发明实施例中,在模型仿真时,不需要改变模态坐标。
(5.2)、取有限元模型输入节点的三维力和三维力矩作为输入变量u,输出节点位移、角位移、速度、角速度、主动指向超静平台与载荷和星体连接节点的距离作为输出变量y;
输入变量u为:
其中,ub为作用于星体执行机构节点的三维力和三维力矩向量,Fa为作用于主动指向超静平台与载荷和星体连接节点的三维力和三维力矩向量,即主动指向超静平台作动杆输出力。
输出变量y为:
其中,Xp为载荷敏感器节点的三维位移和角位移向量,Xb为星体敏感器节点的三维位移和角位移向量,为载荷敏感器节点的三维速度和角速度向量,为星体敏感器节点的三维速度和角速度向量,δl为主动指向超静平台与载荷和星体连接节点的距离。
本发明实施例中,有限元模型的输入节点包括1个上平台质心节点、1个下平台质心节点以及12个主动指向超静平台与载荷和星体连接节点,因此us为(2+R')×R=84维向量,其中,R=6为单个节点的自由度,R'=12为主动指向超静平台与载荷和星体连接节点数。
(5.3)、将us变换成u两次左乘形式,即
us=B1B2u
B2变换矩阵表达式为:
式中,I1为单位阵,为R×R的单位阵,I2为R'×R'单位阵。
B1变换矩阵表达式为:
其中,n3=R′×R+2×R’n4=R+R′,假设第Ni个节点上的控制力向量为其中中的元素分别为作用于该节点,方向沿节点坐标系X、Y、Z轴的力,分别为作用于该节点,方向绕节点坐标系X、Y、Z轴的力矩,I1,R×1=[1,0,0,0,0,...,0]T表示将输入变量u中主动指向超静平台作动杆输出力Fa直接作为主动指向超静平台与载荷和星体连接节点X轴方向上的输入力
本发明实施例中,B2的具体形式为
B1是元素为0和1的坐标分配矩阵,具体形式为
其中,I1,6×1=[1,0,0,0,0,0]T表示将作动杆输出力直接作为节点X轴方向上的输入力,因此要求在有限元模型中,将这12个节点的X轴设定为为沿作动杆轴向方向,正方向均指向上平台。
(5.4)、将ys变换成y两次左乘形式,即:
y=C1C2ys
C2变换矩阵表达式为
其中,I1,1×R=[1,0,0,0,0,...,0],n4=R′+2×R,n5=4×R+R′/2,n6=4×R+2×R′
本发明实施例中,三超平台的输出为1个上平台质心节点、1个下平台质心节点以及12个主动指向超静平台与载荷和星体连接节点的六维位移/角位移和六维速度角速度,因此ys为(2+R')×R×2=168维向量。
其中,yij表示第i个节点的第j个自由度的位移/角位移值,表示第i个节点的第j个自由度的速度/角速度值。
对ys进行两次左乘变换得到y,即
y=C1C2ys
其中,C2是将ys中的节点坐标变换为三超平台坐标系中坐标的变换矩阵,具体形式为:
由此可看出,C2是将ys中的12个主动指向超静平台与载荷和星体连接节点的位移/速度和1个上平台质心节点、1个下平台质心节点的六自由度位移/速度/角位移/角速度取出。
C1是元素为0和1的坐标分配矩阵,具体形式为
综上所述,得到输入输出变换,如图3所示。
(5.5)、将步骤(5.3)和步骤(5.4)变换后的us和ys代入步骤(5.1)得到的物理坐标下的状态空间方程,得到输入节点作用力为输入,以输出节点运动学状态为输出的柔性航天器多级复合控制动力学模型。
由上述模态分析可比较清楚地确定各结构的固有振动特性,识别出模态参数,为基于有限元模型的建模、稳定性分析、控制器设计、振动的传递特性提供了基础。同时主动指向超静平台的引入也带来了同/异位控制、解耦控制等问题。通过细致的动力学分析,把上述真实的动力学对象特性纳入的仿真和设计中来。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (9)
1.一种柔性航天器多级复合控制动力学建模方法,所述柔性航天器包括载荷、星体和主动指向超静平台三部分,所述主动指向超静平台为六自由度并联平台,用于连接星体和载荷,其特征在于包括如下步骤:
(1)、建立包含主动指向超静平台的航天器有限元模型,即航天器一体化有限元模型;
(2)、以主动指向超静平台与载荷和星体连接节点、星体执行机构节点为输入节点,载荷敏感器节点、星体敏感器节点、主动指向超静平台与载荷和星体连接节点为输出节点,建立航天器基于有限元方法的结构动力学方程;
(3)、引入模态坐标变换公式,将航天器基于有限元方法的结构动力学方程变换为航天器模态坐标下的结构动力学方程,并进一步改写成航天器模态坐标下的状态空间方程;
(4)、对模态坐标下的状态空间方程进行输入输出变换,得到以输入节点作用力为输入,以输出节点运动学状态为输出的柔性航天器多级复合控制动力学模型,用于控制系统仿真分析。
2.根据权利要求1所述的一种柔性航天器多级复合控制动力学建模方法,其特征在于:所述航天器基于有限元方法的结构动力学方程为:
其中,x为航天器一体化有限元模型所有节点位移向量,为n×1维向量,n为自由度数,n=R×N,N为航天器一体化有限元模型节点数,R为单个节点的自由度,M为航天器一体化有限元模型质量阵;D为航天器一体化有限元模型阻尼阵,G为航天器一体化模型陀螺阵;K为航天器一体化模型刚度阵;v为航天器一体化有限元模型中输入节点的控制力向量,为nv×1维向量,nv为R×Nv,Nv为输入节点的个数,w为输入扰振力向量,为nv×1维向量;Bv与Bw分别为控制力向量和输入扰振力向量输入变换矩阵,y为输出节点的状态向量,C为输出变换矩阵。
3.根据权利要求2所述的一种柔性航天器多级复合控制动力学建模方法,其特征在于所述航天器模态坐标下的结构动力学方程为:
式中,q为模态坐标,为航天器一体化有限元模型的特征向量振型矩阵,由求解的航天器一体化有限元模型的特征方程得到的m(m≤n)个特征值λi对应的m个特征向量构成,i=1~m,m为模态数。
4.根据权利要求1所述的一种柔性航天器多级复合控制动力学建模方法,其特征在于所述航天器模态坐标下的状态空间方程为:
其中:
其中,xs为模态坐标q表示的状态向量,为xs的导数,q为模态坐标,为q的导数,Λ为求解航天器一体化有限元模型的特征方程得到的m(m<<n)个特征值λi,i=1~m构成的对角矩阵,v为航天器一体化有限元模型中输入节点的控制力向量,为nv×1维向量,nv为R×Nv,Nv为输入节点的个数,w为输入扰振力向量,为nv×1维向量,σ、g、bv、bw、c分别为航天器一体化有限元模型模态坐标下的阻尼矩阵、陀螺矩阵、控制力向量输入变换矩阵、扰振力向量输入变换矩阵和输出变换矩阵。
5.根据权利要求1所述的一种柔性航天器多级复合控制动力学建模方法,其特征在于所述步骤(5)的具体实现为:
(5.1)、忽略航天器一体化有限元模型陀螺阵G,设Bv和C均为单位阵,忽略Bw,设定航天器一体化有限元模型阻尼阵σ=ξωn,并将这些量代入模态坐标下的状态空间方程,得到航天器物理坐标下的状态空间方程:
其中,us为每个输入节点的六维作用力和力矩,Xs为模态坐标q表示的状态向量,ys为每个输出节点的位移和速度;
(5.2)、取有限元模型输入节点的三维力和三维力矩作为输入变量u,输出节点位移、角位移、速度、角速度、主动指向超静平台与载荷和星体连接节点的距离作为输出变量y;
(5.3)、将us变换成u两次左乘形式,即
us=B1B2u
(5.4)、将ys变换成y两次左乘形式,即:
y=C1C2ys
(5.5)、将步骤(5.3)和步骤(5.4)变换后的us和ys代入步骤(5.1)得到的物理坐标下的状态空间方程,得到输入节点作用力为输入,以输出节点运动学状态为输出的柔性航天器多级复合控制动力学模型。
6.根据权利要求5所述的一种柔性航天器多级复合控制动力学建模方法,其特征在于所述步骤(5.2)中输入变量u为:
其中,ub为作用于星体执行机构节点的三维力和三维力矩向量,Fa为作用于主动指向超静平台与载荷和星体连接节点的三维力和三维力矩向量,即主动指向超静平台作动杆输出力。
7.根据权利要求5所述的一种柔性航天器多级复合控制动力学建模方法,其特征在于所述步骤(5.2)中输出变量y为:
其中,Xp为载荷敏感器节点的三维位移和角位移向量,Xb为星体敏感器节点的三维位移和角位移向量,为载荷敏感器节点的三维速度和角速度向量,为星体敏感器节点的三维速度和角速度向量,δl为主动指向超静平台与载荷和星体连接节点的距离。
8.根据权利要求5所述的一种柔性航天器多级复合控制动力学建模方法,其特征在于步骤(5.3)中B2变换矩阵表达式为:
式中,I1为单位阵,为R×R的单位阵,R为单个节点的自由度,n1=R+R′,主动指向超静平台与载荷和星体连接节点数相等,R'为主动指向超静平台与载荷和星体连接节点数之和,I2为单位阵。
B1变换矩阵表达式为:
其中,n3=R′×R+2×R,n4=R+R′,假设第Ni个节点上的控制力向量为其中中的元素分别为作用于该节点,方向沿节点坐标系X、Y、Z轴的力,分别为作用于该节点,方向绕节点坐标系X、Y、Z轴的力矩,I1,R×1=[1,0,0,0,0,...,0]T表示将输入变量u中主动指向超静平台作动杆输出力Fa直接作为主动指向超静平台与载荷和星体连接节点X轴方向上的输入力
9.根据权利要求5所述的一种柔性航天器多级复合控制动力学建模方法,其特征在于:C2变换矩阵表达式为
其中,I1,1×R=[1,0,0,0,0,...,0],n4=R′+2×R,n5=4×R+R′/2,n6=4×R+2×R′。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810634463.1A CN108920789B (zh) | 2018-06-20 | 2018-06-20 | 一种柔性航天器多级复合控制动力学建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810634463.1A CN108920789B (zh) | 2018-06-20 | 2018-06-20 | 一种柔性航天器多级复合控制动力学建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108920789A true CN108920789A (zh) | 2018-11-30 |
CN108920789B CN108920789B (zh) | 2022-07-29 |
Family
ID=64421167
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810634463.1A Active CN108920789B (zh) | 2018-06-20 | 2018-06-20 | 一种柔性航天器多级复合控制动力学建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108920789B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885968A (zh) * | 2019-03-13 | 2019-06-14 | 南京航空航天大学 | 一种复杂运动的运动描述建模方法及系统 |
CN111177859A (zh) * | 2019-12-09 | 2020-05-19 | 上海航天控制技术研究所 | 一种桁架天线的动力学等效连续体建模方法 |
CN111638721A (zh) * | 2020-04-28 | 2020-09-08 | 北京控制工程研究所 | 一种航天器三超控制全链路扰动传递验证系统及验证方法 |
CN111781939A (zh) * | 2020-05-11 | 2020-10-16 | 北京控制工程研究所 | 基于航天器三超相互制约与耦合的姿态控制方法及系统 |
CN113032910A (zh) * | 2021-04-13 | 2021-06-25 | 哈尔滨工业大学 | 一种挠性航天器作动器布局优化方法 |
CN113591220A (zh) * | 2021-08-02 | 2021-11-02 | 中国科学院微小卫星创新研究院 | 一种航天器地面姿态动力学建模仿真方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5203199A (en) * | 1990-10-12 | 1993-04-20 | Teledyne Industries, Inc. | Controlled acceleration platform |
CN102759927A (zh) * | 2012-08-03 | 2012-10-31 | 北京理工大学 | 一种使用多级隔振平台提高光学载荷成像质量的方法 |
CN103605834A (zh) * | 2013-10-31 | 2014-02-26 | 上海卫星工程研究所 | 卫星飞轮被动隔振系统设计方法 |
CN103786901A (zh) * | 2014-02-11 | 2014-05-14 | 北京理工大学 | 一种提高航天器姿态控制性能的方法及隔振平台 |
CN105659888B (zh) * | 2012-08-24 | 2014-07-09 | 北京空间飞行器总体设计部 | 一种基于全柔性卫星模型的控制闭环微振动建模与分析方法 |
CN104158431A (zh) * | 2014-08-05 | 2014-11-19 | 中国空间技术研究院 | 一种静电悬浮式天基超静平台系统 |
CN105203132A (zh) * | 2015-10-21 | 2015-12-30 | 中国矿业大学(北京) | 一种谐振式振动陀螺的输出频率检测方法 |
-
2018
- 2018-06-20 CN CN201810634463.1A patent/CN108920789B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5203199A (en) * | 1990-10-12 | 1993-04-20 | Teledyne Industries, Inc. | Controlled acceleration platform |
CN102759927A (zh) * | 2012-08-03 | 2012-10-31 | 北京理工大学 | 一种使用多级隔振平台提高光学载荷成像质量的方法 |
CN105659888B (zh) * | 2012-08-24 | 2014-07-09 | 北京空间飞行器总体设计部 | 一种基于全柔性卫星模型的控制闭环微振动建模与分析方法 |
CN103605834A (zh) * | 2013-10-31 | 2014-02-26 | 上海卫星工程研究所 | 卫星飞轮被动隔振系统设计方法 |
CN103786901A (zh) * | 2014-02-11 | 2014-05-14 | 北京理工大学 | 一种提高航天器姿态控制性能的方法及隔振平台 |
CN104158431A (zh) * | 2014-08-05 | 2014-11-19 | 中国空间技术研究院 | 一种静电悬浮式天基超静平台系统 |
CN105203132A (zh) * | 2015-10-21 | 2015-12-30 | 中国矿业大学(北京) | 一种谐振式振动陀螺的输出频率检测方法 |
Non-Patent Citations (5)
Title |
---|
CHUNSHI FAN等: "Attitude dynamics modeling and analysis of aggregated modular spacecraft with magnetically levitated spherical momentum actuators", 《2016 35TH CHINESE CONTROL CONFERENCE (CCC)》 * |
廖伯瑜等: "现代机械动力学及其工程应用", 《现代机械动力学及其工程应用》 * |
杨依领等: "综合模态控制力下压电致动器的优化布局", 《浙江大学学报(工学版)》 * |
江洁等: "张拉整体结构振动最优控制的作动器优化配置", 《西安电子科技大学学报》 * |
罗青: "航天器飞轮系统微振动特性及隔振方法研究", 《中国博士学位论文全文数据库 工程科技II辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885968A (zh) * | 2019-03-13 | 2019-06-14 | 南京航空航天大学 | 一种复杂运动的运动描述建模方法及系统 |
CN109885968B (zh) * | 2019-03-13 | 2020-11-20 | 南京航空航天大学 | 一种复杂运动的运动描述建模方法及系统 |
CN111177859A (zh) * | 2019-12-09 | 2020-05-19 | 上海航天控制技术研究所 | 一种桁架天线的动力学等效连续体建模方法 |
CN111177859B (zh) * | 2019-12-09 | 2023-09-12 | 上海航天控制技术研究所 | 一种桁架天线的动力学等效连续体建模方法 |
CN111638721A (zh) * | 2020-04-28 | 2020-09-08 | 北京控制工程研究所 | 一种航天器三超控制全链路扰动传递验证系统及验证方法 |
CN111638721B (zh) * | 2020-04-28 | 2023-08-11 | 北京控制工程研究所 | 一种航天器三超控制全链路扰动传递验证系统及验证方法 |
CN111781939A (zh) * | 2020-05-11 | 2020-10-16 | 北京控制工程研究所 | 基于航天器三超相互制约与耦合的姿态控制方法及系统 |
CN111781939B (zh) * | 2020-05-11 | 2023-06-30 | 北京控制工程研究所 | 基于航天器三超相互制约与耦合的姿态控制方法及系统 |
CN113032910A (zh) * | 2021-04-13 | 2021-06-25 | 哈尔滨工业大学 | 一种挠性航天器作动器布局优化方法 |
CN113591220A (zh) * | 2021-08-02 | 2021-11-02 | 中国科学院微小卫星创新研究院 | 一种航天器地面姿态动力学建模仿真方法 |
CN113591220B (zh) * | 2021-08-02 | 2022-05-17 | 中国科学院微小卫星创新研究院 | 一种航天器地面姿态动力学建模仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108920789B (zh) | 2022-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108920789B (zh) | 一种柔性航天器多级复合控制动力学建模方法 | |
Rui et al. | Transfer matrix method for linear multibody system | |
Jia et al. | Maneuver and active vibration suppression of free-flying space robot | |
Alazard et al. | Two-input two-output port model for mechanical systems | |
Chen et al. | A novel method for the dynamic modeling of Stewart parallel mechanism | |
Banerjee et al. | Efficient simulation of large overall motion of beams undergoing large deflection | |
Hu et al. | Recursive dynamics algorithm for multibody systems with variable-speed control moment gyroscopes | |
Zhang et al. | An integrated control and structural design approach for mesh reflector deployable space antennas | |
Li et al. | A task-space form-finding algorithm for tensegrity robots | |
Senda et al. | Methodology for control of a space robot with flexible links | |
Gong et al. | Suppression in any configuration: A versatile coupling improved multi-objective manipulation framework for modular active vibration isolation system | |
Yang et al. | Influences of space perturbations on robotic assembly process of ultra-large structures | |
Bastos et al. | Analysis of open-loop control design and parallel computation for underactuated manipulators | |
Yu et al. | Active control of a 6-DOF space robot with flexible panels using singular perturbation method | |
Nikravesh et al. | Systematic construction of equations of motion for rigid‐flexible multibody systems containing open and closed kinematic loops | |
PEREZ et al. | A flexible appendage model for use in integrated control/structure spacecraft design | |
Jiang et al. | Orbit-attitude-structure coupled modelling method in local translational coordinate frame for multibody systems | |
CN114115308B (zh) | 一种引力波探测航天器姿态自主规划方法 | |
Senba et al. | Vibration reduction by natural frequency optimization for manipulation of a variable geometry truss | |
Morga et al. | Flexible spacecraft model and robust control techniques for attitude maneuvers | |
Liu et al. | Coordinated control for constrained multiple spacecraft system | |
CN111781939B (zh) | 基于航天器三超相互制约与耦合的姿态控制方法及系统 | |
Tao et al. | Attitude control and vibration suppression of flexible spacecraft based on quintic polynomial path planning | |
Jia et al. | A recursive formulation for open-loop gyroelastic multibody dynamics | |
Zakrzhevskii | Optimal slewing of a flexible spacecraft |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |