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

    基于參數(shù)化CAD模型的船型阻力/耐波性一體化設(shè)計

    2011-03-06 03:06:10邱遼原姜治芳馮佰威劉祖源
    中國艦船研究 2011年1期
    關(guān)鍵詞:學(xué)科融合優(yōu)化

    邱遼原 謝 偉 姜治芳 馮佰威 劉祖源

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

    2武漢理工大學(xué) 交通學(xué)院,湖北 武漢 430063

    基于參數(shù)化CAD模型的船型阻力/耐波性一體化設(shè)計

    邱遼原1謝 偉1姜治芳1馮佰威2劉祖源2

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

    2武漢理工大學(xué) 交通學(xué)院,湖北 武漢 430063

    在艦船概念設(shè)計階段,往往需要快速生成阻力和耐波性能兼優(yōu)的船型。采用基于母型船的船型融合生成方法,實現(xiàn)了參數(shù)化船型自動生成。在此基礎(chǔ)上,采用iSight優(yōu)化平臺,將參數(shù)化船型生成技術(shù)與阻力、耐波性計算模型集成,運用多學(xué)科設(shè)計優(yōu)化技術(shù)實現(xiàn)了船型阻力/耐波性性能一體化設(shè)計。優(yōu)化方法采用多目標(biāo)遺傳算法以獲取Pareto前沿。以一艘46 000 DWT油船的型線優(yōu)化為算例對這個過程進行了具體說明,試驗結(jié)果表明總阻力降低了3%,驗證了這種方法的可行性。

    船舶;參數(shù)化船型;阻力;耐波性;多學(xué)科設(shè)計優(yōu)化

    1 引言

    隨著船舶航行環(huán)境的日益復(fù)雜,耐波性能已成為船型設(shè)計的一個重要性能指標(biāo)。由于耐波性能與船型密切相關(guān),因此船型設(shè)計不僅要考慮船舶的阻力性能,同時也要考慮耐波性能。船舶阻力/耐波性的一體化設(shè)計成了一個重要的研究問題。

    從優(yōu)化設(shè)計角度來看,船型阻力/耐波性一體化設(shè)計是一個多學(xué)科設(shè)計優(yōu)化問題。例如,意大利羅馬水池教授Peri D,利用多學(xué)科設(shè)計優(yōu)化方法,完成了阻力、耐波性兩學(xué)科三目標(biāo)的優(yōu)化問題[1-4];日本大阪大學(xué)的Yusuke Tahara[5]教授利用自主開發(fā)集成框架,完成了阻力和推進的多學(xué)科多目標(biāo)優(yōu)化。

    從實際應(yīng)用角度來看,由于目前CAD軟件在船舶設(shè)計各個階段中獲得了廣泛的應(yīng)用,將CAD軟件融入到多學(xué)科設(shè)計優(yōu)化框架中是必然的趨勢[6-7]。因此,在阻力/耐波性一體化設(shè)計中應(yīng)該集成CAD軟件,從而可帶來以下優(yōu)勢:

    1)CAD軟件提供了豐富的三維造型功能,能準(zhǔn)確地生成船體外形,可為阻力和耐波性能分析提供精確的外形模型;

    2)在以往的阻力、耐波性設(shè)計中,阻力學(xué)科和耐波性學(xué)科通常各自建立一個船體外形模型。而參數(shù)化CAD建??梢詾樽枇W(xué)科和耐波性學(xué)科提供統(tǒng)一的外形模型,避免了不同學(xué)科重復(fù)的外形建模,并消除了不同模型之間的差異;

    3)在船型設(shè)計中,有關(guān)船型的排水量和浮心縱向位置等船型參數(shù)非常重要,CAD軟件可以直接精確計算出各類靜力學(xué)要素;

    4)阻力/耐波性多學(xué)科設(shè)計優(yōu)化結(jié)果本身就是一個三維CAD模型,可被下游設(shè)計階段(詳細設(shè)計階段)直接利用,無需再根據(jù)優(yōu)化結(jié)果重新建立外形CAD模型。

    因此,將CAD軟件與阻力/耐波性多學(xué)科設(shè)計優(yōu)化結(jié)合起來具有實際應(yīng)用價值。本文針對CAD軟件與船型阻力/耐波性多學(xué)科設(shè)計優(yōu)化集成的問題,采用CAD二次開發(fā)技術(shù),提出了基于參數(shù)化CAD模型的阻力/耐波性一體化設(shè)計流程,以46 000 DWT油船阻力/耐波性一體化設(shè)計為算例,詳細闡述流程的實現(xiàn)過程。

    2 基于參數(shù)化CAD的阻力 /耐波性一體化設(shè)計流程

    基于參數(shù)化CAD模型的船型阻力/耐波性一體化設(shè)計流程如圖1所示。

    簡述如下:

    1)設(shè)計參數(shù)即為控制船型變化的優(yōu)化變量;

    2)船型參數(shù)化融合模塊讀入設(shè)計參數(shù),通過母型船的融合,生成新的船型;

    3)自動提取阻力計算及耐波性能計算所需的數(shù)據(jù)文件,同時計算排水量及浮心的縱向位置;

    4)阻力計算軟件及耐波性計算軟件分別讀取相應(yīng)數(shù)據(jù)文件,分別計算阻力及耐波性指標(biāo);

    5)選擇合適的優(yōu)化算法,進行船型阻力及耐波性能綜合優(yōu)化。

    自動評判獲得的船體形狀是否為阻力及耐波性能最優(yōu)的船型,如果不滿足要求,則再自動更改設(shè)計參數(shù),重復(fù)上述步驟的2)~5),直至找到阻力及耐波性能綜合兼優(yōu)的船型。

    上述整個流程可借助多學(xué)科設(shè)計優(yōu)化集成軟件 iSight來實現(xiàn)[8-9]。

    3 應(yīng)用算例

    以一艘46 000 DWT油船的型線優(yōu)化為算例,闡述圖1流程的具體實現(xiàn)過程。其阻力/耐波性一體化設(shè)計問題是:在滿足一定的約束條件下,優(yōu)化船型,改善阻力和耐波性能。

    該船主尺度排水量48 800 t,設(shè)計水線長176.2 m,船寬 32.4 m,吃水 10.5 m,型深 17.9 m,方形系數(shù) 0.794,浮心縱向相對位置 0.467。

    3.1 船型的參數(shù)化融合生成方法

    以型線生成CAD軟件Fastship為二次開發(fā)平臺,基于NURBS的船型描述基礎(chǔ)上,開發(fā)了船型融合自動生成方法。該方法是以現(xiàn)有的多條母型船為基礎(chǔ),通過融合系數(shù)(權(quán)重因子)的調(diào)節(jié),產(chǎn)生一系列光順的船型。而這一融合的過程則是直接修改母型船的NURBS控制頂點,再由融合生成的控制頂點產(chǎn)生船體曲面的網(wǎng)格,進而生成船體曲面。在融合過程中一定要保證融合系數(shù)(權(quán)重因子)的總和為1。其融合過程如下:

    從上面的融合過程可以看到,因融合系數(shù)的和為1,因此無論怎樣調(diào)節(jié)Ci的值,融合后生成的船型則總是在正好容納所有母型船的最小空間內(nèi),如圖2所示。

    另外,如果在優(yōu)化過程中主尺度也作為變量,那么還須要將原主尺度下的融合生成的控制頂點坐標(biāo)按照比例縮放到當(dāng)前的主尺度下的頂點坐標(biāo)。現(xiàn)以船長X方向的縮放為例,得到新船控制頂點X坐標(biāo)的過程可表示為:

    式中,Xnew表示新船的控制頂點X坐標(biāo);Xblending表示融合船的控制頂點X坐標(biāo);Lbasisship表示母型船的船長;Lvariable表示通過優(yōu)化器得到的船長。船寬Y方向和型深Z方向的縮放與此類似。

    3.2 阻力及耐波性能計算

    3.2.1 阻力計算

    船舶總阻力Rt用粘壓阻力Rνp、摩擦阻力Rf和興波阻力Rw三者之和來表達,即Rt=Rνp+Rf+Rw。

    1)摩擦阻力Rf及粘壓阻力Rνp

    采用 Holtrop 方法估算,即 Rf+Rνp=Rf(1+k)。

    式中摩擦阻力系數(shù)采用1957年ITTC公式計算,形狀因子k采用下式估算:

    c2(B /LR)0.92497(0.95-Cp)-0.521448(1-Cp+0.022 5 lcb)0.6906}式中,Cp為棱形系數(shù);lcb為浮心縱向坐標(biāo)占船長的百分數(shù);c1與尾部形狀有關(guān);c2與吃水船長比T/L有關(guān);LR為去流段長度;LR可用下式估算:

    2)粗糙度補貼ΔCf

    根據(jù)1975年第14屆國際船模試驗池會議公式計算粗糙度補貼系數(shù)計算:式中,ks為粗糙度表觀高度,本軟件取150×10-6m。

    3)興波阻力Rw

    采用Shipflow軟件計算,即Rw=CwρV2S。

    Cw的計算通過Shipflow軟件,采用邊界元法、線性自由面邊界條件以及船體表面用高階NURBS表達,源的分布為NURBS表達的高階源分布。Shipflow軟件計算興波阻力所需要的船型曲面上的數(shù)據(jù)均從圖2所示的CAD模型中直接獲得。

    3.2.2 耐波性能計算

    船舶在波浪上的運動計算程序是建立在線性理論基礎(chǔ)上,并假設(shè)流體為均勻、不可壓縮和無粘性的理想流體。在此基礎(chǔ)上建立船舶線性運動微分方程:

    式中,Mij,mij,Nij和 Cij分別為船舶質(zhì)量慣性力系數(shù)、附連質(zhì)量系數(shù)、阻尼系數(shù)和線性恢復(fù)力系數(shù)。

    方程中的水動力系數(shù)采用經(jīng)典的切片理論求解,進而求解運動微分方程得到船舶在波浪上的運動頻率響應(yīng)函數(shù)。耐波性計算所需要的船型各剖面的數(shù)據(jù)均從圖2所示的CAD模型中獲取。圖3為耐波性計算的程序流程。

    3.3 優(yōu)化問題的數(shù)學(xué)表達

    船型阻力及耐波性能一體化設(shè)計問題可用下面的方式來表達。

    目標(biāo)

    1)設(shè)計航速20 kn時的單位排水量總阻力最小,即Min;

    2)固定海況下的垂蕩峰值最小,即min(max-heave),其中maxheave=,ξmax為垂蕩峰值,h為波高;

    3)固定海況下的縱搖峰值最小,即min(maxpitch),其中maxpitch=,φmax為縱搖峰值,k為波數(shù)。

    變量

    船長L、船寬 B、吃水 T、型深 D、融合系數(shù) Ci。約束條件:

    1)浮力平衡約束;

    2)穩(wěn)性約束;

    3) ∑Ci=1 ,且0≤Ci≤1。

    本文利用iSight的過程集成及優(yōu)化策略定制功能,完成了上述各模塊的集成及優(yōu)化問題的表述。其部分界面見圖4。

    3.4 優(yōu)化結(jié)果及分析

    上述的優(yōu)化問題是一個多目標(biāo)優(yōu)化問題。對于這類問題的求解,通常不存在絕對的最優(yōu)解,實為Pareto解集。求解多目標(biāo)優(yōu)化問題的實質(zhì)就是確定Pareto解集[10]。本文采用Shinya提出的多目標(biāo)遺傳算法求解該優(yōu)化問題。優(yōu)化算法中的參數(shù)設(shè)置為:初始種群 20,代數(shù) 10,交叉概率為 1.0,變異概率為0.01。優(yōu)化后得到該問題的Pareto解集,該解集中共有44個優(yōu)化方案,見圖5及表1。

    表1 優(yōu)化后獲得的部分方案Tab.1 Some of solutions obtained by the optimization process

    如果從表中選擇優(yōu)化方案13,那么耐波性能及阻力性能分別提高了 1.3%、1.9%及 12.3%。

    為了驗證該方法的可行性,對優(yōu)化前的母型船和優(yōu)化方案13開展了模型阻力試驗。試驗結(jié)果表明,在設(shè)計航速20 kn附近時,優(yōu)化方案13總阻力較母型船降低約3%,表明該方法是可行的,可應(yīng)用于艦船方案論證和設(shè)計階段方案優(yōu)選,具有一定的工程實用價值。

    4 結(jié)束語

    本文提出了一種采用基于母型船的船型融合生成方法,實現(xiàn)了船型參數(shù)化自動生成。在此基礎(chǔ)上應(yīng)用CAD二次開發(fā)技術(shù),將CAD與阻力計算模型、耐波性計算模型和多學(xué)科設(shè)計優(yōu)化算法集成起來,實現(xiàn)了船型阻力/耐波性一體化設(shè)計。以一艘46 000 DWT油船的型線優(yōu)化為算例,詳細闡述流程的具體實現(xiàn)過程,優(yōu)化結(jié)果和模型試驗驗證了這種方法的可行性,表明該方法具有一定的實用價值。

    在后續(xù)的研究中,將進一步擴展船型母型庫,確定合適的航行性能優(yōu)化目標(biāo),集成更多的船舶航行性能,實現(xiàn)船舶航行性能多學(xué)科設(shè)計優(yōu)化。

    [1]PERI D,ROSSETTI M,CAMPANA E F.Design optimization of ship hulls via CFD techniques [J].Journal of Ship Research,2001,45(2):141-149.

    [2]PERI D,CAMPANA E F.Multidisciplinary design optimzation of a naval combatant [J].Journal of Ship Research,2003,47(1):1-12.

    [3]PERI D,CAMPANA E F.High-fidelity models and multiobjective global optimization algorithms in simulationbased design [J].Journal of Ship Research,2005,49(3):159-175.

    [4]PERI D,CAMPANA E F,DATTOLA R.Multidisciplinary design optimization of a naval frigate [C]//10th AIAA/ISSMO Symposium on Multidisciplinary Analysis and Optimization.Albany,NY,2004:487-517.

    [5]TAHARA Y,TOHYAMA S.CFD-based multi-objective optimization method for ship design[J].International Journal for Numerical Methods in Fluids, 2006,52:499-527.

    [6]HARRIES S,ABT C.Formal hydrodynamic optimization of a fast monohull on the basis of parametric hull design[C]//5th International Conference on Fast Sea Transportation,Seattle,WA,1999:368-398.

    [7]HARRIES S,VALDENAZZI F,ABT C.Investigation on optimization strategies for the hydrodynamic design of fast ferries [C]//6th International Conference on Fast Sea Transportation.Southhampton,UK,2001:255-263.

    [8]馮佰威,劉祖源.多學(xué)科設(shè)計優(yōu)化技術(shù)在船舶設(shè)計中的應(yīng)用[J].中國造船,2009,50(4):15-18.

    [9]馮佰威,劉祖源.船舶CAD&CFD集成優(yōu)化接口開發(fā)及應(yīng)用研究[J].船舶工程,2009,31(增刊 1):32-35.

    [10] ZALEK S F,PARSONS M G,BECK R F.Naval hull form multicriterion hydrodynamic optimization for the conceptual design phase [J].Journal of Ship Reseach,2009,53(4):199-213.

    Integrated Ship Resistence/Seakeeping Optimization Based on Parametric CAD Model

    Qiu Liao-yuan1 Xie wei1 Jiang Zhi-fang1 Feng Bai-wei2 Liu Zu-yuan2
    1 China Ship Development and Design Center, Wuhan 430064,China
    2 School of Transportation,Wuhan University of Technology,Wuhan 430063,China

    The hull forms of ship both with optimal resistance and seakeeping performances are required to create quickly in the conceptual design phase.Parametric ship hulls were achieved automatically with a blending method according to the baseline ship design.Based on that, the resistance and seakeeping performance integrated optimization was implemented with multidisciplinary design optimization method,combining parametric ship hull blending method and the resistance and seakeeping performance procedures through iSight optimization software.A multi-criterion, population-based evolutionary algorithm optimization process was used to provide a more thorough representation of the numerically approximated Pareto front.The resistance and seakeeping integrated optimization of 46 000 DWT oil tanker was presented as an example to illustrate the process.Model test results show that total resistance reduction is up to 3%.The feasibility of this integrated optimization method was also validated.

    ship engineering; parametric CAD model; resistance; seakeeping; multidisciplinary design optimization

    U662.3

    A

    1673-3185(2011)01-18-04

    10.3969/j.issn.1673-3185.2011.01.004

    2010-07-15

    國家自然科學(xué)基金資助項目(10772139);國家863高技術(shù)研究發(fā)展計劃(2006AA04Z124);總裝“十一五”水動力預(yù)研基金

    邱遼原(1977-),男,博士。研究方向:艦船總體設(shè)計、船舶多學(xué)科綜合優(yōu)化技術(shù)。E-mail:Kiuly@tom.com

    猜你喜歡
    學(xué)科融合優(yōu)化
    【學(xué)科新書導(dǎo)覽】
    超限高層建筑結(jié)構(gòu)設(shè)計與優(yōu)化思考
    村企黨建聯(lián)建融合共贏
    融合菜
    民用建筑防煙排煙設(shè)計優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    土木工程學(xué)科簡介
    一道優(yōu)化題的幾何解法
    從創(chuàng)新出發(fā),與高考數(shù)列相遇、融合
    《融合》
    免费看av在线观看网站| 中文字幕免费在线视频6| 69人妻影院| 精品免费久久久久久久清纯| 午夜老司机福利剧场| 一本一本综合久久| 婷婷丁香在线五月| 国产精品电影一区二区三区| 日本熟妇午夜| 成人特级黄色片久久久久久久| 免费av观看视频| 国产在线男女| 精品不卡国产一区二区三区| 国产男靠女视频免费网站| 亚洲乱码一区二区免费版| 亚洲国产精品sss在线观看| 高清日韩中文字幕在线| 少妇裸体淫交视频免费看高清| 亚洲精华国产精华精| 免费高清视频大片| 少妇的逼好多水| 国国产精品蜜臀av免费| 欧美zozozo另类| 欧美色视频一区免费| 亚洲精品456在线播放app | 日韩欧美国产一区二区入口| 国产精品一及| 校园人妻丝袜中文字幕| 两个人视频免费观看高清| 舔av片在线| 免费不卡的大黄色大毛片视频在线观看 | 一级黄片播放器| 中文字幕人妻熟人妻熟丝袜美| 国产高清三级在线| 久久久久久久久久成人| 成人毛片a级毛片在线播放| www日本黄色视频网| 国产精品一区二区三区四区久久| 婷婷精品国产亚洲av在线| 国产爱豆传媒在线观看| 国内揄拍国产精品人妻在线| 国产精品国产三级国产av玫瑰| 国产精品爽爽va在线观看网站| 欧美在线一区亚洲| 免费看日本二区| 国产色婷婷99| 成人特级黄色片久久久久久久| 欧美日韩精品成人综合77777| 日本五十路高清| 男人舔女人下体高潮全视频| 老熟妇仑乱视频hdxx| 3wmmmm亚洲av在线观看| 我要看日韩黄色一级片| bbb黄色大片| 亚洲美女搞黄在线观看 | 国产视频内射| 看黄色毛片网站| 婷婷精品国产亚洲av| 男人舔奶头视频| 亚洲精品在线观看二区| 五月玫瑰六月丁香| 一本久久中文字幕| 亚洲美女搞黄在线观看 | 变态另类成人亚洲欧美熟女| 精品人妻视频免费看| 日韩欧美三级三区| 亚洲av熟女| 精品乱码久久久久久99久播| 国产极品精品免费视频能看的| 男女啪啪激烈高潮av片| 少妇丰满av| 欧美性感艳星| 国产精品一区www在线观看 | 99久久精品一区二区三区| 亚洲色图av天堂| 此物有八面人人有两片| 亚洲图色成人| 日日摸夜夜添夜夜添av毛片 | netflix在线观看网站| av福利片在线观看| 国产老妇女一区| 欧美+亚洲+日韩+国产| 十八禁网站免费在线| 欧美在线一区亚洲| 校园春色视频在线观看| 少妇高潮的动态图| 国产av在哪里看| 欧美日韩综合久久久久久 | 亚洲第一电影网av| 国产伦一二天堂av在线观看| 免费av不卡在线播放| 色在线成人网| 男女之事视频高清在线观看| 干丝袜人妻中文字幕| 国产精品野战在线观看| 亚洲专区中文字幕在线| 亚洲国产欧美人成| 美女 人体艺术 gogo| 淫妇啪啪啪对白视频| 精品久久久噜噜| 简卡轻食公司| 日本爱情动作片www.在线观看 | 久久午夜福利片| 男女之事视频高清在线观看| 免费不卡的大黄色大毛片视频在线观看 | 一级黄片播放器| 美女 人体艺术 gogo| 久久精品综合一区二区三区| 97碰自拍视频| 黄色一级大片看看| 精品久久久久久久末码| 亚洲最大成人手机在线| 毛片一级片免费看久久久久 | 舔av片在线| 天美传媒精品一区二区| 乱人视频在线观看| 午夜福利视频1000在线观看| 神马国产精品三级电影在线观看| АⅤ资源中文在线天堂| 毛片女人毛片| 三级国产精品欧美在线观看| 一级av片app| 长腿黑丝高跟| 国产一区二区三区视频了| 狂野欧美白嫩少妇大欣赏| 最近视频中文字幕2019在线8| 精品久久久久久久久亚洲 | 永久网站在线| 色尼玛亚洲综合影院| 日本欧美国产在线视频| 黄片wwwwww| 婷婷丁香在线五月| 国产精品国产三级国产av玫瑰| 可以在线观看毛片的网站| 一进一出抽搐动态| 久久99热这里只有精品18| 免费观看在线日韩| 国产黄片美女视频| 日韩,欧美,国产一区二区三区 | 亚洲熟妇熟女久久| 一本一本综合久久| 国产 一区 欧美 日韩| 一本一本综合久久| 观看免费一级毛片| 成人综合一区亚洲| 国产精品久久久久久久电影| 国产极品精品免费视频能看的| 亚洲精品亚洲一区二区| 香蕉av资源在线| 床上黄色一级片| 亚洲18禁久久av| 18禁在线播放成人免费| 麻豆av噜噜一区二区三区| 小蜜桃在线观看免费完整版高清| 精品一区二区三区视频在线观看免费| 欧美一区二区精品小视频在线| 热99在线观看视频| 夜夜爽天天搞| 日本三级黄在线观看| 国产亚洲精品av在线| 最近最新中文字幕大全电影3| 99久久精品一区二区三区| 午夜爱爱视频在线播放| 长腿黑丝高跟| 亚洲美女黄片视频| 在线a可以看的网站| 中国美白少妇内射xxxbb| 免费在线观看日本一区| 亚洲精品456在线播放app | 成人毛片a级毛片在线播放| 天堂影院成人在线观看| 好男人在线观看高清免费视频| 亚洲黑人精品在线| 午夜福利欧美成人| 丝袜美腿在线中文| 最近中文字幕高清免费大全6 | 精品人妻1区二区| 老熟妇乱子伦视频在线观看| 国产精品国产高清国产av| 亚洲成人精品中文字幕电影| 欧美成人性av电影在线观看| 男人狂女人下面高潮的视频| 老司机福利观看| a级一级毛片免费在线观看| 在线天堂最新版资源| 国产精品人妻久久久久久| 少妇熟女aⅴ在线视频| 一个人看视频在线观看www免费| 国产又黄又爽又无遮挡在线| 亚洲专区国产一区二区| 少妇高潮的动态图| 男人舔女人下体高潮全视频| 久久亚洲精品不卡| 欧美性感艳星| 亚洲av熟女| 久久午夜福利片| 九九久久精品国产亚洲av麻豆| 久久精品人妻少妇| 国产精品一区二区三区四区久久| 国产高清激情床上av| or卡值多少钱| 亚洲av不卡在线观看| 亚洲三级黄色毛片| 中文资源天堂在线| 久久精品人妻少妇| 在线观看午夜福利视频| 亚洲av成人精品一区久久| 亚洲专区中文字幕在线| av在线天堂中文字幕| 黄色丝袜av网址大全| 国产午夜福利久久久久久| 久久99热这里只有精品18| 欧美日韩乱码在线| 一级黄色大片毛片| 色综合亚洲欧美另类图片| 欧美性感艳星| 亚洲最大成人手机在线| 久久久久久久久久黄片| 最近在线观看免费完整版| 人妻丰满熟妇av一区二区三区| 免费av不卡在线播放| 国产伦精品一区二区三区视频9| 免费高清视频大片| aaaaa片日本免费| 久久久精品大字幕| 在线天堂最新版资源| 国产女主播在线喷水免费视频网站 | 天美传媒精品一区二区| 男人的好看免费观看在线视频| 99久久成人亚洲精品观看| 色综合色国产| 久久精品国产清高在天天线| 亚洲最大成人av| 男人和女人高潮做爰伦理| 一本久久中文字幕| 男人的好看免费观看在线视频| 全区人妻精品视频| 天堂√8在线中文| 久久久久久久午夜电影| 国产一区二区激情短视频| 午夜爱爱视频在线播放| av黄色大香蕉| 亚洲人成网站在线播| 日韩欧美在线乱码| 亚洲精品色激情综合| 啪啪无遮挡十八禁网站| 可以在线观看的亚洲视频| 久久精品国产亚洲网站| 国产亚洲欧美98| 亚洲一区高清亚洲精品| 免费观看精品视频网站| 婷婷精品国产亚洲av在线| 99久久久亚洲精品蜜臀av| 天堂动漫精品| 精品久久国产蜜桃| 亚洲熟妇中文字幕五十中出| 欧美日韩乱码在线| 黄色视频,在线免费观看| 内射极品少妇av片p| 午夜精品在线福利| 国产免费av片在线观看野外av| 精品国内亚洲2022精品成人| 亚洲四区av| 在线观看午夜福利视频| 婷婷六月久久综合丁香| 欧美激情久久久久久爽电影| 国产精品一区二区免费欧美| 中文在线观看免费www的网站| 少妇的逼水好多| 女人被狂操c到高潮| 国产精品一区www在线观看 | 日日撸夜夜添| 免费在线观看影片大全网站| 久久人人精品亚洲av| 国产三级中文精品| 亚洲精品色激情综合| 亚洲aⅴ乱码一区二区在线播放| 精品无人区乱码1区二区| 看十八女毛片水多多多| 国产精品乱码一区二三区的特点| 91av网一区二区| av在线天堂中文字幕| 一级黄片播放器| 日日夜夜操网爽| 精品一区二区三区视频在线观看免费| 在线观看av片永久免费下载| 国产欧美日韩精品亚洲av| 国产一区二区亚洲精品在线观看| 成人无遮挡网站| 国产色婷婷99| av在线蜜桃| 国产精品,欧美在线| 最近视频中文字幕2019在线8| 狂野欧美白嫩少妇大欣赏| av专区在线播放| 久久欧美精品欧美久久欧美| 国产精品人妻久久久影院| 久久久久九九精品影院| 国产久久久一区二区三区| 在线播放无遮挡| 久久国产乱子免费精品| 在线免费观看的www视频| 精品一区二区三区视频在线| 麻豆一二三区av精品| 国内少妇人妻偷人精品xxx网站| 草草在线视频免费看| 亚洲av熟女| 中文字幕免费在线视频6| 亚洲中文字幕日韩| 久久人人爽人人爽人人片va| 国产精品一区www在线观看 | 亚洲乱码一区二区免费版| 色精品久久人妻99蜜桃| 欧美极品一区二区三区四区| 欧美最新免费一区二区三区| 精品人妻1区二区| 日本免费a在线| 亚洲av免费高清在线观看| 免费观看精品视频网站| 欧美又色又爽又黄视频| 国产国拍精品亚洲av在线观看| 中国美女看黄片| 欧美一级a爱片免费观看看| 国产视频一区二区在线看| 亚洲精品亚洲一区二区| 两个人视频免费观看高清| 男插女下体视频免费在线播放| 成人av在线播放网站| 欧美人与善性xxx| 亚洲成人中文字幕在线播放| 最后的刺客免费高清国语| 欧美国产日韩亚洲一区| 特级一级黄色大片| 久久亚洲精品不卡| 夜夜爽天天搞| 亚洲欧美精品综合久久99| 欧美人与善性xxx| 国产不卡一卡二| 欧美在线一区亚洲| 成年女人毛片免费观看观看9| 网址你懂的国产日韩在线| 国内揄拍国产精品人妻在线| 两个人视频免费观看高清| 国产aⅴ精品一区二区三区波| 少妇高潮的动态图| 深夜a级毛片| 欧美黑人巨大hd| 最近最新中文字幕大全电影3| 国产一级毛片七仙女欲春2| 成人国产麻豆网| 日韩 亚洲 欧美在线| 亚洲三级黄色毛片| 嫩草影视91久久| 最近视频中文字幕2019在线8| 成熟少妇高潮喷水视频| 欧美人与善性xxx| 男人舔女人下体高潮全视频| 女人十人毛片免费观看3o分钟| 精品午夜福利视频在线观看一区| 淫妇啪啪啪对白视频| 欧美日本亚洲视频在线播放| 美女xxoo啪啪120秒动态图| 国产麻豆成人av免费视频| 国产精品乱码一区二三区的特点| 精品国产三级普通话版| 国产精品久久视频播放| 国产久久久一区二区三区| 国产精品亚洲一级av第二区| 热99在线观看视频| 中文字幕熟女人妻在线| 桃色一区二区三区在线观看| 一级黄色大片毛片| 精品午夜福利在线看| 欧美3d第一页| 日本黄大片高清| 亚洲狠狠婷婷综合久久图片| 内地一区二区视频在线| av.在线天堂| 97超视频在线观看视频| 两人在一起打扑克的视频| 日韩,欧美,国产一区二区三区 | or卡值多少钱| 亚洲精品456在线播放app | 淫秽高清视频在线观看| 国产 一区 欧美 日韩| 国产精品久久久久久亚洲av鲁大| 露出奶头的视频| 国产大屁股一区二区在线视频| 精品久久久久久,| 毛片女人毛片| 熟妇人妻久久中文字幕3abv| 日本精品一区二区三区蜜桃| 午夜福利成人在线免费观看| 久久久久久九九精品二区国产| 午夜福利在线观看吧| 国产在线男女| 99久久中文字幕三级久久日本| 国产成人aa在线观看| 国语自产精品视频在线第100页| 国产精品爽爽va在线观看网站| 国内精品久久久久精免费| 久久精品国产清高在天天线| 午夜福利成人在线免费观看| 亚洲成人精品中文字幕电影| 国产精品嫩草影院av在线观看 | 一进一出好大好爽视频| 老司机午夜福利在线观看视频| 男女啪啪激烈高潮av片| 国产亚洲精品久久久com| 黄色配什么色好看| 一个人看的www免费观看视频| 狂野欧美激情性xxxx在线观看| 中文字幕久久专区| 国产精品久久久久久精品电影| 国产美女午夜福利| 成人高潮视频无遮挡免费网站| 亚洲精品色激情综合| 一级av片app| 国产高清不卡午夜福利| 欧美日韩乱码在线| 真人做人爱边吃奶动态| 免费电影在线观看免费观看| 精品久久久久久成人av| 免费在线观看影片大全网站| 亚洲精华国产精华精| 国产伦一二天堂av在线观看| 女生性感内裤真人,穿戴方法视频| 久久人妻av系列| 久久亚洲精品不卡| 色视频www国产| 日韩欧美国产一区二区入口| 国产精品嫩草影院av在线观看 | eeuss影院久久| 欧美xxxx性猛交bbbb| 亚洲国产高清在线一区二区三| 久久天躁狠狠躁夜夜2o2o| 校园春色视频在线观看| 日韩中字成人| 国产成人aa在线观看| av在线蜜桃| 亚洲va在线va天堂va国产| 国产亚洲精品av在线| 亚洲性夜色夜夜综合| 日日摸夜夜添夜夜添小说| 亚洲一级一片aⅴ在线观看| 老司机福利观看| 人人妻人人澡欧美一区二区| 99热6这里只有精品| 成人毛片a级毛片在线播放| 天天躁日日操中文字幕| 欧美精品啪啪一区二区三区| 免费看美女性在线毛片视频| 亚洲午夜理论影院| 麻豆国产97在线/欧美| 内地一区二区视频在线| av国产免费在线观看| 女人被狂操c到高潮| 欧美日本亚洲视频在线播放| 国产69精品久久久久777片| av天堂在线播放| 国产淫片久久久久久久久| 亚洲在线自拍视频| 久久久久久久久大av| 欧美色欧美亚洲另类二区| 人妻夜夜爽99麻豆av| 给我免费播放毛片高清在线观看| 日本-黄色视频高清免费观看| 欧美成人一区二区免费高清观看| 欧美日韩国产亚洲二区| 久久中文看片网| 久久国产精品人妻蜜桃| 乱码一卡2卡4卡精品| 三级国产精品欧美在线观看| 精品欧美国产一区二区三| 国产精品1区2区在线观看.| 日日干狠狠操夜夜爽| 国产伦在线观看视频一区| 最近视频中文字幕2019在线8| 大又大粗又爽又黄少妇毛片口| 国产私拍福利视频在线观看| 国产精品久久电影中文字幕| 国产精品无大码| 午夜精品久久久久久毛片777| 免费观看在线日韩| 国模一区二区三区四区视频| 少妇熟女aⅴ在线视频| 亚洲精品粉嫩美女一区| 尤物成人国产欧美一区二区三区| 久久精品国产自在天天线| 国产单亲对白刺激| 此物有八面人人有两片| 国产精品一区www在线观看 | videossex国产| 日本熟妇午夜| 亚洲av一区综合| 中文资源天堂在线| 成年女人永久免费观看视频| 99热精品在线国产| 亚洲无线在线观看| 日韩精品中文字幕看吧| 男女视频在线观看网站免费| 我的老师免费观看完整版| 欧美3d第一页| 99热只有精品国产| 日韩在线高清观看一区二区三区 | 日本一本二区三区精品| 深爱激情五月婷婷| 欧美区成人在线视频| 中文字幕精品亚洲无线码一区| www日本黄色视频网| 久久国产精品人妻蜜桃| 欧美一区二区亚洲| 一级黄色大片毛片| 精品无人区乱码1区二区| 草草在线视频免费看| 欧美潮喷喷水| 嫩草影院新地址| 亚洲中文字幕日韩| 久久久久性生活片| 香蕉av资源在线| 真人做人爱边吃奶动态| 99国产精品一区二区蜜桃av| 91久久精品电影网| 欧美极品一区二区三区四区| 国产成人福利小说| 亚洲国产精品合色在线| 欧美色视频一区免费| 一级a爱片免费观看的视频| 日本一本二区三区精品| 国产伦人伦偷精品视频| 国产高清不卡午夜福利| xxxwww97欧美| 丰满人妻一区二区三区视频av| 精品国内亚洲2022精品成人| 日本爱情动作片www.在线观看 | 十八禁网站免费在线| 国产伦一二天堂av在线观看| 国产亚洲精品久久久com| 少妇人妻一区二区三区视频| 国内精品久久久久精免费| 在线观看美女被高潮喷水网站| 欧美区成人在线视频| 亚洲成a人片在线一区二区| 免费电影在线观看免费观看| 伦精品一区二区三区| 午夜免费激情av| av在线观看视频网站免费| 尤物成人国产欧美一区二区三区| 亚洲性夜色夜夜综合| 成人特级黄色片久久久久久久| 国内精品久久久久久久电影| 色吧在线观看| 91狼人影院| 久久天躁狠狠躁夜夜2o2o| 欧美三级亚洲精品| 麻豆成人午夜福利视频| 五月伊人婷婷丁香| 免费观看精品视频网站| 国内毛片毛片毛片毛片毛片| 欧美精品啪啪一区二区三区| 一边摸一边抽搐一进一小说| 最好的美女福利视频网| 国产又黄又爽又无遮挡在线| 免费看美女性在线毛片视频| 亚洲一区二区三区色噜噜| 女的被弄到高潮叫床怎么办 | 成人高潮视频无遮挡免费网站| 99精品久久久久人妻精品| 欧美bdsm另类| 国产av麻豆久久久久久久| 色精品久久人妻99蜜桃| 国产精品福利在线免费观看| av在线老鸭窝| 久久6这里有精品| 男女啪啪激烈高潮av片| 一本精品99久久精品77| 搡老妇女老女人老熟妇| 精品人妻偷拍中文字幕| 亚洲天堂国产精品一区在线| 美女xxoo啪啪120秒动态图| 身体一侧抽搐| 色av中文字幕| 99在线人妻在线中文字幕| 国产精品1区2区在线观看.| 麻豆精品久久久久久蜜桃| 久久中文看片网| 人妻制服诱惑在线中文字幕| 成人性生交大片免费视频hd| 国产高清激情床上av| 日韩欧美在线乱码| 给我免费播放毛片高清在线观看| www日本黄色视频网| 中文字幕高清在线视频| 成人无遮挡网站| 自拍偷自拍亚洲精品老妇| 精品久久久久久久久久免费视频| 久久久久久久精品吃奶| 午夜福利在线在线| 男女那种视频在线观看| 亚洲国产精品成人综合色| 午夜免费成人在线视频| 亚洲在线观看片| 国产精品久久久久久av不卡| 自拍偷自拍亚洲精品老妇| 99热精品在线国产| 在线观看午夜福利视频| 97人妻精品一区二区三区麻豆| 可以在线观看毛片的网站| 男女之事视频高清在线观看| 日韩人妻高清精品专区| 午夜久久久久精精品| 午夜亚洲福利在线播放| 91午夜精品亚洲一区二区三区 | 欧美性猛交╳xxx乱大交人| 两个人的视频大全免费|