CN113125102B - 内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法 - Google Patents

内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法 Download PDF

Info

Publication number
CN113125102B
CN113125102B CN202110313683.6A CN202110313683A CN113125102B CN 113125102 B CN113125102 B CN 113125102B CN 202110313683 A CN202110313683 A CN 202110313683A CN 113125102 B CN113125102 B CN 113125102B
Authority
CN
China
Prior art keywords
section
ellipse
area
measured
cross
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
CN202110313683.6A
Other languages
English (en)
Other versions
CN113125102A (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 CN202110313683.6A priority Critical patent/CN113125102B/zh
Publication of CN113125102A publication Critical patent/CN113125102A/zh
Application granted granted Critical
Publication of CN113125102B publication Critical patent/CN113125102B/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
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

本发明公开了一种椭圆截面米字耙流量计5×8测点等面积分布方法,具体方案为:包括以下步骤:S1:用4个中心、长轴所在直线分别与待测截面椭圆中心、长轴所在直线重合的相似椭圆将待测截面椭圆沿径向分割成面积相等的1个椭圆和4个椭圆环;S2:以待测截面椭圆中心为端点,向待测截面椭圆圆周线绘制8条第一射线,8条第一射线将待测截面椭圆沿周向分割成8个面积相等的椭圆扇形;S3:结合S1和S2得到面积相等的40个彼此互不重叠、互无间隙的区域单元;S4:求出每个区域单元的面积中心用于布置总压测点。该方法既解决了椭圆截面流量计如何合理分布测点的全新问题,又可以在后处理过程中很方便地按照等面积分布原则求得椭圆截面的流量。

Description

内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法
技术领域
本发明属于实验空气动力学技术领域,更具体地说,它涉及一种内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法及测量方法。
背景技术
流量测量是进气道设计校核、性能评估及数值计算方法检验的重要途径。内转进气道是目前超高速飞行器机体/推进一体化应用越来越广泛的一种进气道,这种进气道的出口通常使用椭圆截面或圆截面。与轴对称进气道或内转进气道出口的圆形截面、二维或侧压进气道出口的矩形截面不同,内转进气道出口的椭圆截面在形状上具有特殊性。在对这类进气道开展流量测量时,如何在其椭圆截面出口布置流量总压测点是一个全新的问题。
目前常见的进气道流量计主要有两种:一种是气动喉道法(即声速截面法);另外一种为总压耙法(也称排架法),即在测量截面上布置多个总压耙,每个耙上布置多个总压探针,每个总压探针对应一个确定的单元面积,且各个总压探针分别位于对应单元的面积中心,同时在测量截面的壁面布置多个静压测点,根据静压总压之间的气动关系求得各个总压测点上的马赫数,与测得的总压、单元面积一起计算单元流量,再求和得到截面流量。前者多用于亚声速流道,后者除了可用于亚声速流道外,更多用于超声速流道。本文针对的是椭圆截面出口的内转进气道,其出口是超声速流,因此流量测量适宜使用总压耙法。
当前,总压耙流量计按照测量截面的形状可分为两种:一种是矩形截面流量耙,在测量截面上布置多个平行的排架,每个排架上布置多个等距分布的总压测点,各个总压测点正好对应着等面积分布的各个区域单元的面积中心;另一种是圆形截面流量耙,在测量截面上沿圆周均匀布置多根排架,每根排架上按照等面积法布置多个测点,例如常见的“米”字耙流量计。
现有总压耙流量计的测点布置方案,不适用于椭圆截面出口流量测量试验,主要表现在以下几个方面:
矩形截面流量耙:如果在椭圆截面出口平面上采用矩形截面流量耙的平行等距测点布置方式,一方面会影响测量结果的准确性,在椭圆截面上除长轴和短轴外,其余区域的测点会出现分布不均匀的情况;另一方面会增加测量数据后处理的难度,等距均匀布置的测点所代表单元面积在边界处难以合理给出,后期插值也存在科学性问题。
圆形截面流量耙:如果在椭圆截面出口平面上采用圆形截面流量耙的测点布置方式,由于椭圆形状的特殊性,同样会降低流量测量数据的可靠性,一方面,除位于椭圆中心附近的几个测点外,其余测点会暴露在椭圆截面出口之外;另一方面,即使按照在椭圆截面内的某一个圆周上均匀设置总压测点,也无法做到相邻排架之间的区域面积相等。
因此,当内转进气道的出口形状为椭圆时,如何合理布置椭圆截面流量计的总压测点,使得每个测点依然能够合理代表一个区域单元,且所有测点代表的区域单元既不重叠又无遗漏地覆盖整个椭圆截面是个全新的技术问题。
发明内容
为解决上述技术问题,本发明提供一种内转进气道椭圆截面米字耙流量计5×8测点的分布方法,充分使用等面积原则分布测点,既解决了椭圆截面流量计如何合理分布测点的全新问题,又可以在后处理过程中能够按照等面积分布原则直接求得椭圆截面的流量。这种求流量的方法还为按流量加权求椭圆截面马赫数与总压恢复系数提供了技术途径。
本发明的上述技术目的是通过以下技术方案得以实现的:
一种内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法,包括以下步骤:
S1:用4个中心、长轴所在直线分别与待测截面椭圆中心、长轴所在直线重合的相似椭圆将待测截面椭圆沿径向分割成面积相等的1个椭圆和4个椭圆环;
S2:以待测截面椭圆中心为端点,向待测截面椭圆圆周线绘制8条第一射线,8条第一射线将待测截面椭圆沿周向分割成8个面积相等的椭圆扇形;
S3:结合S1和S2得到面积相等的40个彼此互不重叠、互无间隙的区域单元;
S4:求出每个区域单元的面积中心用于布置总压测点。
进一步地,S2具体包括以下步骤:
S201:以待测截面椭圆中心为坐标原点O,沿待测截面椭圆长轴方向从左向右为X轴的正方向,沿待测截面椭圆短轴方向从下向上为Y轴的正方向建立平面直角坐标系;
S202:以X轴正方向为始边,逆时针方向旋转到第j条第一射线O l j 的夹角αj为第j条第一射线的方位角,通过公式确定第j条第一射线的方位,即:
Figure 899034DEST_PATH_IMAGE001
]
S203:以坐标原点O为端点,以αj为方位角向待测截面椭圆圆周线绘制8条第一射线Ol 1~Ol 8,8条第一射线将待测截面椭圆沿周向分割成8个面积相等的椭圆扇形。
进一步地,S4具体包括以下步骤:
S401:用9个中心、长轴所在直线分别与待测截面椭圆中心、长轴所在直线重合的相似椭圆将待测截面椭圆沿径向分割成面积相等的1个椭圆和9个椭圆环,从待测截面椭圆中心向外依次形成第1椭圆~第10椭圆;
S402:以X轴正方向为始边,逆时针方向旋转到第j条第二射线OLj的夹角θj为第j条第二射线的方位角,通过公式确定第j条第二射线的方位,即:
Figure 421151DEST_PATH_IMAGE002
S403:以坐标原点O为端点,以θj为方位角向待测截面椭圆圆周线绘制8条第二射线OL1~OL8
S404:S403的第二射线OLj与S401的奇数序号椭圆轮廓线相交的40个交点,即40个区域单元的面积中心布置椭圆截面米字耙流量计的总压测点。
一种内转进气道椭圆截面流量的测量方法,基于上述的椭圆截面米字耙流量计5×8测点等面积分布方法,其特征在于,设置8根测量耙,每根测量耙上设置5个总压测点,得到40个总压测点,每个总压测点对应设置在S404得到的40个区域单元的面积中心上进行测量。
综上所述,本发明具有以下有益效果:
本发明技术直接为椭圆截面提供了一种米字耙流量计5×8测点分布方法,所采用的测点等面积分布原则使得椭圆截面上米字耙以及耙上测点分布科学合理,填补了椭圆截面流量计总压测点分布方法库的空白,也为进一步的椭圆截面平均马赫数、平均总压恢复系数后处理提供了技术途径。
附图说明
图1为本发明实施例中的内转进气道椭圆出口截面的径向分割图;
图2为本发明实施例中的内转进气道椭圆出口截面的周向分割图;
图3为本发明实施例中的内转进气道椭圆出口截面米字耙流量计5×8测点定位图。
具体实施方式
本说明书及权利要求并不以名称的差异来作为区分组件的方式,而是以组件在功能上的差异来作为区分的准则。如在通篇说明书及权利要求当中所提及的“包括”为一开放式用语,故应解释成“包括但不限定于”。“大致”是指在可接受的误差范围内,本领域技术人员能够在一定误差范围内解决所述技术问题,基本达到所述技术效果。
本说明书及权利要求的上下左右等方向名词,是结合附图以便于进一步说明,使得本申请更加方便理解,并不对本申请做出限定,在不同的场景中,上下、左右、里外均是相对而言。
下面结合附图对本发明作进一步的详细说明:
内转进气道椭圆截面出口米字耙流量计5×8测点等面积分布方法概述
在内转进气道的椭圆出口截面上采用辐射分布的8根耙,每根耙分布5个总压测点,共布置40个总压测点。40个总压测点分别位于40个面积相等的不同区域单元的面积中心,这40个区域单元彼此互不重叠、互无间隙的填满整个椭圆截面。则有:
Figure 87756DEST_PATH_IMAGE003
(1)
Figure 120565DEST_PATH_IMAGE004
(2)
每个测点的测量值代表相应区域单元的测量值,40个单元的流量求和就得到了整个椭圆截面的流量。
椭圆截面出口米字耙流量计5×8测点等面积分布方法
(1)椭圆出口截面的等面积分割
根据米字耙流量计每根耙的测点数对椭圆出口截面进行等面积径向分割,在径向分割的基础上,根据米字流量耙的根数对椭圆出口截面进行等面积周向分割。
①椭圆出口截面的等面积径向分割
椭圆出口截面米字耙流量计上每根耙布置5个总压测点,将椭圆出口截面面积平均分为5等分。
方法:用4个中心、长轴所在直线分别与待测截面椭圆中心、长轴所在直线重合的相似椭圆将待测截面椭圆平均分成5个面积相等的部分,由椭圆中心O向外依次得到椭圆2、椭圆环2、椭圆环3、椭圆环4、椭圆环5,见图1。
相似椭圆的长半轴及短半轴确定:
椭圆2、椭圆环2、椭圆环3、椭圆环4、椭圆环5的面积均为原椭圆面积的1/5,则椭圆2、椭圆4、椭圆6、椭圆8、原椭圆的面积比为1:2:3:4:5,即:
Figure 290646DEST_PATH_IMAGE005
(3)
Figure 921347DEST_PATH_IMAGE006
(4)
由此可以得到相似椭圆的长半轴及短半轴长度:
Figure 103279DEST_PATH_IMAGE007
(5)
根据内转进气道出口截面椭圆的长半轴a和短半轴b,结合公式(5)可以计算出相似椭圆的长半轴ai和短半轴bi
②椭圆出口截面的等面积周向分割
如图2所示,根据米字流量耙的根数,将椭圆出口截面沿周向平均分割为8等分。
方法:建立XOY坐标系,椭圆的长轴为X轴,短轴为Y轴,椭圆中心O为坐标原点,以坐标原点为端点向出口截面椭圆作8条射线O l j ,根据椭圆的对称性规定这些射线既关于X轴对称,又关于Y轴对称。这8条射线O l j 将椭圆出口截面沿周向平均分割为8个面积相等的椭圆扇形。
每根射线与椭圆2、椭圆4、椭圆6、椭圆8、原椭圆相交后,就把椭圆出口截面分割为40个互不重叠、互无间隙的等面积区域单元,见图2。
射线O l j 的方位角α j 通过公式给出:
Figure 657888DEST_PATH_IMAGE001
] (6)
(2)椭圆出口截面米字耙流量计总压测点的定位
椭圆出口截面米字耙流量计上40个总压测点,分别位于40个面积相等的不同区域单元的面积中心,这40个区域单元彼此互不重叠、互无间隙的填满整个椭圆截面。
①米字耙流量计总压测点的径向定位
椭圆出口截面米字耙流量计每根耙上的5个测点,分别位于椭圆出口截面径向分割后的5个等面积椭圆或椭圆环的面积中心区域。
测点的径向定位方法:用9个中心、长轴所在直线分别与待测截面椭圆中心、长轴所在直线重合的相似椭圆将待测截面椭圆沿径向分割成面积相等的1个椭圆和9个椭圆环,从中心O向外依次形成椭圆1、椭圆2、…、椭圆9、椭圆10,其中奇数序号的椭圆轮廓线就是米字耙流量计总压测点的径向定位曲线。
用于总压测点定位的相似椭圆长半轴和短半轴长度为:
Figure 869426DEST_PATH_IMAGE008
(7)
根据内转进气道出口截面椭圆的长半轴a和短半轴b,结合公式(7)可以计算出米字耙流量计总压测点径向定位椭圆的长半轴ai和短半轴bi
②米字耙流量计总压测点的周向定位
椭圆出口截面米字耙流量计包括8根耙,将其中的横向与纵向四根耙分别设置在椭圆截面的长轴与短轴上,它们正好分别位于4个面积相等的椭圆扇形的对称轴上;其余四根耙分别位于剩余四个等面积椭圆扇形的面积二等分线上。
总压测点的周向定位方法:建立XOY坐标系,椭圆的长轴为X轴,短轴为Y轴,椭圆中心O为坐标原点,以坐标原点为端点向出口截面椭圆做射线OL j ,射线OL j 的方位角θ j 通过公式给出:
Figure 77554DEST_PATH_IMAGE009
(8)
射线OLj与总压测点径向定位相似椭圆轮廓线形成的交点,即内转进气道椭圆出口截面40个面积相等的不同区域单元的面积中心,就是米字耙流量计上40个总压测点的中心位置,见图3。
实施例公式中相关参数具体说明如表1所示:
Figure 203904DEST_PATH_IMAGE010
表1
本具体实施例仅仅是对本发明的解释,其并不是对本发明的限制,本领域技术人员在阅读完本说明书后可以根据需要对本实施例做出没有创造性贡献的修改,但只要在本发明的权利要求范围内都受到专利法的保护。

Claims (2)

1.内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法,其特征在于:包括以下步骤:
S1:用4个中心、长轴所在直线分别与待测截面椭圆中心、长轴所在直线重合的相似椭圆将待测截面椭圆沿径向分割成面积相等的1个椭圆和4个椭圆环;
S2:以待测截面椭圆中心为端点,向待测截面椭圆圆周线绘制8条第一射线,8条第一射线将待测截面椭圆沿周向分割成8个面积相等的椭圆扇形;
S3:结合S1和S2得到面积相等的40个彼此互不重叠、互无间隙的区域单元;
S4:求出每个区域单元的面积中心用于布置总压测点;
S2具体包括以下步骤:
S201:以待测截面椭圆中心为坐标原点O,沿待测截面椭圆长轴方向从左向右为X轴的正方向,沿待测截面椭圆短轴方向从下向上为Y轴的正方向建立平面直角坐标系;
S202:以X轴正方向为始边,逆时针方向旋转到第j条第一射线O l j 的夹角αj为第j条第一射线的方位角,通过公式确定第j条第一射线的方位,即:
Figure DEST_PATH_IMAGE001
]
S203:以坐标原点O为端点,以αj为方位角向待测截面椭圆圆周线绘制8条第一射线Ol 1~Ol 8,8条第一射线将待测截面椭圆沿周向分割成8个面积相等的椭圆扇形;
S4具体包括以下步骤:
S401:用9个中心、长轴所在直线分别与待测截面椭圆中心、长轴所在直线重合的相似椭圆将待测截面椭圆沿径向分割成面积相等的1个椭圆和9个椭圆环,从待测截面椭圆中心向外依次形成第1椭圆~第10椭圆;
S402:以X轴正方向为始边,逆时针方向旋转到第j条第二射线OLj的夹角θj为第j条第二射线的方位角,通过公式确定第j条第二射线的方位,即:
Figure DEST_PATH_IMAGE003
S403:以坐标原点O为端点,以θj为方位角向待测截面椭圆圆周线绘制8条第二射线OL1~OL8
S404:S403的第二射线OLj与S401的奇数序号椭圆轮廓线相交的40个交点,即40个区域单元的面积中心布置椭圆截面米字耙流量计的总压测点;
其中,参数 a 为内转进气道出口截面椭圆长半轴 的长度,参数 b 为内转进气道出口截面椭圆短半轴的长度。
2.一种内转进气道椭圆截面流量的测量方法,基于权利要求1所述的内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法,其特征在于,设置8根测量耙,每根测量耙上设置5个总压测点,得到5×8个总压测点,每个总压测点对应设置在S404得到的40个区域单元的面积中心上进行测量。
CN202110313683.6A 2021-03-24 2021-03-24 内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法 Active CN113125102B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110313683.6A CN113125102B (zh) 2021-03-24 2021-03-24 内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110313683.6A CN113125102B (zh) 2021-03-24 2021-03-24 内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法

Publications (2)

Publication Number Publication Date
CN113125102A CN113125102A (zh) 2021-07-16
CN113125102B true CN113125102B (zh) 2022-12-02

Family

ID=76774076

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110313683.6A Active CN113125102B (zh) 2021-03-24 2021-03-24 内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法

Country Status (1)

Country Link
CN (1) CN113125102B (zh)

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10160619A (ja) * 1996-11-29 1998-06-19 Ishikawajima Harima Heavy Ind Co Ltd 3次元流れ計測装置とこれを用いた流れ分布計測装置
CN1257390C (zh) * 2003-09-04 2006-05-24 温汉璋 一种具有提速功能的均速流量测量装置
CN101701837B (zh) * 2009-11-20 2011-04-06 中国航空工业空气动力研究院 一种高精度大范围气流流量计
CN202101706U (zh) * 2011-04-22 2012-01-04 燕山大学 一种内管道为椭圆型的电磁流量计传感器
CN104089750B (zh) * 2014-07-15 2016-04-27 首钢总公司 一种顶燃式热风炉流场模拟装置及其模拟方法
CN104848904B (zh) * 2015-06-05 2018-08-07 中国航天空气动力技术研究院 进气道流量测量系统
JP6837026B2 (ja) * 2018-03-19 2021-03-03 株式会社豊田中央研究所 風圧検出位置演算装置、移動体運動制御装置、風圧検出位置演算方法、風圧検出位置演算プログラム、及び移動体運動制御プログラム
IT201700122220A1 (it) * 2017-10-27 2019-04-27 Vertical Dwc Llc Galleria del vento ricircolante trasferibile.
CN107917793B (zh) * 2017-11-14 2020-09-25 中国空气动力研究与发展中心高速空气动力研究所 一种实验用多排测压耙装置
CN111044252B (zh) * 2019-12-19 2021-12-07 中国航天空气动力技术研究院 一种高精准度进气道流量测量方法

Also Published As

Publication number Publication date
CN113125102A (zh) 2021-07-16

Similar Documents

Publication Publication Date Title
Connolly Analytical molecular surface calculation
CN105333855B (zh) 一种测量、计算涡轮导向叶片和导向器喉道面积的方法
CN107228610B (zh) 一种针对筒状曲面的螺旋式扫描测量轨迹规划方法
CN106023156A (zh) 基于检测特征的点云模型与cad模型的配准方法
CN111551339B (zh) 一种翼身融合飞机通气模型畸形流道专用测压耙设计方法
CN102997875B (zh) 一种复杂螺旋曲面高精度测量不确定度的分析方法
CN107917793A (zh) 一种实验用多排测压耙装置
CN105627923A (zh) 一种基于激光测距的叶片稠密点云获取的扫描路径规划方法
CN113125102B (zh) 内转进气道椭圆截面米字耙流量计5×8测点等面积分布方法
Miura et al. Grid drawings of 4-connected plane graphs
CN113125103A (zh) 一种41测点等距分布的椭圆截面流量计数据处理方法
CN113405763B (zh) 一种“3+4+5”米字耙流量计测点布置方法
CN103292654B (zh) 一种计算圆柱体零件作用尺寸的方法
CN106844896B (zh) 一种适用于旋成体外形的来流参数确定方法
Shortley et al. Numerical solution of axisymmetrical problems, with applications to electrostatics and torsion
CN107091603A (zh) 一种测量整体叶盘的方法
CN206756752U (zh) 一种用于超超临界汽轮机低压隔板超声相控阵检测的对比试块
CN114136643A (zh) 一种航空发动机空气流量测点布局方法
CN111487031B (zh) 基于微型扇形探针实现三维流场气流角度实时监测的装置与方法
Smith et al. The use of analytic surfaces for the design of centrifugal impellers by computer graphics
Kayser et al. Neighborhood Analysis of Low Magnification Structures (Glands) in Healthy, Adenomatous, and Carcinomatous Colon Mucos
Dittrich Discharge Coefficients for Combustor-liner Air-entry Holes II: Flush Rectangular Holes, Step Louvers, and Scoops
McNally et al. FORTRAN program for computing coordinates of circular arc single and tandem turbomachinery blade sections on a plane
CN104634294A (zh) 一种曲线槽槽轮几何误差的检测与评定方法
CN114252035B (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