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

    考慮空氣卷入效應(yīng)的甲板上浪數(shù)值模擬

    2018-03-01 02:46:50
    船海工程 2018年1期
    關(guān)鍵詞:壓縮空氣甲板峰值

    ,

    (上海交通大學(xué) a.海洋工程國家重點(diǎn)實(shí)驗(yàn)室;b.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海 200240)

    甲板上浪屬于強(qiáng)非線性問題,早期試驗(yàn)觀察了上浪的整個(gè)過程,并且記錄了上浪過程中甲板上水體的高度和上浪對(duì)上部結(jié)構(gòu)的沖擊載荷[1- 2]。近年有學(xué)者通過在Fluent軟件平臺(tái)上建立起具有造波和消波功能的二維數(shù)值波浪水槽,實(shí)現(xiàn)了對(duì)Greco試驗(yàn)的數(shù)值模擬,最終獲得上浪過程的波高、壓力的模擬結(jié)果,但第一次上浪的水體高度時(shí)歷和壓力時(shí)歷結(jié)果偏差較大[3]。同時(shí),有學(xué)者通過自主開發(fā)的naoe- FOAM- SJTU求解器進(jìn)行二維甲板上浪的數(shù)值模擬,最終結(jié)果和文獻(xiàn)[3]中結(jié)果相似[4]??梢姸S模擬不能完全反映Greco的試驗(yàn),數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果存在一定的偏差,模擬方法存在一定的局限性。

    在入水砰擊的相關(guān)研究中,有學(xué)者發(fā)現(xiàn)由于物體快速進(jìn)入水中,導(dǎo)致一部分空氣未來得及逃逸形成空氣泡,可壓縮空氣泡對(duì)沖擊載荷的緩沖作用稱為空氣墊效應(yīng)[5]。通過數(shù)值方法研究空氣墊對(duì)砰擊壓力的影響,可以發(fā)現(xiàn)砰擊壓力的一些特點(diǎn),但相關(guān)研究中并未與不可壓縮空氣比較[6]。在三體船入水砰擊壓力的研究中發(fā)現(xiàn),空氣墊能夠降低相應(yīng)區(qū)域的壓力峰值[7],但入水砰擊與甲板上浪的特點(diǎn)是不一樣的,空氣墊的影響也不同,故而需要探索新的方法研究上浪空氣墊的影響。為此利用商業(yè)CFD軟件StarCCM+來進(jìn)行三維尺度下的數(shù)值模擬。

    1 甲板上浪數(shù)值模擬

    由于試驗(yàn)水槽較長,獲得穩(wěn)定的波浪的時(shí)間也比較久,如果數(shù)值模擬采用相同的方式會(huì)導(dǎo)致數(shù)值模型太大,數(shù)值計(jì)算需要模擬的時(shí)間也較長。所以,為了提高數(shù)值模擬的效率,通過縮短計(jì)算區(qū)域,設(shè)置合適的參數(shù)減少穩(wěn)定波浪生成的時(shí)間來降低數(shù)值模擬計(jì)算耗時(shí)。同時(shí)通過設(shè)置空氣的屬性,利用真實(shí)氣體模型來模擬可壓縮空氣,探索可壓縮空氣泡對(duì)砰擊載荷的影響。

    1.1 試驗(yàn)描述

    針對(duì)Greco的FPSO甲板上浪試驗(yàn)進(jìn)行模擬計(jì)算,試驗(yàn)用水槽分為造波區(qū)、消波區(qū)、工作區(qū)和后端消波區(qū),整個(gè)計(jì)算區(qū)域長達(dá)13.5 m。在水槽距造波板0.8 m和5.436 m的位置還布置了測(cè)量波高的傳感器WP1和WP2,用來監(jiān)測(cè)造波情況和沖擊甲板的波浪情況。試驗(yàn)用模型甲板上布置有WL1、WL2、WL3三個(gè)用來測(cè)量波高的傳感器,當(dāng)應(yīng)用上部擋板進(jìn)行試驗(yàn)時(shí),在上部擋板距離水面62 mm的位置水平布置pr1、pr3兩個(gè)壓力傳感器。當(dāng)沒有上部擋板時(shí),主要觀察WL1、WL2、WL3三個(gè)測(cè)點(diǎn)上浪時(shí)的波高時(shí)歷以及上浪的水體流動(dòng)過程。當(dāng)有上部擋板時(shí),主要觀察pr1與pr3的壓力時(shí)歷。

    1.2 數(shù)值計(jì)算方法與模型

    Greco的試驗(yàn)為了得到比較好的造波、消波效果,將試驗(yàn)船模布置在距離造波板5.54 m的位置,同時(shí)還設(shè)置了2 m的消波區(qū)域,這使得整個(gè)試驗(yàn)水槽長達(dá)13.5 m。二維數(shù)值模擬采用和試驗(yàn)相同的尺度,但二維數(shù)值模擬網(wǎng)格較少,對(duì)計(jì)算效率影響不大。

    二維數(shù)值模擬和試驗(yàn)在WL1處水體高度時(shí)歷對(duì)比見圖1,整個(gè)試驗(yàn)過程歷時(shí)約9.0 s,其中,前7.0 s的時(shí)間都是在等待穩(wěn)定的波生成并傳播到試驗(yàn)船模附近,在兩個(gè)上浪之間也有約0.8 s左右的時(shí)間沒有任何響應(yīng)。這些時(shí)間的數(shù)值模擬結(jié)果都是不重要的,對(duì)這些時(shí)間的計(jì)算模擬會(huì)造成計(jì)算資源的浪費(fèi)和計(jì)算效率的降低。

    圖1 二維數(shù)值模擬和試驗(yàn)的WL1處水體高度

    進(jìn)行三維尺度模擬網(wǎng)格比二維模擬增加很多,如果采用相同尺度進(jìn)行建模,會(huì)導(dǎo)致計(jì)算時(shí)間成倍增加。同時(shí),在9.0 s的試驗(yàn)中,重點(diǎn)是上浪前后約1 s時(shí)間內(nèi)的響應(yīng)。所以為了減少波浪生成的時(shí)間,縮短模擬計(jì)算的范圍,根據(jù)StarCCM+軟件的功能特性,對(duì)計(jì)算模型進(jìn)行簡化。

    StarCCM+軟件自帶造波和消波功能模塊,不再需要單獨(dú)劃分造波區(qū)和消波區(qū),計(jì)算時(shí)將造波區(qū)、前端消波區(qū)和后端消波區(qū)取消,只保留工作區(qū)。為了有較好的消波效果,工作區(qū)后端增加0.5 m,使整個(gè)建模區(qū)域縮短到8 m,約減小40%。

    由于計(jì)算區(qū)域不一致,如果用和試驗(yàn)值一樣的波高參數(shù)進(jìn)行初始化,會(huì)導(dǎo)致船艏的波浪時(shí)歷和試驗(yàn)值有較大偏差。所以,在進(jìn)行數(shù)值模擬時(shí)相應(yīng)調(diào)整了初始時(shí)刻的波浪參數(shù),保證船首波浪的大小與試驗(yàn)值一致。通過這些調(diào)整,整個(gè)過程只需要模擬1.5~2.0 s便可得到需要的結(jié)果。

    模型造波端設(shè)置為inlet邊界條件,消波端設(shè)置為outlet邊界條件,其他邊界設(shè)置成inlet邊界條件,主要用于模擬無限寬水面,消除壁面的影響[8]。在邊界部分應(yīng)用邊界層網(wǎng)格,在自由表面附近采用了較細(xì)的網(wǎng)格,同時(shí)對(duì)船體表面附近的網(wǎng)格進(jìn)行了細(xì)化。

    1.3 上浪的水體高度和壓力

    在模擬上浪水體高度時(shí),需要應(yīng)用不含上部擋板的模型進(jìn)行計(jì)算。由于第一次上浪時(shí)甲板上都是干的,因此可以保證數(shù)值模擬和試驗(yàn)時(shí)甲板上的流動(dòng)狀態(tài)是一致的,所以此時(shí)只需要控制波浪的大小與試驗(yàn)值一致,即控制wp2處的水體高度時(shí)歷與試驗(yàn)值一致。圖2即為數(shù)值模擬和試驗(yàn)時(shí)無上部擋板wp2處水體高度時(shí)歷對(duì)比圖,可以看出數(shù)值模擬時(shí)上浪的大小和試驗(yàn)是一致的。

    圖2 無上部擋板wp2測(cè)點(diǎn)位置水體高度

    根據(jù)數(shù)值模擬的計(jì)算結(jié)果,繪制如圖3所示的WL1、WL2、WL3三個(gè)測(cè)點(diǎn)上的水體高度時(shí)歷曲線??梢园l(fā)現(xiàn),數(shù)值計(jì)算的結(jié)果和試驗(yàn)的結(jié)果符合較好,而且,相比于文獻(xiàn)[4]二維模型的數(shù)值模擬,三維尺度下的模擬更加精確。

    圖3 WL1、WL2、WL3處上浪水體高度

    在模擬研究上浪時(shí)上部擋板pr1、pr3測(cè)點(diǎn)壓力時(shí),需應(yīng)用含上部擋板的模型進(jìn)行模擬計(jì)算。甲板流動(dòng)狀態(tài)在第一次上浪時(shí)是一致的,只需控制wp2處的水體高度時(shí)歷與試驗(yàn)值一致。圖4為數(shù)值模擬和試驗(yàn)時(shí)有上部擋板wp2處水體高度時(shí)歷對(duì)比圖,可以看出,上浪大小也是一致的。

    圖4 有上部擋板wp2測(cè)點(diǎn)位置水體高度

    根據(jù)數(shù)值模擬的計(jì)算結(jié)果,繪制壓力時(shí)歷曲線見圖5。從圖5可以看出,三維尺度下數(shù)值模擬的結(jié)果和試驗(yàn)值吻合得較好。文獻(xiàn)[3]二維數(shù)值模擬的結(jié)果偏大,但也能反映出大概的趨勢(shì)。同時(shí)可以發(fā)現(xiàn),在上浪的時(shí)候存在兩個(gè)峰值,第一個(gè)峰值是波浪第一次沖擊到上部擋板時(shí)候的響應(yīng)值,第二個(gè)峰值是水流到達(dá)頂端開始從上部擋板回落時(shí)候的響應(yīng)值,而且,由于三維尺度的影響,第二個(gè)峰值的響應(yīng)會(huì)較低。同時(shí),如圖6,在利用真實(shí)可壓縮氣體模型計(jì)算時(shí),得到的結(jié)果與試驗(yàn)結(jié)果更加符合,這也說明了考慮可壓縮性的數(shù)值模型的準(zhǔn)確性。

    圖6 pr1、pr3空氣可壓縮性對(duì)砰擊壓力的影響

    2 卷入空氣泡對(duì)甲板上砰擊壓力的影響與空氣可壓縮性效應(yīng)

    為了研究空氣可壓縮性對(duì)上浪砰擊壓力的影響,在甲板上布置了若干個(gè)壓力測(cè)點(diǎn),測(cè)點(diǎn)FR6、FR5、FR4按20 mm距離布置,其中FR6距船艏20 mm,測(cè)點(diǎn)FR3、FR2、FR1按40 mm距離布置,其中FR3距船艏100 mm。為了保證空氣不會(huì)逃逸,計(jì)算時(shí)將水槽的兩側(cè)邊界設(shè)置成Wall。由于試驗(yàn)中的空氣泡很小,很難研究卷入空氣泡對(duì)甲板上砰擊壓力的影響以及可壓縮和不可壓縮空氣對(duì)砰擊載荷的影響,故而在計(jì)算時(shí)將初始的波高調(diào)整為0.2 m,將干舷減少0.005 m,這樣能夠觀察到更加明顯的空氣泡。圖7即為甲板前端WP2測(cè)點(diǎn)水體高度時(shí)歷曲線。由圖7可見,甲板的上浪大小在考慮和不考慮空氣可壓縮性時(shí)是一致的。上浪水體卷入的空氣泡如圖8所示,可以看到,在空氣泡最初形成的時(shí)候,F(xiàn)R2-6測(cè)點(diǎn)均在空氣泡中,而FR1測(cè)點(diǎn)在其外。

    圖7 甲板前端WP2測(cè)點(diǎn)水體高度變化

    圖8 上浪水體卷入的空氣泡

    圖9繪制出了不同測(cè)點(diǎn)考慮空氣可壓縮性時(shí)的壓力時(shí)歷,并與不考慮可壓縮性時(shí)的結(jié)果繪制在一起進(jìn)行比較。

    圖9 可壓縮與不可壓縮空氣砰擊壓力變化

    在本次數(shù)值模擬中,上浪過程主要經(jīng)歷如下階段:第一階段水體上涌,在上端開始彎曲翻卷,并開始往甲板運(yùn)動(dòng);第二階段水體接觸到甲板,形成最初始的空氣泡;第三階段,接觸到甲板之后,水體往前流動(dòng),并將空氣泡壓縮;第四階段的流動(dòng)為水體的前移和氣泡的運(yùn)動(dòng),氣泡運(yùn)動(dòng)主要表現(xiàn)為前移、破碎和逃逸。

    前兩階段為卷波型甲板上浪模式的水體運(yùn)動(dòng)方式,這兩階段的特點(diǎn)主要取決于上浪模式,不同上浪模式的特點(diǎn)在文獻(xiàn)[2]中有詳細(xì)說明。第三階段時(shí),會(huì)形成第一個(gè)壓力峰值,水體的壓力和砰擊力是造成這個(gè)壓力峰值主要原因。第四階段時(shí),由于水體前移,甲板上的水體高度會(huì)降低,壓力值的總趨勢(shì)是降低的。但是由于此階段空氣泡會(huì)出現(xiàn)破碎和逃逸,這對(duì)空氣泡附近的區(qū)域壓力值會(huì)產(chǎn)生復(fù)雜的影響。

    觀察第三階段的第一個(gè)壓力峰的情況可以發(fā)現(xiàn),F(xiàn)R2~FR6測(cè)點(diǎn)的壓力峰值幾乎是一致的,這是因?yàn)檫@5個(gè)測(cè)點(diǎn)在初始階段都在被上浪水體包圍的空氣泡中,空氣泡中的壓力是一致的。不考慮空氣可壓縮性時(shí),壓力峰值約為920 Pa,而考慮了可壓縮性之后,壓力峰值增加到1 040 Pa。這是由于在上浪水體給空氣的壓力達(dá)到不可壓縮壓力峰值時(shí),可壓縮空氣還能被壓縮,從而氣體內(nèi)的壓力進(jìn)一步增大。而且,在不考慮空氣可壓縮性時(shí),整個(gè)壓力峰維持了約0.013 s;當(dāng)考慮了可壓縮性之后,壓力峰的時(shí)間增加到了0.019 s。

    第一個(gè)壓力峰過后,由于水體前移,甲板上水體高度減少,壓力值逐漸降低。過后,水體繼續(xù)壓縮空氣泡,出現(xiàn)第二個(gè)壓力峰。此時(shí),空氣泡并未發(fā)生破碎和逃逸,壓力峰的值和第一次的壓力峰值是相當(dāng)?shù)?。第二個(gè)峰過后,就會(huì)發(fā)生空氣泡的破碎和逃逸,在這個(gè)階段里,空氣泡內(nèi)、前和后的位置壓力值會(huì)出現(xiàn)不同的表現(xiàn)。

    在0.66~0.68 s之間空氣泡發(fā)生破碎,分成一大一小兩個(gè)空氣泡,前面的空氣泡較小。此時(shí),小空氣泡和大空氣泡都要往前移動(dòng),這就會(huì)在空氣泡后端形成一個(gè)高壓區(qū),可以發(fā)現(xiàn),小空氣泡后的FR1和FR2測(cè)點(diǎn)的壓力值會(huì)變大出現(xiàn)新的峰值,大空氣泡后的FR6測(cè)點(diǎn)也出現(xiàn)了壓力峰值。同樣的特性在0.73 s時(shí)的FR5測(cè)點(diǎn)、0.75 s時(shí)的FR4測(cè)點(diǎn)、0.78 s時(shí)的FR3測(cè)點(diǎn)和0.8 s時(shí)的FR2測(cè)點(diǎn)都有出現(xiàn),這些時(shí)刻都是氣泡逐漸前移,在空氣泡后方形成高壓區(qū)。如圖10,大空氣泡內(nèi)的壓力逐漸降低,這主要是由于甲板上方的水體高度逐漸降低。同時(shí),在氣泡移動(dòng)過程中,后方推動(dòng)空氣泡移動(dòng)的高壓區(qū)的壓力值也在逐漸減小。而空氣泡前的點(diǎn)的壓力主要與水體有關(guān),且較空氣泡內(nèi)的壓力小,這樣空氣泡前后的壓力差就會(huì)推動(dòng)空氣泡移動(dòng)。

    圖10 大空氣泡內(nèi)的壓力變化

    由圖9可見,與不可壓縮空氣計(jì)算的結(jié)果對(duì)比,可壓縮空氣模型在第四階段的壓力峰值要小于不可壓縮情況,且處于空氣泡中部的FR3、FR4、FR5三個(gè)測(cè)點(diǎn)的壓力值均下降10%左右。同時(shí)由圖10可見,考慮可壓縮空氣時(shí)大氣泡內(nèi)的壓力時(shí)歷變化平緩,而不可壓時(shí)候空氣泡內(nèi)的壓力呈現(xiàn)較大的波動(dòng)。這些是可壓縮空氣泡的緩沖效應(yīng)導(dǎo)致的。

    3 結(jié)論

    1)通過數(shù)值模擬研究空氣泡對(duì)甲板載荷的影響,發(fā)現(xiàn)在甲板上浪時(shí)會(huì)有空氣被卷入無法逃逸形成空氣泡,在這些空氣泡的前進(jìn)方向前端出現(xiàn)一個(gè)低壓區(qū),后端出現(xiàn)一個(gè)高壓區(qū),空氣泡也因此被推動(dòng)前進(jìn)。并且隨著上浪時(shí)空氣泡的移動(dòng),當(dāng)空氣泡經(jīng)過甲板某一區(qū)域后,空氣泡后端的高壓區(qū)會(huì)導(dǎo)致此區(qū)域的壓力達(dá)到峰值。

    2)通過引入真實(shí)氣體模型研究可壓縮空氣和不可壓縮空氣在上浪載荷上的差異,最終發(fā)現(xiàn),考慮空氣可壓縮性在不同時(shí)期會(huì)產(chǎn)生不同的影響,在上浪最初不可壓縮空氣相對(duì)于可壓縮空氣作用時(shí)間短,峰值壓力較小,在后續(xù)水體和空氣泡繼續(xù)前進(jìn)時(shí),可壓縮空氣的峰值壓力較小,并且峰值相對(duì)平緩,這充分體現(xiàn)了可壓縮空氣泡的緩沖作用。

    3)在三維尺度下進(jìn)行數(shù)值計(jì)算,數(shù)值計(jì)算結(jié)果和試驗(yàn)值吻合得較好,能夠較好地模擬出甲板上浪現(xiàn)象,而且相對(duì)于二維尺度下的模擬更加精確。

    4)上浪過程中卷入的空氣泡會(huì)形成高壓區(qū),導(dǎo)致甲板結(jié)構(gòu)所受壓力出現(xiàn)一系列壓力峰值,這樣的壓力峰值對(duì)甲板結(jié)構(gòu)具有一定的破壞性。通過引入真實(shí)氣體模型,可以快速預(yù)報(bào)此峰值載荷,對(duì)工程上船舶設(shè)計(jì)進(jìn)行載荷計(jì)算具有一定的實(shí)用價(jià)值。

    5)不同的來波條件會(huì)出現(xiàn)不同的上浪特性和上浪響應(yīng),因此還應(yīng)該進(jìn)一步研究各種來波條件下上浪時(shí)卷入空氣的效應(yīng)。

    [1] FALTINSEN O M, GRECO M, LANDRINI M. Green water loading on FPSO[C]. OMAE. Brail,2001.

    [2] GRECO, M., COLICCHIO, G., FALTINSEN, O. Shipping of water on a two- dimensionalstructure. Part 2[J]. Journal of fluid mechanics.2007,581:371- 399.

    [3] 林兆偉,朱仁傳,繆國平.甲板上浪問題的二維數(shù)值模擬[J].船舶力學(xué),2009,13(1):1- 8.

    [4] YE Haixuan, WAN Decheng. RANS simulations and validations of twodimensional green water on deck[C]. the 3th National CFD Conference on Naval Architecture and Ocean Engineering, Dalian, China, Jul. 2014:191- 197.

    [5] 張林濤.新型連接橋底封板抗砰擊性能研究[D].大連:大連理工大學(xué),2015.

    [6] 陳震,肖熙.空氣墊在平底結(jié)構(gòu)入水砰擊中作用的仿真分析[J].上海交通大學(xué)學(xué)報(bào),2005,39(5):670- 674.

    [7] 曹正林,吳衛(wèi)國.空氣層對(duì)高速三體船連接橋砰擊壓力峰值影響二維仿真研究[J].船舶力學(xué),2008,12(2):237- 242.

    [8] ?STMAN A, SILEO C P L, STANSBERG C T, et al. A fully nonlinear RANS- VOF numerical wave tank applied in the analysis of green water on FPSO in waves[C]. OMAE. 2014:V08BT06A010- V08BT06A010.

    猜你喜歡
    壓縮空氣甲板峰值
    “四單”聯(lián)動(dòng)打造適齡兒童隊(duì)前教育峰值體驗(yàn)
    科考船木甲板安裝工藝
    殲-15艦載戰(zhàn)斗機(jī)在遼寧艦飛行甲板準(zhǔn)備起飛
    國防(2017年12期)2017-04-09 07:52:09
    替代甲板結(jié)合負(fù)壓吸引治療修復(fù)甲床缺損
    300MW循環(huán)流化床機(jī)組壓縮空氣節(jié)能改造及分析
    寬占空比峰值電流型準(zhǔn)PWM/PFM混合控制
    基于峰值反饋的電流型PFM控制方法
    甲板前置治療甲床部分缺損
    發(fā)酵用壓縮空氣余熱再生吸附干燥工藝
    制藥行業(yè)壓縮空氣的凈化系統(tǒng)及驗(yàn)證
    亚洲伊人色综图| 可以免费在线观看a视频的电影网站 | 精品国产国语对白av| 国产乱来视频区| 波多野结衣av一区二区av| 伦理电影大哥的女人| 欧美亚洲 丝袜 人妻 在线| 亚洲精品自拍成人| 午夜久久久在线观看| 一边亲一边摸免费视频| 人妻人人澡人人爽人人| 亚洲专区中文字幕在线 | 日本爱情动作片www.在线观看| 国产精品久久久久久精品电影小说| 久久久国产欧美日韩av| 日本黄色日本黄色录像| 操出白浆在线播放| 好男人视频免费观看在线| 亚洲国产成人一精品久久久| 国产精品二区激情视频| 久久99精品国语久久久| 欧美另类一区| 国产精品一二三区在线看| 人妻 亚洲 视频| 久久99一区二区三区| 国产精品亚洲av一区麻豆 | 国产精品 国内视频| 中文精品一卡2卡3卡4更新| 一级毛片我不卡| 精品亚洲乱码少妇综合久久| 你懂的网址亚洲精品在线观看| 亚洲国产av新网站| 高清欧美精品videossex| 中文字幕制服av| 亚洲精品美女久久久久99蜜臀 | 亚洲情色 制服丝袜| 一二三四在线观看免费中文在| 精品视频人人做人人爽| 麻豆av在线久日| 一级片'在线观看视频| 超碰成人久久| 久久精品熟女亚洲av麻豆精品| 久久精品国产亚洲av高清一级| 最近最新中文字幕免费大全7| 亚洲欧美一区二区三区久久| 女人高潮潮喷娇喘18禁视频| 大话2 男鬼变身卡| 国产亚洲最大av| 最近的中文字幕免费完整| 午夜免费观看性视频| 欧美精品一区二区大全| 成人影院久久| 国产无遮挡羞羞视频在线观看| 午夜福利,免费看| 青草久久国产| 高清不卡的av网站| 中文字幕制服av| 精品酒店卫生间| 亚洲欧洲国产日韩| 免费在线观看视频国产中文字幕亚洲 | 大片免费播放器 马上看| 成年av动漫网址| 日韩熟女老妇一区二区性免费视频| 最黄视频免费看| 国精品久久久久久国模美| 欧美 日韩 精品 国产| 日本黄色日本黄色录像| 国产成人91sexporn| 99九九在线精品视频| 日日啪夜夜爽| 国产精品av久久久久免费| 少妇人妻久久综合中文| 国产精品久久久久久精品电影小说| 久久精品熟女亚洲av麻豆精品| 美女大奶头黄色视频| 国产精品 欧美亚洲| 久热爱精品视频在线9| 热99久久久久精品小说推荐| 亚洲精品日韩在线中文字幕| 亚洲精品av麻豆狂野| 亚洲专区中文字幕在线 | 黄色视频在线播放观看不卡| 国产毛片在线视频| 一区二区av电影网| 2021少妇久久久久久久久久久| 波多野结衣av一区二区av| 亚洲七黄色美女视频| 日韩av免费高清视频| 侵犯人妻中文字幕一二三四区| 另类亚洲欧美激情| 日韩精品有码人妻一区| 国产老妇伦熟女老妇高清| 日韩电影二区| 香蕉国产在线看| 亚洲 欧美一区二区三区| 中文字幕制服av| 波野结衣二区三区在线| 亚洲国产中文字幕在线视频| 久久精品久久久久久久性| 久久久久精品人妻al黑| 免费在线观看黄色视频的| 色播在线永久视频| 日本欧美视频一区| 制服诱惑二区| 国产av国产精品国产| 亚洲av在线观看美女高潮| 一级毛片黄色毛片免费观看视频| 国产片内射在线| 亚洲成人一二三区av| 九草在线视频观看| 精品午夜福利在线看| 黄色视频在线播放观看不卡| 看免费成人av毛片| 亚洲一级一片aⅴ在线观看| 亚洲精品日本国产第一区| 80岁老熟妇乱子伦牲交| 欧美精品一区二区免费开放| av电影中文网址| 日韩 欧美 亚洲 中文字幕| 亚洲一区二区三区欧美精品| xxxhd国产人妻xxx| 韩国精品一区二区三区| a级毛片在线看网站| 久久热在线av| 国产日韩欧美亚洲二区| 丝袜人妻中文字幕| 99九九在线精品视频| 天天影视国产精品| 丝袜人妻中文字幕| 欧美日韩国产mv在线观看视频| 美女午夜性视频免费| 日韩av在线免费看完整版不卡| 久久热在线av| 国产精品麻豆人妻色哟哟久久| 欧美精品一区二区大全| 黄色怎么调成土黄色| 丝瓜视频免费看黄片| 国产成人精品久久久久久| www.熟女人妻精品国产| 嫩草影视91久久| 亚洲七黄色美女视频| 国产视频首页在线观看| 国产乱人偷精品视频| 18禁国产床啪视频网站| 1024视频免费在线观看| 国产男女内射视频| 亚洲一区中文字幕在线| 亚洲国产日韩一区二区| 亚洲精品乱久久久久久| 亚洲精品国产区一区二| 男女下面插进去视频免费观看| 搡老乐熟女国产| 三上悠亚av全集在线观看| 99精国产麻豆久久婷婷| 亚洲情色 制服丝袜| 青春草国产在线视频| 亚洲中文av在线| 亚洲精品乱久久久久久| 欧美 日韩 精品 国产| 香蕉丝袜av| 黑人猛操日本美女一级片| 国产精品久久久久久久久免| 国产精品久久久久久人妻精品电影 | 国产色婷婷99| 免费在线观看黄色视频的| av不卡在线播放| kizo精华| 18禁裸乳无遮挡动漫免费视频| 最近的中文字幕免费完整| 国产成人精品福利久久| 18禁裸乳无遮挡动漫免费视频| 侵犯人妻中文字幕一二三四区| 岛国毛片在线播放| 黑人猛操日本美女一级片| 亚洲精品国产av蜜桃| 肉色欧美久久久久久久蜜桃| 国产视频首页在线观看| 日日啪夜夜爽| 亚洲天堂av无毛| 婷婷色综合www| 国产成人午夜福利电影在线观看| 国产97色在线日韩免费| 18禁动态无遮挡网站| 黑丝袜美女国产一区| 黑人猛操日本美女一级片| 纯流量卡能插随身wifi吗| 天堂中文最新版在线下载| 丰满迷人的少妇在线观看| 国产深夜福利视频在线观看| 只有这里有精品99| 99久国产av精品国产电影| 深夜精品福利| 国产亚洲午夜精品一区二区久久| 亚洲少妇的诱惑av| 免费观看a级毛片全部| 午夜福利影视在线免费观看| 日韩,欧美,国产一区二区三区| 色吧在线观看| 欧美 日韩 精品 国产| 国产免费福利视频在线观看| 国产免费视频播放在线视频| 国产又色又爽无遮挡免| 大码成人一级视频| 日韩制服丝袜自拍偷拍| 国产黄色免费在线视频| 99久国产av精品国产电影| 99re6热这里在线精品视频| 亚洲一级一片aⅴ在线观看| 久久久精品区二区三区| 在线看a的网站| 久久久久久久精品精品| 国产精品久久久人人做人人爽| 在线观看免费午夜福利视频| 亚洲欧美精品综合一区二区三区| 中文字幕人妻丝袜制服| 欧美日韩视频高清一区二区三区二| 精品亚洲成a人片在线观看| 汤姆久久久久久久影院中文字幕| 热99国产精品久久久久久7| 亚洲精品美女久久av网站| 精品久久蜜臀av无| 老司机靠b影院| 高清黄色对白视频在线免费看| 久久久久久久精品精品| kizo精华| 国产精品免费大片| 日本91视频免费播放| 亚洲人成电影观看| 五月天丁香电影| 看免费成人av毛片| 丰满迷人的少妇在线观看| 国产精品麻豆人妻色哟哟久久| 亚洲熟女精品中文字幕| 日韩一区二区三区影片| 午夜福利网站1000一区二区三区| av国产久精品久网站免费入址| 熟女av电影| 成年女人毛片免费观看观看9 | 免费少妇av软件| 欧美人与性动交α欧美精品济南到| 一边摸一边抽搐一进一出视频| 日韩 亚洲 欧美在线| 深夜精品福利| 男女国产视频网站| 精品视频人人做人人爽| 99久久精品国产亚洲精品| 日韩一区二区三区影片| 亚洲国产欧美一区二区综合| 高清欧美精品videossex| 老汉色av国产亚洲站长工具| 亚洲欧洲精品一区二区精品久久久 | 欧美激情 高清一区二区三区| 少妇的丰满在线观看| 黄色 视频免费看| 老鸭窝网址在线观看| 男男h啪啪无遮挡| 久久亚洲国产成人精品v| 777米奇影视久久| 免费黄网站久久成人精品| 这个男人来自地球电影免费观看 | 亚洲精品日韩在线中文字幕| 久久久久精品性色| 免费高清在线观看日韩| √禁漫天堂资源中文www| 欧美日韩精品网址| 久久精品aⅴ一区二区三区四区| 国产精品一国产av| 亚洲五月色婷婷综合| 精品国产露脸久久av麻豆| 国产成人av激情在线播放| 青春草视频在线免费观看| 精品一区二区三区四区五区乱码 | 亚洲av成人不卡在线观看播放网 | 天天躁夜夜躁狠狠久久av| 一级毛片我不卡| 国产欧美亚洲国产| 天堂8中文在线网| 精品人妻在线不人妻| 免费不卡黄色视频| 亚洲图色成人| 人人妻人人澡人人看| 性色av一级| 国产精品久久久久久人妻精品电影 | 国产毛片在线视频| 国产精品一区二区在线不卡| 考比视频在线观看| www.精华液| 免费黄频网站在线观看国产| 少妇人妻久久综合中文| 黑人巨大精品欧美一区二区蜜桃| 婷婷色综合大香蕉| 欧美 日韩 精品 国产| 十八禁人妻一区二区| a级毛片黄视频| 国产熟女欧美一区二区| 国语对白做爰xxxⅹ性视频网站| 看免费成人av毛片| 久久97久久精品| 亚洲国产精品999| 国产成人啪精品午夜网站| 日韩制服丝袜自拍偷拍| bbb黄色大片| 婷婷色av中文字幕| 欧美日韩视频精品一区| 亚洲精品中文字幕在线视频| 制服诱惑二区| 婷婷成人精品国产| 亚洲美女黄色视频免费看| av不卡在线播放| 亚洲人成网站在线观看播放| 国产精品人妻久久久影院| 欧美 日韩 精品 国产| 午夜福利在线免费观看网站| 观看美女的网站| 18禁国产床啪视频网站| 成人漫画全彩无遮挡| av视频免费观看在线观看| 777米奇影视久久| 国产亚洲av片在线观看秒播厂| 亚洲欧美成人综合另类久久久| 亚洲成人一二三区av| 这个男人来自地球电影免费观看 | 美女主播在线视频| 黄色 视频免费看| 国产精品成人在线| 建设人人有责人人尽责人人享有的| 日韩人妻精品一区2区三区| 在现免费观看毛片| 国产女主播在线喷水免费视频网站| 一本大道久久a久久精品| 亚洲av电影在线进入| 国精品久久久久久国模美| 99国产精品免费福利视频| 久久午夜综合久久蜜桃| 1024香蕉在线观看| 中文字幕最新亚洲高清| 国产一卡二卡三卡精品 | 国产亚洲av片在线观看秒播厂| 国产精品秋霞免费鲁丝片| 亚洲专区中文字幕在线 | 久久国产精品男人的天堂亚洲| 国产97色在线日韩免费| 一级爰片在线观看| 一边亲一边摸免费视频| 久久精品aⅴ一区二区三区四区| av网站在线播放免费| 在线观看一区二区三区激情| av网站在线播放免费| 国产又爽黄色视频| 午夜激情久久久久久久| 一区福利在线观看| 亚洲国产中文字幕在线视频| 国产午夜精品一二区理论片| 男女国产视频网站| 国产毛片在线视频| 久久精品久久精品一区二区三区| 国产一区二区激情短视频 | 制服人妻中文乱码| 尾随美女入室| 免费观看性生交大片5| 老汉色av国产亚洲站长工具| 日日摸夜夜添夜夜爱| 看免费成人av毛片| 亚洲婷婷狠狠爱综合网| 中文字幕另类日韩欧美亚洲嫩草| 五月天丁香电影| 男女免费视频国产| 国产精品99久久99久久久不卡 | 18禁动态无遮挡网站| 美女主播在线视频| 国产在线一区二区三区精| 免费在线观看黄色视频的| 999久久久国产精品视频| 国产极品粉嫩免费观看在线| 亚洲少妇的诱惑av| 精品午夜福利在线看| 国产色婷婷99| 黑人猛操日本美女一级片| 欧美变态另类bdsm刘玥| 9色porny在线观看| 岛国毛片在线播放| 亚洲欧美一区二区三区久久| 亚洲欧美精品自产自拍| 最近2019中文字幕mv第一页| 大香蕉久久成人网| 欧美 日韩 精品 国产| 久久久精品国产亚洲av高清涩受| 亚洲专区中文字幕在线 | 日韩精品免费视频一区二区三区| 精品一区在线观看国产| 91精品三级在线观看| 看免费av毛片| 精品少妇黑人巨大在线播放| 中文字幕av电影在线播放| 大香蕉久久网| 午夜福利影视在线免费观看| 纵有疾风起免费观看全集完整版| 一级a爱视频在线免费观看| 1024视频免费在线观看| 超碰97精品在线观看| 成人黄色视频免费在线看| 亚洲国产欧美日韩在线播放| 国产一区二区在线观看av| 男女午夜视频在线观看| 高清av免费在线| 赤兔流量卡办理| 精品第一国产精品| 五月开心婷婷网| 午夜91福利影院| 中文字幕另类日韩欧美亚洲嫩草| 9热在线视频观看99| 亚洲一区二区三区欧美精品| 美女视频免费永久观看网站| 精品久久久久久电影网| 巨乳人妻的诱惑在线观看| 午夜久久久在线观看| 久久久欧美国产精品| 久久热在线av| 日韩熟女老妇一区二区性免费视频| 亚洲欧美一区二区三区久久| 纵有疾风起免费观看全集完整版| 欧美成人午夜精品| 国产熟女午夜一区二区三区| 久久精品国产a三级三级三级| 中国三级夫妇交换| 久热这里只有精品99| 国产日韩一区二区三区精品不卡| 国产在线免费精品| 我要看黄色一级片免费的| 国产成人精品久久二区二区91 | 亚洲七黄色美女视频| 水蜜桃什么品种好| 国产一区二区三区av在线| 三上悠亚av全集在线观看| 99热网站在线观看| 欧美黑人精品巨大| 综合色丁香网| 成人亚洲欧美一区二区av| 欧美国产精品va在线观看不卡| 人成视频在线观看免费观看| 国产日韩欧美视频二区| 别揉我奶头~嗯~啊~动态视频 | 晚上一个人看的免费电影| 亚洲av成人精品一二三区| 少妇精品久久久久久久| 亚洲精品国产色婷婷电影| 亚洲成av片中文字幕在线观看| 国产成人a∨麻豆精品| 中文字幕人妻丝袜制服| 亚洲欧美日韩另类电影网站| av免费观看日本| 黑人巨大精品欧美一区二区蜜桃| 啦啦啦在线观看免费高清www| 日韩视频在线欧美| 精品一区在线观看国产| 国产一区二区三区综合在线观看| 午夜91福利影院| 午夜福利网站1000一区二区三区| 久久久精品区二区三区| 精品一区二区免费观看| 日韩一区二区视频免费看| 天堂8中文在线网| 亚洲av电影在线进入| 亚洲精品日本国产第一区| 亚洲人成电影观看| 99久久人妻综合| 一级片'在线观看视频| 欧美激情 高清一区二区三区| 日日摸夜夜添夜夜爱| 免费女性裸体啪啪无遮挡网站| 国产黄色视频一区二区在线观看| 国产精品一国产av| 免费黄频网站在线观看国产| 中文字幕最新亚洲高清| 亚洲成色77777| 五月天丁香电影| 91国产中文字幕| 国产在线一区二区三区精| 熟女少妇亚洲综合色aaa.| 久久亚洲国产成人精品v| 国产成人a∨麻豆精品| 男女免费视频国产| 亚洲精品国产区一区二| 成人国语在线视频| 宅男免费午夜| netflix在线观看网站| av网站免费在线观看视频| 国产一级毛片在线| 久久久久视频综合| 黄片小视频在线播放| 大片免费播放器 马上看| 欧美精品一区二区大全| 女人被躁到高潮嗷嗷叫费观| xxx大片免费视频| 亚洲第一青青草原| 丰满少妇做爰视频| 亚洲精品久久久久久婷婷小说| 久久人人97超碰香蕉20202| 午夜激情久久久久久久| 老汉色av国产亚洲站长工具| 我要看黄色一级片免费的| 十八禁人妻一区二区| 亚洲图色成人| 激情五月婷婷亚洲| 久久女婷五月综合色啪小说| 九色亚洲精品在线播放| 国产国语露脸激情在线看| 欧美在线黄色| 国产99久久九九免费精品| 成年美女黄网站色视频大全免费| 精品国产超薄肉色丝袜足j| 日本vs欧美在线观看视频| 考比视频在线观看| 大片免费播放器 马上看| 欧美精品人与动牲交sv欧美| 中文字幕制服av| 成人三级做爰电影| 在线天堂中文资源库| 在现免费观看毛片| 日本午夜av视频| 国产又爽黄色视频| 久久青草综合色| 亚洲国产精品一区二区三区在线| 美女视频免费永久观看网站| 在线观看三级黄色| 成年人免费黄色播放视频| 老司机影院成人| 哪个播放器可以免费观看大片| 搡老岳熟女国产| 天美传媒精品一区二区| 精品免费久久久久久久清纯 | 国产探花极品一区二区| 色婷婷久久久亚洲欧美| av福利片在线| 久久婷婷青草| 丝袜美足系列| 丝袜人妻中文字幕| 丁香六月天网| 久久久久久久久免费视频了| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久免费高清国产稀缺| 一本色道久久久久久精品综合| 欧美黄色片欧美黄色片| av不卡在线播放| 波多野结衣一区麻豆| 天天影视国产精品| 啦啦啦在线免费观看视频4| 午夜福利网站1000一区二区三区| 亚洲精品国产av蜜桃| √禁漫天堂资源中文www| 成人亚洲精品一区在线观看| 97精品久久久久久久久久精品| 日韩一卡2卡3卡4卡2021年| 宅男免费午夜| 亚洲自偷自拍图片 自拍| 亚洲色图 男人天堂 中文字幕| 久久国产精品男人的天堂亚洲| 久久精品久久精品一区二区三区| 精品福利永久在线观看| 亚洲精品国产av蜜桃| 丰满乱子伦码专区| 色网站视频免费| 韩国高清视频一区二区三区| av又黄又爽大尺度在线免费看| 成年人免费黄色播放视频| 国产精品99久久99久久久不卡 | 国产av一区二区精品久久| 好男人视频免费观看在线| 免费高清在线观看视频在线观看| 久久久国产欧美日韩av| 最近中文字幕高清免费大全6| 欧美精品一区二区大全| 亚洲精品一区蜜桃| 成人漫画全彩无遮挡| 国产免费福利视频在线观看| 狠狠精品人妻久久久久久综合| 激情五月婷婷亚洲| 免费观看人在逋| 女人被躁到高潮嗷嗷叫费观| 多毛熟女@视频| 99国产综合亚洲精品| 热re99久久精品国产66热6| 午夜激情久久久久久久| 90打野战视频偷拍视频| 韩国高清视频一区二区三区| 色综合欧美亚洲国产小说| 看十八女毛片水多多多| 青春草国产在线视频| 亚洲av福利一区| 国产免费现黄频在线看| 国产乱来视频区| 啦啦啦在线观看免费高清www| videosex国产| 久久久精品区二区三区| 日韩av免费高清视频| 久久国产精品大桥未久av| 在线观看www视频免费| 亚洲精品国产一区二区精华液| 90打野战视频偷拍视频| 精品国产超薄肉色丝袜足j| 2018国产大陆天天弄谢| 久久精品人人爽人人爽视色| 韩国高清视频一区二区三区| 高清av免费在线| 国产av精品麻豆| 美女视频免费永久观看网站| 欧美少妇被猛烈插入视频| 精品国产乱码久久久久久男人| 亚洲精品国产区一区二| svipshipincom国产片| 香蕉丝袜av| 青春草国产在线视频| 菩萨蛮人人尽说江南好唐韦庄| 久久精品熟女亚洲av麻豆精品| 色播在线永久视频|