主权项 |
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>ρ</mi><mi>j</mi></msub><mo>=</mo><msqrt><mn>2</mn><mo>-</mo><mn>2</mn><mi>cos</mi><msub><mi>θ</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>γ</mi><mi>d</mi><mi>k</mi></msubsup><mrow><msubsup><mi>γ</mi><mi>d</mi><mi>k</mi></msubsup><mo>+</mo><mi>β</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>β</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>ρ</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>ω</mi><mo>×</mo></msup><mo>]</mo><mi>Jω</mi><mo>+</mo><mi>J</mi><mrow><mo>(</mo><msub><mover><mi>ω</mi><mo>·</mo></mover><mi>s</mi></msub><mo>+</mo><mfrac><mi>μ</mi><mi>η</mi></mfrac><msub><mi>ω</mi><mi>s</mi></msub><mo>+</mo><mi>s</mi><mrow><mo>(</mo><msub><mi>ω</mi><mi>s</mi></msub><mo>-</mo><mi>ω</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>后并未达到完整的目标姿态,需再绕敏感器矢量方向旋转一定角度,从而实现航天器完整的姿态机动过程。 |