沈 威，陶孟侖，陳定方，劉紅俊，李鵬輝，明廷鑫

(武漢理工大學(xué) 物流工程學(xué)院，湖北 武漢430063)

近年來，無線傳感網(wǎng)絡(luò)被廣泛應(yīng)用于軍事、智能交通、環(huán)境監(jiān)控、醫(yī)療衛(wèi)生等領(lǐng)域，但是網(wǎng)絡(luò)節(jié)點中的傳感器壽命受到了其電源壽命的限制，傳統(tǒng)的化學(xué)電池的壽命有限，電池的更換工作量大，提高了無線網(wǎng)絡(luò)的維護成本，而利用環(huán)境振動能量為無線網(wǎng)絡(luò)傳感器供電成為解決該供電問題的一個有效途徑［1-2］。將振動能量轉(zhuǎn)換為電能有以下3 種常見模式：電磁式、壓電式和靜電式［3］，其中，壓電式能量收集器憑借其尺寸小、易加工、機電能量轉(zhuǎn)換率高和易與電路集成等優(yōu)點，更適合為無線網(wǎng)絡(luò)傳感器供電。

傳統(tǒng)的壓電能量收集器是單懸臂梁形式的，文獻［4 -6］設(shè)計了雙穩(wěn)態(tài)壓電能量收集器，在一定程度上提高了輸出功率。文獻［7 -8］設(shè)計了壓電陣列式能量收集器，與單懸臂壓電能量收集器相比，輸出功率提高了數(shù)倍。但是，以上研究并沒有從機電能量轉(zhuǎn)換率的角度對壓電能量收集器進行分析和改進，且在輸出功率上還存在提高空間。

筆者在線性的壓電陣列式能量收集器基礎(chǔ)上，引入了非線性磁力，建立了非線性壓電陣列式能量收集器的機電耦合模型，并將模型無量綱化，通過Matlab 對線性和非線性的能量收集器模型分別施加白噪聲激勵，對比分析其輸出功率、機電能量轉(zhuǎn)化率。

1 非線性壓電陣列式能量收集器模型

1.1 基本結(jié)構(gòu)

圖1、圖2 所示為筆者設(shè)計的非線性壓電陣列式能量收集器(nonlinear piezoelectric beam array vibration energy harvester，NPBAVEH)，其有3個相同的懸臂梁陣列式排布，懸臂梁根部粘結(jié)有壓電片和電路部分。懸臂梁基底材料為銅，壓電片材料為PZT，PZT 因具有壓電常數(shù)大和靈敏度高的特點［9］，成為壓電能量收集器中應(yīng)用最為廣泛的壓電陶瓷。壓電片工作模式為d31 模式，d31模式下應(yīng)力的方向與電場的方向相互垂直。3 個懸臂梁的端部粘結(jié)有同一個磁鐵塊A，與另外一個固定在基座上的磁鐵塊B產(chǎn)生非線性磁力，B對A的磁力在x軸方向上的分量為FM;磁鐵A在x軸方向上的位移為X;A與B在y軸方向的中心距為d。如圖3 所示，將3 個壓電片串聯(lián)，其輸出電壓比并聯(lián)高。

圖1 非線性壓電陣列式能量收集器結(jié)構(gòu)示意圖1

圖2 非線性壓電陣列式能量收集器結(jié)構(gòu)示意圖2

圖3 非線性壓電陣列式能量收集器的電路示意圖

1.2 機電耦合模型

外部振動激勵通過基座傳遞到懸臂梁根部，由于3 個形狀、結(jié)構(gòu)相同的懸臂梁粘結(jié)同一個磁鐵塊，因此，3 個懸臂梁振動形式一致，振動響應(yīng)相同，只要求出其中一個壓電懸臂梁的電壓輸出，就可以得出3 個壓電懸臂梁的總電壓輸出。在簡諧激勵下，求出第i個懸臂梁在一階模態(tài)下的機電耦合模型(i=1，2，3)：

式中：X為磁鐵A的垂直位移;Vi為第i個壓電懸臂梁的輸出電壓;m、c、k分別為單個懸臂梁的等效質(zhì)量、機械阻尼系數(shù)和等效剛度;m=M/3 +33/140m1，其中M為磁鐵A的質(zhì)量，m1為單個懸臂梁的質(zhì)量;Cp為壓電片的等效電容;R為外部負載;每個懸臂梁所受的磁力在x軸方向上都是相同的，因此Fmi=FM/3;κc、κV均為機電耦合系數(shù)。

為了便于仿真，進行無量綱變換，令x=X/l，則式(1)和式(2)的無量綱形式分別為：

現(xiàn)實中環(huán)境振動往往是隨機的，因此為了模擬實際情況，在研究NPBAVEH 振動響應(yīng)時，輸入的幅值服從標準正態(tài)分布高斯白噪聲，即n(t)～N(0，1)，它是壓電發(fā)電系統(tǒng)研究中最常使用的隨機激勵。則式(3)可變?yōu)椋?/p>

其中，f為無量綱的激勵幅值。

因為3 個壓電片串聯(lián)，所以系統(tǒng)總的輸出電壓為：

表1 所示為部分物理參數(shù)，筆者所取參數(shù)值主要參照文獻［10 -11］。其中，L為懸臂梁長度;hb、hp分別為懸臂梁和壓電片的厚度;le為壓電片的長度;l為磁鐵A中心到梁根部的水平距離;MA、MB分別為磁鐵A、B的磁化強度;VA、VB分別為磁鐵A、B的體積。

表1 物理參數(shù)表

2 數(shù)值計算與仿真分析

2.1 總勢能分析

根據(jù)式(4)和式(5)，計算NPBAVEH 系統(tǒng)中的a、b、γ、μ、ζ、? 的值分別為1.20、0.30、0.05、0.05、0.01、0.50。

能量收集器系統(tǒng)的總勢能包括系統(tǒng)的彈性勢能和磁鐵在磁場中的勢能，可表示為：

a、b取不同值時，系統(tǒng)的總勢能圖如圖4 所示。當(dāng)a=1.2，b=0.3 時，NPBAVEH 系統(tǒng)勢能函數(shù)有兩個穩(wěn)定的奇點和一個不穩(wěn)定的奇點，在兩個穩(wěn)定的奇點附近出現(xiàn)兩個勢阱。當(dāng)FM=0時，相當(dāng)于無磁力耦合作用，即可把能量收集器系統(tǒng)看成線性的壓電陣列式能量收集器(piezoelectric beam array vibration energy harvester，PBAVEH)，此時a= -1.5，b=0，勢能函數(shù)有一個穩(wěn)定的奇點，且只有一個勢阱。

圖4 總勢能圖

2.2 輸出功率的對比分析

算得輸入激勵-fn(t)的方差為σ2，通過Matlab 對式(4)和式(5)進行仿真，得到隨σ2變化的輸出功率圖，如圖5 所示。當(dāng)σ2從小至大取不同值時，對NPBAVEH(a= 1. 2、b= 0. 3)和PBAVEH(a=1.5、b=0)輸出電流的平均功率進行對比，當(dāng)σ2＜1.1 時，NPBAVEH 和PBAVEH 輸出的功率相差較為近似，而當(dāng)σ2＞1. 1 時，NPBAVEH 的輸出功率明顯大于PBAVEH。

圖5 隨σ2 變化的輸出功率圖

2.3 機電能量轉(zhuǎn)換率的對比分析

能量收集器的機電能量轉(zhuǎn)換率可以表示為：

式中：無量綱參數(shù)pe、pm分別為輸出電流的平均功率、輸入機械能的平均功率分別為其相應(yīng)的瞬時功率。根據(jù)式(4)～式(6)的計算式為：

通過Matlab 仿真得到機電能量轉(zhuǎn)換率隨σ2變化的曲線圖，如圖6 所示。PBAVEH 的機電能量轉(zhuǎn)換率基本不隨σ2的改變而改變;而NPBAVEH 的機電能量轉(zhuǎn)換率隨著σ2的變大而逐漸增加，最后趨于平緩。當(dāng)σ2＜1.1 時，PBAVEH 的機電能量轉(zhuǎn)換率較高;而當(dāng)σ2＞1.1 時，NPBAVEH的機電能量轉(zhuǎn)換率較高。

圖6 隨σ2 變化的機電能量轉(zhuǎn)換率圖

取σ2=0.5 時，NPBAVEH 與PBAVEH 無量綱響應(yīng)x的時域圖如圖7 所示。由于激勵的能量較小，NPBAVEH 的響應(yīng)還在單個勢阱中做小幅振蕩。取σ2=2.5 時，NPBAVEH 與PBAVEH 無量綱響應(yīng)x的時域圖如圖8 所示，NPBAVEH 的響應(yīng)不斷在兩個勢阱中躍遷振蕩。無論σ2取多少，PBAVEH 的響應(yīng)都只是在單個勢阱中，而隨著σ2的增大，激勵的能量輸入越來越高，使得NPBAVEH 的響應(yīng)從單個勢阱中做小幅振蕩逐漸變?yōu)樵趦蓚€勢阱中大幅度躍遷振蕩，最終使得NPBAVEH 機電能量轉(zhuǎn)換率逐漸升高。

圖7 σ2 =0.5 時NPBAVEH 與PBAVEH 的響應(yīng)x 時域圖

圖8 σ2 =2.5 時NPBAVEH 與PBAVEH 的響應(yīng)x 時域圖

3 結(jié)論

為了提高壓電能量收集器的輸出性能，筆者設(shè)計了一種新型非線性壓電陣列式能量收集器，該能量收集器將壓電懸臂陣列排布，并將懸臂梁的輸出電壓串聯(lián)。增加兩個永磁鐵對結(jié)構(gòu)進行改進，建立系統(tǒng)的機電耦合模型，通過輸入白噪聲激勵，并從輸出功率和機電能量轉(zhuǎn)換率的角度，與線性壓電陣列式能量收集器進行仿真對比。仿真結(jié)果表明，在激勵強度較大時，相對于線性壓電陣列式能量收集器，所設(shè)計的該壓電能量收集器有效提高了輸出功率和機電能量轉(zhuǎn)換率。

［1］ BOZIDAR M.Piezoelectric wide -bandwidth vibration energy harvester-smart Sand［D］.Yale：Yale University，2011.

［2］ 劉祥建，陳仁文.壓電振動能量收集裝置研究現(xiàn)狀及發(fā)展趨勢［J］.振動與沖擊，2012，31(16)：169-176.

［3］ 王青萍.基于壓電懸臂梁的振動能量收集器的研究［D］.武漢：華中科技大學(xué)，2011.

［4］ ERTURK A，INMAN D J.Broadband piezoelectric power generation on high-energy orbits of the bistable duffing oscillator with electromechanical coupling［J］.Journal of Sound and Vibration，2011，330(10)：2339-2353.

［5］ 崔巖，王飛，董維杰，等.非線性壓電式能量采集器［J］.光學(xué)精密工程，2012，20(12)：2737 -2743.

［6］ STANTON S C，OWENS B A M，MANN B P. Harmonic balance analysis of the bistable piezoelectric inertial generator［J］. Journal of Sound and Vibration，2012，331(15)：3617 -3627.

［7］ 佘引，溫志渝，趙興強，等.MEMS 壓電陣列振動能量收集器［J］.傳感技術(shù)學(xué)報，2014，27(8)：1033-1037.

［8］ WEN Z Y，DENG L C，ZHAO X Q，et al. Improving voltage output with PZT beam array for MEMS-based vibration energy harvester：theory and experiment［J］.Microsystem Technologies，2015，21(2)：331 -339.

［9］ 陳定方，孔建益，楊家軍，等.現(xiàn)代機械設(shè)計師手冊［M］.北京：機械工業(yè)出版社，2014：1326 -1328.

［10］ 孫舒.雙穩(wěn)態(tài)壓電懸臂梁發(fā)電系統(tǒng)的動力學(xué)特性研究［D］.天津：天津大學(xué)，2013.

［11］ 高毓璣，冷永剛，范勝波，等.彈性支撐雙穩(wěn)壓電懸臂梁振動響應(yīng)及能量采集研究［J］. 物理學(xué)報，2014，63(9)：62 -74.