[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN105843232A - 一种飞行器滑翔减速控制方法 - Google Patents

一种飞行器滑翔减速控制方法 Download PDF

Info

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
Application number
CN201610217768.3A
Other languages
English (en)
Other versions
CN105843232B (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.)
Beijing Aerospace Automatic Control Research Institute
Original Assignee
Beijing Aerospace Automatic Control Research Institute
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 Beijing Aerospace Automatic Control Research Institute filed Critical Beijing Aerospace Automatic Control Research Institute
Priority to CN201610217768.3A priority Critical patent/CN105843232B/zh
Publication of CN105843232A publication Critical patent/CN105843232A/zh
Application granted granted Critical
Publication of CN105843232B publication Critical patent/CN105843232B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/04Control of altitude or depth
    • G05D1/06Rate of change of altitude or depth
    • G05D1/0607Rate of change of altitude or depth specially adapted for aircraft
    • G05D1/0653Rate of change of altitude or depth specially adapted for aircraft during a phase of take-off or landing
    • G05D1/0676Rate 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
&alpha; 0 = F n q &times; S m - a x b x ;
其中,ax、bx分别为与飞行器法向力气动特性相关的参数,q为动压,Sm为飞行器特征面积,Fn为无减速控制时的制导力。
较佳的,根据如下的公式计算得到需要耗散掉的速度ΔV:
&Delta; V = V d - V &alpha; T g ;
其中,Tg为距离飞行结束的待飞时间。
较佳的,根据如下的公式计算得到减速需用攻角αn
&alpha; n = &lsqb; ( &alpha; 0 ) N + K V 1 &CenterDot; &Delta; V &CenterDot; m q &CenterDot; S m &CenterDot; b D &rsqb; 1 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
&alpha; 0 = F n q &times; S m - a x b x - - - ( 1 )
其中,ax、bx分别为与飞行器法向力气动特性相关的第一和第二参数,ax、bx可以根据飞行器的气动特性拟合得到;q为动压,Sm为飞行器特征面积,Fn为无减速控制时的制导力。
在本发明的技术方案中,所述无减速控制时的制导力Fn可以根据现有技术中的常用制导算法计算得到,因此,具体计算过程在此不再赘述。
步骤103,计算得到需要耗散掉的速度ΔV;
在本发明的技术方案中,可以根据多种具体的实现方式来计算得到需要耗散掉的速度ΔV。以下将以其中的一种具体实现方式为例,对本发明的技术方案进行说明。
例如,较佳的,在本发明的一个具体实施例中,可以根据Vcx、Vd和距离飞行结束的待飞时间,计算需要耗散掉的速度ΔV。例如,可以根据如下所述的公式计算得到需要耗散掉的速度ΔV:
&Delta; V = V d - V c x T g - - - ( 2 )
其中,Tg为距离飞行结束的待飞时间。在本发明的技术方案中,可以根据实际情况,预先设置Tg的取值,即设Tg为一个常值;或者,也可以根据现有技术中的常用制导算法计算得到该Tg,具体计算过程在此不再赘述。
步骤104,计算得到减速需用攻角αn
在本发明的技术方案中,可以根据多种具体的实现方式来计算得到减速需用攻角αn。以下将以其中的一种具体实现方式为例,对本发明的技术方案进行说明。
例如,较佳的,在本发明的一个具体实施例中,可以根据如下所述的公式计算得到减速需用攻角αn
&alpha; n = &lsqb; ( &alpha; 0 ) N + K V 1 &CenterDot; &Delta; V &CenterDot; m q &CenterDot; S m &CenterDot; b D &rsqb; 1 N - - - ( 3 )
其中,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
&alpha; 0 = F n q &times; S m - a x b x ;
其中,ax、bx分别为与飞行器法向力气动特性相关的参数,q为动压,Sm为飞行器特征面积,Fn为无减速控制时的制导力。
4.根据权利要求3所述的方法,其特征在于,根据如下的公式计算得到需要耗散掉的速度ΔV:
&Delta; V = V d - V c x T g ;
其中,Tg为距离飞行结束的待飞时间。
5.根据权利要求4所述的方法,其特征在于,根据如下的公式计算得到减速需用攻角αn
&alpha; n = &lsqb; ( &alpha; 0 ) N + K V 1 &CenterDot; &Delta; V &CenterDot; m q &CenterDot; S m &CenterDot; b D &rsqb; 1 N ;
其中,bD为与飞行器阻力系数相关的参数;KV1、N为制导参数;m为飞行器质量。
6.根据权利要求5所述的方法,其特征在于,所述计算得到减速控制附加的制导力包括:
计算得到Fn0=(ax+bx·α0)·q·Sm
计算得到Fnn=(ax+bx·αn)·q·Sm
如果Fnn>Fn0,则否则,
CN201610217768.3A 2016-04-08 2016-04-08 一种飞行器滑翔减速控制方法 Expired - Fee Related CN105843232B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 北京航天自动控制研究所 一种升力式飞行器的攻角指令生成方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
刘礅等: "《空间飞行器动力学》", 31 October 2003, 哈尔滨工业大学出版社 *
马宏图: "高可靠性再入轨迹在线设计及实时制导", 《CNKI硕士学位论文》 *

Cited By (6)

* Cited by examiner, † Cited by third party
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