CN105223962B - 一种减小外形跳变时弹体过载的动态攻角限幅方法 - Google Patents
一种减小外形跳变时弹体过载的动态攻角限幅方法 Download PDFInfo
- Publication number
- CN105223962B CN105223962B CN201510706349.1A CN201510706349A CN105223962B CN 105223962 B CN105223962 B CN 105223962B CN 201510706349 A CN201510706349 A CN 201510706349A CN 105223962 B CN105223962 B CN 105223962B
- Authority
- CN
- China
- Prior art keywords
- angle
- attack
- missile
- amplitude
- guided missile
- 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.)
- Active
Links
Landscapes
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明公开了一种减小外形跳变时弹体过载的动态攻角限幅方法,通过对导弹攻角进行动态限幅,以减小气动外形跳变时弹体过载;其中对导弹攻角进行动态限幅,实现方式如下:S1、计算导弹在分离前后稳态过程中的最大攻角幅值和气动外形跳变时导弹的最大攻角幅值,以及导弹攻角和飞行速度的变化范围;S2、确定随导弹飞行速度变化的攻角限幅幅值函数;S3、确定最大升力系数与攻角限幅幅值、飞行速度的函数关系;S4、获取导弹参考面积、质量、以及当前时刻导弹飞行动压、飞行速度和弹道倾角,从而确定当前时刻过载限幅器上限值和下限值。实施本发明可根据导弹飞行实时状态确定最佳的攻角限幅幅值,减小气动外形跳变时弹体过载。
Description
技术领域
本发明属于导弹飞行控制技术领域,更具体地,涉及一种减小外形跳变时弹体过载的动态攻角限幅方法。
背景技术
导弹飞行攻角是表征导弹飞行状态的一个重要参量,它的量值大小、正负号以及正负号的变化规律,对导弹的飞行状态和控制系统设计有着重要作用。目前,为了达到一定的气动设计指标、为制导控制专业提供良好的可控飞行器气动外形,通常会在具体的某一状态点附近,通过变化飞行器的局部结构外形,达到改变气动特性的目的。然而,气动外形的跳变,会带来较大的分离干扰(力矩干扰和力干扰),给结构、控制等系统带来不利影响,例如:当导弹飞行马赫数跨度大、飞行空域广,若采用单一气动外形,很难满足制导与稳定专业的可控性要求,为了提高该型导弹的性能指标,以某一马赫数为区分点,分为两种气动外形。在气动外形跳变时刻,分离干扰会产生较大的过载。传统的攻角限幅器算法是采用固定攻角,以最大固定攻角值进行限幅器设计,这种方法无法抑制气动外形跳变时产生的较大过载,不能满足控制系统的要求。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种减小外形跳变时弹体过载的动态攻角限幅方法,方法能够根据导弹飞行实时状态确定最佳的攻角限幅幅值,避免导弹飞行过程中由于气动外形的跳变而影响导弹飞行性能。
为实现上述目的,按照本发明,提供了一种减小外形跳变时弹体过载的动态攻角限幅方法,所述方法通过对导弹攻角进行动态限幅,以减小气动外形跳变时弹体过载;其中,对导弹攻角进行动态限幅,具体实现方式如下:
S1、计算导弹在分离前后稳态过程中的最大攻角幅值αmax1和气动外形跳变时导弹的最大攻角幅值αmax2,以及导弹攻角α和飞行速度Ma的变化范围;
S2、在所述导弹攻角α和飞行速度Ma的变化范围内,利用导弹飞行过程中实时的速度,经过线性插值得到实时的攻角限幅幅值,得到随导弹飞行速度Ma变化的攻角限幅幅值函数αmax=f(Ma),攻角限幅幅值αmax2≤αmax≤αmax1;
S3、确定导弹升力系数CL与攻角α、侧滑角β、飞行速度Ma、舵偏角的函数关系再利用αmax=f(Ma)得到最大升力系数CLmax与攻角限幅幅值αmax、飞行速度Ma的函数关系CLmax=f(αmax,Ma),其中侧滑角为零,舵偏角取与最大攻角αmax对应的平衡舵偏角;
S4、获取导弹参考面积smax、导弹质量m、以及当前时刻导弹飞行动压q、导弹飞行速度v和导弹弹道倾角θ,从而确定当前时刻过载限幅器上限值和过载限幅器下限值
总体而言,通过本发明所构思的以上技术方案与现有技术相比,主要具备以下的技术优点:本发明对导弹攻角进行了动态限幅,减小了气动外形跳变时的弹体攻角,进而对减小了气动外形跳变时的弹体过载。相比于现有技术,避免气动外形跳变时因弹体过载超过预设值而导致的导弹失稳。本发明方法可以在不改变导弹总体参数的情况下,减小气动外形跳变时的弹体过载,避免导弹飞行过程中由于气动外形的跳变而影响导弹飞行性能,提高导弹的控制精度和性能。
附图说明
图1为本发明减小外形跳变时弹体过载的动态攻角限幅方法流程图;
图2为本发明一个实施例中攻角限幅幅值函数曲线图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明提供一种减小外形跳变时弹体过载的动态攻角限幅方法,方法通过对导弹攻角进行动态限幅,以减气动外形跳变时弹体过载;其中,对导弹攻角进行动态限幅,其实现方式如下:
S1、根据原始气动数据和系统稳定性要求,计算得到导弹在分离前后稳态过程中的最大攻角幅值αmax1;其中原始气动数据由气动吹风试验获得。
S2、根据气动外形跳变过程中的干扰力矩、干扰力数据,确定气动外形跳变时导弹的最大攻角αmax2。
可根据气动外形跳变过程中的干扰力矩、干扰力数据,利用数学仿真的方式,通过调整气动外形跳变时的最大攻角,确保跳变过程中弹体最大过载、弹体角速率不超过预设值;当跳变过程中弹体最大过载、弹体角速率达到预设值的临界值时,输出当前时刻的最大攻角αmax2,即确定了气动外形跳变时导弹的最大攻角αmax2。
S3、根据所述原始气动数据的范围,确定导弹攻角α和飞行速度Ma的变化范围。
S4、在所述导弹攻角α和飞行速度Ma的变化范围内,利用导弹飞行过程中实时的速度,经过线性插值得到实时的攻角限幅幅值,得到随导弹飞行Ma变化的攻角限幅幅值函数αmax=f(Ma),αmax2≤αmax≤αmax1。为了导弹状态的平稳过渡,通常在αmax1与αmax2之间取一个中间值αmax3。
S5、确定导弹升力系数CL与攻角α、侧滑角β、飞行速度Ma、舵偏角的函数关系利用αmax=f(Ma)得到最大升力系数CLmax=f(αmax,Ma),其中侧滑角为零,舵偏角取与最大攻角αmax对应的平衡舵偏角;
S6、获取导弹参考面积smax、导弹质量m、以及当前时刻导弹飞行动压q、导弹飞行速度v和导弹弹道倾角θ,从而确定当前时刻过载限幅器上限值和过载限幅器下限值
以下结合一个具体实施例对本发明方法作进一步说明。
(1)根据原始气动数据、稳定系统要求等条件,得到导弹稳态过程中的最大攻角幅值αmax1=12°;
(2)根据气动外形跳变过程中的干扰力矩、干扰力数据,以及导弹自身的稳定性能,确定气动外形跳变时导弹的最大攻角αmax2=6°(设气动外形跳变时导弹速度为V1=4Ma);
(3)为了导弹状态的平稳过渡,在αmax1与αmax2之间取一中间值,设为αmax3=8.4°;
(4)根据气动数据的范围,确定导弹攻角α、飞行速度Ma的变化范围:α的变化区间为[0~12],Ma的变化范围为[1~7];
(5)在以上步骤基础上,可得到随导弹飞行Ma变化的攻角限幅幅值函数αmax=f(Ma),如:αmax1=12°,αmax2=6°,αmax3=8.4°,V1=4Ma,可得攻角限幅器所用攻角曲线如图2所示;
(6)通常,升力系数CL可表示为攻角α、侧滑角β、马赫数Ma、舵偏角的函数,即对于弹道导弹而言,侧滑角产生的升力较小,可以近似将侧滑角取为零,而舵偏角可取配平舵偏角。设计人员可依据这一关系,找出CL随α和Ma变化的函数,即CL=f(α,Ma),以表1的数据为例;
表1随Ma、攻角α(°)变化的升力系数CL
Ma\|α| | 0.0 | 4.0 | 8.0 | 12.0 |
1.20 | 0.0000 | 0.2732 | 0.5737 | 0.9456 |
1.80 | 0.0000 | 0.3118 | 0.6635 | 1.1446 |
2.50 | 0.0000 | 0.3430 | 0.7454 | 1.3012 |
3.00 | 0.0000 | 0.3569 | 0.7680 | 1.2678 |
3.50 | 0.0000 | 0.3401 | 0.7530 | 1.2135 |
5.00 | 0.0000 | 0.3110 | 0.6666 | 1.0451 |
6.00 | 0.0000 | 0.2888 | 0.6159 | 0.9591 |
7.00 | 0.0000 | 0.2803 | 0.5920 | 0.9191 |
(7)在步骤5和步骤6的基础上,可得最大升力系数CLmax=f(αmax,Ma),当导弹速度V=2Ma时,αmax=12°,CLmax=1.1893,当导弹速度V=4Ma时,αmax=6°,CLmax=0.5273,当导弹速度V=6Ma时,αmax=12°,CLmax=0.9591;
(8)将步骤7与导弹飞行实时参数结合,根据可求得实时过载限幅器上限值,根据可求得实时过载限幅器下限值。通过对导弹攻角进行动态限幅,从而减小了气动外形跳变时弹体过载。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种减小外形跳变时弹体过载的动态攻角限幅方法,其特征在于,所述方法通过对导弹攻角进行动态限幅,以减小气动外形跳变时弹体过载;其中,对导弹攻角进行动态限幅,具体实现方式如下:
S1、计算导弹在分离前后稳态过程中的最大攻角幅值αmax1和气动外形跳变时导弹的最大攻角幅值αmax2,以及导弹攻角α和飞行速度Ma的变化范围;
S2、在所述导弹攻角α和飞行速度Ma的变化范围内,利用导弹飞行过程中实时的速度,经过线性插值得到实时的攻角限幅幅值,得到随导弹飞行速度Ma变化的攻角限幅幅值函数αmax=f(Ma),攻角限幅幅值αmax2≤αmax≤αmax1;
S3、确定导弹升力系数CL与攻角α、侧滑角β、飞行速度Ma、舵偏角的函数关系再利用αmax=f(Ma)得到最大升力系数CLmax与攻角限幅幅值αmax、飞行速度Ma的函数关系CLmax=f(αmax,Ma),其中侧滑角为零,舵偏角取与最大攻角αmax对应的平衡舵偏角;
S4、获取导弹参考面积smax、导弹质量m、以及当前时刻导弹飞行动压q、导弹飞行速度v和导弹弹道倾角θ,从而确定当前时刻过载限幅器上限值和过载限幅器下限值
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510706349.1A CN105223962B (zh) | 2015-10-27 | 2015-10-27 | 一种减小外形跳变时弹体过载的动态攻角限幅方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510706349.1A CN105223962B (zh) | 2015-10-27 | 2015-10-27 | 一种减小外形跳变时弹体过载的动态攻角限幅方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105223962A CN105223962A (zh) | 2016-01-06 |
CN105223962B true CN105223962B (zh) | 2017-09-01 |
Family
ID=54992992
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510706349.1A Active CN105223962B (zh) | 2015-10-27 | 2015-10-27 | 一种减小外形跳变时弹体过载的动态攻角限幅方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105223962B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107121015B (zh) * | 2017-06-16 | 2018-10-16 | 湖北航天技术研究院总体设计所 | 一种快速弹上弹道在线规划方法 |
CN107390718B (zh) * | 2017-08-04 | 2020-09-18 | 中国运载火箭技术研究院 | 一种高速飞行强适应性180度翻转控制方法 |
CN112393873B (zh) * | 2021-01-11 | 2021-04-06 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种攻角匹配方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3963197A (en) * | 1971-12-16 | 1976-06-15 | Messerschmitt-Bolkow-Blohm Gmbh | Control device for avoiding the pitching up of missiles or aircraft |
RU2248304C2 (ru) * | 2002-09-20 | 2005-03-20 | Макеев Роман Владимирович | Способ и устройство ограничения угла атаки и перегрузки самолета |
CN103587681A (zh) * | 2013-10-16 | 2014-02-19 | 中国运载火箭技术研究院 | 抑制侧滑角信号常值偏差影响的高超声速飞行器控制方法 |
CN103863579A (zh) * | 2014-03-31 | 2014-06-18 | 北京控制工程研究所 | 一种深空探测返回过程的预测校正制导方法 |
CN104808492A (zh) * | 2015-03-23 | 2015-07-29 | 北京航天自动控制研究所 | 一种升力式飞行器的攻角指令生成方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100298260B1 (ko) * | 1998-05-27 | 2002-10-19 | 한국항공우주산업 주식회사 | 항공기의받음각제어방법 |
CZ291056B6 (cs) * | 2000-09-26 | 2002-12-11 | Peter Zábranský | Způsob vytváření simulace přetíľení a mezních aerodynamických reľimů pro piloty a zařízení k provádění tohoto způsobu |
-
2015
- 2015-10-27 CN CN201510706349.1A patent/CN105223962B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3963197A (en) * | 1971-12-16 | 1976-06-15 | Messerschmitt-Bolkow-Blohm Gmbh | Control device for avoiding the pitching up of missiles or aircraft |
RU2248304C2 (ru) * | 2002-09-20 | 2005-03-20 | Макеев Роман Владимирович | Способ и устройство ограничения угла атаки и перегрузки самолета |
CN103587681A (zh) * | 2013-10-16 | 2014-02-19 | 中国运载火箭技术研究院 | 抑制侧滑角信号常值偏差影响的高超声速飞行器控制方法 |
CN103863579A (zh) * | 2014-03-31 | 2014-06-18 | 北京控制工程研究所 | 一种深空探测返回过程的预测校正制导方法 |
CN104808492A (zh) * | 2015-03-23 | 2015-07-29 | 北京航天自动控制研究所 | 一种升力式飞行器的攻角指令生成方法 |
Non-Patent Citations (2)
Title |
---|
林献武,等.超远程弹箭降弧段动力平衡角的快速计算方法.《火力与指挥控制》.2010,第35卷(第9期),113-116,123. * |
王宏然,等.模糊反演飞航导弹过载控制系统设计.《弹箭与制导学报》.2006,第26卷(第2期),10-12. * |
Also Published As
Publication number | Publication date |
---|---|
CN105223962A (zh) | 2016-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2974955B1 (en) | Closed loop control of aircraft control surfaces | |
EP3741679B1 (en) | Multi-engine aircraft thrust balancing | |
CN105223962B (zh) | 一种减小外形跳变时弹体过载的动态攻角限幅方法 | |
CN102707723B (zh) | 基于飞行器常规模型的横航向控制器区域设计方法 | |
CN110825115B (zh) | 一种飞机迎角和过载的极限限制控制方法 | |
CN111240212A (zh) | 一种基于优化预测的倾转旋翼无人机控制分配方法 | |
CN106647264B (zh) | 一种基于控制约束的扩展鲁棒h∞的无人机控制方法 | |
CN102707616B (zh) | 基于飞行器三角模型的控制器区域设计方法 | |
CN113848963A (zh) | 一种飞行控制系统的控制律参数设计方法 | |
CN109521673A (zh) | 一种二维机翼颤振问题的区间滑模抑制方法 | |
CN109033487B (zh) | 一种基于蒙特卡洛仿真的飞行器总体概率设计方法 | |
CN109472073B (zh) | 一种飞行器气动布局调整方法、装置及电子设备 | |
CN102707722A (zh) | 基于飞行器常规模型的全维控制器区域设计方法 | |
CN102707629A (zh) | 基于飞行器切换模型的全维控制器区域设计方法 | |
CN109117584A (zh) | 一种低速飞机突风载荷系数计算方法及设备 | |
CN103792848B (zh) | 纵向飞行模型簇人机闭环复合根轨迹多级pid鲁棒控制器设计方法 | |
CN112363521B (zh) | 一种适应滚动通道气动非线性的控制参数调节算法 | |
CN103823367B (zh) | 纵向飞行模型簇颤振抑制复合频率鲁棒控制器设计方法 | |
CN104898685B (zh) | 一种基于概率最优的飞行器俯仰通道攻角跟踪控制方法 | |
CN103809452B (zh) | 纵向飞行模型簇颤振抑制复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103823375B (zh) | 纵向飞行模型簇复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103823374B (zh) | 飞行器多回路模型簇复合根轨迹补偿控制器设计方法 | |
Gao et al. | Three‐Dimensional Integrated Guidance and Control for Near Space Interceptor Based on Robust Adaptive Backstepping Approach | |
CN103809447B (zh) | 飞行器多回路模型簇复合频率控制器设计方法 | |
Pervushina et al. | Designing an adaptive stabilizing system for an unmanned aerial vehicle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |