王曉遠,梁慧靜,李天元,趙曉曉
(天津大學電氣自動化與信息工程學院,天津 300072)
集成起動/發(fā)電機ISG(integrated starter/generator)系統(tǒng)因其快速啟停、質(zhì)量輕、成本低、污染小及能量回收再利用等優(yōu)點,廣泛應(yīng)用于電動汽車、航空航天、工業(yè)機床等領(lǐng)域[1-2]。內(nèi)置式永磁同步電機IPMSM(interior permanent magnet synchronous motor)因其效率高、功率密度高、體積小等優(yōu)點廣泛應(yīng)用于ISG電機系統(tǒng)[3]。
永磁起動/發(fā)電機PMISG(permanent magnet integrated starter/generator)在實際運行時,氣隙中含有不同階次諧波與旋轉(zhuǎn)轉(zhuǎn)子發(fā)生相對運動,在永磁體內(nèi)產(chǎn)生大量渦流損耗[4]。與此同時,混合電動汽車空間有限,PMISG 軸向長度要求盡可能短,為縮短線圈端部長度常采用分數(shù)槽集中繞組形式,相比于整數(shù)槽繞組會產(chǎn)生更多的渦流損耗[5]。渦流損耗會引起永磁體溫度升高,而車內(nèi)空間狹小散熱環(huán)境不良,過高的溫度會造成永磁體退磁的惡劣影響[6]。電動汽車在起動過程中PMISG 工作在電動狀態(tài),帶動發(fā)動機到達臨界轉(zhuǎn)速,此過程中采用了弱磁控制方法獲取大轉(zhuǎn)速,期間磁通會急劇變化,永磁體內(nèi)渦流損耗也因此增大。在弱磁升速過程中,電機磁通飽和程度會發(fā)生變化,從而對dq軸電感及磁鏈等參數(shù)產(chǎn)生影響,同時對渦流損耗也會產(chǎn)生影響。電動狀態(tài)的PMISG 在獲得最大轉(zhuǎn)矩時磁路會達到飽和狀態(tài),使得PMISG 成為一個非線性、多耦合的系統(tǒng),永磁體的渦流損耗分析更加復雜。
永磁體渦流損耗直接影響電機效率和永磁體溫度,為了提高效率、降低永磁體溫度和避免退磁,針對各種類型電機的永磁體渦流損耗研究一直是眾多學者研究的熱點。Yamazaki等[7]研究了在脈寬調(diào)制PWM(pulse width modulation)控制模式下永磁同步電機在二維和三維模型中氣隙磁場的諧波情況,以及各次諧波在轉(zhuǎn)子永磁體所產(chǎn)生的渦流損耗,并進行了實驗驗證。Chen 等[8]探究了表貼式永磁同步電機中永磁體的厚度、充磁方式等因素對渦流損耗的影響。文獻[9]分析了電動汽車用PMISG不同弱磁角與不同轉(zhuǎn)速情況下的永磁體渦流損耗,但是沒有考慮到弱磁條件下轉(zhuǎn)速和超前角的相互聯(lián)系,只進行了獨立分析,沒有對磁場的變化做進一步討論。
綜上所述,PMISG 在電動汽車領(lǐng)域的應(yīng)用一直是研究熱點,其在弱磁過程中的永磁體渦流損耗卻少有研究。在電機實際運行的過程中,弱磁情況下電機磁場分布會發(fā)生變化,導致電感也隨之變化。此外,弱磁運行下磁場分布的變化會造成感應(yīng)在永磁體上的渦流發(fā)生變化,進而影響永磁體渦流損耗。因此,本文針對兩種電機在弱磁運行下的永磁體渦流損耗與磁場進行分析研究,選取兩臺不同轉(zhuǎn)子的35 kW內(nèi)置式PMISG為研究對象,在定子結(jié)構(gòu)完全相同、轉(zhuǎn)子磁鋼用量相同的條件下,比較了V 型和雙V 型兩種轉(zhuǎn)子磁極結(jié)構(gòu)的直交軸電感和永磁體渦流損耗。利用有限元法,選取額定電流和最大電流兩種最具代表性的工況時的電感、永磁體渦流損耗與弱磁程度的關(guān)系進行了具體分析。最后制作樣機進行實驗并對仿真分析結(jié)果進行了驗證。
PMISG 選用了12 極18 槽的極槽配合,每對極下的結(jié)構(gòu)和場量都相同,因此可選擇1對極結(jié)構(gòu)作為二維求解域模型,簡化計算量,模型如圖1 所示。轉(zhuǎn)子磁極結(jié)構(gòu)分別為V 型和雙V 型。在定子完全相同的前提下,兩種轉(zhuǎn)子磁極結(jié)構(gòu)具有相同的磁鋼用量及分塊數(shù)量。

圖1 PMISG 的1/6 轉(zhuǎn)子結(jié)構(gòu)模型Fig.1 Structural model of one sixth of PMISG rotor
混合動力汽車用PMISG的部分基本參數(shù)如表1所示。

表1 PMISG 樣機參數(shù)Tab.1 Parameters of PMISG prototype
對于輕型混合電動汽車,PMISG既作為電動機又作為發(fā)電機。PMISG 的發(fā)電工況已有很多學者進行研究,且在此工況的弱磁程度由所帶負載決定,研究磁場變化與渦流損耗之間的關(guān)系十分重要。因此,本文對PMISG電動工況時弱磁條件下的渦流損耗進行研究。
IPMSM 的dq軸磁路不對稱,dq軸電感相差較大,為弱磁控制帶來了很大優(yōu)勢[10]。IPMSM 運行過程電流、電壓軌跡如圖2所示。

圖2 IPMSM 運行過程電流、電壓軌跡Fig.2 Trajectories of current and voltage of IPMSM under operation
IPMSM弱磁工況時運行在圖2中AC段。在AC段,dq軸電流滿足

式中:irms為電流有效值;imax為最大受限電流;id、iq分別為dq軸電流。
dq軸電壓滿足

式中:urms為電壓有效值;ω為轉(zhuǎn)子旋轉(zhuǎn)電角速度;umax為逆變器由直流母線電壓得到的最大受限電壓;ψf為永磁體磁鏈幅值;Ld、Lq分別為dq軸電感。
直軸電流id為

忽略定子電阻的影響,可推導出弱磁控制下的轉(zhuǎn)速方程,即

式中:n為轉(zhuǎn)子轉(zhuǎn)速;ρ為凸極率,
當電機運行在AC段弱磁控制工況下,轉(zhuǎn)速已經(jīng)超過基速,逆變器的電壓極限輸出為主要限制因素。此時,AC段dq軸電流與轉(zhuǎn)速之間為一一對應(yīng)的關(guān)系,即超前角與轉(zhuǎn)速之間一一對應(yīng)。
對永磁體的渦流損耗進行計算,做如下假設(shè):
(1)諧波磁通密度在空間按照正弦分布,忽略永磁體表面對諧波幅值削弱作用;
(2)永磁材料的磁導率μ為常數(shù);
(3)磁場穿過永磁體是均勻分布的。
圖3 為永磁體渦流損耗計算模型。其中,磁場方向與永磁體充磁方向相同,且在垂直平面Oxy內(nèi)形成渦流,如虛線部分所示;a為永磁體軸向長度;b為寬度;dx、dy分別為渦流在x、y軸方向的長度微分;β為永磁體寬度與渦流回路直角的連線的夾角。

圖3 永磁體渦流損耗計算模型Fig.3 Calculation model of permanent magnet eddy current loss
渦流環(huán)中y方向的有效長度用x表示,即

式中:y0為y方向上的初始長度;ky為正切系數(shù);β為永磁體寬度與渦流回路直角的連線的夾角。
1個渦流環(huán)的磁通φ(x,t)可表示為

式中:B(t)為渦流環(huán)內(nèi)磁通密度;y(x)為渦流環(huán)內(nèi)y方向有效長度。渦流環(huán)內(nèi)產(chǎn)生的感生電動勢可表示為

式中:R(x)為渦流環(huán)內(nèi)阻抗;τ為永磁體電阻率;h為永磁體厚度;ke為端部效應(yīng)有關(guān)的系數(shù),其初始值為1。
永磁體上1個渦流環(huán)的功率損耗Peddy1可表示為

在x方向做積分可得整塊永磁體的渦流損耗Peddy為

由式(11)可以看出,永磁體渦流損耗與其體積、厚度和磁密有關(guān)。在電機磁鋼一定時,永磁體渦流損耗由磁通的交變頻率與諧波含量決定。

式中:B(t,θ)為永磁體內(nèi)磁密;Bm為永磁體內(nèi)m次諧波磁密幅值;ωm為m次諧波電角頻率;θm為m次諧波相角。
本文諧波主要指由于電機結(jié)構(gòu)(例如定子齒槽、轉(zhuǎn)子結(jié)構(gòu)、磁鋼排布等)產(chǎn)生的空間諧波。渦流損耗不僅與弱磁角有關(guān),還與Ld、Lq相關(guān)。
dq軸電感參數(shù)在電機控制中十分關(guān)鍵。PMISG 在較大負載時,磁路飽和、交叉耦合現(xiàn)象更為嚴重。根據(jù)電感定義,dq軸電樞反應(yīng)電感為

式中:Mdq、Mqd為dq軸互感參數(shù);ψd、ψq為dq軸磁鏈。
采用V 型和雙V 型兩種轉(zhuǎn)子磁極結(jié)構(gòu)的PMISG為樣機,對PMISG空載與實際弱磁過程進行仿真分析。
V型和雙V型兩種轉(zhuǎn)子結(jié)構(gòu)電機的空載反電勢如圖4 所示。兩種轉(zhuǎn)子結(jié)構(gòu)空載相反電勢的有效值分別為74.28 V 和68.28 V,諧波畸變率分別為0.32和0.25。雙V型轉(zhuǎn)子結(jié)構(gòu)電機的空載反電動勢波形更接近正弦,諧波含量更少。

圖4 空載反電勢相電壓及其傅里葉分解Fig.4 Phase voltage of no-load back EMF and its Fourier decomposition
兩種轉(zhuǎn)子磁極結(jié)構(gòu)永磁同步電機的空載磁場分布如圖5 所示,其中包含了電機磁力線分布、各部分剖分精度與磁密云圖,證明了磁路設(shè)計與各部分剖分精度的合理性。

圖5 兩種轉(zhuǎn)子磁極結(jié)構(gòu)空載的磁場分布Fig.5 Magnetic field distribution of two types of rotor pole structure under no-load
選取額定電流和最大電流兩種工況進行研究分析,因為這兩種工況最具有代表性。對于恒轉(zhuǎn)矩階段后的弱磁工況,由于此時已到達最大受限電壓,為增大轉(zhuǎn)速需要增大弱磁角度。此后,弱磁角與轉(zhuǎn)速一一對應(yīng)。每隔500 r/min取1個點,將弱磁角設(shè)置為變量,尋找其滿足最大電壓umax的弱磁角,仿真分析轉(zhuǎn)速大于額定轉(zhuǎn)速的運行工況。
圖6 為在額定電流情況下兩種轉(zhuǎn)子結(jié)構(gòu)電機不同弱磁條件下的dq軸電感??梢钥闯?,雙V型電機的交直軸電感差值大,弱磁性能更好。在電流一定的前提下,隨著弱磁階段升速,直軸去磁電流增大,直軸磁路磁密幾乎保持在線性區(qū)域,所采用定子鐵心材料磁化特性曲線的磁導率先增大后減小,因此直軸電感會隨著弱磁程度的增大而略有降低;同時交軸電流減小,交軸磁路由飽和狀態(tài)變?yōu)椴伙柡?,交直軸磁路有重合部分,直軸磁通降低也有利于交軸磁路磁密處于線性區(qū)域,因此交軸電感在弱磁區(qū)域增大,這種影響在交軸磁路較飽和時更顯著。雙V 型PMISG 直軸磁路上的永磁體數(shù)比V 型的多1倍,因此直軸電感更小。

圖6 額定電流時不同弱磁條件下的直交軸電感Fig.6 Direct and quadrature axis inductances under different field weakening conditions at rated current
圖7 為最大電流情況下兩種轉(zhuǎn)子結(jié)構(gòu)電機在不同弱磁條件下的dq軸電感。相比于額定電流工況,由于磁密由飽和變化到不飽和,磁通變化較大,電感變化的速率變大。最大電流工況的電感初始時要小很多,因為此時飽和程度深、電感小。隨著磁密由飽和降到不飽和,交軸電感的變化速率也出現(xiàn)了明顯增加,這是由于去磁電流的增大,削弱了交軸磁場。

圖7 最大電流時不同弱磁條件下的直交軸電感Fig.7 Direct and quadrature axis inductances under different field weakening conditions at maximum current
圖8和圖9分別為額定電流工況與最大電流工況下兩種轉(zhuǎn)子結(jié)構(gòu)在不同弱磁條件下的渦流損耗。

圖8 額定電流時不同弱磁條件下的永磁體渦流損耗Fig.8 Permanent magnet eddy current loss under different field weakening conditions at rated current

圖9 最大電流時不同弱磁條件下的永磁體渦流損耗Fig.9 Permanent magnet eddy current loss under different field weakening conditions at maximum current
從圖8 和圖9 可以看出,雙V 型轉(zhuǎn)子的渦流損耗明顯小于V型轉(zhuǎn)子,這是由于雙V型轉(zhuǎn)子結(jié)構(gòu)諧波含量更少。弱磁控制對交直軸磁場分布產(chǎn)生了影響,相應(yīng)地永磁體內(nèi)感應(yīng)的渦流也發(fā)生了變化。飽和時磁密中諧波含量較多,隨著交軸磁通飽和程度逐漸降低,永磁體渦流損耗變化速率會有所降低。弱磁過程對最大電流工況的雙V 型結(jié)構(gòu)電機永磁體渦流損耗的影響更加顯著,在恒轉(zhuǎn)矩階段的拐點后更加明顯,這是因為在最大電流工況下磁路飽和,而在額定電流工況下磁路是不飽和的。由于飽和強度深,雙V 型弱磁作用更明顯,因此雙V 型永磁體渦流損耗的增速與轉(zhuǎn)速不呈2次函數(shù)關(guān)系,而與磁場、電感、頻率緊密相關(guān),呈現(xiàn)出緩慢的1次增長趨勢。雙V 型PMISG 的交軸磁路磁通飽和程度較低,因此在弱磁階段其渦流損耗的變化速率和恒轉(zhuǎn)矩階段接近。
為了驗證本文仿真分析的有效性,根據(jù)電機設(shè)計參數(shù)對本文所研究的V 型磁極結(jié)構(gòu)轉(zhuǎn)子PMISG進行了空載實驗。圖10為電機轉(zhuǎn)子結(jié)構(gòu)與樣機裝配圖。圖11 為樣機實驗平臺。PMISG 與測功機相連,使用示波器記錄空載反電勢線電壓。

圖10 電機轉(zhuǎn)子結(jié)構(gòu)與電機裝配Fig.10 Motor rotor structure and motor assembly

圖11 樣機實驗平臺Fig.11 Prototype test platform
PMISG永磁體渦流損耗不易直接測量,需要通過間接法計算。永磁材料在不同溫度下會展現(xiàn)出不同的磁性能,根據(jù)其在空載狀態(tài)下定子繞組中感應(yīng)出的反電動勢的變化程度,便可推出轉(zhuǎn)子永磁體的溫升情況,進而根據(jù)熱力學原理計算永磁體渦流損耗。
在溫度為T下定子反電動勢的表達式[11]為

式中:eT為溫度T下的定子反電勢;f為電源頻率;Kdq為繞組系數(shù);φ為通過永磁體的磁通;Kφ為波形系數(shù);N為繞組匝數(shù);Br,T為溫度T下的剩磁;Am為每極下永磁體面積;bmo為永磁體工作點;σ0為空載漏磁系數(shù)。其中,永磁體在溫度為T1和T2時的剩磁關(guān)系為

式中,α為永磁體溫度系數(shù)。
由式(14)和式(15)可得,溫升ΔT與反電勢關(guān)系為

永磁體的渦流損耗表達式為

式中:Pm為永磁體渦流損耗;MPM為永磁體質(zhì)量;CPM為永磁體比熱容;Mrotor為轉(zhuǎn)子鐵心質(zhì)量;Crotor為轉(zhuǎn)子鐵心比熱容;Mshaft為轉(zhuǎn)子轉(zhuǎn)軸質(zhì)量;Cshaft為轉(zhuǎn)子轉(zhuǎn)軸比熱容。將式(16)代入式(17)可求得永磁體渦流損耗。
空載實驗測量PMISG 在不同轉(zhuǎn)速下運行的空載線反電勢隨時間的衰減情況。其中,在額定轉(zhuǎn)速下運行35 min 空載線反電勢隨時間的衰減曲線如圖12所示。

圖12 轉(zhuǎn)速為3 500 r/min 空載反電勢下降法測量結(jié)果Fig.12 Measurement result of no-load back-EMF using drop method at rotation speed of 3 500 r/min
通過計算可得電機轉(zhuǎn)速為3 500~7 000 r/min時的永磁體渦流損耗,并將其與仿真結(jié)果進行對比,如圖13 所示。可見,實驗結(jié)果比仿真結(jié)果大9%左右,原因是有限元仿真并未考慮電機端部效應(yīng),得到的渦流回路電阻較??;電機的制作工藝難以達到高精度,采用間接測量的方法會有誤差。仿真結(jié)果與實驗結(jié)果趨勢相同,誤差可以被接受,驗證了本文分析的合理性。

圖13 渦流損耗實驗結(jié)果與仿真結(jié)果對比Fig.13 Comparison of eddy current loss between experimental and simulation results
本文以兩臺35 kW的內(nèi)置式PMISG為例,在保證其他變量一致的前提下,對V型和雙V型轉(zhuǎn)子結(jié)構(gòu)的PMISG 的空載狀態(tài)和弱磁條件下的負載狀態(tài)進行了有限元分析。分析了空載下磁路的合理性,并對兩種轉(zhuǎn)子結(jié)構(gòu)的額定電流和最大電流情況下的dq軸電感和渦流損耗做了詳細分析。通過樣機進行實驗驗證,得到以下結(jié)論。
(1)在兩臺35 kW的電機永磁體用量與分塊一致時,其空載反電勢相接近,V型轉(zhuǎn)子略大一些,這是由于雙V 型磁極漏磁較大造成的。雙V 型比V型獲得的調(diào)速范圍更廣,這是因為雙V型dq軸電感差值大,凸極率大。在永磁體渦流損耗方面,V 型轉(zhuǎn)子結(jié)構(gòu)電機的永磁體渦流損耗比雙V 型轉(zhuǎn)子大40%左右;在電感參數(shù)方面,交軸電感比直軸電感變化速率更大。
(2)對額定電流與最大電流下的弱磁控制時電機性能進行對比。在電感參數(shù)方面,電感變化率明顯增大,這是因為在最大電流運行時,隨著轉(zhuǎn)速的增加,電機磁密由飽和變?yōu)椴伙柡蜖顟B(tài),造成了電感的變化率變大。在永磁體渦流損耗方面,在最大電流運行時的永磁體渦流損耗明顯增大,且雙V型由于弱磁增長緩慢。在電機設(shè)計時,電機的弱磁能力主要由dq軸電感差值來決定,在速度增大的同時也要考慮電機的磁阻轉(zhuǎn)矩與損耗問題。電機弱磁過程會對電機渦流損耗產(chǎn)生影響,弱磁過程會對諧波磁場產(chǎn)生影響,從而降低永磁體的渦流損耗。