電樞分段供電永磁直線同步電機(jī)的非線性數(shù)學(xué)模型

郭科宇 李耀華 史黎明 周世炯 徐 飛

電樞分段供電永磁直線同步電機(jī)的非線性數(shù)學(xué)模型

郭科宇1,2李耀華1,2史黎明1周世炯1,2徐 飛1

（1. 中國科學(xué)院電工研究所中國科學(xué)院電力電子與電氣驅(qū)動重點實驗室 北京 100190 2. 中國科學(xué)院大學(xué) 北京 100049）

電樞分段供電永磁直線同步電機(jī)（PMLSM）中，相鄰分段鐵心的存在會導(dǎo)致氣隙磁場脈振，動子穿過電樞定子段也會引發(fā)復(fù)雜的暫態(tài)過程。為準(zhǔn)確描述電機(jī)的動態(tài)過程，該文考慮以上兩種特有的電磁特性，建立雙三相分段供電PMLSM的相空間模型。首先，采用解析法分析脈振氣隙磁場的成因，研究電感不平衡以及飽和效應(yīng)問題，將六相電流轉(zhuǎn)化到ab坐標(biāo)系下建立二維電感查詢表，有效地降低了模型的運(yùn)算量。其次，引入易于求導(dǎo)的Sigmoid函數(shù)構(gòu)造隨動子位置變化的耦合系數(shù)，準(zhǔn)確描述動子通過分段過程產(chǎn)生的空載磁鏈及反電動勢。依據(jù)上述方法構(gòu)建電機(jī)相空間模型，與有限元模型進(jìn)行對比，兩者電壓與推力計算結(jié)果一致。最后，在RT-Lab中搭建基于相空間模型的半實物實驗平臺，結(jié)果表明，該模型能夠準(zhǔn)確地體現(xiàn)電機(jī)的各種效應(yīng)，可以替代計算成本較大的有限元模型。

永磁直線同步電機(jī) 分段供電 磁場飽和效應(yīng) 相空間模型

0 引言

雙三相永磁直線同步電機(jī)（Permanent Magnet Linear Synchronous Machine, PMLSM）由于具有響應(yīng)速度高、定位精度高、控制品質(zhì)好、系統(tǒng)可靠性高的優(yōu)點[1-4]，同時兼具多相電機(jī)單相容量低、可靠性高的特點[5-7]，適用于高速磁懸浮交通直線驅(qū)動、高加速的電磁驅(qū)動系統(tǒng)。對于長行程應(yīng)用場合，為了減小電源容量，降低制造和運(yùn)營成本，通常采用長電樞短勵磁的方案，并將長電樞側(cè)分段，當(dāng)動子運(yùn)動到某一個分段下，則為該電樞供電[8-10]。

在研究電機(jī)動態(tài)過程或各種控制方法時，往往需要準(zhǔn)確、全面、高效的動態(tài)數(shù)學(xué)模型[11]。動態(tài)數(shù)學(xué)模型主要分為分布參數(shù)模型和集總參數(shù)模型。有限元模型作為一種分布參數(shù)模型，因可以給出十分精確的計算結(jié)果而被廣泛使用[12-14]；然而有限元模型計算成本大，因而對于行程較長的直線電機(jī)來說，進(jìn)行一次仿真需要的時間過長。而傳統(tǒng)的集總參數(shù)模型因具有計算成本小、計算效率高的特點也得到了廣泛應(yīng)用[15]。電機(jī)的集總參數(shù)模型主要包括dq軸模型和相空間模型[11]，其中dq軸模型是將相空間模型的電感矩陣通過Park正交變換成對角矩陣，從而達(dá)到dq軸解耦與簡化計算的目的。對于電感矩陣循環(huán)對稱且為常數(shù)或者僅為轉(zhuǎn)子或動子位置正弦函數(shù)的情況，dq軸模型易于實施且容易計算[16-18]，但是當(dāng)考慮飽和效應(yīng)等非線性因素以及諧波磁場時，dq軸模型中的電感矩陣便不再是常量，此時相空間模型憑借其靈活性更加適合建立動態(tài)數(shù)學(xué)模型[15, 19-20]。文獻(xiàn)[15]采用相磁鏈直接查表法建立了三相旋轉(zhuǎn)永磁電機(jī)的非線性相空間模型，減小了查表的維度，降低了模型計算量。文獻(xiàn)[19]采用三相電感直接查表法針對星形聯(lián)結(jié)的旋轉(zhuǎn)永磁電機(jī)建立了相空間模型，并通過有限元計算驗證了該模型的正確性。文獻(xiàn)[20]對三相盤式電機(jī)建立了非線性相空間模型，并同時在模型中考慮了鐵心損耗。然而以上文獻(xiàn)集中于旋轉(zhuǎn)電機(jī)的動態(tài)模型，針對鐵心開斷且分段供電的長電樞永磁直線電機(jī)動態(tài)數(shù)學(xué)模型的文獻(xiàn)較少。

與旋轉(zhuǎn)電機(jī)不同的是，長電樞永磁直線電機(jī)存在磁極動子穿過定子電樞分段的過程，同時還存在由相鄰分段電鐵心引起脈振氣隙磁場而導(dǎo)致的電感不平衡的現(xiàn)象[21-28]。文獻(xiàn)[28-29]對分段供電三相永磁直線電機(jī)建立了相空間模型，在模型中同時考慮了動子的通過過程，并分析了電機(jī)運(yùn)行的動態(tài)過程。但是由于該電機(jī)主要工作在持續(xù)運(yùn)行工況，因而飽和程度不明顯，因此該模型假設(shè)動子完全被定子覆蓋時電感為恒定值。

然而由于高加速電磁驅(qū)動應(yīng)用中的直線電機(jī)同時需要滿足輕量化和大推力，往往工作在短時高過載的工況，因而面臨深度飽和效應(yīng)引發(fā)的電感波動問題。在短電樞直線電機(jī)和旋轉(zhuǎn)電機(jī)中，飽和效應(yīng)主要由行波磁場引發(fā)，電感波動循環(huán)對稱；而長電樞分段供電PMLSM中的飽和效應(yīng)主要由脈振氣隙磁場引起，其結(jié)果是各相電感諧波幾乎為同相位。采用精確且計算快速的數(shù)學(xué)模型描述電機(jī)的動態(tài)過程對分段供電PMLSM電磁瞬態(tài)過程研究、電機(jī)實施精確控制都具有重要意義。

本文建立準(zhǔn)確且計算快速的雙三相分段供電PMLSM的相空間模型，針對定子電樞分段造成的脈振磁場，采用解析法對其進(jìn)行深入分析。提出將電機(jī)六相電流轉(zhuǎn)化到ab坐標(biāo)系下建立二維電感查詢表，在保證查表精度的同時，有效提升電感查表速度，降低生成電感查表所需的有限元仿真或?qū)嶒灤螖?shù)。針對動子過分段的過程，提出采用易于求導(dǎo)的Sigmoid函數(shù)構(gòu)造耦合因數(shù)，簡單并準(zhǔn)確地描述了電機(jī)的空載磁鏈和反電動勢特性。最后基于RT-Lab搭建了半實物實驗平臺，驗證了本文建立的相空間模型能夠準(zhǔn)確反映雙三相分段供電PMLSM的電磁特性。

1 電感不平衡和飽和特性分析及電感快速查詢表構(gòu)建方法

雙三相長電樞直線電機(jī)定子由多個分段組成，設(shè)從左向右的分段編號為1～sec，分段供電直線電機(jī)如圖1所示，一個分段下電機(jī)的結(jié)構(gòu)如圖2所示。圖中，PMLSM參數(shù)見表1。為節(jié)約橫向端部空間，方便繞線，定子采用背繞式分布繞組（即克萊姆繞組）。為了減少動子永磁體端部漏磁，動子采用端部半磁極結(jié)構(gòu)。由于電機(jī)工作在短時狀態(tài)，因此電流密度選取較高。

圖1 分段供電直線電機(jī)

圖2 一個分段電機(jī)的結(jié)構(gòu)

表1 PMLSM參數(shù)

Tab.1 The parameters of the PMLSM

1.1 脈振氣隙磁場分析

永磁體的磁導(dǎo)率與空氣相近，且動子中不存在鐵磁材料，可單獨(dú)研究電流在氣隙中產(chǎn)生的磁場。繞組在氣隙中產(chǎn)生的磁場可視為各個導(dǎo)體單獨(dú)產(chǎn)生的磁場疊加的結(jié)果，因此先計算單個導(dǎo)體在氣隙中產(chǎn)生的磁場，單個導(dǎo)體位于雙邊鐵心中的模型如圖3所示。

為簡化計算過程，進(jìn)行如下假設(shè)：①忽略齒槽，設(shè)導(dǎo)體內(nèi)的電流為，將單根載流導(dǎo)體視作一個緊挨鐵心的點電流；②＜-1或＞2處的磁場為0；③氣隙中只含方向磁場；0為真空磁導(dǎo)率；Fe為鐵心磁導(dǎo)率；⑤忽略分段間軌縫。

圖3 單個載流導(dǎo)體位于雙邊鐵心中

分別對區(qū)域2與區(qū)域3列寫關(guān)于法向磁通密度的微分方程為

其中

式中，1y、2y分別為區(qū)域a和區(qū)域b的氣隙法向磁通密度。區(qū)域a與區(qū)域b滿足的邊界條件為

求解式（1）、式（2）得到單個導(dǎo)體在雙邊鐵心下產(chǎn)生的氣隙方向磁通密度分布cony為

一相繞組通電后在氣隙中產(chǎn)生的磁場為繞組中各導(dǎo)體產(chǎn)生磁場的疊加，各相繞組產(chǎn)生的氣隙磁通密度表達(dá)式為

對應(yīng)不同時刻（對應(yīng)不同電角度）下的整體氣隙磁通密度波形如圖4所示。氣隙磁場在向前運(yùn)動的同時，還存在直流脈振分量。當(dāng)電角度為p/2和3p/2時，氣隙磁場的直流分量分別達(dá)到正、負(fù)最大值。由于供電分段的兩側(cè)存在未通電分段的鐵心，根據(jù)磁通連續(xù)性定理，各相繞組產(chǎn)生的磁場在通電分段區(qū)間內(nèi)存在直流分量，因此導(dǎo)致合成磁場中存在直流脈振分量[21]，從而形成脈振氣隙磁場。

圖4 不同電角度下的六相電流產(chǎn)生的氣隙磁通密度

1.2 脈振氣隙磁場導(dǎo)致的電感不對稱現(xiàn)象

雙三相電機(jī)的繞組相量如圖5所示，ABC和DEF分別為兩套互差120°電角度的三相繞組，兩套三相繞組間相差30°電角度。

圖5 雙三相電機(jī)的繞組相量

雙三相電機(jī)的電感可依據(jù)相間電角度分為五類。在旋轉(zhuǎn)電機(jī)中，當(dāng)忽略凸極效應(yīng)及飽和因素時，各個類別中的電感是相等的。

（1）第一類電感：各相自感，包括aa、bb、cc、dd、ee、ff。

（2）第二類電感：互差120°相間互感，包括ab、bc、ac、de、ef、df。

（3）第三類電感：互差30°相間互感，包括ad、be、cf。

（4）第四類電感：互差90°相間互感，包括bd、ce、af。

（5）第五類電感：互差150°相間互感，包括ae、cd、bf。

對于分段供電的雙三相直線電機(jī)，當(dāng)一相繞組通入電流時，繞組產(chǎn)生的氣隙磁通密度存在直流分量。這部分直流分量交鏈到其他相繞組中會導(dǎo)致電感不平衡。直線電機(jī)的自感與互感表達(dá)式均為

式中，1、2為相標(biāo)號，即a～f；為2相的繞組函數(shù)[30]；sw為有效電機(jī)橫向?qū)挾?。?dāng)1=2時，式（5）計算得到1相自感；當(dāng)1≠2時，式（5）計算得到1與2相間互感。

采用式（5）計算得到電機(jī)的五類電感見表2。第一類和第三類電感是內(nèi)部相等的，但第二、四、五類電感內(nèi)部均不相等。在旋轉(zhuǎn)電機(jī)中應(yīng)為0的第四類電感，在長初級直線電機(jī)中不為0。

1.3 脈振氣隙磁場導(dǎo)致的飽和效應(yīng)及電感波動分析

電機(jī)工作于短時高過載狀態(tài)，繞組的電流密度較大（最大達(dá)到44A/mm2），容易導(dǎo)致鐵心飽和。六相電流產(chǎn)生的氣隙磁場存在脈振分量，使得鐵心軛部的飽和程度加深。忽略永磁體，單獨(dú)向供電分段的繞組中通入一個電周期的六相電流，獲得電角度為0、p/2、p、3p/2下的鐵心磁通密度分布如圖6所示。由于脈振磁場的影響，電角度為p/2和3p/2時，鐵心局部磁通密度大幅增加，達(dá)到2T以上；而當(dāng)電角度為0和p時，鐵心磁通密度僅為1.6T。

表2 一個分段的各類電感

Tab.2 The inductances of a power section

圖6 一個電周期內(nèi)不同電角度下的磁通密度分布

一個電周期內(nèi)，脈振磁場兩次達(dá)到峰值，此時由于偏置引發(fā)的軛部磁通密度飽和程度加深則會導(dǎo)致各類電感下降，引發(fā)電感波形中產(chǎn)生2次諧波。且對于各相電感，2次諧波相位基本一致。通入額定六相電流，一個電周期內(nèi)的五類動態(tài)電感波形如圖7所示。

1.4 任意六相電流下的二維電感查詢表

雙三相電機(jī)的動態(tài)電感矩陣階數(shù)=6，共有2個元素，考慮到電感矩陣的對稱性，以及由相鄰未通電鐵心引起的電感矩陣不平衡及飽和效應(yīng)，動態(tài)電感矩陣中含有(2-)/2+=21個獨(dú)立變量。該21個變量均為六相電流（a～f）的函數(shù)，即(a～f)，因此動態(tài)電感矩陣的查詢表為一個6輸入、21輸出系統(tǒng)。

圖7 額定電流下一個電周期內(nèi)的動態(tài)電感

采用拉丁超立方采樣（Latin Hypercube Sampling, LHS）在ab 平面中進(jìn)行電流采樣，獲得各個電感隨ia 與ib 變化的插值曲面，如圖8所示。

2 動子過分段過程空載特性分析及建模

2.1 基于Sigmoid函數(shù)的耦合因數(shù)構(gòu)造

動子采用的是端部半磁極形式，以半個N極與半個S極為一個單元磁路。取分段中的A相繞組為例，動子通過供電分段的過程如圖9所示。將運(yùn)動坐標(biāo)軸原點定于供電分段的左側(cè)端部，表示動子頭部位置。隨著動子進(jìn)入分段，與A相繞組耦合的單元磁路個數(shù)逐漸增多，直到A相繞組耦合的單元磁路達(dá)到最大極數(shù)6，動子穿出過程與穿入過程相反。永磁體在A相中產(chǎn)生的磁鏈pmA可分解為圖9中的A相繞組耦合系數(shù)A與A相單元磁路磁鏈pmuA的乘積，即pmA=ApmuA。

圖9 動子通過一個分段時在A相繞組中產(chǎn)生的磁鏈

A相的耦合因數(shù)A為動子頭部位置的函數(shù)。由于在求取空載反電動勢時，需要對永磁體磁鏈求導(dǎo)，若A采用階躍函數(shù)則表示會在求導(dǎo)過程中產(chǎn)生沖激。為避免該現(xiàn)象，引入Sigmoid函數(shù)（即sig()）柔化A的階躍邊緣，其導(dǎo)數(shù)可用函數(shù)本身的值求出，即sig()=sig()(1-sig())。A相耦合因數(shù)為

式中，sig為Sigmoid函數(shù)的伸縮因子，sig值越大，則階躍邊緣越銳化。其余相的耦合因數(shù)可以通過平移A獲得，即

動子運(yùn)動過程中，耦合因數(shù)隨動子位置的變化如圖10所示。

圖10 耦合因數(shù)隨動子位置的變化

2.2 基于耦合因數(shù)的空載特性建模

由圖9可知，動子在分段繞組中產(chǎn)生的永磁體磁鏈可由單元磁路磁鏈與相應(yīng)的耦合系數(shù)相乘得到。以A相為例，一個單元磁路在一個周期內(nèi)產(chǎn)生的磁鏈及頻譜如圖11所示。

圖11 單元磁路在A相中產(chǎn)生的永磁體磁鏈

A相單元磁路永磁體磁鏈的傅里葉級數(shù)形式為

式中，pmui與pmui分別為第次諧波的幅值和相位。其余各相的單元磁路永磁體磁鏈可通過平移pmuA得到。六相永磁體磁鏈可表示為

動子穿過分段時產(chǎn)生的各相空載反電動勢為永磁體磁鏈對時間的導(dǎo)數(shù)，即

式中，p為各相符號，p=A～F；v為動子速度。將采用式（10）與式（11）計算得到的空載永磁體磁鏈與空載反電動勢分別與有限元計算結(jié)果對比如圖12所示，有限元和空載模型的計算結(jié)果一致。

3 電機(jī)的相空間模型

為描述電機(jī)整體運(yùn)動過程，基于電機(jī)的動態(tài)電感查表方法、空載磁鏈及反電動勢描述方法建立了電機(jī)的狀態(tài)方程。一個分段下，電機(jī)的電壓方程為

式中，abcdef、abcdef和pm分別為六相相電壓、相電流和空載磁鏈?zhǔn)噶?；s為相電阻；、s分別為六相動態(tài)電感矩陣和漏感矩陣。采用虛功原理計算電機(jī)的電磁推力為

選取動子速度、六相電壓abcdef為控制量，abcdef為狀態(tài)量建立一個分段下的相空間模型為

式中，thrust為電機(jī)產(chǎn)生的電磁推力；detent為動子行進(jìn)過程中由端部效應(yīng)、齒槽效應(yīng)引起的定位力。

一個分段下的電機(jī)相空間模型可更直觀地以框圖形式表示，電機(jī)狀態(tài)方程的框圖如圖13所示。

圖13 電機(jī)狀態(tài)方程的框圖

通入六相平衡額定電流，采用式（14）對動子穿過一個分段時的負(fù)載性能進(jìn)行計算，并將端口電壓和推力波形與有限元結(jié)果對比如圖14所示。為檢驗結(jié)果的一致性，計算相空間模型的各相電壓波形與有限元之間的皮爾森相關(guān)系數(shù)（Pearson Correlation Coefficient, PCC）見表3，六相電壓相關(guān)系數(shù)均在99%以上，說明由六相電感矩陣不平衡和飽和效應(yīng)所導(dǎo)致的電壓不平衡均準(zhǔn)確地體現(xiàn)在相空間模型中。采用相空間計算得到的動子穿過分段時的推力波形與有限元計算結(jié)果的相關(guān)系數(shù)為99.6%，表明相空間模型可以精確地描述動子的過分段推力特性。

表3 相空間模型電壓波形與有限元計算的相關(guān)系數(shù)

Tab.3 The correlation coefficient between phase voltage calculated by phase-domain model and FEM-model

4 實驗驗證

為驗證電機(jī)數(shù)學(xué)模型的有效性，在RT-LabOP5700半實物平臺中搭建了電機(jī)相空間模型，并基于該模型實現(xiàn)電機(jī)全過程實時驅(qū)動系統(tǒng)。實驗平臺及驅(qū)動系統(tǒng)整體架構(gòu)如圖15所示。實驗機(jī)箱由CPU板卡和Xilinx公司的Virtex7 FPGA板卡兩部分組成。電機(jī)的驅(qū)動控制器采用CPU板卡實現(xiàn)，控制頻率為10kHz。電機(jī)的相空間模型則編寫在FPGA中，其離散化精度為500ns。

由于分段供電的電機(jī)最多有2個分段同時供電，為節(jié)省計算資源，在模型中僅需兩臺雙三相變流器對兩個電機(jī)相空間模型交替供電，長電樞電機(jī)的分段供電方式如圖16所示。兩個電機(jī)相空間模型分別代表實際系統(tǒng)中的奇數(shù)號分段和偶數(shù)號分段。如圖15b所示，在FPGA中，包含奇數(shù)號分段和偶數(shù)號分段的電機(jī)模型以及動子運(yùn)動方程。電機(jī)模型計算得到電磁推力后作為輸入送入運(yùn)動方程中進(jìn)行動子位置和速度的計算。動子運(yùn)動方程將動子位置和速度反饋到電機(jī)模型完成空載反電動勢與Park變換矩陣的計算。

圖15 基于RT-Lab OP5700的半實物實驗平臺及驅(qū)動系統(tǒng)整體架構(gòu)

圖16 長電樞電機(jī)的分段供電方式

為節(jié)省計算資源，在建立基于FPGA的電機(jī)模型時，所有數(shù)據(jù)均采用定點數(shù)進(jìn)行運(yùn)算。同時，復(fù)雜的函數(shù)（如Sigmoid函數(shù)及其導(dǎo)數(shù)、三角函數(shù)等）均采用查表的方式進(jìn)行計算。兩個電機(jī)方程以及一個動子運(yùn)動方程總共占用FPGA的主要資源見表4。FPGA中的電機(jī)模型離散化精度為500ns，而FPGA的晶振頻率為200MHz。本文建立的電機(jī)模型可以在100個晶振周期之內(nèi)完成計算。

表4 電機(jī)模型以及動子運(yùn)動方程所占用的FPGA資源

Tab.4 The resources in FPGA occupied by PMLSM model and motion equation

采用d=0的控制策略，將280kg的動子以90m/s2的加速度跟蹤最大為30m/s的三角波速度給定。采用數(shù)學(xué)模型計算得到的電機(jī)加速性能如圖17所示。圖17a表明，整個過程中動子成功地跟蹤了參考速度。在加速過程中，不同的供電分段交替為動子提供推力，如圖17b所示，可以看出，當(dāng)動子過分段時，電機(jī)的相空間模型可以準(zhǔn)確體現(xiàn)動子耦合因數(shù)下降引起的推力下降，以及由定位力引發(fā)的波動。電機(jī)dq軸電流給定及各分段的dq軸電流波形如圖17c所示。圖17d為各個分段六相電流波形。電機(jī)運(yùn)動過程中的空載反電動勢如圖17e所示，隨著速度提高，空載反電動勢幅值隨之升高，同時在動子過分段時，動子所離開的分段空載反電動勢逐漸減小為0，因此所建立的模型可以準(zhǔn)確描述動子過分段時的反電動勢特性。

圖17 采用數(shù)學(xué)模型計算得到的電機(jī)加速性能

以上結(jié)果表明，本文建立的相空間模型不僅保證了與有限元模型接近的精度，同時作為集中參數(shù)模型，相空間模型的計算成本小于有限元模型，可以方便地搭建在半實物實驗平臺中，提高控制系統(tǒng)調(diào)試和分析的準(zhǔn)確性。

5 結(jié)論

本文建立了電機(jī)的非線性相空間模型，可以準(zhǔn)確地描述多相分段供電PMLSM中由繞組分段引起的電感不平衡、飽和效應(yīng)和動子通過分段時的暫態(tài)過程。主要結(jié)論如下：

1）由于通電電樞段兩側(cè)存在未通電電樞段的鐵心，從而導(dǎo)致氣隙磁通密度中含有脈振分量，脈振分量導(dǎo)致相電感存在不平衡現(xiàn)象。

2）在高過載工況下，電機(jī)的飽和效應(yīng)主要由脈振磁場引起，從而使得各相自感和互感同步波動。六相電機(jī)的電感矩陣階數(shù)較高，將多相電流從多相坐標(biāo)系映射到ab坐標(biāo)系，通過a和b查詢的二維動態(tài)電感查詢表，可有效地削減電感查詢表構(gòu)建所需的仿真或?qū)嶒灁?shù)量，降低了模型的計算成本。

3）采用Sigmoid函數(shù)構(gòu)造耦合因數(shù)，并基于耦合因數(shù)和單元磁路永磁體磁鏈建立的動子，通過分段時的空載磁鏈及空載反電動勢模型，準(zhǔn)確地描述電機(jī)特性。

4）建立了一個分段下的電機(jī)狀態(tài)方程，對比負(fù)載狀態(tài)下相空間模型和有限元模型的電壓與推力，結(jié)果表明二者計算結(jié)果的一致性為99%以上。

通過實驗驗證的相空間模型可有效地應(yīng)用于相應(yīng)控制系統(tǒng)調(diào)試和分析中，相空間模型不僅保證了近于有限元模型的精度，同時作為集中參數(shù)模型，其計算成本遠(yuǎn)小于有限元模型。

[1] 張邦富, 程明, 王颯颯, 等. 基于改進(jìn)型代理模型優(yōu)化算法的磁通切換永磁直線電機(jī)優(yōu)化設(shè)計[J]. 電工技術(shù)學(xué)報, 2020, 35(5): 1013-1021.

Zhang Bangfu, Cheng Ming, Wang Sasa, et al. Optimal design of flux-switching permanent magnet linear machine based on improved surrogate-based optimization algorithm[J]. Transactions of China Electrotechnical Society, 2020, 35(5): 1013-1021.

[2] 曹瑞武, 蘇恩超, 張學(xué). 軌道交通用次級分段型直線磁通切換永磁電機(jī)研究[J]. 電工技術(shù)學(xué)報, 2020, 35(5): 1001-1012.

Cao Ruiwu, Su Enchao, Zhang Xue. Investigation of linear flux-switching permanent magnet motor with segmented secondary for rail transit[J]. Transactions of China Electrotechnical Society, 2020, 35(5): 1001- 1012.

[3] 王明杰, 徐偉, 楊存祥, 等. 基于精確子域模型的永磁直線同步電機(jī)空載磁場解析計算[J]. 電工技術(shù)學(xué)報, 2020, 35(5): 942-953.

Wang Mingjie, Xu Wei, Yang Cunxiang, et al. Analytical calculation of no-load magnetic field in permanent magnet linear synchronous motors based on an accurate subdomain model[J]. Transactions of China Electrotechnical Society, 2020, 35(5): 942- 953.

[4] 許孝卓, 孫震, 汪旭東, 等. Halbach交替極永磁同步直線電機(jī)特性分析[J]. 電工技術(shù)學(xué)報, 2019, 34(9): 1825-1833.

Xu Xiaozhuo, Sun Zhen, Wang Xudong, et al. Characteristic of a novel Halbach array consequent- pole[J]. Transactions of China Electrotechnical Society, 2019, 34(9): 1825-1833.

[5] 麥志勤, 肖飛, 劉計龍, 等. 基于準(zhǔn)比例諧振級聯(lián)PI的雙三相永磁同步電機(jī)諧波電流抑制策略[J]. 電工技術(shù)學(xué)報, 2018, 33(24): 5751-5759.

Mai Zhiqin, Xiao Fei, Liu Jilong, et al. Harmonic current suppression strategy of dual three-phase permanent magnet synchronous motor based on quasi proportional resonant cascading PI[J]. Transactions of China Electrotechnical Society, 2018, 33(24): 5751-5759.

[6] 高闖, 趙文祥, 吉敬華, 等. 低諧波雙三相永磁同步電機(jī)及其容錯控制[J]. 電工技術(shù)學(xué)報, 2017, 32(增刊1): 124-130.

Gao Chuang, Zhao Wenxiang, Ji Jinghua, et al. Low harmonic dual three-phase permanent magnet synchronous motor and fault-tolerant control[J]. Transactions of China Electrotechnical Society, 2017, 32(S1): 124-130.

[7] 李小慶, 朱景偉, 孫軍浩, 等. 雙繞組永磁容錯電機(jī)矢量控制系統(tǒng)研究[J]. 電工技術(shù)學(xué)報, 2016, 31(5): 26-34.

Li Xiaoqing, Zhu Jingwei, Sun Junhao, et al. Study on the vector control system for dual winding fault- tolerant permanent magnet motors[J]. Transactions of China Electrotechnical Society, 2016, 31(5): 26-34.

[8] 張明元, 馬偉明, 何娜. 長初級直線電動機(jī)分段供電技術(shù)綜述[J]. 中國電機(jī)工程學(xué)報, 2013, 33(27): 96-104, 14.

Zhang Mingyuan, Ma Weiming, He Na. Application of block feeding methods in long primary linear motors[J]. Proceedings of the CSEE, 2013, 33(27): 96-104, 14.

[9] 馬名中, 馬偉明, 范慧麗, 等. 長初級直線感應(yīng)電機(jī)分段供電切換暫態(tài)過程[J]. 電機(jī)與控制學(xué)報, 2015, 19(9): 1-7.

Ma Mingzhong, Ma Weiming, Fan Huili, et al. Switch transient process of section powered long- primary linear induction motor[J]. Electric Machines and Control, 2015, 19(9): 1-7.

[10] 張明元, 馬偉明, 徐興華, 等. 一種考慮電流過零的直線電機(jī)分段供電策略[J]. 海軍工程大學(xué)學(xué)報, 2019, 31(4): 11-16.

Zhang Mingyuan, Ma Weiming, Xu Xinghua, et al. A block feeding strategy for linear motor considering switching at current-crossing point[J]. Journal of Naval University of Engineering, 2019, 31(4): 11-16.

[11] Quéval L, Ohsaki H. Nonlinear abc-model for electrical machines using N-D lookup tables[J]. IEEE Transactions on Energy Conversion, 2015, 30(1): 316-322.

[12] Leonardo L D, Parasiliti F, Tursini M, et al. Transient analysis of PM synchronous motor drives by finite element model co-simulation[C]//IECON 2013-39th Annual Conference of the IEEE Industrial Electronics Society, Vienna, 2013: 6834-6840.

[13] Parasiliti F, Villani M, Tassi A. Dynamic analysis of synchronous reluctance motor drives based on simulink? and finite element model[C]//IECON 2006-32nd Annual Conference on IEEE Industrial Electronics, Paris, 2006: 1516-1520.

[14] Jabbar M A, Liu Zhejie, Dong Jing. Time-stepping finite-element analysis for the dynamic performance of a permanent magnet synchronous motor[J]. IEEE Transactions on Magnetics, 2003, 39(5): 2621-2623.

[15] Tursini M, Villani M, Tullio A D, et al. Nonlinear model suitable for the offline cosimulation of fault- tolerant PM motors drives[J]. IEEE Transactions on Industry Applications, 2017, 53(4): 3719-3729.

[16] Wang L, Jatskevich J, Dommel H W. Re-examination of synchronous machine modeling techniques for electromagnetic transient simulations[J]. IEEE Transa- ctions on Power Systems, 2007, 22(3): 1221-1230.

[17] Fernandez-Bernal F, Garcia-Cerrada A, Faure R. Determination of parameters in interior permanent- magnet synchronous motors with iron losses without torque measurement[J]. IEEE Transactions on Industry Applications, 2001, 37(5): 1265-1272.

[18] Rafaq M S, Mwasilu F, Kim J, et al. Online parameter identification for model-based sensorless control of interior permanent magnet synchronous machine[J]. IEEE Transactions on Power Electronics, 2017, 32(6): 4631-4643.

[19] Mohammed O A, Liu S, Liu Z. Physical modeling of PM synchronous motors for integrated coupling with machine drives[J]. IEEE Transactions on Magnetics, 2005, 41(5): 1628-1631.

[20] Fasil M, Mijatovic N, Jensen B B, et al. Nonlinear dynamic model of PMBLDC motor considering core losses[J]. IEEE Transactions on Industrial Electronics, 2017, 64(12): 9282-9290.

[21] 魯軍勇, 馬偉明, 孫兆龍, 等. 多段初級直線感應(yīng)電機(jī)靜態(tài)縱向邊端效應(yīng)研究[J]. 中國電機(jī)工程學(xué)報, 2009, 29(33): 95-101.

Lu Junyong, Ma Weiming, Sun Zhaolong, et al. Research on static longitudinal end effect of linear induction motor with multi-segment primary[J]. Proceedings of the CSEE, 2009, 29(33): 95-101.

[22] 聶世雄, 付立軍, 許金, 等. 分段供電直線感應(yīng)電機(jī)動子不對稱模型及參數(shù)計算[J]. 電機(jī)與控制學(xué)報, 2017, 21(2): 10-17.

Nie Shixiong, Fu Lijun, Xu Jin, et al. Asymmetrical model and parameter calculation of segment-powered linear inductive motor mover[J]. Electric Machines and Control, 2017, 21(2): 10-17.

[23] 許金, 馬偉明, 魯軍勇, 等. 分段供電直線感應(yīng)電機(jī)氣隙磁場分布和互感不對稱分析[J]. 中國電機(jī)工程學(xué)報, 2011, 31(15): 61-68.

Xu Jin, Ma Weiming, Lu Junyong, et al. Analysis of air-gap magnetic field distribution and mutual inductance asymmetry of sectionally powered linear induction motor[J]. Proceedings of the CSEE, 2011, 31(15): 61-68.

[24] 牟樹君, 柴建云, 孫旭東, 等. 分段供電交流直線電機(jī)三相互感不對稱分析及抑制[J]. 電工技術(shù)學(xué)報, 2015, 30(1): 81-88.

Mu Shujun, Chai Jianyun, Sun Xudong, et al. Analysis and restrain of mutual inductance asymmetry in the sectionally powered AC linear motor[J]. Transactions of China Electrotechnical Society, 2015, 30(1): 81-88.

[25] 牟樹君, 柴建云, 孫旭東, 等. 分段供電交流直線電機(jī)中偏置磁通密度的分析及其消除方法[J]. 電工技術(shù)學(xué)報, 2014, 29(3): 12-20.

Mu Shujun, Chai Jianyun, Sun Xudong, et al. Analysis and elimination of the off-set flux density in the sectionally powered AC linear motors[J]. Transa- ctions of China Electrotechnical Society, 2014, 29(3): 12-20.

[26] 牟樹君, 柴建云, 孫旭東, 等. 長初級分段供電變極距直線感應(yīng)電機(jī)的氣隙磁場分布及電感參數(shù)特性[J]. 清華大學(xué)學(xué)報: 自然科學(xué)版, 2014, 54(9): 1161-1165, 1172.

Mu Shujun, Chai Jianyun, Sun Xudong, et al. Air-gap magnetic field distribution and inductance chara- cteristic in long primary sectional powered variable pole pitch linear induction motor[J]. Journal of Tsinghua University: Science and Technology, 2014, 54(9): 1161-1165, 1172.

[27] 孫兆龍, 馬偉明, 魯軍勇, 等. 長初級雙邊直線感應(yīng)電動機(jī)靜態(tài)縱向邊端效應(yīng)及阻抗矩陣研究[J]. 中國電機(jī)工程學(xué)報, 2010, 30(18): 72-77.

Sun Zhaolong, Ma Weiming, Lu Junyong, et al. Research of static longitudinal end effect and impedance matrix for long primary double-sided linear induction motors[J]. Proceedings of the CSEE，2010, 30(18): 72-77.

[28] 上官璇峰, 勵慶孚, 袁世鷹. 多段初級永磁直線同步電動機(jī)驅(qū)動系統(tǒng)整體建模和仿真[J]. 電工技術(shù)學(xué)報, 2006, 21(3): 52-57.

Shangguan Xuanfeng, Li Qingfu, Yuan Shiying. Integrated modeling and simulation of the systems driven by multi-segment primary permanent linear synchronous motors[J]. Transactions of China Elec- trotechnical Society, 2006, 21(3): 52-57.

[29] 上官璇峰, 勵慶孚, 袁世鷹. 多段初級永磁直線同步電機(jī)驅(qū)動的垂直提升系統(tǒng)[J]. 中國電機(jī)工程學(xué)報, 2007, 27(18): 7-12.

Shangguan Xuanfeng, Li Qingfu, Yuan Shiying. The vertical hoist system driven by OMLSM with multi- segment primary[J]. Proceedings of the CSEE, 2007, 27(18): 7-12.

[30] 劉凱, 張炳義, 馮桂宏. 基于非對稱繞組函數(shù)法永磁同步電機(jī)偏心電感參數(shù)的研究[J]. 電工技術(shù)學(xué)報, 2020, 35(增刊2): 387-394.

Liu Kai, Zhang Bingyi, Feng Guihong. Research on eccentric inductance of permanent magnet machine based on asymmetrical winding function approach[J]. Transactions of China Electrotechnical Society, 2020, 35(S2): 387-394.

Non-Linear Mathematic Model of a Segmented Powered Permanent Magnet Linear Synchronous Machine

1,21,211,21

（1. Key Laboratory of Power Electronics and Electric Drive Institute of Electrical Engineering Chinese Academy of Sciences Beijing 100190 China 2. University of Chinese Academy of Sciences Beijing 100049 China）

Due to the segmented arrangement of windings in long primary dual-three phase permanent magnet linear synchronous machine (PMLSM), the adjacent segments will cause the off-set pulsating magnetic field, and the mover passing through the powered segment will also cause complex transient processes. In this paper, in order to describe the dynamic characteristics of the PMLSM precisely, a phase-domain mathematic model is established and the two fore-mentioned special effects are considered. The causes of off-set pulsation in the air-gap flux distribution are analyzed by the analytical method, and the unbalanced inductance and saturation effects caused by the off-set flux distribution are studied. A 2-D inductance look-up table inabcoordinate system is established to reduce the computation consumption of the model. Then, the easy-to-derivative Sigmoid function is introduced to describe the no load flux linkage and back electromagnetic force. Finally, the mathematic model is established, the voltage and thrust are calculated by the model and verified by the finite element method. Finally, an experiment platform based on the RT-Lab is built. The results show that the model can reflect the various effects of the motor and can replace the finite element model with a large computational cost.

Permanent magnet linear synchronous machine, segmented powered, saturation effect of magnetic field, phase-domain model

TM359.4

10.19595/j.cnki.1000-6753.tces.201078

國家重點研發(fā)計劃高速磁浮交通系統(tǒng)關(guān)鍵技術(shù)研究課題資助項目（2016YFB1200602-20）。

2020-08-30

2020-12-01

郭科宇 男，1993年生，博士研究生，研究方向為永磁直線同步電機(jī)設(shè)計及控制。E-mail: guoky93@mail.iee.ac.cn

史黎明 男，1964年生，研究員，研究方向為直線電機(jī)設(shè)計及其控制、無線電能傳輸技術(shù)。E-mail: limings@mail.iee.ac.cn（通信作者）

（編輯 陳 誠）