龔曙光，蔣自強(qiáng)，謝桂蘭，張建平，盧海山

(湘潭大學(xué) 機(jī)械工程學(xué)院，湘潭 411105)

0 引 言

風(fēng)能是當(dāng)今世界使用最為廣泛、清潔的能源之一，其利用方式也多種多樣，風(fēng)力發(fā)電機(jī)是風(fēng)能最普通的使用方式。近年來，渦激共振壓電的研究是風(fēng)能利用的一個(gè)熱點(diǎn)[2]，但是該方式常用于微型發(fā)電，具有發(fā)電量小的局限性。而傳統(tǒng)帶槳葉的風(fēng)力機(jī)具有運(yùn)輸、制造、維修困難等問題且會(huì)危及飛行的鳥類?；诖?，本文提出一種新型的風(fēng)能利用裝置——無(wú)葉片風(fēng)力發(fā)電機(jī)。與傳統(tǒng)有葉片的風(fēng)力機(jī)相比，無(wú)葉片風(fēng)力機(jī)具有制造成本低、結(jié)構(gòu)簡(jiǎn)單、易加工和維修方便等優(yōu)點(diǎn)，同時(shí)也為未來分布式發(fā)電帶來活力。

永磁發(fā)電機(jī)相對(duì)于勵(lì)磁發(fā)電機(jī)具有更高的發(fā)電效率和節(jié)能效果，其具有結(jié)構(gòu)簡(jiǎn)單且可靠性、穩(wěn)定性高等優(yōu)點(diǎn)[3-4]。如文獻(xiàn)[5]提出了一種新型的無(wú)鐵心的直線發(fā)電機(jī)，探討了其在定常速度下的電磁特性，并提高了系統(tǒng)穩(wěn)定性，但其忽略了發(fā)電機(jī)的電壓幅值和輸出功率；文獻(xiàn)[6-7]研究了永磁直線發(fā)電機(jī)在勻速下的電壓、功率、邊端力等電機(jī)特性。除此以外，許多學(xué)者也致力于研究直線發(fā)電機(jī)在簡(jiǎn)諧運(yùn)動(dòng)下的電磁特性。如文獻(xiàn)[8-9]推導(dǎo)了調(diào)制型永磁直線發(fā)電機(jī)在正弦速度下的電動(dòng)勢(shì)表達(dá)式，并應(yīng)用于波浪能發(fā)電；文獻(xiàn)[10]給出了單相發(fā)電機(jī)在正弦下的運(yùn)動(dòng)方程和發(fā)電機(jī)的輸出方程，通過實(shí)驗(yàn)測(cè)得在最大負(fù)載運(yùn)行下發(fā)電機(jī)的輸出效率可達(dá)83.92 %；文獻(xiàn)[11]通過對(duì)永磁直線發(fā)電機(jī)在簡(jiǎn)諧運(yùn)動(dòng)下的理論推導(dǎo)，找出了適應(yīng)電機(jī)性能的振動(dòng)區(qū)間，并得到了仿真驗(yàn)證；文獻(xiàn)[12]探討了直線發(fā)電機(jī)在自由活塞式運(yùn)動(dòng)時(shí)的電磁特性，比較了矩形截面和工字形截面的永磁體對(duì)發(fā)電機(jī)性能的影響；文獻(xiàn)[13]通過實(shí)驗(yàn)分析了發(fā)電機(jī)在空載和負(fù)載條件下的電機(jī)輸出特性。文獻(xiàn)[14]、文獻(xiàn)[15]研究了永磁直線發(fā)電機(jī)在不規(guī)則正弦運(yùn)動(dòng)時(shí)的電磁特性，推導(dǎo)了波浪對(duì)浮子的作用力，并得到了驗(yàn)證。

目前永磁直線發(fā)電機(jī)主要針對(duì)勻速運(yùn)動(dòng)或簡(jiǎn)諧運(yùn)動(dòng)開展研究，盡管也對(duì)不規(guī)則的正弦運(yùn)動(dòng)有所探討，但其應(yīng)用場(chǎng)景主要是用于波浪能發(fā)電，而波浪能的重力波周期一般為5～10 s(海岸工程)。相反，無(wú)葉片風(fēng)力機(jī)捕能柱的風(fēng)致渦激振動(dòng)是一種高頻低振幅振動(dòng)，且其振動(dòng)也呈現(xiàn)為非周期性?；跓o(wú)葉片風(fēng)力機(jī)捕能柱的振動(dòng)特性，本文首先對(duì)永磁直線發(fā)電機(jī)進(jìn)行了結(jié)構(gòu)設(shè)計(jì)，并對(duì)其空載特性、負(fù)載特性進(jìn)行仿真分析，采用反正切與正弦函數(shù)的乘積逼近發(fā)電機(jī)在渦激共振下的運(yùn)動(dòng)軌跡，以探討捕能柱的來流風(fēng)速和振幅對(duì)發(fā)電機(jī)電磁特性的影響。

1 無(wú)葉片風(fēng)力機(jī)原理

無(wú)葉片風(fēng)力發(fā)電機(jī)主要由捕能柱、轉(zhuǎn)換裝置與發(fā)電機(jī)組成，其原理簡(jiǎn)圖如圖1所示。捕能柱在風(fēng)致渦激振動(dòng)下產(chǎn)生橫向擺動(dòng)，轉(zhuǎn)換裝置則將捕能柱的擺動(dòng)轉(zhuǎn)換成直線運(yùn)動(dòng)再傳給直線發(fā)電機(jī)，永磁直線發(fā)電機(jī)的功能是將機(jī)械能轉(zhuǎn)化為電能，為了方便安裝與維修，其位于風(fēng)力發(fā)電機(jī)的最下端。

圖1 無(wú)葉片風(fēng)力發(fā)電機(jī)原理圖

在風(fēng)致渦激振動(dòng)中，無(wú)葉片風(fēng)力機(jī)捕能柱獲得的風(fēng)能輸入功率：

(1)

式中:ρ為空氣的密度；U為風(fēng)速；D為柱體直徑；H為柱體高度。

設(shè)捕能柱擺動(dòng)的最大角度為θmax，穩(wěn)定擺動(dòng)的角度為θ1，經(jīng)轉(zhuǎn)化后其最大行程為L(zhǎng)max，穩(wěn)定時(shí)行程為L(zhǎng)。本文假定無(wú)葉片風(fēng)力機(jī)的設(shè)計(jì)風(fēng)速為6.0 m/s，其柱體直徑為0.2 m，高為1.3 m。由式(1)即可得到風(fēng)能的輸入功率為36.3 W。

2 圓筒形永磁直線發(fā)電機(jī)

2.1 電機(jī)結(jié)構(gòu)

永磁直線發(fā)電機(jī)的結(jié)構(gòu)如圖2所示，因?yàn)榘l(fā)電機(jī)的次級(jí)是作上下運(yùn)動(dòng)，為減輕次級(jí)質(zhì)量，方便電能傳輸，將發(fā)電機(jī)結(jié)構(gòu)設(shè)計(jì)成外部長(zhǎng)初級(jí)、內(nèi)部短次級(jí)。永磁直線發(fā)電機(jī)由三相初級(jí)鐵心、三相繞組、次級(jí)鐵心和次級(jí)永磁體組成。為了方便制造和獲得較大的氣隙磁密，永磁體采用表貼式，且沿著徑向方向充磁，相鄰永磁體充磁方向相反。三相繞組采用模塊化餅式繞組，單獨(dú)加工成型后嵌入定子槽內(nèi)。

圖2 永磁直線發(fā)電機(jī)二維結(jié)構(gòu)圖

此外，發(fā)電機(jī)的運(yùn)行效率和成本是電機(jī)設(shè)計(jì)過程中需要考慮的重要因素，而行程的長(zhǎng)短又與發(fā)電機(jī)的運(yùn)行效率和成本有直接關(guān)系。因此，結(jié)合額定風(fēng)速下捕能柱在渦激共振時(shí)的振幅，以及其輸入功率為36.3 W，在假定設(shè)計(jì)電壓為12 V，通過磁路設(shè)計(jì)，得到發(fā)電機(jī)的主要參數(shù)，如表1所示。

表1 永磁直線發(fā)電機(jī)的參數(shù)

2.2 勻速下電機(jī)空載特性

通過仿真分析得到發(fā)電機(jī)空載時(shí)的磁通線分布如圖3所示，圖4為發(fā)電機(jī)空載磁感應(yīng)強(qiáng)度的云圖。

圖3 空載磁通

圖4 空載磁密

從圖3可看到，磁力線主回路從N極發(fā)出，經(jīng)過氣隙進(jìn)入定子齒，再經(jīng)過定子軛后從下一個(gè)定子齒流出并再一次進(jìn)入氣隙，然后通過相鄰極性相反的磁極，最后從次級(jí)軛部回到原來的磁極。即除了磁力線主回路之外，存在著另外兩條磁路，一條是從永磁體N極出發(fā)，直接從氣隙進(jìn)入到相鄰的下一個(gè)磁極，再通過次級(jí)軛回到原永磁體；另一條是從N極出發(fā)，經(jīng)過初級(jí)齒后流入到初級(jí)槽，再?gòu)某跫?jí)槽流入氣隙，然后通過相鄰磁極和次級(jí)軛回到該永磁體的S極。

從圖4可知，初級(jí)齒部平均磁感應(yīng)強(qiáng)度為1.09 T，在齒頂處的兩端其磁感應(yīng)強(qiáng)度較大；初級(jí)軛部平均磁感應(yīng)強(qiáng)度為0.55 T，在定子齒尖端處存在極小區(qū)域內(nèi)的磁密過大，但是飽和程度不大，所以從總體來看可忽略磁飽和現(xiàn)象，符合設(shè)計(jì)要求。同時(shí)，由于發(fā)電機(jī)的能量轉(zhuǎn)換是在氣隙中完成，故氣隙的長(zhǎng)度和磁場(chǎng)強(qiáng)度的變化規(guī)律對(duì)發(fā)電機(jī)的性能有著重要影響。

圖5顯示了發(fā)電機(jī)在t=0時(shí)徑向氣隙磁密的分布。由于發(fā)電機(jī)的開斷影響，氣隙磁密在永磁體N極和S極不夠?qū)ΨQ，氣隙磁密的最大值為0.92 T，平均磁密為0.56 T；正對(duì)永磁體兩端的氣隙磁密高于永磁體中間部分，這是因?yàn)樵谝欢ǚ秶鷥?nèi)，永磁體極面面積越大，則極面中心的磁感應(yīng)強(qiáng)度越低于其周邊的磁感應(yīng)強(qiáng)度。

圖5 t=0時(shí)刻空載氣隙徑向磁密

對(duì)氣隙磁密進(jìn)行傅里葉分解，其諧波分布如圖6所示。從圖6可知，氣隙磁密的諧波次數(shù)較多，氣隙磁密的基波幅值為0.788 7 T，波形畸變率為17.5%。

圖6 空載氣隙磁密FFT

發(fā)電機(jī)的感應(yīng)電動(dòng)勢(shì)如圖7所示。從圖7可知，三相電壓幅值分別為18.57 V、18.38 V、18.41 V，與設(shè)計(jì)的電壓值基本吻合，證明所設(shè)計(jì)的電機(jī)是合理的。通過對(duì)電壓諧波分析可知，其三相基波幅值分別為18.50 V、17.82 V、17.90 V，平均波形畸變率為7.5%。

圖7 感應(yīng)電動(dòng)勢(shì)

2.3 負(fù)載特性

發(fā)電機(jī)的負(fù)載特性是衡量發(fā)電機(jī)性能的關(guān)鍵因素。本文中的負(fù)載工況為純電阻電路，在次級(jí)速度為0.5 m/s下計(jì)算了不同負(fù)載工況下的電壓、電流、功率和效率。

圖8顯示了發(fā)電機(jī)輸出功率、效率與負(fù)載的變化曲線。從圖8可知，輸出功率隨負(fù)載的增大呈現(xiàn)先增大后減小的趨勢(shì)，這是因?yàn)殡S著負(fù)載的增大，負(fù)載端電壓逐漸增加，當(dāng)負(fù)載繼續(xù)增大時(shí)，電壓增加會(huì)趨于平緩，而由電路理論可知，當(dāng)電壓保持不變時(shí)，純電阻電路中的輸出功率隨負(fù)載的增大而減小。同時(shí)在純電阻電路中，電機(jī)的損耗主要是銅耗和鐵耗，當(dāng)電阻增大時(shí)，電機(jī)電流下降，繞組銅耗也隨之降低，而鐵耗在低頻電磁場(chǎng)中所占比例不大，因此電機(jī)的效率隨著外電阻電路的阻值增大而增大。圖8的最大輸出功率為114.5 W，最大效率為89.77%。

圖8 輸出功率、效率隨負(fù)載的變化

發(fā)電機(jī)的負(fù)載電壓如圖9所示。從圖9可知，在純阻性負(fù)載狀態(tài)下，其三相電壓幅值分別為17.91 V、17.96 V、17.88 V，相比空載感應(yīng)電動(dòng)勢(shì)分別小了0.58 V、0.42 V、0.53 V，且三相的對(duì)稱性良好。

圖9 負(fù)載電壓

由于外電路采用了純阻性負(fù)載，故其電流波形相位與電壓一致，數(shù)值上符合歐姆定律。為了獲得較好的電機(jī)輸出性能，兼顧功率和效率，現(xiàn)選取外電阻為12 Ω，圖10顯示了在電阻值為12 Ω時(shí)，其電壓、功率和效率隨次級(jí)速度的變化。

圖10 電壓、功率和效率隨次級(jí)速度變化

從圖10可看出，電壓與輸出功率隨次級(jí)速度的增加而增加，這是由于在純電阻電路中阻值不變時(shí)，當(dāng)速度增大其輸出電壓增大，輸出功率也隨之增大，且電壓與速度為線性關(guān)系，而輸出功率與速度的二次方成正比，因此輸出功率的增加趨勢(shì)如同拋物線。同時(shí)，電機(jī)的效率隨著速度的增大逐漸趨近于88%，這是因?yàn)殡S著速度的增大，電機(jī)的損耗也逐漸增大。

通過計(jì)算得到本文所設(shè)計(jì)發(fā)電機(jī)的體積功率密度為59.71 kW/m3，該值與文獻(xiàn)[8]實(shí)驗(yàn)測(cè)試值57.27 kW/m3相接近，這進(jìn)一步說明本文所設(shè)計(jì)的直線發(fā)電機(jī)是可行的。

3 風(fēng)致渦激共振的影響

捕能柱在來流風(fēng)作用下會(huì)產(chǎn)生渦激振動(dòng)，隨著來流風(fēng)速的增大，其渦激頻率會(huì)逐漸接近系統(tǒng)固有頻率，而使捕能柱產(chǎn)生共振，此時(shí)捕能柱的橫向擺動(dòng)角度(也可稱為振幅)達(dá)到最大值，其風(fēng)速區(qū)間也稱為“鎖頻區(qū)間”。然而在自然界中，風(fēng)速一般處于不穩(wěn)定狀態(tài)，捕能柱的橫向擺動(dòng)也處于一個(gè)非穩(wěn)定的簡(jiǎn)諧運(yùn)動(dòng)，因此本文主要針對(duì)來流風(fēng)速和振幅對(duì)發(fā)電機(jī)性能的影響進(jìn)行分析。

3.1 風(fēng)速的影響

本文無(wú)葉片風(fēng)力機(jī)的設(shè)計(jì)風(fēng)速為6 m/s，在 4～8 m/s間選取多個(gè)風(fēng)速，通過計(jì)算得到不同風(fēng)速下次級(jí)速度的變化曲線，如圖11所示。

圖11 不同風(fēng)速下次級(jí)速度變化曲線

從圖11中可看到，隨著來流風(fēng)速的增大，次級(jí)速度幅值出現(xiàn)了先增大后減少的趨勢(shì)，其變化呈現(xiàn)出了不規(guī)則的簡(jiǎn)諧運(yùn)動(dòng)，直到來流風(fēng)速達(dá)到設(shè)計(jì)風(fēng)速時(shí)，次級(jí)速度幅值達(dá)到最大，同時(shí)在“鎖頻”風(fēng)速區(qū)間，次級(jí)的上下振動(dòng)頻率也趨于穩(wěn)定。

圖12給出了永磁直線發(fā)電機(jī)在設(shè)計(jì)風(fēng)速下的瞬時(shí)最大輸出功率和平均輸出功率。從圖12可看到，當(dāng)風(fēng)速為6 m/s時(shí)，瞬時(shí)最大輸出功率可達(dá)35.62 W，平均功率達(dá)到15.95 W。

圖12 設(shè)計(jì)風(fēng)速下發(fā)電機(jī)的輸出功率

按與圖12類似的計(jì)算方法，可得到不同風(fēng)速下發(fā)電機(jī)的平均輸出功率與效率，如圖13所示。

圖13 不同風(fēng)速下發(fā)電機(jī)輸出功率和效率

圖13中，發(fā)電機(jī)的輸出功率與效率均隨風(fēng)速的增大呈現(xiàn)先增大后下降的趨勢(shì)，特別是當(dāng)來流風(fēng)速處于共振區(qū)間時(shí)，發(fā)電機(jī)的輸出功率和效率達(dá)到最大，即其最大平均輸出功率為15.95 W，最大效率為85.44%。這也意味著當(dāng)風(fēng)速處于共振區(qū)間時(shí)，本文所設(shè)計(jì)的發(fā)電機(jī)具有良好的輸出特性。

同時(shí)，發(fā)電機(jī)三相電壓隨風(fēng)速的變化如圖14所示。從圖14可知，發(fā)電機(jī)在捕能柱鎖頻區(qū)間內(nèi)具有較高的電壓，當(dāng)捕能柱處在共振狀態(tài)時(shí)，其具有較大的擺動(dòng)頻率和振幅(擺動(dòng)角度)，使得發(fā)電機(jī)次級(jí)的運(yùn)動(dòng)速度更大。同時(shí)由于次級(jí)運(yùn)動(dòng)速度不恒定，發(fā)電機(jī)的次級(jí)運(yùn)動(dòng)規(guī)律和磁場(chǎng)分布匹配不佳，并且受發(fā)電機(jī)兩端開斷的影響，發(fā)電機(jī)的三相電壓不對(duì)稱，此時(shí)發(fā)電機(jī)的三相電壓的幅值分別是5.03 V、10.30 V、7.75 V。當(dāng)風(fēng)速值離開共振區(qū)間后，捕能柱處于高頻低幅擺動(dòng)，故電壓下降。

圖14 不同風(fēng)速下發(fā)電機(jī)輸出電壓

3.2 振幅(擺動(dòng)角度)的影響

設(shè)振幅倍數(shù)以設(shè)計(jì)風(fēng)速下的次級(jí)運(yùn)動(dòng)行程為基頻，當(dāng)發(fā)電機(jī)的次級(jí)運(yùn)動(dòng)頻率接近系統(tǒng)固有頻率即10 Hz時(shí)，次級(jí)運(yùn)動(dòng)的最大振幅為7.5 mm。

設(shè)次級(jí)運(yùn)動(dòng)位移與時(shí)間的關(guān)系近似采用正弦函數(shù)與反正切函數(shù)的乘積逼近，即有：

x(t)=A(t)sin2(2πft)

(2)

式中:x為位移；t為時(shí)間；A(t)=Aarctan(P(t))為振幅隨時(shí)間變化的反正切函數(shù)；P(t)為t的多項(xiàng)式；f為頻率。

對(duì)位移求導(dǎo)可得到次級(jí)運(yùn)動(dòng)速度與時(shí)間的關(guān)系：

v(t)=A′(t)sin2(2πft)+A(t)2πfsin(4πft)

(3)

由法拉第電磁感應(yīng)定律可知，感應(yīng)電動(dòng)勢(shì)表達(dá)式：

(4)

式中:ψ為磁鏈，且有：

(5)

式中:ψm為最大磁鏈；τp為極距。

將式(2)、式(3)、式(5)代入到式(4)中，可得感應(yīng)電動(dòng)勢(shì)表達(dá)式：

[A′(t)sin2(2πft)+A(t)2πfsin(4πft)]

(6)

由于本文以純電阻電路為研究對(duì)象，故其輸出電壓表達(dá)式：

U=E-Ir

(7)

式中:I為電流；r為內(nèi)阻。

當(dāng)P(t)選取為線性函數(shù)、二次函數(shù)和指數(shù)函數(shù)時(shí)，其在不同振幅倍數(shù)下的電壓值與有限元的仿真結(jié)果對(duì)比，如表2所示。

表2 不同變量函數(shù)下的電壓 (單位：V)

從表2可知，線性函數(shù)在1倍振幅和4倍振幅處、二次函數(shù)在0.5倍振幅和4倍振幅處、指數(shù)函數(shù)在0.5倍和4倍處的結(jié)果與仿真結(jié)果相接近，分析其均方誤差，線性函數(shù)、二次函數(shù)與指數(shù)函數(shù)的均方根差值分別為1.108、1.194、1.106。誤差較大的原因是初始數(shù)據(jù)的正弦度不夠好，并且相鄰振幅之間并不是逐漸上升，中間有少許畸變，綜合比較選取指數(shù)函數(shù)為振幅函數(shù)的變量函數(shù)。

圖15顯示了永磁直線發(fā)電機(jī)電壓、推力隨次級(jí)運(yùn)動(dòng)振幅倍數(shù)的變化曲線。振幅倍數(shù)對(duì)發(fā)電機(jī)輸出功率和效率的影響如圖16所示。

圖15 不同振幅下發(fā)電機(jī)輸出電壓和推力

圖16 不同振幅下發(fā)電機(jī)輸出功率和效率

從圖15可知，當(dāng)振幅增大時(shí)，發(fā)電機(jī)的電壓和電磁推力隨之增大，且電壓的增長(zhǎng)近似為直線，符合電動(dòng)勢(shì)的變化趨勢(shì)。

從圖16可知，發(fā)電機(jī)的輸出功率隨振幅倍數(shù)的增大而增大，且在2.5倍振幅處有最大值97.06 W；發(fā)電機(jī)的效率同樣隨著振幅倍數(shù)的增加而增加，但當(dāng)振幅倍數(shù)由1.0增大到2.5時(shí)，其效率將由85.44%增大到85.76%，即其增幅變緩，這是因?yàn)殡S著電壓的增加，發(fā)電機(jī)損耗增加得更快，其效率逐漸趨于某一定值。

4 結(jié) 論

本文基于無(wú)葉片風(fēng)力機(jī)捕能柱的渦激振動(dòng)特性，完成了永磁直線發(fā)電機(jī)的設(shè)計(jì)，在可行性得到驗(yàn)證的基礎(chǔ)上，探討了捕能柱的來流風(fēng)速和振幅對(duì)發(fā)電機(jī)電磁特性的影響，得到的結(jié)論如下：

(1)所設(shè)計(jì)的永磁直線發(fā)電機(jī)能夠與給定參數(shù)下的無(wú)葉片風(fēng)力機(jī)相匹配。

(2)當(dāng)來流風(fēng)速達(dá)到設(shè)計(jì)風(fēng)速即6 m/s時(shí)，發(fā)電機(jī)的瞬時(shí)功率達(dá)到最大即為35.62 W，其最大效率達(dá)到85.44%，且當(dāng)風(fēng)速處于鎖頻即共振區(qū)間時(shí)，發(fā)電機(jī)具有良好的輸出特性。

(3)通過分析確定了指數(shù)函數(shù)作為振幅函數(shù)的變量函數(shù)，且隨著捕能柱擺動(dòng)幅值的增大，發(fā)電機(jī)的輸出功率和效率均增大，但當(dāng)振幅倍數(shù)大于1時(shí)，發(fā)電機(jī)的效率增幅變緩。

上述所得結(jié)論可為應(yīng)用于無(wú)葉片風(fēng)力機(jī)的永磁直線發(fā)電機(jī)的設(shè)計(jì)提供指導(dǎo)。