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

    利用導(dǎo)流板控制25°傾角Ahmed類車體尾流與氣動(dòng)阻力

    2014-09-27 14:18:59王漢封鄒超張運(yùn)平
    關(guān)鍵詞:尾流

    王漢封+鄒超+張運(yùn)平

    文章編號(hào):16742974(2014)04009407

    收稿日期:20130820

    基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)目(51108468)

    作者簡(jiǎn)介:王漢封(1976-),男,河南開(kāi)封人,中南大學(xué)副教授,博士

    通訊聯(lián)系人,E-mail:wanghfme@gmail.com

    摘 要:通過(guò)風(fēng)洞實(shí)驗(yàn),研究了尾部導(dǎo)流板對(duì)25°傾角Ahmed類車體尾流與氣動(dòng)阻力的影響規(guī)律.對(duì)比了斜面兩側(cè)與斜面上邊緣寬度分別為5 mm,10 mm和15 mm導(dǎo)流板的減阻效果.試驗(yàn)中模型縮尺比為1∶2,基于來(lái)流風(fēng)速與模型長(zhǎng)度的雷諾數(shù)為8.7×105.研究結(jié)果表明,模型尾流中存在一對(duì)規(guī)則的拖曳渦,并伴隨有強(qiáng)烈下掃流,尾部斜面上存在D形流動(dòng)分離區(qū).斜面兩側(cè)5 mm寬導(dǎo)流板對(duì)尾流的影響很小,對(duì)應(yīng)的氣動(dòng)阻力會(huì)增大約2.1%;斜面兩側(cè)10 mm,15 mm寬導(dǎo)流板以及不同寬度的水平導(dǎo)流板可顯著削弱尾流中的拖曳渦.水平導(dǎo)流板能夠消除斜面上的流動(dòng)再附著并破壞D形分離區(qū),其減阻效果明顯高于兩側(cè)導(dǎo)流板,最大減阻率可達(dá)11.8%.

    關(guān)鍵詞:Ahmed模型;尾流;拖曳渦;氣動(dòng)阻力;流動(dòng)控制

    中圖分類號(hào):O355;U461.1 文獻(xiàn)標(biāo)識(shí)碼:A

    Control of the Wake and Aerodynamic Drag of an Ahmed

    Model with 25° Slant Angle by Using Deflectors

    

    WANG Hanfeng1,2, ZOU Chao1, ZHANG Yunping1

    (1.School of Civil Engineering, Central South Univ, Changsha, Hunan 410075,China; 

    2. National Laboratory for HighSpeed Railway Construction,Changsha, Hunan 410075,China)

    Abstract: This paper investigated the effect of rear end deflectors on the near wake and aerodynamic drag of an Ahmed model with 25° slant angle. Drag reduction was compared for deflectors mounted on the two side faces and on the upper edge of the rear slant. The width of deflectors is 5mm, 10mm and 15 mm, respectively. The model scale is 1∶2, and the Reynolds number based oncoming flow velocity and model length is 8.7×105. The results have revealed that there is a pair of organized trailing vortices in the wake, which are accompanied by strong downwash flow. There is a Dshape flow separation zone on the slant face. Deflectors with a width of 5 mm at the two sides of the slant have negligible effect on the near wake, and slightly increase the aerodynamic drag by about 2.1%. On the other hand, all tested horizontal deflectors at the top edge of the slant and deflectors with a width of 10 mm and 15 mm at both sides of the slant considerably weaken the trailing vortices. The horizontal deflectors avoid flow reattachment on the slant and suppress the Dshape separation zone, corresponding to the maximum drag reduction rate of 11.8%, which is higher than that of the deflectors on both sides of the slant.

    Key words: Ahmed model; wakes; tailing vortex; aerodynamic drag; flow control

    

    車輛的氣動(dòng)阻力近似與其行駛速度的平方成正比.當(dāng)時(shí)速為90 km時(shí),發(fā)動(dòng)機(jī)功率的80%左右將用于克服氣動(dòng)阻力[1].通??烧J(rèn)為車輛氣動(dòng)阻力是

    由壓差阻力與摩擦阻力兩部分構(gòu)成,前者在氣動(dòng)阻力中占絕大部分[2].為提高燃油經(jīng)濟(jì)性,圍繞車輛氣動(dòng)阻力的主、被動(dòng)控制方法已開(kāi)展了廣泛的研究,如導(dǎo)流板[2-5]、漩渦發(fā)生器[6-7]、微射流[8-10]和尾部附加隔板[2, 11]等.

    實(shí)際車輛的外形復(fù)雜,不利于相關(guān)研究的對(duì)比.Ahmed模型[12]是目前研究最廣泛的類車體模型之一.該模型頭部由4個(gè)1/4圓柱面過(guò)渡,其尾部?jī)A角α可根據(jù)實(shí)際情況而選擇.研究表明,Ahmed類車體尾流及其氣動(dòng)力特性與尾部?jī)A角α有密切聯(lián)系[12-13].依據(jù)模型尾部斜面上的流動(dòng)特性可分為3個(gè)典型狀態(tài),當(dāng)α< 12.5°時(shí),流動(dòng)在尾部斜面上不會(huì)發(fā)生分離,此時(shí)尾流中會(huì)形成一對(duì)旋向相反的拖曳渦;當(dāng)12.5° <α< 30°時(shí),流動(dòng)在斜面上會(huì)發(fā)生分離與再附著,并形成分離泡,此時(shí)在模型尾流中仍會(huì)出現(xiàn)拖曳渦,但其強(qiáng)度將明顯大于第1種情況,且此時(shí)對(duì)應(yīng)的阻力系數(shù)Cd顯著增大;而在α> 30°時(shí),流動(dòng)在斜面上邊沿發(fā)生分離且無(wú)再附著發(fā)生,尾流中拖曳渦顯著減弱,模型上的壓力分布變得非常均勻,Cd顯著減小.由此可知,斜面上是否出現(xiàn)分離泡、以及拖曳渦強(qiáng)度與模型氣動(dòng)阻力有密切的聯(lián)系,拖曳渦強(qiáng)度越大對(duì)應(yīng)的模型氣動(dòng)阻力也較大[6, 12].

    湖南大學(xué)學(xué)報(bào)(自然科學(xué)版)2014年

    第4期王漢封等:利用導(dǎo)流板控制25°傾角Ahmed類車體尾流與氣動(dòng)阻力

    對(duì)于30°傾角Ahmed模型,斜面兩側(cè)導(dǎo)流板的減阻效果最為顯著,最高可達(dá)17.7%[3].而對(duì)于尾部流動(dòng)狀態(tài)完全不同的25°傾角Ahmed模型,不同位置導(dǎo)流板對(duì)氣動(dòng)阻力的影響仍缺乏系統(tǒng)的研究.依據(jù)文獻(xiàn)[3]中所提出的2種減阻效果較好的導(dǎo)流板布置方式,本文通過(guò)風(fēng)洞實(shí)驗(yàn)系統(tǒng)研究了布置于斜面兩側(cè)和斜面上邊緣的不同寬度的導(dǎo)流板對(duì)25°傾角Ahmed模型尾流與氣動(dòng)力的控制效果.實(shí)驗(yàn)運(yùn)用壓力掃描閥、眼鏡蛇探針與表面油膜流動(dòng)顯示等方法,比較了不同工況下模型氣動(dòng)阻力、尾部壓力分布以及尾流場(chǎng)的變化規(guī)律,揭示了減阻機(jī)理.

    1 試驗(yàn)方法

    1.1 風(fēng)洞模型

    本試驗(yàn)在中南大學(xué)高速鐵路建造技術(shù)國(guó)家工程實(shí)驗(yàn)室的風(fēng)洞高速試驗(yàn)段內(nèi)完成.該風(fēng)洞為回流式風(fēng)洞,具有低速與高速兩個(gè)試驗(yàn)段,其中低速試驗(yàn)段寬12 m,高3.5 m,長(zhǎng)18 m,風(fēng)速范圍為0~18 m/s,湍流度小于2%;高速試驗(yàn)段寬3 m,高3 m,長(zhǎng)15 m,風(fēng)速范圍為5~90 m/s,湍流度小于0.5%.試驗(yàn)裝置如圖1(a)所示.試驗(yàn)中Ahmed模型傾角為25°,縮尺比為1∶2,對(duì)應(yīng)的長(zhǎng)(l)、寬(w)、高(h)分別為522,194.5和144 mm.模型安裝在一個(gè)距風(fēng)洞底面約500 mm的水平板上,避免了風(fēng)洞壁面邊界層的影響.為防止流動(dòng)分離,水平板前邊緣加工成光滑的橢圓形.模型與水平板間隙為25 mm,距水平板前邊緣約550 mm.可以估算模型處平板邊界層厚度約為13.5 mm,即試驗(yàn)中模型完全處于均勻來(lái)流中.本試驗(yàn)裝置與文獻(xiàn)[3, 4, 6, 7]中所述的實(shí)驗(yàn)裝置類似.坐標(biāo)原點(diǎn)定義在水平板上模型尾部中點(diǎn)所對(duì)應(yīng)的位置上,流動(dòng)方向?yàn)楠玿,側(cè)向?yàn)楠珁,高度方向?yàn)楠珃.試驗(yàn)中自由來(lái)流風(fēng)速為U∞ = 25 m/s,對(duì)應(yīng)的基于模型長(zhǎng)度的雷諾數(shù)為8.7×105.模型所造成的風(fēng)洞阻塞率約為0.4%,其影響可忽略不計(jì).

    圖1 試驗(yàn)裝置

    Fig.1 Experimental facility

    

    本文研究了兩類不同位置導(dǎo)流板對(duì)模型尾流與氣動(dòng)阻力的影響.對(duì)無(wú)導(dǎo)流板的工況Case1也進(jìn)行了測(cè)量,以方便結(jié)果的對(duì)比.Case2,Case3和Case4中導(dǎo)流板安裝在模型尾部斜面兩側(cè),導(dǎo)流板寬度分別為5,10和15 mm,約相當(dāng)于模型長(zhǎng)度的1%,2%和3%;Case5,Case6和Case7中導(dǎo)流板安裝于斜面上邊緣處,寬度分別為5,10和15 mm,如圖1(b)所示.

    1.2 測(cè)試方法

    試驗(yàn)中采用眼鏡蛇探針測(cè)量模型尾流中的總壓與速度分布,所用探針響應(yīng)頻率為2.5 kHz,并已成功運(yùn)用于多種湍流場(chǎng)的測(cè)量[14-15].實(shí)驗(yàn)中探針采樣頻率為2 kHz,每一測(cè)點(diǎn)采樣時(shí)間為15 s.測(cè)量分別在模型下游0.5l與l的流向截面內(nèi)進(jìn)行,以觀察尾流中的拖曳渦結(jié)構(gòu).探針固定于計(jì)算機(jī)控制的二維移測(cè)架上在測(cè)量平面內(nèi)逐點(diǎn)進(jìn)行測(cè)量.移測(cè)架位移精度為0.02 mm.考慮到Ahmed模型尾流的對(duì)稱性,測(cè)量?jī)H在y > 0的范圍內(nèi)進(jìn)行.

    為研究不同工況下模型尾部壓力的變化情況,運(yùn)用電子壓力掃描閥對(duì)模型尾部斜面與垂面上壓力分布進(jìn)行了測(cè)量.壓力測(cè)點(diǎn)的布置與文獻(xiàn)[7, 9]相同.試驗(yàn)中每測(cè)點(diǎn)掃描12 000次,以獲得各點(diǎn)平均壓力系數(shù)Cp,其定義式為Cp = -P

    SymboleB@

    /0.5ρU2

    SymboleB@

    ,其中為各測(cè)點(diǎn)平均壓力,P

    SymboleB@

    為風(fēng)洞靜壓力,ρ為空氣密度.本文中上橫線“ˉ”表示時(shí)間平均量.

    試驗(yàn)還采用了表面油膜法對(duì)模型尾部斜面上的流動(dòng)分離情況進(jìn)行了研究.用二甲基硅油、煤油和鈦白粉按一定比例混合拌勻[16-17],并均勻地涂抹在模型尾部斜面上.在25 m/s風(fēng)速下,約10 min,油膜可達(dá)到穩(wěn)定狀態(tài).

    2 結(jié)果及分析

    2.1 流場(chǎng)測(cè)量結(jié)果

    Case4和Case7的對(duì)應(yīng)結(jié)果分別與Case3和Case6非常類似,限于篇幅,它們?cè)诤罄m(xù)討論中未予給出.圖2給出了當(dāng)x=0.5l和l時(shí),各工況以時(shí)均流向渦量ω*x為背景的流線圖.本文中上標(biāo)“*”表示用U∞與l進(jìn)行無(wú)量綱化.由圖2可知,各工況模型尾流中均存在一對(duì)規(guī)則的流向拖曳渦(y < 0沒(méi)有顯示),并總是伴隨著尾流中心線附近的強(qiáng)烈下掃流,這與文獻(xiàn)[3, 4, 6, 18]等報(bào)道的結(jié)果是一致的.

    圖2 x=0.5l和x=l截面內(nèi)的時(shí)均渦量與流線圖

    Fig.2 Timeaveraged vorticity and stremlines in the streamwise planes at x=0.5l and l

    

    在x=0.5l處,Case1的拖曳渦中心ω*x最大值約為13.0.對(duì)于Case2,拖曳渦強(qiáng)度相對(duì)于Case1無(wú)明顯變化.而對(duì)于Case 3, ω*x的最大值僅為7.6,相對(duì)于Case1減弱了約41.5%,且拖曳渦的尺寸也有明顯減小.而對(duì)于導(dǎo)流板水平布置在斜面上邊緣的2種情況Case5和Case6,ω*x最大值分別為7.5和7.4,相對(duì)于Case1分別減弱了42.3%和43.1%.如圖2所示各工況的流線圖也可反映拖曳渦的結(jié)構(gòu)與強(qiáng)度.對(duì)于Case1和Case2,流線在拖曳渦范圍內(nèi)存在強(qiáng)烈的螺旋結(jié)構(gòu);而在Case3,Case5和Case6中,拖曳渦中心附近流線的螺旋結(jié)構(gòu)相對(duì)較弱.與x=0.5l的情況類似,在x=l截面內(nèi),Case1與Case2的拖曳渦強(qiáng)度基本相同,而Case3,Case5和Case6的對(duì)應(yīng)值則明顯較小.5種工況在x=l截面內(nèi)ω*x的最大值分別是0.5l截面內(nèi)對(duì)應(yīng)值的61%,66%,95%, 55%和54%,這說(shuō)明Case3中拖曳渦衰減速率最慢,而Case5與Case6的衰減速率相對(duì)較快.從圖2還可看出,Case1與Case2對(duì)應(yīng)的拖曳渦中心位置也基本相同,而Case3,Case5和Case6的渦團(tuán)中心更靠近尾流中心線,這表明后3種工況下掃流向外側(cè)排開(kāi)拖曳渦的作用相對(duì)較弱.

    為定量比較尾流中下掃流的變化規(guī)律,圖3給出了x=0.5l和l截面內(nèi)的z方向時(shí)均速度W*的分布.各工況拖曳渦中心位置在圖3中用“×”標(biāo)出,以方便對(duì)比.由圖3可知,在各工況下,拖曳渦中心內(nèi)側(cè)均存在著明顯的下掃流,即W*< 0.Case1和Case2對(duì)應(yīng)的W*定性與定量上都非常類似,而Case3相對(duì)于Case1也僅略有減小,這表明斜面兩側(cè)導(dǎo)流板對(duì)下掃流的影響非常有限.相對(duì)于Case1,Case5和Case6中下掃流的強(qiáng)度和其影響范圍都明顯減小了,這與圖2中拖曳渦的變化規(guī)律是一致的.總體來(lái)看,模型尾部的水平導(dǎo)流板對(duì)拖曳渦和下掃流的抑制作用更為顯著.

    圖3 x=0.5l和x=l截面內(nèi)的時(shí)均z方向速度

    Fig.3 Timeaveraged velocity in z direction in the streamwise planes at x=0.5l and l

    

    圖4給出了各工況拖曳渦中心處y方向速度v的能譜.Case1的能譜Ev存在顯著的峰值,其對(duì)應(yīng)的基于l和U

    SymboleB@

    的斯托羅哈數(shù)St=1.55,與文獻(xiàn)[13]的結(jié)果非常吻合.這表明Case1中拖曳渦強(qiáng)度和周期性均較顯著.Case2中,盡管能譜峰值略有減小,但其St數(shù)與Case1相同.與前2種工況不同,Case3,Case5和Case6的能譜中已沒(méi)有明顯峰值出現(xiàn),表明這些工況中拖曳渦已無(wú)顯著的周期性.在x=l截面內(nèi),Ev所表現(xiàn)的規(guī)律與x=0.5l截面內(nèi)完全一致.

    圖4 拖曳渦中心處速度v的能譜

    Fig.4Power spectra density function of v at tailing vortex center

    

    圖5給出用表面油膜法獲得的模型尾部斜面上的流動(dòng)結(jié)構(gòu).由于模型尾流的對(duì)稱性,圖中僅給出了表面油膜流動(dòng)顯示的右半部分,而在左半部分給出了相應(yīng)的流動(dòng)示意圖(Case6大部分區(qū)域流動(dòng)結(jié)構(gòu)已不明顯,故未給出).Case1,Case2和Case3中流動(dòng)分離并非發(fā)生在斜面上邊緣,而是上邊緣略下游

    的實(shí)線所示位置上,如圖5所示.Case1中,上邊緣附近的分離流在斜面上發(fā)生再附,并在斜面上形成一個(gè)D形分離泡,如圖5中流動(dòng)分離線與虛線所圍成范圍,這與文獻(xiàn)[9, 12, 19, 20]中的結(jié)果是完全一致的.Case2中,尾部斜面上的流動(dòng)結(jié)構(gòu)沒(méi)有明顯改變,仍可清晰地觀察到D形流動(dòng)分離區(qū).Case3和Case1相比,D形分離區(qū)仍然存在,但略有減小.對(duì)于Case5,斜面上流動(dòng)結(jié)構(gòu)相對(duì)于Case1發(fā)生了顯著的變化.在水平導(dǎo)流板的作用下,斜面上邊緣附近的流動(dòng)分離線消失了,且分離流在斜面上不會(huì)發(fā)生再附,因而斜面上不再出現(xiàn)封閉的分離泡.Case5所對(duì)應(yīng)的流動(dòng)狀態(tài),非常類似于文獻(xiàn)[13,16,21]中所

    圖5 模型尾部斜面上流動(dòng)顯示結(jié)果

    Fig.5 Surface flow pattern on the slant face

    

    給出的30°或35°傾角Ahmed模型的尾流結(jié)構(gòu).隨著水平導(dǎo)流板的寬度增加到10 mm,Case6中斜面兩側(cè)分離流的影響也基本消失,除斜面左右兩個(gè)角部區(qū)外,整個(gè)斜面基本上都處于分離區(qū)內(nèi),斜面上流動(dòng)較為均勻.

    綜合圖2~圖5可知,對(duì)于斜面兩側(cè)導(dǎo)流板的情況,其尾流特性與Case1是類似的.隨著導(dǎo)流板寬度的增加,斜面上D形分離泡逐漸減小,尾流中拖曳渦強(qiáng)度有所減弱.而斜面上邊緣的水平導(dǎo)流板,可破壞斜面上的D形分離泡,并能夠更為顯著地抑制尾流拖曳渦強(qiáng)度,其作用類似于增大了25°Ahmed模型的尾部?jī)A角.

    2.2 氣動(dòng)阻力試驗(yàn)結(jié)果

    2.2.1 氣動(dòng)阻力

    定義一個(gè)包括模型在內(nèi)的控制體積(如圖6所示)[22],并將動(dòng)量守恒方程應(yīng)用于該控制體積,可以獲得模型氣動(dòng)阻力的精確表達(dá)式[4].當(dāng)控制體足夠大時(shí),可認(rèn)為控制體側(cè)面與頂面上沒(méi)有動(dòng)量輸運(yùn).此外,雷諾應(yīng)力、氣體粘性力等對(duì)氣動(dòng)阻力的貢獻(xiàn)比其他項(xiàng)小一個(gè)數(shù)量級(jí)以上,通常也可忽略[4, 22].模型氣動(dòng)阻力表達(dá)式可簡(jiǎn)化為:

    Fx=-12ρU2

    SymboleB@

    ∫S1-xU

    SymboleB@

    2dS+

    12ρU2

    SymboleB@

    ∫S2yU2

    SymboleB@

    +2zU2

    SymboleB@

    dS+∫S(Pi0-Pi)dS.(1)

    式中:Pi0為來(lái)流總壓;Pi為控制體積出口截面上各點(diǎn)總壓;x,y與z分別為出口截面上3個(gè)方向速度的時(shí)均值;S為控制體的出口面積.式(1)右側(cè)三項(xiàng)分別表示流向速度損失、側(cè)向速度變化以及總壓損失對(duì)氣動(dòng)阻力的貢獻(xiàn).已有文獻(xiàn)[4, 22, 23]成功運(yùn)用式(1)獲得了類車體的氣動(dòng)阻力,本文也將采用此方法估算不同工況下模型氣動(dòng)阻力.

    表1給出了基于x=0.5l和l截面測(cè)量結(jié)果,依據(jù)式(1)估算的模型氣動(dòng)阻力.可以看出,由上述兩個(gè)截面測(cè)量數(shù)據(jù)所得到的氣動(dòng)阻力是非常接近

    的.對(duì)比式(1)右側(cè)三項(xiàng)對(duì)氣動(dòng)阻力的貢獻(xiàn)可發(fā)現(xiàn),尾流中的總壓損失占?xì)鈩?dòng)阻力的絕大部分,而流向速度損失與側(cè)向速度變化對(duì)阻力的貢獻(xiàn)則相對(duì)較小.Case1中模型阻力系數(shù)Cd=0.432,與文獻(xiàn)[19, 24]的結(jié)果非常接近,這也驗(yàn)證了本試驗(yàn)結(jié)果的可靠性.由表1可知,Case2中Cd=0.441,與Case1非常接近,相對(duì)于Case1略微增大約2.1%.這一結(jié)果與3.1節(jié)所述流場(chǎng)變化規(guī)律是吻合的.Case3中Cd=0.415,Case4中Cd=0.399,相對(duì)于Case1的減阻率分別為3.9%和7.6%.而斜面上邊緣水平導(dǎo)流板工況Case5,Case6和Case7,對(duì)應(yīng)的減阻率可達(dá)10.9%,11.6%和11.8%,減阻效果十分顯著,明顯優(yōu)于斜面兩側(cè)導(dǎo)流板各工況.圖7給出了減阻率隨導(dǎo)流板寬度的變化情況.對(duì)于水平導(dǎo)流板,減阻率隨導(dǎo)流板寬度的增加變化很?。欢鴮?duì)于斜面兩側(cè)導(dǎo)流板,減阻率隨著導(dǎo)流板寬度的增加而逐漸增大,但始終低于水平導(dǎo)流板的減阻率.結(jié)合3.1節(jié)中流場(chǎng)測(cè)量

    結(jié)果可知,Cd的減小與尾流中拖曳渦強(qiáng)度的減弱是相關(guān)的.總體來(lái)看,斜面上邊緣導(dǎo)流板對(duì)尾流拖曳渦與氣動(dòng)阻力的抑制作用明顯強(qiáng)于斜面兩側(cè)導(dǎo)流板,這與30°傾角Ahmed模型的對(duì)應(yīng)規(guī)律[3]是截然不同的.

    圖6 以動(dòng)量守恒法計(jì)算模型氣動(dòng)阻力的控制體積[22]

    Fig.6 Control volume for drag estimation using

    momentum conservation[22]

    

    表1 由式(1)計(jì)算的氣動(dòng)阻力

    Tab.1 Aerodynamic drag estimated based on Eq(1)

    工況

    測(cè)量截面

    第1項(xiàng)

    阻力/N

    第2項(xiàng)

    阻力/N

    第3項(xiàng)

    阻力/N

    總阻力/N

    總阻力均值 /N

    阻力系數(shù)

    減阻效果/%

    Case1

    x=0.5l

    -0.972

    1.368

    4.142

    4.538

    x=l

    -0.684

    1.172

    4.040

    4.528

    4.533

    0.432

    -

    Case2

    x=0.5l

    -0.814

    0.988

    4.502

    4.676

    x=l

    -0.612

    1.128

    4.090

    4.606

    4.641

    0.441

    2.1

    Case3

    x=0.5l

    -1.462

    0.548

    5.278

    4.364

    x=l

    -0.786

    0.508

    4.636

    4.358

    4.361

    0.415

    -3.9

    Case4

    x=0.5l

    -1.392

    0.446

    5.214

    4.268

    x=l

    -0.878

    0.454

    4.544

    4.120

    4.194

    0.399

    -7.6

    Case5

    x=0.5l

    -1.702

    0.526

    5.244

    4.068

    x=l

    -1.150

    0.272

    4.928

    4.050

    4.059

    0.385

    -10.9

    Case6

    x=0.5l

    -1.474

    1.032

    4.522

    4.080

    x=l

    -1.094

    0.410

    4.618

    3.934

    4.007

    0.382

    -11.6

    Case7

    x=0.5l

    -1.426

    1.096

    4.342

    4.012

    x=l

    -0.980

    0.406

    4.558

    3.984

    3.998

    0.381

    -11.8

    圖7 各工況的減阻率

    Fig.7 Drag reduction rate

    2.2.2 尾部壓力分布

    圖8給出了5種工況中模型尾部斜面與垂面上的壓力分布.總的來(lái)看,尾部垂面壓力分布受導(dǎo)流板的影響較小,壓力系數(shù)Cp均為-0.25~-0.35.然而,尾部斜面的壓力分布與導(dǎo)流板位置及導(dǎo)流板寬度密切相關(guān).Case1中斜面的上邊緣與右邊緣附近均出現(xiàn)了較強(qiáng)的負(fù)壓,與文獻(xiàn)[7, 9]的測(cè)量結(jié)果是

    一致的.這說(shuō)明Case1中,斜面上邊緣與側(cè)邊緣均存在較強(qiáng)的流動(dòng)分離.Case2中斜面上壓力分布無(wú)明顯變化,僅上邊緣附近的負(fù)壓極值略有增大.這與表1所示Case2中氣動(dòng)阻力的變化規(guī)律是吻合的.Case3的斜面壓力分布與前2種工況截然不同,兩側(cè)與上邊緣附近的負(fù)壓極值明顯減小,Cp的極小值為-0.5左右,僅相當(dāng)于Case1對(duì)應(yīng)值的一半.而對(duì)于水平導(dǎo)流板的工況Case5和Case6,其上邊緣與側(cè)邊緣附近的壓力極值消失,整個(gè)斜面的壓力分布變得非常均勻,且斜面上的負(fù)壓明顯減弱了.這表明Case5和Case6中模型尾部斜面上邊緣與側(cè)邊緣附近的流動(dòng)分離被顯著削弱.Case4和Case7的壓力分布情況分別與Case3和Case6非常類似,圖8中未給出.綜上所述,當(dāng)導(dǎo)流板寬度分別為10和15 mm時(shí),無(wú)論是布置在斜面兩側(cè)還是水平布置在斜面上邊緣處,均能起到減小模型氣動(dòng)阻力的作用,但水平布置在斜面邊上緣處的導(dǎo)流板的減阻效果更優(yōu).

    

    圖8模型尾部斜面和垂面壓力分布

    Fig.8 Pressure distributions on the slant and rear surfaces of the mode

    

    Ahmed模型氣動(dòng)阻力主要由頭部的正壓、尾部斜面與垂面上的負(fù)壓、以及其他各表面的摩擦阻力構(gòu)成,其中尾部斜面和垂面上的負(fù)壓占據(jù)了總氣動(dòng)阻力的絕大部分[2].對(duì)于傾角為25°的Ahmed模型,當(dāng)雷諾數(shù)為7.0×105時(shí),尾部斜面和垂面負(fù)壓對(duì)總氣動(dòng)阻力的貢獻(xiàn)約為80%[25].將斜面與垂面的壓力投影至x方向并積分,可獲得各面對(duì)應(yīng)的阻力系數(shù),如表2所示.對(duì)于Case1,斜面與垂面阻力系數(shù)分別為0.187和0.158,可估算此時(shí)對(duì)應(yīng)的Cd約為0.431,這與表1所示結(jié)果非常吻合.對(duì)比表2中數(shù)據(jù)可知,Case2中,尾部垂面阻力系數(shù)基本沒(méi)有變化.造成Case2中Cd增大的原因主要是斜面阻力系數(shù)增大了,因?yàn)槿鐖D8所示斜面上邊緣附近的負(fù)壓相對(duì)于Case1有所增強(qiáng).對(duì)于Case3~Case7,其尾

    表2 尾部斜面與垂面阻力系數(shù)

    Tab.2 Drag coefficients of slant and rear surfaces

    工況

    尾部斜面

    阻力系數(shù)

    尾部垂面

    阻力系數(shù)

    尾部總

    阻力系數(shù)

    模型總

    阻力系數(shù)

    Case1

    0.187

    0.158

    0.345

    0.431

    Case2

    0.199

    0.159

    0.358

    0.444

    Case3

    0.139

    0.191

    0.330

    0.416

    Case4

    0.121

    0.193

    0.314

    0.400

    Case5

    0.121

    0.184

    0.305

    0.391

    Case6

    0.110

    0.189

    0.299

    0.385

    Case7

    0.109

    0.189

    0.298

    0.384

    部垂面阻力系數(shù)不僅沒(méi)有減小,相對(duì)于Case1反而略有增大,由此可知,各工況減阻效果主要來(lái)源于斜面上負(fù)壓的減弱.即如圖8所示,導(dǎo)流板顯著抑制了斜面上邊緣與側(cè)邊緣附近強(qiáng)烈的流動(dòng)分離,并削弱了斜面上負(fù)壓.

    3 結(jié) 論

    通過(guò)風(fēng)洞試驗(yàn)研究了25°傾角Ahmed類車體尾部斜面兩側(cè)和斜面上邊緣處不同寬度導(dǎo)流板對(duì)模型尾流與氣動(dòng)力的影響規(guī)律.導(dǎo)流板寬度分別為5,10和15 mm,約相當(dāng)于車長(zhǎng)的1%,2%和3%.主要結(jié)論如下:

    1) 當(dāng)模型尾部斜面兩側(cè)導(dǎo)流板寬為5mm時(shí)(Case2),其對(duì)拖曳渦與下掃流的影響可以忽略.斜面兩側(cè)寬分別為10,15 mm導(dǎo)流板(Case3,Case4)和上邊緣寬分別為5,10,15 mm的水平導(dǎo)流板(Case5~Case7)均能夠明顯削弱尾流中拖曳渦與下掃流強(qiáng)度.

    2) 無(wú)導(dǎo)流板時(shí),模型尾部斜面上邊緣附近分離流會(huì)發(fā)生再附著并形成D形流動(dòng)分離泡.對(duì)于斜面兩側(cè)導(dǎo)流板,分離泡隨導(dǎo)流板寬度的增加而有所減小,但不會(huì)消失;而斜面上邊緣導(dǎo)流板能夠抑制斜面上的流動(dòng)再附著,破壞分離泡的形成.隨著上邊緣導(dǎo)流板寬度的增加,斜面上流動(dòng)變得非常均勻,類似于增大了Ahmed模型的尾部?jī)A角.

    3) 基于尾流中總壓與時(shí)均速度的測(cè)量結(jié)果,依據(jù)Onorato等[22]給出的方法估算氣動(dòng)阻力是可行的.斜面兩側(cè)布置5 mm寬導(dǎo)流板(Case2)不僅無(wú)減阻效果,反而使氣動(dòng)阻力增大約2.1%.當(dāng)兩側(cè)導(dǎo)流板寬度分別增加到10和15 mm時(shí)(Case3,Case4),減阻效率分布為3.9%和7.6%.Case5~Case7的減阻效率分別可達(dá)10.9%,11.6%和11.8%,基本不隨導(dǎo)流板寬度而變化.斜面兩側(cè)導(dǎo)流板的減阻效率隨導(dǎo)流板寬度增加而逐漸增大,但始終小于水平導(dǎo)流板的減阻效率.

    參考文獻(xiàn)

    [1] 谷正氣.汽車空氣動(dòng)力學(xué)[M].北京:人民交通出社,2005:10-15.

    GU Zhengqi. Automobile aerodynamics[M].Beijing: China Communications Press,2005:10-15.

    [2] 王勛年,李士偉,陳立.采用尾部隔板降低類客車體阻力的研究[J].實(shí)驗(yàn)流體力學(xué),2011, 25(2): 58-62.

    WANG Xunnian, LI Shiwei, CHEN Li. Investigation on reducing drag of the buslike body using tail clapboards[J]. Journal of Experiments in Fluid Mechanics,2011,25(2):58-62.(In Chinese)

    [3] BEAUDOIN J, AIDER J. Drag and lift reduction of a 3D bluff body using flaps[J].Exp Fluids, 2008, 44:491-501.

    [4] FOURRIE G, KEIRSBULCK L, LABRAGA L,et al. Bluffbody drag reduction using a deflector[J]. Exp Fluids, 2011, 50: 385-395.

    [5] 吳志剛,魏琪,加藤征三,等. 導(dǎo)流翼片的傾角和長(zhǎng)度在降低大后壁車輛氣動(dòng)阻力中的應(yīng)用[J]. 汽車工程,2003,25(6): 634-637.

    WU Zhigang, WEI Qi, SEIZO Kato, et al. The effect of length and inclined angle of air deflectors on reducing drag of large bluff end vehicles[J]. Automotive Engineering,2003,25(6): 634-637.(In Chinese)

    [6] AIDERD J, BEAUDOIN J, WESFREID J. Drag and lift reduction of a 3D bluff body using active vortex generators[J]. Exp Fluids, 2010, 48:771-789.

    [7] PUJALS G, DEPARDON S, COSSU C. Drag reduction of a 3D bluff body using coherent streamwise streaks[J]. Exp Fluids, 2010, 49: 1085-1094.

    [8] 谷正氣,李學(xué)武,何憶斌. 汽車減阻新方法[J]. 汽車工程,2008,30(5):441-443.

    GU Zhengqi LI Xuewu HE Yibin. A new method of reducing aerodynamic drag[J]. Automotive Engineering,2008,30(5):441-443.(In Chinese)

    [9] JOSEPH P, AMANDOLESE X,AIDER J L.Drag reduction on the 25°slant angle Ahmed reference body using pulsed jets[J]. Exp Fluids, 2011, 52:1169-1185.

    [10]AUBRUN S, McNally J, ALVI F,et al. Separation flow control on a generic ground vehicle using steady microjet arrays[J]. Exp Fluids,2011, 51: 1177-1187.

    [11]張攀峰,王晉軍,唐青. 氣動(dòng)附加裝置降低廂式貨車后體阻力[J]. 實(shí)驗(yàn)流體力學(xué),2009,23:12-15.

    ZHANG Panfeng,WANG Jinjun, TANG Qing. Experimental investigation on the aft body drag reduction of the tractortrailer truck by aerodynamic addon device[J].Journal of Experiments in Fluid Mechanics,2009,23:12-15.(In Chinese)

    [12]AHMED S R, RAMM R, FALTIN G. Some salient features ofthe timeaveraged ground vehicle wake[C]// SAE Technical Paper Series 840300,1984.

    [13]VINO G, WATKINS S, MOUSLEY P,et al. Flow structures in the near wake of Ahmed model[J]. J Fluids Struct,2005, 20: 673-695.

    [14]SCHNEIDER G, HOOPER J, MUSGROVE A,et al. Velocity and Reynolds stresses in a precessing jet flow [J]. Exp Fluids, 1997, 22: 489-495.

    [15]CHEN J, HAYNES B, FLETCHER D. Cobra probe measurements of mean velocities, reynolds stresses and highorder velocity correlations in pipe flow [J]. Exp Therm Fluid Sci, 2000, 21: 206-217.

    [16]CONAN B, ANTHOINE J, PLANQUART P. Experimental aerodyanmic study of a cartype bluff body[J]. Exp Fluids, 2011, 50: 1273-1284. 

    [17]THACKER A, AUBRUN S, LEROY A,et al. Effects of suppressing the 3D separation on the rear slant on the flow structures around an Ahmed body[J]. J Wind Eng Ind Aerodyn,2012,107/108:237-243.

    [18]朱暉,楊志剛. 類車體尾跡區(qū)流動(dòng)的實(shí)驗(yàn)研究[J]. 實(shí)驗(yàn)流體力學(xué),2010,24(2):24-27.

    ZHU Hui, YANG Zhigang. Experimental study on the flow field in the wake of Ahmed model[J]. Journal of Experiments in Fluid Mechanics,2010,24(2):24-27.(In Chinese)

    [19]GILLIERON P, CHOMETON F. Modelling of stationary threedimensional separated air flows around an Ahmed reference model[C]//3rd International Workshop on Vortex, ESAIM Proceedings. Amsterdam:Elsevier, 1999,7:173-182.

    [20]KRAJNOVIC S, DAVIDSON L. Flow around a simplified car, part 2: understanding the flow[J]. J Fluids Eng,2005, 127: 919-928.

    [21]GILLIERON P, KOURTA A. Aerodynamic drag control by pulsed jets on simplified car geometry[J]. Exp Fluids, 2013, 54:1-16.

    [22]ONORATO M, COSTELLI A F, GARONNE A. Drag measurement through wake analysis[C]//SAE Technical Paper Series 840302,1984.

    [23]VANDAM C. Recent experience with different methods of drag prediction[J]. Prog Aerosp Sci,1999, 35(8):751-798.

    [24]BRUNN A, WASSEN E, SPERBER D, et al. Active drag control for a generic car model [C]// Notes on Numerical Fluid Mechanics and Multidisciplinary Design.Berlin Heidelberg:Springer, 2007: 247-259. 

    [25]KRAJNOVIC S, BASARA B. LES of the flow around Ahmed body with active flow control [C]//Turbulence and Interactions. Berlin Heidelberg: Springer, 2010: 247-254.

    [5] 吳志剛,魏琪,加藤征三,等. 導(dǎo)流翼片的傾角和長(zhǎng)度在降低大后壁車輛氣動(dòng)阻力中的應(yīng)用[J]. 汽車工程,2003,25(6): 634-637.

    WU Zhigang, WEI Qi, SEIZO Kato, et al. The effect of length and inclined angle of air deflectors on reducing drag of large bluff end vehicles[J]. Automotive Engineering,2003,25(6): 634-637.(In Chinese)

    [6] AIDERD J, BEAUDOIN J, WESFREID J. Drag and lift reduction of a 3D bluff body using active vortex generators[J]. Exp Fluids, 2010, 48:771-789.

    [7] PUJALS G, DEPARDON S, COSSU C. Drag reduction of a 3D bluff body using coherent streamwise streaks[J]. Exp Fluids, 2010, 49: 1085-1094.

    [8] 谷正氣,李學(xué)武,何憶斌. 汽車減阻新方法[J]. 汽車工程,2008,30(5):441-443.

    GU Zhengqi LI Xuewu HE Yibin. A new method of reducing aerodynamic drag[J]. Automotive Engineering,2008,30(5):441-443.(In Chinese)

    [9] JOSEPH P, AMANDOLESE X,AIDER J L.Drag reduction on the 25°slant angle Ahmed reference body using pulsed jets[J]. Exp Fluids, 2011, 52:1169-1185.

    [10]AUBRUN S, McNally J, ALVI F,et al. Separation flow control on a generic ground vehicle using steady microjet arrays[J]. Exp Fluids,2011, 51: 1177-1187.

    [11]張攀峰,王晉軍,唐青. 氣動(dòng)附加裝置降低廂式貨車后體阻力[J]. 實(shí)驗(yàn)流體力學(xué),2009,23:12-15.

    ZHANG Panfeng,WANG Jinjun, TANG Qing. Experimental investigation on the aft body drag reduction of the tractortrailer truck by aerodynamic addon device[J].Journal of Experiments in Fluid Mechanics,2009,23:12-15.(In Chinese)

    [12]AHMED S R, RAMM R, FALTIN G. Some salient features ofthe timeaveraged ground vehicle wake[C]// SAE Technical Paper Series 840300,1984.

    [13]VINO G, WATKINS S, MOUSLEY P,et al. Flow structures in the near wake of Ahmed model[J]. J Fluids Struct,2005, 20: 673-695.

    [14]SCHNEIDER G, HOOPER J, MUSGROVE A,et al. Velocity and Reynolds stresses in a precessing jet flow [J]. Exp Fluids, 1997, 22: 489-495.

    [15]CHEN J, HAYNES B, FLETCHER D. Cobra probe measurements of mean velocities, reynolds stresses and highorder velocity correlations in pipe flow [J]. Exp Therm Fluid Sci, 2000, 21: 206-217.

    [16]CONAN B, ANTHOINE J, PLANQUART P. Experimental aerodyanmic study of a cartype bluff body[J]. Exp Fluids, 2011, 50: 1273-1284. 

    [17]THACKER A, AUBRUN S, LEROY A,et al. Effects of suppressing the 3D separation on the rear slant on the flow structures around an Ahmed body[J]. J Wind Eng Ind Aerodyn,2012,107/108:237-243.

    [18]朱暉,楊志剛. 類車體尾跡區(qū)流動(dòng)的實(shí)驗(yàn)研究[J]. 實(shí)驗(yàn)流體力學(xué),2010,24(2):24-27.

    ZHU Hui, YANG Zhigang. Experimental study on the flow field in the wake of Ahmed model[J]. Journal of Experiments in Fluid Mechanics,2010,24(2):24-27.(In Chinese)

    [19]GILLIERON P, CHOMETON F. Modelling of stationary threedimensional separated air flows around an Ahmed reference model[C]//3rd International Workshop on Vortex, ESAIM Proceedings. Amsterdam:Elsevier, 1999,7:173-182.

    [20]KRAJNOVIC S, DAVIDSON L. Flow around a simplified car, part 2: understanding the flow[J]. J Fluids Eng,2005, 127: 919-928.

    [21]GILLIERON P, KOURTA A. Aerodynamic drag control by pulsed jets on simplified car geometry[J]. Exp Fluids, 2013, 54:1-16.

    [22]ONORATO M, COSTELLI A F, GARONNE A. Drag measurement through wake analysis[C]//SAE Technical Paper Series 840302,1984.

    [23]VANDAM C. Recent experience with different methods of drag prediction[J]. Prog Aerosp Sci,1999, 35(8):751-798.

    [24]BRUNN A, WASSEN E, SPERBER D, et al. Active drag control for a generic car model [C]// Notes on Numerical Fluid Mechanics and Multidisciplinary Design.Berlin Heidelberg:Springer, 2007: 247-259. 

    [25]KRAJNOVIC S, BASARA B. LES of the flow around Ahmed body with active flow control [C]//Turbulence and Interactions. Berlin Heidelberg: Springer, 2010: 247-254.

    [5] 吳志剛,魏琪,加藤征三,等. 導(dǎo)流翼片的傾角和長(zhǎng)度在降低大后壁車輛氣動(dòng)阻力中的應(yīng)用[J]. 汽車工程,2003,25(6): 634-637.

    WU Zhigang, WEI Qi, SEIZO Kato, et al. The effect of length and inclined angle of air deflectors on reducing drag of large bluff end vehicles[J]. Automotive Engineering,2003,25(6): 634-637.(In Chinese)

    [6] AIDERD J, BEAUDOIN J, WESFREID J. Drag and lift reduction of a 3D bluff body using active vortex generators[J]. Exp Fluids, 2010, 48:771-789.

    [7] PUJALS G, DEPARDON S, COSSU C. Drag reduction of a 3D bluff body using coherent streamwise streaks[J]. Exp Fluids, 2010, 49: 1085-1094.

    [8] 谷正氣,李學(xué)武,何憶斌. 汽車減阻新方法[J]. 汽車工程,2008,30(5):441-443.

    GU Zhengqi LI Xuewu HE Yibin. A new method of reducing aerodynamic drag[J]. Automotive Engineering,2008,30(5):441-443.(In Chinese)

    [9] JOSEPH P, AMANDOLESE X,AIDER J L.Drag reduction on the 25°slant angle Ahmed reference body using pulsed jets[J]. Exp Fluids, 2011, 52:1169-1185.

    [10]AUBRUN S, McNally J, ALVI F,et al. Separation flow control on a generic ground vehicle using steady microjet arrays[J]. Exp Fluids,2011, 51: 1177-1187.

    [11]張攀峰,王晉軍,唐青. 氣動(dòng)附加裝置降低廂式貨車后體阻力[J]. 實(shí)驗(yàn)流體力學(xué),2009,23:12-15.

    ZHANG Panfeng,WANG Jinjun, TANG Qing. Experimental investigation on the aft body drag reduction of the tractortrailer truck by aerodynamic addon device[J].Journal of Experiments in Fluid Mechanics,2009,23:12-15.(In Chinese)

    [12]AHMED S R, RAMM R, FALTIN G. Some salient features ofthe timeaveraged ground vehicle wake[C]// SAE Technical Paper Series 840300,1984.

    [13]VINO G, WATKINS S, MOUSLEY P,et al. Flow structures in the near wake of Ahmed model[J]. J Fluids Struct,2005, 20: 673-695.

    [14]SCHNEIDER G, HOOPER J, MUSGROVE A,et al. Velocity and Reynolds stresses in a precessing jet flow [J]. Exp Fluids, 1997, 22: 489-495.

    [15]CHEN J, HAYNES B, FLETCHER D. Cobra probe measurements of mean velocities, reynolds stresses and highorder velocity correlations in pipe flow [J]. Exp Therm Fluid Sci, 2000, 21: 206-217.

    [16]CONAN B, ANTHOINE J, PLANQUART P. Experimental aerodyanmic study of a cartype bluff body[J]. Exp Fluids, 2011, 50: 1273-1284. 

    [17]THACKER A, AUBRUN S, LEROY A,et al. Effects of suppressing the 3D separation on the rear slant on the flow structures around an Ahmed body[J]. J Wind Eng Ind Aerodyn,2012,107/108:237-243.

    [18]朱暉,楊志剛. 類車體尾跡區(qū)流動(dòng)的實(shí)驗(yàn)研究[J]. 實(shí)驗(yàn)流體力學(xué),2010,24(2):24-27.

    ZHU Hui, YANG Zhigang. Experimental study on the flow field in the wake of Ahmed model[J]. Journal of Experiments in Fluid Mechanics,2010,24(2):24-27.(In Chinese)

    [19]GILLIERON P, CHOMETON F. Modelling of stationary threedimensional separated air flows around an Ahmed reference model[C]//3rd International Workshop on Vortex, ESAIM Proceedings. Amsterdam:Elsevier, 1999,7:173-182.

    [20]KRAJNOVIC S, DAVIDSON L. Flow around a simplified car, part 2: understanding the flow[J]. J Fluids Eng,2005, 127: 919-928.

    [21]GILLIERON P, KOURTA A. Aerodynamic drag control by pulsed jets on simplified car geometry[J]. Exp Fluids, 2013, 54:1-16.

    [22]ONORATO M, COSTELLI A F, GARONNE A. Drag measurement through wake analysis[C]//SAE Technical Paper Series 840302,1984.

    [23]VANDAM C. Recent experience with different methods of drag prediction[J]. Prog Aerosp Sci,1999, 35(8):751-798.

    [24]BRUNN A, WASSEN E, SPERBER D, et al. Active drag control for a generic car model [C]// Notes on Numerical Fluid Mechanics and Multidisciplinary Design.Berlin Heidelberg:Springer, 2007: 247-259. 

    [25]KRAJNOVIC S, BASARA B. LES of the flow around Ahmed body with active flow control [C]//Turbulence and Interactions. Berlin Heidelberg: Springer, 2010: 247-254.

    猜你喜歡
    尾流
    考慮風(fēng)切變影響的三維尾流模型風(fēng)場(chǎng)實(shí)驗(yàn)
    基于風(fēng)場(chǎng)實(shí)驗(yàn)的三維疊加尾流模型特性研究
    船舶尾流圖像的數(shù)字化處理和特征描述技術(shù)
    尾流自導(dǎo)魚(yú)雷經(jīng)典三波束彈道導(dǎo)引律設(shè)計(jì)優(yōu)化?
    航空器尾流重新分類(RECAT-CN)國(guó)內(nèi)運(yùn)行現(xiàn)狀分析
    基于單-雙高斯模型擬合法的測(cè)風(fēng)激光雷達(dá)海上風(fēng)電機(jī)組尾流特征分析
    風(fēng)力機(jī)尾流理論模型研究現(xiàn)狀綜述
    飛機(jī)尾流的散射特性與探測(cè)技術(shù)綜述
    錐形流量計(jì)尾流流場(chǎng)分析
    水面艦船風(fēng)尾流效應(yīng)減弱的模擬研究
    av欧美777| 女人十人毛片免费观看3o分钟| 精品无人区乱码1区二区| 国产单亲对白刺激| 18美女黄网站色大片免费观看| 床上黄色一级片| 国产精品爽爽va在线观看网站| 亚洲最大成人手机在线| 高清在线国产一区| 99国产精品一区二区三区| 精品人妻1区二区| 国产一区二区亚洲精品在线观看| 99在线视频只有这里精品首页| 国产av一区在线观看免费| 亚洲不卡免费看| 欧美绝顶高潮抽搐喷水| eeuss影院久久| 乱人视频在线观看| 色视频www国产| 天堂av国产一区二区熟女人妻| 日本黄色视频三级网站网址| 亚洲熟妇中文字幕五十中出| 真人一进一出gif抽搐免费| 欧洲精品卡2卡3卡4卡5卡区| 婷婷亚洲欧美| 欧美在线黄色| 久久久精品欧美日韩精品| 中亚洲国语对白在线视频| 特大巨黑吊av在线直播| 看十八女毛片水多多多| avwww免费| 狠狠狠狠99中文字幕| 首页视频小说图片口味搜索| 成人高潮视频无遮挡免费网站| 久久精品人妻少妇| 久久久精品欧美日韩精品| bbb黄色大片| 国产在线精品亚洲第一网站| 日韩中文字幕欧美一区二区| 亚洲av.av天堂| 一区二区三区免费毛片| av黄色大香蕉| 熟女电影av网| 国产国拍精品亚洲av在线观看| 国产亚洲av嫩草精品影院| av黄色大香蕉| 每晚都被弄得嗷嗷叫到高潮| 露出奶头的视频| 欧美日韩福利视频一区二区| 韩国av一区二区三区四区| 午夜福利在线在线| 日日干狠狠操夜夜爽| 免费看日本二区| 欧美激情国产日韩精品一区| 丁香欧美五月| 日韩av在线大香蕉| 午夜久久久久精精品| 熟女电影av网| 国产国拍精品亚洲av在线观看| 1024手机看黄色片| 国产精品人妻久久久久久| 成人国产一区最新在线观看| 久久久久久大精品| 日本黄色片子视频| 成年版毛片免费区| 午夜福利高清视频| 欧美日韩瑟瑟在线播放| 18美女黄网站色大片免费观看| 亚洲国产精品成人综合色| 18禁黄网站禁片免费观看直播| 婷婷精品国产亚洲av在线| 美女大奶头视频| 国内揄拍国产精品人妻在线| 国产高清视频在线观看网站| 在线免费观看不下载黄p国产 | 欧美极品一区二区三区四区| 熟女电影av网| 国内揄拍国产精品人妻在线| 免费电影在线观看免费观看| 欧美午夜高清在线| 天堂动漫精品| 亚洲av.av天堂| 亚洲,欧美精品.| 我的女老师完整版在线观看| 12—13女人毛片做爰片一| 国产伦一二天堂av在线观看| 国产在线男女| 可以在线观看毛片的网站| 精品人妻偷拍中文字幕| 国内毛片毛片毛片毛片毛片| av黄色大香蕉| 精品一区二区三区视频在线观看免费| 神马国产精品三级电影在线观看| 精品一区二区免费观看| 久久精品影院6| 精品免费久久久久久久清纯| 日本一二三区视频观看| 简卡轻食公司| 国产又黄又爽又无遮挡在线| 久久热精品热| 国产伦人伦偷精品视频| 欧美极品一区二区三区四区| 三级男女做爰猛烈吃奶摸视频| 久久婷婷人人爽人人干人人爱| 性色av乱码一区二区三区2| 床上黄色一级片| 日韩欧美三级三区| 国产伦精品一区二区三区视频9| 草草在线视频免费看| 人人妻,人人澡人人爽秒播| 日本五十路高清| 亚洲精品在线观看二区| 黄色日韩在线| 9191精品国产免费久久| 免费av观看视频| 乱人视频在线观看| 男女做爰动态图高潮gif福利片| 亚洲国产精品999在线| 婷婷色综合大香蕉| 久久久久久久久久黄片| 色综合站精品国产| 90打野战视频偷拍视频| 国产精品98久久久久久宅男小说| 真人一进一出gif抽搐免费| 别揉我奶头~嗯~啊~动态视频| 中出人妻视频一区二区| 人妻丰满熟妇av一区二区三区| 91狼人影院| 久久精品国产亚洲av天美| 日日摸夜夜添夜夜添av毛片 | 日本免费a在线| 一级毛片久久久久久久久女| 精品久久久久久久末码| 99热这里只有是精品在线观看 | 一进一出好大好爽视频| 制服丝袜大香蕉在线| 国产亚洲欧美98| 人人妻人人看人人澡| 亚洲午夜理论影院| 免费av不卡在线播放| 日韩中文字幕欧美一区二区| 国产精品亚洲av一区麻豆| 久久久久久久午夜电影| 亚洲国产欧洲综合997久久,| 成年女人永久免费观看视频| 久久久久九九精品影院| 色av中文字幕| 免费黄网站久久成人精品 | 人妻制服诱惑在线中文字幕| 日本成人三级电影网站| 在线观看av片永久免费下载| 久久人人爽人人爽人人片va | xxxwww97欧美| 亚洲欧美日韩高清专用| 女人被狂操c到高潮| 欧美国产日韩亚洲一区| 搡老熟女国产l中国老女人| 亚洲中文字幕日韩| 中文字幕av成人在线电影| 精品久久久久久,| 精品人妻1区二区| 精品人妻1区二区| 97人妻精品一区二区三区麻豆| 三级男女做爰猛烈吃奶摸视频| 亚洲天堂国产精品一区在线| 啦啦啦观看免费观看视频高清| 亚洲国产精品久久男人天堂| 夜夜爽天天搞| 亚洲欧美日韩高清专用| 午夜激情欧美在线| 午夜激情欧美在线| 黄色丝袜av网址大全| 中文字幕av在线有码专区| 99热这里只有精品一区| 一区福利在线观看| 国产主播在线观看一区二区| 丰满的人妻完整版| 极品教师在线免费播放| 欧美日韩亚洲国产一区二区在线观看| 日韩中文字幕欧美一区二区| 嫩草影院新地址| 国产aⅴ精品一区二区三区波| 亚洲精品亚洲一区二区| 中亚洲国语对白在线视频| 国产精品野战在线观看| 国产欧美日韩一区二区三| 99久久无色码亚洲精品果冻| 好男人在线观看高清免费视频| 婷婷六月久久综合丁香| 高清日韩中文字幕在线| 日日干狠狠操夜夜爽| 中亚洲国语对白在线视频| 一卡2卡三卡四卡精品乱码亚洲| 在线观看免费视频日本深夜| 国产精品精品国产色婷婷| 直男gayav资源| 好男人电影高清在线观看| 69av精品久久久久久| 国产成人a区在线观看| 欧美又色又爽又黄视频| 18美女黄网站色大片免费观看| 午夜视频国产福利| 99热只有精品国产| 一级av片app| 国内精品美女久久久久久| 嫩草影院新地址| 天堂网av新在线| 亚洲精品日韩av片在线观看| 亚洲国产色片| 免费av不卡在线播放| 一二三四社区在线视频社区8| 99热这里只有是精品在线观看 | 岛国在线免费视频观看| 精品久久久久久久人妻蜜臀av| 91av网一区二区| 51午夜福利影视在线观看| 亚洲欧美激情综合另类| 久久伊人香网站| 亚洲三级黄色毛片| 国产毛片a区久久久久| 嫩草影视91久久| 日本 欧美在线| 看黄色毛片网站| 亚洲综合色惰| 精品无人区乱码1区二区| 又紧又爽又黄一区二区| 两个人视频免费观看高清| 啦啦啦观看免费观看视频高清| 自拍偷自拍亚洲精品老妇| 欧美性感艳星| 成人欧美大片| 我的女老师完整版在线观看| 一级av片app| 午夜福利在线观看免费完整高清在 | 日本精品一区二区三区蜜桃| 色尼玛亚洲综合影院| 国产视频一区二区在线看| 欧美bdsm另类| 他把我摸到了高潮在线观看| 国产精品精品国产色婷婷| 免费在线观看成人毛片| 精品一区二区三区人妻视频| 精品免费久久久久久久清纯| 国产伦精品一区二区三区视频9| 久久精品夜夜夜夜夜久久蜜豆| 午夜两性在线视频| 欧美乱色亚洲激情| 色哟哟哟哟哟哟| 久久中文看片网| 亚洲狠狠婷婷综合久久图片| 高清日韩中文字幕在线| 最近中文字幕高清免费大全6 | 欧美+亚洲+日韩+国产| 欧美黄色片欧美黄色片| av视频在线观看入口| 亚洲不卡免费看| 免费看光身美女| 午夜福利在线观看免费完整高清在 | 成人三级黄色视频| 亚洲专区中文字幕在线| 午夜久久久久精精品| 国产熟女xx| 麻豆久久精品国产亚洲av| 国产高清视频在线播放一区| 亚洲三级黄色毛片| 网址你懂的国产日韩在线| 欧美一区二区精品小视频在线| 欧美一区二区精品小视频在线| 国产不卡一卡二| 午夜福利在线观看吧| 特大巨黑吊av在线直播| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 中出人妻视频一区二区| 能在线免费观看的黄片| 亚洲人成网站高清观看| av在线蜜桃| 亚洲,欧美精品.| 超碰av人人做人人爽久久| 久99久视频精品免费| 国产三级黄色录像| 亚洲国产精品成人综合色| 美女 人体艺术 gogo| 亚洲欧美日韩高清在线视频| 国产午夜精品论理片| 嫁个100分男人电影在线观看| 中文字幕精品亚洲无线码一区| 久久久色成人| 日韩欧美在线二视频| 色5月婷婷丁香| 人人妻人人澡欧美一区二区| 久久99热这里只有精品18| 欧美日韩黄片免| 色在线成人网| 99在线人妻在线中文字幕| 日韩 亚洲 欧美在线| 午夜免费男女啪啪视频观看 | 九色成人免费人妻av| 真人一进一出gif抽搐免费| 精品国内亚洲2022精品成人| 日韩精品青青久久久久久| 一个人免费在线观看的高清视频| 亚洲人成电影免费在线| 伦理电影大哥的女人| 舔av片在线| 99久久精品一区二区三区| 日韩欧美三级三区| 国产在线精品亚洲第一网站| 亚洲av免费在线观看| 丰满人妻熟妇乱又伦精品不卡| 嫩草影视91久久| 国产欧美日韩一区二区精品| 麻豆久久精品国产亚洲av| 在线观看午夜福利视频| 成人国产综合亚洲| 91久久精品国产一区二区成人| 亚洲国产日韩欧美精品在线观看| 露出奶头的视频| 窝窝影院91人妻| 国产精品一区二区三区四区久久| 99热这里只有是精品50| 久久精品国产清高在天天线| 蜜桃久久精品国产亚洲av| 18禁黄网站禁片免费观看直播| 亚洲成a人片在线一区二区| 日本免费一区二区三区高清不卡| 久9热在线精品视频| 亚洲成人免费电影在线观看| 免费人成视频x8x8入口观看| h日本视频在线播放| 亚洲七黄色美女视频| 人人妻人人看人人澡| 国产欧美日韩精品亚洲av| av在线老鸭窝| 亚洲,欧美精品.| av黄色大香蕉| 丰满人妻一区二区三区视频av| 亚洲美女搞黄在线观看 | 亚洲一区二区三区色噜噜| 桃红色精品国产亚洲av| 国产一级毛片七仙女欲春2| 性色avwww在线观看| 在线观看66精品国产| 免费人成视频x8x8入口观看| 国产精品1区2区在线观看.| 91在线精品国自产拍蜜月| 午夜老司机福利剧场| 99久久九九国产精品国产免费| 美女xxoo啪啪120秒动态图 | 夜夜爽天天搞| 757午夜福利合集在线观看| 一夜夜www| 国产色婷婷99| 深夜精品福利| www.熟女人妻精品国产| 在线免费观看的www视频| 国产69精品久久久久777片| 麻豆成人av在线观看| 免费高清视频大片| 国产欧美日韩一区二区精品| 18+在线观看网站| 91久久精品电影网| 少妇裸体淫交视频免费看高清| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品在线美女| 国产精品自产拍在线观看55亚洲| 免费看日本二区| 人妻夜夜爽99麻豆av| 丝袜美腿在线中文| av专区在线播放| 少妇丰满av| 久9热在线精品视频| 99久久九九国产精品国产免费| 久久久久久大精品| 国产aⅴ精品一区二区三区波| 欧美一区二区精品小视频在线| 在线观看美女被高潮喷水网站 | 精品不卡国产一区二区三区| 日韩高清综合在线| 欧美日韩黄片免| 亚洲欧美日韩高清专用| 91字幕亚洲| 18禁黄网站禁片午夜丰满| 中文字幕av成人在线电影| 欧美xxxx性猛交bbbb| 国产精品久久久久久亚洲av鲁大| 成年免费大片在线观看| 精品人妻视频免费看| 国产v大片淫在线免费观看| 日本黄大片高清| 色播亚洲综合网| 国产熟女xx| 88av欧美| 精品久久久久久久久久久久久| 波多野结衣高清无吗| 亚洲欧美日韩卡通动漫| 国产在视频线在精品| 久久精品人妻少妇| 免费观看人在逋| 国产久久久一区二区三区| 午夜福利高清视频| 最近中文字幕高清免费大全6 | 一区福利在线观看| 国产精品久久视频播放| 精品免费久久久久久久清纯| 亚洲五月天丁香| 久久精品国产亚洲av天美| 成人国产综合亚洲| 草草在线视频免费看| 精品一区二区三区视频在线观看免费| 97超视频在线观看视频| 久久6这里有精品| 国内精品久久久久精免费| 日本免费一区二区三区高清不卡| 亚洲无线在线观看| 精品久久久久久久久亚洲 | 12—13女人毛片做爰片一| 国产一区二区三区在线臀色熟女| 内地一区二区视频在线| 在线十欧美十亚洲十日本专区| 神马国产精品三级电影在线观看| 我要搜黄色片| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲精品久久久com| 麻豆一二三区av精品| 2021天堂中文幕一二区在线观| 有码 亚洲区| 中文字幕精品亚洲无线码一区| 99精品久久久久人妻精品| 99久久成人亚洲精品观看| 亚洲一区二区三区色噜噜| 观看美女的网站| eeuss影院久久| 国产私拍福利视频在线观看| 搡老妇女老女人老熟妇| 亚洲第一欧美日韩一区二区三区| 三级男女做爰猛烈吃奶摸视频| 欧美精品国产亚洲| 一a级毛片在线观看| 午夜福利视频1000在线观看| 久久久久精品国产欧美久久久| 亚洲国产高清在线一区二区三| 国产主播在线观看一区二区| 成人av在线播放网站| 看十八女毛片水多多多| 琪琪午夜伦伦电影理论片6080| 色av中文字幕| 小说图片视频综合网站| 有码 亚洲区| 真实男女啪啪啪动态图| 亚洲国产日韩欧美精品在线观看| 亚洲va日本ⅴa欧美va伊人久久| 天堂网av新在线| 国产野战对白在线观看| 亚洲欧美清纯卡通| 在线观看舔阴道视频| 亚洲综合色惰| 老司机午夜十八禁免费视频| 一进一出抽搐gif免费好疼| 久久久久免费精品人妻一区二区| 色噜噜av男人的天堂激情| 我要搜黄色片| 国产精品久久久久久久久免 | 亚洲欧美日韩东京热| 国产在线男女| 日韩人妻高清精品专区| 性色avwww在线观看| 天堂av国产一区二区熟女人妻| 色噜噜av男人的天堂激情| 日韩亚洲欧美综合| 国产老妇女一区| 欧美一区二区亚洲| 日日摸夜夜添夜夜添av毛片 | 午夜精品在线福利| 国产av在哪里看| 亚洲人成网站在线播| 亚洲av美国av| 在线观看舔阴道视频| 毛片一级片免费看久久久久 | 国产真实乱freesex| 欧美激情国产日韩精品一区| 宅男免费午夜| 精品无人区乱码1区二区| 网址你懂的国产日韩在线| 欧美成人性av电影在线观看| 亚洲人成网站在线播放欧美日韩| 中文字幕精品亚洲无线码一区| 日韩中字成人| 日韩大尺度精品在线看网址| 日本免费a在线| 在线观看免费视频日本深夜| 在线免费观看不下载黄p国产 | av视频在线观看入口| 一区二区三区激情视频| 国产视频内射| 男女下面进入的视频免费午夜| 精品一区二区三区av网在线观看| 亚洲性夜色夜夜综合| 99久久九九国产精品国产免费| 久久久色成人| 91av网一区二区| 国产精品一区二区免费欧美| 国内毛片毛片毛片毛片毛片| 中文亚洲av片在线观看爽| 精品无人区乱码1区二区| 久久亚洲精品不卡| a级一级毛片免费在线观看| av视频在线观看入口| 欧美日韩黄片免| 亚洲一区高清亚洲精品| 国产精品,欧美在线| 小蜜桃在线观看免费完整版高清| 啪啪无遮挡十八禁网站| 一区二区三区高清视频在线| 老司机福利观看| 一二三四社区在线视频社区8| 欧美激情国产日韩精品一区| 青草久久国产| 久久国产乱子免费精品| 天堂网av新在线| 国产黄色小视频在线观看| 99视频精品全部免费 在线| 很黄的视频免费| 国产伦精品一区二区三区视频9| 国产精品野战在线观看| 夜夜爽天天搞| 91午夜精品亚洲一区二区三区 | 网址你懂的国产日韩在线| 国产人妻一区二区三区在| 91av网一区二区| 国产色婷婷99| 一区二区三区四区激情视频 | 99热这里只有是精品在线观看 | 色综合婷婷激情| 中亚洲国语对白在线视频| 好男人电影高清在线观看| 国产久久久一区二区三区| 亚洲av日韩精品久久久久久密| h日本视频在线播放| 日韩 亚洲 欧美在线| 亚洲美女视频黄频| 午夜日韩欧美国产| 一区二区三区免费毛片| 亚洲精品一区av在线观看| 观看免费一级毛片| 亚洲av一区综合| 不卡一级毛片| 成人亚洲精品av一区二区| 国产av不卡久久| 白带黄色成豆腐渣| 51国产日韩欧美| 热99在线观看视频| 丁香六月欧美| 亚洲精品影视一区二区三区av| 久久久久久九九精品二区国产| 蜜桃久久精品国产亚洲av| 亚洲无线在线观看| 在线观看免费视频日本深夜| 欧美丝袜亚洲另类 | 久久久久国内视频| 自拍偷自拍亚洲精品老妇| 日本一本二区三区精品| 亚洲七黄色美女视频| ponron亚洲| 免费黄网站久久成人精品 | 老熟妇乱子伦视频在线观看| 国产美女午夜福利| 久久久成人免费电影| 91在线观看av| 国产亚洲av嫩草精品影院| 日本与韩国留学比较| 99热这里只有是精品50| 69人妻影院| 老熟妇仑乱视频hdxx| 免费黄网站久久成人精品 | 网址你懂的国产日韩在线| 国产人妻一区二区三区在| 久久这里只有精品中国| 亚洲成av人片在线播放无| 亚洲精品影视一区二区三区av| 久久久久亚洲av毛片大全| 老女人水多毛片| 性色av乱码一区二区三区2| 12—13女人毛片做爰片一| 亚洲美女搞黄在线观看 | 国产精品一区二区三区四区久久| 在线十欧美十亚洲十日本专区| 精品熟女少妇八av免费久了| 丰满的人妻完整版| 精品一区二区三区av网在线观看| 久久久久久久午夜电影| 国产亚洲精品综合一区在线观看| 午夜亚洲福利在线播放| 国产精品一区二区三区四区免费观看 | 国产麻豆成人av免费视频| 久久精品综合一区二区三区| 国产69精品久久久久777片| 黄色日韩在线| 男插女下体视频免费在线播放| 日本免费一区二区三区高清不卡| xxxwww97欧美| 日韩成人在线观看一区二区三区| 蜜桃久久精品国产亚洲av| 宅男免费午夜| 天堂av国产一区二区熟女人妻| 日本免费a在线| 一个人看的www免费观看视频| 99久久无色码亚洲精品果冻| 国产视频一区二区在线看| 桃红色精品国产亚洲av| 两个人的视频大全免费| 亚洲欧美激情综合另类| 亚洲男人的天堂狠狠| 午夜福利欧美成人| av黄色大香蕉| 国产av不卡久久|