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

    非均勻流條件下穿孔管消聲器傳遞損失計算方法與特性分析

    2022-12-05 07:49:28陳志響季振林馬慶鎮(zhèn)陳國強孟祥旗
    聲學技術(shù) 2022年5期
    關(guān)鍵詞:橫流邊界條件聲學

    陳志響,季振林,馬慶鎮(zhèn),陳國強,孟祥旗,張 超

    (1.濰柴動力股份有限公司內(nèi)燃機可靠性國家重點實驗室,山東 濰坊 261061;2.哈爾濱工程大學動力與能源工程學院,黑龍江 哈爾濱 150001)

    0 引言

    消聲器被廣泛應(yīng)用于內(nèi)燃機排氣噪聲控制,聲學性能是其主要的評價指標。頻域方法是消聲器聲學性能計算和分析中普遍使用的方法,主要包括有限元法和邊界元法[1]。這些方法均在頻率域求解線性波動方程,進而計算消聲器的聲學性能。Young等[2]首次使用有限元法計算了無流條件下消聲器的傳遞損失。徐貝貝等[3]利用三維有限元法分析了消聲器的聲學性能。由于當時計算資源的限制,有限元法只能在二維模型下進行計算,但實際的消聲器通常具有三維的幾何特征,因此當時的有限元法預測精度較低。隨著計算機性能的提高和存儲容量的增大,有限元法逐步應(yīng)用于計算和分析三維聲學問題。Mehdizadeh等[4]利用三維聲學有限元法計算了穿孔管消聲器的聲學性能,傳遞損失預測結(jié)果與實驗測量結(jié)果吻合良好。Chaitanya等[5]利用三維有限元法研究進出口管道的厚度對末端修正系數(shù)的影響,并將預測結(jié)果與實驗結(jié)果進行了對比。但是上面的計算過程中并沒有考慮流速對消聲器聲學性能的影響,消聲器作為排氣管路的組件,其消聲性能將不可避免地受到流體流動的影響。1988年,Chang等[6]基于有限差分格式離散了線性化的一維波動方程,在時間域內(nèi)求解了聲學時域信號,利用時頻信號分析和聲學計算得到了消聲器的傳遞損失。徐航手等[7]利用三維時域CFD方法計算了無流條件下單、雙級膨脹腔、直通穿孔管消聲器的傳遞損失。劉晨等[8]詳細研究了三維時域CFD方法用于計算有流狀態(tài)下穿孔管消聲器傳遞損失的改進策略,其計算結(jié)果與實驗測量值吻合較好。雖然時域方法考慮的因素較多,計算精度較高,但是其計算效率較低,并且對網(wǎng)格大小和時間步長的要求較為嚴格,目前該方法還很難用于實際的工程計算與分析[9]。

    在消聲器的管路中,流體的流動會加速聲波向下游的傳播,同時減慢聲波向上游的傳播。一般來說,消聲器不同的腔室中也會存在不同的流動特征。由于流動造成的聲速傳播的差異會影響消聲器的消聲性能。而且消聲器中經(jīng)常使用穿孔結(jié)構(gòu),利用小孔的聲阻抗效應(yīng)來達到消聲的目的,小孔與腔室的組合又形成了共振結(jié)構(gòu),流體的流動會影響共振結(jié)構(gòu)的共振頻率。

    因此,為了考慮流體的運流效應(yīng)對消聲器消聲性能的影響,本文利用計算流體力學(Computational Fluid Dynamics,CFD)數(shù)值方法計算消聲器內(nèi)部流場中的速度分布[10],再將速度分布通過網(wǎng)格映射轉(zhuǎn)移到聲學有限元模型中[11];而流體流動對穿孔聲阻抗的影響通過有流聲阻抗經(jīng)驗公式直接定義,然后進行聲學響應(yīng)的計算。本文使用CFD軟件FLUENT計算消聲器中不同節(jié)點處速度的分布;使用兩種聲阻抗經(jīng)驗公式直接定義小孔的聲阻抗,進而計算有流條件下穿孔管消聲器的傳遞損失。

    1 傳遞損失

    消聲器的傳遞損失只與本體結(jié)構(gòu)有關(guān),不受聲源阻抗與尾管輻射特性的影響,是消聲器聲學性能中最常使用的評價指標。當消聲器的進出口滿足平面波條件時,其傳遞損失可以表示為[12]

    其中:pin為消聲器進口處的入射聲壓,pout為消聲器出口在無反射邊界條件下的透射聲壓,Ain、Aout為消聲器進、出口處的面積。

    2 計算方法

    2.1 消聲器模型

    本文中所計算的消聲器模型為文獻[13]中的橫流式穿孔管消聲器和文獻[14]中的雙穿孔管消聲器。橫流式穿孔管消聲器和雙穿孔管消聲器的幾何模型分別如圖1、2所示,具體尺寸如表1所示。橫流式消聲器中空氣溫度T=347 K,該消聲器為軸對稱結(jié)構(gòu),計算模型的大小取其完整模型的1/16。雙穿孔管消聲器幾何模型中1為消聲器進口管,2為膨脹腔,3為消聲器出口管,空氣溫度T=288 K。

    表1 兩種消聲器的具體尺寸Table 1 Specific dimensions of the two mufflers

    圖1 橫流式穿孔管消聲器的幾何模型Fig.1 Geometric model of cross flow perforated tube muffler

    圖2 雙穿孔管消聲器的幾何模型Fig.2 Geometric model of double perforated tube muffler

    2.2 消聲器內(nèi)部流場的計算

    為獲得消聲器內(nèi)各物理量的分布,首先對幾何模型進行流場網(wǎng)格離散,本文中膨脹腔和進出口管網(wǎng)格的大小為2 mm,穿孔處網(wǎng)格的大小為0.5 mm;然后使用FLUENT[10]軟件對消聲器內(nèi)部流場進行計算。采用壓力基隱式求解器用于求解質(zhì)量、動量和能量守恒方程以及狀態(tài)方程組成的耦合方程組,壓力速度-耦合格式選用壓力關(guān)聯(lián)方程(Semi-Implicit Method for pressure Linked Equation,SIMPLE)算法,空間離散格式為二階精度,可實現(xiàn)k-ε兩方程模型用于模擬湍流流動[7-9,15]。消聲器的介質(zhì)為空氣,密度滿足理想氣體定律,其他物性參數(shù)由氣體溫度決定,邊界條件如下:

    (1)消聲器的進口施加速度進口邊界條件;

    (2)消聲器的出口施加壓力出口邊界條件;

    (3)壁面處設(shè)置為靜止、絕熱和無滑移邊界條件;

    (4)消聲器若為軸對稱結(jié)構(gòu),對其內(nèi)部面采用旋轉(zhuǎn)周期邊界條件。

    邊界條件的設(shè)置如表2中所示。Vinlet為消聲器進口空氣的速度,Poutlet為消聲器出口空氣壓力。

    表2 邊界條件的設(shè)置Table 2 Setting of the boundary conditions

    2.3 消聲器傳遞損失的計算

    對于聲學傳遞損失的計算,首先需要對幾何模型進行網(wǎng)格劃分,一般要求聲學網(wǎng)格的大小不超過最大計算頻率所對應(yīng)波長的1/6,本文中最大計算頻率為3 500 Hz,故離散網(wǎng)格的最大尺寸不能超過16 mm。當消聲器的內(nèi)部流場計算完成后,在LMS Virtual Lab聲學有限元模型中通過網(wǎng)格映射的方式將流體速度分布映射到聲學有限元模型中。在LMS Virtual Lab中映射算法主要包括節(jié)點編號法、最大距離算法、最大單元距離算法和連續(xù)最大單元距離算法。本文使用的是連續(xù)最大單元距離算法,因為該算法適合于流場精細CFD網(wǎng)格向聲學粗糙網(wǎng)格映射,可以保證大范圍的流場尺度的連續(xù)性[11]。消聲器進口設(shè)置為質(zhì)點振速邊界條件,出口為無反射端邊界條件。聲學有限元模型中,不劃分穿孔處的網(wǎng)格,因為在穿孔處由于流體的黏性作用對聲能的耗散尤為嚴重,聲波與流體流動的相互作用在小孔處會產(chǎn)生脫落渦,脫落渦將聲能轉(zhuǎn)化成熱能,加劇了聲能的耗散[9,15-16]。文獻[17]中對直通穿孔管消聲器建出了穿孔處的網(wǎng)格,但是發(fā)現(xiàn)計算值與實驗值相比偏低,這是因為聲學有限元計算過程中不包括流體的黏性作用,無法考慮穿孔處脫落渦對聲能的耗散。因此在消聲器的穿孔處不建立網(wǎng)格,而是通過定義傳遞導納關(guān)系來建立消聲器的穿孔阻抗模型,這樣定義穿孔阻抗就考慮到了脫落渦對聲能的耗散。為了選取最佳的穿孔阻抗模型,消聲器的穿孔阻抗模型分別采用Sullivan模型[13]和康鐘緒模型[18],Sullivan穿孔阻抗模型為

    式中:zS1為無流條件下的穿孔阻抗,zS2為有流條件下的穿孔阻抗,k為聲波數(shù),tw為穿孔管的厚度,dh為小孔的直徑,U0為小孔中的流速,ρ為流體的密度,cc為流體中的聲速,φ為穿孔管的穿孔率。康鐘緒無流條件下穿孔阻抗模型[18]為

    康鐘緒有流條件下的穿孔阻抗模型為

    式(3)和式(4)中,ρ為流體的密度,cc為流體中的聲速,k為聲波數(shù),tw為穿孔管的厚度,α為端部修正系數(shù),dh為小孔的直徑,φ為穿孔管的穿孔率,MB為孔內(nèi)氣流馬赫數(shù)。其余壁面設(shè)置為剛性壁面邊界條件,然后進行聲學響應(yīng)的計算,從而計算得到消聲器的傳遞損失。

    3 計算結(jié)果及分析

    3.1 橫流式穿孔管消聲器結(jié)果及分析

    圖3~5分別為橫流式穿孔管消聲器的壓力、速度、流線分布云圖??梢钥闯?,進口管、膨脹腔、出口管中的壓力依次降低,而消聲器內(nèi)的速度分布并不均勻,小孔附近的速度較高,進口管和出口管中的速度相對較低,膨脹腔中的速度最??;由于小孔的作用,在膨脹腔以及出口管中形成了射流,每個小孔附近的射流速度并不完全相同。由圖5可以看出,進口管處的膨脹腔中流線分布比較復雜,相互交錯,并且與壁面產(chǎn)生了碰撞,形成了很多的漩渦,消耗的能量較多,而出口管處膨脹腔中的流線分布比較光滑,漩渦較少,能量的消耗較少。

    圖3 橫流式穿孔管消聲器靜壓分布云圖Fig.3 Nephogram of static pressure of cross-flow perforatedtube muffler

    圖4 橫流式穿孔管消聲器速度分布云圖Fig.4 Nephogram of velocity magnitude of cross-flow perforated-tube muffler

    圖5 橫流式穿孔管消聲器流線分布云圖Fig.5 Nephogram of streamline of cross-flow perforated-tube muffler

    在得到橫流式穿孔管消聲器流場的計算結(jié)果以后,將流體速度、壓力、密度、溫度等物理量通過網(wǎng)格映射的方法傳遞到聲學有限元模型中,即可計算有流情況下消聲器的傳遞損失。圖6為使用Sullivan穿孔阻抗模型和康鐘緒穿孔阻抗模型計算得到的橫流式穿孔管傳遞損失與實驗測量結(jié)果的比較??梢钥闯?,與實驗值相比較,頻率1 000 Hz以下采用Sullivan穿孔阻抗模型得到的傳遞損較好,頻率1 000~2 800 Hz采用康鐘緒穿孔阻抗模型得到的傳遞損失較好,但在頻率2 800~3 500 Hz范圍內(nèi)兩者與實驗值都有較大差異,可見利用不同的穿孔阻抗模型計算的傳遞損失并不相同。有限元計算值與實驗測量值產(chǎn)生差異的原因可以歸結(jié)如下:(1)實驗測量中可能存在流體與壁面之間的熱傳導,而在有限元計算過程中并沒有考慮;(2)試驗中流速的測量可能存在誤差[8,15]。

    圖6 橫流式穿孔管消聲器的傳遞損失Fig.6 Transmission loss of cross flow perforated tube muffler

    3.2 雙穿孔管消聲器結(jié)果及分析

    圖7~9分別為雙穿孔管消聲器的壓力、速度、流線分布云圖。從圖7~9中的云圖可以看出,進口管、膨脹腔、出口管中的壓力依次降低,而且在進口管中距離進口越遠,進口管中的壓力越大,在出口管中距離出口越遠,壓力越大;雙穿孔管消聲器中速度分布差異較大,膨脹腔中的流體流速較小,出口管中的速度分布不如進口管中速度分布均勻。由圖9可以看出,由進口管中流出的流體需要經(jīng)過多次碰撞、轉(zhuǎn)折才能進入出口管中,由此造成的能量損失較大。

    圖7 雙穿孔管消聲器壓力分布云圖Fig.7 Nephogram of static pressure of double perforated-tube muffler

    圖8 雙穿孔管消聲器速度分布云圖Fig.8 Nephogram of velocity magnitude of double perforatedtube muffler

    圖9 雙穿孔管消聲器流線分布云圖Fig.9 Nephogram of streamline of double perforated-tube muffler

    圖10比較了雙穿孔管消聲器傳遞損失選取兩種不同穿孔阻抗模型的有限元計算結(jié)果和實驗測量結(jié)果。由圖10可以看出,與實驗值相比,采用Sullivan穿孔阻抗模型得到的傳遞損失與實驗測量結(jié)果整體吻合良好,而康鐘緒穿孔阻抗模型在500~820 Hz頻段與實驗值產(chǎn)生的差異較大,這是因為康鐘緒穿孔阻抗模型只能應(yīng)用于穿孔孔徑為2~5 mm的情況,而文中的雙穿孔管消聲器的穿孔孔徑為10 mm。在共振峰680 Hz處兩種阻抗模型的計算結(jié)果與實驗值都略有差異,這可能是由于實驗測量誤差引起的。

    圖10 雙穿孔管消聲器的傳遞損失Fig.10 Transmission loss of double perforated-tube muffler

    3.3 流速對消聲器傳遞損失的影響

    由前兩節(jié)的分析可知,Sullivan穿孔阻抗模型對于有流情況下消聲器傳遞損失的計算結(jié)果比較準確[19]。因此,采用Sullivan穿孔阻抗模型比較流速對消聲器傳遞損失的影響。

    圖11、12比較了在不同流速下橫流式穿孔管消聲器傳遞損失和雙穿孔管消聲器傳遞損失的計算結(jié)果。由于進口的氣流速度最大為34 m·s-1,消聲器中的氣體流速基本都在0.3馬赫數(shù)以下,由于湍流引起的空氣動力性噪聲可以忽略。由圖11、12可以看出,隨著流速的增加,多數(shù)頻率下消聲器的傳遞損失增大,這是由于流速的增加,穿孔的聲阻增大和聲抗減小所致,但在圖12中可以看出共振峰處隨著流速的增加,共振峰值有所下降。通過以上的分析可知,有流情況下,消聲器的傳遞損失不僅沒有減小,而且有所增加。兩種消聲器在不同流速下的壓力損失如表3所示。但是由表3可以看出,在較高的流速條件下,消聲器進、出口的壓力損失也較大。

    表3 兩種消聲器在不同流速下的壓力損失Table 3 Pressure loss of two mufflers with different inlet velocities

    圖11 流速對橫流式消聲器傳遞損失的影響Fig.11 Effect of flow velocity on TL of cross flow perforated tube muffler

    圖12 流速對雙穿孔管消聲器傳遞損失的影響Fig.12 Effect of flow velocity on TL of double perforated pipe muffler

    4 結(jié)論

    為考慮介質(zhì)流動對消聲器消聲性能的影響,本文使用軟件FLUENT和LMS Virtual Lab聯(lián)合仿真計算有流條件下穿孔管消聲器的傳遞損失。仿真計算的結(jié)果與實驗測量結(jié)果吻合良好,說明了此種方法能夠充分考慮流體流動對消聲器傳遞損失的影響。

    文中選取了不同的穿孔阻抗模型來計算消聲器的傳遞損失,檢驗了穿孔阻抗模型的有效性。分析結(jié)果表明,Sullivan聲阻抗經(jīng)驗公式相對康鐘緒聲阻抗模型更為準確,適用穿孔孔徑范圍更廣。文中還比較了不同流速下消聲器傳遞損失的計算結(jié)果,隨著流速的增加,在絕大多數(shù)頻率下消聲器的傳遞損失增大,但是在共振峰處,共振峰的峰值會減小。此種方法過程簡單,計算速度快,具有較高的工程應(yīng)用價值。

    猜你喜歡
    橫流邊界條件聲學
    橫流熱源塔換熱性能研究
    煤氣與熱力(2021年3期)2021-06-09 06:16:20
    一類帶有Stieltjes積分邊界條件的分數(shù)階微分方程邊值問題正解
    帶有積分邊界條件的奇異攝動邊值問題的漸近解
    愛的就是這股Hi-Fi味 Davis Acoustics(戴維斯聲學)Balthus 70
    Acoustical Treatment Primer:Diffusion談?wù)劼晫W處理中的“擴散”
    Acoustical Treatment Primer:Absorption談?wù)劼晫W處理中的“吸聲”(二)
    Acoustical Treatment Primer:Absorption 談?wù)劼晫W處理中的“吸聲”
    基于橫流風扇技術(shù)的直升機反扭驗證
    脊下橫流對PEMFC性能影響的數(shù)值分析
    帶Robin邊界條件的2維隨機Ginzburg-Landau方程的吸引子
    激情五月婷婷亚洲| 欧美日韩国产mv在线观看视频| 国产 精品1| tube8黄色片| 麻豆成人av视频| 久久99一区二区三区| 久久久久精品久久久久真实原创| 男女边吃奶边做爰视频| 夜夜骑夜夜射夜夜干| 免费黄网站久久成人精品| 午夜免费男女啪啪视频观看| 久久人妻熟女aⅴ| 搡女人真爽免费视频火全软件| 搡老乐熟女国产| 人妻少妇偷人精品九色| 精品视频人人做人人爽| 久久久久久久久久久久大奶| 黑人高潮一二区| 日韩欧美精品免费久久| 老女人水多毛片| 国产在线免费精品| 3wmmmm亚洲av在线观看| 91精品伊人久久大香线蕉| 狂野欧美激情性bbbbbb| h视频一区二区三区| 国产成人精品久久久久久| 日本vs欧美在线观看视频 | 91精品一卡2卡3卡4卡| 欧美精品一区二区大全| 欧美国产精品一级二级三级 | 国产亚洲5aaaaa淫片| 一级a做视频免费观看| 热99国产精品久久久久久7| 又爽又黄a免费视频| 成人亚洲精品一区在线观看| 美女国产视频在线观看| 美女cb高潮喷水在线观看| 国产乱人偷精品视频| 成人亚洲精品一区在线观看| 我的女老师完整版在线观看| 午夜91福利影院| 纯流量卡能插随身wifi吗| 亚洲无线观看免费| 久久国产乱子免费精品| 久久韩国三级中文字幕| 精品人妻一区二区三区麻豆| 久久久国产精品麻豆| 美女主播在线视频| 亚洲色图综合在线观看| 伊人久久国产一区二区| 国产免费福利视频在线观看| 国国产精品蜜臀av免费| 日韩 亚洲 欧美在线| 精品亚洲成a人片在线观看| 国产毛片在线视频| 国产毛片在线视频| 大香蕉97超碰在线| 久久久久久人妻| 亚洲精品456在线播放app| 又黄又爽又刺激的免费视频.| 欧美高清成人免费视频www| 色94色欧美一区二区| 亚洲一区二区三区欧美精品| 大码成人一级视频| 国语对白做爰xxxⅹ性视频网站| 亚洲成人一二三区av| 日韩熟女老妇一区二区性免费视频| 国产爽快片一区二区三区| 国产欧美日韩一区二区三区在线 | 亚洲成人一二三区av| 国产一区有黄有色的免费视频| 成人午夜精彩视频在线观看| 欧美变态另类bdsm刘玥| 国产一区二区三区综合在线观看 | 免费黄色在线免费观看| 国产精品久久久久久精品古装| 国产日韩欧美亚洲二区| 国产精品不卡视频一区二区| 国产成人一区二区在线| 免费av不卡在线播放| 中文精品一卡2卡3卡4更新| 久久99热这里只频精品6学生| 悠悠久久av| 成年人免费黄色播放视频| 一区在线观看完整版| 一区二区三区四区激情视频| 无遮挡黄片免费观看| 国产精品国产av在线观看| 五月开心婷婷网| 国产欧美日韩综合在线一区二区| 国产成+人综合+亚洲专区| 不卡一级毛片| 国产精品自产拍在线观看55亚洲 | 性色av乱码一区二区三区2| 亚洲成国产人片在线观看| 五月天丁香电影| 国产精品香港三级国产av潘金莲| a级毛片在线看网站| 欧美大码av| 久久久久久免费高清国产稀缺| 亚洲男人天堂网一区| 欧美日韩视频精品一区| 一区二区三区四区激情视频| 精品亚洲成a人片在线观看| 成人影院久久| 久久久久网色| 一进一出抽搐动态| 日韩制服丝袜自拍偷拍| 久久精品久久久久久噜噜老黄| 大陆偷拍与自拍| 老熟妇仑乱视频hdxx| 多毛熟女@视频| 最新在线观看一区二区三区| 亚洲av国产av综合av卡| 午夜福利一区二区在线看| 黄色视频在线播放观看不卡| 脱女人内裤的视频| 亚洲第一欧美日韩一区二区三区 | 91九色精品人成在线观看| 精品第一国产精品| 人妻久久中文字幕网| av福利片在线| 日韩中文字幕视频在线看片| av福利片在线| 18禁黄网站禁片午夜丰满| 热99re8久久精品国产| 久久中文字幕一级| 波多野结衣一区麻豆| 国产又色又爽无遮挡免| 人人妻,人人澡人人爽秒播| 美女视频免费永久观看网站| 亚洲成人免费电影在线观看| 亚洲精品在线美女| 999久久久国产精品视频| 波多野结衣av一区二区av| 国产精品影院久久| 亚洲少妇的诱惑av| 国产亚洲精品第一综合不卡| 亚洲精品一区蜜桃| 下体分泌物呈黄色| 女性被躁到高潮视频| 女性被躁到高潮视频| 99热全是精品| 午夜精品国产一区二区电影| 亚洲专区中文字幕在线| 国产成人欧美| 成人av一区二区三区在线看 | 岛国毛片在线播放| 蜜桃在线观看..| 日本黄色日本黄色录像| 91成年电影在线观看| 亚洲性夜色夜夜综合| 精品国产一区二区三区久久久樱花| 亚洲av电影在线进入| 不卡av一区二区三区| 飞空精品影院首页| 国产亚洲av片在线观看秒播厂| 亚洲精品美女久久久久99蜜臀| 美女大奶头黄色视频| 色婷婷久久久亚洲欧美| av不卡在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲综合色网址| 正在播放国产对白刺激| 精品欧美一区二区三区在线| 国产高清videossex| 亚洲欧美日韩高清在线视频 | 精品国内亚洲2022精品成人 | 日韩欧美国产一区二区入口| 亚洲专区中文字幕在线| 丁香六月天网| kizo精华| 国产欧美日韩精品亚洲av| 香蕉丝袜av| 十八禁人妻一区二区| 青草久久国产| 国产成人av教育| 亚洲国产中文字幕在线视频| 国产成人免费无遮挡视频| 一本大道久久a久久精品| 一级黄色大片毛片| 激情视频va一区二区三区| 在线亚洲精品国产二区图片欧美| av线在线观看网站| 午夜福利视频在线观看免费| 91字幕亚洲| 午夜免费成人在线视频| av欧美777| 成人免费观看视频高清| 久久中文字幕一级| 最近最新中文字幕大全免费视频| 男女之事视频高清在线观看| 久久精品亚洲av国产电影网| 精品熟女少妇八av免费久了| 日韩制服骚丝袜av| 国产成人啪精品午夜网站| 精品少妇久久久久久888优播| 国产亚洲精品一区二区www | 午夜福利免费观看在线| 国产亚洲精品久久久久5区| 一级片免费观看大全| 国产国语露脸激情在线看| 一级片'在线观看视频| 免费久久久久久久精品成人欧美视频| 淫妇啪啪啪对白视频 | 另类精品久久| 精品国产一区二区三区四区第35| 中文字幕人妻丝袜一区二区| 国产精品.久久久| 久久精品国产亚洲av高清一级| 欧美亚洲 丝袜 人妻 在线| 一区二区av电影网| 欧美另类亚洲清纯唯美| 国产99久久九九免费精品| 国产精品久久久久成人av| 日韩中文字幕欧美一区二区| 超碰97精品在线观看| 久久久久久人人人人人| 操美女的视频在线观看| 俄罗斯特黄特色一大片| 又大又爽又粗| 亚洲色图 男人天堂 中文字幕| 热99re8久久精品国产| 精品久久蜜臀av无| 欧美大码av| 日韩免费高清中文字幕av| 动漫黄色视频在线观看| www.自偷自拍.com| 天天影视国产精品| 亚洲国产精品999| 深夜精品福利| 黑丝袜美女国产一区| 亚洲黑人精品在线| 制服诱惑二区| 国产一区二区在线观看av| 精品一区在线观看国产| 午夜福利影视在线免费观看| 女性被躁到高潮视频| 国产一区有黄有色的免费视频| 男女午夜视频在线观看| 亚洲av美国av| 国产成人精品无人区| 麻豆国产av国片精品| 熟女少妇亚洲综合色aaa.| 黑人操中国人逼视频| 黄网站色视频无遮挡免费观看| 黄色 视频免费看| 男女无遮挡免费网站观看| 精品人妻熟女毛片av久久网站| 色94色欧美一区二区| 日韩人妻精品一区2区三区| 十八禁人妻一区二区| 两人在一起打扑克的视频| 在线看a的网站| 伊人亚洲综合成人网| 亚洲欧美一区二区三区黑人| 青草久久国产| 亚洲avbb在线观看| 亚洲av片天天在线观看| 动漫黄色视频在线观看| 免费少妇av软件| 亚洲精品美女久久av网站| 男女免费视频国产| 各种免费的搞黄视频| 一区二区三区四区激情视频| 麻豆乱淫一区二区| 国产成人系列免费观看| 日韩熟女老妇一区二区性免费视频| 日韩人妻精品一区2区三区| 成年女人毛片免费观看观看9 | 亚洲午夜精品一区,二区,三区| 日韩欧美国产一区二区入口| 精品人妻熟女毛片av久久网站| 永久免费av网站大全| 精品熟女少妇八av免费久了| 中国国产av一级| 国产精品熟女久久久久浪| 丰满少妇做爰视频| 超碰97精品在线观看| 久久国产精品男人的天堂亚洲| 国产成人欧美| 欧美变态另类bdsm刘玥| 黄频高清免费视频| 亚洲人成电影免费在线| 亚洲欧洲精品一区二区精品久久久| 窝窝影院91人妻| 午夜91福利影院| 国产精品国产三级国产专区5o| 欧美+亚洲+日韩+国产| 亚洲成人国产一区在线观看| 国内毛片毛片毛片毛片毛片| 亚洲国产av新网站| 黄色怎么调成土黄色| 狂野欧美激情性bbbbbb| 日本91视频免费播放| 国产免费现黄频在线看| 青草久久国产| 国产成人一区二区三区免费视频网站| 国产高清videossex| 老熟妇乱子伦视频在线观看 | 亚洲精品粉嫩美女一区| 日本av免费视频播放| 国产男人的电影天堂91| av在线播放精品| 一区二区三区激情视频| 中文字幕制服av| 亚洲久久久国产精品| 国产成人影院久久av| 欧美大码av| 又黄又粗又硬又大视频| 韩国精品一区二区三区| www.av在线官网国产| 国产精品偷伦视频观看了| 欧美少妇被猛烈插入视频| 咕卡用的链子| 国产在线一区二区三区精| 精品乱码久久久久久99久播| 80岁老熟妇乱子伦牲交| 黑丝袜美女国产一区| 午夜老司机福利片| 亚洲精品久久午夜乱码| 一级a爱视频在线免费观看| 欧美97在线视频| 欧美日韩中文字幕国产精品一区二区三区 | 黄色视频,在线免费观看| 亚洲国产av影院在线观看| 精品国产乱子伦一区二区三区 | 成年动漫av网址| 亚洲中文av在线| 亚洲人成电影免费在线| 正在播放国产对白刺激| 日韩 亚洲 欧美在线| 日韩精品免费视频一区二区三区| 伊人亚洲综合成人网| 777米奇影视久久| 久久国产亚洲av麻豆专区| 999精品在线视频| 黄色视频不卡| 别揉我奶头~嗯~啊~动态视频 | 亚洲伊人色综图| 国产一区二区激情短视频 | 亚洲一区中文字幕在线| 波多野结衣av一区二区av| 亚洲人成77777在线视频| 性少妇av在线| 黑人操中国人逼视频| 国产精品久久久av美女十八| 午夜福利免费观看在线| 女警被强在线播放| 国产成人免费无遮挡视频| 建设人人有责人人尽责人人享有的| 人妻 亚洲 视频| 在线看a的网站| 国产97色在线日韩免费| netflix在线观看网站| 999久久久国产精品视频| 欧美日韩中文字幕国产精品一区二区三区 | 中国美女看黄片| 国产精品一区二区免费欧美 | 亚洲欧美精品自产自拍| 国产精品亚洲av一区麻豆| 精品国产国语对白av| 精品亚洲成国产av| 免费女性裸体啪啪无遮挡网站| 国产区一区二久久| 久久久久网色| 美女午夜性视频免费| 成年美女黄网站色视频大全免费| 91麻豆av在线| 老司机在亚洲福利影院| 考比视频在线观看| 亚洲国产精品一区二区三区在线| 大码成人一级视频| 亚洲欧美一区二区三区久久| 无遮挡黄片免费观看| 男女免费视频国产| 国产精品欧美亚洲77777| 又紧又爽又黄一区二区| 一二三四社区在线视频社区8| 十八禁人妻一区二区| 制服人妻中文乱码| 老汉色∧v一级毛片| 国产免费现黄频在线看| 中文字幕色久视频| 一本久久精品| 搡老岳熟女国产| 日本a在线网址| 成年女人毛片免费观看观看9 | 久久精品国产a三级三级三级| 亚洲精品美女久久av网站| 亚洲色图 男人天堂 中文字幕| 亚洲色图综合在线观看| 日本91视频免费播放| 免费看十八禁软件| 真人做人爱边吃奶动态| 久久精品国产亚洲av香蕉五月 | 手机成人av网站| 日韩 欧美 亚洲 中文字幕| 永久免费av网站大全| 色视频在线一区二区三区| 极品少妇高潮喷水抽搐| 黄频高清免费视频| 丰满迷人的少妇在线观看| 婷婷色av中文字幕| 亚洲性夜色夜夜综合| 久久久久精品国产欧美久久久 | 免费在线观看完整版高清| 最黄视频免费看| 少妇猛男粗大的猛烈进出视频| 国产深夜福利视频在线观看| 视频区欧美日本亚洲| 日本一区二区免费在线视频| 欧美日韩国产mv在线观看视频| 婷婷色av中文字幕| 精品少妇黑人巨大在线播放| 国产极品粉嫩免费观看在线| 老司机深夜福利视频在线观看 | 欧美日韩视频精品一区| www日本在线高清视频| 91精品三级在线观看| 国产成人免费观看mmmm| 侵犯人妻中文字幕一二三四区| 一二三四社区在线视频社区8| 欧美日韩福利视频一区二区| 午夜日韩欧美国产| 国产精品麻豆人妻色哟哟久久| 免费高清在线观看视频在线观看| 91精品国产国语对白视频| 久久狼人影院| a级毛片在线看网站| 中文字幕高清在线视频| 啦啦啦免费观看视频1| 日韩有码中文字幕| 国产精品久久久久久精品电影小说| 免费在线观看完整版高清| 国产区一区二久久| 久久久久国产精品人妻一区二区| 欧美老熟妇乱子伦牲交| 18禁黄网站禁片午夜丰满| 国产精品一区二区精品视频观看| 亚洲成av片中文字幕在线观看| 亚洲精品乱久久久久久| 12—13女人毛片做爰片一| 少妇裸体淫交视频免费看高清 | 国产精品成人在线| 亚洲精品久久午夜乱码| 宅男免费午夜| 熟女少妇亚洲综合色aaa.| 丝瓜视频免费看黄片| 欧美精品人与动牲交sv欧美| 亚洲色图综合在线观看| 午夜福利乱码中文字幕| 一区在线观看完整版| 黄色视频在线播放观看不卡| 女性生殖器流出的白浆| 国产成人啪精品午夜网站| av线在线观看网站| 成年人黄色毛片网站| 岛国在线观看网站| 精品人妻熟女毛片av久久网站| 中文精品一卡2卡3卡4更新| 国产视频一区二区在线看| 国产成人精品在线电影| 青春草亚洲视频在线观看| 免费少妇av软件| 国产国语露脸激情在线看| 午夜激情av网站| 成在线人永久免费视频| 多毛熟女@视频| 日本五十路高清| 国产av又大| a级毛片黄视频| 亚洲精品粉嫩美女一区| 亚洲欧洲精品一区二区精品久久久| 亚洲成人免费av在线播放| 免费在线观看黄色视频的| 国产亚洲av高清不卡| 天堂中文最新版在线下载| 久久青草综合色| 午夜福利乱码中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 侵犯人妻中文字幕一二三四区| 国产伦人伦偷精品视频| 精品国产乱子伦一区二区三区 | 国产成人一区二区三区免费视频网站| 黑人操中国人逼视频| videos熟女内射| 国产精品久久久久久人妻精品电影 | 80岁老熟妇乱子伦牲交| 国产区一区二久久| 亚洲 国产 在线| 久久性视频一级片| 水蜜桃什么品种好| 黑丝袜美女国产一区| 久久99热这里只频精品6学生| 国产精品香港三级国产av潘金莲| 色婷婷av一区二区三区视频| 亚洲精品国产区一区二| 美女大奶头黄色视频| 国产在视频线精品| videosex国产| 日本av手机在线免费观看| 热99久久久久精品小说推荐| 亚洲精品国产一区二区精华液| 人人妻,人人澡人人爽秒播| 国产熟女午夜一区二区三区| 少妇 在线观看| 国产成人精品无人区| 咕卡用的链子| videosex国产| 日本av手机在线免费观看| 国产免费av片在线观看野外av| 大码成人一级视频| 自线自在国产av| 青草久久国产| 一边摸一边抽搐一进一出视频| 国产男女超爽视频在线观看| 亚洲欧美日韩高清在线视频 | 青春草视频在线免费观看| 国产欧美日韩综合在线一区二区| 亚洲五月色婷婷综合| 男女下面插进去视频免费观看| 无限看片的www在线观看| 性高湖久久久久久久久免费观看| 黄色怎么调成土黄色| 日韩欧美一区视频在线观看| 91麻豆av在线| av线在线观看网站| 国产欧美日韩一区二区精品| 一边摸一边做爽爽视频免费| 纵有疾风起免费观看全集完整版| 91九色精品人成在线观看| 久久久国产精品麻豆| 亚洲 欧美一区二区三区| 宅男免费午夜| 欧美少妇被猛烈插入视频| 美女中出高潮动态图| av又黄又爽大尺度在线免费看| 亚洲成人免费av在线播放| 一区二区av电影网| 亚洲 国产 在线| 美女主播在线视频| 国产亚洲一区二区精品| 男女边摸边吃奶| 午夜激情av网站| 永久免费av网站大全| 欧美国产精品va在线观看不卡| 精品久久久久久久毛片微露脸 | 9热在线视频观看99| 亚洲av日韩在线播放| 国产精品久久久久成人av| 91成年电影在线观看| 丝袜美足系列| 成年人免费黄色播放视频| 日日夜夜操网爽| 女性生殖器流出的白浆| 日日摸夜夜添夜夜添小说| 99re6热这里在线精品视频| 精品第一国产精品| 亚洲av美国av| 国产精品一区二区免费欧美 | 法律面前人人平等表现在哪些方面 | 亚洲国产日韩一区二区| 天堂8中文在线网| 人妻一区二区av| 久久国产精品人妻蜜桃| 欧美午夜高清在线| 久久热在线av| 黄片小视频在线播放| 成人18禁高潮啪啪吃奶动态图| 他把我摸到了高潮在线观看 | 欧美一级毛片孕妇| 美女午夜性视频免费| 国产亚洲欧美精品永久| 久久精品aⅴ一区二区三区四区| 男女高潮啪啪啪动态图| 久久久欧美国产精品| 亚洲欧洲精品一区二区精品久久久| 国产免费av片在线观看野外av| 国产片内射在线| 蜜桃在线观看..| 亚洲国产精品一区三区| 黄片小视频在线播放| 精品第一国产精品| 久久av网站| 人人澡人人妻人| 一二三四社区在线视频社区8| 久久综合国产亚洲精品| 精品人妻熟女毛片av久久网站| 日日夜夜操网爽| 岛国毛片在线播放| 精品少妇内射三级| 久久人妻福利社区极品人妻图片| 咕卡用的链子| 亚洲精品中文字幕一二三四区 | 老司机深夜福利视频在线观看 | 亚洲av成人一区二区三| 成在线人永久免费视频| 精品视频人人做人人爽| 夜夜夜夜夜久久久久| 午夜福利视频在线观看免费| 色综合欧美亚洲国产小说| 久热这里只有精品99| 两个人看的免费小视频| 日本黄色日本黄色录像| 国产免费视频播放在线视频| 亚洲av电影在线进入| 亚洲专区字幕在线| 纵有疾风起免费观看全集完整版| 水蜜桃什么品种好| 91国产中文字幕| 婷婷成人精品国产| 99热国产这里只有精品6| 女人高潮潮喷娇喘18禁视频| 亚洲成国产人片在线观看|