航空学报 > 2014, Vol. 35 Issue (8): 2182-2189   doi: 10.7527/S1000-6893.2013.0520

考虑随机型不确定性的浸入式颤振求解方法

戴玉婷, 杨超   

  1. 北京航空航天大学 航空科学与工程学院, 北京 100191
  • 收稿日期:2013-11-18 修回日期:2014-01-14 出版日期:2014-08-25 发布日期:2014-01-22
  • 通讯作者: 杨超,Tel.:010-82317510,E-mail:yangchao@buaa.edu.cn E-mail:yangchao@buaa.edu.cn
  • 作者简介:戴玉婷女,博士,讲师。主要研究方向:飞行器设计,气动弹性。Tel:010-82317528,E-mail:yutingdai@buaa.edu.cn;杨超男,博士,教授,博士生导师。主要研究方向:飞行器设计,气动弹性。Tel:010-82317510,E-mail:yangchao@buaa.edu.cn
  • 基金资助:

    国家自然科学基金(11302011,11172025);高等学校博士学科点专项科研基金(20131102120051)

Intrusive Flutter Solutions with Stochastic Uncertainty

DAI Yuting, YANG Chao   

  1. School of Aeronautic Science and Engineering, Beihang University, Beijing 100191, China
  • Received:2013-11-18 Revised:2014-01-14 Online:2014-08-25 Published:2014-01-22
  • Supported by:

    National Natural Science Foundation of China (11302011, 11172025); Research Fund for the Doctoral Program of Higher Education of China (20131102120051)

摘要:

气动弹性模型中的参数不确定性一般具有一定的分布规律,为了定量分析随机型参数不确定性对颤振的影响特性,考虑气动弹性系统中广义刚度的随机型不确定性,基于浸入式随机多项式展开(PCE)方法,在传统的颤振求解方法——p-k法的基础上,提出了针对不确定性气动弹性系统稳定性分析的增广p-k法——PCEPK (Polynomial Chaos Expansion with p-k)法,并将该方法应用到某机翼的颤振分析中,研究了均匀分布下的广义刚度不确定性对颤振边界的影响,并同基于结构奇异值μ理论的鲁棒颤振分析的结果和计算效率进行了对比。最后,采用标准的蒙特卡罗模拟(Monte Carlo Simulation, MCS)方法验证了结果的正确性。研究结果表明,PCEPK法计算的颤振边界范围是该分布下的“确定”结果,不因随机样本数而改变,克服了随机方法依赖样本数的缺点。同时,与基于结构奇异值理论的鲁棒颤振分析方法相比,它能够考虑不确定性参数分布对颤振特性的影响,具有更广泛的适用范围。

关键词: 气动弹性, 颤振, p-k法, 不确定性, 随机多项式展开, 浸入式

Abstract:

There is generally a certain distribution for parameter uncertainty in an aeroelastic model. In order to quantify the influence of stochastic uncertain parameters on the flutter boundary, a flutter analysis with generalized stiffness stochastic uncertainty is conducted in this work, based on the intrusive polynomial chaos expansion (PCE) theory. According to the traditional flutter solution method, namely the p-k method, an augmented PCEPK (polynomial chaos expansion with p-k) method for an uncertain aeroelastic system is presented to analyze its aeroelastic stability range. This flutter uncertainty analysis method is applied to a wing model, considering structural uncertainties with uniform distribution. A comparison of accuracy and computational time indicates that the range of flutter velocity by the PCEPK method is consistent with that obtained by the robust μ method. Finally, the standard Monte Carlo simulation (MCS) are employed to validate the results of the PCEPK method. In addition, though the parameter uncertainty is stochastic, the flutter velocity range is deterministic, invariant with the sample number, which resists the dependence on the random samples by the usual stochastic method. Moreover, the influence of the different distribution types of parameter uncertainty on flutter boundary can be obtained by the PCEPK method, which is more applicable when compared with the robust flutter analysis with μ method.

Key words: aeroelasticity, flutter, p-k method, uncertainty, polynomial chaos expansion, intrusive

中图分类号: