飞机带外挂构型的变刚度颤振特性计算方法研究

Study on variable stiffness flutter calculation method of aircraft with stores

摘要

对于飞机带外挂状态的颤振特性需要进行详细研究，包括低速和跨声速的颤振特性。传统基于面元法的颤振计算方法在分析飞机带外挂颤振特性时通常不考虑外挂的气动力作用，且无法考虑跨声速气动力的非线性影响。为了能够考虑带外挂飞机的颤振特性以及跨声速非线性的影响，建立了基于CFD（computational fluid dynamics）的变刚度颤振计算方法。开展了低速和跨声速不同马赫数状态的对称/反对称颤振特性研究，获得了颤振临界点，并通过低速和跨声速缩比模型颤振风洞试验进行验证。通过与低速颤振风洞试验结果对比，相比于面元法计算结果，基于CFD的颤振计算结果规律与低速风洞试验结果更接近。跨声速颤振风洞试验结果与基于CFD颤振计算结果基本一致。基于CFD的颤振计算方法能够较好模拟飞机带外挂颤振特性，可以用于工程计算。

Abstract

The flutter characteristics of aircraft with stores， including the flutter characteristics of low speed and transonic speed， need to be studied in detail. The flutter calculation method based on panel method has some limitations in analyzing flutter characteristics of aircraft with stores without considering the aerodynamic effects and the nonlinear effects of transonic aerodynamic. In order to consider the flutter characteristics of aircraft with stores and the influence of transonic nonlinearity， a variable stiffness flutter calculation method based on computational fluid dynamics（CFD）is generated. Based on the variable stiffness flutter calculation method， the symmetric and antisymmetric flutter characteristics of different Mach states at low speed and transonic speed are studied， and the flutter critical points are obtained and verified by the scaled model wind tunnel test. Compared with the experimental results of low-speed flutter wind tunnel， the results of the CFD-based flutter calculation is closer than those of the panel method. The results of transonic flutter wind tunnel test agree well with the results of CFD flutter calculation. The CFD-based flutter calculation method can effectively simulate the flutter characteristics of the aircraft with stores well and can be used for engineering calculation.

关键词

飞机带外挂状态 / 变刚度颤振计算 / 颤振风洞试验

Key words

aircraft with stores / variable stiffness flutter calculation / flutter wind tunnel test

参考文献

[1] 郭同庆. 复杂组合体跨音速非定常气动力和颤振计算[D]. 南京: 南京航空航天大学， 2006.

[2] 管德. 飞机气动弹性力学手册[M]. 北京: 航空工业出版社， 1994.

[3] BENDIKSEN O O. Review of unsteady transonic aerodynamics:theory and applications[J]. Progress in aerospace sciences, 2011， 47（2）: 135-167.

[4] 张伟伟, 高传强, 叶正寅. 复杂跨声速气动弹性现象及其机理分析[J]. 科学通报, 2018， 63（12）: 1095-1110.

[5] ZHANG Weiwei， GAO Chuanqiang， YE Zhengyin. The complexity and mechanism of transonic aeroelastic problems[J]. Chinese science bulletin, 2018， 63（12）: 1095-1110（in Chinese）.

[6] 许军, 马晓平. 飞机大展弦比外挂机翼弹性颤振特性分析[J]. 计算机仿真, 2015， 32（10）: 19-21.

[7] XU Jun， MA Xiaoping. Aircraft with high aspect ratio store wing elastic flutter characteristic analysis[J]. Computer simulation, 2015， 32（10）: 19-21（in Chinese）.

[8] 刘济科, 张宪民, 赵令诚, 等. 利用不对称外挂提高飞机颤振速度的方法研究[J]. 应用力学学报, 1997， 14（2）: 112-115.

[9] LIU Jike， ZHANG Xianmin， ZHAO Lingcheng， et al. On a method for raising flutter speed by using asymmetric stores[J]. Chinese journal of applied mechanics, 1997， 14（2）: 112-115（in Chinese）.

[10] 解江, 杨青青, 杨飞, 等. 翼吊发动机机翼跨声速颤振特性研究[J]. 飞行力学, 2018， 36（1）: 79-83.

[11] XIE Jiang， YANG Qingqing， YANG Fei， et al. Research on transonic flutter characteristics of a civil aircraft with wing-mounted engine[J]. Flight dynamics, 2018， 36（1）: 79-83（in Chinese）.

[12] 陆志良, 郭同庆, 管德. 跨音速颤振计算方法研究[J]. 航空学报, 2004， 25（3）: 214-217.

[13] LU Zhiliang， GUO Tongqing， GUAN De. A study of calculation method for transonic flutter[J]. Acta aeronautica et astronautica Sinica, 2004， 25（3）: 214-217（in Chinese）.

[14] JAMESON A， SCHMIDT W， TURKEL E. Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes[C]// 14th Fluid and Plasma Dynamics Conference. Reston， VA， USA:AIAA， 1981:AIAA 1981-1259.

[15] BYUN C， GURUSWAMY G P. Aeroelastic computations on wing-body-control configurations on parallel computers[C]// 37th Structure， Structural Dynamics and Materials Conference. Reston， VA， USA:AIAA， 1996:AIAA 1996-1389.

[16] GOODWIN S A， WEED R A， SANKAR L N， et al. Toward cost-effective aeroelastic analysis on advanced parallel computing systems[J]. Journal of aircraft, 1999， 36（4）: 710-715.

[17] YATES E C Jr. AGARD standard aeroelastic configurations for dynamic response. Candidate configuration I. -wing 445.6[EB/OL]. （ 2024-02-11）[ 2024-02-11]. https://ntrs.nasa.gov/citations/19880001820.

基金资助

中国航空工业集团公司沈阳飞机设计研究所博士创新基金资助项目(SBJ004)