• <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

    亚洲av成人不卡在线观看播放网 | 国产日韩欧美视频二区| 自线自在国产av| 国产精品久久久人人做人人爽| 天天添夜夜摸| 中文欧美无线码| 日本猛色少妇xxxxx猛交久久| 国产免费视频播放在线视频| 亚洲专区中文字幕在线 | 国精品久久久久久国模美| 美女国产高潮福利片在线看| 免费黄频网站在线观看国产| 美女高潮到喷水免费观看| 最新的欧美精品一区二区| 国产成人免费观看mmmm| 国产成人a∨麻豆精品| 国产亚洲一区二区精品| 久久天堂一区二区三区四区| 美女大奶头黄色视频| 一本—道久久a久久精品蜜桃钙片| 亚洲 欧美一区二区三区| a级毛片黄视频| 97在线人人人人妻| 90打野战视频偷拍视频| 国产熟女午夜一区二区三区| 久久久久久久精品精品| 精品国产超薄肉色丝袜足j| 热re99久久精品国产66热6| 欧美av亚洲av综合av国产av | 亚洲一区中文字幕在线| 在线观看人妻少妇| 99国产综合亚洲精品| 美女主播在线视频| 91老司机精品| 在线免费观看不下载黄p国产| 国产精品久久久av美女十八| 精品国产一区二区三区久久久樱花| 久久久久精品久久久久真实原创| 人人妻人人添人人爽欧美一区卜| 丰满饥渴人妻一区二区三| 欧美黑人精品巨大| 男女免费视频国产| 日日撸夜夜添| 欧美日韩国产mv在线观看视频| 国产爽快片一区二区三区| 一本一本久久a久久精品综合妖精| 美女大奶头黄色视频| 色综合欧美亚洲国产小说| 亚洲欧美一区二区三区久久| 亚洲图色成人| 岛国毛片在线播放| 中文字幕人妻熟女乱码| 国产成人免费无遮挡视频| 免费看av在线观看网站| 免费在线观看完整版高清| 色精品久久人妻99蜜桃| 黑人猛操日本美女一级片| 99久国产av精品国产电影| av在线app专区| 五月开心婷婷网| 国产 一区精品| 青青草视频在线视频观看| 一级黄片播放器| 在线观看免费视频网站a站| 免费女性裸体啪啪无遮挡网站| 国产麻豆69| 中文字幕av电影在线播放| 亚洲精品国产区一区二| 国产成人一区二区在线| 亚洲一卡2卡3卡4卡5卡精品中文| 啦啦啦中文免费视频观看日本| 青春草国产在线视频| 久久毛片免费看一区二区三区| 永久免费av网站大全| 国产一区亚洲一区在线观看| 2021少妇久久久久久久久久久| 丝袜喷水一区| 亚洲成人av在线免费| 国产精品秋霞免费鲁丝片| 丝袜美足系列| 日本91视频免费播放| 大陆偷拍与自拍| 国产成人91sexporn| 精品少妇内射三级| 极品少妇高潮喷水抽搐| 久久人妻熟女aⅴ| 成人国产av品久久久| 亚洲av日韩精品久久久久久密 | 少妇的丰满在线观看| 一级毛片我不卡| 97在线人人人人妻| 欧美日韩亚洲综合一区二区三区_| 精品久久久久久电影网| 亚洲欧美一区二区三区久久| 国产免费又黄又爽又色| 纵有疾风起免费观看全集完整版| 免费久久久久久久精品成人欧美视频| 777久久人妻少妇嫩草av网站| 久久精品熟女亚洲av麻豆精品| 精品福利永久在线观看| 成人漫画全彩无遮挡| 国产乱来视频区| 国产一卡二卡三卡精品 | 超碰成人久久| 考比视频在线观看| 欧美日韩成人在线一区二区| 日日爽夜夜爽网站| 18禁观看日本| 免费观看人在逋| 国产精品亚洲av一区麻豆 | 国产黄频视频在线观看| 精品少妇久久久久久888优播| 久热爱精品视频在线9| 亚洲精品中文字幕在线视频| 97在线人人人人妻| 亚洲av综合色区一区| 韩国精品一区二区三区| 秋霞伦理黄片| www.av在线官网国产| 久久99一区二区三区| 久久99一区二区三区| 久久99精品国语久久久| 日韩大片免费观看网站| 亚洲av男天堂| 老鸭窝网址在线观看| 婷婷色麻豆天堂久久| 欧美日韩成人在线一区二区| 国产国语露脸激情在线看| 18禁裸乳无遮挡动漫免费视频| 男女国产视频网站| 日本一区二区免费在线视频| 国产成人91sexporn| 欧美另类一区| 精品国产国语对白av| 一级a爱视频在线免费观看| 老汉色∧v一级毛片| 在线观看人妻少妇| 不卡视频在线观看欧美| 色网站视频免费| 只有这里有精品99| 久久久久国产一级毛片高清牌| 久久人人爽av亚洲精品天堂| 一区二区三区乱码不卡18| 最近的中文字幕免费完整| 黄色 视频免费看| 国产精品久久久av美女十八| 老鸭窝网址在线观看| 亚洲专区中文字幕在线 | av卡一久久| 亚洲av成人不卡在线观看播放网 | 日韩熟女老妇一区二区性免费视频| 免费在线观看黄色视频的| 欧美日韩国产mv在线观看视频| 日本欧美视频一区| 久久久亚洲精品成人影院| 亚洲国产精品一区三区| 日日啪夜夜爽| 一级,二级,三级黄色视频| 国产激情久久老熟女| 久久久久精品性色| 久久鲁丝午夜福利片| 国产成人免费观看mmmm| 亚洲熟女毛片儿| 久久久久精品人妻al黑| 操出白浆在线播放| 日韩欧美精品免费久久| 国产淫语在线视频| 只有这里有精品99| 嫩草影院入口| 国产成人午夜福利电影在线观看| 亚洲国产av影院在线观看| 亚洲欧美清纯卡通| 国产伦理片在线播放av一区| 成人国产av品久久久| 中文天堂在线官网| 丰满迷人的少妇在线观看| 男人舔女人的私密视频| 欧美精品一区二区大全| 黄色毛片三级朝国网站| av视频免费观看在线观看| 成人手机av| 国产免费视频播放在线视频| 久久久欧美国产精品| 亚洲天堂av无毛| 精品卡一卡二卡四卡免费| 亚洲欧美色中文字幕在线| 又大又爽又粗| 午夜91福利影院| 久久狼人影院| 99精品久久久久人妻精品| 亚洲欧美清纯卡通| 亚洲在久久综合| 国产一区二区在线观看av| 久久久久精品性色| 国产男女超爽视频在线观看| 成人三级做爰电影| 飞空精品影院首页| av网站免费在线观看视频| 黄片播放在线免费| 精品亚洲乱码少妇综合久久| 男的添女的下面高潮视频| 人人澡人人妻人| 亚洲伊人色综图| 男女下面插进去视频免费观看| 青青草视频在线视频观看| 男女午夜视频在线观看| 久久精品亚洲av国产电影网| 亚洲综合精品二区| 午夜福利,免费看| 久久久精品国产亚洲av高清涩受| 国产日韩欧美视频二区| 国产高清不卡午夜福利| av天堂久久9| 老熟女久久久| 亚洲人成电影观看| 人妻人人澡人人爽人人| 午夜福利,免费看| 操美女的视频在线观看| 久久精品久久久久久久性| 午夜福利免费观看在线| 在线观看国产h片| www.自偷自拍.com| av在线观看视频网站免费| 午夜福利视频精品| a级片在线免费高清观看视频| 少妇人妻 视频| 国产精品99久久99久久久不卡 | a级毛片黄视频| 久久久久精品人妻al黑| 色综合欧美亚洲国产小说| 老汉色∧v一级毛片| 国产成人精品久久二区二区91 | 18禁国产床啪视频网站| 亚洲免费av在线视频| 在线免费观看不下载黄p国产| 亚洲成色77777| 国产成人啪精品午夜网站| 欧美亚洲日本最大视频资源| 一级片免费观看大全| 高清不卡的av网站| 午夜福利一区二区在线看| 久久这里只有精品19| 老司机在亚洲福利影院| 午夜福利在线免费观看网站| 亚洲国产精品一区三区| 久久狼人影院| 亚洲一级一片aⅴ在线观看| 秋霞在线观看毛片| 国产成人精品无人区| 日日摸夜夜添夜夜爱| 日韩欧美一区视频在线观看| 七月丁香在线播放| 亚洲国产欧美日韩在线播放| 91国产中文字幕| 国产在视频线精品| 亚洲精品美女久久av网站| 色综合欧美亚洲国产小说| 18禁动态无遮挡网站| 伊人亚洲综合成人网| 啦啦啦在线免费观看视频4| 精品一区二区三卡| 国产成人91sexporn| 丰满饥渴人妻一区二区三| 男女之事视频高清在线观看 | 久久精品亚洲熟妇少妇任你| 黄色视频不卡| 欧美另类一区| 精品国产一区二区三区四区第35| 亚洲国产av影院在线观看| 黄网站色视频无遮挡免费观看| 你懂的网址亚洲精品在线观看| 国产精品女同一区二区软件| 国产精品久久久人人做人人爽| 中国三级夫妇交换| 狠狠精品人妻久久久久久综合| 欧美精品人与动牲交sv欧美| av不卡在线播放| 国产探花极品一区二区| 激情五月婷婷亚洲| 午夜影院在线不卡| 桃花免费在线播放| 丝瓜视频免费看黄片| 午夜日本视频在线| 日本欧美国产在线视频| 性高湖久久久久久久久免费观看| 国产精品一区二区在线不卡| 一区二区日韩欧美中文字幕| 成人国产av品久久久| 少妇被粗大猛烈的视频| 建设人人有责人人尽责人人享有的| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久久精品精品| 亚洲,欧美精品.| 自拍欧美九色日韩亚洲蝌蚪91| 飞空精品影院首页| 91aial.com中文字幕在线观看| 日日啪夜夜爽| 日本av免费视频播放| a 毛片基地| 99国产综合亚洲精品| 最近中文字幕高清免费大全6| 日韩视频在线欧美| 少妇被粗大猛烈的视频| 亚洲精品国产区一区二| 91aial.com中文字幕在线观看| 老司机在亚洲福利影院| 99国产综合亚洲精品| 亚洲精品av麻豆狂野| 久久毛片免费看一区二区三区| 国产一区亚洲一区在线观看| 涩涩av久久男人的天堂| 国产精品秋霞免费鲁丝片| 亚洲精品中文字幕在线视频| 欧美在线一区亚洲| 一本大道久久a久久精品| 亚洲中文av在线| 国产高清国产精品国产三级| 少妇被粗大的猛进出69影院| av网站在线播放免费| 极品少妇高潮喷水抽搐| 欧美日韩福利视频一区二区| 国产无遮挡羞羞视频在线观看| 精品福利永久在线观看| 高清欧美精品videossex| 亚洲在久久综合| 国产亚洲最大av| 搡老岳熟女国产| 日韩不卡一区二区三区视频在线| 天天躁夜夜躁狠狠躁躁| av在线观看视频网站免费| 一二三四在线观看免费中文在| 欧美日韩亚洲综合一区二区三区_| 中文字幕av电影在线播放| 国产成人一区二区在线| 亚洲情色 制服丝袜| 精品视频人人做人人爽| 777米奇影视久久| 亚洲伊人色综图| 亚洲欧洲国产日韩| 亚洲美女搞黄在线观看| 中文字幕人妻熟女乱码| 国产一级毛片在线| 免费观看a级毛片全部| 亚洲欧美日韩另类电影网站| 亚洲精品自拍成人| 免费黄色在线免费观看| 婷婷色麻豆天堂久久| 黄色视频在线播放观看不卡| 国产福利在线免费观看视频| 男人爽女人下面视频在线观看| 你懂的网址亚洲精品在线观看| 色播在线永久视频| 伦理电影大哥的女人| 啦啦啦中文免费视频观看日本| 亚洲男人天堂网一区| 国产精品欧美亚洲77777| 精品一区二区三卡| 欧美国产精品一级二级三级| 精品人妻在线不人妻| 精品久久蜜臀av无| 波多野结衣av一区二区av| 少妇人妻久久综合中文| 亚洲av成人不卡在线观看播放网 | 亚洲精品,欧美精品| 97人妻天天添夜夜摸| 国产成人精品久久二区二区91 | 夫妻性生交免费视频一级片| 久久精品国产综合久久久| 亚洲av日韩精品久久久久久密 | 五月天丁香电影| 精品国产露脸久久av麻豆| 国产一区二区 视频在线| 久久久精品免费免费高清| 美女扒开内裤让男人捅视频| 国产成人系列免费观看| 亚洲成国产人片在线观看| 男女午夜视频在线观看| 最新在线观看一区二区三区 | 考比视频在线观看| 伦理电影免费视频| 欧美激情 高清一区二区三区| 夫妻午夜视频| 久热这里只有精品99| 国产精品免费大片| 国产午夜精品一二区理论片| 成人18禁高潮啪啪吃奶动态图| 亚洲视频免费观看视频| 黄片播放在线免费| 国产在线视频一区二区| 大片免费播放器 马上看| 九草在线视频观看| 大片免费播放器 马上看| 国产免费现黄频在线看| 日韩欧美精品免费久久| 国产熟女欧美一区二区| 一边摸一边抽搐一进一出视频| 精品免费久久久久久久清纯 | 男人添女人高潮全过程视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品aⅴ在线观看| 午夜免费鲁丝| 日韩,欧美,国产一区二区三区| 91老司机精品| 国产在线免费精品| 亚洲综合精品二区| 国产精品av久久久久免费| 王馨瑶露胸无遮挡在线观看| 欧美激情极品国产一区二区三区| 国产1区2区3区精品| 九草在线视频观看| 99久久99久久久精品蜜桃| 久久精品人人爽人人爽视色| 99九九在线精品视频| 国产精品av久久久久免费| 亚洲欧美清纯卡通| 人成视频在线观看免费观看| 少妇 在线观看| 日韩一本色道免费dvd| 麻豆精品久久久久久蜜桃| 91国产中文字幕| 欧美亚洲日本最大视频资源| 人人妻人人澡人人看| 久久人人97超碰香蕉20202| 免费在线观看黄色视频的| 日韩制服丝袜自拍偷拍| 色婷婷av一区二区三区视频| 精品国产一区二区三区四区第35| 男女国产视频网站| 91精品三级在线观看| 高清不卡的av网站| 人体艺术视频欧美日本| 色婷婷久久久亚洲欧美| 亚洲精品视频女| 国产精品一二三区在线看| 男人操女人黄网站| 老熟女久久久| 久久久久国产一级毛片高清牌| 人妻 亚洲 视频| 成人18禁高潮啪啪吃奶动态图| 欧美黑人欧美精品刺激| 亚洲精品av麻豆狂野| 欧美最新免费一区二区三区| 美女国产高潮福利片在线看| 亚洲美女黄色视频免费看| 国产精品蜜桃在线观看| 免费不卡黄色视频| 久久狼人影院| 人人澡人人妻人| 久久久国产精品麻豆| 亚洲国产最新在线播放| 国产成人精品久久久久久| 亚洲欧美中文字幕日韩二区| 一本色道久久久久久精品综合| 婷婷色av中文字幕| 国产极品粉嫩免费观看在线| 无遮挡黄片免费观看| 熟女少妇亚洲综合色aaa.| 国产午夜精品一二区理论片| 91老司机精品| 波野结衣二区三区在线| 国产精品久久久久久精品电影小说| 亚洲精品日韩在线中文字幕| 国产精品蜜桃在线观看| 国产97色在线日韩免费| 久久久国产欧美日韩av| 亚洲国产中文字幕在线视频| 国产精品国产av在线观看| 精品一区二区三区av网在线观看 | 黑人欧美特级aaaaaa片| 久久精品人人爽人人爽视色| 不卡av一区二区三区| 久久免费观看电影| 黄色 视频免费看| 亚洲国产av影院在线观看| 精品国产一区二区三区四区第35| 国产日韩一区二区三区精品不卡| 国产精品.久久久| 日韩欧美精品免费久久| 电影成人av| 亚洲成人一二三区av| 午夜老司机福利片| 精品一品国产午夜福利视频| 天堂8中文在线网| 国产一区二区激情短视频 | 嫩草影院入口| 永久免费av网站大全| 一区二区三区激情视频| 波多野结衣一区麻豆| 高清在线视频一区二区三区| 亚洲国产中文字幕在线视频| 亚洲精品美女久久av网站| 少妇猛男粗大的猛烈进出视频| 亚洲av日韩在线播放| 国产精品免费视频内射| 一区二区三区精品91| 热re99久久精品国产66热6| 日本vs欧美在线观看视频| 久热这里只有精品99| 一区二区三区四区激情视频| 久久99一区二区三区| 妹子高潮喷水视频| 亚洲av男天堂| 亚洲成人国产一区在线观看 | 乱人伦中国视频| 大片电影免费在线观看免费| 亚洲成av片中文字幕在线观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品视频女| 九九爱精品视频在线观看| 2018国产大陆天天弄谢| 久久精品国产亚洲av涩爱| 欧美乱码精品一区二区三区| 日韩大码丰满熟妇| av在线app专区| 中文字幕最新亚洲高清| 国产av精品麻豆| 一边摸一边做爽爽视频免费| 国产成人91sexporn| 久久毛片免费看一区二区三区| av卡一久久| a级毛片在线看网站| 久久久久久人妻| 男女边吃奶边做爰视频| 色94色欧美一区二区| 亚洲国产精品成人久久小说| 欧美日韩福利视频一区二区| 亚洲精品久久成人aⅴ小说| 亚洲,一卡二卡三卡| 国产成人精品无人区| 晚上一个人看的免费电影| 亚洲综合色网址| 亚洲国产日韩一区二区| 老司机亚洲免费影院| 街头女战士在线观看网站| 无遮挡黄片免费观看| 久久人人爽av亚洲精品天堂| 成年美女黄网站色视频大全免费| 亚洲精品中文字幕在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久精品人妻al黑| a级毛片黄视频| 操出白浆在线播放| 欧美黑人精品巨大| 大片免费播放器 马上看| 1024视频免费在线观看| 精品人妻一区二区三区麻豆| 两个人看的免费小视频| 国产熟女欧美一区二区| 婷婷色av中文字幕| 久久久国产一区二区| 捣出白浆h1v1| 欧美日韩精品网址| 国产一区二区在线观看av| 18禁国产床啪视频网站| 欧美人与性动交α欧美软件| 日本色播在线视频| 国产色婷婷99| av在线播放精品| 母亲3免费完整高清在线观看| 自线自在国产av| 国产人伦9x9x在线观看| 免费观看人在逋| 国产一级毛片在线| 免费久久久久久久精品成人欧美视频| 亚洲av日韩精品久久久久久密 | 国产在线一区二区三区精| 制服人妻中文乱码| 丝袜美腿诱惑在线| 999精品在线视频| 国产男人的电影天堂91| 成人免费观看视频高清| 狂野欧美激情性bbbbbb| 日韩制服骚丝袜av| 亚洲精品一二三| av电影中文网址| 国产黄频视频在线观看| 日韩视频在线欧美| 精品少妇内射三级| 亚洲精品,欧美精品| 久久国产精品男人的天堂亚洲| xxxhd国产人妻xxx| 日韩人妻精品一区2区三区| 日日啪夜夜爽| 精品少妇久久久久久888优播| 国产精品久久久久久人妻精品电影 | 日韩一本色道免费dvd| 女人高潮潮喷娇喘18禁视频| 人人妻人人添人人爽欧美一区卜| 国产精品久久久久久人妻精品电影 | 午夜日本视频在线| 波多野结衣av一区二区av| 我的亚洲天堂| 久久人妻熟女aⅴ| 飞空精品影院首页| 丰满乱子伦码专区| 国产一区二区三区综合在线观看| 一本久久精品| 中文欧美无线码| 三上悠亚av全集在线观看| 国产探花极品一区二区| 亚洲精品aⅴ在线观看| 高清欧美精品videossex| 亚洲国产精品999| 91国产中文字幕| 国产黄色视频一区二区在线观看| 黄片无遮挡物在线观看| 天堂俺去俺来也www色官网| 中文字幕人妻丝袜制服| 欧美精品av麻豆av| 精品国产一区二区三区四区第35| 男女边吃奶边做爰视频| 深夜精品福利| 男女边吃奶边做爰视频| 不卡视频在线观看欧美| 亚洲精品国产av成人精品|