• <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è)計研究
    日韩欧美 国产精品| 国国产精品蜜臀av免费| 日韩免费高清中文字幕av| 嘟嘟电影网在线观看| 欧美变态另类bdsm刘玥| 午夜福利网站1000一区二区三区| 日韩中字成人| 日本爱情动作片www.在线观看| 国产中年淑女户外野战色| 国产色婷婷99| 成人毛片a级毛片在线播放| 久久久久久久久大av| 久久人人爽av亚洲精品天堂 | 午夜免费观看性视频| 亚洲电影在线观看av| 精品国产三级普通话版| 精品少妇久久久久久888优播| 亚洲人与动物交配视频| 欧美高清成人免费视频www| 国产免费一级a男人的天堂| 亚洲av男天堂| 午夜福利视频1000在线观看| 成人美女网站在线观看视频| 国产女主播在线喷水免费视频网站| av免费在线看不卡| 毛片女人毛片| 男女无遮挡免费网站观看| 国产精品一区二区三区四区免费观看| 国产精品av视频在线免费观看| 免费观看性生交大片5| 99热这里只有是精品在线观看| 国产精品蜜桃在线观看| 内射极品少妇av片p| 国产极品天堂在线| 国产精品女同一区二区软件| 久久久久性生活片| 精品久久久久久久人妻蜜臀av| 两个人的视频大全免费| 欧美精品国产亚洲| 最近手机中文字幕大全| 国产大屁股一区二区在线视频| 亚洲经典国产精华液单| 亚洲精品成人久久久久久| 一级毛片久久久久久久久女| 免费看av在线观看网站| 特级一级黄色大片| freevideosex欧美| 国国产精品蜜臀av免费| 婷婷色av中文字幕| 女人久久www免费人成看片| 亚洲熟女精品中文字幕| 久久韩国三级中文字幕| 国产亚洲av片在线观看秒播厂| 国产高清不卡午夜福利| 亚洲欧洲国产日韩| 在线观看免费高清a一片| 亚洲精品456在线播放app| 国产成人一区二区在线| 国产午夜精品久久久久久一区二区三区| av在线天堂中文字幕| 97精品久久久久久久久久精品| 可以在线观看毛片的网站| 国产熟女欧美一区二区| av又黄又爽大尺度在线免费看| 亚洲精品国产av成人精品| 日日啪夜夜爽| 日本黄大片高清| 国产成人精品婷婷| 国产精品一区www在线观看| 免费人成在线观看视频色| 麻豆久久精品国产亚洲av| av又黄又爽大尺度在线免费看| 久久久色成人| 国产色婷婷99| 久久99蜜桃精品久久| 超碰av人人做人人爽久久| 久久精品久久久久久久性| 亚洲无线观看免费| 视频区图区小说| 欧美xxxx黑人xx丫x性爽| 97在线人人人人妻| 大话2 男鬼变身卡| 我要看日韩黄色一级片| 亚洲精品国产av蜜桃| 国产欧美亚洲国产| 一区二区av电影网| 国产精品嫩草影院av在线观看| 欧美一级a爱片免费观看看| 91久久精品电影网| 有码 亚洲区| 熟女人妻精品中文字幕| 亚洲国产日韩一区二区| 国产成年人精品一区二区| 在现免费观看毛片| 欧美老熟妇乱子伦牲交| 成人国产麻豆网| 性插视频无遮挡在线免费观看| 亚洲图色成人| 国产女主播在线喷水免费视频网站| 亚洲内射少妇av| 亚洲成人久久爱视频| 人妻 亚洲 视频| 亚洲精品自拍成人| 欧美激情在线99| 亚洲美女搞黄在线观看| 国产精品一区二区性色av| 国国产精品蜜臀av免费| 亚洲国产成人一精品久久久| 三级男女做爰猛烈吃奶摸视频| 我的女老师完整版在线观看| 久久久久国产精品人妻一区二区| 日韩,欧美,国产一区二区三区| 亚洲精品成人久久久久久| 好男人视频免费观看在线| 涩涩av久久男人的天堂| 精品人妻偷拍中文字幕| 啦啦啦中文免费视频观看日本| 欧美极品一区二区三区四区| 三级经典国产精品| 精品人妻视频免费看| 欧美激情国产日韩精品一区| 日日啪夜夜爽| 久久久a久久爽久久v久久| av在线app专区| 日本熟妇午夜| 人人妻人人看人人澡| 人妻制服诱惑在线中文字幕| 国产男女超爽视频在线观看| 亚洲av日韩在线播放| 国产av不卡久久| 欧美xxxx性猛交bbbb| 熟女av电影| 精华霜和精华液先用哪个| 卡戴珊不雅视频在线播放| 国产日韩欧美亚洲二区| 亚洲av国产av综合av卡| 精品久久久久久久末码| 午夜视频国产福利| 久久99蜜桃精品久久| 亚洲aⅴ乱码一区二区在线播放| 99久久人妻综合| 熟女电影av网| 综合色丁香网| 成人美女网站在线观看视频| 五月伊人婷婷丁香| 亚洲最大成人手机在线| 日本熟妇午夜| 亚洲欧洲日产国产| 亚洲久久久久久中文字幕| 插逼视频在线观看| 欧美性感艳星| 国产高清三级在线| 深夜a级毛片| 最后的刺客免费高清国语| 老师上课跳d突然被开到最大视频| 欧美精品一区二区大全| 一二三四中文在线观看免费高清| 免费人成在线观看视频色| 久久99蜜桃精品久久| 亚洲婷婷狠狠爱综合网| 小蜜桃在线观看免费完整版高清| 欧美日韩在线观看h| 嫩草影院入口| 80岁老熟妇乱子伦牲交| 免费播放大片免费观看视频在线观看| 亚洲精品日本国产第一区| 久久久欧美国产精品| 哪个播放器可以免费观看大片| 色播亚洲综合网| 99热这里只有是精品50| 大陆偷拍与自拍| 高清日韩中文字幕在线| 人妻少妇偷人精品九色| 岛国毛片在线播放| 97超视频在线观看视频| 如何舔出高潮| 亚洲av中文av极速乱| 亚洲,一卡二卡三卡| 天堂网av新在线| 中文欧美无线码| 卡戴珊不雅视频在线播放| 国产色婷婷99| 69人妻影院| 一区二区三区乱码不卡18| 男女下面进入的视频免费午夜| 国产成人a∨麻豆精品| 午夜爱爱视频在线播放| 亚洲精品乱码久久久v下载方式| 伦精品一区二区三区| 午夜福利网站1000一区二区三区| 黄片无遮挡物在线观看| 日韩av不卡免费在线播放| 中文精品一卡2卡3卡4更新| 啦啦啦中文免费视频观看日本| 国产中年淑女户外野战色| 久久久久久久久久久丰满| 十八禁网站网址无遮挡 | 精品少妇久久久久久888优播| 麻豆精品久久久久久蜜桃| 97超碰精品成人国产| 国产乱人视频| 国产av不卡久久| 色婷婷久久久亚洲欧美| 一级爰片在线观看| 国产精品久久久久久久电影| 美女国产视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩亚洲高清精品| 日韩,欧美,国产一区二区三区| 成人黄色视频免费在线看| 欧美日韩视频高清一区二区三区二| 久久精品国产亚洲网站| 婷婷色麻豆天堂久久| 黄片wwwwww| 五月天丁香电影| 国产毛片在线视频| 大片电影免费在线观看免费| av在线app专区| 国产亚洲5aaaaa淫片| 精品久久久久久久末码| 欧美+日韩+精品| 五月玫瑰六月丁香| 高清毛片免费看| 欧美精品一区二区大全| 精品久久久精品久久久| 国产成人精品福利久久| 我要看日韩黄色一级片| 秋霞伦理黄片| 国产精品无大码| 国产免费一区二区三区四区乱码| 国产精品不卡视频一区二区| 亚州av有码| 婷婷色av中文字幕| 久久热精品热| 欧美日韩国产mv在线观看视频 | 国产永久视频网站| 美女xxoo啪啪120秒动态图| 国产日韩欧美亚洲二区| 日韩中字成人| 黄色视频在线播放观看不卡| 日本免费在线观看一区| 丝袜美腿在线中文| 可以在线观看毛片的网站| 少妇人妻 视频| 欧美xxxx性猛交bbbb| 99视频精品全部免费 在线| av在线天堂中文字幕| 国产成人a区在线观看| 日本色播在线视频| 一级黄片播放器| 九九在线视频观看精品| 天天躁日日操中文字幕| 一个人看视频在线观看www免费| 少妇的逼好多水| 精品99又大又爽又粗少妇毛片| 久久精品人妻少妇| 国产精品久久久久久精品电影小说 | 免费黄网站久久成人精品| 精品人妻偷拍中文字幕| 成年人午夜在线观看视频| 国产黄片美女视频| 日本av手机在线免费观看| 夜夜爽夜夜爽视频| 亚洲av男天堂| av线在线观看网站| av国产免费在线观看| 简卡轻食公司| 成人毛片60女人毛片免费| 国产成人免费无遮挡视频| 午夜免费男女啪啪视频观看| 成人国产麻豆网| 精品久久久久久久末码| 中文字幕久久专区| 日产精品乱码卡一卡2卡三| 我的女老师完整版在线观看| 美女cb高潮喷水在线观看| 午夜免费观看性视频| 日韩不卡一区二区三区视频在线| 嘟嘟电影网在线观看| 日韩一本色道免费dvd| 欧美成人a在线观看| 中文资源天堂在线| 久久精品综合一区二区三区| 亚洲丝袜综合中文字幕| 大码成人一级视频| 国产大屁股一区二区在线视频| 3wmmmm亚洲av在线观看| 热re99久久精品国产66热6| 国产免费一级a男人的天堂| 国模一区二区三区四区视频| 免费av不卡在线播放| 一级毛片久久久久久久久女| 亚洲国产精品999| 国产男女内射视频| 国产精品麻豆人妻色哟哟久久| 免费播放大片免费观看视频在线观看| 特大巨黑吊av在线直播| 亚洲欧美日韩无卡精品| 99热这里只有是精品在线观看| 欧美精品一区二区大全| 亚洲国产色片| 寂寞人妻少妇视频99o| 久久午夜福利片| 国产精品99久久久久久久久| 国产精品精品国产色婷婷| 国产探花在线观看一区二区| 啦啦啦在线观看免费高清www| 亚洲图色成人| 亚洲国产av新网站| 蜜桃久久精品国产亚洲av| 亚洲熟女精品中文字幕| 少妇熟女欧美另类| 国产精品国产三级专区第一集| 日韩亚洲欧美综合| av.在线天堂| 成人无遮挡网站| 国产欧美亚洲国产| 69av精品久久久久久| 精品国产露脸久久av麻豆| 欧美性猛交╳xxx乱大交人| 夜夜爽夜夜爽视频| 精品午夜福利在线看| 久久精品国产亚洲av天美| 哪个播放器可以免费观看大片| 免费黄色在线免费观看| 成人毛片a级毛片在线播放| 春色校园在线视频观看| 伦理电影大哥的女人| 观看美女的网站| 国产伦理片在线播放av一区| 久久久精品欧美日韩精品| 男女无遮挡免费网站观看| 黑人高潮一二区| 国产在线一区二区三区精| 精品久久久久久久久亚洲| 成人鲁丝片一二三区免费| 中文字幕av成人在线电影| 国产淫片久久久久久久久| 天堂中文最新版在线下载 | 熟女av电影| 99热这里只有是精品50| 精品午夜福利在线看| av在线老鸭窝| av在线观看视频网站免费| 国产男人的电影天堂91| 人妻夜夜爽99麻豆av| 国产国拍精品亚洲av在线观看| 亚洲精品456在线播放app| 永久网站在线| 有码 亚洲区| 婷婷色av中文字幕| 久久久久久久久久人人人人人人| 亚洲国产欧美在线一区| 亚洲精品456在线播放app| 亚洲av中文av极速乱| 久久久久网色| 日韩国内少妇激情av| 搡老乐熟女国产| 在线观看一区二区三区激情| 久久99精品国语久久久| 在线观看免费高清a一片| 久久久精品免费免费高清| 美女主播在线视频| 国内精品宾馆在线| 亚洲性久久影院| 我的女老师完整版在线观看| 日韩一区二区三区影片| 欧美区成人在线视频| 一级二级三级毛片免费看| 精品亚洲乱码少妇综合久久| 免费高清在线观看视频在线观看| 国产精品久久久久久av不卡| 秋霞伦理黄片| 国产伦精品一区二区三区视频9| 狂野欧美白嫩少妇大欣赏| 老司机影院毛片| 激情 狠狠 欧美| 老女人水多毛片| 婷婷色综合大香蕉| 简卡轻食公司| 99精国产麻豆久久婷婷| 禁无遮挡网站| 麻豆久久精品国产亚洲av| 男人爽女人下面视频在线观看| 边亲边吃奶的免费视频| 亚洲人与动物交配视频| 国产精品成人在线| 久热久热在线精品观看| 又大又黄又爽视频免费| 直男gayav资源| 国产成人免费观看mmmm| 久久久久久久亚洲中文字幕| 国产亚洲91精品色在线| 少妇的逼水好多| 丰满乱子伦码专区| 日韩成人伦理影院| 精品少妇久久久久久888优播| 久久鲁丝午夜福利片| 国产免费视频播放在线视频| 熟女电影av网| 简卡轻食公司| 一区二区三区精品91| 超碰97精品在线观看| 国产精品无大码| 久久久久久久大尺度免费视频| 又大又黄又爽视频免费| 真实男女啪啪啪动态图| 亚洲欧洲日产国产| 人人妻人人爽人人添夜夜欢视频 | 国产高清不卡午夜福利| 久久久色成人| 神马国产精品三级电影在线观看| 成年免费大片在线观看| 精品久久久久久久久亚洲| 亚洲不卡免费看| 少妇人妻久久综合中文| 精品一区二区三卡| 我要看日韩黄色一级片| 岛国毛片在线播放| 一区二区三区免费毛片| av.在线天堂| 简卡轻食公司| 黄片无遮挡物在线观看| 亚洲不卡免费看| 国产精品国产三级国产av玫瑰| 夜夜爽夜夜爽视频| 舔av片在线| 日本熟妇午夜| 一级毛片我不卡| av.在线天堂| 人妻少妇偷人精品九色| 欧美日韩在线观看h| 国产91av在线免费观看| 99久久中文字幕三级久久日本| 成人综合一区亚洲| 91久久精品国产一区二区三区| 日韩成人伦理影院| 男男h啪啪无遮挡| 你懂的网址亚洲精品在线观看| 亚洲精品日本国产第一区| 亚洲精品日韩在线中文字幕| 特大巨黑吊av在线直播| 国产精品国产三级国产av玫瑰| 一级毛片aaaaaa免费看小| av一本久久久久| 99re6热这里在线精品视频| 久久韩国三级中文字幕| 在线看a的网站| 亚洲人与动物交配视频| 欧美日韩精品成人综合77777| 国产乱人偷精品视频| 国产成人精品婷婷| 国产毛片在线视频| 精品久久久久久久末码| 成人毛片a级毛片在线播放| 天堂俺去俺来也www色官网| 26uuu在线亚洲综合色| 亚洲精品乱码久久久v下载方式| 免费播放大片免费观看视频在线观看| 国产精品蜜桃在线观看| 欧美三级亚洲精品| 能在线免费看毛片的网站| 久久久久久久久久人人人人人人| 韩国av在线不卡| 舔av片在线| 岛国毛片在线播放| 午夜福利视频1000在线观看| 国产黄色免费在线视频| 美女高潮的动态| 少妇丰满av| 国产一区有黄有色的免费视频| 久久精品综合一区二区三区| 看免费成人av毛片| 国内精品宾馆在线| 国产精品成人在线| av黄色大香蕉| 色哟哟·www| 人妻一区二区av| 亚洲婷婷狠狠爱综合网| 国内少妇人妻偷人精品xxx网站| 99热全是精品| 最新中文字幕久久久久| 久久久久久久亚洲中文字幕| 99九九线精品视频在线观看视频| av天堂中文字幕网| 久久久久精品性色| 欧美日韩精品成人综合77777| 老女人水多毛片| av一本久久久久| 色综合色国产| 国产精品嫩草影院av在线观看| 蜜臀久久99精品久久宅男| 晚上一个人看的免费电影| 97超视频在线观看视频| 中国美白少妇内射xxxbb| 69av精品久久久久久| 国产成人精品久久久久久| 国产视频内射| 国产精品一及| 久久亚洲国产成人精品v| 一二三四中文在线观看免费高清| 人妻夜夜爽99麻豆av| 又爽又黄无遮挡网站| 99久久中文字幕三级久久日本| 看黄色毛片网站| 国产黄片视频在线免费观看| 亚洲欧美日韩东京热| av在线播放精品| 免费看光身美女| 国产av国产精品国产| 看免费成人av毛片| 好男人在线观看高清免费视频| 国产精品无大码| 久久久a久久爽久久v久久| 嫩草影院精品99| 国产白丝娇喘喷水9色精品| av女优亚洲男人天堂| 欧美日韩在线观看h| 最近最新中文字幕大全电影3| 激情 狠狠 欧美| av天堂中文字幕网| 秋霞在线观看毛片| 亚洲怡红院男人天堂| 少妇的逼水好多| 亚洲高清免费不卡视频| 22中文网久久字幕| 蜜桃久久精品国产亚洲av| 精品久久国产蜜桃| 九草在线视频观看| 色婷婷久久久亚洲欧美| 婷婷色综合大香蕉| 成人一区二区视频在线观看| 国产成人免费无遮挡视频| 男的添女的下面高潮视频| 午夜精品国产一区二区电影 | 亚洲成人一二三区av| 亚洲精品乱久久久久久| 99re6热这里在线精品视频| 毛片女人毛片| 精华霜和精华液先用哪个| 日韩视频在线欧美| 亚洲欧美日韩卡通动漫| 国产精品一区二区性色av| 极品教师在线视频| 中文精品一卡2卡3卡4更新| 亚洲最大成人中文| 日韩亚洲欧美综合| 欧美日韩精品成人综合77777| 大又大粗又爽又黄少妇毛片口| 久久97久久精品| 国产爱豆传媒在线观看| 国产久久久一区二区三区| 欧美国产精品一级二级三级 | 精品少妇黑人巨大在线播放| 99久久精品国产国产毛片| 久久国产乱子免费精品| 大又大粗又爽又黄少妇毛片口| 久久6这里有精品| 在线观看av片永久免费下载| 久久精品国产a三级三级三级| 久久ye,这里只有精品| 九九在线视频观看精品| 婷婷色av中文字幕| 亚洲内射少妇av| 成人亚洲精品一区在线观看 | 国产91av在线免费观看| 丰满乱子伦码专区| 色5月婷婷丁香| 色播亚洲综合网| 日本免费在线观看一区| 亚洲精品456在线播放app| 大片电影免费在线观看免费| 久久久久久久久久久免费av| 久久久久性生活片| 国产欧美另类精品又又久久亚洲欧美| 国产精品99久久久久久久久| 精品一区二区三区视频在线| 精华霜和精华液先用哪个| a级一级毛片免费在线观看| 毛片一级片免费看久久久久| 日产精品乱码卡一卡2卡三| 国产精品麻豆人妻色哟哟久久| 久久精品久久精品一区二区三区| 天天躁夜夜躁狠狠久久av| 极品少妇高潮喷水抽搐| 美女xxoo啪啪120秒动态图| 身体一侧抽搐| 亚洲精品乱码久久久v下载方式| 嫩草影院精品99| 大话2 男鬼变身卡| 亚洲不卡免费看| 偷拍熟女少妇极品色| 午夜福利网站1000一区二区三区| 国产久久久一区二区三区| av在线天堂中文字幕| 人妻夜夜爽99麻豆av| 丝瓜视频免费看黄片| 激情 狠狠 欧美| 国内精品美女久久久久久| 在线亚洲精品国产二区图片欧美 | av天堂中文字幕网| 啦啦啦啦在线视频资源| 国产男女内射视频| 大陆偷拍与自拍| 91精品国产九色| 一级毛片电影观看| 最近中文字幕2019免费版| 国产精品久久久久久久久免| 欧美国产精品一级二级三级 | 亚洲av二区三区四区| 日本黄色片子视频| 一级a做视频免费观看| 老司机影院毛片| 免费看日本二区| 国产亚洲最大av|