CN114234737B - 基于攻角变化的姿态控制系统设计方法 - Google Patents
基于攻角变化的姿态控制系统设计方法 Download PDFInfo
- Publication number
- CN114234737B CN114234737B CN202111467689.5A CN202111467689A CN114234737B CN 114234737 B CN114234737 B CN 114234737B CN 202111467689 A CN202111467689 A CN 202111467689A CN 114234737 B CN114234737 B CN 114234737B
- Authority
- CN
- China
- Prior art keywords
- attack angle
- flight
- rocket
- characteristic
- state
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000013461 design Methods 0.000 title claims abstract description 17
- 238000012937 correction Methods 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 24
- 230000001133 acceleration Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F42—AMMUNITION; BLASTING
- F42B—EXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
- F42B15/00—Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
- F42B15/01—Arrangements 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进行计算，具体计算公式如下：
式中，m为箭体质量，CN为法向力系数，q为动压，Sm为参考面积，P为摆动发动机推力，为发动机俯仰通道摆角，/>为俯仰角加速度，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进行计算，具体计算公式如下：
式中，m为箭体质量，CN为法向力系数，q为动压，Sm为参考面积，P为摆动发动机推力，为发动机俯仰通道摆角，/>为俯仰角加速度，xz为质心位置，xa为火箭惯组安装位置；
在箭体过载ny已知的前提下，根据上述式可以求出时刻tri的法向力系数CN，由于CN为随攻角变化的参数，以CN插值法向力系数矩阵，即可得到时刻tri的飞行攻角αri。
为避免箭体过载ny变化过快造成飞行攻角αri突变，用于攻角辨识的箭体过载ny选取为h拍内的平均值，h的取值可以通过仿真试验进行确定。
步骤五，利用当前飞行时刻及实时攻角插值校正网络矩阵用于姿态控制，具体操作如下：以当前飞行时刻tri和飞行攻角αri插值校正网络矩阵Ejz，得到当前状态下的校正网络参数Gjz_ri用于火箭的姿态控制。
本申请实施例一种基于攻角变化的姿态控制系统设计方法的实施原理为：在飞行环境不确定性较大的情况下，能够实时辨识火箭飞行过程中的攻角变化，根据攻角大小进行校正网络调整，提高姿态控制系统的鲁棒性以及控制性能，姿态控制精度也相对较好，对误差的适应能力较好。
以上所述仅为本申请的较佳实施例而已，并不用以限制本申请，凡在申请的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本申请的保护范围之内。
Claims (2)
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进行计算，具体计算公式如下：
式中，m为箭体质量，CN为法向力系数，q为动压，Sm为参考面积，P为摆动发动机推力，为发动机俯仰通道摆角，/>为俯仰角加速度，xz为质心位置，xa为火箭惯组安装位置；
在箭体过载ny已知的前提下，根据上述式可以求出时刻tri的法向力系数CN，由于CN为随攻角变化的参数，以CN插值法向力系数矩阵，即可得到时刻tri的飞行攻角αri；
步骤五，利用当前飞行时刻及实时攻角插值校正网络矩阵用于姿态控制，具体操作如下：以当前飞行时刻tri和飞行攻角αri插值校正网络矩阵Ejz，得到当前状态下的校正网络参数Gjz_ri用于火箭的姿态控制。
2.根据权利要求1所述的基于攻角变化的姿态控制系统设计方法，其特征在于：用于攻角辨识的箭体过载ny选取为h拍内的平均值，h的取值通过仿真试验进行确定。
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 CN114234737A (zh) | 2022-03-25 |
| CN114234737B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115562314B (zh) * | 2022-10-19 | 2024-06-07 | 航天科工火箭技术有限公司 | 运载火箭子级落区控制方法、系统、介质及计算机设备 |
Citations (7)
| 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 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的攻角制导方法 |
| CN105136422A (zh) * | 2015-09-10 | 2015-12-09 | 中国航天空气动力技术研究院 | 风洞试验中修正飞行器模型侧滑弹性角的方法 |
| CN108983800A (zh) * | 2018-07-27 | 2018-12-11 | 西北工业大学 | 一种基于深度学习的飞机姿态控制方法 |
| CN110377045A (zh) * | 2019-08-22 | 2019-10-25 | 北京航空航天大学 | 一种基于抗干扰技术的飞行器全剖面控制方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9061775B2 (en) * | 2012-06-22 | 2015-06-23 | Isaac M. Ross | Method and apparatus for spacecraft attitude control using polynomial interpolation |
-
2021
- 2021-12-03 CN CN202111467689.5A patent/CN114234737B/zh active Active
Patent Citations (7)
| 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 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的攻角制导方法 |
| 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)
| Title |
|---|
| 基于输出预测的姿控发动机控制律优化设计;贺风华;马克茂;姚郁;;航空学报;第30卷(第06期);第1131页-第1137页 * |
| 带攻角约束的高超声速飞行器自适应反步控制器设计;董朝阳;刘扬;王青;;宇航学报;第41卷(第02期);第174页-第181页 * |
| 运载火箭飞行减载控制技术;宋征宇;;航天控制;20131031;第31卷(第05期);第3-7、18页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114234737A (zh) | 2022-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114234737B (zh) | 基于攻角变化的姿态控制系统设计方法 | |
| CN111220347A (zh) | 一种飞行器气动协调修正方法 | |
| CN109726437B (zh) | 一种舱门气动载荷等效节点力处理方法 | |
| CN114647994B (zh) | 一种爬升性能处理方法 | |
| CN109490855A (zh) | 一种车载雷达标定方法、装置及车辆 | |
| CN112683446B (zh) | 一种飞机实时重心位置估计方法 | |
| CN109631939B (zh) | 一种基于磁强计和加速度计的快速对准方法 | |
| CN109446582A (zh) | 一种考虑地球自转的高精度降阶平稳滑翔动力学建模方法 | |
| CN113392599A (zh) | 一种弹性飞行器动响应的确定方法 | |
| CN103235515B (zh) | 一种利用零运动避免单框架控制力矩陀螺群框架轴转速死区的方法 | |
| CN111460614B (zh) | 一种地月l2点转移轨道中途修正方法 | |
| CN114964259A (zh) | 基于红外焦平面阵列时域信号与红外图像的姿态解算方法及系统 | |
| CN114611420A (zh) | 非定常气动力计算精度评估及修正方法 | |
| CN109492234A (zh) | 一种改进的流固耦合插值方法 | |
| CN116923730B (zh) | 一种具有自调节预设性能约束的航天器姿态主动容错控制方法 | |
| CN111439392A (zh) | 一种航天器编队位置协同控制方法 | |
| CN107962567B (zh) | 浮动柔性机械臂的分布式振动控制系统及其振动控制方法 | |
| CN117171894A (zh) | 一种考虑静稳定裕度约束的飞行器布局气动优化设计方法 | |
| CN110147571B (zh) | 一种组件结构的拓扑优化方法及装置 | |
| CN109443333B (zh) | 一种陀螺阵列反馈加权融合方法 | |
| CN112034879A (zh) | 一种基于高度-射程比的标准轨迹跟踪制导方法 | |
| CN110411707A (zh) | 串联飞行器级间分离气动特性干扰量预测方法 | |
| CN107101693A (zh) | 非稳态液面多传感器加权系数动态分配方法 | |
| CN115114864A (zh) | 一种基于cfd的飞行器全包线气动数据库生成方法 | |
| CN115979567A (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 |