CN105843232A - 一种飞行器滑翔减速控制方法 - Google Patents
一种飞行器滑翔减速控制方法 Download PDFInfo
- Publication number
- CN105843232A CN105843232A CN201610217768.3A CN201610217768A CN105843232A CN 105843232 A CN105843232 A CN 105843232A CN 201610217768 A CN201610217768 A CN 201610217768A CN 105843232 A CN105843232 A CN 105843232A
- Authority
- CN
- China
- Prior art keywords
- aircraft
- calculated
- ctrl
- flag
- control
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000033001 locomotion Effects 0.000 claims abstract description 6
- 238000005553 drilling Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010006 flight Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/04—Control of altitude or depth
- G05D1/06—Rate of change of altitude or depth
- G05D1/0607—Rate of change of altitude or depth specially adapted for aircraft
- G05D1/0653—Rate of change of altitude or depth specially adapted for aircraft during a phase of take-off or landing
- G05D1/0676—Rate of change of altitude or depth specially adapted for aircraft during a phase of take-off or landing specially adapted for landing
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Regulating Braking Force (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种飞行器滑翔减速控制方法,该方法包括:根据制导系统给出的指令速度Vcx与导航系统获取的飞行器当前相对地球运动速度Vd的差值,判断是否需要进行减速控制;当需要进行减速控制时,计算得到基本需用攻角α0;计算得到需要耗散掉的速度ΔV;计算得到减速需用攻角αn;计算得到减速控制附加的制导力根据计算得到的减速控制附加的制导力,对飞行器进行减速控制。通过使用本发明所提供的方法,可以实现对飞行器的精确的速度控制。
Description
技术领域
本发明涉及高超声速飞行器制导技术领域,特别涉及一种飞行器滑翔减速控制方法。
背景技术
对于升力式飞行器,由于其主要依赖气动力飞行,在不确定因素和干扰条件的影响下,飞行结束时飞行器的速度散布大,对后续飞行造成不利,为约束终端点的速度散布,需要对飞行器进行严格的速度控制,使其在干扰条件下速度趋于一致。
然而,现有技术中还尚未能提供一种有效且精度较高的飞行器滑翔减速控制方法,因此该问题亟待得到解决。
发明内容
有鉴于此,本发明提供一种飞行器滑翔减速控制方法,从而可以实现对飞行器的精确的速度控制。
本发明的技术方案具体是这样实现的:
一种飞行器滑翔减速控制方法,该方法包括如下步骤:
根据制导系统给出的指令速度Vcx与导航系统获取的飞行器当前相对地球运动速度Vd的差值,判断是否需要进行减速控制;
当需要进行减速控制时,计算得到基本需用攻角α0;
计算得到需要耗散掉的速度ΔV;
计算得到减速需用攻角αn;
计算得到减速控制附加的制导力
根据计算得到的减速控制附加的制导力,对飞行器进行减速控制。
较佳的,所述根据Vcx与Vd的差值,判断是否需要进行减速控制包括:
预设一个减速控制标志VCtrl_Flag,且该VCtrl_Flag的初始值为0;
当(Vd-Vcx)>Vc1,且VCtrl_Flag的值为0时,置VCtrl_Flag为1;当(Vd-Vcx)<Vc2,且VCtrl_Flag的值为1时,置VCtrl_Flag为0;其中,Vc1为进入减速控制的门槛值,Vc2为退出减速控制的门槛值;
当VCtrl_Flag为0时,判断不需要进行减速控制;当VCtrl_Flag为1时,判断需要进行减速控制。
较佳的,根据如下的公式计算得到基本需用攻角α0:
其中,ax、bx分别为与飞行器法向力气动特性相关的参数,q为动压,Sm为飞行器特征面积,Fn为无减速控制时的制导力。
较佳的,根据如下的公式计算得到需要耗散掉的速度ΔV:
其中,Tg为距离飞行结束的待飞时间。
较佳的,根据如下的公式计算得到减速需用攻角αn:
其中,bD为与飞行器阻力系数相关的参数;KV1、N为制导参数;m为飞行器质量。
较佳的,所述计算得到减速控制附加的制导力包括:
计算得到Fn0=(ax+bx·α0)·q·Sm;
计算得到Fnn=(ax+bx·αn)·q·Sm;
如果Fnn>Fn0,则否则,
如上可见,在本发明中的飞行器滑翔减速控制方法中,由于可以先根据实际速度(即飞行器当前相对地球运动速度)与指令速度的偏差,判断是否需要进行减速控制,然后在需要进行减速控制时,先计算基本需用攻角,然后分别计算需要耗散掉的速度和减速需用攻角,最终计算得到减速控制附加的制导力,并根据减速控制附加的制导力对飞行器进行减速控制,从而可以增大飞行器的飞行攻角,通过增加攻角继而增加阻力,以实现对飞行器的精确的速度控制。上述方法可直接用于滑翔飞行器(例如,以气动力飞行为主的升力式飞行器)的制导,速度控制有效,且速度控制精度高。
附图说明
图1为本发明实施例中的飞行器滑翔减速控制方法的流程示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明进一步详细说明。
本实施例提供了一种飞行器滑翔减速控制方法,该方法适用于以气动力飞行为主的升力式飞行器,从而可以解决飞行器在平衡滑翔飞行状态下的高精度速度控制问题。
图1为本发明实施例中的飞行器滑翔减速控制方法的流程示意图。如图1所示,本发明实施例中的飞行器滑翔减速控制方法包括:。
步骤101,根据制导系统给出的指令速度Vcx与导航系统获取的飞行器当前相对地球运动速度Vd的差值,判断是否需要进行减速控制;
在本发明的技术方案中,首先要判断是否需要进行减速控制。因此,在本步骤中,可以根据Vcx与Vd的差值来判断是否需要进行减速控制,其中,Vcx为制导系统给出的指令速度,Vd为导航系统获取的飞行器当前相对地球运动速度Vd,单位均为m/s。
较佳的,在本发明的具体实施例中,可以根据如下所述的方法来判断是否需要进行减速控制:
预设一个减速控制标志VCtrl_Flag,且该VCtrl_Flag的初始值为0;该VCtrl_Flag用于标志是否需要进行减速控制;
当(Vd-Vcx)>Vc1,且VCtrl_Flag的值为0时,置VCtrl_Flag为1;当(Vd-Vcx)<Vc2,且VCtrl_Flag的值为1时,置VCtrl_Flag为0;其中,Vc1为进入减速控制的门槛值,Vc2为退出减速控制的门槛值;在本发明的技术方案中,可以根据实际应用的需要,预先设置上述Vc1和Vc2的取值;因此,一般情况下,上述Vc1和Vc2均为预设的常值。
当VCtrl_Flag为0时,判断不需要进行减速控制,设减速所附加的侧向力即减速控制所附加的制导力从而退出减速控制流程;
当VCtrl_Flag为1时,判断需要进行减速控制,从而继续执行后续的步骤102。
步骤102,当需要进行减速控制时,计算得到基本需用攻角α0;
对于以气动力飞行为主的升力式飞行器,其飞行所需的力均由飞行器的法向力来提供,而减速则主要是通过增加飞行器的攻角来实现,因此,在本发明的技术方案中,如果通过上述的步骤101确定需要进行减速控制,则可在本步骤中,首先计算基本需用攻角α0。
在本发明的技术方案中,可以根据多种具体的实现方式来计算得到基本需用攻角α0。以下将以其中的一种具体实现方式为例,对本发明的技术方案进行说明。
例如,较佳的,在本发明的一个具体实施例中,可以根据动压、飞行器特征面积和无减速控制时的制导力,计算基本需用攻角α0。例如,可以根据如下所述的公式计算得到基本需用攻角α0:
其中,ax、bx分别为与飞行器法向力气动特性相关的第一和第二参数,ax、bx可以根据飞行器的气动特性拟合得到;q为动压,Sm为飞行器特征面积,Fn为无减速控制时的制导力。
在本发明的技术方案中,所述无减速控制时的制导力Fn可以根据现有技术中的常用制导算法计算得到,因此,具体计算过程在此不再赘述。
步骤103,计算得到需要耗散掉的速度ΔV;
在本发明的技术方案中,可以根据多种具体的实现方式来计算得到需要耗散掉的速度ΔV。以下将以其中的一种具体实现方式为例,对本发明的技术方案进行说明。
例如,较佳的,在本发明的一个具体实施例中,可以根据Vcx、Vd和距离飞行结束的待飞时间,计算需要耗散掉的速度ΔV。例如,可以根据如下所述的公式计算得到需要耗散掉的速度ΔV:
其中,Tg为距离飞行结束的待飞时间。在本发明的技术方案中,可以根据实际情况,预先设置Tg的取值,即设Tg为一个常值;或者,也可以根据现有技术中的常用制导算法计算得到该Tg,具体计算过程在此不再赘述。
步骤104,计算得到减速需用攻角αn;
在本发明的技术方案中,可以根据多种具体的实现方式来计算得到减速需用攻角αn。以下将以其中的一种具体实现方式为例,对本发明的技术方案进行说明。
例如,较佳的,在本发明的一个具体实施例中,可以根据如下所述的公式计算得到减速需用攻角αn:
其中,bD为与飞行器阻力系数相关的参数,其与飞行气动特性相关;KV1为速度比例系数,其值一般可以取为1,、N为由飞行器气动特性拟合得到的攻角函数的阶次;m为飞行器质量。
步骤105,计算得到减速控制附加的制导力
在本发明的技术方案中,可以根据多种具体的实现方式来计算得到减速控制附加的制导力以下将以其中的一种具体实现方式为例,对本发明的技术方案进行说明。
例如,较佳的,在本发明的一个具体实施例中,可以根据如下所述的步骤计算得到减速控制附加的制导力
步骤51,计算得到Fn0=(ax+bx·α0)·q·Sm,其中,Fn0为基本需用攻角产生的制导力;
步骤52,计算得到Fnn=(ax+bx·αn)·q·Sm,其中Fm为基本需用攻角产生的制导力;
步骤53,如果Fnn>Fn0,则否则,
步骤106,根据计算得到的减速控制附加的制导力,对飞行器进行减速控制。
综上可知,在本发明中的飞行器滑翔减速控制方法中,由于可以先根据实际速度(即飞行器当前相对地球运动速度)与指令速度的偏差,判断是否需要进行减速控制,然后在需要进行减速控制时,先计算基本需用攻角,然后分别计算需要耗散掉的速度和减速需用攻角,最终计算得到减速控制附加的制导力,并根据减速控制附加的制导力对飞行器进行减速控制,从而可以增大飞行器的飞行攻角,通过增加攻角继而增加阻力,以实现对飞行器的精确的速度控制。上述方法可直接用于滑翔飞行器(例如,以气动力飞行为主的升力式飞行器)的制导,速度控制有效,且速度控制精度高。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (6)
1.一种飞行器滑翔减速控制方法,其特征在于,该方法包括如下步骤:
根据制导系统给出的指令速度Vcx与导航系统获取的飞行器当前相对地球运动速度Vd的差值,判断是否需要进行减速控制;
当需要进行减速控制时,计算得到基本需用攻角α0;
计算得到需要耗散掉的速度ΔV;
计算得到减速需用攻角αn;
计算得到减速控制附加的制导力
根据计算得到的减速控制附加的制导力,对飞行器进行减速控制。
2.根据权利要求1所述的方法,其特征在于,所述根据Vcx与Vd的差值,判断是否需要进行减速控制包括:
预设一个减速控制标志VCtrl_Flag,且该VCtrl_Flag的初始值为0;
当(Vd-Vcx)>Vc1,且VCtrl_Flag的值为0时,置VCtrl_Flag为1;当(Vd-Vcx)<Vc2,且VCtrl_Flag的值为1时,置VCtrl_Flag为0;其中,Vc1为进入减速控制的门槛值,Vc2为退出减速控制的门槛值;
当VCtrl_Flag为0时,判断不需要进行减速控制;当VCtrl_Flag为1时,判断需要进行减速控制。
3.根据权利要求2所述的方法,其特征在于,根据如下的公式计算得到基本需用攻角α0:
其中,ax、bx分别为与飞行器法向力气动特性相关的参数,q为动压,Sm为飞行器特征面积,Fn为无减速控制时的制导力。
4.根据权利要求3所述的方法,其特征在于,根据如下的公式计算得到需要耗散掉的速度ΔV:
其中,Tg为距离飞行结束的待飞时间。
5.根据权利要求4所述的方法,其特征在于,根据如下的公式计算得到减速需用攻角αn:
其中,bD为与飞行器阻力系数相关的参数;KV1、N为制导参数;m为飞行器质量。
6.根据权利要求5所述的方法,其特征在于,所述计算得到减速控制附加的制导力包括:
计算得到Fn0=(ax+bx·α0)·q·Sm;
计算得到Fnn=(ax+bx·αn)·q·Sm;
如果Fnn>Fn0,则否则,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610217768.3A CN105843232B (zh) | 2016-04-08 | 2016-04-08 | 一种飞行器滑翔减速控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610217768.3A CN105843232B (zh) | 2016-04-08 | 2016-04-08 | 一种飞行器滑翔减速控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105843232A true CN105843232A (zh) | 2016-08-10 |
CN105843232B CN105843232B (zh) | 2017-03-22 |
Family
ID=56597154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610217768.3A Expired - Fee Related CN105843232B (zh) | 2016-04-08 | 2016-04-08 | 一种飞行器滑翔减速控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105843232B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109425754A (zh) * | 2017-08-22 | 2019-03-05 | 成都飞机工业(集团)有限责任公司 | 一种判断无人机减速异常的方法 |
CN114489125A (zh) * | 2022-01-15 | 2022-05-13 | 西北工业大学 | 一种滑翔飞行器高精度临近最优减速控制方法 |
US11358733B2 (en) | 2020-01-21 | 2022-06-14 | Rockwell Collins, Inc. | Optimum idle descent start point display based on current conditions |
CN114812293A (zh) * | 2021-01-27 | 2022-07-29 | 北京理工大学 | 一种末端减速机动控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102139768A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的攻角制导方法 |
US8489258B2 (en) * | 2009-03-27 | 2013-07-16 | The Charles Stark Draper Laboratory, Inc. | Propulsive guidance for atmospheric skip entry trajectories |
CN103587723A (zh) * | 2013-11-07 | 2014-02-19 | 北京临近空间飞行器系统工程研究所 | 一种再入初始段解析式纵向在线轨迹设计及跟踪方法 |
CN104176268A (zh) * | 2014-08-06 | 2014-12-03 | 北京航空航天大学 | 一种滑翔飞行弹道阻尼控制方法 |
CN104731104A (zh) * | 2015-03-09 | 2015-06-24 | 北京航天自动控制研究所 | 一种高超声速飞行器滑翔飞行段的纵向制导方法 |
CN104808492A (zh) * | 2015-03-23 | 2015-07-29 | 北京航天自动控制研究所 | 一种升力式飞行器的攻角指令生成方法 |
-
2016
- 2016-04-08 CN CN201610217768.3A patent/CN105843232B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8489258B2 (en) * | 2009-03-27 | 2013-07-16 | The Charles Stark Draper Laboratory, Inc. | Propulsive guidance for atmospheric skip entry trajectories |
CN102139768A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的攻角制导方法 |
CN103587723A (zh) * | 2013-11-07 | 2014-02-19 | 北京临近空间飞行器系统工程研究所 | 一种再入初始段解析式纵向在线轨迹设计及跟踪方法 |
CN104176268A (zh) * | 2014-08-06 | 2014-12-03 | 北京航空航天大学 | 一种滑翔飞行弹道阻尼控制方法 |
CN104731104A (zh) * | 2015-03-09 | 2015-06-24 | 北京航天自动控制研究所 | 一种高超声速飞行器滑翔飞行段的纵向制导方法 |
CN104808492A (zh) * | 2015-03-23 | 2015-07-29 | 北京航天自动控制研究所 | 一种升力式飞行器的攻角指令生成方法 |
Non-Patent Citations (2)
Title |
---|
刘礅等: "《空间飞行器动力学》", 31 October 2003, 哈尔滨工业大学出版社 * |
马宏图: "高可靠性再入轨迹在线设计及实时制导", 《CNKI硕士学位论文》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109425754A (zh) * | 2017-08-22 | 2019-03-05 | 成都飞机工业(集团)有限责任公司 | 一种判断无人机减速异常的方法 |
US11358733B2 (en) | 2020-01-21 | 2022-06-14 | Rockwell Collins, Inc. | Optimum idle descent start point display based on current conditions |
CN114812293A (zh) * | 2021-01-27 | 2022-07-29 | 北京理工大学 | 一种末端减速机动控制方法 |
CN114812293B (zh) * | 2021-01-27 | 2023-03-24 | 北京理工大学 | 一种末端减速机动控制方法 |
CN114489125A (zh) * | 2022-01-15 | 2022-05-13 | 西北工业大学 | 一种滑翔飞行器高精度临近最优减速控制方法 |
CN114489125B (zh) * | 2022-01-15 | 2024-06-11 | 西北工业大学 | 一种滑翔飞行器高精度临近最优减速控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105843232B (zh) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107941087B (zh) | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 | |
JP4328660B2 (ja) | 航空機の自動離陸装置、自動着陸装置及び自動離着陸装置並びに航空機の自動離陸方法、自動着陸方法及び自動離着陸方法 | |
CN105843232A (zh) | 一种飞行器滑翔减速控制方法 | |
KR20150043392A (ko) | 바람 벡터 추산 | |
US9856032B2 (en) | Method and device for controlling at least one actuator control system of an aircraft, associated computer program product and aircraft | |
CN107085978A (zh) | 一种基于所需到达时间的管制辅助决策指令生成方法 | |
US10261518B2 (en) | Method and apparatus for protecting aircraft maximum lift capability | |
CN103116359B (zh) | 一种无人机着陆下滑初始段的引导方法 | |
JP2015512818A5 (zh) | ||
CN103245257A (zh) | 基于Bezier曲线的多约束飞行器导引方法 | |
CN104019701B (zh) | 一种利用直接力气动力复合控制的前向拦截制导方法 | |
KR20190075770A (ko) | 저속비행상태에서 세로 자세 제어 신호에 기초하여 메인로터의 틸트 각도를 제어하는 방법 및 컴퓨터 프로그램과 수직 이착륙 비행체 | |
US9646505B2 (en) | Method of automatically controlling the descent phase of an aircraft using aircraft avionics executing a descent algorithm | |
US9684309B2 (en) | Aircraft guidance based on partial differential equation for miss distance | |
CN105573340A (zh) | 一种固定翼无人机抗侧风的飞行控制方法 | |
CN106248082A (zh) | 一种飞行器自主导航系统及导航方法 | |
CN104656657A (zh) | 一种常值风干扰平流层飞艇定点控制方法 | |
CN103885453A (zh) | 大角度进场性能改善及优化 | |
CN106372307A (zh) | 一种基于气动模型的民机气流角估计方法 | |
CN105759830B (zh) | 一种升力式飞行器高动态下压段制导方法 | |
WO2018107733A1 (zh) | 一种用于飞艇的控制方法和装置 | |
KR20160036369A (ko) | 동체만의 비행체 형상에 대한 측추력 상호 간섭 현상의 공력 모델링 기법 | |
CN104787361B (zh) | 升力式飞行器再入制导的阻力加速度变化率的确定方法 | |
US9889926B2 (en) | Air vehicles and systems for preemptive turbulence mitigation | |
CN114167885B (zh) | 一种升力式飞行器多模式解析制导方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170322 Termination date: 20200408 |
|
CF01 | Termination of patent right due to non-payment of annual fee |