CN110059863A - 一种基于所需到达时间的航空器四维航迹优化方法 - Google Patents
一种基于所需到达时间的航空器四维航迹优化方法 Download PDFInfo
- Publication number
- CN110059863A CN110059863A CN201910226557.XA CN201910226557A CN110059863A CN 110059863 A CN110059863 A CN 110059863A CN 201910226557 A CN201910226557 A CN 201910226557A CN 110059863 A CN110059863 A CN 110059863A
- Authority
- CN
- China
- Prior art keywords
- aircraft
- air route
- optimized
- arrival time
- node
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000008569 process Effects 0.000 claims abstract description 9
- 230000006870 function Effects 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 3
- 230000007306 turnover Effects 0.000 claims description 3
- 239000000523 sample Substances 0.000 abstract description 2
- 230000007474 system interaction Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 8
- 238000010276 construction Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 239000000446 fuel Substances 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 238000005194 fractionation Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- 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"
-
- 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
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
-
- 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
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种基于所需到达时间的航空器四维航迹优化方法,属于空中航行器交通控制***的技术领域。该方法,创建链表结构的航路数据并设置灵活三边的范围限制,构建四维航迹计算的四个基本模型并依据航迹生成流程对航空器的到达时间、调度范围进行预测,基于该方法的***与进场管理***交互数据,利用进场管理***的冲突探测与调度排序功能得到航空器的所需到达时间,由所需到达时间对航空器的水平轨迹、垂直剖面分别进行优化,在进场管理***中显示航迹优化结果。本发明能够在合理范围内为航空器规划符合到达时间要求的四维航迹,为管制员提供决策辅助和解决方法。
Description
技术领域
本发明公开了一种基于所需到达时间的航空器四维航迹优化方法,涉及空中 交通管制自动化与智能化、空中交通管制决策支持工具,属于空中航行器交通控 制***的技术领域。
背景技术
四维航迹的管理体现在对航空器航迹的生成、预测和优化上。当前研究更多 集中在航迹的初始生成、预测方面以及如何提升它们的精度。准确的航迹初始生 成和预测能为管制员提供良好的态势感知,但对于诸如进场管理***提出的调度 方案却无法提供有效的解决方法。航迹优化即是一种带有目标的航迹生成,以航 迹初始生成和预测为基础,并将进场管理***的调度方案作为优化目标,在满足 运行规则的前提下为航空器制定四维航迹。
目前,有部分航迹优化研究是基于到达时间的,然而,其侧重点主要是机载 航迹,这样带来的问题其一是机载设备缺少全空域的气象、航班等信息,优化结 果顾此失彼;其二是机载优化使得空域管理者空中交通管制员丧失了信息与管制 的主动性。申请号为201710474250的发明申请公开了一种基于所需到达时间的 管制辅助决策指令生成方法,根据管理辅助决策***获得的航路点所需到达时间, 对初始的四维航迹进行拆分与组合确定飞行员的操作时机进而确定管制指令点 及管制方式,提高管制准确性,减轻负荷,然而,该申请从预添加航路中选择水 平航路优化得到的航迹单一,管制指令点及管制方式不灵活,存在调度能力弱的 缺陷。
本申请旨在提出一种对航空器的水平轨迹、垂直剖面分别进行优化的四维航 迹优化方法。
发明内容
本发明的发明目的是针对上述背景技术的不足,提供了一种基于所需到达时 间的航空器四维航迹优化方法,其以准确的四维航迹预测为基础,通过与管制自 动化***的信息交互,为进场航空器规划符合时间要求的四维航迹,解决了从预 添加航路中选择水平航路优化得到的航迹单一进而导致管制不灵活的技术问题。
本发明为实现上述发明目的采用如下技术方案:
一种基于所需到达时间的航空器四维航迹优化方法,包括如下步骤:
步骤1,分析机场终端区进场航路特征,创建链表结构的航路数据,设置灵 活三边的范围限制;
步骤2,构建计算四维航迹的四个基本模型:气象模型、质点模型、性能模 型、意图模型;
步骤3,依据航迹生成流程,运用步骤2中的四个基本模型对航空器的到达 时间、调度范围进行预测;
步骤4,基于本发明的***与进场管理***实现数据交互,利用进场管理系 统冲突探测与调度排序功能得到航空器的所需到达时间;
步骤5,根据步骤4中的所需到达时间对航空器的水平轨迹、垂直剖面分别 进行优化;
步骤6,在管制自动化***中显示航迹优化结果,为管制员提供决策辅助。
上述步骤1的具体实现方法如下:
步骤1a,分析机场终端区的标准进离场航路和历史雷达轨迹,规划实际运行 中航空器的可用航路;
步骤1b,将步骤1a中出现的所有航路点视作链表结点,每一结点代表一个 航路点,结点属性中记录航路点的代码、经纬度、高度/速度限制、特殊飞行限 制;
步骤1c,步骤1b中所提到的每个结点都拥有三个下一级指向结点,采用short 指向、star指向、long指向这三个属性标记每个结点的三个下一级指向结点,short 指向、star指向、long指向分别代表从当前结点到目的结点的最短航路、标准航 路和最长航路指向;
步骤1d,初始化航空器运行目的结点(如跑道结点)的short指向、star指 向、long指向均为NULL;
步骤1e,从步骤1b结点中选择两个恰当的结点,分别作为最短三边结点和 最长三边结点。
上述步骤2中,气象模型的构建方法是:
步骤21a,根据温度偏差和气压高度,确定大气温度T:
T=T0+ΔT+βT·Hp,
其中,T0=288.15K,表示国际标准大气条件下在平均海平面处的温度;ΔT 表示温度偏差;Hp表示气压高度;βT=-0.0065K/m,表示温度垂直递减率;
步骤21b,根据大气温度T,确定大气压力p:
其中,p0=101325Pa,表示国际标准大气条件下的空气压力;g=9.80665m/s2, 表示重力加速度;R=287.05287m2/(K·s2),表示空气常数;
步骤21c,根据温度T与压力p,确定大气密度ρ:
步骤21d,根据气象预报中的风向风速,结合大气温度、大气压力和大气密 度,建立航空器运行的气象模型。
上述步骤2中,质点模型的构建方法是:
步骤22a,建立航空器的位置换算模型:
其中,x、y表示航空器的水平位置,h表示航空器的垂直位置,VTAS表示航 空器真空速,γ为飞行路径角,ψ为航向角,W1和W2分别表示飞行时受到风的 正东向分量和正北向分量;
步骤22b,确定航空器真空速的变化方式:
其中,Thr、D分别表示航空器的推力与受到的阻力,m为航空器质量;
步骤22c,确定航空器航向角的变化方式:
其中,表示转弯坡度;
步骤22d,确定航空器质量变化方程:
其中,fFuel表示燃油流量。
上述步骤2中,性能模型的构建方法是:
步骤23a,计算航空器发动机的推力,其最大起飞推力Thrmax climb如下式:
Thrmax climb=CTc,1·(1-h/CTc,2+CTc,3·h2)·(1-CTc,5·ΔT),
其中,CTc,1、CTc,2、CTc,3和CTc,5均为推力系数,实际推力Thr则根据飞行阶 段由最大起飞推力修正得到;
步骤23b,根据下式计算航空器阻力D:
其中,CD为阻力系数,S为机翼参考面积。
上述步骤2中,意图模型的构建方法是:
步骤24a,结合航空器经纬度与终端区航路设置,判定航空器当前飞行航路 以及在航路中的位置;
步骤24b,根据航空器的当前位置与航向,捕捉航空器即将飞行的下一个航 路点;
步骤24c,依据步骤1中建立的链表航路,选择结点组成最短航路、最长航 路、标准进场航路,初步建立水平意图模型;
步骤24d,由飞行计划与管制移交协议,确定步骤24c中航路的高度/速度约 束,初步建立意图模型。
上述步骤3的具体实现方法如下:
步骤3a,采用步骤24c中的标准进场航路,使用步骤2中的航迹模型,计 算当前航空器的标准进场时间;
步骤3b,采用步骤24c中的最短进场航路,使用步骤2中的航迹模型,计 算当前航空器的最小进场时间;同样的,采用24c中的最长进场航路,计算得到 最大进场时间;
步骤3c,由最小与最大进场时间确定航空器着陆时间窗,即当前航空器的 时间调度范围。
上述步骤4的具体实现方法是:管制自动化***从航迹优化模块获取标准进 场时间及着陆时间窗,并利用其冲突探测与调度排序功能,指定航空器的所需到 达时间,该时间不小于着陆时间窗的最小值。
上述步骤5通过先水平、后垂直的分步规划方法实现整体的航迹优化,水平 航迹优化的实现方法如下:
步骤51a,若航空器未进入三边程序,选择步骤24c中的标准进场航路作为 初始待优化航路;若航空器已进入三边程序,选择标准三边长度为初始待优化航 路,设航空器仍需飞行的三边长度为d;
步骤51b,根据待优化航路的高度/速度限制范围计算当前待优化航路的着陆 时间窗,设计算结果为[tsmall,tlarge];
步骤51c,若tsmall≤RTA≤tlarge,选择当前航路为进场航路,水平航迹优化结 束;
步骤51d,若RTA<tsmall:(1)航空器未进入三边程序:将待优化航路中某 一结点的star指向结点替换为该结点的short指向结点,组成新的待优化航路, 执行步骤51b,逐步替换结点直至水平航迹优化结束;(2)航空器已进入三边 程序:设定航空器的三边飞行长度为d-Δd,其中,Δd为调整步长,逐步减步 长直至水平优化结束;
步骤51e,若RTA>tlarge:(1)航空器未进入三边程序:将待优化航路中某 一结点的star指向结点替换为该结点的long指向结点,组成新的待优化航路, 执行步骤51b,逐步替换结点直至水平航迹优化结束,如果所有结点的star指向 结点已替换为long指向结点后仍不满足步骤51c,则执行步骤51f;(2)航空器 已进入三边程序:设定航空器的三边飞行长度为d+Δd,逐步加步长直至水平优 化结束;
步骤51f,在步骤51a中加入等待程序,执行步骤51b;逐步增加等待次数 直至水平航迹优化结束。
上述步骤5中,垂直剖面的优化方法如下:
步骤52a,根据步骤51c优化得到的水平航迹生成初始四维航迹,按照意图 模型中的一个约束点到下一个约束点,将待优化航路拆分为多个航段;
步骤52b,将每个航段再次拆分为两个子航段:完成减速下降过程的子航段 1,及,完成等速平飞到下一个约束点的子航段2;
步骤52c,在基于所需到达时间的航迹优化过程中,将子航段2再次拆分为 两个等速平飞的子航段:位于子航段1之前等速平飞的子航段2a,及,位于子 航段1之后等速平飞的子航段2b;
步骤52d,设在步骤52a中拆分出的一个航段需要调整的时间为Δ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结束时间;
步骤52e,不考虑转弯过程中的距离损耗,优化前后航空器在每个航段飞过 的水平距离相等:
其中,f1表示优化前的速度-时间函数,f2表示优化后的速度-时间函数;
步骤52f,对每个航段可调整时间Δt作出如下限制:
其中,VTAS0表示子航段1开始点的真空速,VTAS1表示子航段1结束点的真空 速;
步骤52g,比较优化后的航空器过点时间与所需到达时间的差,若在误差范 围内则优化完成;若超出误差范围,则调整优化量并重新进行优化,直到航空器 过点时间与所需到达时间的差在误差范围内。
上述步骤6,点击管制自动化***中航空器标牌,在主界面显示水平航迹, 在辅助界面显示垂直航迹。
本发明采用上述技术方案,具有以下有益效果:
(1)本发明采用链表模式存储航路基础数据,每个链表数据记录了表征当 前航路点到目的航路点构成标准进场航路、最短进场航路、最长进场航路的指向 信息,以航空器所需到达时间位于待优化水平航路着陆时间窗内为目标,利用链 表数据灵活规划航路,再对优化水平轨迹后的四维航迹进行拆分和重组,在不改 变原有的航路限制与飞行操作方式的前提下调节航空器下降与减速时机,进而实 现过航点时间的控制。
(2)本发明使用覆盖面广,从管制自动化***捕捉到进场航空器至脱离进 场三边程序,本发明都可以在合理范围内进行有效航迹优化,水平航迹优化拥有 直飞、迂回、等待、长短三边等多种调节方法,时间调整范围大。
(3)本发明为管制员提供决策的依据,优化结果亦可作为管制员与飞行员 实现调度的方法。
附图说明
图1是本发明优化航空器四维航迹的流程图。
图2是链表化航路及灵活设置三边的示意图。
图3是气象模型生成的示意图。
图4(a)是质点模型中航空器俯仰受力的分析图,图4(b)是质点模型中 航空器偏航受力的分析图,图4(c)是质点模型中航空器横滚受力的分析图。
图5是性能模型生成的示意图。
图6(a)是意图模型中航空器捕捉非star航路点的示意图,图6(b)是意 图模型中航空器捕捉star航路点的示意图,图6(c)是意图模型中航空器捕捉 标准等待程序的示意图,图6(d)是意图模型中航空器捕捉盘旋一周等待的示 意图。
图7是水平航迹优化过程的示意图。
图8是垂直剖面优化过程的示意图。
图9是垂直航迹优化中速度-时间的剖面示意图。
具体实施方式
下面结合附图对发明的技术方案进行详细说明。
本发明提供的一种基于所需到达时间的航空器四维航迹优化方法如图1所示, 包括以下6个步骤:
步骤11,链表化航路点数据:分析机场终端区进场航路特征,创建链表结 构的航路数据,设置灵活三边的范围限制;
步骤12,构建计算四维航迹的四个基本模型:气象模型、质点模型、性能 模型、意图模型;
步骤13,运用步骤12中的四个基本模型对航空器的到达时间、调度范围进 行预测;
步骤14,基于本发明的***与进场管理***实现数据交互,利用进场管理 ***冲突探测与调度排序功能得到航空器的所需到达时间;
步骤15,根据所需到达时间对航空器的水平轨迹、垂直剖面分别进行优化;
步骤16,在管制自动化***中显示航迹优化结果,为管制员提供决策辅助。
如图2所示,链表化航路点及灵活设置三边的具体实现方法如下:
步骤21,设机场终端区的一条标准进场航路为:W1-W3-W4-W6-W8-RWY, 根据机场航图及历史雷达轨迹,分别规划出长航路:W1-W2-W3,短航路:W1-W4, 灵活三边区域:W5-W7;
步骤22,对结点W1而言,其short指向是结点W4,star指向是结点W3,long 指向是结点W2;而结点W3的short指向、star指向、long指向均为结点W4。 W5处是最短三边位置,W7处是最长三边位置;
步骤23,根据航路限制、管制规则及航空器性能,分机型设定结点 Wi(i=1,2,...,8)的高度限制、速度限制、特殊飞行限制,补充结点Wi(i=1,2,...,8)的 航路点代码、经纬度等基础属性;
步骤24,设定终结点RWY的short指向、star指向、long指向均为NULL。
如图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(a)、图4(b)、图4(c)所示,质 点模型的构建方法是:
步骤41,mg,L,Thr,D分别表示航空器飞行时受到的重力、升力、推力 和阻力,VTAS表示航空器的真空速,α,β表示飞行时的攻角与侧滑角,γ为飞行 路径角,Ψ为航向角,表示转弯坡度;在空中交通管理领域,一般认为攻角和 侧滑角为零,以推力、转弯坡度、航径角作为控制变量,以航空器的位置、真空 速、航向和质量作为状态变量;
步骤42,构建质点模型:
式中,x,y表示航空器的水平位置,h表示航空器的垂直位置;W1和W2分 别表示飞行时受到风的正东向分量和正北向分量;fFuel为燃油流量。
如图5所示,性能模型的构建方法如下:
步骤51,计算航空器发动机的推力,其最大起飞推力Thrmax climb如下式:
Thrmax climb=CTc,1·(1-h/CTc,2+CTc,3·h2)·(1-CTc,5·ΔT),
其中,CTc,1、CTc,2、CTc,3和CTc,5均为推力系数,实际推力Thr则根据飞行阶 段最大起飞推力修正得到;
步骤52,根据下式计算航空器阻力D:
其中,CD为阻力系数,S为机翼参考面积。
航空器捕捉非star航路点、star航路点、标准等待程序、盘旋一周等待的意 图如图6(a)、图6(b)、图6(c)、图6(d)所示,意图模型的构建方法如 下:
步骤61,结合航空器经纬度与终端区的航路设置,判定航空器当前飞行航 路以及在航路中的位置;
步骤62,根据航空器的当前位置与航向捕捉航空器的飞行意图,如图6(a)、 图6(b)所示,两架航空器在同一位置,但由于航向不同,故而捕捉到不同的 将飞航路点;如图6(c)、图6(d)所示,根据航向及轨迹判断两架航空器均 在执行转弯程序,图6(c)中的转弯开始点位于标准等待程序起始点,故而判 断图6(c)中的意图是执行标准等待程序,而图6(d)中的意图是盘旋一周等 待;
步骤63,按照链表化数据结构,将后续航路点信息添加至链表结点,初步 建立水平意图;
步骤64,依据飞行计划与管制移交协议确定步骤63中航路的高度和速度约 束,初步建立意图模型。
水平航迹优化过程如图7所示,具体实现方法如下:
步骤71,根据待优化航路的高度/速度约束范围,计算当前待优化航路的着 陆时间窗,设计算结果为[MinTime,MaxTime];
步骤72,若MinTime≤RTA≤MaxTime,选择当前航路为进场航路,水平航迹 优化结束;
步骤73,若RTA<MinTime:(1)航空器未进入三边程序:将待优化航路中 某一结点的star指向结点替换为该结点的short指向结点,更新待优化航路,执 行步骤72,逐步替换结点直至水平航迹优化结束;(2)航空器已进入三边程序: 设定航空器的三边飞行长度为d-Δd,其中,Δd为调整步长,逐步减步长直至 水平优化结束;
步骤74,若RTA>MaxTime:(1)航空器未进入三边程序:将待优化航路中 的某一结点的star指向结点替换为该结点的long指向结点,更新待优化航路, 执行步骤72,逐步替换结点直至水平航迹优化结束,如果所有结点的star指向 结点已替换为long指向结点后仍不满足步骤72,则执行步骤75;(2)航空器 已进入三边程序:设定航空器的三边飞行长度为d+Δd,逐步加步长直至水平优 化结束;
步骤75,在步骤71中加入等待程序,执行步骤72;逐步增加等待次数直至 水平航迹优化结束。
垂直航迹优化过程如图8所示,具体实现方法如下:
步骤81,按照意图模型中的一个约束点(添加了高度/速度约束的链表化航 路点)到下一个约束点,将待优化航路拆分为多个航段,下降航段中可能会采用 减速下降、等速下降、减速平飞、等速平飞四种飞行方式,爬升航段中可能会采 用加速爬升、等速爬升、加速平飞、等速平飞四种飞行方式,然而,两类航段中 等速平飞的飞行方式往往紧随另外三种飞行方式,因此,可将每个航段再次拆分 为两个子航段:完成减速下降过程的子航段1以及完成等速平飞到下一个约束点 的子航段2;
步骤82,在基于所需到达时间的航迹优化过程中,需要将子航段2再次拆 分为两个等速平飞的子航段:位于子航段1之前等速平飞的子航段2a,位于子 航段1之后等速平飞的子航段2b;
步骤83,将拆分后的三个子航段:子航段1,子航段2a,子航段2b按照如 下规则进行重组:
(a)设在步骤81中拆分出的一个航段需要调整的时间为Δt,建立以下方程:
上式中,t0表示优化前子航段1开始时间,t1表示优化前子航段1结束时间 或子航段2开始时间,t2表示优化前子航段2结束时间;t3表示优化后子航段2a 开始时间,t4表示优化后子航段2a结束时间或子航段1开始时间,t5表示优化 后子航段1结束时间或子航段2b开始时间,t6表示优化后子航段2b结束时间,
(b)不考虑转弯过程中的距离损耗,优化前后航空器在每个航段飞过的水 平距离相等:
上式中,f1表示优化前的速度-时间函数,f2表示优化后的速度-时间函数,
(c)对每个航段可调整时间Δt作出如下限制:
上式中,VTAS0表示子航段1开始点的真空速,VTAS1表示子航段1结束点的真 空速;
步骤84,由于优化时改变了航空器经过转弯点的速度,引起转弯半径的变 化,进而导致整个水平距离的变化。为了消除由总距离变化产生的时间误差,比 较优化后的航空器过点时间与RTA的差,若在误差范围内则优化完成;若超出 误差范围,则调整优化量并重新进行优化,直到航空器过点时间与RTA的差在 误差范围内。
图9是垂直航迹优化中的速度-时间剖面示意图,通过对航段的拆分与组合, 航空器从一个航路点到下一航路点的时间缩短了Δt,且不会改变初始四维航迹 在航路点的高度/速度。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围, 凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本 发明保护范围之内。
Claims (7)
1.一种基于所需到达时间的航空器四维航迹优化方法,其特征在于,
创建链表结构的航路数据,链表的每个结点记录一个航路点信息,具体包含:一个航路点的代码及经纬度信息,一个航路点到目的航路点最短航路的指向信息、一个航路点到目的航路点标准航路的指向信息、一个航路点到目的航路点最长航路的指向信息;
根据机场空域情况择标准三边结点、最短三边结点、最长三边结点;
根据计算四维航迹的基本模型预测航空器的标准进场时间及着陆时间窗;
根据航空器的标准进场时间及着陆时间窗指定航空器的所需到达时间;
根据航空器的所需到达时间并基于链表结构的航路数据优化航空器的水平轨迹;
根据航空器的所需到达时间和水平轨迹优化结果优化垂直剖面。
2.根据权利要求1所述一种基于所需到达时间的航空器四维航迹优化方法,其特征在于,根据计算四维航迹的基本模型预测航空器的标准进场时间及着陆时间窗的具体方法为:由计算四维航迹的基本模型确定航空器按照标准进场航路,最短进场航路、最长进场航路飞行时的进场时间,分别记为标准进场时间、最短进场时间、最长进场时间,最短进场时间和最长进场时间确定的时间区间为着陆时间窗。
3.根据权利要求1所述一种基于所需到达时间的航空器四维航迹优化方法,其特征在于,根据航空器的标准进场时间及着陆时间窗指定航空器的所需到达时间的具体方法为:利用进场管理***冲突探测与调度排序的功能指定不小于着陆时间窗最小值的所需到达时间。
4.根据权利要求1所述一种基于所需到达时间的航空器四维航迹优化方法,其特征在于,根据航空器的所需到达时间以及基于链表结构的航路数据优化航空器的水平轨迹的具体方法为:
确定待优化航路:航空器未进入三边程序时,以标准进场航路为待优化航路,反之,以标准三边航路为待优化航路;
结合意图模型,根据航迹预测确定待优化航路的着陆时间窗;
根据航空器的所需到达时间与待优化航路着陆时间窗的关系优化水平航路:
航空器的所需到达时间位于待优化航路着陆时间窗内时,以当前待优化航路为进场航路,结束水平轨迹优化,
航空器的所需到达时间小于待优化航路着陆时间窗的最小值时,对于以标准进场航路为待优化航路的情形,将待优化航路中某一结点的标准航路指向结点替换为该结点的最短航路指向结点后重新确定待优化航路的着陆时间窗,逐步替换待优化航路中的各结点直至航空器的所需到达时间位于待优化航路着陆时间窗内,对于以标准三边长度为待优化航路的情形,逐步减少三边飞行的长度直至航空器的所需到达时间位于待优化航路着陆时间窗内;
航空器的所需到达时间大于待优化航路着陆时间窗的最大值时,对于以标准进场航路为待优化航路的情形,将待优化航路中某一结点的标准航路指向结点替换为该结点的最长航路指向结点后重新确定待优化航路的着陆时间窗,逐步替换待优化航路中的各结点直至航空器的所需到达时间位于待优化航路着陆时间窗内,若在待优化航路所有结点的标准航路指向结点替换为最长航路指向结点后航空器的所需到达时间仍不在待优化航路着陆时间窗内时,增加等待时长重新确定待优化航路,对于以标准三边长度为待优化航路的情形,逐步增加三边飞行的长度直至航空器的所需到达时间位于待优化航路着陆时间窗内。
5.根据权利要求1所述一种基于所需到达时间的航空器四维航迹优化方法,其特征在于,所述意图模型的建立方法为:捕捉航空器的将要飞行结点,添加后续航路结点以完善链表结构的航路,依据飞行计划与管制移交协议确定链表结构航路的高度约束和速度约束。
6.根据权利要求1所述一种基于所需到达时间的航空器四维航迹优化方法,其特征在于,根据航空器的所需到达时间和水平轨迹优化结果优化垂直剖面的具体方法为:根据优化后的水平轨迹更新意图模型,将意图模型中相邻两航路点之间的航路拆分为完成减速下降过程的第一子航段以及完成等速平飞到下一航路点的第二子航段,将第二子航段拆分为位于第一子航段之前的等速平飞的子航段和位于第一子航段之后等速平飞的子航段,调节两个等速平飞子航段的长度,使优化后的过航路点时间与所需到达时间的差值满足阈值要求。
7.一种计算机设备,包括:存储器、处理器及存储在存储器上并在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1至6中任意一项所述基于所需到达时间的航空器四维航迹优化方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910226557.XA CN110059863B (zh) | 2019-03-25 | 2019-03-25 | 一种基于所需到达时间的航空器四维航迹优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910226557.XA CN110059863B (zh) | 2019-03-25 | 2019-03-25 | 一种基于所需到达时间的航空器四维航迹优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110059863A true CN110059863A (zh) | 2019-07-26 |
CN110059863B CN110059863B (zh) | 2023-03-24 |
Family
ID=67316308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910226557.XA Active CN110059863B (zh) | 2019-03-25 | 2019-03-25 | 一种基于所需到达时间的航空器四维航迹优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110059863B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111488848A (zh) * | 2020-04-17 | 2020-08-04 | 中国民航大学 | 一种对轨迹中跑马场形盘旋的特征识别方法 |
CN111898906A (zh) * | 2020-07-29 | 2020-11-06 | 成都民航空管科技发展有限公司 | 基于aman***的计算航班长三边过点的方法及*** |
CN112037583A (zh) * | 2020-09-02 | 2020-12-04 | 中国航空无线电电子研究所 | 一种航空器高度层变更程序优化引导方法 |
CN112215399A (zh) * | 2020-09-11 | 2021-01-12 | 中国航空工业集团公司成都飞机设计研究所 | 一种机载动态电子围栏的高确定数据结构设计方法 |
CN112817330A (zh) * | 2021-01-05 | 2021-05-18 | 北京联合大学 | 一种多无人机四维航迹协同规划方法及*** |
CN113269991A (zh) * | 2021-04-22 | 2021-08-17 | 南京莱斯信息技术股份有限公司 | 基于实时航迹与管制意图的空中交通中期冲突探测方法 |
CN115311905A (zh) * | 2022-10-12 | 2022-11-08 | 珠海翔翼航空技术有限公司 | 基于姿态和速度预测的多跑道飞行进场导航方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070100538A1 (en) * | 2005-10-31 | 2007-05-03 | Honeywell International Inc. | System and method for performing 4-dimensional navigation |
CN102509475A (zh) * | 2011-10-26 | 2012-06-20 | 南京航空航天大学 | 基于4d航迹运行的空中交通管制***及方法 |
CN106846919A (zh) * | 2017-01-16 | 2017-06-13 | 南京航空航天大学 | 一种基于ads‑b信息更新的四维航迹动态预测方法 |
CN107085978A (zh) * | 2017-06-21 | 2017-08-22 | 南京航空航天大学 | 一种基于所需到达时间的管制辅助决策指令生成方法 |
CN109147398A (zh) * | 2018-09-13 | 2019-01-04 | 四川九洲空管科技有限责任公司 | 一种机载防撞***目标监视跟踪性能优化方法 |
-
2019
- 2019-03-25 CN CN201910226557.XA patent/CN110059863B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070100538A1 (en) * | 2005-10-31 | 2007-05-03 | Honeywell International Inc. | System and method for performing 4-dimensional navigation |
CN102509475A (zh) * | 2011-10-26 | 2012-06-20 | 南京航空航天大学 | 基于4d航迹运行的空中交通管制***及方法 |
CN106846919A (zh) * | 2017-01-16 | 2017-06-13 | 南京航空航天大学 | 一种基于ads‑b信息更新的四维航迹动态预测方法 |
CN107085978A (zh) * | 2017-06-21 | 2017-08-22 | 南京航空航天大学 | 一种基于所需到达时间的管制辅助决策指令生成方法 |
CN109147398A (zh) * | 2018-09-13 | 2019-01-04 | 四川九洲空管科技有限责任公司 | 一种机载防撞***目标监视跟踪性能优化方法 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111488848A (zh) * | 2020-04-17 | 2020-08-04 | 中国民航大学 | 一种对轨迹中跑马场形盘旋的特征识别方法 |
CN111898906A (zh) * | 2020-07-29 | 2020-11-06 | 成都民航空管科技发展有限公司 | 基于aman***的计算航班长三边过点的方法及*** |
CN111898906B (zh) * | 2020-07-29 | 2022-06-17 | 成都民航空管科技发展有限公司 | 基于aman***的计算航班长三边过点的方法及*** |
CN112037583A (zh) * | 2020-09-02 | 2020-12-04 | 中国航空无线电电子研究所 | 一种航空器高度层变更程序优化引导方法 |
CN112037583B (zh) * | 2020-09-02 | 2022-08-09 | 中国航空无线电电子研究所 | 一种航空器高度层变更程序优化引导方法 |
CN112215399A (zh) * | 2020-09-11 | 2021-01-12 | 中国航空工业集团公司成都飞机设计研究所 | 一种机载动态电子围栏的高确定数据结构设计方法 |
CN112817330A (zh) * | 2021-01-05 | 2021-05-18 | 北京联合大学 | 一种多无人机四维航迹协同规划方法及*** |
CN112817330B (zh) * | 2021-01-05 | 2023-08-08 | 北京联合大学 | 一种多无人机四维航迹协同规划方法及*** |
CN113269991A (zh) * | 2021-04-22 | 2021-08-17 | 南京莱斯信息技术股份有限公司 | 基于实时航迹与管制意图的空中交通中期冲突探测方法 |
CN115311905A (zh) * | 2022-10-12 | 2022-11-08 | 珠海翔翼航空技术有限公司 | 基于姿态和速度预测的多跑道飞行进场导航方法及*** |
CN115311905B (zh) * | 2022-10-12 | 2022-12-30 | 珠海翔翼航空技术有限公司 | 基于姿态和速度预测的多跑道飞行进场导航方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN110059863B (zh) | 2023-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110059863A (zh) | 一种基于所需到达时间的航空器四维航迹优化方法 | |
EP3598262B1 (en) | Method and system for determining a climb profile | |
CN110930770B (zh) | 一种基于管制意图和飞机性能模型的四维航迹预测方法 | |
EP2631890A2 (en) | Methods for in-flight adjusting of a flight plan | |
KR102242237B1 (ko) | 비행 관리 시스템의 모의 장치 | |
JP6018433B2 (ja) | 航空機軌道に沿った気象データ選択 | |
CN106249592B (zh) | 用于航空器的优化下降和进场剖面的自动确定方法和*** | |
US11094206B2 (en) | Vertical flightpath optimization | |
BR102013004148A2 (pt) | Método para voo de uma aeronave | |
CN112927562B (zh) | 一种基于指定高度限制的航空器飞行高度剖面计算方法 | |
CN106768123A (zh) | 一种无人直升机燃油预估方法 | |
CN110363333A (zh) | 一种基于渐进梯度回归树的天气影响下空中通行能力的预测方法 | |
CN110727286B (zh) | 用于确定下降剖面的方法和*** | |
CN103578299B (zh) | 一种模拟航空器飞行过程的方法 | |
Girardet et al. | Generating optimal aircraft trajectories with respect to weather conditions | |
CN109919391A (zh) | 一种太阳能飞机自主寻优航迹规划方法和*** | |
Robinso et al. | A concurrent sequencing and deconfliction algorithm for terminal area air traffic control | |
CN110986954A (zh) | 基于灰狼优化算法的军用运输机航路规划方法 | |
Ramée et al. | Aircraft flight plan optimization with dynamic weather and airspace constraints | |
CN107610531A (zh) | 基于4d航迹和航路飞行经验数据的航班过点时刻预估方法 | |
CN111445063B (zh) | 基于飞行线路的起降点选取方法及装置 | |
Kiss-Toth et al. | A dynamic programming approach for 4D flight route optimization | |
Xue et al. | Optimal assignment of time of departure under severe weather | |
Drupka et al. | Automated flight planning method to facilitate the route planning process in predicted conditions | |
Matsuno et al. | Robust optimal guidance algorithm for required time of arrival operations using probabilistic weather forecasts |
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 |