力学学报 ›› 2016, Vol. 48 ›› Issue (2): 278-289.doi: 10.6052/0459-1879-14-406

• 流体力学 • 上一篇    下一篇

高超声速舵面热气动弹性不确定性及全局灵敏度分析

叶坤, 叶正寅, 屈展, 邬晓敬, 张伟伟   

  1. 西北工业大学航空学院, 西安 710072
  • 收稿日期:2014-12-17 修回日期:2016-01-04 出版日期:2016-03-18 发布日期:2016-03-03
  • 通讯作者: 叶坤,在读博士生,主要研究方向:高超声速热气动弹性.E-mail:yekun@mail.nwpu.edu.cn E-mail:yekun@mail.nwpu.edu.cn
  • 基金资助:

    国家自然科学基金重点资助项目(91216202).

UNCERTAINTY AND GLOBAL SENSITIVITY ANALYSIS OF HYPERSONIC CONTROL SURFACE AEROTHERMOELASTIC

Ye Kun, Ye Zhengyin, Qu Zhan, Wu Xiaojin, Zhang Weiwei   

  1. School of Aeronautics, Northwestern Polytechnical University, Xi'an 710072, China
  • Received:2014-12-17 Revised:2016-01-04 Online:2016-03-18 Published:2016-03-03

摘要:

鉴于高超声速中气动热预测的不确定性影响热气动弹性分析的可靠性,提出一种温度分布参数化模型,基于此模型,对高超声速舵面热气动弹性中气动热的不确定性及全局灵敏度进行分析,分析方法:求解N-S方程得到物面的温度分布,对此温度分布进行参数化,分别采用蒙特卡罗模拟(Monte Carlo simulation,MCS)方法和稀疏网格数值积分(spare grid numerical integration,SGNI)方法生成不确定性及全局灵敏度分析所需样本,对所有样本都进行热气动弹性分析,热气动弹性分析过程为:由样本得到温度分布,基于此温度分布,考虑热应力和材料属性的影响,对结构进行模态分析,将结构模态插值到气动网格,采用基于CFD的当地流活塞理论进行了气动弹性分析.分别在两种飞行状态下进行分析,计算结果表明:(1) M=5,H=15 km,结构固有频率和颤振分析结果的变异系数约为5.83%;(2) M=6,H=15 km,结构和颤振分析结果的变异系数约为8.84%.两种状态下,两个不确定参数的全局灵敏度都在50%左右,两者耦合作用很小,约为0.与MCS方法相比,SGNI方法显著的提高了不确定性分析效率.

关键词:

高超声速|热气动弹性|不确定性分析|全局灵敏度|稀疏网格方法|当地流活塞理论

Abstract:

Considering that the uncertainty of hypersonic aerothermodynamics prediction affects the reliability of aerothermoelastic analysis, a parameterized model for temperature distribution is therefore proposed.Based on this model, uncertainty and global sensitivity analysis on aerothermodynamics of hypersonic control surface aerothermoelastic are conducted.In the present analysis method, temperature distribution of the control surface is first obtained by solving NS equation and then parameterized.Using Monte Carlo simulation(MCS) method and spare grid numerical integration(SGNI) method to generate samples for analyzing uncertainly and global sensitivity and then analyzing all the samples, aerothermoelastic analysis is carried out as following:To get temperature distribution by the sample, then to analyze structural modal under the effect of structure thermal stress and material property, interpolate structural mode to the aerodynamic grid, and then to analyze aeroelasticity of the control surface in state space based on CFD local piston theory.Under two fly conditions, the calculation results show that:(1) With M=5 and H=15 km, the variation coefficient of natural frequency and flutter analysis is 5.83%,(2) With M=6 and H=15 km, variation coefficient of natural frequency of the structure and flutter analysis is 8.84%, and the global sensitivity of the two uncertainty parameters is about 50% under the two conditions.And the coupling of two parameters is about 0%, which is very small.Comparing with MCS method, SGNI method can be used to improve the efficiency of uncertainty analysis significantly.

Key words:

hypersonic|aerothermoelastic|uncertainty analysis|global sensitivity analysis|spare grid numerical integration|local flow piston theory

中图分类号: 

  • V211.47