发明名称 一种太阳敏感器故障情况下的对日定向姿态控制方法
摘要 一种太阳敏感器故障情况下的对日定向姿态控制方法,利用单轴姿态测量实现双轴的太阳定向稳定控制。当其中一个太阳敏感器发生故障失效时,可利用另外一个健康的太阳敏感器测量信息,结合陀螺角速度测量,利用观测方程重构故障轴上的姿态,从而实现对日定向时的稳定姿态控制。本发明方法具有在太阳敏感器单轴测量故障情况下依然能够实现对日定向的能力,在实现上不需要在卫星上额外增加新的姿态测量设备,只需要在星上控制计算机内通过软件实现本方法即可实现故障情况下的双轴姿态控制,具有实现经济简单的特点,可用作卫星对日定向的一种备份方式。
申请公布号 CN103466103A 申请公布日期 2013.12.25
申请号 CN201310372791.6 申请日期 2013.08.23
申请人 北京控制工程研究所 发明人 林波;曾海波;何英姿;武云丽;范炜;沈莎莎;刘潇翔;周中泽;杜建伟;吴小明;刘江
分类号 B64G1/36(2006.01)I;G01C21/24(2006.01)I 主分类号 B64G1/36(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 陈鹏
主权项 1.一种太阳敏感器故障情况下的对日定向姿态控制方法,其特征在于包括如下步骤:(1)在航天器的-Z面上布局两个太阳敏感器,其中一个太阳敏感器的测量轴与卫星本体的X轴平行,测量角为φ,另一个太阳敏感器的测量轴与卫星本体的Y轴平行,测量角为θ;(2)周期性的采集两个太阳敏感器的数据,计算m<sub>x,n</sub>=tan(φ<sub>n</sub>),m<sub>y,n</sub>=tan(θ<sub>n</sub>),其中φ<sub>n</sub>,θ<sub>n</sub>分别为本周期X轴和Y轴的太阳敏感器测量值,m<sub>x,n</sub>,m<sub>y,n</sub>分别为本周期X轴和Y轴的太阳敏感器测量角度的正切值;(3)建立观测方程,并利用陀螺本周期的测量值,利用观测方程预报下一周期X轴和Y轴的太阳敏感器测量角度的正切估计值,<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>x</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>&omega;</mi><mi>y</mi></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>&omega;</mi><mi>x</mi></msub><msubsup><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&Delta;t</mi></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>y</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>&omega;</mi><mi>x</mi></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>&omega;</mi><mi>y</mi></msub><msubsup><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&Delta;t</mi></mrow></math>]]></maths>其中<img file="FDA0000371084530000013.GIF" wi="124" he="76" /><img file="FDA0000371084530000014.GIF" wi="91" he="77" />分别为本周期X轴和Y轴的太阳敏感器测量角度的正切估计值,初值为0,<img file="FDA0000371084530000015.GIF" wi="154" he="76" /><img file="FDA0000371084530000016.GIF" wi="126" he="77" />为下一周期X轴和Y轴的太阳敏感器测量角度的正切估计值,Δt为控制周期,ω<sub>x</sub>、ω<sub>y</sub>、ω<sub>z</sub>为本周期陀螺测量的卫星本体三轴角速度;(4)判断太阳敏感器的健康状态,利用健康太阳敏感器的测量值对步骤(3)计算得到的<img file="FDA0000371084530000017.GIF" wi="152" he="76" /><img file="FDA0000371084530000018.GIF" wi="126" he="77" />进行修正,当X轴太阳敏感器健康,Y轴太阳敏感器故障时,对X轴太阳敏感器测量角度的正切估计值进行修正,修正方程为<img file="FDA0000371084530000019.GIF" wi="629" he="76" /><maths num="0003"><![CDATA[<math><mrow><msub><mi>K</mi><mi>x</mi></msub><mo>&lt;</mo><mn>0</mn><mo>,</mo><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mi>K</mi><mi>y</mi></msub><mo>+</mo><msubsup><mi>&omega;</mi><mi>z</mi><mn>2</mn></msubsup><mo>></mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>当X轴太阳敏感器故障,Y轴太阳敏感器健康时,对Y轴太阳敏感器测量角度的正切估计值进行修正,修正方程为<img file="FDA00003710845300000111.GIF" wi="637" he="77" /><maths num="0004"><![CDATA[<math><mrow><msub><mi>K</mi><mi>y</mi></msub><mo>&lt;</mo><mn>0</mn><mo>,</mo><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mi>K</mi><mi>x</mi></msub><mo>-</mo><msubsup><mi>&omega;</mi><mi>z</mi><mn>2</mn></msubsup><mo>&lt;</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>其中K<sub>x</sub>,K<sub>y</sub>为修正系数;(5)更新本周期X轴和Y轴的太阳敏感器测量角度的正切估计值,<img file="FDA0000371084530000021.GIF" wi="477" he="76" />由此计算本周期X轴和Y轴的太阳敏感器测量估计值<img file="FDA0000371084530000022.GIF" wi="590" he="88" />其中<img file="FDA0000371084530000023.GIF" wi="48" he="81" />为X轴太阳敏感器测量估计值,<img file="FDA0000371084530000024.GIF" wi="49" he="70" />为Y轴太阳敏感器测量估计值;(6)依据太阳敏感器健康状态和故障太阳敏感器所在轴的测量估计值重构姿态,作为控制的输入量;当X轴的太阳敏感器健康而Y轴的太阳敏感器故障时,令<img file="FDA0000371084530000025.GIF" wi="266" he="82" />当X轴的太阳敏感器故障而Y轴的太阳敏感器健康时,令<img file="FDA0000371084530000026.GIF" wi="269" he="82" />当两轴的太阳敏感器均健康时,令<img file="FDA0000371084530000027.GIF" wi="269" he="71" />当两轴的太阳敏感器均故障时,令φ=0,θ=0;(7)利用步骤(6)得到的计算结果φ和θ,计算对日定向所需的姿态控制量,控制律为:<img file="FDA0000371084530000028.GIF" wi="424" he="75" />U<sub>y</sub>=K<sub>py</sub>θ+K<sub>dy</sub>ω<sub>y</sub>,U<sub>z</sub>=K<sub>dz</sub>(ω<sub>z</sub>-ω<sub>zb</sub>),式中U<sub>m</sub>,m=x,y,z为卫星的三轴控制量,K<sub>dm</sub>为角速度控制系数,ω<sub>m</sub>为卫星陀螺实测三轴角速度,ω<sub>zb</sub>为卫星Z轴目标角速度,K<sub>px</sub>,K<sub>py</sub>分别为滚动角和俯仰角控制系数。
地址 100080 北京市海淀区北京2729信箱