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

通風(fēng)管道彎頭加裝導(dǎo)流片對(duì)流噪聲的影響

2024-12-28 00:00:00王景溫華兵陳聰鄭鏡華
機(jī)械制造與自動(dòng)化 2024年6期

摘要:針對(duì)空調(diào)通風(fēng)系統(tǒng)部件彎頭流噪聲問(wèn)題,以管道出口處的流噪聲為分析對(duì)象,基于有限元方法研究通風(fēng)管道彎頭內(nèi)布置導(dǎo)流片對(duì)管道內(nèi)部流場(chǎng)及出口流噪聲的影響,探明導(dǎo)流片形狀及數(shù)量參數(shù)的作用規(guī)律。模擬分析結(jié)果表明:管道出口處的流噪聲是由于管道中湍流產(chǎn)生的漩渦以及壓力脈動(dòng)所產(chǎn)生的,加裝直導(dǎo)流片反而會(huì)使出口流噪聲增大,而弧形導(dǎo)流片可以起到穩(wěn)流的作用,能減少壁面的漩渦結(jié)構(gòu)和大小,減小脈動(dòng)壓力強(qiáng)度,降低出口流噪聲的聲壓級(jí)水平;隨著弧形導(dǎo)流片數(shù)量增加,出口流噪聲的聲壓級(jí)下降明顯。

關(guān)鍵詞:空調(diào)器;導(dǎo)流片;彎頭;漩渦;湍流;流噪聲

中圖分類(lèi)號(hào):TB535文獻(xiàn)標(biāo)志碼:B文章編號(hào):1671-5276(2024)06-0088-05

Abstract:For the flow noise at the elbow of air conditioning and ventilation system components, with the flow noise at the outlet of the duct taken as the reserch object, the influence of the flow guide sheet arranged inside the elbow of the ventilation duct over the internal flow field and the outlet flow noise is studied based on the finite element method, and the action law of the shape and number parameters of the flow guide sheet is explored. The results show that the flow noise at the outlet of the duct is caused by the vortex and pressure pulsation generated by the turbulent flow in the duct, and the addition of straight deflectors will increase the outlet flow noise, while the curved deflectors can play a role in stabilizing the flow, which reduces the vortex structure and size of the wall, the pulsation pressure intensity, and the sound pressure level of the outlet flow noise. The more the number of curved deflectors increases, the more obviously the sound pressure level of flow noise at the outlet decreases.

Keywords:air conditioner; deflector sheet; elbow; vortex; turbulence; flow noise

0引言

空調(diào)通風(fēng)系統(tǒng)中管道的特點(diǎn)是尺寸大而且復(fù)雜,是由許多直管道和彎管所組成的。流體流過(guò)彎管是一種比較復(fù)雜的現(xiàn)象。在彎管區(qū)域,靠近彎管內(nèi)側(cè)的流體具有較高的速度,而靠近彎管外側(cè)的流體具有較低的速度,從而會(huì)產(chǎn)生比較大的壓力梯度,由于壓力梯度的產(chǎn)生,在流體中產(chǎn)生了一些不平衡的力。比較大的壓力梯度會(huì)在流動(dòng)中產(chǎn)生渦流,這就會(huì)增加湍流的水平,從而會(huì)導(dǎo)致通風(fēng)管道的結(jié)構(gòu)振動(dòng)和流體引起的流噪聲。

目前,國(guó)內(nèi)外一些研究人員對(duì)管道內(nèi)的流場(chǎng)情況進(jìn)行了大量的研究,并且對(duì)管道內(nèi)的流體特性進(jìn)行了數(shù)值模擬和實(shí)驗(yàn)。通過(guò)對(duì)管道流噪聲的實(shí)驗(yàn)研究可以得到其內(nèi)部的傳遞、衰減特性及影響因素的作用規(guī)律。考慮到現(xiàn)場(chǎng)采樣等實(shí)驗(yàn)方法受客觀條件限制,難以得到完整的物理信息,故很難從根源上分析流噪聲的機(jī)制。研究人員從流場(chǎng)的角度對(duì)管道流場(chǎng)進(jìn)行了研究分析,GAO等[1]研究了在空調(diào)通風(fēng)分流三通中安裝導(dǎo)流葉片的減阻方法,提出了安裝導(dǎo)流葉片的合理位置,優(yōu)化了導(dǎo)流葉片的形狀,分析了分流三通在不同流速和長(zhǎng)徑比下的阻力特性,通過(guò)實(shí)驗(yàn)驗(yàn)證了分流三通的優(yōu)化效果。JURAEVA等[2]利用數(shù)值模擬和實(shí)際應(yīng)用研究了導(dǎo)流板對(duì)地鐵隧道內(nèi)風(fēng)道均勻流動(dòng)的影響。結(jié)果表明:在地鐵隧道兩側(cè)安裝導(dǎo)流板可以獲得均勻的流動(dòng)。AI等[3]通過(guò)數(shù)值模擬和實(shí)驗(yàn)研究了5.75m長(zhǎng)直管道內(nèi)部的壓力損失,提出了一種預(yù)測(cè)壓力損失的數(shù)值模擬方法。SUN等[4]通過(guò)數(shù)值模擬研究了加裝和未加裝導(dǎo)流片的90°彎頭對(duì)圓形管道的影響,分析了流場(chǎng)特性,并提出導(dǎo)流葉片應(yīng)設(shè)置在彎管中的最佳位置。CRAWFORD等[5]對(duì)90°彎頭中湍流的壓降進(jìn)行了實(shí)驗(yàn)調(diào)查。所有的測(cè)量數(shù)據(jù)表明流動(dòng)顯示出軸向?qū)ΨQ(chēng)的特點(diǎn),與充分發(fā)展的管道流動(dòng)相比,湍流的法向應(yīng)力分布更加均勻,彎頭曲率的影響是加速管道流動(dòng)中的漩渦衰減。魏立明等[6]采用數(shù)值仿真模擬的方式,對(duì)于直板式節(jié)流閥閥芯進(jìn)行數(shù)值模擬,得到不同閥孔的壓降曲線、流場(chǎng)曲線、速度變化云圖等。經(jīng)過(guò)定量分析得到恒定流速下直板式節(jié)流閥的壓降規(guī)律。數(shù)值模擬結(jié)果表明,在相同的開(kāi)度下,隨著平板式節(jié)流閥閥芯長(zhǎng)軸長(zhǎng)度的增加,管道截面的最大流速有下降的趨勢(shì),隨著開(kāi)度從20%增大到40%,壓降的曲線趨于平緩。MODI等[7]在背景噪聲低于30dB(A)情況下,測(cè)量了不同流速下管口監(jiān)測(cè)點(diǎn)噪聲頻譜特性。在數(shù)值計(jì)算方面,MORI等[8]利用計(jì)算流體軟件和聲學(xué)計(jì)算軟件相結(jié)合的方法,對(duì)空調(diào)系統(tǒng)內(nèi)部流噪聲進(jìn)行了數(shù)值計(jì)算研究。WU等[9]通過(guò)計(jì)算流體動(dòng)力學(xué)(CFD)和計(jì)算流體動(dòng)力聲學(xué)(CAA)軟件來(lái)模擬分體式空調(diào)室內(nèi)機(jī)的空氣動(dòng)力學(xué)和聲學(xué)性能,提出了一種可以提高流速但不增加流噪聲的聲壓級(jí)優(yōu)化方法。

目前已經(jīng)有很多關(guān)于減少空調(diào)通風(fēng)管道局部組件的阻力和減少壓力損失的研究,但是還沒(méi)有明確的管道系統(tǒng)加裝導(dǎo)流片的聲學(xué)優(yōu)化設(shè)計(jì)方案。

本文以空調(diào)通風(fēng)管道彎頭作為研究對(duì)象,運(yùn)用Lighthill方程對(duì)彎頭出口由氣流所引發(fā)的流噪聲源進(jìn)行分析。采用計(jì)算流體力學(xué)和聲學(xué)有限元聯(lián)合仿真的方法,將彎管加裝和未加裝導(dǎo)流片非穩(wěn)態(tài)流場(chǎng)邊界脈動(dòng)壓力導(dǎo)入聲學(xué)有限元網(wǎng)格,轉(zhuǎn)換為在流場(chǎng)中的等效理論聲源,對(duì)彎管內(nèi)部空間是否加裝導(dǎo)流片的流體自輻射聲場(chǎng)進(jìn)行數(shù)值仿真,為空調(diào)通風(fēng)部件彎頭流噪聲的控制提供設(shè)計(jì)參考依據(jù)。

1氣動(dòng)聲學(xué)理論

1.1Navier-Stokes方程

本文主要通過(guò)分析彎管中的壓力分布、速度分布等流場(chǎng)信息研究加裝導(dǎo)流片下彎管的流動(dòng)行為。采用不可壓縮Navier-Stokes方程來(lái)描述流體區(qū)域的連續(xù)性、動(dòng)量和能量:

1.2大渦模擬(LES)湍流模型

大渦模擬湍流模型如下所示。

式中:ui和uj為過(guò)濾后速度分量;μ為湍流黏度系數(shù);τij為亞格子尺度應(yīng)力。

為使控制方程封閉,目前采用較多的亞格子模型是漩渦黏性模型。

式中:δij為克羅內(nèi)克函數(shù);uT為亞格子渦黏性系數(shù);sij為雷諾尺度應(yīng)力張量。

1.3FW-H方程

在Lighthill聲類(lèi)比理論基礎(chǔ)上,推導(dǎo)出FW-H方程:

式中:c0為遠(yuǎn)場(chǎng)聲速;p為觀測(cè)點(diǎn)聲壓;ui為xi方向的流體速度分量;un和vn分別為垂直于積分面的流體速度分量和積分面移動(dòng)速度分量;δ(f)為Dirac函數(shù);H(f)為Heaviside函數(shù);pij和Tij分別為應(yīng)力張量和Lighthill張量。

2彎管流場(chǎng)模擬結(jié)果分析

2.1驗(yàn)證方法的正確性

文獻(xiàn)[10]中對(duì)突擴(kuò)、突縮變截面管道的傳聲特性進(jìn)行了實(shí)驗(yàn)。實(shí)驗(yàn)工況中空氣流經(jīng)突變管道時(shí)由于湍流產(chǎn)生的流噪聲是唯一的噪聲源,可看作一段存在氣流通過(guò)的突擴(kuò)、突縮變截面管道,這是由于局部構(gòu)件引起的管道形式突變對(duì)流場(chǎng)所產(chǎn)生的影響。

圖1所示為突變管道的結(jié)構(gòu)示意圖。參數(shù)為:入口段L1=200mm;出口段L2=500mm;d1=d2=38mm;擴(kuò)張腔的長(zhǎng)度L=280mm;d=200mm。實(shí)驗(yàn)中管道入口流速為45m/s。

對(duì)圖1所示物理模型使用ANSYS ICEM進(jìn)行六面體結(jié)構(gòu)化網(wǎng)格劃分,并且對(duì)邊界層進(jìn)行加密。壁面法向距離的無(wú)量綱值y+約為1,網(wǎng)格總數(shù)約為70萬(wàn)。采用RNG k-ε湍流模型通過(guò)Fluent軟件進(jìn)行管道內(nèi)流場(chǎng)計(jì)算,邊界條件為入口速度、出口壓力,壁面采用絕熱無(wú)滑移邊界。

Virtual.Lab進(jìn)行聲學(xué)計(jì)算時(shí),聲學(xué)網(wǎng)格尺寸的劃分依據(jù)為聲學(xué)模型需要的最大頻率,計(jì)算結(jié)果的精度需要通過(guò)控制最大網(wǎng)格尺寸來(lái)實(shí)現(xiàn),可以參考式(7)。

式中:L為單元?jiǎng)澐殖叽纾琺;

c為當(dāng)?shù)芈曀伲琺/s;

fmax為需要求解的最大噪聲頻率,Hz。

采用本文的數(shù)值模擬方法進(jìn)行仿真復(fù)現(xiàn)。首先在Fluent中進(jìn)行流場(chǎng)計(jì)算,生成格式為CGNS的壓力脈動(dòng)時(shí)域計(jì)算結(jié)果,導(dǎo)入LMS Virtual.Lab中通過(guò)傅里葉變換轉(zhuǎn)換為頻率邊界條件進(jìn)行聲場(chǎng)計(jì)算。圖2為突變管道出口監(jiān)測(cè)點(diǎn)處的聲壓級(jí)頻譜圖和本文通過(guò)數(shù)值模擬方法所得結(jié)果,數(shù)值模擬結(jié)果與文獻(xiàn)結(jié)果在各個(gè)頻段上峰值分布幾乎一致,表明數(shù)值模擬方法計(jì)算管道氣動(dòng)噪聲的可行性。

2.2彎管及導(dǎo)流片模型

空調(diào)通風(fēng)管道系統(tǒng)路徑長(zhǎng),分布廣,形式復(fù)雜,且大多為直管、彎管、T型管和三通等管道形式。管道入口段長(zhǎng)度約為2m;出口段的長(zhǎng)度約為1.5m;管道的半徑約為200mm;彎頭半徑R=150mm。由于所選彎頭位于空調(diào)通風(fēng)管道主管道,流速較大,管道流速為14m/s,導(dǎo)流片的厚度為5mm,形狀分為弧形導(dǎo)流片和直導(dǎo)流片。模型參數(shù)及監(jiān)測(cè)點(diǎn)示意圖如圖3所示。

2.3通風(fēng)管道管內(nèi)流場(chǎng)分析

為了分析流體在彎頭流場(chǎng)中的動(dòng)態(tài)信息,分別計(jì)算無(wú)導(dǎo)流片、弧形導(dǎo)流片和直導(dǎo)流片彎頭的流場(chǎng)分布,如圖4所示。空氣在彎頭中流動(dòng)的時(shí)候,彎頭外側(cè)壓力增大,內(nèi)側(cè)壓力減小,導(dǎo)致邊界層與壁面分離形成回流區(qū)。軸向通道的截面由壓力分布呈現(xiàn)不同的著色,流線表示瞬時(shí)流場(chǎng)的信息。如圖4(a)所示,在無(wú)導(dǎo)流片的情況下,氣體流過(guò)彎頭,彎頭外壁壓力增大,內(nèi)側(cè)壓力減小,形成壓差,靠近彎頭內(nèi)壁的一側(cè)形成大量的漩渦和較強(qiáng)的脈動(dòng)壓力。圖4(b)為弧形導(dǎo)流片數(shù)量N=1彎管的流線圖。由圖可知在弧形導(dǎo)流片數(shù)量N=1的情況下,A截面處沒(méi)有弧形導(dǎo)流片的分流作用,徑向截面的漩渦結(jié)構(gòu)較大;B截面由于弧形導(dǎo)流片的作用,流線分布均勻,流場(chǎng)穩(wěn)定;C截面經(jīng)過(guò)弧形導(dǎo)流片的穩(wěn)流,減渦作用,所以C截面的漩渦結(jié)構(gòu)減小,漩渦數(shù)量降低。由于弧形導(dǎo)流片的作用,彎頭內(nèi)外壁的壓差減小,彎頭處流場(chǎng)更穩(wěn)定,流線更均勻,彎頭內(nèi)壁沒(méi)有過(guò)多的漩渦。由圖4(c)、圖4(d)可以看出,隨著弧形導(dǎo)流片數(shù)量的增多,漩渦的數(shù)量減少,所以弧形導(dǎo)流片的數(shù)量越多,流場(chǎng)就越穩(wěn)定。圖4(e)、圖4(f)和圖4(g)分別為直導(dǎo)流片數(shù)量為1、2和3的彎頭的流線圖。可以看出直導(dǎo)流片會(huì)導(dǎo)致彎管內(nèi)外壁的壓力差、靠近壁面的脈動(dòng)壓力和漩渦數(shù)量的增加。隨著直導(dǎo)流片數(shù)量的增加,漩渦尺度增加,漩渦數(shù)量增多,湍流強(qiáng)度和脈動(dòng)壓力增大。從3種不同數(shù)量直導(dǎo)流片彎頭處均可以看出:C截面處的漩渦結(jié)構(gòu)尺度比A截面大,漩渦數(shù)量比A截面多,說(shuō)明直導(dǎo)流片不能起到減渦和優(yōu)化流場(chǎng)的作用。

3彎管聲場(chǎng)模擬結(jié)果分析

建立兩種不同類(lèi)型的導(dǎo)流片,一種為直導(dǎo)流片,一種為弧形導(dǎo)流片,厚度均為5mm。直導(dǎo)流片的長(zhǎng)度為150mm,弧形導(dǎo)流片半徑為150mm。

為了分析和驗(yàn)證通風(fēng)管路內(nèi)部加裝導(dǎo)流葉片的降噪效果,建立7種不同類(lèi)型的彎頭模型。基于聲學(xué)有限元方法(FEM),利用Virtual.Lab軟件對(duì)彎頭聲場(chǎng)和出口相同位置處監(jiān)測(cè)點(diǎn)場(chǎng)點(diǎn)進(jìn)行聲學(xué)計(jì)算。

3.1導(dǎo)流片的形狀對(duì)流噪聲的影響

圖5所示為不同形狀導(dǎo)流片的出口聲壓級(jí)。3種不同類(lèi)型導(dǎo)流片的聲壓級(jí)頻率峰值出現(xiàn)在400Hz附近,直導(dǎo)流片彎頭出口處總聲壓級(jí)為63dB(A),主要原因?yàn)橹睂?dǎo)流片導(dǎo)致彎頭流場(chǎng)紊亂,湍流強(qiáng)度增大和漩渦增多,彎頭出口處的流噪聲值比無(wú)導(dǎo)流片高,弧形導(dǎo)流片由于增加了一個(gè)弧形結(jié)構(gòu),所以弧形導(dǎo)流片彎頭壁面脈動(dòng)壓力、湍流強(qiáng)度和漩渦大小降低,彎頭出口處的聲壓級(jí)水平降低。

3.2導(dǎo)流片的數(shù)量對(duì)流噪聲的影響

圖6、圖7分別為不同數(shù)量直導(dǎo)流片和弧形導(dǎo)流片彎頭出口聲壓級(jí)。經(jīng)過(guò)A聲級(jí)計(jì)權(quán)和1/3倍頻程處理后,峰值頻率在400Hz左右。無(wú)導(dǎo)流片情況下,彎頭出口總聲壓級(jí)為60dB(A)。如圖8所示,直導(dǎo)流片數(shù)量N=1,彎頭出口總聲壓級(jí)為63dB(A)。總聲壓級(jí)增加,說(shuō)明加裝直導(dǎo)流片會(huì)增加彎頭流場(chǎng)漩渦數(shù)量和湍流強(qiáng)度,出口流噪聲聲壓級(jí)增加;直導(dǎo)流片數(shù)量越增加,出口流噪聲聲壓級(jí)越大,說(shuō)明直導(dǎo)流片不能降低彎頭出口流噪聲水平。加裝一個(gè)弧形導(dǎo)流片后,彎頭出口總聲壓級(jí)為54dB(A),比無(wú)導(dǎo)流片彎頭降低6dB(A),說(shuō)明弧形導(dǎo)流片有降噪效果。由圖7可以看出不同數(shù)量弧形導(dǎo)流片彎頭的出口聲壓級(jí)頻譜曲線的峰值在400~500Hz,弧形導(dǎo)流片數(shù)量N=2時(shí),總聲壓級(jí)降低5dB(A);弧形導(dǎo)流片數(shù)量N=3時(shí),出口氣動(dòng)噪聲聲壓級(jí)比無(wú)導(dǎo)流片的情況時(shí),下降30dB(A)。由此可以得出如下結(jié)論:在管道空間允許的情況,加裝的弧形導(dǎo)流片越多,彎頭出口流噪聲的聲壓級(jí)越低。

4結(jié)語(yǔ)

本文采用有限元方法分析了不同類(lèi)型導(dǎo)流片對(duì)空調(diào)通風(fēng)系統(tǒng)彎頭對(duì)于出口流噪聲的影響,得出如下結(jié)論。

1)在無(wú)導(dǎo)流片的情況下,當(dāng)氣體流經(jīng)彎頭時(shí),彎頭外側(cè)壓力增大,內(nèi)側(cè)壓力減小,從而形成壓差,靠近彎頭內(nèi)壁的一側(cè)形成了較多的渦流,產(chǎn)生流噪聲。流噪聲聲壓級(jí)的頻率峰值與是否加裝導(dǎo)流片以及導(dǎo)流片的數(shù)量無(wú)關(guān)。

2)加裝直導(dǎo)流片時(shí),導(dǎo)致彎頭外壁的壓力差進(jìn)一步增大,形成了較多的渦流,出口流噪聲的聲壓級(jí)反而提高,且隨著直導(dǎo)流片的增加而提高。

3)相較于直導(dǎo)流片,加裝弧形導(dǎo)流片能夠有效地調(diào)節(jié)管道內(nèi)部的流場(chǎng),降低彎頭處的湍流強(qiáng)度和壓力脈動(dòng)水平,減小彎管肘部的漩渦數(shù)量和大小,降低管道流噪聲,且隨著弧形導(dǎo)流片增加,出口流噪聲的聲壓級(jí)降低。

參考文獻(xiàn):

[1] GAO R,LIU K K,LI A G,et al. Biomimetic duct tee for reducing the local resistance of a ventilation and air-conditioning system[J]. Building and Environment,2018,129:130-141.

[2] JURAEVA M,RYU K J,JEONG S H,et al. Influences of the train-wind and air-curtain to reduce the particle concentration inside a subway tunnel[J]. Tunnelling and Underground Space Technology,2016,52:23-29.

[3] AI Z T,MAK C M. Pressure losses across multiple fittings in ventilation ducts[J]. The Scientific World Journal,2013,2013:195763.

[4] SUN T,ZHANG Y,WANG Z Y. Research on flow in 90° curved duct with round section[J]. Physics Procedia,2012,24:692-699.

[5] CRAWFORD N M,CUNNINGHAM G,SPENCE S W T. An experimental investigation into the pressure drop for turbulent flow in 90° elbow bends[J]. Proceedings of the Institution of Mechanical Engineers,Part E:Journal of Process Mechanical Engineering,2007,221(2):77-88.

[6] 魏立明,李偉華,王興義,等. 基于CFD的直板式節(jié)流閥壓降特性研究[J]. 機(jī)械制造與自動(dòng)化,2022,51(4):85-88.

[7] MODI P P,JAYANTI S. Pressure losses and flow maldistribution in ducts with sharp bends[J]. Chemical Engineering Research and Design,2004,82(3):321-331.

[8] MORI M,MASUMOTO T,ISHIHARA K. Study on acoustic,vibration and flow induced noise characteristics of T-shaped pipe with a square cross-section[J]. Applied Acoustics,2017,120:137-147.

[9] WU C J,WANG X J,TANG H B. Transmission loss prediction on a single-inlet/double-outlet cylindrical expansion-chamber muffler by using the modal meshing approach[J]. Applied Acoustics,2008,69(2):173-178.

[10] 楊杰,覃國(guó)周,劉萬(wàn)里. 存在氣流時(shí)消聲器傳聲損失的數(shù)值計(jì)算[J]. 車(chē)輛與動(dòng)力技術(shù),2015(1):42-46.

收稿日期:20230328

第一作者簡(jiǎn)介:王景文(1996—),男,江蘇鎮(zhèn)江人,碩士研究生,研究方向?yàn)榇拜啓C(jī)設(shè)備及系統(tǒng)的振動(dòng)噪聲控制,zachary1631996@163.com。

DOI:10.19344/j.cnki.issn1671-5276.2024.06.017

主站蜘蛛池模板: 国产农村妇女精品一二区| 中国国产高清免费AV片| 999福利激情视频| 亚洲乱强伦| 欧美人在线一区二区三区| 高潮爽到爆的喷水女主播视频| 免费高清自慰一区二区三区| 久久夜色精品国产嚕嚕亚洲av| 色偷偷男人的天堂亚洲av| 国产欧美日韩综合一区在线播放| 国产手机在线小视频免费观看| 香蕉视频在线精品| 日本免费a视频| 91免费观看视频| 精品一区二区三区水蜜桃| 91精品福利自产拍在线观看| 国产午夜一级毛片| 欧美成人午夜在线全部免费| 久久黄色视频影| 国产成年女人特黄特色大片免费| 日本在线亚洲| 国产麻豆福利av在线播放| 国产区免费| 欧美国产日韩在线| 五月婷婷亚洲综合| 日韩欧美一区在线观看| 99成人在线观看| 国产精品三级专区| 久久精品免费国产大片| 日本精品αv中文字幕| 97在线国产视频| 国产精品福利社| 日韩av在线直播| 亚洲成人一区二区| 成人亚洲国产| 久久久国产精品无码专区| 国产91视频免费观看| 亚洲最大福利视频网| 韩日无码在线不卡| 五月天久久综合国产一区二区| 国产一区二区丝袜高跟鞋| 操操操综合网| 久久精品亚洲热综合一区二区| 色哟哟国产成人精品| 97久久人人超碰国产精品| 97se亚洲综合在线| 国产精品专区第1页| 女人18毛片水真多国产| 亚洲国产欧美国产综合久久 | 欧洲亚洲一区| 亚洲第一视频网| 国产成人禁片在线观看| 欧美日韩午夜视频在线观看 | 久久www视频| 日韩黄色大片免费看| 欧美黄网站免费观看| 999国产精品永久免费视频精品久久| 91精品国产91久久久久久三级| 美女国产在线| 亚洲九九视频| 91美女视频在线| 男女性色大片免费网站| 有专无码视频| 免费大黄网站在线观看| 麻豆精品久久久久久久99蜜桃| 欧美一区日韩一区中文字幕页| 中文精品久久久久国产网址 | 亚洲视频免| 国产一区二区免费播放| 国产第一福利影院| 亚洲人成人无码www| 久久精品人人做人人| www.亚洲一区二区三区| 91成人在线观看视频| 久久无码高潮喷水| 日韩无码一二三区| 亚洲69视频| 婷婷亚洲综合五月天在线| 日韩国产一区二区三区无码| 久久香蕉国产线看精品| 欧洲成人在线观看| 影音先锋丝袜制服|