CN113405763B - 一种“3+4+5”米字耙流量计测点布置方法 - Google Patents

一种“3+4+5”米字耙流量计测点布置方法 Download PDF

Info

Publication number
CN113405763B
CN113405763B CN202110654623.0A CN202110654623A CN113405763B CN 113405763 B CN113405763 B CN 113405763B CN 202110654623 A CN202110654623 A CN 202110654623A CN 113405763 B CN113405763 B CN 113405763B
Authority
CN
China
Prior art keywords
measuring
area
bent
points
total pressure
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
CN202110654623.0A
Other languages
English (en)
Other versions
CN113405763A (zh
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.)
Institute of Aerospace Technology of China Aerodynamics Research and Development Center
Original Assignee
Institute of Aerospace Technology of China Aerodynamics Research and Development Center
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 Institute of Aerospace Technology of China Aerodynamics Research and Development Center filed Critical Institute of Aerospace Technology of China Aerodynamics Research and Development Center
Priority to CN202110654623.0A priority Critical patent/CN113405763B/zh
Publication of CN113405763A publication Critical patent/CN113405763A/zh
Application granted granted Critical
Publication of CN113405763B publication Critical patent/CN113405763B/zh
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
    • G01M9/06Measuring arrangements specially adapted for aerodynamic testing

Landscapes

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

Abstract

本发明专利公开了一种“3+4+5“米字耙流量计测点布置方法,具体涉及实验空气动力学技术领域,包括如下步骤:进气道椭圆截面出口布置有米字耙流量计,每根垂直测量排架的3个四等分点上等间距布置3个总压测点,每根倾斜测量排架的4个五等分点上等间距布置4个总压测点,每根水平测量排架的5个六等分点上等间距布置5个总压测点;米字耙流量计的中心布置1个总压测点,测量排架的中心线与椭圆截面出口的交点处分别布置1个静压测点。采用本发明技术方案解决了现有的椭圆截面流量计测点布置方法会导致近壁面测点测量结果不准确、中心测点难以布置的问题,可用于进气道椭圆截面出口流量的精确测量。

Description

一种“3+4+5“米字耙流量计测点布置方法
技术领域
本发明涉及实验空气动力学技术领域,特别涉及一种“3+4+5“米字耙流量计测点布置方法。
背景技术
飞行器的进气道采用椭圆截面出口设计,可使飞行器结构紧致、设计灵活,同时可与前体及燃烧室容易一体化,因此在高超声速飞行器设计中使用越来越频繁。进气道的出口流量是与高超声速飞行器发动机匹配的重要指标,同时也是进气道设计校核、性能评估及数值计算方法检验的重要途径。
对于采用椭圆截面出口设计的进气道,其流量测量通常采用米字耙流量计,目前公开的米字耙流量计测点布置方法有两种:一是在进气道的出口处采用5×8测点等面积分布方法,将进气道出口椭圆截面分割成5×8个彼此互不重叠、互无间隙的区域单元,利用米字耙流量计在每个区域单元的面积中心布置总压测点;二是在进气道的出口处采用41测点等距分布,利用米字耙流量计的8个耙位将椭圆截面出口的面积八等分,8个耙位的中心交点布置1个皮托压测点,每个耙位上沿径向等距布置5个皮托压测点,米字耙流量计共设置41个皮托压测点。
上述两种测点布置方法主要存在以下不足之处:当进气道出口的椭圆截面尺度较小时:一是靠近壁面的测点易与壁面产生激波/边界层干扰,进而影响测量结果的准确性;二是椭圆截面短半轴布置的测压管容易发生气动干扰且易堵塞流道,特别是中心点的固定往往因几何干涉而难以加工与安装。
因此,米字耙流量计的测点布置是进气道椭圆截面出口流量测量的关键所在。
发明内容
本发明意在提供一种“3+4+5“米字耙流量计测点布置方法,解决了现有的测点布置方法会导致测量结果不准确、测点几何干涉、测点易于堵塞流道的问题。
为了达到上述目的,本发明的技术方案如下:一种“3+4+5“米字耙流量计测点布置方法,包括如下步骤:
S1、进气道椭圆截面出口布置有米字耙流量计,米字耙流量计上布置有与椭圆截面短半轴重合的2根垂直测量排架和与椭圆截面长半轴重合的2根水平测量排架,米字耙流量计上还布置有4根分别位于椭圆截面四个象限的倾斜测量排架;
S2、步骤S1的每根垂直测量排架的3个四等分点上等间距布置3个总压测点,每根倾斜测量排架的4个五等分点上等间距布置4个总压测点,每根水平测量排架的5个六等分点上等间距布置5个总压测点;
S3、步骤S1的每根垂直测量排架、倾斜测量排架和水平测量排架的相交处(即进气道椭圆截面出口的中心)布置有1个总压测点;
S4、步骤S1的每根垂直测量排架、倾斜测量排架和水平测量排架的中心线与进气道椭圆截面出口的交点处布置有1个静压测点。
进一步的,倾斜测量排架的中心线与进气道椭圆截面长半轴的夹角通过如下公式确定:
θ=arctan(b/a),0<θ<π/2。
进一步的,测点布置方法的数据处理方法包括如下步骤:
第一步、获得测点数据:
从米字耙流量计测量排架上获得32个测点的总压值,从米字耙流量计的中心测点获得1个总压值,从进气道的椭圆截面出口的静压测点获得8个壁面静压值;
第二步、对进气道的椭圆截面出口进行面积分割:
将椭圆截面出口的面积从周向和径向分割成40个面积相等、且互无间隔互不重叠的区域单元;
第三步、求各区域单元面积中心点的总压:
沿每根测量排架,插值得到该排架上等面积单元区域的面心处总压;
第四步、求各区域单元面积中心点的静压:
对8个壁面静压值进行算数平均得到米字耙流量计中心点的平均静压,由某根测量排架对应的壁面静压与中心点平均静压通过线性插值得到位于该测量排架上的各区域单元面积中心点的静压,从而求得40个区域单元面积中心点的静压;
第五步、求各区域单元上的流量:
根据40个区域单元面积中心点上的总压与静压,以及各区域单元的面积,计算出各区域单元面积中心点上的马赫数和流量气动函数,从而求得各区域单元面上的流量,各区域单元上的流量通过如下公式确定:
Figure BDA0003112129510000031
第六步,求椭圆出口截面的流量:
将40个区域单元面上的流量求和,从而得到椭圆出口截面的流量,椭圆出口截面的流量通过如下公式确定:
Figure BDA0003112129510000032
与现有技术相比,本方案的有益效果:
本方案根据进气道椭圆截面出口的尺寸对米字耙流量计测量排架上的总压测点合理布局,既可以避免流量计边缘测点与椭圆截面出口内壁面的激波边界层相互干扰问题,又可以避免测量排架上总压测点之间的气动干扰和几何干涉,为椭圆截面出口的流量计算和数据处理提供了科学依据。
附图说明
图1是本实施例中进气道椭圆截面出口处的测点分布图。
具体实施方式
下面通过具体实施方式对本发明作进一步详细的说明:
说明书附图中的附图标记包括:水平测量排架1、垂直测量排架2、倾斜测量排架3。
实施例
本实施例中相关符号及术语的说明如下表:
Figure BDA0003112129510000033
如附图1所示:一种“3+4+5“米字耙流量计测点布置方法,包括如下步骤:
S1、进气道椭圆截面出口布置有米字耙流量计,米字耙流量计上布置有与椭圆截面短半轴重合的2根垂直测量排架2和与椭圆截面长半轴重合的2根水平测量排架1,米字耙流量计上还布置有4根分别位于椭圆截面四个象限的倾斜测量排架3。
其中,倾斜测量排架3的中心线与进气道椭圆截面长半轴的夹角通过如下公式确定:
θ=arctan(b/a),0<θ<π/2。
S2、步骤S1的每根垂直测量排架2的3个四等分点上等间距布置3个总压测点,每根倾斜测量排架3的4个五等分点上等间距布置4个总压测点,每根水平测量排架1的5个六等分点上等间距布置5个总压测点。
S3、步骤S1的每根垂直测量排架2、倾斜测量排架3和水平测量排架1的相交处(即进气道椭圆截面出口的中心)布置有1个总压测点。
S4、步骤S1的每根垂直测量排架2、倾斜测量排架3和水平测量排架1的中心线与椭圆截面出口的交点处布置有1个静压测点。
每根垂直测量排架2上分别等间距布置3个总压测点,其总压测点分别位于垂直测量排架2中心线的4等分点上,即总压测点之间的间距为b/4;每根倾斜测量排架3上分别等间距布置的4个总压测点,总压测点分别位于倾斜测量排架3中心线的5等分点上,即总压测点之间的间距为1/5;每根水平测量排架1上分别等间距布置有5个总压测点,总压测点分别位于水平测量排架1中心线的6等分点上,即总压测点之间的间距为a/6;进气道椭圆截面出口的中心布置有1个总压测点,共计33个总压测点。
米字耙流量计各测量排架的中心线与椭圆截面出口的内壁面交点处分别布置1个静压测点,共计8个静压测点。
本方案的数据处理方法包括如下步骤:
第一步、获得测点数据:
从米字耙流量计测量排架上获得32个测点的总压值,从米字耙流量计的中心测点获得1个总压值,从椭圆截面出口壁面的静压测点获得8个壁面静压值。
第二步、对进气道的椭圆截面出口进行面积分割:
将椭圆截面出口的面积从周向和径向分割成40个面积相等、且互无间隔互不重叠的区域单元。
第三步、求各区域单元面积中心点的总压:
沿每根测量排架,插值得到该排架上等面积单元区域的面心处总压。
插值方法:使用Akima或其他插值方法;每根测量排架所在径向均有5~7个总压值(即两根水平测量排架1上各有5个总压测点值、米字耙流量计中心的1个总压测点值、1个壁面总压值(即壁面静压测点值),共7个测点值;同理可知,两根垂直测量排架2上各有3个总压测点值、米字耙流量计中心的1个总压测点值、1个壁面总压值(即壁面静压测点值),共5个测点值;四根倾斜测量排架3上各有4个总压测点值、米字耙流量计中心的1个总压测点值、1个壁面总压值(即壁面静压测点值),共6个测点值),对每一根测量排架的总压值进行插值,可以得到该排架上5个按等环面积分布的区域单元面积中心的总压值,从而求得40个单元面积中心上的总压。
第四步、求各区域单元面积中心点的静压:
对8个壁面测点的静压值进行算数平均得到米字耙流量计中心点的平均静压,由某根测量排架对应的壁面静压与中心点平均静压通过线性插值得到位于该测量排架上的各区域单元面积中心点的静压,从而求得40个区域单元面积中心点的静压。
第五步、求各区域单元上的流量:
根据40个区域单元面积中心点上的总压与静压,以及各区域单元的面积,计算出各区域单元面积中心点上的马赫数和流量气动函数,从而求得各区域单元面上的流量,各区域单元上的流量通过如下公式确定:
Figure BDA0003112129510000051
第六步,求椭圆出口截面的流量:
将40个区域单元面上的流量求和,从而得到椭圆出口截面的流量,椭圆出口截面的流量通过如下公式确定:
Figure BDA0003112129510000052
以上的仅是本发明的实施例,方案中公知的具体结构和/或特性等常识在此未作过多描述。应当指出,对于本领域的技术人员来说,在不脱离本发明结构的前提下,还可以作出若干变形和改进,这些也应该视为本发明的保护范围,这些都不会影响本发明实施的效果和专利的实用性。本申请要求的保护范围应当以其权利要求的内容为准,说明书中的具体实施方式等记载可以用于解释权利要求的内容。

Claims (3)

1.一种“3+4+5“米字耙流量计测点布置方法,其特征在于:包括如下步骤:
S1、进气道椭圆截面出口布置有米字耙流量计,米字耙流量计上布置有与椭圆截面短半轴重合的2根垂直测量排架和与椭圆截面长半轴重合的2根水平测量排架,米字耙流量计上还布置有4根分别位于椭圆截面四个象限的倾斜测量排架;
S2、步骤S1的每根垂直测量排架的3个四等分点上等间距布置3个总压测点,每根倾斜测量排架的4个五等分点上等间距布置4个总压测点,每根水平测量排架的5个六等分点上等间距布置5个总压测点;
S3、步骤S1的每根垂直测量排架、倾斜测量排架和水平测量排架的相交处布置有1个总压测点,其中相交处即进气道椭圆截面出口的中心;
S4、步骤S1的每根垂直测量排架、倾斜测量排架和水平测量排架的中心线与进气道椭圆截面出口的交点处布置有1个静压测点。
2.根据权利要求1所述的一种“3+4+5“米字耙流量计测点布置方法,其特征在于:所述倾斜测量排架的中心线与进气道椭圆截面长半轴的夹角通过如下公式确定:
Figure FDA0003863696180000011
其中,a为进气道出口截面椭圆长半轴的长度,b为进气道出口截面椭圆短半轴的长度。
3.根据权利要求1~2中任意一项所述的一种“3+4+5“米字耙流量计测点布置方法,其特征在于:所述测点布置方法的数据处理方法包括如下步骤:
第一步、获得测点数据:
从米字耙流量计测量排架上获得32个测点的总压值,从米字耙流量计的中心测点获得1个总压值,从进气道的椭圆截面出口的静压测点获得8个壁面静压值;
第二步、对进气道的椭圆截面出口进行面积分割:
将椭圆截面出口的面积从周向和径向分割成40个面积相等、且互无间隔互不重叠的区域单元;
第三步、求各区域单元面积中心点的总压:
沿每根测量排架,插值得到该排架上等面积单元区域的面心处总压;
第四步、求各区域单元面积中心点的静压:
对8个壁面静压值进行算数平均得到米字耙流量计中心点的平均静压,由某根测量排架对应的壁面静压与中心点平均静压通过线性插值得到位于该测量排架上的各区域单元面积中心点的静压,从而求得40个区域单元面积中心点的静压;
第五步、求各区域单元上的流量:
根据40个区域单元面积中心点上的总压与静压,以及各区域单元的面积,计算出各区域单元面积中心点上的马赫数和流量气动函数,从而求得各区域单元面上的流量,各区域单元上的流量通过如下公式确定:
Figure FDA0003863696180000021
其中,Pi*为第i个区域单元面积中心点的总压,T*为流场的总温,K为流量公式中与气体比热比相关的常系数,
Figure FDA0003863696180000022
为第i个区域单元的流量,
Figure FDA0003863696180000023
为椭圆出口截面的流量,qi(Ma)为第i个区域单元面积中心对应的空气动力学流量函数;
第六步,求椭圆出口截面的流量:
将40个区域单元面上的流量求和,从而得到椭圆出口截面的流量,椭圆出口截面的流量通过如下公式确定:
Figure FDA0003863696180000024
CN202110654623.0A 2021-06-11 2021-06-11 一种“3+4+5”米字耙流量计测点布置方法 Active CN113405763B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110654623.0A CN113405763B (zh) 2021-06-11 2021-06-11 一种“3+4+5”米字耙流量计测点布置方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110654623.0A CN113405763B (zh) 2021-06-11 2021-06-11 一种“3+4+5”米字耙流量计测点布置方法

Publications (2)

Publication Number Publication Date
CN113405763A CN113405763A (zh) 2021-09-17
CN113405763B true CN113405763B (zh) 2022-11-29

Family

ID=77683617

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110654623.0A Active CN113405763B (zh) 2021-06-11 2021-06-11 一种“3+4+5”米字耙流量计测点布置方法

Country Status (1)

Country Link
CN (1) CN113405763B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114136643B (zh) * 2021-10-20 2024-01-09 中国航发四川燃气涡轮研究院 一种航空发动机空气流量测点布局方法
CN115979356B (zh) * 2022-12-21 2024-08-02 南京航空航天大学 一种适用于超声速/高超声速流道的流量测试系统及测试方法
CN118376381B (zh) * 2024-06-24 2024-08-30 中国空气动力研究与发展中心低速空气动力研究所 一种测压耙、进气道压力试验装置及试验方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101813554A (zh) * 2010-03-29 2010-08-25 南京航空航天大学 可在同一模型上进行测量的进气道实验装置及工作方法
RU2542791C1 (ru) * 2013-08-29 2015-02-27 Федеральное государственное унитарное предприятие "Центральный аэрогидродинамический институт имени профессора Н.Е. Жуковского" (ФГУП "ЦАГИ") Приемник воздушного давления
CN104848904A (zh) * 2015-06-05 2015-08-19 中国航天空气动力技术研究院 进气道流量测量系统
CA2941899A1 (fr) * 2015-09-15 2017-03-15 Jean-Luc Vialatte Mesure des ecoulements d'air le long d'une paroi
EP3168586A1 (fr) * 2015-11-12 2017-05-17 Safran Aircraft Engines Dispositif de mesure de grandeurs aérodynamiques destiné à être placé dans une veine d'écoulement d'une turbomachine
CN107228690A (zh) * 2017-04-10 2017-10-03 中国计量科学研究院 一种不受非测量方向流速影响的皮托管流速计或流量计
WO2017192976A1 (en) * 2016-05-06 2017-11-09 Frohnapfel Dustin J Improved generalized flow profile production
CN107560822A (zh) * 2017-10-31 2018-01-09 中国商用飞机有限责任公司 发动机进气道前缘防冰冰风洞试验设备及其试验方法
CN111157067A (zh) * 2019-12-31 2020-05-15 姜春林 一种流体流量测量方法
CN111551339A (zh) * 2020-04-28 2020-08-18 中国空气动力研究与发展中心高速空气动力研究所 一种翼身融合飞机通气模型畸形流道专用测压耙设计方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BE1027043B1 (fr) * 2019-02-08 2020-09-08 Safran Aero Boosters Sa Mesure de pression totale et temperature totale dans une turbomachine
CN111044252B (zh) * 2019-12-19 2021-12-07 中国航天空气动力技术研究院 一种高精准度进气道流量测量方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101813554A (zh) * 2010-03-29 2010-08-25 南京航空航天大学 可在同一模型上进行测量的进气道实验装置及工作方法
RU2542791C1 (ru) * 2013-08-29 2015-02-27 Федеральное государственное унитарное предприятие "Центральный аэрогидродинамический институт имени профессора Н.Е. Жуковского" (ФГУП "ЦАГИ") Приемник воздушного давления
CN104848904A (zh) * 2015-06-05 2015-08-19 中国航天空气动力技术研究院 进气道流量测量系统
CA2941899A1 (fr) * 2015-09-15 2017-03-15 Jean-Luc Vialatte Mesure des ecoulements d'air le long d'une paroi
EP3168586A1 (fr) * 2015-11-12 2017-05-17 Safran Aircraft Engines Dispositif de mesure de grandeurs aérodynamiques destiné à être placé dans une veine d'écoulement d'une turbomachine
WO2017192976A1 (en) * 2016-05-06 2017-11-09 Frohnapfel Dustin J Improved generalized flow profile production
CN107228690A (zh) * 2017-04-10 2017-10-03 中国计量科学研究院 一种不受非测量方向流速影响的皮托管流速计或流量计
CN107560822A (zh) * 2017-10-31 2018-01-09 中国商用飞机有限责任公司 发动机进气道前缘防冰冰风洞试验设备及其试验方法
CN111157067A (zh) * 2019-12-31 2020-05-15 姜春林 一种流体流量测量方法
CN111551339A (zh) * 2020-04-28 2020-08-18 中国空气动力研究与发展中心高速空气动力研究所 一种翼身融合飞机通气模型畸形流道专用测压耙设计方法

Also Published As

Publication number Publication date
CN113405763A (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
CN113405763B (zh) 一种“3+4+5”米字耙流量计测点布置方法
CN113125103B (zh) 一种41测点等距分布的椭圆截面流量计数据处理方法
CN111551339B (zh) 一种翼身融合飞机通气模型畸形流道专用测压耙设计方法
US20140064330A1 (en) Apparatus and method for measuring total air temperature within an airflow
Senoo et al. Behavior of severely asymmetric flow in a vaneless diffuser
CN107917793A (zh) 一种实验用多排测压耙装置
CN111780949B (zh) 基于cfd分析的高速进气道前体风洞实验总压修正方法
CN104596769A (zh) 航空发动机在室内试车台试车的真实推力测算方法
CN106545411B (zh) 模拟畸变流场的直连试验装置的设计方法
CN111498141B (zh) 一种基于微型探针实现气流角度实时监测的方法与装置
Poensgen et al. Three-Dimensional Wake Decay Inside of a Compressor Cascade and Its Influence on the Downstream Unsteady Flow Field: Part I—Wake Decay Characteristics in the Flow Passage
CN114528778B (zh) 基于钟形口流量管速度分布模型的空气流量测量方法
CN219736607U (zh) 一种加装外层滞止罩的铂电阻总温探针头部
CN111089669A (zh) 一种用于精准测量涡轮级效率的装置及方法
CN106017582A (zh) 一种基于管径二分原理的皮托管流量测量方法
CN113125800B (zh) 基于皮托管的风速风向测量方法
CN105760647A (zh) 一种质量加权平均值计算方法
Schulze et al. NACA 65-series compressor rotor performance with varying annulus-area ratio, solidity, blade angle, and Reynolds number and comparison with cascade results
CN212458736U (zh) 一种用于精准测量涡轮级效率的装置
Gregory-Smith An investigation of annulus wall boundary layers in axial flow turbomachines
Poensgen et al. Three-Dimensional Wake Decay Inside of a Compressor Cascade and its Influence on the Downstream Unsteady Flow Field: Part I—Wake Decay Characteristics in the Flow Passage
KENDALL Experimental study of laminar boundary layer receptivity to a traveling pressure field
Copp Effects of inlet wall contour on the pressure recovery of a 10 degrees 10-inch-inlet-diameter conical diffuser
Fessler et al. Preliminary survey of compressor rotor-blade wakes and other flow phenomena with a hot-wire anemometer
Allen Experimental investigation of loss in an annular cascade of turbine-nozzle blades of free vortex design

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