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

    中高速船航行姿態(tài)及阻力數(shù)值預(yù)報方法研究*

    2016-08-29 04:35:41周廣禮歐勇鵬
    關(guān)鍵詞:興波船模船體

    周廣禮 歐勇鵬

    (海軍工程大學(xué)艦船工程系 武漢 430033)

    ?

    中高速船航行姿態(tài)及阻力數(shù)值預(yù)報方法研究*

    周廣禮歐勇鵬

    (海軍工程大學(xué)艦船工程系武漢430033)

    針對中高速船不同航速下航態(tài)變化顯著的特點(diǎn),基于RANS方程與重疊網(wǎng)格技術(shù),提出了船舶姿態(tài)及阻力的數(shù)值預(yù)報方法.以DTMB 5415為研究對象,開展了傅氏數(shù)Fr=0.28及0.41下船模航行狀態(tài)的數(shù)值模擬,計算結(jié)果表明,研究航速下,船模升沉及縱傾的預(yù)報誤差可控制在4%及10%以內(nèi),阻力預(yù)報誤差小于3%,同時該方法能夠較為準(zhǔn)確地捕捉船體興波流場的細(xì)節(jié)信息,可為中高速船快速性預(yù)報提供可行途徑.

    中高速船;重疊網(wǎng)格;計算流體力學(xué);航行姿態(tài);阻力

    0 引  言

    排水型船高速航行時,船體表面壓力分布與低航速下相比存在較大不同,易出現(xiàn)顯著的升沉及縱傾,從而增加了船舶流場及阻力的預(yù)報難度.以往船舶航行姿態(tài)主要采用基于Rankine函數(shù)的勢流方法進(jìn)行求解[1-3],該方法忽略了流體粘性,其適用范圍有限.

    隨著計算機(jī)性能的提高及網(wǎng)格技術(shù)的發(fā)展,應(yīng)用粘流CFD方法預(yù)報船舶姿態(tài)已逐步得以實(shí)現(xiàn).倪崇本等[4-5]分別應(yīng)用動態(tài)變形網(wǎng)格技術(shù)對高速三體船及INSEAN 2340標(biāo)模的阻力性能進(jìn)行了計算,探索了計及航態(tài)下船舶阻力的預(yù)報方法,然而該方法基于彈性節(jié)點(diǎn)重構(gòu)技術(shù)對船體周圍的網(wǎng)格進(jìn)行更新,對于姿態(tài)幅度變化較大的工況,往往會由于網(wǎng)格扭曲程度過大而影響計算的收斂性;郭春雨等[6]則以KCS船模為對象,提出了基于整體動網(wǎng)格技術(shù)的船舶阻力預(yù)報方法,應(yīng)用此方法雖較好地解決了網(wǎng)格更新的問題,但求解船舶大幅縱傾升沉運(yùn)動時,水線面加密域較大,增加了網(wǎng)格數(shù)量,降低了計算效率.

    近年來,重疊網(wǎng)格技術(shù)在求解船舶運(yùn)動方面展現(xiàn)出廣闊的應(yīng)用前景[7-9],基于該技術(shù),Wan D等[10-12]實(shí)現(xiàn)了KCS與DTMB 5512等多艘標(biāo)模的自航及復(fù)雜操縱運(yùn)動,極大地推動了船舶計算流體力學(xué)的應(yīng)用與發(fā)展,但目前國內(nèi)應(yīng)用重疊網(wǎng)格開展船舶運(yùn)動計算的相關(guān)研究仍較為滯后,且融合其他模型對船舶運(yùn)動進(jìn)行具體分析時,尚需更多的算例進(jìn)行對比驗(yàn)證.本文基于Star CCM+軟件平臺,應(yīng)用重疊網(wǎng)格技術(shù),開展了不同航速下DTMB 5415船模航行姿態(tài)及阻力性能的數(shù)值預(yù)報,為船舶快速性評定提供了可行方法.

    1 數(shù)值方法

    1.1控制方程及湍流模型

    文中以RANS方程作為求解船舶航行姿態(tài)及阻力的基本方程,其具體形式如下.

    (1)

    式中:ρ為流體密度;μ為流體粘度;p為靜壓;fi為單位質(zhì)量的質(zhì)量力;ui、uf為速度分量.

    湍流模式選取SSTk-ω模型,詳細(xì)推導(dǎo)過程參見文獻(xiàn)[13],下面僅給出該模型的數(shù)學(xué)表達(dá)式.

    (2)

    (3)

    式中:Γk和Γω為k和ω的有效擴(kuò)散率;Gk為平均速度梯度產(chǎn)生的湍流動能;Gω為耗散ω產(chǎn)生的湍流動能;Yk和Yω為湍流k和ω的耗散;Sk和Sω為用戶自定義源項(xiàng).

    1.2重疊網(wǎng)格技術(shù)

    應(yīng)用重疊網(wǎng)格模擬物體復(fù)雜空間運(yùn)動時,將計算域分為多個子域,并分別對子域進(jìn)行網(wǎng)格離散,通過域間空間插值可實(shí)現(xiàn)流場信息的交互,較動態(tài)變形網(wǎng)格技術(shù),其簡化了網(wǎng)格拓?fù)浣Y(jié)構(gòu),增強(qiáng)了網(wǎng)格普適性,見圖1.

    圖1 重疊網(wǎng)格

    在進(jìn)行船舶運(yùn)動求解時,設(shè)定背景域相對大地坐標(biāo)系固定不動,同時在船體周圍建立運(yùn)動域,運(yùn)動域隨船運(yùn)動,通過挖洞計算實(shí)時確定域間邊界,背景域中被挖掉的網(wǎng)格不參與計算,邊界點(diǎn)外兩域重疊部分進(jìn)行數(shù)據(jù)插值完成流場信息的傳遞.

    2 船體姿態(tài)響應(yīng)模型

    為求解船體姿態(tài)的時歷變化過程,見圖2,定義大地坐標(biāo)系E-ξηζ及隨船坐標(biāo)系O-xyz.其中,隨船坐標(biāo)系原點(diǎn)位于船體質(zhì)心處,未釋放船體運(yùn)動時,隨船坐標(biāo)系與大地坐標(biāo)系完全重合,船體處于正浮狀態(tài).

    圖2 坐標(biāo)系定義

    固定來流航速下,船體周圍水動力發(fā)生變化,所受浮力及作用于重心的縱向力矩均偏離正浮平衡狀態(tài),釋放縱向及垂向自由度后,船體開始運(yùn)動,由于流體阻尼的作用,運(yùn)動幅度逐漸衰減,并趨于新的穩(wěn)定平衡狀態(tài).船體的縱向及垂向運(yùn)動方程可表達(dá)為

    (4)

    (5)

    式中:F與M與分別為船體所受浮力及相對質(zhì)心的縱向力矩;w為垂向速度;q為縱搖角速度,B=mw為船體垂向動量;K=Iyq為船體所受縱向動量矩.

    在計算過程中,通過對船體表面壓力進(jìn)行積分,求得不同時刻作用于船體上的力及力矩,結(jié)合上述運(yùn)動模型,即可求解船體運(yùn)動的時歷變化過程,待計算穩(wěn)定后便可獲取船體阻力及流場信息,具體計算流程見圖3.

    圖3 計算流程圖

    3 數(shù)值模擬及結(jié)果分析

    3.1計算對象

    DTMB 5415為ITTC(International Towing Conference)推薦的艦船試驗(yàn)?zāi)P停瑖鴥?nèi)外學(xué)者針對該模型開展了大量的阻力及流場測量試驗(yàn)[14],為CFD計算方法的可靠性驗(yàn)證提供了強(qiáng)有力的支撐,船模主參數(shù)見表1.

    表1 DTMB 5415模型主參數(shù)

    3.2計算域設(shè)置

    由于模型關(guān)于船體中縱剖面對稱,為減小計算量取一半船模開展計算,見圖4,運(yùn)動域設(shè)置為長寬高為1.0L×1.0B×0.3L的方形域,背景域大小及邊界條件設(shè)置如下:

    1) 入口入口距離船體首垂線1.5倍船長,邊界條件均設(shè)置為速度入口.

    2) 出口出口距離船體尾垂線3倍船長,邊界條件設(shè)置為壓力出口,壓力為未擾動時邊界壓力.

    3) 流域的上、下及右邊界設(shè)置無滑移、不可穿透壁面.

    4) 流域左邊界設(shè)置為對稱面.

    5) 船體表面定義為無滑移、不可穿透邊界條件,邊界條件設(shè)置為無滑移壁面.

    圖4 計算域大小及邊界條件定義

    3.3網(wǎng)格劃分

    文中采用切割體網(wǎng)格技術(shù)分別對背景域及運(yùn)動域進(jìn)行網(wǎng)格離散,為確保兩域之間數(shù)據(jù)傳遞的穩(wěn)定性,對域間重疊部分網(wǎng)格尺度進(jìn)行了均一化處理,同時為準(zhǔn)確地捕捉兩相流的自由液面,對船體吃水附近區(qū)域進(jìn)行網(wǎng)格加密,此外,在船體近壁面設(shè)定了局部加密區(qū),確保壁面y+在30~300之間,并對船體首尾處的網(wǎng)格進(jìn)行了細(xì)化處理,船體表面及域剖面網(wǎng)格見圖5,總網(wǎng)格數(shù)為263萬.

    圖5 船體表面及域剖面網(wǎng)格

    3.4結(jié)果與分析

    表2 姿態(tài)計算結(jié)果分析

    本文分別將釋放姿態(tài)前后船模阻力的計算結(jié)果穩(wěn)定值與試驗(yàn)值進(jìn)行了對比(見表3),由誤差分析了可知,釋放姿態(tài)前船模阻力計算誤差均大于計及航行姿態(tài)的工況,計及姿態(tài)下,船模阻力計算誤差均小于3%.

    表3 阻力計算結(jié)果分析

    圖6給出了不同傅氏數(shù)下船模興波波形圖,對比試驗(yàn)值可得,應(yīng)用文中方法可較為準(zhǔn)確地預(yù)報船模興波的演變規(guī)律.船體表面興波波高沿船長方向的變化見圖7,分析可知,F(xiàn)r=0.28時,計算所得船體表面興波與試驗(yàn)值較為吻合,而當(dāng)Fr=0.41時,在后船體表面興波較試驗(yàn)值略低,但仍能反映波高沿船長方向的發(fā)展規(guī)律.

    圖6 不同航速下船體興波圖

    圖7 不同航速下船體表面波高圖

    4 結(jié) 束 語

    基于Star-CCM+內(nèi)嵌的VOF模型及多自由度求解器,應(yīng)用重疊網(wǎng)格技術(shù)開展了DTMB 5415船模不同航速下航行姿態(tài)及阻力的數(shù)值預(yù)報.分析計算結(jié)果可知,文中方法能夠較為準(zhǔn)確地預(yù)報中高速下船模的航行狀態(tài)及流場信息,可為船舶快速性預(yù)報提供可行途徑.

    [1]YASUKAWA H.A rankine panel method to calculate steady wave-making resistance of a ship taking the effect of sinkage and trim into account[J].Transactions of the West-Japan Society of Naval Architects,1993(86):27-35.

    [2]王中,盧曉平,王瑋.考慮升沉和縱傾的方尾船非線性興波阻力計算[J].水動力學(xué)研究與進(jìn)展,2010,25(3):422-428.

    [3]王中,盧曉平,付攀.三體船升沉和縱傾計算及其對興波阻力的影響[J].上海交通大學(xué)學(xué)報,2010,44(10):1388-1392.

    [4]倪崇本,朱仁傳,繆國平,等.計及航行姿態(tài)變化的高

    速多體船阻力預(yù)報[J].水動力學(xué)研究與進(jìn)展,2011,26(1):101-107.

    [5]YAO C B,DONG W Ci.Method to calculate resistance of high-speed displacement ship taking the effect of dynamic sinkage and trim and fluid viscosity into account[J].Journal of Shanghai Jiaotong University,2012,17(4):421-426.

    [6]郭春雨,王戀舟,趙慶新,等.一種計及姿態(tài)變化的船舶阻力預(yù)報方法[J].船舶工程,2015,37(1):31-34.

    [7]CARRICA P M,ISMAIL F,HYMAN M,et al.Turn and zigzag maneuvers of a surface combatant using a URANS approach with dynamic overset grids[J].Journal of Marine Science and Technology,2013(18):166-181.

    [8]CARRICA P M,F(xiàn)U H,STERN F.Computations of self-propulsion free to trim and of motions in waves of the KRISO Container Ship (KCS) model[J].Applied Ocean Research,2011(33):309-320.

    [9] CARRICA P M,CASTRO A M,STERN F.Self-propulsion computations using a speed controller and a discretized propeller with dynamic overset grids[J].Journal of Marine Science and Technology,2010(15):316-330.

    [10]WAN D,SHEN Z,CARRICA P M.Dynamic overset grids in OpenFOAM with application to KCS self-propulsion and maneuvering[J].Ocean Engineering,2015(108):287-306.

    [11] CARRIC P M,WILSON R V,NOACK R W,et al.Ship motions using single-phase level set with dynamic overset grids[J].Computers & Fluids,2007(36):1415-1433.

    [12]TEZDOGAN T,DEMIREL Y K,KELLETT P,et al.Full-scale unsteady RANS CFD simulations of ship behaviour and performance in head seas due to slow steaming[J].Ocean Engineering,2015(97):180-206.

    [13]LARSSON L,RAVEN H.Ship resistance and flow[M].New York:Society of Naval Architects and Marine Engineering,2010.

    [14]OLIVIERI A,PISTANI F,AVANZINI A,et al.Towing tank experiments of resistance,sinkage and trim,boundary layer,wake,and free surface flow around a naval combatant INSEAN 2340 model[R].Iowa,USA:University of Iowa,2001.

    Study on the Numerical Prediction Method of Calculating the Sailing State and Resistance for Medium & High Speed Ship

    ZHOU GuangliOU Yongpeng

    (DepartmentofNavalArchitecture,NavalUniversityofEngineering,Wuhan430033,China)

    The sailing state can be remarkably different when a displacement ship advancing with various high speeds in calm water. A method is presented to calculate the sailing state and resistance for medium & high speed ship based on Reynolds Averaged Navier-Stokes (RANS) equations and overset mesh. The surface combatant DTMB 5415 is chosen as an example to verify the method’s feasibility by simulating the ship’s hydrodynamic performance withFr=0.28 and 0.41. The results show that the errors of the calculated sinkages and trims can be controlled within 4% to 10% and the errors of the calculated resistances are less than 3%. Meanwhile, the current method can more accurately capture the ship wave profiles, which provides a feasible tool for the prediction of sailing states and resistances of medium & high speed ship.

    medium & high speed ship; overset mesh; CFD; sailing state; resistance

    2015-12-28

    U662.3

    10.3963/j.issn.2095-3844.2016.04.014

    周廣禮(1990- ):男,博士生,主要研究領(lǐng)域?yàn)榇八畡恿W(xué)

    *廣東省教育部產(chǎn)學(xué)研結(jié)合項(xiàng)目(2012B091000137)、總裝水動力重點(diǎn)基金資助項(xiàng)目(9140A14030712JB11044)資助

    猜你喜歡
    興波船模船體
    基于模糊PID的船模航向控制研究
    船體行駛過程中的壓力監(jiān)測方法
    船舶興波運(yùn)動的非線性現(xiàn)象算法研究與仿真
    基于Friendship的高速客船特征參數(shù)靈敏度分析
    船模靜水橫搖試驗(yàn)的不確定度分析
    焊接殘余應(yīng)力對船體結(jié)構(gòu)疲勞強(qiáng)度的影響分析
    焊接(2015年9期)2015-07-18 11:03:51
    赴美軍“仁慈”號醫(yī)院船駐船體會
    西洋船模王——童鑑良
    航海(2014年6期)2014-12-12 10:36:03
    水下爆炸氣泡作用下船體總縱強(qiáng)度估算方法
    影響潛航體水面興波因素分析
    欧美日韩乱码在线| 首页视频小说图片口味搜索| 午夜久久久久精精品| 美国免费a级毛片| 久久精品亚洲精品国产色婷小说| 99久久99久久久精品蜜桃| 婷婷六月久久综合丁香| 国产又爽黄色视频| 在线免费观看的www视频| 久久久国产精品麻豆| 淫秽高清视频在线观看| 日韩视频一区二区在线观看| 精品久久久久久久末码| 日韩一卡2卡3卡4卡2021年| 久久久久久亚洲精品国产蜜桃av| 正在播放国产对白刺激| 国产精品电影一区二区三区| 欧美绝顶高潮抽搐喷水| 亚洲熟女毛片儿| 成人永久免费在线观看视频| 国产成人精品久久二区二区免费| 精品卡一卡二卡四卡免费| 看黄色毛片网站| 国产亚洲精品综合一区在线观看 | 亚洲成国产人片在线观看| 国产成人av教育| 精品无人区乱码1区二区| 精品无人区乱码1区二区| 99久久综合精品五月天人人| 男男h啪啪无遮挡| 不卡一级毛片| 哪里可以看免费的av片| 久久久久久大精品| 国产男靠女视频免费网站| 国产精品98久久久久久宅男小说| 色哟哟哟哟哟哟| 欧美三级亚洲精品| 身体一侧抽搐| 午夜视频精品福利| 国产精品综合久久久久久久免费| 成人国产一区最新在线观看| 免费看a级黄色片| 久久欧美精品欧美久久欧美| 天天躁夜夜躁狠狠躁躁| 好男人电影高清在线观看| 成人国产一区最新在线观看| 欧美在线黄色| 在线免费观看的www视频| 午夜福利高清视频| 亚洲va日本ⅴa欧美va伊人久久| 熟妇人妻久久中文字幕3abv| 国产视频一区二区在线看| 久久欧美精品欧美久久欧美| 亚洲片人在线观看| 十八禁网站免费在线| 日韩av在线大香蕉| 一本久久中文字幕| 久99久视频精品免费| 日韩一卡2卡3卡4卡2021年| 两人在一起打扑克的视频| 特大巨黑吊av在线直播 | 国产精华一区二区三区| 国产区一区二久久| 又黄又爽又免费观看的视频| 两个人免费观看高清视频| 欧美成人性av电影在线观看| 亚洲国产精品久久男人天堂| 亚洲人成网站在线播放欧美日韩| 99精品欧美一区二区三区四区| 一本久久中文字幕| 一级毛片精品| 午夜影院日韩av| 午夜福利一区二区在线看| 精华霜和精华液先用哪个| 久久久久久久久中文| 国产亚洲精品一区二区www| 国产精品国产高清国产av| 一进一出抽搐gif免费好疼| 99国产极品粉嫩在线观看| 国产精品香港三级国产av潘金莲| 日韩精品免费视频一区二区三区| 欧美黑人欧美精品刺激| 亚洲av第一区精品v没综合| 99riav亚洲国产免费| 很黄的视频免费| 女人被狂操c到高潮| 成人精品一区二区免费| 国产成+人综合+亚洲专区| 亚洲 欧美 日韩 在线 免费| 神马国产精品三级电影在线观看 | 妹子高潮喷水视频| 午夜福利在线在线| 国产在线观看jvid| 香蕉av资源在线| 日韩视频一区二区在线观看| 久久精品成人免费网站| 亚洲欧美激情综合另类| 狂野欧美激情性xxxx| 91麻豆av在线| 黄片播放在线免费| 午夜日韩欧美国产| 中文字幕精品免费在线观看视频| 亚洲成人免费电影在线观看| 人人妻人人澡欧美一区二区| 久久亚洲精品不卡| 欧美激情极品国产一区二区三区| 人成视频在线观看免费观看| 久久午夜亚洲精品久久| 久久久精品国产亚洲av高清涩受| 亚洲精品av麻豆狂野| 不卡av一区二区三区| 亚洲精品国产区一区二| 午夜影院日韩av| 长腿黑丝高跟| 日韩av在线大香蕉| 亚洲欧美激情综合另类| 亚洲欧洲精品一区二区精品久久久| 中文字幕最新亚洲高清| 侵犯人妻中文字幕一二三四区| 欧美+亚洲+日韩+国产| 亚洲精品久久成人aⅴ小说| 国产av一区二区精品久久| 听说在线观看完整版免费高清| 最近最新免费中文字幕在线| 在线观看免费午夜福利视频| 两人在一起打扑克的视频| 国内精品久久久久精免费| 成年版毛片免费区| 亚洲真实伦在线观看| 天堂动漫精品| а√天堂www在线а√下载| 99精品久久久久人妻精品| 欧美大码av| 曰老女人黄片| 欧美性长视频在线观看| 国产成人欧美| 午夜福利18| 男女床上黄色一级片免费看| 精品不卡国产一区二区三区| 身体一侧抽搐| 香蕉av资源在线| 国产亚洲精品久久久久5区| 给我免费播放毛片高清在线观看| 1024手机看黄色片| 欧美黄色淫秽网站| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美精品综合久久99| 国产高清videossex| 精品高清国产在线一区| 亚洲午夜精品一区,二区,三区| 女同久久另类99精品国产91| 伊人久久大香线蕉亚洲五| 在线永久观看黄色视频| 露出奶头的视频| 亚洲国产欧美一区二区综合| 精品一区二区三区四区五区乱码| 国产精品一区二区免费欧美| 亚洲av五月六月丁香网| 亚洲成人精品中文字幕电影| 麻豆av在线久日| 中文在线观看免费www的网站 | 他把我摸到了高潮在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久免费精品人妻一区二区 | 中文字幕精品亚洲无线码一区 | 亚洲av日韩精品久久久久久密| 欧美丝袜亚洲另类 | 黄色毛片三级朝国网站| av中文乱码字幕在线| 亚洲国产精品成人综合色| 日韩欧美国产在线观看| 动漫黄色视频在线观看| 一个人免费在线观看的高清视频| 国产精品亚洲av一区麻豆| 成人国产综合亚洲| 久久天堂一区二区三区四区| 好男人在线观看高清免费视频 | 18禁美女被吸乳视频| √禁漫天堂资源中文www| 精品国产超薄肉色丝袜足j| 国产精品九九99| 国产日本99.免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日本视频| 国产精品免费一区二区三区在线| 99在线视频只有这里精品首页| 可以在线观看毛片的网站| 日韩欧美国产一区二区入口| 亚洲五月天丁香| 91老司机精品| 丁香欧美五月| 视频在线观看一区二区三区| 啦啦啦韩国在线观看视频| 亚洲成人精品中文字幕电影| 真人一进一出gif抽搐免费| 亚洲中文日韩欧美视频| 欧美黑人精品巨大| 夜夜夜夜夜久久久久| 日本成人三级电影网站| 久久中文字幕人妻熟女| 啦啦啦韩国在线观看视频| 欧美在线黄色| 日韩一卡2卡3卡4卡2021年| 侵犯人妻中文字幕一二三四区| 成人亚洲精品一区在线观看| 最近最新中文字幕大全电影3 | 最近在线观看免费完整版| 不卡av一区二区三区| 亚洲国产欧洲综合997久久, | 亚洲av美国av| 国产成人啪精品午夜网站| 天天躁夜夜躁狠狠躁躁| 禁无遮挡网站| 亚洲色图av天堂| 亚洲中文av在线| 日韩欧美在线二视频| 男人舔女人的私密视频| 亚洲成av人片免费观看| 国语自产精品视频在线第100页| av有码第一页| 精品一区二区三区av网在线观看| 欧美黑人精品巨大| 狂野欧美激情性xxxx| 亚洲一区二区三区色噜噜| 俄罗斯特黄特色一大片| 中文字幕av电影在线播放| 久热这里只有精品99| x7x7x7水蜜桃| 久久久久久人人人人人| 禁无遮挡网站| 少妇熟女aⅴ在线视频| 亚洲精品美女久久久久99蜜臀| 少妇的丰满在线观看| 国产成人欧美| 大香蕉久久成人网| 精品福利观看| 18禁美女被吸乳视频| 国产激情欧美一区二区| 88av欧美| 正在播放国产对白刺激| 久久精品夜夜夜夜夜久久蜜豆 | 美女扒开内裤让男人捅视频| 99国产极品粉嫩在线观看| 大香蕉久久成人网| 国产黄片美女视频| 久久精品国产亚洲av高清一级| 久久 成人 亚洲| 久久婷婷人人爽人人干人人爱| 国产视频内射| 久久精品人妻少妇| 亚洲一卡2卡3卡4卡5卡精品中文| 日本黄色视频三级网站网址| 久久精品91无色码中文字幕| 亚洲国产精品合色在线| 欧美日韩乱码在线| 国产熟女xx| 十八禁网站免费在线| 日韩欧美在线二视频| 亚洲真实伦在线观看| 亚洲色图 男人天堂 中文字幕| 白带黄色成豆腐渣| 亚洲七黄色美女视频| 一级毛片高清免费大全| 男女床上黄色一级片免费看| 最新美女视频免费是黄的| 美女午夜性视频免费| 99久久国产精品久久久| 国产午夜福利久久久久久| cao死你这个sao货| www日本在线高清视频| 国产亚洲欧美在线一区二区| 日韩欧美国产在线观看| 午夜福利一区二区在线看| 国产精品亚洲av一区麻豆| 国产精品综合久久久久久久免费| 真人一进一出gif抽搐免费| 国产精品免费视频内射| 亚洲激情在线av| 精品欧美一区二区三区在线| 村上凉子中文字幕在线| 中文字幕精品亚洲无线码一区 | 日韩视频一区二区在线观看| 国产成人欧美在线观看| 麻豆一二三区av精品| 给我免费播放毛片高清在线观看| 国产成人av激情在线播放| 18禁观看日本| av在线天堂中文字幕| 可以免费在线观看a视频的电影网站| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av香蕉五月| 俺也久久电影网| 国产激情欧美一区二区| 黑丝袜美女国产一区| 最新在线观看一区二区三区| 曰老女人黄片| 国产成年人精品一区二区| 午夜激情福利司机影院| 欧美国产精品va在线观看不卡| 草草在线视频免费看| 国产色视频综合| 亚洲国产欧美网| 国产精品久久久人人做人人爽| 一本一本综合久久| 日日摸夜夜添夜夜添小说| 999久久久国产精品视频| or卡值多少钱| 99热6这里只有精品| 亚洲自偷自拍图片 自拍| 久久99热这里只有精品18| 国产熟女午夜一区二区三区| 国产熟女xx| 免费av毛片视频| 18禁黄网站禁片午夜丰满| 757午夜福利合集在线观看| 在线免费观看的www视频| 99久久国产精品久久久| 国产精品国产高清国产av| 亚洲精品色激情综合| 97超级碰碰碰精品色视频在线观看| 婷婷亚洲欧美| 亚洲av成人av| 国产单亲对白刺激| 亚洲精品在线美女| 99国产综合亚洲精品| 国产精品 欧美亚洲| 少妇粗大呻吟视频| 90打野战视频偷拍视频| 熟女少妇亚洲综合色aaa.| 精品少妇一区二区三区视频日本电影| 成人亚洲精品一区在线观看| 成人三级做爰电影| 国产精品1区2区在线观看.| 十八禁网站免费在线| 99热6这里只有精品| 天天一区二区日本电影三级| 国产精品,欧美在线| 好男人在线观看高清免费视频 | av在线播放免费不卡| 久久天躁狠狠躁夜夜2o2o| 亚洲天堂国产精品一区在线| 丝袜在线中文字幕| 国产精品日韩av在线免费观看| www国产在线视频色| 久久国产精品影院| 一二三四社区在线视频社区8| 欧美色欧美亚洲另类二区| 不卡av一区二区三区| 香蕉av资源在线| 午夜老司机福利片| 好看av亚洲va欧美ⅴa在| 精华霜和精华液先用哪个| 麻豆av在线久日| 欧美日本视频| 91麻豆精品激情在线观看国产| 亚洲国产精品久久男人天堂| 香蕉av资源在线| 欧美色视频一区免费| 日韩中文字幕欧美一区二区| 国产亚洲av嫩草精品影院| 欧美成人午夜精品| 大香蕉久久成人网| 999精品在线视频| 精品欧美国产一区二区三| 精品国产超薄肉色丝袜足j| 黄片大片在线免费观看| 黄色 视频免费看| 成年人黄色毛片网站| 精品国产超薄肉色丝袜足j| 国产成人影院久久av| 欧美精品亚洲一区二区| 日韩大码丰满熟妇| 欧美日韩黄片免| 精品欧美国产一区二区三| 999精品在线视频| 成年人黄色毛片网站| 亚洲一区二区三区色噜噜| 真人一进一出gif抽搐免费| 欧美亚洲日本最大视频资源| 别揉我奶头~嗯~啊~动态视频| 日韩欧美 国产精品| 精品熟女少妇八av免费久了| 日韩欧美国产一区二区入口| 两人在一起打扑克的视频| 2021天堂中文幕一二区在线观 | 久热这里只有精品99| 日本在线视频免费播放| 亚洲男人天堂网一区| 黄色a级毛片大全视频| 国产成年人精品一区二区| 我的亚洲天堂| 亚洲精品国产精品久久久不卡| 免费一级毛片在线播放高清视频| 国产精品精品国产色婷婷| 国产熟女xx| 免费观看人在逋| 岛国在线观看网站| 中国美女看黄片| 一边摸一边做爽爽视频免费| 久久久久九九精品影院| 亚洲一区二区三区色噜噜| 国产熟女午夜一区二区三区| 亚洲专区中文字幕在线| 国产97色在线日韩免费| 亚洲国产精品999在线| 国内精品久久久久久久电影| 一区二区三区激情视频| 国产伦人伦偷精品视频| 亚洲欧洲精品一区二区精品久久久| 国产爱豆传媒在线观看 | 青草久久国产| 老司机在亚洲福利影院| 母亲3免费完整高清在线观看| 国产成人啪精品午夜网站| 久久久久亚洲av毛片大全| 亚洲国产精品成人综合色| 18禁裸乳无遮挡免费网站照片 | 欧美性长视频在线观看| 欧美黑人欧美精品刺激| 亚洲国产精品久久男人天堂| 波多野结衣av一区二区av| 非洲黑人性xxxx精品又粗又长| 不卡一级毛片| 香蕉国产在线看| 精品国产亚洲在线| 日韩av在线大香蕉| 精品久久久久久久人妻蜜臀av| 亚洲成人国产一区在线观看| 亚洲一区高清亚洲精品| 亚洲精品国产区一区二| 最近最新免费中文字幕在线| 中文字幕高清在线视频| 亚洲第一青青草原| 中文字幕人妻熟女乱码| 国产日本99.免费观看| 亚洲一区二区三区色噜噜| 国产高清有码在线观看视频 | 欧美日韩亚洲综合一区二区三区_| 国产aⅴ精品一区二区三区波| 久久国产乱子伦精品免费另类| 成年人黄色毛片网站| 又黄又爽又免费观看的视频| 久久久国产精品麻豆| 两个人看的免费小视频| 在线观看一区二区三区| xxxwww97欧美| 动漫黄色视频在线观看| 18禁美女被吸乳视频| 欧美最黄视频在线播放免费| 一本精品99久久精品77| 天天躁夜夜躁狠狠躁躁| 免费看美女性在线毛片视频| 欧美丝袜亚洲另类 | 国产精品 国内视频| 精品高清国产在线一区| 少妇被粗大的猛进出69影院| 国产单亲对白刺激| 精品人妻1区二区| 99热6这里只有精品| 成人国产一区最新在线观看| 中文字幕人妻熟女乱码| 丝袜在线中文字幕| 成人特级黄色片久久久久久久| 国产精品av久久久久免费| 免费一级毛片在线播放高清视频| 亚洲av成人av| 久久午夜综合久久蜜桃| 日韩精品免费视频一区二区三区| 最新在线观看一区二区三区| 一进一出抽搐gif免费好疼| 亚洲黑人精品在线| 久久精品国产综合久久久| 久久99热这里只有精品18| 国产单亲对白刺激| 欧美又色又爽又黄视频| 俺也久久电影网| 亚洲精品一卡2卡三卡4卡5卡| 香蕉丝袜av| 成人精品一区二区免费| АⅤ资源中文在线天堂| 高清在线国产一区| 日韩 欧美 亚洲 中文字幕| 国产精品亚洲av一区麻豆| 自线自在国产av| 久9热在线精品视频| 757午夜福利合集在线观看| 国产精品精品国产色婷婷| 国产精品野战在线观看| 亚洲av成人av| 观看免费一级毛片| 国产一卡二卡三卡精品| 久久伊人香网站| av在线播放免费不卡| 特大巨黑吊av在线直播 | 久久九九热精品免费| 亚洲专区中文字幕在线| 午夜福利成人在线免费观看| 国产av不卡久久| 久久精品亚洲精品国产色婷小说| a在线观看视频网站| 久久欧美精品欧美久久欧美| 国产精品 欧美亚洲| 国产一区二区激情短视频| 在线十欧美十亚洲十日本专区| 这个男人来自地球电影免费观看| 亚洲av成人不卡在线观看播放网| 欧洲精品卡2卡3卡4卡5卡区| 国产三级黄色录像| 欧美日韩精品网址| 日本一区二区免费在线视频| 午夜福利欧美成人| 看免费av毛片| 久久精品国产亚洲av高清一级| 久久精品91蜜桃| 97超级碰碰碰精品色视频在线观看| 国产免费av片在线观看野外av| 亚洲av中文字字幕乱码综合 | 中文亚洲av片在线观看爽| 黄色片一级片一级黄色片| av在线播放免费不卡| 欧美最黄视频在线播放免费| 亚洲av熟女| 亚洲国产欧美日韩在线播放| 中出人妻视频一区二区| 亚洲精品色激情综合| 国产伦一二天堂av在线观看| 国产亚洲av嫩草精品影院| 婷婷精品国产亚洲av在线| www.熟女人妻精品国产| 嫩草影院精品99| 侵犯人妻中文字幕一二三四区| 精品一区二区三区四区五区乱码| 色综合亚洲欧美另类图片| 精品卡一卡二卡四卡免费| 精品久久久久久久末码| 人人妻,人人澡人人爽秒播| 国内久久婷婷六月综合欲色啪| 欧美日韩中文字幕国产精品一区二区三区| 老熟妇乱子伦视频在线观看| 人人妻人人澡人人看| 亚洲片人在线观看| 精品高清国产在线一区| 这个男人来自地球电影免费观看| 女人高潮潮喷娇喘18禁视频| 观看免费一级毛片| 国产亚洲精品久久久久5区| 不卡一级毛片| 老司机福利观看| 天堂√8在线中文| 国产精品乱码一区二三区的特点| 久久香蕉激情| 国产主播在线观看一区二区| 黄色毛片三级朝国网站| 久久天躁狠狠躁夜夜2o2o| 久久久久亚洲av毛片大全| 91大片在线观看| 欧美日韩乱码在线| 成人av一区二区三区在线看| 亚洲男人的天堂狠狠| 90打野战视频偷拍视频| 午夜两性在线视频| 久久精品国产99精品国产亚洲性色| 国产精品一区二区三区四区久久 | 国产一区二区三区在线臀色熟女| 亚洲熟女毛片儿| 久久精品国产99精品国产亚洲性色| 天堂动漫精品| 亚洲午夜精品一区,二区,三区| 国产国语露脸激情在线看| 国产精品久久久人人做人人爽| 欧美日本视频| 国产色视频综合| 亚洲国产看品久久| 亚洲中文日韩欧美视频| a在线观看视频网站| 最好的美女福利视频网| 黄片播放在线免费| 69av精品久久久久久| 伊人久久大香线蕉亚洲五| 一本精品99久久精品77| 国产亚洲精品久久久久5区| 身体一侧抽搐| 欧美中文综合在线视频| 欧美乱码精品一区二区三区| 一个人观看的视频www高清免费观看 | 性色av乱码一区二区三区2| 欧美日韩黄片免| www.精华液| 欧美性猛交黑人性爽| 中文亚洲av片在线观看爽| 久久草成人影院| 天堂√8在线中文| 亚洲精品一卡2卡三卡4卡5卡| 两人在一起打扑克的视频| 国产欧美日韩精品亚洲av| 777久久人妻少妇嫩草av网站| 亚洲成a人片在线一区二区| 亚洲熟妇中文字幕五十中出| 国语自产精品视频在线第100页| 精品无人区乱码1区二区| 国产成年人精品一区二区| 国产精品国产高清国产av| 精品无人区乱码1区二区| 日韩欧美三级三区| 激情在线观看视频在线高清| 非洲黑人性xxxx精品又粗又长| 亚洲精品色激情综合| 白带黄色成豆腐渣| 久久久久久亚洲精品国产蜜桃av| 亚洲精品美女久久av网站| 午夜激情av网站| 久久久国产成人精品二区| 欧美激情极品国产一区二区三区| 在线观看免费视频日本深夜| 久久久久国产一级毛片高清牌|