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

    基于SHIPFLOW的散貨船阻力預(yù)報分析

    2018-10-09 06:08:18杜敏何珍田中文伍蓉暉
    廣東造船 2018年4期

    杜敏 何珍 田中文 伍蓉暉

    摘 要:本文利用SHIPFLOW軟件對某中低速靈便型散貨船進(jìn)行船模尺度下的阻力預(yù)報分析,提出SHIPFLOW適用于中低速肥大型方尾散貨船阻力預(yù)報的計(jì)算模式及改善尾部網(wǎng)格質(zhì)量、提高計(jì)算精度的方法。將阻力的CFD模擬結(jié)果與船模試驗(yàn)結(jié)果進(jìn)行對比,結(jié)果表明CFD阻力預(yù)報誤差基本控制在±3%范圍以內(nèi),滿足工程精度要求,驗(yàn)證了SHIPFLOW對散貨船阻力預(yù)報的適用性,為后續(xù)方案優(yōu)選論證提供保障。

    關(guān)鍵詞:中低速肥大型散貨船;方尾;SHIPFLOW;阻力預(yù)報;計(jì)算精度

    中圖分類號:U661.32 文獻(xiàn)標(biāo)識碼:A

    Abstract: In this paper, the resistance predictions are carried out at model scale of the medium and low speed handy-size bulk carrier by utilizing the CFD software SHIPFLOW. The appropriate methods are implemented to numerically simulate resistance coefficients of the full block bulk carrier with transom stern at relative low speed and the ways of improving the stern grid quality in terms of transom stern flow to ensure the accuracy of the numerical results are proposed. By comparison with the experimental data, the deviation could be controlled below 3%, satisfying the engineering requirement. Therefore, the feasibility of the numerical prediction method based on SHIPFLOW was validated, which provides a guarantee for the stage of design optimization and demonstration to ensure its results.

    Key words: Medium and low speed full block bulk carrier; Transom stern; SHIPFLOW; Resistance prediction; Computational accuracy

    1 引言

    目前,隨著計(jì)算流體動力學(xué)(CFD)相關(guān)理論的不斷完善和計(jì)算機(jī)技術(shù)的不斷升級,船舶CFD得到突飛猛進(jìn)的發(fā)展。其中,將船舶專用流體計(jì)算軟件SHIPFLOW與船型參數(shù)化軟件CAESES相結(jié)合進(jìn)行多方案阻力預(yù)報,已成為船型設(shè)計(jì)及前期優(yōu)化階段廣泛采用的作法。大量研究表明,SHIPFLOW的勢流/粘流求解器能夠快速有效地預(yù)估船體阻力等性能。但目前的研究大都集中于方形系數(shù)CB小于0.8的中高速細(xì)長船型方面,對于SHIPFLOW在中低速肥大船型阻力預(yù)報方面的適用性研究較少。

    本文利用SHIPFLOW軟件對某中低速靈便型散貨船進(jìn)行船模尺度下多速度點(diǎn)的阻力預(yù)報分析,研究計(jì)算方法的適用性,并通過提升尾部網(wǎng)格質(zhì)量的方法來改善計(jì)算精度,保障阻力預(yù)報的準(zhǔn)確性。

    2 計(jì)算原理

    利用SHIPFLOW進(jìn)行船舶阻力預(yù)報,按RANS粘流計(jì)算區(qū)域的劃分進(jìn)行分類,主要包括兩種計(jì)算模式:軟件默認(rèn)的分區(qū)域ZONAL算法(僅尾部基于RANS計(jì)算粘流)和全船采用RANS計(jì)算粘流的GLOBAL算法。

    如圖1所示:(1)ZONAL分區(qū)域算法,將流場劃分為勢流區(qū)、邊界層區(qū)和粘流區(qū)三個區(qū)域。針對不同的阻力成分采用對應(yīng)的理論模型、計(jì)算原理和網(wǎng)格劃分方式,從而大大縮短了模擬周期。其中,區(qū)域一為外圍勢流區(qū),采用非線性勢流理論計(jì)算興波阻力,船體和自由液面采用面元網(wǎng)格,對應(yīng)軟件的XPAN模塊;區(qū)域二為船首至船舯的邊界層流場,基于邊界層理論計(jì)算船體的摩擦阻力,對應(yīng)軟件的XBOUND模塊;區(qū)域三為船舯至船尾的粘流區(qū)域,采用結(jié)構(gòu)化體網(wǎng)格,通過有限體積法結(jié)合k ? ω SST或EASM湍流模型求解RANS方程,計(jì)算船體的摩擦阻力和粘壓阻力等,對應(yīng)XCHAP模塊;(2)GLOBAL算法,在全船范圍內(nèi)進(jìn)行粘流計(jì)算,求解RANS方程(如圖2所示)。該算法對計(jì)算設(shè)備的要求較高、運(yùn)算周期較長,但某種程度上可以提高計(jì)算精度。

    3 數(shù)值計(jì)算與分析

    3.1 計(jì)算模型

    本文選用某一方形系數(shù)CB=0.87的靈便型散貨船作為計(jì)算模型,分別在設(shè)計(jì)吃水和結(jié)構(gòu)吃水狀態(tài)下進(jìn)行船模尺度下多速度點(diǎn)(Fn=0.137~0.186)的阻力預(yù)報分析。

    3.2 網(wǎng)格劃分設(shè)置

    SHIPFLOW可以自動劃分網(wǎng)格,根據(jù)勢流、粘流計(jì)算的特點(diǎn),采用兩種網(wǎng)格生成器:XMESH面元網(wǎng)格生成器和XGRID結(jié)構(gòu)化體網(wǎng)格生成器.對于網(wǎng)格的尺寸,軟件提供了三種參考模式,便于用戶快速設(shè)置網(wǎng)格。目前,有關(guān)SHIPFLOW網(wǎng)格劃分對計(jì)算精度的影響已有較多研究成果,網(wǎng)格質(zhì)量可以滿足計(jì)算精度要求,因此本文后續(xù)計(jì)算均采用軟件自動生成的fine網(wǎng)格。

    如圖3所示,對于方尾船型的浸水尾封板,在較低速時尾封板浸沒高度較大,此時對應(yīng)的尾流為靜水流;隨著航速增大,當(dāng)相對傅汝德數(shù)Fntr>3時,尾部的水流具有足夠的動能克服粘性影響,從而使自由面脫離尾封板下邊界,此時為干尾封板。Fntr表達(dá)式為:

    通過計(jì)算不同航速對應(yīng)的,發(fā)現(xiàn)僅在設(shè)計(jì)吃水下的相對高速點(diǎn)Fn=0.186時需要設(shè)置干尾封板。計(jì)算發(fā)現(xiàn)設(shè)置Transom后,收斂問題變得復(fù)雜,需要進(jìn)行一系列收斂性測試才能得到較好的結(jié)果。

    3.3 數(shù)值分析

    SHIPFLOW利用XPAN模塊進(jìn)行勢流計(jì)算得到興波阻力,其表達(dá)方式有兩種:壓力積分法和橫切波法。其中,橫切波法對網(wǎng)格質(zhì)量的敏感性相對較低[2],因此本文采用橫切波法表達(dá)興波阻力??偟乃枇ο禂?shù)Ctwtwc等于興波阻力系數(shù)Ctwtwc、粘壓阻力系數(shù)CPV及摩擦阻力系數(shù)Cf之和:

    3.3.1 計(jì)算方法分析

    本文分別選用ZONAL法和GLOBAL法對設(shè)計(jì)吃水下對應(yīng)的多個速度點(diǎn)進(jìn)行阻力預(yù)報,得到總的水阻力系數(shù)及其精度的變化曲線,如圖4所示。其中,精度表達(dá)式為:

    由圖5分析可知:

    (1)兩種算法預(yù)報的阻力系數(shù)曲線與船模試驗(yàn)值的變化趨勢基本一致;

    (2)GLOBAL算法的計(jì)算值略高于試驗(yàn)值,誤差A(yù)Ctwtwc小于2.5%,滿足工程精度要求;ZONAL算法的計(jì)算值遠(yuǎn)低于試驗(yàn)值,誤差超過5%,不滿足工程精度要求;

    (3)相較于ZONAL算法僅對船體后半部進(jìn)行RANS粘流計(jì)算,GLOBAL算法對全船進(jìn)行RANS粘流計(jì)算,粘性阻力預(yù)報更精確,計(jì)算值更接近試驗(yàn)值;

    (4)大量研究表明:對于中高速細(xì)長船型的阻力預(yù)報,SHIPFLOW常用的ZONAL算法能夠滿足精度要求;但對于散貨船等中低速肥大船型,采用ZONAL算法預(yù)報阻力已無法滿足工程精度要求,因此在對阻力預(yù)報精度要求較高時,推薦采用能夠充分考慮全船RANS粘流計(jì)算的GLOBAL算法。

    3.3.2 尾部網(wǎng)格改善設(shè)置

    根據(jù)公式(4):隨著船速增加,興波阻力將急劇增大。因此,對于低速散貨船,由于興波阻力所占比例很小,主要為粘性阻力,故想通過調(diào)整網(wǎng)格等方式來改變摩擦阻力的精度意義不大;粘性阻力計(jì)算精度受浸水尾封板周圍的網(wǎng)格質(zhì)量影響較大。對于低速方尾船浸水尾封板的網(wǎng)格處理,SHIPFLOW默認(rèn)的網(wǎng)格劃分方式為自由面網(wǎng)格會強(qiáng)制接近尾封板的邊界,尾部網(wǎng)格質(zhì)量較差,此時可通過添加濕尾封板(簡稱Wtran)填補(bǔ)尾封板與尾流自由面之間的區(qū)域,保證網(wǎng)格自尾封板邊界至尾流自由液面之間能夠光滑過渡,生成較為真實(shí)的浸水尾封板形狀,同時加密其周圍的貼體結(jié)構(gòu)化網(wǎng)格,從而改善尾部網(wǎng)格質(zhì)量,保證在進(jìn)行粘流計(jì)算XCHAP時能夠較為精確地計(jì)算粘壓阻力。

    基于上文分析,本文后續(xù)研究均采用GLOBAL算法,并且除了Fn=0.186時需要設(shè)置干尾封板,其余速度點(diǎn)均添加濕尾封板。通過計(jì)算得到尾部網(wǎng)格改善前后各阻力系數(shù)及其精度的變化曲線,如圖5、圖6所示。

    通過分析可知:

    (1)通過添加Wtran網(wǎng)格所示,摩擦阻力系數(shù)Cf幾乎沒有變化,粘壓阻力系數(shù)CPv變化很大,減小了約7%左右,粘壓阻力計(jì)算更加合理;

    (2)添加Wtran改善尾部網(wǎng)格后,總水阻力系數(shù)的計(jì)算值仍與試驗(yàn)值的變化趨勢基本一致,且計(jì)算值有所減小,誤差大多在±1%范圍內(nèi),僅在Fn=0.174時誤差絕對值變大,整體優(yōu)化效果明顯;

    (3)當(dāng)Fn=0.174時,添加Wtran導(dǎo)致誤差絕對值變大。分析原因:當(dāng)Fn=0.186時,尾部水流已具備足夠的能量可以完全脫離尾封板,需添加Transom保證尾封板處于干燥狀態(tài);在Fn=0.174的相對高速點(diǎn),出現(xiàn)部分水質(zhì)點(diǎn)脫離尾板的現(xiàn)象,因此尾板的實(shí)際浸沒面積減小,此時若添加Wtran反而會降低計(jì)算精度。

    綜上所述,對于低速方尾船型,當(dāng)尾封板的浸沒面積較大時,可通過添加Wtran網(wǎng)格來改善尾部網(wǎng)格質(zhì)量,從而提高粘壓阻力計(jì)算精度;在相對高速點(diǎn),考慮到尾流分離現(xiàn)象,可適當(dāng)添加Transom。

    3.3.3結(jié)構(gòu)吃水下總阻力系數(shù)計(jì)算結(jié)果

    基于上文結(jié)論,采用GLOBAL算法對結(jié)構(gòu)吃水狀態(tài)下的總阻力系數(shù)進(jìn)行預(yù)報分析,所有速度點(diǎn)均設(shè)置Wtran,得到總阻力系數(shù)及精度的變化曲線,如圖7所示。

    比較Wtran添加前后的計(jì)算結(jié)果可知:在結(jié)構(gòu)吃水狀態(tài)下,添加Wtran后,總阻力系數(shù)的計(jì)算值與試驗(yàn)值的曲線變化趨勢較為一致,且計(jì)算值低于試驗(yàn)值;當(dāng)Fn≤0.161時,誤差<2%;當(dāng)Fn>0.161時,計(jì)算值逐漸變大。

    對尾部網(wǎng)格設(shè)置進(jìn)行修正,僅在Fn≤0.161的低速點(diǎn)添加Wtran,得到總阻力系數(shù)及精度變化曲線,如圖8所示:預(yù)報值更接近試驗(yàn)值,在低速點(diǎn)誤差低于2%,僅在相對高速點(diǎn)誤差較大,但誤差絕對值仍小于3%,整體預(yù)報結(jié)果較穩(wěn)定,滿足精度要求。

    4 結(jié)論

    (1)無論是設(shè)計(jì)吃水還是結(jié)構(gòu)吃水狀態(tài),采用GLOBAL算法得到的總阻力系數(shù)的誤差均在±3%范圍以內(nèi),證明利用SHIPFLOW軟件對中低速肥大型散貨船進(jìn)行阻力預(yù)報是可行的;

    (2)對于中低速方尾船型,可通過以下方法改善阻力預(yù)報精度:在低速點(diǎn)設(shè)置濕尾封板對尾部網(wǎng)格質(zhì)量進(jìn)行改善;在相對高速點(diǎn)考慮尾流分離的影響,適當(dāng)添加干尾封板;

    (3)對于SHIPFLOW提供的兩種阻力預(yù)報算法:ZONAL算法雖然精度欠佳,但模擬速度快、阻力預(yù)報值與試驗(yàn)值的變化趨勢較為一致,可用于方案優(yōu)化初期階段的多方案阻力快速預(yù)報,迅速鎖定優(yōu)化目標(biāo)范圍;GLOBAL算法雖然用時較長,但計(jì)算精度較高,可保障方案優(yōu)選論證結(jié)果的準(zhǔn)確性。二者相結(jié)合,可在保障計(jì)算精度的前提下,提升整體優(yōu)化速度。

    參考文獻(xiàn)

    [1]盧雨,胡安康,尹遜濱等. 基于 SHIPFLOW 的某大型集裝箱船阻力預(yù) 報與試驗(yàn)驗(yàn)證[J].中國造船, 2016,57(2).

    [2]陳偉,許輝, 邱遼原等.基于SHIPFLOW軟件的方尾艦船阻力快速預(yù) 報[J].中國艦船研究, 2012, 7(4).

    [3]陳紅梅,蔡榮泉. 勢流計(jì)算在船舶型線優(yōu)化改型中的適用性研究[J]. 船舶工程, 2012, 34(增2).

    [4] Mierlo K J. Trend Validation of SHIPFLOW Based on the Bare Hull Upright Resistance of the Delft Series [D]. Netherlands: Delft University of Technology, 2006.

    [5]J He, H Chen, H Yu, etc. Resistance Optimization of a Cruise Ship Using a Hybrid Approach [J]. The 27th International Ocean and Polar Engineering Conference, 2017.

    [6] Prasanta K. Sahoo, Lawrence J Doctors, etc. CFD Prediction of the Wave Resistance of a Catamaran [J]. MAHY 2006:International Conference on Marine Hydrodynamics,2006.

    [7]方昭昭,趙丙乾,陳慶任等.基于CFD的雙尾船型阻力數(shù)值預(yù)報[J]. 中國艦船研究,2014, 9(4).

    [8] FLOWTECH. SHIPFLOW6.3 Users Manual[M]. Gothenburg, Sweden, 2017: 8-9,46.

    [9] Xchap Theoretical Manual [M]. Gothenburg, Sweden, 2007.

    [10]H.C. Raven. A Solution Method for the Nonlinear Ship Wave Resistance Problem [D]. Netherlands:Delft University of Technology, 1996.

    [11] SHIPFLOW Design Tutorials Basic [M]. Gothenburg, Sweden, 2017.

    亚洲色图 男人天堂 中文字幕| 热99久久久久精品小说推荐| 亚洲五月色婷婷综合| 午夜日韩欧美国产| 毛片一级片免费看久久久久| 国产免费视频播放在线视频| 久久热在线av| 免费久久久久久久精品成人欧美视频| 亚洲综合色网址| 日韩大片免费观看网站| 99精品久久久久人妻精品| 欧美精品一区二区免费开放| av在线观看视频网站免费| 丝袜人妻中文字幕| 亚洲av欧美aⅴ国产| 久久久久久久久久久久大奶| 亚洲av日韩在线播放| xxxhd国产人妻xxx| 中国三级夫妇交换| 亚洲色图 男人天堂 中文字幕| 国产精品国产三级专区第一集| 亚洲av中文av极速乱| 成年美女黄网站色视频大全免费| 国产成人精品久久久久久| 亚洲七黄色美女视频| 日韩av不卡免费在线播放| 中文天堂在线官网| 男女免费视频国产| 久久人妻熟女aⅴ| 岛国毛片在线播放| 伊人久久国产一区二区| 黑人巨大精品欧美一区二区蜜桃| 嫩草影视91久久| 乱人伦中国视频| 波多野结衣av一区二区av| 蜜桃在线观看..| 免费在线观看完整版高清| 男女床上黄色一级片免费看| 99国产精品免费福利视频| 亚洲精品美女久久久久99蜜臀 | 无限看片的www在线观看| 热99国产精品久久久久久7| 不卡av一区二区三区| 麻豆精品久久久久久蜜桃| 麻豆乱淫一区二区| 国产无遮挡羞羞视频在线观看| 久久久久久久国产电影| 免费黄网站久久成人精品| 美女中出高潮动态图| 国产伦理片在线播放av一区| 成人亚洲精品一区在线观看| 天天操日日干夜夜撸| 免费高清在线观看视频在线观看| 国产日韩欧美视频二区| 在线 av 中文字幕| 久久久久久免费高清国产稀缺| 三上悠亚av全集在线观看| 亚洲成人免费av在线播放| 巨乳人妻的诱惑在线观看| 男女无遮挡免费网站观看| 国产又色又爽无遮挡免| 九草在线视频观看| 99精国产麻豆久久婷婷| 午夜福利网站1000一区二区三区| 亚洲激情五月婷婷啪啪| 日本av免费视频播放| 日韩不卡一区二区三区视频在线| 最近最新中文字幕免费大全7| 又粗又硬又长又爽又黄的视频| 国产片内射在线| av国产精品久久久久影院| 国产免费视频播放在线视频| 无限看片的www在线观看| 少妇 在线观看| 18禁动态无遮挡网站| 欧美精品高潮呻吟av久久| h视频一区二区三区| 国产精品久久久人人做人人爽| 一区二区日韩欧美中文字幕| 国产精品国产av在线观看| 高清视频免费观看一区二区| 国产男女超爽视频在线观看| 午夜福利网站1000一区二区三区| 在线观看免费日韩欧美大片| 一区二区av电影网| 女人精品久久久久毛片| 高清av免费在线| svipshipincom国产片| 高清不卡的av网站| 久久久久精品性色| 国产精品久久久久久精品古装| 熟妇人妻不卡中文字幕| 国产一区有黄有色的免费视频| 男女无遮挡免费网站观看| 在线精品无人区一区二区三| av有码第一页| 国产成人精品无人区| 91国产中文字幕| 亚洲美女视频黄频| 成人三级做爰电影| svipshipincom国产片| 一二三四中文在线观看免费高清| 一级毛片我不卡| 狂野欧美激情性bbbbbb| 成年av动漫网址| 免费少妇av软件| 亚洲成人国产一区在线观看 | 伦理电影大哥的女人| 男女免费视频国产| 欧美精品一区二区免费开放| 看非洲黑人一级黄片| 日韩av免费高清视频| 久久久久久久精品精品| 精品卡一卡二卡四卡免费| 一本久久精品| www日本在线高清视频| 国产精品国产三级国产专区5o| 亚洲精品一区蜜桃| 少妇被粗大的猛进出69影院| 精品人妻熟女毛片av久久网站| 亚洲成色77777| 1024视频免费在线观看| 99香蕉大伊视频| 可以免费在线观看a视频的电影网站 | 91国产中文字幕| 99热国产这里只有精品6| 国产男女内射视频| 卡戴珊不雅视频在线播放| 午夜激情av网站| 亚洲国产成人一精品久久久| 日韩 欧美 亚洲 中文字幕| 亚洲,欧美,日韩| 中文字幕另类日韩欧美亚洲嫩草| 婷婷色麻豆天堂久久| 午夜老司机福利片| 男女边吃奶边做爰视频| 一二三四中文在线观看免费高清| 2021少妇久久久久久久久久久| 天堂俺去俺来也www色官网| 丁香六月天网| 日韩精品免费视频一区二区三区| 黄片小视频在线播放| 国产av码专区亚洲av| 国产 一区精品| 777米奇影视久久| 国产人伦9x9x在线观看| 天天躁日日躁夜夜躁夜夜| 国产精品国产三级国产专区5o| 人人妻,人人澡人人爽秒播 | 色视频在线一区二区三区| 国产在视频线精品| 在线观看www视频免费| 一个人免费看片子| 老汉色av国产亚洲站长工具| 满18在线观看网站| 亚洲精品,欧美精品| 777米奇影视久久| 久久久久精品国产欧美久久久 | 老司机在亚洲福利影院| 99热全是精品| 亚洲国产毛片av蜜桃av| 99国产精品免费福利视频| 亚洲七黄色美女视频| 精品国产一区二区久久| 美女扒开内裤让男人捅视频| 欧美乱码精品一区二区三区| 97人妻天天添夜夜摸| videos熟女内射| 日本黄色日本黄色录像| 天天操日日干夜夜撸| 亚洲美女视频黄频| 丰满迷人的少妇在线观看| 亚洲 欧美一区二区三区| 97人妻天天添夜夜摸| 极品少妇高潮喷水抽搐| 欧美黑人欧美精品刺激| 欧美日韩一区二区视频在线观看视频在线| 大香蕉久久成人网| 一本色道久久久久久精品综合| 中国国产av一级| 无限看片的www在线观看| 嫩草影视91久久| 午夜福利,免费看| 日韩av在线免费看完整版不卡| 精品久久久久久电影网| 黄色怎么调成土黄色| 丝袜脚勾引网站| 精品人妻一区二区三区麻豆| 国产在线视频一区二区| 一级,二级,三级黄色视频| 亚洲一级一片aⅴ在线观看| 男男h啪啪无遮挡| 亚洲精品国产区一区二| 秋霞伦理黄片| 狠狠婷婷综合久久久久久88av| 操出白浆在线播放| 国产一区二区 视频在线| 亚洲精品日本国产第一区| 999久久久国产精品视频| 久久婷婷青草| 久久人妻熟女aⅴ| 日韩一卡2卡3卡4卡2021年| 最近中文字幕高清免费大全6| 欧美日韩精品网址| 在线观看一区二区三区激情| 成人免费观看视频高清| 男人操女人黄网站| 国产视频首页在线观看| 午夜激情av网站| 美女午夜性视频免费| 欧美日韩亚洲综合一区二区三区_| 69精品国产乱码久久久| 久久午夜综合久久蜜桃| 久久久久精品性色| 搡老岳熟女国产| 国产男女超爽视频在线观看| av在线观看视频网站免费| 久久久久精品国产欧美久久久 | avwww免费| 99国产综合亚洲精品| 日韩精品有码人妻一区| 最近中文字幕高清免费大全6| 婷婷成人精品国产| 国产不卡av网站在线观看| 午夜福利网站1000一区二区三区| 久久久久视频综合| 亚洲色图综合在线观看| 亚洲av综合色区一区| 国产午夜精品一二区理论片| 亚洲精品aⅴ在线观看| 国产精品亚洲av一区麻豆 | www.av在线官网国产| 日本午夜av视频| 午夜免费鲁丝| 一级a爱视频在线免费观看| 欧美在线黄色| 在线看a的网站| 大码成人一级视频| 亚洲欧美清纯卡通| 韩国精品一区二区三区| 美女高潮到喷水免费观看| 日本午夜av视频| 人人澡人人妻人| 一二三四在线观看免费中文在| 日韩 欧美 亚洲 中文字幕| av.在线天堂| 一本色道久久久久久精品综合| 色播在线永久视频| 人人妻,人人澡人人爽秒播 | 国产乱人偷精品视频| 黄色视频不卡| 欧美日韩视频精品一区| 日韩av免费高清视频| 99久国产av精品国产电影| 国产日韩欧美亚洲二区| 在线天堂中文资源库| 男女之事视频高清在线观看 | 精品国产露脸久久av麻豆| 国产 一区精品| 日韩成人av中文字幕在线观看| 五月天丁香电影| 国产xxxxx性猛交| videosex国产| 日本午夜av视频| 久久午夜综合久久蜜桃| 精品少妇一区二区三区视频日本电影 | 亚洲av在线观看美女高潮| 欧美激情极品国产一区二区三区| 国产一区二区三区av在线| 9色porny在线观看| 亚洲成人一二三区av| 黄色视频不卡| 老汉色∧v一级毛片| 十八禁高潮呻吟视频| 夫妻性生交免费视频一级片| 久久这里只有精品19| 在线观看免费高清a一片| 狠狠精品人妻久久久久久综合| 国产探花极品一区二区| 如日韩欧美国产精品一区二区三区| 青春草亚洲视频在线观看| 最近2019中文字幕mv第一页| 午夜免费观看性视频| 人成视频在线观看免费观看| 国产精品免费大片| 在线 av 中文字幕| 搡老岳熟女国产| 中文字幕人妻丝袜制服| 国产熟女午夜一区二区三区| 亚洲情色 制服丝袜| 激情视频va一区二区三区| 久久天堂一区二区三区四区| 久久久久久久久久久久大奶| 在线观看三级黄色| 女性被躁到高潮视频| 在线观看免费视频网站a站| 在线看a的网站| 亚洲av福利一区| 免费在线观看视频国产中文字幕亚洲 | www.av在线官网国产| 少妇被粗大猛烈的视频| 久久韩国三级中文字幕| 久久人人97超碰香蕉20202| 国产精品一区二区在线观看99| 日韩精品有码人妻一区| 天天躁日日躁夜夜躁夜夜| 欧美在线黄色| 男女边摸边吃奶| 国产精品 欧美亚洲| 亚洲三区欧美一区| 黄片小视频在线播放| 大话2 男鬼变身卡| 亚洲av福利一区| 日韩 亚洲 欧美在线| 久久毛片免费看一区二区三区| 国产精品偷伦视频观看了| av在线播放精品| 极品少妇高潮喷水抽搐| 51午夜福利影视在线观看| 欧美中文综合在线视频| 另类精品久久| 中文字幕人妻丝袜制服| 亚洲精品国产一区二区精华液| 狂野欧美激情性bbbbbb| 亚洲欧美成人精品一区二区| 亚洲精品,欧美精品| 啦啦啦在线免费观看视频4| 狠狠精品人妻久久久久久综合| 黄色一级大片看看| 最近最新中文字幕免费大全7| 国产精品久久久久久精品电影小说| 亚洲视频免费观看视频| 99香蕉大伊视频| 亚洲欧美清纯卡通| avwww免费| 国产 一区精品| 欧美乱码精品一区二区三区| 国产日韩欧美在线精品| 国产成人精品久久二区二区91 | 国产精品久久久久久精品电影小说| 亚洲视频免费观看视频| 观看av在线不卡| 人妻人人澡人人爽人人| 一区福利在线观看| 两性夫妻黄色片| 一级毛片我不卡| 考比视频在线观看| 国产精品一国产av| 日韩欧美精品免费久久| 精品人妻在线不人妻| 日韩大码丰满熟妇| 久久国产亚洲av麻豆专区| www.自偷自拍.com| 纵有疾风起免费观看全集完整版| 丰满乱子伦码专区| 妹子高潮喷水视频| 国产成人精品在线电影| 成人国产av品久久久| 日韩大片免费观看网站| 亚洲精品国产一区二区精华液| 超碰成人久久| kizo精华| 哪个播放器可以免费观看大片| 午夜日韩欧美国产| 日本欧美视频一区| 亚洲av欧美aⅴ国产| 成人影院久久| 中文字幕制服av| 日韩熟女老妇一区二区性免费视频| 成人国产av品久久久| 中文欧美无线码| 嫩草影院入口| 天天影视国产精品| 亚洲av欧美aⅴ国产| 久久影院123| 熟妇人妻不卡中文字幕| 亚洲欧美一区二区三区黑人| 啦啦啦在线观看免费高清www| 亚洲欧美一区二区三区国产| 亚洲视频免费观看视频| av免费观看日本| 久久精品国产综合久久久| 一级爰片在线观看| 亚洲精品日韩在线中文字幕| 看免费av毛片| 久久久久久久国产电影| 亚洲av电影在线进入| 少妇的丰满在线观看| 97精品久久久久久久久久精品| a级片在线免费高清观看视频| 如日韩欧美国产精品一区二区三区| 日日摸夜夜添夜夜爱| 一个人免费看片子| 51午夜福利影视在线观看| 大片电影免费在线观看免费| 狂野欧美激情性bbbbbb| 亚洲,欧美,日韩| 亚洲成国产人片在线观看| av不卡在线播放| 久久精品熟女亚洲av麻豆精品| 中文欧美无线码| 久久韩国三级中文字幕| 日本欧美视频一区| av在线观看视频网站免费| 亚洲国产最新在线播放| 巨乳人妻的诱惑在线观看| 国产黄色免费在线视频| 亚洲成人手机| 亚洲视频免费观看视频| 日日撸夜夜添| 国产精品一二三区在线看| 男女之事视频高清在线观看 | 亚洲精品日韩在线中文字幕| 嫩草影院入口| 欧美变态另类bdsm刘玥| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品99久久99久久久不卡 | 久久国产亚洲av麻豆专区| 美女午夜性视频免费| 久久精品久久久久久噜噜老黄| 久久久久精品性色| 90打野战视频偷拍视频| 久久影院123| 亚洲欧美一区二区三区久久| 国产日韩欧美亚洲二区| 国产极品天堂在线| 好男人视频免费观看在线| 精品人妻在线不人妻| 国产成人免费观看mmmm| 久久久久久免费高清国产稀缺| 99国产综合亚洲精品| 最近中文字幕高清免费大全6| 午夜免费观看性视频| 亚洲国产av影院在线观看| 亚洲精品成人av观看孕妇| 欧美激情极品国产一区二区三区| 熟女少妇亚洲综合色aaa.| 久久人人97超碰香蕉20202| 一区二区三区精品91| 高清欧美精品videossex| 亚洲久久久国产精品| 一级片'在线观看视频| 亚洲精品国产av蜜桃| 国产精品国产三级专区第一集| 一区二区三区乱码不卡18| 日韩,欧美,国产一区二区三区| 久久影院123| 久久精品aⅴ一区二区三区四区| 国产一区二区在线观看av| 亚洲精品中文字幕在线视频| 天天操日日干夜夜撸| 午夜91福利影院| 亚洲成av片中文字幕在线观看| 亚洲综合色网址| 丝瓜视频免费看黄片| 国产黄频视频在线观看| 国产老妇伦熟女老妇高清| 欧美xxⅹ黑人| 精品人妻熟女毛片av久久网站| 一边亲一边摸免费视频| 久久人人97超碰香蕉20202| 亚洲av中文av极速乱| 热re99久久国产66热| 亚洲,欧美精品.| av又黄又爽大尺度在线免费看| 亚洲成人一二三区av| 日韩av不卡免费在线播放| 国产精品久久久久成人av| 国产精品.久久久| 亚洲国产精品一区三区| 久久狼人影院| 少妇人妻久久综合中文| 国产亚洲一区二区精品| a级毛片在线看网站| 好男人视频免费观看在线| 亚洲在久久综合| 青春草亚洲视频在线观看| 日韩欧美一区视频在线观看| 七月丁香在线播放| 国产爽快片一区二区三区| 91国产中文字幕| 最近2019中文字幕mv第一页| 成年美女黄网站色视频大全免费| 精品少妇内射三级| 黄色一级大片看看| 综合色丁香网| 久久亚洲国产成人精品v| 国产亚洲欧美精品永久| 丁香六月天网| 亚洲欧美成人综合另类久久久| 蜜桃国产av成人99| 1024视频免费在线观看| 99久国产av精品国产电影| 日韩熟女老妇一区二区性免费视频| 精品一区二区三区av网在线观看 | 蜜桃国产av成人99| 日韩av不卡免费在线播放| 日韩免费高清中文字幕av| 五月开心婷婷网| 久久人人爽av亚洲精品天堂| 少妇被粗大的猛进出69影院| 国产一区亚洲一区在线观看| av在线app专区| 精品久久久久久电影网| 晚上一个人看的免费电影| 免费高清在线观看视频在线观看| 国产精品久久久久久精品电影小说| 免费女性裸体啪啪无遮挡网站| 在线观看免费高清a一片| 国产在线免费精品| 亚洲情色 制服丝袜| 欧美激情 高清一区二区三区| √禁漫天堂资源中文www| av在线老鸭窝| 亚洲国产av影院在线观看| 久久久久久人妻| 中文字幕人妻熟女乱码| 精品国产超薄肉色丝袜足j| 天天操日日干夜夜撸| 欧美 亚洲 国产 日韩一| 国产精品av久久久久免费| 亚洲国产欧美网| 老司机影院毛片| 9191精品国产免费久久| 90打野战视频偷拍视频| 啦啦啦在线免费观看视频4| 日本黄色日本黄色录像| av国产精品久久久久影院| 久久精品aⅴ一区二区三区四区| 国产麻豆69| 最近的中文字幕免费完整| 精品酒店卫生间| 久久久精品国产亚洲av高清涩受| 日韩欧美精品免费久久| 欧美av亚洲av综合av国产av | 亚洲国产欧美在线一区| 超碰成人久久| 人妻 亚洲 视频| 七月丁香在线播放| 99久久人妻综合| 日本wwww免费看| 97在线人人人人妻| 十分钟在线观看高清视频www| 丁香六月天网| 各种免费的搞黄视频| 一个人免费看片子| 国产精品亚洲av一区麻豆 | 国产成人91sexporn| 伊人久久国产一区二区| 国产成人免费观看mmmm| 丝袜喷水一区| 国产男人的电影天堂91| 一区福利在线观看| √禁漫天堂资源中文www| 久久精品aⅴ一区二区三区四区| 欧美中文综合在线视频| 久久天堂一区二区三区四区| 搡老乐熟女国产| 99国产精品免费福利视频| 看免费av毛片| 午夜福利免费观看在线| 亚洲综合精品二区| 国产免费现黄频在线看| 丰满少妇做爰视频| 免费女性裸体啪啪无遮挡网站| 欧美国产精品va在线观看不卡| 亚洲av日韩在线播放| 欧美日韩av久久| 亚洲精品久久久久久婷婷小说| 免费看av在线观看网站| 欧美人与性动交α欧美精品济南到| 老司机靠b影院| 久久这里只有精品19| 精品亚洲成国产av| 国产亚洲欧美精品永久| 国产精品国产av在线观看| 在线观看国产h片| 中文字幕人妻丝袜一区二区 | 亚洲国产日韩一区二区| 嫩草影院入口| 黄频高清免费视频| 午夜激情av网站| 久久性视频一级片| 女人被躁到高潮嗷嗷叫费观| 国产欧美日韩综合在线一区二区| 如何舔出高潮| 亚洲av男天堂| 欧美在线黄色| 中国国产av一级| 大香蕉久久网| 精品久久久精品久久久| 男女免费视频国产| 亚洲精品成人av观看孕妇| 免费观看a级毛片全部| 亚洲精品一区蜜桃| av网站在线播放免费| 国产又色又爽无遮挡免| 色视频在线一区二区三区| 国产成人精品福利久久| 亚洲国产欧美网| 亚洲av日韩在线播放| 久久人人爽人人片av| 欧美日韩一级在线毛片| 极品少妇高潮喷水抽搐| 又粗又硬又长又爽又黄的视频| 赤兔流量卡办理| 亚洲精品久久午夜乱码| 久久精品人人爽人人爽视色| 欧美人与性动交α欧美软件| 久久天躁狠狠躁夜夜2o2o | 悠悠久久av| 日本vs欧美在线观看视频| 欧美在线黄色| 久久人人爽av亚洲精品天堂| 国产成人精品在线电影|