• <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ā)電能耗狀況研究及展望
    在线天堂最新版资源| 91九色精品人成在线观看| 能在线免费观看的黄片| 亚洲av.av天堂| 两性午夜刺激爽爽歪歪视频在线观看| 色综合欧美亚洲国产小说| 性欧美人与动物交配| 日韩av在线大香蕉| 丰满人妻一区二区三区视频av| 亚洲精品亚洲一区二区| 特级一级黄色大片| 怎么达到女性高潮| 亚洲av电影不卡..在线观看| 熟女人妻精品中文字幕| 一本久久中文字幕| 国产av不卡久久| 波多野结衣高清无吗| 亚洲人成网站高清观看| 国产淫片久久久久久久久 | 一二三四社区在线视频社区8| 在线观看一区二区三区| 最新中文字幕久久久久| 午夜a级毛片| 九九久久精品国产亚洲av麻豆| 午夜福利在线在线| 午夜福利在线在线| 少妇高潮的动态图| 18+在线观看网站| 高潮久久久久久久久久久不卡| 国产av麻豆久久久久久久| 看黄色毛片网站| 直男gayav资源| 国产成人啪精品午夜网站| 日韩欧美国产在线观看| 国产综合懂色| 十八禁国产超污无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 国产又黄又爽又无遮挡在线| 怎么达到女性高潮| 亚洲av不卡在线观看| 在线免费观看的www视频| 乱人视频在线观看| 久久天躁狠狠躁夜夜2o2o| 热99re8久久精品国产| 国产在视频线在精品| 欧美激情国产日韩精品一区| 网址你懂的国产日韩在线| 国产精品亚洲一级av第二区| 美女大奶头视频| 久久中文看片网| 成年女人永久免费观看视频| 国产精品日韩av在线免费观看| 日韩欧美免费精品| www.www免费av| 久久久国产成人免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 99久久成人亚洲精品观看| 久久久久国产精品人妻aⅴ院| 美女被艹到高潮喷水动态| 在线观看免费视频日本深夜| 欧美激情久久久久久爽电影| 亚洲乱码一区二区免费版| 99久久精品热视频| 中文字幕熟女人妻在线| 91久久精品国产一区二区成人| 一级作爱视频免费观看| 中出人妻视频一区二区| 窝窝影院91人妻| 日韩精品中文字幕看吧| 国产高清视频在线播放一区| 久久精品国产自在天天线| 九色成人免费人妻av| 99国产精品一区二区蜜桃av| 亚洲国产色片| 乱码一卡2卡4卡精品| 在线观看舔阴道视频| 国产毛片a区久久久久| 国产精品一区二区三区四区免费观看 | 此物有八面人人有两片| 丁香欧美五月| 国产精品一区二区三区四区久久| 亚洲性夜色夜夜综合| 哪里可以看免费的av片| 禁无遮挡网站| 国产成人aa在线观看| 99久国产av精品| 亚洲av.av天堂| 国产亚洲欧美在线一区二区| 午夜福利免费观看在线| 嫩草影院入口| 日韩欧美精品免费久久 | 麻豆av噜噜一区二区三区| 两个人视频免费观看高清| 黄色视频,在线免费观看| a级毛片免费高清观看在线播放| 我要搜黄色片| 丰满乱子伦码专区| 综合色av麻豆| 88av欧美| av欧美777| 无遮挡黄片免费观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲专区国产一区二区| 真人做人爱边吃奶动态| 在线播放国产精品三级| 最新在线观看一区二区三区| 亚洲欧美清纯卡通| 在线观看午夜福利视频| 日本 欧美在线| 国产精品女同一区二区软件 | 99久久精品热视频| 丰满的人妻完整版| 成人鲁丝片一二三区免费| 午夜激情欧美在线| 中出人妻视频一区二区| 欧美国产日韩亚洲一区| 国产成人啪精品午夜网站| 观看免费一级毛片| 亚洲中文日韩欧美视频| 一个人看视频在线观看www免费| 成人国产综合亚洲| 成人特级黄色片久久久久久久| 精品人妻偷拍中文字幕| 久久久成人免费电影| 国产亚洲欧美在线一区二区| 最近视频中文字幕2019在线8| 午夜老司机福利剧场| 十八禁人妻一区二区| 成人av在线播放网站| 国产精品嫩草影院av在线观看 | 1000部很黄的大片| 久久国产精品影院| 一区二区三区激情视频| av在线蜜桃| 亚洲精品亚洲一区二区| 最好的美女福利视频网| 国产精品av视频在线免费观看| 国产国拍精品亚洲av在线观看| 丰满的人妻完整版| 久久久久久久久大av| 婷婷六月久久综合丁香| 男女下面进入的视频免费午夜| 国产免费av片在线观看野外av| 我的女老师完整版在线观看| 亚洲国产欧洲综合997久久,| 国产精品三级大全| 欧美在线一区亚洲| 久久久国产成人免费| 国产精品伦人一区二区| 最近最新中文字幕大全电影3| 九色国产91popny在线| 国产伦精品一区二区三区四那| 国产伦一二天堂av在线观看| 国产精品一区二区三区四区免费观看 | 久久精品91蜜桃| 中文字幕av成人在线电影| 最新中文字幕久久久久| 欧美性感艳星| 午夜福利视频1000在线观看| 欧美极品一区二区三区四区| 午夜视频国产福利| 亚洲国产精品合色在线| 亚洲经典国产精华液单 | 哪里可以看免费的av片| 噜噜噜噜噜久久久久久91| 免费人成视频x8x8入口观看| 日韩欧美精品免费久久 | 变态另类成人亚洲欧美熟女| 午夜福利在线观看免费完整高清在 | 欧美xxxx黑人xx丫x性爽| 久久午夜福利片| 国产人妻一区二区三区在| 99久久精品国产亚洲精品| 男人狂女人下面高潮的视频| 亚洲第一电影网av| 成人精品一区二区免费| 免费av观看视频| 亚洲最大成人中文| 久久国产精品影院| 两个人的视频大全免费| 亚洲成人精品中文字幕电影| 一区福利在线观看| 日本精品一区二区三区蜜桃| 欧美在线一区亚洲| 久久性视频一级片| 一级av片app| 国产久久久一区二区三区| 日韩大尺度精品在线看网址| 一本综合久久免费| av在线蜜桃| 久久亚洲精品不卡| 国产精品自产拍在线观看55亚洲| 国产亚洲精品久久久com| 校园春色视频在线观看| 赤兔流量卡办理| 成人高潮视频无遮挡免费网站| 麻豆久久精品国产亚洲av| 丁香欧美五月| 少妇的逼好多水| 欧美成人性av电影在线观看| 国产亚洲欧美98| 午夜日韩欧美国产| 亚洲 国产 在线| 亚洲18禁久久av| 男插女下体视频免费在线播放| 老司机福利观看| 亚洲欧美日韩东京热| 国产野战对白在线观看| 变态另类丝袜制服| 91麻豆精品激情在线观看国产| 国产国拍精品亚洲av在线观看| 色播亚洲综合网| 成年免费大片在线观看| 9191精品国产免费久久| 在线天堂最新版资源| 国产精品一区二区性色av| 午夜两性在线视频| 亚洲一区二区三区色噜噜| 久久精品久久久久久噜噜老黄 | 哪里可以看免费的av片| 免费看a级黄色片| 亚洲av五月六月丁香网| 精品久久久久久久久久久久久| 国产日本99.免费观看| 97超级碰碰碰精品色视频在线观看| 亚洲成av人片免费观看| 可以在线观看的亚洲视频| 亚洲avbb在线观看| 久久国产精品人妻蜜桃| 99热精品在线国产| 精品久久久久久成人av| 老司机福利观看| 99久久精品国产亚洲精品| 亚洲国产色片| 桃红色精品国产亚洲av| 日韩高清综合在线| 亚洲精品久久国产高清桃花| 欧美日韩乱码在线| АⅤ资源中文在线天堂| 亚洲av.av天堂| 日本 av在线| 日韩成人在线观看一区二区三区| 亚洲第一区二区三区不卡| 看十八女毛片水多多多| 伊人久久精品亚洲午夜| 变态另类成人亚洲欧美熟女| 国产精品久久久久久久电影| 999久久久精品免费观看国产| 成人毛片a级毛片在线播放| 精品一区二区三区视频在线观看免费| 99久久九九国产精品国产免费| 免费观看精品视频网站| 老司机福利观看| 亚洲欧美清纯卡通| 亚洲av免费在线观看| 久久人人精品亚洲av| 麻豆成人av在线观看| 日本熟妇午夜| 国产69精品久久久久777片| av女优亚洲男人天堂| 性色avwww在线观看| 免费人成视频x8x8入口观看| 狠狠狠狠99中文字幕| 给我免费播放毛片高清在线观看| 麻豆av噜噜一区二区三区| 免费一级毛片在线播放高清视频| 白带黄色成豆腐渣| 免费av观看视频| 欧美日本视频| 免费在线观看日本一区| 特大巨黑吊av在线直播| 少妇人妻精品综合一区二区 | 噜噜噜噜噜久久久久久91| 亚洲av美国av| 日韩欧美国产一区二区入口| 日本成人三级电影网站| 成人永久免费在线观看视频| 国产伦一二天堂av在线观看| 真人做人爱边吃奶动态| 天堂影院成人在线观看| 国产成人啪精品午夜网站| 国产aⅴ精品一区二区三区波| h日本视频在线播放| 亚洲精品在线观看二区| 精品人妻偷拍中文字幕| 久久久久亚洲av毛片大全| 12—13女人毛片做爰片一| 国产精品久久久久久久电影| 欧美日本视频| 欧美乱色亚洲激情| 美女被艹到高潮喷水动态| 色av中文字幕| 国产亚洲欧美98| 亚洲欧美激情综合另类| 久99久视频精品免费| 国产麻豆成人av免费视频| 熟女人妻精品中文字幕| 亚洲人成电影免费在线| 一级黄片播放器| 日韩 亚洲 欧美在线| 最后的刺客免费高清国语| 91在线精品国自产拍蜜月| avwww免费| 日韩欧美 国产精品| 一个人看的www免费观看视频| 日日摸夜夜添夜夜添av毛片 | 村上凉子中文字幕在线| 少妇人妻一区二区三区视频| 日韩欧美一区二区三区在线观看| 精品午夜福利在线看| 露出奶头的视频| 18禁在线播放成人免费| 观看美女的网站| 日韩av在线大香蕉| 男插女下体视频免费在线播放| 精品无人区乱码1区二区| 夜夜夜夜夜久久久久| 久久伊人香网站| 精品一区二区三区av网在线观看| 午夜免费男女啪啪视频观看 | 亚洲精品亚洲一区二区| 无遮挡黄片免费观看| 久久久成人免费电影| 噜噜噜噜噜久久久久久91| 午夜免费激情av| 国产精品伦人一区二区| 人人妻人人看人人澡| 在线a可以看的网站| 亚洲av成人不卡在线观看播放网| www日本黄色视频网| 色综合亚洲欧美另类图片| 久久久久久久午夜电影| 欧美xxxx黑人xx丫x性爽| 久久性视频一级片| 国产成人aa在线观看| 精品人妻偷拍中文字幕| 欧美三级亚洲精品| 国产精品久久久久久亚洲av鲁大| 日韩欧美国产一区二区入口| 日本 av在线| 九九热线精品视视频播放| 欧美成人a在线观看| 国产精品美女特级片免费视频播放器| 尤物成人国产欧美一区二区三区| 真人做人爱边吃奶动态| 亚洲中文字幕日韩| 91av网一区二区| 黄片小视频在线播放| av在线天堂中文字幕| 99热6这里只有精品| 欧美+日韩+精品| 一a级毛片在线观看| 一本一本综合久久| 日韩欧美国产在线观看| 亚洲av成人精品一区久久| 国产精品野战在线观看| 女同久久另类99精品国产91| 91午夜精品亚洲一区二区三区 | 激情在线观看视频在线高清| 亚洲电影在线观看av| www日本黄色视频网| 变态另类丝袜制服| 一区二区三区四区激情视频 | 亚洲五月天丁香| 亚洲人成伊人成综合网2020| 又黄又爽又刺激的免费视频.| 国产白丝娇喘喷水9色精品| 国产成年人精品一区二区| 国产真实乱freesex| 免费av观看视频| 日韩欧美一区二区三区在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲人成网站高清观看| 动漫黄色视频在线观看| 国产高清视频在线观看网站| 美女 人体艺术 gogo| 在线播放无遮挡| av国产免费在线观看| 久久久久国内视频| 国产精品一区二区免费欧美| 日韩免费av在线播放| 一级毛片久久久久久久久女| 欧美日韩亚洲国产一区二区在线观看| 亚洲色图av天堂| 亚洲av免费在线观看| 在线看三级毛片| 婷婷精品国产亚洲av| 国产91精品成人一区二区三区| 国产真实乱freesex| 老司机深夜福利视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久午夜亚洲精品久久| 伊人久久精品亚洲午夜| 国产91精品成人一区二区三区| 热99re8久久精品国产| 欧美激情国产日韩精品一区| 久久99热这里只有精品18| 亚洲国产色片| 亚洲国产精品合色在线| 他把我摸到了高潮在线观看| 亚洲人成电影免费在线| 日韩欧美在线乱码| 精品乱码久久久久久99久播| 亚洲精品456在线播放app | 国内精品一区二区在线观看| 97超视频在线观看视频| 亚洲人成电影免费在线| 激情在线观看视频在线高清| 女生性感内裤真人,穿戴方法视频| 狂野欧美白嫩少妇大欣赏| 欧美潮喷喷水| 在线观看av片永久免费下载| 亚洲内射少妇av| 一级黄片播放器| 禁无遮挡网站| 日本三级黄在线观看| 全区人妻精品视频| 国产一区二区亚洲精品在线观看| 国产麻豆成人av免费视频| 久久久久免费精品人妻一区二区| 成人三级黄色视频| 一区二区三区免费毛片| 啦啦啦观看免费观看视频高清| 大型黄色视频在线免费观看| 亚洲性夜色夜夜综合| 国内久久婷婷六月综合欲色啪| 免费在线观看影片大全网站| 国产一区二区亚洲精品在线观看| 日本撒尿小便嘘嘘汇集6| 老司机深夜福利视频在线观看| av欧美777| 日本一二三区视频观看| av专区在线播放| 久久久久久久久大av| 精品人妻熟女av久视频| 特级一级黄色大片| 天堂√8在线中文| 精品午夜福利视频在线观看一区| 全区人妻精品视频| 久久久精品欧美日韩精品| 中文字幕人妻熟人妻熟丝袜美| 我的老师免费观看完整版| 亚洲狠狠婷婷综合久久图片| 免费黄网站久久成人精品 | 国产高清有码在线观看视频| 美女被艹到高潮喷水动态| 亚洲av.av天堂| 看十八女毛片水多多多| 欧美午夜高清在线| 精品熟女少妇八av免费久了| 亚洲一区二区三区不卡视频| 久久精品夜夜夜夜夜久久蜜豆| 99riav亚洲国产免费| 精品国内亚洲2022精品成人| 国产三级在线视频| 九九热线精品视视频播放| 午夜激情福利司机影院| 中文字幕av在线有码专区| 美女免费视频网站| 1024手机看黄色片| 成人高潮视频无遮挡免费网站| 国产精品三级大全| 国产aⅴ精品一区二区三区波| av中文乱码字幕在线| 嫁个100分男人电影在线观看| 欧美绝顶高潮抽搐喷水| 亚洲av美国av| 亚洲av不卡在线观看| 国产高清视频在线播放一区| 欧美激情在线99| 在线国产一区二区在线| 日本精品一区二区三区蜜桃| 国产久久久一区二区三区| 免费看a级黄色片| 久久久色成人| 看十八女毛片水多多多| 舔av片在线| 亚洲人成网站高清观看| 亚洲乱码一区二区免费版| 亚洲欧美日韩高清专用| 欧美不卡视频在线免费观看| a级毛片a级免费在线| 丰满人妻熟妇乱又伦精品不卡| 午夜福利18| 国产真实乱freesex| 亚洲人成网站在线播放欧美日韩| 天美传媒精品一区二区| 麻豆一二三区av精品| 男女床上黄色一级片免费看| 成熟少妇高潮喷水视频| 欧美+亚洲+日韩+国产| 最后的刺客免费高清国语| 老司机深夜福利视频在线观看| 有码 亚洲区| 真实男女啪啪啪动态图| 午夜a级毛片| 免费看日本二区| 亚洲av成人不卡在线观看播放网| 在现免费观看毛片| 亚洲专区国产一区二区| 日本a在线网址| 看片在线看免费视频| 在线免费观看的www视频| 99久国产av精品| 神马国产精品三级电影在线观看| 国内毛片毛片毛片毛片毛片| 中文字幕久久专区| 国产主播在线观看一区二区| 免费在线观看成人毛片| 不卡一级毛片| 性插视频无遮挡在线免费观看| 亚洲国产日韩欧美精品在线观看| 宅男免费午夜| 国产精品久久电影中文字幕| 村上凉子中文字幕在线| 成人性生交大片免费视频hd| 熟女电影av网| 精品久久久久久久末码| 在线十欧美十亚洲十日本专区| 丰满乱子伦码专区| 亚洲av成人精品一区久久| 少妇被粗大猛烈的视频| 亚洲avbb在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久精品国产99精品国产亚洲性色| 国产三级在线视频| 在线免费观看不下载黄p国产 | av专区在线播放| 我的女老师完整版在线观看| 国产精品久久电影中文字幕| 免费看a级黄色片| 亚洲人成电影免费在线| 国产精品一区二区性色av| 亚洲五月天丁香| 伊人久久精品亚洲午夜| 国产精品爽爽va在线观看网站| 亚洲欧美日韩高清在线视频| 久久婷婷人人爽人人干人人爱| 欧美丝袜亚洲另类 | 欧美乱色亚洲激情| 成人特级av手机在线观看| 国产精品电影一区二区三区| 欧美黄色片欧美黄色片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 男人舔女人下体高潮全视频| 成年女人毛片免费观看观看9| 亚洲精品日韩av片在线观看| 在线看三级毛片| 久久久久亚洲av毛片大全| 亚洲成a人片在线一区二区| 欧美日本亚洲视频在线播放| 97超级碰碰碰精品色视频在线观看| 成人国产一区最新在线观看| 日韩 亚洲 欧美在线| 91久久精品电影网| 很黄的视频免费| 每晚都被弄得嗷嗷叫到高潮| 久久性视频一级片| 国产精品久久久久久久电影| 欧美xxxx黑人xx丫x性爽| av福利片在线观看| 最后的刺客免费高清国语| 在线观看66精品国产| 一区二区三区免费毛片| 国产极品精品免费视频能看的| 欧美乱妇无乱码| ponron亚洲| 真人一进一出gif抽搐免费| 成年女人看的毛片在线观看| 青草久久国产| 色哟哟哟哟哟哟| 国产精品美女特级片免费视频播放器| 一区二区三区高清视频在线| 男人和女人高潮做爰伦理| 我的老师免费观看完整版| 最近在线观看免费完整版| 成年女人看的毛片在线观看| 国产亚洲av嫩草精品影院| 亚洲av二区三区四区| 久久亚洲真实| 老熟妇乱子伦视频在线观看| 日本 欧美在线| 一级黄色大片毛片| 一本久久中文字幕| 精品人妻偷拍中文字幕| 国产美女午夜福利| 嫩草影视91久久| aaaaa片日本免费| 欧美日本亚洲视频在线播放| 精品午夜福利视频在线观看一区| 一a级毛片在线观看| 色尼玛亚洲综合影院| 国产精品一及| 亚洲在线自拍视频| 亚洲av日韩精品久久久久久密| 亚洲欧美日韩高清在线视频| 国产精品一区二区三区四区久久| 免费在线观看日本一区| 久久久久久大精品| 亚洲国产日韩欧美精品在线观看| 欧美性猛交╳xxx乱大交人| 在线播放无遮挡| 久久久久久久午夜电影| 中文字幕av在线有码专区| 听说在线观看完整版免费高清| 国产精品99久久久久久久久| 丝袜美腿在线中文| 国产老妇女一区| 亚洲av免费高清在线观看| 亚洲人成网站在线播放欧美日韩| 男人的好看免费观看在线视频| 亚洲一区二区三区色噜噜| 欧美日韩福利视频一区二区| 自拍偷自拍亚洲精品老妇| 亚洲第一欧美日韩一区二区三区|