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

    基于差分進(jìn)化算法的散貨船“最后分艙裝載”計(jì)算

    2016-10-12 01:32:47,,,
    中國航海 2016年4期
    關(guān)鍵詞:貨艙散貨船差分

    , , ,

    (大連海事大學(xué) 航海動(dòng)態(tài)仿真和控制交通行業(yè)重點(diǎn)實(shí)驗(yàn)室, 遼寧 大連 116026)

    基于差分進(jìn)化算法的散貨船“最后分艙裝載”計(jì)算

    劉春雷,孫霄峰,吳海波,尹勇

    (大連海事大學(xué)航海動(dòng)態(tài)仿真和控制交通行業(yè)重點(diǎn)實(shí)驗(yàn)室,遼寧大連116026)

    為實(shí)現(xiàn)散貨船“最分分艙裝載”計(jì)算,根據(jù)其特點(diǎn)建立“最后分艙裝載”的優(yōu)化模型,將其歸結(jié)為單目標(biāo)約束優(yōu)化問題并采用差分進(jìn)化算法進(jìn)行求解。以某25萬t散貨船為例進(jìn)行實(shí)例計(jì)算:實(shí)例一為平均吃水變化較大時(shí),利用傳統(tǒng)“縱傾力矩法”與上述方法計(jì)算得到的平均吃水和吃水差的誤差分別為1.026%及8.33%,0%及0%;實(shí)例二為平均吃水變化較小時(shí),利用傳統(tǒng)“縱傾力矩法”與上述方法計(jì)算得到的平均吃水和吃水差的誤差分別為0%和2.17%,0%和0%。計(jì)算結(jié)果表明:上述方法不僅適用于平均吃水變化較小時(shí)的“最后分艙裝載”計(jì)算,而且適用于平均吃水變化較大時(shí)的計(jì)算,尤其是在平均吃水變化較大時(shí)可有效減小誤差。

    水路運(yùn)輸; 散貨船; 配載儀; 最后分艙裝載; 差分進(jìn)化算法

    Abstract: An optimization model for bulk carrier “Final Trimming” calculation is built according to the characteristics of the process. The trimming process is regarded as a single-objective optimization with constraints. The differential evolution algorithm is used to solve the problem. Two cases of trimming for a 250 000 tons bulk carrier are carried out for illustration. The first case, associated with large mean draft changes, shows that the deviations of mean draft and trim of draft are 1.026% and 8.33% by traditional “Trimming Moment Method”. The second case, representing the final trimming with only small mean draft change, shows that the deviations of mean draft and trim of draft are 0% and 2.17% by traditional method. On the contrary, the suggested method makes all the deviations zero in both cases, which means that the method can be effective for both general loading and final trimming.

    Keywords: waterway transportation; bulk carrier; loading computer; final trimming; differential evolution algorithm

    隨著航運(yùn)業(yè)不斷發(fā)展,散貨船逐漸趨于大型化,這使得其進(jìn)出港口時(shí)的吃水常受泊位水深的限制。散貨船裝載貨物時(shí)通常按預(yù)定的配載方案裝載,理論上可達(dá)到要求的浮態(tài),但實(shí)際裝載時(shí)因裝載量和裝載位置與理論值不同,導(dǎo)致裝貨后的浮態(tài)與所要求的浮態(tài)有較大差別。因此,散貨船裝貨時(shí)通常會(huì)預(yù)留一部分貨物用于“最后分艙裝載”[1],使船舶最終的浮態(tài)滿足進(jìn)出港要求。

    目前針對(duì)散貨船“最后分艙裝載”的計(jì)算大多是使用Excel根據(jù)“加載100 t吃水法”[2]計(jì)算或根據(jù)“縱傾力矩法”[3]計(jì)算。這2種方法存在每厘米吃水噸數(shù)TPC和每厘米縱傾力矩MTC一直不變的假設(shè),當(dāng)平均吃水和吃水差變化較大時(shí)會(huì)產(chǎn)生較大的誤差。此外,文獻(xiàn)[4]采用“內(nèi)點(diǎn)法”解決質(zhì)量分配問題,將艙室簡(jiǎn)化成方體,認(rèn)為其重心是不依賴于載重量的,而實(shí)際上艙室是不規(guī)則體,其重心和載重量存在著一一對(duì)應(yīng)的關(guān)系,這樣的簡(jiǎn)化也會(huì)導(dǎo)致“最后分艙裝載”計(jì)算產(chǎn)生誤差。

    針對(duì)以上不足,將散貨船“最后分艙裝載”歸結(jié)為單目標(biāo)約束優(yōu)化問題,建立散貨船“最后分艙裝載”的優(yōu)化模型,并將差分進(jìn)化算法引入到該優(yōu)化模型的求解中,通過建立離線的艙容表實(shí)現(xiàn)散貨船的“最后分艙裝載”計(jì)算,提高計(jì)算精度。

    1 建立“最后分艙裝載”優(yōu)化模型

    “最后分艙裝載”結(jié)束后,船舶在重力和浮力的作用下能平衡漂浮于靜水面,此時(shí)可能同時(shí)存在橫傾和縱傾,但總滿足任意浮態(tài)下的平衡方程。船舶任意浮態(tài)平衡時(shí)所確定的平衡方程組[5]為

    (1)

    1.1 設(shè)計(jì)變量

    1.2 目標(biāo)函數(shù)

    通?!白詈蠓峙撗b載”前船舶的浮態(tài)和各個(gè)貨艙的裝載情況是已知的;而“最后分艙裝載”結(jié)束后船舶處在一個(gè)平衡狀態(tài),滿足船舶任意狀態(tài)下的平衡方程組(該平衡方程組為一個(gè)隱式非線性方程組)。這里主要使用差分進(jìn)化算法求解該方程組,進(jìn)而求出“最后分艙裝載”的貨物分配量,因此平衡方程組中的任何一個(gè)方程都可作為目標(biāo)函數(shù)。但是,出于計(jì)算簡(jiǎn)單考慮,根據(jù)平衡方程組中的第一式,以用于“最后分艙裝載”的貨艙總載貨量與總排水量變化量差值的絕對(duì)值作為目標(biāo)函數(shù),即

    (2)

    式(2)中:Mlast為分艙裝載后船舶總的排水量;M0為分艙裝載前除調(diào)整艙之外船舶總的排水量。

    1.3約束條件

    根據(jù)船舶平衡時(shí)所確定的浮態(tài)方程、每個(gè)艙室允許裝載的最大裝貨量和最小裝貨量,得到以下約束條件。

    (3)

    式(3)中:pimax為要調(diào)整的艙中第i個(gè)貨艙的最大裝貨量;xG0,yG0和zG0分別為調(diào)整吃水前除調(diào)整艙之外船舶重心的縱向坐標(biāo)、橫向坐標(biāo)及垂向坐標(biāo);xGi,yGi和zGi分別為要調(diào)整的艙中第i個(gè)貨艙重心的縱向坐標(biāo)、橫向坐標(biāo)及垂向坐標(biāo);xBlast,yBlast和zBlast分別為調(diào)整后所需的船舶浮心的縱向坐標(biāo)、橫向坐標(biāo)及垂向坐標(biāo)。

    2 用于“最后分艙裝載”的差分進(jìn)化算法

    差分進(jìn)化算法具有結(jié)構(gòu)簡(jiǎn)單、可調(diào)參數(shù)少、收斂速度快和魯棒性強(qiáng)等特點(diǎn)。[6-8]該算法作為一類新興的優(yōu)化算法,在處理單目標(biāo)優(yōu)化問題方面具有很強(qiáng)的優(yōu)勢(shì),已廣泛運(yùn)用于科學(xué)研究和工程實(shí)踐當(dāng)中,如控制器參數(shù)優(yōu)化、濾波器設(shè)計(jì)、數(shù)據(jù)挖掘和石油化工等。[9]

    差分進(jìn)化算法是從一個(gè)隨機(jī)初始化的種群開始搜索,經(jīng)過變異、交叉和選擇操作產(chǎn)生下一時(shí)刻的主種群;該過程重復(fù)進(jìn)行,直到滿足停止條件。

    2.1 編碼及初始化過程

    差分進(jìn)化算法中的個(gè)體即優(yōu)化設(shè)計(jì)中的設(shè)計(jì)變量,個(gè)體上的每一維變量即為每個(gè)設(shè)計(jì)變量。每個(gè)個(gè)體的維數(shù)與設(shè)計(jì)變量相同。

    1) 設(shè)計(jì)變量

    (4)

    2) 種群個(gè)體

    (5)

    式(4)和式(5)中:n為個(gè)體的維數(shù);NP為種群規(guī)模。

    初始種群應(yīng)盡可能地覆蓋整個(gè)搜索域,用隨機(jī)方法產(chǎn)生,每個(gè)個(gè)體Xi的每一維生成方法為

    xi,j=pjmin+rand(0,1)×(pjmax-pjmin)

    (6)

    式(6)中:rand(0,1)為一個(gè)位于0~1的滿足均勻分布的隨機(jī)數(shù);pjmin,pjmax分別為初始化用于“最后分艙裝載”的貨艙中第j個(gè)貨艙載重量的最大值和最小值,取pjmin=-0.3pjmax,pjmax=1.3pjmax。

    2.2 評(píng)價(jià)過程

    由于約束的存在,在計(jì)算每個(gè)個(gè)體目標(biāo)函數(shù)值的同時(shí)還要計(jì)算個(gè)體的違反約束值[10],以區(qū)分可行個(gè)體與不可行個(gè)體。

    計(jì)算每個(gè)個(gè)體的目標(biāo)函數(shù)值為

    (7)

    計(jì)算每個(gè)個(gè)體的違反約束值,將約束條件改寫成式(8),其中g(shù)r(Xi)及hr(Xi)分別為不等式約束和等式約束;q及m-q分別為不等式約束的數(shù)量和等式約束的數(shù)量,這里q=2n,m=2n+2。

    (8)

    個(gè)體違反第r個(gè)約束的值表示為

    (9)

    式(9)中:δ為等式約束的容許誤差;δ的取值范圍為0.001~0.010。由此可得個(gè)體的違反約束值為

    (10)

    當(dāng)G(Xi)=0時(shí),Xi為可行個(gè)體;否則Xi為不可行個(gè)體。

    2.3 變異過程

    對(duì)于種群中每個(gè)個(gè)體Xi,隨機(jī)生成3個(gè)互不相同的整數(shù)r1,r2,r3∈{1,2,…,NP},且要求r1,r2,r3,i互不相同,按照式(11)產(chǎn)生變異個(gè)體Vi。

    Vi=Xr1+F·(Xr2-Xr3)

    (11)

    式(11)中:縮放因子F為正數(shù),通常被限定在區(qū)間[0,1]內(nèi)。

    2.4 交叉過程

    交叉操作按每個(gè)個(gè)體向量的分量進(jìn)行。變異產(chǎn)生的變異個(gè)體Vi和目標(biāo)個(gè)體Xi按式(12)進(jìn)行交叉操作產(chǎn)生試驗(yàn)個(gè)體Ui。

    (12)

    式(12)中:sn為隨機(jī)整數(shù),滿足sn∈{1,2,…,n};對(duì)每個(gè)解個(gè)體需要產(chǎn)生一個(gè)sn;Cr為交叉概率因子,Cr∈[0,1]。

    2.5 選擇過程

    差分進(jìn)化算法釆用的是一對(duì)一的貪婪選擇。這里使用自適應(yīng)懲罰函數(shù)法[11]對(duì)目標(biāo)個(gè)體和試驗(yàn)個(gè)體的目標(biāo)函數(shù)值及違反約束程度值進(jìn)行處理,求得目標(biāo)個(gè)體和試驗(yàn)個(gè)體的修正目標(biāo)函數(shù)值分別為F(Xi)及F(Ui),計(jì)算方法為:

    3) 計(jì)算目標(biāo)個(gè)體和試驗(yàn)個(gè)體的距離分別為d(Xi)及d(Ui)。

    4) 計(jì)算目標(biāo)個(gè)體和試驗(yàn)個(gè)體的懲罰項(xiàng)分別為p(Xi)及p(Ui)。

    5) 按照式(13)和式(14)計(jì)算目標(biāo)個(gè)體及試驗(yàn)個(gè)體的修正目標(biāo)函數(shù)值為

    則選擇過程為

    (15)

    式(15)中:Xinext為下一代種群中第i個(gè)個(gè)體。

    2.6 終止條件

    當(dāng)問題達(dá)到所需精度時(shí)即可停止,這樣能明顯縮短計(jì)算時(shí)間。同時(shí),進(jìn)化代數(shù)T為另一停止準(zhǔn)則,當(dāng)所設(shè)精度不能達(dá)到所需時(shí),只要進(jìn)化代數(shù)超過給定的最大進(jìn)化代數(shù),進(jìn)化即停止。

    3 程序設(shè)計(jì)

    基于Visual Studio 2013開發(fā)平臺(tái),采用C#語言編寫“最后分艙裝載”計(jì)算程序,其流程圖見圖1。

    4 算例及分析

    以某25萬噸級(jí)礦砂船為例,選擇1艙和9艙進(jìn)行“最后分艙裝載”的實(shí)例計(jì)算。船舶主尺度為:垂線間長(zhǎng)LPP=319.5 m;型寬B=57 m;型深D=25 m;設(shè)計(jì)吃水T=18 m;全船共有9個(gè)貨艙。

    圖1 “最后分艙裝載”計(jì)算流程圖

    4.1 算例一

    貨物密度ρRHO=1 t/m3,分艙裝載前后平均吃水變化為ΔT=1.1 m。

    該船“最后分艙裝載”前后的浮態(tài)和分艙前各貨艙的裝載量分別見表1及表2。使用傳統(tǒng)的“縱傾力矩法”和文中方法計(jì)算1艙及9艙的分貨量見表3。將此結(jié)果輸入到Onboard-Napa裝載計(jì)算機(jī)中得到平均吃水、吃水差和由此計(jì)算的誤差見表4。

    表1 船舶初始浮態(tài)和預(yù)期最終浮態(tài) m

    表2 分艙前各貨艙的裝載量

    表3 傳統(tǒng)方法、文中方法的計(jì)算結(jié)果 t

    表4 傳統(tǒng)方法與文中方法計(jì)算結(jié)果對(duì)比

    利用傳統(tǒng)“縱傾力矩法”得到的平均吃水和吃水差分別為3.86 m及-4.07 m,平均吃水和吃水差的誤差分別為1.026%及8.33%;使用文中方法得到的平均吃水和吃水差分別為3.90 m及-4.44 m,平均吃水和吃水差的誤差均為0%。

    4.2 算例二

    貨物密度ρRHO=1.584 t/m3,分艙裝載前后平均吃水變化為ΔT=0.23 m。

    該船“最后分艙裝載”前后的浮態(tài)和分艙前各貨艙的裝載量分別見表5及表6。使用傳統(tǒng)的“縱傾力矩法”和文中方法計(jì)算1艙和9艙的分貨量見表7。將此結(jié)果輸入到Onboard-Napa裝載計(jì)算機(jī)中得到平均吃水、吃水差和由此計(jì)算的誤差見表8。

    表5 船舶初始浮態(tài)和預(yù)期最終浮態(tài) m

    表6 分艙前各貨艙的裝載 t

    表7 傳統(tǒng)方法、文中方法計(jì)算的結(jié)果

    表8 傳統(tǒng)方法與文中方法計(jì)算結(jié)果對(duì)比

    利用傳統(tǒng)“縱傾力矩法”得到的平均吃水和吃水差分別為18.64 m及-0.45 m,平均吃水和吃水差的誤差分別為0%及2.17%;使用文中方法得到的平均吃水和吃水差分別為18.64 m及-0.46 m,平均吃水和吃水差的誤差均為0%。

    由上述2個(gè)算例可知:傳統(tǒng)“縱傾力矩法”僅在平均吃水變化小時(shí)誤差較小,在平均吃水變化較大時(shí)誤差會(huì)變大;文中方法不僅在平均吃水變化小時(shí)誤差較小,且在平均吃水變化較大時(shí)誤差也較小。

    5 結(jié)束語

    根據(jù)散貨船“最后分艙裝載”的特點(diǎn)建立散貨船“最后分艙裝載”的優(yōu)化模型;基于散貨船的艙容表,使用差分進(jìn)化算法對(duì)該模型進(jìn)行求解。以某25萬噸級(jí)礦砂船為例進(jìn)行實(shí)例計(jì)算,計(jì)算結(jié)果表明,文中方法不僅適合使用少量貨物進(jìn)行“最后分艙裝載”計(jì)算,而且適合使用大量貨物進(jìn)行“最后分艙裝載”計(jì)算,可有效減小誤差,使船舶達(dá)到預(yù)期的平均吃水和吃水差值,具有通用性。該方法已成功應(yīng)用到自主研發(fā)的裝載計(jì)算中,具有一定的工程實(shí)用價(jià)值。

    這里僅對(duì)散貨船的艏艙和艉艙進(jìn)行“最后分艙裝載”計(jì)算,在以后的工作中將研究使用多艙進(jìn)行“最后分艙裝載”計(jì)算。

    [1] 東昉. 散裝船吃水差調(diào)整的一種新方法[J].大連海事大學(xué)學(xué)報(bào),1996,22(1):15-19.

    [2] 欒法敏. 散貨船“最后分艙裝載”的一種計(jì)算方法[J].中國航海,2013,36(2):135-137.

    [3] 邢向輝. 散貨船裝載計(jì)算機(jī)數(shù)學(xué)模型及應(yīng)用的研究[D].大連:大連海事大學(xué),2000.

    [4] 雷坤,林焰,紀(jì)卓尚. 基于內(nèi)點(diǎn)法的駁船接載配載方案優(yōu)化模型[J].中國造船, 2010, 51(2):149-154.

    [5] 杜嘉立,姜華. 船舶原理[M].大連:大連海事大學(xué)出版社, 2011:13-17.

    [6] STOM R, PRICE K. Differential Evolution-A Simple and Efficient Adaptive Scheme for Global Optimization over Continuous Spaces[C]. International Computer Science Institute, Berkeley: ICSI, 1995.

    [7] STOM R,PRICE K. Differential Evolution-A Simple and Efficient Heuristic for Global Optimization over Continuous Spaces [J]. Journal of Global Optimization. 1997, 11(4):341-359.

    [8] STOM R, PRICE K, LAMPINEN J A. Differential Evolution: A Practical Approach to Global Optimization[M]. Berlin, Germany: Springer, 2006:1-6.

    [9] 徐斌. 基于差分進(jìn)化算法的多目標(biāo)優(yōu)化方法[D]. 上海:華東理工大學(xué),2013:24-32.

    [10] 王勇,蔡自興,周育人,等. 約束優(yōu)化進(jìn)化算法[J]. 軟件學(xué)報(bào), 2009,20(1):11-29.

    [11] TESSEMA B, YEN G G. A Self Adaptive Penalty Function Based Algorithm for Constrained Optimization [C]. IEEE Congress on Evolutionary Computation, 2006: 246-253.

    Calculation Method of “Final Trimming” for Bulk Carriers Based on Differential Evolution Algorithm

    LIUChunlei,SUNXiaofeng,WUHaibo,YINYong

    (Key Laboratory of Marine Dynamic Simulation & Control for Ministry of Transport,Dalian Maritime University, Dalian 116026, China)

    1000-4653(2016)04-0052-04

    U693.6

    A

    2016-07-12

    國家高技術(shù)研究發(fā)展計(jì)劃(“八六三”計(jì)劃)課題(2015AA016404)

    劉春雷(1987—),男,河北保定人,博士生,研究方向?yàn)樯⒇洿悄芑漭d。E-mail:18940930042@qq.com 孫霄峰(1978—),男,山東濟(jì)寧人,副教授,博士生,研究方向?yàn)榇安倏v性及航海動(dòng)態(tài)仿真。E-mail:xfsun_dlmu@163.com

    猜你喜歡
    貨艙散貨船差分
    數(shù)列與差分
    20艘!今年好望角型散貨船拆解量遠(yuǎn)超去年同期
    中國船檢(2020年3期)2020-11-27 09:14:20
    空中客車A330 飛機(jī)貨艙通風(fēng)系統(tǒng)
    基于VECM的國際干散貨船投資
    中國航海(2019年2期)2019-07-24 08:26:58
    5萬噸長(zhǎng)貨艙多用途船舶開發(fā)設(shè)計(jì)
    民用飛機(jī)貨艙煙霧探測(cè)系統(tǒng)研究
    民用飛機(jī)貨艙地板設(shè)計(jì)研究
    科技視界(2019年3期)2019-04-20 02:23:20
    79800 DWT雙燃料散貨船LNG儲(chǔ)氣罐位置確定
    船海工程(2015年4期)2016-01-05 15:53:50
    5 000 t散貨船清潔推進(jìn)系統(tǒng)配置論證與控制設(shè)計(jì)
    基于差分隱私的大數(shù)據(jù)隱私保護(hù)
    国产精品人妻久久久影院| 国产精品99久久久久久久久| 熟女电影av网| 久久精品国产亚洲av香蕉五月| 久久精品综合一区二区三区| 国产淫片久久久久久久久| 亚洲图色成人| 欧美不卡视频在线免费观看| 26uuu在线亚洲综合色| 丰满乱子伦码专区| 99热全是精品| 免费一级毛片在线播放高清视频| 我的老师免费观看完整版| 国产精品一及| 国产亚洲精品久久久com| 国产久久久一区二区三区| 国产精品99久久久久久久久| 午夜a级毛片| 亚洲熟妇中文字幕五十中出| 最近视频中文字幕2019在线8| 日本av手机在线免费观看| 我要看日韩黄色一级片| 亚洲av一区综合| 听说在线观看完整版免费高清| 日韩欧美三级三区| 国产av在哪里看| 国内少妇人妻偷人精品xxx网站| 国产极品精品免费视频能看的| 国产成人freesex在线| 亚洲色图av天堂| 中文资源天堂在线| 大香蕉久久网| 国产精品99久久久久久久久| 精品不卡国产一区二区三区| 久久久成人免费电影| 亚洲aⅴ乱码一区二区在线播放| 99热6这里只有精品| 国产伦精品一区二区三区四那| 中文亚洲av片在线观看爽| 在线观看一区二区三区| 丝袜美腿在线中文| 日韩欧美国产在线观看| 亚洲成人av在线免费| 男女做爰动态图高潮gif福利片| 99九九线精品视频在线观看视频| 桃色一区二区三区在线观看| 亚洲性久久影院| 日产精品乱码卡一卡2卡三| 久久久久久久午夜电影| 国产真实乱freesex| 日韩视频在线欧美| 国产高潮美女av| 欧美精品国产亚洲| 午夜爱爱视频在线播放| 欧美+日韩+精品| 欧美日韩在线观看h| 18禁在线播放成人免费| 成年av动漫网址| 亚洲色图av天堂| 国产一区二区三区av在线 | a级毛色黄片| 乱系列少妇在线播放| 精品人妻视频免费看| 久久精品国产自在天天线| 我的女老师完整版在线观看| 丰满乱子伦码专区| 国产精品爽爽va在线观看网站| 婷婷色综合大香蕉| 日韩中字成人| 色噜噜av男人的天堂激情| 99久久精品热视频| 日韩欧美三级三区| 国产精品av视频在线免费观看| www日本黄色视频网| 久久国产乱子免费精品| 精华霜和精华液先用哪个| 久久精品久久久久久噜噜老黄 | 91麻豆精品激情在线观看国产| 亚洲av熟女| 最近的中文字幕免费完整| 麻豆一二三区av精品| 免费观看在线日韩| 国产成年人精品一区二区| 一区二区三区免费毛片| 成人鲁丝片一二三区免费| 国产av麻豆久久久久久久| 看黄色毛片网站| 最好的美女福利视频网| 欧美bdsm另类| 最近视频中文字幕2019在线8| 一本一本综合久久| 菩萨蛮人人尽说江南好唐韦庄 | 色哟哟哟哟哟哟| 91av网一区二区| 91久久精品国产一区二区三区| 亚洲国产精品成人综合色| 精品一区二区三区视频在线| 国产69精品久久久久777片| 51国产日韩欧美| 成人二区视频| 久久久国产成人精品二区| 女的被弄到高潮叫床怎么办| 亚洲国产精品国产精品| 美女被艹到高潮喷水动态| 看黄色毛片网站| 少妇高潮的动态图| 国产精品一区二区在线观看99 | 老司机影院成人| 美女被艹到高潮喷水动态| 亚洲精品亚洲一区二区| 欧美zozozo另类| 国产乱人偷精品视频| 国产成人精品久久久久久| 亚洲内射少妇av| 99热全是精品| 插逼视频在线观看| 中文字幕熟女人妻在线| 直男gayav资源| 在线观看免费视频日本深夜| 国产一区二区亚洲精品在线观看| 能在线免费观看的黄片| 久久精品国产亚洲av涩爱 | 91精品国产九色| av福利片在线观看| or卡值多少钱| 99热这里只有精品一区| 成人毛片a级毛片在线播放| 国产一区二区三区在线臀色熟女| 日本撒尿小便嘘嘘汇集6| 欧美日韩国产亚洲二区| 亚洲欧美精品自产自拍| 麻豆精品久久久久久蜜桃| 成熟少妇高潮喷水视频| 桃色一区二区三区在线观看| 国产av一区在线观看免费| 午夜福利成人在线免费观看| 欧美成人精品欧美一级黄| 国产精品麻豆人妻色哟哟久久 | 丰满人妻一区二区三区视频av| 丰满乱子伦码专区| 日本免费a在线| 我要搜黄色片| 午夜福利在线观看吧| 啦啦啦啦在线视频资源| 在线观看66精品国产| 久久中文看片网| 一进一出抽搐动态| 国产精品.久久久| 国产成人91sexporn| 老司机影院成人| 国产精品人妻久久久久久| 蜜臀久久99精品久久宅男| 一区二区三区四区激情视频 | 国产精品一及| 日韩一区二区三区影片| www.av在线官网国产| 亚洲av男天堂| 精品无人区乱码1区二区| 中文字幕av在线有码专区| 国产三级在线视频| 少妇裸体淫交视频免费看高清| 性欧美人与动物交配| 搞女人的毛片| 免费不卡的大黄色大毛片视频在线观看 | 国内少妇人妻偷人精品xxx网站| 综合色av麻豆| 亚洲第一电影网av| 身体一侧抽搐| 校园人妻丝袜中文字幕| 亚洲欧美日韩高清在线视频| 国产精品福利在线免费观看| videossex国产| 成人综合一区亚洲| 国产一区二区在线观看日韩| 久久久色成人| 亚洲色图av天堂| 国产精品人妻久久久影院| 桃色一区二区三区在线观看| 人人妻人人澡人人爽人人夜夜 | 久久精品91蜜桃| 有码 亚洲区| 老司机福利观看| 成人特级黄色片久久久久久久| 哪里可以看免费的av片| 亚洲一区高清亚洲精品| av福利片在线观看| 中国国产av一级| 又爽又黄a免费视频| 国产乱人视频| 欧美xxxx黑人xx丫x性爽| 成人美女网站在线观看视频| 色噜噜av男人的天堂激情| 久久久色成人| 亚洲经典国产精华液单| 国产一区二区三区在线臀色熟女| 高清在线视频一区二区三区 | 国产精品福利在线免费观看| 欧美色视频一区免费| 久久久久免费精品人妻一区二区| 国产私拍福利视频在线观看| 99热全是精品| 国产一级毛片七仙女欲春2| 一级二级三级毛片免费看| 国产在视频线在精品| 蜜臀久久99精品久久宅男| 性色avwww在线观看| 中文字幕av在线有码专区| 99久久精品热视频| 国产精品伦人一区二区| 边亲边吃奶的免费视频| videossex国产| 18+在线观看网站| 色综合亚洲欧美另类图片| 午夜精品在线福利| 国产精品国产高清国产av| 久久久欧美国产精品| 此物有八面人人有两片| 欧美成人一区二区免费高清观看| 欧美变态另类bdsm刘玥| 精品午夜福利在线看| 国产精品久久久久久久久免| 2022亚洲国产成人精品| 禁无遮挡网站| 伦理电影大哥的女人| 精品久久久久久久久av| 国产 一区 欧美 日韩| 亚洲欧美精品专区久久| 乱系列少妇在线播放| 网址你懂的国产日韩在线| 麻豆国产97在线/欧美| 国产91av在线免费观看| 黄片无遮挡物在线观看| 一进一出抽搐动态| 国产一级毛片七仙女欲春2| 狠狠狠狠99中文字幕| 男女那种视频在线观看| 男女做爰动态图高潮gif福利片| 亚洲最大成人手机在线| 尾随美女入室| 精品久久国产蜜桃| av在线老鸭窝| 男女做爰动态图高潮gif福利片| 国产美女午夜福利| 日本成人三级电影网站| 69av精品久久久久久| 青春草视频在线免费观看| av国产免费在线观看| 麻豆成人午夜福利视频| 亚洲国产高清在线一区二区三| av在线天堂中文字幕| 最近最新中文字幕大全电影3| 国产av在哪里看| 一本久久精品| av天堂中文字幕网| 女人十人毛片免费观看3o分钟| 日韩大尺度精品在线看网址| 久久精品国产自在天天线| 午夜激情欧美在线| 国产精品一区二区性色av| 嫩草影院入口| 色综合色国产| 日本免费a在线| 久久久精品大字幕| 99热网站在线观看| 高清毛片免费看| 亚洲乱码一区二区免费版| 好男人在线观看高清免费视频| 亚洲欧美日韩东京热| 亚洲欧洲日产国产| 联通29元200g的流量卡| 国产高清视频在线观看网站| 日韩欧美国产在线观看| 在线观看一区二区三区| 色视频www国产| 国产成人精品久久久久久| 久久久精品欧美日韩精品| 桃色一区二区三区在线观看| 久久亚洲国产成人精品v| 日本免费一区二区三区高清不卡| 中文精品一卡2卡3卡4更新| 国产淫片久久久久久久久| 国产高清三级在线| 美女脱内裤让男人舔精品视频 | 美女cb高潮喷水在线观看| 国产乱人偷精品视频| av天堂中文字幕网| 国产高清视频在线观看网站| 成人性生交大片免费视频hd| 天堂av国产一区二区熟女人妻| 九九在线视频观看精品| 国产一区二区在线av高清观看| 久久精品国产亚洲网站| 亚洲精品亚洲一区二区| 国产成人aa在线观看| 中文亚洲av片在线观看爽| 中文字幕熟女人妻在线| kizo精华| 欧美性猛交╳xxx乱大交人| 大又大粗又爽又黄少妇毛片口| 亚洲欧洲国产日韩| 亚洲久久久久久中文字幕| 欧美日韩在线观看h| 午夜激情福利司机影院| 亚洲精品日韩在线中文字幕 | 最好的美女福利视频网| 国产精品美女特级片免费视频播放器| 国产极品精品免费视频能看的| 99国产精品一区二区蜜桃av| kizo精华| 欧美性猛交╳xxx乱大交人| 免费一级毛片在线播放高清视频| 亚洲18禁久久av| 国产av一区在线观看免费| 少妇人妻一区二区三区视频| 干丝袜人妻中文字幕| 国产精品久久久久久亚洲av鲁大| 美女cb高潮喷水在线观看| 欧美最新免费一区二区三区| 美女内射精品一级片tv| av福利片在线观看| 午夜老司机福利剧场| 欧美不卡视频在线免费观看| 少妇猛男粗大的猛烈进出视频 | 亚洲精品影视一区二区三区av| 欧美日韩国产亚洲二区| 超碰av人人做人人爽久久| 欧美变态另类bdsm刘玥| 狂野欧美激情性xxxx在线观看| 久久久精品大字幕| 国产成人精品婷婷| 啦啦啦观看免费观看视频高清| 国产高清激情床上av| 亚洲一区高清亚洲精品| 日本五十路高清| 国内精品宾馆在线| 国产一级毛片七仙女欲春2| 插阴视频在线观看视频| 变态另类成人亚洲欧美熟女| 国模一区二区三区四区视频| 人人妻人人澡人人爽人人夜夜 | av国产免费在线观看| 亚洲一区高清亚洲精品| 自拍偷自拍亚洲精品老妇| 最近的中文字幕免费完整| 久久久久久九九精品二区国产| 亚洲精品国产av成人精品| 免费av不卡在线播放| 精品免费久久久久久久清纯| 91狼人影院| 日韩强制内射视频| 中文字幕av在线有码专区| 久久久久久久午夜电影| 亚洲av中文字字幕乱码综合| 噜噜噜噜噜久久久久久91| 精品久久久久久久人妻蜜臀av| 97超视频在线观看视频| 久久久久久久亚洲中文字幕| 国产av不卡久久| 国产成人精品一,二区 | 精品久久久久久久人妻蜜臀av| 波多野结衣高清作品| 亚洲av成人精品一区久久| 日韩三级伦理在线观看| 国产色婷婷99| 一级毛片电影观看 | 天天一区二区日本电影三级| 丝袜美腿在线中文| 变态另类成人亚洲欧美熟女| 禁无遮挡网站| 亚洲中文字幕一区二区三区有码在线看| 中国美女看黄片| 婷婷精品国产亚洲av| 亚洲婷婷狠狠爱综合网| 亚洲电影在线观看av| 一边摸一边抽搐一进一小说| 人体艺术视频欧美日本| 日韩亚洲欧美综合| 在线观看av片永久免费下载| 国产日本99.免费观看| 波多野结衣高清无吗| 在线免费十八禁| 日本免费一区二区三区高清不卡| 国产成人aa在线观看| 欧美+亚洲+日韩+国产| 午夜激情福利司机影院| 91aial.com中文字幕在线观看| 国产黄a三级三级三级人| 国产成人精品一,二区 | 欧美高清成人免费视频www| 久久久久久久午夜电影| 久久久久久久久中文| 亚洲国产欧美人成| av在线亚洲专区| 国产精品女同一区二区软件| 国产毛片a区久久久久| 欧美在线一区亚洲| 插逼视频在线观看| 精品一区二区三区视频在线| 神马国产精品三级电影在线观看| 乱系列少妇在线播放| 亚洲精品自拍成人| 日本黄色片子视频| 2021天堂中文幕一二区在线观| 久久久久性生活片| 婷婷色综合大香蕉| 级片在线观看| 黑人高潮一二区| 久久久精品大字幕| 淫秽高清视频在线观看| 久久久久免费精品人妻一区二区| 少妇裸体淫交视频免费看高清| 国产午夜精品论理片| 日本免费a在线| 国产免费男女视频| 欧美激情在线99| 夫妻性生交免费视频一级片| 一本一本综合久久| 99国产极品粉嫩在线观看| 又爽又黄无遮挡网站| 26uuu在线亚洲综合色| 国产乱人偷精品视频| 黑人高潮一二区| 国产av麻豆久久久久久久| 日韩强制内射视频| www.av在线官网国产| 一区二区三区免费毛片| 我的老师免费观看完整版| 51国产日韩欧美| 欧美日本亚洲视频在线播放| 国产精品伦人一区二区| 亚洲综合色惰| 美女被艹到高潮喷水动态| 舔av片在线| 日韩国内少妇激情av| av视频在线观看入口| 黄色一级大片看看| 日韩高清综合在线| 国内精品美女久久久久久| 成人性生交大片免费视频hd| 亚洲自拍偷在线| 婷婷色av中文字幕| 亚洲美女视频黄频| 久久久精品大字幕| 99九九线精品视频在线观看视频| 国产精品乱码一区二三区的特点| av卡一久久| 天堂中文最新版在线下载 | 看十八女毛片水多多多| 亚洲av免费高清在线观看| 久久精品人妻少妇| 亚洲精品456在线播放app| 九九久久精品国产亚洲av麻豆| 别揉我奶头 嗯啊视频| 五月玫瑰六月丁香| 日日摸夜夜添夜夜添av毛片| 69av精品久久久久久| 三级男女做爰猛烈吃奶摸视频| 少妇裸体淫交视频免费看高清| 青春草视频在线免费观看| 91麻豆精品激情在线观看国产| 午夜激情欧美在线| 亚洲欧美精品自产自拍| 性插视频无遮挡在线免费观看| 欧美变态另类bdsm刘玥| 老熟妇乱子伦视频在线观看| 有码 亚洲区| 国产又黄又爽又无遮挡在线| 久久精品久久久久久噜噜老黄 | 一本精品99久久精品77| 少妇的逼好多水| 久久亚洲国产成人精品v| 欧美在线一区亚洲| 麻豆久久精品国产亚洲av| 久99久视频精品免费| 亚洲精品乱码久久久久久按摩| 国产亚洲欧美98| 看非洲黑人一级黄片| 十八禁国产超污无遮挡网站| 国产一区二区三区在线臀色熟女| 少妇熟女aⅴ在线视频| 欧美成人a在线观看| 国产免费一级a男人的天堂| 青春草国产在线视频 | 精品人妻熟女av久视频| 亚洲精品成人久久久久久| av黄色大香蕉| 狂野欧美激情性xxxx在线观看| 久久精品国产亚洲av涩爱 | 寂寞人妻少妇视频99o| 麻豆精品久久久久久蜜桃| 亚洲美女视频黄频| 美女被艹到高潮喷水动态| 嫩草影院新地址| 99视频精品全部免费 在线| 狂野欧美白嫩少妇大欣赏| 黄片wwwwww| 国产成人a区在线观看| 久久久精品欧美日韩精品| 久久鲁丝午夜福利片| 国产亚洲av嫩草精品影院| 久久久久久九九精品二区国产| 日本黄色片子视频| 久久中文看片网| 日韩在线高清观看一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 国产熟女欧美一区二区| 午夜激情福利司机影院| 国产精品av视频在线免费观看| 久久久久久久亚洲中文字幕| 超碰av人人做人人爽久久| 人体艺术视频欧美日本| 日韩一区二区三区影片| 日日撸夜夜添| 日韩成人伦理影院| 午夜福利在线观看吧| 久久精品人妻少妇| 国产久久久一区二区三区| 男女边吃奶边做爰视频| www.av在线官网国产| 麻豆国产av国片精品| 免费观看a级毛片全部| 久久人妻av系列| 成人特级av手机在线观看| 国产精品.久久久| 简卡轻食公司| 国产一区二区在线观看日韩| 青春草视频在线免费观看| 最近中文字幕高清免费大全6| 99视频精品全部免费 在线| 免费在线观看成人毛片| 高清日韩中文字幕在线| 91狼人影院| 美女cb高潮喷水在线观看| 久久人人精品亚洲av| 亚洲图色成人| 久久久精品欧美日韩精品| 看十八女毛片水多多多| 51国产日韩欧美| 男女做爰动态图高潮gif福利片| 99热精品在线国产| 一区二区三区四区激情视频 | 日本黄色片子视频| 黑人高潮一二区| 伊人久久精品亚洲午夜| 午夜精品在线福利| 国产精品无大码| 精品久久久噜噜| 国产高清有码在线观看视频| 欧美精品一区二区大全| 欧美高清性xxxxhd video| 色综合站精品国产| 99久久九九国产精品国产免费| 久久久色成人| 在线观看一区二区三区| 亚洲三级黄色毛片| 久久韩国三级中文字幕| 亚洲精品日韩在线中文字幕 | 精品久久久久久久久av| 网址你懂的国产日韩在线| 热99re8久久精品国产| 美女国产视频在线观看| 91狼人影院| 在现免费观看毛片| 少妇丰满av| 在线播放国产精品三级| 99久久精品热视频| 国产色婷婷99| 国产久久久一区二区三区| 26uuu在线亚洲综合色| 国内久久婷婷六月综合欲色啪| 国产乱人视频| 亚洲人与动物交配视频| 哪个播放器可以免费观看大片| 精品久久久久久成人av| 日日撸夜夜添| 天堂av国产一区二区熟女人妻| 我要看日韩黄色一级片| 少妇熟女aⅴ在线视频| 国产精品女同一区二区软件| 欧美+日韩+精品| 毛片一级片免费看久久久久| 日本成人三级电影网站| 国产淫片久久久久久久久| 波多野结衣高清无吗| 三级毛片av免费| 大型黄色视频在线免费观看| 男人和女人高潮做爰伦理| 两个人的视频大全免费| 午夜激情欧美在线| 久久精品国产自在天天线| 国产色婷婷99| 深爱激情五月婷婷| 欧美不卡视频在线免费观看| 国产探花极品一区二区| 欧美xxxx黑人xx丫x性爽| 国产免费一级a男人的天堂| 亚洲成a人片在线一区二区| 国产 一区精品| 中文字幕人妻熟人妻熟丝袜美| 国产精品人妻久久久影院| 久久久欧美国产精品| 日日摸夜夜添夜夜爱| 国产一级毛片七仙女欲春2| 国产日本99.免费观看| 26uuu在线亚洲综合色| 亚洲av成人精品一区久久| 嫩草影院精品99| 国产成人freesex在线| 午夜免费激情av| 99久久久亚洲精品蜜臀av| 99热全是精品| 男女做爰动态图高潮gif福利片| 最近最新中文字幕大全电影3| 精品午夜福利在线看| 久久亚洲国产成人精品v|