• <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又大| 纵有疾风起免费观看全集完整版| 两个人免费观看高清视频| 亚洲精品中文字幕在线视频| 一区二区三区乱码不卡18| 久久国产精品人妻蜜桃| 成人国产av品久久久| 精品人妻1区二区| a级毛片黄视频| 成人影院久久| av福利片在线| 欧美精品一区二区免费开放| 欧美在线一区亚洲| 久久国产精品大桥未久av| 欧美黑人欧美精品刺激| 欧美av亚洲av综合av国产av| 婷婷成人精品国产| 亚洲精品第二区| 色婷婷久久久亚洲欧美| 黄色毛片三级朝国网站| 我的亚洲天堂| 国产片内射在线| 日本猛色少妇xxxxx猛交久久| 成人手机av| 免费人妻精品一区二区三区视频| 久久免费观看电影| 国产福利在线免费观看视频| 香蕉丝袜av| 午夜福利在线免费观看网站| 色94色欧美一区二区| 国产精品国产av在线观看| 亚洲成人免费av在线播放| 男人舔女人的私密视频| 久久久久久久大尺度免费视频| 九色亚洲精品在线播放| 欧美人与性动交α欧美软件| 国产av国产精品国产| 日本av手机在线免费观看| 制服人妻中文乱码| av网站在线播放免费| 欧美日韩亚洲国产一区二区在线观看 | 日韩 亚洲 欧美在线| 王馨瑶露胸无遮挡在线观看| 91精品伊人久久大香线蕉| 成人18禁高潮啪啪吃奶动态图| 久久精品久久久久久噜噜老黄| 一区在线观看完整版| 熟女少妇亚洲综合色aaa.| 伊人亚洲综合成人网| 99国产精品一区二区蜜桃av | 亚洲欧洲精品一区二区精品久久久| 日韩大片免费观看网站| av片东京热男人的天堂| 一本—道久久a久久精品蜜桃钙片| 久9热在线精品视频| 波多野结衣一区麻豆| 69精品国产乱码久久久| 午夜福利在线观看吧| 老鸭窝网址在线观看| 欧美在线一区亚洲| av网站免费在线观看视频| 高清黄色对白视频在线免费看| 男女下面插进去视频免费观看| 久久中文字幕一级| 女人久久www免费人成看片| 黄色视频不卡| 久久国产精品男人的天堂亚洲| 性色av一级| 丝袜美足系列| 精品久久久久久电影网| 国产极品粉嫩免费观看在线| 日韩有码中文字幕| 国产精品秋霞免费鲁丝片| 性少妇av在线| av片东京热男人的天堂| 自拍欧美九色日韩亚洲蝌蚪91| 丰满人妻熟妇乱又伦精品不卡| www.精华液| 久久久久精品国产欧美久久久 | 精品久久久精品久久久| av不卡在线播放| 亚洲va日本ⅴa欧美va伊人久久 | 日本91视频免费播放| 亚洲第一av免费看| 极品少妇高潮喷水抽搐| 日韩制服丝袜自拍偷拍| 亚洲 国产 在线| 国产一级毛片在线| 嫩草影视91久久| 国产精品影院久久| tube8黄色片| 男女无遮挡免费网站观看| 一区二区日韩欧美中文字幕| av视频免费观看在线观看| 宅男免费午夜| 欧美黑人欧美精品刺激| 女性生殖器流出的白浆| 美女脱内裤让男人舔精品视频| 成在线人永久免费视频| 亚洲欧美精品自产自拍| 欧美xxⅹ黑人| 天天操日日干夜夜撸| 91国产中文字幕| 日韩制服骚丝袜av| 黄片大片在线免费观看| 久9热在线精品视频| 中文欧美无线码| 色婷婷av一区二区三区视频| 首页视频小说图片口味搜索| 久久久久视频综合| 老司机午夜福利在线观看视频 | 国产xxxxx性猛交| av天堂在线播放| 国产一区有黄有色的免费视频| 无遮挡黄片免费观看| 新久久久久国产一级毛片| av天堂在线播放| 91字幕亚洲| 国产一级毛片在线| 日韩一区二区三区影片| 欧美+亚洲+日韩+国产| 国产高清视频在线播放一区 | 亚洲av国产av综合av卡| 免费人妻精品一区二区三区视频| 一级,二级,三级黄色视频| a级毛片在线看网站| 亚洲精品日韩在线中文字幕| 亚洲熟女精品中文字幕| 国产成人av教育| 久久久久久亚洲精品国产蜜桃av| 91麻豆av在线| 久久久久精品国产欧美久久久 | 免费日韩欧美在线观看| 国产成+人综合+亚洲专区| 久久久久网色| 男人爽女人下面视频在线观看| 亚洲午夜精品一区,二区,三区| 9热在线视频观看99| 色老头精品视频在线观看| 亚洲va日本ⅴa欧美va伊人久久 | 精品国产国语对白av| 日日夜夜操网爽| 日日夜夜操网爽| 亚洲国产精品成人久久小说| 国产视频一区二区在线看| 另类精品久久| 黄色 视频免费看| 无遮挡黄片免费观看| 欧美黄色片欧美黄色片| 亚洲国产日韩一区二区| 国产三级黄色录像| 中文字幕人妻丝袜制服| 国产成人影院久久av| 三级毛片av免费| 免费久久久久久久精品成人欧美视频| 久久久精品94久久精品| 久久青草综合色| 岛国在线观看网站| 999久久久精品免费观看国产| 亚洲全国av大片| 国产精品麻豆人妻色哟哟久久| 午夜精品久久久久久毛片777| 人人妻人人澡人人看| 久久久久国产一级毛片高清牌| 亚洲欧洲日产国产| 女性被躁到高潮视频| 日韩制服丝袜自拍偷拍| 少妇被粗大的猛进出69影院| 19禁男女啪啪无遮挡网站| 免费高清在线观看日韩| 汤姆久久久久久久影院中文字幕| 亚洲国产av新网站| 大片电影免费在线观看免费| 欧美另类一区| 国产在视频线精品| 成年人午夜在线观看视频| 97精品久久久久久久久久精品| 老熟女久久久| 宅男免费午夜| 大型av网站在线播放| 美女主播在线视频| 亚洲国产中文字幕在线视频| 久久影院123| 欧美激情久久久久久爽电影 | 国产精品.久久久| 大香蕉久久网| tocl精华| 国产成人啪精品午夜网站| 国产免费一区二区三区四区乱码| 国产精品99久久99久久久不卡| 啦啦啦啦在线视频资源| 新久久久久国产一级毛片| 亚洲成人国产一区在线观看| 欧美日韩亚洲综合一区二区三区_| 国产黄频视频在线观看| 午夜激情久久久久久久| 天天添夜夜摸| 亚洲伊人色综图| 久久久欧美国产精品| 99国产综合亚洲精品| a在线观看视频网站| 老司机午夜十八禁免费视频| 久久ye,这里只有精品| 99九九在线精品视频| 亚洲精品国产区一区二| 亚洲欧美日韩另类电影网站| 亚洲熟女毛片儿| av免费在线观看网站| 午夜免费鲁丝| 一本一本久久a久久精品综合妖精| 精品亚洲成a人片在线观看| 国产xxxxx性猛交| 日韩熟女老妇一区二区性免费视频| 午夜福利在线观看吧| 最新在线观看一区二区三区| 久久99热这里只频精品6学生| 久久久精品区二区三区| 日韩电影二区| 99国产综合亚洲精品| 国产一区二区在线观看av| 久久久久久久久免费视频了| 欧美国产精品va在线观看不卡| 国产精品国产三级国产专区5o| tube8黄色片| 国产亚洲欧美在线一区二区| 免费高清在线观看视频在线观看| 人妻久久中文字幕网| 五月天丁香电影| 午夜福利视频在线观看免费| 人人妻人人澡人人爽人人夜夜| 麻豆国产av国片精品| 亚洲欧洲精品一区二区精品久久久| 99国产精品一区二区三区| 国产在线一区二区三区精| 欧美日本中文国产一区发布| 啪啪无遮挡十八禁网站| 久久久久国产一级毛片高清牌| 精品卡一卡二卡四卡免费| 三上悠亚av全集在线观看| 午夜福利,免费看| 黄片大片在线免费观看| 欧美日韩亚洲国产一区二区在线观看 | 人妻一区二区av| 男女边摸边吃奶| 国产精品免费大片| 亚洲欧美清纯卡通| av线在线观看网站| 女性生殖器流出的白浆| 免费人妻精品一区二区三区视频| 亚洲综合色网址| 免费日韩欧美在线观看| 97在线人人人人妻| 亚洲精品美女久久久久99蜜臀| av网站免费在线观看视频| 日韩一卡2卡3卡4卡2021年| 精品一区二区三区av网在线观看 | 欧美久久黑人一区二区| 久久久久国产精品人妻一区二区| 精品国产一区二区三区四区第35| 俄罗斯特黄特色一大片| 一进一出抽搐动态| 50天的宝宝边吃奶边哭怎么回事| av免费在线观看网站| 亚洲精品久久成人aⅴ小说| 性色av乱码一区二区三区2| 免费不卡黄色视频| 一级毛片女人18水好多| 黄色 视频免费看| 国产在线免费精品| 免费高清在线观看视频在线观看| 色综合欧美亚洲国产小说| 国产成+人综合+亚洲专区| 韩国精品一区二区三区| 性高湖久久久久久久久免费观看| 桃红色精品国产亚洲av| 法律面前人人平等表现在哪些方面 | 免费在线观看日本一区| 我的亚洲天堂| 99久久99久久久精品蜜桃| 人妻人人澡人人爽人人| www.自偷自拍.com| 91精品三级在线观看| 国产黄频视频在线观看| 多毛熟女@视频| 美女国产高潮福利片在线看| 青春草视频在线免费观看| 亚洲精品中文字幕一二三四区 | 成人黄色视频免费在线看| 国产精品 欧美亚洲| 亚洲专区国产一区二区| 国产免费av片在线观看野外av| 久久免费观看电影| 国产精品免费大片| 欧美人与性动交α欧美软件| 国产又爽黄色视频| 欧美日韩中文字幕国产精品一区二区三区 | 国产伦人伦偷精品视频| 日韩大片免费观看网站| 美女福利国产在线| 日韩欧美国产一区二区入口| videosex国产| bbb黄色大片| 亚洲av成人一区二区三| 亚洲伊人色综图| 性色av乱码一区二区三区2| 亚洲综合色网址| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩视频一区二区在线观看| 一区在线观看完整版| 人人妻人人添人人爽欧美一区卜| 18禁观看日本| 国产男女超爽视频在线观看| 极品人妻少妇av视频| 亚洲欧美清纯卡通| 国产又爽黄色视频| cao死你这个sao货| 美女扒开内裤让男人捅视频| 秋霞在线观看毛片| 啦啦啦免费观看视频1| 一级,二级,三级黄色视频| 正在播放国产对白刺激| 亚洲熟女毛片儿| 性高湖久久久久久久久免费观看| 国产熟女午夜一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产深夜福利视频在线观看| 黑丝袜美女国产一区| 国产成人欧美在线观看 | 欧美在线一区亚洲| 国产成人啪精品午夜网站| 亚洲午夜精品一区,二区,三区| 久久亚洲国产成人精品v| 亚洲色图 男人天堂 中文字幕| 看免费av毛片| 丝袜在线中文字幕| 日韩熟女老妇一区二区性免费视频| 最新在线观看一区二区三区| 黑人欧美特级aaaaaa片| 久久人人爽人人片av| 亚洲男人天堂网一区| 青草久久国产| 国产高清国产精品国产三级| 久久久精品区二区三区| videosex国产| 日韩大片免费观看网站| 欧美日韩av久久| 亚洲,欧美精品.| 女警被强在线播放| 人妻人人澡人人爽人人| 亚洲国产看品久久| 成人手机av| 免费一级毛片在线播放高清视频 | 在线永久观看黄色视频| 美女大奶头黄色视频| 丰满人妻熟妇乱又伦精品不卡| 国产成人欧美| 视频区欧美日本亚洲| 各种免费的搞黄视频| 一级黄色大片毛片| 五月开心婷婷网| 成人国产av品久久久| 男女国产视频网站| 一边摸一边做爽爽视频免费| 高清视频免费观看一区二区| 一区在线观看完整版| 一区二区三区乱码不卡18| 在线观看免费视频网站a站| 青草久久国产| 真人做人爱边吃奶动态| 国产真人三级小视频在线观看| 久久久久久久久久久久大奶| 成年女人毛片免费观看观看9 | 午夜影院在线不卡| 永久免费av网站大全| av天堂久久9| 国产精品久久久久久精品古装| 人人妻人人澡人人爽人人夜夜| 欧美一级毛片孕妇| 亚洲专区国产一区二区| 精品人妻熟女毛片av久久网站| 精品熟女少妇八av免费久了| 老熟妇乱子伦视频在线观看 | 欧美精品人与动牲交sv欧美| 搡老乐熟女国产| 999精品在线视频| 老司机影院毛片| 高清在线国产一区| 中文字幕制服av| 肉色欧美久久久久久久蜜桃| 久久综合国产亚洲精品| 亚洲av成人不卡在线观看播放网 | 男男h啪啪无遮挡| 999精品在线视频| 美女脱内裤让男人舔精品视频| 欧美少妇被猛烈插入视频| 飞空精品影院首页| 一个人免费在线观看的高清视频 | 国产免费一区二区三区四区乱码| 12—13女人毛片做爰片一| 99国产综合亚洲精品| 亚洲综合色网址| 老汉色∧v一级毛片| 中文字幕另类日韩欧美亚洲嫩草| 久久国产精品人妻蜜桃| 国产精品影院久久| 人成视频在线观看免费观看| 国产1区2区3区精品| 久久亚洲精品不卡| 精品国产乱码久久久久久小说| 欧美激情久久久久久爽电影 | 黄色毛片三级朝国网站| 免费少妇av软件| tocl精华| www日本在线高清视频| 亚洲国产精品999| 中文字幕色久视频| 少妇粗大呻吟视频| 亚洲国产精品一区二区三区在线| 久久综合国产亚洲精品| 丁香六月天网| 欧美日韩视频精品一区| 下体分泌物呈黄色| 久久久国产精品麻豆| 亚洲第一欧美日韩一区二区三区 | 丰满迷人的少妇在线观看| 涩涩av久久男人的天堂| 亚洲精华国产精华精| 中文字幕高清在线视频| 蜜桃在线观看..| 久久久久久亚洲精品国产蜜桃av| 9热在线视频观看99| 日韩三级视频一区二区三区| 欧美中文综合在线视频| 国产又爽黄色视频| 日韩大片免费观看网站| 亚洲专区中文字幕在线| 最新的欧美精品一区二区| 少妇 在线观看| 性色av乱码一区二区三区2| 丁香六月欧美| 考比视频在线观看| 少妇的丰满在线观看| 亚洲精品在线美女| 狠狠狠狠99中文字幕| av在线播放精品| 久久精品熟女亚洲av麻豆精品| 欧美成狂野欧美在线观看| 久久精品亚洲熟妇少妇任你| 高清视频免费观看一区二区| 亚洲成国产人片在线观看| av视频免费观看在线观看| 一二三四在线观看免费中文在| 色94色欧美一区二区| 久久国产精品人妻蜜桃| 欧美在线黄色| 国产极品粉嫩免费观看在线| 建设人人有责人人尽责人人享有的| 国产成人精品久久二区二区91| netflix在线观看网站| 国产在线免费精品| 亚洲情色 制服丝袜| 日韩欧美一区视频在线观看| 国产精品亚洲av一区麻豆| 妹子高潮喷水视频| 色精品久久人妻99蜜桃| 不卡av一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 九色亚洲精品在线播放| 欧美日韩亚洲综合一区二区三区_| 黄片播放在线免费| av视频免费观看在线观看| 免费久久久久久久精品成人欧美视频| 亚洲av成人不卡在线观看播放网 | 妹子高潮喷水视频| 亚洲激情五月婷婷啪啪| 亚洲一区中文字幕在线| 丝袜美腿诱惑在线| 久久久久久久精品精品| 成年动漫av网址| 大码成人一级视频| 欧美+亚洲+日韩+国产| 男人操女人黄网站| 亚洲五月色婷婷综合| 国产福利在线免费观看视频| 亚洲av成人一区二区三| 岛国毛片在线播放| 久久久久久久久免费视频了| 欧美变态另类bdsm刘玥| 亚洲国产欧美在线一区| 精品亚洲乱码少妇综合久久| 国产精品秋霞免费鲁丝片| 美女视频免费永久观看网站| 久久精品成人免费网站| 久久久久网色| 狠狠婷婷综合久久久久久88av| 黄色毛片三级朝国网站| 国产精品二区激情视频| 多毛熟女@视频| h视频一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 国产xxxxx性猛交| 99热网站在线观看| 欧美 日韩 精品 国产| 欧美精品一区二区免费开放| 日日爽夜夜爽网站| 午夜免费观看性视频| 亚洲伊人久久精品综合| 男人舔女人的私密视频| 午夜精品久久久久久毛片777| 亚洲欧美激情在线| 麻豆av在线久日| 免费日韩欧美在线观看| 黄色视频,在线免费观看| 国产1区2区3区精品| 一二三四社区在线视频社区8| 国产片内射在线| 另类亚洲欧美激情| 国产男人的电影天堂91| av在线老鸭窝| 久久这里只有精品19| 久久人人97超碰香蕉20202| 欧美日韩视频精品一区| 性色av一级| 大香蕉久久成人网| 国产亚洲av高清不卡| 久久免费观看电影| 黄色视频在线播放观看不卡| 国产成人精品久久二区二区91| 国产老妇伦熟女老妇高清| 国产精品av久久久久免费| 午夜福利视频精品| 日韩三级视频一区二区三区| 啦啦啦在线免费观看视频4| 日本黄色日本黄色录像| 色视频在线一区二区三区| 久久av网站| 高清在线国产一区| 亚洲视频免费观看视频| 搡老乐熟女国产| 午夜福利,免费看| 亚洲国产精品一区三区| 国产一卡二卡三卡精品| 亚洲国产精品一区三区| 国产欧美日韩精品亚洲av| 欧美大码av| 深夜精品福利| 久久久久久久精品精品| 久久久久精品人妻al黑| 亚洲色图综合在线观看| 久久久国产精品麻豆| 王馨瑶露胸无遮挡在线观看| 午夜福利免费观看在线| 黑人操中国人逼视频| 叶爱在线成人免费视频播放| 夜夜骑夜夜射夜夜干| 最近最新中文字幕大全免费视频| 老司机亚洲免费影院| 69av精品久久久久久 | 亚洲精品一卡2卡三卡4卡5卡 | 久久综合国产亚洲精品| 美女高潮喷水抽搐中文字幕| 国产精品一区二区免费欧美 | 久久ye,这里只有精品| 精品少妇黑人巨大在线播放| av在线播放精品| 久久久国产精品麻豆| 又黄又粗又硬又大视频| 精品国产乱码久久久久久男人| 久久精品国产综合久久久| 香蕉丝袜av| 欧美 亚洲 国产 日韩一| 欧美日韩亚洲综合一区二区三区_| 老司机深夜福利视频在线观看 | 黄色视频在线播放观看不卡| 老鸭窝网址在线观看| 啦啦啦啦在线视频资源| 69精品国产乱码久久久| 成人国产一区最新在线观看| 热99re8久久精品国产| av在线老鸭窝| 免费久久久久久久精品成人欧美视频| 国产免费av片在线观看野外av| 人人澡人人妻人| 女性生殖器流出的白浆| 99九九在线精品视频| 免费观看av网站的网址| 欧美日韩av久久| tocl精华| 人人妻,人人澡人人爽秒播| 中文精品一卡2卡3卡4更新| 午夜福利乱码中文字幕| 久久国产精品人妻蜜桃| 亚洲av日韩精品久久久久久密| 中文字幕精品免费在线观看视频| 亚洲av欧美aⅴ国产| 妹子高潮喷水视频| 人人妻人人澡人人看| 欧美av亚洲av综合av国产av| 高清av免费在线| 在线观看免费视频网站a站| 色综合欧美亚洲国产小说| a级片在线免费高清观看视频| 脱女人内裤的视频| 亚洲人成77777在线视频| 免费在线观看视频国产中文字幕亚洲 | av福利片在线| 人人妻,人人澡人人爽秒播| 欧美另类亚洲清纯唯美| 国产一卡二卡三卡精品| av国产精品久久久久影院| 国产精品九九99| 欧美乱码精品一区二区三区|