一种固体捆绑火箭助推器落点计算方法技术

技术编号:39520847 阅读:22 留言:0更新日期:2023-11-25 19:00
本发明专利技术公开了一种固体捆绑火箭助推器落点计算方法,包括建立固体助推器的六自由度动力学模型;六自由度动力学模型为被动段状态变量关于箭体系上的合力矩的函数;建立来流攻角模型和来流方位角模型;根据来流攻角模型和来流方位角模型得到助推器下落过程中箭体系上的合力矩;以助推器分离点的数据作为仿真初始值,将所述箭体系上的合力矩代入六自由度动力学模型进行仿真,得到助推器落点信息

【技术实现步骤摘要】
一种固体捆绑火箭助推器落点计算方法


[0001]本专利技术涉及一种固体捆绑火箭助推器落点计算方法,适用于固体捆绑运载火箭固体助推器落点计算,属于运载火箭弹道设计领域


技术介绍

[0002]固体捆绑运载火箭利用固体助推器的大推重比特点,有效的增加火箭的运载能力,同时由于固体助推器分离后将无控的返回地面,将面临落点安全性要求,固体助推器相比液体助推器在被动段下落过程中有着不一样的规律,固体助推器分离点质量要远大于液体助推器分离点质量,相比于液体助推器质心集中在尾部发动机处,固体助推器质心基本处于助推器中心偏下处,在无控下落过程中,液体助推器在气动力矩的作用下将最终稳定为尾部发动机朝下的状态,而固体助推器由于压心的反复变动,使得固体助推器在下落过程中姿态会一直处于旋转状态,因此需要建立精细化的气动力及气动力矩模型,来分析固体助推器下落过程中的弹道特性,同时由于固体助推器质量较大,落点速度较大,因此对固体助推器的落点计算精度提出了更高的要求,现有的三自由度落点计算方法仅考虑落体的质心运动,不适应固体助推的运动状态,无法满足高精度需求


技术实现思路

[0003]本专利技术的目的在于克服上述缺陷,提供一种固体捆绑火箭助推器落点计算方法,解决了传统三自由度落点计算方法无法满足高精度需求的技术问题,本专利技术可快速分析固体助推器在下落段的各弹道特征参数,发现其被动段运动规律,使固体助推器落点计算结果更加准确

[0004]为实现上述专利技术目的,本专利技术提供如下技术方案:<br/>[0005]本专利技术涉及一种固体捆绑火箭助推器落点计算方法,属于运载火箭弹道设计领域,包括如下步骤:首先建立固体助推器的6自由度动力学模型,定义来流攻角,通过气动仿真建模,生成不同来流攻角

不同马赫数下的助推器被动段气动数据,然后建立助推器下落过程中气动力及气动力矩公式;根据助推器分离点的位置

速度

高度

姿态

质量,进行动力学仿真,计算固体助推器在无控状态下的被动段落点,分析固体助推器在下落过程中的气动力变化曲线及姿态运动曲线,为较好的预报助推器落点提供理论支持

本设计方法对固体助推器下落过程中所受到的气动力及力矩进行精细化建模,可快速分析固体助推器在下落段的各弹道特征参数,发现其被动段运动规律,使固体助推器落点计算结果更加准确

[0006]本专利技术一种固体捆绑火箭助推器落点计算方法,包括:
[0007]建立固体助推器的六自由度动力学模型;六自由度动力学模型为被动段状态变量关于箭体系上的合力矩的函数;
[0008]建立来流攻角模型和来流方位角模型;根据来流攻角模型和来流方位角模型得到助推器下落过程中箭体系上的合力矩;
[0009]获取助推器分离点的数据;
[0010]以助推器分离点的数据作为仿真初始值,将所述箭体系上的合力矩代入六自由度动力学模型进行仿真,得到助推器落点信息

[0011]进一步的,六自由度动力学模型包括:
[0012][0013]其中,
X
t
为被动段状态变量,为固体捆绑火箭助推器的发惯系位置,为固体捆绑火箭助推器的发惯系速度,
[q
0 q
1 q
2 q3]T
为固体捆绑火箭助推器的姿态四元数,
[
ω
bx ω
by ω
bz
]T
为火箭的角速度,
G
为固体捆绑火箭助推器的质量;
[0014][
ω
bx ω
by ω
bz
]T
为箭体系上的合力矩
[M
bx M
by M
bz
]T
的函数

[0015]进一步的,六自由度动力学模型还包括:
[0016][0017]其中
,

[q
0 q
1 q
2 q3]T
的导数,为
[
ω
bx ω
by ω
bz
]T
的导数,为视加速度在箭体系中坐标分量,为地心发射惯性系中的地球引力加速度,为被动段的总秒耗量,
[MI]为惯量矩阵,
[MI]‑1为惯量矩阵的逆矩阵,为箭体系到发惯系的转换矩阵

[0018]进一步的,来流攻角模型为:
[0019][0020]其中,
α
Q
为来流攻角,具体为火箭相对来流的速度矢量与箭体坐标系
X1
轴夹角,
α
Q
∈[0

π
];
|V
Q
|
为火箭相对于来流的速度矢量
V
Q
的模;为火箭相对于来流的速度矢量
V
Q
在箭体坐标系
X1
轴的分量

[0021]进一步的,方位角模型为:
[0022][0023]其中,
Φ
Q
为来流方位角,具体为按箭体坐标系
X1
轴正向右手法则旋转,箭体
I
基准到的夹角,
Φ
Q
∈[0
,2π
)
;为
V
Q
在箭体坐标系
Y1
轴的分量;为
V
Q
在箭体坐标系
Z1
轴的分量,为
V
Q
在箭体坐标系
Y1Z1
平面的分量

[0024]进一步的,根据来流攻角模型和方位角模型建立助推器下落过程中箭体系上的合力矩的方法包括:
[0025]在不同来流攻角

不同马赫数的条件下,通过气动仿真建模的方法得到气动阻力系数
C
x
和法向力系数
C
n
关于来流攻角和马赫数的二维被动段气动插值数表;
[0026]根据当前助推器的状态
X
t
得到
V
Q

[0027]根据
V
Q
和来流攻角模型得到当前的来流攻角
α
Q
和马赫数;基于二维被动段气动插值数表,将当前的来流攻角
α
Q
和马赫数进行二维插值得到当前的气动阻力系数
C
x
和法向力系数
C
n

[0028]根据
V
Q
和方位角模型得到当前的来流方位角
Φ
Q

[0029]根据当前的气动阻力系数
C
x

法向力系数<本文档来自技高网
...

【技术保护点】

【技术特征摘要】
1.
一种固体捆绑火箭助推器落点计算方法,其特征在于,包括:建立固体助推器的六自由度动力学模型;六自由度动力学模型为被动段状态变量关于箭体系上的合力矩的函数;建立来流攻角模型和来流方位角模型;根据来流攻角模型和来流方位角模型得到助推器下落过程中箭体系上的合力矩;获取助推器分离点的数据;以助推器分离点的数据作为仿真初始值,将所述箭体系上的合力矩代入六自由度动力学模型进行仿真,得到助推器落点信息
。2.
根据权利要求1所述的一种固体捆绑火箭助推器落点计算方法,其特征在于,六自由度动力学模型包括:其中,
X
t
为被动段状态变量,为固体捆绑火箭助推器的发惯系位置,为固体捆绑火箭助推器的发惯系速度,
[q
0 q
1 q
2 q3]
T
为固体捆绑火箭助推器的姿态四元数,
[
ω
bx ω
by ω
bz
]
T
为火箭的角速度,
G
为固体捆绑火箭助推器的质量;
[
ω
bx ω
by ω
bz
]
T
为箭体系上的合力矩
[M
bx M
by M
bz
]
T
的函数
。3.
根据权利要求2所述的一种固体捆绑火箭助推器落点计算方法,其特征在于,六自由度动力学模型还包括:
其中
,

[q
0 q
1 q
2 q3]
T
的导数,为
[
ω
bx ω
by ω
bz
]
T
的导数,为视加速度在箭体系中坐标分量,为地心发射惯性系中的地球引力加速度,为被动段的总秒耗量,
[MI]
为惯量矩阵,
[MI]
‑1为惯量矩阵的逆矩阵,为箭体系到发惯系的转换矩阵
。4.
根据权利要求3所述的一种固体捆绑火箭助推器落点计算方法,其特征在于,来流攻角模型为:
其中,
α
Q
为来流攻角,具体为火箭相对来流的速度矢量与箭体坐标系
X1
轴夹角,
α
Q
∈[0

π
]

|V
Q
|
为火箭相对于来流的速度矢量
V
Q
的模;为火箭相对于来流的速度矢量
V
Q
在箭体坐标系
X1
轴的分量
。5.
根据权利要求4所述的一种固体捆绑火箭助推器落点计算方法,其特征在于,方位角模型为:其中,
Φ
Q
为来流方位角,具体为按箭体坐标系
X1
轴正向右手法则旋转,箭体
I
基准到的夹角,
Φ
Q
∈[0
,2π
)
;为
V
Q
在箭体坐标系
Y1
轴的分量;为
V
Q
在箭体坐标系
Z1
轴的...

【专利技术属性】
技术研发人员:陈凯张哲玮冯继航王颖程川唐明亮匡东政潘文俊张学功张卫东洪刚张亮
申请(专利权)人:上海宇航系统工程研究所
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1