做直線電機(jī)設(shè)計繞不開有限元仿真這一步。而 Maxwell仿真直線電機(jī)最核心、最常碰到的就是空載反電動勢、推力輸出和磁阻力這三類分析?？蛰d反電動勢反映磁場設(shè)計和繞組配合的水平推力輸出決定電機(jī)能不能達(dá)到額定出力要求磁阻力由端部力和齒槽力構(gòu)成則直接決定了推力波動和定位精度。這三塊跑通了一臺平板型直線電機(jī)的基本電磁性能就完全清楚了。這篇文章的目標(biāo)讀者有兩類一類是剛接觸直線電機(jī)電磁仿真模型能搭出來但后面不知道看什么結(jié)果的工程師和研究生另一類是已經(jīng)跑過一些仿真但波形讀得不夠透、力曲線里的波動成分分不清來源的人。文章按“先想清楚仿什么再動手建模最后把結(jié)果讀透”的順序來寫把平板型直線電機(jī)的三種典型仿真完整走一遍同時把我在實際調(diào)試中踩過的坑一并說明。1. 整體設(shè)計與仿真思路拆解1.1 三類仿真分別回答什么問題直線電機(jī)本質(zhì)上是一臺展開的旋轉(zhuǎn)電機(jī)。把旋轉(zhuǎn)電機(jī)的定子沿圓周展開旋轉(zhuǎn)磁場變成了沿直線運動的行波磁場推力取代轉(zhuǎn)矩成為核心輸出。但展開帶來的不只是結(jié)構(gòu)變化電磁特性也跟著改變了旋轉(zhuǎn)電機(jī)磁路首尾閉合完全對稱直線電機(jī)有頭有尾鐵心兩端磁路被截斷于是多出了旋轉(zhuǎn)電機(jī)沒有的端部效應(yīng)。因此直線電機(jī)的電磁仿真天然要回答三個層次的問題磁場設(shè)計對不對空載反電動勢的波形、幅值、對稱性直接反映永磁體極對數(shù)、繞組極相組配合、氣隙尺寸這些參數(shù)選得合不合理。出力夠不夠通電后能產(chǎn)生多大的穩(wěn)定推力推力波動有多大這決定了電機(jī)帶負(fù)載能力和運動平穩(wěn)性。波動從哪里來磁阻力也就是定位力由端部力和齒槽力兩部分組成。端部力來自鐵心長度有限齒槽力來自永磁體與齒槽的相互吸引必須分開看清楚才能針對性地去抑制。這三類仿真不是并列關(guān)系而是遞進(jìn)關(guān)系先確認(rèn)磁場沒問題再看出力表現(xiàn)最后分析波動來源并驗證抑制方案。1.2 仿真順序與模型復(fù)用邏輯實際項目里我建議的順序是先做空載反電動勢仿真再做零電流磁阻力仿真最后做負(fù)載推力仿真。空載和磁阻力仿真的共同點是繞組不用通電流模型幾乎完全復(fù)用區(qū)別只在運動設(shè)置和后處理提取角度。負(fù)載推力仿真是在空載模型基礎(chǔ)上把電流激勵加進(jìn)去同時還能把磁阻力影響疊加進(jìn)來觀察。這套流程最大的價值在模型復(fù)用。模型只建一次通過切換激勵方式、運動方式和后處理表達(dá)式就能覆蓋三種仿真場景。所以在建模階段就要為復(fù)用留好余地運動域和邊界條件按統(tǒng)一規(guī)則設(shè)定繞組激勵用電流源形式這樣后面切換工況只需要改激勵值不需要動幾何結(jié)構(gòu)。1.3 工具選型解析Maxwell瞬態(tài)場的適用邊界這一類仿真普遍選擇Maxwell的瞬態(tài)場求解器而不是靜態(tài)求解器關(guān)鍵原因在空載反電動勢的本質(zhì)。它是一個隨時間變化的感應(yīng)電壓必須讓動子實際運動起來在時間軸上采樣才能得到。瞬態(tài)場求解器天然支持物體運動與網(wǎng)格重構(gòu)適合這種“動起來看結(jié)果”的場景。磁阻力的仿真也同理。磁阻力雖然與電流無關(guān)但它隨動子位置周期性變化必須在多個位置采樣才能畫出完整的力-位移曲線。用靜態(tài)場掃描位置也能做到但每個位置要重新剖分網(wǎng)格、重新求解效率太低。瞬態(tài)求解讓動子以恒定速度滑過一段距離一次性把所有位置的力都算出來而且同一段結(jié)果里還能同時提取反電動勢和磁阻力。有一點需要明確二維模型是平板型直線電機(jī)仿真的主力。平板結(jié)構(gòu)沿寬度方向基本不變用縱向剖面做二維求解精度足夠且計算量小。只有電機(jī)寬度較窄、橫向端部漏磁明顯或者需要評估疊片端部渦流時才值得上三維。二維模型的力輸出默認(rèn)按單位深度計算如果設(shè)置了疊厚軟件會自動乘上這個細(xì)節(jié)后面再展開。2. 平板型直線電機(jī)建模與邊界條件配置2.1 平板型結(jié)構(gòu)的兩種拓?fù)湎确智迤桨逍椭本€電機(jī)有兩種常見拓?fù)浣Ｇ耙欢ㄒ确智?。第一種是動圈式動子是線圈繞組定子是鋪在軌道上的永磁體陣列常見于短行程高動態(tài)平臺。第二種是動鐵式動子是帶齒槽的鐵心加繞組定子同樣是永磁體陣列。兩者電磁本質(zhì)相同但運動域、網(wǎng)格和力的提取差別很明顯動鐵式要把整個鐵心連同繞組包進(jìn)運動域動圈式則只包線圈。兩種拓?fù)湓贛axwell里的處理差別主要落在運動域覆蓋范圍上。動圈式運動體幾乎不含鐵磁材料運動域內(nèi)磁場梯度相對平緩網(wǎng)格要求寬松動鐵式運動體包含硅鋼片鐵心運動域要跟隨齒槽的幾何特征網(wǎng)格必須精細(xì)很多。在動手建模前確定拓?fù)淠鼙苊夂竺娲罅糠倒び绕涫沁\動域尺寸和網(wǎng)格剖分策略的選擇。2.2 幾何建模與模型簡化平板型直線電機(jī)的幾何建模我通常直接在Maxwell的建模環(huán)境里畫。結(jié)構(gòu)本身不復(fù)雜鐵心、齒、繞組、永磁體、背鐵都是一個個矩形塊用坐標(biāo)和布爾運算就能搭出來。模型簡化有幾個原則要記住第一倒角、圓角、螺絲孔這類細(xì)節(jié)一律去掉它們對電磁結(jié)果影響極小卻會明顯增加網(wǎng)格剖分負(fù)擔(dān)。第二端部繞組在二維模型里簡化成槽內(nèi)導(dǎo)體的截面不需要畫出繞組實際走向因為二維模型里本來也畫不出。第三如果電機(jī)極對數(shù)較多先建包含完整端部的整機(jī)模型確保端部效應(yīng)被真實還原不要直接用半模型除非你清楚周期邊界條件的適用范圍。特別強(qiáng)調(diào)一點直線電機(jī)的端部效應(yīng)是固有特性模型必須保留真實兩端結(jié)構(gòu)。這與旋轉(zhuǎn)電機(jī)用周期邊界節(jié)約計算量不同直線電機(jī)的邊界是“真實的截斷”不能為了簡化人為把端部抹掉。2.3 運動域、邊界條件與網(wǎng)格處理的三個關(guān)鍵點運動域的設(shè)置直接決定仿真能不能跑起來。運動域需要完全包住運動體同時留足間隙避免運動過程中運動體與靜止部件干涉。運動域形狀建議與運動體外形一致或略大這樣運動域內(nèi)部空氣網(wǎng)格才能跟隨運動體一起移動氣隙區(qū)域的網(wǎng)格不會被反復(fù)重構(gòu)。邊界條件方面平板型直線電機(jī)二維模型的外部輪廓一般設(shè)置自然邊界或第一類邊界條件。氣隙和運動域之間不需要額外設(shè)邊界關(guān)鍵是保證運動域外的靜止氣隙也參與求解不能出現(xiàn)重疊或縫隙。網(wǎng)格剖分的核心就三個氣隙網(wǎng)格要加密沿厚度方向至少兩層以上齒頂附近和永磁體邊緣要局部細(xì)化這些位置磁場梯度大運動域內(nèi)部空氣網(wǎng)格應(yīng)與氣隙網(wǎng)格在交界處保持尺寸一致。如果交界面網(wǎng)格尺寸差異過大運動域滑移邊界上的插值誤差會被放大直接體現(xiàn)在反電動勢波形的光滑度和力曲線的噪聲水平上。3. 空載反電動勢仿真實操3.1 激勵與運動設(shè)置讓電機(jī)空載滑行空載反電動勢仿真的目標(biāo)是考察永磁體磁場在繞組中產(chǎn)生的感應(yīng)電動勢因此繞組不需要通電流。在Maxwell瞬態(tài)場里把三相繞組的激勵都設(shè)成零安培電流源即可相當(dāng)于開路狀態(tài)。這里用零電流源而不用零電壓源是為了避免電壓源帶來的額外電氣方程仿真過程更穩(wěn)定。運動設(shè)置是這一步的關(guān)鍵。空載反電動勢只與速度和磁鏈有關(guān)不需要考慮動力學(xué)所以把運動模式設(shè)為恒速運動而不是帶負(fù)載的機(jī)械瞬態(tài)。速度值取額定速度或任意合適速度都可以但要注意反電動勢幅值隨速度線性變化讀取波形后要換算成單位速度下的反電動勢常數(shù)也就是伏秒每米的量綱這樣不同速度下的結(jié)果才有可比性也才能客觀評估磁場設(shè)計水平。3.2 反電動勢的提取與波形解讀反電動勢的讀取方式主要有兩種。第一種直接查看繞組激勵的感應(yīng)電壓結(jié)果Maxwell后處理里能列出每相繞組的感應(yīng)電動勢隨時間變化的波形第二種是提取每相磁鏈然后對時間求導(dǎo)得到反電動勢。兩種方法各有適用場景感應(yīng)電壓結(jié)果方便快捷但時間步長偏大時波形容易出毛刺磁鏈求導(dǎo)更光滑適合對波形細(xì)節(jié)要求較高的分析。拿到三相反電動勢波形后重點看四件事幅值是否一致、相位是否互差120度、波形正弦性好不好、反電動勢系數(shù)大小是否落在設(shè)計預(yù)期。反電動勢幅值偏小說明永磁體用量不足或氣隙過大波形畸變明顯說明繞組極相配合或齒形設(shè)計有問題三相對稱性差要檢查繞組接線和幾何建模是否出錯。我習(xí)慣順手對波形做一次快速傅里葉分析看高次諧波含量。反電動勢的諧波直接決定推力波動的諧波成分這個數(shù)據(jù)在后面的推力仿真里很有參考價值。3.3 常見波形異常與調(diào)整方向空載反電動勢仿真最容易暴露建模問題。波形不對稱、幅值忽大忽小、相位錯位這些問題在空載仿真里通常一眼就能發(fā)現(xiàn)。如果三相反電動勢不對稱優(yōu)先檢查繞組相序和槽內(nèi)導(dǎo)體的激勵方向。Maxwell里每個繞組需要設(shè)置正確的激勵方向某個槽的導(dǎo)體電流方向反了該相反電動勢波形就會明顯異常。如果三相都對稱但幅值與理論值差很多檢查永磁體充磁方向是否設(shè)置正確以及氣隙尺寸與實際設(shè)計是否一致。波形毛刺多一般不是物理問題而是數(shù)值問題。要么時間步長太大要么氣隙網(wǎng)格不夠密。把每個電周期的采樣點數(shù)增加到一百個以上同時加密氣隙網(wǎng)格毛刺通常大幅改善。4. 推力輸出仿真實操4.1 負(fù)載工況配置電流激勵與電流角推力輸出仿真是在空載模型基礎(chǔ)上加電流激勵。三相繞組通入幅值和頻率按工況設(shè)定的正弦電流三相互差120度相位。電流幅值取額定電流頻率和速度要匹配電機(jī)工作點直線電機(jī)運行速度與電頻率的關(guān)系由極距決定電頻率等于速度除以兩倍極距。這個關(guān)系要在激勵設(shè)置里算準(zhǔn)確否則電流行波速度與動子運動速度不匹配推力會明顯下降。電流角的選擇同樣影響推力大小。對表貼式永磁直線電機(jī)通常采用直軸電流為零的控制方式也就是電流矢量與反電動勢矢量同相位此時單位電流產(chǎn)生的推力最大。如果電機(jī)有凸極效應(yīng)則需要掃描電流角找最大推力點這屬于后期優(yōu)化內(nèi)容初次驗證額定推力時直軸電流為零的假設(shè)已經(jīng)夠用。4.2 推力結(jié)果提取與推力波動系數(shù)推力的提取很簡單。瞬態(tài)仿真結(jié)束后選中所有運動部件查看力分量隨時間的變化曲線。運動方向如果是橫軸方向就提取對應(yīng)的分力。注意不要只選鐵心要把運動域內(nèi)的所有物體包括鐵心、繞組和隨動的永磁體全部選中作為一個整體提取合力否則得到的是不完整的分力。推力波形最有價值的信息有兩個平均推力和推力波動。平均推力要跟設(shè)計值對照誤差在幾個百分點內(nèi)正常。推力波動用推力波動系數(shù)量化定義有多種常見的是最大值減最小值除以兩倍平均推力也有按峰峰值除以平均值的寫法寫報告時注明定義即可。推力波形的形狀同樣值得觀察。正常情況下的推力由恒定分量加周期性波動組成波動頻率與電頻率存在固定倍數(shù)關(guān)系。如果出現(xiàn)緩慢變化的非周期分量多半是仿真還沒達(dá)到穩(wěn)態(tài)需要繼續(xù)跑更長時間讓磁鏈建立過程的影響消退。4.3 推力曲線背后的電機(jī)設(shè)計信息推力波動大不代表電機(jī)一定有問題關(guān)鍵在于分清波動來源。推力波動至少有三個來源電流產(chǎn)生的電磁推力自身的諧波、磁阻力中的齒槽力、以及端部力。齒槽力周期與齒距相關(guān)端部力周期與極距相關(guān)三者疊加在一起才是總推力波動。實際分析時可以做一次減法實驗把負(fù)載工況的總推力減去空載時的磁阻力剩下的就是純電磁推力部分。如果減完之后波動基本消失說明推力波動主要是磁阻力貢獻(xiàn)下一步該優(yōu)化端部結(jié)構(gòu)和齒槽形狀如果減完之后波動仍然很大問題就出在電流與反電動勢的配合上比如反電動勢諧波與電流諧波相互作用產(chǎn)生額外脈動。這一步分析對定位設(shè)計問題非常有效建議每次推力仿真后都做一遍。5. 磁阻力、端部力與齒槽力仿真分離5.1 三種力的物理來源先分清再仿真直線電機(jī)的磁阻力在不同文獻(xiàn)里有不同叫法定位力、齒槽效應(yīng)力都有工程上習(xí)慣按來源進(jìn)一步拆成端部力和齒槽力。力的類型物理來源位置周期主要影響因素抑制方向端部力鐵心兩端磁路截斷端部磁不對稱一個極距鐵心長度、端部齒形輔助齒、端部修正齒槽力永磁體與鐵心齒之間切向磁拉力與齒距相關(guān)具體取決于極槽配合槽口寬度、齒形、磁極形狀半閉口槽、磁極偏移磁阻力端部力與齒槽力之和繞組不通電時的合力兩者周期的最小公倍數(shù)上述因素疊加綜合優(yōu)化理解這個層次關(guān)系仿真分離才有意義。如果說“仿真磁阻力”其實得到的是三者合力要區(qū)分端部力和齒槽力各自的貢獻(xiàn)需要特殊的模型或后處理手段。5.2 磁阻力仿真電流為零時的推力成分磁阻力仿真直接復(fù)用空載模型條件繞組激勵設(shè)為零電流動子以恒定速度滑過一段距離。這樣得到的推力時變曲線就是磁阻力隨位置變化的曲線。為了獲得完整周期動子滑行距離至少要覆蓋一個位置周期也就是端部力周期和齒槽力周期的最小公倍數(shù)。實際仿真時多跑兩個周期更穩(wěn)妥既能確認(rèn)周期性又方便截取穩(wěn)定段。仿真結(jié)束后把橫軸從時間換算成位移就得到磁阻力-位移曲線。曲線峰峰值是磁阻力的核心指標(biāo)它直接影響定位精度和低速推力波動水平。如果峰值過大就要考慮抑制措施。注意一點磁阻力仿真電流雖然為零永磁體始終全程參與所以永磁體材料屬性和充磁方向必須保證正確。磁阻力計算結(jié)果對氣隙特別敏感建模時氣隙尺寸不能有偏差否則磁阻力會被系統(tǒng)性放大或縮小。5.3 端部力與齒槽力的分離方法要單獨提取端部力和齒槽力通用辦法是構(gòu)造一個無端部參考模型。思路是這樣的齒槽力本質(zhì)上是周期現(xiàn)象不依賴鐵心兩端截斷。如果構(gòu)造一臺鐵心長度兩端周期延伸、采用主從周期邊界條件的模型端部效應(yīng)被消除這臺模型里算出的定位力就是純齒槽力。然后拿完整模型的磁阻力減去純齒槽力差值就是端部力。這個操作在項目里實現(xiàn)起來不復(fù)雜在完整模型基礎(chǔ)上復(fù)制出一個周期模型修改邊界條件和幾何長度跑一遍同樣的零電流運動仿真然后做減法。周期模型的主從邊界需要保證端部匹配幾何長度必須是周期長度的整數(shù)倍否則邊界對不上算出來的齒槽力同樣失真。直線電機(jī)的周期長度最典型的情況是每個極距重復(fù)一次或者每兩個極距重復(fù)一次與繞組排布有關(guān)。也有更精細(xì)的做法通過改變鐵心長度做一組掃描利用長度變化時端部力相位偏移、齒槽力相位不變的特性用數(shù)學(xué)方法把兩者分離。這個方法理論上更嚴(yán)謹(jǐn)?shù)僮鲝?fù)雜偏科研向工程上用參考模型做減法已經(jīng)足夠準(zhǔn)確。5.4 抑制手段的仿真驗證思路分離出端部力和齒槽力之后優(yōu)化方向就清晰了齒槽力靠調(diào)整齒形、槽口寬度、永磁體磁極形狀來抑制端部力靠輔助齒、端部齒形修正、磁極錯位來抑制。在Maxwell里驗證這些抑制方案非常方便。每改一次幾何或材料參數(shù)重跑一遍零電流仿真對比磁阻力峰峰值和波形形狀的變化就能快速判斷措施是否有效。比如半閉口槽設(shè)計在二維模型里就是改一個槽口寬度參數(shù)磁極錯位就是在永磁體陣列里把某幾個磁極端部位置偏移一小段距離。這些結(jié)構(gòu)改動在參數(shù)化設(shè)計里都可以用變量控制跑一組參數(shù)掃描就能得到最優(yōu)值。需要提醒的是抑制端部力的措施和抑制齒槽力的措施有時候會互相影響。輔助齒改動端部鐵心形狀可能也會影響附近的齒槽力分布。所以最終方案要在完整模型里統(tǒng)一驗證不能只信任單獨齒槽模型的結(jié)果。這也是我一直建議把完整模型和參考模型都保留優(yōu)化階段交替使用的原因。6. 常見問題與排查技巧實錄6.1 反電動勢波形不對稱或幅值異常這個問題幾乎算是入門必踩坑。最常見的原因有三個。第一繞組激勵方向設(shè)置錯誤。Maxwell中每個繞組組成部分需要指定電流流入流出方向某個槽的導(dǎo)體方向設(shè)置反了該相反電動勢就會明顯異常。檢查方法是把所有槽導(dǎo)體的方向高亮出來逐個核對。第二永磁體充磁方向錯誤。平板型直線電機(jī)常用垂直充磁加輔助水平充磁的Halbach陣列某個磁極的充磁分量方向?qū)懛创艌龇植疾粚ΨQ反電動勢自然也不對稱。第三幾何模型不對稱。繞組端部或鐵心在建模時因布爾運算失誤出現(xiàn)微小偏差肉眼很難發(fā)現(xiàn)但電磁結(jié)果會暴露。碰到波形不對稱我排查順序是先看繞組分相和方向再看永磁體充磁最后查幾何前兩個基本覆蓋了絕大多數(shù)情況。6.2 推力波動異常、波形有毛刺或鋸齒感推力波形有毛刺通常不是電機(jī)問題而是網(wǎng)格和時間步的問題。氣隙網(wǎng)格太粗齒槽附近的磁導(dǎo)率變化無法準(zhǔn)確捕捉力曲線上就會疊加高頻毛刺。時間步長太大運動過程中的磁場變化被跨越力波形出現(xiàn)臺階狀。處理辦法是把氣隙網(wǎng)格沿厚度方向分到三層以上同時讓每個電周期采樣點不少于一百個。毛刺只在某一小段位移區(qū)間出現(xiàn)時檢查運動域與靜止氣隙邊界上網(wǎng)格尺寸差異統(tǒng)一交界區(qū)域網(wǎng)格通常能解決。還有一種鋸齒感來自磁阻力和電磁推力頻率不同。磁阻力以極距為周期電磁推力以電周期為周期兩者疊加后如果頻率接近波形上會出現(xiàn)明顯拍頻特征。這不是錯誤而是物理現(xiàn)象遇到這種情況先做零電流仿真分離磁阻力就會看得清楚。6.3 仿真速度慢、夭折或不收斂直線電機(jī)瞬態(tài)仿真的計算量集中在運動域網(wǎng)格重構(gòu)和每個時間步的非線性迭代。仿真慢或中途夭折優(yōu)先檢查三件事。一是運動域過大。運動域包住運動體即可包太多周圍空氣會讓運動域內(nèi)部網(wǎng)格數(shù)量猛增。二是氣隙網(wǎng)格過度細(xì)化。氣隙網(wǎng)格要密但不需要無節(jié)制加密三層左右通常就是收益上限。三是非線性迭代步數(shù)不夠。鐵心飽和程度高時默認(rèn)迭代次數(shù)可能不滿足收斂要求適當(dāng)調(diào)大最大迭代步數(shù)并觀察殘差曲線。如果仍然太慢可以用平移周期對稱性減少模型范圍前提是確認(rèn)不破壞端部效應(yīng)這在整機(jī)模型和參考模型配合下可以實現(xiàn)。6.4 力分量方向與實際不一致提取推力時發(fā)現(xiàn)力的符號相反或者力一直很小首先檢查運動方向與電流產(chǎn)生的行波磁場方向是否一致。直線電機(jī)推力方向由電流相序決定三相電流相序反了推力方向就反了。解決辦法是對調(diào)任意兩相電流的相位或者反過來設(shè)置動子運動方向。力結(jié)果分量的正負(fù)取決于模型坐標(biāo)。電機(jī)沿橫軸正方向運動時選中運動體提該方向分力即可符號正負(fù)代表力與坐標(biāo)正方向相同或相反。判斷推力方向最簡單的辦法是看電機(jī)實際應(yīng)該往哪個方向推動負(fù)載再跟仿真里的運動方向?qū)φ铡，F(xiàn)象可能原因排查方向反電動勢不對稱繞組方向錯、充磁錯、幾何不對稱方向高亮核對再查充磁反電動勢毛刺多時間步長過大、氣隙網(wǎng)格太粗采樣點加到100以上加密氣隙推力波動帶毛刺網(wǎng)格粗、時間步大、滑移邊界網(wǎng)格不匹配統(tǒng)一交界面網(wǎng)格尺寸加密時間步仿真不收斂運動域過大、網(wǎng)格過密、非線性迭代不足縮小運動域控制網(wǎng)格數(shù)量調(diào)大迭代次數(shù)推力方向反電流相序反交換任意兩相電流相位要說體會的話我最大的感受是直線電機(jī)的仿真建模只占三成功夫結(jié)果解讀占七成。模型建得再漂亮讀不對波形和力曲線等于白算。每次跑完仿真先問自己三個問題波形對不對數(shù)值是否符合物理直覺波動成分能不能用端部力、齒槽力、電磁諧波拆解出來。三條都對上了這個仿真才算真正做完。對于新手我更建議把空載反電動勢、零電流磁阻力和負(fù)載推力這三套仿真當(dāng)成固定套餐每個模型都完整跑一遍、對比一遍直線電機(jī)電磁設(shè)計的功底就是這樣慢慢攢出來的。