?

平流層螺旋槳等離子體流動控制地面實驗方法

2015-06-21 15:08陳慶亞田希暉車學科聶萬勝周思引
實驗流體力學 2015年5期
關鍵詞:平流層雷諾數螺旋槳

陳慶亞, 田希暉, 車學科, 聶萬勝, 周思引

(裝備學院 航天裝備系, 北京 101416)

平流層螺旋槳等離子體流動控制地面實驗方法

陳慶亞, 田希暉*, 車學科, 聶萬勝, 周思引

(裝備學院 航天裝備系, 北京 101416)

根據螺旋槳雷諾相似準則和等離子體射流相似準則,提出了一種基于螺旋槳葉素理論,利用地面實驗設備開展平流層螺旋槳等離子體流動控制研究的實驗方法。首先根據螺旋槳幾何參數和運動參數計算葉素微段來流速度和迎角,然后根據螺旋槳雷諾相似準則確定常壓翼型風洞模擬平流層葉素流動的吹風參數,最后根據等離子體射流雷諾相似準則,確定激勵器和激勵電源參數模擬平流層等離子體射流并評估其流動控制效果。利用該方法研究了20km高度S1223翼型螺旋槳的等離子體流動控制效果,實驗表明:飛艇以5~20m/s的速度前進時,SDBD激勵電壓峰-峰值13.6kV,頻率10kHz時,誘導的等離子體射流使螺旋槳300r/min時推力最大可提高10.9%,600r/min時推力反而減小了0.52%~1.7%。

平流層螺旋槳;等離子體;流動控制;葉素理論;實驗方法

0 引 言

平流層飛艇具有駐留時間長、效費比高、生存能力強和高度適中等特點,作為搭載偵察、預警、通信等電子設備的平臺可以發揮其特有優勢,在民用和軍用領域均有很大現實意義[1]。但平流層飛艇螺旋槳在低雷諾數飛行環境下,葉素邊界層轉捩、分離等問題造成螺旋槳性能和效率降低[2],提高螺旋槳低雷諾數下飛行效率和寬工況飛行能力已成為重要研究課題。當前國內外對平流層螺旋槳增效的研究主要集中在優化翼型氣動外形、采用槳梢小翼等領域,獲得了有價值的研究成果,但是螺旋槳增效能力有限[3-7]。

表面介質阻擋放電(Surface Dielectric Barrier Discharge-SDBD)是等離子體流動控制技術的重要

研究方向,能夠有效抑制翼型流動分離,發揮減阻增升作用,提高螺旋槳低雷諾數下飛行效率和飛艇寬工況飛行能力;同時具有結構簡單、質量輕、功耗低和響應快等特點,已成為流動控制技術的熱點研究領域[8-9]。 目前,在翼型等離子體流動控制方面,主要通過優化激勵器幾何參數和激勵電源參數,探討翼型在不同的來流速度和迎角以及低雷諾數條件下抑制流動分離的控制效果[10-12],其中有美國圣母大學和我國南京航空航天大學等科研單位開發了等離子體流動控制的實驗無人機。

針對平流層高效螺旋槳,基于螺旋槳葉素理論,國防科技大學雷光新、劉巍等人[13]給出了平流層螺旋槳設計流程;中國電子科技集團二十七研究所姜波、翟彬等人[14]建立了計算NACA-0015翼型螺旋槳氣動性能的模型,并結合實驗數據驗證了模型的正確性。平流層螺旋槳等離子體流動控制技術方面,程鈺鋒等人[15-16]通過數值仿真開展了螺旋槳等離子體流動控制增效的研究,但是當前還沒有見到通過實驗方法開展平流層螺旋槳等離子體流動控制研究的公開文獻,難以推廣螺旋槳在平流層的應用。實驗開展高空平流層等離子體螺旋槳的研究需要低密度風洞,雖然可以實現,但是存在技術難度高、效率低以及成本高等問題。因此,迫切需要綜合采用現有地面實驗手段模擬平流層螺旋槳等離子體流動控制,定量評估效果。

作者根據劉沛清等[17-18]開展的地面常規風洞進行平流層螺旋槳實驗的可行性研究,在車學科等人[19]發展的地面條件下模擬高空等離子體流動控制效果的方法基礎上,基于螺旋槳葉素理論,采用實驗手段開展平流層螺旋槳等離子體流動控制增效的研究。

1 平流層螺旋槳地面實驗原理

1.1 平流層螺旋槳葉素運動參數確定

葉素理論又稱孤立槳葉理論[20],如圖1所示,沿切線方向將槳葉切成若干個稱之為葉素的微小段,認為繞過每個葉素的運動氣流是二維的,所以葉素之間、槳葉之間互不影響。

圖1 螺旋槳葉素示意圖

圖1中,r為葉素展向中點距槳心距離;R為螺旋槳半徑;ξ為葉素相對半徑,即ξ=r/R;br為葉素弦長;θr為葉素安裝角,dr為葉素微小段寬度。螺旋槳在平流層以轉速ns、相對來流前進速度V0運動時,根據運動相對性,其r處葉素運動分析如圖2所示。

圖2 葉素運動分析

在不考慮螺旋槳軸向干涉速度和環向干涉速度的情況下,忽略其干涉合成速度V1影響,幾何合成速度W0可以近似代替實際入流速度W1,則幾何入流角度Φ0可以近似等于實際入流角度Φ1。根據葉素運動分析,其周向速度U0、幾何合成速度W0、幾何入流角度Φ0和葉素迎角α可以表示為:

(1)

(2)

(3)

(4)

不同半徑r處葉素的幾何合成速度W0和迎角α分別由Wr、αr表示,根據式(1)~(4)可知,槳葉葉素的運動狀態可以由來流速度Wr、迎角αr2個參數確定。因此,在平流層以轉速ns、相對來流速度V0運動的螺旋槳可以轉化成不同半徑r處迎角為αr的葉素在速度為Wr的來流中的運動。

1.2 地面實驗參數確定

螺旋槳相同轉速在不同高度工作時,葉素繞流的雷諾數也會隨之變化。在平流層20km低氣壓、低密度條件下,螺旋槳雷諾數相當于低空飛行時的1/10,葉素近壁面層流難以克服逆壓梯度而容易發生流動分離造成螺旋槳性能下降。文獻[18]認為,密度基本不會對螺旋槳的相似參數產生影響。根據螺旋槳雷諾相似準則,地面二元翼型風洞實驗模擬平流層各葉素微段流動時對應雷諾數相等,即:

(5)

式中:ρh、uh為平流層大氣密度、動力粘度;ρg、ug為實驗氣體密度、動力粘度;b為實驗翼型弦長;Wrg為地面實驗風速。根據式(5)可得到地面實驗風速Wrg,實驗翼型迎角與平流層葉素迎角αr相同。

文獻[19]認為,根據等離子體射流雷諾相似準則,通過調整激勵器參數和激勵電源參數,可以由地面SDBD等離子體放電模擬不同高度下放電的流動控制效果。因此,實驗基于射流雷諾相似準則,根據高空平流層等離子體射流雷諾數確定地面激勵電源頻率和電壓。

1.3 螺旋槳推力評估

通過二元翼型風洞實驗系統,將采集的實驗數據經過編程計算得到分別垂直和平行于來流方向的升力dL系數CL、阻力dD系數CD,如圖3所示。

圖3 受力分析

根據翼型受力分析可知,葉素沿槳盤豎直向上的分量dT即螺旋槳產生的推力,則dL、dD和dT可表示為:

(6)

(7)

(8)

式中:pr為實驗測得來流空氣動壓。對槳葉葉素推力dT積分可得螺旋槳單槳葉推力T:

(9)

式中:r0為螺旋槳槳盤半徑。本實驗主要探討施加等離子體流動控制對螺旋槳推力的變化效果,采用計算得出的推力指標定量評估平流層螺旋槳性能的改變情況。

2 實驗方案

實驗模擬的平流層飛艇螺旋槳基于文獻[18]的設計,葉素剖面為S1223翼型,螺旋槳半徑R=3 250mm。實驗在裝備學院的低湍流度風洞中進行,湍流度0.02%,滿足平流層低湍流度要求,風洞試驗段尺寸為0.8m×1.0m×4.0m,實驗S1223翼型弦長200mm,展長790mm,采用環氧樹脂加工。實驗艙段和翼型如圖4所示,翼型上下表面分別布置18和17個測壓孔,測壓設備采用西北工業大學研制的壓力采集系統,由DS-32電子壓力掃描器、DS-104微壓校準單元和氣壓傳感器3部分組成,測量誤差0.1%。

分別在翼型上表面前緣、下表面后緣處布置1個SDBD激勵器,激勵器暴露電極和植入電極寬度分別為5mm、10mm,電極長度500mm,重合區長度230mm,電極間隙0mm,兩電極間通過粘貼3層Kapton膠帶作為介質阻擋層,其等效模擬的平流層20km螺旋槳激勵器的暴露電極、植入電極寬度為10mm和60mm,電極間隙0mm[19]。翼型上表面激勵器暴露電極靠近翼型前緣,距離前緣5.0mm;下表面激勵器暴露電極靠近翼型后緣,距離后緣5.0mm。根據上述激勵器鋪設方式可知上表面等離子體射流方向與來流方向相同,下表面等離子體射流與來流方向相反[20],翼型上表面1~5號測壓孔由于激勵器的遮擋無法正常測壓。

激勵電源采用HFHV30-1高頻高壓交流電源,輸出電壓±15kV,輸出頻率1~50kHz,上、下表面激勵器采用并聯方式,使得激勵器工作狀態時刻相同。電壓采用安捷倫N2771B高壓探頭測量,電流采用皮爾森電流線圈6595測量,電壓、電流測量結果使用安捷倫DSO3024A示波器進行顯示和記錄。

圖4 S1223翼型風洞實驗

實驗方案設定平流層20km處飛艇螺旋槳以300和600r/min 2種轉速工作,由于平流層平均風速為10m/s[21],同時又要滿足飛艇駐空懸停和機動能力,螺旋槳的前進速度V0分別取5、10和20m/s,根據前文所述實驗原理可得地面二元翼型風洞試驗參數,如表1所示。

文獻[19]通過低壓艙平板射流實驗,得出平流層20km、激勵電壓峰-峰值13.6kV時,等離子體沿壁面切向的射流雷諾數為36~40。在本實驗中,進一步考慮了等離子體對低氣壓稀薄氣體的壓縮作用導致

表1 風洞試驗參數

(b)螺旋槳轉速600r/min

密度的變化,沿壁面射流的雷諾數約為79~112,由于激勵電源頻率對雷諾數影響很小,實驗頻率采用10kHz,主要是為了更好模擬平流層激勵器放電誘導的流場,實驗電壓峰-峰值根據射流雷諾數分別為6和7kV。

3 實驗結果分析

3.1 300r/min時螺旋槳推力分析

300r/min(ns=5r/s)時不同相對半徑ξ處葉素小微段推力dT如圖5所示。

圖5中無控制指激勵器關閉狀態,在不同前進速度V0下,峰-峰值6和7kV 2種模式與無控制時的葉素推力進行比較可以看到,SDBD誘導的等離子體射流均有明顯增效作用。

相對半徑ξ介于0.30~0.45區間時,施加激勵與無控制的推力曲線重合度較高,這主要是此區間葉素迎角相對較大,出現了嚴重的流動分離,超出了等離子體抑制流動分離的能力范疇,因此該區間等離子體對葉素推力作用效果不明顯;ξ介于0.45~0.90區間段,該區間推力dT明顯提高,是等離子體增效的主要區間,具體在ξ=0.45處開始出現增效,其中ξ~(0.50,0.80)推力dT增加幅度最為顯著,ξ~(0.80,0.90)處的增加幅度相對減??;ξ介于0.90~0.975區間段,此范圍葉素已接近螺旋槳槳尖部位,迎角和弦長都比較小,幾乎沒有發生流動分離,等離子體不發揮增效作用,因此推力曲線重合。

隨著激勵電壓增大,等離子體射流對螺旋槳增效幅值的分布規律大致相同。7kV時增效作用略好,主要體現在ξ介于0.50~0.65區間段,說明該區間等離子體體積力和射流速度的增大可以更有效地抑制流動分離。

圖6給出了螺旋槳在實驗電壓6和7kV時等離子體射流對螺旋槳推力的增效情況??梢钥闯?,300r/min轉速時在3種前進速度V0下等離子體射流均可顯著增大螺旋槳推力,不同V0時其增效情況略有差異,其中7kV時的流動控制效果要好于6kV,螺旋槳推力分別增加了6.6%~10.8%、9.2%~10.9%。

3.2 600r/min時螺旋槳推力分析

螺旋槳轉速600r/min(ns=10r/s)時,在2種激勵電壓下的各葉素微段推力dT如圖7所示。

可以看出,600r/min時等離子體作用范圍主要介于ξ~(0.35~0.80)區間。其中,僅在V0=10、20m/s2種速度時,ξ~(0.45,0.55)區間的等離子體增加了葉素微小段推力,而其余工況和葉素區間均出現了使推力減小的負作用。其中,ξ介于0.65~0.8葉素區間螺旋槳推力減小幅度最大,此區間同時也是葉素推力下降最為劇烈的一段,說明該處流場處于臨界狀態,等離子體射流對流場的改變存在不確定性。峰-峰值6和7kV2種激勵電壓下槳葉推力分布曲線重合程度很高,表明螺旋槳轉速600r/min時,流場氣體動量很大,等離子體射流對流場的影響能力減小。

與螺旋槳300r/min時施加等離子體射流的效果比較,兩種轉速下等離子體對螺旋槳推力產生影響的葉素分布區間長度大致相同,相對半徑比區間長度均為0.45左右。但600r/min時等離子體在ξ=0.35

處開始產生氣動影響,相對于300r/min時ξ=0.45,其氣動作用區間整體向槳根移動了大約0.1的相對半徑長度。

圖8為螺旋槳600r/min時,施加等離子體射流后對螺旋槳推力的作用情況。結合葉素推力分布圖7可知:激勵器在6和7kV誘導等離子體射流時,均只在不同V0下的一小段葉素區間略微增加了葉素推力,但總作用效果減小了螺旋槳推力,減小幅值0.52%~2.1%。其中,6kV比7kV放電時使推力減小的幅度更大,對應分別減小了1.4%~2.1%、0.52%~1.7%。

圖7 600r/min時葉素推力dT

圖8 600r/min時槳葉推力T比較

綜上可知,螺旋槳轉速不同時等離子體射流產生的作用效果也會不同,由于不同轉速螺旋槳流場的不同,等離子體射流對推力存在增效和減效2種情況;在增效和減效2種情況時的等離子體射流,當提高激勵電壓,相比之前較低電壓均有利于提高螺旋槳性能;螺旋槳前進速度變化會使等離子體的控制作用發生輕微變化,但其影響程度與轉速相比差距較大,不是主要影響因素,說明飛艇在以較低的相對來流速度前進時,應主要針對螺旋槳轉速這一參量尋求最佳的等離子體增效方案。

4 結 論

(1) 提出了一種模擬平流層螺旋槳等離子體流動控制效果的地面實驗方法,該方法基于等離子體射流雷諾相似準則和螺旋槳雷諾相似準則,根據螺旋槳葉素理論,通過地面實驗分析平流層螺旋槳在低雷諾數條件下的氣動特性。

(2) 平流層螺旋槳SDBD激勵器采用微秒正弦連續式激勵模式時,等離子體射流對平流層螺旋槳推力的作用效果主要與轉速有關;在同一轉速下,螺旋槳不同前進速度對等離子體控制效果的影響很??;在螺旋槳300r/min可增效轉速,提高激勵器電壓有助于增大螺旋槳推力。

(3) 本文實驗條件下,模擬了平流層20km SDBD激勵器在峰-峰值13.6kV,頻率10kHz時的等離子體射流對螺旋槳的氣動效果,具體當螺旋槳以5,10和20m/s的速度前進時,在轉速300和600r/min 2種情況下,等離子體放電時使螺旋槳推力分別提高了6.6%~10.9%和減小了0.52%~1.7%。

[1] 郭建國, 周軍. 臨近空間低動態飛行器控制研究綜述[J]. 航空學報, 2013, 35(2): 320-331.

Guo Jianguo, Zhou Jun. Review of the control of low dynamic vehicles in near space[J]. Acta Aeronautica et Astronautica Sinica, 2013, 35(2): 320-331.

[2] Gambler D E, Arena A. Automated dynamic propeller testing at low Reynolds numbers[R]. AIAA-2010-853, 2010.

[3] Li Y H, Wu Y, Zhou M, et al. Control of the corner separation in a compressor cascade by steady and unsteady plasma aerodynamic actuation[J]. Experiment in Fluids, 2010, 48: 1015-1023.

[4] Ma R, Liu P Q. Optimization design study of low-Reynolds-number high-lift airfoils for the high-efficiency propeller of low-dynamic vehicles in stratosphere[J]. Science China Technological Sciences, 2010, 53(10): 2792-2807.

[5] Selig M S. Wind tunnel aerodynamic tests of six airfoils for use on small wind turbines[R]. AIAA-2004-1188, 2004.

[6] Benard N, Moreau E. Effects of altitude on the electromechanical characteristics of dielectric barrier discharge plasma actuators[R], AIAA-2010-4633, 2010.

[7] 許成杰, 楊旭東, 朱敏臨. 近空間槳梢小翼螺旋槳布局氣動增效研究[J]. 航空計算技術, 2011, 41(5): 61-64.

Xu Chengjie, Yang Xudong, Zhu Minlin. Aerodynamic synergistic mechanism of high altitude proplet propeller configurations[J]. Aeronautical Computing Technique, 2011, 41(5): 61-64.

[8] Wang J J, Choi K, Feng L H, et al. Recent development in DBD plasma flow control[J]. Progress in Aerospace Sciences, 2013, 62: 52-79.

[9] 王勛年, 王萬波, 黃勇, 等. 介質阻擋放電等離子體對翼型流動分離控制的實驗研究[J]. 實驗流體力學, 2011, 25(4): 9-14.

Wang Xunnian, Wang Wanbo, Huang Yong, et al. Investigation of flow separation control on an airfoil using DBD plasma actuators [J]. Journal of Experiments in Fluid Mechanics, 2011, 25(4): 9-14.

[10] Rizzetta D, Visbal M . Plasma flow control simulations of a low-Reynolds number low-aspect-ratio wing[J]. Computers and Fluids, 2012, 70: 95-114.

[11] 張攀峰, 王晉軍, 施威毅, 等. 等離子體激勵低速分離流動控制實驗研究[J]. 實驗流體力學, 2007, 21(2): 35-39.

Zhang Panfeng, Wang Jinjun, Shi Weiyi, et al. Experimental study on the separation control by plasma actuator in subsonic flow[J]. Journal of Experiments in Fluid Mechanics, 2007, 21(2): 35-39.

[12] Benard N, Jolibois J, Moreau E. Lift and drag performances of an axisymmetric airfoil controlled by plasma actuator[J]. Journal of Electrostatics, 2009, 67: 133-139.

[13] 雷光新, 劉巍, 楊濤. 平流層飛艇螺旋槳初步設計方法[J]. 導彈與航天運載技術, 2014, (2): 1-4.

Lei Guangxin, Liu Wei, Yang Tao. Primary design method of stratosphere airship propeller[J]. Missiles and Space Vehicles, 2014, (2): 1-4.

[14] 姜波, 翟彬, 靳磊. 基于葉素和動量組合理論的直升機性能計算[J]. 電光系統, 2011, (1): 45-48.

Jiang Bo, Zhai Bin, Jin Lei. Performance calculation of helicopter based on blade element and momentum theory[J]. Electronic and Electro-optical Systems, 2011, (1): 45-48.

[15] 程鈺鋒, 聶萬勝, 車學科. 等離子體提高螺旋槳槳根翼型氣動性能的仿真研究[J]. 核聚變與等離子體物理, 2012, 32(3): 265-270.

Cheng Yufeng, Nie Wansheng, Che Xueke. Numerical study on plasma flow control to enhance the aerodynamic characteristic of the aerofoil on propeller root region[J]. Nuclear Fusion and Plasma Physics, 2012, 32(3): 265-270.

[16] 程鈺鋒, 聶萬勝. 離子體提高平流層螺旋槳氣動性能的數值分析[J]. 核聚變與等離子體物理, 2012, 32(4): 372-378.

Cheng Yufeng, Nie Wansheng. Numerical analysis of the effect of plasma flow control on enhancing the aerodynamic characteristics of stratospheric screw propeller[J]. Nuclear Fusion and Plasma Physics, 2012, 32(4): 372-378.

[17] 劉沛清, 馬蓉, 段中喆, 等. 平流層飛艇螺旋槳地面風洞試驗[J]. 航空動力學報, 2011, 26(8): 1775-1781.

Liu Peiqing, Ma Rong, Duan Zhongzhe, et al. Ground wind tunnel test study of the propeller of stratospheric airships[J]. Journal of Aerospace Power, 2011, 26(8): 1775-1781.

[18] 劉沛清, 馬利川, 段中喆, 等. 平流層飛艇螺旋槳相似準則分析與驗證[J]. 北京航空航天大學學報, 2012, 38(7): 957-961.

Liu Peiqing, Ma Lichuan, Duan Zhongzhe, et al. Study and verification on similarity theory for propellers of stratospheric airships[J]. Journal of Beijing University of Aeronautics and Astronautic, 2012, 38(7): 957-961.

[19] 車學科, 聶萬勝, 侯志勇, 等. 地面實驗模擬高空等離子體流動控制效果研究[J]. 航空學報, 2015, 36(2): 441-448.

Che Xueke, Nie Wansheng, Hou Zhiyong, et al. Study on high altitude plasma flow control simulation through ground experiment[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(2): 441-448.

[20] 馮立好, 王晉軍, Choi K S, 等. 離子體環量控制翼型增升的實驗研究[J]. 力學學報, 2013, 45(6): 815-821.

Feng Lihao, Wang Jinjun, Choi K S, et al. Experimental investigation on lift of a plasma circulation control airfoil[J]. Chinese Journal of Theoretical and Applied Mechanics, 2013, 45(6): 815-821.

[21] 劉沛清. 空氣螺旋槳理論及其應用[M]. 北京: 北京航空航天大學出版社, 2006: 59-62.

[22] Liu P Q, Ma L C, Duan Z Z. Aerodynamics properties and design method of high efficiency-light propeller of stratospheric airships[R]. IEEE-2011-8041, 2011.

(編輯:張巧蕓)

Ground experimental method for stratospheric propeller plasma flow control

Chen Qingya, Tian Xihui*, Che Xueke, Nie Wansheng, Zhou Siyin

(Department of Space Equipment,Equipment Academy, Beijing 101416, China)

According to the propeller Reynolds-similarity theory and the plasma induced jet similarity theory, the method which is used to study the performance of plasma flow control of high-attitude propeller by ground experiment is put forward based on the blade element theory. First, the propeller geometry and motion parameters are converted to the inflow velocity and angle of attack on the blade element. Second, the induced plasma voltage and inflow parameters are determined based on the similarity theory in ground experiment. Finally, the collected experimental data is processed to assess the aerodynamic performance of propeller. A high-attitude 20km S1223 airfoil propeller plasma flow control experiment is conducted by this method. It is found that when the propeller revolves at the speed of 300r/min under forward velocity of 5~20m/s, the thrust of propeller is increased by 6.6%~10.9% with the high voltage AC input, but at the speed of 600r/min the thrust is reduced by 0.52%~1.7%.

high-attitude propeller;plasma;flow control;blade element theory;experimental method

1672-9897(2015)05-0090-07

10.11729/syltlx20140140

2014-12-03;

2015-03-11

國家自然科學基金(11205244),高分專項青年創新基金(GFZX04060103-5)

ChenQY,TianXH,CheXK,etal.Groundexperimentalmethodforstratosphericpropellerplasmaflowcontrol.JournalofExperimentsinFluidMechanics, 2015, 29(5): 90-96. 陳慶亞, 田希暉, 車學科, 等. 平流層螺旋槳等離子體流動控制地面實驗方法. 實驗流體力學, 2015, 29(5): 90-96.

V211.44

A

陳慶亞(1991-),男,河北滄州人,碩士研究生。研究方向:等離子體流動控制技術。通信地址:裝備學院航天裝備系(101416)。E-mail:chenqingya123@163.com

*通信作者 E-mail: tianxihui@sina.com

猜你喜歡
平流層雷諾數螺旋槳
青藏高原上空平流層水汽的時空演變特征
基于副氣囊的平流層浮空器高度控制
安第斯山上空平流層山地波統計特征的AIRS觀測研究
船用螺旋槳研究進展
基于CFD的螺旋槳拉力確定方法
非接觸機械密封端面間流體膜流動狀態臨界雷諾數的討論*
基于Transition SST模型的高雷諾數圓柱繞流數值研究
船模螺旋槳
亞臨界雷諾數圓柱繞流遠場氣動噪聲實驗研究
“平流層”是個啥——話說飛行高度
91香蕉高清国产线观看免费-97夜夜澡人人爽人人喊a-99久久久无码国产精品9-国产亚洲日韩欧美综合