摘  要： 詳細敘述了基于ANSYS有限元軟件的APDL語言的三相12/8結(jié)構(gòu)的開關(guān)磁阻電機幾何建模、網(wǎng)格剖分、載荷施加、邊界條件、求解及后處理等步驟，建立了二維開關(guān)磁阻電機的有限元模型，分析隨電機轉(zhuǎn)子位置變化的開關(guān)磁阻電機靜態(tài)磁場，獲得了電感、磁鏈及轉(zhuǎn)矩等數(shù)據(jù)，為開關(guān)磁阻電機的本體優(yōu)化及控制系統(tǒng)設(shè)計提供了理論依據(jù)。
關(guān)鍵詞： 開關(guān)磁阻電機；有限元；電磁分析；APDL語言

    開關(guān)磁阻電機SRM（Switched Reluctance Machine）具有結(jié)構(gòu)簡單堅固、可靠性強、起動轉(zhuǎn)矩大、容錯性好等優(yōu)點，在礦山機械、電動汽車及航空航天等領(lǐng)域得到廣泛應(yīng)用[1-4]。但SRM工作過程是高度非線性的，無法嚴格列寫相電感和相電流的具體解析表達式，這給前期設(shè)計帶來很大困難。因此需要借助有限元軟件對SRM的電磁特性進行分析。
    本文采用ANSYS軟件的APDL語言建立三相12/8結(jié)構(gòu)SRM的二維有限元模型并對其靜態(tài)磁場特性進行分析[5-7]。詳細敘述了幾何建模、網(wǎng)格剖分、載荷施加、邊界條件、求解及后處理等步驟，分析隨電機轉(zhuǎn)子位置變化的靜態(tài)磁場，獲得了電感、磁鏈及轉(zhuǎn)矩等數(shù)據(jù)。
1 有限元模型的建立及求解
    整個電機建模過程分為前處理、求解及后處理3部分，如圖1所示。

1.1 模型建立
    前處理過程是有限元分析電機磁場的一個重要過程。由于本文采用ANSYS軟件的APDL語言建立SRM的二維有限元模型，其相對菜單圖形用戶界面GUI（Graphical User Interface）操作的優(yōu)點在于項目的可移植性強。在建立模型時，把項目中使用的電機結(jié)構(gòu)尺寸、電機的材料屬性等參數(shù)變量化，并根據(jù)不通過開發(fā)項目設(shè)定賦具體的值。另外把項目文件名以及建模用的單元等信息放在程序前面部分。這樣，如果再開發(fā)另外一個類似項目時，只需要修改項目文件名、建模單元信息、電機尺寸、材料屬性等參數(shù)變量值即可，項目程序的可移植性非常強。
    建立待建模電機的幾何模型是前處理過程的一種重要環(huán)節(jié)。建模的原則是由點生線，由線生面。由于SRM電機截面的對稱性，在建立模型時，定子部分可以先建立一個定子齒和繞在該齒上的線圈，然后利用Agen命令語句進行12份30°軸中心陣列，即可得到完整的12/8結(jié)構(gòu)SRM定子部分模型。同理，轉(zhuǎn)子部分亦可先建立一個轉(zhuǎn)子齒及其齒中心對稱線45°范圍內(nèi)模型，然后利用Agen命令語句進行8份45°軸中心陣列，即可得到完整的12/8結(jié)構(gòu)SRM轉(zhuǎn)子部分模型。電機定子部分包括定子、繞組和靠近定子的一半定轉(zhuǎn)子間氣隙面積。轉(zhuǎn)子部分包括轉(zhuǎn)子、轉(zhuǎn)軸、轉(zhuǎn)子齒間槽面積和靠近轉(zhuǎn)子的一半定轉(zhuǎn)子間氣隙面積。圖2所示為標有AREA單元標號的12/8結(jié)構(gòu)SRM的幾何模型。

    模型建立后，即可按照實際電機特性對各AREA單元進行賦材料屬性。完成此步驟后，即可進行剖分操作。剖分操作時應(yīng)避免模型中出現(xiàn)尖角(如倒角部分)；避免模型中出現(xiàn)面積比過大的區(qū)域，如氣隙與周圍區(qū)域；避免模型中出現(xiàn)極不規(guī)則的形狀，如包括氣隙在內(nèi)的所有空氣區(qū)域。對電機而言，應(yīng)由氣隙開始剖分，即先進行電機定轉(zhuǎn)子間的氣隙的剖分。圖3為電機剖分圖，其中圖4為定轉(zhuǎn)子氣隙部分剖分圖。

1.2 求解
    若求解m種不同位置時n種電流下的電機磁場數(shù)據(jù)，需要分別建立m×n次模型一步一步地求解，非常繁瑣。為了提高計算效率，建立模型時，將電機轉(zhuǎn)子部分作為一個旋轉(zhuǎn)組件，并設(shè)定一個旋轉(zhuǎn)角度變量，這樣每次程序執(zhí)行時只需要為該旋轉(zhuǎn)角度變量賦值，將此旋轉(zhuǎn)組件旋轉(zhuǎn)至相應(yīng)位置，而不需要重新編寫程序建立電機模型，最后將屬于轉(zhuǎn)子部分旋轉(zhuǎn)組件的氣隙圓環(huán)與屬于定子部件的氣隙圓環(huán)臨近單元合并，同時利用一個循環(huán)在繞組中加載不同電流即可求解出相應(yīng)位置時不同電流下的電感、磁鏈等電磁數(shù)據(jù)。這樣SRM在m種不同位置時n種電流下的磁場數(shù)據(jù)只需要經(jīng)過m次的程序計算即可。按照以上思想確定好位置循環(huán)和電流循環(huán)之間的關(guān)系后，一般把電流環(huán)嵌套于位置環(huán)中為宜。
    另外在定子軛部最外層的所有單元施加矢量磁位Az=0條件，確保電機模型外沒有磁場分布。模型的求解器采用直接波前法求解器，其命令語句為：“Eqslv, Front”。為了提高求解精度，每一個加載計算步的收斂迭代次數(shù)設(shè)為40。
2 仿真結(jié)果及分析
    求解結(jié)束后，可以用菜單操作或者命令行顯示SRM在相應(yīng)位置及相應(yīng)相電流勵磁時的靜態(tài)磁場分布，如圖5所示為12/8結(jié)構(gòu)SRM在定轉(zhuǎn)子對齊位置和不對齊位置時的磁力線分布圖。
    圖6所示為對應(yīng)圖5（a）定子齒和轉(zhuǎn)子齒對齊位置和圖5（b）定子齒和轉(zhuǎn)子槽對齊位置時的磁通密度矢量圖。

    模型程序中，把在不同位置時不同電流的SRM相電感數(shù)據(jù)、相磁鏈及轉(zhuǎn)矩數(shù)據(jù)分別寫入L.dat、ψ.dat和T.dat文件中，打開后提取數(shù)據(jù)到專業(yè)繪圖軟件(如Origin)即可得到圖7～圖9所示的電機電感L(θ，i)、磁鏈ψ(θ，i)和轉(zhuǎn)矩T(θ，i)特性曲線簇。電感L(θ，i)、磁鏈ψ(θ，i)和轉(zhuǎn)矩T(θ，i)的數(shù)據(jù)可為MATLAB/Simulink等方法搭建SRM的動態(tài)性能分析模型提供參考數(shù)據(jù)。

    本文詳細敘述了利用ANSYS軟件的APDL語言建立三相12/8結(jié)構(gòu)SRM的二維有限元模型進行靜態(tài)電磁分析的幾何建模、網(wǎng)格剖分、載荷施加、邊界條件、求解及后處理等步驟，實現(xiàn)了轉(zhuǎn)子組件的自動旋轉(zhuǎn)并與定子組件單元耦合的功能。利用ANSYS軟件的APDL語言建立電機模型的程序可移植性強、求解效率高，計算結(jié)果對電機的靜態(tài)性能分析具有參考價值。
參考文獻
[1] 陳昊．開關(guān)磁阻調(diào)速電動機的原理·設(shè)計·應(yīng)用[M].徐州：中國礦業(yè)大學出版社，2000．
[2] GOPALAKRISHNAN S，OMEKANDA A M,LEQUESNE B．Classification and remediation of electrical faults in the switched reluctance drive[J]．IEEE Transactions on Industry Applications，2006，42(2)：479-486．
[3] 陳小元，鄧智泉，許培林，等．整距繞組分塊轉(zhuǎn)子開關(guān)磁阻電機的電磁設(shè)計[J].中國電機工程學報，2011，31(36)：109-115．
[4] 陳小元，鄧智泉，范娜，等.雙極性分塊轉(zhuǎn)子開關(guān)磁阻電機[J].電機與控制學報，2010，14(10)：1-7.
[5] 嚴云.基于ANSYS參數(shù)化設(shè)計語言的結(jié)構(gòu)優(yōu)化設(shè)計[J].華東交通大學學報，2004，21(4)：52-55.
[6] 丁文，周會軍，魚振民.基于ANSYS二次開發(fā)的開關(guān)磁阻電機電磁場分析軟件[J].微電機，2006，39(2)：19-21，55.
[7] 劉蕓蕓，吳建華，張式勤.基于ANSYS平臺開發(fā)開關(guān)磁阻電機電磁分析軟件[J].中小型電機，2005，32(3)：1-4.