CN114278460B - 一种轴对称天线罩包罩试验喷管型线设计方法 - Google Patents

一种轴对称天线罩包罩试验喷管型线设计方法 Download PDF

Info

Publication number
CN114278460B
CN114278460B CN202111594629.XA CN202111594629A CN114278460B CN 114278460 B CN114278460 B CN 114278460B CN 202111594629 A CN202111594629 A CN 202111594629A CN 114278460 B CN114278460 B CN 114278460B
Authority
CN
China
Prior art keywords
control point
flow
calculating
ith
gap
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
CN202111594629.XA
Other languages
English (en)
Other versions
CN114278460A (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.)
Beijing Electromechanical Engineering Research Institute
Original Assignee
Beijing Electromechanical Engineering Research Institute
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 Beijing Electromechanical Engineering Research Institute filed Critical Beijing Electromechanical Engineering Research Institute
Priority to CN202111594629.XA priority Critical patent/CN114278460B/zh
Publication of CN114278460A publication Critical patent/CN114278460A/zh
Application granted granted Critical
Publication of CN114278460B publication Critical patent/CN114278460B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

Landscapes

  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
  • Details Of Aerials (AREA)

Abstract

本发明涉及一种轴对称天线罩包罩试验喷管型线设计方法,属于飞行器结构试验技术领域,解决了现有技术中天线罩包罩试验中能量利用率低、成本较高的问题。该方法包括:设置喉道间隙尺寸以及喷管上各控制点的轴向位置,计算喉道面积,并根据各控制点压力计算来流总压及各控制点的流动马赫数;根据所述各控制点的流动马赫数和所述喉道面积,计算各控制点的无粘流道间隙尺寸;根据所述来流总压及各控制点的流动马赫数,计算各控制点的流道间隙边界层厚度;根据所述各控制点的轴向位置、所述无粘流道间隙尺寸和所述流道间隙边界层厚度,计算各控制点坐标,将得到的所述控制点拟合成连续光滑曲线,生成喷管型线。

Description

一种轴对称天线罩包罩试验喷管型线设计方法
技术领域
本发明涉及飞行器结构试验技术领域,尤其涉及一种轴对称天线罩包罩试验喷管型线设计方法。
背景技术
随着现代飞行器飞行速度越来越快,位于飞行器头部的天线罩面临着严酷的力热载荷环境。为确保天线罩结构安全,保证天线罩内部设备正常工作,需要对天线罩结构进行地面考核试验。传统的天线罩力热联合试验中,力载荷加载方式主要为皮囊或者卡箍,热载荷加载方式主要为石英灯,无法模拟天线罩表面压力分布,也没有气流的剪切、冲刷效应;射流试验能够真实模拟天线罩表面压力,但能量利用率低、成本较高。
因此,现有技术中缺少一种能够节省气流流量和加热功率,提升能量利用率的轴对称天线罩包罩试验喷管型线设计方法。
发明内容
鉴于上述的分析,本发明实施例旨在提供一种轴对称天线罩包罩试验喷管型线设计方法,用以解决现有天线罩包罩试验中能量利用率低、成本较高的问题。
一方面,本发明实施例提供了一种轴对称天线罩包罩试验喷管型线设计方法,包括:
设置喉道间隙尺寸以及喷管上各控制点的轴向位置,计算喉道面积,并根据各控制点压力计算来流总压及各控制点的流动马赫数;
根据所述各控制点的流动马赫数和所述喉道面积,计算各控制点的无粘流道间隙尺寸;
根据所述来流总压及各控制点的流动马赫数,计算各控制点的流道间隙边界层厚度;
根据所述各控制点的轴向位置、所述无粘流道间隙尺寸和所述流道间隙边界层厚度,计算各控制点坐标,将得到的所述控制点拟合成连续光滑曲线,生成喷管型线。
进一步地,所述根据所述各控制点的轴向位置、所述无粘流道间隙尺寸和所述流道间隙边界层厚度,计算各控制点坐标,包括:
根据所述无粘流道间隙尺寸,以及所述流道间隙边界层厚度计算各控制点的流道间隙尺寸;
根据所述各控制点的轴向位置、所述流道间隙尺寸和所述控制点对应的天线罩外表面半径,计算各控制点坐标。
进一步地,所述各控制点坐标,表达为:
(xi,yi)=(xi,Rii)
δi=δi'+2δi
其中,xi为第i个控制点与天线罩顶点的轴向距离,yi为第i个控制点与天线罩旋转轴间的纵向距离,Ri为第i个控制点对应的天线罩外表面半径,δi为第i个控制点的流道间隙尺寸,δi'为第i个控制点的无粘流道间隙尺寸,δi”为第i个控制点的流道间隙边界层厚度。
进一步地,所述根据所述各控制点的流动马赫数和所述喉道面积,计算各控制点无粘流道间隙尺寸,包括:
根据第i个控制点的马赫数、喉道面积和气体比热比计算第i个控制点处的无粘流道面积;
根据所述第i个控制点处的无粘流道面积和第i个控制点对应的天线罩外表面半径计算所述无粘流道间隙尺寸。
进一步地,第i个控制点处的所述无粘流道间隙尺寸δi',表达为:
Figure BDA0003430173720000031
Figure BDA0003430173720000032
其中,σi为第i个控制点处的无粘流道面积,Ri为第i个控制点对应的天线罩外表面半径,σ*为喉道面积,Mai为第i个控制点的马赫数,γ为气体比热比。
进一步地,所述根据所述来流总压及各控制点的流动马赫数,计算各控制点的流道间隙边界层厚度,包括:
根据第i个控制点的马赫数、来流总温和气体比热比计算第i个控制点处的气流温度;
根据所述第i个控制点处的气流温度、第i个控制点压力和理想气体常数计算第i个控制点处的气流密度;
根据所述第i个控制点处的气流密度和气流温度计算第i个控制点处的气流雷诺数;
根据所述第i个控制点处的气流雷诺数和第i个控制点与天线罩顶点的轴向距离,计算第i个控制点的流道间隙边界层厚度。
进一步地,所述第i个控制点的流道间隙边界层厚度δi”,表达为:
Figure BDA0003430173720000033
Figure BDA0003430173720000041
Figure BDA0003430173720000042
Figure BDA0003430173720000043
其中,xi为第i个控制点与天线罩顶点的轴向距离,Rexi为第i个控制点处的气流雷诺数,Ti为第i个控制点处的气流温度,ρi为第i个控制点处的气流密度,pi为第i个控制点压力,R为理想气体常数,T0为来流总温,γ为气体比热比,Mai为第i个控制点的马赫数。
进一步地,设置喉道间隙尺寸以及各控制点的轴向位置,计算喉道面积,并根据各控制点压力计算来流总压及各控制点的流动马赫数,包括:
根据待测天线罩尺寸和电弧加热器的质量流量设置喉道位置、喉道间隙尺寸初始值以及各控制点的轴向位置;
根据所述喉道间隙尺寸计算喉道面积;
根据所述各控制点压力,确定来流总压及各控制点的流动马赫数,并满足天线罩顶点和喉道之间控制点的马赫数小于1,喉道外远离顶点方向的控制点的马赫数大于1;
根据所述喉道面积及来流总压,计算空气质量流量,判断所述空气质量流量是否小于电弧加热器所能提供的最大空气质量流量Gmax
当所述空气质量流量大于所述最大空气质量流量Gmax时,减小喉道间隙尺寸,重新计算空气质量流量,直至所述空气质量流量小于所述最大空气质量流量Gmax
其中,所述各控制点压力为天线罩上与各控制点轴向位置相同的位置处的考核压力。
进一步地,所述根据所述各控制点压力,确定来流总压及各控制点的流动马赫数包括:根据如下公式确定来流总压及各控制点的流动马赫数:
Figure BDA0003430173720000051
其中,p0为来流总压,γ为气体比热比,pi为第i个控制点压力,Mai为第i个控制点的马赫数。
进一步地,所述空气质量流量G,表达为:
Figure BDA0003430173720000052
Figure BDA0003430173720000053
其中,p0为来流总压,σ*为喉道面积,Cd为声速流量系数,T*为喉道处气流温度,T0为来流总温,γ为气体比热比,pi为第i个控制点压力,Mai为第i个控制点的马赫数。
与现有技术相比,本发明至少可实现如下有益效果之一:
1、本发明利用一维变截面可压缩流动和边界层理论,设计包罩喷管型线,在电弧加热器产生的气流条件下可准确模拟天线罩表面压力分布,对天线罩结构承载能力验证更为真实;
2、本发明通过天线罩外包覆一个随型包罩喷管,让电弧加热器产生的热气流在喷管和天线罩之间的缝隙内流动,可以在较小的加热功率和空气流量条件下,实现较高的压力和温度参数,能够显著节省气流流量和加热功率,能量利用率提升几十倍以上。
本发明中,上述各技术方案之间还可以相互组合,以实现更多的优选组合方案。本发明的其他特征和优点将在随后的说明书中阐述,并且,部分优点可从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过说明书以及附图中所特别指出的内容中来实现和获得。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制,在整个附图中,相同的参考符号表示相同的部件。
图1为本申请一个实施例所示的轴对称天线罩包罩试验喷管型线设计方法流程图;
图2为本申请一个实施例所示的天线罩包罩喷管型线示意图;
图3为本申请一个实施例所示的天线罩与喷管构成的流道间隙示意图;
图4为本申请一个实施例所示的流道间隙中无粘流道及边界层示意图;
图5为本申请一个实施例所示的某天线罩包罩喷管型线及压力仿真云图;
图6为本申请一个实施例所示的喷管实现压力与考核要求对比图。
具体实施方式
下面结合附图来具体描述本发明的优选实施例,其中,附图构成本申请一部分,并与本发明的实施例一起用于阐释本发明的原理,并非用于限定本发明的范围。
本发明的一个具体实施例,公开了一种轴对称天线罩包罩试验喷管型线设计方法,如图1所示,包括:
S10、设置喉道间隙尺寸以及喷管上各控制点的轴向位置,计算喉道面积,并根据各控制点压力计算来流总压及各控制点的流动马赫数;
具体地,天线罩的表面压力一般从天线罩顶点开始,往后压力逐渐降低。如距天线罩顶点x=300mm截面压力为0.2MPa,距天线罩顶点x=400mm截面压力为0.15MPa,如图2所示,包罩喷管是包覆在天线罩外面的,二者之间形成环形的缝隙流道,电弧加热器产生的热气流在喷管和天线罩之间的缝隙内流动,可以在较小的加热功率和空气流量条件下,实现较高的压力和温度参数,既能模拟天线罩表面压力分布,也能显著降低试验成本。本实施例就是以已知截面的压力为输入,来进行轴对称天线罩包罩试验喷管型线设计。
具体地,设置喉道间隙尺寸以及各控制点的轴向位置,计算喉道面积,并根据各控制点压力计算来流总压及各控制点的流动马赫数,包括:
S101、根据待测天线罩尺寸及电弧加热器质量流量设置喉道位置、喉道间隙尺寸初始值以及各控制点的轴向位置;
具体地,对于亚声速包罩试验,喉道一般位于天线罩根部的金属框上,喉道间隙尺寸一般为5~15mm,喉道间隙尺寸主要取决于天线罩的尺寸以及电弧加热器所能提供的质量流量,电弧加热器所能提供的流量是一定的,质量流量与天线罩半径、喉道间隙尺寸成正比,因此如果天线罩半径较大,流道间隙尺寸须小一些。流道面积应在试验有效考核区域内最小,以保证气流在喉道位置达到声速,可选地,喉道尺寸可先根据电弧加热器的典型流量(即电弧加热器所能提供的最大流量的一半)进行估算。
轴对称天线罩母线上的压力变化较为平缓,设置5~10个点即能满足设计要求,如图3所示,喉道附近可适当加密,在设置控制点时应注意控制点的设置应能够体现天线罩几何外形和表面压力的变化,同时方便拟合成光滑的喷管型线。可选地,本实施例中选用8个控制点。
S102、根据所述喉道间隙尺寸计算喉道面积;
具体地,喉道在天线罩根部金属环上,喉道面积与喉道位置对应的天线罩外表面半径和喉道间隙尺寸有关。喉道面积为喉道位置对应的天线罩外表面形成的圆沿喉道向圆外扩展喉道间隙长度形成的圆环的面积;所述喉道位置对应的天线罩外表面形成的圆是指天线罩上轴向位置与喉道轴向位置相同的位置处形成的圆,当天线罩外形和喉道间隙尺寸确定后,喉道面积就是一定的。
S103、根据所述各控制点压力,确定来流总压及各控制点的流动马赫数,并满足天线罩顶点和喉道之间控制点的马赫数小于1,喉道外远离顶点方向的控制点的马赫数大于1;
具体地,对于亚声速包罩试验,流道面积应在试验有效考核区域(即天线罩非金属部分)内气流为亚声速(即马赫数小于1),在喉道外远离顶点方向的控制点处气流为超声速(即马赫数大于1),具体地,使用一维变截面可压缩流动关系式,调整来流总压以及天线罩母线控制点处的流动马赫数,使天线罩母线控制点处的压力满足上述试验考核要求;设计过程中可根据需求设置多组数据进行对比,一方面模拟尽可能多的点压力,另一方面要保证有效考核区域内压力最低的点马赫数仍小于1。更具体地,一维变截面可压缩流动关系式中的总压~静压关系式,表达为:
Figure BDA0003430173720000091
其中,p0为来流总压,单位为帕(Pa);γ为气体比热比;pi为第i个控制点压力,单位为帕(Pa);Mai为第i个控制点的马赫数。
可选地,来流总压p0依经验设置初始值为260kPa;气体比热比γ在来流总温不高于1000℃时取1.4。
需要说明的是,上述控制点压力是指天线罩上与各控制点轴向位置相同的位置处的考核压力,是根据试验需要预先设置好的。
S104、根据所述喉道面积及来流总压,计算空气质量流量,判断所述空气质量流量是否小于电弧加热器所能提供的最大空气质量流量Gmax;最大空气质量流量Gmax是电弧加热器的重要性能参数,其受管道直径、气源压力、结构强度限制。国内大多数电弧加热器最大流量不超过10kg/s。
所述根据所述喉道面积及来流总压,计算空气质量流量,包括:根据来流总温和气体比热比计算喉道处气流温度;根据所述喉道处气流温度、来流总温、喉道面积和声速流量系数计算空气流量质量。所述空气质量流量G,表达为:
Figure BDA0003430173720000092
Figure BDA0003430173720000093
其中,p0为来流总压;σ*为喉道面积,单位为平方米(m2);Cd为声速流量系数;T*为喉道处气流温度,单位为开尔文(K);T0为来流总温,单位为开尔文(K);γ为气体比热比;pi为第i个控制点压力;Mai为第i个控制点的马赫数。
可选地,声速流量系数Cd,取0.96~0.98。
S105、当所述空气质量流量大于所述最大空气质量流量Gmax时,减小喉道间隙尺寸,重复步骤S102至S104重新计算空气质量流量,直至所述空气质量流量小于所述最大空气质量流量Gmax
S20、根据所述各控制点的流动马赫数和所述喉道面积,计算各控制点的无粘流道间隙尺寸,具体包括:
根据第i个控制点的马赫数、喉道面积和气体比热比计算第i个控制点处的无粘流道面积,第i个控制点处的无粘流道面积σi,单位为平方米(m2),表达为:
Figure BDA0003430173720000101
其中,σ*为喉道面积,Mai为第i个控制点的马赫数,γ为气体比热比。
根据所述第i个控制点处的无粘流道面积和第i个控制点对应的天线罩外表面半径计算所述无粘流道间隙尺寸。所述第i个控制点对应的天线罩外表面半径是指天线罩上轴向位置与第i个控制点轴向位置相同的位置处的半径。
具体地,第i个控制点处的所述无粘流道间隙尺寸δi',单位为米(m),表达为:
Figure BDA0003430173720000102
其中,σi为第i个控制点处的无粘流道面积;Ri为第i个控制点对应的天线罩外表面半径,单位为米(m);σ*为喉道面积,Mai为第i个控制点的马赫数,γ为气体比热比。
S30、根据所述来流总压及各控制点的流动马赫数,计算各控制点的流道间隙边界层厚度,包括:
S301、根据第i个控制点的马赫数、来流总温和气体比热比计算第i个控制点处的气流温度;具体地,第i个控制点处的气流温度,表达为:
Figure BDA0003430173720000111
其中,Ti为第i个控制点处的气流温度,pi为第i个控制点压力,T0为来流总温,γ为气体比热比,Mai为第i个控制点的马赫数。
S302、根据所述第i个控制点处的气流温度、第i个控制点压力和理想气体常数计算第i个控制点处的气流密度;具体地,第i个控制点处的气流密度ρi,表达为:
Figure BDA0003430173720000112
其中,Ti为第i个控制点处的气流温度,pi为第i个控制点压力,R为理想气体常数,理想气体常数R=287J/(kg*K)。
S303、根据所述第i个控制点处的气流密度和气流温度计算第i个控制点处的气流雷诺数;具体地,第i个控制点处的气流雷诺数Rexi,表达为:
Figure BDA0003430173720000113
其中,Ti为第i个控制点处的气流温度,ρi为第i个控制点处的气流密度。
S304、根据所述第i个控制点处的气流雷诺数和第i个控制点与天线罩顶点的轴向距离,计算第i个控制点的流道间隙边界层厚度。具体地,所述第i个控制点的流道间隙边界层厚度δi”,表达为:
Figure BDA0003430173720000121
其中,xi为第i个控制点与天线罩顶点的轴向距离,单位为米(m),Rexi为第i个控制点处的气流雷诺数。
S40、根据所述各控制点的轴向位置、所述无粘流道间隙尺寸和所述流道间隙边界层厚度,计算各控制点坐标,将得到的所述控制点拟合成连续光滑曲线,生成喷管型线。
具体地,本发明中利用一维变截面可压缩流动方程确定无粘流道尺寸,利用管流边界层厚度公式对流道进行粘性修正。
具体地,所述根据所述各控制点的轴向位置、所述无粘流道间隙尺寸和所述流道间隙边界层厚度,计算各控制点坐标,包括:根据所述无粘流道间隙尺寸,以及所述流道间隙边界层厚度计算各控制点的流道间隙尺寸;更具体地,如图4所示,流道间隙尺寸中包括喷管表面的流道间隙边界层厚度和天线罩表面的流道间隙边界层厚度,所述各控制点的流道间隙尺寸,表达为:
δi=δi'+2δi”  (10)
其中,δi为第i个控制点的流道间隙尺寸,δi'为第i个控制点的无粘流道间隙尺寸,δi”为第i个控制点的流道间隙边界层厚度。
根据所述各控制点的轴向位置、所述流道间隙尺寸和所述控制点对应的天线罩外表面半径,计算各控制点坐标,所述各控制点坐标,表达为:
(xi,yi)=(xi,Rii)
其中,xi为第i个控制点与天线罩顶点的轴向距离,yi为第i个控制点与天线罩旋转轴间的纵向距离,Ri为第i个控制点对应的天线罩外表面半径,δi为第i个控制点的流道间隙尺寸,δi'为第i个控制点的无粘流道间隙尺寸,δi”为第i个控制点的流道间隙边界层厚度。
具体地,通过下文说明一种轴对称天线罩包罩试验喷管型线设计方法的工作流程:
(1)针对某天线罩,将流场喉道布置于天线罩根部的金属框上,距离天线罩顶点轴向距离为600mm,喉道间隙尺寸初始值为10mm。喉道面积为9424.8mm2,或9.4248×10-3m2。
(2)根据上述公式(1),调整来流总压以及天线罩母线控制点(母线上8~10个轴向位置)处的流动马赫数,使天线罩母线控制点处的压力满足试验考核要求。针对本实例中的天线罩,来流总压初始值设定为260kPa,如表1所示,在天线罩母线上x=0.1m~0.6m范围内布置了8个控制点,天线罩考核区域(天线罩非金属部分)内控制点对应的最大压力为239kPa,最小压力为143kPa,来流总压为260kPa时,满足控制点1~6处的流动马赫数均小于1,喉道位置达到声速,即马赫数等于1,保证包罩喷管内考核区域内气流为亚声速,控制点8处为超声速,因此,本实施例的计算过程中来流总压为260kPa。
表1
Figure BDA0003430173720000131
(3)根据天线罩测试要求来流总温取785K,总压取260kPa,Cd取0.96,根据上述公式(2)计算空气质量流量为3.69kg/s,小于电弧加热器所能提供的最大空气流量,因此喉道面积满足电弧加热器流量要求,喉道间隙尺寸取初始值10mm。
(4)根据喉道面积和表1中各控制点处的流动马赫数,使用上述公式(4)计算控制点处的无粘流道面积,再使用上述公式(5)计算无粘流道间隙尺寸,本实例中各控制点对应的无粘流道面积和无粘流道间隙尺寸如表2所示:
表2
Figure BDA0003430173720000141
(5)根据上述公式(6)~(9)计算边界层厚度,得到各控制点的边界层厚度参数。如表3所示,本实例中喉道处单边边界层厚度最厚为1.133mm,约为无粘流道间隙尺寸的11.3%。
表3
Figure BDA0003430173720000142
(6)如表4所示根据上述表2和表3中各控制点对应的所述无粘流道间隙尺寸和所述流道间隙边界层厚度,通过上述公式(10)计算控制点流道间隙尺寸设计值。本实例中,流道间隙尺寸最小为喉道处,为12.27mm。
表4
控制点编号 <![CDATA[x<sub>i</sub>(m)]]> <![CDATA[δ<sub>i</sub>(mm)]]>
1 0.1 41.83
2 0.15 30.41
3 0.2 23.88
4 0.3 18.51
5 0.4 15.49
6 0.5 13.75
7 0.6 12.27
8 0.655 17.07
(7)将表4中各控制点流道间隙尺寸加上对应位置的天线罩外表面半径,得到喷管型线上8个控制点坐标,具体如表5所示。
表5
控制点编号 <![CDATA[x<sub>i</sub>(mm)]]> <![CDATA[y<sub>i</sub>(mm)]]>
1 100 86.03
2 150 88.01
3 200 93.68
4 300 110.31
5 400 127.39
6 500 144.05
7 600 157.28
8 655 162.07
(8)将表5中型线各控制点拟合成连续、光滑的曲线,形成最终喷管型线。如图5所示,图5为根据本发明设计的某天线罩包罩喷管型线及压力仿真云图;图6为根据本发明设计的喷管实现的压力与考核要求对比,喷管可准确模拟天线罩表面压力分布。
与现有技术相比,本发明中一种轴对称天线罩包罩试验喷管型线设计方法,利用一维变截面可压缩流动和边界层理论,设计包罩喷管型线,在电弧加热器产生的气流条件下可准确模拟天线罩表面压力分布,对天线罩结构承载能力验证更为真实;通过天线罩外包覆一个随型包罩喷管,让电弧加热器产生的热气流在喷管和天线罩之间的缝隙内流动,可以在较小的加热功率和空气流量条件下,实现较高的压力和温度参数,能够显著节省气流流量和加热功率,能量利用率提升几十倍以上。
本领域技术人员可以理解,实现上述实施例方法的全部或部分流程,可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于计算机可读存储介质中。其中,所述计算机可读存储介质为磁盘、光盘、只读存储记忆体或随机存储记忆体等。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。

Claims (9)

1.一种轴对称天线罩包罩试验喷管型线设计方法,其特征在于,包括:
设置喉道间隙尺寸以及喷管上各控制点的轴向位置,计算喉道面积,并根据各控制点压力计算来流总压及各控制点的流动马赫数,其中,喉道为天线罩外表面与喷管内表面之间形成的流道的面积最小处;各控制点压力为天线罩上与各控制点轴向位置相同的位置处的考核压力;
根据所述各控制点的流动马赫数和所述喉道面积,计算各控制点的无粘流道间隙尺寸;
根据所述来流总压及各控制点的流动马赫数,计算各控制点的流道间隙边界层厚度;
根据所述各控制点的轴向位置、所述无粘流道间隙尺寸和所述流道间隙边界层厚度,计算各控制点坐标,将得到的所述控制点拟合成连续光滑曲线,生成喷管型线;
设置喉道间隙尺寸以及喷管上各控制点的轴向位置,计算喉道面积,并根据各控制点压力计算来流总压及各控制点的流动马赫数,包括:
根据待测天线罩尺寸和电弧加热器的质量流量设置喉道位置、喉道间隙尺寸初始值以及各控制点的轴向位置;
根据所述喉道间隙尺寸计算喉道面积;
根据所述各控制点压力,确定来流总压及各控制点的流动马赫数,并满足天线罩顶点和喉道之间控制点的马赫数小于1,喉道外远离顶点方向的控制点的马赫数大于1;
根据所述喉道面积及来流总压,计算空气质量流量,判断所述空气质量流量是否小于电弧加热器所能提供的最大空气质量流量Gmax
当所述空气质量流量大于所述最大空气质量流量Gmax时,减小喉道间隙尺寸,重新计算空气质量流量,直至所述空气质量流量小于所述最大空气质量流量Gmax
2.根据权利要求1所述的轴对称天线罩包罩试验喷管型线设计方法,其特征在于,
所述根据所述各控制点的轴向位置、所述无粘流道间隙尺寸和所述流道间隙边界层厚度,计算各控制点坐标,包括:
根据所述无粘流道间隙尺寸,以及所述流道间隙边界层厚度计算各控制点的流道间隙尺寸;
根据所述各控制点的轴向位置、所述流道间隙尺寸和所述控制点对应的天线罩外表面半径,计算各控制点坐标。
3.根据权利要求2所述的轴对称天线罩包罩试验喷管型线设计方法,其特征在于,
所述各控制点坐标,表达为:
(xi,yi)=(xi,Rii)
δi=δi'+2δi
其中,xi为第i个控制点与天线罩顶点的轴向距离,yi为第i个控制点与天线罩旋转轴间的纵向距离,Ri为第i个控制点对应的天线罩外表面半径,δi为第i个控制点的流道间隙尺寸,δi'为第i个控制点的无粘流道间隙尺寸,δi”为第i个控制点的流道间隙边界层厚度。
4.根据权利要求3所述的轴对称天线罩包罩试验喷管型线设计方法,其特征在于,
所述根据所述各控制点的流动马赫数和所述喉道面积,计算各控制点无粘流道间隙尺寸,包括:
根据第i个控制点的马赫数、喉道面积和气体比热比计算第i个控制点处的无粘流道面积;
根据所述第i个控制点处的无粘流道面积和第i个控制点对应的天线罩外表面半径计算所述无粘流道间隙尺寸。
5.根据权利要求4所述的轴对称天线罩包罩试验喷管型线设计方法,其特征在于,
第i个控制点处的所述无粘流道间隙尺寸δi',表达为:
Figure FDA0004129658330000031
Figure FDA0004129658330000032
其中,σi为第i个控制点处的无粘流道面积,Ri为第i个控制点对应的天线罩外表面半径,σ*为喉道面积,Mai为第i个控制点的马赫数,γ为气体比热比。
6.根据权利要求3所述的轴对称天线罩包罩试验喷管型线设计方法,其特征在于,
所述根据所述来流总压及各控制点的流动马赫数,计算各控制点的流道间隙边界层厚度,包括:
根据第i个控制点的马赫数、来流总温和气体比热比计算第i个控制点处的气流温度;
根据所述第i个控制点处的气流温度、第i个控制点压力和理想气体常数计算第i个控制点处的气流密度;
根据所述第i个控制点处的气流密度和气流温度计算第i个控制点处的气流雷诺数;
根据所述第i个控制点处的气流雷诺数和第i个控制点与天线罩顶点的轴向距离,计算第i个控制点的流道间隙边界层厚度。
7.根据权利要求6所述的轴对称天线罩包罩试验喷管型线设计方法,其特征在于,
所述第i个控制点的流道间隙边界层厚度δi”,表达为:
Figure FDA0004129658330000041
Figure FDA0004129658330000042
Figure FDA0004129658330000043
Figure FDA0004129658330000044
其中,xi为第i个控制点与天线罩顶点的轴向距离,Rexi为第i个控制点处的气流雷诺数,Ti为第i个控制点处的气流温度,ρi为第i个控制点处的气流密度,pi为第i个控制点压力,R为理想气体常数,T0为来流总温,γ为气体比热比,Mai为第i个控制点的马赫数。
8.根据权利要求1所述的轴对称天线罩包罩试验喷管型线设计方法,其特征在于,
所述根据所述各控制点压力,确定来流总压及各控制点的流动马赫数包括:根据如下公式确定来流总压及各控制点的流动马赫数:
Figure FDA0004129658330000045
其中,P0为来流总压,γ为气体比热比,pi为第i个控制点压力,Mai为第i个控制点的马赫数。
9.根据权利要求8所述的轴对称天线罩包罩试验喷管型线设计方法,其特征在于,
所述空气质量流量G,表达为:
Figure FDA0004129658330000051
Figure FDA0004129658330000052
其中,P0为来流总压,σ*为喉道面积,Cd为声速流量系数,T*为喉道处气流温度,T0为来流总温,γ为气体比热比。
CN202111594629.XA 2021-12-23 2021-12-23 一种轴对称天线罩包罩试验喷管型线设计方法 Active CN114278460B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111594629.XA CN114278460B (zh) 2021-12-23 2021-12-23 一种轴对称天线罩包罩试验喷管型线设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111594629.XA CN114278460B (zh) 2021-12-23 2021-12-23 一种轴对称天线罩包罩试验喷管型线设计方法

Publications (2)

Publication Number Publication Date
CN114278460A CN114278460A (zh) 2022-04-05
CN114278460B true CN114278460B (zh) 2023-05-12

Family

ID=80875370

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111594629.XA Active CN114278460B (zh) 2021-12-23 2021-12-23 一种轴对称天线罩包罩试验喷管型线设计方法

Country Status (1)

Country Link
CN (1) CN114278460B (zh)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3258236B1 (en) * 2016-06-14 2019-10-16 Haute école du paysage, d'ingénierie et d'architecture de Genève Wind generation means and wind test facility comprising the same
US20190265125A1 (en) * 2018-02-23 2019-08-29 Global Reach Aerospace LLC Large test area compressed air wind tunnel
CN208128550U (zh) * 2018-03-02 2018-11-20 中国航空工业集团公司沈阳飞机设计研究所 一种石英灯辐射加热器
CN210487644U (zh) * 2019-05-27 2020-05-08 中国航天空气动力技术研究院 电弧风洞烧蚀透波联合试验装置
CN110712764A (zh) * 2019-10-21 2020-01-21 中国航天空气动力技术研究院 一种用于高焓条件下的亚声速包罩烧蚀试验装置
CN112945507B (zh) * 2021-02-03 2023-05-19 中国空气动力研究与发展中心高速空气动力研究所 一种高超声速风洞轴对称喷管

Also Published As

Publication number Publication date
CN114278460A (zh) 2022-04-05

Similar Documents

Publication Publication Date Title
Louw et al. Investigation of the flow field in the vicinity of an axial flow fan during low flow rates
Wilkinson et al. Performance testing of an axial flow fan designed for air-cooled heat exchanger applications
CN115879396B (zh) 高空模拟试车台进气前室流程化一维气动设计方法
CN109815624A (zh) 一种考虑进气总压畸变影响的压气机稳定边界判断方法
CN109583147A (zh) 一种离心叶轮预旋转盘心梯度模拟件设计方法
Lee et al. Performance assessment of various fan ribs inside a centrifugal blower
Brossman et al. Sensitivity of multistage compressor performance to inlet boundary conditions
Zhao et al. Film cooling of showerhead holes from the twisted leading edge of a gas turbine blade: Complex mainstream characteristics and reasonable angle arrangement
CN114278460B (zh) 一种轴对称天线罩包罩试验喷管型线设计方法
Xu et al. Effect of rotor axial blade loading distribution on compressor stability
Li et al. Effect of non-axisymmetric endwall profiling on heat transfer and film cooling effectiveness of a transonic rotor blade
CN111780949A (zh) 基于cfd分析的高速进气道前体风洞实验总压修正方法
Yang et al. Prediction and optimization of aerodynamic noise in an automotive air conditioning centrifugal fan
CN112883575A (zh) 一种考虑表面粗糙度的叶轮机械边界层转捩模型修正方法
Zhao et al. A novel clocking effect between inlet bend and volute in an automotive turbocharging system
Hu et al. Performance prediction of transonic axial compressor based on streamline curvature method
Dong et al. Experimental investigation on SPS casing treatment with bias flow
Scribben et al. Effectiveness of a serpentine inlet duct flow control technique at design and off-design simulated flight conditions
Zenz et al. Aerodynamical and aeroelastic investigations of a riblet design applied on the surface of turbine exit guide vanes of a low pressure turbine
CN114048697B (zh) 一种飞机试验气流组织设计方法
CN113959691A (zh) 一种中介机匣气动性能试验件的设计方法
Li et al. Experimental investigations on cooling hole diameters of an impingement-effusion cooling system
Abbasi et al. Investigation of inlet distortion effects on axial compressor performance based on streamline curvature method
Guo et al. Design of a continuously variable Mach-number nozzle
Zhou et al. Calculation of Cooling Fan Blade Deformation and Aerodynamic Performance Based on Fluid-Structure Model

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