发明名称 一种航空遥感惯性稳定平台耦合力矩估计与补偿方法
摘要 一种航空遥感惯性稳定平台耦合力矩估计与补偿方法,通过建立三轴惯性稳定平台动力学模型,根据POS、横滚/俯仰陀螺、方位码盘实时测出的信息,估算出基座及框架间耦合力矩,然后计算出抵消此耦合力矩所需的力矩电机电流值,并将此电流值补偿到电流环的电流给定输入值中,使电机输出与耦合力矩大小相等、方向相反的力矩,补偿平台耦合力矩干扰,提高平台稳定精度。本发明实现了耦合力矩的实时估算与补偿,提高了稳定精度,适用于具有基座及框架间耦合力矩的航空遥感惯性稳定平台。
申请公布号 CN103175530B 申请公布日期 2015.04.08
申请号 CN201310068184.0 申请日期 2013.03.04
申请人 北京航空航天大学 发明人 周向阳;张宏燕;李建平;房建成
分类号 G01C21/18(2006.01)I;G01C25/00(2006.01)I 主分类号 G01C21/18(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 成金玉;卢纪
主权项 一种航空遥感惯性稳定平台耦合力矩估计与补偿方法,其特征在于:基于建立的三轴动力学方程,根据POS、陀螺、码盘测量值估算耦合力矩值,进一步计算补偿电流值,具体包括以下步骤:(1)针对航空遥感三轴惯性稳定平台建立动基座下简化动力学模型,如下:①方位框绕方位轴的动力学方程为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>J</mi><mi>az</mi></msub><mrow><mo>(</mo><msubsup><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>ifz</mi><mi>f</mi></msubsup><mo>+</mo><msub><mover><mi>&theta;</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>M</mi><mi>z</mi></msub></mrow>]]></math><img file="FDA0000662428730000011.GIF" wi="376" he="89" /></maths>②俯仰框组件绕俯仰轴的动力学方程简化为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mrow><mo>[</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>ax</mi></msub><mo>+</mo><msub><mi>J</mi><mi>fx</mi></msub><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mi>ifx</mi><mi>f</mi></msubsup><mo>]</mo></mrow><mo>&prime;</mo></msup><mo>+</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>az</mi></msub><mo>+</mo><msub><mi>J</mi><mi>fz</mi></msub><mo>-</mo><msub><mi>J</mi><mi>ay</mi></msub><mo>-</mo><msub><mi>J</mi><mi>fy</mi></msub><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mi>ify</mi><mi>f</mi></msubsup><msubsup><mi>&omega;</mi><mi>ibz</mi><mi>b</mi></msubsup><mo>-</mo><msub><mi>J</mi><mi>az</mi></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>a</mi></msub><msubsup><mi>&omega;</mi><mi>ify</mi><mi>f</mi></msubsup><mo>=</mo><msub><mi>M</mi><mi>x</mi></msub></mrow>]]></math><img file="FDA0000662428730000012.GIF" wi="1249" he="101" /></maths>③横滚框组件绕横滚轴方向动力学方程简化为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mo>[</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>ry</mi></msub><mo>+</mo><msub><mi>J</mi><mi>ay</mi></msub><mo>+</mo><msub><mi>J</mi><mi>fy</mi></msub><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mi>ify</mi><mi>f</mi></msubsup><msup><mo>]</mo><mo>&prime;</mo></msup><mo>+</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>rx</mi></msub><mo>-</mo><msub><mi>J</mi><mi>rz</mi></msub><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mi>ibx</mi><mi>b</mi></msubsup><msubsup><mi>&omega;</mi><mi>ibz</mi><mi>b</mi></msubsup><mo>+</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>ax</mi></msub><mo>+</mo><msub><mi>J</mi><mi>fx</mi></msub><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mi>ifx</mi><mi>f</mi></msubsup><msubsup><mi>&omega;</mi><mi>ibz</mi><mi>b</mi></msubsup></mtd></mtr><mtr><mtd><mo>-</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>az</mi></msub><mo>+</mo><msub><mi>J</mi><mi>fz</mi></msub><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mi>ibz</mi><mi>b</mi></msubsup><msubsup><mi>&omega;</mi><mi>ibz</mi><mi>b</mi></msubsup><mo>+</mo><msub><mi>J</mi><mi>az</mi></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>a</mi></msub><msubsup><mi>&omega;</mi><mi>ibx</mi><mi>b</mi></msubsup><mo>=</mo><msub><mi>M</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000662428730000013.GIF" wi="1178" he="195" /></maths>其中定义:J<sub>ax</sub>、J<sub>ay</sub>、J<sub>az</sub>分别为方位框架在x、y、z轴的转动惯量;J<sub>fx</sub>、J<sub>fy</sub>、J<sub>fz</sub>分别为俯仰框架在x、y、z轴的转动惯量;J<sub>rx</sub>、J<sub>ry</sub>、J<sub>rz</sub>分别为横滚框架在x、y、z轴的转动惯量;<img file="FDA0000662428730000014.GIF" wi="73" he="84" />为基座相对惯性空间的角速度在基座坐标系的表达,其在x、y、z轴上的投影分别记为<img file="FDA0000662428730000015.GIF" wi="481" he="92" />为俯仰框架相对惯性空间的角速度在俯仰框坐标系的表达,其在x、y、z轴上的投影分别记为<img file="FDA0000662428730000016.GIF" wi="232" he="92" /><img file="FDA0000662428730000017.GIF" wi="216" he="92" />为俯仰框架相对惯性空间的角加速度在俯仰框坐标系的表达在z轴的投影;θ<sub>a</sub>为方位框架相对俯仰框架转角,<img file="FDA0000662428730000021.GIF" wi="58" he="86" />为方位框架相对俯仰框架转角角速度;<img file="FDA0000662428730000022.GIF" wi="53" he="87" />为方位框架相对俯仰框架转角角加速度;M<sub>x</sub>为俯仰框组件沿x轴受到的综合力矩,M<sub>y</sub>为横滚框组件沿y轴受到的综合力矩,M<sub>z</sub>为方位框架沿z轴受到的综合力矩;(2)利用安装于基座的POS测量<img file="FDA0000662428730000023.GIF" wi="96" he="83" />和<img file="FDA0000662428730000024.GIF" wi="131" he="83" />利用安装于俯仰框架的横滚和俯仰陀螺测量<img file="FDA0000662428730000025.GIF" wi="83" he="86" />和<img file="FDA0000662428730000026.GIF" wi="129" he="86" />利用安装于方位框架的码盘测量θ<sub>a</sub>,间接计算<img file="FDA0000662428730000027.GIF" wi="102" he="85" />(3)根据步骤(2)中直接测量值或间接计算值<img file="FDA0000662428730000028.GIF" wi="533" he="91" />θ<sub>a</sub>和<img file="FDA0000662428730000029.GIF" wi="56" he="84" />计算横滚框耦合力矩值M<sub>couple_y</sub>和俯仰框耦合力矩值M<sub>couple_x</sub>;(4)根据步骤(3)中计算得到的横滚框耦合力矩值M<sub>couple_y</sub>和俯仰框耦合力矩值M<sub>couple_x</sub>,进一步计算出抵消横滚框耦合力矩M<sub>couple_y</sub>所需的力矩电机电流值I<sub>com_y</sub>,抵消俯仰框耦合力矩M<sub>couple_x</sub>所需的力矩电机电流值I<sub>com_x</sub>,并将计算出的电流值补偿到电流环的电流给定输入值中,电机输出与耦合力矩大小相等、方向相反的力矩,补偿惯性平台耦合干扰力矩;所述步骤(3)中根据步骤(2)中直接测量值或间接计算值计算俯仰框耦合力矩值M<sub>couple_x</sub>、横滚框耦合力矩值M<sub>couple_y</sub>,具体如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mi>couple</mi><mo>_</mo><mi>x</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>az</mi></msub><mo>+</mo><msub><mi>J</mi><mi>fz</mi></msub><mo>-</mo><msub><mi>J</mi><mi>ay</mi></msub><mo>-</mo><msub><mi>J</mi><mi>fy</mi></msub><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mi>ify</mi><mi>f</mi></msubsup><msubsup><mi>&omega;</mi><mi>ibz</mi><mi>b</mi></msubsup><mo>-</mo><msub><mi>J</mi><mi>az</mi></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>a</mi></msub><msubsup><mi>&omega;</mi><mi>ify</mi><mi>f</mi></msubsup></mrow>]]></math><img file="FDA00006624287300000210.GIF" wi="994" he="93" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mi>couple</mi><mo>_</mo><mi>y</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>rx</mi></msub><mo>-</mo><msub><mi>J</mi><mi>rz</mi></msub><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mi>ibx</mi><mi>b</mi></msubsup><msubsup><mi>&omega;</mi><mi>ibz</mi><mi>b</mi></msubsup><mo>+</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>ax</mi></msub><mo>+</mo><msub><mi>J</mi><mi>fx</mi></msub><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mi>ifx</mi><mi>f</mi></msubsup><msubsup><mi>&omega;</mi><mi>ibz</mi><mi>b</mi></msubsup><mo>-</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>az</mi></msub><mo>+</mo><msub><mi>J</mi><mi>fz</mi></msub><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mi>ibz</mi><mi>b</mi></msubsup><msubsup><mi>&omega;</mi><mi>ibz</mi><mi>b</mi></msubsup><mo>+</mo><msub><mi>J</mi><mi>az</mi></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>a</mi></msub><msubsup><mi>&omega;</mi><mi>ibx</mi><mi>b</mi></msubsup><mo>.</mo></mrow>]]></math><img file="FDA00006624287300000211.GIF" wi="1588" he="87" /></maths>
地址 100191 北京市海淀区学院路37号