柱棒式超磁致伸縮振動(dòng)能量收集裝置建模與實(shí)驗(yàn)

何忠波 柏 果 鄭佳偉 周景濤 劉國(guó)平

(1.陸軍工程大學(xué)車(chē)輛與電氣工程系， 石家莊 050003； 2.32184部隊(duì)， 北京 100072)

0 引言

超磁致伸縮材料(Giant magnetostrictive material， GMM)是一種新型磁功能智能材料，可實(shí)現(xiàn)機(jī)械能-磁場(chǎng)能-電能之間的高效能量轉(zhuǎn)換，與其他智能材料相比，其具有輸出力大、應(yīng)變大、滯環(huán)小、響應(yīng)速度快、能量密度高等優(yōu)點(diǎn)[1-3]。目前GMM在致動(dòng)器、換能器和傳感器等領(lǐng)域已經(jīng)得到了廣泛應(yīng)用，為超磁致伸縮能量收集相關(guān)研究奠定了一定的理論基礎(chǔ)[4-6]。振動(dòng)能量收集是指利用特定能量的技術(shù)，對(duì)于緩解能源壓力以及提升能源的利用率具有一定意義[7-9]。超磁致伸縮式能量收集裝置(Giant magnetostrictive energy harvesting device，GMEHD)基于GMM的維拉里效應(yīng)，將外界的機(jī)械振動(dòng)能轉(zhuǎn)換為可儲(chǔ)存使用的電能，具有輸出功率高、工作穩(wěn)定性強(qiáng)、應(yīng)用適應(yīng)性廣等優(yōu)勢(shì)，可為微機(jī)電系統(tǒng)甚至是常用的中小型功率器件供電，具有較好的發(fā)展?jié)摿10-13]。從20世紀(jì)80年代開(kāi)始，國(guó)外研究機(jī)構(gòu)嘗試采集振動(dòng)能量，并且論證了振動(dòng)能量回收應(yīng)用于微功耗器件驅(qū)動(dòng)的可行性。

目前，針對(duì)超磁致伸縮致動(dòng)器的研究很全面，體系也已經(jīng)很成熟，但對(duì)能量收集方面研究并不多。國(guó)內(nèi)外學(xué)者致力于研究實(shí)用性更強(qiáng)的應(yīng)用型新結(jié)構(gòu)，然而對(duì)材料非線性、機(jī)電耦合特性等考慮較少，現(xiàn)有模型的研究大多集中于GMM棒磁化模型，并未遵循能量收集的相關(guān)工作過(guò)程；對(duì)GMEHD的有限元仿真分析也十分有限，部分仿真的假設(shè)過(guò)于簡(jiǎn)單，與材料本征不夠吻合，仿真結(jié)論不可靠。本文根據(jù)機(jī)械振動(dòng)通常具有周期性往復(fù)的實(shí)際特點(diǎn)，設(shè)計(jì)一種柱棒式GMEHD，理論上棒式GMEHD進(jìn)行能量回收可采集任意頻率下的振動(dòng)能，而且其對(duì)體積要求低，可做成閉合磁路、避免磁泄漏，有效提高材料的機(jī)電轉(zhuǎn)換效率；并且，棒式能量回收結(jié)構(gòu)與磁致伸縮棒式致動(dòng)器結(jié)構(gòu)幾乎一樣，材料的磁機(jī)耦合特性也完全相同，由于致動(dòng)器研究的成熟，該類(lèi)型結(jié)構(gòu)具有較好的研究基礎(chǔ)[14]。為準(zhǔn)確計(jì)算能量收集產(chǎn)生的電動(dòng)勢(shì)，結(jié)合其實(shí)際工作過(guò)程，建立系統(tǒng)的輸入輸出模型，并通過(guò)有限元仿真分析以及實(shí)驗(yàn)研究對(duì)模型準(zhǔn)確性進(jìn)行驗(yàn)證及對(duì)比分析，以期為GMEHD結(jié)構(gòu)優(yōu)化和實(shí)際應(yīng)用奠定基礎(chǔ)。

1 基本結(jié)構(gòu)與原理

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

設(shè)計(jì)了圖1所示的柱棒式GMEHD，核心元件為GMM棒；永磁體作為偏置磁鐵為GMM棒提供偏置磁場(chǎng)，能夠使GMEHD工作于機(jī)電耦合系數(shù)較高的線性區(qū)間[15-16]，同時(shí)將偏置磁鐵上下對(duì)置，相對(duì)周向放置減小了一定體積；由于GMM棒抗拉強(qiáng)度低，脆性較大，同時(shí)為了便于GMEHD接受外界激勵(lì)，在實(shí)際應(yīng)用中需要設(shè)計(jì)承力桿[17]；預(yù)緊彈簧為GMM棒提供預(yù)壓力，通過(guò)調(diào)整蓋帽和套筒之間的螺紋可調(diào)整預(yù)緊力；蓋帽、套筒和支座均采用高導(dǎo)磁材料，整個(gè)裝置構(gòu)成閉合磁路，實(shí)現(xiàn)輸出效率最大化。

圖1 能量收集裝置結(jié)構(gòu)圖Fig.1 Structure diagram of energy harvesting device1.承力桿 2.預(yù)緊彈簧 3.偏置磁鐵 4.采集線圈 5.力傳感器預(yù)留位置 6.底部引線槽 7.固定螺栓孔 8.支座 9.超磁致伸縮棒 10.側(cè)方引線槽 11.套筒 12.蓋帽

1.2 工作原理

在機(jī)械應(yīng)力作用下材料產(chǎn)生應(yīng)變，會(huì)引起GMM棒內(nèi)部磁化狀態(tài)發(fā)生變化，這種現(xiàn)象稱(chēng)為磁致伸縮逆效應(yīng)，又稱(chēng)維拉里效應(yīng)。

根據(jù)維拉里效應(yīng)，如圖2所示，GMM棒一端固定，一端受力，外界壓力通過(guò)承力桿和預(yù)緊彈簧作用到GMM棒上。GMM棒形狀發(fā)生變化，其縱向長(zhǎng)度減小使其內(nèi)部磁疇會(huì)發(fā)生偏轉(zhuǎn)，從而使GMM棒的內(nèi)部磁場(chǎng)分布發(fā)生改變，其磁通量發(fā)生了變化。由法拉第電磁感應(yīng)定律，如果GMM棒外部繞有感應(yīng)線圈，就會(huì)產(chǎn)生感應(yīng)電動(dòng)勢(shì)并可通過(guò)線圈進(jìn)行收集[18-19]。

圖2 GMEHD工作原理圖Fig.2 Working principles of GMEHD

2 系統(tǒng)建模

能量收集是將機(jī)械能轉(zhuǎn)換為磁場(chǎng)能最終轉(zhuǎn)換為電勢(shì)能的過(guò)程，故模型需建立輸入力和輸出電動(dòng)勢(shì)之間的關(guān)系。相較于大多數(shù)現(xiàn)有模型不同，本模型充分考慮了能量收集工作整個(gè)過(guò)程，從GMEHD受到外界激勵(lì)到GMM棒產(chǎn)生形變，再到GMM棒內(nèi)部產(chǎn)生變化的磁場(chǎng)，最后通過(guò)線圈產(chǎn)生感應(yīng)電動(dòng)勢(shì)，模型更貼近裝置的實(shí)際工作過(guò)程。建模過(guò)程如圖3所示。

圖3 GMEHD建模流程圖Fig.3 Flow chart of GMEHD modeling

GMM棒在偏置磁場(chǎng)作用下的振動(dòng)過(guò)程和磁化情況比較復(fù)雜，所以對(duì)裝置進(jìn)行建模時(shí)作如下假設(shè)：①GMM棒的橫截面在振動(dòng)時(shí)始終保持為平面作整體運(yùn)動(dòng)，并略去其縱向伸縮引起的橫向變形。②GMM棒受到的應(yīng)力σ、產(chǎn)生的壓應(yīng)變?chǔ)?、偏置磁?chǎng)強(qiáng)度H0和產(chǎn)生的磁場(chǎng)強(qiáng)度H均勻分布。③將GMM棒等效為質(zhì)量-彈簧-阻尼模型，并假設(shè)其等效阻尼系數(shù)c和彈性模量E不隨外界力的變化而改變。④忽略預(yù)緊彈簧對(duì)系統(tǒng)輸入力的影響，GMM棒能夠完全得到外界力的輸入。

2.1 動(dòng)力學(xué)模型

基于上述假設(shè)，將GMEHD簡(jiǎn)化為單自由度振動(dòng)系統(tǒng)，如圖4所示。

圖4 GMEHD等效模型Fig.4 GMEHD equivalent model

取GMM棒中心平衡點(diǎn)為坐標(biāo)原點(diǎn)，x軸正方向豎直向下，由圖4的二階線性有阻尼質(zhì)量-彈簧系統(tǒng)，得到系統(tǒng)的動(dòng)力學(xué)微分方程

(1)

式中m——GMM棒等效質(zhì)量

k——GMM棒等效剛度

x——GMM棒位移

由GMM棒的位移可得到其應(yīng)變?yōu)?/p>

(2)

式中l(wèi)——GMM棒長(zhǎng)度

當(dāng)對(duì)GMM棒施加壓應(yīng)力時(shí)，其磁疇將主要沿垂直軸線的易磁化方向分布。此時(shí)，GMM棒的應(yīng)變?chǔ)伺c磁化強(qiáng)度M的關(guān)系近似為基于能量基礎(chǔ)的二次疇轉(zhuǎn)模型

(3)

式中λs——GMM棒的飽和磁致伸縮應(yīng)變

Ms——GMM棒的飽和磁化強(qiáng)度

2.2 磁場(chǎng)模型

在外力作用下，GMM產(chǎn)生伸縮應(yīng)變，進(jìn)而導(dǎo)致其內(nèi)部有效磁場(chǎng)發(fā)生變化。由熱動(dòng)力學(xué)基本定理可得棒內(nèi)的Helmholtz自由能密度為

(4)

式中μ0——真空磁導(dǎo)率

H0——偏置磁場(chǎng)強(qiáng)度

α——疇壁相互作用系數(shù)

σ——GMM棒在外力作用下的應(yīng)力

T——內(nèi)部溫度S——熵

Helmholtz自由能密度J關(guān)于磁化強(qiáng)度M的微分即為有效磁場(chǎng)，式(4)對(duì)M微分，得

(5)

式(3)兩邊對(duì)M微分，得

(6)

將式(6)代入式(5)中，得

(7)

(8)

式中A——橫截面積

2.3 電動(dòng)勢(shì)模型

考慮GMM固有的磁滯非線性，由電磁學(xué)原理可得棒內(nèi)磁感應(yīng)強(qiáng)度為

B=μ0(M+H)

(9)

在得到GMM棒內(nèi)磁感應(yīng)強(qiáng)度B后，由法拉第電磁感應(yīng)定律可計(jì)算產(chǎn)生的感應(yīng)電動(dòng)勢(shì)為

(10)

式中Φ——GMM棒內(nèi)部的磁通量

n——線圈匝數(shù)

將式(9)代入式(10)，可得

(11)

將式(3)、(8)代入式(11)，得

(12)

將式(2)代入式(12)，得

(13)

聯(lián)立式(1)、(13)即可得到在一定輸入力作用下GMEHD所產(chǎn)生的電動(dòng)勢(shì)。

2.4 模型求解

在外界環(huán)境中，振動(dòng)往往具有周期性或階段周期性，通過(guò)傅里葉變換，可將非正弦的周期性外力分解為一系列周期性正弦力的線性組合，因此以正弦力為輸入，借助式(1)、(13)，計(jì)算系統(tǒng)的輸出電動(dòng)勢(shì)。

為了確保正弦力在不同時(shí)刻均能夠完全作用到GMEHD上，即為了使正弦力始終為正，預(yù)應(yīng)力應(yīng)大于或等于正弦力的幅值。在預(yù)應(yīng)力和正弦力的共同作用下，令外力F=F0+F1sin(ωt)，則式(1)為

(14)

式中F0——預(yù)應(yīng)力

F1——正弦力幅值

ω——正弦力頻率

由振動(dòng)理論可知，單自由度強(qiáng)迫振動(dòng)微分方程的全部解包括穩(wěn)態(tài)振動(dòng)與瞬態(tài)振動(dòng)。穩(wěn)態(tài)振動(dòng)是外界激勵(lì)下的強(qiáng)迫振動(dòng)，在正弦激勵(lì)下為持續(xù)等幅振動(dòng)，瞬態(tài)振動(dòng)在阻尼情況下為衰減振動(dòng)，只在振動(dòng)開(kāi)始后的一段時(shí)間才有意義，一般情況下可以不考慮。分開(kāi)討論穩(wěn)態(tài)振動(dòng)預(yù)應(yīng)力和正弦力。

預(yù)應(yīng)力F0為施加在GMM棒上的靜力，其產(chǎn)生的位移x1為靜位移，不會(huì)有速度和加速度的產(chǎn)生，所以x1為

(15)

正弦力F1sin(ωt)為施加在GMM棒上的簡(jiǎn)諧力，由單自由度強(qiáng)迫振動(dòng)對(duì)簡(jiǎn)諧激勵(lì)的響應(yīng)結(jié)果可以得到動(dòng)位移x2為

x2=Xsin(ωt-θ)

(16)

(17)

(18)

故GMM棒的總位移x為

(19)

將式(19)代入式(13)中，得到GMEHD在正弦力下的輸出電動(dòng)勢(shì)為

(20)

(21)

通過(guò)式(21)可以看出，在一定預(yù)應(yīng)力和正弦力共同作用下，GMEHD的輸出電動(dòng)勢(shì)幅值和輸入正弦力的頻率呈線性關(guān)系；在外力及GMEHD一定的情況下，除了簡(jiǎn)諧量外均為常數(shù)，其中正弦量不影響輸出電動(dòng)勢(shì)u的相位，其相位與輸入正弦力相比，超前π/2。

2.5 線圈優(yōu)化

2.5.1線圈結(jié)構(gòu)

GMEHD采用多層繞組圓柱形線圈，其結(jié)構(gòu)如圖5所示。線圈骨架為空心圓筒，r1為內(nèi)徑，r2為外徑，h為骨架高度，銅漆包線由銅芯構(gòu)成，半徑為r，繞組高度為D(即r2-r1)。

圖5 多層圓柱形線圈示意圖Fig.5 Schematic of multilayered cylindrical coil

結(jié)合多層繞組圓柱形線圈的結(jié)構(gòu)以及電阻公式，每一匝線圈的平均長(zhǎng)度為π(r1+r2)，則整個(gè)線圈的電阻R為

(22)

式中ρ——銅芯電阻率

在線圈骨架結(jié)構(gòu)一定的情況下，線圈匝數(shù)N與銅芯半徑r會(huì)受到結(jié)構(gòu)尺寸的約束，即

(23)

在線圈的實(shí)際纏繞過(guò)程中，不可避免會(huì)在漆包線之間產(chǎn)生間隔，實(shí)際的線圈匝數(shù)達(dá)不到理想計(jì)算結(jié)果，所以對(duì)式(23)添加修正系數(shù)η(0<η<1)，即

(24)

消去r，得

(25)

在線圈骨架一定的情況下，r1、r2、D、h均為定值，同時(shí)在銅漆包線材質(zhì)一定的情況下，電阻率ρ也為定值，因此線圈電阻R與匝數(shù)平方成正比，令

(26)

則

R=C1n2

(27)

在線圈高h(yuǎn)大于繞組高度D時(shí)，內(nèi)含磁心的多層繞組圓柱形線圈的電感L為[20]

(28)

式中μr——GMM棒相對(duì)磁導(dǎo)率

W——(r1+r2)/h的函數(shù)

V——h/D的函數(shù)

在線圈骨架一定的情況下，這三者均可視為常數(shù)，因此線圈自感L與匝數(shù)平方成正比，令

(29)

則

L=C2n2

(30)

2.5.2電路分析

圖6 線圈等效電路Fig.6 Coil equivalent circuit

由于GMEHD產(chǎn)生交流電動(dòng)勢(shì)，在分析線圈對(duì)輸出電動(dòng)勢(shì)的影響時(shí)，除了考慮線圈的直流阻抗外，還需考慮其感抗對(duì)輸出電動(dòng)勢(shì)的影響。圖6是線圈的等效電路。其中，u是GMEHD產(chǎn)生的交流電動(dòng)勢(shì)，i為電路中的電流，R2為外接負(fù)載電阻，u2為負(fù)載電阻上的電壓。將線圈等效為電阻R和電感L的串聯(lián)，u1為線圈上的電壓，ur、uL分別是電阻R和電感L上的電壓。為簡(jiǎn)化求解過(guò)程，對(duì)電路進(jìn)行分析時(shí)采用相量圖作為輔助分析工具，根據(jù)相量圖的幾何關(guān)系便可對(duì)電路進(jìn)行相關(guān)運(yùn)算。得到圖6的相量圖如圖7所示。

圖7 電路相量圖Fig.7 Phasor diagram of circuit

以I為參考相量，外接負(fù)載電阻R兩端電壓U2與I同相，線圈上的電壓U1含有兩個(gè)分量：電阻r上的電壓Ur也與I同相；電感L上的電壓UL較I超前90°。由相量圖得(U2，U1，U)和(Ur，UL，U1)2個(gè)電壓三角形，由相量圖應(yīng)用余弦定理計(jì)算

(31)

其中

XL=ωL=2πfL

(32)

式中φ——線圈電路相位角

XL——線圈感抗

由式(27)、(30)～(32)可以解得

(33)

從式(33)可以看出，φ不受線圈匝數(shù)n的影響。由式(21)可知，在外界輸入力一定以及GMM棒結(jié)構(gòu)參數(shù)確定的情況下，GMEHD產(chǎn)生的電動(dòng)勢(shì)與匝數(shù)n呈線性關(guān)系

U=c1n

(34)

式中c1——比例常數(shù)

則整個(gè)電路的有功功率P為

(35)

在外接負(fù)載電阻R2一定時(shí)，P只與線圈匝數(shù)n相關(guān)，將式(35)變形得

(36)

(37)

2.5.3線圈匝數(shù)

根據(jù)GMEHD以及使用的GMM棒的結(jié)構(gòu)尺寸，確定線圈各參數(shù)如表1所示。

表1 線圈各參數(shù)Tab.1 Parameter in coil

線圈各項(xiàng)參數(shù)確定后，即可確定出W與V取值為5與0.87。由式(37)可知，除了線圈結(jié)構(gòu)參數(shù)與外接電阻之外，輸入力頻率f也會(huì)影響線圈匝數(shù)。圖8是匝數(shù)n與輸入力頻率f的關(guān)系曲線。

圖8 線圈匝數(shù)與輸入力頻率的關(guān)系曲線Fig.8 Relationship curve of coil turns and input force frequency

3 仿真

利用有限元軟件COMSOL Multiphysics建立有限元模型，設(shè)置模型相關(guān)組件，對(duì)GMEHD進(jìn)行仿真分析。

3.1 有限元模型

圖9 幾何模型Fig.9 Geometric model

建立GMEHD幾何模型，如圖9所示。由于GMEHD為軸對(duì)稱(chēng)結(jié)構(gòu)，無(wú)需建立GMEHD的實(shí)際三維空間結(jié)構(gòu)，僅建立了GMEHD的二維平面模型，在保證不影響計(jì)算結(jié)果精度的情況下，簡(jiǎn)化了模型，減小了計(jì)算量。

3.2 模型參數(shù)設(shè)置

對(duì)幾何模型各部分添加不同的材料，材料參數(shù)如表2所示。

表2 有限元模型各參數(shù)Tab.2 Parameter in finite element model

物理場(chǎng)選擇包含固體力學(xué)場(chǎng)和磁場(chǎng)的耦合場(chǎng)，在GMEHD底端設(shè)置固定約束，頂桿上端設(shè)置預(yù)應(yīng)力以及相應(yīng)載荷。網(wǎng)格剖分選擇由用戶控制網(wǎng)格，對(duì)不同區(qū)域選擇不同網(wǎng)格形狀和大小，在保證精度的情況下減小計(jì)算量，網(wǎng)格剖分結(jié)果如圖10所示。在完成相關(guān)設(shè)置后進(jìn)行計(jì)算，即可得到仿真分析的結(jié)果。

3.3 仿真結(jié)果

將1 N預(yù)應(yīng)力以及幅值1 N、頻率1 Hz的正弦力共同作用到承力桿上，可得GMEHD的磁感應(yīng)強(qiáng)度分布圖，圖11為一個(gè)正弦周期內(nèi)不同時(shí)間節(jié)點(diǎn)的磁感應(yīng)強(qiáng)度分布圖。

圖10 網(wǎng)格剖分圖Fig.10 Mesh split diagram

圖11 磁感應(yīng)強(qiáng)度分布圖Fig.11 Magnetic induction distribution diagrams

通過(guò)圖11可以看出，隨著時(shí)間的變化，在GMEHD中磁感應(yīng)強(qiáng)度變化率最大的是GMM棒，基于此可通過(guò)GMM的逆磁效應(yīng)得到感生電動(dòng)勢(shì)。

4 實(shí)驗(yàn)

為驗(yàn)證模型的有效性，制作了GMEHD樣機(jī)(力傳感器在試驗(yàn)時(shí)置于GMEHD內(nèi)部)，GMM棒采用定向凝固法制備，輸出桿、端蓋、外殼和底座均采用高磁導(dǎo)率的15號(hào)低碳鋼，各結(jié)構(gòu)參數(shù)均與有限元仿真設(shè)定值一致，樣機(jī)實(shí)物如圖12所示。

圖12 GMEHD和力傳感器Fig.12 GMEHD and force sensor

4.1 實(shí)驗(yàn)系統(tǒng)

圖13 實(shí)驗(yàn)系統(tǒng)Fig.13 Experiment system1.函數(shù)信號(hào)發(fā)生器 2.模態(tài)激振器 3.固定支架 4.電荷放大器 5.數(shù)據(jù)采集器 6.功率放大器 7.計(jì)算機(jī)

用于測(cè)試GMEHD輸出電動(dòng)勢(shì)實(shí)驗(yàn)系統(tǒng)如圖13所示。實(shí)驗(yàn)過(guò)程為：DC1022U型函數(shù)信號(hào)發(fā)生器產(chǎn)生輸入信號(hào)，信號(hào)經(jīng)PA020型功率放大器放大后輸入到JZ020型模態(tài)激振器中，模態(tài)激振器通過(guò)輸出頂桿將力傳遞到GMEHD，GMM棒受到外力后產(chǎn)生電動(dòng)勢(shì)；將GMEHD中力傳感器測(cè)得的輸入力信號(hào)經(jīng)PE01AD型電荷放大器后接入到SA1604型數(shù)據(jù)采集器中的一個(gè)通道中，同時(shí)將GMEHD產(chǎn)生的電動(dòng)勢(shì)接入到數(shù)據(jù)采集器的另一通道中，再將數(shù)據(jù)采集器接入到計(jì)算機(jī)中，計(jì)算機(jī)端便可通過(guò)數(shù)據(jù)采集器的配套軟件進(jìn)行實(shí)驗(yàn)測(cè)量與數(shù)據(jù)記錄。

4.2 實(shí)驗(yàn)測(cè)量

由于GMM材料抗壓不抗拉，同時(shí)為了使GMM棒在正弦力輸入下能夠不出現(xiàn)伸長(zhǎng)的情況，在實(shí)驗(yàn)前通過(guò)調(diào)整蓋帽和套筒之間的螺紋給予GMEHD一定的預(yù)壓力，保證GMM棒工作在合適的變形范圍內(nèi)，并防止其受到拉應(yīng)力而發(fā)生斷裂。在實(shí)驗(yàn)中以不同頻率單位正弦力為激勵(lì)，通過(guò)置于GMM棒下的力傳感器及數(shù)據(jù)采集器，可以得到輸入力-時(shí)間曲線，如圖14所示。

圖14 不同輸入力頻率的力幅值曲線Fig.14 Time curves of different frequency input force

將線圈得到的電動(dòng)勢(shì)信號(hào)接入到數(shù)據(jù)采集器的另一通道中，可得到在不同頻率正弦力下的電動(dòng)勢(shì)隨時(shí)間變化的曲線，如圖15所示。

圖15 不同輸入力頻率下的輸出電動(dòng)勢(shì)曲線Fig.15 Output electromotive force curves under different frequency forces

圖17 不同輸入力頻率下的輸出電動(dòng)勢(shì)曲線對(duì)比Fig.17 Comparison of output electromotive force curves under different frequency forces

由圖15可知，在輸入正弦力頻率只有1 Hz時(shí)，電動(dòng)勢(shì)幅值在2 mV左右，因幅值較小，所以受振動(dòng)干擾影響較大，輸出波形噪聲較大，但仍能夠看出波形形狀；當(dāng)頻率達(dá)到10 Hz時(shí)，輸出波形較為穩(wěn)定且清晰，電動(dòng)勢(shì)幅值在15 mV左右；當(dāng)頻率達(dá)到50 Hz時(shí)，電動(dòng)勢(shì)幅值能夠達(dá)到75 mV，輸出電動(dòng)勢(shì)幅值與輸入正弦力頻率基本呈線性關(guān)系；圖16是輸入力10 Hz時(shí)與輸出電動(dòng)勢(shì)的對(duì)比，可以看出輸入力頻率不影響GMEHD輸出電動(dòng)勢(shì)的相位，其始終超前于輸入力相位π/2，原因是線圈為電感元件，而在電感元件電路中，在相位上電流始終比電壓滯后π/2，且電流與磁場(chǎng)強(qiáng)度、輸入力的相位一致，會(huì)使得輸出電動(dòng)勢(shì)相位超前輸入力π/2。這均與理論分析結(jié)果一致。

圖16 10 Hz時(shí)輸入正弦力與輸出電動(dòng)勢(shì)對(duì)比Fig.16 Comparison of input sinusoidal force and output electromotive force at 10 Hz

4.3 對(duì)比分析

將實(shí)驗(yàn)所得電動(dòng)勢(shì)曲線分別與模型計(jì)算值和有限元仿真結(jié)果進(jìn)行對(duì)比，如圖17所示。由圖17可以看出，在頻率為1 Hz時(shí)因受振動(dòng)噪聲影響較大，造成實(shí)驗(yàn)結(jié)果比模型計(jì)算值和仿真結(jié)果大，在頻率為10、50 Hz時(shí)，實(shí)驗(yàn)結(jié)果比模型計(jì)算值和仿真結(jié)果略小，因?yàn)樵谀Ｐ陀?jì)算和有限元仿真分析時(shí)，沒(méi)有考慮從GMEHD輸入頂桿到GMM棒這一碰撞過(guò)程中的輸入力能量損耗，而在實(shí)驗(yàn)時(shí)這一損耗實(shí)際存在，從而在一定程度上影響輸出電動(dòng)勢(shì)的幅值，造成實(shí)驗(yàn)結(jié)果比模型計(jì)算值和有限元仿真結(jié)果略小。以輸入力頻率10 Hz時(shí)為例，模型計(jì)算與仿真得到的輸出電動(dòng)勢(shì)幅值為17 mV，實(shí)驗(yàn)所得輸出電動(dòng)勢(shì)幅值為15 mV，相對(duì)誤差為11.7%，由于碰撞損耗不可避免，誤差在合理范圍之內(nèi)，不影響實(shí)驗(yàn)對(duì)模型準(zhǔn)確性的驗(yàn)證。

5 結(jié)論

(1)設(shè)計(jì)了一種柱棒式超磁致伸縮振動(dòng)能量收集裝置(GMEHD)，整個(gè)裝置磁路高度閉合，輸出效率較大，能量密度高。

(2)利用振動(dòng)理論、電磁場(chǎng)理論建立了GMEHD的輸出電動(dòng)勢(shì)模型，推導(dǎo)了輸出電動(dòng)勢(shì)的計(jì)算式；利用仿真軟件COMSOL Multiphysics，建立了GMEHD有限元模型，研究了其內(nèi)磁感應(yīng)強(qiáng)度以及輸出電動(dòng)勢(shì)，得到了兩者的變化規(guī)律。

(3)通過(guò)電路分析，對(duì)一定尺寸骨架下的線圈匝數(shù)進(jìn)行了優(yōu)化分析，得到了最優(yōu)線圈匝數(shù)與輸入力頻率的關(guān)系曲線。

(4)搭建了實(shí)驗(yàn)測(cè)試系統(tǒng)，對(duì)GMEHD進(jìn)行了不同頻率的正弦激勵(lì)實(shí)驗(yàn)。結(jié)果表明，GMEHD輸出電動(dòng)勢(shì)的相位不受輸入力頻率的影響，其始終超前于輸入力相位π/2；在一定頻率范圍內(nèi)，輸入力頻率與輸出電動(dòng)勢(shì)幅值呈線性關(guān)系。輸出波形與模型計(jì)算值、有限元仿真結(jié)果吻合較好，誤差在合理范圍內(nèi)，驗(yàn)證了模型的準(zhǔn)確性。