永磁步進(jìn)電動機(jī)磁場的計(jì)算
郭 宏 王宗培 郭慶吉 劉桂芝(哈爾濱工業(yè)大學(xué))
l 引 言
為了保證步進(jìn)電動機(jī)合理設(shè)計(jì)和可靠運(yùn)行，需運(yùn)用數(shù)值解法計(jì)算電機(jī)的磁場。電機(jī)電磁場的數(shù)值解法中，有限元法最常用。在有元法的整個計(jì)算過程中，非線性方程組的求解占有關(guān)鍵的地位。由于電機(jī)磁路中鐵磁材料部分的磁飽和特性，使得所建立的非線性方程組具有較復(fù)雜的結(jié)構(gòu)形式，這就常常引起非線性方程組求解過程的不收斂性。本文通過對一臺二相永磁步進(jìn)電動機(jī)空載磁場的實(shí)際解算，具體分析非線性方程組求解過程中出現(xiàn)的不收斂性，并采用可變步長的牛頓一拉斐森法解決牛頓一拉斐森法(定步長)大范圍不收斂，使得迭代初值的選取具有較大的任意性。
2二相永磁步進(jìn)電動機(jī)空載磁場有限元計(jì)算的不收斂性
二相永磁步進(jìn)電動機(jī)的子具有24個齒(均布)，轉(zhuǎn)子具有18個磁極(****磁鋼，N、S均布)。據(jù)其結(jié)構(gòu)特點(diǎn)，可以選擇圓周60o范圍為求解區(qū)域，其余區(qū)域內(nèi)磁場分布與其相同或?qū)ΨQ。圖l是這一區(qū)域的剖分圖，共有344個節(jié)點(diǎn)，630個單元。
利用矢量位A求解，邊界I、Ⅱ具有第一類邊界條件(A—O、，邊界Ⅲ、Ⅳ具有半周期邊界條件  在半周期邊界條件的處理上采用了不降階的處理方法，即在  邊界Ⅲ的節(jié)點(diǎn)上建立方程時，考慮邊界Ⅳ內(nèi)側(cè)單元的作用，而在邊界Ⅳ的節(jié)點(diǎn)上建立方程時，亦考慮邊界Ⅲ內(nèi)側(cè)單元的作用[1]。從而得非線性方程組：
剛度矩陣[K]中的元素是通過A—B
 這一順序過程求得的，K值與A值之間具有較復(fù)雜的關(guān)系。 對式(1)采用傳統(tǒng)的牛頓一拉斐森法(定步長)求解，迭代公式是：
[A]的初值[A]（0）的選取是認(rèn)為規(guī)定媒質(zhì)系數(shù)全部為線性（對于非線性材料，取磁化曲線線性段斜率），將非線性方程組直接作為線性議程組求解一次獲得。
利用這種傳統(tǒng)的方法解算二相永磁步進(jìn)電動機(jī)時，發(fā)現(xiàn)迭代是不收斂的，選擇取了不同的初值，仍得不到收斂結(jié)果，具體表現(xiàn)有兩種，一是迭代進(jìn)行到某步以后，A值出現(xiàn)振蕩，達(dá)不到預(yù)定的精度；二是迭代最初幾步就出現(xiàn)齒磁密大于2.0T的情況，計(jì)算無法繼續(xù)。
根據(jù)牛頓一拉非森法的局部收斂性定理，要解決不收斂性，需進(jìn)一步選擇非常接近
真值的初值[2][A](0)。但是，要選擇這樣的初值，需經(jīng)過反復(fù)的試算，計(jì)算量太大，而且，不同的區(qū)域需選擇不同的初值才行，這樣做是很困難的。
通過反復(fù)的摸索，采用了一種稱為可變步長的牛頓一拉非森法，解決了磁場計(jì)算中出現(xiàn)的不收斂性，使得初值的選取是大范圍的，提高了收斂性。
3 可變步長的牛頓-拉非森法的實(shí)際應(yīng)用
可變步長的牛頓-拉非森法迭代公式為：
與傳統(tǒng)的牛頓一拉非森法（定步長）比較，主要區(qū)別在于增加了變步長因子   變步長因子[ω]<k>可隨迭代步數(shù)的不同而自動調(diào)整。[ω](k)的調(diào)整是依據(jù)于每步迭代的誤差[R]-[F](k)而定的。 在對這臺二相永磁步進(jìn)電動機(jī)空載磁場進(jìn)行解算時，發(fā)現(xiàn)位于邊界Ⅲ、Ⅳ上的點(diǎn)191,213,...,301;212,234,...3 |
