[发明专利]一种轴对称无翼无舵导弹主动段零升阻力系数辨识方法

权利要求书

1.一种轴对称无翼无舵导弹主动段零升阻力系数辨识方法，其特征在于，其步骤为：

第一步、预处理实测数据，对试验数据进行各种预处理，数据预处理包括采样时刻校正、异常值剔除；

第二步、预处理理论数据，针对每次飞行试验获取质量、室压、推力数据，并对上述数据进行插值加密处理；

第三步、获取推导参数，按照海拔高度获取动压、马赫数推导参数；

第四步、导弹飞行推力修正，给出导弹主动段飞行过程中发动机推力修正值；

第五步、数据截取，按照如下两个准则截取数据：

1)截取主动段数据，以发动机点火时刻开始到轴向视加速度出现负值结束作为主动段，从实测数据中截取该时间段加速度；

2)截取小攻角飞行段数据，截取弹体攻角α_k≤2°时间段的数据；

第六步、计算零升阻力系数，选取导弹底部面积作为参考面积，重力加速度g，计算零升阻力系数；

第七步、零升阻力系数平滑与数据表构建；

第八步、多次飞行试验辨识结果融合；

第九步、获取其他数据源；

获取地面计算或风洞试验给出的主动段宽域Ma-H-C_A数据表；采用同时模拟外流与底部喷流的计算流体力学方法获取由不同攻角、不同高度、不同马赫数构成的宽域工况下的主动段轴向力系数表Ma-h-α_k-C_A；

第十步、与其他数据源融合，形成易用的气动参数模型：

在Ma-h-α_k-C_A数据表插值获取与数据表中相对应的轴向力系数数据C'_A0，按照C'_A0与差值修正Ma-h-α_k-C_A数据表中对应马赫数下其他攻角和高度轴向力系数数据，以修正后的宽域Ma-h-α_k-C_A数据表作为与现有气动模型相容的轴向力气动参数模型提供总体、控制使用。

2.根据权利要求1所述的方法，其特征在于，所述第一步、预处理实测数据包括：

1)轴向视加速度A_x，单位m/s²；

2)飞行海拔高度h，单位取km，如果没有该数据，需由遥测数据中导航基准系测量的导弹位置x、y、z并结合发射点海拔高度h₀、地球半径R按照公式(1)计算给出，

3)飞行速率v，单位取m/s。如果没有该数据，可由导航基准系导弹速度v的分量v_x、v_y、v_z按照公式(2)计算给出，

4)弹体攻角α_k，单位取“°”，由导航基准系导弹速度分量按照公式(3)计算给出，

5)摆动喷管θ，单位取“°”，由喷管Y、Z两个方向的摆角θ_y、θ_z按照公式(4)计算给出，

6)发动机贮室压强P_zf，单位取MPa。

3.根据权利要求1所述的方法，其特征在于，所述第二步、预处理理论数据包括：

1)导弹飞行理论质量m，对理论质量进行线性插值加密获得导弹飞行过程每一采样时刻的导弹质量；

2)导弹飞行理论室压P_zt，对燃烧室压力理论值进行线性插值加密获得导弹飞行过程每一采样时刻的室压；

3)导弹飞行理论推力F_t，对理论推力进行线性插值加密获得导弹飞行过程每一采样时刻的导弹推力。

4.根据权利要求1所述的方法，其特征在于，所述第三步、获取推导参数包括：

1)导弹主动段飞行动压Q，所述动压用于计算气动系数，标准大气参数表中大气密度与海拔高度对应，按照第一步中的高度h查表获取当地密度ρ，按照公式(5)计算给出动压，

2)飞行马赫数Ma，气动模型按照飞行马赫数给出，由高度h在大气参数表中查询获取当地音速a后按照公式(6)计算给出，