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

    地面效應(yīng)作用下翼尖渦特性的PIV實(shí)驗(yàn)研究

    2015-04-14 08:42:30曠,代欽,2,*
    關(guān)鍵詞:翼尖翼面渦量

    章 曠,代 欽,2,*

    (1.上海大學(xué)上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所,上海 200072; 2.上海市力學(xué)在能源工程中的應(yīng)用重點(diǎn)實(shí)驗(yàn)室,上海 200072)

    地面效應(yīng)作用下翼尖渦特性的PIV實(shí)驗(yàn)研究

    章 曠1,代 欽1,2,*

    (1.上海大學(xué)上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所,上海 200072; 2.上海市力學(xué)在能源工程中的應(yīng)用重點(diǎn)實(shí)驗(yàn)室,上海 200072)

    完成了NACA23012機(jī)翼地面效應(yīng)條件下翼尖渦結(jié)構(gòu)及升阻力特性實(shí)驗(yàn)。實(shí)驗(yàn)在拖曳水槽中模擬機(jī)翼的飛行狀態(tài),獲得了在多種飛行高度、0°攻角時(shí)機(jī)翼在水平地面和正弦波浪地面附近的升/阻力、翼尖渦流場(chǎng)的變化規(guī)律,對(duì)比分析了水平地面和波浪地面附近翼尖渦速度、渦量分布的區(qū)別及其可能對(duì)機(jī)翼升/阻力造成的影響。實(shí)驗(yàn)結(jié)果表明:即使在正弦波浪地面附近,隨著機(jī)翼逐漸靠近地面,升力逐漸減小至負(fù)升力,翼尖渦的強(qiáng)度亦發(fā)生相應(yīng)變化;尤其是在小間隙比、負(fù)升力情況下,翼尖渦的旋轉(zhuǎn)方向產(chǎn)生了改變;流場(chǎng)結(jié)構(gòu)不僅受機(jī)翼距地面高度影響,也隨著波浪地面與機(jī)翼瞬時(shí)所處位置構(gòu)成的相對(duì)相位關(guān)系的不同而變化,并且渦量沿波浪地面運(yùn)動(dòng)的變化呈現(xiàn)周期性,但變化規(guī)律并不符合正弦周期,主要原因在于波浪地形與下翼面所構(gòu)成的流道形狀及狹窄程度的周期性變化對(duì)翼尖渦流場(chǎng)結(jié)構(gòu)的發(fā)展和演化產(chǎn)生不同程度的抑制。

    地面效應(yīng);翼尖渦;PIV;波浪地面;拖曳水槽;實(shí)驗(yàn)研究

    0 引 言

    機(jī)翼作為地效飛行器的主要升力面,其空氣動(dòng)力學(xué)特性與流動(dòng)結(jié)構(gòu)之間具有密切的相關(guān)性,并對(duì)飛行器性能產(chǎn)生最直接的影響。深入了解機(jī)翼近地面的流動(dòng)結(jié)構(gòu)及其演化規(guī)律,對(duì)于揭示地面效應(yīng)的物理本質(zhì)具有重要意義。眾所周知,當(dāng)機(jī)翼靠近地面或水面飛行時(shí),相比在自由空間中飛行,往往能獲得更大的升阻比。因此,地效飛行器比普通飛行器具有更高的能源利用率。為此已經(jīng)在相關(guān)領(lǐng)域開(kāi)展了許多研究工作,對(duì)于地面效應(yīng)的空氣動(dòng)力學(xué)特性及流動(dòng)基本規(guī)律有了相當(dāng)多的認(rèn)識(shí)。

    李盾[1]等用數(shù)值模擬的方法發(fā)現(xiàn)在間隙比足夠小、攻角足夠大時(shí),氣流受地面和機(jī)翼下翼面的擠壓,很難立即流向外側(cè),這會(huì)使飛行器底部的壓力增高,升力增加,產(chǎn)生正地面效應(yīng)。但是當(dāng)小于某攻角飛行時(shí),升阻比急劇下降,產(chǎn)生負(fù)效應(yīng)。Ahmed[2]等進(jìn)行的NACA4412機(jī)翼動(dòng)地板地效實(shí)驗(yàn)研究也發(fā)現(xiàn)——在低高度小攻角情況下機(jī)翼下翼面產(chǎn)生負(fù)地面效應(yīng),他們認(rèn)為下翼面與地面形成的先收縮后擴(kuò)張的狹窄流道是產(chǎn)生負(fù)地面效應(yīng)的原因。此外還有更多的針對(duì)波浪地面地效規(guī)律的研究報(bào)告。陳新[3]等的數(shù)值計(jì)算結(jié)果表明:機(jī)翼擾動(dòng)在自由表面上的興波幅度非常小,可以忽略不計(jì);當(dāng)考慮波浪水面對(duì)機(jī)翼氣動(dòng)力的影響時(shí),可用“剛性”波浪地面代替“柔性”的波浪水面。Byelinskyy V G[4]等計(jì)算了使用剛性波浪地面替換柔性波浪水面所造成的誤差,其結(jié)果不超過(guò)1%;測(cè)量了機(jī)翼在硬質(zhì)波浪地面附近水動(dòng)力特性隨間隙及攻角變化的規(guī)律,發(fā)現(xiàn)在正弦波浪地面上方運(yùn)動(dòng)的機(jī)翼其受力情況的變化曲線并不符合正弦規(guī)律,其受力極值的出現(xiàn)位置比機(jī)翼前緣到達(dá)波峰的時(shí)刻略超前。Zhang Xin[5]等使用PIV及LDA對(duì)機(jī)翼地面效應(yīng)進(jìn)行實(shí)驗(yàn)研究,研究指出翼尖渦特性和上翼面的流動(dòng)分離會(huì)改變地面效應(yīng)下機(jī)翼的氣動(dòng)力特性,在機(jī)翼裝有端板時(shí)升力增加,翼尖渦誘導(dǎo)的下洗能減弱翼尖附近的流動(dòng)分離。Jung Kwang Hyo[6]等利用煙線裝置對(duì)NACA6409機(jī)翼進(jìn)行了地面效應(yīng)實(shí)驗(yàn)研究,實(shí)驗(yàn)直觀地觀察了翼尖渦的產(chǎn)生和演化過(guò)程,結(jié)果顯示受地面效應(yīng)影響,間隙比越小翼尖渦越被抑制。Hiemcke[7]在風(fēng)洞中測(cè)量在地面效應(yīng)下NACA5312翼型表面的壓力分布,發(fā)現(xiàn)在大攻角、小間隙高度飛行時(shí)上翼面尾緣流動(dòng)分離被抑制。屈秋林[8]等數(shù)值模擬了地效飛行器近波浪表面飛行的分離流場(chǎng),結(jié)果顯示波浪表面飛行的氣動(dòng)力平均值和水平地面飛行氣動(dòng)力平均值隨飛行高度變化的規(guī)律一致,在波浪表面附近飛行,上翼面流動(dòng)分離的范圍和翼剖面壓力分布均呈周期性變化。

    以往的多數(shù)機(jī)翼地面效應(yīng)風(fēng)洞實(shí)驗(yàn)均采用固定地板或動(dòng)地板裝置,難以避免地面邊界層對(duì)流場(chǎng)的干擾,使流場(chǎng)結(jié)果與實(shí)際流動(dòng)情況不符。對(duì)于固定地板實(shí)驗(yàn),機(jī)翼相對(duì)于地面靜止,這使得實(shí)驗(yàn)?zāi)M的運(yùn)動(dòng)關(guān)系并不真實(shí),同時(shí)在機(jī)翼距離地面較低的情況下,機(jī)翼很可能已經(jīng)處于地面產(chǎn)生的邊界層中。動(dòng)地面雖然能模擬真實(shí)的機(jī)翼相對(duì)地面的運(yùn)動(dòng)關(guān)系,但是由于地面存在震動(dòng),長(zhǎng)時(shí)間運(yùn)動(dòng)地板會(huì)發(fā)熱等,均對(duì)地面附近的流場(chǎng)產(chǎn)生影響。為避免上述問(wèn)題,本研究采用水槽拖曳實(shí)驗(yàn)?zāi)M機(jī)翼地面效應(yīng)。實(shí)驗(yàn)以文獻(xiàn)[3-4]的研究結(jié)果為依據(jù),采用文獻(xiàn)[4]的實(shí)驗(yàn)方法,將機(jī)翼模型在水平底面或正弦波形底面上拖曳,以水為介質(zhì)模擬空氣、硬質(zhì)水平地板和波浪狀地板模擬平靜水面和波浪水面,使得實(shí)驗(yàn)所模擬的相對(duì)運(yùn)動(dòng)關(guān)系與實(shí)際飛行狀況相同。同時(shí),由于使用了水模擬空氣介質(zhì),可以在較小的拖曳速度下獲得一個(gè)相對(duì)大的雷諾數(shù)Re=2.0×105。實(shí)驗(yàn)使用高速PIV系統(tǒng)測(cè)量機(jī)翼翼尖渦速度矢量分布,以期從流場(chǎng)結(jié)構(gòu)與翼載荷的相關(guān)性來(lái)探究地面效應(yīng)的機(jī)理。需要指出的是,文獻(xiàn)[4]實(shí)驗(yàn)側(cè)重于對(duì)機(jī)翼在水平和波浪地面附近的升力阻力特性及規(guī)律的研究,尤其是受力的相位與波浪地面相位之間的關(guān)系,但沒(méi)有研究機(jī)翼繞流的流場(chǎng)結(jié)構(gòu),且實(shí)驗(yàn)雷諾數(shù)達(dá)Re=1×106。而本實(shí)驗(yàn)雷諾數(shù)略低,而且與文獻(xiàn)[4]的最大區(qū)別在于更關(guān)注的是機(jī)翼流場(chǎng)結(jié)構(gòu)的測(cè)量及其與機(jī)翼升力阻力之間的關(guān)系。

    1 實(shí)驗(yàn)裝置

    實(shí)驗(yàn)用拖曳水槽尺寸為6 m×0.5 m×0.5 m,側(cè)壁均為平面玻璃。水平地面效應(yīng)實(shí)驗(yàn)即以水槽玻璃底面為實(shí)驗(yàn)地面。波浪地面效應(yīng)實(shí)驗(yàn)采用波高0.02 m、波長(zhǎng)0.16 m的特制橡膠正弦波浪板拼接鋪滿水槽底,橡膠板與水槽底面同寬,表面光潔。拖曳裝置為一根長(zhǎng)6 m的直線拖曳軌道,架設(shè)于水槽一側(cè)上方,拖曳軌道上裝有滑塊?;瑝K由伺服電機(jī)驅(qū)動(dòng)同步帶使之運(yùn)動(dòng),運(yùn)動(dòng)速度及加減速時(shí)間可調(diào),空載最大速度可達(dá)4.0 m/s,啟動(dòng)最大加速度可達(dá)8.0 m/s2。本實(shí)驗(yàn)中機(jī)翼拖曳運(yùn)動(dòng)速度為2.0 m/s,使用鋁制NACA23012機(jī)翼,弦長(zhǎng)c=0.1 m,展長(zhǎng)L=0.2 m,表面光潔,以弦長(zhǎng)為參考長(zhǎng)度的雷諾數(shù)為Re=2.0×105。機(jī)翼升降、攻角調(diào)節(jié)及測(cè)力裝置均安裝于軌道滑塊之上。拖曳裝置如圖1(a)所示,滑塊上方架設(shè)一個(gè)三角懸臂梁伸至水槽中,懸臂梁上垂直安裝有行程0.2 m直線電控滑軌,用于調(diào)節(jié)機(jī)翼距水槽底面的高度,滑軌平臺(tái)上裝有手動(dòng)攻角調(diào)節(jié)機(jī)構(gòu),同時(shí)裝有夾持器,用于夾持測(cè)力天平。機(jī)翼通過(guò)導(dǎo)板與天平連接,導(dǎo)板具有對(duì)稱機(jī)翼剖面,以減小對(duì)水流的擾動(dòng)。機(jī)翼高度的調(diào)節(jié)精度為1 μm,攻角的調(diào)節(jié)精度為1',高度和攻角調(diào)節(jié)裝置均可鎖止。

    圖1 實(shí)驗(yàn)裝置及布局示意圖Fig.1 Schematic diagram of the experimental setup

    拖曳水槽實(shí)驗(yàn)布局如圖1(b)所示。PIV實(shí)驗(yàn)光源為半導(dǎo)體連續(xù)激光器,輸出功率3 W,波長(zhǎng)532 nm。激光經(jīng)柱透鏡擴(kuò)束后,形成的片光分別從水槽外側(cè)(對(duì)波浪地形)或底部(對(duì)水平地形)照明流體內(nèi)豎直截面,截面與機(jī)翼運(yùn)動(dòng)方向垂直。片光截面位于實(shí)驗(yàn)段中部距離機(jī)翼起始位置約2 m處,位置固定不變。拍攝波浪地面不同相位流場(chǎng)時(shí),僅前后整體移動(dòng)波浪地板,而不調(diào)節(jié)光路。PIV高速相機(jī)布置于水槽上游外0.5 m處,采用105 mm定焦鏡頭,拍攝方向與機(jī)翼前進(jìn)方向相同。相機(jī)空間分辨率1024×1024,全畫(huà)幅幀頻為3000 fps。本實(shí)驗(yàn)中畫(huà)幅寬度與片光照明區(qū)域同寬,采集頻率設(shè)為500幀/s。示蹤粒子采用跟隨性極好且密度與水相近的尼龍微珠,對(duì)激光有良好的散射性,拍攝圖像清晰,信噪比高,且在靜止水中可懸浮三個(gè)小時(shí)以上保持不沉降或上浮。每組拖曳實(shí)驗(yàn)開(kāi)始前在光面位置注入示蹤粒子,待擴(kuò)散均勻后,保持水體靜止一小時(shí),以保證流場(chǎng)湍流度足夠低且每次實(shí)驗(yàn)的湍流度相近。PIV實(shí)驗(yàn)與測(cè)力實(shí)驗(yàn)可同時(shí)進(jìn)行。定義間隙比h*為機(jī)翼后緣到水平地面高度h (或到波浪地面波浪節(jié)點(diǎn)高度h)與機(jī)翼弦長(zhǎng)c的比值(h*=h/c)。本實(shí)驗(yàn)中水平地面間隙比工況為0.1、0.15、0.2、0.3、0.5、1.0。波浪地面間隙比工況為0.2、0.3、0.5(波浪地面時(shí),下翼面最低點(diǎn)低于后緣6 mm,間隙比h*=0.1~0.15時(shí)下翼面最低點(diǎn)會(huì)觸碰波峰表面,因此波浪地面間隙比工況最低為h*=0.2),實(shí)驗(yàn)幾何攻角α=0°。

    2 實(shí)驗(yàn)結(jié)果與討論

    2.1 水平地面測(cè)力結(jié)果

    圖2所示為0°攻角時(shí),機(jī)翼在水平地面附近,不同間隙比h*下的升力/阻力系數(shù)曲線。可以看出水平地面時(shí),隨著間隙比的逐漸降低,升力系數(shù)CL-flat隨之減小;首先在h*=1.0~0.5段CL-flat的下降十分緩慢,升力方向向上;h*<0.5后CL-flat下降的程度加劇但機(jī)翼仍承受正向升力,并在h*=0.3附近升力接近于零;在h*<0.2時(shí)CL-flat劇烈減小并呈負(fù)值,這意味著0°攻角間隙比h*小于0.3后,機(jī)翼開(kāi)始受到負(fù)升力的作用,且越靠近地面,負(fù)升力的作用越強(qiáng)。而阻力系數(shù)CD-flat曲線顯示,隨著間隙比的降低,阻力系數(shù)在較大h*范圍內(nèi)(h*=0.3~1.0)基本保持不變,即該h*范圍內(nèi)阻力系數(shù)CD-flat受地面效應(yīng)的影響非常弱;只是在h*<0.2之后,CD-flat有略有增大。該測(cè)力結(jié)果與以往文獻(xiàn)中所報(bào)道的實(shí)驗(yàn)或數(shù)值計(jì)算結(jié)果相吻合。波浪地面時(shí),對(duì)測(cè)得的瞬時(shí)升力阻力曲線取周期平均值,可以看到,得到的升力系數(shù)周期均值CL-wavy和阻力系數(shù)周期均值CD-wavy整體趨勢(shì)與水平地面時(shí)相近,數(shù)值上均略大于水平地面結(jié)果。

    圖2 水平地面和波浪地面升力系數(shù)與阻力系數(shù)(α=0°)Fig.2 Variation of the lift and drag coefficient CLand CD(α=0°)

    2.2 翼尖渦瞬時(shí)流場(chǎng)結(jié)構(gòu)

    本節(jié)討論各工況翼尖渦結(jié)構(gòu)特性及差異,觀察方向與機(jī)翼運(yùn)動(dòng)方向相同。速度矢量圖中坐標(biāo)參數(shù)均已無(wú)量綱化,橫坐標(biāo)x/c為實(shí)際空間寬度與弦長(zhǎng)之比,機(jī)翼正中位于x/c=0處,右、左側(cè)翼梢水平坐標(biāo)分別為x/c=±1處??v坐標(biāo)y/c為實(shí)際空間高度與弦長(zhǎng)之比,y/c=0標(biāo)識(shí)位置為水平地面,機(jī)翼垂直坐標(biāo)y/c在各間隙工況有所不同。坐標(biāo)z/c為下游截面離開(kāi)機(jī)翼后緣的距離與弦長(zhǎng)之比,z/c=0為機(jī)翼后緣所在截面。矢量圖中部空白處是由于受到垂直導(dǎo)板繞流對(duì)流場(chǎng)產(chǎn)生的干擾數(shù)據(jù),因此覆蓋該部分?jǐn)?shù)據(jù)應(yīng)不會(huì)對(duì)翼尖渦結(jié)構(gòu)的分析討論造成影響。

    2.2.1 不同間隙比翼尖渦速度場(chǎng)對(duì)比

    圖3為0°攻角,間隙比h*=1.0時(shí)機(jī)翼后緣下游z/c=1.0截面的速度矢量圖。圖中虛線標(biāo)識(shí)了機(jī)翼后緣寬度及所在高度??梢杂^察到在機(jī)翼兩端有兩個(gè)旋轉(zhuǎn)方向相反的旋渦,左側(cè)為順時(shí)針旋轉(zhuǎn),右側(cè)為逆時(shí)針旋轉(zhuǎn)。機(jī)翼尾流下洗明顯,說(shuō)明由于下翼面的壓強(qiáng)大于上翼面,使下翼面流體繞過(guò)兩側(cè)翼尖向上翼面翻卷形成翼尖渦,從而兩旋渦之間的流體誘導(dǎo)出強(qiáng)烈的下洗流動(dòng),因此機(jī)翼有效攻角小于0°,但仍大于零升攻角,使得機(jī)翼的升力為正值。因距離地面較遠(yuǎn),旋渦的生成和發(fā)展沒(méi)有受到明顯抑制,其形狀大致為較規(guī)則的圓形,兩渦核坐標(biāo)基本與兩翼尖坐標(biāo)重合;另外,由于誘導(dǎo)作用,兩旋渦之間的流體速度明顯大于旋渦外側(cè)的流速。對(duì)比測(cè)力結(jié)果,也說(shuō)明此時(shí)機(jī)翼受到正向升力的作用。

    圖4為0°攻角,間隙比h*=0.3時(shí)機(jī)翼后緣下游z/c=1.0截面的速度矢量圖??梢杂^察到該截面全場(chǎng)速度矢量相對(duì)于h*=1.0的工況顯著減小,機(jī)翼的兩端仍然有翼尖渦生成,但強(qiáng)度已大為減弱。兩側(cè)旋渦的旋轉(zhuǎn)方向與h*=1.0時(shí)相同,其誘導(dǎo)的下洗運(yùn)動(dòng)明顯削弱,同時(shí)旋渦的形狀因受地面擠壓,其輪廓已不再是規(guī)則的圓形,且兩渦之間伴有若干小尺度二次渦。由于流體從下翼面繞過(guò)翼尖向上翼面翻轉(zhuǎn),說(shuō)明下翼面壓力仍大于上翼面,但是翼尖渦強(qiáng)度的降低說(shuō)明上下翼面間的壓差減小,升力降低。這與測(cè)力結(jié)果中h*=0.3時(shí)升力較h*=1.0時(shí)減小的變化趨勢(shì)相吻合。

    圖5為0°攻角,間隙比h*=0.1時(shí)機(jī)翼后緣下游z/c=1.0截面的速度矢量圖。我們注意到該間隙比時(shí),流場(chǎng)結(jié)構(gòu)與前兩工況完全不同。此時(shí)機(jī)翼尾流有明顯上洗,流體從上翼面繞過(guò)兩側(cè)翼尖向下翼面翻卷形成翼尖渦。即左側(cè)旋渦逆時(shí)針旋轉(zhuǎn),右側(cè)旋渦順時(shí)針旋轉(zhuǎn),與h*=1.0、h*=0.3兩工況時(shí)的旋渦旋轉(zhuǎn)方向恰好相反。同時(shí)兩渦核位置已經(jīng)從翼尖向機(jī)翼內(nèi)側(cè)移動(dòng)了約0.2倍弦長(zhǎng)距離,其間距明顯小于h*=1.0時(shí)的情形。

    圖3 α=0°,h*=1.0,z/c=1.0時(shí)的速度矢量圖Fig.3 Tip vortex velocity fields at α=0°,h*=1.0,z/c=1.0

    圖4 α=0°,h*=0.3,z/c=1.0時(shí)的速度矢量圖Fig.4 Tip vortex velocity fields at α=0°,h*=0.3,z/c=1.0

    圖5 α=0°,h*=0.1,z/c=1.0時(shí)的速度矢量圖Fig.5 Tip vortex velocity fields at α=0°,h*=0.1,z/c=1.0

    圖5進(jìn)一步與圖3比較,兩工況的流場(chǎng)結(jié)構(gòu)有明顯區(qū)別:1)旋渦的形狀不再是圓形,由于受到地面的擠壓,旋渦的下半部分變扁;2)同側(cè)旋渦的旋轉(zhuǎn)方向相反,說(shuō)明兩種工況下上下翼面的壓力分布發(fā)生了相應(yīng)的變化;小間隙比時(shí),上翼面壓力大于下翼面,機(jī)翼受到負(fù)升力的作用,升力方向向下;3)雖然是處于下游的同一距離截面,但渦的移動(dòng)方向及流場(chǎng)速度也有著明顯的差異;在高間隙比h*=1.0時(shí),翼尖渦向下游延伸的過(guò)程中,渦核不產(chǎn)生水平或垂向的移動(dòng),渦核間距基本保持恒定;而低間隙比h*=0.1時(shí),旋渦在翼尖生成后,兩渦核分別向機(jī)翼內(nèi)側(cè)產(chǎn)生明顯移動(dòng),越向下游,兩渦核間距越近。

    對(duì)比測(cè)力曲線,流場(chǎng)的速度分布狀況印證了機(jī)翼受到負(fù)升力的合理性,下翼面流速大于上翼面,導(dǎo)致下翼面壓力變?yōu)槲?,機(jī)翼整體受負(fù)升力,與測(cè)力結(jié)果吻合。另外,在h*=0.1時(shí),由于展向的流動(dòng)較強(qiáng),在地面表面分離形成小尺度二次渦。

    2.2.2 不同間隙比翼尖渦渦量演化過(guò)程對(duì)比

    圖6~圖11為0°攻角h*=0.1~1.0各間隙比時(shí),機(jī)翼下游各截面翼尖渦瞬時(shí)渦量分布演化過(guò)程的云圖,每個(gè)間隙比均截取了機(jī)翼下游z/c=0.08、0.48、1.0、2.0、3.0、4.0六個(gè)截面的渦量分布,每個(gè)截面上均有一對(duì)渦量相反的集中渦分布,圖中坐標(biāo)z/c=0與機(jī)翼后緣重合。

    圖6 α=0°,h*=1.0時(shí)的不同截面渦核渦量分布Fig.6 Tip vortex vorticity contours at α=0°,h*=1.0

    圖7 α=0°,h*=0.5時(shí)的不同截面渦核渦量分布Fig.7 Tip vortex vorticity contours at α=0°,h*=0.5

    圖8 α=0°,h*=0.3時(shí)的不同截面渦核渦量分布Fig.8 Tip vortex vorticity contours at α=0°,h*=0.3

    圖9 α=0°,h*=0.2時(shí)的不同截面渦核渦量分布Fig.9 Tip vortex vorticity contours at α=0°,h*=0.2

    圖10 α=0°,h*=0.15時(shí)的不同截面渦核渦量分布Fig.10 Tip vortex vorticity contours at α=0°,h*=0.15

    圖11 α=0°,h*=0.1時(shí)的不同截面渦核渦量分布Fig.11 Tip vortex vorticity contours at α=0°,h*=0.1

    圖6~圖11可以看出翼尖渦渦量Ω總的變化趨勢(shì)是:隨著間隙比的減小,渦核處的最大渦量也相應(yīng)下降,表明機(jī)翼越接近地面,翼尖渦的生成和發(fā)展所受的抑制越強(qiáng),并且在向下游延伸的過(guò)程中,旋渦的耗散速度也越快。h*=1.0時(shí)翼尖渦渦核處的最大渦量Ωmax衰減到z/c=4.0截面時(shí)的變化量ΔΩ=7;h*=0.5時(shí)為ΔΩ=9;h*=0.3時(shí)衰減量為ΔΩ=15。而h*=0.2時(shí)的情況較為特殊,渦量衰減量ΔΩ僅為6,這是由于該小間隙比時(shí),初始時(shí)刻渦的生成已經(jīng)受到強(qiáng)烈抑制,渦核處最大渦量絕對(duì)值很小,機(jī)翼兩端幾乎沒(méi)有輪廓清晰的集中渦形成,當(dāng)衰減到z/c=4.0截面時(shí),渦量幾乎與靜止水體的渦量相同,翼尖渦被完全耗散,因此ΔΩ較小。h*=0.15時(shí)衰減量ΔΩ=12,h*=0.1時(shí)衰減量為ΔΩ=15。分析原因如下:從圖6~圖9中看出,隨著間隙比的減小,旋渦強(qiáng)度只是表現(xiàn)出渦量上的變化和耗散速度的區(qū)別,并無(wú)流動(dòng)結(jié)構(gòu)上的根本改變,機(jī)翼兩端的集中渦渦量都是左側(cè)為負(fù)、右側(cè)為正,說(shuō)明在機(jī)翼的左端旋渦為順時(shí)針、右端為逆時(shí)針,即水流從下翼面經(jīng)翼尖翻轉(zhuǎn)至上翼面,下翼面的壓力大于上翼面,兩側(cè)翼尖渦渦核軸線均平行向下延伸,兩渦間距保持不變。然而當(dāng)間隙比減小到h*=0.15和h*=0.1時(shí),如圖10、圖11所示,流動(dòng)結(jié)構(gòu)及渦量分布發(fā)生了本質(zhì)的變化,可以發(fā)現(xiàn)渦量在左側(cè)變?yōu)檎?,右?cè)為負(fù)值,說(shuō)明在機(jī)翼的左端旋渦為逆時(shí)針旋轉(zhuǎn),即上翼面水流繞過(guò)翼尖向下翼面旋轉(zhuǎn)形成翼尖渦。同時(shí)在向下游發(fā)展的過(guò)程中,兩端渦核不斷向機(jī)翼內(nèi)側(cè)運(yùn)動(dòng),渦量也逐步耗散。另外從這兩個(gè)工況的結(jié)果看出,機(jī)翼越靠近地面,翼尖渦的強(qiáng)度越大,渦核向機(jī)翼內(nèi)側(cè)的位移也越大,兩渦越靠近。

    對(duì)于機(jī)翼旋渦變化的現(xiàn)象,做如下解釋:圖12所示為地面附近的0°攻角機(jī)翼剖面圖,x軸代表水平地面??梢钥闯鰴C(jī)翼下翼面與地面構(gòu)成的流道先收縮后擴(kuò)張。在流道形狀固定的情況下,隨著間隙比的減小,流道逐漸變得狹窄,使得下翼面的流速增加,間隙比越小流速增加越大,當(dāng)間隙比降為h*=0.2時(shí),上下翼面流速接近,壓力差趨于零,機(jī)翼兩端沒(méi)有明顯的旋渦生成。h*<0.2后,下翼面的流速即超過(guò)上翼面。同時(shí),由于機(jī)翼逐漸靠近地面,產(chǎn)生的翼尖渦受到地面的擠壓變形,在一定程度上抑制了翼尖渦的生成并且也加快了翼尖渦的耗散速度。

    圖12 機(jī)翼及與地面形成的通道示意圖Fig.12 Side view of the wing and floor

    2.2.3 波浪地面不同間隙比渦量對(duì)比

    在正弦波浪地形不同相位的截面上,翼尖渦流場(chǎng)結(jié)構(gòu)及其隨時(shí)間的演化過(guò)程均有所區(qū)別,因此本實(shí)驗(yàn)研究了位于波浪地形的波峰、波谷和前后波節(jié)等典型截面處的的翼尖渦渦量隨時(shí)間的變化規(guī)律,以機(jī)翼后緣掠過(guò)該截面的時(shí)刻為計(jì)時(shí)零點(diǎn)。首先給出每個(gè)截面計(jì)時(shí)起始時(shí)刻機(jī)翼后緣相對(duì)于波浪地面相位的示意圖,如圖13所示。波浪地形波長(zhǎng)為1.6 c,波幅0.1 c。此時(shí)定義平均間隙比h*為機(jī)翼后緣至波節(jié)的高度與弦長(zhǎng)之比,因此后緣在波谷處的間隙比H*為H*=h*+0.1;后緣在波峰處的間隙比H*為H*=h*-0.1。

    圖13 機(jī)翼后緣截面與波浪相對(duì)位置關(guān)系示意圖(α=0°)Fig.13 Schematic diagram of the relative position of the trailing edge and the wavesurface(α=0°)

    圖14~圖16為0°攻角,h*=0.2~0.5時(shí)后緣下游各截面(z/c=0.08,0.48,1.0,2.0,3.0,4.0)渦核渦量Ωcore在不同相位對(duì)比圖。該渦量是將機(jī)翼兩側(cè)翼尖渦核極值取絕對(duì)值后的平均值結(jié)果。首先比較z/c= 0.08截面各間隙比時(shí)結(jié)果;圖14中可以看到當(dāng)后緣經(jīng)過(guò)波谷時(shí),z/c=0.08截面翼尖渦渦量為四個(gè)相位中的最大值,渦量大約為相同間隙比時(shí)水平地面值的4倍,后緣此時(shí)離開(kāi)曲面實(shí)際高度為h*=0.3,是h*= 0.3時(shí)水平地面渦量值的2倍,說(shuō)明在機(jī)翼后緣掠過(guò)波谷的瞬間,機(jī)翼后半段與波浪地形構(gòu)成的擴(kuò)張通道不但沒(méi)有抑制翼尖渦的生成,反而增大了其渦量值。當(dāng)后緣掠過(guò)波峰后,z/c=0.08截面翼尖渦渦量最大值是同間隙比時(shí)水平地面結(jié)果的3倍;此時(shí)后緣與曲面的實(shí)際高度差為H*=h*-0.1=0.1,但是此時(shí)并沒(méi)有出現(xiàn)水平地面時(shí)翼尖渦反向旋轉(zhuǎn)的結(jié)果。再對(duì)比后緣經(jīng)過(guò)后波節(jié)時(shí)的結(jié)果,得到與前兩個(gè)相位類似的結(jié)果。說(shuō)明這三個(gè)相位時(shí)下翼面流道在小間隙比時(shí)對(duì)地面效應(yīng)的影響都是削弱的。再看圖13中前節(jié)點(diǎn)渦量變化曲線,其與水平地面相同間隙比時(shí)的曲線幾乎重合,因此也可認(rèn)為當(dāng)機(jī)翼經(jīng)過(guò)前波節(jié)的瞬間,下翼面流道的影響與水平地面時(shí)構(gòu)成的流道對(duì)地面效應(yīng)的影響相同,均為先收縮后擴(kuò)張的通道。圖15、圖16中,除了后緣位于前波節(jié)時(shí)的渦量結(jié)果與水平地面相同間隙比時(shí)的情況吻合外,后緣位于其他三個(gè)相位(波峰、波谷、后波節(jié))的結(jié)果均遠(yuǎn)大于水平地面相同間隙比時(shí)的結(jié)果,且該三個(gè)相位的最大渦量之間的差距比較小,說(shuō)明這三個(gè)相位在較大間隙比時(shí)下翼面與波浪地面形成的流道雖然形狀不盡相似,但地面效應(yīng)的影響均已削弱且效果相近。值得注意的是,這三個(gè)相位在較大間隙比時(shí)的渦量極值均明顯大于水平地面h*=1.0時(shí)的結(jié)果。圖14中,波峰、波谷、后節(jié)點(diǎn)三個(gè)相位處的渦量變化如曲線所示,可以看出在平均間隙比h*=0.2時(shí),隨著機(jī)翼遠(yuǎn)離該三截面,翼尖渦渦量的耗散明顯加快。圖15、圖16中,各相位處渦量變化曲線下降并不明顯,說(shuō)明在間隙比稍大時(shí),翼尖渦渦量的耗散速度沒(méi)有加快,波浪地面對(duì)地面效應(yīng)的影響亦削弱。

    取相對(duì)于波浪地面不同相位的各z/c=0.08截面渦核處渦量值Ωcore并繪成曲線,同時(shí)加入機(jī)翼波浪地面瞬時(shí)升力系數(shù)與水平地面升力系數(shù)比值的變化曲線,如圖17所示。圖17(a)為機(jī)翼波浪地面瞬時(shí)升力系數(shù)比CL-wavy/CD-flat變化曲線,圖17(b)為渦量變化曲線,橫坐標(biāo)表示機(jī)翼后緣下方波浪地面的相位,橫坐標(biāo)0值處為后波節(jié)。由于在三個(gè)不同間隙比下,瞬時(shí)升力系數(shù)趨勢(shì)相同,此處僅討論升力與波浪地面相位關(guān)系,故只畫(huà)出一條升力變化趨勢(shì)線。

    圖14 α=0°,h*=0.2時(shí)的不同截面渦核渦量Fig.14 Variation of the vortex cores vorticity at α=0°,h*=0.2

    圖15 α=0°,h*=0.3時(shí)的不同截面渦核渦量Fig.15 Variation of the vortex cores vorticity at α=0°,h*=0.3

    圖16 α=0°,h*=0.5時(shí)的不同截面渦核渦量Fig.16 Variation of the vortex cores vorticity at α=0°,h*=0.5

    圖17 α=0°,不同間隙比下渦量和升力系數(shù)比變化與波浪相對(duì)位置關(guān)系Fig.17 Variation of the vortex cores vorticityand ratio of CLwith the phase of the wave surface(α=0°)

    可以看到在h*=0.2時(shí),渦量在機(jī)翼后緣經(jīng)過(guò)前節(jié)點(diǎn)時(shí)達(dá)到最小值,隨著機(jī)翼后緣不斷向前運(yùn)動(dòng)掠過(guò)波峰、后節(jié)點(diǎn)和波谷,渦量不斷上升至波谷時(shí)的最大值,然后再次迅速跌落至前節(jié)點(diǎn)處的最小值。當(dāng)h*=0.3~0.5時(shí),渦量亦在前節(jié)點(diǎn)達(dá)到最小值,但是隨即在波峰處恢復(fù)到接近最大值的水平,在機(jī)翼后緣掠過(guò)后節(jié)點(diǎn)時(shí)達(dá)到最大值,在機(jī)翼后緣越過(guò)波谷后迅速下降,運(yùn)動(dòng)到前節(jié)點(diǎn)時(shí)渦量達(dá)到最小值。渦量隨波浪地面相對(duì)位置的變化結(jié)果顯示,在一個(gè)地形波浪周期內(nèi),渦核處渦量變化表現(xiàn)出明顯周期性變化,但是不符合波浪地面的正弦變化規(guī)律。因此,從波浪地面與水平地面的渦核渦量Ωcore對(duì)比結(jié)果來(lái)看,除了在上坡節(jié)點(diǎn)(前節(jié)點(diǎn))處兩者相近外,波浪地面對(duì)于機(jī)翼的地面效應(yīng)的影響較水平地面情況是削弱的。從波浪地面的不同相位結(jié)果來(lái)比較,上坡節(jié)點(diǎn)處受地面影響最大,此時(shí)流道形狀最狹窄。在小間隙比條件下,其余節(jié)點(diǎn)處受地面影響程度不同,隨著間隙比增大,受地面影響程度的差異迅速減小,這也是渦量測(cè)量結(jié)果不符合正弦規(guī)律的原因。而在瞬時(shí)升力系數(shù)比曲線上,可以發(fā)現(xiàn)當(dāng)機(jī)翼后緣在后節(jié)點(diǎn)與波谷之間時(shí)升力達(dá)到最大值,后緣在前節(jié)點(diǎn)與波峰之間時(shí)升力達(dá)最小值,同時(shí),在機(jī)翼升力增加的過(guò)程中,幅值變化明顯,表現(xiàn)為在最大值附近曲率大而在最小值附近曲率很小。此時(shí)可以認(rèn)為在后緣上坡的過(guò)程中,下翼面與地面形成的流道造成了下翼面流速的增加,導(dǎo)致上下翼面壓力差減小,升力損失,旋渦強(qiáng)度減弱。當(dāng)機(jī)翼后緣下坡時(shí),后緣與地面間隙較小,下翼面流體較難迅速?gòu)牧鞯懒鞒?,造成下翼面流體受到阻滯,速度減小,升力增加,旋渦強(qiáng)度增加。

    3 結(jié) 論

    利用拖曳水槽模擬了機(jī)翼在水平及波浪表面附近的地面效應(yīng),采用二維高速PIV系統(tǒng)測(cè)量了各工況下尾流場(chǎng)的瞬時(shí)變化過(guò)程。研究認(rèn)為在0°攻角下翼尖渦流場(chǎng)結(jié)構(gòu)受到地面與下翼面所構(gòu)成的流道形狀及狹窄程度的影響;并發(fā)現(xiàn)在水平地面附近不同飛行高度時(shí),機(jī)翼所受正升力和負(fù)升力所對(duì)應(yīng)的流場(chǎng)結(jié)構(gòu)差異,并分析了機(jī)翼受負(fù)升力的原因。在波浪地面的測(cè)量中,獲得了機(jī)翼在飛躍正弦波浪地形的過(guò)程中流場(chǎng)結(jié)構(gòu)的變化,地面效應(yīng)由于波浪地面與機(jī)翼下翼面瞬間相對(duì)相位關(guān)系的不同而不同,僅在機(jī)翼經(jīng)過(guò)上坡節(jié)點(diǎn)時(shí),地面效應(yīng)與水平地面時(shí)相似;而在其他波浪相位時(shí)地面效應(yīng)均受到削弱。同時(shí)發(fā)現(xiàn)正弦波浪上各相位截面翼尖渦核處渦量的變化呈現(xiàn)出明顯的周期性,但該周期性不符合正弦變化規(guī)律,且其極值與正弦波型有相位差。這些特性對(duì)于地面效應(yīng)飛行器縱向穩(wěn)定性的控制具有一定的參考價(jià)值。

    [1] Li Dun,WangYining.Three-dimension numerical simulation of the ground effect on the ground effect vehicle[J].Flight Dynamics,2006,24(2):18-21.(in Chinese)

    李盾,王義寧.地效飛行器三維地面效應(yīng)的數(shù)值模擬[J].飛行力學(xué),2006,24(2):18-21.

    [2] Ahmed M R,Takasaki T,Kohama Y.Aerodynamics of a NACA4412 airfoil in ground effect[J].AIAA Journal,2007,45(1): 37-47.

    [3] Chen Xin,Shan Xuexiong.Numerical simulation on unsteady aerodynamics and wace making of a wing flying over offing[C]//2003 Symposium on Research Frontier of Aerodynamics.Beijing:China Astronautic Publishing House,2003:132-137.(in Chinese)陳新,單雪雄.三維機(jī)翼掠海飛行時(shí)非定常氣動(dòng)力和興波的數(shù)值計(jì)算[C]//2003空氣動(dòng)力學(xué)前沿研究論文集.北京:中國(guó)宇航出版社,2003:132-137.

    [4] Byelinskyy V G,Zinchuk P I.Hydrodynamical characteristics of an ekranoplane wing flying near the wavy sea surface[C]//RTO Meeting proceedings,1999.

    [5] Zhang Xin,Zeriban J,Ruhrmann A,et al.Tip vortices generated by a wing in ground effect[C]//Proceedings of the 11thInternational Symposium on Applications of laser Techniques to Fluid Mechanics.University of Southampton,2007.

    [6] Jung K H,Chun H H,Kim H J.Experimental investigation of wingin-ground effect with a NACA6409 section[J].Journal of Marine Science and Technology,2008,13:317-327.

    [7] Hiemcke C.NACA5312 in ground effect:wind tunnel and panel code studies[R].AIAA 97-2320.

    [8] Qu Qiulin,Liu Peiqing,Qin Xuguo.Numerical research on separated flow around a WIG craft flying over wavy ground at high incidence angle[J].Acta Aerodynamica Sinica,2008,26(2):221-226.(in Chinese)

    屈秋林,劉沛清,秦緒國(guó).地效飛行器近波浪地面大迎角飛行分離流動(dòng)數(shù)值研究[J].空氣動(dòng)力學(xué)學(xué)報(bào),2008,26(2):221-226.

    [9] Qu Qiulin,Liu Peiqing,Qin Xuguo.Numerical research on aerodynamic characteristics of WIG craft flying over wavy ground[J].Acta Aeronautica et Astronauticasinica,2007,28(6):1327-1333.(in Chinese)

    屈秋林,劉沛清,秦緒國(guó).地效飛行器波浪地面飛行氣動(dòng)性能數(shù)值研究[J].航空學(xué)報(bào),2007,28(6):1327-1333.

    [10]Im Y H,Chang K S.Unsteady aerodynamics of a wing in ground effect airfoil flying over a wavy wall[J].Journal of Aircraft,2000,37(4):690-696.

    [11]Sun Ruimin,Daichin.Experimental investigation on tip vortices and aerodynamics of a wing in ground effect[J].Theoretical and Applied Mechanics Letters,2011,1(3):1-6.

    [12]Yang Wei,Yang Zhigang.Aerodynamic investigation of a 2D wing and flows in ground effect[J].Chinese Journal of Computational Physics,2009,26(2):231-240.

    [13]Rozhdestvensky K V.Wing-in ground effect vehicles[J].Progress in Aerospace Sciences,2006,42:211-283.

    [14]Zhang Xin,Zeriban J.Aerodynamics of a single element wing in ground effect[J].Journal of Aircraft,2000,37(6):1058-1064.

    [15]Harvey J K,Perry F J.Flowfield produced by trailing vortices in the vicinity of the ground[J].AIAA Journal,1971,9(8):1659-1660.

    Experimental study on the tip vortices of a wing close to a flat and a wavy surface using PIV

    Zhang Kuang1,Daichin1,2,*
    (1.Shanghai Institute of Applied Mathematics and Mechanics,Shanghai University,Shanghai 200072,China; 2.Shanghai Key Laboratory of Mechanics in Energy Engineering,Shanghai 200072,China)

    An experimental study on the tip vortices and lift/drag forces of a NACA23012 wing in adjacent to the flat and the sine wavy surface at zero angle of attack is introduced.The experiment is conducted in a towing tank.The study focuses on the hydrodynamics and evolution of the tip vortices of the wing close to the flat and the sine wavy surface with different ground height,velocity and vorticity distribution,as well as their influence on the lift/drag force of the wing are compared.The results show that,the lift force decreases from positive to negative as the wing gets closer to the wavy surface,the intensity of the tip vortices change as well,and the shape of tip vorticesis is deformed correspondingly.Especially,for the case of small ground height and negative lift force,the rotation direction of the tip vortices is changed.The flow structure is not only affected by the wing ground height,but also influenced by the phase relation of the wing and the ground,i.e.the instantaneous shape of the flow channel between the lower surface of the wing and the wavy ground surface.And the vorticity of the tip vortices along the wavy surface presents some periodic variation,whereas the variation does not comply with sine regularity.The reason is that the periodical variation of the shape of the flow channel and the stenosis degree of the channel beneath the wing brings a different level of inhabitation to the evolution of the tip vortices.

    wing in ground effect;tip vortices;PIV;sine wavy surface;towing tank;experimental study

    V211.7

    A

    10.7638/kqdlxxb-2013.0049

    0258-1825(2015)03-0367-09

    2013-04-27;

    2013-06-07

    國(guó)家自然科學(xué)基金(11472169;11072142)

    章曠(1987-),男,浙江臺(tái)州人,碩士研究生,研究方向:實(shí)驗(yàn)流體力學(xué)和機(jī)翼地面效應(yīng)空氣動(dòng)力學(xué).E-mail:toughstone3@gmail.com

    代欽*(1966-),男,蒙古族,內(nèi)蒙古錫林郭勒盟人,工學(xué)博士,研究員,研究方向:實(shí)驗(yàn)流體力學(xué).E-mail:daichin@staff.shu.edu.cn

    章曠,代欽.地面效應(yīng)作用下翼尖渦特性的PIV實(shí)驗(yàn)研究[J].空氣動(dòng)力學(xué)學(xué)報(bào),2015,33(3):367-374.

    10.7638/kqdlxxb-2013.0049 Zhang K,Daichin.Experimental study on the tip vortices of a wing close to a flat and a wavy surface using PIV[J].Acta Aerodynamica Sinica,2015,33(3):367-374.

    猜你喜歡
    翼尖翼面渦量
    基于拔銷器鎖定的飛行器氣動(dòng)控制面解鎖控制方法
    中高速條件下不同翼尖小翼的數(shù)值模擬分析
    含沙空化對(duì)軸流泵內(nèi)渦量分布的影響
    固定翼二維彈道修正引信升力翼面位置的影響
    自由表面渦流動(dòng)現(xiàn)象的數(shù)值模擬
    基于Aerobook平臺(tái)的復(fù)合材料翼面結(jié)構(gòu)設(shè)計(jì)流程
    基于翼尖渦物理特征的誘導(dǎo)阻力減阻機(jī)制實(shí)驗(yàn)研究
    基于流動(dòng)顯示的翼尖渦不穩(wěn)定頻率測(cè)量
    基于 FFD 技術(shù)的民用運(yùn)輸機(jī)翼尖裝置設(shè)計(jì)研究
    航態(tài)對(duì)大型船舶甲板氣流場(chǎng)的影響
    亚洲av成人一区二区三| 亚洲成人精品中文字幕电影 | 国产一区二区三区在线臀色熟女 | 亚洲自拍偷在线| 久久精品影院6| 国产精品久久视频播放| 久久久国产一区二区| 久久精品国产综合久久久| 久久久久久免费高清国产稀缺| 美女午夜性视频免费| 国产av精品麻豆| 别揉我奶头~嗯~啊~动态视频| 级片在线观看| 一级a爱视频在线免费观看| 日韩精品青青久久久久久| 国产xxxxx性猛交| 搡老岳熟女国产| 男女做爰动态图高潮gif福利片 | 欧美日韩黄片免| 天天躁夜夜躁狠狠躁躁| 国产主播在线观看一区二区| 天堂动漫精品| www.999成人在线观看| 大陆偷拍与自拍| 日韩欧美免费精品| 男女之事视频高清在线观看| 久久亚洲真实| 午夜免费激情av| 久久久国产成人免费| 日韩一卡2卡3卡4卡2021年| 黄色丝袜av网址大全| 亚洲自偷自拍图片 自拍| 99国产极品粉嫩在线观看| 久久香蕉国产精品| 精品国产一区二区久久| 91九色精品人成在线观看| 自线自在国产av| 天堂√8在线中文| 国产精华一区二区三区| 亚洲色图综合在线观看| 自线自在国产av| 亚洲精品粉嫩美女一区| 人人妻人人澡人人看| 一级毛片女人18水好多| 欧美+亚洲+日韩+国产| 免费高清视频大片| 欧美日韩视频精品一区| 最新美女视频免费是黄的| 天堂影院成人在线观看| 他把我摸到了高潮在线观看| 精品人妻在线不人妻| 久久精品人人爽人人爽视色| 高清黄色对白视频在线免费看| 欧美日韩视频精品一区| 成熟少妇高潮喷水视频| 搡老乐熟女国产| 黑人操中国人逼视频| 国产精品国产高清国产av| 成年版毛片免费区| 中文字幕精品免费在线观看视频| 久久久久国产精品人妻aⅴ院| 国产精品亚洲一级av第二区| 久久久久久人人人人人| 国产成年人精品一区二区 | www日本在线高清视频| 涩涩av久久男人的天堂| 国产亚洲精品久久久久5区| 99久久国产精品久久久| 欧美乱色亚洲激情| 精品国产一区二区三区四区第35| 极品教师在线免费播放| 国产精品美女特级片免费视频播放器 | 美女扒开内裤让男人捅视频| 女同久久另类99精品国产91| 99热只有精品国产| 亚洲人成电影观看| 欧美中文日本在线观看视频| 中文亚洲av片在线观看爽| 精品乱码久久久久久99久播| 国产精品成人在线| 成人精品一区二区免费| 一级a爱片免费观看的视频| 亚洲欧美一区二区三区黑人| 欧美日韩亚洲国产一区二区在线观看| 在线观看免费视频日本深夜| 国产一卡二卡三卡精品| 精品久久久久久,| 成人三级黄色视频| 交换朋友夫妻互换小说| 久久久久久久久免费视频了| 精品电影一区二区在线| 日本精品一区二区三区蜜桃| 精品一区二区三区四区五区乱码| 色综合婷婷激情| 男女午夜视频在线观看| 啪啪无遮挡十八禁网站| 精品国产乱子伦一区二区三区| 成人av一区二区三区在线看| 欧美日韩乱码在线| 欧美一级毛片孕妇| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久人人做人人爽| 又大又爽又粗| 看黄色毛片网站| 免费看a级黄色片| 日本免费a在线| 99久久综合精品五月天人人| 久久香蕉精品热| 一二三四社区在线视频社区8| 亚洲人成电影免费在线| 不卡av一区二区三区| 亚洲国产看品久久| 后天国语完整版免费观看| 亚洲精品av麻豆狂野| 日韩欧美国产一区二区入口| 日日干狠狠操夜夜爽| 欧美中文综合在线视频| 精品电影一区二区在线| 最近最新中文字幕大全免费视频| 国产精品永久免费网站| 精品一区二区三区视频在线观看免费 | 两性夫妻黄色片| 两个人看的免费小视频| 涩涩av久久男人的天堂| 大型av网站在线播放| 成人影院久久| 精品电影一区二区在线| 亚洲精华国产精华精| 久久人人爽av亚洲精品天堂| 神马国产精品三级电影在线观看 | 国产99白浆流出| 国产一区二区激情短视频| 热99re8久久精品国产| 国产成+人综合+亚洲专区| 日韩成人在线观看一区二区三区| 一区二区三区精品91| 日日爽夜夜爽网站| 午夜a级毛片| 这个男人来自地球电影免费观看| 精品卡一卡二卡四卡免费| 欧美av亚洲av综合av国产av| 久久精品成人免费网站| 美女福利国产在线| 亚洲国产毛片av蜜桃av| av在线天堂中文字幕 | 91精品国产国语对白视频| 色老头精品视频在线观看| ponron亚洲| 国产不卡一卡二| 久久香蕉国产精品| 18禁裸乳无遮挡免费网站照片 | 欧美+亚洲+日韩+国产| 一进一出好大好爽视频| 黄频高清免费视频| 欧美日韩一级在线毛片| 成人18禁在线播放| 精品福利观看| 一a级毛片在线观看| 淫秽高清视频在线观看| 亚洲色图 男人天堂 中文字幕| 欧美一级毛片孕妇| 岛国在线观看网站| 亚洲精品在线观看二区| 色综合站精品国产| 国产精品美女特级片免费视频播放器 | 高清欧美精品videossex| 日本 av在线| 男人操女人黄网站| 无限看片的www在线观看| 国产精品国产av在线观看| 久9热在线精品视频| 精品久久久精品久久久| 国产欧美日韩一区二区三| 久久久国产成人免费| 成人永久免费在线观看视频| 另类亚洲欧美激情| 99久久精品国产亚洲精品| 手机成人av网站| 人人妻人人爽人人添夜夜欢视频| 狂野欧美激情性xxxx| 一级黄色大片毛片| 午夜福利在线免费观看网站| 正在播放国产对白刺激| 亚洲全国av大片| 满18在线观看网站| 国产精品野战在线观看 | 亚洲情色 制服丝袜| 久久久久久免费高清国产稀缺| 精品一品国产午夜福利视频| 满18在线观看网站| 久久热在线av| 美女大奶头视频| 亚洲五月婷婷丁香| 欧美激情久久久久久爽电影 | 国产欧美日韩精品亚洲av| 午夜老司机福利片| 天堂俺去俺来也www色官网| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看免费午夜福利视频| 手机成人av网站| 在线天堂中文资源库| 婷婷丁香在线五月| 波多野结衣高清无吗| 啪啪无遮挡十八禁网站| 精品久久久精品久久久| 国产一区二区三区综合在线观看| 国产精品久久电影中文字幕| 国产野战对白在线观看| 亚洲五月婷婷丁香| 9热在线视频观看99| 国产一区二区三区综合在线观看| 黑人猛操日本美女一级片| 日韩精品青青久久久久久| a级毛片在线看网站| 99国产精品99久久久久| 99热只有精品国产| 91麻豆av在线| 亚洲伊人色综图| 最近最新中文字幕大全免费视频| 色播在线永久视频| 99国产精品免费福利视频| 久久久国产成人精品二区 | 国产麻豆69| 老司机在亚洲福利影院| 久久中文字幕人妻熟女| 人妻丰满熟妇av一区二区三区| 亚洲国产精品999在线| 制服诱惑二区| 亚洲精品一卡2卡三卡4卡5卡| 90打野战视频偷拍视频| 久久热在线av| 精品人妻在线不人妻| 精品国产一区二区久久| 麻豆久久精品国产亚洲av | 欧美日韩黄片免| 在线看a的网站| 啦啦啦在线免费观看视频4| 一级a爱片免费观看的视频| 99在线视频只有这里精品首页| 757午夜福利合集在线观看| 19禁男女啪啪无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 国产成人系列免费观看| 国产麻豆69| 日韩中文字幕欧美一区二区| 亚洲黑人精品在线| 女人爽到高潮嗷嗷叫在线视频| 久久天堂一区二区三区四区| 一级片'在线观看视频| 91字幕亚洲| 日韩高清综合在线| 岛国在线观看网站| 叶爱在线成人免费视频播放| 啦啦啦在线免费观看视频4| 精品人妻在线不人妻| 日本vs欧美在线观看视频| 亚洲三区欧美一区| 国产精品 国内视频| 一区二区三区国产精品乱码| 日本精品一区二区三区蜜桃| 搡老乐熟女国产| 美女大奶头视频| 国产极品粉嫩免费观看在线| 久久精品91蜜桃| а√天堂www在线а√下载| 丁香六月欧美| 欧美午夜高清在线| 色综合欧美亚洲国产小说| 欧美日韩亚洲高清精品| 如日韩欧美国产精品一区二区三区| 热99国产精品久久久久久7| 天堂√8在线中文| 黑人猛操日本美女一级片| 亚洲美女黄片视频| 真人一进一出gif抽搐免费| 水蜜桃什么品种好| 黑丝袜美女国产一区| 最新美女视频免费是黄的| 国产xxxxx性猛交| 亚洲人成伊人成综合网2020| 91字幕亚洲| 成人手机av| www日本在线高清视频| 成人亚洲精品一区在线观看| 日日摸夜夜添夜夜添小说| 丰满的人妻完整版| 国产精品亚洲av一区麻豆| 12—13女人毛片做爰片一| 午夜福利在线观看吧| 美国免费a级毛片| 免费高清视频大片| 国产成人av激情在线播放| 精品一区二区三卡| 可以免费在线观看a视频的电影网站| 999久久久国产精品视频| 在线观看午夜福利视频| 亚洲精品一卡2卡三卡4卡5卡| 黄频高清免费视频| 国产成人精品久久二区二区免费| 亚洲中文av在线| 久久国产精品人妻蜜桃| 18禁国产床啪视频网站| 久久久久精品国产欧美久久久| 亚洲自拍偷在线| 国产成人av教育| 999久久久精品免费观看国产| 国产精品 国内视频| 国产免费男女视频| 精品国产乱码久久久久久男人| 91精品国产国语对白视频| 岛国视频午夜一区免费看| av网站在线播放免费| 黑人巨大精品欧美一区二区蜜桃| 日韩有码中文字幕| 18禁观看日本| 亚洲人成电影免费在线| ponron亚洲| 欧美成狂野欧美在线观看| 久热爱精品视频在线9| 亚洲精品中文字幕一二三四区| 久久精品国产综合久久久| 亚洲精品粉嫩美女一区| 日日爽夜夜爽网站| av有码第一页| 精品高清国产在线一区| 99在线人妻在线中文字幕| av天堂久久9| 久久亚洲精品不卡| svipshipincom国产片| 国产又爽黄色视频| 婷婷丁香在线五月| av在线播放免费不卡| 国产亚洲欧美98| ponron亚洲| 久9热在线精品视频| 激情视频va一区二区三区| 亚洲精华国产精华精| 久久久久国产精品人妻aⅴ院| 欧美激情高清一区二区三区| 老司机深夜福利视频在线观看| 欧美日本中文国产一区发布| 午夜免费鲁丝| www.999成人在线观看| 国产日韩一区二区三区精品不卡| 欧美黄色片欧美黄色片| www.熟女人妻精品国产| 久久精品国产清高在天天线| 精品午夜福利视频在线观看一区| 在线观看www视频免费| 这个男人来自地球电影免费观看| 亚洲熟女毛片儿| 18美女黄网站色大片免费观看| a级毛片黄视频| 狠狠狠狠99中文字幕| 欧美乱码精品一区二区三区| 大陆偷拍与自拍| 亚洲av成人一区二区三| 在线天堂中文资源库| 一级毛片女人18水好多| 国产一区二区激情短视频| 日日夜夜操网爽| 日韩三级视频一区二区三区| 麻豆一二三区av精品| 日本黄色日本黄色录像| 中文字幕av电影在线播放| 99精品欧美一区二区三区四区| 国产高清videossex| 精品国产一区二区三区四区第35| 国产成人啪精品午夜网站| 国产xxxxx性猛交| 99久久综合精品五月天人人| 不卡一级毛片| 亚洲黑人精品在线| 91九色精品人成在线观看| a级片在线免费高清观看视频| 很黄的视频免费| 亚洲色图av天堂| 五月开心婷婷网| 欧美中文综合在线视频| 村上凉子中文字幕在线| 满18在线观看网站| 人妻久久中文字幕网| 91老司机精品| 97人妻天天添夜夜摸| 男人舔女人下体高潮全视频| 亚洲一码二码三码区别大吗| 亚洲欧美激情综合另类| 国产成人av激情在线播放| 精品一区二区三区视频在线观看免费 | 97人妻天天添夜夜摸| 少妇被粗大的猛进出69影院| av福利片在线| 国产精品免费视频内射| 日本免费a在线| 麻豆一二三区av精品| 18禁国产床啪视频网站| 少妇的丰满在线观看| 丝袜在线中文字幕| 亚洲五月天丁香| 美国免费a级毛片| 黄色视频不卡| 国产av一区在线观看免费| 天天躁狠狠躁夜夜躁狠狠躁| av福利片在线| 中文亚洲av片在线观看爽| а√天堂www在线а√下载| 一级毛片精品| 在线免费观看的www视频| www.精华液| 久久精品国产亚洲av香蕉五月| 中文字幕最新亚洲高清| 国产一区二区三区综合在线观看| netflix在线观看网站| 久久中文字幕人妻熟女| 成年版毛片免费区| 成人手机av| 精品电影一区二区在线| 操美女的视频在线观看| 午夜a级毛片| 成人影院久久| 国产精品九九99| 欧美日韩国产mv在线观看视频| 国产有黄有色有爽视频| 国产精品一区二区精品视频观看| 午夜精品在线福利| 亚洲精品在线观看二区| 欧美另类亚洲清纯唯美| 88av欧美| 成人精品一区二区免费| 亚洲国产精品一区二区三区在线| 一二三四在线观看免费中文在| 在线观看免费日韩欧美大片| 自线自在国产av| 欧美中文日本在线观看视频| 另类亚洲欧美激情| 国产精品电影一区二区三区| 黑丝袜美女国产一区| www日本在线高清视频| 视频在线观看一区二区三区| 老熟妇仑乱视频hdxx| 免费不卡黄色视频| 一区二区日韩欧美中文字幕| 人人妻人人添人人爽欧美一区卜| 免费高清在线观看日韩| 亚洲av五月六月丁香网| 欧美乱妇无乱码| 韩国精品一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 欧美黑人精品巨大| 50天的宝宝边吃奶边哭怎么回事| 亚洲性夜色夜夜综合| 欧美性长视频在线观看| 亚洲一区二区三区欧美精品| 午夜视频精品福利| 亚洲免费av在线视频| av欧美777| 天堂影院成人在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 免费在线观看日本一区| 精品国产乱码久久久久久男人| 777久久人妻少妇嫩草av网站| 高清毛片免费观看视频网站 | 亚洲欧美一区二区三区久久| 国产成人精品久久二区二区91| 亚洲精品成人av观看孕妇| 欧美乱妇无乱码| 神马国产精品三级电影在线观看 | 黄片小视频在线播放| 欧美中文综合在线视频| 在线观看一区二区三区| 我的亚洲天堂| 久久中文字幕一级| 侵犯人妻中文字幕一二三四区| 波多野结衣一区麻豆| 国产精品影院久久| 久久久国产成人精品二区 | 亚洲成a人片在线一区二区| 久久人人精品亚洲av| 人妻久久中文字幕网| 欧美av亚洲av综合av国产av| 少妇的丰满在线观看| 午夜福利一区二区在线看| 丰满迷人的少妇在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 人妻丰满熟妇av一区二区三区| 一个人免费在线观看的高清视频| 成年版毛片免费区| 一区在线观看完整版| 韩国精品一区二区三区| 免费久久久久久久精品成人欧美视频| 国产精品一区二区三区四区久久 | 国产区一区二久久| 80岁老熟妇乱子伦牲交| 中文字幕av电影在线播放| 精品免费久久久久久久清纯| 黄网站色视频无遮挡免费观看| 视频区图区小说| 亚洲国产精品sss在线观看 | 久久亚洲精品不卡| 亚洲 欧美一区二区三区| 俄罗斯特黄特色一大片| 在线观看免费视频日本深夜| 三级毛片av免费| 男人舔女人下体高潮全视频| 欧美最黄视频在线播放免费 | 亚洲精品av麻豆狂野| 性色av乱码一区二区三区2| 精品少妇一区二区三区视频日本电影| 久99久视频精品免费| 啦啦啦免费观看视频1| 两个人看的免费小视频| 久久中文字幕人妻熟女| 中文亚洲av片在线观看爽| 国产麻豆69| 婷婷六月久久综合丁香| 啦啦啦在线免费观看视频4| 亚洲av第一区精品v没综合| 久久精品91蜜桃| 国产亚洲欧美98| 久久精品国产清高在天天线| 91在线观看av| 国产成人影院久久av| 亚洲 国产 在线| 一级片'在线观看视频| 天堂√8在线中文| 欧美国产精品va在线观看不卡| 成人手机av| 无遮挡黄片免费观看| 丰满迷人的少妇在线观看| 俄罗斯特黄特色一大片| 亚洲三区欧美一区| 欧美黄色淫秽网站| 国产精品二区激情视频| 欧美国产精品va在线观看不卡| tocl精华| 美女大奶头视频| 人人妻人人添人人爽欧美一区卜| 18禁美女被吸乳视频| 老汉色∧v一级毛片| 大码成人一级视频| 老汉色av国产亚洲站长工具| 岛国在线观看网站| 亚洲人成电影免费在线| 久热这里只有精品99| 热re99久久国产66热| 国产精品久久视频播放| 人人妻,人人澡人人爽秒播| 伦理电影免费视频| 一个人免费在线观看的高清视频| 亚洲五月色婷婷综合| 99久久精品国产亚洲精品| 18禁国产床啪视频网站| 成人永久免费在线观看视频| 午夜精品在线福利| 国产xxxxx性猛交| 国产精品电影一区二区三区| 日韩av在线大香蕉| 大型黄色视频在线免费观看| 黑人猛操日本美女一级片| 精品久久久久久久毛片微露脸| 在线观看免费视频日本深夜| 桃色一区二区三区在线观看| 色综合站精品国产| 国产精品一区二区精品视频观看| 黄色片一级片一级黄色片| 少妇被粗大的猛进出69影院| 黑人操中国人逼视频| 桃红色精品国产亚洲av| 男女午夜视频在线观看| 国产在线精品亚洲第一网站| 国产亚洲精品久久久久5区| 老汉色∧v一级毛片| 欧美另类亚洲清纯唯美| videosex国产| 欧美成人午夜精品| 国产单亲对白刺激| 中文字幕人妻丝袜一区二区| www日本在线高清视频| 日本一区二区免费在线视频| 国产高清激情床上av| 精品国产国语对白av| 在线观看免费日韩欧美大片| 在线观看一区二区三区激情| 日日爽夜夜爽网站| 中国美女看黄片| 亚洲自拍偷在线| 国产99久久九九免费精品| 69av精品久久久久久| 日韩精品青青久久久久久| 久久影院123| 中文字幕另类日韩欧美亚洲嫩草| 国产99白浆流出| 亚洲男人的天堂狠狠| 黑丝袜美女国产一区| 999久久久精品免费观看国产| 高清av免费在线| 久久精品人人爽人人爽视色| 亚洲三区欧美一区| 色综合欧美亚洲国产小说| videosex国产| 在线观看www视频免费| 欧美另类亚洲清纯唯美| 亚洲av电影在线进入| 亚洲欧美激情在线| 精品日产1卡2卡| 午夜福利在线免费观看网站| 999久久久精品免费观看国产| 欧美中文日本在线观看视频| 久久香蕉国产精品| 久久香蕉精品热| 自拍欧美九色日韩亚洲蝌蚪91| 久久久国产精品麻豆| 在线观看www视频免费| 国产成人啪精品午夜网站| 亚洲九九香蕉|