999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

受電弓艙對(duì)受電弓區(qū)域氣動(dòng)噪聲的影響

2021-12-20 08:48:46史佳偉圣小珍
振動(dòng)與沖擊 2021年23期
關(guān)鍵詞:區(qū)域

史佳偉,葛 帥,圣小珍

(西南交通大學(xué) 牽引動(dòng)力國家重點(diǎn)實(shí)驗(yàn)室,成都 610031)

隨著列車運(yùn)行速度的不斷提高,列車噪聲對(duì)周邊環(huán)境的影響愈加顯著,降噪已成為高速列車設(shè)計(jì)中一個(gè)重要的考慮因素。高速列車的降噪設(shè)計(jì)最終可以歸結(jié)為對(duì)主要噪聲源的降噪。其中,受電弓就是列車重要的噪聲源之一,作為列車車頂?shù)耐怀霾考渲苯邮艿礁咚贇饬鳑_擊,噪聲強(qiáng)度很大,是高速列車局部最強(qiáng)聲源[1]。此外,對(duì)于兩側(cè)安裝了聲屏障的線路而言,受電弓的噪聲會(huì)更加突顯,因?yàn)榱熊囅虏康脑肼暱梢员宦暺琳险趽酰図數(shù)穆曉磩t是暴露的[2]。綜上所述,控制受電弓的氣動(dòng)噪聲對(duì)降低高速列車的噪聲水平至關(guān)重要。

控制受電弓氣動(dòng)噪聲的一類重要措施是引入導(dǎo)流罩、隔聲板等附屬設(shè)備。這些設(shè)備可以使受電弓避免直接暴露于高速氣流中或隔斷受電弓噪聲向周圍傳播的路徑。對(duì)于這類降噪方法,國內(nèi)外學(xué)者進(jìn)行了大量研究。Noger等[3]在低噪聲風(fēng)洞中測(cè)試了TGV列車受電弓系統(tǒng)的氣動(dòng)噪聲,發(fā)現(xiàn)受電弓落弓槽后方的豎直面是主要的聲源部位,受電弓的存在會(huì)增大寬頻噪聲的幅值;與一般的空腔流動(dòng)不同,受電弓艙不會(huì)出現(xiàn)空腔共振的現(xiàn)象。Kurita等[4]通過聲源識(shí)別,發(fā)現(xiàn) FASTECH360S 試驗(yàn)列車車頂?shù)腪形隔聲板的前緣和后緣是主要噪聲源;與傳統(tǒng)隔聲板相比,Z型隔聲板的隔聲效果更好,但其本身產(chǎn)生的氣動(dòng)噪聲更大,因而總體上降噪效果并不明顯。Yu等[5-6]通過數(shù)值模擬對(duì)比分析了4種受電弓導(dǎo)流罩的降噪效果,發(fā)現(xiàn)僅有單獨(dú)使用兩側(cè)隔聲板的方案可以降低受電弓氣動(dòng)噪聲,且4種方案都會(huì)使受電弓區(qū)域的氣動(dòng)阻力增大。Zhang等[7-8]從降噪的角度研究了受電弓開閉口方式、導(dǎo)流罩結(jié)構(gòu)及導(dǎo)流罩安裝位置的影響。Kim等[9]將受電弓簡(jiǎn)化為單個(gè)圓柱,研究了受電弓艙附近的空氣流動(dòng)特性和噪聲特性。Kim等[10]還比較了3種不同的受電弓附屬設(shè)備對(duì)受電弓區(qū)域流場(chǎng)和噪聲的影響,認(rèn)為受電弓隔聲板的方案降噪效果最好。陳羽等[11]通過風(fēng)洞試驗(yàn)研究了受電弓導(dǎo)流罩對(duì)受電弓氣動(dòng)噪聲的影響,發(fā)現(xiàn)降弓狀態(tài)下導(dǎo)流罩有比較明顯的降噪效果,而升弓狀態(tài)下導(dǎo)流罩的降噪效果有限。

在中國標(biāo)準(zhǔn)動(dòng)車組上采用了受電弓艙的方案[12],即受電弓位置的車頂下沉一定高度形成一個(gè)凹腔,受電弓安裝于凹腔內(nèi),這也被稱為受電弓下沉式安裝。孔學(xué)舟等[13]和林鵬等[14]分別研究了這種安裝方式對(duì)受電弓區(qū)域氣動(dòng)阻力的影響,結(jié)果表明:受電弓艙的方案可以使受電弓區(qū)域?qū)崿F(xiàn)近20%的減阻,但目前關(guān)于受電弓艙對(duì)受電弓區(qū)域氣動(dòng)噪聲影響的研究還比較少。為此,本文基于大渦模擬 (large eddy simulation,LES)方法和FW-H (Ffowcs Williams-Hawkings)方程,通過數(shù)值模擬研究了受電弓艙對(duì)受電弓區(qū)域氣流流動(dòng)特性、表面氣動(dòng)聲源以及遠(yuǎn)場(chǎng)氣動(dòng)噪聲的影響,相關(guān)結(jié)果可為高速列車受電弓氣動(dòng)噪聲控制提供參考。

1 數(shù)值計(jì)算模型

1.1 幾何模型

本文采用三車模型,受電弓安裝在中間車二位端處,分別建立受電弓外置和下沉式(受電弓艙)安裝的列車模型。由于研究重點(diǎn)為受電弓區(qū)域的氣動(dòng)噪聲,故所用模型忽略了轉(zhuǎn)向架等對(duì)車頂流場(chǎng)影響較小的部件并對(duì)列車底部做了平整化處理,將車端風(fēng)擋簡(jiǎn)化為光滑曲面,以減少數(shù)值計(jì)算所需的網(wǎng)格數(shù)量。幾何模型,如圖1所示。

(a)整車模型

1.2 計(jì)算域及邊界條件

本文模擬列車以350 km/h的速度在明線上運(yùn)行的工況,相應(yīng)的馬赫數(shù)小于0.3,故將空氣視作不可壓縮流體[15]。圖2給出了列車外流場(chǎng)計(jì)算域和邊界條件的示意圖。其中,列車車頭鼻尖距離計(jì)算域入口100 m,車尾鼻尖距離計(jì)算域出口200 m,計(jì)算域?qū)挾群透叨确謩e為60 m和30 m,列車底部距地面0.376 m。計(jì)算域入口設(shè)定為速度入口,輸入列車運(yùn)行速度。出口設(shè)定為壓力出口,靜壓為0。計(jì)算域兩側(cè)及頂部邊界設(shè)定為對(duì)稱邊界。計(jì)算域底部設(shè)置為滑移壁面以消除地面附面層的影響,滑移速度為來流速度。

圖2 計(jì)算域及邊界條件 (m)Fig.2 Computational domain and boundary conditions (m)

1.3 網(wǎng)格劃分

劃分切割體網(wǎng)格離散計(jì)算域,為避免網(wǎng)格差異對(duì)計(jì)算結(jié)果的影響,所有計(jì)算模型采用相同的網(wǎng)格尺寸進(jìn)行網(wǎng)格劃分。其中,受電弓表面網(wǎng)格尺寸控制在5~10 mm,車身部位網(wǎng)格的最大尺寸為60 mm,流線型車頭和車尾最大網(wǎng)格尺寸為30 mm,外流場(chǎng)區(qū)域最大網(wǎng)格尺寸為2 000 mm。對(duì)列車周圍和受電弓區(qū)域進(jìn)行了局部加密,列車周圍加密區(qū)域加密尺寸為60 mm,受電弓周圍兩個(gè)加密區(qū)的網(wǎng)格尺寸分別為15 mm和30 mm。在列車表面拉伸出初始高度0.1 mm,增長(zhǎng)率為1.2的8層棱柱層網(wǎng)格。各模型劃分后的網(wǎng)格數(shù)量均在2 400萬左右。圖3給出了受電弓下沉式安裝時(shí)受電弓區(qū)域和整個(gè)計(jì)算域的網(wǎng)格分布示意圖。

(a)受電弓區(qū)域面網(wǎng)格

1.4 求解過程

為了加快收斂速度,首先進(jìn)行穩(wěn)態(tài)模擬,之后以穩(wěn)態(tài)模擬結(jié)果作為初值開始瞬態(tài)模擬。穩(wěn)態(tài)模擬采用SST-kω模型,共迭代4 000步以確保計(jì)算收斂。瞬態(tài)模擬采用基于WALE亞格子應(yīng)力的LES模型,時(shí)間離散格式為二階隱式。瞬態(tài)計(jì)算時(shí)間步長(zhǎng)取為0.000 1 s,相應(yīng)的最大分析頻率為5 000 Hz,每個(gè)時(shí)間步內(nèi)迭代步數(shù)設(shè)為10步。瞬態(tài)計(jì)算總時(shí)長(zhǎng)設(shè)為0.4 s。對(duì)遠(yuǎn)場(chǎng)聲學(xué)信號(hào)的采集在瞬態(tài)模擬進(jìn)行到0.15 s時(shí)開始。

本文采用FW-H方程計(jì)算受電弓區(qū)域的遠(yuǎn)場(chǎng)噪聲。已有研究表明,受電弓的遠(yuǎn)場(chǎng)噪聲主要由受電弓表面的偶極子聲源貢獻(xiàn)[16-17],故本文將積分面取在固體表面,忽略了四極子源,同時(shí)由于列車靜止,故單極子源項(xiàng)為0,即相當(dāng)于僅考慮了偶極子源的影響。需要說明的是,使用不可壓縮流體模型結(jié)合FW-H方程計(jì)算受電弓區(qū)域的氣動(dòng)噪聲相當(dāng)于僅考慮了受電弓艙通過改變受電弓區(qū)域流場(chǎng)對(duì)受電弓區(qū)域氣動(dòng)噪聲產(chǎn)生的影響,而沒有考慮艙壁的反射作用。

2 流場(chǎng)分析

2.1 受電弓區(qū)域的流動(dòng)特征

圖4為兩種不同安裝方式下受電弓區(qū)域縱對(duì)稱面時(shí)均流場(chǎng)的速度矢量圖。由圖4可見,當(dāng)受電弓外置安裝在車頂時(shí),來流直接沖擊受電弓的所有部件并在滑板、平衡臂、底座支架及絕緣子等處發(fā)生分離,在這些部件的背風(fēng)側(cè)形成明顯的分離渦。當(dāng)受電弓下沉式安裝時(shí),流入受電弓艙的氣流在受電弓艙前部形成向下翻卷的漩渦,受電弓艙內(nèi)氣流流速明顯降低,受電弓艙內(nèi)部件周圍的氣流流動(dòng)狀態(tài)與外置安裝時(shí)有很大不同。受漩渦影響,位于受電弓艙上方距離受電弓艙較近的氣流會(huì)比較明顯的向下偏轉(zhuǎn),而距離受電弓艙較遠(yuǎn)的氣流則基本不受影響。

(a)外置-升弓

2.2 受電弓區(qū)域的尾渦結(jié)構(gòu)

圖5給出了當(dāng)使用Q準(zhǔn)則等值面所顯示的t=0.4 s時(shí)受電弓區(qū)域的尾渦結(jié)構(gòu),采用渦量值渲染。Q準(zhǔn)則將渦定義為渦量量級(jí)大于剪切應(yīng)變率量級(jí)的區(qū)域,Q值的計(jì)算式為[18]

(1)

式中:xi(i=1,2,3)為直角坐標(biāo)分量;ui為流體速度在xi方向的分量。

由圖5可見,Q值越大,渦區(qū)范圍越小。當(dāng)受電弓外置安裝于車頂時(shí),流經(jīng)受電弓時(shí)的氣流首先在受電弓表面分裂成泡沫狀的小渦,而后逐漸發(fā)展成蠕蟲狀或發(fā)夾狀的大渦。受電弓表面剛剛形成的渦具有較大的渦量,隨著渦向后運(yùn)動(dòng)遠(yuǎn)離受電弓,其渦量值逐漸減小。當(dāng)存在受電弓艙時(shí),氣流在受電弓艙前緣即發(fā)生分離。形成的分離渦向下游運(yùn)動(dòng),作用于受電弓表面。進(jìn)入受電弓艙的氣流流經(jīng)受電弓時(shí),會(huì)因受電弓的干擾再次發(fā)生分離,但分離渦的形態(tài)與外置安裝時(shí)有所不同。在升弓狀態(tài)下,可以觀察到受電弓底架和下臂桿附近渦的數(shù)量及渦量大小有一定程度的減小。在降弓狀態(tài)下,可以觀察到滑板 和平衡臂附近渦的數(shù)量及渦量大小有一定程度的減小。在受電弓的后方,流出受電弓艙的氣流與艙外的氣流發(fā)生匯聚,形成更大尺度的漩渦,同時(shí),還可以觀察到受電弓艙后緣拐角處流動(dòng)分離形成的尾渦,這些渦在受電弓后方很大的范圍內(nèi)存在,但渦量很小。

(a)外置-升弓,Q=500

3 受電弓表面氣動(dòng)聲源分析

圖6給出了兩種不同安裝方式下受電弓表面脈動(dòng)壓力標(biāo)準(zhǔn)差sp的大小和分布,該值反映了受電弓表面一點(diǎn)處壓力波動(dòng)的劇烈程度,可以用來表征受電弓表面一點(diǎn)處聲源的強(qiáng)度。sp的定義為

(2)

(a)外置-升弓

由圖6可見,受電弓艙的存在改變了受電弓表面氣動(dòng)噪聲源的強(qiáng)度和分布。受電弓外置安裝時(shí),平衡臂的上表面、迎風(fēng)側(cè)和背風(fēng)側(cè),后方滑板的迎風(fēng)側(cè)和背風(fēng)側(cè)、底座支架和絕緣子的背風(fēng)側(cè),底座支架前方橫梁的上下表面以及后方絕緣子和底座支架與下臂桿相連位置的迎風(fēng)側(cè)壓力波動(dòng)劇烈,是受電弓表面的主要?dú)鈩?dòng)聲源位置。當(dāng)受電弓下沉式安裝時(shí),底座支架和絕緣子的背風(fēng)側(cè),底座支架前方橫梁的上下表面,后方絕緣子以及底座支架與下臂桿相連位置的迎風(fēng)側(cè)壓力波動(dòng)明顯減弱,而底架前方橫梁和前排絕緣子的迎風(fēng)側(cè)、以及降弓狀態(tài)下前方滑板和平衡臂的上表面壓力波動(dòng)加強(qiáng)。此外,還可以發(fā)現(xiàn),在降弓狀態(tài)下,平衡臂表面的壓力波動(dòng)有所加強(qiáng),而后方滑板表面的壓力波動(dòng)有所減弱。

為了解釋不同安裝方式下受電弓表面聲源強(qiáng)度和分布的差異,圖7給出了不同安裝方式下當(dāng)t=0.4 s時(shí)受電弓區(qū)域縱對(duì)稱面的瞬時(shí)渦量圖,這些云圖更加清晰的顯示了受電弓的尾渦形態(tài)及其與受電弓各部件之間的相互作用。由圖7可見,當(dāng)存在受電弓艙時(shí),氣流在受電弓艙前緣即發(fā)生分離,分離渦沖擊底架、絕緣子和前方滑板(降弓狀態(tài))的迎風(fēng)側(cè),使這些位置的壓力脈動(dòng)加劇。但受電弓前緣形成的渦向下翻卷,艙內(nèi)氣流速度降低,這減弱了底架和絕緣子本身引起的流動(dòng)分離,從而使底架和絕緣子背風(fēng)側(cè)分離區(qū)的壓力脈動(dòng)較外置安裝時(shí)顯著降低。此外,由于底架和絕緣子處分離減弱,分離渦的尺度變小,渦量降低,其與下游部件的相互作用也較外置時(shí)變?nèi)酰@使得受電弓底架和絕緣子尾渦影響的區(qū)域(如底架與下臂桿相連位置的迎風(fēng)側(cè)和后方絕緣子迎風(fēng)側(cè))壓力脈動(dòng)減弱。此外,在降弓狀態(tài)下,可以觀察到滑板和平衡臂處的流動(dòng)分離也有一定程度的減弱。值得注意的是,降弓狀態(tài)下受電弓艙的存在使得前方滑板的尾渦向下偏轉(zhuǎn)撞擊平衡臂上表面,這可能是降弓狀態(tài)下平衡臂表面壓力波動(dòng)加劇的原因。

(a)外置-升弓

4 遠(yuǎn)場(chǎng)氣動(dòng)噪聲分析

第2章和第3章分析了受電弓艙對(duì)受電弓區(qū)域流場(chǎng)及受電弓表面氣動(dòng)聲源的影響,本節(jié)進(jìn)一步研究受電弓艙對(duì)受電弓區(qū)域遠(yuǎn)場(chǎng)氣動(dòng)噪聲的影響。用于評(píng)價(jià)受電弓艙對(duì)遠(yuǎn)場(chǎng)噪聲影響的遠(yuǎn)場(chǎng)聲學(xué)測(cè)點(diǎn),布置如圖8所示。其中,縱向測(cè)點(diǎn)沿z軸負(fù)向依次命名為q1~q5,測(cè)點(diǎn)間距均為3 m,q3測(cè)點(diǎn)與受電弓底部中心縱向坐標(biāo)一致。垂向測(cè)點(diǎn)沿y軸正向(y軸正向垂直紙面向外)命名為r1~r8,垂向間距均為1 m,其中r1測(cè)點(diǎn)距離地面高度為0.5 m,r4測(cè)點(diǎn)與q3測(cè)點(diǎn)重合,垂向測(cè)點(diǎn)距列車縱對(duì)稱面橫向距離均為7.5 m。橫向測(cè)點(diǎn)沿x軸正向依次命名為p1~p5,距計(jì)算域縱對(duì)稱面距離依次為7.5 m、15.0 m、25.0 m、30.0 m和50.0 m,p1測(cè)點(diǎn)與q3測(cè)點(diǎn)重合。

圖8 遠(yuǎn)場(chǎng)噪聲測(cè)點(diǎn)布置示意圖(不按比例)Fig.8 Layout of far field noise measurement points (not to scale)

圖9給出了分別以受電弓和受電弓艙(或外置狀態(tài)下的受電弓底板)為聲源,各測(cè)點(diǎn)聲壓級(jí)的計(jì)算結(jié)果。由圖9可見,各工況下受電弓輻射噪聲的聲壓級(jí)均較受電弓艙(或受電弓底板)高6~20 dB(A)。這說明無論受電弓外置安裝還是下沉式安裝,無論升弓還是降弓,受電弓區(qū)域向軌道兩側(cè)輻射的氣動(dòng)噪聲均由受電弓本身所主導(dǎo)。此外,對(duì)于本文所建立的受電弓模型,受電弓外置安裝且為降弓狀態(tài)時(shí)遠(yuǎn)場(chǎng)測(cè)點(diǎn)的聲壓級(jí)最大,受電弓下沉式安裝且為降弓狀態(tài)時(shí)遠(yuǎn)場(chǎng)測(cè)點(diǎn)的聲壓級(jí)最小。

統(tǒng)計(jì)兩種不同安裝方式下受電弓遠(yuǎn)場(chǎng)噪聲聲壓級(jí)在3個(gè)方向的平均變化量(下沉安裝的結(jié)果減去外置安裝時(shí)的結(jié)果),如圖10所示。

分析圖10中的數(shù)據(jù)可得:升弓狀態(tài)下受電弓下沉式的安裝方案在縱向和橫向輻射的噪聲較外置安裝時(shí)分別平均低約0.2 dB(A)和0.48 dB(A),垂向測(cè)點(diǎn)還略有增大。而降弓狀態(tài)下,受電弓下沉式的安裝與受電弓外置安裝時(shí)相比,3個(gè)方向測(cè)點(diǎn)的聲壓級(jí)都平均降低了約5 dB(A),降噪效果明顯。圖11給出了兩種安裝方式下測(cè)點(diǎn)q3的A計(jì)權(quán)聲壓級(jí)頻譜圖。由圖11可見,升弓狀態(tài)下,受電弓艙的存在主要降低了受電弓315 Hz以下頻段的噪聲。而在降弓狀態(tài)下,受電弓艙的存在幾乎使得整個(gè)分析頻段的噪聲都有所降低。由于A計(jì)權(quán)弱化了低頻的影響,故受電弓艙對(duì)升弓狀態(tài)下受電弓的降噪效果并不明顯,而對(duì)降弓狀態(tài)下受電弓的氣動(dòng)噪聲則有明顯的降噪效果。

(a)縱向測(cè)點(diǎn)

圖10 降噪量統(tǒng)計(jì)Fig.10 Noise reduction statistics

(a)升弓狀態(tài)

5 結(jié) 論

本文結(jié)合LES方法與FW-H方程通過數(shù)值模擬研究了受電弓艙對(duì)受電弓區(qū)域流場(chǎng)和氣動(dòng)噪聲的影響。分析計(jì)算結(jié)果,得到如下結(jié)論:

(1)受電弓艙會(huì)顯著影響受電弓區(qū)域的氣流流動(dòng)狀態(tài)進(jìn)而改變受電弓表面的氣動(dòng)聲源強(qiáng)度和分布。當(dāng)存在受電弓艙時(shí),來流在受電弓艙前緣發(fā)生分離,分離渦與受電弓相互作用,加劇受電弓表面的壓力波動(dòng)。但進(jìn)入受電弓艙的氣流向下偏轉(zhuǎn),流速降低,這使得受電弓本身所引起的流動(dòng)分離減弱,受電弓表面分離區(qū)及受到上游部件尾渦干擾的區(qū)域壓力波動(dòng)明顯減弱。

(2)無論受電弓是外置式安裝還是下沉式安裝,無論升弓還是降弓,受電弓區(qū)域向軌道兩側(cè)輻射的氣動(dòng)噪聲均由受電弓本身所主導(dǎo)。

(3)對(duì)于本文所建立的模型,在升弓狀態(tài)下,受電弓艙的存在主要降低了受電弓315 Hz以下頻段的噪聲。而在降弓狀態(tài)下,受電弓艙的存在幾乎使整個(gè)分析頻段的噪聲都有所降低。總體而言,受電弓艙對(duì)降弓狀態(tài)下受電弓的氣動(dòng)噪聲的降噪效果更好。

猜你喜歡
區(qū)域
分割區(qū)域
探尋區(qū)域創(chuàng)新的密碼
科學(xué)(2020年5期)2020-11-26 08:19:22
基于BM3D的復(fù)雜紋理區(qū)域圖像去噪
軟件(2020年3期)2020-04-20 01:45:18
小區(qū)域、大發(fā)展
商周刊(2018年15期)2018-07-27 01:41:20
論“戎”的活動(dòng)區(qū)域
區(qū)域發(fā)展篇
區(qū)域經(jīng)濟(jì)
關(guān)于四色猜想
分區(qū)域
公司治理與技術(shù)創(chuàng)新:分區(qū)域比較
主站蜘蛛池模板: 国产精品久久国产精麻豆99网站| 中国国产A一级毛片| 91免费观看视频| 亚洲资源在线视频| 色婷婷狠狠干| 精品無碼一區在線觀看 | 欧美a在线视频| 国产真实乱人视频| 久久国语对白| 国产精品亚洲欧美日韩久久| 久久国产精品电影| 九九这里只有精品视频| 亚洲一区二区三区国产精品| 亚洲精品国产成人7777| 亚洲一区二区成人| 亚洲国产精品久久久久秋霞影院| 免费国产高清精品一区在线| 免费啪啪网址| 欧美日韩午夜| 久久青草热| 中美日韩在线网免费毛片视频 | 国产极品美女在线观看| 久久一色本道亚洲| 三上悠亚一区二区| 激情综合图区| 免费精品一区二区h| 国产成人91精品| 国产超薄肉色丝袜网站| 98超碰在线观看| 波多野结衣一区二区三区AV| 18禁黄无遮挡网站| 精品伊人久久久香线蕉| 国产三级毛片| 女同国产精品一区二区| 午夜精品一区二区蜜桃| 丁香婷婷综合激情| 东京热一区二区三区无码视频| 美女一级毛片无遮挡内谢| 亚洲色偷偷偷鲁综合| 91久久夜色精品国产网站| 国产视频只有无码精品| 精品一区二区三区视频免费观看| 99热这里只有精品免费国产| 国产xx在线观看| 日本a级免费| 亚洲无码一区在线观看| 美女被操黄色视频网站| 国产成人综合久久精品下载| 久久黄色免费电影| www.99精品视频在线播放| 91精品久久久无码中文字幕vr| 男女男免费视频网站国产| 午夜福利亚洲精品| 欧美日韩北条麻妃一区二区| 色综合激情网| 强奷白丝美女在线观看| 国产色伊人| 亚洲欧美自拍一区| 欧美一级高清免费a| 一本大道香蕉久中文在线播放 | 亚洲无码A视频在线| 日韩成人在线视频| 98超碰在线观看| 国产成年女人特黄特色毛片免| 日韩精品免费在线视频| 欧美成人国产| 亚洲欧洲日韩久久狠狠爱| 91精品视频在线播放| 日韩区欧美区| 91年精品国产福利线观看久久| 中文字幕在线视频免费| 久久精品国产免费观看频道| 亚洲天堂首页| 亚洲第一区精品日韩在线播放| 青青热久麻豆精品视频在线观看| 成人午夜免费观看| 特级精品毛片免费观看| 波多野结衣爽到高潮漏水大喷| 精品视频91| 国产精品吹潮在线观看中文| 欧美色视频日本| 国产成人a在线观看视频|