发明名称 火星着陆器喷气推力器和质量矩复合控制系统
摘要 本发明公开了一种火星着陆器喷气推力器和质量矩复合控制系统,该系统包括有姿态控制系统、复合执行机构、着陆器动力学和运动学模型;其中,姿态控制系统包括有姿态控制器和控制分配;其中,复合执行机构包括有RCS和MM。在考虑多种环境干扰因素情况下,建立基于RCS/MM的着陆器复合动力学模型;根据着陆器系统姿态误差产生控制系统所需要的控制力矩;控制分配则将总的控制力矩单独分配到两个执行机构,产生控制指令;RCS/MM系统根据各自的输入指令产生实际的控制力矩,调整着陆器姿态。本发明可有效降低着陆器系统燃料消耗,具有较强的机动性能,能够产生连续的控制力矩,改善着陆器的姿态控制精度和控制余度,为火星的精确着陆提供保障。
申请公布号 CN103869823B 申请公布日期 2016.06.29
申请号 CN201410138869.2 申请日期 2014.04.08
申请人 北京航空航天大学 发明人 吴忠;王振;郭雷
分类号 G05D1/10(2006.01)I 主分类号 G05D1/10(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 李有浩
主权项 一种火星着陆器喷气推力器和质量矩复合控制系统,火星着陆器至少包括有位置姿态控制器、火星着陆器动力学和运动学模型(3)和执行机构;其特征在于:所述执行机构是指喷气反作用控制系统RCS和变质心MM的复合执行机构(2);所述变质心MM中包括有第一质量滑块p和第二质量滑块q;第一质量滑块p安装在着陆器的本体体坐标系O<sub>b</sub>‑X<sub>b</sub>Y<sub>b</sub>Z<sub>b</sub>的Y<sub>b</sub>轴上,且第一质量滑块p是沿Y<sub>b</sub>轴往复运动;第二质量滑块q安装在着陆器的本体体坐标系O<sub>b</sub>‑X<sub>b</sub>Y<sub>b</sub>Z<sub>b</sub>的Z<sub>b</sub>轴上,且第二质量滑块q是沿Z<sub>b</sub>轴往复运动;所述火星着陆器动力学和运动学模型(3)是在喷气反作用控制系统RCS和变质心MM的复合执行机构(2)下构建的;所述位置姿态控制器中包括有姿态控制模块(11)和控制分配模块(12)构成的姿态控制系统(1);姿态控制系统(1)中的姿态控制模块(11)根据接收到的姿态误差<img file="FDA0000962891590000011.GIF" wi="53" he="79" />产生姿态控制所需的期望姿态控制力矩<img file="FDA0000962891590000012.GIF" wi="85" he="63" />并将所述<img file="FDA0000962891590000013.GIF" wi="50" he="61" />输出给控制分配模块(12)中的控制力矩分配模块(121);控制力矩分配模块(121)将所述的期望姿态控制力矩<img file="FDA0000962891590000014.GIF" wi="51" he="62" />进行分配处理后一方面输出喷气反作用控制系统RCS力矩指令<img file="FDA0000962891590000015.GIF" wi="517" he="95" />给喷气反作用控制系统RCS指令分解模块(122),另一方面输出变质心MM力矩指令<img file="FDA0000962891590000016.GIF" wi="564" he="94" />给变质心MM指令分解模块(123);τ<sub>rdx</sub>为控制力矩分配模块(121)分配给喷气反作用控制系统RCS的X轴的姿态控制力矩;τ<sub>rdy</sub>为控制力矩分配模块(121)分配给喷气反作用控制系统RCS的Y轴的姿态控制力矩;τ<sub>rdz</sub>为控制力矩分配模块(121)分配给喷气反作用控制系统RCS的Z轴的姿态控制力矩;τ<sub>mdx</sub>为控制力矩分配模块(121)分配给变质心MM的X轴的姿态控制力矩;τ<sub>mdy</sub>为控制力矩分配模块(121)分配给变质心MM的Y轴的姿态控制力矩;τ<sub>mdz</sub>为控制力矩分配模块(121)分配给变质心MM的Z轴的姿态控制力矩;喷气反作用控制系统RCS指令分解模块(122)对接收到的所述<img file="FDA0000962891590000021.GIF" wi="516" he="95" />进行分解处理,输出喷气反作用控制系统RCS启动指令PWM<sub>RCS</sub>给喷气反作用控制系统RCS和变质心MM的复合执行机构(2)中的喷气反作用控制系统RCS;变质心MM指令分解模块(123)对接收到的所述<img file="FDA0000962891590000022.GIF" wi="564" he="95" />进行分解处理,输出变质心MM的位置指令s<sub>p</sub>、s<sub>q</sub>给喷气反作用控制系统RCS和变质心MM的复合执行机构(2)中的变质心MM;s<sub>p</sub>表示第一质量滑块p的位置指令,s<sub>q</sub>表示第二质量滑块q的位置指令;喷气反作用控制系统RCS和变质心MM的复合执行机构(2)中的喷气反作用控制系统RCS依据PWM<sub>RCS</sub>指令产生控制力矩<img file="FDA0000962891590000023.GIF" wi="82" he="62" />变质心MM依据s<sub>p</sub>、s<sub>q</sub>指令产生控制力矩<img file="FDA0000962891590000024.GIF" wi="91" he="63" />喷气反作用控制系统RCS和变质心MM的复合执行机构(2)的力矩总和<img file="FDA0000962891590000025.GIF" wi="43" he="46" />为作用于火星着陆器动力学和运动学模型(3)的复合控制力矩;火星着陆器动力学和运动学模型(3)是在喷气反作用控制系统RCS和变质心MM的复合执行机构(2)的力矩总和<img file="FDA0000962891590000026.GIF" wi="38" he="48" />与环境因素产生的空气干扰力矩<img file="FDA0000962891590000027.GIF" wi="90" he="71" />共同作用下,改变火星着陆器姿态,进而改变火星着陆器的着陆轨迹,最终完成火星着陆器姿态和位置的解算;所述火星着陆器动力学和运动学模型(3)中火星着陆器平动时的动力学模型为<img file="FDA0000962891590000028.GIF" wi="1278" he="167" />M为火星着陆器的总质量;<img file="FDA0000962891590000029.GIF" wi="102" he="70" />为火星着陆器所受的空气动力;<img file="FDA00009628915900000210.GIF" wi="109" he="71" />为火星引力势能;<img file="FDA00009628915900000211.GIF" wi="108" he="71" />为喷气反作用控制系统RCS的推力矢量;<img file="FDA00009628915900000212.GIF" wi="158" he="167" />为在惯性坐标系O<sub>a</sub>‑X<sub>a</sub>Y<sub>a</sub>Z<sub>a</sub>下的火星着陆器质心点的二次微分矢量;<img file="FDA00009628915900000213.GIF" wi="45" he="63" />为火星着陆器质心O<sub>b</sub>点的绝对矢量;d表示微分标识;t为火星着陆器的着陆时刻;a为惯性坐标系的标识;<img file="FDA0000962891590000031.GIF" wi="149" he="79" />为建模时动态引入的火星着陆器燃耗、RCS推力或者MM运动的作用力不确定项;<img file="FDA0000962891590000032.GIF" wi="142" he="87" />为摄动加速度;<img file="FDA0000962891590000033.GIF" wi="886" he="183" />其中:μ<sub>s</sub>为太阳引力常数;μ<sub>e</sub>为地球引力常数;<img file="FDA0000962891590000034.GIF" wi="54" he="62" />为太阳相对于着陆器的位置矢量;<img file="FDA0000962891590000035.GIF" wi="54" he="62" />为地球相对于着陆器的位置矢量;<img file="FDA0000962891590000036.GIF" wi="60" he="69" />为太阳在火星惯性系的位置矢量;<img file="FDA0000962891590000037.GIF" wi="66" he="66" />为地球在火星惯性系的位置矢量;<img file="FDA0000962891590000038.GIF" wi="134" he="71" />为质量滑块相对着陆器的运动产生的惯性力;<img file="FDA0000962891590000039.GIF" wi="1259" he="419" />其中:p为第一质量滑块的标识;q为第二质量滑块的标识;<img file="FDA00009628915900000310.GIF" wi="46" he="55" />为第一质量滑块p的位置矢量;<img file="FDA00009628915900000311.GIF" wi="44" he="55" />为第二质量滑块q的位置矢量;m<sub>p</sub>为第一质量滑块p的质量;m<sub>q</sub>为第二质量滑块q的质量;<img file="FDA00009628915900000312.GIF" wi="141" he="175" />为在火星着陆器的本体体坐标系O<sub>b</sub>‑X<sub>b</sub>Y<sub>b</sub>Z<sub>b</sub>下的第一质量滑块的二次微分位置矢量;<img file="FDA00009628915900000313.GIF" wi="133" he="174" />为在火星着陆器的本体体坐标系O<sub>b</sub>‑X<sub>b</sub>Y<sub>b</sub>Z<sub>b</sub>下的第二质量滑块的二次微分位置矢量;<img file="FDA00009628915900000314.GIF" wi="102" he="158" />为在火星着陆器的本体体坐标系O<sub>b</sub>‑X<sub>b</sub>Y<sub>b</sub>Z<sub>b</sub>下的第一质量滑块的一次微分位置矢量;<img file="FDA00009628915900000315.GIF" wi="109" he="159" />为在火星着陆器的本体体坐标系O<sub>b</sub>‑X<sub>b</sub>Y<sub>b</sub>Z<sub>b</sub>下的第二质量滑块的一次微分位置矢量;<img file="FDA00009628915900000316.GIF" wi="78" he="70" />为火星着陆器本体坐标系相对于惯性坐标系的角加速度矢量的反对称矩阵;<img file="FDA0000962891590000041.GIF" wi="78" he="79" />为火星着陆器本体坐标系相对于惯性坐标系的角加速度矢量的反对称矩阵的一阶微分,下角标ab为两个坐标系的标识,上角标“×”表示反对称矩阵形式,<img file="FDA0000962891590000042.GIF" wi="713" he="271" /><img file="FDA0000962891590000043.GIF" wi="93" he="63" />为火星着陆器本体坐标系相对于惯性坐标系的X轴角加速度矢量,<img file="FDA0000962891590000044.GIF" wi="94" he="70" />为火星着陆器本体坐标系相对于惯性坐标系的Y轴角加速度矢量,<img file="FDA0000962891590000045.GIF" wi="93" he="67" />为火星着陆器本体坐标系相对于惯性坐标系的Z轴角加速度矢量;所述火星着陆器动力学和运动学模型(3)中火星着陆器转动时的动力学模型为<img file="FDA0000962891590000046.GIF" wi="1475" he="103" />J<sub>O</sub>为火星着陆器的转动惯量;<img file="FDA0000962891590000047.GIF" wi="78" he="63" />为火星着陆器相对于惯性坐标系的角速度矢量;<img file="FDA0000962891590000048.GIF" wi="82" he="79" />为角速度矢量<img file="FDA0000962891590000049.GIF" wi="77" he="63" />的一阶微分;<img file="FDA00009628915900000410.GIF" wi="85" he="71" />为空气干扰力矩;<img file="FDA00009628915900000411.GIF" wi="97" he="63" />为建模时动态引入的喷气反作用控制系统RCS控制力矩、空气动力矩、空气阻尼力矩或者质量滑块运动的力矩不确定项;<img file="FDA00009628915900000412.GIF" wi="78" he="71" />为变质心MM运动产生的惯性力矩;且:<img file="FDA00009628915900000413.GIF" wi="1206" he="819" />μ<sub>p</sub>为第一质量滑块的质量比,<img file="FDA00009628915900000414.GIF" wi="230" he="127" />μ<sub>q</sub>为第二质量滑块的质量比,<img file="FDA00009628915900000416.GIF" wi="230" he="127" />ΔJ为变质心MM运动引起的附加转动惯量;且:<img file="FDA0000962891590000051.GIF" wi="837" he="207" />式中:<img file="FDA0000962891590000052.GIF" wi="63" he="63" />为第一质量块p位置矢量的反对称矩阵,上角标“×”表示反对称矩阵形式,<img file="FDA0000962891590000053.GIF" wi="477" he="269" />s<sub>p</sub>表示第一质量滑块p的位置指令;<img file="FDA0000962891590000054.GIF" wi="59" he="63" />为第二质量块q位置矢量的反对称矩阵,上角标“×”表示反对称矩阵形式,<img file="FDA0000962891590000055.GIF" wi="469" he="271" />s<sub>q</sub>表示第二质量滑块q的位置指令。
地址 100191 北京市海淀区学院路37号