CN107958102B - 一种用于高超声速气动热预测的偏差大气参数确定方法 - Google Patents
一种用于高超声速气动热预测的偏差大气参数确定方法 Download PDFInfo
- Publication number
- CN107958102B CN107958102B CN201711086208.XA CN201711086208A CN107958102B CN 107958102 B CN107958102 B CN 107958102B CN 201711086208 A CN201711086208 A CN 201711086208A CN 107958102 B CN107958102 B CN 107958102B
- Authority
- CN
- China
- Prior art keywords
- atmospheric
- deviation
- height
- density
- rho
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
本发明提供了一种用于高超声速气动热预测的偏差大气参数确定方法,属于高超声速飞行器气动热环境预示技术领域。该方法包括如下步骤:(1)、根据飞行弹道点的飞行高度H,由标准大气方程组,得出该弹道点对应的标准大气密度ρ;(2)、根据飞行弹道点的飞行高度H,由大气密度偏差Δρ与高度的关系,得出对应的大气密度偏差量Δρ,由标准大气密度ρ和大气密度偏差量Δρ,得出该飞行高度H对应的偏差大气密度ρ';(3)、根据偏差大气密度ρ',由标准大气方程组,反查出与偏差大气密度ρ'对应的偏差大气高度H';(4)、根据偏差大气高度H',由标准大气方程组,分别计算得到偏差大气压力P'和偏差大气温度T'。本发明相对其它方法来确定偏差大气参数,具有方便快速的特点。
Description
技术领域
本发明涉及一种用于高超声速气动热预测的偏差大气参数确定方法,属于高超声速飞行器气动热环境预示技术领域。
背景技术
对于高超声速飞行器气动热环境的预测来说,其输入条件是大气环境参数(压力、温度、密度)、飞行速度和攻角。临近空间高超声速滑翔飞行器在大气层内长时间飞行,实际大气环境参数与标准大气参数并不一致、存在一定的偏差,且偏差量随着高度的增加而增加,高度60km以上大气密度偏差超过了40%。大气环境参数的偏差会影响飞行器气动力的预测,进而影响飞行器的实际飞行弹道。大气环境参数和飞行弹道的偏差均会使飞行器气动热环境的预测产生一定的不确定度,进而影响飞行器的热防护设计。因此有必要在高超声速飞行器气动热环境的预测中考虑大气环境参数偏差的影响,以此来提高气动热环境预测准确度,为飞行器的热防护设计提供更为准确和合理的热环境输入。
表征大气环境的参数满足气体状态方程:
p=ρRT
其中p表示压力,ρ表示密度,T表示温度,R表示气体常数。由于目前工程实际中大气参数的偏差一般仅给出大气密度的偏差,大气压力和温度的偏差即使给出的话,由于大气参数偏差是针对多年大气参数进行统计分析后而给出,大气密度、压力和温度三个参数的偏差目前不满足气体状态方程,难以做到具有相关性,因此难以合理确定偏差大气的三个参数。由于高超声速气动热环境的输入条件包含三个参数,在仅有大气密度偏差或者密度、压力和温度偏差不具有相关性的情况下,如何预测考虑大气参数偏差情况下的高超声速飞行器气动热环境成为了一个需要解决的问题。
发明内容
本发明的技术解决问题是:针对目前仅有大气密度偏差或者密度、压力和温度偏差不具有相关性的情况,提供了一种用于高超声速气动热预测的偏差大气参数确定方法,可仅根据密度偏差数据通过反查标准大气函数得出符合气体状态方程的偏差大气参数,为偏差大气下的高超声速气动热预测提供合理的大气环境参数输入。
本发明的技术解决方案是:一种用于高超声速气动热预测的偏差大气参数确定方法,所述偏差大气参数包括偏差大气密度、偏差大气温度和偏差大气压力,该方法包括如下步骤:
(1)、根据飞行弹道点的飞行高度H,由标准大气方程组,得出该弹道点对应的标准大气密度ρ;
(2)、根据飞行弹道点的飞行高度H,由大气密度偏差Δρ与高度的关系,得出对应的大气密度偏差量Δρ,由标准大气密度ρ和大气密度偏差量Δρ,得出该飞行高度H对应的偏差大气密度ρ';
(3)、根据偏差大气密度ρ',由标准大气方程组,反查出与偏差大气密度ρ'对应的偏差大气高度H';
(4)、根据偏差大气高度H',由标准大气方程组,分别计算得到偏差大气压力P'和偏差大气温度T'。
所述标准大气方程组为:
其中,T为标准大气温度、P为标准大气压力,Hg为位势高度,rb为地球半径,R为气体常数和g为重力加速度;
Hb、Tb和pb分别是大气相应层的位势高度、温度和压力的下界值,b为温度梯度,它们随Hg的变化如下:
大气密度偏差Δρ与高度H的关系Δρ(H):
步骤(3)的具体实现为:
(3.1)、针对飞行弹道点的飞行高度H,设置搜索范围为:[H-10km,H+10km],初始化高度Hi为H-10km;
(3.2)、根据高度Hi,由标准大气压力和温度方程,分别计算得到高度Hi对应的大气压力Pi和大气温度Ti;
(3.3)、将大气压力Pi和大气温度Ti,代入气体状态方程,计算得到高度Hi对应的大气密度ρi;
(3.4)、判断大气密度ρi的大小,如果大气密度ρi与偏差大气密度ρ'相对差在预设的范围内,则确定高度Hi为偏差大气密度ρ'对应的偏差大气高度H',结束;否则,在高度Hi的基础上加上步长,更新高度Hi,重复执行步骤(3.1)~(3.4)。
所述预设的范围为:1%~2%。
所述步长范围为:0.01km~0.1km。
本发明与现有技术相比的有益效果是:
(1)、针对仅有大气密度偏差或密度、压力和温度偏差不具有相关性的情况,本发明提供了合理确定偏差大气参数的方法,可仅根据密度偏差数据通过反查标准大气函数得出符合气体状态方程的偏差大气参数,相对其它方法如复杂的统计分析方法或通过探空试验来确定偏差大气参数,具有方便快速的特点,并且节约成本;
(2)、本发明根据飞行高度得到偏差大气密度ρ',再由偏差大气密度ρ'反查偏差大气高度H',采用该偏差大气高度H'得到对应的偏差大气压力P'和偏差大气温度T',在保证了作为高超声速气动热预测的偏差大气参数满足气体状态方程的前提下,更加接近实际情况,提高后续高超声速气动热预测精度。
(3)、本发明采用遍历的手段反查偏差大气密度ρ'所对应的偏差大气高度H',计算准确性高。
(4)、所确定的偏差大气参数可以为高超声速气动热预测提供合理的参数输入,从而得出准确和合理的偏差大气下的飞行器气动热环境。
附图说明
图1为标准/偏差大气状态下表面热流沿x方向分布;
图2为本发明高超声速气动热预测中大气参数偏差的处理方法流程图。
具体实施方式
以下结合附图和具体实施例对本发明进行详细说明。
表征大气环境的参数满足气体状态方程:
p=ρRT
其中p表示压力,ρ表示密度,T表示温度,R表示气体常数。由于目前工程实际中大气参数的偏差一般仅给出大气密度的偏差,大气压力和温度的偏差即使给出的话,由于大气参数偏差是针对多年大气参数进行统计分析后而给出,大气密度、压力和温度三个参数的偏差目前不满足气体状态方程,难以做到具有相关性,因此难以合理确定偏差大气的三个参数。由于高超声速气动热环境的输入条件包含三个参数,在仅有大气密度偏差或者密度、压力和温度偏差不具有相关性的情况下,如何预测考虑大气参数偏差情况下的高超声速飞行器气动热环境成为了一个需要解决的问题。
为了解决该问题,进行了高超声速气动热环境与来流环境参数关系的理论分析和数值计算。首先对驻点热流与来流环境参数的关系进行分析。下式给出了高超声速飞行条件下的驻点热流理论公式:
qw0,FR=0.763Pr-0.6(h0-hw0)(ρ0μ0)0.4(ρw0μw0)0.1
其中Pr为普朗特数,h为焓值,μ为气体粘性系数,V为飞行速度,下标“0”表示驻点条件,下标“∞”表示来流条件,下标“w”表示壁面条件,hD为空气离解焓;hw0为驻点壁面处的焓;ρw0为驻点壁面处的密度,μw0为驻点壁面处的气体粘性系数;R0为端头半径。
由上式可以看出驻点热流公式非常复杂,它与大气参数(密度、压力和温度)并没有特别明显的直接关系。为了使驻点热流公式与大气参数具有较明晰的关系,对上式进行了理论推导,推导结果如下:
上式表明驻点热流只近似地与来流大气参数的密度相关。这样以来,上面提到的问题就可以得到解决:既然驻点热流只与大气密度相关,那么在大气参数有偏差的情况下,只需要知道大气密度的偏差就可以得到偏差状态下的驻点热流,而不需要知道大气压力和温度的偏差。
接下来对飞行器大面积上的热流公式进行了推导和分析,可以得出同样的结论—热流近似只与大气密度相关。利用高超声速气动热环境工程计算方法对标准大气和偏差大气状态下的驻点热流进行了计算,标准大气和偏差大气下的压力、温度和密度均满足气体状态方程,表1给出了标准大气状态和偏差大气状态下的大气参数和驻点热流关系,由表可以看出:只要密度偏差一定,无论压力和温度的偏差如何变化,热流偏差基本一致。利用高超声速气动热环境数值方法对标准大气和偏差大气状态下的表面热流进行了计算,标准大气和偏差大气下的压力、温度和密度均满足气体状态方程。图1给出了标准/偏差大气状态下的表面热流沿X向分布,图中:Standard为密度、压力和温度参数均为标准大气状态,rho_cons.P-30%T_-30%为密度无偏差,压力和温度取-30%偏差,rho_-30%P_-30%T_cons为温度无偏差,压力和密度取-30%偏差,rho_+42.9%P_consT_-30%为压力无偏差,密度取+42.9%偏差,温度取-30%偏差。从图中可以看出,在密度没有偏差的情况下,即使压力和温度有偏差,表面热流也基本没有偏差。上述两个计算印证了上面的理论推导和分析。
表1标准大气状态和偏差大气状态下的大气参数与驻点热流
根据上面的分析可知:对于偏差大气下的高超声速飞行器的气动热环境,在仅知大气密度偏差的情况下,无论压力和温度的偏差如何变化,只要三者(压力、温度和密度)满足气体状态方程,则气动热环境的大小与压力和温度的偏差基本无关。
基于上面的认识,本发明提出了一种用于高超声速气动热预测的偏差大气参数确定方法,所述偏差大气参数包括偏差大气密度、偏差大气温度和偏差大气压力,该方法包括如下步骤:
(1)、根据飞行弹道点的飞行高度H(单位为km),由标准大气方程组,得出该弹道点对应的标准大气密度ρ;
所述标准大气方程组为:
其中,T为标准大气温度、P为标准大气压力,Hg为位势高度;
rb为地球半径,R为气体常数和g为重力加速度,分别为:
rb=6356.766(km)
R=287.05281(J/kgK)
g=9.80665(m/s2)
Hb、Tb和pb分别是Hg相应大气层的位势高度、温度和压力的下界值,b为温度梯度,它们随Hg的变化如下:
(2)、根据飞行弹道点的飞行高度H,由大气密度偏差Δρ与高度的关系,得出对应的大气密度偏差量Δρ,由标准大气密度ρ和大气密度偏差量Δρ,得出该飞行高度H对应的偏差大气密度ρ':
ρ'=ρ+Δρ
大气密度偏差Δρ与高度的关系可以用函数表示,也可以用表格表示。表2列出了一种大气密度偏差Δρ与高度的关系Δρ(H)实施例:
大气密度偏差Δρ与高度的关系Δρ(H)
(3)、根据偏差大气密度ρ',由标准大气方程组,反查出与偏差大气密度ρ'对应的偏差大气高度H';
(3.1)、针对飞行弹道点的飞行高度H,设置搜索范围为:[H-10km,H+10km],初始化高度Hi为H-10km;
(3.2)、根据高度Hi,由标准大气压力方程,分别计算得到高度Hi对应的标准大气压力Pi和标准大气温度Ti;
(3.3)、将标准大气压力Pi和标准大气温度Ti,代入气体状态方程,计算得到高度Hi对应的标准大气密度ρi;
(3.4)、判断偏差大气密度ρi的大小,如果偏差大气密度ρi与偏差大气密度ρ'相对差在预设的范围(1%~2%)内,则确定高度Hi为与偏差大气密度ρ'对应的偏差大气高度H',结束;否则,在高度Hi的基础上加上步长(0.01km~0.1km),更新高度Hi,重复执行步骤(3.1)~(3.4)。
该方法可以在仅有密度偏差数据或密度、压力和温度偏差数据不具有相关性的情况下,利用密度偏差数据通过反查标准大气函数得出符合气体状态方程的偏差大气参数,从而可以为偏差大气下的高超声速气动热预测提供合理的大气环境参数输入。偏差大气参数与飞行速度V和飞行姿态角(攻角α、侧滑角β)作为输入条件,由高超声速气动热预测程序计算得出飞行器典型部位在偏差大气环境下的气动热环境。
上述实施例子只是对本发明的解释,而不能作为对本发明的限制,因此凡是与本发明思路类似的实施方式或大气参数偏差处理方法均在本发明的保护范围内。
本发明未进行详细描述部分属于本领域技术公知常识。
Claims (6)
1.一种用于高超声速气动热预测的偏差大气参数确定方法,所述偏差大气参数包括偏差大气密度、偏差大气温度和偏差大气压力,其特征在于包括如下步骤:
(1)、根据飞行弹道点的飞行高度H,由标准大气方程组,得出该弹道点对应的标准大气密度ρ;
(2)、根据飞行弹道点的飞行高度H,由大气密度偏差Δρ与高度的关系,得出对应的大气密度偏差量Δρ,由标准大气密度ρ和大气密度偏差量Δρ,得出该飞行高度H对应的偏差大气密度ρ';
(3)、根据偏差大气密度ρ',由标准大气方程组,反查出与偏差大气密度ρ'对应的偏差大气高度H';
步骤(3)的具体实现为:
(3.1)、针对飞行弹道点的飞行高度H,设置搜索范围为:[H-10km,H+10km],初始化高度Hi为H-10km;
(3.2)、根据高度Hi,由标准大气压力和温度方程,分别计算得到高度Hi对应的大气压力Pi和大气温度Ti;
(3.3)、将大气压力Pi和大气温度Ti,代入气体状态方程,计算得到高度Hi对应的大气密度ρi;
(3.4)、判断大气密度ρi的大小,如果大气密度ρi与偏差大气密度ρ'相对差在预设的范围内,则确定高度Hi为与偏差大气密度ρ'对应的偏差大气高度H',结束;否则,在高度Hi的基础上加上步长,更新高度Hi,重复执行步骤(3.1)~(3.4);
(4)、根据偏差大气高度H',由标准大气方程组,分别计算得到偏差大气压力P'和偏差大气温度T'。
2.根据权利要求1所述的一种用于高超声速气动热预测的偏差大气参数确定方法,其特征在于:所述标准大气方程组为:
其中,T为标准大气温度、P为标准大气压力,Hg为位势高度,rb为地球半径,R为气体常数和g为重力加速度;
Hb、Tb和pb分别是大气相应层的位势高度、温度和压力的下界值,b为温度梯度,它们随Hg的变化如下:
。
3.根据权利要求1所述的一种用于高超声速气动热预测的偏差大气参数确定方法,其特征在于大气密度偏差Δρ与高度H的关系Δρ(H):
。
5.根据权利要求1所述的一种用于高超声速气动热预测的偏差大气参数确定方法,其特征在于所述预设的范围为:1%~2%。
6.根据权利要求1所述的一种用于高超声速气动热预测的偏差大气参数确定方法,其特征在于所述步长范围为:0.01km~0.1km。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711086208.XA CN107958102B (zh) | 2017-11-07 | 2017-11-07 | 一种用于高超声速气动热预测的偏差大气参数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711086208.XA CN107958102B (zh) | 2017-11-07 | 2017-11-07 | 一种用于高超声速气动热预测的偏差大气参数确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107958102A CN107958102A (zh) | 2018-04-24 |
CN107958102B true CN107958102B (zh) | 2021-06-11 |
Family
ID=61964503
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711086208.XA Active CN107958102B (zh) | 2017-11-07 | 2017-11-07 | 一种用于高超声速气动热预测的偏差大气参数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107958102B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103708045A (zh) * | 2014-01-16 | 2014-04-09 | 中国人民解放军国防科学技术大学 | 一种探月飞船跳跃式再入的在线参数辨识方法 |
CN104568652A (zh) * | 2015-01-08 | 2015-04-29 | 浙江大学 | 一种高精度测量临近空间大气密度的方法及测量装置 |
CN105628051A (zh) * | 2014-10-29 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种嵌入式大气测量装置性能评估方法 |
CN106705996A (zh) * | 2016-11-25 | 2017-05-24 | 北京航天自动控制研究所 | 一种基于大气特征参数的飞行器导航信息修正方法 |
CN107132155A (zh) * | 2017-05-02 | 2017-09-05 | 中国西安卫星测控中心 | 基于实测大气密度的短时大气模型修正方法 |
-
2017
- 2017-11-07 CN CN201711086208.XA patent/CN107958102B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103708045A (zh) * | 2014-01-16 | 2014-04-09 | 中国人民解放军国防科学技术大学 | 一种探月飞船跳跃式再入的在线参数辨识方法 |
CN105628051A (zh) * | 2014-10-29 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种嵌入式大气测量装置性能评估方法 |
CN104568652A (zh) * | 2015-01-08 | 2015-04-29 | 浙江大学 | 一种高精度测量临近空间大气密度的方法及测量装置 |
CN106705996A (zh) * | 2016-11-25 | 2017-05-24 | 北京航天自动控制研究所 | 一种基于大气特征参数的飞行器导航信息修正方法 |
CN107132155A (zh) * | 2017-05-02 | 2017-09-05 | 中国西安卫星测控中心 | 基于实测大气密度的短时大气模型修正方法 |
Non-Patent Citations (5)
Title |
---|
"Atmospheric Environments for Entry, Descent and Landing (EDL)";C. G. Justus et al.;《https://ntrs.nasa.gov/search.jsp?R=20070032693 2020-08-04T08:55:34+00:00Z》;20070630;第1-37页 * |
"临近空间大气温度和密度特性分析";韩丁等;《遥感信息》;20170630;第32卷(第3期);第17-24页 * |
"北半球高空大气参数波动对临近空间飞行热环境的影响";陈闽慷等;《科学通报》;20170531;第62卷(第13期);第1402-1409页 * |
"基于特征模型的高超声速飞行器自适应控制研究进展";孟斌;《控制理论与应用》;20141231;第31卷(第12期);第1640-1649页 * |
"弹道导弹再入段大气参数修正";杨涛等;《弹箭与制导学报》;20090630;第29卷(第3期);第166-168页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107958102A (zh) | 2018-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108304601B (zh) | 一种高超声速飞行器边界层转捩的判断方法 | |
CN107103117B (zh) | 一种高超声速飞行器控制舵缝隙的热环境设计方法 | |
CN108287054B (zh) | 一种飞行条件下的转捩雷诺数获取方法 | |
CN107367941B (zh) | 高超声速飞行器攻角观测方法 | |
CN112668104B (zh) | 一种高超声速飞行器气动参数在线辨识方法 | |
CN109033493B (zh) | 基于无迹卡尔曼滤波的辨识高速旋转弹气动参数滤波方法 | |
CN105466660B (zh) | 一种采用加速度观测器的飞机气动参数分区分步辨识方法 | |
CN103471803A (zh) | 一种模型自由飞试验的气动参数确定方法 | |
CN109612676B (zh) | 基于飞行试验数据的气动参数反算方法 | |
RU2615220C1 (ru) | Способ определения управляющего сигнала по углу крена модели гиперзвукового летательного аппарата (ГЛА) для контроля аэродинамической идентичности по числам Рейнольдса траекторий полёта модели и натурного изделия при проведении опережающих лётных исследований | |
CN104077469A (zh) | 基于速度预测的分段迭代剩余时间估计方法 | |
CN110823216B (zh) | 一种用于惯导的自适应加速度计温度补偿方法 | |
CN105628051A (zh) | 一种嵌入式大气测量装置性能评估方法 | |
CN110826168B (zh) | 一种飞行器气动辨识修正方法及介质 | |
CN108303228B (zh) | 一种风洞试验的边界层流态控制装置及方法 | |
CN106096091B (zh) | 一种飞机运动模拟方法 | |
CN110187713A (zh) | 一种基于气动参数在线辨识的高超声速飞行器纵向控制方法 | |
CN107844128A (zh) | 一种基于复合比例导引的高超声速飞行器巡航段制导方法 | |
Rasuo | The influence of Reynolds and Mach numbers on two-dimensional wind-tunnel testing: An experience | |
CN107831653B (zh) | 一种抑制参数摄动的高超声速飞行器指令跟踪控制方法 | |
CN103034125A (zh) | 一种飞船返回舱气动与喷流控制力矩参数辨识方法 | |
CN107958102B (zh) | 一种用于高超声速气动热预测的偏差大气参数确定方法 | |
CN110765669B (zh) | 一种轴对称无翼无舵导弹主动段零升阻力系数辨识方法 | |
CN111412887B (zh) | 一种基于卡尔曼滤波的攻角、侧滑角辨识方法 | |
KR101944596B1 (ko) | 기상 데이터 보정 방법 및 장치 |
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 |