• <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久久人妻综合| kizo精华| 久久99蜜桃精品久久| 日韩大片免费观看网站| 能在线免费看毛片的网站| 18禁动态无遮挡网站| 国产成人a区在线观看| 国产成人一区二区在线| 免费av不卡在线播放| 日韩人妻高清精品专区| 欧美日韩视频精品一区| 一区二区av电影网| www.色视频.com| 日韩大片免费观看网站| 一级片'在线观看视频| 男女边吃奶边做爰视频| 久久影院123| 免费大片18禁| 亚洲aⅴ乱码一区二区在线播放| 国产成人a∨麻豆精品| 一区二区三区乱码不卡18| 久久国内精品自在自线图片| 久久久久国产精品人妻一区二区| 特大巨黑吊av在线直播| 人妻一区二区av| 亚洲精品乱码久久久久久按摩| 午夜激情久久久久久久| 日韩 亚洲 欧美在线| 中文字幕久久专区| 日韩人妻高清精品专区| 久久综合国产亚洲精品| 美女视频免费永久观看网站| 五月开心婷婷网| 五月伊人婷婷丁香| 嫩草影院入口| 日本一本二区三区精品| 男的添女的下面高潮视频| 精品午夜福利在线看| 婷婷色综合www| 观看免费一级毛片| 国产成人精品久久久久久| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产最新在线播放| 午夜精品一区二区三区免费看| av.在线天堂| 亚洲综合色惰| 99久久精品国产国产毛片| tube8黄色片| 在线天堂最新版资源| 99久久中文字幕三级久久日本| 成人无遮挡网站| 亚洲精品乱码久久久v下载方式| 黄色一级大片看看| 少妇人妻 视频| 国产成人a∨麻豆精品| 久久精品国产亚洲av天美| 伦理电影大哥的女人| 久久精品国产自在天天线| 69人妻影院| 晚上一个人看的免费电影| 久久鲁丝午夜福利片| 免费观看a级毛片全部| 亚洲av一区综合| 亚洲精品中文字幕在线视频 | 国产成人免费观看mmmm| 好男人在线观看高清免费视频| 成人综合一区亚洲| 纵有疾风起免费观看全集完整版| 我的女老师完整版在线观看| 禁无遮挡网站| 亚洲人成网站高清观看| 成年女人看的毛片在线观看| 激情五月婷婷亚洲| 免费看日本二区| 久久鲁丝午夜福利片| 又大又黄又爽视频免费| 99久国产av精品国产电影| 成人美女网站在线观看视频| 在线观看av片永久免费下载| 亚洲av一区综合| 晚上一个人看的免费电影| 99久久九九国产精品国产免费| 日韩不卡一区二区三区视频在线| 男人和女人高潮做爰伦理| 国产成人精品一,二区| 身体一侧抽搐| 18禁在线播放成人免费| 亚洲久久久久久中文字幕| 岛国毛片在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 伦精品一区二区三区| 深夜a级毛片| av国产精品久久久久影院| 国产欧美日韩一区二区三区在线 | 亚洲欧美清纯卡通| 国产精品国产av在线观看| 2022亚洲国产成人精品| av在线播放精品| 亚洲久久久久久中文字幕| 色网站视频免费| 国产精品一及| 国产精品嫩草影院av在线观看| 联通29元200g的流量卡| 亚洲av国产av综合av卡| 91久久精品国产一区二区成人| 欧美高清性xxxxhd video| 久久久久久久久久人人人人人人| 男人添女人高潮全过程视频| 国产精品偷伦视频观看了| 色视频在线一区二区三区| 久久久久久久大尺度免费视频| 国产一区二区三区综合在线观看 | 国产成人a∨麻豆精品| 中文字幕人妻熟人妻熟丝袜美| 日本-黄色视频高清免费观看| 国产亚洲91精品色在线| 性色av一级| 国产成人精品福利久久| 各种免费的搞黄视频| 午夜爱爱视频在线播放| 日韩成人伦理影院| www.色视频.com| 色婷婷久久久亚洲欧美| 国产一区二区三区综合在线观看 | videossex国产| 99热这里只有是精品50| 91午夜精品亚洲一区二区三区| 欧美性猛交╳xxx乱大交人| 国产精品伦人一区二区| 亚洲欧美一区二区三区黑人 | 成人国产av品久久久| 欧美日韩综合久久久久久| 一级毛片 在线播放| 在线 av 中文字幕| 亚洲,欧美,日韩| 亚洲欧美成人综合另类久久久| 久久久久九九精品影院| 精品久久久久久久人妻蜜臀av| 欧美xxⅹ黑人| 日韩av不卡免费在线播放| 联通29元200g的流量卡| 亚洲无线观看免费| 欧美日韩综合久久久久久| 黄片无遮挡物在线观看| 97超视频在线观看视频| 亚洲最大成人av| 在线 av 中文字幕| av线在线观看网站| 尾随美女入室| 日本黄大片高清| 婷婷色综合大香蕉| 亚洲国产欧美人成| 男插女下体视频免费在线播放| 免费看av在线观看网站| 国产又色又爽无遮挡免| 三级国产精品欧美在线观看| 99热网站在线观看| 欧美成人a在线观看| 免费观看性生交大片5| 欧美 日韩 精品 国产| 国产午夜精品一二区理论片| 国产爱豆传媒在线观看| 亚洲,一卡二卡三卡| 久久久a久久爽久久v久久| 国产永久视频网站| 18禁在线无遮挡免费观看视频| 在线观看一区二区三区激情| 2021天堂中文幕一二区在线观| av卡一久久| 秋霞伦理黄片| 日韩制服骚丝袜av| 青春草国产在线视频| 国产黄片美女视频| 久久综合国产亚洲精品| videos熟女内射| 日韩人妻高清精品专区| 午夜老司机福利剧场| 99久久精品一区二区三区| 日本午夜av视频| 久久久久九九精品影院| 岛国毛片在线播放| 激情 狠狠 欧美| 久久久久久久久久久丰满| 午夜福利高清视频| 综合色丁香网| 一级毛片aaaaaa免费看小| 久久国产乱子免费精品| 国产精品99久久99久久久不卡 | 综合色av麻豆| 国产 一区 欧美 日韩| 2022亚洲国产成人精品| 日本一二三区视频观看| 国产毛片a区久久久久| 亚洲av.av天堂| 男插女下体视频免费在线播放| 国产综合懂色| 日本免费在线观看一区| 国产男女内射视频| 人妻系列 视频| 不卡视频在线观看欧美| 看非洲黑人一级黄片| 18禁裸乳无遮挡免费网站照片| 联通29元200g的流量卡| 成人毛片60女人毛片免费| 国产男女内射视频| 99热6这里只有精品| 可以在线观看毛片的网站| 欧美最新免费一区二区三区| 26uuu在线亚洲综合色| 亚洲精品影视一区二区三区av| 久久99热6这里只有精品| 免费看不卡的av| 麻豆久久精品国产亚洲av| xxx大片免费视频| 插阴视频在线观看视频| 亚洲国产成人一精品久久久| 国产伦理片在线播放av一区| 男女边摸边吃奶| 国产精品三级大全| 国产精品熟女久久久久浪| 中文字幕久久专区| 午夜福利高清视频| 狂野欧美激情性xxxx在线观看| 三级国产精品片| 搞女人的毛片| 欧美日韩精品成人综合77777| 国产精品女同一区二区软件| 国产69精品久久久久777片| 亚洲无线观看免费| 日本午夜av视频| 亚洲欧美日韩卡通动漫| 人妻 亚洲 视频| 久久久久精品久久久久真实原创| 51国产日韩欧美| 久久久久久久久久人人人人人人| 69人妻影院| 午夜免费观看性视频| 亚洲美女视频黄频| 韩国av在线不卡| 人体艺术视频欧美日本| 欧美一区二区亚洲| 久久久欧美国产精品| 最近最新中文字幕免费大全7| 一级黄片播放器| 在线a可以看的网站| 激情五月婷婷亚洲| 一级毛片aaaaaa免费看小| 亚洲怡红院男人天堂| 天堂俺去俺来也www色官网| 国产 一区精品| 春色校园在线视频观看| 成年版毛片免费区| 天美传媒精品一区二区| 18+在线观看网站| 亚洲精品乱久久久久久| 高清毛片免费看| 久久久午夜欧美精品| 亚洲最大成人手机在线| 国产成人精品福利久久| 亚洲精品久久午夜乱码| 欧美变态另类bdsm刘玥| 精品久久久久久久久av| 免费少妇av软件| 亚洲av成人精品一二三区| 天堂俺去俺来也www色官网| 人妻夜夜爽99麻豆av| 王馨瑶露胸无遮挡在线观看| 国产成人精品福利久久| 久久99热这里只频精品6学生| 80岁老熟妇乱子伦牲交| xxx大片免费视频| av线在线观看网站| 精品国产三级普通话版| 国语对白做爰xxxⅹ性视频网站| 天堂网av新在线| 边亲边吃奶的免费视频| 亚洲精品国产色婷婷电影| 青春草国产在线视频| 99久久精品一区二区三区| 日本与韩国留学比较| 大片免费播放器 马上看| 男女无遮挡免费网站观看| 国产黄色免费在线视频| 日韩av免费高清视频| 久久久久国产精品人妻一区二区| 久久久精品免费免费高清| 日韩人妻高清精品专区| 亚洲最大成人av| 欧美日韩亚洲高清精品| 三级男女做爰猛烈吃奶摸视频| 精品午夜福利在线看| 韩国av在线不卡| 免费看光身美女| 亚洲美女视频黄频| 高清欧美精品videossex| 可以在线观看毛片的网站| 亚洲精品亚洲一区二区| 亚洲精品久久午夜乱码| 亚洲久久久久久中文字幕| 天堂中文最新版在线下载 | 久久久久久久久久久丰满| 亚洲成人中文字幕在线播放| a级毛色黄片| 日韩人妻高清精品专区| 国产美女午夜福利| 国产爱豆传媒在线观看| 精品人妻熟女av久视频| 中文天堂在线官网| 黄色视频在线播放观看不卡| 国产永久视频网站| 亚洲欧美一区二区三区国产| 在线观看av片永久免费下载| 狂野欧美激情性bbbbbb| 久久热精品热| 黄色视频在线播放观看不卡| 22中文网久久字幕| 亚洲四区av| 国产熟女欧美一区二区| 男插女下体视频免费在线播放| 亚洲av二区三区四区| 久久精品久久精品一区二区三区| 天天一区二区日本电影三级| 日本爱情动作片www.在线观看| 国产午夜精品一二区理论片| 综合色av麻豆| 中国美白少妇内射xxxbb| 有码 亚洲区| 在线观看一区二区三区激情| 久久国产乱子免费精品| 在线免费十八禁| 久久久久久久午夜电影| 99久国产av精品国产电影| 亚洲最大成人中文| 亚洲欧美成人综合另类久久久| 简卡轻食公司| 人体艺术视频欧美日本| 亚洲欧洲日产国产| 日本三级黄在线观看| 成人亚洲欧美一区二区av| 久久久久久久午夜电影| 香蕉精品网在线| 欧美精品一区二区大全| 亚洲最大成人av| 亚洲天堂国产精品一区在线| 精品人妻偷拍中文字幕| 欧美激情在线99| 黄色欧美视频在线观看| 精华霜和精华液先用哪个| 热99国产精品久久久久久7| 欧美日本视频| 久久久久久伊人网av| 亚洲欧美日韩另类电影网站 | 久久久久久久国产电影| 国产高潮美女av| 色视频在线一区二区三区| 国产 一区 欧美 日韩| 少妇人妻一区二区三区视频| 婷婷色综合www| 91狼人影院| 亚洲av免费在线观看| 亚洲成人一二三区av| 日本黄大片高清| 青春草亚洲视频在线观看| 亚洲精品456在线播放app| 免费观看av网站的网址| 在线观看国产h片| tube8黄色片| 爱豆传媒免费全集在线观看| av在线老鸭窝| 国产男女内射视频| 欧美高清成人免费视频www| 国产高清国产精品国产三级 | 只有这里有精品99| 男人狂女人下面高潮的视频| 97人妻精品一区二区三区麻豆| 91精品国产九色| 亚洲精品成人久久久久久| 亚洲精品456在线播放app| 婷婷色综合大香蕉| 亚洲国产精品国产精品| 国产成年人精品一区二区| 欧美一级a爱片免费观看看| 久久久国产一区二区| 成人欧美大片| av在线天堂中文字幕| 国产伦理片在线播放av一区| 亚洲欧美日韩卡通动漫| 蜜桃亚洲精品一区二区三区| 又大又黄又爽视频免费| 乱码一卡2卡4卡精品| 亚洲精品色激情综合| 欧美一级a爱片免费观看看| 少妇的逼好多水| 波多野结衣巨乳人妻| 成人特级av手机在线观看| 在线播放无遮挡| 视频中文字幕在线观看| h日本视频在线播放| 一区二区三区免费毛片| 国产精品.久久久| 色播亚洲综合网| 在线a可以看的网站| 毛片一级片免费看久久久久| 久久久久久久久久久免费av| 搡老乐熟女国产| 亚洲精品一区蜜桃| 国内揄拍国产精品人妻在线| 中文字幕制服av| 插逼视频在线观看| 日韩成人av中文字幕在线观看| 中国美白少妇内射xxxbb| 中国三级夫妇交换| 观看免费一级毛片| 黄片wwwwww| 91精品国产九色| 97人妻精品一区二区三区麻豆| 建设人人有责人人尽责人人享有的 | 嫩草影院新地址| 国产欧美亚洲国产| 性插视频无遮挡在线免费观看| 亚洲一级一片aⅴ在线观看| 禁无遮挡网站| 亚洲精品乱码久久久久久按摩| 国产精品一二三区在线看| 成人毛片60女人毛片免费| 五月天丁香电影| 身体一侧抽搐| 男人和女人高潮做爰伦理| 国产免费福利视频在线观看| 男人爽女人下面视频在线观看| 寂寞人妻少妇视频99o| 男人和女人高潮做爰伦理| 成人毛片a级毛片在线播放| 日韩亚洲欧美综合| av网站免费在线观看视频| av女优亚洲男人天堂| 国产高潮美女av| 久久韩国三级中文字幕| 亚洲人成网站在线观看播放| 中国美白少妇内射xxxbb| 国产精品.久久久| 中国国产av一级| 老女人水多毛片| 男女边摸边吃奶| 日韩欧美精品免费久久| 国产又色又爽无遮挡免| 亚洲伊人久久精品综合| 99热这里只有是精品在线观看| 国产午夜福利久久久久久| 亚洲国产精品成人综合色| 另类亚洲欧美激情| 国内揄拍国产精品人妻在线| 99热6这里只有精品| 亚洲电影在线观看av| 亚洲av电影在线观看一区二区三区 | 天堂俺去俺来也www色官网| 国内精品美女久久久久久| 午夜福利视频精品| 国产精品不卡视频一区二区| 爱豆传媒免费全集在线观看| 全区人妻精品视频| 一级片'在线观看视频| 日韩电影二区| 色哟哟·www| 国产伦精品一区二区三区四那| 我要看日韩黄色一级片| 人妻制服诱惑在线中文字幕| 久久午夜福利片| 熟女电影av网| av女优亚洲男人天堂| 日产精品乱码卡一卡2卡三| 亚洲精品456在线播放app| 久久精品国产自在天天线| 99久久人妻综合| 日本午夜av视频| 青春草国产在线视频| 国产高清不卡午夜福利| 国产高潮美女av| 我要看日韩黄色一级片| 啦啦啦啦在线视频资源| 亚洲av欧美aⅴ国产| 一级毛片我不卡| 成年版毛片免费区| 国产精品人妻久久久影院| 亚洲不卡免费看| 婷婷色综合大香蕉| 美女内射精品一级片tv| 国产成人午夜福利电影在线观看| 亚洲av福利一区| 又黄又爽又刺激的免费视频.| 亚洲美女视频黄频| 中文字幕免费在线视频6| 成人黄色视频免费在线看| 国产精品蜜桃在线观看| 看非洲黑人一级黄片| 亚洲人成网站在线播| 精品一区二区免费观看| 亚洲精品乱码久久久v下载方式| av在线观看视频网站免费| 热re99久久精品国产66热6| 国产一区二区三区av在线| 久久99热6这里只有精品| 蜜臀久久99精品久久宅男| 欧美xxxx性猛交bbbb| 中文字幕av成人在线电影| 久久久久九九精品影院| 精品国产一区二区三区久久久樱花 | 国产在线男女| 少妇的逼好多水| 秋霞伦理黄片| 免费电影在线观看免费观看| 国产大屁股一区二区在线视频| 欧美bdsm另类| 欧美 日韩 精品 国产| 国产一区二区在线观看日韩| 极品教师在线视频| 97人妻精品一区二区三区麻豆| 嘟嘟电影网在线观看| 日韩中字成人| 日本三级黄在线观看| 精品一区二区免费观看| 青青草视频在线视频观看| 天堂网av新在线| 久久久久久久午夜电影| 内射极品少妇av片p| 亚洲,一卡二卡三卡| 黄色一级大片看看| 嫩草影院入口| 99精国产麻豆久久婷婷| 美女内射精品一级片tv| 国产高清不卡午夜福利| 亚洲综合精品二区| 高清午夜精品一区二区三区| 久久国内精品自在自线图片| 在线看a的网站| 丝袜脚勾引网站| 亚洲激情五月婷婷啪啪| a级毛色黄片| 少妇人妻一区二区三区视频| 18禁动态无遮挡网站| 少妇人妻 视频| 国产日韩欧美亚洲二区| 五月玫瑰六月丁香| 国产探花在线观看一区二区| 国产亚洲午夜精品一区二区久久 | 欧美zozozo另类| 欧美精品人与动牲交sv欧美| 中文字幕亚洲精品专区| 97在线视频观看| 国产男女超爽视频在线观看| 国产视频内射| 亚州av有码| 亚洲欧美日韩东京热| 好男人在线观看高清免费视频| 少妇被粗大猛烈的视频| 一级毛片aaaaaa免费看小| 毛片一级片免费看久久久久| freevideosex欧美| 免费黄频网站在线观看国产| 国产综合精华液| 大话2 男鬼变身卡| 中文字幕制服av| 欧美另类一区| 在线a可以看的网站| 性色av一级| 啦啦啦在线观看免费高清www| 国产精品一区二区三区四区免费观看| 日韩 亚洲 欧美在线| 免费av毛片视频| 欧美+日韩+精品| 亚洲国产精品成人久久小说| 国产伦精品一区二区三区视频9| 在线播放无遮挡| 亚洲四区av| 又大又黄又爽视频免费| 国产精品人妻久久久久久| 午夜福利视频1000在线观看| 欧美日韩精品成人综合77777| 亚洲av中文字字幕乱码综合| 一级毛片我不卡| 久久国产乱子免费精品| 午夜免费观看性视频| 日韩欧美精品免费久久| 成人黄色视频免费在线看| 秋霞伦理黄片| 性色avwww在线观看| av免费观看日本| 午夜视频国产福利| 亚洲国产最新在线播放| 久久国产乱子免费精品| 人妻 亚洲 视频| 最近手机中文字幕大全| 少妇被粗大猛烈的视频| 亚洲精华国产精华液的使用体验| 2021少妇久久久久久久久久久| 日韩伦理黄色片| kizo精华| 久久久久网色| 国内揄拍国产精品人妻在线| 尾随美女入室| 在线免费十八禁| av免费在线看不卡| 国产成人freesex在线| 97精品久久久久久久久久精品| 各种免费的搞黄视频| 女人久久www免费人成看片|