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

    三體船阻力回歸分析及預(yù)報(bào)方法研究

    2010-06-21 01:02:26盧曉平
    船舶力學(xué) 2010年4期
    關(guān)鍵詞:體船實(shí)船船模

    王 中,盧曉平,王 瑋

    (1海軍工程大學(xué)船舶與動(dòng)力學(xué)院,武漢 430033;2海軍裝備研究院艦船所,北京 100073)

    三體船阻力回歸分析及預(yù)報(bào)方法研究

    王 中1,盧曉平1,王 瑋2

    (1海軍工程大學(xué)船舶與動(dòng)力學(xué)院,武漢 430033;2海軍裝備研究院艦船所,北京 100073)

    基于三體船船模系列阻力試驗(yàn)結(jié)果,以非均勻B樣條擬合船模阻力試驗(yàn)數(shù)據(jù)隨航速變化關(guān)系曲線,采用模長(zhǎng)界限控制的Gram-Schmidt正交化方法建立給定排水量和航速三體船阻力與側(cè)體位置關(guān)系的回歸模型,該方法能夠快速而有效地獲取有效范圍內(nèi)任意給定航速下三體船阻力隨側(cè)體位置的變化關(guān)系,具有一定的理論意義和較高的工程實(shí)用價(jià)值。

    三體船;阻力;回歸分析;非均勻B樣條

    1 引 言

    三體船成為近年研究熱點(diǎn)[1-5],各種三體船阻力算法紛紛被推出,但理論方法還不足以準(zhǔn)確預(yù)報(bào)三體船阻力,系列模型試驗(yàn)仍然是必要的手段。采用模型試驗(yàn)的方法只能獲取離散的速度點(diǎn)及側(cè)體位置情況下的阻力數(shù)據(jù),無法連續(xù)地反映阻力隨速度和側(cè)體位置變化的情況。而三體船阻力變化曲線較常規(guī)單體船復(fù)雜得多,將離散的模型試驗(yàn)數(shù)據(jù)用于阻力預(yù)報(bào)、阻力特性分析以及三體船工程設(shè)計(jì)較不方便。采用回歸分析的方法可以將離散信息連續(xù)化,從而全面地獲取阻力值隨航速、側(cè)體位置變化的信息,求出給定航速下阻力性能最優(yōu)的側(cè)體位置。本文基于三體船船模系列阻力試驗(yàn)結(jié)果,以非均勻B樣條擬合船模阻力試驗(yàn)數(shù)據(jù)隨航速變化關(guān)系曲線,采用模長(zhǎng)界限控制的Gram-Schmidt正交化方法[6]建立給定排水量和航速下的三體船阻力與側(cè)體位置關(guān)系回歸方程,以便快速而有效地獲取有效范圍內(nèi)給定航速下三體船阻力隨側(cè)體位置的變化關(guān)系。

    2 模長(zhǎng)界限控制的Gram-Schmidt正交化方法

    特別適用于多元多項(xiàng)式回歸的模長(zhǎng)界限控制的Gram-Schmidt正交化方法在船舶科研設(shè)計(jì)領(lǐng)域得到了廣泛應(yīng)用,它的主要特點(diǎn)是在數(shù)學(xué)上熟知的Gram-Schmidt正交化過程中加上模長(zhǎng)界限控制技巧。因此,本文選擇該方法建立回歸模型。

    三元多項(xiàng)式回歸方程可轉(zhuǎn)化為線性回歸方程如下:式中n0=(n1+1)(n2+1)(n3+1)-1。采用模長(zhǎng)界限控制的Gram-Schmidt正交化方法可以將(1)式進(jìn)一步化作正交多項(xiàng)式為基函數(shù)的多元線性回歸方程。模長(zhǎng)界限控制的Gram-Schmidt正交化是一個(gè)遞推過程,寫成矩陣形式

    現(xiàn)在如果把pk(k= 1 ,…,n )作為自變量,欲從有限次試驗(yàn)觀察值 (ηr, xr,yr,zr)r=1,…,m 來求 pk所對(duì)應(yīng)的偏回歸系數(shù)bk,則確定bk估計(jì)量k的正規(guī)方程為

    由于pk(k= 1 ,…,n )是正交多項(xiàng)式,所以P′P=In(In為n階單位矩陣)。 故由(4)式可以立即求得對(duì)應(yīng)于正交多項(xiàng)式pk的偏回歸系數(shù)的估計(jì)量為:

    在計(jì)算過程中可以通過手工或是自動(dòng)方式剔除不顯著變量。

    3 三體船阻力回歸分析模型的建立

    三體船坐標(biāo)系如圖1所示,橫向偏距為p,縱向偏距為a。給出三體船阻力隨p、a變化的圖譜或回歸公式對(duì)三體船的工程設(shè)計(jì)具有十分重要的意義。以下即采用上節(jié)介紹的方法來建立這類回歸方程。在以往的研究中發(fā)現(xiàn)剩余阻力系數(shù)隨p、a變化的規(guī)律十分復(fù)雜,剩余阻力系數(shù)隨p、a變化的曲線拐點(diǎn)較多,不平坦,簡(jiǎn)單的回歸模型不一定能很好地表示其變化規(guī)律;而總阻力隨p、a變化的趨勢(shì)相對(duì)平緩些,簡(jiǎn)單的回歸公式就能夠比較好地表示總阻力隨p、a的變化。由于船模與實(shí)船粘性不相似,對(duì)船??傋枇⒒貧w公式是沒有意義的,為此作者提出針對(duì)實(shí)船(設(shè)計(jì)目標(biāo)船)建立總阻力(或有效功率)回歸公式。

    在由試驗(yàn)數(shù)據(jù)換算實(shí)船阻力過程中采用傅汝德方法,即總阻力為摩擦阻力與剩余阻力之和。記高速三體船摩擦阻力為Rf=0.5ρU2S·Cf,其中S=S1+2S2,S為三體船浸濕面積,S1、S2分別為中體和側(cè)體的浸濕面積[7-8]。Cf按下式計(jì)算:

    其中Cf1、Cf2分別為中體和側(cè)體的摩擦阻力系數(shù),摩擦系數(shù)的計(jì)算采用ITTC-57公式或者Prandtl-Schlichting公式。實(shí)船粗糙度補(bǔ)貼系數(shù)為ΔCf=0.4×10-3。

    本文三體船阻力回歸分析及阻力預(yù)報(bào)的主要步驟歸納如下:

    (1)建立模型試驗(yàn)圖譜;

    (2)根據(jù)給定實(shí)船排水量及航速,基于非均勻B樣條[9-10]插值求取相應(yīng)速度下各個(gè)側(cè)體位置的船模阻力值,將模型圖譜換算到待求實(shí)船;

    (3)對(duì)實(shí)船圖譜回歸分析;

    (4)根據(jù)回歸公式計(jì)算實(shí)船總阻力值,或分析三體船排水量、航速給定條件下實(shí)船阻力隨p、a變化規(guī)律,以確定最佳、較佳或可行的p、a設(shè)計(jì)值。

    4 實(shí)例分析

    作者基于Acess數(shù)據(jù)庫,采用VC#開發(fā)了基于系列船模阻力試驗(yàn)數(shù)據(jù)的三體船阻力回歸及預(yù)報(bào)系統(tǒng)。該系統(tǒng)具有良好的人機(jī)操作界面,只需用戶輸入給定的實(shí)船航速及排水量,系統(tǒng)便可計(jì)算出該航速及排水量下的阻力與側(cè)體位置關(guān)系的回歸方程,同時(shí)可以計(jì)算出有效范圍內(nèi)給定側(cè)體位置的三體船阻力值。

    本文基于中體和側(cè)體均為Wigley船型的高速三體船模在Fn=0.1~0.8時(shí)3個(gè)橫向偏距、5個(gè)縱向偏距共15個(gè)狀態(tài)的阻力試驗(yàn)數(shù)據(jù)進(jìn)行回歸分析,建立三體實(shí)船總阻力隨 p、a變化的回歸方程,船模阻力試驗(yàn)于2006年6月在中國特種飛行器研究所船池實(shí)施。Wigley數(shù)學(xué)三體船主體和側(cè)體船型表達(dá)式如(7)式所示。其中主體 L=5m,B=0.4m,T=0.178m,側(cè)體 L=1.842m,B=0.147 4m,T=0.066m,主體與單個(gè)小側(cè)體排水量之比為19.86。

    采用上述方法,我們給出4 000噸級(jí)三體船對(duì)應(yīng)航速32Kn,5 000噸級(jí)三體船對(duì)應(yīng)航速20Kn的阻力回歸結(jié)果如下:

    (8)式為三體船阻力與排水量無因次化比值隨側(cè)體位置的回歸表達(dá)式,相應(yīng)系數(shù)如表1所示?;貧w方程通常只能在側(cè)體位置參數(shù)變化范圍內(nèi)有效(即只適用于內(nèi)插情況),表2給出了相應(yīng)回歸模型表達(dá)式的有效范圍。某些情況下可能需要應(yīng)用回歸公式進(jìn)行外推,此時(shí)應(yīng)該慎重,外推可能產(chǎn)生較大誤差,甚至出現(xiàn)完全無效的結(jié)果。

    表1 多項(xiàng)式系數(shù)表Tab.1 Polynomial coefficient

    表2 回歸公式適用范圍(單位:m)Tab.2 Range of the regressive expressions(unit:m)

    圖2、3分別給出了兩種情況下實(shí)船阻力與側(cè)體位置關(guān)系的等值線云圖,每個(gè)圖中(a)圖為由試驗(yàn)數(shù)據(jù)直接繪制的等值線云圖,(b)圖為根據(jù)回歸模型繪制的等值線云圖。從圖2、3可以看出,阻力回歸模型能夠很逼真地反應(yīng)模型試驗(yàn)實(shí)際測(cè)量預(yù)報(bào)結(jié)果,同時(shí),更為清晰明確地反應(yīng)了三體船阻力與側(cè)體的位置關(guān)系,上述回歸圖譜可直接用于三體船工程設(shè)計(jì)。

    5 基于回歸模型分析側(cè)體位置對(duì)三體船阻力的影響

    根據(jù)本文方法,作者對(duì)5 000噸級(jí)三體船在18~40Kn航速范圍內(nèi)的總阻力與側(cè)體位置關(guān)系進(jìn)行回歸,系統(tǒng)研究了側(cè)體位置對(duì)三體船阻力的影響,如圖4所示。

    從圖4可以看出,對(duì)本船型從總阻力最小進(jìn)行側(cè)體布局優(yōu)化的角度看,20Kn中低航速時(shí)側(cè)體縱向位置應(yīng)布置在主體中部稍后,30Kn以上高航速時(shí)側(cè)體應(yīng)布置在主體后側(cè),隨著航速的增大橫向位置向主船體靠近,三體船阻力的這些特性與文獻(xiàn)[2]得到的結(jié)論一致;此外上述5 000噸級(jí)三體船阻力回歸圖譜系列中曲線趨勢(shì)協(xié)調(diào),其規(guī)律與模型試驗(yàn)原始數(shù)據(jù)預(yù)報(bào)結(jié)果吻合很好,該圖譜可直接應(yīng)用于工程,指導(dǎo)三體船型減阻優(yōu)化設(shè)計(jì)。

    6 結(jié) 語

    本文從三體船船模系列阻力試驗(yàn)結(jié)果出發(fā),以非均勻B樣條擬合船模阻力試驗(yàn)數(shù)據(jù)隨航速變化關(guān)系曲線,采用模長(zhǎng)界限控制的正交化方法建立給定航速下三體船阻力與側(cè)體位置關(guān)系的回歸方程,具有一定的理論意義和較高的工程實(shí)用價(jià)值。

    所采用的模長(zhǎng)界限控制的Gram-Schmidt正交化方法在船舶科研設(shè)計(jì)領(lǐng)域被廣泛應(yīng)用,具有較好的可靠性和穩(wěn)定性。

    對(duì)由船模阻力試驗(yàn)結(jié)果換算得到的給定排水量的實(shí)船總阻力進(jìn)行回歸分析,能更真實(shí)、清晰、明確反映實(shí)船阻力特性。該方法針對(duì)任意給定實(shí)船排水量以及航速進(jìn)行回歸分析,運(yùn)算量非常大,采用傳統(tǒng)的各個(gè)步驟分散的計(jì)算方法效率比較低,而作者將整套方法集成到一個(gè)界面友好的軟件系統(tǒng)中則可以大大簡(jiǎn)化操作,任意改變排水量和航速可以快速得到相應(yīng)的回歸方程。

    建立的整套回歸分析方法適用于任意三體船模型試驗(yàn)數(shù)據(jù)到實(shí)船的換算,所開發(fā)的基于系列船模阻力試驗(yàn)數(shù)據(jù)的三體船阻力回歸及預(yù)報(bào)系統(tǒng)可以直接用于各種三體船型,只需加入新的三體船船模阻力試驗(yàn)數(shù)據(jù)即可進(jìn)行實(shí)船阻力性能的回歸分析和預(yù)報(bào),具有較好的通用性。

    在試驗(yàn)?zāi)P妥銐蚨嗟臈l件下,采用本文的回歸分析方法和軟件系統(tǒng)還可望建立系統(tǒng)的三體船阻力優(yōu)化設(shè)計(jì)圖譜。

    [1]Robert P,Saunders Jr.An investigation of the resistance properties of a modern trimaran combatant ship base on Taylor Standard Series and Series64[D].Monterey,CA:Master’s Thesis,Naval Postgraduate School,1995.

    [2]Elcin Z.Wave making resistance characteristics of trimaran hulls[D].Monterey,CA:Master’s Thesis,Naval Postgraduate School,2003.

    [3]盧曉平,潘雨村.高速三體船興波阻力與片體布局優(yōu)化研究[J].水動(dòng)力學(xué)研究與進(jìn)展,A輯,2004,19(3):347-359.

    [4]李云波,陳 康,黃德波.三體船粘性阻力計(jì)算與計(jì)算方法比較[J].水動(dòng)力學(xué)研究與進(jìn)展,A輯,2005,20(4):452-457.

    [5]Wang Hu,Zou Zaojian.Numerical research on wave-making resistance of trimaran[J].J Shanghai Jiaotong Univ.(Sci.),2008,13(3):348-351.

    [6]周連第.多元回歸分析方法及在船舶科研設(shè)計(jì)中的應(yīng)用[M].北京:國防工業(yè)出版社,1979.

    [7]酈 云,盧曉平.高速三體船阻力性能研究[J].船舶力學(xué),2007,11(2):191-198.

    [8]李培勇,裘泳銘,顧敏童等.三體船阻力模型試驗(yàn)[J].中國造船,2002,43(4):6-12.

    [9]Piegl L.On NURBS:A survey[J].IEEE Computer Graphics and Applications,1991,11(1):55-71.

    [10]施法中.計(jì)算機(jī)輔助幾何設(shè)計(jì)與非均勻有理B樣條[M].北京:高等教育出版社,2001.

    Research on regression and prediction of the trimaran resistance

    WANG Zhong1,LU Xiao-ping1,WANG Wei2

    (1 Ship and Power Eng.College,Naval Univ.of Engineering,Wuhan 430033,China;
    2 Naval Academy of Armament,Beijing 100073,China)

    By using Non-Uniform-B-Spline to express the resistance curve and adopting the Gram-Schmidt method with finite condition,the regression of the relationship between the trimaran resistance and the side hull position was made,which is based on a series trimaran model experiment results.With the method put forward,the resistance of the trimaran with certain displacement and speed could be calculated quickly,and it is useful to solve such kind of problem of engineering design.

    trimaran;resistance;regression;Non-Uniform-B-Spline

    U661.31

    A

    1007-7294(2010)04-0355-07

    2009-09-10

    王 中(1981-),男,博士研究生。主要從事船舶與海洋工程流體力學(xué)研究。

    猜你喜歡
    體船實(shí)船船模
    船舶阻力試驗(yàn)實(shí)船換算的蒙特卡羅仿真
    基于模糊PID的船模航向控制研究
    某雙體客滾船軸系安裝工藝
    廣東造船(2021年2期)2021-05-23 12:16:39
    長(zhǎng)江1 000車滾裝船EEDI驗(yàn)證分析
    船海工程(2019年6期)2019-12-25 06:12:50
    5 668 TEU集裝箱船球艏改型與浮態(tài)節(jié)能的實(shí)船驗(yàn)證
    中國航海(2017年3期)2017-11-03 02:46:54
    船模靜水橫搖試驗(yàn)的不確定度分析
    基于三維水動(dòng)力方法的三體船連接橋波浪載荷計(jì)算分析
    西洋船模王——童鑑良
    航海(2014年6期)2014-12-12 10:36:03
    三體船側(cè)體位置優(yōu)化設(shè)計(jì)研究
    船舶(2012年4期)2012-08-11 01:48:36
    基于VDMS系統(tǒng)的通航船模操縱性率定研究
    久久中文看片网| 欧美黄色片欧美黄色片| 欧美日韩国产mv在线观看视频| 国产精品亚洲av一区麻豆| a在线观看视频网站| 国产深夜福利视频在线观看| netflix在线观看网站| 久久精品久久久久久噜噜老黄| 欧美日韩中文字幕国产精品一区二区三区 | 91精品三级在线观看| 老司机靠b影院| 丰满少妇做爰视频| 国产精品 欧美亚洲| 欧美精品一区二区免费开放| 中文字幕人妻丝袜一区二区| 丝瓜视频免费看黄片| 不卡av一区二区三区| 搡老乐熟女国产| 亚洲国产av新网站| 久久免费观看电影| 国产主播在线观看一区二区| 热re99久久精品国产66热6| 五月天丁香电影| 国产男女内射视频| 黄色视频,在线免费观看| 纯流量卡能插随身wifi吗| 亚洲三区欧美一区| 97人妻天天添夜夜摸| 久久国产亚洲av麻豆专区| 日本猛色少妇xxxxx猛交久久| 免费不卡黄色视频| 亚洲av男天堂| 欧美在线一区亚洲| 啦啦啦免费观看视频1| 午夜福利影视在线免费观看| 人人妻人人澡人人看| 久久久欧美国产精品| 亚洲少妇的诱惑av| 91字幕亚洲| 咕卡用的链子| 亚洲伊人久久精品综合| 亚洲人成电影观看| 老汉色∧v一级毛片| 亚洲av片天天在线观看| 国产片内射在线| 久9热在线精品视频| 免费人妻精品一区二区三区视频| a在线观看视频网站| 最近最新免费中文字幕在线| kizo精华| 极品少妇高潮喷水抽搐| 欧美日韩黄片免| av不卡在线播放| 欧美日韩一级在线毛片| 亚洲国产精品一区二区三区在线| 亚洲精品美女久久av网站| 日韩视频在线欧美| 好男人电影高清在线观看| 成年人午夜在线观看视频| 亚洲情色 制服丝袜| 菩萨蛮人人尽说江南好唐韦庄| 一个人免费在线观看的高清视频 | 啦啦啦免费观看视频1| 少妇粗大呻吟视频| 亚洲国产精品成人久久小说| 亚洲中文日韩欧美视频| 国产国语露脸激情在线看| 人人妻人人澡人人爽人人夜夜| 在线精品无人区一区二区三| tube8黄色片| 午夜视频精品福利| 黄色怎么调成土黄色| 国产精品一区二区免费欧美 | 十八禁高潮呻吟视频| 黄色视频在线播放观看不卡| 日日摸夜夜添夜夜添小说| 黄色怎么调成土黄色| 精品少妇内射三级| 青春草亚洲视频在线观看| 欧美激情极品国产一区二区三区| 热99国产精品久久久久久7| 人妻久久中文字幕网| 三级毛片av免费| 精品久久蜜臀av无| 一区二区av电影网| 久久中文字幕一级| 国产一区二区三区在线臀色熟女 | 日韩一区二区三区影片| 欧美日韩国产mv在线观看视频| 一级毛片精品| 国产免费一区二区三区四区乱码| 久久99热这里只频精品6学生| 免费观看a级毛片全部| 亚洲av电影在线进入| 女人久久www免费人成看片| 国内毛片毛片毛片毛片毛片| 精品亚洲成a人片在线观看| 国产精品一二三区在线看| 水蜜桃什么品种好| 五月天丁香电影| 欧美性长视频在线观看| 高清欧美精品videossex| 高潮久久久久久久久久久不卡| 亚洲欧洲日产国产| 天天操日日干夜夜撸| 亚洲精品美女久久av网站| 18禁国产床啪视频网站| 一区二区日韩欧美中文字幕| 国产野战对白在线观看| 飞空精品影院首页| 另类精品久久| 又紧又爽又黄一区二区| 777久久人妻少妇嫩草av网站| 91精品国产国语对白视频| 欧美另类亚洲清纯唯美| 老司机午夜福利在线观看视频 | 性色av一级| 亚洲专区字幕在线| 国产成+人综合+亚洲专区| 老司机亚洲免费影院| 久久精品国产综合久久久| 丁香六月欧美| netflix在线观看网站| 精品第一国产精品| 日日夜夜操网爽| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产欧美一区二区综合| 99精国产麻豆久久婷婷| 男女免费视频国产| 黑丝袜美女国产一区| 80岁老熟妇乱子伦牲交| 国产欧美日韩精品亚洲av| 国产欧美日韩一区二区三 | 波多野结衣av一区二区av| 999久久久国产精品视频| 欧美日韩亚洲高清精品| av在线app专区| 蜜桃在线观看..| 国产一区二区激情短视频 | 中国美女看黄片| 满18在线观看网站| 黄频高清免费视频| 国产精品成人在线| 90打野战视频偷拍视频| 亚洲av美国av| 十八禁网站网址无遮挡| 国产一区二区激情短视频 | 久久人人爽av亚洲精品天堂| 91大片在线观看| 亚洲国产av新网站| 日韩一卡2卡3卡4卡2021年| 三级毛片av免费| 中文欧美无线码| 亚洲精品日韩在线中文字幕| 美女大奶头黄色视频| a级毛片在线看网站| 亚洲九九香蕉| 国产伦人伦偷精品视频| 满18在线观看网站| 69av精品久久久久久 | 亚洲国产日韩一区二区| 日韩三级视频一区二区三区| 亚洲精品美女久久久久99蜜臀| 十八禁网站免费在线| 99国产精品免费福利视频| 日韩欧美一区视频在线观看| 老熟女久久久| 天天躁狠狠躁夜夜躁狠狠躁| 精品一品国产午夜福利视频| 捣出白浆h1v1| 亚洲第一av免费看| 一个人免费看片子| 日韩精品免费视频一区二区三区| 在线精品无人区一区二区三| 免费观看人在逋| 麻豆av在线久日| 人成视频在线观看免费观看| 美女大奶头黄色视频| 国产精品秋霞免费鲁丝片| 免费久久久久久久精品成人欧美视频| 日韩大片免费观看网站| 欧美xxⅹ黑人| 黑丝袜美女国产一区| 91成年电影在线观看| 久久精品成人免费网站| 最黄视频免费看| 天天添夜夜摸| 俄罗斯特黄特色一大片| 亚洲精品自拍成人| 久久99一区二区三区| 高清黄色对白视频在线免费看| 天堂中文最新版在线下载| 亚洲性夜色夜夜综合| 视频在线观看一区二区三区| 首页视频小说图片口味搜索| 啦啦啦在线免费观看视频4| 真人做人爱边吃奶动态| 2018国产大陆天天弄谢| 伦理电影免费视频| 99九九在线精品视频| 激情视频va一区二区三区| 国产成+人综合+亚洲专区| 老司机亚洲免费影院| 秋霞在线观看毛片| 成年av动漫网址| 免费黄频网站在线观看国产| 午夜福利一区二区在线看| 欧美日本中文国产一区发布| 亚洲久久久国产精品| 叶爱在线成人免费视频播放| 俄罗斯特黄特色一大片| 国产一区二区三区在线臀色熟女 | videosex国产| 亚洲欧美一区二区三区黑人| 亚洲七黄色美女视频| cao死你这个sao货| 日韩大码丰满熟妇| 欧美日韩福利视频一区二区| 亚洲av男天堂| svipshipincom国产片| 国产高清国产精品国产三级| 一级毛片女人18水好多| 亚洲国产欧美网| 欧美国产精品va在线观看不卡| 97精品久久久久久久久久精品| 动漫黄色视频在线观看| 乱人伦中国视频| 日韩 欧美 亚洲 中文字幕| a级片在线免费高清观看视频| 精品亚洲乱码少妇综合久久| 久久国产精品男人的天堂亚洲| 男女床上黄色一级片免费看| 久久精品成人免费网站| 久久性视频一级片| 中文字幕人妻熟女乱码| 国产麻豆69| 制服人妻中文乱码| 天天躁日日躁夜夜躁夜夜| 麻豆乱淫一区二区| 国产精品久久久久久精品古装| 久久精品国产a三级三级三级| 男女床上黄色一级片免费看| 国产福利在线免费观看视频| 一进一出抽搐动态| 日韩有码中文字幕| 亚洲色图 男人天堂 中文字幕| 人人妻人人澡人人看| 国产精品欧美亚洲77777| 午夜福利视频精品| 另类精品久久| 精品视频人人做人人爽| 91字幕亚洲| 美女中出高潮动态图| 亚洲av电影在线进入| 免费高清在线观看日韩| 人成视频在线观看免费观看| 久久香蕉激情| 真人做人爱边吃奶动态| 国产精品一区二区在线不卡| 夜夜骑夜夜射夜夜干| 久久久久久久大尺度免费视频| 新久久久久国产一级毛片| 久久久精品免费免费高清| 中文字幕人妻丝袜一区二区| 伊人亚洲综合成人网| 亚洲综合色网址| 欧美人与性动交α欧美软件| 婷婷成人精品国产| 黑人巨大精品欧美一区二区蜜桃| 另类亚洲欧美激情| 国内毛片毛片毛片毛片毛片| 别揉我奶头~嗯~啊~动态视频 | 中文字幕制服av| 男人添女人高潮全过程视频| 一区二区日韩欧美中文字幕| 日韩视频一区二区在线观看| 国产一区二区在线观看av| 丰满迷人的少妇在线观看| 王馨瑶露胸无遮挡在线观看| 国产亚洲av高清不卡| 午夜福利视频在线观看免费| 99精品欧美一区二区三区四区| 男女之事视频高清在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜久久久在线观看| 在线观看免费高清a一片| 一级黄色大片毛片| 在线观看人妻少妇| 99re6热这里在线精品视频| 国内毛片毛片毛片毛片毛片| 91国产中文字幕| av网站在线播放免费| 亚洲国产欧美网| 亚洲成人免费电影在线观看| 国产无遮挡羞羞视频在线观看| 国产xxxxx性猛交| 精品乱码久久久久久99久播| 天天操日日干夜夜撸| 爱豆传媒免费全集在线观看| 啦啦啦视频在线资源免费观看| 国产日韩欧美视频二区| 国产精品香港三级国产av潘金莲| 啪啪无遮挡十八禁网站| 美女扒开内裤让男人捅视频| 五月天丁香电影| 男人舔女人的私密视频| 国产精品久久久av美女十八| kizo精华| 免费观看av网站的网址| 精品国产一区二区久久| 大片免费播放器 马上看| 99国产极品粉嫩在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品麻豆人妻色哟哟久久| 成年动漫av网址| 国产亚洲精品久久久久5区| 免费观看av网站的网址| 99久久国产精品久久久| 丝袜美足系列| www.自偷自拍.com| 性色av乱码一区二区三区2| 首页视频小说图片口味搜索| 色视频在线一区二区三区| 9热在线视频观看99| 在线观看免费视频网站a站| 国产黄频视频在线观看| 成人影院久久| 国产成人精品久久二区二区免费| 18禁观看日本| 精品久久久久久久毛片微露脸 | 桃花免费在线播放| 亚洲第一欧美日韩一区二区三区 | 久久九九热精品免费| 国产日韩一区二区三区精品不卡| 久久久久国产精品人妻一区二区| 人妻 亚洲 视频| a级毛片在线看网站| 国产亚洲午夜精品一区二区久久| 欧美日韩av久久| av在线app专区| 51午夜福利影视在线观看| 狂野欧美激情性xxxx| 国产av国产精品国产| 999久久久精品免费观看国产| 国产免费av片在线观看野外av| 国产福利在线免费观看视频| 亚洲精品乱久久久久久| 丰满少妇做爰视频| 男人舔女人的私密视频| 国产欧美亚洲国产| 制服诱惑二区| 人人澡人人妻人| 久久久精品区二区三区| 欧美午夜高清在线| 精品福利永久在线观看| 大片电影免费在线观看免费| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲欧美色中文字幕在线| 啦啦啦 在线观看视频| 少妇人妻久久综合中文| 色老头精品视频在线观看| 久久久久久人人人人人| 欧美黑人精品巨大| 久久人妻福利社区极品人妻图片| 国产一区二区三区av在线| 日本av手机在线免费观看| 国产一级毛片在线| 叶爱在线成人免费视频播放| 日本wwww免费看| 大陆偷拍与自拍| 熟女少妇亚洲综合色aaa.| 亚洲五月色婷婷综合| 美女高潮到喷水免费观看| 午夜精品久久久久久毛片777| 9色porny在线观看| 曰老女人黄片| 12—13女人毛片做爰片一| avwww免费| 亚洲av日韩在线播放| 国产高清国产精品国产三级| 超碰97精品在线观看| 国产男女超爽视频在线观看| 免费久久久久久久精品成人欧美视频| 中文字幕精品免费在线观看视频| 女人被躁到高潮嗷嗷叫费观| 国产一区二区激情短视频 | 亚洲精品国产一区二区精华液| av网站在线播放免费| 国产成+人综合+亚洲专区| 欧美+亚洲+日韩+国产| 免费少妇av软件| 国产福利在线免费观看视频| 50天的宝宝边吃奶边哭怎么回事| 亚洲伊人色综图| cao死你这个sao货| 精品少妇久久久久久888优播| 午夜精品国产一区二区电影| 国产在线视频一区二区| 在线观看免费日韩欧美大片| 国产精品一区二区在线不卡| 91字幕亚洲| 久久国产精品大桥未久av| 多毛熟女@视频| 久久久久国产精品人妻一区二区| 国产不卡av网站在线观看| 嫁个100分男人电影在线观看| 久久人人爽av亚洲精品天堂| 成年av动漫网址| 男女之事视频高清在线观看| 久久久久精品国产欧美久久久 | 免费黄频网站在线观看国产| 一级片免费观看大全| 男女国产视频网站| 在线观看www视频免费| 丝袜美腿诱惑在线| 成人18禁高潮啪啪吃奶动态图| 狠狠精品人妻久久久久久综合| 欧美日韩亚洲综合一区二区三区_| 在线av久久热| 美国免费a级毛片| 精品视频人人做人人爽| 丝袜人妻中文字幕| 免费在线观看日本一区| 91字幕亚洲| 51午夜福利影视在线观看| 久久久精品国产亚洲av高清涩受| 啦啦啦在线免费观看视频4| 啦啦啦免费观看视频1| 国产欧美日韩一区二区精品| 国产在视频线精品| 国产主播在线观看一区二区| 色精品久久人妻99蜜桃| 黑人欧美特级aaaaaa片| 在线亚洲精品国产二区图片欧美| 90打野战视频偷拍视频| 国产成人啪精品午夜网站| 亚洲国产欧美日韩在线播放| 最近最新中文字幕大全免费视频| 脱女人内裤的视频| 女人被躁到高潮嗷嗷叫费观| 黄色视频,在线免费观看| 777米奇影视久久| 啦啦啦啦在线视频资源| 日韩视频一区二区在线观看| 成人影院久久| 波多野结衣av一区二区av| 久久免费观看电影| 久久久久久久久免费视频了| 亚洲精品第二区| 国产主播在线观看一区二区| 久久精品亚洲av国产电影网| 日韩欧美免费精品| 免费高清在线观看视频在线观看| 天天躁日日躁夜夜躁夜夜| 国产真人三级小视频在线观看| 麻豆国产av国片精品| 国产精品久久久久成人av| 另类亚洲欧美激情| 国产精品亚洲av一区麻豆| 成人av一区二区三区在线看 | 国产一区二区 视频在线| 亚洲第一av免费看| 99久久国产精品久久久| 久久久久网色| 亚洲精品美女久久av网站| 自拍欧美九色日韩亚洲蝌蚪91| 美女视频免费永久观看网站| 在线看a的网站| 丁香六月欧美| 少妇被粗大的猛进出69影院| 老司机午夜十八禁免费视频| 亚洲精品第二区| 久久中文字幕一级| 高潮久久久久久久久久久不卡| 国产精品麻豆人妻色哟哟久久| 男人操女人黄网站| 欧美人与性动交α欧美软件| 中文字幕人妻丝袜一区二区| 欧美日韩亚洲国产一区二区在线观看 | 高潮久久久久久久久久久不卡| 午夜福利视频在线观看免费| 欧美日韩成人在线一区二区| 人妻人人澡人人爽人人| 另类亚洲欧美激情| 菩萨蛮人人尽说江南好唐韦庄| 一级a爱视频在线免费观看| kizo精华| 国产欧美日韩综合在线一区二区| 王馨瑶露胸无遮挡在线观看| av一本久久久久| 水蜜桃什么品种好| 亚洲成av片中文字幕在线观看| 视频区图区小说| 亚洲欧美精品综合一区二区三区| 久久久水蜜桃国产精品网| 亚洲伊人久久精品综合| 一级黄色大片毛片| 午夜福利在线免费观看网站| 黑人猛操日本美女一级片| 丝瓜视频免费看黄片| 欧美精品人与动牲交sv欧美| 国产精品一区二区免费欧美 | 十八禁网站免费在线| 十分钟在线观看高清视频www| a 毛片基地| 美女高潮到喷水免费观看| 精品一区在线观看国产| 99久久人妻综合| 亚洲美女黄色视频免费看| 老司机在亚洲福利影院| 人人妻人人澡人人爽人人夜夜| 久久久久精品国产欧美久久久 | 亚洲国产精品999| 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利视频在线观看免费| 黄片大片在线免费观看| 免费日韩欧美在线观看| 欧美黑人欧美精品刺激| 亚洲avbb在线观看| av线在线观看网站| 国产成人影院久久av| 欧美精品亚洲一区二区| 美女脱内裤让男人舔精品视频| av国产精品久久久久影院| 99精品欧美一区二区三区四区| 精品一区在线观看国产| 成年动漫av网址| 可以免费在线观看a视频的电影网站| 日韩欧美国产一区二区入口| 亚洲国产欧美一区二区综合| av不卡在线播放| 国内毛片毛片毛片毛片毛片| 国产成人a∨麻豆精品| 国产91精品成人一区二区三区 | 汤姆久久久久久久影院中文字幕| 久久久久网色| av有码第一页| 高潮久久久久久久久久久不卡| 久久狼人影院| 婷婷丁香在线五月| 在线观看人妻少妇| 日本av免费视频播放| 精品福利观看| 久久99一区二区三区| 欧美性长视频在线观看| 99国产综合亚洲精品| 国产成人a∨麻豆精品| 国产在线视频一区二区| 2018国产大陆天天弄谢| 日韩 欧美 亚洲 中文字幕| 最新在线观看一区二区三区| 国产精品影院久久| 国产一区二区在线观看av| 亚洲全国av大片| 成人手机av| 久久久久国产精品人妻一区二区| 最近最新中文字幕大全免费视频| 亚洲欧美清纯卡通| 老司机影院毛片| 91精品伊人久久大香线蕉| 亚洲精品美女久久久久99蜜臀| 久久久国产一区二区| 亚洲国产毛片av蜜桃av| 午夜老司机福利片| 成人国产一区最新在线观看| 女警被强在线播放| 午夜视频精品福利| 美女主播在线视频| 一级a爱视频在线免费观看| 亚洲精品中文字幕一二三四区 | 午夜老司机福利片| 悠悠久久av| 青草久久国产| 激情视频va一区二区三区| 丝袜人妻中文字幕| 亚洲五月色婷婷综合| 麻豆av在线久日| 一二三四在线观看免费中文在| 纵有疾风起免费观看全集完整版| 亚洲va日本ⅴa欧美va伊人久久 | 美女主播在线视频| 每晚都被弄得嗷嗷叫到高潮| 制服诱惑二区| 久久久久久久国产电影| 性少妇av在线| 国产高清视频在线播放一区 | 久久久久久人人人人人| 国产熟女午夜一区二区三区| 中文字幕制服av| 国产精品秋霞免费鲁丝片| 日本五十路高清| 可以免费在线观看a视频的电影网站| 国产在线观看jvid| 亚洲avbb在线观看| 热99久久久久精品小说推荐| 久久久久精品国产欧美久久久 | 美女脱内裤让男人舔精品视频| 王馨瑶露胸无遮挡在线观看| 国产精品av久久久久免费| 欧美人与性动交α欧美软件| 欧美在线一区亚洲| 国产视频一区二区在线看| 成年人黄色毛片网站| 悠悠久久av| 老司机午夜十八禁免费视频| av在线老鸭窝| 国产男女内射视频| 欧美变态另类bdsm刘玥| 国产人伦9x9x在线观看| 国产91精品成人一区二区三区 | 制服人妻中文乱码| 日韩三级视频一区二区三区| 精品卡一卡二卡四卡免费| 国产av精品麻豆| 老司机亚洲免费影院|