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

    基于翼尖渦物理特征的誘導(dǎo)阻力減阻機(jī)制實驗研究

    2017-11-07 10:06:35黃文濤顧定一
    實驗流體力學(xué) 2017年5期
    關(guān)鍵詞:環(huán)量渦的翼尖

    黃文濤,向 陽,王 笑,劉 洪,顧定一

    (上海交通大學(xué) 航空航天學(xué)院,上海 200240)

    基于翼尖渦物理特征的誘導(dǎo)阻力減阻機(jī)制實驗研究

    黃文濤,向 陽,王 笑,劉 洪*,顧定一

    (上海交通大學(xué) 航空航天學(xué)院,上海 200240)

    本文通過風(fēng)洞實驗研究了翼尖渦的物理特征以及誘導(dǎo)阻力的減阻機(jī)制。實驗中利用3DPIV(三維粒子圖像測速技術(shù))技術(shù)得到了翼尖渦的物理特征,并基于本文提出并設(shè)計的翼尖氣動力測量裝置,得到了機(jī)翼翼尖處的誘導(dǎo)阻力。實驗結(jié)果表明,機(jī)翼翼尖渦的無量綱環(huán)量會隨機(jī)翼迎角及風(fēng)速的增大而增大。翼尖渦無量綱環(huán)量的減小以及翼尖渦與機(jī)翼之間距離的增大都會引起誘導(dǎo)阻力的減小。具體而言,通過抑制翼尖渦的無量綱環(huán)量,增加翼尖渦與主機(jī)翼之間的距離,減小翼尖渦與機(jī)翼之間的相互作用,實現(xiàn)機(jī)翼翼尖誘導(dǎo)阻力的減阻。

    翼尖渦;物理特征;誘導(dǎo)阻力;環(huán)量;減阻

    0 引 言

    在我國大型客機(jī)氣動設(shè)計中,民機(jī)的氣動減阻一直是目前面臨的主要挑戰(zhàn)之一,具有很重要的意義。在飛機(jī)翼尖處產(chǎn)生的誘導(dǎo)阻力是飛機(jī)飛行時受到阻力的重要組成部分,在巡航狀態(tài)時約占飛機(jī)總阻力的40%,有很大的減阻空間[1]。根據(jù)誘導(dǎo)阻力的定義可知,飛機(jī)翼尖處的誘導(dǎo)阻力是由于飛機(jī)翼尖處的翼尖渦誘導(dǎo)引起的[2-3]。阻礙或者破壞翼尖渦結(jié)構(gòu)形成,設(shè)法消除翼尖處的集中渦,使洗流具有均勻的分布是目前飛機(jī)誘導(dǎo)阻力減阻的主要措施[4-7]。因此,本文將結(jié)合翼尖渦的物理特征和機(jī)翼翼尖的誘導(dǎo)阻力研究翼尖誘導(dǎo)阻力的減阻方法,并揭示相應(yīng)的減阻機(jī)制,進(jìn)而指導(dǎo)飛機(jī)機(jī)翼的翼尖設(shè)計。

    實驗和數(shù)值計算[8-10]是現(xiàn)階段研究翼尖渦引起的誘導(dǎo)阻力減阻的2種主要方法。C.S.Lee等人[11-12]通過在翼尖處增加噴流裝置進(jìn)行實驗,發(fā)現(xiàn)翼尖渦的渦核會隨著噴流方向的不同而改變,從而導(dǎo)致翼尖處誘導(dǎo)阻力的變化。Margaris和Gursul[13]等人研究了噴流裝置的噴口位置對翼尖渦的影響,發(fā)現(xiàn)不同的噴口位置對翼尖渦結(jié)構(gòu)的改變是不同的。F.M.Catalano等人利用熱線風(fēng)速儀進(jìn)行實驗研究,發(fā)現(xiàn)降低翼尖渦強(qiáng)度會削弱翼尖渦引起的下洗流動,從而降低翼尖渦引起的誘導(dǎo)阻力[14]。Luis Falcao和德國Saarland大學(xué)的Christian Boiler等人通過數(shù)值研究發(fā)現(xiàn),增強(qiáng)翼尖渦的耗散能力,會提高相應(yīng)的翼尖誘導(dǎo)阻力的減阻效果[15-16]。由此,根據(jù)前人的研究可知,翼尖渦的物理特征與誘導(dǎo)阻力之間存在緊密的聯(lián)系[17-22]。本文將通過風(fēng)洞實驗進(jìn)一步研究它們之間的關(guān)系,揭示誘導(dǎo)阻力減阻機(jī)制。

    實驗中,為了測量翼尖處的誘導(dǎo)阻力,本文提出并設(shè)計了一套翼尖氣動阻力的測量裝置。裝置中傳感器測量的阻力方向始終與風(fēng)洞軸線平行,不會隨機(jī)翼迎角的變化而變化,即始終平行于來流方向,這種測量阻力的方法可以不通過坐標(biāo)變換直接測量阻力,從而避免引入由角度測量偏差引起的阻力誤差。同時,與翼尖渦的物理特征相結(jié)合,進(jìn)而分析基于翼尖渦物理特征的翼尖減阻原則,為以后機(jī)翼翼尖的設(shè)計提供基礎(chǔ)數(shù)據(jù)和理論指導(dǎo)。

    綜上所述,本文的工作具體分為以下幾個部分。第一部分是研究概述,第二部分是實驗裝置的介紹,第三部分是翼尖渦物理特征與誘導(dǎo)阻力的測量與計算,第四部分是基于翼尖渦物理特征的誘導(dǎo)阻力減阻機(jī)制研究,最后一部分是對本文工作的總結(jié)。

    1 實驗裝置

    本實驗是在上海交通大學(xué)的低速回流式風(fēng)洞中進(jìn)行的。風(fēng)洞收縮比為6.25,試驗段截面形狀為矩形,截面尺寸為1.2m×0.9m,長2m。實驗中的最低風(fēng)速為10m/s,最高風(fēng)速為70m/s。在試驗段上游,通過收縮段前面的蜂窩網(wǎng)狀結(jié)構(gòu)整流,以降低實驗段中的湍流度。為了提高流場品質(zhì),實驗前需要對風(fēng)洞流場進(jìn)行校驗。通過熱線風(fēng)速儀測量得到風(fēng)洞自由來流的湍流度為0.05%,圖1為風(fēng)洞實驗裝置簡圖。

    實驗中采用簡化的機(jī)翼模型,機(jī)翼翼型選為NACA0020。機(jī)翼采用鋁合金材料加工,弦長為25cm,如圖2所示。通過改變機(jī)翼的展長高度來滿足實驗所需的展弦比。展向長度范圍為25~65cm,即展弦比為1~2.69。機(jī)翼模型安裝在有刻度的水平轉(zhuǎn)盤上,水平轉(zhuǎn)盤固定在風(fēng)洞底面上,轉(zhuǎn)動水平轉(zhuǎn)盤即可變換機(jī)翼模型的迎角。為了減小前緣渦的分離,實驗中采用的迎角范圍為0°~12°[23]。其中,0°迎角機(jī)翼模型的翼尖位于風(fēng)洞中心軸線上,如圖2所示。進(jìn)行翼尖渦的流場測試時,風(fēng)速V∞=20m/s,以模型弦長為參考長度的雷諾數(shù)Re=3×105。

    圖1 風(fēng)洞實驗裝置示意圖Fig.1 Experimental setup in wind tunnel

    圖2 機(jī)翼模型的風(fēng)洞安裝示意圖Fig.2 Installation of wing model in the wind tunnel

    1.1三維PIV(3DPIV)系統(tǒng)

    實驗中采用三維PIV系統(tǒng)進(jìn)行圖像采集。其中,三維PIV系統(tǒng)是利用2臺型號為Imager-Pco的高分辨率CCD相機(jī)進(jìn)行拍攝,相機(jī)的最高分辨率可達(dá)4032pixel×2688pixel,圖像采集速率為1幀/s。相機(jī)鏡頭采用NIKKOR 100mm 1∶1.4D定焦鏡頭,并通過532nm波長的帶通濾光鏡來提高圖像信噪比。2臺相機(jī)并不正對著激光面,而是與之成一定的夾角。實驗中,相機(jī)與鏡頭與激光面之間的傾斜角必須滿足Scheimpflug條件[24]。在本實驗中,激光片與CCD相機(jī)記錄的平面相距1000mm,2個相機(jī)之間的夾角為60°。此時,三維PIV的拍攝窗口大約200mm×200mm。激光器則選用雙脈沖Nd:YAG激光器,激光面與試驗段高透明度玻璃垂直,激光光片厚度為0.5~3mm可調(diào)。CCD相機(jī)與雙脈沖Nd:YAG激光共同連接到同步器的延遲系統(tǒng)上,以此來實現(xiàn)激光發(fā)射與圖像采集之間的同步。通過Insight軟件對3DPIV拍攝的流場進(jìn)行處理,得到流場中的翼尖渦結(jié)構(gòu)信息。

    風(fēng)洞中,示蹤粒子選為乙二醇加熱后產(chǎn)生的煙霧,發(fā)煙器置于機(jī)翼模型下游的風(fēng)洞管道內(nèi),粒徑為1μm≤d≤5μm。為了保證粒子布撒的均勻性,粒子將在風(fēng)洞中循環(huán)幾次,從而可以減輕粒子不均勻?qū)α鲌龅挠绊憽y量時,片光垂直來流方向且距離翼尖后緣100mm,2個相機(jī)與片光成30°夾角,記錄片光內(nèi)示蹤粒子的運動軌跡。每個實驗狀態(tài)均采集300幀圖像,進(jìn)行時均處理后,可得到該狀態(tài)下的流場信息。

    1.2誘導(dǎo)阻力測量系統(tǒng)

    Gharib等人[19]在研究平板翼起動過程中翼尖渦對阻力的影響時,根據(jù)平板翼表面的流動結(jié)構(gòu)將平板翼分成2個區(qū)域,即準(zhǔn)二維區(qū)域和翼尖渦影響的三維區(qū)域。本文基于該思想,設(shè)計了一套機(jī)翼翼尖氣動力測量裝置,通過該裝置測量機(jī)翼翼尖處的翼尖渦引起的氣動阻力。裝置模型如圖3所示,該氣動力測量裝置將機(jī)翼分為2部分,主機(jī)翼部分與翼尖部分。傳感器不與主機(jī)翼部分相連,而是與翼尖部分直接連接,由此可測量翼尖處的氣動阻力。

    圖3 測力裝置模型圖Fig.3 Model of force measuring system

    實驗中,主機(jī)翼和傳感器分別安裝在底板上,兩者不相連。在改變主機(jī)翼迎角時,同時改變傳感器與翼尖的相對位置,使機(jī)翼翼尖與主機(jī)翼迎角相同。由于傳感器固定在底板上并不與主機(jī)翼相連接,傳感器始終保持不動,即其測量方向始終不變,可以保證測量得到的氣動力為機(jī)翼受到的氣動阻力,此時,該氣動阻力為機(jī)翼受到的總阻力。在來流速度較低的情況下,機(jī)翼的總阻力可以分為翼型阻力和誘導(dǎo)阻力。其中,翼型阻力可以通過測量不同展弦比下準(zhǔn)二維區(qū)域的阻力計算得到。進(jìn)行一個假設(shè):翼型阻力在三維翼尖渦影響區(qū)域不發(fā)生變化。因此,通過傳感器測量的總阻力減去上面計算求得的準(zhǔn)二維翼型阻力即可以獲得翼尖處翼尖渦引起的誘導(dǎo)阻力。

    為了能夠精確測量翼尖處的氣動力,在實驗中選用HBM的高精度的SP4M單分量傳感器測量與來流平行的氣動力,即氣動力中的阻力。傳感器量程為5kg,測量精度可達(dá)萬分之二,所以可以很好地滿足實驗的需要。單分量傳感器測量得到的翼尖處的氣動力信號通過NI公司型號為NI PXIe-1071的數(shù)據(jù)采集系統(tǒng)來采集,采樣頻率為200Hz。數(shù)據(jù)采集系統(tǒng)與電腦相連,將數(shù)據(jù)讀入電腦中。之后將采集得到的數(shù)據(jù)通過MATLAB程序進(jìn)行再一次的濾波處理。

    2 機(jī)翼翼尖渦的物理特征與誘導(dǎo)阻力

    在實驗中,風(fēng)洞的風(fēng)速范圍設(shè)定在10~70m/s,機(jī)翼的弦長C=25cm,展弦比為2,即機(jī)翼展長為50cm,機(jī)翼的迎角為4°、6°、8°、10°和12°。為了分析近壁面翼尖渦結(jié)構(gòu)的生長演化過程,將用3DPIV來進(jìn)行實驗,記錄尾跡中距離機(jī)翼尾緣不同弦長處的近壁面翼尖渦結(jié)構(gòu)。

    機(jī)翼翼尖渦的無量綱環(huán)量:

    式中:Γ為翼尖渦結(jié)構(gòu)的環(huán)量;U為風(fēng)洞的風(fēng)速;L為翼尖渦的特征尺寸。

    2.1機(jī)翼翼尖渦的物理特征

    圖4給出了迎角為6°時距離機(jī)翼尾緣不同截面處的翼尖渦結(jié)構(gòu)。圖中的截面從右至左分別為距離機(jī)翼尾緣1倍弦長、2倍弦長、3倍弦長和4倍弦長截面處的翼尖渦結(jié)構(gòu)。圖4左上角中的rc為翼尖渦渦核與機(jī)翼尾緣翼尖處之間的無量綱距離。

    圖4 6°迎角下機(jī)翼尖渦云圖Fig.4 Contour of the tip vortex of wing with the angle of attack at 6°

    從圖4可以看到,機(jī)翼翼尖渦向下游傳輸?shù)倪^程中也會不斷地向翼根處移動,同時翼尖渦的強(qiáng)度在脫落之后基本保持不變,即翼尖渦已經(jīng)處于穩(wěn)定階段。

    圖5給出了風(fēng)速U=15m/s時,在D/C=4截面處,不同迎角下的翼尖渦結(jié)構(gòu)的速度場和渦量云圖。

    圖5 不同機(jī)翼迎角下的翼尖渦的瞬時速度場和渦量云圖Fig.5 Instantaneous velocity (arrows) and vorticity (color) fields of the tip vortex of wing with different angles of attack

    從圖5中可以看到,翼尖渦會隨著機(jī)翼迎角的增大而增大。圖6給出了翼尖渦對應(yīng)的無量綱環(huán)量變化曲線圖。從圖6中可以看出,翼尖渦的無量綱環(huán)量隨著機(jī)翼迎角的增大而不斷地增大。在迎角較小時(在6°~10°變化),翼尖渦的無量綱環(huán)量隨迎角的增加較緩。在迎角較大時(從10°到12°),翼尖渦的無量綱環(huán)量隨迎角的增加變得較大。

    圖6 不同迎角對翼尖渦的無量綱環(huán)量的影響Fig.6 Non-dimensional circulation of vortex structure on wingtip influenced by the angle of attack

    圖7給出了迎角為6°時,在D/C=4截面處,不同風(fēng)速下翼尖渦結(jié)構(gòu)的速度場和渦量云圖。

    從圖7中可以看到,翼尖渦隨著迎角的增大不斷地增大。圖8給出了翼尖渦的無量綱環(huán)量隨雷諾數(shù)的變化規(guī)律。從圖中可以看到,隨著雷諾數(shù)的增大,翼尖渦的無量綱環(huán)量也會隨之增大。這主要是由于隨著雷諾數(shù)的增大,翼尖處脫落的渦結(jié)構(gòu)更強(qiáng)。

    圖7 不同風(fēng)速下的翼尖渦的瞬時速度場和渦量云圖Fig.7 Instantaneous velocity (arrows) and vorticity (color) fields of the tip vortex of the plate wing with different velocities

    圖8 不同雷諾數(shù)對翼尖渦的無量綱環(huán)量的影響Fig.8 Non-dimensional circulation of vortex structure on wingtip influenced by the Reynolds number

    2.2翼尖誘導(dǎo)阻力的測量與計算

    為了得到翼尖部分的翼尖阻力,將機(jī)翼根據(jù)其表面流動分為準(zhǔn)二維區(qū)域和翼尖三維區(qū)域[19,25]。準(zhǔn)二維區(qū)域氣動力的改變呈線性變化,翼尖三維區(qū)域的氣動力呈非線性變化。實驗中將改變翼尖部分的高度,測量不同高度翼尖的氣動阻力,由此確定翼尖處受到的氣動阻力。不過,此時測量得到的翼尖處的氣動阻力包括誘導(dǎo)阻力和型阻2部分。此外,根據(jù)文獻(xiàn)[25],型阻會隨著測量高度的增大逐漸增大,誘導(dǎo)阻力則不會隨高度增加而增大。那么,可根據(jù)測量得到的不同高度翼尖處的氣動阻力估算出相應(yīng)的誘導(dǎo)阻力。

    圖9給出了在不同迎角下翼尖渦引起的誘導(dǎo)阻力的變化規(guī)律。

    圖9 迎角對翼尖誘導(dǎo)阻力的影響Fig.9 Induced drag of wingtip influenced by the angle of attack

    圖9中橫坐標(biāo)為機(jī)翼迎角,縱坐標(biāo)為機(jī)翼翼尖處的誘導(dǎo)阻力。圖中紅色實線、綠色虛線和藍(lán)色點畫線則分別給出了風(fēng)速U=20、25和35m/s時翼尖渦引起的誘導(dǎo)阻力的變化曲線。從圖中可以看到,不同風(fēng)速下的翼尖渦引起的誘導(dǎo)阻力隨機(jī)翼迎角的變化規(guī)律一致,均是隨著迎角的增大而不斷地增大。

    圖10給出了在不同風(fēng)速下,翼尖渦引起的誘導(dǎo)阻力的變化規(guī)律。圖中紅色實線、綠色虛線、藍(lán)色虛點線以及黑色點畫線則分別給出了迎角AOA=6°、8°、10°和12°時翼尖渦引起的誘導(dǎo)阻力的變化曲線。

    圖10 風(fēng)速對翼尖誘導(dǎo)阻力的影響Fig.10 Induced drag of wingtip influenced by airflow speed

    圖10中橫坐標(biāo)為風(fēng)速,縱坐標(biāo)為機(jī)翼翼尖處的誘導(dǎo)阻力。從圖中可以看到,不同機(jī)翼迎角下的翼尖渦引起的誘導(dǎo)阻力隨風(fēng)速的變化規(guī)律一致,均是隨著風(fēng)速的增大而不斷地增大,這主要是由于隨著風(fēng)速的增大,翼尖渦不斷增強(qiáng),從而使得翼尖渦與主機(jī)翼間的相互作用變大,導(dǎo)致了翼尖處誘導(dǎo)阻力的增大。

    3 基于翼尖渦物理特征的誘導(dǎo)阻力減阻機(jī)制

    根據(jù)上面的研究,通過3DPIV技術(shù)獲得了翼尖渦的物理特征,此外,利用傳感器測量和計算得到了翼尖處的誘導(dǎo)阻力。基于這2部分的內(nèi)容將進(jìn)一步分析研究翼尖渦的物理特征對誘導(dǎo)阻力的影響以及相應(yīng)的誘導(dǎo)阻力減阻機(jī)制。

    圖11給出了翼尖處的誘導(dǎo)阻力系數(shù)CDinduce與翼尖渦的無量綱環(huán)量Γ*的關(guān)系曲線圖。其中,誘導(dǎo)阻力系數(shù)CDinduce可表示為:

    式中:F為翼尖處的誘導(dǎo)阻力;U為風(fēng)洞的風(fēng)速;ρ為風(fēng)洞內(nèi)的空氣密度。

    圖11中可以看到,機(jī)翼翼尖誘導(dǎo)阻力系數(shù)隨著翼尖渦的無量綱環(huán)量Γ*增大而不斷地增大。在近壁面處,隨著翼尖渦無量綱環(huán)量的增大,翼尖渦與主機(jī)翼之間的相互作用變強(qiáng),引起的誘導(dǎo)阻力也會隨之增大。即,翼尖處的誘導(dǎo)阻力系數(shù)與翼尖渦的無量綱環(huán)量呈正相關(guān),即

    圖11 翼尖渦無量綱環(huán)量對誘導(dǎo)阻力系數(shù)的影響Fig.11 Coefficient of induced drag of wingtip influenced by non-dimensional circulation

    由此,與傳統(tǒng)的減阻方法相同,翼尖渦的無量綱環(huán)量的減小,即抑制翼尖渦的生長也將作為機(jī)翼翼尖誘導(dǎo)阻力的減阻目標(biāo)之一。

    式中:rc為翼尖渦渦核與機(jī)翼尾緣翼尖處之間的距離,如圖4所示。rx、ry和rz為翼尖渦渦核與機(jī)翼尾緣翼尖處在3個方向上的距離,C為機(jī)翼弦長。

    除了抑制翼尖渦的強(qiáng)度以外,控制翼尖渦與主機(jī)翼之間的距離也將作為機(jī)翼翼尖誘導(dǎo)阻力的減阻目標(biāo)之一。

    不過,當(dāng)翼尖渦與主機(jī)翼之間的距離過大時,由于翼尖渦此時遠(yuǎn)離主機(jī)翼,它對誘導(dǎo)阻力系數(shù)CDinduce的影響不敏感。那么

    4 結(jié) 論

    本文通過風(fēng)洞實驗,研究了翼尖渦的物理特征以及機(jī)翼誘導(dǎo)阻力的減阻原理。通過對實驗結(jié)果的分析得到以下的結(jié)論:

    (1) 實驗中利用3DPIV技術(shù)拍攝并進(jìn)行流場處理,得到了翼尖渦的物理特征。通過實驗結(jié)果表明,翼尖渦的無量綱環(huán)量以及渦核的尺寸會隨著機(jī)翼迎角以及風(fēng)速的增大而增大。

    (2) 基于高精度單分量傳感器,作者提出并設(shè)計了一套翼尖氣動阻力測量實驗裝置,并根據(jù)該實驗裝置測量得到的氣動阻力進(jìn)行計算求得翼尖渦引起的誘導(dǎo)阻力。

    [1]馬漢東,崔爾杰.大型飛機(jī)阻力預(yù)示與減阻研究[J].力學(xué)與實踐,2007,29(2):1-8.Ma H D,Cui E J.Drag prediction and reduction for civil transportation aircraft[J].Mechanics in Engineering,2007,29(2):1-8.

    [2]Birch D,Lee T.Structure and induced drag of a tip vortex[J].Journal of Aircraft,2004,41(5):1138-1145.

    [3]Devenport W J,Rife M C,Liapis S I,et al.The structure and development of a wing-tip vortex[J].Journal of Fluid Mechanics,1996,312:67-106.

    [4]江永泉.翼梢小翼的空氣動力機(jī)理[J].民用飛機(jī)設(shè)計與研究,1993,(3):16-22.

    [5]顧蘊(yùn)松,程克明,鄭新軍.翼尖渦流場特性及其控制[J].空氣動力學(xué)學(xué)報,2008,26(4):446-451.Gu Y S,Cheng K M,Zheng X J.Flow field characteristics of wing tips vortex and its control[J].Acta Aerodynamica Sinica,2008,26(4):446-451.

    [6]Gatto A,Mattioni F,Friswell M I.Experimental investigation of bistable winglets to enhance wing lift takeoff capability[J].Journal of Aircraft,2009,46(2):647-655.

    [7]Sun R M,Daichin.Experimental investigation on tip vortices and aerodynamics[J].Theoretical and Applied Mechanics Letters,2011,1(3):1-6.

    [8]Ashurst W T,Meiburg E.Three-dimensional shear layers via vortex dynamics[J].J Fluid Mech,1988,189:87-116.

    [9]Knio O M,Ghoniem A F.Numerical study of a 3-dimensional vortex method[J].Journal of Computational Physics,1990,86(1):75-106.

    [10]Cai J.LES for wing tip vortex around an airfoil[D].Texas:University of Texas,2006:74-75.

    [11]Lee C S.Flow structure and scaling laws in lateral wing-tip blowing[J].AIAA J,1989,27:1002-1007.

    [12]Lee C S,Tavella D.Flow structure of lateral wing-tip blowing[R].AIAA-1986-1810,1986.

    [13]Margaris P,Gursul I.Effect of steady blowing on wing tip flowfield[R].AIAA-2004-2619,2004.

    [14]Catalano F M,Ceron-Mufioz H D.Experimental analysis of the aerodynamic characteristics of adaptive multi-winglets[C].43rd AIAA Aerospace Sciences Meeting and Exhibit,Reno:AIAA,2005.

    [15]Falclio L,Gomes A A,Suleman A.Aero-structural design optimization of a morphing wingtip[J].Journal of Intelligent Material Systems and Structures,2011,22(10):1113-1124.

    [16]Boiler C,Kuo C M.Demonstration of adaptive structure performance on modular micro airvehicle[C].51st AIAA/ASME/ASCE/AHS/ASC Structures,Structural Dynamics,and Materials Conference,Orlando:AIAA,2010.

    [17]Kroo I.Drag due to lift concepts for prediction and reduction[J].Annual Review of Fluid Mechanics,2001,33:587-617.

    [18]Chao D D,Dam C P.Wing drag prediction and de-composition[R].AIAA-2004-5074,2004.

    [19]Ringuette M J,Milano M,Gharib M.Role of the tip vortex in the force generation of low-aspect-ratio normal flat plates[J].Journal of Fluid Mechanics,2007,581:453-468.

    [20]Dacles-Mariani J,Kwak D,Zilliac G G.On numerical errors and turbulence modeling in tip vortex flow prediction[J].International Journal for Numerical Methods in Fluids,1999,30:65-82.

    [21]Takahashi R K,Mcalister P K W.Study of a wing-tip vortex using laser velocimetry[R].NASA-Technical Memorandum-88343,Washington,D.C.:NASA,1987.

    [22]Li J,Cai J,Liu C.Large eddy simulation of wing tip vortex in the near field[R].Mathematics Preprint Series.The University of Texas at Arlington,2007.

    [23]Bailey S C C,Tavoularis S.Measurements of the velocity field of a wing-tip vortex,wandering in grid turbulence[J].J Fluid Mech,2008,601:281-315.

    [24]Prasad A K,Jensen K.Scheimpflug stereo camera for particle image velocimetry in liquid flows[J].Appl Opt,1995,34:7092.

    [25]Jardin T,Farcy A,David L.Three-dimensional effects in hovering flapping flight[J].J Fluid Mech,2012,702:102-125.

    黃文濤(1985-),男,山東臨沂人,博士研究生。研究方向:非定??諝鈩恿W(xué)。通信地址:上海市閔行區(qū)東川路800號上海交通大學(xué)航空航天學(xué)院2406室(200240)。E-mail:herosleeper2000@163.com

    (編輯:張巧蕓)

    Experimentalstudyondrag-reductionmechanismsbasedonthephysicalcharacteristicoftipvortex

    Huang Wentao,Xiang Yang,Wang Xiao,Liu Hong*,Gu Dingyi

    (School of Aeronautics and Astronautics,Shanghai Jiao Tong University,Shanghai 200240,China)

    The drag-reduction mechanisms based on the physical characteristic of the tip vortex are investigated through experiments in the wind tunnel.Through 3DPIV (3D Particle Image Velocimetry) experiments,the physical characteristic of tip vortex is obtained,and the induced drag of wing is calculated based on the aerodynamic force measurement setup of wingtip proposed in this paper.Experimental results show that the non-dimensional circulation of the tip vortex increases with the increasing angle of attack and the wind speed.Meanwhile,with the decrease of the non-dimension circulation of the tip vortex or the increase of the spacing between the wing and the tip vortex,the induced drag becomes smaller and smaller.Specifically,the induced drag reduction can be achieved by inhibiting the non-dimensional intensity of the tip vortex,which weakens the interaction between the main wing and the tip vortex.

    tip vortex;physical characteristic;induced drag;circulation;drag reduction

    O355

    A

    1672-9897(2017)05-0053-07

    10.11729/syltlx20160194

    2016-12-11;

    2017-08-23

    973計劃項目“大型客機(jī)減阻機(jī)理和方法研究”(2014CB744802)

    *通信作者 E-mail:hongliu@sjtu.edu.cn

    HuangWT,XiangY,WangX,etal.Experimentalstudyondrag-reductionmechanismsbasedonthephysicalcharacteristicoftipvortex.JournalofExperimentsinFluidMechanics,2017,31(5):53-59.黃文濤,向 陽,王 笑,等.基于翼尖渦物理特征的誘導(dǎo)阻力減阻機(jī)制實驗研究.實驗流體力學(xué),2017,31(5):53-59.

    猜你喜歡
    環(huán)量渦的翼尖
    中高速條件下不同翼尖小翼的數(shù)值模擬分析
    南海中尺度渦的形轉(zhuǎn)、內(nèi)轉(zhuǎn)及平移運動研究
    海洋通報(2020年5期)2021-01-14 09:26:52
    葉輪出口環(huán)量非線性分布條件下混流泵性能研究
    等-變環(huán)量設(shè)計葉片軸流風(fēng)機(jī)性能研究
    基于模式函數(shù)和變分法的螺旋槳最佳環(huán)量計算方法
    翼尖尖點散射特性分析*
    溝槽對湍流邊界層中展向渦影響的實驗研究
    開縫圓柱縫隙傾斜角對脫落渦的影響
    基于流動顯示的翼尖渦不穩(wěn)定頻率測量
    基于 FFD 技術(shù)的民用運輸機(jī)翼尖裝置設(shè)計研究
    中文字幕熟女人妻在线| 国产爱豆传媒在线观看| 久久久a久久爽久久v久久| 国产片特级美女逼逼视频| 91久久精品国产一区二区三区| 精品人妻偷拍中文字幕| 亚洲,欧美,日韩| av在线老鸭窝| 干丝袜人妻中文字幕| 亚洲精品影视一区二区三区av| 97超视频在线观看视频| 寂寞人妻少妇视频99o| 亚洲中文日韩欧美视频| 日本黄色片子视频| 久久精品影院6| 亚洲国产精品sss在线观看| 白带黄色成豆腐渣| 亚洲国产精品国产精品| 给我免费播放毛片高清在线观看| 久久精品综合一区二区三区| 99久久精品热视频| 亚洲中文字幕一区二区三区有码在线看| 日韩成人伦理影院| 老师上课跳d突然被开到最大视频| 午夜激情福利司机影院| 色播亚洲综合网| 国产一区二区在线av高清观看| 亚洲成人精品中文字幕电影| 成年女人看的毛片在线观看| 伦精品一区二区三区| av黄色大香蕉| 在线看三级毛片| 欧美xxxx黑人xx丫x性爽| 欧美性猛交╳xxx乱大交人| 久久久久久久久大av| 国产高清视频在线播放一区| 国产在线男女| 亚洲欧美精品综合久久99| 插阴视频在线观看视频| 国产国拍精品亚洲av在线观看| 我的老师免费观看完整版| 少妇高潮的动态图| 国产精品人妻久久久久久| 国产v大片淫在线免费观看| 久久韩国三级中文字幕| 国产精品精品国产色婷婷| 亚洲国产高清在线一区二区三| 亚洲在线观看片| 亚洲欧美精品综合久久99| 在现免费观看毛片| 乱人视频在线观看| 国产成人91sexporn| 麻豆久久精品国产亚洲av| 久久久久免费精品人妻一区二区| 老司机福利观看| 中国美女看黄片| 91久久精品国产一区二区三区| 看非洲黑人一级黄片| 国产极品精品免费视频能看的| 国产在线男女| 中文字幕人妻熟人妻熟丝袜美| 国产人妻一区二区三区在| 在线观看美女被高潮喷水网站| 日本 av在线| 精品人妻偷拍中文字幕| 日韩,欧美,国产一区二区三区 | 久久精品人妻少妇| 国产探花极品一区二区| 国产精品永久免费网站| 国产亚洲欧美98| 国产黄色视频一区二区在线观看 | 久久久久国产精品人妻aⅴ院| 免费av观看视频| 波多野结衣巨乳人妻| 99久久九九国产精品国产免费| 麻豆一二三区av精品| 日韩欧美国产在线观看| 久久久成人免费电影| 国产真实伦视频高清在线观看| 99国产极品粉嫩在线观看| 国产女主播在线喷水免费视频网站 | 亚洲综合色惰| 在线观看一区二区三区| 最新中文字幕久久久久| 色噜噜av男人的天堂激情| av在线观看视频网站免费| 久久欧美精品欧美久久欧美| 久久精品人妻少妇| 一进一出抽搐gif免费好疼| 人人妻人人澡人人爽人人夜夜 | 给我免费播放毛片高清在线观看| av视频在线观看入口| 国产高清有码在线观看视频| 免费观看人在逋| 深爱激情五月婷婷| 一级毛片久久久久久久久女| 97在线视频观看| 久久精品影院6| 丰满人妻一区二区三区视频av| 搡老岳熟女国产| 国产av在哪里看| 不卡一级毛片| 欧美日韩一区二区视频在线观看视频在线 | 一本精品99久久精品77| 久久久国产成人精品二区| 少妇熟女欧美另类| 成人av在线播放网站| 欧美成人一区二区免费高清观看| 日本免费a在线| 黄色视频,在线免费观看| 国产精品免费一区二区三区在线| 丰满的人妻完整版| 香蕉av资源在线| 日日撸夜夜添| 国产精品久久久久久精品电影| 最新中文字幕久久久久| 国产精品无大码| 日本与韩国留学比较| 欧美三级亚洲精品| 九九在线视频观看精品| 久久久久国产网址| 91在线观看av| 大又大粗又爽又黄少妇毛片口| 十八禁国产超污无遮挡网站| 免费黄网站久久成人精品| 成年女人毛片免费观看观看9| 在线看三级毛片| 久久久久九九精品影院| 99国产精品一区二区蜜桃av| 久久午夜福利片| 亚洲五月天丁香| 午夜老司机福利剧场| 俺也久久电影网| 亚洲成人久久爱视频| 夜夜夜夜夜久久久久| 欧美成人精品欧美一级黄| 校园春色视频在线观看| 免费一级毛片在线播放高清视频| 久久鲁丝午夜福利片| 亚洲专区国产一区二区| 午夜福利高清视频| 亚洲精品一区av在线观看| 亚洲自拍偷在线| 国产午夜精品久久久久久一区二区三区 | 日本成人三级电影网站| 亚洲人成网站在线播放欧美日韩| 三级国产精品欧美在线观看| 国产高潮美女av| 91午夜精品亚洲一区二区三区| av天堂在线播放| 亚洲成人精品中文字幕电影| 久久精品夜夜夜夜夜久久蜜豆| 精品乱码久久久久久99久播| АⅤ资源中文在线天堂| 久久久久久国产a免费观看| 中文字幕熟女人妻在线| 天美传媒精品一区二区| 18禁在线无遮挡免费观看视频 | 国产精品99久久久久久久久| 国产伦一二天堂av在线观看| 国产午夜精品论理片| 免费看光身美女| 免费在线观看影片大全网站| 悠悠久久av| 日韩av不卡免费在线播放| 久久精品人妻少妇| 美女大奶头视频| 国产伦精品一区二区三区视频9| 97超视频在线观看视频| 亚洲性夜色夜夜综合| 在现免费观看毛片| 女人被狂操c到高潮| 又爽又黄a免费视频| 精品久久久久久久人妻蜜臀av| 伊人久久精品亚洲午夜| 国产精品久久久久久久电影| 国产女主播在线喷水免费视频网站 | 国产探花在线观看一区二区| 日韩欧美在线乱码| 看免费成人av毛片| 三级国产精品欧美在线观看| 一a级毛片在线观看| 男女边吃奶边做爰视频| 乱人视频在线观看| 少妇猛男粗大的猛烈进出视频 | 国产伦精品一区二区三区四那| 国产午夜精品论理片| 精品福利观看| 久久精品91蜜桃| 男人舔奶头视频| 亚洲欧美成人综合另类久久久 | 校园春色视频在线观看| 精品久久久噜噜| 欧美人与善性xxx| av.在线天堂| 国产精品福利在线免费观看| 免费电影在线观看免费观看| 国产国拍精品亚洲av在线观看| 成人漫画全彩无遮挡| 午夜福利在线观看吧| 51国产日韩欧美| 亚洲天堂国产精品一区在线| 在线国产一区二区在线| 在线a可以看的网站| 精品人妻偷拍中文字幕| 日韩中字成人| 久久人人爽人人爽人人片va| 观看美女的网站| 久久久色成人| 久久草成人影院| 老师上课跳d突然被开到最大视频| 亚洲av中文字字幕乱码综合| 黄色配什么色好看| ponron亚洲| 18禁黄网站禁片免费观看直播| www.色视频.com| 深夜精品福利| 超碰av人人做人人爽久久| 嫩草影院入口| 美女高潮的动态| 最近最新中文字幕大全电影3| 亚洲中文日韩欧美视频| 成人亚洲欧美一区二区av| 22中文网久久字幕| 久久婷婷人人爽人人干人人爱| a级毛片a级免费在线| 日韩av在线大香蕉| 中文在线观看免费www的网站| av在线老鸭窝| 色5月婷婷丁香| 日本黄色视频三级网站网址| 国产一区二区三区av在线 | 亚洲美女视频黄频| 中文字幕熟女人妻在线| 最后的刺客免费高清国语| 一级av片app| 日韩大尺度精品在线看网址| 久久久久国产网址| 日韩av不卡免费在线播放| 国产蜜桃级精品一区二区三区| 亚洲精品国产成人久久av| 亚洲欧美日韩卡通动漫| 日本欧美国产在线视频| 亚洲av不卡在线观看| 丝袜喷水一区| 日日摸夜夜添夜夜添小说| 一级黄片播放器| 日韩高清综合在线| 色哟哟·www| 九九久久精品国产亚洲av麻豆| 国产精品亚洲一级av第二区| 久久久久久九九精品二区国产| 蜜臀久久99精品久久宅男| 国产高清激情床上av| 黄色日韩在线| 国产成人freesex在线 | 两个人视频免费观看高清| 国产真实乱freesex| 九九爱精品视频在线观看| 日韩成人伦理影院| 三级国产精品欧美在线观看| 丰满的人妻完整版| 性色avwww在线观看| 国产黄色视频一区二区在线观看 | 国产高清视频在线播放一区| eeuss影院久久| 亚洲国产欧美人成| 在线观看午夜福利视频| 欧美高清性xxxxhd video| 18+在线观看网站| 国产精品久久久久久久电影| 免费看a级黄色片| 国产欧美日韩精品亚洲av| 在线观看午夜福利视频| 免费看日本二区| av在线老鸭窝| 伦理电影大哥的女人| 日日啪夜夜撸| 精品久久久久久久久av| 69av精品久久久久久| a级毛色黄片| 真人做人爱边吃奶动态| 夜夜看夜夜爽夜夜摸| www日本黄色视频网| 亚洲精品国产成人久久av| 性色avwww在线观看| 欧洲精品卡2卡3卡4卡5卡区| 日韩国内少妇激情av| 亚洲一区二区三区色噜噜| 国产精品一及| 卡戴珊不雅视频在线播放| 秋霞在线观看毛片| 大又大粗又爽又黄少妇毛片口| 黄色日韩在线| 久久人人精品亚洲av| 丰满乱子伦码专区| 麻豆国产av国片精品| 男女那种视频在线观看| 精品福利观看| 国产 一区精品| 亚洲无线在线观看| 欧美色视频一区免费| 久久久久精品国产欧美久久久| 国产精品一区二区三区四区免费观看 | 欧美又色又爽又黄视频| 十八禁网站免费在线| 亚洲va在线va天堂va国产| 级片在线观看| 亚洲av二区三区四区| avwww免费| 欧美激情在线99| 欧美中文日本在线观看视频| 看黄色毛片网站| 成人毛片a级毛片在线播放| 国产精品人妻久久久影院| 国产亚洲精品久久久com| 老女人水多毛片| 你懂的网址亚洲精品在线观看 | 色av中文字幕| 内射极品少妇av片p| 村上凉子中文字幕在线| 中文字幕熟女人妻在线| 成人av在线播放网站| 亚洲人成网站高清观看| 白带黄色成豆腐渣| 亚洲成人精品中文字幕电影| 露出奶头的视频| 亚洲欧美日韩高清专用| 国产乱人偷精品视频| 日韩精品中文字幕看吧| av天堂在线播放| or卡值多少钱| 久久国内精品自在自线图片| 亚洲图色成人| 国产白丝娇喘喷水9色精品| 国产一区二区在线观看日韩| 久久6这里有精品| 如何舔出高潮| 一区福利在线观看| 最近手机中文字幕大全| 日本成人三级电影网站| 国产成人a∨麻豆精品| 欧美激情久久久久久爽电影| 日韩av在线大香蕉| 亚洲经典国产精华液单| 中国美女看黄片| 一个人看的www免费观看视频| 亚洲欧美日韩卡通动漫| 精品福利观看| 中国美女看黄片| 国语自产精品视频在线第100页| 亚洲av熟女| 搞女人的毛片| 亚洲欧美日韩卡通动漫| 网址你懂的国产日韩在线| 69av精品久久久久久| 成人高潮视频无遮挡免费网站| 欧美性猛交╳xxx乱大交人| 国产亚洲精品综合一区在线观看| 深爱激情五月婷婷| 久久精品国产亚洲av涩爱 | 欧美3d第一页| 亚洲精品成人久久久久久| 国产中年淑女户外野战色| 日韩av不卡免费在线播放| 精品人妻熟女av久视频| 欧美又色又爽又黄视频| 啦啦啦啦在线视频资源| www日本黄色视频网| 精品午夜福利视频在线观看一区| 成人国产麻豆网| 九九爱精品视频在线观看| 一区福利在线观看| 精品人妻一区二区三区麻豆 | 亚洲av不卡在线观看| 婷婷色综合大香蕉| 亚洲人成网站在线播| 成人国产麻豆网| 亚洲不卡免费看| 国内揄拍国产精品人妻在线| 欧美国产日韩亚洲一区| 国产精品免费一区二区三区在线| 中文字幕久久专区| 露出奶头的视频| 97超碰精品成人国产| 成人漫画全彩无遮挡| 欧美+日韩+精品| 嫩草影院新地址| 国产毛片a区久久久久| 少妇猛男粗大的猛烈进出视频 | 日日撸夜夜添| 午夜精品一区二区三区免费看| 久久精品影院6| 一级毛片aaaaaa免费看小| 麻豆国产97在线/欧美| 赤兔流量卡办理| 久久久精品欧美日韩精品| 中国美白少妇内射xxxbb| 偷拍熟女少妇极品色| 国产伦在线观看视频一区| 天美传媒精品一区二区| 国语自产精品视频在线第100页| 久久午夜福利片| 国产精品久久久久久久电影| 卡戴珊不雅视频在线播放| 婷婷精品国产亚洲av| 亚洲精品久久国产高清桃花| 国产精品一区二区性色av| 久99久视频精品免费| 亚洲av中文av极速乱| 久久久久性生活片| 少妇被粗大猛烈的视频| 国产在视频线在精品| 亚洲av美国av| 国产激情偷乱视频一区二区| 国产精品女同一区二区软件| 国内精品一区二区在线观看| 性色avwww在线观看| 国产一级毛片七仙女欲春2| 国产淫片久久久久久久久| 久久精品人妻少妇| 蜜桃亚洲精品一区二区三区| 丰满的人妻完整版| 精品国内亚洲2022精品成人| 国产精品精品国产色婷婷| 亚洲欧美日韩东京热| 丝袜美腿在线中文| 免费观看在线日韩| 午夜激情福利司机影院| 18禁黄网站禁片免费观看直播| 高清毛片免费看| 国产私拍福利视频在线观看| 免费一级毛片在线播放高清视频| 黄色视频,在线免费观看| 人妻制服诱惑在线中文字幕| 久久久久国内视频| 午夜福利高清视频| 亚洲精品456在线播放app| 亚洲欧美成人综合另类久久久 | 中文字幕免费在线视频6| 久久久久久伊人网av| 亚洲欧美日韩高清在线视频| 成人漫画全彩无遮挡| 国产精品国产高清国产av| 久久久久国内视频| 日韩人妻高清精品专区| 亚洲欧美精品综合久久99| 一进一出好大好爽视频| 国产激情偷乱视频一区二区| 嫩草影院新地址| 国产高清不卡午夜福利| 久99久视频精品免费| 亚洲成人精品中文字幕电影| 精品久久久久久久久久久久久| 亚洲熟妇中文字幕五十中出| 97超碰精品成人国产| 亚洲av熟女| 夜夜看夜夜爽夜夜摸| 日韩 亚洲 欧美在线| 欧美日韩在线观看h| 日本三级黄在线观看| 亚洲国产欧洲综合997久久,| 91狼人影院| 欧美成人a在线观看| 亚洲成人精品中文字幕电影| 欧美另类亚洲清纯唯美| 亚洲熟妇中文字幕五十中出| 国产成人一区二区在线| 在线观看午夜福利视频| 91午夜精品亚洲一区二区三区| 亚洲精品456在线播放app| 丰满的人妻完整版| 天堂av国产一区二区熟女人妻| 欧美日韩乱码在线| 中文亚洲av片在线观看爽| 日韩 亚洲 欧美在线| 搡老岳熟女国产| 国产大屁股一区二区在线视频| 亚洲熟妇熟女久久| 国产成人freesex在线 | 最好的美女福利视频网| 久久精品国产鲁丝片午夜精品| 天堂动漫精品| 欧美激情国产日韩精品一区| 俺也久久电影网| 国产三级中文精品| 高清毛片免费看| 少妇的逼好多水| 久久精品国产亚洲av涩爱 | 99热这里只有是精品50| 黄色视频,在线免费观看| 亚洲经典国产精华液单| 可以在线观看的亚洲视频| 美女免费视频网站| 淫妇啪啪啪对白视频| 18禁在线播放成人免费| 插逼视频在线观看| 在线看三级毛片| 国产伦精品一区二区三区视频9| 精品国内亚洲2022精品成人| а√天堂www在线а√下载| 亚洲国产高清在线一区二区三| 久久人人爽人人片av| 热99在线观看视频| 天堂av国产一区二区熟女人妻| 两个人视频免费观看高清| 午夜福利在线观看免费完整高清在 | 一级av片app| 丰满人妻一区二区三区视频av| 欧美高清成人免费视频www| av卡一久久| 国产高潮美女av| 精品午夜福利视频在线观看一区| 国产亚洲91精品色在线| 亚洲最大成人中文| 精品一区二区三区视频在线| 3wmmmm亚洲av在线观看| 成人一区二区视频在线观看| av视频在线观看入口| 色哟哟哟哟哟哟| 亚洲精品456在线播放app| 黑人高潮一二区| 欧美+日韩+精品| 午夜福利视频1000在线观看| 亚洲成人精品中文字幕电影| 99热这里只有是精品50| 精品免费久久久久久久清纯| 身体一侧抽搐| 人妻久久中文字幕网| 国产国拍精品亚洲av在线观看| 日韩,欧美,国产一区二区三区 | 熟女电影av网| 桃色一区二区三区在线观看| 在线免费观看的www视频| 中国美白少妇内射xxxbb| 一级毛片我不卡| 国产男靠女视频免费网站| 国产大屁股一区二区在线视频| 听说在线观看完整版免费高清| 欧美区成人在线视频| 99热全是精品| 美女黄网站色视频| 搞女人的毛片| 麻豆成人午夜福利视频| 午夜a级毛片| 真人做人爱边吃奶动态| 午夜福利在线观看吧| av天堂在线播放| 观看免费一级毛片| 少妇裸体淫交视频免费看高清| 日本-黄色视频高清免费观看| 亚洲精品456在线播放app| 香蕉av资源在线| 色5月婷婷丁香| 国内精品宾馆在线| 国产一区二区在线观看日韩| 天堂av国产一区二区熟女人妻| 18禁在线无遮挡免费观看视频 | 色尼玛亚洲综合影院| 亚洲精品国产成人久久av| 女同久久另类99精品国产91| 啦啦啦啦在线视频资源| 亚洲av成人av| 黄色一级大片看看| 51国产日韩欧美| 一进一出好大好爽视频| 国产精品亚洲美女久久久| 久久久色成人| 3wmmmm亚洲av在线观看| 麻豆一二三区av精品| 久久久久国内视频| 国产真实伦视频高清在线观看| 淫妇啪啪啪对白视频| 欧美zozozo另类| 91午夜精品亚洲一区二区三区| 天堂影院成人在线观看| 日韩成人av中文字幕在线观看 | 晚上一个人看的免费电影| 亚洲熟妇熟女久久| av福利片在线观看| 干丝袜人妻中文字幕| 免费高清视频大片| 亚洲经典国产精华液单| 亚洲成a人片在线一区二区| 精品午夜福利在线看| 成年版毛片免费区| 99热全是精品| 国产精品一区二区三区四区久久| 大香蕉久久网| 国产日本99.免费观看| 亚洲国产欧洲综合997久久,| 精品福利观看| 夜夜夜夜夜久久久久| 亚洲av成人av| 欧美最新免费一区二区三区| 欧美3d第一页| 男女边吃奶边做爰视频| 99视频精品全部免费 在线| 国产精品久久久久久久久免| 色综合站精品国产| 欧美日韩一区二区视频在线观看视频在线 | 人妻少妇偷人精品九色| 国产精品久久久久久精品电影| 国产免费男女视频| 国产美女午夜福利| 亚洲一级一片aⅴ在线观看| 欧美激情久久久久久爽电影| 韩国av在线不卡| 亚洲av成人av| 在线a可以看的网站| 蜜桃亚洲精品一区二区三区| 久久久精品94久久精品| 日本-黄色视频高清免费观看| 免费黄网站久久成人精品| 免费看美女性在线毛片视频|