亚轨道运载火箭级间比推力弹道一体化优化方法技术

亚轨道运载火箭级间比推力弹道一体化优化方法，通过遗传算法给出若干组级间比、推力以及滑行时间的参数组合，作为初始种群的个体；分别将每一组参数组合作为输入条件，以交班速度作为目标函数，通过内层优化算法进行求解，得到每组级间比、推力下的最优飞行弹道和对应的交班速度；以交班速度作为每一组参数的适应度，利用遗传算法得到最优运载器级间比、推力组合和对应的飞行弹道。本发明专利技术在满足亚轨道运载火箭飞行过程中的动压、过载、热流约束的情况下，在运载器总装药量不变情况下，尽可能的提高亚轨道运载火箭的交班速度以使有效载荷拥有更大的交班速度。

【技术实现步骤摘要】
亚轨道运载火箭级间比推力弹道一体化优化方法
本专利技术涉及一种应用于亚轨道运载火箭的级间比推力弹道一体化优化方法，属于运载火箭总体参数设计领域。
技术介绍
亚轨道飞行器作为一种全新的空中平台，开拓了飞机与卫星之间存在的广阔真空地带，具有传统飞行器无法比拟的得天独厚的优势。相对于传统的运载火箭，亚轨道运载火箭弹道形式上具有以下特点：1)较低的交班高度，且助推段全程在大气层内飞行。2)接近0°的交班弹道倾角，需要运载器在大气层内迅速转弯。使用针对空间轨道任务设计的传统运载火箭进行亚轨道飞行器的发射任务势必带来成本增加与可靠性降低的问题，随着各国对亚轨道空间开发力度的加大，有必要对现有的运载火箭进行优化设计以适应于亚轨道发射任务，使得在满足亚轨道运载火箭飞行过程中的动压、过载、热流约束的情况下，在运载器总装药量不变情况下，尽可能的提高亚轨道运载火箭的交班速度以使有效载荷拥有更大的交班速度。运载火箭的设计工作涉及到总体、弹道、气动、动力、控制等诸多学科，是一项复杂的系统工程。总体设计及优化作为当中的牵引项目，不仅将为后续的诸多分系统设计提供指导的方向，更为重要的是将决定飞行器所能达到的最优性能。此外，亚轨道运载火箭的弹道设计也同样对有效载荷的性能有着十分重要的意义。任何涉及亚轨道运载火箭的总体优化都离不开经过优化弹道的检验，能否将总体优化设计所赋予飞行器的性能提升完全地发挥出来取决于弹道设计的优劣。弹道设计与总体的设计的关系就如同计算机的软件与硬件的关系一样，二者的最优组合才能带来最佳的性能。目前还没有对传统运载火箭进行优化设计以适应于亚轨道发射任务的方法。专利技术内容本专利技术的技术解决问题是：克服现有技术的不足，提供亚轨道运载火箭级间比推力弹道一体化优化方法，在满足亚轨道运载火箭飞行过程中的动压、过载、热流约束的情况下，在运载器总装药量不变情况下，使有效载荷拥有更大的交班速度。本专利技术的技术解决方案是：亚轨道运载火箭级间比推力弹道一体化优化方法，包括如下步骤：第一步：通过遗传算法给出若干组级间比、推力以及滑行时间的参数组合，作为初始种群的个体；第二步：分别将每一组参数组合作为输入条件，以交班速度作为目标函数，通过内层优化算法进行求解，得到每组级间比、推力下的最优飞行弹道和对应的交班速度；第三步：以第二步得到的交班速度作为每一组参数的适应度；第四步：判断是否满足迭代的终止条件，如果满足终止条件，则将最后一代种群中适应度最大的个体作为优化的最优结果，输出最优运载器级间比、推力组合和对应的飞行弹道；否则执行第五步；第五步：对每一组参数进行选择、交叉、变异遗传操作，然后返回第二步继续执行。所述内层优化算法的实现步骤如下：(2.1)对亚轨道运载火箭动力学模型进行时域变换，得到时域为[-1,1]的亚轨道运载火箭动力学模型；(2.2)将步骤(2.1)得到的亚轨道运载火箭动力学方程中的控制量和状态量进行离散；(2.3)将步骤(2.1)得到的亚轨道运载火箭动力学方程进行离散；(2.4)将过程约束和终端约束进行离散；(2.5)将目标函数进行离散；(2.6)将每一组参数组合作为输入条件，采用序列二次规划方法进行求解，得到每组级间比、推力下的最优飞行弹道和对应的交班速度。所述步骤(2.1)的实现方法为：(3.1)利用如下变换公式对求解的时间区域进行时域变换：其中，τ为[-1,1]之间的变量，t0为亚轨道运载火箭起飞时刻，tf为亚轨道运载火箭交班时刻，t为亚轨道运载火箭从起飞到交班之间的任一时刻；(3.2)对亚轨道运载火箭动力学方程进行时域变换，得到时域为[-1,1]的亚轨道运载火箭动力学方程：式中：v为亚轨道运载火箭速度，θ为当地弹道倾角，φ为射程角，r为地心距，P为发动机推力，L和D分别为亚轨道运载火箭的升力和推力，α为攻角，m为亚轨道运载火箭质量，g为重力加速度,为亚轨道运载火箭发动机质量流量；其中r、φ、v、θ、m为状态量，α为控制量；(3.3)对过程约束条件进行时域变换，得到时域为[-1,1]的动压约束条件和法向过载约束条件如下：动压约束：法向过载约束：其中，为时域[-1,1]的动压约束函数，qmax为亚轨道运载火箭起飞到交班过程中的最大动压，为时域[-1,1]的过载约束函数，nymax为亚轨道运载火箭起飞到交班过程中的最大法向过载；(3.4)对边界约束条件进行时域变换，得到时域为[-1,1]的终端高度约束条件和终端倾角约束条件如下：终端高度约束：终端倾角约束：其中，为时域[-1,1]的高度约束函数，hf为亚轨道运载火箭交班时的高度，为时域[-1,1]的倾角约束函数，θf为亚轨道运载火箭交班时的弹道倾角。所述步骤(2.2)的实现方式为：(4.1)在区间[-1,1]上配置N+1个配点,配点位置为N阶多项式的所有一阶导数的零点以及-1和1，约定τi的值为所有这些配点中由小至大排序的第i个数；(4.2)利用如下公式对亚轨道运载火箭动力学方程中的状态量进行离散：利用如下公式对亚轨道运载火箭动力学方程中的控制量进行离散：式中Lagi(τ)为Lagrange多项式式中，R(τ)为地心距离散后的函数，Φ(τ)为射程角离散后的函数，V(τ)为亚轨道运载火箭速度离散后的函数，Θ(τ)为当地弹道倾角离散后的函数，M(τ)为亚轨道运载火箭质量离散后的函数，A(τ)为攻角离散后的函数。所述步骤(2.3)中，离散后的亚轨道运载火箭动力学方程为：其中，k∈[0,N]，所述步骤(2.4)中，过程约束和终端约束离散后的公式如下：(1)动压约束：(2)法向过载约束：边界约束条件：(1)终端高度约束：(2)终端倾角约束：所述步骤(2.5)的实现过程如下：(7.1)优化目标为使交班点速度最大，所以目标函数为V(t)为亚轨道运载火箭速度随时间变化的函数；(7.2)将目标函数进行时域变换，得到时域为[-1,1]的目标函数：(7.3)对(7.2)的目标函数进行离散，得到：与现有技术相比，本专利技术具有如下有益效果：(1)本专利技术外层采用遗传算法对级间比进行优化，使得运载火箭的理想速度得到提升，增大了运载器自身的潜力。(2)本专利技术通过外层采用遗传算法对推力进行优化，通过内层优化算法对飞行弹道进行优化，使得运载火箭的速度损失减小，提高了运载器的效率。(3)本专利技术通过增大理想速度与减小速度损失两个手段可以使运载火箭的交班速度得到有效地提升，在满足亚轨道运载火箭飞行过程中的动压、过载、热流约束的情况下，在运载器总装药量不变情况下，增大运载火箭的射程能力。(4)本专利技术通过飞行弹道优化和级间比、推力的优化，使得运载火箭能够达到全局最优的性能。附图说明图1为本专利技术流程图；图2为内层优化算法流程图。具体实施方式运载火箭优化的性能指标通常可以归纳为以下几类：1.最大射程或最大交班速度；2.最小起飞质量；3.最大有效载荷；4.最小成本。对于不同类型的运载火箭，它们的重要程度并不相同，根据问题的性质和设计经验，应选择对目标函数影响显著的参数作为设计变量。本专利技术所选择的目标函数是运载火箭总装药量不变情况下的最大交班速度。影响运载火箭性能指标的参数很多。本专利技术以各级发动机装药量推力以及运载火箭的飞行弹道作为优化变量。本专利技术所涉及到的优化问题是一个总体参数与飞行轨迹结合的优化问题，总体参数的优本文档来自技高网...

【技术保护点】
亚轨道运载火箭级间比推力弹道一体化优化方法，其特征在于包括如下步骤：第一步：通过遗传算法给出若干组级间比、推力以及滑行时间的参数组合，作为初始种群的个体；第二步：分别将每一组参数组合作为输入条件，以交班速度作为目标函数，通过内层优化算法进行求解，得到每组级间比、推力下的最优飞行弹道和对应的交班速度；第三步：以第二步得到的交班速度作为每一组参数的适应度；第四步：判断是否满足迭代的终止条件，如果满足终止条件，则将最后一代种群中适应度最大的个体作为优化的最优结果，输出最优运载器级间比、推力组合和对应的飞行弹道；否则执行第五步；第五步：对每一组参数进行选择、交叉、变异遗传操作，然后返回第二步继续执行。

【技术特征摘要】
1.亚轨道运载火箭级间比推力弹道一体化优化方法，其特征在于包括如下步骤：第一步：通过遗传算法给出若干组级间比、推力以及滑行时间的参数组合，作为初始种群的个体；第二步：分别将每一组参数组合作为输入条件，以交班速度作为目标函数，通过内层优化算法进行求解，得到每组级间比、推力下的最优飞行弹道和对应的交班速度；第三步：以第二步得到的交班速度作为每一组参数的适应度；第四步：判断是否满足迭代的终止条件，如果满足终止条件，则将最后一代种群中适应度最大的个体作为优化的最优结果，输出最优运载器级间比、推力组合和对应的飞行弹道；否则执行第五步；第五步：对每一组参数进行选择、交叉、变异遗传操作，然后返回第二步继续执行。2.根据权利要求1所述的亚轨道运载火箭级间比推力弹道一体化优化方法，其特征在于：所述内层优化算法的实现步骤如下：(2.1)对亚轨道运载火箭动力学模型进行时域变换，得到时域为[-1,1]的亚轨道运载火箭动力学模型；(2.2)将步骤(2.1)得到的亚轨道运载火箭动力学方程中的控制量和状态量进行离散；(2.3)将步骤(2.1)得到的亚轨道运载火箭动力学方程进行离散；(2.4)将过程约束和终端约束进行离散；(2.5)将目标函数进行离散；(2.6)将每一组参数组合作为输入条件，采用序列二次规划方法进行求解，得到每组级间比、推力下的最优飞行弹道和对应的交班速度。3.根据权利要求1所述的亚轨道运载火箭级间比推力弹道一体化优化方法，其特征在于：所述步骤(2.1)的实现方法为：(3.1)利用如下变换公式对求解的时间区域进行时域变换：其中，τ为[-1,1]之间的变量，t0为亚轨道运载火箭起飞时刻，tf为亚轨道运载火箭交班时刻，t为亚轨道运载火箭从起飞到交班之间的任一时刻；(3.2)对亚轨道运载火箭动力学方程进行时域变换，得到时域为[-1,1]的亚轨道运载火箭动力学方程：式中：v为亚轨道运载火箭速度，θ为当地弹道倾角，φ为射程角，r为地心距，P为发动机推力，L和D分别为亚轨道运载火箭的升力和推力，α为攻角，m为亚轨道运载火箭质量，g为重力加速度,为亚轨道运载火箭发动机质量流量；其中r、φ、v、θ、m为状态量，α为控制量；(3.3)对过程约束条件进行时域变换，得到时域为[-1,1]的...

【专利技术属性】
技术研发人员：孙石杰，张烨琛，闵昌万，于明星，汪雷，李加玉，徐春铃，袁江波，张宁宁，杨凌霄，
申请(专利权)人：北京临近空间飞行器系统工程研究所，中国运载火箭技术研究院，
类型：发明
国别省市：北京,11