• <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)
    人妻少妇偷人精品九色| 听说在线观看完整版免费高清| 亚洲18禁久久av| 91在线精品国自产拍蜜月| 麻豆国产av国片精品| 亚洲美女搞黄在线观看 | 精品久久久久久久人妻蜜臀av| 亚洲经典国产精华液单| 午夜免费成人在线视频| 免费观看人在逋| 亚洲欧美激情综合另类| av专区在线播放| 亚洲人成网站在线播放欧美日韩| 国产精品一及| 亚洲精华国产精华精| 干丝袜人妻中文字幕| 窝窝影院91人妻| 永久网站在线| 亚洲精品色激情综合| 最后的刺客免费高清国语| 蜜桃久久精品国产亚洲av| 久久精品国产亚洲av涩爱 | 麻豆成人av在线观看| 午夜福利18| 国产一区二区三区在线臀色熟女| 乱系列少妇在线播放| 国产大屁股一区二区在线视频| 亚洲欧美激情综合另类| 午夜福利在线观看免费完整高清在 | 久久久久久伊人网av| 1000部很黄的大片| 国产精品无大码| 精品人妻视频免费看| 91麻豆精品激情在线观看国产| 欧美bdsm另类| 精品人妻1区二区| 自拍偷自拍亚洲精品老妇| 嫩草影视91久久| 亚洲成人久久性| 91久久精品电影网| ponron亚洲| 欧美又色又爽又黄视频| 18禁黄网站禁片免费观看直播| 免费电影在线观看免费观看| 狂野欧美激情性xxxx在线观看| 99精品久久久久人妻精品| 日韩大尺度精品在线看网址| 波多野结衣高清无吗| 欧美性猛交黑人性爽| 久99久视频精品免费| 精品久久久噜噜| 午夜福利在线观看免费完整高清在 | 麻豆一二三区av精品| 国产精品一区二区性色av| 天堂动漫精品| 亚洲人与动物交配视频| 天堂av国产一区二区熟女人妻| 国产高清视频在线播放一区| 久久久精品欧美日韩精品| 午夜激情欧美在线| 日韩国内少妇激情av| 亚洲美女黄片视频| 免费人成视频x8x8入口观看| 国产大屁股一区二区在线视频| 亚洲成人免费电影在线观看| 有码 亚洲区| 熟女人妻精品中文字幕| 午夜日韩欧美国产| 少妇人妻一区二区三区视频| 夜夜爽天天搞| 国产高清三级在线| 亚洲 国产 在线| 精品久久国产蜜桃| 国产精品伦人一区二区| 伦理电影大哥的女人| 热99在线观看视频| 给我免费播放毛片高清在线观看| 人妻丰满熟妇av一区二区三区| 夜夜看夜夜爽夜夜摸| 亚洲av免费在线观看| 亚洲四区av| 日韩大尺度精品在线看网址| 国产精品福利在线免费观看| 亚洲精品影视一区二区三区av| 直男gayav资源| 夜夜看夜夜爽夜夜摸| 欧美精品国产亚洲| 欧美色视频一区免费| 毛片一级片免费看久久久久 | 动漫黄色视频在线观看| 国产爱豆传媒在线观看| 日韩精品青青久久久久久| 欧美日韩精品成人综合77777| 国内少妇人妻偷人精品xxx网站| 亚洲 国产 在线| 久久精品久久久久久噜噜老黄 | 69av精品久久久久久| 久久九九热精品免费| 女人被狂操c到高潮| 国产在视频线在精品| 亚洲人成网站在线播放欧美日韩| 琪琪午夜伦伦电影理论片6080| 亚洲欧美激情综合另类| 看十八女毛片水多多多| 国产黄a三级三级三级人| 热99re8久久精品国产| 91在线观看av| 久久99热这里只有精品18| 久久久久九九精品影院| a级一级毛片免费在线观看| 波多野结衣高清无吗| 免费观看人在逋| 精品99又大又爽又粗少妇毛片 | 人妻少妇偷人精品九色| 人妻制服诱惑在线中文字幕| 五月伊人婷婷丁香| 色精品久久人妻99蜜桃| 亚洲av中文字字幕乱码综合| 3wmmmm亚洲av在线观看| 久久午夜亚洲精品久久| 色综合婷婷激情| 亚洲精品日韩av片在线观看| 国产精华一区二区三区| 欧美+日韩+精品| 国产不卡一卡二| 久久久久国内视频| 看片在线看免费视频| 午夜福利18| 免费搜索国产男女视频| 成人欧美大片| 国产不卡一卡二| 久久久久九九精品影院| 91狼人影院| 一a级毛片在线观看| 中国美女看黄片| 国内精品久久久久精免费| 亚洲 国产 在线| 亚洲av第一区精品v没综合| 国产久久久一区二区三区| 欧美三级亚洲精品| 国产精品女同一区二区软件 | 国产精品国产三级国产av玫瑰| 亚洲最大成人手机在线| 免费av观看视频| 欧美日韩精品成人综合77777| 日韩欧美精品免费久久| 国产伦精品一区二区三区四那| av在线观看视频网站免费| 狠狠狠狠99中文字幕| av黄色大香蕉| 精品人妻1区二区| 亚洲国产精品成人综合色| 老女人水多毛片| 97碰自拍视频| 久久久久国内视频| 熟妇人妻久久中文字幕3abv| 久久精品影院6| 亚洲成人久久爱视频| 一卡2卡三卡四卡精品乱码亚洲| 国产高清激情床上av| 精品久久久久久,| 亚洲av美国av| 亚洲av免费高清在线观看| 精品人妻一区二区三区麻豆 | 亚洲无线在线观看| 国产亚洲精品久久久com| 国产av一区在线观看免费| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品一区av在线观看| 久久人人爽人人爽人人片va| 在线观看午夜福利视频| 亚洲美女搞黄在线观看 | 精品福利观看| 搡老岳熟女国产| 免费大片18禁| 精品无人区乱码1区二区| 小蜜桃在线观看免费完整版高清| 99视频精品全部免费 在线| 此物有八面人人有两片| 亚洲av日韩精品久久久久久密| www.www免费av| 淫妇啪啪啪对白视频| 久久欧美精品欧美久久欧美| 俺也久久电影网| 日韩欧美在线乱码| 99热只有精品国产| 成人亚洲精品av一区二区| 日本黄色视频三级网站网址| 99国产精品一区二区蜜桃av| 小蜜桃在线观看免费完整版高清| 久久婷婷人人爽人人干人人爱| 久9热在线精品视频| 特级一级黄色大片| 97超视频在线观看视频| 99久国产av精品| 成人特级黄色片久久久久久久| 日韩强制内射视频| 99久久成人亚洲精品观看| 男人舔奶头视频| 国产精品久久久久久av不卡| 亚洲无线观看免费| 久久精品影院6| 如何舔出高潮| 1000部很黄的大片| 日韩高清综合在线| 人妻制服诱惑在线中文字幕| 日韩强制内射视频| 国产色爽女视频免费观看| 亚洲真实伦在线观看| 亚洲国产欧美人成| 成人二区视频| 久久久久久九九精品二区国产| 欧美黑人欧美精品刺激| 狂野欧美白嫩少妇大欣赏| 老熟妇乱子伦视频在线观看| 免费大片18禁| 小蜜桃在线观看免费完整版高清| 亚洲18禁久久av| 中出人妻视频一区二区| 日韩欧美在线乱码| 久久久久免费精品人妻一区二区| 午夜精品久久久久久毛片777| 真实男女啪啪啪动态图| 两个人的视频大全免费| 88av欧美| 男人的好看免费观看在线视频| 国产探花在线观看一区二区| 亚洲国产精品成人综合色| 亚洲第一电影网av| 成人三级黄色视频| 简卡轻食公司| 日本-黄色视频高清免费观看| 我要搜黄色片| 日日夜夜操网爽| 色噜噜av男人的天堂激情| 能在线免费观看的黄片| av.在线天堂| 偷拍熟女少妇极品色| 老女人水多毛片| 91av网一区二区| 精品人妻熟女av久视频| 国产精品亚洲一级av第二区| 欧美激情在线99| 赤兔流量卡办理| 99热网站在线观看| 性欧美人与动物交配| 12—13女人毛片做爰片一| 欧美在线一区亚洲| 内射极品少妇av片p| 亚洲av美国av| 精品人妻一区二区三区麻豆 | 亚洲 国产 在线| 日韩高清综合在线| h日本视频在线播放| 极品教师在线免费播放| 人妻丰满熟妇av一区二区三区| 亚洲精品成人久久久久久| 久久久国产成人免费| 网址你懂的国产日韩在线| 乱人视频在线观看| 中文字幕av在线有码专区| 香蕉av资源在线| 免费看日本二区| 国内少妇人妻偷人精品xxx网站| 国模一区二区三区四区视频| 精品福利观看| 丝袜美腿在线中文| 国产精品乱码一区二三区的特点| 丰满乱子伦码专区| 久久久久久久亚洲中文字幕| 日日撸夜夜添| 久久久久久久久久久丰满 | 别揉我奶头~嗯~啊~动态视频| 国产乱人伦免费视频| 18禁在线播放成人免费| 亚洲精华国产精华液的使用体验 | 亚洲成人中文字幕在线播放| 三级毛片av免费| 日本黄色视频三级网站网址| 亚洲va在线va天堂va国产| 欧美日韩瑟瑟在线播放| 少妇裸体淫交视频免费看高清| 日日夜夜操网爽| 国产91精品成人一区二区三区| 国产v大片淫在线免费观看| 嫩草影视91久久| 亚洲四区av| 波野结衣二区三区在线| 一个人观看的视频www高清免费观看| 韩国av一区二区三区四区| 哪里可以看免费的av片| 精品一区二区三区av网在线观看| 免费不卡的大黄色大毛片视频在线观看 | 欧美xxxx性猛交bbbb| 国语自产精品视频在线第100页| 俄罗斯特黄特色一大片| 日日摸夜夜添夜夜添av毛片 | 日韩国内少妇激情av| a级毛片a级免费在线| 亚洲在线观看片| 九九久久精品国产亚洲av麻豆| 久久人人精品亚洲av| 一本久久中文字幕| 日本成人三级电影网站| 毛片一级片免费看久久久久 | 久久久久久大精品| 日日摸夜夜添夜夜添小说| 亚洲av美国av| 精品久久久久久久久久免费视频| 美女高潮喷水抽搐中文字幕| 久久国产精品人妻蜜桃| 久久精品国产亚洲av香蕉五月| 久久久国产成人精品二区| 男人舔奶头视频| 成人性生交大片免费视频hd| 亚洲一区二区三区色噜噜| 午夜老司机福利剧场| 国产精品国产高清国产av| 制服丝袜大香蕉在线| 成人欧美大片| 日韩精品有码人妻一区| 岛国在线免费视频观看| 亚洲内射少妇av| 久久午夜亚洲精品久久| 看黄色毛片网站| 成人国产一区最新在线观看| 亚洲国产精品成人综合色| 精品久久久久久久久久久久久| 国产一区二区三区视频了| 天堂√8在线中文| 韩国av在线不卡| 人妻制服诱惑在线中文字幕| 国产精品免费一区二区三区在线| 乱系列少妇在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 美女免费视频网站| 午夜老司机福利剧场| 国产精品野战在线观看| 亚洲欧美日韩无卡精品| 一边摸一边抽搐一进一小说| 麻豆一二三区av精品| 午夜福利欧美成人| 99在线视频只有这里精品首页| 熟妇人妻久久中文字幕3abv| 国产成年人精品一区二区| 久久久久性生活片| 麻豆成人av在线观看| 亚洲性夜色夜夜综合| 久久国产精品人妻蜜桃| 亚洲欧美日韩无卡精品| 国产精品嫩草影院av在线观看 | 美女高潮喷水抽搐中文字幕| 天堂网av新在线| 97超视频在线观看视频| 人妻丰满熟妇av一区二区三区| 88av欧美| 窝窝影院91人妻| 国产成人福利小说| 午夜日韩欧美国产| 我要搜黄色片| 中文字幕av在线有码专区| 少妇被粗大猛烈的视频| 成人性生交大片免费视频hd| av在线蜜桃| 亚洲欧美日韩高清专用| 国产乱人视频| 成年免费大片在线观看| 国产真实伦视频高清在线观看 | 久久婷婷人人爽人人干人人爱| 九九爱精品视频在线观看| 国内揄拍国产精品人妻在线| 尤物成人国产欧美一区二区三区| 亚洲一级一片aⅴ在线观看| 国产高清三级在线| 哪里可以看免费的av片| 午夜日韩欧美国产| 国产伦人伦偷精品视频| 黄色配什么色好看| 亚洲五月天丁香| 三级男女做爰猛烈吃奶摸视频| 欧美xxxx黑人xx丫x性爽| 在线观看66精品国产| 97热精品久久久久久| 九九在线视频观看精品| 久久九九热精品免费| 国内精品美女久久久久久| 日韩亚洲欧美综合| 嫩草影视91久久| 国产日本99.免费观看| 97热精品久久久久久| 黄色一级大片看看| 性色avwww在线观看| av.在线天堂| av福利片在线观看| 日本 欧美在线| 人妻制服诱惑在线中文字幕| 午夜福利在线观看吧| 久久香蕉精品热| 国产亚洲精品久久久com| 国产极品精品免费视频能看的| 日本欧美国产在线视频| 亚洲精品亚洲一区二区| 成人鲁丝片一二三区免费| 亚洲精品成人久久久久久| 亚洲男人的天堂狠狠| 色综合亚洲欧美另类图片| 国产色爽女视频免费观看| 十八禁网站免费在线| 国产 一区 欧美 日韩| 婷婷色综合大香蕉| 能在线免费观看的黄片| 国产视频一区二区在线看| 狂野欧美激情性xxxx在线观看| 国产在视频线在精品| 神马国产精品三级电影在线观看| 高清毛片免费观看视频网站| 不卡视频在线观看欧美| 五月玫瑰六月丁香| 亚洲avbb在线观看| 国产男靠女视频免费网站| 麻豆国产av国片精品| 在线观看66精品国产| 制服丝袜大香蕉在线| 亚洲av中文字字幕乱码综合| 免费电影在线观看免费观看| 九九爱精品视频在线观看| 午夜影院日韩av| a在线观看视频网站| 国产精品乱码一区二三区的特点| 网址你懂的国产日韩在线| 长腿黑丝高跟| 乱码一卡2卡4卡精品| 一区福利在线观看| 中文在线观看免费www的网站| 亚州av有码| 窝窝影院91人妻| 欧美性猛交╳xxx乱大交人| 免费看光身美女| 观看美女的网站| 一个人看的www免费观看视频| 日本a在线网址| 亚洲av二区三区四区| 久久九九热精品免费| 久久精品人妻少妇| 女的被弄到高潮叫床怎么办 | 亚洲av电影不卡..在线观看| 91麻豆精品激情在线观看国产| 成人美女网站在线观看视频| 中文字幕熟女人妻在线| 国内精品一区二区在线观看| 亚洲va在线va天堂va国产| 91久久精品电影网| 男人狂女人下面高潮的视频| 色综合亚洲欧美另类图片| 亚洲天堂国产精品一区在线| 欧美zozozo另类| 18+在线观看网站| 九九热线精品视视频播放| 午夜免费成人在线视频| 国产精品乱码一区二三区的特点| 男女之事视频高清在线观看| 成年人黄色毛片网站| 免费黄网站久久成人精品| 成年人黄色毛片网站| 欧美潮喷喷水| 搡女人真爽免费视频火全软件 | 欧美绝顶高潮抽搐喷水| 欧美+亚洲+日韩+国产| 精品午夜福利在线看| 久99久视频精品免费| 国产白丝娇喘喷水9色精品| 久久精品综合一区二区三区| 在线播放国产精品三级| 91午夜精品亚洲一区二区三区 | 3wmmmm亚洲av在线观看| 一级a爱片免费观看的视频| 最新中文字幕久久久久| 亚洲最大成人手机在线| 在线天堂最新版资源| 一区二区三区四区激情视频 | 自拍偷自拍亚洲精品老妇| 男人的好看免费观看在线视频| 女同久久另类99精品国产91| 国产成人影院久久av| 亚洲成人久久性| 美女xxoo啪啪120秒动态图| 韩国av一区二区三区四区| 校园人妻丝袜中文字幕| 久久精品国产亚洲av涩爱 | 久久久久久久午夜电影| 亚洲午夜理论影院| 国产黄色小视频在线观看| 国产色爽女视频免费观看| 久久精品国产自在天天线| 免费看美女性在线毛片视频| 成熟少妇高潮喷水视频| 午夜影院日韩av| 一边摸一边抽搐一进一小说| 免费观看精品视频网站| 免费看a级黄色片| 欧美黑人巨大hd| av黄色大香蕉| 精品一区二区三区人妻视频| 亚洲av免费在线观看| 午夜精品久久久久久毛片777| 久久久久久久久久成人| 日日撸夜夜添| 成熟少妇高潮喷水视频| 中文字幕av成人在线电影| 欧美激情在线99| 日日夜夜操网爽| 亚洲精品日韩av片在线观看| 午夜精品久久久久久毛片777| 99视频精品全部免费 在线| a级毛片a级免费在线| 级片在线观看| 精品人妻视频免费看| 一本久久中文字幕| 国产精品98久久久久久宅男小说| 精品久久国产蜜桃| 精品一区二区三区av网在线观看| 日韩欧美一区二区三区在线观看| 亚洲欧美精品综合久久99| 久久久午夜欧美精品| av在线观看视频网站免费| 久久国内精品自在自线图片| 三级男女做爰猛烈吃奶摸视频| 狠狠狠狠99中文字幕| 欧美人与善性xxx| 偷拍熟女少妇极品色| 91在线观看av| 舔av片在线| 在线播放无遮挡| 国产精品自产拍在线观看55亚洲| 免费一级毛片在线播放高清视频| 亚洲熟妇中文字幕五十中出| 日本一本二区三区精品| 亚洲av.av天堂| 日本欧美国产在线视频| 国产麻豆成人av免费视频| 两个人的视频大全免费| 成人无遮挡网站| 日韩精品中文字幕看吧| 国产黄色小视频在线观看| 国产伦在线观看视频一区| 欧美性猛交黑人性爽| 天堂网av新在线| 国产私拍福利视频在线观看| 亚洲美女搞黄在线观看 | 亚洲精品一卡2卡三卡4卡5卡| 观看美女的网站| 午夜免费男女啪啪视频观看 | 日韩欧美免费精品| 日韩大尺度精品在线看网址| 亚洲av美国av| av在线老鸭窝| 久久久成人免费电影| 精品日产1卡2卡| 成人三级黄色视频| 老司机午夜福利在线观看视频| 午夜福利高清视频| 给我免费播放毛片高清在线观看| 精品不卡国产一区二区三区| 99国产精品一区二区蜜桃av| 久久久久国产精品人妻aⅴ院| 99国产精品一区二区蜜桃av| 一个人看的www免费观看视频| 中亚洲国语对白在线视频| 女同久久另类99精品国产91| 亚洲va在线va天堂va国产| 亚洲成人免费电影在线观看| 给我免费播放毛片高清在线观看| 免费观看在线日韩| 毛片女人毛片| 91久久精品电影网| 一进一出抽搐动态| 国产高清视频在线播放一区| 1024手机看黄色片| 亚洲一区二区三区色噜噜| 好男人在线观看高清免费视频| 日韩欧美 国产精品| 在线免费观看不下载黄p国产 | 精品一区二区三区av网在线观看| 国产伦人伦偷精品视频| av在线蜜桃| 99久国产av精品| 国产一区二区亚洲精品在线观看| .国产精品久久| 亚洲一区高清亚洲精品| 日本精品一区二区三区蜜桃| 国内精品一区二区在线观看| 熟女人妻精品中文字幕| 久久99热这里只有精品18| 中文资源天堂在线| 亚洲四区av| 日韩中字成人| a级一级毛片免费在线观看| 99riav亚洲国产免费| 国产一区二区三区av在线 | 观看免费一级毛片| 精品免费久久久久久久清纯| 国产高潮美女av| 国产精品久久电影中文字幕| av在线天堂中文字幕| 午夜福利在线观看免费完整高清在 | 身体一侧抽搐| 毛片一级片免费看久久久久 | h日本视频在线播放| www.www免费av| 午夜福利在线观看免费完整高清在 | 国产又黄又爽又无遮挡在线| 亚洲av成人av| 色综合婷婷激情| 国产黄色小视频在线观看| 免费无遮挡裸体视频|