梅俊偉 劉保國(guó) 馮 偉

(①河南工業(yè)大學(xué)機(jī)電工程學(xué)院，河南 鄭州450001;②河南省超硬磨料磨削裝備重點(diǎn)實(shí)驗(yàn)室，河南 鄭州450001)

與普通磨削相比，超高速磨削(砂輪線速度>150 m/s)具有效率高、加工精度高、表面質(zhì)量好等優(yōu)點(diǎn)[1]。由于其內(nèi)置電機(jī)電磁損耗和軸承摩擦產(chǎn)生熱量，導(dǎo)致主軸產(chǎn)生較大的熱變形。研究表明，電機(jī)是電主軸的主要熱源之一[2-3]。隨著轉(zhuǎn)速的增加，主軸電機(jī)產(chǎn)生的溫升進(jìn)一步加劇，嚴(yán)重影響了高速電主軸的工作壽命和精度。因此，必須對(duì)電主軸內(nèi)置電機(jī)的電磁發(fā)熱進(jìn)行深入分析，以提高超高速磨削電主軸的工作性能。

目前，國(guó)內(nèi)外學(xué)者對(duì)電主軸內(nèi)置電機(jī)溫度場(chǎng)的分析進(jìn)行了大量研究。文獻(xiàn)[4]提出了一種電主軸三維有限元模型來(lái)預(yù)測(cè)高速電主軸熱性能。將內(nèi)置電機(jī)發(fā)熱量加載到溫度場(chǎng)模型中，但沒(méi)有考慮溫度變化對(duì)電機(jī)材料屬性的影響。文獻(xiàn)[5]通過(guò)使用有限元分析軟件ANSYS Workbench分析了電主軸系統(tǒng)在不同轉(zhuǎn)速下的溫度場(chǎng)分布和熱變形。文獻(xiàn)[6]分析了電主軸電機(jī)定、轉(zhuǎn)子間的氣隙厚度受到熱載荷和離心力作用的減小程度，計(jì)算了定、轉(zhuǎn)子間的對(duì)流換熱系數(shù)隨氣隙長(zhǎng)度的變化量，得到定、轉(zhuǎn)子間的熱分布規(guī)律。文獻(xiàn)[7]通過(guò)實(shí)驗(yàn)較為精確的獲得電主軸電動(dòng)機(jī)損耗，并將損耗以生熱量的方式加載到有限元模型中，但沒(méi)有將電主軸電動(dòng)機(jī)材料的溫度特性引起電磁損耗變化這一因素進(jìn)行分析。文獻(xiàn)[8]提出了一種新型軸芯冷卻結(jié)構(gòu)，建立了電主軸熱-結(jié)構(gòu)耦合分析模型，分析得到有、無(wú)軸心冷卻轉(zhuǎn)子的瞬態(tài)溫升規(guī)律。文獻(xiàn)[9]將電主軸內(nèi)置電機(jī)的損耗以生熱率的方式導(dǎo)入到溫度場(chǎng)模型中，運(yùn)用有限元法研究了電主軸熱穩(wěn)態(tài)下的溫度場(chǎng)分布，但沒(méi)有考慮內(nèi)置電機(jī)材料的溫度特性。文獻(xiàn)[10]采用有限元法以同步電機(jī)電主軸為研究對(duì)象，分析了電主軸電機(jī)定子和轉(zhuǎn)子的溫度場(chǎng)分布以及熱變形。文獻(xiàn)[11]將計(jì)算得到的電主軸內(nèi)置永磁電機(jī)的電磁損耗加載到有限元熱分析模型中，分析了永磁同步電主軸在啟動(dòng)、穩(wěn)定運(yùn)行以及不同轉(zhuǎn)速下各關(guān)鍵部位的溫度變化，但忽略了內(nèi)置電機(jī)溫度變化對(duì)材料屬性的反饋?zhàn)饔?。文獻(xiàn)[12]采用熱阻節(jié)點(diǎn)網(wǎng)絡(luò)法分析了高速磨削電主軸在不同轉(zhuǎn)速、冷卻水流量及冷卻水入口溫度條件下電主軸電機(jī)定子的溫升規(guī)律。

通過(guò)以上分析可以看出，現(xiàn)有文獻(xiàn)對(duì)電主軸內(nèi)置電機(jī)溫度場(chǎng)的研究多集中于電磁損耗影響溫度場(chǎng)的單向耦合，沒(méi)有考慮溫度升高對(duì)電主軸內(nèi)置電機(jī)材料屬性的影響。隨著內(nèi)置電機(jī)溫度的升高，繞組電阻率和永磁體磁性會(huì)產(chǎn)生相應(yīng)的變化，進(jìn)而會(huì)引起繞組銅耗和永磁體渦流損耗的變化，采用單向耦合分析可能造成內(nèi)置電機(jī)的溫度場(chǎng)產(chǎn)生較大的誤差。因此，本文采用電磁場(chǎng)與溫度場(chǎng)雙向耦合的方法，考慮了溫度變化對(duì)電主軸內(nèi)置電機(jī)繞組銅線電阻率和永磁體磁性的影響。首先，通過(guò)電磁場(chǎng)分析得到內(nèi)置電機(jī)的電磁損耗，其次，將電磁損耗加載到溫度場(chǎng)模型中得到內(nèi)置電機(jī)繞組的溫度，然后，根據(jù)溫度值改變繞組和永磁體的材料參數(shù)，再次進(jìn)行電磁場(chǎng)分析以更新電磁損耗，依次循環(huán)，直到溫度值達(dá)到穩(wěn)定。

1 有限元分析模型建立

本文以超高速磨削電主軸為研究對(duì)象進(jìn)行仿真分析，電主軸內(nèi)置電機(jī)的主要參數(shù)如表1所示。

表1 電主軸內(nèi)置電機(jī)參數(shù)

為便于建立溫度場(chǎng)求解模型，將電主軸內(nèi)置電機(jī)雙層繞組分別等效為等長(zhǎng)的直導(dǎo)體,由于該電主軸為對(duì)稱結(jié)構(gòu)，為了減少計(jì)算時(shí)間和節(jié)約計(jì)算機(jī)資源，取磨削電主軸的1/6模型進(jìn)行分析。模型省去了連接螺釘、倒角、軸承預(yù)緊裝置等，電主軸簡(jiǎn)化分析模型如圖1所示。

2 溫度對(duì)材料的影響

電主軸磁熱耦合分析主要基于內(nèi)置電機(jī)的材料屬性隨溫度變化而建立的。隨著溫度的升高，內(nèi)置電機(jī)的繞組銅線電阻率以及永磁體的磁性會(huì)發(fā)生變化，這些會(huì)改變電主軸內(nèi)置電機(jī)的電磁損耗，進(jìn)而影響內(nèi)置電機(jī)的溫度場(chǎng)分布[13]。

2.1 繞組銅線電阻率的溫度特性

電主軸內(nèi)置電機(jī)繞組材料為銅。電主軸在工作時(shí)，繞組電阻率直接影響內(nèi)置電機(jī)銅耗的產(chǎn)生，進(jìn)而影響內(nèi)置電機(jī)的溫升。繞組電阻率隨溫度的變化采用下式計(jì)算[14]:

(1)

式中：ρt=ρ15分別為溫度在t、15 ℃時(shí)的繞組的電阻率值，取ρ15=1.7×10-8Ω·m。

2.2 永磁體的溫度特性

永磁體的磁性隨著溫度的升高而發(fā)生變化，從而改變內(nèi)置電機(jī)的整體磁場(chǎng)，與磁場(chǎng)相關(guān)的電磁損耗也會(huì)發(fā)生變化。溫升對(duì)永磁體磁性的影響如下式[15]：

(2)

式中:Brt0：t0℃時(shí)的剩磁密度；Hct0：t0℃時(shí)的計(jì)算矯頑力；αBr：Br的可逆溫度系數(shù)；αHr：Hc的可逆溫度系數(shù)；IL：Br的不可逆損失率；t1：永磁體工作溫度，℃。

內(nèi)置電機(jī)永磁體的電導(dǎo)率對(duì)渦流損耗影響較大，通過(guò)計(jì)算永磁體上的電阻率來(lái)間接計(jì)算電導(dǎo)率。溫度對(duì)永磁體電阻率ρm的影響為[16]

ρm=c·T2+b·T+a

(3)

式中：a、b、c為與材料有關(guān)的系數(shù)，本文中電主軸內(nèi)置電機(jī)的永磁體采用釹鐵硼永磁材料，其中c=-5.468×10-6,b=-1.765×10-3,a=1.520。

3 磁熱雙向耦合計(jì)算流程

本文通過(guò)Maxwell有限元分析軟件對(duì)電主軸內(nèi)置電機(jī)進(jìn)行電磁場(chǎng)分析，通過(guò)ANSYS Workbench中的Transient Thermal進(jìn)行瞬態(tài)溫度場(chǎng)分析。將電磁場(chǎng)分析得到的電磁損耗導(dǎo)入到溫度場(chǎng)中計(jì)算溫度場(chǎng)分布，然后根據(jù)溫度場(chǎng)改變電磁場(chǎng)材料屬性參數(shù)以更新電磁損耗，使電磁場(chǎng)與溫度場(chǎng)的數(shù)據(jù)實(shí)現(xiàn)雙向傳遞，從而實(shí)現(xiàn)磁熱雙向耦合的目的。電主軸磁熱雙向耦合分析流程圖如圖2所示。具體步驟如下：

(1)通過(guò)內(nèi)置電機(jī)的初始參數(shù)分別建立電磁場(chǎng)和溫度場(chǎng)有限元分析模型。

(2)在電磁場(chǎng)分析模型中分析得到電磁損耗，并將損耗作為熱源導(dǎo)入到溫度場(chǎng)分析模型。

(3)在溫度場(chǎng)分析模型中分析得到內(nèi)置電機(jī)的溫度場(chǎng)分布，并將溫度值反饋到電磁場(chǎng)模型中。

(4)在電磁場(chǎng)分析模型中根據(jù)相應(yīng)的溫度變化值改變材料的屬性(繞組電阻率和永磁體磁性)，通過(guò)電磁場(chǎng)分析，將更新的電磁損耗導(dǎo)入到溫度場(chǎng)分析模型中重新計(jì)算溫度場(chǎng)分布。

(5)重復(fù)以上步驟，直到溫度值收斂。

4 仿真分析

通過(guò)采用磁熱雙向耦合的方法分析磨削電主軸在工作轉(zhuǎn)速為10 000 r/min、軸向載荷為400 N、徑向載荷為1 000 N時(shí)的溫度場(chǎng)分布。本文所用電主軸的軸承采用油氣潤(rùn)滑，殼體采用循環(huán)水冷的方式。經(jīng)過(guò)計(jì)算，電主軸的熱載荷和邊界條件如表2所示。

表2 電主軸熱載荷和邊界條件

圖3和圖4分別為通過(guò)單向耦合和雙向耦合計(jì)算得到的繞組的溫度場(chǎng)仿真云圖。由溫度場(chǎng)云圖可以看出，繞組的最高溫度出現(xiàn)在端部，單向耦合計(jì)算的最高溫度為42.385 ℃，雙向耦合計(jì)算的最高溫度為44.853 ℃。

圖5所示為繞組耦合分析的溫度變化曲線。由于繞組的電阻率隨著溫度的升高而逐漸增大，繞組銅耗與電阻率成正相關(guān)。因此，經(jīng)過(guò)雙向耦合計(jì)算，繞組溫度呈現(xiàn)上升趨勢(shì)。由圖中數(shù)據(jù)可以看出，經(jīng)過(guò)5次耦合計(jì)算后，繞組的溫度值趨于平穩(wěn)，溫度逐漸實(shí)現(xiàn)收斂。

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

為了驗(yàn)證上述電主軸磁熱耦合方法的正確性，對(duì)磨削電主軸進(jìn)行溫升試驗(yàn)，溫升測(cè)試試驗(yàn)平臺(tái)如圖6所示。

試驗(yàn)前電主軸處于不工作狀態(tài)24 h以上，以使電主軸各部位溫度與外界環(huán)境溫度盡量保持一致。使用預(yù)埋在繞組端部的溫度傳感器PT100進(jìn)行溫度測(cè)量。如圖7所示分別為單向耦合、雙向耦合和實(shí)驗(yàn)測(cè)量得到的繞組端部的溫度隨時(shí)間變化的曲線圖。由圖中可以得出，考慮溫度對(duì)材料性能影響的雙向耦合得到的結(jié)果更加接近于實(shí)際值。

由圖7可以看出，試驗(yàn)值在穩(wěn)定階段呈現(xiàn)上下波動(dòng)的趨勢(shì)，主要是由于當(dāng)循環(huán)冷卻水溫度達(dá)到水冷機(jī)設(shè)定溫度(28 ℃)時(shí)，水冷機(jī)開(kāi)始工作，使循環(huán)冷卻水溫度降低到25 ℃，之后依此循環(huán)。因此，繞組端部的溫度在穩(wěn)定階段呈現(xiàn)周期性的變化。由于本次試驗(yàn)的負(fù)載較小，繞組電流較小，繞組溫升較小。因此，通過(guò)單向耦合與雙向耦合分析得到的繞組端部的溫度相差較小。

6 結(jié)語(yǔ)

本文考慮溫度變化對(duì)繞組和永磁體材料屬性的影響，采用磁熱雙向耦合的方法，對(duì)磨削電主軸進(jìn)行瞬態(tài)溫度場(chǎng)分析?？梢缘贸觯翰牧系臏囟忍匦詴?huì)對(duì)磨削電主軸內(nèi)置電機(jī)的溫度場(chǎng)產(chǎn)生影響。通過(guò)仿真與試驗(yàn)結(jié)果分析對(duì)比也說(shuō)明，通過(guò)磁熱雙向耦合的分析方法，可以提高磨削電主軸溫度場(chǎng)分析的準(zhǔn)確性。