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

    船舶阻力換算中形狀因子的計(jì)算研究

    2015-07-01 07:57:34李彥強(qiáng)
    兵器裝備工程學(xué)報(bào) 2015年4期
    關(guān)鍵詞:模型

    江 杰,李彥強(qiáng)

    (海軍裝備技術(shù)研究所動(dòng)力技術(shù)研究室,北京 102400)

    船模阻力試驗(yàn)換算方法最早由傅汝德提出,該方法將船體總阻力分成摩擦阻力Rf和剩余阻力Rr,并認(rèn)為摩擦阻力僅與雷諾數(shù)Re 有關(guān),剩余阻力僅與傅汝德數(shù)Fr 有關(guān),且兩者互不干擾。在此基礎(chǔ)上,休斯提出了三因次換算方法,引入形狀因子1 +k 來(lái)表達(dá)船體粘性阻力和相當(dāng)平板摩擦阻力之間的關(guān)系。引入形狀因子的三因次法較二因次法合理,其中形狀因子可由低速拖曳或疊模試驗(yàn)確定。目前,三因次換算方法被廣泛應(yīng)用于各大水池阻力試驗(yàn)中,但是休斯提出的1 +k 為常數(shù)的假定是否正確有待進(jìn)一步探討[1]。

    有關(guān)尺度效應(yīng)的試驗(yàn)研究有幾何相似船模組試驗(yàn),通過(guò)幾何相似船模組試驗(yàn)不僅可以比較不同換算方法的合理性,還能分析形狀因子及推進(jìn)效率中各種成分的尺度效應(yīng)等問(wèn)題。但是船模試驗(yàn)往往受到各種限制,幾何相似船模組試驗(yàn)亦不能給出實(shí)船雷諾數(shù)下的信息。隨著現(xiàn)代船舶流體力學(xué)的發(fā)展,有關(guān)船模實(shí)船阻力換算中尺度效應(yīng)的計(jì)算和驗(yàn)證方面的研究不斷增多。歐洲幾大水池合作的VIRTUE[2],EFFORT 等項(xiàng)目均包含實(shí)船數(shù)據(jù)測(cè)量及全尺度外推技術(shù)和尺度效應(yīng)相關(guān)的研究。Raven[3]等對(duì)HTC 船模進(jìn)行了模型和實(shí)尺雷諾數(shù)下的數(shù)值計(jì)算,計(jì)算結(jié)果表明實(shí)尺雷諾數(shù)下形狀因子比模型雷諾數(shù)下大7%。倪崇本[4]采用k -ε 湍流模型計(jì)算了S60 船模在傅汝德數(shù)在0.1 到0.2 各航速下形狀因子,計(jì)算結(jié)果表明形狀因子1 +k 隨航速的增加而增大。

    船模實(shí)船阻力換算常常通過(guò)低速時(shí)船模阻力數(shù)據(jù)或者低速時(shí)疊模計(jì)算來(lái)求得船體的形狀因子。但是本文計(jì)算表明雷諾數(shù)及航行姿態(tài)對(duì)某些船型的形狀因子產(chǎn)生顯著的影響,若采用正浮姿態(tài)下船模尺度的疊模計(jì)算求取形狀因子將導(dǎo)致較大的誤差。本文基于數(shù)值計(jì)算比較了航速,姿態(tài),實(shí)船雷諾數(shù)對(duì)DTMB5415 及KCS 船模形狀因子的影響。計(jì)算表明:對(duì)DTMB5415 船型,考慮航速、姿態(tài)、實(shí)船雷諾數(shù)的影響,在Fr=0.41 時(shí)實(shí)船雷諾數(shù)下的形狀因子比Fr=0.1 模型雷諾數(shù)下形狀因子大44%。對(duì)KCS 船型,考慮航速,實(shí)船雷諾數(shù)的影響,其在Fr =0.271 下實(shí)船形狀因子比其在Fr =0.086 下船模形狀因子大7%?;诘退傧炉B模計(jì)算或者拖曳試驗(yàn)得到的形狀因子必然導(dǎo)致?lián)Q算得到的實(shí)船阻力產(chǎn)生較大的誤差。

    1 數(shù)值計(jì)算基礎(chǔ)理論

    1.1 控制方程及計(jì)算模型

    采用有限體積法離散控制方程,不可壓縮粘性流體的控制方程如下:

    式(1)、(2)分別為連續(xù)性方程和動(dòng)量方程。ui分別為流體質(zhì)點(diǎn)在i 方向的速度分量,fi是質(zhì)量力;p 為流體的壓力;μ 是相體積分?jǐn)?shù)平均的動(dòng)力黏度系數(shù);μt為湍流動(dòng)力黏性系數(shù)。

    上述方程需要結(jié)合湍流模型封閉方程組,常見的兩方程湍流模型有k-ε、RNG k-ε、k-ω、SST k-ω、RSM 等。由于船體繞流帶有復(fù)雜的分離,伴流現(xiàn)象,而湍流模式對(duì)準(zhǔn)確預(yù)報(bào)伴流及流動(dòng)分離有重要的作用。本文參考文獻(xiàn)[5],選用SST k-ω 湍流模型,基本方程參考文獻(xiàn)[6]。

    SST k-ω 湍流模型集合了k -ε 模型和k -ω 模型的優(yōu)點(diǎn):在近壁面區(qū)域,采用k -ω 模型,湍流耗散較小,收斂性好;在湍流充分發(fā)展的區(qū)域,采用k -ε 模型數(shù)值模擬,計(jì)算效率高。SST k-ω 湍流模型對(duì)壁面和尾流場(chǎng)模擬效果好,對(duì)自由來(lái)流的湍流度不敏感[7],具有較好的穩(wěn)定性。

    1.2 船舶阻力換算方法

    船模阻力試驗(yàn)換算方法最早由傅汝德提出,該方法將船體總阻力分成摩擦阻力Rf 和剩余阻力Rr,并認(rèn)為摩擦阻力僅與雷諾數(shù)Rn 有關(guān),剩余阻力僅與傅汝德數(shù)Fn 有關(guān),且兩者互不干擾[1]。因此總阻力系數(shù)可以表示為

    因此實(shí)船尺度下的阻力系數(shù)可由船模阻力試驗(yàn)中測(cè)得的阻力系數(shù)表示為

    其中Cf0代表相當(dāng)平板摩擦阻力系數(shù),Cr 為剩余阻力系數(shù),下標(biāo)s 和m 分別代表變量對(duì)應(yīng)于實(shí)船和模型。在此基礎(chǔ)上,休斯提出了三因次換算方法,引入形狀因子1 +k 來(lái)表達(dá)船體粘性阻力和相當(dāng)平板摩擦阻力之間的關(guān)系:

    此時(shí)有

    引入形狀因子的三因次法較二因次法合理,其中形狀因子可由低速拖曳或疊模試驗(yàn)確定??紤]到實(shí)船船體表面粗糙度的作用以及船模實(shí)船阻力換算過(guò)程中由于雷諾數(shù)不相等導(dǎo)致的尺度效應(yīng),需要補(bǔ)貼一定的經(jīng)驗(yàn)系數(shù),即船模實(shí)船換算補(bǔ)貼Ca,Ca可由式(7)確定:

    Lpp為垂線間長(zhǎng),ks 為粗糙度表觀高度,一般可取ks=0.15mm。

    2 數(shù)值計(jì)算及分析

    2.1 DTMB5415 船型數(shù)值計(jì)算

    DTMB5415 是ITTC 推薦的軍艦類型中唯一有大量公開試驗(yàn)數(shù)據(jù)的研究對(duì)象(圖1)。美國(guó)的DTMB,IIHR 以及意大利的INSEAN 三家研究機(jī)構(gòu)對(duì)該船型進(jìn)行了全面的試驗(yàn)。DTMB 5415 船??s尺比為λ=24.8,其實(shí)船與船模的主要數(shù)據(jù)如表1 所示。

    表1 DTMB5415 船體幾何尺寸

    圖1 DTMB5415 疊模模型幾何,表面網(wǎng)格及網(wǎng)格質(zhì)量

    在實(shí)際船模試驗(yàn)時(shí),通常忽略低航速下船舶的興波阻力,采用傅汝德數(shù)低于0.15 的船模試驗(yàn)數(shù)據(jù)來(lái)計(jì)算船模的形狀因子。但是低航速下流動(dòng)不穩(wěn)定,導(dǎo)致測(cè)量阻力也不穩(wěn)定。本文采用疊模計(jì)算可以在較大的范圍內(nèi)求取船舶的粘壓阻力及形狀因子。疊模模型尺寸與相關(guān)模型試驗(yàn)一致,以船模中縱剖面與靜水面為對(duì)稱面。數(shù)值計(jì)算區(qū)域計(jì)算區(qū)域入口取船艏向前延伸1 倍船長(zhǎng)處,出口取船艉向后延伸2 倍船長(zhǎng)處,側(cè)邊界及下方邊界均取1 倍船長(zhǎng)。

    2.1.1 模型雷諾數(shù)下形狀因子隨航速變化

    鑒于船模試驗(yàn)通常采用傅汝德數(shù)低于0.15 的船模阻力數(shù)據(jù)來(lái)計(jì)算形狀因子,本文首先計(jì)算了傅汝德數(shù)0.1 ~0.15各航速下的粘壓阻力系數(shù)及形狀因子。然后計(jì)算了傅汝德數(shù)在0.28,0.36,0.41 下船模的形狀因子,以此分析粘壓阻力系數(shù)及形狀因子在較大航速范圍內(nèi)變化規(guī)律。為表示粘壓阻力系數(shù)與摩擦阻力系數(shù)變化的關(guān)系,定義k1為粘壓阻力與摩擦阻力的比值。表2 列出了Fr 從0.1 ~0.41 航速范圍內(nèi)粘壓阻力系數(shù),形狀因子1 +k 及k1的變化。

    表2 不同航速下形狀因子的變化

    從表1 可以看出,隨著航速增加(即雷諾數(shù)增大),粘壓阻力系數(shù)下降,但是粘壓阻力系數(shù)相對(duì)于摩擦阻力系數(shù)不斷增大,即k1和形狀因子增大。當(dāng)Fr 從0.1 增加到0.41 時(shí),k1增大了13%。在較低的航速范圍內(nèi),形狀因子變化不大,F(xiàn)r 從0.1 ~0.15 航速范圍內(nèi),形狀因子增大了0.7%。但是Fr 由0.1 增大到0.41,1 +k 值增大了5%。因此若采用低航速下計(jì)算得到的形狀因子來(lái)外推實(shí)船在高航速時(shí)的阻力,導(dǎo)致粘性阻力低估5%。

    2.1.2 模型雷諾數(shù)下形狀因子隨姿態(tài)的變化

    當(dāng)船舶航速較大時(shí),船體的姿態(tài)對(duì)阻力產(chǎn)生顯著的影響,此時(shí)主要體現(xiàn)在濕表面積及各項(xiàng)阻力系數(shù)的變化。考慮船體姿態(tài)變化,針對(duì)不同航速下姿態(tài)分別進(jìn)行建模與網(wǎng)格劃分,此時(shí)疊模模型為該航速下首尾吃水的連線以下的船體部分。鑒于低航速下船模姿態(tài)變化較小,對(duì)阻力的影響可以忽略不計(jì),本文僅計(jì)算Fr 大于0.28 的3 個(gè)航速。不同航速下船體姿態(tài)如表3,不同姿態(tài)下形狀因子的變化如表4。

    表3 不同航速下船模的航行姿態(tài)

    表4 不同姿態(tài)下形狀因子的變化

    計(jì)算表明,考慮船體吃水和縱傾變化之后,濕表面積增大,摩擦阻力增大,但是由于雷諾數(shù)不變,摩擦阻力系數(shù)不變。在較高航速下船模產(chǎn)生較大的平行下沉及縱傾,粘壓阻力系數(shù)及形狀因子k 增大,3 個(gè)航速下形狀因子分別增大了1%,2%,12%。Fr=0.28 和0.36 時(shí),船體姿態(tài)變化不大,粘壓阻力系數(shù)變化不明顯;當(dāng)Fr =0.41 時(shí),船體姿態(tài)變化較大,粘壓阻力系數(shù)增大了117%。

    2.1.3 實(shí)船雷諾數(shù)下形狀因子隨航速變化

    不考慮船體姿態(tài),實(shí)船雷諾數(shù)下的粘壓阻力系數(shù)及形狀因子變化如表5。

    表5 不考慮姿態(tài)實(shí)船雷諾數(shù)下形狀因子的變化

    在實(shí)船雷諾數(shù)下,隨著航速的增大,粘壓阻力系數(shù)略微減小,摩擦阻力系數(shù)不斷降低,故形狀因子增大。實(shí)尺雷諾數(shù)的形狀因子要比模型雷諾數(shù)下形狀因子分別大5.7%,6.4%,6.4%,6.2%。

    考慮船體姿態(tài),實(shí)船雷諾數(shù)下粘壓阻力系數(shù)及形狀因子在各航速變化如表6??紤]船體姿態(tài)變化和實(shí)船雷諾數(shù)的影響之后,船舶的形狀因子變化比較明顯。其中在Fr 為0.41時(shí),形狀因子較該航速模型雷諾數(shù)下不考慮姿態(tài)計(jì)算值增加了37%,較Fr=0.1 模型雷諾數(shù)下形狀因子計(jì)算值增加了44%。若按照形狀因子為一常數(shù)的假定,采用低速正浮下疊模計(jì)算得到的形狀因子來(lái)外推高航速的粘性阻力時(shí),將會(huì)導(dǎo)致粘性阻力低估44%(表6)。

    表6 考慮姿態(tài)實(shí)船雷諾數(shù)下形狀因子的變化

    2.2 KCS 船模計(jì)算

    KCS 作為國(guó)際船舶操縱性比較研究的貨船船型,其各種試驗(yàn)數(shù)據(jù)在Gothenburg 2000 會(huì)議資料中公布,表7 給出了其主要尺寸。本文對(duì)KCS 進(jìn)行疊模計(jì)算的疊模模型均未考慮各航速下船體姿態(tài)的變化(圖2)。

    表7 KCS 船體幾何尺寸

    圖2 KCS 疊模模型幾何,表面網(wǎng)格及網(wǎng)格質(zhì)量

    2.2.1 模型雷諾數(shù)下形狀因子隨航速變化

    模型雷諾數(shù)下KCS 形狀因子隨航速變化如下。

    表8 不同航速下形狀因子的變化

    從表8 可以看出,F(xiàn)r 從0.087 增大到0.271,粘壓阻力系數(shù)與摩擦阻力系數(shù)下降,形狀因子增大了3. 4%,k1增大了6%。

    2.2.2 實(shí)船雷諾數(shù)下形狀因子隨航速變化

    實(shí)船雷諾數(shù)下KCS 隨航速變化如下(表9)。

    表9 不考慮姿態(tài)實(shí)船雷諾數(shù)下形狀因子的變化

    實(shí)船雷諾數(shù)下,形狀因子隨航速變化較小,在Fr 從0.087 到0.271 的航速范圍內(nèi),形狀因子變化不超過(guò)0.5%,k1隨航速變化不超過(guò)3.5%。相比模型雷諾數(shù)下疊模計(jì)算值,形狀因子分別增大了6.3%,5.4%,3.6%,3.4%,k1的變化不超過(guò)1.3%。

    從上述計(jì)算可以看出,KCS 在Fr =0.271 時(shí)實(shí)船雷諾數(shù)下形狀因子比Fr=0.087 模型雷諾數(shù)下形狀因子大7%,若采用模型雷諾數(shù)下低航速疊模計(jì)算得到形狀因子會(huì)產(chǎn)生較大誤差。KCS 的形狀因子變化沒有5415 形狀因子變化大,主要原因可能是:其最高航速較小且未考慮航行姿態(tài)的變化,水線以下線型較為光順,流動(dòng)分離較少。

    3 結(jié)論

    1)本文通過(guò)對(duì)DTMB 5415 和KCS 的疊模計(jì)算證明了船舶的形狀因子并不是一個(gè)常數(shù),航速,航行姿態(tài)及實(shí)船雷諾數(shù)均對(duì)形狀因子產(chǎn)生顯著的影響。

    2)考慮航速、姿態(tài)、實(shí)船雷諾數(shù)的影響,DTMB 5415 在Fr=0.41 時(shí)實(shí)船雷諾數(shù)下的形狀因子比Fr =0.1 模型雷諾數(shù)下形狀因子大44%。若采用低速正浮疊模計(jì)算得到形狀因子進(jìn)行換算來(lái)求取實(shí)船高速時(shí)的阻力將導(dǎo)致明顯的誤差。

    3)對(duì)KCS 船型而言,航速及實(shí)船雷諾數(shù)的影響使形狀因子變化了7%。

    現(xiàn)有的船模實(shí)船換算方法之所以沒有產(chǎn)生明顯的誤差,主要是采用低速拖曳試驗(yàn)得到的阻力中含有一定的興波阻力成分,此時(shí)得到的形狀因子較疊模計(jì)算值偏大。且換算過(guò)程中結(jié)合船模實(shí)船換算補(bǔ)貼,使其在一定航速范圍內(nèi)具有較高的精度。為精確預(yù)報(bào)實(shí)船在較高航速的阻力,必須要考慮形狀因子的變化,此時(shí)形狀因子可以通過(guò)數(shù)值計(jì)算或者試驗(yàn)來(lái)求取。

    [1]盛振邦,劉應(yīng)中.船舶原理[M].上海:上海交通大學(xué)出版社,2003.

    [2]李勝忠,李斌,趙峰,李力楓.VIRTUE 計(jì)劃研究進(jìn)展綜述[J].船舶力學(xué).2009:13(4).

    [3]Raven H C,van der Ploeg A,Starke A R,et al.Towards a CFD-based prediction of ship performance-Progress in predicting full-scale resistance and scale effects[J]. International Journal of Maritime Engineering,Transactions RINA,2009,150(3):54-57.

    [4]倪崇本,朱仁傳,繆國(guó)平,等.一種基于CFD 的船舶總阻力預(yù)報(bào)方法[J].水動(dòng)力學(xué)研究與進(jìn)展(A 輯),2010,25(5):579-586.

    [5]張志榮,趙峰,李百奇.k -ω 湍流模式在船舶粘性流場(chǎng)計(jì)算中的應(yīng)用[J].船舶力學(xué),2003,7(1):33-37.

    [6]STERN F.Final report and recommendations of the resistance committee[C]//Proceedings of 23rd international towing tank conference. Seoul/Shanghai: ITTC,2002,173-246.

    [7]Menter F R.Two-equation Eddy-viscosity turbulence models for engineering applications[J].AIAA Journal,1994(32):1598-1605.

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    永久网站在线| 久久精品国产亚洲av天美| 色综合色国产| 国产精品一二三区在线看| 亚洲内射少妇av| 国产亚洲精品久久久com| 亚洲国产日韩欧美精品在线观看| 日韩欧美三级三区| 欧美日韩亚洲高清精品| 亚洲精品视频女| 最近手机中文字幕大全| 中文资源天堂在线| 舔av片在线| 伊人久久精品亚洲午夜| 97精品久久久久久久久久精品| 一夜夜www| 国产黄色免费在线视频| 精品国内亚洲2022精品成人| 精品人妻熟女av久视频| 欧美另类一区| 免费大片黄手机在线观看| 国产精品女同一区二区软件| 欧美潮喷喷水| 2018国产大陆天天弄谢| 十八禁国产超污无遮挡网站| 欧美日韩精品成人综合77777| 亚洲一区高清亚洲精品| 91久久精品电影网| 69av精品久久久久久| 成人高潮视频无遮挡免费网站| 18禁在线播放成人免费| 国产伦在线观看视频一区| 男女边吃奶边做爰视频| 精品久久久噜噜| 美女黄网站色视频| 真实男女啪啪啪动态图| 亚州av有码| 99热这里只有是精品50| 高清日韩中文字幕在线| 成人美女网站在线观看视频| 午夜激情久久久久久久| 亚洲欧美成人精品一区二区| 91在线精品国自产拍蜜月| 女人久久www免费人成看片| 国产精品人妻久久久久久| 国产精品精品国产色婷婷| 色5月婷婷丁香| 日日摸夜夜添夜夜爱| 国产欧美另类精品又又久久亚洲欧美| 成人高潮视频无遮挡免费网站| 乱码一卡2卡4卡精品| 伦理电影大哥的女人| av在线天堂中文字幕| 午夜激情久久久久久久| 熟妇人妻不卡中文字幕| 久久久久久久久久黄片| 18+在线观看网站| 美女高潮的动态| 日韩 亚洲 欧美在线| 日韩 亚洲 欧美在线| 91在线精品国自产拍蜜月| 午夜激情久久久久久久| 3wmmmm亚洲av在线观看| 纵有疾风起免费观看全集完整版 | 美女大奶头视频| 老师上课跳d突然被开到最大视频| 男女国产视频网站| av福利片在线观看| 亚洲人与动物交配视频| 精品久久久久久久久亚洲| 成人综合一区亚洲| 搡老妇女老女人老熟妇| 三级国产精品欧美在线观看| 秋霞在线观看毛片| 亚洲伊人久久精品综合| 午夜视频国产福利| 好男人视频免费观看在线| 在线观看av片永久免费下载| 国产v大片淫在线免费观看| 久久久久久久久久黄片| 内射极品少妇av片p| 国内精品宾馆在线| 亚洲成人精品中文字幕电影| 亚洲av成人精品一区久久| 亚洲最大成人av| 国产男女超爽视频在线观看| 非洲黑人性xxxx精品又粗又长| 哪个播放器可以免费观看大片| 又爽又黄a免费视频| 听说在线观看完整版免费高清| 免费av不卡在线播放| 久久精品夜色国产| 99久久精品热视频| 久久精品国产亚洲av涩爱| 一级黄片播放器| 大香蕉97超碰在线| 在线观看美女被高潮喷水网站| 亚洲av福利一区| 1000部很黄的大片| 直男gayav资源| 国产黄色小视频在线观看| 22中文网久久字幕| 亚洲不卡免费看| 噜噜噜噜噜久久久久久91| 美女高潮的动态| 国产黄频视频在线观看| 免费观看av网站的网址| 中文天堂在线官网| 亚洲精品,欧美精品| 国产成人91sexporn| 久久精品综合一区二区三区| 一夜夜www| 91精品伊人久久大香线蕉| 国产精品久久久久久精品电影小说 | 国产综合懂色| 午夜福利成人在线免费观看| 高清在线视频一区二区三区| av线在线观看网站| 亚洲精华国产精华液的使用体验| 日韩成人伦理影院| 老司机影院毛片| 亚洲av中文av极速乱| 国产精品一区二区性色av| 三级男女做爰猛烈吃奶摸视频| 国产伦理片在线播放av一区| 午夜久久久久精精品| 亚洲图色成人| 看黄色毛片网站| 国产成人freesex在线| 亚洲精品视频女| 国产亚洲5aaaaa淫片| 日本色播在线视频| 日韩欧美三级三区| 久久久精品94久久精品| 大片免费播放器 马上看| 亚洲欧美中文字幕日韩二区| 婷婷色综合www| 国产单亲对白刺激| 少妇的逼水好多| 免费看光身美女| 国产精品久久久久久久电影| 久久久午夜欧美精品| 国产精品人妻久久久影院| 在线观看美女被高潮喷水网站| 亚洲国产精品成人久久小说| 99久久人妻综合| 欧美丝袜亚洲另类| 老司机影院成人| 成人鲁丝片一二三区免费| 免费少妇av软件| 十八禁网站网址无遮挡 | 国产 一区精品| 综合色av麻豆| av在线播放精品| 男女啪啪激烈高潮av片| 精品国产三级普通话版| 日本av手机在线免费观看| 日日摸夜夜添夜夜爱| 又爽又黄无遮挡网站| 亚洲欧美成人精品一区二区| 男女啪啪激烈高潮av片| 一级片'在线观看视频| 91久久精品国产一区二区成人| 熟女人妻精品中文字幕| 欧美高清性xxxxhd video| 欧美激情久久久久久爽电影| 久久国产乱子免费精品| 少妇裸体淫交视频免费看高清| 嫩草影院入口| 国产精品蜜桃在线观看| kizo精华| 黄片无遮挡物在线观看| 欧美性猛交╳xxx乱大交人| 欧美丝袜亚洲另类| 国产成人一区二区在线| 国产不卡一卡二| 啦啦啦中文免费视频观看日本| 如何舔出高潮| 十八禁国产超污无遮挡网站| 亚洲第一区二区三区不卡| 我的女老师完整版在线观看| 少妇猛男粗大的猛烈进出视频 | 精品欧美国产一区二区三| 又黄又爽又刺激的免费视频.| 午夜免费男女啪啪视频观看| 色综合色国产| 在线观看一区二区三区| 18禁在线无遮挡免费观看视频| 天堂影院成人在线观看| 色网站视频免费| 日韩欧美精品v在线| 亚洲av一区综合| 久99久视频精品免费| 国产成人午夜福利电影在线观看| 男女边摸边吃奶| 在线观看一区二区三区| 卡戴珊不雅视频在线播放| 亚洲精品中文字幕在线视频 | 国产色婷婷99| 天天躁夜夜躁狠狠久久av| 97在线视频观看| 久久久久性生活片| 亚洲,欧美,日韩| 精品久久久久久久久亚洲| 日日干狠狠操夜夜爽| 99久国产av精品| 国产伦理片在线播放av一区| 男女边摸边吃奶| 日韩精品有码人妻一区| 99九九线精品视频在线观看视频| 久99久视频精品免费| 三级男女做爰猛烈吃奶摸视频| 国产极品天堂在线| 少妇裸体淫交视频免费看高清| 免费看日本二区| 大香蕉久久网| 少妇人妻一区二区三区视频| 欧美一区二区亚洲| 韩国av在线不卡| 一级毛片aaaaaa免费看小| 久久久久久久久久久免费av| 热99在线观看视频| 好男人视频免费观看在线| 91精品国产九色| 国产午夜精品论理片| 中文欧美无线码| 啦啦啦中文免费视频观看日本| 国产乱人偷精品视频| av天堂中文字幕网| 黑人高潮一二区| 美女cb高潮喷水在线观看| 精品人妻一区二区三区麻豆| 女人久久www免费人成看片| 日韩伦理黄色片| 国产黄色小视频在线观看| 三级国产精品片| a级毛片免费高清观看在线播放| 激情 狠狠 欧美| 韩国高清视频一区二区三区| 中文字幕av在线有码专区| 男的添女的下面高潮视频| 日本猛色少妇xxxxx猛交久久| 亚洲人与动物交配视频| 91精品一卡2卡3卡4卡| 国产乱人视频| 日韩av免费高清视频| 少妇人妻精品综合一区二区| 亚洲人成网站在线观看播放| 日本与韩国留学比较| 偷拍熟女少妇极品色| 久久精品国产亚洲av天美| 免费av不卡在线播放| 亚洲精品乱久久久久久| av女优亚洲男人天堂| 久久精品国产鲁丝片午夜精品| 精品不卡国产一区二区三区| 国产精品av视频在线免费观看| 国产黄色小视频在线观看| 久久久久性生活片| 国产伦一二天堂av在线观看| 久久热精品热| 亚洲国产欧美人成| 国语对白做爰xxxⅹ性视频网站| 99热这里只有精品一区| 十八禁网站网址无遮挡 | 午夜免费男女啪啪视频观看| 亚洲精品日韩在线中文字幕| 夫妻午夜视频| 亚洲精品影视一区二区三区av| 超碰97精品在线观看| 久久久欧美国产精品| 91午夜精品亚洲一区二区三区| 亚洲欧美日韩东京热| 天堂网av新在线| 两个人的视频大全免费| 嫩草影院精品99| 韩国高清视频一区二区三区| 国产精品一区二区在线观看99 | 久久韩国三级中文字幕| 久久久久久久久久久丰满| 亚洲欧美一区二区三区黑人 | 国产男人的电影天堂91| 亚洲精品国产av成人精品| 成人综合一区亚洲| 日日啪夜夜撸| 一级毛片久久久久久久久女| 少妇熟女欧美另类| 国产午夜精品一二区理论片| 男人舔奶头视频| 最近视频中文字幕2019在线8| 国产人妻一区二区三区在| 色5月婷婷丁香| 在线观看一区二区三区| 久久午夜福利片| 亚洲欧美成人精品一区二区| 日本熟妇午夜| 五月伊人婷婷丁香| 一区二区三区四区激情视频| 婷婷色av中文字幕| 国产一区二区三区综合在线观看 | 淫秽高清视频在线观看| a级毛色黄片| 人妻一区二区av| 亚洲成人精品中文字幕电影| 国产高清有码在线观看视频| 亚洲人成网站在线播| 99视频精品全部免费 在线| 爱豆传媒免费全集在线观看| 国产精品久久久久久精品电影小说 | 少妇人妻一区二区三区视频| 嫩草影院入口| 亚洲av成人精品一区久久| 三级经典国产精品| 黑人高潮一二区| 亚洲国产精品成人久久小说| 久久久久久久久久久免费av| 免费不卡的大黄色大毛片视频在线观看 | 亚洲一级一片aⅴ在线观看| 插阴视频在线观看视频| 观看免费一级毛片| 日韩在线高清观看一区二区三区| 亚洲精品自拍成人| 最近中文字幕高清免费大全6| 久久草成人影院| 国产精品一区二区在线观看99 | 女人十人毛片免费观看3o分钟| a级毛色黄片| 一个人看视频在线观看www免费| 美女大奶头视频| 亚洲高清免费不卡视频| 免费黄网站久久成人精品| 少妇熟女aⅴ在线视频| 亚洲欧洲国产日韩| 女人十人毛片免费观看3o分钟| 欧美极品一区二区三区四区| 三级经典国产精品| 久久人人爽人人爽人人片va| 最近2019中文字幕mv第一页| av在线亚洲专区| 亚洲综合色惰| 成年人午夜在线观看视频 | 国产精品一区二区性色av| 国产精品精品国产色婷婷| 狠狠精品人妻久久久久久综合| 亚洲成色77777| 激情五月婷婷亚洲| 熟妇人妻不卡中文字幕| 国内精品美女久久久久久| 国产高清国产精品国产三级 | 九九在线视频观看精品| 亚洲国产欧美人成| 国产69精品久久久久777片| 亚洲天堂国产精品一区在线| 高清午夜精品一区二区三区| 精品午夜福利在线看| 国产老妇伦熟女老妇高清| .国产精品久久| 97在线视频观看| 午夜久久久久精精品| 深夜a级毛片| 午夜精品在线福利| 亚洲无线观看免费| 婷婷六月久久综合丁香| 毛片女人毛片| 国产午夜精品论理片| 亚州av有码| 亚洲av电影不卡..在线观看| 国产一级毛片七仙女欲春2| 搡老乐熟女国产| 老女人水多毛片| 亚洲自拍偷在线| 麻豆成人午夜福利视频| 亚洲色图av天堂| 亚洲av电影在线观看一区二区三区 | 日本免费a在线| 精品久久久噜噜| 国产亚洲5aaaaa淫片| 最后的刺客免费高清国语| 色尼玛亚洲综合影院| 成人亚洲欧美一区二区av| av在线天堂中文字幕| 精品一区二区三卡| 免费av毛片视频| 成人二区视频| 婷婷色麻豆天堂久久| 国产片特级美女逼逼视频| 岛国毛片在线播放| 又大又黄又爽视频免费| 亚洲成人久久爱视频| 熟女电影av网| 国产 亚洲一区二区三区 | 亚洲国产色片| 久久久欧美国产精品| 黄色日韩在线| 看非洲黑人一级黄片| 午夜免费激情av| 久久鲁丝午夜福利片| 男女边吃奶边做爰视频| 18禁动态无遮挡网站| 亚洲av国产av综合av卡| 美女内射精品一级片tv| 日韩成人伦理影院| 天堂影院成人在线观看| 久久久久久九九精品二区国产| 日韩强制内射视频| 又黄又爽又刺激的免费视频.| 亚洲国产精品sss在线观看| 午夜精品在线福利| 国产成人精品福利久久| 91久久精品国产一区二区三区| 亚洲av在线观看美女高潮| 国产精品福利在线免费观看| 中文乱码字字幕精品一区二区三区 | 日产精品乱码卡一卡2卡三| 禁无遮挡网站| 国产成人精品福利久久| 国产精品爽爽va在线观看网站| 真实男女啪啪啪动态图| 黄片wwwwww| 日韩一区二区三区影片| 我的老师免费观看完整版| 在线观看免费高清a一片| 日韩在线高清观看一区二区三区| 又粗又硬又长又爽又黄的视频| 亚洲精品一区蜜桃| 色综合亚洲欧美另类图片| 天堂影院成人在线观看| 精品久久久久久久久av| 国产伦在线观看视频一区| 免费观看性生交大片5| 一个人看的www免费观看视频| 国语对白做爰xxxⅹ性视频网站| 国产免费福利视频在线观看| 秋霞在线观看毛片| 亚洲精品日韩av片在线观看| 伊人久久国产一区二区| 国产av不卡久久| 中文字幕av在线有码专区| 深爱激情五月婷婷| 国产一区二区三区综合在线观看 | 日韩强制内射视频| 国内精品一区二区在线观看| 久久久久九九精品影院| 日韩人妻高清精品专区| 欧美丝袜亚洲另类| 日本午夜av视频| 丝瓜视频免费看黄片| 日日撸夜夜添| 国产精品国产三级专区第一集| 亚洲最大成人手机在线| 纵有疾风起免费观看全集完整版 | 日日摸夜夜添夜夜添av毛片| av一本久久久久| 特级一级黄色大片| 欧美xxxx黑人xx丫x性爽| 国产一级毛片七仙女欲春2| 最近的中文字幕免费完整| 中文欧美无线码| 如何舔出高潮| 亚洲av男天堂| 免费观看精品视频网站| 久久综合国产亚洲精品| h日本视频在线播放| 日本欧美国产在线视频| freevideosex欧美| 内地一区二区视频在线| 亚洲熟妇中文字幕五十中出| 国产精品人妻久久久久久| av又黄又爽大尺度在线免费看| 欧美变态另类bdsm刘玥| 国产高清有码在线观看视频| 亚洲国产高清在线一区二区三| 国产精品久久久久久精品电影小说 | 可以在线观看毛片的网站| 国产视频首页在线观看| 免费看不卡的av| 欧美日韩视频高清一区二区三区二| 丝袜喷水一区| 日日啪夜夜撸| 九九爱精品视频在线观看| 日韩欧美三级三区| 久久久久久久久中文| 男女边吃奶边做爰视频| 内射极品少妇av片p| 亚洲精品影视一区二区三区av| 插逼视频在线观看| 国产精品久久视频播放| 不卡视频在线观看欧美| 久久精品久久精品一区二区三区| 人妻制服诱惑在线中文字幕| 欧美zozozo另类| 午夜福利在线观看免费完整高清在| 欧美xxxx性猛交bbbb| 婷婷色综合大香蕉| 精品久久久噜噜| 99九九线精品视频在线观看视频| 国产男女超爽视频在线观看| 亚洲欧美一区二区三区国产| 日韩伦理黄色片| 亚洲人成网站高清观看| 深爱激情五月婷婷| 国产老妇女一区| 久久精品夜色国产| 热99在线观看视频| 99久久精品一区二区三区| 久久精品久久精品一区二区三区| 午夜福利在线观看免费完整高清在| 人妻系列 视频| 亚洲国产av新网站| 欧美变态另类bdsm刘玥| 一个人观看的视频www高清免费观看| 亚洲欧美一区二区三区黑人 | 亚洲欧美精品专区久久| 国产精品久久久久久久久免| eeuss影院久久| 99热这里只有是精品在线观看| 亚洲真实伦在线观看| 国产一级毛片在线| 日韩伦理黄色片| 午夜精品国产一区二区电影 | 国产亚洲av片在线观看秒播厂 | 国产亚洲av嫩草精品影院| 最近中文字幕2019免费版| 成人无遮挡网站| 国产老妇女一区| 午夜福利视频1000在线观看| 欧美xxⅹ黑人| 亚洲av免费高清在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久久久久久电影| 欧美日韩国产mv在线观看视频 | 亚洲伊人久久精品综合| 国产乱人偷精品视频| 亚洲av成人精品一二三区| 亚洲精品456在线播放app| 精品人妻视频免费看| 中文字幕亚洲精品专区| 最新中文字幕久久久久| 精品99又大又爽又粗少妇毛片| 日韩精品有码人妻一区| 欧美成人午夜免费资源| 麻豆国产97在线/欧美| 日日摸夜夜添夜夜添av毛片| 久久久久九九精品影院| 亚洲自拍偷在线| 男女那种视频在线观看| 校园人妻丝袜中文字幕| 人妻夜夜爽99麻豆av| 晚上一个人看的免费电影| 精品午夜福利在线看| 免费看光身美女| 中文乱码字字幕精品一区二区三区 | 国产精品嫩草影院av在线观看| 91av网一区二区| 大又大粗又爽又黄少妇毛片口| 校园人妻丝袜中文字幕| 1000部很黄的大片| 99热这里只有精品一区| av一本久久久久| 最近中文字幕2019免费版| 建设人人有责人人尽责人人享有的 | 亚洲av成人精品一区久久| 天美传媒精品一区二区| 纵有疾风起免费观看全集完整版 | 国产亚洲5aaaaa淫片| 在线观看人妻少妇| 18禁在线播放成人免费| 国产免费福利视频在线观看| 成人亚洲精品av一区二区| 有码 亚洲区| .国产精品久久| 麻豆精品久久久久久蜜桃| 爱豆传媒免费全集在线观看| 国产成人freesex在线| 成人av在线播放网站| 国内揄拍国产精品人妻在线| 国产又色又爽无遮挡免| 欧美zozozo另类| 亚洲精品亚洲一区二区| 国产亚洲av片在线观看秒播厂 | 国产在线男女| 免费观看性生交大片5| 国产爱豆传媒在线观看| 亚洲美女搞黄在线观看| 国产精品精品国产色婷婷| 亚洲欧美精品自产自拍| 久久草成人影院| 日韩av免费高清视频| 欧美最新免费一区二区三区| 成年免费大片在线观看| 亚洲av福利一区| 91精品一卡2卡3卡4卡| 91久久精品国产一区二区成人| 国产午夜福利久久久久久| 狠狠精品人妻久久久久久综合| 国产视频首页在线观看| 欧美xxⅹ黑人| 熟妇人妻不卡中文字幕| 欧美性猛交╳xxx乱大交人| 大话2 男鬼变身卡| 丝瓜视频免费看黄片| 国产精品女同一区二区软件| 国产女主播在线喷水免费视频网站 | 肉色欧美久久久久久久蜜桃 | 国内精品一区二区在线观看| 午夜激情欧美在线| 久久久久久久久久黄片| 色综合色国产| 国产视频内射| 一边亲一边摸免费视频| 国产高清不卡午夜福利| 久热久热在线精品观看| 成人特级av手机在线观看| 免费观看在线日韩| 亚洲久久久久久中文字幕|