发明名称 一种航天器指向约束姿态机动控制方法
摘要 本发明涉及一种受指向约束的自主姿态机动控制的方法,属于航天器姿态控制技术领域。该方法采用在单位球面S上,以敏感器目标指向矢量r<sub>d</sub>的末端点为目标点位置,敏感器当前的指向矢量r的末端点为当前位置,指向约束所形成的球冠面为障碍区域,构建关于r的末端点在球面上运动的导航函数V,然后根据导航函设计控制力矩表达式,通过改变控制力矩参数来调整控制力矩的幅值,从而驱动航天器使得敏感器指向目标矢量r<sub>d</sub>,航天器完成敏感器指向目标矢量r<sub>d</sub>后,需再绕敏感器矢量方向旋转角度θ,从而实现航天器完整的姿态机动过程。本发明方法能够明确地对障碍区域的指向规避进行处理,同时对于多个障碍约束时能避免局部极小值,确保航天器安全机动到目标姿态,满足了执行机构有界性的要求,实现了自主航天器指向约束姿态机动的控制。
申请公布号 CN102331785B 申请公布日期 2014.04.16
申请号 CN201110199698.0 申请日期 2011.07.16
申请人 北京理工大学 发明人 崔祜涛;程小军;崔平远;朱圣英;徐瑞
分类号 G05D1/08(2006.01)I 主分类号 G05D1/08(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 李爱英;郭德忠
主权项 1.一种航天器指向约束姿态机动控制的方法,其特征在于:该方法实现的具体过程如下:步骤一、根据航天器本身的结构安装信息,得到指向约束的敏感器方向矢量在本体系下的表示为r<sub>b</sub>,在惯性系下的表示为r;根据姿态敏感器信息和星历信息,得到航天器质心到相关的n个天体的矢量在惯性系下的表示分别为r<sub>oj</sub>,j=1,...,n;根据姿态敏感器信息,得到航天器在惯性系下相对于在本体系下的姿态矩阵C<sub>Ib</sub>,下标表示该姿态矩阵是从本体系b到惯性系I的转换;根据期望姿态矩阵<img file="FDA0000388661810000011.GIF" wi="112" he="87" />得到惯性系下敏感器指向目标矢量r<sub>d</sub>;步骤二、以航天器质心为中心,建立单位球面S;根据敏感器的视场顶角θ<sub>0</sub>以及相关的n个天体的视角θ<sub>oj</sub>,得到敏感器方向矢量r与第n个天体矢量r<sub>oj</sub>之间的约束关系为r<sup>T</sup>r<sub>oj</sub><cosθ<sub>j</sub>,其中θ<sub>j</sub>=θ<sub>0</sub>+θ<sub>oj</sub>,即敏感器方向矢量r在姿态机动过程中不能进入到由第j个天体矢量r<sub>oj</sub>为对称中心,顶点在航天器质心,锥顶角为θ<sub>j</sub>的空间锥内,n个天体的指向将形成n个空间锥,这些空间锥与单位球面S相交后,将单位球面S截出n个球冠面,其中,第j个球冠面的中心与球冠边缘的距离为ρ<sub>j</sub>,<maths num="0001"><![CDATA[<math><mrow><msub><mi>&rho;</mi><mi>j</mi></msub><mo>=</mo><msqrt><mn>2</mn><mo>-</mo><mn>2</mn><mi>cos</mi><msub><mi>&theta;</mi><mi>j</mi></msub></msqrt><mo>;</mo></mrow></math>]]></maths>步骤三、在单位球面S上,以敏感器指向目标矢量r<sub>d</sub>的末端点为目标点位置,敏感器当前的方向矢量r的末端点为当前位置,指向约束所形成的球冠面为障碍区域,构建关于r的末端点在球面上运动的导航函数V;导航函数V为:<maths num="0002"><![CDATA[<math><mrow><msup><mrow><mi>V</mi><mo>=</mo><mrow><mo>(</mo><mfrac><msubsup><mi>&gamma;</mi><mi>d</mi><mi>k</mi></msubsup><mrow><msubsup><mi>&gamma;</mi><mi>d</mi><mi>k</mi></msubsup><mo>+</mo><mi>&beta;</mi></mrow></mfrac><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mi>k</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></math>]]></maths>其中,k为常数,k≥2,γ<sub>d</sub>=||r-r<sub>d</sub>||<sup>2</sup>,β和β<sub>j</sub>为中间变量,<img file="FDA0000388661810000014.GIF" wi="272" he="155" /><maths num="0003"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mi>j</mi></msub><mo>=</mo><msup><mrow><mo>|</mo><mo>|</mo><mi>r</mi><mo>-</mo><msub><mi>r</mi><mi>oj</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>&rho;</mi><mi>j</mi><mn>2</mn></msubsup><mo>;</mo></mrow></math>]]></maths>步骤四、将导航函数融合到控制律的设计过程中,结合姿态动力学和运动学,并利用反演法backstepping设计控制力矩u表达式,航天器在控制力矩驱动下旋转;控制力矩u为:<maths num="0004"><![CDATA[<math><mrow><mi>u</mi><mo>=</mo><mo>[</mo><msup><mi>&omega;</mi><mo>&times;</mo></msup><mo>]</mo><mi>J&omega;</mi><mo>+</mo><mi>J</mi><mrow><mo>(</mo><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>s</mi></msub><mo>+</mo><mfrac><mi>&mu;</mi><mi>&eta;</mi></mfrac><msub><mi>&omega;</mi><mi>s</mi></msub><mo>+</mo><mi>s</mi><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>s</mi></msub><mo>-</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,μ,η,s为调节参数,J为航天器的惯量矩阵,ω为航天器的姿态角速度,ω<sub>s</sub>为期望姿态角速度和<img file="FDA0000388661810000022.GIF" wi="70" he="77" />为期望姿态角速度的导数,[ω<sup>×</sup>]为航天器姿态角速度的叉乘矩阵;步骤五、根据步骤四中确定的控制力矩表达式,通过改变控制参数μ,η和s,来调整控制力矩的幅值,从而满足航天器执行机构的输出要求;步骤六、根据步骤四和步骤五所确定的控制力矩来驱动航天器,从而使敏感器指向目标矢量r<sub>d</sub>;步骤七、航天器完成敏感器指向目标矢量r<sub>d</sub>后并未达到完整的目标姿态,需再绕敏感器矢量方向旋转一定角度,从而实现航天器完整的姿态机动过程。
地址 100081 北京市海淀区中关村南大街5号