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

    散貨船自由浮態(tài)計(jì)算簡(jiǎn)化方法

    2017-02-09 05:23:17雷,勇,峰,鳳,
    關(guān)鍵詞:水線散貨船傾角

    劉 春 雷, 尹 勇, 孫 霄 峰, 張 秀 鳳, 神 和 龍

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

    散貨船自由浮態(tài)計(jì)算簡(jiǎn)化方法

    劉 春 雷, 尹 勇*, 孫 霄 峰, 張 秀 鳳, 神 和 龍

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

    根據(jù)散貨船特點(diǎn)探討了一種面向工程應(yīng)用的船舶自由浮態(tài)計(jì)算簡(jiǎn)化方法.在矩陣法的基礎(chǔ)上,將7個(gè)水線面要素簡(jiǎn)化到4個(gè),簡(jiǎn)化雅克比矩陣求解.采用全選主元高斯消去法求解線性方程組,以船舶浮態(tài)平衡方程組作為迭代終止條件,保證計(jì)算精度.以散貨船“DOLCE VITA”及“RUI AN CHENG”為例進(jìn)行了實(shí)例計(jì)算,對(duì)20個(gè)壓載艙進(jìn)行了20萬(wàn)次隨機(jī)裝載實(shí)驗(yàn),驗(yàn)證了算法在有橫傾下的收斂性和穩(wěn)定性;對(duì)5個(gè)貨艙進(jìn)行了20萬(wàn)次隨機(jī)裝載實(shí)驗(yàn),驗(yàn)證了算法在有縱傾下的收斂性和穩(wěn)定性;對(duì)典型載況進(jìn)行計(jì)算,和裝載手冊(cè)給定值相比,吃水差誤差都在0.01 m以下.結(jié)果表明:該算法只需計(jì)算任意傾斜水線面下的排水體積和浮心坐標(biāo),計(jì)算量較小,程序?qū)崿F(xiàn)簡(jiǎn)單,實(shí)時(shí)性、魯棒性較好,具有一定的工程實(shí)用價(jià)值.

    船舶工程;自由浮態(tài);散貨船;配載儀;吃水差

    0 引 言

    船舶浮態(tài)計(jì)算是指計(jì)算船舯吃水Tm、橫傾角θ及縱傾角φ,在船體坐標(biāo)系下多采用艏吃水Tf、艉吃水Ta及橫傾角θ表示.船舶浮態(tài)計(jì)算是保障船舶航行安全的重要前提,也是船舶穩(wěn)性計(jì)算、船舶剪力和彎矩計(jì)算、自動(dòng)配載及裝卸貨過(guò)程優(yōu)化的基礎(chǔ).

    目前船舶浮態(tài)計(jì)算方法大體可分為4類:第1類是依靠船舶靜水力資料計(jì)算的常規(guī)法[1-2];第2類是采用牛頓迭代的矩陣法[3-6];第3類計(jì)算方法可歸結(jié)為優(yōu)化法,主要有文獻(xiàn)[7-8]提出的非線性規(guī)劃法及文獻(xiàn)[9-11]提出的遺傳算法;第4類是有限元法,通過(guò)有限元分析計(jì)算船舶浮態(tài)[12-14].

    但是這4類算法存在一些不足:第1類常規(guī)法以船舶靜水力數(shù)據(jù)為基礎(chǔ),將縱傾和橫傾分開計(jì)算,適用于小縱傾計(jì)算,縱傾和橫傾較大時(shí)誤差會(huì)增大;第2類矩陣法在每一次迭代計(jì)算時(shí)都要計(jì)算包含水線面面積、漂心、慣性矩、排水體積及浮心等多項(xiàng)要素的雅克比矩陣,需實(shí)時(shí)計(jì)算船舶傾斜水線面,計(jì)算量較大,程序處理困難;第3類優(yōu)化法只需要計(jì)算排水體積和浮心,避免了雅克比矩陣的計(jì)算,減少了計(jì)算量,但是迭代次數(shù)較多,收斂速度有時(shí)比較緩慢,影響程序?qū)崟r(shí)性,不太適用于工程應(yīng)用;第4類算法采用有限元法,需建立船舶有限元模型,計(jì)算量大,實(shí)時(shí)性不好.

    針對(duì)以上不足,本文根據(jù)散貨船方形系數(shù)較大的特點(diǎn),探討一種面向工程應(yīng)用的浮態(tài)迭代算法,在求解雅克比矩陣時(shí),用過(guò)船艉吃水處且平行于基平面的水線面代替傾斜水線面在基平面的投影面,將求解雅克比矩陣時(shí)的7個(gè)水線面要素簡(jiǎn)化到4個(gè),簡(jiǎn)化雅克比矩陣求解,程序只需計(jì)算任意傾斜水線面下的排水體積及浮心坐標(biāo),以使編程實(shí)現(xiàn)簡(jiǎn)單,穩(wěn)定可靠.

    1 船舶自由浮態(tài)方程組求解

    本文采用固定在船上的Oxyz坐標(biāo)系統(tǒng):原點(diǎn)O為基平面、中站面和中線面的交點(diǎn).x軸為基平面和中線面的交線,指向船艏為正.y軸是基平面和中站面的交線,指向右舷為正.z軸是中線面和中站面的交線,向上為正.船舶任意浮態(tài)平衡方程組為[5]

    (1)

    式中:ρ為海水密度;V為排水體積;Δ為船舶總裝載量;xb、yb、zb為船舶浮心縱向、橫向及垂向坐標(biāo);xg、yg、zg為船舶重心縱向、橫向及垂向坐標(biāo);θ為船舶橫傾角;φ為船舶縱傾角.

    在求解時(shí),一般通過(guò)牛頓迭代求解非線性方程組[5],其雅克比矩陣為

    (2)

    在求解雅克比矩陣時(shí),其難點(diǎn)是7個(gè)水線面要素S、xf、yf、zf、Ixf、Iyf、Ixyf的求解,需要傾斜水線面和船體實(shí)時(shí)求交計(jì)算水線面的輪廓數(shù)據(jù),文獻(xiàn)[6]提出一種獲得水線面輪廓的排序方法.本文根據(jù)散貨船的特點(diǎn)提出一種簡(jiǎn)化方法,將7個(gè)水線面要素簡(jiǎn)化到4個(gè),且不需要實(shí)時(shí)計(jì)算.如圖1所示,散貨船的方形系數(shù)較大,可用過(guò)船艉吃水Ta處的平行于基平面的水線面Aw1近似代替傾斜水線面在基平面的投影面Aw.

    根據(jù)水平水線面S1特點(diǎn)可得到如下性質(zhì):yf=0,Ixy=0,zf=Tm-dtanφ≈Tm,d為漂心縱向坐標(biāo)距船舯的距離,另外根據(jù)平行移軸原理則可得到Ix、Iy、Ixy與Ixf、Iyf、Ixyf的對(duì)應(yīng)轉(zhuǎn)換關(guān)系(Ix、Iy、Ixy分別為水線面對(duì)x、y軸的慣性矩和慣性積):

    (3)

    圖1 水線面投影

    則原雅克比矩陣可簡(jiǎn)化為

    (4)

    由式(4)可知,求解方程組只需要計(jì)算4個(gè)水線面要素S、xf、Ix及Iyf.在編寫程序時(shí)可離線計(jì)算出各個(gè)吃水下的水線面參數(shù)值,建立離線數(shù)據(jù)庫(kù),根據(jù)船艉吃水進(jìn)行離線插值,提高程序計(jì)算效率.在建立離線水線面數(shù)據(jù)庫(kù)時(shí)可由船舶設(shè)計(jì)部門提供的船舶靜水力表獲得,也可采用格林公式進(jìn)行計(jì)算,計(jì)算公式如式(5)所示:

    (5)

    其中My是水線面對(duì)y軸的靜矩.

    采用全選主元高斯消去法求解三元一次線性方程組,以船舶任意浮態(tài)平衡條件作為迭代終止判斷條件,當(dāng)滿足船舶平衡條件時(shí)輸出浮態(tài)計(jì)算結(jié)果:

    (6)

    式中:ε1、ε2、ε3為迭代精度控制,本文取ε1=5,ε2=0.001,ε3=0.001.

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

    首先采用常規(guī)法計(jì)算船舶浮態(tài)初始值,然后采用迭代法進(jìn)行迭代,直到滿足船舶浮態(tài)平衡條件,程序設(shè)計(jì)流程圖如圖2所示,其具體步驟如下:

    圖2 程序設(shè)計(jì)流程圖

    步驟1 計(jì)算船舶總裝載量Δ及重心坐標(biāo)(xg,yg,zg).

    步驟2 常規(guī)法計(jì)算船舶Tm、tanθ、tanφ.

    步驟3 計(jì)算船舶傾斜水線面方程xtanφ+ytanθ-z+Tm=0,傾斜水線面和船舶外殼求交計(jì)算水線面下排水體積V和浮心坐標(biāo)(xb,yb,zb).

    步驟4 根據(jù)式(6)判斷是否滿足平衡條件,即判斷總裝載量和排水量是否相等、重心和浮心是否在同一條垂線上,如果不滿足,則轉(zhuǎn)步驟6.

    步驟5 浮態(tài)程序計(jì)算結(jié)束,輸出當(dāng)前浮態(tài)結(jié)果Tm、tanθ、tanφ.

    步驟6 根據(jù)Ta插值靜水力表計(jì)算Ix、Iyf、xf、S.

    步驟7 計(jì)算雅克比矩陣.

    步驟8 全選主元高斯消去法解三元一次方程組得到δTm、δtanθ和δtanφ.

    步驟9 計(jì)算船舶Tm、tanθ、tanφ:

    (7)

    轉(zhuǎn)步驟3.

    3 算 例

    本文以38 000 t散貨船“DOLCE VITA”和“RUI AN CHENG”為例進(jìn)行算法收斂性和準(zhǔn)確性測(cè)試:(1)對(duì)其20個(gè)壓載艙進(jìn)行20萬(wàn)次隨機(jī)裝載實(shí)驗(yàn)測(cè)試算法在有橫傾情況下的收斂性和穩(wěn)定性;(2)對(duì)其5個(gè)貨艙進(jìn)行20萬(wàn)次隨機(jī)裝載實(shí)驗(yàn)測(cè)試算法在船舶有吃水差下的收斂性和穩(wěn)定性;(3)對(duì)幾種典型載況進(jìn)行計(jì)算,和裝載手冊(cè)給定值相比較驗(yàn)證算法準(zhǔn)確性.

    算例1 20個(gè)壓載艙20萬(wàn)次隨機(jī)裝載實(shí)驗(yàn)

    在初始裝載狀況為空載的情況下進(jìn)行了20萬(wàn)次隨機(jī)裝載實(shí)驗(yàn),以每個(gè)壓載艙的最大艙容為限制,在20個(gè)壓載艙隨機(jī)產(chǎn)生裝貨量進(jìn)行浮態(tài)計(jì)算.計(jì)算結(jié)果如圖3所示,船舶橫傾角計(jì)算結(jié)果的分布范圍從左傾12.160°到右傾12.856°,最大迭代次數(shù)為5次,最小迭代次數(shù)為1次,全部20萬(wàn)次裝載方案的平均迭代次數(shù)為2.29次.1次迭代即可收斂的載況為10 737個(gè),2、3、4、5次迭代收斂的載況分別為124 402、60 813、4 013及35個(gè),迭代失敗的次數(shù)為0.

    圖4所示為各橫傾角下的平均迭代次數(shù),從圖中可看出隨著橫傾角絕對(duì)值的減小,迭代次數(shù)逐漸減少:船舶橫傾角位于-13°~-12°時(shí)平均迭代次數(shù)為5次;橫傾角位于0°~1°時(shí)平均迭代次數(shù)為1.9次;橫傾角位于11°~12°時(shí)平均迭代次數(shù)為5次.

    圖3 各迭代次數(shù)下的載況數(shù)(橫傾)

    圖4 各橫傾角下的平均迭代次數(shù)

    從實(shí)驗(yàn)結(jié)果可以看出算法在橫傾角較小時(shí)一般需要1~2次迭代即可收斂,當(dāng)橫傾角較大時(shí)一般不超過(guò)5次迭代即可收斂,驗(yàn)證了算法在船舶存在橫傾時(shí)的收斂性和穩(wěn)定性.

    算例2 5個(gè)貨艙20萬(wàn)次隨機(jī)裝載實(shí)驗(yàn)

    在初始裝載狀況為空載的情況下進(jìn)行了20萬(wàn)次隨機(jī)裝載實(shí)驗(yàn),以每個(gè)貨艙的最大艙容為限制,在5個(gè)貨艙隨機(jī)產(chǎn)生裝貨量進(jìn)行浮態(tài)計(jì)算.計(jì)算結(jié)果如圖5所示,船舶吃水差計(jì)算結(jié)果從-13.040 m到13.508 m,最大迭代次數(shù)為6次,最小迭代次數(shù)為1次,全部20萬(wàn)次裝載方案的平均迭代次數(shù)為1.74次.1次迭代即可收斂的載況為89 949個(gè),2、3、4、5、6次迭代收斂的載況分別為80 759、22 143、5 430、1 567及152個(gè),迭代失敗的次數(shù)為0.

    圖6所示為各吃水差下的平均迭代次數(shù),從圖中可看出隨著吃水差絕對(duì)值的減小,迭代次數(shù)逐漸減?。捍俺运钗挥?14~-13 m時(shí)平均迭代次數(shù)為4次;吃水差位于0~1 m時(shí)平均迭代次數(shù)為1.27次;吃水差位于12~13 m時(shí)平均迭代次數(shù)為5次.

    圖5 各迭代次數(shù)下的載況數(shù)(縱傾)

    圖6 各吃水差下平均迭代次數(shù)

    從實(shí)驗(yàn)結(jié)果可以看出算法在吃水差較小時(shí)一般需要1~2次迭代即可收斂,當(dāng)吃水差較大時(shí)一般不超過(guò)6次迭代即可收斂,驗(yàn)證了算法在船舶存在縱傾時(shí)的收斂性和穩(wěn)定性.

    算例3 算法準(zhǔn)確性驗(yàn)證

    為驗(yàn)證算法準(zhǔn)確性,本文對(duì) “DOLCE VITA”4種典型載況進(jìn)行了實(shí)例計(jì)算,如表1所示.表2為“RUI AN CHENG”8種典型載況的浮態(tài)計(jì)算結(jié)果.

    分析表1、2中數(shù)據(jù)可知:(1)吃水差誤差都在0.010 m以下,最大誤差為0.009 m,可見采用本文方法計(jì)算的船舶浮態(tài)計(jì)算結(jié)果和裝載手冊(cè)給定值誤差較小,具有較高的計(jì)算精度.(2)船舶空船載況下,船舶吃水差較大,本文方法需要進(jìn)行4次迭代,其他載況基本需要1次迭代即可獲得較精確的結(jié)果.

    表1 散貨船“DOLCE VITA”浮態(tài)計(jì)算結(jié)果

    表2 散貨船“RUI AN CHENG”浮態(tài)計(jì)算結(jié)果

    可見,本文算法以船舶浮態(tài)方程組作為迭代終止條件,只有滿足船舶平衡條件時(shí)才輸出計(jì)算結(jié)果,保證了計(jì)算精度,和裝載手冊(cè)給定值相比誤差較小,驗(yàn)證了算法的可行性與準(zhǔn)確性.

    基于文中所述浮態(tài)算法,大連海事大學(xué)航海研究所自主研發(fā)出一套基于STL模型散貨船裝載計(jì)算機(jī)程序,如圖7所示,目前該軟件正在申請(qǐng)CCS及DNV船級(jí)社的型式認(rèn)可.

    圖7 散貨船配載儀

    4 結(jié) 論

    (1)針對(duì)散貨船的特點(diǎn),在矩陣法的基礎(chǔ)上探討了一種面向工程應(yīng)用的船舶浮態(tài)迭代算法.相比于常規(guī)法,本文算法適用于任意浮態(tài);相比于矩陣法減少了水線面面積、漂心、慣性矩等要素的實(shí)時(shí)求解,減少了計(jì)算量;相比于優(yōu)化法及有限元法減少了迭代次數(shù).

    (2)算法以船舶浮態(tài)平衡方程組作為迭代終止條件,保證了計(jì)算精度.

    (3)算法簡(jiǎn)單,容易編程實(shí)現(xiàn),依賴于離線靜水力表,只需要計(jì)算出任意傾斜水線面下的排水體積和浮心,具有實(shí)時(shí)性、魯棒性及準(zhǔn)確性,工程實(shí)用價(jià)值較高.

    本文算法已經(jīng)成功應(yīng)用在自主研發(fā)的散貨船裝載計(jì)算機(jī)系統(tǒng)中,計(jì)算船型覆蓋了38 000、64 000 及250 000 t級(jí)別的主流散貨船型,未來(lái)將在此基礎(chǔ)上研究船舶自動(dòng)配載、貨物裝載過(guò)程中優(yōu)化等方面內(nèi)容.

    [1] 段興鋒,任鴻翔,東 昉. 縱傾狀態(tài)船舶浮態(tài)的研究與計(jì)算 [J]. 中國(guó)航海, 2015, 38(1):54-57.

    DUAN Xingfeng, REN Hongxiang, DONG Fang. Calculation of floating condition of a longitudinal inclined ship [J]. Navigation of China, 2015, 38(1):54-57. (in Chinese)

    [2] 劉春雷,尹 勇,孫霄峰,等. 船舶浮態(tài)計(jì)算的一種修正方法[J]. 大連海事大學(xué)學(xué)報(bào), 2014, 40(4):1-6.

    LIU Chunlei, YIN Yong, SUN Xiaofeng,etal. A modified method to ship floating calculation [J]. Journal of Dalian Maritime University, 2014, 40(4):1-6. (in Chinese)

    [3] CHEN Jing, LIN Yan, HUO Junzhou,etal. Optimization of ship′s subdivision arrangement for offshore sequential ballast water exchange using a non-dominated sorting genetic algorithm [J]. Ocean Engineering, 2010, 37(11/12):978-988.

    [4] CHEN Jing, LIN Yan, HUO Junzhou,etal. Optimal ballast water exchange sequence design using symmetrical multitank strategy [J]. Journal of Marine Science and Technology, 2010, 15(3):280-293.

    [5] 趙曉非,林 焰. 關(guān)于解船舶浮態(tài)問(wèn)題的矩陣方法[J]. 中國(guó)造船, 1985(3):55-64.

    ZHAO Xiaofei, LIN Yan. Matrix methods for solving ship floating state problem [J]. Shipbuilding of China, 1985(3):55-64. (in Chinese)

    [6] 林 焰,李鐵驪,紀(jì)卓尚. 破損船舶自由浮態(tài)計(jì)算[J]. 大連理工大學(xué)學(xué)報(bào), 2001, 41(1):85-89.

    LIN Yan, LI Tieli, JI Zhuoshang. Ship damage floating calculation [J]. Journal of Dalian University of Technology, 2001, 41(1):85-89. (in Chinese)

    [7] 馬 坤,張明霞,紀(jì)卓尚. 基于非線性規(guī)劃法的船舶浮態(tài)計(jì)算[J]. 大連理工大學(xué)學(xué)報(bào), 2003, 43(3):329-331.

    MA Kun, ZHANG Mingxia, JI Zhuoshang. Ship floating calculation based on nonlinear programming [J]. Journal of Dalian University of Technology, 2003, 43(3):329-331. (in Chinese)

    [8] 孫承猛,劉寅東. 船舶浮態(tài)計(jì)算的一種優(yōu)化方法[J]. 大連海事大學(xué)學(xué)報(bào), 2006, 32(2):39-41.

    SUN Chengmeng, LIU Yindong. Optimization method for ship′s floating condition calculation [J]. Journal of Dalian Maritime University, 2006, 32(2):39-41. (in Chinese)

    [9] 陸叢紅,林 焰,紀(jì)卓尚. 遺傳算法在船舶自由浮態(tài)計(jì)算中的應(yīng)用[J]. 上海交通大學(xué)學(xué)報(bào), 2005, 39(5):701-705,710.

    LU Conghong, LIN Yan, JI Zhuoshang. Application of genetic algorithm in ship free floatation calculation [J]. Journal of Shanghai Jiaotong University, 2005, 39(5):701-705, 710. (in Chinese)

    [10] LU Conghong, LIN Yan, JI Zhuoshang. Free trim calculation using genetic algorithm based on NURBS shipform [J]. International Shipbuilding Progress, 2007, 54(1):45-62.

    [11] LU Conghong, LIN Yan, JI Zhuoshang. An algorithm for plane-surface intersection and its application to shipbuilding [J]. Ship Technology Research, 2005, 52(3):103-106.

    [12] CALABRESE F, MANCARELLA L, ZIZZARI A A,etal. A multidisciplinary method for evaluating ship stability [J]. Journal of Shipping and Ocean Engineering, 2012, 2(6):321-326.

    [13] CALABRESE F, CATALDO M, DE PASCALIS A,etal. A method for evaluating the ship equilibrium point using 3D meshes developed for embedded device [C] // MTS/IEEE OCEANS 2015-Genova: Discovering Sustainable Ocean Energy for a New World. Piscataway:IEEE, 2015:7271763.

    [14] AZCUETA R. Computation of turbulent free-surface flows around ships and floating bodies [J]. Ship Technology Research, 2002, 49(2):46-69.

    Simplified method of bulk carrier′s free floatation calculation

    LIU Chunlei, YIN Yong*, SUN Xiaofeng, ZHANG Xiufeng, SHEN Helong

    ( Key Laboratory of Marine Simulation and Control Ministry of Communication, Dalian Maritime University, Dalian 116026, China )

    According to the characteristics of bulk carrier, a simplified method of ship free floatation calculation for engineering application is studied. Based on matrix method, 7 waterline plane elements are reduced to 4, and the solution of Jacobi matrix is simplified. Complete Gaussian pivoting elimination method is adopted to solve the linear equations. In order to ensure the accuracy of the calculation, the floating state balance equations are used as the iteration stop conditions. Bulk carrier “DOLCE VITA” and “RUI AN CHENG” are taken as examples, 200 000 random loading experiments are carried out on 20 ballast tanks to verify the convergence and robustness of the algorithm in heeling condition. 200 000 random loading experiments are carried out on 5 cargo holds to verify the convergence and robustness of the algorithm in trim condition. Typical load conditions are calculated, and compared with the given value in the loading manual the trim error is less than 0.01 m. The calculation results show that in this algorithm the calculating amount is decreased because only displacement volume under any tilt waterline plane and center of buoyancy need to be computed, the calculation amount is small, the program is very simple, real-time, robust and has a certain engineering application value.

    ship engineering; free floatation; bulk carrier; loading computer; trim

    1000-8608(2017)01-0055-06

    2016-05-21;

    2016-11-24.

    “八六三”國(guó)家高技術(shù)研究發(fā)展計(jì)劃資助項(xiàng)目(2015AA016404);海洋公益性行業(yè)科研專項(xiàng)(201505017-4);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(3132016310).

    劉春雷(1987-),男,博士生,E-mail:18940930042@163.com;尹 勇*(1969-),男,教授,E-mail:bushyin@163.com.

    U661.21

    A

    10.7511/dllgxb201701008

    猜你喜歡
    水線散貨船傾角
    地球軸傾角的改斜歸正
    激光傾角儀在CT引導(dǎo)下經(jīng)皮肺穿刺活檢中的應(yīng)用
    車輪外傾角和前束角匹配研究
    北京汽車(2021年3期)2021-07-17 08:25:42
    20艘!今年好望角型散貨船拆解量遠(yuǎn)超去年同期
    基于VECM的國(guó)際干散貨船投資
    水線光纜
    水線光纜
    規(guī)?;怆u養(yǎng)殖場(chǎng)水線管理
    79800 DWT雙燃料散貨船LNG儲(chǔ)氣罐位置確定
    船海工程(2015年4期)2016-01-05 15:53:50
    5 000 t散貨船清潔推進(jìn)系統(tǒng)配置論證與控制設(shè)計(jì)
    人人妻,人人澡人人爽秒播| 三级毛片av免费| 成人av一区二区三区在线看| 亚洲图色成人| 国产高清有码在线观看视频| 色综合亚洲欧美另类图片| 日产精品乱码卡一卡2卡三| 寂寞人妻少妇视频99o| 亚洲色图av天堂| 联通29元200g的流量卡| 中出人妻视频一区二区| 国产亚洲av嫩草精品影院| 成人漫画全彩无遮挡| eeuss影院久久| 成年女人看的毛片在线观看| 国产伦一二天堂av在线观看| 直男gayav资源| 在现免费观看毛片| 我要看日韩黄色一级片| www.色视频.com| 日韩av不卡免费在线播放| 免费人成在线观看视频色| 春色校园在线视频观看| 我要看日韩黄色一级片| 一a级毛片在线观看| av中文乱码字幕在线| 午夜福利高清视频| 精品一区二区三区视频在线观看免费| 国产精品一区二区三区四区免费观看 | 高清日韩中文字幕在线| 男人舔奶头视频| 久久久午夜欧美精品| 国产免费一级a男人的天堂| 国产精品不卡视频一区二区| 国产精品野战在线观看| 麻豆精品久久久久久蜜桃| 一级毛片我不卡| 超碰av人人做人人爽久久| 国产亚洲91精品色在线| 国产男靠女视频免费网站| 一进一出抽搐gif免费好疼| 国内久久婷婷六月综合欲色啪| 国产成人a∨麻豆精品| 婷婷六月久久综合丁香| 国产精品久久久久久久电影| 中文在线观看免费www的网站| 国模一区二区三区四区视频| 国产私拍福利视频在线观看| av在线蜜桃| 日韩欧美在线乱码| 日本三级黄在线观看| 久久午夜亚洲精品久久| 综合色av麻豆| 欧美日韩在线观看h| 看十八女毛片水多多多| 国产色爽女视频免费观看| 亚洲欧美日韩高清专用| 亚洲自偷自拍三级| 精品人妻偷拍中文字幕| 亚洲无线观看免费| 在线播放国产精品三级| 久久久国产成人免费| 一夜夜www| av免费在线看不卡| 国产午夜福利久久久久久| 免费av观看视频| 51国产日韩欧美| or卡值多少钱| 在线观看一区二区三区| 国产精品亚洲一级av第二区| 亚洲av成人精品一区久久| 国产日本99.免费观看| 午夜福利在线观看吧| 精品久久国产蜜桃| 神马国产精品三级电影在线观看| 亚洲国产日韩欧美精品在线观看| 啦啦啦啦在线视频资源| 男插女下体视频免费在线播放| 久久久精品94久久精品| 国内揄拍国产精品人妻在线| 91麻豆精品激情在线观看国产| 国产精品一区二区三区四区久久| 中国国产av一级| а√天堂www在线а√下载| 亚洲av免费高清在线观看| 国产片特级美女逼逼视频| 在线免费观看不下载黄p国产| 精品久久久噜噜| 国产单亲对白刺激| 午夜激情福利司机影院| 九色成人免费人妻av| 国产男人的电影天堂91| 精品午夜福利在线看| 久久婷婷人人爽人人干人人爱| 国产高清激情床上av| 久久这里只有精品中国| 搡女人真爽免费视频火全软件 | 最新中文字幕久久久久| 97人妻精品一区二区三区麻豆| 亚洲精品一区av在线观看| 久久久午夜欧美精品| 日本五十路高清| av中文乱码字幕在线| 国产三级在线视频| 日本熟妇午夜| 免费看光身美女| 精品久久久久久久久久免费视频| 中文在线观看免费www的网站| 精品日产1卡2卡| 中文字幕免费在线视频6| 九九爱精品视频在线观看| 午夜福利视频1000在线观看| 午夜免费激情av| 日韩亚洲欧美综合| 精品欧美国产一区二区三| 97超级碰碰碰精品色视频在线观看| 国产乱人偷精品视频| 日韩欧美在线乱码| 久久精品91蜜桃| 老熟妇仑乱视频hdxx| www.色视频.com| 亚洲自拍偷在线| 网址你懂的国产日韩在线| 午夜日韩欧美国产| 国内精品久久久久精免费| 99热全是精品| 国产亚洲精品久久久com| 给我免费播放毛片高清在线观看| 国产精品久久久久久av不卡| 偷拍熟女少妇极品色| 听说在线观看完整版免费高清| 亚洲最大成人av| 亚洲熟妇熟女久久| 国产一区亚洲一区在线观看| 久久国产乱子免费精品| 欧美+亚洲+日韩+国产| 色哟哟哟哟哟哟| 精品一区二区三区人妻视频| 国内少妇人妻偷人精品xxx网站| 午夜福利在线观看免费完整高清在 | 女人十人毛片免费观看3o分钟| 亚洲七黄色美女视频| 免费搜索国产男女视频| 亚洲精品影视一区二区三区av| 国产精品久久久久久av不卡| 国产成人a∨麻豆精品| 男插女下体视频免费在线播放| 日本 av在线| av在线观看视频网站免费| 男女啪啪激烈高潮av片| 免费看a级黄色片| 国产精品福利在线免费观看| 2021天堂中文幕一二区在线观| 2021天堂中文幕一二区在线观| 国产 一区精品| 国产爱豆传媒在线观看| 国产亚洲精品久久久久久毛片| 欧美一级a爱片免费观看看| 又爽又黄a免费视频| 亚洲av不卡在线观看| 国产精品一二三区在线看| 久久精品夜色国产| 亚洲成人久久爱视频| 麻豆成人午夜福利视频| 99久久精品国产国产毛片| 观看美女的网站| 日本免费a在线| 一级毛片久久久久久久久女| 欧美一区二区亚洲| 欧美性猛交╳xxx乱大交人| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品久久国产高清桃花| 熟女人妻精品中文字幕| 成人二区视频| 国产在视频线在精品| 美女被艹到高潮喷水动态| 不卡视频在线观看欧美| 国产熟女欧美一区二区| 在线观看av片永久免费下载| 国产 一区精品| 丰满人妻一区二区三区视频av| 亚洲18禁久久av| 国产真实伦视频高清在线观看| 极品教师在线视频| 99久国产av精品| 亚洲av成人精品一区久久| 99精品在免费线老司机午夜| 国产探花极品一区二区| 夜夜看夜夜爽夜夜摸| 久久欧美精品欧美久久欧美| 热99在线观看视频| 亚洲av免费在线观看| 我要看日韩黄色一级片| 欧美zozozo另类| 1024手机看黄色片| 亚洲精品456在线播放app| 国产aⅴ精品一区二区三区波| 国产精品电影一区二区三区| 日本熟妇午夜| 免费观看在线日韩| 久久人人精品亚洲av| 亚洲精品粉嫩美女一区| 美女内射精品一级片tv| 精品99又大又爽又粗少妇毛片| 欧美人与善性xxx| 人人妻人人澡欧美一区二区| 成年女人永久免费观看视频| 国产亚洲av嫩草精品影院| 夜夜爽天天搞| 男人的好看免费观看在线视频| 国产美女午夜福利| 人妻久久中文字幕网| 中文字幕人妻熟人妻熟丝袜美| 天堂影院成人在线观看| 看黄色毛片网站| 久久鲁丝午夜福利片| 久久精品国产亚洲av涩爱 | 中出人妻视频一区二区| 真实男女啪啪啪动态图| 午夜免费男女啪啪视频观看 | 男女边吃奶边做爰视频| 日本爱情动作片www.在线观看 | 国产精品国产高清国产av| 国内精品美女久久久久久| 亚洲欧美日韩无卡精品| 国产精品一及| 狂野欧美激情性xxxx在线观看| 嫩草影院入口| 一进一出抽搐动态| 老司机午夜福利在线观看视频| 真实男女啪啪啪动态图| 久久久久免费精品人妻一区二区| 99视频精品全部免费 在线| av天堂在线播放| 九九久久精品国产亚洲av麻豆| 哪里可以看免费的av片| 国产一区二区三区av在线 | 精华霜和精华液先用哪个| 亚洲精品色激情综合| 搡老妇女老女人老熟妇| 久久国内精品自在自线图片| 久久久国产成人精品二区| 十八禁国产超污无遮挡网站| 国产精品99久久久久久久久| 午夜福利在线观看免费完整高清在 | 一本精品99久久精品77| 久久亚洲国产成人精品v| 免费观看在线日韩| 精华霜和精华液先用哪个| 波多野结衣高清作品| 欧美在线一区亚洲| 国产精品久久视频播放| 99热全是精品| 国产又黄又爽又无遮挡在线| 久久久久免费精品人妻一区二区| 成人特级黄色片久久久久久久| 亚洲欧美日韩无卡精品| АⅤ资源中文在线天堂| 毛片一级片免费看久久久久| 日韩 亚洲 欧美在线| 亚洲av不卡在线观看| 国产探花在线观看一区二区| 少妇高潮的动态图| 日韩欧美 国产精品| 成年女人看的毛片在线观看| 老熟妇仑乱视频hdxx| 精品人妻视频免费看| 丰满人妻一区二区三区视频av| 亚洲熟妇中文字幕五十中出| 久久久精品欧美日韩精品| av在线播放精品| 日本欧美国产在线视频| 亚洲成人久久爱视频| 一卡2卡三卡四卡精品乱码亚洲| 国产在线精品亚洲第一网站| 亚洲欧美成人精品一区二区| 黄色配什么色好看| av在线蜜桃| 国产av在哪里看| 美女大奶头视频| 国产精品久久久久久精品电影| 麻豆国产97在线/欧美| 亚洲av熟女| 色视频www国产| 成人精品一区二区免费| 嫩草影视91久久| 99热网站在线观看| 亚洲aⅴ乱码一区二区在线播放| 麻豆乱淫一区二区| 久久久a久久爽久久v久久| 欧美激情久久久久久爽电影| 日本 av在线| 噜噜噜噜噜久久久久久91| 一本一本综合久久| 欧美国产日韩亚洲一区| 国产精品亚洲美女久久久| 男人舔奶头视频| 精品久久国产蜜桃| 亚洲欧美日韩高清专用| 国产在视频线在精品| 看十八女毛片水多多多| 黄片wwwwww| 在线观看美女被高潮喷水网站| 在线免费观看不下载黄p国产| 啦啦啦韩国在线观看视频| 高清午夜精品一区二区三区 | 亚洲成人久久性| 色吧在线观看| 在线观看午夜福利视频| 色尼玛亚洲综合影院| 国内精品一区二区在线观看| 亚洲自拍偷在线| 欧美性感艳星| 99热这里只有是精品在线观看| 国产激情偷乱视频一区二区| 色哟哟哟哟哟哟| 天堂动漫精品| 毛片女人毛片| 乱码一卡2卡4卡精品| 99视频精品全部免费 在线| 黄片wwwwww| 日日摸夜夜添夜夜爱| 亚洲人成网站在线播| 亚洲,欧美,日韩| 18禁黄网站禁片免费观看直播| av在线天堂中文字幕| 一级黄片播放器| 国产亚洲精品综合一区在线观看| av.在线天堂| 99久久中文字幕三级久久日本| 免费搜索国产男女视频| 老司机福利观看| 天天一区二区日本电影三级| 国产一区二区三区av在线 | 精品午夜福利在线看| 最后的刺客免费高清国语| 一卡2卡三卡四卡精品乱码亚洲| 不卡视频在线观看欧美| 成人三级黄色视频| 欧美日韩乱码在线| 欧美3d第一页| 一个人看视频在线观看www免费| 久久久午夜欧美精品| 97在线视频观看| 久99久视频精品免费| 啦啦啦韩国在线观看视频| 国产成人精品久久久久久| 国产成年人精品一区二区| 俺也久久电影网| 天堂影院成人在线观看| 国产精品亚洲一级av第二区| 亚洲av.av天堂| 最近手机中文字幕大全| 国产综合懂色| 免费观看在线日韩| 国产高清视频在线观看网站| 日韩亚洲欧美综合| 免费人成在线观看视频色| 91久久精品电影网| 国产成人福利小说| av女优亚洲男人天堂| 可以在线观看毛片的网站| 22中文网久久字幕| www.色视频.com| 18禁在线无遮挡免费观看视频 | 又粗又爽又猛毛片免费看| 一级毛片电影观看 | 中文字幕久久专区| 国产高潮美女av| 国产精品日韩av在线免费观看| 国产精品美女特级片免费视频播放器| 可以在线观看毛片的网站| 91午夜精品亚洲一区二区三区| 激情 狠狠 欧美| 女人被狂操c到高潮| 特级一级黄色大片| 久久精品国产自在天天线| 国产精品久久久久久av不卡| 天堂影院成人在线观看| 三级国产精品欧美在线观看| 国产一区二区激情短视频| 中文字幕精品亚洲无线码一区| 欧美日本视频| 精品人妻视频免费看| 丝袜美腿在线中文| 尤物成人国产欧美一区二区三区| 夜夜看夜夜爽夜夜摸| 哪里可以看免费的av片| 最近的中文字幕免费完整| 少妇的逼水好多| 国产精品福利在线免费观看| 嫩草影院新地址| 校园春色视频在线观看| 婷婷精品国产亚洲av| 亚洲综合色惰| 亚洲四区av| 国产色爽女视频免费观看| 在线a可以看的网站| 久久婷婷人人爽人人干人人爱| 天天躁日日操中文字幕| 99热这里只有精品一区| 久久久久国产网址| 色噜噜av男人的天堂激情| 日日摸夜夜添夜夜爱| 精品99又大又爽又粗少妇毛片| 亚洲国产精品合色在线| 精华霜和精华液先用哪个| 国产高清有码在线观看视频| 晚上一个人看的免费电影| 大香蕉久久网| 草草在线视频免费看| 国模一区二区三区四区视频| 午夜久久久久精精品| 亚洲中文字幕日韩| 亚洲自偷自拍三级| 一本久久中文字幕| 亚洲,欧美,日韩| 人人妻人人澡欧美一区二区| 禁无遮挡网站| 午夜精品在线福利| 精品一区二区三区av网在线观看| 不卡一级毛片| 国产一区二区激情短视频| 色av中文字幕| 22中文网久久字幕| 日本免费a在线| 成人美女网站在线观看视频| 久久久午夜欧美精品| 日本爱情动作片www.在线观看 | 午夜福利在线在线| 欧美在线一区亚洲| 国产在线精品亚洲第一网站| 午夜福利在线在线| 卡戴珊不雅视频在线播放| 黄色欧美视频在线观看| av在线观看视频网站免费| 女人十人毛片免费观看3o分钟| 一级av片app| 日韩成人伦理影院| 亚洲欧美成人综合另类久久久 | 直男gayav资源| 国产高清视频在线播放一区| 综合色丁香网| 99热这里只有是精品50| 最好的美女福利视频网| 午夜福利在线观看吧| 麻豆一二三区av精品| 久久精品夜夜夜夜夜久久蜜豆| 亚洲国产精品成人久久小说 | 国产激情偷乱视频一区二区| 欧美日韩精品成人综合77777| 毛片女人毛片| 国产探花极品一区二区| 亚洲av美国av| 男人狂女人下面高潮的视频| 全区人妻精品视频| 亚洲最大成人av| 国产在视频线在精品| 国产av在哪里看| 俺也久久电影网| 色5月婷婷丁香| 亚洲自偷自拍三级| 夜夜夜夜夜久久久久| 黄色一级大片看看| 国产三级中文精品| 日日摸夜夜添夜夜添av毛片| 国产精品无大码| 日本欧美国产在线视频| 日韩国内少妇激情av| 亚洲自偷自拍三级| 国产精品久久久久久久久免| 精品一区二区三区视频在线观看免费| 我要搜黄色片| 天堂网av新在线| www日本黄色视频网| 简卡轻食公司| h日本视频在线播放| 久久久精品94久久精品| 99久久成人亚洲精品观看| 日本五十路高清| 给我免费播放毛片高清在线观看| 国产精品99久久久久久久久| 国产高清视频在线观看网站| 国产91av在线免费观看| 91午夜精品亚洲一区二区三区| 18禁在线无遮挡免费观看视频 | 亚洲一级一片aⅴ在线观看| 欧美色欧美亚洲另类二区| 久久精品91蜜桃| 美女内射精品一级片tv| 欧美日韩国产亚洲二区| 日日撸夜夜添| 晚上一个人看的免费电影| 在线观看av片永久免费下载| 午夜福利在线观看吧| 国产成人一区二区在线| 成人特级av手机在线观看| 一本精品99久久精品77| 一级a爱片免费观看的视频| 国产精品一区www在线观看| 精品不卡国产一区二区三区| 观看美女的网站| 国产精品野战在线观看| 我的女老师完整版在线观看| 性欧美人与动物交配| 一进一出抽搐gif免费好疼| 成人无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看| 97在线视频观看| 日本一本二区三区精品| 男女下面进入的视频免费午夜| 99热全是精品| 1000部很黄的大片| 插阴视频在线观看视频| 精品福利观看| av在线播放精品| 国产成人一区二区在线| 久久精品综合一区二区三区| 国产高清三级在线| 国产高清激情床上av| 成人二区视频| 波多野结衣高清无吗| 亚洲欧美日韩高清专用| 成人亚洲欧美一区二区av| 最近中文字幕高清免费大全6| 在线国产一区二区在线| 亚洲av成人av| 乱码一卡2卡4卡精品| 久久这里只有精品中国| 国产精品久久久久久久久免| 人人妻,人人澡人人爽秒播| 麻豆成人午夜福利视频| 在线天堂最新版资源| 久久午夜福利片| 日本熟妇午夜| 偷拍熟女少妇极品色| 免费观看精品视频网站| 99热网站在线观看| 老司机福利观看| av国产免费在线观看| 日本黄色片子视频| 亚洲真实伦在线观看| 最近2019中文字幕mv第一页| 麻豆久久精品国产亚洲av| 99riav亚洲国产免费| 伦精品一区二区三区| 最近视频中文字幕2019在线8| 免费一级毛片在线播放高清视频| 天堂动漫精品| 精品一区二区三区av网在线观看| 国产成人91sexporn| 高清毛片免费看| 最新中文字幕久久久久| 中文字幕免费在线视频6| av天堂在线播放| 99热全是精品| av在线观看视频网站免费| 久久九九热精品免费| 亚洲中文字幕日韩| 一级毛片我不卡| 三级男女做爰猛烈吃奶摸视频| 免费人成视频x8x8入口观看| 欧美激情久久久久久爽电影| 日本a在线网址| 久久鲁丝午夜福利片| 亚洲av熟女| 黄色一级大片看看| 欧美激情在线99| 99久久九九国产精品国产免费| 日本爱情动作片www.在线观看 | 国产成人91sexporn| 国产av一区在线观看免费| 久久99热这里只有精品18| 内地一区二区视频在线| 精品国产三级普通话版| 成年av动漫网址| 国内揄拍国产精品人妻在线| 欧美另类亚洲清纯唯美| 少妇被粗大猛烈的视频| 一卡2卡三卡四卡精品乱码亚洲| eeuss影院久久| 精品99又大又爽又粗少妇毛片| 男人舔女人下体高潮全视频| 国产v大片淫在线免费观看| 搡老岳熟女国产| 晚上一个人看的免费电影| 3wmmmm亚洲av在线观看| 精品午夜福利视频在线观看一区| 精品久久久久久久久亚洲| 久久草成人影院| eeuss影院久久| 久久精品国产亚洲av天美| 菩萨蛮人人尽说江南好唐韦庄 | 欧美日韩一区二区视频在线观看视频在线 | 可以在线观看的亚洲视频| 精品人妻偷拍中文字幕| 最近的中文字幕免费完整| 日产精品乱码卡一卡2卡三| 色视频www国产| 美女免费视频网站| 国产精品久久电影中文字幕| 日本 av在线| 中文在线观看免费www的网站| 久久久久久久久久黄片| 久久久久久久亚洲中文字幕| 国产国拍精品亚洲av在线观看| 国产在线男女| 欧美性感艳星| 伊人久久精品亚洲午夜| 色播亚洲综合网| 一级黄色大片毛片| 日韩亚洲欧美综合| 亚洲精品亚洲一区二区| 免费在线观看成人毛片| 亚洲精品日韩av片在线观看|