菱形前体飞行器无驻点压力FADS技术研究

张宗源, 顾蕴松, 李琳恺, 赵冬凯

张宗源, 顾蕴松, 李琳恺, 等. 菱形前体飞行器无驻点压力FADS技术研究[J]. 实验流体力学, doi: 10.11729/syltlx20230125.
引用本文: 张宗源, 顾蕴松, 李琳恺, 等. 菱形前体飞行器无驻点压力FADS技术研究[J]. 实验流体力学, doi: 10.11729/syltlx20230125.
ZHANG Z Y, GU Y S, LI L K, et al. Research on FADS technology of diamond-nosed aircraft without stagnation pressure[J]. Journal of Experiments in Fluid Mechanics, doi: 10.11729/syltlx20230125.
Citation: ZHANG Z Y, GU Y S, LI L K, et al. Research on FADS technology of diamond-nosed aircraft without stagnation pressure[J]. Journal of Experiments in Fluid Mechanics, doi: 10.11729/syltlx20230125.

菱形前体飞行器无驻点压力FADS技术研究

基金项目: 中央高校基本科研业务费专项经费资助项目(NS2022013);江苏高校优势学科建设工程资助项目
详细信息
    作者简介:

    张宗源: (1999—),男,江苏盐城人,硕士研究生。研究方向:FADS技术。通信地址:江苏省南京市秦淮区御道街29号南京航空航天大学明故宫校区飞行测控创新实验室(210016)。E-mail:zzy083117@163.com

    通讯作者:

    顾蕴松: E-mail:yunsonggu@nuaa.edu.cn

  • 中图分类号: V448

Research on FADS technology of diamond-nosed aircraft without stagnation pressure

  • 摘要: 现代战斗机受头部火控雷达等设备限制,无法在机头驻点附近开设测压孔,缺乏驻点压力会导致常规嵌入式大气数据传感系统测试精度大大下降。针对典型战斗机常用的菱形前体外形,对无驻点压力FADS(Flush Air Data Sensing, FADS)系统解算算法及精度进行研究。通过亚跨声速风洞校准试验,获得了测压孔压力分布特性,基于卡尔曼滤波算法构建了无驻点压力FADS技术。通过引入差压数据改进了算法,改进算法实现了部分解耦,提高了解算精度且迭代次数较少,计算量较小。风洞试验结果表明:无驻点压力解算算法可以在外插车次下较好地解算大气参数,其中迎角测量精度为0.33°,侧滑角测量精度为0.30°,静压测量精度为0.67%,马赫数测量精度为0.011。
    Abstract: Limited by equipment such as fire control radars near the nose, modern fighters cannot set pressure ports near the stagnation point, which causes the test accuracy of the conventional flush air data sensing system to greatly decline. For the diamond nose used in the typical fighter, the algorithm and accuracy of FADS (Flush Air Data Sensing) without stagnation pressure are studied. Through subsonic and transonic wind tunnel calibration experiments, the pressure distribution characteristics of the pressure ports were obtained and the FADS technology without stagnation pressure was constructed based on the Kalman filtering algorithm. The algorithm has been improved by importing differential pressure. The improved algorithm is partly decoupled, resulting in improved accuracy with few iterations and low computational complexity. The experimental results show that the algorithm without stagnation pressure can effectively work out air data in external experiment, with an accuracy of 0.33° for angle of attack, 0.30° for angle of sideslip, 0.67% for static pressure, and 0.011 for Mach number.
  • 为精确控制飞行器姿态和轨迹,需时刻掌握飞行器周围的大气数据(迎角、侧滑角、马赫数、静压、动压等)。通过惯性导航系统间接获得的大气参数精度难以满足要求,传统探针式大气数据测量系统则会给飞行器带来隐身性能下降、气动热等问题[1],两者都难以满足现代飞行器的测量需求。高精度、高可靠性、低成本的嵌入式大气数据传感(Flush Air Data Sensing, FADS)系统是新型战斗机及高超声速飞行器大气数据测量的优选。

    与探针式大气数据传感系统不同,FADS系统是在飞行器表面布置测压孔,通过内嵌压力传感器阵列测得飞行器表面压力分布,利用表面压力分布与飞行状态之间的对应关系,采用解算算法获取大气参数[2]。FADS系统可实时测量飞行器大气参数,为大气层内飞行器的飞行控制提供实时高精度飞行参数,目前已应用于多种飞行器,如F–18战斗机[3],X–31[4]、X–33[5]、X–43A[6]等技术验证机,B–2、X–47B等飞翼飞行器[7]及多种航天飞行器[8-9]。国内FADS技术研究起步较晚,飞行试验公开报道较少,目前多为理论研究和风洞试验研究(包括解算算法、补偿算法[10]、故障识别方法[11]、组合导航技术[12]等研究),也有部分团队完成了工程样机研制[13]

    常见的FADS技术一般需在飞行器头部驻点开设测压孔测量总压,但是受限于飞行器头部外形结构和内部空间,以及火控雷达等设备的影响,工程上难以在飞行器头部配置测压孔[14]。部分尖前体飞行器前缘半径过小,也无法在头部驻点配置测压孔[15]。驻点压力对于马赫数解算至关重要[16],缺乏驻点压力信息会导致测压点接收的动压大幅减小,进而导致大气数据解算精度大幅下降,因此,亟需建立高精度、无驻点压力的FADS技术,以满足工程需求。

    Whitmore等[17]提出将FADS系统安装于机翼前缘(LE-FADS),在算法中使用加权最小二乘法。低速风洞试验结果表明该方案具有一定可行性,但受风洞堵塞效应影响,大迎角下的测量值与理论值相差较大。王岩和郑伟[14]针对某型战斗机研究了基于查表法的分布式FADS系统(DFADS)。该系统测压点全部位于雷达罩后方,避免了对雷达系统的干扰,但系统精度尚未得到验证。王鹏等[18]针对某尖楔前体飞行器研究了驻点压力对BP神经网络算法精度的影响,在驻点压力不作为训练输入时,迎角测试误差从0.15°升高至0.25°。王禹等[7]针对飞翼飞行器设计了基于机翼上下或左右压力系数差的解算算法,仿真结果表明该算法精度较高且稳定收敛。

    本文针对菱形前体外形开展亚跨声速风洞试验,对基于卡尔曼滤波算法的无驻点压力FADS技术进行研究及精度分析。

    本文工作基于现代战斗机常用的菱形机头模型。模型长L = 280 mm,宽b = 88 mm,高H = 89 mm。沿机体轴向在5个横截面上开设14个测压孔(以P1~P14表示,孔径均为1.2 mm),其中11个测压孔位于机体下表面,呈十字形分布,3个测压孔位于上表面,如图1所示。将机头驻点定义为坐标系原点,x轴沿机体轴向向后为正。模型支杆设置了基准块,利用角度象限仪可将模型迎角及滚转角调整至0°。迎角机构底座设有弯槽,可供支杆安装的侧滑角至多可达到 ± 6°,模型侧滑角依靠支杆安装角度确定,安装误差≤0.02°。迎角机构控制误差≤0.06°。

    图  1  模型及测压孔
    Fig.  1  Model and FADS points
    表  1  测压孔位置信息
    Table  1  Location information of pressure ports
    测压孔编号 轴向位置/mm 圆周角/(°)
    1 35.0 0
    2 52.5 0
    3 70.0 0
    4 87.5 0
    5 105.0 0
    6 70.0 −30
    7 70.0 −60
    8 70.0 −90
    9 70.0 30
    10 70.0 60
    11 70.0 90
    12 35.0 180
    13 70.0 180
    14 105.0 180
    下载: 导出CSV 
    | 显示表格

    试验在南京航空航天大学空气动力学系NH–1亚跨超声速风洞进行。NH–1风洞是一座间歇型直流下吹式亚跨超三声速风洞,运行马赫数范围为0.2~3.5,试验段口径为0.6 m × 0.6 m。经过标校,亚跨声速下NH–1风洞试验段核心流区气流偏角<0.13°,马赫数误差<0.0037。在本文试验中,校准马赫数为0.3、0.6、0.9、1.05,模拟高度分别为300 m、1.5 km、3.0 km、3.8 km;迎角范围为−10°~10°;侧滑角范围为−6°~3°。此外,还开展了马赫数0.8、侧滑角6°的试验,试验结果用于检验算法模型精度。模型最大横截面积为5422 mm2,风洞试验最大阻塞度约为4.6%。

    图  2  风洞试验布局
    Fig.  2  Experimental Setup

    表面压力测量系统为PSI公司的DTC Initium电子压力扫描阀系统。该系统由主机及多个ESP–64HD模块组成,主机可同时连接多个模块同步采集,每个模块可提供64个通道压力数据,且有多种量程可供选择。模块使用了差压传感器,64个通道共用1个参考端,实验中,参考端一般连接常温静止大气。测量精度为 ± 0.05% FS,试验过程中采样频率为333 Hz,量程为15 psi(约103.425 kPa)。

    目前,FADS技术已发展出较多解算算法,例如三点法(应用案例:X–33[5])、最小二乘法(应用案例:F–18[19])、神经网络法[20]、卡尔曼滤波算法(应用案例:HFLEX[21]、“火星科学实验室MSL”[22])等。三点法、最小二乘法基于钝头体空气动力学模型推导而来,一般认为不适用于尖头体外形。神经网络法拟合精度通常较高,但需要大量训练样本数据。卡尔曼滤波算法已经在HFLEX、“火星科学实验室”中得到了应用,测量精度较好(测压孔在非机头锥布局情况下,该方法也能够保持较好的测量精度[23])。因此,本文基于该算法建立亚跨声速小迎角稳态飞行FADS解算算法。

    大气参数解算是复杂的非线性问题,应用卡尔曼滤波算法需利用泰勒展开将系统线性化。飞行器表面某点压力可以下式表示:

    $$ p_{i} = p_{∞} + qC_{pi} $$ (1)

    式中:pi为第i个测压点的压力值;p为来流静压,q为来流动压;Cpi为与马赫数Ma、迎角α、侧滑角β有关的压力系数函数[14],该函数是FADS算法的核心。pq和来流总压p0、马赫数Ma中只有2个独立变量,因此决定飞行器表面压力的独立变量仅有4个。本文选取$X $= [α, β, p, q]T作为独立变量预测压力。通过风洞试验获得各测压点压力与大气参数变量的关系,即:

    $$ {p_i} = {f_i}\left( {\boldsymbol{X}} \right) + {e_i} $$ (2)

    式中,ei为测量误差。任取n个测压点用于解算,有:

    $$ \left[ {\begin{array}{*{20}{c}} {{p_1}} \\ {{p_2}} \\ \vdots \\ {{p_n}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{f_1}(X)} \\ {{f_2}(X)} \\ \vdots \\ {{f_n}(X)} \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{e_1}} \\ {{e_2}} \\ \vdots \\ {{e_n}} \end{array}} \right] $$ (3)

    将上式简写为向量形式:

    $$ {\boldsymbol{P}} = {\boldsymbol{F}}\left( {\boldsymbol{X}} \right) + {\boldsymbol{e}} $$ (4)

    式(4)为非线性方程,通过级数展开将其线性化,可得压力近似值:

    $$ {\boldsymbol{P}} = {\boldsymbol{F}}\left( {{{\boldsymbol{X}}_0}} \right) + {\left[ {\frac{{\partial {\boldsymbol{F}}}}{{\partial {\boldsymbol{X}}}}} \right]_{{\boldsymbol{X}} = {{\boldsymbol{X}}_0}}}\Delta {\boldsymbol{X}} + {\boldsymbol{e}} $$ (5)

    式中,${\boldsymbol{X}}_0$为先验状态估计值,Δ${\boldsymbol{X}} $为递归过程中的更新向量,即Δ${\boldsymbol{X}} $= ${\boldsymbol{X}} $ − ${\boldsymbol{X}}_0$。定义${\boldsymbol{H}} $为雅可比矩阵:

    $$ {\boldsymbol{H}} = {\left[ {\frac{{\partial {\boldsymbol{F}}}}{{\partial {\boldsymbol{X}}}}} \right]_{{\boldsymbol{X}} = {{\boldsymbol{X}}_0}}} $$ (6)

    定义${\boldsymbol{y}} $为压力残差:

    $$ {\boldsymbol{y}} = {\boldsymbol{P}} - {\boldsymbol{F}}\left( {{{\boldsymbol{X}}_0}} \right) $$ (7)

    式(5)可简化为:

    $$ {\boldsymbol{y}} = {\boldsymbol{H}}\Delta {\boldsymbol{X}} + {\boldsymbol{e}} $$ (8)

    至此,式(5)已从非线性问题转化为线性问题。

    假设测量误差${\boldsymbol{e}} $均值为0,协方差为${{S}} $,根据Gauss–Markov定理,可得Δ${\boldsymbol{X}} $最佳线性无偏估计为:

    $$ \Delta {\boldsymbol{X}} = {\left( {{{\boldsymbol{H}}^{\mathrm{T}}}{{\boldsymbol{S}}^{ - 1}}{\boldsymbol{H}}} \right)^{ - 1}}{{\boldsymbol{H}}^{\mathrm{T}}}{{\boldsymbol{S}}^{ - 1}}{\boldsymbol{y}} $$ (9)

    假设测量误差不相关,则测量误差协方差矩阵${\boldsymbol{S}} $为:

    $${\boldsymbol{ S }}= E\left[ {{\boldsymbol{e}}{{\boldsymbol{e}}^{\mathrm{T}}}} \right] = \left[ {\begin{array}{*{20}{c}} {\sigma _1^2}&{}&{}&{} \\ {}&{\sigma _2^2}&{}&{} \\ {}&{}& \ddots &{} \\ {}&{}&{}&{\sigma _n^2} \end{array}} \right] $$ (10)

    式中,σi为压力测量误差的标准差。

    解算出Δ${\boldsymbol{X}} $后,需进行递归运算,令:

    $$ {\boldsymbol{X}}_{0} = {\boldsymbol{X}}_{0} + \Delta {\boldsymbol{X }} $$ (11)

    经多次迭代运算,当压力残差${\boldsymbol{y}} $趋于0时,估计值${\boldsymbol{X}}_0 $已趋于真实值${\boldsymbol{X}} $。

    得到动压、静压后,马赫数Ma和总压p0可分别根据以下公式求出:

    $$ q = 0.5\gamma {p_{{\infty}} }M{a^2} $$ (12)
    $$ \frac{{{p_0}}}{{{p_{{\infty}} }}} = {\left( {1 + \frac{{\gamma - 1}}{2}M{a^2}} \right)^{\frac{\gamma }{{\gamma - 1}}}} $$ (13)

    式中,γ为空气比热比,值为1.4。

    卡尔曼滤波算法目前应用较少,其缺点主要在于:1)需进行大量矩阵运算,且需要迭代,计算量较大;2)使用绝压进行计算,为满足飞行器全速域飞行要求,需使用大量程传感器,导致在起飞等低速阶段传感器测量精度不高;3)算法不解耦,当某一变量计算误差较大时,其他变量的计算误差也会同时增大。

    为弥补这些缺点,本文将纵截面或横截面测压点的差压引入预测变量${\boldsymbol{P}} $。此时,预测压力的拟合式从pi = p + qCpi简化为pij = qCpijpij = pi pj, Cpij = Cpi Cpj),计算量减小;两点之间的差压远小于每个点的绝压,因此可以使用小量程差压型传感器,提高测量精度。一般纵截面内差压与侧滑角关联较小,横截面内差压与迎角关联较小,引入这部分差压可将预测变量部分解耦,其形式如下:

    $$ \begin{gathered} \left[ \begin{gathered} {p_i} \\ {p_j} \\ {p_k} \\ \vdots \\ \end{gathered} \right] = \left[ \begin{gathered} {p_{{\infty}} } + q{C_{pi}}(\alpha ,\beta ,Ma) \\ {p_{{\infty}} } + q{C_{pj}}(\alpha ,\beta ,Ma) \\ {p_{{\infty}} } + q{C_{pk}}(\alpha ,\beta ,Ma) \\ \vdots \\ \end{gathered} \right] \\ \Downarrow \\ \left[ \begin{gathered} {p_i} \\ {p_i} - {p_j} \\ {p_i} - {p_k} \\ \vdots \\ \end{gathered} \right] = \left[ \begin{gathered} {p_{{\infty}} } + q{C_{pi}}(\alpha ,\beta ,Ma) \\ q{C_{pij}}(\alpha ,Ma) \\ q{C_{pik}}(\beta ,Ma) \\ \vdots \\ \end{gathered} \right] \\ \end{gathered} $$ (14)

    通过风洞试验建立了各测压点的压力数据库,对各测压点的试验数据进行分析,图3展示了Ma = 0.9、β = 0时部分测压点的压力值(图中pi表示第i个测压点的压力值)。可以看出:纵截面上P1、P3、P5、P12、P13和P14的压力随迎角变化显著,即关于迎角敏感;下表面P1、P3和P5的压力随迎角增大单调上升,上表面P12、P13和P14的压力随迎角增大单调下降,符合压力变化规律;P1、P3、P12和P14的压力变化幅度大、线性度好,更有利于解算。侧滑角为0时,对称位置测压点的压力值应相等:P6与P9、P7与P10呈现出较好的对称性,更有利于解算;P8与P11呈现出一定非对称性,这可能是由模型加工误差所导致。解算4个大气参数至少需4个测压孔,增加测压孔数量可以提升精度,但数量过多会导致计算量大幅增加且易引入故障,本文使用了6个测压孔,初步选用P3、P6、P7、P10、P12、P14展示解算过程。

    图  3  部分测压点的风洞试验数据
    Fig.  3  Wind tunnel test data of pressure ports

    关于马赫数Ma、迎角α、侧滑角βCp函数是非常复杂的非线性函数,通常难以直接获得。根据泰勒公式,可将复杂函数近似为形式简单的多项式函数。多项式方程计算简单,易于实现,适合工程应用。

    采用多项式形式对Cp函数进行了拟合。图4展示了不同马赫数下P1的压力系数空间曲面拟合结果。从图中可以看出Cp与迎角、侧滑角大致呈线性关系,与马赫数呈非线性关系。由于马赫数仅有4个,因此关于马赫数采用了二次拟合。

    图  4  不同马赫数下Cp1随迎角和侧滑角变化的空间曲面
    Fig.  4  Cp1 with angle of attack and sideslip at different Mach numbers

    首先,仅使用单一马赫数下的试验数据进行拟合,此时拟合的Cp函数仅涉及迎角和侧滑角,暂不涉及马赫数,便于确定算法精度上限。在不融合其他马赫数数据时,Cp函数拟合度极高,R2达到0.97以上。多次试验的结论基本一致,本文仅介绍其中部分数据结果。以马赫数0.6、侧滑角3°工况的试验数据为例,拟合精度如图5所示(本文以均方根误差衡量精度)。由图5可以得出:FADS测量迎角精度0.11°,侧滑角精度0.15°,静压精度0.18%,马赫数精度0.004;迎角最大误差<0.37°,侧滑角最大误差<0.61°,静压最大误差<0.67%,马赫数最大误差<0.011。可以看出,在精确拟合出Cp函数的前提下,无驻点压力解算算法可以较好地还原马赫数。在本试验的校准范围内,初始值选择上一轮迭代收敛值或恒定值,算法均可收敛,前者一般迭代3~5次达到收敛,后者一般迭代7~10次达到收敛。

    图  5  FADS测量值与真实值对比
    Fig.  5  Comparison of FADS measurement with real values

    风洞校准试验仅能获得离散马赫数下的数据,而飞行器实际飞行过程是一个连续调速调姿过程,需基于离散校准数据建立连续算法,使FADS系统在飞行器整个速域内均可解算大气数据。将所有校准马赫数下的试验数据全部用于拟合,得到较宽速域内的连续Cp函数。为兼容多个马赫数下的数据,连续Cp函数的拟合度有所下降。

    以马赫数0.8、侧滑角6°工况的试验数据检验算法模型的解算能力。图6展示了使用连续Cp函数解算马赫数0.8、侧滑角6°工况试验数据的精度:迎角精度0.27°,侧滑角精度0.44°,静压精度1.14%,马赫数精度0.017;迎角最大误差<0.42°,侧滑角最大误差<1.05°,静压最大误差<1.93%,马赫数最大误差<0.041。在非校准车次里,算法解算误差增大,这是因为校准马赫数仅有4个且间隔较大,导致算法存在欠拟合问题,整体拟合度有所下降。此外,侧滑角6°在校准侧滑角范围之外,这证明算法具有一定外插预测能力。

    图  6  FADS测量值与真实值对比
    Fig.  6  Comparison of FADS measurement with real values

    从以上不同工况的试验结果可以看出:无驻点压力解算算法具有较高精度,且压力系数函数拟合度越高,解算精度越高。若要投入实际应用,还需根据实际需求合理规划校准马赫数,以建立较为精确完整的表面压力模型。

    图7展示了马赫数为0.9时的P14与P3压力系数差值。在校准范围内,Cp14Cp3随迎角变化显著,且几乎与侧滑角无关,引入此差压可以部分地解耦,提高整体解算精度。本文预测变量P最终调整如式(15)。

    图  7  Ma = 0.9,Cp14Cp3随侧滑角变化曲线
    Fig.  7  Cp14Cp3 with sideslip angle at Ma = 0.9
    $$ P{\text{ = }}\left[ {\begin{array}{*{20}{c}} {{p_3}} \\ {{p_6}} \\ \begin{gathered} {p_7} \\ {p_{10}} \\ {p_{12}} \\ \end{gathered} \\ {{p_{14}}} \end{array}} \right] \to P{\text{ = }}\left[ {\begin{array}{*{20}{c}} \begin{gathered} {p_{14}} - {p_3} \\ {p_6} \\ {p_7} \\ \end{gathered} \\ {{p_{10}}} \\ {{p_{12}} - {p_3}} \\ {{p_{14}}} \end{array}} \right] $$ (15)

    图8展示了以优化算法解算马赫数0.8、侧滑角6°工况试验数据的精度:迎角精度0.33°,侧滑角精度0.30°,静压精度0.67%,马赫数精度0.011;迎角最大误差<0.55°,侧滑角最大误差<0.78°,静压最大误差<1.70%,马赫数最大误差<0.036。对比原算法的解算结果,侧滑角、静压、马赫数的解算误差均有所减小,迎角解算误差略增大。在非校准马赫数下,算法迭代次数仍然少于10次即可达到收敛。

    图  8  FADS优化算法测量值与真实值对比
    Fig.  8  Comparison of improved FADS measurement with real values

    表2列出了多种测压孔方案下的解算精度提升幅度(以均方根误差的下降幅度衡量,正值代表误差减小,精度提升)。引入差压数据后,大部分参数的解算精度获得提升,少数参数的解算精度下降但降幅较小,迭代所需次数无明显增加。因此可以得出结论:引入差压数据可以减小计算量并提升解算精度。若应用于其他构型,还需根据气动外形特点,尝试不同测压孔方案和差压绝压组合方案,获得最佳解算精度。

    表  2  精度提升幅度
    Table  2  The increase of accuracy
    测压孔编号 迎角精度
    /(°)
    侧滑角精度
    /(°)
    静压精度
    /%
    马赫数
    精度
    3、6、7、10、12、14 −0.07 0.14 0.47 0.006
    1、6、7、10、12、14 0.89 1.02 −0.08 0.059
    1、3、6、7、10、14 0.59 0.14 0.37 0.018
    1、3、6、7、10、12 0.75 0.61 0.05 0.035
    3、7、9、10、12、14 1.13 0.25 −0.02 0.060
    3、6、9、10、12、14 −0.22 0.13 0.37 0.004
    3、6、7、9、12、14 −0.20 0.19 0.25 0.003
    下载: 导出CSV 
    | 显示表格

    以战斗机常用菱形机头为研究对象,开展了亚跨声速风洞试验,对基于卡尔曼滤波算法的无驻点压力FADS技术进行了研究和改进。研究结果表明:

    1)无驻点压力FADS技术可以在缺乏驻点压力信息的条件下解算出大气数据,其中迎角精度为0.33°,侧滑角精度为0.30°,静压精度为0.67%,马赫数精度为0.011。

    2)合适地引入差压可以减小计算量并提升无驻点压力FADS技术的精度。

    3)无驻点压力FADS技术迭代次数少,计算较快,具有一定工程应用前景。

    本文仅研究了小迎角稳态飞行状态下的测量方法,而真实战斗机常需进行大迎角飞行和高动态机动飞行。较大的运动角速率会对基于表面压力的FADS技术应用产生不利影响。未来还需针对大迎角测量方法及运动角速率的影响开展深入研究。

  • 图  4   不同马赫数下Cp1随迎角和侧滑角变化的空间曲面

    Fig.  4   Cp1 with angle of attack and sideslip at different Mach numbers

    图  1   模型及测压孔

    Fig.  1   Model and FADS points

    图  2   风洞试验布局

    Fig.  2   Experimental Setup

    图  3   部分测压点的风洞试验数据

    Fig.  3   Wind tunnel test data of pressure ports

    图  5   FADS测量值与真实值对比

    Fig.  5   Comparison of FADS measurement with real values

    图  6   FADS测量值与真实值对比

    Fig.  6   Comparison of FADS measurement with real values

    图  7   Ma = 0.9,Cp14Cp3随侧滑角变化曲线

    Fig.  7   Cp14Cp3 with sideslip angle at Ma = 0.9

    图  8   FADS优化算法测量值与真实值对比

    Fig.  8   Comparison of improved FADS measurement with real values

    表  1   测压孔位置信息

    Table  1   Location information of pressure ports

    测压孔编号 轴向位置/mm 圆周角/(°)
    1 35.0 0
    2 52.5 0
    3 70.0 0
    4 87.5 0
    5 105.0 0
    6 70.0 −30
    7 70.0 −60
    8 70.0 −90
    9 70.0 30
    10 70.0 60
    11 70.0 90
    12 35.0 180
    13 70.0 180
    14 105.0 180
    下载: 导出CSV

    表  2   精度提升幅度

    Table  2   The increase of accuracy

    测压孔编号 迎角精度
    /(°)
    侧滑角精度
    /(°)
    静压精度
    /%
    马赫数
    精度
    3、6、7、10、12、14 −0.07 0.14 0.47 0.006
    1、6、7、10、12、14 0.89 1.02 −0.08 0.059
    1、3、6、7、10、14 0.59 0.14 0.37 0.018
    1、3、6、7、10、12 0.75 0.61 0.05 0.035
    3、7、9、10、12、14 1.13 0.25 −0.02 0.060
    3、6、9、10、12、14 −0.22 0.13 0.37 0.004
    3、6、7、9、12、14 −0.20 0.19 0.25 0.003
    下载: 导出CSV
  • [1]

    JOST M, SCHWEGMANN F, KOHLER T. Flush air data system – an advanced air data system for the aerospace industry[C]//Proc of the AIAA Guidance, Navigation, and Control Conference and Exhibit. 2004. .

    [2] 丁智坚, 周欢, 吴东升, 等. 嵌入式大气数据测量系统技术研究进展[J]. 宇航学报, 2019, 40(3): 247–257. DOI: 10.3873/j.issn.1000-1328.2019.03.001

    DING Z J, ZHOU H, WU D S, et al. Review of flush air data sensing system[J]. Journal of Astronautics, 2019, 40(3): 247–257. doi: 10.3873/j.issn.1000-1328.2019.03.001

    [3]

    WHITMORE S, MOES T, LARSON T. Preliminary results from a subsonic high-angle-of-attack flush airdata sensing (HI-FADS) system: Design, calibration, algorithm development, and flight test evaluation[C]//Proc of the 28th Aerospace Sciences Meeting. 1990. .

    [4]

    WEISS S. Comparing three algorithms for modeling flush air data systems[C]//Proc of the 40th AIAA Aerospace Sciences Meeting & Exhibit. 2002. .

    [5]

    WHITMORE S, COBLEIGH B, HAERING E Jr. Design and calibration of the X-33 flush airdata sensing (FADS) system[C]//Proc of the 36th AIAA Aerospace Sciences Meeting and Exhibit. 1998. .

    [6]

    BAUMANN E, PAHLE J W, DAVIS M C, et al. X-43A flush airdata sensing system flight-test results[J]. Journal of Spacecraft and Rockets, 2010, 47(1): 48–61. doi: 10.2514/1.41163

    [7] 王禹, 郑伟, 童建忠, 等. 飞翼飞机嵌入式大气数据系统算法研究[J]. 测控技术, 2022, 41(9): 101–106. DOI: 10.19708/j.ckjs.2021.08.263

    WANG Y, ZHENG W, TONG J Z, et al. Flush air data sensing system algorithm of flying wing aircraft[J]. Measurement & Control Technology, 2022, 41(9): 101–106. doi: 10.19708/j.ckjs.2021.08.263

    [8]

    HENRY M, WOLF H, SIEMERS P. An evaluation of Shuttle Entry Air Data System (SEADS) flight pressures - Comparisons with wind tunnel and theoretical predictions[C]//Proc of the 15th Aerodynamic Testing Conference. 1988. .

    [9]

    STEPHAN T, MARKUS S, MICHAEL C, et al. Hybrid navigation system for the SHEFEX II mission[C]//Proc of the AIAA Guidance, Navigation and Control Sensor Technologies in Europe. 2008. .

    [10] 王希洋, 柏楠, 苑景春, 等. FADS系统变管径引气管路压力延迟误差补偿方法[J]. 战术导弹技术, 2015(2): 37–42. DOI: 10.16358/j.issn.1009-1300.2015.02.08

    WANG X Y, BAI N, YUAN J C, et al. Compensation method for pressure delay error of variable diameter exhaust pipe in FADS system[J]. Tactical Missile Technology, 2015(2): 37–42. doi: 10.16358/j.issn.1009-1300.2015.02.08

    [11] 李清东, 陈璐璐, 张孝功, 等. FADS快速智能故障检测和诊断技术[J]. 系统工程与电子技术, 2009, 31(10): 2544–2546. , 2009, 31(10): 2544–2546. doi: 10.3321/j.issn:1001-506X.2009.10.058

    LI Q D, CHEN L L, ZHANG X G, et al. Flush airdata sensing system fast intelligent fault detection and diagnosis technology[J]. Systems Engineering and Electronics, 2009, 31(10): 25442546. DOI: 10.3321/j.issn:1001-506X.2009.10.058

    [12] 杨胜江, 赵景朝, 杨志红. 嵌入式大气数据传感与惯性导航信息融合方法研究[J]. 战术导弹技术, 2016(2): 95–100. DOI: 10.16358/j.issn.1009-1300.2016.02.17

    YANG S J, ZHAO J C, YANG Z H. Research on information fusion method of embedded atmospheric data sensing and inertial navigation[J]. Tactical Missile Technology, 2016(2): 95–100. doi: 10.16358/j.issn.1009-1300.2016.02.17

    [13] 陈广强, 刘吴月, 豆修鑫, 等. 吸气式空空导弹FADS系统设计[J]. 中国科学:技术科学, 2016, 46(11): 1193–1206. DOI: 10.1360/N092016-00258

    CHEN G Q, LIU W Y, DOU X X, et al. Flush air data sensing system design for air breathing air-to-air missile[J]. Scientia Sinica:Technologica, 2016, 46(11): 1193–1206. doi: 10.1360/N092016-00258

    [14] 王岩, 郑伟. 分布嵌入式大气数据系统算法的初步研究[J]. 飞机设计, 2008, 28(6): 5–11,26. DOI: 10.3969/j.issn.1673-4599.2008.06.002

    WANG Y, ZHENG W. Elementary study on the distributed flush air data system arithmetic[J]. Aircraft Design, 2008, 28(6): 5–11,26. doi: 10.3969/j.issn.1673-4599.2008.06.002

    [15] 王鹏, 金鑫. 尖锥前体飞行器FADS系统的人工神经网络建模及风洞试验研究[J]. 实验流体力学, 2019, 33(5): 57–63. DOI: 10.11729/syltlx20180125

    WANG P, JIN X. Study on artificial neural network modeling and wind tunnel test for the FADS system applied to the vehicle with sharp nosed fore-bodies[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(5): 57–63. doi: 10.11729/syltlx20180125

    [16] 陈广强, 刘吴月, 豆修鑫, 等. 超声速飞行器FADS系统实时解算设计与验证[J]. 空气动力学学报, 2018, 36(4): 561–570. DOI: 10.7638/kqdlxxb-2016.0073

    CHEN G Q, LIU W Y, DOU X X, et al. Flush air data sensing system real-time solving design and verification for supersonic vehicle[J]. Acta Aerodynamica Sinica, 2018, 36(4): 561–570. doi: 10.7638/kqdlxxb-2016.0073

    [17]

    WHITMORE S, MOES T, CZERNIEJEWSKI M, et al. Application of a flush airdata sensing system to a wing leading edge(LE-FADS)[C]//Proc of the 31st Aerospace Sciences Meeting. 1993. .

    [18] 王鹏, 胡远思, 金鑫. 尖楔前体飞行器FADS系统驻点压力对神经网络算法精度的影响[J]. 宇航学报, 2016, 37(9): 1072–1079. DOI: 10.3873/j.issn.1000-1328.2016.09.006

    WANG P, HU Y S, JIN X. Effect of stagnation pressure on the neural network algorithm accuracy for FADS system applied to the vehicle with sharp wedged fore-bodies[J]. Journal of Astronautics, 2016, 37(9): 1072–1079. doi: 10.3873/j.issn.1000-1328.2016.09.006

    [19]

    WHITMORE S A, MOES T R, LEONDES C T. Development of a pneumatic high-angle-of-attack flush airdata sensing system[J]. Control and Dynamic Systems, 1992: 453-511. .

    [20]

    ROHLOFF T J, WHITMORE S A, CATTON I. Fault-tolerant neural network algorithm for flush air data sensing[J]. Journal of Aircraft, 1999, 36(3): 541–549. doi: 10.2514/2.2489

    [21]

    TAKAKI R, TAKIZAWA M, et al. ADS measurement of HYFLEX (HYpersonic FLight EXperiment)[C]//Proc of the 35th Aerospace Sciences Meeting and Exhibit. 1997. .

    [22]

    KARLGAARD C D, KUTTY P, SCHOENENBERGER M, et al. Mars entry atmospheric data system trajectory reconstruction algorithms and flight results[C]//Proc of the 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition. 2013. .

    [23] 黄喜元, 陈洪波, 朱如意. 高超声速飞行器嵌入式大气数据获取技术研究[J]. 导弹与航天运载技术, 2017(3): 58–64. DOI: 10.7654/j.issn.1004-7182.20170313

    HUANG X Y, CHEN H B, ZHU R Y. Research on air data acquisition technology of hypersonic vehicles[J]. Missiles and Space Vehicles, 2017(3): 58–64. doi: 10.7654/j.issn.1004-7182.20170313

图(8)  /  表(2)
计量
  • 文章访问数:  146
  • HTML全文浏览量:  61
  • PDF下载量:  10
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-09-26
  • 修回日期:  2023-11-18
  • 录用日期:  2023-11-21
  • 网络出版日期:  2024-01-09

目录

/

返回文章
返回
x 关闭 永久关闭