武志勇，郭宏，錢浩

（北京航空航天大學(xué)自動化科學(xué)與電氣工程學(xué)院，北京 100191）

考慮噪聲因子的永磁同步電機轉(zhuǎn)矩脈動分布特性分析

武志勇，郭宏，錢浩

（北京航空航天大學(xué)自動化科學(xué)與電氣工程學(xué)院，北京 100191）

針對永磁同步電機在批量生產(chǎn)過程中，受到電機尺寸加工偏差及材料屬性偏差等噪聲因子的影響，存在批量生產(chǎn)電機的轉(zhuǎn)矩脈動峰值會產(chǎn)生分散性的問題。為了研究存在噪聲因子條件下轉(zhuǎn)矩脈動峰值的分布特性，通過采用敏感性分析的方法得到了影響轉(zhuǎn)矩脈動峰值的主要噪聲因子包括磁鋼形位公差和材料屬性偏差，槽口尺寸偏差及氣隙尺寸偏差。通過理論分析和解析計算得到了存在噪聲因子的情況下，批量生產(chǎn)的電機產(chǎn)品轉(zhuǎn)矩脈動峰值的統(tǒng)計分布規(guī)律。得出在電機可行設(shè)計域內(nèi)，不同電機設(shè)計方案，其批量產(chǎn)品轉(zhuǎn)矩脈動峰值的均值和標(biāo)準(zhǔn)差呈正比關(guān)系，并通過粒子群優(yōu)化結(jié)果證明了結(jié)論的有效性。

永磁;同步電機;敏感性分析;統(tǒng)計分析;轉(zhuǎn)矩

0 引言

永磁同步電動機具有反電勢正弦度好、控制性能優(yōu)良、等諸多優(yōu)點，廣泛應(yīng)用于高精度電機伺服系統(tǒng)內(nèi)。限制永磁同步電機整體性能提升的關(guān)鍵問題在于電機的轉(zhuǎn)矩脈動［1－4］。轉(zhuǎn)矩脈動作為高精密伺服系統(tǒng)的寄生特性，一直受到人們的重視。國內(nèi)外學(xué)者對轉(zhuǎn)矩脈動開展了深入的機理分析并研究從電機設(shè)計層面對轉(zhuǎn)矩脈動峰值進行有效的抑制。

但是，在實際生產(chǎn)過程中，產(chǎn)生了如下兩個問題，其一，由于加工過程存在不可抑制的結(jié)構(gòu)尺寸公差和材料屬性偏差等噪聲因子，而轉(zhuǎn)矩波動又對這些噪聲因子異常敏感，進而造成了實際加工之后，電機轉(zhuǎn)矩脈動值偏離了設(shè)計期望值。其二，由于在批量生產(chǎn)的電機中，每臺電機的噪聲因子情況均不一致，從而造成了在批量生產(chǎn)電機內(nèi)部，轉(zhuǎn)矩脈動峰值呈現(xiàn)一定的分散性。這就使得初始的電機設(shè)計期望值無法正確的反映出批量電機內(nèi)部轉(zhuǎn)矩脈動的特性。進而影響了伺服系統(tǒng)的控制性能。

為了研究轉(zhuǎn)矩脈動峰值在批量生產(chǎn)電機內(nèi)部的分布規(guī)律，一些文獻針對此問題開展了分析。文獻［5－6］重點分析了磁鋼加工偏差條件下，齒槽轉(zhuǎn)矩峰值在批量電機內(nèi)部的分布情況，文獻［7］對由磁鋼偏差引起的反電勢有效值的分散性進行統(tǒng)計分析。文獻［8］針對轉(zhuǎn)矩脈動在存在加工偏差的條件下，進行電機結(jié)構(gòu)尺寸優(yōu)化，提升電機系統(tǒng)的穩(wěn)健性。

從上述文獻可以總結(jié)出，首先，由于電機性能的分散性由加工公差引起，因此在分析批量電機關(guān)鍵特性的分散性時，文獻均針對影響性能分散性最甚的偏差進行處理。比如磁鋼偏差等。其次，這些文獻均采用了統(tǒng)計的手段對分散性進行描述。這種統(tǒng)計學(xué)規(guī)律往往可以由樣本的均值和標(biāo)準(zhǔn)差來表示。均值反映了批次電機整體性能與設(shè)計目標(biāo)值所接近的程度，標(biāo)準(zhǔn)差反映的為批次電機的性能分散性。通過這兩個統(tǒng)計學(xué)變量，能夠較為直觀的反映出批量電機特性的整體分布規(guī)律。

在解決轉(zhuǎn)矩脈動分布問題中，均值和標(biāo)準(zhǔn)差之間的關(guān)系值得進一步深入研究，其反映了批量電機關(guān)鍵特性的分布情況，對兩個變量關(guān)系的研究也是對電機性能分布的研究。通過對均值和標(biāo)準(zhǔn)差的統(tǒng)計分析，可以得到在一定的公差條件下，批量電機在制造過程中所呈現(xiàn)的分布規(guī)律。目前，文獻已對齒槽定位轉(zhuǎn)矩和反電勢波形進行了有效的分析，還沒有針對轉(zhuǎn)矩脈動分散性的有效統(tǒng)計分析。本文通過理論分析和解析計算討論了各種加工偏差對均值和標(biāo)準(zhǔn)差的影響程度，研究了在可行域內(nèi)考慮加工偏差等噪聲因子后電機轉(zhuǎn)矩脈動特性的均值和標(biāo)準(zhǔn)差之間的關(guān)系。

本文首先針對電機制造過程中產(chǎn)生的電機尺寸和形位公差進行敏感性分析，得到對于電機關(guān)鍵特性較為敏感的噪聲因子。之后針對幾種典型的電機制造公差，對批次電機在可行域內(nèi)均值和標(biāo)準(zhǔn)差的分布規(guī)律進行理論推導(dǎo)和解析計算，得到了均值和標(biāo)準(zhǔn)差之間的數(shù)值關(guān)系。最后采用粒子群算法分別將均值和標(biāo)準(zhǔn)差作為優(yōu)化目標(biāo)進行優(yōu)化，通過優(yōu)化結(jié)果證明了結(jié)論的有效性。

1 噪聲因子敏感性分析

對于噪聲因子的敏感性分析是轉(zhuǎn)矩脈動統(tǒng)計分析的前提。在電機加工過程中，噪聲因子主要包括以下幾種:

1）定子和轉(zhuǎn)子部分:槽開口誤差（Δbo），氣隙偏差（Δa－g），轉(zhuǎn)子偏心（ΔInn）。

2）磁鋼部分:磁鋼厚度（Δw），磁鋼貼偏（Δδ），剩磁偏差（Δbr），充磁角度偏差（Δθ）。

噪聲因子示意圖如圖1所示。

圖1 電機各種偏差示意圖Fig.1The tolerances of the motors

敏感性分析主要通過計算含噪聲因子情況下與理想無噪聲因子條件下轉(zhuǎn)矩脈動峰值的比值，得到轉(zhuǎn)矩脈動對于各噪聲因子的敏感程度。

本文選擇兩種較為典型的分?jǐn)?shù)槽集中繞組永磁同步電機作為敏感性分析的對象，分別是14極18槽電機以及22極24槽電機。這兩種分?jǐn)?shù)槽集中繞組永磁同步電機的反電勢正弦度好，能夠抑制由于反電勢非正弦造成的轉(zhuǎn)矩脈動，同時兩種電機的槽數(shù)和磁鋼數(shù)的公倍數(shù)較大也相對削弱了電機的齒槽轉(zhuǎn)矩。由于這兩種極槽配合的電機轉(zhuǎn)矩脈動值較小，噪聲因子會導(dǎo)致實際生產(chǎn)值嚴(yán)重偏離設(shè)計期望值，造成了設(shè)計和生產(chǎn)結(jié)果之間的很大的誤差。因此采用這種低轉(zhuǎn)矩脈動的永磁同步電機極槽配合結(jié)構(gòu)進行敏感性分析是十分必要的。

電機的初始設(shè)計方案已經(jīng)完成。具體電動機參數(shù)見表1。

表1 電機初始化設(shè)計方案Table 1Motor initial design parameter

由于電機轉(zhuǎn)矩脈動可分為電磁轉(zhuǎn)矩脈動（Tem）以及齒槽轉(zhuǎn)矩（Tcogg）兩部分，因此對轉(zhuǎn)矩脈動這兩個部分分別進行敏感性分析。

在敏感性分析中噪聲因子不同尺度的公差帶也會對敏感性產(chǎn)生影響，為了避免由于公差帶選擇過大或過小干擾敏感性分析的結(jié)果，14極18槽電機的部分噪聲因子選擇2～3個公差帶，綜合分析在不同公差標(biāo)準(zhǔn)的情況下，噪聲因子的敏感程度，所有公差帶的選取參考國標(biāo)GB/T 1800.3要求。敏感性分析采用有限元方法進行。在有限元計算時，針對每種電機模型均進行3次網(wǎng)格剖分，并檢查計算結(jié)果誤差在1%范圍之內(nèi)，從而避免了由于有限元時空離散對計算結(jié)果所造成的影響。有限元結(jié)果如表2，3所示，表中的百分比是指含噪聲因子時與理想無噪聲因子時相比，轉(zhuǎn)矩脈動峰值相對增加的百分比。

根據(jù)敏感性分析結(jié)果可以得出:Tem和Tcogg對于Δδ，ΔBr，Δθ三個噪聲因子異常敏感，噪聲因子的微小偏差會造成轉(zhuǎn)矩脈動特別是齒槽轉(zhuǎn)矩的強烈變化。同時Δbo以及Δa_g對于齒槽轉(zhuǎn)矩也較為敏感，因此根據(jù)以上結(jié)論，在統(tǒng)計分析時針對如下幾個噪聲因子進行:Δδ，ΔBr，Δθ，Δbo以及Δa－g。以上分析是依據(jù)設(shè)計資料和計算結(jié)果得到。

表2 電機定子和轉(zhuǎn)子部分尺寸偏差的敏感性分析Table 2Sensitive analysis of the size errors of the stators and rotors

表3 電機磁鋼部分尺寸偏差的敏感性分析Table 3Sensitive analysis of the size errors of the magnets

2 批量生產(chǎn)電機轉(zhuǎn)矩脈動峰值均值和標(biāo)準(zhǔn)差分析

本節(jié)分析電機在各種噪聲因子作用下，批量生產(chǎn)電機轉(zhuǎn)矩脈動峰值的統(tǒng)計規(guī)律以及均值和標(biāo)準(zhǔn)差之間的關(guān)系。根據(jù)敏感性分析的結(jié)果，轉(zhuǎn)矩脈動峰值對磁鋼類噪聲因子最為敏感，因此主要對磁鋼類噪聲因子造成的影響進行分析。首先通過理論推導(dǎo)對批量電機轉(zhuǎn)矩脈動均值和標(biāo)準(zhǔn)差之間的關(guān)系進行說明，隨后采用解析計算方法對批量生產(chǎn)電機轉(zhuǎn)矩脈動的統(tǒng)計分布進行解析計算。本節(jié)計算所使用的電機案例為上節(jié)所述的14極18槽電機。

2.1 磁鋼類噪聲因子

主要針對Δδ，ΔBr，Δθ這三種磁鋼類噪聲因子展開研究。

根據(jù)文獻［5－6］可知，電機的轉(zhuǎn)矩脈動可分為固有的轉(zhuǎn)矩脈動（TNHC）和附加轉(zhuǎn)矩脈動（TAHC）。其中，TNHC為電機不存在噪聲因子時轉(zhuǎn)矩脈動的固有部分。而TAHC為考慮噪聲因子作用后，附加的轉(zhuǎn)矩脈動部分。磁鋼偏差如圖2所示;

圖214 p18s電機磁鋼偏差示意圖，每塊磁鋼包含磁鋼貼偏，剩磁偏差，以及充磁角度偏差.Fig.2Tolerances of the 14 poles and 18 slots motors，including the magnet misplacements，Br errors and magnetization directions errors

當(dāng)磁鋼的偏差遠(yuǎn)小于磁鋼名義尺寸時，計算由磁鋼偏差所產(chǎn)生的TAHC可進行兩種近似處理。

1）由各噪聲因子作用產(chǎn)生的TAHC與噪聲因子數(shù)值之間呈線性關(guān)系;

2）各塊磁鋼產(chǎn)生的TAHC之間滿足相移關(guān)系;

文獻［5，7］針對以上兩個近似處理有詳細(xì)的理論證明。當(dāng)不考慮磁路飽和時，上面的結(jié)論嚴(yán)格成立，當(dāng)考慮磁路飽和時，結(jié)論近似成立。

根據(jù)以上的近似，首先得到一塊磁鋼單位偏差產(chǎn)生的TAHC的波形，通過相移和比例縮放及線性加權(quán)可得到各塊磁鋼噪聲因子總作用下的TAHC波形。

在電機加工制造過程中，上述三種磁鋼類噪聲因子在實際工程加工中常按標(biāo)準(zhǔn)正態(tài)分布處理［5，7］。其分布規(guī)律如下所示:

其中，σerr為三種磁鋼誤差的標(biāo)準(zhǔn)差。將σerr數(shù)值作為單位偏差數(shù)值。根據(jù)單位偏差數(shù)值對噪聲因子進行歸一化處理。

其中，σδ，σθ，σbr分別為三種公差的標(biāo)準(zhǔn)差。因此，Δδ'，Δθ'，Δbr'∈［～，N（0，1）］。

由于三種噪聲因子的推導(dǎo)過程一致，因此在計算時由Δerr'統(tǒng)一代表Δδ'，Δθ'，Δbr'不再單獨列寫表達(dá)式。

通過以上的分析，當(dāng)各磁鋼存在噪聲因子時，噪聲因子產(chǎn)生的TAHC可由傅里葉級數(shù)展開式表示為:

在分析過程中，各個頻次的幅值是考慮的核心問題。其表達(dá)式可以寫為:

當(dāng)h≠nNm，n=1，2，3…時為一閉合矢量和;其向量大約有一半的長度在實軸上，一半在虛軸上。因此

由于上式的實部和虛部兩部分為標(biāo)準(zhǔn)差相同的正態(tài)分布，因此向量的模值服從瑞利分布。瑞利分布的期望為，標(biāo)準(zhǔn)差為

因此，當(dāng)h≠nNm，n=1，2，3…

可以看到，磁鋼誤差導(dǎo)致的轉(zhuǎn)矩脈動分散性，其期望和標(biāo)準(zhǔn)差都與電機本身的電磁環(huán)境（h）相關(guān)，并且還與磁鋼的塊數(shù)以及磁鋼偏差的分布相關(guān)。根據(jù)上述推導(dǎo)可知:

從式11中可以看到，附加轉(zhuǎn)矩脈動峰值在大部分頻次的期望和標(biāo)準(zhǔn)差均成正比，而在一些頻次上均值和標(biāo)準(zhǔn)差的規(guī)律則有變動。轉(zhuǎn)矩脈動峰值是對各頻次的峰值進行矢量求和的結(jié)果。因此轉(zhuǎn)矩脈動的峰值會由一些含有較大轉(zhuǎn)矩脈動均值和標(biāo)準(zhǔn)差的頻次所主導(dǎo)。雖然期望和標(biāo)準(zhǔn)差的比值會有所改變，但其整體上應(yīng)仍呈現(xiàn)正比關(guān)系。

上述是針對一種噪聲因子獨立作用時的分析結(jié)論，當(dāng)三種噪聲因子共同作用時，由于三種偏差的幅值遠(yuǎn)小于磁鋼尺寸和材料屬性名義值，其各自作用產(chǎn)生的TAHC可看做相互獨立的隨機變量，當(dāng)加權(quán)三種噪聲因子的計算結(jié)果后，不改變均值和標(biāo)準(zhǔn)差之間的比例關(guān)系。

因此可得出當(dāng)電機存在磁鋼類噪聲因子時，其轉(zhuǎn)矩脈動的分布規(guī)律服從均值和標(biāo)準(zhǔn)差呈正比的正態(tài)分布。

為了進一步驗證得到的結(jié)論，通過仿真模擬實際批量生產(chǎn)的電機轉(zhuǎn)矩脈動的分布情況對上述理論推導(dǎo)進行驗證。為了符合實際工況的需要，對電機的可行設(shè)計域進行限定，其可行設(shè)計域的上下限如表4所示。在電機可行設(shè)計域內(nèi)對電機的各尺寸進行隨機組合，選取50個電機設(shè)計方案作為實驗對象，每個實驗對象在磁鋼公差帶范圍內(nèi)隨機生成1 000臺電機樣本，其中的每塊磁鋼三種噪聲因子均嚴(yán)格服從正態(tài)分布。

盡管有限元計算能夠有效的得到精確的轉(zhuǎn)矩脈動波形，但采用有限元方法對此問題進行計算時會消耗大量的計算時間。因此，本文采用建立解析模型方法代替有限元方法用以計算轉(zhuǎn)矩脈動特性。文獻［9－10］提出的Exact Subdomain方法在電機磁路弱飽和時能夠精確的得出電機的電磁分布規(guī)律，本文針對Exact Subdomain方法對其進行改良，使其能夠有效的計算加工偏差對轉(zhuǎn)矩脈動的影響。雖然Exact Subdomain方法是基于電機開路的電磁解析計算方法，在處理帶載問題時會受到磁路飽和的影響產(chǎn)生一定的計算誤差，但針對統(tǒng)計學(xué)的計算，由于噪聲因子對電機磁路飽和程度的影響非常小，在含有不同噪聲因子的各電機內(nèi)，由飽和所造成的計算偏差程度基本一致，因此飽和對于均值和方差的影響可近似視為相同且對均值和方差的整體關(guān)系不構(gòu)成顛覆性的影響。因此采用Exact Subdomain方法能夠有效的反映兩者之間的數(shù)值關(guān)系。在計算時重點考慮了磁鋼貼偏，磁鋼的剩磁誤差以及充磁方向偏差情況下的磁場分布。

表4 電機設(shè)計變量可行域范圍及噪聲因子范圍Table 4Range of the parameters and the noise factors

采用解析方法對50個電機設(shè)計方案的1 000臺含噪聲因子的樣本進行轉(zhuǎn)矩脈動計算，分別得到各實驗對象轉(zhuǎn)矩脈動峰值均值和標(biāo)準(zhǔn)差值。

通過解析計算，50個電機設(shè)計方案轉(zhuǎn)矩脈動峰值的均值和標(biāo)準(zhǔn)差分布規(guī)律如下圖所示:根據(jù)圖3所示，各電機設(shè)計方案的均值和標(biāo)準(zhǔn)差基本呈現(xiàn)正比關(guān)系。

圖3 50個電機設(shè)計方案轉(zhuǎn)矩脈動均值和標(biāo)準(zhǔn)差Fig.3Expected values and the standard deviation of the torque ripple peak values in 50 motors

2.2 定子噪聲因子

關(guān)于定子類噪聲因子，主要針對Δbo和Δa_g展開研究，如圖4所示。根據(jù)敏感性分析的結(jié)論，兩種噪聲因子對轉(zhuǎn)矩脈動峰值的影響要遠(yuǎn)小于磁鋼類噪聲因子的影響。因此這里不做理論推導(dǎo)，而是通過解析計算的方法得到存在兩種偏差情況下轉(zhuǎn)矩脈動的分布特性。

圖4 14p18s電機槽口偏差以及氣隙偏差示意圖Fig.4Slots and air gap errors of the 14p18s

在電機設(shè)計可行域之內(nèi)，選擇50個電機設(shè)計方案，每個實驗對象在定子偏差公差帶范圍內(nèi)隨機生成1 000臺樣本，得到50個電機設(shè)計方案均值和標(biāo)準(zhǔn)差的分布圖如圖6、圖7所示。圖5為槽口偏差時的分布圖，圖6為氣隙誤差的分布圖。從圖中可以看到槽口偏差所引起的性能分散性，轉(zhuǎn)矩脈動峰值的均值和標(biāo)準(zhǔn)差沒有任何聯(lián)系，是呈散亂分布的。而氣隙偏差近似呈現(xiàn)正比關(guān)系，但與磁鋼相比分布較散亂。

圖5 50個電機設(shè)計方案存在槽開口偏差時轉(zhuǎn)矩脈動均值和標(biāo)準(zhǔn)差Fig.5Expected values and standard deviations of the torque ripple in 50 motors

圖6 50個電機設(shè)計方案存在氣隙偏差時轉(zhuǎn)矩脈動均值和標(biāo)準(zhǔn)差Fig.6Expected values and standard deviations of the torque ripple in 50 motors

2.3 同時存在定子和磁鋼噪音因子

當(dāng)電機同時存在磁鋼和定子噪聲因子時，同樣在可行設(shè)計域內(nèi)隨機組合50個電機設(shè)計方案，電機的各種噪聲因子均在所給定的公差帶內(nèi)變動，服從正態(tài)分布規(guī)律。每個電機設(shè)計方案進行1 000臺樣本的轉(zhuǎn)矩脈動計算，得到50個電機設(shè)計方案轉(zhuǎn)矩脈動的均值和標(biāo)準(zhǔn)差。其分布可由圖7所示。

圖7 50個電機設(shè)計方案在定子和磁鋼主要噪音因子存在的條件下轉(zhuǎn)矩脈動的均值和標(biāo)準(zhǔn)差。選擇1#，2#，3#電機進一步的樣本統(tǒng)計分布規(guī)律分析Fig.7Expected value and standard deviations of the torque ripple values in 50 motors designs with main noise factors.Choosing 1#，2#，3#motor to accomplish the statistical analysis

從圖7中可以看到，50個電機設(shè)計方案的轉(zhuǎn)矩脈動均值和標(biāo)準(zhǔn)差依然呈現(xiàn)線性分布。這是由于磁鋼偏差對電機性能分散性所造成的影響比槽口和氣隙偏差要大的多，因此，磁鋼偏差的影響主導(dǎo)了性能分散性的分布。在圖7中選擇三個不同的電機設(shè)計方案。圖8顯示了三臺電機在可行設(shè)計域范圍內(nèi)的轉(zhuǎn)矩脈動分布特性，當(dāng)轉(zhuǎn)矩脈動均值變大時，其分散性也變大，如2#電機;當(dāng)轉(zhuǎn)矩脈動均值變小，分散性也變小，轉(zhuǎn)矩脈動分布，如1#電機。

圖8 1#，2#，3#電機樣本中1 000臺樣機的轉(zhuǎn)矩脈動的分布情況Fig.8Torque ripple peak values distributions of 1 000 1#，2#and 3#motors

綜上所述，可以得到如下結(jié)論，永磁同步電機在設(shè)計可行設(shè)計域內(nèi)，各種加工偏差影響下，批次電機的均值和標(biāo)準(zhǔn)差是成正比的。其分布規(guī)律服從均值和標(biāo)準(zhǔn)差呈正比的正態(tài)分布。

3 粒子群算法優(yōu)化驗證

采用多目標(biāo)優(yōu)化的方式驗證批量電機轉(zhuǎn)矩脈動的均值標(biāo)準(zhǔn)差之間的正比關(guān)系。

針對批量生產(chǎn)電機的分散性進行優(yōu)化。在優(yōu)化過程中，將電機的輸出平均轉(zhuǎn)矩（Tav）的倒數(shù)作為一個目標(biāo)，將1 000臺的電機樣本轉(zhuǎn)矩脈動的均值（Me1000）或標(biāo)準(zhǔn)差（std1000）作為另一個目標(biāo)，分別開展兩目標(biāo)優(yōu)化。優(yōu)化目標(biāo)分別為

如果兩個優(yōu)化過程所得到的最優(yōu)解和最劣解是一致的，則證明了使用均值或者使用標(biāo)準(zhǔn)差作為電機優(yōu)化的目標(biāo)，兩者是等效的，兩個目標(biāo)之間可進行相互替代，從而驗證了兩個目標(biāo)的線性關(guān)系。

優(yōu)化采用粒子群（PSO）優(yōu)化算法，其優(yōu)化變量為［boa，Rm，Wpm，A－p，leq，αp］;電機的可行設(shè)計域如表4所示。

采用含精英歸檔策略的PSO算法進行優(yōu)化。優(yōu)化中粒子群個數(shù)為50，精英集合為40個，進化代數(shù)為20。C1=C2=0.8，其流程如圖9所示。

圖9 PSO優(yōu)化算法流程圖Fig.9Flow chart for the PSO optimization

圖10顯示了兩個優(yōu)化過程的最終結(jié)果。

圖10 兩個不同優(yōu)化目標(biāo)的優(yōu)化結(jié)果Fig.10Optimization results for two different objectives

表5顯示了兩個優(yōu)化結(jié)果得到的最穩(wěn)健解和最不穩(wěn)健解。從表中可以看到，最穩(wěn)健解完全一致，最不穩(wěn)健解僅在槽口變量上有0.3 mm的差異。從兩個優(yōu)化過程得到的優(yōu)化結(jié)果中可以看到，采用均值或標(biāo)準(zhǔn)差進行的兩目標(biāo)優(yōu)化，優(yōu)化結(jié)果是一致的，因此兩個目標(biāo)在優(yōu)化過程中是可以等效替換的，從統(tǒng)計角度分析，均值和標(biāo)準(zhǔn)差是成正比關(guān)系的。

表5 兩目標(biāo)優(yōu)化結(jié)果Table 5The results of the two different optimization

從表5中可以看到通過進行永磁同步電機的穩(wěn)健設(shè)計，電機產(chǎn)品的穩(wěn)健性顯著提高，從而使得電機產(chǎn)品轉(zhuǎn)矩脈動峰值的均值和標(biāo)準(zhǔn)差均有大幅的下降，這就保證了電機產(chǎn)品兼?zhèn)涞娃D(zhuǎn)矩脈動及低性能分散性。

4 結(jié)論

1）通過噪聲因子敏感性分析得到對轉(zhuǎn)矩脈動較為敏感的噪聲因子為磁鋼的尺寸和材料屬性偏差，以及氣隙和槽口的尺寸偏差。

2）當(dāng)存在電機加工制造偏差時，批量生產(chǎn)電機轉(zhuǎn)矩脈動峰值的均值和標(biāo)準(zhǔn)差近似呈正比關(guān)系，其分布規(guī)律服從均值和標(biāo)準(zhǔn)差呈正比的正態(tài)分布。

3）通過進行穩(wěn)健設(shè)計，提高了電機產(chǎn)品的穩(wěn)健性，抑制了產(chǎn)品整體的轉(zhuǎn)矩脈動峰值及產(chǎn)品的性能分散性。

［1］齊蓉，林輝，陳明.無刷直流電機換向轉(zhuǎn)矩脈動分析與抑制［J］.電機與控制學(xué)報，2006，10（3）:287－290.

QI Rong，LIN Hui，CHEN Ming.Analysis and attenuation on commutation torque ripples of brushless DC motor［J］.Electric Machines and Control，2006，10（3）:287－290.

［2］王淑紅，熊光煜.無刷直流電機換相轉(zhuǎn)矩脈動減小及動態(tài)仿真［J］.電機與控制學(xué)報，2008，12（2）:169－173.

WANG Shuhong，XIONG Guangyu.Commutation torque ripple minimization and dynamic simulation of brushless DC motor［J］.Electric Machines and Control，2008，12（2）:169－173.

［3］李耀華，劉衛(wèi)國.永磁同步電機直接轉(zhuǎn)矩控制不合理轉(zhuǎn)矩脈動［J］.電機與控制學(xué)報，2007，11（2）:148－152.

LI Yaohua，LIU Weiguo.Unreasonable torque ripple of direct torque control in permanent magnet synchronous motor［J］.Electric Machines and Control，2007，11（2）:148－152.

［4］郭宏，錢浩.永磁同步電機低轉(zhuǎn)矩脈動的穩(wěn)健設(shè)計［J］.中國電機工程學(xué)報，2012，32（24）:88－95.

GUO Hong，QIAN Hao.Robust design for reducing torque ripple in permanent magnet synchronous motor［J］.Proceedings of the CSEE，2012，32（24）:88－95.

［5］HEINS G，BROWN T，THIELE M.Statistical analysis of the effect of magnet placement on cogging torque in fractional pitch permanent magnet motors［J］.IEEE Transactions on Magnetics，2011，47（8）:2142－2148.

［6］GASPARIN L，CERNIGOJ A，MARKIC S，et al.Additional cogging torque components in permanent-magnet motors due to manufacturing imperfections［J］.IEEE Transactions on Magnetics，2009，45（3）:1210－1213.

［7］SIMON－SEMPERE V，BURGOS－PAYAN M，CERQUIDESBUENO J R.Influence of manufacturing tolerances on the electromotive force in permanent-magnet motors［J］.IEEE Transactions on Magnetics，2013，49（11）:5522－5532.

［8］PONOMAREV P，PETROV I，PYRHONEN J.Influence of travelling current linkage harmonics on inductance variation，torque ripple and sensorless capability of tooth-coil permanent-magnet synchronous machines［J］.IEEE Transactions on Magnetics，2014，50（1）:1－8.

［9］ZHU Z Q，WU L J，XIA Z P.An accurate subdomain model for magnetic field computation in slotted surface-mounted permanentmagnet machines［J］.IEEE Transactions on Magnetics，2010，46 （4）:1100－1115.

［10］WU L J，ZHU Z Q，STATON D A，et al.Comparison of analytical models of cogging torque in surface-mounted pm machines［J］.IEEE Transactions on Industrial Electronics，2012，59（6）: 2414－2425.

（編輯:張詩閣）

Distribution characteristics analysis of torque ripple in permanent magnet synchronous motors considering the noise factors

WU Zhi-yong，GUO Hong，QIAN Hao
（School of Automation Science and Electrical Engineering，Beijing University of Aeronautics and Astronautics，Beijing 100191，China）

In mass productions of permanent magnet synchronous motors（PMSM），the torque ripple dispersion is caused by the sizes errors and the material properties variations，which can also be named as noise factors.In order to analyze the distribution of torque ripple in mass productions of PMSM，the sensitive analysis was operated and the most sensitive noise factors which greatly affect the torque ripple were obtained that are the magnets misplacement properties error，the slot opening error and the air-gap size error.Then the distribution of the torque ripple in mass productions was estimated by the firstly the theoretical analysis and by the analytical calculation.The relationship of the mean value and standard derivation of the torque ripple in mass production are linear which is then proved by the optimum design results.

permanent magnets;synchronous motors;sensitivity analysis;statistical analysis;torque

10.15938/j.emc.2015.01.002

TM 34

A

1007－449X（2015）01－0008－08

2014－01－25

國家自然科學(xué)基金（51277007）

武志勇（1985—），男，博士，研究方向為特種電機設(shè)計及驅(qū)動控制器設(shè)計;

郭宏（1969—），男，博士，教授，博士生導(dǎo)師，研究方向為特種電機及其控制;電機多物理場綜合設(shè)計與優(yōu)化等;

錢浩（1985—），男，博士，研究方向為高精度伺服電機設(shè)計與驅(qū)動控制。

武志勇