CN106644362B - 一种超音速空腔流动的马赫数敏感性分析方法 - Google Patents
一种超音速空腔流动的马赫数敏感性分析方法 Download PDFInfo
- Publication number
- CN106644362B CN106644362B CN201611061983.5A CN201611061983A CN106644362B CN 106644362 B CN106644362 B CN 106644362B CN 201611061983 A CN201611061983 A CN 201611061983A CN 106644362 B CN106644362 B CN 106644362B
- Authority
- CN
- China
- Prior art keywords
- mach number
- attack
- sensitivity analysis
- angle
- cavity
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种超音速空腔流动的马赫数敏感性分析方法,该方案通过调节模型迎角,在平板前缘形成不同角度的膨胀波,气流经过膨胀波后马赫数将升高,通过对比迎角变化前后的脉动压力数据,开展马赫数的敏感性分析。能够实现在超声速条件下开展马赫数敏感性分析,解决了风洞敏感性分析试验的关键问题。
Description
技术领域
本发明涉及的是实验流体力学领域,尤其是一种适用于超音速空腔流动的马赫数敏感性分析方法。
背景技术
随着我国武器装备的自主创新发展,越来越多的非常规气动布局在先进战机上得到了应用。然而,新气动布局方案在提升飞行器机敏性、隐身效果的同时,也带来了许多新问题,如大攻角飞行时涡破裂导致垂尾抖振、内埋武器舱开舱时产生的剧烈压力脉动等问题。这些问题导致的后果之一就是剧烈的动载荷作用下,垂尾蒙皮、内埋武器舱薄壁结构表面裂纹加速形成和扩展,严重威胁飞行器的飞行安全。为了解决上述动载荷作用下的结构疲劳失效问题,飞机设计单位已经构建了较为成熟的结构设计和疲劳分析方法,存在的主要问题是缺乏有效精确的动载荷数据作为飞机结构设计和疲劳分析的输入。
现阶段,设计人员主要通过风洞试验来获取动载荷数据。然而,由于动载荷预报的风洞试验相似模拟理论还不完善,动载荷数据的不确定性分析方法尚未建立,导致当前动载荷预测偏差较大。为了安全起见,设计人员往往采用较高的安全系数,导致结构重量增加,严重制约飞行器综合性能的提升。
构建相似模拟理论的关键,在于找到影响主要物理现象的相似参数。郑哲敏、谈庆明指出,对于不能完全模拟的问题,采用局部模拟是可行的,即对控制主要物理现象和过程的相似参数进行完全模拟,而放松对其他次要因素的要求。周济福、李家春提出了一种计算相似参数的敏感度因子的方法,指出敏感度因子可用于表征相似参数的重要程度,并将该方法应用于油藏多孔介质渗流研究,获得了与理论分析一致的相似准则。对于复杂强非线性流动,理论分析构建的相似准则往往考虑到所有的影响因素,在试验模拟时几乎难以实现。而周济福、李家春提出的计算敏感度因子的方法,有助于辨识各相似参数的重要程度,从而构建基于主要相似参数的风洞试验相似准则。
开展流场参数的敏感度分析,不仅有助于风洞试验相似准则的构建,同时敏感性分析数据还可用于静载荷、动载荷等试验数据的不确定性分析,解决目前动载荷预报试验技术中相似准则构建和预测精度估计的两大难题,为飞行器设计提供准确的动载荷数据。此外,基于模型几何参数的敏感性分析数据可用于气动外形优化。基于流动控制参数的敏感性分析结果,不仅对于控制参数的选取和优化具有重要意义,在开展流场主动控制研究、揭示控制机理等方面也将发挥重要作用。
目前,敏感性分析方法在总体参数设计、工程结构设计、石油开采等方面应用较为广泛。罗鹏程等开展了武器装备敏感性分析研究,指出敏感性分析结果中可以提炼出哪些装备是重要武器装备的结论。金镭等通过敏感性分析找出了性能较差飞机的各个组成部分与先进飞机之间的差距,定量的给出了这些部分具体能够提升的空间,为战斗机改进改型和作战使用提供理论指导。唐冕针对大跨度自锚式悬索桥结构,基于桥梁多振型耦合的气弹理论,建立了气动参数敏感性分析的理论和方法。白玉湖等通过相似参数的数值实验,定量分析了各个相似参数对于水驱油计算结果的影响程度,对各个相似参数的敏感因子进行了比较,从而确定了各个相似参数的主次关系,为实现部分相似提供了理论依据和设计原则。
敏感性分析方法在诸多工程研究领域应用十分广泛,但是在空气动力研究方面的应用较少,仅有西北工业大学的徐林程等采用数值计算的方法开展了翼型的敏感性分析研究。数值计算在对单目标的稳态流动进行敏感性分析方面具有一定的优势,但是对于多目标的非定常流动,则计算量显著上升,计算结果的收敛性也急剧下降。
随着风洞流场品质和控制精度的不断提升,目前,风洞试验基本实现了对流场参数的连续高精度调节,这对于开展流场参数的敏感性研究十分有利。然而,在超音速风洞试验中,雷诺数、模型姿态等参数都能实现连续变化,但是由于马赫数由喷管型面决定,不能像亚声速一样实现连续变化。超声速条件下,如何开展马赫数敏感性分析,是风洞敏感性分析试验亟待解决的关键问题。
发明内容
本发明的目的,就是针对现有技术所存在的不足,而提供一种超音速空腔流动的马赫数敏感性分析方法,该方案通过调节模型迎角,在平板前缘形成不同角度的膨胀波,气流经过膨胀波后马赫数将升高,通过对比迎角变化前后的脉动压力数据,开展马赫数的敏感性分析。
本方案是通过如下技术措施来实现的:
一种超音速空腔流动的马赫数敏感性分析方法,包括有以下步骤:
a.使平板起始迎角为0度,设置每步迎角增加Δα,依次测量3步;
b.在空腔前缘安装总压耙板和静压测孔,测量获得在不同迎角下的空腔入口马赫数;
c.控制每步间的马赫数差量小于0.01,若步与步之间马赫数差量大于0.01,则下调每步的迎角差量,直至每步间马赫数差量在0.01以内;
d.根据步骤c确定的迎角差量Δα,进行0度、Δα和2Δα迎角下的空腔脉动压力测量试验;
e.采用差分方法,计算不同测点的脉动压力系数,开展空腔流动马赫数敏感性分析。
作为本方案的优选:步骤b中测量空腔入口马赫数的方法为:根据雷列公式(1),计算得到i测点处马赫数Mi,选择最外层5个测点马赫数的平均值作为空腔入口马赫数M;
式中,Pi为第i个总压测点压力值,Ps为静压测孔的压力值。
作为本方案的优选:步骤e中,马赫数敏感性分析方法为:
0度、Δα和2Δα迎角下,空腔入口马赫数依次记为M0,M1,M2,测点K处脉动压力系数依次记为Cp0,Cp1,Cp2,则测点K处的马赫数敏感性导数可用式(3)计算得到:
本方案的有益效果可根据对上述方案的叙述得知,由于在该方案通过调节模型迎角,在平板前缘形成不同角度的膨胀波,气流经过膨胀波后马赫数将升高,通过对比迎角变化前后的脉动压力数据,开展马赫数的敏感性分析。
由此可见,本发明与现有技术相比,具有实质性特点和进步,其实施的有益效果也是显而易见的。
附图说明
图1为本方案实施方式的示意图。
图中,1为平板,2为空腔,3为静压测孔,4为总压耙板。
具体实施方式
本说明书中公开的所有特征,或公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合。
本说明书(包括任何附加权利要求、摘要和附图)中公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换。即,除非特别叙述,每个特征只是一系列等效或类似特征中的一个例子而已。
本方案包括有以下步骤:
a.使平板起始迎角为0度,设置每步迎角增加Δα,依次测量3步;
b.在空腔前缘安装总压耙板和静压测孔,测量获得在不同迎角下的空腔入口马赫数;
c.控制每步间的马赫数差量小于0.01,若步与步之间马赫数差量大于0.01,则下调每步的迎角差量,直至每步间马赫数差量在0.01以内;
d.根据步骤c确定的迎角差量Δα,进行0度、Δα和2Δα迎角下的空腔脉动压力测量试验;
e.采用差分方法,计算不同测点的脉动压力系数,开展空腔流动马赫数敏感性分析。
步骤b中测量空腔入口马赫数的方法为:根据雷列公式(1),计算得到i测点处马赫数Mi,选择最外层5个测点马赫数的平均值作为空腔入口马赫数M;
式中,Pi为第i个总压测点压力值,Ps为静压测孔的压力值。
步骤e中,马赫数敏感性分析方法为:
0度、Δα和2Δα迎角下,空腔入口马赫数依次记为M0,M1,M2,测点K处脉动压力系数依次记为Cp0,Cp1,Cp2,则测点K处的马赫数敏感性导数可用式(3)计算得到:
采用本方法,通过调节模型迎角,在平板前缘形成不同角度的膨胀波,气流经过膨胀波后马赫数将升高,通过对比迎角变化前后的脉动压力数据,开展马赫数的敏感性分析。能够有效解决在超声速条件下,如何开展的马赫数敏感性分析的问题,具有显著的技术贡献。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合,以及披露的任一新的方法或过程的步骤或任何新的组合。
Claims (3)
1.一种超音速空腔流动的马赫数敏感性分析方法,其特征是:包括有以下步骤:
a.使平板起始迎角为0度,设置每步迎角增加Δα,依次测量3步;
b.在空腔前缘安装总压耙板和静压测孔,测量获得在不同迎角下的空腔入口马赫数;
c.控制每步间的马赫数差量小于0.01,若步与步之间马赫数差量大于0.01,则下调每步的迎角差量,直至每步间马赫数差量在0.01以内;
d.根据步骤c确定的迎角差量Δα,进行0度、Δα和2Δα迎角下的空腔脉动压力测量试验;
e.采用差分方法,计算不同测点的脉动压力系数,开展空腔流动马赫数敏感性分析。
2.根据权利要求1所述的方法,其特征是:所述步骤b中测量空腔入口马赫数的方法为:根据雷列公式(1),计算得到i测点处马赫数Mi,选择最外层5个测点马赫数的平均值作为空腔入口马赫数M;
式中,Pi为第i个总压测点压力值,Ps为静压测孔测的压力值。
3.根据权利要求1所述的方法,其特征是:所述步骤e中,马赫数敏感性分析方法为:
0度、Δα和2Δα迎角下,空腔入口马赫数依次记为M0,M1,M2,测点K处脉动压力系数依次记为Cp0,Cp1,Cp2,则测点K处的马赫数敏感性导数可用式(3)计算得到:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611061983.5A CN106644362B (zh) | 2016-11-24 | 2016-11-24 | 一种超音速空腔流动的马赫数敏感性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611061983.5A CN106644362B (zh) | 2016-11-24 | 2016-11-24 | 一种超音速空腔流动的马赫数敏感性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106644362A CN106644362A (zh) | 2017-05-10 |
CN106644362B true CN106644362B (zh) | 2018-11-02 |
Family
ID=58811621
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611061983.5A Active CN106644362B (zh) | 2016-11-24 | 2016-11-24 | 一种超音速空腔流动的马赫数敏感性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106644362B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107917793B (zh) * | 2017-11-14 | 2020-09-25 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种实验用多排测压耙装置 |
CN108491572B (zh) * | 2018-02-11 | 2021-11-23 | 合肥市太泽透平技术有限公司 | 一种基于两面流理论的透平机械二维流动仿真方法 |
CN113029508B (zh) * | 2021-03-24 | 2023-06-02 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种用于风洞模型底部压力测量的微型组合式压力传感器 |
CN115200826B (zh) * | 2022-09-16 | 2022-11-15 | 中国空气动力研究与发展中心高速空气动力研究所 | 背负式进气小展弦比飞翼布局的背风涡破裂迎角修正方法 |
CN117076956B (zh) * | 2023-10-16 | 2024-01-26 | 西安石油大学 | 缝洞油藏物理模型相似准数优化方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5612492A (en) * | 1995-06-07 | 1997-03-18 | Mcdonnell Douglas Corporation | Formulations and method of use of pressure sensitive paint |
JP2005157816A (ja) * | 2003-11-26 | 2005-06-16 | Japan Aerospace Exploration Agency | 磁場のフィードバック制御を行う磁力支持装置 |
CN103150423A (zh) * | 2013-02-05 | 2013-06-12 | 中国人民解放军国防科学技术大学 | 轴对称连续变马赫数喷管及其壁面确定方法 |
CN103969022A (zh) * | 2014-05-23 | 2014-08-06 | 厦门大学 | 一种高超声速风洞湍流度间接测量方法 |
CN105157948A (zh) * | 2015-09-14 | 2015-12-16 | 南京航空航天大学 | 一种适用于超声速/高超声速流道的流量测试系统及测试方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6696690B2 (en) * | 2001-12-17 | 2004-02-24 | The Boeing Company | Method and apparatus to correct for the temperature sensitivity of pressure sensitive paint |
JP5747343B2 (ja) * | 2010-12-14 | 2015-07-15 | 国立研究開発法人宇宙航空研究開発機構 | 超音速航空機の実機相当レイノルズ数における自然層流翼の設計方法 |
-
2016
- 2016-11-24 CN CN201611061983.5A patent/CN106644362B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5612492A (en) * | 1995-06-07 | 1997-03-18 | Mcdonnell Douglas Corporation | Formulations and method of use of pressure sensitive paint |
JP2005157816A (ja) * | 2003-11-26 | 2005-06-16 | Japan Aerospace Exploration Agency | 磁場のフィードバック制御を行う磁力支持装置 |
CN103150423A (zh) * | 2013-02-05 | 2013-06-12 | 中国人民解放军国防科学技术大学 | 轴对称连续变马赫数喷管及其壁面确定方法 |
CN103969022A (zh) * | 2014-05-23 | 2014-08-06 | 厦门大学 | 一种高超声速风洞湍流度间接测量方法 |
CN105157948A (zh) * | 2015-09-14 | 2015-12-16 | 南京航空航天大学 | 一种适用于超声速/高超声速流道的流量测试系统及测试方法 |
Non-Patent Citations (2)
Title |
---|
基于全局灵敏度分析的侧向气动导数不确定性对侧向飞行载荷的影响;阮文斌 等;《航空学报》;20160625;第37卷(第6期);第1827-1832页 * |
翼型风洞试验中不确定性分析的自动微分方法;徐林程 等;《航空学报》;20140825;第35卷(第8期);第2102-2111页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106644362A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106644362B (zh) | 一种超音速空腔流动的马赫数敏感性分析方法 | |
CN106650046B (zh) | 一种舰船空气流场的非定常特性获取方法 | |
CN104298826A (zh) | 一种反推力状态下的航空发动机气动稳定性预测与评估方法 | |
El Khoury et al. | Crossflow past a prolate spheroid at Reynolds number of 10000 | |
Cenko | Store separation lessons learned during the last 30 years | |
Ramanujam et al. | Improving airfoil lift prediction | |
CN110702363A (zh) | 针对雷诺数影响的高空螺旋桨风洞试验数据修正方法 | |
CN109033548A (zh) | 一种计算槽壁边界条件主要系数的拟合方法 | |
Bobke et al. | Large-eddy simulations of adverse pressure gradient turbulent boundary layers | |
Aiguo et al. | Experimental study of a controlled variable double-baffle distortion generator engine test rig | |
Chen et al. | Aerodynamic optimization design for high pressure turbines based on the adjoint approach | |
Hashimoto et al. | Transonic wind tunnel simulation with porous wall and support devices | |
Valger et al. | Structure of turbulent separated flow in the neighborhood of a plate-mounted prism of square section | |
Cenko | Store separation overview | |
CN110207937B (zh) | 一种考虑粗糙效应的飞行器湍流确定方法及系统 | |
Mahmuddina | The effect of flat plate theory assumption in post-stall lift and drag coefficients extrapolation with Viterna method | |
Yu et al. | Mission segment division of the whole aeroengine loading spectrum based on flight actions | |
Murty et al. | Numerical simulation of nozzle flow field with jet vane thrust vector control | |
Tucker et al. | Laminar flow control flight experiment design and execution | |
Benderskii et al. | The use of the RANS/ILES method to study the influence of coflow wind on the flow in a hot, nonisobaric, supersonic airdrome jet during its interaction with the jet blast deflector | |
Lee et al. | Computational simulations of transitional flows around turbulence stimulators at low speeds | |
Halila et al. | Inflow turbulence effects on transition prediction using a correlation-based transition model | |
Huang et al. | Experimental and numerical studies of geometry effects on UCAV's aerodynamics | |
Zamiri et al. | Numerical evaluation of wind direction effects on the turbulence aerodynamics of a ship airwake | |
Kumar et al. | Modelling of cascade fin aerodynamics near stall using Kirchhoff’s steady-state stall 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 |