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

    城門洞型隧洞明滿流水力瞬變計(jì)算研究

    2018-11-01 06:17:24鞠小明孫立成
    中國(guó)農(nóng)村水利水電 2018年10期
    關(guān)鍵詞:水深水力隧洞

    鞠 鋮,鞠小明,2,孫立成

    (1.四川大學(xué)水利水電學(xué)院,成都 610065;2.四川大學(xué)水力學(xué)與山區(qū)河流開(kāi)發(fā)保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,成都 610065)

    0 引 言

    城門洞型尾水隧洞在水電站輸水系統(tǒng)中應(yīng)用廣泛,受到下游運(yùn)行水位變化的影響,有些水電站尾水隧洞中可能會(huì)出現(xiàn)名滿交替水流現(xiàn)象,例如當(dāng)下游水位接近洞頂高程,機(jī)組增加負(fù)荷,原本是無(wú)壓流狀態(tài)的尾水隧洞中可能出現(xiàn)涌波,從而發(fā)生有壓流和無(wú)壓流(也稱明流)交替流動(dòng),對(duì)于頂部通氣條件不良的城門洞型隧洞,明滿流發(fā)生時(shí)可能還會(huì)封存部分氣腔在隧洞頂部[1]。實(shí)驗(yàn)研究發(fā)現(xiàn)明滿流發(fā)生時(shí),隧洞頂部壓力會(huì)突然升高或降低,壓力變化劇烈, 若設(shè)計(jì)控制不良可能會(huì)影響輸水隧洞運(yùn)行安全[1,2]。對(duì)輸水隧洞明滿流進(jìn)行水力瞬變計(jì)算研究,目的是確定洞頂壓力變化值以及變化過(guò)程,為輸水隧洞設(shè)計(jì)和運(yùn)行提供技術(shù)依據(jù)。由于有壓流的水力瞬變方程組和無(wú)壓流水力瞬變的方程組不同,因此在明滿流水力瞬變計(jì)算方法上,通常采用兩種計(jì)算途徑,一是將無(wú)壓流和有壓流分開(kāi)計(jì)算,各自采用對(duì)應(yīng)的瞬變流方程組進(jìn)行求解,將無(wú)壓流和有壓流交界面作為一個(gè)滑動(dòng)的激波來(lái)處理,這種方法也稱為激波擬合法[3,4]。另一種方法是基于普雷斯曼(Preissmann)窄縫概念,也就是當(dāng)無(wú)壓流變成滿流時(shí),假想隧洞頂部有一個(gè)既不增加過(guò)流斷面積也不改變過(guò)水?dāng)嗝鏉裰艿恼p,此時(shí)的有壓流斷面隧洞頂部壓力就是窄縫中的水深,這樣就可以采用無(wú)壓流瞬變流方程組(圣維南方程組)來(lái)求解有壓流的斷面壓力和流量,這樣的方法也稱為窄縫法或激波捕捉法[5,3]。比較而言,窄縫法比激波擬合法更加實(shí)用,在城市污水管網(wǎng)和水電站輸水管道的水力瞬變模擬中得到了應(yīng)用[6~10]。對(duì)于城門洞型輸水隧洞,即使采用窄縫法進(jìn)行計(jì)算,由于水位處于拱頂部位時(shí)的計(jì)算變量如面積和濕周等不僅與水深有關(guān),還與自由水面形成的拱頂夾角有關(guān),這使得城門洞型隧洞水力瞬變計(jì)算更加復(fù)雜和不容易穩(wěn)定,因此研究城門洞型輸水隧洞的水力瞬變方法既有理論指導(dǎo)意義又有現(xiàn)實(shí)應(yīng)用價(jià)值。

    1 數(shù)學(xué)模型和基本方程

    1.1 有壓流水力瞬變計(jì)算基本方程和求解方法

    有壓流水力瞬變計(jì)算的基本方程見(jiàn)下式。

    運(yùn)動(dòng)方程:

    (1)

    連續(xù)方程:

    (2)

    式中:H為壓頭;V為流速;x為從管段左端起算的距離;g為重力加速度;f為沿程損失系數(shù);D為管徑;a為水錘波速;t為時(shí)間。

    將基本方程沿正特征線和負(fù)特征線分別積分可得正負(fù)特征方程,利用邊界條件就可以求解有壓流隧洞各斷面節(jié)點(diǎn)的壓力值和流量值。

    1.2 無(wú)壓流基本方程和求解方法

    無(wú)壓流水力瞬變計(jì)算的基本方程。

    運(yùn)動(dòng)方程:

    (3)

    連續(xù)方程:

    (4)

    式中:y為無(wú)壓流水深;V和t意義同前;A為無(wú)壓流過(guò)水?dāng)嗝娣e;B為無(wú)壓流水面寬;S0為無(wú)壓流隧洞底坡;Sf為無(wú)壓流阻力損失引起的能量坡度, 由下式計(jì)算:

    方程(3)和(4)也稱為無(wú)壓流瞬變流計(jì)算圣維南(Saint-Venanf)方程組,式中Sf由公式(5)計(jì)算。

    (5)

    式中:n為曼寧粗糙度系數(shù);R為水力半徑,由公式(6)計(jì)算:

    (6)

    式中:X為無(wú)壓流過(guò)水?dāng)嗝鏉裰堋?/p>

    明流水力瞬變計(jì)算方法主要有差分法、特征線法、瞬時(shí)流態(tài)法、微幅波理論法等,其中差分法又分為顯式差分法和隱式差分法。也可以將差分法和特征線法結(jié)合起來(lái)使用,即在邊界上采用特征線法,內(nèi)部斷面采用差分法,這樣的方法Chaudhry M H稱之為擴(kuò)散法(diffusive scheme)[5]。無(wú)論采用哪種計(jì)算方法,在無(wú)壓流瞬變過(guò)程中的邊界條件方程和計(jì)算穩(wěn)定性都是無(wú)法回避和必須考慮的重要問(wèn)題,有關(guān)明流計(jì)算方法、差分格式、2維計(jì)算方法等可參考有關(guān)文獻(xiàn)資料[11~14]。

    在城門洞型隧洞明流計(jì)算中,當(dāng)瞬時(shí)水深超過(guò)城門洞直墻高度進(jìn)入弧形頂拱時(shí),瞬時(shí)水面寬度、濕周和過(guò)水面積、無(wú)壓流波速等不僅是瞬時(shí)水深yp的函數(shù),同時(shí)還是水面扇形瞬時(shí)夾角θ的函數(shù),如圖1所示,相應(yīng)的偏導(dǎo)數(shù)項(xiàng)將增加一個(gè)水面夾角變量,這也是城門洞型隧洞水力瞬變過(guò)程計(jì)算更加復(fù)雜的原因之一。

    圖1 窄縫法計(jì)算城門洞隧洞斷面圖Fig.1 Section of rectangular tunnel with arch crown and a slot

    下列公式(7)~(10)用于計(jì)算無(wú)壓流水位超過(guò)直墻高度位于弧形頂拱部位時(shí)的明流瞬時(shí)變量。

    扇形水面夾角θ由下式計(jì)算:

    (7)

    式中:yP為隧洞瞬時(shí)水深;Hc為城門洞隧洞直墻高度;Rc為城門洞隧洞拱頂半徑;Bc為城門洞隧洞底寬。

    瞬時(shí)水面寬由下式計(jì)算:

    (8)

    式中:θ為扇形水面夾角;其他符號(hào)意義同前。

    瞬時(shí)濕周由下式計(jì)算:

    X=Bc+2Hc+Rc(α-θ)

    (9)

    式中:α為城門洞拱頂夾角(弧度);其他符號(hào)意義同前。

    瞬時(shí)過(guò)水面積由下式計(jì)算:

    (10)

    式中:A為城門洞隧洞瞬時(shí)過(guò)水?dāng)嗝娣e;其他符號(hào)意義同前。

    城門洞型隧洞水力瞬變計(jì)算時(shí)需要知道系統(tǒng)內(nèi)所有斷面上的初始恒定狀態(tài),即需要知道每一節(jié)點(diǎn)處的水深和流速,無(wú)壓明流隧洞初始水面曲線可由明渠漸變流常微分方程計(jì)算確定[5]。

    (11)

    式中:Q為隧洞明流過(guò)流量;其他符號(hào)意義同前。

    用四階龍格-庫(kù)塔法求解上述方程,已知X=Xi處水深為yi,那么X=Xi+1處水深為:

    (12)

    式中:

    a1=Δxf(xi,yi)

    a4=Δxf(xi+Δx,yi+a3)

    隧洞內(nèi)明流各斷面的初始水位可由上述方程(11)和(12)計(jì)算確定。

    2 普雷斯曼窄縫法計(jì)算明滿流水力瞬變方法

    2.1 普雷斯曼(Preissmann)窄縫法計(jì)算原理

    當(dāng)城門洞隧洞出口水位接近洞頂高程時(shí),若隧洞內(nèi)水流發(fā)生水力瞬變,此時(shí)就可能發(fā)生明滿流現(xiàn)象(也稱明滿交替水流)。普雷斯曼(Preissmann)提出了一個(gè)有趣的計(jì)算方法,該方法假想隧洞內(nèi)水流由無(wú)壓流轉(zhuǎn)變?yōu)橛袎毫鲿r(shí),隧洞頂部通氣良好,不存在封存氣囊或氣腔,并假想隧洞頂部存在一個(gè)窄縫(參見(jiàn)圖1),該窄縫既不增加過(guò)流斷面積,也不增加水力半徑,當(dāng)發(fā)生有壓流時(shí),只要滿足無(wú)壓流波速(celerity of surface wave)與有壓流水錘波速(velocity of water-hammer)相等的方法來(lái)選擇窄縫寬度[15,,5],這樣窄縫中的水深變量就是有壓流的壓力水頭變量,并且物理概念也完全一致,如此無(wú)壓流和有壓流就可以用同一形式的圣維南方程組來(lái)描述,求解這個(gè)方程就可以進(jìn)行明滿流水力瞬變的計(jì)算,這就是明滿流水力瞬變計(jì)算的基本原理和方法,該方法在城市下水道的水力瞬變計(jì)算中得到了應(yīng)用。由于該方法沒(méi)有考慮城門洞洞頂可能存在的氣囊或氣腔內(nèi)空氣體積的變化,僅適用于洞頂通氣條件良好的隧洞,本文采用窄縫法的基本思路和原理對(duì)城門洞型輸水隧洞進(jìn)行水力瞬變計(jì)算研究。

    2.2 城門洞隧洞明滿流算例

    某大型發(fā)電輸水尾水隧洞采用城門洞型斷面,隧洞底寬17.5 m,直墻高17.04 m,頂拱夾角118.15°,隧洞斷面最大凈空高度22.0 m,隧洞長(zhǎng)300.0 m,底坡0.00,如圖2所示。由于隧洞斷面尺寸大,不具備采用變頂高的通氣和排氣隧洞條件,當(dāng)電站尾水水位接近洞頂高程附近時(shí),此時(shí)上游斷面流量增加,隧洞內(nèi)將產(chǎn)生明滿交替水流;或者隧洞內(nèi)原為有壓流動(dòng),當(dāng)上游斷面流量減小,也將產(chǎn)生明滿交替水流。

    圖2 某輸水隧洞布置圖Fig.2 Profile of an engineering tunnel with arch crown

    采用前述的計(jì)算原理和計(jì)算方法,隧洞全長(zhǎng)300 m,共分為50段,也即6 m一小段,包括首尾共有51個(gè)計(jì)算節(jié)點(diǎn)。為了能夠觀察隧洞全長(zhǎng)范圍內(nèi)不同斷面的水力瞬變情況,并且兩個(gè)觀察斷面的水力瞬變量能有比較大的區(qū)別,方便用不同曲線來(lái)表示,從上游0 m斷面開(kāi)始每間隔60 m選擇一個(gè)觀察特征斷面,因此分別選取距離上游0、120、240 m三個(gè)斷面作為計(jì)算觀察特征斷面,研究這些特征斷面發(fā)生明滿流時(shí)的水深變化過(guò)程,而隧洞下游出口斷面始終保持初始恒定狀態(tài)的水深。計(jì)算三種發(fā)生明滿流的工況,如表1所示。

    表1 發(fā)生明滿流的計(jì)算工況Tab.1 Free-surface-pressure flow transient condition of calculation

    工況1是隧洞內(nèi)初始流量為0,下游斷面水深21.0 m,距離洞頂1.0 m,隧洞內(nèi)初始狀態(tài)為無(wú)壓流,隧洞上游斷面流量在20 s內(nèi)線性增加至1 600 m3/s,隧洞內(nèi)各斷面水深變化過(guò)程如圖3所示。工況2是下游斷面水深27.862 m,超過(guò)洞頂5.862 m,初始工況是有壓流,初始流量1 600 m3/s,上游斷面在10 s內(nèi)流量線性減小至零,隧洞內(nèi)各特征斷面的水深變化如圖4所示。工況3是將第一種工況的流量變化時(shí)間縮短為10 s,隧洞內(nèi)各特征斷面的水深變化如圖5所示。

    圖3 上游斷面流量增加特征斷面水深變化過(guò)程(工況1)Fig.3 Process of water depth variation on particular section with flow discharge increasing upstream (condition 1)

    圖4 上游斷面流量減小特征斷面水深變化過(guò)程(工況2)Fig.4 Process of water depth variation on particular section with flow discharge decreasing upstream (condition 2)

    圖5 上游斷面流量增加特征斷面水深變化過(guò)程(工況3)Fig.5 Process of water depth variation on particular section with flow discharge increasing upstream (condition 3)

    2.3 算例計(jì)算成果分析

    從工況1的計(jì)算成果圖3可以看出,隧洞初始狀態(tài)為無(wú)壓流,隨著上游斷面流量增加,在上游流量開(kāi)始變化后的1.74 s首先是上游第一個(gè)斷面(距離上游0 m斷面)變成有壓流,隨后隧洞處于有壓流狀態(tài),并且當(dāng)從無(wú)壓流變?yōu)橛袎毫鲿r(shí),隧洞水深(拱頂壓力)急劇變化,之后維持一段時(shí)間的有壓流狀態(tài),當(dāng)上游斷面流量增加結(jié)束,不再有補(bǔ)充流量時(shí),隧洞內(nèi)水流恢復(fù)無(wú)壓狀態(tài),從圖中可以看出距離流量變化斷面越近,壓力變化越大,圖3中距離上游0 m斷面窄縫內(nèi)最大水深43.755 m,算例隧洞最大凈高22.0 m,相當(dāng)于頂拱在明滿交替時(shí)壓力劇變21.755 m,這種劇變壓力類似于水錘壓力,是瞬間產(chǎn)生的,它不同于靜水壓強(qiáng)分布規(guī)律,對(duì)隧洞頂部的襯砌會(huì)造成不利影響,嚴(yán)重時(shí)可能會(huì)影響襯砌結(jié)構(gòu)安全。圖4為初始狀態(tài)是有壓流,上游斷面流量減小,可以看出最先是距離上游0 m斷面的水深減小,然后是后續(xù)斷面水深跟著減小,從有壓流變?yōu)闊o(wú)壓流時(shí)沒(méi)有發(fā)生壓力劇變(模型中不考慮洞頂空氣體積變化),距離上游0 m斷面最小水深17.238 m,基本處于城門洞直墻頂高程附近(算例直墻高17.04 m),也就是水力瞬變處于拱頂變水面寬的穹頂內(nèi),說(shuō)明前述計(jì)算方法中關(guān)于水位處于穹頂內(nèi)的計(jì)算原理是可行的,隨著無(wú)壓流過(guò)渡過(guò)程的發(fā)展,因?yàn)樗矶醋钕掠蔚臄嗝媸怯袎旱?窄縫內(nèi)水深27.862 m,隧洞凈高22.0 m),因此在流量變化后的11.38 s時(shí)刻,首先是靠近下游斷面的特征斷面出現(xiàn)有壓流,如圖4中的紅色線(圖中更靠近下游的虛線最小水位均高于22.0 m,該斷面是有壓狀態(tài)),然后上游斷面也開(kāi)始出現(xiàn)有壓流,但是依然是最靠近流量變化斷面的節(jié)點(diǎn)出現(xiàn)最大水深,如圖4中的距離上游0 m斷面窄縫內(nèi)最大水深38.729 m,即明滿流交替時(shí)壓力劇變16.729 m。表2是初始工況是有壓流,上游斷面流量減小后各特征斷面的最大和最小水深。

    表2 上游斷面流量減小各特征斷面最大和最小水深 m

    圖5與圖3類似,區(qū)別是上游斷面流量變化比圖3快1倍,窄縫內(nèi)最大水深達(dá)到53.896 m,即明滿流交替時(shí)壓力劇變達(dá)到31.896 m,說(shuō)明流量變化越快,明滿流交替產(chǎn)生的壓力變化越大,兩種情況的拱頂最大壓力數(shù)據(jù)比較如表3所示。

    表3 不同流量變化時(shí)間各特征斷面水深和拱頂壓力比較(工況1/工況3) m

    3 結(jié) 語(yǔ)

    算例研究表明采用窄縫法計(jì)算城門洞型隧洞內(nèi)的明滿交替水力瞬變過(guò)程是可行的;在不考慮洞頂可能存在的氣囊或氣腔內(nèi)空氣體積變化的情形下,隧洞內(nèi)水流由無(wú)壓流變?yōu)橛袎毫鲿r(shí)將產(chǎn)生較大的壓力變化,由有壓流變?yōu)闊o(wú)壓流時(shí)壓力變化較小;隧洞內(nèi)斷面流量變化越快,明滿流交替產(chǎn)生的壓力變化也越大;靠近流量變化越近的斷面水力瞬變程度越急劇。

    猜你喜歡
    水深水力隧洞
    水力全開(kāi)
    書(shū)法靜水深流
    河北水利(2022年10期)2022-12-29 11:48:12
    基于水深分段選擇因子的多光譜影像反演水深
    隧洞止水帶安裝質(zhì)量控制探討
    滇中引水工程大理段首條隧洞順利貫通
    自密實(shí)混凝土在水工隧洞襯砌中的應(yīng)用
    球墨鑄鐵管的水力計(jì)算
    水力噴射壓裂中環(huán)空水力封隔全尺寸實(shí)驗(yàn)
    GPS RTK技術(shù)在水深測(cè)量中的應(yīng)用
    浸入式水深監(jiān)測(cè)儀器的設(shè)計(jì)
    99热网站在线观看| 欧美一级a爱片免费观看看| a级一级毛片免费在线观看| 久久鲁丝午夜福利片| 日韩欧美 国产精品| 色吧在线观看| 久久午夜福利片| 麻豆乱淫一区二区| 精品久久久噜噜| 欧美区成人在线视频| 久久99蜜桃精品久久| 一级a做视频免费观看| 91久久精品国产一区二区成人| 亚洲精品456在线播放app| 国产色婷婷99| 欧美xxⅹ黑人| 一级毛片 在线播放| 亚洲,一卡二卡三卡| 99热全是精品| 国产成人午夜福利电影在线观看| 日产精品乱码卡一卡2卡三| 国产国拍精品亚洲av在线观看| 午夜免费鲁丝| 欧美3d第一页| 久久久久国产网址| 久久久久精品久久久久真实原创| 国产免费视频播放在线视频| 纯流量卡能插随身wifi吗| 欧美人与善性xxx| 国产精品国产三级专区第一集| 色婷婷av一区二区三区视频| 精品国产一区二区久久| 欧美日韩一区二区视频在线观看视频在线| 精品一品国产午夜福利视频| h日本视频在线播放| 日产精品乱码卡一卡2卡三| 国产精品人妻久久久影院| 在线观看国产h片| 毛片一级片免费看久久久久| 乱系列少妇在线播放| 亚洲av在线观看美女高潮| 亚洲欧美成人精品一区二区| 亚洲国产毛片av蜜桃av| 青春草视频在线免费观看| a级片在线免费高清观看视频| 91精品一卡2卡3卡4卡| 免费大片黄手机在线观看| 亚洲第一区二区三区不卡| 男女免费视频国产| a级毛片在线看网站| 久久99蜜桃精品久久| 少妇人妻 视频| 99热网站在线观看| 色吧在线观看| 97超视频在线观看视频| 日韩中文字幕视频在线看片| 亚洲av福利一区| 一级片'在线观看视频| 成人特级av手机在线观看| 男男h啪啪无遮挡| 国产成人aa在线观看| 国产淫语在线视频| 婷婷色麻豆天堂久久| 91久久精品国产一区二区三区| 水蜜桃什么品种好| 内射极品少妇av片p| 精华霜和精华液先用哪个| 国产成人a∨麻豆精品| 最近的中文字幕免费完整| 成人午夜精彩视频在线观看| 国产精品福利在线免费观看| 91精品伊人久久大香线蕉| 黑人巨大精品欧美一区二区蜜桃 | 精品人妻一区二区三区麻豆| 成人毛片60女人毛片免费| 国产男人的电影天堂91| 久久99热6这里只有精品| 欧美+日韩+精品| 国产 一区精品| 人人妻人人添人人爽欧美一区卜| 精品亚洲乱码少妇综合久久| 欧美日韩一区二区视频在线观看视频在线| 精品一区在线观看国产| 亚洲av不卡在线观看| www.av在线官网国产| 日韩伦理黄色片| 18禁在线无遮挡免费观看视频| 久久精品国产a三级三级三级| 黑丝袜美女国产一区| 国产极品粉嫩免费观看在线 | 亚洲国产精品成人久久小说| 国产伦理片在线播放av一区| 久久国产乱子免费精品| 一本久久精品| 少妇被粗大的猛进出69影院 | 国产毛片在线视频| 日本wwww免费看| 国模一区二区三区四区视频| 人妻 亚洲 视频| 亚洲精品亚洲一区二区| 免费观看无遮挡的男女| 久久久a久久爽久久v久久| 国产成人精品久久久久久| 少妇的逼好多水| 亚洲欧美日韩另类电影网站| 黄色配什么色好看| 中国三级夫妇交换| 一级黄片播放器| 欧美+日韩+精品| 国产欧美日韩综合在线一区二区 | 国产欧美亚洲国产| 99久久精品热视频| 国产美女午夜福利| 免费大片黄手机在线观看| 国产伦精品一区二区三区视频9| 国模一区二区三区四区视频| 国产成人午夜福利电影在线观看| 亚洲国产av新网站| 51国产日韩欧美| 国产毛片在线视频| 草草在线视频免费看| 午夜老司机福利剧场| 国产成人精品一,二区| 天天操日日干夜夜撸| 久久久久人妻精品一区果冻| 久久免费观看电影| 久久人人爽av亚洲精品天堂| 国产色爽女视频免费观看| 天天操日日干夜夜撸| 成人18禁高潮啪啪吃奶动态图 | 免费黄频网站在线观看国产| 亚洲精品久久久久久婷婷小说| 校园人妻丝袜中文字幕| 欧美97在线视频| 亚洲av中文av极速乱| 日本爱情动作片www.在线观看| 欧美成人午夜免费资源| 永久网站在线| 黄色怎么调成土黄色| 老熟女久久久| 亚洲av中文av极速乱| 国产高清有码在线观看视频| 精品一区二区三卡| 只有这里有精品99| 最近手机中文字幕大全| 日韩一区二区视频免费看| 久久99蜜桃精品久久| 美女福利国产在线| 国产精品国产av在线观看| 如日韩欧美国产精品一区二区三区 | 视频区图区小说| 婷婷色综合www| 午夜视频国产福利| 国产一区二区在线观看日韩| a 毛片基地| 久久av网站| 日韩不卡一区二区三区视频在线| 一级黄片播放器| 亚洲中文av在线| 五月玫瑰六月丁香| 亚洲精品国产av蜜桃| 国产熟女欧美一区二区| 丰满迷人的少妇在线观看| 国产91av在线免费观看| 国产成人精品无人区| 亚洲国产精品一区三区| 国产日韩欧美在线精品| av免费观看日本| 最后的刺客免费高清国语| 久久精品熟女亚洲av麻豆精品| 街头女战士在线观看网站| 午夜福利视频精品| 国产又色又爽无遮挡免| 色吧在线观看| 特大巨黑吊av在线直播| 男人舔奶头视频| 久久 成人 亚洲| 夜夜爽夜夜爽视频| 国产欧美日韩综合在线一区二区 | 久久久久久久久久人人人人人人| 欧美另类一区| 亚洲精品一区蜜桃| 中文欧美无线码| 好男人视频免费观看在线| 99久久精品一区二区三区| 亚洲精品日韩av片在线观看| 狂野欧美激情性bbbbbb| 国产精品一区二区三区四区免费观看| 亚洲内射少妇av| 亚洲av不卡在线观看| 永久网站在线| 美女福利国产在线| 伦精品一区二区三区| 精品久久久久久电影网| 99久久精品国产国产毛片| 国产亚洲av片在线观看秒播厂| 大陆偷拍与自拍| 美女主播在线视频| 国产成人免费无遮挡视频| 极品教师在线视频| 欧美老熟妇乱子伦牲交| 黄色日韩在线| 男人舔奶头视频| 国产成人91sexporn| av在线播放精品| 蜜桃久久精品国产亚洲av| 欧美变态另类bdsm刘玥| 国产av码专区亚洲av| 亚洲真实伦在线观看| 中文精品一卡2卡3卡4更新| 国产免费视频播放在线视频| 久久99热6这里只有精品| 亚洲精品一二三| 亚洲av中文av极速乱| 欧美另类一区| 水蜜桃什么品种好| 欧美日韩综合久久久久久| 成人特级av手机在线观看| 久久久久久久久久久丰满| 91午夜精品亚洲一区二区三区| 大片免费播放器 马上看| 午夜免费男女啪啪视频观看| 国产精品人妻久久久影院| 欧美成人精品欧美一级黄| 日韩欧美 国产精品| 日日摸夜夜添夜夜爱| 在线观看免费视频网站a站| 日韩不卡一区二区三区视频在线| 婷婷色av中文字幕| 男女啪啪激烈高潮av片| 亚洲欧美清纯卡通| 日韩熟女老妇一区二区性免费视频| 成人特级av手机在线观看| 十八禁高潮呻吟视频 | 自线自在国产av| 男的添女的下面高潮视频| 黄色毛片三级朝国网站 | 91久久精品国产一区二区成人| 两个人免费观看高清视频 | 午夜视频国产福利| 国产精品国产三级国产av玫瑰| 精品人妻熟女av久视频| 日韩欧美一区视频在线观看 | 成人无遮挡网站| 大又大粗又爽又黄少妇毛片口| 日本黄大片高清| 男女国产视频网站| 亚洲精品国产av蜜桃| av福利片在线观看| 人人妻人人澡人人看| 欧美另类一区| 国产 精品1| 久久久国产欧美日韩av| 久久久久久久久久久免费av| 免费观看a级毛片全部| 激情五月婷婷亚洲| 女性被躁到高潮视频| tube8黄色片| 夫妻性生交免费视频一级片| 欧美日韩视频精品一区| 97在线视频观看| 2021少妇久久久久久久久久久| 国产乱人偷精品视频| 你懂的网址亚洲精品在线观看| 免费不卡的大黄色大毛片视频在线观看| 老司机影院成人| 成人亚洲精品一区在线观看| 国精品久久久久久国模美| 亚洲人成网站在线观看播放| 免费久久久久久久精品成人欧美视频 | 菩萨蛮人人尽说江南好唐韦庄| 人妻系列 视频| 久久人妻熟女aⅴ| 国产免费视频播放在线视频| 久久鲁丝午夜福利片| 激情五月婷婷亚洲| 男女无遮挡免费网站观看| 国产成人91sexporn| 免费看日本二区| 亚洲精品aⅴ在线观看| av国产久精品久网站免费入址| 麻豆成人av视频| 国产精品免费大片| 国产淫语在线视频| 啦啦啦在线观看免费高清www| 国产爽快片一区二区三区| 国产精品福利在线免费观看| 麻豆精品久久久久久蜜桃| 国产免费一级a男人的天堂| av卡一久久| 麻豆成人av视频| 成人特级av手机在线观看| 免费少妇av软件| 国内揄拍国产精品人妻在线| 国产精品蜜桃在线观看| 日本午夜av视频| 国产成人aa在线观看| 制服丝袜香蕉在线| 免费观看的影片在线观看| av有码第一页| 高清不卡的av网站| 国产在线免费精品| 亚洲精品久久午夜乱码| 国产成人午夜福利电影在线观看| 男女国产视频网站| 18+在线观看网站| 啦啦啦在线观看免费高清www| 精品人妻一区二区三区麻豆| 美女国产视频在线观看| 日日啪夜夜爽| 精品一区二区免费观看| 亚洲人成网站在线播| 少妇人妻一区二区三区视频| 秋霞在线观看毛片| 汤姆久久久久久久影院中文字幕| 人妻系列 视频| 精品人妻偷拍中文字幕| 下体分泌物呈黄色| 午夜福利网站1000一区二区三区| 亚洲av免费高清在线观看| 亚洲成色77777| 高清不卡的av网站| 大香蕉久久网| 中文字幕av电影在线播放| 日日摸夜夜添夜夜添av毛片| 亚洲情色 制服丝袜| 人人妻人人看人人澡| 一边亲一边摸免费视频| 国产伦精品一区二区三区四那| 在线观看三级黄色| 成人综合一区亚洲| 国产国拍精品亚洲av在线观看| 一级爰片在线观看| 日韩三级伦理在线观看| 国产高清三级在线| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品国产精品| 日日撸夜夜添| 涩涩av久久男人的天堂| 我的女老师完整版在线观看| 老女人水多毛片| 99热网站在线观看| 少妇精品久久久久久久| 各种免费的搞黄视频| 极品少妇高潮喷水抽搐| 亚洲一级一片aⅴ在线观看| av黄色大香蕉| 国产成人91sexporn| 亚洲精品乱码久久久v下载方式| 如日韩欧美国产精品一区二区三区 | 2021少妇久久久久久久久久久| 99视频精品全部免费 在线| 91精品伊人久久大香线蕉| 亚洲中文av在线| 免费观看无遮挡的男女| 成人亚洲精品一区在线观看| 久久午夜福利片| 伦理电影大哥的女人| 久久6这里有精品| 日本黄色日本黄色录像| 免费观看无遮挡的男女| 国产成人精品福利久久| 91精品国产国语对白视频| 十八禁网站网址无遮挡 | 免费观看在线日韩| 国产精品嫩草影院av在线观看| 精品国产国语对白av| 久久久国产精品麻豆| 午夜福利视频精品| 美女脱内裤让男人舔精品视频| av国产久精品久网站免费入址| 国产高清不卡午夜福利| 美女cb高潮喷水在线观看| 夫妻性生交免费视频一级片| 美女cb高潮喷水在线观看| 亚洲电影在线观看av| 国语对白做爰xxxⅹ性视频网站| av天堂中文字幕网| 美女cb高潮喷水在线观看| 日本-黄色视频高清免费观看| 亚洲欧美清纯卡通| 亚洲精品aⅴ在线观看| 下体分泌物呈黄色| 国产精品国产三级专区第一集| 三级国产精品片| 街头女战士在线观看网站| 亚洲av综合色区一区| 哪个播放器可以免费观看大片| 亚洲精品中文字幕在线视频 | 韩国av在线不卡| 亚洲精品自拍成人| 久久 成人 亚洲| 亚洲精品一二三| 美女xxoo啪啪120秒动态图| 国产成人精品一,二区| 日韩中字成人| 亚洲精品国产av蜜桃| 精品国产一区二区久久| 中文字幕亚洲精品专区| 菩萨蛮人人尽说江南好唐韦庄| 极品少妇高潮喷水抽搐| 国产成人精品婷婷| 午夜免费男女啪啪视频观看| 天美传媒精品一区二区| 亚洲成色77777| 国语对白做爰xxxⅹ性视频网站| 肉色欧美久久久久久久蜜桃| 成年人午夜在线观看视频| 精品人妻偷拍中文字幕| 啦啦啦啦在线视频资源| 精品国产一区二区久久| 亚洲第一av免费看| 日本猛色少妇xxxxx猛交久久| 在线观看免费视频网站a站| 亚洲精品国产成人久久av| 熟女av电影| 视频中文字幕在线观看| 免费av不卡在线播放| 一本久久精品| 天天躁夜夜躁狠狠久久av| 国产乱来视频区| 日本色播在线视频| 老熟女久久久| 国产亚洲91精品色在线| 久久毛片免费看一区二区三区| 我要看日韩黄色一级片| av有码第一页| 国产亚洲精品久久久com| 91久久精品电影网| 欧美3d第一页| 精品一区二区免费观看| 精品视频人人做人人爽| 国产日韩一区二区三区精品不卡 | 夫妻性生交免费视频一级片| 老司机影院成人| 亚洲精品乱码久久久久久按摩| 日韩一区二区三区影片| 免费播放大片免费观看视频在线观看| 国产亚洲av片在线观看秒播厂| 大片电影免费在线观看免费| 中国美白少妇内射xxxbb| 亚洲色图综合在线观看| 美女脱内裤让男人舔精品视频| 老司机影院成人| 男的添女的下面高潮视频| 老司机影院成人| 制服丝袜香蕉在线| 日本av免费视频播放| 亚洲精品国产成人久久av| 国产免费又黄又爽又色| 亚洲av.av天堂| 十八禁高潮呻吟视频 | 男人狂女人下面高潮的视频| 欧美精品一区二区免费开放| 狂野欧美白嫩少妇大欣赏| 精品熟女少妇av免费看| 好男人视频免费观看在线| 色94色欧美一区二区| 亚洲性久久影院| 97在线视频观看| 久久精品熟女亚洲av麻豆精品| 亚州av有码| 又爽又黄a免费视频| 成人国产麻豆网| 特大巨黑吊av在线直播| 日本色播在线视频| 在线观看www视频免费| 亚洲av中文av极速乱| 久久久久国产精品人妻一区二区| 寂寞人妻少妇视频99o| 亚洲,一卡二卡三卡| 欧美人与善性xxx| 偷拍熟女少妇极品色| 国产又色又爽无遮挡免| 国产亚洲午夜精品一区二区久久| 18+在线观看网站| 日日摸夜夜添夜夜爱| 亚洲第一区二区三区不卡| 一级黄片播放器| 久久久久精品久久久久真实原创| 日韩视频在线欧美| 色视频在线一区二区三区| 色网站视频免费| 国产精品99久久99久久久不卡 | 精品少妇久久久久久888优播| 黑丝袜美女国产一区| 乱系列少妇在线播放| 哪个播放器可以免费观看大片| 在线观看人妻少妇| 中国国产av一级| 波野结衣二区三区在线| 狂野欧美激情性bbbbbb| 又大又黄又爽视频免费| 美女中出高潮动态图| 久久综合国产亚洲精品| 高清午夜精品一区二区三区| 久久人人爽人人片av| 搡女人真爽免费视频火全软件| 在线观看一区二区三区激情| 久久影院123| 女性生殖器流出的白浆| 在线 av 中文字幕| 精华霜和精华液先用哪个| 久久久久久久亚洲中文字幕| 天堂8中文在线网| 精品国产露脸久久av麻豆| 亚洲欧美一区二区三区国产| 日韩不卡一区二区三区视频在线| 精品一品国产午夜福利视频| 高清不卡的av网站| av网站免费在线观看视频| 欧美成人精品欧美一级黄| 大码成人一级视频| 99国产精品免费福利视频| 亚洲va在线va天堂va国产| 国产极品天堂在线| 精品亚洲成a人片在线观看| 狂野欧美激情性xxxx在线观看| 性高湖久久久久久久久免费观看| 久久久久精品性色| 免费播放大片免费观看视频在线观看| 久久这里有精品视频免费| 免费观看无遮挡的男女| 久久久久久久久久久久大奶| 特大巨黑吊av在线直播| 国产男女内射视频| 成人免费观看视频高清| 51国产日韩欧美| 亚洲怡红院男人天堂| h视频一区二区三区| 一二三四中文在线观看免费高清| 成人毛片a级毛片在线播放| 嘟嘟电影网在线观看| 黑人高潮一二区| 欧美+日韩+精品| 热99国产精品久久久久久7| 亚洲熟女精品中文字幕| 中文字幕久久专区| 亚洲,一卡二卡三卡| 天天操日日干夜夜撸| 久久久久久久国产电影| 九色成人免费人妻av| 欧美激情极品国产一区二区三区 | 亚洲国产最新在线播放| 久久久午夜欧美精品| 精品亚洲乱码少妇综合久久| 在线观看国产h片| 国产精品久久久久久久电影| 成人二区视频| 2018国产大陆天天弄谢| 九草在线视频观看| 最近2019中文字幕mv第一页| 精品国产一区二区久久| 亚洲欧美成人精品一区二区| 中文字幕精品免费在线观看视频 | 亚洲av中文av极速乱| 97在线人人人人妻| 99九九线精品视频在线观看视频| 国产精品久久久久久久电影| 一本—道久久a久久精品蜜桃钙片| 中文资源天堂在线| 亚洲综合精品二区| 街头女战士在线观看网站| 亚洲第一av免费看| 日日摸夜夜添夜夜爱| 91午夜精品亚洲一区二区三区| 又爽又黄a免费视频| 久久人人爽av亚洲精品天堂| 大话2 男鬼变身卡| 中文字幕制服av| 欧美少妇被猛烈插入视频| 久久女婷五月综合色啪小说| 免费在线观看成人毛片| 国产一区亚洲一区在线观看| 精品酒店卫生间| 久久精品夜色国产| 十分钟在线观看高清视频www | 亚洲欧洲精品一区二区精品久久久 | 99国产精品免费福利视频| 国产老妇伦熟女老妇高清| 人人妻人人澡人人看| 少妇被粗大的猛进出69影院 | 各种免费的搞黄视频| 久久鲁丝午夜福利片| 51国产日韩欧美| 中国三级夫妇交换| 免费大片黄手机在线观看| 97精品久久久久久久久久精品| 日日撸夜夜添| 欧美日本中文国产一区发布| 51国产日韩欧美| 18禁裸乳无遮挡动漫免费视频| 日韩,欧美,国产一区二区三区| 麻豆成人午夜福利视频| 亚洲伊人久久精品综合| 国产精品不卡视频一区二区| 最后的刺客免费高清国语| 亚洲天堂av无毛| 欧美变态另类bdsm刘玥| 高清欧美精品videossex| 观看av在线不卡| 赤兔流量卡办理| 2021少妇久久久久久久久久久| 国产极品粉嫩免费观看在线 | 人体艺术视频欧美日本| 婷婷色av中文字幕| 久久久精品94久久精品| 熟女电影av网| 在线亚洲精品国产二区图片欧美 | 亚洲成人一二三区av| 国产精品三级大全| 一级毛片黄色毛片免费观看视频| 香蕉精品网在线| 91久久精品国产一区二区三区| 欧美日韩一区二区视频在线观看视频在线|