• <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)孔精加工工藝分析
    久久久久久人人人人人| 日韩 亚洲 欧美在线| 中文字幕人妻熟女乱码| 日韩,欧美,国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 一区二区三区精品91| 日韩精品免费视频一区二区三区| 黄色a级毛片大全视频| 中文欧美无线码| 国产不卡av网站在线观看| 丰满少妇做爰视频| 国产成人免费观看mmmm| 啦啦啦啦在线视频资源| 女警被强在线播放| 国产女主播在线喷水免费视频网站| av欧美777| 久久久久国产精品人妻一区二区| 亚洲精品久久午夜乱码| 亚洲欧美精品自产自拍| avwww免费| 一二三四在线观看免费中文在| 国产精品秋霞免费鲁丝片| 大陆偷拍与自拍| 亚洲国产毛片av蜜桃av| 97人妻天天添夜夜摸| 亚洲欧洲日产国产| 1024视频免费在线观看| 十八禁网站网址无遮挡| 成人国产一区最新在线观看 | 免费在线观看视频国产中文字幕亚洲 | 女人被躁到高潮嗷嗷叫费观| 国产精品一区二区免费欧美 | av有码第一页| tube8黄色片| 丝袜在线中文字幕| 99国产精品99久久久久| 观看av在线不卡| 在线看a的网站| 人妻一区二区av| 久久久久久久国产电影| 免费在线观看视频国产中文字幕亚洲 | 国产av一区二区精品久久| 久久久久久久大尺度免费视频| 国产成人a∨麻豆精品| 国产成人免费观看mmmm| 国产成人精品久久久久久| 国产欧美亚洲国产| 亚洲国产看品久久| 热re99久久精品国产66热6| 一级毛片我不卡| 一本色道久久久久久精品综合| 国产精品 国内视频| www.av在线官网国产| 欧美日韩黄片免| 久久精品国产综合久久久| 一本—道久久a久久精品蜜桃钙片| 51午夜福利影视在线观看| 精品一品国产午夜福利视频| 性色av一级| 丝袜脚勾引网站| 在线观看免费高清a一片| 国产片内射在线| 操美女的视频在线观看| 男人操女人黄网站| 精品福利永久在线观看| 九草在线视频观看| 日韩精品免费视频一区二区三区| 蜜桃国产av成人99| 亚洲国产精品成人久久小说| 国产成人精品久久二区二区免费| svipshipincom国产片| 波多野结衣一区麻豆| 大型av网站在线播放| 久久精品人人爽人人爽视色| 精品视频人人做人人爽| 晚上一个人看的免费电影| 我要看黄色一级片免费的| 久久综合国产亚洲精品| 亚洲欧洲日产国产| 亚洲成国产人片在线观看| 午夜老司机福利片| 亚洲av成人精品一二三区| 美女大奶头黄色视频| 亚洲av日韩精品久久久久久密 | 一二三四社区在线视频社区8| 久久久欧美国产精品| 国产在线观看jvid| 亚洲国产毛片av蜜桃av| 老汉色av国产亚洲站长工具| 香蕉丝袜av| 久久精品人人爽人人爽视色| 1024香蕉在线观看| avwww免费| 免费黄频网站在线观看国产| 男女边摸边吃奶| 久久狼人影院| 自线自在国产av| 国产精品久久久久久精品古装| 欧美日韩av久久| 超碰成人久久| 午夜91福利影院| 国产免费一区二区三区四区乱码| 啦啦啦视频在线资源免费观看| 久久久久久久久免费视频了| 中文字幕最新亚洲高清| √禁漫天堂资源中文www| 午夜激情久久久久久久| netflix在线观看网站| 水蜜桃什么品种好| 亚洲欧美日韩另类电影网站| 热re99久久精品国产66热6| 18禁观看日本| 久久久久久久大尺度免费视频| 一级黄片播放器| 婷婷成人精品国产| 亚洲五月色婷婷综合| 一边摸一边做爽爽视频免费| 国产精品国产三级国产专区5o| 97精品久久久久久久久久精品| 下体分泌物呈黄色| 夜夜骑夜夜射夜夜干| 亚洲国产最新在线播放| 丁香六月欧美| 高潮久久久久久久久久久不卡| 午夜免费鲁丝| 咕卡用的链子| 久久国产精品人妻蜜桃| 青春草亚洲视频在线观看| 国产高清videossex| 国产亚洲午夜精品一区二区久久| av网站在线播放免费| 国产亚洲午夜精品一区二区久久| 2018国产大陆天天弄谢| 国产成人欧美在线观看 | 亚洲九九香蕉| 久久女婷五月综合色啪小说| 成年人免费黄色播放视频| 国产欧美日韩精品亚洲av| 精品久久久久久久毛片微露脸 | 国产精品一二三区在线看| 国产真人三级小视频在线观看| 青春草视频在线免费观看| 亚洲国产毛片av蜜桃av| 国产成人系列免费观看| 男女无遮挡免费网站观看| 少妇被粗大的猛进出69影院| 在线观看免费视频网站a站| 在线观看免费视频网站a站| avwww免费| kizo精华| 一二三四社区在线视频社区8| 在线精品无人区一区二区三| 高清视频免费观看一区二区| 亚洲,欧美精品.| 精品人妻熟女毛片av久久网站| 黄色视频不卡| 国产三级黄色录像| 精品久久久久久久毛片微露脸 | 黄色 视频免费看| 国产亚洲一区二区精品| 一边摸一边做爽爽视频免费| 久久久久久免费高清国产稀缺| 午夜两性在线视频| 国产成人精品久久二区二区免费| 精品少妇久久久久久888优播| 日本一区二区免费在线视频| 国产成人av教育| 久久久精品国产亚洲av高清涩受| 一区二区av电影网| 菩萨蛮人人尽说江南好唐韦庄| 成人国产av品久久久| 国产精品久久久久成人av| 99九九在线精品视频| 一区二区三区乱码不卡18| 一级毛片我不卡| 国产精品熟女久久久久浪| 亚洲国产精品成人久久小说| 一二三四在线观看免费中文在| 久久久精品区二区三区| 晚上一个人看的免费电影| 嫩草影视91久久| 涩涩av久久男人的天堂| 亚洲欧美色中文字幕在线| 久久狼人影院| 两个人免费观看高清视频| 久久av网站| 大香蕉久久网| 久久青草综合色| 免费一级毛片在线播放高清视频 | 亚洲 欧美一区二区三区| 亚洲国产精品成人久久小说| 美女视频免费永久观看网站| 国产精品香港三级国产av潘金莲 | 久久国产精品人妻蜜桃| 国产成人免费无遮挡视频| 曰老女人黄片| 老熟女久久久| 国产日韩欧美在线精品| 在线精品无人区一区二区三| 满18在线观看网站| 日本午夜av视频| 国产男人的电影天堂91| av一本久久久久| 啦啦啦在线免费观看视频4| 久久精品国产综合久久久| 欧美精品啪啪一区二区三区 | e午夜精品久久久久久久| www.熟女人妻精品国产| 国产成人a∨麻豆精品| 国产精品免费视频内射| 国产成人系列免费观看| 国产成人一区二区在线| 久久精品国产亚洲av高清一级| 69精品国产乱码久久久| 校园人妻丝袜中文字幕| 2021少妇久久久久久久久久久| 久久人人爽人人片av| 亚洲图色成人| 涩涩av久久男人的天堂| 亚洲男人天堂网一区| 欧美日韩亚洲综合一区二区三区_| 国产女主播在线喷水免费视频网站| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻人人澡人人看| 高清av免费在线| 成人18禁高潮啪啪吃奶动态图| 精品卡一卡二卡四卡免费| 水蜜桃什么品种好| www.自偷自拍.com| 两个人免费观看高清视频| 亚洲中文av在线| 亚洲av欧美aⅴ国产| 又大又黄又爽视频免费| 一边亲一边摸免费视频| 国产1区2区3区精品| 中文乱码字字幕精品一区二区三区| 嫩草影视91久久| av视频免费观看在线观看| 成人亚洲欧美一区二区av| 亚洲五月婷婷丁香| 亚洲综合色网址| 亚洲国产最新在线播放| 亚洲精品一卡2卡三卡4卡5卡 | 成年女人毛片免费观看观看9 | 国产成人精品无人区| 国产欧美日韩综合在线一区二区| 一级片'在线观看视频| 高潮久久久久久久久久久不卡| 色播在线永久视频| 少妇 在线观看| 国产野战对白在线观看| 成年人免费黄色播放视频| 1024视频免费在线观看| 国产精品三级大全| 国产日韩欧美视频二区| 亚洲成人国产一区在线观看 | 国精品久久久久久国模美| 老司机影院成人| 三上悠亚av全集在线观看| 亚洲精品一二三| 久久免费观看电影| 国产精品亚洲av一区麻豆| 精品一区在线观看国产| 99久久99久久久精品蜜桃| 午夜两性在线视频| 亚洲成人免费av在线播放| 欧美人与善性xxx| 亚洲欧美一区二区三区黑人| 亚洲人成电影观看| 99热网站在线观看| 国产99久久九九免费精品| 亚洲午夜精品一区,二区,三区| 国产在线一区二区三区精| 一本一本久久a久久精品综合妖精| 久久99一区二区三区| 国产成人一区二区三区免费视频网站 | 宅男免费午夜| 亚洲一码二码三码区别大吗| 国产精品人妻久久久影院| 午夜91福利影院| 午夜视频精品福利| 黄色一级大片看看| 18禁黄网站禁片午夜丰满| 青青草视频在线视频观看| 丰满饥渴人妻一区二区三| 国产精品香港三级国产av潘金莲 | 欧美久久黑人一区二区| av不卡在线播放| 久久久精品免费免费高清| 免费观看a级毛片全部| 国产淫语在线视频| 日韩,欧美,国产一区二区三区| 日韩制服骚丝袜av| 美女福利国产在线| 久久人人爽人人片av| 日韩电影二区| 美女大奶头黄色视频| 丁香六月欧美| 午夜福利乱码中文字幕| 日韩,欧美,国产一区二区三区| 老司机午夜十八禁免费视频| 老鸭窝网址在线观看| 黄色片一级片一级黄色片| 人人妻人人添人人爽欧美一区卜| 亚洲欧美日韩另类电影网站| 新久久久久国产一级毛片| 日韩,欧美,国产一区二区三区| 日韩大码丰满熟妇| 啦啦啦在线观看免费高清www| 精品久久久精品久久久| 男人舔女人的私密视频| 亚洲成av片中文字幕在线观看| 中国美女看黄片| 亚洲精品国产一区二区精华液| 夜夜骑夜夜射夜夜干| 好男人电影高清在线观看| 人人妻人人爽人人添夜夜欢视频| 老鸭窝网址在线观看| 18禁裸乳无遮挡动漫免费视频| 真人做人爱边吃奶动态| 在线观看www视频免费| 在线观看免费视频网站a站| 99精国产麻豆久久婷婷| 观看av在线不卡| 亚洲情色 制服丝袜| 亚洲av综合色区一区| 少妇精品久久久久久久| 国产亚洲欧美在线一区二区| 免费高清在线观看视频在线观看| 久久久国产一区二区| 高清欧美精品videossex| 人人妻,人人澡人人爽秒播 | 国产精品av久久久久免费| av片东京热男人的天堂| 欧美日韩视频高清一区二区三区二| 一二三四在线观看免费中文在| 国产精品久久久av美女十八| 成在线人永久免费视频| 婷婷色综合大香蕉| 久久综合国产亚洲精品| 午夜激情av网站| 美女午夜性视频免费| 香蕉丝袜av| 色婷婷av一区二区三区视频| 日日爽夜夜爽网站| 蜜桃国产av成人99| 欧美人与善性xxx| 国产黄频视频在线观看| 搡老乐熟女国产| 一级黄片播放器| 国产精品香港三级国产av潘金莲 | 色视频在线一区二区三区| 国产亚洲av高清不卡| 成人午夜精彩视频在线观看| 交换朋友夫妻互换小说| 夫妻午夜视频| 天堂俺去俺来也www色官网| 国产91精品成人一区二区三区 | 曰老女人黄片| 国产av一区二区精品久久| 亚洲国产欧美一区二区综合| 大片电影免费在线观看免费| 丰满少妇做爰视频| 精品亚洲成国产av| 成人亚洲欧美一区二区av| 国产极品粉嫩免费观看在线| 久久精品久久精品一区二区三区| 晚上一个人看的免费电影| 久久精品国产亚洲av涩爱| 一区在线观看完整版| 中文字幕av电影在线播放| 成年av动漫网址| 亚洲av电影在线进入| e午夜精品久久久久久久| 高清欧美精品videossex| 日本欧美视频一区| 欧美在线一区亚洲| 国产视频首页在线观看| 婷婷色av中文字幕| 岛国毛片在线播放| 亚洲国产av新网站| 欧美乱码精品一区二区三区| 一个人免费看片子| 高清黄色对白视频在线免费看| www.精华液| 国产视频一区二区在线看| 蜜桃在线观看..| 国产色视频综合| 91精品国产国语对白视频| 亚洲伊人色综图| 少妇裸体淫交视频免费看高清 | 汤姆久久久久久久影院中文字幕| 在线观看国产h片| 午夜免费成人在线视频| 黄频高清免费视频| 欧美精品av麻豆av| 国产高清视频在线播放一区 | av有码第一页| 国产精品人妻久久久影院| 久久精品亚洲熟妇少妇任你| 国产真人三级小视频在线观看| 亚洲国产欧美一区二区综合| 午夜免费鲁丝| 日本五十路高清| 色网站视频免费| av网站在线播放免费| 男人爽女人下面视频在线观看| 日韩精品免费视频一区二区三区| 美女午夜性视频免费| 午夜福利免费观看在线| 99久久精品国产亚洲精品| 99国产精品一区二区三区| 亚洲成人国产一区在线观看 | 精品福利观看| 人人妻人人爽人人添夜夜欢视频| 观看av在线不卡| 王馨瑶露胸无遮挡在线观看| 国产一区有黄有色的免费视频| 国产女主播在线喷水免费视频网站| 国产精品欧美亚洲77777| 国产日韩欧美视频二区| 美女高潮到喷水免费观看| 国产一区亚洲一区在线观看| 欧美日韩精品网址| 亚洲av美国av| 日韩一卡2卡3卡4卡2021年| 欧美精品啪啪一区二区三区 | 亚洲自偷自拍图片 自拍| 亚洲国产看品久久| 电影成人av| 色视频在线一区二区三区| 老熟女久久久| 满18在线观看网站| 精品第一国产精品| 嫩草影视91久久| 欧美精品人与动牲交sv欧美| 久久久久久久精品精品| 新久久久久国产一级毛片| 久久免费观看电影| 精品卡一卡二卡四卡免费| 午夜影院在线不卡| 国产精品成人在线| 热re99久久国产66热| 精品视频人人做人人爽| 别揉我奶头~嗯~啊~动态视频 | 亚洲中文日韩欧美视频| 美女中出高潮动态图| 后天国语完整版免费观看| 中文字幕精品免费在线观看视频| 久久天堂一区二区三区四区| 国产一区二区在线观看av| 9热在线视频观看99| 九草在线视频观看| 久久久国产一区二区| 国产精品久久久人人做人人爽| 激情视频va一区二区三区| 精品国产一区二区三区久久久樱花| 丰满人妻熟妇乱又伦精品不卡| 国产一级毛片在线| 久久精品国产亚洲av高清一级| 黄片播放在线免费| 国产又爽黄色视频| 婷婷色综合大香蕉| 欧美日韩av久久| 欧美激情极品国产一区二区三区| 一边亲一边摸免费视频| 狠狠精品人妻久久久久久综合| 一级片免费观看大全| 大片免费播放器 马上看| kizo精华| 欧美国产精品一级二级三级| 99久久人妻综合| 如日韩欧美国产精品一区二区三区| 777久久人妻少妇嫩草av网站| 中文字幕人妻丝袜制服| 久久青草综合色| 一区福利在线观看| 99热国产这里只有精品6| 97人妻天天添夜夜摸| 欧美激情高清一区二区三区| 久久精品久久精品一区二区三区| av在线老鸭窝| 亚洲色图 男人天堂 中文字幕| 99国产精品一区二区三区| 中文字幕人妻丝袜一区二区| 日韩伦理黄色片| 国产精品一区二区精品视频观看| 少妇裸体淫交视频免费看高清 | 亚洲五月婷婷丁香| 国语对白做爰xxxⅹ性视频网站| 亚洲国产日韩一区二区| 精品一区在线观看国产| 新久久久久国产一级毛片| 欧美成狂野欧美在线观看| 在现免费观看毛片| 大香蕉久久网| www.熟女人妻精品国产| 女人高潮潮喷娇喘18禁视频| 日本av免费视频播放| 日韩 亚洲 欧美在线| 亚洲国产精品一区二区三区在线| 日本欧美国产在线视频| 国产免费现黄频在线看| 国产人伦9x9x在线观看| 久久亚洲国产成人精品v| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久大尺度免费视频| 男女下面插进去视频免费观看| 午夜福利视频精品| 国产一卡二卡三卡精品| 一边亲一边摸免费视频| 国产成人系列免费观看| 色婷婷av一区二区三区视频| 婷婷色av中文字幕| 亚洲激情五月婷婷啪啪| 欧美精品啪啪一区二区三区 | 亚洲精品成人av观看孕妇| 美女脱内裤让男人舔精品视频| 精品少妇久久久久久888优播| 国产淫语在线视频| 操美女的视频在线观看| 激情视频va一区二区三区| 亚洲一码二码三码区别大吗| 99久久综合免费| 亚洲中文av在线| av国产久精品久网站免费入址| 你懂的网址亚洲精品在线观看| 丰满少妇做爰视频| 日本a在线网址| 狂野欧美激情性bbbbbb| 丝袜人妻中文字幕| 老汉色∧v一级毛片| 欧美日韩福利视频一区二区| 老司机在亚洲福利影院| 叶爱在线成人免费视频播放| xxxhd国产人妻xxx| 久久毛片免费看一区二区三区| 亚洲国产欧美在线一区| 亚洲自偷自拍图片 自拍| 看免费成人av毛片| 国产日韩一区二区三区精品不卡| 国产淫语在线视频| 日韩大片免费观看网站| 三上悠亚av全集在线观看| 一区二区三区精品91| 9热在线视频观看99| 亚洲九九香蕉| 99九九在线精品视频| 黄色视频不卡| 丰满迷人的少妇在线观看| 精品人妻在线不人妻| 日韩中文字幕视频在线看片| 亚洲国产精品国产精品| 中文字幕人妻熟女乱码| 亚洲三区欧美一区| 一级,二级,三级黄色视频| 亚洲自偷自拍图片 自拍| 99国产综合亚洲精品| 国产成人免费无遮挡视频| 国产精品国产av在线观看| 国产男女内射视频| 在线观看一区二区三区激情| 亚洲av男天堂| 老熟女久久久| 在线观看免费午夜福利视频| 免费在线观看影片大全网站 | 国产精品人妻久久久影院| 人人妻,人人澡人人爽秒播 | 欧美久久黑人一区二区| 国产午夜精品一二区理论片| 日韩av免费高清视频| 国产又爽黄色视频| 中文字幕色久视频| 无遮挡黄片免费观看| 十分钟在线观看高清视频www| 日韩制服丝袜自拍偷拍| 亚洲人成电影免费在线| 每晚都被弄得嗷嗷叫到高潮| 欧美精品一区二区免费开放| 亚洲熟女毛片儿| 最新的欧美精品一区二区| 精品少妇一区二区三区视频日本电影| 久久99一区二区三区| 我要看黄色一级片免费的| 亚洲av日韩精品久久久久久密 | 亚洲精品成人av观看孕妇| 久久天躁狠狠躁夜夜2o2o | 免费观看人在逋| 午夜福利影视在线免费观看| 国产av精品麻豆| 高清欧美精品videossex| 狠狠婷婷综合久久久久久88av| 亚洲欧美一区二区三区久久| 波野结衣二区三区在线| 精品少妇内射三级| 精品国产一区二区三区久久久樱花| 亚洲欧美清纯卡通| 大香蕉久久网| 人妻人人澡人人爽人人| 欧美xxⅹ黑人| 日本欧美国产在线视频| 亚洲激情五月婷婷啪啪| 在线观看www视频免费| 中文欧美无线码| 亚洲中文字幕日韩| 日韩av在线免费看完整版不卡| 欧美成人午夜精品| 久久九九热精品免费| 亚洲中文字幕日韩| 亚洲人成电影免费在线| 99久久人妻综合| 日韩中文字幕视频在线看片| 欧美亚洲日本最大视频资源| 国产成人av教育| 男女午夜视频在线观看|