CN116991171A - 基于能量的无动力滑降无人机着陆下滑道动态生成方法 - Google Patents

基于能量的无动力滑降无人机着陆下滑道动态生成方法 Download PDF

Info

Publication number
CN116991171A
CN116991171A CN202310850244.8A CN202310850244A CN116991171A CN 116991171 A CN116991171 A CN 116991171A CN 202310850244 A CN202310850244 A CN 202310850244A CN 116991171 A CN116991171 A CN 116991171A
Authority
CN
China
Prior art keywords
energy
aircraft
landing
downslide
section
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
CN202310850244.8A
Other languages
English (en)
Other versions
CN116991171B (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.)
Institute of Mechanics of CAS
Original Assignee
Institute of Mechanics of CAS
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 Institute of Mechanics of CAS filed Critical Institute of Mechanics of CAS
Priority to CN202310850244.8A priority Critical patent/CN116991171B/zh
Publication of CN116991171A publication Critical patent/CN116991171A/zh
Application granted granted Critical
Publication of CN116991171B publication Critical patent/CN116991171B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明提供了一种基于能量的无动力滑降无人机着陆下滑道动态生成方法，是从飞机飞至跑道后方对准跑道准备着陆开始，将飞机下滑道分为能量管理段、陡下滑段、圆弧过渡段和浅下滑段，着陆开始后，在能量管理段开始实时计算飞机当前能量，根据能量剖面计算浅下滑段起始高度，然后根据剩余航程分别计算飞机的浅下滑段、圆弧过渡段和陡下滑段对应航程与高度的关系，生成整个下滑道高度与航程曲线。本发明动态生成的下滑道可有效解决无人机降落初始状态不同导致的降落状态的不确定性，可在无减速板情况下仅通过动态生成的下滑道使飞机在降落时达到理想的速度与下沉率，对高度跟踪效果好，着陆精度高。

Description

基于能量的无动力滑降无人机着陆下滑道动态生成方法
技术领域
本发明涉及无人机技术领域，具体涉及一种基于能量的无动力滑降无人机着陆下滑道动态生成方法。
背景技术
固定翼无人机自动着陆是整个飞行的最后阶段，接地时飞机需处于理想的飞行速度及下沉率，否则易对飞机造成损坏。
固定翼无人机飞行前将下滑道装订至飞机，降落阶段严格按装订的下滑道飞行，无法做到根据飞机状态调整；对于无减速板的无动力滑降无人机，使用俯仰角控制飞机着陆时的速度及下沉率，除此之外俯仰角还用于控制高度，两种方式耦合导致飞机对高度的跟踪效果不佳，着陆精度差。
发明内容
针对上述背景技术中存在的技术问题，本发明提出了一种基于能量的无动力滑降无人机着陆下滑道动态生成方法，其构思合理，动态生成的下滑道可有效解决无人机降落初始状态不同导致的降落状态的不确定性，可在无减速板情况下仅通过动态生成的下滑道使飞机在降落时达到理想的速度与下沉率，无需改变飞机纵向姿态控制着陆速度，对高度跟踪效果好，着陆精度高。
为解决上述技术问题，本发明提供的一种基于能量的无动力滑降无人机着陆下滑道动态生成方法，是从飞机飞至跑道后方对准跑道准备着陆开始，将飞机下滑道分为能量管理段、陡下滑段、圆弧过渡段和浅下滑段，在能量管理段实时计算飞机当前能量，根据能量剖面计算浅下滑段起始高度，然后根据剩余航程分别计算飞机的浅下滑段、圆弧过渡段和陡下滑段对应航程与高度的关系，生成下滑道高度与航程曲线。
所述基于能量的无动力滑降无人机着陆下滑道动态生成方法，其中，应用于无人机进入着陆段的情况,具体包括如下步骤：
S100、轨迹仿真确定能量剖面边界
能量剖面边界是最陡下滑能量剖面和最大升阻比下滑能量剖面构成，标称下滑能量剖面位于最陡下滑和最大升租比下滑之间，是飞机沿理想下滑道飞行得到的下滑轨迹，能量剖面是待飞距离x的函数；最陡下滑是无人机进入陡下滑段前按照不超过最大允许动压限制达到的最陡下滑轨迹；最大升阻比下滑是无人机陡下滑前保持最大升阻比飞行所得到的下滑轨迹；根据飞机沿最陡下滑轨迹与最大升阻比轨迹来计算飞机的能量剖面，能量为单位质量的动能和势能,且飞机能量的计算公式如下：
上式(1)中，Eow为能量，g为重力加速度，V为飞行速度，H为飞机高度；以待飞距离x为自变量，能量Eow为因变量，即可得到飞机的能量剖面；
着陆段能量剖面以沿跑道方向落点后方L公里开始计算，在该位置处，飞机分别以最陡飞行、最大升阻比飞行以及标称飞行，以低于陡下滑起始高度Hd结束；根据飞机以最陡飞行、最大升阻比飞行以及标称飞行时的能量以及对应待飞距离x便可得到飞机的能量剖面；
能量剖面确定后，使用二级傅里叶级数拟合能量曲线，分别得到最陡下滑、标称下滑、最大升阻比下滑能量与待飞距x的函数：
上式(2)-(4)中，Eh、Em、El分别为最陡下滑、标称下滑、最大升阻比下滑时的能量，x为飞机当前位置与落点位置的水平距离即待飞距；a0h，a1h，b0h，b1h，wh，a0m，a1m，b0m，b1m，wm，a0l，a1l，a2l，b0l，b1l，wl为曲线拟合的系数；
S200、设计标准着陆下滑道
根据已知的飞机理想的触地速度与下沉率，通过多次仿真迭代，确定标准着陆下滑道；
S300、设计最陡下滑及最大升阻比下滑着陆下滑道
分别以最陡下滑及最大升阻比下滑结束时的状态作为初始条件，使用标准下滑轨迹的陡下滑角r1以及圆弧半径R，仅通过改变浅下滑起始高度Hs，使飞机触地时满足理想速度与下沉率，分别得到最陡下滑及最大升阻比浅下滑的起始高度H1,H2；
S400、根据无人机能量实时计算浅下滑段起始高度；
无人机下滑时的能量剖面确定后，通过所述步骤S100-S300得到满足理想速度及下沉率的标称、最陡下滑和最大升阻比的着陆轨迹；无人机在实际飞行时，根据标称、最陡下滑和最大升阻比的浅下滑起始高度及能量剖面，实时计算飞机当前能量在能量剖面中的位置，通过一维线性插值方式计算当前能量下浅下滑开始时应有初始高度，进而动态生成不同的下滑轨迹；
S500、根据能量计算的浅下滑段起始高度动态生成下滑道
陡下滑角r1、圆弧半径R以及浅下滑角r2均为定值，根据步骤S400实时计算更新的浅下滑起始高度Hs，即可动态生成完整的着陆下滑道。
所述基于能量的无动力滑降无人机着陆下滑道动态生成方法，其中，所述步骤S200的具体过程为：定义浅下滑起始高度Hs，圆弧段半径R，陡下滑轨迹倾角r1；r1的选取原则是使动压在陡下滑段保持恒定；初始状态选择标称下滑结束时的状态，轨迹设计的过程如下：
S201、给定标称下滑结束时飞机状态；
S202、选取飞机的陡下滑角r1、圆弧半径R和浅下滑起始高度Hs；
S203、计算飞机的浅下滑、圆弧段和深下滑段坐标；
S204、飞机以高度跟踪进行飞行仿真；
S205、飞机落地速度下沉率若满足条件，则确定标准下滑道陡下滑角r1、圆弧半径R、浅下滑起始高度H2，结束轨迹设计；飞机落地速度下沉率若不满足条件，则返回到步骤S202继续进行；最终通过多次仿真迭代，确定一组陡下滑角r1、圆弧半径R以及浅下滑起始高度Hs作为标准着陆下滑道。
所述基于能量的无动力滑降无人机着陆下滑道动态生成方法，其中，所述步骤S300中获得最陡下滑的起始高度H1的具体过程为：
S301、给定最陡下滑结束时飞机状态；
S302、使用所述步骤S200得到的飞机的陡下滑角r1、圆弧半径R；
S303、选取一个起始高度Hs；
S304、计算飞机的浅下滑、圆弧段和深下滑段坐标；
S305、飞机以高度跟踪进行飞行仿真；
S306、飞机落地速度下沉率若满足条件，则确定最陡下滑浅下滑起始高度H1，结束轨迹设计；飞机落地速度下沉率若不满足条件，则返回到步骤S303继续进行。
所述基于能量的无动力滑降无人机着陆下滑道动态生成方法，其中，所述步骤S300中获得最大升阻比的浅下滑起始高度H2的具体过程为:
S311、给定最大升阻比下滑结束时飞机状态；
S312、使用所述步骤S200得到的飞机的陡下滑角r1、圆弧半径R；
S313、选取一个起始高度Hs；
S314、计算飞机的浅下滑、圆弧段和深下滑段坐标；
S315、飞机以高度跟踪进行飞行仿真；
S316、飞机落地速度下沉率若满足条件，则确定最陡下滑浅下滑起始高度H2，不满足条件落地速度下沉率，则返回到步骤S313继续进行。
所述基于能量的无动力滑降无人机着陆下滑道动态生成方法，其中，所述步骤S400中当前能量下浅下滑开始时应有初始高度的计算公式如下：
飞机当前能量处于标称能量与最陡下滑能量之间时：
Hs＝(H1-H2)×(Ec-Em)/(Eh-Em)+H2 (5)；
飞机当前能量处于标称能量与最大升阻比下滑能量之间时：
Hs＝(H2-H3)×(Ec-El)/(Em-El)+H3 (6)；
上式(5)-(6)中，Hs为计算的浅下滑起始高度，Ec为飞机当前能量，Em为当前待飞距离下标准下滑时应该具备的能量，Eh为当前待飞距离下，最陡下滑时飞机具备的能量；El为当前待飞距离下，最大升阻比飞行时飞机具备的能量。
所述基于能量的无动力滑降无人机着陆下滑道动态生成方法，其中：所述下滑道高度在计算时，首先根据飞机能量动态生成浅下滑段高度，再反向推导圆弧过渡段以及陡下滑段不同待飞距离x飞机应处的不同高度给定Hcmd；
所述浅下滑段高度的计算方法如下：
根据所述步骤400实时计算得到浅下滑应有起始高度Hs，可得浅下滑开始时与着陆点的距离为：Ls＝Hs/tan(r2)；
定义判断条件为：x<Ls，其中，x为无人机距离与落点的距离，即待飞距；
高度给定为：Hcmd＝tan(r2)·x。
所述基于能量的无动力滑降无人机着陆下滑道动态生成方法，其中，所述圆弧段高度的计算方法如下：
定义圆弧段起始高度为：
Hc＝R·(-cos(r1)+cos(r2))+Hs；
定义圆弧段起始待飞距为：
Lc＝R·sin(r1)-R·sin(r2)+Ls；
定义判断条件为:x<Lc；
高度给定为：
所述基于能量的无动力滑降无人机着陆下滑道动态生成方法，其中，所述陡下滑段高度的计算方法如下：
定义陡下滑起始高度为：Hd；
定义陡下滑起始待飞距为：Ld＝(Hd-Hc)/tan(r1)+Lc；
定位判断条件：x<Ld；
高度给定为：
Hcmd＝(x-Ld)/tan(r1)+Hc。
采用上述技术方案，本发明具有如下有益效果：
本发明基于能量的无动力滑降无人机着陆下滑道动态生成方法构思合理，动态生成的下滑道可有效解决无人机降落初始状态不同导致的降落状态的不确定性，可在无减速板情况下仅通过动态生成的下滑道使飞机在降落时达到理想的速度与下沉率；无需改变飞机纵向姿态控制着陆速度，对高度跟踪效果好，着陆精度高。
本发明具体优点还体现在以下几方面：
(1)动态生成着陆下滑道，通过计算飞机当前能量实时调整飞机浅下滑初始时刻高度，进而更改整个飞行高度剖面，达到对着陆速度的有效控制；
(2)将纵向通道解耦，俯仰角仅用于控制飞机高度，提升了高度的跟踪效果及着陆时的精度；
(3)动态生成的下滑道可有效解决无人机降落初始状态不同导致的降落状态的不确定性，可在无减速板情况下仅通过动态生成的下滑道使飞机在降落时达到理想的速度与下沉率。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明基于能量的无动力滑降无人机着陆下滑道动态生成方法的流程图；
图2为本发明基于能量的无动力滑降无人机着陆下滑道动态生成方法中的无人机滑降高度剖面图；
图3为本发明基于能量的无动力滑降无人机着陆下滑道动态生成方法中的无人机能量剖面图；
图4为本发明基于能量的无动力滑降无人机着陆下滑道动态生成方法中的标准下滑道设计流程图；
图5为本发明基于能量的无动力滑降无人机着陆下滑道动态生成方法中的最陡及最大升阻比下滑道设计流程图；
图6为本发明基于能量的无动力滑降无人机着陆下滑道动态生成方法中的三条下滑道高度轨迹图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
下面结合具体的实施方式对本发明做进一步的解释说明。
如图1所示,本发明基于能量的无动力滑降无人机着陆下滑道动态生成方法,是从飞机飞至跑道后方L公里对准跑道准备着陆开始，先将飞机下滑道分为能量管理段、陡下滑段、圆弧过渡段和浅下滑段，从能量管理段实时计算飞机当前能量，根据能量剖面计算浅下滑段起始高度，然后根据剩余航程分别计算飞机的浅下滑段、圆弧过渡段和陡下滑段对应航程与高度的关系，生成下滑道高度与航程曲线。
其中，上述的陡下滑段、圆弧过渡段和浅下滑段的高度均为待飞距离x的函数。陡下滑阶段开始飞机起落架下放且飞机高度快速下降，轨迹倾斜角较大。在浅下滑与陡下滑之间，有一个过渡段：圆弧过渡段。圆弧过渡段的目的是使飞机倾角逐渐减小，使轨迹平滑过渡至浅下滑段。浅下滑为飞机着陆的最后阶段，该阶段轨迹角小，可有效对飞机进行减速，该阶段越长，飞机减速效果越明显，本发明主要通过动态改变此段初始高度，用于对飞机不同的初始能量状态减速，使飞机降落时刻处于理想的速度与下沉率。
下滑道高度计算时(即根据能量确定浅下滑起始高度Hs,进而反向推导圆弧段及浅下滑段不同距离对应的高度轨迹)，首先根据飞机能量动态生成浅下滑段高度，再反向推导圆弧过渡段以及陡下滑段不同待飞距离x飞机应处的不同高度给定Hcmd。如图2所示，r1为陡下滑轨迹倾角，r2为浅下滑轨迹倾角，R为圆弧过渡段半径，这几个参数在设计阶段均已确定数值。Hs为浅下滑段起始高度，Ls为浅下滑段起始距离；Hc为圆弧过渡段起始高度，Lc为圆弧过渡段起始距离；Hd为陡下滑段起始高度，Ld为陡下滑段起始距离。
各阶段下滑道高度给定计算方法如下：
1)浅下滑段高度计算：
根据下述步骤S400实时计算得到浅下滑应有起始高度Hs，可得浅下滑开始时与着陆点的距离：
Ls＝Hs/tan(r2)；
判断条件：x<Ls，其中，x为无人机距离与落点的距离，即待飞距
高度给定：
Hcmd＝tan(r2)·x；
2)圆弧段高度计算：
圆弧段起始高度：
Hc＝R·(-cos(r1)+cos(r2))+Hs；
圆弧段起始待飞距：
Lc＝R·sin(r1)-R·sin(r2)+Ls；
判断条件:x<Lc；
高度给定：
3)陡下滑段计算：
陡下滑起始高度：Hd；
陡下滑起始待飞距：Ld＝(Hd-Hc)/tan(r1)+Lc；
判断条件：x<Ld；
高度给定：
Hcmd＝(x-Ld)/tan(r1)+Hc。
如图1所示，本实施例提供的基于能量的无动力滑降无人机着陆下滑道动态生成方法，主要应用于无人机进入着陆段(着陆点后L公里开始，整个流程都是着陆段，如图2所示)的情况，包括如下步骤：
S100、轨迹仿真确定能量剖面边界
能量剖面边界是最陡下滑能量剖面和最大升阻比下滑能量剖面构成，标称下滑能量剖面位于最陡下滑和最大升租比下滑之间，是飞机沿理想下滑道飞行得到的下滑轨迹，纵向制导(是指使飞机沿着预定的高度速度飞行)的目的就是使飞机状态尽量保持在标称下滑能量剖面，能量剖面是待飞距离x的函数；最陡下滑是无人机进入陡下滑段前按照不超过最大允许动压限制达到的最陡下滑轨迹；最大升阻比下滑是无人机陡下滑前保持最大升阻比飞行所得到的下滑轨迹；根据飞机沿最陡下滑轨迹与最大升阻比轨迹计算得到飞机的能量剖面，能量为单位质量的动能和势能,且飞机能量的计算公式如下：
上式(1)中，Eow为能量，g为重力加速度，V为飞行速度，H为飞机高度。
着陆段能量剖面以沿跑道方向落点后方L公里开始计算，本发明以跑道落点后方40km处开始，即L＝40km。在该位置处，飞机分别以最陡飞行，最大升阻比飞行以及标称飞行，以低于陡下滑起始高度Hd结束。根据三组飞行时的能量以及对应待飞距离x便可得到飞机的能量剖面，如图3所示。
能量剖面确定后，使用二级傅里叶级数拟合能量曲线，分别得到最陡下滑、标称下滑、最大升阻比下滑能量与待飞距x的函数。
上式(2)-(4)中，Eh，Em，El分别为最陡下滑、标称下滑、最大升阻比下滑时的能量，x为飞机当前位置与落点位置的水平距离即待飞距。a0h，a1h，b0h，b1h，wh，a0m，a1m，b0m，b1m，wm，a0l，a1l，a2l，b0l，b1l，wl为曲线拟合的系数。
S200、设计标准着陆下滑道
根据已知的飞机理想的触地速度与下沉率，通过多次仿真迭代，确定标准着陆下滑道；
着陆轨迹设计是一个两端约束问题，已知起始及结束条件，在高度剖面上规划速度剖面，使其满足着陆时的约束条件。
由于飞机触地速度与下沉率在设计阶段已知，由它便可确定浅下滑轨迹倾角r2；因此主要的设计参数为浅下滑起始高度Hs，圆弧段半径R，陡下滑轨迹倾角r1；r1的选取原则是使动压在陡下滑段保持恒定；初始状态选择标称下滑结束时的状态，轨迹设计的过程如下(图4所示)：
S201、给定标称下滑结束时飞机状态(高度、速度、姿态)；
S202、选取飞机的陡下滑角r1、圆弧半径R和浅下滑起始高度Hs；
S203、计算飞机的浅下滑、圆弧段和深下滑段坐标；
S204、飞机以高度跟踪进行飞行仿真；
S205、飞机落地速度下沉率若满足条件，则确定标准下滑道陡下滑角r1、圆弧半径R、浅下滑起始高度H2，结束轨迹设计；飞机落地速度下沉率若不满足条件，则返回到步骤S202继续进行；最终通过多次仿真迭代，确定一组陡下滑角r1、圆弧半径R以及浅下滑起始高度Hs作为标准着陆下滑道。
S300、设计最陡下滑及最大升阻比下滑着陆下滑道
分别以最陡下滑及最大升阻比下滑结束时的状态作为初始条件，使用标准下滑轨迹的陡下滑角r1以及圆弧半径R，仅通过改变浅下滑起始高度Hs，使飞机触地时满足理想速度与下沉率，分别得到最陡下滑及最大升阻比浅下滑的起始高度H1,H2，具体过程如图5所示；
其中，获得最陡下滑的起始高度H1的具体过程为：
S301、给定最陡下滑结束时飞机状态(高度、速度、姿态)；
S302、使用上述步骤S200得到的飞机的陡下滑角r1、圆弧半径R；
S303、选取一个起始高度Hs；
S304、计算飞机的浅下滑、圆弧段和深下滑段坐标；
S305、飞机以高度跟踪进行飞行仿真；
S306、飞机落地速度下沉率若满足条件，则确定最陡下滑浅下滑起始高度H1，结束轨迹设计；飞机落地速度下沉率若不满足条件，则返回到步骤S303继续进行。
同时，获得最大升阻比浅下滑的起始高度H2的具体过程为：
S311、给定最大升阻比下滑结束时飞机状态(高度、速度、姿态)；
S312、使用上述步骤S200得到的飞机的陡下滑角r1、圆弧半径R；
S313、选取一个起始高度Hs；
S314、计算飞机的浅下滑、圆弧段和深下滑段坐标；
S315、飞机以高度跟踪进行飞行仿真；
S316、飞机落地速度下沉率若满足条件，则确定最陡下滑浅下滑起始高度H2，不满足条件落地速度下沉率，则返回到步骤S313继续进行。
S400、根据无人机能量实时计算浅下滑段起始高度
无人机下滑时的能量剖面确定后，通过前面步骤S100-S300得到标称、最陡下滑、最大升阻比的三条着陆轨迹。这三条着陆轨迹均满足理想速度及下沉率，且陡下滑角r1、圆弧半径R以及浅下滑角r2均相同，仅浅下滑起始高度Hs不同。三条着陆轨迹的浅下滑(即标称、最陡下滑和最大升阻比的浅下滑)的起始高度分别为H1，H2，H3，如图6所示。H1为飞机以最陡下滑飞行的浅下滑起始高度；H2为飞机机以标准下滑道飞行浅下滑起始高度；H3为飞机以最大升阻比飞行时浅下滑起始高度。无人机在实际飞行时，可根据三条着陆轨迹的起始高度(即标称、最陡下滑和最大升阻比轨迹的起始高度)及能量剖面，实时计算飞机当前能量在能量剖面中的位置，通过一维线性插值方式计算当前能量下浅下滑开始时应有初始高度，进而动态生成不同的下滑轨迹。
当前能量下浅下滑开始时应有初始高度的计算公式如下：
飞机当前能量处于标称能量与最陡下滑能量之间时：
Hc＝(H1-H2)×(Ec-Em)/(Eh-Em)+H2 (5)；
飞机当前能量处于标称能量与最大升阻比下滑能量之间时：
Hc＝(H2-H3)×(Ec-El)/(Em-El)+H3 (6)；
式中，Hc为计算的浅下滑起始高度，Ec为飞机当前能量，Em为当前待飞距离下标准下滑时应该具备的能量，Eh为当前待飞距离下，最陡下滑时飞机具备的能量；El为当前待飞距离下，最大升阻比飞行时飞机具备的能量。
S500、根据能量计算的浅下滑段起始高度动态生成下滑道
陡下滑角r1、圆弧半径R以及浅下滑角r2均为定值，再根据步骤S400实时计算更新的浅下滑起始高度Hs，即可动态生成完整的着陆下滑道。
本发明动态生成的下滑道可有效解决无人机降落初始状态不同导致的降落状态的不确定性，可在无减速板情况下仅通过动态生成的下滑道使飞机在降落时达到理想的速度与下沉率；无需改变飞机纵向姿态控制着陆速度，对高度跟踪效果好，着陆精度高。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (9)

1.一种基于能量的无动力滑降无人机着陆下滑道动态生成方法，其特征在于：从飞机飞至跑道后方对准跑道准备着陆开始，将飞机下滑道分为能量管理段、陡下滑段、圆弧过渡段和浅下滑段，在能量管理段实时计算飞机当前能量，根据能量剖面计算浅下滑段起始高度，然后根据剩余航程分别计算飞机的浅下滑段、圆弧过渡段和陡下滑段对应航程与高度的关系，生成下滑道高度与航程曲线。
2.如权利要求1所述的基于能量的无动力滑降无人机着陆下滑道动态生成方法，其特征在于，应用于无人机进入着陆段的情况,具体包括如下步骤：
S100、轨迹仿真确定能量剖面边界
能量剖面边界是最陡下滑能量剖面和最大升阻比下滑能量剖面构成，标称下滑能量剖面位于最陡下滑和最大升租比下滑之间，是飞机沿理想下滑道飞行得到的下滑轨迹，能量剖面是待飞距离x的函数；最陡下滑是无人机进入陡下滑段前按照不超过最大允许动压限制达到的最陡下滑轨迹；最大升阻比下滑是无人机陡下滑前保持最大升阻比飞行所得到的下滑轨迹；根据飞机沿最陡下滑轨迹与最大升阻比轨迹来计算飞机的能量剖面，能量为单位质量的动能和势能,且飞机能量的计算公式如下：
上式(1)中，Eow为能量，g为重力加速度，V为飞行速度，H为飞机高度；以待飞距离x为自变量，能量Eow为因变量，即可得到飞机的能量剖面；
着陆段能量剖面以沿跑道方向落点后方L公里开始计算，在该位置处，飞机分别以最陡飞行、最大升阻比飞行以及标称飞行，以低于陡下滑起始高度Hd结束；根据飞机以最陡飞行、最大升阻比飞行以及标称飞行时的能量以及对应待飞距离x便可得到飞机的能量剖面；
能量剖面确定后，使用二级傅里叶级数拟合能量曲线，分别得到最陡下滑、标称下滑、最大升阻比下滑能量与待飞距x的函数：
上式(2)-(4)中，Eh、Em、El分别为最陡下滑、标称下滑、最大升阻比下滑时的能量，x为飞机当前位置与落点位置的水平距离即待飞距；a0h，a1h，b0h，b1h，wh，a0m，a1m，b0m，b1m，wm，a0l，a1l，a2l，b0l，b1l，wl为曲线拟合的系数；
S200、设计标准着陆下滑道
根据已知的飞机理想的触地速度与下沉率，通过多次仿真迭代，确定标准着陆下滑道；
S300、设计最陡下滑及最大升阻比下滑着陆下滑道
分别以最陡下滑及最大升阻比下滑结束时的状态作为初始条件，使用标准下滑轨迹的陡下滑角r1以及圆弧半径R，仅通过改变浅下滑起始高度Hs，使飞机触地时满足理想速度与下沉率，分别得到最陡下滑及最大升阻比浅下滑的起始高度H1,H2；
S400、根据无人机能量实时计算浅下滑段起始高度；
无人机下滑时的能量剖面确定后，通过所述步骤S100-S300得到满足理想速度及下沉率的标称、最陡下滑和最大升阻比的着陆轨迹；无人机在实际飞行时，根据标称、最陡下滑和最大升阻比的浅下滑起始高度及能量剖面，实时计算飞机当前能量在能量剖面中的位置，通过一维线性插值方式计算当前能量下浅下滑开始时应有初始高度，进而动态生成不同的下滑轨迹；
S500、根据能量计算的浅下滑段起始高度动态生成下滑道
陡下滑角r1、圆弧半径R以及浅下滑角r2均为定值，根据步骤S400实时计算更新的浅下滑起始高度Hs，即可动态生成完整的着陆下滑道。
3.如权利要求2所述的基于能量的无动力滑降无人机着陆下滑道动态生成方法，其特征在于，所述步骤S200的具体过程为：定义浅下滑起始高度Hs，圆弧段半径R，陡下滑轨迹倾角r1；r1的选取原则是使动压在陡下滑段保持恒定；初始状态选择标称下滑结束时的状态，轨迹设计的过程如下：
S201、给定标称下滑结束时飞机状态；
S202、选取飞机的陡下滑角r1、圆弧半径R和浅下滑起始高度Hs；
S203、计算飞机的浅下滑、圆弧段和深下滑段坐标；
S204、飞机以高度跟踪进行飞行仿真；
S205、飞机落地速度下沉率若满足条件，则确定标准下滑道陡下滑角r1、圆弧半径R、浅下滑起始高度H2，结束轨迹设计；飞机落地速度下沉率若不满足条件，则返回到步骤S202继续进行；最终通过多次仿真迭代，确定一组陡下滑角r1、圆弧半径R以及浅下滑起始高度Hs作为标准着陆下滑道。
4.如权利要求2所述的基于能量的无动力滑降无人机着陆下滑道动态生成方法，其特征在于，所述步骤S300中获得最陡下滑的起始高度H1的具体过程为：
S301、给定最陡下滑结束时飞机状态；
S302、使用所述步骤S200得到的飞机的陡下滑角r1、圆弧半径R；
S303、选取一个起始高度Hs；
S304、计算飞机的浅下滑、圆弧段和深下滑段坐标；
S305、飞机以高度跟踪进行飞行仿真；
S306、飞机落地速度下沉率若满足条件，则确定最陡下滑浅下滑起始高度H1，结束轨迹设计；飞机落地速度下沉率若不满足条件，则返回到步骤S303继续进行。
5.如权利要求2所述的基于能量的无动力滑降无人机着陆下滑道动态生成方法，其特征在于，所述步骤S300中获得最大升阻比的浅下滑起始高度H2的具体过程为:
S311、给定最大升阻比下滑结束时飞机状态；
S312、使用所述步骤S200得到的飞机的陡下滑角r1、圆弧半径R；
S313、选取一个起始高度Hs；
S314、计算飞机的浅下滑、圆弧段和深下滑段坐标；
S315、飞机以高度跟踪进行飞行仿真；
S316、飞机落地速度下沉率若满足条件，则确定最陡下滑浅下滑起始高度H2，不满足条件落地速度下沉率，则返回到步骤S313继续进行。
6.如权利要求2所述的基于能量的无动力滑降无人机着陆下滑道动态生成方法，其特征在于，所述步骤S400中当前能量下浅下滑开始时应有初始高度的计算公式如下：
飞机当前能量处于标称能量与最陡下滑能量之间时：
Hs＝(H1-H2)×(Ec-Em)/(Eh-Em)+H2 (5)；
飞机当前能量处于标称能量与最大升阻比下滑能量之间时：
Hs＝(H2-H3)×(Ec-El)/(Em-El)+H3 (6)；
上式(5)-(6)中，Hs为计算的浅下滑起始高度，Ec为飞机当前能量，Em为当前待飞距离下标准下滑时应该具备的能量，Eh为当前待飞距离下，最陡下滑时飞机具备的能量；El为当前待飞距离下，最大升阻比飞行时飞机具备的能量。
7.如权利要求2所述的基于能量的无动力滑降无人机着陆下滑道动态生成方法，其特征在于：所述下滑道高度在计算时，首先根据飞机能量动态生成浅下滑段高度，再反向推导圆弧过渡段以及陡下滑段不同待飞距离x飞机应处的不同高度给定Hcmd；
所述浅下滑段高度的计算方法如下：
根据所述步骤400实时计算得到浅下滑应有起始高度Hs，可得浅下滑开始时与着陆点的距离为：Ls＝Hs/tan(r2)；
定义判断条件为：x<Ls，其中，x为无人机距离与落点的距离，即待飞距；
高度给定为：Hcmd＝tan(r2)·x。
8.如权利要求2所述的基于能量的无动力滑降无人机着陆下滑道动态生成方法，其特征在于，所述圆弧段高度的计算方法如下：
定义圆弧段起始高度为：
Hc＝R·(-cos(r1)+cos(r2))+Hs；
定义圆弧段起始待飞距为：
Lc＝R·sin(r1)-R·sin(r2)+Ls；
定义判断条件为:x<Lc；
高度给定为：
9.如权利要求2所述的基于能量的无动力滑降无人机着陆下滑道动态生成方法，其特征在于，所述陡下滑段高度的计算方法如下：
定义陡下滑起始高度为：Hd；
定义陡下滑起始待飞距为：Ld＝(Hd-Hc)/tan(r1)+Lc；
定位判断条件：x<Ld；
高度给定为：
Hcmd＝(x-Ld)/tan(r1)+Hc。
CN202310850244.8A 2023-07-12 2023-07-12 基于能量的无动力滑降无人机着陆下滑道动态生成方法 Active CN116991171B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310850244.8A CN116991171B (zh) 2023-07-12 2023-07-12 基于能量的无动力滑降无人机着陆下滑道动态生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310850244.8A CN116991171B (zh) 2023-07-12 2023-07-12 基于能量的无动力滑降无人机着陆下滑道动态生成方法

Publications (2)

Publication Number Publication Date
CN116991171A true CN116991171A (zh) 2023-11-03
CN116991171B CN116991171B (zh) 2024-04-16

Family

ID=88529300

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310850244.8A Active CN116991171B (zh) 2023-07-12 2023-07-12 基于能量的无动力滑降无人机着陆下滑道动态生成方法

Country Status (1)

Country Link
CN (1) CN116991171B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107941087A (zh) * 2017-10-18 2018-04-20 北京航空航天大学 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法
CN109253730A (zh) * 2018-11-07 2019-01-22 广州大学 可重复使用运载器末端能量管理段的三维轨迹在线规划方法及系统
US20190371187A1 (en) * 2018-06-05 2019-12-05 Honeywell International Inc. Methods and systems for stabilized approach energy management
CN113534847A (zh) * 2021-08-23 2021-10-22 哈尔滨工业大学 一种有动力可重复使用飞行器能量管理轨迹设计方法
US20220121223A1 (en) * 2019-04-25 2022-04-21 Aerovironment, Inc. Methods of Climb and Glide Operations of a High Altitude Long Endurance Aircraft
CN114647249A (zh) * 2020-12-18 2022-06-21 英西图公司(波音公司子公司) 垂直降落载具的降落
EP4050585A1 (en) * 2021-02-24 2022-08-31 Honeywell International Inc. Methods, and systems for enabling effective modeling of landing gear for energy management

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107941087A (zh) * 2017-10-18 2018-04-20 北京航空航天大学 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法
US20190371187A1 (en) * 2018-06-05 2019-12-05 Honeywell International Inc. Methods and systems for stabilized approach energy management
CN109253730A (zh) * 2018-11-07 2019-01-22 广州大学 可重复使用运载器末端能量管理段的三维轨迹在线规划方法及系统
US20220121223A1 (en) * 2019-04-25 2022-04-21 Aerovironment, Inc. Methods of Climb and Glide Operations of a High Altitude Long Endurance Aircraft
CN114647249A (zh) * 2020-12-18 2022-06-21 英西图公司(波音公司子公司) 垂直降落载具的降落
EP4050585A1 (en) * 2021-02-24 2022-08-31 Honeywell International Inc. Methods, and systems for enabling effective modeling of landing gear for energy management
CN113534847A (zh) * 2021-08-23 2021-10-22 哈尔滨工业大学 一种有动力可重复使用飞行器能量管理轨迹设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
宋辉 等: "《大展弦比无人机无动力返场控制技术》", 《科学技术创新》, no. 18, pages 1 - 5 *

Also Published As

Publication number Publication date
CN116991171B (zh) 2024-04-16

Similar Documents

Publication Publication Date Title
CN110908396B (zh) 可重复使用运载器全阶段再入返回制导方法
CN110497943B (zh) 一种基于强化学习的城轨列车节能运行策略在线优化方法
CN102927851B (zh) 一种基于轨迹在线规划的末制导方法
CA2674732C (en) Implementing continuous descent approaches for maximum predictability in aircraft
US8660722B2 (en) Method for optimizing aircraft landing on a runway
CN103983143A (zh) 包含速度过程约束和多终端约束的空地导弹投放下滑段制导方法
CN109240335B (zh) 一种空天飞行器进场着陆制导方法
CN111580547B (zh) 一种高超声速飞行器编队控制方法
CN106021628A (zh) 一种运载火箭垂直返回弹道设计方法
CA2599474A1 (en) Method for implementing required navigational performance procedures
CN101893892B (zh) 一种无人机自动伞降回收控制方法
CN105159308A (zh) 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法
CN110031003B (zh) 一种火箭上面级最优可达轨道快速规划计算方法
CN109240323A (zh) 一种实时解析构造的空天飞行器再入制导方法
CN109657256A (zh) 一种高精度弹道式再入标称返回轨道仿真方法
CN105717937A (zh) 用飞行器航空电子设备自动控制飞行器的下降阶段的方法
CN108298110A (zh) 一种两级入轨空天飞行器上升段轨迹及设计方法
CN116991171B (zh) 基于能量的无动力滑降无人机着陆下滑道动态生成方法
CN104656659B (zh) 舰载机滑跃起飞自动飞行控制方法
CN104656661B (zh) 一种公务机下降着陆控制方法
CN106873615A (zh) 应急返场着陆速度指令集设计方法
CN110262558B (zh) 一种无人机定点着陆的控制方法
Wu et al. Optimum flight trajectory guidance based on total energy control of aircraft
CN116301028B (zh) 基于吸气式高超声速平台的多约束在线飞行轨迹规划中段导引方法
CN113190024B (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