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

    基于SPH法的正弦形消力池底板水躍現(xiàn)象數(shù)值模擬

    2019-03-26 08:56:16薄夫萍吳海濤張?jiān)圃?/span>顧聲龍
    水利水電科技進(jìn)展 2019年2期
    關(guān)鍵詞:消能率水躍消力池

    薄夫萍,吳海濤,王 星,張?jiān)圃?,顧聲?/p>

    (青海大學(xué)水利電力學(xué)院,青海 西寧 810016)

    水躍是底流消能工中常用的消能形式。經(jīng)過水躍這一局部水力現(xiàn)象,水流流態(tài)發(fā)生急劇變化,水躍區(qū)段內(nèi)底部主流區(qū)和上部旋滾區(qū)相互摻混、碰撞、摩擦,下泄水流消耗大部分動(dòng)能從而使流速銳減、流態(tài)趨于平穩(wěn),以達(dá)到消能效果。但當(dāng)消能率較低時(shí),下泄水流經(jīng)消力池后仍攜帶大量余能,致使下游河床遭受嚴(yán)重沖刷甚至造成壩體損壞[1]。因此,提高消力池的消能率是目前水工結(jié)構(gòu)設(shè)計(jì)中迫切需要解決的問題。

    正弦形消力池底板最初由Ead等[2]提出,其本質(zhì)是將光滑的消力池底板改為正弦形以增加壁面粗糙度,使水流之間的紊動(dòng)更為劇烈,水躍特性急劇改變,從而增加消能率。目前對(duì)此種消力池的研究方法主要有試驗(yàn)研究、理論分析、數(shù)值模擬。Hughes等[3]通過對(duì)消力池底板加設(shè)不同礫石的試驗(yàn)研究表明,隨著底板粗糙度的增加,水躍躍長明顯縮短,消能率有所提高;Tokyay[4]、Izadjoo等[5]、Abbaspour等[6]、Elsebaie等[7]、Samadi等[8]分別通過大量試驗(yàn)研究了正弦形、梯形、三角形底板上的水躍,指明共軛水深是與弗勞德數(shù)(Fr)有關(guān)的函數(shù),且正弦形底板在消能方面具有很大潛力。傅銘煥等[9]、張志昌等[10-11]結(jié)合斷面流速分布對(duì)數(shù)規(guī)律及動(dòng)量方程給出了共軛水深、水躍長度等理論計(jì)算方法,并指出壁面阻力對(duì)波狀底板水躍特性有很大的影響,是解決實(shí)際問題中不可忽略的重要因素。由于試驗(yàn)研究受到儀器、測(cè)量方法等因素的限制,很難準(zhǔn)確測(cè)量復(fù)雜水力現(xiàn)象中急劇變化的參數(shù);而理論分析需對(duì)研究問題做大量假設(shè)以簡(jiǎn)化控制方程便于求解,但與實(shí)際現(xiàn)象偏離且求解過程繁瑣,且參數(shù)選取主要基于經(jīng)驗(yàn)公式。數(shù)值模擬則可以彌補(bǔ)試驗(yàn)研究及理論分析的不足,記錄現(xiàn)象演變的全過程,為所研究問題提供完整、具體的數(shù)據(jù),便于細(xì)部結(jié)構(gòu)的觀察,有利于揭示現(xiàn)象變化的本質(zhì)。目前研究正弦形底板水躍現(xiàn)象采用的數(shù)值模擬方法主要為網(wǎng)格方法:魏文禮等[12-13]、程香菊等[14]、Abbaspour等[15]采用網(wǎng)格法建立不同形式的湍流模型,對(duì)波狀底板上的水躍進(jìn)行數(shù)值模擬,得到水躍的演變過程及近壁面處存在的小渦分布。但網(wǎng)格法由于受到網(wǎng)格節(jié)點(diǎn)的限制,在模擬自由表面、運(yùn)動(dòng)交界面、大變形等引起網(wǎng)格嚴(yán)重變形的多尺度問題時(shí)精度和準(zhǔn)確性差強(qiáng)人意。而水躍現(xiàn)象中水流紊動(dòng)劇烈并伴隨空氣不斷地?fù)饺隱16],自由表面存在旋滾,因而水面波動(dòng)大且水體破碎嚴(yán)重,因此考慮采用新的數(shù)值模擬方法即光滑質(zhì)點(diǎn)水動(dòng)力學(xué)法(SPH方法)進(jìn)行研究。SPH方法是純拉格朗日形式的無網(wǎng)格粒子法之一,對(duì)需要計(jì)算的問題域采用任意分布的粒子框架進(jìn)行描述,在處理大變形等問題上具有獨(dú)特優(yōu)勢(shì),目前已廣泛地應(yīng)用于滑坡涌浪[17]、多相流[18]、潰壩[19]、溢洪道[20]等工程流體動(dòng)力學(xué)問題的研究,并取得了重要的成果。

    結(jié)合水躍的特點(diǎn)及SPH方法的優(yōu)勢(shì),筆者將使用SPHYSICS_2D軟件對(duì)正弦形消力池底板上水躍現(xiàn)象建模計(jì)算,以分析其水面線、流速分布、躍長、共軛水深、消能率等水躍特性相較于光滑底板的變化趨勢(shì)。

    1 經(jīng)SPH方法離散后的N-S方程

    流體動(dòng)力學(xué)中通常使用邊界條件、初始狀態(tài)結(jié)合質(zhì)量、動(dòng)量、能量三大守恒定律建立控制方程,以確定流體系統(tǒng)的運(yùn)動(dòng)情況。對(duì)無窮小流體單元采用拉格朗日描述法建立控制方程可得偏微分形式下的連續(xù)性方程、動(dòng)量方程、能量方程,即N-S方程:

    (1)

    (2)

    (3)

    式中:α、β表示坐標(biāo)方向;ρ為流體密度;t為時(shí)間;V為流體速度向量;X為流體位置向量;g為重力加速度;e為流體能量;σαβ為總應(yīng)力張量;P各向同性壓力;θ為黏性項(xiàng)。

    SPH方法中將任意函數(shù)改寫為SPH形式需要通過核近似、粒子近似2個(gè)關(guān)鍵步驟。經(jīng)核近似后得場(chǎng)函數(shù)積分表示(式(4)),經(jīng)粒子近似后得函數(shù)的離散形式(式(5)):

    (4)

    (5)

    其中Wij=W(ri-rj,h)

    式中:A(r)為任意函數(shù);W(r-r′,h)為光滑函數(shù);h為光滑長度,用以控制計(jì)算域;r和r′代表粒子位置;N為粒子總數(shù);mj和ρj分別為j粒子質(zhì)量和密度;Wij為i、j粒子間的光滑函數(shù)。

    SPH方法中采用亞粒子模型(SPS)建立湍流模型,該模型由大渦模型(LES)和亞網(wǎng)格尺寸模型(SGS)組合而成。對(duì)黏性項(xiàng)θ的處理采用SPS方法,可得動(dòng)量方程(式(6));在SPH方法中流體被認(rèn)為是弱可壓縮性的,因此為控制流體的壓縮性引入了狀態(tài)方程(7),通過改變聲速調(diào)節(jié)數(shù)值模型時(shí)間步長以滿足壓縮性調(diào)節(jié)要求。但是聲速也不能過小,至少是流體最大流速的10倍,這樣才能確保粒子密度變化率小于1%,控制流體為弱可壓縮狀態(tài)。

    (6)

    (7)

    式中:ρ0為預(yù)設(shè)密度,取值為1 000 kg/m3;c0為預(yù)設(shè)密度對(duì)應(yīng)的聲速;?0為層流動(dòng)力黏度,取10-6m2/s;τ為黏性應(yīng)力。

    綜上所述,經(jīng)SPH方法離散后的N-S方程分別為以下形式:

    (8)

    (9)

    (10)

    式中:下標(biāo)a、b表示計(jì)算域內(nèi)任意2個(gè)粒子;mb為b粒子的質(zhì)量;Vab為a、b粒子之間的相對(duì)速度向量;Wab為a、b粒子間核函數(shù)梯度;ψab為經(jīng)SPH方法離散后的黏性項(xiàng)。

    水躍消能率的計(jì)算通常使用躍前斷面能量E1與躍后斷面能量E2之差作為能量損失值ΔE,因此可以得到躍后能量損失為

    (11)

    式中:Z1、Z2分別為躍前斷面、躍后斷面水位高度;α1、α2分別為躍前斷面、躍后斷面的動(dòng)能修正系數(shù);U1、U2分別為躍前斷面、躍后斷面平均流速;P1、P2分別為躍前斷面、躍后斷面壓強(qiáng)。

    當(dāng)渠底水平時(shí),根據(jù)總流的動(dòng)量方程可以推導(dǎo)出單寬流量q與躍前水深y1、躍后水深y2之間的關(guān)系,從而得到光滑底板上能量損失ΔE及消能率η:

    (12)

    (13)

    2 水躍現(xiàn)象數(shù)值模型

    為驗(yàn)證SPH方法研究正弦形底板上水躍現(xiàn)象的可行性及適用性,選取文獻(xiàn)[2]中的試驗(yàn)數(shù)據(jù)與SPH模擬值作對(duì)比。文獻(xiàn)[2]中試驗(yàn)?zāi)P褪褂玫木匦嗡鄄捎糜袡C(jī)玻璃制作,寬0.446 m、高0.6 m、長7.6 m,并在水槽底部鋪設(shè)鋁制正弦形波狀底板。試驗(yàn)布置如圖1所示,圖中Lj為躍長;Lrj為旋滾區(qū)長度;Um為斷面最大流速;S為波長,T為波幅。研究的波形有2種:S均為68 mm,T分別為13 mm、22 mm。由于正弦形底板凹凸不平,易引發(fā)空蝕問題,Ead等[21]研究發(fā)現(xiàn)將正弦形底板設(shè)置在光滑水平面以下,波峰與光滑底板位于同一高度時(shí),能有效避免壁面凹凸不平引起的空蝕問題。

    圖1 試驗(yàn)布置示意圖[2]

    根據(jù)試驗(yàn)?zāi)P筒捎肧PHYSICS_2D軟件建立相應(yīng)的數(shù)值模型,如圖2所示。模型主要結(jié)構(gòu)包括蓄水池、水槽、尾水收集池。正弦邊界采用公式y(tǒng)=Asin(π/0.034x)建立,當(dāng)參數(shù)A(A為波幅的1/2)分別為0.006 5和0.011時(shí)可得T分別為13 mm、22 mm,S為68 mm的正弦形底板。模型由2種粒子構(gòu)成,分別為可流動(dòng)的水粒子和固定邊界粒子。各工況閘孔自由出流流量調(diào)節(jié)通過蓄水池的水位調(diào)節(jié)實(shí)現(xiàn);流經(jīng)水槽的水流暫存在收集池以防水粒子運(yùn)動(dòng)到計(jì)算域以外,對(duì)數(shù)值計(jì)算產(chǎn)生干擾。建模時(shí)為減小運(yùn)行時(shí)間預(yù)先在整個(gè)水槽段添加一定深度的水域。

    圖2 數(shù)值模型示意圖

    對(duì)文獻(xiàn)[2]中5組工況進(jìn)行數(shù)值模擬,具體參數(shù)如表1所示,表中y2*為光滑壁面躍后水深。數(shù)值運(yùn)算時(shí)選取的主要參數(shù)如下:核函數(shù)為5次樣條函數(shù);黏性項(xiàng)處理采用亞粒子湍流模型;邊界采用動(dòng)力邊界條件。模型粒子總數(shù)為96 476個(gè),其中邊界粒子4 610個(gè),初始粒子間距為5×10-3m,光滑長度為3.68×10-3m,時(shí)間步長為0.01 s,模擬時(shí)長為10 s。

    表1 工況參數(shù)

    3 結(jié)果分析

    3.1 模型驗(yàn)證

    由于本數(shù)值模型只能在初始狀態(tài)添加水域進(jìn)行計(jì)算,故為滿足試驗(yàn)工況所需的恒定入流條件,數(shù)值建模時(shí)在水池中添加一個(gè)推板,即推板模型[22]。通過設(shè)定推板的速度對(duì)閘口附近水域進(jìn)行補(bǔ)給,以保證蓄水池中水位恒定。分析得知各試驗(yàn)工況閘口處流量的SPH模擬值與物理模型試驗(yàn)值間的誤差分別為3.9%、6.6%。由流量對(duì)比(圖3)可得:推板模型能夠在短時(shí)間內(nèi)達(dá)到各工況所需的穩(wěn)定流量且流量值出現(xiàn)波動(dòng)的情況較少,表明推板模型能夠有效解決正弦形底板水躍研究中的流量恒定問題。

    將工況A的SPH水面線模擬值和已有文獻(xiàn)試驗(yàn)值及模擬值作對(duì)比,如圖4所示。由圖4可得各工況水面變化趨勢(shì)基本保持一致,SPH數(shù)值模擬結(jié)果能夠充分展現(xiàn)水面線形態(tài)的變化。

    圖3 閘口處流量對(duì)比

    圖4 工況A水面線對(duì)比

    將5種工況水面線的SPH模擬值與文獻(xiàn)[2]試驗(yàn)值作對(duì)比,如圖5所示,可知兩者基本吻合,僅在某些工況,如工況B、C旋滾區(qū)段水面有一定的偏差??紤]到水躍是一種較為復(fù)雜的流態(tài),在水躍區(qū)段內(nèi)主流區(qū)和旋滾區(qū)之間的水質(zhì)點(diǎn)不斷摻混,伴隨大量的氣泡產(chǎn)生或摻入,導(dǎo)致水面劇烈波動(dòng),部分粒子飛濺,造成試驗(yàn)值存在較大誤差。工況D和E流量相同、波形不同時(shí),水面線變化趨勢(shì)基本相同,水面線基本重合,說明波形的改變對(duì)水面線高度的影響較小,分析其原因可能是由于波狀底板位于水平面下方,沒有伸入流體內(nèi)部;且水面的高度遠(yuǎn)大于波狀底板波幅值,故導(dǎo)致正弦形底板上波形的變化對(duì)水面線影響相對(duì)較小。

    圖5 各工況水面線對(duì)比

    從5種工況水面線的整體情況來看,SPH數(shù)值模擬結(jié)果和試驗(yàn)結(jié)果吻合較好,說明SPH方法能夠較為準(zhǔn)確地反映水躍過程中水面的變化趨勢(shì)及形態(tài),在處理波狀消力池底板上的水躍特性方面具有適用性及可靠性。

    3.2 速度云圖、矢量圖

    由工況A速度分布(圖6)得:水躍段包括速度為正的順流及速度為負(fù)的逆流,其中逆流主要聚集在水流表面和底部的主流處,表面由于部分水流回流形成逆流,底部因正弦形底板的波狀結(jié)構(gòu)形成波內(nèi)旋渦導(dǎo)致流速為負(fù);水躍段存在表面旋滾,且有大量氣體摻入,速度分布波動(dòng)較大,當(dāng)?shù)装鍨檎倚螘r(shí)波內(nèi)形成旋渦,相當(dāng)于增加了邊界層的厚度和粗糙度,近壁面及表面流速較小,中間層流速較大;當(dāng)水流流經(jīng)明渠段時(shí)受邊界層影響,表面流速最大,近壁面處最小。

    圖6 工況A速度云圖

    圖7 工況A速度矢量

    圖8 工況B速度矢量

    結(jié)合工況A、B的速度矢量(圖7、圖8)可得:水躍段內(nèi)有明顯水流破碎和摻氣現(xiàn)象;在閘口開度相同時(shí),隨著Fr的增大,水流破碎、摻氣現(xiàn)象加劇并且存在明顯渦結(jié)構(gòu),分析其原因主要是水躍段流速分布不均勻,流層間的水質(zhì)點(diǎn)有相對(duì)運(yùn)動(dòng),導(dǎo)致流層間產(chǎn)生的內(nèi)摩擦切應(yīng)力在水流表面為順流方向,而在水流底部為逆流方向,故形成了順時(shí)針流向的旋滾。對(duì)比兩工況旋滾得:同等條件下旋滾傾斜程度隨Fr的增大逐漸增加,即形成旋滾的長度逐漸縮短,旋滾內(nèi)的渦結(jié)構(gòu)愈發(fā)明顯,渦團(tuán)尺寸逐漸增加,這將意味著旋滾的影響范圍逐漸變大。

    3.3 躍長

    水躍段水流劇烈翻滾摻混使水面形態(tài)急劇變化,旋滾的尾部水流不斷匯入主流并被逐漸淹沒,直到水流表面不會(huì)因?yàn)樗鞯姆瓭L而產(chǎn)生較大波動(dòng)時(shí)流態(tài)演變?yōu)槊髑鲃?dòng)。根據(jù)圖7可得,水躍前部1—1斷面至形成旋滾的2—2斷面之間的距離為旋滾區(qū)長度Lrj;1—1斷面至水流經(jīng)過旋滾又回到較為穩(wěn)定的自由水面3—3斷面之間的距離為Lj。

    將5種工況Lrj、Lj模擬值與試驗(yàn)值及光滑底板躍長試驗(yàn)值Lj*對(duì)比得出,Lrj、Lj模擬值相比對(duì)應(yīng)試驗(yàn)值的誤差分別在7%、8%以內(nèi);Lj試驗(yàn)值較光滑底板減小30%左右,表明正弦形底板上的躍長較光滑底板明顯縮短,若應(yīng)用在實(shí)際工程中能夠有效縮短消力池長度。結(jié)合表1對(duì)比工況A、B、C的Lrj、Lj如圖9所示,表明試驗(yàn)值和模擬值變化趨勢(shì)是一致的;且同等條件下,Lrj、Lj均隨著Fr的遞增呈遞增趨勢(shì)。

    圖9 躍長、旋滾區(qū)長度對(duì)比

    3.4 躍后水深

    圖10 共軛水深、躍后水深對(duì)比

    將圖7中水面變化較穩(wěn)定的3—3斷面處水深值作為躍后水深y2。為比較共軛水深與躍后水深的關(guān)系,將A、B、C這3種工況的共軛水深值統(tǒng)一縮小至原來的1/10與其躍后水深y2繪制在圖10中,對(duì)比y2模擬值與試驗(yàn)值及光滑底板躍后水深y2*得:共軛水深、躍后水深的變化趨勢(shì)相同;y2模擬值與試驗(yàn)值之間的誤差在15%以內(nèi),y2試驗(yàn)值較y2*減少30%左右;比較工況A、B、C得同等條件下y2、y2*隨Fr的增加而增加;結(jié)合表1中工況D、C可發(fā)現(xiàn)波形對(duì)躍后水深影響較小。采用無量綱分析法探究躍后水深、共軛水深與躍長間的相互關(guān)系。結(jié)果表明:y2和Lj變化趨勢(shì)呈正比關(guān)系;當(dāng)共軛水深呈遞增趨勢(shì)時(shí),躍長呈減小趨勢(shì),兩者呈反比關(guān)系,表明波狀底板雖然能夠減小躍長及躍后水深,但是不能同時(shí)減小共軛水深和躍長。

    3.5 消能率

    水躍消能是泄水建筑物中常用的消能形式之一,高速水流經(jīng)過消力池后形成水躍,躍后流速的急劇衰減消耗了上游水流的大部分能量,從而起到消能的目的。因此,在工程造價(jià)預(yù)算范圍內(nèi)對(duì)消能工進(jìn)行合理設(shè)計(jì),以盡可能地提高消能率。

    圖11 流速分析

    將5組工況正弦形底板消能率模擬值、試驗(yàn)值及光滑底板消能率進(jìn)行對(duì)比,見表2,可知當(dāng)Fr=4~8時(shí),光滑消力池底板的消能率范圍為39.36%~66.12%,正弦形消力池底板的消能率范圍為47.08%~73.12%,因此正弦形底板較光滑底板消能率明顯提高了10%左右,且5種工況的消能率均在47%以上。由表2還可得知,消能率試驗(yàn)值和模擬值吻合度較高,變化規(guī)律相似。在所研究的Fr范圍里,光滑底板和正弦形底板的消能率均隨Fr的增加而增加;對(duì)比工況A、B、C可得,當(dāng)波形相同時(shí)消能率隨Fr的增加而增加;對(duì)比D、E兩種工況可得當(dāng)Fr相同時(shí),正弦波形不同時(shí)消能率較為接近,說明其他條件相同時(shí),正弦波形對(duì)消能率的影響不明顯,正弦形底板消能效果明顯。

    表2 各工況不同F(xiàn)r時(shí)的消能率

    4 結(jié) 論

    a. 推板模型能夠獲取試驗(yàn)所需的恒定流量條件,使模擬流量誤差在7%以內(nèi);流量能在較短的時(shí)間內(nèi)達(dá)到穩(wěn)定且很少出現(xiàn)數(shù)值波動(dòng)。

    b. 數(shù)值模擬水面線與試驗(yàn)值基本吻合,僅存在較小波動(dòng),驗(yàn)證了SPH方法研究正弦形底板水躍現(xiàn)象的可行性及適用性。

    c. 水躍段流速分布不均勻且波動(dòng)較大;相同條件下,隨著Fr的增加摻氣量及水面破碎現(xiàn)象越來越明顯,旋滾的影響范圍逐漸變大。

    d. 正弦形底板躍長、躍后水深較光滑底板均較小30%左右,表明波狀底板能夠有效降低躍后水深,減小消力池護(hù)坦高度;對(duì)比模擬工況發(fā)現(xiàn)水面線、躍長受波形影響較??;躍長、躍后水深的變化趨勢(shì)一致,但共軛水深與躍長變化趨勢(shì)相反,故不能同時(shí)減小躍長和共軛水深。

    e. 所模擬工況斷面平均流速、最大流速分布規(guī)律具有相似性,即在水躍段存在較大波動(dòng)、衰減,明渠段基本保持穩(wěn)定;正弦形底板消能率較光滑底板提高10%左右,且5種工況消能率均在47%以上。

    猜你喜歡
    消能率水躍消力池
    套筒閥消能與氣蝕研究及結(jié)構(gòu)改進(jìn)
    那吉航運(yùn)樞紐大壩消力池?fù)p毀成因和搶修方案分析
    一種新型消能結(jié)構(gòu)水力試驗(yàn)研究
    突然擴(kuò)散水躍方程的改進(jìn)與比較
    漸擴(kuò)式消力池水躍長度的分析與計(jì)算
    U型鋼板樁在干河子壩二級(jí)消力池修復(fù)工程中的應(yīng)用
    某攔河大壩下游消力池結(jié)構(gòu)對(duì)消能防沖的影響
    沖沙閘滲透穩(wěn)定計(jì)算及消力池防沖設(shè)計(jì)
    臺(tái)階式溢洪道純臺(tái)階消能率變化規(guī)律研究
    底流消力池內(nèi)懸柵消能工數(shù)值模擬研究
    国产亚洲精品一区二区www| 男女那种视频在线观看| 看片在线看免费视频| 日本成人三级电影网站| 91麻豆精品激情在线观看国产| 一个人免费在线观看电影 | 亚洲自偷自拍图片 自拍| 国产一级毛片七仙女欲春2| avwww免费| 国产黄色小视频在线观看| 亚洲人成电影免费在线| 看黄色毛片网站| 亚洲一区二区三区色噜噜| 久久精品国产清高在天天线| 午夜福利视频1000在线观看| 淫妇啪啪啪对白视频| 国产精品1区2区在线观看.| 少妇人妻一区二区三区视频| 免费大片18禁| 男女之事视频高清在线观看| 亚洲欧美一区二区三区黑人| 丝袜人妻中文字幕| 亚洲av成人精品一区久久| 精品不卡国产一区二区三区| 在线永久观看黄色视频| 麻豆av在线久日| www.精华液| 欧美日韩瑟瑟在线播放| 国产免费av片在线观看野外av| 欧美zozozo另类| 国产v大片淫在线免费观看| 特级一级黄色大片| 久久亚洲精品不卡| 黄片大片在线免费观看| 99热这里只有是精品50| 久久久久久人人人人人| 亚洲av电影在线进入| 国产美女午夜福利| 国产午夜精品论理片| 国产av在哪里看| 国产一区二区三区在线臀色熟女| 亚洲成a人片在线一区二区| 99热精品在线国产| 不卡一级毛片| 午夜福利免费观看在线| www国产在线视频色| 欧美乱色亚洲激情| 国内揄拍国产精品人妻在线| 桃红色精品国产亚洲av| 久久久精品大字幕| 日韩欧美国产一区二区入口| 最近视频中文字幕2019在线8| 女警被强在线播放| 欧美日韩国产亚洲二区| 亚洲真实伦在线观看| 国产精品,欧美在线| 亚洲成人精品中文字幕电影| 黄片小视频在线播放| 欧美黄色片欧美黄色片| 熟女电影av网| 色哟哟哟哟哟哟| av中文乱码字幕在线| 在线视频色国产色| 亚洲自拍偷在线| 国产精品av视频在线免费观看| 美女黄网站色视频| 90打野战视频偷拍视频| 久久久久久久精品吃奶| 国产精华一区二区三区| 欧美3d第一页| 国产又黄又爽又无遮挡在线| 脱女人内裤的视频| 757午夜福利合集在线观看| 老汉色∧v一级毛片| 国产一区在线观看成人免费| 五月玫瑰六月丁香| 亚洲av第一区精品v没综合| 午夜精品在线福利| 99久久精品国产亚洲精品| 国产97色在线日韩免费| 久久草成人影院| 人妻丰满熟妇av一区二区三区| 欧美乱码精品一区二区三区| 女同久久另类99精品国产91| 日韩人妻高清精品专区| 久久精品aⅴ一区二区三区四区| 国产成人福利小说| 制服丝袜大香蕉在线| 国产亚洲精品久久久com| 桃色一区二区三区在线观看| 九色国产91popny在线| 黑人操中国人逼视频| 曰老女人黄片| 每晚都被弄得嗷嗷叫到高潮| 亚洲成人精品中文字幕电影| 免费电影在线观看免费观看| 免费大片18禁| 亚洲国产色片| 法律面前人人平等表现在哪些方面| 精品久久久久久久久久免费视频| 特级一级黄色大片| 悠悠久久av| 国产主播在线观看一区二区| 国产一区二区激情短视频| 欧美在线一区亚洲| 热99re8久久精品国产| 午夜久久久久精精品| 国产亚洲精品综合一区在线观看| 精品不卡国产一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 天堂动漫精品| 天天躁日日操中文字幕| 午夜激情福利司机影院| 99久久国产精品久久久| 亚洲avbb在线观看| 亚洲午夜理论影院| 精品一区二区三区视频在线观看免费| 美女午夜性视频免费| 女生性感内裤真人,穿戴方法视频| 一个人观看的视频www高清免费观看 | 久久人人精品亚洲av| 哪里可以看免费的av片| 亚洲黑人精品在线| 亚洲人成伊人成综合网2020| 制服人妻中文乱码| 老司机深夜福利视频在线观看| 日本成人三级电影网站| 国产成人系列免费观看| 国产一区二区在线av高清观看| 日韩av在线大香蕉| 中文字幕精品亚洲无线码一区| 中文字幕最新亚洲高清| 亚洲最大成人中文| 久久欧美精品欧美久久欧美| 99热只有精品国产| 亚洲aⅴ乱码一区二区在线播放| 午夜激情福利司机影院| 久久国产精品人妻蜜桃| 欧美精品啪啪一区二区三区| 国产乱人视频| 在线免费观看的www视频| 日韩av在线大香蕉| 欧美丝袜亚洲另类 | 99久久久亚洲精品蜜臀av| 婷婷六月久久综合丁香| 午夜亚洲福利在线播放| bbb黄色大片| 精品久久久久久久久久免费视频| av天堂中文字幕网| 亚洲成人中文字幕在线播放| 女生性感内裤真人,穿戴方法视频| 99riav亚洲国产免费| 欧美另类亚洲清纯唯美| 一个人看视频在线观看www免费 | 国产一区在线观看成人免费| 天天添夜夜摸| 欧美日韩一级在线毛片| 欧美大码av| 亚洲色图av天堂| 日韩成人在线观看一区二区三区| 色综合站精品国产| 午夜精品在线福利| 久久久久精品国产欧美久久久| 757午夜福利合集在线观看| av中文乱码字幕在线| 亚洲中文av在线| 青草久久国产| 最近最新中文字幕大全电影3| 国产乱人伦免费视频| 美女免费视频网站| 99久国产av精品| 久久久久久久久免费视频了| 亚洲国产高清在线一区二区三| 丰满的人妻完整版| 亚洲五月天丁香| 欧美成人免费av一区二区三区| 欧美日本视频| 国产高清视频在线播放一区| 丰满人妻一区二区三区视频av | www.999成人在线观看| 18禁裸乳无遮挡免费网站照片| 国产精品久久电影中文字幕| 不卡av一区二区三区| 99国产极品粉嫩在线观看| 国产成人福利小说| 国产欧美日韩一区二区精品| 国产又黄又爽又无遮挡在线| av天堂中文字幕网| 男女做爰动态图高潮gif福利片| 亚洲av熟女| 狂野欧美激情性xxxx| 亚洲国产日韩欧美精品在线观看 | 午夜福利在线观看免费完整高清在 | 美女大奶头视频| 麻豆国产97在线/欧美| 黄色片一级片一级黄色片| 日韩精品中文字幕看吧| avwww免费| 岛国在线免费视频观看| 首页视频小说图片口味搜索| 特级一级黄色大片| 神马国产精品三级电影在线观看| 亚洲18禁久久av| 欧美成狂野欧美在线观看| 九色国产91popny在线| 国产野战对白在线观看| 亚洲成av人片免费观看| 757午夜福利合集在线观看| 免费高清视频大片| 久久伊人香网站| 好男人电影高清在线观看| 男女视频在线观看网站免费| 成人精品一区二区免费| 美女高潮喷水抽搐中文字幕| 少妇的丰满在线观看| 亚洲精品在线美女| 最近最新免费中文字幕在线| 国产免费男女视频| 中文字幕人妻丝袜一区二区| 国内久久婷婷六月综合欲色啪| 久久欧美精品欧美久久欧美| 男人的好看免费观看在线视频| 91九色精品人成在线观看| 国产私拍福利视频在线观看| e午夜精品久久久久久久| 国产视频一区二区在线看| av在线天堂中文字幕| 久久人妻av系列| 久久人人精品亚洲av| 欧美乱妇无乱码| 757午夜福利合集在线观看| 男女床上黄色一级片免费看| 美女高潮的动态| 日韩中文字幕欧美一区二区| 免费搜索国产男女视频| 无限看片的www在线观看| 亚洲熟妇中文字幕五十中出| 校园春色视频在线观看| 嫩草影视91久久| 18美女黄网站色大片免费观看| 看黄色毛片网站| 91麻豆精品激情在线观看国产| 嫁个100分男人电影在线观看| 国产 一区 欧美 日韩| av福利片在线观看| 国产高清视频在线观看网站| 真人一进一出gif抽搐免费| 麻豆一二三区av精品| 大型黄色视频在线免费观看| 99国产综合亚洲精品| 日韩有码中文字幕| 国产亚洲精品综合一区在线观看| 巨乳人妻的诱惑在线观看| 欧美zozozo另类| 欧美日韩亚洲国产一区二区在线观看| 亚洲一区二区三区不卡视频| 精品久久久久久,| 欧美日韩乱码在线| 国产v大片淫在线免费观看| 欧美精品啪啪一区二区三区| 久久天堂一区二区三区四区| 俺也久久电影网| 成人欧美大片| 国产av在哪里看| 亚洲五月天丁香| 中文字幕久久专区| 国产男靠女视频免费网站| 免费看美女性在线毛片视频| 亚洲成a人片在线一区二区| 国产伦一二天堂av在线观看| 亚洲国产精品999在线| 18禁美女被吸乳视频| av黄色大香蕉| 九九久久精品国产亚洲av麻豆 | 一二三四在线观看免费中文在| 国产精品精品国产色婷婷| 中文字幕久久专区| 国产高清视频在线观看网站| 久久国产精品人妻蜜桃| 午夜免费观看网址| 久久久久精品国产欧美久久久| 欧美av亚洲av综合av国产av| 久久久久性生活片| 91字幕亚洲| 给我免费播放毛片高清在线观看| 老司机在亚洲福利影院| 免费大片18禁| 国产午夜精品论理片| 女生性感内裤真人,穿戴方法视频| 中文在线观看免费www的网站| 美女免费视频网站| 免费看日本二区| 嫩草影视91久久| 真人一进一出gif抽搐免费| 麻豆成人av在线观看| 最近最新中文字幕大全免费视频| 又黄又粗又硬又大视频| 19禁男女啪啪无遮挡网站| 亚洲欧美精品综合久久99| 变态另类成人亚洲欧美熟女| 欧美在线黄色| 欧美日韩乱码在线| 亚洲第一电影网av| 黄色 视频免费看| a级毛片a级免费在线| 亚洲欧美精品综合久久99| 成人特级黄色片久久久久久久| 午夜福利免费观看在线| 成人精品一区二区免费| 欧美一区二区国产精品久久精品| 一级毛片精品| 日韩三级视频一区二区三区| 国产激情欧美一区二区| 精品国产超薄肉色丝袜足j| 五月伊人婷婷丁香| 亚洲五月天丁香| 免费观看精品视频网站| 无遮挡黄片免费观看| 欧美成人一区二区免费高清观看 | 国产免费av片在线观看野外av| 国产美女午夜福利| 白带黄色成豆腐渣| 亚洲熟妇熟女久久| 亚洲一区高清亚洲精品| h日本视频在线播放| 久久久久久国产a免费观看| 精品久久久久久成人av| 12—13女人毛片做爰片一| 国产伦人伦偷精品视频| 黄色丝袜av网址大全| 大型黄色视频在线免费观看| 视频区欧美日本亚洲| 99久国产av精品| 精品一区二区三区视频在线观看免费| 日本免费a在线| 99热这里只有精品一区 | 国产欧美日韩精品亚洲av| 成年女人毛片免费观看观看9| 婷婷精品国产亚洲av| 欧美黑人欧美精品刺激| 国产美女午夜福利| 变态另类成人亚洲欧美熟女| 国产亚洲欧美在线一区二区| www.精华液| 国产高潮美女av| cao死你这个sao货| 国产在线精品亚洲第一网站| 国产av在哪里看| 亚洲精品久久国产高清桃花| 亚洲国产欧美网| 久久99热这里只有精品18| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩东京热| 哪里可以看免费的av片| 成年女人毛片免费观看观看9| 免费在线观看影片大全网站| 99国产精品一区二区蜜桃av| 夜夜躁狠狠躁天天躁| 成人三级黄色视频| 麻豆国产av国片精品| 亚洲午夜精品一区,二区,三区| 19禁男女啪啪无遮挡网站| 色综合欧美亚洲国产小说| 日本黄色视频三级网站网址| 午夜亚洲福利在线播放| 亚洲人成网站高清观看| 五月伊人婷婷丁香| 麻豆国产av国片精品| 女同久久另类99精品国产91| 99国产综合亚洲精品| 亚洲在线自拍视频| 亚洲国产精品999在线| 可以在线观看的亚洲视频| 国产精品乱码一区二三区的特点| 长腿黑丝高跟| 法律面前人人平等表现在哪些方面| av天堂在线播放| 午夜福利在线观看吧| 午夜日韩欧美国产| 99热只有精品国产| 久久精品国产99精品国产亚洲性色| 成年女人看的毛片在线观看| 91av网站免费观看| 精品一区二区三区av网在线观看| av中文乱码字幕在线| 精品欧美国产一区二区三| 曰老女人黄片| 久久久久久国产a免费观看| 欧美国产日韩亚洲一区| 久久久久久九九精品二区国产| 日本免费a在线| 啦啦啦韩国在线观看视频| 18禁黄网站禁片免费观看直播| 男插女下体视频免费在线播放| 特级一级黄色大片| 国产精品久久久久久精品电影| 女警被强在线播放| 小说图片视频综合网站| 欧美性猛交╳xxx乱大交人| 两个人看的免费小视频| 日韩 欧美 亚洲 中文字幕| 无人区码免费观看不卡| 在线观看免费视频日本深夜| 在线观看美女被高潮喷水网站 | 热99在线观看视频| 亚洲成a人片在线一区二区| 欧美日韩黄片免| 两个人视频免费观看高清| 中文资源天堂在线| 99在线人妻在线中文字幕| 欧美国产日韩亚洲一区| 在线视频色国产色| 毛片女人毛片| av天堂中文字幕网| 性色av乱码一区二区三区2| 窝窝影院91人妻| 国产精品,欧美在线| svipshipincom国产片| 最好的美女福利视频网| 午夜精品久久久久久毛片777| 亚洲av五月六月丁香网| 久久久久久大精品| 精品一区二区三区av网在线观看| 久久精品aⅴ一区二区三区四区| 国产单亲对白刺激| 欧美+亚洲+日韩+国产| 久久九九热精品免费| 欧美不卡视频在线免费观看| 黄色视频,在线免费观看| 啦啦啦韩国在线观看视频| 男女做爰动态图高潮gif福利片| 欧美色欧美亚洲另类二区| 国产一区二区三区视频了| 99久久无色码亚洲精品果冻| 中文资源天堂在线| 床上黄色一级片| 最近最新中文字幕大全电影3| 国产一区二区在线观看日韩 | 午夜免费激情av| 国产精品乱码一区二三区的特点| 黄色视频,在线免费观看| 国产精品永久免费网站| 亚洲自拍偷在线| 热99re8久久精品国产| 亚洲一区二区三区色噜噜| 午夜福利在线观看吧| 少妇熟女aⅴ在线视频| 亚洲性夜色夜夜综合| 免费观看精品视频网站| 宅男免费午夜| 欧美日韩亚洲国产一区二区在线观看| 最好的美女福利视频网| 国产美女午夜福利| 国产精品亚洲美女久久久| 久久精品综合一区二区三区| 99热这里只有是精品50| 色综合亚洲欧美另类图片| 精品日产1卡2卡| 人人妻,人人澡人人爽秒播| 亚洲五月婷婷丁香| 黄色日韩在线| 久久久久久久精品吃奶| 特级一级黄色大片| 中文字幕熟女人妻在线| 成人欧美大片| www.999成人在线观看| 久久热在线av| 久久天堂一区二区三区四区| 国产精品久久电影中文字幕| 美女免费视频网站| 色视频www国产| 大型黄色视频在线免费观看| 欧美色视频一区免费| 亚洲国产高清在线一区二区三| 精品久久久久久久末码| 在线十欧美十亚洲十日本专区| 黄色视频,在线免费观看| 久久久久久久久免费视频了| 婷婷六月久久综合丁香| 最新在线观看一区二区三区| 日本黄色视频三级网站网址| 这个男人来自地球电影免费观看| h日本视频在线播放| 91久久精品国产一区二区成人 | 两人在一起打扑克的视频| 亚洲 国产 在线| 99在线视频只有这里精品首页| 亚洲中文字幕一区二区三区有码在线看 | 欧美黄色淫秽网站| 亚洲国产高清在线一区二区三| 久久天躁狠狠躁夜夜2o2o| 国产精品国产高清国产av| а√天堂www在线а√下载| 国产免费男女视频| 91av网站免费观看| 国产免费男女视频| 午夜亚洲福利在线播放| 搞女人的毛片| 日本a在线网址| 亚洲自偷自拍图片 自拍| 国产伦人伦偷精品视频| 欧美绝顶高潮抽搐喷水| 搡老妇女老女人老熟妇| 在线十欧美十亚洲十日本专区| 精品免费久久久久久久清纯| 免费在线观看视频国产中文字幕亚洲| 免费大片18禁| 国产免费男女视频| 51午夜福利影视在线观看| 搞女人的毛片| 欧美成人性av电影在线观看| 又黄又爽又免费观看的视频| 亚洲欧美一区二区三区黑人| 日本熟妇午夜| 国产成人精品久久二区二区91| 99国产精品一区二区三区| 亚洲av第一区精品v没综合| 日本熟妇午夜| 亚洲国产精品成人综合色| 精品福利观看| 精品久久久久久久末码| 黑人巨大精品欧美一区二区mp4| 亚洲精品久久国产高清桃花| 中文亚洲av片在线观看爽| 男插女下体视频免费在线播放| 成年女人看的毛片在线观看| 久久午夜亚洲精品久久| 悠悠久久av| 可以在线观看的亚洲视频| 国产三级中文精品| e午夜精品久久久久久久| 亚洲欧美日韩无卡精品| 91麻豆av在线| 久久亚洲精品不卡| 黄色丝袜av网址大全| 黄片小视频在线播放| 最近最新中文字幕大全电影3| 欧美黄色片欧美黄色片| 国内精品久久久久精免费| 亚洲精品粉嫩美女一区| 麻豆一二三区av精品| 亚洲狠狠婷婷综合久久图片| 好看av亚洲va欧美ⅴa在| 亚洲精品中文字幕一二三四区| 国产毛片a区久久久久| 国产精品久久久久久人妻精品电影| 网址你懂的国产日韩在线| 18美女黄网站色大片免费观看| 变态另类成人亚洲欧美熟女| 俺也久久电影网| av在线蜜桃| 三级毛片av免费| 亚洲美女黄片视频| 老司机福利观看| 热99在线观看视频| 最新美女视频免费是黄的| 长腿黑丝高跟| 人人妻人人看人人澡| 在线观看美女被高潮喷水网站 | 日韩欧美免费精品| 亚洲va日本ⅴa欧美va伊人久久| 狠狠狠狠99中文字幕| 婷婷精品国产亚洲av| 男人的好看免费观看在线视频| 人妻夜夜爽99麻豆av| 在线免费观看不下载黄p国产 | 国产精品精品国产色婷婷| 不卡av一区二区三区| 最近最新免费中文字幕在线| 夜夜看夜夜爽夜夜摸| 成人av在线播放网站| 精品人妻1区二区| 欧美一级毛片孕妇| 国产免费av片在线观看野外av| 两个人看的免费小视频| 亚洲精品国产精品久久久不卡| 亚洲av成人一区二区三| 亚洲精品中文字幕一二三四区| 婷婷精品国产亚洲av| 毛片女人毛片| 久久欧美精品欧美久久欧美| 亚洲精品456在线播放app | 俺也久久电影网| 午夜视频精品福利| 国内精品久久久久久久电影| 国产精品久久久av美女十八| svipshipincom国产片| 中文字幕人妻丝袜一区二区| 亚洲中文日韩欧美视频| 搡老妇女老女人老熟妇| 精品福利观看| 精品欧美国产一区二区三| 神马国产精品三级电影在线观看| 免费看日本二区| 日韩欧美一区二区三区在线观看| 国产在线精品亚洲第一网站| a级毛片a级免费在线| 757午夜福利合集在线观看| 午夜免费激情av| 亚洲 国产 在线| 欧美一级毛片孕妇| 亚洲,欧美精品.| 亚洲国产欧美网| 母亲3免费完整高清在线观看| 精品不卡国产一区二区三区| 国产真实乱freesex| 亚洲美女黄片视频| 高清在线国产一区| 国产三级黄色录像| 99国产精品一区二区三区| 久久久久九九精品影院| 国产乱人视频| 精品午夜福利视频在线观看一区| 国产91精品成人一区二区三区| 99精品在免费线老司机午夜|