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

    液柱分離狀態(tài)下水錘二次升壓現(xiàn)象計算與研究

    2022-05-09 05:38:08解朝蓬
    三峽大學學報(自然科學版) 2022年2期
    關鍵詞:液柱水錘波速

    富 友 解朝蓬 蔣 勁

    (1.蘭州理工大學 能源與動力工程學院,蘭州 730050;2.安徽三聯(lián)泵業(yè)股份有限公司,安徽 馬鞍山 243000;3.武漢大學 動力與機械學院,武漢 430072)

    在液柱分離水錘現(xiàn)象中,二次升壓現(xiàn)象一直是研究的重點內(nèi)容.這種壓力激升往往伴隨著空泡的生成與潰滅,在閥門處尤為明顯,其壓力激升甚至比根據(jù)Joukowsky公式計算的最大壓力還要高得多[1].有學者分析了二次升壓過程中閥門處空泡的初生—發(fā)育—潰滅的過程,從而探究空泡潰滅對流動特性的影響[2];也有試驗結果表明,二次升壓現(xiàn)象在低瞬變強度下會尤為突出[3],但都對其內(nèi)部壓力波的變化規(guī)律研究較少.

    在伴有液柱分離的瞬變流動分析中,空泡的生成使得水錘壓力波的變化程度與單相水錘不同.探索更加精確的瞬變流動摩阻模型對水錘壓力波進行捕捉,有助于分析二次升壓壓力波的內(nèi)部變化規(guī)律[4].此外,這種二次升壓的產(chǎn)生會對機組、閥門和儀表產(chǎn)生較大的危害,所以對其進行研究是非常有必要的.本文基于雙流體數(shù)學模型,采用數(shù)值模擬與試驗對比相結合的研究方法,構建出一種改良的水錘數(shù)值解法,目的是為了探討二次升壓的內(nèi)部流動機理,分析邊界水頭、邊界壓力、水錘波速、管道長度和管道直徑對二次升壓的影響.

    1 計算方法

    1.1 雙流體數(shù)學模型

    本文采用歐拉-歐拉雙流體模型[5-7],該模型中兩相均為可壓縮.模型包含一個空泡輸運方程、兩個連續(xù)性方程、兩個動量方程和兩個能量方程.表達式如下:

    方程中:下標1代表氣相;下標2代表液相;α為空穴率;u為速度(m/s);ρ為密度(kg/m3);p為壓力(Pa);J u為非恒定摩阻;J s=fu|u|/2D;D為管道直徑(m);J g為重力項,J g=gsinθ;μ為壓力弛豫系數(shù);λ為速度弛豫系數(shù);E=e+u2/2;pI為兩相相間壓力,pI=α1p1+α2p2;uI為兩相相間速度,uI=(α1ρ1u1+α2ρ2u2)/(α1ρ1+α2ρ2);α1+α2=1.

    為使上述方程組封閉,引入狀態(tài)方程:

    其中:γk為比熱容比,γk=C p,k/C v,k;C p,k為定壓比熱容;C v,k為定容比熱容.π1通常為0,π2表達式見式(9):

    其中:Tref為試驗環(huán)境溫度;psat為飽和蒸汽壓.

    波速方程為:

    方程(8~10)中:k=1代表氣相;k=2代表液相.

    1.1.1 相間弛豫系數(shù)

    在使用雙流體模型時,必須考慮瞬態(tài)流動中的壓力和速度弛豫問題.其決定了模型的計算穩(wěn)定性[8-9].相間弛豫系數(shù)包括速度弛豫系數(shù)和壓力弛豫系數(shù).

    速度弛豫系數(shù)可表達為:

    在壓力弛豫系數(shù)中引入空泡靜態(tài)平衡條件:

    其中:S為表面張力;R為泡粒半徑.

    進而改進的壓力弛豫系數(shù)表示為:

    1.1.2 非恒定摩阻模型

    瞬態(tài)流動模型可采用單系數(shù)和雙系數(shù)摩阻模型來計算壓力波的衰減.單系數(shù)和雙系數(shù)摩阻模型由于時間和空間偏微分項的存在,使其在雙流體模型中的應用變得復雜.加權函數(shù)非恒定摩阻模型[10]是對Zielke模型和Trikha模型[11]的簡化.基于該模型的數(shù)值編碼,使用加權函數(shù),可以節(jié)省大量計算所需的時間.

    對Trikha模型簡化,通過加權平均,將瞬態(tài)流動阻力模型轉化為如下形式:

    式中:

    其中:υ為運動粘度(m2/s);g為重力加速度(m/s2);y i為速度變量(m/s);y i中i取1、2、3,對應m i和n i中3個數(shù).

    1.2 求解雙流體模型的數(shù)值方法

    1.2.1 Godunov-Rusanov格式

    將方程組(1~7)進行重新整理,忽略空泡輸運方程中的弛豫源項.整理后如下:

    其中,U=(αρ,αρu,αρE)T,F=(αρ,αρu2+αP,αρE+αP)T,H=(0,pI,pIuI).

    礦床內(nèi)發(fā)育了兩組斷裂,一組為NNW向張扭性斷裂,該組斷裂以扭性為主,因扭力是順時針方向的,故多為右斜列式出現(xiàn),同時還伴生有小角度的扭裂隙和張裂隙。在這些構造裂隙帶中充填有閃鋅礦、輝銻礦、脆硫銻鉛礦和石英、含錳菱鐵礦、含鐵菱錳礦、鐵白云石等金屬礦物和脈石。該組斷裂多成群平行排列出現(xiàn),地表出露30余條。它是本礦床的主要控礦構造,礦脈的產(chǎn)狀及形態(tài),嚴格受其控制。第二組為NNE向扭性斷裂,分布在礦床西側,多充填有乳白色石英脈,局部有輝銻礦化,有切穿NNW向斷裂的現(xiàn)象,為第二期形成的斷裂(圖2)。在礦床主要地段未發(fā)現(xiàn),對礦脈無破壞性。礦區(qū)內(nèi)主要控礦構造及性質見表1。

    求解格式如下:

    式中:Δ=?α1/?x,其離散格式為

    式中:S=max{|u i±c i|,|u i+1-c i+1|,|u i-1-c i-1|}.

    因為空泡輸運方程(16)為非守恒型方程,不能采用上述守恒形式格式求解.因此,對其采用如下求解格式:

    1.2.2 邊界條件

    如圖1所示,邊界條件的構造采用特征線法(圖中橫坐標為網(wǎng)格長度上的網(wǎng)格節(jié)點,縱坐標為時間層).

    圖1 特征值及特征線在空間坐標的分布趨勢d x/d t=u±c

    相容性方程可以表示為:

    式中:φR代表當特征波速為負值時,向上游傳遞的信息節(jié)點;φL代表當特征波速為正值時,向下游傳遞的信息節(jié)點;φA,φB,φC是不同網(wǎng)格節(jié)點上的值.

    2 數(shù)值計算結果與試驗結果對比

    管線長度分為5個部分,管道長度在0~0.65 m之間,管道標高為0.47 m;管道長度在0.65~24.83 m 之間,標高為0.26 m;管道長度在24.83~49.17 m之間,標高為0.13 m;管道長度在在73.52~98.56 m之間,標高為0 m.計算時間步長為5×10-5s,管道波速為1 240~1 291.8 m/s.在管線沿線有4處壓力檢測點,分別位于管線1/4、2/4、3/4 和4/4 處.文獻[12]給出了相應的流體熱力學參數(shù).使得該模型較好地保證了計算的穩(wěn)定性和精度.

    試驗系統(tǒng)原理圖如圖2所示.

    圖2 試驗系統(tǒng)示意圖

    表1 試驗具體參數(shù)

    數(shù)值模擬結果與試驗結果對比如圖3~8所示.

    圖3 管線P3處流速為2.74 m/s的水頭隨時間的變化特征

    圖4 管線P1處流速為2.74 m/s的水頭隨時間的變化特征

    圖5 管線P3處流速為1.43 m/s的水頭隨時間的變化特征

    圖6 管線P1處流速為1.43 m/s的水頭隨時間的變化特征

    圖7 管線P3處流速為0.47 m/s的水頭隨時間的變化特征

    無論是管道末端還是管道中點的模擬結果都表明:采用非恒定摩阻的新模型,不僅可以捕捉到液柱分離過程中水擊的最大升壓和降壓值,而且還可以捕捉最大升壓之后的壓力衰減變化.

    3 內(nèi)部流動特性分析

    3.1 水錘二次升壓分析

    從圖8可以看出,管道末端的水錘壓力波的二次壓力升高值明顯高于首級波峰.這種現(xiàn)象是由液柱分離再彌合產(chǎn)生的壓力波反射疊加所導致.

    圖8 管線P1處流速為0.47m/s的水頭隨時間的變化特征

    為探究上述二次升壓現(xiàn)象的流動機理,以“閥門-管道-穩(wěn)壓水箱”系統(tǒng)為研究對象,忽略摩阻影響.選取兩種瞬變流動現(xiàn)象,第1種是初始流速為0.47 m/s,初始水頭為70 m 無二次升壓的水錘現(xiàn)象,其壓力分布如圖9(a)所示;第2種是初始流速為0.47 m/s,初始水頭為50 m 有二次升壓的水錘現(xiàn)象,其壓力分布如圖9(b)所示.

    圖9 閥門處兩種瞬態(tài)流動壓力分布

    圖10和圖11為壓力波在多相瞬變流中的傳播過程,其中實線為升壓波,虛線為降壓波.壓力波的傳遞過程遵循兩個規(guī)律:(1)壓力波與壁面、閥門等邊界條件相接觸時,其反射波與原波性質相同;(2)壓力波與穩(wěn)壓邊界相作用時,其反射波與原波性質相異.

    圖10 無二次升壓的水錘現(xiàn)象

    圖11 伴有二次升壓的水錘現(xiàn)象

    如圖10所示,壓力波1、4、5為升壓波,2、3、6為降壓波.在無液柱分離情況下,壓力變化按照理論計算公式(Hmax/min=H0±,H0為邊界水頭)下降與上升,壓力下降與上升的變化范圍如圖12所示(橫坐標0代表上游,1代表下游),該條件下并未出現(xiàn)二次升壓現(xiàn)象.

    圖12 管線壓力波分布

    如圖11所示,當末尾閥門快速關閉,閥門處產(chǎn)生升壓波1向上游傳播,升壓波在上游與穩(wěn)壓罐作用形成降壓波2反射至閥門處,傳播時間為2L/c.閥門處的低壓波2和3反射導致壓力降低至飽和蒸汽壓,閥門處發(fā)生空化.降壓波3從上游穩(wěn)壓罐處變化為升壓波4反射,升壓波4在管道末端與空泡相遇,此時波速產(chǎn)生突變,如圖13所示.

    圖13 波速隨管線變化

    升壓波4在兩相混合區(qū)域產(chǎn)生折射與反射兩種壓力波:第1種壓力波與空泡作用反射形成降壓波5向上游傳遞;第2種壓力波在兩相區(qū)域以升壓波4繼續(xù)向閥門邊界折射傳遞.當壓力波4傳遞至閥門處,反射為升壓波6向上游傳遞.在降壓波5和升壓波6的傳遞期間,相同管線位置處降壓波5 早于升壓波6,其變化如圖14(a)所示.

    升壓波6傳遞至上游處變化為降壓波8向下游傳遞,降壓波5傳遞至上游變化為升壓波7向下游傳遞.在相同管線位置升壓波7早于降壓波8,其變化如圖14(b)所示.升壓波7在上游處與升壓波6相遇,兩種相同性質的正壓波在N 點疊加傳遞至末端閥門處,最終使得末尾閥門處壓力峰值高于理論計算值.

    圖14 管線壓力波分布

    3.2 水錘二次升壓隨系統(tǒng)參數(shù)變化

    為了判定水錘二次升壓值隨系統(tǒng)參數(shù)的變化,本文選取了邊界水頭、邊界流速、管道直徑、管道長度、和水錘波速5個參數(shù)來進行內(nèi)部流動特性分析.

    3.2.1 邊界水頭對二次升壓的影響

    圖15統(tǒng)計了二次升壓值(二次峰值)在特定初始流速下隨邊界水頭的變化,圖中橫坐標為邊界水頭(H0)與Joukowsky升壓(ΔH=(cv0)/g)的比值,縱坐標表示二次升壓值(H2max)與Joukowsky 升壓值(ΔH)的比值.通過分析發(fā)現(xiàn),在特定的初始流速下,二次升壓值隨著初始水頭的增大呈現(xiàn)出先增大后減小的趨勢,并且在H0≈0.8ΔH處會取得極值.此結論應用在水錘泵的輸水過程中,可以使其獲得更高的揚程.

    圖15 二次升壓值隨邊界水頭的變化

    3.2.2 邊界流速對二次升壓的影響

    表2統(tǒng)計了二次升壓值在特定初始水頭下隨初始流速的變化.可以看出,在初始水頭為24.8 m 條件下,初始流速在0.30~0.82 m/s之間,二次升壓值隨著初始流速的增大而升高,當初始流速大于1.02 m/s,二次升壓現(xiàn)象便不再發(fā)生.在初始水頭為41.8 m的條件下,初始流速為0.30 m/s,并沒有發(fā)生二次升壓現(xiàn)象,因為在此條件下的初始邊界壓力大于Joukowsky升壓值,不能發(fā)生液柱分離,從而使得一次峰值與二次峰值相等.而初始流速在0.47~1.14 m/s之間,二次升壓值隨著初始流速的增大而增大,當初始流速大于1.35 m/s,二次升壓現(xiàn)象便不再發(fā)生.在初始水頭為71.8 m 條件下,初始流速在0.64~1.35 m/s之間,二次升壓值隨著初始流速的增大而增大,當初始流速大于1.70 m/s,二次升壓現(xiàn)象便不再發(fā)生.以上說明,在液柱分離過程中,二次升壓值隨著初始流速的增大而增大,當初始流速大于臨界流速,二次升壓便不會發(fā)生.因此,在工程應用中,如果要避免因液柱分離而產(chǎn)生的二次升壓,可以通過數(shù)值模擬來選擇最優(yōu)工況.

    表2 二次升壓值隨初始流速變化(單位:Pa)

    3.2.3 管道直徑、管道長度和水錘波速對二次升壓的影響

    表3統(tǒng)計了在初始邊界H0=0.8ΔH條件下,二次升壓相對值隨管道直徑、管道長度和管道波速的變化.在理論情況下,壓力波在水中的傳播速度不超過1 450 m/s,實際波速值要小于理論值,并且為使得該條件下發(fā)生液柱分離,選取波速最小值為1 100 m/s.

    表3 二次升壓相對值H2max/ΔH 的變化情況

    可以看出,隨著管道直徑的增加,二次升壓相對值保持不變;隨著管道長度的增加,二次升壓值相應增大;隨著水錘波速的增大,二次升壓值隨之增大.

    3.2.4 二次升壓試驗驗證

    為了驗證二次升壓值隨邊界水頭的變化規(guī)律,本節(jié)對圖15的模擬結果進行了相應的實驗驗證.

    表4 試驗具體數(shù)據(jù)

    首先進行了初始流速為0.348 3 m/s,邊界水頭為25.428 m 的水錘試驗.在此基礎上,保持初始流速基本不變,取邊界水頭為H0=0.6ΔH、H0=0.7ΔH、H0=0.8ΔH和H0=0.9ΔH,分別對相應試驗2、3、4、5進行驗證.試驗結果如圖16所示,可以看出,在H0=0.8ΔH處,二次升壓值最大.

    圖16 二次升壓值隨邊界水頭變化

    圖17 管線中點處水頭隨時間的變化特征

    上述試驗結果對圖15的模擬結果進行了驗證.從水錘壓力波隨管線分布可以看出,由于摩阻的存在,升壓波在閥門處反射向上游傳遞的過程中持續(xù)衰減,并且在管線中點處二次升壓現(xiàn)象不再發(fā)生.

    4 結論

    本文通過對瞬態(tài)流動中的二次升壓現(xiàn)象進行分析,得到了以下結論:

    1)在伴有液柱分離的水錘中,空泡的生成與潰滅會使得管道內(nèi)部發(fā)生壓力波的疊加.管道中兩種升壓波疊加沿管線傳遞至閥門處,使得閥門處流體進一步被壓縮,從而產(chǎn)生二次升壓現(xiàn)象,這種二次升壓值會高于理論計算值.

    2)二次升壓值在H0≈0.8ΔH處取得極值,隨初始水頭增大呈現(xiàn)先增大后減小的趨勢.二次升壓值隨著初始流速的增大而增大,當初始流速大于臨界流速,二次升壓便不會再發(fā)生.二次升壓值隨水錘波速和管道長度的增大而增大,管道直徑對二次升壓值無影響.

    3)通過5 組試驗結果驗證了二次升壓值在H0≈0.8ΔH處取得極值.在升壓波從閥門處向上游傳遞的過程中,管線阻力對其衰減影響顯著.

    猜你喜歡
    液柱水錘波速
    巧用“形象思維”,速解液柱(活塞)移動問題
    高水頭短距離泵站水錘計算分析
    基于實測波速探討地震反射波法超前預報解譯標志
    培養(yǎng)科學思維 落實核心素養(yǎng)
    水錘防護措施在某長距離供水系統(tǒng)中的應用
    輸水管道水錘防護方案對比分析
    吉林地區(qū)波速比分布特征及構造意義
    地震地質(2015年3期)2015-12-25 03:29:41
    豎直液柱與水平液面作用激起毛細波探究
    液柱在激波沖擊下RM不穩(wěn)定性和破裂過程的數(shù)值計算
    基于分位數(shù)回歸的剪切波速變化規(guī)律
    午夜视频国产福利| 亚洲美女视频黄频| 国产精品一区二区在线观看99 | 久久久久久伊人网av| 午夜福利高清视频| 久久精品国产亚洲av涩爱| 最近最新中文字幕大全电影3| 色视频www国产| 老司机影院成人| 嫩草影院入口| 天天躁日日操中文字幕| 日韩成人av中文字幕在线观看| 一级毛片电影观看 | 丰满少妇做爰视频| 深爱激情五月婷婷| 国国产精品蜜臀av免费| 国产高清视频在线观看网站| 变态另类丝袜制服| 久久99热这里只有精品18| 日韩欧美精品v在线| 神马国产精品三级电影在线观看| 日本免费在线观看一区| 亚洲,欧美,日韩| 免费观看精品视频网站| 十八禁国产超污无遮挡网站| 一卡2卡三卡四卡精品乱码亚洲| 桃色一区二区三区在线观看| 欧美xxxx性猛交bbbb| 午夜福利视频1000在线观看| 亚洲欧美日韩卡通动漫| av在线天堂中文字幕| 高清av免费在线| 国产成人精品婷婷| 日本wwww免费看| 在线观看一区二区三区| 中文欧美无线码| 国产爱豆传媒在线观看| 亚洲成色77777| 国产伦理片在线播放av一区| 一级毛片aaaaaa免费看小| 国产成人精品一,二区| 成年av动漫网址| 日本三级黄在线观看| 热99re8久久精品国产| 日本一二三区视频观看| 国产精品熟女久久久久浪| av国产久精品久网站免费入址| 免费观看人在逋| 在线观看美女被高潮喷水网站| 男女国产视频网站| 精品久久国产蜜桃| 久久久久久伊人网av| 免费av毛片视频| 男人舔女人下体高潮全视频| 亚洲精品亚洲一区二区| 国产爱豆传媒在线观看| 国产免费一级a男人的天堂| 国产高潮美女av| 国产精品.久久久| 国产亚洲精品av在线| 久久精品国产99精品国产亚洲性色| 国产在视频线精品| 搡老妇女老女人老熟妇| 搡老妇女老女人老熟妇| 一级二级三级毛片免费看| 亚洲自偷自拍三级| 亚州av有码| 狂野欧美激情性xxxx在线观看| 偷拍熟女少妇极品色| 一级二级三级毛片免费看| 男女国产视频网站| 亚洲精品国产成人久久av| 少妇的逼好多水| 有码 亚洲区| 噜噜噜噜噜久久久久久91| 色综合站精品国产| 长腿黑丝高跟| 黑人高潮一二区| 黄片wwwwww| 日韩中字成人| 国产 一区精品| 久久99热这里只有精品18| 久久人人爽人人爽人人片va| 亚洲综合色惰| 哪个播放器可以免费观看大片| 在线观看一区二区三区| 久久这里只有精品中国| 26uuu在线亚洲综合色| 精品久久久久久电影网 | av在线老鸭窝| 大话2 男鬼变身卡| 亚洲欧洲国产日韩| 神马国产精品三级电影在线观看| 国产午夜精品久久久久久一区二区三区| kizo精华| 国产色爽女视频免费观看| a级毛色黄片| 一区二区三区免费毛片| 久久精品久久精品一区二区三区| 亚洲综合色惰| 少妇丰满av| av免费观看日本| 久久精品91蜜桃| 99视频精品全部免费 在线| 伊人久久精品亚洲午夜| 中文字幕av在线有码专区| 国产一区有黄有色的免费视频 | 日韩一本色道免费dvd| 小说图片视频综合网站| 人人妻人人澡人人爽人人夜夜 | 午夜激情福利司机影院| 人妻少妇偷人精品九色| 亚洲成人久久爱视频| 亚洲av免费高清在线观看| 国产亚洲精品久久久com| 免费观看的影片在线观看| 能在线免费看毛片的网站| 桃色一区二区三区在线观看| 两个人视频免费观看高清| 青青草视频在线视频观看| 99九九线精品视频在线观看视频| 中文字幕熟女人妻在线| 亚洲av一区综合| 男女视频在线观看网站免费| 好男人视频免费观看在线| 日韩欧美精品免费久久| 最近手机中文字幕大全| 成人漫画全彩无遮挡| 深夜a级毛片| 成年版毛片免费区| 特级一级黄色大片| 欧美激情久久久久久爽电影| 九草在线视频观看| 最近中文字幕2019免费版| 寂寞人妻少妇视频99o| 天天躁夜夜躁狠狠久久av| 国产av在哪里看| 26uuu在线亚洲综合色| 99久久中文字幕三级久久日本| 精品久久久久久久久av| 国产精品精品国产色婷婷| 日日干狠狠操夜夜爽| 97热精品久久久久久| 中文乱码字字幕精品一区二区三区 | 老司机影院毛片| 亚洲成人av在线免费| 亚洲国产精品合色在线| 欧美三级亚洲精品| av在线观看视频网站免费| 亚洲欧美精品专区久久| 日产精品乱码卡一卡2卡三| av播播在线观看一区| av在线观看视频网站免费| 26uuu在线亚洲综合色| 精品久久久久久久久av| 一级二级三级毛片免费看| av在线播放精品| 99久久无色码亚洲精品果冻| 亚洲欧美日韩高清专用| 三级毛片av免费| 婷婷色av中文字幕| 夜夜爽夜夜爽视频| 国产成人午夜福利电影在线观看| 亚洲,欧美,日韩| 九草在线视频观看| 欧美+日韩+精品| 在线观看66精品国产| 1000部很黄的大片| 亚洲成人中文字幕在线播放| 亚洲av男天堂| 免费大片18禁| 国语对白做爰xxxⅹ性视频网站| av女优亚洲男人天堂| 三级国产精品片| 国产 一区 欧美 日韩| 国产一区有黄有色的免费视频 | 日本一本二区三区精品| 日韩大片免费观看网站 | 69av精品久久久久久| 老司机影院成人| 国产亚洲一区二区精品| 欧美极品一区二区三区四区| 精品久久久久久成人av| 国产大屁股一区二区在线视频| 国产高清三级在线| 一级毛片aaaaaa免费看小| 久久鲁丝午夜福利片| 小说图片视频综合网站| 日本与韩国留学比较| 七月丁香在线播放| 久久久久久大精品| 最近中文字幕高清免费大全6| 一本一本综合久久| 久久精品久久久久久久性| 尤物成人国产欧美一区二区三区| 欧美日本视频| 免费av不卡在线播放| 亚洲精品久久久久久婷婷小说 | 丝袜美腿在线中文| 熟女人妻精品中文字幕| 午夜精品国产一区二区电影 | 精品人妻偷拍中文字幕| 哪个播放器可以免费观看大片| 亚洲国产欧美人成| av卡一久久| 韩国av在线不卡| 国产亚洲av片在线观看秒播厂 | 91精品国产九色| 人人妻人人看人人澡| .国产精品久久| 久久久午夜欧美精品| 久久久久精品久久久久真实原创| 免费不卡的大黄色大毛片视频在线观看 | 国产精品综合久久久久久久免费| 26uuu在线亚洲综合色| 丰满乱子伦码专区| 国产精品久久久久久av不卡| 国产人妻一区二区三区在| 国产在线一区二区三区精 | 男人舔奶头视频| 看非洲黑人一级黄片| 日韩大片免费观看网站 | 国产视频内射| 国产高清三级在线| 看十八女毛片水多多多| 一区二区三区高清视频在线| 成人国产麻豆网| 99国产精品一区二区蜜桃av| 免费黄色在线免费观看| 人妻系列 视频| 人妻少妇偷人精品九色| 免费看av在线观看网站| 成人亚洲精品av一区二区| 久久久久久久久中文| 一个人看的www免费观看视频| 日韩欧美精品v在线| 国产av不卡久久| 一区二区三区乱码不卡18| 精品无人区乱码1区二区| 免费人成在线观看视频色| 国产成人精品久久久久久| 美女黄网站色视频| 狠狠狠狠99中文字幕| 国产69精品久久久久777片| 国产激情偷乱视频一区二区| 插逼视频在线观看| 2021天堂中文幕一二区在线观| 国产成人福利小说| 在现免费观看毛片| 久久久精品大字幕| 在线a可以看的网站| 亚洲美女视频黄频| 丝袜美腿在线中文| 久久久国产成人精品二区| 在线播放无遮挡| 欧美日本亚洲视频在线播放| 国产精品日韩av在线免费观看| 成人av在线播放网站| 又爽又黄a免费视频| 三级经典国产精品| 黄片wwwwww| 免费观看人在逋| 国产精品一二三区在线看| 青春草视频在线免费观看| 寂寞人妻少妇视频99o| 国产成年人精品一区二区| 午夜老司机福利剧场| 干丝袜人妻中文字幕| 少妇的逼水好多| 男人舔奶头视频| 丰满乱子伦码专区| 国产69精品久久久久777片| 日韩,欧美,国产一区二区三区 | 99在线视频只有这里精品首页| 成人美女网站在线观看视频| 国产成人精品久久久久久| 久久久久久久久久黄片| 亚洲婷婷狠狠爱综合网| 九九爱精品视频在线观看| 伦理电影大哥的女人| 在线观看一区二区三区| 成人午夜高清在线视频| 久久久久久久国产电影| 黄色配什么色好看| 日韩国内少妇激情av| 天堂中文最新版在线下载 | 99热6这里只有精品| 啦啦啦韩国在线观看视频| 亚洲第一区二区三区不卡| 日韩三级伦理在线观看| 亚洲四区av| 国产色爽女视频免费观看| 成年版毛片免费区| 亚洲av男天堂| 毛片一级片免费看久久久久| 高清在线视频一区二区三区 | 少妇丰满av| 亚洲五月天丁香| 国产精品永久免费网站| 日本黄色视频三级网站网址| 美女大奶头视频| 99久久无色码亚洲精品果冻| 日本三级黄在线观看| 欧美日韩综合久久久久久| 美女黄网站色视频| 午夜激情欧美在线| 日韩成人伦理影院| 国产 一区精品| 视频中文字幕在线观看| 中文字幕免费在线视频6| 国产探花在线观看一区二区| 99久久九九国产精品国产免费| 波多野结衣高清无吗| 高清av免费在线| 国产成人a∨麻豆精品| 99热这里只有精品一区| 国产探花极品一区二区| 亚洲精品久久久久久婷婷小说 | 国产精品久久久久久精品电影小说 | 嫩草影院新地址| 如何舔出高潮| 网址你懂的国产日韩在线| 国内揄拍国产精品人妻在线| 午夜亚洲福利在线播放| 色视频www国产| 亚洲无线观看免费| 蜜臀久久99精品久久宅男| 综合色av麻豆| 国产成人aa在线观看| 少妇的逼水好多| 国产午夜精品久久久久久一区二区三区| 日日干狠狠操夜夜爽| 中文字幕制服av| 伦理电影大哥的女人| 美女黄网站色视频| 最后的刺客免费高清国语| 联通29元200g的流量卡| ponron亚洲| 99在线人妻在线中文字幕| 狂野欧美白嫩少妇大欣赏| 男女那种视频在线观看| 国产乱人视频| 日本爱情动作片www.在线观看| 亚洲av不卡在线观看| 亚洲国产精品合色在线| 精品一区二区免费观看| 日本一二三区视频观看| 精品国内亚洲2022精品成人| 国产av在哪里看| 一个人免费在线观看电影| 亚洲国产精品久久男人天堂| 精品国产露脸久久av麻豆 | 国产毛片a区久久久久| 日日撸夜夜添| 中文乱码字字幕精品一区二区三区 | 亚洲一区高清亚洲精品| 白带黄色成豆腐渣| 丰满人妻一区二区三区视频av| 热99re8久久精品国产| 麻豆成人av视频| 精品一区二区免费观看| 久久久国产成人免费| 麻豆一二三区av精品| 亚洲成人中文字幕在线播放| 91精品伊人久久大香线蕉| 国产在视频线在精品| 成人毛片a级毛片在线播放| 天美传媒精品一区二区| 久久精品国产亚洲av天美| 国产熟女欧美一区二区| 一级毛片电影观看 | h日本视频在线播放| 亚洲18禁久久av| 秋霞在线观看毛片| 国产精品国产高清国产av| 国产精品一区二区三区四区久久| 一级黄片播放器| 国语自产精品视频在线第100页| 久久精品影院6| 男人狂女人下面高潮的视频| 熟女电影av网| 国产成人福利小说| 亚洲久久久久久中文字幕| 成人毛片60女人毛片免费| 搡老妇女老女人老熟妇| 欧美97在线视频| 久久99蜜桃精品久久| 美女内射精品一级片tv| 国产色爽女视频免费观看| 自拍偷自拍亚洲精品老妇| 亚洲精品影视一区二区三区av| av专区在线播放| 欧美日韩在线观看h| 午夜福利在线观看吧| 中文在线观看免费www的网站| 亚洲精品久久久久久婷婷小说 | 中文字幕人妻熟人妻熟丝袜美| 一区二区三区高清视频在线| 69av精品久久久久久| 麻豆一二三区av精品| 国产精品一及| 国产精品国产三级国产专区5o | 卡戴珊不雅视频在线播放| a级一级毛片免费在线观看| 免费播放大片免费观看视频在线观看 | 亚洲va在线va天堂va国产| 日韩成人伦理影院| 欧美日韩精品成人综合77777| 午夜精品在线福利| 久久鲁丝午夜福利片| 精品久久久久久久久av| 在线播放无遮挡| 国产午夜精品久久久久久一区二区三区| 午夜免费激情av| av在线亚洲专区| 国产亚洲午夜精品一区二区久久 | 欧美性猛交黑人性爽| 尤物成人国产欧美一区二区三区| 国产三级在线视频| 日韩高清综合在线| 日本一本二区三区精品| 美女大奶头视频| 国产美女午夜福利| 最近中文字幕2019免费版| 校园人妻丝袜中文字幕| 久久精品熟女亚洲av麻豆精品 | .国产精品久久| 中文字幕熟女人妻在线| 欧美成人午夜免费资源| 中文欧美无线码| 亚洲欧美一区二区三区国产| 久久久久久久久久久免费av| 18禁在线播放成人免费| 国产 一区 欧美 日韩| av国产免费在线观看| 狂野欧美白嫩少妇大欣赏| 久久精品91蜜桃| 亚洲aⅴ乱码一区二区在线播放| 免费一级毛片在线播放高清视频| 日韩精品有码人妻一区| 国产精品美女特级片免费视频播放器| 精品不卡国产一区二区三区| 女人十人毛片免费观看3o分钟| 国产午夜精品一二区理论片| 最近最新中文字幕免费大全7| 超碰97精品在线观看| 边亲边吃奶的免费视频| 久久久久九九精品影院| 国产毛片a区久久久久| 免费看美女性在线毛片视频| 久久精品熟女亚洲av麻豆精品 | 国产免费福利视频在线观看| 毛片一级片免费看久久久久| 日韩成人伦理影院| 国产一区二区在线av高清观看| 内地一区二区视频在线| 日本色播在线视频| 成人性生交大片免费视频hd| 欧美日韩一区二区视频在线观看视频在线 | 国产乱人视频| 午夜福利视频1000在线观看| 中国美白少妇内射xxxbb| 国产男人的电影天堂91| 亚洲国产色片| 欧美不卡视频在线免费观看| 日韩三级伦理在线观看| 精品一区二区三区视频在线| 赤兔流量卡办理| 午夜福利在线在线| 春色校园在线视频观看| 老女人水多毛片| 亚洲无线观看免费| 欧美不卡视频在线免费观看| 村上凉子中文字幕在线| 午夜福利在线观看免费完整高清在| 日韩强制内射视频| 麻豆乱淫一区二区| av在线蜜桃| 色噜噜av男人的天堂激情| 欧美xxxx黑人xx丫x性爽| 国产大屁股一区二区在线视频| 少妇人妻一区二区三区视频| 18禁裸乳无遮挡免费网站照片| 欧美丝袜亚洲另类| 嫩草影院精品99| 狂野欧美白嫩少妇大欣赏| 18禁在线播放成人免费| 一边亲一边摸免费视频| 别揉我奶头 嗯啊视频| 日本免费在线观看一区| 一级毛片久久久久久久久女| 欧美zozozo另类| 男女那种视频在线观看| 男人舔女人下体高潮全视频| 国产不卡一卡二| 亚洲成人中文字幕在线播放| 国产黄色小视频在线观看| 乱码一卡2卡4卡精品| 午夜日本视频在线| 97在线视频观看| 99热6这里只有精品| 久久久久久久久久成人| 日本与韩国留学比较| 成人漫画全彩无遮挡| 欧美极品一区二区三区四区| 亚洲欧美精品专区久久| av卡一久久| 成人二区视频| 免费av不卡在线播放| 最近视频中文字幕2019在线8| 中文字幕制服av| 国产午夜福利久久久久久| 国产高清国产精品国产三级 | 亚洲国产精品久久男人天堂| 亚洲精品亚洲一区二区| 成年女人永久免费观看视频| 国产老妇伦熟女老妇高清| 联通29元200g的流量卡| 波多野结衣高清无吗| 亚洲自拍偷在线| 中文字幕精品亚洲无线码一区| 国产成人freesex在线| 精品午夜福利在线看| 亚洲精品国产av成人精品| 美女高潮的动态| 中文天堂在线官网| 国产精品精品国产色婷婷| 国产久久久一区二区三区| 亚洲美女视频黄频| 蜜桃亚洲精品一区二区三区| 国产午夜精品久久久久久一区二区三区| 美女内射精品一级片tv| 欧美高清成人免费视频www| 成人国产麻豆网| 色播亚洲综合网| 国产乱来视频区| 久热久热在线精品观看| 欧美日韩综合久久久久久| 淫秽高清视频在线观看| 99久久人妻综合| 永久免费av网站大全| 中国国产av一级| 亚洲久久久久久中文字幕| 亚洲av中文字字幕乱码综合| 亚洲av免费在线观看| 国产一级毛片在线| 久久欧美精品欧美久久欧美| 欧美xxxx黑人xx丫x性爽| www.av在线官网国产| 在线免费观看的www视频| 久久久久网色| 三级国产精品欧美在线观看| 麻豆国产97在线/欧美| 99久久精品一区二区三区| 亚洲伊人久久精品综合 | 国内精品宾馆在线| 国产精品电影一区二区三区| 亚洲欧美成人综合另类久久久 | 少妇猛男粗大的猛烈进出视频 | 亚洲一区高清亚洲精品| av播播在线观看一区| 亚洲国产精品成人久久小说| 麻豆久久精品国产亚洲av| av在线播放精品| 大香蕉久久网| 国产成年人精品一区二区| 夜夜爽夜夜爽视频| 3wmmmm亚洲av在线观看| 久久韩国三级中文字幕| 亚洲欧美日韩东京热| a级毛片免费高清观看在线播放| 午夜视频国产福利| 亚洲欧美精品自产自拍| 国产成人aa在线观看| 亚洲中文字幕日韩| 色播亚洲综合网| 国产精品一区二区在线观看99 | 午夜激情福利司机影院| 日产精品乱码卡一卡2卡三| 欧美丝袜亚洲另类| 久久久久国产网址| 国产真实伦视频高清在线观看| 一本一本综合久久| 国产免费一级a男人的天堂| 五月玫瑰六月丁香| 啦啦啦啦在线视频资源| 哪个播放器可以免费观看大片| 国产高清视频在线观看网站| 国产在线一区二区三区精 | 波多野结衣巨乳人妻| 成年女人永久免费观看视频| 又黄又爽又刺激的免费视频.| 久久久久久伊人网av| 欧美3d第一页| 日韩欧美精品v在线| 国产片特级美女逼逼视频| 国产伦在线观看视频一区| 综合色丁香网| 国产淫语在线视频| 欧美性猛交黑人性爽| 极品教师在线视频| 亚洲怡红院男人天堂| 晚上一个人看的免费电影| 日韩欧美精品v在线| 51国产日韩欧美| 午夜福利在线在线| 精品欧美国产一区二区三| 我的老师免费观看完整版| 欧美激情在线99| 国语自产精品视频在线第100页| 91在线精品国自产拍蜜月| 国产精品乱码一区二三区的特点| 日韩三级伦理在线观看| 国产精品久久久久久久久免| 国产伦精品一区二区三区视频9|