CN113589847A - 一种柔性飞行器飞行半径确定方法 - Google Patents

一种柔性飞行器飞行半径确定方法 Download PDF

Info

Publication number
CN113589847A
CN113589847A CN202111042267.3A CN202111042267A CN113589847A CN 113589847 A CN113589847 A CN 113589847A CN 202111042267 A CN202111042267 A CN 202111042267A CN 113589847 A CN113589847 A CN 113589847A
Authority
CN
China
Prior art keywords
aircraft
climbing
section
engine
oil consumption
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111042267.3A
Other languages
English (en)
Other versions
CN113589847B (zh
Inventor
李道春
邵浩原
赵仕伟
阚梓
申童
姚卓尔
向锦武
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202111042267.3A priority Critical patent/CN113589847B/zh
Publication of CN113589847A publication Critical patent/CN113589847A/zh
Application granted granted Critical
Publication of CN113589847B publication Critical patent/CN113589847B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10—Simultaneous control of position or course in three dimensions
    • G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • G05D1/106—Change initiated in response to external conditions, e.g. avoidance of elevated terrain or of no-fly zones

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

本发明公开了一种柔性飞行器飞行半径确定方法，针对飞机燃油消耗问题，建立燃油消耗量与平均高度、平均速度、飞行时间、油箱平均油量以及发动机工作状态的多元回归模型，通过对飞行器任务剖面的划分，基于柔性飞机动力学，修正飞机模型气动数据，该方法对柔性飞行器飞行半径的计算快捷简单，可以保证柔性飞行器在任务剖面内作战半径计算的准确性。该方法可为柔性飞行器的飞行性能计算提供参考，为解决柔性飞行器的航程计算，确定飞行性能等提供解决方案。

Description

一种柔性飞行器飞行半径确定方法
技术领域
本发明涉及航空技术领域，尤其涉及一种柔性飞行器飞行半径确定方法。
背景技术
飞行半径是衡量飞机性能的重要指标。对于飞行器而言，飞行半径计算精度与飞行器重量、载油量等因素有关。精确的飞行半径计算，是飞行性能评估和任务规划的重要环节。在设计之初由于缺乏完整详尽的计算分析数据，需要对飞行器对飞行器起飞前，需要精确估算全航线航程，以便计算飞行器的飞行时间、机动突防距离等战术参量；飞行中，导航系统必须实时精确的计算当前飞行航程，才能保证任务剖面的正确性和时效性。
目前，在作战半径计算时，往往受限于飞行力学建模对象的刚体假设等因素，对于具有一定柔性的飞行器(如运输机、大展弦比无人机)而言，以此为基础的作战半径计算往往产生较大误差。因此本发明提供了一种计算快捷简单，满足一定精度要求的柔性飞行器飞行半径计算方法。该方法无需专业性能软件或与仿真软件所用的大量原始数据。
发明内容
有鉴于此，本发明提供了一种柔性飞行器飞行半径确定方法，包括如下步骤：
S1：建立发动机燃油消耗率与平均高度、平均速度、发动机推力、发动机设定工作状态的多元线性回归模型，即为发动机燃油消耗率模型；
发动机推力与高度、飞行马赫数、发动机设定工作状态的线性回归模型，即为发动机推力模型；
一般情况下，发动机工作状态可以分为慢车状态、加力状态和中间推力状态。数据可由实际飞行试验测得，也可由发动机模拟试验测得，一般为离散的数据点，需要进行多元线性回归。
S2：根据高度，建立温度、声速、压力、大气密度的数学模型；
S3：根据建立考虑柔性影响修正的升力系数模型，根据公式
Figure BDA0003249762450000011
其中ηi为i阶弹性模态，
Figure BDA0003249762450000012
为弹性模态对应的升力增量。确定飞行器升力系数，并修正飞行器极曲线的气动数据，用于校正后续计算过程。
S4：根据已知飞行器起飞重量、停机角、起飞机场高度、温度、离地迎角和跑道摩擦系数，计算起飞耗油量。起飞耗油量包括起飞前暖机耗油量和起飞滑跑耗油量；具体步骤为：
S41.根据起飞暖机时间和发动机暖机耗油率，计算起飞暖机耗油；根据起飞机场高度，代入大气密度的数学模型，修正大气密度。
S42.根据设置的离地迎角和飞行器的升力系数曲线，求得离地升力系数，进而由飞行器极曲线求得离地阻力系数。
S43.迭代计算飞行器离地速度与发动机离地推力。方法如下：计算飞行器离地速度。通过发动机推力模型，求得发动机离地推力。进而更新得到新的离地速度。若两次计算的离地速度满足误差范围，结果为飞行器离地速度，若不满足误差范围，则用新的离地速度通过发动机推力模型，迭代求解，直至满足误差范围。
S44.根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器初始速度、发动机初始推力，求得起飞初始耗油率。根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器离地速度、发动机离地推力，求得起飞离地耗油率。计算起飞时间
Figure BDA0003249762450000021
其中g为当地重力加速度，f为跑道摩擦系数，
Figure BDA0003249762450000022
为离地升阻比，Gtakeoff为飞行器起飞重量，T0为初始推力，Tqf为起飞离地推力，VLD为离地速度。取初始耗油率和起飞离地耗油率的平均值作为起飞过程耗油率，结合起飞时间，计算得到起飞耗油量。
S5：根据场高等条件，计算爬升至任务剖面规定巡航高度所需爬升高度，估算爬升过程耗油量，分为加速爬升和等速爬升两个阶段，分别对应加速爬升耗油量和等速爬升耗油量；具体步骤为：
S51.确定起飞爬升初始重量、爬升开始高度和爬升终点高度。一般的，爬升开始速度为飞行器起飞离地速度，爬升终点高度为飞行器巡航速度，爬升开始高度为机场高度，爬升终点高度为飞行器出航巡航高度。
S52.根据爬升高度和爬升起止速度，将爬升段平均分为N段进行计算。根据每段爬升的平均高度，代入大气密度与声速数学模型，修正该段高度的大气密度、声速。
S53.迭代计算每段爬升的飞行器速度、马赫数、爬升速度、飞行器重量、耗油量。计算该爬升段飞行器升力系数，并通过飞机极曲线求得对应阻力系数。由马赫数、该爬升段平均高度和发动机工作状态，代入所述发动机推力模型，得到发动机推力。根据该段发动机设定的工作状态、飞行速度、发动机推力和平均高度，计算燃油消耗率。计算动能变化修正系数
Figure BDA0003249762450000023
爬升速度
Figure BDA0003249762450000024
该段爬升时间
Figure BDA0003249762450000025
该爬升段燃油消耗量Wi＝sfci*DT，该爬升段水平距离，其中Vi为计算段飞行速度，DV为计算段速度变化量，DH为计算段高度变化量，Gps为爬升段重量，Ti为爬升段推力，Di为爬升段阻力，sfci为计算段燃油消耗率。更新飞行器速度、马赫数、爬升速度、飞行器重量，进行下一段爬升耗油量计算。
S54.累加每一段爬升的耗油量，得爬升段耗油量，累加每一段爬升的水平距离，得到爬升段水平距离。
额外的，如爬升段有发动机加力工作情况，需按照发动机加力工作状态，则需要计算给定加力时长、飞行器高度、飞行器推力的发动机耗油量，累加到爬升段总耗油量中。
S6：根据任务剖面规定，结合发动机工作状态、任务要求推力、飞行高度和飞行速度，结合S1所述发动机燃油消耗率模型，结合巡逻、作战等任务段时间，计算任务段耗油量；
S7：出航巡航计算。确定飞行器安全着陆余油，根据前述计算的起飞耗油量、爬升过程耗油量、任务段耗油量，确定巡航可用余油，并以巡航可用余油的一半确定出航可用油量；由巡航段所需推力、设定巡航速度、设定巡航高度推算巡航耗油率。具体步骤如下：
S71.将计算段油量分为N个计算段，确定出航巡航开始前飞行器重量，根据出航巡航段高度，修正出航巡航段马赫数、大气密度、大气压力，求得动压。
S72.确定出航巡航段的计算段发动机设定的工作状态、飞行器升力系数，由飞行器极曲线得到阻力系数、飞行器升阻比。结合发动机推力损失计算巡航所需推力。
S73.根据该巡航段发动机设定的工作状态、飞行速度、发动机推力和平均高度，计算燃油消耗率，计算段航程
Figure BDA0003249762450000031
其中sfci为计算段燃油消耗率，DW为计算段燃油消耗量，V为巡航速度。更新飞行器重量，进行下一计算段巡航耗油量计算。
S74.累加每一段的航程，得到出航巡航距离。
特别的，计算段数量N的确定由飞行器设计航程确定，一般航程越远，所设置N的值越大。当飞行器有额外挂载时，第S72步的飞行器阻力系数需要增加飞行器挂载的阻力系数增量。第S73步，由计算段耗油量和耗油率，累加可以计算得到航时，同时也可以根据计算段航程与耗油量，得到单位航程耗油量，进行远航速度计算。
S8：根据具体任务需要，确定飞行器返航构型、返航巡航高度，修正飞机返航重量、马赫数、大气密度、大气压力和升阻力系数，结合S1所述发动机燃油消耗率模型，计算方法同S7所述出航巡航计算，最终得到返航巡航距离；
S9：出航巡航距离与返航巡航距离比较，若误差较大，则适当增加或减少出航油量，迭代计算。若误差合理，则爬升水平距离与出/返航距离累加得作战半径。
附图说明
图1为本发明的流程图；
图2为本发明实施例中的大展弦比固定翼无人机作战半径示意图。
具体实施方式
下面将结合本发明实施方式中的附图，对本发明实施方式中的技术方案进行清楚、完整的描述，显然，所描述的实施方式仅仅是作为例示，并非用于限制本发明。
下面通过一个具体的实施例对本发明提供的上述柔性飞行器飞行半径计算方法的具体实施进行详细说明。
实施例1为某大展弦比柔性飞行器飞行半径计算：
如图1所示，包括如下步骤：
S1：建立发动机燃油消耗率与平均高度、平均速度、发动机推力、发动机设定工作状态的多元线性回归模型，即为发动机燃油消耗率模型；
发动机推力与高度、飞行马赫数、发动机设定工作状态的线性回归模型，即为发动机推力模型；
所述发动机燃油消耗率模型表示为：
Sfc＝f(H,V,T,state)
发动机推力模型表示为：
T＝T(H,Ma,state)
其中H为发动机/飞行器平均高度，V为飞行器平均速度，Ma为飞行器飞行马赫数，T为发动机推力，state为发动机工作状态，一般情况下，发动机工作状态可以分为慢车状态、加力状态和中间推力状态。数据可由实际飞行试验测得，也可由发动机模拟试验测得，一般为离散的数据点，需要进行多元线性回归。
S2：根据气压高度，建立温度、声速、压力、大气密度的数学模型；所述计算数学模型，分别如下：
温度：
Figure BDA0003249762450000041
声速：
Figure BDA0003249762450000042
压力：
Figure BDA0003249762450000051
大气密度：
Figure BDA0003249762450000052
S3：根据建立考虑柔性影响修正的升力系数模型，根据公式
Figure BDA0003249762450000053
其中ηi为i阶弹性模态，
Figure BDA0003249762450000054
为弹性模态对应的升力增量。确定并修正飞行器升力系数，修正飞行器气动数据，用于校正后续计算过程。
S4：根据已知飞行器起飞重量、停机角、起飞机场高度、温度、离地迎角和跑道摩擦系数，计算起飞耗油量。起飞耗油量Wqf包括起飞前暖机耗油量W0和起飞滑跑耗油量Whp；具体步骤为：
S41.根据起飞暖机时间T0和发动机暖机耗油率sfc0，计算起飞暖机耗油W0＝T0*sfc0；根据起飞机场高度，修正大气密度ρ。
S42.根据设置的离地迎角和飞行器的升力系数曲线，求得离地升力系数CYLD，进而由极曲线求得离地阻力系数CXLD。
S43.迭代计算飞行器离地速度VLD与发动机离地推力TLD。方法如下：计算飞行器离地速度
Figure BDA0003249762450000055
其中Gtakeoff为飞行器器起飞重量。根据该方法求得的离地速度，通过发动机推力模型，求得发动机离地推力TLD0。进而更新得到
Figure BDA0003249762450000056
其中αLD为离地迎角。若VLD0与VLD1满足误差范围，则VLD1即为飞行器离地速度，若VLD0与VLD1，不满足误差范围，则令VLD0＝VLD1，通过发动机推力模型，迭代求解，直至VLD0与VLD1满足误差范围，求得飞行器离地速度。
S44.计算起飞耗油量。步骤如下：发动机初始推力
Figure BDA0003249762450000057
其中A为停机角。根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器初始速度、发动机初始推力，求得起飞初始耗油率sfcqf0。根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器离地速度、发动机离地推力，求得起飞离地耗油率sfcqf。起飞时间
Figure BDA0003249762450000061
其中g为当地重力加速度，f为跑道摩擦系数，
Figure BDA0003249762450000062
为离地升阻比。起飞滑跑耗油量
Figure BDA0003249762450000063
所以起飞耗油量Wqf＝W0+Whp。
S5：根据场高等条件，计算爬升至任务剖面规定巡航高度所需爬升高度，估算爬升过程耗油量Wps，分为加速爬升和等速爬升两个阶段，分别对应加速爬升耗油量和等速爬升耗油量；具体步骤为：
S51.确定起飞爬升初始重量Gps＝Gtakeoff-Wqf，确定爬升开始高度VI和爬升终点高度VF。爬升高度Hps＝HI-HF，其中根据爬升初始高度HI和爬升终点高度HF。一般的，爬升开始速度为飞行器起飞离地速度，爬升终点高度为飞行器巡航速度，爬升开始高度为机场高度，爬升终点高度为飞行器出航巡航高度。
S52.根据爬升高度和爬升起止速度，将爬升段平均分为N段进行计算。根据每段爬升的平均高度，修正该段高度的大气密度、声速。计算辅助变量，
Figure BDA0003249762450000064
Figure BDA0003249762450000065
S53.迭代计算每段爬升的飞行器速度、马赫数、爬升速度、飞行器重量、耗油量。由
Figure BDA0003249762450000066
计算该爬升段飞行器升力系数，并通过飞机极曲线求得对应阻力系数。由马赫数、该爬升段平均高度和发动机工作状态，代入S1所述发动机推力模型，得到发动机推力Ti。计算动能变化修正系数
Figure BDA0003249762450000067
爬升速度
Figure BDA0003249762450000068
该段爬升时间
Figure BDA0003249762450000069
根据该段发动机设定的工作状态、飞行速度、发动机推力和平均高度，计算燃油消耗率sfci。该段燃油消耗量Wi＝sfci*DT。该爬升段水平距离Li＝VicosθiDT，其中
Figure BDA00032497624500000610
更新飞行器速度、马赫数、爬升速度、飞行器重量，进行下一段爬升耗油量计算。
S54.累加每一段爬升的耗油量，得爬升段耗油量
Figure BDA0003249762450000071
额外的，如爬升段有发动机加力工作情况，需按照发动机加力工作状态，则需要计算给定加力时长、飞行器高度、飞行器推力的发动机耗油量，累加到爬升段总耗油量中。
S6：根据任务剖面规定，结合发动机工作状态、任务要求推力、飞行高度和飞行速度，结合S1所述发动机燃油消耗率模型，结合巡逻、作战等任务段时间，计算任务段耗油量；
S7：确定飞行器安全着陆余油，根据前述计算的起飞耗油量、爬升过程耗油量、任务段耗油量，确定巡航可用余油Wcruise，并以巡航可用余油的一半确定出航可用油量Wsetout；由巡航段所需推力、设定巡航速度、设定巡航高度推算巡航耗油率。具体步骤如下：
S71.将计算段油量分为N个计算段，得到辅助变量
Figure BDA0003249762450000072
确定出航巡航开始前飞行器重量GSI，根据出航巡航段高度，修正出航巡航段马赫数Ss、大气密度ρ、大气压力P，求得动压q＝0.7PMa2，其中
Figure BDA0003249762450000073
V为飞行器飞行速度。
S72.确定出航巡航段的计算段发动机设定的工作状态、飞行器升力系数
Figure BDA0003249762450000074
由飞行器极曲线得到阻力系数CXi，飞行器升阻比
Figure BDA0003249762450000075
巡航所需推力
Figure BDA0003249762450000076
其中IT为发动机推力损失。
S73.根据该巡航段发动机设定的工作状态、飞行速度、发动机推力和平均高度，代入发动机燃油消耗率模型，计算燃油消耗率sfci，计算段航程
Figure BDA0003249762450000077
更新飞行器重量，进行下一计算段巡航耗油量计算。
S74.累加每一段的航程，得到出航巡航距离
Figure BDA0003249762450000078
特别的，计算段数量N的确定由飞行器设计航程确定，一般航程越远，所设置N的值越大。当飞行器有额外挂载时，第2步的飞行器阻力系数需要增加飞行器挂载的阻力系数增量。第3步，由计算段耗油量和耗油率，累加可以计算得到航时，同时也可以根据计算段航程与耗油量，得到单位航程耗油量，进行远航速度计算。
S8：根据具体任务需要，确定飞行器返航构型、返航巡航高度，修正飞机返航重量、马赫数、大气密度、大气压力和升阻力系数，结合S1所述发动机燃油消耗率模型，计算方法同S7所述出航巡航计算，最终得到返航巡航距离
Figure BDA0003249762450000079
S9：出航巡航距离与出航巡航距离比较，若误差较大，则适当增加或减少出航油量，迭代计算。若误差合理，则得出作战半径R＝Lps+max{Lsc,Lrc}。
最终得出飞行半径剖面如图2所示。
本发明提供的柔性飞行器飞行半径计算方法，通过对任务剖面不同阶段的计算，基于柔性飞行器飞行动力学理论，对飞行器模型、巡航段耗油率进行修正，进而完成作战半径的计算。本发明中所述的柔性飞行器飞行半径计算方法，无需专业性能软件或仿真软件所用的大量原始数据，计算快捷简单，可实现大展弦比无人机等飞行器的作战半径估算，为解决其总体性能论证、制定作战方案点等问题提供高可靠性数据支撑，具有重要的工程应用价值。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (8)

1.一种柔性飞行器飞行半径确定方法，其特征在于，包括以下步骤：
S1：建立发动机燃油消耗率与平均高度、平均速度、发动机推力、发动机设定工作状态的多元线性回归模型，即为发动机燃油消耗率模型；
发动机推力与高度、飞行马赫数、发动机设定工作状态的线性回归模型，即为发动机推力模型；
S2：根据高度，建立温度、声速、压力、大气密度的数学模型；
S3：根据建立考虑柔性影响修正的升力系数模型，根据公式
Figure FDA0003249762440000011
其中ηi为i阶弹性模态，
Figure FDA0003249762440000012
为弹性模态对应的升力增量；确定飞行器升力系数，并修正飞行器极曲线的气动数据，用于校正后续计算过程；
S4：根据已知飞行器起飞重量、停机角、起飞机场高度、温度、离地迎角和跑道摩擦系数，计算起飞耗油量；
S5：根据场高等条件，计算爬升至任务剖面规定巡航高度所需爬升高度，估算爬升过程耗油量；
S6：根据任务剖面规定，结合发动机工作状态、任务要求推力、飞行高度和飞行速度，结合S1所述发动机燃油消耗率模型，结合巡逻、作战等任务段时间，计算任务段耗油量；
S7：出航巡航计算；确定飞行器安全着陆余油，根据前述计算的起飞耗油量、爬升过程耗油量、任务段耗油量，确定巡航可用余油，并以巡航可用余油的一半确定出航可用油量；由巡航段所需推力、设定巡航速度、设定巡航高度推算巡航耗油率；
S8：根据具体任务需要，确定飞行器返航构型、返航巡航高度，修正飞机返航重量、马赫数、大气密度、大气压力和升阻力系数，结合S1所述发动机燃油消耗率模型，计算方法同S7所述出航巡航计算，最终得到返航巡航距离；
S9：出航巡航距离与返航巡航距离比较，若误差较大，则适当增加或减少出航油量，迭代计算；若误差合理，则爬升水平距离与出/返航距离累加得作战半径。
2.根据权利要求1所述的一种柔性飞行器飞行半径确定方法，其特征在于，步骤S4中，所述的起飞耗油量包括起飞前暖机耗油量和起飞滑跑耗油量。
3.根据权利要求2所述的一种柔性飞行器飞行半径确定方法，其特征在于，步骤S4的具体步骤为：
S41.根据起飞暖机时间和发动机暖机耗油率，计算起飞暖机耗油；根据起飞机场高度，代入大气密度的数学模型，修正大气密度；
S42.根据设置的离地迎角和飞行器的升力系数曲线，求得离地升力系数，进而由飞行器极曲线求得离地阻力系数；
S43.迭代计算飞行器离地速度与发动机离地推力；方法如下：计算飞行器离地速度；通过发动机推力模型，求得发动机离地推力；进而更新得到新的离地速度；若两次计算的离地速度满足误差范围，结果为飞行器离地速度，若不满足误差范围，则用新的离地速度通过发动机推力模型，迭代求解，直至满足误差范围；
S44.根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器初始速度、发动机初始推力，求得起飞初始耗油率；根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器离地速度、发动机离地推力，求得起飞离地耗油率；计算起飞时间
Figure FDA0003249762440000021
其中g为当地重力加速度，f为跑道摩擦系数，
Figure FDA0003249762440000022
为离地升阻比，Gtakeoff为飞行器起飞重量，T0为初始推力，Tqf为起飞离地推力，VLD为离地速度；取初始耗油率和起飞离地耗油率的平均值作为起飞过程耗油率，结合起飞时间，计算得到起飞耗油量。
4.根据权利要求1所述的一种柔性飞行器飞行半径确定方法，其特征在于，步骤S5中，估算爬升过程耗油量，分为加速爬升和等速爬升两个阶段，分别对应加速爬升耗油量和等速爬升耗油量。
5.根据权利要求4所述的一种柔性飞行器飞行半径确定方法，其特征在于，步骤S5中具体步骤为：
S51.确定起飞爬升初始重量、爬升开始高度和爬升终点高度；爬升开始速度为飞行器起飞离地速度，爬升终点高度为飞行器巡航速度，爬升开始高度为机场高度，爬升终点高度为飞行器出航巡航高度；
S52.根据爬升高度和爬升起止速度，将爬升段平均分为N段进行计算；根据每段爬升的平均高度，代入大气密度与声速数学模型，修正该段高度的大气密度、声速；
S53.迭代计算每段爬升的飞行器速度、马赫数、爬升速度、飞行器重量、耗油量；计算该爬升段飞行器升力系数，并通过飞机极曲线求得对应阻力系数；由马赫数、该爬升段平均高度和发动机工作状态，代入所述发动机推力模型，得到发动机推力；根据该段发动机设定的工作状态、飞行速度、发动机推力和平均高度，计算燃油消耗率；计算动能变化修正系数
Figure FDA0003249762440000023
爬升速度
Figure FDA0003249762440000024
该段爬升时间
Figure FDA0003249762440000025
该爬升段燃油消耗量Wi＝sfci*DT，该爬升段水平距离，其中Vi为计算段飞行速度，DV为计算段速度变化量，DH为计算段高度变化量，Gps为爬升段重量，Ti为爬升段推力，Di为爬升段阻力，sfci为计算段燃油消耗率；更新飞行器速度、马赫数、爬升速度、飞行器重量，进行下一段爬升耗油量计算；
S54.累加每一段爬升的耗油量，得爬升段耗油量，累加每一段爬升的水平距离，得到爬升段水平距离。
6.根据权利要求5所述的一种柔性飞行器飞行半径确定方法，其特征在于，步骤S5中，如爬升段有发动机加力工作情况，需按照发动机加力工作状态，则需要计算给定加力时长、飞行器高度、飞行器推力的发动机耗油量，累加到爬升段总耗油量中。
7.根据权利要求1所述的一种柔性飞行器飞行半径确定方法，其特征在于，步骤S7中，具体步骤如下：
S71.将计算段油量分为N个计算段，确定出航巡航开始前飞行器重量，根据出航巡航段高度，修正出航巡航段马赫数、大气密度、大气压力，求得动压；
S72.确定出航巡航段的计算段发动机设定的工作状态、飞行器升力系数，由飞行器极曲线得到阻力系数、飞行器升阻比；结合发动机推力损失计算巡航所需推力；
S73.根据该巡航段发动机设定的工作状态、飞行速度、发动机推力和平均高度，计算燃油消耗率，计算段航程
Figure FDA0003249762440000031
其中sfci为计算段燃油消耗率，DW为计算段燃油消耗量，V为巡航速度；更新飞行器重量，进行下一计算段巡航耗油量计算；
S74.累加每一段的航程，得到出航巡航距离。
8.根据权利要求7所述的一种柔性飞行器飞行半径确定方法，其特征在于，步骤S5中，计算段数量N的确定由飞行器设计航程确定，航程越远，所设置N的值越大；当飞行器有额外挂载时，第S72步的飞行器阻力系数需要增加飞行器挂载的阻力系数增量；第S73步，由计算段耗油量和耗油率，累加计算得到航时，同时根据计算段航程与耗油量，得到单位航程耗油量，进行远航速度计算。
CN202111042267.3A 2021-09-07 2021-09-07 一种柔性飞行器飞行半径确定方法 Active CN113589847B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111042267.3A CN113589847B (zh) 2021-09-07 2021-09-07 一种柔性飞行器飞行半径确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111042267.3A CN113589847B (zh) 2021-09-07 2021-09-07 一种柔性飞行器飞行半径确定方法

Publications (2)

Publication Number Publication Date
CN113589847A true CN113589847A (zh) 2021-11-02
CN113589847B CN113589847B (zh) 2022-03-01

Family

ID=78241461

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111042267.3A Active CN113589847B (zh) 2021-09-07 2021-09-07 一种柔性飞行器飞行半径确定方法

Country Status (1)

Country Link
CN (1) CN113589847B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114065399A (zh) * 2021-11-23 2022-02-18 北京航空航天大学 一种考虑复杂气象条件下的无人飞行器飞行性能计算方法
CN114065398A (zh) * 2021-11-23 2022-02-18 北京航空航天大学 一种大展弦比柔性飞行器飞行性能计算方法
CN114217634A (zh) * 2021-12-06 2022-03-22 中国人民解放军国防科技大学 一种基于物联网的无人机蓝牙通信系统及方法
CN117807916A (zh) * 2024-02-29 2024-04-02 中国航空工业集团公司西安飞机设计研究所 一种喷气式飞机给定推力最大航程巡航设计方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160257420A1 (en) * 2015-03-04 2016-09-08 The Boeing Company Method for optimum maximum range cruise speed in an aircraft
CN110794866A (zh) * 2019-10-17 2020-02-14 成都飞机工业（集团）有限责任公司 一种集爬升-巡航-下降为整体的航时性能优化方法
CN112179350A (zh) * 2020-09-29 2021-01-05 中国航空工业集团公司沈阳飞机设计研究所 一种基于高效排序算法的巡航解算方法
CN112254734A (zh) * 2020-09-29 2021-01-22 中国航空工业集团公司沈阳飞机设计研究所 一种基于高效排序算法的经济导航解算方法
CN112287459A (zh) * 2020-10-30 2021-01-29 中国航发沈阳发动机研究所 一种飞机发动机燃油效率计算方法
CN112528407A (zh) * 2020-12-10 2021-03-19 中国航空工业集团公司沈阳飞机设计研究所 一种固定翼飞机亚音速巡航航程优化设计方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160257420A1 (en) * 2015-03-04 2016-09-08 The Boeing Company Method for optimum maximum range cruise speed in an aircraft
CN110794866A (zh) * 2019-10-17 2020-02-14 成都飞机工业（集团）有限责任公司 一种集爬升-巡航-下降为整体的航时性能优化方法
CN112179350A (zh) * 2020-09-29 2021-01-05 中国航空工业集团公司沈阳飞机设计研究所 一种基于高效排序算法的巡航解算方法
CN112254734A (zh) * 2020-09-29 2021-01-22 中国航空工业集团公司沈阳飞机设计研究所 一种基于高效排序算法的经济导航解算方法
CN112287459A (zh) * 2020-10-30 2021-01-29 中国航发沈阳发动机研究所 一种飞机发动机燃油效率计算方法
CN112528407A (zh) * 2020-12-10 2021-03-19 中国航空工业集团公司沈阳飞机设计研究所 一种固定翼飞机亚音速巡航航程优化设计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LI YI 等: "Optimization of Range and Endurance of a Propeller UAV Based on SQP Algorithm", 《PROCEEDINGS OF 2017 IEEE INTERNATIONAL CONFERENCE ON UNMANNED SYSTEMS (ICUS)》 *
张云飞 等: "螺旋桨飞机等高度飞行的航程与航时计算", 《飞行力学》 *
惠辉辉 等: "基于飞机模型的性能计算与管理", 《飞机设计》 *
李亮明: "基于参考重量的直升机航程性能试飞研究", 《工程与试验》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114065399A (zh) * 2021-11-23 2022-02-18 北京航空航天大学 一种考虑复杂气象条件下的无人飞行器飞行性能计算方法
CN114065398A (zh) * 2021-11-23 2022-02-18 北京航空航天大学 一种大展弦比柔性飞行器飞行性能计算方法
CN114065398B (zh) * 2021-11-23 2022-08-16 北京航空航天大学 一种大展弦比柔性飞行器飞行性能计算方法
CN114065399B (zh) * 2021-11-23 2024-04-16 北京航空航天大学 一种考虑复杂气象条件下的无人飞行器飞行性能计算方法
CN114217634A (zh) * 2021-12-06 2022-03-22 中国人民解放军国防科技大学 一种基于物联网的无人机蓝牙通信系统及方法
CN114217634B (zh) * 2021-12-06 2023-10-10 中国人民解放军国防科技大学 一种基于物联网的无人机蓝牙通信系统及方法
CN117807916A (zh) * 2024-02-29 2024-04-02 中国航空工业集团公司西安飞机设计研究所 一种喷气式飞机给定推力最大航程巡航设计方法
CN117807916B (zh) * 2024-02-29 2024-05-17 中国航空工业集团公司西安飞机设计研究所 一种喷气式飞机给定推力最大航程巡航设计方法

Also Published As

Publication number Publication date
CN113589847B (zh) 2022-03-01

Similar Documents

Publication Publication Date Title
CN113589847B (zh) 一种柔性飞行器飞行半径确定方法
CN114065398B (zh) 一种大展弦比柔性飞行器飞行性能计算方法
CA2940044C (en) Flight path optimization using nonlinear programming
CN114065399B (zh) 一种考虑复杂气象条件下的无人飞行器飞行性能计算方法
CN101441478B (zh) 一种小型无人机自动着陆拉平控制方法及其装置
CN105701552B (zh) 一种飞行航路垂直剖面的确定方法
EP2691297B1 (en) Method and appartus for operating an aircraft propeller
CN106768123A (zh) 一种无人直升机燃油预估方法
CN105279290A (zh) 一种四发螺旋桨飞机续航性能计算方法
EP3598261A1 (en) Method and system for determining a descent profile
CN114329766A (zh) 面向深度强化学习的飞行动力学模型可信度评估方法
CN114004021B (zh) 用于飞行管理系统性能管理的巡航燃油流量计算方法
Hansen et al. Induced moment effects of formation flight using two F/A-18 aircraft
Donovan High speed problems of aircraft and experimental methods
CN113094938B (zh) 一种面向海上搜救任务仿真的直升机油耗模型构建方法
CN113625768B (zh) 一种火星直升机航迹规划方法、系统、设备及存储介质
US4445179A (en) Aircraft minimum drag speed system
CN112035947B (zh) 一种带整体油箱的机翼剖面载荷计算方法
Labonté Formulas for the fuel of climbing propeller driven airplanes
CN112441253A (zh) 一种多旋翼无人机实时尾迹的仿真分析系统及方法
CN115856359B (zh) 一种直升机空速在线修正方法
CN115344068B (zh) 一种飞行模拟机快速稳定方法、电子设备及介质
Bridges Flight performance envelope for an aircraft with a fixed-pitch propeller
Howard et al. Flight test of a half-scale unmanned air vehicle
Papež et al. Experimental method for unmanned aerial vehicles

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant