• <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)孔精加工工藝分析
    亚洲男人天堂网一区| 国产又色又爽无遮挡免| 欧美日韩国产mv在线观看视频| 国产精品三级大全| 一区二区三区四区激情视频| 韩国av在线不卡| 女的被弄到高潮叫床怎么办| 免费在线观看完整版高清| 成人亚洲精品一区在线观看| 日韩伦理黄色片| 国产片内射在线| 国产成人a∨麻豆精品| 亚洲欧美激情在线| 亚洲精品久久成人aⅴ小说| 午夜福利视频精品| 欧美成人午夜精品| 久久天堂一区二区三区四区| 亚洲av国产av综合av卡| 欧美精品av麻豆av| 国产精品一二三区在线看| 男女无遮挡免费网站观看| 成人亚洲欧美一区二区av| 午夜日本视频在线| 一本色道久久久久久精品综合| 国产免费视频播放在线视频| 亚洲专区中文字幕在线 | 国产 精品1| 日本猛色少妇xxxxx猛交久久| 亚洲精品久久午夜乱码| 亚洲国产欧美日韩在线播放| 久久久欧美国产精品| 亚洲av日韩在线播放| 亚洲av男天堂| 美国免费a级毛片| 国产精品一二三区在线看| 考比视频在线观看| 日韩大码丰满熟妇| 精品第一国产精品| 热re99久久精品国产66热6| 日韩电影二区| 老鸭窝网址在线观看| 成年女人毛片免费观看观看9 | 2021少妇久久久久久久久久久| 黑人巨大精品欧美一区二区蜜桃| 久久女婷五月综合色啪小说| 亚洲免费av在线视频| 热re99久久国产66热| 如何舔出高潮| 五月开心婷婷网| 成人黄色视频免费在线看| 久久久久人妻精品一区果冻| 久久久久久免费高清国产稀缺| 赤兔流量卡办理| 久久久亚洲精品成人影院| 亚洲综合色网址| 婷婷色综合www| av线在线观看网站| 99国产综合亚洲精品| 精品久久久久久电影网| 婷婷色综合www| 在线观看人妻少妇| 在线看a的网站| 日韩中文字幕视频在线看片| 久久久久精品性色| 国产成人午夜福利电影在线观看| 大香蕉久久成人网| 精品一区二区三区av网在线观看 | 亚洲欧美一区二区三区黑人| 国产激情久久老熟女| 五月天丁香电影| 在线观看一区二区三区激情| 欧美精品人与动牲交sv欧美| 国产人伦9x9x在线观看| 蜜桃国产av成人99| www.av在线官网国产| 午夜av观看不卡| 国产精品女同一区二区软件| 午夜精品国产一区二区电影| 叶爱在线成人免费视频播放| 午夜激情久久久久久久| 一区二区三区四区激情视频| 婷婷色av中文字幕| 18禁国产床啪视频网站| 国产熟女欧美一区二区| 日韩av免费高清视频| 曰老女人黄片| 亚洲第一青青草原| 欧美日韩福利视频一区二区| av国产久精品久网站免费入址| 国产亚洲精品第一综合不卡| 午夜福利网站1000一区二区三区| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美激情在线| 久久久久久久久久久免费av| 国产男人的电影天堂91| 国产亚洲av高清不卡| 欧美老熟妇乱子伦牲交| 中文字幕人妻熟女乱码| 一本久久精品| 考比视频在线观看| 在线观看免费午夜福利视频| 亚洲四区av| 国产色婷婷99| 欧美亚洲 丝袜 人妻 在线| 久久热在线av| 天天躁夜夜躁狠狠躁躁| 国产伦人伦偷精品视频| 久久久久精品国产欧美久久久 | 欧美日韩视频高清一区二区三区二| 在线观看www视频免费| 亚洲成av片中文字幕在线观看| 一区在线观看完整版| 高清不卡的av网站| 一级爰片在线观看| 国产精品国产三级专区第一集| 我的亚洲天堂| 女人高潮潮喷娇喘18禁视频| 男女边摸边吃奶| 国产高清国产精品国产三级| 成年av动漫网址| 99精国产麻豆久久婷婷| 国产亚洲午夜精品一区二区久久| 蜜桃国产av成人99| 欧美日韩亚洲综合一区二区三区_| 在线观看三级黄色| 热re99久久精品国产66热6| 国产欧美亚洲国产| 成年女人毛片免费观看观看9 | 国产欧美日韩一区二区三区在线| 肉色欧美久久久久久久蜜桃| 国产亚洲最大av| 亚洲精品中文字幕在线视频| 香蕉丝袜av| 亚洲伊人色综图| 久久久精品区二区三区| 麻豆乱淫一区二区| 男人舔女人的私密视频| 男人添女人高潮全过程视频| 日韩,欧美,国产一区二区三区| 国产黄色免费在线视频| 国产精品人妻久久久影院| 黑人欧美特级aaaaaa片| 一级毛片我不卡| 亚洲熟女毛片儿| 波多野结衣av一区二区av| 久久鲁丝午夜福利片| 久久久亚洲精品成人影院| 日韩一卡2卡3卡4卡2021年| 男男h啪啪无遮挡| 精品亚洲成国产av| 这个男人来自地球电影免费观看 | 亚洲成人av在线免费| 欧美日韩亚洲国产一区二区在线观看 | 国产精品成人在线| 国语对白做爰xxxⅹ性视频网站| 亚洲第一av免费看| 亚洲色图 男人天堂 中文字幕| 99久久99久久久精品蜜桃| 欧美精品高潮呻吟av久久| 伊人亚洲综合成人网| 精品第一国产精品| 国产一区有黄有色的免费视频| 国产精品 国内视频| 国产一区二区 视频在线| 晚上一个人看的免费电影| 国产男女超爽视频在线观看| 中国国产av一级| 丝瓜视频免费看黄片| 黄频高清免费视频| 人人妻人人添人人爽欧美一区卜| 在线观看一区二区三区激情| av福利片在线| 欧美精品亚洲一区二区| 国产精品.久久久| 久久久久视频综合| 老司机靠b影院| 最近中文字幕2019免费版| 一级片'在线观看视频| 热re99久久精品国产66热6| 成人亚洲欧美一区二区av| √禁漫天堂资源中文www| 亚洲国产欧美日韩在线播放| 老鸭窝网址在线观看| 久久99热这里只频精品6学生| 91成人精品电影| 视频在线观看一区二区三区| 婷婷色综合www| 男女无遮挡免费网站观看| 制服诱惑二区| 国产精品久久久久久精品电影小说| 亚洲精品美女久久久久99蜜臀 | 美国免费a级毛片| 赤兔流量卡办理| 视频区图区小说| 亚洲欧美中文字幕日韩二区| 成人亚洲精品一区在线观看| 国产一区二区 视频在线| 女人高潮潮喷娇喘18禁视频| 亚洲熟女精品中文字幕| 一级a爱视频在线免费观看| 19禁男女啪啪无遮挡网站| 爱豆传媒免费全集在线观看| 亚洲专区中文字幕在线 | 老汉色∧v一级毛片| 欧美日韩亚洲综合一区二区三区_| 国产精品 欧美亚洲| 久久久久久人人人人人| 亚洲国产欧美一区二区综合| 亚洲综合精品二区| 亚洲精品美女久久av网站| 国产成人欧美在线观看 | 国产一区二区三区综合在线观看| 日韩av免费高清视频| 久久久久久久大尺度免费视频| 叶爱在线成人免费视频播放| 亚洲,一卡二卡三卡| 国产免费又黄又爽又色| 日韩不卡一区二区三区视频在线| 日韩中文字幕欧美一区二区 | 狂野欧美激情性xxxx| 丰满饥渴人妻一区二区三| 色网站视频免费| av一本久久久久| 极品少妇高潮喷水抽搐| 欧美人与性动交α欧美精品济南到| 亚洲精品国产av成人精品| 免费在线观看完整版高清| 一本—道久久a久久精品蜜桃钙片| 久久久久久免费高清国产稀缺| 免费观看av网站的网址| 精品少妇内射三级| 久久精品国产亚洲av涩爱| 日韩中文字幕视频在线看片| 久久精品国产a三级三级三级| 午夜日本视频在线| 国产精品嫩草影院av在线观看| 日韩,欧美,国产一区二区三区| 日韩av免费高清视频| 伊人久久大香线蕉亚洲五| 一级片免费观看大全| 色播在线永久视频| 十八禁网站网址无遮挡| 制服人妻中文乱码| 青春草视频在线免费观看| 视频区图区小说| 精品少妇一区二区三区视频日本电影 | 制服诱惑二区| 在线免费观看不下载黄p国产| 欧美最新免费一区二区三区| 久久久久久久久免费视频了| 久热这里只有精品99| 亚洲欧美成人综合另类久久久| 97精品久久久久久久久久精品| 免费高清在线观看日韩| 国产99久久九九免费精品| 久热这里只有精品99| 日日撸夜夜添| 久久精品国产a三级三级三级| 一区二区三区精品91| 亚洲自偷自拍图片 自拍| 不卡视频在线观看欧美| 亚洲精品国产色婷婷电影| videos熟女内射| 欧美日韩视频精品一区| 麻豆精品久久久久久蜜桃| 亚洲七黄色美女视频| 卡戴珊不雅视频在线播放| 亚洲成人国产一区在线观看 | 国产 精品1| 大陆偷拍与自拍| 亚洲精品久久午夜乱码| 午夜福利,免费看| 美女扒开内裤让男人捅视频| 国产成人系列免费观看| 中文字幕人妻熟女乱码| 亚洲精品国产区一区二| 91国产中文字幕| 亚洲精品日本国产第一区| 国产深夜福利视频在线观看| 少妇被粗大的猛进出69影院| 欧美最新免费一区二区三区| 久久这里只有精品19| 久久久久久久精品精品| 国产精品蜜桃在线观看| 亚洲成av片中文字幕在线观看| 免费观看a级毛片全部| av国产精品久久久久影院| 欧美成人精品欧美一级黄| 日韩熟女老妇一区二区性免费视频| 在线观看免费午夜福利视频| netflix在线观看网站| tube8黄色片| 亚洲欧美一区二区三区黑人| 多毛熟女@视频| 国产精品一二三区在线看| 国产一区二区激情短视频 | 久久毛片免费看一区二区三区| 男人操女人黄网站| 女性被躁到高潮视频| 19禁男女啪啪无遮挡网站| 久久久久视频综合| 麻豆精品久久久久久蜜桃| 成人午夜精彩视频在线观看| 成人手机av| 在线观看国产h片| 91精品伊人久久大香线蕉| 久久精品国产a三级三级三级| 日韩av在线免费看完整版不卡| 国产成人a∨麻豆精品| 亚洲伊人色综图| 视频在线观看一区二区三区| 日韩欧美精品免费久久| 国产毛片在线视频| 视频区图区小说| 国产精品免费视频内射| 色网站视频免费| xxxhd国产人妻xxx| 青春草视频在线免费观看| 国产一区二区在线观看av| 久久久精品区二区三区| 纵有疾风起免费观看全集完整版| 少妇猛男粗大的猛烈进出视频| 国产精品久久久人人做人人爽| 一本—道久久a久久精品蜜桃钙片| 老司机靠b影院| 人成视频在线观看免费观看| 久久影院123| 亚洲成人免费av在线播放| 午夜av观看不卡| 最近2019中文字幕mv第一页| 人体艺术视频欧美日本| 日本av手机在线免费观看| h视频一区二区三区| 日日撸夜夜添| 在线亚洲精品国产二区图片欧美| 亚洲av国产av综合av卡| 国产精品嫩草影院av在线观看| 性高湖久久久久久久久免费观看| 伦理电影免费视频| 久久国产精品男人的天堂亚洲| 日韩精品有码人妻一区| 亚洲精品成人av观看孕妇| 在线看a的网站| 一级毛片我不卡| 十八禁网站网址无遮挡| 成人三级做爰电影| 国产爽快片一区二区三区| 国产精品久久久久久人妻精品电影 | 天美传媒精品一区二区| 国产精品麻豆人妻色哟哟久久| 青草久久国产| 久久精品国产综合久久久| 1024香蕉在线观看| 国产精品人妻久久久影院| 天天躁狠狠躁夜夜躁狠狠躁| 可以免费在线观看a视频的电影网站 | 午夜av观看不卡| 一区二区三区激情视频| 黄片无遮挡物在线观看| 亚洲av综合色区一区| 汤姆久久久久久久影院中文字幕| 另类精品久久| 国产熟女午夜一区二区三区| 中文字幕制服av| 波多野结衣一区麻豆| 美女视频免费永久观看网站| 日韩精品有码人妻一区| av在线app专区| 桃花免费在线播放| 高清在线视频一区二区三区| 可以免费在线观看a视频的电影网站 | av福利片在线| 久久精品国产a三级三级三级| 亚洲国产av新网站| 丝袜喷水一区| 成人三级做爰电影| 国产午夜精品一二区理论片| 青春草亚洲视频在线观看| 国产伦理片在线播放av一区| av在线播放精品| 大香蕉久久网| 国产免费一区二区三区四区乱码| 黄频高清免费视频| 9热在线视频观看99| 欧美在线一区亚洲| 99热全是精品| 国产熟女午夜一区二区三区| 国产深夜福利视频在线观看| 国产精品亚洲av一区麻豆 | 久久人人爽av亚洲精品天堂| 不卡视频在线观看欧美| 午夜日韩欧美国产| 天堂中文最新版在线下载| 五月天丁香电影| 男的添女的下面高潮视频| 美女国产高潮福利片在线看| 一级,二级,三级黄色视频| 狂野欧美激情性bbbbbb| 香蕉国产在线看| 亚洲精品中文字幕在线视频| 亚洲av成人精品一二三区| 亚洲av在线观看美女高潮| 免费观看a级毛片全部| 亚洲专区中文字幕在线 | av视频免费观看在线观看| av天堂久久9| 99re6热这里在线精品视频| 亚洲 欧美一区二区三区| 性高湖久久久久久久久免费观看| 国产在视频线精品| 天天操日日干夜夜撸| 少妇人妻精品综合一区二区| 国产一区二区三区av在线| 老汉色∧v一级毛片| 久久精品久久久久久噜噜老黄| 人人妻,人人澡人人爽秒播 | 黄片播放在线免费| 成人免费观看视频高清| 美国免费a级毛片| 国产免费福利视频在线观看| 九色亚洲精品在线播放| 一区二区三区精品91| 欧美成人精品欧美一级黄| 亚洲欧洲国产日韩| 亚洲第一av免费看| 国产高清国产精品国产三级| 黄色视频不卡| 观看美女的网站| 99久久人妻综合| 少妇人妻精品综合一区二区| 人妻人人澡人人爽人人| 日韩av免费高清视频| 人人妻,人人澡人人爽秒播 | 97在线人人人人妻| 99精品久久久久人妻精品| 国产精品国产三级国产专区5o| 中文天堂在线官网| 超碰成人久久| 国产精品.久久久| 人人澡人人妻人| 考比视频在线观看| 色婷婷av一区二区三区视频| 最近的中文字幕免费完整| 香蕉丝袜av| av在线播放精品| 免费看不卡的av| 国产97色在线日韩免费| 三上悠亚av全集在线观看| 亚洲精品视频女| 男人爽女人下面视频在线观看| 中文字幕人妻丝袜制服| 久久人人爽av亚洲精品天堂| 国产又爽黄色视频| 午夜日本视频在线| 青春草国产在线视频| 久久精品亚洲av国产电影网| 国产 精品1| 婷婷色麻豆天堂久久| 在线免费观看不下载黄p国产| 亚洲精品自拍成人| 亚洲欧美精品综合一区二区三区| 人成视频在线观看免费观看| 777久久人妻少妇嫩草av网站| 高清黄色对白视频在线免费看| 不卡av一区二区三区| 肉色欧美久久久久久久蜜桃| 欧美人与性动交α欧美软件| 51午夜福利影视在线观看| 超碰成人久久| 满18在线观看网站| 丝袜脚勾引网站| 亚洲激情五月婷婷啪啪| 麻豆乱淫一区二区| 黄片无遮挡物在线观看| 国产精品久久久久久久久免| 中文字幕人妻熟女乱码| 九九爱精品视频在线观看| 亚洲美女视频黄频| 国产精品一二三区在线看| 电影成人av| 国产一级毛片在线| 天天操日日干夜夜撸| 亚洲国产精品成人久久小说| 女性生殖器流出的白浆| 两个人看的免费小视频| 日韩精品免费视频一区二区三区| 国产精品秋霞免费鲁丝片| 天堂俺去俺来也www色官网| 高清黄色对白视频在线免费看| 嫩草影院入口| 99re6热这里在线精品视频| 男女边吃奶边做爰视频| 男女午夜视频在线观看| 一区在线观看完整版| 久久人人爽人人片av| 免费观看人在逋| 中文字幕高清在线视频| 99国产精品免费福利视频| 欧美另类一区| 久久国产精品男人的天堂亚洲| 日韩一卡2卡3卡4卡2021年| 亚洲国产欧美日韩在线播放| 亚洲人成电影观看| 久久久久久久久免费视频了| 久久久国产欧美日韩av| 国产极品粉嫩免费观看在线| 日韩制服丝袜自拍偷拍| 成人手机av| 一级黄片播放器| 亚洲精品自拍成人| 青春草国产在线视频| 精品久久久久久电影网| 久久久欧美国产精品| av卡一久久| 国产精品 国内视频| 中文字幕亚洲精品专区| 国产男女内射视频| 黄色一级大片看看| 久久久久精品性色| 日韩精品有码人妻一区| 成人午夜精彩视频在线观看| 丝袜美足系列| 日韩av不卡免费在线播放| 亚洲国产成人一精品久久久| 亚洲综合色网址| 欧美日本中文国产一区发布| 中文乱码字字幕精品一区二区三区| 黄片无遮挡物在线观看| 一级爰片在线观看| 亚洲精品久久午夜乱码| 亚洲一码二码三码区别大吗| 中文字幕亚洲精品专区| 国产精品免费大片| videosex国产| a 毛片基地| 国产国语露脸激情在线看| 秋霞伦理黄片| 丝袜美足系列| 婷婷色麻豆天堂久久| 国产乱人偷精品视频| 国产亚洲最大av| 欧美日韩综合久久久久久| 老司机在亚洲福利影院| 日韩伦理黄色片| 午夜免费观看性视频| 在线观看免费高清a一片| 中文字幕另类日韩欧美亚洲嫩草| 中文乱码字字幕精品一区二区三区| 母亲3免费完整高清在线观看| 中文字幕制服av| 国产精品三级大全| 欧美激情 高清一区二区三区| 最近中文字幕高清免费大全6| 夫妻性生交免费视频一级片| 亚洲欧美激情在线| 亚洲成国产人片在线观看| 99国产精品免费福利视频| 老司机亚洲免费影院| 免费高清在线观看视频在线观看| 波多野结衣av一区二区av| 黄色怎么调成土黄色| 51午夜福利影视在线观看| 飞空精品影院首页| 国产精品免费大片| 国产精品一区二区在线观看99| 精品少妇一区二区三区视频日本电影 | 免费av中文字幕在线| 夫妻性生交免费视频一级片| 国产精品一二三区在线看| 精品久久蜜臀av无| 国产在线免费精品| 少妇人妻精品综合一区二区| 成人国产av品久久久| 激情五月婷婷亚洲| 亚洲成色77777| 亚洲欧美激情在线| 精品久久蜜臀av无| 女人久久www免费人成看片| 少妇 在线观看| 成人国语在线视频| 国产在线视频一区二区| 美女国产高潮福利片在线看| 一区二区av电影网| kizo精华| 久久精品亚洲熟妇少妇任你| 日本一区二区免费在线视频| kizo精华| 色婷婷久久久亚洲欧美| 日韩不卡一区二区三区视频在线| 亚洲成人国产一区在线观看 | 国产精品香港三级国产av潘金莲 | 丁香六月天网| 精品少妇久久久久久888优播| 国产成人精品久久二区二区91 | 国产精品偷伦视频观看了| 国产又爽黄色视频| 热re99久久国产66热| 青草久久国产| avwww免费| 好男人视频免费观看在线| 一区二区三区精品91| 久久狼人影院| 欧美日韩国产mv在线观看视频| 亚洲色图 男人天堂 中文字幕| 日韩一区二区视频免费看| 一区二区三区乱码不卡18| 国产 一区精品| 日韩一区二区视频免费看| 人妻人人澡人人爽人人| 99香蕉大伊视频| 1024香蕉在线观看| 久久久久网色| 日韩一卡2卡3卡4卡2021年| 精品少妇黑人巨大在线播放| 日韩 亚洲 欧美在线|