一种考虑反作用飞轮特性及功率受限的航天器姿态控制方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN110502025A(43)申请公布日2019.11.26(21)申请号201910707553.3(22)申请日2019.08.01(71)申请人北京航空航天大学地址100191北京市海淀区学院路37号(72)发明人胡庆雷张鑫鑫董宏洋郑建英(51)Int.Cl.G05D1/08(2006.01)权利要求书6页说明书13页附图1页(54)发明名称一种考虑反作用飞轮特性及功率受限的航天器姿态控制方法(57)摘要本发明公开了一种考虑反作用飞轮特性及功率受限的航天器姿态控制方法，包括以下步骤：建立四元数描述方式下的航天器运动学和动力学模型；建立航天器体装式太阳能帆板输出功率与太阳入射角的关系；考虑反作用飞轮特性，建立角动量约束、控制力矩约束以及体装式太阳能帆板故障时反作用飞轮的实际可用功率约束；根据任务需求设计目标性能函数；基于模型预测控制策略完成上述约束条件下的航天器姿态控制任务。本方法能够很好的解决航天器在考虑反作用飞轮特性以及体装式太阳能帆板发生故障时的姿态控制问题，通过对目标函数的设计达到控制精度和能量消耗的综合最优，保证航天器在反作用飞轮特性及功率受限的情况下仍能完成高精度姿态控制任务。CN110502025ACN110502025A权利要求书1/6页1.一种考虑反作用飞轮特性及功率受限的航天器姿态控制方法，其特征在于，包括以下步骤：S1：针对航天器姿态控制问题，建立四元数描述方式下的航天器姿态运动学和动力学模型，将所述航天器姿态运动学与动力学模型离散化，建立模型预测控制方法的预测模型；S2：建立航天器体装式太阳能帆板输出功率与太阳入射角的关系，并将其转化为与航天器姿态角有关的表达式；S3：基于航天器执行机构反作用飞轮特性，建立反作用飞轮的角动量约束和控制力矩约束，同时建立航天器体装式太阳能帆板故障时反作用飞轮的实际可用功率约束；S4：根据任务需求设计目标性能函数；S5：基于步骤S1中建立的模型预测控制方法的预测模型，完成步骤S3中建立的约束下的航天器姿态控制任务。2.根据权利要求1所述的方法，其特征在于，所述步骤S1中的四元数描述方式下的航天器姿态运动学和动力学模型如下：3×3T3其中，J∈R为航天器惯量矩阵；ω＝(ω1,ω2,ω3)为航天器本体三轴角速度；H∈R3为航天器总角动量，其具体形式为H＝Jω+HRWA，HRWA∈R为反作用飞轮组合角动量在航天器3T惯量主轴上的投影；uc∈R为作用于航天器的姿态控制力矩；q＝(q0,q1,q2,q3)为描述航天器姿态的四元数，且满足Ξ(ω)是斜反对称矩阵，其具体形式如下：对于任意一个向量x∈R3，其对应的斜反对称阵x×形式如下：在忽略航天器轨道角速度对反作用飞轮角动量的影响时，反作用飞轮组合的动力学模型如下：LL其中，τRWA∈R为由L个反作用飞轮产生的控制力矩；hRWA∈R为反作用飞轮组合的角动L×L量，为反作用飞轮组合角动量关于时间的导数；JRWA∈R为反作用飞轮组合的转动惯量矩阵，其形式为为各个反作用飞轮的转动惯量；为各个反作用飞轮的2CN110502025A权利要求书2/6页转速，为各个反作用飞轮转速的关于时间的导数；基于反作用飞轮组合的动力学模型及反作用飞轮组合构型，得到反作用飞轮组合作用在航天器上的姿态控制力矩为：其中，C∈R3×L为反作用飞轮组合安装矩阵；根据角动量守恒定理，作用于航天器的姿态控制力矩uc为将式(5)代入式(1)，得到由反作用飞轮组合作为执行机构的航天器姿态运动学与动力学模型：将所述航天器姿态运动学与动力学模型离散化，取采样时间间隔为Δt，以将所述航天器姿态运动学与动力学模型作为模型预测控制方法的预测模型，离散化后的航天器姿态运动学与动力学模型为：其中，下标k和k+1分别表示k时刻和k+1时刻；Hk为k时刻的航天器总角动量；I4×4为四阶单位矩阵；上标x为斜反对称矩阵；T令xk＝(qk,ωk)为第k时刻航天器姿态控制系统状态量，表征航天器的姿态和角速度，为第k时刻系统控制量，将离散化后的航天器姿态运动学与动力学模型(7)简化为如下形式：定义则xk+1＝F(xk,uk)(10)。3.根据权利要求2所述的方法，其特征在于，步骤S2中建立的航天器体装式太阳能帆板输出功率与太阳入射角的关系为：其中，Ptotal为航天器上m个体装式太阳能帆板的输出总功率；为第i个太阳能帆板3CN110502025A权利要求书3/6页的输出功率；为第i个太阳能帆板的有效面积；θi为太阳入射角，表示第i个太阳能帆板法线与太阳光的夹角；PEOL(θi)为太阳入射角θi的函数，表示第i个太阳能帆板的单位面积功率，其具体关系式如下：其中，Ps为入射太阳辐射；η为太阳能帆板光电转换率；Id为太阳能帆板固有退化