• <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ī)律
    国产精品国产三级国产av玫瑰| 嘟嘟电影网在线观看| 777米奇影视久久| 日本黄大片高清| 亚洲av欧美aⅴ国产| 日本wwww免费看| 日韩av不卡免费在线播放| 在线看a的网站| 亚洲欧美一区二区三区黑人 | 国产成人aa在线观看| 插阴视频在线观看视频| 最近手机中文字幕大全| 久久6这里有精品| 久久久久视频综合| 在线观看免费日韩欧美大片 | 国产又色又爽无遮挡免| 九九爱精品视频在线观看| 妹子高潮喷水视频| 91在线精品国自产拍蜜月| 亚洲av成人精品一二三区| 看十八女毛片水多多多| 国产日韩欧美亚洲二区| 老司机影院成人| av天堂久久9| 69精品国产乱码久久久| 国产午夜精品一二区理论片| 下体分泌物呈黄色| 国产亚洲最大av| 久久99一区二区三区| 国产av一区二区精品久久| 韩国av在线不卡| 亚洲精品国产色婷婷电影| 街头女战士在线观看网站| 人妻制服诱惑在线中文字幕| 国产成人一区二区在线| 国产精品人妻久久久久久| 国产日韩欧美在线精品| 午夜老司机福利剧场| 免费不卡的大黄色大毛片视频在线观看| 精品国产一区二区三区久久久樱花| 纯流量卡能插随身wifi吗| 综合色丁香网| 亚洲综合色惰| 免费观看在线日韩| 久久毛片免费看一区二区三区| 午夜激情福利司机影院| 9色porny在线观看| 国产欧美日韩综合在线一区二区 | 久久99蜜桃精品久久| 国产精品99久久久久久久久| 国产精品久久久久久精品古装| 国产高清有码在线观看视频| 亚洲av男天堂| 国产精品一区二区性色av| 国产高清三级在线| 国产成人aa在线观看| 欧美日韩视频高清一区二区三区二| 免费人妻精品一区二区三区视频| 午夜福利,免费看| 亚洲人成网站在线播| 最近手机中文字幕大全| av不卡在线播放| 中国国产av一级| 男女边吃奶边做爰视频| 国产精品福利在线免费观看| 免费看光身美女| 极品少妇高潮喷水抽搐| 国产美女午夜福利| 亚洲情色 制服丝袜| 欧美bdsm另类| 欧美激情极品国产一区二区三区 | 街头女战士在线观看网站| 五月开心婷婷网| 一级毛片aaaaaa免费看小| av女优亚洲男人天堂| 性色avwww在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲第一区二区三区不卡| 99久久精品一区二区三区| 欧美成人精品欧美一级黄| 亚洲精品日本国产第一区| 国产精品女同一区二区软件| 亚洲欧美成人综合另类久久久| 夜夜看夜夜爽夜夜摸| 亚洲精品一二三| 久久久a久久爽久久v久久| 狂野欧美白嫩少妇大欣赏| 视频中文字幕在线观看| 一本色道久久久久久精品综合| 大话2 男鬼变身卡| 国产精品99久久久久久久久| 国产高清三级在线| 亚洲av综合色区一区| 亚洲va在线va天堂va国产| 亚洲情色 制服丝袜| 18禁裸乳无遮挡动漫免费视频| 国产黄色免费在线视频| 人体艺术视频欧美日本| 熟女av电影| 久久久国产欧美日韩av| 国产精品国产三级国产专区5o| 少妇猛男粗大的猛烈进出视频| 精品午夜福利在线看| 久久99精品国语久久久| 免费av不卡在线播放| 在线观看免费视频网站a站| 国产无遮挡羞羞视频在线观看| 黄色视频在线播放观看不卡| 王馨瑶露胸无遮挡在线观看| 欧美性感艳星| 成人影院久久| 久久99热这里只频精品6学生| 一本色道久久久久久精品综合| 天天躁夜夜躁狠狠久久av| 久久毛片免费看一区二区三区| 久热这里只有精品99| 两个人免费观看高清视频 | 国产av精品麻豆| 精品视频人人做人人爽| 国产精品三级大全| 国产亚洲5aaaaa淫片| 国产精品一区二区三区四区免费观看| 国产在线一区二区三区精| 18禁动态无遮挡网站| 精品人妻熟女av久视频| 插阴视频在线观看视频| 又大又黄又爽视频免费| 国内少妇人妻偷人精品xxx网站| 黄色一级大片看看| 如何舔出高潮| 国产精品久久久久久精品电影小说| 最后的刺客免费高清国语| 国产成人免费观看mmmm| 最近最新中文字幕免费大全7| 国产视频首页在线观看| 99久国产av精品国产电影| 涩涩av久久男人的天堂| 久久久久久久亚洲中文字幕| 国产欧美亚洲国产| 久久99一区二区三区| 日韩大片免费观看网站| av在线播放精品| 亚洲精品日本国产第一区| 热re99久久精品国产66热6| 久久精品国产a三级三级三级| 在线观看免费视频网站a站| 久久久久久久久久久免费av| 久久99一区二区三区| 精品国产露脸久久av麻豆| 欧美精品人与动牲交sv欧美| 视频区图区小说| 各种免费的搞黄视频| 黄色怎么调成土黄色| 女性被躁到高潮视频| 国产欧美日韩一区二区三区在线 | 99九九在线精品视频 | 国产精品一区二区三区四区免费观看| 尾随美女入室| 国产伦理片在线播放av一区| 亚洲精品成人av观看孕妇| 国产91av在线免费观看| 少妇被粗大的猛进出69影院 | 亚洲av福利一区| 日韩三级伦理在线观看| 国产91av在线免费观看| a 毛片基地| 人妻系列 视频| 亚洲国产成人一精品久久久| 亚洲欧美日韩另类电影网站| 国产乱来视频区| 狂野欧美白嫩少妇大欣赏| 热99国产精品久久久久久7| 欧美3d第一页| 亚洲第一av免费看| 精品久久久噜噜| 国产伦精品一区二区三区四那| 国产欧美亚洲国产| 精品人妻偷拍中文字幕| 国产精品国产三级专区第一集| 如日韩欧美国产精品一区二区三区 | 日本免费在线观看一区| 久久久午夜欧美精品| 色婷婷av一区二区三区视频| 国产欧美日韩综合在线一区二区 | av国产精品久久久久影院| 99久国产av精品国产电影| 国产91av在线免费观看| 两个人免费观看高清视频 | 精品一区在线观看国产| 亚洲欧洲日产国产| 日本av手机在线免费观看| 国产成人a∨麻豆精品| 久久青草综合色| 国产精品三级大全| 亚洲欧洲日产国产| 国产精品女同一区二区软件| 一区二区三区四区激情视频| 国产精品久久久久久av不卡| av国产精品久久久久影院| 大片免费播放器 马上看| 亚洲成人手机| 一本久久精品| 国产精品一区www在线观看| 建设人人有责人人尽责人人享有的| 国产69精品久久久久777片| 久久久久久久大尺度免费视频| av线在线观看网站| 精品一区二区三卡| 大香蕉久久网| 大陆偷拍与自拍| 天天操日日干夜夜撸| 亚洲av中文av极速乱| 只有这里有精品99| 久久精品国产自在天天线| 日本av免费视频播放| 亚洲美女视频黄频| 国产免费福利视频在线观看| 亚洲av二区三区四区| 秋霞伦理黄片| 国产成人freesex在线| 免费av不卡在线播放| 国产精品.久久久| 午夜福利在线观看免费完整高清在| 久久久亚洲精品成人影院| 国产亚洲5aaaaa淫片| 国产熟女午夜一区二区三区 | 国产成人freesex在线| 欧美日韩综合久久久久久| 婷婷色麻豆天堂久久| 91aial.com中文字幕在线观看| 国产精品国产三级国产av玫瑰| 国产成人91sexporn| 一级片'在线观看视频| 人人妻人人添人人爽欧美一区卜| 男男h啪啪无遮挡| 最近手机中文字幕大全| 免费看av在线观看网站| 新久久久久国产一级毛片| 人妻 亚洲 视频| 亚洲经典国产精华液单| 人妻夜夜爽99麻豆av| 久久99精品国语久久久| 亚洲国产精品专区欧美| 偷拍熟女少妇极品色| 天美传媒精品一区二区| 一区二区三区乱码不卡18| 欧美日韩av久久| 另类精品久久| av在线app专区| 午夜视频国产福利| 国产精品一区二区三区四区免费观看| 寂寞人妻少妇视频99o| 国产 一区精品| 男女边吃奶边做爰视频| 国精品久久久久久国模美| 国产老妇伦熟女老妇高清| 国产欧美日韩一区二区三区在线 | 国产精品一二三区在线看| 新久久久久国产一级毛片| 亚洲,欧美,日韩| 国产高清不卡午夜福利| 三级经典国产精品| 欧美精品人与动牲交sv欧美| 3wmmmm亚洲av在线观看| √禁漫天堂资源中文www| 免费av中文字幕在线| 一二三四中文在线观看免费高清| 中文乱码字字幕精品一区二区三区| 国产69精品久久久久777片| 国产乱人偷精品视频| 久久av网站| 国产成人精品久久久久久| 一区二区av电影网| 九色成人免费人妻av| 久久韩国三级中文字幕| 狂野欧美激情性bbbbbb| 久久精品国产亚洲av天美| 国产午夜精品一二区理论片| 色94色欧美一区二区| 欧美亚洲 丝袜 人妻 在线| 日本猛色少妇xxxxx猛交久久| 久久6这里有精品| 国产中年淑女户外野战色| 少妇精品久久久久久久| 欧美日本中文国产一区发布| 亚洲av欧美aⅴ国产| 搡女人真爽免费视频火全软件| 噜噜噜噜噜久久久久久91| 日韩欧美精品免费久久| 在线观看www视频免费| 久久狼人影院| 99视频精品全部免费 在线| 亚洲欧美精品专区久久| 国产精品三级大全| 久久韩国三级中文字幕| 99热这里只有是精品在线观看| 又黄又爽又刺激的免费视频.| 成年人免费黄色播放视频 | 菩萨蛮人人尽说江南好唐韦庄| 人人澡人人妻人| 中文乱码字字幕精品一区二区三区| 国产在线一区二区三区精| 国产毛片在线视频| 天天躁夜夜躁狠狠久久av| 一级二级三级毛片免费看| 三上悠亚av全集在线观看 | 国产成人精品婷婷| 美女内射精品一级片tv| 在线亚洲精品国产二区图片欧美 | 精品亚洲成a人片在线观看| 日本黄色日本黄色录像| 2022亚洲国产成人精品| 国产欧美另类精品又又久久亚洲欧美| 高清不卡的av网站| 久久99热6这里只有精品| 夜夜爽夜夜爽视频| .国产精品久久| 69精品国产乱码久久久| 91久久精品国产一区二区三区| 欧美精品亚洲一区二区| 丝袜喷水一区| 色视频在线一区二区三区| 亚洲精品色激情综合| 精品卡一卡二卡四卡免费| 成年av动漫网址| 高清视频免费观看一区二区| 国产免费福利视频在线观看| 亚洲怡红院男人天堂| 国产av国产精品国产| 青春草国产在线视频| 大又大粗又爽又黄少妇毛片口| 日韩欧美 国产精品| 丰满人妻一区二区三区视频av| 国产成人a∨麻豆精品| 国产精品久久久久久久久免| 91精品一卡2卡3卡4卡| 国产永久视频网站| 亚洲精品乱久久久久久| 六月丁香七月| 午夜精品国产一区二区电影| 中文资源天堂在线| 亚洲激情五月婷婷啪啪| 精品久久久久久久久av| 99久久精品国产国产毛片| tube8黄色片| av福利片在线观看| 欧美bdsm另类| 91在线精品国自产拍蜜月| 婷婷色av中文字幕| 在线精品无人区一区二区三| 精品人妻偷拍中文字幕| 日本猛色少妇xxxxx猛交久久| 男人爽女人下面视频在线观看| 一个人免费看片子| 成人亚洲欧美一区二区av| 99热国产这里只有精品6| 亚洲精品久久午夜乱码| 久久99热6这里只有精品| 精品少妇黑人巨大在线播放| 这个男人来自地球电影免费观看 | 少妇精品久久久久久久| 国产黄片美女视频| 国产欧美日韩综合在线一区二区 | 建设人人有责人人尽责人人享有的| 亚洲丝袜综合中文字幕| 91久久精品电影网| 777米奇影视久久| 热99国产精品久久久久久7| 国产成人精品一,二区| 欧美三级亚洲精品| 丰满人妻一区二区三区视频av| 亚洲国产最新在线播放| 中文在线观看免费www的网站| 国产高清不卡午夜福利| 成人国产av品久久久| 噜噜噜噜噜久久久久久91| 97在线人人人人妻| 街头女战士在线观看网站| 日韩亚洲欧美综合| 日本黄色日本黄色录像| 日韩av免费高清视频| 三上悠亚av全集在线观看 | a 毛片基地| av有码第一页| 国产欧美日韩综合在线一区二区 | 啦啦啦啦在线视频资源| 麻豆成人av视频| 欧美日韩精品成人综合77777| 2021少妇久久久久久久久久久| 99久久精品一区二区三区| 日韩人妻高清精品专区| 伦理电影免费视频| a级毛色黄片| 国产在线视频一区二区| 国内揄拍国产精品人妻在线| av卡一久久| 久久久a久久爽久久v久久| 大话2 男鬼变身卡| 国产男女内射视频| 欧美日韩av久久| 人妻 亚洲 视频| 国产极品天堂在线| 天美传媒精品一区二区| 日日啪夜夜爽| 天堂中文最新版在线下载| 一级片'在线观看视频| 99热6这里只有精品| 人妻一区二区av| 国产极品粉嫩免费观看在线 | 国产 一区精品| 国产成人免费无遮挡视频| 国产乱人偷精品视频| 国产免费福利视频在线观看| 高清av免费在线| 久久精品国产亚洲网站| 内射极品少妇av片p| 99久久精品热视频| 午夜福利影视在线免费观看| 精品视频人人做人人爽| 亚洲熟女精品中文字幕| 久久国产亚洲av麻豆专区| 丁香六月天网| 中文天堂在线官网| 国产日韩一区二区三区精品不卡 | 少妇的逼好多水| a级毛色黄片| 久久久国产精品麻豆| 国产精品国产三级国产av玫瑰| av女优亚洲男人天堂| 免费看不卡的av| 3wmmmm亚洲av在线观看| 亚洲va在线va天堂va国产| 两个人免费观看高清视频 | 久久精品夜色国产| 看十八女毛片水多多多| 七月丁香在线播放| 麻豆成人午夜福利视频| 久久av网站| 国产一区二区在线观看av| 亚洲av免费高清在线观看| 久热这里只有精品99| 日韩伦理黄色片| 能在线免费看毛片的网站| 国产色婷婷99| 亚洲欧美精品专区久久| 亚洲成人手机| 人妻 亚洲 视频| 一区在线观看完整版| 亚洲欧美清纯卡通| 国产精品国产av在线观看| 26uuu在线亚洲综合色| 性色avwww在线观看| 国产又色又爽无遮挡免| 亚洲av成人精品一二三区| 国产成人免费无遮挡视频| 一级爰片在线观看| 亚洲真实伦在线观看| 国产免费一区二区三区四区乱码| 青春草视频在线免费观看| 日韩中字成人| 男人和女人高潮做爰伦理| 国产伦理片在线播放av一区| 久久久久精品久久久久真实原创| 在现免费观看毛片| 欧美+日韩+精品| 一级毛片 在线播放| 在线免费观看不下载黄p国产| 大香蕉久久网| 精品久久国产蜜桃| 99热这里只有精品一区| 久久久久精品久久久久真实原创| 久久99热6这里只有精品| 中文欧美无线码| 午夜激情福利司机影院| 成人毛片a级毛片在线播放| 国产淫片久久久久久久久| 亚洲经典国产精华液单| 日本猛色少妇xxxxx猛交久久| 婷婷色av中文字幕| 日本-黄色视频高清免费观看| 免费大片18禁| 一级爰片在线观看| 日韩伦理黄色片| 人妻少妇偷人精品九色| 尾随美女入室| 十分钟在线观看高清视频www | 丁香六月天网| 只有这里有精品99| 亚洲va在线va天堂va国产| 尾随美女入室| 日韩三级伦理在线观看| 免费看日本二区| 中文天堂在线官网| 国产免费视频播放在线视频| 日韩电影二区| 一区二区av电影网| 国产视频内射| 午夜福利在线观看免费完整高清在| 天美传媒精品一区二区| 一级爰片在线观看| 精华霜和精华液先用哪个| 曰老女人黄片| 亚洲成人av在线免费| 国产精品女同一区二区软件| 又黄又爽又刺激的免费视频.| 日韩中字成人| 国产欧美日韩一区二区三区在线 | 久久久久久久久久人人人人人人| 成人18禁高潮啪啪吃奶动态图 | 久久国内精品自在自线图片| 亚洲国产日韩一区二区| 亚洲av福利一区| 亚洲美女搞黄在线观看| 欧美亚洲 丝袜 人妻 在线| 国产亚洲精品久久久com| 久久久久国产网址| a级毛片免费高清观看在线播放| 97在线人人人人妻| 久久99蜜桃精品久久| 老司机影院成人| 久久女婷五月综合色啪小说| 国产精品女同一区二区软件| 免费少妇av软件| 99久久精品热视频| 欧美成人精品欧美一级黄| 91午夜精品亚洲一区二区三区| 最后的刺客免费高清国语| 国产成人freesex在线| 黑人猛操日本美女一级片| 丝袜喷水一区| 亚洲欧美成人综合另类久久久| 欧美精品亚洲一区二区| 老司机影院毛片| 久久精品久久精品一区二区三区| 国产av一区二区精品久久| av免费在线看不卡| 欧美成人精品欧美一级黄| 亚洲精品乱码久久久久久按摩| 久久国内精品自在自线图片| 国产欧美亚洲国产| 亚洲国产精品一区三区| 青春草视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 晚上一个人看的免费电影| 只有这里有精品99| 国产精品99久久99久久久不卡 | 久久人人爽av亚洲精品天堂| 国产欧美日韩精品一区二区| 99久久中文字幕三级久久日本| 高清毛片免费看| 国产欧美日韩综合在线一区二区 | 青春草视频在线免费观看| 日韩中文字幕视频在线看片| 久久精品国产a三级三级三级| 日本午夜av视频| 国产视频内射| a级毛片在线看网站| 免费观看的影片在线观看| 亚洲av二区三区四区| 久久久久久久精品精品| 少妇被粗大的猛进出69影院 | 欧美日韩国产mv在线观看视频| 欧美人与善性xxx| 国产视频首页在线观看| 国产淫语在线视频| 我的女老师完整版在线观看| 一级av片app| 国产精品偷伦视频观看了| 亚洲国产色片| 精品熟女少妇av免费看| 美女内射精品一级片tv| 黄色配什么色好看| 综合色丁香网| 色5月婷婷丁香| 伦理电影免费视频| 两个人免费观看高清视频 | 国产深夜福利视频在线观看| 深夜a级毛片| 在线免费观看不下载黄p国产| 天天操日日干夜夜撸| 国产精品久久久久成人av| 大陆偷拍与自拍| 国模一区二区三区四区视频| 女人久久www免费人成看片| 免费观看性生交大片5| 麻豆成人av视频| 搡老乐熟女国产| 久久国产精品男人的天堂亚洲 | 乱人伦中国视频| 26uuu在线亚洲综合色| 18禁动态无遮挡网站| av播播在线观看一区| 三级国产精品片| 最近中文字幕2019免费版| 亚洲欧洲日产国产| 美女中出高潮动态图| 精品国产一区二区久久| 亚洲欧美中文字幕日韩二区| 欧美区成人在线视频| 啦啦啦中文免费视频观看日本| 亚洲,欧美,日韩| 午夜免费鲁丝| 久久精品国产亚洲av涩爱| 青青草视频在线视频观看| 国产毛片在线视频| 青青草视频在线视频观看| 高清视频免费观看一区二区| 国产精品女同一区二区软件| 亚洲性久久影院| 综合色丁香网| 国产在线视频一区二区| 最近中文字幕高清免费大全6| 亚洲久久久国产精品| 色视频在线一区二区三区| 国产成人freesex在线|