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

    圓管壁粗糙度對(duì)圓管內(nèi)流態(tài)影響研究

    2016-05-17 03:33:05曾亞武方海龍

    李 歡,曾亞武,夏 磊,方海龍

    (武漢大學(xué) 土木建筑工程學(xué)院, 湖北 武漢 430072)

    ?

    圓管壁粗糙度對(duì)圓管內(nèi)流態(tài)影響研究

    李歡,曾亞武,夏磊,方海龍

    (武漢大學(xué) 土木建筑工程學(xué)院, 湖北 武漢 430072)

    摘要:基于尼古拉茲圓管實(shí)驗(yàn)思想,利用顆粒流程序中的流體計(jì)算模塊,實(shí)現(xiàn)了對(duì)圓管壁粗糙度的模擬,進(jìn)而分析了不同粗糙度對(duì)圓管斷面流態(tài)分布的影響以及在不同壓力差作用下平均流速與圓管壁粗糙度之間的關(guān)系。研究結(jié)果表明:圓管內(nèi)流體流速受管壁的擾動(dòng)影響,在橫截面上呈U型分布,而非拋物線型分布;圓管壁對(duì)管內(nèi)流體流態(tài)的影響范圍隨管壁粗糙度變化而變化,圓管管壁相對(duì)粗糙度越大,管內(nèi)受擾動(dòng)流體范圍越大,即管內(nèi)層流范圍越??;在管壁粗糙度一定的情況下,圓管斷面流量隨上下端面壓力差的增大而增大;在相同壓力差作用下,圓管斷面流量隨管壁相對(duì)粗糙度增大而逐漸減小,顆粒的阻流作用明顯;隨著管壁粗糙度的增大,圓管上下端面受壓力差作用的影響范圍也隨之增大,使管內(nèi)流體流態(tài)變得更加復(fù)雜。

    關(guān)鍵詞:尼古拉茲實(shí)驗(yàn);粗糙圓管;顆粒流流體計(jì)算;相對(duì)粗糙度;流體流態(tài)

    1885年雷諾(Reynolds O)曾用試驗(yàn)揭示了流體運(yùn)動(dòng)存在著兩種流態(tài):層流和紊流,其試驗(yàn)結(jié)果表明:流態(tài)不同,沿程水頭損失的規(guī)律也不同。達(dá)西[1](Darcy H)根據(jù)不同尺寸的圓管和不同類(lèi)型、不同長(zhǎng)度的土樣進(jìn)行試驗(yàn)得到了線性滲透定律:在層流時(shí),沿程水頭損失與流速的一次方成比例。1901年福希海梅(Forchheimer)首先提出滲流的高雷諾數(shù)非線性關(guān)系:在充分紊流時(shí),沿程水頭損失與流速的二次方成比例。在過(guò)渡流態(tài),其沿程阻力水頭損失與流速的1.75~2次方成比例。由此可見(jiàn),處于層流、過(guò)渡(層流與紊流間過(guò)渡態(tài))和充分紊流流態(tài)的水流遵循不同的沿程阻力規(guī)律[2]。

    為了探討紊流存在時(shí)流體的沿程阻力規(guī)律,尼古拉茲[3](Nikuradse J)在圓管內(nèi)壁粘貼上經(jīng)過(guò)篩分具有不同粒徑的砂礫,制成人工均勻顆粒粗糙圓管;然后在不同粗糙度的圓管上進(jìn)行系統(tǒng)實(shí)驗(yàn),實(shí)驗(yàn)結(jié)果認(rèn)為沿程阻力系數(shù)取決于雷諾數(shù)和相對(duì)粗糙度。黃琨、萬(wàn)軍偉等[4-5]對(duì)尼古拉茲實(shí)驗(yàn)的數(shù)據(jù)進(jìn)行整理時(shí)指出尼古拉茲實(shí)驗(yàn)存在誤差,即尼古拉茲實(shí)驗(yàn)二次刷膠引起的實(shí)驗(yàn)誤差。喬昌凱等[6]用玻璃球在矩形水槽中鋪設(shè)不同層數(shù)的透水層模擬卵石河床,得到不同厚度透水層阻力系數(shù)隨雷諾數(shù)變化的關(guān)系。Herwig H等[7]研究表明,在層流流動(dòng)中壁面粗糙度的影響也是值得重視的并對(duì)粗糙度參數(shù)進(jìn)行了定義。郝鵬飛等[8]研究了粗糙度對(duì)矩形截面微管道內(nèi)液體流動(dòng)阻力特性的影響,采用微觀粒子圖像測(cè)速技術(shù)測(cè)量了粗糙微管道內(nèi)的流場(chǎng)結(jié)構(gòu)。隨著計(jì)算機(jī)技術(shù)和數(shù)值計(jì)算方法的發(fā)展,使得利用數(shù)值方法研究流體滲流特性、流場(chǎng)結(jié)構(gòu)等成為可能。周健[9],劉洋[10-12]等利用顆粒流程序研究了砂土液化、管涌、砂沸等多孔介質(zhì)流體滲流現(xiàn)象的微觀機(jī)理,結(jié)果表明顆粒流程序能較好地模擬流體和顆粒的相互作用,并重點(diǎn)研究了流體對(duì)顆粒的影響,但對(duì)流體與流場(chǎng)分布研究較少。

    尼古拉茲(Nikuradse)實(shí)驗(yàn)提出了粗糙度對(duì)流態(tài)的影響,進(jìn)而得出流體沿程阻力與雷諾數(shù)和相對(duì)粗糙度的關(guān)系,但受試驗(yàn)觀測(cè)條件的限制,試驗(yàn)中并沒(méi)有開(kāi)展圓管內(nèi)流場(chǎng)結(jié)構(gòu)的研究。本文基于尼古拉茲(Nikuradse)圓管實(shí)驗(yàn)思想,利用顆粒流程序中的流體計(jì)算模塊,實(shí)現(xiàn)對(duì)圓管壁粗糙度的模擬,進(jìn)而分析不同粗糙度對(duì)圓管斷面流態(tài)分布的影響以及在不同壓力梯度下平均流速與圓管壁粗糙度之間的關(guān)系。

    1顆粒流CFD計(jì)算原理

    1.1流體運(yùn)動(dòng)方程

    對(duì)于單位體積固液兩相系統(tǒng)的常密度液相不可壓縮流的Navier-Stokes方程可表示為:

    (1)

    (2)

    式中:u為流速矢量;τ為黏性應(yīng)力張量;g為重力加速度;fint為單位體積內(nèi)顆粒與流體的相互作用力;n為孔隙率;ρf為流體的密度。為拉普拉斯(Laplace)算子。

    1.2流體-顆粒相互作用

    流體-顆粒相互作用拽曳力可以用一些半經(jīng)驗(yàn)公式計(jì)算。本文的研究即采用Ergun S[13]的孔隙率小于0.8(n≤0.8)的計(jì)算公式:

    (3)

    Ergun方程涵蓋了很大范圍的流體條件,適用于能量非線性耗散的高雷諾系數(shù)的孔隙流,對(duì)于孔隙率大于0.8(n≥0.8)的情況,可以應(yīng)用Wen C Y等[14]的公式:

    (4)

    其中球形顆粒的拖拽力系數(shù)CD是雷諾數(shù)Re的方程,如下

    (5)

    通過(guò)式(3)和式(4)求得顆粒與流體之間的相互作用力,然后代入式(2)求解Navier-Stokes方程。在顆粒流流固耦合計(jì)算內(nèi)核是采用由Patanker和Spalding提出SIMPLE算法,即壓力耦合方程組的半隱式方法,其核心是采用“猜測(cè)—修正的方法在交錯(cuò)網(wǎng)格的基礎(chǔ)上計(jì)算流體壓力場(chǎng),進(jìn)而求解流體動(dòng)量方程(N-S方程)。

    2粗糙圓管模擬實(shí)驗(yàn)

    2.1建立圓管粗糙度模型

    圖1尼古拉茲(NikuradseJ)表面粗糙度圓管示意圖

    本文基于尼古拉茲(NikuradseJ)人工粗糙管實(shí)驗(yàn)的思想,在圓管中不斷增加沙粒的粒徑直至極限狀態(tài),從而建立不同相對(duì)粗糙度圓管模型。

    首先在顆粒流程序中建立光滑圓管模型,然后在圓管內(nèi)壁粘結(jié)上(通過(guò)算法固定)一定直徑的顆粒,模擬一定的粗糙度,然后增大顆粒的半徑而達(dá)到極限粗糙度狀態(tài),并對(duì)一系列的模型進(jìn)行流體模擬實(shí)驗(yàn)。

    對(duì)于人工粗糙管的模擬,本文嘗試了兩種不同的方法:① 隨機(jī)粗糙度模型,② 規(guī)則粗糙度模型。由于隨機(jī)法生成的模型的相對(duì)粗糙度難以定量化,故最終采用規(guī)則粗糙度模型。

    規(guī)則粗糙度模型建模分為以下兩步:

    (1) 建立圓管直徑,以及絕對(duì)粗糙度,還有圓管壁上顆粒半徑等函數(shù)關(guān)系表達(dá)式;

    (2) 推導(dǎo)出規(guī)則陣列顆粒的圓心表達(dá),在顆粒流程序中直接用FISH函數(shù)生成相應(yīng)的顆粒。

    (6)

    式中:θ為環(huán)向弧度增量。

    圖2相對(duì)粗糙度與圓管直徑關(guān)系圖

    表1 不同粗糙度圓管模型建模參數(shù)

    注:參數(shù)推導(dǎo)時(shí),取圓管的直徑為70mm,軸線長(zhǎng)為200mm。其中模型1為光滑圓管模型,模型5為極限粗糙度模型。

    通過(guò)顆粒流程序編制Fish函數(shù)并且按照表1的參數(shù)建立的圓管粗糙度模型如圖3所示。圖3中RRT為轉(zhuǎn)置粗糙度(RelativeRoughnessTranspose,RRT),與相對(duì)粗糙度RR為倒數(shù)關(guān)系。在Δ/d=0時(shí),對(duì)應(yīng)于d/Δ=∞,為表示方便,仍以d/Δ=0表示。由圖3可見(jiàn),模型1的相對(duì)粗糙度為零,為光滑圓管;模型5為極限粗糙圓管,其轉(zhuǎn)置粗糙度d/Δ=2.4;模型4、模型3、和模型2的轉(zhuǎn)置粗糙度d/Δ分別為模型5(極限轉(zhuǎn)置粗糙度)的3倍,6倍,12.5倍。其中d/Δ=0與d/Δ=30是尼古拉茲實(shí)驗(yàn)人工粗糙圓管采用的模型。

    2.2實(shí)驗(yàn)?zāi)M過(guò)程

    模擬實(shí)驗(yàn)采用滲透實(shí)驗(yàn)中達(dá)西實(shí)驗(yàn)?zāi)P?即在圓管上頂面和下底面加壓)。對(duì)不同粗糙度的圓管分別在頂部和底部施加流體壓力邊界條件,使之產(chǎn)生壓力差,在不同的水力梯度條件下進(jìn)行PFC-CFD耦合運(yùn)算,使試樣最終達(dá)到穩(wěn)定的流態(tài),平衡穩(wěn)定的標(biāo)準(zhǔn)是計(jì)算所得的通過(guò)圓管截面的流量趨于某一定值。如圖所示其中流體的密度ρf=1 000kg/m3,流體動(dòng)力黏滯系數(shù)μf=10-3Pa.s,圓管和顆粒的摩擦系數(shù)為0.5,底部和頂部的壓力差依次為P1=2 000Pa,P2=4 000Pa,P3=6 000Pa。圓管上下頂面為壓力邊界,而圓管側(cè)邊界為滑動(dòng)邊界。

    3模擬結(jié)果分析

    模擬實(shí)驗(yàn)中測(cè)定各模型在圓管頂面(Z=200mm)處的橫截面流量和XOZ縱剖面的流速分布。其中流量通過(guò)顆粒流程序編制的FISH函數(shù)直接計(jì)算,用以分析相對(duì)粗糙度與流量之間的關(guān)系;XOZ縱剖面流速分布在三維模型中不易觀察,故利用FISH函數(shù)導(dǎo)出數(shù)據(jù),在MATLAB中繪制縱剖面速度云圖,觀察縱剖面內(nèi)的流速分布,用以分析粗糙度對(duì)圓管內(nèi)流速分布的影響。

    圖3不同相對(duì)粗糙度圓管模型

    3.1流量與時(shí)間關(guān)系

    首先計(jì)算不同壓力差下,圓管頂面截面流量與時(shí)間的關(guān)系,如圖4所示。為了便于比較,本文統(tǒng)一設(shè)置流體時(shí)步長(zhǎng)為1×10-5s,設(shè)置循環(huán)至平衡時(shí)間為0.6s,流體的計(jì)算時(shí)步為60 000步,而力學(xué)計(jì)算步長(zhǎng)由程序自動(dòng)控制,一般小于或等于流體計(jì)算時(shí)步的步長(zhǎng)。

    首先,從圖4(a)~圖4(c)可以看出,圓管頂面截面流量都是先隨著時(shí)間緩慢增加,然后逐漸趨于某一定值而達(dá)到平衡,說(shuō)明流體運(yùn)動(dòng)從加壓到平衡需要一定的時(shí)間,體現(xiàn)在數(shù)值模擬中,就是需要經(jīng)過(guò)若干時(shí)步的計(jì)算,才能達(dá)到平衡。

    其次,從圖4(a)~圖4(c)還可以看出,對(duì)于不同的圓管粗糙度,流量曲線達(dá)到平衡后的穩(wěn)定流量與管壁相對(duì)粗糙度成反比例關(guān)系,即隨著管壁相對(duì)粗糙度的增大,圓管截面流量逐漸減小。這一規(guī)律在不同的壓力差下均成立。分析原因,主要有二:其一是粗糙顆粒的存在,阻塞了流體的流動(dòng)路徑,特別是極限粗糙度的情況下,顆粒的阻塞作用最明顯,流量也最??;其二,顆粒的存在,使過(guò)流凈斷面面積減小,從而使通過(guò)圓管截面的流量減小。另外從圖中也可以觀測(cè)到,在截面流量達(dá)到平衡之前,流量-時(shí)間曲線的斜率是隨著相對(duì)粗糙度的增加而減小的,說(shuō)明管壁粗糙度對(duì)管內(nèi)流體的流態(tài)存在影響,使得相對(duì)粗糙度越大,圓管截面流量達(dá)到平衡時(shí)需要經(jīng)過(guò)更多的計(jì)算時(shí)步。

    最后,從圖4(a)~圖4(c)中可以看到,壓力差不斷增大時(shí),縱軸坐標(biāo)的數(shù)值在增大,那么反映通過(guò)相同粗糙度圓管的流量也不斷增大。為了說(shuō)明這種現(xiàn)象;特別在圖4(d)中列出了光滑圓管在不同壓力差下的流量-時(shí)間關(guān)系曲線,從圖4(d)中可以明顯看到隨著壓力差的增大,通過(guò)圓管的流量是逐漸增大的。

    3.2圓管粗糙度、壓力差與流量之間的關(guān)系

    圖4 不同壓力差下流體流量與時(shí)間的關(guān)系曲線

    圖5圓管粗糙度、壓力差與流量之間的關(guān)系

    由圖5可見(jiàn),當(dāng)圓管粗糙度一定時(shí),隨著圓管上下斷面水力壓差增大,圓管斷面流量(流速)也相應(yīng)增大,與達(dá)西定律表現(xiàn)出一致性;當(dāng)圓管上下斷面壓差一定時(shí),隨著圓管相對(duì)粗糙度逐步增大,斷面流量(流速)則逐漸減小,說(shuō)明圓管壁粗糙度通過(guò)影響斷面過(guò)流面積和流體流態(tài),進(jìn)而影響斷面流量(流速)。

    3.3斷面流速分布圖

    圖6(a)~圖6(e)分別表示在相同壓力差P=2000Pa作用下,不同粗糙度圓管XOZ平面的穩(wěn)定流速分布圖。

    由圖6(a)~圖6(e)可以看出:

    (1) 在XOZ平面中,流速是關(guān)于圓管軸線對(duì)稱(chēng)分布的,這一分布符合模型邊界條件和幾何條件軸對(duì)稱(chēng)的特點(diǎn)。

    (2) 在XOZ平面中,靠近軸線的區(qū)域流速分布均勻,而靠近管壁部分流速分布則不均勻,且隨著圓管壁相對(duì)粗糙度的增大,圓管壁附近受影響的區(qū)域逐步擴(kuò)大,即圓管軸線附近層流區(qū)域逐步減小。如圖6(a)所示的光滑圓管流速分布圖,圓管壁對(duì)流體的流速影響范圍很小,而圖6(e)所示的極限粗糙度圓管流速分布圖,管壁附近流速不均勻分布區(qū)則要大得多。

    (3) 從圖中可以看出,靠近圓管壁的流體流速普遍小于圓管軸線中心附近的流體流速,在圓管壁設(shè)置滑動(dòng)邊界的條件下,斷面流速呈現(xiàn)U型分布,而非拋物線型分布,該結(jié)論和文獻(xiàn)[15]是一致的。

    (4) 從圖6中還可以看出,由于加壓差的作用,圓管頂部和底部出現(xiàn)了流速擾動(dòng)區(qū),一般情況下該擾動(dòng)區(qū)域相對(duì)較小,對(duì)結(jié)果分析影響較小。但在極限粗糙度情況下(如圖6(e)所示),管道端部影響使得管道中流體的流動(dòng)非常復(fù)雜。

    圖6壓力差P=2000Pa時(shí)不同粗糙度圓管XOZ平面的穩(wěn)定流速分布圖

    4結(jié)語(yǔ)

    本文通過(guò)顆粒流流體計(jì)算程序模擬尼古拉茲人工粗糙圓管實(shí)驗(yàn),說(shuō)明顆粒流流體計(jì)算模塊不僅可以實(shí)現(xiàn)對(duì)流體和顆粒的相互作用分析,而且對(duì)于粗糙圓管模擬可以較為真實(shí)的反應(yīng)流體在圓管中的流動(dòng)狀態(tài),對(duì)分析圓管內(nèi)流體流態(tài)有一定的借鑒意義和參考價(jià)值。

    參考文獻(xiàn):

    [1]陳仲頤,周景星,王洪瑾編著.土力學(xué)[M].北京:清華大學(xué)出版社,1994:36-46.

    [2]李大美,楊小亭主編.水力學(xué)[M].武漢:武漢大學(xué)出版社,2004:261-285.

    [3]Nikuradse,Johann.LawsofFlowinRoughPipes[M].Washington:NationalAdvisoryCommitteeforAeronautics, 1950:1-26.

    [4]黃琨.孔隙介質(zhì)滲流基本方程的探索[D].武漢:中國(guó)地質(zhì)大學(xué),2012:1-18.

    [5]萬(wàn)軍偉,黃琨,李銀泉,等.尼古拉茲圓管沿程水流阻力實(shí)驗(yàn)的新發(fā)現(xiàn)[J].地質(zhì)科技情報(bào),2012,31(5):157-161.

    [6]喬昌凱,劉興年,陳興偉,等.卵石透水河床滲流阻力研究[J].四川大學(xué)學(xué)報(bào):工程科學(xué)版,2013,45(S2):36-40.

    [7]HerwigH,GlossD,WenterodtT. 粗糙壁面通道內(nèi)的流動(dòng)研究[J].中國(guó)計(jì)量學(xué)院學(xué)報(bào),2008,19(4):296-303.

    [8]郝鵬飛,姚朝暉,何楓.粗糙微管道內(nèi)液體流動(dòng)特性的實(shí)驗(yàn)研究[J].物理學(xué)報(bào),2007,56(8):4728-4732.

    [9]周健,姚志雄,張剛.基于散體介質(zhì)理論的砂土管涌機(jī)制研究[J].巖石力學(xué)與工程學(xué)報(bào),2008,27(4):749-756.

    [10]劉洋,周健,付建新.飽和砂土流固耦合細(xì)觀數(shù)值模型及其在液化分析中的應(yīng)用[J].水利學(xué)報(bào),2009,40(2):250-256.

    [11]劉洋,李飛,柴小慶,等.滲流的PFC-CFD耦合細(xì)觀數(shù)值模擬[J].水文地質(zhì)工程地質(zhì),2008,35(2):66-70.

    [12]劉洋,吳順川,周健.砂沸現(xiàn)象的細(xì)觀數(shù)值模擬及影響因素分析[J].北京科技大學(xué)學(xué)報(bào),2008,30(5):473-477.

    [13]ErgunS.Fluidflowthroughpackedcolumns[J].JournalofMaterialsScience&ChemicalEngineering, 1952,48(2):89-94.

    [14]WenCY,Y.HoYu.Mechanicsoffluidization[J].Chem.engngProg.symp.ser, 1966,62:100-111.

    [15]WhiteFM.Fluidmechanics[M].NewYork:McGraw-Hill, 2011:1-49.

    Effects of Pipe Roughness on Flow Pattern Within Circular Pipe

    LI Huan,ZENG Yawu, XIA Lei, FANG Hailong

    (SchoolofCivilEngineering,WuhanUniversity,Wuhan,Hubei430072,China)

    Abstract:Based on the idea of Nikuradse’s experiment, this paper simulated the roughness of pipes by using Particle Flow Code (PFC) FLUID SCHEME, and analyzed the effects of different roughness of pipes to the flow pattern and the relationship between the mean velocity and roughness under different pressures. The results showed that flow velocity within pipe was influenced by pipe wall, and presented U-shaped in cross-section profile, rather than parabolic distribution; The flow pattern was affected by the pipe roughness, the larger relative roughness of the pipe, the greater range fluid was perturbed, namely laminar flow range become smaller. Under certain circumstances of pipe roughness, pipe flow increases with the increased pressure; while under the same pressure difference, pipe flow decreases with the increased roughness of pipe and the particles damping action is significant. While increasing roughness, the affected area within pipe also increased under different pressure, which make fluid flow patterns within pipe become more complex.

    Keywords:Nikuradse’s experiment; roughness pipe; PFC Fluid Scheme; relative roughness; flow pattern

    中圖分類(lèi)號(hào):TV134.1

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1672—1144(2016)02—0045—06

    作者簡(jiǎn)介:李歡(1990—),男,湖北黃梅人,碩士研究生,研究方向?yàn)殡x散元模擬及水土相互作用。 E-mail:llihhuan@whu.edu.cn

    收稿日期:2015-12-08修稿日期:2015-12-28

    DOI:10.3969/j.issn.1672-1144.2016.02.009

    亚洲国产成人一精品久久久| 亚洲精品乱码久久久v下载方式| 成人美女网站在线观看视频| 亚洲经典国产精华液单| 国内少妇人妻偷人精品xxx网站| 亚洲精品,欧美精品| 九九久久精品国产亚洲av麻豆| 精品一区在线观看国产| 亚洲国产色片| 国产色爽女视频免费观看| 伦精品一区二区三区| 精品熟女少妇av免费看| 亚州av有码| 校园人妻丝袜中文字幕| 久久久久性生活片| 午夜视频国产福利| 午夜福利网站1000一区二区三区| 日韩制服骚丝袜av| 女人十人毛片免费观看3o分钟| 简卡轻食公司| 丝袜喷水一区| 国产老妇女一区| 日本色播在线视频| 最近视频中文字幕2019在线8| 嫩草影院入口| 高清毛片免费看| 久久久久免费精品人妻一区二区| 日韩欧美 国产精品| 亚洲美女视频黄频| 亚洲精品久久午夜乱码| 亚洲成人一二三区av| 成人欧美大片| 国产女主播在线喷水免费视频网站 | 午夜福利在线在线| 久久久久久国产a免费观看| 亚洲综合精品二区| 国产在视频线在精品| 特级一级黄色大片| 亚洲自拍偷在线| 亚洲熟妇中文字幕五十中出| 国语对白做爰xxxⅹ性视频网站| 国产精品伦人一区二区| 日韩成人伦理影院| 国产av在哪里看| 午夜激情久久久久久久| 插阴视频在线观看视频| 日本三级黄在线观看| 国产成人freesex在线| 日日撸夜夜添| 搡老妇女老女人老熟妇| 午夜精品国产一区二区电影 | www.av在线官网国产| 免费高清在线观看视频在线观看| 国产成人精品婷婷| 亚洲最大成人手机在线| 精品国内亚洲2022精品成人| 国产女主播在线喷水免费视频网站 | 春色校园在线视频观看| 日日啪夜夜爽| 最近中文字幕2019免费版| 亚洲精品日韩av片在线观看| 一区二区三区四区激情视频| 晚上一个人看的免费电影| 看黄色毛片网站| 日本黄大片高清| 精品一区二区免费观看| 最后的刺客免费高清国语| 亚洲一区高清亚洲精品| 国产精品人妻久久久久久| 亚洲国产高清在线一区二区三| 国产爱豆传媒在线观看| 亚洲人与动物交配视频| 免费大片18禁| 岛国毛片在线播放| 国产免费福利视频在线观看| 国产精品伦人一区二区| 国产精品人妻久久久影院| 日韩制服骚丝袜av| 韩国av在线不卡| av一本久久久久| 亚洲激情五月婷婷啪啪| 亚洲国产日韩欧美精品在线观看| 国产精品久久久久久精品电影小说 | 国产精品嫩草影院av在线观看| 2021天堂中文幕一二区在线观| 日韩视频在线欧美| 熟女人妻精品中文字幕| 夫妻午夜视频| 街头女战士在线观看网站| 亚洲欧美清纯卡通| 黑人高潮一二区| 国产av在哪里看| 大香蕉久久网| 我的女老师完整版在线观看| 丝瓜视频免费看黄片| 亚洲欧美清纯卡通| 激情五月婷婷亚洲| 黄片无遮挡物在线观看| 欧美最新免费一区二区三区| 国产激情偷乱视频一区二区| 天堂俺去俺来也www色官网 | 丰满乱子伦码专区| 婷婷色麻豆天堂久久| 亚洲欧美日韩东京热| videos熟女内射| 一区二区三区乱码不卡18| 国产av不卡久久| 亚洲av中文av极速乱| 久久久久国产网址| 晚上一个人看的免费电影| 男插女下体视频免费在线播放| 国产精品国产三级国产av玫瑰| 日韩人妻高清精品专区| 尾随美女入室| 国产乱人偷精品视频| 一夜夜www| 99热网站在线观看| 欧美最新免费一区二区三区| 国产淫语在线视频| 久久久久久久国产电影| 最近手机中文字幕大全| 九色成人免费人妻av| av免费观看日本| 天堂中文最新版在线下载 | 乱人视频在线观看| 国产精品女同一区二区软件| 中文天堂在线官网| 国产成人精品福利久久| 国产淫片久久久久久久久| 街头女战士在线观看网站| 赤兔流量卡办理| 欧美成人a在线观看| 中文欧美无线码| 亚洲av在线观看美女高潮| av在线播放精品| 国产精品综合久久久久久久免费| 别揉我奶头 嗯啊视频| 亚洲欧洲日产国产| 一级毛片黄色毛片免费观看视频| 秋霞在线观看毛片| 最近中文字幕高清免费大全6| 国产精品一及| 国产成人91sexporn| 亚洲熟女精品中文字幕| 亚洲精品自拍成人| 免费av不卡在线播放| 日韩欧美一区视频在线观看 | 又黄又爽又刺激的免费视频.| 国产色爽女视频免费观看| 亚洲在线自拍视频| 国产欧美日韩精品一区二区| 久久久久久久久久成人| 成年版毛片免费区| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久久电影| 精品人妻偷拍中文字幕| 国产不卡一卡二| 日韩三级伦理在线观看| 中国国产av一级| 日韩中字成人| 夜夜看夜夜爽夜夜摸| 白带黄色成豆腐渣| 亚洲欧美精品自产自拍| 免费高清在线观看视频在线观看| 高清日韩中文字幕在线| av卡一久久| 麻豆精品久久久久久蜜桃| 午夜精品一区二区三区免费看| 欧美日韩综合久久久久久| 黄片wwwwww| 综合色av麻豆| 嫩草影院入口| av一本久久久久| 国产人妻一区二区三区在| 亚洲成人一二三区av| 久久久a久久爽久久v久久| 亚洲国产日韩欧美精品在线观看| 欧美三级亚洲精品| 国产成人91sexporn| 69人妻影院| 亚洲欧美一区二区三区国产| 国产单亲对白刺激| 最近最新中文字幕大全电影3| 亚洲,欧美,日韩| 神马国产精品三级电影在线观看| 国产亚洲午夜精品一区二区久久 | 天天躁日日操中文字幕| 亚洲国产色片| a级毛色黄片| 黄色欧美视频在线观看| 国产黄片美女视频| 欧美bdsm另类| 免费观看性生交大片5| 久久人人爽人人片av| 秋霞伦理黄片| 身体一侧抽搐| 亚洲18禁久久av| 男人爽女人下面视频在线观看| 亚洲图色成人| 赤兔流量卡办理| 少妇高潮的动态图| 国产精品福利在线免费观看| 亚洲熟妇中文字幕五十中出| 少妇熟女欧美另类| 一个人看视频在线观看www免费| 亚洲精品成人久久久久久| 在线观看免费高清a一片| 久久久久久久国产电影| 国产精品一区二区性色av| 精品久久国产蜜桃| 91精品国产九色| 国产精品一区二区在线观看99 | 国语对白做爰xxxⅹ性视频网站| 18禁动态无遮挡网站| 黄色一级大片看看| 国产成人免费观看mmmm| 国产69精品久久久久777片| 日韩,欧美,国产一区二区三区| 精品国产三级普通话版| 久久99热6这里只有精品| 一级毛片黄色毛片免费观看视频| 69av精品久久久久久| 日韩欧美国产在线观看| 亚洲激情五月婷婷啪啪| 午夜福利视频精品| 可以在线观看毛片的网站| 久久精品久久久久久噜噜老黄| 18禁在线无遮挡免费观看视频| 国产高清三级在线| 只有这里有精品99| 最后的刺客免费高清国语| 男女那种视频在线观看| 亚洲真实伦在线观看| 一夜夜www| 久久精品国产亚洲网站| 大话2 男鬼变身卡| 久久97久久精品| av黄色大香蕉| 男女啪啪激烈高潮av片| 欧美一区二区亚洲| 蜜桃亚洲精品一区二区三区| 日本免费在线观看一区| 欧美另类一区| 女的被弄到高潮叫床怎么办| 免费观看在线日韩| 日本免费a在线| 亚洲人成网站在线播| 肉色欧美久久久久久久蜜桃 | 毛片一级片免费看久久久久| 国产成人a区在线观看| 亚洲精品亚洲一区二区| 能在线免费观看的黄片| 成年人午夜在线观看视频 | 久久精品夜色国产| 色综合色国产| 五月伊人婷婷丁香| 国产69精品久久久久777片| 免费黄网站久久成人精品| 日韩,欧美,国产一区二区三区| 亚洲av.av天堂| 日韩av不卡免费在线播放| 18禁在线播放成人免费| 国产精品蜜桃在线观看| 国产一级毛片七仙女欲春2| 久久久久久久午夜电影| 午夜久久久久精精品| 日本-黄色视频高清免费观看| 又大又黄又爽视频免费| 亚洲婷婷狠狠爱综合网| 乱人视频在线观看| 高清在线视频一区二区三区| 国产亚洲5aaaaa淫片| 国产精品一区二区三区四区久久| 久久久久精品久久久久真实原创| 成人av在线播放网站| 特级一级黄色大片| 亚洲国产欧美在线一区| 国产精品国产三级国产专区5o| 中文字幕人妻熟人妻熟丝袜美| 日本黄大片高清| 深爱激情五月婷婷| 国产精品美女特级片免费视频播放器| 亚洲最大成人中文| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区乱码不卡18| 午夜视频国产福利| 九九久久精品国产亚洲av麻豆| 国产老妇伦熟女老妇高清| 色5月婷婷丁香| 亚洲国产av新网站| 能在线免费观看的黄片| 国产白丝娇喘喷水9色精品| 91av网一区二区| 国产黄a三级三级三级人| 韩国高清视频一区二区三区| 国产黄片美女视频| 亚洲av免费在线观看| 好男人在线观看高清免费视频| 80岁老熟妇乱子伦牲交| 亚洲av.av天堂| 在线 av 中文字幕| 熟妇人妻久久中文字幕3abv| 自拍偷自拍亚洲精品老妇| 精品久久久久久久末码| 久99久视频精品免费| 黄色一级大片看看| 国产av在哪里看| 午夜福利在线观看吧| 欧美性感艳星| 精品久久久噜噜| 亚洲欧美日韩无卡精品| 禁无遮挡网站| av在线亚洲专区| 蜜臀久久99精品久久宅男| 一个人免费在线观看电影| 久久精品国产鲁丝片午夜精品| 日本午夜av视频| 麻豆av噜噜一区二区三区| 老司机影院毛片| av一本久久久久| 亚洲一区高清亚洲精品| 身体一侧抽搐| 国产精品爽爽va在线观看网站| 亚洲图色成人| 免费av毛片视频| 国产美女午夜福利| 色网站视频免费| 男女国产视频网站| 狠狠精品人妻久久久久久综合| 免费av不卡在线播放| 91久久精品国产一区二区三区| 日本黄色片子视频| 亚洲av日韩在线播放| 26uuu在线亚洲综合色| 久久久久精品性色| 中文精品一卡2卡3卡4更新| 久久精品久久久久久久性| 中文资源天堂在线| 高清av免费在线| 又爽又黄无遮挡网站| 永久免费av网站大全| 一区二区三区四区激情视频| 久久久精品免费免费高清| 日本与韩国留学比较| 婷婷色综合www| 97超碰精品成人国产| 亚洲av男天堂| 嘟嘟电影网在线观看| 三级国产精品欧美在线观看| 亚洲国产精品成人综合色| 中文字幕av在线有码专区| 国产精品av视频在线免费观看| 国产精品嫩草影院av在线观看| 国产亚洲精品av在线| 日韩制服骚丝袜av| 欧美成人精品欧美一级黄| 国产免费一级a男人的天堂| 少妇猛男粗大的猛烈进出视频 | 女人久久www免费人成看片| 国产亚洲91精品色在线| 国产在视频线在精品| 亚洲激情五月婷婷啪啪| 老司机影院成人| 国产精品99久久久久久久久| 亚洲国产欧美人成| 麻豆av噜噜一区二区三区| 我的老师免费观看完整版| 国产国拍精品亚洲av在线观看| 欧美 日韩 精品 国产| 日韩欧美 国产精品| 深爱激情五月婷婷| 别揉我奶头 嗯啊视频| 2022亚洲国产成人精品| 两个人视频免费观看高清| 别揉我奶头 嗯啊视频| 2021天堂中文幕一二区在线观| 国产亚洲精品久久久com| 免费av不卡在线播放| 一个人观看的视频www高清免费观看| 成人一区二区视频在线观看| 亚洲自偷自拍三级| 亚洲av免费高清在线观看| 三级经典国产精品| 亚洲欧美清纯卡通| 欧美人与善性xxx| 成人特级av手机在线观看| 毛片女人毛片| 亚洲第一区二区三区不卡| 国产黄色免费在线视频| 能在线免费看毛片的网站| 九草在线视频观看| or卡值多少钱| 天天一区二区日本电影三级| 色综合站精品国产| 精品一区二区三卡| 亚洲欧洲日产国产| 国产淫语在线视频| 丰满人妻一区二区三区视频av| 国产色爽女视频免费观看| 中文字幕av成人在线电影| 亚洲欧美清纯卡通| 免费av不卡在线播放| 日本爱情动作片www.在线观看| 午夜免费男女啪啪视频观看| 五月伊人婷婷丁香| 日本午夜av视频| 亚洲精品国产av蜜桃| 2021天堂中文幕一二区在线观| 午夜精品在线福利| 在线天堂最新版资源| freevideosex欧美| 亚洲最大成人手机在线| 国产成人freesex在线| 亚洲内射少妇av| 国产色爽女视频免费观看| 亚洲国产精品专区欧美| 亚洲国产色片| 2021天堂中文幕一二区在线观| 久久久色成人| 亚洲成人av在线免费| 啦啦啦中文免费视频观看日本| 精品久久久久久久久久久久久| 国产三级在线视频| 日韩欧美精品免费久久| 成人亚洲精品一区在线观看 | 精品国产三级普通话版| 欧美极品一区二区三区四区| 亚洲精品成人久久久久久| 三级国产精品片| 国产高清不卡午夜福利| 久久久色成人| 欧美区成人在线视频| 国精品久久久久久国模美| 久久久久免费精品人妻一区二区| 国产精品av视频在线免费观看| 免费黄频网站在线观看国产| 国产永久视频网站| 99久国产av精品国产电影| 三级男女做爰猛烈吃奶摸视频| 大话2 男鬼变身卡| 一本久久精品| 99视频精品全部免费 在线| 国产一级毛片七仙女欲春2| 日韩 亚洲 欧美在线| 男的添女的下面高潮视频| 日本一二三区视频观看| av.在线天堂| 欧美激情在线99| 尾随美女入室| 久久综合国产亚洲精品| 十八禁国产超污无遮挡网站| 床上黄色一级片| 亚洲三级黄色毛片| 床上黄色一级片| 欧美潮喷喷水| av天堂中文字幕网| freevideosex欧美| 免费观看无遮挡的男女| 精品少妇黑人巨大在线播放| 国产精品女同一区二区软件| 中国美白少妇内射xxxbb| 激情 狠狠 欧美| 久久精品夜色国产| 蜜臀久久99精品久久宅男| 能在线免费看毛片的网站| 女性被躁到高潮视频| 看十八女毛片水多多多| 伦理电影免费视频| 精品视频人人做人人爽| 啦啦啦在线观看免费高清www| 欧美 日韩 精品 国产| 亚洲熟女精品中文字幕| 国产野战对白在线观看| av在线观看视频网站免费| 国产精品无大码| av网站在线播放免费| 久久精品久久精品一区二区三区| 欧美+日韩+精品| 黄片播放在线免费| 1024香蕉在线观看| 亚洲欧美成人综合另类久久久| 国产一级毛片在线| 欧美精品一区二区免费开放| 青春草国产在线视频| 亚洲精品自拍成人| 最近2019中文字幕mv第一页| 永久免费av网站大全| 2021少妇久久久久久久久久久| 三级国产精品片| 丝袜人妻中文字幕| 下体分泌物呈黄色| 亚洲av欧美aⅴ国产| 一区二区三区激情视频| 午夜福利网站1000一区二区三区| 水蜜桃什么品种好| 99精国产麻豆久久婷婷| 午夜福利影视在线免费观看| 国产av码专区亚洲av| 欧美中文综合在线视频| 精品人妻在线不人妻| 最近中文字幕2019免费版| 丰满少妇做爰视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产男女超爽视频在线观看| 亚洲av电影在线观看一区二区三区| 国产深夜福利视频在线观看| 哪个播放器可以免费观看大片| 婷婷成人精品国产| 婷婷色麻豆天堂久久| 青春草国产在线视频| 欧美日韩视频精品一区| 国产日韩欧美视频二区| 丝袜人妻中文字幕| 国产精品久久久久成人av| 黄频高清免费视频| 夜夜骑夜夜射夜夜干| 亚洲第一区二区三区不卡| 啦啦啦视频在线资源免费观看| 国产成人精品一,二区| av在线观看视频网站免费| 有码 亚洲区| 亚洲成人手机| 欧美精品高潮呻吟av久久| 久久人妻熟女aⅴ| 国产老妇伦熟女老妇高清| 国产无遮挡羞羞视频在线观看| 亚洲激情五月婷婷啪啪| 亚洲欧美日韩另类电影网站| 欧美精品人与动牲交sv欧美| 久久久久国产精品人妻一区二区| 视频在线观看一区二区三区| 97精品久久久久久久久久精品| 欧美精品一区二区免费开放| 精品第一国产精品| 美女午夜性视频免费| 午夜久久久在线观看| 国产色婷婷99| 成人影院久久| 欧美黄色片欧美黄色片| 国产 一区精品| 久久人人爽人人片av| 一级片'在线观看视频| 一区福利在线观看| 国产人伦9x9x在线观看 | 在线观看三级黄色| 久久久久久久亚洲中文字幕| 亚洲美女黄色视频免费看| 国产精品久久久久久精品电影小说| 国产精品国产三级专区第一集| 新久久久久国产一级毛片| 欧美人与性动交α欧美精品济南到 | 秋霞在线观看毛片| 精品一区二区免费观看| 男人舔女人的私密视频| 亚洲精品国产一区二区精华液| 九草在线视频观看| 国产1区2区3区精品| 91午夜精品亚洲一区二区三区| 成年av动漫网址| 久久久久久伊人网av| 18禁国产床啪视频网站| 亚洲精品国产av成人精品| 男女下面插进去视频免费观看| 日韩免费高清中文字幕av| 精品亚洲成a人片在线观看| 亚洲美女视频黄频| 亚洲欧美色中文字幕在线| 国产一区二区三区av在线| 丝袜在线中文字幕| av免费在线看不卡| 蜜桃国产av成人99| 可以免费在线观看a视频的电影网站 | 国产成人a∨麻豆精品| 精品国产乱码久久久久久男人| 在线 av 中文字幕| 婷婷色麻豆天堂久久| 女的被弄到高潮叫床怎么办| 亚洲伊人久久精品综合| 日本欧美视频一区| 久久影院123| 老汉色av国产亚洲站长工具| 综合色丁香网| 日韩熟女老妇一区二区性免费视频| 国语对白做爰xxxⅹ性视频网站| 亚洲av日韩在线播放| 在线免费观看不下载黄p国产| 五月伊人婷婷丁香| 在线观看人妻少妇| 欧美精品国产亚洲| 性高湖久久久久久久久免费观看| 王馨瑶露胸无遮挡在线观看| 国产有黄有色有爽视频| 婷婷色综合大香蕉| 国产成人欧美| 日韩伦理黄色片| 国产野战对白在线观看| 一级毛片电影观看| 精品酒店卫生间| 狂野欧美激情性bbbbbb| 99香蕉大伊视频| 亚洲av在线观看美女高潮| 2018国产大陆天天弄谢| 考比视频在线观看| 大片免费播放器 马上看| 亚洲欧美一区二区三区黑人 | 日本黄色日本黄色录像| 免费观看无遮挡的男女| 国产成人精品无人区| 国产精品久久久久久精品古装| 免费观看av网站的网址| 女的被弄到高潮叫床怎么办| 桃花免费在线播放| 日韩精品免费视频一区二区三区| 丰满迷人的少妇在线观看| 90打野战视频偷拍视频|