連軍強，孫宏坤

1 海軍駐桂林地區(qū)軍事代表室，廣西桂林 541000

2 海軍駐黃埔造船有限公司軍事代表室，廣東廣州 510715

0 引 言

隨著艦船綜合電力和稀土永磁技術(shù)的發(fā)展，永磁電機已被廣泛應(yīng)用于艦船推進領(lǐng)域。與傳統(tǒng)的勵磁電機相比，永磁電機具有低損耗、高功率密度等優(yōu)點。然而，齒槽轉(zhuǎn)矩作為永磁電機的固有特點會引起電機的振動、噪聲和控制困難；因此，削弱永磁電機齒槽轉(zhuǎn)矩成為永磁電機研究的重要方面。永磁電機齒槽轉(zhuǎn)矩的產(chǎn)生來自永磁體與電樞齒之間的切向力，是永磁體與電樞齒槽相互作用的結(jié)果［1-3］，它受電機極弧系數(shù)、永磁體厚度、氣隙長度、定子槽口寬度等多個設(shè)計參數(shù)的影響。這些設(shè)計參數(shù)與齒槽轉(zhuǎn)矩之間呈現(xiàn)出既相互影響又相互制約的關(guān)系。所以在齒槽轉(zhuǎn)矩參數(shù)優(yōu)化時，需要在允許的參數(shù)空間內(nèi)尋找一組最優(yōu)的參數(shù)。采用有限元方法可以精確分析計算永磁電機的齒槽轉(zhuǎn)矩，但該方法計算耗時長，且無法給出永磁電機設(shè)計參數(shù)與齒槽轉(zhuǎn)矩之間的顯性關(guān)系。本文擬提出一種基于粒子群算法的參數(shù)優(yōu)化方法。該方法將建立齒槽轉(zhuǎn)矩解析模型，并根據(jù)參數(shù)實際的取值范圍確定約束條件，以齒槽轉(zhuǎn)矩最小為目標(biāo)函數(shù)，使用粒子群算法對該解析模型在設(shè)計約束的多維參數(shù)空間內(nèi)進行參數(shù)尋優(yōu)，從而得到最優(yōu)的參數(shù)組合。

1 齒槽轉(zhuǎn)矩解析模型及改進

目前，齒槽轉(zhuǎn)矩的計算方法可以分為兩類：解析方法和數(shù)值方法。其中，解析法的物理概念清晰。齒槽轉(zhuǎn)矩的解析模型很多，主要分為基于能量法的齒槽轉(zhuǎn)矩解析模型和基于麥克斯韋張量法的齒槽轉(zhuǎn)矩解析模型兩大類［4-7］。齒槽轉(zhuǎn)矩參數(shù)優(yōu)化對齒槽轉(zhuǎn)矩解析模型有以下要求：

1）解析模型中要包含所有需要優(yōu)化的設(shè)計參數(shù)，這樣才能求得全局最優(yōu)解；

2）解析模型要能精確計算齒槽轉(zhuǎn)矩，且計算過程耗時要短，以減少粒子群尋優(yōu)的時間。

由于基于能量法的齒槽轉(zhuǎn)矩解析模型能直觀地體現(xiàn)各設(shè)計參數(shù)與齒槽轉(zhuǎn)矩的關(guān)系，且求解簡單，所以選取基于能量法的解析模型，在參數(shù)選擇時，取對齒槽轉(zhuǎn)矩影響最大的極弧系數(shù)ap、永磁體厚度hm和定子槽口寬度b等參數(shù)為優(yōu)化設(shè)計參數(shù)。齒槽轉(zhuǎn)矩參數(shù)優(yōu)化可描述為：

由于采用基于能量法的齒槽轉(zhuǎn)矩解析模型，為便于計算，根據(jù)能量法［8］，齒槽轉(zhuǎn)矩解析模型可進行以下化簡：

式中：θ為永磁體相對某一指定齒的中心線旋轉(zhuǎn)的角度；Wg(θ)為電機在角 θ時的儲能；G(α)為定轉(zhuǎn)子相對機械角為α的氣隙磁導(dǎo)；B(α，θ)為氣隙磁密；NL為磁極數(shù)與定子齒數(shù)的最小公倍數(shù)；n=1，2，…，N ；GnNL為磁導(dǎo)的NL的整數(shù)次諧波；為磁密的NL的整數(shù)次諧波。

從上式可以看出，建立齒槽轉(zhuǎn)矩解析模型的關(guān)鍵是推導(dǎo)出氣隙相對磁導(dǎo)GnNL和氣隙磁密BnNL的表達式，氣隙相對磁導(dǎo)GnNL的表達式為

現(xiàn)對BnNL進行推導(dǎo)。在齒槽轉(zhuǎn)矩模型中，一般都是將氣隙磁密理想化為矩形，根據(jù)拉普拉斯和準(zhǔn)泊松方程可得到精確的氣隙磁密表達式［9-12］。本文選取徑向充磁表貼式永磁電機為研究對象，在圖1所示的表貼式永磁電機模型中，Rs為定子半徑，Rr為轉(zhuǎn)子鐵心半徑，Rm為轉(zhuǎn)子外半徑，永磁體厚度表示為：hm=Rm-Rr，氣隙軸線的位置選在r==，對氣隙磁密徑向分量表達式進行變形，得

式中：Mn為磁化強度矢量；μ0為真空磁導(dǎo)率；μr為永磁體相對磁導(dǎo)率；Np為轉(zhuǎn)子極對數(shù)。

在一個磁極范圍內(nèi)對式（4）積分，得

將式（5）和式（3）代入式（2），得到齒槽轉(zhuǎn)矩的解析模型。解析模型中只有極弧系數(shù)ap、永磁體厚度hm和定子槽口寬度b為未知參數(shù)，其余均為已知量。

2 粒子群優(yōu)化算法

粒子群算法（Particle Swarm Optimization，PSO）最早由Eberhart和Kennedy于1995年提出，基本概念源于對鳥群覓食行為的研究［13-15］。

假設(shè)在一個D維的目標(biāo)搜索空間中，有N個粒子組成一個群落，其中第i個粒子表示為一個D維的向量：

第i個粒子的“飛行”速度也是一個D維的向量，記為

第i個粒子迄今為止搜索到的最優(yōu)位置稱為個體極值，記為

整個粒子群迄今為止搜索到的最優(yōu)位置為全局極值，記為

在找到這兩個最優(yōu)值時，粒子根據(jù)式（10）和式（11）來更新自己的速度和位置。

式中：c1和c2為學(xué)習(xí)因子，也稱加速常數(shù)（acceleration constant）；r1和 r2為［0，1］范圍內(nèi)的均勻隨機數(shù)［16］。算法流程如圖2所示。

圖2 PSO算法流程圖Fig.2 The flow chart of PSO

3 齒槽轉(zhuǎn)矩解析模型參數(shù)尋優(yōu)及優(yōu)化結(jié)果分析

以圖1所示的表貼式永磁電機為研究對象，電機軸長94 mm，電機參數(shù)見表1。永磁電機極弧系數(shù)ap的取值范圍[0 .6，0.9]，永磁體厚度hm（單位：mm）的取值范圍[3 .5，6.5]，定子槽口寬度b（單位：mm）的取值范圍[2 .25，2.75] ，解析模型中傅里葉級數(shù)迭代次數(shù)取20。將電機已知參數(shù)代入式（4）并進行化簡，可得到電機徑向磁場強度表達式：

表1 電機參數(shù)Tab.1 Parameters of permanent magnet motor

由式（2）、式（3）、式（4）和式（12），可得到齒槽轉(zhuǎn)矩的解析表達式

由式（13）所得的計算結(jié)果為沿氣隙軸線分布的多個點的齒槽轉(zhuǎn)矩值，而在電機設(shè)計過程中是將齒槽轉(zhuǎn)矩的峰值作為性能指標(biāo)，所以優(yōu)化目標(biāo)函數(shù)為

使用Matlab軟件對齒槽轉(zhuǎn)矩解析模型（式（14））進行編程，粒子群優(yōu)化算法的參數(shù)選取為：種群數(shù)為50；最大迭代次數(shù)分別取50，100和200。圖3所示為參數(shù)優(yōu)化過程和優(yōu)化函數(shù)收斂曲線。其中：圖3（a）和圖3（c）為迭代次數(shù)分別取50和200時，表貼式永磁電機3個對齒槽轉(zhuǎn)矩影響較大的參數(shù)，即極弧系數(shù)ap、永磁體厚度hm和定子槽口寬度b在粒子群尋優(yōu)過程中的變化；圖3（b）和圖3（d）為優(yōu)化目標(biāo)函數(shù)，即電機齒槽轉(zhuǎn)矩峰值在迭代過程中的收斂曲線。迭代50次時，目標(biāo)函數(shù)值為0.23 N·m；粒子群算法迭代100次時，目標(biāo)函數(shù)值為0.13 N·m；粒子群算法迭代200次時，目標(biāo)函數(shù)值仍為0.13 N·m；所以，粒子群算法在迭代100次已尋到全局最優(yōu)值。

4 優(yōu)化方法的有限元驗證

Ansoft軟件為專門的電機有限元計算軟件，使用Ansoft軟件自帶的優(yōu)化功能對文中所提出的優(yōu)化算法進行驗證，選擇Optimetrics選項中的GA選項，即遺傳算法優(yōu)化選項，Ansoft軟件中GA選項的數(shù)據(jù)樣本是大量的有限元分析數(shù)據(jù)，因此優(yōu)化結(jié)果可信，可作為本文所提出的優(yōu)化算法的驗證手段。GA選項的參數(shù)設(shè)置如圖4所示，遺傳算法尋優(yōu)圖和優(yōu)化結(jié)果分別如圖5和圖6所示，迭代52次時取得全局最優(yōu)值。

圖3 參數(shù)優(yōu)化過程和優(yōu)化函數(shù)收斂曲線Fig.3 Progress of optimization of parameters and optimization function convergence

粒子群優(yōu)化結(jié)果與Ansoft軟件中GA選項優(yōu)化結(jié)果的對比見表2。從中可見，粒子群優(yōu)化后齒槽轉(zhuǎn)矩值為0.13 N·m，遺傳算法優(yōu)化后齒槽轉(zhuǎn)矩值為0.132 N·m，兩者非常接近，但本文所提優(yōu)化算法的優(yōu)化時間要比Ansoft軟件中GA選項優(yōu)化時間少得多。這是由于基于粒子群算法的齒槽轉(zhuǎn)矩優(yōu)化方法是對齒槽轉(zhuǎn)矩解析模型進行優(yōu)化，它的樣本數(shù)據(jù)由解析模型計算得到，而Ansoft軟件中GA選項的樣本數(shù)據(jù)是由有限元計算得到的，耗時較長。

圖4 遺傳算法參數(shù)設(shè)置Fig.4 The parameter setup of Genetic Algorithm

圖5 遺傳算法尋優(yōu)圖Fig.5 The progress of Genetic Algorithm search for optimistic result

圖6 遺傳算法優(yōu)化結(jié)果Fig.6 The optimization result by using Genetic Algorithm

表2 優(yōu)化結(jié)果對比Tab.2 Comparison of the optimization results

使用粒子群優(yōu)化后的參數(shù)建立表貼式永磁電機的有限元模型，進行有限元分析，結(jié)果如圖7所示。在進行仿真時，設(shè)置電機角速度為1（°）/s，故每秒對應(yīng)機械度1°。電機參數(shù)優(yōu)化前齒槽轉(zhuǎn)矩值為1.23 N·m，使用粒子群算法對電機參數(shù)優(yōu)化后，齒槽轉(zhuǎn)矩值降為0.13 N·m，降低了近90%，說明基于粒子群算法的齒槽轉(zhuǎn)矩優(yōu)化方法能夠有效對電機參數(shù)進行優(yōu)化，達到大幅度削弱齒槽轉(zhuǎn)矩的目的。

5 結(jié)語

圖7 PSO算法優(yōu)化前后齒槽轉(zhuǎn)矩對比圖Fig.7 Cogging torque before optimization and after optimization

本文在對齒槽轉(zhuǎn)矩解析模型進行變形和簡化的基礎(chǔ)上，選取對齒槽轉(zhuǎn)矩影響較大的3個參數(shù)作為優(yōu)化變量，以齒槽轉(zhuǎn)矩峰值為優(yōu)化目標(biāo)量，使用粒子群算法進行優(yōu)化，經(jīng)有限元驗證，該方法取得了較精確的優(yōu)化參數(shù)，且計算時間大大縮短，說明基于粒子群算法的齒槽轉(zhuǎn)矩優(yōu)化方法能夠有效對電機參數(shù)進行優(yōu)化，達到大幅度削弱齒槽轉(zhuǎn)矩的目的。

［1］SWANN S A.Effect of rotor eccentricity on the magnetic field in the air-gap of a non-salient-pole machine［J］.Proceedings of the Institution of Electrical Engineers，1963，110（5）：903-915.

［2］LORIMER W，HARTMAN A.Magnetically induced vibrations in brushless DC magnet motors［C］//28th Annual Symposium Proceedings of Incremental Motion Control Systems and Devices，1999：299-308.

［3］冀溥，王秀和，王道涵，等.轉(zhuǎn)子靜態(tài)偏心的表面永磁電機齒槽轉(zhuǎn)矩研究［J］.中國電機工程學(xué)報，2004，24（9）：188-191.JI Pu，WANG Xiuhe，WANG Daohan，et al.Study of cogging torque in surface-mounted permanent magnet motors with static eccentricity［J］.Proceedings of the CSEE，2004，24（9）：188-191.

［4］張冉，王秀和，楊玉波，等.基于等效剩磁法的永磁電動機轉(zhuǎn)子偏心磁場解析計算［J］.電工技術(shù)學(xué)報，2009，24（5）：7-12.ZHANG Ran，WANG Xiuhe，YANG Yubo，et al.Analytical prediction of magnetic field in permanent magnet motors with rotor eccentricity based on the method of equivalent remanence［J］.Transactions of China Electronics Society，2009，24（5）：7-12.

［5］HARTMAN A，LORIMER W.Undriven vibrations in brushless DC motors［J］.IEEE Transactions on Magnetics，2001，37（2）：789-792.

［6］章名濤，肖如鴻.電機的電磁場［M］.北京：機械工業(yè)出版社，1988：182-201.

［7］HWANG S M，KIM K T，JEONG W B，et al.Comparison of vibration sources between symmetric and asymmetric HDD spindle motors with rotor eccentricity［J］.IEEE Transactions on Industry Applications，2001，37（6）：1727-1731.

［8］KIM U，LIEU D K.Effects of magnetically induced vibration force in brushless permanent-magnet motors［J］.IEEE Transactions on Magnetics，2005，41（6）：2164-2172.

［9］LI J T，LIU Z J，NAY L H A.Effect of radial magnetic forces in permanent magnet motors with rotor eccentricity［J］.IEEE Transactions on Magnetics，2007，43（6）：2525-2527.

［10］NIAZI P，TOLIYAT H A.Online parameter estimation of permanent magnet assisted synchronous reluctance motor［J］.IEEE Transactions on Industry Applications，2007，43（2）：609-615.

［11］BURTH M，VéLEZ-REYES G C.Subset selection for improved parameter estimation on on-line identification on a synchronous generator［J］.IEEE Transactions on Power Systems，1999，14（1）：218-225.

［12］程樹康，宮海龍，柴鳳，等.永磁輪轂電機齒槽轉(zhuǎn)矩研究［J］.中國電機 工程學(xué)報，2009，29（30）：47-51.CHENG Shukang，GONG Hailong，CHAI Feng，et al.Study on cogging torque of permanent magnet in-wheel motor［J］.Proceedings of the CSEE，2009，29（30）：47-51.

［13］KATHIRAVAN R，GANGULI R.Strength design of composite beam using gradient and particle swarm optimization［J］.Composite Structures，2007，81（4）：471-479

［14］FAN S K S，ZAHARA E.A hybrid simplex search and particle swarm optimization for unconstrained optimization［J］.European Journal of Operational Research，2007，181（2）：527-548.

［15］DE FALCO I，DELLA C A，TARANTINO E.Facing classification problems with particle swarm optimization［J］.Applied Soft Computing，2007，7（3）：652-658.

［16］阮旻智，李慶民，王紅軍，等.人工免疫粒子群算法在系統(tǒng)可靠性優(yōu)化中的應(yīng)用［J］.控制理論與應(yīng)用，2010，27（9）：1253-1258.RUAN Minzhi，LI Qingmin，WANG Hongjun，et al.Application of artificial immune particle swarm optimization algorithm to system-reliability optimization［J］.Control Theory and Application，2010，27（9）：1253-1258.