CN111257591A - 一种七孔探针气流场风速风向测量方法及装置 - Google Patents

一种七孔探针气流场风速风向测量方法及装置 Download PDF

Info

Publication number
CN111257591A
CN111257591A CN202010242738.4A CN202010242738A CN111257591A CN 111257591 A CN111257591 A CN 111257591A CN 202010242738 A CN202010242738 A CN 202010242738A CN 111257591 A CN111257591 A CN 111257591A
Authority
CN
China
Prior art keywords
wind speed
partition
determining
airflow field
matrix
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.)
Pending
Application number
CN202010242738.4A
Other languages
English (en)
Inventor
朱传焕
唐君
张恒萍
郎波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing HKC Optoelectronics Technology Co Ltd
China Ship Development and Design Centre
Original Assignee
Chongqing HKC Optoelectronics Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing HKC Optoelectronics Technology Co Ltd filed Critical Chongqing HKC Optoelectronics Technology Co Ltd
Priority to CN202010242738.4A priority Critical patent/CN111257591A/zh
Publication of CN111257591A publication Critical patent/CN111257591A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/08—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring variation of an electric variable directly affected by the flow, e.g. by using dynamo-electric effect
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P13/00—Indicating or recording presence, absence, or direction, of movement

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Abstract

本发明公开了一种七孔探针气流场风速风向测量方法及装置，属于气流场测试领域，其中，方法的实现包括：获取七孔探针输出的七个孔的瞬时压力值；确认气流场来流对应的七孔探针的分区，根据各孔的瞬时压力值及分区对应的校准系数得到俯仰角、方位角及和风速相关的总压系数和静压系数；由俯仰角和方位角得到风向角，由风速相关的总压系数和静压系数得到气流场的风速。本发明可以用于高效快速准确测量气流场风速风向，规范测量流程，提高工作效率，减小不确定度影响因数。

Description

一种七孔探针气流场风速风向测量方法及装置
技术领域
本发明属于气流场测试领域，更具体地，涉及一种七孔探针气流场风速风向测量方法及装置。
背景技术
多孔探针测量速度场的方式，根据孔的数目不同，又分为三孔、五孔。并且随着加工工艺的进步，探针直径也越来越小，目前七孔探针孔直径达到4mm，后续发展直径将到达2mm。
七孔探针是头部具有七个压力测量孔的一种方向探针，图1展示的为七孔探针的实物图，其中，图1中(a)为七孔探针实物图，(b)为探针输出孔示意图，(c)为从迎风面观察的七孔位置关系。如图1所示，探针在头部共7个测压孔，分别定义为1～7号孔，其中，中心孔定义为7号测压孔。使用时，根据探针头部7个测压孔的压力值，可以进行精确计算获得来流的风速、风向的信息。
但是由于七孔探针压力方向孔的孔数多，其测量范围逐渐增加，但是其校准和测量原理也变得十分复杂。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提出了一种七孔探针气流场风速风向测量方法及装置，由此解决由于七孔探针压力方向孔的孔数多，导致的测量范围逐渐增加，进而使得校准和测量原理也变得十分复杂的技术问题。
为实现上述目的，按照本发明的一个方面，提供了一种七孔探针气流场风速风向测量方法，包括：
获取七孔探针输出的七个孔的瞬时压力值；
确认气流场来流对应的七孔探针的分区，根据各孔的瞬时压力值及所述分区对应的校准系数得到俯仰角、方位角及和风速相关的总压系数和静压系数；
由所述俯仰角和所述方位角得到风向角，由所述风速相关的总压系数和静压系数得到气流场的风速。
优选地，由Pj＝max{P1、P2、P3、P4、P5、P6、P7}确认气流场来流对应的七孔探针的分区，其中，P1表示第一个孔的瞬时压力值，P2表示第二个孔的瞬时压力值，P3表示第三个孔的瞬时压力值，P4表示第四个孔的瞬时压力值，P5表示第五个孔的瞬时压力值，P6表示第六个孔的瞬时压力值，P7表示第七个孔的瞬时压力值，j＝1或2或3或4或5或6或7。
优选地，若j＝1，则迎风区在第一分区，由
Figure BDA0002433095650000021
得到Cθ1、
Figure BDA0002433095650000022
由Cθ1、
Figure BDA0002433095650000023
确定矩阵
Figure BDA0002433095650000024
由
Figure BDA0002433095650000025
Figure BDA0002433095650000026
Figure BDA0002433095650000027
确定θ1、
Figure BDA0002433095650000028
CO1和Cq1，
Figure BDA0002433095650000029
Figure BDA00024330956500000210
表示第一分区对应的校准系数矩阵，由
Figure BDA00024330956500000211
确定风向角α和β，由
Figure BDA00024330956500000212
确定气流场的风速V，ρ为空气密度，第一分区为外区。
优选地，若j＝2，则迎风区在第二分区，由
Figure BDA0002433095650000031
得到Cθ2、
Figure BDA0002433095650000032
由Cθ2、
Figure BDA0002433095650000033
确定矩阵
Figure BDA0002433095650000034
由
Figure BDA0002433095650000035
Figure BDA0002433095650000036
Figure BDA0002433095650000037
确定θ2、
Figure BDA0002433095650000038
CO2和Cq2，
Figure BDA0002433095650000039
Figure BDA00024330956500000310
表示第二分区对应的校准系数矩阵，由
Figure BDA00024330956500000311
确定风向角α和β，由
Figure BDA00024330956500000312
确定气流场的风速V，ρ为空气密度，第二分区为外区。
优选地，若j＝3，则迎风区在第三分区，由
Figure BDA00024330956500000313
得到Cθ3、
Figure BDA00024330956500000314
由Cθ3、
Figure BDA00024330956500000315
确定矩阵
Figure BDA00024330956500000316
由
Figure BDA00024330956500000317
Figure BDA00024330956500000318
Figure BDA00024330956500000319
确定θ3、
Figure BDA00024330956500000320
CO3和Cq3，
Figure BDA00024330956500000321
Figure BDA00024330956500000322
表示第三分区对应的校准系数矩阵，由
Figure BDA00024330956500000323
确定风向角α和β，由
Figure BDA0002433095650000041
确定气流场的风速V，ρ为空气密度，第三分区为外区。
优选地，若j＝4，则迎风区在第四分区，由
Figure BDA0002433095650000042
得到Cθ4、
Figure BDA0002433095650000043
由Cθ4、
Figure BDA0002433095650000044
确定矩阵
Figure BDA0002433095650000045
由
Figure BDA0002433095650000046
Figure BDA0002433095650000047
Figure BDA0002433095650000048
确定θ4、
Figure BDA0002433095650000049
CO4和Cq4，
Figure BDA00024330956500000410
Figure BDA00024330956500000411
表示第四分区对应的校准系数矩阵，由
Figure BDA00024330956500000412
确定风向角α和β，由
Figure BDA00024330956500000413
确定气流场的风速V，ρ为空气密度，第四分区为外区。
优选地，若j＝5，则迎风区在第五分区，由
Figure BDA00024330956500000414
得到Cθ5、
Figure BDA00024330956500000415
由Cθ5、
Figure BDA00024330956500000416
确定矩阵
Figure BDA00024330956500000417
由
Figure BDA00024330956500000418
Figure BDA0002433095650000051
Figure BDA0002433095650000052
确定θ5、
Figure BDA0002433095650000053
CO5和Cq5，
Figure BDA0002433095650000054
Figure BDA0002433095650000055
表示第五分区对应的校准系数矩阵，由
Figure BDA0002433095650000056
确定风向角α和β，由
Figure BDA0002433095650000057
确定气流场的风速V，ρ为空气密度，第五分区为外区。
优选地，若j＝6，则迎风区在第六分区，由
Figure BDA0002433095650000058
得到Cθ6、
Figure BDA0002433095650000059
由Cθ6、
Figure BDA00024330956500000510
确定矩阵
Figure BDA00024330956500000511
由
Figure BDA00024330956500000512
Figure BDA00024330956500000513
Figure BDA00024330956500000514
确定θ6、
Figure BDA00024330956500000515
CO6和Cq6，
Figure BDA00024330956500000516
Figure BDA00024330956500000517
表示第六分区对应的校准系数矩阵，由
Figure BDA00024330956500000518
确定风向角α和β，由
Figure BDA00024330956500000519
确定场的风速V，ρ为空气密度，第六分区为外区。
优选地，若j＝7，则迎风区在第七分区，由
Figure BDA0002433095650000061
得到Cα7、Cβ7，其中
Figure BDA0002433095650000062
Figure BDA0002433095650000063
为P1～P6的算数平均值，由Cα7、Cβ7确定矩阵
Figure BDA0002433095650000064
由θ7＝[C7][Kα]15×1，
Figure BDA0002433095650000065
Figure BDA0002433095650000066
确定θ7、
Figure BDA0002433095650000067
Co7和Cq7，[Kα]15×1、[Kβ]15×1、
Figure BDA0002433095650000068
表示第七分区对应的校准系数矩阵，由
Figure BDA0002433095650000069
确定风向角α和β，由
Figure BDA00024330956500000610
确定气流场的风速V，ρ为空气密度，第七分区为内区。
优选地，所述气孔探针为4mm气孔探针。
按照本发明的另一个方面，提供了一种七孔探针气流场风速风向测量装置，包括：
传感器模块，用于将气孔探针的七个孔的瞬时压力值转换为对应的七个直流电压信号；
压力信号采集模块，用于将各所述直流电压信号转换为对应的七个数字信号，其中，每个数字信号对应一个孔的瞬时压力值；
数据处理模块，用于确认气流场来流对应的七孔探针的分区，根据各孔的瞬时压力值及所述分区对应的校准系数得到俯仰角、方位角及和风速相关的总压系数和静压系数；由所述俯仰角和所述方位角得到风向角，由所述风速相关的总压系数和静压系数得到气流场的风速。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：
气流场测试技术经过了多年的发展，产生了多种可用于流场测试的技术，然而综合比较而言，工程应用最为广泛的依然是七孔探针。4mm七孔探针是一种高性能的流场测试仪器，可以得到流场中三维速度及压力信息，最大可测流动角可达80°，为复杂气流场的测量提供了有力的实验手段。通过测量七个孔的压力值，能够精确地确定空间三维气流场的速度大小和方向、总压、静压等参数。与热线风速仪、激光多普勒测速仪相比，4mm七孔探针测量风速风向具有设备费用低、使用方便、不易损坏、精度高等优点，测量方法科学，可提高测量效率和测量结果可靠性，可用于高效快速准确测量气流场风速风向，规范测量流程，提高工作效率，减小不确定度影响因数。
附图说明
图1是本发明实施例提供的一种七孔探针模型示意图，其中，(a)为七孔探针实物图，(b)为探针输出孔示意图，(c)为从迎风面观察的七孔位置关系；
图2是本发明实施例提供的一种小流动角(内区)与大流动角(外区)分区图；
图3是本发明实施例提供的一种小流动角(内区)与大流动角(外区)流动特性示意图，其中，(a)为小流动迎角，(b)为大流动迎角；
图4是本发明实施例提供的一种七孔探针气流场风速风向测量方法流程示意图；
图5是本发明实施例提供的一种七孔探针气流场风速风向测量装置的结构示意图。
图6是本发明实施例提供的一种七孔探针在坐标系中α、β、θ、
Figure BDA0002433095650000071
角度关系示意图，其中，(a)为内区，(b)为外区。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
在本发明实例中，“第一”、“第二”、“第三”等是用于区别不同的对象，而不必用于描述特定的顺序或先后次序。
尽管流体流场测量已步入复杂而昂贵的光、电设备测量时代，例如激光多普勒测速仪(Laser Doppler Velocimetry，LDV)，粒子图像测速仪(Particle ImageVelocimetry，PIV)、热线风速仪(Hot Wire Anemometer，HWA)等，但在某些测量场合下，特别是测量平均速度，多孔压力探头仍是主要选择之一。在三维测量情况下，一个经过仔细设计、精良制造、严格校准的七孔探针，其测速精度可超过三维热线，在可以忽略接触式测量对流场的干扰的场合(例如外场实物全尺寸测量)，往往优先考虑简单实用的七孔探针，来测量速度压力，而不选用调试操作复杂的三维激光测速仪。
七孔探针可以得到流场中三维速度及压力信息。其主要应用有：
1)钝体尾迹复杂流动测量；
2)前体涡、边条涡、大后掠角机翼的前缘涡；翼尖涡、后体涡的定量测量；
3)进气道内流测量、叶栅流场测量；
4)风洞流场品质校测。
在钝体尾流场的测量中，由于流动的复杂性，气流有较大的偏角。尽管五孔探针、三维热线也可以测量尾流场局部点的三维速度分量，但都不能测量大于45°流动角的气流(气流相对于探头体轴线)，特别是三维旋涡复杂流动时，局部气流偏角往往超过45°。因此采用五孔探针、三维热线进行测量时误差很大，并且这种误差由于无法甄别是否为系统误差而无法进行修正。而七孔探针可测量气流偏角为78°的大偏角流动，测试精度为1％，并且可以得到空间流场某点的总压、静压。为提高测量效率，利用多根七孔探头制成耙，可以同时测量空间多个点的气动参数。而利用数控坐标架也可以提高测量效率。
七孔探针能够实现大角度的测量与其构造和测量原理是分不开的。七孔探针头部的7个测压孔将测量区域划分为7个区域，如图2所示。同时，定义中心孔所在的为内区，边孔所在的区域为外区。其中，第七孔所在的流动空间成为内区，其它六个区成为外区，每个区各占60°的扇形空间。关于探针的孔位定义如图1所示，其中，图1中的孔顺序及标号还可以按照其它方式进行标号，图1中的孔标号方式不应理解为对本发明的唯一性限定，相应的，图2中的分区号也不应理解为对本发明的唯一性限定。
内区与外区的定义与气流绕探针头部的流动状态密切相关，如图3所示，其中，(a)为小流动迎角，(b)为大流动迎角。
当气流在与探针体轴呈小角度范围(不大于24°)时，探针头部七个测压孔对气流角度的测量均有贡献，因此，对于这种小角度流动，定义为内区流动；在较大角度绕流(大于24°)时，气流在探针头部出现流动分离现象，仅有迎风面的压力孔对角度测量有贡献，因此，定义边孔所在的区域为大角度外区流动。
七孔探针对于角度的测量依赖于由探针头部所得压力数据，压力数据处理后的无量纲参数即敏感系数，与气流角度有对应的函数关系，与探针内区与外区的流动特性有关。
其中，七孔探针压力校准系数矩阵通过实验室校准得到，在利用七孔探针测量风速风向时，测量过程步骤和校准系数矩阵数据利用是准确测量气流场风速风向的关键和难点，在本发明实施例中，优选通过4mm七孔探针测量气流场的压力值，利用自适应算法准确得到气流场的风速风向大小。
如图4所示是本发明实施例提供的一种七孔探针气流场风速风向测量方法流程示意图，具体包括以下步骤：
S1：获取七孔探针输出的七个孔的瞬时压力值；
S2：确认气流场来流对应的七孔探针的分区，根据各孔的瞬时压力值及所述分区对应的校准系数得到俯仰角、方位角及和风速相关的总压系数和静压系数；
S3：由俯仰角和方位角得到风向角，由风速相关的总压系数和静压系数得到气流场的风速。
以下对各步骤的具体实施方式进行详细说明：
1)七孔探针测试原理是通过测量七个孔的压力值，利用校准系数矩阵和特定测量模型解算出流场方向、总压和静压，最后利用总压和静压求解流场速度。
2)测量采集
将4mm七孔探针安装在支架上，调整位置，测量未知的气流场，七孔探针输出七个瞬时压力值(P1、P2、P3、P4、P5、P6、P7)至7个相应的压力传感器，压力传感器输出相应电信号至压力信号采集模块和数据处理模块。
3)按自适应算法进行综合计算：
数据处理模块进行综合数据处理，根据自适应算法确认气流场来流对应的七孔探针的分区并提取校准系数计算矩阵对应系数，通过分区按特定的测量模型计算风速、风向；
自适应算法如公式(1)所示，确认气流场来流对应的七孔探针的分区：
Pj＝max{P1、P2、P3、P4、P5、P6、P7} (1)
其中，j＝1或2或3或4或5或6或7。
Figure BDA0002433095650000101
符号说明：
Figure BDA0002433095650000102
Figure BDA0002433095650000111
如图6所示为七孔探针在坐标系中α、β、θ、
Figure BDA0002433095650000112
角度关系示意图，其中，(a)为内区，(b)为外区。
定义：
Figure BDA0002433095650000113
若j＝1，说明迎风区在1区：
利用P1、P2、P6、P7(此时压力孔3、4、5在背风区，不予采用)的压力数据带入公式(4)，即可计算此时的Cθ1、
Figure BDA0002433095650000114
Figure BDA0002433095650000115
由于此时风向在1区，所以此时使用1区的校准系数
Figure BDA0002433095650000116
Figure BDA0002433095650000117
根据公式(3)可获得θ1和
Figure BDA00024330956500001110
以及和风速相关的系数Co1和Cq1。
将θ1和
Figure BDA00024330956500001111
代入到公式(2)即可获得风向角α和β。
Figure BDA0002433095650000118
Figure BDA0002433095650000119
Figure BDA0002433095650000121
由公式(5)、(6)、(7)即可计算出来流风速V，其中，ρ表示空气密度，
未知气流场的风速V和风向角α和β即可准确记录显示。
若j＝2，说明迎风区在2区：
利用P1、P2、P3、P7的压力数据带入公式(8)，即可计算此时的Cθ2、
Figure BDA0002433095650000122
Figure BDA0002433095650000123
由于此时迎风区是2区，所以此时使用2区的校准系数
Figure BDA0002433095650000124
Figure BDA0002433095650000125
根据公式(3)即可获得θ2和
Figure BDA0002433095650000129
以及和风速相关的系数Co2和Cq2。
将θ2和
Figure BDA00024330956500001210
代入到公式(2)即可获得风向角α和β。
Figure BDA0002433095650000126
Figure BDA0002433095650000127
将P∝L、PoL代入公式(7)即可计算出来流风速V。
未知气流场风速V和风向角α和β即可准确记录显示。
若j＝3，说明迎风区是3区：
利用P2、P3、P4、P7的压力数据带入公式(11)，即可计算此时的Cθ3、
Figure BDA0002433095650000128
Figure BDA0002433095650000131
由于此时迎风区是3区，所以此时使用3区的校准系数
Figure BDA0002433095650000132
Figure BDA0002433095650000133
根据公式(3)即可获得θ3和
Figure BDA00024330956500001310
以及和风速相关的系数Co3和Cq3。
将θ3和
Figure BDA00024330956500001311
代入到公式(2)即可获得风向角α和β。
Figure BDA0002433095650000134
Figure BDA0002433095650000135
将P∝L、PoL代入公式(7)即可计算出来流风速V。
未知气流场风速V和风向角α和β即可准确记录显示。
若j＝4，说明迎风区是4区：
利用P3、P4、P5、P7的压力数据带入公式(14)，即可计算此时的Cθ4、
Figure BDA0002433095650000136
Figure BDA0002433095650000137
由于此时迎风区是4区，所以此时使用4区的校准系数
Figure BDA0002433095650000138
Figure BDA0002433095650000139
根据公式(3)可得到θ4和
Figure BDA00024330956500001312
以及和风速相关的系数Co4和Cq4。
将θ4和
Figure BDA00024330956500001313
代入到公式(2)即可获得风向角α和β。
Figure BDA0002433095650000141
Figure BDA0002433095650000142
将P∝L、PoL代入公式(7)即可计算出来流风速。
未知气流场风速V和风向角α和β即可准确记录显示。
若j＝5，说明迎风区是5区：
利用P4、P5、P6、P7的压力数据带入公式(17)，即可计算此时的Cθ5、
Figure BDA0002433095650000143
Figure BDA0002433095650000144
由于此时迎风区是5区，所以此时使用5区的校准系数
Figure BDA0002433095650000145
Figure BDA0002433095650000146
根据公式(3)即可获得θ5和
Figure BDA0002433095650000147
以及和风速相关的系数Co5和Cq5。
将θ5和
Figure BDA0002433095650000148
代入到公式(2)即可获得风向角α和β。
Figure BDA0002433095650000149
Figure BDA00024330956500001410
将P∝L、PoL代入公式(7)即可计算出来流风速V。
未知气流场风速V和风向角α和β即可准确记录显示。
若j＝6，说明迎风区是6区：
利用P1、P5、P6、P7的压力数据带入公式(20)，即可计算此时的Cθ6、
Figure BDA00024330956500001411
Figure BDA0002433095650000151
由于此时迎风区是6区，此时使用6区的校准系数
Figure BDA0002433095650000152
Figure BDA0002433095650000153
根据公式(3)即可获得θ6和
Figure BDA0002433095650000154
以及和风速相关的系数Co6和Cq6。
将θ6和
Figure BDA0002433095650000155
代入到公式(2)即可获得风向角α和β。
Figure BDA0002433095650000156
Figure BDA0002433095650000157
将P∝L、PoL代入公式(7)即可计算出来流风速V。
未知气流场风速V和风向角α和β即可准确记录显示。
若j＝7，说明迎风区是7区：
利用P1～P7的压力数据带入公式(23)、(24)，即可计算此时的Cα7、Cβ7；
Figure BDA0002433095650000158
Figure BDA0002433095650000159
其中，P1、P2、P3、P4、P5、P6、P7分别表示第1区至第7区的压力值，
Figure BDA00024330956500001510
为P1～P6的算数平均值。由公式(23)、(24)得到Cα7、Cβ7。
由于此时迎风区是7区，此时使用7区的校准系数[Kα]15×1、[Kβ]15×1、
Figure BDA00024330956500001511
具体风速风向根据以下公式计算：
Figure BDA0002433095650000161
θ7＝[C7][Kα]15×1，
Figure BDA0002433095650000162
Figure BDA0002433095650000163
由此，可得到θ7和
Figure BDA0002433095650000164
以及和风速相关的系数Co7和Cq7。
将θ7和
Figure BDA0002433095650000165
代入到公式(2)即可获得α和β。
Figure BDA0002433095650000166
Figure BDA0002433095650000167
将P∝L、PoL代入公式(7)即可计算出来流风速V。
未知气流场风速V和风向角α和β即可准确记录显示。
如图5所示，“七孔探针”接收7路风速、风向信号后输出至压力传感器模块，压力传感器模块将信号转换为直流电压信号，输出至压力信号采集模块，压力信号采集模块将电压信号采集后编码成数字信号，数据处理模块将数字信号进行处理、存储、显示等。
具体地，传感器模块，用于将气孔探针的七个孔的瞬时压力值转换为对应的七个直流电压信号；
压力信号采集模块，用于将各所述直流电压信号转换为对应的七个数字信号，其中，每个数字信号对应一个孔的瞬时压力值；
数据处理模块，用于确认气流场来流对应的七孔探针的分区，根据各孔的瞬时压力值及所述分区对应的校准系数得到俯仰角、方位角及和风速相关的总压系数和静压系数；由所述俯仰角和所述方位角得到风向角，由所述风速相关的总压系数和静压系数得到气流场的风速。
其中，各模块的具体实施方式可以参考方法实施例的描述，本发明实施例将不再复述。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种七孔探针气流场风速风向测量方法，其特征在于，包括：
获取七孔探针输出的七个孔的瞬时压力值；
确认气流场来流对应的七孔探针的分区，根据各孔的瞬时压力值及所述分区对应的校准系数得到俯仰角、方位角及和风速相关的总压系数和静压系数；
由所述俯仰角和所述方位角得到风向角，由所述风速相关的总压系数和静压系数得到气流场的风速。
2.根据权利要求1所述的方法，其特征在于，由Pj＝max{P1、P2、P3、P4、P5、P6、P7}确认气流场来流对应的七孔探针的分区，其中，P1表示第一个孔的瞬时压力值，P2表示第二个孔的瞬时压力值，P3表示第三个孔的瞬时压力值，P4表示第四个孔的瞬时压力值，P5表示第五个孔的瞬时压力值，P6表示第六个孔的瞬时压力值，P7表示第七个孔的瞬时压力值，j＝1或2或3或4或5或6或7。
3.根据权利要求2所述的方法，其特征在于，若j＝1，则迎风区在第一分区，由
Figure FDA0002433095640000011
得到Cθ1、
Figure FDA0002433095640000012
由Cθ1、
Figure FDA0002433095640000013
确定矩阵
Figure FDA0002433095640000014
Figure FDA0002433095640000015
由
Figure FDA0002433095640000016
Figure FDA0002433095640000017
确定θ1、
Figure FDA0002433095640000018
CO1和Cq1，
Figure FDA0002433095640000019
表示第一分区对应的校准系数矩阵，由
Figure FDA00024330956400000110
确定风向角α和β，由
Figure FDA0002433095640000021
确定气流场的风速V，ρ为空气密度，第一分区为外区。
4.根据权利要求2所述的方法，其特征在于，若j＝2，则迎风区在第二分区，由
Figure FDA0002433095640000022
得到Cθ2、
Figure FDA0002433095640000023
由Cθ2、
Figure FDA0002433095640000024
确定矩阵
Figure FDA0002433095640000025
Figure FDA0002433095640000026
由
Figure FDA0002433095640000027
Figure FDA0002433095640000028
确定θ2、
Figure FDA0002433095640000029
CO2和Cq2，
Figure FDA00024330956400000210
表示第二分区对应的校准系数矩阵，由
Figure FDA00024330956400000211
确定风向角α和β，由
Figure FDA00024330956400000212
确定气流场的风速V，ρ为空气密度，第二分区为外区。
5.根据权利要求2所述的方法，其特征在于，若j＝3，则迎风区在第三分区，由
Figure FDA00024330956400000213
得到Cθ3、
Figure FDA00024330956400000214
由Cθ3、
Figure FDA00024330956400000215
确定矩阵
Figure FDA00024330956400000216
Figure FDA00024330956400000217
由
Figure FDA00024330956400000218
Figure FDA0002433095640000031
确定θ3、
Figure FDA0002433095640000032
CO3和Cq3，
Figure FDA0002433095640000033
表示第三分区对应的校准系数矩阵，由
Figure FDA0002433095640000034
确定风向角α和β，由
Figure FDA0002433095640000035
确定气流场的风速V，ρ为空气密度，第三分区为外区。
6.根据权利要求2所述的方法，其特征在于，若j＝4，则迎风区在第四分区，由
Figure FDA0002433095640000036
得到Cθ4、
Figure FDA0002433095640000037
由Cθ4、
Figure FDA0002433095640000038
确定矩阵
Figure FDA0002433095640000039
Figure FDA00024330956400000310
由
Figure FDA00024330956400000311
Figure FDA00024330956400000312
确定θ4、
Figure FDA00024330956400000313
CO4和Cq4，
Figure FDA00024330956400000314
表示第四分区对应的校准系数矩阵，由
Figure FDA00024330956400000315
确定风向角α和β，由
Figure FDA00024330956400000316
确定气流场的风速V，ρ为空气密度，第四分区为外区。
7.根据权利要求2所述的方法，其特征在于，若j＝5，则迎风区在第五分区，由
Figure FDA0002433095640000041
得到Cθ5、
Figure FDA0002433095640000042
由Cθ5、
Figure FDA0002433095640000043
确定矩阵
Figure FDA0002433095640000044
Figure FDA0002433095640000045
由
Figure FDA00024330956400000422
Figure FDA0002433095640000047
确定θ5、
Figure FDA0002433095640000048
CO5和Cq5，
Figure FDA0002433095640000049
表示第五分区对应的校准系数矩阵，由
Figure FDA00024330956400000410
确定风向角α和β，由
Figure FDA00024330956400000411
确定气流场的风速V，ρ为空气密度，第五分区为外区。
8.根据权利要求2所述的方法，其特征在于，若j＝6，则迎风区在第六分区，由
Figure FDA00024330956400000412
得到Cθ6、
Figure FDA00024330956400000413
由Cθ6、
Figure FDA00024330956400000414
确定矩阵
Figure FDA00024330956400000415
Figure FDA00024330956400000416
由
Figure FDA00024330956400000417
Figure FDA00024330956400000418
确定θ6、
Figure FDA00024330956400000419
CO6和Cq6，
Figure FDA00024330956400000420
表示第六分区对应的校准系数矩阵，由
Figure FDA00024330956400000421
确定风向角α和β，由
Figure FDA0002433095640000051
确定气流场的风速V，ρ为空气密度，第六分区为外区。
9.根据权利要求2所述的方法，其特征在于，若j＝7，则迎风区在第七分区，由
Figure FDA0002433095640000052
得到Cα7、Cβ7，其中，
Figure FDA0002433095640000053
Figure FDA0002433095640000054
为P1～P6的算数平均值，由Cα7、Cβ7确定矩阵
Figure FDA0002433095640000055
Figure FDA0002433095640000056
由θ7＝[C7][Kα]15×1，
Figure FDA0002433095640000057
Figure FDA0002433095640000058
确定θ7、
Figure FDA0002433095640000059
Co7和Cq7，[Kα]15×1、[Kβ]15×1、
Figure FDA00024330956400000510
表示第七分区对应的校准系数矩阵，由
Figure FDA00024330956400000511
确定风向角α和β，由
Figure FDA00024330956400000512
确定气流场的风速V，ρ为空气密度，第七分区为内区。
10.一种七孔探针气流场风速风向测量装置，其特征在于，包括：
传感器模块，用于将气孔探针的七个孔的瞬时压力值转换为对应的七个直流电压信号；
压力信号采集模块，用于将各所述直流电压信号转换为对应的七个数字信号，其中，每个数字信号对应一个孔的瞬时压力值；
数据处理模块，用于确认气流场来流对应的七孔探针的分区，根据各孔的瞬时压力值及所述分区对应的校准系数得到俯仰角、方位角及和风速相关的总压系数和静压系数；由所述俯仰角和所述方位角得到风向角，由所述风速相关的总压系数和静压系数得到气流场的风速。
CN202010242738.4A 2020-03-31 2020-03-31 一种七孔探针气流场风速风向测量方法及装置 Pending CN111257591A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010242738.4A CN111257591A (zh) 2020-03-31 2020-03-31 一种七孔探针气流场风速风向测量方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010242738.4A CN111257591A (zh) 2020-03-31 2020-03-31 一种七孔探针气流场风速风向测量方法及装置

Publications (1)

Publication Number Publication Date
CN111257591A true CN111257591A (zh) 2020-06-09

Family

ID=70949941

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010242738.4A Pending CN111257591A (zh) 2020-03-31 2020-03-31 一种七孔探针气流场风速风向测量方法及装置

Country Status (1)

Country Link
CN (1) CN111257591A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856072A (zh) * 2020-07-29 2020-10-30 中国汽车工程研究院股份有限公司 一种气流速度计算方法、系统、设备及一种存储介质
CN111896211A (zh) * 2020-07-29 2020-11-06 大连海事大学 一种大来流角度下五孔探针的分区插值方法
CN112113740A (zh) * 2020-11-18 2020-12-22 中国空气动力研究与发展中心低速空气动力研究所 一种测压装置、测压方法及风洞局部气流偏角测量方法
CN113390601A (zh) * 2021-08-16 2021-09-14 中国空气动力研究与发展中心低速空气动力研究所 一种动态流场测量方法、系统及设备
CN113753262A (zh) * 2021-11-09 2021-12-07 中国空气动力研究与发展中心低速空气动力研究所 一种直升机平尾区域流场速度的测量装置及方法
CN114062710A (zh) * 2021-11-18 2022-02-18 西安西热锅炉环保工程有限公司 一种3d流速测定探头及测定方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014036701A1 (zh) * 2012-09-06 2014-03-13 Lu Ming 测量旋流场的压力和速度的工具
CN207703890U (zh) * 2017-12-26 2018-08-07 北京航天益森风洞工程技术有限公司 一种基于五孔探针的高速列车侧风数据测量系统
CN108593964A (zh) * 2018-03-22 2018-09-28 中国科学院工程热物理研究所 一种七孔气流探针的标定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014036701A1 (zh) * 2012-09-06 2014-03-13 Lu Ming 测量旋流场的压力和速度的工具
CN207703890U (zh) * 2017-12-26 2018-08-07 北京航天益森风洞工程技术有限公司 一种基于五孔探针的高速列车侧风数据测量系统
CN108593964A (zh) * 2018-03-22 2018-09-28 中国科学院工程热物理研究所 一种七孔气流探针的标定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
丁涛: "七孔探针实时测试系统及其可压缩流场测量问题研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856072A (zh) * 2020-07-29 2020-10-30 中国汽车工程研究院股份有限公司 一种气流速度计算方法、系统、设备及一种存储介质
CN111896211A (zh) * 2020-07-29 2020-11-06 大连海事大学 一种大来流角度下五孔探针的分区插值方法
CN111856072B (zh) * 2020-07-29 2021-10-15 中国汽车工程研究院股份有限公司 一种气流速度计算方法、系统、设备及一种存储介质
CN111896211B (zh) * 2020-07-29 2022-06-21 大连海事大学 一种大来流角度下五孔探针的分区插值方法
CN112113740A (zh) * 2020-11-18 2020-12-22 中国空气动力研究与发展中心低速空气动力研究所 一种测压装置、测压方法及风洞局部气流偏角测量方法
CN113390601A (zh) * 2021-08-16 2021-09-14 中国空气动力研究与发展中心低速空气动力研究所 一种动态流场测量方法、系统及设备
CN113753262A (zh) * 2021-11-09 2021-12-07 中国空气动力研究与发展中心低速空气动力研究所 一种直升机平尾区域流场速度的测量装置及方法
CN114062710A (zh) * 2021-11-18 2022-02-18 西安西热锅炉环保工程有限公司 一种3d流速测定探头及测定方法

Similar Documents

Publication Publication Date Title
CN111257591A (zh) 一种七孔探针气流场风速风向测量方法及装置
US11422057B2 (en) Dynamic five-hole probe
WO2021174681A1 (zh) 一种复合型五孔压力-温度探针
JP3234894B2 (ja) 非定常流速の測定方法および装置
EP0255056B1 (en) Method for measuring the speed of a gas flow
US7010970B2 (en) Embedded-sensor multi-hole probes
CN102112850A (zh) 流场压力和速度的组合测量工具和方法
CN212254182U (zh) 一种复合型压力-温度探针
CN211234909U (zh) 一种复合型五孔压力-温度探针
CN111551215A (zh) 一种复合型压力-温度探针及其气流速度计算方法
CN115435929B (zh) 一种高频总温、总压探针
CN105588703A (zh) 一种用于亚声速复杂三维流场测量的12孔全向探针
CN108195510A (zh) 一种热膜式剪应力传感器的热风洞校准方法
JP2010518370A (ja) ガス状流体の流れの速さおよび方向を測定する装置および方法
Haneda et al. Neural network-based airflow vector sensor using multiple MEMS differential pressure sensors
CN112945500A (zh) 一种跨声速风洞堵塞度影响的模拟测量装置及方法
CN111498141A (zh) 一种基于微型探针实现气流角度实时监测的方法与装置
CN218822484U (zh) 一种测量级间三维全参数高频探针
CN113125800B (zh) 基于皮托管的风速风向测量方法
WO2014036701A1 (zh) 测量旋流场的压力和速度的工具
Efremova et al. Models of Signals, Algorithms and Errors of Instrumental Channels of the Air Data System Based on the Vortex Method
Rupesh Ms et al. Experimental and Computational Evaluation of Five Hole Five Probe Flow Analyzer for Subsonic Wind Calibration
CN111487031A (zh) 基于微型扇形探针实现三维流场气流角度实时监测的装置与方法
Yao et al. Theoretical Calibration and Experimental Test of Hemispherical 7-hole Probe Under Three-dimensional Flow Field in Wind Tunnel
Chernoray et al. Improving the accuracy of multihole probe measurements in velocity gradients

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200609

RJ01 Rejection of invention patent application after publication