陳清華，林漢毅，許斌

（1.安徽理工大學(xué) 省部共建深部煤礦采動響應(yīng)與災(zāi)害防控國家重點(diǎn)實(shí)驗(yàn)室，安徽淮南 232001；2.安徽理工大學(xué) 礦山智能裝備與技術(shù)安徽省重點(diǎn)實(shí)驗(yàn)室，安徽淮南 232001）

超磁致伸縮材料（Giant magnetostrictive material,GMM）在外部磁場的作用下會按一定規(guī)律發(fā)生形變，即磁致伸縮現(xiàn)象，利用這個變化特性可以制成超磁 致 伸 縮 致 動 器 （Giant magnetostrictive actuator,GMA）[1-3]。與傳統(tǒng)微定位工作平臺的驅(qū)動器相比，GMA具有響應(yīng)速度快、能量轉(zhuǎn)換效率高、輸出力大等優(yōu)點(diǎn)[4-6]。當(dāng)GMA采用勵磁線圈作為磁場源驅(qū)使GMM棒產(chǎn)生位移時會產(chǎn)生大量的熱，從而導(dǎo)致GMM棒磁變能力衰減及自身熱變形[7-9]。因此為了保證GMA的輸出精度，必須采取一定的熱變形控制措施以降低或消除溫度變化所帶來的影響。

目前，控制GMA溫升的主要方法有兩大類，一類是被動補(bǔ)償法，一類是主動溫控法。被動補(bǔ)償法能補(bǔ)償熱膨脹所引起的誤差，但補(bǔ)償控制算法十分復(fù)雜，而且不能補(bǔ)償GMM棒因磁致伸縮系數(shù)變化所引起的誤差[10-12]。主動溫控法需要增加溫度控制裝置，但是卻能夠控制GMM棒的溫升，不僅能夠抑制GMM棒自身的熱變形，還能保證磁致伸縮系數(shù)的穩(wěn)定性[13-15]。因此，主動溫控法逐漸受到國內(nèi)外學(xué)者的關(guān)注，并取得廣泛的應(yīng)用[16-18]。

本文設(shè)計(jì)了一種內(nèi)管外腔式的GMA冷卻系統(tǒng)，根據(jù)歐姆定律，得出驅(qū)動線圈的熱阻功耗。另外，基于傅里葉定律、牛頓冷卻定律以及等效熱阻定律[19-20]，建立了GMA的傳熱模型。在對不同工況下GMA的溫度場進(jìn)行有限元仿真后，通過比較同等工況自然對流和強(qiáng)制對流情況下的GMA溫度場，對冷卻系統(tǒng)的溫控特性進(jìn)行評估，分析驅(qū)動電流強(qiáng)度和冷卻液流速對GMA的溫度分布以及GMM熱變形的影響。

1 GMA結(jié)構(gòu)及工作原理

內(nèi)管外腔式GMA結(jié)構(gòu)如圖1所示。套筒、前導(dǎo)磁端蓋、前墊片推桿、GMM棒、后墊片、冷卻腔骨架、后導(dǎo)磁端蓋共同構(gòu)成了GMA的閉合磁路。前導(dǎo)磁端蓋和前墊片推桿之間設(shè)置了碟簧，通過調(diào)整后導(dǎo)磁端蓋上的預(yù)緊螺釘，可實(shí)現(xiàn)GMM棒預(yù)應(yīng)力的調(diào)節(jié)。進(jìn)水銅管、GMM棒上纏繞的銅管與線圈外側(cè)設(shè)置的冷卻腔骨架形成了GMA的冷卻結(jié)構(gòu)，可及時帶走線圈產(chǎn)生的熱量。另外，為降低線圈向GMM棒的傳熱量，在銅管和線圈之間設(shè)置具有絕熱性能的尼龍線圈骨架。

圖1 GMA 結(jié)構(gòu)示意圖

GMA的工作原理為：當(dāng)線圈中通入驅(qū)動電流時，在GMM棒周圍會產(chǎn)生驅(qū)動磁場，同時為了消除倍頻現(xiàn)象，需要在線圈中通入偏執(zhí)電流以產(chǎn)生勵磁磁場，GMM棒在兩種磁場作用下伸長，從而驅(qū)使推桿輸出位移。根據(jù)GMA對輸出位移大小的要求，通過調(diào)節(jié)驅(qū)動電流的大小，GMM棒可以獲得不同強(qiáng)度的驅(qū)動磁場而產(chǎn)生不同大小的伸長量。

2 GMA傳熱分析

2.1 熱源分析

GMA工作時的熱量來源主要有3部分：線圈通入驅(qū)動電流時產(chǎn)生的熱量Qs1，通入偏置電流時產(chǎn)生的熱量Qs2，GMM棒的渦流和磁滯損耗Qs3，由此可以得到GMA產(chǎn)生的總熱量

當(dāng)驅(qū)動線圈工作在低頻或直流狀態(tài)時，GMM棒上產(chǎn)生的渦流損耗及磁滯損耗與線圈上的歐姆損耗相比很小，可以忽略不計(jì)，所以GMA的主要熱源為線圈歐姆損耗產(chǎn)生的焦耳熱， GMA的總產(chǎn)熱量為

線圈的歐姆損耗為：

式中：Qi為線圈的歐姆損耗；I為控制電流；Rc線圈的總電阻；ρc為導(dǎo)線電阻率；Ac為導(dǎo)線截面積；Lc為導(dǎo)線總長度；nc為導(dǎo)線總匝數(shù)；do為線圈外徑；di為線圈內(nèi)徑。

2.2 傳熱理論分析

GMA內(nèi)部發(fā)熱，傳熱以及與恒溫冷卻水的換熱過程復(fù)雜，假設(shè)材料各項(xiàng)同性，以帶內(nèi)熱源的非穩(wěn)態(tài)導(dǎo)熱微分方程來描述GMA的溫度分布[21]，即

式中：t為溫度；τ為時間；a為熱擴(kuò)散率；r、z、θ為圓柱坐標(biāo)系的坐標(biāo)；為單位體積產(chǎn)熱量；ρ為物體的密度；C為物體比熱。

忽略碟簧、墊片、推桿、螺釘?shù)任⑿〗Y(jié)構(gòu)的影響，對GMA傳熱過程進(jìn)行分析，其傳熱模型如圖2所示。線圈作為GMA的內(nèi)熱源，其熱量主要沿兩個方向傳遞：其一，通過線圈骨架向內(nèi)傳遞，向內(nèi)傳遞的熱量一部分在傳向銅管時被管內(nèi)冷卻水帶走，另一部分則通過銅管及冷卻水傳如GMM棒；其二，通過冷卻腔骨架向外傳遞，向外傳遞的熱量一部分被腔內(nèi)冷卻水帶走，另一部分通過冷卻水和套筒傳到空氣中。

圖2 GMA 傳熱模型

通過對GMA傳熱過程分析，GMA的傳熱模型可簡化為一維導(dǎo)熱問題，由此建立GMA等效熱阻模型如圖3所示。線圈產(chǎn)生的熱量Qs近似為一個熱源環(huán)[22]，Qi、Qo分別為自熱源環(huán)向內(nèi)和向外傳遞的熱量，Qwi為管內(nèi)冷卻水帶走的熱量，Qwo為腔內(nèi)冷卻水帶走的熱量，QM為銅管傳向GMM棒的熱量，Qa為套筒傳向空氣的熱量。Rci、Rco分別為熱源環(huán)內(nèi)側(cè)和外側(cè)的線圈熱阻，Rf、Rq、Rt分別為線圈骨架、冷卻腔骨架、套筒的熱阻，Rp1、Rp2分別為與線圈骨架接觸側(cè)及與GMM棒接觸側(cè)銅管熱阻，Rwi、Rwo分別為冷卻水與銅管對流和冷卻水與腔體對流的熱阻，Ra為與空氣對流的熱阻。

圖3 GMA 等效熱阻模型

根據(jù)式（6）和熱阻的定義可得：

式中：λc、λf、λq、λt分別為線圈（銅管）、線圈骨架、冷卻腔骨架、套筒的導(dǎo)熱系數(shù)；lc、lp分別為線圈軸向長度和銅管長度；hi為銅管和水的對流換熱系數(shù)；ho為冷卻腔骨架和水的對流換熱系數(shù)；ha為套筒和空氣的對流換熱系數(shù)。

當(dāng)GMA溫度未達(dá)到穩(wěn)態(tài)時，線圈產(chǎn)生的熱量一部分被GMM棒所吸收，另一部分通過對流換熱帶走，即

當(dāng)GMA溫度達(dá)到穩(wěn)態(tài)時，所有熱量都將通過對流換熱帶走，此時QM= 0，TM = Twi，Qi = Qwi，其中TM為GMM棒的溫度，Twi為管內(nèi)冷卻水的平均溫度。

類比電學(xué)中的歐姆定律得：

式中：Ts、TA、Two分別為熱源面溫度、環(huán)境溫度、腔內(nèi)冷卻水平均溫度。

銅管和冷卻腔內(nèi)冷卻水帶走的熱量為：

式中：To、Tout_i分別為銅管進(jìn)、出口冷卻水溫度；Tin_o、Tout_o分別為冷卻腔進(jìn)、出口冷卻水溫度；Cp為水的比熱；qm冷卻水的質(zhì)量流量。

由于流出銅管的冷卻水直接流入冷卻腔，因此可得Tin_o= Tin_i。

結(jié)合式（10）~式（12），各支路的等效熱阻分別為：

因此GMM棒在穩(wěn)態(tài)時溫度為

由式（18）可知，GMM棒的溫度與GMA的產(chǎn)熱功率成正比，與冷卻水的流量成反比。因此，可以通過調(diào)節(jié)冷卻水的流量對GMM棒的溫度進(jìn)行控制。

3 仿真分析

3.1 條件設(shè)定

對GMA的結(jié)構(gòu)進(jìn)行簡化，忽略螺釘?shù)任⑿〗Y(jié)構(gòu)的影響，利用FloEFD熱分析軟件對GMA進(jìn)行仿真分析主要步驟如下：

1）建立模型：根據(jù)圖1所示GMA結(jié)構(gòu)示意圖，按照表1所示關(guān)鍵參數(shù)建立三維模型，并為零件添加材料屬性以及劃分網(wǎng)格。

表1 GMA結(jié)構(gòu)參數(shù)表

2）確定初始條件和邊界條件：線圈中電流大小與產(chǎn)熱功率的關(guān)系如表2所示，根據(jù)線圈中電流的大小，對其施加相應(yīng)的功率載荷；GMA各部件和環(huán)境的初始溫度設(shè)置為25 ℃，冷卻水入口溫度為25 ℃；GMA殼體外表面與空氣進(jìn)行自然對流換熱，參照推薦范圍[23]，設(shè)置換熱系數(shù)為 5 W/（m2·℃）。

表2 線圈電流大小與產(chǎn)熱功率表

3.2 結(jié)果分析

線圈中分別通入1 A和6 A的電流時，自然對流條件下即冷卻水靜止?fàn)顟B(tài)下的GMA溫度分布云圖如圖4a）、圖4b）所示，強(qiáng)制對流條件下即冷卻水以0.5 m/s狀態(tài)下的GMA溫度分布云圖如圖4c）、圖4d）所示。在對應(yīng)的不同條件下，GMA各部分的最高溫度如表3所示。

圖4 GMA 溫度分布云圖

表3 兩種入口流速下GMA溫度比較表

從圖4和表3可以看出，當(dāng)線圈中分別通入1 A和6 A的電流時，在自然對流條件下，GMA中各部件溫度劇烈上升，溫度最高分別達(dá)到了28.77 ℃和120.64 ℃，GMM棒溫度分別達(dá)到了 28.71 ℃和118.88 ℃；而在強(qiáng)制對流條件下，雖然GMA中各部件的溫度都有上升，但其溫升都得到了有效地控制，溫度最高分別為25.14 ℃和29.85 ℃，GMM棒溫度分別為 25.03 ℃ 和 26.12 ℃。

隨著工作電流的增大，線圈的產(chǎn)熱功率也逐漸增大，在自然對流和強(qiáng)制對流條件下，GMM棒的溫度隨電流變化如圖5所示。

圖5 電流大小對 GMA 棒溫度的影響

從圖5中可以看出，在冷卻水為0 時，GMM棒的最高溫度隨電流的增大而急劇上升，最高溫度可達(dá)到118.98 ℃。當(dāng)冷卻水在以0.5 m/s的入口速度下流動時，GMM棒的最高溫度雖然會隨著電流的增大而上升，但GMM棒的溫升卻得到了明顯的控制，最高溫度僅達(dá)到了26.19 ℃。

根據(jù)理論分析可知，GMM棒的溫度與冷卻水入口流速成正比，因此為了得到冷卻水流速對GMM棒溫度的影響，對不同冷卻水入口速度進(jìn)行了仿真，仿真結(jié)果如圖6所示。

圖6 冷卻水流速對GMM棒溫度的影響

由圖6可知，當(dāng)冷卻水流速一定時，線圈中的電流值越大，GMM的溫度越高。隨著冷卻水流速的增加，GMM棒的溫度逐漸降低，當(dāng)冷卻水流速大于2.5 m/s時，各電流值下的GMM棒溫度都小于25.2 ℃；當(dāng)冷卻水流速大于3 m/s時，GMM棒的溫度變化趨于平緩。

為了得到在采取強(qiáng)制水冷時，GMM棒的熱變形誤差，作電流強(qiáng)度為6 A、冷卻水流速為3.0 m/s下GMM棒軸向溫度分布曲線，如圖7所示。

圖7 GMM 棒軸向溫度分布

從圖7中可以看出，由于GMM棒軸向溫度存在差異，因此在計(jì)算熱變形量時需要對其進(jìn)行積分

式中：αl為GMM棒熱膨脹系數(shù)；l為GMM棒長度；Δt為GMM棒沿軸向的溫升。

根據(jù)圖7對曲線進(jìn)行反向擬合，通過式（19）可得GMM棒因溫度升高而引起的熱變形量為0.086 μm。

4 結(jié)論

1）針對內(nèi)管外腔式GMA水冷結(jié)構(gòu)確定了熱源功率，并基于熱阻定律建立了GMA的傳熱數(shù)學(xué)模型，分析了影響GMM棒溫度的因素。

2）GMA的溫度隨線圈電流強(qiáng)度的增大而增大，而通過強(qiáng)制對流換熱可以有效地控制GMA的溫升。

3）GMM棒溫升隨著冷卻水入口速度的增大而減小，當(dāng)冷卻水入口流速大于3 m/s時，GMM棒的升變化趨于平緩。當(dāng)電流為6 A時，冷卻水入口流速為 3 m/s時，GMM 棒的熱變形可控在 0.086 μm。