CN114234737A - 基于攻角变化的姿态控制系统设计方法 - Google Patents

基于攻角变化的姿态控制系统设计方法 Download PDF

Info

Publication number
CN114234737A
CN114234737A CN202111467689.5A CN202111467689A CN114234737A CN 114234737 A CN114234737 A CN 114234737A CN 202111467689 A CN202111467689 A CN 202111467689A CN 114234737 A CN114234737 A CN 114234737A
Authority
CN
China
Prior art keywords
angle
attack
flight
attack angle
rocket
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
CN202111467689.5A
Other languages
English (en)
Other versions
CN114234737B (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.)
CASIC Rocket Technology Co
Original Assignee
CASIC Rocket Technology Co
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 CASIC Rocket Technology Co filed Critical CASIC Rocket Technology Co
Priority to CN202111467689.5A priority Critical patent/CN114234737B/zh
Publication of CN114234737A publication Critical patent/CN114234737A/zh
Application granted granted Critical
Publication of CN114234737B publication Critical patent/CN114234737B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42AMMUNITION; BLASTING
    • F42BEXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
    • F42B15/00Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
    • F42B15/01Arrangements thereon for guidance or control

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Combustion & Propulsion (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本申请涉及一种基于攻角变化的姿态控制系统设计方法,包括以下步骤:计算飞行攻角包络;设计各时刻点下各攻角状态的校正网络;对每一特征时刻点重复步骤二,将各状态下的校正网络组合为参数矩阵装订于飞控软件;火箭飞行过程中实时辨识飞行攻角;利用当前飞行时刻及实时攻角插值校正网络矩阵用于姿态控制。本申请具有以下可预期的技术效果:在飞行环境不确定性较大的情况下,能够实时辨识火箭飞行过程中的攻角变化,根据攻角大小进行校正网络调整,提高姿态控制系统的鲁棒性以及控制性能,姿态控制精度也相对较好,对误差的适应能力较好。

Description

基于攻角变化的姿态控制系统设计方法
技术领域
本申请涉及火箭制导技术领域,尤其是涉及一种基于攻角变化的姿态控制系统设计方法。
背景技术
随着商业发射需求的增加,由于可回收液体运载火箭具有提高利用效率、降低飞行成本的优势,国内外研究机构及商业航天公司均开展了相应的研制工作。
可回收火箭在返回过程中发动机朝下,与传统的低阻力流线体飞行器气动特性差异较大,存在气动线性度差的问题,即不同攻角状态下气动力和力矩系数对攻角的导数变化较大。在进行姿态控制系统设计时,这一气动特性将带来较大的参数误差。
采用传统设计方法对误差的适应能力有限,在飞行环境不确定性较大的情况下,控制性能和姿态控制精度相对较差,故而有待改进。
发明内容
为了改善传统设计方法的控制性能和姿态控制精度相对较差等问题,本申请提供一种基于攻角变化的姿态控制系统设计方法。
本申请提供一种基于攻角变化的姿态控制系统设计方法,采用如下的技术方案:
一种基于攻角变化的姿态控制系统设计方法,包括以下步骤:
步骤一,计算飞行攻角包络;
步骤二,设计各时刻点下各攻角状态的校正网络;
步骤三,对每一特征时刻点重复步骤二,将各状态下的校正网络组合为参数矩阵装订于飞控软件;
步骤四,火箭飞行过程中实时辨识飞行攻角;
步骤五,利用当前飞行时刻及实时攻角插值校正网络矩阵用于姿态控制。
进一步地,步骤一的具体操作如下:
根据总体提供的风场数据使用要求,计算火箭飞行过程中的攻角包络,对于任一特征时刻点ti,存在攻角变化范围[0,αimax],其中i取值为1到n,其中n表示特征时刻点个数,n的取值为正整数,特征时刻点ti的选取特征依据包括但不限于:火箭的飞行弹道和气动特性。
进一步地,步骤二的具体操作如下:
对特征时刻点ti,针对攻角变化范围在[0,αimax]内的每一个状态各设计一组校正网络参数Gij,其中j取值为1到m,m的取值为正整数,Gij的下角标i表示第i个特征时刻点,Gij的下角标j表示第j个攻角状态点,对每一个特征时刻点,攻角变化范围不同,选取的攻角状态点个数应该相同,即对每个时间点m为定值。
进一步地,步骤三的具体操作如下:
对每一个特征时刻点重复步骤二直至完成该飞行时段内所有特征点校正网络设计,将每一个特征时刻点的每一个攻角状态下的校正网络参数组合为一个mxkxn的矩阵Ejz,将矩阵Ejz装订至飞控软件中,其中k表示任一状态下的校正网络参数的长度,且对每一个状态k值保持不变。
进一步地,步骤四的具体操作如下:
飞行过程中,从火箭启控开始实时辨识飞行攻角,时刻tri的攻角辨识方法如下:通过火箭惯组反馈的箭体过载ny进行计算,具体计算公式如下:
Figure BDA0003392227000000031
式中,m为箭体质量,CN为法向力系数,q为动压,Sm为参考面积,P为摆动发动机推力,
Figure BDA0003392227000000032
为发动机俯仰通道摆角,
Figure BDA0003392227000000033
为俯仰角加速度,xz为质心位置,xa为火箭惯组安装位置;
在箭体过载ny已知的前提下,根据上述式可以求出时刻tri的法向力系数CN,由于CN为随攻角变化的参数,以CN插值法向力系数矩阵,即可得到时刻tri的飞行攻角αri
进一步地,用于攻角辨识的箭体过载ny选取为h拍内的平均值,h的取值可以通过仿真试验进行确定。
进一步地,步骤五的具体操作如下:
以当前飞行时刻tri和飞行攻角αri插值校正网络矩阵Ejz,得到当前状态下的校正网络参数Gjz_ri用于火箭的姿态控制。
综上所述,本申请包括以下至少一种有益技术效果:
在飞行环境不确定性较大的情况下,能够实时辨识火箭飞行过程中的攻角变化,根据攻角大小进行校正网络调整,提高姿态控制系统的鲁棒性以及控制性能,姿态控制精度也相对较好,对误差的适应能力较好。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例的基于攻角变化的姿态控制系统设计方法的流程示意图。
具体实施方式
为了使本申请所要解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
需要说明的是,当元件被称为“固定于”或“设置于”另一个元件,它可以直接在另一个元件上或者间接在该另一个元件上。当一个元件被称为是“连接于”另一个元件,它可以是直接连接到另一个元件或间接连接至该另一个元件上。
需要理解的是,术语“长度”、“宽度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请的限制。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本申请的描述中,“多个”的含义是两个或两个以上,除非另有明确具体的限定。
以下结合附图1对本申请作进一步详细说明。
本申请实施例公开一种基于攻角变化的姿态控制系统设计方法。参照图1,基于攻角变化的姿态控制系统设计方法包括以下步骤:
步骤一,计算飞行攻角包络,具体操作如下:根据总体提供的风场数据使用要求,计算火箭飞行过程中的攻角包络,对于任一特征时刻点ti,存在攻角变化范围[0,αimax],其中i取值为1到n,其中n表示特征时刻点个数,n的取值为正整数,特征时刻点ti的选取特征依据包括但不限于:火箭的飞行弹道和气动特性。
步骤二,设计各时刻点下各攻角状态的校正网络,具体操作如下:对特征时刻点ti,针对攻角变化范围在[0,αimax]内的每一个状态各设计一组校正网络参数Gij,其中j取值为1到m,m的取值为正整数,Gij的下角标i表示第i个特征时刻点,Gij的下角标j表示第j个攻角状态点,对每一个特征时刻点,攻角变化范围不同,选取的攻角状态点个数应该相同,即对每个时间点m为定值。
步骤三,对每一特征时刻点重复步骤二,将各状态下的校正网络组合为参数矩阵装订于飞控软件,具体操作如下:对每一个特征时刻点重复步骤二直至完成该飞行时段内所有特征点校正网络设计,将每一个特征时刻点的每一个攻角状态下的校正网络参数组合为一个mxkxn的矩阵Ejz,将矩阵Ejz装订至飞控软件中,其中k表示任一状态下的校正网络参数的长度,且对每一个状态k值保持不变,飞控软件即火箭的飞行控制软件。
步骤四,火箭飞行过程中实时辨识飞行攻角,具体操作如下:飞行过程中,从火箭启控开始实时辨识飞行攻角,时刻tri的攻角辨识方法如下:通过火箭惯组反馈的箭体过载ny进行计算,具体计算公式如下:
Figure BDA0003392227000000051
式中,m为箭体质量,CN为法向力系数,q为动压,Sm为参考面积,P为摆动发动机推力,
Figure BDA0003392227000000052
为发动机俯仰通道摆角,
Figure BDA0003392227000000053
为俯仰角加速度,xz为质心位置,xa为火箭惯组安装位置;
在箭体过载ny已知的前提下,根据上述式可以求出时刻tri的法向力系数CN,由于CN为随攻角变化的参数,以CN插值法向力系数矩阵,即可得到时刻tri的飞行攻角αri
为避免箭体过载ny变化过快造成飞行攻角αri突变,用于攻角辨识的箭体过载ny选取为h拍内的平均值,h的取值可以通过仿真试验进行确定。
步骤五,利用当前飞行时刻及实时攻角插值校正网络矩阵用于姿态控制,具体操作如下:以当前飞行时刻tri和飞行攻角αri插值校正网络矩阵Ejz,得到当前状态下的校正网络参数Gjz_ri用于火箭的姿态控制。
本申请实施例一种基于攻角变化的姿态控制系统设计方法的实施原理为:在飞行环境不确定性较大的情况下,能够实时辨识火箭飞行过程中的攻角变化,根据攻角大小进行校正网络调整,提高姿态控制系统的鲁棒性以及控制性能,姿态控制精度也相对较好,对误差的适应能力较好。
以上所述仅为本申请的较佳实施例而已,并不用以限制本申请,凡在申请的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本申请的保护范围之内。

Claims (7)

1.一种基于攻角变化的姿态控制系统设计方法,其特征在于,包括以下步骤:
步骤一,计算飞行攻角包络;
步骤二,设计各时刻点下各攻角状态的校正网络;
步骤三,对每一特征时刻点重复步骤二,将各状态下的校正网络组合为参数矩阵装订于飞控软件;
步骤四,火箭飞行过程中实时辨识飞行攻角;
步骤五,利用当前飞行时刻及实时攻角插值校正网络矩阵用于姿态控制。
2.根据权利要求1所述的基于攻角变化的姿态控制系统设计方法,其特征在于,步骤一的具体操作如下:
根据总体提供的风场数据使用要求,计算火箭飞行过程中的攻角包络,对于任一特征时刻点ti,存在攻角变化范围[0,αimax],其中i取值为1到n,其中n表示特征时刻点个数,n的取值为正整数,特征时刻点ti的选取特征依据包括但不限于:火箭的飞行弹道和气动特性。
3.根据权利要求2所述的基于攻角变化的姿态控制系统设计方法,其特征在于,步骤二的具体操作如下:
对特征时刻点ti,针对攻角变化范围在[0,αimax]内的每一个状态各设计一组校正网络参数Gij,其中j取值为1到m,m的取值为正整数,Gij的下角标i表示第i个特征时刻点,Gij的下角标j表示第j个攻角状态点,对每一个特征时刻点,攻角变化范围不同,选取的攻角状态点个数应该相同,即对每个时间点m为定值。
4.根据权利要求3所述的基于攻角变化的姿态控制系统设计方法,其特征在于,步骤三的具体操作如下:
对每一个特征时刻点重复步骤二直至完成该飞行时段内所有特征点校正网络设计,将每一个特征时刻点的每一个攻角状态下的校正网络参数组合为一个m×k×n的矩阵Ejz,将矩阵Ejz装订至飞控软件中,其中k表示任一状态下的校正网络参数的长度,且对每一个状态k值保持不变。
5.根据权利要求4所述的基于攻角变化的姿态控制系统设计方法,其特征在于,步骤四的具体操作如下:
飞行过程中,从火箭启控开始实时辨识飞行攻角,时刻tri的攻角辨识方法如下:通过火箭惯组反馈的箭体过载ny进行计算,具体计算公式如下:
Figure FDA0003392226990000021
式中,m为箭体质量,CN为法向力系数,q为动压,Sm为参考面积,P为摆动发动机推力,
Figure FDA0003392226990000022
为发动机俯仰通道摆角,
Figure FDA0003392226990000023
为俯仰角加速度,xz为质心位置,xa为火箭惯组安装位置;
在箭体过载ny已知的前提下,根据上述式可以求出时刻tri的法向力系数CN,由于CN为随攻角变化的参数,以CN插值法向力系数矩阵,即可得到时刻tri的飞行攻角αri
6.根据权利要求5所述的基于攻角变化的姿态控制系统设计方法,其特征在于:用于攻角辨识的箭体过载ny选取为h拍内的平均值,h的取值可以通过仿真试验进行确定。
7.根据权利要求5所述的基于攻角变化的姿态控制系统设计方法,其特征在于,步骤五的具体操作如下:
以当前飞行时刻tri和飞行攻角αri插值校正网络矩阵Ejz,得到当前状态下的校正网络参数Gjz_ri用于火箭的姿态控制。
CN202111467689.5A 2021-12-03 2021-12-03 基于攻角变化的姿态控制系统设计方法 Active CN114234737B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111467689.5A CN114234737B (zh) 2021-12-03 2021-12-03 基于攻角变化的姿态控制系统设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111467689.5A CN114234737B (zh) 2021-12-03 2021-12-03 基于攻角变化的姿态控制系统设计方法

Publications (2)

Publication Number Publication Date
CN114234737A true CN114234737A (zh) 2022-03-25
CN114234737B CN114234737B (zh) 2024-05-28

Family

ID=80753076

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111467689.5A Active CN114234737B (zh) 2021-12-03 2021-12-03 基于攻角变化的姿态控制系统设计方法

Country Status (1)

Country Link
CN (1) CN114234737B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115562314A (zh) * 2022-10-19 2023-01-03 航天科工火箭技术有限公司 运载火箭子级落区控制方法、系统、介质及计算机设备

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2206044C2 (ru) * 2001-08-20 2003-06-10 Григорьев Владимир Григорьевич Способ определения границ области возможного пуска ракеты
RU2216707C1 (ru) * 2002-04-19 2003-11-20 Федеральное государственное унитарное предприятие "Конструкторское бюро машиностроения" Вращающаяся самонаводящаяся ракета
CN102139766A (zh) * 2010-10-28 2011-08-03 中国科学院力学研究所 一种亚轨道飞行器再入飞行的设计攻角获取方法
CN102139768A (zh) * 2010-10-28 2011-08-03 中国科学院力学研究所 一种亚轨道飞行器再入飞行的攻角制导方法
US20140379176A1 (en) * 2012-06-22 2014-12-25 Isaac M. Ross Method and apparatus for spacecraft attitude control using polynomial interpolation
CN105136422A (zh) * 2015-09-10 2015-12-09 中国航天空气动力技术研究院 风洞试验中修正飞行器模型侧滑弹性角的方法
CN108983800A (zh) * 2018-07-27 2018-12-11 西北工业大学 一种基于深度学习的飞机姿态控制方法
CN110377045A (zh) * 2019-08-22 2019-10-25 北京航空航天大学 一种基于抗干扰技术的飞行器全剖面控制方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2206044C2 (ru) * 2001-08-20 2003-06-10 Григорьев Владимир Григорьевич Способ определения границ области возможного пуска ракеты
RU2216707C1 (ru) * 2002-04-19 2003-11-20 Федеральное государственное унитарное предприятие "Конструкторское бюро машиностроения" Вращающаяся самонаводящаяся ракета
CN102139766A (zh) * 2010-10-28 2011-08-03 中国科学院力学研究所 一种亚轨道飞行器再入飞行的设计攻角获取方法
CN102139768A (zh) * 2010-10-28 2011-08-03 中国科学院力学研究所 一种亚轨道飞行器再入飞行的攻角制导方法
US20140379176A1 (en) * 2012-06-22 2014-12-25 Isaac M. Ross Method and apparatus for spacecraft attitude control using polynomial interpolation
CN105136422A (zh) * 2015-09-10 2015-12-09 中国航天空气动力技术研究院 风洞试验中修正飞行器模型侧滑弹性角的方法
CN108983800A (zh) * 2018-07-27 2018-12-11 西北工业大学 一种基于深度学习的飞机姿态控制方法
CN110377045A (zh) * 2019-08-22 2019-10-25 北京航空航天大学 一种基于抗干扰技术的飞行器全剖面控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
宋征宇;: "运载火箭飞行减载控制技术", 航天控制, vol. 31, no. 05, 31 October 2013 (2013-10-31), pages 3 - 7 *
董朝阳;刘扬;王青;: "带攻角约束的高超声速飞行器自适应反步控制器设计", 宇航学报, vol. 41, no. 02, pages 174 *
贺风华;马克茂;姚郁;: "基于输出预测的姿控发动机控制律优化设计", 航空学报, vol. 30, no. 06, pages 1131 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115562314A (zh) * 2022-10-19 2023-01-03 航天科工火箭技术有限公司 运载火箭子级落区控制方法、系统、介质及计算机设备
CN115562314B (zh) * 2022-10-19 2024-06-07 航天科工火箭技术有限公司 运载火箭子级落区控制方法、系统、介质及计算机设备

Also Published As

Publication number Publication date
CN114234737B (zh) 2024-05-28

Similar Documents

Publication Publication Date Title
US20100211294A1 (en) Control device for internal combustion engine
CN105136422A (zh) 风洞试验中修正飞行器模型侧滑弹性角的方法
CN114234737A (zh) 基于攻角变化的姿态控制系统设计方法
CN111220347A (zh) 一种飞行器气动协调修正方法
CN109540459B (zh) 一种气动特性数值计算结果修正方法
CN110702106B (zh) 一种无人机及其航向对准方法、装置和存储介质
CN114735239B (zh) 航天器轨道机动控制方法、装置和电子设备
CN110155363A (zh) 基于cfd方法的弹性气动数据精确获取方法
CN114608794A (zh) 一种模型风洞虚拟飞行试验气动力系数测量方法
CN114647994B (zh) 一种爬升性能处理方法
CN109631939B (zh) 一种基于磁强计和加速度计的快速对准方法
CN114611420A (zh) 非定常气动力计算精度评估及修正方法
CN116820101A (zh) 一种距离信息缺失下的欠驱动无人艇编队控制方法
CN111460614A (zh) 一种地月l2点转移轨道中途修正方法
CN112034879B (zh) 一种基于高度-射程比的标准轨迹跟踪制导方法
CN111397579B (zh) 一种分段修正的气压高度估算方法
CN113238572A (zh) 基于预设性能控制的预设时间四旋翼无人机姿态跟踪方法
CN110069842B (zh) 舵效快速估算方法
CN110411707A (zh) 串联飞行器级间分离气动特性干扰量预测方法
CN114675054B (zh) 基于风力发电机组塔基载荷的风向识别方法与系统
CN115327906A (zh) 一种四旋翼无人机故障容错控制器的设计方法及系统
CN114234736A (zh) 固体运载火箭转移轨道控制方法
CN113108779B (zh) 一种基于1553b总线三冗余的自主箭测系统
CN108873862A (zh) 一种针对飞行器的控制系统稳定性的综合评估方法
CN114519232B (zh) 一种箭体运动方程系数计算方法及系统

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