CN114065398A - 一种大展弦比柔性飞行器飞行性能计算方法 - Google Patents
一种大展弦比柔性飞行器飞行性能计算方法 Download PDFInfo
- Publication number
- CN114065398A CN114065398A CN202111395385.2A CN202111395385A CN114065398A CN 114065398 A CN114065398 A CN 114065398A CN 202111395385 A CN202111395385 A CN 202111395385A CN 114065398 A CN114065398 A CN 114065398A
- Authority
- CN
- China
- Prior art keywords
- aircraft
- calculating
- speed
- climbing
- takeoff
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
本发明公开了一种大展弦比柔性飞行器飞行性能计算方法,结合平均高度、平均速度、飞行时间、油箱平均油量以及发动机工作状态的燃油消耗率的多元回归模型,建立考虑柔性飞行动力学修正的无人飞行器飞行力学模型,并在此基础上建立飞行性能计算方法,具体飞行性能主要包括:最大爬升率、最大盘旋角速度、最小平飞速度、最大平飞速度、起飞性能、着陆性能和最大航程计算。该方法无需专业性能软件与仿真软件所用的大量原始数据,对飞行性能的计算高效快捷,可以提高大展弦比柔性飞行器在给定任务情况的飞行性能计算的准确性。
Description
技术领域
本发明涉及航空技术领域,尤其涉及一种大展弦比柔性飞行器飞行性能计算方法。
背景技术
飞行性能评估作为无人机飞行任务规划的重要环节,精确的飞行性能计算可以保证飞行任务剖面的正确性和时效性。同时在飞行器设计之初由于缺乏完整详尽的计算分析数据,通过飞行性能评估,可以对飞行器起飞、着陆、飞行器的飞行时间、航程、机动、巡逻等参量进行评估,进而规划制定合适的飞行任务。
目前,在上述起飞性能、最大航程和着陆性能等飞行性能计算缺少具体系统的方法,对于具有一定柔性的大展弦比飞行器(如运输机、大展弦比无人机)而言,以此为基础的飞行性能计算往往产生较大误差,缺少一个具备计算精度的飞行器性能模型。因此本发明提供了一种满足一定精度要求,同时计算快捷简单的飞行性能计算方法,该方法无需专业性能软件与仿真软件所用的大量原始数据,对飞行性能的计算高效快捷,可以提高大展弦比柔性飞行器在给定任务情况的飞行性能计算的准确性。
发明内容
有鉴于此,本发明提供了一种大展弦比柔性飞行器飞行性能计算方法,包括飞行器起飞性能、最大爬升率、最大盘旋角速度、着陆性能和最大航程6个部分,可分为如下步骤:
S1:建立发动机燃油消耗率与平均高度、平均速度、发动机推力、发动机设定工作状态的多元线性回归模型;发动机推力与高度、飞行马赫数、发动机设定工作状态的线性回归模型与发动机推力回归模型。
一般情况下,发动机工作状态可以分为慢车状态、加力状态和中间推力状态。数据可由实际飞行试验测得,也可由发动机模拟试验测得,一般为离散的数据点,需要进行多元线性回归。
S2:建立环境模型,根据气压高度,建立温度、声速、压力、空气密度和大展弦比柔性影响的数学模型。大展弦比柔性影响主要通过考虑模态修正升力系数计算,根据建立考虑大展弦比柔性条件下的升力系数模型,根据公式其中ηi为i阶弹性模态,为弹性模态对应的升力系数增量,需要通过已知飞行器机翼翼型与尺寸进行结构模态计算。确定并计算飞行器各项升力系数,修正飞行器气动数据,用于校正后续计算过程。
S3:进行起飞性能计算,计算起飞滑跑距离以及相应的起飞需要时间以及消耗燃料量。
具体步骤为:
S31.已知飞行器起飞重量、停机角、起飞机场高度、温度、离地迎角和跑道摩擦系数。根据起飞暖机时间和发动机暖机耗油率,计算起飞暖机耗油;根据起飞机场高度,修正大气密度。
S32.根据设置的离地迎角和飞行器的升力系数曲线,求得离地升力系数,进而由极曲线求得阻力系数,求离地阻力系数。
S33.根据大展弦比柔性飞行器升力系数模型和起飞重量求得的离地速度,通过发动机推力回归模型,求得发动机离地推力。进而更新得到离地速度。迭代求解,直至满足误差范围,得到飞行器离地速度,并推得发动机离地推力。
S34.按照公式根据发动机初始推力、起飞离地推力和起飞离地耗油率计算起飞滑跑距离以及相应的起飞需要时间以及消耗燃料量。起飞时间起飞滑跑距离其中g为当地重力加速度,f为跑道摩擦系数,为离地升阻比。起飞滑跑耗油量所以起飞耗油量Wqf=W0+Whp。
S4:计算飞行器给定的飞行高度和环境条件下的最大爬升率。根据场高等条件,计算爬升至任务剖面规定巡航高度所需爬升高度,估算爬升过程耗油量,分为加速爬升和等速爬升两个阶段,分别对应加速爬升耗油量和等速爬升耗油量,具体步骤为:
S41.根据已知高度条件,修正空气密度、声速。确定起飞爬升初始重量,迭代计算飞行器速度和爬升轨迹角。判断飞行器速度是否满足最大表速限制,满足情况下进行下一步计算。
S42.根据升力系数曲线,确定最大升力系数。确定起飞爬升初始重量Gps,爬升升力系数CYp,判断爬升升力系数与最大升力系数关系。如果CYp>CYmax,则要增大飞行器速度,返回第一步进行迭代计算;如果CYp<CYmax,进行下一步迭代计算。
S43.根据升力系数和极曲线求得对应阻力系数,得到飞行器阻力系数,进而得到飞行器阻力。由马赫数、该爬升段平均高度和发动机工作状态,代入S1所述发动机推力模型,得到发动机推力并计算切向过载。如果切向过载小于1,则增大飞行器速度,重复上述迭代计算。如果切向过载大于等于1,则进行下一步计算。
S44.计算最大爬升速度,并更新临时上升轨迹角,如果临时上升轨迹角与该迭代步飞行器角的差值满足一定范围,则此时最大爬升率与上升轨迹角即为该状态下的飞行器的最大爬升率与上升轨迹角。如果不满足误差范围,则需要增大迭代飞行器速度,重复上述步骤。
S5:计算最小平飞速度,即飞机在某高度条件下能够保持等速水平直线飞行的最小速度。根据任务剖面规定,确定飞行器重量与飞行高度。根据高度修正空气密度与大气压强。根据公式求得最小平飞速度,最小平飞马赫数。
S6:计算最大盘旋角速度,计算以飞机强度或升力特性所允许的最大过载作非定常盘旋所能达到的最大盘旋角速度和盘旋速度。具体步骤如下:
S61.根据已知盘旋高度和气象条件,修正空气密度、声速。确定飞行器盘旋重量,迭代飞行器速度,计算该速度下飞行器最大升力系数CYmax。
S62.确定盘旋升力系数,以及对应的升力线斜率CYA。判断升力系数与最大升力系数关系。如果CYp>CYmax,则要增大迭代飞行器速度Vi,返回第一步进行迭代计算;如果CYp<CYmax,进行下一步迭代计算。
S63.更新飞行器攻角,由马赫数、盘旋高度和发动机工作状态,代入S1所述发动机推力模型,得到发动机推力,计入推力分量的盘旋升力系数。计算飞行器的法向过载ny,判断过载是否符合要求,包括结构过载限制和超音速限制,如果ny>1,且满足过载限制,进行下一步迭代计算。如果过载不满足要求,则增加迭代飞行器速度,返回第一步。
S7:计算飞行器着陆性能,主要包括着陆总距离、着陆滑跑距离以及相应的着陆需要时间和消耗燃料量的计算。具体步骤如下:
S71.已知飞行器着陆重量、给定接地迎角、机场高度、温度、离地迎角和刹车阻力系数(一般取0.2-0.3)。根据起飞机场高度,修正大气密度ρ。
S72.根据设置的接地迎角和飞行器的升力系数曲线,求得接地升力系数,进而由极曲线求得阻力系数。
S73.迭代计算飞行器离地速度。计算飞行器接地速度Vjd。根据该方法求得的离地速度。预设安全系数Q(一般取1.2-1.3之间),求得飞行器进场速度Vjc。
S74.按照公式计算着陆距离和对应的着陆时间。估算着陆空中段距离其中为接地空中段升阻比Kjd和接地滑跑段升阻比Kjc的平均值。着陆空中段时间着陆滑跑距离其中g为当地重力加速度,为跑道摩擦系数,f2为刹车阻力系数,Kjd,为接地时刻升阻比。估算着陆滑跑时间计算着陆滑跑距离Lzl=Lzh+Lzk,着陆滑跑时间tzl=tzh+tzk。
S8:计算飞行器转场航程,包含起飞段、起飞后加速上升段、等速上升段、等高等速巡航段、下滑段叠加组成。主要步骤包括:
S81.起飞段油量计算。根据起飞暖机时间和发动机暖机耗油率,计算起飞暖机耗油;根据起飞机场高度,修正大气密度。根据设置的离地迎角和飞行器的升力系数曲线,求得离地升力系数,进而由极曲线求得离地阻力系数。计算飞行器离地速度。根据该方法求得的离地速度,通过发动机推力回归模型,求得发动机离地推力。进而更新得到离地速度。若两次计算离地速度满足误差范围,则此时计算的离地速度即为飞行器离地速度,若不满足误差范围,则将更新得到的离地速度代入发动机推力回归模型,迭代求解,直至满足误差范围,求得飞行器离地速度。根据发动机初始推力、起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器初始速度,求得起飞初始耗油率。根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器离地速度、发动机离地推力,求得起飞离地耗油率。进而依据公式计算起飞与起飞过程耗油量。
S82.起飞爬升段油量与航程计算。确定起飞爬升初始重量,确定爬升高度。一般的,爬升开始速度为飞行器起飞离地速度,爬升终点高度为飞行器巡航速度,爬升开始高度为机场高度,爬升终点高度为飞行器出航巡航高度。根据爬升高度和爬升起止速度,将爬升段平均分为N段进行计算。根据每段爬升的平均高度,修正该段高度的空气密度、声速。迭代计算每段爬升的飞行器速度、马赫数、爬升速度、飞行器重量、耗油量。累加每一段爬升的耗油量,得爬升段耗油量以及爬升段水平距离。
S83.等速巡航段油量与航程计算。将计算段油量分为N个计算段,确定巡航开始前飞行器重量,根据巡航段高度,修正巡航段马赫数、大气密度和大气压力。确定巡航段的计算段发动机设定的工作状态、飞行器升力系数,由飞行器极曲线得到阻力系数,计算得到飞行器升阻比。依据巡航段飞行器重量和升阻比,结合发动机推力损失,计算巡航所需推力。根据该巡航段发动机设定的工作状态、飞行速度、发动机推力和平均高度,计算燃油消耗率与计算段航程。更新飞行器重量,进行下一计算段巡航耗油量计算。累加每一段的航程,得到出航巡航距离。
S84.下滑段航程设定与转场航程计算。按照不同高度设置不同的下滑段航程,下滑段一般按照设定下滑程序进行,一般为给定的高度和下滑段航。累计爬升段、巡航段和下滑段水平距离,得到转场航程;累计各段时间,得到转场任务时间。
本发明提供的柔性飞行器飞行半径计算方法,通过对任务剖面不同阶段的计算,基于柔性飞行器飞行动力学理论,对飞行器模型、巡航段耗油率进行修正,进而完成作战半径的计算。本发明中所述的柔性飞行器飞行半径计算方法,无需专业性能软件或仿真软件所用的大量原始数据,计算快捷简单,可实现大展弦比无人机等飞行器的作战半径估算,为解决其总体性能论证、制定作战方案点等问题提供高可靠性数据支撑,具有重要的工程应用价值。
附图说明
图1为本发明实施例中的大展弦比固定翼飞行器起飞性能计算流程图;
图2为本发明实施例中的大展弦比固定翼飞行器爬升性能计算流程图;
图3为本发明实施例中的大展弦比固定翼飞行器最小平飞速度计算流程图;
图4为本发明实施例中的大展弦比固定翼飞行器盘旋性能计算流程图;
图5为本发明实施例中的大展弦比固定翼飞行器着陆性能计算流程图;
图6为本发明实施例中的大展弦比固定翼飞行器转场航程计算流程图。
具体实施方式
下面将结合本发明实施方式中的附图,对本发明实施方式中的技术方案进行清楚、完整的描述,显然,所描述的实施方式仅仅是作为例示,并非用于限制本发明。
下面通过一个具体的实施例对本发明提供的上述柔性飞行器飞行性能计算方法的具体实施进行详细说明。
步骤分别如图所1示,包括如下步骤:
S1:建立发动机燃油消耗率与平均高度、平均速度、发动机推力、发动机设定工作状态的多元线性回归模型;发动机推力与高度、飞行马赫数、发动机设定工作状态的线性回归模型。S1所述燃油消耗率回归模型可表示为:
Sfc=f(H,V,T,state)
发动机推力回归模型可表示为:
T=T(H,Ma,state)
其中H为发动机/飞行器平均高度,V为飞行器平均速度,Ma为飞行器飞行马赫数,T为发动机推力,state为发动机工作状态,一般情况下,发动机工作状态可以分为慢车状态、加力状态和中间推力状态。数据可由实际飞行试验测得,也可由发动机模拟试验测得,一般为离散的数据点,需要进行多元线性回归。
S2:建立环境模型,根据气压高度,建立温度、声速、压力、空气密度和大展弦比柔性影响的数学模型,分别如下:
温度:
声速:
压力:
空气密度:
大展弦比柔性影响:
根据建立考虑大展弦比柔性条件下的升力系数模型,根据公式其中ηi为i阶弹性模态,为弹性模态对应的升力系数增量,需要通过已知飞行器机翼翼型与尺寸进行结构模态计算。确定并计算飞行器各项升力系数,修正飞行器气动数据,用于校正后续计算过程。
S3:进行起飞性能计算,计算起飞滑跑距离以及相应的起飞需要时间以及消耗燃料量。
流程如图1所示,具体步骤为:
S31.已知飞行器起飞重量、停机角、起飞机场高度、温度、离地迎角和跑道摩擦系数。根据起飞暖机时间T0和发动机暖机耗油率sfc0,计算起飞暖机耗油W0=T0*sfc0;根据起飞机场高度,修正大气密度ρ。
S32.根据设置的离地迎角和结合前述大展弦比柔性飞行器飞行器的升力系数曲线,求得离地升力系数CYLD,进而由极曲线求得阻力系数CXL。
S33.迭代计算飞行器离地速度VLD与发动机离地推力TLD。计算飞行器离地速度其中Gtakeoff为飞行器器起飞重量。根据该方法求得的离地速度,通过发动机推力回归模型,求得发动机离地推力TLD0。进而更新得到其中αLD为离地迎角。若VLD0与VLD1满足误差范围,则VLD1即为飞行器离地速度,若VLD0与VLD1,不满足误差范围,则令VLD0=VLD1,通过发动机推力回归模型,迭代求解,直至VLD0与VLD1满足误差范围,求得飞行器离地速度。
S34.计算起飞滑跑距离以及相应的起飞需要时间以及消耗燃料量。发动机初始推力其中A为停机角。起飞时间起飞滑跑距离其中g为当地重力加速度,f为跑道摩擦系数,为离地升阻比。根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器初始速度、发动机初始推力,求得起飞初始耗油率sfcqf0。根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器离地速度、发动机离地推力,求得起飞离地耗油率sfcqf。起飞滑跑耗油量所以起飞耗油量Wqf=W0+Whp。
本计算例中,得到计算结果为:起飞距离448.7m,起飞离地速度105.28m/s,起飞时间25.71s。
S4:飞行器爬升性能计算,计算流程如图2所示。计算飞行器给定的飞行高度和环境条件下的最大爬升率。根据场高等条件,计算爬升至任务剖面规定巡航高度所需爬升高度,估算爬升过程耗油量Wps,分为加速爬升和等速爬升两个阶段,分别对应加速爬升耗油量和等速爬升耗油量,具体步骤为:
S42.根据升力系数曲线,确定最大升力系数CYmax。确定起飞爬升初始重量Gps,爬升升力系数判断爬升升力系数与最大升力系数关系。如果CYp>CYmax,则要增大迭代飞行器速度Vi,返回第一步进行迭代计算;如果CYp<CYmax,进行下一步迭代计算。
S43.根据升力系数和极曲线求得对应阻力系数CX,计算飞行器阻力系数CXp=CX+CX0,进而阻力由马赫数、该爬升段平均高度和发动机工作状态,代入S1所述发动机推力模型,得到发动机推力T。计算切向过载如果Nx<1,则增大迭代飞行器速度Vi,重复上述计算。如果Nx≥1,则进行下一步计算。
S44.计算最大爬升速度Vymax=NxV,更新临时上升轨迹角如果θi与θi+1的差值满足一定范围,则此时最大爬升率与上升轨迹角即为该状态下的飞行器的最大爬升率与上升轨迹角。如果不满足误差范围,则需要增大迭代飞行器速度Vi,重复上述步骤。
在本实施例中,飞行器在高度5000m下最大爬升率为87.67m/s,上升轨迹角为18.62°。
S5:计算最小平飞速度,即飞机在某高度条件下能够保持等速水平直线飞行的最小速度,具体计算流程如图3所示。根据任务剖面规定,确定飞行器重量与飞行高度。根据高度修正空气密度与大气压强。根据公式求得最小平飞速度最小平飞马赫数其中G为飞行器重量,Ry为,ρ为空气密度,P为大气压强S为飞行器升力面积,Cymax为飞行器最大升力系数。
本实施例中得到计算结果:13000米高度下,飞行器最小平飞速度为307km/h。
S6:计算最大盘旋角速度,计算以飞机强度或升力特性所允许的最大过载作非定常盘旋所能达到的最大盘旋角速度和盘旋速度,流程如图4所示。具体步骤如下:
S61.根据已知盘旋高度和气象条件,修正空气密度、声速。确定飞行器盘旋重量Gturn,迭代飞行器速度Vi,计算该速度下飞行器最大升力系数CYmax。
S62.确定盘旋升力系数以及对应的升力线斜率CYA。判断升力系数与最大升力系数关系。如果CYp>CYmax,则要增大迭代飞行器速度Vi,返回第一步进行迭代计算;如果CYp<CYmax,进行下一步迭代计算。
S63.更新飞行器攻角由马赫数、盘旋高度和发动机工作状态,代入S1所述发动机推力模型,得到发动机推力T,计入推力分量的盘旋升力系数为飞行器的法向过载判断过载是否符合要求,包括结构过载限制和超音速限制,如果ny>1,且满足过载限制,进行下一步迭代计算。如果过载不满足要求,则增加迭代飞行器速度,返回第一步。
本实施例中计算结果为在高度10000m情况下,飞行器最大盘旋角速度为0.147rad/s,盘旋速度573.97km/h,盘旋半径1084.6m。
S7:计算飞行器着陆性能,流程如图5所示。主要包括着陆总距离、着陆滑跑距离以及相应的着陆需要时间和消耗燃料量的计算。具体步骤如下:
S71.已知飞行器着陆重量、给定接地迎角、机场高度、温度、离地迎角和刹车阻力系数(一般取0.2-0.3)。根据起飞机场高度,修正大气密度ρ。
S72.根据设置的接地迎角和根据大展弦比柔性飞行器的升力系数,求得接地升力系数CYjd,进而由极曲线求得阻力系数CXL,求得接地阻力系数CXjd=CXL+ΔCD。
S73.迭代计算飞行器离地速度VLD。计算飞行器接地速度其中Glanding为飞行器着陆重量。根据该方法求得的离地速度。预设安全系数Q(一般取1.2-1.3之间),求得飞行器进场速度Vjc=VjkQ。
S74.计算着陆距离和对应的着陆时间。估算着陆空中段距离其中为接地空中段升阻比Kjd和接地滑跑段升阻比Kjc的平均值。着陆空中段时间着陆滑跑距离其中g为当地重力加速度,为跑道摩擦系数,f2为刹车阻力系数,为接地时刻升阻比。估算着陆滑跑时间计算着陆距离Lzl=Lzh+Lzk,着陆时间tzl=tzh+tzk。
本实施例中计算结果为:进场速度为210km/s,着陆迎角8°情况下的着陆滑跑距离为657.4m,着陆总距离为1819.17m,滑跑时间为13.1min,着陆总时间为14.7min。
S8:计算飞行器转场航程,流程如图6所示。转场航程包含起飞段、起飞后加速上升段、等速上升段、等高等速巡航段、下滑段叠加组成。主要步骤包括:
S81.起飞段油量计算。根据起飞暖机时间T0和发动机暖机耗油率sfc0,计算起飞暖机耗油W0=T0*sfc0;根据起飞机场高度,修正大气密度ρ。根据设置的离地迎角和飞行器的升力系数曲线,求得离地升力系数CYLD,进而由极曲线求得离地阻力系数CXLD。计算飞行器离地速度其中Gtakeoff为飞行器器起飞重量。根据该方法求得的离地速度,通过发动机推力回归模型,求得发动机离地推力TLD0。进而更新得到其中αLD为离地迎角。若VLD0与VLD1满足误差范围,则VLD1即为飞行器离地速度,若VLD0与VLD1,不满足误差范围,则令VLD0=VLD1,通过发动机推力回归模型,迭代求解,直至VLD0与VLD1满足误差范围,求得飞行器离地速度。发动机初始推力其中A为停机角。根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器初始速度、发动机初始推力,求得起飞初始耗油率sfcqf0。根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器离地速度、发动机离地推力,求得起飞离地耗油率sfcqf。起飞时间其中g为当地重力加速度,f为跑道摩擦系数,为离地升阻比。起飞滑跑耗油量所以起飞耗油量Wqf=W0+Whp。
S82.起飞爬升段油量与航程计算。确定起飞爬升初始重量Gps=Gtakeoff-Wqf,确定爬升开始高度VI和爬升终点高度VF。爬升高度Hps=HI-HF,其中根据爬升初始高度HI和爬升终点高度HF。一般的,爬升开始速度为飞行器起飞离地速度,爬升终点高度为飞行器巡航速度,爬升开始高度为机场高度,爬升终点高度为飞行器出航巡航高度。根据爬升高度和爬升起止速度,将爬升段平均分为N段进行计算。根据每段爬升的平均高度,修正该段高度的空气密度、声速。计算辅助变量,迭代计算每段爬升的飞行器速度、马赫数、爬升速度、飞行器重量、耗油量。由计算该爬升段飞行器升力系数,并通过飞机极曲线求得对应阻力系数。由马赫数、该爬升段平均高度和发动机工作状态,代入S1所述发动机推力模型,得到发动机推力Ti。计算动能变化修正系数爬升速度该段爬升时间根据该段发动机设定的工作状态、飞行速度、发动机推力和平均高度,计算燃油消耗率sfci。该段燃油消耗量Wi=sfci*DT。该爬升段水平距离Li=VicosθiDT,其中更新飞行器速度、马赫数、爬升速度、飞行器重量,进行下一段爬升耗油量计算。累加每一段爬升的耗油量,得爬升段耗油量爬升段水平距离
S83.等速巡航段油量与航程计算。将巡航段油量分为N个计算段,计算辅助变量确定巡航开始前飞行器重量GSI,根据巡航段高度,修正巡航段马赫数Ss、大气密度ρ、大气压力P,求得动压q=0.7PMa2,其中V为飞行器飞行速度。确定巡航段的计算段发动机设定的工作状态、飞行器升力系数由飞行器极曲线得到阻力系数CXi,飞行器升阻比巡航所需推力其中IT为发动机推力损失。根据该巡航段发动机设定的工作状态、飞行速度、发动机推力和平均高度,计算燃油消耗率sfci,计算段航程更新飞行器重量,进行下一计算段巡航耗油量计算。累加每一段的航程,得到巡航距离
S84.下滑段航程设定与转场航程计算。按照不同高度设置不同的下滑段航程Ldescent,下滑段一般按照设定下滑程序进行,本实施例按照13000m高度,设定下滑段航程20.5km进行计算。转场航程L=Lps+Lcruise+Ldescent。
本实施例转场航程计算结果为:在巡航高度为10000m,巡航马赫数0.6的情况下,转场航程2162.51km,时间为3.79h。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (10)
1.一种大展弦比柔性飞行器飞行性能计算方法,其特征在于,包括以下步骤:
S1:建立发动机燃油消耗率与平均高度、平均速度、发动机推力、发动机设定工作状态的多元线性回归模型;发动机推力与高度、飞行马赫数、发动机设定工作状态的线性回归模型与发动机推力回归模型;
S2:建立环境模型,根据气压高度,建立温度、声速、压力、空气密度和大展弦比柔性影响的数学模型;
大展弦比柔性影响通过考虑模态修正升力系数计算,根据建立考虑大展弦比柔性条件下的升力系数模型,根据公式其中ηi为i阶弹性模态,为弹性模态对应的升力系数增量,需要通过已知飞行器机翼翼型与尺寸进行结构模态计算;
确定并计算飞行器各项升力系数,修正飞行器气动数据,用于校正后续计算过程;
S3:进行起飞性能计算,计算起飞滑跑距离以及相应的起飞需要时间以及消耗燃料量;
S4:计算飞行器给定的飞行高度和环境条件下的最大爬升率;根据场高等条件,计算爬升至任务剖面规定巡航高度所需爬升高度,估算爬升过程耗油量,分为加速爬升和等速爬升两个阶段,分别对应加速爬升耗油量和等速爬升耗油量;
S5:计算最小平飞速度,即飞机在某高度条件下能够保持等速水平直线飞行的最小速度;根据任务剖面规定,确定飞行器重量与飞行高度;根据高度修正空气密度与大气压强;根据公式求得最小平飞速度,最小平飞马赫数;
S6:计算最大盘旋角速度,计算以飞机强度或升力特性所允许的最大过载作非定常盘旋所能达到的最大盘旋角速度和盘旋速度;
S7:计算飞行器着陆性能,主要包括着陆总距离、着陆滑跑距离以及相应的着陆需要时间和消耗燃料量的计算;
S8:计算飞行器转场航程,包含起飞段、起飞后加速上升段、等速上升段、等高等速巡航段、下滑段叠加组成。
2.根据权利要求1所述的一种大展弦比柔性飞行器飞行性能计算方法,其特征在于,S3具体包括以下子步骤:
S31.已知飞行器起飞重量、停机角、起飞机场高度、温度、离地迎角和跑道摩擦系数;根据起飞暖机时间和发动机暖机耗油率,计算起飞暖机耗油;根据起飞机场高度,修正大气密度;
S32.根据设置的离地迎角和飞行器的升力系数曲线,求得离地升力系数,进而由极曲线求得阻力系数,求离地阻力系数;
S33.根据大展弦比柔性飞行器升力系数模型和起飞重量求得的离地速度,通过发动机推力回归模型,求得发动机离地推力;进而更新得到离地速度;迭代求解,直至满足误差范围,得到飞行器离地速度,并推得发动机离地推力;
3.根据权利要求1所述的一种大展弦比柔性飞行器飞行性能计算方法,其特征在于,S4具体包括以下子步骤:
S41.根据已知高度条件,修正空气密度、声速;确定起飞爬升初始重量,迭代计算飞行器速度和爬升轨迹角;判断飞行器速度是否满足最大表速限制,满足情况下进行下一步计算;
S42.根据升力系数曲线,确定最大升力系数;确定起飞爬升初始重量Gps,爬升升力系数CYp,判断爬升升力系数与最大升力系数关系;如果CYp>CYmax,则要增大飞行器速度,返回第一步进行迭代计算;如果CYp<CYmax,进行下一步迭代计算;
S43.根据升力系数和极曲线求得对应阻力系数,得到飞行器阻力系数,进而得到飞行器阻力;由马赫数、该爬升段平均高度和发动机工作状态,代入S1所述发动机推力模型,得到发动机推力并计算切向过载;如果切向过载小于1,则增大飞行器速度,重复上述迭代计算;如果切向过载大于等于1,则进行下一步计算;
S44.计算最大爬升速度,并更新临时上升轨迹角,如果临时上升轨迹角与该迭代步飞行器角的差值满足一定范围,则此时最大爬升率与上升轨迹角即为该状态下的飞行器的最大爬升率与上升轨迹角;如果不满足误差范围,则需要增大迭代飞行器速度,重复上述步骤。
4.根据权利要求1所述的一种大展弦比柔性飞行器飞行性能计算方法,其特征在于,S6具体包括以下子步骤:
S61.根据已知盘旋高度和气象条件,修正空气密度、声速;确定飞行器盘旋重量,迭代飞行器速度,计算该速度下飞行器最大升力系数CYmax;
S62.确定盘旋升力系数,以及对应的升力线斜率CYA;判断升力系数与最大升力系数关系;如果CYp>CYmax,则要增大迭代飞行器速度Vi,返回第一步进行迭代计算;如果CYp<CYmax,进行下一步迭代计算;
S63.更新飞行器攻角,由马赫数、盘旋高度和发动机工作状态,代入S1所述发动机推力模型,得到发动机推力,计入推力分量的盘旋升力系数;计算飞行器的法向过载ny,判断过载是否符合要求,包括结构过载限制和超音速限制,如果ny>1,且满足过载限制,进行下一步迭代计算;如果过载不满足要求,则增加迭代飞行器速度,返回第一步;
5.根据权利要求1所述的一种大展弦比柔性飞行器飞行性能计算方法,其特征在于,S7具体包括以下子步骤:
S71.已知飞行器着陆重量、给定接地迎角、机场高度、温度、离地迎角和刹车阻力系数;根据起飞机场高度,修正大气密度ρ;
S72.根据设置的接地迎角和飞行器的升力系数曲线,求得接地升力系数,进而由极曲线求得阻力系数;
S73.迭代计算飞行器离地速度;计算飞行器接地速度Vjd;根据该方法求得的离地速度;预设安全系数Q,求得飞行器进场速度Vjc;
6.根据权利要求1所述的一种大展弦比柔性飞行器飞行性能计算方法,其特征在于,S8具体包括以下子步骤:
S81.起飞段油量计算;
S82.起飞爬升段油量与航程计算;
S83.等速巡航段油量与航程计算;
S84.下滑段航程设定与转场航程计算。
7.根据权利要求6所述的一种大展弦比柔性飞行器飞行性能计算方法,其特征在于,S81.起飞段油量计算,包括:
根据起飞暖机时间和发动机暖机耗油率,计算起飞暖机耗油;根据起飞机场高度,修正大气密度;根据设置的离地迎角和飞行器的升力系数曲线,求得离地升力系数,进而由极曲线求得离地阻力系数;计算飞行器离地速度;根据该方法求得的离地速度,通过发动机推力回归模型,求得发动机离地推力;进而更新得到离地速度;若两次计算离地速度满足误差范围,则此时计算的离地速度即为飞行器离地速度,若不满足误差范围,则将更新得到的离地速度代入发动机推力回归模型,迭代求解,直至满足误差范围,求得飞行器离地速度;根据发动机初始推力、起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器初始速度,求得起飞初始耗油率;根据起飞滑跑过程发动机设定的工作状态、起飞机场高度、飞行器离地速度、发动机离地推力,求得起飞离地耗油率;进而依据公式计算起飞与起飞过程耗油量。
8.根据权利要求6所述的一种大展弦比柔性飞行器飞行性能计算方法,其特征在于,S82.起飞爬升段油量与航程计算,包括:
确定起飞爬升初始重量,确定爬升高度;爬升开始速度为飞行器起飞离地速度,爬升终点高度为飞行器巡航速度,爬升开始高度为机场高度,爬升终点高度为飞行器出航巡航高度;根据爬升高度和爬升起止速度,将爬升段平均分为N段进行计算;根据每段爬升的平均高度,修正该段高度的空气密度、声速;迭代计算每段爬升的飞行器速度、马赫数、爬升速度、飞行器重量、耗油量;累加每一段爬升的耗油量,得爬升段耗油量以及爬升段水平距离。
9.根据权利要求6所述的一种大展弦比柔性飞行器飞行性能计算方法,其特征在于,S83.等速巡航段油量与航程计算,包括:
将计算段油量分为N个计算段,确定巡航开始前飞行器重量,根据巡航段高度,修正巡航段马赫数、大气密度和大气压力;确定巡航段的计算段发动机设定的工作状态、飞行器升力系数,由飞行器极曲线得到阻力系数,计算得到飞行器升阻比;依据巡航段飞行器重量和升阻比,结合发动机推力损失,计算巡航所需推力;根据该巡航段发动机设定的工作状态、飞行速度、发动机推力和平均高度,计算燃油消耗率与计算段航程;更新飞行器重量,进行下一计算段巡航耗油量计算;累加每一段的航程,得到出航巡航距离。
10.根据权利要求6所述的一种大展弦比柔性飞行器飞行性能计算方法,其特征在于,S84.下滑段航程设定与转场航程计算,包括:按照不同高度设置不同的下滑段航程,下滑段按照设定下滑程序进行,为给定的高度和下滑段航;累计爬升段、巡航段和下滑段水平距离,得到转场航程;累计各段时间,得到转场任务时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111395385.2A CN114065398B (zh) | 2021-11-23 | 2021-11-23 | 一种大展弦比柔性飞行器飞行性能计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111395385.2A CN114065398B (zh) | 2021-11-23 | 2021-11-23 | 一种大展弦比柔性飞行器飞行性能计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114065398A true CN114065398A (zh) | 2022-02-18 |
CN114065398B CN114065398B (zh) | 2022-08-16 |
Family
ID=80279685
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111395385.2A Active CN114065398B (zh) | 2021-11-23 | 2021-11-23 | 一种大展弦比柔性飞行器飞行性能计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114065398B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114638129A (zh) * | 2022-05-23 | 2022-06-17 | 中国航空工业集团公司西安飞机设计研究所 | 一种基于发动机状态变化的爬升能力数据获得方法 |
CN114647994A (zh) * | 2022-05-24 | 2022-06-21 | 中国航空工业集团公司西安飞机设计研究所 | 一种爬升性能快速处理方法 |
CN115593647A (zh) * | 2022-11-03 | 2023-01-13 | 清华大学(Cn) | 用于垂直起降飞机的串联混合动力系统航程最优设计方法 |
CN117807916A (zh) * | 2024-02-29 | 2024-04-02 | 中国航空工业集团公司西安飞机设计研究所 | 一种喷气式飞机给定推力最大航程巡航设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120158220A1 (en) * | 2010-12-20 | 2012-06-21 | Selex Systemi Integrati S.P.A. | Fast vertical trajectory prediction method for air traffic management, and relevant atm system |
CN112528407A (zh) * | 2020-12-10 | 2021-03-19 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种固定翼飞机亚音速巡航航程优化设计方法 |
CN112613122A (zh) * | 2020-12-24 | 2021-04-06 | 中国航空工业集团公司西安飞机设计研究所 | 一种层流控制飞机顶层参数设计方法 |
CN113589847A (zh) * | 2021-09-07 | 2021-11-02 | 北京航空航天大学 | 一种柔性飞行器飞行半径确定方法 |
-
2021
- 2021-11-23 CN CN202111395385.2A patent/CN114065398B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120158220A1 (en) * | 2010-12-20 | 2012-06-21 | Selex Systemi Integrati S.P.A. | Fast vertical trajectory prediction method for air traffic management, and relevant atm system |
CN112528407A (zh) * | 2020-12-10 | 2021-03-19 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种固定翼飞机亚音速巡航航程优化设计方法 |
CN112613122A (zh) * | 2020-12-24 | 2021-04-06 | 中国航空工业集团公司西安飞机设计研究所 | 一种层流控制飞机顶层参数设计方法 |
CN113589847A (zh) * | 2021-09-07 | 2021-11-02 | 北京航空航天大学 | 一种柔性飞行器飞行半径确定方法 |
Non-Patent Citations (4)
Title |
---|
YUN CHENG 等: "Efficiency evaluation of energy harvesting from discrete gust response for plate wing", 《JOURNAL OF MECHANICAL ENGINEERING SCIENCE》 * |
向锦武 等: "长航时无人机关键技术研究进展", 《哈尔滨工业大学学报》 * |
王敏 等: "机种机型运动识别方法的研究", 《计算机与数字工程》 * |
石文 等: "高空超长航时太阳能无人机主要技术问题分析", 《飞航导弹》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114638129A (zh) * | 2022-05-23 | 2022-06-17 | 中国航空工业集团公司西安飞机设计研究所 | 一种基于发动机状态变化的爬升能力数据获得方法 |
CN114647994A (zh) * | 2022-05-24 | 2022-06-21 | 中国航空工业集团公司西安飞机设计研究所 | 一种爬升性能快速处理方法 |
CN115593647A (zh) * | 2022-11-03 | 2023-01-13 | 清华大学(Cn) | 用于垂直起降飞机的串联混合动力系统航程最优设计方法 |
CN117807916A (zh) * | 2024-02-29 | 2024-04-02 | 中国航空工业集团公司西安飞机设计研究所 | 一种喷气式飞机给定推力最大航程巡航设计方法 |
CN117807916B (zh) * | 2024-02-29 | 2024-05-17 | 中国航空工业集团公司西安飞机设计研究所 | 一种喷气式飞机给定推力最大航程巡航设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114065398B (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114065398B (zh) | 一种大展弦比柔性飞行器飞行性能计算方法 | |
CN114065399B (zh) | 一种考虑复杂气象条件下的无人飞行器飞行性能计算方法 | |
CN113589847B (zh) | 一种柔性飞行器飞行半径确定方法 | |
CN109446582B (zh) | 一种考虑地球自转的高精度降阶平稳滑翔动力学建模方法 | |
CN110598294B (zh) | 常规布局飞机飞行轨迹稳定性确定方法及装置 | |
CN111717411A (zh) | 一种基于试飞数据标准重量下巡航推力增量的修正方法 | |
CN114004021B (zh) | 用于飞行管理系统性能管理的巡航燃油流量计算方法 | |
Hansen et al. | Induced moment effects of formation flight using two F/A-18 aircraft | |
Gong et al. | Using flight manual data to derive aero-propulsive models for predicting aircraft trajectories | |
CN113468828A (zh) | 一种飞机空中飞行颠簸强度指数计算方法 | |
CN116203600B (zh) | 一种无人机通信信号丢失后带动力运动轨迹循迹的方法 | |
CN108398883B (zh) | 一种rlv进场着陆轨迹快速推演及确定方法 | |
Tian et al. | Wake encounter simulation and flight validation with UAV close formation flight | |
CN113094938B (zh) | 一种面向海上搜救任务仿真的直升机油耗模型构建方法 | |
CN114912284A (zh) | 基于第一法则的飞行管理系统爬升性能预测算法 | |
CN109117584A (zh) | 一种低速飞机突风载荷系数计算方法及设备 | |
Tian et al. | UAS-based wind estimation using sinusoidal gust model | |
CN111859540B (zh) | 一种大气扰动中飞机颠簸响应的计算方法 | |
Smith et al. | Effect of competing V-tail models on a UAS 6-DOF nonlinear simulation | |
Lan et al. | Analysis of the crash of a transport aircraft and assessment of fuzzy-logic stall recovery | |
Curry et al. | An in-flight investigation of ground effect on a forward-swept wing airplane | |
CN116700358B (zh) | 一种无人机转弯阶段非线性定高补偿控制方法 | |
CN117421825B (zh) | 一种近地层风环境影响下的大型飞机cfd模拟方法及系统 | |
Wang et al. | Automatic Landing Control Design of Gyroplane | |
Howard et al. | Flight test of a half-scale unmanned air vehicle |
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 |