CN106648673B - 一种连续下降运行程序的分析和设计方法 - Google Patents

一种连续下降运行程序的分析和设计方法 Download PDF

Info

Publication number
CN106648673B
CN106648673B CN201611233019.6A CN201611233019A CN106648673B CN 106648673 B CN106648673 B CN 106648673B CN 201611233019 A CN201611233019 A CN 201611233019A CN 106648673 B CN106648673 B CN 106648673B
Authority
CN
China
Prior art keywords
aircraft
continuous descending
operation program
model
determining
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.)
Active
Application number
CN201611233019.6A
Other languages
English (en)
Other versions
CN106648673A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201611233019.6A priority Critical patent/CN106648673B/zh
Publication of CN106648673A publication Critical patent/CN106648673A/zh
Application granted granted Critical
Publication of CN106648673B publication Critical patent/CN106648673B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F8/00—Arrangements for software engineering
    • G06F8/20—Software design
    • G06F8/22—Procedural
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F11/00—Error detection; Error correction; Monitoring
    • G06F11/36—Preventing errors by testing or debugging software
    • G06F11/3604—Software analysis for verifying properties of programs
    • G06F11/3608—Software analysis for verifying properties of programs using formal methods, e.g. model checking, abstract interpretation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Quality & Reliability (AREA)
  • Traffic Control Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种连续下降运行程序的分析和设计方法，根据初始设计的连续下降运行程序，联合导航点坐标与数据库编码表，建立航空器意图模型；选择典型机型，考虑航空器质量以及风速风向的不确定性，基于航空器质点模型，结合上述航空器意图模型，采用航空器基本性能数据库，生成航空器四维航迹；基于上述生成的一系列航空器四维航迹，分析不确定条件下，各类航空器过航路点高度与速度分布，以及不同机型配比时，前后机的间隔分布，判断CDO程序初始设计是否恰当，现行管制移交间隔是否合理的结论，从而指导CDO程序的设计与修改。本发明克服了现有CDO程序存在的缺点，设计出在节能、减排、降噪方面具备一定优势的CDO程序。

Description

一种连续下降运行程序的分析和设计方法
技术领域
本发明属于民航技术领域，特别涉及了一种连续下降运行程序的分析和设计方法。
背景技术
近年来，随着能源危机的加剧，环境保护意识的加强，燃油消耗、机场噪声、废气排放等问题日益凸显。如何在安全运行的前提下，尽可能地降低燃油消耗、缓解机场噪声、减少废气排放等已成为民航业关注的焦点。连续下降运行程序(Continuous DescentOperation,CDO)正是基于上述背景，率先在航空发达国家提出并应用。目前，美国、荷兰、新加坡等国家的CDO技术已经相对成熟，并且已经在大部分机场运行，带来了巨大的经济与环境效益。
我国在《国务院关于促进民航业发展的若干意见》(国发〔2012〕24号)中明确指出——“到2020年，我国初步形成安全、便捷、高效、绿色的现代化民用航空体系”、“要切实打造绿色低碳航空”。同时，中国民航局在近期发布的多个规划性文件中，也着重强调要加速行业节能减排工作。这些诉求必将有力地推进我国民航CDO程序的研究与应用工作。
我国在CDO程序方面的研究刚刚起步，与上述航空发达国家的差距较大。同时，CDO程序的设计与应用，均离不开CDO程序高度/速度设计的分析以及管制移交间隔的分析，四维航迹生成与预测方面的研究能够为此类分析提供基础。目前四维航迹预测的方法，主要有全能量方程法和质点模型法，均采用是正向计算生成四维航迹。然而，CDO程序的关键问题是如何确定TOD点的位置，并在TOD点之后生成连续下降剖面。上述采用正向计算的航迹预测方法都不能解决该问题。同时，目前对于飞行程序的分析主要从安全性、节能减排、经济性等方面入手，鲜考虑飞行过程中不确定性因素与管制移交习惯等方面，从而不能够模拟出更适合终端空域实际运行情况的飞行程序。
发明内容
为了解决上述背景技术提出的技术问题，本发明旨在提供一种连续下降运行程序的分析和设计方法，克服现有CDO程序存在的缺点，设计在节能、减排、降噪方面具备一定优势的CDO程序。
为了实现上述技术目的，本发明的技术方案为：
一种连续下降运行程序的分析和设计方法，包括以下步骤：
(1)根据初始设计的连续下降运行程序，联合导航点坐标与数据库编码表，考虑水平航迹与垂直剖面，建立航空器意图模型；
(2)选择典型机型，考虑航空器质量和运行环境的不确定性，结合步骤(1)建立的航空器意图模型，建立航空器质点模型，生成航空器四维航迹；
(3)基于生成的航空器四维航迹，分析不确定条件下各类航空器过航路点高度与速度分布，以及不同机型配比时前后机的间隔分布，从而判断连续下降运行程序设计是否合理，以及现行管制移交间隔是否合理；
(4)根据步骤(3)的判断结论，修改连续下降运行程序，更新连续下降运行程序的水平航迹、高度/速度限制以及管制移交间隔，并将修改后的连续下降运行程序返回步骤(1)，重新分析，直至连续下降运行程序安全合理。
进一步地，步骤(1)的具体步骤如下：
(A)根据初始设计的连续下降运行程序，确定连续下降运行程序中包含的航路点；
(B)根据导航点坐标和数据库编码表，确定航空器所要经过的航路点的经纬度、航向以及是否为转弯点，建立连续下降运行程序的水平航迹；
(C)根据导航点坐标和数据库编码表，确定航空器所要经过的航路点的速度限制和高度限制信息，建立连续下降运行程序的垂直剖面；
(D)根据连续下降运行程序的水平航迹与垂直剖面，以及连续下降运行程序的运行特点，建立航空器意图模型。
进一步地，步骤(2)的具体步骤如下：
(a)基于航空器基本性能数据库，确定航空器的性能参数；
(b)确定运行期间的环境模型，包括确定随高度变化的温度、压力与空气密度参数，以及确定随高度与位置变化的风场数据；
(c)确定航空器的质量参数，以及质量参数的不确定性；
(d)基于步骤(1)建立的航空器意图模型，建立航空器运动的质点模型；
(e)基于航空器运动质点模型，生成航空器四维航迹。
进一步地，在步骤(b)中，确定随高度变化的温度、压力与空气密度参数的步骤如下：
(a1)根据气压高度，确定温度T：
T＝T0+ΔT+βT·Hp
上式中，T0＝288.15K，表示国际标准大气条件下在平均海平面处的温度；ΔT表示温度偏差；Hp表示气压高度；βT＝-0.0065K/m，表示温度垂直递减率；
(b1)根据温度T，确定压力p：
Figure BDA0001194819190000031
上式中，p0＝101325Pa，表示国际标准大气条件下的空气压力；g0＝9.80665m/s2，表示重力加速度；R＝287.05287m2/(K·s2)，表示空气常数。
(c1)根据温度T与压力p，确定空气密度r：
Figure BDA0001194819190000032
进一步地，在步骤(b)中，确定随高度与位置变化的风场数据的步骤如下：
(a2)在天气预报中心上选取所需的各气压层上的风场数据，包含日期、时刻、风分量和高度层；
(b2)根据连续下降运行程序的水平航迹的范围，裁剪区域，由纬度范围确定南北区域、经度范围确定东西区域；
(c2)根据精度要求确定单位网格跨度；
(d2)根据设定下载风场数据；
(e2)针对高度层的u风分量与v风分量，计算出风速和风向：
Figure BDA0001194819190000041
Figure BDA0001194819190000042
上式中，Vwind为风速，
Figure BDA0001194819190000043
为风向。
进一步地，在步骤(d)中，建立航空器运动的质点模型的步骤如下：
(a3)根据航空器意图模型中连续下降运行程序的水平航迹，分为直线航段与转弯航段，进而构建水平方向上的质点运动模型，其中直线航段：
Figure BDA0001194819190000044
式中，s为飞行距离；VGS为地速；VTAS为真空速；
Figure BDA0001194819190000045
为风角，
Figure BDA0001194819190000046
MC为航线角；
Figure BDA0001194819190000047
为偏流，
Figure BDA0001194819190000048
其中转弯航段：
Figure BDA0001194819190000049
上式中，ROT为转弯率；φ为转弯坡度；
(b3)根据航空器意图模型中连续下降运行程序的垂直剖面，建立高度与速度变化的质点方程：
Figure BDA00011948191900000410
上式中，h为高度；γ为航径角；D为航空器阻力；Thr为航空器推力；m为航空器质量；g为重力加速度；
(c3)考虑航空器在飞行过程中的燃油消耗，建立航空器质量变化方程。
进一步地，在步骤(c)中，确定航空器质量参数以及质量参数不确定性的步骤如下：
(a4)根据性能参数，以对应机型的参考质量作为该航空器的质量参数；
(b4)设执行连续下降运行程序的航空器质量符合正态分布，选择合适的均值与方差构建航空器质量参数的不确定性模型。
进一步地，在步骤(e)中，生成航空器四维航迹的步骤如下：
(a5)采用4阶龙格库塔法对航空器运动的质点模型进行反向积分，生成航空器的垂直剖面；
(b5)对航空器所要经过的航路点采用等距离散的方式，生成航空器的水平航迹；
(c5)将垂直剖面与水平航迹融合，得到航空器四维航迹。
进一步地，步骤(3)的具体过程如下：
确定执飞连续下降运行程序的主力机型，基于各机型执飞连续下降运行程序的四维航迹，以高度/速度剖面以及在各航路点高度/速度分布图的形式展示，分析下降顶点，分析是否满足国际民航组织提出的垂直剖面边界要求，从而判断连续下降运行程序设计是否合理；
确定执飞连续下降运行程序的前后机型配比情况，基于各机型配比执飞连续下降运行程序的四维航迹，依据现行的移交间隔，推算连续下降运行程序结束处的前后机间隔，分析不同机型配比时，现行的移交间隔是否满足安全间隔要求。
进一步地，在步骤(4)中，修改连续下降运行程序的过程如下：
如果航空器的过点高度、速度不满足国际民航组织提出的垂直剖面边界要求，则修改航路点的高度、速度限制；
如果连续下降运行程序结束处的前后机间隔不满足雷达管制的间隔要求，则增大管制移交间隔；如果满足要求，则适量减小管制移交间隔。
采用上述技术方案带来的有益效果：
(1)综合考虑影响因素，使得本发明具有准确性特点：
本发明考虑了航空器质量以及风速风向的不确定性，基于航空器质点模型，结合航空器意图模型，采用航空器基本性能数据库(BADA)，生成航空器四维航迹，因此确保了生成航迹的准确性与全面性。
(2)能够生成CDO运行下的航空器的四维航迹，使得本发明具有国内领先与国外新技术接轨特点：
本发明提出了一种利用4阶龙格库塔法(ODE45)可以反向计算的特点，结合特殊的航空器意图，反向生成连续下降的航空器高度剖面并确定TOD点位置的方法，使得本发明与正向生成航空器四维航迹的方法不同。
(3)利用生成的航空器四维航迹进行多角度的全面分析，然后根据分析结果进行程序设计，使得本发明具有可靠性特点：
本发明提出了分析各类航空器过航路点高度与速度分布、不同机型配比时前后机间隔分布，然后根据分析结果修改CDO程序的方法，使得本发明具有可靠性高的特点，同时对于管制移交间隔的分析，使得本发明更能贴近终端空域实际的运行情况。
(4)技术解决方案简单可靠，使得本发明便于应用：
本发明在设计各个模块时，通过深入研究欧美的四维航迹预测和飞行程序设计与分析的结构、功能，为满足实时性、可靠性的需求，采用了简单可靠的技术解决方案。
附图说明
图1为本发明的流程示意图；
图2为本发明中建立航空器意图模型流程示意图；
图3为本发明中生成航空器四维航迹流程示意图；
图4为本发明环境模型中确定随高度变化的温度、压力与空气密度参数流程示意图；
图5为本发明环境模型中确定随高度与位置变化的风场数据流程示意图；
图6为本发明建立航空器运动质点模型的流程示意图；
图7为本发明中航空器四维航迹生成算法流程示意图；
图8为典型的CDO程序下的航空器进场垂直飞行航迹示意图；
图9为本发明分析各类航空器过航路点高度与速度分布流程示意图；
图10为本发明分析不同机型配比时前后机间隔分布流程示意图；
图11为本发明修改CDO程序流程示意图。
具体实施方式
以下将结合附图，对本发明的技术方案进行详细说明。
本发明提出的一种连续下降运行程序的分析和设计方法，流程图如图1所示，包括以下步骤：
步骤11，根据初始设计的CDO程序，联合导航点坐标与数据库编码表，考虑水平航迹与垂直剖面，建立航空器意图模型；
步骤12，选择典型机型，考虑航空器质量以及风速风向的不确定性，基于航空器质点模型，结合上述航空器意图模型，采用航空器基本性能数据库(BADA)，生成航空器四维航迹；
步骤13，基于上述生成的一系列航空器四维航迹，分析不确定条件下，各类航空器过航路点高度与速度分布，以及不同机型配比时，前后机的间隔分布，获得：CDO程序初始设计是否恰当，现行管制移交间隔是否合理的结论。
步骤14，根据上述结论，指导CDO程序的设计与修改，更新CDO程序的水平航迹、高度/速度限制以及管制移交间隔，并重新分析，直至CDO程序安全合理。
图2为建立航空器意图模型流程示意图，具体包括以下步骤：
步骤21，根据初始设计的CDO程序，确定CDO程序中包含的航路点；
步骤22，根据导航点坐标和数据库编码表，确定航空器所要经过的航路点的经纬度、航向、是否转弯点等信息，建立CDO程序的水平航迹；
步骤23，根据导航点坐标和数据库编码表，确定航空器所要经过的航路点的速度限制、高度限制等信息，建立CDO程序的垂直剖面；
步骤24，根据上述CDO程序的水平航迹与垂直剖面，以及CDO程序的运行特点，建立航空器意图模型，包括：直线飞行、转弯飞行、等马赫数(Mach)下降、等校正空速(CAS,Calibrated Air Speed)下降、等下降率(ROD,Rate of Descent)减速下降等方式。本发明中CDO程序运行时的航空器意图与使用阶段对应关系在表1中记录。
表1
图3为生成航空器四维航迹流程示意图，具体包括以下步骤：
步骤31，基于欧控实验中心发布航空器基础资料(BADA)，确定航空器的性能参数，包括：航空器机型参数(含发动机数目、发动机类型、尾流等级)；质量参数(含最大/最小/参考质量，以及最大配载质量)；飞行包络参数(含最大飞行速度、最大运行高度等)；空气动力学参数(含机翼参考面积、航空器各类构型的失速速度、航空器各类构型的附加/诱导阻力系数)、发动机推力参数(含最大爬升/下降/进近/着陆推力系数)、燃油流量参数(含推力相关、下降与巡航的燃油流量系数)；
步骤32，确定运行期间的环境模型，包括：一方面确定随高度变化的温度、压力与空气密度参数；另一方面利用欧洲中期天气预报中心(ECMWF)给出的数据，确定随高度与位置变化的风场数据；
步骤33，确定航空器的质量参数，以及质量参数的不确定性；
步骤34，基于航空器的意图模型，建立航空器运动的质点模型；
步骤35，基于航空器运动的质点模型，确定航空器四维航迹生成的求解算法。
图4为环境模型中确定随高度变化的温度、压力与空气密度参数流程示图，具体包括以下步骤：
步骤41，根据气压高度(对流层顶以下)确定温度T：
T＝T0+ΔT+βT·Hp
上式中，T0＝288.15K，表示国际标准大气条件下在平均海平面处的温度；ΔT表示温度偏差；Hp表示气压高度；βT＝-0.0065K/m，表示温度垂直递减率；
步骤42，根据温度T确定压力p：
上式中，p0＝101325Pa，表示国际标准大气条件下的空气压力；g0＝9.80665m/s2，表示重力加速度；R＝287.05287m2/(K·s2)，表示空气常数；
步骤43，根据温度T与压力p确定空气密度：
Figure BDA0001194819190000101
图5为环境模型中确定随高度与位置变化的风场数据流程示意图，具体包括以下步骤：
步骤51，从欧洲中期天气预报中心(ECMWF)网站上选取所需的各气压层上风的数据，包含日期、时刻(0h 6h 12h 18h，UTC时间)、风分量、高度层等；
步骤52，根据初始设计CDO程序的水平航迹的范围，裁剪区域，分别由纬度范围确定南北区域，经度范围确定东西区域；
步骤53，确定单位网格跨度，数据可以定义在0.75°×0.75°经纬度跨度的网格上，也可根据所需精度自定义网格跨度；
步骤54，各选项确定后，下载文件，文件下载格式为GRIB格式；
步骤55，将所获取的以GRIB格式存储的风场信息，进行解码：首先，利用工具grib2ctl.exe生成整个文件的描述文件.ctl；然后，利用工具gribmap.exe生成映射文件.idx；最后，通过其描述文件，提取该文件存储记录的形式与结构信息，处理并保存数据(各高度层的u风分量与v风分量)；
步骤56，针对高度层的u风分量与v风分量，根据下述公式计算出风速和风向：
Figure BDA0001194819190000103
上式中，Vwind为风速，
Figure BDA0001194819190000104
为风向，南北方向，南风为正；东西方向，西风为正。
步骤57，将所得的风速、风向进行概率统计，分析其分布函数与系数，从而完成风场数据的估算。
确定航空器质量参数以及不确定性，具体包括以下步骤：
步骤61，根据航空器基础资料(BADA)提供的性能参数，确定以对应机型的参考质量作为该航空器的质量参数；
步骤62，假设执行CDO程序的航空器质量符合正态分布，选择合适的均值与方差构建航空器质量参数的不确定性模型。
图6为航空器运动的质点模型的建立流程示意图，具体包括以下步骤：
步骤71，根据建立的航空器意图模型中CDO程序的水平航迹，将其分为直线航段与转弯航段，进而构建水平方向上的质点运动模型，其中直线航段为：
Figure BDA0001194819190000111
上式中，s为飞行距离；VGS为地速；VTAS为真空速；
Figure BDA0001194819190000112
为风角，
Figure BDA0001194819190000113
MC为航线角；
Figure BDA0001194819190000114
为偏流，
Figure BDA0001194819190000115
其中转弯航段为：
Figure BDA0001194819190000116
上式中，ROT为转弯率；φ为转弯坡度；
步骤72，根据建立的航空器意图模型中的CDO程序的垂直剖面，建立高度与速度变化的质点方程：
Figure BDA0001194819190000117
Figure BDA0001194819190000118
上式中，γ为航径角；D为航空器阻力；Thr为航空器推力；m为航空器质量；g为重力加速度；
步骤73，计算航空器阻力，如下式：
上式中，CD为阻力系数，CD＝CD0+CD2·(CL)2，其中CL为升力系数；S为机翼参考面积；各系数参见航空器基础资料(BADA)；
步骤74，计算航空器推力，其最大起飞推力如下式：
Thrmax climb＝CTc,1·(1-h/CTc,2+CTc,3·h2)·(1-CTc,5·ΔT)
上式中，CTc,1、CTc,2、CTc,3和CTc,5均为推力系数，参见航空器基础资料(BADA)，且下降/进近/着陆的推力可视作最大爬升推力的函数，但与所处的高度以及飞行阶段相关；
步骤75，考虑航空器在飞行过程中的燃油消耗，建立航空器质量变化方程，如下式：
Figure BDA0001194819190000121
上式中，Cf1、Cf2、Cf3和Cf4均为燃油消耗系数，参见航空器基础资料(BADA)。
图7为航空器四维航迹生成中确定求解算法流程示意图，具体包括以下步骤：
步骤81，利用4阶龙格库塔法(ODE45)对航空器运动的质点模型进行反向积分，生成航空器的垂直剖面；图8为典型的CDO程序下的航空器进场垂直飞行航迹示意图；
步骤82，对航空器所要经过的航路点采用等距离散的方式，生成航空器的水平航迹；
步骤83，将垂直剖面与水平航迹融合，得到执飞CDO程序的四维航迹。
图9为分析各类航空器过航路点高度与速度分布流程示意图，具体包括以下步骤：
步骤91，确定执飞CDO程序的主力机型；
步骤92，确定仿真次数，以及确定CDO程序区域的风场信息，确定各机型的航空器质量分布；
步骤93，建立各机型执飞CDO程序的航空器意图；
步骤94，生成各机型执飞CDO程序的四维航迹；
步骤95，基于各机型执飞CDO程序的四维航迹，以高度/速度剖面以及其在各航路点高度/速度分布图的形式展示，分析下降顶点(TOD)的位置，分析是否满足国际民航组织(ICAO)提出的垂直剖面边界要求。
图10为分析不同机型配比时前后机间隔分布流程示意图，具体包括以下步骤：
步骤101，确定执飞CDO程序的前后机型配比情况；
步骤102，确定仿真次数，以及确定CDO程序区域的风场信息，确定各机型的航空器质量分布；
步骤103，建立各机型配比执飞CDO程序的航空器意图；
步骤104，生成各机型配比执飞CDO程序的四维航迹；
步骤105，基于各机型配比执飞CDO程序的四维航迹，依据现行的移交间隔，推算CDO程序结束处的前后机间隔，分析不同机型配比时，现行的移交间隔是否满足安全间隔要求。
图11为修改CDO程序流程示意图，具体包括以下步骤：
步骤111，根据分析各类航空器过航路点高度与速度分布，如果航空器的过点高度、速度不满足国际民航组织(ICAO)提出的垂直剖面边界要求，则修改航路点的高度、速度限制。限制主要分为4类：无限制、上界限制、下界限制与窗口限制；
步骤112，分析不同机型配比时前后机间隔分布，如果航空器对之间在CDO程序结束处，不满足雷达管制的间隔要求，则增大管制移交间隔；如果满足要求，则可以适当减小管制移交间隔；
步骤113，针对修改结果，重新构建航空器意图模型、生成航空器四维航迹，进行再次分析。
以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。

Claims (8)

1.一种连续下降运行程序的分析和设计方法，其特征在于，包括以下步骤：
(1)根据初始设计的连续下降运行程序，联合导航点坐标与数据库编码表，考虑水平航迹与垂直剖面，建立航空器意图模型；
(2)选择典型机型，考虑航空器质量和运行环境的不确定性，结合步骤(1)建立的航空器意图模型，建立航空器质点模型，生成航空器四维航迹；步骤(2)的具体过程如下：
(a)基于航空器基本性能数据库，确定航空器的性能参数；
(b)确定运行期间的环境模型，包括确定随高度变化的温度、压力与空气密度参数，以及确定随高度与位置变化的风场数据；
(c)确定航空器的质量参数，以及质量参数的不确定性；
(d)基于步骤(1)建立的航空器意图模型，建立航空器运动的质点模型；
(e)基于航空器运动质点模型，生成航空器四维航迹；所述生成航空器四维航迹的过程如下：
(a5)采用4阶龙格库塔法对航空器运动的质点模型进行反向积分，生成航空器的垂直剖面；
(b5)对航空器所要经过的航路点采用等距离散的方式，生成航空器的水平航迹；
(c5)将垂直剖面与水平航迹融合，得到航空器四维航迹；
(3)基于生成的航空器四维航迹，分析不确定条件下各类航空器过航路点高度与速度分布，以及不同机型配比时前后机的间隔分布，从而判断连续下降运行程序设计是否合理，以及现行管制移交间隔是否合理；
(4)根据步骤(3)的判断结论，修改连续下降运行程序，更新连续下降运行程序的水平航迹、高度/速度限制以及管制移交间隔，并将修改后的连续下降运行程序返回步骤(1)，重新分析，直至连续下降运行程序安全合理。
2.根据权利要求1所述一种连续下降运行程序的分析和设计方法，其特征在于：步骤(1)的具体步骤如下：
(A)根据初始设计的连续下降运行程序，确定连续下降运行程序中包含的航路点；
(B)根据导航点坐标和数据库编码表，确定航空器所要经过的航路点的经纬度、航向以及是否为转弯点，建立连续下降运行程序的水平航迹；
(C)根据导航点坐标和数据库编码表，确定航空器所要经过的航路点的速度限制和高度限制信息，建立连续下降运行程序的垂直剖面；
(D)根据连续下降运行程序的水平航迹与垂直剖面，以及连续下降运行程序的运行特点，建立航空器意图模型。
3.根据权利要求1所述一种连续下降运行程序的分析和设计方法，其特征在于：在步骤(b)中，确定随高度变化的温度、压力与空气密度参数的步骤如下：
(a1)根据气压高度，确定温度T：
T＝T0+ΔT+βT·Hp
上式中，T0＝288.15K，表示国际标准大气条件下在平均海平面处的温度；ΔT表示温度偏差；Hp表示气压高度；βT＝-0.0065K/m，表示温度垂直递减率；
(b1)根据温度T，确定压力p：
Figure FDA0002145015960000021
上式中，p0＝101325Pa，表示国际标准大气条件下的空气压力；g0＝9.80665m/s2，表示重力加速度；R＝287.05287m2/(K·s2)，表示空气常数；
(c1)根据温度T与压力p，确定空气密度r：
4.根据权利要求1所述一种连续下降运行程序的分析和设计方法，其特征在于：在步骤(b)中，确定随高度与位置变化的风场数据的步骤如下：
(a2)在天气预报中心上选取所需的各气压层上的风场数据，包含日期、时刻、风分量和高度层；
(b2)根据连续下降运行程序的水平航迹的范围，裁剪区域，由纬度范围确定南北区域、经度范围确定东西区域；
(c2)根据精度要求确定单位网格跨度；
(d2)根据设定下载风场数据；
(e2)针对高度层的u风分量与v风分量，计算出风速和风向：
Figure FDA0002145015960000032
Figure FDA0002145015960000033
上式中，Vwind为风速，
Figure FDA0002145015960000034
为风向。
5.根据权利要求4所述一种连续下降运行程序的分析和设计方法，其特征在于：在步骤(d)中，建立航空器运动的质点模型的步骤如下：
(a3)根据航空器意图模型中连续下降运行程序的水平航迹，分为直线航段与转弯航段，进而构建水平方向上的质点运动模型，其中直线航段：
Figure FDA0002145015960000035
式中，s为飞行距离；VGS为地速；VTAS为真空速；
Figure FDA0002145015960000036
为风角，
Figure FDA0002145015960000037
MC为航线角；
Figure FDA0002145015960000038
为偏流，
Figure FDA0002145015960000039
其中转弯航段：
Figure FDA0002145015960000041
上式中，ROT为转弯率；φ为转弯坡度；
(b3)根据航空器意图模型中连续下降运行程序的垂直剖面，建立高度与速度变化的质点方程：
Figure FDA0002145015960000042
Figure FDA0002145015960000043
上式中，h为高度；γ为航径角；D为航空器阻力；Thr为航空器推力；m为航空器质量；g为重力加速度；
(c3)考虑航空器在飞行过程中的燃油消耗，建立航空器质量变化方程。
6.根据权利要求1所述一种连续下降运行程序的分析和设计方法，其特征在于：在步骤(c)中，确定航空器质量参数以及质量参数不确定性的步骤如下：
(a4)根据性能参数，以对应机型的参考质量作为该航空器的质量参数；
(b4)设执行连续下降运行程序的航空器质量符合正态分布，选择合适的均值与方差构建航空器质量参数的不确定性模型。
7.根据权利要求1所述一种连续下降运行程序的分析和设计方法，其特征在于：步骤(3)的具体过程如下：
确定执飞连续下降运行程序的主力机型，基于各机型执飞连续下降运行程序的四维航迹，以高度/速度剖面以及在各航路点高度/速度分布图的形式展示，分析下降顶点，分析是否满足国际民航组织提出的垂直剖面边界要求，从而判断连续下降运行程序设计是否合理；
确定执飞连续下降运行程序的前后机型配比情况，基于各机型配比执飞连续下降运行程序的四维航迹，依据现行的移交间隔，推算连续下降运行程序结束处的前后机间隔，分析不同机型配比时，现行的移交间隔是否满足安全间隔要求。
8.根据权利要求1所述一种连续下降运行程序的分析和设计方法，其特征在于：在步骤(4)中，修改连续下降运行程序的过程如下：
如果航空器的过点高度、速度不满足国际民航组织提出的垂直剖面边界要求，则修改航路点的高度、速度限制；
如果连续下降运行程序结束处的前后机间隔不满足雷达管制的间隔要求，则增大管制移交间隔；如果满足要求，则适量减小管制移交间隔。
CN201611233019.6A 2016-12-28 2016-12-28 一种连续下降运行程序的分析和设计方法 Active CN106648673B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611233019.6A CN106648673B (zh) 2016-12-28 2016-12-28 一种连续下降运行程序的分析和设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611233019.6A CN106648673B (zh) 2016-12-28 2016-12-28 一种连续下降运行程序的分析和设计方法

Publications (2)

Publication Number Publication Date
CN106648673A CN106648673A (zh) 2017-05-10
CN106648673B true CN106648673B (zh) 2020-02-14

Family

ID=58832551

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611233019.6A Active CN106648673B (zh) 2016-12-28 2016-12-28 一种连续下降运行程序的分析和设计方法

Country Status (1)

Country Link
CN (1) CN106648673B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107085978B (zh) * 2017-06-21 2020-02-18 南京航空航天大学 一种基于所需到达时间的管制辅助决策指令生成方法
CN110415379A (zh) * 2019-07-31 2019-11-05 深圳市瑞达飞行科技有限公司 飞行数据处理方法、装置、计算机系统和可读存储介质
CN110765640B (zh) * 2019-11-05 2022-05-10 中国船舶重工集团海装风电股份有限公司 一种转子有效风速的计算方法、系统及设备
DE102020107456A1 (de) * 2020-03-18 2021-09-23 Volocopter Gmbh Verfahren und Steuergerät zur Kurvenkoordinierung eines Fluggerätes sowie ein Fluggerät mit Kurvenkoordinierung

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101692315A (zh) * 2009-09-25 2010-04-07 民航总局空管局技术中心 基于实时雷达数据的高精度4d飞机航迹分析方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101692315A (zh) * 2009-09-25 2010-04-07 民航总局空管局技术中心 基于实时雷达数据的高精度4d飞机航迹分析方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"基于计划到达时刻的四维航迹规划";刘杰 等;《航空计算技术》;20160731;第46卷(第4期);全文 *
"基于BADA及航空器意图的四维航迹预测";张军峰 等;《西南交通大学学报》;20140630;第49卷(第3期);前言、第1-3节 *
"支持连续下降运行的飞行管理技术研究";齐林 等;《2016第五届民用飞机航电系统国际论坛论文集》;20160421;第3节 *

Also Published As

Publication number Publication date
CN106648673A (zh) 2017-05-10

Similar Documents

Publication Publication Date Title
CN106648673B (zh) 一种连续下降运行程序的分析和设计方法
CN110930770B (zh) 一种基于管制意图和飞机性能模型的四维航迹预测方法
CN110046735B (zh) 基于飞行数据分析的飞机离场燃油消耗评估方法
CN101692315B (zh) 基于实时雷达数据的高精度4d飞机航迹分析方法
CN106529093A (zh) 一种针对大展弦比机翼的气动/结构/静气弹耦合优化方法
CN106768123A (zh) 一种无人直升机燃油预估方法
CN105303896A (zh) 航班预计到达时间精确预估方法
CN107016880B (zh) 一种基于模拟退火算法的连续下降进近油耗优化方法
JP2020077387A (ja) 垂直飛行経路の最適化
CN109738035A (zh) 基于ads-b航迹数据的航空器油耗计算方法
CN112669651B (zh) 一种基于飞行动态情报中eet值修正预测过点时间的方法
EP2899685A1 (en) Systems and methods for estimating net solar energy production for airborne photovoltaic systems
Porretta et al. Performance evaluation of a novel 4D trajectory prediction model for civil aircraft
Fernandes de Oliveira et al. Benefits of optimal flight planning on noise and emissions abatement at the Frankfurt airport
Jiang et al. Research on method of trajectory prediction in aircraft flight based on aircraft performance and historical track data
Bronsvoort et al. Enhanced descent wind forecast for aircraft
Schuster et al. High-performance trajectory prediction for civil aircraft
CN105759630A (zh) 基于模糊pid控制的飞机4d航迹仿真系统及仿真方法
CN107688689B (zh) 一种基于分层加权的飞行程序噪声评估方法
Ellerbroek et al. Fuel and emission benefits for continuous descent approaches at Schiphol
Pisani et al. City-pair trajectory optimization in the presence of winds using the GATAC framework
de Oliveira et al. On-board trajectory optimization of rnav departure and arrival procedures concerning emissions and population annoyance
Liu et al. Research on Trajectory Generation and Optimization in Continuous Descent Operations
Huang et al. Aircraft Departure Fuel Consumption Improvement Model Based On Flight Data
CN104881559B (zh) 微下冲气流工程化仿真建模方法

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