王 欽，楊慶東

（北京信息科技大學 機電工程學院，北京 100192）

0 引言

齒槽轉(zhuǎn)矩是產(chǎn)生轉(zhuǎn)矩波動的主要因素之一，數(shù)控磨床轉(zhuǎn)臺對電機轉(zhuǎn)矩波動有較高要求，因此削弱電機的齒槽轉(zhuǎn)矩對提高磨床轉(zhuǎn)臺的回轉(zhuǎn)精度具有重要意義。國內(nèi)外學者為削弱齒槽轉(zhuǎn)矩總結(jié)了許多方法［1－4］，本文基于能量法和傅里葉級數(shù)解析法提出了一種盤式力矩電機齒槽轉(zhuǎn)矩波動的計算方法，并利用有限元分析工具Ansoft Maxwell建立盤式電機模型和外轉(zhuǎn)子電機模型，通過對比分析驗證盤式力矩電機齒槽轉(zhuǎn)矩波動峰值較小的優(yōu)勢。

1 齒槽轉(zhuǎn)矩計算表達式

1.1 齒槽轉(zhuǎn)矩產(chǎn)生機理

盤式力矩電機主要磁場結(jié)構(gòu)如圖1所示，它由定子和轉(zhuǎn)子構(gòu)成。定子由沖卷機一邊沖孔一邊沿Z軸卷制而成，由于現(xiàn)有繞線技術的限制，定子必須要有齒槽，方便導線嵌入；轉(zhuǎn)子由永磁體和磁軛組成，每一塊磁鐵均為扇形，N極、S極交替排列。盤式力矩電機磁場方向如圖2所示，沿Z軸方向，即電機軸向方向。根據(jù)相互作用力原理，轉(zhuǎn)子和定子所受齒槽力大小相等、方向相反，因此分析定子鐵芯的受力情況即可得到電機的轉(zhuǎn)子受力情況，進而得到電機的齒槽轉(zhuǎn)矩。

轉(zhuǎn)子旋轉(zhuǎn)過程中，永磁體與定子齒槽產(chǎn)生非線性吸引，力F均勻分布在齒上，轉(zhuǎn)軸垂直于X－Y平面，單齒輸出齒槽轉(zhuǎn)矩T是由受力大小相同、半徑不同的扭矩疊加形成，電機的齒槽轉(zhuǎn)矩是由所有齒的齒槽轉(zhuǎn)矩疊加而成。齒槽轉(zhuǎn)矩波動峰值在高速運轉(zhuǎn)時表現(xiàn)為振動，當磨床轉(zhuǎn)臺小切削量精加工外圓面，齒槽轉(zhuǎn)矩波動峰值產(chǎn)生的振動與磨削顫振共振時，就會在磨削表面產(chǎn)生明顯波紋。

圖1 盤式力矩電機磁場結(jié)構(gòu)

圖2 盤式力矩電機磁場方向

1.2 盤式力矩電機齒槽轉(zhuǎn)矩表達式

根據(jù)能量法［5］，齒槽轉(zhuǎn)矩是由轉(zhuǎn)子旋轉(zhuǎn)時電機能量的變化引起的。忽略永磁體和鐵芯的能量變化，假定盤式力矩電機永磁體磁感線穿過的面為矩形，建立的齒槽轉(zhuǎn)矩表達式為：

其中：α為轉(zhuǎn)子位置，即永磁體中心線與齒中心線夾角；Wairgap為氣隙磁場能量；μ0為空氣磁導率；θ為轉(zhuǎn)子角度；G（θ）為氣隙相對磁導率；B（θ，α）為無槽時氣隙磁密分布；V為求解區(qū)域。

其中：Vairgap為氣隙體積；BanNL和GanNL均為傅里葉系數(shù)；n為一個整數(shù)；NL為Ns和Np的最小公倍數(shù)。

氣隙體積由下式計算：

其中：R1為定子內(nèi)半徑；R2為定子外半徑；g為氣隙。

根據(jù)磁路定律，計算式（5）中的氣隙磁密Bδ的公式為：

其中：Br為剩磁密度；kml為考慮永磁體間漏磁后的系數(shù)。計算kml的方程為：

其中：Pg為一個極下氣隙磁導；Pm為永磁體磁導；Pml為永磁體間的漏磁。上述幾個量的計算方程為：

其中：wm為永磁體寬度；LFe為永磁體長度；wf為永磁體間的寬度。

將式（9）～（11）代入式（8）得到kml，將式（8）代入式（7）得到Bδ，將式（7）代入式（5）得到BanNL，將式（4）～式（6）代入式（2）得到徑向磁場盤式力矩電機齒槽轉(zhuǎn)矩Tcog（α）的表達式：

從以上分析可以看出，氣隙體積Vairgap、氣隙磁密Bδ是影響齒槽轉(zhuǎn)矩波動幅值的重要因素；極數(shù)和槽數(shù)的最小公倍數(shù)NL、極弧系數(shù)αp是影響齒槽轉(zhuǎn)矩波動頻率的重要因素。在設計電機時應盡量減小氣隙g以降低齒槽轉(zhuǎn)矩波動最大值；在選取極槽配合參數(shù)時，應考慮電機工況，避開共振頻率。氣隙磁密Bδ是永磁體能量的表現(xiàn)形式，氣隙磁密越大，電機磁能越大，輸出扭矩的能力越大，通過優(yōu)化極弧系數(shù)在保證氣隙磁密正弦分布的同時也可改善齒槽轉(zhuǎn)矩波動。

由于本文采用的能量法解析計算模型將平底槽型簡化成矩形，假定氣隙磁密為矩形，磁通集中存在于齒部，實際上槽口處也存在磁通，使得計算結(jié)果偏小。

2 計算與仿真

為證明盤式力矩電機齒槽轉(zhuǎn)矩較小的優(yōu)勢，本文以外轉(zhuǎn)子力矩電機作為對比對象進行分析。電機主要參數(shù)如表1所示。

表1 電機主要參數(shù)表

將兩個電機參數(shù)代入式（12），計算結(jié)果如圖3所示。因盤式力矩電機的氣隙體積Vairgap較小，所以齒槽轉(zhuǎn)矩波動峰值較小。

使用有限元工具Ansofe Maxwell建立電機仿真模型。根據(jù)電磁場的對稱分布原理，計算一個對稱單元即可，即1/6個電機。因磁場方向為軸向，2D切面無法完整反映電磁分布，所以盤式力矩電機要用3D仿真模型進行計算，但外轉(zhuǎn)子電機使用2D仿真模型即可。有限元模型參數(shù)設定如表2所示。

圖3 齒槽轉(zhuǎn)矩表達式計算對比

表2 有限元模型參數(shù)設定表

通過多次仿真發(fā)現(xiàn)，將外轉(zhuǎn)子電機2D模型Max Length（最大長度）設定為0.5mm，Normal Deviation（正態(tài)偏差）設定為0.2°是最佳值，此時增大網(wǎng)格數(shù)量仿真波形變化不明顯，已經(jīng)達到了有限元計算的最佳精度。

3D運算量比2D大很多，為達到3D最佳運算精度，需找到電腦性能所能承受的最大值，經(jīng)過多次嘗試發(fā)現(xiàn)，盤式力矩電機Max Length（最大長度）設定為2mm，Normal Deviation（正態(tài)偏差）設定為5°～10°是極限值，如果再增大網(wǎng)格數(shù)量，電腦就會內(nèi)存不足無法運算，如果減小網(wǎng)格數(shù)量，波形趨勢基本不變，能夠滿足運算精度。齒槽轉(zhuǎn)矩有限元計算結(jié)果如圖4所示。

理論計算與有限元仿真結(jié)果對比如表3所示。從表3中可以看出，兩種計算方法得出的結(jié)果都顯示盤式力矩電機齒槽轉(zhuǎn)矩較小，理論計算結(jié)果顯示盤式力矩電機齒槽轉(zhuǎn)矩比外轉(zhuǎn)子力矩電機減小24.5%。有限元計算結(jié)果較為精確，結(jié)果顯示盤式力矩電機齒槽轉(zhuǎn)矩比外轉(zhuǎn)子力矩電機減小28.9%，η1與η2相差4.4%，由此證明式（12）能夠在一定程度上反映齒槽轉(zhuǎn)矩波動的規(guī)律，但是由于模型簡化較多，幅值誤差較大。

圖4 齒槽轉(zhuǎn)矩有限元計算結(jié)果

表3 理論計算與仿真結(jié)果對比表

3 結(jié)論

本文推導了盤式力矩電機齒槽轉(zhuǎn)矩波動表達式，建立有限元模型對盤式力矩電機和外轉(zhuǎn)子力矩電機的齒槽轉(zhuǎn)矩波動進行了對比仿真計算。仿真結(jié)果表明，在槽型結(jié)構(gòu)、永磁體間隙、定子沖片外形尺寸相同的條件下，盤式力矩電機齒槽轉(zhuǎn)矩波動幅值減小25%，氣隙體積Vairgap是影響齒槽轉(zhuǎn)矩波動峰值的重要因素。為提高轉(zhuǎn)臺加工精度，減小齒槽轉(zhuǎn)矩波動，可將盤式力矩電機嵌入轉(zhuǎn)臺內(nèi)，設計電機時應盡量減小氣隙，以減小轉(zhuǎn)矩波動峰值。

［1］劉婷.表貼式永磁同步電機齒槽轉(zhuǎn)矩削弱方法研究［D］.長沙：湖南大學，2013：7-12.

［2］鄧秋玲，黃守道，劉婷，等.永磁電機齒槽轉(zhuǎn)矩的研究分析［J］.湖南大學學報（自然科學版），2011（3）：56-59.

［3］楊玉波，王秀和，朱常青.基于分塊永磁磁極的永磁電機齒槽轉(zhuǎn)矩削弱方法［J］.電工技術學報，2012（3）：73-78.

［4］黃守道，劉婷，歐陽紅林，等.基于槽口偏移的永磁電機齒槽轉(zhuǎn)矩削弱方法［J］.電工技術學報，2013（3）：99-106.

［5］朱莉，姜淑忠，諸自強，等.表面式永磁電機齒槽轉(zhuǎn)矩解析模型比較［J］.微電機，2010（1）：10-15.