直線型磁通切換電機(jī)特性分析的時(shí)步有限元模型

摘" 要： 針對直線型雙邊磁通切換永磁電機(jī)（DLFSPM）結(jié)構(gòu)相對復(fù)雜，現(xiàn)階段運(yùn)行特性研究不充分的問題，以時(shí)步有限元為工具，建立了一種研究DLFSPM運(yùn)行特性的場路耦合分析模型。闡述了DLFSPM電動(dòng)運(yùn)行原理，并推導(dǎo)出ABC坐標(biāo)系下的運(yùn)行特性方程。為解決直線電機(jī)仿真環(huán)境下的開域處理問題，建立圓弧形電機(jī)仿真拓?fù)洌砸慌_額定輸出電磁推力為400 N的DLFSPM樣機(jī)為例，得到可以同時(shí)考慮DLFSPM內(nèi)部磁場和外電路特性的場路耦合時(shí)步有限元模型。最后，利用該時(shí)步有限元分析模型研究了其空載特性和力角特性，仿真與實(shí)測的相關(guān)對比曲線結(jié)果表明二者吻合較好，說明該模型具有良好的計(jì)算精度，進(jìn)一步驗(yàn)證了其正確性和有效性。

關(guān)鍵詞： 直線型雙邊磁通切換永磁電機(jī)（DLFSPM）； 場路耦合； 運(yùn)行原理； 時(shí)步有限元分析模型； 空載特性； 力角特性

中圖分類號： TN92?34； TM351" " " " " " " " " " "文獻(xiàn)標(biāo)識碼： A" " " " " " " " " " "文章編號： 1004?373X（2024）24?0025?08

Time?step finite element model for operational characteristics analysis of LFSPM

DU Wei1， 2， WU Yulei1， 2， WEN Cheng1， 2， MA Zengqiang1， 2

（1. School of Electrical and Electronic Engineering， Shijiazhuang Tiedao University， Shijiazhuang 050043， China;

2. Hebei Collaborative Innovation Center of Transportation Power Grid Intelligent Integration Technology and Equipment， Shijiazhuang 050043， China）

Abstract： In response to the relatively complex structure and insufficient research on the operating characteristics of double?sided linear flux switched permanent magnet motors （DLFSPM）， a field?circuit coupling analysis model for studying the operating characteristics of DLFSPM is established by using time step finite element method as a tool. The motor?operation principle of DLFSPM is described， and the operating characteristic equation in ABC coordinate system is derived. A simulation topology for circular?arc motors is established to address the issue of open domain processing in the simulation environment of linear motors. Taking a DLFSPM prototype with a rated electromagnetic thrust of 400 N as an example， the field?circuit coupling time?step finite element model is obtained that can simultaneously consider the internal magnetic field and external circuit characteristics of the DLFSPM. The time step finite element analysis model is used to study its no?load characteristics and force angle characteristics. The comparison curve results between simulation and actual measurement show good agreement between the two， indicating that the model has good computational accuracy and further verifying its correctness and effectiveness.

Keywords： double?sided linear flux switched permanent magnet motor; field?circuit coupling; operating principle; time step finite element analysis model; empty load characteristic; force angle characteristic

0" 引" 言

直線型磁通切換永磁電機(jī)（Linear Flux Switching Permanent Magnet， LFSPM）集成了直線感應(yīng)電機(jī)成本低與傳統(tǒng)永磁直線電機(jī)功率及密度高的特點(diǎn)，其在長行程、隔磁條件下的直線運(yùn)動(dòng)領(lǐng)域，如城市軌道交通、數(shù)控機(jī)床系統(tǒng)和物流系統(tǒng)等方面，具有較高的研究價(jià)值和廣闊的應(yīng)用前景[1?3]。

LFSPM包括單邊和雙邊兩種拓?fù)浣Y(jié)構(gòu)，文獻(xiàn)[4?6]列舉了LFSPM的一般設(shè)計(jì)原則。文獻(xiàn)[7?8]基于城市軌道交通牽引系統(tǒng)，分別設(shè)計(jì)了次級分段型和次級齒槽型的單邊LFSPM，并通過與直線感應(yīng)電機(jī)進(jìn)行對比，得出其在推力密度、效率、成本和能耗等方面具有一定的優(yōu)勢。文獻(xiàn)[9]提出了一種雙邊LFSPM，與單邊結(jié)構(gòu)相比，雙邊結(jié)構(gòu)在降低定位力和推力波動(dòng)、減少漏磁及提高推力密度等方面更具優(yōu)勢，但其安裝較為復(fù)雜，應(yīng)用場景較少。文獻(xiàn)[10?11]在dq坐標(biāo)系下建立了LFSPM的數(shù)學(xué)模型，分析了該電機(jī)的速度閉環(huán)控制原理及性能。仿真及試驗(yàn)結(jié)果表明，當(dāng)采用[id=0]的控制方式時(shí)，電機(jī)具有較好的動(dòng)態(tài)性能以及較高的功率因數(shù)和效率。作為一種新型拓?fù)潆姍C(jī)，盡管國內(nèi)外學(xué)者在基礎(chǔ)理論、設(shè)計(jì)方法和控制策略等多個(gè)方面取得了一定的研究成果，但相比于傳統(tǒng)永磁直線電機(jī)，其理論尚不完善，研究的深度和廣度均顯不足。此外，大部分學(xué)者基于LFSPM的靜態(tài)特性[12]開展了研究，但在實(shí)際運(yùn)行中，LFSPM產(chǎn)生氣隙磁場和電磁推力的機(jī)理與傳統(tǒng)永磁直線電機(jī)明顯不同，然而目前鮮有文獻(xiàn)對LFSPM的運(yùn)行特性進(jìn)行深入研究。

鑒于此，本文以低速直線驅(qū)動(dòng)為應(yīng)用背景，以直線型雙邊磁通切換永磁電機(jī)（Double?sided Linear Flux Switching Permanent Magnet， DLFSPM）為研究對象，提出了面向其運(yùn)行特性分析的時(shí)步有限元模型。針對DLFSPM特殊拓?fù)浣Y(jié)構(gòu)，文中推導(dǎo)出abc坐標(biāo)系下運(yùn)行的特性方程，通過建立圓弧形電機(jī)仿真拓?fù)?，解決了直線電機(jī)仿真環(huán)境下的開域處理問題；并以一臺額定推力為400 N的DLFSPM樣機(jī)為例，得出可以同時(shí)考慮DLFSPM內(nèi)部磁場和外電路特性的場路耦合時(shí)步有限元模型。本文利用該模型研究了其空載特性和力角特性；且為驗(yàn)證模型的正確性和有效性，還進(jìn)行了仿真試驗(yàn)，發(fā)現(xiàn)仿真結(jié)果與實(shí)測結(jié)果吻合較優(yōu)，證明了該時(shí)步有限元模型具有良好的計(jì)算精度。

1" DLFSPM電機(jī)電動(dòng)運(yùn)行原理

1.1" 拓?fù)浣Y(jié)構(gòu)

圖1為本文提出的DLFSPM結(jié)構(gòu)及參數(shù)示意圖。表1為樣機(jī)基本參數(shù)。由圖表可知，該電機(jī)包括1個(gè)短初級動(dòng)子和2個(gè)相同的長次級定子，動(dòng)子導(dǎo)磁齒與次級導(dǎo)磁齒均為凸極結(jié)構(gòu)。其中，初級動(dòng)子由2片“H”型磁鋼和5片雙“H”型磁鋼組成。接著，將永磁體嵌于磁鋼之間，采用交替平行充磁形式對6塊永磁體進(jìn)行充磁，而電樞繞組環(huán)繞于導(dǎo)磁齒上。全機(jī)本體結(jié)構(gòu)總計(jì)24個(gè)線槽，線圈按照120°的相位差間隔安裝在上下對稱的線槽中，以構(gòu)成電機(jī)的三相繞組。

1.2" 工作原理

DLFSPM電機(jī)電動(dòng)運(yùn)行原理如圖2所示。動(dòng)子處于位置1時(shí)，初級動(dòng)子鐵芯內(nèi)永磁體與電樞繞組形成合成磁場，m1齒電樞反應(yīng)為去磁，m2齒電樞反應(yīng)為增磁，依據(jù)磁阻最小原理，以初級齒m2的運(yùn)動(dòng)趨勢為主，m2齒向左運(yùn)動(dòng)與次級齒重合對齊；當(dāng)動(dòng)子位置處于位置2時(shí)，m1齒電樞反應(yīng)為增磁，m2齒電樞反應(yīng)為去磁，以初級齒m1的運(yùn)動(dòng)趨勢為主，m1齒向左運(yùn)動(dòng)與次級齒重合對齊?？傊?，通過控制電樞繞組中電流的大小、方向及頻率，即可控制DLFSPM的運(yùn)動(dòng)速度和方向。

2" 運(yùn)行特性方程

2.1" 假設(shè)條件和正方向規(guī)定

本文DLFSPM電機(jī)的運(yùn)行特性指在一定外部條件下，當(dāng)電機(jī)運(yùn)行速度處于穩(wěn)定狀態(tài)時(shí)，其所輸出的電磁推力與力角（功角）間的關(guān)系。為簡化分析模型，做如下假設(shè)：

1） 忽略電機(jī)鐵芯的磁飽和；

2） 忽略不計(jì)鐵芯中的磁滯和渦流損耗；

3） 電樞繞組基波電勢按正弦規(guī)律變化，不考慮電機(jī)磁場中的高次諧波影響；

4） 各相繞組參數(shù)與結(jié)構(gòu)相同且對稱。

為便于分析，規(guī)定A相繞組匝鏈永磁勵(lì)磁磁鏈為最大值時(shí)的定動(dòng)子相對位置為直d軸，具體見圖3a）；順著運(yùn)動(dòng)方向，超前d軸90°（電角度）的位置為交q軸，當(dāng)動(dòng)子運(yùn)動(dòng)到此位置時(shí)，永磁體磁鏈被定子鐵芯短路，如圖3b）所示。正方向的規(guī)定如下：

1） 電樞繞組電壓、電流方向采用電動(dòng)機(jī)慣例；

2） 動(dòng)子向左運(yùn)動(dòng)為正方向；

3） 在相量分析時(shí)，d軸逆時(shí)針旋轉(zhuǎn)為正，且q軸超前d軸90°電角度。

2.2" 力角特性方程

DLFSPM電機(jī)電動(dòng)運(yùn)行的具體分析如圖4所示。圖中，[U]為電樞繞組端電壓；[I]為電樞電流；[R]為電樞電阻；[Xσ]為電樞端部漏抗；[Xs]為同步電抗，包括電樞反應(yīng)電抗和電樞漏電抗（端部漏抗除外）；[E0]為空載反電勢；[Eδ]為氣隙內(nèi)電勢；電樞電壓超前電樞電流的角度（外功率因數(shù)角）為[φ]；電樞電壓超前空載反電勢的角度（力角）為[δ]；電樞電流超前d軸的角度為[λ]。

DLFSPM電機(jī)電壓平衡方程為：

[U=E0+jIXs+jIXσ+IR =Eδ+jIXσ+IR]" （1）

考慮到DLFSPM電機(jī)采用了集中繞組，因此可忽略其端部漏抗，則：

[U=E0+jIXs+IR =Eδ+IR]" " （2）

選[E0]作為參考相量，根據(jù)相量圖可得：

[U=Ucosδ+jsinδ] " （3）

由式（2）和式（3）可得：

[IP=Ucosδ-E0R+XsUsinδR2+X2sIQ=RUsinδ-Ucosδ-E0XsR2+X2s]" "（4）

式中：[IP]為電樞電流的有功分量；[IQ]為電樞電流的無功分量。

DLFSPM電機(jī)氣隙電磁功率為：

[Pm=mE0IP =mE0Ucosδ-E20R+mXsE0UsinδR2+X2s]" （5）

當(dāng)[R]遠(yuǎn)小于[Xs]時(shí)，可忽略[R]的影響，則：

[Pm=mE0IP=mE0UsinδXs] " "（6）

如果不考慮次級定子的渦流損耗，則[P輸出≈Pm]。因此，電磁推力可以表示為：

[F輸出=P輸出v=P輸出τPf=mE0UsinδvXs]" （7）

式中：m為電機(jī)相數(shù)；v為動(dòng)子機(jī)械線速度；[P輸出]為實(shí)際輸出的機(jī)械功率；f為電源頻率；[τP]為次級極距。

此外，采用負(fù)載角[λ]表達(dá)氣隙功率：

[Pm=mE0IP=mE0Icos90°-λ]" （8）

則DLFSPM電機(jī)輸出的電磁推力可表示為：

[F輸出=P輸出v=mE0Icos90°-λv]" （9）

在電機(jī)帶載運(yùn)行時(shí)，由于電樞反應(yīng)使得永磁體工作點(diǎn)發(fā)生了變化，其電樞鐵芯的非線性也會導(dǎo)致電樞反應(yīng)電抗發(fā)生變化。因此，為獲得更為精確的電機(jī)運(yùn)行特性，本文借助有限元法展開研究。

3" 時(shí)步有限元分析模型

為計(jì)及電樞反應(yīng)電抗的變化及氣隙諧波的影響，本文以時(shí)步有限元法[13]為工具，建立了面向DLFSPM電機(jī)運(yùn)行特性分析的場路耦合模型。

3.1" 電磁場有限元方程

考慮到DLFSPM電機(jī)定子、動(dòng)子的實(shí)際結(jié)構(gòu)，在對電機(jī)內(nèi)磁場進(jìn)行分析時(shí)，做如下假設(shè)：

1） 磁場沿電機(jī)縱向分布不變，將磁場作為二維場處理，并采用直角坐標(biāo)系進(jìn)行計(jì)算；

2） 電機(jī)鐵芯沖片材料為各向同性，忽略磁滯的影響，將任一瞬間的磁場作為非線性恒定磁場處理；

3） 動(dòng)子電樞采用了集中繞組，可忽略各繞組的端部漏抗。

綜上，在求解區(qū)域內(nèi)可獲得電磁場方程的邊值形式：

[Ω：??x1μ?A?x+??y1μ?A?y=-Js+JpmΓ：A=0] （10）

式中：[A]為矢量磁位；[μ]為磁導(dǎo)率；[Js]為電樞電流密度；[Jpm]為永磁體的等效面電流密度；[Γ]為第一類邊界條件。

用加權(quán)余量法將電磁場方程離散化，其有限元方程為：

[KA = CI+Q] （11）

式中：[K]為剛度矩陣；[C]為電樞電流激勵(lì)的關(guān)聯(lián)矩陣；[I]為動(dòng)子電樞電流，[I=iA，iB，iCT]；[Q]為永磁體激勵(lì)列向量。其中，[K]、[C]、[Q]矩陣各元素的表達(dá)式參見文獻(xiàn)[6]。

3.2" 電樞繞組回路方程

電樞繞組采用集中參數(shù)形式的電路方程，并根據(jù)圖4將電樞繞組回路方程表達(dá)為：

[Us=Eδ+IRs] （12）

式中：[Us]為電樞繞組端電壓，[Us=uA，uB，uCT]；[Eδ]為氣隙電勢，[Eδ=eA，eB，eCT]；[Rs]為電樞繞組電阻，[Rs=diagra，rb，rcT]。將[Eδ=lef?CT?A]代入式（12）可得：

[Us=lef?CT?A+IRs] （13）

3.3" 方程的耦合與求解

將磁場有限元方程式（11）與外電路方程式（13）耦合到同一個(gè)方程中，能夠得到DLFSPM電機(jī)場路耦合時(shí)步有限元方程：

[K-Clef?CTRsAI=QUs] " （14）

式中，[A]、[I]、[Us]均為時(shí)間的函數(shù)，可采用后差分歐拉法對其進(jìn)行時(shí)間離散，計(jì)算電樞反應(yīng)電抗變化及氣隙諧波后，再采用牛頓?拉夫遜法對該非線性微分方程進(jìn)行求解。

3.4" 仿真環(huán)境下開域問題的處理

與旋轉(zhuǎn)電機(jī)不同，直線電機(jī)的運(yùn)行路徑具有無循環(huán)特性，其動(dòng)子鐵芯的開斷會使得氣隙處于開域環(huán)境，而動(dòng)子移動(dòng)則會帶來網(wǎng)格重剖分問題，降低計(jì)算效率。為解決該問題，本文建立了DLFSPM電機(jī)的圓弧形運(yùn)行路徑，具體如圖5所示。此次采用移動(dòng)網(wǎng)格的分析方法[14]，對動(dòng)子移動(dòng)時(shí)的網(wǎng)格剖分問題進(jìn)行分析。

DLFSPM電機(jī)有兩個(gè)定子位于動(dòng)子兩側(cè)，且在動(dòng)子兩側(cè)均存在有氣隙。為了便于表述，本文將其稱為內(nèi)定子側(cè)氣隙和外定子側(cè)氣隙，如圖5a）所示。以動(dòng)子向左移動(dòng)即逆時(shí)針旋轉(zhuǎn)為例說明具體的處理方法。

1） 為減少剖分工作量及網(wǎng)格畸變對計(jì)算結(jié)果的影響，將氣隙至少分為3層，中間層為一個(gè)規(guī)則的圓環(huán)，圓環(huán)靠近動(dòng)子側(cè)的節(jié)點(diǎn)屬于動(dòng)子區(qū)域，靠近定子側(cè)的節(jié)點(diǎn)屬于定子區(qū)域。此外，圓環(huán)兩側(cè)節(jié)點(diǎn)數(shù)量需相等，以保證節(jié)點(diǎn)旋轉(zhuǎn)時(shí)的網(wǎng)格質(zhì)量。

2） DLFSPM電機(jī)動(dòng)子初始位置如圖5b）所示；首先將全部動(dòng)子節(jié)點(diǎn)、內(nèi)定子側(cè)氣隙節(jié)點(diǎn)和內(nèi)定子節(jié)點(diǎn)逆時(shí)針旋轉(zhuǎn)一個(gè)角度，結(jié)果如圖5c）所示；然后將內(nèi)定子區(qū)域順時(shí)針旋轉(zhuǎn)相同角度，結(jié)果如圖5d）所示，從而完成動(dòng)子的前移。

該方法在處理動(dòng)子移動(dòng)時(shí)，只需處理氣隙中間層的網(wǎng)格連接，無需重復(fù)剖分網(wǎng)格，因此大幅提高了仿真效率。

4" 運(yùn)行特性研究

4.1" 空載特性

DLFSPM電機(jī)的一個(gè)重要指標(biāo)就是空載反電勢。分析其空載反電勢的波形和諧波分布，對于驗(yàn)證電機(jī)設(shè)計(jì)的合理性和判斷有限元基本解的正確性具有重要意義。

以永磁體作為單獨(dú)激勵(lì)源，電樞電流為零，繞組處于開路狀態(tài)。在電動(dòng)機(jī)慣例下，A相繞組的空載反電勢為：

[eA=dΨAdt=ddtN1lefaSbe=1NaaΔeAi+Aj+Ak3] （15）

式中：[ΨA]為A相繞組永磁磁鏈；[N1]為A相繞組串聯(lián)匝數(shù)；[lef]為電機(jī)縱向深度；[aa]為A相繞組并聯(lián)支路數(shù)；[N]為繞組單元總數(shù)；[Δe]為單元面積；[Sb]為A相繞組面積。當(dāng)單元屬于A相繞組流出區(qū)域時(shí)，[aa=1]；當(dāng)單元屬于A相繞組流入?yún)^(qū)域時(shí)，[aa=-1]；當(dāng)單元不屬于A相繞組區(qū)域時(shí)，[aa=0]。

保存瞬態(tài)計(jì)算每個(gè)時(shí)步下的繞組端電壓，可以得到空載反電勢波形，具體如圖6所示。

計(jì)算得到DLFSPM電機(jī)的空載反電勢幅值分別為E01=22.47 V，E03=0.57 V，E05=0.32 V。其中，反電勢最大諧波分量分別為3次和5次諧波，諧波畸變率（THD）較小，為2.9%。由此說明所設(shè)計(jì)電機(jī)的反電動(dòng)勢正弦度適用于BLAC控制。

4.2" 力角特性

本文所提出的場路耦合時(shí)步有限元模型具有較優(yōu)的通用性，基于電樞繞組電阻及電流或電壓幅值與相位信息，即可進(jìn)行運(yùn)行特性的仿真驗(yàn)證。此外，考慮電機(jī)運(yùn)行的穩(wěn)定性，還計(jì)算了單側(cè)法向吸力與力角的關(guān)系，計(jì)算流程如下。

1） 給定A相電樞電流的幅值和相位[λ]，從而確定三相電流的瞬時(shí)值。

2） 根據(jù)時(shí)步有限元法計(jì)算電機(jī)在動(dòng)子移動(dòng)過程中不同時(shí)刻的二維恒定非線性磁場。

3） 根據(jù)磁場計(jì)算結(jié)果，計(jì)算電磁推力、法向力及力角。

4） 保存每個(gè)電樞電流相位[λ]下的輸出電磁推力、法向力與力角，得到穩(wěn)態(tài)運(yùn)行特性。

電磁推力和法向力的計(jì)算采用文獻(xiàn)[15]提出的改進(jìn)麥克斯韋應(yīng)力法，電磁推力[Fe]和法向力[Fn]的具體計(jì)算公式分別為：

[Fe=lefμ0δe=1NBτBnΔe] （16）

[Fn=lef2μ0δe=1NB2n-B2τΔe] （17）

式中：[μ0]為氣隙磁導(dǎo)率；[Bτ]為切向磁密；[Bn]為法向磁密。

按照上述計(jì)算流程，可以得到DLFSPM電機(jī)力角特性，具體如圖7所示。

圖7a）給出了電磁推力隨力角變化的曲線。由圖可以看出：在額定電樞電流條件下，該DLFSPM電機(jī)最大輸出電磁推力為680 N，對應(yīng)的功角約為85°（B點(diǎn)位置），對應(yīng)的輸入電流相位[λ]為95°；額定電磁推力為400 N，對應(yīng)的功角約為35°（A點(diǎn)位置），對應(yīng)的輸入電流相位[λ]為40°，最大過載倍數(shù)約為1.7。

圖7b）給出了單側(cè)法向力隨力角變化的曲線?？梢钥闯?，隨著力角從0°～180°逐漸增大，法向力為減小的趨勢。在初始位置δ=0°時(shí)，為最大法向力4 285 N，在δ=180°時(shí)，為最小法向力2 788 N。圖中C點(diǎn)為額定電磁推力時(shí)的法向力，為4 156 N，單側(cè)法向力約為電磁推力的10倍；圖中D點(diǎn)為最大電磁推力時(shí)的法向力，為3 569 N，此時(shí)法向力約為電磁推力的5倍。

可借助相量圖對力角特性進(jìn)行原理分析。由式（7）和式（9）可知：在忽略電樞電阻和鐵芯磁飽和的假設(shè)下，當(dāng)δ=90°即λ=90°時(shí)，電磁推力達(dá)到最大；考慮電樞電阻和鐵芯磁飽和后，經(jīng)有限元計(jì)算可知，δ=85°即λ=95°時(shí)，出現(xiàn)最大電磁力，與相量分析吻合。

法向力隨[δ]的變化趨勢可以借助圖8進(jìn)行分析。圖中，[ψ0]為永磁體勵(lì)磁磁場，[ψs]為電樞反應(yīng)磁場，[ψδ]為氣隙合成磁場。在忽略電樞電阻和鐵芯磁飽和的假設(shè)下，可以認(rèn)為[ψ0]和[ψs]保持幅值不變，電流相位[λ]即為[ψ0]和[ψs]之間的夾角，隨著[λ]在0°～180°變化，此時(shí)力角[δ]也在0°～180°變化，氣隙合成磁場[ψδ]呈現(xiàn)減小的趨勢。而對于本文研究的DLFSPM電機(jī)而言，其法向磁場分量要比切向磁場分量大得多，因此，其法向力必然會隨之呈現(xiàn)減小的趨勢。

5" 試驗(yàn)驗(yàn)證

DLFSPM樣機(jī)實(shí)物圖如圖9所示。限于試驗(yàn)條件，本文所研究的DLFSPM電機(jī)尚未做全面測試，為了驗(yàn)證所提方法的可靠性，從方便對比又不失一般性角度考慮，測試了該樣機(jī)在反拖空載運(yùn)行時(shí)的感應(yīng)電勢波形。圖10a）和圖10b）分別為樣機(jī)的實(shí)測和仿真空載感應(yīng)電勢波形。

限于動(dòng)子的運(yùn)動(dòng)行程，本文在動(dòng)子運(yùn)動(dòng)速度為0.076 m/s的情況下進(jìn)行了空載感應(yīng)電勢的測試，通過比較可以看出，仿真與試驗(yàn)結(jié)果較為相似，二者均具有較好的正弦特性。其中，有限元仿真得到的感應(yīng)電勢幅值為EA=3.12 V，EB=3.13 V，EC=3.18 V，實(shí)測得到的感應(yīng)電勢幅值為EA=3.21 V，比仿真結(jié)果大2.88%；EB=3.25 V，比仿真結(jié)果大3.83%；EC=3.19 V，比仿真結(jié)果大0.31%。兩者之間的誤差主要是由電機(jī)端部效應(yīng)和測量誤差所引起的。

空載感應(yīng)電勢仿真與實(shí)測結(jié)果的對比證明了有限元基本解的正確性和有效性，驗(yàn)證了所提模型的可靠性。

6" 結(jié)" 論

針對直線型磁通切換永磁電機(jī)運(yùn)行特性研究不充分的問題，本文提出了一種面向DLFSPM電機(jī)運(yùn)行特性分析的時(shí)步有限元模型。該模型將外電路方程和磁場有限元方程聯(lián)立求解，可得出電磁推力、法向力隨力角的變化規(guī)律。通過對比仿真與實(shí)測的相關(guān)曲線，證明了該方法的有效性。主要結(jié)論如下。

1） 通過推導(dǎo)DLFSPM電機(jī)運(yùn)行特性方程，得到了力角[δ]和負(fù)載角[λ]表征的電磁推力表達(dá)式。考慮到DLFSPM電機(jī)結(jié)構(gòu)的特殊性，文中還建立了面向其運(yùn)行特性分析的時(shí)步有限元模型，針對直線電機(jī)運(yùn)行路徑的無循環(huán)特性，采用圓弧形運(yùn)行路徑和移動(dòng)網(wǎng)格剖分技術(shù)，避免了動(dòng)子移動(dòng)時(shí)的重復(fù)剖分。

2） 使用所提模型計(jì)算得到了樣機(jī)的空載感應(yīng)電勢波形，通過研究該樣機(jī)在額定電流條件下的力角特性曲線，得到了其過載能力。

3） 為驗(yàn)證該模型的有效性，以空載感應(yīng)電勢響應(yīng)為標(biāo)準(zhǔn)，將上述模型的計(jì)算結(jié)果和試驗(yàn)測量的結(jié)果進(jìn)行對比。由對比結(jié)果可知，兩者的數(shù)據(jù)高度吻合，能夠證明本文模型的可靠性。

參考文獻(xiàn)

[1] 沈燚明，盧琴芬.初級勵(lì)磁型永磁直線電機(jī)研究現(xiàn)狀與展望[J].電工技術(shù)學(xué)報(bào)，2021，36（11）：2325?2343.

[2] NOMAN U， FAISAL K， ABDUL B， et al. Analytical airgap field model and experimental validation of double sided hybrid excited linear flux switching machine [J]. IEEE access， 2021（9）： 117120?117131.

[3] SHAHID H， FAISAL K， WASIQ U， et al. Development of a low?cost modular structure fault tolerant field excited flux switching linear machine for urban rail transit [J]. IEEE access， 2021（9）： 165854?165864.

[4] HAO C， AYMAN M E， NABEEL A O. Flux?switching permanent magnet machines： a review of opportunities and challenges?part I： fundamentals and topologies [J]. IEEE transactions on energy conversion， 2020， 35（2）： 684?698.

[5] YUE Y， JIA S F， LIANG D L. New topologies of high torque density machine based on magnetic field modulation principle [J]. CES transactions on electrical machines and systems， 2023， 7（1）： 1?10.

[6] 程明，花為.定子永磁無刷電機(jī)理論、設(shè)計(jì)與控制[M].北京：科學(xué)出版社，2021：35?56.

[7] 曹瑞武，蘇恩超，張學(xué).軌道交通用次級分段型直線磁通切換永磁電機(jī)研究[J].電工技術(shù)學(xué)報(bào)，2020，35（5）：1001?1012.

[8] CAO R， LU M， JIANG N， et al. Comparison between linear induction motor and linear flux?switching permanent?magnet motor for railway transportation [J]. IEEE transactions on industrial electronics， 2019， 66（12）： 9394?9405.

[9] WEN C， ZHAO Q K， LI M W， et al. Multi?objective optimization based on hyper? parameter random forest regression. for linear motor design [J]. International journal of machine learning and cybernetics， 2022， 13（10）： 2929?2942.

[10] CAO R W， CHENG M， MI C， et al. Modeling of a complementary and modular linear fluxswitching permanent magnet motor [J]. IEEE transactions on energy conversion， 2012， 27（2）： 489?497.

[11] CAO R W， CHENG M， ZHANG B F. Speed control of complementary and modular linear flux?switching permanent magnet motor [J]. IEEE transactions on industrial electronics， 2015， 62（7）： 4056?4064.

[12] 姚新元.初級錯(cuò)位雙邊直線磁通切換永磁電機(jī)控制系統(tǒng)研究[D].鎮(zhèn)江：江蘇大學(xué)，2022.

[13] 楊曉霞，王翠珍，張芳，等.計(jì)及勵(lì)磁系統(tǒng)的同步電機(jī)時(shí)步有限元模型研究[J].大電機(jī)技術(shù)，2023（1）：21?29.

[14] 周琴.二維拋物型奇異攝動(dòng)問題的移動(dòng)網(wǎng)格方法[J].工程數(shù)學(xué)學(xué)報(bào)，2021，38（6）：869?878.

[15] 康錦萍，李志強(qiáng)，劉曉芳，等.汽輪發(fā)電機(jī)電磁轉(zhuǎn)矩計(jì)算方法的對比[J].電工技術(shù)學(xué)報(bào)，2009，24（9）：38?43.

作者簡介：杜" ?。?984—），男，河北保定人，博士研究生，講師，研究方向?yàn)橹本€電機(jī)及其控制。

吳玉磊（1997—），男，河北唐山人，在讀碩士研究生，研究方向?yàn)榇磐ㄇ袚Q永磁直線電機(jī)設(shè)計(jì)及其控制。

聞" 程（1985—），女，湖北武漢人，博士研究生，副教授，研究方向?yàn)橹本€電機(jī)及其控制。

馬增強(qiáng)（1975—），男，河北石家莊人，博士研究生，教授，研究方向?yàn)殡姎庠O(shè)備的在線檢測與故障診斷。