[go: up one dir, main page]

CN111896211B - A partition interpolation method for five-hole probes under large inflow angles - Google Patents

A partition interpolation method for five-hole probes under large inflow angles Download PDF

Info

Publication number
CN111896211B
CN111896211B CN202010744265.8A CN202010744265A CN111896211B CN 111896211 B CN111896211 B CN 111896211B CN 202010744265 A CN202010744265 A CN 202010744265A CN 111896211 B CN111896211 B CN 111896211B
Authority
CN
China
Prior art keywords
pressure
angle
hole
coefficient
interpolation
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
Application number
CN202010744265.8A
Other languages
Chinese (zh)
Other versions
CN111896211A (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.)
Dalian Maritime University
Original Assignee
Dalian Maritime University
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 Dalian Maritime University filed Critical Dalian Maritime University
Priority to CN202010744265.8A priority Critical patent/CN111896211B/en
Publication of CN111896211A publication Critical patent/CN111896211A/en
Application granted granted Critical
Publication of CN111896211B publication Critical patent/CN111896211B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/06Measuring arrangements specially adapted for aerodynamic testing

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明提供一种大来流角度下五孔探针的分区插值方法,包括获取内区插值文件和外区插值文件;获取五孔探针采集到的五孔压力值,根据五孔压力值确定插值分区模式,所述插值分区模式分为内区模式和外区模式;当插值分区模式为内区模式时,基于所述内区插值文件对所述五孔压力值进行处理,得到流场的总压和静压;当插值分区模式为外区模式时,基于外区插值文件对所述五孔压力值进行处理,得到流场的总压和静压;基于流场的总压和静压计算得到流场的流速。当探头来流角度过大时,探头部分区域的气流会发生附面层分离,处在分离区中的感压孔将失效,从而导致测量不准。采用本发明方法可以克服上述问题,将探头适应的来流角度范围由±30°提升至±45°。

Figure 202010744265

The invention provides a subarea interpolation method for a five-hole probe under a large inflow angle, which includes acquiring an inner-area interpolation file and an outer-area interpolation file; acquiring the five-hole pressure value collected by the five-hole probe, and determining according to the five-hole pressure value Interpolation partition mode, the interpolation partition mode is divided into inner region mode and outer region mode; when the interpolation partition mode is inner region mode, the pressure value of the five holes is processed based on the inner region interpolation file, and the flow field is obtained. Total pressure and static pressure; when the interpolation partition mode is the outer zone mode, the five-hole pressure values are processed based on the outer zone interpolation file to obtain the total pressure and static pressure of the flow field; the total pressure and static pressure based on the flow field Calculate the flow velocity of the flow field. When the incoming flow angle of the probe is too large, the air flow in the part of the probe will be separated from the boundary layer, and the pressure-sensing hole in the separation area will fail, resulting in inaccurate measurement. By adopting the method of the present invention, the above problems can be overcome, and the range of the incoming flow angle adapted to the probe is increased from ±30° to ±45°.

Figure 202010744265

Description

一种大来流角度下五孔探针的分区插值方法A partition interpolation method for five-hole probes under large inflow angles

技术领域technical field

本发明涉及风洞工程测量领域,具体而言,尤其涉及一种大来流角度下五孔探针的分区插值方法。The invention relates to the field of wind tunnel engineering measurement, in particular, to a subarea interpolation method of a five-hole probe under a large incoming flow angle.

背景技术Background technique

五孔探针作为一种流场测量工具,在平面叶栅风洞实验中受到广泛的应用。使用传统插值方法只能测量到来流±30°内的气流,当来流角度过大,探头部分区域的气流会发生附面层分离,处在分离区中的感压孔都将失效,从而导致测量不准。在平面叶栅的相关实验中的大迎角工况下,由于柵后的附面层分离严重,气流矢量的变化大。特别是柵后流场角区中的流动非常复杂,除马赫数的剧烈变化外,气流的角度变化幅度也非常大,这已经超过了传统插值方法的测量范围。针对这种情况,提出分区插值法。As a flow field measurement tool, the five-hole probe is widely used in the plane cascade wind tunnel experiment. Using the traditional interpolation method can only measure the air flow within ±30° of the incoming flow. When the incoming flow angle is too large, the air flow in the part of the probe will be separated from the boundary layer, and the pressure-sensing holes in the separation area will fail, resulting in Inaccurate measurement. Under the condition of high angle of attack in the related experiments of the plane cascade, due to the severe separation of the boundary layer behind the cascade, the change of the airflow vector is large. In particular, the flow in the angular region of the flow field behind the grid is very complex. In addition to the drastic change of Mach number, the angular change of the airflow is also very large, which has exceeded the measurement range of traditional interpolation methods. In view of this situation, a partition interpolation method is proposed.

发明内容SUMMARY OF THE INVENTION

根据上述提出的大迎角工况超出现有插值测量范围造成测量不准的技术问题,而提供一种大来流角度下五孔探针的分区插值方法,通过分区插值方法,扩大探头适应的来流角度。According to the technical problem of inaccurate measurement caused by the above-mentioned large angle of attack working conditions exceeding the existing interpolation measurement range, a partitioned interpolation method of five-hole probe under large incoming flow angle is provided. inflow angle.

本发明采用的技术手段如下:The technical means adopted in the present invention are as follows:

一种大来流角度下五孔探针的分区插值方法,包括:A partition interpolation method for a five-hole probe under a large inflow angle, comprising:

获取内区插值文件和外区插值文件;Get the inner area interpolation file and the outer area interpolation file;

获取五孔探针采集到的五孔压力值,根据所述五孔压力值确定插值分区模式,所述插值分区模式分为内区模式和外区模式;acquiring the five-hole pressure value collected by the five-hole probe, and determining an interpolation partition mode according to the five-hole pressure value, and the interpolation partition mode is divided into an inner zone mode and an outer zone mode;

当插值分区模式为内区模式时,基于所述内区插值文件对所述五孔压力值进行处理,得到测点的气流方向角、总压和静压;When the interpolation partition mode is the inner zone mode, the pressure values of the five holes are processed based on the inner zone interpolation file to obtain the airflow direction angle, total pressure and static pressure of the measuring point;

当插值分区模式为外区模式时,基于所述外区插值文件对所述五孔压力值进行处理,得到测点的气流方向角、总压和静压;When the interpolation partition mode is the outer zone mode, the pressure values of the five holes are processed based on the outer zone interpolation file to obtain the airflow direction angle, total pressure and static pressure of the measuring point;

基于所述流场的总压和静压计算得到流场的流速。The flow velocity of the flow field is calculated based on the total pressure and static pressure of the flow field.

进一步地,所述内区插值文件包括内区方向校准系数曲线、内区总压系数校准曲线以及内区静压系数校准曲线;所述外区插值文件包括外区方向校准系数曲线、外区总压系数校准曲线以及外区静压系数校准曲线。Further, the inner zone interpolation file includes the inner zone direction calibration coefficient curve, the inner zone total pressure coefficient calibration curve and the inner zone static pressure coefficient calibration curve; the outer zone interpolation file includes the outer zone direction calibration coefficient curve, the outer zone total pressure coefficient curve Pressure coefficient calibration curve and outer zone static pressure coefficient calibration curve.

进一步地,所述根据所述五孔压力值确定插值分区模式,包括:Further, determining the interpolation partition mode according to the five-hole pressure value includes:

根据所述五孔压力值中压力值最大的孔的编号确定插值分区模式:若2号孔探针压力值最大时为内区模式,其余孔探针压力值最大时为外区模式。The interpolation partition mode is determined according to the number of the hole with the largest pressure value among the five hole pressure values: if the pressure value of the No. 2 hole probe is the largest, it is the inner zone mode, and the other holes are the outer zone mode when the probe pressure value is the largest.

进一步地,所述基于所述内区插值文件对所述五孔压力值进行处理,得到测点的气流方向角、总压和静压,包括:Further, the pressure value of the five holes is processed based on the interpolation file in the inner region to obtain the airflow direction angle, total pressure and static pressure of the measuring point, including:

根据所述五孔压力值计算偏航角校准系数和俯仰角校准系数;Calculate the yaw angle calibration coefficient and the pitch angle calibration coefficient according to the pressure value of the five holes;

基于所述偏航角校准系数和俯仰角校准系数通过线性插值方法查询内区方向校准系数曲线,从而获得偏航角与俯仰角;Based on the yaw angle calibration coefficient and the pitch angle calibration coefficient, the inner zone direction calibration coefficient curve is inquired by the linear interpolation method, so as to obtain the yaw angle and the pitch angle;

基于所述偏航角与俯仰角通过线性插值方法查询内区总压系数校准曲线,从而获得总压校准系数,根据所述五孔压力值和总压校准系数计算流场的总压;Based on the yaw angle and pitch angle, the calibration curve of the total pressure coefficient in the inner area is inquired by the linear interpolation method, so as to obtain the total pressure calibration coefficient, and the total pressure of the flow field is calculated according to the pressure value of the five holes and the total pressure calibration coefficient;

基于所述偏航角与俯仰角通过线性插值方法查询内区静压系数校准曲线从而获得静压校准系数,根据所述五孔压力值和静压校准系数计算流场的静压。Based on the yaw angle and pitch angle, the static pressure coefficient calibration curve of the inner region is queried by the linear interpolation method to obtain the static pressure calibration coefficient, and the static pressure of the flow field is calculated according to the five-hole pressure value and the static pressure calibration coefficient.

进一步地,所述基于所述外区插值文件对所述五孔压力值进行处理,得到测点的气流方向角、总压和静压,包括:Further, the pressure value of the five holes is processed based on the outer area interpolation file, and the airflow direction angle, total pressure and static pressure of the measuring point are obtained, including:

根据所述五孔压力值计算圆周角校准系数和锥角校准系数;Calculate the calibration coefficient of the circumference angle and the calibration coefficient of the cone angle according to the pressure value of the five holes;

基于所述圆周角校准系数和锥角校准系数通过线性插值方法查询外区方向校准系数曲线,获得圆周角与锥角;Based on the circumference angle calibration coefficient and the cone angle calibration coefficient, the outer zone direction calibration coefficient curve is inquired through a linear interpolation method to obtain the circumference angle and the cone angle;

基于所述圆周角与锥角通过线性插值方法查询外区总压系数校准曲线,从而获得总压校准系数,根据所述五孔压力值和总压校准系数计算流场的总压;Based on the circumference angle and the cone angle, the total pressure coefficient calibration curve of the outer area is inquired by the linear interpolation method, so as to obtain the total pressure calibration coefficient, and the total pressure of the flow field is calculated according to the five-hole pressure value and the total pressure calibration coefficient;

基于所述圆周角与锥角通过线性插值方法查询外区静压系数校准曲线从而获得静压校准系数,根据所述五孔压力值和静压校准系数计算流场的静压。Based on the circumference angle and the cone angle, the static pressure coefficient calibration curve of the outer area is queried by the linear interpolation method to obtain the static pressure calibration coefficient, and the static pressure of the flow field is calculated according to the five-hole pressure value and the static pressure calibration coefficient.

较现有技术相比,本发明具有以下优点:Compared with the prior art, the present invention has the following advantages:

1、本发明通过分区的方式,对来流角度处在内区的测点采用内区插值文件进行处理,对来流角度处在外区的测点采用外区插值文件进行测量,将探头适应的来流角度范围由±30°提升至±45°。1. The present invention adopts the interpolation file in the inner area to process the measuring points whose inflow angle is in the inner area, and uses the interpolation file in the outer area to measure the measuring points whose inflow angle is in the outer area. The incoming flow angle range is increased from ±30° to ±45°.

2、本发明克服了来流角度过大时,探头部分区域的气流会发生附面层分离,处在分离区中的感压孔将失效,从而导致测量不准的问题。2. The present invention overcomes the problem of inaccurate measurement caused by the separation of the boundary layer of the air flow in some areas of the probe when the incoming flow angle is too large, and the pressure-sensing holes in the separation area will fail.

基于上述理由本发明可在流场测量领域广泛推广。Based on the above reasons, the present invention can be widely promoted in the field of flow field measurement.

附图说明Description of drawings

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做以简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to illustrate the embodiments of the present invention or the technical solutions in the prior art more clearly, the following briefly introduces the accompanying drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description These are some embodiments of the present invention, and for those of ordinary skill in the art, other drawings can also be obtained from these drawings without any creative effort.

图1为本发明分区插值方法流程图。FIG. 1 is a flow chart of the method of partition interpolation according to the present invention.

图2为本发明插值分区模式示意图。FIG. 2 is a schematic diagram of an interpolation partition mode of the present invention.

图3为本发明线性插值运算流程图。FIG. 3 is a flow chart of the linear interpolation operation of the present invention.

图4为本发明内区方向校准系数曲线示意图。FIG. 4 is a schematic diagram of the calibration coefficient curve of the inner region direction according to the present invention.

图5为本发明内区总压系数校准曲线示意图。FIG. 5 is a schematic diagram of the calibration curve of the total pressure coefficient in the inner region of the present invention.

图6为本发明内区静压系数校准曲线示意图。FIG. 6 is a schematic diagram of the calibration curve of the static pressure coefficient in the inner zone of the present invention.

图7为本发明内区来流角度示意图。FIG. 7 is a schematic diagram of the inflow angle of the inner zone of the present invention.

图8为本发明外区来流角度示意图。FIG. 8 is a schematic diagram of the angle of incoming flow from the outer zone of the present invention.

具体实施方式Detailed ways

为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。In order to make those skilled in the art better understand the solutions of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only Embodiments are part of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。It should be noted that the terms "first", "second" and the like in the description and claims of the present invention and the above drawings are used to distinguish similar objects, and are not necessarily used to describe a specific sequence or sequence. It is to be understood that the data so used may be interchanged under appropriate circumstances such that the embodiments of the invention described herein can be practiced in sequences other than those illustrated or described herein. Furthermore, the terms "comprising" and "having" and any variations thereof, are intended to cover non-exclusive inclusion, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those expressly listed Rather, those steps or units may include other steps or units not expressly listed or inherent to these processes, methods, products or devices.

如图1所示,本发明提供了一种基于大来流角度下流场的五孔探针的分区插值法,包括:As shown in FIG. 1 , the present invention provides a subarea interpolation method of a five-hole probe based on a flow field under a large incoming flow angle, including:

获取内区插值文件和外区插值文件。其中,内区插值文件包括内区方向校准系数曲线、内区总压系数校准曲线以及内区静压系数校准曲线。内区插值文件在五孔探针出厂前,即已经由厂家标定完成,如图4到图6。外区插值文件包括外区方向校准系数曲线、外区总压系数校准曲线以及外区静压系数校准曲线。外区插值文件由使用者在使用前标定完成。Get the inner and outer interpolation files. The inner region interpolation file includes the inner region direction calibration coefficient curve, the inner region total pressure coefficient calibration curve, and the inner region static pressure coefficient calibration curve. The inner area interpolation file has been calibrated by the manufacturer before the five-hole probe leaves the factory, as shown in Figure 4 to Figure 6. The outer zone interpolation file includes the outer zone direction calibration coefficient curve, the outer zone total pressure coefficient calibration curve and the outer zone static pressure coefficient calibration curve. The outer area interpolation file is calibrated by the user before use.

获取五孔探针采集到的五孔压力值,根据五孔压力值确定插值分区模式,所述插值分区模式分为内区模式和外区模式。具体包括根据五孔压力值中压力值最大的孔的编号确定插值分区模式:若2号孔探针压力值最大时为内区模式,其余孔探针压力值最大为外区模式。区域划分如图2所示。在本发明中,外区模式包括1号孔探针压力值最大的1区模式、3号孔探针压力值最大的3区模式、4号孔探针压力值最大的4区模式、5号孔探针压力值最大的5区模式。且每一个外区模式下均包括独立的外区校准文件,即外区方向校准系数曲线、外区总压系数校准曲线以及外区静压系数校准曲线。The five-hole pressure value collected by the five-hole probe is acquired, and the interpolation partition mode is determined according to the five-hole pressure value, and the interpolation partition mode is divided into an inner zone mode and an outer zone mode. Specifically, the interpolation partition mode is determined according to the number of the hole with the largest pressure value among the five hole pressure values: if the pressure value of the No. 2 hole probe is the largest, it is the inner zone mode, and the rest of the holes with the highest pressure value are the outer zone mode. The area division is shown in Figure 2. In the present invention, the outer zone modes include zone 1 mode with the highest probe pressure value in hole 1, zone 3 mode with the highest probe pressure value in hole 3, zone 4 mode with the highest probe pressure value in hole 4, and zone 4 mode with the highest probe pressure value in hole 4. The 5-zone mode with the maximum value of the hole probe pressure. And each outer zone mode includes independent outer zone calibration files, namely outer zone direction calibration coefficient curve, outer zone total pressure coefficient calibration curve and outer zone static pressure coefficient calibration curve.

具体地,当来流的气流角较小时,探头上的每个感压孔均有效,都能正确地感受到压力,利用5个感压孔测得的压力数据可以对测点参数进行插值计算。但当气流方向大于探头角度(>25°)时,气流通过探针头部时会分离,使部分感压孔处于分离区中而失效。这时,选择对着来流方向的感压孔(压力值最大的感压孔)、紧挨着该孔两侧的感压孔以及2号孔获得的压力数据来对测点参数进行插值计算。这时五孔针变成4个四孔针。五孔探针在使用时通过测压管连接五个压力传感器,获得各孔压力数值。Specifically, when the airflow angle of the incoming flow is small, each pressure-sensing hole on the probe is effective, and the pressure can be sensed correctly. Using the pressure data measured by the five pressure-sensing holes, the parameters of the measuring point can be interpolated. . However, when the airflow direction is greater than the probe angle (>25°), the airflow will be separated when passing through the probe head, so that some of the pressure-sensing holes are in the separation area and fail. At this time, the pressure-sensing hole facing the incoming flow direction (the pressure-sensing hole with the largest pressure value), the pressure-sensing hole on both sides of the hole, and the pressure data obtained from the No. 2 hole are selected to interpolate the parameters of the measuring point. . At this time, the five-hole needle becomes four four-hole needles. When the five-hole probe is in use, five pressure sensors are connected through a pressure measuring tube to obtain the pressure value of each hole.

当来流角度较小时,中心感压孔P2具有最大压力。随着来流角度变大,探头部分区域发生附面层分离,导致部分感压孔失效时,压力值最大的孔变为新的中心孔。根据五个感压孔采集到的压力值,探针头部被分为五个区。定义P2孔所对的区域为内区,其它四个感压孔所对应的区域按编号排列,为外区。每个外区占据45°的扇形空间,如图2所示。When the incoming flow angle is small, the central pressure-sensing hole P2 has the maximum pressure. As the incoming flow angle increases, the boundary layer is separated in some areas of the probe, and when some pressure-sensing holes fail, the hole with the highest pressure value becomes a new central hole. According to the pressure values collected by the five pressure-sensing holes, the probe head is divided into five zones. The area corresponding to the P2 hole is defined as the inner area, and the areas corresponding to the other four pressure-sensing holes are arranged in numbers and are the outer area. Each outer zone occupies a 45° sector of space, as shown in Figure 2.

当插值分区模式为内区模式时,基于所述内区插值文件对所述五孔压力值进行处理,得到测点的气流方向角、总压和静压。如图3所示,包括:根据所述五孔压力值计算偏航角校准系数和俯仰角校准系数;基于所述偏航角校准系数和俯仰角校准系数通过线性插值方法查询内区方向校准系数曲线,从而获得偏航角与俯仰角,如图4所示,图中每一个点均对应一对已知的偏航角和俯仰角,其他非已知位置的偏航角和俯仰角可以通过其与已知偏航角与俯仰角的相对位置插值求得。基于所述偏航角与俯仰角通过线性插值方法查询内区总压系数校准曲线,从而获得总压校准系数,如图5所示,根据所述五孔压力值和总压校准系数计算流场的总压;基于所述偏航角与俯仰角通过线性插值方法查询内区静压系数校准曲线从而获得静压校准系数,如图6所示,根据所述五孔压力值和静压校准系数计算流场的静压。When the interpolation partition mode is the inner zone mode, the pressure values of the five holes are processed based on the inner zone interpolation file to obtain the airflow direction angle, total pressure and static pressure of the measuring point. As shown in FIG. 3 , it includes: calculating the yaw angle calibration coefficient and the pitch angle calibration coefficient according to the pressure value of the five holes; based on the yaw angle calibration coefficient and the pitch angle calibration coefficient, querying the inner zone direction calibration coefficient through a linear interpolation method curve to obtain the yaw and pitch angles, as shown in Figure 4, each point in the figure corresponds to a pair of known yaw and pitch angles, and the yaw and pitch angles of other non-known positions can be obtained through It is obtained by interpolating its relative position to the known yaw and pitch angles. Based on the yaw angle and pitch angle, the calibration curve of the total pressure coefficient in the inner area is queried by the linear interpolation method, so as to obtain the total pressure calibration coefficient. The total pressure of Calculate the static pressure of the flow field.

具体来说,内区模式下,内区插值文件定义如下参数:Specifically, in the inner region mode, the inner region interpolation file defines the following parameters:

Figure BDA0002607810890000051
Figure BDA0002607810890000051

Figure BDA0002607810890000052
Figure BDA0002607810890000052

Figure BDA0002607810890000053
Figure BDA0002607810890000053

Figure BDA0002607810890000054
Figure BDA0002607810890000054

Figure BDA0002607810890000055
Figure BDA0002607810890000055

其中,Kα为α方向校准系数,Kβ为β方向校准系数,α与β分别为偏航角与俯仰角,CPt为总压校准系数,CPs为静压校准系数,P1-P5为1-5孔的压力值,Pt为总压,Ps为静压,

Figure BDA0002607810890000061
为孔1、孔3、孔4、孔5压力平均值。探针头部的五个感压孔会分别采集到不同的压力。利用这五个压力值,通过公式(1)-(5)可以计算出对应来流角度下的Kα和Kβ,在通过线性插值方式,查询内区方向校准系数曲线,从而可以确定偏航角α与俯仰角β。具体地如图4所示,在方向校准曲线中找到点(Kα,Kβ)的位置,计算出点(Kα,Kβ)与其所在的最小四边形各边之间的距离,用线性插值法插值得到偏航角α与俯仰角β。Among them, K α is the α direction calibration coefficient, K β is the β direction calibration coefficient, α and β are the yaw angle and pitch angle respectively, C Pt is the total pressure calibration coefficient, C Ps is the static pressure calibration coefficient, P 1 -P 5 is the pressure value of holes 1-5, P t is the total pressure, P s is the static pressure,
Figure BDA0002607810890000061
It is the average pressure of hole 1, hole 3, hole 4 and hole 5. The five pressure-sensing holes on the probe head will collect different pressures respectively. Using these five pressure values, K α and K β at the corresponding inflow angle can be calculated by formulas (1)-(5), and the yaw can be determined by querying the calibration coefficient curve in the direction of the inner area by means of linear interpolation. angle α and pitch angle β. Specifically, as shown in Figure 4, find the position of the point (K α , K β ) in the direction calibration curve, calculate the distance between the point (K α , K β ) and each side of the smallest quadrilateral where it is located, and use linear interpolation The yaw angle α and the pitch angle β are obtained by normal interpolation.

在根据偏航角α与俯仰角β通过线性插值方式查询内区总压系数校准曲线和内区静压系数校准曲线,从而获取总压校准系数CPt和静压校准系数CPs。再根据公式(3)和(4)反推求出流程的总压Pt和静压PsAccording to the yaw angle α and the pitch angle β, the calibration curve of the total pressure coefficient in the inner area and the calibration curve of the static pressure coefficient in the inner area are inquired by linear interpolation, so as to obtain the total pressure calibration coefficient C Pt and the static pressure calibration coefficient C Ps . Then, according to formulas (3) and (4), the total pressure P t and static pressure P s of the process can be obtained by inversion.

当插值分区模式为外区模式时,基于所述外区插值文件对所述五孔压力值进行处理,得到流场的总压和静压。包括:根据所述五孔压力值计算圆周角校准系数和锥角校准系数;基于所述圆周角校准系数和锥角校准系数通过线性插值方法查询外区方向校准系数曲线,从而获得圆周角与锥角;基于所述圆周角与锥角通过线性插值方法查询外区总压系数校准曲线,从而获得总压校准系数,根据所述五孔压力值和总压校准系数计算流场的总压;基于所述圆周角与锥角通过线性插值方法查询外区静压系数校准曲线从而获得静压校准系数,根据所述五孔压力值和静压校准系数计算流场的静压。When the interpolation partition mode is the outer zone mode, the pressure values of the five holes are processed based on the outer zone interpolation file to obtain the total pressure and static pressure of the flow field. The method includes: calculating the calibration coefficient of the circumference angle and the calibration coefficient of the cone angle according to the pressure value of the five holes; based on the calibration coefficient of the circumference angle and the calibration coefficient of the cone angle, querying the calibration coefficient curve of the direction of the outer area through the linear interpolation method, so as to obtain the circumference angle and the cone angle Based on the circumference angle and the cone angle, the total pressure coefficient calibration curve of the outer area is queried by the linear interpolation method, so as to obtain the total pressure calibration coefficient, and the total pressure of the flow field is calculated according to the five-hole pressure value and the total pressure calibration coefficient; based on The circumference angle and the cone angle are obtained by querying the static pressure coefficient calibration curve of the outer area through a linear interpolation method to obtain the static pressure calibration coefficient, and the static pressure of the flow field is calculated according to the five-hole pressure value and the static pressure calibration coefficient.

具体来说,当探针测量大方向来流时,气流通过探针头部时会分离,使部分感压孔处于分离区中而失效。这时,选择对着来流方向的感压孔(压力值最大的感压孔)、紧挨着该孔两侧的感压孔以及2号孔获得的压力数据来计算锥角与圆周角系数及总压、静压系数,即此时五孔针相当于四孔针使用。Specifically, when the probe measures the incoming flow in a general direction, the air flow will be separated when passing through the probe head, so that part of the pressure-sensing hole is in the separation area and fails. At this time, select the pressure-sensitive hole facing the incoming flow direction (the pressure-sensitive hole with the largest pressure value), the pressure-sensitive hole on both sides of the hole, and the pressure data obtained from the No. 2 hole to calculate the taper angle and the circumference angle coefficient And the total pressure and static pressure coefficient, that is, the five-hole needle is equivalent to the four-hole needle at this time.

在外区模式下,对应1号孔探针压力值最大的1区模式、3号孔探针压力值最大的3区模式、4号孔探针压力值最大的4区模式、5号孔探针压力值最大的5区模式。分别预先标定了各模式下的外区方向校准系数曲线、外区总压系数校准曲线以及外区静压系数校准曲线。标定时,对于小来流角度区域,先固定1个α角度,然后旋转β角度,逐个扫描校准,直到扫过整个内区;对于大来流角度区域,先通过圆周角θ和锥角

Figure BDA0002607810890000062
在确定所要校准的布点,再换算得到对应的偏航角α和俯仰角β。校准的探针转角范围为α≤45°,β≤45°。小来流角度下,不同校准点间的偏航角α和俯仰角β相距5°;大来流角度下,不同校准点间的圆周角θ和锥角
Figure BDA0002607810890000071
相距5°。标定后的总压系数校准曲线示意图和静压系数校准曲线示意图形式与图5-图6中的曲线类似,此处不再赘述。In the outer zone mode, it corresponds to zone 1 mode with the highest probe pressure value in hole 1, zone 3 mode with the highest probe pressure value in hole 3, zone 4 mode with the highest probe pressure value in hole 4, and zone 4 mode with probe pressure value in hole 5. 5-zone mode with maximum pressure value. The calibration coefficient curve of the outer zone direction, the calibration curve of the total pressure coefficient of the outer zone and the calibration curve of the static pressure coefficient of the outer zone in each mode are pre-calibrated respectively. When calibrating, for the small inflow angle area, first fix an α angle, then rotate the β angle, and scan and calibrate one by one until the entire inner area is swept; for the large inflow angle area, first pass the circumference angle θ and cone angle.
Figure BDA0002607810890000062
After determining the points to be calibrated, the corresponding yaw angle α and pitch angle β are obtained by conversion. The calibrated probe rotation angle range is α≤45°, β≤45°. Under small incoming flow angle, the yaw angle α and pitch angle β between different calibration points are 5° apart; under large incoming flow angle, the circumference angle θ and cone angle between different calibration points
Figure BDA0002607810890000071
5° apart. The schematic diagram of the calibration curve of the total pressure coefficient and the schematic diagram of the calibration curve of the static pressure coefficient after calibration are similar to the curves in Figs.

在小来流角度(内区)下,来流角度使用偏航角α和俯仰角β来定义。迎着来流方向,右偏航为正。在大来流角度(外区)下,来流角度由圆周角θ和锥角φ来定义,其中φ始终为正,是气流方向与探针头部的夹角,θ定义为测点来流角度在XY平面内的投影与X轴之间的夹角,逆时针方向为正。如图7给出了内区模式下偏航角α和俯仰角β的示意图,图8给出了外区模式下圆周角θ和锥角φ的定义。各角度之间的换算关系为:At a small incoming flow angle (inner zone), the incoming flow angle is defined using the yaw angle α and the pitch angle β. Facing the direction of oncoming flow, the right yaw is positive. Under the large incoming flow angle (outer area), the incoming flow angle is defined by the circumferential angle θ and the cone angle φ, where φ is always positive and is the angle between the airflow direction and the probe head, and θ is defined as the incoming flow from the measuring point The angle between the projection of the angle in the XY plane and the X axis, the counterclockwise direction is positive. Figure 7 shows a schematic diagram of the yaw angle α and pitch angle β in the inner zone mode, and Figure 8 shows the definitions of the circumference angle θ and the cone angle φ in the outer zone mode. The conversion relationship between each angle is:

Figure BDA0002607810890000072
Figure BDA0002607810890000072

进一步地,外区模式下,外区插值文件定义如下参数:Further, in the outer zone mode, the outer zone interpolation file defines the following parameters:

Figure BDA0002607810890000073
Figure BDA0002607810890000073

Figure BDA0002607810890000074
Figure BDA0002607810890000074

Figure BDA0002607810890000075
Figure BDA0002607810890000075

Figure BDA0002607810890000081
Figure BDA0002607810890000081

上述公式中i表示受力最大的孔的编号,相应的公式代表该中工况下,五孔探针各系数的计算方法。Kθi,

Figure BDA0002607810890000082
为第i区方向系数,CPti,CPsi为第i区总压系数和静压系数。In the above formula, i represents the number of the hole with the largest force, and the corresponding formula represents the calculation method of each coefficient of the five-hole probe under the middle working condition. K θi ,
Figure BDA0002607810890000082
is the direction coefficient of the i-th zone, C Pti , C Psi are the total pressure coefficient and the static pressure coefficient of the i-th zone.

通过特定工况下的Kθi,和

Figure BDA0002607810890000083
再通过线性插值方式,查询对应的第i外区方向校准系数曲线,从而可以确定圆周角θ和锥角φ。在第i外区方向校准曲线中找到点
Figure BDA0002607810890000084
的位置,用线性插值法插值得到圆周角θ和锥角φ。再根据圆周角θ和锥角φ通过线性插值方式查询第i外区总压系数校准曲线和第i外区静压系数校准曲线,从而获取总压校准系数CPti和静压校准系数CPsi。再根据第i区的公式反推求出流程的总压Pt和静压Ps。By K θi under specific conditions, and
Figure BDA0002607810890000083
Then through the linear interpolation method, query the corresponding ith outer zone direction calibration coefficient curve, so that the circumference angle θ and the cone angle φ can be determined. Find the point in the ith outer zone orientation calibration curve
Figure BDA0002607810890000084
The position of , use linear interpolation to obtain the circumference angle θ and cone angle φ. Then according to the circumference angle θ and cone angle φ, the total pressure coefficient calibration curve of the i-th outer zone and the static pressure coefficient calibration curve of the i-th outer zone are inquired by linear interpolation, so as to obtain the total pressure calibration coefficient C Pti and the static pressure calibration coefficient C Psi . The total pressure P t and static pressure P s of the process can be obtained by inversion according to the formula in the i-th zone.

基于所述流场的总压和静压计算得到流场的流速。包括通过以下公式反推求得流场速度:The flow velocity of the flow field is calculated based on the total pressure and static pressure of the flow field. Including the back-calculation of the flow field velocity by the following formula:

Pt-Ps=(1/2)ρv2 P t -P s =(1/2)ρv 2

其中ρ表示空气密度,v表示流场速度。where ρ is the air density and v is the velocity of the flow field.

最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, but not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: The technical solutions described in the foregoing embodiments can still be modified, or some or all of the technical features thereof can be equivalently replaced; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the technical solutions of the embodiments of the present invention. scope.

Claims (5)

1. A partition interpolation method of a five-hole probe under a large incoming flow angle is characterized by comprising the following steps:
acquiring an inner region interpolation file and an outer region interpolation file;
acquiring a five-hole pressure value acquired by a five-hole probe, and determining an interpolation partition mode according to the five-hole pressure value, wherein the interpolation partition mode is divided into an inner partition mode and an outer partition mode;
when the interpolation partition mode is an inner region mode, processing the five-hole pressure value based on the inner region interpolation file to obtain an airflow direction angle, a total pressure and a static pressure of a measuring point, wherein the inner region interpolation file defines the following parameters:
Figure FDA0003566159620000011
Figure FDA0003566159620000012
Figure FDA0003566159620000013
Figure FDA0003566159620000014
Figure FDA0003566159620000015
wherein, KαIs a calibration coefficient of alpha direction, KβIs a beta direction calibration coefficient, alpha and beta are yaw angle and pitch angle, respectively, CPtFor total pressure calibration factor, CPsFor static pressure calibration factor, P1-P5At a pressure of 1-5 holes, PtIs the total pressure, PsIn order to be at a static pressure,
Figure FDA0003566159620000016
the average values of the pressures of the holes 1, 3, 4 and 5 are shown;
when the interpolation partition mode is an outer area mode, processing the five-hole pressure value based on the outer area interpolation file to obtain an airflow direction angle, total pressure and static pressure of a measuring point, wherein the outer area interpolation file defines the following parameters:
Figure FDA0003566159620000017
Figure FDA0003566159620000021
Figure FDA0003566159620000022
Figure FDA0003566159620000023
in the above formula, i represents the number of the hole with the largest stress, the corresponding formula represents the calculation method of each coefficient of the five-hole probe under the working condition, and Kθi,
Figure FDA0003566159620000024
Is the i-th zone directional coefficient, CPti,CPsiThe total pressure coefficient and the static pressure coefficient of the ith area are shown;
and calculating the flow velocity of the flow field based on the total pressure and the static pressure of the airflow at the measuring point.
2. The zonal interpolation method of a five-hole probe at a high inflow angle of claim 1, wherein the inner zone interpolation file comprises an inner zone direction calibration coefficient curve, an inner zone total pressure coefficient calibration curve, and an inner zone static pressure coefficient calibration curve; the outer area interpolation file comprises an outer area direction calibration coefficient curve, an outer area total pressure coefficient calibration curve and an outer area static pressure coefficient calibration curve.
3. The partition interpolation method for the five-hole probe under the large incoming flow angle according to claim 1 or 2, wherein the determining of the interpolation partition mode according to the five-hole pressure value comprises:
determining an interpolation partition mode according to the serial number of the hole with the maximum pressure value in the five-hole pressure values: if be the district mode when No. 2 hole probe pressure values are the biggest, the district mode is outer to the biggest to the probe pressure value in the remaining hole.
4. The method for performing zone interpolation on a five-hole probe under a large incoming flow angle according to claim 3, wherein the step of processing the five-hole pressure value based on the inner zone interpolation file to obtain the airflow direction angle, the total pressure and the static pressure of the measuring point comprises the following steps:
calculating a yaw angle calibration coefficient and a pitch angle calibration coefficient according to the five-hole pressure value;
inquiring an inner area direction calibration coefficient curve by a linear interpolation method based on the yaw angle calibration coefficient and the pitch angle calibration coefficient so as to obtain a yaw angle and a pitch angle;
inquiring an inner area total pressure coefficient calibration curve by a linear interpolation method based on the yaw angle and the pitch angle so as to obtain a total pressure calibration coefficient, and calculating the total pressure of the flow field according to the five-hole pressure value and the total pressure calibration coefficient;
and inquiring an inner region static pressure coefficient calibration curve through a linear interpolation method based on the yaw angle and the pitch angle so as to obtain a static pressure calibration coefficient, and calculating the static pressure of the flow field according to the five-hole pressure value and the static pressure calibration coefficient.
5. The zonal interpolation method of the five-hole probe under the incoming flow angle of claim 3, wherein the processing of the five-hole pressure value based on the outer zone interpolation file to obtain the airflow direction angle, total pressure and static pressure of the measurement point comprises:
calculating a circumferential angle calibration coefficient and a cone angle calibration coefficient according to the five-hole pressure value;
inquiring an outer region direction calibration coefficient curve by a linear interpolation method based on the circumferential angle calibration coefficient and the cone angle calibration coefficient so as to obtain a circumferential angle and a cone angle;
inquiring an outer region total pressure coefficient calibration curve by a linear interpolation method based on the circumferential angle and the cone angle so as to obtain a total pressure calibration coefficient, and calculating the total pressure of the flow field according to the five-hole pressure value and the total pressure calibration coefficient;
and inquiring an outer zone static pressure coefficient calibration curve by a linear interpolation method based on the circumferential angle and the cone angle so as to obtain a static pressure calibration coefficient, and calculating the static pressure of the flow field according to the five-hole pressure value and the static pressure calibration coefficient.
CN202010744265.8A 2020-07-29 2020-07-29 A partition interpolation method for five-hole probes under large inflow angles Active CN111896211B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010744265.8A CN111896211B (en) 2020-07-29 2020-07-29 A partition interpolation method for five-hole probes under large inflow angles

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010744265.8A CN111896211B (en) 2020-07-29 2020-07-29 A partition interpolation method for five-hole probes under large inflow angles

Publications (2)

Publication Number Publication Date
CN111896211A CN111896211A (en) 2020-11-06
CN111896211B true CN111896211B (en) 2022-06-21

Family

ID=73182517

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010744265.8A Active CN111896211B (en) 2020-07-29 2020-07-29 A partition interpolation method for five-hole probes under large inflow angles

Country Status (1)

Country Link
CN (1) CN111896211B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113155401B (en) * 2021-02-24 2022-12-20 中国空气动力研究与发展中心空天技术研究所 Rotatable linear flow measuring device and measuring method
CN114323536B (en) * 2021-12-30 2023-05-26 大连海事大学 Interpolation method for improving measurement accuracy of five-hole probe
CN115615612B (en) * 2022-09-28 2025-05-13 哈尔滨工业大学 A supersonic probe and calibration method thereof

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102636662A (en) * 2012-04-20 2012-08-15 南京航空航天大学 Full-flow direction flow velocity measurement probe and measurement method thereof
CN104713693A (en) * 2014-12-15 2015-06-17 中国燃气涡轮研究院 Pressure-leading type supersonic velocity five-hole probe with orthogonal auxiliary holes
CN106096223A (en) * 2016-05-10 2016-11-09 中国科学院工程热物理研究所 A kind of five-hole probe data processing method
CN108562420A (en) * 2017-12-26 2018-09-21 北京航天益森风洞工程技术有限公司 A kind of automatic calibration system and method for porous probe
CN110514391A (en) * 2019-08-12 2019-11-29 北京航空航天大学 A Method for Uncertainty Evaluation of Three-dimensional Flow Field Measurement by Five-hole Pressure Probe
CN111257591A (en) * 2020-03-31 2020-06-09 中国舰船研究设计中心 Method and device for measuring wind speed and direction of airflow field with seven-hole probe

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102636662A (en) * 2012-04-20 2012-08-15 南京航空航天大学 Full-flow direction flow velocity measurement probe and measurement method thereof
CN104713693A (en) * 2014-12-15 2015-06-17 中国燃气涡轮研究院 Pressure-leading type supersonic velocity five-hole probe with orthogonal auxiliary holes
CN106096223A (en) * 2016-05-10 2016-11-09 中国科学院工程热物理研究所 A kind of five-hole probe data processing method
CN108562420A (en) * 2017-12-26 2018-09-21 北京航天益森风洞工程技术有限公司 A kind of automatic calibration system and method for porous probe
CN110514391A (en) * 2019-08-12 2019-11-29 北京航空航天大学 A Method for Uncertainty Evaluation of Three-dimensional Flow Field Measurement by Five-hole Pressure Probe
CN111257591A (en) * 2020-03-31 2020-06-09 中国舰船研究设计中心 Method and device for measuring wind speed and direction of airflow field with seven-hole probe

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《1种副孔正交型超声速5孔探针的设计与应用》;张有 等;《航空发动机》;20181031;第44卷(第5期);第65-72页 *
《可调导叶尾缘吹气对离心压缩机动静干涉与叶轮气动载荷影响的研究》;辛建池;《中国优秀博硕士学位论文全文数据库(博士) 工程科技II辑》;20200115(第1期);第21-24页 *
《非对向法和混合编程在叶栅流场测量中的应用》;陈浮 等;《哈尔滨工业大学学报》;20031231;第35卷(第12期);第1426-1428页 *

Also Published As

Publication number Publication date
CN111896211A (en) 2020-11-06

Similar Documents

Publication Publication Date Title
CN111896211B (en) A partition interpolation method for five-hole probes under large inflow angles
WO2021174681A1 (en) Composite five-hole pressure-temperature probe
CN101046422B (en) Apparatus and method for detecting blockage of impulse piping
CN111551215A (en) A composite pressure-temperature probe and its air velocity calculation method
CN113405805B (en) Aero-engine inlet flow obtaining method considering boundary layer of air inlet channel
CN113916486B (en) Low-speed wind tunnel rapid pressure measurement and control method based on multi-parameter pressure regression algorithm
CN211234909U (en) Compound five-hole pressure-temperature probe
CN102112850A (en) Combined measurement instrument for flow field pressure and velocity and method thereof
CN212254182U (en) Composite pressure-temperature probe
WO2012028007A1 (en) Thickness measuring method for wafer marginal membrane
CN111487440A (en) Calibration method of five-hole probe
CN108593964B (en) A kind of calibration method of seven-hole airflow probe
CN110618293B (en) Calibration method of cross hot-wire anemometer
CN113405763B (en) Measuring point arrangement method of '3 +4+ 5' rice-shaped harrow flow meter
CN114459554A (en) Method for improving instantaneous flow metering precision of ultrasonic gas meter based on pressure data
CN115435930A (en) Three-dimensional full-parameter high-frequency probe for measuring interstage
CN107907250A (en) A kind of temperature-compensation method and device of silicon on sapphire pressure sensor
CN111089668B (en) A device and method for accurately measuring compressor stage temperature rise efficiency
CN212082824U (en) A probe for measuring full parameters of transonic three-dimensional steady flow field
CN212082771U (en) A full-parameter probe for measuring the two-dimensional steady flow field of the outer wall of the wheel hub
CN116182997A (en) Method for acquiring internal airflow speed of compressor and testing device
CN113125800B (en) Wind speed and direction measuring method based on pitot tube
CN111498141B (en) A method and device for real-time monitoring of airflow angle based on micro-probe
CN112964368B (en) Turbine blade radiation temperature measurement correction method
CN117214456B (en) An integrated wind speed and direction measuring instrument

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