发明名称 一种用于控制力矩陀螺奇异规避的航天器姿态控制方法
摘要 本发明公开了一种用于控制力矩陀螺奇异规避的航天器姿态控制方法,其具体步骤为:第一步,SGCMG群角动量集合的求取;第二步,SGCMG群与航天器整体系统角动量幅值判定;第三步,确定航天器姿控系统最优控制性能指标;第四步,确定SGCMG群与航天器整体系统线性化模型;第五步,确定SGCMG群与航天器所组成整体系统的线性化范围;第六步,确定满足约束条件的线性化范围;第七步,利用非线性预测控制实现航天器姿态控制。本发明完成对航天器姿态的精确控制,达到了降低航天器发射成本、提高航天器在轨运行寿命的目的。
申请公布号 CN105223961A 申请公布日期 2016.01.06
申请号 CN201510666229.3 申请日期 2015.10.16
申请人 北京机械设备研究所 发明人 张佳为;许诺;张丽华
分类号 G05D1/08(2006.01)I;G05D1/10(2006.01)I 主分类号 G05D1/08(2006.01)I
代理机构 中国航天科工集团公司专利中心 11024 代理人 岳洁菱;姜中英
主权项 一种用于控制力矩陀螺奇异规避的航天器姿态控制方法,其特征在于具体步骤为:第一步  确定SGCMG群角动量集合依据SGCMG群的奇异角动量<img file="617613dest_path_image002.GIF" wi="20" he="25" />绘制SGCMG群冗余构型奇异动量面;SGCMG群奇异角动量表达式为:<img file="dest_path_image004.GIF" wi="206" he="57" />(1)其中<img file="dest_path_image006.GIF" wi="14" he="16" />为SGCMG个数,<img file="dest_path_image008.GIF" wi="17" he="25" />为第<img file="dest_path_image010.GIF" wi="10" he="18" />个SGCMG的角动量,<img file="dest_path_image012.GIF" wi="21" he="25" />为SGCMG群奇异方向单位矢量,<img file="dest_path_image014.GIF" wi="16" he="25" />为第<img file="891907dest_path_image010.GIF" wi="10" he="18" />个SGCMG框架轴方向单位矢量;依据所绘制的奇异动量面得到SGCMG群角动量集合包含最大的以原点为球心的球域半径<img file="dest_path_image016.GIF" wi="13" he="14" />为:<img file="dest_path_image018.GIF" wi="154" he="32" />(2)其中<img file="dest_path_image020.GIF" wi="13" he="17" />为SGCMG群框架角,<img file="dest_path_image022.GIF" wi="17" he="16" />表示由坐标原点指向动量包络面的方向矢量;<img file="dest_path_image024.GIF" wi="85" he="32" />表示SGCMG群动量包络面上的最小角动量,即求取各个方向矢量<img file="978679dest_path_image022.GIF" wi="17" he="16" />上的最大角动量,然后在所求得的各个方向上的最大角动量中求得最小值;第二步  判定SGCMG群与航天器组成的整体系统角动量幅值当SGCMG群与航天器组成的整体系统初始角动量幅值<img file="dest_path_image026.GIF" wi="25" he="25" />满足<img file="dest_path_image028.GIF" wi="60" he="28" />时,确定航天器姿控系统最优控制性能指标;第三步  确定航天器姿控系统最优控制性能指标由传感器实时采集航天器姿态四元数<img file="dest_path_image030.GIF" wi="13" he="17" />、角速度<img file="dest_path_image032.GIF" wi="16" he="14" />、SGCMG群框架角<img file="41270dest_path_image020.GIF" wi="13" he="17" />及框架角速度<img file="dest_path_image034.GIF" wi="14" he="28" />作为状态变量,反馈给航天器姿控系统,确定航天器姿控系统初始时刻<img file="dest_path_image036.GIF" wi="16" he="25" />、控制时域<img file="dest_path_image038.GIF" wi="21" he="20" />、预测时域<img file="dest_path_image040.GIF" wi="16" he="18" />及正定矩阵<img file="dest_path_image042.GIF" wi="22" he="18" />,正定矩阵<img file="dest_path_image044.GIF" wi="21" he="18" />及正定矩阵<img file="dest_path_image046.GIF" wi="18" he="18" />,建立SGCMG群与航天器所组成整体系统的最优控制性能指标:<img file="dest_path_image048.GIF" wi="428" he="38" />(3)其中<img file="dest_path_image050.GIF" wi="137" he="49" />,<img file="dest_path_image052.GIF" wi="113" he="49" />,输入<img file="dest_path_image054.GIF" wi="41" he="30" />,<img file="dest_path_image056.GIF" wi="22" he="26" />为系统平衡点;<img file="dest_path_image058.GIF" wi="89" he="26" />,<img file="dest_path_image060.GIF" wi="70" he="25" />,<img file="dest_path_image062.GIF" wi="52" he="25" />,其中<img file="dest_path_image064.GIF" wi="33" he="25" />为<img file="dest_path_image066.GIF" wi="36" he="20" />维单位矩阵,<img file="dest_path_image068.GIF" wi="38" he="25" />为<img file="dest_path_image070.GIF" wi="45" he="20" />维单位矩阵,<img file="dest_path_image072.GIF" wi="22" he="25" />为<img file="705029dest_path_image006.GIF" wi="14" he="16" />维单位矩阵;第四步  确定SGCMG群与航天器整体系统线性化模型SGCMG群与航天器所组成整体系统在平衡点邻域的线性化模型为:<img file="dest_path_image074.GIF" wi="429" he="113" />(4)其中<img file="dest_path_image076.GIF" wi="17" he="18" />为航天器指令姿态四元数对应的欧拉角,<img file="dest_path_image078.GIF" wi="18" he="25" />为SGCMG群的框架角平衡点,<img file="dest_path_image080.GIF" wi="104" he="76" />,<img file="dest_path_image082.GIF" wi="33" he="25" />为<img file="184331dest_path_image006.GIF" wi="14" he="16" />维单位矩阵,<img file="dest_path_image084.GIF" wi="16" he="20" />为航天器转动惯量,<img file="dest_path_image086.GIF" wi="36" he="28" />表示三维向量对应的反对称矩阵,<img file="dest_path_image088.GIF" wi="21" he="18" />为SGCMG角动量,<img file="dest_path_image090.GIF" wi="14" he="20" />为飞轮角动量幅值,<img file="dest_path_image092.GIF" wi="18" he="18" />为SGCMG群的安装矩阵,<img file="dest_path_image094.GIF" wi="20" he="25" />为系统外干扰力矩;第五步 确定SGCMG群与航天器所组成整体系统的线性化范围依据SGCMG群与航天器所组成整体系统的线性化模型确定状态反馈控制器<img file="dest_path_image096.GIF" wi="18" he="18" />,依据<img file="991970dest_path_image094.GIF" wi="20" he="25" />最大幅值<img file="dest_path_image098.GIF" wi="14" he="20" />、控制输入约束<img file="dest_path_image100.GIF" wi="66" he="25" />确定系统线性化范围;状态变量系数矩阵为<img file="dest_path_image102.GIF" wi="20" he="25" />,干扰力矩系数矩阵为<img file="dest_path_image104.GIF" wi="20" he="25" />,控制输入系数矩阵为<img file="dest_path_image106.GIF" wi="22" he="25" />,依据不等式约束公式(5)~公式(7),得到<img file="dest_path_image108.GIF" wi="14" he="18" />的下确界及正定矩阵<img file="dest_path_image110.GIF" wi="18" he="18" />:<img file="dest_path_image112.GIF" wi="49" he="22" /><img file="dest_path_image114.GIF" wi="222" he="52" />(5)<img file="dest_path_image116.GIF" wi="500" he="48" />(6)<img file="dest_path_image118.GIF" wi="238" he="57" />(7)其中<img file="dest_path_image120.GIF" wi="18" he="25" />为<img file="dest_path_image122.GIF" wi="10" he="20" />个临近系统平衡点的已知状态点,<img file="dest_path_image124.GIF" wi="36" he="22" />为待求正定矩阵,<img file="dest_path_image126.GIF" wi="38" he="22" />,<img file="dest_path_image128.GIF" wi="18" he="25" />为<img file="dest_path_image130.GIF" wi="18" he="20" />的第<img file="244089dest_path_image010.GIF" wi="10" he="18" />行;<img file="dest_path_image132.GIF" wi="432" he="78" />(8)其中<img file="dest_path_image134.GIF" wi="58" he="28" />为<img file="dest_path_image136.GIF" wi="37" he="21" />的第<img file="114699dest_path_image010.GIF" wi="10" he="18" />行;第六步 确定满足约束条件的线性化范围满足系统状态约束及控制输入约束的线性化范围:<img file="dest_path_image138.GIF" wi="236" he="30" />(9)其中<img file="dest_path_image140.GIF" wi="146" he="28" />,<img file="dest_path_image142.GIF" wi="102" he="28" />;<img file="dest_path_image144.GIF" wi="329" he="196" />,<img file="dest_path_image146.GIF" wi="93" he="78" />,<img file="dest_path_image148.GIF" wi="26" he="59" />为<img file="930297dest_path_image110.GIF" wi="18" he="64" />的第<img file="672119dest_path_image010.GIF" wi="10" he="55" />行第<img file="dest_path_image150.GIF" wi="14" he="90" />列元素,<img file="dest_path_image152.GIF" wi="16" he="103" />为系统状态约束最大幅值;第七步 利用非线性预测控制实现航天器姿态控制在航天器机动准备阶段,采用序列二次规划离线确定性能指标<img file="dest_path_image154.GIF" wi="16" he="20" />对应的最优控制输入<img file="dest_path_image156.GIF" wi="20" he="25" />,并将最优控制输入<img file="945188dest_path_image156.GIF" wi="20" he="25" />作为下次序列二次规划实时迭代的初值,采用非线性预测控制进行循环运算,并实时判定系统当前状态变量是否满足<img file="dest_path_image158.GIF" wi="69" he="25" />,满足<img file="989629dest_path_image158.GIF" wi="69" he="25" />时将系统控制器切换为线性控制器<img file="dest_path_image160.GIF" wi="17" he="16" />,直至系统状态被控制到平衡点。
地址 100854 北京市海淀区北京142信箱208分箱