CN103576554B - 基于控制需求的飞行器气动误差模型分量、分级设计方法 - Google Patents

基于控制需求的飞行器气动误差模型分量、分级设计方法 Download PDF

Info

Publication number
CN103576554B
CN103576554B CN201310549486.XA CN201310549486A CN103576554B CN 103576554 B CN103576554 B CN 103576554B CN 201310549486 A CN201310549486 A CN 201310549486A CN 103576554 B CN103576554 B CN 103576554B
Authority
CN
China
Prior art keywords
deviation
rudder
rolling
pitching
moment coefficient
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
Application number
CN201310549486.XA
Other languages
English (en)
Other versions
CN103576554A (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.)
China Academy of Launch Vehicle Technology CALT
Beijing Institute of Near Space Vehicles System Engineering
Original Assignee
China Academy of Launch Vehicle Technology CALT
Beijing Institute of Near Space Vehicles System Engineering
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 China Academy of Launch Vehicle Technology CALT, Beijing Institute of Near Space Vehicles System Engineering filed Critical China Academy of Launch Vehicle Technology CALT
Priority to CN201310549486.XA priority Critical patent/CN103576554B/zh
Publication of CN103576554A publication Critical patent/CN103576554A/zh
Application granted granted Critical
Publication of CN103576554B publication Critical patent/CN103576554B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明涉及一种基于控制需求的飞行器气动误差模型分量、分级设计方法,依次包括:一、实时采集参数;二、获得俯仰舵偏产生的俯仰力矩系数增量;获得滚动舵偏产生的俯仰力矩系数增量、偏航力矩系数增量、滚动力矩系数增量;三、获得俯仰常值项偏差;俯仰舵效项偏差;滚动舵偏诱发的俯仰力矩系数偏差;偏航常值项偏差;偏航稳定项偏差;滚动舵偏诱发的偏航力矩系数偏差;滚动常值项偏差;滚动稳定项偏差;滚动舵效项偏差;四、建立俯仰力矩系数偏差;偏航力矩系数偏差;滚动力矩系数偏差;五、划分三级气动误差模型。本发明可平衡总体设计中各分系统的设计难度,加快总体方案闭合,充分预示飞行试验风险点,为飞行试验后气动误差量值缩小提供参考。

Description

基于控制需求的飞行器气动误差模型分量、分级设计方法
技术领域
本发明涉及一种基于控制需求的飞行器气动误差模型分量、分级设计方法,特别是涉及一种实现控制特性恶劣模式拉偏组合,加强地面预示及控制能力的基于控制需求的飞行器气动误差模型分量、分级设计方法。
背景技术
气动误差模型是控制系统分析、设计和仿真中误差状态选取的最重要组成部分,气动误差模型的构架、量值及使用方式决定了控制系统设计结果的可靠性,因此需要基于控制需求进行气动误差模型设计及使用的研究。
对于高超声速类滑翔飞行器,主要依靠气动力控制实现滑翔飞行,需要对气动误差模型进行精细刻画。传统的气动误差模型形式为“常值项+比例项”的形式,以法向力系数偏差为例说明,一般为|ΔCN|≤(0.05+15%·|CN|)形式所示。此种形式存在一定的问题,即无法满足控制系统设计对稳定性、舵效、交叉导数等气动增量项的拉偏需求,无法实现控制系统关键参数误差极限状态的考核。为充分激发高超声速滑翔飞行器操纵性与稳定性更为恶劣的模式,以加强地面预示及状态控制的能力。根据飞行器本体特性创新性提出一种分量式气动误差模型,此种模型形式将偏差量划分为“常值项+稳定项+舵效项+控制耦合项”的分量形式进行表述,可实现各分量不同极性的拉偏,形成恶劣模式的组合状态。
在气动误差模型已建立的基础上,合理有效的使用及考核方式,一方面可平衡总体设计中各分系统的设计难度,加快总体方案闭合,另外,可充分预示飞行试验风险点,并为飞行试验后气动误差量值的缩小提供参考。对于气动误差模型的使用及考核方式,一般的做法是针对一级误差带,进行极限拉偏组合开展仿真,要求全偏差组合状态下均能实现稳定飞行,此种做法经常面临的问题是无法实现误差源设计与控制设计的平衡,方案难以闭合。为实现方案的闭合,气动误差模型需要极尽所能的缩小,控制能力需要极尽所能的增强,其后果为加大了气动误差地面预示不足的风险,而对于主要依靠空气舵实现滑翔飞行的飞行器,则存在空气舵控制能力挖掘不充分的情况(这里尤指配有姿控动力系统的飞行器,即为适应较大的气动误差,在空气舵控制能力尚未充分发挥作用的情况下,即采用了RCS进行辅助控制)。
基于以上气动误差模型在控制应用中存在的两方面问题,因此亟需提供一种新型的基于控制需求的飞行器气动误差模型分量、分级设计方法。
发明内容
本发明要解决的技术问题是提供一种实现控制特性恶劣模式拉偏组合,加强地面预示及控制能力的,为飞行试验后气动误差大小的收敛建立攻关目标,平衡气动误差模型实现难度和姿控设计难度,提前预示气动误差不确定性带来的稳定飞行风险的基于控制需求的飞行器气动误差模型分量、分级设计方法。
为解决上述技术问题,本发明一种基于控制需求的飞行器气动误差模型分量、分级设计方法,依次包括以下步骤:
第一步、实时采集飞行高度Hi,飞行马赫数Ma,飞行攻角α,飞行侧滑角β,飞行器俯仰舵偏飞行器偏航舵偏δψ,飞行器滚动舵偏δγ
根据误差带包络确定常值系数a、b、k;
第二步、获得俯仰舵偏产生的力矩系数获得0°俯仰舵偏下的力矩系数相减得到俯仰舵偏产生的俯仰力矩系数增量
获得滚动舵偏产生的力矩系数获得0°滚动舵偏下的力矩系数相减得到滚动舵偏产生的俯仰力矩系数增量滚动舵偏产生的偏航力矩系数增量滚动舵偏产生的滚动力矩系数增量
第三步、获得俯仰常值项偏差ΔCmz0=±(a*α+b);
俯仰舵效项偏差
滚动舵偏诱发的俯仰力矩系数偏差
偏航常值项偏差ΔCmy0=±(a*α+b);
偏航稳定项偏差 ΔC my β = ± ( a * α + b ) ;
滚动舵偏诱发的偏航力矩系数偏差
滚动常值项偏差ΔCmx0=±(a*α+b);
滚动稳定项偏差 ΔC mx β = ± ( a * α + b ) ;
滚动舵效项偏差
第四步、建立俯仰力矩系数偏差
建立偏航力矩系数偏差 ΔC my = ΔC my 0 + ΔC my β · β + ΔC my _ δ γ ;
建立滚动力矩系数偏差 ΔC mx = ΔC mx 0 + ΔC mx β · β + ΔC mx _ δ γ ;
第五步、当空气舵起控后,依靠空气舵独立控制即可实现全偏差状态稳定,此情况下可适应的气动误差模型为第一级气动误差模型;
在第一级误差模型基础上对各分量项进行不同程度的放大,当空气舵起控后,依靠RCS与空气舵复合控制实现全偏差状态稳定,此情况下可适应的气动误差模型为第二级气动误差模型;
在第二级误差模型基础上对各分量项进行不同程度的放大,当空气舵起控后,依靠RCS与空气舵复合控制,对全偏差组合状态进行稳定概率摸底。
采用本发明可以挖掘出高超声速滑翔飞行器布局评价或控制设计更为恶劣的状态,实现控制特性恶劣拉偏组合,充分挖掘地面设计风险和控制能力;
本发明突破传统姿控设计中偏差使用的思路,使设计不过于保守和冒进,为气动误差的设计和姿控设计均提出了攻关目标,并平衡了气动误差模型实现难度和姿控设计难度,一方面可平衡总体设计中各分系统的设计难度,加快总体方案闭合,另外,可充分预示飞行试验风险点,并为飞行试验后气动误差量值的缩小提供参考。
具体实施方式
对于俯仰通道,常值项反应纵向压心预示偏差。对于偏航、滚动通道,常值项反应结构制造、安装、飞行中烧蚀等因素产生的外形不对称偏差。对于俯仰通道,稳定性项反应俯仰静稳定性偏差。对于偏航、滚动通道,稳定性项偏差分别反应偏航静稳定性偏差、滚动静稳定性偏差。对于舵效项偏差,分别反应三通道气动舵能力偏差。对于控制耦合项偏差,分别反应某一通道气动舵诱发的其它通道气动力矩大小的偏差。
对于俯仰通道,一般需考虑纵向压心偏差,舵效偏差,因俯仰舵偏与滚动舵偏共用升降舵偏,滚动舵偏诱发的俯仰力矩系数变化不可忽略,所以需特别考虑滚动诱发的俯仰力矩系数偏差项。
对于偏航通道,一般需考虑常值项偏差、稳定项偏差、舵效偏差,对于具有较大后缘舵布局的飞行器,因滚动舵偏诱发的偏航力矩系数变化不可忽略,所以需特别考虑滚动诱发的偏航力矩系数偏差项。
具体来说,本发明依次包括以下步骤:
第一步、实时采集飞行高度Hi,飞行马赫数Ma,飞行攻角α,飞行侧滑角β,飞行器俯仰舵偏飞行器偏航舵偏δψ,飞行器滚动舵偏δγ
根据误差带包络确定常值系数a、b、k;
第二步、获得俯仰舵偏产生的力矩系数获得0°俯仰舵偏下的力矩系数相减得到俯仰舵偏产生的俯仰力矩系数增量
获得滚动舵偏产生的力矩系数获得0°滚动舵偏下的力矩系数相减得到滚动舵偏产生的俯仰力矩系数增量滚动舵偏产生的偏航力矩系数增量滚动舵偏产生的滚动力矩系数增量
第三步、对于常值项偏差ΔCmz0,其偏差的变化特性表现为攻角的一次线性函数,在压心前后移动量一致的情况下,获得俯仰常值项偏差ΔCmz0=±(a*α+b);
对于舵效项偏差,可用舵偏增量百分比的形式,获得俯仰舵效项偏差
对于滚动舵偏诱发的俯仰力矩系数偏差项,根据对飞行器本体特性的认识,其诱发的偏差大小与俯仰舵偏位置和滚动舵产生的俯仰力矩系数增量值相关,获得滚动舵偏诱发的俯仰力矩系数偏差
偏航常值项偏差ΔCmy0=±(a*α+b);
偏航稳定项偏差 ΔC my β = ± ( a * α + b ) ;
滚动舵偏诱发的偏航力矩系数偏差
滚动常值项偏差ΔCmx0=±(a*α+b);
滚动稳定项偏差 ΔC mx β = ± ( a * α + b ) ;
滚动舵效项偏差
第四步、建立俯仰力矩系数偏差
在偏航通道无方向舵控制的情况下,建立偏航力矩系数偏差 ΔC my = ΔC my 0 + ΔC my β · β + ΔC my _ δ γ ;
建立滚动力矩系数偏差 ΔC mx = ΔC mx 0 + ΔC mx β · β + Δ C mx _ δ γ ;
第五步、在误差模型构架确定的基础上,综合控制能力可适应的气动误差量值及仿真方式,确定相应级别的误差带量值,将误差带分为三级进行设计和考核。
一:首先确定第一级气动误差模型。综合考虑气动偏差及其它设计偏差,进行全偏差极限拉偏组合进行仿真,当空气舵起控后,依靠空气舵独立控制即可实现全偏差状态稳定,此情况下可适应的气动误差模型为第一级气动误差模型;
二:在第一级气动误差模型基础上,确定第二级气动误差模型,此级误差模型在第一级误差模型基础上对各分量项进行不同程度的放大得到。综合考虑此级气动偏差及其它设计偏差,进行全偏差极限拉偏组合进行仿真,当空气舵起控后,依靠RCS与空气舵复合控制实现全偏差状态稳定,此情况下可适应的气动误差模型为第二级气动误差模型;
三:在第二级气动误差模型基础上,确定第三级气动误差模型,此级误差模型在第二级误差模型基础上对各分量项进行不同程度的放大得到。综合考虑此级气动偏差及其它设计偏差,对气动偏差进行抽样取值,对其它偏差进行极限取值,当空气舵起控后,依靠RCS与空气舵复合控制,对全偏差组合状态进行稳定概率摸底。

Claims (1)

1.一种基于控制需求的飞行器气动误差模型分量、分级设计方法,依次包括以下步骤:
第一步、实时采集飞行高度Hi,飞行马赫数Ma,飞行攻角α,飞行侧滑角β,飞行器俯仰舵偏飞行器偏航舵偏δψ,飞行器滚动舵偏δγ
根据误差带包络确定常值系数a、b、k;
第二步、获得俯仰舵偏产生的力矩系数获得0°俯仰舵偏下的力矩系数相减得到俯仰舵偏产生的俯仰力矩系数增量
获得滚动舵偏产生的力矩系数获得0°滚动舵偏下的力矩系数相减得到滚动舵偏产生的俯仰力矩系数增量滚动舵偏产生的偏航力矩系数增量滚动舵偏产生的滚动力矩系数增量
第三步、获得俯仰常值项偏差ΔCmz0=±(a*α+b);
俯仰舵效项偏差
滚动舵偏诱发的俯仰力矩系数偏差
偏航常值项偏差ΔCmy0=±(a*α+b);
偏航稳定项偏差 ΔC m y β = ± ( a * α + b ) ;
滚动舵偏诱发的偏航力矩系数偏差
滚动常值项偏差ΔCmx0=±(a*α+b);
滚动稳定项偏差 ΔC m x β = ± ( a * α + b ) ;
滚动舵效项偏差
第四步、建立俯仰力矩系数偏差
建立偏航力矩系数偏差 ΔC m y = ΔC m y 0 + ΔC m y β · β + ΔC m y _ δ γ ;
建立滚动力矩系数偏差 ΔC m x = ΔC m x 0 + ΔC m x β · β + ΔC m x _ δ γ ;
第五步、当空气舵起控后,依靠空气舵独立控制即可实现全偏差状态稳定,此情况下可适应的气动误差模型为第一级气动误差模型;
在第一级误差模型基础上对各分量项进行不同程度的放大,当空气舵起控后,依靠RCS与空气舵复合控制实现全偏差状态稳定,此情况下可适应的气动误差模型为第二级气动误差模型;
在第二级误差模型基础上对各分量项进行不同程度的放大,当空气舵起控后,依靠RCS与空气舵复合控制,对全偏差组合状态进行稳定概率摸底。
CN201310549486.XA 2013-11-07 2013-11-07 基于控制需求的飞行器气动误差模型分量、分级设计方法 Active CN103576554B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310549486.XA CN103576554B (zh) 2013-11-07 2013-11-07 基于控制需求的飞行器气动误差模型分量、分级设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310549486.XA CN103576554B (zh) 2013-11-07 2013-11-07 基于控制需求的飞行器气动误差模型分量、分级设计方法

Publications (2)

Publication Number Publication Date
CN103576554A CN103576554A (zh) 2014-02-12
CN103576554B true CN103576554B (zh) 2016-05-18

Family

ID=50048571

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310549486.XA Active CN103576554B (zh) 2013-11-07 2013-11-07 基于控制需求的飞行器气动误差模型分量、分级设计方法

Country Status (1)

Country Link
CN (1) CN103576554B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104155983B (zh) * 2014-08-08 2015-05-20 北京航天自动控制研究所 飞行器姿态运动通道间气动耦合特性的交联影响确定方法
CN104155985B (zh) * 2014-08-08 2015-05-20 北京航天自动控制研究所 飞行器姿态运动通道间惯性耦合特性的交联影响确定方法
CN105628333B (zh) * 2014-10-29 2018-06-22 北京临近空间飞行器系统工程研究所 一种高空高马赫数条件下的气动误差确定方法
CN104914736A (zh) * 2015-05-07 2015-09-16 南京航空航天大学 一种高超声速飞行器通用仿真模型的构建方法
CN105136422B (zh) * 2015-09-10 2017-10-13 中国航天空气动力技术研究院 风洞试验中修正飞行器模型侧滑弹性角的方法
CN111985169B (zh) * 2020-08-26 2024-04-02 北京中科宇航技术有限公司 一种近轴对称飞行器气动特性建模方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4970638A (en) * 1988-10-31 1990-11-13 The United States Of America As Represented By The Secretary Of The Air Force Control of unknown systems via deconvolution
CN101660914A (zh) * 2009-08-19 2010-03-03 南京航空航天大学 耦合惯性位置误差的机载星光和惯性组合的自主导航方法
CN103063216A (zh) * 2013-01-06 2013-04-24 南京航空航天大学 一种基于星像坐标建模的惯性与天文组合导航方法
CN103245359A (zh) * 2013-04-23 2013-08-14 南京航空航天大学 一种惯性导航系统中惯性传感器固定误差实时标定方法
CN103267531A (zh) * 2013-04-23 2013-08-28 上海卫星工程研究所 一种光纤陀螺随机误差的高精度补偿方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005015324A2 (en) * 2002-12-09 2005-02-17 Georgia Tech Research Corp. Adaptive output feedback apparatuses and methods capable of controlling a non­minimum phase system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4970638A (en) * 1988-10-31 1990-11-13 The United States Of America As Represented By The Secretary Of The Air Force Control of unknown systems via deconvolution
CN101660914A (zh) * 2009-08-19 2010-03-03 南京航空航天大学 耦合惯性位置误差的机载星光和惯性组合的自主导航方法
CN103063216A (zh) * 2013-01-06 2013-04-24 南京航空航天大学 一种基于星像坐标建模的惯性与天文组合导航方法
CN103245359A (zh) * 2013-04-23 2013-08-14 南京航空航天大学 一种惯性导航系统中惯性传感器固定误差实时标定方法
CN103267531A (zh) * 2013-04-23 2013-08-28 上海卫星工程研究所 一种光纤陀螺随机误差的高精度补偿方法

Also Published As

Publication number Publication date
CN103576554A (zh) 2014-02-12

Similar Documents

Publication Publication Date Title
CN103576554B (zh) 基于控制需求的飞行器气动误差模型分量、分级设计方法
CN105752314B (zh) 一种高空低速自然层流高升力翼型
CN103587680B (zh) 一种飞行器侧滑转弯控制方法
CN103587681A (zh) 抑制侧滑角信号常值偏差影响的高超声速飞行器控制方法
CN103625637B (zh) 一种大型飞机侧向阵风缓和方法
CN102707624A (zh) 基于飞行器常规模型的纵向控制器区域设计方法
CN103287574A (zh) 一种飞机增升装置控制方法
CN102117362B (zh) 滑流影响下的轻型飞机水平尾翼设计载荷确定方法
CN107330152A (zh) 一种适用于旋翼飞行器的高效气动配平方法
CN102052266B (zh) 基于尖尾缘翼型设计的后加载钝尾缘翼型
CN104401503B (zh) 一种适用于多舵面再入飞行器的配平舵偏角确定方法
Hufford et al. Viscous flow around marine propellers using boundary-layer strip theory
Shindo Simplified tunnel correction method
Lamar Prediction of vortex flow characteristics of wings at subsonic and supersonic speeds
CN104155985B (zh) 飞行器姿态运动通道间惯性耦合特性的交联影响确定方法
RU2352914C1 (ru) Аэромеханический способ измерения параметров состояния вертолета в полете и устройство для его осуществления
CN106168529B (zh) 一种修正飞行器刚性体升力系数的风洞试验方法
CN106950982A (zh) 再入飞行器姿控动力系统高空力矩特性辨识方法
CN112541305A (zh) 一种基于全局变量求导的飞机气动力特性分析方法
CN110160737A (zh) 基于工程面元法的弹性气动数据精确获取方法
Huffman et al. Subsonic longitudinal and lateral-directional static aerodynamic characteristics for a model with swept back and swept forward wings
Nichols Analysis and compilation of missile aerodynamic data. Volume 1: Data presentation and analysis
Rein Measurements of aerodynamic forces and moments on the DLR-F17E model in low-and high-speed flows
Ferris Static aerodynamic characteristics of a model with a 17 percent thick supercritical wing
Spreiter Downwash and sidewash fields behind cruciform wings

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant