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

    柱體流激振動能量俘獲理論與技術(shù)研究若干進(jìn)展1)

    2024-04-15 02:52:34徐萬海馬燁璇
    力學(xué)學(xué)報(bào) 2024年3期
    關(guān)鍵詞:振動

    徐萬海 馬燁璇

    * (天津大學(xué)水利工程智能建設(shè)與運(yùn)維全國重點(diǎn)實(shí)驗(yàn)室,天津 300072)

    ? (中國科學(xué)院力學(xué)研究所,北京 100190)

    引言

    大力發(fā)展可再生能源是推動能源降碳減排、實(shí)現(xiàn)“雙碳”目標(biāo)的必由之路.潮流能因其分布廣泛且儲量巨大,是備受關(guān)注的可再生能源之一.相較于波浪能、風(fēng)能等海洋綠色能源,潮流能相對穩(wěn)定且可提前預(yù)報(bào).全世界范圍內(nèi)潮流能的理論估計(jì)值為0.5~1.0 TW 量級[1-2].我國擁有18000 km 的海岸線,可利用的水道超過130 個(gè),擁有豐富的潮流能儲量,理論可開采量約1.4×107kW[3],具備規(guī)?;_發(fā)利用潮流能的條件.

    目前,潮流能轉(zhuǎn)化裝置按照能量俘獲原理可分為渦輪型和非渦輪型.渦輪型潮流能獲取裝置主要通過流動介質(zhì)推動水輪機(jī)槳葉轉(zhuǎn)動,將流體動能轉(zhuǎn)化為獲能裝置的機(jī)械能,再通過電力變化裝置轉(zhuǎn)換成電能.渦輪型潮流能獲取裝置包括:水平軸水輪機(jī)和垂直軸水輪機(jī),其中水平軸水輪機(jī)分為風(fēng)車式、空心貫穿式和導(dǎo)流罩式等,垂直軸水輪機(jī)分為直葉片式和螺旋葉片式等.非渦輪型潮流能獲取裝置采用新型原理進(jìn)行能量轉(zhuǎn)換,主要包括:升力-阻力式、振蕩水翼式、擺動翼式、壓電式、帆式和文丘里式等.相較于非渦輪型潮流能獲取裝置,渦輪型潮流能獲取裝置因其設(shè)備的經(jīng)濟(jì)性與靈活性,在潮流能發(fā)電工程當(dāng)中得到了廣泛的應(yīng)用[4].與世界上潮流能資源最好的英國、挪威等國家相比,我國沿海的潮流總體流速偏低、水深較淺且潮流能分布很不均勻.在東海主要以浙江省沿海潮流能資源最為豐富并且開發(fā)條件最好.而在北方沿海地區(qū),有相當(dāng)多潮流流速較低約在1.0 m/s 左右的海域.由于潮流能的能量密度與流速的立方成正比,因此在潮流流速較低的海域利用水輪機(jī)開發(fā)潮流能勢必需要增大潮流能發(fā)電裝置的尺寸,對水深也提出了更高要求[5].同時(shí),由于水輪機(jī)需要一定的啟動流速,在流速較低時(shí)很難實(shí)現(xiàn)有效的能量轉(zhuǎn)換.這樣不僅造成了潮流能資源的浪費(fèi),也阻礙了潮流能利用技術(shù)在更大范圍內(nèi)推廣.近年來,隨著海洋工程技術(shù)迅猛發(fā)展,研究人員發(fā)現(xiàn)可以利用結(jié)構(gòu)物的流激振動進(jìn)行潮流能發(fā)電.當(dāng)流體流過結(jié)構(gòu)物時(shí),會對結(jié)構(gòu)物表面施加周期性變化的流體力,使其發(fā)生往復(fù)運(yùn)動,而結(jié)構(gòu)物的往復(fù)運(yùn)動又會改變流體的流動模式,進(jìn)而改變作用于結(jié)構(gòu)物表面的流體力,這種流固耦合現(xiàn)象稱為流激振動(flow-induced vibration,FIV)[6-10].流激振動是造成海洋工程結(jié)構(gòu)物疲勞損傷的主要原因之一.從另外方面去看,流激振動具有自激且自持的特性,非常適合用于潮流能的開發(fā)利用[11-13].

    流激振動能量俘獲裝置的研制始于20 世紀(jì)90 年代,按照能量轉(zhuǎn)化方式可分為壓電式[14-16]和電磁式兩種類型[17-19].壓電式流激振動發(fā)電裝置是利用流體流過壓電材料時(shí)使其發(fā)生變形而產(chǎn)生的電壓變化實(shí)現(xiàn)能量轉(zhuǎn)換.比較經(jīng)典的壓電式流激振動發(fā)電裝置包括懸臂梁式[20]、鰻魚式[21-22]和樹式[23-25]等.懸臂梁式能量俘獲裝置利用壓電梁材料流激振動作用發(fā)電,可應(yīng)用于微型、小型傳感器等電子設(shè)備的供電系統(tǒng),但整體電能輸出量低,僅在μW 或mW 量級.鰻魚式能量俘獲裝置利用柔性壓電材料的擺動或扭轉(zhuǎn)發(fā)電,最初被應(yīng)用于潮流能利用[21-22,26],經(jīng)改進(jìn)后進(jìn)一步應(yīng)用于風(fēng)能利用.樹式能量俘獲裝置利用了流體旋渦脫落與結(jié)構(gòu)的共振作用,能量轉(zhuǎn)化能力相對較強(qiáng),該類裝置形式多樣,改進(jìn)型裝置層出不窮[23-24,27].整體而言,壓電式能量俘獲裝置的能量轉(zhuǎn)換能力較弱,現(xiàn)階段難以實(shí)現(xiàn)大規(guī)模開發(fā).

    電磁式能量俘獲裝置通過振子流激振動帶動永磁體切割磁場,從而實(shí)現(xiàn)發(fā)電.Wang 等[17]設(shè)計(jì)了一種電磁式能量俘獲裝置,利用旋渦脫落使頂梁結(jié)構(gòu)形變,從而推動永磁體切割磁感線進(jìn)行發(fā)電,不過該裝置的能量俘獲功率較低,僅為μW 量級.Vortex Bladeless 公司借鑒了樹式能量俘獲裝置的理念,設(shè)計(jì)了電磁式流激振動能量俘獲裝置.該裝置內(nèi)部不存在齒輪、螺栓或其他機(jī)械的運(yùn)動部件,成本極大降低、發(fā)電效率顯著提升,但僅停留在樣機(jī)研究階段,能量轉(zhuǎn)換能力是否能達(dá)到原型機(jī)的效果仍有待進(jìn)一步驗(yàn)證[13].2005 年,美國密歇根大學(xué)的Bernitsas教授團(tuán)隊(duì)研制了VIVACE (vortex-induced vibration aquatic clean energy)裝置(如圖1 所示),通過柱體流激振動將流體能量轉(zhuǎn)換為機(jī)械能,再利用發(fā)電機(jī)將機(jī)械能轉(zhuǎn)換為電能[18-19].該裝置中柱體振子受到的流體力遠(yuǎn)大于壓電式裝置振子受到的流體力,對阻尼的克服強(qiáng)度遠(yuǎn)高于壓電式裝置,能量轉(zhuǎn)換能力顯著增強(qiáng),也有利于大規(guī)模開發(fā),成為電磁式流激振動能量俘獲研究的突破點(diǎn)和里程碑.

    圖1 流激振動能量俘獲裝置示意圖[18]Fig.1 Schematic of energy harvesting device based on FIV[18]

    相較于波浪能俘獲裝置和其他潮流能俘獲裝置,VIVACE 這種基于柱體結(jié)構(gòu)流激振動的能量俘獲裝置具有較高的能量密度[5](如圖2 所示),極具工程應(yīng)用前景.眾多學(xué)者針對單個(gè)柱體結(jié)構(gòu)、柱群結(jié)構(gòu)的流激振動能量俘獲理論和技術(shù)開展了大量研究,重點(diǎn)圍繞雷諾數(shù)條件、結(jié)構(gòu)截面形狀、系統(tǒng)阻尼特性和多柱體排布方式等多種影響因素,探討了流激振動能量轉(zhuǎn)換效率.本文將對國內(nèi)外流激振動能量俘獲理論和技術(shù)方面的研究進(jìn)展進(jìn)行全方位綜述,重點(diǎn)對現(xiàn)有流激振動能量轉(zhuǎn)化裝置的響應(yīng)機(jī)理、能量俘獲功率和能量轉(zhuǎn)換效率的影響因素進(jìn)行總結(jié)闡述,并對未來的研究工作進(jìn)行展望,以期促進(jìn)流激振動能量俘獲理論的發(fā)展和流激振動能量轉(zhuǎn)換裝置的工程應(yīng)用.

    圖2 VIVACE 與波浪能及潮流能發(fā)電裝置的能量密度[18]Fig.2 The energy density of VIVACE,wave and other tidal energy harvesting devices[18]

    1 流激振動能量俘獲評估理論模型

    流激振動能量俘獲評估理論模型可用于計(jì)算流激振動能量俘獲功率和轉(zhuǎn)換效率,目前流激振動能量俘獲理論模型的推導(dǎo)多以Bernitsas 等[18]提出的數(shù)學(xué)模型為參考.

    假設(shè)來流沿x方向,柱體振子在y方向振動.那么柱體振子的動力學(xué)響應(yīng)可以由如下二階線性微分方程表示

    其中mosc為振蕩系統(tǒng)的質(zhì)量,ctotal為總體阻尼系數(shù),K為系統(tǒng)剛度系數(shù),y為柱體振子的振動位移,符號(?)表示對時(shí)間t求導(dǎo).FFluid為柱體在振動方向上受到的流體力,可分解為兩部分,即黏性力部分Fviscous和非黏性力部分Finciscid

    黏性力Fviscous和非黏性力Finciscid可分別由式(3)和式(4)表示

    其中cy(t)為柱體振子振動方向的流體力系數(shù),ρ為流體密度,U為來流速度,D為柱體截面的特征尺寸,L為柱體振子的長度,ma為柱體振子在流體中的附加質(zhì)量.將式(2)~式(4) 代入式(1),進(jìn)一步化簡可得

    在一個(gè)振動周期Tosc內(nèi),流體對柱體振子系統(tǒng)所做的功,即柱體振子可轉(zhuǎn)化的機(jī)械能功率可表示為

    若柱體流激振動過程中可近似簡化為簡諧振動,式(6)中最后一項(xiàng)中唯一非0 項(xiàng)為流激振動能量的轉(zhuǎn)化過程為:先由流體動能轉(zhuǎn)換為柱體系統(tǒng)機(jī)械能,再由系統(tǒng)機(jī)械能轉(zhuǎn)化為電能.機(jī)械能向電能轉(zhuǎn)化的過程中可引入獲能阻尼系數(shù)charness,獲能阻尼和系統(tǒng)損耗阻尼共同構(gòu)成總體阻尼

    其中csystem為系統(tǒng)損耗阻尼系數(shù).柱體流激振動俘獲的能量輸出功率可表示為

    能量耗散功率可表示為

    柱體振子掃過區(qū)域內(nèi)流體所蘊(yùn)含的能量可表述為

    式中Amax為柱體振子振動過程中的最大振幅.柱體流激振動的能量轉(zhuǎn)換效率可表示為

    式中BetzLimit為貝茨極限,是理想狀態(tài)下流體動能轉(zhuǎn)換成機(jī)械能的極限比值,其值為16/27 (約為59.26%)[28].可基于式(11)計(jì)算柱體流激振動能量轉(zhuǎn)換效率的理論值,對流激振動能量俘獲裝置的效果進(jìn)行評估.

    2 單柱結(jié)構(gòu)流激振動能量俘獲

    根據(jù)流激振動能量俘獲評估理論模型可知,單個(gè)柱體流激振動能量轉(zhuǎn)換效率主要與振幅和頻率相關(guān),而柱體的振幅和頻率受到雷諾數(shù)、質(zhì)量比、阻尼比和截面形狀等多種因素影響.學(xué)者們重點(diǎn)關(guān)注了圓柱和非圓截面柱體在不同雷諾數(shù)、質(zhì)量比、阻尼比等條件下的振動響應(yīng)特性和能量俘獲性能.

    2.1 單圓柱流激振動能量俘獲

    早期VIVACE 能量俘獲裝置的柱體振子設(shè)計(jì)為單個(gè)光滑圓柱結(jié)構(gòu),基于單個(gè)圓柱結(jié)構(gòu)流激振動響應(yīng)機(jī)理推導(dǎo)了能量轉(zhuǎn)換效率理論模型,并通過模型實(shí)驗(yàn)初步證明了VIVACE 能量俘獲裝置具有較高的能量轉(zhuǎn)換效率.隨后,Bernitsas 等[18-19]研究了雷諾數(shù)、阻尼比和邊界條件等多種因素對單根圓柱流激振動能量轉(zhuǎn)換效率的影響.

    光滑圓柱的流激振動幾乎發(fā)生在整個(gè)雷諾數(shù)范圍內(nèi),僅在3 個(gè)雷諾數(shù)范圍內(nèi)不發(fā)生振動,即Re<40 (邊界層不分離和形成一對穩(wěn)定旋渦的區(qū)域),150

    考慮到VIVACE 能量俘獲裝置的實(shí)際應(yīng)用場景,自由液面和水底邊壁可能會對柱體振子流激振動響應(yīng)和能量俘獲功率產(chǎn)生影響.VIVACE 裝置在淺水區(qū)域(河流或溪流)中作業(yè)時(shí),會受到自由表面的不利影響.柱體靠近自由液面時(shí),流激振動的共振范圍減小、頻率鎖定特性產(chǎn)生顯著改變[33].主要原因在于:圓柱振動對自由液面產(chǎn)生強(qiáng)烈調(diào)制作用,改變了柱體頂部的渦量和渦脫落周期,進(jìn)而導(dǎo)致在渦脫落頻率、振動頻率和固有頻率之間的鎖頻帶寬減小,流激振動難以維持大范圍的共振.在自由液面的影響下柱體振幅降低,且在振動啟動階段并未像遠(yuǎn)離自由液面那樣跳躍增長,而是隨流速增大緩慢增長.由于圓柱尾流的模式切換,柱體流激振動共振開始和結(jié)束階段均存在滯回現(xiàn)象[31].VIVACE 裝置靠近水底邊壁時(shí),邊壁效應(yīng)會對柱體振子的響應(yīng)特性產(chǎn)生影響.在雷諾數(shù)Re=8000~150000 范圍內(nèi),邊壁效應(yīng)的影響隨間距比發(fā)生變化.間距比小于0.65時(shí),流激振動響應(yīng)的共振范圍向高流速區(qū)域轉(zhuǎn)移,并在共振結(jié)束階段出現(xiàn)滯回現(xiàn)象,無量綱正向振幅因圓柱后僅出現(xiàn)頂部尾渦可增大至2.0.間距比在0.65~3.0 之間時(shí),流激振動響應(yīng)的共振范圍未發(fā)生顯著變化,正向振幅和負(fù)向振幅整體對稱,但振幅因尾渦強(qiáng)度的減弱而減小.間距比大于3.0 后,邊壁效應(yīng)的影響可忽略[31].

    圓柱表面粗糙度通過影響邊界層分離點(diǎn)、湍流水平、尾渦強(qiáng)度,進(jìn)而對流激振動響應(yīng)和流固耦合作用產(chǎn)生重要影響[34-36].為了進(jìn)一步提高能量轉(zhuǎn)換效率,可將由直條粗糙帶構(gòu)成的被動湍流控制器(passive turbulence control,PTC)安裝于VIVACE 系統(tǒng)中圓柱振子的外表面(如圖3 所示).PTC 對VIVACE裝置產(chǎn)生3 種有利作用:一是通過局部增大圓柱的表面粗糙度,使邊界層由層流轉(zhuǎn)變?yōu)橥牧?二是有助于控制旋渦的脫落,使圓柱展向不同位置的流動參數(shù)同步,增大尾流的展向相關(guān)長度;三是將流動的分離點(diǎn)固定,促使圓柱發(fā)生馳振[37-39].整體而言,PTC觸發(fā)了流動分離并給邊界層充能,從而顯著提升尾流強(qiáng)度,產(chǎn)生更高水平的升力.安裝PTC 后的圓柱能夠在光滑圓柱不發(fā)生振動的臨界雷諾數(shù)范圍內(nèi)產(chǎn)生大幅振動.PTC 的最佳安裝角度范圍為57°~80°,在這一范圍內(nèi)圓柱流激振動的幅值顯著增大,最高無量綱振幅可達(dá)2.7.PTC 使旋渦脫落的軸向相關(guān)性增強(qiáng),能夠增大流激振動共振范圍,并提前進(jìn)入共振階段[40].這些影響能夠最大限度地將流體能量轉(zhuǎn)化為機(jī)械能,從而增強(qiáng)裝置的能量俘獲功率,見圖4.徐萬海等[41-42]進(jìn)一步關(guān)注了自由液面對安裝PTC 后圓柱流激振動的影響.近自由液面PTC 圓柱隨流速增大可觀測到3 個(gè)區(qū)域,即渦激區(qū)、渦激向馳振轉(zhuǎn)化區(qū)和馳振區(qū).在渦激區(qū),自由液面的影響不明顯;在馳振區(qū),自由液面效應(yīng)使振動幅值和流體力系數(shù)出現(xiàn)顯著差異.

    圖3 安裝PTC 的圓柱模型[35-36]Fig.3 Cylinder model with PTC[35-36]

    圖4 近自由液面PTC 圓柱的能量俘獲功率[42]Fig.4 Harnessed power of a single cylinder with PTC near the free surface[42]

    為了增強(qiáng)VIVACE 裝置的性能,Lee 等[43-44]開發(fā)了第一代虛擬阻尼-彈簧系統(tǒng) (virtual damper–spring system,VCK)替代VIVACE 裝置中的實(shí)體彈簧.VCK 系統(tǒng)由圓柱振子、同步齒帶和滑輪傳動、電動機(jī)/發(fā)電機(jī)和控制器組成.控制器通過測量的速度和位移反饋向系統(tǒng)提供阻尼力和回復(fù)力.因此,在此過程中不會使力與位移之間產(chǎn)生人為的相位差,導(dǎo)致能量轉(zhuǎn)換的偏移.第一代VCK 系統(tǒng)采用數(shù)據(jù)采集系統(tǒng)讀取圓柱的位置并完成力反饋,成功提升了VIVACE 裝置的運(yùn)行效率.但這套數(shù)據(jù)采集系統(tǒng)中的ADC (analog digital conversion) 模塊和DAC(digital analog conversion)模塊需要進(jìn)行數(shù)字信號濾波,導(dǎo)致力反饋和位移之間出現(xiàn)相位差,造成一定的遲滯效應(yīng).Sun 等[45]在第一代VCK 系統(tǒng)基礎(chǔ)上研制了第二代VCK 系統(tǒng).第二代VCK 系統(tǒng)采用數(shù)字信號,使得系統(tǒng)響應(yīng)速度更快,整個(gè)系統(tǒng)只在中央處理器運(yùn)行部分存在10 μs 的延遲,實(shí)現(xiàn)了對位移和速度的實(shí)時(shí)測量,如圖5 所示.

    圖5 虛擬阻尼-彈簧系統(tǒng)[45]Fig.5 Virtual damper-spring system[45]

    目前,基于單個(gè)圓柱振子的VIVACE 裝置初步完成了理論和技術(shù)積累,已申請多項(xiàng)專利[46-48],并多次在美國的圣克萊爾河及荷蘭的多條運(yùn)河中進(jìn)行了現(xiàn)場樣機(jī)測試.

    2.2 非圓截面單柱流激振動能量俘獲

    對于不安裝PTC 湍流控制器的圓柱結(jié)構(gòu),通常在來流作用下主要激發(fā)渦激振動,渦激振動具有自限性,能量俘獲有限.非圓截面柱體結(jié)構(gòu),特別是具有尖角截面的柱體,通常存在渦激振動和馳振兩種模式.馳振為非自限制振動,隨著流速增大振幅逐漸增大,有利于能量俘獲.因此,非圓截面柱體結(jié)構(gòu)(如圖6 所示)流激振動的響應(yīng)特性和能量俘獲亦引起了諸多關(guān)注.

    圖6 非圓截面柱體示意圖[49-55]Fig.6 Schematic of the cylinder with non-circular cross-section[49-55]

    學(xué)術(shù)界對正方形截面柱體結(jié)構(gòu)的流激振動關(guān)注較早,并已初步揭示了其馳振特征[56-60].Tamimi 等[49]實(shí)驗(yàn)研究了正方形截面柱體的流激振動特性和能量俘獲能力.實(shí)驗(yàn)觀測的雷諾數(shù)范圍為Re=2300~40000,來流攻角工況為0°和45°.攻角為0°時(shí),正方形截面柱體表現(xiàn)為典型的馳振特征,振幅隨流速持續(xù)增長,升力系數(shù)和響應(yīng)頻率隨流速保持緩慢變化,柱體獲取的機(jī)械能隨流速持續(xù)增長,

    在實(shí)驗(yàn)觀測的流速范圍內(nèi)并未出現(xiàn)極值;攻角為45°時(shí),柱體表現(xiàn)為典型的渦激振動特征,振幅隨流速先增大后減小,且同步出現(xiàn)鎖頻現(xiàn)象和相位角跳變現(xiàn)象,柱體獲取的機(jī)械能先增大后減小,在共振范圍內(nèi)存在極值.雖然渦激振動時(shí)的獲能效率高于馳振,但渦激振動獲取的機(jī)械能存在極值,而發(fā)生馳振時(shí)獲取的機(jī)械能會隨流速增大而持續(xù)增長.Han等[61]在更低的質(zhì)量比條件下(m*=1.45)開展了正方形截面柱體流激振動的數(shù)值仿真分析.攻角為45°時(shí),低質(zhì)量比柱體的渦激作用顯著增強(qiáng),隨著流速增大出現(xiàn)了持續(xù)鎖頻現(xiàn)象;攻角為22.5°時(shí),正方形截面柱體振動中存在渦激作用的高頻成分和馳振的低頻成分,高流速工況下的大振幅是由渦激作用和馳振共同導(dǎo)致的.由于持續(xù)鎖頻,攻角為45°條件下的低質(zhì)量比柱體具備更高的能量俘獲能力,且在阻尼比小于0.3 時(shí)能量轉(zhuǎn)換效率對阻尼比變化不敏感,而阻尼比為0.7 時(shí)能量轉(zhuǎn)換效率最高.Zhang 等[62]通過數(shù)值仿真在更高的雷諾數(shù)范圍內(nèi)(Re=15500~232000)研究了攻角為0°時(shí)正方形截面柱體的流激振動響應(yīng)和能量俘獲能力.隨著雷諾數(shù)增大,柱體振動響應(yīng)依次出現(xiàn)渦激振動初始分支區(qū)(Re<46000)、渦激振動上分支區(qū)(46000 159000).渦激振動上分支區(qū)能量轉(zhuǎn)換效率出現(xiàn)最大值,而在馳振區(qū)能量轉(zhuǎn)換效率降低,但能量俘獲功率隨流速增大持續(xù)增長.柱體質(zhì)量比和阻尼比的影響顯著,阻尼比小于0.6 時(shí),隨著阻尼比的增大,柱體振幅減小,但能量俘獲功率和轉(zhuǎn)換效率增加;阻尼比達(dá)到0.8 時(shí),馳振響應(yīng)消失.增加質(zhì)量比會縮短渦激振動響應(yīng)區(qū)和渦激振動向馳振過渡區(qū),柱體在更低的雷諾數(shù)條件下發(fā)生馳振,且能量俘獲功率和能量轉(zhuǎn)換效率隨質(zhì)量比的增加而降低(如圖7 所示).

    圖7 不同質(zhì)量比下方柱的振動幅值[62]Fig.7 Vibration amplitude of a square cylinder with different mass ratio[62]

    對于正方形截面柱體,渦激振動的能量轉(zhuǎn)換效率高于馳振的能量轉(zhuǎn)換效率,主要原因是流體所蘊(yùn)含的能量與流速呈三次方關(guān)系,能量俘獲功率在由渦激振動向馳振轉(zhuǎn)變的增長量小于流體蘊(yùn)含能量的增長幅度.但由于馳振的不穩(wěn)定性,能量俘獲功率隨流速增大持續(xù)增長而沒有極限,表明利用馳振進(jìn)行流激振動能量俘獲的潛力巨大[63].Zhang 等[64]進(jìn)一步研究了長方形截面柱體流激振動能量俘獲能力.邊長比(流動方向邊長與垂直來流方向邊長之比)對柱體流激振動響應(yīng)和能量俘獲功率存在顯著影響.當(dāng)邊長比從1/6 增加到1.5 時(shí),柱體振動逐漸受到抑制,當(dāng)邊長比為2.0 時(shí),在任何流速下都不會出現(xiàn)渦激振動和馳振.對于能量收集,當(dāng)邊長比從2.0 減小到1/4 時(shí),能量俘獲功率呈增大趨勢,最大能量轉(zhuǎn)換效率可同步提高至15.5%;當(dāng)邊長比減小到1/6 時(shí),能量俘獲功率保持增加趨勢,但最大能量轉(zhuǎn)換效率并未進(jìn)一步提高.

    三角形截面柱體在外部流場作用下也會產(chǎn)生大振幅馳振[50,65-66].對于正三角形截面柱體,當(dāng)來流方向與其某一邊長垂直(來流攻角為60°)時(shí),渦激振動響應(yīng)和馳振會隨流速增大依次激發(fā),并且在渦激振動向馳振的轉(zhuǎn)化區(qū)存在響應(yīng)頻率降低的現(xiàn)象.隨著系統(tǒng)阻尼比增大,正三角形截面柱體會發(fā)生軟馳振和硬馳振.軟馳振是指柱體由渦激振動轉(zhuǎn)變而成的馳振,具有自激性;硬馳振是指柱體不可通過自激勵(lì)由渦激振動轉(zhuǎn)變?yōu)轳Y振,而是需要較大的外部激勵(lì)產(chǎn)生馳振[67].系統(tǒng)存在臨界阻尼比,超過臨界阻尼比后柱體僅發(fā)生渦激振動,系統(tǒng)的臨界阻尼比為柱體能量俘獲的最優(yōu)阻尼比.阻尼比超過臨界阻尼比時(shí),正三角形截面柱體的能量轉(zhuǎn)換效率在渦激振動階段對阻尼比的變化不敏感,在馳振階段隨阻尼比的增大而升高[68].Zou 等[69]研究了來流攻角對正三角形截面柱體流激振動響應(yīng)和能量俘獲性能的影響.來流攻角在0°~30°之間時(shí),柱體流激振動幅值較小,能量俘獲功率近乎為0;來流攻角在30°~60°之間時(shí),能量俘獲功率隨流速增大單調(diào)增長.攻角為60°時(shí),能量俘獲功率最高,此時(shí)無量綱振幅可高達(dá)3.0.通過分析柱體的尾渦模式可以發(fā)現(xiàn),攻角較高的條件下旋渦脫落模式由2S 型(2 個(gè)獨(dú)立旋渦)過渡到P+S 型(1 對旋渦+1 個(gè)獨(dú)立旋渦),導(dǎo)致作用于柱體的非定常流體載荷增強(qiáng),從而提高了柱體振幅和能量俘獲功率.Yan 等[70]基于正三角形截面柱體流激振動能量俘獲應(yīng)用場景,考慮發(fā)電機(jī)磁通密度的影響建立能量俘獲模型,研究了能量轉(zhuǎn)換效率和能量俘獲能力.系統(tǒng)阻尼隨磁通密度的增加而增大,正三角形截面柱體的振動響應(yīng)由軟馳振轉(zhuǎn)變?yōu)橛柴Y振.通過合理的控制策略調(diào)整發(fā)電機(jī)的勵(lì)磁電壓,進(jìn)而控制磁通密度,能夠使柱體保持較高的能量俘獲能力.為了進(jìn)一步提高三角形截面柱體流激振動能量俘獲能力,Shao 等[71]研究了不同邊高比(等腰三角形底邊與底邊高之比)的等腰三角形截面柱體流激振動響應(yīng)特性.邊高比低于1.0 時(shí),柱體產(chǎn)生渦激振動,僅在外部激勵(lì)下發(fā)生硬馳振,在硬馳振過程中能量俘獲功率較高;邊高比為1.0 和1.5 時(shí),隨著流速增大柱體由渦激振動轉(zhuǎn)變?yōu)轳Y振,即產(chǎn)生軟馳振,能量俘獲功率隨流速持續(xù)增長,但軟馳振過程中的能量俘獲功率低于硬馳振過程.

    除了圓形、三角形和方形截面柱體,梯形、菱形、六邊形、八邊形等其他截面形式柱體的流激振動特性也獲得了一定關(guān)注,并對不同截面形式柱體流激振動的能量俘獲能力進(jìn)行了綜合對比[51-52,72].Ding 等[51]通過數(shù)值仿真分析了PTC 圓柱、正三角形截面柱體、正方形截面柱體、準(zhǔn)梯形截面柱體流激振動的能量轉(zhuǎn)換效率.準(zhǔn)梯形截面柱體長邊迎流時(shí),隨流速增大柱體依次產(chǎn)生渦激振動和馳振;短邊迎流時(shí),柱體出現(xiàn)低幅高頻振動響應(yīng).對比能量轉(zhuǎn)換效率結(jié)果,長邊迎流的準(zhǔn)梯形截面柱體略高于PTC 圓柱,顯著高于正三角形截面柱體和方形截面柱體.李恒[72]實(shí)驗(yàn)觀測了梯形和菱形截面柱體的流激振動特性.當(dāng)梯形截面柱體長邊迎流時(shí),柱體能夠產(chǎn)生渦激振動和馳振,其馳振階段的振幅整體略高于正三角形截面柱體的振幅.菱形截面柱體的流激振動特征與攻角為45°時(shí)的正方形截面柱體的振動特征類似,但其振幅存在較大差異.若菱形的迎流向?qū)蔷€長于垂直流向?qū)蔷€時(shí),菱形截面柱體振幅高于正方形截面柱體的振幅,反之振幅結(jié)果降低.Zhang 等[52]對6 種不同截面形狀(三角形、正方形、六邊形、八邊形、二十四邊形和圓形)的柱體流激振動能量俘獲能力進(jìn)行了對比(如圖8 所示).三角形、正方形、六邊形和八邊形截面柱體在來流作用下先發(fā)生渦激振動,隨著流速增大后發(fā)生馳振;而二十四邊形和圓形截面柱體僅發(fā)生渦激振動,隨著流速增大柱體振幅保持在極低水平.在渦激振動區(qū)域,圓柱的流激振動振幅和能量轉(zhuǎn)換效率最高;在馳振區(qū)域,柱體流激振動振幅和能量俘獲功率隨橫截面形狀邊數(shù)減小而增大,即三角形 >正方形 >六邊形 >八邊形.從柱體的尾流模式可以進(jìn)一步闡述不同截面形式產(chǎn)生不同獲能效果的機(jī)理.在低流速的渦激振動區(qū)域,正方形、六邊形、八邊形和圓形截面柱體的尾渦模式為2S 型,而三角形截面柱體后的尾渦模式為2P 型,導(dǎo)致三角形截面柱體受到更強(qiáng)的流體載荷作用,從而提高能量俘獲功率.在高流速的馳振區(qū)域,三角形截面柱體后的尾渦模式為2S +2P 型,相較于其他截面形式的柱體,三角形截面柱體后的尾渦密度更高、強(qiáng)度更大,因此獲能功率和能量轉(zhuǎn)換效率最高.

    圖8 不同非圓截面柱體的能量俘獲功率[52]Fig.8 Harnessed power of a single cylinder with non-circular crosssection[52]

    通過圓柱流激振動抑制研究發(fā)現(xiàn),柱體結(jié)構(gòu)通過安裝附屬裝置或改變截面形狀后并非總是能有效抑制振動,某些工況下反而增強(qiáng)振動,這為柱體流激振動能量俘獲提供了新思路.Wang 等[53]發(fā)現(xiàn)在圓柱迎流面開不同形狀的凹槽會對其流激振動特性產(chǎn)生顯著影響.對于方形凹槽圓柱,開槽角度為30°,60°和150°時(shí)柱體振幅增大、共振范圍增加、能量俘獲功率提升;開槽角度為90°時(shí),柱體振動被抑制,共振范圍降低50%.對于三角形和圓形凹槽圓柱,開槽角度為90°時(shí),柱體振幅減小;其他開槽角度下共振范圍增大、能量俘獲功率略有提升.對于丁字形凹槽圓柱,開槽角度為120°時(shí),柱體振幅增大、能量俘獲功率提升;開槽角度為90°時(shí),柱體振幅減小、能量俘獲功率降低;其他開槽角度下柱體振幅和能量俘獲功率與光滑圓柱結(jié)果基本一致.Zhu 等[54]基于仿生原理提出在圓柱外表面安裝鰭狀條帶的方式提高流激振動能量俘獲能力.鰭狀條帶的安裝角度為20°和45°時(shí),柱體產(chǎn)生I 型硬馳振,能量轉(zhuǎn)換效率和俘獲功率均顯著提升;安裝角度為0°和60°時(shí),柱體產(chǎn)生II 型硬馳振,能量轉(zhuǎn)換效率和俘獲功率提升幅度較小;安裝角度為90°時(shí),柱體僅發(fā)生渦激振動;安裝角度為120°時(shí),柱體振動被抑制.隨后,Zhu 等[55]發(fā)明了一種安裝自由旋轉(zhuǎn)五角星葉輪的圓柱振子,能夠同時(shí)俘獲柱體振動能量和葉輪旋轉(zhuǎn)能量.葉輪的旋轉(zhuǎn)方向不斷變化,會導(dǎo)致流體流動不穩(wěn)定和水動力不穩(wěn)定,能夠使柱體振子尾流寬度增加,并產(chǎn)生不穩(wěn)定的旋渦脫落,從而增強(qiáng)柱體振子振幅和能量俘獲功率.Barati 等[73]通過數(shù)值仿真研究了弓形(圓弧與其對應(yīng)的弦圍成的圖形)截面柱體流激振動能量俘獲能力.柱體弧面迎流,關(guān)注了5 種長寬比(柱體沿流向長度與柱體垂直流向?qū)挾戎?的影響.隨著長寬比的增大,柱體振子的旋渦脫落頻率降低,但振幅和能量俘獲功率增大.Shao 等[74]和燕翔等[75]實(shí)驗(yàn)研究了T 形截面柱體振子的流激振動特性與能量俘獲能力.柱體T 形截面的頂部迎流,關(guān)注了5 種長寬比(柱體沿流向長度與柱體垂直流向?qū)挾戎?的影響.隨著截面長寬比的增大,振動響應(yīng)逐漸由硬馳振轉(zhuǎn)化為軟馳振,且馳振階段的振幅與頻率均有所降低,對應(yīng)的能量轉(zhuǎn)換效率減小.Lian 等[76]和冉聃頡等[77]將圓形和多邊形組合形成了多種復(fù)雜截面形狀,并實(shí)驗(yàn)研究了復(fù)雜截面形狀柱體的能量俘獲能力.弓形、三角形和T 形組合而成的復(fù)雜截面形狀柱體較正三角形截面柱體具有更高的能量俘獲能力.從眾多復(fù)雜截面形狀柱體流激振動能量俘獲結(jié)果總結(jié)發(fā)現(xiàn),當(dāng)柱體截面迎流面具有對稱尖銳凸起、被流面無旋渦再附著的特點(diǎn)時(shí),有利于柱體振子在大阻尼下發(fā)生自激馳振,更有利于能量俘獲.

    3 柱群流激振動能量俘獲

    對于大規(guī)模的流激振動能量開發(fā),各柱體振子之間不可避免地存在流場干涉.在尾流作用下,柱體結(jié)構(gòu)的振動可能增強(qiáng)[78-83],有利于流激振動能量俘獲.

    3.1 圓柱群流激振動能量俘獲

    雙圓柱系統(tǒng)通常被認(rèn)為是柱群系統(tǒng)的基本組成單元,具有并列、串列和交錯(cuò)3 種排布方式[84].對于串列和交錯(cuò)排布的雙圓柱,在特定間距比范圍內(nèi)因尾流干涉而產(chǎn)生大幅振動.低雷諾數(shù)下的光滑雙圓柱系統(tǒng)流激振動的能量轉(zhuǎn)換效率偏低[85],仍然有必要為圓柱振子安裝PTC 裝置以提高能量俘獲能力.Sun 等[86]通過模型實(shí)驗(yàn)研究了間距和剛度對串列排布帶有PTC 雙圓柱振子VIVACE 裝置的流體能量俘獲效果.對于所關(guān)注的3 種間距比(1.57,2.0 和2.57),雙圓柱的能量俘獲功率是單圓柱的2.56~13.49 倍;能量轉(zhuǎn)換效率為單圓柱的2.0~6.68 倍.受到下游圓柱的影響,上游圓柱的能量轉(zhuǎn)換功率增長多達(dá)100%.下游圓柱的能量轉(zhuǎn)換功率因干涉作用的影響在不同流速下發(fā)生顯著變化.串列排布雙圓柱振子的VIVACE裝置能夠利用流速低至0.4 m/s 的流動能量,并且無流速上限,最大能量俘獲功率出現(xiàn)在間距較小、彈簧剛度較大時(shí).Xu 等[87]實(shí)驗(yàn)研究了阻尼比和剛度對交錯(cuò)排布帶有PTC 雙圓柱振子VIVACE 裝置的流體能量俘獲效果.雙圓柱在垂直來流方向的間距比為1.0,在沿來流方向的間距比為2.57.雙圓柱的能量俘獲功率隨阻尼比的增大而升高,且高于兩個(gè)單體圓柱流激振動的能量俘獲功率(如圖9 所示).較低的剛度有利于雙圓柱的啟動,能夠在較低的流速下實(shí)現(xiàn)流體能量俘獲.Ding 等[88]通過數(shù)值仿真分析了下游圓柱流激振動能量俘獲功率低于上游圓柱的原因.在Re為60000 的渦激振動響應(yīng)區(qū)域,下游圓柱在小間距比為2.0 工況下受上游圓柱遮蔽作用的影響,振動被顯著抑制,能量俘獲功率和轉(zhuǎn)換效率降低;在Re為110000 的馳振響應(yīng)區(qū)域,由上游圓柱脫落的旋渦在不同的時(shí)機(jī)沖擊下游圓柱,導(dǎo)致下游圓柱的能量俘獲功率和轉(zhuǎn)換效率低于上游圓柱的結(jié)果.通過尾流場發(fā)現(xiàn),在低流速的渦激振動區(qū)域,上游圓柱和下游圓柱之間的間隙流方向與下游圓柱振動方向相反,下游圓柱的振動響應(yīng)被抑制,能量轉(zhuǎn)換功率降低;在高流速的馳振區(qū)域,間隙流抵達(dá)下游圓柱的下側(cè)區(qū)域,推動下游圓柱朝其運(yùn)動方向移動,進(jìn)一步增大振幅,提升能量轉(zhuǎn)換功率[86].通過分析柱體系統(tǒng)參數(shù)影響規(guī)律發(fā)現(xiàn),在圓柱發(fā)生渦激振動區(qū)域,隨著系統(tǒng)剛度和固有頻率的增加,遮蔽效應(yīng)大幅減弱,但若下游圓柱上、下兩側(cè)脫落的旋渦由于上游圓柱脫落的旋渦吸引而積累不夠時(shí),會導(dǎo)致下游圓柱振動受到一定程度的抑制.在圓柱發(fā)生馳振區(qū)域,遮蔽效應(yīng)的強(qiáng)弱取決于下游圓柱附近的渦量是被上游圓柱剪切層產(chǎn)生的渦流增強(qiáng)還是減弱.

    圖9 交錯(cuò)排布圓柱的能量俘獲功率[87]Fig.9 Harnessed power of two staggered cylinders[87]

    為了提高下游圓柱流激振動能量俘獲能力,研究上游圓柱對下游圓柱的遮蔽效應(yīng)十分必要.Sun等[89]系統(tǒng)地總結(jié)了遮蔽效應(yīng)的產(chǎn)生原因及規(guī)避措施.下游圓柱受到的遮蔽效應(yīng)與阻尼比密切相關(guān),隨著阻尼比的增加,遮蔽效應(yīng)減弱,下游圓柱的振動恢復(fù)較為穩(wěn)定的幅值和頻率.為了改善下游圓柱的遮蔽效應(yīng),Sun 等[90]系統(tǒng)地研究了自適應(yīng)阻尼系數(shù)下串列排布雙圓柱VIVACE 裝置的流激振動響應(yīng)特性及能量俘獲性能.通過使用與速度成正比的自適應(yīng)阻尼系數(shù),VIVACE 裝置的能量俘獲能力有所提升.對于串列排布的雙圓柱振子,自適應(yīng)阻尼能夠減小上游圓柱對下游圓柱的遮蔽效應(yīng),使下游圓柱的能量俘獲功率翻倍.自適應(yīng)阻尼能降低VIVACE 裝置的啟動流速,提高其高效運(yùn)行的流速范圍,并提升能力輸出的穩(wěn)定性,使能量俘獲功率在渦激振動和馳振的過渡區(qū)不產(chǎn)生顯著下滑.Zhang 等[91]基于串列雙圓柱的4 種尾流干涉區(qū),即剪切層再附著區(qū)(間距比1.5~2.5)、共同脫渦區(qū)(間距比2.5~6.5)、尾流干擾區(qū)(間距比6.5~13.0)和不干涉區(qū)(間距比13.0~20.0),提出了線性和非線性的分段變阻尼模型,能夠增強(qiáng)下游圓柱振動響應(yīng).在剪切層再附著區(qū)和不干涉區(qū),非線性變阻尼模型能夠提升下游圓柱能量俘獲功率和能量轉(zhuǎn)換效率;而在共同脫渦區(qū)和尾流干擾區(qū),線性變阻尼模型更具優(yōu)勢(如圖10 所示).

    圖10 分段阻尼下串列圓柱的能量俘獲功率[91]Fig.10 Harnessed power of two tandem cylinders with segmental damping[91]

    兩個(gè)圓柱外徑不一致的情形也得到一定關(guān)注.Chen 等[92]實(shí)驗(yàn)研究了不等徑串列圓柱流激振動響應(yīng)特性,上游圓柱外徑與下游圓柱外徑的比值為0.4.阻尼比較小時(shí),下游圓柱隨流速增大依次發(fā)生渦激振動和馳振,且存在兩者之間的過渡區(qū);阻尼比較大時(shí),低流速下的渦激振動被抑制,渦激振動幅值和范圍隨阻尼增大而減小,高流速下發(fā)生馳振,能量俘獲功率較高.下游圓柱馳振發(fā)生的臨界速度隨阻尼比增大而升高,但隨間距比增大而降低.Bai 等[93]研究了雙圓柱外徑比和間距比對流激振動能量俘獲的影響.在較小的間距比下,等徑雙圓柱的能量俘獲能力較低,隨著間距比增大,等徑雙圓柱的能量俘獲功率增大,但整體低于不等徑雙圓柱的結(jié)果.間距比3~6、直徑比0.6 的雙圓柱能量俘獲功率更高且更穩(wěn)定.Tang 等[94]研究了上游圓柱外徑高于下游圓柱外徑時(shí)(直徑比大于1.0)串列雙圓柱的流激振動能量俘獲能力.相較于等徑雙圓柱,直徑比為2.0 時(shí),能量俘獲功率的極值略有增大,但極值對應(yīng)的流速也同步升高;直徑比為3.0 時(shí),能量俘獲功率在高流速區(qū)域出現(xiàn)顯著增大的極值.

    隨著圓柱數(shù)量的增多,圓柱間的尾流干涉更為復(fù)雜.因此,一些學(xué)者關(guān)注了多根圓柱系統(tǒng)的流激振動能量轉(zhuǎn)換特性.Ding 等[95]研究了串列排布帶有PTC 三圓柱振子VIVACE 裝置的響應(yīng)特性和流體能量俘獲效果.圓柱之間的間距比固定為2.5.隨著流速增大三圓柱的振動響應(yīng)可劃分為渦激振動初始分支區(qū)、渦激振動上分支區(qū)和渦激振動向馳振轉(zhuǎn)換區(qū)和馳振區(qū).下游圓柱受到上游圓柱尾流的強(qiáng)烈干擾,最下游圓柱的渦激振動初始分支響應(yīng)被抑制,中間圓柱的尾流模式在渦激振動向馳振轉(zhuǎn)換區(qū)與上游圓柱發(fā)生變化(如圖11 所示).三圓柱系統(tǒng)的能量俘獲功率高于3 個(gè)孤立單圓柱的能量俘獲功率,能量轉(zhuǎn)換效率在渦激振動上分支較單圓柱高35%,在馳振區(qū)與單圓柱相當(dāng).Han 等[96]通過數(shù)值仿真研究了正三角形排布的三圓柱振子流體能量俘獲效果.3 個(gè)圓柱的振動完全同步,其中兩個(gè)圓柱連線形成的正三角形底邊與來流方向垂直,圓柱間距比的變化范圍為1.01~2.0.隨著間距比的減小,三圓柱系統(tǒng)易產(chǎn)生低頻的大幅值振動,間距比為1.5 時(shí)三圓柱系統(tǒng)產(chǎn)生一個(gè)高于渦激振動上分支的新響應(yīng)分支.間距比低于1.5 時(shí)三圓柱系統(tǒng)在高流速范圍內(nèi)存在更好的能量俘獲能力,間距比為2.0 時(shí)三圓柱系統(tǒng)在低流速范圍的能量轉(zhuǎn)換效率提升,綜合而言1.2 為最佳間距比.Wang 等[97]重點(diǎn)關(guān)注了三角形排布的三圓柱系統(tǒng)中上游雙圓柱尾流對下游圓柱流激振動能量俘獲性能的影響.在上游雙圓柱的尾流影響下,下游圓柱發(fā)生馳振,能量俘獲功率和轉(zhuǎn)換效率均高于孤立單圓柱的結(jié)果.隨著上游雙圓柱與下游圓柱的距離增大,尾流干擾作用減弱,下游圓柱的能量俘獲功率下降.Rabiee 等[98]研究了來流攻角對三角形排布三圓柱系統(tǒng)流體能量俘獲效果的影響.來流攻角為0°和90°時(shí),三圓柱系統(tǒng)的能量俘獲功率較三個(gè)孤立單圓柱能量俘獲功率總和提升173%;攻角為45°時(shí),三圓柱系統(tǒng)的能量俘獲功率較3 個(gè)孤立單圓柱能量俘獲功率總和提升288%.

    圖11 串列三圓柱的尾渦模式[95]Fig.11 Wake modes of three tandem cylinders[95]

    目前關(guān)于大規(guī)模圓柱群流激振動能量俘獲的相關(guān)研究相對偏少.羅竹梅等[99-100]發(fā)現(xiàn)耦合連接的四圓柱系統(tǒng)和五圓柱系統(tǒng)更有利于俘獲在低流速下的流體能量.Zhang 等[101]通過數(shù)值仿真揭示了交錯(cuò)排布四圓柱系統(tǒng)的流激振動能量俘獲性能.交錯(cuò)排布四圓柱之間的干涉作用可根據(jù)間距比劃分為3 個(gè)區(qū)域.間距比小于2.0 時(shí)為四圓柱同步振動區(qū)域,4 個(gè)圓柱的振動相位和頻率保持一致;間距比2.0~4.0 時(shí)為振動抑制區(qū)域,上游3 個(gè)圓柱的振動因旋渦脫落被下游圓柱干擾而被顯著抑制,但下游圓柱的振動因上游不穩(wěn)定流場作用被顯著增強(qiáng),尾流中的低壓渦流在上游圓柱的推動下,循環(huán)作用于下游圓柱,增強(qiáng)了下游圓柱受到的流體載荷,進(jìn)一步提高能量轉(zhuǎn)換功率;間距比大于4.0 時(shí)為振動恢復(fù)區(qū)域,四圓柱之間的相互干涉作用減弱,振動恢復(fù)至孤立單圓柱的振動水平.整體而言,交錯(cuò)排布四圓柱的能量俘獲功率隨間距增大呈下降趨勢.Kim 等[102]研究了由串列排布的4 臺單圓柱VIVACE 模組組成的柱群結(jié)構(gòu)流激振動能量轉(zhuǎn)換性能.相比于孤立的單圓柱VIVACE 裝置,多臺VIVACE 裝置組成的柱群結(jié)構(gòu)能夠利用更多的能量.四圓柱振子的VIVACE裝置能夠有效地利用流速低至0.8~1.5 m/s 的河流或潮流中蘊(yùn)含的流動能量,且對最高流速沒有限制.為了推動VIVACE 裝置的工程應(yīng)用,美國VHE(Vortex Hydro Energy)公司獲取了VIVACE 裝置的專利轉(zhuǎn)化許可,并開展了多次現(xiàn)場應(yīng)用測試.2010 年,VHE 公司在美國密歇根州休倫港的圣克萊爾河對VIVACE 裝置進(jìn)行了首次公開測試.2012 年,VHE 公司又在圣克萊爾河開展了VIVACE 裝置的第2 次測試.2013 年,VHE 公司在美國新澤西州溢油反應(yīng)研究中心開展了VIVACE 裝置的拖曳測試.2016 年,VHE 公司對最新型的VIVACE 裝置在圣克萊爾河開展了第3 次測試,重點(diǎn)關(guān)注了VIVACE裝置在海洋中的生存能力.這些現(xiàn)場測試工作為商業(yè)化應(yīng)用提供了經(jīng)驗(yàn)參考,但目前為止,VIVACE 裝置仍未開始商業(yè)化應(yīng)用.

    3.2 非圓截面柱群流激振動能量俘獲

    由于三角形、四邊形等非圓截面單柱在某些工況下具有較高的能量俘獲能力,非圓截面柱群流激振動的能量俘獲效果也引起了一定關(guān)注.

    兩個(gè)正三角形截面柱體組成的能量收集系統(tǒng)設(shè)計(jì)應(yīng)重點(diǎn)考慮排布方式和間距比.并列排布的正三角形截面柱體不利于能量利用,主要原因是兩個(gè)柱體振子之間存在強(qiáng)烈的相互干擾,振動極不穩(wěn)定,振幅均被明顯抑制.在小間距比為1.0 工況下,兩個(gè)柱體可視為同一阻流體,尾流模式為S+P 型,尾渦的尾跡較長;間距比為2.0 工況下,兩個(gè)柱體分別脫渦,由于尾渦的相互擠壓,一個(gè)柱體尾渦為2P+S 型,另一個(gè)柱體尾渦為3P 型;在間距比為2.0 工況下,兩個(gè)柱體的尾渦均為3P 型,柱體的尾渦同步脫落現(xiàn)象更為明顯,尾渦更為清晰,相互擠壓現(xiàn)象減弱[103].串列排布的正三角形截面柱體響應(yīng)耦合機(jī)制受流速、振子間距等因素的共同作用,最佳的能量俘獲間距比約為5.0.低流速時(shí)因上游柱體的遮蔽效應(yīng),下游柱體的振動被抑制,隨著流速增大下游柱體振幅逐漸增大.串列柱體能量俘獲的啟動流速為0.4 m/s,能量俘獲功率隨速度增大持續(xù)增加.整體而言,下游柱體對上游柱體的能量俘獲具有促進(jìn)作用,下游柱體的振動和能量轉(zhuǎn)換在一定程度上受到上游柱體的不利影響.交錯(cuò)排布的正三角形截面柱體在間距比4.0~6.0 時(shí),上游柱體的旋渦脫落對下游柱體的振動具有促進(jìn)作用,下游柱體具有較高的能量俘獲功率,兩個(gè)柱體的振動聯(lián)動性強(qiáng),相位差維持在穩(wěn)定范圍,有利于能量俘獲,能量俘獲功率高于串列雙柱體[104],如圖12 所示.

    圖12 交錯(cuò)三角形截面柱體中下游柱體的能量俘獲功率[104]Fig.12 Harnessed power of downstream cylinder of two staggered cylinders with triangular section[104]

    方形截面柱體在上游固定方形截面柱體尾流的作用下能量俘獲功率呈上升趨勢[105].攻角為0°時(shí),下游柱體不再出現(xiàn)馳振響應(yīng),而是在較高的流速條件下持續(xù)保持較高振幅,最高能量俘獲功率較單個(gè)方形截面柱體提升3.8 倍,但隨間距比的增大而降低.攻角為45°時(shí),由于高流速條件下尾流模式的變化,下游柱體的振幅被上游尾流顯著增強(qiáng),能量俘獲功率較單個(gè)方形截面柱體提升約3~5 倍.隨后,Tamimi 等[106]進(jìn)一步研究了可自由振動的兩個(gè)方形截面柱體的能量俘獲性能.在上游振動柱體的尾流作用下,雖然下游柱體的能量俘獲功率略有下降,但上下游兩個(gè)柱體的整體能量俘獲功率和轉(zhuǎn)換效率顯著提升.

    對于串列排布的T 形截面柱體系統(tǒng),應(yīng)當(dāng)謹(jǐn)慎設(shè)計(jì)柱體的間距[107].在大多數(shù)間距條件下,由于兩個(gè)柱體的相互干擾,上、下游柱體的振幅和能量俘獲功率均低于單個(gè)T 形截面柱體的結(jié)果.但在某些特定工況下,兩個(gè)柱體之間的相互干擾有利于流體能量俘獲.間距比3.0~6.0 時(shí),上游柱體在下游柱體的影響下發(fā)生軟馳振響應(yīng),而不是單個(gè)T 形截面柱體發(fā)生的渦激振動響應(yīng),上游圓柱的能量俘獲功率隨流速增大而持續(xù)增長,同時(shí)由于上游柱體尾流的促進(jìn)作用,下游柱體的振幅隨流速增大而增長,最終高于單個(gè)T 形截面柱體的振幅,下游柱體的能量俘獲功率也同步提升.雙T 形截面柱體系統(tǒng)能量轉(zhuǎn)換的最佳功率出現(xiàn)在間距比為6.0 的工況下.

    為了最大程度地實(shí)現(xiàn)能量俘獲,一些學(xué)者提出了不同截面形式柱體組合形成的流激振動能量俘獲系統(tǒng).Liu 等[108]提出了一種由方形截面柱體和圓形截面柱體組合構(gòu)成的流激振動能量俘獲系統(tǒng)(如圖13所示).方形截面柱體和圓形截面柱體通過彈簧-阻尼系統(tǒng)耦合連接,其中圓形截面柱體在靠近自由液面的下方布置,用于俘獲水流中的能量,方形截面柱體在靠近自由液面的上方布置,用于俘獲空氣中的能量.不同截面形式柱體的組合應(yīng)用實(shí)現(xiàn)了能量俘獲功率、能量轉(zhuǎn)換效率和流速適應(yīng)性的改善.Tamimi等[109]提出將方形截面柱體和圓形截面柱體串列排布,發(fā)現(xiàn)在上游方形截面柱體尾流作用下,下游圓形截面柱體的能量俘獲功率提高了2 倍.

    圖13 組合形式的流激振動能量俘獲系統(tǒng)示意圖[108]Fig.13 Schematic diagram of a combined FIV energy harvesting system[108]

    4 結(jié)語與展望

    流激振動是工程中十分常見的流固耦合現(xiàn)象,利用流激振動進(jìn)行潮流能開發(fā)極具工程應(yīng)用前景.本文回顧了近年來關(guān)于柱體流激振動能量俘獲的研究工作,總結(jié)了多種截面形式的單個(gè)柱體、柱群結(jié)構(gòu)流激振動能量俘獲理論與技術(shù)方面的研究進(jìn)展.

    (1)單個(gè)光滑圓柱在某些雷諾數(shù)范圍內(nèi)無法激發(fā)高幅值振動,難以實(shí)現(xiàn)能量俘獲,通過安裝被動湍流控制器PTC 能夠顯著提升能量俘獲性能.以美國密歇根大學(xué)Bernitsas 教授團(tuán)隊(duì)為代表的學(xué)者已基本揭示了PTC 參數(shù)、系統(tǒng)阻尼、系統(tǒng)剛度、雷諾數(shù)、自由液面效應(yīng)和底邊壁效應(yīng)等因素對能量俘獲性能的影響規(guī)律,初步完成了理論和技術(shù)積累.

    (2)非圓截面柱體易激發(fā)大幅值的馳振,其流激振動的響應(yīng)特性和能量俘獲被廣泛關(guān)注.某些特定來流攻角、系統(tǒng)質(zhì)量比、系統(tǒng)阻尼、系統(tǒng)剛度、雷諾數(shù)條件下三角形、四邊形、多邊形和異形等多種截面形式柱體的流激振動作用機(jī)理和能量俘獲能力已初步明確.當(dāng)柱體截面迎流面具有對稱尖銳凸起、被流面無旋渦再附著的特點(diǎn)時(shí),有利于在大阻尼下發(fā)生自激馳振,能量俘獲能力更強(qiáng).

    (3)對于大規(guī)模的流激振動能量開發(fā),各柱體振子之間存在流場干涉,需要合理設(shè)計(jì)柱體排布形式、柱體間距和系統(tǒng)阻尼等參數(shù),實(shí)現(xiàn)流體能量俘獲最大化.對于圓形截面柱群,下游圓柱在上游圓柱的遮蔽效應(yīng)下能量俘獲性能下降,需要通過系統(tǒng)自適應(yīng)阻尼控制,提高振動幅值和穩(wěn)定性,從而提升下游圓柱的能量俘獲能力.合理設(shè)計(jì)的圓柱群系統(tǒng)相較于同等數(shù)量的單個(gè)圓柱具有更高的能量俘獲功率和轉(zhuǎn)換效率.

    (4)三角形、方形和T 形截面雙柱體的流激振動能量俘獲的相關(guān)研究較少.通常上游柱體的流激振動幅值和能量俘獲性能會受到下游柱體的促進(jìn)作用.下游柱體的能量俘獲性能受到間距比的顯著影響,小間距比時(shí)下游柱體的能量俘獲性能因上游柱體的遮蔽效應(yīng)下降,隨著間距比增大,下游柱體的能量俘獲性能逐漸提升.

    本文在回顧與總結(jié)中發(fā)現(xiàn):柱體流激振動能量俘獲研究及工程應(yīng)用推廣中存在一些問題有待進(jìn)一步解決:

    (1)以雙柱體、多柱體為對象的流激振動能量俘獲性能研究成果如何拓展至大規(guī)模柱群系統(tǒng)流激振動能量開發(fā)利用有待于進(jìn)一步論證;

    (2)影響流激振動能量俘獲性能的因素眾多,包括來流速度、雷諾數(shù)、排布形式、柱體截面形式、系統(tǒng)阻尼和剛度參數(shù)等,如何綜合考慮多因素影響實(shí)現(xiàn)能量俘獲性能最優(yōu)設(shè)計(jì)是有待解決的難題;

    (3)目前關(guān)于柱體結(jié)構(gòu)流激振動能量俘獲功率和轉(zhuǎn)換效率的結(jié)果大多基于理論、數(shù)值分析和室內(nèi)模型實(shí)驗(yàn),仍有待通過樣機(jī)測試和示范性工程應(yīng)用加以檢驗(yàn).

    猜你喜歡
    振動
    振動的思考
    某調(diào)相機(jī)振動異常診斷分析與處理
    振動與頻率
    This “Singing Highway”plays music
    具非線性中立項(xiàng)的廣義Emden-Fowler微分方程的振動性
    中立型Emden-Fowler微分方程的振動性
    基于ANSYS的高速艇艉軸架軸系振動響應(yīng)分析
    船海工程(2015年4期)2016-01-05 15:53:26
    主回路泵致聲振動分析
    UF6振動激發(fā)態(tài)分子的振動-振動馳豫
    帶有強(qiáng)迫項(xiàng)的高階差分方程解的振動性
    亚洲欧美日韩卡通动漫| 岛国在线免费视频观看| av女优亚洲男人天堂| 欧美国产日韩亚洲一区| 无限看片的www在线观看| 久久久久久九九精品二区国产| 午夜a级毛片| 日韩欧美免费精品| 三级毛片av免费| 精品久久久久久久久久久久久| 日本熟妇午夜| av国产免费在线观看| 国产精品爽爽va在线观看网站| 国产亚洲av嫩草精品影院| 欧美日韩亚洲国产一区二区在线观看| 老司机深夜福利视频在线观看| 嫩草影视91久久| 久久精品国产清高在天天线| 国产成人影院久久av| 国产精品久久电影中文字幕| 国产乱人伦免费视频| 2021天堂中文幕一二区在线观| 亚洲精品国产精品久久久不卡| 亚洲最大成人手机在线| 国产成人av激情在线播放| 国产老妇女一区| 最新中文字幕久久久久| 丝袜美腿在线中文| 亚洲国产中文字幕在线视频| 亚洲av二区三区四区| 国产真人三级小视频在线观看| 听说在线观看完整版免费高清| 亚洲人成电影免费在线| 国产激情偷乱视频一区二区| 久久精品人妻少妇| 嫩草影院精品99| 日韩有码中文字幕| 中文字幕av在线有码专区| 悠悠久久av| 国产成年人精品一区二区| 一区二区三区免费毛片| 在线播放国产精品三级| 亚洲成人免费电影在线观看| 亚洲av免费在线观看| 免费看十八禁软件| 国产高清videossex| 老熟妇乱子伦视频在线观看| 国产免费av片在线观看野外av| 亚洲欧美日韩卡通动漫| 99热这里只有精品一区| 欧美不卡视频在线免费观看| 淫妇啪啪啪对白视频| 色老头精品视频在线观看| 国产野战对白在线观看| 偷拍熟女少妇极品色| 一进一出抽搐动态| 精品久久久久久久人妻蜜臀av| 两个人的视频大全免费| 欧美性猛交黑人性爽| 熟女人妻精品中文字幕| 免费一级毛片在线播放高清视频| 神马国产精品三级电影在线观看| 国产三级黄色录像| 日本撒尿小便嘘嘘汇集6| 日韩欧美在线乱码| 在线视频色国产色| 99久久无色码亚洲精品果冻| 一卡2卡三卡四卡精品乱码亚洲| 一个人看视频在线观看www免费 | 亚洲欧美日韩高清专用| 国产激情偷乱视频一区二区| 日韩欧美精品v在线| 深爱激情五月婷婷| 最新在线观看一区二区三区| 日本熟妇午夜| 久久久久久久久大av| 欧美成人一区二区免费高清观看| 观看免费一级毛片| 欧美一区二区精品小视频在线| 国产精品日韩av在线免费观看| bbb黄色大片| 一区二区三区国产精品乱码| 国产av在哪里看| 嫩草影院精品99| 久久久国产成人精品二区| 特级一级黄色大片| 美女大奶头视频| 51午夜福利影视在线观看| 中出人妻视频一区二区| 久久国产精品影院| 长腿黑丝高跟| 日本撒尿小便嘘嘘汇集6| 国产伦一二天堂av在线观看| 欧美av亚洲av综合av国产av| 亚洲熟妇熟女久久| 国产伦一二天堂av在线观看| 欧美大码av| 国产精品99久久久久久久久| 中文字幕av成人在线电影| 欧美+亚洲+日韩+国产| 神马国产精品三级电影在线观看| 日本免费一区二区三区高清不卡| 亚洲精品亚洲一区二区| 国产男靠女视频免费网站| 色综合亚洲欧美另类图片| 极品教师在线免费播放| 国产精华一区二区三区| 精品日产1卡2卡| 色尼玛亚洲综合影院| 国产久久久一区二区三区| 国产成人a区在线观看| 午夜两性在线视频| 婷婷精品国产亚洲av在线| 日韩欧美一区二区三区在线观看| 国产精品免费一区二区三区在线| 亚洲精品色激情综合| a级毛片a级免费在线| 在线观看午夜福利视频| 女人被狂操c到高潮| 怎么达到女性高潮| 99久久九九国产精品国产免费| 国产精品久久久久久亚洲av鲁大| 亚洲,欧美精品.| 亚洲七黄色美女视频| 可以在线观看的亚洲视频| 日本在线视频免费播放| 亚洲一区二区三区不卡视频| 国产成人系列免费观看| 青草久久国产| 国产精品99久久久久久久久| 精品乱码久久久久久99久播| 51国产日韩欧美| 韩国av一区二区三区四区| 美女免费视频网站| 国产亚洲av嫩草精品影院| 免费一级毛片在线播放高清视频| 欧美在线黄色| 午夜精品在线福利| 欧美乱色亚洲激情| 丰满人妻一区二区三区视频av | 看免费av毛片| 久久国产乱子伦精品免费另类| 岛国在线免费视频观看| 欧美一级a爱片免费观看看| 午夜福利高清视频| 欧美成人性av电影在线观看| 亚洲人成网站高清观看| 男女午夜视频在线观看| 一个人观看的视频www高清免费观看| 中文字幕人妻熟人妻熟丝袜美 | 久久精品国产99精品国产亚洲性色| 日韩精品中文字幕看吧| 天堂√8在线中文| 欧美色欧美亚洲另类二区| 欧美精品啪啪一区二区三区| www.999成人在线观看| 51国产日韩欧美| www日本在线高清视频| 久久久成人免费电影| 久久午夜亚洲精品久久| 欧美乱色亚洲激情| 亚洲国产欧洲综合997久久,| 少妇的丰满在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲熟妇中文字幕五十中出| 色av中文字幕| 一进一出抽搐gif免费好疼| 欧美性猛交╳xxx乱大交人| 成人性生交大片免费视频hd| 99热这里只有是精品50| 欧美又色又爽又黄视频| 久久久久久久午夜电影| 亚洲av免费高清在线观看| 亚洲人成网站在线播| 在线播放国产精品三级| 精品一区二区三区av网在线观看| 国产精品亚洲美女久久久| 国产黄片美女视频| 日韩成人在线观看一区二区三区| 久久精品国产亚洲av涩爱 | 99久久综合精品五月天人人| 国产av麻豆久久久久久久| 国产午夜福利久久久久久| 18禁国产床啪视频网站| 天美传媒精品一区二区| www国产在线视频色| 熟女电影av网| 精品久久久久久久毛片微露脸| 午夜激情福利司机影院| 久久精品91蜜桃| 18禁在线播放成人免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 又粗又爽又猛毛片免费看| 一夜夜www| 亚洲天堂国产精品一区在线| 大型黄色视频在线免费观看| 国产高清有码在线观看视频| 国产伦精品一区二区三区四那| 3wmmmm亚洲av在线观看| 国产成人aa在线观看| 国产精品乱码一区二三区的特点| 亚洲内射少妇av| 精品日产1卡2卡| 每晚都被弄得嗷嗷叫到高潮| 嫩草影院入口| 亚洲午夜理论影院| 久久精品国产清高在天天线| av片东京热男人的天堂| 国产一区二区三区视频了| 天美传媒精品一区二区| 亚洲在线观看片| 久久精品人妻少妇| 欧美成狂野欧美在线观看| 搡老熟女国产l中国老女人| 婷婷六月久久综合丁香| aaaaa片日本免费| 亚洲色图av天堂| 国产成人福利小说| 九色国产91popny在线| 日韩成人在线观看一区二区三区| 18禁黄网站禁片午夜丰满| 日韩精品中文字幕看吧| 18禁在线播放成人免费| 天堂av国产一区二区熟女人妻| 亚洲美女视频黄频| 中出人妻视频一区二区| 每晚都被弄得嗷嗷叫到高潮| 日本 欧美在线| 一区二区三区国产精品乱码| 久久精品国产亚洲av香蕉五月| 亚洲成人久久爱视频| 激情在线观看视频在线高清| 欧美性猛交╳xxx乱大交人| 日本一本二区三区精品| 精品国产三级普通话版| 最新在线观看一区二区三区| 亚洲 国产 在线| 日韩欧美精品免费久久 | 国产av在哪里看| 美女cb高潮喷水在线观看| 欧美绝顶高潮抽搐喷水| 亚洲熟妇熟女久久| xxx96com| 国产主播在线观看一区二区| 女人被狂操c到高潮| 我的老师免费观看完整版| 免费搜索国产男女视频| 可以在线观看毛片的网站| 欧美日韩国产亚洲二区| 久久久久国内视频| 午夜福利在线观看吧| 欧美成狂野欧美在线观看| 欧美极品一区二区三区四区| 免费人成在线观看视频色| 一级a爱片免费观看的视频| 男女视频在线观看网站免费| 淫妇啪啪啪对白视频| 国产一级毛片七仙女欲春2| 亚洲精品粉嫩美女一区| 亚洲一区高清亚洲精品| 三级男女做爰猛烈吃奶摸视频| 女生性感内裤真人,穿戴方法视频| 亚洲欧美激情综合另类| 国内精品一区二区在线观看| 久久99热这里只有精品18| 一个人观看的视频www高清免费观看| 成熟少妇高潮喷水视频| 欧美日韩黄片免| 琪琪午夜伦伦电影理论片6080| 免费在线观看成人毛片| 国产乱人伦免费视频| 91麻豆精品激情在线观看国产| 18禁裸乳无遮挡免费网站照片| 精品99又大又爽又粗少妇毛片 | 欧美成人性av电影在线观看| 精品久久久久久成人av| 亚洲欧美日韩东京热| 中文字幕人成人乱码亚洲影| 午夜福利免费观看在线| 日本 欧美在线| 丰满人妻熟妇乱又伦精品不卡| 在线观看66精品国产| 床上黄色一级片| 少妇人妻精品综合一区二区 | 国产野战对白在线观看| 国产一区二区在线av高清观看| 国产精华一区二区三区| 国产色爽女视频免费观看| 老汉色∧v一级毛片| 久久久久久久午夜电影| tocl精华| 午夜精品一区二区三区免费看| 国产高潮美女av| 高清毛片免费观看视频网站| 黄色日韩在线| 日韩av在线大香蕉| 国产精品一区二区三区四区免费观看 | 日日干狠狠操夜夜爽| 欧美另类亚洲清纯唯美| 亚洲最大成人中文| 又粗又爽又猛毛片免费看| 99精品欧美一区二区三区四区| 久9热在线精品视频| 日日摸夜夜添夜夜添小说| 男女午夜视频在线观看| 免费看十八禁软件| 国产精品 欧美亚洲| 精品人妻偷拍中文字幕| 性色av乱码一区二区三区2| 精品久久久久久成人av| 日本黄大片高清| 免费电影在线观看免费观看| 欧美一区二区国产精品久久精品| 女警被强在线播放| 日韩欧美在线乱码| 村上凉子中文字幕在线| 99久久无色码亚洲精品果冻| 国产色爽女视频免费观看| 日韩人妻高清精品专区| 免费观看人在逋| 国产成人a区在线观看| 日韩欧美 国产精品| 国产高清有码在线观看视频| 免费人成在线观看视频色| 美女大奶头视频| 午夜老司机福利剧场| 国产又黄又爽又无遮挡在线| 国产精品99久久99久久久不卡| 国产亚洲欧美在线一区二区| 国产单亲对白刺激| 欧美绝顶高潮抽搐喷水| 搡老岳熟女国产| 精品免费久久久久久久清纯| 又爽又黄无遮挡网站| 变态另类成人亚洲欧美熟女| 最近视频中文字幕2019在线8| 久久午夜亚洲精品久久| 脱女人内裤的视频| 网址你懂的国产日韩在线| av天堂在线播放| 极品教师在线免费播放| 国产淫片久久久久久久久 | 日韩欧美三级三区| 亚洲七黄色美女视频| 国产麻豆成人av免费视频| 亚洲在线自拍视频| 成人av一区二区三区在线看| 偷拍熟女少妇极品色| 欧美色欧美亚洲另类二区| 搞女人的毛片| 日韩 欧美 亚洲 中文字幕| 嫩草影院精品99| 非洲黑人性xxxx精品又粗又长| 在线十欧美十亚洲十日本专区| 一a级毛片在线观看| 操出白浆在线播放| 国语自产精品视频在线第100页| 宅男免费午夜| 在线免费观看不下载黄p国产 | 国产精品国产高清国产av| 色哟哟哟哟哟哟| 女生性感内裤真人,穿戴方法视频| 国产精品精品国产色婷婷| 色av中文字幕| 欧美三级亚洲精品| 欧美性猛交╳xxx乱大交人| bbb黄色大片| 国产伦一二天堂av在线观看| 叶爱在线成人免费视频播放| 制服人妻中文乱码| 丁香欧美五月| 高清毛片免费观看视频网站| 中文字幕人妻熟人妻熟丝袜美 | 美女黄网站色视频| 69av精品久久久久久| 精品午夜福利视频在线观看一区| 一个人看视频在线观看www免费 | 狂野欧美白嫩少妇大欣赏| 欧美精品啪啪一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 欧美激情久久久久久爽电影| 国产精品久久视频播放| 动漫黄色视频在线观看| 国产真实乱freesex| 国产爱豆传媒在线观看| 久久久久性生活片| 亚洲自拍偷在线| 色视频www国产| 亚洲成人精品中文字幕电影| 99久久综合精品五月天人人| 婷婷六月久久综合丁香| 真实男女啪啪啪动态图| 一二三四社区在线视频社区8| 脱女人内裤的视频| or卡值多少钱| 成人永久免费在线观看视频| 91在线精品国自产拍蜜月 | 毛片女人毛片| 亚洲专区中文字幕在线| 美女高潮的动态| 色在线成人网| 精品日产1卡2卡| 国产主播在线观看一区二区| 真实男女啪啪啪动态图| 国产一区二区激情短视频| 午夜视频国产福利| 欧美日韩乱码在线| 国产精品98久久久久久宅男小说| 午夜福利视频1000在线观看| 国产精品综合久久久久久久免费| 国产成人av教育| av中文乱码字幕在线| 性欧美人与动物交配| av天堂在线播放| 成人性生交大片免费视频hd| 中文亚洲av片在线观看爽| 少妇人妻一区二区三区视频| 中文资源天堂在线| 亚洲精品成人久久久久久| 午夜激情欧美在线| 在线天堂最新版资源| 欧美日韩乱码在线| 最新美女视频免费是黄的| 手机成人av网站| 丰满的人妻完整版| 国内精品久久久久久久电影| 三级男女做爰猛烈吃奶摸视频| 露出奶头的视频| 亚洲国产精品久久男人天堂| 中文字幕av成人在线电影| 最近最新免费中文字幕在线| 真实男女啪啪啪动态图| 亚洲乱码一区二区免费版| 国产亚洲av嫩草精品影院| 国产毛片a区久久久久| 黄色视频,在线免费观看| 久久久久久久亚洲中文字幕 | 1000部很黄的大片| 最后的刺客免费高清国语| 天天一区二区日本电影三级| 国产在视频线在精品| 国内精品久久久久精免费| 欧美日韩中文字幕国产精品一区二区三区| 国产精品亚洲美女久久久| 亚洲欧美日韩卡通动漫| 日本黄色视频三级网站网址| 午夜激情欧美在线| eeuss影院久久| 成人18禁在线播放| 国产精品一区二区三区四区久久| 9191精品国产免费久久| 成人午夜高清在线视频| 九九在线视频观看精品| 成人av在线播放网站| 久久久久亚洲av毛片大全| 又爽又黄无遮挡网站| 久久婷婷人人爽人人干人人爱| 亚洲av成人精品一区久久| 亚洲内射少妇av| 久久久久免费精品人妻一区二区| 天天躁日日操中文字幕| 亚洲成人久久性| 亚洲成a人片在线一区二区| 美女免费视频网站| 男女那种视频在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲色图av天堂| 成年人黄色毛片网站| 波多野结衣高清无吗| www国产在线视频色| 久久精品夜夜夜夜夜久久蜜豆| 国产高清三级在线| 中文字幕熟女人妻在线| 一个人看的www免费观看视频| 国产精华一区二区三区| 国产av在哪里看| 亚洲国产中文字幕在线视频| 国产av不卡久久| 亚洲一区二区三区不卡视频| 久久久久国内视频| 97碰自拍视频| 成人亚洲精品av一区二区| 99热只有精品国产| 搡老熟女国产l中国老女人| 欧美成人a在线观看| 我的老师免费观看完整版| 美女高潮的动态| 国产伦精品一区二区三区视频9 | 精品国产三级普通话版| 中文字幕久久专区| 香蕉丝袜av| 欧美日韩综合久久久久久 | 国产一区二区在线av高清观看| 国产成+人综合+亚洲专区| 嫩草影院入口| 中文字幕av成人在线电影| 老司机福利观看| 日本精品一区二区三区蜜桃| 亚洲一区二区三区不卡视频| 高清毛片免费观看视频网站| 国产av在哪里看| 欧美国产日韩亚洲一区| 丰满乱子伦码专区| h日本视频在线播放| 久久午夜亚洲精品久久| 成人三级黄色视频| 亚洲在线自拍视频| 天美传媒精品一区二区| 国产极品精品免费视频能看的| 很黄的视频免费| 久久精品国产亚洲av香蕉五月| 日日摸夜夜添夜夜添小说| 两个人视频免费观看高清| 成人特级黄色片久久久久久久| 国内精品久久久久久久电影| 色综合亚洲欧美另类图片| 成人精品一区二区免费| 日本 欧美在线| 亚洲专区中文字幕在线| 日韩欧美在线二视频| 亚洲精品国产精品久久久不卡| 757午夜福利合集在线观看| 精品欧美国产一区二区三| 国产精品影院久久| 级片在线观看| 精品久久久久久久毛片微露脸| 国产精品1区2区在线观看.| 国产麻豆成人av免费视频| 日本 av在线| 麻豆一二三区av精品| 色播亚洲综合网| 我的老师免费观看完整版| 99久久综合精品五月天人人| 国产av不卡久久| 国内精品一区二区在线观看| av女优亚洲男人天堂| 欧美bdsm另类| 99精品在免费线老司机午夜| 丝袜美腿在线中文| 日韩人妻高清精品专区| 国产精品永久免费网站| 亚洲av电影在线进入| 国产精品1区2区在线观看.| 少妇人妻精品综合一区二区 | 久久香蕉国产精品| 国产乱人伦免费视频| 热99re8久久精品国产| 久久亚洲精品不卡| 日韩欧美在线乱码| 高清毛片免费观看视频网站| 日韩欧美精品v在线| 欧美一级a爱片免费观看看| 狠狠狠狠99中文字幕| 国产高清视频在线播放一区| 国产成人啪精品午夜网站| 久久久国产精品麻豆| 小说图片视频综合网站| 国产97色在线日韩免费| 国产精品一区二区三区四区久久| 757午夜福利合集在线观看| 国产真实乱freesex| 久久亚洲精品不卡| 99久久99久久久精品蜜桃| 午夜免费激情av| 成人三级黄色视频| 日日夜夜操网爽| 国产精品 国内视频| 又粗又爽又猛毛片免费看| 级片在线观看| 日本免费a在线| 亚洲中文字幕一区二区三区有码在线看| 亚洲第一欧美日韩一区二区三区| 十八禁网站免费在线| 国产精品乱码一区二三区的特点| 精品熟女少妇八av免费久了| 日韩欧美一区二区三区在线观看| 国产亚洲精品一区二区www| 一本一本综合久久| 少妇丰满av| 亚洲av中文字字幕乱码综合| 午夜老司机福利剧场| 欧美日本亚洲视频在线播放| 国内精品久久久久精免费| 免费搜索国产男女视频| av天堂在线播放| 精品久久久久久久久久免费视频| 搡老妇女老女人老熟妇| 两个人的视频大全免费| 日本在线视频免费播放| 免费大片18禁| 久久久久久人人人人人| av黄色大香蕉| 亚洲七黄色美女视频| 好男人电影高清在线观看| 黄色丝袜av网址大全| av在线天堂中文字幕| 国产午夜福利久久久久久| 色av中文字幕| 精品久久久久久久毛片微露脸| 国产午夜福利久久久久久| 99热这里只有精品一区| 极品教师在线免费播放| 国产免费一级a男人的天堂| 熟女少妇亚洲综合色aaa.| 嫩草影院精品99| 色哟哟哟哟哟哟| 欧美国产日韩亚洲一区| 制服丝袜大香蕉在线| 精品电影一区二区在线| 国产综合懂色| 18禁美女被吸乳视频| 男人舔奶头视频| xxx96com| 国产伦人伦偷精品视频| 麻豆成人av在线观看|