CN105628333A - 一种高空高马赫数条件下的气动误差确定方法 - Google Patents
一种高空高马赫数条件下的气动误差确定方法 Download PDFInfo
- Publication number
- CN105628333A CN105628333A CN201410594551.5A CN201410594551A CN105628333A CN 105628333 A CN105628333 A CN 105628333A CN 201410594551 A CN201410594551 A CN 201410594551A CN 105628333 A CN105628333 A CN 105628333A
- Authority
- CN
- China
- Prior art keywords
- exp
- mach number
- aerodynamic
- deviation
- altitude
- 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
本发明属于飞行器结构技术领域,具体涉及一种高空高马赫数条件下的气动误差确定方法。本发明提出的高空高马赫数条件下的气动误差确定方法,从气动预示影响因素出发,解决了地面风洞试验能力不足问题,反映了高马赫数低雷诺数流动下的特殊物理现象,将飞行器理论外形气动预示结果和受实际加工制造水平及飞行环境的影响进行了联系,全面综合了气动误差影响源,弥补了传统方法的不足,所获结果具有可信性,可解决工程实际问题。本发明给出超出风洞能力条件下的高空高马赫数状态飞行器气动系数误差值,为飞行器控制***设计提供依据,提高成功飞行概率。
Description
技术领域
本发明属于飞行器结构技术领域,具体涉及一种高空高马赫数条件下的气动误差确定方法。
背景技术
为了对飞行器气动特性合理预示,需要给出相应的气动误差,明确预示偏差,指导气动特性预示结果的使用。现有公开文献对气动误差带确定方法介绍较少。由于传统飞行器剖面可被地面风洞试验能力所覆盖,因此传统气动误差带主要基于地面风洞试验给出,具体包括预示方法与试验结果的偏差和试验测量偏差两部分。而在高空高马赫数条件下(高度大于50km,马赫数大于15),现有国内试验能力已经不能覆盖飞行工况,仅靠风洞试验无法给出高空高马赫数条件下气动特性预示结果对应的气动误差。
发明内容
本发明的目的在于针对地面风洞试验能力的不足,提出可适用于高空高马赫数条件下的气动误差确定方法,以满足工程研制需求。
为达到上述目的,本发明所采取的技术方案为:
一种高空高马赫数条件下的气动误差确定方法,包括如下步骤:
(1)基于完全气体NS方程给出实际高空高马赫数飞行条件下飞行器六分量气动特性预示设计结果Ccp;
(2)选取典型飞行攻角、侧滑角和飞行器控制舵舵偏状态,在不少于2座的多座高超声速风洞中,针对马赫数8的风洞试验能力上限,开展相同状态风洞试验,且在单座风洞同状态重复性试验不少于3次,获取多座风洞试验数据;
(3)获得风洞试验数据均值Cexp_ave和上偏差Δexp_up、下偏差Δexp_down;
(4)基于气动特性预示方法,计算给出风洞状态下的气动数据Cexp_cp;
(5)计算获得飞行器设计工况条件下相同状态的气动特性计算校验结果Ccp_1~Ccp_n;以Ccp为基准,给出高空高马赫数条件下计算偏差Δcp_up和Δcp_down;
(6)将风洞条件下预示方法偏差Δexp_cp和高空条件下预示方法Δcp偏差取均方和,获得完全气体条件下预示方法偏差
(7)若飞行马赫数大于10以上,计算获得考虑真实气体效应的气动系数Crg,并与对应状态完全气体结果相减,获得真实气体效应对气动特性影响量Δrg;
(8)若飞行高度达于70km以上,计算获得考虑稀薄气体效应的气动系数Cld,并与对应状态完全气体结果相减,获得真实气体效应对气动特性影响量Δld;
(9)针对典型偏差外形,采用设计使用的完全气体NS方程,计算获得各类外形偏差条件下的气动系数Cpc1~Cpcn,并与对应状态理论外形气动系数结果相减,获得各类外形偏差对气动特性影响量Δpc1~Δpcn。将Δpc1~Δpcn取均方和,获得外形偏差对气动特性影响量
(10)将完全气体预示方法偏差Δpg、真实气体效应对气动特性影响量Δrg、稀薄气体效应对气动特性影响量Δld和外形偏差对气动特性影响量Δpc取均方和,获得气动系数C的误差值
所述步骤(3)具体为:将在多座风洞中开展的相同状态试验数据按t分布进行数学统计,获得风洞试验数据均值Cexp_ave和上偏差Δexp_up、下偏差Δexp_down。
所述步骤(4)具体为:若Cexp_cp与Cexp_ave量值接近,则以Cexp_cp为基准,根据Δexp_up、Δexp_down获得风洞条件下的气动预示方法误差上、下限Δexp_cp_up和Δexp_cp_down;否则应对气动特性预示方法进行确认和修正,直至Cexp_cp与Cexp_ave量值接近为止。
所述步骤(5)具体为:针对实际飞行剖面下的高度和马赫数,选取与风洞试验状态相同的攻角、侧滑角和飞行器控制舵舵偏状态,采用不同于设计使用的经过校验、有一定可信度的不少于7的多套完全气体NS方程气动数值模拟方法,计算获得飞行器设计工况条件下相同状态的气动特性计算校验结果Ccp_1~Ccp_n;以Ccp为基准,给出高空高马赫数条件下计算偏差Δcp_up和Δcp_down。
所述步骤(7)具体为:若飞行马赫数大于10以上,则针对实际飞行剖面下的高度和马赫数,选取与风洞试验状态相同的攻角、侧滑角和飞行器控制舵舵偏状态,采用平衡/非平衡气体模型,计算获得考虑真实气体效应的气动系数Crg,并与对应状态完全气体结果相减,获得真实气体效应对气动特性影响量Δrg。
所述步骤(8)具体为:若飞行高度达于70km以上,则针对实际飞行剖面下的高度和马赫数,选取与风洞试验状态相同的攻角、侧滑角和飞行器控制舵舵偏状态,采用完全气体NS方程滑移边界条件(70km~80km)/桥函数、DSMC方法(大于80km),计算获得考虑稀薄气体效应的气动系数Cld,并与对应状态完全气体结果相减,获得真实气体效应对气动特性影响量Δld。
所述步骤(9)具体为:根据飞行器实际加工制造偏差、烧蚀可能性和结构弹性变形情况,确定典型偏差外形,外形偏差的选取应具有独立性。
本发明所取得的有益效果为:
本发明提出的高空高马赫数条件下的气动误差确定方法,从气动预示影响因素出发,解决了地面风洞试验能力不足问题,反映了高马赫数低雷诺数流动下的特殊物理现象,将飞行器理论外形气动预示结果和受实际加工制造水平及飞行环境的影响进行了联系,全面综合了气动误差影响源,弥补了传统方法的不足,所获结果具有可信性,可解决工程实际问题。本发明给出超出风洞能力条件下的高空高马赫数状态飞行器气动系数误差值,为飞行器控制***设计提供依据,提高成功飞行概率。
具体实施方式
下面结合具体实施例对本发明进行详细说明。
本发明所述高空高马赫数条件下的气动误差确定方法如下:
(1)基于完全气体NS方程给出实际高空高马赫数飞行条件下飞行器六分量气动特性预示设计结果Ccp。
(2)选取典型飞行攻角、侧滑角和飞行器控制舵舵偏状态,在不少于2座的多座高超声速风洞中,针对马赫数8的风洞试验能力上限,开展相同状态风洞试验,且在单座风洞同状态重复性试验不少于3次,获取多座风洞试验数据。
(3)将在多座风洞中开展的相同状态试验数据按t分布进行数学统计,获得风洞试验数据均值Cexp_ave和上偏差Δexp_up、下偏差Δexp_down;
(4)基于气动特性预示方法,计算给出风洞状态下的气动数据Cexp_cp。若Cexp_cp与Cexp_ave量值接近,则以Cexp_cp为基准,根据Δexp_up、Δexp_down获得风洞条件下的气动预示方法误差上、下限Δexp_cp_up和Δexp_cp_down。否则应对气动特性预示方法进行确认和修正,直至Cexp_cp与Cexp_ave量值接近为止。
(5)针对实际飞行剖面下的高度和马赫数,选取与风洞试验状态相同的攻角、侧滑角和飞行器控制舵舵偏状态,采用其他(不同于设计使用的)经过校验、有一定可信度的多套(不少于7)完全气体NS方程气动数值模拟方法,计算获得飞行器设计工况条件下相同状态的气动特性计算校验结果Ccp_1~Ccp_n。以Ccp为基准,给出高空高马赫数条件下计算偏差Δcp_up和Δcp_down。
(6)将风洞条件下预示方法偏差Δexp_cp和高空条件下预示方法Δcp偏差取均方和,获得完全气体条件下预示方法偏差
(7)若飞行马赫数大于10以上,则针对实际飞行剖面下的高度和马赫数,选取与风洞试验状态相同的攻角、侧滑角和飞行器控制舵舵偏状态,采用平衡/非平衡气体模型,计算获得考虑真实气体效应的气动系数Crg,并与对应状态完全气体结果相减,获得真实气体效应对气动特性影响量Δrg。
(8)若飞行高度达于70km以上,则针对实际飞行剖面下的高度和马赫数,选取与风洞试验状态相同的攻角、侧滑角和飞行器控制舵舵偏状态,采用完全气体NS方程滑移边界条件(70km~80km)/桥函数、DSMC方法(大于80km),计算获得考虑稀薄气体效应的气动系数Cld,并与对应状态完全气体结果相减,获得真实气体效应对气动特性影响量Δld。
(9)根据飞行器实际加工制造偏差、烧蚀可能性和结构弹性变形情况,确定典型偏差外形。外形偏差的选取应具有独立性。针对典型偏差外形,采用设计使用的完全气体NS方程,计算获得各类外形偏差条件下的气动系数Cpc1~Cpcn,并与对应状态理论外形气动系数结果相减,获得各类外形偏差对气动特性影响量Δpc1~Δpcn。将Δpc1~Δpcn取均方和,获得外形偏差对气动特性影响量
(10)将完全气体预示方法偏差Δpg、真实气体效应对气动特性影响量Δrg、稀薄气体效应对气动特性影响量Δld和外形偏差对气动特性影响量Δpc取均方和,获得气动系数C的误差值
Claims (7)
1.一种高空高马赫数条件下的气动误差确定方法,其特征在于:包括如下步骤:
(1)基于完全气体NS方程给出实际高空高马赫数飞行条件下飞行器六分量气动特性预示设计结果Ccp;
(2)选取典型飞行攻角、侧滑角和飞行器控制舵舵偏状态,在不少于2座的多座高超声速风洞中,针对马赫数8的风洞试验能力上限,开展相同状态风洞试验,且在单座风洞同状态重复性试验不少于3次,获取多座风洞试验数据;
(3)获得风洞试验数据均值Cexp_ave和上偏差Δexp_up、下偏差Δexp_down;
(4)基于气动特性预示方法,计算给出风洞状态下的气动数据Cexp_cp;
(5)计算获得飞行器设计工况条件下相同状态的气动特性计算校验结果Ccp_1~Ccp_n;以Ccp为基准,给出高空高马赫数条件下计算偏差Δcp_up和Δcp_down;
(6)将风洞条件下预示方法偏差Δexp_cp和高空条件下预示方法Δcp偏差取均方和,获得完全气体条件下预示方法偏差
(7)若飞行马赫数大于10以上,计算获得考虑真实气体效应的气动系数Crg,并与对应状态完全气体结果相减,获得真实气体效应对气动特性影响量Δrg;
(8)若飞行高度达于70km以上,计算获得考虑稀薄气体效应的气动系数Cld,并与对应状态完全气体结果相减,获得真实气体效应对气动特性影响量Δld;
(9)针对典型偏差外形,采用设计使用的完全气体NS方程,计算获得各类外形偏差条件下的气动系数Cpc1~Cpcn,并与对应状态理论外形气动系数结果相减,获得各类外形偏差对气动特性影响量Δpc1~Δpcn;将Δpc1~Δpcn取均方和,获得外形偏差对气动特性影响量
(10)将完全气体预示方法偏差Δpg、真实气体效应对气动特性影响量Δrg、稀薄气体效应对气动特性影响量Δld和外形偏差对气动特性影响量Δpc取均方和,获得气动系数C的误差值
2.根据权利要求1所述的高空高马赫数条件下的气动误差确定方法,其特征在于:所述步骤(3)具体为:将在多座风洞中开展的相同状态试验数据按t分布进行数学统计,获得风洞试验数据均值Cexp_ave和上偏差Δexp_up、下偏差Δexp_down。
3.根据权利要求1所述的高空高马赫数条件下的气动误差确定方法,其特征在于:所述步骤(4)具体为:若Cexp_cp与Cexp_ave量值接近,则以Cexp_cp为基准,根据Δexp_up、Δexp_down获得风洞条件下的气动预示方法误差上、下限Δexp_cp_up和Δexp_cp_down;否则应对气动特性预示方法进行确认和修正,直至Cexp_cp与Cexp_ave量值接近为止。
4.根据权利要求1所述的高空高马赫数条件下的气动误差确定方法,其特征在于:所述步骤(5)具体为:针对实际飞行剖面下的高度和马赫数,选取与风洞试验状态相同的攻角、侧滑角和飞行器控制舵舵偏状态,采用不同于设计使用的经过校验、有一定可信度的不少于7的多套完全气体NS方程气动数值模拟方法,计算获得飞行器设计工况条件下相同状态的气动特性计算校验结果Ccp_1~Ccp_n;以Ccp为基准,给出高空高马赫数条件下计算偏差Δcp_up和Δcp_down。
5.根据权利要求1所述的高空高马赫数条件下的气动误差确定方法,其特征在于:所述步骤(7)具体为:若飞行马赫数大于10以上,则针对实际飞行剖面下的高度和马赫数,选取与风洞试验状态相同的攻角、侧滑角和飞行器控制舵舵偏状态,采用平衡/非平衡气体模型,计算获得考虑真实气体效应的气动系数Crg,并与对应状态完全气体结果相减,获得真实气体效应对气动特性影响量Δrg。
6.根据权利要求1所述的高空高马赫数条件下的气动误差确定方法,其特征在于:所述步骤(8)具体为:若飞行高度达于70km以上,则针对实际飞行剖面下的高度和马赫数,选取与风洞试验状态相同的攻角、侧滑角和飞行器控制舵舵偏状态,采用完全气体NS方程滑移边界条件(70km~80km)/桥函数、DSMC方法(大于80km),计算获得考虑稀薄气体效应的气动系数Cld,并与对应状态完全气体结果相减,获得真实气体效应对气动特性影响量Δld。
7.根据权利要求1所述的高空高马赫数条件下的气动误差确定方法,其特征在于:所述步骤(9)具体为:根据飞行器实际加工制造偏差、烧蚀可能性和结构弹性变形情况,确定典型偏差外形,外形偏差的选取应具有独立性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410594551.5A CN105628333B (zh) | 2014-10-29 | 2014-10-29 | 一种高空高马赫数条件下的气动误差确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410594551.5A CN105628333B (zh) | 2014-10-29 | 2014-10-29 | 一种高空高马赫数条件下的气动误差确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105628333A true CN105628333A (zh) | 2016-06-01 |
CN105628333B CN105628333B (zh) | 2018-06-22 |
Family
ID=56043472
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410594551.5A Active CN105628333B (zh) | 2014-10-29 | 2014-10-29 | 一种高空高马赫数条件下的气动误差确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105628333B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106932164A (zh) * | 2017-02-16 | 2017-07-07 | 北京临近空间飞行器***工程研究所 | 一种基于气动导数辨识结果的气动数据修正方法 |
CN110702364A (zh) * | 2019-10-22 | 2020-01-17 | 西北工业大学 | 针对桨尖马赫数影响的高空螺旋桨风洞试验数据修正方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3876755B2 (ja) * | 2002-04-18 | 2007-02-07 | 石川島播磨重工業株式会社 | 風洞装置 |
CN102521482A (zh) * | 2011-11-15 | 2012-06-27 | 中国航天空气动力技术研究院 | 一种粘性干扰效应气动力天地换算方法 |
CN103576554A (zh) * | 2013-11-07 | 2014-02-12 | 北京临近空间飞行器***工程研究所 | 基于控制需求的飞行器气动误差模型分量、分级设计方法 |
-
2014
- 2014-10-29 CN CN201410594551.5A patent/CN105628333B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3876755B2 (ja) * | 2002-04-18 | 2007-02-07 | 石川島播磨重工業株式会社 | 風洞装置 |
CN102521482A (zh) * | 2011-11-15 | 2012-06-27 | 中国航天空气动力技术研究院 | 一种粘性干扰效应气动力天地换算方法 |
CN103576554A (zh) * | 2013-11-07 | 2014-02-12 | 北京临近空间飞行器***工程研究所 | 基于控制需求的飞行器气动误差模型分量、分级设计方法 |
Non-Patent Citations (1)
Title |
---|
梁德旺等: "热完全气体的热力学特性及其N—S方程的求解", 《南京航空航天大学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106932164A (zh) * | 2017-02-16 | 2017-07-07 | 北京临近空间飞行器***工程研究所 | 一种基于气动导数辨识结果的气动数据修正方法 |
CN106932164B (zh) * | 2017-02-16 | 2019-02-19 | 北京临近空间飞行器***工程研究所 | 一种基于气动导数辨识结果的气动数据修正方法 |
CN110702364A (zh) * | 2019-10-22 | 2020-01-17 | 西北工业大学 | 针对桨尖马赫数影响的高空螺旋桨风洞试验数据修正方法 |
CN110702364B (zh) * | 2019-10-22 | 2020-06-30 | 西北工业大学 | 针对桨尖马赫数影响的高空螺旋桨风洞试验数据修正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105628333B (zh) | 2018-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104792633B (zh) | 一种飞机机体裂纹扩展寿命预测方法 | |
CN102722624B (zh) | 一种用于飞机适航审定及设计评估的飞行场景开发方法 | |
CN106840572A (zh) | 一种临近空间大展弦比柔性飞行器风洞试验数据修正方法 | |
CN103869383B (zh) | 基于卡尔曼滤波算法的大气数据计算机及其实现方法 | |
CN105115692A (zh) | 一种cfd数值模拟与风洞试验相结合的气动预测方法 | |
CN111159817A (zh) | 一种结冰风洞试验用的混合缩比机翼翼型设计方法 | |
Woodard et al. | Summary of ice shape geometric fidelity studies on an iced swept wing | |
CN108984909B (zh) | 一种基于Mar-Lin模型的含大损伤飞机复合材料结构剩余强度分析方法 | |
CN103914623A (zh) | 用于飞机结构腐蚀疲劳关键件延寿的寿命包线扩展方法 | |
Pérez et al. | Survey of aircraft structural dynamics non-linear problems and some recent solutions | |
CN109540459A (zh) | 一种气动特性数值计算结果修正方法 | |
CN105373647A (zh) | 一种通过地面滑跑试验辨识气动焦点的方法 | |
CN110702364B (zh) | 针对桨尖马赫数影响的高空螺旋桨风洞试验数据修正方法 | |
CN111274648B (zh) | 一种民用飞机前缘襟翼的分布式飞行载荷设计方法 | |
CN105628333A (zh) | 一种高空高马赫数条件下的气动误差确定方法 | |
CN103544402A (zh) | 一种疲劳开裂结构等效分析谱的构造方法 | |
Woodard et al. | Experimental Aerodynamic Simulation of a Scallop Ice Accretion on a Swept Wing | |
Sitaraman et al. | Coupled CFD/CSD prediction of rotor aerodynamic and structural dynamic loads for three critical flight conditions | |
CN104598745B (zh) | 一种月地高速再入返回飞行品质评估方法 | |
CN105334032A (zh) | 一种机型结构寿命的计算方法 | |
Biedron et al. | An examination of unsteady airloads on a UH-60A Rotor: Computation versus measurement | |
CN104677634A (zh) | 航空发动机机上侧风地面试验方法 | |
Ball et al. | F-35 full scale durability modeling and test | |
CN103983413A (zh) | 一种飞机起落架缓冲性能验证的落震试验方法 | |
Gebbink et al. | High-speed wind tunnel test of the CAE-AVM in DNW-HST for CFD validation purposes |
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 |