CN107977491A - 一种非稳态情况下飞行器空气舵缝隙的气动热评估方法 - Google Patents
一种非稳态情况下飞行器空气舵缝隙的气动热评估方法 Download PDFInfo
- Publication number
- CN107977491A CN107977491A CN201711117139.4A CN201711117139A CN107977491A CN 107977491 A CN107977491 A CN 107977491A CN 201711117139 A CN201711117139 A CN 201711117139A CN 107977491 A CN107977491 A CN 107977491A
- Authority
- CN
- China
- Prior art keywords
- airvane
- gap
- aircraft
- gap area
- rudder
- 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
Classifications
-
- 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
- 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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
一种非稳态情况下飞行器空气舵缝隙的气动热评估方法,包括步骤如下:一、通过数值求解飞行器流场的N-S方程,获得飞行器外壁表面热流;二、对舵缝隙内是否存在非定常效应进行判断并相应处理;三、获得若干周期内舵缝隙区域特征点处定常方法的平均热流;四、获得若干周期内特征点处非定常方法的平均热流并和定常方法结果比较,根据情况相应处理;五、获得舵缝隙区域干扰因子,利用曲线拟合方法获得干扰因子随舵偏变化的分段解析函数曲线;六、将干扰因子的分段函数曲线嵌入到气动热工程计算程序,获得飞行器在设定弹道时间段的舵缝隙区域热环境结果。本发明在保证空气舵缝隙气动热评估结果可靠性的同时能够有效减小评估结果的冗余度。
Description
技术领域
本发明涉及一种气动热评估方法。
背景技术
良好的控制特性和机动性对飞行器性能十分重要,而空气舵是高超声速飞行器进行姿态控制和维持稳定的重要部件,使得高超声速飞行器离不开空气舵的设计以及在飞行过程中空气舵的频繁摆动。由于结构特点及功能要求,使得空气舵和飞行器舱体之间不可避免地会存在缝隙,如图1所示,空气舵缝隙区域一般包括空气舵缝隙中与来流空气接触部分和空气舵附近受空气舵缝隙流动影响显著的区域。空气舵缝隙往往是飞行器流动规律最复杂、气动加热最严酷的区域之一。
在空气舵缝隙流动复杂性的众多表现中,非稳态是其特性之一。一方面,空气舵缝隙内往往会存在不稳定的分离涡流动结构,这种流动结构可能形成非定常流动,而使得空气舵缝隙内流动处于不稳定状态;另一方面,出于飞行控制或姿态维持的需要,空气舵需要频繁摆动而使得空气舵处于一种非稳态状态中,这两种情况下的空气舵缝隙气动热评估统称为非稳态情况下的空气舵缝隙气动热评估。
通常情况下,高超声速飞行器空气舵区域的气动热环境评估基本方式是采用气动热数值计算和地面试验相结合的方法进行,其中尤以数值计算能够评估真实状态下的气动热环境而被更广泛地采用。从现有技术来看,定舵偏情况下的定常气动热模拟已相当普遍并在工程上广泛采用,常用的手段是基于定常假设数值求解Navier-Stokes(N-S)方程而获得空气舵区域的热流,但对计算过程中可能出现的非定常流动现象则关注不足。若按传统方法,如果空气舵缝隙内出现分离涡流动结构,涡的蹿动导致空气舵舵缝隙内干扰区模拟的热流忽高忽低时,则根据数值模拟结果去取值则可能会取到极高值或极低值,得到的结果具有很大的不确定性,难以保证可靠;而如若采用非定常方法计算,则因其对计算资源的需求巨大而难以在工程上广泛使用。
对于空气舵频繁摆动的情况,因气动热评估时数值模拟工作量巨大而通常采用极限舵偏包络的方法进行弹道设计,极限舵偏情况下的热环境往往非常严酷,给防隔热方案设计带来极大困难,经常因为空气舵缝隙区域气动热评估结果过于严酷而使得方案难以闭合,也使得飞行器在研制过程中过多地受制于飞行器的姿控设计方案和空气舵缝隙区域的防隔热设计方案。
综上,对于高超声速飞行器空气舵缝隙区域的气动热评估,是影响飞行器姿控系统设计和局部防隔热设计方案的重点问题,也是当前设计水平下的普遍难点,如何准确、快速地评估空气舵缝隙区域的热环境,使之能够为飞行器方案的设计提供准确的输入,满足飞行试验的要求,十分迫切。
发明内容
本发明要解决的技术问题是:克服现有技术的不足,针对空气舵缝隙内产生的非定常流动现象和空气舵工作时频繁摆动这两类非稳态问题而导致空气舵缝隙区域气动热评估难以把握并且计算量巨大的问题,本发明提出了一种非稳态情况下空气舵缝隙干扰区的气动热环境设计方法,能够有效地考虑空气舵缝隙内存在非定常流动时的气动热变化特征,建立空气舵频繁摆动情况下的空气舵缝隙气动热评估方法,在可接受的计算资源耗费情况下,保证评估可靠性的同时有效减小评估的冗余度。
本发明所采用的技术方案是:一种非稳态情况下飞行器空气舵缝隙的气动热评估方法,包括如下步骤:
步骤一、利用有限体积法或有限差分法数值求解飞行器流场的N-S方程,获得飞行器外壁表面热流,并获取残差收敛曲线;
步骤二、对飞行器空气舵与舱体之间的缝隙内是否存在非定常效应进行判断,具体判断方法如下:
根据残差收敛曲线进行流场收敛情况判断,如果残差曲线下降4~5个量级,则表明飞行器表面的气动热计算结果收敛,则进入步骤五;
如果残差曲线下降小于3个量级,则结合热流计算结果进行判断,判断方法包括如下两种:
方法一:按照设定的步数输出空气舵缝隙区域的热流分布云图,观察热流分布云图是否随计算步数的变化呈周期性的变化,如果呈周期性变化则认为空气舵缝隙区域流动存在非定常效应,则进入步骤三;如果不呈周期性变化则认为空气舵缝隙区域流动不存在非定常效应,则进入步骤五;
方法二:通过输出空气舵缝隙区域固定位置处点的随计算步数变化的热流曲线,观察热流曲线是否随计算步数呈周期性变化,若呈周期性变化,则认为空气舵缝隙区域流动存在非定常效应,则进入步骤三;如果不呈周期性变化则认为空气舵缝隙区域流动不存在非定常效应,则进入步骤五;
步骤三、在空气舵缝隙内选取特征点,对若干周期内定常计算方法获得的热流数据进行平均处理得到空气舵缝隙区域在特征点处的平均热流结果;
步骤四、采用非定常计算方法获得步骤三中选取的特征点处若干周期内的平均热流结果,对比特征点处非定常计算方法和定常计算方法的热流结果:如果非定常计算方法的平均热流结果小于定常计算方法的平均热流结果或非定常计算方法的平均热流结果与定常计算方法的平均热流结果的偏差在10%以内,则进入步骤五;
步骤五、获得空气舵缝隙区域的干扰因子,对舵偏变化范围进行分段,利用曲线拟合方法获得干扰因子关于舵偏变化的分段解析函数曲线;干扰因子指步骤三中获得的空气舵缝隙区域的平均热流结果与非干扰热流结果的比值;
步骤六、将步骤五中获得的干扰因子关于随舵偏变化的分段解析函数曲线嵌入到气动热工程计算程序中,计算获得飞行器在设定的弹道时间段内的空气舵缝隙区域热环境曲线。
本发明与现有技术相比的优点在于:
(1)传统方法进行高超声速空气舵缝隙区域气动热评估时,往往忽略可能存在的非定常效应,若出现非定常效应时,则评估结果可能大幅偏大或大幅偏小,难以保证结果的可靠性;本发明采用理论分析和数值计算相结合的手段,能够有效判断高超声速飞行时空气舵缝隙区域内是否存在明显的非定常效应,并提出了一种合理解决存在非定常效应时空气舵缝隙的气动热评估方法,在保证评估结果可靠性的同时有效减小评估结果的冗余度;
(2)传统方法对于空气舵频繁摆动的弹道,由于工作量巨大而无法按真实舵偏进行气动热评估;而采用极限舵偏包络设计,则会对防隔热系统方案造成很大压力甚至影响飞行方案闭合;本发明基于提取不同舵偏情况下的空气舵缝隙区域数值结果,采用分段曲线拟合干扰因子的方法实现了空气舵频繁变化情况下的气动热评估,能够在可接受的较短时间内实现对飞行方案真实舵偏情况下的气动热评估,从而减小防隔热系统方案设计的压力。
附图说明
图1为空气舵安装结构及空气舵缝隙区域示意图;
图2为本发明的空气舵缝隙区域热环境设计流程图;
图3为残差收敛曲线示意;
图4为典型残差收敛曲线振荡情况;
图5(a)~图5(h)为典型周期内空气舵缝隙舵轴附近区域热流分布;
图6为典型空气舵缝隙位置点热流随计算步数的变化;
图7(a)~图7(l)为空气舵缝隙区域热流随时间的变化;
图8(a)~图8(d)为空气舵舵轴干扰因子随舵偏的变化及拟合曲线;
图9为典型弹道时间段空气舵舵轴气动热评估结果。
具体实施方式
下面结合附图对本发明的具体实施方式进行进一步的详细描述。
参照图2,示出了本发明实施例中一种控制舵缝隙的热环境设计方法的步骤流程图。一种非稳态情况下飞行器空气舵缝隙的气动热评估方法,包括步骤如下:
1.飞行器热环境数值计算
针对滑翔段典型飞行状态,通过求解N-S方程(式1)获得真实飞行器外壁表面热流。
式中,x、y、z为笛卡尔坐标系下坐标,为求解矢量,为无粘通量,为粘性通量,t为物理流动时间。采用有限体积法或有限差分法数值离散N-S方程,并采用一定的数值方法对空间项和时间项进行离散,然后迭代求解可获得飞行器外壁表面热流。湍流流态时,采用构建湍流模型(常用的有S-A湍流模型、SST湍流模型等)的方法进行求解。计算过程中保存流场的残差收敛情况,如图3所示,一般情况下残差随计算步数的增加向机器零趋近。
2.判断空气舵缝隙区域是否存在非定常效应
根据残差收敛曲线大致判断流场收敛情况,一般在网格合理的情况下,如果残差曲线能够下降4~5个量级,则基本表明气动热计算结果收敛,如图 3所示;如果残差曲线下降量级小于3个量级,并在较高的数值上振荡,则空气舵缝隙区域流场结果可能不收敛,如图4所示,这种情况进一步通过热流计算结果进行判断。按一定的步数输出空气舵缝隙区域的气动热计算结果,观察其结果是否随计算步数的变化呈周期性的变化,若呈周期性变化则认为空气舵缝隙区域流动存在非定常效应,如图5(a)~图5(h)所示,图中展示了计算步数86200步~87600步之间每隔200步空气舵缝隙舵轴附近区域热流分布及量值的变化,图5(a)中计算步数=86200,图5(b)中计算步数=86400,图5(c)中计算步数=86600,图5(d)中计算步数=86800,图5(e)中计算步数=87000,图5(f)中计算步数=87200,图5(g)中计算步数=87400,图5(h) 中计算步数=87600,可以看出该区域热流经历了一个多周期的变化,其中图 5(f)的热流分布和图5(a)相近,之后热流变化便继续延续之前的规律;或通过输出空气舵缝隙固定位置点的热流观察其是否随计算步数呈周期性变化,若是则认为空气舵缝隙区域流动存在非定常效应,如图6所示,图中展示了空气舵缝隙区域特征点处的热流变化,热流随计算步数的增加呈现出明显的周期性变化,最大值结果约为最小值结果的4倍。
3.定常方法的非定常效应计算
相对于空气舵缝隙尺寸的小尺度(mm量级),缝隙内的高速(km/s量级)流动可以说明上述时间周期非常短,相对于飞行器姿态变化时间来说,这个时间周期可以短到忽略不计。因而可以采用若干周期内的定常方法进行数据的平均处理得到空气舵区域的热环境结果,具体处理方法如式(2)。
式中,n为保存的热流文件个数,为统计平均热流。每个计算步数周期内建议保存5~9个热流文件,总周期数应不小于5个周期。
采用定常方法计算空气舵缝隙区域热流,并在若干周期内取多个样本的平均值,可以有效减小随机造成的计算结果取值大幅偏大或大幅偏小的概率,一定程度上保证气动热评估结果的可靠性。当然,该平均值还需通过非定常方法计算的结果验证才能更为充分地保证其数值的可靠性和有效性。
4.非定常效应计算结果的确定
以稳态计算结果作为初场,采用非定常方法计算了约2.5ms物理时间的流动结果,图7(a)~图7(l)给出了物理时间=0.02ms~0.46ms时间段内空气舵缝隙舵轴附近区域热流随时间的演化,图7(a)中物理时间=0.02ms,图7(b) 中物理时间=0.06ms,图7(c)中物理时间=0.10ms,图7(d)中物理时间=0.14ms,图7(e)中物理时间=0.18ms,图7(f)中物理时间=0.22ms,图7(g)中物理时间=0.26ms,图7(h)中物理时间=0.30ms,图7(i)中物理时间=0.34ms,图7(j)中物理时间=0.38ms,图7(k)中物理时间=0.42ms,图7(l)中物理时间=0.46ms,可以看到该区域热流随物理时间增加呈现近周期性变化,高热流区位置沿着舵轴上下蹿动。采用舵轴干扰区特征点的热流进行统计,发现舵轴特征点的热流也随周期性的变化,呈现明显的非定常流动的现象;对比特征点非定常方法和定常方法计算的热流结果,如表1,表明定常方法计算得到的统计平均结果接近非定常方法计算结果并且能够包络,表明了定常方法统计平均处理的有效性。
表1平均后的特征点热流(MW/m2)
5.获得舵偏摆动情况下的热流干扰因子拟合曲线
采用定常方法统计平均处理得到空气舵缝隙区域的热流结果,并和非干扰热流结果进行比较,得到热流干扰因子。以空气舵轴热流为例,表2给出了不同空气舵舵偏角情况下的空气舵舵轴干扰区热流变化情况和空气舵舵轴干扰区热流干扰因子(以驻点热流为非干扰参考热流)随空气舵舵偏的变化情况。
表2典型状态空气舵舵轴干扰区热流随空气舵舵偏变化
基于获得的空气舵缝隙区域干扰因子,按照其数值变化规律将空气舵舵偏分为四段,为(-10°~-1°)、(-1°~0°)、(0°~2°)、(2°~10°),并分别用解析函数曲线对干扰因子进行拟合,如图8(a)~图8(d)所示,拟合结果如下:
如图8(a)所示,拟合曲线1:y=0.003x2-0.225x-0.092;
如图8(b)所示,拟合曲线2:y=0.180x2+0.170x+0.150;
如图8(c)所示,拟合曲线3:y=0.252x2+1.036x+0.128;
如图8(d)所示,拟合曲线4:y=-0.046x2+0.823x+1.766。
6.沿弹道对舵缝隙区域进行热环境设计
根据获得的空气舵缝隙区域随舵偏变化的分段解析化的干扰因子,将其嵌入到气动热工程计算程序中,用于空气舵缝隙干扰区气动热环境的评估,典型弹道时间段空气舵缝隙干扰区的气动热评估结果如图9所示,可以看出,空气舵缝隙区域热流随舵偏变化敏感,热流峰值对应舵偏的正极值时刻;图中用直线a 标示出了极限舵偏,直线b标示出了极限舵偏情况下对应的热流,相比于极限舵偏下的气动热评估结果,基于真实舵偏情况下的气动热评估结果明显大幅下降。图9所得基于真实舵偏情况下的气动热评估结果,可用于后续飞行器防热、隔热方案的评估及设计。
上述实施例子只是对本发明的解释,而不能作为对本发明的限制,因此凡是与本发明思路类似的实施方式或空气舵缝隙区域的气动热评估方法均在本发明的保护范围内。
本发明未详细说明部分属于本领域技术人员公知技术。
Claims (4)
1.一种非稳态情况下飞行器空气舵缝隙的气动热评估方法,其特征在于,包括如下步骤:
步骤一、利用有限体积法或有限差分法数值求解飞行器流场的N-S方程,获得飞行器外壁表面热流,并获取残差收敛曲线;
步骤二、对飞行器空气舵与舱体之间的缝隙内是否存在非定常效应进行判断,如果空气舵缝隙区域流动存在非定常效应,则进入步骤三;如果空气舵缝隙区域流动不存在非定常效应,则进入步骤五;
步骤三、在空气舵缝隙内选取特征点,对若干周期内定常计算方法获得的热流数据进行平均处理得到空气舵缝隙区域在特征点处的平均热流结果;
步骤四、采用非定常计算方法获得步骤三中选取的特征点处若干周期内的平均热流结果,对比特征点处非定常计算方法和定常计算方法的热流结果:如果非定常计算方法的平均热流结果小于定常计算方法的平均热流结果或非定常计算方法的平均热流结果与定常计算方法的平均热流结果的偏差在10%以内,则进入步骤五;
步骤五、获得空气舵缝隙区域的干扰因子,对舵偏变化范围进行分段,利用曲线拟合方法获得干扰因子关于舵偏变化的分段解析函数曲线;
步骤六、将步骤五中获得的干扰因子关于随舵偏变化的分段解析函数曲线嵌入到气动热工程计算程序中,计算获得飞行器在设定的弹道时间段内的空气舵缝隙区域热环境曲线。
2.根据权利要求1所述的一种非稳态情况下飞行器空气舵缝隙的气动热评估方法,其特征在于:所述步骤二对飞行器空气舵与舱体之间的缝隙内是否存在非定常效应进行判断的具体判断方法如下:
根据残差收敛曲线进行流场收敛情况判断,如果残差曲线下降4~5个量级,则表明飞行器表面的气动热计算结果收敛,则空气舵缝隙区域流动不存在非定常效应;
如果残差曲线下降小于3个量级,则结合热流计算结果进行判断,判断方法如下:按照设定的步数输出空气舵缝隙区域的热流分布云图,观察热流分布云图是否随计算步数的变化呈周期性的变化,如果呈周期性变化则认为空气舵缝隙区域流动存在非定常效应;如果不呈周期性变化则认为空气舵缝隙区域流动不存在非定常效应。
3.根据权利要求1所述的一种非稳态情况下飞行器空气舵缝隙的气动热评估方法,其特征在于:所述步骤二对飞行器空气舵与舱体之间的缝隙内是否存在非定常效应进行判断的具体判断方法如下:
根据残差收敛曲线进行流场收敛情况判断,如果残差曲线下降4~5个量级,则表明飞行器表面的气动热计算结果收敛,则空气舵缝隙区域流动不存在非定常效应;
如果残差曲线下降小于3个量级,则结合热流计算结果进行判断,判断方法如下:通过输出空气舵缝隙区域固定位置处点的随计算步数变化的热流曲线,观察热流曲线是否随计算步数呈周期性变化,若呈周期性变化,则认为空气舵缝隙区域流动存在非定常效应;如果不呈周期性变化则认为空气舵缝隙区域流动不存在非定常效应。
4.根据权利要求2或3所述的一种非稳态情况下飞行器空气舵缝隙的气动热评估方法,其特征在于:所述干扰因子指步骤三中获得的空气舵缝隙区域的平均热流结果与非干扰热流结果的比值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711117139.4A CN107977491B (zh) | 2017-11-13 | 2017-11-13 | 一种非稳态情况下飞行器空气舵缝隙的气动热评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711117139.4A CN107977491B (zh) | 2017-11-13 | 2017-11-13 | 一种非稳态情况下飞行器空气舵缝隙的气动热评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107977491A true CN107977491A (zh) | 2018-05-01 |
CN107977491B CN107977491B (zh) | 2021-09-03 |
Family
ID=62013360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711117139.4A Active CN107977491B (zh) | 2017-11-13 | 2017-11-13 | 一种非稳态情况下飞行器空气舵缝隙的气动热评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107977491B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108995792A (zh) * | 2018-07-30 | 2018-12-14 | 上海机电工程研究所 | 复合材料结构的空气舵 |
CN109470374A (zh) * | 2018-11-13 | 2019-03-15 | 中国航天空气动力技术研究院 | 一种针对3-4mm缝隙内舵轴热流测量装置 |
CN110567669A (zh) * | 2019-08-06 | 2019-12-13 | 北京空天技术研究所 | 一种风洞试验的高速飞行器翼舵缝隙热流测量方法及装置 |
CN112287611A (zh) * | 2020-09-25 | 2021-01-29 | 北京临近空间飞行器系统工程研究所 | 一种降低凸起物气动热干扰的局部外形优化方法 |
CN112507630A (zh) * | 2020-11-16 | 2021-03-16 | 中国运载火箭技术研究院 | 热环境设计流场数据处理方法以及装置、存储介质 |
CN112798229A (zh) * | 2020-12-25 | 2021-05-14 | 中国空气动力研究与发展中心设备设计及测试技术研究所 | 一种气动光学效应自适应矫正系统的评估方法 |
CN113867381A (zh) * | 2021-12-02 | 2021-12-31 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种飞行器姿态控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120166148A1 (en) * | 2010-12-14 | 2012-06-28 | Japan Aerospace Exploration Agency | Method of designing natural laminar flow wing for reynolds numbers equivalent to actual supersonic aircraft |
CN104133933A (zh) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | 一种高超声速飞行器热环境下气动弹性力学特性分析方法 |
CN105808954A (zh) * | 2016-03-11 | 2016-07-27 | 中国航天空气动力技术研究院 | 一种适用于cfd数值模拟的周期非定常流场的预测方法 |
CN106202804A (zh) * | 2016-07-22 | 2016-12-07 | 北京临近空间飞行器系统工程研究所 | 基于数据库的复杂外形飞行器分布式热环境参数预测方法 |
CN107103117A (zh) * | 2017-03-27 | 2017-08-29 | 北京临近空间飞行器系统工程研究所 | 一种高超声速飞行器控制舵缝隙的热环境设计方法 |
CN107247839A (zh) * | 2017-06-08 | 2017-10-13 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种低速风洞虚拟飞行试验飞机模型设计方法 |
-
2017
- 2017-11-13 CN CN201711117139.4A patent/CN107977491B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120166148A1 (en) * | 2010-12-14 | 2012-06-28 | Japan Aerospace Exploration Agency | Method of designing natural laminar flow wing for reynolds numbers equivalent to actual supersonic aircraft |
CN104133933A (zh) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | 一种高超声速飞行器热环境下气动弹性力学特性分析方法 |
CN105808954A (zh) * | 2016-03-11 | 2016-07-27 | 中国航天空气动力技术研究院 | 一种适用于cfd数值模拟的周期非定常流场的预测方法 |
CN106202804A (zh) * | 2016-07-22 | 2016-12-07 | 北京临近空间飞行器系统工程研究所 | 基于数据库的复杂外形飞行器分布式热环境参数预测方法 |
CN107103117A (zh) * | 2017-03-27 | 2017-08-29 | 北京临近空间飞行器系统工程研究所 | 一种高超声速飞行器控制舵缝隙的热环境设计方法 |
CN107247839A (zh) * | 2017-06-08 | 2017-10-13 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种低速风洞虚拟飞行试验飞机模型设计方法 |
Non-Patent Citations (3)
Title |
---|
SHIN HYUNG RHEE ET AL.: ""A suggestion of gap flow control devices for the suppression of rudder cavitation"", 《JOURNAL OF MARINE ENCE & TECHNOLOGY》 * |
宁万春: ""某空空导弹燃气舵气动特性的非定常研究"", 《中国优秀硕士学位论文全文数据库工程科技辑》 * |
贾忠湖等: ""翼型舵面偏转非定常流动数值模拟"", 《航空计算技术》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108995792A (zh) * | 2018-07-30 | 2018-12-14 | 上海机电工程研究所 | 复合材料结构的空气舵 |
CN109470374A (zh) * | 2018-11-13 | 2019-03-15 | 中国航天空气动力技术研究院 | 一种针对3-4mm缝隙内舵轴热流测量装置 |
CN109470374B (zh) * | 2018-11-13 | 2020-11-10 | 中国航天空气动力技术研究院 | 一种针对3-4mm缝隙内舵轴热流测量装置 |
CN110567669A (zh) * | 2019-08-06 | 2019-12-13 | 北京空天技术研究所 | 一种风洞试验的高速飞行器翼舵缝隙热流测量方法及装置 |
CN112287611A (zh) * | 2020-09-25 | 2021-01-29 | 北京临近空间飞行器系统工程研究所 | 一种降低凸起物气动热干扰的局部外形优化方法 |
CN112287611B (zh) * | 2020-09-25 | 2023-05-12 | 北京临近空间飞行器系统工程研究所 | 一种降低凸起物气动热干扰的局部外形优化方法 |
CN112507630A (zh) * | 2020-11-16 | 2021-03-16 | 中国运载火箭技术研究院 | 热环境设计流场数据处理方法以及装置、存储介质 |
CN112507630B (zh) * | 2020-11-16 | 2022-12-13 | 中国运载火箭技术研究院 | 热环境设计流场数据处理方法以及装置、存储介质 |
CN112798229A (zh) * | 2020-12-25 | 2021-05-14 | 中国空气动力研究与发展中心设备设计及测试技术研究所 | 一种气动光学效应自适应矫正系统的评估方法 |
CN113867381A (zh) * | 2021-12-02 | 2021-12-31 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种飞行器姿态控制方法 |
CN113867381B (zh) * | 2021-12-02 | 2022-02-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种飞行器姿态控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107977491B (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107977491A (zh) | 一种非稳态情况下飞行器空气舵缝隙的气动热评估方法 | |
CN104714553B (zh) | 基于几何规划的滑翔飞行器末端能量管理轨迹规划方法 | |
Müller et al. | Modelling of crossflow-induced transition based on local variables | |
CN109446582B (zh) | 一种考虑地球自转的高精度降阶平稳滑翔动力学建模方法 | |
CN104699947B (zh) | 一种采用rans/les混合技术模拟飞行器摇滚运动的方法 | |
CN105069311B (zh) | 一种远程火箭发射初态误差传播估计方法 | |
Yongjie et al. | A new single-blade based hybrid CFD method for hovering and forward-flight rotor computation | |
CN109751204A (zh) | 一种风力机结冰数值模拟方法 | |
CN109871586A (zh) | 适用于共轴旋翼地面效应的气动干扰数值确定方法和系统 | |
CN110083968A (zh) | 基于修正气封泄露量影响数值模型的压气机特性预测方法 | |
CN106484980B (zh) | 一种固定舵二维弹道修正弹气动参数计算方法 | |
CN106708065A (zh) | 无人机降落控制方法、装置与无人机 | |
CN107357976A (zh) | 一种飞行器的动导数的计算方法 | |
Zhu et al. | Numerical study of a variable camber plunge airfoil under wind gust condition | |
CN102566446A (zh) | 基于线性模型组的无人直升机全包线数学模型构建方法 | |
CN114611416A (zh) | 导弹非线性非定常气动特性ls-svm建模方法 | |
Venkatachari et al. | Assessment of transition modeling capability in OVERFLOW with emphasis on swept-wing configurations | |
CN105224726B (zh) | 结构网格动网格技术用于非结构网格流场求解器的方法 | |
CN109885598A (zh) | 故障识别方法、装置、计算机可读存储介质及电子设备 | |
CN115542746A (zh) | 高超声速飞行器的能量管控再入制导方法及装置 | |
CN105955029B (zh) | 一种保鲁棒性的pid控制参数优化方法 | |
CN105628333B (zh) | 一种高空高马赫数条件下的气动误差确定方法 | |
Gaunaa et al. | Investigation of droplet path in a rain erosion tester | |
CN106991209A (zh) | 一种火星大气真实气体环境气动特性预测方法 | |
Chen et al. | Research on simulation method of missile adapter’s separation based on combined calculation |
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 |