CN113281001A - 基于集成式微型大气数据模块的全速域大气数据解算方法 - Google Patents
基于集成式微型大气数据模块的全速域大气数据解算方法 Download PDFInfo
- Publication number
- CN113281001A CN113281001A CN202110408075.3A CN202110408075A CN113281001A CN 113281001 A CN113281001 A CN 113281001A CN 202110408075 A CN202110408075 A CN 202110408075A CN 113281001 A CN113281001 A CN 113281001A
- Authority
- CN
- China
- Prior art keywords
- pressure
- angle
- atmospheric data
- data
- formula
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/08—Aerodynamic models
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/02—Wind tunnels
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了基于集成式微型大气数据模块的全速域大气数据解算方法,涉及嵌入式大气数据测量技术领域,该算法通过线性插值提供了一种将离散的风洞校准数据扩展为连续的数据的方法;以马赫数为判别依据,通过残差收敛的方法提供了一种提高大气参数解算精度的方式。该算法流程简洁有效,计算效率高,解算结果精度较高,适用于低功耗的机载集成式大气数据模块。
Description
技术领域
本发明涉及大气参数测量领域,尤其涉及了基于集成式微型大气数据模块的全速域大气数据解算方法。
背景技术
随着航空航天飞行器不断向高速高空发展,高超声速航空航天飞行器对大气层内飞行时,迎角、侧滑角、空速和动压等大气参数的精准测量提出了更高的要求。传统的航空飞行器所使用的风标式传感器,存在机械部件多、动态特性差的缺点;传统的航天飞行器使用的惯导系统,则存在测量累计误差逐渐增大的缺点。基于嵌入式大气数据传感系统技术(Flush Air Data Sensing System)的测量方式近年来成为国内外大气参数测量研究的热点,其具有不破坏飞行器表面气动外形、受气动热效应影响小、精度高以及响应速度快的特点。
目前,国外搭载FADS技术的验证机已进行飞行验证,同时对嵌入式大气数据传感系统技术的机理研究已具备一定基础,国内对于FADS技术的研究基本停留在理论阶段和原理样机试制阶段。普遍使用“压力传感器——电脑”的搭配方式,由传感器进行数据采集,电脑进行大气数据的解算与传输,该方案体积大、部件多,不适合在飞行器上搭载,解算算法对解算单元的性能要求较高。
目前,国内并未见到可适用于集成化、小型化、运算能力相对有限的集成式大气数据模块的可测量高超声速飞行器全速域大气数据解算方法。
发明内容
本发明提供基于集成式微型大气数据模块的全速域大气数据解算方法,能够结合风洞校准试验结果,为高超声速飞行器提供可行的、连续的全速域数据解算算法,计算效率高,涉及的计算量有限。
基于集成式微型大气数据模块的全速域大气数据解算方法,包括:
S1、将测量压力数据带入无量纲系数CMa的经验公式,求得无量纲系数CMa的值,再根据CMa和飞行器马赫数的关系得到飞行器马赫数,即空速的粗略值Ma0,将空速的粗略值Ma0标记为计算初值Mai-1;
S2、以所述计算初值Mai-1为输入值,根据气动数据库中搭载的风洞试验给定马赫数下的迎角αA、侧滑角βA以及马赫数MaA等大气数据,与气动数据库中风洞试验得到的无量纲系数Cα A、Cβ A间的映射关系,利用线性差值算法以及多项式最小二乘法,拟合得到空速为Mai时,飞行器迎角α、侧滑角β的解算公式;
S3、将所述空速为Mai时,测量得到的压力数据P1、P2、P3、P4、P5根据测压孔的分布规律组合为无量纲系数Cα、Cβ,利用Cα、Cβ解算迎角α、侧滑角β;
S4、利用所述压力数据和所述气动数据库中搭载的风洞校准试验给定的迎角αA、侧滑角βA与风洞校准试验得到的无量纲修正系数Co A和Cq A间的对应关系,利用线性差值算法以及多项式最小二乘法,拟合得到空速为Mai时,无量纲修正系数Co和Cq的解算公式;
S5、利用S3中计算得到的迎角α、侧滑角β,利用S4中公式计算得到修正系数Co和Cq,联合计算得到静压P∞、动压q∞和修正后的空速Mai;
S6、判断Mai和Mai-1间的差值是否满足精度要求,满足则输出迎角α、侧滑角β、动压q∞和空速Mai,否则返回S2,循环执行S2-S6。
进一步的,当所述测量压力数据包含五个飞行器模型头部测压孔压力值时,所述经验公式为:
其中,P1、P2、P3、P4、P5,分别为所述飞行器模型头部测压孔测得的压力值。
进一步的,所述无量纲系数CMa和飞行器马赫数Ma的关系为:
进一步的,所述线性差值算法为:
式中,Aj、Bj、Cj和Dj为计算得到的空间直线系数;
进一步的,当所述测压孔数量为5个,且呈十字形分布时,无量纲系数Cα、Cβ和所述压力数据之间的关系是:
P1、P2、P3、P4、P5,分别为所述飞行器模型头部测压孔测得的压力值。
进一步的,在所述S3中,所述迎角α、侧滑角β、无量纲修正系数Co、Cq的计算方法为:
当马赫数Ma给定时,无量纲系数与大气数据的函数关系可由多项式最小二乘法拟合的曲线公式表征,公式如下:
计算角度时,Ai表征的含义为迎角α、侧滑角β;Bi表征的含义为无量纲系数Cα、Cβ;
修正马赫数时,Ai表征的含义为无量纲系数Co、Cq;
Bi表征的含义为迎角α、侧滑角β;
n为多项式公式的次数。
进一步的,在S5中,所述动压q∞和空速Mai的计算方法为:
SS1、根据所述S4计算得到的无量纲修正系数Co、Cq,计算总压Po和静压P∞,其中总压Po的计算公式为:
静压P∞的计算公式为:
P1、P2、P3、P4、P5,分别为所述飞行器模型头部测压孔测得的压力值,
动压q∞的计算公式为:
q∞=Po-P∞
SS2、空速Mai由上一循环Mai-1得到的计算公式得到:
若Mai-1≤0.3,马赫数的计算服从伯努利方程,空速Mai计算公式如下:
式中,q∞为动压,ρ为当地大气密度,a为当地音速;
若Mai-1>0.3,马赫数的计算服从一维等熵绝热流公式,空速Mai计算公式如下:
式中,γ为气体比热比常数。
进一步的,所述气动数据库中搭载的数据包括:
风洞校准实验时得到的迎角αA、侧滑角βA、动压q∞ A、静压p∞ A以及马赫数MaA,与模型表面压力Pi A特定组合得到的无量纲系数Cα A、Cβ A、Co A和Cq A,式中上标A表示为风洞标定试验结果。
本发明的有益效果是:
(1)本发明结合风洞校准实验特点,使用线性插值方法提供了适用于集成式大气数据模块的解算公式扩展方法;
(2)本发明通过以马赫数为判别参数,通过残差迭代的方法优化了大气数据的解算精度。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是集成式微型大气数据模块的示意图;
图2是实施例中使用的某运载火箭头部大气数据测量方案内部框架图;
图3是实施例中使用的某运载火箭头部测压孔示意图;
图4是线性插值方法的图示;
图5是一种基于集成式大气数据模块的全速域大气数据解算方法的设计流程图。
具体实施方式
为使本领域技术人员更好地理解本发明的技术方案,下面结合具体实施例对本发明作进一步详细描述。
一种基于集成式大气数据模块的全速域大气数据解算方法,基于一种集成式大气数据模块,集成式大气数据模块是一种基于stm32微控制器的一种包括传感器、存储器、stm32微控制器、串口通信端口,如图1、2所示。集成式大气数据模块内嵌于飞行器内部,飞行器表面测压孔和传感器通过引压管路连接。存储器内部存储风洞校准实验结果得到的气动数据库以及解算程序,可以用来解算实时大气数据。通过RS422串口向上位机传输包括迎角、侧滑角、空速和动压等大气数据。一种基于集成式大气数据模块的全速域大气数据解算方法,其使用的气动数据库来源于风洞校准实验,解算算法对数据库的使用有如下特点:风洞试验是若干个离散的、特定的空速下的校准,得到的是特定空速下的大气数据解算公式。为获取其它未经过校准的空速下的大气数据解算公式,通过线性插值的方法得到Ma≤2.0连续的解算公式。
本实施例中,采用模型为1:40的某运载火箭缩比模型,模型头部测压孔布局形式如图3所示,风洞校准实验风速如下表所示:
表1、实施例风洞校准实验指标
(a)解算公式
模型头部测压孔分布如图3所示,压力值分别取为P1、P2、P3、P4及P5。通过如下公式,经风洞实验后,计算各马赫数下的无量纲系数并建立Cα、Cβ、Co和Cq的气动数据库。
无量纲系数CMa只与空速相关,可以近似的计算得到高超声速飞行器的所处马赫数的粗略值。
设定无量纲系数Cα、Cβ如下所示:
式中,无量纲系数Cα、Cβ分别与迎角α、侧滑角β呈对应的函数关系,用于计算迎角α、侧滑角β。
上述公式适用于0.3≤Ma≤2.0时计算马赫数。
(b)线性插值方法
将表1风洞校准实验中给定的离散的马赫数通过线性插值扩展到Ma≤2.0的连续域。其方法如图4所示,(α1,CαCα1,Ma1)和(α2,CαCα2,Ma2)为风洞标定试验中得到的气动数据库中的求解迎角的两组数据,(α3,CαCα3,Ma3)介于(α1,CαCα1,Ma1)和(α2,CαCα2,Ma2)之间,根据Ma3的值使用线性插值方法求出中间点(α3,CαCα3)的值,即:
根据(α1,CαCα1,Ma1)和(α2,CαCα2,Ma2)数据组中的每对数据可求解空间直线公式为:
式中,Aj、Bj、Cj和Dj为计算得到的空间直线系数;
将Ma3代入空间直线公式,即可得到数据组(α3,CαCα3,Ma3)。
(C)算法流程如图5所示:
(1)由CMa计算Ma0;
(2)由Ma0,线性插值得到马赫数为Ma0情况下的无量纲系数Cα、Cβ及其解算公式;
(3)计算迎角α、侧滑角β,修正Co、Cq,并根据Co、Cq计算动压和修正后的马赫数;
(4)以ΔMaΔMa=|MaiMai-Mai-1Mai-1|<δ<δ(δ数值根据精度要求选定)为精度要求,满足时输出迎角、侧滑角、空速和动压等大气数据;若不满足返回步骤2,以新的马赫数插值得到新的解算公式,迭代直至满足马赫数精度要求。
经过实验,在风洞试验中,测量范围内,该方法测量的迎角、侧滑角的精度优于1°,空速精度优于5%,数据输出频率为100hz。
上述算法由基于stm32微控制器的集成式微型大气数据模块运行。目前,精度很高的神经网络算法由于计算量大,难以在算力有限的单片机上运行。本发明实现了连续的大气数据实时计算,且精度较高。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (8)
1.一种基于集成式大气数据模块的全速域大气数据解算方法,其特征在于,包括:
S1、将测量压力数据带入无量纲系数CMa的经验公式,求得无量纲系数CMa的值,再根据CMa和飞行器马赫数的关系得到飞行器马赫数,即空速的粗略值Ma0,将空速的粗略值Ma0标记为计算初值Mai-1;
S2、以所述计算初值Mai-1为输入值,根据气动数据库中搭载的风洞试验给定马赫数下的迎角αA、侧滑角βA以及马赫数MaA等大气数据,与气动数据库中风洞试验得到的无量纲系数Cα A、Cβ A间的映射关系,利用线性差值算法以及多项式最小二乘法,拟合得到空速为Mai时,飞行器迎角α、侧滑角β的解算公式;
S3、将所述空速为Mai时,测量得到的压力数据P1、P2、P3、P4、P5根据测压孔的分布规律组合为无量纲系数Cα、Cβ,利用Cα、Cβ解算迎角α、侧滑角β;
S4、利用所述压力数据和所述气动数据库中搭载的风洞校准试验给定的迎角αA、侧滑角βA与风洞校准试验得到的无量纲修正系数Co A和Cq A间的对应关系,利用线性差值算法以及多项式最小二乘法,拟合得到空速为Mai时,无量纲修正系数Co和Cq的解算公式;
S5、利用S3中计算得到的迎角α、侧滑角β,利用S4中公式计算得到修正系数Co和Cq,并由此计算得到静压P∞、动压q∞和修正后的空速Mai;
S6、判断Mai和Mai-1间的差值是否满足精度要求,满足则输出迎角α、侧滑角β、动压q∞和空速Mai,否则返回S2,循环执行S2-S6。
7.根据权利要求1所述的基于集成式大气数据模块的全速域大气数据解算方法,其特征在于,在S5中,所述动压q∞和空速Mai的计算方法为:
a)由步骤S4计算得到的无量纲修正系数Co、Cq,计算总压Po和静压P∞:
总压Po的计算公式为:
静压P∞的计算公式为:
动压q∞的计算公式为:
q∞=Po-P∞
b)空速Mai由上一循环Mai-1得到的计算公式得到:
若Mai-1≤0.3,马赫数的计算服从伯努利方程,空速Mai计算公式如下:
式中,q∞为动压,ρ为当地大气密度,a为当地音速。
若Mai-1>0.3,马赫数的计算服从一维等熵绝热流公式,空速Mai计算公式如下:
式中,γ为气体比热比常数。
8.根据权利要求1所述的基于集成式大气数据模块的全速域大气数据解算方法,其特征在于,所述气动数据库中搭载的数据包括:
风洞校准实验时得到的迎角αA、侧滑角βA、动压q∞ A、静压p∞ A以及马赫数MaA等大气数据,与模型表面压力Pi A特定组合得到的无量纲系数Cα A、Cβ A、Co A和Cq A,式中上标A表示为风洞标定试验结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110408075.3A CN113281001B (zh) | 2021-04-15 | 2021-04-15 | 基于集成式微型大气数据模块的全速域大气数据解算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110408075.3A CN113281001B (zh) | 2021-04-15 | 2021-04-15 | 基于集成式微型大气数据模块的全速域大气数据解算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113281001A true CN113281001A (zh) | 2021-08-20 |
CN113281001B CN113281001B (zh) | 2022-04-22 |
Family
ID=77277060
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110408075.3A Active CN113281001B (zh) | 2021-04-15 | 2021-04-15 | 基于集成式微型大气数据模块的全速域大气数据解算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113281001B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116380396A (zh) * | 2023-05-30 | 2023-07-04 | 中国航空工业集团公司沈阳空气动力研究所 | 一种连续式风洞飞行器大气数据系统鉴定试验系统及方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040011124A1 (en) * | 2000-11-17 | 2004-01-22 | Joel Choisnet | Method for determining aerodynamic parameters and method for detecting failure of a probe used for determining aerodynamic parameters |
CN105066994A (zh) * | 2015-08-21 | 2015-11-18 | 中国运载火箭技术研究院 | 一种嵌入式大气数据系统与惯性导航系统的数据融合方法 |
CN105301275A (zh) * | 2014-07-23 | 2016-02-03 | 空中客车运营简化股份公司 | 估算飞行器的马赫数的方法和装置 |
US10006928B1 (en) * | 2016-03-31 | 2018-06-26 | Textron Innovations Inc. | Airspeed determination for aircraft |
CN112461432A (zh) * | 2020-10-30 | 2021-03-09 | 成都凯天电子股份有限公司 | 一种实现大气数据测量长期稳定的在线修正方法 |
CN112580275A (zh) * | 2020-12-29 | 2021-03-30 | 成都信息工程大学 | 一种简单高效的高超声速飞行器fads系统测压孔布局方法 |
-
2021
- 2021-04-15 CN CN202110408075.3A patent/CN113281001B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040011124A1 (en) * | 2000-11-17 | 2004-01-22 | Joel Choisnet | Method for determining aerodynamic parameters and method for detecting failure of a probe used for determining aerodynamic parameters |
CN105301275A (zh) * | 2014-07-23 | 2016-02-03 | 空中客车运营简化股份公司 | 估算飞行器的马赫数的方法和装置 |
CN105066994A (zh) * | 2015-08-21 | 2015-11-18 | 中国运载火箭技术研究院 | 一种嵌入式大气数据系统与惯性导航系统的数据融合方法 |
US10006928B1 (en) * | 2016-03-31 | 2018-06-26 | Textron Innovations Inc. | Airspeed determination for aircraft |
CN112461432A (zh) * | 2020-10-30 | 2021-03-09 | 成都凯天电子股份有限公司 | 一种实现大气数据测量长期稳定的在线修正方法 |
CN112580275A (zh) * | 2020-12-29 | 2021-03-30 | 成都信息工程大学 | 一种简单高效的高超声速飞行器fads系统测压孔布局方法 |
Non-Patent Citations (2)
Title |
---|
朱亮聪 等: "基于部分遥测数据的飞行器力学环境评估方法", 《上海航天》 * |
陆辰: "嵌入式大气数据系统算法及大气数据传感信息融合关键技术研究", 《中国博士学位论文全文数据库(工程科技Ⅱ辑)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116380396A (zh) * | 2023-05-30 | 2023-07-04 | 中国航空工业集团公司沈阳空气动力研究所 | 一种连续式风洞飞行器大气数据系统鉴定试验系统及方法 |
CN116380396B (zh) * | 2023-05-30 | 2023-08-15 | 中国航空工业集团公司沈阳空气动力研究所 | 一种连续式风洞飞行器大气数据系统鉴定试验系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113281001B (zh) | 2022-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104318107B (zh) | 一种跨大气层飞行飞行器的高精度大气数据获取方法 | |
CN104048808B (zh) | 一种动态熵探针 | |
CN105122248A (zh) | 用于对部件尤其是涡轮叶片进行建模的方法 | |
CN111076940A (zh) | 一种测量亚音三维流场动态全参数的棱台五孔探针 | |
CN106545411B (zh) | 模拟畸变流场的直连试验装置的设计方法 | |
CN113281001B (zh) | 基于集成式微型大气数据模块的全速域大气数据解算方法 | |
CN114486164B (zh) | 一种基于神经网络的三孔探针梯度流场测试标定方法 | |
CN112685972A (zh) | 一种高效率的高超声速飞行器压力传感器分布优化方法 | |
CN114444216A (zh) | 基于数值模拟的高空条件下飞行器姿态控制方法及系统 | |
CN113051661A (zh) | 基于微小容腔多动力学的高温气流动态总压智能软测量方法 | |
CN111220341A (zh) | 一种风洞高马赫数低雷诺数轴对称型面喷管设计方法 | |
CN111498141B (zh) | 一种基于微型探针实现气流角度实时监测的方法与装置 | |
Peyvan et al. | Axial-flow compressor performance prediction in design and off-design conditions through 1-D and 3-D modeling and experimental study | |
Yang et al. | Aerodynamic characteristics of transonic and supersonic flow over rectangular cavities | |
CN109595040A (zh) | 一种燃气轮机叶片扭曲规律设计方法 | |
CN110836713B (zh) | 一种考虑校准箱气体质量变化的文氏管流量系数标定方法 | |
CN112883575A (zh) | 一种考虑表面粗糙度的叶轮机械边界层转捩模型修正方法 | |
CN111487031B (zh) | 基于微型扇形探针实现三维流场气流角度实时监测的装置与方法 | |
Sears et al. | Interference-free wind-tunnel flows by adaptive-wall technology | |
Torre et al. | The effect of turning angle on the loss generation of LP turbines | |
CN107202664A (zh) | 一种用于嵌入式大气数据系统的大气参数解算方法 | |
CN112597613B (zh) | 确定压气机叶片所受气动合力及其作用点的方法 | |
CN117744539B (zh) | 基于流动特征的高超飞行器初始流场生成方法 | |
CN117332512B (zh) | 一种飞机等校正空速爬升性能及加速因子的处理方法 | |
CN114077775B (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 |