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

    波浪滑翔器柔性水翼推進(jìn)性能數(shù)值仿真分析

    2021-09-10 08:25:50孫秀軍桑宏強(qiáng)
    關(guān)鍵詞:水翼滑翔波浪

    劉 芬,彭 彬,孫秀軍,桑宏強(qiáng),2

    波浪滑翔器柔性水翼推進(jìn)性能數(shù)值仿真分析

    劉 芬1,彭 彬1,孫秀軍3,4*,桑宏強(qiáng)1,2

    (1.天津工業(yè)大學(xué)機(jī)械工程學(xué)院,天津,300387; 2.天津市現(xiàn)代機(jī)電裝備技術(shù)重點(diǎn)實(shí)驗(yàn)室,天津,300387; 3.中國海洋大學(xué)物理海洋教育部重點(diǎn)實(shí)驗(yàn)室,山東 青島,266100; 4.青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室,山東 青島,266237)

    水翼是波浪滑翔器波浪動力轉(zhuǎn)換的關(guān)鍵部件,直接決定了波浪滑翔器的航行性能。文中研究了波浪滑翔器柔性水翼的推進(jìn)性能,根據(jù)波浪滑翔器工作時水翼的運(yùn)動特點(diǎn),基于COMSOL Multiphysics軟件仿真出柔性水翼主動升沉運(yùn)動下被動變形的過程,重點(diǎn)研究波浪參數(shù)、俯仰運(yùn)動參數(shù)和斯特哈勞爾數(shù)對推進(jìn)性能的影響。仿真結(jié)果表明,在一定范圍內(nèi),水翼推進(jìn)性能隨波高、斯特哈勞爾數(shù)的增大而增大,隨波周期的增大而減小; 當(dāng)波周期較大時,波高和波周期的變化對水翼推進(jìn)性能的影響較小; 當(dāng)俯仰周期與波周期相等時,俯仰運(yùn)動對水翼推進(jìn)性能有積極作用。文中結(jié)果可為波浪滑翔器被動柔性水翼后續(xù)研究提供參考。

    波浪滑翔器; 柔性水翼; 推進(jìn)性能; 數(shù)值仿真

    0 引言

    波浪滑翔器是一種新型海洋移動觀測平臺,可以將波浪能轉(zhuǎn)換為前向驅(qū)動動能,具有長期自主觀測能力,目前已廣泛應(yīng)用于水文氣象數(shù)據(jù)觀測、環(huán)境監(jiān)測、生物追蹤及遠(yuǎn)程預(yù)警等領(lǐng)域[1-3]。

    水翼作為決定波浪滑翔器波浪動力轉(zhuǎn)換的關(guān)鍵部件,國內(nèi)外學(xué)者已經(jīng)進(jìn)行了大量研究。桑宏強(qiáng)等[4]通過仿真分析提出了振蕩水翼的最優(yōu)運(yùn)動方法,研究了水翼運(yùn)動的極限位置等對波浪滑翔器推進(jìn)速度的影響; 呂博元等[5]應(yīng)用STAR-CCM+軟件對水翼的擺動進(jìn)行仿真,研究了被動擺動、主動擺動、限位角及波浪參數(shù)等對水翼推進(jìn)性能的影響,可知主動擺動的推進(jìn)性能優(yōu)于被動擺動,水翼推進(jìn)性能受波浪參數(shù)影響較大,且隨波高和波浪頻率的增加而增加,在波浪頻率較小時,波浪頻率的變化對推進(jìn)性能的影響較小; Liu等[6]對規(guī)則波作用下振蕩柔性水翼的推進(jìn)性能進(jìn)行了研究,可知在規(guī)則波波動頻率與水翼振蕩頻率相等的情況下,水翼的推進(jìn)性能最佳; Liu等[7]針對后緣柔性的主動振蕩水翼能量汲取性能進(jìn)行研究,發(fā)現(xiàn)水翼柔性激勵了前緣渦的產(chǎn)生,提升了能量汲取能力; Alben等[8]發(fā)現(xiàn)自主推進(jìn)柔性體的推進(jìn)速度存在共振峰值,提出理論模型,并通過了實(shí)驗(yàn)驗(yàn)證; Prempraneerach等[9]進(jìn)行了二維柔性擺動水翼實(shí)驗(yàn),發(fā)現(xiàn)合理的弦向變形使得水翼效率提高了36%; Katz等[10]通過勢流理論分析得出,當(dāng)柔性尾鰭具有弦向變形時,其推進(jìn)效率比剛性尾鰭高; 田寶強(qiáng)等[11]在傳統(tǒng)波浪滑翔器的基礎(chǔ)上,引入了柔性蹼翼,搭建了原理樣機(jī),通過水池實(shí)驗(yàn)驗(yàn)證了結(jié)構(gòu)的可行性; 周凱[12]研究了典型仿生水翼的推進(jìn)及能量吸收機(jī)制,結(jié)果表明對稱運(yùn)動模式下水翼的推進(jìn)性能隨斯特哈勞爾數(shù)的增大而增大,斯特哈勞爾數(shù)過大,會導(dǎo)致流場的不穩(wěn)定。

    文中在傳統(tǒng)波浪滑翔器的基礎(chǔ)上,引入柔性水翼進(jìn)行研究,以期進(jìn)一步提升波浪滑翔器的航行性能。根據(jù)波浪滑翔器及其水翼機(jī)構(gòu)特殊的運(yùn)動方式與特點(diǎn),利用COMSOL Multiphysics軟件中的流體模塊、多體動力學(xué)模塊和多物理場耦合模塊,仿真柔性水翼隨波浪滑翔器運(yùn)動過程中產(chǎn)生的升沉運(yùn)動、俯仰運(yùn)動及形變的耦合過程。研究運(yùn)動參數(shù)對柔性水翼推進(jìn)性能的影響,主要包括波浪參數(shù)、俯仰運(yùn)動參數(shù)以及斯特哈勞數(shù),可為柔性水翼的開發(fā)提供參考。

    1 柔性水翼物理模型

    1.1 工作原理

    為進(jìn)一步提升波浪滑翔器的航行性能,引入柔性水翼。柔性水翼在結(jié)構(gòu)上與傳統(tǒng)波浪滑翔器的水翼機(jī)構(gòu)有很大的差別,其主要通過在水動力下的彈性變形實(shí)現(xiàn)波浪能的吸收與利用。柔性水翼的優(yōu)點(diǎn)主要表現(xiàn)在: 較剛性水翼結(jié)構(gòu)更簡單,可靠性強(qiáng); 較剛性水翼更易輕質(zhì),同等水動力條件下具有較小的慣性; 可綜合利用波浪滑翔器的升沉、俯仰和橫滾運(yùn)動產(chǎn)生前向驅(qū)動力。

    根據(jù)以往的實(shí)驗(yàn)可知,波浪滑翔器在實(shí)際運(yùn)動中的橫滾較小,一般不超過±5°,因此文中不考慮波浪滑翔器橫滾運(yùn)動對柔性水翼的影響。假定波浪滑翔器以恒定的速度前向運(yùn)動,當(dāng)波峰來臨時,波浪滑翔器在波浪的作用下產(chǎn)生上升運(yùn)動與俯仰運(yùn)動,同時柔性水翼受水動力的作用沿中心線向下彎曲,與水流形成一定攻角,柔性水翼因此受到升力和阻力的共同作用而產(chǎn)生前向動力,進(jìn)而作用于波浪滑翔器,如圖1(a)所示; 當(dāng)波谷來臨時,波浪滑翔器在波浪的作用下產(chǎn)生下沉運(yùn)動與俯仰運(yùn)動,同時柔性水翼在水動力的作用下沿中心線向上彎曲,與水流形成一定的攻角,柔性水翼因此受到升力和阻力的共同作用產(chǎn)生前向運(yùn)動,進(jìn)而作用于波浪滑翔器,如圖1(b)所示。因此,柔性水翼在波浪中始終通過前向驅(qū)動力完成對滑翔機(jī)的推進(jìn)。

    1.2 運(yùn)動規(guī)律

    實(shí)際海況中,波浪形式非常復(fù)雜,導(dǎo)致水翼在波浪中的運(yùn)動也極為復(fù)雜,通過數(shù)值仿真極難實(shí)現(xiàn)。為方便后續(xù)研究,結(jié)合線性波理論與柔性水翼的工作原理,對運(yùn)動過程進(jìn)行簡化,賦予柔性水翼升沉運(yùn)動、俯仰運(yùn)動以及流場作用下的變形,將波浪的運(yùn)動、柔性水翼的升沉運(yùn)動與俯仰運(yùn)動均設(shè)定為正弦函數(shù),具體表達(dá)式為

    波浪滑翔器能夠充分利用波浪的起伏將水動力轉(zhuǎn)化為前向推力,在升沉方向具有良好的隨波性,因此認(rèn)為波浪運(yùn)動與柔性水翼升沉運(yùn)動變化一致; 柔性水翼由初始狀態(tài)下沉?xí)r,在水動力的作用下,俯仰運(yùn)動角由平衡位逐漸向負(fù)向增大; 翼前緣相對于波浪滑翔器是固定的,所以翼俯仰運(yùn)動中心相對位置也是固定的,但為了后續(xù)研究的典型性,文中簡化翼前緣與俯仰運(yùn)動中心一致。根據(jù)上述表述,可將式(1)中的參數(shù)設(shè)定為

    1.3 動力學(xué)參數(shù)

    對于柔性水翼而言,無量綱參數(shù)斯特哈勞爾數(shù)對推進(jìn)性能有很大的影響,且

    為更方便和直觀地衡量柔性水翼推進(jìn)性能,將其在流場中前向方向的受力轉(zhuǎn)化為無量綱形式[5],推力系數(shù)為

    為更直觀地對比柔性水翼推進(jìn)性能,計(jì)算得到1個周期內(nèi)的平均推力系數(shù)[5]為

    式中,為波周期。

    2 仿真模型建立

    2.1 計(jì)算模型

    柔性水翼采用NACA0006翼型剖面,如圖2所示。文中取翼弦長度為180 mm,翼后緣長度C為140 mm,翼前緣長度C為40 mm。

    圖2 柔性水翼模型

    翼后緣呈彈性,其柔度系數(shù)E[12]定義為

    式中: 取E為5.7×103;為翼后緣的楊氏模量。

    2.2 計(jì)算域及網(wǎng)格劃分

    采用動網(wǎng)格技術(shù)研究柔性水翼在流場中的水動力特性。如圖3所示,柔性水翼的計(jì)算域?yàn)?0×(+20),翼表面設(shè)置為無滑移邊界條件,其中心與計(jì)算域縱向中心線的距離為/2,與計(jì)算域橫向中心線重合,右側(cè)邊界設(shè)置為速度入口,左側(cè)邊界設(shè)置為壓力出口,上下邊界設(shè)置為壁面。為確保計(jì)算的準(zhǔn)確性,翼表面增設(shè)大小為5×5的變形域,其中心與翼中心重合,變形域頂點(diǎn)與整體計(jì)算域頂點(diǎn)連線作輔助邊界。

    采用湍流模型和COMSOL Multiphysics軟件對計(jì)算域進(jìn)行非結(jié)構(gòu)網(wǎng)格劃分,圖4為網(wǎng)格劃分示意圖。

    圖3 計(jì)算域及邊界條件示意圖

    圖4 計(jì)算域網(wǎng)格劃分

    2.3 模型驗(yàn)證

    為驗(yàn)證所建立模型的準(zhǔn)確性,對網(wǎng)格尺寸以及計(jì)算可靠性進(jìn)行驗(yàn)證。

    2.3.1 網(wǎng)格無關(guān)驗(yàn)證

    圖5 不同網(wǎng)格數(shù)量下翼推力系數(shù)變化曲線

    2.3.2 可靠性驗(yàn)證

    表1 計(jì)算結(jié)果對比

    3 仿真結(jié)果與分析

    3.1 波浪參數(shù)對翼推進(jìn)性能的影響

    重點(diǎn)分析波高和波周期對翼推進(jìn)性能的影響,其中波高為0.2~0.8 m,以0.2 m間隔選取,波周期為2~8 s,以2 s間隔選取,具體參數(shù)如表2所示。

    表2 數(shù)值仿真波浪參數(shù)

    3.1.1 波高

    圖6是波周期分別為=2,4 ,6,8 s時,不同波高下翼推力系數(shù)C()隨時間變化曲線。由于第1個周期處于啟動階段,流場處于不穩(wěn)定狀態(tài),而第2個周期穩(wěn)定程度相對較好,因此選擇第2個周期進(jìn)行分析。

    由圖6(a)可知,推力系數(shù)隨波高的增大而增大,且在波高=0.4,0.6,0.8 m時,推力系數(shù)隨時間的變化波動較大,沒有呈現(xiàn)出明顯的周期性,說明柔性水翼運(yùn)動中周圍的流場并未趨于穩(wěn)定; 而在波高為0.2 m時,推力系數(shù)隨時間的變化明顯較小,且呈現(xiàn)出明顯的周期性。由此可知,當(dāng)波周期=2 s時,較小的波高可使流場和推力系數(shù)趨于穩(wěn)定,較大的波高會導(dǎo)致流場與推力系數(shù)的不穩(wěn)定。

    由圖6(b)~(d)可知,波周期=4,6,8 s時仍具有上述規(guī)律,且隨著波周期的增大,流場和推力系數(shù)在更大波高下穩(wěn)定變化,其中: 波周期=4 s、波高=0.4 m; 波周期=6 s、波高=0.8 m; 波周期=8 s、波高=0.8 m時,推力系數(shù)趨于穩(wěn)定。

    筆者相信,以村落民俗志為基礎(chǔ)的鄉(xiāng)村研究,在當(dāng)今中國重述“亞洲”、重寫“世界史”和“全球史”的學(xué)術(shù)熱潮中,不僅沒有過時,而且不可或缺。畢竟,體察中國國情,理解中國道路,提煉中國話語,仍要以“在村落里研究”[注][美]克利福德·格爾茨:《文化的解釋》,韓莉譯,譯林出版社,1999年,第29頁。 的鄉(xiāng)土中國學(xué)術(shù)實(shí)踐為立足點(diǎn)。

    圖6 不同波高下翼推力系數(shù)隨時間變化曲線

    圖7 不同波高下翼平均推力系數(shù)變化曲線

    3.1.2 波周期

    圖8為翼推力系數(shù)在2個波周期隨時間變化的曲線,波高為=0.2,0.4,0.6,0.8 m?;谏瞎?jié)所述原因,選擇第2個周期進(jìn)行分析。

    由圖8(a)可知,推力系數(shù)隨波周期的增大而減小,且在波周期=4,6,8 s時,推力系數(shù)隨時間變化波動較小,呈現(xiàn)出明顯的周期性,流場趨于穩(wěn)定,而在波周期=2 s時,推力系數(shù)隨時間的變化波動較大,沒有呈現(xiàn)出明顯的周期性,流場并未趨于穩(wěn)定??梢?當(dāng)波高為0.2 m時,較大的波周期可使流場和推力系數(shù)趨于穩(wěn)定,較小的波周期會導(dǎo)致流場與推力系數(shù)的不穩(wěn)定。

    由圖8(b)~(d)可知,波高=0.4,0.6,0.8 m時仍具有上述規(guī)律,且隨著波高的增大,流場和推力系數(shù)在更大的波周期下穩(wěn)定變化,其中波高=0.4 m、波周期=4 s,波高=0.6 m、波周期=4 s,波高=0.8 m、波周期=6 s時,推力系數(shù)趨于穩(wěn)定。

    圖9為不同波周期下柔性水翼平均推力系數(shù)的變化曲線。由圖可知,不同波高下,平均推力系數(shù)均隨波周期的增大而減小,且當(dāng)波周期較小時,變化幅度較大,說明波周期對翼推進(jìn)性能的影響程度較大; 波周期較大時,變化幅度較小,說明波周期的變化對翼推進(jìn)性能的影響較小。該結(jié)果與文獻(xiàn)[5]相近。

    3.2 俯仰運(yùn)動參數(shù)對翼推進(jìn)性能的影響

    為分析俯仰運(yùn)動參數(shù)對柔性水翼推進(jìn)性能的影響,重點(diǎn)分析俯仰角幅值0和俯仰運(yùn)動周期2對推進(jìn)性能的影響,其中俯仰角幅值為0~15°,以5°間隔選取,俯仰運(yùn)動周期為0~4 s,以2 s的間隔選取,設(shè)計(jì)的具體工況如表3所示。

    圖8 不同波周期下翼推力系數(shù)變化曲線

    圖9 不同波周期下翼平均推力系數(shù)變化曲線

    表3 數(shù)值仿真俯仰運(yùn)動參數(shù)

    3.2.1 俯仰角幅值

    圖10為柔性水翼推力系數(shù)在俯仰運(yùn)動周期2=2,4 s時隨時間變化的曲線,選擇計(jì)算中的第2個周期進(jìn)行分析。由圖10(a)可知,在周期2=2 s時,隨俯仰角幅值的增大,推力系數(shù)出現(xiàn)負(fù)值的區(qū)域逐漸增大,隨時間變化的波動明顯增大,沒有呈現(xiàn)出明顯的周期性,流場并未趨于穩(wěn)定; 由圖10(b)可知,在周期2=4 s時,隨俯仰角幅值的增大,推力系數(shù)呈現(xiàn)良好的周期性,流場均趨于穩(wěn)定。

    圖10 不同俯仰角幅值下翼推力系數(shù)變化曲線

    圖11 不同俯仰角幅值下翼平均推力系數(shù)變化曲線

    圖11中,當(dāng)俯仰運(yùn)動周期2與波周期相等時,俯仰角幅值在一定范圍內(nèi)增大,平均推力系數(shù)隨之增大,變化幅度較小; 在俯仰運(yùn)動周期為波周期的1/2時,平均推力系數(shù)隨俯仰角幅值的增大而減小,變化幅度較大??梢?當(dāng)俯仰運(yùn)動周期與波周期相等時,在一定范圍內(nèi),俯仰角幅值越大,推進(jìn)性能越佳; 當(dāng)俯仰運(yùn)動周期為波周期的1/2時,在一定范圍內(nèi),俯仰角幅值越小,推進(jìn)性能越佳。

    3.2.2 俯仰運(yùn)動周期

    圖12為柔性水翼推力系數(shù)在2個波周期隨時間變化的曲線,俯仰角幅值0=5°,10°,15°時,選擇分析第2個周期。由圖可知,俯仰運(yùn)動周期的減小使得推力系數(shù)隨著時間變化的波動明顯增大,且周期性明顯變差,流場由較穩(wěn)定向較不穩(wěn)定轉(zhuǎn)變。

    圖12 不同俯仰運(yùn)動周期下翼推力系數(shù)變化曲線

    圖13為不同俯仰運(yùn)動周期下柔性水翼的平均推力系數(shù)變化曲線,由圖可知,平均推力系數(shù)均隨俯仰運(yùn)動周期的增大呈先減小后增大的趨勢,且在俯仰角幅值較小時,俯仰運(yùn)動周期對平均推力系數(shù)的影響較小。

    圖13 不同俯仰運(yùn)動周期下翼平均推力系數(shù)變化曲線

    3.3 斯特哈勞爾數(shù)對翼推進(jìn)性能的影響

    圖14為柔性水翼推力系數(shù)在2個波周期內(nèi)隨時間變化曲線,斯特哈勞爾數(shù)=0.450,0.225,0.150,0.1125,0.090,選擇計(jì)算中的第2個周期進(jìn)行分析。

    表4 數(shù)值仿真斯特哈勞爾數(shù)

    圖14 不同斯特哈勞爾數(shù)下翼推力系數(shù)變化曲線

    由圖可知,隨著斯特哈勞爾數(shù)的減小,推力系數(shù)隨時間變化的波動明顯減小,流場逐漸趨于穩(wěn)定,呈現(xiàn)出明顯的周期性,且當(dāng)斯特哈勞爾數(shù)較小時,對推力系數(shù)的影響較小。

    圖15為不同斯特哈勞爾數(shù)下柔性水翼的平均推力系數(shù)變化曲線,由圖可知,平均推力系數(shù)基本符合隨斯特哈勞爾數(shù)增大而增大,斯特哈勞爾數(shù)較小時,其改變對平均推力系數(shù)的影響較小,該結(jié)果與文獻(xiàn)[12]研究結(jié)果相近。

    圖15 不同斯特哈勞爾數(shù)下翼平均推力系數(shù)變化曲線

    4 水池實(shí)驗(yàn)驗(yàn)證

    4.1 實(shí)驗(yàn)平臺搭建

    水翼推進(jìn)實(shí)驗(yàn)平臺能夠?qū)崿F(xiàn)柔性水翼的升沉振蕩運(yùn)動,使得柔性水翼在流場作用下產(chǎn)生形變,從而產(chǎn)生前向推力。實(shí)驗(yàn)平臺主要由氣動裝置、主體框架、滑動導(dǎo)軌、拉力測試儀和柔性水翼機(jī)構(gòu)組成,通過調(diào)節(jié)氣動裝置的氣壓來調(diào)節(jié)柔性水翼升沉運(yùn)動的行程和周期,主體框架搭載氣動裝置和柔性水翼機(jī)構(gòu),在與柔性水翼前向相反的方向安裝拉力測試儀,實(shí)時接收柔性水翼機(jī)構(gòu)的前向推力數(shù)據(jù),實(shí)驗(yàn)平臺如圖16所示。

    圖16 實(shí)驗(yàn)平臺

    圖17 柔性水翼幾何參數(shù)

    圖18 柔性水翼實(shí)物圖

    4.2 實(shí)驗(yàn)結(jié)果分析

    圖19 波周期T=4 s下仿真與實(shí)驗(yàn)平均推力對比曲線

    盡管存在一定偏差,但仿真數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)趨勢性一致,符合前文仿真分析得出的規(guī)律,水池實(shí)驗(yàn)在一定程度上驗(yàn)證了文中的仿真結(jié)果。

    圖20 波周期T=6 s下仿真與實(shí)驗(yàn)平均推力對比曲線

    5 結(jié)論

    文中以波浪滑翔器中的柔性水翼為研究對象,對柔性水翼在波浪中的運(yùn)動進(jìn)行仿真分析,重點(diǎn)分析了運(yùn)動參數(shù)對翼推進(jìn)性能的影響,得到如下結(jié)論: 1) 通過COMSOL Multiphysics軟件的流體場與固體場相結(jié)合的模型,實(shí)現(xiàn)了柔性水翼在升沉運(yùn)動和俯仰運(yùn)動耦合中產(chǎn)生被動變形的仿真; 2) 在一定范圍內(nèi),翼推進(jìn)性能隨波高的增大而增大,且較大的波高易導(dǎo)致翼周圍的流場不穩(wěn)定,當(dāng)波周期較大時,波高對翼推進(jìn)性能的影響較小; 3) 在一定范圍內(nèi),翼推進(jìn)性能隨波周期的增大而減小,較小的波周期易導(dǎo)致翼周圍流場不穩(wěn)定,對翼推進(jìn)性能的影響較大; 4) 在一定范圍內(nèi),當(dāng)俯仰運(yùn)動周期與波周期相等時,俯仰運(yùn)動對翼推進(jìn)性能有積極作用,且俯仰角幅值越大,翼推進(jìn)性能越佳,當(dāng)俯仰運(yùn)動周期與波周期不相等時,俯仰運(yùn)動對翼推進(jìn)性能有消極作用,且俯仰角幅值越小,翼推進(jìn)性能越佳; 5) 在一定范圍內(nèi),翼推進(jìn)性能隨斯特哈勞爾數(shù)的減小而減小,且在低斯特哈勞爾數(shù)下對翼推進(jìn)性能影響較小; 6) 仿真數(shù)據(jù)和實(shí)驗(yàn)數(shù)據(jù)存在一定偏差,但二者趨勢性相近,一定程度上驗(yàn)證了文中的仿真結(jié)果。

    [1] Daniel T,Manley J,Trenaman N. The Wave Glider: Enabling a New Approach to Persistent Ocean Observation and Research[J]. Ocean Dynamics,2011,61(10): 1509-1520.

    [2] Hine R,Willcox S,Hine G,et al. The Wave Glider: A Wave-Powered Autonomous Marine Vehicle[C]//Procee- dings of MTS/IEEE Biloxi-Marine Technology for Our Future: Global and Local Challenges. Biloxi,USA: IEEE,2009: 296-301.

    [3] 廖煜雷,李曄,劉濤,等. 波浪滑翔器技術(shù)的回顧與展望[J]. 哈爾濱工程大學(xué)學(xué)報(bào),2016,37(9): 1227-1236.

    Liao Yu-lei,Li Ye,Liu Tao,et al. Unmanned Wave Glider Technology: State of Art and Perspective[J]. Journal of Harbin Engineering University,2016,37(9): 1227-1236.

    [4] 桑宏強(qiáng),李燦,孫秀軍. 波浪滑翔器縱向速度與波浪參數(shù)定量分析[J]. 水下無人系統(tǒng)學(xué)報(bào),2018,26(1): 16-22.

    Sang Hong-qiang,Li Can,Sun Xiu-jun. Quantitative Analysis on Longitudinal Velocity and Wave Parameter of Wave Glider[J]. Journal of Unmanned Undersea Systems,2018,26(1): 16-22.

    [5] 呂元博,田新亮,李欣,等. NACA 0012擺動撲翼水動力特性的二維數(shù)值模擬[J]. 中國艦船研究,2018,13(2): 7-15.

    Lü Yuan-bo,Tian Xin-liang,Li Xin,et al. Two-dimen- sional Numerical Simulation of NACA 0012 Flapping Foil Hydrodynamics[J]. Chinese Journal of Ship Research,2018,13(2): 7-15.

    [6] Liu P,Liu Y B,Huang S L,et al. Effects of Regular Waves on Propulsion Performance of Flexible Flapping Foil[J]. Applied Sciences,2018,8(6): 934.

    [7] Liu W,Xiao Q,Zhu Q. Passive Flexibility Effect on Oscillating Foil Energy Harvester[J]. AIAA Journal,2016,54(4): 1-16.

    [8] Alben S,Witt C,Baker T V,et al. Dynamics of Freely Swimming Flexible Foils[J]. Physics of Fluids,2012,24(5): 109-133.

    [9] Prempraneerach P,Hober F S,Triantfyllou M S. The Effect of Chordwise Flexibility on the Trust and Efficiency of a Flapping foil[C]//Proceedings of the 13th International Symposium on Unmanned Untethered Submersible Technology. Durham,USA: IEEE,2003: 130-128.

    [10] Katz J,Weihs D. Hydrodynamic Propulsion by Large Amplitude Oscillation of an Airfoil with Chordwise Flexibility[J]. Journal of Fluid Mechanics,1978,88(3): 485-497.

    [11] 田寶強(qiáng),李玲瓏. 蹼翼型波浪滑翔機(jī)結(jié)構(gòu)設(shè)計(jì)和運(yùn)動原理分析[J]. 中國機(jī)械工程,2017,28(24): 2939-2943.

    Tian Bao-qiang,Li Ling-long. Structure Design and Movement Principle of Wave Gliders with Webbed Wings[J]. China Mechanical Engineering,2017,28(24): 2939-2943.

    [12] 周凱. 典型仿生水翼的推進(jìn)及能量吸收機(jī)制研究[D]. 哈爾濱: 哈爾濱工程大學(xué),2018.

    [13] Kinsey T,Dumas G. Computational Fluid Dynamics Analysis of a Hydrokinetic Turbine Based on Oscillating Hydrofoils[J]. Journal of Fluids Engineering,2012,134(2): 021104.

    Numerical Simulation Analysis of Flexible Hydrofoil Propulsion Performance for Wave Glider

    1,1,3,4*,1,2

    (1. School of Mechanical Engineering,Tianjin Polytechnic University,Tianjin 300387,China; 2.Tianjin Key Laboratory of Advanced Mechatronic Equipment Technology,Tianjin 30087,China; 3. Physical Oceanorgraphy Laboratory,Ocean University of China,Qingdao 266100,China; 4. Laboratory of Marine Dynamics and Climate Function,Pilot National Laboratory for Marine Science and Technology(Qingdao),Qingdao 266237,China)

    The hydrofoil is a key component that determines the capability of wave power conversions and the navigation performance of wave gliders. The propulsion performance of a flexible hydrofoil is investigated in this study. The passive deformation process of a flexible hydrofoil under active heave motion is simulated using the COMSOL Multiphysics software according to the motion characteristics of a hydrofoil when the wave glider is working. The effects of wave and pitching motion parameters and the Strouhal number on the propulsion performance of the hydrofoil are investigated. The simulation results show that the propulsion performance of the hydrofoil increases with wave height and Strouhal number and decreases with wave period within a certain range. When the wave period is high,the change in wave height and wave period has little effect on the propulsion performance of the hydrofoil. In contrast,when the pitch and wave periods are equal,the pitching motion has a positive effect on the performance.

    wave glider; flexible hydrofoil; propulsion performance; numerical simulation

    TJ630; U674.941; O353

    A

    2096-3920(2021)04-0374-09

    10.11993/j.issn.2096-3920.2021.04.002

    劉芬,彭彬,孫秀軍,等. 波浪滑翔器柔性水翼推進(jìn)性能數(shù)值仿真分析[J]. 水下無人系統(tǒng)學(xué)報(bào),2021,29(4): 374-382.

    2020-08-28;

    2020-11-11.

    山東省重大科技創(chuàng)新項(xiàng)目(2019JZZY020701);天津市自然科學(xué)基金重點(diǎn)基金(18JCZDJC40100); 國家重點(diǎn)研發(fā)計(jì)劃重點(diǎn)專項(xiàng)(2017YFC0305902); 青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室“問海計(jì)劃”項(xiàng)目(2017WHZZB0101).

    通信作者簡介:孫秀軍(1981-),男,博士,教授,主要研究方向?yàn)椴ɡ嘶杵骷夹g(shù)及其海洋觀測應(yīng)用.

    (責(zé)任編輯: 楊力軍)

    猜你喜歡
    水翼滑翔波浪
    蘭花螳螂會滑翔
    波浪谷和波浪巖
    攻天掠地的先鋒武器——滑翔導(dǎo)彈
    波浪滑翔機(jī)橢圓形后緣水翼動力特性研究
    袖珍水翼突防潛艇的設(shè)計(jì)構(gòu)想及運(yùn)用研究
    波浪谷隨想
    去看神奇波浪谷
    讓你心跳加速的“滑翔飛板”
    三維扭曲水翼空化現(xiàn)象CFD模擬
    空中滑翔大比拼(下)——滑翔傘
    国产av一区在线观看免费| 我要看日韩黄色一级片| 国产三级中文精品| 人人妻人人澡人人爽人人夜夜 | 国产精品久久久久久亚洲av鲁大| 亚洲精品亚洲一区二区| 日韩中字成人| 亚洲一级一片aⅴ在线观看| 精品一区二区免费观看| 99久久精品热视频| 亚洲自偷自拍三级| 色哟哟·www| 欧美性猛交╳xxx乱大交人| 欧美日韩乱码在线| 成人三级黄色视频| 日韩欧美 国产精品| 久久久久国产网址| 岛国毛片在线播放| 伊人久久精品亚洲午夜| 国产亚洲精品久久久com| 亚洲成人av在线免费| 晚上一个人看的免费电影| 亚洲五月天丁香| 久久久成人免费电影| 不卡一级毛片| 联通29元200g的流量卡| 亚洲一区高清亚洲精品| 卡戴珊不雅视频在线播放| www.色视频.com| 亚洲精品久久久久久婷婷小说 | 亚洲自偷自拍三级| 日韩人妻高清精品专区| 黄色配什么色好看| 欧美最黄视频在线播放免费| 一级黄片播放器| 亚洲婷婷狠狠爱综合网| 国产一区二区激情短视频| 亚洲精品自拍成人| 亚洲在线观看片| 欧美丝袜亚洲另类| 少妇猛男粗大的猛烈进出视频 | 一个人看的www免费观看视频| 亚洲人成网站高清观看| 国产淫片久久久久久久久| 国产精品.久久久| 麻豆乱淫一区二区| 少妇的逼好多水| 小说图片视频综合网站| 国产成人a区在线观看| 久久久久网色| av又黄又爽大尺度在线免费看 | 97热精品久久久久久| 国产亚洲精品久久久久久毛片| 99久久人妻综合| 亚洲av免费高清在线观看| 天天躁夜夜躁狠狠久久av| 中国美白少妇内射xxxbb| 青春草国产在线视频 | 亚洲国产欧洲综合997久久,| 免费av不卡在线播放| 亚洲第一电影网av| 日韩一区二区三区影片| 亚洲国产色片| 亚洲七黄色美女视频| 国产一区二区在线观看日韩| 两性午夜刺激爽爽歪歪视频在线观看| 久久鲁丝午夜福利片| 1024手机看黄色片| 黑人高潮一二区| 亚洲国产高清在线一区二区三| 精品99又大又爽又粗少妇毛片| av视频在线观看入口| av在线蜜桃| 18+在线观看网站| 久久草成人影院| 51国产日韩欧美| 久久久久久久久大av| 变态另类成人亚洲欧美熟女| 国产av一区在线观看免费| 日韩亚洲欧美综合| 国产一区二区激情短视频| 日韩一本色道免费dvd| 中国美白少妇内射xxxbb| 中文字幕制服av| 国产精品女同一区二区软件| 91aial.com中文字幕在线观看| 99久久九九国产精品国产免费| 欧洲精品卡2卡3卡4卡5卡区| 能在线免费看毛片的网站| 国产蜜桃级精品一区二区三区| 人体艺术视频欧美日本| 天堂网av新在线| 久久久国产成人免费| 日韩强制内射视频| 欧美日本亚洲视频在线播放| 久久久久久久久久成人| 内地一区二区视频在线| 久久国内精品自在自线图片| 亚洲国产欧美人成| 日韩欧美精品v在线| 欧美激情国产日韩精品一区| 亚洲不卡免费看| 久久午夜亚洲精品久久| 日韩精品有码人妻一区| 在线观看66精品国产| 久久精品久久久久久久性| 国产精品女同一区二区软件| 国产三级在线视频| 婷婷亚洲欧美| 国产成人精品一,二区 | 亚洲人成网站高清观看| 禁无遮挡网站| 欧美一区二区国产精品久久精品| 精品少妇黑人巨大在线播放 | 美女黄网站色视频| 国产精品国产高清国产av| 99久久成人亚洲精品观看| 一本一本综合久久| 久久精品影院6| 亚洲人成网站在线播放欧美日韩| 精品国内亚洲2022精品成人| 亚洲欧美清纯卡通| 久久精品国产清高在天天线| 国产精品一区二区三区四区免费观看| 成人二区视频| 成人午夜高清在线视频| 日韩中字成人| 天堂中文最新版在线下载 | 亚洲欧美精品综合久久99| 六月丁香七月| 一个人免费在线观看电影| 久久综合国产亚洲精品| 在线观看午夜福利视频| 成人av在线播放网站| 91精品国产九色| 免费黄网站久久成人精品| 高清毛片免费观看视频网站| 九色成人免费人妻av| 久久国产乱子免费精品| 日韩av不卡免费在线播放| 伦理电影大哥的女人| 亚洲在线自拍视频| 99久久成人亚洲精品观看| 国产精品一二三区在线看| 毛片一级片免费看久久久久| 我要看日韩黄色一级片| 日韩中字成人| 日韩人妻高清精品专区| .国产精品久久| 天堂影院成人在线观看| 热99re8久久精品国产| 99国产精品一区二区蜜桃av| 欧美性猛交黑人性爽| www日本黄色视频网| av在线蜜桃| 国产极品精品免费视频能看的| 免费观看精品视频网站| 99久久精品国产国产毛片| 国产一区二区在线av高清观看| 99国产极品粉嫩在线观看| 欧美不卡视频在线免费观看| 蜜桃亚洲精品一区二区三区| 国产成人a∨麻豆精品| 国产一区亚洲一区在线观看| 日产精品乱码卡一卡2卡三| 少妇熟女aⅴ在线视频| 午夜久久久久精精品| 男人的好看免费观看在线视频| 国产成人精品婷婷| 你懂的网址亚洲精品在线观看 | 欧美成人一区二区免费高清观看| 日韩人妻高清精品专区| 日韩国内少妇激情av| 午夜久久久久精精品| 国产亚洲精品久久久久久毛片| 久久久久性生活片| 久久久久久久午夜电影| 男女做爰动态图高潮gif福利片| 91午夜精品亚洲一区二区三区| 精品不卡国产一区二区三区| 亚洲欧美精品综合久久99| 美女cb高潮喷水在线观看| 久久精品国产自在天天线| 变态另类成人亚洲欧美熟女| 久久久午夜欧美精品| 欧美激情久久久久久爽电影| 欧美另类亚洲清纯唯美| 日韩欧美 国产精品| 国内久久婷婷六月综合欲色啪| 日日干狠狠操夜夜爽| 丰满的人妻完整版| 99国产极品粉嫩在线观看| 欧美日韩一区二区视频在线观看视频在线 | 在线a可以看的网站| 美女高潮的动态| 一夜夜www| 亚洲18禁久久av| 最好的美女福利视频网| 色哟哟哟哟哟哟| 夫妻性生交免费视频一级片| 色尼玛亚洲综合影院| 久久久久九九精品影院| 久久九九热精品免费| 中文字幕精品亚洲无线码一区| av视频在线观看入口| 日韩成人av中文字幕在线观看| 欧美成人a在线观看| 国产视频首页在线观看| 日韩制服骚丝袜av| 2021天堂中文幕一二区在线观| 国产av麻豆久久久久久久| 国产单亲对白刺激| 国产日本99.免费观看| 午夜激情欧美在线| 成人av在线播放网站| 成人无遮挡网站| 国产黄片视频在线免费观看| 亚洲最大成人av| 99热只有精品国产| 国内揄拍国产精品人妻在线| 久久99蜜桃精品久久| eeuss影院久久| 国产精品久久久久久精品电影| 国产成人91sexporn| 搡老妇女老女人老熟妇| 亚洲精品色激情综合| 精华霜和精华液先用哪个| 一区二区三区高清视频在线| 亚洲精品久久久久久婷婷小说 | 有码 亚洲区| 久久久久久久久久黄片| 午夜a级毛片| 黑人高潮一二区| 亚洲精品影视一区二区三区av| 免费观看人在逋| 噜噜噜噜噜久久久久久91| 五月玫瑰六月丁香| 日本-黄色视频高清免费观看| av在线老鸭窝| 成人鲁丝片一二三区免费| 亚洲七黄色美女视频| 免费观看a级毛片全部| 国产激情偷乱视频一区二区| 热99re8久久精品国产| 天美传媒精品一区二区| 99久久中文字幕三级久久日本| 亚洲激情五月婷婷啪啪| 黄色日韩在线| 看黄色毛片网站| 成年免费大片在线观看| 欧美xxxx黑人xx丫x性爽| 美女大奶头视频| 婷婷精品国产亚洲av| 欧美潮喷喷水| 丝袜喷水一区| 中文在线观看免费www的网站| 国产黄色视频一区二区在线观看 | 欧美一级a爱片免费观看看| 亚洲精品456在线播放app| 最近的中文字幕免费完整| 岛国在线免费视频观看| 嫩草影院精品99| 中文字幕熟女人妻在线| 国产精品精品国产色婷婷| 超碰av人人做人人爽久久| 久久精品国产自在天天线| 久久精品久久久久久噜噜老黄 | 成人一区二区视频在线观看| 亚洲不卡免费看| 丝袜喷水一区| 久久久久久久久中文| 亚洲无线观看免费| 成人美女网站在线观看视频| 亚洲在线自拍视频| 真实男女啪啪啪动态图| 最近最新中文字幕大全电影3| 国产成人a区在线观看| 亚洲国产精品国产精品| 日韩视频在线欧美| 成人综合一区亚洲| 国产精品电影一区二区三区| 久久精品国产亚洲av天美| 亚洲av电影不卡..在线观看| 亚洲精品乱码久久久v下载方式| 久久久久久久久久久丰满| 尤物成人国产欧美一区二区三区| 乱码一卡2卡4卡精品| 成人午夜精彩视频在线观看| 国产毛片a区久久久久| 国产熟女欧美一区二区| 内射极品少妇av片p| 国产精品.久久久| 熟女电影av网| 欧美区成人在线视频| 国产精品久久电影中文字幕| 夜夜夜夜夜久久久久| 狂野欧美白嫩少妇大欣赏| 久99久视频精品免费| 99riav亚洲国产免费| 欧美色视频一区免费| 久99久视频精品免费| 高清毛片免费观看视频网站| 午夜福利视频1000在线观看| 日本一二三区视频观看| 午夜精品在线福利| 午夜爱爱视频在线播放| 亚洲国产色片| 亚洲国产精品成人综合色| 观看美女的网站| 国产一区亚洲一区在线观看| 男女做爰动态图高潮gif福利片| 97超视频在线观看视频| 色哟哟哟哟哟哟| ponron亚洲| 久久午夜亚洲精品久久| 岛国在线免费视频观看| 久久午夜福利片| 亚洲va在线va天堂va国产| 一本久久精品| 99在线视频只有这里精品首页| 99热6这里只有精品| 校园春色视频在线观看| 校园春色视频在线观看| 久久精品国产鲁丝片午夜精品| 最近最新中文字幕大全电影3| 午夜久久久久精精品| 成人一区二区视频在线观看| 尤物成人国产欧美一区二区三区| 啦啦啦啦在线视频资源| 日本黄色片子视频| 亚洲人成网站在线播| 国产在线精品亚洲第一网站| 又黄又爽又刺激的免费视频.| 欧美潮喷喷水| 熟妇人妻久久中文字幕3abv| 日本熟妇午夜| 网址你懂的国产日韩在线| 99视频精品全部免费 在线| 不卡一级毛片| 只有这里有精品99| 国产精品无大码| 久久精品国产亚洲av香蕉五月| 免费电影在线观看免费观看| 久久精品影院6| 哪里可以看免费的av片| 九色成人免费人妻av| 精品久久国产蜜桃| 成人高潮视频无遮挡免费网站| 99久久精品国产国产毛片| 一本久久中文字幕| 成人av在线播放网站| 亚洲av熟女| 国产成人影院久久av| 久久人人爽人人片av| 亚洲av不卡在线观看| 又爽又黄a免费视频| 国产伦精品一区二区三区四那| 亚洲无线观看免费| 亚洲国产精品合色在线| 人体艺术视频欧美日本| 国产不卡一卡二| 69av精品久久久久久| 成人二区视频| 亚洲av电影不卡..在线观看| 青春草视频在线免费观看| 亚洲性久久影院| 不卡一级毛片| 人人妻人人澡欧美一区二区| 久久婷婷人人爽人人干人人爱| 亚洲av成人av| 天美传媒精品一区二区| 久久婷婷人人爽人人干人人爱| 亚洲国产欧美在线一区| 国产黄色小视频在线观看| 精品久久久久久久末码| 简卡轻食公司| 嘟嘟电影网在线观看| 成人二区视频| 99精品在免费线老司机午夜| www.av在线官网国产| 国产不卡一卡二| 一边摸一边抽搐一进一小说| 国产高清三级在线| 三级男女做爰猛烈吃奶摸视频| av免费观看日本| 亚洲aⅴ乱码一区二区在线播放| 国产黄色视频一区二区在线观看 | 亚洲av电影不卡..在线观看| 日本-黄色视频高清免费观看| 一边亲一边摸免费视频| 免费看av在线观看网站| 欧美最新免费一区二区三区| 爱豆传媒免费全集在线观看| 久久精品夜色国产| 国产精品.久久久| 色哟哟·www| 亚洲欧美精品综合久久99| 国产一级毛片七仙女欲春2| 激情 狠狠 欧美| 麻豆乱淫一区二区| 欧美最黄视频在线播放免费| 亚洲欧美日韩东京热| 国产精品1区2区在线观看.| 国产午夜福利久久久久久| 蜜臀久久99精品久久宅男| 九九久久精品国产亚洲av麻豆| 麻豆国产97在线/欧美| 成人午夜高清在线视频| 99热6这里只有精品| 十八禁国产超污无遮挡网站| 国产日韩欧美在线精品| 亚洲成人av在线免费| 色综合亚洲欧美另类图片| 亚洲婷婷狠狠爱综合网| 简卡轻食公司| 午夜精品国产一区二区电影 | 午夜精品国产一区二区电影 | 毛片一级片免费看久久久久| 一本一本综合久久| 天美传媒精品一区二区| 久久久久久国产a免费观看| 亚洲av第一区精品v没综合| 可以在线观看的亚洲视频| 天堂影院成人在线观看| 久久这里有精品视频免费| 成人综合一区亚洲| h日本视频在线播放| 日韩一本色道免费dvd| 91麻豆精品激情在线观看国产| 欧美性感艳星| 国产亚洲av嫩草精品影院| 精品久久久久久成人av| 亚洲国产日韩欧美精品在线观看| 网址你懂的国产日韩在线| 欧美又色又爽又黄视频| 国产亚洲91精品色在线| а√天堂www在线а√下载| 久久人人爽人人片av| 国产精品电影一区二区三区| 在线天堂最新版资源| 中文欧美无线码| 男女边吃奶边做爰视频| 亚洲精品国产av成人精品| 校园春色视频在线观看| 91精品一卡2卡3卡4卡| 久久热精品热| 我的女老师完整版在线观看| 丝袜美腿在线中文| 亚洲国产色片| 久久久久久久亚洲中文字幕| 国产毛片a区久久久久| 午夜福利高清视频| 国产精品国产三级国产av玫瑰| 岛国在线免费视频观看| 亚洲欧美中文字幕日韩二区| 精品日产1卡2卡| 日韩,欧美,国产一区二区三区 | 欧美日韩在线观看h| 黑人高潮一二区| 国产黄片视频在线免费观看| 啦啦啦啦在线视频资源| 麻豆一二三区av精品| 国产成人精品久久久久久| 丰满的人妻完整版| 69人妻影院| 欧美3d第一页| 国产精品野战在线观看| 亚洲精品亚洲一区二区| 国产精品,欧美在线| 一区二区三区免费毛片| 内地一区二区视频在线| 三级男女做爰猛烈吃奶摸视频| avwww免费| 最近中文字幕高清免费大全6| 久久久久性生活片| av在线老鸭窝| 熟女电影av网| 成人国产麻豆网| 欧美在线一区亚洲| 精品日产1卡2卡| 欧美激情久久久久久爽电影| 国产一区二区三区av在线 | 又爽又黄无遮挡网站| 赤兔流量卡办理| 2021天堂中文幕一二区在线观| 九九热线精品视视频播放| 国产成人aa在线观看| 大又大粗又爽又黄少妇毛片口| 12—13女人毛片做爰片一| 丰满乱子伦码专区| 嘟嘟电影网在线观看| 国产成人精品久久久久久| 99热网站在线观看| 精品欧美国产一区二区三| 人妻制服诱惑在线中文字幕| 国产精品久久久久久精品电影| 亚洲第一电影网av| 亚洲三级黄色毛片| 日韩一本色道免费dvd| 国产成人a∨麻豆精品| 伊人久久精品亚洲午夜| 国产伦精品一区二区三区视频9| 亚洲av免费高清在线观看| 成人综合一区亚洲| 男人舔女人下体高潮全视频| 国内精品宾馆在线| 搞女人的毛片| 亚洲国产欧美在线一区| 啦啦啦啦在线视频资源| 不卡一级毛片| 国产高清有码在线观看视频| 亚洲在线观看片| 久久99精品国语久久久| 一本久久中文字幕| 伦精品一区二区三区| 久久久久久久久久黄片| 51国产日韩欧美| 国产毛片a区久久久久| 黄色视频,在线免费观看| 国产色婷婷99| 深夜a级毛片| 欧美丝袜亚洲另类| 欧美精品一区二区大全| 美女 人体艺术 gogo| av在线蜜桃| 性插视频无遮挡在线免费观看| 精品久久久噜噜| 国产精品久久久久久亚洲av鲁大| 成人漫画全彩无遮挡| 国产 一区 欧美 日韩| 男的添女的下面高潮视频| 六月丁香七月| 国产精品人妻久久久久久| 波多野结衣高清无吗| 午夜亚洲福利在线播放| 亚洲丝袜综合中文字幕| 又粗又爽又猛毛片免费看| a级毛片a级免费在线| 天堂网av新在线| 久久99热这里只有精品18| 啦啦啦啦在线视频资源| 亚洲精品日韩在线中文字幕 | 国产精品人妻久久久久久| 波多野结衣高清作品| 特大巨黑吊av在线直播| 亚洲国产高清在线一区二区三| 色播亚洲综合网| 伦精品一区二区三区| 九九热线精品视视频播放| 精品久久久噜噜| 热99在线观看视频| 色视频www国产| 身体一侧抽搐| av卡一久久| 中国美白少妇内射xxxbb| 国产高清视频在线观看网站| 欧美丝袜亚洲另类| 男的添女的下面高潮视频| 最近2019中文字幕mv第一页| 天天躁夜夜躁狠狠久久av| 三级毛片av免费| 天天躁夜夜躁狠狠久久av| 日韩欧美精品免费久久| 日日摸夜夜添夜夜添av毛片| 亚洲精品粉嫩美女一区| 久久久久久久久久久丰满| 村上凉子中文字幕在线| 国产精品一区www在线观看| 变态另类成人亚洲欧美熟女| av在线天堂中文字幕| 99久久精品一区二区三区| 亚洲av男天堂| 校园春色视频在线观看| 日本与韩国留学比较| 人妻少妇偷人精品九色| 岛国毛片在线播放| 观看美女的网站| 久久亚洲国产成人精品v| 丰满人妻一区二区三区视频av| 午夜精品国产一区二区电影 | 国产在线精品亚洲第一网站| 亚洲欧洲国产日韩| 午夜精品国产一区二区电影 | 国产精品电影一区二区三区| 深夜精品福利| 国产色爽女视频免费观看| eeuss影院久久| 又爽又黄无遮挡网站| 久久99蜜桃精品久久| 99热6这里只有精品| 午夜福利在线在线| 午夜精品在线福利| 国产亚洲精品久久久久久毛片| 欧美zozozo另类| 亚洲av二区三区四区| av视频在线观看入口| 色综合色国产| а√天堂www在线а√下载| www.av在线官网国产| 亚洲国产精品成人综合色| 久久久久国产网址| 亚洲七黄色美女视频| 岛国毛片在线播放| av视频在线观看入口| 一个人看的www免费观看视频| 天天躁日日操中文字幕| 国产午夜精品论理片| 99九九线精品视频在线观看视频| 在线国产一区二区在线| a级毛片免费高清观看在线播放| 成人欧美大片| 欧美一区二区国产精品久久精品| 精品久久久久久久久久久久久| 久久婷婷人人爽人人干人人爱| 99久久人妻综合|