CN114088335A - 基于气动探针正反向连续运动的定常流场快速测量方法 - Google Patents
基于气动探针正反向连续运动的定常流场快速测量方法 Download PDFInfo
- Publication number
- CN114088335A CN114088335A CN202210029277.1A CN202210029277A CN114088335A CN 114088335 A CN114088335 A CN 114088335A CN 202210029277 A CN202210029277 A CN 202210029277A CN 114088335 A CN114088335 A CN 114088335A
- Authority
- CN
- China
- Prior art keywords
- pneumatic probe
- pressure
- pneumatic
- probe
- pressure data
- 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
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/06—Measuring arrangements specially adapted for aerodynamic testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L11/00—Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L27/00—Testing or calibrating of apparatus for measuring fluid pressure
- G01L27/002—Calibrating, i.e. establishing true relation between transducer output value and value to be measured, zeroing, linearising or span error determination
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
本发明公开了一种基于气动探针正反向连续运动的定常流场快速测量方法。本发明属于气动测试技术领域。该测量方法由气动探针正向和反向连续运动获得两组压力测量数据,首先进行无相位偏移的低通滤波处理,再设定离散***压力递推公式的传递函数的系数进行修正,获取两组修正后的压力测量数据,随后利用连续同步采集的气动探针位置信息插值,获取相同的空间位置处的修正后的压力测量数据,最后以相同的空间位置的正向修正后的压力测量数据和反向修正后的压力测量数据的平均值作为气动探针的最终的压力测量数据。该测量方法能够显著提高气动探针的定常流场气动参数的测试效率,降低试验成本,无新增硬件设备,实施便捷且具有较强的环境适应性。
Description
技术领域
本发明属于气动测试技术领域,具体涉及一种基于气动探针正反向连续运动的定常流场快速测量方法。
背景技术
标准的气动探针测量***由气动探针、测压软管和压力传感器(压力扫描阀)等组成,一般用于稳态定常流场气动参数的测量。由于气动探针和测压软管形成的管腔效应会造成连接在测压软管末端的压力传感器测量到的变化压力信号发生滞后,因此在使用气动探针测量不同空间位置的流场参数时,通常采用离散点测量方式,即气动探针到达测量点位后,需等待一定时间,待测压管路压力达到平衡后再进行测量。然而当测量所需的空间测点数较多或测压管路压力的稳定时间较长时,这种测量方式会大大增加试验测量时间。
当前,亟需发展一种基于气动探针正反向连续运动的定常流场快速测量方法。
发明内容
本发明所要解决的技术问题是提供一种基于气动探针正反向连续运动的定常流场快速测量方法。
本发明的基于气动探针正反向连续运动的定常流场快速测量方法由气动探针正向和反向连续运动获得两组压力测量数据进行无相位偏移的低通滤波处理。通过设定阶离散***压力递推公式的传递函数的系数,经过修正后可获取两组修正后的压力测量数据,两组修正后的压力测量数据利用连续同步采集的气动探针位置信息插值能够获取相同的空间位置处的修正后的压力测量数据。以相同的空间位置的正向修正后的压力测量数据和反向修正后的压力测量数据的平均值作为气动探针的最终的修正后的压力测量数据。
本发明的基于气动探针正反向连续运动的定常流场快速测量方法,包括以下步骤:
S10.如果已经获得定常流场条件下的气动探针管路***压力传递的动态特性参数,转入步骤S100,否则转入步骤S20;
S20.在风洞的位移机构上安装气动探针,气动探针连接测压软管,测压软管末端连接压力信号测量装置;风洞启动,待流场稳定后,位移机构沿正向从静止加速到匀速速度再减速到静止,使气动探针从起点连续运动至终点;随后,位移机构沿反向从静止加速到与正向相同的匀速速度再减速到静止,使气动探针从终点运动至起点;运动过程中,数采***以相同的采样率同步采集气动探针位置数据和测量压力数据;
S30.对测量压力数据进行无相位偏移的低通滤波处理;
S60.根据公式(1)计算获取修正后的压力数据;
S70.根据同步采集的气动探针位置数据,对修正后的压力数据进行插值处理,分别获得相同位置处正向运动修正后压力数据和反向运动修正后压力数据的空间分布;
S80.计算相同的空间位置的正向修正后压力数据和反向修正后压力数据偏差以及数据序列y的极值点处压力数值和对应的修正后压力数值偏差,计算所有数据偏差的平方和;
S90.通过迭代优化的方法不断调整步骤S50设定的动态特性参数、,重复步骤S50~S80,直至所有数据偏差的平方和达到最小;以相同的空间位置的正向修正后压力数据和反向修正后压力数据的平均值作为气动探针的最终的修正后压力数据;转入步骤S110;
S100.在风洞的位移机构上安装气动探针,气动探针连接测压软管,测压软管末端连接压力信号测量装置;风洞启动,待流场稳定后,位移机构在按照预先设置的路线使气动探针由起点从静止加速到匀速速度再减速到静止,连续运动至终点,同时数采***以气动探针管路***的动态传递特性所标识的采样率同步采集探针位置和测量压力数据;对测量压力数据进行无相位偏移的低通滤波处理;由公式(1)获取修正后的压力数据;
S110.输出位置坐标和对应的修正后压力数据。
进一步地,所述的步骤S20和步骤S100中压力信号测量装置的固有频率大于等于100Hz。
进一步地,所述的步骤S20和步骤S100中数采***的采样率范围为10Hz~100Hz。
简而言之,本发明的基于气动探针正反向连续运动的定常流场快速测量方法不需要提前进行探针管路***动态特性标定试验,而是直接由位移机构带动气动探针在流场中分别进行正向和反向连续运动,对探针正反向运动测量数据进行修正,通过迭代优化方法使探针正向和反向运动修正后压力数据的偏差最小化,最终获得气动探针连续运动测孔位置感受到的“真实”压力变化。
具体而言,本发明的基于气动探针正反向连续运动的定常流场快速测量方法由位移机构带动气动探针在被测流场中分别沿正向和反向连续运动,数采***同步采集探针位置数据和探针管路末端压力测量数据。为避免高频脉动噪声信号的影响,对测量压力数据进行无相位偏移的低通滤波处理。依据线性离散***理论建立修正压力数据序列的递推公式,并给定气动探针管路***动态特性参数的预估初值,对气动探针正反向连续运动测量数据进行修正。利用同步采集的探针位置信息对修正后的正反向压力数据在相同探针位置坐标处进行插值处理,获取相同位置坐标下的修正后的正反向压力数据的空间分布。计算两组修正后的压力数据的偏差,同时计算修正前压力极值点位置数据与修正后对应数据的偏差。利用迭代优化方法不断调整管路***动态特性参数,使所有数据偏差的平方和最小化。理想情况下两组修正后的测量压力数据在相同的空间位置应该具有相同的数值,而且,气动探针管路***一般是一个过阻尼***,气动探针管路***的输入信号应该经过输出信号的极值点,可以通过迭代优化的方法,不断优化气动探针管路***的动态特性参数,使得两组修正后的测量压力数据在若干个相同的空间位置处的压力数值偏差平方和以及输出信号的极值点处压力数值和相应的修正后的测量压力数据的偏差平方和达到最小。最后以优化迭代收敛后的相同的空间位置处探针正向和反向运动修正后的压力数据的平均值作为最终气动探针连续运动修正后的压力数据,即为气动探针测孔位置实际感受到的压力空间分布。
本发明的基于气动探针正反向连续运动的定常流场快速测量方法可以方便有效地适应气动探针的各种定常流场环境,且在定常流场分布特征变化不大时也可以只进行一个方向上的气动探针连续运动测量,利用近似的定常流场的气动探针管路***的动态传递特性的阶离散***压力递推公式获取修正后的测量压力数据。
本发明的基于气动探针正反向连续运动的定常流场快速测量方法能够显著提高气动探针的定常流场气动参数的测试效率,降低试验成本,同时该方法相对传统离散点测量方式无新增测量硬件设备,实施便捷且具有较强的环境适应性。
具体实施方式
下面结合实例对本发明进行详细说明。
实施例1
本实施例的气动探针在进行风洞试验前,已经获得近似的定常流场的气动探针管路***的动态传递特性。具体工作步骤如下:
a1.在风洞的位移机构上安装气动探针;
b1.风洞启动,流场稳定后,位移机构在定常流场中按照预先设置的运动路线匀速运动气动探针,并以气动探针管路***的动态传递特性所标识的采样率采集位置坐标和测量压力数据;
c1.对测量压力数据进行无相位偏移的低通滤波处理;
e1.输出位置坐标和最终的修正后的测量压力数据;
实施例2
本实施例的气动探针在进行风洞试验前,未获得近似的定常流场的气动探针管路***的动态传递特性。具体工作步骤如下:
a2.在风洞的位移机构上安装气动探针;风洞启动,流场稳定后,位移机构沿正向,从静止加速到匀速再减速到静止,将气动探针从起点运动至终点;随后,位移机构沿反向,从静止加速到匀速再减速到静止,将气动探针从终点运动至起点;运动过程中,压力采集测试***以相同的采样率同步采集空间位置和测量压力数据;
b2.对测量压力数据进行无相位偏移的低通滤波处理;
对气动探针连续运动测量得到的测量压力数据进行无相位偏移的低通滤波处理;
展开后:
e2.将修正后的测量压力数据进行插值处理,获得相同的空间位置的正向修正后的测量压力数据和反向修正后的测量压力数据;
f2.计算相同的空间位置的正向修正后的测量压力数据和反向修正后的测量压力数据的偏差的平方和,以及输出信号的极值点处压力数值和相应的修正后的压力数值的偏差的平方和,不断调整步骤c2的传递函数的系数、,重复步骤c2~e2,直至两种偏差的平方和达到最小或者处于给定的偏差范围;以相同的空间位置的正向修正后的测量压力数据和反向修正后的测量压力数据的平均值作为气动探针的最终的修正后的测量压力数据;
g2.输出位置坐标和最终的修正后的测量压力数据。
第二种,假定气动探针在开始运动的起点位置处于压力平衡状态,此时压力传感器的测值即为起点位置气动探针测实际感受的压力值;由于气动探针在压力平衡状态时的输入值等于输出值,即频率为0时,***的增益值为1,此时动态特性参数有以下约束关系:
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域。对于熟悉本领域的人员而言,在不脱离本发明原理的前提下,可容易地实现另外的改进和润饰,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节。
Claims (4)
1.基于气动探针正反向连续运动的定常流场快速测量方法,其特征在于,包括以下步骤:
S10.如果已经获得定常流场条件下的气动探针管路***压力传递的动态特性参数,转入步骤S100,否则转入步骤S20;
S20.在风洞的位移机构上安装气动探针,气动探针连接测压软管,测压软管末端连接压力信号测量装置;风洞启动,待流场稳定后,位移机构沿正向从静止加速到匀速速度再减速到静止,使气动探针从起点连续运动至终点;随后,位移机构沿反向从静止加速到与正向相同的匀速速度再减速到静止,使气动探针从终点运动至起点;运动过程中,数采***以相同的采样率同步采集气动探针位置数据和测量压力数据;
S30.对测量压力数据进行无相位偏移的低通滤波处理;
S60.根据公式(1)计算获取修正后的压力数据;
S70.根据同步采集的气动探针位置数据,对修正后的压力数据进行插值处理,分别获得相同位置处正向运动修正后压力数据和反向运动修正后压力数据的空间分布;
S80.计算相同的空间位置的正向修正后压力数据和反向修正后压力数据偏差以及数据序列y的极值点处压力数值和对应的修正后压力数值偏差,计算所有数据偏差的平方和;
S90.通过迭代优化的方法不断调整步骤S50设定的动态特性参数、,重复步骤S50~S80,直至所有数据偏差的平方和达到最小;以相同的空间位置的正向修正后压力数据和反向修正后压力数据的平均值作为气动探针的最终的修正后压力数据;转入步骤S110;
S100.在风洞的位移机构上安装气动探针,气动探针连接测压软管,测压软管末端连接压力信号测量装置;风洞启动,待流场稳定后,位移机构在按照预先设置的路线使气动探针由起点从静止加速到匀速速度再减速到静止,连续运动至终点,同时数采***以气动探针管路***的动态传递特性所标识的采样率同步采集探针位置和测量压力数据;对测量压力数据进行无相位偏移的低通滤波处理;由公式(1)获取修正后的压力数据;
S110.输出位置坐标和对应的修正后压力数据。
2.根据权利要求1所述的基于气动探针正反向连续运动的定常流场快速测量方法,其特征在于,所述的步骤S20和步骤S100中压力信号测量装置的固有频率大于等于100Hz。
3.根据权利要求1所述的基于气动探针正反向连续运动的定常流场快速测量方法,其特征在于,所述的步骤S20和步骤S100中数采***的采样率范围为10Hz~100Hz。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210029277.1A CN114088335B (zh) | 2022-01-12 | 2022-01-12 | 基于气动探针正反向连续运动的定常流场快速测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210029277.1A CN114088335B (zh) | 2022-01-12 | 2022-01-12 | 基于气动探针正反向连续运动的定常流场快速测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114088335A true CN114088335A (zh) | 2022-02-25 |
CN114088335B CN114088335B (zh) | 2022-04-08 |
Family
ID=80308610
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210029277.1A Active CN114088335B (zh) | 2022-01-12 | 2022-01-12 | 基于气动探针正反向连续运动的定常流场快速测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114088335B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115452313A (zh) * | 2022-11-14 | 2022-12-09 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种声爆试验探针角度敏感性快速标定方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5929331A (en) * | 1997-01-14 | 1999-07-27 | The Texas A&M University System | Multi-directional, three component velocity measurement pressure probe |
CN101655406A (zh) * | 2009-09-11 | 2010-02-24 | 中国铁道科学研究院机车车辆研究所 | 气体压力传感器零点标定方法及装置 |
DE202012103700U1 (de) * | 2012-09-26 | 2012-10-22 | Ming Lu | Messgerät zum Messen des Drucks und der Geschwindigkeit eines Wirbelströmungsfeldes |
CN104807508A (zh) * | 2015-04-13 | 2015-07-29 | 浙江大学 | 一种实验教学用带测压管显示的数显流量计及测量方法 |
CN107101798A (zh) * | 2017-05-12 | 2017-08-29 | 中国科学院工程热物理研究所 | 一种动态五孔探针 |
CN110441027A (zh) * | 2019-09-10 | 2019-11-12 | 中国航发沈阳发动机研究所 | 一种控制探针自动跟踪***零点漂移的修正方法 |
CN110793745A (zh) * | 2019-12-04 | 2020-02-14 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种超音速风洞流校测压软管保护装置 |
CN111498141A (zh) * | 2020-04-21 | 2020-08-07 | 中国人民解放军空军工程大学 | 一种基于微型探针实现气流角度实时监测的方法与装置 |
WO2021174681A1 (zh) * | 2020-03-06 | 2021-09-10 | 上海海事大学 | 一种复合型五孔压力-温度探针 |
-
2022
- 2022-01-12 CN CN202210029277.1A patent/CN114088335B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5929331A (en) * | 1997-01-14 | 1999-07-27 | The Texas A&M University System | Multi-directional, three component velocity measurement pressure probe |
CN101655406A (zh) * | 2009-09-11 | 2010-02-24 | 中国铁道科学研究院机车车辆研究所 | 气体压力传感器零点标定方法及装置 |
DE202012103700U1 (de) * | 2012-09-26 | 2012-10-22 | Ming Lu | Messgerät zum Messen des Drucks und der Geschwindigkeit eines Wirbelströmungsfeldes |
CN104807508A (zh) * | 2015-04-13 | 2015-07-29 | 浙江大学 | 一种实验教学用带测压管显示的数显流量计及测量方法 |
CN107101798A (zh) * | 2017-05-12 | 2017-08-29 | 中国科学院工程热物理研究所 | 一种动态五孔探针 |
CN110441027A (zh) * | 2019-09-10 | 2019-11-12 | 中国航发沈阳发动机研究所 | 一种控制探针自动跟踪***零点漂移的修正方法 |
CN110793745A (zh) * | 2019-12-04 | 2020-02-14 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种超音速风洞流校测压软管保护装置 |
WO2021174681A1 (zh) * | 2020-03-06 | 2021-09-10 | 上海海事大学 | 一种复合型五孔压力-温度探针 |
CN111498141A (zh) * | 2020-04-21 | 2020-08-07 | 中国人民解放军空军工程大学 | 一种基于微型探针实现气流角度实时监测的方法与装置 |
Non-Patent Citations (2)
Title |
---|
A. R. PORRO: "Pressure probe designs for dynamic pressure measurements in a supersonic flow field", 《ICIASF 2001 RECORD, 19TH INTERNATIONAL CONGRESS ON INSTRUMENTATION IN AEROSPACE SIMULATION FACILITIES 》 * |
陈峰 等: "基于探针管路动态修正的压气机动态总压测试", 《实验流体力学》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115452313A (zh) * | 2022-11-14 | 2022-12-09 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种声爆试验探针角度敏感性快速标定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114088335B (zh) | 2022-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110553786B (zh) | 压力传感器的补偿方法和*** | |
CN114088335B (zh) | 基于气动探针正反向连续运动的定常流场快速测量方法 | |
CN107271081B (zh) | 基于两阶段最小二乘拟合的硅压阻式压力变送器温度补偿方法及装置 | |
CN108362911B (zh) | 一种基于最小二乘法和插值法的流速仪标定方法 | |
CN105510706B (zh) | 一种高精度欠采样测频方法 | |
CN111624594B (zh) | 一种基于转换量测重构的组网雷达跟踪方法及*** | |
CN114526796A (zh) | 一种提高量程比的流量计仪表系数动态修正方法 | |
JP3670900B2 (ja) | 伝送器自動校正方法及び伝送器自動校正装置 | |
US20240178856A1 (en) | Propagation delay compensation and interpolation filter | |
CN114046960B (zh) | 基于预先动态标定的气动探针定常流场连续测试方法 | |
CN109683524B (zh) | 一种对无采样保持的各采样信号进行采样同步的处理方法 | |
CN111917413A (zh) | 一种ti-adc通道间时序偏差校准方法 | |
CN109655775B (zh) | 一种任意波形发生器幅度扫频多尺度校准方法和装置 | |
WO2018161377A1 (zh) | 用于蓝牙发射功率校准的方法及装置 | |
CN112461489B (zh) | 用于低压测量的电子扫描阀参考压力控制***及应用方法 | |
CN113765517A (zh) | 一种温度自适应校准的高精度多通道同步采集*** | |
CN108848447B (zh) | 一种采用未知节点修正的差分DV_Distance节点定位方法 | |
CN114553334B (zh) | 相控阵天线指向误差测量方法、***、终端及装置 | |
US20210054799A1 (en) | Gas flow rate measurement device and gas flow rate measurement method | |
CN115389090A (zh) | 基于参考端控制的压力测量***及测量方法 | |
CN114710383A (zh) | 一种利用ofdm信号对仪表射频通道频响进行校准的方法 | |
CN114691437A (zh) | 一种测试校正装置及方法 | |
CN110207730B (zh) | 一种电阻式位移传感器温度自补偿方法 | |
CN1272603C (zh) | 交叉敏感情况下多输入检测***的动态补偿方法 | |
CN108152527B (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 |