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

    基于改進(jìn)MPS方法的彈性支撐側(cè)壁對潰壩現(xiàn)象的影響分析

    2019-05-17 02:30:12吳巧瑞謝永和白興蘭
    關(guān)鍵詞:水柱潰壩側(cè)壁

    吳巧瑞,王 磊,謝永和,白興蘭

    (浙江海洋大學(xué)船舶與機(jī)電工程學(xué)院,浙江舟山 316022)

    潰壩現(xiàn)象的物理模型是使用擋板在水箱內(nèi)維持一定高度的水柱,在初始時(shí)刻迅速撤去擋板,在重力的作用下,水柱會(huì)坍塌形成一個(gè)自由表面波,并對周圍的障礙物進(jìn)行砰擊作用的一種現(xiàn)象,見圖1。此現(xiàn)象廣泛存在于各個(gè)領(lǐng)域,尤其在船舶與海洋結(jié)構(gòu)物設(shè)計(jì)制造、海洋水動(dòng)力學(xué)領(lǐng)域具有廣泛的應(yīng)用背景,經(jīng)常用作甲板上浪、液艙晃蕩、入水沖擊、海洋平臺(tái)遭受波浪載荷等工程實(shí)際問題的簡化來研究,其共同點(diǎn)是在砰擊過程中具有流體的大變形和強(qiáng)間斷性。因此,盡管潰壩現(xiàn)象的物理模型非常簡單,卻一直吸引著大量學(xué)者對其進(jìn)行研究,包括實(shí)驗(yàn)研究[1-3]和數(shù)值研究[4-6],但潰壩問題仍然是大變形自由液面問題處理中的難點(diǎn)之一。

    鑒于潰壩現(xiàn)象的特殊性,無網(wǎng)格方法越來越被廣泛地應(yīng)用來處理這類具有大變形自由液面的流動(dòng)問題,MPS方法就是其中之一。MPS方法最早是由日本東京大學(xué)的KOSHIZUKA S,et al[7]于1995年提出的,相對于其它無網(wǎng)格方法,MPS方法采用預(yù)估-修正的半隱式方法來求解不可壓縮流體的控制方程,從而可以有效地避免離散對流項(xiàng)所引起的數(shù)值發(fā)散問題,另外,粒子具有運(yùn)動(dòng)性并具有非常便捷的自由液面判別方法,這些都使得MPS方法在處理自由液面問題尤其是大變形問題中有其獨(dú)特的優(yōu)勢。但同時(shí),由于該方法提出不久,MPS方法本身在穩(wěn)定性、精確性方面還存在一定的問題。因此,筆者對原始MPS方法進(jìn)行了改進(jìn)[8],改進(jìn)后的MPS方法采用弧度判定法和粒子碰撞模型,從而可以有效地避免粒子聚集現(xiàn)象和粒子的誤判現(xiàn)象,大大提高了數(shù)值模擬的穩(wěn)定性和準(zhǔn)確性。本文采用改進(jìn)的MPS方法,通過模擬固定側(cè)壁的潰壩現(xiàn)象驗(yàn)證了該方法在潰壩問題中的有效性,探究了彈性支撐側(cè)壁對潰壩現(xiàn)象中的流動(dòng)形態(tài)、側(cè)壁受力以及表面波對側(cè)壁的砰擊壓力的影響規(guī)律,旨在為解決實(shí)際工程問題提供一定的借鑒。

    圖1 潰壩物理模型圖Fig.1 Physical model of the dam break

    1 改進(jìn)的MPS方法數(shù)值模型

    對于不可壓縮粘性流體,其控制方程為連續(xù)方程和N-S方程,如下所示:

    其中,ρ為流體密度;v為流體質(zhì)點(diǎn)的速度;f為外部的作用力;ρ為流體質(zhì)點(diǎn)的壓力;υ為流體的運(yùn)動(dòng)學(xué)粘性系數(shù)。

    MPS方法是通過預(yù)估-修正求解不可壓縮流體的控制方程,具體實(shí)現(xiàn)方法如下。

    1.1 粒子相互作用模型

    在MPS方法中,采用可以自由移動(dòng)的粒子來離散問題域,相鄰之間的相互作用是通過一個(gè)核函數(shù)來衡量,本文采用的是Koshizuka和Oka[9]在1996年提出的核函數(shù),如下:

    在該式中,r是相鄰兩個(gè)粒子之間的距離,re是核函數(shù)的控制半徑。某個(gè)粒子的相鄰粒子即在其核函數(shù)控制半徑內(nèi)的粒子。

    粒子數(shù)密度表征的是流場中粒子的分布狀態(tài),保持粒子數(shù)密度恒定就是保證流體的不可壓縮性。該參數(shù)可以通過粒子在核函數(shù)的控制域內(nèi)與其周圍相鄰粒子的核函數(shù)數(shù)值的疊加表示,即:

    式中,N是粒子i的相鄰粒子總數(shù),ri和rj分別是粒子i和粒子j的位置矢量。

    梯度算子和拉普拉斯算子可以分別表示為:

    1.2 邊界條件

    邊界條件可分為兩種類型,即自由表面邊界條件和固壁邊界條件。

    對于自由表面邊界條件,首先要進(jìn)行自由表面粒子的判定,原始的判定方法是根據(jù)粒子數(shù)密度進(jìn)行判定,具體判定公式如下[9]:

    式中,β是一個(gè)小于1的系數(shù),一般取β=0.97。

    上面根據(jù)粒子數(shù)判定自由液面粒子的方法常常導(dǎo)致誤判現(xiàn)象的發(fā)生,為了克服這個(gè)問題,改進(jìn)的MPS方法中采用了一種弧度判定方法,這種方法通過對某個(gè)粒子的相鄰粒子相對該粒子的弧度進(jìn)行疊加,如果疊加后的弧度大于等于2 π并且在(0,2 π)區(qū)間內(nèi)無間斷,則該粒子判定為流體內(nèi)部粒子,否則為自由表面粒子,這種方法已證明可以有效避免粒子的誤判現(xiàn)象[8]。

    對于自由表面粒子,其壓力賦為大氣壓即p=patm=0。

    至于固壁條件,首先,固壁用三層邊界粒子來表示,靠近流體的那一層為流體邊界粒子,參與壓力的計(jì)算;外層的兩層粒子稱為固體邊界粒子,不參與壓力計(jì)算。另外,為了防止內(nèi)部的流體粒子穿透固體邊界,當(dāng)流體粒子靠近邊界時(shí),會(huì)收到固壁粒子的一個(gè)排斥力,具體計(jì)算如下[10]:

    其中,r0為初始粒子間距;D=5 gh,h為水柱的初始高度;P1=12,P2=4。

    1.3 計(jì)算流程

    (1)根據(jù)第n個(gè)時(shí)間步的信息,通過粘性項(xiàng)和外力項(xiàng)顯式預(yù)估粒子的中間速度和中間位置信息,即:

    (2)通過式(3)計(jì)算粒子的中間粒子數(shù)密度信息n*,然后隱式求解Poisson方程,得到pn+1:

    (3)根據(jù)pn+1對中間速度和中間位置信息進(jìn)行修正,得到第n+1時(shí)間步的速度和位置信息:

    (4)為了避免內(nèi)部流體粒子的聚集現(xiàn)象,引入了一種碰撞模型,這是改進(jìn)MPS方法的另外一個(gè)改進(jìn)的措施。即當(dāng)兩個(gè)粒子之間的間距滿足d<βL0時(shí),應(yīng)用的碰撞模型如下[8]:

    2 固定側(cè)壁潰壩現(xiàn)象

    選取的固定側(cè)壁潰壩模型如圖2所示,水箱的長度為3.22 m,高度為2 m,水柱初始高度0.6 m,位于水柱對面墻壁上的一個(gè)點(diǎn)P(3.22 m,0.16 m)作為壓力監(jiān)測點(diǎn)。離散模型的數(shù)據(jù)為:流體粒子7 200個(gè),流體邊界粒子723個(gè),固體邊界粒子2 196個(gè),共計(jì)10 116個(gè)粒子。

    圖2 潰壩模型及離散圖Fig.2 Model of the dam break and the corresponding dicretisizing model

    圖3 潰壩過程及壓力分布情況Fig.3 The process of dam break and the pressure distribution

    圖3給出了改進(jìn)MPS方法模擬的潰壩過程中的六個(gè)典型時(shí)刻的流動(dòng)現(xiàn)象及對應(yīng)的壓力分布。從圖中可以看出,當(dāng)潰壩現(xiàn)象發(fā)生后,水柱在重力的作用下發(fā)生坍塌形成了一個(gè)自由表面波或稱為潰壩波,該表面波在慣性的作用下向右傳播,然后與右側(cè)壁面發(fā)生砰擊現(xiàn)象,之后出現(xiàn)了自由表面的翻卷、破碎和飛濺現(xiàn)象,之后接著反向向左傳播,能量逐漸減弱直至達(dá)到新的靜平衡狀態(tài)。

    圖4(a)給出的是潰壩波前端位置隨時(shí)間的變化曲線,圖4(b)給出了P點(diǎn)的壓力隨時(shí)間的變化情況,這兩幅圖均與對應(yīng)的實(shí)驗(yàn)值[11]進(jìn)行了比較。從圖中可以看出,潰壩波大約在T=2.5附近到達(dá)右面的側(cè)壁并緊接著發(fā)生了砰擊作用,整個(gè)過程共發(fā)生了兩次比較明顯的砰擊現(xiàn)象。另外,改進(jìn)的MPS方法的計(jì)算結(jié)果較為平滑,與實(shí)驗(yàn)值吻合更好。但是,改進(jìn)的MPS方法的計(jì)算結(jié)果,無論是潰壩波前端位置的變化還是P點(diǎn)的壓力時(shí)歷值,相對于實(shí)驗(yàn)值存在一定的延時(shí),這種情況可能是由于數(shù)值模擬中運(yùn)動(dòng)學(xué)粘性系數(shù)的取值較大造成的。

    圖4 潰壩波前端位置和P點(diǎn)壓力的時(shí)歷曲線Fig.4 Time history of the leading edge of collapsed water column wave and the pressure at the point of P

    3 彈性支撐側(cè)壁潰壩現(xiàn)象

    實(shí)際工程中,當(dāng)發(fā)生潰壩時(shí),由于潰壩波的巨大砰擊力的作用,會(huì)導(dǎo)致一些建筑物的移動(dòng)或旋轉(zhuǎn)。因此,本文選取如圖5所示的具有彈性支撐側(cè)壁的潰壩模型,水柱右面的墻壁為彈性支撐,在潰壩波的砰擊力作用下將圍繞彈性支撐點(diǎn)做旋轉(zhuǎn)運(yùn)動(dòng),其轉(zhuǎn)動(dòng)方程為:

    式中,I為慣性矩;M為轉(zhuǎn)動(dòng)體的質(zhì)量;H為轉(zhuǎn)動(dòng)體質(zhì)心到旋轉(zhuǎn)點(diǎn)的距離;K為彈性支撐的剛度系數(shù);θ為彈性支撐側(cè)壁的轉(zhuǎn)角;Mf為流體壓力作用在彈性支撐側(cè)壁上的合力矩,可按照下式計(jì)算:

    其中,pi為彈性支撐側(cè)壁粒子i處的壓力;N為彈性支撐側(cè)壁離散的粒子總數(shù);l為粒子i距彈性支撐點(diǎn)的距離;dl可取為粒子初始分布間距。

    具體尺寸選取如圖2所示的潰壩現(xiàn)象的物理模型,只是右側(cè)固壁變成了彈性支撐側(cè)壁。在整個(gè)數(shù)值模擬中運(yùn)動(dòng)性粘性系數(shù)v取0.001 m/s2,計(jì)算時(shí)間步長dt為0.000 5 s。圖6給出了具有彈性支撐側(cè)壁的潰壩現(xiàn)象模擬結(jié)果,包括6個(gè)在整個(gè)潰壩現(xiàn)象中比較典型時(shí)刻點(diǎn)的圖像,分布對應(yīng)T=3.859 9、6.147 3、7.576 9、9.649 8、11.293 8和21.444 0。圖7顯示了在整個(gè)潰壩過程中彈性支撐側(cè)壁的旋轉(zhuǎn)角度θ隨時(shí)間的變化情況,圖8中給出了P點(diǎn)壓力隨時(shí)間的變化規(guī)律,圖9比較了右側(cè)側(cè)壁在固定和彈性支撐兩種情況下的所受流體合力矩隨時(shí)間的變化情況。

    圖5 彈性支撐側(cè)壁的潰壩模型示意圖Fig.5 Illustration of dam break with elastic-supported rigid wall

    圖6 具有彈性支撐側(cè)壁的潰壩現(xiàn)象Fig.6 Phenomenon of dam break with elastic-supported side wall

    結(jié)合圖6-9可以看出,當(dāng)潰壩現(xiàn)象發(fā)生之后,水柱坍塌并形成一個(gè)潰壩波,在T=3.8附近到達(dá)右側(cè)的墻壁,右側(cè)側(cè)壁所受流體的合力矩從0開始增長,在該力矩的作用下,彈性支撐側(cè)壁開始逐漸向右側(cè)旋轉(zhuǎn);緊接著潰壩波與右側(cè)側(cè)壁發(fā)生了比較激烈的抨擊現(xiàn)象,右側(cè)側(cè)壁所受流體的合力矩在T=8左右達(dá)到最大值,與此同時(shí)其轉(zhuǎn)角也達(dá)到最大值0.25左右;此后,潰壩波在左右兩個(gè)墻壁之間反復(fù)的傳播直至靜止?fàn)顟B(tài),右側(cè)側(cè)壁的旋轉(zhuǎn)角度也在這個(gè)往復(fù)震蕩過程中逐漸減小至0。

    圖7 轉(zhuǎn)角θ隨時(shí)間的變化曲線Fig.7 The time history of the rotation angle

    圖8 P點(diǎn)的壓力隨時(shí)間的變化曲線Fig.8 The time history of the pressure at the point of P

    另外,結(jié)合前節(jié)固定側(cè)壁的潰壩現(xiàn)象,彈性支撐側(cè)壁的潰壩模型主要有以下特征:首先,具有彈性支撐側(cè)壁的潰壩現(xiàn)象從潰壩發(fā)生到重新恢復(fù)平衡狀態(tài)所消耗的時(shí)間要大于固定側(cè)壁所需的時(shí)間,前者是在T大于20以后才達(dá)到靜平衡狀態(tài),后者是在T=11左右。與固定側(cè)壁的潰壩模型相比,彈性支撐側(cè)壁的潰壩模型的整個(gè)潰壩過程和側(cè)壁上流體壓力合力矩的峰值的出現(xiàn)均具有一定的延遲性。其次,這兩種潰壩模型右側(cè)側(cè)壁的壓力時(shí)歷曲線有明顯的不同。雖然這兩種潰壩模型的最大抨擊壓力峰值均在1.2左右,但固定側(cè)壁模型中有兩次比較明顯的砰擊現(xiàn)象,而對于彈性支撐側(cè)壁的潰壩模型,當(dāng)潰壩波到達(dá)右側(cè)側(cè)壁發(fā)生砰擊現(xiàn)象時(shí),右側(cè)墻壁在砰擊力矩的作用下向右側(cè)旋轉(zhuǎn),彈性支撐起到了明顯的延遲用,從而沒有明顯的第一次砰擊現(xiàn)象。

    圖9 右側(cè)側(cè)壁的流體合力矩隨時(shí)間的變化曲線Fig.9 The fluid force momentum on the right-side wall

    4 小結(jié)

    本文基于改進(jìn)的MPS方法,采用Fortran語言進(jìn)行自主編程,分別模擬了固定側(cè)壁潰壩現(xiàn)象和彈性支撐側(cè)壁潰壩現(xiàn)象。研究證明,改進(jìn)的MPS方法可以成功地模擬潰壩現(xiàn)象中自由液面的翻卷、破碎和飛濺現(xiàn)象,驗(yàn)證了該方法在處理這種大變形自由液面流動(dòng)問題中的優(yōu)勢。另一方面,通過這兩種潰壩模型的對比發(fā)現(xiàn),雖然兩種模型的最大砰擊壓力相近,但固定側(cè)壁的潰壩模型有兩次明顯的砰擊現(xiàn)象,彈性支撐側(cè)壁的潰壩現(xiàn)象只有一次明顯的砰擊現(xiàn)象。其次,由于彈性支撐側(cè)壁的緩沖作用,整個(gè)潰壩過程的進(jìn)行出現(xiàn)了顯著的延遲效應(yīng)。

    猜你喜歡
    水柱潰壩側(cè)壁
    探探鯨的水柱
    高速公路隧道洞口淺埋段施工中的雙側(cè)壁導(dǎo)坑法
    柳鋼6號(hào)高爐應(yīng)對爐缸側(cè)壁溫度超高的護(hù)爐實(shí)踐
    四川冶金(2019年4期)2019-11-18 09:30:50
    越野車側(cè)壁加強(qiáng)筋板工藝改進(jìn)
    模具制造(2019年3期)2019-06-06 02:10:52
    Run through the rain
    水柱有“魔力”
    徐家河尾礦庫潰壩分析
    潰壩涌浪及其對重力壩影響的數(shù)值模擬
    水柱測量中的水下滑翔機(jī)轉(zhuǎn)向性能
    潰壩波對單橋墩作用水力特性研究
    麻豆成人午夜福利视频| 激情五月婷婷亚洲| 日本黄大片高清| 在线播放无遮挡| 亚洲美女搞黄在线观看| 91在线精品国自产拍蜜月| 高清日韩中文字幕在线| 国产一区二区在线观看日韩| 国产乱人视频| 国产精品久久久久久久电影| 免费人成在线观看视频色| 久久久久久久久久久丰满| 国产在视频线精品| 美女福利国产在线 | 国产精品蜜桃在线观看| 亚洲人与动物交配视频| 高清欧美精品videossex| 夜夜骑夜夜射夜夜干| 国产v大片淫在线免费观看| 夜夜爽夜夜爽视频| 国产成人午夜福利电影在线观看| 国产精品国产av在线观看| 成人一区二区视频在线观看| 国产一区二区三区综合在线观看 | 中文欧美无线码| 亚洲激情五月婷婷啪啪| 七月丁香在线播放| 成人18禁高潮啪啪吃奶动态图 | 男女无遮挡免费网站观看| 国产一区二区三区综合在线观看 | 国产精品一及| 国产精品欧美亚洲77777| 久久青草综合色| 一级毛片久久久久久久久女| 午夜激情久久久久久久| 日本免费在线观看一区| 丰满迷人的少妇在线观看| h日本视频在线播放| 成年av动漫网址| 久久99热这里只频精品6学生| 日韩精品有码人妻一区| 网址你懂的国产日韩在线| 日韩,欧美,国产一区二区三区| 在线观看一区二区三区激情| 成人美女网站在线观看视频| 日韩电影二区| 丝袜脚勾引网站| 欧美丝袜亚洲另类| 99热6这里只有精品| 国语对白做爰xxxⅹ性视频网站| 在线观看美女被高潮喷水网站| 亚洲精品乱久久久久久| 观看av在线不卡| 亚洲av欧美aⅴ国产| 亚洲中文av在线| 久久女婷五月综合色啪小说| 美女脱内裤让男人舔精品视频| 色网站视频免费| 亚洲成色77777| 五月开心婷婷网| 交换朋友夫妻互换小说| 亚洲,欧美,日韩| 国产成人一区二区在线| 亚洲无线观看免费| 久久精品国产亚洲av天美| 尤物成人国产欧美一区二区三区| 久久精品国产a三级三级三级| 超碰97精品在线观看| 欧美日韩国产mv在线观看视频 | 交换朋友夫妻互换小说| 一个人看视频在线观看www免费| 欧美+日韩+精品| 欧美激情极品国产一区二区三区 | 中文字幕亚洲精品专区| 美女视频免费永久观看网站| 日韩视频在线欧美| 精品一品国产午夜福利视频| 亚洲精华国产精华液的使用体验| 国产成人a区在线观看| 直男gayav资源| 亚洲精品一区蜜桃| 国产亚洲av片在线观看秒播厂| 国产成人精品一,二区| 国产av一区二区精品久久 | 亚洲精品视频女| 全区人妻精品视频| 国产有黄有色有爽视频| 人妻制服诱惑在线中文字幕| 国产精品不卡视频一区二区| 亚洲精品一二三| 另类亚洲欧美激情| 少妇的逼水好多| 久久婷婷青草| 久久久欧美国产精品| 一个人看视频在线观看www免费| 久久鲁丝午夜福利片| 亚洲无线观看免费| 伊人久久精品亚洲午夜| 国产白丝娇喘喷水9色精品| 欧美高清成人免费视频www| 亚洲精品一二三| 亚洲欧美成人精品一区二区| 在线观看免费视频网站a站| 国产精品不卡视频一区二区| 赤兔流量卡办理| 欧美老熟妇乱子伦牲交| av视频免费观看在线观看| 99久久精品一区二区三区| 久久国产精品男人的天堂亚洲 | 伊人久久国产一区二区| 欧美 日韩 精品 国产| 国产69精品久久久久777片| 狂野欧美激情性bbbbbb| 国模一区二区三区四区视频| 欧美成人一区二区免费高清观看| 一区二区三区免费毛片| 在线观看免费高清a一片| 午夜福利视频精品| 午夜福利视频精品| 97热精品久久久久久| 一级毛片我不卡| 最近手机中文字幕大全| videossex国产| 午夜老司机福利剧场| av不卡在线播放| 在现免费观看毛片| 插阴视频在线观看视频| 日韩欧美精品免费久久| 最近手机中文字幕大全| 美女cb高潮喷水在线观看| 高清不卡的av网站| 久久国产亚洲av麻豆专区| 卡戴珊不雅视频在线播放| 日日摸夜夜添夜夜爱| 亚洲国产最新在线播放| 久久综合国产亚洲精品| 偷拍熟女少妇极品色| 超碰av人人做人人爽久久| 我要看日韩黄色一级片| 丰满乱子伦码专区| 国精品久久久久久国模美| 91在线精品国自产拍蜜月| 亚洲aⅴ乱码一区二区在线播放| 一区二区三区免费毛片| 好男人视频免费观看在线| 美女福利国产在线 | 在线免费观看不下载黄p国产| 久久99热6这里只有精品| 欧美成人精品欧美一级黄| 老熟女久久久| 国产大屁股一区二区在线视频| 欧美成人a在线观看| 91精品国产国语对白视频| 亚洲欧美日韩卡通动漫| 亚洲高清免费不卡视频| 18禁裸乳无遮挡免费网站照片| 夜夜骑夜夜射夜夜干| 国产片特级美女逼逼视频| 麻豆成人午夜福利视频| 亚洲欧美一区二区三区国产| 水蜜桃什么品种好| 亚洲四区av| 亚洲欧美一区二区三区国产| 韩国av在线不卡| 亚洲精品中文字幕在线视频 | 日韩中文字幕视频在线看片 | av在线老鸭窝| 国产欧美日韩精品一区二区| 亚洲国产精品999| 国产成人精品福利久久| 国产高潮美女av| 26uuu在线亚洲综合色| 一区二区av电影网| 成人国产麻豆网| 亚洲精品乱久久久久久| 爱豆传媒免费全集在线观看| 亚洲自偷自拍三级| 舔av片在线| 免费观看无遮挡的男女| 啦啦啦视频在线资源免费观看| 精品久久久精品久久久| 简卡轻食公司| 超碰97精品在线观看| 丝袜喷水一区| 人妻 亚洲 视频| 大香蕉久久网| 日韩欧美 国产精品| 日本午夜av视频| 各种免费的搞黄视频| 九色成人免费人妻av| 91在线精品国自产拍蜜月| 国产精品嫩草影院av在线观看| 久热久热在线精品观看| 青春草视频在线免费观看| 亚洲,欧美,日韩| av.在线天堂| 午夜福利在线在线| 又粗又硬又长又爽又黄的视频| 日韩强制内射视频| 午夜福利视频精品| 在线免费十八禁| 国产91av在线免费观看| 免费人成在线观看视频色| 赤兔流量卡办理| 亚洲色图av天堂| 日韩人妻高清精品专区| 26uuu在线亚洲综合色| 精品一区在线观看国产| 高清不卡的av网站| 国产 精品1| tube8黄色片| 日韩在线高清观看一区二区三区| 人体艺术视频欧美日本| av在线观看视频网站免费| 各种免费的搞黄视频| 国产免费视频播放在线视频| 国内精品宾馆在线| 99热全是精品| 国产一区有黄有色的免费视频| 18+在线观看网站| 久久久久久久久久久免费av| 国产免费福利视频在线观看| 亚洲一区二区三区欧美精品| 日韩不卡一区二区三区视频在线| 搡老乐熟女国产| 亚洲婷婷狠狠爱综合网| 亚洲人与动物交配视频| 亚洲不卡免费看| 久久这里有精品视频免费| 黄色一级大片看看| 国产高潮美女av| 黄片wwwwww| 欧美少妇被猛烈插入视频| 成人漫画全彩无遮挡| 熟女人妻精品中文字幕| av线在线观看网站| 80岁老熟妇乱子伦牲交| 三级经典国产精品| 午夜福利在线在线| 国产精品一区二区在线不卡| 日本欧美国产在线视频| 极品教师在线视频| 国产色婷婷99| 精品酒店卫生间| 婷婷色麻豆天堂久久| 最近最新中文字幕免费大全7| 成人影院久久| 国产中年淑女户外野战色| 少妇熟女欧美另类| 极品少妇高潮喷水抽搐| 尾随美女入室| 在线天堂最新版资源| 亚洲欧美精品专区久久| 精品久久久久久久久亚洲| 亚洲欧洲日产国产| 丰满少妇做爰视频| 国产无遮挡羞羞视频在线观看| 免费少妇av软件| 女人久久www免费人成看片| 毛片一级片免费看久久久久| 亚洲国产色片| 欧美丝袜亚洲另类| 少妇精品久久久久久久| 人人妻人人爽人人添夜夜欢视频 | 丝袜喷水一区| 大香蕉久久网| av在线观看视频网站免费| 少妇的逼好多水| 欧美日韩国产mv在线观看视频 | 亚洲欧洲日产国产| 国产美女午夜福利| 精品酒店卫生间| av又黄又爽大尺度在线免费看| 国产中年淑女户外野战色| 亚洲久久久国产精品| 欧美精品人与动牲交sv欧美| 国产亚洲最大av| 18禁动态无遮挡网站| 精品国产露脸久久av麻豆| 色哟哟·www| 久热久热在线精品观看| 美女脱内裤让男人舔精品视频| 亚洲自偷自拍三级| 日韩 亚洲 欧美在线| 亚洲精华国产精华液的使用体验| 久久久午夜欧美精品| av国产免费在线观看| 欧美zozozo另类| 国产成人a区在线观看| 亚洲欧美一区二区三区黑人 | 欧美最新免费一区二区三区| 国产精品一及| 秋霞在线观看毛片| 欧美+日韩+精品| 亚洲成人中文字幕在线播放| 亚洲aⅴ乱码一区二区在线播放| 国产欧美日韩一区二区三区在线 | 一本一本综合久久| 国产精品麻豆人妻色哟哟久久| .国产精品久久| 亚洲精品日韩在线中文字幕| 国产精品免费大片| 亚洲国产色片| 午夜视频国产福利| av播播在线观看一区| 偷拍熟女少妇极品色| 天堂俺去俺来也www色官网| 国产午夜精品一二区理论片| 一区二区三区精品91| 亚洲欧美日韩卡通动漫| 蜜桃久久精品国产亚洲av| 熟女电影av网| 看非洲黑人一级黄片| 精品久久久久久久末码| 精品亚洲成a人片在线观看 | 欧美成人a在线观看| 国产视频内射| 国产成人一区二区在线| 国产精品偷伦视频观看了| 国精品久久久久久国模美| 女的被弄到高潮叫床怎么办| 日韩成人伦理影院| 亚洲欧美精品专区久久| 女人十人毛片免费观看3o分钟| 久久99热6这里只有精品| 日韩,欧美,国产一区二区三区| 街头女战士在线观看网站| 啦啦啦中文免费视频观看日本| 如何舔出高潮| 黄色一级大片看看| 婷婷色综合大香蕉| 丰满少妇做爰视频| 日韩电影二区| 五月伊人婷婷丁香| 丰满少妇做爰视频| 欧美97在线视频| 日本黄大片高清| 国国产精品蜜臀av免费| 春色校园在线视频观看| 日韩电影二区| 交换朋友夫妻互换小说| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲欧美清纯卡通| 少妇精品久久久久久久| 91狼人影院| 黄片wwwwww| 国国产精品蜜臀av免费| 好男人视频免费观看在线| 欧美成人午夜免费资源| h日本视频在线播放| 九九久久精品国产亚洲av麻豆| 一区二区三区免费毛片| 中文在线观看免费www的网站| 肉色欧美久久久久久久蜜桃| 亚州av有码| 欧美xxxx性猛交bbbb| 91狼人影院| 一区在线观看完整版| 在线观看免费高清a一片| 亚洲人成网站在线播| 亚洲国产日韩一区二区| 亚洲欧美成人综合另类久久久| 97热精品久久久久久| av视频免费观看在线观看| 22中文网久久字幕| 欧美一区二区亚洲| 亚洲成人一二三区av| 久久久久国产网址| 一区二区av电影网| 性色avwww在线观看| 亚洲人成网站高清观看| 永久免费av网站大全| 激情五月婷婷亚洲| 国产精品一二三区在线看| 亚洲精华国产精华液的使用体验| av免费观看日本| 精品久久久久久久末码| 免费黄频网站在线观看国产| 亚洲精品乱码久久久久久按摩| 嫩草影院新地址| 亚洲av在线观看美女高潮| 人妻夜夜爽99麻豆av| av一本久久久久| 性色av一级| 亚洲成色77777| 免费久久久久久久精品成人欧美视频 | 午夜精品国产一区二区电影| 亚洲精品亚洲一区二区| 六月丁香七月| 久热这里只有精品99| 99热网站在线观看| 97精品久久久久久久久久精品| 在线精品无人区一区二区三 | 国产精品麻豆人妻色哟哟久久| 国产男女超爽视频在线观看| 最近最新中文字幕免费大全7| 中文乱码字字幕精品一区二区三区| 午夜免费男女啪啪视频观看| 欧美日韩在线观看h| 日本猛色少妇xxxxx猛交久久| 久久精品久久久久久久性| 黄色日韩在线| 亚洲精品乱久久久久久| 国产精品一区www在线观看| 久久亚洲国产成人精品v| 免费播放大片免费观看视频在线观看| 美女cb高潮喷水在线观看| 秋霞伦理黄片| 少妇人妻一区二区三区视频| 欧美成人a在线观看| 联通29元200g的流量卡| 日本午夜av视频| 免费观看的影片在线观看| 国产欧美日韩精品一区二区| 欧美日韩精品成人综合77777| 能在线免费看毛片的网站| 久热这里只有精品99| 亚洲欧美日韩无卡精品| 在线观看美女被高潮喷水网站| 新久久久久国产一级毛片| 亚洲精品日韩在线中文字幕| 啦啦啦中文免费视频观看日本| 精品人妻一区二区三区麻豆| 亚洲精品亚洲一区二区| 少妇猛男粗大的猛烈进出视频| 国产 一区 欧美 日韩| 国产永久视频网站| 午夜免费鲁丝| 亚洲第一av免费看| 亚洲真实伦在线观看| 毛片一级片免费看久久久久| 欧美一级a爱片免费观看看| 欧美日韩视频精品一区| 一级毛片我不卡| 成人高潮视频无遮挡免费网站| 国产男女内射视频| 国产爽快片一区二区三区| 一区二区三区精品91| 亚洲国产色片| 99精国产麻豆久久婷婷| 校园人妻丝袜中文字幕| 国产欧美日韩精品一区二区| 18禁裸乳无遮挡动漫免费视频| 18禁动态无遮挡网站| 国产女主播在线喷水免费视频网站| 91精品伊人久久大香线蕉| 日日摸夜夜添夜夜添av毛片| 爱豆传媒免费全集在线观看| 亚洲精品国产av蜜桃| av在线观看视频网站免费| 五月天丁香电影| 老司机影院毛片| 色视频在线一区二区三区| 久久婷婷青草| 一级毛片 在线播放| 女的被弄到高潮叫床怎么办| 午夜日本视频在线| 夜夜骑夜夜射夜夜干| 成人亚洲精品一区在线观看 | 久久久色成人| 国产成人精品福利久久| 男女国产视频网站| 免费黄频网站在线观看国产| 免费观看a级毛片全部| 亚洲精品第二区| 亚洲精品aⅴ在线观看| 国产精品久久久久久av不卡| 久久精品久久久久久久性| a级毛色黄片| av.在线天堂| 亚洲欧美日韩卡通动漫| 亚洲人成网站在线观看播放| 久久久久久伊人网av| 国产精品伦人一区二区| 午夜福利网站1000一区二区三区| 一二三四中文在线观看免费高清| 人妻 亚洲 视频| 18禁裸乳无遮挡动漫免费视频| 欧美成人精品欧美一级黄| 国产免费一级a男人的天堂| 99久久精品国产国产毛片| 亚洲成人av在线免费| 毛片一级片免费看久久久久| 欧美老熟妇乱子伦牲交| 男女国产视频网站| 在线观看免费视频网站a站| 色哟哟·www| 美女中出高潮动态图| 国产一区二区三区av在线| 国产一级毛片在线| 成人黄色视频免费在线看| 久久久久久久大尺度免费视频| 日本av免费视频播放| 少妇 在线观看| 亚洲精品自拍成人| 免费黄色在线免费观看| 交换朋友夫妻互换小说| 黄色视频在线播放观看不卡| 人体艺术视频欧美日本| 人人妻人人添人人爽欧美一区卜 | 日本黄色日本黄色录像| 激情五月婷婷亚洲| 精品一区二区三卡| 国产亚洲91精品色在线| 青春草视频在线免费观看| 亚州av有码| av在线app专区| 搡女人真爽免费视频火全软件| 国产精品久久久久久精品古装| 国产极品天堂在线| 简卡轻食公司| 久久热精品热| 大码成人一级视频| 亚洲图色成人| 成人黄色视频免费在线看| 久久久久国产精品人妻一区二区| 亚洲国产高清在线一区二区三| 国产成人a∨麻豆精品| 欧美激情国产日韩精品一区| 黑人猛操日本美女一级片| 国产日韩欧美在线精品| 久久久久久伊人网av| 久久久久久久国产电影| 国产精品国产三级专区第一集| 国产精品嫩草影院av在线观看| 人人妻人人澡人人爽人人夜夜| 我要看日韩黄色一级片| 成人一区二区视频在线观看| 日韩制服骚丝袜av| 国产精品久久久久久久久免| 亚洲av中文av极速乱| 大香蕉97超碰在线| 99精国产麻豆久久婷婷| 91久久精品国产一区二区成人| kizo精华| 男女边吃奶边做爰视频| 亚洲国产毛片av蜜桃av| 免费播放大片免费观看视频在线观看| 久久久久久久大尺度免费视频| 少妇人妻 视频| 乱码一卡2卡4卡精品| 黑人高潮一二区| 男女无遮挡免费网站观看| 免费av中文字幕在线| 国产伦理片在线播放av一区| 一级二级三级毛片免费看| 日本黄色日本黄色录像| 久久国产亚洲av麻豆专区| 肉色欧美久久久久久久蜜桃| 亚洲av免费高清在线观看| 久久久久久九九精品二区国产| 六月丁香七月| 高清午夜精品一区二区三区| 亚洲欧美中文字幕日韩二区| 午夜福利在线在线| 啦啦啦啦在线视频资源| 婷婷色综合大香蕉| 晚上一个人看的免费电影| 精品久久久久久久久亚洲| 久久久久久久大尺度免费视频| av免费观看日本| 在线亚洲精品国产二区图片欧美 | 国产日韩欧美在线精品| 2018国产大陆天天弄谢| 中文字幕免费在线视频6| 草草在线视频免费看| 国产精品麻豆人妻色哟哟久久| 亚洲内射少妇av| 午夜免费男女啪啪视频观看| 国产精品爽爽va在线观看网站| 夜夜爽夜夜爽视频| 高清av免费在线| 一级毛片aaaaaa免费看小| 在线观看三级黄色| 在线观看免费视频网站a站| 亚洲精品国产色婷婷电影| 成人高潮视频无遮挡免费网站| 大片免费播放器 马上看| 成人亚洲欧美一区二区av| 2021少妇久久久久久久久久久| 亚洲av欧美aⅴ国产| 最近2019中文字幕mv第一页| 亚洲精品久久久久久婷婷小说| 国产精品一二三区在线看| 国产精品福利在线免费观看| 国产精品久久久久久精品古装| 亚洲欧美日韩无卡精品| 女的被弄到高潮叫床怎么办| 国产一区二区在线观看日韩| 免费黄网站久久成人精品| 亚洲av.av天堂| 国产精品麻豆人妻色哟哟久久| 国产精品人妻久久久影院| 一个人免费看片子| 亚洲av福利一区| 国产午夜精品一二区理论片| 在线观看美女被高潮喷水网站| 免费看光身美女| 女人十人毛片免费观看3o分钟| 久久久久久久久久久免费av| av国产久精品久网站免费入址| 九九在线视频观看精品| 久久国内精品自在自线图片| 永久免费av网站大全| 久久精品夜色国产| 亚洲美女黄色视频免费看| 精品亚洲乱码少妇综合久久| 91久久精品国产一区二区三区| 国产一区有黄有色的免费视频| 日本猛色少妇xxxxx猛交久久| 亚洲av中文字字幕乱码综合| 久久人人爽人人爽人人片va| 另类亚洲欧美激情| 观看美女的网站| 欧美日韩综合久久久久久|