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

    波浪與水囊潛堤相互作用的數(shù)值模擬*

    2022-03-07 09:14:08楊志堅(jiān)趙西增
    水運(yùn)工程 2022年2期
    關(guān)鍵詞:消波渦量水囊

    王 臣,楊志堅(jiān),趙西增,邱 俞

    (1.舟山市港航工程規(guī)劃設(shè)計(jì)院有限公司,浙江 舟山316000;2.浙江大學(xué) 海洋學(xué)院,浙江 舟山,316021;3.舟山市普陀交通工程有限公司,浙江 舟山316100)

    防波堤是港口、海岸最常見的防護(hù)工程,具有防御波浪入侵、圍護(hù)港池平穩(wěn)和減少海岸侵蝕等功能。傳統(tǒng)防波堤存在建造成本隨水深大幅增加、安裝和維護(hù)費(fèi)用較高、影響近岸自然景觀等問題,為此,新型防波堤的研究和設(shè)計(jì)應(yīng)運(yùn)而生,其中以橡膠袋裝水類的潛堤為代表的新型柔性防波堤具有快速布置、建造成本低、后期維護(hù)方便和不影響水體交換等優(yōu)點(diǎn),引起了廣泛的關(guān)注[1]。

    波浪與水囊結(jié)構(gòu)相互作用是一個(gè)復(fù)雜的、涉及多相耦合的過程,存在結(jié)構(gòu)的非線性變形和波浪破碎和漩渦脫落現(xiàn)象,因此波浪與水囊結(jié)構(gòu)相互作用的理論推導(dǎo)研究不多見,主要集中在物理模型試驗(yàn)和基于勢流理論的數(shù)值模擬方面。Cho等[2]在二維線性水彈性理論中,采用多域邊界元方法得到了數(shù)值解,發(fā)現(xiàn)水下水平柔性薄膜防波堤對后方的掩護(hù)性能通常優(yōu)于相應(yīng)的剛性結(jié)構(gòu)。Lee等[3]基于線性波理論和結(jié)構(gòu)小幅度響應(yīng)的前提下,利用最小二乘法建立了波浪與多排出水垂直薄膜防波堤相互作用的水彈性分析模型,發(fā)現(xiàn)增大膜的張力會降低波浪的透射率。Broderick等[4]將水囊結(jié)構(gòu)視為一個(gè)內(nèi)部壓力為常數(shù)、集中質(zhì)量的系統(tǒng),利用邊界元理論耦合有限元的方法來模擬非線性波與充液膜相互作用問題。Phadke等[5]從圓柱殼理論出發(fā),建立了一個(gè)分區(qū)邊界元與有限元的耦合模型,研究了波浪與水囊潛堤相互作用的水動(dòng)力特性。Phadke等[6-7]關(guān)注了波浪與水囊潛堤相互作用中的線性結(jié)構(gòu)響應(yīng)和共振現(xiàn)象,并考慮了膜結(jié)構(gòu)的幾何非線性變形,發(fā)現(xiàn)即使是相對較小的變形,幾何非線性也很重要,膜幾何非線性增加了系統(tǒng)的剛度,并在某些頻率下引起遲滯響應(yīng)。Das等[8]將分區(qū)邊界元和有限元耦合模型擴(kuò)展到三維,并與二維數(shù)值結(jié)果和Ohyama等[9]的試驗(yàn)結(jié)果相對比,發(fā)現(xiàn)膜系統(tǒng)的共振現(xiàn)象抵消了入射波、降低了波的透射率,合適的膜內(nèi)壓力可以降低入射波頻率范圍內(nèi)的透射波高。Li等[10]采用多極展開法和變量分離法得到了波浪與半圓形水囊潛堤相互作用的理論解,研究水囊潛堤的形狀和初始內(nèi)壓對水囊潛堤消波能力的影響,發(fā)現(xiàn)設(shè)計(jì)合理的水囊潛堤可作為一種有效的防波堤。Liu等[11]基于混合歐拉-拉格朗日方法建立了水囊潛堤與非線性波的相互作用耦合模型。

    綜上可知,相關(guān)研究主要基于勢流理論,由于忽略了流體黏性的影響,對波浪破碎和漩渦脫落現(xiàn)象都未涉及,同時(shí)無法精確模擬強(qiáng)非線性波浪運(yùn)動(dòng)和結(jié)構(gòu)的大幅運(yùn)動(dòng)響應(yīng),水囊潛堤消波特性的機(jī)理尚未清晰。本文將基于高階有限差分法求解黏性流方程,采用有限元方法求解結(jié)構(gòu)動(dòng)力學(xué)方程,采用浸入邊界法實(shí)現(xiàn)流體和固體求解器的耦合求解,構(gòu)建黏性流數(shù)值波浪水槽,開展波浪與水囊潛堤相互作用的數(shù)值模擬研究。

    1 數(shù)值模型

    1.1 流場求解器

    考慮流體為二維不可壓縮的黏性流體,其控制方程由連續(xù)性方程和動(dòng)量方程組成,具體形式如下:

    ?·u=0

    (1)

    (2)

    式中:u為速度矢量;p為壓強(qiáng);ρ為流體密度;μ為動(dòng)力黏性系數(shù);F為質(zhì)量力。

    本模型求解控制方程采用分步算法:1)采用緊致插值曲線求解對流項(xiàng),得到中間速度u*;2)采用中心差分法求解非對流項(xiàng)Ⅰ得到第2個(gè)中間速度u**;3)最后利用SOR(Successive Over-Relation)算法進(jìn)行壓力-速度耦合,求解非對流項(xiàng)Ⅱ,得到速度解u和壓強(qiáng)解p,具體求解流程可參考文獻(xiàn)[12-13]。

    1.2 結(jié)構(gòu)求解器

    結(jié)構(gòu)運(yùn)動(dòng)是以位移為基本未知量的有限元法,待求解對象經(jīng)過有限元離散、單元分析、系統(tǒng)組集和引入邊界條件后,可得到結(jié)構(gòu)運(yùn)動(dòng)控制方程[14],表達(dá)式為:

    (3)

    1.3 流固耦合

    流固耦合采用改進(jìn)的虛擬單元浸入邊界法,通過高階插值賦予虛擬單元流體屬性,進(jìn)而在流固邊界設(shè)置不可穿透與無滑移條件,實(shí)現(xiàn)流體和固體高精度耦合求解。改進(jìn)的虛擬單元法插值方案見圖1,該方法詳細(xì)的推導(dǎo)和驗(yàn)證見文獻(xiàn)[16-17]。

    圖1 改進(jìn)的虛擬單元法插值方案

    圖2 流固耦合求解流程

    2 數(shù)值設(shè)置及驗(yàn)證

    2.1 數(shù)值設(shè)置

    波浪與水囊潛堤相互作用的試驗(yàn)在浙江大學(xué)紫金港校區(qū)的斷面波浪數(shù)值水槽中進(jìn)行,參照試驗(yàn)布置,數(shù)值波浪水槽設(shè)置見圖3,其水囊潛堤結(jié)構(gòu)模型位置、波面監(jiān)測點(diǎn)的設(shè)置均與試驗(yàn)設(shè)置相似。數(shù)值水槽總長34 m,采用源造波法產(chǎn)生波浪,源點(diǎn)設(shè)置在x=0 m處。在水槽左側(cè)x=-8 m至x=0 m處和右側(cè)x=18 m至x=26 m處為海綿層消波區(qū),可防止邊界處的波浪反射;x=0 m處為源項(xiàng)造波點(diǎn),可產(chǎn)生向左右兩個(gè)方向傳播的規(guī)則波,水囊潛堤模型中心位于x=9 m處。水囊潛堤使用氯丁橡膠薄膜,厚度為1 mm,彈性模量為3 MPa,密度為1.24 g/cm3,泊松比為0.47。試驗(yàn)中采用高速相機(jī)和浪高儀,分別布置了4個(gè)浪高儀,具體位置如圖3中G1~G4所示,分別觀測水囊結(jié)構(gòu)的變形及波浪經(jīng)過潛堤的變形問題。

    圖3 水囊潛堤的數(shù)值波浪水槽(單位:m)

    2.2 數(shù)值驗(yàn)證

    圖4表示水深0.4 m、波浪頻率f=0.7 Hz時(shí)預(yù)測得到的結(jié)構(gòu)響應(yīng)及波面與試驗(yàn)結(jié)果的定性比較驗(yàn)證,分別考慮了剛體、完全充盈和未完全充盈3種狀態(tài)。可發(fā)現(xiàn)剛性水囊潛堤的數(shù)值計(jì)算結(jié)果在波面形態(tài)與試驗(yàn)吻合。充盈的柔性水囊潛堤的數(shù)值計(jì)算結(jié)果在波面形態(tài)和結(jié)構(gòu)響應(yīng)上,幾乎完全與試驗(yàn)吻合。充盈度85%的柔性水囊潛堤數(shù)值計(jì)算結(jié)果在結(jié)構(gòu)響應(yīng)形狀和幅度整體上與試驗(yàn)結(jié)果相比基本一致。

    圖4 水深0.4 m波浪頻率f=0.7 Hz試驗(yàn)與數(shù)模對比的結(jié)構(gòu)響應(yīng)

    圖5給出了試驗(yàn)中的波浪透射系數(shù)與數(shù)值計(jì)算結(jié)果的對比。對于波浪與不同類型的水囊潛堤相互作用后的透射系數(shù),數(shù)模計(jì)算結(jié)果和試驗(yàn)結(jié)果略有差異,特別是未充盈的柔性水囊潛堤,但趨勢上吻合良好。這是由于波浪與未充盈的柔性水囊潛堤相互作用產(chǎn)生的強(qiáng)非線性耦合作用所致。通過對比水深0.4 m和水深0.335 m的透浪系數(shù),發(fā)現(xiàn)水深0.335 m的數(shù)值對比結(jié)果的差異更大。這主要是由于水深0.335 m時(shí)產(chǎn)生了波浪破碎現(xiàn)象,且未充盈的柔性水囊潛堤的復(fù)雜的結(jié)構(gòu)運(yùn)動(dòng)響應(yīng),強(qiáng)非線性的耦合作用進(jìn)一步增強(qiáng)。

    圖5 試驗(yàn)與數(shù)模對比的透射系數(shù)CT隨無因次水深kh的變化

    3 結(jié)果分析

    3.1 充盈度變化對水囊潛堤消波性能的影響

    圖6為水深0.4 m、波浪頻率0.7 Hz不同充盈度下柔性水囊潛堤的結(jié)構(gòu)變形。由圖6可知,隨著充盈度的減小,充盈的柔性水囊潛堤轉(zhuǎn)變?yōu)槲闯溆娜嵝运覞摰?,柔性水囊潛堤的結(jié)構(gòu)響應(yīng)增大,充盈度85%、80%出現(xiàn)了大幅度的結(jié)構(gòu)響應(yīng);當(dāng)充盈度降低到75%時(shí),未充盈的柔性水囊潛堤結(jié)構(gòu)響應(yīng)減小,不再出現(xiàn)大幅度的擺動(dòng),只在頂部左右搖擺;當(dāng)充盈度降低到70%時(shí),未充盈的柔性水囊潛堤塌陷成充盈的柔性水囊潛堤,柔性水囊潛堤的結(jié)構(gòu)響應(yīng)與充盈度100%的水囊相似,此時(shí)冠層高度(靜水面到結(jié)構(gòu)物頂點(diǎn)的豎向距離)增大,這表明在一定范圍內(nèi)減小充盈度可增大柔性水囊潛堤的結(jié)構(gòu)響應(yīng),未充盈的柔性水囊潛堤的結(jié)構(gòu)響應(yīng)大于充盈的柔性水囊潛堤,但充盈度減小過多會導(dǎo)致未充盈的柔性水囊潛堤塌陷成充盈的柔性水囊潛堤,反而降低了水囊潛堤的結(jié)構(gòu)高度。

    圖6 水深0.4 m不同充盈度下柔性水囊潛堤的結(jié)構(gòu)響應(yīng)

    圖7為不同充盈度下水囊潛堤的波浪透射系數(shù),可發(fā)現(xiàn)充盈的柔性水囊潛堤的透射系數(shù)最大,消波性能最差;充盈度85%、80%的柔性水囊潛堤透射系數(shù)相接近且最小,消波性能最好;充盈度75%的柔性水囊潛堤的透射系數(shù)高于充盈度80%的柔性水囊潛堤,消波性能變差,這是由于充盈度75%的柔性水囊潛堤結(jié)構(gòu)響應(yīng)減小,不再出現(xiàn)大幅度的擺動(dòng),只在頂部左右搖擺,而水囊潛堤的結(jié)構(gòu)響應(yīng)與其消波性能相關(guān),因此其透射性能降低。充盈度70%的柔性水囊潛堤的透射系數(shù)與充盈度100%的水囊潛堤相接近,這是由于充盈度70%的柔性水囊潛堤塌陷成充盈的柔性水囊潛堤的緣故,冠層高度增加,進(jìn)一步降低了充盈的柔性水囊潛堤的消波性能。

    3.2 水囊潛堤消波機(jī)理

    圖8為波浪經(jīng)過不同整體剛度水囊潛堤過程中的渦量場。由圖8可知,剛性水囊潛堤只在結(jié)構(gòu)物附近產(chǎn)生少量渦量,看不到尾渦的演化,這些渦量是由于結(jié)構(gòu)物的阻礙作用和淺水效應(yīng)所致。相較于剛性水囊潛堤,對于充盈的(充盈度100%和充盈度70%)柔性水囊潛堤而言,由于受迫微幅振動(dòng)在結(jié)構(gòu)物附近產(chǎn)生的渦量增加,在結(jié)構(gòu)物的背浪側(cè)能明顯看到一個(gè)完全發(fā)育的正渦,但是漩渦基本不脫落,充盈度70%的柔性水囊潛堤由于冠層高度增加,漩渦演化位置相較于充盈度100%的柔性水囊潛堤出現(xiàn)上移現(xiàn)象。對于未充盈的柔性水囊潛堤而言,相較于充盈度100%的柔性水囊潛堤,充盈度90%的柔性水囊潛堤的尾渦發(fā)生了脫落,表明其結(jié)構(gòu)響應(yīng)增強(qiáng)促進(jìn)了漩渦的脫落;充盈度85%的柔性水囊潛堤出現(xiàn)了大幅度的結(jié)構(gòu)響應(yīng),產(chǎn)生的渦量最多,在膜結(jié)構(gòu)頂部能夠觀察到明顯正渦和負(fù)渦交替脫落的現(xiàn)象,并在結(jié)構(gòu)物背浪側(cè)能明顯看到一個(gè)發(fā)育完全的正渦發(fā)生脫落;充盈度75%的柔性水囊潛堤只在頂部左右擺動(dòng),沒有出現(xiàn)大幅度的運(yùn)動(dòng),其頂部會出現(xiàn)交替脫落的漩渦,漩渦強(qiáng)度較小,這表明渦量的演化與結(jié)構(gòu)物的結(jié)構(gòu)響應(yīng)程度相關(guān)。此外,由于柔性水囊潛堤受迫振動(dòng)擾動(dòng)水體,在膜結(jié)構(gòu)附近處不斷產(chǎn)生渦量,并傳播到水面附近與入射波碰撞,消散一部分波能。

    圖8 水深0.4 m渦量場的數(shù)值計(jì)算結(jié)果

    4 結(jié)論

    1)基于有限差分法與有限元法的耦合模型可有效模擬水囊潛堤的消波過程。

    2)合適的充盈度范圍內(nèi),未充盈的柔性水囊潛堤存在大幅度振動(dòng)促進(jìn)了尾渦的演化與脫落,并且柔性水囊潛堤受迫振動(dòng)擾動(dòng)水體,在膜結(jié)構(gòu)附近處不斷產(chǎn)生渦量,消散了部分波能,具有更好的消波效果。

    3)未完全充盈的潛堤可增大柔性水囊潛堤的結(jié)構(gòu)響應(yīng),未充盈的柔性水囊潛堤的結(jié)構(gòu)響應(yīng)大于完全充盈的柔性水囊潛堤,但充盈度減小過多會導(dǎo)致未充盈的柔性水囊潛堤塌陷、降低消波效果。

    猜你喜歡
    消波渦量水囊
    主體結(jié)構(gòu)荷載可控的新型組合式防護(hù)結(jié)構(gòu)(Ⅱ):影響因素及設(shè)計(jì)理念
    爆炸與沖擊(2025年1期)2025-02-19 00:00:00
    水囊聯(lián)合小劑量催產(chǎn)素應(yīng)用于足月引產(chǎn)的臨床療效觀察
    含沙空化對軸流泵內(nèi)渦量分布的影響
    循環(huán)水槽多層孔板消波裝置開發(fā)及消波特性數(shù)值模擬
    自由表面渦流動(dòng)現(xiàn)象的數(shù)值模擬
    足月妊娠不同規(guī)格水囊引產(chǎn)促宮頸成熟的療效分析
    紗布填塞與水囊壓迫治療產(chǎn)后大出血患者的臨床效果觀察
    多層變孔徑傾斜孔板式消波裝置的試驗(yàn)研究
    航態(tài)對大型船舶甲板氣流場的影響
    T型透空式防波堤消波性能的理論分析
    中出人妻视频一区二区| 亚洲精品国产一区二区精华液| 国产欧美日韩综合在线一区二区| 操出白浆在线播放| 中文字幕人妻熟女乱码| 日韩高清综合在线| 香蕉国产在线看| 午夜精品久久久久久毛片777| 成年人黄色毛片网站| 18禁黄网站禁片午夜丰满| e午夜精品久久久久久久| 国产色视频综合| 麻豆国产av国片精品| 亚洲一卡2卡3卡4卡5卡精品中文| 十八禁人妻一区二区| 国产熟女午夜一区二区三区| 日本一区二区免费在线视频| 亚洲人成电影免费在线| 女人精品久久久久毛片| 亚洲男人的天堂狠狠| 亚洲国产精品久久男人天堂| 人人澡人人妻人| 真人做人爱边吃奶动态| 中文字幕高清在线视频| 亚洲av片天天在线观看| 欧美性长视频在线观看| 国产亚洲精品一区二区www| 757午夜福利合集在线观看| 9191精品国产免费久久| 欧美日韩瑟瑟在线播放| 夜夜看夜夜爽夜夜摸| 成人av一区二区三区在线看| av天堂在线播放| 精品一品国产午夜福利视频| 亚洲av美国av| 99热只有精品国产| 国产精品免费视频内射| 精品一区二区三区视频在线观看免费| 黑人巨大精品欧美一区二区蜜桃| 村上凉子中文字幕在线| 一区二区日韩欧美中文字幕| 中文字幕人妻丝袜一区二区| 免费高清视频大片| 国产精品免费一区二区三区在线| 久久香蕉国产精品| 欧美激情高清一区二区三区| 欧美在线一区亚洲| 亚洲色图 男人天堂 中文字幕| 精品国产一区二区久久| 性少妇av在线| 最新在线观看一区二区三区| 精品日产1卡2卡| 最新美女视频免费是黄的| 精品卡一卡二卡四卡免费| www日本在线高清视频| 国产精品99久久99久久久不卡| 丝袜在线中文字幕| 一级a爱片免费观看的视频| 精品国产一区二区三区四区第35| 国产精品,欧美在线| 久久精品人人爽人人爽视色| 嫩草影院精品99| 亚洲av第一区精品v没综合| 日韩高清综合在线| 国产麻豆69| 国产激情久久老熟女| 亚洲天堂国产精品一区在线| 波多野结衣巨乳人妻| 亚洲午夜理论影院| 亚洲最大成人中文| 国产一级毛片七仙女欲春2 | 精品国产超薄肉色丝袜足j| 国产av在哪里看| 国产私拍福利视频在线观看| 久久精品影院6| 久久人妻熟女aⅴ| 久久九九热精品免费| 国产97色在线日韩免费| 美女扒开内裤让男人捅视频| 人妻丰满熟妇av一区二区三区| 日韩国内少妇激情av| 麻豆成人av在线观看| 99国产极品粉嫩在线观看| 国产免费av片在线观看野外av| 久久久国产精品麻豆| av视频在线观看入口| 欧美人与性动交α欧美精品济南到| 操美女的视频在线观看| 久久人人爽av亚洲精品天堂| 麻豆一二三区av精品| av电影中文网址| 国产亚洲av高清不卡| av中文乱码字幕在线| 91麻豆精品激情在线观看国产| 视频区欧美日本亚洲| 午夜亚洲福利在线播放| 国产精华一区二区三区| 精品久久蜜臀av无| 国内久久婷婷六月综合欲色啪| 精品国产国语对白av| 成人特级黄色片久久久久久久| 国产精品久久久久久亚洲av鲁大| 97人妻精品一区二区三区麻豆 | 久久人人97超碰香蕉20202| 日本免费一区二区三区高清不卡 | 亚洲成国产人片在线观看| www.熟女人妻精品国产| 真人做人爱边吃奶动态| 黑人巨大精品欧美一区二区mp4| 亚洲专区中文字幕在线| 午夜影院日韩av| 国产精品二区激情视频| 黄色视频,在线免费观看| 国产高清有码在线观看视频 | 国产伦一二天堂av在线观看| 动漫黄色视频在线观看| 午夜精品国产一区二区电影| 日韩欧美国产一区二区入口| 精品国产美女av久久久久小说| 欧美日韩亚洲综合一区二区三区_| 亚洲国产欧美一区二区综合| 精品国产一区二区久久| 伦理电影免费视频| 波多野结衣巨乳人妻| 精品少妇一区二区三区视频日本电影| 欧美激情 高清一区二区三区| 久久狼人影院| 19禁男女啪啪无遮挡网站| 满18在线观看网站| 可以在线观看的亚洲视频| 亚洲成a人片在线一区二区| 啦啦啦观看免费观看视频高清 | 成在线人永久免费视频| 夜夜爽天天搞| 女人精品久久久久毛片| 两个人看的免费小视频| 看黄色毛片网站| 久久精品91蜜桃| 在线观看一区二区三区| 亚洲av电影在线进入| 91国产中文字幕| 亚洲精品国产一区二区精华液| 91九色精品人成在线观看| 久久国产精品影院| 久久久国产欧美日韩av| 亚洲久久久国产精品| 午夜福利在线观看吧| 国产精品美女特级片免费视频播放器 | 国产精品久久久久久亚洲av鲁大| 亚洲国产精品999在线| 99在线人妻在线中文字幕| 18禁裸乳无遮挡免费网站照片 | 精品欧美国产一区二区三| 日本在线视频免费播放| 成人永久免费在线观看视频| 午夜福利高清视频| 色综合亚洲欧美另类图片| 日韩视频一区二区在线观看| 久久性视频一级片| 亚洲第一电影网av| 国产精品乱码一区二三区的特点 | 国产精华一区二区三区| 欧美另类亚洲清纯唯美| 亚洲一区二区三区不卡视频| 在线观看66精品国产| 亚洲国产毛片av蜜桃av| 久久久国产精品麻豆| 99香蕉大伊视频| 欧美国产日韩亚洲一区| 在线十欧美十亚洲十日本专区| 精品人妻1区二区| 国产午夜精品久久久久久| а√天堂www在线а√下载| 日本一区二区免费在线视频| 我的亚洲天堂| 手机成人av网站| 欧美亚洲日本最大视频资源| 亚洲国产毛片av蜜桃av| 久久久国产精品麻豆| 这个男人来自地球电影免费观看| 91麻豆av在线| 国产精品亚洲美女久久久| 乱人伦中国视频| 免费av毛片视频| 免费人成视频x8x8入口观看| 人人妻人人澡人人看| 制服人妻中文乱码| 亚洲最大成人中文| 午夜影院日韩av| 嫩草影院精品99| 非洲黑人性xxxx精品又粗又长| 多毛熟女@视频| 搞女人的毛片| 国产野战对白在线观看| 久久久久久大精品| 国产成人影院久久av| 女人爽到高潮嗷嗷叫在线视频| 在线观看舔阴道视频| 国产视频一区二区在线看| 日韩精品免费视频一区二区三区| 精品国产乱码久久久久久男人| 久久这里只有精品19| 国产高清videossex| 色婷婷久久久亚洲欧美| 亚洲成人精品中文字幕电影| 91成年电影在线观看| 精品欧美国产一区二区三| 黑人欧美特级aaaaaa片| 老司机在亚洲福利影院| 国产又爽黄色视频| 国产成人精品在线电影| 三级毛片av免费| 国产99久久九九免费精品| 一边摸一边做爽爽视频免费| 亚洲成人精品中文字幕电影| 国产精品自产拍在线观看55亚洲| 国产麻豆成人av免费视频| 中文字幕久久专区| 不卡一级毛片| 日韩欧美一区视频在线观看| 两个人视频免费观看高清| 免费不卡黄色视频| 欧美一级a爱片免费观看看 | 一区二区三区精品91| 国产精品永久免费网站| 97人妻天天添夜夜摸| 两个人免费观看高清视频| 性色av乱码一区二区三区2| 欧美黄色淫秽网站| 美国免费a级毛片| 亚洲自拍偷在线| 天堂√8在线中文| 久久久久九九精品影院| 亚洲第一电影网av| 777久久人妻少妇嫩草av网站| av天堂在线播放| av在线天堂中文字幕| 99久久国产精品久久久| 99精品在免费线老司机午夜| 欧美不卡视频在线免费观看 | 久久国产亚洲av麻豆专区| 又紧又爽又黄一区二区| 少妇被粗大的猛进出69影院| 久久久久久久久久久久大奶| 午夜久久久在线观看| 99久久99久久久精品蜜桃| 国产亚洲精品第一综合不卡| av在线天堂中文字幕| 首页视频小说图片口味搜索| 精品第一国产精品| 9热在线视频观看99| 一a级毛片在线观看| 免费搜索国产男女视频| 亚洲成国产人片在线观看| 欧美性长视频在线观看| 精品国产一区二区三区四区第35| 国产精品爽爽va在线观看网站 | 叶爱在线成人免费视频播放| 精品欧美国产一区二区三| 啦啦啦观看免费观看视频高清 | 国产熟女午夜一区二区三区| 国产成人av激情在线播放| 妹子高潮喷水视频| 亚洲熟妇熟女久久| 国产成人av教育| 精品久久久久久久久久免费视频| 亚洲精品在线美女| 91精品三级在线观看| 亚洲精华国产精华精| 午夜久久久在线观看| 一级a爱视频在线免费观看| 亚洲少妇的诱惑av| 极品教师在线免费播放| 免费不卡黄色视频| 亚洲专区字幕在线| 97碰自拍视频| 大陆偷拍与自拍| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久久中文| 亚洲少妇的诱惑av| 精品国产乱码久久久久久男人| 精品无人区乱码1区二区| 色播在线永久视频| 午夜影院日韩av| 51午夜福利影视在线观看| 少妇的丰满在线观看| 97人妻天天添夜夜摸| 18禁国产床啪视频网站| 欧美人与性动交α欧美精品济南到| 人人澡人人妻人| 国产精品综合久久久久久久免费 | 国产97色在线日韩免费| 久久精品亚洲熟妇少妇任你| 国产野战对白在线观看| cao死你这个sao货| 国产aⅴ精品一区二区三区波| 曰老女人黄片| 久热这里只有精品99| 黄片大片在线免费观看| 怎么达到女性高潮| 午夜精品国产一区二区电影| 久久人妻福利社区极品人妻图片| 美女扒开内裤让男人捅视频| 欧美国产精品va在线观看不卡| 亚洲少妇的诱惑av| 欧美日韩乱码在线| 老司机福利观看| 国产精品亚洲一级av第二区| 精品国产美女av久久久久小说| 午夜福利视频1000在线观看 | 亚洲精品国产精品久久久不卡| 极品人妻少妇av视频| 国产亚洲精品综合一区在线观看 | 国产精品电影一区二区三区| 动漫黄色视频在线观看| АⅤ资源中文在线天堂| 日本一区二区免费在线视频| 久久精品亚洲熟妇少妇任你| 91大片在线观看| 久久人人97超碰香蕉20202| 一夜夜www| 国产精品一区二区在线不卡| 久久婷婷成人综合色麻豆| 国产午夜精品久久久久久| 成人国语在线视频| 国产精品免费视频内射| 99精品欧美一区二区三区四区| 国产区一区二久久| 亚洲一区中文字幕在线| 亚洲男人天堂网一区| 精品国产一区二区久久| 欧美不卡视频在线免费观看 | 757午夜福利合集在线观看| 欧美老熟妇乱子伦牲交| 国产成人av教育| 免费无遮挡裸体视频| 男女做爰动态图高潮gif福利片 | 久久国产亚洲av麻豆专区| 一a级毛片在线观看| 老鸭窝网址在线观看| 欧美最黄视频在线播放免费| 亚洲成人精品中文字幕电影| 1024视频免费在线观看| 精品久久久久久久人妻蜜臀av | 叶爱在线成人免费视频播放| 级片在线观看| www.自偷自拍.com| 色精品久久人妻99蜜桃| 日日摸夜夜添夜夜添小说| 久久人人97超碰香蕉20202| 亚洲七黄色美女视频| 免费久久久久久久精品成人欧美视频| 黄网站色视频无遮挡免费观看| 99国产精品免费福利视频| 免费在线观看亚洲国产| 777久久人妻少妇嫩草av网站| 天堂√8在线中文| 久久久精品国产亚洲av高清涩受| 咕卡用的链子| 国产野战对白在线观看| 久久精品成人免费网站| 久久国产亚洲av麻豆专区| 亚洲国产看品久久| 大香蕉久久成人网| 欧美日本视频| 国产精品免费一区二区三区在线| 欧美色欧美亚洲另类二区 | 亚洲五月天丁香| 九色国产91popny在线| 丝袜在线中文字幕| 亚洲国产精品成人综合色| 国产成人免费无遮挡视频| 午夜a级毛片| 无人区码免费观看不卡| 精品福利观看| 亚洲av日韩精品久久久久久密| av片东京热男人的天堂| 亚洲三区欧美一区| 国产伦人伦偷精品视频| 精品一区二区三区av网在线观看| 欧美av亚洲av综合av国产av| 国产精品久久久久久人妻精品电影| 啦啦啦韩国在线观看视频| 国产在线观看jvid| 一区二区三区精品91| 国产伦一二天堂av在线观看| 国产主播在线观看一区二区| 欧美色视频一区免费| 中文字幕精品免费在线观看视频| 亚洲久久久国产精品| 91麻豆av在线| 999久久久精品免费观看国产| 亚洲狠狠婷婷综合久久图片| 午夜成年电影在线免费观看| 91字幕亚洲| 亚洲黑人精品在线| 丝袜美足系列| 精品国产美女av久久久久小说| 亚洲欧美激情综合另类| 久久香蕉精品热| 成人亚洲精品一区在线观看| 久久伊人香网站| 日韩精品青青久久久久久| 久久天躁狠狠躁夜夜2o2o| 精品卡一卡二卡四卡免费| 亚洲精品久久成人aⅴ小说| 很黄的视频免费| 精品国产超薄肉色丝袜足j| 男女下面插进去视频免费观看| 91精品国产国语对白视频| avwww免费| 国产极品粉嫩免费观看在线| 亚洲专区国产一区二区| 91av网站免费观看| 成人欧美大片| 一区二区日韩欧美中文字幕| 精品熟女少妇八av免费久了| 激情在线观看视频在线高清| 亚洲狠狠婷婷综合久久图片| 正在播放国产对白刺激| 日韩免费av在线播放| 视频在线观看一区二区三区| 搡老熟女国产l中国老女人| 在线免费观看的www视频| 亚洲成a人片在线一区二区| 动漫黄色视频在线观看| 成年版毛片免费区| 男女下面进入的视频免费午夜 | 久久中文字幕人妻熟女| 欧美最黄视频在线播放免费| 精品人妻在线不人妻| 18禁黄网站禁片午夜丰满| 亚洲中文字幕日韩| 亚洲黑人精品在线| 韩国av一区二区三区四区| 国产一区二区三区综合在线观看| 老司机午夜福利在线观看视频| 精品电影一区二区在线| 国产色视频综合| 一级黄色大片毛片| 18禁美女被吸乳视频| 精品卡一卡二卡四卡免费| 精品乱码久久久久久99久播| 欧美一级毛片孕妇| 亚洲精品在线观看二区| 色精品久久人妻99蜜桃| 欧美大码av| 美女大奶头视频| 啪啪无遮挡十八禁网站| 亚洲男人的天堂狠狠| 久久国产乱子伦精品免费另类| 69精品国产乱码久久久| 成人永久免费在线观看视频| АⅤ资源中文在线天堂| 午夜福利成人在线免费观看| 最好的美女福利视频网| 国产亚洲精品av在线| 日韩欧美一区视频在线观看| 操美女的视频在线观看| 国产午夜福利久久久久久| 这个男人来自地球电影免费观看| 美国免费a级毛片| 亚洲熟妇熟女久久| 亚洲男人的天堂狠狠| 国产av在哪里看| 女警被强在线播放| 亚洲人成网站在线播放欧美日韩| av视频在线观看入口| 香蕉久久夜色| 国产精品日韩av在线免费观看 | 国产成年人精品一区二区| 国产在线精品亚洲第一网站| 丁香欧美五月| 99久久综合精品五月天人人| 国产在线观看jvid| 在线av久久热| 亚洲精品在线观看二区| 两性夫妻黄色片| 一边摸一边做爽爽视频免费| 啪啪无遮挡十八禁网站| 欧美不卡视频在线免费观看 | 黄色 视频免费看| 视频在线观看一区二区三区| 午夜激情av网站| av视频在线观看入口| 亚洲午夜理论影院| 久久久久久大精品| 在线国产一区二区在线| 亚洲久久久国产精品| 九色国产91popny在线| 国产私拍福利视频在线观看| 伦理电影免费视频| 午夜福利一区二区在线看| 自拍欧美九色日韩亚洲蝌蚪91| 欧美av亚洲av综合av国产av| 国产色视频综合| 亚洲中文字幕日韩| 久久久久国内视频| 一级黄色大片毛片| 一边摸一边抽搐一进一小说| 国产精品1区2区在线观看.| 日本a在线网址| 9191精品国产免费久久| 国产一区二区在线av高清观看| 黄色a级毛片大全视频| 国产精品免费一区二区三区在线| 久久影院123| 久久精品国产亚洲av香蕉五月| 国产亚洲欧美98| 波多野结衣高清无吗| 国产极品粉嫩免费观看在线| 18禁观看日本| 无人区码免费观看不卡| 亚洲av成人不卡在线观看播放网| 两性夫妻黄色片| 亚洲精品国产区一区二| 国产熟女午夜一区二区三区| 99久久久亚洲精品蜜臀av| 免费在线观看黄色视频的| 久久久久久人人人人人| 国产亚洲av嫩草精品影院| 欧美成人午夜精品| 十八禁人妻一区二区| 亚洲五月色婷婷综合| 午夜免费激情av| 国产一级毛片七仙女欲春2 | 男人操女人黄网站| 亚洲精品在线美女| 一边摸一边抽搐一进一小说| 麻豆成人av在线观看| 久久精品91无色码中文字幕| 色老头精品视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 久久精品亚洲精品国产色婷小说| 90打野战视频偷拍视频| 亚洲 欧美一区二区三区| 欧美日本中文国产一区发布| 欧美一级a爱片免费观看看 | 99国产精品一区二区三区| 国产伦人伦偷精品视频| 曰老女人黄片| 这个男人来自地球电影免费观看| 国产亚洲av高清不卡| 制服诱惑二区| 精品一区二区三区视频在线观看免费| 精品欧美国产一区二区三| 精品电影一区二区在线| 国产麻豆69| 村上凉子中文字幕在线| www国产在线视频色| 男女午夜视频在线观看| 精品国产一区二区久久| 久久人人97超碰香蕉20202| 午夜老司机福利片| 18禁观看日本| 亚洲人成伊人成综合网2020| 999精品在线视频| 香蕉国产在线看| 99久久久亚洲精品蜜臀av| 久99久视频精品免费| 韩国精品一区二区三区| 午夜老司机福利片| av中文乱码字幕在线| 午夜老司机福利片| 欧美一区二区精品小视频在线| 午夜福利视频1000在线观看 | 亚洲五月婷婷丁香| 97人妻天天添夜夜摸| a级毛片在线看网站| 亚洲五月色婷婷综合| 91国产中文字幕| 亚洲国产欧美日韩在线播放| 黄色片一级片一级黄色片| 夜夜看夜夜爽夜夜摸| 午夜精品在线福利| 999久久久国产精品视频| 亚洲成a人片在线一区二区| 午夜影院日韩av| 精品欧美一区二区三区在线| 男女下面进入的视频免费午夜 | 后天国语完整版免费观看| 国产一区二区三区视频了| tocl精华| 伊人久久大香线蕉亚洲五| 电影成人av| 在线观看66精品国产| 桃红色精品国产亚洲av| 身体一侧抽搐| 午夜老司机福利片| 亚洲少妇的诱惑av| 自拍欧美九色日韩亚洲蝌蚪91| 在线永久观看黄色视频| 成人18禁在线播放| 岛国视频午夜一区免费看| 日日夜夜操网爽| 黄色视频,在线免费观看| 在线十欧美十亚洲十日本专区| 中国美女看黄片| 纯流量卡能插随身wifi吗| 日本撒尿小便嘘嘘汇集6| 国产男靠女视频免费网站| 成人亚洲精品av一区二区| 岛国视频午夜一区免费看| 99国产精品99久久久久| 欧美成人性av电影在线观看| 夜夜看夜夜爽夜夜摸| 后天国语完整版免费观看| 成人国产一区最新在线观看| 狂野欧美激情性xxxx| 国产一卡二卡三卡精品| 亚洲 欧美一区二区三区| 搡老妇女老女人老熟妇| 国语自产精品视频在线第100页| 午夜福利视频1000在线观看 | 精品无人区乱码1区二区|