• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于環(huán)量控制的虛擬舵面機翼氣動特性計算研究

    2020-01-21 09:36:33付志杰許和勇杜海王宇航徐悅
    航空科學技術(shù) 2020年5期
    關(guān)鍵詞:能耗

    付志杰 許和勇 杜海 王宇航 徐悅

    摘要:在機翼后緣應(yīng)用環(huán)量控制技術(shù)可以改變機翼的氣動力,為了研究環(huán)量控制技術(shù)在虛擬舵面飛行器上的控制效果,分別對不同舵偏角的機械舵面模型和不同射流動量系數(shù)的虛擬舵面模型進行了數(shù)值模擬。通過對比二者的升阻力系數(shù)和力矩系數(shù)曲線,發(fā)現(xiàn)前者在舵偏角θ= 0°、10°、20°和30°時的舵效分別與后者在Cμ= 0、0.005、0.009和0.012時相同,且θ與Cμ為二次多項式關(guān)系。為進一步研究環(huán)量控制系統(tǒng)在虛擬舵面上的氣動效率和能耗,對不同噴口高度的模型進行數(shù)值模擬,發(fā)現(xiàn)在射流速度相同的情況下,不同噴口高度的虛擬舵面的等效升阻比相同,但是大噴口的虛擬舵面需要耗費更大的功率。

    關(guān)鍵詞:虛擬舵面;環(huán)量控制;等舵效;能耗;等效升阻比

    中圖分類號:V11文獻標識碼:ADOI:10.19452/j.issn1007-5453.2020.05.002

    基金項目:國家自然科學基金(11972306);中央高?;究蒲袠I(yè)務(wù)費(310201901A004);翼型葉柵空氣動力學重點實驗室基金

    舵面是固定翼飛行器的重要組成部分,偏轉(zhuǎn)舵面可以改變自由來流對飛行器的作用力,從而實現(xiàn)對飛行器的控制,如升降舵控制俯仰運動、方向舵控制偏航運動和副翼控制滾轉(zhuǎn)運動。機械舵面已經(jīng)伴隨飛行器存在了百余年時間,且應(yīng)用在幾乎所有的固定翼飛行器上。存在時間長和應(yīng)用范圍廣似乎說明了機械舵面具有不可替代性。但近年來出現(xiàn)了一種新概念無舵面飛行器,這種飛行器通過環(huán)量控制技術(shù)和射流推力矢量技術(shù)來完成飛行器的俯仰、偏航和滾轉(zhuǎn)運動。環(huán)量控制技術(shù)應(yīng)用在機翼后緣上,高壓氣流從科恩達(Coanda)曲面上下兩端的槽道噴出形成射流,射流沿Coanda曲面繼續(xù)流動直至分離,從而改變自由來流的方向,起到類似副翼和襟翼的作用。射流推力矢量技術(shù)應(yīng)用在發(fā)動機尾噴口上,借助引入的二次噴射氣流去改變發(fā)動機尾流的流向,進而產(chǎn)生特定方向的推力。

    對環(huán)量控制技術(shù)的研究大致分為兩個階段,第一個階段主要發(fā)生在20世紀后30年,Englar[1-3]和Abramson[4]等主要研究Coanda的曲面形狀和噴口高度等幾何參數(shù)對機翼增升效果的影響,也有Loth[5]等研究將環(huán)量控制技術(shù)作為機械舵面的補充來產(chǎn)生更大的升力。這一階段的研究目的主要是利用環(huán)量控制技術(shù)提高機翼升力,實現(xiàn)短距起降,但Grumman A6等驗證機都存在諸如較高的發(fā)動機引氣量和阻力等問題,很難將這一技術(shù)真正應(yīng)用到實際中。第二個階段始于2000年左右,Warsop[6]、Frith[7]等利用環(huán)量控制技術(shù)產(chǎn)生類似舵面的控制力來控制飛行器。因為相比于達到明顯增升效果所需的引氣量,產(chǎn)生有效控制力所需的引氣量更小,僅占發(fā)動機總進氣量的1%~2%。近幾年,英國BAE系統(tǒng)公司聯(lián)合多所高校研發(fā)的Demon[8-9]和MAGMA[10]無舵面(虛擬舵面)無人機相繼試飛成功,說明環(huán)量控制技術(shù)有很大潛力取代機械舵面,為飛行器提供控制力。

    比起機械舵面,虛擬舵面機翼一體性強,沒有明顯的尖銳邊緣、縫隙和凸起物等雷達散射源,提高了隱身性能。虛擬舵面還去除了機械舵面所需的液壓作動器、傳力關(guān)鍵件和傳動接頭等組件,減輕機身重量的同時提高了機翼的維修性[11]。

    機械舵面依靠控制舵面偏轉(zhuǎn)角度來控制飛行,虛擬舵面依靠調(diào)節(jié)射流動量系數(shù)大小來控制飛行,為了研究這兩種完全不同的控制方式能否產(chǎn)生相同的控制效果,本文分別對不同舵偏角的機械舵面模型和不同射流動量系數(shù)的虛擬舵面模型進行了數(shù)值模擬,并與試驗結(jié)果進行對比,從機翼的受力等方面分析二者的異同。為了進一步研究虛擬舵面的環(huán)量控制系統(tǒng)的效率,本文對不同噴口高度的虛擬舵面進行了能耗計算和對比,并定義了等效升阻比的概念來比較不同噴口高度的虛擬舵面的氣動效率。

    1模型和數(shù)值方法

    1.1模型和網(wǎng)格

    本文使用的是半展長機翼模型[12],基準模型的俯視圖如圖1所示。機翼由內(nèi)、中、外三段組成,三段翼的4個端面A、B、C、D處的翼型如圖2所示,各個端面的弦長和展向分布在表1中列出,其中展向分布是指各端面到翼根的距離與展長的比值。機翼的前緣后掠角Λ= 35°,半展長b/2 = 688mm,平均氣動弦長MAC = 328.7mm。

    機械舵面和虛擬舵面模型均由基準模型得到。圖3是兩種模型的簡化示意圖,機械舵面的后緣是傳統(tǒng)的可上下偏轉(zhuǎn)的襟翼(見圖3(a)紅色部分),虛擬舵面的后緣是切向吹氣的Coanda曲面(見圖3(b)紅色部分)。襟翼和Coanda曲面只分布在外段翼上,后文的弦長c指外段翼的弦長,即c = 262.6mm。

    機械舵面模型的襟翼繞位于x = 0.772c處的旋轉(zhuǎn)軸轉(zhuǎn)動。為了便于生成高質(zhì)量結(jié)構(gòu)網(wǎng)格,對基準模型進行如下處理:(1)在襟翼兩端各留出一個3mm的縫隙,如圖4所示;(2)襟翼與主翼的翼面用樣條線光順連接,如圖5所示。前者避免了偏轉(zhuǎn)襟翼的兩端與主翼形成的“剪刀差”幾何不連續(xù)問題,有利于提升網(wǎng)格質(zhì)量;后者提高了機翼的一體性,便于生成O形網(wǎng)格。

    機械舵面整體為O形網(wǎng)格拓撲,網(wǎng)格總量約為4.3×106,第一層網(wǎng)格高度0.01mm,保證了第一層網(wǎng)格的y+< 1。機械舵面模型共有4個,分別為襟翼舵偏角θ= 0°,10°,20°和30°。

    虛擬舵面模型由基準模型修型得到。將傳統(tǒng)翼型修型為環(huán)量控制翼型通常有兩種方法,第一種方法是適當增大后半段翼型的厚度,變尖后緣為鈍后緣,這樣不僅滿足Coanda曲面的鈍后緣要求,也保證了充足的內(nèi)部空間來布置管道和閥等組件。此方法不改變弦長,但在加厚翼型時要保證翼面足夠光順,難度大。第二種方法是直接切去尖后緣,變成鈍后緣翼型,這種方法操作簡單,但缺點是弦長變小。本文的機械舵面本身為鈍后緣,可直接將后緣修形為Coanda曲面,如圖6所示。圖7是Coanda曲面的放大圖,藍色區(qū)域表示Coanda曲面,紅色區(qū)域表示射流噴口。

    Coanda曲面為半圓形,如圖8所示,半徑r = 5.36mm,噴口在x = 0.976c的位置,噴口高度h = 0.35mm,噴口寬度bj= 476.8mm。虛擬舵面整體為O形網(wǎng)格拓撲,圖9是翼型示意圖以及噴口位置放大圖。

    1.2數(shù)值方法和網(wǎng)格無關(guān)性驗證

    本文的所有算例使用ANSYS CFX求解器計算,湍流模型采用基于RANS方程的SST模型。美國國家航空航天局(NASA)[13]建立了環(huán)量控制翼型的試驗數(shù)據(jù)庫,目的是幫助驗證計算流體力學(CFD)程序的可靠性和準確性。

    本文用CC020-010EJ翼型的風洞試驗數(shù)據(jù)來驗證本文中使用的數(shù)值模擬方法。計算條件是Ma = 0.1,Cμ= 0.047,α= 0°,基于弦長的雷諾數(shù)Re = 5×105。射流入口的邊界條件為質(zhì)量流入口。對于給定的Cμ,通過以下公式可計算出對應(yīng)的質(zhì)量流率:

    對比CC020-010EJ翼型的表面壓力系數(shù)的計算值與試驗值[11],如圖11所示,發(fā)現(xiàn)計算結(jié)果與試驗結(jié)果高度吻合,三條不同網(wǎng)格密度的壓力系數(shù)曲線幾乎重合,說明本文中使用的數(shù)值方法有較高的可靠性和準確性。

    對不同網(wǎng)格密度的CC020-010EJ翼型的升阻力系數(shù)的計算結(jié)果與試驗結(jié)果進行比較(見表3),可以發(fā)現(xiàn)數(shù)值模擬的升阻力系數(shù)比試驗結(jié)果稍大,但總體而言計算值與試驗值很接近。隨著網(wǎng)格量增大,計算結(jié)果的精度更高,但加密網(wǎng)格所帶來的精度收益在逐漸減小,同時消耗更多的計算資源和時間。所以平衡計算精度和效率,本文選擇中等網(wǎng)格的節(jié)點分布作為兩種模型網(wǎng)格劃分的參考。

    2結(jié)果與討論

    2.1試驗結(jié)果與計算結(jié)果對比

    兩種模型的計算條件是自由來流速度V∞=20m/s,基于弦長的雷諾數(shù)Re = 5×105,迎角范圍從-4°到30°。機械舵面的舵偏角θ分別為0°,10°,20°和30°;虛擬舵面的射流動量系數(shù)分別為0,0.001,0.003,0.005,0.009,0.013,0.015,0.02,0.025,0.03,0.035,0.04。

    本文的試驗數(shù)據(jù)來自于西華大學的風洞試驗。圖12是機械舵面升阻力系數(shù)的計算與試驗結(jié)果對比。從升力系數(shù)曲線對比圖看,迎角在-4°~8°的線性段內(nèi)試驗結(jié)果與計算結(jié)果吻合程度很好,尤其是舵偏角為0°和10°時,二者幾乎重合,舵偏角增大到20°和30°時,試驗結(jié)果比計算結(jié)果偏大,但升力系數(shù)曲線的斜率保持相同;在大迎角狀態(tài)時,試驗結(jié)果表明機械舵面有較好的失速特性,沒有出現(xiàn)升力突然下降的情況,保證了飛機的失速改出能力,計算結(jié)果的最大升力系數(shù)和失速迎角均比試驗結(jié)果大,原因可能是使用的RANS湍流模型在模擬大分離流動時精度不高。從阻力系數(shù)曲線對比圖看,由于CFX軟件對阻力的計算精度不高,計算結(jié)果與試驗結(jié)果相差較大。

    計算結(jié)果顯示Cμ> 0.015時虛擬舵面的升阻力系數(shù)比機械舵面舵最大偏角θ= 30°的大很多,說明Cμ超過0.015后,其產(chǎn)生的控制力超過了機械舵面最大舵偏角產(chǎn)生的控制力,二者無法作對比,所以本節(jié)僅針對具有對比意義的0≤Cμ≤0.013區(qū)間進行分析。

    圖13是虛擬舵面升阻力系數(shù)曲線的計算與試驗結(jié)果對比。從升力系數(shù)曲線對比圖看,迎角在-4°~8°的線性段內(nèi)Cμ= 0,0.005和0.009的試驗結(jié)果與計算結(jié)果吻合程度很好,Cμ= 0.013的計算結(jié)果稍大于試驗結(jié)果;大迎角時,虛擬舵面與機械舵面的情況類似,計算結(jié)果的最大升力系數(shù)和失速迎角均比試驗結(jié)果大。從阻力系數(shù)曲線看,在較小迎角(6°以下)及失速以后(20°以上),計算值與試驗值的變化趨勢和吻合度較好,但是在6°至失速前的范圍內(nèi),二者相差較大,特別是試驗值呈下降趨勢,而計算值呈單調(diào)遞增趨勢。

    總體來看,CFD方法能準確計算出機械舵面和虛擬舵面的氣動力,二者的對比具有可靠性和準確性。

    2.2機械舵面與虛擬舵面對比

    2.2.1兩種舵面的等舵效關(guān)系

    通過調(diào)節(jié)射流動量系數(shù),虛擬舵面能否產(chǎn)生與機械舵面相同的控制力和力矩是本文研究的主要內(nèi)容。本文將4個不同舵偏角的機械舵面的升阻力系數(shù)和俯仰力矩系數(shù)曲線作為目標曲線,用不同Cμ的虛擬舵面的力和俯仰力矩系數(shù)曲線和目標曲線作對比,若三條曲線均吻合,則認為該Cμ具有與目標曲線舵偏角相同的舵效。為得到最好的匹配結(jié)果,增加計算Cμ在0~0.013范圍內(nèi)以0.001為間隔的之前沒有計算的狀態(tài)。

    數(shù)值模擬的匹配結(jié)果如圖14(a)~圖14(c)所示,圖14(d)~圖14(f)是對應(yīng)的試驗結(jié)果的匹配對比。從圖14(a)~圖14(c)中的匹配情況看,θ=0°與Cμ=0、θ=10°與Cμ=0.005、θ=20°與Cμ= 0.009和θ=30°與Cμ=0.012具有相同的舵效。除了虛擬舵面的失速迎角比機械舵面大一點,4組等舵效的升力系數(shù)和俯仰力矩系數(shù)曲線幾乎重合。從圖14(b)還發(fā)現(xiàn)θ=0°的阻力系數(shù)比Cμ=0的大,這是因為修形后虛擬舵面后緣更鈍,壓差阻力增大,其余三組的虛擬舵面的阻力系數(shù)均小于機械舵面,這說明虛擬舵面的環(huán)量控制技術(shù)具有減阻的作用。

    圖14(d)~圖14(f)的試驗結(jié)果對比也證明了上述的等舵效關(guān)系。從圖14(d)還發(fā)現(xiàn),機械舵面失速后沒有出現(xiàn)升力突然下降的情況,而虛擬舵面的升力下降非常明顯,說明前者的失速特性要優(yōu)于后者。

    2.2.2兩種舵面構(gòu)型的流場對比

    雖然等舵效的機械舵面和虛擬舵面具有相同的控制力和力矩,但是兩種舵面存在結(jié)構(gòu)差異,機翼表面的壓力分布和流線必然存在一些不同之處。圖16對比了兩種舵面的極限流線和表面壓力分布。從圖16(a)機械舵面上翼面的極限流線可以發(fā)現(xiàn),氣流在主翼面上為附著流動,在襟翼上的展向流動十分明顯,產(chǎn)生了流動分離現(xiàn)象,而圖16(b)虛擬舵面的對應(yīng)位置展向流動幾乎沒有,整個翼面均為附著流動。對比二者極限流線,說明虛擬舵面在根本上避免了襟翼大角度偏轉(zhuǎn)所帶來的流動分離現(xiàn)象,機翼弦向的載荷分布更好。

    機械舵面在襟翼上表面發(fā)生了流動分離,而虛擬舵面為附著流動,這兩種控制方式對機翼下游的流場是否產(chǎn)生不同的影響?圖17給出了機械舵面和虛擬舵面的三維流線圖。對比發(fā)現(xiàn),兩種舵面的下游流場比較類似,即外段翼兩端均有兩個大小不同的渦,內(nèi)側(cè)的渦小但向下偏折角度大。這更直觀地說明了兩種舵面的不同控制方式對下游的流場產(chǎn)生了基本等效的影響。

    2.3噴口高度對氣動效率和能耗的影響

    2.3.1功率系數(shù)和等效升阻比的定義

    機械舵面飛行器靠活動操縱面改變流場,活動操縱面由液壓作動器驅(qū)動,作動器所消耗的能量占比很小,甚至可忽略不計。與機械舵面飛行器不同,虛擬舵面飛行器需要消耗額外的能量產(chǎn)生所需的射流,射流向流場注入能量,起到控制作用。虛擬舵面飛行器在機動飛行或巡航時,氣壓泵持續(xù)不斷引氣所消耗的能量不可忽略,對其進行能耗分析十分必要。

    假設(shè)環(huán)量控制系統(tǒng)的氣壓泵入口連接自由流動的氣體,出口連接高壓腔入口。自由流動的氣體經(jīng)氣壓泵壓縮后進入高壓腔(即射流通道),高壓氣體從噴口噴出形成射流。氣壓泵所消耗的功率計算公式如下:

    圖19給出了不同噴口高度的虛擬舵面升阻力系數(shù)、俯仰力矩系數(shù)和等效升阻比曲線。從圖19(a)~圖19(c)可以看出,相同Cμ下噴口高度h的虛擬舵面的升阻力系數(shù)和俯仰力矩系數(shù)均比噴口高度2h的大,說明相同Cμ下噴口高度h的控制力強于噴口高度2h的。這是因為噴口高度越小,射流速度越大,射流附著Coanda曲面流動的距離越遠,機翼繞流環(huán)量越大,升力系數(shù)越大。但從圖19(d)發(fā)現(xiàn),噴口高度h的虛擬舵面的等效升阻比反而要比噴口高度2h的低,這是因為雖然前者的升力大于后者,但是前者的阻力和能耗也大于后者且所占比重大,導(dǎo)致前者的等效升阻比比后者小。

    為保證射流的噴口處速度相同,現(xiàn)對噴口高度2h的虛擬舵面的升阻力系數(shù)、俯仰力矩系數(shù)和等效升阻比曲線的橫坐標進行坐標變換,即Cμ= 1/2 Cμ。變化后的2h-scaled曲線如各圖中紅色虛線所示。從圖19(a)~圖19(c)可以看出,h與2h-scaled兩條曲線在Cμ< 0.01時的升阻力系數(shù)和俯仰力矩系數(shù)幾乎相同;Cμ> 0.01時2h-scaled曲線的值小于h曲線。但是圖19(d)表明h與2h-scaled的等效升阻比曲線幾乎重合,說明在射流速度相同的情況下,不同噴口高度的虛擬舵面的等效升阻比相同。

    2.3.3噴口高度對能耗的影響

    圖20給出了不同噴口高度的虛擬舵面所消耗的pc, jet隨Cμ的變化趨勢。從圖中可以發(fā)現(xiàn),隨著Cμ增大,pc, jet—Cμ曲線的斜率逐漸增大。這說明在大Cμ時,增加相同的Cμ需要更大的pc, jet增量,意味著更明顯的能量消耗提升。圖中兩條實線為噴口高度h和2h的虛擬舵面的原始pc, jet—Cμ曲線,發(fā)現(xiàn)相同Cμ下,小噴口需要耗費更大的功率。

    前文圖18已經(jīng)給出了相同Cμ情況下,噴口高度大,射流速度小。為了分析在相同射流速度下不同噴口高度的pc, jet變化,對噴口高度2h的虛擬舵面的pc, jet—Cμ曲線進行坐標變換,即Cμ =(1/2)Cμ,得到圖中2h-scaled的紅色虛線。比較噴口高度h和2h-scaled這兩條曲線,發(fā)現(xiàn)噴口高度2h的pc, jet比噴口高度h大,且差值隨Cμ增大而增大。說明在射流速度相同的情況下,噴口高度較大的氣壓泵比噴口高度較小的消耗更多的能量。

    圖19和圖20證明了相同Cμ下,噴口高度h的虛擬舵面控制力強但氣壓泵所需功率高,噴口高度2h的虛擬舵面控制力弱但氣壓泵所需功率低,說明噴口高度的設(shè)計受到機翼的控制力大小和氣壓泵所需功率大小的雙重限制。噴口高度的最優(yōu)化設(shè)計仍需要繼續(xù)深入研究。

    3結(jié)論

    本文對機械舵面和虛擬舵面兩種機翼進行數(shù)值模擬,并與試驗結(jié)果進行對比,進一步探尋了機械舵面舵偏角θ和虛擬舵面射流動量系數(shù)Cμ之間的等舵效關(guān)系。為了研究虛擬舵面的噴口高度對氣動效率和能耗的影響,又對噴口高度2h的模型進行計算,并且定義了一個新的計算虛擬舵面氣動效率的公式,對比了射流動能相同的情況下不同噴口高度的虛擬舵面的升阻力系數(shù)、俯仰力矩系數(shù)和等效升阻比曲線,還對比了不同噴口高度下的環(huán)量控制系統(tǒng)的能量消耗??梢缘贸鋈缦陆Y(jié)論:

    (1)通過對比試驗結(jié)果,發(fā)現(xiàn)用本文中的CFD方法能較準確模擬出機械舵面和虛擬舵面的真實氣動力。

    (3)機械舵面失速后沒有出現(xiàn)升力突然下降的情況,而虛擬舵面失速后升力下降非常明顯,說明機械舵面的失速特性要優(yōu)于虛擬舵面。

    (4)機械舵面的襟翼上表面會發(fā)生流動分離現(xiàn)象,而虛擬舵面上全部為附著流動,射流控制技術(shù)能避免機械舵面上的襟翼大角度偏折情況下的流動分離現(xiàn)象,且兩種控制方式對下游氣流的影響幾乎相同。

    (5)在射流速度相同的情況下,不同噴口高度的虛擬舵面的等效升阻比相同,但噴口高度h的氣壓泵所消耗的功率小于噴口高度2h的氣壓泵所消耗功率。

    (6)在Cμ相同時,噴口高度h的虛擬舵面控制力強但氣壓泵所需功率高,噴口高度2h的虛擬舵面控制力弱但氣壓泵所需功率低。噴口高度的設(shè)計受到機翼的控制力大小和氣壓泵所需功率大小的雙重限制。

    參考文獻

    [1]Englar R J,Williams R M. Design of a circulation control stern plane for submarine applications[R]. Bethesda:David W. Taylor Naval Ship Research and Development Center,1971.

    [2]Englar R J. Two-dimensional subsonic wind tunnel tests of two 15-percent thick circulation control airfoils[R]. Bethesda:David W. Taylor Naval Ship Research and Development Center,1971.

    [3]Englar R J. Experimental investigation of the high velocity Coanda wall jet applied to bluff trailing edge circulation control airfoils[R]. Bethesda:David W. Taylor Naval Ship Research and Development Center,1975.

    [4]AbramsonJ,RogersE.High-speedcharacteristicsof circulationcontrolairfoils[C]//21stAerospaceSciences Meeting,1983:265.

    [5]Loth J L,F(xiàn)anucci J B,Roberts S C. Flight performance of a circulation controlled STOL aircraft[J]. Journal of Aircraft,1976,13(3):169-173.

    [6]Warsop C,Crowther W J. Fluidic flow control effectors for flight control[J].AIAAJournal,2018,56(10):3808-3824.

    [7]Frith S,Wood N. Effect of trailing edge geometry on a circulationcontroldeltawing[C]//21stAIAAApplied Aerodynamics Conference,2003:3797.

    [8]Buonanno A. Aerodynamic circulation control for flapless flight control of an unmanned air vehicle[D]. Bedfordshire:Cranfield University,2009.

    [9]Fielding J P,Mills A,Smith H. Design and manufacture of the DEMONunmannedairvehicledemonstratorvehicle[J]. Proceedings of the Institution of Mechanical Engineers,Part G:Journal of Aerospace Engineering,2010,224(4):365-372.

    [10]Warsop C,Crowther W. NATO AVT-239 task group:Flight demonstration of fluidic flight controls on the magma subscale demonstrator aircraft[C]//AIAA Scitech 2019 Forum,2019:0282.

    [11]蔡琰.國外射流飛行控制技術(shù)發(fā)展及前景分析[J].航空科學技術(shù), 2019, 31(1): 85-86. Cai Yan. Development and prospect analysis of jet flight controltechnologyabroad.[J].AeronauticalScience&Technology, 2020, 31(1):85-86.(in Chinese)

    [12]徐悅,杜海,李巖,等.基于射流飛控技術(shù)的無操縱面飛行器研究進展[J].航空科學技術(shù), 2019, 30(4): 1-7. Xu Yue,Du Hai,Li Yan,et al.Progress research of fluidic flight control technology for flapless aircraft [J]. Aeronautical Science &Technology, 2019, 30(4): 1-7.(in Chinese)

    [13]Englar R,Jones G,Allan B,et al. 2-D circulation control airfoilbenchmarkexperimentsintendedforCFDcode validation[C]//47thAIAAAerospaceSciencesMeeting Including the New Horizons Forum and Aerospace Exposition,2009:902.

    [14]Lefebvre A M,Zha G C. Design of high wing loading compact electric airplane utilizing co-flow jet flow control[C]//53rd AIAAAerospace Sciences Meeting,2015:0772.

    [15]Xu H Y,Qiao C L,Yang H Q,et al. Active circulation control on the blunt trailing edge wind turbine airfoil[J]. AIAA Journal,2018,56(2):554-570.

    (責任編輯陳東曉)

    作者簡介

    付志杰(1994-)男,碩士研究生。主要研究方向:計算流體力學、流動控制。

    Tel:15620032372

    E-mail:zjfu@mail.nwpu.edu.cn

    許和勇(1980-)男,博士,教授。主要研究方向:計算流體力學、流動控制。

    Tel:15802935215

    E-mail:xuheyong@nwpu.edu.cn

    杜海(1985-)男,博士,講師。主要研究方向:空氣空力學、流動控制。

    Tel:15196686983

    E-mail:duhai2017@163.com

    王宇航(1991-)男,碩士,工程師。主要研究方向:飛行器設(shè)計。

    Tel:15810113662

    E-mail:yunmengjingtian@163.com

    徐悅(1979-)男,博士,研究員。主要研究方向:空氣動力學、流動控制。

    Tel:010-84929359

    E-mail:xuyue@cae.ac.cn

    Investigation on Flapless Wing Based on Circulation Control

    Fu Zhijie1,Xu Heyong1,*,Du Hai2,Wang Yuhang3,Xu Yue3

    1. National Key Laboratory of Science and Technology on Aerodynamic Design and Research,Northwestern

    Polytechnical University,Xian 700072,China

    2. Key Laboratory of Fluid and Power Machinery,Ministry of Education,Xihua University,Chengdu 610039,China 3. Chinese Aeronautical Establishment,Beijing 100012,China

    Abstract: Applying circulation control at the wing trailing edge could change the aerodynamics of the wing. The numerical simulations of the flap wing with different flap deflection anglesθand the flapless wing with different jet momentum coefficients Cμare conducted to investigate the control effect of the circulation control applied on the flapless aircraft. It is found that the control authority of the flap wing atθ=0°, 10°, 20°, 30°are equivalent to that of the flapless wing with Cμ= 0, 0.005, 0.009, 0.012 after comparing the lift, drag and moment coefficient curves between them. Andθand Cμare quadratic polynomial relations. Further, the numerical simulations of the flapless wing with different slot heights are conducted to access the aerodynamic efficiency and the energy expenditure for the flapless wing. It is found that the equivalent lift-to-drag ratios of the flapless wing with different slot heights are equal when they have the same jet velocity, however, the flapless wing with lager slot height needs relatively more power.

    Key Words: flapless wing; circulation control; equivalent control authority; energy expenditure; equivalent lift-todrag ratio

    猜你喜歡
    能耗
    EnMS在航空發(fā)動機試驗?zāi)芎目刂浦械膽?yīng)用實踐
    從能耗“雙控”向碳排放“雙控”轉(zhuǎn)變
    120t轉(zhuǎn)爐降低工序能耗生產(chǎn)實踐
    昆鋼科技(2022年2期)2022-07-08 06:36:14
    能耗雙控下,漲價潮再度來襲!
    寧夏議能耗“雙控”三年行動計劃
    中國氯堿(2021年10期)2021-12-21 06:18:02
    探討如何設(shè)計零能耗住宅
    水下飛起滑翔機
    日本先進的“零能耗住宅”
    華人時刊(2018年15期)2018-11-10 03:25:26
    基于ZigBee的大型公共建筑能耗采集器設(shè)計
    我國火力發(fā)電能耗狀況研究及展望
    两个人看的免费小视频| 午夜福利影视在线免费观看| 成人18禁高潮啪啪吃奶动态图| 亚洲精品av麻豆狂野| 午夜91福利影院| 一区二区三区国产精品乱码| 精品久久久精品久久久| 麻豆成人av在线观看| 国产一区在线观看成人免费| 国产午夜精品久久久久久| 亚洲av日韩精品久久久久久密| 久久精品亚洲av国产电影网| 狂野欧美激情性xxxx| 精品卡一卡二卡四卡免费| 每晚都被弄得嗷嗷叫到高潮| 女生性感内裤真人,穿戴方法视频| 黄片大片在线免费观看| 欧美久久黑人一区二区| 亚洲欧美激情在线| 长腿黑丝高跟| bbb黄色大片| xxxhd国产人妻xxx| 纯流量卡能插随身wifi吗| 岛国视频午夜一区免费看| 国产黄色免费在线视频| 男女高潮啪啪啪动态图| 变态另类成人亚洲欧美熟女 | 视频区欧美日本亚洲| 欧美色视频一区免费| 精品午夜福利视频在线观看一区| 亚洲国产精品合色在线| 欧美在线黄色| 夜夜爽天天搞| 久久午夜综合久久蜜桃| 身体一侧抽搐| 亚洲av日韩精品久久久久久密| 99精品久久久久人妻精品| 夜夜躁狠狠躁天天躁| 亚洲人成网站在线播放欧美日韩| 天堂动漫精品| 一进一出抽搐动态| svipshipincom国产片| 99久久精品国产亚洲精品| 男女下面进入的视频免费午夜 | 日韩视频一区二区在线观看| 天堂中文最新版在线下载| 电影成人av| a在线观看视频网站| 精品一区二区三区四区五区乱码| 亚洲av片天天在线观看| 亚洲av成人不卡在线观看播放网| 精品国产乱码久久久久久男人| 久久久精品国产亚洲av高清涩受| 欧美最黄视频在线播放免费 | 欧美激情 高清一区二区三区| 天天添夜夜摸| 又紧又爽又黄一区二区| 国产精品成人在线| 免费在线观看黄色视频的| 欧美中文日本在线观看视频| av天堂久久9| 精品久久久久久久久久免费视频 | 亚洲精品av麻豆狂野| 看黄色毛片网站| 黄色毛片三级朝国网站| 窝窝影院91人妻| 级片在线观看| 精品国内亚洲2022精品成人| 一本综合久久免费| 国产精品免费视频内射| 每晚都被弄得嗷嗷叫到高潮| 好男人电影高清在线观看| 欧美在线黄色| 视频区欧美日本亚洲| 欧美成人性av电影在线观看| 这个男人来自地球电影免费观看| 国产伦人伦偷精品视频| 欧美成人午夜精品| 日韩精品青青久久久久久| 国产在线观看jvid| 欧美日韩亚洲综合一区二区三区_| 国产欧美日韩一区二区精品| 国产亚洲欧美精品永久| 精品久久久久久久久久免费视频 | 变态另类成人亚洲欧美熟女 | 日韩成人在线观看一区二区三区| 18禁国产床啪视频网站| 国产一区二区三区在线臀色熟女 | 日韩精品中文字幕看吧| 一级a爱视频在线免费观看| 亚洲精品在线美女| 日本五十路高清| 国产精品国产av在线观看| 一二三四社区在线视频社区8| 电影成人av| av天堂久久9| 亚洲精品国产区一区二| 俄罗斯特黄特色一大片| 可以免费在线观看a视频的电影网站| 欧美国产精品va在线观看不卡| av欧美777| 久久天躁狠狠躁夜夜2o2o| 亚洲欧洲精品一区二区精品久久久| 亚洲色图av天堂| 欧美精品亚洲一区二区| 亚洲人成伊人成综合网2020| 久久久久国内视频| tocl精华| 亚洲欧美激情综合另类| 色综合欧美亚洲国产小说| 黄片播放在线免费| 色综合婷婷激情| 一二三四在线观看免费中文在| 久久婷婷成人综合色麻豆| 国产三级在线视频| av视频免费观看在线观看| 亚洲精品国产精品久久久不卡| 欧美成人性av电影在线观看| bbb黄色大片| 国产高清国产精品国产三级| 男男h啪啪无遮挡| 女性被躁到高潮视频| 日韩欧美国产一区二区入口| 欧美日韩亚洲综合一区二区三区_| 成人三级做爰电影| 亚洲欧美日韩另类电影网站| 19禁男女啪啪无遮挡网站| 精品国产一区二区久久| av免费在线观看网站| 午夜福利,免费看| 在线观看日韩欧美| 欧美日韩瑟瑟在线播放| 啦啦啦 在线观看视频| 色婷婷av一区二区三区视频| 午夜激情av网站| a级片在线免费高清观看视频| 后天国语完整版免费观看| 中文字幕高清在线视频| 精品人妻在线不人妻| 午夜精品在线福利| 精品国产国语对白av| 国产精品久久久av美女十八| 两性夫妻黄色片| 久久久久久久久中文| 巨乳人妻的诱惑在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美日韩高清在线视频| 一级a爱片免费观看的视频| av超薄肉色丝袜交足视频| 亚洲情色 制服丝袜| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产欧美一区二区综合| 亚洲av片天天在线观看| av在线播放免费不卡| 十八禁人妻一区二区| 国产成人av激情在线播放| 人人妻,人人澡人人爽秒播| 久久精品国产亚洲av高清一级| 搡老乐熟女国产| 久久久久久人人人人人| 日韩人妻精品一区2区三区| 久久久久久久久免费视频了| 国产xxxxx性猛交| a在线观看视频网站| 麻豆一二三区av精品| 无遮挡黄片免费观看| 中亚洲国语对白在线视频| 夜夜躁狠狠躁天天躁| 国产三级在线视频| 午夜福利在线观看吧| 亚洲 欧美一区二区三区| 久久人人精品亚洲av| 最近最新免费中文字幕在线| 丝袜美足系列| 国产欧美日韩一区二区三| 99香蕉大伊视频| 无限看片的www在线观看| 国产欧美日韩一区二区三| 人人妻人人澡人人看| 女人爽到高潮嗷嗷叫在线视频| 亚洲片人在线观看| 日本wwww免费看| 欧美日韩瑟瑟在线播放| 国产三级黄色录像| 黄色丝袜av网址大全| 日韩免费av在线播放| а√天堂www在线а√下载| 国产aⅴ精品一区二区三区波| 丰满饥渴人妻一区二区三| 国产精品免费一区二区三区在线| 中文字幕人妻丝袜制服| 一进一出抽搐动态| 搡老乐熟女国产| 性欧美人与动物交配| 中文字幕最新亚洲高清| 香蕉久久夜色| 欧美国产精品va在线观看不卡| 亚洲专区中文字幕在线| 超色免费av| 两性夫妻黄色片| 久久人妻av系列| 纯流量卡能插随身wifi吗| 露出奶头的视频| 国产成人av教育| 高潮久久久久久久久久久不卡| 国产伦一二天堂av在线观看| 久久欧美精品欧美久久欧美| 国产高清激情床上av| 久久99一区二区三区| 精品无人区乱码1区二区| 最近最新中文字幕大全电影3 | 国产99白浆流出| 国产主播在线观看一区二区| 性色av乱码一区二区三区2| 热99re8久久精品国产| 男女高潮啪啪啪动态图| 欧美性长视频在线观看| 韩国精品一区二区三区| 天堂动漫精品| 真人一进一出gif抽搐免费| 校园春色视频在线观看| 9191精品国产免费久久| 精品国产超薄肉色丝袜足j| 久久天堂一区二区三区四区| 99在线视频只有这里精品首页| 亚洲国产欧美网| 神马国产精品三级电影在线观看 | 亚洲 国产 在线| cao死你这个sao货| 成人永久免费在线观看视频| 午夜免费观看网址| 久久久久国产一级毛片高清牌| 麻豆久久精品国产亚洲av | 一区二区三区激情视频| 一个人观看的视频www高清免费观看 | 熟女少妇亚洲综合色aaa.| 午夜精品国产一区二区电影| 大香蕉久久成人网| 一边摸一边抽搐一进一出视频| 国产精品日韩av在线免费观看 | 一二三四在线观看免费中文在| 午夜福利在线免费观看网站| 亚洲精品久久午夜乱码| 另类亚洲欧美激情| 亚洲成a人片在线一区二区| 欧美国产精品va在线观看不卡| 欧美黑人精品巨大| 丝袜美腿诱惑在线| 欧美成狂野欧美在线观看| 美女扒开内裤让男人捅视频| 91精品国产国语对白视频| 欧美日韩亚洲高清精品| 亚洲精品美女久久久久99蜜臀| 久久精品成人免费网站| 欧美一级毛片孕妇| 岛国在线观看网站| 一区二区日韩欧美中文字幕| 中国美女看黄片| 亚洲情色 制服丝袜| 丰满饥渴人妻一区二区三| 午夜福利一区二区在线看| 久久精品成人免费网站| 久久热在线av| 不卡一级毛片| 国产精品二区激情视频| 18禁黄网站禁片午夜丰满| 精品卡一卡二卡四卡免费| 久久国产精品男人的天堂亚洲| 久久亚洲真实| 亚洲va日本ⅴa欧美va伊人久久| 国产三级黄色录像| 免费在线观看日本一区| 日本精品一区二区三区蜜桃| 法律面前人人平等表现在哪些方面| 亚洲av成人不卡在线观看播放网| 两性夫妻黄色片| 777久久人妻少妇嫩草av网站| 黑丝袜美女国产一区| 亚洲情色 制服丝袜| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美98| 黄色 视频免费看| av网站免费在线观看视频| 韩国av一区二区三区四区| 一二三四社区在线视频社区8| 一级片'在线观看视频| 99久久人妻综合| 久久午夜亚洲精品久久| 亚洲熟女毛片儿| 一夜夜www| 久久 成人 亚洲| 日韩欧美三级三区| 丰满人妻熟妇乱又伦精品不卡| 精品久久久久久,| 中文欧美无线码| 亚洲av五月六月丁香网| 精品免费久久久久久久清纯| 亚洲自偷自拍图片 自拍| 国产三级黄色录像| 久久国产亚洲av麻豆专区| 亚洲成人国产一区在线观看| 欧美另类亚洲清纯唯美| 国产亚洲精品第一综合不卡| 两人在一起打扑克的视频| 少妇粗大呻吟视频| 99精品久久久久人妻精品| 一区二区三区激情视频| 大香蕉久久成人网| 一区二区三区激情视频| 在线观看免费日韩欧美大片| 欧美成人午夜精品| 嫁个100分男人电影在线观看| 欧美日韩精品网址| 一级片免费观看大全| 色婷婷久久久亚洲欧美| 日日干狠狠操夜夜爽| 欧美日韩国产mv在线观看视频| 动漫黄色视频在线观看| 亚洲av第一区精品v没综合| 色综合婷婷激情| 在线天堂中文资源库| 久久久久久久久免费视频了| 欧美黑人欧美精品刺激| 欧美久久黑人一区二区| 男女床上黄色一级片免费看| 久久久久久久久免费视频了| 91精品三级在线观看| 91大片在线观看| 黄色女人牲交| av天堂久久9| 国产免费av片在线观看野外av| 自拍欧美九色日韩亚洲蝌蚪91| 男女午夜视频在线观看| 久久久国产成人精品二区 | 精品午夜福利视频在线观看一区| 国产一区在线观看成人免费| 熟女少妇亚洲综合色aaa.| 在线观看一区二区三区激情| 老司机午夜福利在线观看视频| 欧美日本亚洲视频在线播放| 国产成人欧美在线观看| 一级a爱片免费观看的视频| e午夜精品久久久久久久| 神马国产精品三级电影在线观看 | 精品国产一区二区三区四区第35| 91麻豆av在线| 在线观看免费视频网站a站| 中文字幕色久视频| 美女午夜性视频免费| 女人被狂操c到高潮| 亚洲伊人色综图| 日韩有码中文字幕| 男人舔女人的私密视频| 国产成人一区二区三区免费视频网站| 人人妻,人人澡人人爽秒播| 热99re8久久精品国产| 91字幕亚洲| av福利片在线| 好男人电影高清在线观看| 日本一区二区免费在线视频| 国产精品影院久久| 99久久国产精品久久久| av中文乱码字幕在线| 亚洲欧美激情综合另类| 欧美av亚洲av综合av国产av| 精品午夜福利视频在线观看一区| 国产高清激情床上av| 欧美中文综合在线视频| av视频免费观看在线观看| 国内久久婷婷六月综合欲色啪| 久久精品国产亚洲av香蕉五月| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品一二三| 亚洲精品国产色婷婷电影| 另类亚洲欧美激情| 久久久久久久久免费视频了| 19禁男女啪啪无遮挡网站| 亚洲精品av麻豆狂野| 女生性感内裤真人,穿戴方法视频| 中文字幕av电影在线播放| 国产亚洲欧美在线一区二区| 色播在线永久视频| 国产伦人伦偷精品视频| 水蜜桃什么品种好| 久久久久九九精品影院| 色哟哟哟哟哟哟| 别揉我奶头~嗯~啊~动态视频| 999精品在线视频| 午夜精品久久久久久毛片777| 12—13女人毛片做爰片一| 日韩有码中文字幕| 男女午夜视频在线观看| av网站在线播放免费| 一区福利在线观看| 久久影院123| 亚洲男人天堂网一区| xxxhd国产人妻xxx| 亚洲人成77777在线视频| 搡老熟女国产l中国老女人| 91av网站免费观看| 欧美日韩av久久| 一级片免费观看大全| 国产亚洲av高清不卡| 超色免费av| 欧美亚洲日本最大视频资源| 美女 人体艺术 gogo| 精品人妻1区二区| 亚洲一区高清亚洲精品| 国产高清激情床上av| 国产高清国产精品国产三级| 成年人免费黄色播放视频| 男男h啪啪无遮挡| 黄色视频,在线免费观看| 国产成+人综合+亚洲专区| 黄频高清免费视频| 黄色成人免费大全| 青草久久国产| 亚洲 国产 在线| e午夜精品久久久久久久| 美女大奶头视频| www.精华液| 日韩欧美一区二区三区在线观看| 夜夜躁狠狠躁天天躁| 琪琪午夜伦伦电影理论片6080| 国产1区2区3区精品| 欧美黄色片欧美黄色片| 亚洲国产中文字幕在线视频| 午夜老司机福利片| 超碰成人久久| 国产蜜桃级精品一区二区三区| 国产成人一区二区三区免费视频网站| 丰满饥渴人妻一区二区三| 日韩av在线大香蕉| 精品国产国语对白av| 成年人免费黄色播放视频| 亚洲第一欧美日韩一区二区三区| 伊人久久大香线蕉亚洲五| 久久久久久久午夜电影 | 日韩免费av在线播放| 丝袜人妻中文字幕| 亚洲欧美一区二区三区久久| 十八禁网站免费在线| 成人18禁在线播放| 在线视频色国产色| 精品一区二区三区四区五区乱码| 色哟哟哟哟哟哟| 久久久久九九精品影院| 午夜激情av网站| 男女下面插进去视频免费观看| 精品久久久精品久久久| 视频区欧美日本亚洲| 怎么达到女性高潮| 99热只有精品国产| svipshipincom国产片| 国产亚洲精品第一综合不卡| 在线观看免费日韩欧美大片| 美女国产高潮福利片在线看| 日本免费a在线| 又黄又粗又硬又大视频| 久久精品影院6| 欧美一级毛片孕妇| 久久久久亚洲av毛片大全| 成年人黄色毛片网站| 黑人巨大精品欧美一区二区mp4| 免费搜索国产男女视频| 久久精品亚洲精品国产色婷小说| 90打野战视频偷拍视频| 亚洲精品av麻豆狂野| 免费高清视频大片| xxx96com| 国产精品免费一区二区三区在线| 成人18禁在线播放| 一区福利在线观看| 免费在线观看亚洲国产| 老熟妇乱子伦视频在线观看| 国产精品自产拍在线观看55亚洲| 高潮久久久久久久久久久不卡| 国产av在哪里看| 91国产中文字幕| 精品国产超薄肉色丝袜足j| 桃红色精品国产亚洲av| 欧美成人午夜精品| 这个男人来自地球电影免费观看| 日韩中文字幕欧美一区二区| x7x7x7水蜜桃| av欧美777| 人人妻人人添人人爽欧美一区卜| 日韩视频一区二区在线观看| 国产精品1区2区在线观看.| 国产熟女午夜一区二区三区| 村上凉子中文字幕在线| 成人18禁在线播放| 视频区欧美日本亚洲| 又黄又爽又免费观看的视频| 成人黄色视频免费在线看| 村上凉子中文字幕在线| 最近最新中文字幕大全电影3 | 99国产极品粉嫩在线观看| 久久久国产成人精品二区 | 久久久久久久久中文| 久久久久久久精品吃奶| 女人被躁到高潮嗷嗷叫费观| 欧美成人免费av一区二区三区| 男女床上黄色一级片免费看| 一本大道久久a久久精品| 99riav亚洲国产免费| 99久久人妻综合| svipshipincom国产片| 女生性感内裤真人,穿戴方法视频| 免费日韩欧美在线观看| 男人操女人黄网站| 国产xxxxx性猛交| 夜夜躁狠狠躁天天躁| 日本免费a在线| 黄网站色视频无遮挡免费观看| 无人区码免费观看不卡| 亚洲片人在线观看| 91精品国产国语对白视频| 亚洲黑人精品在线| 欧美一区二区精品小视频在线| 成人三级做爰电影| 国产一区二区在线av高清观看| 女性被躁到高潮视频| 99在线视频只有这里精品首页| 天天影视国产精品| 久久青草综合色| 18禁美女被吸乳视频| 丰满饥渴人妻一区二区三| 免费不卡黄色视频| 亚洲视频免费观看视频| 91九色精品人成在线观看| 99riav亚洲国产免费| 新久久久久国产一级毛片| 在线播放国产精品三级| 成人亚洲精品一区在线观看| 少妇 在线观看| 极品人妻少妇av视频| 亚洲视频免费观看视频| 女人爽到高潮嗷嗷叫在线视频| 十八禁人妻一区二区| 欧美日本亚洲视频在线播放| 久久人人爽av亚洲精品天堂| 少妇裸体淫交视频免费看高清 | 久久天躁狠狠躁夜夜2o2o| 好男人电影高清在线观看| svipshipincom国产片| 亚洲av电影在线进入| 久热爱精品视频在线9| 性少妇av在线| 欧美黄色片欧美黄色片| 久久精品亚洲熟妇少妇任你| av福利片在线| 久热这里只有精品99| 色播在线永久视频| 桃红色精品国产亚洲av| 亚洲国产毛片av蜜桃av| 欧美最黄视频在线播放免费 | 一本大道久久a久久精品| 欧美乱妇无乱码| 久久人人精品亚洲av| 91老司机精品| 狠狠狠狠99中文字幕| 国产高清视频在线播放一区| √禁漫天堂资源中文www| 久久久精品欧美日韩精品| 夜夜看夜夜爽夜夜摸 | 亚洲全国av大片| www.熟女人妻精品国产| 欧美日韩亚洲综合一区二区三区_| 国产99白浆流出| 久久久国产一区二区| 国产欧美日韩精品亚洲av| 日韩精品中文字幕看吧| 三级毛片av免费| a级片在线免费高清观看视频| 老司机深夜福利视频在线观看| 男人操女人黄网站| 亚洲精品久久成人aⅴ小说| 精品卡一卡二卡四卡免费| 久久性视频一级片| 国产不卡一卡二| 91精品国产国语对白视频| 日韩大码丰满熟妇| 一级a爱片免费观看的视频| 国产1区2区3区精品| 又大又爽又粗| 无遮挡黄片免费观看| 99精品在免费线老司机午夜| 一级,二级,三级黄色视频| 国产av在哪里看| 精品国产乱子伦一区二区三区| 国产精品久久久久久人妻精品电影| 99久久人妻综合| 免费av中文字幕在线| 国产有黄有色有爽视频| 12—13女人毛片做爰片一| 黄色丝袜av网址大全| 久久人妻熟女aⅴ| 午夜视频精品福利| 黄色丝袜av网址大全| 性色av乱码一区二区三区2| 日本免费一区二区三区高清不卡 | 在线观看www视频免费| 国产亚洲欧美精品永久| 天天添夜夜摸| 亚洲国产欧美日韩在线播放| 高清在线国产一区| 精品免费久久久久久久清纯| 欧美+亚洲+日韩+国产| 亚洲情色 制服丝袜| 免费搜索国产男女视频| 新久久久久国产一级毛片| 国产精品电影一区二区三区| 午夜福利在线观看吧| 亚洲第一av免费看|