发明名称 一种深空探测器动力下降初始关键参数优化方法
摘要 本发明提供一种深空探测器动力下降初始关键参数优化方法,具体过程为:获取目标近月点高度h<sub>p</sub>与计算的近月点高度h之间的差量Δh,获取目标近月点月理纬度φ<sub>p</sub>与计算的月理纬度φ之间的差量Δφ;并以Δh和Δφ小于设定阈值为目标,通过迭代优化获取深空探测器动力下降初始关键参数,所述参数包括开始点火时刻T<sub>ib</sub>,点火速度增量ΔV<sub>f</sub>和点火方向<img file="DDA0000509615660000011.GIF" wi="256" he="88" />利用该方法能够得到一次环月降轨的参数,使深空探测器到达动力下降初始点时,月理纬度、高度、速度、测控弧段等关键参数满足任务要求。
申请公布号 CN104035333B 申请公布日期 2015.06.10
申请号 CN201410220160.7 申请日期 2014.05.23
申请人 北京空间飞行器总体设计部 发明人 董捷;周文艳;马继楠;杨眉;吴学英;张熇;孙泽洲
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 李爱英;杨志兵
主权项 一种深空探测器动力下降初始关键参数优化方法,其特征在于,具体步骤为:获取目标近月点高度h<sub>p</sub>与计算的近月点高度h之间的差量Δh,获取目标近月点月理纬度<img file="FDA0000685000980000011.GIF" wi="64" he="69" />与计算的月理纬度<img file="FDA0000685000980000012.GIF" wi="50" he="65" />之间的差量<img file="FDA0000685000980000013.GIF" wi="104" he="77" />并以Δh和<img file="FDA0000685000980000014.GIF" wi="85" he="76" />小于设定阈值为目标,通过迭代优化获取深空探测器动力下降初始关键参数,所述参数包括开始点火时刻T<sub>ib</sub>,点火速度增量ΔV<sub>f</sub>和点火方向<img file="FDA0000685000980000015.GIF" wi="256" he="109" />具体过程为:步骤一、设定参量:用于环月轨道外推的时刻T<sub>0</sub>,在T<sub>0</sub>时刻月心J2000惯性坐标系下的瞬时轨道参数<img file="FDA0000685000980000016.GIF" wi="260" he="98" />环月降轨脉冲点火时刻T<sub>i</sub>距时刻T<sub>0</sub>的时间间隔ΔT<sub>0</sub>,环月降轨点火速度增量<img file="FDA0000685000980000017.GIF" wi="656" he="90" />轨道飞行周期数N<sub>T</sub>;步骤二、从时刻T<sub>0</sub>开始,在月心J2000惯性坐标系下根据所述<img file="FDA0000685000980000018.GIF" wi="226" he="87" />将轨道外推至T<sub>i</sub>时刻,计算T<sub>i</sub>时刻对应的瞬时轨道参数<img file="FDA0000685000980000019.GIF" wi="258" he="101" />将所述<img file="FDA00006850009800000110.GIF" wi="229" he="101" />外推N<sub>T</sub>周期后计算对应近月点高度h和月理纬度<img file="FDA00006850009800000111.GIF" wi="70" he="66" />步骤三、判断|Δh|&lt;0.1km、<img file="FDA00006850009800000112.GIF" wi="260" he="85" />且深空探测器的累计时间达到预设量是否都成立,若是则将当前的ΔT<sub>0</sub>和ΔVx<sub>0</sub>作为脉冲点火环月降轨参数,进入步骤五,若否则进入步骤四,其中Δh=h‑h<sub>p</sub>,<img file="FDA00006850009800000113.GIF" wi="236" he="81" />h<sub>p</sub>为目标近月点高度,<img file="FDA00006850009800000114.GIF" wi="60" he="69" />为目标近月点月理纬度;步骤四、令N<sub>T</sub>值加一,以h<sub>p</sub>、<img file="FDA00006850009800000115.GIF" wi="61" he="70" />为目标量,并采用微分校正法进行迭代计算,计算ΔT<sub>0</sub>的迭代量Δt,点火速度增量<img file="FDA00006850009800000116.GIF" wi="118" he="90" />的迭代量<img file="FDA00006850009800000117.GIF" wi="110" he="76" />然后更新ΔT<sub>0</sub>=ΔT<sub>0</sub>+Δt,更新<img file="FDA00006850009800000118.GIF" wi="359" he="84" />返回步骤二;步骤五、设置点火时间长度初值为ΔT<sub>f0</sub>;步骤六、根据步骤三得到的ΔT<sub>0</sub>确定脉冲点火时刻T<sub>i</sub>,并将其作为有限推力点火时间区间的中间时刻;从时刻T<sub>0</sub>开始,在月心J2000惯性坐标系下根据所述<img file="FDA0000685000980000021.GIF" wi="232" he="97" />将轨道外推至时刻T<sub>ib</sub>=T<sub>i</sub>‑1/2ΔT<sub>f0</sub>,计算至该时刻对应的瞬时轨道参数<img file="FDA00006850009800000210.GIF" wi="290" he="90" />基于所述<img file="FDA0000685000980000022.GIF" wi="253" he="87" />计算点火过程的轨道参数至时刻T<sub>ie</sub>=T<sub>i</sub>+1/2ΔT<sub>f0</sub>得到对应的瞬时轨道参数<img file="FDA0000685000980000023.GIF" wi="260" he="87" />其中点火过程中点火方向为根据所述步骤三计算的ΔVx<sub>0</sub>所确定的方向,并将所述<img file="FDA0000685000980000024.GIF" wi="235" he="91" />外推N<sub>T</sub>周期后计算对应近月点高度h和月理纬度<img file="FDA0000685000980000025.GIF" wi="69" he="67" />步骤七、判断|Δh|&lt;0.1km、<img file="FDA0000685000980000026.GIF" wi="260" he="85" />且深空探测器的累计时间达到预设量是否都成立,若是则将当前的T<sub>ib</sub>和ΔT<sub>f0</sub>作为有限推力点火环月降轨参数,进入步骤九,若否则进入步骤八,其中Δh=h‑h<sub>p</sub>,<img file="FDA0000685000980000027.GIF" wi="237" he="84" />步骤八、令N<sub>T</sub>值加一,以h<sub>p</sub>、<img file="FDA0000685000980000028.GIF" wi="62" he="67" />为目标量,采用微分校正法进行迭代计算,计算点火时间ΔT<sub>f0</sub>的迭代量Δt<sub>f</sub>,更新ΔT<sub>f0</sub>=ΔT<sub>f0</sub>+Δt<sub>f</sub>,进而更新T<sub>ib</sub>=T<sub>i</sub>‑1/2ΔT<sub>f0</sub>,返回步骤六;步骤九、根据计算的ΔT<sub>f0</sub>计算点火速度增量ΔV<sub>f</sub>,获得最终的环月降轨参数包括:T<sub>ib</sub>、ΔV<sub>f</sub>和<img file="FDA0000685000980000029.GIF" wi="259" he="106" />
地址 100094 北京市海淀区友谊路104号
您可能感兴趣的专利