CN107085978A - 一种基于所需到达时间的管制辅助决策指令生成方法 - Google Patents
一种基于所需到达时间的管制辅助决策指令生成方法 Download PDFInfo
- Publication number
- CN107085978A CN107085978A CN201710474250.2A CN201710474250A CN107085978A CN 107085978 A CN107085978 A CN 107085978A CN 201710474250 A CN201710474250 A CN 201710474250A CN 107085978 A CN107085978 A CN 107085978A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- leg
- airborne vehicle
- sub
- 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
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0073—Surveillance aids
- G08G5/0091—Surveillance aids for monitoring atmospheric conditions
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Quality & Reliability (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Development Economics (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开一种基于所需到达时间的管制辅助决策指令生成方法,步骤是:构建四维航迹预测所需的航空器初始意图模型、外部环境模型与质点模型;对航空器进行航迹预测,得到初始的四维航迹以及航路点的预计到达时间;利用管制自动化***或管制辅助决策***中航班排序、冲突解脱的功能,获得航路点的所需到达时间;对航空器初始的四维航迹进行优化;在空中交通管制辅助***中显示航空器预计的水平轨迹与垂直剖面,并直接生成管制员应当向飞行员下达的管制指令。此种方法通过航迹预测与航迹优化快速生成管制指令,其规划的指令与航迹对管制员和飞行员而言操作简单,可控性强,并能使航空器按照所需到达时间通过航路点。
Description
技术领域
本发明属于民用航空器空中交通管理技术领域,涉及空中交通管制自动化与智能化、空中交通管制决策支持工具领域,特别涉及一种基于所需到达时间的管制辅助决策指令生成方法。
背景技术
随着空中交通流量的迅猛增长,空域拥挤状况愈发严重,因此而引起的航班延误率升高、冲突调配难度加大已成为目前管制工作中的一大难题。与此同时,管制的工作压力也与日俱增,以往凭借经验判断的管制方法已很难满足当前运行的需求,必须在管制设备和管制手段上寻求革新。对于空中交通管理部门来说,管制辅助决策***是一种提高效率、减轻负荷与保障安全的有效工具。
管制辅助决策***是以航迹预测为基础的:根据当前航空器状态、对飞行员和管制员意图的估计、结合气象信息与航空器性能,推算航空器在未来一段时间内的飞行状态;管制辅助决策***根据航迹推算的结果,使用计算机快速有效地对航空器进行排序、冲突探测,进而依据一定的规则进行调度和冲突解脱。鉴于此,辅助决策***的核心和关键是控制过点时间,即控制到达某一航路点的所需到达时间(Required Time of Arrival,RTA)。
以下文章或专利是关于基于所需到达时间的管制辅助决策***方面的:
The Boeing Company.Method and System of Controlling A Flight of anAircraft Subjected to a Required Time of Arrival Constraint:United States,US20160379500A1[P].2016-05-24;
SHIH-YIH YOUNG M I,KRISTEN M JEROME M I.Predictable And Required Timeof Arrival Compliant Optimized Profile Descents With Four Dimensional FlightManagement System And Related Method:United States,US9193442 B1[P].2015-11-24;
MACWILLIAMS P V,ZAKRZEWSKI E.Terminal Area Required Time of Arrival(RTA)Concept of operations and Automation Prototype,AIAA-2008-8930[R],2008;
SCHARL J,HARALDSDOTTIR A,KING J,et al.A Fast-Time required time ofarrival(RTA)model for analysis of 4D arrival management concepts,AIAA-2008-7027[R],2008.
然而值得注意的是,上述研究目前尚停留在运行概念设计与原型***验证阶段,离实际应用存在以下两个缺陷:一是管制辅助***仅仅是为管制员提供“应该达到的效果”,并没有具体的管制方法,这对管制员而言实施困难,仍然难以解决使航空器准时过点的问题;二是基于所需到达时间的航迹规划是由机载来完成的,相比于地面机载气象探测范围小、不能综合考虑空域内的航空器,同时该策略有待于地空数据通信技术取得突破性进展。
发明内容
本发明的目的,在于提供一种基于所需到达时间的管制辅助决策指令生成方法,其通过航迹预测与航迹优化快速生成管制指令,其规划的指令与航迹对管制员和飞行员而言操作简单,可控性强,并能使航空器按照所需到达时间通过航路点。
为了达成上述目的,本发明的解决方案是:
一种基于所需到达时间的管制辅助决策指令生成方法,包括如下步骤:
步骤1,构建四维航迹预测所需的航空器初始意图模型、外部环境模型与质点模型;
步骤2,按照步骤1中建立的意图模型、环境模型、质点模型,对航空器进行航迹预测,得到初始的四维航迹以及航路点的预计到达时间;
步骤3,利用管制自动化***或管制辅助决策***中航班排序、冲突解脱的功能,获得航路点的所需到达时间;
步骤4,由步骤3中的所需到达时间,对航空器初始的四维航迹进行优化;
步骤5,在空中交通管制辅助***中显示航空器预计的水平轨迹与垂直剖面,并直接生成管制员应当向飞行员下达的管制指令。
上述步骤1中,航空器的初始意图模型的构建方法是:根据飞行计划与航路点坐标,建立航空器的水平航迹;结合数据库编码表与管制移交协议,确定航空器所要经过的航路点的速度限制和高度限制;由此得到航空器的初始意图模型。
上述步骤1中,外部环境模型的构建方法是:
步骤1a,根据温度偏差和气压高度,确定大气温度T:
T=T0+ΔT+βT·Hp
其中,T0=288.15K,表示国际标准大气条件下在平均海平面处的温度;ΔT表示温度偏差;Hp表示气压高度;βT=-0.0065K/m,表示温度垂直递减率;
步骤1b,根据大气温度T,确定大气压力p:
其中,p0=101325Pa,表示国际标准大气条件下的空气压力;g0=9.80665m/s2,表示重力加速度;R=287.05287m2/(K·s2),表示空气常数;
步骤1c,根据温度T与压力p,确定大气密度ρ:
步骤1d,根据气象预报中的风向风速,结合大气温度、大气压力和大气密度,建立航空器运行的环境模型。
上述步骤1中,质点模型的构建方法是:
步骤10a,计算航空器发动机的推力,其最大起飞推力Thrmax climb如下式:
Thrmax climb=CTc,1·(1-h/CTc,2+CTc,3·h2)·(1-CTc,5·ΔT)
其中,h为大地高,CTc,1、CTc,2、CTc,3和CTc,5均为推力系数;
步骤10b,根据下式计算航空器阻力:
其中,VTAS为航空器真空速;CD为阻力系数,S为机翼参考面积;
步骤10c,航空器运动学方程如下:
其中,m表示航空器质量,d/dt表示时间微分;
该方程推导为:
其中,定义f{M}为能量分配系数;
步骤10d,由航空器推力、阻力及运动学方程建立航空器质点模型。
上述步骤3中,管制自动化***或管制辅助决策***在指定航空器所需到达时间时,需在时间窗内;若所需到达时间超出时间窗范围,则通过调整步骤1中的意图模型来进行调整。
上述步骤4的具体步骤如下:
步骤41,由步骤2生成的初始四维航迹,按照意图模型中的一个约束点到下一个约束点,拆分为多个航段;
步骤42,将每个航段再次拆分为两个子航段:完成减速下降或加速爬升过程的定义为子航段1,完成等速平飞到下一个约束点的定义为子航段2;
步骤43,在基于所需到达时间的航迹优化过程中,将子航段2再次拆分为两个等速平飞的子航段:位于子航段1之前等速平飞的子航段2a,位于子航段1之后等速平飞的子航段2b;
步骤44,设在步骤41中拆分出的一个航段需要调整的时间为Δt,建立以下方程:
t0=t3
Δt=t2-t6
t1-t0=t5-t4
其中,t0表示优化前子航段1开始时间,t1表示优化前子航段1结束时间或子航段2开始时间,t2表示优化前子航段2结束时间;t3表示优化后子航段2a开始时间,t4表示优化后子航段2a结束时间或子航段1开始时间,t5表示优化后子航段1结束时间或子航段2b开始时间,t6表示优化后子航段2b结束时间;
步骤45,不考虑转弯过程中的距离损耗,优化前后航空器在每个航段飞过的水平距离相等:
其中,f1表示优化前的速度-时间函数,f2表示优化后的速度-时间函数;
步骤46,对每个航段可调整时间Δt作出如下限制:
下降航段:
爬升航段:
其中,VTAS0表示子航段1开始点的真空速(True Airspeed,TAS),VTAS1表示子航段1结束点的真空速;
步骤47,比较优化后的航空器过点时间与所需到达时间的差,若在误差范围内则优化完成;若超出误差范围,则调整优化量并重新进行优化,直到航空器过点时间与所需到达时间的差在误差范围内。
采用上述方案后,本发明提出了一种新的航迹生成与优化方式,在地面的管制辅助***中完成基于所需到达时间的航迹规划,并直接向管制员提供应当下达的管制指令,从而辅助管制员决策,达到提高管制情景意识、降低管制工作负荷的目的。具体而言,本发明具有如下有益效果:
(1)本发明采用的航迹优化方式是轨迹的拆分与组合,因此在优化的过程中不会改变航空器原来的下降减速和加速爬升方式,航空器的性能不会受到影响;对于飞行员来说,优化后的航迹并不改变原来的操作方式,仅涉及操作的时机;而对于管制员来说,明确的管制指令点以及管制方式,可以很方便的下达管制指令,同时提高管制准确性,减轻负荷;
(2)本发明实践性强,支持实时的轨迹优化与修正,即使由于人为失误导致适当可控的偏差,本发明的航迹优化方法可以在后期进行修正,消除前序偏差;
(3)本发明可以使航空器在所需到达时间到达指定航路点;
(4)本发明向管制员提供直接的管制建议,有助于提高管制情景意识、降低管制工作负荷。
附图说明
图1是本发明的流程图;
图2是航空器初始意图模型生成示意图;
图3是环境模型生成示意图;
图4是质点模型生成示意图;
图5是航迹优化过程示意图;
图6是下降航段航迹优化中的速度-时间剖面示意图;
图7是管制辅助决策***中管制指令生成及显示示意图。
具体实施方式
以下将结合附图,对本发明的技术方案及有益效果进行详细说明。
如图1所示,本发明提供一种基于所需到达时间的管制辅助决策指令生成方法,包括以下步骤:
步骤11,根据飞行计划与航路点坐标,建立航空器的水平航迹,结合数据库编码表与管制移交协议,确定航空器所要经过的航路点的速度限制和高度限制,由此得到航空器的初始意图模型;根据气象观测数据和气象预报,建立航空器运行期间的环境模型,包括不同高度与不同位置的大气温度、大气压力、大气密度和风场等数据;根据航空器动力学模型与运动学方程,建立航空器质点运动方程,基于航空器基本性能数据库(BADA)中的性能参数,建立航空器的质点模型;
步骤12,结合步骤11中的意图模型、环境模型、质点模型,对航空器进行航迹预测,并获取航路点的预计到达时间(Estimated Time of Arrival,ETA);初始的四维航迹通过步长预测累加递推而成,在步长预测前先比较航空器的当前状态与意图模型中下一个约束点的速度/高度关系,应当采用的飞行方式如表1所示:
表1
表1中,VCAS1表示航空器当前点的校正空速(Calibrated Air Speed,CAS),VCAS2表示意图模型中下一个约束点的校正空速;h1表示航空器当前点的高度(Altitude,ALT),h2表示意图模型中下一个约束点的高度。
步骤13,利用管制自动化***或管制辅助决策***中航班排序、冲突解脱的功能,获得航路点的所需到达时间(Required Time of Arrival,RTA);在步骤12中对一条水平航路进行航迹预测时会得到航路点的预计到达时间,同时得到航路点的最早和最晚到达时间,最早和最晚到达时间的这个时间区间称为时间窗;管制自动化***或管制辅助决策***在指定航空器RTA时,应当在时间窗内;若由于冲突等原因RTA超出时间窗范围,则应通过改变水平路径或等待策略来进行调整,即调整步骤11中的意图模型;
步骤14,基于步骤13中的所需到达时间,按照本发明中基于所需到达时间的四维航迹优化方法,对初始四维航迹进行优化;
步骤15,在空中交通管制辅助***中显示航空器预计的水平轨迹与垂直剖面,并直接生成管制员应当向飞行员下达的管制指令。
图2为航空器初始四维航迹生成示意图,具体包括以下步骤:
步骤21,根据飞行计划与航路点坐标,建立航空器的水平航迹;
步骤22,结合数据库编码表与管制移交协议,确定航空器所要经过的航路点的速度限制和高度限制;
步骤23,联合步骤21、步骤22建立航空器的初始意图模型。
图3为环境模型生成示意图,具体包括以下步骤:
步骤31,根据温度偏差和气压高度,确定大气温度T:
T=T0+ΔT+βT·Hp
上式中,T0=288.15K,表示国际标准大气条件下在平均海平面处的温度;ΔT表示温度偏差;Hp表示气压高度;βT=-0.0065K/m,表示温度垂直递减率;
步骤32,根据大气温度T,确定大气压力p:
上式中,p0=101325Pa,表示国际标准大气条件下的空气压力;g0=9.80665m/s2,表示重力加速度;R=287.05287m2/(K·s2),表示空气常数;
步骤33,根据温度T与压力p,确定大气密度ρ:
步骤34,根据气象预报中的风向风速,结合大气温度、大气压力和大气密度,建立航空器运行的环境模型;
图4为质点模型生成示意图,具体步骤如下:
步骤41,计算航空器发动机的推力,其最大起飞推力Thrmax climb如下式:
Thrmax climb=CTc,1·(1-h/CTc,2+CTc,3·h2)·(1-CTc,5·ΔT)
上式中,h为大地高,CTc,1、CTc,2、CTc,3和CTc,5均为推力系数,参见航空器基础资料(BADA),且起飞/爬升/巡航/下降/进近/着陆的推力可视作最大爬升推力的函数,与所处的高度以及飞行阶段相关;
步骤42,计算航空器阻力:
上式中,VTAS为航空器真空速;CD为阻力系数,S为机翼参考面积,各系数参见航空器基础资料(BADA);且航空器起飞/爬升/巡航/进近/着陆形态下的阻力系数均为CD的相关函数;
步骤43,航空器运动学方程:
上式中,m表示航空器质量,d/dt表示时间微分;
进一步地,该方程可推导为:
上式中,定义f{M}为能量分配系数;f{M}在各阶段的取值及计算参见航空器基础资料(BADA);
步骤44,由航空器推力、阻力及运动学方程建立航空器质点模型;
图5为航迹优化过程示意图,具体步骤如下:
步骤51,按照意图模型中的一个约束点到下一个约束点,拆分为多个航段;在下降航段中,可能会采用减速下降、等速下降、减速平飞、等速平飞四种飞行方式;在爬升航段中,可能会采用加速爬升、等速爬升、加速平飞、等速平飞四种飞行方式;但在两种类型的航段中,等速平飞的飞行方式往往紧随另外三种飞行方式;因此可将每个航段再次拆分为两个子航段:完成减速下降过程的定义为子航段1,完成等速平飞到下一个约束点的定义为子航段2;
步骤52,在基于所需到达时间的航迹优化过程中,需要将子航段2再次拆分为两个等速平飞的子航段:位于子航段1之前等速平飞的子航段2a,位于子航段1之后等速平飞的子航段2b;
步骤53,将拆分后的三个子航段:子航段1,子航段2a,子航段2b按照如下规则进行重组:
(a)设在步骤51中拆分出的一个航段需要调整的时间为Δt,建立以下方程:
t0=t3
Δt=t2-t6
t1-t0=t5-t4
上式中,t0表示优化前子航段1开始时间,t1表示优化前子航段1结束时间或子航段2开始时间,t2表示优化前子航段2结束时间;t3表示优化后子航段2a开始时间,t4表示优化后子航段2a结束时间或子航段1开始时间,t5表示优化后子航段1结束时间或子航段2b开始时间,t6表示优化后子航段2b结束时间;
(b)不考虑转弯过程中的距离损耗,优化前后航空器在每个航段飞过的水平距离相等:
上式中,f1表示优化前的速度-时间函数,f2表示优化后的速度-时间函数;
(c)对每个航段可调整时间Δt作出如下限制:
下降航段:
爬升航段:
上式中,VTAS0表示子航段1开始点的真空速(True Airspeed,TAS),VTAS1表示子航段1结束点的真空速。
步骤54,由于优化时改变了航空器经过转弯点的速度,引起转弯半径的变化,进而导致整个水平距离的变化。为了消除由总距离变化产生的时间误差,比较优化后的航空器过点时间与RTA的差,若在误差范围内则优化完成;若超出误差范围,则调整优化量并重新进行优化,直到航空器过点时间与RTA的差在误差范围内。
图6为下降航段航迹优化中的速度-时间剖面示意图,通过对航段的拆分与组合,航空器从一个航路点到下一航路点的时间缩短了Δt,且不会改变初始四维航迹在航路点的高度/速度;航空器按照剖面飞行,当到达剖面中高度/速度变化点时,管制员向飞行员下达一次指令。
图7为管制辅助***中管制指令生成及显示示意图。图7中(a)部分为指令提示区,时间轴下方黑色竖线表示指令执行线;每架航空器对应各自的进度条,进度条会随着时间轴一起向左移动:进度条共分为两种状态,黑色的表示执行段(速度变化或者高度变化),白色的表示等速平飞段;右边“指令区”表示管制指令;当黑色的执行段接近指令执行线时,管制员即可向飞行员下达管制指令。图7中(b)部分为航迹优化后的航空器水平进场轨迹及高度、速度剖面示意图。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (6)
1.一种基于所需到达时间的管制辅助决策指令生成方法,其特征在于包括如下步骤:
步骤1,构建四维航迹预测所需的航空器初始意图模型、外部环境模型与质点模型;
步骤2,按照步骤1中建立的意图模型、环境模型、质点模型,对航空器进行航迹预测,得到初始的四维航迹以及航路点的预计到达时间;
步骤3,利用管制自动化***或管制辅助决策***中航班排序、冲突解脱的功能,获得航路点的所需到达时间;
步骤4,由步骤3中的所需到达时间,对航空器初始的四维航迹进行优化;
步骤5,在空中交通管制辅助***中显示航空器预计的水平轨迹与垂直剖面,并直接生成管制员应当向飞行员下达的管制指令。
2.如权利要求1所述的一种基于所需到达时间的管制辅助决策指令生成方法,其特征在于:所述步骤1中,航空器的初始意图模型的构建方法是:根据飞行计划与航路点坐标,建立航空器的水平航迹;结合数据库编码表与管制移交协议,确定航空器所要经过的航路点的速度限制和高度限制;由此得到航空器的初始意图模型。
3.如权利要求1所述的一种基于所需到达时间的管制辅助决策指令生成方法,其特征在于:所述步骤1中,外部环境模型的构建方法是:
步骤1a,根据温度偏差和气压高度,确定大气温度T:
T=T0+ΔT+βT·Hp
其中,T0=288.15K,表示国际标准大气条件下在平均海平面处的温度;ΔT表示温度偏差;Hp表示气压高度;βT=-0.0065K/m,表示温度垂直递减率;
步骤1b,根据大气温度T,确定大气压力p:
<mrow>
<mi>p</mi>
<mo>=</mo>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>T</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>T</mi>
</mrow>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>g</mi>
<mn>0</mn>
</msub>
<mrow>
<msub>
<mi>&beta;</mi>
<mi>T</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>R</mi>
</mrow>
</mfrac>
</mrow>
</msup>
</mrow>
其中,p0=101325Pa,表示国际标准大气条件下的空气压力;g0=9.80665m/s2,表示重力加速度;R=287.05287m2/(K·s2),表示空气常数;
步骤1c,根据温度T与压力p,确定大气密度ρ:
<mrow>
<mi>&rho;</mi>
<mo>=</mo>
<mfrac>
<mi>p</mi>
<mrow>
<mi>R</mi>
<mo>&CenterDot;</mo>
<mi>T</mi>
</mrow>
</mfrac>
</mrow>
步骤1d,根据气象预报中的风向风速,结合大气温度、大气压力和大气密度,建立航空器运行的环境模型。
4.如权利要求1所述的一种基于所需到达时间的管制辅助决策指令生成方法,其特征在于:所述步骤1中,质点模型的构建方法是:
步骤10a,计算航空器发动机的推力,其最大起飞推力Thrmax climb如下式:
Thrmax climb=CTc,1·(1-h/CTc,2+CTc,3·h2)·(1-CTc,5·ΔT)
其中,h为大地高,CTc,1、CTc,2、CTc,3和CTc,5均为推力系数;
步骤10b,根据下式计算航空器阻力:
<mrow>
<mi>D</mi>
<mo>=</mo>
<mn>0.5</mn>
<mo>&CenterDot;</mo>
<msub>
<mi>C</mi>
<mi>D</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>&rho;</mi>
<mo>&CenterDot;</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>T</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>&CenterDot;</mo>
<mi>S</mi>
</mrow>
其中,VTAS为航空器真空速;CD为阻力系数,S为机翼参考面积;
步骤10c,航空器运动学方程如下:
<mrow>
<mo>(</mo>
<mi>T</mi>
<mi>h</mi>
<mi>r</mi>
<mo>-</mo>
<mi>D</mi>
<mo>)</mo>
<mo>&CenterDot;</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>T</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>mg</mi>
<mn>0</mn>
</msub>
<mfrac>
<mrow>
<mi>d</mi>
<mi>h</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>mV</mi>
<mrow>
<mi>T</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msub>
<mi>dV</mi>
<mrow>
<mi>T</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
其中,m表示航空器质量,d/dt表示时间微分;
该方程推导为:
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>h</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mi>h</mi>
<mi>r</mi>
<mo>-</mo>
<mi>D</mi>
<mo>)</mo>
<mo>&CenterDot;</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>T</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>mg</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>V</mi>
<mrow>
<mi>T</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
<msub>
<mi>g</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>dV</mi>
<mrow>
<mi>T</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>h</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mi>h</mi>
<mi>r</mi>
<mo>-</mo>
<mi>D</mi>
<mo>)</mo>
<mo>&CenterDot;</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>T</mi>
<mi>A</mi>
<mi>S</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>mg</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mi>f</mi>
<mo>{</mo>
<mi>M</mi>
<mo>}</mo>
</mrow>
其中,定义f{M}为能量分配系数;
步骤10d,由航空器推力、阻力及运动学方程建立航空器质点模型。
5.如权利要求1所述的一种基于所需到达时间的管制辅助决策指令生成方法,其特征在于:所述步骤3中,管制自动化***或管制辅助决策***在指定航空器所需到达时间时,需在时间窗内;若所需到达时间超出时间窗范围,则通过调整步骤1中的意图模型来进行调整。
6.如权利要求1所述的一种基于所需到达时间的管制辅助决策指令生成方法,其特征在于:所述步骤4的具体步骤如下:
步骤41,由步骤2生成的初始四维航迹,按照意图模型中的一个约束点到下一个约束点,拆分为多个航段;
步骤42,将每个航段再次拆分为两个子航段:完成减速下降或加速爬升过程的定义为子航段1,完成等速平飞到下一个约束点的定义为子航段2;
步骤43,在基于所需到达时间的航迹优化过程中,将子航段2再次拆分为两个等速平飞的子航段:位于子航段1之前等速平飞的子航段2a,位于子航段1之后等速平飞的子航段2b;
步骤44,设在步骤41中拆分出的一个航段需要调整的时间为Δt,建立以下方程:
t0=t3
Δt=t2-t6
t1-t0=t5-t4
其中,t0表示优化前子航段1开始时间,t1表示优化前子航段1结束时间或子航段2开始时间,t2表示优化前子航段2结束时间;t3表示优化后子航段2a开始时间,t4表示优化后子航段2a结束时间或子航段1开始时间,t5表示优化后子航段1结束时间或子航段2b开始时间,t6表示优化后子航段2b结束时间;
步骤45,不考虑转弯过程中的距离损耗,优化前后航空器在每个航段飞过的水平距离相等:
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msubsup>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>4</mn>
</msub>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msubsup>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msubsup>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>3</mn>
</msub>
<msub>
<mi>t</mi>
<mn>6</mn>
</msub>
</msubsup>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
其中,f1表示优化前的速度-时间函数,f2表示优化后的速度-时间函数;
步骤46,对每个航段可调整时间Δt作出如下限制:
下降航段:
爬升航段:
其中,VTAS0表示子航段1开始点的真空速(True Airspeed,TAS),VTAS1表示子航段1结束点的真空速;
步骤47,比较优化后的航空器过点时间与所需到达时间的差,若在误差范围内则优化完成;若超出误差范围,则调整优化量并重新进行优化,直到航空器过点时间与所需到达时间的差在误差范围内。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710474250.2A CN107085978B (zh) | 2017-06-21 | 2017-06-21 | 一种基于所需到达时间的管制辅助决策指令生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710474250.2A CN107085978B (zh) | 2017-06-21 | 2017-06-21 | 一种基于所需到达时间的管制辅助决策指令生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107085978A true CN107085978A (zh) | 2017-08-22 |
CN107085978B CN107085978B (zh) | 2020-02-18 |
Family
ID=59606208
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710474250.2A Active CN107085978B (zh) | 2017-06-21 | 2017-06-21 | 一种基于所需到达时间的管制辅助决策指令生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107085978B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109191919A (zh) * | 2018-08-01 | 2019-01-11 | 清华大学 | 基于所需到达时间限制的四维轨迹速度调节方法及装置 |
CN109255981A (zh) * | 2018-08-31 | 2019-01-22 | 温州云航信息科技有限公司 | 基于关注焦点感知的增强型管制指令支持***及其方法 |
CN109491405A (zh) * | 2018-11-23 | 2019-03-19 | 中国航空工业集团公司沈阳飞机设计研究所 | 加油机和受油机的自主会和方法及装置 |
CN109615937A (zh) * | 2019-01-18 | 2019-04-12 | 南京航空航天大学 | 航段运行时间优化计算装置 |
CN109859531A (zh) * | 2018-12-04 | 2019-06-07 | 中国航空无线电电子研究所 | 一种面向飞行员不完整输入计算预测点预报风的方法 |
CN110059863A (zh) * | 2019-03-25 | 2019-07-26 | 南京航空航天大学 | 一种基于所需到达时间的航空器四维航迹优化方法 |
CN110111792A (zh) * | 2018-02-01 | 2019-08-09 | 雷明凯 | 空中交通管制自动指挥*** |
CN110276993A (zh) * | 2019-07-26 | 2019-09-24 | 郭静方 | 一种基于4d航迹的航空管制方法及*** |
CN110349445A (zh) * | 2018-04-05 | 2019-10-18 | 通用电气航空***有限责任公司 | 具有多个rta约束的有效飞行剖面 |
CN110853335A (zh) * | 2019-11-14 | 2020-02-28 | 东南大学 | 高速公路常发性瓶颈路段协作车队冲突避险自主决策方法 |
CN110889543A (zh) * | 2019-11-20 | 2020-03-17 | 西北工业大学 | 一种基于四维航迹运行的rta序列分配方法 |
CN111009154A (zh) * | 2019-10-18 | 2020-04-14 | 中国电子科技集团公司第二十八研究所 | 基于航班流预排序的终端区空域冲突探测与解脱算法 |
CN111344645A (zh) * | 2017-11-14 | 2020-06-26 | 湾流航空航天公司 | 轨迹建模中校准空速和真实空速之间的转换 |
CN112037583A (zh) * | 2020-09-02 | 2020-12-04 | 中国航空无线电电子研究所 | 一种航空器高度层变更程序优化引导方法 |
CN113269991A (zh) * | 2021-04-22 | 2021-08-17 | 南京莱斯信息技术股份有限公司 | 基于实时航迹与管制意图的空中交通中期冲突探测方法 |
CN115292557A (zh) * | 2022-07-29 | 2022-11-04 | 深圳微品致远信息科技有限公司 | 滑跑起飞的推算方法、装置、计算机设备和存储介质 |
CN116029570A (zh) * | 2023-03-24 | 2023-04-28 | 中国科学院心理研究所 | 一种塔台管制员辅助决策***的构建方法及装置 |
CN117978916A (zh) * | 2024-04-01 | 2024-05-03 | 中国民用航空飞行学院 | 一种管制员通话负荷预测方法及设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101692315A (zh) * | 2009-09-25 | 2010-04-07 | 民航总局空管局技术中心 | 基于实时雷达数据的高精度4d飞机航迹分析方法 |
US20160069688A1 (en) * | 2014-09-04 | 2016-03-10 | Honeywell International Inc. | System and method for managing speed constraints during required time of arrival operations |
CN106340208A (zh) * | 2015-01-07 | 2017-01-18 | 江苏理工学院 | 基于4d航迹运行的空中交通管制***的工作方法 |
CN106648673A (zh) * | 2016-12-28 | 2017-05-10 | 南京航空航天大学 | 一种连续下降运行程序的分析和设计方法 |
CN106846919A (zh) * | 2017-01-16 | 2017-06-13 | 南京航空航天大学 | 一种基于ads‑b信息更新的四维航迹动态预测方法 |
-
2017
- 2017-06-21 CN CN201710474250.2A patent/CN107085978B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101692315A (zh) * | 2009-09-25 | 2010-04-07 | 民航总局空管局技术中心 | 基于实时雷达数据的高精度4d飞机航迹分析方法 |
US20160069688A1 (en) * | 2014-09-04 | 2016-03-10 | Honeywell International Inc. | System and method for managing speed constraints during required time of arrival operations |
CN106340208A (zh) * | 2015-01-07 | 2017-01-18 | 江苏理工学院 | 基于4d航迹运行的空中交通管制***的工作方法 |
CN106648673A (zh) * | 2016-12-28 | 2017-05-10 | 南京航空航天大学 | 一种连续下降运行程序的分析和设计方法 |
CN106846919A (zh) * | 2017-01-16 | 2017-06-13 | 南京航空航天大学 | 一种基于ads‑b信息更新的四维航迹动态预测方法 |
Non-Patent Citations (5)
Title |
---|
OSTWALD PAUL: "Impacts of ATC related maneuvers on meeting a required time of arrival", 《DIGITAL AVIONICS SYSTEMS CONFERENCE》 * |
VILARDAGA S: "Operating cost sensitivity to required time of arrival commands to ensure separation in optimal aircraft 4D trajectories", 《TRANSPORTATION RESEARCH PART C》 * |
刘杰: "基于计划到达时刻的四维航迹规划", 《航空计算技术》 * |
武晓光: "进场航空器动态优化排序研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
马广辉: "终端空域进场航空器航迹规划研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111344645A (zh) * | 2017-11-14 | 2020-06-26 | 湾流航空航天公司 | 轨迹建模中校准空速和真实空速之间的转换 |
CN111344645B (zh) * | 2017-11-14 | 2023-12-12 | 湾流航空航天公司 | 轨迹建模中校准空速和真实空速之间的转换 |
CN110111792A (zh) * | 2018-02-01 | 2019-08-09 | 雷明凯 | 空中交通管制自动指挥*** |
CN110111792B (zh) * | 2018-02-01 | 2021-02-19 | 雷明凯 | 空中交通管制自动指挥*** |
CN110349445A (zh) * | 2018-04-05 | 2019-10-18 | 通用电气航空***有限责任公司 | 具有多个rta约束的有效飞行剖面 |
CN110349445B (zh) * | 2018-04-05 | 2022-04-29 | 通用电气航空***有限责任公司 | 具有多个rta约束的有效飞行剖面 |
CN109191919A (zh) * | 2018-08-01 | 2019-01-11 | 清华大学 | 基于所需到达时间限制的四维轨迹速度调节方法及装置 |
CN109255981B (zh) * | 2018-08-31 | 2022-09-13 | 温州云航信息科技有限公司 | 基于关注焦点感知的增强型管制指令支持***及其方法 |
CN109255981A (zh) * | 2018-08-31 | 2019-01-22 | 温州云航信息科技有限公司 | 基于关注焦点感知的增强型管制指令支持***及其方法 |
CN109491405A (zh) * | 2018-11-23 | 2019-03-19 | 中国航空工业集团公司沈阳飞机设计研究所 | 加油机和受油机的自主会和方法及装置 |
CN109859531A (zh) * | 2018-12-04 | 2019-06-07 | 中国航空无线电电子研究所 | 一种面向飞行员不完整输入计算预测点预报风的方法 |
CN109615937A (zh) * | 2019-01-18 | 2019-04-12 | 南京航空航天大学 | 航段运行时间优化计算装置 |
CN110059863A (zh) * | 2019-03-25 | 2019-07-26 | 南京航空航天大学 | 一种基于所需到达时间的航空器四维航迹优化方法 |
CN110059863B (zh) * | 2019-03-25 | 2023-03-24 | 南京航空航天大学 | 一种基于所需到达时间的航空器四维航迹优化方法 |
CN110276993A (zh) * | 2019-07-26 | 2019-09-24 | 郭静方 | 一种基于4d航迹的航空管制方法及*** |
CN111009154A (zh) * | 2019-10-18 | 2020-04-14 | 中国电子科技集团公司第二十八研究所 | 基于航班流预排序的终端区空域冲突探测与解脱算法 |
CN110853335A (zh) * | 2019-11-14 | 2020-02-28 | 东南大学 | 高速公路常发性瓶颈路段协作车队冲突避险自主决策方法 |
CN110889543B (zh) * | 2019-11-20 | 2023-06-09 | 西北工业大学 | 一种基于四维航迹运行的rta序列分配方法 |
CN110889543A (zh) * | 2019-11-20 | 2020-03-17 | 西北工业大学 | 一种基于四维航迹运行的rta序列分配方法 |
CN112037583B (zh) * | 2020-09-02 | 2022-08-09 | 中国航空无线电电子研究所 | 一种航空器高度层变更程序优化引导方法 |
CN112037583A (zh) * | 2020-09-02 | 2020-12-04 | 中国航空无线电电子研究所 | 一种航空器高度层变更程序优化引导方法 |
CN113269991A (zh) * | 2021-04-22 | 2021-08-17 | 南京莱斯信息技术股份有限公司 | 基于实时航迹与管制意图的空中交通中期冲突探测方法 |
CN115292557A (zh) * | 2022-07-29 | 2022-11-04 | 深圳微品致远信息科技有限公司 | 滑跑起飞的推算方法、装置、计算机设备和存储介质 |
CN115292557B (zh) * | 2022-07-29 | 2023-08-25 | 深圳微品致远信息科技有限公司 | 滑跑起飞的推算方法、装置、计算机设备和存储介质 |
CN116029570A (zh) * | 2023-03-24 | 2023-04-28 | 中国科学院心理研究所 | 一种塔台管制员辅助决策***的构建方法及装置 |
CN117978916A (zh) * | 2024-04-01 | 2024-05-03 | 中国民用航空飞行学院 | 一种管制员通话负荷预测方法及设备 |
CN117978916B (zh) * | 2024-04-01 | 2024-05-28 | 中国民用航空飞行学院 | 一种管制员通话负荷预测方法及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN107085978B (zh) | 2020-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107085978A (zh) | 一种基于所需到达时间的管制辅助决策指令生成方法 | |
JP5473672B2 (ja) | 連続降下アプローチによる航空機の予測可能性の最大化 | |
JP6033549B2 (ja) | 航空機の飛行経路を予測するためのデータの供給 | |
Patron et al. | Flight trajectories optimization under the influence of winds using genetic algorithms | |
AU2012350148A1 (en) | A flight prediction system | |
CN110059863B (zh) | 一种基于所需到达时间的航空器四维航迹优化方法 | |
EP3598261B1 (en) | Method and system for determining a descent profile | |
Erzberger et al. | A time-based concept for terminal-area traffic management | |
No et al. | Control and simulation of arbitrary flight trajectory-tracking | |
Ghazi et al. | Method to calculate Cessna Citation X aircraft climb and cruise trajectory using an aero-propulsive model | |
Murrieta-Mendoza et al. | Method to calculate aircraft VNAV trajectory cost using a performance database | |
Fernandes de Oliveira et al. | Benefits of optimal flight planning on noise and emissions abatement at the Frankfurt airport | |
Guan et al. | A study of 4D trajectory prediction based on machine deep learning | |
Dalmau et al. | Performance comparison between TEMO and a typical FMS in presence of CTA and wind uncertainties | |
Gu et al. | Real-time optimization of short-term flight profiles to control time of arrival | |
Dalmau et al. | Performance comparison of guidance strategies to accomplish RTAs during a CDO | |
CN105759630B (zh) | 基于模糊pid控制的飞机4d航迹仿真***及仿真方法 | |
Prats et al. | Evaluation of in-flight trajectory optimisation with time constraints in a moving base flight simulator | |
Chan et al. | Improving and validating ctas performance models | |
Gierszewski et al. | Toward Onboard Trajectory Optimization for Fuel-Saving Climb of Aircraft with Automatic Flight Control | |
Battipede et al. | Analysis of the impact of performance model accuracy on 4d trajectory optimization | |
Royer | Design of an automatic landing system for the meridian UAV using fuzzy Logic | |
Tobias et al. | Mixing four-dimensional equipped and unequipped aircraft in the terminal area | |
Jiang et al. | Study on 4D track prediction and optimization of aircraft climbing phase | |
van der Plas et al. | Trajectory optimisation of an aerobatic air race |
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 |