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

    鈍體表面附著物對低速水流壓電俘能器性能影響研究1)

    2023-08-06 08:46:08孫維鵬劉宸涵郁小彬鐘可欣趙道利
    力學(xué)學(xué)報(bào) 2023年7期
    關(guān)鍵詞:圓柱型俘能器附著物

    孫維鵬 劉宸涵 郁小彬 胡 珅 鐘可欣 趙道利,2)

    * (西安理工大學(xué)省部共建西北旱區(qū)生態(tài)水利國家重點(diǎn)實(shí)驗(yàn)室,西安 710048)

    ? (華東天荒坪抽水蓄能有限責(zé)任公司,杭州 310012)

    引言

    近年來,能源危機(jī)愈演愈烈[1],為了應(yīng)對新環(huán)境下全球能源格局發(fā)生的深刻調(diào)整[2],高效開發(fā)利用可再生能源迫在眉睫[3-6].水能作為占比最高的可再生能源,如何高效利用一直備受關(guān)注.天然河道中的高速水流能可通過水輪發(fā)電機(jī)組加以利用[7-8],這一技術(shù)目前已相對成熟;而湖泊與小型河流中的低速水流能雖可作為新能源的重要補(bǔ)充,但因其分布廣、能量低而存在利用困難的問題.馳振[9]作為一種能夠發(fā)生在低速水流中的現(xiàn)象,有流速響應(yīng)范圍寬和無鎖定限制的特點(diǎn).因此基于馳振設(shè)計(jì)且已在流場中展現(xiàn)出游優(yōu)異性能的馳振式壓電俘能器為低速水流能的利用提供了解決方案[10].利用馳振式壓電俘能器將低速水流能量轉(zhuǎn)化為電能,可以為微機(jī)電系統(tǒng)[11-12]和無線傳感器[13-14]供電.

    目前已經(jīng)存在一些關(guān)于利用低速水流能的壓電俘能器研究.曹東興等[15]設(shè)計(jì)一種利用永磁鐵增強(qiáng)流致振動的壓電俘能器.實(shí)驗(yàn)結(jié)果表明: 水流速度為0.5 m/s 時,輸出功率的均方根值為120 μW,與無附加磁體的壓電俘能器相比提高了57.8%.單小彪等[16]提出一種在水中完全淹沒的尾流激振式壓電俘能器.結(jié)果表明,當(dāng)鈍體的直徑為15 mm 時,在上游距其中心20 mm 處放置直徑5 mm 擾流柱能夠使壓電俘能器輸出23.35 V 的電壓,與沒有擾流柱相比輸出電壓提升了51%.Liu 等[17]基于馳振的觸發(fā)方式,提出一種新的時間變流型方法.該方法解釋了傳統(tǒng)的時間平均流型方法難以解釋的橫向馳振機(jī)制.Lian 等[18]研究鈍體截面為正方形、三角形和圓形的壓電俘能器在水中的性能.實(shí)驗(yàn)表明: 截面為三角形的鈍體,其振幅比大于1.8;圓柱型鈍體的能量轉(zhuǎn)化率最高為80%.燕翔等[19]實(shí)驗(yàn)研究T 型截面鈍體的馳振特性,并進(jìn)一步探究其在不同外力條件及不同截面寬高比和阻尼比條件下的振動響應(yīng)差異.He 等[20]通過研究證實(shí)改變鈍體的幾何形狀,可以實(shí)現(xiàn)耦合或獨(dú)立的馳振和渦激振動.Zhao等[21]采用圓柱形擾流柱,對鈍體截面形狀為圓形、D 形和倒D 形的壓電俘能器進(jìn)行水槽實(shí)驗(yàn)測試.結(jié)果表明,3 種壓電俘能器的采集功率分別達(dá)到570.3,596.4和1074 W/m3,是無尾流激勵情況的43.2,25.3 和31 倍.張軍等[22]對三棱柱的流致振動進(jìn)行了實(shí)驗(yàn)研究,并評估其能量轉(zhuǎn)化效率.Song 等[23]研究圓柱型鈍體壓電俘能器在不同直徑、質(zhì)量和阻尼比條件下的輸出性能.研究證實(shí)在一定范圍內(nèi),降低質(zhì)量和增加鈍體直徑都可以提高壓電俘能器的輸出.以上研究主要通過改變鈍體的形狀來優(yōu)化壓電俘能器,而在鈍體表面添加小型附著物的研究還很少.

    存在于鈍體表面的附著物能夠影響壓電俘能器系統(tǒng)的振動與輸出[24-25],目前已存在一些將附著物與壓電俘能器相結(jié)合的應(yīng)用: Zhu 等[26]提出一種由圓柱體和五角星組成的新型壓電俘能器,雙向流固耦合仿真的結(jié)果表明其功率密度可達(dá)到885 W/m3.Hu 等[27]通過調(diào)整圓形、三角形和矩形桿件在圓柱形鈍體上的位置來提升設(shè)備的風(fēng)能量采集效率.實(shí)驗(yàn)證實(shí)當(dāng)安裝的三棱柱夾角為60°時,能量采集器的性能最佳.Zhu 等[28]設(shè)計(jì)了在圓柱體表面安裝對稱鰭狀條帶的新鈍體,數(shù)值模擬的結(jié)果表明收集到的最大功率超過60 W/m.目前關(guān)于附著物在壓電俘能器上的研究都只局限于其形狀及安裝形式,缺少在鈍體上進(jìn)行切除特定形狀的嘗試.附著物的凸起高度和凹陷深度值對壓電俘能器性能的影響有待研究.

    綜上所述,本文主要研究不同凸起高度/凹陷深度和形狀的附著物對馳振式壓電俘能器輸出特性的影響.引入水的附加質(zhì)量力改進(jìn)現(xiàn)有的分段參數(shù)模型.基于準(zhǔn)穩(wěn)態(tài)假設(shè),建立馳振的水動力模型并計(jì)算出馳振力.水槽實(shí)驗(yàn)驗(yàn)證了數(shù)學(xué)模型的準(zhǔn)確性,進(jìn)一步通過模型分析壓電俘能器的輸出功率和懸臂梁振幅與凸起高度/凹陷深度之間的關(guān)系.最后基于數(shù)值模擬從流體動力學(xué)的角度闡述超表面對系統(tǒng)影響的根本原因.

    1 數(shù)學(xué)模型

    本節(jié)利用拓展的哈密頓原理推導(dǎo)出壓電俘能器的機(jī)械控制方程,引入高斯定理建立電場-位移控制方程,并基于準(zhǔn)穩(wěn)態(tài)假設(shè)計(jì)算馳振的水動力及力矩,利用伽遼金法離散懸臂梁的位移,并在此基礎(chǔ)上解耦控制方程,得出輸出功率和懸臂梁位移的近似解.

    1.1 俘能器機(jī)械控制方程

    本文研究的壓電俘能器由鈍體、懸臂梁和MFC壓電片組成,如圖1 所示.懸臂梁的上端固定,中部黏貼有壓電片,下端與鈍體連接.x軸為懸臂梁的長度方向,y軸為水流方向,受到水流沖擊的鈍體在z方向上發(fā)生振動,考慮到振幅較大,選擇脆性低、韌性高的宏觀纖維復(fù)合材料(macro fiber composites,MFC)作為壓電材料.

    圖1 壓電俘能器Fig.1 Piezoelectric energy harvester

    橢圓柱鈍體的表面裝配了人工設(shè)計(jì)的附著物,其形狀采用橢球型和圓柱型兩種,各自的最大高度被定義為特征值,分別為Eell和Ecy,如圖2 所示.橢球型和圓柱型附著物的特征值均設(shè)置為6,4,2,-1 和-2 mm,其中正值為凸起,負(fù)值為凹陷.

    圖2 附著物的結(jié)構(gòu)圖Fig.2 Structural diagram of the attachments

    基于拓展的哈密頓原理[29-30],壓電俘能器的控制方程可表示為

    式中,ν(x,t) 表示t時刻懸臂梁x處的位置在z方向上的位移;ca和cs分別表示懸臂梁的氣動阻尼系數(shù)和黏性應(yīng)變系數(shù);δ表示克羅內(nèi)克符號,其值在x=l時取1,反之取0;?P為壓電耦合系數(shù),?P=-e31bp·(z1+z2)/2,e31為壓電應(yīng)力系數(shù),e31=Epd31,d31為壓電層的應(yīng)變系數(shù);V(t)壓電材料的電壓;Ftip和Mtip分別表示懸臂梁的馳振力和力矩.

    1.2 電場-位移控制方程

    引入高斯定理[31]

    系統(tǒng)電場-位移控制方程可表示為

    式中Cp是壓電俘能器的電容,表示為

    1.3 馳振的水動力計(jì)算

    為了建立馳振的水動力模型,引入準(zhǔn)穩(wěn)態(tài)假設(shè).而假設(shè)中的馳振力會受到攻角α的影響.攻角可表示為

    而攻角與馳振系數(shù)Cy存在如下關(guān)系

    式中a1和a3可通過實(shí)驗(yàn)獲得.

    考慮到附加質(zhì)量力對馳振力的影響,鈍體單位長度的馳振力可表示為[32]

    式中CM表示附加質(zhì)量系數(shù).

    懸臂梁所受到的馳振水動力和力矩可由下式計(jì)算

    1.4 馳振模型的機(jī)電解耦

    應(yīng)用伽遼金法[33-34]離散懸臂梁的位移ν(x,t),將其劃分為空間變量和時間變量

    式中,Qi(t)和φi(t)分別為懸臂梁的第i階模態(tài)坐標(biāo)和模態(tài).

    馳振一階模態(tài)的水動力可以根據(jù)f(t)=φ(l)Ftip+φ′(l)Mtip來計(jì)算.依據(jù)式 (7),得出一階模態(tài)的水動力

    式中,k1,k2和k3分別為

    Abdelkefi 等[35]的研究證實(shí),懸臂梁的振動主要受一階模態(tài)支配,Bibo 等[36]也發(fā)現(xiàn)一階模態(tài)滿足被捕獲的條件.因此本文采用一階模態(tài)來簡化懸臂梁的機(jī)電耦合模型

    式中,ζ表示阻尼系數(shù);ω表示壓電俘能器的一階固有頻率;θp是壓電耦合項(xiàng),表示為 θP=?P[φ′(l2)-φ′(l1)];f(t)表示一階振型的水動力.解耦的控制方程[37]由式 (9) 和式 (10) 得

    式中,?和Ce分別為壓電俘能器的修正頻率和阻尼,表示為

    進(jìn)一步可得位移模態(tài)坐標(biāo)與電壓之間的關(guān)系

    模態(tài)坐標(biāo)的振幅為

    由式 (14) 確定模態(tài)振幅Q0后,壓電俘能器的采集功率可表示為

    2 實(shí)驗(yàn)設(shè)置

    本研究的實(shí)驗(yàn)部分在圖3 所示的U 型水槽中進(jìn)行.水從地下水庫通過泵與進(jìn)水管被引入矩形薄壁堰中,薄壁堰上游的蜂窩網(wǎng)可以降低水的湍流,從而減少水流不穩(wěn)定對上游水位探針讀數(shù)的影響;水流離開水堰流經(jīng)第一段長直段,在這里被充分穩(wěn)定以滿足實(shí)驗(yàn)段對穩(wěn)定水流的需求;隨后水流流入實(shí)驗(yàn)段,實(shí)驗(yàn)段的中部安裝了壓電俘能器;最后水流由尾水口流回地下水庫.實(shí)驗(yàn)段的長度為6.2 m,其斷面寬度為0.4 m.實(shí)驗(yàn)中的流速是通過泵閥門及尾水口上方的電動閥門控制的,轉(zhuǎn)動泵閥門能夠調(diào)節(jié)進(jìn)入水槽內(nèi)的流量;調(diào)整電動閥門的高度可控制水槽內(nèi)的水位,然后利用流量計(jì)算公式得出流速.

    圖3 U 型水槽Fig.3 U-shape water channel

    壓電俘能器的鈍體、懸臂梁及MFC 壓電片的參數(shù)如表1 所示.本實(shí)驗(yàn)的流速范圍為0.26~0.51 m/s,取0.05 m/s 為一流速步長.實(shí)驗(yàn)中利用示波器測試出與壓電片串聯(lián)的負(fù)載電阻的電壓并將數(shù)據(jù)儲存在計(jì)算機(jī),然后計(jì)算出電壓的均方根值(RMS voltage).實(shí)驗(yàn)的負(fù)載電阻值為3.5×105Ω.

    表1 鈍體、懸臂梁和壓電片的物理參數(shù)Table 1 Physical parameters of bluff body,cantilever beam and piezoelectric sheet

    3 模型驗(yàn)證

    為了驗(yàn)證所建立數(shù)學(xué)模型的準(zhǔn)確性,分別繪制0.41,0.46 和0.51 m/s 流速下橢球型和圓柱型附著物在2 mm 凹陷條件下電壓與預(yù)測值的對照,如圖4所示.圖中的實(shí)線為實(shí)驗(yàn)測得的電壓曲線,是從單次140 s 的實(shí)驗(yàn)電壓中選取相對穩(wěn)定的4 s 繪制的;虛線為電壓的數(shù)學(xué)模型計(jì)算值.圖中實(shí)線和虛線基本重合,二者之間的微小差異是由水流的輕微波動引起的.所建立數(shù)學(xué)模型的計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果十分吻合,由此驗(yàn)證了數(shù)學(xué)模型的準(zhǔn)確性.

    圖4 模型計(jì)算和實(shí)驗(yàn)結(jié)果的時域曲線對比Fig.4 Comparison of time histories of the model calculations and experimental results

    4 壓電俘能器輸出特性和懸臂梁振幅的分析

    結(jié)合實(shí)驗(yàn)值和數(shù)學(xué)模型的近似解,分析橢球型和圓柱型附著物在不同的特征值和流速條件下對馳振式壓電俘能器RMS 輸出功率和鈍體振幅的影響.

    4.1 流速對壓電俘能器輸出特性影響的分析

    橢球型附著物的壓電俘能器RMS 輸出功率隨流速的變化如圖5(a) 所示,可以看出在所有特征值下,壓電俘能器的輸出均隨著流速的增加而增加.這是由于鈍體的振幅隨流速而增大,從而帶動壓電片發(fā)生更大程度的形變進(jìn)而產(chǎn)生更高的電壓,符合馳振的特性.此外附著物的特征值對壓電俘能器輸出功率的影響也不同.在U=0.51 m/s 的流速下,Eell=-2 mm 的壓電俘能器RMS 輸出功率為1.411 mW,相比于橢圓柱增幅為69.88%;而部分特征值的附著物則抑制壓電俘能器的輸出: 當(dāng)Eell=6,4 mm 時,壓電俘能器的RMS 輸出功率在各個流速下始終低于未包裹附著物的橢圓柱,其中U=0.51 m/s 時,6 mm 特征值的橢球型附著物使輸出功率降低了80.02%.

    圖5 流速對壓電俘能器RMS 輸出功率的影響Fig.5 Effect of flow velocity on the RMS output power of the piezoelectric energy harvester

    圖5(b) 展示了包裹圓柱型附著物對壓電俘能器性能的影響.與橢球型附著物相比,圓柱型附著物的施加使得壓電俘能器的出輸出整體偏低.當(dāng)Ecy=4,6 mm 時,壓電俘能器在所有流速下的輸出均被大幅度抑制,其中Ecy=6 mm 的壓電俘能器在U=0.51 m/s 的輸出為0.059 mW,降低了92.9%.此外隨著流速的不斷增加,Ecy=-1,-2 mm 在流速超過0.41 m/s 后能夠提高輸出功率,其中Ecy=-2 mm 在0.51 m/s 的流速下能夠產(chǎn)生1.308 mW 的RMS 功率.

    4.2 特征值對壓電俘能器輸出特性的影響

    4.1 中比較了不同特征值下流速對于壓電俘能器RMS 輸出功率的影響,接下來結(jié)合數(shù)學(xué)模型的計(jì)算結(jié)果,進(jìn)一步分析特征值的變化對于壓電俘能器性能的影響.圖6(a) 表明隨著特征值從-2 mm 逐漸增大到6 mm,橢球型附著物壓電俘能器的輸出先減小再增大后減小的趨勢,最大值分別在Eell=-2,2 mm 處.這說明當(dāng)特征值絕對值在0 附近,即凸起或凹陷相對較小時,附著物的施加能夠增強(qiáng)鈍體的振動,提高壓電俘能器的輸出功率,并且最大值的出現(xiàn)證明實(shí)現(xiàn)這種增強(qiáng)效果的附著物存在“最優(yōu)特征值”.在0.51 m/s 流速下,未包裹附著物的原始鈍體對應(yīng)的壓電俘能器RMS 輸出功率為0.83 mW,而分別應(yīng)用和-2 mm 特征值的橢球型與圓柱型附著物均有最大輸出,分別對應(yīng)于1.41 和1.31 mW,相比于原始俘能器分別提高了69.88%和57.83%.此外當(dāng)特征值大于2 mm 后,壓電俘能器的輸出產(chǎn)生不同程度的降低,這表明較大高度的附著物會抑制鈍體的振動.

    圖6 不同形狀附著物的壓電俘能器RMS 輸出功率隨特征值的變化Fig.6 Variation of RMS output power for energy harvesters versus eigenvalues with different shaped attachment applied

    對于圓柱型附著物如圖6(b) 所示,隨著附著物高度從最大凹陷增加到最大凸起,壓電俘能器的RMS 輸出功率在Ecy< 0 的范圍內(nèi)不斷降低,并在E=0 處達(dá)到極小值;而在Ecy> 0 時又呈現(xiàn)先增大后減小的趨勢.這與橢球型附著物的趨勢基本一致,說明兩種形狀附著物的特征值變化對壓電俘能器輸出的影響規(guī)律是基本相同的.

    此外,附著物的恰當(dāng)應(yīng)用能夠增加壓電俘能器的輸出.在橢球型附著物特征值為-2 mm、流速為0.51 m/s 的條件下,俘能器的最大RMS 輸出功率為1.411 mW,對應(yīng)的輸出功率密度ρp可被計(jì)算為11998.3 W/m3,與部分研究的俘能器輸出功率密度對比如表2 所示.

    表2 與同類研究對比Table 2 Comparison with similar studies

    4.3 特征值對懸臂梁振幅影響的分析

    為了進(jìn)一步分析附著物特征值的變化對壓電俘能器輸出的影響,通過機(jī)電解耦得到懸臂梁的底端振幅,如圖7 所示.通過與圖6 中相應(yīng)特征值條件下RMS 輸出功率的對比可以發(fā)現(xiàn),懸臂梁的振幅與壓電俘能器輸出功率的變化趨勢基本相同,二者呈正相關(guān).此外當(dāng)U=0.51 m/s 時,沒有施加附著物的懸臂梁振幅為20.25 mm,相同流速下鈍體包裹了Ecy=6 mm 附著物的懸臂梁振幅則為3.07 mm,降低了84.83%,此時懸臂梁幾乎靜止,鈍體的振動得到大幅度抑制.

    圖7 圓柱型附著物懸臂梁的振幅隨特征值的變化Fig.7 Variation of amplitude with eigenvalue of cantilever beam for cylindrical attachment

    4.4 包裹附著物的鈍體擾流數(shù)值模擬

    適當(dāng)特征值的附著物對壓電俘能器功率輸出具有增幅效應(yīng),并且無論是橢球型還是圓柱型附著物,凹陷或者凸起高度為2 mm 時,增幅效應(yīng)均最大.因此選取Eell=-2 mm,Ecy=-2 mm,Eell=2 mm,Ecy=2 mm 和E=0 的5 個鈍體進(jìn)行三維數(shù)值模擬研究.D為鈍體的迎水面寬度,計(jì)算域設(shè)置為20D×6D×2D,鈍體距上游面6D.圍繞鈍體的周圍設(shè)置加密的四面體網(wǎng)格,其他區(qū)域采用六面體網(wǎng)格.對網(wǎng)格進(jìn)行無關(guān)性驗(yàn)證,均衡了效率和精度,最終網(wǎng)格數(shù)目控制在60 萬左右,整體網(wǎng)格質(zhì)量控制在0.5 以上.為了定性地研究附著物對鈍體流致振動的影響,并考慮到計(jì)算效率及收斂性,將整個計(jì)算域的網(wǎng)格設(shè)為固定網(wǎng)格,湍流模型采用SSTk-ω,設(shè)置0.51 m/s 流速入口和環(huán)境出口,模擬時間步長為0.00125 s,計(jì)算步數(shù)為8000 步.

    圖8 繪制了5 個鈍體的壓力和渦量云圖,從左側(cè)流入的水流在受到鈍體阻礙后流速降低并在其上游形成高壓區(qū);隨后從側(cè)邊繞過加速流向下游并產(chǎn)生低壓區(qū),加速過程中水流受黏性的影響形成脫落渦.對比圖8(a) 所示的壓力云圖發(fā)現(xiàn): 包裹附著物的鈍體周圍產(chǎn)生的高、低壓區(qū)之間的壓力差明顯大于未包裹附著物情況,這也對應(yīng)更大的升力系數(shù),如表3所示.圖8(b) 所示的三維渦量云圖也表明包裹附著物后,從鈍體兩側(cè)脫落的尾渦明顯更大.絕對值為2 mm 附著物的存在,使在鈍體上形成的壓差升高,加劇流致振動,最終提高了壓電俘能器的輸出,這與試驗(yàn)及模型得出的結(jié)果是一致的.

    表3 升力系數(shù)Table 3 Lift coefficient

    圖8 U=0.51 m/s 時,Eell=2,-2 mm,Ecy=2,-2 mm 和E=0 的云圖Fig.8 Contour of Eell=2,-2 mm,Ecy=2,-2 mm and E=0 with U=0.51 m/s

    此外進(jìn)一步對比包裹兩種類型附著物鈍體的脫落尾渦,發(fā)現(xiàn)在相同特征值下,包裹橢球型附著物的鈍體尾流渦街要比包裹圓柱型附著物的情況略微紊亂,推測是因?yàn)楦街锏慕Y(jié)構(gòu)干擾了流體分層,這種干擾作用在橢球型附著物上要大于圓柱型附著物.

    5 總結(jié)

    本文采用數(shù)學(xué)模型、實(shí)驗(yàn)及數(shù)值模擬相結(jié)合的方式研究了附著物的形狀和附著物的特征值對馳振式壓電俘能器輸出的影響.基于拓展的哈密頓原理建立壓電俘能器的機(jī)械控制方程,利用高斯定理建立電場-位移控制方程,通過準(zhǔn)穩(wěn)態(tài)假設(shè)計(jì)算出馳振力,在利用伽遼金法離散懸臂梁位移的基礎(chǔ)上,實(shí)現(xiàn)對機(jī)電耦合控制模型的化簡和解耦,最終得到系統(tǒng)輸出及懸臂梁振動位移的近似解.模型計(jì)算的準(zhǔn)確性得到了水槽實(shí)驗(yàn)的驗(yàn)證.利用數(shù)值模擬研究進(jìn)一步分析了附著物在鈍體的振動過程中發(fā)揮的作用.本文主要結(jié)論如下.

    (1) 在U=0.51 m/s 和Eell=-2 mm 的條件下,壓電俘能器最大RMS 輸出功率為1.411 mW,相比于E=0 的情況增加了69.88%.

    (2) 隨著附著物的特征值從-2 mm 增加到6 mm,壓電俘能器的輸出功率均出現(xiàn)先減小再增大后減小的趨勢,凸起和凹陷分別在-2 和 2 mm 取得最大值.對于增強(qiáng)壓電俘能器的輸出而言,附著物的凸起和凹陷分別存在“最優(yōu)特征值”.

    (3) 當(dāng)U=0.51 m/s 時,Ecy=6 mm 條件下懸臂梁末端振幅為3.07 mm,與E=0 的情況相比降低了84.83%,附著物的減振效果為大跨度柔性結(jié)構(gòu)的減振提供了參考.

    (4) 數(shù)值模擬的結(jié)果表明: 特征值絕對值為2 mm 的附著物能夠令在鈍體上形成的壓差升高、升力變大,引起鈍體愈發(fā)劇烈的流致振動;并且附著物的應(yīng)用會干擾鈍體兩側(cè)的流體分層,進(jìn)而對脫落尾渦的大小產(chǎn)生影響.

    猜你喜歡
    圓柱型俘能器附著物
    汽車懸架的振動俘能與汽車平順性仿真研究
    壓電疊堆—蜂鳴片復(fù)合俘能器設(shè)計(jì)*
    板式換熱器進(jìn)口處液體分布器的數(shù)值模擬*
    一種質(zhì)量塊-彈簧自參數(shù)共振壓電俘能器研究
    壓電與聲光(2020年5期)2020-10-28 05:05:26
    一種壓電式微型俘能器
    船底涂層有什么作用?
    發(fā)明解讀刀具夾持裝置以及包括這種刀具夾持裝置的手持式電動工具
    電動工具(2018年4期)2018-08-22 02:02:28
    圓柱型功能梯度雙材料界面裂紋問題研究
    輸電線路附著物測算系統(tǒng)測算功能模塊的研究
    圓柱型正交各向異性圓板的自由振動分析
    欧美日韩乱码在线| netflix在线观看网站| av福利片在线观看| 久久人人精品亚洲av| 五月玫瑰六月丁香| 日本一二三区视频观看| 精品人妻一区二区三区麻豆 | 亚洲黑人精品在线| 免费黄网站久久成人精品 | 亚洲黑人精品在线| 又爽又黄a免费视频| 久久久久国产精品人妻aⅴ院| 国产乱人视频| 麻豆国产97在线/欧美| a在线观看视频网站| 亚洲一区高清亚洲精品| 国产熟女xx| 三级国产精品欧美在线观看| 免费电影在线观看免费观看| 脱女人内裤的视频| 国产亚洲精品综合一区在线观看| 亚洲经典国产精华液单 | 最后的刺客免费高清国语| 久久草成人影院| 亚洲午夜理论影院| 午夜福利18| 国产精品爽爽va在线观看网站| 精品人妻熟女av久视频| 久久久久久久久中文| av天堂中文字幕网| 美女黄网站色视频| 欧美色视频一区免费| 久久人妻av系列| 国产精品乱码一区二三区的特点| 一a级毛片在线观看| 一区二区三区四区激情视频 | 精品久久久久久久久亚洲 | 亚洲欧美精品综合久久99| 搡老熟女国产l中国老女人| 亚洲 国产 在线| 国产黄a三级三级三级人| 欧美成人a在线观看| 特大巨黑吊av在线直播| 国产精品1区2区在线观看.| 中文字幕人成人乱码亚洲影| 精品久久久久久久末码| 男女那种视频在线观看| 国产高清视频在线播放一区| 综合色av麻豆| 久久精品国产99精品国产亚洲性色| 国产亚洲精品综合一区在线观看| 美女免费视频网站| 啪啪无遮挡十八禁网站| 国产综合懂色| 久久热精品热| av在线天堂中文字幕| 亚洲人成伊人成综合网2020| 久久久久久久久中文| 91午夜精品亚洲一区二区三区 | 欧美xxxx性猛交bbbb| 99热精品在线国产| 午夜福利在线在线| 中文字幕熟女人妻在线| 极品教师在线视频| 在线a可以看的网站| 男女做爰动态图高潮gif福利片| 丝袜美腿在线中文| 别揉我奶头 嗯啊视频| 老熟妇乱子伦视频在线观看| 最新中文字幕久久久久| 亚洲欧美清纯卡通| 日日摸夜夜添夜夜添小说| 日韩欧美 国产精品| 国产欧美日韩精品亚洲av| 在现免费观看毛片| 免费av观看视频| 99国产极品粉嫩在线观看| 噜噜噜噜噜久久久久久91| x7x7x7水蜜桃| 亚洲国产精品合色在线| 宅男免费午夜| 精品人妻熟女av久视频| 久久久久久久久中文| 最好的美女福利视频网| 免费人成视频x8x8入口观看| 观看免费一级毛片| 欧美+亚洲+日韩+国产| 久久久精品欧美日韩精品| 国产日本99.免费观看| 在线播放无遮挡| 日本五十路高清| 99国产综合亚洲精品| 欧美黄色淫秽网站| 亚洲专区国产一区二区| 欧美黑人欧美精品刺激| 很黄的视频免费| 91九色精品人成在线观看| 自拍偷自拍亚洲精品老妇| 国产黄色小视频在线观看| 国产精品伦人一区二区| 国产精品亚洲一级av第二区| 国产一区二区亚洲精品在线观看| 欧美中文日本在线观看视频| 日本成人三级电影网站| 床上黄色一级片| 国产成人av教育| 国产午夜精品论理片| 欧美在线黄色| 在线观看午夜福利视频| 啦啦啦韩国在线观看视频| 国产免费一级a男人的天堂| 色哟哟·www| 国产成人av教育| 国产精品99久久久久久久久| 久久中文看片网| 国产淫片久久久久久久久 | 88av欧美| 欧美色欧美亚洲另类二区| 国产精品99久久久久久久久| 国产91精品成人一区二区三区| 国产亚洲av嫩草精品影院| 亚洲三级黄色毛片| 久久国产乱子伦精品免费另类| 亚洲av第一区精品v没综合| 中文字幕免费在线视频6| 波野结衣二区三区在线| 亚洲内射少妇av| x7x7x7水蜜桃| 日本成人三级电影网站| 欧美一区二区国产精品久久精品| 国产精品亚洲av一区麻豆| 欧美日韩国产亚洲二区| 他把我摸到了高潮在线观看| 真人一进一出gif抽搐免费| 高潮久久久久久久久久久不卡| 嫩草影视91久久| 哪里可以看免费的av片| 欧美3d第一页| av福利片在线观看| 国产一区二区三区在线臀色熟女| 特大巨黑吊av在线直播| 日韩中文字幕欧美一区二区| 精品午夜福利在线看| 亚洲人与动物交配视频| 欧美不卡视频在线免费观看| 亚洲av免费在线观看| 能在线免费观看的黄片| 黄色配什么色好看| 亚洲av美国av| 老司机午夜福利在线观看视频| 国产黄片美女视频| 久久人人爽人人爽人人片va | 亚洲成人中文字幕在线播放| 给我免费播放毛片高清在线观看| 亚洲一区二区三区不卡视频| 色吧在线观看| 国产69精品久久久久777片| 757午夜福利合集在线观看| 丰满人妻一区二区三区视频av| 成人av在线播放网站| 91九色精品人成在线观看| 九九久久精品国产亚洲av麻豆| 亚洲五月婷婷丁香| 国产伦在线观看视频一区| 18禁黄网站禁片免费观看直播| 欧美xxxx黑人xx丫x性爽| 国内久久婷婷六月综合欲色啪| 全区人妻精品视频| 偷拍熟女少妇极品色| 女人被狂操c到高潮| 国产精品国产高清国产av| 国产一区二区在线观看日韩| 日韩亚洲欧美综合| 日韩精品中文字幕看吧| 亚洲av免费在线观看| bbb黄色大片| 99国产精品一区二区三区| 久久中文看片网| 午夜激情福利司机影院| 精品不卡国产一区二区三区| 久久久精品欧美日韩精品| 亚洲精品乱码久久久v下载方式| 午夜两性在线视频| 亚洲av一区综合| 老司机福利观看| 三级国产精品欧美在线观看| 日本a在线网址| 久久久久精品国产欧美久久久| 12—13女人毛片做爰片一| 俄罗斯特黄特色一大片| 精品久久久久久久久久免费视频| 无遮挡黄片免费观看| 中国美女看黄片| 亚洲av中文字字幕乱码综合| 国产欧美日韩一区二区三| 午夜精品久久久久久毛片777| 一本综合久久免费| 一二三四社区在线视频社区8| 精品无人区乱码1区二区| 激情在线观看视频在线高清| 色哟哟·www| 如何舔出高潮| 久久精品人妻少妇| 99riav亚洲国产免费| 国产精品电影一区二区三区| 99久久精品热视频| 日本黄色片子视频| ponron亚洲| 免费搜索国产男女视频| netflix在线观看网站| 欧美午夜高清在线| 亚洲精品456在线播放app | 丰满人妻一区二区三区视频av| 午夜激情欧美在线| 99视频精品全部免费 在线| 国产精品三级大全| 亚洲国产精品久久男人天堂| 欧美又色又爽又黄视频| 欧美乱妇无乱码| 欧美乱色亚洲激情| 精品乱码久久久久久99久播| 91麻豆av在线| 亚洲av成人不卡在线观看播放网| 一个人观看的视频www高清免费观看| 看免费av毛片| www.色视频.com| 欧美xxxx性猛交bbbb| 变态另类成人亚洲欧美熟女| 超碰av人人做人人爽久久| 老女人水多毛片| 亚洲精品日韩av片在线观看| 亚洲人成电影免费在线| 国产精品国产高清国产av| 色哟哟·www| 别揉我奶头 嗯啊视频| 91午夜精品亚洲一区二区三区 | 五月玫瑰六月丁香| 国产人妻一区二区三区在| 免费一级毛片在线播放高清视频| 国产爱豆传媒在线观看| 天美传媒精品一区二区| 99久久99久久久精品蜜桃| 一进一出好大好爽视频| 蜜桃亚洲精品一区二区三区| 亚洲 国产 在线| 午夜老司机福利剧场| 在线观看舔阴道视频| 亚洲精品色激情综合| 人妻久久中文字幕网| 久久热精品热| 真人一进一出gif抽搐免费| 国产精品国产高清国产av| 九九热线精品视视频播放| 亚洲成a人片在线一区二区| 免费一级毛片在线播放高清视频| 久久热精品热| 国产精品精品国产色婷婷| 亚洲欧美精品综合久久99| 亚洲一区二区三区色噜噜| 我要搜黄色片| 91在线观看av| 国产黄a三级三级三级人| www.熟女人妻精品国产| 婷婷色综合大香蕉| 精品无人区乱码1区二区| 久久香蕉精品热| 亚洲av第一区精品v没综合| 免费一级毛片在线播放高清视频| www.999成人在线观看| 老司机福利观看| 久久欧美精品欧美久久欧美| 国产真实乱freesex| 麻豆成人av在线观看| 国产精品女同一区二区软件 | 全区人妻精品视频| 日本免费a在线| 99精品在免费线老司机午夜| 精品乱码久久久久久99久播| 精品人妻一区二区三区麻豆 | 欧美不卡视频在线免费观看| 18禁在线播放成人免费| 中文字幕高清在线视频| 在线免费观看的www视频| 精品午夜福利视频在线观看一区| 亚洲精品一区av在线观看| 欧美成狂野欧美在线观看| 激情在线观看视频在线高清| 亚洲七黄色美女视频| 赤兔流量卡办理| 国产伦精品一区二区三区视频9| 在现免费观看毛片| 国产午夜精品论理片| 色吧在线观看| 麻豆av噜噜一区二区三区| 伊人久久精品亚洲午夜| 日韩中字成人| 老鸭窝网址在线观看| 欧美极品一区二区三区四区| 国产乱人伦免费视频| 黄片小视频在线播放| 能在线免费观看的黄片| 国产伦在线观看视频一区| 亚洲美女视频黄频| 给我免费播放毛片高清在线观看| 天堂av国产一区二区熟女人妻| 色5月婷婷丁香| 亚洲精品成人久久久久久| www.www免费av| 成人美女网站在线观看视频| 亚洲精品在线观看二区| 亚洲成人久久爱视频| 国产精品久久久久久久电影| 一本久久中文字幕| 一区二区三区激情视频| 青草久久国产| 亚洲中文字幕一区二区三区有码在线看| 色播亚洲综合网| 日韩国内少妇激情av| 综合色av麻豆| 国产在线男女| 久久精品综合一区二区三区| 小蜜桃在线观看免费完整版高清| 国产一区二区在线av高清观看| 淫妇啪啪啪对白视频| 黄色女人牲交| 中文字幕人妻熟人妻熟丝袜美| 精品国产三级普通话版| 亚洲av电影不卡..在线观看| 在线观看一区二区三区| 亚洲不卡免费看| 少妇丰满av| 国产精品亚洲av一区麻豆| 人人妻,人人澡人人爽秒播| www.www免费av| 欧美日韩国产亚洲二区| 97超视频在线观看视频| 成人性生交大片免费视频hd| 久久久久久久久久成人| 国产免费av片在线观看野外av| 好男人在线观看高清免费视频| 国产免费男女视频| 91麻豆精品激情在线观看国产| 三级毛片av免费| 亚洲av二区三区四区| 九色成人免费人妻av| 亚洲第一欧美日韩一区二区三区| 全区人妻精品视频| 深爱激情五月婷婷| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久亚洲精品不卡| 国产真实乱freesex| 成年版毛片免费区| 成人特级av手机在线观看| 99久国产av精品| 国产精品综合久久久久久久免费| 国产精品美女特级片免费视频播放器| 悠悠久久av| 久久久色成人| 免费在线观看影片大全网站| 久久九九热精品免费| .国产精品久久| 制服丝袜大香蕉在线| 国产在线男女| 色综合站精品国产| 国产熟女xx| 免费大片18禁| 两性午夜刺激爽爽歪歪视频在线观看| 欧美激情在线99| 少妇高潮的动态图| 亚洲国产精品成人综合色| 国产精品一区二区三区四区免费观看 | 国产精品一区二区三区四区免费观看 | 国产中年淑女户外野战色| 久久久久久久久久成人| 免费电影在线观看免费观看| 男人狂女人下面高潮的视频| 全区人妻精品视频| 一个人免费在线观看电影| 窝窝影院91人妻| 久久精品久久久久久噜噜老黄 | 欧美乱色亚洲激情| 两个人视频免费观看高清| 一个人看视频在线观看www免费| 搡老妇女老女人老熟妇| 国产精品国产高清国产av| 91久久精品国产一区二区成人| 亚洲色图av天堂| 亚洲黑人精品在线| 午夜两性在线视频| 别揉我奶头~嗯~啊~动态视频| 国内少妇人妻偷人精品xxx网站| 国产真实乱freesex| 国内少妇人妻偷人精品xxx网站| 国产美女午夜福利| 亚洲乱码一区二区免费版| 亚洲一区高清亚洲精品| 在线观看一区二区三区| 一本精品99久久精品77| 天堂动漫精品| 少妇的逼水好多| 90打野战视频偷拍视频| 老司机深夜福利视频在线观看| 岛国在线免费视频观看| 亚洲经典国产精华液单 | 亚洲精品亚洲一区二区| 成人鲁丝片一二三区免费| 一级a爱片免费观看的视频| 欧美另类亚洲清纯唯美| 宅男免费午夜| 免费人成在线观看视频色| 亚洲av日韩精品久久久久久密| 亚洲国产日韩欧美精品在线观看| 国产黄色小视频在线观看| 国产毛片a区久久久久| 国产精品久久久久久人妻精品电影| 老熟妇乱子伦视频在线观看| 午夜福利视频1000在线观看| 国产乱人视频| 成年女人毛片免费观看观看9| 欧美在线黄色| 老女人水多毛片| 国产国拍精品亚洲av在线观看| 成人精品一区二区免费| 深爱激情五月婷婷| 亚洲av二区三区四区| 亚洲欧美清纯卡通| 嫩草影视91久久| 欧美黄色淫秽网站| 久久婷婷人人爽人人干人人爱| 真实男女啪啪啪动态图| 麻豆国产97在线/欧美| 中文资源天堂在线| 99久久无色码亚洲精品果冻| 亚洲中文字幕一区二区三区有码在线看| 亚洲人成电影免费在线| 日韩成人在线观看一区二区三区| 亚洲精品乱码久久久v下载方式| 国产精品免费一区二区三区在线| 天堂av国产一区二区熟女人妻| av在线天堂中文字幕| 又黄又爽又免费观看的视频| 日韩高清综合在线| 91久久精品电影网| 国产单亲对白刺激| 久久久久久久久大av| 国产成人a区在线观看| 丰满乱子伦码专区| 欧美+日韩+精品| 69人妻影院| 国产精品国产高清国产av| 日韩高清综合在线| 亚洲 欧美 日韩 在线 免费| 久久久国产成人免费| 国产乱人视频| 成人av在线播放网站| 国产一级毛片七仙女欲春2| 女人十人毛片免费观看3o分钟| 国产 一区 欧美 日韩| 国产精品一区二区三区四区久久| netflix在线观看网站| 国产伦精品一区二区三区四那| 亚洲熟妇熟女久久| 好看av亚洲va欧美ⅴa在| 亚洲国产精品成人综合色| 日本黄色视频三级网站网址| 国产中年淑女户外野战色| 日韩中字成人| 亚洲在线观看片| 久久久精品欧美日韩精品| 成人国产综合亚洲| 久久性视频一级片| 人人妻人人澡欧美一区二区| 色尼玛亚洲综合影院| 国产av不卡久久| 国产伦在线观看视频一区| 日本 欧美在线| 欧美日韩综合久久久久久 | 久久精品国产99精品国产亚洲性色| 国产伦精品一区二区三区视频9| 午夜精品一区二区三区免费看| 中文字幕熟女人妻在线| 婷婷精品国产亚洲av| 久9热在线精品视频| 人人妻,人人澡人人爽秒播| 欧美另类亚洲清纯唯美| 久久精品影院6| 亚洲,欧美,日韩| 成人无遮挡网站| 99久久久亚洲精品蜜臀av| 国产精品亚洲av一区麻豆| 免费av毛片视频| 看片在线看免费视频| av视频在线观看入口| 99视频精品全部免费 在线| bbb黄色大片| 国产一区二区在线观看日韩| 高潮久久久久久久久久久不卡| 极品教师在线免费播放| 国产精品亚洲一级av第二区| x7x7x7水蜜桃| 久久久成人免费电影| 国产伦精品一区二区三区四那| 69人妻影院| 亚洲中文字幕一区二区三区有码在线看| 欧美色欧美亚洲另类二区| 国产伦精品一区二区三区视频9| 欧美潮喷喷水| 91在线精品国自产拍蜜月| 动漫黄色视频在线观看| 性欧美人与动物交配| 宅男免费午夜| 91狼人影院| 国产精品久久视频播放| 此物有八面人人有两片| 亚洲男人的天堂狠狠| 十八禁网站免费在线| 久久久久久久久中文| 性插视频无遮挡在线免费观看| 我的老师免费观看完整版| 色视频www国产| 色综合亚洲欧美另类图片| 精品久久久久久久久久免费视频| 成年女人永久免费观看视频| 看十八女毛片水多多多| 可以在线观看的亚洲视频| 亚洲欧美日韩无卡精品| 真实男女啪啪啪动态图| 乱人视频在线观看| 亚洲经典国产精华液单 | 婷婷色综合大香蕉| 亚洲三级黄色毛片| 久久6这里有精品| 成人av一区二区三区在线看| 欧美极品一区二区三区四区| 亚洲精品粉嫩美女一区| 网址你懂的国产日韩在线| 国产麻豆成人av免费视频| 搡女人真爽免费视频火全软件 | 精品无人区乱码1区二区| 国产伦精品一区二区三区视频9| 精品久久久久久久人妻蜜臀av| 毛片一级片免费看久久久久 | 日本 av在线| 成人一区二区视频在线观看| 亚洲黑人精品在线| 97碰自拍视频| 午夜久久久久精精品| 九九热线精品视视频播放| 在线观看66精品国产| 日本撒尿小便嘘嘘汇集6| 高清日韩中文字幕在线| 亚洲国产日韩欧美精品在线观看| 国产成人影院久久av| 亚洲美女黄片视频| 日本一本二区三区精品| 少妇裸体淫交视频免费看高清| 国产一级毛片七仙女欲春2| 91麻豆av在线| 亚洲国产欧美人成| 一级黄色大片毛片| aaaaa片日本免费| 一进一出好大好爽视频| 一本综合久久免费| 亚洲欧美日韩高清专用| 精华霜和精华液先用哪个| 日韩大尺度精品在线看网址| а√天堂www在线а√下载| 国产精品综合久久久久久久免费| 久久精品国产亚洲av天美| or卡值多少钱| 成年女人看的毛片在线观看| 色综合欧美亚洲国产小说| 久久久久久国产a免费观看| av天堂中文字幕网| 国产69精品久久久久777片| 国产精品一区二区三区四区久久| 黄色一级大片看看| 国产精品美女特级片免费视频播放器| 久久久成人免费电影| 好看av亚洲va欧美ⅴa在| 一本一本综合久久| 久久精品夜夜夜夜夜久久蜜豆| 免费一级毛片在线播放高清视频| 99久久无色码亚洲精品果冻| 嫩草影视91久久| av黄色大香蕉| 偷拍熟女少妇极品色| 麻豆久久精品国产亚洲av| 免费黄网站久久成人精品 | 又黄又爽又免费观看的视频| 亚洲国产色片| www.色视频.com| 91九色精品人成在线观看| 好男人在线观看高清免费视频| 精品国产三级普通话版| 久久久久免费精品人妻一区二区| 国产精华一区二区三区| 97人妻精品一区二区三区麻豆| 欧美日韩亚洲国产一区二区在线观看| 国产国拍精品亚洲av在线观看| av福利片在线观看| 国产精品伦人一区二区| 亚洲av日韩精品久久久久久密| 亚洲av二区三区四区| 在线看三级毛片| 久久人人爽人人爽人人片va | 国产精品自产拍在线观看55亚洲| 亚洲中文字幕日韩| 极品教师在线视频| 少妇的逼好多水| 一个人看视频在线观看www免费| 国产精华一区二区三区| 国内揄拍国产精品人妻在线| 国产亚洲精品综合一区在线观看| 亚洲欧美日韩高清专用|