張 欣， 宋子晗， 王文斌， 韓 瑜， 祝麗花

(天津工業(yè)大學(xué) 天津市電工電能新技術(shù)重點(diǎn)實(shí)驗(yàn)室，天津 300387)

電機(jī)被廣泛應(yīng)用于人類(lèi)生活中，尤其應(yīng)用于艦船驅(qū)動(dòng)等軍工領(lǐng)域以及重工業(yè)行業(yè)等對(duì)大型電機(jī)的低噪聲要求更加嚴(yán)格，永磁同步電機(jī)的平穩(wěn)和安靜運(yùn)行性能在高性能應(yīng)用中非常重要。

電機(jī)的振動(dòng)噪聲研究起源于20世紀(jì)30年代，而磁致伸縮效應(yīng)是在1842年由英國(guó)學(xué)者發(fā)現(xiàn)的。在2004年美國(guó)學(xué)者M(jìn)ohammed等[1]測(cè)量了逆向磁致伸縮對(duì)電工鋼的磁化特性的影響，并用虛功原理計(jì)算出了磁致伸縮力。Zhu等[2]建立了包括磁致伸縮的磁彈耦合模型，用于疊片鐵心的振動(dòng)分析，然后應(yīng)用數(shù)值模型并測(cè)量磁化強(qiáng)度和MS曲線，分析發(fā)現(xiàn)定子中的變形不僅由磁阻力引起，還由磁致伸縮效應(yīng)和諧波引起，而且磁致伸縮和諧波非常重要，這對(duì)于設(shè)計(jì)更低振動(dòng)和噪聲的電動(dòng)機(jī)很有意義。Ghalamestani等[3]同時(shí)考慮電磁力和磁致伸縮力對(duì)電機(jī)定子的影響，用等效力表示出了電機(jī)定子的形變，分析發(fā)現(xiàn)電磁力和磁致伸縮力在一定情況下會(huì)互相增強(qiáng)和抵消。

韓雪巖對(duì)電工硅鋼片的磁致伸縮特性和不同電流下的電機(jī)振動(dòng)噪聲進(jìn)行了測(cè)量和分析，發(fā)現(xiàn)當(dāng)考慮磁致伸縮效應(yīng)時(shí)，測(cè)試結(jié)果更接近實(shí)際值[4-5]。祝麗花[6]測(cè)量了硅鋼片在有無(wú)絕緣層時(shí)磁致伸縮效應(yīng)對(duì)應(yīng)力的影響以及變壓器鐵心硅鋼片鐵軛和中柱在軋制方向與垂直方向的應(yīng)力與應(yīng)變的關(guān)系。張欣等[7]通過(guò)對(duì)不同類(lèi)型的電機(jī)進(jìn)行定子振動(dòng)響應(yīng)測(cè)試，發(fā)現(xiàn)電機(jī)定子硅鋼片的主要振動(dòng)頻率為供電頻率的一倍和二倍電源頻率處，而且以二倍供電頻率為主。吳勝男等[8]通過(guò)建立電機(jī)定子解析模型對(duì)定子振動(dòng)與磁致伸縮的關(guān)系進(jìn)行了解析計(jì)算，發(fā)現(xiàn)由磁致伸縮效應(yīng)引起的定子上的振動(dòng)與硅鋼片的磁致伸縮系數(shù)成正比。閆榮格等[9-10]測(cè)試了不同含量諧波情況下的電機(jī)振動(dòng)加速度與電機(jī)的應(yīng)力分布情況，分析了不同諧波對(duì)電機(jī)振動(dòng)噪聲的影響。馮大軍等[11]對(duì)不同種類(lèi)的硅鋼片在不同退火工藝下進(jìn)行了磁致伸縮性能測(cè)試，發(fā)現(xiàn)磁致伸縮與退火時(shí)間和磁化方向等因素有關(guān)。Kuroishi等[12]通過(guò)進(jìn)行靜磁耦合和結(jié)構(gòu)有限元分析來(lái)估計(jì)定子的振動(dòng)形狀，并與在每個(gè)邊帶頻率上計(jì)算出的振動(dòng)形狀進(jìn)行了對(duì)比試驗(yàn)。Yamagashira等[13]在旋轉(zhuǎn)通量條件下，測(cè)量了非定向電工鋼板二維磁致伸縮特性，明確了磁通密度向量，磁場(chǎng)強(qiáng)度向量和二維磁致伸縮之間的關(guān)系。

Belahcen[14]研究發(fā)現(xiàn)在考慮耦合時(shí)，同步發(fā)電機(jī)定子的計(jì)算振動(dòng)是不同的。取決于所考慮的振動(dòng)頻率，由于磁彈性耦合，觀察到振動(dòng)幅度的增大或減小大約20%。Zhang等[15]測(cè)量并分析了交變磁化加直流偏置下硅鋼片的磁致伸縮各向異性和應(yīng)力依賴性，提出了一種利用直流偏置下的磁致伸縮曲線的磁性有限元分析方法。孫小光[16]通過(guò)Matlab繪制了電機(jī)定子中磁致伸縮力與外加磁場(chǎng)方向和大小的插值曲面圖，發(fā)現(xiàn)磁致伸縮力的方向與磁場(chǎng)方向相同且大小相關(guān)。王園弟等[17]對(duì)旋轉(zhuǎn)磁化下的無(wú)取向硅鋼片進(jìn)行了應(yīng)力分析，發(fā)現(xiàn)磁致伸縮效應(yīng)下的鐵軛應(yīng)變主要是平行軋制方向，而定子齒部最大應(yīng)變則垂直與鐵軛應(yīng)變的方向。郝清亮等[18]通過(guò)試驗(yàn)對(duì)比了磁致伸縮力和麥克斯韋力對(duì)電機(jī)振動(dòng)噪聲的貢獻(xiàn)分析，發(fā)現(xiàn)雖然磁致伸縮力的貢獻(xiàn)小于麥克斯韋力，但對(duì)磁致伸縮力的研究也是非常重要的。馮旭[19]對(duì)無(wú)取向硅鋼片在交變磁場(chǎng)中的磁致伸縮效應(yīng)進(jìn)行分析，又對(duì)永磁同步電機(jī)的振動(dòng)進(jìn)行了試驗(yàn)測(cè)試，發(fā)現(xiàn)當(dāng)考慮磁致伸縮效應(yīng)時(shí)電機(jī)的振動(dòng)噪聲大于不考慮磁致伸縮效應(yīng)時(shí)的振動(dòng)噪聲。Xiao等[20]提出了一種解耦方法來(lái)處理鐵磁材料中的磁致伸縮效應(yīng)和反磁致伸縮效應(yīng)的計(jì)算。將從磁場(chǎng)溶液中導(dǎo)出的磁致伸縮轉(zhuǎn)換為體積力密度，并將其用作外部載荷，通過(guò)沿每個(gè)局部主軸重構(gòu)非線性BH曲線，可以考慮反磁致伸縮對(duì)磁場(chǎng)的影響。趙小軍等[21]中構(gòu)建了硅鋼片的直流偏磁動(dòng)態(tài)磁致模型，用于模擬硅鋼片在直流偏磁下的損耗特性。翁玲等[22]測(cè)量了多種磁致伸縮材料在不同磁感應(yīng)強(qiáng)度和勵(lì)磁頻率下的磁導(dǎo)率以及電磁損耗情況。

電磁力導(dǎo)致的振動(dòng)噪聲和磁致伸縮力導(dǎo)致的振動(dòng)噪聲可以分別采取一定方法來(lái)降低，最終降低電機(jī)整體振動(dòng)噪聲。目前電磁力部分的噪聲問(wèn)題已經(jīng)有很多成熟的降噪方案，而由磁致伸縮力導(dǎo)致的振動(dòng)噪聲的抑制方法還鮮有研究。本文提出通過(guò)填充具有負(fù)磁致伸縮特性的軟磁復(fù)合材料來(lái)抑制電機(jī)振動(dòng)噪聲的方法，使擬填充材料的負(fù)磁致伸縮與硅鋼片的正磁致伸縮效應(yīng)引起的形變基本相互抵消，從而減小電機(jī)的振動(dòng)噪聲。首先建立電機(jī)定子的電磁-機(jī)械耦合數(shù)值模型進(jìn)行有限元計(jì)算仿真，對(duì)電機(jī)定子填充負(fù)磁致伸縮材料，同時(shí)比較了不同填充位置和距離對(duì)定子面應(yīng)力的影響情況。然后根據(jù)有限元仿真結(jié)果搭建了電機(jī)定子應(yīng)變測(cè)量試驗(yàn)平臺(tái)，通過(guò)對(duì)實(shí)際電機(jī)定子打孔并填充負(fù)磁致伸縮材料鎳進(jìn)行試驗(yàn)，驗(yàn)證了有限元計(jì)算的結(jié)果。

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

首先建立考慮磁致伸縮效應(yīng)的電磁-機(jī)械的耦合模型，為了測(cè)試負(fù)磁致伸縮材料對(duì)電機(jī)定子振動(dòng)噪聲的影響，本文選取有限元仿真軟件COMSOL Multiphysics建立了二維模型，對(duì)電機(jī)定子進(jìn)行了仿真試驗(yàn)。其磁-機(jī)械耦合數(shù)值模型可用以下公式表示

[S][A]=[Je]

(1)

[K][U]=[F]

(2)

式中：K表示機(jī)械剛度矩陣；S表示電磁剛度矩陣；Je是外部電流密度雅可比矩陣；F是作用力；通過(guò)求解得到振動(dòng)位移U；矢量磁勢(shì)A。

固體力學(xué)中線彈性材料的研究方程

(3)

F=I+?u

(4)

式中：u表示位移矢量；F是載荷量；ρ表示密度；S是應(yīng)變；υ代表泊松比。

對(duì)于彈性張量D，在三維線彈性材料中可表示為

(5)

在有限元中固體力學(xué)模塊的求解方程如下所示

(6)

T-Ti=C(S-Si)

(7)

式中：ρ是密度；μ是位移矢量；F是載荷量；T是所受應(yīng)力；Ti是所給出的應(yīng)力的初值；C是剛度矩陣；Si是初始應(yīng)變?cè)O(shè)定值。

對(duì)于電機(jī)定子硅鋼片來(lái)說(shuō)，其磁致伸縮面的應(yīng)力求解公式為

σ=Dξ

(8)

式中：σ為磁致伸縮應(yīng)力；D為彈性張量；ξ為磁致伸縮應(yīng)變。

2 仿真結(jié)果與分析

本文采用金屬鎳作為填充材料，首先鎳容易成型，比較容易獲得，而且有著很好的機(jī)械強(qiáng)度和延展性，同時(shí)空氣中不易氧化，具備負(fù)磁致伸縮性能，所以本次試驗(yàn)選用鎳(Ni)進(jìn)行了試驗(yàn)。Ni的參數(shù)如表1所示。而后基于以上的數(shù)學(xué)模型對(duì)永磁同步電機(jī)的定子模型進(jìn)行了限元仿真，包括打孔填充負(fù)磁致伸縮材料前后電機(jī)定子的面應(yīng)力變化情況以及點(diǎn)應(yīng)力變化情況對(duì)比。

表1 Ni的屬性值

2.1 初始條件

在有限元仿真軟件COMSOLMultiphysics構(gòu)建電機(jī)定子模型，為了減少轉(zhuǎn)子磁場(chǎng)對(duì)定子的影響，本試驗(yàn)去掉了電機(jī)轉(zhuǎn)子，這樣可以減小外部氣隙中的漏磁，從而使磁通主要集中在定子鐵心內(nèi)部，減小電機(jī)定子齒表面的電磁力，使電機(jī)定子鐵心的主要應(yīng)力為磁致伸縮力。

2.2 電機(jī)定子面應(yīng)力仿真結(jié)果

圖1為電機(jī)定子填充負(fù)磁致伸縮材料前的局部面應(yīng)力分布情況，圖2為打孔填充負(fù)磁致伸縮材料后的局部面應(yīng)力分布。兩圖對(duì)比發(fā)現(xiàn)，在定子填充負(fù)磁致伸縮材料后最大面應(yīng)力由1.69×107N/m2減小到1.67×107N/m2，變化量很少，但是可以明顯地看出圖2中左上方深色區(qū)域(應(yīng)力較大區(qū)域)明顯減少了。為了更直接的看出負(fù)磁致伸縮材料填充前后的應(yīng)力變化情況，于是對(duì)圖3中三個(gè)位置A、B、C處的點(diǎn)應(yīng)力分別進(jìn)行了測(cè)量。

圖1 未打孔時(shí)局部應(yīng)力

圖2 打孔填充后局部應(yīng)力

2.3 電機(jī)定子點(diǎn)應(yīng)力仿真結(jié)果

如圖3所示，選取的三個(gè)測(cè)量點(diǎn)為A、B、C，點(diǎn)A位于打孔附近區(qū)域，點(diǎn)B和點(diǎn)C距離打孔位置較遠(yuǎn)，由圖4和圖5可以看出在打孔添加負(fù)磁致伸縮材料后測(cè)量點(diǎn)的應(yīng)力變化情況。

圖3 測(cè)量點(diǎn)

在本次仿真中孔3為填充負(fù)磁致伸縮材料的一個(gè)較優(yōu)點(diǎn)，孔1和孔2為隨機(jī)兩個(gè)位置。

較優(yōu)點(diǎn)尋找方法為：首先在COMSOL中設(shè)定打孔的形狀為圓形，打孔半徑為r，孔圓心到定子模型中心距離為R，孔圓心與定子模型中心連線與水平方向夾角為angle，如圖4所示。

圖4 測(cè)量點(diǎn)

然后在COMSOL中進(jìn)行參數(shù)化掃描，對(duì)三項(xiàng)數(shù)據(jù)設(shè)定好步長(zhǎng)以及掃描范圍，其中r范圍設(shè)定為1～3 mm，R范圍設(shè)定為55～90 mm，angle范圍設(shè)定為0°～90°。對(duì)多組填充方案下的電機(jī)應(yīng)力進(jìn)行有限元參數(shù)化掃描計(jì)算，然后對(duì)每組參數(shù)下的應(yīng)力進(jìn)行比較，選出最小應(yīng)力下的填充方案作為最較優(yōu)填充方案。

對(duì)比圖5和圖6(c)，當(dāng)在孔3處填充負(fù)磁致伸縮材料后，點(diǎn)A應(yīng)力由0.95×106N/m2減小為0.38×106N/m2，減小了60%；點(diǎn)B由1.50×106N/m2變?yōu)?.52×106N/m2，點(diǎn)C處應(yīng)力由1.25×106N/m2改變?yōu)?.26×106N/m2，點(diǎn)B和點(diǎn)C應(yīng)力變化很小。

圖5 測(cè)量點(diǎn)原應(yīng)力

(a) 孔1添加Ni

對(duì)比圖5和圖6(a)，當(dāng)在孔1填充負(fù)磁致伸縮材料后，點(diǎn)A應(yīng)力由0.95×106N/m2變?yōu)?.8×106N/m2，與原應(yīng)力相比減小了15.7%；點(diǎn)B應(yīng)力由1.50×106N/m2變?yōu)?.52×106N/m2，點(diǎn)C應(yīng)力由1.25×106N/m2變?yōu)?.29×106N/m2，點(diǎn)B和C應(yīng)力變化很小。

對(duì)比圖5和圖6(b)當(dāng)在孔2填充負(fù)磁致伸縮材料后，點(diǎn)A應(yīng)力變?yōu)?.20×106N/m2，與原應(yīng)力相比增大了26%；點(diǎn)B應(yīng)力為1.55×106N/m2，點(diǎn)C變?yōu)?1.29×106N/m2，點(diǎn)B和C應(yīng)力在孔2填充負(fù)磁致伸縮材料前后變化不大。

由以上三組對(duì)比可以發(fā)現(xiàn)，打孔填充負(fù)磁致伸縮材料后，對(duì)填充點(diǎn)附近區(qū)域的影響較大，而隨著距離變大后，對(duì)較遠(yuǎn)距離的影響較??；并且只有在合適的位置打孔并填充負(fù)磁致伸縮材料后才會(huì)使應(yīng)力減小，隨意打孔填充負(fù)磁致伸縮材料后可能會(huì)使應(yīng)力增大，產(chǎn)生負(fù)面影響。

3 試驗(yàn)結(jié)果與分析

根據(jù)仿真結(jié)果，制定了試驗(yàn)方案，分別測(cè)量了不填充材料和打孔填充負(fù)磁致伸縮材料后的電機(jī)定子面應(yīng)力，并選出幾個(gè)點(diǎn)對(duì)比其填充材料前后的應(yīng)力變化。

3.1 試驗(yàn)儀器與準(zhǔn)備

如圖7所示，本試驗(yàn)采用的為三片直角型應(yīng)變片，可以對(duì)測(cè)量位置的0°、45°、90°三個(gè)方向進(jìn)行應(yīng)變測(cè)量。

圖7 應(yīng)變片

本文根據(jù)仿真思路設(shè)計(jì)了試驗(yàn)裝置，如圖8所示，去掉一個(gè)永磁同步電機(jī)的轉(zhuǎn)子，這樣磁通主要集中約束在定子鐵心內(nèi)部，外部氣隙中的漏磁通很小，定子齒表面的電磁力很小，定子鐵心的主要應(yīng)力為磁致伸縮力。

圖8 電機(jī)定子

在電機(jī)4個(gè)位置粘貼應(yīng)變片，并分別命名為A、B、C、D，應(yīng)變片A和B周?chē)鷽](méi)有孔，C右側(cè)為直徑4 mm孔，D左右兩側(cè)都有4 mm孔。從左到右分別為編號(hào)1、2、3、4、5、6的直徑4 mm孔。其中孔1為添加負(fù)磁致伸縮材料的適合位置之一。

圖9中利用變頻器給電機(jī)繞組中通入50 Hz的工頻電產(chǎn)生交變旋轉(zhuǎn)磁場(chǎng)，利用JHDY 動(dòng)態(tài)應(yīng)變測(cè)量系統(tǒng)測(cè)量電機(jī)定子上應(yīng)變片粘貼位置的應(yīng)變，采樣頻率設(shè)定為1 000 Hz，在電機(jī)定子線圈未通電時(shí)對(duì)數(shù)據(jù)進(jìn)行調(diào)零，消除環(huán)境的影響，再通電進(jìn)行測(cè)量。每次改變硅鋼片片數(shù)后通電測(cè)量前都要重新調(diào)零。之后將數(shù)據(jù)在電腦中進(jìn)行保存。

圖9 電機(jī)應(yīng)變測(cè)量裝置

3.2 孔1添加鎳后對(duì)應(yīng)變片B和C的影響

圖10中兩圖為應(yīng)變片B在孔1打孔填充負(fù)磁致伸縮材料前后的應(yīng)變對(duì)比?？?添加負(fù)磁致伸縮材料鎳后應(yīng)變峰值由165.913 με變?yōu)榱?63.217 με，變化量為1.6%，谷值由-111.283 με變?yōu)?108.878 με，變化量為2.1%。由于孔1與應(yīng)變片B距離較遠(yuǎn)，大約70°，所以變化量很小，與仿真規(guī)律相同。

(a) 添加材料前B應(yīng)變

圖11(a)、(b)、(c)三圖分別為應(yīng)變片C在未添加負(fù)磁致伸縮材料時(shí)、僅孔1填加Ni時(shí)、孔1和4添加Ni時(shí)的應(yīng)變對(duì)比?？?與應(yīng)變片C相差約43.6°、孔5與應(yīng)變片C相差約60°，兩孔都距離應(yīng)變片C較遠(yuǎn)；而孔1緊挨著應(yīng)變片C，距離非常近。

圖11(a)和圖11(b)對(duì)比，發(fā)現(xiàn)當(dāng)在孔1填充負(fù)磁致伸縮材料后，應(yīng)變片C峰值由148.424 με減小為124.110 με，應(yīng)變峰值減小了16%；圖11(b)和圖11(c)對(duì)比，應(yīng)變峰值由124.110 με變?yōu)?24.054 με，峰值變化量為0.05%；谷值由-96.699 με變?yōu)?90.659 με，變化量為6.2%，總體考慮誤差后變化較小。已知孔4距離應(yīng)變片C較遠(yuǎn)，而孔1距離應(yīng)變片C較近，所以孔1填充負(fù)磁致伸縮材料后對(duì)應(yīng)變片C的影響較大，而孔4填充負(fù)磁致伸縮材料則影響較小。

(a) 添加材料前C應(yīng)變

3.3 孔1添加鎳后應(yīng)變片B、C、D應(yīng)變對(duì)比

圖12為應(yīng)變片D在孔1添加Ni后的應(yīng)變，結(jié)合圖10(b)孔1添加Ni后B應(yīng)變圖，和圖11(b)孔1添加Ni后C應(yīng)變圖。發(fā)現(xiàn)應(yīng)變片C的應(yīng)變峰值為124.110 με，遠(yuǎn)小于應(yīng)變片B的163.217 με和D的179.540 με?？芍獞?yīng)變片C距離孔1較近，而且孔1未添加Ni時(shí)與僅孔1添加Ni時(shí)應(yīng)變片C最大應(yīng)變減小了16%。

圖12 孔1添加Ni后D應(yīng)變

以上對(duì)比可發(fā)現(xiàn)當(dāng)在電機(jī)定子合適的位置填充負(fù)磁致伸縮材料后，可以使填充材料附近的應(yīng)變減小。磁致伸縮效應(yīng)會(huì)使電機(jī)定子伴隨電源頻率產(chǎn)生規(guī)律性的伸長(zhǎng)和縮短，從而引起電機(jī)本身的振動(dòng)和內(nèi)部應(yīng)力，同時(shí)硅鋼片表面的伸長(zhǎng)和縮短正是應(yīng)變的來(lái)源，所以說(shuō)應(yīng)變量減小了，對(duì)應(yīng)的內(nèi)部應(yīng)力也會(huì)減小，從而電機(jī)的振動(dòng)噪聲也就減小了。

還可發(fā)現(xiàn)打孔填充負(fù)磁致伸縮材料后對(duì)填充材料位置附近區(qū)域的影響較大，而距離填充負(fù)磁致伸縮材料位置較遠(yuǎn)的區(qū)域，則影響較小。這是因?yàn)樗钊胴?fù)磁致伸縮材料產(chǎn)生的形變改變了電機(jī)定子內(nèi)部的應(yīng)力，而由形變產(chǎn)生的內(nèi)部應(yīng)力又是隨著距離增大而遞減的，所以隨著測(cè)量點(diǎn)與填充位置處距離的增大，測(cè)量點(diǎn)的應(yīng)變所受到的影響也就越小。

4 結(jié) 論

本文首先建立了永磁同步電機(jī)的磁-機(jī)械耦合數(shù)值模型，在有限元中對(duì)電機(jī)定子應(yīng)力進(jìn)行了系統(tǒng)的分析。并根據(jù)仿真模型制作真實(shí)的電機(jī)模型，進(jìn)行電機(jī)定子打孔填充負(fù)磁致伸縮材料的仿真試驗(yàn)進(jìn)行驗(yàn)證，試驗(yàn)結(jié)果與仿真基本相符。得出結(jié)論如下：

(1) 根據(jù)擬填充負(fù)磁致伸縮材料的特性找到合適的打孔填充位置和尺寸，然后在電機(jī)定子上打孔并填充負(fù)磁致伸縮材料，可以使電機(jī)定子的應(yīng)變和面應(yīng)力減小，從而達(dá)到減振降噪的目的。

(2) 在電機(jī)定子上打孔填充負(fù)磁致伸縮材料后，會(huì)使定子的面應(yīng)力發(fā)生變化，而應(yīng)力改變的幅度與測(cè)量點(diǎn)到填充位置處的距離成反比。距離越近，應(yīng)力變化的幅度越大；距離越遠(yuǎn)，應(yīng)力變化的幅度越小。本試驗(yàn)驗(yàn)證了局部應(yīng)力的變化情況，在電機(jī)其他定子齒對(duì)稱位置同樣打孔填加材料可以減小電機(jī)定子整體應(yīng)力。

此結(jié)論對(duì)低噪聲永磁同步電機(jī)的設(shè)計(jì)有重要的參考價(jià)值。