CN111650958A - 一种固定翼无人机起飞段切入航路点的在线路径规划方法 - Google Patents

一种固定翼无人机起飞段切入航路点的在线路径规划方法 Download PDF

Info

Publication number
CN111650958A
CN111650958A CN201911287849.0A CN201911287849A CN111650958A CN 111650958 A CN111650958 A CN 111650958A CN 201911287849 A CN201911287849 A CN 201911287849A CN 111650958 A CN111650958 A CN 111650958A
Authority
CN
China
Prior art keywords
path
waypoint
aerial vehicle
unmanned aerial
arc
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
CN201911287849.0A
Other languages
English (en)
Other versions
CN111650958B (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.)
Hubei Institute Of Aerospacecraft
Original Assignee
Hubei Institute Of Aerospacecraft
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 Hubei Institute Of Aerospacecraft filed Critical Hubei Institute Of Aerospacecraft
Priority to CN201911287849.0A priority Critical patent/CN111650958B/zh
Publication of CN111650958A publication Critical patent/CN111650958A/zh
Application granted granted Critical
Publication of CN111650958B publication Critical patent/CN111650958B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft

Abstract

本发明公开了固定翼无人机起飞段切入航路点的在线路径规划方法:在起飞段,无人机通过导航系统器件和卫星接收机测量自身的飞行姿态、速度和地理位置三维坐标信息;将飞行器自身三维位置高度和飞行速度方向等信息与期望第一个航路点三维位置高度和期望飞行速度方向等信息进行比较,以便于规划爬升和转弯路径;纵向上规划出到第一个航路点高度的可行爬升路径;水平方向上规划出满足第一个航路点期望方向的转弯路径。针对固定翼无人机起飞段,增强了无人机自主化程度,设计合理可行并且计算量小、在线执行的从起飞点切入第一个航路点的飞行路径规划策略,使无人机在高度和飞行方向上都能够平滑的切入到第一个预设航路点,提高系统的可靠性和实用性。

Description

一种固定翼无人机起飞段切入航路点的在线路径规划方法
技术领域
本发明属于飞行器导航、制导和控制的技术领域,具体涉及一种固定翼无人机起飞段切入航路点的在线路径规划方法。
背景技术
无人机具有体积小、机动性高、能有效减少驾驶人员伤亡及全寿命周期成本低等方面的优势,可以替代有人驾驶飞机执行各种枯燥的、危险的、恶劣的任务,例如远程监视、侦查和打击等军事领域,以及航拍、植保、快递运输、灾难救援和测绘等民事领域,具有广阔的军用和民用前景。
随着信息化、自动化、智能化技术的不断发展,无人机技术也发生了很大的改变,无人飞行器自主化程度要求的不断提高。目前无人机的任务规划系统为了圆满完成飞行任务,一般预先设置飞行途中必须经过的一系列关键目标飞行航路点,这些航路点通过一系列直线和圆弧连接而组合成飞行航路,这样对于每一个航点就包含了三维位置和期望飞行方向的信息。通常由于起飞场地的限制,如起飞位置和跑道方向或无人机弹射起飞方向与第一个航点的位置和方向不一致;又由于无人机起飞条件受到现场风向影响较大,通常期望固定翼无人机初始逆风起飞来增加相对空气的速度来提高升力,这样起飞点的位置和方向与第一个航点的位置和方向也不一致。同时,对于固定翼无人机或垂直起降复合式固定翼无人机来说,由于每次起飞位置和第一个航路点相对高度位置不确定,同时起飞的机头方向和预设的第一个飞行航点的期望飞行方向也不能每次预先固定,传统的起飞策略都需要有经验的地面站操作员规划出一种合理的起飞路线策略,从三维位置和方向上导引飞行器起飞到第一个航路点,倘若每次起飞前都要由人工规划出飞行器切入航路点的飞行路径,会大大降低飞行器的自主化程度。
固定翼无人机在起飞阶段,飞行器的爬升角同样也不会设置过大,这是为了防止俯仰角过大使飞行器产生失速,失速会导致机翼升力不足而坠机,这对于如何规划起飞阶段的高度策略也是一种挑战。
发明内容
针对固定翼无人机起飞段路径规划这一特殊飞行阶段,如何增强无人机自主化程度,能够在线规划出从起飞点切入第一个航路点的合理可行路径亟待解决。
实现本发明目的所采用的技术方案具体如下:
一种固定翼无人机起飞段切入航路点的在线路径规划方法,包括如下步骤:
S1、在起飞段,无人机通过导航系统(如惯性传感器)器件和卫星接收机测量自身的飞行姿态、速度和地理位置三维坐标信息;
S2、将飞行器自身三维位置高度和飞行速度方向等信息与期望第一个航路点三维位置高度和期望飞行速度方向等信息进行比较,以便于规划爬升和转弯路径;
S3、纵向上规划出到第一个航路点高度的可行爬升路径;
S4、水平方向上规划出满足第一个航路点期望方向的转弯路径。
进一步地,所述步骤S1中无人机测量自身的飞行姿态、速度和地理位置三维坐标信息的具体方法为:无人机相对地面惯性坐标系的姿态角信息由惯性传感器(陀螺仪、加速度计、惯性测量单元)来估计,而地理位置三维坐标、水平飞行速度和速度方向等信息由卫星全球定位导航传感器 (GPS)来估计。
再进一步地,所述步骤S3中纵向上规划出到第一个航路点高度的可行爬升路径时,充分考虑无人机爬升角过大可能会导致失速的的工程实践风险,设计随实际水平待飞距离变化的直线爬升和盘旋爬升两种路径规划策略。
更进一步地,所述步骤S4中水平方向上规划出满足第一个航路点期望方向的转弯路径时,是一种飞行路径最短的优化路径。
优选地,所述飞行路径最短的优化路径采用基于Dubins路径的规划转弯策略予以规划。
具体地,所述纵向随实际水平待飞距离变化的爬升路径规划策略具体为:
设起飞点位置向量
Figure RE-GDA0002615882700000031
处的当前高度为H1,由于第一个期望航路点A处的期望高度为HA,为了让飞机有一个较为平稳的爬升过程,高度指令一般以直线斜坡信号的方式给出,即按照起飞方向纵向平面内选择一个与航路点A高度相同的中间航路点p2,其位置向量为
Figure RE-GDA00026158827000000314
其高度H2与A处的高度 HA相同,则
Figure RE-GDA0002615882700000032
Figure RE-GDA0002615882700000033
的期望高度偏差ΔH12
ΔH12=H2-H1 (1)
定义
Figure RE-GDA0002615882700000034
Figure RE-GDA0002615882700000035
的位置向量为
Figure RE-GDA0002615882700000036
定义
Figure RE-GDA0002615882700000037
Figure RE-GDA0002615882700000038
的水平投影距离为|p1p2|,这也是实际水平待飞距离。爬升指令即根据飞机的当前位置,计算出当前的高度指令,通常给定飞机一个固定的爬升角,按照固定的爬升角飞行到达期望的高度。定义期望爬升角即期望航迹倾斜角γ
Figure RE-GDA0002615882700000039
Figure RE-GDA00026158827000000310
Figure RE-GDA00026158827000000311
的水平投影距离|p1p2|即实际水平待飞距离较短时,同时高度偏差ΔH12较大时,飞机的期望航迹倾斜角γ会变大,当该角过大时,俯仰角过大使飞行器产生失速,失速会导致机翼升力不足而坠机。因此这种直线爬升的方式仅适合γ较小时的情况,考虑到无人飞机的失速迎角一般在15 度左右,考虑预留一定的裕度,以期望航迹倾斜角10°为分界线(如果航迹倾斜角大于10°则飞行器产生失速):
1)当
Figure RE-GDA00026158827000000312
时,即γ≤10°时采用直线爬升的方式;
直线爬升的高度指令Hg如下:
Figure RE-GDA00026158827000000313
式中p代表无人机当前的位置,定义飞行点到当前无人机的水平投影距离为|p1p|,变量tanγ应该满足tanγ<tan10°。
2)当
Figure RE-GDA0002615882700000041
时,即γ>10°时采用直线爬升和盘旋爬升结合的方式,直线爬升的期望航迹倾斜角为10°,先按期望航迹倾斜角10°直线爬升到中间航路点p2下方的水平投影p3点处,再按盘旋的半径Rpan盘旋爬升到中间航路点p2,盘旋的半径Rpan为预先设定,盘旋爬升时的飞行器滚转角φ满足飞行时最大滚转角φmax的约束,即φ≤φmax
具体地,所述飞行路径最短的优化路径采用基于Dubins路径的规划转弯策略予以规划的策略具体为:
从一个航点配置(位置和航向)过渡到另一个航点配置时,两个不同配置之间的时间最优路径由一段圆弧跟着一段直线,结尾再由一段圆弧到达最总配置,这种圆弧—直线—圆弧的路径对于固定高度和恒定速度的情况是时间最优路径;
假设无人机在执行水平飞行时飞控系统能够保持地速V绝对值不变,定义圆弧路径的半径用R表示,半径与无人机的转弯半径一样大;对于任何给定的开始和结束配置,由一段圆弧接一段直线再加上一段圆弧的路径有四种情况,
第一种情况:右圆弧—直线—右圆弧;
第二种情况:右圆弧—直线—左圆弧;
第三种情况:左圆弧—直线—右圆弧;
第四种情况:左圆弧—直线—左圆弧;
为了确定路径最短的路径,需要计算四种情况的路径长度;给定无人机的位置向量信息
Figure RE-GDA0002615882700000042
航向χ和转弯半径R,右转和左转圆弧的中心位置向量
Figure RE-GDA0002615882700000043
Figure RE-GDA0002615882700000044
分别由下式给出:
Figure RE-GDA0002615882700000045
Figure RE-GDA0002615882700000051
式中符号T代表向量的转置;
为了计算不同情况轨迹的长度,需要分别针对顺时针圆和逆时针圆确定角距离,假设角度θ1和θ2分别在0和2π之间,对于顺时针圆,θ1和θ2之间的角距离如下:
21|=<2π+θ21> (6)
式中,符号<θ>代表角θ对2π取模,即
<θ>=θmod 2π (7)
同样,对于逆时针圆,θ1和θ2之间的角距离如下:
21|=<2π-θ21> (8)
假设Dubins路径开始圆心
Figure RE-GDA0002615882700000052
下标为s,结束圆心
Figure RE-GDA0002615882700000053
下标为e,开始圆心
Figure RE-GDA0002615882700000054
与结束圆心
Figure RE-GDA0002615882700000055
连线与正北方向的夹角定义为
Figure RE-GDA0002615882700000056
则开始圆心到结束圆心的距离l如下:
Figure RE-GDA0002615882700000057
因此,四种情况的路径长度计算公式分别如下:
第一种情况中右圆弧—直线—右圆弧的路径长度L1如下:
Figure RE-GDA0002615882700000058
第二种情况中右圆弧—直线—左圆弧的路径长度L2如下:
Figure RE-GDA0002615882700000059
其中,
Figure RE-GDA0002615882700000061
第三种情况中左圆弧—直线—右圆弧的路径长度L3如下:
Figure RE-GDA0002615882700000062
其中,
Figure RE-GDA0002615882700000063
第四种情况中左圆弧—直线—左圆弧的路径长度L4如下:
Figure RE-GDA0002615882700000064
具体地,所述无人机从起飞点位置P1纵向随实际水平待飞距离变化的爬升路径规划策略爬升到中间航路点p2后,按照所述飞行路径最短的优化路径采用基于Dubins路径的规划转弯策略予以规划的策略具体为:
从一个航点配置(位置和航向)过渡到另一个航点配置时,两个不同配置之间的时间最优路径由一段圆弧跟着一段直线,结尾再由一段圆弧到达最总配置,这种圆弧—直线—圆弧的路径对于固定高度和恒定速度的情况是时间最优路径;
假设无人机在执行水平飞行时飞控系统能够保持地速V绝对值不变,定义圆弧路径的半径用R表示,半径与无人机的转弯半径一样大;对于任何给定的开始和结束配置,航点连接时由一段圆弧接一段直线再加上一段圆弧的路径有四种情况,第一种情况:右圆弧—直线—右圆弧;
第二种情况:右圆弧—直线—左圆弧;
第三种情况:左圆弧—直线—右圆弧;
第四种情况:左圆弧—直线—左圆弧;
为了确定路径最短的路径,需要计算四种情况的路径长度;给定无人机的位置向量信息
Figure RE-GDA0002615882700000071
航向χ和转弯半径R,右转和左转圆弧的中心位置向量
Figure RE-GDA0002615882700000072
Figure RE-GDA0002615882700000073
分别由下式给出:
Figure RE-GDA0002615882700000074
Figure RE-GDA0002615882700000075
式中符号T代表向量的转置;
为了计算不同情况轨迹的长度,需要分别针对顺时针圆和逆时针圆确定角距离,假设角度θ1和θ2分别在0和2π之间,对于顺时针圆,θ1和θ2之间的角距离如下:
21|=<2π+θ21> (6)
式中,符号<θ>代表角θ对2π取模,即
<θ>=θmod 2π (7)
同样,对于逆时针圆,θ1和θ2之间的角距离如下:
21|=<2π-θ21> (8)
假设Dubins路径开始圆心
Figure RE-GDA0002615882700000076
下标为s,结束圆心
Figure RE-GDA0002615882700000077
下标为e,开始圆心
Figure RE-GDA0002615882700000078
与结束圆心
Figure RE-GDA0002615882700000079
连线与正北方向的夹角定义为
Figure RE-GDA00026158827000000710
则开始圆心到结束圆心的距离l如下:
Figure RE-GDA00026158827000000711
因此,四种情况的路径长度计算公式分别如下:
第一种情况中右圆弧—直线—右圆弧的路径长度L1如下:
Figure RE-GDA0002615882700000081
第二种情况中右圆弧—直线—左圆弧的路径长度L2如下:
Figure RE-GDA0002615882700000082
其中,
Figure RE-GDA0002615882700000083
第三种情况中左圆弧—直线—右圆弧的路径长度L3如下:
Figure RE-GDA0002615882700000084
其中,
Figure RE-GDA0002615882700000085
第四种情况中左圆弧—直线—左圆弧的路径长度L4如下:
Figure RE-GDA0002615882700000086
当无人机到达中间航路点p2后,按照以上四种情况分别计算路径长度,由于仅需要位置信息,计算量很小,可以在线完成路径规划,最终选择一条在水平面路径长度最短,并且满足从当前飞行位置和飞行方向到第一个航路点A期望方向的转弯路径,再由无人机制导控制系统来导引跟踪这条规划好的路径。
本发明的固定翼无人机起飞段切入航路点的在线路径规划方法,在起飞段,飞行器通过导航系统如惯性传感器器件和卫星接收机等设备测量自身的飞行姿态、速度和地理位置三维坐标等信息,将飞行器自身三维位置和飞行方向信息与期望第一个航路点三维位置和期望飞行方向比较,纵向上规划出到第一个航路点高度的可行爬升路径,充分考虑了无人机爬升角过大可能会导致失速的的工程实践风险,设计了随实际水平待飞距离变化的直线爬升和盘旋爬升两种路径规划策略,满足了无人机起飞过程中爬升路径不会过于陡峭的要求,提高了系统的可靠性;设计的纵向上到第一个航路点水平面的转弯路径,采用基于Dubins路径的规划转弯策略,是一种飞行路径最短的优化路径,同时通过待飞航点间“曲线圆弧—直线—曲线圆弧”的平滑路径连接,满足了使无人机的飞行方向往期望飞行方向平滑过渡的要求,提高了系统的实用性。
总之,针对固定翼无人机起飞段,增强了无人机自主化程度,设计合理可行并且计算量小、在线执行的从起飞点切入第一个航路点的飞行路径规划策略,使无人机在高度和飞行方向上都能够平滑的切入到第一个预设航路点,为后续无人机可以按照任务航路点飞行,并且为飞行路径的制导跟踪系统提供参考路径,提高系统的可靠性和实用性。
附图说明
图1为本发明无人机机体、期望飞行航路点、提出的在线规划从起飞段切入航路点的路径和三维北东地坐标系之间关系图;
图2为本发明提出的在线规划从起飞点切入第一个航路点的合理可行路径方案流程图;
图3为本发明提出的纵向随实际水平待飞距离变化的爬升路径规划策略示意图;
图4为本发明提出的水平方向上规划出满足第一个航路点期望方向的转弯路径策略示意图,图中情况1即第一种情况、情况2即第二种情况、情况3即第三种情况、情况4即第四种情况。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及具体实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明无人机机体、期望飞行航路点、提出的在线规划从起飞段切入航路点的路径和三维北东地坐标系之间关系图如图1所示。具体步骤如下:
无人机相对地面惯性坐标系的姿态角信息通常由惯性传感器(陀螺仪、加速度计、惯性测量单元)来估计,而三维世界位置和平移速度及飞行方向等信息通常由卫星全球定位导航传感器(GPS)来估计。定义飞行轨迹状态包括北东地惯性系的当前无人机位置p的向量信息
Figure RE-GDA0002615882700000101
地速V,飞行速度方向与正北方向的夹角航线角χ。
无人机的任务规划系统为了圆满完成飞行任务,预先设置飞行途中必须经过的一系列关键目标飞行航路点A点、B点、C点和D点,期望按照航路点A-B-C-D的顺序飞行,这样对于每一个航点就包含了三维位置和期望飞行方向的信息,如飞行时需要满足航路点A到航路点B的直线期望飞行方向;然而,通常由于起飞场地的限制,如北东地坐标系中起飞点位置向量
Figure RE-GDA0002615882700000102
和跑道方向或无人机弹射起飞方向与第一个航路点A的位置和方向不一致;又由于无人机起飞条件受到现场风向影响较大,通常期望固定翼无人机初始逆风起飞来增加相对空气的速度来提高升力,这样起飞点位置向量
Figure RE-GDA0002615882700000103
的位置和起飞方向与第一个航路点A的位置和方向也不一致。针对固定翼无人机起飞段路径规划这一特殊飞行阶段,为增强无人机自主化程度,本发明提出的在线规划从起飞点切入第一个航路点的合理可行路径方案流程图如图2所示,具体分别在纵向平面和水平平面两个维度进行路径规划:
(1)纵向随实际水平待飞距离变化的爬升路径规划策略
如图3所示,设起飞点位置向量
Figure RE-GDA0002615882700000111
处的当前高度为H1,由于第一个期望航路点A处的期望高度为HA,为了让飞机有一个较为平稳的爬升过程,高度指令一般以直线斜坡信号的方式给出,即按照起飞方向纵向平面内选择一个与航路点A高度相同的中间航路点p2,其位置向量为
Figure RE-GDA0002615882700000112
其高度H2与A处的高度HA相同,则
Figure RE-GDA0002615882700000113
Figure RE-GDA0002615882700000114
的期望高度偏差ΔH12
ΔH12=H2-H1 (1)
定义
Figure RE-GDA0002615882700000115
Figure RE-GDA0002615882700000116
的位置向量为
Figure RE-GDA0002615882700000117
定义
Figure RE-GDA0002615882700000118
Figure RE-GDA0002615882700000119
的水平投影距离为|p1p2|,这也是实际水平待飞距离。爬升指令即根据飞机的当前位置,计算出当前的高度指令,通常给定飞机一个固定的爬升角,按照固定的爬升角飞行到达期望的高度。定义期望爬升角即期望航迹倾斜角γ
Figure RE-GDA00026158827000001110
Figure RE-GDA00026158827000001111
Figure RE-GDA00026158827000001112
的水平投影距离|p1p2|即实际水平待飞距离较短时,同时高度偏差ΔH12较大时,飞机的期望航迹倾斜角γ会变大,当该角过大时,俯仰角过大使飞行器产生失速,失速会导致机翼升力不足而坠机。因此这种直线爬升的方式仅适合γ较小时的情况,考虑到无人飞机的失速迎角一般在15 度左右,考虑预留一定的裕度,以期望航迹倾斜角10°为分界线(如果航迹倾斜角大于10°则飞行器产生失速),本发明纵向随实际水平待飞距离变化的爬升路径规划策略如下:
1)当
Figure RE-GDA0002615882700000121
时,即γ≤10°时采用直线爬升的方式;
直线爬升的高度指令Hg如下:
Figure RE-GDA0002615882700000122
式中p代表无人机当前的位置,定义飞行点到当前无人机的水平投影距离为|p1p|,变量tanγ应该满足tanγ<tan10°。
2)当
Figure RE-GDA0002615882700000123
时,即γ>10°时采用直线爬升和盘旋爬升结合的方式,直线爬升的期望航迹倾斜角为10°,先按期望航迹倾斜角10°直线爬升到中间航路点p2下方的水平投影p3点处,再按盘旋的半径Rpan盘旋爬升到中间航路点p2,盘旋的半径Rpan为预先设定,盘旋爬升时的飞行器滚转角φ满足飞行时最大滚转角φmax的约束,即φ≤φmax
(2)水平方向上规划出满足第一个航路点期望方向的转弯路径策略
当无人机到达中间航路点p2后,此时在高度上和第一个指令航路点A 的高度相同,但由于按照直线方向爬升,到p2点时的飞行方向一般和航路点A的期望方向不同,到达中间航点后需要进行水平转弯,本发明水平转弯路径规化策略如下:
采用基于Dubins路径的规划转弯策略,即从一个航点配置(位置和航向)过渡到另一个航点配置时,两个不同配置之间的时间最优路径由一段圆弧跟着一段直线,结尾再由一段圆弧到达最总配置,这种圆弧—直线—圆弧的路径对于固定高度和恒定速度的情况是时间最优路径。
假设无人机在执行水平飞行时飞控系统能够保持地速V绝对值不变,定义圆弧路径的半径用R表示,半径与无
人机的转弯半径一样大。对于任何给定的开始和结束配置,如图4所示,由一段圆弧接一段直线再加上一段圆弧的路径有四种可能。
第一种情况:右圆弧—直线—右圆弧;
第二种情况:右圆弧—直线—左圆弧;
第三种情况:左圆弧—直线—右圆弧;
第四种情况:左圆弧—直线—左圆弧;
为了确定路径最短的路径,需要计算四种情况的路径长度。给定无人机的位置向量信息
Figure RE-GDA0002615882700000131
航向χ和转弯半径R,右转和左转圆弧的中心位置向量
Figure RE-GDA0002615882700000132
Figure RE-GDA0002615882700000133
分别由下式给出:
Figure RE-GDA0002615882700000134
Figure RE-GDA0002615882700000135
式中符号T代表向量的转置。
为了计算不同情况轨迹的长度,需要分别针对顺时针圆和逆时针圆确定角距离,假设角度θ1和θ2分别在0和2π之间,对于顺时针圆,θ1和θ2之间的角距离如下:
21|=<2π+θ21> (6)
式中,符号<θ>代表角θ对2π取模,即
<θ>=θmod 2π (7)
同样,对于逆时针圆,θ1和θ2之间的角距离如下:
21|=<2π-θ21> (8)
假设Dubins路径开始圆心
Figure RE-GDA0002615882700000136
下标为s,结束圆心
Figure RE-GDA0002615882700000137
下标为e,开始圆心
Figure RE-GDA0002615882700000138
与结束圆心
Figure RE-GDA0002615882700000139
连线与正北方向的夹角定义为
Figure RE-GDA00026158827000001310
则开始圆心到结束圆心的距离l如下:
Figure RE-GDA00026158827000001311
因此,四种情况的路径长度计算公式分别如下:
第一种情况中右圆弧—直线—右圆弧的路径长度L1如下:
Figure RE-GDA0002615882700000141
第二种情况中右圆弧—直线—左圆弧的路径长度L2如下:
Figure RE-GDA0002615882700000142
其中,
Figure RE-GDA0002615882700000143
第三种情况中左圆弧—直线—右圆弧的路径长度L3如下:
Figure RE-GDA0002615882700000144
其中,
Figure RE-GDA0002615882700000145
第四种情况中左圆弧—直线—左圆弧的路径长度L4如下:
Figure RE-GDA0002615882700000146
当无人机到达中间航路点p2后,按照以上四种情况分别计算路径长度,由于仅需要位置信息,计算量很小,可以在线完成路径规划,最终选择一条在水平面路径长度最短,并且满足从当前飞行位置和飞行方向到第一个航路点期望方向的转弯路径,再由无人机制导控制系统来导引跟踪这条规划好的路径。
综上所述,将纵向路径规划与水平路径规划连接在一起,满足了在线规划无人机从起飞点切入第一个航路点的合理可行路径方案的要求。
本发明不仅局限于上述具体实施方式,本领域一般技术人员根据本发明公开的内容,可以采用其它多种具体实施方式实施本发明,因此,凡是采用本发明的设计结构和思路,做一些简单的变化或更改的设计,都落入本发明保护的范围。

Claims (8)

1.一种固定翼无人机起飞段切入航路点的在线路径规划方法,其特征在于,包括如下步骤:
S1、在起飞段,无人机通过导航系统器件和卫星接收机测量自身的飞行姿态、速度和地理位置三维坐标信息;
S2、将飞行器自身三维位置高度和飞行速度方向等信息与期望第一个航路点三维位置高度和期望飞行速度方向等信息进行比较,以便于规划爬升和转弯路径;
S3、纵向上规划出到第一个航路点高度的可行爬升路径;
S4、水平方向上规划出满足第一个航路点期望方向的转弯路径。
2.根据权利要求1所述的固定翼无人机起飞段切入航路点的在线路径规划,其特征在于,所述步骤S1中无人机测量自身的飞行姿态、速度和地理位置三维坐标信息的具体方法为:无人机相对地面惯性坐标系的姿态角信息由惯性传感器来估计,而地理位置三维坐标、水平飞行速度和速度方向信息由卫星全球定位导航传感器GPS来估计。
3.根据权利要求1所述的固定翼无人机起飞段切入航路点的在线路径规划,其特征在于,所述步骤S3中纵向上规划出到第一个航路点高度的可行爬升路径时,充分考虑无人机爬升角过大可能会导致失速的的工程实践风险,设计随实际水平待飞距离变化的直线爬升和盘旋爬升两种路径规划策略。
4.根据权利要求1所述的固定翼无人机起飞段切入航路点的在线路径规划,其特征在于,所述步骤S4中水平方向上规划出满足第一个航路点期望方向的转弯路径时,是一种飞行路径最短的优化路径。
5.根据权利要求4所述的的固定翼无人机起飞段切入航路点的在线路径规划,其特征在于,所述飞行路径最短的优化路径采用基于Dubins路径的规划转弯策略予以规划。
6.根据权利要求3所述的固定翼无人机起飞段切入航路点的在线路径规划,其特征在于,所述纵向随实际水平待飞距离变化的爬升路径规划策略具体为:
设起飞点位置向量
Figure RE-FDA0002615882690000011
处的当前高度为H1,由于第一个期望航路点A处的期望高度为HA,为了让飞机有一个较为平稳的爬升过程,高度指令一般以直线斜坡信号的方式给出,即按照起飞方向纵向平面内选择一个与航路点A高度相同的中间航路点p2,其位置向量为
Figure RE-FDA0002615882690000021
其高度H2与A处的高度HA相同,则
Figure RE-FDA0002615882690000022
Figure RE-FDA0002615882690000023
的期望高度偏差ΔH12
ΔH12=H2-H1 (1)
定义
Figure RE-FDA0002615882690000024
Figure RE-FDA0002615882690000025
的位置向量为
Figure RE-FDA0002615882690000026
定义
Figure RE-FDA0002615882690000027
Figure RE-FDA0002615882690000028
的水平投影距离为|p1p2|,这也是实际水平待飞距离。爬升指令即根据飞机的当前位置,计算出当前的高度指令,通常给定飞机一个固定的爬升角,按照固定的爬升角飞行到达期望的高度。定义期望爬升角即期望航迹倾斜角γ
Figure RE-FDA0002615882690000029
Figure RE-FDA00026158826900000210
Figure RE-FDA00026158826900000211
的水平投影距离|p1p2|即实际水平待飞距离较短时,同时高度偏差ΔH12较大时,飞机的期望航迹倾斜角γ会变大,当该角过大时,俯仰角过大使飞行器产生失速,失速会导致机翼升力不足而坠机。因此这种直线爬升的方式仅适合γ较小时的情况,考虑到无人飞机的失速迎角一般在15度左右,考虑预留一定的裕度,以期望航迹倾斜角10°为分界线:
1)当
Figure RE-FDA00026158826900000212
时,即γ≤10°时采用直线爬升的方式;
直线爬升的高度指令Hg如下:
Figure RE-FDA00026158826900000213
式中p代表无人机当前的位置,定义飞行点到当前无人机的水平投影距离为|p1p|,变量tanγ应该满足tanγ<tan10°。
2)当
Figure RE-FDA00026158826900000214
时,即γ>10°时采用直线爬升和盘旋爬升结合的方式,直线爬升的期望航迹倾斜角为10°,先按期望航迹倾斜角10°直线爬升到中间航路点p2下方的水平投影p3点处,再按盘旋的半径Rpan盘旋爬升到中间航路点p2,盘旋的半径Rpan为预先设定,盘旋爬升时的飞行器滚转角φ满足飞行时最大滚转角φmax的约束,即φ≤φmax
7.根据权利要求5所述的固定翼无人机起飞段切入航路点的在线路径规划,其特征在于,所述飞行路径最短的优化路径采用基于Dubins路径的规划转弯策略予以规划的策略具体为:
从一个航点配置(位置和航向)过渡到另一个航点配置时,两个不同配置之间的时间最优路径由一段圆弧跟着一段直线,结尾再由一段圆弧到达最总配置,这种圆弧—直线—圆弧的路径对于固定高度和恒定速度的情况是时间最优路径;
假设无人机在执行水平飞行时飞控系统能够保持地速V绝对值不变,定义圆弧路径的半径用R表示,半径与无人机的转弯半径一样大;对于任何给定的开始和结束配置,由一段圆弧接一段直线再加上一段圆弧的路径有四种情况,
第一种情况:右圆弧—直线—右圆弧;
第二种情况:右圆弧—直线—左圆弧;
第三种情况:左圆弧—直线—右圆弧;
第四种情况:左圆弧—直线—左圆弧;
为了确定路径最短的路径,需要计算四种情况的路径长度;给定无人机的位置向量信息
Figure RE-FDA0002615882690000031
航向χ和转弯半径R,右转和左转圆弧的中心位置向量
Figure RE-FDA0002615882690000032
Figure RE-FDA0002615882690000033
分别由下式给出:
Figure RE-FDA0002615882690000034
Figure RE-FDA0002615882690000035
式中符号T代表向量的转置;
为了计算不同情况轨迹的长度,需要分别针对顺时针圆和逆时针圆确定角距离,假设角度θ1和θ2分别在0和2π之间,对于顺时针圆,θ1和θ2之间的角距离如下:
21|=<2π+θ21> (6)
式中,符号<θ>代表角θ对2π取模,即
<θ>=θmod2π (7)
同样,对于逆时针圆,θ1和θ2之间的角距离如下:
21|=<2π-θ21> (8)
假设Dubins路径开始圆心
Figure RE-FDA0002615882690000041
下标为s,结束圆心
Figure RE-FDA0002615882690000042
下标为e,开始圆心
Figure RE-FDA0002615882690000043
与结束圆心
Figure RE-FDA0002615882690000044
连线与正北方向的夹角定义为
Figure RE-FDA0002615882690000045
则开始圆心到结束圆心的距离l如下:
Figure RE-FDA0002615882690000046
因此,四种情况的路径长度计算公式分别如下:
第一种情况中右圆弧—直线—右圆弧的路径长度L1如下:
Figure RE-FDA0002615882690000047
第二种情况中右圆弧—直线—左圆弧的路径长度L2如下:
Figure RE-FDA0002615882690000048
其中,
Figure RE-FDA0002615882690000049
第三种情况中左圆弧—直线—右圆弧的路径长度L3如下:
Figure RE-FDA0002615882690000051
其中,
Figure RE-FDA0002615882690000052
第四种情况中左圆弧—直线—左圆弧的路径长度L4如下:
Figure RE-FDA0002615882690000053
8.根据权利要求6所述的固定翼无人机起飞段切入航路点的在线路径规划,其特征在于,所述无人机从起飞点位置P1纵向随实际水平待飞距离变化的爬升路径规划策略爬升到中间航路点p2后,按照所述飞行路径最短的优化路径采用基于Dubins路径的规划转弯策略,规划其到达第一个航路点A期望方向的转弯路径的策略具体为:
从一个航点配置包括位置和航向,过渡到另一个航点配置时,两个不同配置之间的时间最优路径由一段圆弧跟着一段直线,结尾再由一段圆弧到达最终配置,这种圆弧—直线—圆弧的路径对于固定高度和恒定速度的情况是时间最优路径;
假设无人机在执行水平飞行时飞控系统能够保持地速V绝对值不变,定义圆弧路径的半径用R表示,半径与无人机的转弯半径一样大;对于任何给定的开始和结束配置,由一段圆弧接一段直线再加上一段圆弧的路径有四种情况,
第一种情况:右圆弧—直线—右圆弧;
第二种情况:右圆弧—直线—左圆弧;
第三种情况:左圆弧—直线—右圆弧;
第四种情况:左圆弧—直线—左圆弧;
为了确定路径最短的路径,需要计算四种情况的路径长度;给定无人机的位置向量信息
Figure RE-FDA0002615882690000061
航向χ和转弯半径R,右转和左转圆弧的中心位置向量
Figure RE-FDA0002615882690000062
Figure RE-FDA0002615882690000063
分别由下式给出:
Figure RE-FDA0002615882690000064
Figure RE-FDA0002615882690000065
式中符号T代表向量的转置;
为了计算不同情况轨迹的长度,需要分别针对顺时针圆和逆时针圆确定角距离,假设角度θ1和θ2分别在0和2π之间,对于顺时针圆,θ1和θ2之间的角距离如下:
21|=<2π+θ21> (6)
式中,符号<θ>代表角θ对2π取模,即
<θ>=θmod2π (7)
同样,对于逆时针圆,θ1和θ2之间的角距离如下:
21|=<2π-θ21> (8)
假设Dubins路径开始圆心
Figure RE-FDA0002615882690000066
下标为s,结束圆心
Figure RE-FDA0002615882690000067
下标为e,开始圆心
Figure RE-FDA0002615882690000068
与结束圆心
Figure RE-FDA0002615882690000069
连线与正北方向的夹角定义为
Figure RE-FDA00026158826900000610
则开始圆心到结束圆心的距离l如下:
Figure RE-FDA00026158826900000611
因此,四种情况的路径长度计算公式分别如下:
第一种情况中右圆弧—直线—右圆弧的路径长度L1如下:
Figure RE-FDA0002615882690000071
第二种情况中右圆弧—直线—左圆弧的路径长度L2如下:
Figure RE-FDA0002615882690000072
其中,
Figure RE-FDA0002615882690000073
第三种情况中左圆弧—直线—右圆弧的路径长度L3如下:
Figure RE-FDA0002615882690000074
其中,
Figure RE-FDA0002615882690000075
第四种情况中左圆弧—直线—左圆弧的路径长度L4如下:
Figure RE-FDA0002615882690000076
当无人机到达中间航路点p2后,按照以上四种情况分别计算路径长度,由于仅需要位置信息,计算量很小,可以在线完成路径规划,最终选择一条在水平面路径长度最短,并且满足从当前飞行位置和飞行方向到第一个航路点A期望方向的转弯路径,再由无人机制导控制系统来导引跟踪这条规划好的路径。
CN201911287849.0A 2019-12-15 2019-12-15 一种固定翼无人机起飞段切入航路点的在线路径规划方法 Active CN111650958B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911287849.0A CN111650958B (zh) 2019-12-15 2019-12-15 一种固定翼无人机起飞段切入航路点的在线路径规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911287849.0A CN111650958B (zh) 2019-12-15 2019-12-15 一种固定翼无人机起飞段切入航路点的在线路径规划方法

Publications (2)

Publication Number Publication Date
CN111650958A true CN111650958A (zh) 2020-09-11
CN111650958B CN111650958B (zh) 2023-11-03

Family

ID=72342531

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911287849.0A Active CN111650958B (zh) 2019-12-15 2019-12-15 一种固定翼无人机起飞段切入航路点的在线路径规划方法

Country Status (1)

Country Link
CN (1) CN111650958B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112306077A (zh) * 2020-11-06 2021-02-02 广州极飞科技有限公司 飞行控制方法、装置、飞行器及存储介质
CN112306091A (zh) * 2020-10-28 2021-02-02 深圳市领峰电动智能科技有限公司 一种三维轨迹规划方法、装置、设备及存储介质
CN112684810A (zh) * 2020-12-18 2021-04-20 彩虹无人机科技有限公司 一种固定翼无人机环绕飞行的导航方法
CN113625741A (zh) * 2021-09-01 2021-11-09 北京北航天宇长鹰无人机科技有限公司 飞行段控制方法及装置
CN114111797A (zh) * 2021-11-30 2022-03-01 北京信息科技大学 基于fpga的卡尔曼滤波器、ip核及导航用芯片
CN114111773A (zh) * 2021-11-30 2022-03-01 北京信息科技大学 组合导航方法、装置、系统及存储介质
CN115097863A (zh) * 2022-06-24 2022-09-23 西北工业大学 一种基于Dubins方法的水平引导过渡路径构建方法
CN116736747A (zh) * 2023-08-15 2023-09-12 陕西德鑫智能科技有限公司 一种无人机应急处理方法及装置
CN117275292A (zh) * 2023-11-08 2023-12-22 北京中兵数字科技集团有限公司 用于航空器单发离场的航径规划方法、装置和计算设备
CN117492474A (zh) * 2022-07-22 2024-02-02 海鹰航空通用装备有限责任公司 无人机航迹自主导航获取方法及无人机侦察抵近控制方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090112535A1 (en) * 2007-10-26 2009-04-30 Sensis Corporation Method of integrating point mass equations to include vertical and horizontal profiles
EP2921823A1 (en) * 2014-03-18 2015-09-23 BAE Systems PLC Path planning
US20160364988A1 (en) * 2014-03-18 2016-12-15 Bae Systems Plc Path planning
CN108332753A (zh) * 2018-01-30 2018-07-27 北京航空航天大学 一种无人机电力巡检路径规划方法
CN108388263A (zh) * 2018-03-01 2018-08-10 中国航空无线电电子研究所 圆形航线自动飞行系统
CN109708639A (zh) * 2018-12-07 2019-05-03 湖北航天飞行器研究所 飞行器平飞跟踪直线和圆弧路径的侧向制导指令生成方法
CN109947123A (zh) * 2019-02-27 2019-06-28 南京航空航天大学 一种基于视线导引律的无人机路径跟踪与自主避障方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090112535A1 (en) * 2007-10-26 2009-04-30 Sensis Corporation Method of integrating point mass equations to include vertical and horizontal profiles
EP2921823A1 (en) * 2014-03-18 2015-09-23 BAE Systems PLC Path planning
US20160364988A1 (en) * 2014-03-18 2016-12-15 Bae Systems Plc Path planning
CN108332753A (zh) * 2018-01-30 2018-07-27 北京航空航天大学 一种无人机电力巡检路径规划方法
CN108388263A (zh) * 2018-03-01 2018-08-10 中国航空无线电电子研究所 圆形航线自动飞行系统
CN109708639A (zh) * 2018-12-07 2019-05-03 湖北航天飞行器研究所 飞行器平飞跟踪直线和圆弧路径的侧向制导指令生成方法
CN109947123A (zh) * 2019-02-27 2019-06-28 南京航空航天大学 一种基于视线导引律的无人机路径跟踪与自主避障方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
赵锋;杨伟;杨朝旭;孙绍山;: "无人机三维航路动态规划及导引控制研究" *
项梅;段刘刚;: "一种无人机转弯航路平滑算法及实现" *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022089506A1 (zh) * 2020-10-28 2022-05-05 深圳市领峰电动智能科技有限公司 一种三维轨迹规划方法、装置、设备及存储介质
CN112306091A (zh) * 2020-10-28 2021-02-02 深圳市领峰电动智能科技有限公司 一种三维轨迹规划方法、装置、设备及存储介质
CN112306091B (zh) * 2020-10-28 2023-06-20 峰飞航空科技(昆山)有限公司 一种三维轨迹规划方法、装置、设备及存储介质
CN112306077A (zh) * 2020-11-06 2021-02-02 广州极飞科技有限公司 飞行控制方法、装置、飞行器及存储介质
CN112684810A (zh) * 2020-12-18 2021-04-20 彩虹无人机科技有限公司 一种固定翼无人机环绕飞行的导航方法
CN113625741A (zh) * 2021-09-01 2021-11-09 北京北航天宇长鹰无人机科技有限公司 飞行段控制方法及装置
CN113625741B (zh) * 2021-09-01 2022-03-08 北京北航天宇长鹰无人机科技有限公司 飞行段控制方法及装置
CN114111773B (zh) * 2021-11-30 2023-11-17 北京信息科技大学 组合导航方法、装置、系统及存储介质
CN114111773A (zh) * 2021-11-30 2022-03-01 北京信息科技大学 组合导航方法、装置、系统及存储介质
CN114111797A (zh) * 2021-11-30 2022-03-01 北京信息科技大学 基于fpga的卡尔曼滤波器、ip核及导航用芯片
CN114111797B (zh) * 2021-11-30 2024-02-20 北京信息科技大学 基于fpga的卡尔曼滤波器、ip核及导航用芯片
CN115097863A (zh) * 2022-06-24 2022-09-23 西北工业大学 一种基于Dubins方法的水平引导过渡路径构建方法
CN117492474A (zh) * 2022-07-22 2024-02-02 海鹰航空通用装备有限责任公司 无人机航迹自主导航获取方法及无人机侦察抵近控制方法
CN116736747A (zh) * 2023-08-15 2023-09-12 陕西德鑫智能科技有限公司 一种无人机应急处理方法及装置
CN116736747B (zh) * 2023-08-15 2023-11-14 陕西德鑫智能科技有限公司 一种无人机应急处理方法及装置
CN117275292A (zh) * 2023-11-08 2023-12-22 北京中兵数字科技集团有限公司 用于航空器单发离场的航径规划方法、装置和计算设备
CN117275292B (zh) * 2023-11-08 2024-04-09 北京中兵数字科技集团有限公司 用于航空器单发离场的航径规划方法、装置和计算设备

Also Published As

Publication number Publication date
CN111650958B (zh) 2023-11-03

Similar Documents

Publication Publication Date Title
CN111650958B (zh) 一种固定翼无人机起飞段切入航路点的在线路径规划方法
US9274529B2 (en) Safe emergency landing of a UAV
US9257048B1 (en) Aircraft emergency landing route system
US20180102058A1 (en) High-precision autonomous obstacle-avoidance flying method for unmanned aerial vehicle
US9520066B2 (en) Determining landing sites for aircraft
CN102582826B (zh) 一种四旋翼无人飞行器的驾驶方法和系统
US7689328B2 (en) Determining suitable areas for off-airport landings
US6963795B2 (en) Vehicle position keeping system
CN101095090B (zh) 用于自动绕圈飞行的控制系统
US20130035810A1 (en) Method and system to autonomously direct aircraft to emergency-contingency landing sites using on-board sensors
CN109708639B (zh) 飞行器平飞跟踪直线和圆弧路径的侧向制导指令生成方法
CN102620736A (zh) 一种无人机的导航方法
US10502584B1 (en) Mission monitor and controller for autonomous unmanned vehicles
CN110764527A (zh) 一种无人机无动力应急返场在线航迹规划方法
CN112748743A (zh) 空中运载器导航系统
US11535394B2 (en) Aircraft landing assistance method and memory storage device including instructions for performing an aircraft landing assistance method
Harbick et al. Planar spline trajectory following for an autonomous helicopter
JP5166349B2 (ja) 固定翼機、固定翼機システムおよび固定翼機の着陸方法
Jantawong et al. Automatic landing control based on GPS for fixed-wing aircraft
Sedlmair et al. Design and experimental validation of UAV control laws-3D spline-path-following and easy-handling remote control
CN109213186B (zh) 无人机喷洒作业的方法、装置以及无人机
Eng et al. Simulation of a fixed-wing UAV forced landing with dynamic path planning
EP4181104A1 (en) Systems and methods for providing safe landing assistance for a vehicle
CN114724413B (zh) 一种固定翼飞行器水平航线切换的引导方法
US20230076554A1 (en) Automatic aircraft taxiing

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