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

    平整輸沙床面水流阻力研究

    2023-07-20 09:27:34賀承明張根廣許曉陽(yáng)
    中國(guó)農(nóng)村水利水電 2023年7期
    關(guān)鍵詞:粗化輸沙床面

    賀承明,張根廣,許曉陽(yáng)

    (西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院,陜西 楊凌 712100)

    0 引 言

    水流阻力是河流運(yùn)動(dòng)力學(xué)領(lǐng)域重點(diǎn)研究的問題之一,對(duì)于研究水流結(jié)構(gòu)和泥沙運(yùn)動(dòng)問題具有重要意義。影響水流阻力的因素較多,主要包括泥沙級(jí)配、水流條件、床面形態(tài)及邊界條件等。通常學(xué)者將水流阻力劃分為沙粒阻力、沙波阻力和邊壁阻力[1],并建立了相應(yīng)的阻力計(jì)算公式。對(duì)于明渠水流,床面阻力是水流阻力的主要組成部分,包括沙粒阻力和沙波阻力。當(dāng)床面平整或是沙波發(fā)育不明顯時(shí),沙波阻力可忽略不計(jì)。

    在研究床面阻力時(shí),許多學(xué)者常常采用床面粗糙高度Ks來表示床面粗糙程度。例如,尼古拉茲在研究管道壁面阻力時(shí),通過在管道內(nèi)壁黏附均勻沙粒得到了圓管中水流阻力系數(shù)與雷諾數(shù)的關(guān)系[1];Engelund[2]取2.5 倍的均勻沙粒徑作為床面粗糙高度,即Ks=2.5d;Van Rijn[3]取3倍的d90作為非均勻床面的粗糙高度。趙連白[4]認(rèn)為,對(duì)于非均勻沙而言,粗細(xì)泥沙顆粒對(duì)阻力的貢獻(xiàn)不同,粗顆粒的貢獻(xiàn)大于細(xì)顆粒。當(dāng)d50≤0.2 mm時(shí),取Ks=0.5d50;當(dāng)0.2 mm<d50<6 mm 時(shí),取Ks=d50;當(dāng)d50≥6 mm 時(shí),取Ks=2d50。喻國(guó)良[5]認(rèn)為床面粗糙高度與泥沙不均勻性有關(guān),并分析得到Ks=σg2d50,σg為床沙級(jí)配的非均勻性。還有一些學(xué)者研究了沙粒阻力與水流條件的關(guān)系。例如,黃才安[6]利用泥沙起動(dòng)的Shields曲線與明渠阻力特性之間的相似性,結(jié)合實(shí)測(cè)泥沙起動(dòng)資料得到了沙粒阻力系數(shù)只與床面坡度有關(guān)的結(jié)論。何文社[7]通過理論分析,提出了與泥沙起動(dòng)水力條件及泥沙重度有關(guān)的泥沙等效粒徑的概念及其計(jì)算方法,并在此基礎(chǔ)上得到了粗化穩(wěn)定床面的粗糙高度計(jì)算公式。陳學(xué)彪[8]基于阻力特性,綜合考慮水沙條件,建立了以泥沙中值粒徑、非均勻系數(shù)、水力半徑和床面坡度為變量的床面粗糙度計(jì)算公式。

    在可動(dòng)床面上,水流運(yùn)動(dòng)往往會(huì)引起泥沙輸移,因此將消耗水流的能量。推移質(zhì)輸沙率越大,水流搬運(yùn)泥沙所消耗的能量越多,水流受到的阻力也越大。Song[9]、Gao[10]通過定床清水和定床輸沙試驗(yàn)發(fā)現(xiàn),在推移質(zhì)輸沙率較大時(shí),輸沙阻力系數(shù)與床面總阻力系數(shù)的比值超過了0.3,可見推移質(zhì)輸移引起的水流阻力增大比較顯著,不可忽略。bagnold[11,12]基于能量平衡原理在推導(dǎo)推移質(zhì)輸沙率時(shí)認(rèn)為,泥沙顆粒在運(yùn)動(dòng)過程中,因其重力作用及相互碰撞消耗一定能量,將使泥沙保持下沉,因此需要水流提供一定作用力維持其運(yùn)動(dòng);并通過旋轉(zhuǎn)圓筒的顆粒碰撞剪切試驗(yàn),得到了水流搬運(yùn)推移質(zhì)泥沙所需做功的計(jì)算公式。黃才安[13]從挾沙水流的能量平衡原理出發(fā),考慮泥沙運(yùn)動(dòng)對(duì)水流阻力的影響,結(jié)合bagnold 水流功率理論,給出了動(dòng)平整床面挾沙水流阻力的計(jì)算公式。

    綜上分析,平整床面水流阻力與床面泥沙條件、水流條件及輸沙情況等因素有關(guān),前人的阻力計(jì)算公式較少考慮床面結(jié)構(gòu)及輸沙阻力,且輸沙阻力計(jì)算公式大多是基于實(shí)驗(yàn)資料的經(jīng)驗(yàn)公式或未對(duì)水流阻力進(jìn)行劃分。因此,本文將動(dòng)床阻力進(jìn)行分割,擬討論床面粗化過程中床面粗化對(duì)床面沙粒阻力的影響,結(jié)合水流條件分析計(jì)算床面沙粒阻力,并通過分析滾動(dòng)顆粒的運(yùn)動(dòng)特征,探討推移質(zhì)能量消耗規(guī)律,最后建立平整輸沙床面的水流阻力公式。

    1 阻力系數(shù)

    床面阻力常常通過床面阻力系數(shù)或床面粗糙高度來反映[1],其與水流結(jié)構(gòu)的關(guān)系可表示為:

    式中:f為床面阻力系數(shù);U為斷面平均流速,m∕s;u*為摩阻流速,m∕s;R為水力半徑,m;χ為流態(tài)校正系數(shù),在床面粗糙時(shí),取χ= 1;Ks為床面粗糙高度,m。

    在平整或沙波發(fā)育不明顯河床上,床面阻力可以劃分為沙粒阻力與輸沙阻力,沙粒阻力為床面上的泥沙顆粒對(duì)水流產(chǎn)生的表面阻力;輸沙阻力為水流為維持推移質(zhì)運(yùn)動(dòng)所消耗能量對(duì)應(yīng)的阻力,即:

    式中:f0為沙粒阻力系數(shù);fs為輸沙阻力系數(shù)。

    1.1 沙粒阻力系數(shù)

    由式(1)和式(2)可得床面沙粒阻力系數(shù)與床面沙粒阻力對(duì)應(yīng)的床面粗糙高度的關(guān)系為:

    式中:R'為床面沙粒阻力對(duì)應(yīng)的水力半徑,m;Ks0為床面沙粒阻力對(duì)應(yīng)的床面粗糙高度,m。

    何文社[7]認(rèn)為非均勻沙床面粗化穩(wěn)定后,其粗糙高度可由起動(dòng)粒徑dc計(jì)算:

    式中:m為粗糙高度系數(shù),當(dāng)床面粗化穩(wěn)定時(shí),;γs為泥沙顆粒的容重,N∕m3;γ為水的容重,N∕m3;ρ為水的密度,kg∕m3;u'*為床面沙粒阻力對(duì)應(yīng)的摩阻流速,;g為重力加速度,取9.8 m∕s2;J為床面坡度;θc為臨界Shields數(shù)。

    黃才安[6]利用泥沙起動(dòng)的Shields 曲線與明渠阻力特性之間的相似性,得到了均勻沙床面沙粒阻力系數(shù)計(jì)算公式,為:

    聯(lián)立式(4)和式(7)可得均勻沙床面沙粒阻力對(duì)應(yīng)的粗糙高度為:

    式中:當(dāng)過水?dāng)嗝鎸捝畋容^大時(shí),可取R'=h'。

    由式(5)與式(8)可見,床面粗糙高度可由床面上泥沙的起動(dòng)粒徑dc表示,起動(dòng)粒徑反映了水流作用在床面上的剪切力。兩式不同之處僅為系數(shù)的差異,究其原因是兩式適用的床面條件不同,公式(5)適用于粗化完成后的非均勻沙床面,公式(8)適用于均勻沙床面。

    在非均勻沙床面粗化完成前,床面粗化層伴隨著泥沙輸移而逐漸發(fā)展變化;在完成粗化后,床面粗細(xì)泥沙顆粒通過有規(guī)律的排列形成尺寸較大的局部穩(wěn)定結(jié)構(gòu)或顆粒簇[14],增大了床面的粗糙度。而在均勻沙床面上,床面泥沙均勻,不能形成粗化結(jié)構(gòu)。因此可見,在非均勻沙床面粗化過程中,床面泥沙粗化層結(jié)構(gòu)伴隨著泥沙輸移逐漸形成,相應(yīng)床面的粗糙高度是變化的,因此假定床面粗化過程中的床面粗糙高度可表示為:

    式中:c為床面粗化程度系數(shù);c與m的區(qū)別是,c為變量而非常數(shù);其他符號(hào)意義同前。

    為了確定床面粗糙高度系數(shù)c,作者對(duì)床面粗化過程進(jìn)行分析認(rèn)為,天然河流的床面往往是由非均勻沙組成,非均勻沙床面的輸沙過程也是床面的粗化過程,不同的水流條件塑造的床面粗化層也各不相同,其具體過程可簡(jiǎn)述為,在水流條件變化之前,河道流量較小,床面表層已完成粗化,床面泥沙基本上處于個(gè)別起動(dòng)狀態(tài)或靜止?fàn)顟B(tài),河流輸沙率接近于0;當(dāng)河道流量增大,床面切應(yīng)力也隨之增大,床面原有的粗化層被破壞,床面上可動(dòng)泥沙逐漸增多,當(dāng)粗化層被完全破壞時(shí),床面泥沙近似松散排列,輸沙率達(dá)到最大,可以認(rèn)為此時(shí)的河流輸沙率處于平衡狀態(tài),其單寬推移質(zhì)輸沙率為g'b0;隨著時(shí)間的推移,床面逐漸粗化,床面上聚集體結(jié)構(gòu)增多,并對(duì)細(xì)顆粒泥沙產(chǎn)生隱蔽作用[15],床面上可動(dòng)泥沙減少,床面單寬推移質(zhì)輸沙率g'b隨之減少;當(dāng)床面條件與來流量完全匹配時(shí),床面粗化完成,床面上基本沒有泥沙起動(dòng)輸移,輸沙率接近于0。由此可見,床面粗化程度與輸沙率高度匹配,因此可由輸沙率相對(duì)變化過程來反映床面粗化程度,則床面粗化程度系數(shù)c可表示為:

    式中:g'b0為以水下重量計(jì)的單寬推移質(zhì)平衡輸沙率,N∕(m·s);g'b為以水下重量計(jì)的單寬推移質(zhì)輸沙率,N∕(m·s)。

    g'b0可采用Meyer-Peter和Muller公式[1]計(jì)算

    式中:θ為Shields 數(shù);d50為床沙的中值粒徑,m;其他符號(hào)意義同前。

    有關(guān)泥沙起動(dòng)條件θc的取值。何文社[16]取顆粒完全暴露時(shí)的θc為0.02,適合于床面上沒有泥沙起動(dòng)的臨界情況。孟震[17]在研究推移質(zhì)輸沙率時(shí)建議,泥沙顆粒處于運(yùn)動(dòng)狀態(tài),泥沙起動(dòng)條件θc應(yīng)取對(duì)應(yīng)泥沙“個(gè)別動(dòng)”起動(dòng)條件的值,即θc=0.031。此外,其他一些學(xué)者采用的“個(gè)別”起動(dòng)條件θc基本都在0.03 左右,如樂培九[18]取0.03、周雙[19]取0.032 等。因此,本文選取泥沙“個(gè)別”起動(dòng)條件θc= 0.031。

    聯(lián)立式(4)、式(9)及式(10)可得床面沙粒阻力系數(shù)f0為:

    1.2 輸沙阻力系數(shù)

    推移質(zhì)泥沙在運(yùn)動(dòng)過程中會(huì)受到水流的剪切作用,使水流的部分能量轉(zhuǎn)化為推移質(zhì)顆粒的動(dòng)能,推移質(zhì)輸沙率越大,消耗的水流能量越多。推移質(zhì)受到的水流作用力與推移質(zhì)泥沙顆粒在水流中的運(yùn)動(dòng)特征有關(guān),而根據(jù)推移質(zhì)運(yùn)動(dòng)形式的不同,推移質(zhì)可以分為滾動(dòng)推移質(zhì)和躍移推移質(zhì)[1],其運(yùn)動(dòng)軌跡如圖1所示。

    圖1 滾動(dòng)和躍移推移質(zhì)運(yùn)動(dòng)過程Fig.1 The movement of rolling and saltating bedload

    由圖1可見,滾動(dòng)泥沙顆粒始終在床面附近運(yùn)動(dòng),軌跡幾乎為一條與床面平行的直線;躍移顆粒的運(yùn)動(dòng)軌跡為曲線,在垂直于水流方向,顆粒的運(yùn)動(dòng)范圍相對(duì)較大,躍高可達(dá)(4~5)d[20]。推移質(zhì)顆粒在滾動(dòng)和躍移運(yùn)動(dòng)時(shí)都會(huì)受到水下有效重力、水流拖曳力及上舉力:

    泥沙顆粒水下有效重力G0:

    水流拖曳力FD:

    上舉力FL:

    式中:CD為拖曳力系數(shù);CL為上舉力系數(shù);d為泥沙顆粒粒徑,m;u為水流作用流速,m∕s;us為推移質(zhì)運(yùn)動(dòng)速度,m∕s。

    根據(jù)周雙[19]的研究,水流作用在滾動(dòng)推移質(zhì)顆粒的作用流速可取為5.58u*,滾動(dòng)推移質(zhì)的運(yùn)動(dòng)速度為:

    躍移推移質(zhì)顆粒的運(yùn)動(dòng)軌跡為曲線,在垂直于水流方向,水流作用在顆粒上的流速與顆粒運(yùn)動(dòng)高度有關(guān),因此,顆粒所受的上舉力和拖曳力不斷變化,bagnold[12]由躍移顆粒運(yùn)動(dòng)過程中的動(dòng)量守恒得到躍移顆粒在沿水流運(yùn)動(dòng)方向受到的平均作用力為:

    式中:α為摩擦角;對(duì)于天然泥沙,tanα可取為0.63[11];其他符號(hào)意義同前。

    然后基于水流功率理論推導(dǎo)得到單位時(shí)間水流搬運(yùn)躍移推移質(zhì)泥沙所需做的功為:

    而滾動(dòng)推移質(zhì)顆粒在運(yùn)動(dòng)過程中還受到床面的滾動(dòng)摩擦力:

    式中:η為滾動(dòng)摩擦系數(shù),對(duì)于天然泥沙,通常取為0.63。

    假定在水流運(yùn)動(dòng)方向,滾動(dòng)推移質(zhì)顆粒做勻速運(yùn)動(dòng),則其在水流方向的受力平衡條件為:

    綜上可知,由于滾動(dòng)推移質(zhì)和躍移推移質(zhì)的運(yùn)動(dòng)軌跡不同,導(dǎo)致其受到的水流拖曳力和上舉力不同,在沿水流運(yùn)動(dòng)方向,滾動(dòng)和躍移推移質(zhì)顆粒受到的水流作用力分別為ˉFx和FD,因此,水流搬運(yùn)滾動(dòng)推移質(zhì)和躍移推移質(zhì)消耗的能量也應(yīng)有所差異。

    惠遇甲[20]利用高速攝影,通過水槽泥沙運(yùn)動(dòng)試驗(yàn),得到了推移質(zhì)中滾動(dòng)和躍移顆粒所占百分比與水流強(qiáng)度θ的關(guān)系式:

    式中:P1、P2分別為滾動(dòng)和躍移顆粒百分?jǐn)?shù);其他符號(hào)意義同前。

    一些學(xué)者將躍移作為推移質(zhì)運(yùn)動(dòng)的基本形式[1],但由式(21)可見,在水流強(qiáng)度較低時(shí),水流強(qiáng)度不足以使顆粒起躍,推移質(zhì)以滾動(dòng)為主;隨著水流強(qiáng)度的增大,滾動(dòng)顆粒逐漸減少;當(dāng)θ=0.6 時(shí),推移質(zhì)中滾動(dòng)顆粒所占百分比為25%,滾動(dòng)推移質(zhì)對(duì)總的推移質(zhì)能量消耗的影響不可忽略。因此,在計(jì)算水流搬運(yùn)推移質(zhì)所做的功時(shí),需要分別考慮水流搬運(yùn)滾動(dòng)推移質(zhì)和躍移推移質(zhì)所做的功,則單位時(shí)間水流搬運(yùn)推移質(zhì)消耗的總能量可表示為:

    式中:Es2為水流搬運(yùn)滾動(dòng)推移質(zhì)所消耗的能量。

    在水流泥沙運(yùn)動(dòng)方向,水流通過水流拖曳力對(duì)滾動(dòng)泥沙做功,則單位時(shí)間單位面積床面上水流搬運(yùn)滾動(dòng)推移質(zhì)泥沙所做的功可表示為:

    式中:N為單位床面面積上滾動(dòng)推移質(zhì)顆粒數(shù)量;FDi、usi分別為第i顆滾動(dòng)泥沙受到的水流拖曳力和顆粒的運(yùn)動(dòng)速度;其他符號(hào)意義同前。

    將式(14)及式(16)代入式(24)可得:

    因式中:

    則式(25)可簡(jiǎn)化為:

    式中:滾動(dòng)推移質(zhì)在運(yùn)動(dòng)過程中受到床面的遮蔽作用較小,其拖曳力系數(shù)CD可取為1.306[21];其他符號(hào)意義同前。

    將式(15)、式(16)及式(20)代入式(24)可得單位時(shí)間單位面積床面上水流搬運(yùn)滾動(dòng)推移質(zhì)泥沙所做的功也可表示為:

    對(duì)比式(18)和式(28)可知,搬運(yùn)滾動(dòng)推移質(zhì)消耗的能量還與水流上舉力有關(guān),水流搬運(yùn)滾動(dòng)推移質(zhì)和躍移推移質(zhì)消耗的能量有所差異,因此,式(18)并不適用于計(jì)算搬運(yùn)滾動(dòng)推移質(zhì)消耗的能量,水流搬運(yùn)滾動(dòng)推移質(zhì)消耗的能量可由式(27)或式(28)計(jì)算。

    將式(18)和式(27)代入式(23)可得水流搬運(yùn)推移質(zhì)消耗的能量為:

    黃才安[13]根據(jù)能量平衡原理得到輸沙阻力系數(shù)為:

    將式(29)代入式(30)可得輸沙阻力系數(shù)為:

    將式(12)和式(31)代入式(3)可得床面總阻力系數(shù)為:

    采用陸永軍[22-26]等人共136組非均勻沙推移質(zhì)輸移試驗(yàn)資料對(duì)式(32)做擬合處理,得到f(g'b∕g'b0) = 2.04g'b∕g'b0+ 0.42,其相關(guān)系數(shù)為0.83,則床面總阻力系數(shù)為:

    2 公式驗(yàn)證

    采用William[27]整理的Casey、Chyn S D、Ho P 等人的286 組非均勻沙試驗(yàn)資料分別按式(33)、黃才安公式[13]及何文社公式[7]計(jì)算床面阻力系數(shù),并將計(jì)算結(jié)果代入阻力變換公式[式(1)]計(jì)算斷面平均流速,結(jié)果如圖2所示??梢钥闯?,相比黃才安公式和何文社公式,本文公式與實(shí)測(cè)資料符合的更好。

    圖2 計(jì)算結(jié)果驗(yàn)證與對(duì)比Fig.2 Verification and comparison of calculation results

    為分析上述公式的計(jì)算精度,采用均方根相對(duì)誤差(RMSRE)來評(píng)價(jià)公式的計(jì)算精度:

    式中:RMSRE為實(shí)測(cè)值與計(jì)算值的均方根相對(duì)誤差;n為流速的樣本個(gè)數(shù);Xi為流速的計(jì)算值;Yi為流速的實(shí)測(cè)值;RMSRE的值越小,說明計(jì)算公式的精度越高。

    經(jīng)過計(jì)算,本文公式的RMSRE為0.094,黃才安公式的RMSRE為0.151,何文社公式的RMSRE為0.124,本文公式的計(jì)算精度明顯更高。

    3 結(jié) 論

    研究綜合考慮床面泥沙、水流條件和輸沙情況,將水流阻力劃分為沙粒阻力和輸沙阻力,基于床面粗化過程中床面表層結(jié)構(gòu)變化對(duì)水流阻力的影響,分析了推移質(zhì)運(yùn)動(dòng)特征,得到了床面粗化程度與相對(duì)輸沙率的關(guān)系,完善了推移質(zhì)輸移所消耗能量的計(jì)算公式,建立了平整輸沙床面的水流阻力計(jì)算公式。經(jīng)實(shí)測(cè)資料驗(yàn)證,并與其他學(xué)者的阻力計(jì)算公式分析對(duì)比表明,本文公式具有較高的精度,適用于有泥沙輸移的平整或沙波發(fā)育不明顯的天然沙非均勻床面。

    猜你喜歡
    粗化輸沙床面
    魚鱗狀床面粗糙特性
    對(duì)瓦里安碳纖維治療床面模型的評(píng)估
    淹沒植物明渠床面沖淤及其對(duì)水流運(yùn)動(dòng)的影響
    分段平移相滲曲線方法校準(zhǔn)網(wǎng)格粗化效果
    改進(jìn)的投影覆蓋方法對(duì)遼河河道粗糙床面分維量化研究
    油藏地質(zhì)模型粗化的方法及其適用性分析
    非均勻多孔介質(zhì)滲透率粗化的有限分析算法
    湖南省四水流域森林、徑流、輸沙變化特性
    基于分布式水文模型的流域輸沙過程模擬
    塔克拉瑪干沙漠腹地塔中地區(qū)的風(fēng)動(dòng)力環(huán)境
    久久久久精品国产欧美久久久| 亚洲成人精品中文字幕电影| 午夜福利免费观看在线| a在线观看视频网站| 色综合欧美亚洲国产小说| 男男h啪啪无遮挡| 亚洲中文字幕日韩| 宅男免费午夜| 最近最新免费中文字幕在线| 可以免费在线观看a视频的电影网站| 国产成人免费无遮挡视频| 精品国产超薄肉色丝袜足j| 国产免费av片在线观看野外av| 国产伦人伦偷精品视频| 精品午夜福利视频在线观看一区| 国产精品亚洲美女久久久| 国产一区二区三区视频了| 午夜a级毛片| 国产伦人伦偷精品视频| 一卡2卡三卡四卡精品乱码亚洲| 无限看片的www在线观看| 波多野结衣高清无吗| 他把我摸到了高潮在线观看| 午夜福利在线观看吧| 99久久综合精品五月天人人| 亚洲最大成人中文| 精品人妻在线不人妻| 最近最新中文字幕大全免费视频| 国产精品久久电影中文字幕| 女性被躁到高潮视频| 人人妻人人澡欧美一区二区 | 97超级碰碰碰精品色视频在线观看| 成人欧美大片| 老熟妇仑乱视频hdxx| 色综合亚洲欧美另类图片| 人成视频在线观看免费观看| 国内久久婷婷六月综合欲色啪| 在线十欧美十亚洲十日本专区| 午夜影院日韩av| 亚洲av电影在线进入| 99精品久久久久人妻精品| 91老司机精品| 亚洲中文字幕一区二区三区有码在线看 | 一区二区三区精品91| 精品久久久久久久毛片微露脸| 两个人免费观看高清视频| 国产成人精品无人区| www国产在线视频色| 韩国av一区二区三区四区| 亚洲av成人不卡在线观看播放网| 99国产综合亚洲精品| 国内精品久久久久久久电影| 老司机深夜福利视频在线观看| 亚洲美女黄片视频| 69av精品久久久久久| 自拍欧美九色日韩亚洲蝌蚪91| 老司机深夜福利视频在线观看| 成人三级黄色视频| 国产精品久久久久久人妻精品电影| 午夜精品久久久久久毛片777| 精品欧美一区二区三区在线| 免费观看精品视频网站| 国产精品野战在线观看| 欧美最黄视频在线播放免费| 一级毛片高清免费大全| 99国产精品一区二区三区| 侵犯人妻中文字幕一二三四区| 日本免费一区二区三区高清不卡 | 久9热在线精品视频| 很黄的视频免费| 国产激情久久老熟女| 亚洲天堂国产精品一区在线| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲精品综合一区在线观看 | 性少妇av在线| 女人爽到高潮嗷嗷叫在线视频| 精品人妻1区二区| 国产一区二区三区在线臀色熟女| 中出人妻视频一区二区| 不卡av一区二区三区| 免费少妇av软件| 亚洲一码二码三码区别大吗| 午夜久久久在线观看| 国产乱人伦免费视频| 欧美黑人精品巨大| 久久性视频一级片| 脱女人内裤的视频| 精品电影一区二区在线| 欧美日韩乱码在线| 亚洲精品国产一区二区精华液| 亚洲av五月六月丁香网| 亚洲av美国av| 极品人妻少妇av视频| 香蕉久久夜色| 美女免费视频网站| 日韩av在线大香蕉| 天堂影院成人在线观看| 日韩免费av在线播放| 亚洲av电影在线进入| 欧美黄色淫秽网站| 亚洲专区国产一区二区| 国产一区二区在线av高清观看| bbb黄色大片| 丝袜美足系列| 久久狼人影院| 男人操女人黄网站| 久久久久久久久中文| 亚洲少妇的诱惑av| 老鸭窝网址在线观看| 国产精品电影一区二区三区| 伊人久久大香线蕉亚洲五| 国产一区二区三区在线臀色熟女| 多毛熟女@视频| 亚洲熟女毛片儿| 十八禁网站免费在线| 欧美日韩福利视频一区二区| 亚洲av成人一区二区三| 亚洲国产毛片av蜜桃av| 一边摸一边抽搐一进一小说| 看片在线看免费视频| 乱人伦中国视频| 中文字幕另类日韩欧美亚洲嫩草| 正在播放国产对白刺激| 天天躁夜夜躁狠狠躁躁| 久久亚洲真实| 黄色毛片三级朝国网站| 搡老岳熟女国产| 母亲3免费完整高清在线观看| 国产人伦9x9x在线观看| 久久人人97超碰香蕉20202| 极品教师在线免费播放| 最近最新中文字幕大全电影3 | 啦啦啦免费观看视频1| 午夜福利在线观看吧| 亚洲欧美日韩高清在线视频| 丁香欧美五月| 给我免费播放毛片高清在线观看| 香蕉久久夜色| 亚洲久久久国产精品| 日本一区二区免费在线视频| 亚洲自偷自拍图片 自拍| 日本 av在线| 亚洲五月婷婷丁香| 成人永久免费在线观看视频| 国产精品免费一区二区三区在线| 日本撒尿小便嘘嘘汇集6| 女人爽到高潮嗷嗷叫在线视频| 韩国精品一区二区三区| √禁漫天堂资源中文www| videosex国产| 91精品三级在线观看| 久久久水蜜桃国产精品网| 久久中文字幕人妻熟女| 黄片大片在线免费观看| 亚洲男人的天堂狠狠| 国产伦一二天堂av在线观看| 纯流量卡能插随身wifi吗| 欧美黄色片欧美黄色片| 18美女黄网站色大片免费观看| 人妻久久中文字幕网| 精品国产国语对白av| 国产高清视频在线播放一区| 日韩欧美三级三区| 中亚洲国语对白在线视频| 欧美成人免费av一区二区三区| 亚洲国产精品合色在线| 国产黄a三级三级三级人| 十分钟在线观看高清视频www| 国产乱人伦免费视频| 国产一区二区三区综合在线观看| 校园春色视频在线观看| 国产精品久久久久久亚洲av鲁大| 日本 欧美在线| 啦啦啦韩国在线观看视频| 人人妻人人澡欧美一区二区 | 欧美激情久久久久久爽电影 | 亚洲av片天天在线观看| 国产一卡二卡三卡精品| 在线av久久热| 老汉色∧v一级毛片| 久久中文字幕人妻熟女| 欧美一级a爱片免费观看看 | 国产三级在线视频| 给我免费播放毛片高清在线观看| 亚洲国产精品合色在线| 啦啦啦 在线观看视频| 可以免费在线观看a视频的电影网站| 亚洲三区欧美一区| 香蕉丝袜av| 亚洲精品国产色婷婷电影| 日韩免费av在线播放| 大型av网站在线播放| 免费观看精品视频网站| av视频免费观看在线观看| 18禁黄网站禁片午夜丰满| 757午夜福利合集在线观看| 亚洲男人天堂网一区| 日韩一卡2卡3卡4卡2021年| 国产乱人伦免费视频| 黄色视频不卡| 国产日韩一区二区三区精品不卡| 无遮挡黄片免费观看| 久久性视频一级片| 又黄又粗又硬又大视频| 黄色视频不卡| 两个人看的免费小视频| 国产不卡一卡二| 国产熟女午夜一区二区三区| 九色亚洲精品在线播放| 免费看美女性在线毛片视频| 啦啦啦韩国在线观看视频| 美女免费视频网站| aaaaa片日本免费| 国产亚洲av高清不卡| 大陆偷拍与自拍| 国产高清videossex| 97碰自拍视频| 日韩高清综合在线| 一边摸一边抽搐一进一出视频| 日韩三级视频一区二区三区| 欧美中文综合在线视频| 乱人伦中国视频| 日韩高清综合在线| 国产亚洲精品一区二区www| 亚洲激情在线av| 99久久国产精品久久久| 国产一级毛片七仙女欲春2 | 午夜久久久久精精品| 国产精品久久久久久亚洲av鲁大| 好男人在线观看高清免费视频 | 亚洲片人在线观看| 国产麻豆成人av免费视频| 欧美精品亚洲一区二区| 欧美激情高清一区二区三区| 热re99久久国产66热| 日韩欧美一区二区三区在线观看| www国产在线视频色| 每晚都被弄得嗷嗷叫到高潮| 色播亚洲综合网| 国产1区2区3区精品| 国产精品久久久人人做人人爽| 久久这里只有精品19| 中文亚洲av片在线观看爽| 黄频高清免费视频| 9191精品国产免费久久| 99国产精品免费福利视频| 色综合欧美亚洲国产小说| 天天躁夜夜躁狠狠躁躁| 一区在线观看完整版| 久久久久国产精品人妻aⅴ院| 99在线人妻在线中文字幕| 日本三级黄在线观看| 婷婷精品国产亚洲av在线| 久久久久九九精品影院| 亚洲色图av天堂| 午夜福利免费观看在线| 可以在线观看毛片的网站| 免费无遮挡裸体视频| 久久精品影院6| 亚洲国产中文字幕在线视频| 亚洲成人精品中文字幕电影| 18禁黄网站禁片午夜丰满| 国产亚洲精品久久久久久毛片| 久久人妻熟女aⅴ| 窝窝影院91人妻| 两个人免费观看高清视频| 最新在线观看一区二区三区| 日韩精品青青久久久久久| 男人的好看免费观看在线视频 | 大陆偷拍与自拍| 欧美中文日本在线观看视频| 999久久久国产精品视频| 午夜福利在线观看吧| 一本大道久久a久久精品| 一二三四社区在线视频社区8| 亚洲av电影不卡..在线观看| 丁香欧美五月| 国产成人欧美在线观看| 午夜精品国产一区二区电影| 成年人黄色毛片网站| 国产精品98久久久久久宅男小说| 狠狠狠狠99中文字幕| 亚洲午夜理论影院| 99在线人妻在线中文字幕| 亚洲天堂国产精品一区在线| av在线播放免费不卡| 精品无人区乱码1区二区| 亚洲成人免费电影在线观看| 又大又爽又粗| 国产精品99久久99久久久不卡| or卡值多少钱| 午夜福利成人在线免费观看| 免费观看精品视频网站| 变态另类成人亚洲欧美熟女 | 亚洲男人的天堂狠狠| 夜夜看夜夜爽夜夜摸| 久久香蕉国产精品| 在线观看免费日韩欧美大片| 中文字幕最新亚洲高清| 久久国产精品影院| 18禁国产床啪视频网站| 一区在线观看完整版| 韩国av一区二区三区四区| 亚洲精品久久成人aⅴ小说| 午夜福利视频1000在线观看 | or卡值多少钱| 亚洲avbb在线观看| 国产99白浆流出| 90打野战视频偷拍视频| 校园春色视频在线观看| 亚洲中文字幕一区二区三区有码在线看 | 美女高潮喷水抽搐中文字幕| e午夜精品久久久久久久| 女性被躁到高潮视频| 90打野战视频偷拍视频| 亚洲av片天天在线观看| 天天添夜夜摸| 99久久99久久久精品蜜桃| 女生性感内裤真人,穿戴方法视频| 欧美大码av| 男女做爰动态图高潮gif福利片 | 亚洲国产毛片av蜜桃av| 久久精品国产99精品国产亚洲性色 | 久久久水蜜桃国产精品网| 久久热在线av| 精品乱码久久久久久99久播| 最新美女视频免费是黄的| 黄色女人牲交| 丰满人妻熟妇乱又伦精品不卡| 男女下面进入的视频免费午夜 | 制服人妻中文乱码| 日本 av在线| 老司机福利观看| 91成人精品电影| 精品国产一区二区久久| 精品人妻1区二区| 热99re8久久精品国产| 免费在线观看影片大全网站| 给我免费播放毛片高清在线观看| 午夜两性在线视频| 高清毛片免费观看视频网站| 美女国产高潮福利片在线看| av网站免费在线观看视频| 少妇熟女aⅴ在线视频| 久久中文字幕人妻熟女| 熟妇人妻久久中文字幕3abv| 免费女性裸体啪啪无遮挡网站| 欧美乱码精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 国产成+人综合+亚洲专区| 国产成人精品在线电影| 国产亚洲精品一区二区www| 黑丝袜美女国产一区| 99久久99久久久精品蜜桃| 91成年电影在线观看| 99精品在免费线老司机午夜| 视频区欧美日本亚洲| 免费一级毛片在线播放高清视频 | 搞女人的毛片| 男女之事视频高清在线观看| 欧美国产精品va在线观看不卡| 国产精品av久久久久免费| 国产熟女xx| 黑丝袜美女国产一区| 欧美日本视频| 男女之事视频高清在线观看| 美女国产高潮福利片在线看| 亚洲aⅴ乱码一区二区在线播放 | 超碰成人久久| 精品国产乱码久久久久久男人| 黑人欧美特级aaaaaa片| 日韩精品青青久久久久久| 久久久久国内视频| 夜夜看夜夜爽夜夜摸| 大香蕉久久成人网| 女人被躁到高潮嗷嗷叫费观| 精品国产一区二区久久| 亚洲自拍偷在线| 国产区一区二久久| 精品国产一区二区三区四区第35| 午夜久久久在线观看| 两个人免费观看高清视频| 国产精品永久免费网站| 国产精品,欧美在线| 国产日韩一区二区三区精品不卡| АⅤ资源中文在线天堂| 精品卡一卡二卡四卡免费| 大陆偷拍与自拍| 成人永久免费在线观看视频| 男女下面插进去视频免费观看| 久久久久国产精品人妻aⅴ院| 亚洲aⅴ乱码一区二区在线播放 | 亚洲av五月六月丁香网| 日韩大码丰满熟妇| 久久久久久免费高清国产稀缺| 青草久久国产| 91字幕亚洲| 性少妇av在线| 亚洲国产看品久久| 欧美中文综合在线视频| 99热只有精品国产| 亚洲色图综合在线观看| 国产成人精品在线电影| 欧美性长视频在线观看| 18美女黄网站色大片免费观看| 国产亚洲精品久久久久久毛片| 两性夫妻黄色片| 国产高清视频在线播放一区| 免费不卡黄色视频| 国产又爽黄色视频| 黄色丝袜av网址大全| 91麻豆精品激情在线观看国产| 国产精品一区二区免费欧美| 国产黄a三级三级三级人| 丝袜美足系列| 满18在线观看网站| 午夜两性在线视频| 国产精品美女特级片免费视频播放器 | 亚洲人成电影观看| 每晚都被弄得嗷嗷叫到高潮| 97超级碰碰碰精品色视频在线观看| 日韩成人在线观看一区二区三区| 精品欧美一区二区三区在线| 国产又爽黄色视频| 校园春色视频在线观看| 精品乱码久久久久久99久播| 成人欧美大片| 中文字幕另类日韩欧美亚洲嫩草| 九色亚洲精品在线播放| 真人一进一出gif抽搐免费| 天天躁夜夜躁狠狠躁躁| 中文字幕色久视频| 精品人妻在线不人妻| 9色porny在线观看| 色综合亚洲欧美另类图片| 老司机福利观看| 亚洲国产精品sss在线观看| 精品欧美一区二区三区在线| 国产精品98久久久久久宅男小说| 丝袜美腿诱惑在线| 人人澡人人妻人| 丰满的人妻完整版| 性欧美人与动物交配| 欧美+亚洲+日韩+国产| 亚洲一卡2卡3卡4卡5卡精品中文| 黑人巨大精品欧美一区二区mp4| 国内毛片毛片毛片毛片毛片| 91精品国产国语对白视频| 免费在线观看完整版高清| 精品久久久久久成人av| 精品久久久久久,| 日韩一卡2卡3卡4卡2021年| 老司机靠b影院| 久久久久国内视频| 日韩国内少妇激情av| 国产成人精品在线电影| 中文字幕久久专区| 久久天堂一区二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91| 波多野结衣高清无吗| 国产欧美日韩一区二区精品| 精品不卡国产一区二区三区| 亚洲七黄色美女视频| 丝袜美足系列| 亚洲熟妇熟女久久| 欧美另类亚洲清纯唯美| 村上凉子中文字幕在线| 两性夫妻黄色片| 侵犯人妻中文字幕一二三四区| 国产成人精品无人区| 亚洲一区二区三区不卡视频| 亚洲中文av在线| 免费观看精品视频网站| 精品久久久久久久毛片微露脸| 欧美日韩瑟瑟在线播放| 给我免费播放毛片高清在线观看| 国产男靠女视频免费网站| 成人亚洲精品av一区二区| 亚洲av美国av| 一本综合久久免费| 18禁裸乳无遮挡免费网站照片 | 国产在线观看jvid| 亚洲情色 制服丝袜| 免费高清在线观看日韩| 一二三四在线观看免费中文在| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产精品999在线| svipshipincom国产片| 变态另类丝袜制服| 亚洲国产精品合色在线| 波多野结衣一区麻豆| 久久久国产欧美日韩av| 国产日韩一区二区三区精品不卡| 久久久精品欧美日韩精品| 久久精品亚洲熟妇少妇任你| 亚洲成av人片免费观看| 日韩成人在线观看一区二区三区| 久久午夜综合久久蜜桃| 国产精品爽爽va在线观看网站 | 亚洲成人久久性| 亚洲成国产人片在线观看| 1024香蕉在线观看| 欧美日本视频| 国产成人影院久久av| 波多野结衣高清无吗| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品美女久久av网站| 两性午夜刺激爽爽歪歪视频在线观看 | av欧美777| 我的亚洲天堂| 亚洲欧美一区二区三区黑人| 999精品在线视频| 久久久久久亚洲精品国产蜜桃av| 午夜免费观看网址| 91九色精品人成在线观看| av网站免费在线观看视频| 一进一出抽搐gif免费好疼| 亚洲成国产人片在线观看| 精品乱码久久久久久99久播| 亚洲一区二区三区不卡视频| 嫁个100分男人电影在线观看| 日韩精品青青久久久久久| 午夜成年电影在线免费观看| 99在线人妻在线中文字幕| 91精品国产国语对白视频| 成在线人永久免费视频| 熟妇人妻久久中文字幕3abv| 69av精品久久久久久| 国产亚洲av高清不卡| 精品无人区乱码1区二区| 日韩欧美一区二区三区在线观看| 精品卡一卡二卡四卡免费| 免费看十八禁软件| 亚洲 国产 在线| 欧美国产日韩亚洲一区| 欧美日韩中文字幕国产精品一区二区三区 | 日韩 欧美 亚洲 中文字幕| 黄色a级毛片大全视频| 欧美日韩乱码在线| 老司机靠b影院| 午夜福利成人在线免费观看| 国产三级黄色录像| 亚洲精品一卡2卡三卡4卡5卡| 午夜老司机福利片| 日本 av在线| 中文字幕高清在线视频| 这个男人来自地球电影免费观看| 国产三级在线视频| 黄色a级毛片大全视频| 久久精品国产亚洲av高清一级| 亚洲国产欧美日韩在线播放| 99久久综合精品五月天人人| 欧美激情极品国产一区二区三区| 国产精品99久久99久久久不卡| 非洲黑人性xxxx精品又粗又长| 两个人免费观看高清视频| 精品国产国语对白av| 一本久久中文字幕| 国产在线观看jvid| 日日爽夜夜爽网站| 麻豆av在线久日| 亚洲精品一区av在线观看| 国产高清videossex| 免费看美女性在线毛片视频| 不卡一级毛片| 精品午夜福利视频在线观看一区| 亚洲精品国产精品久久久不卡| 欧美精品亚洲一区二区| 欧美色欧美亚洲另类二区 | 欧美乱码精品一区二区三区| 国产色视频综合| 美女午夜性视频免费| 久久人妻av系列| 亚洲精品粉嫩美女一区| 91在线观看av| av网站免费在线观看视频| 老司机靠b影院| 欧洲精品卡2卡3卡4卡5卡区| 嫩草影院精品99| 国产亚洲精品久久久久久毛片| 性欧美人与动物交配| 亚洲精品美女久久av网站| 国产精品一区二区精品视频观看| 国产精品秋霞免费鲁丝片| 亚洲视频免费观看视频| 国产精品亚洲美女久久久| 欧美人与性动交α欧美精品济南到| 久热这里只有精品99| 国产亚洲精品一区二区www| 麻豆国产av国片精品| 久热这里只有精品99| 亚洲一区二区三区不卡视频| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色女人牲交| 亚洲自偷自拍图片 自拍| 免费在线观看黄色视频的| 91麻豆精品激情在线观看国产| 免费在线观看视频国产中文字幕亚洲| 亚洲人成电影免费在线| 国产精品久久久人人做人人爽| 国产成+人综合+亚洲专区| 精品熟女少妇八av免费久了| 1024香蕉在线观看| 又黄又爽又免费观看的视频| 老司机午夜十八禁免费视频| 欧美黄色片欧美黄色片| 成年版毛片免费区| 日本vs欧美在线观看视频| 99精品欧美一区二区三区四区| 成人国产综合亚洲| av片东京热男人的天堂| 亚洲熟妇中文字幕五十中出| 国产精品免费视频内射| 成年女人毛片免费观看观看9| 黄色视频不卡|