发明名称 一种基于指令力矩矢量调节的控制力矩陀螺奇异规避方法
摘要 一种基于指令力矩矢量调节的控制力矩陀螺奇异规避方法,首先采集控制力矩陀螺群框架角向量,进而计算控制力矩陀螺运动方程的雅克比矩阵、奇异度量值、指令力矩矢量调节系数,然后根据姿态控制器给出的控制力矩指令得到控制力矩指令调节矢量、零运动奇异规避强度系数,最后根据控制力矩指令调节矢量、零运动奇异规避强度系数得到控制力矩陀螺框架角速度指令向量来控制力矩陀螺框架角速度。本发明克服了在框架“锁死”时力矩指令与其特定方向重合而无法脱离的情况,解决了在现有技术在奇异规避过程中所存在的框架角“锁死”问题,实现了对控制力矩陀螺奇异的有效规避。
申请公布号 CN105388902A 申请公布日期 2016.03.09
申请号 CN201510860445.1 申请日期 2015.11.30
申请人 北京控制工程研究所 发明人 雷拥军;姚宁;刘洁;赵江涛;朱琦;何海锋;李晶心;曹荣向
分类号 G05D1/08(2006.01)I;G05B13/02(2006.01)I 主分类号 G05D1/08(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 陈鹏
主权项 一种基于指令力矩矢量调节的控制力矩陀螺奇异规避方法,其特征在于包括如下步骤:(1)采集控制力矩陀螺群框架角向量δ,然后根据控制力矩陀螺运动方程计算框架角运动的雅克比矩阵J;(2)计算奇异度量值Sv为Sv=det(J·J<sup>T</sup>)进而得到指令力矩矢量调节系数υ为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&upsi;</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><mi>S</mi><mi>v</mi><mo>&gt;</mo><msub><mi>D</mi><mrow><mi>s</mi><mn>3</mn></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>k</mi><mrow><mi>t</mi><mi>m</mi><mi>p</mi><mi>&upsi;</mi></mrow></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>S</mi><mi>v</mi><mo>-</mo><msub><mi>D</mi><mrow><mi>s</mi><mn>3</mn></mrow></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><mi>S</mi><mi>v</mi><mo>&le;</mo><msub><mi>D</mi><mrow><mi>s</mi><mn>3</mn></mrow></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000863239770000011.GIF" wi="798" he="159" /></maths>当|υ|&gt;υ<sub>Limt</sub>时,υ=sgn(υ)·υ<sub>Limt</sub>,其中,k<sub>tmpυ</sub>为调节增益,D<sub>s3</sub>为矢量调节阈值,υ<sub>Limt</sub>为矢量调节系数限幅值,sgn(υ)为取υ的符号运算;(3)根据姿态控制器给出的控制力矩指令τ<sub>c</sub>得到控制力矩指令调节矢量τ<sub>cAd</sub>为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&tau;</mi><mrow><mi>c</mi><mi>A</mi><mi>d</mi></mrow></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mi>&upsi;</mi></mrow></mtd><mtd><mi>&upsi;</mi></mtd></mtr><mtr><mtd><mi>&upsi;</mi></mtd><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mi>&upsi;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>&upsi;</mi></mrow></mtd><mtd><mi>&upsi;</mi></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>&CenterDot;</mo><msub><mi>&tau;</mi><mi>c</mi></msub></mrow>]]></math><img file="FDA0000863239770000012.GIF" wi="602" he="263" /></maths>(4)根据奇异度量值Sv确定零运动奇异规避强度系数α<sub>s1</sub>、鲁棒奇异规避强度系数α<sub>s2</sub>,对于零运动奇异规避强度系数α<sub>s1</sub>,当Sv&lt;D<sub>s1</sub>时,零运动奇异规避强度系数<img file="FDA0000863239770000013.GIF" wi="526" he="182" />当D<sub>s1</sub>≤Sv&lt;D<sub>s2</sub>时,零运动奇异规避强度系数<img file="FDA0000863239770000014.GIF" wi="654" he="183" />当Sv≥D<sub>s2</sub>时,零运动奇异规避强度系数α<sub>s1</sub>=0;其中,D<sub>s1</sub>为鲁棒奇异规避启动阈值,D<sub>s2</sub>为零运动奇异规避启动阈值,C<sub>1</sub>为零运动奇异规避强度增益系数,α<sub>s10</sub>为偏置量;(5)对于鲁棒奇异规避强度系数α<sub>s2</sub>,当Sv&lt;D<sub>s1</sub>时,鲁棒奇异规避强度系数<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>=</mo><msub><mi>C</mi><mn>2</mn></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><mi>S</mi><mi>v</mi><mo>+</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>6</mn></mrow></msup></mrow></mfrac><mo>-</mo><mfrac><mn>1</mn><mrow><msub><mi>D</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>+</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>6</mn></mrow></msup></mrow></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000863239770000021.GIF" wi="887" he="189" /></maths>当Sv≥D<sub>s1</sub>时,鲁棒奇异规避强度系数α<sub>s2</sub>=0,其中,C<sub>2</sub>为鲁棒奇异规避强度增益系数;(6)根据控制力矩指令调节矢量τ<sub>cAd</sub>、奇异规避强度系数α<sub>s1</sub>及α<sub>s2</sub>得到控制力矩陀螺框架角速度指令向量<img file="FDA0000863239770000022.GIF" wi="63" he="78" />为<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mover><mi>&delta;</mi><mo>&CenterDot;</mo></mover><mi>d</mi></msub><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><msub><mi>H</mi><mrow><mi>c</mi><mi>m</mi><mi>g</mi><mn>0</mn></mrow></msub></mfrac><msup><mi>J</mi><mi>T</mi></msup><msup><mrow><mo>(</mo><mi>J</mi><mo>&CenterDot;</mo><msup><mi>J</mi><mi>T</mi></msup><mo>+</mo><msub><mi>&alpha;</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><msub><mi>I</mi><mn>3</mn></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&tau;</mi><mrow><mi>c</mi><mi>A</mi><mi>d</mi></mrow></msub><mo>+</mo><msub><mi>&alpha;</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><msub><mi>I</mi><mn>3</mn></msub><mo>-</mo><msup><mi>J</mi><mi>T</mi></msup><msup><mrow><mo>(</mo><mrow><mi>J</mi><mo>&CenterDot;</mo><msup><mi>J</mi><mi>T</mi></msup><mo>+</mo><msub><mi>&alpha;</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><msub><mi>I</mi><mn>3</mn></msub></mrow><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>J</mi><mo>)</mo></mrow><mfrac><mrow><mo>&part;</mo><mi>S</mi><mi>v</mi><mrow><mo>(</mo><mi>&delta;</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><mi>&delta;</mi></mrow></mfrac><mo>.</mo></mrow>]]></math><img file="FDA0000863239770000023.GIF" wi="1940" he="177" /></maths>使用控制力矩陀螺框架角速度指令向量<img file="FDA0000863239770000025.GIF" wi="68" he="78" />调整控制力矩陀螺群中控制力矩陀螺框架角速度,其中,H<sub>cmg0</sub>为控制力矩陀螺角动量,I<sub>3</sub>为3阶单位矩阵,<img file="FDA0000863239770000024.GIF" wi="109" he="143" />为Sv对框架角向量δ的偏导数,‑1为矩阵逆运算符号。
地址 100080 北京市海淀区北京2729信箱