彭志遠，杜長虹，陳 健，馬永泉，周安健，任 勇

(1.重慶長安新能源汽車科技有限公司，重慶 401120；2.長安福特馬自達發(fā)動機公司，南京 211100)

在永磁同步電機的控制中，獲取轉子溫度具有非常重要的意義，主要表現在以下3方面[1-3]：① 電機轉子溫度過高極有可能造成永磁磁鋼不可逆的退磁，造成電機系統(tǒng)運行安全風險；② 電機轉子磁鋼退磁會進一步影響電驅動系統(tǒng)扭矩的計算精度，造成整車性能控制的偏差；③ 電機轉子溫度的精確估算可以有效增加電機峰值轉矩持續(xù)時間，在提升電機潛能的同時，減少電機定子磁材料的使用。因此，電機轉子溫度精確估算不僅能確保電機系統(tǒng)安全運行，還可以提升電機系統(tǒng)驅動性能以及降低其制造成本。

目前的研究中，電機轉子溫度估算方法主要分為“經驗公式間接計算”“熱網絡分析計算”“反電勢法計算”三大類[4-9]，第1類方法是間接測量與電機轉子溫度相關的變量，通過經驗公式計算獲得；第2類方法是將電機細分成單元體，建立單元體之間的熱阻并形成熱網絡模型，通過在熱網絡模型中加入補償單元，結合有限元法對計算結果進行驗證；第3類方法是測量電機反電動勢，計算得到電機剩余磁通密度，通過查詢剩余磁通密度與轉子溫度對應關系獲得轉子實際溫度。綜合分析目前現有的研究方法：首先，沒有考慮整車系統(tǒng)下電后，自然冷卻條件下的電機轉子溫度變化情況，無法做到系統(tǒng)再次上電時對轉子初始溫度進行賦值計算；其次，只考慮常溫條件下的轉子溫度特性，沒有考慮環(huán)境溫度對轉子溫度特性的影響，造成算法適應性差且精度有限；再次，利用反電勢法測量轉子溫度時，需要將電機電流卸載，不適用于整車實際工況下的估算。因此，本文將綜合考慮影響電機轉子溫度的熱節(jié)點，通過試驗方法獲得電機轉子溫度特性的規(guī)律應不同環(huán)境溫度、不同運行工況下的轉子溫度算法，以提高電機系統(tǒng)性能與運行安全，實現低制造成本的目的。

1 轉子溫度熱模型分析

為實現電機轉子溫度的精確估算，需要從影響轉子溫度的根源入手，重點分析電機實際運行工況下各主要溫度節(jié)點生熱及熱傳導機理。圖1給出了電機系統(tǒng)剖面結構示意圖，現分別對電機在帶載運行、空載運行、停機工況下，電機轉子與各主要熱節(jié)點的熱交換進行定性分析，建立各工況下電機轉子的等效熱模型。

1.冷卻水道；2.定子繞組；3.定子；4.轉子；5.轉子磁石鋼；6.軸承；7.轉軸

1.1 電機運行狀態(tài)熱分析

電機帶載運行過程中，必然要產生相應的損耗，這些損耗轉換成熱與周圍環(huán)境進行熱交換，最終達到熱平衡。此過程中，電機轉子溫升主要受到銅損耗、鐵損耗、機械損耗以及冷卻損耗的影響[10]。

1.1.1電機銅損耗計算

銅損耗占總損耗的絕大部分，是主要的熱源。假定電流在定子繞組截面上均勻分布，且忽略溫度對定子繞組相電阻的影響，銅損耗可以由下面公式計算[11]：

Pcu=nI2R

(1)

式中：n為電機相數；I為相電流；R為相電阻。

1.1.2電機鐵損耗計算

鐵損耗包括磁滯損耗和渦流損耗。磁滯損耗由定子繞組中交變電流引起交變磁場變化而產生的，渦流損耗由鐵心中的磁場變化，產生感應電流而引起。鐵損耗[12]：

(2)

式中：kh為磁滯損耗系數；ke為渦流損耗系數；kexc為附加損耗系數；f為電樞磁場交變頻率；Bm為定子鐵心磁通密度幅值。

1.1.3機械損耗計算

機械損耗包括軸承轉動引起的摩擦損耗和通風引起的損耗。機械損耗[13]：

(3)

式中：kc為表面粗糙系數；Cf為摩擦因數；ρair為空氣密度；ωm為電機角速度；l為轉子長度；r為轉子半徑。

1.1.4冷卻損耗計算

冷卻損耗包括冷卻回路帶走的定子功率損耗、定子及其繞組對周圍環(huán)境的耗散功率損耗。冷卻液吸收的定子功率損耗[14]：

Pw=ρwCwAwv(Tin-Tout)

(4)

式中：ρw為冷卻液密度；Cw為冷卻液比熱容；Aw為冷卻管路截面積；v為冷卻液流動速度；Tin、Tout分別為進、出水溫度。

定子及其繞組對周圍環(huán)境的耗散功率損耗：

(5)

式中：δ為對流換熱系數；As為定子及繞組表面對流換熱面積；Ts為定子及繞組表面溫度；T0為環(huán)境溫度；t2、t1分別為單位時間內的起止時刻點。

對流換熱系數采用如下經驗公式計算[15]：

δ=9.73+14V0.62

(6)

式中V為散熱表面空氣流通速度。

1.1.5轉子溫度模型建立

電機帶載運行單位時間內消耗的能量應滿足能量守恒定律，因此有：

Pcu+PFe+Pfr=Pw+Pair+Pr

(7)

式中：Pr為單位時間內轉子發(fā)熱損耗；Pair為單位時間內系統(tǒng)總耗散功率。

轉子熱損耗帶來其溫度變化，那么電機帶載運行工況下，任何時刻轉子溫度可由下式計算：

(8)

式中：Cr為轉子材質比熱容；Mr為轉子質量；Tr1、Tr2分別為單位時間內起止時刻點的轉子溫度。

1.2 電機停機狀態(tài)熱分析

電機系統(tǒng)處于停機狀態(tài)時，轉子僅與周圍環(huán)境存在自然冷卻條件下的熱交換并達到最終熱平衡。此狀態(tài)下的熱模型可以等效為電機轉子作為一個熱源向環(huán)境釋放熱量，直到溫度下降到與環(huán)境溫度保持一致時，熱交換完成。

轉子及磁鋼對周圍環(huán)境的耗散功率損耗由下式計算：

(9)

式中：Ar為轉子表面對流換熱面積；Tr1與式(8)中變量保持一致；δ計算見式(6)。

2 轉子溫度算法構建

電動車運行環(huán)境與工況十分復雜，為精確估算電機在不同環(huán)境溫度與不同載荷工況下轉子溫度，構建了電機轉子溫度算法，如圖2所示。該算法首先判斷系統(tǒng)是否上電，然后根據電機實際轉速、轉矩判斷電機實時工作狀態(tài)，針對電機運行狀態(tài)，采用電機系統(tǒng)損耗功率守恒的原則對轉子溫度熱模型進行簡化，最后通過修正模塊對模型計算的轉子溫度進行校準，在系統(tǒng)下電時對修正后的轉子溫度進行非遺失性存儲，以便系統(tǒng)再次上電后對轉子初始溫度的計算。

2.1 停機狀態(tài)溫度計算

整車上電后，系統(tǒng)首先進入停機狀態(tài)(狀態(tài)標志SF=2)，此時需要對電機轉子初始溫度進行計算。為減小控制系統(tǒng)運算量，通過臺架試驗分別記錄環(huán)境溫度點(0～60 ℃，相鄰間隔為10 ℃)自然冷卻條件下，電機轉子由最高穩(wěn)定溫度下降到環(huán)境溫度的變化規(guī)律(如圖3所示)，建立相關數值模型，并簡化停機狀態(tài)下的轉子熱模型。

圖2 電機轉子溫度算法流程

圖3 自然冷卻條件下轉子溫度變化

由圖3可以看出：隨著環(huán)境溫度的升高，電機轉子由初始溫度自然冷卻到環(huán)境溫度的熱平衡時間逐漸增大。由此可見，環(huán)境溫度對電機轉子的耗散功率影響較大，為有效簡化停機狀態(tài)下的轉子熱模型，通過以下步驟實現系統(tǒng)上電后的轉子初始溫度計算：

1)讀取電機控制系統(tǒng)記錄的上次下電時存儲的轉子溫度Tr_pre、電池控制系統(tǒng)發(fā)送的停機時長tstop以及整車控制器發(fā)送的環(huán)境溫度T0；

2)查詢對應環(huán)境溫度自然冷卻條件下的轉子溫度數值曲線(如環(huán)境溫度30 ℃)，根據系統(tǒng)上次下電時記錄的轉子溫度Tr_pre，找到該溫度對應時間點t0并加上停機時長tstop，查詢時間點(t0+tstop)時間點的溫度值，即獲得系統(tǒng)本次上電時的轉子初始溫度Tr_init(若系統(tǒng)停機時長tstop大于等于自然冷卻到環(huán)境溫度的熱平衡時長tb-t0，則視轉子初始溫度為環(huán)境溫度)。

3)對于其他環(huán)境溫度點下的轉子初始溫度，可以通過式(10)插值擬合計算，其中擬合系數由試驗標定確定，見表1所示。

Tr_init=T1_2=ξT1+(1-ξ)T2

(10)

式中：T1、T2分別為相鄰環(huán)境溫度點下獲得的轉子初始溫度；T1_2為其他環(huán)境溫度點下插值擬合計算出的轉子初始溫度；ξ為擬合系數，見表1。

表1 轉子初始溫度擬合系數

2.2 運行狀態(tài)溫度計算

電機系統(tǒng)運行過程中產生的銅損耗、鐵損耗、機械損耗、冷卻損耗將對定子與轉子溫度產生影響。由建立的系統(tǒng)熱模型可知，電機運行產生的銅損耗、鐵損耗與機械損耗主要使定子升溫，并作為一個熱源對轉子進行加熱，根據采樣時間內電機系統(tǒng)損耗功率守恒的原則對運行狀態(tài)轉子溫度模型進行簡化，如下所示：

(11)

式中：Cs為定子材料比熱容；Ms為定子質量；Ts1、Ts2分別為前后采樣時間點的定子溫度；Pw、Ps_air、Pr_air分別由式(4)(5)及式(9)計算獲得。

系統(tǒng)進入運行狀態(tài)時(狀態(tài)標志SF=1)，首先調用停機狀態(tài)計算的轉子初始溫度Tr_init，然后將其代入式(8)與式(9)中的變量Tr1，并結合式(11)中實時計算的Pr求出單位時間內終止時刻的轉子溫度Tr2，最后將當前終止時刻獲得的Tr2不斷循環(huán)迭代到式(8)中變量Tr1，計算出實時運行狀態(tài)下的轉子溫度。

2.3 轉子溫度修正計算

利用式(8)不斷迭代計算的實時轉子溫度在一定累計時長后會產生相應誤差，為確保轉子溫度的精確計算，通過試驗數據擬合出電機反電勢與轉子溫度對應關系(如圖4所示)，采用以下步驟與修正條件實現電機轉子溫度修正：

1)當電機轉速在1 000～5 000 r/min且轉矩為0 N·m狀態(tài)時，通過三相合成電壓Us推算電機反電勢Uphi。

圖4 電機反電勢與轉子溫度對應關系

3 試驗驗證

為了驗證本文所建立的基于系統(tǒng)熱節(jié)點等效模型的轉子溫度實時控制算法，利用AVL測試系統(tǒng)搭建電機轉子溫度精度驗證臺架，如圖5所示。該臺架系統(tǒng)由電力測功機、電池模擬器、高低溫箱、冷卻系統(tǒng)、電機總成及控制器、電參數測試儀、可調低壓電源及相關傳感器組成。

圖5 電機轉子溫度試驗臺架實物圖

由于環(huán)境溫度、運行工況對電機轉子溫度影響較大，本文在不同環(huán)境溫度、不同運行工況條件下，將模型計算值與實測值進行對比，通過不斷標定修正系統(tǒng)冷卻損耗數值模型來優(yōu)化算法精度。圖6給出環(huán)境溫度30 ℃且電機處于停機冷卻條件下模型計算值與實測值隨冷卻時間的變化關系，由于系統(tǒng)無運行功率，電機定子作為唯一熱源對外界耗散熱量，模型計算值與實測值隨冷卻時間的增加而階梯減小，最大誤差在±3 ℃以內。圖7、圖8分別給出環(huán)境溫度30 ℃與60 ℃且變功率條件下模型計算值與實測值隨電機運行時間的變化關系，當電機運行功率由高向低切換過程中，由于定子發(fā)熱功耗減小使得其傳遞給轉子的加熱功率低于冷卻損耗功率，轉子溫度逐漸下降，反之轉子溫度逐漸上升。整個變功率運行工況過程中，模型計算值與實測值動態(tài)變化跟隨性較好，最大誤差在±10 ℃以內。

圖6 停機冷卻轉子溫度估算驗證結果

圖7 環(huán)境溫度30 ℃轉子溫度估算驗證結果

圖8 環(huán)境溫度60 ℃轉子溫度估算驗證結果

4 結論

1)通過對電機系統(tǒng)各個熱節(jié)點在不同工作狀態(tài)下的發(fā)熱機理進行分析，建立了轉子溫度熱模型。

2)根據建立的轉子溫度熱模型，結合電機系統(tǒng)實際運行工況，提出了基于系統(tǒng)熱節(jié)點等效模型的轉子溫度實時控制算法，對電機不同工況下的轉子溫度計算模型進行了簡化，并利用測試數據構建起算法中相應的數值模型。

3)搭建電機轉子溫度控制算法驗證臺架，在不同環(huán)境溫度條件下，采用變功率運行工況對電機轉子溫度模型計算值與實測值進行比較，通過反復優(yōu)化系統(tǒng)相關熱節(jié)點功率損耗數值模型，實現停機冷卻條件下模型計算值與實測值最大誤差在±3 ℃以內，變功率動態(tài)工況下模型計算值與實測值最大誤差在±10 ℃以內，滿足電機轉子溫度估算的使用需求。