CN113532722B - 一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法 - Google Patents
一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法 Download PDFInfo
- Publication number
- CN113532722B CN113532722B CN202110571918.1A CN202110571918A CN113532722B CN 113532722 B CN113532722 B CN 113532722B CN 202110571918 A CN202110571918 A CN 202110571918A CN 113532722 B CN113532722 B CN 113532722B
- Authority
- CN
- China
- Prior art keywords
- pressure
- pulsating pressure
- double
- time period
- aircraft
- 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
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L9/00—Measuring steady of quasi-steady pressure of fluid or fluent solid material by electric or magnetic pressure-sensitive elements; Transmitting or indicating the displacement of mechanical pressure-sensitive elements, used to measure the steady or quasi-steady pressure of a fluid or fluent solid material, by electric or magnetic means
- G01L9/02—Measuring steady of quasi-steady pressure of fluid or fluent solid material by electric or magnetic pressure-sensitive elements; Transmitting or indicating the displacement of mechanical pressure-sensitive elements, used to measure the steady or quasi-steady pressure of a fluid or fluent solid material, by electric or magnetic means by making use of variations in ohmic resistance, e.g. of potentiometers, electric circuits therefor, e.g. bridges, amplifiers or signal conditioning
- G01L9/04—Measuring steady of quasi-steady pressure of fluid or fluent solid material by electric or magnetic pressure-sensitive elements; Transmitting or indicating the displacement of mechanical pressure-sensitive elements, used to measure the steady or quasi-steady pressure of a fluid or fluent solid material, by electric or magnetic means by making use of variations in ohmic resistance, e.g. of potentiometers, electric circuits therefor, e.g. bridges, amplifiers or signal conditioning of resistance-strain gauges
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L19/00—Details of, or accessories for, apparatus for measuring steady or quasi-steady pressure of a fluent medium insofar as such details or accessories are not special to particular types of pressure gauges
- G01L19/06—Means for preventing overload or deleterious influence of the measured medium on the measuring device or vice versa
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
Description
技术领域
本发明涉及一种双谱分析转捩辨识方法。
背景技术
脉动压力是飞行器设计中的一种重要物理量。高速飞行器容易遭受强烈的压力脉动,这些压力脉动可以引起舱内元件的震动并且导致结构上的问题,从而影响着飞行器的性能。根据本单位的飞行试验结果,飞行器随着飞行高度降低,来流雷诺数增大,飞行器表面的边界层由层流过渡到湍流过程中,发现脉动压力幅值并非连续增大。针对飞行试验脉动压力数据的转捩辨识,以往的公开文献鲜有报道。然而,单纯地通过压力脉动幅值也并不能准确地判读转捩发生时刻。
现有的飞行试验转捩测量技术多数使用热流传感器与温度传感器来判别转捩,但是热流与温度传感器频响较低(一般在几十Hz到几百Hz),无法兼顾测量边界层内扰动波(一般在kHz)的特性。
发明内容
本发明要解决的技术问题是:本发明提出一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法,能够较为准确地辨识出飞行试验中飞行器边界层的转捩时刻,对后续飞行器的弹道优化、热防护设计提供数据支撑。
本发明所采用的技术方案是:一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法,包括步骤如下:
S2:对脉动压力在时间域上的分布曲线进行划分,截取若干个Δt时间段内的脉动压力数据;Δt的取值范围根据采样频率选取;
S3:对Δt时间段内的脉动压力数据进行双谱分析,得到每个Δt时间段内的脉动压力双谱值;
S4:根据不同时间段内的脉动压力双谱值大小,辨识飞行器转捩发生时刻。
飞行器上的压力传感器选取压阻式Kulite传感器或压电式PCB传感器。
S1中,利用matlab软件中的smooth函数进行移动平均滤波,其中点数范围SPAN根据实际测量中的采样频率来确定,滤波方法采用Lowess。
S3的具体方法为:
(a)利用matlab软件中的离散小波函数wavedec对每个Δt时间段内的脉动压力数据进行小波分解;其中,wavedec函数中的小波函数wname根据样本数据具体形式进行选择;
(b)对每个尺度频率下的脉动压力数据进行互相关分析,获取每个Δt时间段内的脉动压力的相关性大小,即双谱值大小。
一种基于飞行试验脉动压力数据的双谱分析转捩辨识系统,包括:
第一模块,采用移动平均方法对飞行器上压力传感器测得的瞬时压力p进行滤波,得到平均压力获取飞行器的脉动压力在时间域上的分布曲线;对脉动压力在时间域上的分布曲线进行划分,截取若干个Δt时间段内的脉动压力数据;Δt的取值范围根据采样频率选取;
第二模块,对Δt时间段内的脉动压力数据进行双谱分析,得到每个Δt时间段内的脉动压力双谱值;
第三模块,根据不同时间段内的脉动压力双谱值大小,辨识飞行器转捩发生时刻。
飞行器上的压力传感器选取压阻式Kulite传感器或压电式PCB传感器。
第一模块中,利用matlab软件中的smooth函数进行移动平均滤波,其中,点数范围SPAN根据实际测量中的采样频率来确定,滤波方法采用Lowess。
对Δt时间段内的脉动压力数据进行双谱分析的具体方法为:
(a)利用matlab软件中的离散小波函数wavedec对每个Δt时间段内的脉动压力数据进行小波分解;其中,wavedec函数中的小波函数wname根据样本数据具体形式进行选择;
(b)对每个尺度频率下的脉动压力数据进行互相关分析,获取每个Δt时间段内的脉动压力的相关性大小,即双谱值大小。
本发明与现有技术相比的优点在于:
本发明提供一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法,其特征在于,针对飞行器脉动压力在时域上的分布特征,不是直观的从脉动压力幅值来判读转捩时刻。根据边界层发展演化过程中的扰动波特征,利用双谱方法的最大相关系数对边界层转捩时刻进行辨识。该方法已被飞行试验的热流与温度数据佐证,已应用于多个外形飞行器的转捩时刻辨识。
现有的飞行试验转捩测量技术多数使用热流传感器与温度传感器来判别转捩,但是热流与温度传感器频响较低(一般在几十Hz到几百Hz),无法兼顾测量边界层内扰动波(一般在kHz)的特性。脉动压力传感器的高频响特性不仅可以直观地测量边界层内压力的幅值变化,同时也可反映边界层内扰动波的发展演化特征。
附图说明
图1为本发明实施例公开的基于飞行试验脉动压力数据的转捩辨识方法的流程图;
图2为本发明实施例公开的飞行器试验不同时间段的脉动压力相关性脉动压力在时间域上的分布曲线示意图。
图3为本发明实施例公开的飞行器脉动压力在时间域上的分布曲线划分示意图。
图4为脉动压力数据双谱分析结果图。
具体实施方式
结合附图对本发明进行说明。
根据飞行试验测量的瞬时压力数据,对其进行滤波,获得脉动压力在时域上的分布。针对脉动压力数据,对其进行不同时间段的划分,利用小波变换方法对不同时间段的脉动压力数据进行多尺度分解,然后对分解后的每个尺度脉动压力进行互相关分析。根据不同时间段的双谱值大小,判断边界层转捩发生时刻。
如图1所示,一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法,包括步骤如下:
飞行试验转捩测量所用的压力传感器一般选取压阻式Kulite传感器和压电式PCB传感器。
S101中,利用matlab软件中的smooth函数进行移动平均滤波,其中点数范围SPAN根据实际测量中的采样频率来确定,滤波方法采用Lowess。
S102:对脉动压力在时间域上的分布曲线进行划分,截取若干个Δt时间段内的脉动压力数据。Δt的取值范围根据采样频率选取(这里采样频率为200kHz,Δt=0.1s);
S103:对Δt时间段内的脉动压力数据进行双谱分析;
小波变换为一种具有多分辨分析特点的时频局部化分析方法,其在低频部分具有较低的时间分辨率和较高的频率分辨率,在高频部分具有较高的时间分辨率和较低的频率分辨率。基于傅立叶变换的双谱方法无法对飞行试验中测量的非稳态信号进行处理,小波变换适合分析非平稳的信号和提取信号的局部特征,且基于小波变换的双谱方法可根据测量数据的分布特征选取与之匹配的小波基函数。
具体方法为:
(a)利用matlab软件中的离散小波函数wavedec对每个Δt时间段内的脉动压力数据进行小波分解。其中,wavedec函数中的小波函数wname无明确要求,可根据样本数据具体形式进行选择。
(b)对每个尺度频率下的脉动压力数据进行互相关分析,获取每个Δt时间段内的脉动压力的相关性大小,即双谱值大小。
S104:根据不同时间段内的脉动压力双谱值大小,辨识飞行器转捩发生时刻。
一种基于飞行试验脉动压力数据的双谱分析转捩辨识系统,包括:
第一模块,采用移动平均方法对飞行器上压力传感器测得的瞬时压力p进行滤波,得到平均压力获取飞行器的脉动压力在时间域上的分布曲线;对脉动压力在时间域上的分布曲线进行划分,截取若干个Δt时间段内的脉动压力数据;Δt的取值范围根据采样频率选取;
第二模块,对Δt时间段内的脉动压力数据进行双谱分析,得到每个Δt时间段内的脉动压力双谱值;
第三模块,根据不同时间段内的脉动压力双谱值大小,辨识飞行器转捩发生时刻。
飞行器上的压力传感器选取压阻式Kulite传感器或压电式PCB传感器。
第一模块中,利用matlab软件中的smooth函数进行移动平均滤波,其中,点数范围SPAN根据实际测量中的采样频率来确定,滤波方法采用Lowess。
对Δt时间段内的脉动压力数据进行双谱分析的具体方法为:
(a)利用matlab软件中的离散小波函数wavedec对每个Δt时间段内的脉动压力数据进行小波分解;其中,wavedec函数中的小波函数wname根据样本数据具体形式进行选择;
(b)对每个尺度频率下的脉动压力数据进行互相关分析,获取每个Δt时间段内的脉动压力的相关性大小,即双谱值大小。
实施例:
一种基于飞行试验脉动压力数据的转捩辨识方法,从工程应用上能够较为准确地辨识出飞行试验中飞行器边界层的转捩时刻,为工程转捩预示方法精度的提升提供重要技术支撑。
本发明基于某外形飞行试验的脉动压力压力测量结果,获取不同时间段的双谱结果,根据双谱值大小,辨识判读转捩发生时刻。
一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法,包括步骤如下:
S101:采用中值滤波对瞬时压力进行滤波得到平均压力,获取压力的脉动在时间域上的分布,如图2所示。
S102:对脉动压力在时间域上进行划分,截取0.1s内的脉动压力数据(采用频率为200kHz),如图3所示。
S103:对划分时间段的脉动压力数据进行双谱分析,即对每个时间段内的脉动压力数据进行小波变换,对每个尺度频率下的数据进行互相关分析,获取每个时间段内的脉动压力双谱值大小,如图4所示。
S104:转捩时刻判读,即根据不同时间段内的脉动压力双谱值大小,辨识转捩发生时刻。
本发明未详细说明部分属于本领域技术人员公知技术。
Claims (4)
1.一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法,其特征在于,包括步骤如下:
S1中,利用matlab软件中的smooth函数进行移动平均滤波,其中点数范围SPAN根据实际测量中的采样频率来确定,滤波方法采用Lowess;
S2:对脉动压力在时间域上的分布曲线进行划分,截取若干个Δt时间段内的脉动压力数据;Δt的取值范围根据采样频率选取;
S3:对Δt时间段内的脉动压力数据进行双谱分析,得到每个Δt时间段内的脉动压力双谱值;
S3的具体方法为:
(a)利用matlab软件中的离散小波函数wavedec对每个Δt时间段内的脉动压力数据进行小波分解;其中,wavedec函数中的小波函数wname根据样本数据具体形式进行选择;
(b)对每个尺度频率下的脉动压力数据进行互相关分析,获取每个Δt时间段内的脉动压力的相关性大小,即双谱值大小;
S4:根据不同时间段内的脉动压力双谱值大小,辨识飞行器转捩发生时刻。
2.根据权利要求1所述的一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法,其特征在于,飞行器上的压力传感器选取压阻式Kulite传感器或压电式PCB传感器。
3.一种基于飞行试验脉动压力数据的双谱分析转捩辨识系统,其特征在于,包括:
第一模块,采用移动平均方法对飞行器上压力传感器测得的瞬时压力p进行滤波,得到平均压力获取飞行器的脉动压力在时间域上的分布曲线;对脉动压力在时间域上的分布曲线进行划分,截取若干个Δt时间段内的脉动压力数据;Δt的取值范围根据采样频率选取;
第一模块中,利用matlab软件中的smooth函数进行移动平均滤波,其中,点数范围SPAN根据实际测量中的采样频率来确定,滤波方法采用Lowess;第二模块,对Δt时间段内的脉动压力数据进行双谱分析,得到每个Δt时间段内的脉动压力双谱值;
对Δt时间段内的脉动压力数据进行双谱分析的具体方法为:
(a)利用matlab软件中的离散小波函数wavedec对每个Δt时间段内的脉动压力数据进行小波分解;其中,wavedec函数中的小波函数wname根据样本数据具体形式进行选择;
(b)对每个尺度频率下的脉动压力数据进行互相关分析,获取每个Δt时间段内的脉动压力的相关性大小,即双谱值大小;
第三模块,根据不同时间段内的脉动压力双谱值大小,辨识飞行器转捩发生时刻。
4.根据权利要求3所述的一种基于飞行试验脉动压力数据的双谱分析转捩辨识系统,其特征在于,飞行器上的压力传感器选取压阻式Kulite传感器或压电式PCB传感器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110571918.1A CN113532722B (zh) | 2021-05-25 | 2021-05-25 | 一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110571918.1A CN113532722B (zh) | 2021-05-25 | 2021-05-25 | 一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113532722A CN113532722A (zh) | 2021-10-22 |
CN113532722B true CN113532722B (zh) | 2023-04-14 |
Family
ID=78094722
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110571918.1A Active CN113532722B (zh) | 2021-05-25 | 2021-05-25 | 一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113532722B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115840889B (zh) * | 2023-02-17 | 2023-05-26 | 中国空气动力研究与发展中心计算空气动力研究所 | 用于转捩预测的特征值的处理方法、装置、设备及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110806300A (zh) * | 2019-10-12 | 2020-02-18 | 北京临近空间飞行器系统工程研究所 | 一种适用于高超声速飞行试验转捩研究的测点布置方法 |
CN111207903A (zh) * | 2020-03-02 | 2020-05-29 | 北京空天技术研究所 | 适用于亚跨超音速风洞的转捩测量方法 |
CN112304563A (zh) * | 2020-10-30 | 2021-02-02 | 中国空气动力研究与发展中心超高速空气动力研究所 | 研究转捩对高超声速飞行器气动特性影响的风洞试验方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2648567B1 (fr) * | 1989-05-24 | 1991-10-04 | Inst Nat Sante Rech Med | Procede de traitement numerique d'un signal par transformation reversible en ondelettes |
US7251592B1 (en) * | 2003-08-18 | 2007-07-31 | United Technologies Corporation | Boundary layer transition model |
CN104730518B (zh) * | 2015-03-30 | 2017-03-15 | 北京空间飞行器总体设计部 | 一种基于高斯拟合的雷达多普勒谱估计海面流场的方法 |
US10953979B2 (en) * | 2015-11-11 | 2021-03-23 | The Arizona Board Of Regents On Behalf Of The University Of Arizona | Control of hypersonic boundary layer transition |
CN109141731A (zh) * | 2018-07-26 | 2019-01-04 | 西北工业大学 | 一种可用于水下湍流边界层壁面脉动压力测试的柔性基微传感器及其制造方法 |
CN111090907B (zh) * | 2019-10-30 | 2023-09-29 | 中国航天空气动力技术研究院 | 一种飞行试验转捩判断方法 |
CN111780948B (zh) * | 2020-06-10 | 2022-01-04 | 北京临近空间飞行器系统工程研究所 | 一种高超声速飞行试验中飞行器边界层转捩过程特性的测量方法 |
CN111832159B (zh) * | 2020-06-23 | 2023-08-29 | 北京临近空间飞行器系统工程研究所 | 一种基于飞行试验数据的边界层转捩阵面动态演化过程确定方法 |
CN111766039B (zh) * | 2020-07-10 | 2022-04-15 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种亚声速风洞可压缩流体扰动模态测量结果解算方法 |
CN112484954A (zh) * | 2020-11-24 | 2021-03-12 | 中国航天空气动力技术研究院 | 一种判断流场状态信息的方法、系统及存储介质 |
-
2021
- 2021-05-25 CN CN202110571918.1A patent/CN113532722B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110806300A (zh) * | 2019-10-12 | 2020-02-18 | 北京临近空间飞行器系统工程研究所 | 一种适用于高超声速飞行试验转捩研究的测点布置方法 |
CN111207903A (zh) * | 2020-03-02 | 2020-05-29 | 北京空天技术研究所 | 适用于亚跨超音速风洞的转捩测量方法 |
CN112304563A (zh) * | 2020-10-30 | 2021-02-02 | 中国空气动力研究与发展中心超高速空气动力研究所 | 研究转捩对高超声速飞行器气动特性影响的风洞试验方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113532722A (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110455490B (zh) | 超声速和高超声速风洞流场湍流度的计算方法及装置 | |
Thomas et al. | On the mechanism of unsteady shock oscillation in shock wave/turbulent boundary layer interactions | |
Brusniak et al. | Physics of unsteady blunt-fin-induced shock wave/turbulent boundary layer interactions | |
CN113532722B (zh) | 一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法 | |
CN106247173B (zh) | 管道泄漏检测的方法及装置 | |
CN107014668A (zh) | 一种基于压电和智能涂层传感器的疲劳裂纹综合监测方法 | |
Nolan et al. | Particle image velocimetry measurements of a transitional boundary layer under free stream turbulence | |
CN105953080B (zh) | 基于同侧传感器布置的声波泄漏定位方法 | |
CN104595729A (zh) | 一种基于声波幅值的油气管道泄漏定位方法 | |
CN107963239A (zh) | 一种基于音频的运载火箭故障检测装置及检测方法 | |
CN113494990A (zh) | 一种风洞扰动对超声速层流边界层厚度影响的分析方法 | |
CN104545870A (zh) | 用于心率检测的床垫和心率检测方法 | |
CN110953488B (zh) | 基于堆栈自编码的气液两相流管道泄漏声发射检测方法 | |
CN109916434B (zh) | 一种光纤法珀传感器腔长多倍频相关解调方法 | |
CN108992053A (zh) | 一种实时的无束缚检测心率和心跳间隔的方法 | |
CN108107086A (zh) | 一种基于阵列气体传感器的气体检测方法及气体传感器 | |
CN104698080A (zh) | 一种利用兰姆波对结构损伤进行状态监测的方法 | |
CN108646248A (zh) | 一种针对低速运动声源的被动声学测速测距方法 | |
CN104696711A (zh) | 管道泄漏点快速准确定位的方法 | |
CN112857730B (zh) | 一种高超声速脉动压力试验数据分析处理方法 | |
CN110186521A (zh) | 基于小波脊特征提取的涡街湿气过读补偿与流量测量方法 | |
CN103940575B (zh) | 一种基于应变信号能量的风洞试验天平的评估方法 | |
EP1668320A2 (en) | Method and system for locating critical flow feature indicators in three dimensional flow regimes | |
CN109163894B (zh) | 一种基于摩擦温度信号的磨合状态识别方法 | |
CN116738152A (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 |