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

    三體船阻力預(yù)報(bào)與側(cè)體布局的優(yōu)化?

    2016-06-16 02:54:46陳飛宇余建星馬維林劉見德趙羿羽

    陳飛宇, 余建星, 馬維林, 劉見德, 趙羿羽

    (天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室, 天津 300072)

    三體船阻力預(yù)報(bào)與側(cè)體布局的優(yōu)化?

    陳飛宇, 余建星, 馬維林, 劉見德??, 趙羿羽

    (天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室, 天津 300072)

    摘要:三體船具有一系列突出的優(yōu)點(diǎn),它優(yōu)良的快速性、耐波性和穩(wěn)性等特點(diǎn)為其應(yīng)用提供了廣闊的前景。側(cè)體布局的設(shè)計(jì)是優(yōu)化其性能的關(guān)鍵,為了探究其對(duì)阻力性能的影響,本文首先在GAMBIT中創(chuàng)建三體船數(shù)值模型,再采用FLUENT軟件對(duì)其進(jìn)行數(shù)值計(jì)算。與船模試驗(yàn)相比,數(shù)值計(jì)算方法的相對(duì)誤差小于2%,同時(shí)能準(zhǔn)確高效地捕捉到自由液面波形等流場(chǎng)數(shù)據(jù),為三體船阻力預(yù)報(bào)提供了一種高效可行的方法。在此基礎(chǔ)上改變側(cè)體縱橫偏距,計(jì)算并分析了25種側(cè)體布局方案的興波干擾情況和阻力性能,得到了側(cè)體布置在主體興波波峰中時(shí)阻力較大,布置在主體興波波谷中時(shí)阻力較小的結(jié)論。最后根據(jù)阻力最優(yōu)原則,確定了側(cè)體的最佳布置區(qū)域。不僅揭示了側(cè)體布局對(duì)三體船阻力的影響規(guī)律,還為三體實(shí)船的側(cè)體布局設(shè)計(jì)提供了新的思路。

    關(guān)鍵詞:三體船; CFD; 阻力預(yù)報(bào); 側(cè)體布局優(yōu)化

    引用格式:陳飛宇, 余建星, 馬維林, 等. 三體船阻力預(yù)報(bào)與側(cè)體布局的優(yōu)化[J].中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),2016, 46(5): 125-131.

    CHEN Fei-Yu, YU Jian-Xing, LIU Jian-De, et al. Prediction on resistance of trimaran and side hull optimization[J].Periodical of Ocean University of China, 2016, 46(5): 125-131.

    21世紀(jì)是海洋的世紀(jì),海洋考察、開發(fā)利用、海運(yùn)、海防等都對(duì)船舶設(shè)計(jì)提出了更高要求。高速三體船作為一種嶄新的船型,受到了廣泛關(guān)注。

    國(guó)外對(duì)三體船的研究已經(jīng)持續(xù)了三十幾年,軍用及民用實(shí)船日益增多,各方面研究成果領(lǐng)先于中國(guó)。Suzuki和Larsson等使用三維Rankine面元法對(duì)高速三體船的興波阻力進(jìn)行了數(shù)值計(jì)算[1-2]。Brizzolara等對(duì)圓舭型三體船和主體為折角線型的三體船進(jìn)行了船模試驗(yàn)和數(shù)值計(jì)算,對(duì)三體船的剩余阻力進(jìn)行了研究[3]。Javanmardi等通過自行編制的CFD程序計(jì)算并分析了不同側(cè)體位置對(duì)三體船阻力性能和操縱性能的影響[4]。

    中國(guó)對(duì)三體船的研究還處于起步階段。閆宏生等人提出了一種可調(diào)長(zhǎng)度高速三體船的概念,并初步分析了不同側(cè)體布局下的興波干擾情況和阻力性能[5]。酈云等對(duì)主體和側(cè)體均為Wigley船型的高速三體船進(jìn)行了模型試驗(yàn)研究[6]。陳京普等采用改進(jìn)的Dawson方法,開發(fā)了雙/三體船興波阻力數(shù)值預(yù)報(bào)程序,并考慮了方尾的處理[7]。蔡新功等利用Michell薄船理論結(jié)合模型試驗(yàn)對(duì)三體船的側(cè)體布局規(guī)律進(jìn)行了研究[8]。韓開佳等對(duì)新細(xì)長(zhǎng)體理論進(jìn)行改進(jìn),使之適用于卵型船等多體船型的興波阻力計(jì)算[9]。

    目前的阻力研究方法中,船模試驗(yàn)結(jié)果最為可靠,但船模制作耗時(shí)且費(fèi)用較高,改變型線進(jìn)行船型優(yōu)化效率低;理論計(jì)算方法則往往要求對(duì)計(jì)算對(duì)象進(jìn)行抽對(duì)船模試驗(yàn)三體船型進(jìn)行阻力數(shù)值模擬,將模擬結(jié)果將與船模試驗(yàn)結(jié)果進(jìn)行對(duì)比,驗(yàn)證FLUENT軟件在三體船黏性流場(chǎng)數(shù)值模擬方面的可靠性;在此基礎(chǔ)上設(shè)計(jì)25種側(cè)體布置方案,充分研究側(cè)體布局對(duì)三體船阻力性能的影響,為側(cè)體布局優(yōu)化提供指導(dǎo)。

    1基礎(chǔ)理論

    流體流動(dòng)要受物理守恒定律的支配,控制方程就是這些守恒定律的數(shù)學(xué)描述。對(duì)于船舶阻力計(jì)算,由于不考慮熱傳導(dǎo)和熱交換,因此能量守恒方程自動(dòng)成立。下面給出質(zhì)量守恒定律和動(dòng)量守恒定律對(duì)應(yīng)的控制方程:

    (1)質(zhì)量守恒方程又稱為連續(xù)性方程,任何流動(dòng)問題都必須滿足質(zhì)量守恒定律,該定律可以表述為:?jiǎn)挝粫r(shí)間內(nèi)流場(chǎng)中微元體質(zhì)量的增加,等于同一時(shí)間間隔內(nèi)流入該微元體的凈質(zhì)量,由此可導(dǎo)出流體流動(dòng)連續(xù)性方程的積分形式為:

    (1)

    (2)動(dòng)量守恒方程也稱作運(yùn)動(dòng)方程或N-S方程。象和簡(jiǎn)化,復(fù)雜形體的計(jì)算難以實(shí)現(xiàn)[10-12]。因此本文動(dòng)量守恒定律也是任何流動(dòng)系統(tǒng)都必須滿足的基本定律,其本質(zhì)是牛頓第二定律。該定律可表述為:微元體中流體的動(dòng)量對(duì)時(shí)間的變化率等于外界作用在該微元體上的各種力之和。按照這一定律,可以得出在慣性(非加速)坐標(biāo)系中的動(dòng)量守恒方程:

    (2)

    式中:p是靜壓;τij是應(yīng)力張量;ρgi,F(xiàn)i是重力體積力和其它體積力(如離散相相互作用產(chǎn)生的升力);Fi還可以包括其它模型源項(xiàng)或者用戶自定義源項(xiàng)。

    2計(jì)算模型與網(wǎng)格參數(shù)

    采用某課題的設(shè)計(jì)船型進(jìn)行計(jì)算,縮尺比取0.1,即采用船模尺度進(jìn)行計(jì)算,船?;緟?shù)如表1所示。其中主體為方尾船型,長(zhǎng)寬比為12.93;側(cè)體左右對(duì)稱,單個(gè)側(cè)體排水量占三體船總排水量的4.47%。

    表1 船?;緟?shù)

    將船模橫剖面型線導(dǎo)入Gambit創(chuàng)建船體曲面及計(jì)算域,創(chuàng)建的船體模型如圖1所示,圖中a表示側(cè)體船舯與主體船舯的縱向距離,p表示側(cè)體中縱剖面與主體中縱剖面的橫向距離。

    考慮到模型左右對(duì)稱,取一半模型進(jìn)行計(jì)算。計(jì)算域入口距主體艏部0.5倍船長(zhǎng),出口距主體艉部2.5倍船長(zhǎng);船寬方向取1倍船長(zhǎng)。計(jì)算域根據(jù)設(shè)計(jì)水線面分為空氣和水2部分,空氣層高度取4倍吃水,水深取6倍吃水。

    圖1 船體模型

    在Gambit中分塊劃分流場(chǎng),以便采用結(jié)構(gòu)網(wǎng)格離散空間計(jì)算域。為了準(zhǔn)確捕捉自由液面波形,在船體表面附近及設(shè)計(jì)水線面附近加密流場(chǎng)。網(wǎng)格劃分結(jié)果如圖2所示。

    圖2 流場(chǎng)網(wǎng)格

    邊界條件:空氣入口和水入口均設(shè)為速度入口;出口為自由出流;對(duì)稱面采用對(duì)稱邊界條件;船體表面采用無(wú)滑移壁面條件,其余采用滑移壁面邊界條件;采用VOF模擬自由液面;采用標(biāo)準(zhǔn)壁面函數(shù)計(jì)算摩擦力。

    3阻力預(yù)報(bào)方法驗(yàn)證

    采用上述三體船模型及參數(shù)對(duì)以下4個(gè)三體船側(cè)體布局方案進(jìn)行數(shù)值模擬,航速取4.88m/s,F(xiàn)r=0.677,相當(dāng)于實(shí)船航速30 kn。

    方案1:a=13.25m,p=7.38m;

    方案2:a=13.25m,p=9.02m;

    方案3:a=14.84m,p=7.38m;

    方案4:a=14.84m,p=9.02m。

    計(jì)算收斂后,打印阻力報(bào)告,得到各方案三體船的靜水阻力如表2所示。由于數(shù)值計(jì)算以主體中縱剖面為對(duì)稱面,只取一半模型進(jìn)行模擬,因此三體船總阻力為數(shù)值計(jì)算結(jié)果的2倍。表中還列出了船模試驗(yàn)相同航速下的阻力數(shù)值,同時(shí)計(jì)算出數(shù)值計(jì)算結(jié)果的相對(duì)誤差??梢钥吹綌?shù)值模擬結(jié)果的精度很高,相對(duì)誤差最大值僅為1.65%。

    表2 數(shù)值計(jì)算結(jié)果及船模試驗(yàn)結(jié)果

    圖3為方案3三體船自由液面波形,圖中Y表示自由液面高度。與圖4方案3船模試驗(yàn)波形對(duì)比可以看到數(shù)值計(jì)算得到的波形規(guī)則平整、逼真。我們知道三體船的阻力性能在很大程度上取決于側(cè)體的位置,合理布置主側(cè)體的相對(duì)位置,能夠充分獲取片體間的有利興波干擾,使興波阻力得到最大程度的減小。而興波阻力是由船體興波引起的壓力分布的改變所產(chǎn)生的,因此準(zhǔn)確模擬三體船自由液面波形是三體船阻力預(yù)報(bào)的關(guān)鍵點(diǎn)。

    圖3 方案3波形圖

    圖4 方案3船模試驗(yàn)波形圖

    根據(jù)上述阻力及波形對(duì)比結(jié)果,可以知道FLUENT軟件在三體船阻力數(shù)值預(yù)報(bào)中具有一定的可靠性,能夠較為準(zhǔn)確的計(jì)算出三體船阻力,并捕捉到三體船自由液面波形等流場(chǎng)數(shù)據(jù),可以作為三體船側(cè)體布局優(yōu)化的一種工具。

    4側(cè)體布局優(yōu)化研究

    下面對(duì)設(shè)計(jì)船型進(jìn)行側(cè)體布局優(yōu)化分析,以30 kn航速(Fr=0.677)下的總阻力為優(yōu)化目標(biāo)。為充分研究側(cè)體布局對(duì)三體船阻力性能的影響,分別設(shè)計(jì)了5個(gè)縱向偏距和5個(gè)橫向偏距,a/L=0~0.5,間隔為0.125,p/B=1~5,間隔為1,即a=0、6.625、13.25、19.875、26.5,p=4.1、8.2、12.3、16.4、20.5m,其中L和B分別為主體的設(shè)計(jì)水線長(zhǎng)和設(shè)計(jì)水線寬。每個(gè)橫向偏距對(duì)應(yīng)5種縱向布局,因此共有25種方案。

    為了更好地研究高速三體船主側(cè)體間的興波干擾規(guī)律,本文還對(duì)三體船主、側(cè)體單獨(dú)進(jìn)行數(shù)值模擬,以便和三體船進(jìn)行對(duì)比。模型的網(wǎng)格密度、參數(shù)設(shè)定等與之前的算例相同。

    計(jì)算收斂后,輸出各方案的自由液面波形圖,限于篇幅,下面給出部分波形圖(見圖5)。

    圖5 自由液面波形圖

    分析以上自由面波形圖,可以看出:

    (1) 三體船首興波波峰的高度及位置與主體波形沒有明顯差異,舯部以前的波形也基本相同。這是因?yàn)槿w船舯部以前的波形主要由主體產(chǎn)生,且側(cè)體布置在主體側(cè)后方,因此側(cè)體布局的變化對(duì)其影響不大。

    (2) 在p=4.1m的5個(gè)方案中,側(cè)體布局縱向變化對(duì)主體舯部以后的水線產(chǎn)生了明顯的影響;但當(dāng)p≥8.2m時(shí),側(cè)體布局變化對(duì)主體的水線影響不大,但仍會(huì)對(duì)主體側(cè)后方波形產(chǎn)生較大影響。

    (3) 當(dāng)側(cè)體處于主體首興波的波峰中時(shí),側(cè)體艏部?jī)?nèi)側(cè)水線急劇升高,甚至可以觀測(cè)到飛濺現(xiàn)象;當(dāng)側(cè)體處于主體側(cè)后方的波谷中時(shí),側(cè)體周圍的自由液面高度一般較小。

    (4) 總體上,側(cè)體內(nèi)側(cè)自由液面高度高于外側(cè);側(cè)體遠(yuǎn)離主體時(shí),側(cè)體整體水線較低。

    根據(jù)各方案計(jì)算得到的阻力數(shù)據(jù),繪制主側(cè)體興波阻力、黏性阻力、總阻力曲線及三體船總阻力等值線,如圖6~11所示。

    圖6 主體興波阻力曲線

    圖7 主體粘性阻力曲線

    分析上述阻力曲線,可以看出:

    (1)p=4.1m時(shí),受側(cè)體布局影響,主體阻力變化很大;p=8.2m時(shí),側(cè)體布局對(duì)主體的影響明顯減弱;p≥12.3m時(shí),主體阻力基本保持不變;這與圖5中的波形是一致的。對(duì)比圖6和7可以發(fā)現(xiàn),主體黏性阻力與興波阻力受側(cè)體布局影響的程度差異很大,主體黏性阻力與興波阻力最大值與最小值之差分別為2.4 N和33.2 N,因此側(cè)體布局對(duì)主體阻力的影響主要體現(xiàn)在興波阻力上。

    圖8 主體總阻力曲線

    圖9 側(cè)體興波阻力曲線

    圖10 側(cè)體粘性阻力曲線

    (2) 側(cè)體布局變化時(shí),側(cè)體黏性阻力與興波阻力均受到較大影響。總體上看,側(cè)體遠(yuǎn)離主體時(shí),主側(cè)體間興波干擾減弱,側(cè)體阻力變化幅度也相應(yīng)減小。將阻力曲線與圖5的波形圖相結(jié)合進(jìn)行分析,可以發(fā)現(xiàn)在側(cè)體橫向位置p相同的情況下,側(cè)體處于主體興波的波峰中時(shí),其阻力也處于峰值。

    圖11 側(cè)體總阻力曲線

    圖12 三體船總阻力等值線

    5結(jié)論

    根據(jù)上述研究,可以得到以下結(jié)論:

    (1) FLUENT軟件在三體船阻力數(shù)值預(yù)報(bào)中具有一定的可靠性,能夠較為準(zhǔn)確的計(jì)算出三體船阻力,并捕捉到三體船自由液面波形等流場(chǎng)數(shù)據(jù),可以作為三體船側(cè)體布局優(yōu)化的工具;

    (2) 側(cè)體布局對(duì)主體阻力的影響主要體現(xiàn)在興波阻力上,且僅當(dāng)側(cè)體橫向位置p/B<2時(shí)有較大影響;

    (3) 側(cè)體布局變化時(shí),側(cè)體黏性阻力與興波阻力均受到較大影響。側(cè)體縱向位置不變,橫向遠(yuǎn)離主體時(shí),主側(cè)體間興波干擾減弱,側(cè)體阻力變化幅度也相應(yīng)減??;側(cè)體橫向位置不變,則當(dāng)側(cè)體布置于主體興波的波峰中時(shí),側(cè)體阻力最大,布置于主體興波的波谷中時(shí),側(cè)體阻力最??;

    (4) 僅考慮阻力最優(yōu),則三體船側(cè)體的最佳布置位置為0.25≤a/L≤0.5,1≤p/B≤2的區(qū)域。

    參考文獻(xiàn):

    [1]Suzuki K, Nakata Y, Ikehata M. Numerical prediction on wave making resistance of high speed trimaran [C]. //Sydney: Proceedings of FAST’97: Conference, 1997: 611-621.

    [2]Larsson L, Janson C E, Brun P. A numerical investigation on trimaran configurations [C]. Sydney: Proceedings of FAST' 97 Conference, 1997: 537-544.

    [3]Brizzolara S, Bruzzone D. Hydrodynamic optimization of high-speed trimaran hull forms [C]. Sydney: Proceedings of the Jnternational Offshore and Polar Engineering Conference, 2008: 547-554.

    [4]Javanmardi M R, Jahanbakhsh E, Seif M S. Hydrodynamic analysis of trimaran vessels[J]. Polish Maritime Research, 2008, 15(55): 11-18.

    [5]閆宏生, 徐小穎, 張英晟, 等.可調(diào)長(zhǎng)度高速三體船的概念及水動(dòng)力研究[J].船舶工程, 2012, 34(5): 9-12.

    Yan Hong-Sheng,Xu Xiao-Ying,Zhang Ying-Sheng, et al. Concept and hydrodynamic analysis of high-speed adjustable-length trimaran[J]. Ship Engineering, 2012, 34(5): 9-12.

    [6]酈云, 盧曉平. 高速三體船阻力性能研究[J]. 船舶力學(xué), 2007, 11(2): 191-198.

    Li Yun, LU Xiao-Ping.An investigation on the resistance of high speed trimarans[J]. Journal of Ship Mechanics, 2007, 11(2):191-198.

    [7]陳京普, 朱德祥, 何術(shù)龍.雙體船/三體船興波阻力數(shù)值預(yù)報(bào)方法研究[J]. 船舶力學(xué), 2006, 10(2): 23-29.

    Chen Jing-Pu,Zhu De-Xiang,He Shu-Long.Research on numerical prediction method for wavemaking resistance of catamaran/trimaran[J].Journal of Ship Mechanics, 2006, 10(2): 23-29.

    [8]蔡新功, 王平, 謝小敏. 三體船方案優(yōu)化布局的阻力計(jì)算與試驗(yàn)研究[J]. 水動(dòng)力學(xué)研究與進(jìn)展, 2007, 22(2): 202-207.

    Cai Xin-Gong, Wang Ping, Xie Xiao-Min. Resistance study on alterative layouts of the trmiaran hull configuration[J]. Journal of Hydrodynamics 2007, 22(2): 202-207.

    [9]韓開佳, 黃德波. 三體船的興波阻力計(jì)算[J]. 哈爾濱工程大學(xué)學(xué)報(bào), 2000, 21(1): 6-10.

    Han Kai-Jia,Huang De-Bo. Wavemaking Resistance Calculation of Trimaran[J]. Journal of Harbin Engineering University, 2000, 21(1): 6-10.

    [10]盧曉平, 王中, 孫永華, 等. Rankine源Dawson型方法求解三體船興波阻力[J]. 華中科技大學(xué)學(xué)報(bào), 2008, 36(11): 103-107.

    Lu Xiao-Ping,Wang Zhong,Sun Yong-Hua, et al. Solution to wave resistance to trimanan using-Rankine source method of Dawson type[J]. Journal of Huazhong University of Science and Technology(Nature Science Edition), 2008, 36(11): 103-107.

    [11]謝玲玲, 徐海祥.深淺水中高速三體船興波阻力數(shù)值計(jì)算[J]. 華中科技大學(xué)學(xué)報(bào), 2009, 37(4): 130-132.

    Xie Ling-Ling,Xu Hai-Xiang. Numerical calculation of wave-making resistance of a high-speed trimaran in deep or shallow water[J]. Journal of Huazhong University of Science and Technology(Nature Science Edition), 2009, 37(4): 130-132.

    [12]劉昌明, 李云波, 李欲龍.三體船興阻力計(jì)算方法比較及興波干擾研究[J].中國(guó)艦船研究, 2009, 4(5): 8-32.

    Liu Chang-Ming,Li Yun-Bo, Li Yu-Long. The comparison of calculation methods to wave-making resistance and wave interference in trimaran[J]. Chinese Journal of Ship Research, 2009, 4(5): 8-32.

    責(zé)任編輯陳呈超

    Prediction on Resistance of Trimaran and Side Hull Optimization

    CHEN Fei-Yu, YU Jian-Xing, MA Wei-Lin, LIU Jian-De, ZHAO Yi-Yu

    (State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China)

    Abstract:Trimaran has a series of outstanding merits, its excellent performance of speediness, sea-keeping and stability have made it a wide application. Although the trimaran, a new type of hull form with excellent performance, has received more and more attention, problems still exist in studying its hull form and the arrangement of side hulls is the key factor to optimize its performance. Here, viscous computational fluid dynamics (CFD) technology was applied to research resistance performance of a trimaran. In order to explore the impact on the performance of resistance, the numerical trimaran model was designed first in GAMBIT, then conduct numerical calculation in FLUENT software, while taking free surfaces and flow viscosity into consideration. To validate the simulation results, conduct ship model test to calculate the resistances of the target trimaran with different layouts. And convert the experiment results to the practice ship by the two-dimensional and ‘1+k’ resistance calculation method, then compare the experiment results with the simulation results. Compared with model test results, the relative error of numerical resistance results is less than 2%, at the same time the information of the flow field was accurately captured, such as the waveform of the free surface, proves the reliability of FLUENT in predicting trimaran resistance. Then the wave-making interference and resistance performance for 25 different side hull configurations were computed and analyzed. The results show that the total resistance is higher when the side hulls are arranged in the peak of main hull wave and smaller when arranged in wave trough. In conclusion, the best area for side hull arrangement is signed out based on the principle of total resistance optimization. Not only reveals the influence rules of side hulls arrangement on trimaran resistances, but also offers a new method for the design of trimaran side hulls arrangement. At the same time, an apparent discrepancy in calculated and test results at high Froude numbers shows that the application or the code itself should be improved.

    Key words:Trimaran; CFD; resistance prediction; side hull optimization

    基金項(xiàng)目:? 國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(2014CB046803);國(guó)家自然科學(xué)基金創(chuàng)新研究群體科學(xué)基金項(xiàng)目(51021004);國(guó)家自然科學(xué)基金項(xiàng)目(51239008);海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室(上海交通大學(xué))開放課題項(xiàng)目(課題編號(hào)1201)資助

    收稿日期:2014-06-20;

    修訂日期:2015-01-15

    作者簡(jiǎn)介:陳飛宇(1987-),男,博士生,主要研究方向?yàn)榇傲黧w力學(xué)、結(jié)構(gòu)系統(tǒng)可靠性分析與優(yōu)化等。 ??通訊作者:E-mail:thisismrliu@163.com

    中圖法分類號(hào):U674.951

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1672-5174(2016)05-125-07

    DOI:10.16441/j.cnki.hdxb.20140236

    Supported by the Major State Basic Research Development Program of China,(2014CB046803);Project supported by the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (51021004);Project supported by the State Key Program of National Natural Science of China(51239008);Project supported by the State Key Laboratory of Ocean Engineering (Shanghai Jiao Tong University)

    听说在线观看完整版免费高清| x7x7x7水蜜桃| 国产精品永久免费网站| 特级一级黄色大片| 无遮挡黄片免费观看| 欧美一区二区国产精品久久精品| 这个男人来自地球电影免费观看| 欧美日韩亚洲国产一区二区在线观看| 日本熟妇午夜| 国产不卡一卡二| 国内毛片毛片毛片毛片毛片| 少妇的丰满在线观看| 麻豆一二三区av精品| 亚洲中文字幕一区二区三区有码在线看 | 欧美午夜高清在线| 日本 欧美在线| 91麻豆av在线| 国产乱人伦免费视频| 91av网一区二区| 亚洲性夜色夜夜综合| 久99久视频精品免费| 高潮久久久久久久久久久不卡| 国产精品98久久久久久宅男小说| 精品国产三级普通话版| 国产高清有码在线观看视频| 亚洲男人的天堂狠狠| 一个人看视频在线观看www免费 | 亚洲七黄色美女视频| 99久久国产精品久久久| 久久精品国产综合久久久| 男人舔女人的私密视频| 国产一区二区三区在线臀色熟女| 午夜激情福利司机影院| 人人妻人人看人人澡| 老熟妇仑乱视频hdxx| 中文字幕人妻丝袜一区二区| 别揉我奶头~嗯~啊~动态视频| 国产在线精品亚洲第一网站| 悠悠久久av| 国产精品一区二区三区四区久久| 久久九九热精品免费| 午夜久久久久精精品| h日本视频在线播放| 欧美在线黄色| 亚洲精华国产精华精| 亚洲avbb在线观看| 国产精品98久久久久久宅男小说| 他把我摸到了高潮在线观看| 亚洲成av人片免费观看| 可以在线观看毛片的网站| 亚洲成人久久性| 国产欧美日韩精品亚洲av| 桃红色精品国产亚洲av| 精品国产乱子伦一区二区三区| 一区二区三区激情视频| 美女黄网站色视频| 精品熟女少妇八av免费久了| 亚洲av成人一区二区三| 欧美av亚洲av综合av国产av| 亚洲黑人精品在线| 亚洲天堂国产精品一区在线| 亚洲 欧美 日韩 在线 免费| 久久精品国产清高在天天线| 禁无遮挡网站| 免费在线观看日本一区| 亚洲av成人不卡在线观看播放网| av在线天堂中文字幕| 女警被强在线播放| 国产精品99久久久久久久久| 久久亚洲精品不卡| 两性夫妻黄色片| 亚洲美女黄片视频| 亚洲片人在线观看| 国产蜜桃级精品一区二区三区| 日本精品一区二区三区蜜桃| 最近视频中文字幕2019在线8| 色播亚洲综合网| 99精品久久久久人妻精品| 成年女人毛片免费观看观看9| 精品久久久久久成人av| 亚洲第一欧美日韩一区二区三区| 欧美色视频一区免费| 91久久精品国产一区二区成人 | 久久久久久大精品| 午夜亚洲福利在线播放| 国产成人精品久久二区二区91| 午夜福利免费观看在线| 亚洲激情在线av| 无遮挡黄片免费观看| 日韩免费av在线播放| 美女免费视频网站| 黄片小视频在线播放| 日韩国内少妇激情av| 午夜激情欧美在线| 久久久久久久久久黄片| 又黄又粗又硬又大视频| 国产真实乱freesex| 精品国产三级普通话版| 亚洲一区二区三区不卡视频| 欧美色视频一区免费| 国产高清激情床上av| 亚洲成av人片在线播放无| 午夜福利视频1000在线观看| 巨乳人妻的诱惑在线观看| 99热精品在线国产| 最近最新免费中文字幕在线| 女人高潮潮喷娇喘18禁视频| 国产综合懂色| 欧美zozozo另类| 人人妻,人人澡人人爽秒播| 黑人欧美特级aaaaaa片| 国产精品久久久久久人妻精品电影| 999久久久国产精品视频| 午夜久久久久精精品| 久久久久久久精品吃奶| 亚洲精品一区av在线观看| 国内精品久久久久精免费| 老鸭窝网址在线观看| 在线视频色国产色| 小蜜桃在线观看免费完整版高清| 国产成人影院久久av| 欧美在线黄色| 91九色精品人成在线观看| 欧美日韩一级在线毛片| 亚洲激情在线av| 成年版毛片免费区| 国产97色在线日韩免费| 99在线人妻在线中文字幕| 国产视频一区二区在线看| 精品午夜福利视频在线观看一区| 国产一区二区激情短视频| 99视频精品全部免费 在线 | 三级国产精品欧美在线观看 | 99国产极品粉嫩在线观看| 91老司机精品| 欧美成人一区二区免费高清观看 | 后天国语完整版免费观看| 一进一出好大好爽视频| 色老头精品视频在线观看| 色视频www国产| 久久久久久大精品| 日韩免费av在线播放| 亚洲 欧美一区二区三区| 欧美成狂野欧美在线观看| 九九久久精品国产亚洲av麻豆 | 国产精品久久久久久久电影 | 在线观看免费午夜福利视频| 99热只有精品国产| 全区人妻精品视频| 国产精品亚洲一级av第二区| 久久久久久久久中文| 99在线视频只有这里精品首页| 免费电影在线观看免费观看| 网址你懂的国产日韩在线| 欧美中文综合在线视频| 日韩欧美免费精品| 日韩成人在线观看一区二区三区| 9191精品国产免费久久| 午夜激情福利司机影院| 欧美日韩中文字幕国产精品一区二区三区| 国产精品亚洲一级av第二区| 丝袜人妻中文字幕| 嫩草影视91久久| 亚洲第一电影网av| 欧美xxxx黑人xx丫x性爽| 精品国产美女av久久久久小说| 欧美日韩精品网址| 婷婷丁香在线五月| 两个人视频免费观看高清| 国产一区二区在线观看日韩 | 在线观看午夜福利视频| 国产精品久久久久久精品电影| av欧美777| 一个人免费在线观看的高清视频| 90打野战视频偷拍视频| 中文字幕熟女人妻在线| 偷拍熟女少妇极品色| 免费大片18禁| 亚洲人成电影免费在线| 亚洲七黄色美女视频| 久久天堂一区二区三区四区| 舔av片在线| 麻豆国产av国片精品| 国产精品女同一区二区软件 | 欧美日韩国产亚洲二区| 久久久久国产一级毛片高清牌| 国产精品久久久久久人妻精品电影| 欧美午夜高清在线| 12—13女人毛片做爰片一| 色综合婷婷激情| 婷婷精品国产亚洲av| 午夜福利在线在线| 人人妻人人澡欧美一区二区| 亚洲国产精品sss在线观看| 美女午夜性视频免费| 午夜a级毛片| 91麻豆精品激情在线观看国产| 免费高清视频大片| 一级a爱片免费观看的视频| 国产亚洲精品久久久久久毛片| 脱女人内裤的视频| 国产精品亚洲av一区麻豆| 免费搜索国产男女视频| 亚洲av成人av| 2021天堂中文幕一二区在线观| 手机成人av网站| 国产在线精品亚洲第一网站| 99热6这里只有精品| cao死你这个sao货| 国产高清视频在线观看网站| 国产精品久久久久久亚洲av鲁大| 国产精品99久久久久久久久| 国产人伦9x9x在线观看| 精品熟女少妇八av免费久了| 在线观看美女被高潮喷水网站 | 精品乱码久久久久久99久播| 两个人的视频大全免费| 搞女人的毛片| 最近最新中文字幕大全电影3| 三级毛片av免费| cao死你这个sao货| АⅤ资源中文在线天堂| 日韩欧美精品v在线| 亚洲人成网站在线播放欧美日韩| 国产伦精品一区二区三区四那| 99久久综合精品五月天人人| 日本一二三区视频观看| 最近视频中文字幕2019在线8| 草草在线视频免费看| 99久久久亚洲精品蜜臀av| 超碰成人久久| 精品一区二区三区av网在线观看| 最新在线观看一区二区三区| 亚洲欧美激情综合另类| 18美女黄网站色大片免费观看| 国产精品九九99| av国产免费在线观看| 999精品在线视频| 天天一区二区日本电影三级| 午夜精品一区二区三区免费看| 亚洲成a人片在线一区二区| www.www免费av| 免费看美女性在线毛片视频| 国产精品香港三级国产av潘金莲| 一区二区三区高清视频在线| 天天一区二区日本电影三级| 免费在线观看日本一区| 熟女少妇亚洲综合色aaa.| 在线观看免费午夜福利视频| 免费人成视频x8x8入口观看| 国产成人福利小说| 国内揄拍国产精品人妻在线| 91字幕亚洲| 99视频精品全部免费 在线 | 丰满人妻一区二区三区视频av | 午夜福利成人在线免费观看| 麻豆成人午夜福利视频| 亚洲国产欧美一区二区综合| 成在线人永久免费视频| 男女床上黄色一级片免费看| 色噜噜av男人的天堂激情| 此物有八面人人有两片| 欧美日韩综合久久久久久 | 99国产综合亚洲精品| 美女高潮喷水抽搐中文字幕| 色噜噜av男人的天堂激情| 波多野结衣高清无吗| 日韩中文字幕欧美一区二区| 又黄又粗又硬又大视频| 99久久精品一区二区三区| 国产aⅴ精品一区二区三区波| 少妇的丰满在线观看| 国产精品美女特级片免费视频播放器 | 性色avwww在线观看| 最近最新中文字幕大全免费视频| 国产精品日韩av在线免费观看| 他把我摸到了高潮在线观看| 少妇人妻一区二区三区视频| 久久久久久九九精品二区国产| 国产精品一区二区精品视频观看| 国产伦在线观看视频一区| 99久久精品一区二区三区| 黑人欧美特级aaaaaa片| 日韩人妻高清精品专区| or卡值多少钱| 老司机午夜福利在线观看视频| 国产精品美女特级片免费视频播放器 | 最近在线观看免费完整版| 老鸭窝网址在线观看| 香蕉丝袜av| 午夜免费激情av| 可以在线观看毛片的网站| 国语自产精品视频在线第100页| 欧美+亚洲+日韩+国产| 国产精品香港三级国产av潘金莲| 九九热线精品视视频播放| 国产精华一区二区三区| 美女cb高潮喷水在线观看 | 欧美乱码精品一区二区三区| 一区二区三区国产精品乱码| 国产精品野战在线观看| 美女被艹到高潮喷水动态| 国产探花在线观看一区二区| 国产又黄又爽又无遮挡在线| 国产高清videossex| 日本在线视频免费播放| 成熟少妇高潮喷水视频| 欧美成人免费av一区二区三区| 中出人妻视频一区二区| 久久久国产成人精品二区| 国产免费男女视频| 嫩草影视91久久| 十八禁网站免费在线| а√天堂www在线а√下载| 在线观看日韩欧美| 久久久久精品国产欧美久久久| 国产精品99久久久久久久久| 日本与韩国留学比较| 国产精品久久久av美女十八| 国产成人系列免费观看| 一级a爱片免费观看的视频| 午夜福利欧美成人| 国内揄拍国产精品人妻在线| 草草在线视频免费看| 亚洲欧洲精品一区二区精品久久久| 亚洲欧美激情综合另类| 国产97色在线日韩免费| 老司机午夜福利在线观看视频| 亚洲第一欧美日韩一区二区三区| 亚洲欧美精品综合久久99| 日韩欧美在线二视频| 天堂影院成人在线观看| 精品久久久久久久久久免费视频| 午夜精品一区二区三区免费看| 中文字幕久久专区| 少妇的逼水好多| 国产熟女xx| 久久欧美精品欧美久久欧美| 精品免费久久久久久久清纯| 免费在线观看成人毛片| 中文字幕熟女人妻在线| 国产一区二区在线av高清观看| 在线视频色国产色| 日本 av在线| 丰满人妻一区二区三区视频av | netflix在线观看网站| 国产伦在线观看视频一区| 欧美激情久久久久久爽电影| 一级黄色大片毛片| 国产一级毛片七仙女欲春2| 两个人的视频大全免费| 99国产极品粉嫩在线观看| 五月伊人婷婷丁香| 99久久精品热视频| www.999成人在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国产精品av久久久久免费| 91久久精品国产一区二区成人 | 久久久久久九九精品二区国产| 国产av在哪里看| 日本 av在线| 五月玫瑰六月丁香| 一个人免费在线观看的高清视频| 男女之事视频高清在线观看| 天天躁日日操中文字幕| 欧美乱妇无乱码| 美女 人体艺术 gogo| 欧美一级毛片孕妇| 最新美女视频免费是黄的| 一级毛片高清免费大全| 欧美绝顶高潮抽搐喷水| 在线a可以看的网站| 国产乱人视频| 97人妻精品一区二区三区麻豆| 国产蜜桃级精品一区二区三区| АⅤ资源中文在线天堂| 91av网站免费观看| 久久久久九九精品影院| x7x7x7水蜜桃| 久久久久精品国产欧美久久久| 可以在线观看毛片的网站| 超碰成人久久| 日本黄大片高清| 国产野战对白在线观看| 久久精品国产清高在天天线| 90打野战视频偷拍视频| 啪啪无遮挡十八禁网站| 国产欧美日韩一区二区精品| 少妇的丰满在线观看| 亚洲在线自拍视频| 无限看片的www在线观看| 中文亚洲av片在线观看爽| 国产三级中文精品| 亚洲国产精品999在线| 中文字幕人妻丝袜一区二区| 国产伦人伦偷精品视频| 桃红色精品国产亚洲av| 欧美丝袜亚洲另类 | 久久精品91无色码中文字幕| 在线观看免费午夜福利视频| 国产激情久久老熟女| 99国产精品99久久久久| 国语自产精品视频在线第100页| 亚洲 国产 在线| 香蕉av资源在线| 制服丝袜大香蕉在线| 免费观看的影片在线观看| 久久草成人影院| 看免费av毛片| 国产精品国产高清国产av| 1000部很黄的大片| 草草在线视频免费看| 久久久精品欧美日韩精品| 国产精华一区二区三区| 91麻豆av在线| 久久久久久久久中文| 亚洲中文字幕日韩| 啦啦啦免费观看视频1| 国产午夜精品久久久久久| 国产精品 国内视频| 两人在一起打扑克的视频| 欧美一级a爱片免费观看看| 午夜亚洲福利在线播放| 国内精品久久久久久久电影| 国产精品99久久久久久久久| 午夜a级毛片| 禁无遮挡网站| 国产1区2区3区精品| 国产一级毛片七仙女欲春2| 在线永久观看黄色视频| 18禁美女被吸乳视频| 免费在线观看成人毛片| 日日夜夜操网爽| 黑人巨大精品欧美一区二区mp4| 美女大奶头视频| 综合色av麻豆| 久9热在线精品视频| 亚洲欧洲精品一区二区精品久久久| 国产伦一二天堂av在线观看| 午夜视频精品福利| 色吧在线观看| 男女视频在线观看网站免费| 两人在一起打扑克的视频| 丁香六月欧美| av在线蜜桃| 51午夜福利影视在线观看| 成人18禁在线播放| 悠悠久久av| а√天堂www在线а√下载| 天堂影院成人在线观看| 国产乱人视频| 搡老妇女老女人老熟妇| 久久国产精品人妻蜜桃| 97碰自拍视频| 嫁个100分男人电影在线观看| 久久午夜综合久久蜜桃| 女同久久另类99精品国产91| 久久久成人免费电影| 日本 欧美在线| 法律面前人人平等表现在哪些方面| 成人高潮视频无遮挡免费网站| 久久国产乱子伦精品免费另类| or卡值多少钱| 亚洲人成伊人成综合网2020| 日日摸夜夜添夜夜添小说| 国产高潮美女av| 床上黄色一级片| 国产精品美女特级片免费视频播放器 | 国产成+人综合+亚洲专区| 成人特级av手机在线观看| 国产精品电影一区二区三区| 欧美黄色片欧美黄色片| 久久久久精品国产欧美久久久| 男女午夜视频在线观看| 欧美性猛交╳xxx乱大交人| x7x7x7水蜜桃| 国产亚洲精品一区二区www| 亚洲精品在线观看二区| 淫秽高清视频在线观看| 在线免费观看不下载黄p国产 | xxxwww97欧美| 国产亚洲精品一区二区www| 午夜福利18| 嫩草影院精品99| 欧美黄色淫秽网站| 欧美日韩瑟瑟在线播放| 国产一区在线观看成人免费| 亚洲国产精品成人综合色| 亚洲av电影不卡..在线观看| 亚洲av五月六月丁香网| 日韩欧美一区二区三区在线观看| 丝袜人妻中文字幕| 91老司机精品| 国内少妇人妻偷人精品xxx网站 | 国产成人啪精品午夜网站| 国产乱人伦免费视频| 久久精品国产清高在天天线| 禁无遮挡网站| 美女被艹到高潮喷水动态| 国产精品久久久久久亚洲av鲁大| 黄色视频,在线免费观看| 中出人妻视频一区二区| h日本视频在线播放| 丰满人妻一区二区三区视频av | 人人妻,人人澡人人爽秒播| 免费看美女性在线毛片视频| 精品午夜福利视频在线观看一区| 黄色成人免费大全| 欧美黑人欧美精品刺激| 熟女人妻精品中文字幕| 韩国av一区二区三区四区| 亚洲成人精品中文字幕电影| 亚洲国产精品成人综合色| 国产aⅴ精品一区二区三区波| 免费在线观看影片大全网站| 国产精品日韩av在线免费观看| 日本黄色片子视频| 精品国内亚洲2022精品成人| 久久香蕉精品热| 色吧在线观看| 国产黄色小视频在线观看| 日本在线视频免费播放| 999精品在线视频| a级毛片a级免费在线| 婷婷丁香在线五月| 一区二区三区激情视频| 九九热线精品视视频播放| 91麻豆av在线| 99视频精品全部免费 在线 | 美女扒开内裤让男人捅视频| 女同久久另类99精品国产91| 国产一区二区激情短视频| 老司机午夜福利在线观看视频| 老熟妇乱子伦视频在线观看| 在线免费观看的www视频| 一进一出抽搐动态| 日韩精品中文字幕看吧| 亚洲熟妇中文字幕五十中出| 99久久国产精品久久久| 欧美激情久久久久久爽电影| 亚洲av美国av| 国产69精品久久久久777片 | 亚洲国产欧洲综合997久久,| 桃色一区二区三区在线观看| 三级男女做爰猛烈吃奶摸视频| 国产午夜精品久久久久久| 国产v大片淫在线免费观看| 欧美日韩黄片免| 成人国产一区最新在线观看| av视频在线观看入口| 午夜精品一区二区三区免费看| 亚洲美女黄片视频| 男插女下体视频免费在线播放| 中文资源天堂在线| 一进一出好大好爽视频| 久久久久性生活片| 国产亚洲精品久久久久久毛片| 中文字幕熟女人妻在线| 在线播放国产精品三级| 午夜福利在线在线| 俺也久久电影网| 久久久久免费精品人妻一区二区| 国产精品一及| 日本黄色片子视频| 国产极品精品免费视频能看的| 国产黄色小视频在线观看| 免费在线观看日本一区| 在线观看一区二区三区| 欧美日韩一级在线毛片| 精品福利观看| 亚洲性夜色夜夜综合| 亚洲av成人精品一区久久| av国产免费在线观看| 熟女人妻精品中文字幕| 免费av不卡在线播放| 俄罗斯特黄特色一大片| 久久久精品大字幕| 一a级毛片在线观看| 18禁黄网站禁片午夜丰满| 小说图片视频综合网站| 黄频高清免费视频| 国产伦精品一区二区三区四那| 丰满的人妻完整版| 久久草成人影院| 日韩精品青青久久久久久| 嫁个100分男人电影在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产欧美人成| 中文字幕最新亚洲高清| 亚洲成人中文字幕在线播放| 在线观看免费午夜福利视频| 一个人看的www免费观看视频| 偷拍熟女少妇极品色| 国产伦精品一区二区三区四那| 91麻豆精品激情在线观看国产| 欧美成人免费av一区二区三区| 亚洲av成人不卡在线观看播放网| 免费观看的影片在线观看| h日本视频在线播放| 日本免费a在线| 亚洲国产欧美网| 麻豆国产97在线/欧美| 九色国产91popny在线| 国产伦精品一区二区三区视频9 | 久久亚洲精品不卡| 中出人妻视频一区二区| 国产精品久久久久久精品电影| 亚洲av中文字字幕乱码综合| 亚洲七黄色美女视频| 日本成人三级电影网站| 亚洲成人久久爱视频| 色噜噜av男人的天堂激情| 久久久久久久久免费视频了| 亚洲无线在线观看| 欧美精品啪啪一区二区三区|