CN107958206B - 一种飞行器表面热流辨识装置温度测量数据预处理方法 - Google Patents
一种飞行器表面热流辨识装置温度测量数据预处理方法 Download PDFInfo
- Publication number
- CN107958206B CN107958206B CN201711086206.0A CN201711086206A CN107958206B CN 107958206 B CN107958206 B CN 107958206B CN 201711086206 A CN201711086206 A CN 201711086206A CN 107958206 B CN107958206 B CN 107958206B
- Authority
- CN
- China
- Prior art keywords
- temperature
- data
- measurement data
- heat flow
- flow identification
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K7/00—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
- G01K7/02—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using thermoelectric elements, e.g. thermocouples
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明提供了一种飞行器表面热流辨识装置温度测量数据预处理方法,属于航空航天飞行试验热学参数测量及处理技术领域。该方法首先对热流辨识装置的温升测量数据进行局部失真点(局部跳点)进行剔除的处理,然后利用N个相邻数据点平均的光滑处理方法对测量数据进行平滑处理,最后得到满足热流辨识要求的温度测量数据。所述N值根据温度传感器相关参数和温度曲线特征进行确定。使用本发明完成预处理后的温度测量数据进行热流辨识,可以有效改善温度阶跃和局部跳点对热流辨识结果的影响,提高热流辨识结果的准确度和可靠性。
Description
技术领域
本发明涉及一种飞行器表面热流辨识装置温度测量数据预处理方法,适用于目前高超声速飞行器飞行试验的测量数据的预处理,属于航空航天飞行试验热学参数测量及处理技术领域。
背景技术
高超声速飞行器表面气动加热环境严酷,飞行器防隔热***面临严峻的挑战,对飞行器防隔热***设计精度和保守量提出了苛刻的要求。而准确预测气动加热环境及分布则是实现防隔热***精细化设计的前提和根本。目前受制于气动热环境预测技术的现状,热环境预示方法尚不能达到很好的精度,而地面试验设备的能力目前不能复现真实飞行环境。在这种背景下,通过飞行试验开展飞行器表面热流测量进而完善热环境预示方法就变得日益迫切。
然而由于高超声速飞行器表面热环境非常严酷,且飞行时间长,受制于传感器耐温极限、产品尺寸和测量精度,目前利用传统热流传感器直接进行测量的难度很大,逐渐难以适应长时间高热流的测量需求。气动热辨识是一种热流间接测量方法,其利用测量防热材料层间温升或通过测量热流辨识装置的温升,通过数学方法反演得到表面热流信息。由于环境适应性更好、测量结构相对简单、产品尺寸相对较小,气动热辨识方法和相应的热流辨识装置得到了越来越多的研究和应用。
由热流辨识的原理可知,敏感元件温升测量数据是进行热流辨识的前提和关键。敏感元件温升信息测量过程为:首先通过安装于敏感元件内部的热电偶获得电信号信息,然后通过放大器和变换器转换成为温度信息。由于变换器的采编位数有限,例如目前航天飞行器使用的传感器的采编位数一般为8位,即在量程范围内分为255个分度,在测量值明显小于量程最大值的情况下,温度测量数据会产生明显的阶跃现象,如图2(a)所示,温度传感器量程为0~200℃,每个分度值约为0.8℃,由于测量值最大不超过30℃,温升测量结果存在明显阶跃现象。另外,电信号在放大和变换的过程中由于受到电磁干扰、电源稳定性和飞行器瞬时振动和冲击的影响,温度测量结果可能存在局部失真情况,图2(b)给出了某飞行试验温度测量结果,可以看出温度测量结果局部存在跳点。温度测量结果的阶跃现象和局部跳点均会对热流辨识过程产生影响,导致热流辨识结果失真,因此需要在开展热流辨识工作前对温度测量结果进行预处理,为热流辨识提供合理的输入条件,提高热流辨识结果的准确性。
目前对于高超声速飞行试验具有阶跃特征的非光滑温度测量数据还缺乏有效的预处理方法。采用常用的数据拟合方法和数据平滑方法开展数据预处理研究,结果表明:数据拟合方法如多项式拟合方法,其处理结果曲线非常光滑,从宏观上看与原始数据曲线吻合较好,但会抹平或忽略温升曲线的细小变化,导致热流辨识结果失真;数据平滑方法如各类平滑滤波方法,由于难以确定滤波方法的一些特征参数(如所用的数据点个数,滤波的数据宽度等),其对具有阶跃特征的测量数据的处理效果并不令人满意。
发明内容
本发明的技术解决问题是:针对目前热流辨识组件测量数据质量不高、难以作为热流辨识有效输入条件的问题,提供一种飞行器表面热流辨识装置温度测量数据预处理方法,可以为热流辨识提供合理的温度数据输入,有效改善温度阶跃和局部跳点对热流辨识结果的影响,提高热流辨识结果的准确度和可靠性。
其中,nmax为所有数据点个数,N为预设的平滑处理点数,为正整数。
作为优选方案,还可以在所述步骤之前增加如下处理:
分析飞行器表面热流辨识装置测量所得到的温度曲线变化规律和趋势,剔除温度测量数据中的跳点。
所述数据点个数N通过如下方法确定:
(1)、对温度测量数据曲线进行统计分析,确定温度测量数据的温度变化范围ΔT和测量时长为Δt;
(2)、根据温度测量数据测量所采用的温度传感器的量程TR和采编位数为Nd,计算温度传感器测量分度值ΔT;
(3)、根据温度变化范围ΔT、温度传感器测量分度值ΔT,计算温度测量数据的台阶个数mstep;
(4)、根据测量时长Δt、温度测量数据的台阶个数mstep和采编器采样率f,计算预处理平滑数据点个数N。
本发明与现有技术相比的有益效果是:
(1)、本发明提出对飞行试验热流辨识装置的温度测量数据采用相邻数据点平均的方法进行平滑处理,提高其作为热流辨识输入条件的数据质量;
(2)、本发明进一步提出了平滑数据点数N的一种确定方法,利用根据该平滑数据点数N完成预处理后的温度测量数据进行热流辨识,可以有效改善温度阶跃和局部跳点对热流辨识结果的影响,提高热流辨识结果的准确度和可靠性。
附图说明
图1为本发明温度测量数据预处理方法流程图;
图2(a)为本发明实施例温度测量数据曲线;
图2(b)为本发明实施例存在跳点的温度测量数据曲线;
图3为本发明实施例温度变化范围ΔT和测量时长Δt定义示意图;
图4(a)为本发明实施例原始数据、多项式拟合数据和平滑数据的整体图;
图4(b)为本发明实施例原始数据、多项式拟合数据和平滑数据的局部图;
图5(a)为本发明实施例原始数据、不同点数数据平滑数据的整体图;
图5(b)为本发明实施例原始数据、不同点数平滑数据的局部放大图;
图6为本发明实施例基于不同点数平滑数据的辨识结果图;
图7为本发明基于原始数据、多项式拟合数据和平滑数据的辨识结果与飞行攻角的对比。
具体实施方式
以下结合附图和具体实施方式对本发明进行详细说明。
针对目前热流辨识装置测量数据预处理方面存在的问题,根据上述分析和研究,本发明提供了一种飞行器表面热流辨识装置温度测量数据预处理方法。该方法的技术特点为:首先对热流辨识装置的温升测量数据进行局部失真点(局部跳点)剔除的处理,然后利用N个相邻数据点平均的光滑处理方法对测量数据进行平滑处理(N值根据传感器相关参数和温度曲线特征进行确定),最后得到满足热流辨识要求的温度测量数据。具体过程如下:
s1、分析飞行器表面热流辨识装置测量所得到的温度曲线变化规律和趋势,剔除温度测量数据的中的跳点。
其中,nmax为所有数据点个数。
作为优选方案,本发明进一步提出根据温度传感器量程TR、采编器采编位数Nd、采编器采样率f以及温度变化范围ΔT和测量时长Δt,确定预处理平滑数据点个数N,具体为:
(1)、对温度测量数据曲线进行分析和统计,确定温度测量数据的温度变化范围ΔT和测量时长为Δt;
图3为典型飞行试验热流辨识装置温度测量数据曲线温度变化范围ΔT和测量时长Δt定义示意图,其中横坐标为时间t,纵坐标为温度T,图中曲线的温升为ΔT,测量时长为Δt。
温度测量数据的起始时间为t1,结束时间为t2,起始温度为T1,结束时间为T2,则温升为ΔT=|T2-T1|,测量时长为Δt=t2-t1。
(2)、根据温度传感器量程TR,计算温度传感器测量分度值ΔT:
(3)、根据温度变化范围ΔT、温度传感器测量分度值ΔT,计算温度测量数据的台阶个数mstep:
(4)、根据温度测量数据的测量时长Δt、温度测量数据的台阶个数mstep和采编器采样率f,计算预处理平滑数据点个数N:
实施例:
图4(a)和图4(b)分别给出了典型飞行试验热流辨识装置原始数据、多项式拟合数据和平滑数据整体和局部曲线图,从图中可以看出,从宏观上看,两种方法处理结果均与原始数据曲线吻合较好,但从局部来看,本发明预处理结果与原始数据更吻合,更能精确地反应原始数据的变化。
图5(a)和图5(b)分别给出了典型飞行试验热流辨识装置原始温度测量数据和不同点数数据平滑数据的整体和局部放大图。其中,给出了采用N=200、400、800个相邻数据点进行光滑处理的温度结果,从宏观上看不同点数光滑处理结果几乎一致,从局部放大图可看出三者存在细微的差别。图6给出了不同点数相应的热流辨识结果的对比,由对比可知,采用不同数量的数据点进行数据光滑处理,得到热流辨识结果有明显差异。若光滑处理所使用的数据点数不够,则不足以改善热流辨识结果的剧烈跳变现象,如N=200的热流辨识结果所示;而光滑处理所使用的数据点数越多,其对热流剧烈变化的抹平效果也越明显,当N=800时,热流的阶梯变化已被抹平。所以对于数据点光滑的处理方法来说,采用多少数据点(即N值的选取)进行光滑是一个关键的问题。由此可见,对温升曲线进行光滑预处理时,要合理选取数据光滑预处理时所采用的数据点个数。
图7为基于原始数据、多项式拟合数据和平滑数据的热流辨识结果的对比,图中还给出了飞行攻角的分布,可以看出,原始数据由于温度阶跃的存在,热流辨识结果呈现出剧烈的跳变现象、严重失真,预处理后数据的辨识结果则有效消除了原始测量数据辨识结果的跳变现象,而且热流辨识结果的幅值和变化趋势得到了较好的保留,热流随攻角有明显的阶梯状变化,符合热流随攻角改变而变化这一真实物理现象;多项式拟合数据的热流辨识结果过于光滑,热流幅值和变化趋势均被抹平,存在较大的失真现象。
由上述分析可知,本发明数据平滑的方法相对能更好地反映温升曲线的变化,热流辨识结果的准确度和可靠性更高。
本发明提供的测量数据预处理方法可以为热流辨识提供合理的温度数据输入,有效改善温度阶跃和局部跳点对热流辨识结果的影响,提高热流辨识结果的准确度和可靠性。
本发明提供的数据预处理方法不仅仅适用于热流辨识装置的测量数据,还适用于其它类型测温装置获得的温度测量数据以及其它热相关测量数据。上述实施例子只是对本发明的解释,而不能作为对本发明的限制,因此凡是与本发明思路类似的实施方式均在本发明的保护范围内。
本说明书中未进行详细描述部分属于本领域技术人员公知常识。
Claims (4)
其中,nmax为所有数据点个数,N为预设的平滑处理点数,为正整数;
所述数据点个数N通过如下方法确定:
(1)、对温度测量数据曲线进行统计分析,确定温度测量数据的温度变化范围ΔT和测量时长为Δt;
(2)、根据温度测量数据测量所采用的温度传感器的量程TR和采编位数为Nd,计算温度传感器测量分度值ΔT;
(3)、根据温度变化范围ΔT、温度传感器测量分度值ΔT,计算温度测量数据的台阶个数mstep;
(4)、根据测量时长Δt、温度测量数据的台阶个数mstep和采编器采样率f,计算预处理平滑数据点个数N:
所述预处理平滑数据点个数N具体的计算公式为:
2.根据权利要求1所述的一种飞行器表面热流辨识装置温度测量数据预处理方法,其特征在于在所述步骤之前增加如下处理:
分析飞行器表面热流辨识装置测量所得到的温度曲线变化规律和趋势,剔除温度测量数据中的跳点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711086206.0A CN107958206B (zh) | 2017-11-07 | 2017-11-07 | 一种飞行器表面热流辨识装置温度测量数据预处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711086206.0A CN107958206B (zh) | 2017-11-07 | 2017-11-07 | 一种飞行器表面热流辨识装置温度测量数据预处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107958206A CN107958206A (zh) | 2018-04-24 |
CN107958206B true CN107958206B (zh) | 2021-02-09 |
Family
ID=61964466
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711086206.0A Active CN107958206B (zh) | 2017-11-07 | 2017-11-07 | 一种飞行器表面热流辨识装置温度测量数据预处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107958206B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109598342B (zh) * | 2018-11-23 | 2021-07-13 | 中国运载火箭技术研究院 | 一种决策网络模型自博弈训练方法及*** |
CN109556762B (zh) * | 2018-12-07 | 2023-08-29 | 中国航天空气动力技术研究院 | 一种用于测量凹腔热流分布的探头 |
CN117451217B (zh) * | 2023-12-25 | 2024-03-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种航天热流传感器及基于双温差补偿的热流修正方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105068032A (zh) * | 2015-08-13 | 2015-11-18 | 湖南大学 | 一种光伏汇流箱电流采集通道温度漂移系数的校准方法 |
CN105548250A (zh) * | 2016-01-13 | 2016-05-04 | 北京机电工程研究所 | 飞行器气动热试验的热流控制方法、装置及*** |
CN106021181A (zh) * | 2016-05-17 | 2016-10-12 | 中国电子科技集团公司第四十研究所 | 一种对fft数据实施的相邻点相关性均值降噪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9435694B2 (en) * | 2013-03-14 | 2016-09-06 | Kelsey-Hayes Company | Outside air temperature measurement device and method |
-
2017
- 2017-11-07 CN CN201711086206.0A patent/CN107958206B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105068032A (zh) * | 2015-08-13 | 2015-11-18 | 湖南大学 | 一种光伏汇流箱电流采集通道温度漂移系数的校准方法 |
CN105548250A (zh) * | 2016-01-13 | 2016-05-04 | 北京机电工程研究所 | 飞行器气动热试验的热流控制方法、装置及*** |
CN106021181A (zh) * | 2016-05-17 | 2016-10-12 | 中国电子科技集团公司第四十研究所 | 一种对fft数据实施的相邻点相关性均值降噪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107958206A (zh) | 2018-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107958206B (zh) | 一种飞行器表面热流辨识装置温度测量数据预处理方法 | |
CN102128990B (zh) | 一种电磁辐射敏感度自动测试中判断功率放大器饱和的方法 | |
CN109668506B (zh) | 一种基于涡流稳态特性的磁性金属材料厚度检测方法 | |
CN103344174A (zh) | 一种不均匀导电材料表面覆盖层厚度的涡流测厚方法 | |
CN105548917A (zh) | 一种非屏蔽环境下磁传感器阵列的校准方法 | |
CN201653844U (zh) | 热障涂层抗高温氧化性能测试装置 | |
Sazhin | Novel mass air flow meter for automobile industry based on thermal flow microsensor. II. Flow meter, test procedures and results | |
US20130088222A1 (en) | System and method for measuring wrinkle depth in a composite structure | |
CN108169565B (zh) | 一种用于电导率测量的非线性温度补偿方法 | |
CN113587794A (zh) | 电磁轴承转子位移测量方法、装置以及电磁轴承*** | |
CN103698724A (zh) | 一种多铁性材料的磁电性能测量***和测量方法 | |
CN110618332A (zh) | 基于热流测量的电容器温升测量方法及*** | |
CN103293506A (zh) | 一种无需检测前试块标定的电导率测量仪器的实现方法 | |
CN104006903A (zh) | 发动机进气口温度传感器的检测方法 | |
CN105319444A (zh) | 一种导电材料电导率均匀程度评估方法 | |
CN113296045B (zh) | 基于传感阵列的微型电流传感器的误差校正方法 | |
CN113156349B (zh) | 一种材料磁致机械特性的测量方法和装置 | |
CN111999377B (zh) | 漏磁检测表征缺陷宽度的方法 | |
CN112666395B (zh) | 非接触式的金属材质电导率测量方法及*** | |
CN108051648A (zh) | 一种基于直流电位和涡流检测法的材料电磁属性测量方法 | |
Demetriades | Turbulence correlations in a compressible wake | |
CN107563100B (zh) | 基于黑盒理论分析接触电阻变化规律的方法 | |
Bakken et al. | A velocity dependent effective angle method for calibration of X-probes at low velocities | |
CN111044604A (zh) | 一种acfm单轴磁信号评估方法 | |
Wang et al. | Vibration measurement using a low-cost imaging sensor and image processing techniques |
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 |