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

    雙斷級滑行艇的阻力回歸公式

    2010-03-06 03:22:52沈小紅吳啟銳鄧愛民
    中國艦船研究 2010年2期
    關(guān)鍵詞:計算公式回歸方程阻力

    沈小紅 吳啟銳 肖 笛 鄧愛民

    中國艦船研究設(shè)計中心,湖北 武 漢 430064

    雙斷級滑行艇的阻力回歸公式

    沈小紅 吳啟銳 肖 笛 鄧愛民

    中國艦船研究設(shè)計中心,湖北 武 漢 430064

    本文以5艘雙斷級滑行艇模型的試驗(yàn)數(shù)據(jù)為基礎(chǔ),采用多元線性回歸的方法對試驗(yàn)數(shù)據(jù)進(jìn)行回歸分析,從而導(dǎo)出在不同斷級參數(shù)情況下雙斷級滑行艇的阻力回歸公式。通過相對誤差分析法驗(yàn)證了回歸公式的精確性,并根據(jù)回歸的結(jié)果繪制成圖以供設(shè)計參考。

    雙斷級滑行艇;多元線性回歸;阻力回歸公式;相對誤差分析法

    1 引言

    雙斷級滑行艇的主要優(yōu)點(diǎn)是航速高、阻力小,而兩個斷級的幾何參數(shù)(斷級的縱向位置和高度)對其阻力性能的影響至關(guān)重要,由于各滑行面之間的干擾比較復(fù)雜[1],因而斷級的確定主要還是依靠模型試驗(yàn),至今尚無一種有效的理論計算方法。本文將根據(jù)模型試驗(yàn)所得到的試驗(yàn)數(shù)據(jù),用多項(xiàng)式回歸的方法[2-3]得到在不同斷級參數(shù)情況下雙斷級滑行艇的阻力公式。

    本文將取5艘雙斷級滑行艇模型的試驗(yàn)數(shù)據(jù)進(jìn)行回歸分析,除了斷級參數(shù)以外,這些模型的主尺度完全相同。這些模型的編號及斷級參數(shù)值如表1所示。

    表中,l1為斷級總長度;θ1為斷級角度 (任何縱剖線上兩個斷級下緣點(diǎn)與艇尾封板下緣點(diǎn)的連線在一條直線上[4],此直線與基線的夾角即為斷級角度);l2為前斷級長度;l3為后斷級長度;h為斷級高度(兩個斷級高度相等)。通過斷級總長度l1和斷級角度θ1這兩個因素應(yīng)用幾何關(guān)系可求出前斷級長度l2、后斷級長度l3和斷級高度h的值,這樣就能將雙斷級艇的斷級參數(shù)從3個因素的變化轉(zhuǎn)化為兩個因素的變化。因此,本文所選取的試驗(yàn)?zāi)P偷臄嗉墔?shù)就是由斷級總長度l1和斷級角度θ1這兩個因素來確定的。

    表1 模型的編號及對應(yīng)的斷級參數(shù)值

    2 多元線性回歸

    設(shè)隨機(jī)變量y與非隨機(jī)變量之間有如下的p元線性回歸模型:

    式中,β0、 β1、 …、 βp是p+1 個待估參數(shù);x1、 x2、…、xp是p個可以精確測量或可控制的一般變量;ε1、 ε2、…、 εN是 N 個相互獨(dú)立且服從同一正態(tài)分布N(0,σ)的隨機(jī)變量。這就是多元線性回歸的數(shù)學(xué)模型[5]。

    采用最小二乘估計法可求出系數(shù)β0、β1、…和βp,而后可確定 y 與 x1、 x2、…、 xp之間的關(guān)系為:

    y=β0+β1x1+ β2x2+ …+ βpxp(2)

    在上述過程中,尚需討論兩個方面的問題:一方面是線性相關(guān)關(guān)系的假設(shè)檢驗(yàn),以驗(yàn)證y與x1、x2、…、xp之間確實(shí)存在線性關(guān)系;另一方面是回歸系數(shù)的假檢驗(yàn),以驗(yàn)證每個因子x1、x2、…、xp對y的影響程度,剔除那些影響不顯著的因子,這就需要選擇“最優(yōu)”的回歸方程[6]。

    所謂“最優(yōu)”回歸方程,就是包含所有對y影響顯著的變量而不包含對y影響不顯著的變量的回歸方程。選擇“最優(yōu)”回歸方程的方法有幾種,本文將采用逐步回歸分析的方法來建立“最優(yōu)”回歸方程。

    逐步回歸分析的基本思想是將因子逐個引入。引入因子的條件是,該因子的偏回歸平方和經(jīng)檢驗(yàn)是顯著的。同時,每引入一個新因子后,要對舊因子逐個檢驗(yàn),將偏回歸平方和變?yōu)椴伙@著的因子剔除。采用上述方法,由于每一步都進(jìn)行檢驗(yàn),因而保證了最后所得的方程中所有的因子都是顯著的。

    3 雙斷級滑行艇的阻力回歸公式

    下面將模型A1、A2和A3(這些模型具有相同的斷級總長度l1)作為一組進(jìn)行回歸,再將模型A2、B2和C2(這些模型具有相同的斷級夾角θ1)作為一組進(jìn)行回歸。

    3.1 A1、A2 和 A3 阻力回歸公式

    在回歸分析前,首先要將因變量和所有自變量進(jìn)行無因次化[7],各變量無因次化結(jié)果如下:

    阻力無因次化:阻升比R/Δ;重心無因次化:xg/B(B為折角線平均寬度);斷級長度 l1無因次化:l1/B。

    斷級夾角θ1在實(shí)艇和模型中具有相同的角度,它們的幾何關(guān)系不會發(fā)生變化,故斷級夾角θ1無需進(jìn)行無因次化。

    以斷級夾角θ1、xg/B和動負(fù)荷系數(shù)CB為自變量,阻升比R/Δ為因變量,進(jìn)行三元多項(xiàng)式回歸分析[8],這樣可以把阻升比R/Δ的試驗(yàn)數(shù)據(jù)在斷級夾角θ1、xg/B和動負(fù)荷系數(shù)CB3個方向上進(jìn)行光順。故可設(shè)回歸多項(xiàng)式的表達(dá)式如下:

    從繪出的試驗(yàn)點(diǎn)草圖可以看出試驗(yàn)曲線比較平緩,沒有出現(xiàn)拐點(diǎn),則可取 n1=n2=n3=2,即可將式(3)轉(zhuǎn)化為:

    顯然,式(4)為一非線性表達(dá)式,為了能應(yīng)用回歸的方法進(jìn)行處理,可以采用變量替換的方法將其轉(zhuǎn)化為多元線性回歸表達(dá)式[9]。

    本文所取試驗(yàn)數(shù)據(jù)為:斷級夾角θ1取0.81°、1.11°、1.41°等 3 個值; xg/B 取 1.062 5、1.125 0、1.187 5 等 3 個值;動負(fù)荷系數(shù) CB取 0.030 16、0.032 18、0.033 42、0.034 19、0.035 65、0.037 88、0.039 89、0.042 55、0.045 21、0.047 48、0.050 64、0.053 81、0.057 44、0.061 27、0.065 10 等 15 個值。形成的阻升比數(shù)據(jù)如表2所示。

    表2 A1、A2和A3阻升比R/Δ 數(shù)據(jù)

    根據(jù)逐步回歸的計算方法,應(yīng)用Matlab編程計算可算得顯著變量回歸系數(shù)Lijk的值,其具體結(jié)果如表3。

    表3 A1、A2和A3阻升比R/Δ 計算公式的Lijk值

    3.2 A2、B2 和 C2 阻力回歸公式

    按照回歸A1、A2和A3阻力公式的方法,同樣可以對A2、B2和C2阻力公式進(jìn)行回歸,故可設(shè)其回歸多項(xiàng)式的表達(dá)式為:

    采用變量替換的方法將式(6)轉(zhuǎn)化為多元線性回歸表達(dá)式的方法同上節(jié)內(nèi)容,故在此不再贅述。 在本節(jié)所取的試驗(yàn)數(shù)據(jù)中,l1/B 取 1.481 25、1.575 0、1.668 75 等 3 個值;xg/B 和動負(fù)荷系數(shù) CB所取的值同上一節(jié)。同樣也形成了個數(shù)據(jù),其阻升比R/Δ的數(shù)據(jù)如表4所示,應(yīng)用Matlab編程計算所得的阻升比R/Δ計算公式的Lijk值如表5所示。

    4 計算公式精確性分析

    對于以上回歸的阻升比計算公式,還應(yīng)分析其精確性。本文將采用相對誤差分析法[10],從兩個方面來驗(yàn)證以上回歸公式的精確性。首先從進(jìn)行回歸的試驗(yàn)數(shù)據(jù)中取出原始數(shù)據(jù)進(jìn)行比較 (由于試驗(yàn)數(shù)據(jù)太多,故只取小部分?jǐn)?shù)據(jù)進(jìn)行比較驗(yàn)證),其比較結(jié)果如表6和表7所示;其次,由于模型 A3 在重心縱向位置 xg=640 mm(即 xg/B=1.0)時的試驗(yàn)數(shù)據(jù)并沒有參加回歸,則可用此試驗(yàn)數(shù)據(jù)與回歸公式計算結(jié)果進(jìn)行比較,其比較結(jié)果如表8所示。

    由以上分析結(jié)果可見,根據(jù)回歸計算公式所求得結(jié)果與原始試驗(yàn)值非常接近。此外,如圖1~5所示,回歸曲線與試驗(yàn)數(shù)值點(diǎn)具有良好的貼近性,由此,可以看出本文的回歸公式計算精度比較高。

    表4 A2、B2和C2阻升比R/Δ 數(shù)據(jù)

    表5 A2、B2和C2阻升比R/Δ 計算公式的Lijk值

    表 6 A1、A2 和 A3 在 xg/B = 1.575 時比較結(jié)果

    表 7 A2、B2 和 C2 在 xg/B = 1.575 時比較結(jié)果

    表 8 A3 在 xg/B=1.0 時比較結(jié)果

    5 結(jié) 論

    本文采用多元線性回歸中逐步回歸的計算方法,得出雙斷級滑行艇在不同斷級參數(shù)下的 “最優(yōu)”阻力回歸公式,通過對計算值與試驗(yàn)值的相對誤差分析,可以看出回歸所得的計算公式精度比較高,從而證明了根據(jù)試驗(yàn)數(shù)據(jù),采用回歸分析的方法得出雙斷級滑行艇的阻力計算公式是有效而可信的。本文所得的回歸公式可以為設(shè)計雙斷級滑行艇的斷級參數(shù)選擇提供一些理論依據(jù)。

    [1]董祖舜.快艇動力學(xué)[M].武漢:華中理工大學(xué)出版社出版,1991,4.

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

    [3]譚仲楷.應(yīng)用回歸分析方法估算圓艇舭阻力[J].中國造船,1979,67(2):14-29.

    [4]夏翔.超高速斷級滑行艇推進(jìn)系統(tǒng)工況配合特性的研究[J].中國艦船研究,2008,3(2):39-42.

    [5]上海師范大學(xué)數(shù)學(xué)系.回歸分析及其試驗(yàn)設(shè)計 [M].上海:上海教育出版社,1978.

    [6]何志飛,邢圣德.滑行艇之最佳阻力[J].江蘇船舶,1996,13(4):1-4.

    [7]郭亨翔,倪家俊,張蓓玉.當(dāng)今世界優(yōu)秀船型主要參數(shù)的回歸分析[J].船舶設(shè)計通訊,2002(4):12-31.

    [8] LAMMEREN W P A VAN,MANEN J D VAN,OOSTERVELD M W C.The wageningen B-screw series Trans[J].Society of Naval Architects and Marine Engineer,1969(77).

    [9]陳賓康,趙成璧,洪安祥.船用組合舵水動力系數(shù)的回歸公式[J].中國造船,1999,147(4)20-27.

    [10]李紅.數(shù)值分析[M].武漢:華中科技大學(xué)出版社,2003.

    A Resistance Regression Formula for Double-stepped Planing Craft

    Shen Xiao-h(huán)ong Wu Qi-ruiXiao DiDeng Ai-min
    China Ship Development and Design Center, Wuhan 430064, China

    Based on the model test data of five double-stepped planing crafts, the multiple linear regression method was used to treat the test data and thus a resistance regression formula was obtained for the craft with different parameters of its double steps.Then, the accuracy of the regression formula was verified through the relative error analysis.Finally,the regression results were plotted in charts so as to be references for design.

    double-stepped planing craft; multi-factor linear regression; resistance regression formula;relative-error analysis method

    U661.3

    A

    1673-3185(2010)02-59-05

    2009-04-16

    沈小紅(1980-),女,碩士研究生。研究方向:船舶與海洋結(jié)構(gòu)物設(shè)計制造。E-mail:shenxiaohong304@163.com

    吳啟銳(1965-),男,研究員,碩士生導(dǎo)師。研究方向:船舶與海洋結(jié)構(gòu)物設(shè)計制造

    猜你喜歡
    計算公式回歸方程阻力
    電機(jī)溫升計算公式的推導(dǎo)和應(yīng)用
    鼻阻力測定在兒童OSA診療中的臨床作用
    采用直線回歸方程預(yù)測桑癭蚊防治適期
    線性回歸方程的求解與應(yīng)用
    線性回歸方程要點(diǎn)導(dǎo)學(xué)
    零阻力
    英語文摘(2020年10期)2020-11-26 08:12:12
    2019離職補(bǔ)償金計算公式一覽表
    走進(jìn)回歸分析,讓回歸方程不再是你高考的絆腳石
    別讓摩擦成為學(xué)習(xí)的阻力
    阻力不小 推進(jìn)當(dāng)循序漸進(jìn)
    一本色道久久久久久精品综合| 99久国产av精品国产电影| 国产片特级美女逼逼视频| 韩国高清视频一区二区三区| 国产成人系列免费观看| 在线精品无人区一区二区三| 在线看a的网站| 一区二区三区精品91| 欧美日韩福利视频一区二区| 啦啦啦在线观看免费高清www| av国产精品久久久久影院| 又黄又粗又硬又大视频| 熟妇人妻不卡中文字幕| 欧美激情 高清一区二区三区| 操美女的视频在线观看| 男女边吃奶边做爰视频| 午夜91福利影院| 亚洲欧美一区二区三区国产| av在线老鸭窝| 国产又爽黄色视频| 伦理电影大哥的女人| 国产女主播在线喷水免费视频网站| 国产熟女午夜一区二区三区| 成人国产麻豆网| 啦啦啦在线观看免费高清www| 熟女少妇亚洲综合色aaa.| 中文字幕最新亚洲高清| 曰老女人黄片| 午夜福利视频在线观看免费| 午夜福利免费观看在线| 两性夫妻黄色片| 婷婷色综合大香蕉| 青春草亚洲视频在线观看| 国产精品嫩草影院av在线观看| 亚洲国产欧美日韩在线播放| 久久精品aⅴ一区二区三区四区| 成人免费观看视频高清| 观看av在线不卡| 日韩精品免费视频一区二区三区| 在线观看免费视频网站a站| 九草在线视频观看| 午夜福利免费观看在线| 别揉我奶头~嗯~啊~动态视频 | 成年动漫av网址| 美女高潮到喷水免费观看| 欧美精品人与动牲交sv欧美| 久久精品久久久久久噜噜老黄| 夫妻午夜视频| 五月天丁香电影| 18禁观看日本| 熟女av电影| 国产有黄有色有爽视频| 又大又黄又爽视频免费| 91老司机精品| 亚洲精品第二区| 亚洲精品一区蜜桃| 国产一卡二卡三卡精品 | 操出白浆在线播放| 午夜久久久在线观看| 精品卡一卡二卡四卡免费| 99热全是精品| av视频免费观看在线观看| 欧美在线一区亚洲| 国产一区二区在线观看av| 亚洲国产精品一区三区| a 毛片基地| 亚洲av在线观看美女高潮| 亚洲人成网站在线观看播放| 成年美女黄网站色视频大全免费| 欧美日韩福利视频一区二区| 人人妻人人添人人爽欧美一区卜| 欧美成人午夜精品| 免费黄频网站在线观看国产| 日韩人妻精品一区2区三区| 色吧在线观看| 男男h啪啪无遮挡| 亚洲五月色婷婷综合| 涩涩av久久男人的天堂| 午夜福利免费观看在线| 99热国产这里只有精品6| 男女国产视频网站| 水蜜桃什么品种好| 两个人免费观看高清视频| 欧美日韩综合久久久久久| 欧美精品一区二区免费开放| 精品第一国产精品| 99热全是精品| 国产伦理片在线播放av一区| 又大又黄又爽视频免费| 亚洲成人手机| 精品亚洲乱码少妇综合久久| 九草在线视频观看| 免费在线观看黄色视频的| 成人影院久久| 久久久久网色| 国产色婷婷99| 国产伦理片在线播放av一区| 中文天堂在线官网| xxx大片免费视频| 精品国产国语对白av| 国产精品久久久久久久久免| 久久亚洲国产成人精品v| 天天躁夜夜躁狠狠久久av| 亚洲色图 男人天堂 中文字幕| www.熟女人妻精品国产| 国产精品熟女久久久久浪| 亚洲精品久久成人aⅴ小说| 美女脱内裤让男人舔精品视频| 精品国产乱码久久久久久小说| 国产男女内射视频| 男男h啪啪无遮挡| 老汉色∧v一级毛片| 国产成人午夜福利电影在线观看| 亚洲欧美精品自产自拍| 久久久久国产一级毛片高清牌| 精品卡一卡二卡四卡免费| 如何舔出高潮| 欧美激情极品国产一区二区三区| 国产精品一国产av| netflix在线观看网站| 中文字幕制服av| 免费在线观看完整版高清| 色精品久久人妻99蜜桃| av天堂久久9| 美女国产高潮福利片在线看| 欧美97在线视频| 啦啦啦 在线观看视频| 街头女战士在线观看网站| 久久久精品94久久精品| 久久久精品区二区三区| 久久精品国产亚洲av涩爱| 午夜日本视频在线| 国产精品 欧美亚洲| 亚洲四区av| 1024香蕉在线观看| 叶爱在线成人免费视频播放| 91精品三级在线观看| a级毛片在线看网站| 成人国语在线视频| 午夜久久久在线观看| 在现免费观看毛片| 日本欧美视频一区| 满18在线观看网站| 国产男人的电影天堂91| 热re99久久国产66热| 日本91视频免费播放| 欧美日韩综合久久久久久| 精品国产一区二区久久| 91精品伊人久久大香线蕉| 国产精品久久久久久久久免| 五月开心婷婷网| 男人操女人黄网站| 精品一区二区三区av网在线观看 | 极品人妻少妇av视频| 老司机深夜福利视频在线观看 | 亚洲国产精品成人久久小说| 一区二区日韩欧美中文字幕| 男女边摸边吃奶| 看非洲黑人一级黄片| 青青草视频在线视频观看| 麻豆乱淫一区二区| 亚洲欧美成人精品一区二区| 制服诱惑二区| 精品一区二区三区av网在线观看 | netflix在线观看网站| 天天躁夜夜躁狠狠躁躁| 午夜福利在线免费观看网站| 波多野结衣一区麻豆| av国产精品久久久久影院| 如何舔出高潮| 叶爱在线成人免费视频播放| 久久久久久久久免费视频了| 久久久久久久大尺度免费视频| 国产欧美亚洲国产| 视频区图区小说| 9色porny在线观看| 伊人久久国产一区二区| 亚洲第一区二区三区不卡| 国产成人欧美在线观看 | 可以免费在线观看a视频的电影网站 | 狂野欧美激情性bbbbbb| 久久久精品免费免费高清| 日本av免费视频播放| 色视频在线一区二区三区| 飞空精品影院首页| 免费高清在线观看日韩| 亚洲欧洲日产国产| 九色亚洲精品在线播放| 高清视频免费观看一区二区| 两性夫妻黄色片| 人妻人人澡人人爽人人| 欧美97在线视频| 久久久久久久久久久久大奶| 校园人妻丝袜中文字幕| av线在线观看网站| 亚洲国产日韩一区二区| 桃花免费在线播放| 免费在线观看视频国产中文字幕亚洲 | 国产精品 国内视频| 欧美国产精品一级二级三级| 久热爱精品视频在线9| 国产又爽黄色视频| 大陆偷拍与自拍| 纵有疾风起免费观看全集完整版| 青草久久国产| 亚洲精品久久午夜乱码| 日日啪夜夜爽| 母亲3免费完整高清在线观看| www.精华液| 国产 精品1| 考比视频在线观看| 制服诱惑二区| 丁香六月欧美| 波多野结衣av一区二区av| 免费女性裸体啪啪无遮挡网站| 可以免费在线观看a视频的电影网站 | 欧美精品高潮呻吟av久久| 99国产综合亚洲精品| 色吧在线观看| 色婷婷久久久亚洲欧美| 宅男免费午夜| 亚洲熟女精品中文字幕| 一本久久精品| av视频免费观看在线观看| 久久精品亚洲av国产电影网| 久久 成人 亚洲| 超碰97精品在线观看| 亚洲第一av免费看| 亚洲av国产av综合av卡| 日日爽夜夜爽网站| 亚洲成人免费av在线播放| 国产成人精品久久久久久| 亚洲国产av新网站| 午夜福利影视在线免费观看| av在线老鸭窝| 母亲3免费完整高清在线观看| 精品一区二区免费观看| 97人妻天天添夜夜摸| 黄片小视频在线播放| 91aial.com中文字幕在线观看| 丁香六月欧美| 国产又爽黄色视频| 在线观看免费日韩欧美大片| netflix在线观看网站| 久久久国产欧美日韩av| 一级,二级,三级黄色视频| 黄片播放在线免费| 亚洲图色成人| 国产精品 国内视频| 青青草视频在线视频观看| 亚洲美女搞黄在线观看| 交换朋友夫妻互换小说| 亚洲精品av麻豆狂野| 在线亚洲精品国产二区图片欧美| 久久女婷五月综合色啪小说| 欧美中文综合在线视频| 天堂俺去俺来也www色官网| 国产一区二区激情短视频 | 老司机靠b影院| 精品人妻一区二区三区麻豆| 日韩一本色道免费dvd| 久久人人爽av亚洲精品天堂| 国产一级毛片在线| 亚洲av电影在线进入| 中文字幕色久视频| 亚洲国产精品成人久久小说| 亚洲精品一二三| 99香蕉大伊视频| av卡一久久| 晚上一个人看的免费电影| 国产精品av久久久久免费| 在线观看三级黄色| 久久久久国产精品人妻一区二区| www.av在线官网国产| 十分钟在线观看高清视频www| 黄色视频在线播放观看不卡| 欧美黑人欧美精品刺激| 色婷婷av一区二区三区视频| 成年人午夜在线观看视频| 欧美老熟妇乱子伦牲交| 日日撸夜夜添| 青春草国产在线视频| 性高湖久久久久久久久免费观看| 免费观看a级毛片全部| 国产一区二区 视频在线| 亚洲av成人不卡在线观看播放网 | 婷婷色av中文字幕| 你懂的网址亚洲精品在线观看| 日本猛色少妇xxxxx猛交久久| 国产一级毛片在线| 欧美xxⅹ黑人| 中国国产av一级| 成年av动漫网址| 午夜福利网站1000一区二区三区| 校园人妻丝袜中文字幕| 午夜影院在线不卡| 99精品久久久久人妻精品| 亚洲自偷自拍图片 自拍| 18禁国产床啪视频网站| 色婷婷av一区二区三区视频| 另类精品久久| 精品免费久久久久久久清纯 | 欧美日韩亚洲国产一区二区在线观看 | 热99国产精品久久久久久7| 十八禁高潮呻吟视频| 欧美xxⅹ黑人| av一本久久久久| 国产精品蜜桃在线观看| 国产精品三级大全| 熟妇人妻不卡中文字幕| 交换朋友夫妻互换小说| 91精品伊人久久大香线蕉| 欧美黑人精品巨大| 视频区图区小说| 国产 精品1| 女性生殖器流出的白浆| 欧美另类一区| 久久久久久久大尺度免费视频| 高清黄色对白视频在线免费看| 精品一区二区三区av网在线观看 | 国产成人午夜福利电影在线观看| 欧美中文综合在线视频| 欧美 日韩 精品 国产| 51午夜福利影视在线观看| 国产精品一区二区在线不卡| 精品人妻一区二区三区麻豆| 精品免费久久久久久久清纯 | 七月丁香在线播放| 韩国高清视频一区二区三区| 欧美日韩亚洲综合一区二区三区_| 人人妻人人澡人人看| 在现免费观看毛片| 成人午夜精彩视频在线观看| 久久精品久久久久久噜噜老黄| 日韩 亚洲 欧美在线| 看免费成人av毛片| 美女扒开内裤让男人捅视频| 高清在线视频一区二区三区| 国产黄频视频在线观看| 欧美 亚洲 国产 日韩一| 亚洲,欧美精品.| 国产成人91sexporn| 一级黄片播放器| 成人亚洲欧美一区二区av| 国产成人啪精品午夜网站| 久久久久久人人人人人| 国产女主播在线喷水免费视频网站| 国精品久久久久久国模美| 美女中出高潮动态图| 观看av在线不卡| 国产 精品1| 久久女婷五月综合色啪小说| 国产日韩一区二区三区精品不卡| 少妇精品久久久久久久| 99香蕉大伊视频| 水蜜桃什么品种好| 国产av码专区亚洲av| 国产日韩一区二区三区精品不卡| 国产成人系列免费观看| 亚洲欧洲国产日韩| 啦啦啦啦在线视频资源| 亚洲一码二码三码区别大吗| 丁香六月欧美| 亚洲成人国产一区在线观看 | 制服丝袜香蕉在线| 亚洲 欧美一区二区三区| 成年动漫av网址| 精品一区二区免费观看| 久久久久久人人人人人| 国产精品成人在线| 国产片特级美女逼逼视频| 久久久久久久久久久免费av| 欧美日韩视频高清一区二区三区二| 日韩视频在线欧美| 欧美日韩福利视频一区二区| 久久久久国产精品人妻一区二区| 99国产综合亚洲精品| 天美传媒精品一区二区| 成人午夜精彩视频在线观看| 各种免费的搞黄视频| 妹子高潮喷水视频| 国产成人精品无人区| 在线观看免费高清a一片| 精品人妻在线不人妻| 少妇人妻久久综合中文| 激情视频va一区二区三区| 精品免费久久久久久久清纯 | 亚洲国产看品久久| 黑人猛操日本美女一级片| 无遮挡黄片免费观看| 最黄视频免费看| 狠狠精品人妻久久久久久综合| 看非洲黑人一级黄片| 自线自在国产av| 日本vs欧美在线观看视频| 国产乱人偷精品视频| 日韩一卡2卡3卡4卡2021年| 乱人伦中国视频| 日本wwww免费看| 中文字幕制服av| 亚洲欧美成人综合另类久久久| 自线自在国产av| 色婷婷av一区二区三区视频| 国产xxxxx性猛交| www.精华液| 亚洲精华国产精华液的使用体验| xxxhd国产人妻xxx| 久久女婷五月综合色啪小说| 99久国产av精品国产电影| 人成视频在线观看免费观看| 少妇人妻精品综合一区二区| 久久性视频一级片| 成年av动漫网址| av片东京热男人的天堂| 国产精品久久久av美女十八| 国产熟女午夜一区二区三区| 超碰97精品在线观看| 男女边摸边吃奶| www日本在线高清视频| 90打野战视频偷拍视频| 中文字幕av电影在线播放| 宅男免费午夜| 国产人伦9x9x在线观看| 国产av一区二区精品久久| 七月丁香在线播放| 自线自在国产av| 一区在线观看完整版| 麻豆精品久久久久久蜜桃| av网站免费在线观看视频| 人妻人人澡人人爽人人| 国产av精品麻豆| 国产1区2区3区精品| 欧美人与善性xxx| 激情五月婷婷亚洲| 精品国产一区二区三区四区第35| 十分钟在线观看高清视频www| 天天躁夜夜躁狠狠躁躁| 69精品国产乱码久久久| 悠悠久久av| 久久精品亚洲av国产电影网| 韩国高清视频一区二区三区| 亚洲精品日韩在线中文字幕| 777久久人妻少妇嫩草av网站| 天堂中文最新版在线下载| 欧美乱码精品一区二区三区| 免费高清在线观看日韩| 国产免费一区二区三区四区乱码| 波野结衣二区三区在线| 毛片一级片免费看久久久久| 国产欧美亚洲国产| 亚洲精品av麻豆狂野| 深夜精品福利| 亚洲国产精品一区三区| 久久久欧美国产精品| 国产男人的电影天堂91| 免费黄频网站在线观看国产| 日本vs欧美在线观看视频| 蜜桃国产av成人99| 久久99一区二区三区| 国产免费视频播放在线视频| 亚洲在久久综合| 久久ye,这里只有精品| 一区二区日韩欧美中文字幕| 2018国产大陆天天弄谢| 免费黄频网站在线观看国产| 免费不卡黄色视频| 极品少妇高潮喷水抽搐| 精品一区二区免费观看| 99久久人妻综合| 日韩精品有码人妻一区| 美女视频免费永久观看网站| 制服诱惑二区| 亚洲精品国产区一区二| 中文欧美无线码| 欧美精品亚洲一区二区| 国产成人精品久久二区二区91 | 国产爽快片一区二区三区| 69精品国产乱码久久久| 久久国产亚洲av麻豆专区| 在线天堂中文资源库| xxxhd国产人妻xxx| 九草在线视频观看| 毛片一级片免费看久久久久| 国产男女超爽视频在线观看| 人妻人人澡人人爽人人| 黄频高清免费视频| 欧美国产精品va在线观看不卡| 欧美日韩综合久久久久久| 精品卡一卡二卡四卡免费| 女人久久www免费人成看片| 午夜免费鲁丝| 欧美老熟妇乱子伦牲交| 91老司机精品| 捣出白浆h1v1| 热99国产精品久久久久久7| 亚洲精品日韩在线中文字幕| 一区二区三区激情视频| 汤姆久久久久久久影院中文字幕| 最近的中文字幕免费完整| 久久天躁狠狠躁夜夜2o2o | 人人澡人人妻人| 亚洲成人国产一区在线观看 | 亚洲情色 制服丝袜| 考比视频在线观看| 在线观看国产h片| 黑人巨大精品欧美一区二区蜜桃| 视频区图区小说| 伊人久久大香线蕉亚洲五| 又大又黄又爽视频免费| 久久人人97超碰香蕉20202| 日本av免费视频播放| 亚洲精品中文字幕在线视频| 免费观看a级毛片全部| 999精品在线视频| 黄片无遮挡物在线观看| 99热全是精品| 久久99精品国语久久久| 丰满饥渴人妻一区二区三| 亚洲精品,欧美精品| 韩国av在线不卡| 国产黄色免费在线视频| 侵犯人妻中文字幕一二三四区| 超色免费av| 婷婷成人精品国产| 欧美精品一区二区免费开放| 在线观看www视频免费| 国产精品无大码| 亚洲精品成人av观看孕妇| 免费人妻精品一区二区三区视频| 亚洲精品在线美女| 别揉我奶头~嗯~啊~动态视频 | 成人三级做爰电影| 在线观看免费视频网站a站| 亚洲精品国产色婷婷电影| 一本一本久久a久久精品综合妖精| 人人澡人人妻人| 久久久国产精品麻豆| 午夜福利视频精品| 国产高清不卡午夜福利| 91精品伊人久久大香线蕉| 国产精品av久久久久免费| 免费日韩欧美在线观看| 母亲3免费完整高清在线观看| 一边摸一边做爽爽视频免费| 日本色播在线视频| 欧美日韩亚洲高清精品| 99热国产这里只有精品6| 亚洲精品一区蜜桃| 中文字幕色久视频| 久久婷婷青草| 国产精品 欧美亚洲| 尾随美女入室| 在线观看三级黄色| 日韩一区二区视频免费看| av一本久久久久| 国产不卡av网站在线观看| 天美传媒精品一区二区| 亚洲av电影在线进入| 国产精品熟女久久久久浪| 久久天躁狠狠躁夜夜2o2o | 精品卡一卡二卡四卡免费| 亚洲国产精品一区二区三区在线| 国产精品久久久久久久久免| 日韩av免费高清视频| 考比视频在线观看| 久久久久久久久免费视频了| 国产精品.久久久| 日韩免费高清中文字幕av| 国产一卡二卡三卡精品 | 亚洲三区欧美一区| 亚洲精品,欧美精品| 中文字幕制服av| av片东京热男人的天堂| 青春草视频在线免费观看| 欧美精品av麻豆av| 亚洲精品美女久久av网站| 色播在线永久视频| 天堂中文最新版在线下载| 美女国产高潮福利片在线看| 又大又爽又粗| 亚洲欧美清纯卡通| 久久精品亚洲av国产电影网| 午夜激情av网站| 亚洲精品,欧美精品| 国产xxxxx性猛交| 美女福利国产在线| 免费人妻精品一区二区三区视频| 国产又爽黄色视频| av卡一久久| 久久性视频一级片| 国产淫语在线视频| 亚洲欧美一区二区三区久久| 青草久久国产| 男人添女人高潮全过程视频| 久久99热这里只频精品6学生| 电影成人av| 丁香六月天网| 亚洲色图综合在线观看| 亚洲一区二区三区欧美精品| 国产爽快片一区二区三区| 青春草亚洲视频在线观看| 男女免费视频国产| 亚洲av成人不卡在线观看播放网 | 狂野欧美激情性xxxx| 精品酒店卫生间| 一级黄片播放器| 中文精品一卡2卡3卡4更新| 国产精品一区二区在线观看99| 日韩一区二区三区影片| 啦啦啦在线观看免费高清www| 90打野战视频偷拍视频| 性色av一级| 熟女少妇亚洲综合色aaa.| 最近中文字幕2019免费版| 亚洲成人一二三区av|