歡迎參加第十一屆產(chǎn)品創(chuàng)新數(shù)字化國際峰會 12月3-4日 杭州
0 引言
現(xiàn)代高性能永磁電機主要由永磁電機本體和驅(qū)動控制器構(gòu)成,永磁同步電動機具有結(jié)構(gòu)簡單、體積小、重量輕和高效節(jié)能等一系列優(yōu)點,近年來得到了迅速發(fā)展。高性能稀土永磁材料的出現(xiàn),其優(yōu)越的磁性能和相對較低的價格,使得高性能永磁同步電機的開發(fā)和研究成為世界各國的熱點,并在國防、工農(nóng)業(yè)生產(chǎn)和日常生活等方面獲得越來越廣泛的應(yīng)用。
永磁電機,特別是內(nèi)嵌磁鋼轉(zhuǎn)子永磁電機(IPM),其結(jié)構(gòu)復雜,傳統(tǒng)的磁路法已經(jīng)無法準確計算磁路和電機性能,需要借助高性能的有限元磁場求解工具提升仿真精度和效率。同時,高性能永磁電機應(yīng)用在調(diào)速、伺服控制等場合,在電機產(chǎn)品開發(fā)初期需要結(jié)合控制電路和逆變器系統(tǒng)考慮設(shè)計和優(yōu)化電機性能。
1 永磁體渦流損耗高效精確計算
永磁電機轉(zhuǎn)子與定子基波磁勢是同步旋轉(zhuǎn)的,因此通常計算中忽略轉(zhuǎn)子內(nèi)的永磁體渦流損耗。而實際永磁電機中,由于存在齒槽效應(yīng),且繞組磁動勢的非正弦分布,或者由PWM逆變器引入高次諧波電流等,均會產(chǎn)生諧波磁勢,從而導致基波轉(zhuǎn)速下的轉(zhuǎn)子永磁體及固定永磁體的金屬護套中引起渦流損耗。
通常情況下,與電機定子的繞組銅線損耗和鐵損相比,轉(zhuǎn)子磁鋼渦流損耗占比很小。但是,由于轉(zhuǎn)子散熱條件相對封閉,熱量不容易散發(fā),磁鋼渦流損耗可能會引起磁鋼內(nèi)局部高溫升,從而引起永磁體局部熱退磁。特別是燒結(jié)釹鐵硼(NdFeB)具有較大電導率和較低的居里溫度,更需要特別計算渦流損耗和校核溫升。因此,設(shè)計永磁電機初期,就需要精確計算磁鋼渦流損耗,保證電機磁鋼穩(wěn)定的熱性能。
1.1 渦流效應(yīng)趨膚深度計算
永磁體的損耗,主要由氣隙高次諧波含量在磁鋼中感應(yīng)渦流導致,因此,準確計算高次諧波在永磁體中的趨膚深度是精確計算磁鋼渦流損耗的前提保證。本文采用Maxwell有限元計算直接計算磁鋼損耗,而電磁場問題實現(xiàn)數(shù)值求解的一個關(guān)鍵步驟就是網(wǎng)格剖分。因此網(wǎng)格剖分的規(guī)模和質(zhì)量直接影響數(shù)值計算的精度、效率和經(jīng)濟性。
對于較高頻率的高次諧波電磁場有限元分析,計算磁鋼渦流損耗首要考慮到其趨膚深度的影響。當磁場進入磁鋼表面以內(nèi)時,其最大值按指數(shù)規(guī)律衰減,因此,隨著深度的增加,磁場也將按指數(shù)規(guī)律衰減。
磁場在導體內(nèi)的趨膚深度,可由公式1進行計算:
式中,δ為趨膚深度,ω為角頻率=2*π*f,σ為磁鋼電導率,μ為磁鋼相對磁導率。
計算不同頻率下的釹鐵硼磁鋼的趨膚深度,結(jié)果如表1所示:
表1 不同頻率下釹鐵硼永磁體的趨膚深度
ANSYS Maxwell提供了軟件的易用性和方便性,能夠自動依據(jù)磁鋼材料、銅線等導電材料的電導率,在指定頻率下自動計算趨膚深度,如圖1所示銅線的Skin Depth的數(shù)據(jù)由軟件自動功能自動計算出在1kHz下趨膚深度為2.0898mm,并分為4層。
圖1 ANSYS Maxwell自動計算趨膚深度
經(jīng)典有限元算法為了能夠達到較高的計算精度,在趨膚深度內(nèi)劃分的有限元網(wǎng)格一般來說必須足夠小。但是太小的網(wǎng)格會增加計算機運算開銷,因此ANSYS Maxwell R16中改進了趨膚深度的TAU網(wǎng)格剖分技術(shù),優(yōu)化后的趨膚深度網(wǎng)格能夠在法向方向上體現(xiàn)網(wǎng)格極小尺寸,而垂向量方向網(wǎng)格尺寸依然可以足夠大,既可以明顯降低網(wǎng)格數(shù)量,又節(jié)約運算資源,數(shù)值計算結(jié)果精度與傳統(tǒng)網(wǎng)格技術(shù)相當,非常適合磁鋼渦流損耗仿真計算。
如圖2所示,在二維有限元中,與傳統(tǒng)有限元算法比較,網(wǎng)格數(shù)量減少到二十分之一,同樣能獲得幾乎等值的計算精度,計算效率卻成倍顯著提高,所以ANSYS Maxwell R16.0新的TAU網(wǎng)格剖分技術(shù)優(yōu)勢非常明顯。
圖2 R16.0新TAU網(wǎng)格剖分技術(shù)優(yōu)勢對比分析
1.2 周期模型網(wǎng)格克隆技術(shù)
永磁電機往往由多個磁極組成,且具有周向?qū)ΨQ性,且電機定子鐵芯、轉(zhuǎn)子鐵芯、磁鋼等,軸向具有拉伸對稱性,為了更進一步減小網(wǎng)格數(shù)量,ANSYS Maxwell R16在2D和3D中分別引進了全新的網(wǎng)格克隆技術(shù)(Clone Mesh),從周向和軸向分別減少網(wǎng)格剖分數(shù)量,提高網(wǎng)格復用性和對稱性。
如圖3所示,打開Clone Mesh選項,定義4極電機父磁鋼幾何位置的扇形區(qū)域,設(shè)置此區(qū)域內(nèi)的Skin Depth選項,生成網(wǎng)格后的效果如圖3右側(cè)所示,4個磁鋼中除扇形區(qū)域之外的3個磁鋼網(wǎng)格直接網(wǎng)格克隆而來,所以4個磁鋼的網(wǎng)格完全相同。值得注意的是,盡管磁鋼的網(wǎng)格克隆而來,但是僅針對初始網(wǎng)格,如果使用Maxwell自動自適應(yīng)網(wǎng)格剖分技術(shù)迭代網(wǎng)格,則迭代后的4個磁鋼網(wǎng)格不再是網(wǎng)格克隆關(guān)系,且Maxwell有限元矩陣運行時所有網(wǎng)格節(jié)點均參與計算。
圖3 Maxwell 2D Clone mesh設(shè)置和初始網(wǎng)格效果
當采用Maxwell 3D Clone Mesh技術(shù)后,電機模型軸向網(wǎng)格更趨于統(tǒng)一,能大規(guī)模降低3D模型的網(wǎng)格量,如圖4所示,Maxwell 3D Clone Mesh使用前后網(wǎng)格對比圖,改善后的網(wǎng)格數(shù)量只有傳統(tǒng)網(wǎng)格數(shù)量的63%。
圖4 Maxwell 3D Clone mesh結(jié)果對比圖
1.3 計算效率提升
在網(wǎng)格技術(shù)改進的同時,Maxwell求解算法也致力與提升求解效率,不斷加強的高性能并行求解技術(shù)(ANSYS HPC)為大型三維模型的求解提供了高效運算技術(shù),相對于單核技術(shù),多核并行求解,可以獲得線性加速效果。即使不使用多核并行求解技術(shù),Maxwell新算法也能顯著提升運算效率,表2統(tǒng)計了Maxwell 十年來計算效率提升水平,幾乎提升了60倍。
表2 Maxwell 3D運算效率提升對比表
2 有限元降階模型提取
永磁電機調(diào)速時需要結(jié)合控制電路考慮逆變電流對電機的控制和性能影響,為了將Maxwell 2D/3D電磁模型和控制電路模型耦合分析,ANSYS Simplorer提供了無與倫比的實時耦合-場路協(xié)同仿真技術(shù)。然而,Maxwell 3D直接與Simplorer實時耦合仿真時,由于時間步長由Simplorer控制,實時計算時電腦CPU和內(nèi)存開銷大,因此為了快速求解場路耦合問題,ANSYS提供了獨有的Maxwell 2D/3D有限元模型降階輸出功能。
ANSYS MaxwellR16軟件中開發(fā)的獨特降階ECE模型輸出功能,直接在瞬態(tài)磁場外電路中引入ECE控制模塊,通過定義Maxwell circuit editor中三相電機的ECE3_model和ECER_model,快捷抽取電機的E-Model模型。
ECE抽取工作直接由Maxwell軟件自動完成,僅需要分別對ECE3_model和ECER_model定義參數(shù)。如圖5所示,其中Windings輸入電流掃描的繞組,以逗號相隔;CurrentSweeps定義電流掃描值和范圍,圖示意義為以0電流為中間值,正負各掃描5個點,電流值間隔2A。圖6顯示了ECER_model的定義和使用,其中RotAngMax是最大掃描角度,對于三相電機而言只需要60度電角度,這里需要填入機械角度;RotAngIntervals定義角度掃描的點數(shù)。在本例中,角度周期是30度、掃描30個點;2相電流各掃描11個點,所以總的計算點數(shù)是11*11*30=3630點,即在Maxwell瞬態(tài)場中會計算3630個時間點,啟動仿真后自動計算并輸出降階模型,不需要手動操作,十分方便。
圖5 Maxwell 3D不同版本軟件性能改進對比柱狀圖
圖6 Maxwell外電路中的ECE3_model模型參數(shù)
圖7 Maxwell外電路中的ECER_model模型參數(shù)
3 場路耦合協(xié)同仿真
求解完成后,在Maxwell對應(yīng)的result結(jié)果文件夾下會自動生成一個ece_model.sml文件,直接輸入到ANSYS Simplorer控制電路系統(tǒng)工具中。本論文案例永磁電機為圖3所示電機模型,在Maxwell中定義5*sin(2*pi*50*time)電流源激勵輸入。
ANSYS Simplorer中構(gòu)建與Maxwell輸入相同電流源控制電路,如圖8所示,仿真相同電流激勵下Maxwell直接輸出結(jié)果,以及Simplorer中ECE模型輸出結(jié)果,對比驗證ECE模型計算精度。
圖8 Simplorer三相電流源控制電路
圖9和圖10分別輸出轉(zhuǎn)矩和C相反電動勢波形對比,實線為Simplorer耦合仿真波形,虛線為Maxwell仿真結(jié)果的波形導入到Simplorer中,從波形對比可知,二者吻合度極高,從而驗證了ECE降階模型的精確度。而此狀態(tài)下Simplorer的仿真速度非常快,數(shù)秒時間便已經(jīng)完成,僅為Maxwell仿真時間的幾十分之一,大大提高了永磁電機仿真分析中場路耦合協(xié)同仿真的時效性和可行性,以便研究高保真永磁電機模型下的各種復雜控制電路和控制策略。
圖9 Simplorer+ECE降階模型vs Maxwell時域仿真負載轉(zhuǎn)矩波形
圖10 Simplorer+ECE降階模型vs Maxwell時域仿真繞組感應(yīng)電勢波形
4 結(jié)論
本文闡述了基于高速發(fā)展的CAE工程仿真領(lǐng)域領(lǐng)導企業(yè)ANSYS有限元電磁場數(shù)值計算工具,運用Maxwell和Simplorer 2015年R16新版本,對高性能永磁電動機磁鋼渦流損耗計算、與控制電路場合協(xié)同等研究熱點應(yīng)用進行分析,提供高效仿真分析方法和手段,以便提升永磁電機設(shè)計性能和提高產(chǎn)品仿真能力,最終加速產(chǎn)品上市時間。
核心關(guān)注:拓步ERP系統(tǒng)平臺是覆蓋了眾多的業(yè)務(wù)領(lǐng)域、行業(yè)應(yīng)用,蘊涵了豐富的ERP管理思想,集成了ERP軟件業(yè)務(wù)管理理念,功能涉及供應(yīng)鏈、成本、制造、CRM、HR等眾多業(yè)務(wù)領(lǐng)域的管理,全面涵蓋了企業(yè)關(guān)注ERP管理系統(tǒng)的核心領(lǐng)域,是眾多中小企業(yè)信息化建設(shè)首選的ERP管理軟件信賴品牌。
轉(zhuǎn)載請注明出處:拓步ERP資訊網(wǎng)http://www.ezxoed.cn/
本文標題:2015產(chǎn)品創(chuàng)新數(shù)字化征文:基于Maxwell的高性能永磁電機和控制系統(tǒng)仿真分析
本文網(wǎng)址:http://www.ezxoed.cn/html/solutions/14019318872.html