• <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手机在线免费观看| 午夜精品久久久久久毛片777| 亚洲人成77777在线视频| 久久久精品免费免费高清| 免费日韩欧美在线观看| 亚洲自偷自拍图片 自拍| 欧美精品av麻豆av| 久久免费观看电影| 汤姆久久久久久久影院中文字幕| 欧美在线一区亚洲| 这个男人来自地球电影免费观看| 又大又爽又粗| 热99re8久久精品国产| 俄罗斯特黄特色一大片| 少妇裸体淫交视频免费看高清 | 91成年电影在线观看| 大陆偷拍与自拍| 成年动漫av网址| 国产主播在线观看一区二区| 国产伦理片在线播放av一区| 露出奶头的视频| 久久狼人影院| 国产av一区二区精品久久| 精品人妻1区二区| 在线观看免费视频网站a站| 亚洲第一青青草原| 成人免费观看视频高清| 欧美日韩中文字幕国产精品一区二区三区 | 青草久久国产| 国产亚洲精品久久久久5区| 亚洲专区中文字幕在线| 老司机靠b影院| 精品国产一区二区三区四区第35| 日韩欧美三级三区| 亚洲欧美激情在线| 九色亚洲精品在线播放| 成年人午夜在线观看视频| av不卡在线播放| 麻豆成人av在线观看| 亚洲专区国产一区二区| 九色亚洲精品在线播放| 国产高清视频在线播放一区| 免费日韩欧美在线观看| 老司机靠b影院| 中文亚洲av片在线观看爽 | 中文字幕精品免费在线观看视频| 日韩免费高清中文字幕av| 国产精品1区2区在线观看. | 午夜免费鲁丝| 欧美中文综合在线视频| 久久99热这里只频精品6学生| 亚洲色图综合在线观看| 欧美另类亚洲清纯唯美| 日日夜夜操网爽| 欧美日韩av久久| 悠悠久久av| 午夜视频精品福利| 嫩草影视91久久| 亚洲精品久久午夜乱码| 99国产精品一区二区三区| 久久久久久亚洲精品国产蜜桃av| 亚洲情色 制服丝袜| av国产精品久久久久影院| 一本色道久久久久久精品综合| 久久久久精品国产欧美久久久| 91精品国产国语对白视频| 99riav亚洲国产免费| 老鸭窝网址在线观看| 久久久精品免费免费高清| 亚洲精品中文字幕在线视频| 精品高清国产在线一区| 这个男人来自地球电影免费观看| 午夜福利在线免费观看网站| 丰满饥渴人妻一区二区三| 热re99久久国产66热| 桃花免费在线播放| 人成视频在线观看免费观看| 在线观看66精品国产| 精品国产乱子伦一区二区三区| 露出奶头的视频| 国产精品久久久久久人妻精品电影 | 国产精品成人在线| 亚洲欧美一区二区三区黑人| 老鸭窝网址在线观看| 日本一区二区免费在线视频| 桃红色精品国产亚洲av| 一级,二级,三级黄色视频| 久久久久国内视频| 人人妻人人澡人人看| 国产精品秋霞免费鲁丝片| 久久这里只有精品19| av一本久久久久| 一级毛片精品| 国产成人啪精品午夜网站| 99九九在线精品视频| 精品少妇内射三级| 成人特级黄色片久久久久久久 | 国产伦人伦偷精品视频| videosex国产| 人成视频在线观看免费观看| 极品教师在线免费播放| tocl精华| 亚洲国产欧美在线一区| 精品少妇内射三级| 老司机靠b影院| 啦啦啦免费观看视频1| 99在线人妻在线中文字幕 | 国产真人三级小视频在线观看| 国产精品.久久久| 久久亚洲精品不卡| 极品教师在线免费播放| 女人久久www免费人成看片| 亚洲人成电影观看| 久久久精品国产亚洲av高清涩受| 日本a在线网址| 亚洲伊人久久精品综合| 亚洲国产精品一区二区三区在线| 另类亚洲欧美激情| 大香蕉久久网| 日韩视频在线欧美| 夜夜爽天天搞| 极品教师在线免费播放| 欧美激情高清一区二区三区| 久久国产精品大桥未久av| 香蕉国产在线看| 亚洲av日韩在线播放| bbb黄色大片| 一区在线观看完整版| 国产成人一区二区三区免费视频网站| 91麻豆av在线| 午夜福利视频精品| av免费在线观看网站| 精品国产国语对白av| 亚洲国产欧美日韩在线播放| 国产不卡av网站在线观看| 老熟妇仑乱视频hdxx| 国产又爽黄色视频| 国产精品久久久av美女十八| 久久精品国产亚洲av香蕉五月 | 啦啦啦 在线观看视频| 十八禁人妻一区二区| 久久久国产一区二区| 咕卡用的链子| 色播在线永久视频| 国产精品国产高清国产av | 啦啦啦中文免费视频观看日本| 精品少妇一区二区三区视频日本电影| 欧美日韩av久久| 中文亚洲av片在线观看爽 | 午夜福利一区二区在线看| 日本精品一区二区三区蜜桃| 波多野结衣av一区二区av| 国产精品九九99| 大香蕉久久网| 超色免费av| 国产成人精品在线电影| 1024视频免费在线观看| 桃花免费在线播放| 久久久久久久久免费视频了| 国产成人精品久久二区二区免费| 亚洲午夜精品一区,二区,三区| 丁香六月天网| 国产一区二区三区视频了| 国产免费av片在线观看野外av| 精品人妻在线不人妻| 亚洲一区二区三区欧美精品| 19禁男女啪啪无遮挡网站| 亚洲天堂av无毛| xxxhd国产人妻xxx| 麻豆成人av在线观看| 日韩精品免费视频一区二区三区| 久久久国产一区二区| 国产成人av激情在线播放| 色婷婷久久久亚洲欧美| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利影视在线免费观看| 自线自在国产av| 91成人精品电影| 免费在线观看日本一区| 真人做人爱边吃奶动态| 日日摸夜夜添夜夜添小说| 不卡一级毛片| 亚洲成人免费av在线播放| 亚洲成av片中文字幕在线观看| 五月开心婷婷网| 亚洲av日韩在线播放| 亚洲色图 男人天堂 中文字幕| 男女免费视频国产| aaaaa片日本免费| 亚洲欧洲精品一区二区精品久久久| 亚洲av成人一区二区三| 亚洲全国av大片| 91成年电影在线观看| 高清视频免费观看一区二区| 麻豆av在线久日| 麻豆国产av国片精品| 少妇的丰满在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产精品一区二区在线不卡| av不卡在线播放| 久久精品国产a三级三级三级| 欧美日韩中文字幕国产精品一区二区三区 | 日日摸夜夜添夜夜添小说| 亚洲五月色婷婷综合| 国产深夜福利视频在线观看| 国产一区二区激情短视频| 成年人免费黄色播放视频| 亚洲熟女精品中文字幕| 成人手机av| 国产成人av激情在线播放| 久久久久网色| 日本欧美视频一区| 日韩人妻精品一区2区三区| 夫妻午夜视频| 一本久久精品| av欧美777| 高清毛片免费观看视频网站 | 搡老岳熟女国产| videosex国产| 久久这里只有精品19| 亚洲午夜理论影院| av线在线观看网站| 日韩熟女老妇一区二区性免费视频| 69精品国产乱码久久久| 女人高潮潮喷娇喘18禁视频| 高清视频免费观看一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产aⅴ精品一区二区三区波| 精品一区二区三区视频在线观看免费 | 在线天堂中文资源库| 亚洲精华国产精华精| www.熟女人妻精品国产| 变态另类成人亚洲欧美熟女 | 母亲3免费完整高清在线观看| www.精华液| 亚洲精品一卡2卡三卡4卡5卡| 国产不卡av网站在线观看| 精品免费久久久久久久清纯 | 啦啦啦在线免费观看视频4| 亚洲精品中文字幕在线视频| 女性生殖器流出的白浆| 操出白浆在线播放| 成人精品一区二区免费| 嫩草影视91久久| 狠狠狠狠99中文字幕| av天堂久久9| 欧美日韩黄片免| 女性生殖器流出的白浆| 亚洲成国产人片在线观看| 黄色视频在线播放观看不卡| 日韩欧美三级三区| 久久久国产成人免费| 久久 成人 亚洲| 午夜免费成人在线视频| 免费观看人在逋| 国产精品免费大片| 中文亚洲av片在线观看爽 | 免费女性裸体啪啪无遮挡网站| 一级片'在线观看视频| 欧美激情久久久久久爽电影 | 1024视频免费在线观看| 国产精品一区二区在线不卡| 超碰97精品在线观看| 午夜福利免费观看在线| 精品国产乱码久久久久久小说| 99香蕉大伊视频| 少妇猛男粗大的猛烈进出视频| tube8黄色片| 欧美精品人与动牲交sv欧美| 黑人猛操日本美女一级片| 亚洲专区国产一区二区| 国产欧美日韩精品亚洲av| 亚洲天堂av无毛| 欧美大码av| 亚洲国产欧美日韩在线播放| 精品少妇黑人巨大在线播放| 在线观看免费午夜福利视频| 中文字幕高清在线视频| 国产成人系列免费观看| 视频在线观看一区二区三区| 精品一区二区三区av网在线观看 | 久久人妻熟女aⅴ| 国产成人啪精品午夜网站| 999久久久精品免费观看国产| 久久久水蜜桃国产精品网| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久精品电影小说| 亚洲精品一二三| 欧美日韩一级在线毛片| 亚洲专区字幕在线| 另类亚洲欧美激情| 中文字幕人妻熟女乱码| 成人永久免费在线观看视频 | 色婷婷av一区二区三区视频| 亚洲va日本ⅴa欧美va伊人久久| av在线播放免费不卡| 一级,二级,三级黄色视频| 成人精品一区二区免费| 国产又爽黄色视频| 又大又爽又粗| 色播在线永久视频| 啦啦啦免费观看视频1| 国产精品偷伦视频观看了| 91精品国产国语对白视频| 色在线成人网| 99国产精品一区二区蜜桃av | 精品少妇内射三级| 天天躁狠狠躁夜夜躁狠狠躁| 大香蕉久久成人网| 91精品国产国语对白视频| 精品国产一区二区三区四区第35| 欧美亚洲 丝袜 人妻 在线| 丝袜美足系列| 成年动漫av网址| 91国产中文字幕| 午夜91福利影院| 国产高清激情床上av| 成年人黄色毛片网站| 色播在线永久视频| 他把我摸到了高潮在线观看 | 中文字幕制服av| 日韩视频一区二区在线观看| 久久久久视频综合| 国产亚洲一区二区精品| 国产免费现黄频在线看| 精品午夜福利视频在线观看一区 | 亚洲成av片中文字幕在线观看| 首页视频小说图片口味搜索| 又黄又粗又硬又大视频| 久久av网站| 视频区图区小说| 国产精品久久久久久精品电影小说| kizo精华| 成人国产av品久久久| 高清视频免费观看一区二区| 在线观看免费午夜福利视频| 精品卡一卡二卡四卡免费| 另类亚洲欧美激情| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩国产mv在线观看视频| 国产一区有黄有色的免费视频| 两性夫妻黄色片| 欧美国产精品va在线观看不卡| 欧美黑人欧美精品刺激| 久久久精品国产亚洲av高清涩受| 一区二区三区乱码不卡18| 91麻豆精品激情在线观看国产 | 亚洲精品国产精品久久久不卡| 免费黄频网站在线观看国产| 18在线观看网站| 亚洲av片天天在线观看| 久久亚洲精品不卡| 一区二区av电影网| 亚洲午夜精品一区,二区,三区| 国产精品二区激情视频| 18禁裸乳无遮挡动漫免费视频| 日韩一区二区三区影片| 极品教师在线免费播放| 亚洲情色 制服丝袜| 涩涩av久久男人的天堂| 桃红色精品国产亚洲av| 丰满少妇做爰视频| 一本一本久久a久久精品综合妖精| 久久久久视频综合| 久久婷婷成人综合色麻豆| 99国产精品一区二区蜜桃av | 国产精品国产高清国产av | 美女高潮喷水抽搐中文字幕| 亚洲成国产人片在线观看| 12—13女人毛片做爰片一| 精品免费久久久久久久清纯 | 日日摸夜夜添夜夜添小说| 搡老岳熟女国产| 国产国语露脸激情在线看| 欧美日韩国产mv在线观看视频| 国产伦理片在线播放av一区| 亚洲第一欧美日韩一区二区三区 | 欧美精品高潮呻吟av久久| 国产黄频视频在线观看| 国产高清videossex| 999久久久国产精品视频| 男女边摸边吃奶| 亚洲中文日韩欧美视频| 国产不卡av网站在线观看| 在线av久久热| 国产日韩欧美视频二区| 国产av精品麻豆| 啦啦啦视频在线资源免费观看| 亚洲 国产 在线| 一级,二级,三级黄色视频| 亚洲国产欧美网| 天堂动漫精品| 中文字幕色久视频| 成人av一区二区三区在线看| 免费观看人在逋| 国产在线精品亚洲第一网站| 啦啦啦免费观看视频1| 精品一区二区三区四区五区乱码| 午夜两性在线视频| 色94色欧美一区二区| 男女高潮啪啪啪动态图| 乱人伦中国视频| 久久九九热精品免费| 欧美性长视频在线观看| 欧美激情极品国产一区二区三区| 成人影院久久| 欧美av亚洲av综合av国产av| 久久香蕉激情| 一边摸一边抽搐一进一出视频| 国产精品.久久久| 精品第一国产精品| h视频一区二区三区| 无人区码免费观看不卡 | 欧美在线一区亚洲| 欧美日韩国产mv在线观看视频| 国产在线视频一区二区| 欧美日韩亚洲国产一区二区在线观看 | 中文字幕色久视频| 国产深夜福利视频在线观看| 国产av又大| 变态另类成人亚洲欧美熟女 | 搡老乐熟女国产| 啦啦啦免费观看视频1| 久久久精品国产亚洲av高清涩受| 妹子高潮喷水视频| 男女边摸边吃奶| 99国产精品一区二区三区| 成人亚洲精品一区在线观看| 成人18禁高潮啪啪吃奶动态图| 美女午夜性视频免费| 久久精品91无色码中文字幕| 国产老妇伦熟女老妇高清| 久久人妻福利社区极品人妻图片| 窝窝影院91人妻| 欧美黑人欧美精品刺激| 国产片内射在线| 日韩精品免费视频一区二区三区| 日韩大片免费观看网站| 好男人电影高清在线观看| 啪啪无遮挡十八禁网站| 国产亚洲精品一区二区www | 国产福利在线免费观看视频| 成人精品一区二区免费| 日韩三级视频一区二区三区| 美女扒开内裤让男人捅视频| 在线观看一区二区三区激情| 日本vs欧美在线观看视频| 男人舔女人的私密视频| 制服诱惑二区| 在线观看免费视频日本深夜| 18禁黄网站禁片午夜丰满| 日韩免费高清中文字幕av| 国产精品自产拍在线观看55亚洲 | 美国免费a级毛片| 日本黄色日本黄色录像| 成人影院久久| 国产免费av片在线观看野外av| 欧美另类亚洲清纯唯美| 亚洲成人免费av在线播放| 两人在一起打扑克的视频| 日韩制服丝袜自拍偷拍| 一区在线观看完整版| 不卡一级毛片| 成年人黄色毛片网站| 动漫黄色视频在线观看| 久久精品成人免费网站| 欧美国产精品va在线观看不卡| 少妇裸体淫交视频免费看高清 | 国产淫语在线视频| 亚洲综合色网址| 人人妻,人人澡人人爽秒播| 在线观看免费午夜福利视频| 欧美日韩亚洲高清精品| 久久香蕉激情| 欧美日韩亚洲高清精品| 美女主播在线视频| 999久久久精品免费观看国产| 亚洲五月婷婷丁香| 一级黄色大片毛片| 1024视频免费在线观看| 国产精品一区二区精品视频观看| 国产精品 国内视频| 97人妻天天添夜夜摸| a级毛片在线看网站| 国产精品一区二区精品视频观看| av网站免费在线观看视频| 97人妻天天添夜夜摸| 色精品久久人妻99蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 中国美女看黄片| 下体分泌物呈黄色| 亚洲七黄色美女视频| 免费日韩欧美在线观看| 亚洲性夜色夜夜综合| 国产成人精品久久二区二区免费| 色94色欧美一区二区| 亚洲欧美一区二区三区久久| 99国产综合亚洲精品| 国产1区2区3区精品| 国产精品一区二区在线观看99| 亚洲成人手机| 国产一区二区三区在线臀色熟女 | 亚洲一卡2卡3卡4卡5卡精品中文| 免费在线观看影片大全网站| 99riav亚洲国产免费| 热re99久久精品国产66热6| 男人操女人黄网站| 国产一区二区三区视频了| 亚洲一卡2卡3卡4卡5卡精品中文| 男人舔女人的私密视频| 精品一区二区三区四区五区乱码| 老熟女久久久| 亚洲成a人片在线一区二区| 18禁观看日本| 女警被强在线播放| 韩国精品一区二区三区| 亚洲人成伊人成综合网2020| 极品少妇高潮喷水抽搐| 夜夜夜夜夜久久久久| 精品亚洲成国产av| 777米奇影视久久| 在线观看免费视频网站a站| 国产片内射在线| 人人妻人人爽人人添夜夜欢视频| 我的亚洲天堂| 中文字幕人妻丝袜一区二区| 欧美中文综合在线视频| 两个人免费观看高清视频| 夫妻午夜视频| 欧美激情 高清一区二区三区| 黄色a级毛片大全视频| 亚洲第一青青草原| 看免费av毛片| 国产一区二区三区在线臀色熟女 | 日韩有码中文字幕| 久久精品国产a三级三级三级| 99久久精品国产亚洲精品| 涩涩av久久男人的天堂| 色综合欧美亚洲国产小说| 国产免费av片在线观看野外av| 国产aⅴ精品一区二区三区波| 国产午夜精品久久久久久| 色综合婷婷激情| 日韩欧美免费精品| 欧美精品啪啪一区二区三区| 国产精品 国内视频| 欧美成人午夜精品| 久9热在线精品视频| av天堂在线播放| 午夜激情久久久久久久| 久久久精品免费免费高清| 黑人欧美特级aaaaaa片| 国产亚洲欧美在线一区二区| 中文字幕av电影在线播放| 亚洲九九香蕉| 国产老妇伦熟女老妇高清| 亚洲精华国产精华精| 激情视频va一区二区三区| 欧美成狂野欧美在线观看| www.自偷自拍.com| 国产成人免费观看mmmm| 国产精品香港三级国产av潘金莲| 国产亚洲欧美精品永久| 人妻 亚洲 视频| 9色porny在线观看| 国产亚洲欧美在线一区二区| 一本综合久久免费| 国产精品 国内视频| 欧美日韩黄片免| 日韩一卡2卡3卡4卡2021年| 成年动漫av网址| 精品一区二区三区av网在线观看 | 自拍欧美九色日韩亚洲蝌蚪91| 国产免费av片在线观看野外av| 亚洲欧美日韩高清在线视频 | 一本综合久久免费| 欧美精品人与动牲交sv欧美| 超碰97精品在线观看| 亚洲人成伊人成综合网2020| 日韩精品免费视频一区二区三区| 国产一区二区 视频在线| 国产精品一区二区在线观看99| 757午夜福利合集在线观看| 日韩人妻精品一区2区三区| 国产成人欧美| 国产单亲对白刺激| 国产91精品成人一区二区三区 | 成人免费观看视频高清| 老司机深夜福利视频在线观看| 汤姆久久久久久久影院中文字幕| 亚洲全国av大片| 大码成人一级视频| 黄片小视频在线播放| 99久久99久久久精品蜜桃| 老鸭窝网址在线观看| 欧美亚洲日本最大视频资源| 国产av又大| 日韩成人在线观看一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 黄色视频不卡| 国产精品熟女久久久久浪| 久久久久国内视频| 好男人电影高清在线观看| 亚洲avbb在线观看| www.999成人在线观看| 夜夜骑夜夜射夜夜干| 久久精品国产亚洲av高清一级| 精品一区二区三区视频在线观看免费 | 亚洲国产欧美网| 日韩人妻精品一区2区三区| 亚洲全国av大片| 国产男女内射视频| 岛国毛片在线播放|