CN105989205A - 飞行器表面脉动压力的确定方法 - Google Patents
飞行器表面脉动压力的确定方法 Download PDFInfo
- Publication number
- CN105989205A CN105989205A CN201510076087.5A CN201510076087A CN105989205A CN 105989205 A CN105989205 A CN 105989205A CN 201510076087 A CN201510076087 A CN 201510076087A CN 105989205 A CN105989205 A CN 105989205A
- Authority
- CN
- China
- Prior art keywords
- fluctuation pressure
- pulsating pressure
- pressure
- aircraft surface
- aircraft
- 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.)
- Granted
Links
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供一种飞行器表面脉动压力的确定方法,包括:对飞行器进行几何建模并划分单元,并计算表面脉动压力;按照预设倍数对所述单元进行N次加密,并在每一次加密后计算表面脉动压力,其中N为大于或等于1的整数;检测第N次加密后计算得到的脉动压力与上一次计算得到的脉动压力是否收敛在预设范围内;如果是,则将第N次加密后计算得到的脉动压力与典型风洞试验结果进行比较,根据比较结果对模型进行局部加密,之后计算获得脉动压力基础值;根据典型风洞试验结果与所述脉动压力基础值建立脉动压力修正曲线,获得修正差值,根据所述修正差值计算得到对应的飞行器表面脉动压力结果值;该方法能够有效提高飞行器表面脉动压力测量的准确性和普适性。
Description
技术领域
本发明涉及飞行器环境工程技术领域,尤其涉及一种飞行器表面脉动压力的确定方法。
背景技术
飞行器在飞行过程中会受到表面脉动压力的作用而产生振动,而飞行器振动是导致飞行器故障的主要因素,准确确定脉动压力对于开展飞行器抗振动设计和试验验证,保证飞行安全具有重要意义。目前脉动压力预测仍以经验公式为主,工程预测精度低,无法描述流场细节,无法预测压力随时间的响应,主要用来预估脉动压力的变化趋势。脉动压力的数值计算仍以雷诺平均法(RANS)为主,RANS方法可以用现有计算机资源实现高雷诺数复杂流动的数值模拟,但是雷诺平均法只能提供湍流的平均信息,且湍流模型没有普适性,所以计算准确性较差,而另一种数值计算方法——直接数值模拟,由于计算机资源的限制,主要用于湍流理论的研究,不适合用于工程应用。对飞行器表面脉动压力进行精细研究的风洞试验因本底噪声较高,头部区域难以测得有效数据,同时,受试验能力或试验成本的限制,风洞试验只能进行少数几个角度、速度状态的试验,所能得到的状态结果有限。因此,如何用较小的代价实现飞行器表面脉动压力计算成为现在亟待需要解决的问题。
发明内容
在下文中给出关于本发明的简要概述,以便提供关于本发明的某些方面的基本理解。应当理解,这个概述并不是关于本发明的穷举性概述。它并不是意图确定本发明的关键或重要部分,也不是意图限定本发明的范围。其目的仅仅是以简化的形式给出某些概念,以此作为稍后论述的更详细描述的前序。
为解决上述问题,本发明提出一种飞行器表面脉动压力的确定方法,包括:
对飞行器进行几何建模并划分单元,并计算表面脉动压力;
按照预设倍数对所述单元进行N次加密,并在每一次加密后计算表面脉动压力,其中N为大于或等于1的整数;
检测第N次加密后计算得到的脉动压力与上一次计算得到的脉动压力是否收敛在预设范围内;
如果是,则将第N次加密后计算得到的脉动压力与典型风洞试验结果进行比较,根据比较结果对模型进行局部加密,之后计算获得脉动压力基础值;
根据典型风洞试验结果与所述脉动压力基础值建立脉动压力修正曲线,获得修正差值,根据所述修正差值计算得到对应的飞行器表面脉动压力结果值。
本发明提供的飞行器表面脉动压力的确定方法,能够有效提高飞行器表面脉动压力测量的准确性和普适性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的飞行器表面脉动压力的确定方法一种实施例的流程图。
图2为本发明提供的飞行器表面脉动压力的确定方法中对飞行器进行几何建模和单元划分的示意图。
图3为本发明提供的飞行器表面脉动压力的确定方法中对飞行器模型单元进行加密的示意图。
图4为为本发明提供的飞行器表面脉动压力的确定方法中脉动压力修正曲线一种实施例的示意图。
具体实施方式
下面参照附图来说明本发明的实施例。在本发明的一个附图或一种实施方式中描述的元素和特征可以与一个或者更多个其他附图或实施方式中示出的元素和特征相结合。应当注意,为了清楚目的,附图和说明中省略了与本发明无关的、本领域普通技术人员已知的部件和处理的表示和描述。
参考图1,本实施例提供一种飞行器表面脉动压力的确定方法,包括:
步骤S101,对飞行器进行几何建模并划分单元,并计算表面脉动压力;
步骤S102,按照预设倍数对所述单元进行N次加密,并在每一次加密后计算表面脉动压力,其中N为大于或等于1的整数;
步骤S103,检测第N次加密后计算得到的脉动压力与上一次计算得到的脉动压力是否收敛在预设范围内;
步骤S104,如果是,则将第N次加密后计算得到的脉动压力与典型风洞试验结果进行比较,根据比较结果对模型进行局部加密,之后计算获得脉动压力基础值;
步骤S105,根据典型风洞试验结果与所述脉动压力基础值建立脉动压力修正曲线,获得修正差值,根据所述修正差值计算得到对应的飞行器表面脉动压力结果值。
本实施例提供的飞行器表面脉动压力的确定方法,能够有效提高飞行器表面脉动压力测量的准确性和普适性。
具体地,参考图2,首先对行器进行几何建模并划分单元,并采用大涡模拟方法进行第一轮的脉动压力计算。
通过大涡模拟方法计算脉动压力是本领域的公知技术,包括:计算亚格子粘性系数;确定空间离散格式(迎风面推荐AUSM类格式);确定时间推进方法(由于大涡数值模拟是非定常计算,为了保证计算的精度,常常采用的是双时间法的LU-SGS。双时间法在隐式时间离散的公式中加入伪时间项,并进行与单时间法相似的推导);确定脉动压力时间步长;计算脉动压力。
参考图3,之后在第一轮计算模型的基础上,按照预设倍数为划分的单元进行N次加密,N为大于或等于1的整数,在每一次加密后均采用大涡模拟方法进行脉动压力的计算,并将当前得到的脉动压力与上一次加密计算得到的脉动压力进行比较,直到二者的差值收敛在预设范围内,作为一种优选的实施方式,预设倍数为2,预设范围为15%。
当二者的差值收敛在预设范围内,则将当次得到的脉动压力与典型风洞试验结果进行比较,根据比较结果对模型进行局部加密,具体地,对比较结果误差大于20%的局部进行加密,直到比较结果误差小于20%,并在模型几何外形突变的地方和计算结果变化剧烈的地方进行加密。
局部加密之后再进行脉动压力的计算,获得脉动压力基础值。
进一步地,参考图4,根据典型风洞试验结果与脉动压力基础值建立脉动压力修正曲线,获得相应的修正差值,该修正差值为所述典型风洞试验结果与所述脉动压力基础值之间的差值,将该修正差值与脉动压力基础值进行叠加,即可以获得飞行器表面脉动压力结果值。
典型风洞试验是将飞行器放置在风洞中,通过改变飞行器迎风角度、空气来流速度来模拟飞行器不同攻角、不同速度的飞行情况,并在飞行器表面安装脉动压力传感器进行脉动压力测量,从而获得典型飞行状态的脉动压力值。
由于典型风洞试验对应的状态比较少,在对脉动压力基础值进行修正时,还需在典型风洞试验脉动压力基础上对其他速度状态进行插值,得到其他速度状态的脉动压力修正差值。
虽然已经详细说明了本发明及其优点,但是应当理解在不超出由所附的权利要求所限定的本发明的精神和范围的情况下可以进行各种改变、替代和变换。而且,本申请的范围不仅限于说明书所描述的过程、设备、手段、方法和步骤的具体实施例。本领域内的普通技术人员从本发明的公开内容将容易理解,根据本发明可以使用执行与在此所述的相应实施例基本相同的功能或者获得与其基本相同的结果的、现有和将来要被开发的过程、设备、手段、方法或者步骤。因此,所附的权利要求旨在它们的范围内包括这样的过程、设备、手段、方法或者步骤。
Claims (7)
1.一种飞行器表面脉动压力的确定方法,其特征在于,包括:
对飞行器进行几何建模并划分单元,并计算表面脉动压力;
按照预设倍数对所述单元进行N次加密,并在每一次加密后计算表面脉动压力,其中N为大于或等于1的整数;
检测第N次加密后计算得到的脉动压力与上一次计算得到的脉动压力是否收敛在预设范围内;
如果是,则将第N次加密后计算得到的脉动压力与典型风洞试验结果进行比较,根据比较结果对模型进行局部加密,之后计算获得脉动压力基础值;
根据典型风洞试验结果与所述脉动压力基础值建立脉动压力修正曲线,获得修正差值,根据所述修正差值计算得到对应的飞行器表面脉动压力结果值。
2.根据权利要求1所述的飞行器表面脉动压力的确定方法,其特征在于,采用大涡模拟方法计算所述表面脉动压力。
3.根据权利要求1所述的飞行器表面脉动压力的确定方法,其特征在于,所述预设倍数为2,所述预设范围为15%。
4.根据权利要求1所述的飞行器表面脉动压力的确定方法,其特征在于,根据比较结果对模型进行局部加密,包括:
对比较结果误差大于20%的局部进行加密,直到比较结果误差小于20%。
5.根据权利要求4所述的飞行器表面脉动压力的确定方法,其特征在于,根据比较结果对模型进行局部加密之后,还包括:
在所述模型几何外形突变的地方进行加密;
在计算结果变化剧烈的地方进行加密。
6.根据权利要求1所述的飞行器表面脉动压力的确定方法,其特征在于,所述修正差值为所述典型风洞试验结果与所述脉动压力基础值之间的差值。
7.根据权利要求1所述的飞行器表面脉动压力的确定方法,其特征在于,所述飞行器表面脉动压力结果值为所述脉动压力基础值和所述修正差值的叠加。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510076087.5A CN105989205B (zh) | 2015-02-13 | 2015-02-13 | 飞行器表面脉动压力的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510076087.5A CN105989205B (zh) | 2015-02-13 | 2015-02-13 | 飞行器表面脉动压力的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105989205A true CN105989205A (zh) | 2016-10-05 |
CN105989205B CN105989205B (zh) | 2019-04-23 |
Family
ID=57042040
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510076087.5A Active CN105989205B (zh) | 2015-02-13 | 2015-02-13 | 飞行器表面脉动压力的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105989205B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107657109A (zh) * | 2017-09-22 | 2018-02-02 | 哈尔滨工业大学 | 一种超声速飞行器跨声速段脉动压力频域识别方法 |
CN109341945A (zh) * | 2018-11-23 | 2019-02-15 | 中国航空工业集团公司北京长城计量测试技术研究所 | 基于双管腔的脉动压力测量修正方法 |
CN111504595A (zh) * | 2020-03-31 | 2020-08-07 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种高速列车模型在声学风洞内的表面声载荷测量方法 |
CN112100895A (zh) * | 2020-10-28 | 2020-12-18 | 北京机电工程研究所 | 一种针对螺栓连接组合结构的分层模型修正方法和系统 |
CN112446091A (zh) * | 2020-11-19 | 2021-03-05 | 中国航天空气动力技术研究院 | 一种基于人工神经网络的脉动压力预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1125152A (en) * | 1965-01-19 | 1968-08-28 | Marconi Co Ltd | Improvements in or relating to sub-aqueous pressure wave sonar systems |
CN101881287A (zh) * | 2010-06-03 | 2010-11-10 | 西北工业大学 | 飞行器高压液压脉动试验系统 |
-
2015
- 2015-02-13 CN CN201510076087.5A patent/CN105989205B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1125152A (en) * | 1965-01-19 | 1968-08-28 | Marconi Co Ltd | Improvements in or relating to sub-aqueous pressure wave sonar systems |
CN101881287A (zh) * | 2010-06-03 | 2010-11-10 | 西北工业大学 | 飞行器高压液压脉动试验系统 |
Non-Patent Citations (1)
Title |
---|
刘振皓: "航天飞行器脉动压力数值计算方法综述", 《强度与环境》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107657109A (zh) * | 2017-09-22 | 2018-02-02 | 哈尔滨工业大学 | 一种超声速飞行器跨声速段脉动压力频域识别方法 |
CN107657109B (zh) * | 2017-09-22 | 2020-02-11 | 哈尔滨工业大学 | 一种超声速飞行器跨声速段脉动压力频域识别方法 |
CN109341945A (zh) * | 2018-11-23 | 2019-02-15 | 中国航空工业集团公司北京长城计量测试技术研究所 | 基于双管腔的脉动压力测量修正方法 |
CN109341945B (zh) * | 2018-11-23 | 2020-08-21 | 中国航空工业集团公司北京长城计量测试技术研究所 | 基于双管腔的脉动压力测量修正方法 |
CN111504595A (zh) * | 2020-03-31 | 2020-08-07 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种高速列车模型在声学风洞内的表面声载荷测量方法 |
CN111504595B (zh) * | 2020-03-31 | 2022-03-08 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种高速列车模型在声学风洞内的表面声载荷测量方法 |
CN112100895A (zh) * | 2020-10-28 | 2020-12-18 | 北京机电工程研究所 | 一种针对螺栓连接组合结构的分层模型修正方法和系统 |
CN112100895B (zh) * | 2020-10-28 | 2023-09-26 | 北京机电工程研究所 | 一种针对螺栓连接组合结构的分层模型修正方法和系统 |
CN112446091A (zh) * | 2020-11-19 | 2021-03-05 | 中国航天空气动力技术研究院 | 一种基于人工神经网络的脉动压力预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105989205B (zh) | 2019-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Numerical simulation of 3-D wing flutter with fully coupled fluid–structural interaction | |
CN105989205A (zh) | 飞行器表面脉动压力的确定方法 | |
Kim et al. | Numerical analysis on added resistance of ships | |
Chen et al. | Fully coupled fluid–structural interactions using an efficient high resolution upwind scheme | |
CN114662425B (zh) | 一种水轮机启停工况流场仿真预测方法及系统 | |
Ye et al. | Benchmark computations for flows around a stationary cylinder with high Reynolds numbers by RANS-overset grid approach | |
CN103400035B (zh) | 一种高可信度快速预测飞行器滚转动导数的方法 | |
CN105955928A (zh) | 基于cfd预报船舶阻力的计算方法 | |
Zhang et al. | A flutter prediction method with low cost and low risk from test data | |
Winter et al. | CFD-based aeroelastic reduced-order modeling robust to structural parameter variations | |
CN111428424A (zh) | 一种基于螺旋度的亚格子涡粘模型的验证方法及装置 | |
Ishihara et al. | A numerical study of the aerodynamic characteristics of ice-accreted transmission lines | |
Eguchi et al. | A finite element LES for high-Re flow in a short-elbow pipe with undisturbed inlet velocity | |
Wang et al. | Extraction of aerodynamic damping and prediction of vortex-induced vibration of bridge deck using CFD simulation of forced vibration | |
Kang et al. | A numerical investigation of two-degree-of-freedom VIV of a circular cylinder using the modified turbulence model | |
Chwalowski et al. | Fun3d analyses in support of the second aeroelastic prediction workshop | |
Wang et al. | Numerical simulation of transonic limit cycle oscillations using high-order low-diffusion schemes | |
Hine et al. | Fundamental Physics Validation Using CREATE-AV KestrelTM Part I | |
CN102270249B (zh) | 识别零部件特征频率的方法 | |
CN105224726A (zh) | 结构网格动网格技术用于非结构网格流场求解器的方法 | |
Ashcroft et al. | On the Simulation of Unsteady Turbulence and Transition Effects in a Multistage Low Pressure Turbine: Part I—Verification and Validation | |
Palacios et al. | Shape sensitivity of free-surface interfaces using a level set methodology | |
Orang | Genuinely characteristic-based scheme for the incompressible turbulent flows | |
Tatar | Global nonlinear aerodynamic reduced-order modeling and parameter estimation by radial basis functions | |
Jimbo et al. | Studies on unsteady flow characteristics in a high pressure turbine cascade based on a high-order large eddy simulation turbulence model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |