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

    基于CFD的46 000 t油船型線優(yōu)化設(shè)計(jì)

    2018-01-13 00:17:30鄭強(qiáng)馮佰威
    計(jì)算機(jī)輔助工程 2017年6期
    關(guān)鍵詞:優(yōu)化

    鄭強(qiáng)++馮佰威

    摘要: 為適應(yīng)國(guó)際海事組織船舶能效設(shè)計(jì)指數(shù)對(duì)綠色船舶的要求,以46 000 t油船為例,采用CAD軟件FRIENDSHIP對(duì)船舶進(jìn)行參數(shù)化建模,同時(shí)集成CFD軟件Shipflow,搭建船型優(yōu)化平臺(tái),實(shí)現(xiàn)參數(shù)化船型自動(dòng)生成、性能計(jì)算及優(yōu)化,獲得阻力性能較好的船型。研究結(jié)果表明:所采用的46 000 t油船型線優(yōu)化方法具有可行性,可為46 000 t油船型線設(shè)計(jì)提供參考。該研究方法可實(shí)現(xiàn)仿真驅(qū)動(dòng)設(shè)計(jì),對(duì)縮短設(shè)計(jì)周期具有重要意義。

    關(guān)鍵詞: 參數(shù)化建模; 油船型線; 優(yōu)化; 總阻力

    中圖分類號(hào): U441文獻(xiàn)標(biāo)志碼: B

    Optimization design of 46 000 t oil tanker line based on CFD

    ZHENG Qianga, FENG Baiweib

    (a.School of Transportation; b.Key Laboratory of High Performance Ship Technology

    (Wuhan University of Technology), Ministry of Education, Wuhan University of Technology, Wuhan 430063, China)

    Abstract: To meet the requirements of International Maritime Organization(IMO) energy efficiency design index for green ship, the parametric modeling of a 46 000 t oil tanker is carried out using CAD software FRIENDSHIP. The hull form optimization platform is built by integrating the CFD software Shipflow. The automatic generation of parametric hull form, the performance calculation and optimization are realized. The hull form with better resistance performance is obtained. The study results show that the line optimization method on the 46 000 t oil tanker is feasible. The simulationdriven design can be achieved using this research method. It is important for the time saving of the design cycle.

    Key words: parametric modeling; oil tanker line; optimization; total resistance

    收稿日期: 2017[KG*9〗07[KG*9〗20修回日期: 2017[KG*9〗09[KG*9〗08

    基金項(xiàng)目: 國(guó)家自然科學(xué)基金重點(diǎn)國(guó)際(地區(qū))合作研究項(xiàng)目(51720105011);國(guó)家自然科學(xué)基金(51479150,51709213);

    自主創(chuàng)新基金(173102007)

    作者簡(jiǎn)介: 鄭強(qiáng)(1994—),男,福建寧德人,碩士研究生,研究方向?yàn)榇岸鄬W(xué)科設(shè)計(jì)優(yōu)化技術(shù),(Email)451535308@qq.com

    馮佰威(1974—),男,遼寧北鎮(zhèn)人,副教授,博士,研究方向?yàn)榇皵?shù)字化和船舶多學(xué)科設(shè)計(jì)優(yōu)化,(Email)275951659@qq.com0引言

    傳統(tǒng)的船舶型線設(shè)計(jì)受限于船型數(shù)據(jù)和船舶設(shè)計(jì)者的經(jīng)驗(yàn),設(shè)計(jì)過(guò)程復(fù)雜繁瑣,資源和人力浪費(fèi)嚴(yán)重,得到的方案僅能滿足既定要求而非最優(yōu)設(shè)計(jì)。在科技發(fā)達(dá)的今天,傳統(tǒng)方法已經(jīng)不適應(yīng)船舶工業(yè)的發(fā)展,所以對(duì)船體型線的研究顯得尤為重要。當(dāng)前,計(jì)算機(jī)技術(shù)蓬勃發(fā)展,各行各業(yè)紛紛應(yīng)用計(jì)算機(jī)技術(shù)簡(jiǎn)化工作,船舶工業(yè)也不例外。計(jì)算機(jī)輔助系統(tǒng)在船舶型線設(shè)計(jì)中的作用越來(lái)越重要,常用的計(jì)算工具,如CAD和CFD等,已經(jīng)日漸成熟,越來(lái)越多的學(xué)者應(yīng)用這些軟件進(jìn)行優(yōu)化計(jì)算。但是,由于CAD和CFD之間缺少有效的集成,無(wú)法更好地發(fā)揮計(jì)算機(jī)的優(yōu)勢(shì),使得CFD仍不能很好地應(yīng)用于船舶設(shè)計(jì)開(kāi)發(fā)。

    本文圍繞船體的水動(dòng)力性能,利用參數(shù)化建模軟件FRIENDSHIP構(gòu)建船舶參數(shù)化模型,同時(shí)集成CFD計(jì)算程序Shipflow搭建船型優(yōu)化平臺(tái),實(shí)現(xiàn)參數(shù)化船型生成、性能計(jì)算和船型優(yōu)化。

    1船型參數(shù)化設(shè)計(jì)

    在參數(shù)化建模中,船體的曲面形狀由船體特征參數(shù)控制,特征參數(shù)對(duì)船舶水動(dòng)力性能影響很大。各類船舶有不同的特征參數(shù)和曲線,其中特征參數(shù)有船長(zhǎng)、型寬、設(shè)計(jì)吃水、主甲板高、棱形系數(shù)、方形系數(shù)等,特征曲線有龍骨線、平邊線、平底線、舭部升高線等。特征參數(shù)和特征曲線數(shù)據(jù)量要遠(yuǎn)遠(yuǎn)小于確定船體曲面的型值數(shù)據(jù)量。用特征參數(shù)、特征曲線表達(dá)船型,然后進(jìn)行初步的性能估算,可以減少計(jì)算量,提高計(jì)算速度。[13]船型參數(shù)化設(shè)計(jì)方法分為4步:(1)確定特征參數(shù);(2)設(shè)計(jì)縱向特征曲線;(3)生成橫剖面曲線;(4)得到船體曲面。[45]

    1.1特征參數(shù)

    特征參數(shù)控制特征曲線的形狀,進(jìn)而影響船體曲面形狀。特征參數(shù)可以分為2類:一類為全局參數(shù),如船體的垂線間長(zhǎng)度、船寬、型深、設(shè)計(jì)吃水、主甲板高度、棱形系數(shù)、方形系數(shù)等船體主要要素;另一類為局部參數(shù),主要控制船體縱向和橫向的特征曲線,如控制平邊線、平底線、舭部升高線和球艏等的參數(shù)。[6]endprint

    1.2特征曲線

    參數(shù)特征曲線可以分為縱向特征曲線和橫剖面曲線??v向特征曲線,如橫剖面面積曲線、設(shè)計(jì)水線、船體輪廓線、平邊線、平底線等,都是生成船體曲面必要的曲線。設(shè)計(jì)縱向特征曲線時(shí),可根據(jù)具體情況決定是設(shè)計(jì)全部縱向曲線還是僅設(shè)計(jì)部分縱向曲線??v向曲線主要組成見(jiàn)表1。縱向特征曲線通過(guò)控制橫剖面曲線上點(diǎn)的坐標(biāo),實(shí)現(xiàn)橫剖面曲線的設(shè)計(jì)。橫剖面曲線的創(chuàng)建一般用FRIENDSHIP中的Feature實(shí)現(xiàn)。

    1.3船體曲面的生成

    FRIENDSHIP軟件以形狀參數(shù)為基礎(chǔ),整個(gè)建模過(guò)程細(xì)分為3個(gè)連續(xù)步驟:首先,根據(jù)所選的特征參數(shù)建立參數(shù)點(diǎn),用軟件中的各種線型設(shè)計(jì)和生成光順的縱向特征曲線:然后,運(yùn)用軟件的相應(yīng)模塊生成一系列光順的橫剖面曲線;最后,利用蒙面法生成船體曲面。[710]船體曲面生成過(guò)程見(jiàn)圖1。

    2油船船體曲面全參數(shù)化建模

    以某46 000 t油船為例進(jìn)行船體曲面全參數(shù)化建模。由于后續(xù)黏流計(jì)算的需要,所建參數(shù)化模型按1∶27的比例進(jìn)行縮小。

    船體中部曲面的曲率變化較小,建模通常都從船體中部位置開(kāi)始,然后依次向艏、艉延伸。46 000 t油船在船體中部位置設(shè)有一定長(zhǎng)度的平行舯體,并且艏部帶有球鼻,因此以平行舯體段作為起始段,再以平行舯體兩端的橫剖面向艏、艉進(jìn)行延伸建模,主要將船體劃分為平行舯體段,船體前部,船體后部3個(gè)部分。根據(jù)不同部分曲面的特點(diǎn),對(duì)曲面進(jìn)行進(jìn)一步劃分,見(jiàn)圖2。

    2.1主要特征曲線創(chuàng)建

    根據(jù)船舶主體尺寸要素等數(shù)據(jù)建立船體的平邊線、平底線、中剖面輪廓線、船的設(shè)計(jì)吃水、橫剖面面積、設(shè)計(jì)水線、球艏等特征曲線。下面以球艏為例介紹特征曲線建立過(guò)程。

    球艏建模相關(guān)的主要特征參數(shù)見(jiàn)表2。

    2.2橫剖面構(gòu)造

    在特征曲線的基礎(chǔ)上用NURBS曲線表達(dá)光順的橫剖面曲線。由于船體橫剖面曲線隨船長(zhǎng)變化,需要將其劃分為幾個(gè)部分單獨(dú)建模,此處僅以船體中部曲面的創(chuàng)建為例詳細(xì)闡述,其他部位方法類似。

    船體中部曲面橫剖線分為底部、舭部和舷側(cè)3個(gè)部分,需要4個(gè)點(diǎn)控制橫剖面曲線形狀,見(jiàn)圖4。

    2.3曲面生成

    將Feature導(dǎo)入到Curve Engine中,設(shè)定好特征參數(shù)和曲線的配置,選擇橫剖面曲線生成的坐標(biāo)系統(tǒng),再將該Curve Engine添加到曲面生成器Meta Surface中,輸入曲面生成的起始和終止位置,就會(huì)在該范圍內(nèi)會(huì)生成無(wú)數(shù)的橫剖面從而得到光順的船體曲面,見(jiàn)圖5。

    3油船阻力性能優(yōu)化

    以參數(shù)化46 000 t油船為基礎(chǔ),選擇船體的首、尾部型線為優(yōu)化對(duì)象,建立相應(yīng)的數(shù)學(xué)模型。采用Sobol算法和Tsearch算法對(duì)數(shù)學(xué)模型進(jìn)行求解,完成46 000 t油船總阻力性能的數(shù)值優(yōu)化。

    3.1型值數(shù)據(jù)的自動(dòng)提取

    為實(shí)現(xiàn)FRIENDSHIP和Shipflow之間的集成,需要將參數(shù)化模型轉(zhuǎn)化為Shipflow所能讀取的格式,實(shí)現(xiàn)CFD計(jì)算和后續(xù)的船型自動(dòng)優(yōu)化。在研究過(guò)程中,定義Feature實(shí)現(xiàn)船體參數(shù)化曲面型值的自動(dòng)提取。根據(jù)曲面特點(diǎn)將曲面型值分為4個(gè)部分:艉軸、主體、艉部、球艏。采用做橫剖線的形式,對(duì)船體提取4部分型值,設(shè)置好Feature中的相關(guān)參數(shù),即可得到全船型值,見(jiàn)圖6。將該型值數(shù)據(jù)將輸入到SHIPFLOW中完成總阻力的計(jì)算。

    3.2計(jì)算網(wǎng)格

    對(duì)該46 000 t油船在設(shè)計(jì)航速(弗勞德數(shù)Fr=0.183)時(shí)阻力和流場(chǎng)進(jìn)行數(shù)值計(jì)算,興波阻力采用勢(shì)流面元法計(jì)算,在船體和自由表面進(jìn)行網(wǎng)格劃分,船體計(jì)算面元總數(shù)約為3 800個(gè)。自由液面大小取值:船前取0.5LPP(LPP為垂線間長(zhǎng)度),船后取1.6LPP,船側(cè)取0.6LPP,計(jì)算面元數(shù)約為3 500個(gè),見(jiàn)圖7。計(jì)算黏性阻力時(shí)采用RANS方法,自動(dòng)劃分網(wǎng)格,軸向×周向×徑向網(wǎng)格劃分為113×49×71,網(wǎng)格總數(shù)約為39萬(wàn)個(gè)。

    3.3優(yōu)化變量設(shè)定

    選擇船體的首、尾部型線為優(yōu)化對(duì)象,在船體的首部選擇10個(gè)參數(shù)、尾部選擇9個(gè)參數(shù),其相應(yīng)的變化范圍見(jiàn)表4。

    3.4優(yōu)化目標(biāo)定義

    重點(diǎn)對(duì)船體的首、尾部型線進(jìn)行優(yōu)化,以減小船舶的總阻力,因此將總阻力作為優(yōu)化的目標(biāo)函數(shù)。編寫(xiě)CFD計(jì)算所需的配置文件,型值數(shù)據(jù)按第3.1節(jié)所得到的型值。在FRIENDSHIP中集成Shipflow,設(shè)置好計(jì)算路徑進(jìn)行總阻力數(shù)值計(jì)算。

    3.5約束條件定義

    在保證46 000 t油船的阻力性能得到改善的同時(shí),要保持初始船型的營(yíng)運(yùn)性能。為保證船型在優(yōu)化前后的浮性和穩(wěn)性不發(fā)生太大的變化,將排水量和浮心的縱向位置作為優(yōu)化的約束條件,要求船體型線在約束條件范圍內(nèi)變化。選擇約束變化范圍為:(1)排水體積在母型船的基礎(chǔ)上浮動(dòng)范圍小于1%;(2)浮心縱向位置在母型船的基礎(chǔ)上浮動(dòng)范圍小于1%。

    3.6優(yōu)化算法設(shè)置

    選擇Sobol算法在全局范圍內(nèi)進(jìn)行較優(yōu)船型的探索,生成總計(jì)200組方案,從中選出總阻力最小的方案,并在此方案的基礎(chǔ)上進(jìn)行Tsearch算法優(yōu)化,最大迭代次數(shù)設(shè)置為200次,最終得到阻力最小的方案。

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

    經(jīng)過(guò)130 h的優(yōu)化,優(yōu)化結(jié)果見(jiàn)表5。優(yōu)化前后自由面波形輪廓圖及舷側(cè)縱向波切比較見(jiàn)圖8,優(yōu)化前后型線比較見(jiàn)圖9。

    優(yōu)化后總阻力下降5.60%。將優(yōu)化船與母型船的型線進(jìn)行比較可以發(fā)現(xiàn),排水體積的減少使得船體前部的橫剖面形狀由“U”形逐漸變?yōu)椤癡”形?!癡”形橫剖面形狀會(huì)使得船體的濕表面積減小,因此會(huì)使得船體的摩擦阻力相應(yīng)減小,同時(shí)由于“V”形橫剖面形狀使得船體的舭部漩渦減少,雖然“V”形剖面可能會(huì)使得船體的興波阻力增大,但由于該46 000 t油船屬于低速船,興波阻力所占的比重較小,因此興波阻力增大的幅度小于船體摩擦阻力增大的幅度。船體的尾部曲線略微收縮有向“V”形變化的趨勢(shì),船體的濕表面積減小從而使得摩擦阻力減小。endprint

    5結(jié)束語(yǔ)

    利用FRIENDSHIP和Shipflow成功實(shí)現(xiàn)46 000 t油船船型全參數(shù)化建模和總阻力自動(dòng)優(yōu)化,獲得阻力性能較好的船型,這對(duì)實(shí)現(xiàn)仿真驅(qū)動(dòng)設(shè)計(jì)、縮短設(shè)計(jì)周期具有重要意義。

    參考文獻(xiàn):

    [1]張萍. 船型參數(shù)化設(shè)計(jì)[D]. 無(wú)錫: 江南大學(xué), 2009.

    [2]張萍, 朱德祥, 何術(shù)龍. 參數(shù)化的船型設(shè)計(jì)方法[J]. 中國(guó)造船, 2008, 49(4): 2635.DOI: 10.3969/j.issn.10004882.2008.04.004.

    ZHANG P, ZHU D X, HE S L. Parametric design method of hull form[J]. Shipbuilding of China, 2008, 49(4): 2635. DOI: 10.3969/j.issn.10004882.2008.04.004.

    [3]黃金鋒, 萬(wàn)松林. 基于設(shè)計(jì)特征的FRIENDSHIP船型參數(shù)化方法及實(shí)現(xiàn)[J]. 中國(guó)艦船研究, 2012, 7(2): 7985. DOI: 10.3969/j.issn.16733185.2012.02.015.

    HUANG J F, WAN S L. Parametric method of FRIENDSHIP hull form based on design feature[J]. Chinese Journal of Ship Research, 2012, 7(2): 7985. DOI: 10.3969/j.issn.16733185.2012.02.015.

    [4]蔣國(guó)偉, 常海超, 馮佰威. 基于FRIENDSHIP的船型優(yōu)化方法研究[C]//2015年CAD/CAM學(xué)術(shù)交流會(huì)議論文集. 南京, 2015.

    [5]沈通, 馮佰威, 劉祖源, 等. 基于徑向基函數(shù)插值的船體曲面修改方法研究[J]. 中國(guó)造船, 2013, 54(4): 4554. DOI: 10.3969/j.issn.10004882.2013.04.006.

    SHEN T, FENG B W, LIU Z Y, et al. Research of hull form modification method based on interpolation of radial basis function[J]. Shipbuilding of China, 2013, 54(4): 4554. DOI: 10.3969/j.issn.10004882.2013.04.006.

    [6]胡春平. 基于FRIENDSHIP的船體型線優(yōu)化技術(shù)研究[D]. 武漢: 武漢理工大學(xué), 2012.

    [7]劉祖源, 馮佰威, 詹成勝. 船體型線多學(xué)科設(shè)計(jì)優(yōu)化[M]. 北京: 國(guó)防工業(yè)出版社, 2010.

    [8]馮佰威, 劉祖源, 詹成勝, 等. 基于船型修改融合方法的參數(shù)化建模技術(shù)[J]. 計(jì)算機(jī)輔助工程, 2010, 19(4): 37. DOI: 10.3969/j.issn.10060871.2010.04.002.

    FENG B W, LIU Z Y, ZHAN C S, et al. Parametrization modeling based on hull form morphing method[J]. Computer Aided Engineering, 2010, 19(4): 37. DOI: 10.3969/j.issn.10060871.2010.04.002.

    [9]YANG C, HUANG F X, KIM H. Hydrodynamic optimization of a triswach[J]. Journal of Hydrodynamics, 2014, 26(6): 856864.

    [10]馮佰威, 劉祖源, 常海超. 多學(xué)科設(shè)計(jì)優(yōu)化技術(shù)在船舶初步設(shè)計(jì)中的應(yīng)用[J]. 中國(guó)造船, 2009, 50(4): 109116. DOI: 10.3969/j.issn.10004882.2009.04.015.

    FENG B W, LIU Z Y, CHANG H C. Application of multidisciplinary design optimization techniques in ships preliminary design[J]. Shipbuilding of China, 2009, 50(4): 109116. DOI: 10.3969/j.issn.10004882.2009.04.015.(編輯武曉英)第26卷 第6期2017年12月計(jì) 算 機(jī) 輔 助 工 程Computer Aided EngineeringVol.26 No.6Dec. 2017endprint

    猜你喜歡
    優(yōu)化
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    PEMFC流道的多目標(biāo)優(yōu)化
    能源工程(2022年1期)2022-03-29 01:06:28
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    圍繞“地、業(yè)、人”優(yōu)化產(chǎn)業(yè)扶貧
    事業(yè)單位中固定資產(chǎn)會(huì)計(jì)處理的優(yōu)化
    4K HDR性能大幅度優(yōu)化 JVC DLA-X8 18 BC
    幾種常見(jiàn)的負(fù)載均衡算法的優(yōu)化
    電子制作(2017年20期)2017-04-26 06:57:45
    免费女性裸体啪啪无遮挡网站| 人人妻人人爽人人添夜夜欢视频| 欧美+日韩+精品| 成年美女黄网站色视频大全免费| 美女国产视频在线观看| 国产精品女同一区二区软件| 午夜日韩欧美国产| 亚洲精品国产av蜜桃| 在线观看人妻少妇| 久久久国产欧美日韩av| 国产精品久久久av美女十八| 一级爰片在线观看| 人妻 亚洲 视频| 日本-黄色视频高清免费观看| 亚洲美女视频黄频| 啦啦啦中文免费视频观看日本| 侵犯人妻中文字幕一二三四区| 国产成人精品一,二区| 在线天堂中文资源库| 免费人妻精品一区二区三区视频| 人人妻人人添人人爽欧美一区卜| 国产精品蜜桃在线观看| 五月开心婷婷网| 亚洲色图 男人天堂 中文字幕| 国产高清国产精品国产三级| 超碰97精品在线观看| 成年女人在线观看亚洲视频| 一级,二级,三级黄色视频| 中文精品一卡2卡3卡4更新| 午夜老司机福利剧场| 亚洲国产欧美日韩在线播放| 久久久国产精品麻豆| 美女主播在线视频| 少妇猛男粗大的猛烈进出视频| 亚洲精品中文字幕在线视频| 国产一区二区激情短视频 | 性色av一级| 一区二区日韩欧美中文字幕| 91国产中文字幕| av天堂久久9| 免费黄色在线免费观看| 国产伦理片在线播放av一区| 日本午夜av视频| 人人澡人人妻人| 国产亚洲一区二区精品| 国产一区二区三区综合在线观看| 国产麻豆69| 毛片一级片免费看久久久久| 国产亚洲av片在线观看秒播厂| 国产成人精品福利久久| av片东京热男人的天堂| 在线 av 中文字幕| 国产精品不卡视频一区二区| 69精品国产乱码久久久| 大片免费播放器 马上看| 丝瓜视频免费看黄片| 亚洲中文av在线| 美女主播在线视频| 亚洲成人手机| 亚洲国产精品国产精品| 国产免费福利视频在线观看| 丰满乱子伦码专区| 亚洲av综合色区一区| a级片在线免费高清观看视频| 国产黄色免费在线视频| 亚洲精品久久午夜乱码| 久久精品国产鲁丝片午夜精品| 午夜福利在线免费观看网站| 久久精品国产自在天天线| 少妇人妻久久综合中文| 免费高清在线观看日韩| av在线观看视频网站免费| 亚洲欧美成人综合另类久久久| 国产人伦9x9x在线观看 | 精品少妇黑人巨大在线播放| 国产精品av久久久久免费| 久久99一区二区三区| 观看av在线不卡| 国产男人的电影天堂91| 久久国产亚洲av麻豆专区| 欧美日韩av久久| 久久久久久久国产电影| 久久精品aⅴ一区二区三区四区 | 精品国产一区二区久久| 亚洲精品国产一区二区精华液| 一区二区三区四区激情视频| 人妻少妇偷人精品九色| 亚洲美女搞黄在线观看| 18禁国产床啪视频网站| 天堂中文最新版在线下载| 国产一区亚洲一区在线观看| 人妻系列 视频| 人妻一区二区av| 日本午夜av视频| 久久国产亚洲av麻豆专区| 乱人伦中国视频| 美女脱内裤让男人舔精品视频| 少妇的丰满在线观看| 在线 av 中文字幕| 欧美人与善性xxx| 国产熟女午夜一区二区三区| av电影中文网址| 久久国产精品男人的天堂亚洲| 成人毛片60女人毛片免费| 国产精品欧美亚洲77777| 侵犯人妻中文字幕一二三四区| 久久av网站| 哪个播放器可以免费观看大片| 国产日韩一区二区三区精品不卡| 免费日韩欧美在线观看| 国产亚洲精品第一综合不卡| 久久精品人人爽人人爽视色| 天天操日日干夜夜撸| 女人精品久久久久毛片| 亚洲欧美一区二区三区国产| 九草在线视频观看| 大码成人一级视频| 亚洲成人一二三区av| 亚洲国产日韩一区二区| 最近最新中文字幕免费大全7| 欧美 亚洲 国产 日韩一| 日本欧美视频一区| 91国产中文字幕| 国产免费视频播放在线视频| 亚洲国产看品久久| 亚洲av中文av极速乱| 青青草视频在线视频观看| 国产欧美日韩一区二区三区在线| 最近最新中文字幕大全免费视频 | 少妇人妻 视频| 一本久久精品| 91精品三级在线观看| av福利片在线| 午夜老司机福利剧场| 日韩电影二区| 日本午夜av视频| 日韩欧美精品免费久久| 久久毛片免费看一区二区三区| 成人手机av| 免费观看性生交大片5| 精品国产超薄肉色丝袜足j| 国产又色又爽无遮挡免| 国产有黄有色有爽视频| 日韩伦理黄色片| 亚洲精品久久午夜乱码| 老司机影院成人| 国产在线一区二区三区精| 亚洲精品美女久久久久99蜜臀 | 亚洲美女搞黄在线观看| 欧美+日韩+精品| freevideosex欧美| 男女边吃奶边做爰视频| 春色校园在线视频观看| www.自偷自拍.com| 国产探花极品一区二区| 美女大奶头黄色视频| 日韩三级伦理在线观看| 日韩欧美精品免费久久| 亚洲美女搞黄在线观看| 黄色一级大片看看| 爱豆传媒免费全集在线观看| 黄色 视频免费看| 久久国产精品大桥未久av| 又粗又硬又长又爽又黄的视频| 国产精品国产三级国产专区5o| 日韩熟女老妇一区二区性免费视频| 日本vs欧美在线观看视频| 中文乱码字字幕精品一区二区三区| 亚洲精华国产精华液的使用体验| 美女国产视频在线观看| 久久精品久久久久久久性| 久久国产精品大桥未久av| 男女边吃奶边做爰视频| 十八禁高潮呻吟视频| 男的添女的下面高潮视频| 综合色丁香网| 精品国产一区二区三区四区第35| av福利片在线| 亚洲国产精品成人久久小说| 国产欧美亚洲国产| 激情五月婷婷亚洲| 久久久a久久爽久久v久久| 中文字幕人妻丝袜制服| 国产97色在线日韩免费| 亚洲婷婷狠狠爱综合网| 另类精品久久| 亚洲精品一二三| 欧美精品一区二区大全| 成年女人毛片免费观看观看9 | 欧美成人精品欧美一级黄| 熟女av电影| 在线观看免费日韩欧美大片| 国产亚洲最大av| 久久久久国产精品人妻一区二区| 亚洲国产看品久久| 久久99一区二区三区| 国产成人av激情在线播放| 精品亚洲乱码少妇综合久久| av在线观看视频网站免费| 欧美日韩国产mv在线观看视频| 久久久欧美国产精品| 搡女人真爽免费视频火全软件| 久久久久精品性色| 久久女婷五月综合色啪小说| 在线观看美女被高潮喷水网站| 久久久久久久亚洲中文字幕| 99香蕉大伊视频| 99久久精品国产国产毛片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品久久久久久精品电影小说| 十八禁网站网址无遮挡| 日韩精品有码人妻一区| 久久国产亚洲av麻豆专区| 激情视频va一区二区三区| h视频一区二区三区| 永久网站在线| 搡女人真爽免费视频火全软件| 国产成人免费无遮挡视频| 婷婷成人精品国产| 亚洲在久久综合| av免费观看日本| 啦啦啦中文免费视频观看日本| 久久久久久久大尺度免费视频| 在线天堂中文资源库| 视频区图区小说| 国产精品久久久久久av不卡| 久久精品国产亚洲av天美| 国产精品一区二区在线观看99| 亚洲国产av影院在线观看| 久久精品国产综合久久久| 精品亚洲乱码少妇综合久久| 国产av精品麻豆| 少妇的丰满在线观看| 蜜桃国产av成人99| 亚洲精品久久久久久婷婷小说| 人妻人人澡人人爽人人| 久久亚洲国产成人精品v| 国产成人免费观看mmmm| 亚洲欧美日韩另类电影网站| 国产精品二区激情视频| 男人操女人黄网站| 91在线精品国自产拍蜜月| 中文字幕制服av| 丝袜美腿诱惑在线| 菩萨蛮人人尽说江南好唐韦庄| av电影中文网址| 在线观看人妻少妇| 成人毛片60女人毛片免费| 国产精品三级大全| 精品少妇一区二区三区视频日本电影 | 美女主播在线视频| 人妻系列 视频| 日本-黄色视频高清免费观看| 啦啦啦中文免费视频观看日本| 久久精品aⅴ一区二区三区四区 | 久久精品国产a三级三级三级| 国产极品天堂在线| 王馨瑶露胸无遮挡在线观看| 国产xxxxx性猛交| a级毛片在线看网站| 色视频在线一区二区三区| 天天躁夜夜躁狠狠躁躁| 欧美最新免费一区二区三区| 一边亲一边摸免费视频| 色播在线永久视频| 成年av动漫网址| 成人亚洲欧美一区二区av| 亚洲av福利一区| 天天躁夜夜躁狠狠躁躁| 国产熟女午夜一区二区三区| 在线观看www视频免费| 久久精品国产亚洲av高清一级| 91在线精品国自产拍蜜月| 少妇被粗大的猛进出69影院| 啦啦啦在线观看免费高清www| 自线自在国产av| 亚洲精品久久久久久婷婷小说| av不卡在线播放| 国产一区二区激情短视频 | 老汉色av国产亚洲站长工具| 国产精品女同一区二区软件| 亚洲第一青青草原| 多毛熟女@视频| 国产精品蜜桃在线观看| 亚洲欧美精品自产自拍| 中国国产av一级| 午夜精品国产一区二区电影| 精品国产乱码久久久久久男人| 各种免费的搞黄视频| 99久国产av精品国产电影| 80岁老熟妇乱子伦牲交| 如日韩欧美国产精品一区二区三区| 大片免费播放器 马上看| 男女无遮挡免费网站观看| 1024视频免费在线观看| 久久久久精品性色| 永久免费av网站大全| 一级毛片我不卡| 亚洲中文av在线| 欧美日韩视频高清一区二区三区二| 精品一品国产午夜福利视频| 午夜激情久久久久久久| 亚洲第一av免费看| 精品久久蜜臀av无| 亚洲第一av免费看| 亚洲国产欧美在线一区| tube8黄色片| 精品一区二区三区四区五区乱码 | 久久99热这里只频精品6学生| 久久毛片免费看一区二区三区| 一级片'在线观看视频| 精品亚洲成a人片在线观看| 中文字幕人妻熟女乱码| 国产一级毛片在线| 哪个播放器可以免费观看大片| 黄网站色视频无遮挡免费观看| 中文天堂在线官网| 午夜激情av网站| 亚洲av日韩在线播放| 性色avwww在线观看| 国产一区二区激情短视频 | 免费观看av网站的网址| 美国免费a级毛片| 久久久国产精品麻豆| 看免费av毛片| 捣出白浆h1v1| av在线播放精品| 中文字幕人妻丝袜一区二区 | 激情视频va一区二区三区| 视频在线观看一区二区三区| 99热全是精品| 亚洲精品,欧美精品| 亚洲欧美日韩另类电影网站| 汤姆久久久久久久影院中文字幕| 青春草国产在线视频| 1024视频免费在线观看| 99热全是精品| 街头女战士在线观看网站| 日韩一区二区三区影片| 成人手机av| 建设人人有责人人尽责人人享有的| 母亲3免费完整高清在线观看 | 亚洲成人av在线免费| 少妇猛男粗大的猛烈进出视频| 岛国毛片在线播放| 啦啦啦视频在线资源免费观看| 国产探花极品一区二区| 高清不卡的av网站| 久久久久久久久久久免费av| 丰满少妇做爰视频| 免费女性裸体啪啪无遮挡网站| av电影中文网址| 美女脱内裤让男人舔精品视频| 王馨瑶露胸无遮挡在线观看| 国产av码专区亚洲av| 亚洲国产欧美网| 高清视频免费观看一区二区| 精品少妇久久久久久888优播| 在线精品无人区一区二区三| 国精品久久久久久国模美| 捣出白浆h1v1| 毛片一级片免费看久久久久| 美女脱内裤让男人舔精品视频| 久久精品国产综合久久久| 免费高清在线观看日韩| 一区二区三区乱码不卡18| 色网站视频免费| 丰满少妇做爰视频| 国产成人精品久久久久久| 人妻人人澡人人爽人人| 大话2 男鬼变身卡| 丝袜脚勾引网站| 亚洲一码二码三码区别大吗| 亚洲欧美精品综合一区二区三区 | 99久久精品国产国产毛片| 狠狠婷婷综合久久久久久88av| 日韩欧美精品免费久久| 亚洲欧美一区二区三区黑人 | 免费人妻精品一区二区三区视频| 成人二区视频| 亚洲欧洲国产日韩| av在线app专区| 色吧在线观看| 久久精品熟女亚洲av麻豆精品| 一边亲一边摸免费视频| 蜜桃在线观看..| 国产精品一区二区在线不卡| 不卡视频在线观看欧美| 亚洲精品久久午夜乱码| 人人妻人人澡人人爽人人夜夜| 国产av国产精品国产| 午夜免费观看性视频| 日日撸夜夜添| 国产成人精品一,二区| 成人国产麻豆网| 麻豆乱淫一区二区| 18在线观看网站| 热re99久久精品国产66热6| 波多野结衣一区麻豆| 欧美日韩亚洲高清精品| 日本vs欧美在线观看视频| 丁香六月天网| 熟女少妇亚洲综合色aaa.| 久久婷婷青草| 美女xxoo啪啪120秒动态图| 在线精品无人区一区二区三| av卡一久久| 丝袜美腿诱惑在线| 下体分泌物呈黄色| 欧美日韩成人在线一区二区| 中文字幕制服av| 亚洲国产av新网站| 在线观看一区二区三区激情| 国产精品三级大全| 高清视频免费观看一区二区| av卡一久久| 欧美激情高清一区二区三区 | 国产精品香港三级国产av潘金莲 | 欧美日韩亚洲国产一区二区在线观看 | 久久综合国产亚洲精品| 看免费成人av毛片| 日韩欧美精品免费久久| av网站免费在线观看视频| 少妇人妻精品综合一区二区| 波多野结衣一区麻豆| 自线自在国产av| 十分钟在线观看高清视频www| 精品一区二区三卡| 在线观看国产h片| 国产精品二区激情视频| 亚洲婷婷狠狠爱综合网| 男女免费视频国产| 秋霞在线观看毛片| 最黄视频免费看| 日日啪夜夜爽| 日本猛色少妇xxxxx猛交久久| 久久毛片免费看一区二区三区| 在线观看美女被高潮喷水网站| 伦理电影大哥的女人| 婷婷成人精品国产| 成人国产麻豆网| 波多野结衣av一区二区av| 免费久久久久久久精品成人欧美视频| 人妻一区二区av| 久久 成人 亚洲| 精品国产一区二区三区久久久樱花| 2021少妇久久久久久久久久久| 2018国产大陆天天弄谢| 高清在线视频一区二区三区| 一级片'在线观看视频| 亚洲精华国产精华液的使用体验| 啦啦啦视频在线资源免费观看| 国产有黄有色有爽视频| 少妇 在线观看| 日本色播在线视频| 国产日韩欧美亚洲二区| 黑人猛操日本美女一级片| 亚洲,欧美精品.| 国产免费视频播放在线视频| 大话2 男鬼变身卡| 精品国产超薄肉色丝袜足j| 又大又黄又爽视频免费| 久久国产精品男人的天堂亚洲| 亚洲欧洲精品一区二区精品久久久 | 久久国产精品男人的天堂亚洲| 18禁裸乳无遮挡动漫免费视频| 国产日韩欧美亚洲二区| 少妇人妻久久综合中文| 久久人人97超碰香蕉20202| 国产精品一国产av| 97人妻天天添夜夜摸| 五月开心婷婷网| 亚洲精品久久午夜乱码| 久久av网站| 久久影院123| 日韩在线高清观看一区二区三区| 国产精品欧美亚洲77777| 人妻 亚洲 视频| 2018国产大陆天天弄谢| 精品一区二区三卡| 黄频高清免费视频| 多毛熟女@视频| 亚洲人成网站在线观看播放| 国产精品欧美亚洲77777| 热99久久久久精品小说推荐| √禁漫天堂资源中文www| 国产又爽黄色视频| 日韩大片免费观看网站| 精品99又大又爽又粗少妇毛片| 久久久久精品久久久久真实原创| 成年av动漫网址| 色婷婷av一区二区三区视频| 极品少妇高潮喷水抽搐| 又粗又硬又长又爽又黄的视频| 两个人免费观看高清视频| 亚洲精品国产色婷婷电影| 国产又色又爽无遮挡免| 午夜免费鲁丝| 亚洲精品美女久久av网站| 免费观看无遮挡的男女| 女性生殖器流出的白浆| 好男人视频免费观看在线| 波野结衣二区三区在线| 男男h啪啪无遮挡| 99国产综合亚洲精品| av不卡在线播放| 日韩av不卡免费在线播放| 久久久精品区二区三区| 欧美bdsm另类| 免费黄网站久久成人精品| 婷婷色综合大香蕉| 女人被躁到高潮嗷嗷叫费观| 丝袜人妻中文字幕| 欧美日韩av久久| a级毛片在线看网站| 在线免费观看不下载黄p国产| 国产极品天堂在线| 人妻系列 视频| 如日韩欧美国产精品一区二区三区| 精品99又大又爽又粗少妇毛片| 69精品国产乱码久久久| 亚洲国产精品国产精品| 色吧在线观看| 国产淫语在线视频| 亚洲精品国产一区二区精华液| 国产av一区二区精品久久| 在线观看美女被高潮喷水网站| 女人久久www免费人成看片| 9191精品国产免费久久| 国产免费现黄频在线看| 大片免费播放器 马上看| 黄频高清免费视频| 狂野欧美激情性bbbbbb| 国产精品成人在线| 天天躁夜夜躁狠狠躁躁| 一本—道久久a久久精品蜜桃钙片| 精品国产一区二区三区四区第35| 久久久久人妻精品一区果冻| 人妻一区二区av| 色播在线永久视频| 中文字幕另类日韩欧美亚洲嫩草| 性色avwww在线观看| 亚洲一码二码三码区别大吗| 亚洲欧美日韩另类电影网站| 91aial.com中文字幕在线观看| 久久人人97超碰香蕉20202| 桃花免费在线播放| 美女午夜性视频免费| 国产精品 欧美亚洲| 免费在线观看完整版高清| 看免费av毛片| 色94色欧美一区二区| 不卡av一区二区三区| 麻豆av在线久日| 超碰97精品在线观看| 少妇人妻久久综合中文| 日韩精品免费视频一区二区三区| 色哟哟·www| 国产精品国产av在线观看| 最近中文字幕高清免费大全6| 人体艺术视频欧美日本| 亚洲国产色片| freevideosex欧美| 亚洲精品中文字幕在线视频| 99香蕉大伊视频| 亚洲精品一区蜜桃| 国产一区亚洲一区在线观看| 少妇的逼水好多| 日本av手机在线免费观看| 免费看不卡的av| av免费在线看不卡| 久久国产精品男人的天堂亚洲| 婷婷色综合大香蕉| 最近的中文字幕免费完整| 亚洲av欧美aⅴ国产| 一区二区三区激情视频| 婷婷色综合大香蕉| 丝袜喷水一区| av网站在线播放免费| 亚洲视频免费观看视频| 黄色怎么调成土黄色| 欧美激情极品国产一区二区三区| 亚洲四区av| 欧美日韩成人在线一区二区| 18禁动态无遮挡网站| 极品人妻少妇av视频| 男人添女人高潮全过程视频| 黑人猛操日本美女一级片| 欧美黄色片欧美黄色片| 国产人伦9x9x在线观看 | 视频区图区小说| 天天躁狠狠躁夜夜躁狠狠躁| 一级毛片电影观看| 亚洲精品美女久久久久99蜜臀 | 街头女战士在线观看网站| 亚洲精品国产av蜜桃| 中文字幕色久视频| 国产又色又爽无遮挡免| 免费观看性生交大片5| 欧美日韩综合久久久久久| 国产精品久久久久久精品电影小说| 亚洲精品,欧美精品| 国产在线一区二区三区精| 久久久国产一区二区| 18禁国产床啪视频网站| 久久 成人 亚洲| 制服丝袜香蕉在线| 亚洲三区欧美一区| 一区二区av电影网| 18禁国产床啪视频网站| 女性被躁到高潮视频| 精品一区二区免费观看| 亚洲精品美女久久久久99蜜臀 | 日本欧美国产在线视频| 搡女人真爽免费视频火全软件|