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

    基于間斷有限元模型的給排水干濕分離研究?

    2017-10-16 09:07:40楊金英
    計算機(jī)與數(shù)字工程 2017年9期
    關(guān)鍵詞:淺水水深水體

    楊金英

    基于間斷有限元模型的給排水干濕分離研究?

    楊金英

    (秦皇島市山海關(guān)排水管理處 秦皇島 066200)

    針對給排水干濕分離難以確定運(yùn)動交界的問題,該研究將利用離散格式有限元建立淺水方程模式,結(jié)合細(xì)分化單元內(nèi)物理量進(jìn)行線性重構(gòu),提出了一種新的給排水干濕分離方法。該方法保證了干濕交界處水體靜止?fàn)顟B(tài)不被破壞,并且在水體非物理流動條件下的干濕處的過渡區(qū)域之間可連續(xù)。經(jīng)過數(shù)值實(shí)驗(yàn)結(jié)果表明,采用間斷式的有限元模式能夠準(zhǔn)確地模擬給排水的水體情形,計算過程中可精確計算間斷間運(yùn)動過程。

    有限元模型;給排水;干濕分離;淺水方程

    AbstractIn order to solve the problem that the movement boundary is difficult to be determ ined by the dry and wet separation of water supply and drainage,a new method of water and wet separation is proposed by using discrete finite element method to estab?lish the shallow water equation model and linear reconstruction with the physical quantities in the subdividing unit.The method can ensure that the static state of the water and wet junction is not destroyed,and can be continuous between the transitional areas of the wet and dry place under the non-physical flow condition of the water body.The results of numerical experiments show that the inter?mittent finite element model can accurately simulate the situation of water supply and drainage,and the process of intermittent mo?tion can be accurately calculated.

    Key Wordsfnite element model,water supply and drainage,wet and dry separation,shallow water equation

    Class NumberTN915.6

    1 引言

    有限元方法被廣泛地應(yīng)用于求解水體流動模型中,利用改進(jìn)數(shù)值算法的形式將有限元劃分與空間數(shù)階相結(jié)合可以得到高精度的數(shù)值解[1~3]。在面對給排水工程項(xiàng)目時,間斷有限元模型可通過淺水域模擬得到非結(jié)構(gòu)化網(wǎng)格單元,確保了流體運(yùn)動中的初始條件與流動邊界的有機(jī)結(jié)合[4]。在模擬相流問題和激流問題過程中,傳統(tǒng)的數(shù)值求解在處理強(qiáng)對流問題時易出現(xiàn)數(shù)值求解的耗散并導(dǎo)致模擬失真的結(jié)果[5]。間斷有限元模型在模擬給排水這類間斷性流體可避免連續(xù)性空間單元數(shù)值求解,因此本文為了研究給排水過程中的干濕分離技術(shù),利用間斷有限元模型分析水體淺水方程組,采用HLL函數(shù)的形式來對邊界處數(shù)值的通量加以計算,結(jié)合水深非負(fù)算子對間斷有限元模型進(jìn)行Runge-Kutta格式化以提高算法精度,通過水體體積守恒定律對干濕處過渡進(jìn)行數(shù)值分析。

    2 間斷有限元模型

    2.1 淺水方程

    在計算水體質(zhì)量時采用固定網(wǎng)格的方法并且需要保證單元的拓?fù)潢P(guān)系穩(wěn)定,由于基于水體質(zhì)量守恒定律[6],因此整個過程牽涉到的計算量較小。通常會采用一種線性近似方法來對單元內(nèi)的水深進(jìn)行重新構(gòu)建[7],以此達(dá)到干濕交接單元內(nèi)水深為負(fù)的目的。但是采用上述的方法可能存在弊端,可能會造成水體自由表面計算錯誤,這樣的情況就使得單元內(nèi)的水體向周圍鄰近的單元內(nèi)形成非物理性的流動,同時會導(dǎo)致單元內(nèi)干濕交界面位置的計算錯誤[8]。采用上述的方法還會導(dǎo)致水體干濕界面捕捉困難的問題。存在這樣的原因主要在于真實(shí)的界面有可能會存在于水體單元的內(nèi)部而并非水體單元的表面或者交界處[9]。除上述情況之外,采用近似線性的薄層方法來處理,可能會使得淺水方程和諧性質(zhì)遭到破壞,這樣就無法對水體靜止的情況進(jìn)行模擬仿真,加大了研究的難度。即水體靜止形式如下[10]

    盡管列出了上述的一些局限性,大多數(shù)的研究仍然利用薄層方法對水體加以計算,因?yàn)槠溆嬎愕牧枯^小,且使用的條件和范圍比較廣泛,因此仍然能夠廣泛的應(yīng)用于各種水體的淺水方程模型當(dāng)中。

    關(guān)于一維的水體淺水方程組,一般由兩個方程所組成,一個方程用于描述水體的質(zhì)量守恒,另一個方程用于描述物理動量恒等。用于描述能量守恒的方程形式如式(2),其中U代表能量守恒,F(xiàn)(U)代表流量項(xiàng),S(U)代表源項(xiàng)等,對于方程中每一項(xiàng)的定義如式(3)。

    在上式的描述中,關(guān)于水體的總水深用h表示,水體的水流速度用u表示,q=hu則表示能量守恒,動量用x表示,關(guān)于水體中的底坡的高程則利用b來表示。在水中的重力加速度用g來表示。為了方便研究的進(jìn)行,在研究中我們將只考慮一種情況,即無粘情形,這樣的情況適用于更加復(fù)雜的附加源項(xiàng),類似于科氏力、底部摩擦項(xiàng)以及風(fēng)應(yīng)力等因素都將被考察在內(nèi)。

    2.2 間斷有限元格式

    為了對間斷有限元格式加以應(yīng)用,我們首先需要對域離散互不重合的區(qū)間Ωi=[xi-1/2,xi+1/2]加以計算。利用上述的域離散互不重合控制方程(4)乘以實(shí)驗(yàn)函數(shù)φ∈Vp(Ωi),同時在不重合區(qū)間Ωi內(nèi)進(jìn)行積分的運(yùn)算,則可以得到方程的弱解:

    在上述的弱解形式中,Uh表示物理能量守恒變量的U的近似解。與實(shí)驗(yàn)函數(shù)φ類似,它們同屬于解空間Vp(Ωi)內(nèi)。在Ωi內(nèi),存在不超過 p的多項(xiàng)式共同組成了解空間Vp(Ωi),這樣的解空間使得多項(xiàng)式在單元的邊界不連續(xù)成為可能。同時為了使得解空間內(nèi)的單元間能夠?qū)崿F(xiàn)信息的交換和處理[11],將利用數(shù)值通量 F*(Uh)來代替通量項(xiàng)F(Uh)。在本文的研究中,將采用HLL函數(shù)的形式來對邊界處數(shù)值的通量加以計算[12],同時為了獲取到半離散的格式,本文中選用的基函數(shù)為Lagrang?ian函數(shù)[13]。利用上述提出的方法,我們將得到水體一維淺水方程的半離散化的形式,同時將其表示為常微分的形式,如下

    一般情況下關(guān)于半離散形式的常微分方程,則能夠?qū)崿F(xiàn)二階顯式TVD Runge-Kutta格式化的計算,詳細(xì)的推理過程如下

    在推理的式中,Runge-Kutta的各個分步的系數(shù)一般用 βil以及αil來表示。Un以及Un+1則用于表示n以及n+1時刻的數(shù)值的解。通過對上述的模型采用TVB限制器,實(shí)現(xiàn)對于每一步數(shù)值的解的計算,同時還需要對數(shù)值振蕩進(jìn)行控制以保證數(shù)據(jù)的格式穩(wěn)定。通過上述的公式與模型,我們最終能夠得到空間精度的Runge-Kutta以及具有二階時間的有限元間斷淺水方程模型。

    3 給排水干濕分離

    3.1 水深非負(fù)算子

    對于水深我們需要設(shè)定其為正值,因此本文將采用水深非負(fù)算子,并將其作用于每個Runge-Kut?ta的分步過程中[14]。為了方便研究的進(jìn)行,我們需要對一些變量進(jìn)行標(biāo)識。Ui+=(hi+,q+i)用于表示非重復(fù)區(qū)間Ωi內(nèi)的變量均值,Ui=(hi,qi)用于表示節(jié)點(diǎn)i∈{1,2 }處的解。其中非負(fù)算子的具體表達(dá)式如下

    上式中,hmini表示非重復(fù)區(qū)間Ωi內(nèi)的最小水深。根據(jù)上式的內(nèi)容不難發(fā)現(xiàn),對于所有的水深進(jìn)行修正后得到的水深值都為非負(fù)數(shù)值。同時還能夠保持水體單元內(nèi)的總體積以及總動量不發(fā)生變化。關(guān)于水深的平均值為負(fù)的干單元,則需要將其平均水深重新修正為h+i=0。

    3.2 干單元處理方法

    3.3 干濕處過渡

    為了避免Gibbs現(xiàn)象的產(chǎn)生,本部分將采用對水體單元進(jìn)行細(xì)分的方法實(shí)現(xiàn)對過渡單元的重新構(gòu)建。完成干濕邊界分區(qū)之后,將過渡的水體單元以及為節(jié)點(diǎn)分為兩個子單元。同時為了使得干濕交界點(diǎn)的流速較小,也將以動量q為對象使用線性重構(gòu)方法在子單元實(shí)現(xiàn)對動量的重新分布。細(xì)分后h與q在子單元的分布將如下所示。

    對于干濕過渡區(qū)域處于全濕狀態(tài)的單元,這樣的重新構(gòu)建的過程要較為簡單。全濕狀態(tài)的單元底坡以及水體自由面所組成的形狀為梯形形狀。水深連續(xù)以及參照物理水體體積守恒定律,單元內(nèi)水體深度的分布由以下公式確定:

    4 實(shí)驗(yàn)分析

    在本部分的試驗(yàn)中,將不均勻底坡上水面靜止的情況考慮在內(nèi),同時需要驗(yàn)證在干濕交界時模型對于水體數(shù)值的處理能夠保持和諧的性質(zhì)。設(shè)定[0,1]為計算域,同時需要將區(qū)間劃分為100個均勻的小單元。底坡的高程如下式計算:

    其中r=|x-0.5|,a的系數(shù)設(shè)定為1.3,rm的系數(shù)設(shè)為0.5。關(guān)于水體的初始水位高程則設(shè)定為η(x)=max(b(x),1),水體靜止即速度為0。根據(jù)CFL的條件,時間步的長度取為0.0001s,而模擬2s后則計算過程結(jié)束,結(jié)果如圖1所示。

    圖1 給排水干濕分離計算結(jié)果

    由圖1可見,對于干濕的處理方式有不同的方法,(a)與(b)為進(jìn)行比較的結(jié)果。CONV為針對傳統(tǒng)的線性重新構(gòu)建方法得出的計算結(jié)果,HREF為針對滿足和諧性質(zhì)的干濕處理線性重新構(gòu)建方法得出的計算結(jié)果。通過結(jié)果的對比顯示,傳統(tǒng)干濕處理方法無法對于水體的狀態(tài)穩(wěn)定下來,即無法保證水體的靜止。如圖1(b)所示,模擬2s后,在濕的單元內(nèi)水體流入能夠過渡到過渡單元內(nèi),從而使得動量q由干濕交界的單元向水深值較大的區(qū)域流動。同時對于新的干濕區(qū)域處理方法,通過將過渡的單元劃分為干和濕兩個部分,并使其保持連續(xù),這樣在動量q的計算過程中就能避免大量的舍入誤差的產(chǎn)生。

    5 結(jié)語

    本研究主要圍繞給排水過程中的干濕分離技術(shù),以水體靜止情形為出發(fā)點(diǎn),利用間斷有限元模型對離散互不重合的區(qū)間進(jìn)行計算。采用HLL函數(shù)的形式對邊界處數(shù)值的通量加以計算,同時,選用的基函數(shù)為Lagrangian函數(shù)得到水體一維淺水方程的半離散化的形式和二階顯式TVD Run?ge-Kutta格式化。運(yùn)用水深非負(fù)算子和Lagrangian基函數(shù)作用于每個Runge-Kutta的分步。參照Lewy的限制條件將過渡的水體單元以及為節(jié)點(diǎn)分為兩個子單元,使用線性重構(gòu)方法在子單元實(shí)現(xiàn)對動量的重新分布。

    [1]唐亮,凌賢長,艾哈邁德·艾格瑪.液化側(cè)向流動場地樁基動力反應(yīng)振動臺試驗(yàn)三維有限元數(shù)值模擬方法[J].土木工程學(xué)報,2013(1):180-184.

    TANG Liang,LIN Xianchang,AIGema.Three Dimension?al Finite Element Numerical Simulation of Shaking Table Test for Dynamic Response of Pile Foundation in Lique?fied Lateral Flowing[J].China Civil,2013(1):180-184.

    [2]魯陽,鄒建鋒,鄭耀.基于非結(jié)構(gòu)網(wǎng)格的TTGC有限元格式的實(shí)現(xiàn)及在超聲速流動中的應(yīng)用[J].計算力學(xué)學(xué)報,2013(5):712-716.

    LU Yang,ZOU Jianfeng,ZHANG Yao.Realization of TT?GC Finite Element Scheme Based on Unstructured Grid and Its Application in Supersonic Flow[J].Chinese Jour?nal of Computational Mechanics,2013(5):712-716.

    [3]魏匡民,蘇佩珍.有限元法中非關(guān)聯(lián)流動條件下剛度矩陣的存儲方式及平衡方程求解[J].三峽大學(xué)學(xué)報(自然科學(xué)版),2011,33(5):10-14.

    WEI Kuangmin,SU Peizhen.Storing Method of Stiffness Matrix and Solving Equilibrium Equation in Non-associat?ed Flow Condition[J].Journal of China Three Gorges Uni?versity Natural Sciences,2011,33(5):10-14.

    [4]張忠宇,姚熊亮,張阿漫.基于間斷有限元方法的并列圓柱層流流動特性[J].物理學(xué)報,2016,65(8):1-11.

    ZHANG Zhongyu,YAO Xiongliang,ZHANG Aman.Flow characteristics of parallel laminar flow based on discontin?uous finite element method[J].Acta Physica Sinica,2016,65(8):1-11.

    [5]楊青遠(yuǎn),周蘇芬,盧偉真,等.明渠交匯水流阻力對平面二維數(shù)值模擬的影響[J].水科學(xué)進(jìn)展,2012,23(2):236-242.

    YANG Qingyuan,ZHOU Sufen,LU Weizhen,et al.Influ?ence of Flow Resistance in Open Channel Crossing on Two-dimensional Numerical Simulation[J].Advances in Water Science.2012,23(2):236-242.

    [6]馬鑫磊,任光明,夏敏.滑坡涌浪預(yù)測評價方法綜述[J].水利水電科技進(jìn)展,2015(3):89-98.

    MA Xinlei,REN Guangming,XIA Min.Summarization of forecasting and evaluating methods of landslide surge[J].Advances in Science and Technology of Water Resources,2015(3):89-98.

    [7]李艷華,丁建麗,閆人華.基于國產(chǎn)GF-1遙感影像的山區(qū)細(xì)小水體提取方法研究[J].資源科學(xué),2015,37(2):408-416.

    LI Yanhua,DING Jianli,YAN Renhua.Study on Extrac?tion of Small Waters in Mountainous Area Based on Do?mestic GF-1 Remote Sensing Image[J].Resources Sci?ence,2015,37(2):408-416.

    [8]吳鋒,鐘萬勰.淺水動邊界問題的位移法模擬[J].計算機(jī)輔助工程,2016,25(2):5-13.

    WU Feng,ZHONGWanxie.Displacement method for shal?low water dynamic boundary[J].Computer Aided Engi?neering,2016,25(2):5-13.

    [9]王增輝,夏軍強(qiáng),李濤.兩步計算模式下的三門峽水庫異重流模擬[J].四川大學(xué)學(xué)報(工程科學(xué)版),2016,25(2):5-13.

    WANG Zenghui,XIA Junqiang,LI Tao.Simulations of three-grade calculation of three-gorges reservoirs with density flow[J].Journal of Sichuan University Engineer?ing Science Edition,2016,25(2):5-13.

    [10]程文龍,史英標(biāo),吳修廣.基于AUSM+-UP格式的有限體積模型及其在錢塘江河口的應(yīng)用[J].水力發(fā)電學(xué)報,2013,32(3):82-86.

    CHENG Wenlong,SHI Yingbiao,WU Xiuguang.Finite Volume Model Based on AUSM-UPScheme and Its Ap?plication in Qiantang Estuary[J].Journal of Hydroelec?tric Engineering,2013,32(3):82-86.

    [11]徐立云,榮巨,郭昆吾.基于改進(jìn)粒子群算法的制造單元構(gòu)建方法研究[J].制造技術(shù)與機(jī)床,2014(3):47-53.

    XU Liyun,RONG Ju,GUO Kunwu.Research on Manu?facturing Cell Construction Based on Improved Particle Swarm Optimization[J].Manufacturing Technology&Machine Tool,2014(3):47-53.

    [12]楊金波,李訂芳,陳華.和諧WAF格式在帶干河床淺水波方程中的應(yīng)用[J].華中科技大學(xué)學(xué)報(自然科學(xué)版),2012,40(2):54-57.

    YANG Jinbo,LI Dingfang,CHEN Hua.Application of Harmonic WAF Scheme in Shallow Water Wave Equa?tions with Dry Riverbed[J].Journal of Huazhong Univer?sity of Science and Technology Natural Science Edition,2012,40(2):54-57.

    [13]李俊杰,張穎.高分子稀溶液的多尺度模型研究進(jìn)展[J].廈門大學(xué)學(xué)報(自然科學(xué)版),2013,52(2):289-296.

    LI Junjie,ZHANG Ying.Progress in Research on Multi?scale Models of Polymer Dilute Solution[J].Journal of Xiamen University Natural Science,2013,52(2):289-296.

    [14]趙張益,張慶河,李世森.Runge-Kutta間斷有限元格式在一維淺水方程中的應(yīng)用[J].水科學(xué)進(jìn)展,2012,23(5):695-701.

    ZHAO Zhangyi,ZHANG Qinghe,LI Shisen.Application of Runge-Kutta Discontinuous Finite Element Scheme in One Dimensional Shallow Water Equation[J].Advances in Water Science,2012,23(5):695-701.

    [15]劉存良,潘振寬,鄭永果.兩種保持符號距離函數(shù)的水平集分割方法[J].吉林大學(xué)學(xué)報(工學(xué)版),2013,43(3):115-119.

    LIU Cunliang,PAN Zhenkuan,ZHENG Yongguo.Two Level Set Segmenting Methods for Preserving Symbol Dis?tance Function[J].Journal of Jilin University Engineer?ing and Technology Edition,2013,43(3):115-119.

    Study on Dry and W et Separation of W ater Supp ly and Drainage Based on Discontinuous Finite Elem ent

    YANG Jinying
    (Drainage Administrative Department of Shanhaiguan District of Qinhuangdao City,Qinhuangdao 066200)

    10.3969/j.issn.1672-9722.2017.09.043

    2017年3月8日,

    2017年4月23日

    楊金英,女,碩士,高級工程師,研究方向:給排水方面。

    猜你喜歡
    淺水水深水體
    書法靜水深流
    河北水利(2022年10期)2022-12-29 11:48:12
    基于水深分段選擇因子的多光譜影像反演水深
    海洋通報(2022年4期)2022-10-10 07:40:32
    農(nóng)村黑臭水體治理和污水處理淺探
    新型淺水浮托導(dǎo)管架的應(yīng)用介紹
    云南化工(2021年10期)2021-12-21 07:33:40
    生態(tài)修復(fù)理念在河道水體治理中的應(yīng)用
    廣元:治理黑臭水體 再還水清岸美
    帶阻尼的隨機(jī)淺水波方程的隨機(jī)吸引子
    GPS RTK技術(shù)在水深測量中的應(yīng)用
    (2+1)維廣義淺水波方程的Backlund變換和新精確解的構(gòu)建
    找不同
    免费在线观看影片大全网站| 男人舔奶头视频| 最近最新中文字幕大全电影3| 国产成人精品无人区| 老熟妇仑乱视频hdxx| 俺也久久电影网| 国产亚洲精品久久久久久毛片| 欧美3d第一页| 好男人电影高清在线观看| 精品一区二区三区av网在线观看| 国产欧美日韩精品亚洲av| 久久久久久九九精品二区国产 | 中文字幕最新亚洲高清| 亚洲av电影在线进入| 大型黄色视频在线免费观看| 淫妇啪啪啪对白视频| 久久午夜亚洲精品久久| 少妇被粗大的猛进出69影院| 成年版毛片免费区| 91麻豆av在线| 欧美日韩乱码在线| 国产精品98久久久久久宅男小说| 美女扒开内裤让男人捅视频| 俄罗斯特黄特色一大片| 国产精品久久久久久亚洲av鲁大| 亚洲av成人精品一区久久| 亚洲无线在线观看| 亚洲中文日韩欧美视频| 1024手机看黄色片| 免费av毛片视频| 麻豆av在线久日| 国产又色又爽无遮挡免费看| 一区二区三区国产精品乱码| 亚洲色图 男人天堂 中文字幕| 精品乱码久久久久久99久播| 熟女电影av网| 热99re8久久精品国产| 老汉色∧v一级毛片| 两个人视频免费观看高清| 午夜免费成人在线视频| 黄色毛片三级朝国网站| 亚洲av第一区精品v没综合| 嫩草影视91久久| 日韩欧美在线二视频| 狠狠狠狠99中文字幕| 国内揄拍国产精品人妻在线| 日本一本二区三区精品| 国产麻豆成人av免费视频| 成在线人永久免费视频| 午夜福利在线在线| √禁漫天堂资源中文www| АⅤ资源中文在线天堂| 999久久久国产精品视频| 国产99久久九九免费精品| 老鸭窝网址在线观看| 此物有八面人人有两片| 伦理电影免费视频| 久9热在线精品视频| 女生性感内裤真人,穿戴方法视频| www.www免费av| 少妇人妻一区二区三区视频| 激情在线观看视频在线高清| 婷婷丁香在线五月| 国产精品久久久久久人妻精品电影| 欧美在线黄色| 免费在线观看视频国产中文字幕亚洲| 又爽又黄无遮挡网站| 亚洲狠狠婷婷综合久久图片| 欧美精品啪啪一区二区三区| 精品一区二区三区四区五区乱码| 露出奶头的视频| 亚洲真实伦在线观看| 国产人伦9x9x在线观看| 欧美成人一区二区免费高清观看 | 亚洲精品色激情综合| 亚洲中文字幕日韩| 日韩欧美国产在线观看| 黄片大片在线免费观看| 久久久水蜜桃国产精品网| 啦啦啦免费观看视频1| 精品国产亚洲在线| 黑人操中国人逼视频| 精品无人区乱码1区二区| 午夜精品久久久久久毛片777| 天堂√8在线中文| 婷婷精品国产亚洲av在线| 国语自产精品视频在线第100页| 国产精品一区二区精品视频观看| 午夜久久久久精精品| 国模一区二区三区四区视频 | 中亚洲国语对白在线视频| 天堂影院成人在线观看| 亚洲国产精品999在线| 看黄色毛片网站| 欧美人与性动交α欧美精品济南到| 麻豆国产97在线/欧美 | 最近在线观看免费完整版| 国产一区二区三区在线臀色熟女| 国产在线观看jvid| 日日摸夜夜添夜夜添小说| 国产区一区二久久| 国产精品 国内视频| 99精品欧美一区二区三区四区| 成人av一区二区三区在线看| 久久精品亚洲精品国产色婷小说| 一本精品99久久精品77| 热99re8久久精品国产| av在线天堂中文字幕| 亚洲自拍偷在线| 香蕉久久夜色| 蜜桃久久精品国产亚洲av| 两性夫妻黄色片| 男女之事视频高清在线观看| 99久久精品热视频| 三级国产精品欧美在线观看 | av欧美777| 妹子高潮喷水视频| 国产精品国产高清国产av| 国产成人啪精品午夜网站| 露出奶头的视频| 此物有八面人人有两片| 欧美激情久久久久久爽电影| 黄色a级毛片大全视频| 麻豆国产av国片精品| 99国产精品一区二区蜜桃av| 亚洲真实伦在线观看| 欧洲精品卡2卡3卡4卡5卡区| 丝袜美腿诱惑在线| 午夜a级毛片| 久久久久免费精品人妻一区二区| 欧美在线一区亚洲| 久久久久久国产a免费观看| 色噜噜av男人的天堂激情| 国产成人一区二区三区免费视频网站| 一夜夜www| 深夜精品福利| 搡老熟女国产l中国老女人| 黄色a级毛片大全视频| 真人一进一出gif抽搐免费| 美女免费视频网站| 老司机靠b影院| 欧美激情久久久久久爽电影| 窝窝影院91人妻| 欧美乱妇无乱码| 欧美性猛交黑人性爽| 中文字幕人成人乱码亚洲影| 精品一区二区三区视频在线观看免费| 国产1区2区3区精品| 亚洲天堂国产精品一区在线| 两人在一起打扑克的视频| 国产又色又爽无遮挡免费看| 亚洲精品中文字幕在线视频| 丁香六月欧美| 色尼玛亚洲综合影院| 一进一出抽搐动态| 欧美成狂野欧美在线观看| 国产精品久久久久久精品电影| 高清在线国产一区| av超薄肉色丝袜交足视频| 最近视频中文字幕2019在线8| www日本黄色视频网| 日本熟妇午夜| 中文字幕久久专区| 亚洲aⅴ乱码一区二区在线播放 | 麻豆成人av在线观看| 国产伦在线观看视频一区| 国产成人av激情在线播放| 一卡2卡三卡四卡精品乱码亚洲| 在线十欧美十亚洲十日本专区| 久久九九热精品免费| 91国产中文字幕| 欧美人与性动交α欧美精品济南到| 亚洲成人久久爱视频| 国产区一区二久久| 成人三级做爰电影| 日韩欧美在线二视频| 久久热在线av| 男女那种视频在线观看| 国产精品野战在线观看| 好看av亚洲va欧美ⅴa在| 国产亚洲精品av在线| 久久久久久久久中文| 久久久久久久久中文| 色综合婷婷激情| 久久久精品大字幕| 国内精品久久久久精免费| 国产熟女xx| netflix在线观看网站| 久久香蕉精品热| 香蕉久久夜色| 99久久99久久久精品蜜桃| 观看免费一级毛片| 亚洲精品国产一区二区精华液| 午夜福利视频1000在线观看| 久久九九热精品免费| 香蕉丝袜av| 亚洲精品美女久久av网站| 一级a爱片免费观看的视频| 国产99白浆流出| 欧美乱码精品一区二区三区| 亚洲av成人精品一区久久| 性欧美人与动物交配| 免费在线观看亚洲国产| 亚洲熟妇熟女久久| 巨乳人妻的诱惑在线观看| 日本黄色视频三级网站网址| 欧美日韩精品网址| 在线观看66精品国产| 成人高潮视频无遮挡免费网站| 一级毛片高清免费大全| 99精品在免费线老司机午夜| 淫妇啪啪啪对白视频| www.熟女人妻精品国产| 久久久久国产精品人妻aⅴ院| 精品国产超薄肉色丝袜足j| 老司机靠b影院| 亚洲自拍偷在线| 国产av一区二区精品久久| 国产精品久久久人人做人人爽| 国产蜜桃级精品一区二区三区| 国产精品久久久久久亚洲av鲁大| 国产三级中文精品| 成人国产一区最新在线观看| 嫩草影视91久久| 精品人妻1区二区| 国产精品亚洲一级av第二区| 精品国产超薄肉色丝袜足j| 精品国产亚洲在线| 国产一区二区三区视频了| 又粗又爽又猛毛片免费看| 男人舔奶头视频| 久久人人精品亚洲av| 国产三级在线视频| 午夜福利欧美成人| 久久中文看片网| 亚洲成人中文字幕在线播放| 巨乳人妻的诱惑在线观看| av在线播放免费不卡| 99re在线观看精品视频| 久久精品国产清高在天天线| 欧美黑人欧美精品刺激| 国产欧美日韩精品亚洲av| 久久久久久大精品| 国产一区二区激情短视频| 午夜福利免费观看在线| 午夜亚洲福利在线播放| 久久天躁狠狠躁夜夜2o2o| 国产伦人伦偷精品视频| 一个人观看的视频www高清免费观看 | 国产伦在线观看视频一区| 亚洲av五月六月丁香网| 在线a可以看的网站| 亚洲 欧美 日韩 在线 免费| 少妇粗大呻吟视频| 国产三级在线视频| 18禁观看日本| 人人妻,人人澡人人爽秒播| 国产精品久久久久久亚洲av鲁大| 国产野战对白在线观看| 国产精品久久久久久精品电影| 久久精品影院6| 国产亚洲精品久久久久5区| 麻豆av在线久日| 亚洲成人久久爱视频| 欧美成人免费av一区二区三区| 国产aⅴ精品一区二区三区波| 99国产精品99久久久久| 欧美日韩国产亚洲二区| 国产黄片美女视频| xxxwww97欧美| 日本免费a在线| 精品欧美一区二区三区在线| 日本三级黄在线观看| 成人欧美大片| 免费在线观看视频国产中文字幕亚洲| 久热爱精品视频在线9| 亚洲 国产 在线| 午夜福利欧美成人| 精华霜和精华液先用哪个| 免费在线观看完整版高清| 亚洲性夜色夜夜综合| 99久久无色码亚洲精品果冻| 香蕉国产在线看| 两人在一起打扑克的视频| 国产探花在线观看一区二区| 老司机午夜福利在线观看视频| 免费看日本二区| 50天的宝宝边吃奶边哭怎么回事| 亚洲男人的天堂狠狠| 色综合欧美亚洲国产小说| 午夜精品一区二区三区免费看| 一区二区三区国产精品乱码| 手机成人av网站| 亚洲国产欧洲综合997久久,| 欧美+亚洲+日韩+国产| 国产精品av久久久久免费| 男女床上黄色一级片免费看| 国产精品免费视频内射| 美女高潮喷水抽搐中文字幕| 午夜老司机福利片| 性欧美人与动物交配| 免费在线观看成人毛片| 十八禁网站免费在线| 欧美国产日韩亚洲一区| 最新美女视频免费是黄的| 亚洲无线在线观看| 男女做爰动态图高潮gif福利片| 国产一级毛片七仙女欲春2| 麻豆久久精品国产亚洲av| 免费看a级黄色片| 久久中文看片网| 国产精品免费一区二区三区在线| 嫩草影视91久久| 国产亚洲精品第一综合不卡| 一进一出好大好爽视频| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲av片天天在线观看| 欧美日韩亚洲国产一区二区在线观看| 久久国产精品人妻蜜桃| 成人国产一区最新在线观看| 99精品欧美一区二区三区四区| 在线永久观看黄色视频| 19禁男女啪啪无遮挡网站| 日本成人三级电影网站| 国产黄片美女视频| 动漫黄色视频在线观看| 精品国产亚洲在线| 久热爱精品视频在线9| 亚洲精品在线观看二区| 日本精品一区二区三区蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 哪里可以看免费的av片| 欧美成人一区二区免费高清观看 | 一级毛片高清免费大全| 91九色精品人成在线观看| 国产午夜精品论理片| 精品午夜福利视频在线观看一区| 亚洲av电影不卡..在线观看| 国产人伦9x9x在线观看| 特级一级黄色大片| 特大巨黑吊av在线直播| 亚洲国产中文字幕在线视频| 可以在线观看的亚洲视频| av国产免费在线观看| 久久这里只有精品中国| 99久久精品热视频| 免费看a级黄色片| 亚洲男人天堂网一区| АⅤ资源中文在线天堂| 国产精品av久久久久免费| 黄片大片在线免费观看| 午夜激情福利司机影院| 欧美中文综合在线视频| 亚洲av五月六月丁香网| 久久欧美精品欧美久久欧美| 丰满的人妻完整版| 极品教师在线免费播放| 午夜福利视频1000在线观看| 女人被狂操c到高潮| 91老司机精品| 激情在线观看视频在线高清| 在线十欧美十亚洲十日本专区| 亚洲中文字幕日韩| 757午夜福利合集在线观看| 在线观看www视频免费| 国产精品影院久久| 黄片小视频在线播放| а√天堂www在线а√下载| tocl精华| 法律面前人人平等表现在哪些方面| 法律面前人人平等表现在哪些方面| 色尼玛亚洲综合影院| 亚洲熟女毛片儿| 天天躁狠狠躁夜夜躁狠狠躁| 天堂动漫精品| 成人av在线播放网站| 国产1区2区3区精品| 国产精品99久久99久久久不卡| 亚洲国产精品合色在线| 国产精品免费视频内射| 亚洲欧美日韩无卡精品| 日韩欧美精品v在线| 久久久水蜜桃国产精品网| 51午夜福利影视在线观看| 欧美丝袜亚洲另类 | bbb黄色大片| 9191精品国产免费久久| 国产视频一区二区在线看| 村上凉子中文字幕在线| 国产精品av久久久久免费| 亚洲欧美日韩高清专用| 午夜福利成人在线免费观看| 国内精品久久久久精免费| 999久久久国产精品视频| 国产免费av片在线观看野外av| 在线观看美女被高潮喷水网站 | 最近最新中文字幕大全电影3| 中文亚洲av片在线观看爽| 美女午夜性视频免费| 国产精品日韩av在线免费观看| www.www免费av| 亚洲av成人精品一区久久| 国产精品98久久久久久宅男小说| 男人舔奶头视频| 黄色视频,在线免费观看| 久久久久九九精品影院| 午夜福利成人在线免费观看| 在线a可以看的网站| 国产精品免费一区二区三区在线| 欧美日韩中文字幕国产精品一区二区三区| 麻豆av在线久日| 欧美另类亚洲清纯唯美| 在线a可以看的网站| 搡老岳熟女国产| 在线观看一区二区三区| 母亲3免费完整高清在线观看| 很黄的视频免费| 久久久精品大字幕| 嫩草影视91久久| 黄色视频不卡| 夜夜夜夜夜久久久久| 一区二区三区国产精品乱码| 日韩成人在线观看一区二区三区| 欧美乱妇无乱码| 成熟少妇高潮喷水视频| 国产伦在线观看视频一区| 18禁黄网站禁片免费观看直播| 亚洲性夜色夜夜综合| 黄片大片在线免费观看| 欧美最黄视频在线播放免费| 久久国产精品影院| 日本精品一区二区三区蜜桃| 免费看十八禁软件| 制服诱惑二区| 国内久久婷婷六月综合欲色啪| 又黄又爽又免费观看的视频| 丰满的人妻完整版| cao死你这个sao货| 亚洲欧美激情综合另类| 亚洲中文av在线| 亚洲国产精品成人综合色| 又爽又黄无遮挡网站| 国产精品免费视频内射| 欧美在线黄色| 国产精品爽爽va在线观看网站| 欧美另类亚洲清纯唯美| 欧美性长视频在线观看| 精品人妻1区二区| 麻豆av在线久日| 国产人伦9x9x在线观看| 国语自产精品视频在线第100页| 亚洲熟妇熟女久久| 别揉我奶头~嗯~啊~动态视频| 高潮久久久久久久久久久不卡| 亚洲成人精品中文字幕电影| 欧美性猛交╳xxx乱大交人| 国产视频内射| 国产黄片美女视频| 国产精品久久久久久精品电影| 最近最新免费中文字幕在线| 日本三级黄在线观看| 久久久久久九九精品二区国产 | 少妇裸体淫交视频免费看高清 | 大型av网站在线播放| 99热这里只有是精品50| 一区二区三区激情视频| 欧美黄色片欧美黄色片| 这个男人来自地球电影免费观看| 午夜日韩欧美国产| 国产精品影院久久| 18禁裸乳无遮挡免费网站照片| 中文字幕最新亚洲高清| 十八禁人妻一区二区| 国产精品一区二区免费欧美| 久久婷婷成人综合色麻豆| 51午夜福利影视在线观看| 特大巨黑吊av在线直播| 成人国产综合亚洲| 久久精品国产亚洲av香蕉五月| 久久这里只有精品19| 搡老熟女国产l中国老女人| 久久国产乱子伦精品免费另类| 九九热线精品视视频播放| 女警被强在线播放| 全区人妻精品视频| 午夜福利免费观看在线| 亚洲国产精品久久男人天堂| 欧美成人午夜精品| 男女床上黄色一级片免费看| 又爽又黄无遮挡网站| 日本熟妇午夜| 国产又色又爽无遮挡免费看| 大型av网站在线播放| 国产av一区在线观看免费| 又黄又爽又免费观看的视频| 午夜福利18| 国产亚洲精品久久久久5区| 国产成人啪精品午夜网站| 国产黄片美女视频| 国产精品九九99| 成人亚洲精品av一区二区| 老熟妇仑乱视频hdxx| 亚洲黑人精品在线| 成人18禁在线播放| 国产成人欧美在线观看| 国产午夜福利久久久久久| 51午夜福利影视在线观看| 久久热在线av| 色在线成人网| 可以在线观看的亚洲视频| 夜夜看夜夜爽夜夜摸| 日本a在线网址| 亚洲国产欧美网| 免费高清视频大片| 午夜成年电影在线免费观看| 日本五十路高清| 欧美日韩福利视频一区二区| 欧美三级亚洲精品| 中国美女看黄片| 欧美不卡视频在线免费观看 | 欧美中文综合在线视频| 亚洲国产精品999在线| 国产成人av教育| 一区二区三区激情视频| 18美女黄网站色大片免费观看| 又黄又爽又免费观看的视频| 午夜亚洲福利在线播放| 在线观看舔阴道视频| 亚洲国产欧美人成| 91大片在线观看| 9191精品国产免费久久| 国产精品日韩av在线免费观看| www国产在线视频色| 亚洲中文字幕日韩| 亚洲av中文字字幕乱码综合| 99国产精品一区二区三区| 亚洲欧美日韩无卡精品| 久久久久国内视频| 在线国产一区二区在线| 搡老岳熟女国产| 久热爱精品视频在线9| 极品教师在线免费播放| 少妇被粗大的猛进出69影院| 亚洲精品中文字幕一二三四区| 一进一出抽搐动态| 精品免费久久久久久久清纯| 国产精品 欧美亚洲| 最好的美女福利视频网| 69av精品久久久久久| 欧美中文综合在线视频| 午夜视频精品福利| a级毛片a级免费在线| 亚洲成av人片免费观看| 国产成人啪精品午夜网站| 国产精品99久久99久久久不卡| 男女下面进入的视频免费午夜| 精品熟女少妇八av免费久了| 国产av麻豆久久久久久久| 亚洲人成网站在线播放欧美日韩| 国产精品亚洲一级av第二区| ponron亚洲| 一个人免费在线观看的高清视频| 老司机午夜福利在线观看视频| 精品国产亚洲在线| 日韩欧美三级三区| 亚洲欧美激情综合另类| 人人妻人人澡欧美一区二区| 亚洲欧美激情综合另类| 中文资源天堂在线| 日韩 欧美 亚洲 中文字幕| 久久九九热精品免费| 亚洲国产欧美网| 国产精品久久久久久亚洲av鲁大| 露出奶头的视频| 中文字幕久久专区| 男人舔女人下体高潮全视频| 久久久国产成人精品二区| 国产亚洲av高清不卡| 黄频高清免费视频| 日本三级黄在线观看| 国产精品99久久99久久久不卡| 搡老妇女老女人老熟妇| 久久久久精品国产欧美久久久| 国产av一区二区精品久久| 久久久久九九精品影院| 国产成人精品久久二区二区免费| 人人妻,人人澡人人爽秒播| 午夜福利免费观看在线| 久久精品国产99精品国产亚洲性色| 日韩欧美 国产精品| 99国产精品一区二区蜜桃av| 在线十欧美十亚洲十日本专区| 黄片小视频在线播放| 久久性视频一级片| 亚洲av片天天在线观看| 日韩成人在线观看一区二区三区| 免费在线观看亚洲国产| 亚洲欧美精品综合一区二区三区| 成人午夜高清在线视频| 国产精品免费视频内射| 99国产精品一区二区蜜桃av| 夜夜躁狠狠躁天天躁| 国产男靠女视频免费网站| 俺也久久电影网| 欧美不卡视频在线免费观看 | 看免费av毛片| 欧美日本亚洲视频在线播放| 99久久国产精品久久久| 亚洲午夜精品一区,二区,三区| 亚洲18禁久久av| 成年女人毛片免费观看观看9| 亚洲 欧美一区二区三区| 欧美+亚洲+日韩+国产| 日本一区二区免费在线视频| 特级一级黄色大片|