• <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方程的吸引子
    亚洲四区av| 国内精品宾馆在线| 久久精品国产自在天天线| 少妇高潮的动态图| 亚洲精品456在线播放app| 亚洲欧洲日产国产| 欧美精品人与动牲交sv欧美| 国产 一区 欧美 日韩| 亚洲国产精品专区欧美| 成人午夜精彩视频在线观看| 久久 成人 亚洲| 婷婷色综合大香蕉| 久久久久久久久久久免费av| 国产日韩欧美亚洲二区| 国产高清国产精品国产三级 | 亚洲精品一二三| 男人爽女人下面视频在线观看| 黑人高潮一二区| 这个男人来自地球电影免费观看 | 国产精品嫩草影院av在线观看| 久久久久精品久久久久真实原创| 人人妻人人澡人人爽人人夜夜| 纯流量卡能插随身wifi吗| 一区二区三区四区激情视频| 婷婷色综合大香蕉| 久久女婷五月综合色啪小说| 中文字幕亚洲精品专区| 国产精品精品国产色婷婷| 最近最新中文字幕免费大全7| 日韩人妻高清精品专区| 直男gayav资源| 久久久久网色| 国产视频首页在线观看| 日韩av不卡免费在线播放| 97超碰精品成人国产| 国产成人免费无遮挡视频| 日韩伦理黄色片| 看免费成人av毛片| 免费看av在线观看网站| 国产免费视频播放在线视频| 国产永久视频网站| 欧美成人精品欧美一级黄| 综合色丁香网| 亚洲国产精品专区欧美| av网站免费在线观看视频| 蜜桃亚洲精品一区二区三区| 亚洲色图综合在线观看| 美女脱内裤让男人舔精品视频| 日韩一区二区视频免费看| 最近2019中文字幕mv第一页| 国产男人的电影天堂91| 成人免费观看视频高清| 在线亚洲精品国产二区图片欧美 | 久久精品熟女亚洲av麻豆精品| 免费观看的影片在线观看| 国产男人的电影天堂91| 久久久亚洲精品成人影院| 日韩成人伦理影院| 建设人人有责人人尽责人人享有的 | 日韩,欧美,国产一区二区三区| 午夜激情福利司机影院| 亚洲欧洲日产国产| 国产真实伦视频高清在线观看| 国产探花极品一区二区| 亚洲三级黄色毛片| 韩国高清视频一区二区三区| 我要看日韩黄色一级片| 国产精品偷伦视频观看了| 99精国产麻豆久久婷婷| 国产日韩欧美在线精品| 久久精品久久久久久久性| 久久人人爽av亚洲精品天堂 | 国产真实伦视频高清在线观看| 高清视频免费观看一区二区| 最近最新中文字幕大全电影3| 国产精品伦人一区二区| 超碰97精品在线观看| 欧美3d第一页| 晚上一个人看的免费电影| 欧美 日韩 精品 国产| 欧美日韩亚洲高清精品| 国产亚洲一区二区精品| 不卡视频在线观看欧美| 久久精品久久久久久久性| 五月天丁香电影| 午夜免费鲁丝| 99九九线精品视频在线观看视频| 少妇裸体淫交视频免费看高清| 亚洲一级一片aⅴ在线观看| 亚洲熟女精品中文字幕| 天天躁夜夜躁狠狠久久av| 国产黄片美女视频| 亚洲第一av免费看| 免费黄频网站在线观看国产| 极品教师在线视频| 丝袜脚勾引网站| 春色校园在线视频观看| 亚洲av电影在线观看一区二区三区| av又黄又爽大尺度在线免费看| 亚洲av中文av极速乱| 国产精品久久久久久久电影| 看免费成人av毛片| 欧美区成人在线视频| 国产精品一区www在线观看| 女性被躁到高潮视频| av国产久精品久网站免费入址| 日韩成人av中文字幕在线观看| 一级爰片在线观看| 中文字幕av成人在线电影| 精品一区二区三卡| 男人爽女人下面视频在线观看| 国产乱人视频| 大香蕉久久网| 熟女av电影| 亚洲三级黄色毛片| 最近最新中文字幕大全电影3| 亚洲人成网站在线播| 建设人人有责人人尽责人人享有的 | 女人十人毛片免费观看3o分钟| 最近最新中文字幕免费大全7| 亚洲欧美日韩无卡精品| 精品熟女少妇av免费看| 老女人水多毛片| 熟女人妻精品中文字幕| 久久6这里有精品| 少妇的逼好多水| 欧美精品一区二区免费开放| 在线天堂最新版资源| 少妇的逼水好多| 国产大屁股一区二区在线视频| 成人二区视频| 91精品国产国语对白视频| 我要看日韩黄色一级片| 男女边摸边吃奶| 亚洲成人av在线免费| 内射极品少妇av片p| 色5月婷婷丁香| 欧美老熟妇乱子伦牲交| 国产色婷婷99| 久久久久久久大尺度免费视频| 亚洲欧美中文字幕日韩二区| 亚洲欧美中文字幕日韩二区| 亚洲av免费高清在线观看| 国产淫片久久久久久久久| 99久久精品国产国产毛片| 亚洲精品日本国产第一区| 插阴视频在线观看视频| 91久久精品国产一区二区成人| 2018国产大陆天天弄谢| 欧美bdsm另类| 女性生殖器流出的白浆| 欧美人与善性xxx| av在线老鸭窝| 国产精品国产av在线观看| 亚洲精华国产精华液的使用体验| 国产一区二区三区综合在线观看 | 99久久人妻综合| 亚洲人成网站高清观看| 尾随美女入室| 黄色日韩在线| 亚洲激情五月婷婷啪啪| 日韩精品有码人妻一区| av又黄又爽大尺度在线免费看| 综合色丁香网| 在现免费观看毛片| 久久久久国产精品人妻一区二区| 中文字幕亚洲精品专区| 久久婷婷青草| 精品少妇黑人巨大在线播放| 久久人妻熟女aⅴ| 亚洲欧美精品专区久久| 最近最新中文字幕免费大全7| 久久久久久久精品精品| 国产精品嫩草影院av在线观看| 99热6这里只有精品| 女的被弄到高潮叫床怎么办| videos熟女内射| 日韩一区二区视频免费看| 女性生殖器流出的白浆| 欧美精品一区二区大全| 深夜a级毛片| 在线观看三级黄色| 一二三四中文在线观看免费高清| 男人爽女人下面视频在线观看| 美女视频免费永久观看网站| 亚洲人成网站高清观看| 亚洲性久久影院| 最近中文字幕2019免费版| 日本与韩国留学比较| 久久久久久久久大av| 美女xxoo啪啪120秒动态图| 熟女人妻精品中文字幕| 特大巨黑吊av在线直播| 亚洲丝袜综合中文字幕| 美女脱内裤让男人舔精品视频| av在线观看视频网站免费| 亚洲国产欧美在线一区| 久久综合国产亚洲精品| 亚洲精品国产av成人精品| 晚上一个人看的免费电影| 国产淫语在线视频| 网址你懂的国产日韩在线| 精品人妻偷拍中文字幕| 亚州av有码| av国产精品久久久久影院| 国产淫片久久久久久久久| 亚洲av综合色区一区| 男人添女人高潮全过程视频| 黄色配什么色好看| 欧美日韩综合久久久久久| 人妻夜夜爽99麻豆av| 精品国产乱码久久久久久小说| 亚洲国产欧美在线一区| 女人十人毛片免费观看3o分钟| 久久久久国产精品人妻一区二区| 晚上一个人看的免费电影| 成人午夜精彩视频在线观看| 人体艺术视频欧美日本| 91aial.com中文字幕在线观看| 亚洲欧洲国产日韩| 黑人猛操日本美女一级片| 热99国产精品久久久久久7| 国产v大片淫在线免费观看| 亚洲人成网站高清观看| 亚洲性久久影院| 国产成人精品久久久久久| 欧美日韩在线观看h| 亚洲精品乱久久久久久| 国产综合精华液| 亚洲国产最新在线播放| 妹子高潮喷水视频| 日韩电影二区| 亚洲成色77777| 汤姆久久久久久久影院中文字幕| 成人国产av品久久久| 日日摸夜夜添夜夜添av毛片| 性高湖久久久久久久久免费观看| 亚洲激情五月婷婷啪啪| 亚洲av男天堂| 黄色怎么调成土黄色| 99久久精品国产国产毛片| 在线亚洲精品国产二区图片欧美 | 少妇丰满av| 一级片'在线观看视频| 亚洲高清免费不卡视频| 久久99热6这里只有精品| 亚洲av成人精品一区久久| 人妻一区二区av| 天堂中文最新版在线下载| 亚洲av男天堂| 国产成人a∨麻豆精品| 妹子高潮喷水视频| 人人妻人人看人人澡| 亚洲av日韩在线播放| 菩萨蛮人人尽说江南好唐韦庄| 校园人妻丝袜中文字幕| 欧美精品亚洲一区二区| 国产淫语在线视频| 亚洲欧美精品自产自拍| 爱豆传媒免费全集在线观看| 天美传媒精品一区二区| 干丝袜人妻中文字幕| 国产精品国产av在线观看| 亚洲电影在线观看av| av免费在线看不卡| 国产精品久久久久久av不卡| 欧美3d第一页| 人人妻人人爽人人添夜夜欢视频 | 欧美老熟妇乱子伦牲交| 欧美三级亚洲精品| 久久久久久久久大av| 国模一区二区三区四区视频| 狠狠精品人妻久久久久久综合| 麻豆成人av视频| 自拍偷自拍亚洲精品老妇| 亚洲欧洲日产国产| 亚洲国产色片| 欧美一级a爱片免费观看看| 国产成人精品久久久久久| 七月丁香在线播放| 六月丁香七月| 国产乱来视频区| 久热这里只有精品99| 国产精品一区二区性色av| 色吧在线观看| 美女xxoo啪啪120秒动态图| 国产片特级美女逼逼视频| 久久久久久久亚洲中文字幕| 永久网站在线| 少妇人妻久久综合中文| av在线播放精品| 天堂俺去俺来也www色官网| 人人妻人人爽人人添夜夜欢视频 | 免费看不卡的av| 亚洲av日韩在线播放| 国产男女内射视频| 久久久色成人| 中文在线观看免费www的网站| 全区人妻精品视频| 99久久中文字幕三级久久日本| 91在线精品国自产拍蜜月| 如何舔出高潮| 十八禁网站网址无遮挡 | 亚洲精品乱码久久久久久按摩| 久久久成人免费电影| 欧美zozozo另类| 内地一区二区视频在线| 成人亚洲欧美一区二区av| 久久久久久久久大av| 18禁裸乳无遮挡免费网站照片| 一个人免费看片子| 国产av码专区亚洲av| 一本色道久久久久久精品综合| 天天躁日日操中文字幕| 狠狠精品人妻久久久久久综合| 岛国毛片在线播放| 亚洲怡红院男人天堂| 成人无遮挡网站| 99热这里只有是精品50| 一级毛片黄色毛片免费观看视频| 亚洲av成人精品一二三区| 嫩草影院入口| 新久久久久国产一级毛片| 国产精品一区二区三区四区免费观看| 精品久久久久久久久亚洲| 日本av手机在线免费观看| 99热这里只有精品一区| 久久国产精品大桥未久av | 免费观看a级毛片全部| 欧美成人a在线观看| 久久久久久久国产电影| 色综合色国产| 亚洲天堂av无毛| 亚洲精品久久午夜乱码| 男人狂女人下面高潮的视频| 亚洲国产日韩一区二区| 99热国产这里只有精品6| 91精品伊人久久大香线蕉| 国产日韩欧美亚洲二区| 我的老师免费观看完整版| 在线观看三级黄色| 麻豆精品久久久久久蜜桃| 国产一区二区三区综合在线观看 | av国产免费在线观看| 免费大片18禁| 免费黄网站久久成人精品| 亚洲成人av在线免费| 我的女老师完整版在线观看| 国产精品久久久久成人av| 国产深夜福利视频在线观看| 国产av精品麻豆| 亚洲av成人精品一二三区| 国语对白做爰xxxⅹ性视频网站| 亚洲美女搞黄在线观看| .国产精品久久| 中国三级夫妇交换| 美女脱内裤让男人舔精品视频| 亚洲av国产av综合av卡| 亚洲色图综合在线观看| 亚洲精品乱码久久久v下载方式| av福利片在线观看| 国产成人91sexporn| 黄色视频在线播放观看不卡| 小蜜桃在线观看免费完整版高清| 国产永久视频网站| 丰满迷人的少妇在线观看| 国产精品伦人一区二区| 看十八女毛片水多多多| 久久久久久久久大av| 日本欧美国产在线视频| 成年女人在线观看亚洲视频| 中文天堂在线官网| 噜噜噜噜噜久久久久久91| 婷婷色av中文字幕| 黄色一级大片看看| 美女国产视频在线观看| 成人无遮挡网站| 亚洲精品一区蜜桃| 亚洲国产高清在线一区二区三| 五月玫瑰六月丁香| 国产片特级美女逼逼视频| 直男gayav资源| 精品熟女少妇av免费看| 2018国产大陆天天弄谢| 色视频www国产| 黄色一级大片看看| 亚洲国产精品专区欧美| 国产片特级美女逼逼视频| 亚洲精品一二三| 国产一区二区三区av在线| 青春草亚洲视频在线观看| 免费黄网站久久成人精品| 亚洲电影在线观看av| 国产亚洲91精品色在线| 最近的中文字幕免费完整| 欧美激情国产日韩精品一区| 又大又黄又爽视频免费| 日本一二三区视频观看| h视频一区二区三区| 亚洲精品亚洲一区二区| 色视频在线一区二区三区| 精品99又大又爽又粗少妇毛片| 免费观看av网站的网址| 精品久久国产蜜桃| 一级黄片播放器| 日韩,欧美,国产一区二区三区| 成人国产麻豆网| 最近手机中文字幕大全| 国产av一区二区精品久久 | 国产精品免费大片| 少妇 在线观看| 国产成人freesex在线| 亚洲自偷自拍三级| 色视频www国产| 高清午夜精品一区二区三区| 久久ye,这里只有精品| 中文字幕制服av| 国产熟女欧美一区二区| 一级毛片电影观看| 晚上一个人看的免费电影| 一级毛片黄色毛片免费观看视频| 日韩强制内射视频| 在线亚洲精品国产二区图片欧美 | 欧美少妇被猛烈插入视频| 2021少妇久久久久久久久久久| 国产成人精品一,二区| 久久久亚洲精品成人影院| 久久6这里有精品| 久久久久久久久久久丰满| www.av在线官网国产| 七月丁香在线播放| 欧美xxxx性猛交bbbb| 国产成人精品福利久久| 国模一区二区三区四区视频| 又黄又爽又刺激的免费视频.| 国产黄频视频在线观看| 丰满少妇做爰视频| 国产精品不卡视频一区二区| 视频区图区小说| 欧美精品亚洲一区二区| 成人高潮视频无遮挡免费网站| 久久久久视频综合| 国产精品一二三区在线看| 国产av精品麻豆| 亚洲av福利一区| 亚洲精品中文字幕在线视频 | 亚洲图色成人| 成年美女黄网站色视频大全免费 | 亚洲最大成人中文| 成人国产av品久久久| 黄片无遮挡物在线观看| 美女内射精品一级片tv| 久久精品熟女亚洲av麻豆精品| 久久99蜜桃精品久久| 五月天丁香电影| 一个人看的www免费观看视频| 18禁在线播放成人免费| 国产av国产精品国产| 亚洲人成网站高清观看| 国产精品无大码| 久久久色成人| 亚洲最大成人中文| 色哟哟·www| 观看免费一级毛片| 免费高清在线观看视频在线观看| 18禁在线无遮挡免费观看视频| 久久av网站| 亚洲精品乱久久久久久| 国产欧美日韩一区二区三区在线 | 久久人人爽人人片av| 国产黄片美女视频| 久久久久人妻精品一区果冻| 伊人久久国产一区二区| 国产亚洲5aaaaa淫片| 国产 精品1| 在线天堂最新版资源| 少妇猛男粗大的猛烈进出视频| 80岁老熟妇乱子伦牲交| 国产精品人妻久久久影院| 精品99又大又爽又粗少妇毛片| 午夜福利在线在线| 在线 av 中文字幕| 国产高清不卡午夜福利| 欧美性感艳星| av视频免费观看在线观看| 99久久中文字幕三级久久日本| 亚洲成人一二三区av| 国产乱人视频| 2021少妇久久久久久久久久久| 久久国产精品大桥未久av | 久久人人爽av亚洲精品天堂 | 国产一区亚洲一区在线观看| 伊人久久国产一区二区| 免费观看a级毛片全部| 国产成人aa在线观看| 在线天堂最新版资源| 一二三四中文在线观看免费高清| 日韩成人伦理影院| 大片电影免费在线观看免费| 国产在线一区二区三区精| 日韩人妻高清精品专区| 丝袜脚勾引网站| 免费观看av网站的网址| 欧美成人a在线观看| 亚洲久久久国产精品| av一本久久久久| 免费看日本二区| 亚洲第一区二区三区不卡| 男男h啪啪无遮挡| 久久精品国产亚洲av涩爱| 国产精品爽爽va在线观看网站| 成人美女网站在线观看视频| 一本—道久久a久久精品蜜桃钙片| 特大巨黑吊av在线直播| av在线观看视频网站免费| 最近的中文字幕免费完整| 国产免费一级a男人的天堂| 日本欧美国产在线视频| 高清午夜精品一区二区三区| 精品久久久久久久末码| 国产成人精品福利久久| 久久综合国产亚洲精品| 在线播放无遮挡| 成人漫画全彩无遮挡| 另类亚洲欧美激情| 日本猛色少妇xxxxx猛交久久| 免费黄频网站在线观看国产| 免费播放大片免费观看视频在线观看| 美女高潮的动态| 亚洲av中文av极速乱| 免费黄色在线免费观看| 欧美xxⅹ黑人| 午夜老司机福利剧场| 午夜免费观看性视频| 啦啦啦在线观看免费高清www| 精品一区二区三卡| 免费不卡的大黄色大毛片视频在线观看| 99热这里只有精品一区| 深夜a级毛片| 91精品国产九色| 精品国产乱码久久久久久小说| 国产精品伦人一区二区| av在线观看视频网站免费| 国产69精品久久久久777片| 欧美亚洲 丝袜 人妻 在线| 国产中年淑女户外野战色| 纵有疾风起免费观看全集完整版| 99久久精品国产国产毛片| av.在线天堂| 久久久色成人| 成年美女黄网站色视频大全免费 | 尾随美女入室| 日韩不卡一区二区三区视频在线| 国产免费一区二区三区四区乱码| 美女高潮的动态| 人妻 亚洲 视频| 最近中文字幕2019免费版| 欧美精品一区二区免费开放| 日韩强制内射视频| 99热这里只有是精品在线观看| 成人毛片60女人毛片免费| 美女cb高潮喷水在线观看| 久久久久久久久久人人人人人人| 亚洲aⅴ乱码一区二区在线播放| 伊人久久精品亚洲午夜| 青春草亚洲视频在线观看| 久久毛片免费看一区二区三区| 精品久久久久久久久av| 久久精品熟女亚洲av麻豆精品| 最近最新中文字幕免费大全7| 久久久久久久亚洲中文字幕| 久久国产乱子免费精品| 丰满少妇做爰视频| 99视频精品全部免费 在线| 日韩成人av中文字幕在线观看| 最近最新中文字幕大全电影3| 亚洲欧美成人综合另类久久久| 久久久久久久国产电影| 日韩中字成人| 久久精品国产a三级三级三级| 视频区图区小说| 99热这里只有精品一区| 国产精品秋霞免费鲁丝片| 成年av动漫网址| 亚洲av男天堂| 国产成人a区在线观看| 国产成人精品婷婷| 青春草视频在线免费观看| 一级a做视频免费观看| 国产淫语在线视频| 国产亚洲午夜精品一区二区久久| 大又大粗又爽又黄少妇毛片口| 少妇的逼好多水| a级一级毛片免费在线观看| 激情 狠狠 欧美| 欧美激情国产日韩精品一区| 久久精品久久精品一区二区三区| 亚洲在久久综合| 亚洲婷婷狠狠爱综合网| 欧美区成人在线视频| 女性被躁到高潮视频| 久久热精品热| 亚洲精品日本国产第一区| 国产乱人视频| 日本黄色片子视频| 观看av在线不卡| 99久久人妻综合| 97超视频在线观看视频| 成人美女网站在线观看视频| 欧美另类一区| 高清视频免费观看一区二区| 国精品久久久久久国模美| 免费大片黄手机在线观看| 久久午夜福利片| 午夜福利在线观看免费完整高清在|