• <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
    基于最佳逼近理論的主動吸收造波算法研究
    多層變孔徑傾斜孔板式消波裝置的試驗研究
    18+在线观看网站| 99国产综合亚洲精品| 久久久久精品国产欧美久久久| 日本五十路高清| 国产精品嫩草影院av在线观看 | 国产黄a三级三级三级人| 五月玫瑰六月丁香| 亚洲av成人av| 日韩欧美三级三区| 国产精品,欧美在线| 国产精品,欧美在线| АⅤ资源中文在线天堂| 麻豆国产97在线/欧美| 国产一区二区三区在线臀色熟女| 欧美性感艳星| 午夜老司机福利剧场| e午夜精品久久久久久久| 国产激情偷乱视频一区二区| 91字幕亚洲| 亚洲熟妇中文字幕五十中出| 久9热在线精品视频| 两个人看的免费小视频| 国内精品美女久久久久久| 欧美午夜高清在线| 日本一二三区视频观看| 怎么达到女性高潮| 久久精品人妻少妇| 叶爱在线成人免费视频播放| www国产在线视频色| 中文资源天堂在线| 女同久久另类99精品国产91| 夜夜看夜夜爽夜夜摸| 又黄又粗又硬又大视频| 好男人在线观看高清免费视频| 亚洲欧美日韩卡通动漫| 伊人久久大香线蕉亚洲五| 国产69精品久久久久777片| 久久久久久久久大av| av福利片在线观看| 日本黄色片子视频| 国产伦精品一区二区三区视频9 | 中文字幕av成人在线电影| 国产午夜福利久久久久久| 国产高清视频在线播放一区| 又粗又爽又猛毛片免费看| 国产真实伦视频高清在线观看 | 少妇的逼水好多| 久久精品国产亚洲av涩爱 | 99热精品在线国产| 国产午夜福利久久久久久| 在线免费观看的www视频| 中国美女看黄片| 内地一区二区视频在线| 一区二区三区激情视频| 琪琪午夜伦伦电影理论片6080| 精品熟女少妇八av免费久了| 搡老熟女国产l中国老女人| 日本与韩国留学比较| 99国产综合亚洲精品| 亚洲成a人片在线一区二区| 国产在视频线在精品| 久久人妻av系列| 国内毛片毛片毛片毛片毛片| 久久精品夜夜夜夜夜久久蜜豆| 内射极品少妇av片p| 亚洲精品456在线播放app | 免费人成在线观看视频色| 18美女黄网站色大片免费观看| 精华霜和精华液先用哪个| 日日摸夜夜添夜夜添小说| 亚洲av成人精品一区久久| 日韩大尺度精品在线看网址| 日韩欧美在线乱码| 网址你懂的国产日韩在线| 亚洲专区国产一区二区| 久久欧美精品欧美久久欧美| www日本黄色视频网| 精品乱码久久久久久99久播| 国产亚洲欧美在线一区二区| 天美传媒精品一区二区| 午夜福利在线观看吧| 国产亚洲av嫩草精品影院| 亚洲欧美精品综合久久99| 久久久久国产精品人妻aⅴ院| 国产精品香港三级国产av潘金莲| 欧美成人一区二区免费高清观看| 天堂网av新在线| 欧美一区二区精品小视频在线| 国产一区二区在线av高清观看| 国产成人aa在线观看| 日韩欧美在线乱码| 久久这里只有精品中国| 亚洲欧美日韩高清专用| 亚洲国产日韩欧美精品在线观看 | 国产精品av视频在线免费观看| 欧美成人a在线观看| 亚洲第一电影网av| 欧美一区二区精品小视频在线| 色综合站精品国产| 校园春色视频在线观看| 亚洲精品在线美女| 成年女人永久免费观看视频| 久久精品国产99精品国产亚洲性色| 久久精品国产亚洲av涩爱 | 黄色女人牲交| 一区二区三区激情视频| 亚洲欧美日韩高清在线视频| 99riav亚洲国产免费| www.999成人在线观看| 国产乱人视频| 国产午夜精品论理片| 亚洲精品456在线播放app | 人人妻人人澡欧美一区二区| 久久欧美精品欧美久久欧美| 丰满人妻一区二区三区视频av | 免费无遮挡裸体视频| 国产精品av视频在线免费观看| 精品一区二区三区视频在线 | 极品教师在线免费播放| av黄色大香蕉| 久久久久久人人人人人| 噜噜噜噜噜久久久久久91| 国产精品亚洲av一区麻豆| 亚洲最大成人中文| 日韩欧美三级三区| 久久久国产成人精品二区| 色老头精品视频在线观看| 他把我摸到了高潮在线观看| 久久性视频一级片| 久久精品国产亚洲av涩爱 | 日本一二三区视频观看| 午夜视频国产福利| 国产真人三级小视频在线观看| 欧美成人一区二区免费高清观看| 久久久久久久久中文| 啦啦啦免费观看视频1| 免费看十八禁软件| 日韩大尺度精品在线看网址| 美女大奶头视频| 日本 欧美在线| 国产精品一区二区免费欧美| 精品电影一区二区在线| 国产黄色小视频在线观看| 亚洲男人的天堂狠狠| 亚洲成人免费电影在线观看| 亚洲av电影在线进入| 婷婷丁香在线五月| 午夜福利高清视频| 夜夜爽天天搞| 中文字幕久久专区| 色噜噜av男人的天堂激情| 嫁个100分男人电影在线观看| 窝窝影院91人妻| 亚洲乱码一区二区免费版| 亚洲精华国产精华精| 精品一区二区三区av网在线观看| 久久欧美精品欧美久久欧美| 夜夜爽天天搞| 亚洲中文日韩欧美视频| 免费av毛片视频| 欧美色视频一区免费| 亚洲成av人片在线播放无| 国产一区在线观看成人免费| 麻豆成人午夜福利视频| 久久久久国产精品人妻aⅴ院| 久久久久久久久久黄片| 又爽又黄无遮挡网站| 又紧又爽又黄一区二区| 九色成人免费人妻av| 美女cb高潮喷水在线观看| 熟女人妻精品中文字幕| 在线观看美女被高潮喷水网站 | 久久久久久久久久黄片| 日韩精品青青久久久久久| 成年版毛片免费区| 人人妻人人看人人澡| 成人精品一区二区免费| 国产高清三级在线| 特大巨黑吊av在线直播| 在线观看日韩欧美| 国产亚洲精品久久久com| 国产v大片淫在线免费观看| 免费在线观看日本一区| 亚洲专区国产一区二区| 男人舔女人下体高潮全视频| 少妇裸体淫交视频免费看高清| 狂野欧美白嫩少妇大欣赏| 级片在线观看| 国产毛片a区久久久久| 97碰自拍视频| 少妇的逼好多水| 性色avwww在线观看| 首页视频小说图片口味搜索| 欧美一级a爱片免费观看看| 精品一区二区三区视频在线观看免费| 少妇裸体淫交视频免费看高清| 最新美女视频免费是黄的| 人人妻,人人澡人人爽秒播| 日本 欧美在线| 真人做人爱边吃奶动态| 精品国产美女av久久久久小说| 国产美女午夜福利| 99视频精品全部免费 在线| 久久久国产精品麻豆| 日本在线视频免费播放| 欧美av亚洲av综合av国产av| 国产真实乱freesex| 成人国产综合亚洲| 极品教师在线免费播放| 99热这里只有是精品50| 婷婷六月久久综合丁香| 国产69精品久久久久777片| 国产爱豆传媒在线观看| 十八禁人妻一区二区| 在线国产一区二区在线| 日本精品一区二区三区蜜桃| svipshipincom国产片| 好男人电影高清在线观看| 男女视频在线观看网站免费| 在线观看舔阴道视频| 亚洲片人在线观看| 中国美女看黄片| 成年女人毛片免费观看观看9| 国内久久婷婷六月综合欲色啪| 少妇人妻精品综合一区二区 | 日韩国内少妇激情av| 欧美日韩国产亚洲二区| 九色国产91popny在线| 人妻久久中文字幕网| 中文字幕精品亚洲无线码一区| 久久午夜亚洲精品久久| 狂野欧美激情性xxxx| 波多野结衣巨乳人妻| 久久精品人妻少妇| 亚洲性夜色夜夜综合| 国产69精品久久久久777片| 搞女人的毛片| 午夜a级毛片| 国产精品 国内视频| 精品久久久久久久久久久久久| 日韩中文字幕欧美一区二区| 色吧在线观看| 最后的刺客免费高清国语| 淫妇啪啪啪对白视频| 亚洲无线观看免费| 成人av一区二区三区在线看| 国产欧美日韩一区二区三| 母亲3免费完整高清在线观看| 在线国产一区二区在线| 19禁男女啪啪无遮挡网站| 校园春色视频在线观看| 国产真实伦视频高清在线观看 | 欧美xxxx黑人xx丫x性爽| 国产亚洲精品一区二区www| 国产精品一区二区三区四区久久| 久久久久精品国产欧美久久久| 亚洲真实伦在线观看| 欧美乱妇无乱码| 少妇高潮的动态图| 51午夜福利影视在线观看| 国产av麻豆久久久久久久| 久久这里只有精品中国| 亚洲av不卡在线观看| 大型黄色视频在线免费观看| 国产淫片久久久久久久久 | 成人国产一区最新在线观看| av女优亚洲男人天堂| 又紧又爽又黄一区二区| 男女视频在线观看网站免费| 久久久久国产精品人妻aⅴ院| 午夜福利免费观看在线| 精品国内亚洲2022精品成人| 日韩欧美三级三区| 内地一区二区视频在线| 国产不卡一卡二| 午夜激情福利司机影院| 日韩人妻高清精品专区| 欧美日韩一级在线毛片| 老汉色∧v一级毛片| 国产精品久久视频播放| 国产精品99久久99久久久不卡| 俄罗斯特黄特色一大片| 欧美一级a爱片免费观看看| 国产免费av片在线观看野外av| 首页视频小说图片口味搜索| 免费看a级黄色片| 高清毛片免费观看视频网站| 日日夜夜操网爽| 国产99白浆流出| 久久精品影院6| 十八禁网站免费在线| 亚洲av一区综合| 成人18禁在线播放| 久久伊人香网站| 亚洲在线自拍视频| 午夜免费男女啪啪视频观看 | 国产精品国产高清国产av| 日韩精品中文字幕看吧| 国语自产精品视频在线第100页| 波野结衣二区三区在线 | 长腿黑丝高跟| 少妇人妻精品综合一区二区 | tocl精华| 看黄色毛片网站| 国产成人系列免费观看| 中国美女看黄片| 日日摸夜夜添夜夜添小说| 亚洲国产精品sss在线观看| 色精品久久人妻99蜜桃| 窝窝影院91人妻| 五月伊人婷婷丁香| 亚洲人成电影免费在线| 欧美一区二区精品小视频在线| 国产亚洲av嫩草精品影院| 亚洲专区国产一区二区| 一区二区三区免费毛片| av天堂中文字幕网| 中国美女看黄片| 白带黄色成豆腐渣| 老鸭窝网址在线观看| 国产成+人综合+亚洲专区| 欧美一区二区国产精品久久精品| 波多野结衣巨乳人妻| 哪里可以看免费的av片| 看免费av毛片| 免费高清视频大片| 色综合欧美亚洲国产小说| 国产精品久久久久久久电影 | 蜜桃亚洲精品一区二区三区| 91麻豆av在线| 久久精品国产亚洲av涩爱 | 国产亚洲精品久久久com| 国产真实乱freesex| 国产亚洲欧美98| 亚洲成a人片在线一区二区| 51国产日韩欧美| 人妻丰满熟妇av一区二区三区| 国产高潮美女av| 在线十欧美十亚洲十日本专区| 亚洲乱码一区二区免费版| 欧美成人一区二区免费高清观看| 国产欧美日韩精品一区二区| 国产一区在线观看成人免费| 99在线视频只有这里精品首页| 国产免费男女视频| 久久精品国产清高在天天线| 午夜激情福利司机影院| 亚洲av二区三区四区| 天堂动漫精品| 内地一区二区视频在线| 久久99热这里只有精品18| 99热这里只有精品一区| 久久这里只有精品中国| 午夜久久久久精精品| 久久婷婷人人爽人人干人人爱| 亚洲av电影不卡..在线观看| 搡老岳熟女国产| 成年女人毛片免费观看观看9| 国产午夜精品久久久久久一区二区三区 | 国产成人av教育| 久久久久久人人人人人| 精品无人区乱码1区二区| 日韩欧美国产一区二区入口| 最新在线观看一区二区三区| 久久国产乱子伦精品免费另类| 夜夜夜夜夜久久久久| 最好的美女福利视频网| 成人鲁丝片一二三区免费| 亚洲自拍偷在线| 一级作爱视频免费观看| 成人高潮视频无遮挡免费网站| 国产一区在线观看成人免费| 久久久久国内视频| 老司机午夜福利在线观看视频| 99久久99久久久精品蜜桃| 午夜免费观看网址| 久久6这里有精品| 夜夜夜夜夜久久久久| 脱女人内裤的视频| 久久久久久久久中文| 色老头精品视频在线观看| 午夜免费男女啪啪视频观看 | 欧美日韩黄片免| 啦啦啦免费观看视频1| 免费搜索国产男女视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 18+在线观看网站| 91麻豆精品激情在线观看国产| 最新在线观看一区二区三区| 国内精品美女久久久久久| 日韩av在线大香蕉| 可以在线观看毛片的网站| 国产主播在线观看一区二区| 国内精品久久久久久久电影| 精品不卡国产一区二区三区| 叶爱在线成人免费视频播放| 人妻丰满熟妇av一区二区三区| 少妇丰满av| 老司机午夜十八禁免费视频| 国产99白浆流出| 日本 欧美在线| 精品不卡国产一区二区三区| 日韩亚洲欧美综合| 精品人妻偷拍中文字幕| 两人在一起打扑克的视频| 亚洲乱码一区二区免费版| 蜜桃亚洲精品一区二区三区| 黄色女人牲交| 亚洲久久久久久中文字幕| 日韩 欧美 亚洲 中文字幕| 欧美另类亚洲清纯唯美| 亚洲,欧美精品.| 国产又黄又爽又无遮挡在线| 亚洲欧美日韩无卡精品| 18美女黄网站色大片免费观看| 亚洲欧美日韩高清在线视频| av福利片在线观看| 久久久久久久久大av| 熟女人妻精品中文字幕| 高潮久久久久久久久久久不卡| 亚洲 欧美 日韩 在线 免费| 久久精品国产自在天天线| 国产黄色小视频在线观看| 免费在线观看成人毛片| 夜夜躁狠狠躁天天躁| 亚洲成人免费电影在线观看| 精品久久久久久,| 少妇的逼好多水| 一进一出抽搐动态| 真人做人爱边吃奶动态| 亚洲成人精品中文字幕电影| 日韩亚洲欧美综合| 女同久久另类99精品国产91| 久9热在线精品视频| 亚洲精品亚洲一区二区| 欧美成人免费av一区二区三区| 男女做爰动态图高潮gif福利片| 国产老妇女一区| 最近最新免费中文字幕在线| 身体一侧抽搐| 免费在线观看成人毛片| 少妇丰满av| 制服人妻中文乱码| 欧美一区二区国产精品久久精品| 老司机深夜福利视频在线观看| 成人午夜高清在线视频| 亚洲国产高清在线一区二区三| 亚洲国产色片| 最好的美女福利视频网| 国产久久久一区二区三区| 岛国视频午夜一区免费看| 麻豆一二三区av精品| 免费大片18禁| 婷婷精品国产亚洲av| 美女cb高潮喷水在线观看| ponron亚洲| 国产97色在线日韩免费| 丰满的人妻完整版| 精品久久久久久久人妻蜜臀av| 少妇裸体淫交视频免费看高清| 国语自产精品视频在线第100页| 禁无遮挡网站| 动漫黄色视频在线观看| 性色avwww在线观看| www日本在线高清视频| 19禁男女啪啪无遮挡网站| 日本三级黄在线观看| 久久香蕉精品热| 舔av片在线| 在线播放无遮挡| 黄片小视频在线播放| 国产中年淑女户外野战色| 亚洲午夜理论影院| 国产蜜桃级精品一区二区三区| 别揉我奶头~嗯~啊~动态视频| 香蕉久久夜色| 国产91精品成人一区二区三区| 三级毛片av免费| 国语自产精品视频在线第100页| 长腿黑丝高跟| 老司机午夜福利在线观看视频| 亚洲精品在线美女| 亚洲性夜色夜夜综合| 国产毛片a区久久久久| 亚洲熟妇中文字幕五十中出| 久久久久精品国产欧美久久久| 国产一区二区在线观看日韩 | 国产精品亚洲美女久久久| 听说在线观看完整版免费高清| 看片在线看免费视频| 19禁男女啪啪无遮挡网站| 一级毛片女人18水好多| 一二三四社区在线视频社区8| 夜夜躁狠狠躁天天躁| 欧美黄色片欧美黄色片| 午夜久久久久精精品| 欧美午夜高清在线| 欧美大码av| 在线免费观看的www视频| 欧美大码av| 97超视频在线观看视频| xxx96com| 亚洲熟妇中文字幕五十中出| 国产免费男女视频| 青草久久国产| 国产av在哪里看| 欧洲精品卡2卡3卡4卡5卡区| 国产精品三级大全| 国产欧美日韩一区二区三| 真人做人爱边吃奶动态| 夜夜爽天天搞| 又粗又爽又猛毛片免费看| 精品久久久久久成人av| 亚洲国产欧美人成| 国产乱人伦免费视频| 操出白浆在线播放| 国产主播在线观看一区二区| 内地一区二区视频在线| 久久婷婷人人爽人人干人人爱| 亚洲第一欧美日韩一区二区三区| 神马国产精品三级电影在线观看| 欧美日韩一级在线毛片| 欧美激情久久久久久爽电影| 国产黄a三级三级三级人| 夜夜躁狠狠躁天天躁| 3wmmmm亚洲av在线观看| 国产麻豆成人av免费视频| 熟女电影av网| 欧美zozozo另类| 看黄色毛片网站| 天堂动漫精品| 亚洲国产中文字幕在线视频| 我的老师免费观看完整版| 色精品久久人妻99蜜桃| 波多野结衣巨乳人妻| 久久婷婷人人爽人人干人人爱| 久久精品91无色码中文字幕| 国产伦在线观看视频一区| 久久久久久久久大av| 日韩有码中文字幕| 欧美日韩黄片免| 亚洲国产色片| 亚洲av成人不卡在线观看播放网| 亚洲成人久久爱视频| 欧美又色又爽又黄视频| 首页视频小说图片口味搜索| 美女大奶头视频| 熟女少妇亚洲综合色aaa.| 国产免费av片在线观看野外av| 国产色爽女视频免费观看| 波多野结衣高清作品| 国产伦精品一区二区三区四那| 精品久久久久久久毛片微露脸| 一本综合久久免费| 欧美在线一区亚洲| 欧美日韩综合久久久久久 | 亚洲国产精品合色在线| 99久久无色码亚洲精品果冻| 欧美区成人在线视频| 精品欧美国产一区二区三| 精品电影一区二区在线| 成年女人毛片免费观看观看9| 亚洲国产欧洲综合997久久,| 色老头精品视频在线观看| 女警被强在线播放| 午夜免费观看网址| 又黄又粗又硬又大视频| 一级作爱视频免费观看| 国产蜜桃级精品一区二区三区| 在线观看舔阴道视频| 亚洲成av人片免费观看| 中文字幕人成人乱码亚洲影| 国产伦精品一区二区三区四那| 欧美成人一区二区免费高清观看| 午夜影院日韩av| 怎么达到女性高潮| 啪啪无遮挡十八禁网站| 午夜免费成人在线视频| 国产成人av激情在线播放| 久久久久久久精品吃奶| 欧美又色又爽又黄视频| 久久久久久久亚洲中文字幕 | 亚洲国产精品久久男人天堂| 久久久精品欧美日韩精品| 国产午夜福利久久久久久| 狂野欧美白嫩少妇大欣赏| 在线视频色国产色| svipshipincom国产片| 热99在线观看视频| 午夜久久久久精精品| 99精品欧美一区二区三区四区| 久久人妻av系列| АⅤ资源中文在线天堂| 小蜜桃在线观看免费完整版高清| 深爱激情五月婷婷| www日本在线高清视频| 嫩草影院入口| 日本黄色视频三级网站网址| 一本一本综合久久| 伊人久久大香线蕉亚洲五| 夜夜夜夜夜久久久久| 少妇丰满av| 最近最新中文字幕大全电影3| 国产欧美日韩精品亚洲av| 精品久久久久久久人妻蜜臀av| 黄色视频,在线免费观看| 日本熟妇午夜| 久久天躁狠狠躁夜夜2o2o| 精品一区二区三区视频在线观看免费| 成人av在线播放网站| 最好的美女福利视频网| 成熟少妇高潮喷水视频| 久久精品国产清高在天天线| 成人一区二区视频在线观看| 少妇裸体淫交视频免费看高清|