姜 濤，曲詩(shī)健，孫佳偉

(中車(chē)大連電力牽引研發(fā)中心有限公司，遼寧 大連 116052)

0 引 言

近年來(lái)，隨著稀土永磁材料的發(fā)展，永磁電機(jī)應(yīng)用越來(lái)越廣泛。永磁電機(jī)具有高轉(zhuǎn)矩密度、高效率、高功率因數(shù)等優(yōu)點(diǎn)[1]，滿足了現(xiàn)代高性能伺服系統(tǒng)對(duì)寬調(diào)速范圍、低速大轉(zhuǎn)矩、穩(wěn)定性好與快速動(dòng)態(tài)響應(yīng)的要求。但齒槽轉(zhuǎn)矩是導(dǎo)致永磁電機(jī)轉(zhuǎn)矩脈動(dòng)的主要原因之一[2]，影響伺服電機(jī)控制精度，所以永磁伺服電機(jī)齒槽轉(zhuǎn)矩的抑制尤為重要。

齒槽轉(zhuǎn)矩是永磁體與定子齒之間相互吸引產(chǎn)生的[3]?，F(xiàn)有文獻(xiàn)提出了各種各樣齒槽轉(zhuǎn)矩抑制方法，但絕大多數(shù)是以理想模型為基礎(chǔ)的理論分析與仿真計(jì)算，并未計(jì)及批量生產(chǎn)中可能存在的加工誤差。包括定子參數(shù)優(yōu)化的定子斜槽、不等槽口寬、定子齒開(kāi)虛槽等方法，轉(zhuǎn)子參數(shù)優(yōu)化的優(yōu)化極靴深度、轉(zhuǎn)子斜極、優(yōu)化極弧系數(shù)、磁極分段移位、磁極采用不對(duì)稱(chēng)角度、采用不等極弧系數(shù)等方法[4-11]，都是基于理想模型的。

計(jì)及加工誤差后，齒槽轉(zhuǎn)矩為兩部分的疊加，第一部分為理想模型的齒槽轉(zhuǎn)矩TNHC(Native Harmonic Components, NHC),另一部分為加工誤差產(chǎn)生的額外齒槽轉(zhuǎn)矩TAHC(Additional Harmonic Components, AHC)[12]，TAHC的產(chǎn)生導(dǎo)致文獻(xiàn)[4-11]中的齒槽轉(zhuǎn)矩抑制方法理論與實(shí)際區(qū)別較大。文獻(xiàn)[13]定義了一個(gè)與槽數(shù)和極數(shù)有關(guān)的參數(shù)F，用來(lái)衡量TAHC對(duì)加工誤差的敏感度。文獻(xiàn)[14-15]研究了相同極性永磁體沿同一方向位移誤差、單個(gè)定子齒誤差和單個(gè)槽口寬誤差對(duì)齒槽轉(zhuǎn)矩的影響規(guī)律，但是沒(méi)有考慮加工誤差的隨機(jī)性。文獻(xiàn)[16]為了加快有限元的計(jì)算速度，將N極和S極永磁體(Permanent Magnet，PM)分別取相鄰的一半作為一個(gè)單元，并對(duì)只考慮這個(gè)單元永磁體的電機(jī)進(jìn)行有限元仿真，將齒槽轉(zhuǎn)矩計(jì)算結(jié)果進(jìn)行位移疊加得到最終的齒槽轉(zhuǎn)矩，但只能針對(duì)理想模型。文獻(xiàn)[17]從泊松方程和拉普拉斯方程出發(fā)，推導(dǎo)考慮加工誤差的齒槽轉(zhuǎn)矩解析模型，模型計(jì)算結(jié)果與有限元算法接近，但計(jì)算推導(dǎo)過(guò)程較為復(fù)雜。然而，以上文獻(xiàn)均未分析齒槽轉(zhuǎn)矩對(duì)各加工誤差的敏感性以及加工誤差對(duì)齒槽轉(zhuǎn)矩概率分布的影響。

本文首先采用能量法建立解析模型，分別計(jì)算單一加工誤差對(duì)齒槽轉(zhuǎn)矩頻譜的影響規(guī)律，以此判斷對(duì)齒槽轉(zhuǎn)矩影響最大的誤差種類(lèi)，指出了齒槽轉(zhuǎn)矩實(shí)際與理論差距大的本質(zhì)原因。再基于蒙特卡羅方法，結(jié)合有限元仿真軟件，計(jì)算同時(shí)考慮所有加工誤差后的齒槽轉(zhuǎn)矩概率密度分布，以此預(yù)測(cè)量產(chǎn)電機(jī)齒槽轉(zhuǎn)矩合格率。并且采用定子沖片循環(huán)旋轉(zhuǎn)方法較好地抑制了定子側(cè)誤差產(chǎn)生的TAHC，同時(shí)進(jìn)行了理論論證，得到了沖片依次旋轉(zhuǎn)一個(gè)槽時(shí)齒槽轉(zhuǎn)矩抑制效果最好的結(jié)論。最后，用實(shí)驗(yàn)驗(yàn)證了本文齒槽轉(zhuǎn)矩分析抑制方法的可行性。

1 能量法模型

1.1 理想模型

根據(jù)能量法，齒槽轉(zhuǎn)矩等于電樞繞組不通電時(shí)，永磁電機(jī)內(nèi)部能量對(duì)定轉(zhuǎn)子相對(duì)位置角的負(fù)導(dǎo)數(shù)。忽略定轉(zhuǎn)子鐵心中的能量，只考慮氣隙和永磁體中的能量[18]，齒槽轉(zhuǎn)矩為

(1)

(2)

如圖1所示，假定θ=0位置為指定的兩個(gè)永磁體中心線，α為某一指定的定子槽中心線和θ=0位置所夾角度，即定轉(zhuǎn)子相對(duì)位置角。

圖1 永磁體與電樞相對(duì)位置

圖2 的分布

圖3 定子槽中磁力線路徑

永磁體對(duì)著定子槽的位置，假定槽中磁力線以四分之一圓弧路徑進(jìn)入定子齒，圓弧半徑Rs等于磁力線進(jìn)入位置到定子齒頂點(diǎn)的距離，氣隙長(zhǎng)度δ等于物理氣隙δ0加上以Rs為半徑的四分之一圓弧長(zhǎng)度[19]。電機(jī)齒槽分布均勻，G2(θ,α)在周向上的分布如圖4所示，So為槽口寬對(duì)應(yīng)的機(jī)械角度。將式(4)帶入式(2)，可得G2(θ,α)傅里葉展開(kāi)如式:

圖4 G2(q,α)的分布

(3)

(4)

(5)

將式(3)、式(5)帶入式(1)可得：

(6)

式中，La為電機(jī)鐵心軸向長(zhǎng)度，R1和R2分別為轉(zhuǎn)子外半徑和定子內(nèi)半徑，z為槽數(shù)，p為極對(duì)數(shù)，q為使得qz/(2p)為整數(shù)的整數(shù)

1.2 非理想模型

電機(jī)批量生產(chǎn)中，即使微小的誤差，也可能造成齒槽轉(zhuǎn)矩成倍增長(zhǎng)。后文采用定子沖片循環(huán)旋轉(zhuǎn)方法，能較好地抑制定子側(cè)誤差對(duì)TAHC的影響，故此處僅分析轉(zhuǎn)子側(cè)永磁體充磁、厚度、寬度、周向位置等誤差對(duì)齒槽轉(zhuǎn)矩的影響。

(7)

(8)

(9)

(10)

式中,m=n=q′z，q′為正整數(shù)。相比于式(6)新增了槽數(shù)倍數(shù)次齒槽轉(zhuǎn)矩諧波。

1.3 加工誤差敏感性分析

假定每塊永磁體都各不相同，對(duì)于12槽10極電機(jī)來(lái)說(shuō)，每塊永磁體都將產(chǎn)生12次齒槽轉(zhuǎn)矩諧波?；?2次諧波和相鄰永磁體機(jī)械角度(36°)可以得到永磁體誤差導(dǎo)致的額外齒槽轉(zhuǎn)矩矢量圖[19]，如圖5所示。相鄰兩塊永磁體產(chǎn)生的12次齒槽轉(zhuǎn)矩諧波相位依次相差72°。根據(jù)該矢量圖與表5所列加工誤差的分布，可以直接得到各類(lèi)誤差對(duì)齒槽轉(zhuǎn)矩影響最?lèi)毫拥那闆r并列于表1中。

圖5 齒槽轉(zhuǎn)矩矢量圖

表1 計(jì)算模型永磁體參數(shù)

分別考慮永磁體4種誤差，并結(jié)合表2中的電機(jī)標(biāo)稱(chēng)參數(shù)，代入能量法模型，算得各次諧波幅值如圖6所示。永磁體各種加工誤差都將引入槽數(shù)倍數(shù)次的齒槽轉(zhuǎn)矩諧波，并且對(duì)齒槽轉(zhuǎn)矩影響最大的誤差為永磁體周向位置誤差。即永磁體周向位置是導(dǎo)致齒槽轉(zhuǎn)矩理論與實(shí)際相差大的首要原因。因此，加強(qiáng)對(duì)永磁體周向位置加工精度，將產(chǎn)生較好的齒槽轉(zhuǎn)矩抑制效果。

圖6 各次齒槽轉(zhuǎn)矩諧波幅值

表2 電機(jī)標(biāo)稱(chēng)參數(shù)

2 有限元模型

2.1 理想模型

本文針對(duì)一個(gè)12槽10極的表貼式永磁伺服電機(jī)進(jìn)行研究分析，電機(jī)標(biāo)稱(chēng)參數(shù)列在表2中，拓?fù)浣Y(jié)構(gòu)如圖7所示。有限元計(jì)算的理想模型齒槽轉(zhuǎn)矩波形及頻譜如圖8所示。在轉(zhuǎn)子轉(zhuǎn)動(dòng)的360°機(jī)械角度內(nèi)，齒槽轉(zhuǎn)矩周期性變化，變化周期數(shù)為槽數(shù)與極數(shù)的最小公倍數(shù)[15]。

圖7 電機(jī)拓?fù)浣Y(jié)構(gòu)

圖8 理想模型齒槽轉(zhuǎn)矩

2.2 齒槽轉(zhuǎn)矩概率分布

如果隨機(jī)變量X服從一個(gè)以μ為數(shù)學(xué)期望、σ2為方差的正態(tài)分布，即N(μ，σ2)。那么X落在區(qū)間(μ-3σ，μ+3σ)以外的概率小于千分之三，基本上可以把區(qū)間(μ-3σ，μ+3σ)看作是隨機(jī)變量X實(shí)際可能的取值區(qū)間，這稱(chēng)為正態(tài)分布的“3σ”原則。

假定實(shí)際加工中電機(jī)各參數(shù)都服從正態(tài)分布，數(shù)學(xué)期望為對(duì)應(yīng)參數(shù)標(biāo)稱(chēng)值，方差通過(guò)對(duì)樣機(jī)各參數(shù)測(cè)量得到。對(duì)一臺(tái)樣機(jī)的定子齒直徑、PM寬度、PM厚度、PM表面磁密等參數(shù)進(jìn)行測(cè)量，測(cè)量結(jié)果分別如表3、表4所示，其中，定子齒直徑定義為兩個(gè)相對(duì)定子齒的距離。

表3 定子齒直徑誤差

表4 永磁體誤差

根據(jù)正態(tài)分布的“3σ”原則，各參數(shù)的標(biāo)準(zhǔn)差σ應(yīng)等于相應(yīng)最大誤差的三分之一。但考慮到本次測(cè)量只是針對(duì)一個(gè)樣機(jī)，樣本數(shù)量有限，各參數(shù)的標(biāo)準(zhǔn)差取得稍大。定子齒直徑最大誤差為0.04 mm，建立有限元模型時(shí)，參數(shù)化定子齒半徑，取其標(biāo)準(zhǔn)差為0.025/3 mm；PM寬度最大誤差為0.05°，取標(biāo)準(zhǔn)差0.060°/3；PM厚度的最大誤差為0.03 mm，取標(biāo)準(zhǔn)差0.040/3 mm；PM表面磁密最大誤差為1.96%，取剩磁標(biāo)準(zhǔn)差為0.025/3T。

另外，有限元軟件建模時(shí)，參數(shù)化定子槽口寬的一半，如圖9所示，這樣能同時(shí)考慮槽口寬度與位置誤差。且PM周向位置與槽口寬誤差由供應(yīng)商提供，取PM周向位置標(biāo)準(zhǔn)差為0.600°/3，取半個(gè)槽口寬標(biāo)準(zhǔn)差為0.009/3 mm。具體參數(shù)列在表5中。

圖9 槽口寬的一半

表5 加工誤差分布

根據(jù)各類(lèi)誤差的概率分布，隨機(jī)生成300組考慮加工誤差的參數(shù)樣本，如圖10所示，并將各組誤差樣本進(jìn)行隨機(jī)組合，建立相應(yīng)的不對(duì)稱(chēng)有限元模型，計(jì)算各模型的齒槽轉(zhuǎn)矩。

圖10 誤差樣本

對(duì)有限元仿真結(jié)果進(jìn)行統(tǒng)計(jì)分析，齒槽轉(zhuǎn)矩幅值與諧波分布分別如圖11、圖12所示?？紤]加工誤差后，產(chǎn)生了TAHC，使得所有樣本齒槽轉(zhuǎn)矩幅值均大于理想模型的19 mNm，且齒槽轉(zhuǎn)矩幅值分布范圍很廣，從25 mNm到58 mNm。TAHC諧波主要包括定子側(cè)誤差產(chǎn)生的10次諧波和轉(zhuǎn)子側(cè)誤差產(chǎn)生的12次、24次諧波。其中12次諧波的分布范圍最廣，從0.8 mNm到35.5 mNm，是導(dǎo)致齒槽轉(zhuǎn)矩幅值分布范圍廣的主要原因。并且根據(jù)能量法模型中的分析，齒槽轉(zhuǎn)矩12次諧波分布范圍廣是由PM周向位置誤差導(dǎo)致的，再次印證了提高永磁體周向位置加工精度是抑制齒槽轉(zhuǎn)矩的有效措施。60次齒槽轉(zhuǎn)矩諧波幅值分布在16.5 mNm到19 mNm，即加工誤差對(duì)60次齒槽轉(zhuǎn)矩諧波即TNHC影響很小。

圖11 齒槽轉(zhuǎn)矩幅值分布

圖12 齒槽轉(zhuǎn)矩諧波分布

3 定子沖片循環(huán)旋轉(zhuǎn)

在定子沖片疊壓過(guò)程中，下一片沖片總是相對(duì)于前一片定向旋轉(zhuǎn)一個(gè)固定的槽數(shù)，如圖13所示。這樣，每一個(gè)定子沖片與轉(zhuǎn)子相互作用產(chǎn)生的齒槽轉(zhuǎn)矩都有一個(gè)相位差，通過(guò)旋轉(zhuǎn)合適的槽數(shù)，就能削弱特定極數(shù)倍數(shù)次的齒槽轉(zhuǎn)矩諧波。

圖13 定子沖片循環(huán)旋轉(zhuǎn)

對(duì)于任意槽極配合的電機(jī)，設(shè)其定子槽數(shù)為Z，轉(zhuǎn)子極數(shù)為P，定子側(cè)加工誤差導(dǎo)致的齒槽轉(zhuǎn)矩諧波為P的倍數(shù)次，P次諧波對(duì)應(yīng)幅值為KP。設(shè)單元電機(jī)數(shù)為Nu=GCD(P,Z)，即槽數(shù)與極數(shù)的最大公約數(shù)。記Z′=Z/Nu，P′=P/Nu，采用定子沖片循環(huán)旋轉(zhuǎn)方法，設(shè)每次旋轉(zhuǎn)一個(gè)槽對(duì)應(yīng)的角度，記第一片沖片的齒槽轉(zhuǎn)矩P次諧波為

TP次-第一片=KPsin(Pα)

(11)

根據(jù)正弦函數(shù)的性質(zhì)，相鄰Z′片定子沖片的齒槽轉(zhuǎn)矩P次諧波之和為

(12)

若Z′不是質(zhì)數(shù)，設(shè)a為Z′的一個(gè)因數(shù)，且Z′=a*b，定子沖片旋轉(zhuǎn)時(shí)，設(shè)每次旋轉(zhuǎn)a個(gè)槽對(duì)應(yīng)的角度，當(dāng)b不為1時(shí)，相鄰b片定子沖片的齒槽轉(zhuǎn)矩P次諧波之和為

(13)

由上可知，當(dāng)單元電機(jī)槽數(shù)不是質(zhì)數(shù)時(shí)，旋轉(zhuǎn)的槽數(shù)將有多種選擇，每片定子沖片旋轉(zhuǎn)a個(gè)槽，相鄰b片沖片的齒槽轉(zhuǎn)矩P次諧波就能相互抵消。同理，仍設(shè)定子沖片依次旋轉(zhuǎn)a個(gè)槽，c為正整數(shù)，相鄰b片定子沖片齒槽轉(zhuǎn)矩cP次諧波之和為

(14)

當(dāng)cP′/b不為整數(shù)時(shí)，相鄰b片定子沖片齒槽轉(zhuǎn)矩cP次諧波才能抵消。因此，b應(yīng)取最大值Z′，才能使得更多極數(shù)倍數(shù)次齒槽轉(zhuǎn)矩諧波被抵消，得到最佳的齒槽轉(zhuǎn)矩抑制效果，對(duì)應(yīng)定子沖片依次旋轉(zhuǎn)1個(gè)槽。

綜上所述，任意槽極配合的電機(jī)均能采用定子沖片循環(huán)旋轉(zhuǎn)方法抑制極數(shù)倍數(shù)次齒槽轉(zhuǎn)矩諧波，且定子沖片依次旋轉(zhuǎn)1個(gè)槽能取得最好的齒槽轉(zhuǎn)矩抑制效果。實(shí)際加工時(shí)也采用這種方案。

對(duì)于本文研究的12槽10極電機(jī)，定子側(cè)誤差產(chǎn)生的是10的倍數(shù)次齒槽轉(zhuǎn)矩諧波。單元電機(jī)槽數(shù)為6，極數(shù)為5，定子沖片依次旋轉(zhuǎn)1個(gè)槽，那么c可取值1、2、3、4和5，即相鄰6片定子沖片齒槽轉(zhuǎn)矩10、20、30、40和50次諧波之和都等于0。

對(duì)考慮加工誤差的有限元模型應(yīng)用定子沖片循環(huán)旋轉(zhuǎn)方法，并忽略沖片間的不一致性，得到齒槽轉(zhuǎn)矩幅值分布如圖14所示，齒槽轉(zhuǎn)矩諧波分布如圖15所示。定子沖片循環(huán)旋轉(zhuǎn)之后，齒槽轉(zhuǎn)矩幅值分布從21 mNm到52 mNm，數(shù)學(xué)期望為33 mNm，比普通定子電機(jī)的40 mNm減小了17.5%。齒槽轉(zhuǎn)矩10次諧波基本上消除了，12次、24次和60次沒(méi)有變化。

圖14 齒槽轉(zhuǎn)矩幅值分布(定子沖片旋轉(zhuǎn))

圖15 齒槽轉(zhuǎn)矩諧波分布(定子沖片旋轉(zhuǎn))

采用定子循環(huán)換位后和未采用循環(huán)換位的齒槽轉(zhuǎn)矩概率分布對(duì)比如圖16所示，可以看出通過(guò)循環(huán)換位后齒槽轉(zhuǎn)矩的分布整體向較小的幅值區(qū)間偏移，其數(shù)學(xué)期望相比換位前減小了17.5%，且在齒槽轉(zhuǎn)矩幅值較高的分布區(qū)間的頻數(shù)均有所減小，說(shuō)明采用定向旋轉(zhuǎn)沖片的方法可以在不提升加工精度的情況下，較好地抑制定子側(cè)加工誤差引入的齒槽轉(zhuǎn)矩諧波，進(jìn)而抑制總的齒槽轉(zhuǎn)矩；另外，循環(huán)旋轉(zhuǎn)疊壓的定子沖片內(nèi)圓呈波浪紋，清晰地顯示了沖片內(nèi)徑加工誤差。普通疊壓定子沖片內(nèi)圓光滑，即沖片旋轉(zhuǎn)之前相同定子齒下的沖片內(nèi)徑一致，這樣才能保證旋轉(zhuǎn)之后有較好的齒槽轉(zhuǎn)矩抑制效果。

圖16 優(yōu)化前后齒槽轉(zhuǎn)矩幅值對(duì)比

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

從生產(chǎn)的一批電機(jī)中隨機(jī)抽取樣本進(jìn)行齒槽轉(zhuǎn)矩測(cè)試，測(cè)試平臺(tái)如圖17所示。齒槽轉(zhuǎn)矩測(cè)試結(jié)果如圖18所示，齒槽轉(zhuǎn)矩幅值約為40.4 mNm，與定子沖片循環(huán)旋轉(zhuǎn)有限元模型計(jì)算的21 mNm至52 mNm相符。傅里葉分解得到的額外齒槽轉(zhuǎn)矩諧波主要是12次，與能量法模型和有限元模型的分析結(jié)果一致。極數(shù)倍數(shù)次齒槽轉(zhuǎn)矩諧波幅值均較小，驗(yàn)證了定子沖片循環(huán)旋轉(zhuǎn)方法的可行性。

圖17 齒槽轉(zhuǎn)矩測(cè)試平臺(tái)

圖18 實(shí)測(cè)齒槽轉(zhuǎn)矩

5 結(jié) 論

本文針對(duì)有加工誤差的表貼式永磁伺服電機(jī)進(jìn)行研究分析，主要得到如下結(jié)論：

(1)基于能量法，建立的計(jì)及加工誤差的表貼式永磁伺服電機(jī)齒槽轉(zhuǎn)矩解析模型，能分別評(píng)估齒槽轉(zhuǎn)矩對(duì)永磁體剩磁、厚度、寬度及周向位置等誤差的敏感性大小。解析計(jì)算結(jié)果表明，永磁體周向位置誤差對(duì)齒槽轉(zhuǎn)矩影響最大。實(shí)際生產(chǎn)中，若加強(qiáng)對(duì)該參數(shù)的加工精度，可取得較好的齒槽轉(zhuǎn)矩抑制效果。

(2)對(duì)任意槽極配合電機(jī)，均能采用定子沖片循環(huán)旋轉(zhuǎn)方法削弱定子側(cè)誤差產(chǎn)生的額外齒槽轉(zhuǎn)矩，并且相鄰沖片依次旋轉(zhuǎn)一個(gè)槽距能取得最好的齒槽轉(zhuǎn)矩抑制效果。

(3)基于蒙特卡羅方法，利用有限元軟件，建立的一套計(jì)及實(shí)際加工誤差的分析體系，能較好地預(yù)測(cè)批量生產(chǎn)永磁伺服電機(jī)齒槽轉(zhuǎn)矩及其頻譜分布規(guī)律，用于指導(dǎo)電機(jī)生產(chǎn)。