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

    基于Fluent的頂導(dǎo)向套筒單座調(diào)節(jié)閥流場(chǎng)數(shù)值模擬

    2022-08-24 12:19:34張希恒張耀壬
    化工機(jī)械 2022年3期
    關(guān)鍵詞:節(jié)流調(diào)節(jié)閥套筒

    張希恒 張耀壬 王 宇 張 超

    (蘭州理工大學(xué)石油化工學(xué)院)

    調(diào)節(jié)閥作為一種常見的閥門在工業(yè)領(lǐng)域有著廣泛的應(yīng)用,其性能的改進(jìn)對(duì)提高工藝效率和降低能耗具有重要意義。 因此,研究調(diào)節(jié)閥的流量特性具有重要的工業(yè)應(yīng)用價(jià)值。

    近年來,CFD技術(shù)在閥門內(nèi)流場(chǎng)的數(shù)值模擬方 面[1~3]和 結(jié) 構(gòu) 優(yōu) 化 方 面[4~6]均 有 諸 多 應(yīng) 用。 An Young Joon等研究了截止閥的性能, 應(yīng)用流體力學(xué)方法在定常流狀態(tài)下分析截止閥的特性并取得一定的突破[7]。 何世權(quán)等運(yùn)用CFD技術(shù)對(duì)角式調(diào)節(jié)閥內(nèi)部流場(chǎng)進(jìn)行三維模擬,通過對(duì)比分析不同流向下的速度場(chǎng)與壓力場(chǎng)分布,得出流閉型對(duì)提高流量系數(shù)有著積極作用[8]。 陶正良等運(yùn)用數(shù)值模擬和試驗(yàn)方法對(duì)電站調(diào)節(jié)閥內(nèi)流場(chǎng)進(jìn)行研究,該研究可以判斷出流體在短時(shí)間內(nèi)的壓力損失,并且對(duì)其流動(dòng)特性可以進(jìn)行預(yù)測(cè),為調(diào)節(jié)閥的改型提供方向[9]。 潘永成等建立調(diào)節(jié)閥的三維模型, 利用CFD技術(shù)對(duì)調(diào)節(jié)閥內(nèi)流體進(jìn)行數(shù)值模擬,得到不同開度下閥芯所對(duì)應(yīng)的理想流量特性曲線以及流量與壓降的關(guān)系曲線[10]。 在此,筆者以頂導(dǎo)向套筒單座調(diào)節(jié)閥為研究對(duì)象,對(duì)其內(nèi)部三維流場(chǎng)進(jìn)行仿真模擬,根據(jù)壓力云圖、速度云圖和流線圖對(duì)不同開度下內(nèi)流場(chǎng)流動(dòng)特性進(jìn)行分析,最后通過試驗(yàn)驗(yàn)證數(shù)值模擬的準(zhǔn)確性。

    1 研究對(duì)象

    頂導(dǎo)向套筒單座調(diào)節(jié)閥在管道系統(tǒng)中具有重要的控制作用,用于調(diào)節(jié)管路系統(tǒng)重流體的壓力和流量。 流體流經(jīng)套筒閥時(shí),其流場(chǎng)中存在渦流、氣穴等一些復(fù)雜的流動(dòng)現(xiàn)象,這對(duì)套筒閥流場(chǎng)特性的預(yù)估帶來很大困擾。 在以往工程實(shí)踐中,試驗(yàn)是解決這個(gè)問題的主要方法,但試驗(yàn)往往受到尺寸模型、流場(chǎng)干擾、測(cè)量精度及周期過長(zhǎng)等的限制。 為此,筆者以頂導(dǎo)向套筒單座調(diào)節(jié)閥為研究對(duì)象, 運(yùn)用Fluent對(duì)其內(nèi)流道流場(chǎng)進(jìn)行模擬,通過數(shù)值模擬得到其內(nèi)流道的壓力、速度和流線分布規(guī)律,并進(jìn)行流量流阻試驗(yàn),將模擬結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比,以驗(yàn)證該數(shù)值模擬方法的可行性。

    2 數(shù)值模擬

    2.1 標(biāo)準(zhǔn)k-ε模型

    標(biāo)準(zhǔn)k-ε模型是半經(jīng)驗(yàn)公式,其中k方程是精確方程,ε方程是由經(jīng)驗(yàn)公式導(dǎo)出的方程。 k-ε模型假定流場(chǎng)完全是湍流,可忽略分子間的粘性,因而標(biāo)準(zhǔn)k-ε模型適用于完全是湍流的流場(chǎng)。 標(biāo)準(zhǔn)k-ε模型的方程為湍動(dòng)能方程k和擴(kuò)散方程ε[11]:

    其中,ρ是流體密度,ui是時(shí)均速度,μ是動(dòng)量,μt是湍流粘度,Gk是由平均速度梯度引起的湍動(dòng)能k的產(chǎn)生項(xiàng),Gb是由浮力引起的湍動(dòng)能k的產(chǎn)生項(xiàng),YM是可壓縮湍流中脈動(dòng)擴(kuò)張的貢獻(xiàn), 水作為不可壓縮流體YM=0,σk=1.0和σε=1.3分別是與湍動(dòng)能k和耗散率ε對(duì)應(yīng)的Prandtl數(shù),Sk和Sε是用戶定義的源項(xiàng)。 針對(duì)空氣、水等常見介質(zhì)的基本湍流,通過大量的試驗(yàn)總結(jié)出標(biāo)準(zhǔn)k-ε模型經(jīng)驗(yàn)常數(shù)C1ε=1.44、C2ε=1.92、C3ε=0.09。

    2.2 建模與網(wǎng)格劃分

    通過三維建模軟件SolidWorks2018建立頂導(dǎo)向套筒單座調(diào)節(jié)閥的三維模型(圖1),該閥公稱尺寸為40 mm,結(jié)構(gòu)長(zhǎng)度為251 mm,開啟行程為25 mm。

    圖1 調(diào)節(jié)閥結(jié)構(gòu)示意圖

    為了使閥內(nèi)流場(chǎng)流動(dòng)更穩(wěn)定,對(duì)閥兩端進(jìn)行加管處理,進(jìn)口端延長(zhǎng)2D、出口端延長(zhǎng)6D(D為閥門公稱直徑);由于閥門的實(shí)際結(jié)構(gòu)復(fù)雜,對(duì)Fluent模擬計(jì)算結(jié)果的收斂有一定程度的影響,為了便于計(jì)算,將一些對(duì)流體流動(dòng)產(chǎn)生影響較小的結(jié)構(gòu)做適當(dāng)簡(jiǎn)化處理, 建立相對(duì)開度為10%、20%、30%、40%、50%、60%、70%、80%、90%和100%的分析模型, 在ANSYS Workench軟件中運(yùn)用DM(Design Modeler)模塊,將不同開度的分析模塊導(dǎo)入, 利用反向建模的方法抽取其內(nèi)部流道模型(圖2)。

    圖2 頂導(dǎo)向套筒單座調(diào)節(jié)閥100%內(nèi)部流道模型

    運(yùn)用預(yù)處理軟件ICEM對(duì)內(nèi)部流道模型進(jìn)行網(wǎng)格劃分(圖3)。 進(jìn)出口直管段采用六面體結(jié)構(gòu)網(wǎng)格,中間流道部分由于結(jié)構(gòu)不規(guī)則,結(jié)構(gòu)化網(wǎng)格可能會(huì)存在較大的偽擴(kuò)散,造成極大的計(jì)算誤差,甚至是非物理解,而利用非結(jié)構(gòu)網(wǎng)格則不會(huì)存在此類問題;由于閥芯和閥座之間的流動(dòng)間隙較小,因此對(duì)閥芯和閥座之間的流動(dòng)區(qū)域進(jìn)行局部加密,以提高計(jì)算效率的同時(shí)能更好地捕捉流動(dòng)細(xì)節(jié)。

    圖3 調(diào)節(jié)閥內(nèi)部流道網(wǎng)格劃分

    數(shù)值計(jì)算結(jié)果的準(zhǔn)確性不僅與網(wǎng)格質(zhì)量有關(guān),還與網(wǎng)格的數(shù)量密切相關(guān),選擇合理的網(wǎng)格數(shù)量可以在滿足計(jì)算精度的前提下提高計(jì)算效率,節(jié)約計(jì)算資源。 因此,進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證對(duì)選擇合理的網(wǎng)格數(shù)量尤為重要。 現(xiàn)對(duì)頂導(dǎo)向套筒單座調(diào)節(jié)閥50%開度下的模型進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證,結(jié)果見表1、圖4。

    表1 網(wǎng)格無關(guān)性驗(yàn)證結(jié)果

    圖4 網(wǎng)格無關(guān)性驗(yàn)證曲線

    從圖4可以看出, 網(wǎng)格數(shù)量小于3.0×106時(shí)質(zhì)量流量變化較大, 當(dāng)網(wǎng)格數(shù)量增加到3.0×106后,繼續(xù)增加網(wǎng)格數(shù)量質(zhì)量流量基本保持不變,因此可以認(rèn)為網(wǎng)格數(shù)量已經(jīng)達(dá)到了網(wǎng)格無關(guān)性,故網(wǎng)格數(shù)量選擇為3.0×106。 在后續(xù)的計(jì)算過程中,隨著開度的變化, 閥芯與閥座的間隙也在變化,網(wǎng)格數(shù)量也應(yīng)做適當(dāng)?shù)恼{(diào)整,保持網(wǎng)格數(shù)量與上述網(wǎng)格數(shù)量相差不大,這樣既能加快計(jì)算的速度又能保證計(jì)算的準(zhǔn)確度。

    2.3 邊界條件設(shè)置

    調(diào)節(jié)閥工作介質(zhì)為水, 采用標(biāo)準(zhǔn)k-ε湍流模型,用SIMPLE算法和具有二階精度的迎風(fēng)格式進(jìn)行流場(chǎng)計(jì)算求解,邊界條件為壓力100 kPa,出口壓力0,壁面設(shè)置為光滑無滑移,殘差曲線精度設(shè)置為10-4;對(duì)入口流量和出口流量進(jìn)行檢測(cè),當(dāng)殘差曲線基本保持不變且所檢測(cè)的物理量保持不變時(shí)可認(rèn)為模擬收斂。

    3 結(jié)果分析

    對(duì)頂導(dǎo)向套筒單座調(diào)節(jié)閥20%、40%、60%、70%、80%、100%相對(duì)開度下的閥內(nèi)部流體對(duì)稱面上的壓力、速度及流線等進(jìn)行流動(dòng)規(guī)律分析。

    3.1 對(duì)稱面壓力云圖

    圖5是頂導(dǎo)向套筒單座調(diào)節(jié)閥不同開度下的對(duì)稱面壓力云圖。

    圖5 頂導(dǎo)向套筒單座調(diào)節(jié)閥不同開度下的壓力云圖

    從圖5可以看出: 閥門相鄰開度的壓力分布具有相似性, 當(dāng)調(diào)節(jié)閥處于20%~40%開度時(shí),流體在閥門前部的延伸流場(chǎng)內(nèi)的壓降很小,調(diào)節(jié)閥的整個(gè)壓降主要集中在節(jié)流孔部位,在閥芯靠近出口位置出現(xiàn)負(fù)壓區(qū);當(dāng)達(dá)到70%開度時(shí),節(jié)流孔處發(fā)生比較明顯的壓力梯度; 當(dāng)達(dá)到80%開度以后,不止節(jié)流孔周圍,節(jié)流孔下游處也有明顯的壓力變化。 說明調(diào)節(jié)閥在啟閉過程中閥芯和閥座之間的節(jié)流是產(chǎn)生壓降的重要因素之一。

    3.2 對(duì)稱面速度云圖

    圖6是頂導(dǎo)向套筒單座調(diào)節(jié)閥不同開度下的對(duì)稱面速度云圖。

    圖6 頂導(dǎo)向套筒單座調(diào)節(jié)閥不同開度下的速度云圖

    由圖6可知:當(dāng)調(diào)節(jié)閥處于20%~60%開度時(shí),進(jìn)出口速度基本均勻,閥芯與閥座之間的節(jié)流區(qū)域會(huì)出現(xiàn)明顯的速度梯度,這是由于流體在通過節(jié)流孔部位時(shí)流動(dòng)面積突然減小, 出現(xiàn)射流,在節(jié)流口后出現(xiàn)高速區(qū);當(dāng)開度達(dá)到70%以后,高速區(qū)出現(xiàn)在閥芯與套筒間,最大流速發(fā)生在閥芯出口側(cè),隨著開度呈先增大后減小的趨勢(shì);80%開度時(shí),流速達(dá)到最大(16.85 m/s)。

    閥體內(nèi)流體的最大速度出現(xiàn)在節(jié)流孔附近,高速流體對(duì)節(jié)流孔周圍沖擊會(huì)形成明顯的侵蝕區(qū), 沖刷區(qū)附近的流速高于附近區(qū)域的流速,高速度和低速度在閥芯和套筒之間的區(qū)域混合,從而產(chǎn)生強(qiáng)烈的剪切作用, 導(dǎo)致低速區(qū)域形成渦流, 并在對(duì)沖區(qū)域附近形成一個(gè)流速較低的區(qū)域,即低速空穴區(qū)。

    3.3 對(duì)稱面流線圖

    圖7是頂導(dǎo)向套筒單座調(diào)節(jié)閥不同開度下的對(duì)稱面流線圖。

    圖7 頂導(dǎo)向套筒單座調(diào)節(jié)閥不同開度下的流線圖

    從圖7可以看出, 在閥體中腔的右下側(cè)以及閥芯與套筒間、套筒與閥體間均有不同程度的漩渦產(chǎn)生,開度越小閥體中腔右下側(cè)產(chǎn)生漩渦的區(qū)域越大, 此區(qū)域只有部分流體參與主流流動(dòng),其余流體相互碰撞造成能量損失。 流體流經(jīng)節(jié)流孔后,一部分流體形成漩渦在閥芯與套筒間繞流后從套筒窗口流出, 另一部分直接從套筒窗口流出,在套筒與閥體間產(chǎn)生漩渦。 在閥門出口位置由于過流斷面面積的增大, 流體流動(dòng)不穩(wěn)定,出現(xiàn)二次流, 隨著流道的延長(zhǎng)流體流動(dòng)逐漸平穩(wěn),但此時(shí)流體呈螺旋態(tài)。

    對(duì)比壓力云圖、 速度云圖和流線圖分析可知:壓力云圖中的低壓區(qū)、速度云圖中的高速區(qū)和流線圖中密集區(qū)的位置相互對(duì)應(yīng),漩渦出現(xiàn)在高速區(qū)末端,表明流場(chǎng)中出現(xiàn)的漩渦區(qū)域?qū)毫退俣鹊淖兓泻艽笥绊懀瑫?huì)產(chǎn)生較大的能量損失,流阻系數(shù)增大;流體通過套筒后,流動(dòng)相對(duì)穩(wěn)定,壓力和速度分布也比較均勻,在管道出口內(nèi)出現(xiàn)二次流, 呈螺旋態(tài)向下游流動(dòng)且逐漸消失;流體流經(jīng)節(jié)流孔時(shí),能量損失主要是由流動(dòng)收縮產(chǎn)生的。

    3.4 流量系數(shù)和流阻系數(shù)

    計(jì)算流量系數(shù)和流阻系數(shù)的結(jié)果見表2。 由表2可知,在壓差近似不變時(shí),頂導(dǎo)向套筒單座調(diào)節(jié)閥的流量、流速和流量系數(shù)隨著相對(duì)開度的增大而增大,流阻系數(shù)隨著開度的增大而減小。

    表2 頂導(dǎo)向套筒單座調(diào)節(jié)閥流量系數(shù)和流阻系數(shù)

    4 閥門流量流阻試驗(yàn)

    4.1 試驗(yàn)設(shè)備

    試驗(yàn)是在閥門流量流阻測(cè)試裝置中完成的,該試驗(yàn)裝置測(cè)量范圍廣,適用于各類通用閥門的流量系數(shù)、流阻系數(shù)和流量特性測(cè)試。 試驗(yàn)采用變頻泵并聯(lián)和穩(wěn)壓罐供水技術(shù), 其穩(wěn)壓效果好。裝置配有高精度傳感器和先進(jìn)的數(shù)據(jù)采集器,可以同時(shí)采集流量、壓力、壓差及溫度等數(shù)據(jù)。 其測(cè)試軟件功能強(qiáng)大,操作簡(jiǎn)單,能夠提高閥門流阻測(cè)試的效率和自動(dòng)化水平。

    4.2 試驗(yàn)與模擬結(jié)果對(duì)比

    按 照GB/T 17213.9—2005/IEC 60534-2-3:1997《工業(yè)過程控制閥 第2-3部分:流通能力 試驗(yàn)程序》的要求安裝試驗(yàn)閥和完成試驗(yàn)程序。 在閥前與閥后2D和6D的位置(D為被測(cè)閥門的公稱直徑)設(shè)置取壓孔,測(cè)試介質(zhì)為水[12,13]。

    流量系數(shù)、流阻系數(shù)的模擬值與試驗(yàn)值之間的相對(duì)偏差計(jì)算式如下:

    將模擬所得的流量系數(shù)、流阻系數(shù)與試驗(yàn)值進(jìn)行對(duì)比,結(jié)果如圖8、9所示。

    由圖8、9可得, 數(shù)值模擬與試驗(yàn)所得的流量系數(shù)在不同開度下的變化趨勢(shì)基本相近,最大偏差值為-4.94%,小于標(biāo)準(zhǔn)要求的±5%;數(shù)值模擬與試驗(yàn)所得的流阻系數(shù)在不同開度下的變化趨勢(shì)基本相近,最大偏差值為-5.18%。 說明數(shù)值模擬得到的流量系數(shù)與流阻系數(shù)符合調(diào)節(jié)閥的實(shí)際調(diào)節(jié)特性,同時(shí)也驗(yàn)證了數(shù)值模擬計(jì)算方法的可靠性。

    圖9 頂導(dǎo)向套筒單座調(diào)節(jié)閥流阻系數(shù)試驗(yàn)值與模擬值對(duì)比

    5 結(jié)論

    5.1 通過閥門壓降分析看出,流體在閥門前部的延伸流場(chǎng)內(nèi)的壓降很小,調(diào)節(jié)閥的整個(gè)壓降主要集中在節(jié)流孔部位,節(jié)流孔處發(fā)生比較明顯的壓力梯度。 由此可得,調(diào)節(jié)閥在啟閉過程中閥芯和閥座之間的節(jié)流是產(chǎn)生壓降的重要因素之一。

    5.2 由不同開度下速度分析結(jié)果可以看出,進(jìn)出口速度基本均勻,閥芯與閥座之間的節(jié)流區(qū)域出現(xiàn)明顯的速度梯度,閥體內(nèi)流體的最大速度出現(xiàn)在節(jié)流孔附近,高速流動(dòng)對(duì)節(jié)流孔周圍沖擊形成明顯的侵蝕區(qū),高速流和低速流在閥芯和套筒之間混合,從而產(chǎn)生強(qiáng)烈的剪切作用,導(dǎo)致低速區(qū)域形成渦流,在對(duì)沖區(qū)域附近形成低速空穴區(qū)。

    5.3 在閥體中腔的右下側(cè)以及閥芯與套筒間、套筒與閥體間均有不同程度的漩渦產(chǎn)生,開度越小閥體中腔右下側(cè)產(chǎn)生漩渦的區(qū)域越大。 在閥門出口位置出現(xiàn)二次流,隨著流道的延長(zhǎng)流體流動(dòng)逐漸平穩(wěn),但此時(shí)流體呈螺旋態(tài)流動(dòng)。

    5.4 對(duì)比壓力云圖、速度云圖和流線圖可知,低壓區(qū)、高速區(qū)和流線密集區(qū)的位置相互對(duì)應(yīng),漩渦出現(xiàn)在高速區(qū)末端, 表明流場(chǎng)中出現(xiàn)的漩渦區(qū)域?qū)毫退俣鹊淖兓泻艽笥绊懀?導(dǎo)致能量損失增大,流阻系數(shù)增大,能量損失主要由流動(dòng)收縮產(chǎn)生。

    5.5 將流量系數(shù)和流阻系數(shù)的數(shù)值模擬值與試驗(yàn)值分別進(jìn)行對(duì)比,其相對(duì)偏差滿足標(biāo)準(zhǔn)要求的±5%,說明數(shù)值模擬計(jì)算方法可行、計(jì)算結(jié)果準(zhǔn)確。

    猜你喜歡
    節(jié)流調(diào)節(jié)閥套筒
    天然氣井井下節(jié)流器研究現(xiàn)狀及應(yīng)用前景
    重慶川儀調(diào)節(jié)閥有限公司
    重慶川儀調(diào)節(jié)閥有限公司
    套筒灌漿連接密實(shí)性檢測(cè)研究現(xiàn)狀及展望
    重慶川儀調(diào)節(jié)閥有限公司
    超高壓氣井井下節(jié)流技術(shù)應(yīng)用和設(shè)計(jì)方法
    PR方程模擬節(jié)流效應(yīng)的數(shù)值研究
    山東化工(2019年12期)2019-07-05 08:44:26
    一種調(diào)節(jié)閥口徑計(jì)算的方法
    一種尾架套筒自動(dòng)機(jī)械鎖緊機(jī)構(gòu)
    套筒類零件內(nèi)孔精加工工藝分析
    18禁在线播放成人免费| 黄色一级大片看看| 亚洲精品久久久久久婷婷小说| 26uuu在线亚洲综合色| 美女xxoo啪啪120秒动态图| 女人十人毛片免费观看3o分钟| 日本猛色少妇xxxxx猛交久久| 老司机影院毛片| 黄色配什么色好看| 99久久九九国产精品国产免费| 亚洲精华国产精华液的使用体验| 欧美变态另类bdsm刘玥| 18+在线观看网站| 久久鲁丝午夜福利片| 成人亚洲精品一区在线观看 | 久久99热这里只有精品18| 最近最新中文字幕免费大全7| 在线免费观看的www视频| 床上黄色一级片| 日本爱情动作片www.在线观看| 性色avwww在线观看| 蜜桃亚洲精品一区二区三区| 麻豆成人av视频| 国产黄色免费在线视频| 青青草视频在线视频观看| 国产精品国产三级国产专区5o| 午夜视频国产福利| 高清午夜精品一区二区三区| 亚洲成人av在线免费| 国产片特级美女逼逼视频| 国产人妻一区二区三区在| 99热6这里只有精品| 精品久久久久久久久久久久久| 国产毛片a区久久久久| 五月伊人婷婷丁香| 18禁在线无遮挡免费观看视频| 欧美成人精品欧美一级黄| 国产伦一二天堂av在线观看| 日韩中字成人| 一级毛片我不卡| 又爽又黄a免费视频| 网址你懂的国产日韩在线| 男人舔女人下体高潮全视频| 国语对白做爰xxxⅹ性视频网站| 精品人妻一区二区三区麻豆| 永久网站在线| 国产视频首页在线观看| 欧美3d第一页| av播播在线观看一区| 国产亚洲午夜精品一区二区久久 | 久久久久久久大尺度免费视频| 80岁老熟妇乱子伦牲交| freevideosex欧美| 啦啦啦中文免费视频观看日本| 如何舔出高潮| 亚洲色图av天堂| 又大又黄又爽视频免费| 一区二区三区乱码不卡18| 日韩欧美国产在线观看| 寂寞人妻少妇视频99o| 精品午夜福利在线看| 99久久精品国产国产毛片| 在线观看人妻少妇| 精品熟女少妇av免费看| 天堂√8在线中文| 好男人在线观看高清免费视频| 亚洲在线观看片| 午夜福利网站1000一区二区三区| a级毛片免费高清观看在线播放| 久久精品夜色国产| 国产探花极品一区二区| 菩萨蛮人人尽说江南好唐韦庄| 六月丁香七月| 黄色欧美视频在线观看| 汤姆久久久久久久影院中文字幕 | 欧美一级a爱片免费观看看| av.在线天堂| 97超视频在线观看视频| 亚洲国产色片| 熟妇人妻久久中文字幕3abv| 久久这里只有精品中国| 成人无遮挡网站| 国产 一区 欧美 日韩| 欧美日韩一区二区视频在线观看视频在线 | 黄色日韩在线| 国产亚洲91精品色在线| 麻豆精品久久久久久蜜桃| 搞女人的毛片| 国产爱豆传媒在线观看| 精品久久久久久久人妻蜜臀av| 人妻系列 视频| 非洲黑人性xxxx精品又粗又长| 熟妇人妻不卡中文字幕| 国产又色又爽无遮挡免| 最近中文字幕2019免费版| 日韩精品青青久久久久久| 3wmmmm亚洲av在线观看| 亚洲国产欧美在线一区| 两个人视频免费观看高清| 一本一本综合久久| 久久草成人影院| 午夜精品一区二区三区免费看| 99久久九九国产精品国产免费| 人人妻人人澡欧美一区二区| 91精品伊人久久大香线蕉| 国产久久久一区二区三区| 美女脱内裤让男人舔精品视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 18禁在线无遮挡免费观看视频| 国产大屁股一区二区在线视频| 亚洲第一区二区三区不卡| 一级片'在线观看视频| 人妻少妇偷人精品九色| 人妻制服诱惑在线中文字幕| 国产精品一及| 国产片特级美女逼逼视频| 午夜免费男女啪啪视频观看| 国产三级在线视频| 午夜福利在线观看吧| 黄色欧美视频在线观看| 免费看a级黄色片| 天美传媒精品一区二区| 精品亚洲乱码少妇综合久久| 高清午夜精品一区二区三区| 草草在线视频免费看| 一级毛片 在线播放| 国产伦一二天堂av在线观看| 少妇熟女aⅴ在线视频| 国内少妇人妻偷人精品xxx网站| 噜噜噜噜噜久久久久久91| 婷婷色av中文字幕| 精品一区二区三区人妻视频| 国产av国产精品国产| 欧美成人午夜免费资源| 国产精品久久久久久久电影| 舔av片在线| 久久草成人影院| 老司机影院毛片| 一边亲一边摸免费视频| 国产一级毛片在线| 国产精品久久久久久久电影| 日本免费a在线| 欧美高清性xxxxhd video| 欧美区成人在线视频| 国产成人精品婷婷| 国产单亲对白刺激| 两个人的视频大全免费| 一夜夜www| 伊人久久国产一区二区| 22中文网久久字幕| av又黄又爽大尺度在线免费看| 久久精品久久久久久噜噜老黄| 亚洲美女视频黄频| 一夜夜www| 久久精品久久久久久久性| 国产精品人妻久久久久久| 成人毛片a级毛片在线播放| 国产免费一级a男人的天堂| 亚洲综合色惰| 久久99热这里只频精品6学生| 成人二区视频| 波多野结衣巨乳人妻| 国产精品综合久久久久久久免费| 黑人高潮一二区| 日日干狠狠操夜夜爽| 在线 av 中文字幕| 亚洲美女搞黄在线观看| 国产又色又爽无遮挡免| 男的添女的下面高潮视频| 久久久久久久久久黄片| 97热精品久久久久久| 狠狠精品人妻久久久久久综合| 色5月婷婷丁香| 午夜日本视频在线| 国产亚洲午夜精品一区二区久久 | 一个人免费在线观看电影| 免费在线观看成人毛片| 免费观看精品视频网站| av在线蜜桃| 777米奇影视久久| 午夜精品在线福利| 欧美激情国产日韩精品一区| 国产 亚洲一区二区三区 | 亚洲精品日本国产第一区| 国产伦理片在线播放av一区| 久久人人爽人人爽人人片va| 特大巨黑吊av在线直播| 国内少妇人妻偷人精品xxx网站| 免费观看性生交大片5| 丝袜喷水一区| 婷婷色麻豆天堂久久| 色播亚洲综合网| 国产免费又黄又爽又色| 中文欧美无线码| 国产成人a∨麻豆精品| 可以在线观看毛片的网站| 国产极品天堂在线| 禁无遮挡网站| 欧美激情久久久久久爽电影| 久久精品综合一区二区三区| 九九爱精品视频在线观看| 久久久久久九九精品二区国产| 一区二区三区四区激情视频| 天堂av国产一区二区熟女人妻| 男女下面进入的视频免费午夜| 亚洲av免费高清在线观看| kizo精华| 美女黄网站色视频| 99热这里只有精品一区| 黄片无遮挡物在线观看| 亚洲国产色片| 人体艺术视频欧美日本| 日本爱情动作片www.在线观看| 国产日韩欧美在线精品| 国产不卡一卡二| 简卡轻食公司| 天堂中文最新版在线下载 | 精品久久久精品久久久| 日本熟妇午夜| 男女边吃奶边做爰视频| 亚洲人成网站在线播| 日本免费在线观看一区| 日韩大片免费观看网站| 亚洲国产精品成人久久小说| 80岁老熟妇乱子伦牲交| 免费观看a级毛片全部| 18禁在线无遮挡免费观看视频| 久久久a久久爽久久v久久| 国产老妇女一区| 国产国拍精品亚洲av在线观看| 国产成人aa在线观看| 中国美白少妇内射xxxbb| 又大又黄又爽视频免费| 干丝袜人妻中文字幕| 欧美xxⅹ黑人| 日日摸夜夜添夜夜爱| 国产淫语在线视频| 国产精品国产三级国产av玫瑰| 啦啦啦啦在线视频资源| 天天躁夜夜躁狠狠久久av| 六月丁香七月| 九草在线视频观看| 超碰av人人做人人爽久久| 亚洲在线观看片| 国产不卡一卡二| 国产中年淑女户外野战色| 国产男女超爽视频在线观看| 噜噜噜噜噜久久久久久91| 免费黄频网站在线观看国产| 久久精品久久久久久噜噜老黄| 狂野欧美白嫩少妇大欣赏| 免费av不卡在线播放| 国产精品爽爽va在线观看网站| 久久久a久久爽久久v久久| 欧美日韩精品成人综合77777| 乱系列少妇在线播放| 三级男女做爰猛烈吃奶摸视频| a级毛片免费高清观看在线播放| videos熟女内射| 国产午夜精品久久久久久一区二区三区| 久久久色成人| 亚洲电影在线观看av| 女人十人毛片免费观看3o分钟| 赤兔流量卡办理| 亚洲精品乱久久久久久| 极品少妇高潮喷水抽搐| 亚洲国产精品sss在线观看| 男女视频在线观看网站免费| 中文在线观看免费www的网站| 人人妻人人澡人人爽人人夜夜 | 国产视频首页在线观看| 午夜激情福利司机影院| 久久久久久久久中文| 亚洲精品久久久久久婷婷小说| 青春草亚洲视频在线观看| 超碰97精品在线观看| 日韩,欧美,国产一区二区三区| 日韩中字成人| 18禁裸乳无遮挡免费网站照片| 午夜精品一区二区三区免费看| 成人毛片60女人毛片免费| 日韩欧美三级三区| 久99久视频精品免费| 97超碰精品成人国产| 日韩欧美精品v在线| 日日摸夜夜添夜夜爱| 又爽又黄a免费视频| 美女大奶头视频| 亚洲精品第二区| 国产免费一级a男人的天堂| 成人无遮挡网站| 内射极品少妇av片p| 床上黄色一级片| 真实男女啪啪啪动态图| 色吧在线观看| 老师上课跳d突然被开到最大视频| 在线观看一区二区三区| 嘟嘟电影网在线观看| 男女啪啪激烈高潮av片| 成人性生交大片免费视频hd| 成人毛片60女人毛片免费| 亚洲va在线va天堂va国产| 九草在线视频观看| 日日摸夜夜添夜夜爱| 欧美人与善性xxx| 男女国产视频网站| 久久久欧美国产精品| 少妇的逼好多水| av在线老鸭窝| 亚洲自偷自拍三级| 欧美一级a爱片免费观看看| 久久久久久久久大av| 狠狠精品人妻久久久久久综合| 欧美日韩在线观看h| 99热全是精品| 亚洲精品亚洲一区二区| 最近的中文字幕免费完整| 亚洲第一区二区三区不卡| av国产久精品久网站免费入址| 岛国毛片在线播放| 啦啦啦中文免费视频观看日本| 色5月婷婷丁香| 国产亚洲最大av| 久久这里有精品视频免费| 亚洲熟妇中文字幕五十中出| 联通29元200g的流量卡| 麻豆国产97在线/欧美| av在线老鸭窝| 五月玫瑰六月丁香| 大又大粗又爽又黄少妇毛片口| 国产精品国产三级国产av玫瑰| 亚洲精品一区蜜桃| 热99在线观看视频| 一级毛片 在线播放| 九九爱精品视频在线观看| eeuss影院久久| 男人狂女人下面高潮的视频| 欧美日韩在线观看h| 国产视频内射| 亚洲乱码一区二区免费版| 欧美xxxx性猛交bbbb| av在线蜜桃| 男人和女人高潮做爰伦理| 国产有黄有色有爽视频| 一本久久精品| 人体艺术视频欧美日本| 国产亚洲精品久久久com| 在线天堂最新版资源| 97热精品久久久久久| 夫妻性生交免费视频一级片| 亚洲婷婷狠狠爱综合网| 熟女电影av网| 在线 av 中文字幕| 精品人妻偷拍中文字幕| 99热全是精品| 亚洲成人精品中文字幕电影| 成人二区视频| 又大又黄又爽视频免费| av黄色大香蕉| 日韩大片免费观看网站| 成人综合一区亚洲| 久久久久久久午夜电影| 97超碰精品成人国产| 国产高清国产精品国产三级 | 亚洲av成人精品一区久久| 免费不卡的大黄色大毛片视频在线观看 | 久久精品人妻少妇| 亚洲精品成人av观看孕妇| 久久热精品热| 99热这里只有是精品在线观看| 国产亚洲5aaaaa淫片| 天堂俺去俺来也www色官网 | 亚洲精品日本国产第一区| 国产精品一区二区性色av| 五月天丁香电影| 国产精品一区二区性色av| 成年人午夜在线观看视频 | 一本久久精品| 我的老师免费观看完整版| 亚洲最大成人手机在线| 国产爱豆传媒在线观看| a级毛色黄片| av在线观看视频网站免费| 免费在线观看成人毛片| 美女大奶头视频| 亚洲av不卡在线观看| 男女下面进入的视频免费午夜| 免费黄网站久久成人精品| 免费播放大片免费观看视频在线观看| 亚洲精品乱码久久久v下载方式| 日韩不卡一区二区三区视频在线| 中文在线观看免费www的网站| 在线免费观看的www视频| 午夜精品一区二区三区免费看| 九九久久精品国产亚洲av麻豆| 毛片女人毛片| 插阴视频在线观看视频| 免费观看精品视频网站| 少妇丰满av| 亚洲人与动物交配视频| 欧美不卡视频在线免费观看| 美女内射精品一级片tv| 免费观看av网站的网址| 久久久久网色| 亚洲精品乱码久久久久久按摩| 亚洲成人久久爱视频| 日韩一本色道免费dvd| 一级黄片播放器| 亚洲精品,欧美精品| 超碰97精品在线观看| 精华霜和精华液先用哪个| 嘟嘟电影网在线观看| 日韩伦理黄色片| av线在线观看网站| 亚洲自偷自拍三级| 亚洲人成网站在线观看播放| av在线天堂中文字幕| 亚洲性久久影院| 天美传媒精品一区二区| 国产精品1区2区在线观看.| 18禁裸乳无遮挡免费网站照片| 日本-黄色视频高清免费观看| 禁无遮挡网站| 男女边摸边吃奶| 天堂俺去俺来也www色官网 | 免费观看的影片在线观看| 亚洲av日韩在线播放| 简卡轻食公司| 国产一区二区三区av在线| 国产91av在线免费观看| 热99在线观看视频| 欧美一级a爱片免费观看看| 色综合站精品国产| 亚洲最大成人中文| 日韩欧美精品免费久久| 成人毛片a级毛片在线播放| 国产乱来视频区| 亚洲精品日韩av片在线观看| 日韩欧美国产在线观看| 国产一区二区亚洲精品在线观看| 一个人看视频在线观看www免费| 欧美成人a在线观看| 天堂中文最新版在线下载 | 亚洲人成网站在线播| 高清日韩中文字幕在线| 天堂av国产一区二区熟女人妻| 久久综合国产亚洲精品| h日本视频在线播放| 亚洲最大成人手机在线| 日产精品乱码卡一卡2卡三| 久久99热6这里只有精品| 日产精品乱码卡一卡2卡三| 国产黄色视频一区二区在线观看| 日韩,欧美,国产一区二区三区| 亚洲国产高清在线一区二区三| 国产精品三级大全| 久久久久性生活片| 国产精品一区二区性色av| 18禁裸乳无遮挡免费网站照片| 国产白丝娇喘喷水9色精品| 青青草视频在线视频观看| 中文字幕av在线有码专区| 日韩大片免费观看网站| 99久久人妻综合| 最新中文字幕久久久久| 干丝袜人妻中文字幕| 婷婷色麻豆天堂久久| 人妻少妇偷人精品九色| 亚洲欧美精品自产自拍| 一区二区三区乱码不卡18| 亚洲经典国产精华液单| 超碰av人人做人人爽久久| 中国美白少妇内射xxxbb| 噜噜噜噜噜久久久久久91| 大又大粗又爽又黄少妇毛片口| 国产一区亚洲一区在线观看| 亚洲久久久久久中文字幕| 日韩不卡一区二区三区视频在线| 亚洲国产成人一精品久久久| 黄片wwwwww| 熟妇人妻不卡中文字幕| 哪个播放器可以免费观看大片| 啦啦啦啦在线视频资源| 大话2 男鬼变身卡| 国产成人免费观看mmmm| 亚洲天堂国产精品一区在线| 又大又黄又爽视频免费| 成人亚洲精品av一区二区| 国产一区二区三区av在线| 中文字幕亚洲精品专区| 国产 一区精品| 麻豆av噜噜一区二区三区| 中文欧美无线码| 久久久久久久大尺度免费视频| 久久这里只有精品中国| 色尼玛亚洲综合影院| 日韩不卡一区二区三区视频在线| 亚洲不卡免费看| 久久精品久久久久久噜噜老黄| 亚洲不卡免费看| 99久久精品热视频| 老师上课跳d突然被开到最大视频| 久久久久久久久久黄片| 成年女人在线观看亚洲视频 | 久久久久久久久久黄片| 国产成人a区在线观看| 床上黄色一级片| 国产亚洲精品av在线| 啦啦啦中文免费视频观看日本| 亚洲av二区三区四区| av在线播放精品| 国产午夜精品一二区理论片| 听说在线观看完整版免费高清| 不卡视频在线观看欧美| 18禁在线播放成人免费| 亚洲最大成人中文| 18禁在线无遮挡免费观看视频| 在线a可以看的网站| 欧美成人a在线观看| 精品熟女少妇av免费看| 国产麻豆成人av免费视频| 久久久亚洲精品成人影院| 秋霞在线观看毛片| 亚洲欧美日韩无卡精品| 亚洲国产色片| 自拍偷自拍亚洲精品老妇| 99热这里只有是精品在线观看| 哪个播放器可以免费观看大片| 亚洲成人一二三区av| 汤姆久久久久久久影院中文字幕 | 国产伦精品一区二区三区四那| 99久久精品热视频| 一级av片app| 亚洲欧美成人综合另类久久久| 国产大屁股一区二区在线视频| 又粗又硬又长又爽又黄的视频| 中文字幕av在线有码专区| 视频中文字幕在线观看| 伊人久久精品亚洲午夜| 国产黄色小视频在线观看| .国产精品久久| 亚洲国产色片| 自拍偷自拍亚洲精品老妇| 欧美高清性xxxxhd video| 国产成年人精品一区二区| 80岁老熟妇乱子伦牲交| 久久99热这里只频精品6学生| 国产白丝娇喘喷水9色精品| 美女内射精品一级片tv| av.在线天堂| 乱系列少妇在线播放| 街头女战士在线观看网站| 中文字幕亚洲精品专区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av成人精品一区久久| 欧美另类一区| 免费看不卡的av| 日本爱情动作片www.在线观看| 大陆偷拍与自拍| 精品人妻熟女av久视频| 日韩欧美精品v在线| 欧美高清性xxxxhd video| 国产精品一区二区三区四区久久| 欧美97在线视频| 偷拍熟女少妇极品色| 成人亚洲精品一区在线观看 | 少妇熟女aⅴ在线视频| 久久午夜福利片| 中文天堂在线官网| 国产视频内射| 日韩制服骚丝袜av| 国产成人精品一,二区| 久久人人爽人人爽人人片va| av在线亚洲专区| 一二三四中文在线观看免费高清| 国产精品一区二区三区四区免费观看| 美女内射精品一级片tv| 男的添女的下面高潮视频| 久热久热在线精品观看| 一级毛片 在线播放| 大又大粗又爽又黄少妇毛片口| 最近视频中文字幕2019在线8| 国内精品一区二区在线观看| 97热精品久久久久久| 日韩欧美 国产精品| 99热全是精品| 亚洲精品日韩av片在线观看| 久久久久久久国产电影| 99久久精品热视频| 国产成人精品一,二区| 91久久精品国产一区二区三区| 在线播放无遮挡| 天天躁日日操中文字幕| 亚洲av日韩在线播放| 久久人人爽人人片av| 亚洲欧洲国产日韩| 中文字幕久久专区| 一级二级三级毛片免费看| .国产精品久久| 国产成人aa在线观看| 看非洲黑人一级黄片| 国产在线一区二区三区精| 精品久久久久久久久av| 天天一区二区日本电影三级| 亚洲欧美清纯卡通| 午夜日本视频在线| 狠狠精品人妻久久久久久综合| 国产成人免费观看mmmm| 国产激情偷乱视频一区二区| 一个人免费在线观看电影| 全区人妻精品视频| 建设人人有责人人尽责人人享有的 | 黄片wwwwww| 高清在线视频一区二区三区| 男女国产视频网站|