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

    2 200箱集裝箱船球鼻首優(yōu)化

    2018-07-02 10:24:52田中文白秀琴
    艦船科學(xué)技術(shù) 2018年6期
    關(guān)鍵詞:興波集裝箱船船型

    董 鵬,田中文,白秀琴,何 珍

    (1.武漢理工大學(xué) 國(guó)家水運(yùn)安全工程技術(shù)研究中心可靠性工程研究所,湖北 武漢 430063;2.武漢理工大學(xué) 船舶動(dòng)力工程技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430063;3.中船黃埔文沖船舶有限公司,廣東 廣州 510725)

    0 引 言

    海運(yùn)貿(mào)易對(duì)全球經(jīng)濟(jì)發(fā)展非常重要,目前,全球90%的貿(mào)易運(yùn)輸量由船舶航運(yùn)承擔(dān)[1]。同時(shí),海運(yùn)業(yè)會(huì)消耗大量的能源,并產(chǎn)生大量廢氣等污染物,因此保護(hù)環(huán)境、發(fā)展綠色船舶技術(shù)顯得尤為重要。船舶減阻是一種重要的船舶節(jié)能、減排方法。球鼻首作為船體的特殊部位,可以有效地改善水面波系的分布,減小船體興波阻力,對(duì)節(jié)能型船舶的發(fā)展有很重要的意義。

    傳統(tǒng)的球鼻首型線設(shè)計(jì)方法是針對(duì)具體的船型特點(diǎn),根據(jù)一些基本規(guī)律和設(shè)計(jì)經(jīng)驗(yàn),確定幾個(gè)型線方案,通過(guò)多組模型試驗(yàn)分析計(jì)算來(lái)改進(jìn)型線并最終確定,這種設(shè)計(jì)過(guò)程非常繁瑣,需要反復(fù)修改,消耗大量的人力物力,效率不高[2]。國(guó)外方面,C.Yang和H.Kim基于徑向基函數(shù)修改船體曲面的參數(shù)化方法,完成了KCS船型的興波阻力性能優(yōu)化[3],但是由于船體曲面的復(fù)雜性,故很難用一個(gè)或者多個(gè)函數(shù)將曲面形狀精確表達(dá)出來(lái),這種方法也一直沒(méi)有得到廣泛應(yīng)用。

    為了快速獲得球鼻首阻力最小最優(yōu)型線,本文利用CAESES軟件對(duì)2 200箱集裝箱船球鼻首進(jìn)行參數(shù)化建模,并集成Shipflow軟件,利用遺傳算法對(duì)球鼻首型線進(jìn)行自動(dòng)優(yōu)化,獲得具有最小阻力的該船型球鼻首型線。

    1 球鼻首參數(shù)化模型的建立

    2 200箱集裝箱船是廣州文沖船廠建造的淺吃水支線集裝箱船,適合全球航行,并滿足巴拿馬、蘇伊士運(yùn)河等航道的通航要求。該船尺度靈便,航速適中,節(jié)能環(huán)保,具有優(yōu)秀的裝載能力和油耗性能,滿足最新的國(guó)際規(guī)范規(guī)則;布置合理緊湊,貨艙利用率低,振動(dòng)噪音性能佳,技術(shù)性能已經(jīng)達(dá)到國(guó)際先進(jìn)水平行列。該船的主要參數(shù)及數(shù)據(jù)如表1所示。

    表1 2 200箱集裝箱船主要參數(shù)

    1.1 球鼻首建模

    球鼻首參數(shù)化是球鼻首進(jìn)行優(yōu)化設(shè)計(jì)的基礎(chǔ),采用CAESES軟件對(duì)球鼻首進(jìn)行參數(shù)化建模。建模思想為:首先構(gòu)建控制船體面元的特征曲線以及曲面生成器Curve Engine,然后根據(jù)這些特征曲線和曲面生成器Curve Engine通過(guò)特征Feature生成橫剖線,最后由曲面生成器Meta Surface將這些創(chuàng)建好的橫剖線平滑連接起來(lái),生成光順曲面。

    由于本文只對(duì)球鼻首部分進(jìn)行優(yōu)化研究,故只對(duì)球鼻首進(jìn)行參數(shù)化建模。球鼻首創(chuàng)建的大致過(guò)程如下:

    1)創(chuàng)建構(gòu)成球鼻首曲面的特征輪廓線,并創(chuàng)建好相應(yīng)的 Feature,然后利用 Curve Engine 和 Meta Surface功能生成球鼻首前段曲面bulb。

    2)用Fillet Surface功能實(shí)現(xiàn)前端曲面到船體的光滑過(guò)渡,并將過(guò)渡面元分為兩部分S1和S2。

    3)為了實(shí)現(xiàn)前段和上部分船體的連接,需要?jiǎng)?chuàng)建球鼻首尖端bulbtip。和球鼻首前端曲面一樣,創(chuàng)建好曲面控制線和Feature之后,啟動(dòng)Curve Engine,并利用Meta Surface生成球鼻首尖端面元bulbtip。

    4)采用Coons Patch功能連接bulbtip、S2面元以及相鄰船體2塊面元的4條邊,即可形成最后一塊面元S4,并且與周?chē)?個(gè)曲面保持相切。

    球鼻首參數(shù)化之后的模型如圖1所示。

    1.2 模型精確度分析

    對(duì)于建立的船體模型,需要對(duì)船體進(jìn)行靜水力數(shù)據(jù)計(jì)算。獲得球鼻首建模之后的船體模型的排水量、浮心位置、濕表面積等性能之后,分析和初始船型之間的誤差,從而保證球鼻首建模之后船體型線精確度[4]。

    通過(guò)在CAESES軟件中添加靜水力計(jì)算模塊進(jìn)行計(jì)算,并對(duì)數(shù)據(jù)進(jìn)行分析之后,發(fā)現(xiàn)設(shè)計(jì)吃水時(shí)的排水體積誤差僅為0.033%;浮心縱向位置誤差為0.020%;一些其他靜水力數(shù)據(jù)也符合設(shè)計(jì)要求,建模后船型的靜水力數(shù)據(jù)誤差保持在很小得到范圍內(nèi),如表2所示。

    其中誤差計(jì)算公式如下式:

    表2 設(shè)計(jì)吃水下的靜水力數(shù)據(jù)對(duì)比

    2 優(yōu)化模型的建立

    在優(yōu)化之前需要建立恰當(dāng)?shù)膬?yōu)化數(shù)學(xué)模型[5],對(duì)船型進(jìn)行優(yōu)化獲得最小阻力是典型的優(yōu)化問(wèn)題。對(duì)于一個(gè)完整的優(yōu)化問(wèn)題,應(yīng)該要包括設(shè)計(jì)變量、約束條件、優(yōu)化目標(biāo)和合適的優(yōu)化算法。

    2.1 優(yōu)化變量的選擇

    設(shè)計(jì)變量的選擇取決于要研究的問(wèn)題,選取對(duì)優(yōu)化問(wèn)題有關(guān)系、影響優(yōu)化目標(biāo)的參數(shù)作為優(yōu)化變量。球鼻首是一個(gè)復(fù)雜的曲面空間,對(duì)球鼻首線型優(yōu)化來(lái)說(shuō),需要選擇一些能夠表達(dá)和改變球鼻首線型的參數(shù)作為設(shè)計(jì)變量。在2 200箱集裝箱船的參數(shù)化球鼻首中,選擇具有代表性的特征參數(shù)作為設(shè)計(jì)變量,并根據(jù)實(shí)際情況確定設(shè)計(jì)變量的范圍。本文選取球鼻首建模時(shí)創(chuàng)建的26個(gè)特征參數(shù)作為優(yōu)化變量,并確定參數(shù)變化范圍。這些參數(shù)包括控制輪廓線上點(diǎn)的坐標(biāo)值,如長(zhǎng)度、最前點(diǎn)高度、最大寬度、最大寬度線的高度等參數(shù)和控制面元豐滿度的角度等參數(shù)。

    2.2 約束條件的定義

    對(duì)于獲得最小阻力球鼻首的問(wèn)題,其約束條件必須滿足船體的需要,在獲得最小阻力的同時(shí),要保證船體優(yōu)化前后某些參數(shù)保持不變,例如排水體積的限制、浸水面的限制、吃水匹配和寬度的限制、耐波性要求等。本文選擇無(wú)因次量排水體積V′和浮心縱向位置XCOF作為球鼻首優(yōu)化的約束條件,控制其在1%以內(nèi)變化[6]。

    其中:;;V為排水體積;LPP為首尾柱之間長(zhǎng)度;LCB為浮心縱向位置。

    2.3 優(yōu)化目標(biāo)的定義

    船型優(yōu)化是一個(gè)綜合權(quán)衡的過(guò)程,優(yōu)化過(guò)程中要考慮到船舶的各項(xiàng)性能和經(jīng)濟(jì)指標(biāo),本文的研究只限于船舶型線的阻力性能優(yōu)化,即在約束條件下以船體阻力為優(yōu)化對(duì)象進(jìn)行優(yōu)化。在優(yōu)化問(wèn)題中,優(yōu)化目標(biāo)是核心,優(yōu)化目標(biāo)一般用目標(biāo)函數(shù)來(lái)表示。根據(jù)優(yōu)化設(shè)計(jì)變量的輸入,按照某種計(jì)算規(guī)則獲得一個(gè)優(yōu)化指標(biāo),獲得目標(biāo)函數(shù)最優(yōu)解即為優(yōu)化問(wèn)題的結(jié)果。在船型優(yōu)化中,目標(biāo)函數(shù)可以是耐波性、阻力、操縱性等。對(duì)球鼻首進(jìn)行優(yōu)化研究,文章中輸出興波阻力作為優(yōu)化問(wèn)題的目標(biāo)函數(shù)。

    式中:Rwm,vm,Sm,Cwm分別為興波阻力,速度,濕表面積,興波阻力系數(shù)。

    2.4 優(yōu)化算法的選擇

    當(dāng)一個(gè)實(shí)際的問(wèn)題轉(zhuǎn)化為數(shù)學(xué)問(wèn)題時(shí),需要采用數(shù)學(xué)優(yōu)化算法來(lái)求解,在數(shù)學(xué)理論中有很多的優(yōu)化算法,大部分傳統(tǒng)算法是從單個(gè)初始值迭代求最優(yōu)解,往往誤入局部最優(yōu)解,而且難以解決復(fù)雜的非線性的問(wèn)題。遺傳算法基于生物遺傳和進(jìn)化機(jī)制,適合于復(fù)雜系統(tǒng)優(yōu)化問(wèn)題的求解,是一種高效、全局的自適應(yīng)概率優(yōu)化技術(shù),這里選取遺傳算法作為2 200箱集裝箱船球鼻首的優(yōu)化算法。

    3 基于遺傳算法的球鼻首優(yōu)化分析

    在自然界的生物進(jìn)化中,每一種物種的子代都來(lái)源于父代,但其遺傳物質(zhì)又不完全同于父代。個(gè)體的性征由染色體決定,染色體上面又是由有序排列的基因組成的,染色體上的基因?qū)Νh(huán)境有不同的適應(yīng)性,通過(guò)基因雜交和變異也可以產(chǎn)生適應(yīng)性強(qiáng)的個(gè)體被保留下來(lái)。遺傳算法(GA)[7–8]是基于自然選擇的生物進(jìn)化機(jī)制的優(yōu)化算法,把優(yōu)化問(wèn)題和自然遺傳一一對(duì)應(yīng)起來(lái),優(yōu)化解的搜索空間就是遺傳空間,將優(yōu)化問(wèn)題的每一個(gè)可能出現(xiàn)的解編碼為染色體的二進(jìn)制串(二進(jìn)制編碼只是一種編碼方法,還有其他的編碼方法,比如浮點(diǎn)編碼),每個(gè)染色體對(duì)應(yīng)一個(gè)解,染色體二進(jìn)制串的每一位稱(chēng)為基因;這里將每個(gè)染色體對(duì)應(yīng)于個(gè)體,一定數(shù)量的個(gè)體(即一定數(shù)量的染色體)組成種群。

    遺傳算法首先隨機(jī)產(chǎn)生一定數(shù)量的種群(即優(yōu)化問(wèn)題的解的集合),按照提前設(shè)置的優(yōu)化問(wèn)題的適應(yīng)度函數(shù)計(jì)算每個(gè)個(gè)體(每個(gè)解)的適應(yīng)度,再根據(jù)適應(yīng)度值對(duì)每個(gè)個(gè)體對(duì)應(yīng)的染色體進(jìn)行選擇,淘汰掉適應(yīng)度值低的染色體,適應(yīng)度值高的染色體就會(huì)保留下來(lái),然后對(duì)染色體進(jìn)行交叉、變異等遺傳算法,進(jìn)而產(chǎn)生新一代的種群,這樣一代又一代的計(jì)算,最終達(dá)到滿足某種優(yōu)化約束條件而生成最優(yōu)的群體,即是種群中最優(yōu)解。

    3.1 球鼻首優(yōu)化的具體步驟

    遺傳算法是一個(gè)以適應(yīng)度函數(shù)(目標(biāo)函數(shù))為依據(jù),通過(guò)對(duì)種群個(gè)體施加遺傳操作實(shí)現(xiàn)種群內(nèi)個(gè)體結(jié)構(gòu)重組的迭代處理過(guò)程。在這一過(guò)程中,群體個(gè)體(問(wèn)題的解)逐代優(yōu)化并逐漸逼近最優(yōu)解。使用遺傳算法NSGAII對(duì)球鼻首型線進(jìn)行優(yōu)化分析,應(yīng)用于球鼻首的型線優(yōu)化中的具體步驟如圖2所示。

    3.2 優(yōu)化結(jié)果分析

    在Shipflow進(jìn)行計(jì)算時(shí),需要導(dǎo)入相應(yīng)的實(shí)船或船模尺度下的模型格式,接著設(shè)置好相關(guān)參數(shù)進(jìn)行計(jì)算,其結(jié)果可以輸出船模尺度下的,也可以輸出實(shí)船尺度下的。文章中在進(jìn)行計(jì)算結(jié)果分析時(shí),導(dǎo)入實(shí)船模型對(duì)球鼻首進(jìn)行參數(shù)化建模并分析靜水力數(shù)據(jù),并使用CAESES集成Shipflow軟件計(jì)算,輸出的優(yōu)化結(jié)果為船模尺度下的解。

    3.2.1 設(shè)計(jì)速度下優(yōu)化結(jié)果分析

    使用CAESES軟件集成Shipflow軟件,運(yùn)用遺傳算法對(duì)參數(shù)化之后的2 200箱集裝箱船球鼻首在設(shè)計(jì)吃水、設(shè)計(jì)航速下進(jìn)行優(yōu)化計(jì)算。使用建立好的優(yōu)化模型,選擇球鼻首建模時(shí)的26個(gè)參數(shù)作為設(shè)計(jì)變量,無(wú)因次量排水體積V′和浮心縱向位置XCOF作為約束條件,興波阻力作為目標(biāo)函數(shù)。利用遺傳算法NSGAII進(jìn)行優(yōu)化,設(shè)置種群數(shù)量為28,遺傳代數(shù)為20,交叉概率為0.8,變異概率為0.008,生成560組船型設(shè)計(jì)方案,剔除一部分出錯(cuò)的結(jié)果,得到模型優(yōu)化之后的最優(yōu)興波阻力方案。計(jì)算并對(duì)比優(yōu)化前后的目標(biāo)和約束條件如表3所示。

    由表3可以看出,和初始方案相比,在船型設(shè)計(jì)吃水,航速為19 kn的情況下,經(jīng)過(guò)NSGAII算法優(yōu)化后的船型興波阻力下降了4.137%,興波阻力系數(shù)降低3.822%,而浮心縱向位置和排水體積的變化都維持在1%以內(nèi),誤差計(jì)算同式(1)。圖3為優(yōu)化前后球鼻首面元對(duì)比圖,圖4為優(yōu)化前后球鼻首部分型值組橫剖線對(duì)比圖,由圖3和圖4可以看出,優(yōu)化之后的球鼻首相對(duì)初始球鼻首,長(zhǎng)度有明顯的增加,并且有上翹的趨勢(shì)。圖5和圖6分別為航速19kn時(shí),初始船型和優(yōu)化后船型沿船體表面自由面波高對(duì)比圖和波形對(duì)比圖,比較圖5和圖6可以看出,優(yōu)化之后的船型在球鼻首處的波高明顯降低,并且船中后部的散波也有所改善。

    表3 優(yōu)化前船型和優(yōu)化后船型數(shù)據(jù)對(duì)比

    3.2.2 優(yōu)化船型低速航行性能分析

    集裝箱船逐漸向低速方向發(fā)展,一方面由于集裝箱船在低速時(shí)能夠減少運(yùn)營(yíng)成本;另一方面,在全球海運(yùn)低迷期,低速行駛可以增加船舶的使用量,讓更多閑置或者新增運(yùn)力投入運(yùn)營(yíng)[9, 10]。對(duì)優(yōu)化之后的船型在15 kn進(jìn)行計(jì)算,并和原始船型進(jìn)行對(duì)比。

    計(jì)算結(jié)果顯示,航速為15 kn時(shí),優(yōu)化船型的興波阻力系數(shù)相對(duì)初始船型下降了9.208%;總阻力系數(shù)也下降了1.010%,如表4所示。圖7和圖8分別為航速15 kn時(shí),初始船型和優(yōu)化后船型沿船體表面自由面波高對(duì)比圖和波形對(duì)比圖。由圖7可以看出,在球鼻首優(yōu)化之后船體中前部的波高有明顯降低,由圖8也可以看出優(yōu)化之后船型的波形有明顯的改善。

    表4 初始船型和優(yōu)化后船型數(shù)據(jù)對(duì)比(v=15 kn)

    4 結(jié) 語(yǔ)

    使用CAESES軟件對(duì)2 200箱集裝箱船球鼻首進(jìn)行了參數(shù)化建模,獲得參數(shù)化之后的模型,基于Shipflow和CAESES軟件集成平臺(tái),結(jié)合遺傳算法對(duì)2 200箱集裝箱船球鼻首在設(shè)計(jì)吃水,航速為19 kn下的型線進(jìn)行了優(yōu)化,并將優(yōu)化好的船型在低速航行狀態(tài)15 kn下進(jìn)行了阻力驗(yàn)證。結(jié)果證明,最終獲得阻力性能良好的球鼻首型線,提高該船型的水動(dòng)力性能。

    [1]郝金鳳, 強(qiáng)兆新, 石俊令, 等.船舶設(shè)計(jì)節(jié)能減排技術(shù)策略[J].艦船科學(xué)技術(shù), 2012, 09: 3–10.

    [2]孫存樓, 徐炳香, 韓文玉.基于母型船的船體線型優(yōu)化設(shè)計(jì)研究[J].船舶工程, 2015, 01: 20–23.

    [3]HYUNYUL KIM, CHI Yang.A new surface modification approach for CFD-based hull form optimization [C]// 9th International Conferences on Hydrodynamics, 2010, 10: 11-15.

    [4]胡春平.基于Friendship的船體型線優(yōu)化技術(shù)研究[D].武漢:武漢理工大學(xué), 2012.

    [5]黃金森.基于最小阻力的船型優(yōu)化研究[D].天津: 天津大學(xué),2006.

    [6]鄧賢輝, 方昭昭, 趙丙乾, 等.基于計(jì)算流體動(dòng)力學(xué)的最小阻力船型自動(dòng)優(yōu)化[J].中國(guó)艦船研究, 2015, 10(3): 19–25.

    [7]謝宜, 王長(zhǎng)雷.基于小生境遺傳算法的單體船船型參數(shù)優(yōu)化[J].艦船科學(xué)技術(shù), 2011, 33(3): 35–38.

    [8]李華昌, 謝淑蘭, 易忠勝.遺傳算法的原理與應(yīng)用[J].礦冶,2005(01): 87–90.

    [9]辛吉誠(chéng).未來(lái)大型集 裝 箱 船 發(fā)展趨勢(shì)[J].中國(guó)船檢,2015(02): 92–95.

    [10]張華, 王超峰.超低速航行及其經(jīng)濟(jì)效益淺析[J].航海,2010(02): 44–46.

    猜你喜歡
    興波集裝箱船船型
    巨大的集裝箱船
    全球最大集裝箱船首航青島港
    軍事文摘(2023年14期)2023-07-28 08:39:46
    船舶興波運(yùn)動(dòng)的非線性現(xiàn)象算法研究與仿真
    基于Friendship的高速客船特征參數(shù)靈敏度分析
    世界最大級(jí)別集裝箱船“宇宙號(hào)”
    軍事文摘(2018年24期)2018-12-26 00:57:56
    中國(guó)古代著名的海船船型
    ——福船
    中國(guó)古代著名的海船船型—廣船
    Anti-ageing effects of a new Dimethylaminoethanol-based formulation on DGalactose induced skin ageing model of rat
    內(nèi)傾船型的非線性橫搖仿真分析
    影響潛航體水面興波因素分析
    一区二区三区高清视频在线| 午夜精品国产一区二区电影 | 国产成人a∨麻豆精品| 春色校园在线视频观看| 久久久精品94久久精品| 亚洲无线在线观看| 亚洲国产欧美人成| 国产精品亚洲一级av第二区| 美女高潮的动态| 精品无人区乱码1区二区| 国产成人aa在线观看| 国产久久久一区二区三区| 一级av片app| 人妻久久中文字幕网| 老司机午夜福利在线观看视频| 亚洲精品乱码久久久v下载方式| 国产91av在线免费观看| 亚洲av电影不卡..在线观看| 亚洲国产精品成人久久小说 | 少妇猛男粗大的猛烈进出视频 | 成人一区二区视频在线观看| 欧美3d第一页| 亚洲七黄色美女视频| 两个人的视频大全免费| 最近手机中文字幕大全| 亚洲自偷自拍三级| 22中文网久久字幕| 亚洲四区av| 一个人观看的视频www高清免费观看| 亚洲精品色激情综合| 12—13女人毛片做爰片一| 搞女人的毛片| 成人一区二区视频在线观看| 欧美日本视频| 国产伦精品一区二区三区四那| 久99久视频精品免费| 在现免费观看毛片| 国产成人aa在线观看| 精品欧美国产一区二区三| 深爱激情五月婷婷| 一a级毛片在线观看| 色综合亚洲欧美另类图片| 成人特级黄色片久久久久久久| а√天堂www在线а√下载| 99热网站在线观看| 一区福利在线观看| 亚洲av中文av极速乱| 国产 一区精品| 男人舔奶头视频| 老司机福利观看| 毛片一级片免费看久久久久| 久久久精品94久久精品| 精品午夜福利在线看| 国产精品福利在线免费观看| 中文字幕av成人在线电影| 看非洲黑人一级黄片| 免费在线观看影片大全网站| www日本黄色视频网| 国产三级中文精品| 欧美三级亚洲精品| 全区人妻精品视频| 国产午夜福利久久久久久| 老熟妇仑乱视频hdxx| 国产精品久久久久久亚洲av鲁大| 亚洲成人精品中文字幕电影| eeuss影院久久| 国产精品久久久久久精品电影| 久久久久久久久久黄片| 麻豆国产97在线/欧美| 俄罗斯特黄特色一大片| 舔av片在线| 欧美日韩乱码在线| 久久久久久国产a免费观看| 日韩精品青青久久久久久| 我的老师免费观看完整版| 欧美+亚洲+日韩+国产| 日本熟妇午夜| 精品少妇黑人巨大在线播放 | 欧美三级亚洲精品| 亚洲欧美日韩无卡精品| 99久久精品热视频| 深爱激情五月婷婷| 国内精品宾馆在线| 亚洲欧美精品综合久久99| 麻豆一二三区av精品| 婷婷精品国产亚洲av| 小说图片视频综合网站| 精品人妻偷拍中文字幕| 最新中文字幕久久久久| 久久久久久国产a免费观看| 97热精品久久久久久| 免费观看的影片在线观看| 最近在线观看免费完整版| 成年版毛片免费区| 精品少妇黑人巨大在线播放 | a级毛片免费高清观看在线播放| 久久精品影院6| 不卡视频在线观看欧美| 午夜视频国产福利| 日本熟妇午夜| 亚洲精品一卡2卡三卡4卡5卡| 日韩精品中文字幕看吧| av天堂中文字幕网| 午夜福利在线观看免费完整高清在 | 成年女人看的毛片在线观看| 中出人妻视频一区二区| 免费看av在线观看网站| 如何舔出高潮| 亚洲国产精品久久男人天堂| 国产精品一区二区免费欧美| 天堂网av新在线| 赤兔流量卡办理| 伦理电影大哥的女人| 国产精品乱码一区二三区的特点| 国产成人一区二区在线| av在线亚洲专区| 国产亚洲精品综合一区在线观看| 99精品在免费线老司机午夜| 精品一区二区三区av网在线观看| 国产精品电影一区二区三区| 午夜免费男女啪啪视频观看 | 国产欧美日韩精品亚洲av| 成年女人永久免费观看视频| 久久综合国产亚洲精品| 欧美zozozo另类| 亚洲精品国产av成人精品 | 国产精品一区二区性色av| 男人舔奶头视频| 美女免费视频网站| 久久精品国产清高在天天线| 女生性感内裤真人,穿戴方法视频| 国产黄a三级三级三级人| 乱人视频在线观看| 国产一区二区三区在线臀色熟女| 亚洲最大成人av| 18禁在线无遮挡免费观看视频 | 天天一区二区日本电影三级| .国产精品久久| 此物有八面人人有两片| 国产精品久久久久久精品电影| 九色成人免费人妻av| 69人妻影院| 美女cb高潮喷水在线观看| 一级黄片播放器| 久久久久精品国产欧美久久久| 精品久久久久久久久久久久久| 久久久久久久亚洲中文字幕| 人人妻人人看人人澡| 伊人久久精品亚洲午夜| 中文字幕av成人在线电影| 一a级毛片在线观看| 国产av不卡久久| 亚洲av第一区精品v没综合| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品久久久com| 午夜激情欧美在线| 亚洲第一电影网av| 午夜免费男女啪啪视频观看 | 国产精品人妻久久久影院| 色哟哟·www| 日本a在线网址| 色播亚洲综合网| 日韩亚洲欧美综合| 亚洲精品456在线播放app| 在线播放国产精品三级| 国产精品野战在线观看| 尾随美女入室| 在线a可以看的网站| 国产女主播在线喷水免费视频网站 | 美女内射精品一级片tv| 亚洲av第一区精品v没综合| 性色avwww在线观看| 国内精品宾馆在线| 乱系列少妇在线播放| 桃色一区二区三区在线观看| 乱系列少妇在线播放| 高清毛片免费看| 久久久精品大字幕| 免费人成在线观看视频色| 亚洲av五月六月丁香网| 国产精品综合久久久久久久免费| 午夜福利在线在线| 久久国产乱子免费精品| 国内揄拍国产精品人妻在线| 精品乱码久久久久久99久播| 久久精品影院6| 偷拍熟女少妇极品色| 我要搜黄色片| 久久九九热精品免费| 欧美极品一区二区三区四区| 可以在线观看毛片的网站| 国产伦在线观看视频一区| 成人亚洲欧美一区二区av| 午夜福利在线观看免费完整高清在 | 国产又黄又爽又无遮挡在线| 亚洲国产精品合色在线| 99视频精品全部免费 在线| 欧美又色又爽又黄视频| 高清午夜精品一区二区三区 | 国产午夜精品论理片| 一区二区三区高清视频在线| 赤兔流量卡办理| 免费av不卡在线播放| 午夜福利在线在线| 97超视频在线观看视频| 国产欧美日韩精品亚洲av| 色播亚洲综合网| 国产麻豆成人av免费视频| 人妻少妇偷人精品九色| 久久99热6这里只有精品| 国产精品久久久久久久久免| 最近最新中文字幕大全电影3| 精品久久久久久久人妻蜜臀av| 国产精品国产高清国产av| 小说图片视频综合网站| 亚洲乱码一区二区免费版| 噜噜噜噜噜久久久久久91| 一个人免费在线观看电影| 中国美女看黄片| 久99久视频精品免费| 久久精品国产亚洲av天美| 日韩,欧美,国产一区二区三区 | 午夜爱爱视频在线播放| 色av中文字幕| 免费一级毛片在线播放高清视频| 一级黄色大片毛片| 色哟哟哟哟哟哟| 午夜福利18| 日韩欧美三级三区| 日韩强制内射视频| 18禁裸乳无遮挡免费网站照片| 国产日本99.免费观看| 国产高清三级在线| 欧美三级亚洲精品| 亚洲精品日韩av片在线观看| 久久人人爽人人爽人人片va| 丝袜美腿在线中文| 性欧美人与动物交配| 看免费成人av毛片| 日本五十路高清| 亚洲精品国产成人久久av| 99久国产av精品国产电影| 又黄又爽又刺激的免费视频.| 精品不卡国产一区二区三区| 国产精品国产三级国产av玫瑰| 六月丁香七月| 日本三级黄在线观看| 网址你懂的国产日韩在线| 国产av在哪里看| 久久久久久大精品| 亚洲人与动物交配视频| 午夜免费男女啪啪视频观看 | 成人无遮挡网站| 久久久国产成人精品二区| 少妇的逼水好多| 在线观看66精品国产| 2021天堂中文幕一二区在线观| 国产精品人妻久久久影院| 国产精品久久久久久av不卡| 性欧美人与动物交配| 搞女人的毛片| 在线观看一区二区三区| a级毛片a级免费在线| 午夜老司机福利剧场| 一级黄片播放器| 久久久久久久亚洲中文字幕| 国产探花极品一区二区| 12—13女人毛片做爰片一| 久久99热这里只有精品18| 超碰av人人做人人爽久久| 1024手机看黄色片| 99热这里只有精品一区| 国产午夜精品论理片| 91久久精品电影网| 久久中文看片网| 天美传媒精品一区二区| 精品人妻熟女av久视频| 亚洲第一电影网av| 国产精品99久久久久久久久| 国产真实伦视频高清在线观看| 亚洲欧美日韩卡通动漫| 日本撒尿小便嘘嘘汇集6| 少妇的逼水好多| 国产精品综合久久久久久久免费| 免费不卡的大黄色大毛片视频在线观看 | 午夜福利在线在线| 麻豆国产97在线/欧美| 男人的好看免费观看在线视频| 别揉我奶头 嗯啊视频| 亚洲专区国产一区二区| 村上凉子中文字幕在线| 日韩强制内射视频| 99热6这里只有精品| 国产av麻豆久久久久久久| 久久久久久九九精品二区国产| 高清日韩中文字幕在线| 69av精品久久久久久| 欧美中文日本在线观看视频| 久久亚洲国产成人精品v| 久久久久久九九精品二区国产| av免费在线看不卡| 嫩草影院入口| 国产91av在线免费观看| 国产精品国产高清国产av| 国产精品无大码| 亚洲性夜色夜夜综合| 麻豆乱淫一区二区| 在线播放无遮挡| 22中文网久久字幕| 变态另类丝袜制服| 午夜精品在线福利| 亚洲人成网站高清观看| 国产高潮美女av| 精品久久国产蜜桃| 国产淫片久久久久久久久| 成人美女网站在线观看视频| 欧美激情久久久久久爽电影| 国国产精品蜜臀av免费| 在线免费十八禁| 人人妻人人看人人澡| 最好的美女福利视频网| 亚洲自偷自拍三级| 国模一区二区三区四区视频| 麻豆av噜噜一区二区三区| 欧美人与善性xxx| 深夜a级毛片| 日日干狠狠操夜夜爽| 又爽又黄无遮挡网站| 亚洲av五月六月丁香网| 色综合亚洲欧美另类图片| 欧美日韩乱码在线| 国产日本99.免费观看| 日韩一区二区视频免费看| 成人亚洲欧美一区二区av| 久久精品久久久久久噜噜老黄 | 亚洲天堂国产精品一区在线| 一个人看视频在线观看www免费| 日产精品乱码卡一卡2卡三| 丰满乱子伦码专区| 亚洲人与动物交配视频| 在线播放无遮挡| 日韩精品青青久久久久久| 亚洲精品粉嫩美女一区| 欧美xxxx性猛交bbbb| 日韩一区二区视频免费看| 国产高清视频在线播放一区| 国产精品三级大全| 国产精品美女特级片免费视频播放器| 欧美潮喷喷水| 国产女主播在线喷水免费视频网站 | 精品一区二区三区视频在线观看免费| 亚洲av美国av| 少妇人妻精品综合一区二区 | 精品人妻熟女av久视频| 看黄色毛片网站| 欧美又色又爽又黄视频| 搡老妇女老女人老熟妇| 亚洲精品乱码久久久v下载方式| 中文亚洲av片在线观看爽| 国产精品日韩av在线免费观看| 99热6这里只有精品| 国产成人一区二区在线| 美女大奶头视频| 免费搜索国产男女视频| 老师上课跳d突然被开到最大视频| 国产极品精品免费视频能看的| 天堂√8在线中文| 日韩 亚洲 欧美在线| 国产亚洲精品av在线| 少妇被粗大猛烈的视频| 成人午夜高清在线视频| 丝袜喷水一区| av中文乱码字幕在线| 日韩av在线大香蕉| 日韩欧美国产在线观看| 国产精品一区二区性色av| 成人漫画全彩无遮挡| 亚洲欧美精品综合久久99| 黄色一级大片看看| 国产成年人精品一区二区| 久久国内精品自在自线图片| 中文在线观看免费www的网站| 日韩精品青青久久久久久| 中国国产av一级| 尤物成人国产欧美一区二区三区| 在线天堂最新版资源| 精品久久国产蜜桃| 国产日本99.免费观看| 欧美潮喷喷水| 日日啪夜夜撸| 欧美日韩乱码在线| 91精品国产九色| 欧美日韩综合久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 夜夜看夜夜爽夜夜摸| 波多野结衣高清作品| 亚洲国产色片| 在线播放无遮挡| eeuss影院久久| 国产成人影院久久av| 国产在线精品亚洲第一网站| 日韩欧美 国产精品| 亚洲无线在线观看| 国产69精品久久久久777片| 亚洲精品一卡2卡三卡4卡5卡| 亚州av有码| 久久精品夜色国产| 亚洲av五月六月丁香网| 哪里可以看免费的av片| aaaaa片日本免费| 欧美最新免费一区二区三区| 久久久精品94久久精品| 久久草成人影院| 一个人看的www免费观看视频| 国产男人的电影天堂91| 麻豆精品久久久久久蜜桃| 国语自产精品视频在线第100页| 亚洲国产欧洲综合997久久,| 欧美另类亚洲清纯唯美| 免费在线观看影片大全网站| 欧美一区二区亚洲| 欧美高清性xxxxhd video| av中文乱码字幕在线| 麻豆国产av国片精品| 国产免费男女视频| 国产精品人妻久久久影院| 可以在线观看的亚洲视频| 在线国产一区二区在线| 精品久久久久久久久亚洲| 亚洲高清免费不卡视频| 国产综合懂色| 熟女电影av网| 久久久精品欧美日韩精品| 最近最新中文字幕大全电影3| ponron亚洲| 欧美成人a在线观看| 久久久久久久久久久丰满| 22中文网久久字幕| 亚洲aⅴ乱码一区二区在线播放| av专区在线播放| 亚洲乱码一区二区免费版| 日韩欧美精品免费久久| 国产精品野战在线观看| 3wmmmm亚洲av在线观看| 春色校园在线视频观看| 真人做人爱边吃奶动态| 久久精品国产亚洲av香蕉五月| 麻豆乱淫一区二区| 老司机福利观看| 你懂的网址亚洲精品在线观看 | 天堂网av新在线| 搡老岳熟女国产| 99久国产av精品国产电影| 不卡视频在线观看欧美| 国产黄a三级三级三级人| 内地一区二区视频在线| 精品一区二区三区人妻视频| 91在线精品国自产拍蜜月| 欧美最黄视频在线播放免费| 久久99热这里只有精品18| 亚洲av免费高清在线观看| 欧美日本视频| 久久这里只有精品中国| 51国产日韩欧美| 在线天堂最新版资源| 国产精品精品国产色婷婷| 中文字幕精品亚洲无线码一区| 国产白丝娇喘喷水9色精品| 精品久久久久久久久av| 日本黄色片子视频| 欧美日韩一区二区视频在线观看视频在线 | 又粗又爽又猛毛片免费看| 两个人视频免费观看高清| 国产精品永久免费网站| 又黄又爽又免费观看的视频| 精品久久久久久久久亚洲| 欧美激情在线99| 一级毛片aaaaaa免费看小| 人人妻人人澡欧美一区二区| 一个人免费在线观看电影| 精品一区二区三区av网在线观看| 亚洲精品日韩av片在线观看| 熟女人妻精品中文字幕| 国产高清不卡午夜福利| 长腿黑丝高跟| 在线播放国产精品三级| 婷婷精品国产亚洲av在线| 禁无遮挡网站| 久久久a久久爽久久v久久| av天堂在线播放| 日韩高清综合在线| 国产成人一区二区在线| 久久久久久久久久黄片| 精品午夜福利视频在线观看一区| 精品一区二区三区视频在线观看免费| 中文字幕熟女人妻在线| 一本精品99久久精品77| 国产午夜精品久久久久久一区二区三区 | 国产精品人妻久久久影院| 成人亚洲精品av一区二区| 亚洲av中文av极速乱| 久久久国产成人免费| 香蕉av资源在线| 91精品国产九色| 日本免费一区二区三区高清不卡| 搡老妇女老女人老熟妇| 国产精品人妻久久久久久| 国产亚洲91精品色在线| 麻豆av噜噜一区二区三区| 国产乱人视频| 成人性生交大片免费视频hd| 成人高潮视频无遮挡免费网站| 黄片wwwwww| 久久久久性生活片| 国产爱豆传媒在线观看| 亚洲精品国产av成人精品 | 日本 av在线| 99热6这里只有精品| 国产极品精品免费视频能看的| 免费观看人在逋| 啦啦啦观看免费观看视频高清| 亚洲精品成人久久久久久| 国产麻豆成人av免费视频| 欧美三级亚洲精品| 国产精品久久视频播放| 亚洲,欧美,日韩| 久久精品影院6| 天堂影院成人在线观看| 久久99热这里只有精品18| 日本撒尿小便嘘嘘汇集6| 国产精品综合久久久久久久免费| 日韩欧美精品免费久久| 国内精品一区二区在线观看| 舔av片在线| 美女被艹到高潮喷水动态| 香蕉av资源在线| 精品99又大又爽又粗少妇毛片| 国产精品无大码| 在线看三级毛片| 此物有八面人人有两片| 国产片特级美女逼逼视频| 亚洲第一区二区三区不卡| 一本久久中文字幕| 三级男女做爰猛烈吃奶摸视频| 午夜亚洲福利在线播放| 内射极品少妇av片p| 免费观看人在逋| 亚洲,欧美,日韩| 日韩强制内射视频| 美女被艹到高潮喷水动态| 久久久久久国产a免费观看| 成人美女网站在线观看视频| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品不卡视频一区二区| 国产av不卡久久| 亚洲av成人精品一区久久| 国产成年人精品一区二区| 亚洲在线自拍视频| 精品午夜福利在线看| 国产色婷婷99| 一进一出好大好爽视频| 亚洲国产精品国产精品| 欧美最新免费一区二区三区| 99国产精品一区二区蜜桃av| 久久精品国产亚洲av涩爱 | 久久久久久久久久黄片| 亚洲国产精品国产精品| 少妇丰满av| 欧美成人精品欧美一级黄| 亚洲乱码一区二区免费版| 韩国av在线不卡| 在线观看av片永久免费下载| 麻豆成人午夜福利视频| 欧美三级亚洲精品| 国产成人影院久久av| 日本黄大片高清| 99久久精品国产国产毛片| 插逼视频在线观看| 亚洲高清免费不卡视频| 国产黄片美女视频| 免费人成视频x8x8入口观看| 熟女电影av网| 久久国内精品自在自线图片| 国产视频一区二区在线看| 一边摸一边抽搐一进一小说| 欧美日韩综合久久久久久| 美女内射精品一级片tv| 亚洲av二区三区四区| av卡一久久| 亚洲色图av天堂| 亚洲五月天丁香| av福利片在线观看| 亚洲欧美清纯卡通| 国产成人精品久久久久久| 亚洲国产色片| 午夜福利在线在线| 免费在线观看成人毛片| 一进一出抽搐gif免费好疼| 国产美女午夜福利| 此物有八面人人有两片| 国产日本99.免费观看| 国产三级在线视频| 99热全是精品| 美女免费视频网站| 又爽又黄无遮挡网站| 精品久久久久久久久亚洲| 99热这里只有精品一区| 一区福利在线观看| 亚洲图色成人| 亚洲av免费在线观看| 欧美最黄视频在线播放免费| 国产精品爽爽va在线观看网站| 成人高潮视频无遮挡免费网站| 日本爱情动作片www.在线观看 | 床上黄色一级片|