許焰 譚建平 王軍鋒 劉云龍 祝忠彥

(中南大學 機電工程學院，湖南 長沙 410083)

磁力傳動屬于耦合傳動，是以替代機械式傳動為目的，應用永磁材料或電磁機構所產(chǎn)生的磁力作用，實現(xiàn)力或力矩非接觸式傳遞的一種新技術［1-2］，因其可以實現(xiàn)非接觸的動力傳遞，近年來被廣泛應用于血泵驅(qū)動研究領域［3-5］.譚建平等［6］提出了一種大間隙磁力傳動方法，并將其應用于可植入式永磁軸流式血泵的外磁場驅(qū)動，實驗表明，在耦合間隙40 mm 的情況下，可輸出0.7 N·mm 的扭矩，驅(qū)動軸流式血泵達到5 000 r/min，證明了這種驅(qū)動方式的可行性，并為軸流式血泵的臨床應用提供了一個新的思路［7-11］.

然而，電磁能量的轉(zhuǎn)換與傳輸，必然伴隨著熱量的產(chǎn)生和傳遞，倘若不能使系統(tǒng)內(nèi)部熱量及時散發(fā)，會導致工作部件持續(xù)溫升，進而影響材料強度及系統(tǒng)性能，因此，為保證系統(tǒng)的可靠工作，必須對系統(tǒng)溫度場進行了分析計算.Yuan 等［12］對微型發(fā)電機的溫度場分布進行了分析;Higuchi 等［13］分析了圓柱導體中瞬態(tài)渦流場、溫度場的變化情況;王艷武等［14］研究了電機轉(zhuǎn)子的溫度場和熱應力場;丁樹業(yè)等［15］分析計算了大型發(fā)電機定子的三維溫度場;李志強等［16］研究了熱沖擊載荷下鋁蜂窩夾芯結(jié)構的溫度場及變形情況.

文中以軸流式血泵的大間隙磁力驅(qū)動系統(tǒng)為研究對象，對電磁體的溫度分布以及熱量損失、熱流密度、熱梯度等熱物理參數(shù)進行分析，通過有限元仿真計算，得到了系統(tǒng)溫度場的分布和變化規(guī)律，并對電磁驅(qū)動系統(tǒng)的電磁參數(shù)、運行環(huán)境參數(shù)等進行了優(yōu)化設計.

1 系統(tǒng)傳熱機制分析

對于大間隙磁力傳動系統(tǒng)，電磁體內(nèi)部的熱源主要來自線圈的銅損和鐵心的鐵損，通過熱傳導和熱對流兩種方式散發(fā).圖1 所示為“雙T”型電磁發(fā)生裝置熱交換機制的模型圖(文獻［9］詳細介紹了“雙T”型驅(qū)動裝置的磁傳動機理)，在無強制散熱措施時，線圈產(chǎn)生的熱量，一部分通過熱傳導，經(jīng)線圈和鐵心之間的骨架絕緣層傳給鐵心，再通過鐵心表面按對流換熱方式散發(fā)到空氣中;另一部分則直接靠對流換熱，通過暴露在空氣中的線圈表面?zhèn)髦量諝庵校?7］.在電磁體線圈通電開始一段時間內(nèi)，電磁體溫度分布隨時間變化，屬于瞬態(tài)導熱問題;連續(xù)工作一段時間后，其內(nèi)部溫度分布將保持穩(wěn)定，溫度場可視為穩(wěn)態(tài)場.

圖1 電磁體熱交換機制的模型Fig.1 Model of heat exchange mechanism of electromagnet

2 電磁體溫升的仿真計算

2.1 仿真條件及建模假設

線圈匝數(shù)均為650 匝、線徑為0.31 mm;單個線圈通電電流0.8 A，環(huán)境溫度為25 ℃，電磁體初始溫度為25 ℃，空氣流速為0 m/s.電磁體由多種材料組成，不是均勻發(fā)熱體，假設其厚度方向各個截面的溫度相等，從而將電磁體溫度場計算模型簡化為二維導熱問題.

2.2 仿真步驟

(1)模型的構建

設置單元類型:Type1 PLANE55，四邊形4 節(jié)點數(shù).劃分網(wǎng)格時設置單元大小為0.001，電磁體的結(jié)構參數(shù)如圖1 所示.材料屬性為:硅鋼片、線圈和線圈骨架的密度分別為7650.0、8429.3、1700.0kg/m，比熱容分別為446.0、395.7、1360.0 J/(kg·K).

(2)施加載荷及邊界條件

設置環(huán)境溫度為25 ℃;第三類邊界條件需要施加對流載荷于模型周圍，空氣流速為0 m/s，故定義模型邊界上的對流換熱系數(shù)為α=14.2W/(m2·K);給鐵心施加生熱率QFe，給線圈單元體施加生熱率QCu.

2.3 仿真結(jié)果及分析

用POST26 進行后處理，得到通電開始12 000 s時的電磁體溫度場分布，如圖2 所示.過高的溫度必然會對系統(tǒng)的長時間運行產(chǎn)生嚴重影響，為研究電磁發(fā)生裝置中的溫度場變化規(guī)律，選取圖2 中A、B、D 點為代表，查看電磁體通電開始12000 s 內(nèi)不同位置的溫度變化情況如圖3 所示.

圖2 電磁體穩(wěn)態(tài)溫度場分布Fig.2 Distribution of steady-state thermal field of electromagnet

圖3 電磁體3 個不同位置點溫度隨時間的變化曲線Fig.3 Change curves of temperature with time at three different positions of electromagnet

由圖2 可以看出:①電磁體溫度場分布左右對稱，單個電磁體溫度場分布也呈左右對稱;②當電磁體通電12000s 時，電磁體最低溫度出現(xiàn)在骨架邊緣(D 點附近)，為79.780 ℃;最高溫度出現(xiàn)在通電線圈(B、C 所示位置)，達到108.860 ℃.

由圖3 可以看出:①在前2400 s 內(nèi)，電磁體溫度上升較快;隨著時間的推移，電磁體溫度上升趨勢減弱并最終達到穩(wěn)定;②當t =9600 s 時，溫度最大值與穩(wěn)態(tài)時108.860 ℃相差不到1 ℃，故認為此時系統(tǒng)溫度場達到穩(wěn)定.

3 電磁體不同位置的溫度測定實驗

3.1 實驗條件及設備

基本實驗條件:電磁體鐵芯采用142 片厚度為0.35 mm 的硅鋼片壓制而成，四組線圈匝數(shù)均為650 匝，線徑為0.31mm，單個線圈通電電流為0.8A，保持室內(nèi)溫度25℃.利用DS18B20 數(shù)字溫度傳感器采集溫度值并發(fā)送到上位機，測溫范圍為-55 ℃到125 ℃.

為避免在實驗過程中燒壞線圈，本實驗中僅對線圈通電開始300 s 內(nèi)電磁體不同位置的溫度進行測量，測點A-F 位置布置同樣如圖2 所示.

3.2 實驗結(jié)果及分析

將各測點在線圈通電開始300 s 內(nèi)不同時刻的溫度仿真值及實驗數(shù)據(jù)擬合成曲線，得到各個測點位置的溫度值隨時間變化的曲線如圖4 所示.

圖4 各測點溫度隨時間的變化曲線Fig.4 Change curves of temperature with time at each measuring point

由圖4 可知:電磁體測試點溫度隨通電時間線性升高，仿真值與實驗值相比較，兩組曲線重合度較高，總體上是相符的，仿真結(jié)果是可信的.但由于幾何建模、實驗過程本身及材料屬性假設等原因產(chǎn)生了一定偏差.

4 降低電磁體溫升的措施分析

由于普通人體可承受的溫度為46～55 ℃，須保證電磁體最高溫度不超過50 ℃.通過對電磁體內(nèi)部傳熱機制的分析可知，通過減少電磁體損耗、改善電磁體導熱條件和散熱條件等，可以有效降低電磁體溫升.

(1)減少電磁體損耗

電磁體線圈通電電流是影響電磁體損耗的主要因素是，通過仿真可以得到，電磁體最高溫度Tm隨單個線圈通電電流I 大小變化的曲線如圖5 所示.

圖5 不同通電電流下的電磁體最高溫度Fig.5 Maximum temperature of electromagnet with different currents

由溫度隨通電電流的變化曲線，擬合出溫度和通電電流的函數(shù)關系式:

式中，θm為電磁體最高溫度，I 為單個線圈通電電流.要保證電磁體最高溫度θm不超過50℃，則電流I需要滿足I ＜0.56 A.

(2)改善電磁體導熱條件

由于電磁體線圈繞組導熱能力最差，因此主要考慮從增大線圈繞組的導熱系數(shù)方面來改善電磁體導熱條件.設置單個線圈通電電流為0.55 A，環(huán)境溫度為10 ℃.仿真得到不同導熱系數(shù)下的電磁體最高溫度變化情況如圖6 所示.可以看出，隨著線圈繞組導熱系數(shù)的增大，電磁體最高溫度逐漸減小，但是減小量很少.因此，通過改善電磁體導熱條件并不能很好地解決電磁體溫升問題.

(3)改善電磁體表面散熱條件

忽略散熱表面幾何尺寸等因素，則可近似認為電磁體表面的換熱系數(shù)α 僅決定于空氣的流速v，電流、環(huán)境溫度等參數(shù)設置同前，仿真得到不同空氣流速下的電磁體最高溫度數(shù)值和變化情況，如圖7所示.

圖6 不同導熱系數(shù)下電磁體最高溫度Fig.6 Maximum temperature of electromagnet with different thermal conductivities

圖7 不同對流換熱系數(shù)下電磁體最高溫度Fig.7 Maximum temperature of electromagnet with different convective heat transfer coefficients

可以看出，當空氣流速v =3.5 m/s 時，電磁體最高溫度θm=50.103 ℃，所以要保證電磁體最高溫度θm不超過50℃，則空氣流速v 需要滿足v ＞3.5m/s.

5 結(jié)論

通過對軸流式血泵的大間隙磁力驅(qū)動系統(tǒng)溫度分布以及熱量損失、熱流密度、熱梯度等熱物理參數(shù)進行分析，以及有限元仿真計算，得到了系統(tǒng)溫度場的分布和變化規(guī)律，電磁體最低溫度出現(xiàn)在骨架邊緣位置，為79.780 ℃;最高溫度出現(xiàn)在通電線圈中央附近，達到108.860℃;當通電時間t=170min 時，系統(tǒng)溫度場達到穩(wěn)定.

針對電磁體溫升過高的問題，通過對不同條件下電磁體穩(wěn)態(tài)溫度場進行仿真分析可知:減少電磁體損耗和改善電磁體散熱條件能夠有效降低電磁體溫升，當電磁體單個線圈通電電流I ＜0.56 A、空氣流速v ＞3.5 m/s 時，可以保證電磁體最高溫度不超過50 ℃，為降低電磁體溫升、保證系統(tǒng)穩(wěn)定可靠地工作提供了理論依據(jù).

［1］Alscher H，Iguchi M，Eastham A R，et al.Non-contact suspension and propulsion technology［J］.Vehicle System Dynamics，1983，12(4/5):259-289.

［2］徐先懂，龔中良，譚建平.基于外磁場耦合的血泵驅(qū)動系統(tǒng)［J］.中南大學學報:自然科學版，2007，38(4):711-714.Xu Xian-dong，Gong Zhong-liang，Tan Jian-ping.Blood pump driven system based on extracorporeal magnetic field couple［J］.Journal of Central South University:Science and Technology，2007，38(4):711-714.

［3］Ohashi T，Kuyama H，Suzuki K，et al.Control of aqueous droplets using magnetic and electrostatic forces［J］.Analytica Chimica Acta，2008，612(2):218-225.

［4］Ren Z H，Jahanmir S，Heshmat H.Design analysis and performance assessment of hybrid magnetic bearings for a rotary centrifugal blood pump［J］.ASAIO Journal，2009，55(4):340-347.

［5］Hoshi H，Katakoa K，Ohuchi K，et al.Magnetically suspended centrifugal blood pump with a radial magnetic driver［J］.ASAIO Journal，2005，51(1):60-64.

［6］譚建平，許焰，廖平，等.一種非接觸式大間隙磁力驅(qū)動方法:中國，200810030545.1［P］.2008-01-25.

［7］Tan Jian-ping，Liu Yun-long，Xu Yan，et al.Dynamic characteristics of a large gap magnetic driving blood pump during start-up process ［J］.Magnetohydrodynamics，2011，47(3):283-294.

［8］Tan Jian-ping，Liu Yun-long，Xu Yan，et al.Study on phase angle of blood pump control system driven by large gap magnetic field［J］.Advanced Science Letters，2011，4(4/5):1357-1360.

［9］劉云龍，譚建平，許焰，等.大間隙磁力傳動系統(tǒng)能量傳遞效率［J］.中南大學學報:自然科學版.2011，42(5):1296-1302.Liu Yun-long，Tan Jian-ping，Xu Yan，et al.Energy transfer efficiency of large gap magnetic drive system ［J］.Journal of Central South University:Science and Technology，2011，42(5):1296-1302.

［10］許焰，譚建平，李譚喜，等.行波磁場驅(qū)動的磁力傳動系統(tǒng)空間磁場數(shù)學模型［J］.中南大學學報:自然科學版.2009，40(6):1573-1579.Xu Yan，Tan Jian-ping，Li Tan-xi，et al.Space magneticfield's mathematical model of magnetic drive system driven by traveling wave magnetic field ［J］.Journal of Central South University:Science and Technology，2009，40(6):1573-1579.

［11］劉云龍，譚建平，祝忠彥，等.大間隙磁力傳動系統(tǒng)電磁場切換相位角研究［J］.湖南大學學報:自然科學版.2012，39(2):41-46.Liu Yun-long，Tan Jian-ping，Zhu Zhong-yan，et al.Phase angle of electromagnetic field in large gap magnetic drive system［J］.Journal of Hunan University:Natural Sciences，2012，39(2):41-46.

［12］Yuan Chun，Wei Ke-yin，Wang Xiao-nian.Coupled-field thermal analysis of high-speed permanent magnetic generator applied in micro-turbine generator［C］∥Proceedings of the Eighth International Conference on Electrical Machines and Systems.Nanjing:IEEE，2005:2458-2461.

［13］Higuchi M，Kawamura R，Tanigawa Y.Magneto-thermoelastic stresses induced by a transient magnetic field in a conducting solid circular cylinder ［J］.International Journal of Solids and Structures，2007，44(16):5316-5335.

［14］王艷武，楊立，陳翾，等.異步電機轉(zhuǎn)子三維溫度場及熱應力場研究［J］.電機與控制學報，2010，14(6):27-32.Wang Yan-wu，Yang Li，Chen Xuan，et al.Study on 3D thermal field and thermal stress field of the induction motor rotor［J］.Electric Machines and Control，2010，14(6):27-32.

［15］丁樹業(yè)，孫兆瓊，苗立杰，等.大型發(fā)電機定子主絕緣溫度場數(shù)值研究［J］.電機與控制學報，2010，14(7):53-58.Ding Shu-ye，Sun Zhao-qiong，Miao Li-jie，et al.Numerical investigation of stator main insulation temperature field for large generator［J］.Electric Machines and Control，2010，14(7):53-58.

［16］李志強，王志華，趙隆茂.熱沖擊載荷下鋁蜂窩夾芯結(jié)構的溫度場及變形［J］.華南理工大學學報:自然科學版，2011，39(9):97-102，127.Li Zhi-qiang，Wang Zhi-hua，Zhao Long-mao.Temperature field and deformation of sandwich structure with aluminum honeycomb cores under thermal loading ［J］.Journal of South China University of Technology:Natural Science Edition，2011，39(9):97-102，127.

［17］祝忠彥.血泵大間隙磁力傳動系統(tǒng)磁力矩相位角及電磁體溫升研究［D］.長沙:中南大學機電工程學院，2011:40-51.