本发明公开的一种基于路径点规划的火星大气进入段预测制导方法，涉及一种火星大气进入段预测制导方法，属于深空探测技术领域。本发明包括如下步骤：步骤1：建立探测器动力学模型；步骤2：引入路径点规划方法，通过规划路径点(e，sinγ)，使预测制导生成的进入轨迹能够满足过载约束条件；步骤3：求取制导指令σcmd；步骤4：每个制导周期重复步骤2、3，直至满足开伞条件。本发明首次将路径点规划引入预测制导中，以保证进入过程中飞行器始终满足路径约束，从而保证飞行器及乘员安全。本发明可用于火星大气进入段预测制导。
1.一种基于路径点规划的火星大气进入段预测制导方法，包括步骤1，步骤1、建立探测器动力学模型；具体方法为，考虑火星自转影响的探测器对无量纲时间 的三自由度无量纲进入动力学模型为，
其中，s为剩余纵程，表征从探测器当前位置到目标开伞位置的火星表面大圆弧的距离，r为火星质心到探测器质心的距离，无量纲参数为火星半径R0，v为探测器相对于火星的速度，无量纲参数为 其中g0为火星表面重力加速度，γ为航迹角，σ为倾侧角，g为当地重力加速度，其无量纲参数为g0；D和L分别为阻力加速度和升力加速度：其无量纲参数均为g0，CD和CL分别为阻力系数和升力系数，S为探测器参考面积，m为探测器质量，q＝ρv2/2为动压，β＝m/SCD为探测器弹道系数，L/D为探测器升阻比；火星大气密度ρ采用指数模型如公式(3)，其中ρ0为参考密度，h0为参考高度，hs为大气密度标高；
其特征在于：具体实现方法还包括步骤2、3、4，步骤2：引入路径点规划方法，所述的路径点规划方法是指通过规划路径点(e，sinγ)，使预测制导生成的进入轨迹能够满足过载约束条件，进而保证探测器及乘员安全；
在确定sinγ-e剖面(航迹角正弦-能量剖面)后，制导指令σcmd由动力学方程(1)反解得到，其中，参考航迹角的导数 由下述关系式(5)确定，求得σcmd即为要求取的制导指令；
步骤4：每个制导周期重复步骤2、3，直至满足开伞条件。
2.根据权利要求1所述的一种基于路径点规划的火星大气进入段预测制导方法，其特征在于：所述的步骤2具体实现方法如下，进入段所需要满足的过载约束如下：
其中，gE为地球表面重力加速度，式(6)可以表达为下述形式式(7)，其中， C1＝CnDρ0/2β；
步骤2.1：根据路径约束，确定路径点(e，sinγ)所需满足的物理条件；
为保证在进入段中飞行器满足路径约束，则在规划的路径点处需满足，或
其中，(·)tra|(r,v)和(·)lf|(r,v)分别为在r-v空间内，飞行轨迹以及过载边界上在(r,v)点处的导数；根据进入动力学方程(1)，有，根据式(7)，在过载边界上有，
式(14)给出了轨迹与过载约束边界曲线曲率相同时的数学关系，但要使过载约束得到满足，还需保证在曲率相同时，轨迹不能位于过载边界下方；
在利用式(14)进行路径点规划时，必须首先确定比能量e与航迹角正弦sinγ之间的关系，所述的比能量e与航迹角正弦sinγ之间的关系通过步骤2.2加以确定；
步骤2.2：建立路径点(e，sinγ)所需满足的比能量e与航迹角正弦sinγ的数学关系；
所述的比能量e与航迹角正弦sinγ的数学关系根据工程需要确定，可以为公式(15)或二次曲线关系；
其中，(e0,sinγ0)和(ef,sinγf)为sinγ-e空间(航迹角正弦-能量空间)内初始和末端点，(e1,sinγ1)和(e2,sinγ2)为路径规划算法需要确定的剖面控制点；为了保证在转折点处可导，在转折点处采用哈密尔顿插值算法；
由于在进入过程中，需要满足末端开伞位置约束式(16)s(ef)＝0 (16)
因此，sinγ-e剖面(航迹角正弦-能量剖面)需有单一可变量；为此，令其中△γ为常数，为路径规划算法需要确定的变量；
步骤2.3：联立式(14)与(15)，即可求得路径点(e，sinγ)。
一种基于路径点规划的火星大气进入段预测制导方法
技术领域
[0001] 本发明涉及火星大气进入段预测制导方法，尤其涉及一种基于路径点规划的火星大气进入段预测制导方法，属于深空探测技术领域。
背景技术
[0002] 在火星着陆探测任务中，不仅要满足着陆的精确性，还要保证着陆安全性。在火星进入、下降与着陆过程中，火星大气进入段是历时最久，探测器热防护系统工作条件最恶劣的一个阶段，保障其安全实施尤为重要。该阶段中，由于受到结构强度以及热防护系统的制约，进入段制导律的设计需要在保证着陆精度的同时，满足过载、热流等路径约束，从而保证任务的安全实施。
[0003] 火星大气进入段现有制导方法主要分为标称轨迹法和预测制导法。其中，标称轨迹法可以通过设计标称轨迹来满足路径约束，进而保证进入段的安全性。预测制导法根据所预测的开伞位置偏差产生制导指令，从而对偏差加以校正。该方法能够实现较高的开伞位置精度，并且对初始再入条件不敏感，在火星进入制导方法研究领域受到广泛关注。现有火星大气进入段预测制导的相关文献对于保证进入段末端位置精度这一问题进行了较为充分的研究，而对满足路径约束这一问题需要做进一步研究。
[0004] 为了保证火星大气进入段飞行的安全性，有必要针对预测制导法满足路径约束这一问题，设计一种能够同时保证着陆精度和相关约束的制导律。
发明内容
[0005] 本发明公开的一种基于路径点规划的火星大气进入段预测制导方法，要解决的技术问题是使火星大气进入段预测制导方法在保障着陆精度的同时，能够在进入过程中满足路径约束，从而保证探测器及乘员安全。
[0006] 本发明是通过下述技术方案实现的：
[0007] 本发明公开的一种基于路径点规划的火星大气进入段预测制导方法，在已有火星大气进入段预测制导方法的基础上，引入路径点规划方法，从而使得该制导方法能够同时满足着陆精度与路径约束，从而保证探测器及乘员安全。
[0008] 本发明公开的一种基于路径点规划的火星大气进入段预测制导方法，包括如下步骤：
[0009] 步骤1、建立探测器动力学模型。具体方法为，
[0010] 考虑火星自转影响的探测器对无量纲时间 的三自由度无量纲进入动力学模型为，
[0011]
[0012] 其中，s为剩余纵程，表征从探测器当前位置到目标开伞位置的火星表面大圆弧的距离，r为火星质心到探测器质心的距离，无量纲参数为火星半径R0，v探测器相对于火星的速度，无量纲参数为 其中g0为火星表面重力加速度，γ为航迹角，σ为倾侧角，g为当地重力加速度，其无量纲参数为g0。D和L分别阻力加速度和升力加速度：
[0013]
[0014] 其无量纲参数均为g0，CD和CL分别为阻力系数和升力系数，S为探测器参考面积，m为探测器质量，q＝ρv2/2为动压，β＝m/SCD为探测器弹道系数，L/D为探测器升阻比。火星大气密度ρ采用指数模型如公式(3)，
[0015]
[0016] 其中ρ0为参考密度，h0为参考高度，hs为大气密度标高。
[0017] 步骤2：引入路径点规划方法，所述的路径点规划方法是指通过规划路径点(e，sinγ)，使预测制导生成的进入轨迹能够满足过载约束条件，进入段所需要满足的过载约束如下：
[0018]
[0019] 其中，gE为地球表面重力加速度，式(4)可以表达为下述形式式(5)，[0020]
[0021] 其中， C1＝CnDρ0/2β。
[0022] 步骤2.1：根据路径约束，确定路径点(e，sinγ)所需满足的物理条件。
[0023] 为保证在进入段中飞行器满足路径约束，则在规划的路径点处需满足，[0024]
[0025] 或
[0026]
[0027] 其中，(·)tra|(r,v)和(·)lf|(r,v)分别为在r-v空间内，飞行轨迹以及过载边界上在(r,v)点处的导数。根据进入动力学式(1)，有
[0028]
[0029] 根据式(5)，在过载边界上有
[0030]
[0031] 则根据式(7)，有
[0032]
[0033] 定义进入段飞行器的比能量为，
[0034]
[0035] 则路径点可根据式(12)进行规划，
[0036]
[0037] 式(12)给出了轨迹与过载约束边界曲线曲率相同时的数学关系，但要使过载约束得到满足，还需保证在曲率相同时，轨迹不能位于过载边界下方。
[0038] 在利用式(12)进行路径点规划时，必须首先确定能量e与航迹角正弦sinγ之间的关系，所述的能量e与航迹角正弦sinγ之间的关系通过步骤2.2加以确定。
[0039] 步骤2.2：建立路径点(e，sinγ)所需满足的能量e与航迹角正弦sinγ的数学关系。
[0040] 所述的能量e与航迹角正弦sinγ的数学关系根据工程需要确定，可以为公式(13)或二次曲线关系。
[0041]
[0042] 其中，(e0,sinγ0)和(ef,sinγf)为sinγ-e空间(航迹角正弦-能量空间)内初始和末端点，(e1,sinγ1)和(e2,sinγ2)为路径规划算法需要确定的剖面控制点。为了保证在转折点处可导，在转折点处采用哈密尔顿插值算法。
[0043] 由于在进入过程中，需要满足末端开伞位置约束式(14)
[0044] s(ef)＝0 (14)
[0045] 因此，sinγ-e剖面(航迹角正弦-能量剖面)需有单一可变量。为此，令[0046]
[0047] 其中△γ为常数， 为路径规划算法需要确定的变量。
[0048] 步骤2.3：联立式(12)与(13)，即可求得路径点(e，sinγ)。
[0049] 步骤3：求取制导指令σcmd。
[0050] 在确定sinγ-e剖面(航迹角正弦-能量剖面)后，制导指令σcmd可由动力学方程(1)反解得到，
[0051]
[0052] 其中，参考航迹角的导数 由下述关系式(17)确定，
[0053]
[0054] 求得σcmd即为要求取的制导指令。
[0055] 步骤4：每个制导周期重复步骤2、3，直至满足开伞条件。
[0056] 有益效果：
[0057] 本发明公开的一种基于路径点规划的火星大气进入段预测制导方法，首次将路径点规划引入预测制导中，以保证进入过程中飞行器始终满足路径约束，从而保证飞行器及乘员安全。
附图说明
[0058] 图1为基于路径点规划的预测制导指令生成流程图。
[0059] 图2为火星大气进入段r-v空间(火心距-速度空间)内路径约束条件下的可行域。
[0060] 图3给出了采用路径点规划的预测制导方法在sinγ-e空间(航迹角正弦-能量空间)上的轨迹。
具体实施方式
[0061] 为了更好的说明本发明的目的和优点，下面结合附图和实例对发明内容做进一步说明。
[0062] 实施例1：
[0063] 本实施例公开的一种基于路径点规划的火星大气进入段预测制导方法，包括如下步骤：
[0064] 步骤1：建立探测器动力学模型。具体方法为，
[0065] 考虑火星自转影响的探测器对无量纲时间 的三自由度无量纲进入动力学模型为：
[0066]
[0067] 其中，s为剩余纵程，表征从探测器当前位置到目标开伞位置的火星表面大圆弧的距离，r为火星质心到探测器质心的距离，无量纲参数为火星半径R0，v探测器相对于火星的速度，无量纲参数为 其中g0为火星表面重力加速度，γ为航迹角，σ为倾侧角，g为当地重力加速度，其无量纲参数为g0。D和L分别阻力加速度和升力加速度：
[0068]
[0069] 其无量纲参数均为g0，CD和CL分别为阻力系数和升力系数，S为探测器参考面积，m为探测器质量，q＝ρv2/2为动压，β＝m/SCD为探测器弹道系数，L/D为探测器升阻比。火星大气密度ρ采用指数模型如公式(20)，
[0070]
[0071] 其中ρ0＝2×10-4kg/m3为参考密度，h0＝40000m为参考高度，hs＝7500km为大气密度标高。
[0072] 步骤2：引入路径点规划方法，所述的路径点规划方法是指通过规划路径点(e，sinγ)使预测制导生成的进入轨迹能够满足过载约束条件，进入段所需要满足的过载约束如下：
[0073]
[0074] 其中，gE＝9.81m/s2为地球表面重力加速度，(21)式可以表达为下述形式，[0075]
[0076] 其中， C1＝CnDρ0/2β。
[0077] 步骤2.1：根据路径约束，确定路径点(e，sinγ)所需满足的物理条件。
[0078] 为保证在进入段中飞行器满足路径约束，则在规划的路径点处需满足，[0079]
[0080] 或
[0081]
[0082] 其中，(·)tra|(r,v)和(·)lf|(r,v)分别为在r-v空间(火心距-速度空间)内，飞行轨迹以及过载边界上在(r,v)点处的导数。根据进入动力学式(18)，有
[0083]
[0084] 根据式(22)，在过载边界上有
[0085]
[0086] 则根据式(24)，有
[0087]
[0088] 定义进入段飞行器的比能量
[0089]
[0090] 则路径点可根据下式进行规划
[0091]
[0092] 其中，1/r在路径点处可以近似看做常量，取为0.991。式(29)给出了轨迹与过载约束边界曲线曲率相同时的数学关系，但要使过载约束得到满足，还需保证在曲率相同时，轨迹不能位于过载边界下方。
[0093] 在利用式(29)进行路径点规划时，必须首先确定能量e与航迹角正弦sinγ之间的关系，该关系通过步骤2.2加以确定。
[0094] 步骤2.2：建立路径点(e，sinγ)所需满足的能量e与航迹角正弦sinγ的数学关系。
[0095] 所述的能量e与航迹角正弦sinγ的数学关系根据工程需要确定，可以为公式(30)或二次曲线关系。
[0096]
[0097] 其中，(e0,sinγ0)和(ef,sinγf)为sinγ-e空间(航迹角正弦-能量空间)内初始和末端点，(e1,sinγ1)和(e2,sinγ2)为路径规划算法需要确定的剖面控制点。为了保证在转折点处可导，在转折点处采用哈密尔顿插值算法。
[0098] 由于在进入过程中，数值算法需要满足末端位置式(31)这一约束。
[0099] s(ef)＝0 (31)
[0100] 因此，sinγ-e剖面(航迹角正弦-能量剖面)需有单一可变量。为此，令，[0101]
[0102] 其中△γ＝0.1为常数， 为路径规划算法需要确定的变量。
[0103] 步骤2.3：联立式(29)与(30)，即可求得路径点(e，sinγ)。
[0104] 步骤3：求取制导指令σcmd。
[0105] 在确定sinγ-e剖面(航迹角正弦-能量剖面)后，制导指令σcmd可有动力学方程(18)反解得到
[0106]
[0107] 其中，参考航迹角的导数 由下述关系确定
[0108]
[0109] 求得σcmd即为要求取的制导指令。
[0110] 步骤4：每个制导周期重复步骤2、3，直至满足开伞条件。
[0111] 图3给出了采用路径点规划的预测制导方法在sinγ-e空间(航迹角正弦-能量空间)上的轨迹。可以看出经过进入过程中的一次路径点规划，可以使飞行器在整个进入过程不违背路径约束。
[0112] 本发明保护范围不仅局限于实施例，实施例用于解释本发明，凡与本发明在相同原理和构思条件下的变更或修改均在本发明公开的保护范围之内。
