展弦比對(duì)NACA0012彈性平直機(jī)翼顫振臨界特性影響研究

摘 要：展弦比是飛行器機(jī)翼設(shè)計(jì)的重要參數(shù)，研究其對(duì)機(jī)翼顫振臨界特性的影響具有重要的工程應(yīng)用價(jià)值。本文采用三維雷諾平均N-S方程和S-A湍流模型求解NACA0012彈性平直機(jī)翼的非定常氣動(dòng)力，結(jié)合由模態(tài)法建立的機(jī)翼結(jié)構(gòu)動(dòng)力學(xué)方程，發(fā)展了計(jì)算流體力學(xué)（CFD）/計(jì)算結(jié)構(gòu)力學(xué)（CSD）顫振時(shí)域數(shù)值模擬方法。利用該方法研究了不同展弦比條件下機(jī)翼的顫振臨界速度和頻率。結(jié)果表明，由于展弦比較大且材質(zhì)輕柔，機(jī)翼會(huì)出現(xiàn)一定幅度的垂向靜彎曲變形。展弦比的增大不會(huì)影響機(jī)翼模態(tài)固有振型，但是會(huì)導(dǎo)致機(jī)翼顫振臨界速度、顫振頻率和模態(tài)固有頻率的降低。

關(guān)鍵詞：展弦比； 顫振臨界特性； 彈性平直機(jī)翼； 模態(tài)； 數(shù)值模擬

中圖分類號(hào)：V215.34 文獻(xiàn)標(biāo)識(shí)碼：A DOI：10.19452/j.issn1007-5453.2024.06.006

基金項(xiàng)目： 航空科學(xué)基金（201958058001）；中國空氣動(dòng)力研究與發(fā)展中心旋翼空氣動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室研究開放課題（RAL202204-9）

機(jī)翼展弦比的定義是展長b與平均弦長cmean的比值，也可以表示為展長的平方b2與參考面積S的比值[1]。為了降低誘導(dǎo)阻力，提高升阻比和任務(wù)載荷，高空長航時(shí)無人機(jī)、大型運(yùn)輸機(jī)和大型客機(jī)等飛行器在追求結(jié)構(gòu)輕量化的同時(shí)，普遍采用展弦比在10以上的大展弦比機(jī)翼[2]。然而，輕量化結(jié)構(gòu)和大展弦比設(shè)計(jì)使得機(jī)翼容易與非定常流動(dòng)進(jìn)行耦合，誘發(fā)出復(fù)雜的氣動(dòng)彈性問題，如顫振、彎曲和扭轉(zhuǎn)發(fā)散等[3]，不僅影響飛行器的運(yùn)行性能，還可能導(dǎo)致機(jī)翼產(chǎn)生疲勞損壞甚至結(jié)構(gòu)破壞，產(chǎn)生嚴(yán)重后果[4]。因此，大展弦比彈性機(jī)翼顫振分析對(duì)于提高飛行器的穩(wěn)定性和安全性具有重要意義。

顫振臨界特性（即顫振臨界速度和頻率）的預(yù)測是機(jī)翼顫振分析的重要部分。對(duì)于飛行器機(jī)翼而言，隨著來流速度U￥的增大，在一定U￥范圍內(nèi)，由氣流擾動(dòng)引起的振動(dòng)會(huì)衰減。當(dāng)U￥達(dá)到某一臨界值時(shí)，機(jī)翼振動(dòng)幅值剛好維持不變，這一速度在顫振分析中被稱為顫振臨界速度Uf，等幅振動(dòng)的頻率即為顫振頻率ωf。當(dāng)U￥超過Uf后，振動(dòng)會(huì)發(fā)散，機(jī)翼會(huì)發(fā)生顫振[5]。當(dāng)前，大展弦比彈性機(jī)翼顫振臨界特性的研究主要涉及兩個(gè)科學(xué)問題：（1）顫振分析方法的構(gòu)造和驗(yàn)證；（2）外掛載荷分布[6]、根部迎角[7]、操縱面及其偏轉(zhuǎn)角[8]、材料剛度[9]和密度[10]等對(duì)機(jī)翼顫振臨界特性的影響。速度-阻尼比外推法[11]、顫振裕度法[12]和p-k法[13]均為學(xué)術(shù)界和工程界中常用的顫振分析方法。然而，這些方法采用的氣動(dòng)力計(jì)算手段（如西奧道森模型、渦格法和偶極子格網(wǎng)法等）大多基于勢流理論，未考慮氣體的黏性效應(yīng)。近年來，隨著計(jì)算機(jī)性能的發(fā)展和大規(guī)模并行計(jì)算的普及，以N-S方程為基礎(chǔ)的計(jì)算流體力學(xué)（CFD）/計(jì)算結(jié)構(gòu)力學(xué)（CSD）時(shí)域數(shù)值模擬方法[14]已廣泛應(yīng)用于飛行器機(jī)翼的設(shè)計(jì)。CFD/CSD方法直接從流體和結(jié)構(gòu)基本方程出發(fā)，合理性和可靠性較高，而且除了顫振臨界速度和頻率以外，還能夠獲得振動(dòng)響應(yīng)、瞬時(shí)氣動(dòng)力和流場等信息。

綜上所述，現(xiàn)有的大展弦比彈性機(jī)翼顫振臨界特性研究大多在特定展弦比條件下開展，鮮有涉及展弦比對(duì)機(jī)翼顫振臨界速度和頻率的影響。本文以大展弦比NACA0012彈性平直機(jī)翼為研究對(duì)象，采用模態(tài)法進(jìn)行機(jī)翼的結(jié)構(gòu)動(dòng)力學(xué)分析，結(jié)合三維雷諾平均N-S方程和S-A湍流模型[15]，建立了CFD/CSD顫振時(shí)域數(shù)值模擬方法。利用該方法對(duì)機(jī)翼顫振臨界速度和頻率隨展弦比的變化規(guī)律進(jìn)行研究。

1 CFD/CSD顫振時(shí)域數(shù)值模擬方法

1.1 流動(dòng)控制方程

本文采用自研CFD軟件求解大展弦比NACA0012彈性平直機(jī)翼的非定常氣動(dòng)力，其控制方程為積分形式的守恒型三維非定常雷諾平均N-S方程，直角坐標(biāo)系下的形式為

采用基于預(yù)估-校正技術(shù)的四階雜交線性多步法[17]對(duì)式（4）進(jìn)行時(shí)域推進(jìn)求解。

圖1給出了大展弦比NACA0012彈性平直機(jī)翼顫振時(shí)域數(shù)值模擬流程圖。可以看出，機(jī)翼的結(jié)構(gòu)參數(shù)矩陣M、G和K和固有振型由有限元模態(tài)分析獲得，氣動(dòng)力矩陣F由CFD數(shù)值模擬確定。一般而言，翼面的CFD網(wǎng)格點(diǎn)和模態(tài)分析控制點(diǎn)是不重合的，因此在每一個(gè)時(shí)間步內(nèi)，需要通過結(jié)構(gòu)固有振型在翼面處進(jìn)行數(shù)據(jù)傳遞：壓力分布通過振型轉(zhuǎn)化為廣義氣動(dòng)力，進(jìn)而求解機(jī)翼結(jié)構(gòu)力學(xué)方程獲得廣義位移，廣義位移再通過振型轉(zhuǎn)化為翼面物理位移，作為下一個(gè)時(shí)間步CFD求解的邊界條件。模態(tài)分析得到的結(jié)構(gòu)固有振型通過徑向基函數(shù)插值法[18]轉(zhuǎn)化為翼面CFD網(wǎng)格點(diǎn)描述下的振型分布。機(jī)翼振動(dòng)后CFD網(wǎng)格的變形和運(yùn)動(dòng)也采用基于徑向基函數(shù)插值的動(dòng)網(wǎng)格技術(shù)[19]實(shí)現(xiàn)。

1.3 算例驗(yàn)證

AGARD445.6機(jī)翼顫振特性風(fēng)洞試驗(yàn)是由美國國家航空航天局（NASA）蘭利研究中心完成的。目前，該機(jī)翼模型已成為國際上顫振程序考核的標(biāo)準(zhǔn)三維算例。

該機(jī)翼模型的幾何參數(shù)為：展弦比為1.644，根梢比為1.517，1/4弦線后掠角為45°，翼型為NACA65A004。機(jī)翼試驗(yàn)?zāi)Ｐ偷牟牧蠟樘一ㄐ哪尽C(jī)翼顫振分析取前四階振動(dòng)模態(tài)，其固有振型依次為一階垂向彎曲、一階扭轉(zhuǎn)、二階垂向彎曲和二階扭轉(zhuǎn)。

2 結(jié)果與分析

本文的研究重點(diǎn)在于展弦比對(duì)NACA0012彈性平直機(jī)翼顫振臨界特性（即顫振臨界速度和頻率）的影響。機(jī)翼模型的幾何參數(shù)、材料參數(shù)和飛行環(huán)境見表1，考慮10、12、14和16的4個(gè)展弦比AR工況。顫振數(shù)值模擬取前五階振動(dòng)模態(tài)，圖3給出了不同展弦比條件下機(jī)翼模態(tài)固有頻率。從圖3中可以看出，隨著展弦比的增大，機(jī)翼模態(tài)固有振型不變，固有頻率逐漸降低。機(jī)翼的第一、二和五階振動(dòng)模態(tài)為垂向彎曲模態(tài)，第三階為扭轉(zhuǎn)模態(tài)，第四階則為弦向彎曲模態(tài)。

以展弦比AR=12工況為例，圖4給出了不同來流速度U￥時(shí)NACA0012彈性平直機(jī)翼的廣義位移時(shí)間響應(yīng)曲線。可以看出，當(dāng)U￥=30m/s和38.6m/s時(shí)，第一階振動(dòng)模態(tài)廣義位移（ξ1）存在一定幅度的靜漂移，說明由于展弦比較大，且材質(zhì)輕柔，機(jī)翼發(fā)生了垂向靜彎曲變形。當(dāng)U￥=30m/s時(shí)，各階模態(tài)廣義位移響應(yīng)表現(xiàn)為振動(dòng)衰減。當(dāng)U￥提高至38.6m/s時(shí)，廣義位移均表現(xiàn)為近似等幅響應(yīng)。進(jìn)一步提高U￥至50m/s時(shí)，廣義位移表現(xiàn)出明顯的振動(dòng)發(fā)散趨勢。綜合以上數(shù)值模擬結(jié)果可以得出，機(jī)翼在U￥=38.6m/s時(shí)處于臨界穩(wěn)定狀態(tài)，即AR=12時(shí)，機(jī)翼的顫振臨界速度Uf= 38.6m/s。

利用結(jié)構(gòu)固有振型將機(jī)翼的廣義位移轉(zhuǎn)化為翼面物理位移。圖5給出了AR=12的NACA0012彈性平直機(jī)翼在顫振速度（Uf=38.6m/s）時(shí)翼尖位置的垂向彎曲位移（以向上為正）和扭轉(zhuǎn)角響應(yīng)。圖6則給出了翼尖響應(yīng)的快速傅立葉變換（FFT）頻譜。可以看出，垂向彎曲位移的平衡位置不為0，而扭轉(zhuǎn)角的平衡位置為0°，表明機(jī)翼在垂向出現(xiàn)了靜彎曲變形。機(jī)翼的振動(dòng)幅值隨時(shí)間t基本保持不變，翼尖彎曲位移和扭轉(zhuǎn)角的主頻率均為31.04rad/s。因此，可以得出，AR=12時(shí)，機(jī)翼的顫振頻率ωf=31.04rad/s。

采用顫振數(shù)值模擬方法對(duì)其他AR條件下NACA0012彈性平直機(jī)翼的顫振臨界特性進(jìn)行計(jì)算，圖7給出了機(jī)翼顫振臨界速度和頻率隨展弦比的變化趨勢。可以看出，當(dāng)AR從10增大至16，機(jī)翼的顫振臨界速度會(huì)減小約45.71%，而顫振頻率會(huì)減小約33.16%。說明了AR的增大會(huì)引起機(jī)翼顫振臨界速度和頻率的降低。

3 結(jié)論

本文采用以三維雷諾平均N-S方程和S-A湍流模型為基礎(chǔ)的CFD數(shù)值模擬方法求解NACA0012彈性平直機(jī)翼的非定常氣動(dòng)力，結(jié)合模態(tài)法構(gòu)造的機(jī)翼結(jié)構(gòu)動(dòng)力學(xué)方程，發(fā)展了CFD/CSD顫振時(shí)域數(shù)值模擬方法，進(jìn)而研究了展弦比對(duì)機(jī)翼顫振臨界特性的影響，主要得到以下結(jié)論：

（1）由于展弦比較大，且材質(zhì)輕柔，機(jī)翼在垂向上會(huì)發(fā)生一定幅度的靜彎曲變形。

（2）機(jī)翼展弦比的增大不會(huì)帶來結(jié)構(gòu)模態(tài)固有振型的變化，但是會(huì)降低顫振臨界速度、顫振頻率和結(jié)構(gòu)模態(tài)固有頻率。

參考文獻(xiàn)

[1]吳子牛， 白晨媛， 李娟， 等. 空氣動(dòng)力學(xué)[M]. 北京： 北京航空航天大學(xué)出版社， 2016. Wu Ziniu， Bai Chenyuan， Li Juan， et al. Aerodynamics[M]. Beijing： Beihang University Press， 2016.（in Chinese）

[2]劉國春. 大展弦比機(jī)翼氣動(dòng)外形設(shè)計(jì)方案研究[J]. 飛機(jī)設(shè)計(jì)， 2011， 31（3）： 9-12+36. Liu Guochun. Research on design of high aspect wing shape[J]. Aircraft Design， 2011， 31（3）： 9-12+36. （in Chinese）

[3]Afonso F， Vale J， Oliveira é， et al. A review on non-linear aero‐elasticity of high aspect-ratio wings[J].Progress in Aerospace Sciences， 2017， 89： 40-57.

[4]Arena A， Lacarbonara W， Marzocca P. Nonlinear aeroelastic formulation and postflutter analysis of flexible high-aspectratio wings[J]. Journal of Aircraft， 2013， 50（6）： 1748-1764.

[5]楊超. 飛行器氣動(dòng)彈性原理[M]. 2版. 北京： 北京航空航天大學(xué)出版社， 2016. Yang Chao. Aeroelastic principle of aircraft[M]. 2ed edition. Beijing： Beihang University Press， 2016.（in Chinese）

[6]祁武超， 張賀銘， 田素梅， 等. 外掛物對(duì)大展弦比直機(jī)翼顫振特性的影響[J]. 計(jì)算力學(xué)學(xué)報(bào)， 2022， 39（5）： 598-607. Qi Wuchao， Zhang Heming， Tian Sumei， et al. Effects of external stores on flutter characteristics of a straight wing with high aspect ratio[J]. Chinese Journal of Computational Mechanics， 2022， 39（5）： 598-607. （in Chinese）

[7]付志超， 陳占軍， 劉子強(qiáng). 大展弦比機(jī)翼氣動(dòng)彈性的幾何非線性效應(yīng)[J]. 工程力學(xué)， 2017， 34（4）： 231-240. Fu Zhichao， Chen Zhanjun， Liu Ziqiang. Geometric nonlinear aeroelastic behavior of high aspect ratio wings[J]. Engineering Mechanics， 2017， 34（4）： 231-240. （in Chinese）

[8]胡志勇， 黨云卿， 陳海. 副翼偏轉(zhuǎn)角對(duì)大展弦比機(jī)翼顫振特性的影響研究[J]. 航空科學(xué)技術(shù)， 2016， 27（9）： 12-16. Hu Zhiyong， Dang Yunqing， Chen Hai. Influence of aileron de‐flection on the flutter of high-aspect-ratio wings[J]. Aeronauti‐cal Science & Technology， 2016， 27（9）： 12-16. （in Chinese）

[9]Qiao Shengjun， Jin Jiao， Ni Yangge， et al. Effect of stiffness on flutter of composite wings with high aspect ratio[J]. Mathemati‐cal Problems in Engineering， 2021， 2021： 1-14.

[10]余秋陽， 李興華， 李新濤， 等. 大展弦比NACA0012彈性平直機(jī)翼顫振特性研究[J]. 西北工業(yè)大學(xué)學(xué)報(bào)， 2022， 40（5）： 1071-1079. Yu Qiuyang， Li Xinghua， Li Xintao， et al. Study on flutter char‐acteristics of NACA0012 flexible straight wing with high-as‐pect-ratio[J]. Journal of Northwestern Polytechnical Universi‐ty， 2022， 40（5）： 1071-1079. （in Chinese）

[11]胡家亮， 吳江鵬， 脫朝智， 等. 3D打印飛機(jī)顫振風(fēng)洞試驗(yàn)?zāi)Ｐ驮O(shè)計(jì)與應(yīng)用[J]. 航空科學(xué)技術(shù)， 2022， 33（9）： 43-50. Hu Jialiang， Wu Jiangpeng， Tuo Chaozhi， et al. Design and application of 3D printed aircraft flutter wind tunnel test model[J]. Aeronautical Science & Technology， 2022， 33（9）： 43-50.（in Chinese）

[12]Zimmerman N H， Weissenburger J T. Prediction of flutter onset speed based on flight testing at subcritical speeds[J]. Journal of Aircraft， 1964， 1（4）： 190-202.

[13]Li Xing， Bai Yuguang， Xiao Ling， et al. Vibration control studies of a high-aspect-ratio wing with geometric nonlinearity[J]. Aerospace Science and Technology， 2022， 123： 107461.

[14]張科施， 凌圣博， 韓忠華. 跨聲速運(yùn)輸機(jī)機(jī)翼氣動(dòng)/結(jié)構(gòu)優(yōu)化平臺(tái)AeroStruct的發(fā)展及應(yīng)用[J]. 航空科學(xué)技術(shù)， 2022， 33（4）： 47-56. Zhang Keshi， Ling Shengbo， Han Zhonghua. Development and application of Aerostruct，an aerodynamic/structural optimiza‐tion platform for transonic transport aircraft wings[J]. Aero‐nautical Science & Technology， 2022， 33（4）： 47-56.（in Chinese）

[15]Spalart P R， Allmaras S R. A one-equation turbulence model for aerodynamic flows[J]. La Recherche Aérospatiale， 1994， 1： 5-21.

[16]葉正寅， 王剛， 張偉偉. 流固耦合力學(xué)基礎(chǔ)及其應(yīng)用[M]. 2版. 哈爾濱： 哈爾濱工業(yè)大學(xué)出版社， 2016. Ye Zhengyin， Wang Gang， Zhang Weiwei. Fundamentals of fluid-structure coupling and its application[M]. 2nd edition. Harbin： Harbin Institute of Technology Press， 2016. （in Chinese）

[17]Zhang Weiwei， Jiang Yuewen， Ye Zhengyin. Two better loosely coupled solution algorithms of CFD based aeroelastic simula‐tion[J]. Engineering Applications of Computational Fluid Me‐chanics， 2007， 1（4）： 253-262.

[18]Wright G B. Radial basis function interpolation： numerical and analytical developments[D]. Boulder： University of Colorado at Boulder， 2003.

[19]De Boer A， Van der Schoot M S， Bijl H. Mesh deformation based on radial basis function interpolation[J]. Computers & structures， 2007， 85（11-14）： 784-795.

[20]Yates E C J. AGARD standard aeroelastic configurations for dynamic response I：wing 445.6[R]. NASA-TM-100492， 1987.

[21]Silva W A， Lii B P， Chwalowski P. Evaluation of linear， inviscid， viscous， and reduced-order modelling aeroelastic solutions of the AGARD 445.6 wing using root locus analysis[J]. International Journal of Computational Fluid Dynamics， 2014， 28（3）： 122-139.

Effect of Aspect-ratio on Flutter Critical Characteristics of NACA0012 Flexible Straight Wing

Yu Qiuyang， Xu Shengjin

Tsinghua University， Beijing 100084， China

Abstract： Aspect-ratio is an important parameter in the design of aircraft wings， and studying its influence on the flutter critical characteristics of wings has significant engineering application value. This paper uses the threedimensional Reynolds Averaged N-S equations and the S-A turbulence model to solve the unsteady aerodynamic forces of high-aspect-ratio NACA0012 flexible straight wings. Combining with the wing dynamic equation established by the modal method， a CFD/CSD flutter time-domain numerical simulation method is developed. This method is employed to study the influence of aspect-ratio on the flutter critical characteristics of the wing. The results indicate that due to the high-aspect-ratio and soft material， the wing has a certain range of vertical static bending deformation. The increase in aspect-ratio does not affect the natural mode shape of the wing， but reduces the modal natural frequency， flutter critical velocity， and flutter frequency.

Key Words： aspect-ratio； flutter critical characteristics； flexible straight wing； mode； numerical simulation