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

    數(shù)值水槽模擬波浪運(yùn)動及合理性研究

    2021-12-03 11:34:00楊彤麟甘賢備
    水利規(guī)劃與設(shè)計 2021年12期
    關(guān)鍵詞:造波消波動量

    楊彤麟,羅 浩,晏 亮,甘賢備,郭 輝

    (1.湖南科技大學(xué)土木工程學(xué)院,湖南 湘潭 411201;2.中國鐵道科學(xué)研究院集團(tuán)有限公司鐵道建筑研究所,北京 100081)

    環(huán)境條件一直是橋梁等結(jié)構(gòu)設(shè)計需要考慮的一大關(guān)鍵因素,在某些情況下要面對跨越水域的需求及波浪沖擊橋梁的問題,由于在現(xiàn)場進(jìn)行測量難度大,成本高,又為了能方便準(zhǔn)確地進(jìn)行波浪對結(jié)構(gòu)產(chǎn)生的沖擊作用方面的研究,國內(nèi)外學(xué)者在波浪模擬上做了大量地相關(guān)研究。

    Peter Troch[1]首次實(shí)現(xiàn)了推板造波法及邊界消波,減少了水槽長度,大幅降低了邊緣反射效應(yīng);劉毅[2]指出現(xiàn)有造波機(jī)在一定頻率范圍下振幅保持和精確控制等方面的問題,提出應(yīng)用于造波機(jī)的新型液壓驅(qū)動技術(shù);葉瀟瀟[3]研究了不同造波條件下生成的內(nèi)孤立波運(yùn)動的流場特征;張可心[4]采用二階造波理論,分別進(jìn)行了物理水槽和數(shù)值水槽造波實(shí)驗;虞佳穎[5]利用fluent探究了預(yù)留沖刷坑對挑流水舌在下游水墊塘中消殺能量的效果。

    由于造波實(shí)驗水池的成本較高,試驗結(jié)果受影響因素多,本文總結(jié)了指導(dǎo)造波的波浪理論,并從數(shù)值水槽方面,總結(jié)了其造波方法和消波的原理,并以一個工程實(shí)例驗證了數(shù)值造波的模擬效果。

    1 波浪理論

    當(dāng)前學(xué)者們一般把波浪模擬分為3種:理論推導(dǎo)、物理模擬及數(shù)值模擬,所以波浪研究常采用理論及物理模擬或者理論及數(shù)值模擬。在進(jìn)行數(shù)值模擬波浪運(yùn)動時,假定波浪為一種不可壓縮的無黏性的理想液體,其流體運(yùn)動的速度勢函數(shù)Ф(x,y,z,t)滿足的連續(xù)性拉普拉斯控制方程為[6]:

    ?2φ=0

    (1)

    在描述波浪運(yùn)動的同時,給定其水深為d的初始邊界條件為:

    (2)

    其自由水面z=η(x,t)時的運(yùn)動方程可由伯努利方程積分得到:

    (3)

    式中,g—重力加速度。

    波浪沿x方向傳播則滿足周期性邊界條件[7]:

    φ(x,t)=φ(x+L,t)

    (4)

    φ(x,t)=φ(x,t+T)

    (5)

    式中,L—波長;T—波浪周期。

    1.1 線性波浪理論

    研究人員最初假定波浪運(yùn)動的幅值很小,僅對邊界條件取理想的線性化處理,將重力作為唯一外力,從而水質(zhì)點(diǎn)保持不變的頻率進(jìn)行往復(fù)運(yùn)動,波面呈簡諧波式并以一定波速傳遞,這便是線性波浪理論,又稱為小振幅波浪理論,其波面方程為:

    (6)

    水質(zhì)點(diǎn)的水平與豎直方向速度分別為:

    (7)

    (8)

    波速C為:

    (9)

    式中,k—波數(shù);z—距水底距離;η—瞬時波面高度;H—波高;ω—波的圓頻率。

    1.2 斯托克斯波浪理論

    由于線性波忽略了波高對波浪作用的影響,研究人員進(jìn)而把波陡(波高/波長)作為參數(shù),得到了非線性的stokes波浪理論,其波面方程為:

    (10)

    水質(zhì)點(diǎn)水平方向的運(yùn)動速度為:

    (11)

    水質(zhì)點(diǎn)豎直方向的運(yùn)動速度為:

    (12)

    二階stokes波波速與線性波相同,其波長為:

    (13)

    式中,L—波長;T—波浪周期;H—波高;η—瞬時波面高度;z—距水底距離;k—波數(shù);d—水深;ω—波的圓頻率。

    1.3 橢圓余弦波浪理論

    相對于線性波及斯托克斯波來說,橢圓余弦波是一種適用于近海岸區(qū)的非線性波浪理論,其具有更陡的波峰和更平坦的波谷,能在淺水條件下更明顯地反應(yīng)波浪運(yùn)動的主要因素,波面方程用雅可比橢圓余弦函數(shù)表示為:

    (14)

    (15)

    Zt=ZC-H

    (16)

    式中,Zt—地面與波谷底端的垂直距離;ZC—地面與波峰頂?shù)拇怪本嚯x;k—橢圓積分模數(shù);Kk—k的第一類完全橢圓積分;Ek—k的第二類橢圓積分;Cn—橢圓余弦函數(shù),當(dāng)0

    1.4 孤立波浪理論

    孤立波是一種在長距離傳播中可以較好地保持波浪形態(tài)以及波浪能量的非線性波浪理論,它一般用于在淺水條件下,波浪沖擊船舶等物體破碎或產(chǎn)生極端波浪的情況。當(dāng)式(14)、(15)中k=1時,波浪為孤立波;當(dāng)k=0時,波浪為線性波。

    2 數(shù)值造波

    2.1 速度入口造波方法

    速度入口造波是提供一個目標(biāo)波的水平及垂向運(yùn)動速度給入口邊界處的水質(zhì)點(diǎn),促使水質(zhì)點(diǎn)周期性運(yùn)動進(jìn)行造波。數(shù)值波浪水槽的造波區(qū)入口處需要將邊界設(shè)為速度入口邊界,邊界入口定義的入射波形式為:

    H=Acos(kx+ky-ωεt+ε)

    (17)

    式中,kx=kcosθ;ky=ksinθ;A—振幅;ε—相位差;k—波數(shù)。

    有效波頻率定義為:

    ωε=ω+kU

    (18)

    式中,U—入射波速度;ω—波頻率。

    2.2 動量源造波方法

    該方法是通過fluent自帶的二次編譯功能在造波源區(qū)域給動量方程添加動量源項s(x,y,t),使水槽入口處水質(zhì)點(diǎn)擁有運(yùn)動速度,此時水面將在重力作用下進(jìn)行升降,從而形成波浪。其動量源s(x,y,t)可用以下公式確定:

    (19)

    (20)

    式中,u(x,y,t)—t時刻目標(biāo)波浪在坐標(biāo)(x,y)處的水質(zhì)點(diǎn)水平速度,根據(jù)線性波理論或Stokes波理論取值;u、v—x、y方向速度;Δx—水平方向網(wǎng)格長度。

    2.3 質(zhì)量源造波方法

    質(zhì)量源造波與物理實(shí)驗中的推板造波機(jī)的原理相同,將入口處固壁邊界處理成動網(wǎng)格來模擬推波板運(yùn)動。其水平運(yùn)動速度函數(shù)u(t)如下:

    (21)

    (22)

    式中,η—波面升高,根據(jù)Stokes波理論或線性波理論確定其值;ω—角頻率;C0—造波板運(yùn)動與給定波浪間的傳遞函數(shù);k—波數(shù);d—水深。

    3 數(shù)值消波

    3.1 多孔介質(zhì)消波方法

    在物理水池實(shí)驗中,經(jīng)常需要在水槽末端布置海綿等多孔介質(zhì)吸收波浪,這是為了避免反射波對造波效果產(chǎn)生影響。所以在數(shù)值水槽的末端可以同樣設(shè)定一定長度的多孔介質(zhì)區(qū)域,用以達(dá)到吸收反射波的效果,保證數(shù)值消波的合理性。

    其原理是在消波區(qū)域的動量方程中添加動量衰減的源項,從而達(dá)到消除波浪的作用,該源項表達(dá)式為:

    (23)

    式中,第一項為黏性耗散項;第二項為慣性耗散項;v—速度;1/α、C2—常數(shù)。

    考慮到該動量源過大會對工作區(qū)域來波產(chǎn)生影響,所以僅設(shè)置一個合適的慣性耗散系數(shù)C2,并將此系數(shù)沿波浪傳播方向線性遞增,則在消波區(qū)x=xi處有[8]:

    (C2)i=KP[(xj-x0)/(xe-x0)]

    (24)

    式中,x0、xe—消波區(qū)前端和尾端的橫坐標(biāo),x0

    3.2 動量阻尼消波方法

    動量阻尼消波方法的消波原理是在數(shù)值水槽末端的消波區(qū)域設(shè)定一個動量源,將運(yùn)動方向設(shè)為波浪傳播反方向,相當(dāng)于提供一個與波浪運(yùn)動方向相反的阻力項。阻尼消波的動量方程為:

    (25)

    (26)

    式中,ρ—流體相對密度;t—波浪傳播時間;u—波浪傳播的水平方向的速度分量;v—波浪傳播的豎向速度分量;p—劃分流體網(wǎng)格單元上的壓力;g—地表重力加速度;μ—消波動力系數(shù)。

    4 數(shù)值模擬

    本文基于FLUENT計算流體軟件,采用VOF法追蹤水氣交界面[9],以不可壓縮黏性流體的N-S為時均控制方程,采用k-ε湍流模型,并以PISO算法的二階迎風(fēng)格模式進(jìn)行求解計算,從而得到了合理的數(shù)值結(jié)果。

    4.1 工況條件及網(wǎng)格劃分

    根據(jù)某一大跨斜拉橋橋塔樣式建立數(shù)值水槽,并根據(jù)造波區(qū)與消波區(qū)長度[10]設(shè)定長為368m,寬度為125m,高為250m,水深為40m,選取有效波高Hs為6.4m,有效波周期為7.38s,波長為84.50m。數(shù)值水槽網(wǎng)格劃分采用六面體結(jié)構(gòu)化網(wǎng)格,在波浪傳播方向的網(wǎng)格均采用0.15m加密網(wǎng)格處理,在水槽豎向以自由水面為界采用漸變網(wǎng)格加密處理,網(wǎng)格尺寸在0.25~1m之間,為了更加精確反應(yīng)橋墩壁面上所受到波浪力作用大小,對橋墩流域處的縱向網(wǎng)格設(shè)置邊界層加密到0.1m。最后劃分了2140000左右個網(wǎng)格單元,網(wǎng)格質(zhì)量均在0.6以上符合計算要求,網(wǎng)格劃分如圖1所示?;?種不同波浪理論的造波結(jié)果,線性波與stokes波的波面時程對比如圖2所示,水質(zhì)點(diǎn)速度對比如圖3所示。

    圖1 網(wǎng)格劃分圖

    圖2 波面時程

    圖3 水質(zhì)點(diǎn)速度

    從圖2中不難看出,由于Stokes波在加入波陡參數(shù)后,其理論下的波浪波峰較線性波更陡,波谷較線性波更平坦。同時從圖3中看出線性波與Stokes波水質(zhì)點(diǎn)水平速度與豎向速度基本保持一致。

    4.2 數(shù)值結(jié)果對比

    基于fluent造波方法中,源造波及推板造波需要使用動網(wǎng)格技術(shù),并且波高隨水域衰減快[11],本文采用造波效果更好的速度入口造波,消波區(qū)通過在fluent軟件設(shè)置自定義UDF,編寫對應(yīng)數(shù)值阻尼源項宏的消波代碼,最后得到數(shù)值結(jié)果與理論值的對比結(jié)果如圖4所示。

    圖4 數(shù)值模擬與理論計算對比

    由圖4可知,由于采用速度入口造波,在距入射邊界處100m的前15s,波浪尚在傳遞,此處水面呈發(fā)展?fàn)顟B(tài),在23s處和理論值得到了較好的擬合,并保持穩(wěn)定,誤差較小,波高衰弱較小,得證數(shù)值模擬效果良好。

    5 結(jié)論

    本文首先介紹當(dāng)前主要的4種波浪理論,再系統(tǒng)地總結(jié)了當(dāng)下數(shù)值模擬波浪研究主要采用的造波及消波方法,通過建立一個實(shí)況斜拉橋橋塔的數(shù)值水槽,得到了波浪數(shù)值解,并與理論界對比,分析得到以下結(jié)論:

    (1)當(dāng)前主要波浪理論為線性波浪理論,非線性斯托克斯波浪理論,橢圓余弦波及孤立波。通過fluent軟件發(fā)現(xiàn)stokes波波峰較線性波更陡,波谷較線性波更平坦,但線性波與Stokes波水質(zhì)點(diǎn)水平速度與豎向速度基本保持一致。

    (2)當(dāng)前數(shù)組模擬波浪研究主要采用的造波方法為:速度入口造波方法、動量源造波方法、質(zhì)量源造波方法等;消波方法有:多孔介質(zhì)消波方法、動量阻尼消波方法等。

    (3)基于fluent計算流體軟件,采用stokes波,速度入口造波方法及動量源消波方法得到的數(shù)值模擬解與理論計算解擬合較好,驗證了數(shù)值造波的可靠性。

    猜你喜歡
    造波消波動量
    主體結(jié)構(gòu)荷載可控的新型組合式防護(hù)結(jié)構(gòu)(Ⅱ):影響因素及設(shè)計理念
    爆炸與沖擊(2025年1期)2025-02-19 00:00:00
    動量守恒定律在三個物體系中的應(yīng)用
    基于模糊PID控制的主動造波系統(tǒng)研究與應(yīng)用
    分層流水槽箱內(nèi)垂蕩板式內(nèi)波造波模擬研究
    應(yīng)用動量守恒定律解題之秘訣
    動量相關(guān)知識的理解和應(yīng)用
    循環(huán)水槽多層孔板消波裝置開發(fā)及消波特性數(shù)值模擬
    基于潛堤地形上的波浪傳播模擬
    科技資訊(2017年19期)2017-08-08 08:39:37
    基于最佳逼近理論的主動吸收造波算法研究
    多層變孔徑傾斜孔板式消波裝置的試驗研究
    日本av手机在线免费观看| 18禁在线播放成人免费| 国产免费福利视频在线观看| 久久久久久久久久成人| 亚洲av中文字字幕乱码综合| 日日啪夜夜爽| 在线亚洲精品国产二区图片欧美 | 熟妇人妻不卡中文字幕| 国产视频内射| 97精品久久久久久久久久精品| 午夜免费鲁丝| 免费不卡的大黄色大毛片视频在线观看| av免费观看日本| 久久久午夜欧美精品| 精品人妻偷拍中文字幕| 纵有疾风起免费观看全集完整版| 亚洲第一区二区三区不卡| 在线天堂最新版资源| 欧美日韩综合久久久久久| 人人妻人人澡人人爽人人夜夜| 99国产精品免费福利视频| 免费看不卡的av| 成年女人在线观看亚洲视频| 少妇精品久久久久久久| 丝袜喷水一区| 尤物成人国产欧美一区二区三区| 国产色爽女视频免费观看| 中文字幕制服av| 成年免费大片在线观看| 成年免费大片在线观看| 国产无遮挡羞羞视频在线观看| 妹子高潮喷水视频| 国产淫片久久久久久久久| 国产精品99久久99久久久不卡 | 国产乱来视频区| 免费在线观看成人毛片| 99久国产av精品国产电影| 热99国产精品久久久久久7| 国内精品宾馆在线| 美女高潮的动态| 少妇人妻精品综合一区二区| 天天躁日日操中文字幕| 好男人视频免费观看在线| 欧美日韩国产mv在线观看视频 | 免费大片黄手机在线观看| 国产探花极品一区二区| 赤兔流量卡办理| 内射极品少妇av片p| 色综合色国产| 国产在线免费精品| 国产黄片美女视频| 久久国产精品大桥未久av | 人妻制服诱惑在线中文字幕| 久久精品国产亚洲av涩爱| 成人高潮视频无遮挡免费网站| 色综合色国产| 精品国产一区二区三区久久久樱花 | 久久韩国三级中文字幕| 偷拍熟女少妇极品色| 免费观看无遮挡的男女| 国产精品久久久久久av不卡| 免费观看无遮挡的男女| 亚洲欧美清纯卡通| 一本一本综合久久| 亚洲国产精品成人久久小说| 麻豆精品久久久久久蜜桃| 成人二区视频| 国产亚洲av片在线观看秒播厂| 亚洲三级黄色毛片| 日韩国内少妇激情av| 欧美日韩一区二区视频在线观看视频在线| 欧美高清性xxxxhd video| 亚洲精品一区蜜桃| 国产色婷婷99| 久久久精品免费免费高清| 久久鲁丝午夜福利片| 在线精品无人区一区二区三 | 97精品久久久久久久久久精品| 久久久久久人妻| 18禁在线播放成人免费| 人人妻人人添人人爽欧美一区卜 | 人妻系列 视频| 久久久久久人妻| 少妇的逼好多水| 国产真实伦视频高清在线观看| 肉色欧美久久久久久久蜜桃| 精品国产一区二区三区久久久樱花 | av免费在线看不卡| 美女国产视频在线观看| av又黄又爽大尺度在线免费看| 成人国产av品久久久| 美女中出高潮动态图| 精品一区二区三卡| 国产成人午夜福利电影在线观看| 国产毛片在线视频| 国产精品.久久久| 色网站视频免费| 大话2 男鬼变身卡| 国产爱豆传媒在线观看| 日本与韩国留学比较| 免费播放大片免费观看视频在线观看| 国产精品一及| 九色成人免费人妻av| 青春草亚洲视频在线观看| 一二三四中文在线观看免费高清| 国产亚洲精品久久久com| 熟女人妻精品中文字幕| 亚洲av综合色区一区| 中文字幕免费在线视频6| 亚洲第一av免费看| tube8黄色片| 久久婷婷青草| 国产成人一区二区在线| 天堂中文最新版在线下载| 亚洲美女黄色视频免费看| 自拍偷自拍亚洲精品老妇| av国产免费在线观看| 最近最新中文字幕免费大全7| 国产亚洲91精品色在线| 精品久久久久久久久av| av在线播放精品| 18禁裸乳无遮挡免费网站照片| 简卡轻食公司| 三级国产精品片| 久热久热在线精品观看| 国产久久久一区二区三区| 中文字幕免费在线视频6| 国产色婷婷99| 久久久久久久久久成人| 亚洲电影在线观看av| 亚洲av中文av极速乱| 欧美日韩视频精品一区| 亚洲av成人精品一二三区| 亚洲国产av新网站| 亚洲内射少妇av| 久久6这里有精品| 在线免费十八禁| 干丝袜人妻中文字幕| 最近最新中文字幕大全电影3| 大又大粗又爽又黄少妇毛片口| 97超碰精品成人国产| 九草在线视频观看| 大话2 男鬼变身卡| 亚洲真实伦在线观看| 国产精品秋霞免费鲁丝片| 亚洲精品色激情综合| av又黄又爽大尺度在线免费看| 一级毛片久久久久久久久女| 一本久久精品| 建设人人有责人人尽责人人享有的 | 毛片一级片免费看久久久久| 51国产日韩欧美| 中文字幕亚洲精品专区| 午夜免费男女啪啪视频观看| 少妇高潮的动态图| 亚洲国产高清在线一区二区三| 亚洲天堂av无毛| 免费久久久久久久精品成人欧美视频 | 久久久久国产网址| 人体艺术视频欧美日本| 国产一区二区三区av在线| 又粗又硬又长又爽又黄的视频| 久久6这里有精品| 亚洲人与动物交配视频| 老女人水多毛片| 亚洲国产精品999| 人体艺术视频欧美日本| 日本av免费视频播放| av一本久久久久| 日本欧美国产在线视频| a级毛色黄片| av在线app专区| 99re6热这里在线精品视频| 欧美日韩在线观看h| 亚洲国产av新网站| 久久久久人妻精品一区果冻| 夜夜骑夜夜射夜夜干| 久久久久久久大尺度免费视频| 久久久精品免费免费高清| 熟女电影av网| 亚洲成人中文字幕在线播放| 久久久成人免费电影| 国产精品三级大全| 久久6这里有精品| 建设人人有责人人尽责人人享有的 | 九九爱精品视频在线观看| 久久亚洲国产成人精品v| 免费观看在线日韩| 哪个播放器可以免费观看大片| 一本色道久久久久久精品综合| 18禁在线播放成人免费| 在线 av 中文字幕| 亚洲av在线观看美女高潮| 亚洲国产精品999| 寂寞人妻少妇视频99o| 亚洲av国产av综合av卡| 亚洲伊人久久精品综合| 青春草亚洲视频在线观看| 男人和女人高潮做爰伦理| av在线老鸭窝| 伦理电影大哥的女人| 久久热精品热| 少妇人妻久久综合中文| 天美传媒精品一区二区| 狂野欧美白嫩少妇大欣赏| 美女中出高潮动态图| 久久国产精品大桥未久av | 22中文网久久字幕| 中国三级夫妇交换| 国产av国产精品国产| 久久久久性生活片| 国产一区有黄有色的免费视频| 欧美人与善性xxx| 国产乱人视频| 黑人高潮一二区| 久久国产亚洲av麻豆专区| 亚洲av免费高清在线观看| 亚洲美女搞黄在线观看| 一个人看视频在线观看www免费| 国产爱豆传媒在线观看| 大片电影免费在线观看免费| 国产精品爽爽va在线观看网站| 我的女老师完整版在线观看| av国产精品久久久久影院| 午夜激情福利司机影院| 午夜精品国产一区二区电影| 另类亚洲欧美激情| 纵有疾风起免费观看全集完整版| 爱豆传媒免费全集在线观看| 好男人视频免费观看在线| 免费看日本二区| 秋霞伦理黄片| 色网站视频免费| 亚洲欧美成人综合另类久久久| 综合色丁香网| 青春草国产在线视频| 日韩视频在线欧美| 亚洲av二区三区四区| 久久精品久久久久久噜噜老黄| 国产精品女同一区二区软件| 亚洲欧美日韩另类电影网站 | 久久毛片免费看一区二区三区| 国产女主播在线喷水免费视频网站| 亚州av有码| 女性被躁到高潮视频| 成人综合一区亚洲| 免费观看在线日韩| 各种免费的搞黄视频| 国产免费一区二区三区四区乱码| 99久久人妻综合| 九草在线视频观看| 欧美日韩精品成人综合77777| 亚洲精品久久午夜乱码| 国产无遮挡羞羞视频在线观看| 国产毛片在线视频| 午夜福利在线在线| 亚洲精品国产色婷婷电影| 久久精品国产亚洲网站| 国产伦精品一区二区三区四那| 高清黄色对白视频在线免费看 | av国产免费在线观看| 嫩草影院新地址| 国产乱来视频区| 99re6热这里在线精品视频| 欧美日韩视频精品一区| av国产精品久久久久影院| 色视频www国产| 欧美高清成人免费视频www| 久久久久久九九精品二区国产| 精品一区二区免费观看| 六月丁香七月| av黄色大香蕉| 精品国产露脸久久av麻豆| 国产成人免费观看mmmm| 国产免费一区二区三区四区乱码| 乱码一卡2卡4卡精品| 亚洲av日韩在线播放| xxx大片免费视频| 美女cb高潮喷水在线观看| 日日啪夜夜爽| 亚洲国产高清在线一区二区三| 精品久久久久久电影网| 能在线免费看毛片的网站| 亚洲av成人精品一区久久| 日韩人妻高清精品专区| 久久影院123| 免费看日本二区| 色综合色国产| 少妇人妻久久综合中文| 性色avwww在线观看| 国产精品麻豆人妻色哟哟久久| 我的老师免费观看完整版| 免费少妇av软件| 久久av网站| 最黄视频免费看| 亚洲综合精品二区| 国产日韩欧美在线精品| 在线看a的网站| 精品久久国产蜜桃| 国产成人aa在线观看| 亚洲欧美清纯卡通| 麻豆成人午夜福利视频| 亚洲人成网站在线观看播放| av国产免费在线观看| 日韩免费高清中文字幕av| 成人一区二区视频在线观看| 男人舔奶头视频| 一级片'在线观看视频| 亚洲av免费高清在线观看| 日韩电影二区| 亚洲精品视频女| av国产精品久久久久影院| 国产无遮挡羞羞视频在线观看| 尾随美女入室| 婷婷色综合大香蕉| videossex国产| 久久久精品免费免费高清| 2022亚洲国产成人精品| 日韩成人av中文字幕在线观看| 亚洲国产成人一精品久久久| 看免费成人av毛片| 日日啪夜夜撸| 黄色怎么调成土黄色| 极品教师在线视频| 伦理电影大哥的女人| 一区二区av电影网| 久久午夜福利片| 国产高清三级在线| 亚洲电影在线观看av| 高清日韩中文字幕在线| 亚洲av免费高清在线观看| 亚洲美女搞黄在线观看| 性高湖久久久久久久久免费观看| 国产精品av视频在线免费观看| 丰满少妇做爰视频| 国产色爽女视频免费观看| 男的添女的下面高潮视频| 亚洲av福利一区| 中文乱码字字幕精品一区二区三区| 在线观看免费视频网站a站| 日韩一本色道免费dvd| 久久久久人妻精品一区果冻| 免费看日本二区| 久久鲁丝午夜福利片| 精华霜和精华液先用哪个| 国产免费一区二区三区四区乱码| 国产乱人偷精品视频| av在线蜜桃| 黑人猛操日本美女一级片| 超碰97精品在线观看| 国产精品一区二区在线观看99| 亚洲精品,欧美精品| 伦理电影免费视频| 性色av一级| 久久鲁丝午夜福利片| 亚洲国产av新网站| 高清欧美精品videossex| 男女无遮挡免费网站观看| 永久免费av网站大全| 国产精品三级大全| 在线免费观看不下载黄p国产| 国产在线男女| 日产精品乱码卡一卡2卡三| 亚洲精品久久午夜乱码| 丝瓜视频免费看黄片| 99精国产麻豆久久婷婷| 国产精品国产三级国产av玫瑰| 热99国产精品久久久久久7| 日韩视频在线欧美| 亚洲精品自拍成人| 搡老乐熟女国产| 极品教师在线视频| 夫妻性生交免费视频一级片| 我的女老师完整版在线观看| 男人和女人高潮做爰伦理| 成人综合一区亚洲| freevideosex欧美| 99九九线精品视频在线观看视频| 午夜精品国产一区二区电影| 国产高清不卡午夜福利| 亚洲精品色激情综合| 丝袜脚勾引网站| 伦精品一区二区三区| 久久精品国产自在天天线| 女人十人毛片免费观看3o分钟| 亚洲精品一区蜜桃| 国产爱豆传媒在线观看| 国国产精品蜜臀av免费| 一级黄片播放器| 国产成人freesex在线| 日韩大片免费观看网站| 亚洲国产精品成人久久小说| 熟女av电影| 日韩不卡一区二区三区视频在线| 插阴视频在线观看视频| 免费少妇av软件| 国内少妇人妻偷人精品xxx网站| 18禁在线播放成人免费| 亚洲成人一二三区av| 夫妻性生交免费视频一级片| 爱豆传媒免费全集在线观看| 99热这里只有是精品在线观看| 九草在线视频观看| 国产成人免费观看mmmm| 在现免费观看毛片| 免费看不卡的av| 国产一区二区三区av在线| 一二三四中文在线观看免费高清| 国产乱人偷精品视频| 日韩欧美 国产精品| 精品少妇久久久久久888优播| 亚洲国产精品国产精品| 亚洲精品自拍成人| 高清不卡的av网站| 直男gayav资源| 欧美三级亚洲精品| 国产亚洲91精品色在线| 美女高潮的动态| 伊人久久精品亚洲午夜| 日本vs欧美在线观看视频 | 亚洲国产最新在线播放| 国产伦精品一区二区三区四那| 在线观看一区二区三区| 午夜老司机福利剧场| av专区在线播放| 最新中文字幕久久久久| 国语对白做爰xxxⅹ性视频网站| 亚洲av中文字字幕乱码综合| 哪个播放器可以免费观看大片| 视频中文字幕在线观看| 高清毛片免费看| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产a三级三级三级| 免费观看在线日韩| 久久综合国产亚洲精品| 极品少妇高潮喷水抽搐| 国产成人精品久久久久久| 亚洲国产精品一区三区| 极品教师在线视频| 肉色欧美久久久久久久蜜桃| 精品久久国产蜜桃| 欧美xxxx性猛交bbbb| 麻豆成人av视频| 精品国产三级普通话版| 麻豆成人午夜福利视频| 国产在线免费精品| 麻豆国产97在线/欧美| 中文字幕免费在线视频6| 91久久精品电影网| 综合色丁香网| 尤物成人国产欧美一区二区三区| 汤姆久久久久久久影院中文字幕| 亚洲va在线va天堂va国产| 欧美激情国产日韩精品一区| 日韩大片免费观看网站| 青春草国产在线视频| 国产精品一区二区性色av| 欧美一区二区亚洲| 国产久久久一区二区三区| 中文字幕av成人在线电影| 国语对白做爰xxxⅹ性视频网站| 亚洲国产精品999| 99热网站在线观看| 免费黄色在线免费观看| 国产伦理片在线播放av一区| 美女主播在线视频| 亚洲精品国产av成人精品| 在线观看人妻少妇| 日本黄大片高清| 成人特级av手机在线观看| 狠狠精品人妻久久久久久综合| 亚洲自偷自拍三级| 男男h啪啪无遮挡| 97在线人人人人妻| 亚洲av日韩在线播放| 丰满乱子伦码专区| 欧美丝袜亚洲另类| 国产免费福利视频在线观看| 久久韩国三级中文字幕| 午夜免费观看性视频| 日日撸夜夜添| 国产精品一二三区在线看| 亚洲欧美一区二区三区黑人 | 国产成人免费无遮挡视频| 久久亚洲国产成人精品v| av又黄又爽大尺度在线免费看| 国产美女午夜福利| 大香蕉久久网| 国国产精品蜜臀av免费| 观看免费一级毛片| 日本av免费视频播放| 性色avwww在线观看| 日韩欧美一区视频在线观看 | 夜夜骑夜夜射夜夜干| 中文字幕人妻熟人妻熟丝袜美| 三级国产精品欧美在线观看| 亚洲精品久久午夜乱码| 极品少妇高潮喷水抽搐| 精品人妻偷拍中文字幕| 欧美成人a在线观看| 赤兔流量卡办理| 成人国产麻豆网| 国产片特级美女逼逼视频| 日本黄色日本黄色录像| 有码 亚洲区| 欧美成人午夜免费资源| 久久久午夜欧美精品| 大香蕉97超碰在线| 女人久久www免费人成看片| 色哟哟·www| 尤物成人国产欧美一区二区三区| 欧美精品一区二区免费开放| 国产精品久久久久久精品电影小说 | 国产av码专区亚洲av| 校园人妻丝袜中文字幕| 欧美高清成人免费视频www| 在线观看三级黄色| 国产久久久一区二区三区| 男女边摸边吃奶| 大香蕉97超碰在线| 大香蕉久久网| 18禁在线播放成人免费| 日本欧美视频一区| 97超视频在线观看视频| 久久久久久人妻| 久久国产乱子免费精品| 少妇 在线观看| 欧美精品人与动牲交sv欧美| 欧美日韩国产mv在线观看视频 | a级毛片免费高清观看在线播放| 老女人水多毛片| 少妇人妻 视频| 最近2019中文字幕mv第一页| 肉色欧美久久久久久久蜜桃| 综合色丁香网| 18+在线观看网站| 亚州av有码| 国产乱人偷精品视频| 久久精品久久久久久噜噜老黄| 国产视频内射| 国内少妇人妻偷人精品xxx网站| 2021少妇久久久久久久久久久| 少妇高潮的动态图| 亚洲av电影在线观看一区二区三区| 男人添女人高潮全过程视频| 亚洲综合色惰| 女的被弄到高潮叫床怎么办| 国产精品偷伦视频观看了| 亚洲精品日本国产第一区| 全区人妻精品视频| 卡戴珊不雅视频在线播放| 国产伦理片在线播放av一区| 大片电影免费在线观看免费| 涩涩av久久男人的天堂| 色婷婷av一区二区三区视频| 熟妇人妻不卡中文字幕| av免费观看日本| 嘟嘟电影网在线观看| 少妇被粗大猛烈的视频| 国产精品伦人一区二区| 亚洲真实伦在线观看| 欧美xxxx性猛交bbbb| 欧美xxⅹ黑人| 亚洲国产最新在线播放| 性高湖久久久久久久久免费观看| 免费观看性生交大片5| 国产伦在线观看视频一区| 国产精品一区二区在线不卡| 日本色播在线视频| 人妻系列 视频| 国产乱人视频| 久久精品夜色国产| 午夜福利高清视频| 成年人午夜在线观看视频| 国产免费福利视频在线观看| 2021少妇久久久久久久久久久| 全区人妻精品视频| 亚洲精品国产av成人精品| 肉色欧美久久久久久久蜜桃| 插阴视频在线观看视频| 欧美区成人在线视频| 王馨瑶露胸无遮挡在线观看| 人人妻人人澡人人爽人人夜夜| 国产成人免费无遮挡视频| 在线精品无人区一区二区三 | 蜜桃在线观看..| 精品一区在线观看国产| 人人妻人人看人人澡| 久久影院123| 麻豆成人av视频| 人人妻人人添人人爽欧美一区卜 | 自拍偷自拍亚洲精品老妇| 欧美3d第一页| 日本色播在线视频| 极品教师在线视频| 国产一区二区在线观看日韩| 国产欧美日韩一区二区三区在线 | 在线亚洲精品国产二区图片欧美 | 久久久久久久亚洲中文字幕| 日日撸夜夜添| 又黄又爽又刺激的免费视频.| 一本色道久久久久久精品综合| 国产高清三级在线| 亚洲内射少妇av| 亚洲人成网站在线观看播放| 国产视频首页在线观看| 亚洲,欧美,日韩| 高清不卡的av网站| 夫妻性生交免费视频一级片| 好男人视频免费观看在线| 99精国产麻豆久久婷婷| 大香蕉97超碰在线| 亚洲国产高清在线一区二区三| 亚洲va在线va天堂va国产| 免费观看在线日韩|