发明名称 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法
摘要 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法,本发明涉及异面交叉快变轨道稳定姿态指向跟踪控制方法。本发明为了解决现有技术未考虑航天器的惯量不确定性,依赖于状态初值,无法自由调整收敛时间,以及飞轮在奇异方向产生的补偿力矩需要人为设计的问题。具体是按照以下步骤进行的:步骤一、设追踪星与目标星位于异面交叉轨道上,需要确定期望姿态;步骤二、期望姿态跟踪控制律的设计;步骤三、消除期望姿态跟踪控制律的抖振;步骤四、追踪星与目标星轨道交叉点的期望姿态随追踪星与目标星轨道交叉点间的距离而变化,根据期望姿态跟踪控制律来确定执行机构的配置方案,求解期望姿态控制力矩。本发明应用于卫星控制领域。
申请公布号 CN104527994B 申请公布日期 2016.08.24
申请号 CN201510030643.5 申请日期 2015.01.21
申请人 哈尔滨工业大学 发明人 孙延超;马广富;刘昱晗;赵文锐;朱津津;李传江;李卓
分类号 B64G1/24(2006.01)I;B64G3/00(2006.01)I 主分类号 B64G1/24(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法,其特征在于:异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法具体是按照以下步骤进行的:步骤一、设追踪星与目标星位于异面交叉轨道上,追踪星需要自主探测目标星的位置,需要确定期望姿态;具体过程为:期望姿态坐标系各坐标轴单位矢量在惯性系中由下式确定:<img file="FDA0000906078420000011.GIF" wi="238" he="128" /><img file="FDA0000906078420000012.GIF" wi="921" he="139" />k=i×j其中r<sub>t</sub>与r<sub>c</sub>分别为追踪星与目标星在地心惯性坐标系中的位置矢量;若期望姿态坐标系与质心轨道坐标系相对质心惯性坐标系的余弦转换矩阵分别为Rdi和Roi,则期望姿态坐标系相对于质心轨道坐标系的转换矩阵为<img file="FDA0000906078420000013.GIF" wi="931" he="75" />步骤二、期望姿态跟踪控制律的设计;具体步骤为:在期望姿态坐标系各坐标轴单位矢量下对期望姿态角和本体姿态角进行计算:将期望姿态角和本体姿态角分别记为a<sub>d</sub>和a<sub>b</sub>,它们的差为欧拉姿态角误差,记为e,即e=a<sub>b</sub>‑a<sub>d</sub>,令系统状态为<img file="FDA0000906078420000014.GIF" wi="227" he="74" />可写出线性化状态方程<img file="FDA0000906078420000015.GIF" wi="1162" he="217" />其中I为对应轴主惯量,在姿态跟踪过程中期望姿态角加速度是有界的,将<img file="FDA0000906078420000016.GIF" wi="118" he="67" />项当成有界干扰γ,设<img file="2.GIF" wi="203" he="63" />γ<sub>0</sub>为干扰上界,<img file="FDA0000906078420000018.GIF" wi="42" he="57" />为欧拉姿态角的一阶导数,<img file="FDA0000906078420000019.GIF" wi="40" he="57" />为欧拉姿态角的二阶导数,u为输入力矩,<img file="FDA00009060784200000110.GIF" wi="59" he="67" />为本体姿态角的二阶导数;线性化状态方程<img file="FDA00009060784200000111.GIF" wi="647" he="217" />满足算法假设条件rank[b,Ab]=n,式中,n=2为系统阶数,<img file="FDA00009060784200000112.GIF" wi="277" he="155" /><img file="FDA00009060784200000113.GIF" wi="183" he="214" />u为输入力矩;并对线性化状态方程<img file="FDA0000906078420000021.GIF" wi="646" he="216" />进行线性变换,x=Ge,G=[Ab,b]<sup>‑1</sup>,可得到<img file="FDA0000906078420000022.GIF" wi="1166" he="158" />x<sub>1</sub>和x<sub>2</sub>是x=Ge中的x,是个矢量<img file="FDA0000906078420000023.GIF" wi="191" he="155" /><img file="FDA0000906078420000024.GIF" wi="54" he="67" />是x<sub>1</sub>的一阶导数,<img file="FDA0000906078420000025.GIF" wi="55" he="67" />是x<sub>2</sub>的二阶导数,γ为有界干扰,u为输入力矩;引入不丢失正负号信息的乘方运算符z<sup>[q]</sup>=|z|<sup>q</sup>sign(z),z,q∈R;滑模面的形式为:<img file="FDA0000906078420000026.GIF" wi="1155" he="78" />式中,α<sub>1</sub>β<sub>1</sub>是可变的常值系数;期望姿态跟踪控制律中期望姿态跟踪滑模控制律具有如下形式:<img file="FDA0000906078420000027.GIF" wi="1327" he="135" />式中,u<sub>1</sub>为期望力矩,<img file="FDA0000906078420000028.GIF" wi="810" he="76" />sign(s)为符号函数;α<sub>1</sub>、α<sub>2</sub>、β<sub>1</sub>和β<sub>2</sub>都是可变常值系数,是可以任意赋值的变量;<img file="FDA0000906078420000029.GIF" wi="641" he="75" />由期望姿态跟踪滑模控制律(6)形式可知需要调节的参数为T<sub>max</sub>和γ<sub>0</sub>;x在时间T<sub>max</sub>内收敛到原点,即欧拉姿态角误差e在有限时间内收敛到原点;步骤三、消除期望姿态跟踪控制律的抖振;具体步骤为:消除期望姿态跟踪控制律中期望姿态跟踪滑模控制律的抖振,在姿态跟踪滑模控制律中用饱和函数sat(s)代替符号函数sign(s),<img file="FDA00009060784200000210.GIF" wi="1248" he="135" /><img file="FDA00009060784200000211.GIF" wi="450" he="291" />其中Δ&gt;0,为边界层;<img file="FDA00009060784200000212.GIF" wi="640" he="94" />步骤四、追踪星与目标星轨道交叉点的期望姿态随追踪星与目标星轨道交叉点间的距离而变化,根据期望姿态跟踪控制律来确定执行机构的配置方案,求解期望姿态控制力矩;具体步骤为:选用反作用飞轮和控制力矩陀螺组合配置;当偏航轴用两个平行放置的单框架控制力矩陀螺控制;整个卫星角动量为零,即在卫星本体角速度为零时,使h<sub>1</sub>和h<sub>2</sub>方向相反,框架角为零;当两个陀螺转子角动量大小都为h,框架角分别为δ<sub>1</sub>和δ<sub>2</sub>,陀螺群在卫星本体系中总角动量为<img file="FDA0000906078420000031.GIF" wi="1147" he="158" />h<sub>1</sub>和h<sub>2</sub>为二维角动量,有方向有大小,只能输出x、z方向力矩;h<sub>1</sub>和h<sub>2</sub>是矢量,h是在某个方向上的h<sub>1</sub>和h<sub>2</sub>的标量;各陀螺框架转动产生的合成陀螺力矩T可表示为:<img file="FDA0000906078420000032.GIF" wi="1408" he="166" />式中,δ为陀螺框架角,<img file="FDA0000906078420000033.GIF" wi="50" he="83" />是陀螺的框架角速度,C(δ)为陀螺群的力矩矩阵,<img file="FDA0000906078420000034.GIF" wi="67" he="67" />是H的一阶导数,<img file="FDA0000906078420000035.GIF" wi="57" he="81" />是δ<sub>1</sub>的一阶导数,<img file="FDA0000906078420000036.GIF" wi="58" he="81" />是δ<sub>2</sub>的一阶导数;控制力矩陀螺产生奇异,会陷入奇异状态,不能有效输出控制力矩,设计合适的操纵律来回避或脱离奇异状态,奇异值度量为:D=det(CC<sup>T</sup>)    (10)奇异时D=0,非奇异时D&gt;0,且该值越大表明奇异程度越小,框架角速度采用鲁棒伪逆操纵律进行计算,具有如下形式:<img file="FDA0000906078420000037.GIF" wi="1023" he="127" />其中,<img file="FDA0000906078420000038.GIF" wi="40" he="67" />为框架角速度,T为指令力矩,只需输出z轴力矩,即[0,T<sub>z</sub>]<sup>T</sup>;α为权系数,可根据D的大小实时调整;E是单位矩阵;C为陀螺群的力矩矩阵;h为陀螺转子角动量大小;而滚动轴和俯仰轴分别用一个飞轮控制,期望姿态控制力矩工作模式下飞轮作为单纯一阶惯性环节,对<img file="FDA0000906078420000039.GIF" wi="1056" he="139" />进行分析与控制器设计;式中,t<sub>s</sub>为一阶惯性系统的时间常数;采用PD控制器进行设计,从e到飞轮实际输出期望姿态控制力矩u<sub>w</sub>的传递函数为<img file="FDA0000906078420000041.GIF" wi="1131" he="139" />式中,k<sub>p</sub>为比例环节系数,k<sub>d</sub>为微分环节系数,<img file="FDA0000906078420000042.GIF" wi="87" he="72" />为系统的传递函数;期望姿态控制力矩u<sub>w</sub>为<img file="1.GIF" wi="99" he="63" />。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号