張 戟，王洪武，蔡浩雄

(1.同濟大學(xué)新能源汽車工程中心，上海 201804； 2.上海汽車集團股份有限公司，上海 201804)

前言

根據(jù)動力來源的不同，電動汽車可分為純電動汽車、混合動力汽車和燃料電池汽車等。它們與傳統(tǒng)汽車驅(qū)動系統(tǒng)的很大一個不同點是，集成了大功率的電力電子裝置構(gòu)成的電氣動力系統(tǒng)，包括動力電池、大功率DC/DC轉(zhuǎn)換器或DC/AC逆變器、驅(qū)動電動機及其之間相連的高壓電纜等[1]。這些電力電子系統(tǒng)工作時，將發(fā)出強烈的電磁干擾。開關(guān)器件在開通和關(guān)斷中，由于電壓和電流在短時間內(nèi)發(fā)生跳變，從而形成輻射和傳導(dǎo)電磁干擾。另一方面，連接動力驅(qū)動部件的電纜為高壓電纜，工作時有大電流(>100A)流過，電流的波動也會產(chǎn)生強烈的電磁輻射[2]。上述的零部件在工作中，供電情況、負(fù)載情況和在車輛中的布置和連接方式與傳統(tǒng)車輛零部件有較大不同。目前國內(nèi)外對這方面的研究很少，本文中針對電動汽車驅(qū)動系統(tǒng)的電磁兼容(EMC)性能進行研究具有理論和實際意義。

1 基于RMxprt的電機模型建立

應(yīng)用RMxprt模塊可分析12大類、15種常用電機[3]。交流變頻永磁同步電機轉(zhuǎn)子采用永磁體，不需要勵磁，因此功率因數(shù)大，電機具有較高的功率密度和效率[4]。目前多數(shù)驅(qū)動系統(tǒng)均采用永磁同步電機，本文中永磁同步電機設(shè)定為內(nèi)轉(zhuǎn)子；電機的極數(shù)為8；機械摩擦損耗估算為300W；風(fēng)磨損耗給定為600W；基準(zhǔn)轉(zhuǎn)速為4 000r/min；控制類型為交流AC；驅(qū)動電路類型為Y型三相。

定子參數(shù)設(shè)定：鐵芯外徑為260mm；鐵芯內(nèi)徑為160mm；鐵芯實際軸向長度為200mm；鐵芯疊壓系數(shù)為0.95；鐵芯沖片材料為DW315_50；定子槽數(shù)為24個；定子槽型設(shè)置為2號槽型——梨形槽；斜槽數(shù)為1，即斜過1個定子槽。

定子繞組參數(shù)設(shè)定：繞組層數(shù)為2層；匝間連接方式為全極式；并聯(lián)支路數(shù)為2；每槽匝數(shù)為4；線圈節(jié)距為3；線圈不并繞；不用漆絕緣；線規(guī)為0.724mm。

轉(zhuǎn)子參數(shù)設(shè)定：轉(zhuǎn)子外徑為159mm；轉(zhuǎn)子內(nèi)徑為60mm；轉(zhuǎn)子軸向長度為200mm；轉(zhuǎn)子沖片材料為DW315_50；鐵芯疊壓系數(shù)為0.95；電極類型為5。

設(shè)置轉(zhuǎn)軸導(dǎo)磁，通過定子、定子繞組、轉(zhuǎn)子和轉(zhuǎn)軸等參數(shù)的設(shè)定，可得電機的RMxprt模型，如圖1所示。

仿真參數(shù)設(shè)置：電機的負(fù)載特性直接影響電機驅(qū)動系統(tǒng)的電氣特性，但由于仿真軟件的限制，本文中對穩(wěn)態(tài)工況進行研究，負(fù)載類型為恒功率負(fù)載；額定輸出功率為90kW；額定電壓為250V；額定轉(zhuǎn)速為4 000r/min；工作溫度為75℃。

仿真結(jié)果如圖2和圖3所示。圖2為三相電流曲線，峰值電流約520A；圖3為空氣隙磁通密度曲線，最大磁密約0.73T。

2 Maxwell 2D下的電磁場分析

將建立的電機模型導(dǎo)入至Maxwell 2D模塊進行后續(xù)的有限元仿真。

邊界條件設(shè)置：本文模型在瞬態(tài)磁場求解器下，可將模型的邊線設(shè)置為矢量磁位邊界條件(Vector Potential Boundary或稱為狄利克萊邊界條件)，定義該邊線上的所有點都滿足：

XY坐標(biāo)系A(chǔ)Z=Const

(1)

RZ坐標(biāo)系rAθ=Const

(2)

式中：Const為給定常數(shù)；AZ和Aθ分別為XY坐標(biāo)系下Z方向和RZ坐標(biāo)系下θ方向的矢量磁位。

狄利克萊邊界條件設(shè)置Const=0表示磁力線平行于此邊界線，為理想磁絕緣情況[5]。

將1/8模型的兩條最外側(cè)半徑分別定義為主、從邊界條件(Master/Slave Boundary)，箭頭方向指向無窮遠(yuǎn)處，設(shè)置從邊界條件反對稱于主邊界條件。

激勵源設(shè)置：瞬態(tài)磁場的激勵源比較豐富，有電流源和電流密度源，還可以將導(dǎo)條形成線圈后對線圈施加電流源、電壓源和復(fù)雜控制的外電路源。本文中激勵源選擇繞組電壓源激勵，以Phase A為例(Phase B和Phase C的設(shè)置和Phase A類似，只是電壓源相位有所差異)設(shè)置：多匝繞制，在此為2匝；初始電流值為0；分別設(shè)上整體電阻值、附加電感值(非繞組自身電感值)和電壓激勵源的值，其中電壓激勵源為時間的函數(shù)。

網(wǎng)格剖分：Maxwell 2D網(wǎng)格剖分采用了金字塔型剖分設(shè)置，無需用戶過多的參與，直接利用內(nèi)置的自適應(yīng)剖分也可以得到正確的計算結(jié)果。

系統(tǒng)自帶的網(wǎng)格剖分設(shè)置共有三大項，分別是On Selection、Inside Selection和Surface Approximation，其各自的意義分別為對物體邊界內(nèi)指定剖分規(guī)則、對物體內(nèi)部指定剖分規(guī)則和對物體表層指定剖分規(guī)則。

線圈、電極和主體部分的內(nèi)部剖分均采用Length Based Refinement(基于單元變長的剖分設(shè)置)，設(shè)置單元最大邊長分別為5.5、4和10.1mm。

電極和主體部分的邊界剖分均采用Surface Approximation(表層細(xì)致剖分設(shè)置)，設(shè)置最大弦長分別為0.079 5和0.13mm，弦所對應(yīng)的最大角均為15°。

求解設(shè)置：采用瞬時磁場求解器，設(shè)置停止時刻為0.04s；時間步長為0.04ms；非線性的殘差為0.000 1。

通過邊界條件設(shè)置、激勵源設(shè)置、網(wǎng)格剖分和求解設(shè)置后進行仿真，結(jié)果如圖4所示。當(dāng)前計算時刻為0.04s，轉(zhuǎn)速為4 000r/min，轉(zhuǎn)子位置處于248.46°。由圖可見：運行后電機定子槽內(nèi)漏磁增加；電機在軛部磁密較高。

3 Simplorer平臺下的聯(lián)合仿真

本文中采用異步正弦脈寬調(diào)制方式，調(diào)制時通常保持輸入的三角載波信號頻率ft不變，只改變正弦調(diào)制波信號頻率fr，這樣可提高低頻時的載波比，從而使SPWM逆變器每個輸出電壓周期內(nèi)的SPWM矩形脈沖數(shù)在輸出頻率fo降低情況下增加，兩脈沖間的間距縮小，諧波受到抑制而減小，負(fù)載電機產(chǎn)生的脈動轉(zhuǎn)矩和噪聲減小，有效改善低頻工作時的性能。

3.1 Simplorer與RMxprt的聯(lián)合仿真

理想情況下的聯(lián)合仿真模型如圖5所示：電源采用理想電壓源E1，電壓為250V；濾波電容C1容值選取為9 900μF；功率開關(guān)管采用IGBT加反向續(xù)流二極管，IGBT導(dǎo)通電壓為0.8V；RMxprt模型采用第1節(jié)所設(shè)計的電機模型；采用定轉(zhuǎn)速控制，轉(zhuǎn)速設(shè)定為3 750r/min；控制策略采用正弦脈寬調(diào)制(SPWM)中的異步調(diào)制方式；電流表AM1～AM3測量三相電流，電壓表VM1～VM3測量線電壓，電流表AM4測量直流母線電流，電壓表VM4測量三相相電壓。

三相電流曲線如圖6所示。由圖可見：三相電流收斂穩(wěn)定后峰值電流約為520A，與圖2所示RMxprt模型下的三相電流曲線基本一致。

電機輸出轉(zhuǎn)矩曲線如圖7所示。由圖可見：穩(wěn)定后電機輸出轉(zhuǎn)矩恒定在約200N·m。

但在實際情況下，電池、逆變器、電機和高壓線纜等都存在一些非線性的寄生電容和電感、分布電容和電感、雜散電容和電感等不可見的隱性元素，這些元素給干擾的傳播提供有效途徑，使得該系統(tǒng)的電磁兼容研究變得非常復(fù)雜。

本文中建立了驅(qū)動系統(tǒng)的一個非理想模型：電源采用共模電流模型，存在對地分布電容Cg1；高壓線纜采用T型集總參數(shù)模型；增加電機對地分布電容Cg2；AM4、AM5用于測量共模電流，如圖8所示。

電機驅(qū)動系統(tǒng)共模電流的傳播主要通過以下幾條途徑：電力電子器件與散熱器之間的寄生電容耦合；電機的繞組和定子機殼之間的分布電容耦合；電池和高壓線纜與地之間的分布電容耦合。由于仿真軟件的限制，本文中主要研究流過電池和電機與地之間分布電容的共模電流以及流過高壓線纜與地之間分布電容的共模電流。

通過AM4、AM5測量的直流端高壓線纜和交流端高壓線纜到地的共模電流曲線如圖9所示。由圖可見：直流端高壓線纜到地的共模電流約為16mA，峰值約為36mA；交流端高壓線纜到地的共模電流約為24mA，峰值約為48mA。

流過電池和電機的共模電流曲線如圖10所示。由圖可見：流過電池的共模電流約為45mA，峰值約為90mA；流過電機的共模電流約為38mA，峰值約為80mA。由此可得出另外有一部分共模電流通過逆變器的分布參數(shù)由逆變器到大地形成共?；芈贰?/p>

3.2 Simplorer與Maxwell 2D的聯(lián)合仿真

非理想情況下Simplorer與Maxwell 2D的聯(lián)合仿真模型如圖11所示。電源同樣采用共模電流模型，存在對地分布電容Cg1；高壓線纜同樣采用T型集總參數(shù)模型；功率開關(guān)管采用IGBT加反向續(xù)流二極管，IGBT導(dǎo)通電壓為0.8V；Maxwell 2D模型采用第2節(jié)所設(shè)計的電機模型，電機對地分布電容為Cg2；轉(zhuǎn)速同樣設(shè)定為3 750r/min；控制策略同樣采用正弦脈寬調(diào)制(SPWM)中的異步調(diào)制方式；AM1～AM3用于測量三相電流，AM4和AM5用于測量共模電流。

三相電流曲線如圖12所示。由圖可見：與圖11所示Simplorer與RMxprt聯(lián)合仿真的三相電流曲線基本一致，峰值約為450A。

電機輸出轉(zhuǎn)矩曲線如圖13所示。由圖可見：穩(wěn)定后峰值約為170N·m。

通過AM4和AM5測量的直流端高壓線纜和交流端高壓線纜到地的共模電流曲線如圖14所示。由圖可見：直流端高壓線纜到地的共模電流約為20mA，峰值約為30mA；交流端高壓線纜到地的共模電流約為25mA，峰值約為60mA。

流過電池和電機的共模電流曲線如圖15所示。由圖可見：流過電池的共模電流約為250mA，峰值約為750mA；流過電機的共模電流約為250mA，峰值約為700mA。

下面分析3種EMC設(shè)計對流過電池和電機的共模電流的影響。

(1) 減小分布電容容值

圖16為將電池和電機的分布電容由10nF改為1nF后流過電池和電機的共模電流曲線。由圖可見：與圖15相比，分布電容減小1個數(shù)量級，共模電流就有明顯減少；此時流過電池的共模電流約為100mA，峰值約為175mA；流過電機的共模電流約為100mA，峰值約為140mA。從而驗證了分布電容的大小對共模干擾的大小影響非常大。

(2) 緩沖電路設(shè)計

在直流母線上、IGBT邊分別并聯(lián)一個單電容緩沖電路，如圖17所示，共3個。流過電池和電機的共模電流曲線如圖18所示。由圖可見：流過電池的共模電流約為54mA，峰值約為90mA；流過電機的共模電流也約為54mA，峰值約為78mA?？梢娕c圖16相比均有明顯降低，證實了緩沖電路對改善系統(tǒng)EMC性能有明顯作用。

(3) 濾波器設(shè)計

在直流電源輸入端加入了LISN阻抗穩(wěn)定網(wǎng)絡(luò)，在三相交流輸出端加入LC濾波電路，流過電池和電機的共模電流曲線如圖19所示。由圖可見：流過電池的共模電流幾乎降到0；流過電機的共模電流約為60mA，峰值約為120mA?？梢娕c圖20相比均有明顯降低，證實了濾波器設(shè)計對改善系統(tǒng)EMC性能有明顯作用。

4 結(jié)論

針對電動汽車驅(qū)動系統(tǒng)進行了EMC建模仿真研究，采用磁路法建立RMxprt交流永磁同步電機模

型，運用有限元法分析Maxwell 2D電機模型，并在Simplorer平臺下進行整個驅(qū)動系統(tǒng)的聯(lián)合建模仿真。通過多個平臺下的仿真結(jié)果對比分析，驗證了電動汽車驅(qū)動系統(tǒng)的EMC仿真真實可信；同時聯(lián)合仿真下的非理想?yún)?shù)設(shè)計和EMC設(shè)計對流過電池和電機的共模電流影響的分析，均與理論較符合，證實了驅(qū)動系統(tǒng)EMC仿真設(shè)計的有效性。

參考文獻

[1] 王洪武,張戟,楊騰飛.電動汽車驅(qū)動系統(tǒng)EMC問題研究現(xiàn)狀[J].電子測量技術(shù),2011(34):18-22.

[2] Kremer Frank, Frei Stephan. Simulation of Emissions of Power Electronic Devices in Electrical and Hybrid Electrical Vehicles[C]. 2010 Asia-Pacific Symposium on Electromagnetic Compatibility, APEMC 2010:691-694.

[3] 趙博,張洪亮.Ansoft 12在工程電磁場中的應(yīng)用[M].北京:中國水利水電出版社,2010.

[4] Maki Kohji, Funato Hiroki. Motor Modeling for EMC Simulation by 3-D Electromagnetic Field Analysis[C]. 2009 IEEE International Electric Machines and Drives Conference, IEMDC 2009:103-108.

[5] 張戟,孫澤昌.現(xiàn)代汽車電磁兼容理論與設(shè)計基礎(chǔ)[M].北京:清華大學(xué)出版社,2009.