首页 / 专利库 / 飞机类型 / 巡航高度 / 一种天地往返飞行器全轨迹的优化方法

一种天地往返飞行器全轨迹的优化方法

阅读:172发布:2020-05-21

专利汇可以提供一种天地往返飞行器全轨迹的优化方法专利检索,专利查询,专利分析的服务。并且本 发明 涉及一种天地往返 飞行器 全轨迹的优化方法,包括:S1.基于第一模型参数信息建立飞行器巡航轨迹段的单目标优化模型;S2.基于第二模型参数信息建立飞行器剩余轨迹段的多目标优化模型;S3.基于飞行器参数信息建立所述飞行器的飞行器模型;S4.基于所述飞行器模型并采用单目标优化 算法 对所述单目标优化模型进行求解,并获取第一优化结果;S5.基于所述飞行器模型和所述第一优化结果,并采用多目标 优化算法 对所述多目标优化模型进行求解,并获取第二优化结果;S6.汇总所述第一优化结果和第二优化结果生成优化解集库,根据所述解集库生成所述飞行器的最优飞行轨迹方案。,下面是一种天地往返飞行器全轨迹的优化方法专利的具体信息内容。

1.一种天地往返飞行器全轨迹的优化方法,包括:
S1.基于第一模型参数信息建立飞行器巡航轨迹段的单目标优化模型;
S2.基于第二模型参数信息建立飞行器剩余轨迹段的多目标优化模型;
S3.基于飞行器参数信息建立所述飞行器的飞行器模型;
S4.基于所述飞行器模型并采用单目标优化算法对所述单目标优化模型进行求解,并获取第一优化结果;
S5.基于所述飞行器模型和所述第一优化结果,并采用多目标优化算法对所述多目标优化模型进行求解,并获取第二优化结果;
S6.汇总所述第一优化结果和第二优化结果生成优化解集库,根据所述解集库生成所述飞行器的最优飞行轨迹方案。
2.根据权利要求1所述的优化方法,其特征在于,所述第一模型参数信息包括第一优化目标、第一设计变量和第一约束条件;
所述第一优化目标包括巡航时间tcr、巡航高度hcr和巡航赫数Macr;
所述第一设计变量包括巡航初始时刻的高度h0cr、巡航初始时刻的马赫数Ma0cr和巡航期间的飞行αcr;
所述第一约束条件包括巡航马赫数Macr约束、巡航航迹倾斜角γcr约束、动压qd约束和驻点热流Qs峰值约束。
3.根据权利要求2所述的优化方法,其特征在于,所述巡航马赫数Macr约束和所述巡航航迹倾斜角γcr约束满足:
其中, 为最大巡航马赫数, 为最小巡航马赫数, 为最大巡航航迹倾斜角,为最小巡航航迹倾斜角;
所述动压qd约束满足:
其中, 为动压上边界;
所述驻点热流Qs峰值约束满足:
其中,Qsm为驻点热流Qs的峰值, 为驻点热流Qs的峰值上边界。
4.根据权利要求2所述的优化方法,其特征在于,所述第二模型参数信息包括第二优化目标、第二设计变量和第二约束条件;
所述第二优化目标包括驻点热流Qs和轨迹振荡Δγsum,其中,
其中,tf为所述飞行器全轨迹的总飞行时间;
所述第二设计变量包括爬升滑行段飞行攻角αcg和返回段飞行攻角αrt;
所述第二约束条件包括动压qd约束、过载nload约束、驻点热流Qs峰值约束、爬升滑行段终端飞行状态约束和进场着陆初始飞行状态约束。
5.根据权利要求4所述的优化方法,其特征在于,所述动压qd约束满足:
其中, 为动压上边界;
所述过载nload约束满足:
其中, 为过载上边界;
所述驻点热流Qs峰值约束满足:
其中,Qsm为驻点热流Qs的峰值, 为驻点热流Qs的峰值上边界;
所述爬升滑行段终端飞行状态约束满足:
其中,hfcg为所述飞行器爬升滑行段结束时的高度,Mafcg为所述飞行器爬升滑行段结束时的马赫数,h0cr为巡航初始时刻的高度,Ma0cr为巡航初始时刻的马赫数,γfcg为所述飞行器爬升滑行段结束时的爬升航迹倾斜角;
所述进场着陆初始飞行状态约束满足:
其中,hf为飞行器的飞行高度,hfr为飞行器的飞行高度的要求值,vf为飞行器的飞行速度, 为飞行速度vf的上边界, 为飞行速度vf的下边界,γf为所述飞行器返回段结束时的返回航迹倾斜角, 为返回航迹倾斜角γf的上边界, 为返回航迹倾斜角γf的下边界,xf为飞行器的航程, 为航程xf的上边界, 为航程xf的下边界。
6.根据权利要求1或4所述的优化方法,其特征在于,所述飞行器参数信息包括:飞行器初始质量m0、发动机模型参数和气动模型参数;
所述飞行器初始质量m0满足:
m0=m1dry+m1fuel+m2dry+m2fuel
其中,m1dry为固体助推器干重,m1fuel为固体火箭推进剂的质量,m2dry为飞行器本体干重,m2fuel为飞行器本体的燃料质量;
所述发动机模型参数包括发动机推Ft和质量流量m;
所述气动模型参数包括阻力系数CD和升力系数CL。
7.根据权利要求2所述的优化方法,其特征在于,步骤S4中,选取所述巡航时间tcr、所述巡航高度hcr、所述巡航马赫数Macr中的任一进行单目标优化求解,获取所述第一优化结果;
或者,分别对所述巡航时间tcr、所述巡航高度hcr、所述巡航马赫数Macr进行单目标优化求解,获取所述第一优化结果。
8.根据权利要求6所述的优化方法,其特征在于,步骤S5中包括:
S501.目标函数归一化;
S502.搜寻锚点
S503.生成理想超平面,过各所述锚点生成理想超平面;
S504.在所述理想超平面上生成优化参考点;
S505.压缩目标搜索域;
S506.对于每一个所述优化参考点,结合目标搜索域压缩策略,生成其对应的单目标优化问题,并采用单目标优化算法进行求解;
S507.若多目标优化模型为三维以下的多目标优化问题,则直接滤除局部多目标优化解,获取所述第二优化结果。
9.根据权利要求8所述的优化方法,其特征在于,步骤S5中还包括:
S508.若多目标优化模型为三维以上的多目标优化问题,则旋转目标搜索域;
S509.搜寻边界优化点;
S510.增添优化点;
S511.滤除局部多目标优化解,获取所述第二优化结果。
10.根据权利要求9所述的优化方法,其特征在于,步骤S509中,基于二分法,根据旋转后的所述目标搜索域,搜寻目标空间内位于边界的边界优化点。

说明书全文

一种天地往返飞行器全轨迹的优化方法

技术领域

[0001] 本发明涉及一种天地往返飞行器全轨迹的优化方法,尤其涉及一种高超声速天地往返飞行器全轨迹的优化方法。

背景技术

[0002] 随着超燃冲压发动机技术的兴起和发展,高超声速的巡航飞行器逐渐进入科研人员的研究视线。其推进效率相比于传统火箭得到大幅提高,可在20km以上的高度进行
数大于5的高超声速巡航飞行,并返回地面的规定机场进行着陆。这种类型的高超声速巡航飞行器,具有发动机效率高、巡航速度更快、可重复使用等优势,在未来的空基发射入轨领域具有很好的发展前景。
[0003] 对于高超声速巡航飞行器,其整个飞行轨迹可划分为上升段、巡航段和下降段等多个飞行阶段。各飞行阶段均包含了众多不同类型的约束条件,且飞行器的状态不尽相同,发动机、气动与轨迹的耦合效应较为明显。在飞行器的整个飞行轨迹中,上升段和下降段分别以巡航飞行状态作为终端飞行状态约束和飞行状态初值,来进行各自轨迹的优化。因此,巡航飞行状态对全轨迹的影响十分直接和重要。目前已有对飞行器飞行轨迹优化的方法大
都事先人为确定了一个固定的巡航高度、马赫数等巡航飞行状态,这很大程度的限制了整
个飞行轨迹的优化设计空间。

发明内容

[0004] 本发明的目的在于提供一种天地往返飞行器全轨迹的优化方法,可以在满足一定轨迹约束条件的要求下,给出巡航飞行任务目标最优和飞行技术难度指标矛盾权衡优化时
的优化解集库。用于天地往返巡航飞行器的轨迹优化设计和飞行方案设计。
[0005] 为实现上述发明目的,本发明提供一种天地往返飞行器全轨迹的优化方法,包括:
[0006] S1.基于第一模型参数信息建立飞行器巡航轨迹段的单目标优化模型;
[0007] S2.基于第二模型参数信息建立飞行器剩余轨迹段的多目标优化模型;
[0008] S3.基于飞行器参数信息建立所述飞行器的飞行器模型;
[0009] S4.基于所述飞行器模型并采用单目标优化算法对所述单目标优化模型进行求解,并获取第一优化结果;
[0010] S5.基于所述飞行器模型和所述第一优化结果,并采用多目标优化算法对所述多目标优化模型进行求解,并获取第二优化结果;
[0011] S6.汇总所述第一优化结果和第二优化结果生成优化解集库,根据所述解集库生成所述飞行器的最优飞行轨迹方案。
[0012] 根据本发明的一个方面,所述第一模型参数信息包括第一优化目标、第一设计变量和第一约束条件;
[0013] 所述第一优化目标包括巡航时间tcr、巡航高度hcr和巡航马赫数Macr;
[0014] 所述第一设计变量包括巡航初始时刻的高度h0cr、巡航初始时刻的马赫数Ma0cr和巡航期间的飞行αcr;
[0015] 所述第一约束条件包括巡航马赫数Macr约束、巡航航迹倾斜角γcr约束、动压qd约束和驻点热流Qs峰值约束。
[0016] 根据本发明的一个方面,所述巡航马赫数Macr约束和所述巡航航迹倾斜角γcr约束满足:
[0017]
[0018] 其中, 为最大巡航马赫数, 为最小巡航马赫数, 为最大巡航航迹倾斜角, 为最小巡航航迹倾斜角;
[0019] 所述动压qd约束满足:
[0020]
[0021] 其中, 为动压上边界;
[0022] 所述驻点热流Qs峰值约束满足:
[0023]
[0024] 其中,Qsm为驻点热流Qs的峰值, 为驻点热流Qs的峰值上边界。
[0025] 根据本发明的一个方面,所述第二模型参数信息包括第二优化目标、第二设计变量和第二约束条件;
[0026] 所述第二优化目标包括驻点热流Qs和轨迹振荡Δγsum,其中,
[0027]
[0028] 其中,tf为所述飞行器全轨迹的总飞行时间;
[0029] 所述第二设计变量包括爬升滑行段飞行攻角αcg和返回段飞行攻角αrt;
[0030] 所述第二约束条件包括动压qd约束、过载nload约束、驻点热流Qs峰值约束、爬升滑行段终端飞行状态约束和进场着陆初始飞行状态约束。
[0031] 根据本发明的一个方面,所述动压qd约束满足:
[0032]
[0033] 其中, 为动压上边界;
[0034] 所述过载nload约束满足:
[0035]
[0036] 其中, 为过载上边界;
[0037] 所述驻点热流Qs峰值约束满足:
[0038]
[0039] 其中,Qsm为驻点热流Qs的峰值, 为驻点热流Qs的峰值上边界;
[0040] 所述爬升滑行段终端飞行状态约束满足:
[0041]
[0042] 其中,hfcg为所述飞行器爬升滑行段结束时的高度,Mafcg为所述飞行器爬升滑行段结束时的马赫数,h0cr为巡航初始时刻的高度,Ma0cr为巡航初始时刻的马赫数,γfcg为所述飞行器爬升滑行段结束时的爬升航迹倾斜角;
[0043] 所述进场着陆初始飞行状态约束满足:
[0044]
[0045] 其中,hf为飞行器的飞行高度,hfr为飞行器的飞行高度的要求值,vf为飞行器的飞行速度, 为飞行速度vf的上边界, 为飞行速度vf的下边界,γf为所述飞行器返回段结束时的返回航迹倾斜角, 为返回航迹倾斜角γf的上边界, 为返回航迹倾斜角γf的下
边界,xf为飞行器的航程, 为航程xf的上边界, 为航程xf的下边界。
[0046] 根据本发明的一个方面,所述飞行器参数信息包括:飞行器初始质量m0、发动机模型参数和气动模型参数;
[0047] 所述飞行器初始质量m0满足:
[0048] m0=m1dry+m1fuel+m2dry+m2fuel
[0049] 其中,m1dry为固体助推器干重,m1fuel为固体火箭推进剂的质量,m2dry为飞行器本体干重,m2fuel为飞行器本体的燃料质量;
[0050] 所述发动机模型参数包括发动机推Ft和质量流量m;
[0051] 所述气动模型参数包括阻力系数CD和升力系数CL。
[0052] 根据本发明的一个方面,步骤S4中,选取所述巡航时间tcr、所述巡航高度hcr、所述巡航马赫数Macr中的任一进行单目标优化求解,获取所述第一优化结果;或者,分别对所述巡航时间tcr、所述巡航高度hcr、所述巡航马赫数Macr进行单目标优化求解,获取所述第一优化结果。
[0053] 根据本发明的一个方面,步骤S5中包括:
[0054] S501.目标函数归一化;
[0055] S502.搜寻锚点
[0056] S503.生成理想超平面,过各所述锚点生成理想超平面;
[0057] S504.在所述理想超平面上生成优化参考点;
[0058] S505.压缩目标搜索域;
[0059] S506.对于每一个所述优化参考点,结合目标搜索域压缩策略,生成其对应的单目标优化问题,并采用单目标优化算法进行求解;
[0060] S507.若多目标优化模型为三维以下的多目标优化问题,则直接滤除局部多目标优化解,获取所述第二优化结果。
[0061] 根据本发明的一个方面,步骤S5中还包括:
[0062] S508.若多目标优化模型为三维以上的多目标优化问题,则旋转目标搜索域;
[0063] S509.搜寻边界优化点;
[0064] S510.增添优化点;
[0065] S511.滤除局部多目标优化解,获取所述第二优化结果。
[0066] 根据本发明的一个方面,步骤S509中,基于二分法,根据旋转后的所述目标搜索域,搜寻目标空间内位于边界的边界优化点。
[0067] 根据本发明的一种方案,在满足设计约束要求的前提下,设计优化得到令巡航飞行任务目标最优的巡航轨迹,并得到以此最优巡航状态为基础的能够尽量降低整个飞行技
术难度的剩余段飞行轨迹库。
[0068] 根据本发明的一种方案,巡航段作为高超声速巡航飞行器的主要飞行阶段,其对应的巡航飞行状态应当根据飞行任务需求和约束,采用优化的手段设计确定,从而更全面
的挖掘巡航段乃至整个飞行轨迹的设计优化空间,得到更具意义的优化解。同时,采用本发明的方法对飞行器全轨迹分步与多目标优化,实现了巡航飞行任务目标的优化和飞行技术
难度优化。
附图说明
[0069] 图1示意性表示根据本发明的一种实施方式的优化方法的步骤框图
[0070] 图2示意性表示根据本发明的一种实施方式的优化方法的飞行器全轨迹分步优化关系图;
[0071] 图3示意性表示根据本发明的一种实施方式优化方法的多目标优化方法步骤流程图;
[0072] 图4示意性表示根据本发明的一种实施方式的优化方法的剩余段轨迹多目标优化点集图;
[0073] 图5示意性表示根据本发明的一种实施方式的优化方法的高度-时间曲线对比图。

具体实施方式

[0074] 为了更清楚地说明本发明实施方式或现有技术中的技术方案,下面将对实施方式中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是本发明的一些
实施方式,对于本领域普通技术人员而言,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
[0075] 在针对本发明的实施方式进行描述时,术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“平”、“顶”、“底”“内”、“外”所表达的方位或位置关系是基于相关附图所示的方位或位置关系,其仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此上述术语不能
理解为对本发明的限制。
[0076] 下面结合附图和具体实施方式对本发明作详细地描述,实施方式不能在此一一赘述,但本发明的实施方式并不因此限定于以下实施方式。
[0077] 如图1所示,根据本发明的一种实施方式,本发明的一种天地往返飞行器全轨迹的优化方法,包括:
[0078] S1.基于第一模型参数信息建立飞行器巡航轨迹段的单目标优化模型;
[0079] S2.基于第二模型参数信息建立飞行器剩余轨迹段的多目标优化模型;
[0080] S3.基于飞行器参数信息建立飞行器的飞行器模型;
[0081] S4.基于飞行器模型并采用单目标优化算法对单目标优化模型进行求解,并获取第一优化结果;
[0082] S5.基于飞行器模型和第一优化结果,并采用多目标优化算法对多目标优化模型进行求解,并获取第二优化结果;
[0083] S6.汇总第一优化结果和第二优化结果生成优化解集库,根据解集库生成飞行器的最优飞行轨迹方案。
[0084] 根据本发明的一种实施方式,在步骤S1中建立飞行器巡航轨迹段的单目标优化模型。在本实施方式中,通过第一模型参数信息建立飞行器巡航轨迹段的单目标优化模型。其中,第一模型参数信息包括第一优化目标、第一设计变量和第一约束条件。在本实施方式
中,首先设置巡航段的第一优化目标,第一优化目标包括巡航时间tcr、巡航高度hcr和巡航马赫数Macr;然后设置飞行器巡航段的第一设计变量,第一设计变量包括巡航初始时刻的高度h0cr、巡航初始时刻的马赫数Ma0cr和巡航期间的飞行攻角αcr;最后整理并设置飞行器巡航段优化的第一约束条件,第一约束条件包括巡航马赫数Macr约束、巡航航迹倾斜角γcr约
束、动压qd约束和驻点热流Qs峰值约束。在本实施方式中,巡航飞行器采用稳态巡航,期间高度、马赫数变化不大。因此,第一优化目标中的巡航高度hcr和巡航马赫数Macr可采用巡航开始时刻的巡航初始时刻的高度h0cr、巡航初始时刻的马赫数Ma0cr等表征整个巡航期间的相应飞行状态信息。在本实施方式中,第一优化目标需要实现巡航时间tcr最长,巡航初始时刻的高度h0cr最高,巡航初始时刻的马赫数Ma0cr最大。在实际应用中,应从上述第一优化目标中确定一个进行单目标优化。或这分别对第一优化目标中的3个分别进行单目标优化。
[0085] 在本实施方式中,在本步骤S1中,设置第一步优化的设计变量为:巡航初始时刻的高度h0cr、巡航初始时刻的马赫数Ma0cr和巡航期间的飞行攻角αcr。由于不同飞行器的任务状态不同,因此需要根据实际需求确定各变量的取值范围。
[0086] 在本实施方式中,在整理并设置飞行器巡航段的约束条件时,巡航马赫数Macr约束和巡航航迹倾斜角γcr约束根据稳态巡航的推阻平衡和升重平衡条件,可知飞行器在巡航期间巡航马赫数Macr、巡航航迹倾斜角γcr的理论变化量应为零,但在实际飞行中很难保证上述变化量绝对为零,需要给定一个很小的容许变化范围,因此,巡航马赫数Macr约束和巡航航迹倾斜角γcr约束满足:
[0087]
[0088] 其中, 为最大巡航马赫数, 为最小巡航马赫数, 为最大巡航航迹倾斜角, 为最小巡航航迹倾斜角;
[0089] 在本实施方式中,在实际巡航飞行过程中,飞行器的飞行动压应小于动压约束上边界,因此动压qd约束满足:
[0090]
[0091] 其中, 为动压上边界;
[0092] 在本实施方式中,在实际巡航飞行过程中,飞行器的驻点热流Qs的最大值,即驻点峰值热流Qsm应小于上边界,驻点热流Qs峰值约束满足:
[0093]
[0094] 其中,Qsm为驻点热流Qs的峰值, 为驻点热流Qs的峰值上边界。
[0095] 根据本发明的一种实施方式,步骤S2中建立飞行器剩余轨迹段的多目标优化模型。在本实施方式中,通过第二模型参数信息建立飞行器剩余轨迹段的多目标优化模型。剩余轨迹段包括飞行器的爬升滑行轨迹段和返回轨迹段。第二模型参数信息包括第二优化目
标、第二设计变量和第二约束条件。在本实施方式中,首先设置剩余轨迹段的第二优化目
标,第二优化目标包括驻点热流Qs和轨迹振荡Δγsum;然后设置飞行器剩余轨迹段的第二设计变量,第二设计变量包括爬升滑行段飞行攻角αcg和返回段飞行攻角αrt;最后整理并设置飞行器剩余轨迹段优化的第二约束条件,第二约束条件包括动压qd约束、过载nload约束、驻点热流Qs峰值约束、爬升滑行段终端飞行状态约束和进场着陆初始飞行状态约束。在本
实施方式中,第二优化目标目的是减小飞行器的爬升滑行轨迹段和返回轨迹段的飞行难
度,在保证达到规定巡航飞行状态要求的前提下,减小整个飞行任务执行的难度。在本实施方式中,从驻点热流Qs的驻点热流峰值Qsm和轨迹振荡Δγsum程度两方面评价飞行难度,并作为优化目标考虑。因此,在剩余轨迹段第二优化目标的优化为二维多目标优化,并且使第二优化目标中的驻点热流Qs的驻点热流峰值Qsm和轨迹振荡Δγsum满足:
[0096] 驻点热流峰值Qsm最低,Min Qsm;
[0097] 轨迹振荡Δγsum最小,用航迹倾斜角的全程累积变化量表示,MinΔγsum,其满足:
[0098]
[0099] 其中,tf为所述飞行器全轨迹的总飞行时间。
[0100] 在本实施方式中,在本步骤S2中,设置第二步优化的设计变量为:爬升滑行段飞行攻角αcg和返回段飞行攻角αrt。由于不同飞行器的任务状态不同,需要根据实际需求确定各变量的取值范围。
[0101] 在本实施方式中,在整理并设置飞行器剩余轨迹段的约束条件时,分别对动压qd约束、过载nload约束、驻点热流Qs峰值约束、爬升滑行段终端飞行状态约束和进场着陆初始飞行状态约束进行整理并设置。
[0102] 在本实施方式中,在实际剩余轨迹段飞行过程中,飞行器的飞行动压应小于动压约束上边界,因此动压qd约束满足:
[0103]
[0104] 其中, 为动压上边界;
[0105] 在本实施方式中,在实际剩余轨迹段飞行过程中,飞行器的过载nload应小于过载约束上边界,因此过载nload约束满足:
[0106]
[0107] 其中, 为过载上边界;
[0108] 在本实施方式中,在实际剩余轨迹端飞行过程中,飞行器的驻点热流Qs的最大值,即驻点峰值热流Qsm应小于上边界,驻点热流Qs峰值约束满足:
[0109]
[0110] 其中,Qsm为驻点热流Qs的峰值, 为驻点热流Qs的峰值上边界。
[0111] 在本实施方式中,飞行器爬升滑行段结束时,其高度和马赫数必须与前述步骤S1中得到的巡航段入口位置的巡航初始时刻的高度h0cr和巡航初始时刻的马赫数Ma0cr相等,且此时航迹倾斜角应基本为0,但在实际飞行中很难保证上述航迹倾斜角绝对为零,这里同样设置一个很小的容许范围。因此,爬升滑行段终端飞行状态约束满足:
[0112]
[0113] 其中,hfcg为所述飞行器爬升滑行段结束时的高度,Mafcg为所述飞行器爬升滑行段结束时的马赫数,h0cr为巡航初始时刻的高度,Ma0cr为巡航初始时刻的马赫数,γfcg为所述飞行器爬升滑行段结束时的爬升航迹倾斜角;
[0114] 在本实施方式中,飞行器在返回段结束时,须满足一定的关于高度、速度、航迹倾斜角和航程的进场着陆条件,因此进场着陆初始飞行状态约束满足:
[0115]
[0116] 其中,hf为飞行器的飞行高度,hfr为飞行器的飞行高度的要求值,vf为飞行器的飞行速度, 为飞行速度vf的上边界, 为飞行速度vf的下边界,γf为所述飞行器返回段结束时的返回航迹倾斜角, 为返回航迹倾斜角γf的上边界, 为返回航迹倾斜角γf的下
边界,xf为飞行器的航程, 为航程xf的上边界, 为航程xf的下边界。
[0117] 根据本发明的一种实施方式,步骤S3中,建立飞行器的飞行器模型。在本实施方式中,根据飞行器参数信息建立飞行器的飞行器模型。其中,飞行器参数信息包括:飞行器初始质量m0、发动机模型参数和气动模型参数。在本实施方式中,由于飞行器包含固体火箭助推器(包含过渡段)和天地往返飞行器本体两部分,整个组合体在发射时刻的质量,即飞行器初始质量m0满足:
[0118] m0=m1dry+m1fuel+m2dry+m2fuel
[0119] 其中,m1dry为固体助推器干重(包含过渡段质量),m1fuel为固体火箭推进剂的质量,m2dry为飞行器本体干重,m2fuel为飞行器本体的燃料质量。
[0120] 在本实施方式中,发动机模型参数包括发动机推力Ft和质量流量m。在本实施方式中,飞行器为高超声速飞行器,其发送机采用冲压发动机,因此需要建立发动机推力Ft和质量流量m随马赫数Ma、高度h、攻角α等参数的变化关系:
[0121] Ft=ft(Ma,h,α)
[0122] m=fm(Ma,h,α)
[0123] 其中,ft和fm分别为根据已知发动机数据,得到发动机推力和质量流量的插值函数。
[0124] 在本实施方式中,气动模型参数包括阻力系数CD和升力系数CL。在本实施方式中,飞行器在高超声速飞行过程中需要建立阻力系数CD和升力系数CL随马赫数Ma和攻角α的变化关系:
[0125] CD=fD(Ma,α)
[0126] CL=fL(Ma,α)
[0127] 其中,fD和fL分别为根据已知气动参数二维插值表,得到阻力系数CD和升力系数CL的插值函数。
[0128] 根据本发明的一种实施方式,步骤S4中,基于飞行器模型并采用单目标优化算法对单目标优化模型进行求解,并获取第一优化结果。在本实施方式中,采用单目标优化算法序列二次规划法(Sequential Quadratic Programming,SQP)求解单目标优化模型中飞行
器巡航轨迹端的优化问题。在本实施方式中,可选取巡航时间tcr、巡航高度hcr、巡航马赫数Macr中的任一一项通过单目标优化算法进行单目标优化求解,获取第一优化结果。或者,分别对巡航时间tcr、巡航高度hcr、巡航马赫数Macr通过单目标优化算法进行单目标优化求解,获取所述第一优化结果。
[0129] 如图2所示,根据本发明的一种实施方式,步骤S5中,根据步骤S4中获取的第一优化结果对飞行器巡航轨迹段进行优化,基于优化后的巡航轨迹段对飞行器全轨迹巡航端进
行分步优化,即剩余轨迹段进行优化,参见图2所示。在本实施方式中,需要针对驻点热流峰值Qsm和轨迹振荡总量Δγsum两个优化目标进行多目标优化。在本实施方式中,采用多目标优化算法第三代直接搜索域法(Directed Search Domain III,DSD-III),求解多目标优化
模型剩余轨迹段优化问题。在本实施方式中,采用优化求解策略将多目标优化问题转换为
单目标优化问题进行求解,对所有优化解进行汇总和处理后,得到多目标优化解,即第二优化结果。
[0130] 如图3所示,根据本发明的一种实施方式,步骤S5中,基于飞行器模型和第一优化结果,并采用多目标优化算法对多目标优化模型进行求解,并获取第二优化结果。在本实施方式中,步骤S5中包括:
[0131] S501.目标函数归一化。
[0132] S502.搜寻锚点,在本实施方式中,锚点即为单目标最优时的优化点。
[0133] S503.生成理想超平面,过各所述锚点生成理想超平面。
[0134] S504.在所述理想超平面上生成优化参考点。
[0135] S505.压缩目标搜索域,在本实施方式中,多目标优化算法第三代直接搜索域法集成了原始第一代直接搜索域法(DSD)和第二代直接搜索域法(DSD-II)算法的相应策略,供
针对具体问题进行选用。在本实施方式中,优先选用第二代直接搜索域法(DSD-II)的压缩
策略,但优化问题中存在初值敏感性问题时,可选择第一代直接搜索域法(DSD)算法的压缩策略。
[0136] S506.对于每一个所述优化参考点,结合前述步骤S505中相应的目标搜索域压缩策略,生成其对应的单目标优化问题,并采用单目标优化算法进行求解。
[0137] S507.若多目标优化模型为三维以下的多目标优化问题,则直接滤除局部多目标优化解,获取第二优化结果。
[0138] 如图3所示,在本实施方式中,步骤S5中还包括:
[0139] S508.若多目标优化模型为三维以上的多目标优化问题,则旋转目标搜索域。在本实施方式中,通过调用旋转搜索策略,对前述步骤中压缩后的搜索域进行旋转,从而保证优化搜索能够较全面的覆盖整个优化空间。
[0140] S509.搜寻边界优化点,在本实施方式中,第三代直接搜索域法基于二分法,根据旋转后的目标搜索域,搜寻目标空间内位于边界的边界优化点。
[0141] S510.增添优化点,在本实施方式中,根据第三代直接搜索域法关于优化点是否需要增加的判断条件,按照一定的计算和优化方法,新增优化问题和得到的优化点。
[0142] S511.滤除局部多目标优化解,获取第二优化结果,在本实施方式中,第三代直接搜索域法集成了两种方法滤除局部的优化解,以保证优化解的全局最优性。
[0143] 根据本发明的一种实施方式,步骤S6中对第一优化结果和第二优化结果进行汇总生成优化解集库。在本实施方式中,根据解集库生成飞行器的最优飞行轨迹方案。例如,参见表1和表2,表1中数据为以一个算例得到的飞行器在巡航轨迹段巡航时间最长的优化轨
迹飞行状态。表2中数据为通过汇总步骤S507和步骤S511中获取的第二优化结果从而生成
飞行器剩余轨迹段多目标优化解集库。同时,根据表2的飞行器剩余轨迹段多目标优化解集库生成剩余段轨迹多目标优化点集如图4所示。如图5所示,在本实施中,可选取表2中优化解标号为1和19的两组数据生成图5中高度-时间曲线对比图。通过表1和表2中的数据进行
飞行器巡航轨迹段最优解(即第一步优化结果)和剩余轨迹段最优解(即第二步优化结果)
的组合从而构成了飞行器全轨迹的优化结果。
[0144]
[0145] 表1最优巡航轨迹算例结果
[0146]
[0147] 表2剩余段多目标优化解集库
[0148] 上述内容仅为本发明的具体方案的例举,对于其中未详尽描述的设备和结构,应当理解为采取本领域已有的通用设备及通用方法来予以实施。
[0149] 以上所述仅为本发明的一个方案而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修
改、等同替换、改进等,均应包含在本发明的保护范围之内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈