程志遠(yuǎn) 宋曉逸 吳曉婷 貴子航 李東東

無線充電系統(tǒng)旋轉(zhuǎn)式電磁耦合器損耗計算及熱點溫度研究

程志遠(yuǎn) 宋曉逸 吳曉婷 貴子航 李東東

（上海電力大學(xué)電氣工程學(xué)院 上海 200090）

電磁耦合器是實現(xiàn)無線電能傳輸（WPT）技術(shù)的核心元件，其工作溫度決定了系統(tǒng)運行的穩(wěn)定性與使用壽命。為此該文以抗偏移能力強的旋轉(zhuǎn)式電磁耦合器為研究對象，針對該電磁耦合器在非正弦激勵下?lián)p耗及溫度精確計算的問題，同時為避免大量實驗數(shù)據(jù)繁瑣采集過程，利用瞬態(tài)電磁-熱流場雙向間接耦合方法，建立考慮絕緣引起的匝間距及溫度對材料屬性影響的三維電磁場仿真模型，而后將電磁損耗計算結(jié)果作為載荷耦合至熱流場中并建立對應(yīng)關(guān)系，實現(xiàn)對熱點溫度計算。仿真與實驗結(jié)果表明，物理模型的熱點溫度與測量值最大誤差為5.15%，仿真得出的溫升分布規(guī)律能較準(zhǔn)確地展示實際情況，驗證了物理仿真模型的合理性，可為旋轉(zhuǎn)式電磁耦合器的散熱設(shè)計提供依據(jù)，在工程應(yīng)用中具有一定指導(dǎo)意義。

無線電能傳輸 旋轉(zhuǎn)式電磁耦合器 多物理場耦合 溫升特性 電磁損耗

0 引言

自主式水下航行器（Autonomous Underwater Vehicles, AUV）是目前人類探索海洋資源的重要設(shè)備之一，而采用有線輸電方式是限制其應(yīng)用范圍和發(fā)展的關(guān)鍵問題[1-2]。無線電能傳輸（Wireless Power Transfer, WPT）技術(shù)采用非物理接觸傳能方式，同時具有高自動化程度以及安全性、強靈活性和適用性等特點，可作為AUV可靠的供電方式[3-5]。典型的磁耦合諧振式WPT系統(tǒng)示意圖如圖1所示，其中電磁耦合器是實現(xiàn)WPT的核心組件，用于進(jìn)行原、副邊間“電-磁-電”能量轉(zhuǎn)換傳輸。

圖1 磁耦合諧振式WTP系統(tǒng)示意圖

為減小暗流及水壓對電磁耦合器原、副邊對位的影響，AUV大多采用籠狀對接系統(tǒng)。本文設(shè)計的籠狀對接無線充電裝置如圖2所示。其中電磁耦合器原邊安裝于籠狀對接裝置中，副邊安裝于AUV機身處。隨著系統(tǒng)頻率和容量的不斷提升、整機體積的不斷縮小，電磁耦合器的損耗和溫升問題逐漸明顯，其溫升過高不僅會加快整機系統(tǒng)的絕緣老化速度、增大組件的熱應(yīng)力，影響系統(tǒng)使用壽命，同時系統(tǒng)中電子器件的參數(shù)也會因為周圍溫度變化發(fā)生漂移，進(jìn)一步影響系統(tǒng)效率和運行穩(wěn)定性[6-9]。因此，精確預(yù)估高頻電磁耦合器線圈繞組及磁心損耗、研究正常模態(tài)下電磁耦合器溫升分布，對AUV的散熱設(shè)計及提高WPT系統(tǒng)的可靠性至關(guān)重要。

圖2 AUV對接系統(tǒng)示意圖

為減小高頻環(huán)境下電磁耦合器的磁心損耗，WPT系統(tǒng)通常采用磁滯回線狹長、磁導(dǎo)率較高的軟磁材料，同時利用多股絞合Litz線繞制線圈繞組來降低由趨膚效應(yīng)引起的渦流損耗[10-11]。將線圈繞組簡化為整塊體積相同的導(dǎo)體是目前常用的建模方法，該方法能極大降低建模難度、縮短計算時間，但在高頻電磁場中的計算精度往往難以保證。文獻(xiàn)[12-13, 16-18]均采用簡化模型，忽略由于絕緣層厚度引起的匝間距對電磁耦合器電阻等參數(shù)的影響。文獻(xiàn)[12]使用頻域求解器對軸式松耦合變壓器進(jìn)行電磁損耗計算，忽略經(jīng)高頻逆變器之后高次諧波對損耗計算結(jié)果的影響；文獻(xiàn)[13]采用高磁導(dǎo)率、低損耗的磁條增大原、副邊之間耦合系數(shù)，但在計算損耗時默認(rèn)Litz線圈交流電阻等于直流電阻。目前，電磁耦合器溫升計算方法主要包括等效熱網(wǎng)絡(luò)法和數(shù)值計算法。等效熱網(wǎng)絡(luò)法通常是根據(jù)研究對象的熱流路徑建立熱回路模型，對各組件發(fā)熱損耗功率進(jìn)行計算及設(shè)置。如果建立正確的熱回路模型并設(shè)置合理的功率添加方法，電磁耦合器的溫升分布規(guī)律就比較容易得到。文獻(xiàn)[14]研究空氣中磁耦合器的熱回路模型和加熱機理，計算各部件發(fā)熱功率并定性地給出溫度分布，實驗結(jié)果與理論分析吻合度較高；文獻(xiàn)[15]為提高熱模型精度，利用仿真與實驗相結(jié)合的方法確定研究對象的表面?zhèn)鳠嵯禂?shù)，結(jié)果表明，根據(jù)實驗測算建立的熱模型有足夠高的工程精度。由于等效熱網(wǎng)絡(luò)法基于大量實驗數(shù)據(jù)測算，同時建模過程繁瑣，因此不適合多參數(shù)復(fù)雜模型計算。有限元法考慮磁心材料的非線性及導(dǎo)體各向異性參數(shù)的影響，可極大地提高計算精度且能夠直觀地得到空間內(nèi)磁場與溫升分布，在工程計算中逐漸得到應(yīng)用。文獻(xiàn)[16]對通入相反電流的兩個銅導(dǎo)體提出電-磁-熱分析模型，但該結(jié)果僅適用于二維模型計算；文獻(xiàn)[17]提出一種精確測量磁耦合器損耗的實驗方法，但僅在特定溫度下進(jìn)行測量對比，并且缺少在系統(tǒng)運行時磁耦合器各部件溫度分布；文獻(xiàn)[18]通過多物理場耦合研究了工作頻率及溫度對車用平板式電磁耦合器材料特性的影響規(guī)律。雖然在文章中討論模型各部分的發(fā)熱量，但并未給出各部件的溫升分布。

針對上述問題，本文首先建立了考慮旋轉(zhuǎn)式電磁耦合器線圈絕緣層的三維電磁場仿真模型，同時考慮電磁材料屬性的溫度效應(yīng)，基于三維瞬態(tài)場仿真得到電磁組件損耗非均勻分布。采用電磁-熱流場間接耦合計算方法求解自然對流下電磁耦合器的穩(wěn)態(tài)溫升分布，最終搭建一套帶有紅外熱成像儀的WPT系統(tǒng)樣機對電磁耦合器進(jìn)行溫升測試，驗證仿真模型的合理性與準(zhǔn)確性。

1 電磁場計算模型

結(jié)合籠狀對接無線充電系統(tǒng)外形，為減小由于對位偏差對系統(tǒng)傳輸功率帶來的影響，本文研究對象為環(huán)型磁心結(jié)構(gòu)的旋轉(zhuǎn)式電磁耦合器，如圖3所示。其主要由環(huán)形磁心及多芯Litz線圈繞組構(gòu)成，原邊線圈緊貼原邊環(huán)形磁心的內(nèi)側(cè)，副邊線圈緊貼副邊環(huán)形磁心外側(cè)。副邊線圈嵌套于原邊線圈內(nèi)，原、副邊磁心由TDK PC_40鐵氧體構(gòu)成。為使得原邊線圈與磁心緊密貼合，利用3D打印定制原邊緊固件?？紤]到用于磁心與線圈間的緊固件、墊板等部件由樹脂材料構(gòu)成，對磁場分布基本沒有影響，故在電磁場中建模時可忽略。實驗樣機具體參數(shù)見表1。

圖3 旋轉(zhuǎn)式電磁耦合器示意圖

表1 實驗樣機具體參數(shù)

Tab.1 Specific parameters of experimental prototype

針對旋轉(zhuǎn)式電磁耦合器電磁損耗精確計算問題，本文考慮了由絕緣引起的匝間距對計算電磁耦合器高頻交流電阻的影響，同時考慮到多芯Litz線難以精細(xì)化建模，根據(jù)等效電阻不變的原則，利用孔隙率系數(shù)將多芯Litz線圈簡化為截面積相同的方形導(dǎo)體[19]。

2 多物理場耦合計算原理

圖4為求解旋轉(zhuǎn)式電磁耦合器三維磁-熱-流耦合計算流程。首先對三維電磁模型設(shè)置與溫度相關(guān)的材料屬性，在初始溫度下計算電磁耦合器的損耗密度；其次采用間接耦合方式將電磁損耗作為熱流場激勵，利用網(wǎng)格-節(jié)點單元對應(yīng)關(guān)系對流熱場模型的流體、傳熱及能量方程進(jìn)行耦合分析，計算出三維模型的表面?zhèn)鳠嵯禂?shù)；然后，根據(jù)溫度結(jié)果反饋至電磁場中得到對應(yīng)的材料系數(shù)后，再反復(fù)迭代計算；最后，當(dāng)?shù)綌?shù)內(nèi)監(jiān)測點的相鄰兩步間溫度差值小于1 K時，終止迭代并得到電磁耦合器穩(wěn)態(tài)溫升分布。

圖4 多場耦合計算流程

2.1 瞬態(tài)電磁場計算

2.2 電磁損耗計算

電磁耦合器的電磁損耗是產(chǎn)生溫升的主要原因，其主要包括線圈損耗和磁心損耗。

2.2.1 線圈損耗計算

線圈損耗是由于線圈中通入電流引起的損耗，通常是電磁耦合器損耗產(chǎn)生的主要來源，高頻情況下可分解為趨膚效應(yīng)損耗和鄰近效應(yīng)損耗。為減小趨膚效應(yīng)損耗，WPT系統(tǒng)常采用多芯絞合型Litz線繞制線圈。根據(jù)Litz線的空間結(jié)構(gòu)分析可知，當(dāng)流經(jīng)頻率較高的電流時，單匝Litz線電阻con包括因自身電流產(chǎn)生的趨膚效應(yīng)電阻skin和組成Litz線的細(xì)線間由于鄰近效應(yīng)產(chǎn)生的內(nèi)部感應(yīng)電阻in_pro。WPT中諧振線圈的交流電阻還應(yīng)包括由相鄰其他線匝流經(jīng)交流電，使本線匝中產(chǎn)生的匝間感應(yīng)電阻turn_pro，以及系統(tǒng)中其他線圈流經(jīng)交流電產(chǎn)生的交變磁場在本線圈中產(chǎn)生的線圈間感應(yīng)電阻coil_pro。因此單個諧振線圈的總電阻coil的表達(dá)式為

假設(shè)除了該待求線圈外，系統(tǒng)中還存在個其他線圈，同時認(rèn)為Litz線中每股細(xì)線導(dǎo)通電流及物理參數(shù)均相等，忽略線圈彎曲部分的內(nèi)外半徑延展導(dǎo)致的輕微變化，根據(jù)文獻(xiàn)[20]式（3），skin、in_pro、turn_pro和coil_pro分別為

當(dāng)系統(tǒng)工作在諧振點附近且系統(tǒng)中僅存在一對諧振線圈時，副邊線圈中感應(yīng)電流與原邊線圈中的發(fā)射電流的相位差為p/2，故可推導(dǎo)出諧振線圈間的感應(yīng)耦合電阻為零。至此，在本文研究的系統(tǒng)中單個線圈的損耗計算公式為

2.2.2 磁心損耗計算

為提高WPT系統(tǒng)的傳輸能力，通常會在電磁耦合器中加入鐵氧體等磁心材料用以配合諧振線圈來提高系統(tǒng)互感，但與此同時也會引入磁心損耗。磁心損耗通常與磁心的物理參數(shù)、工作頻率以及磁感應(yīng)強度等相關(guān)，可用帶溫度系數(shù)的修正Steinmetz公式進(jìn)行計算。

2.3 熱流場計算方程

由于電磁耦合器的磁心、線圈、緊固件表面與周圍空氣間均存在溫度差，熱量會從溫度高的部分流向溫度低的部分以達(dá)到熱穩(wěn)定。在自然對流冷卻方式下，電磁耦合器的傳熱方式主要包括熱傳導(dǎo)、熱對流及熱輻射三種方式?？紤]到原邊線圈與原邊磁心、副邊線圈與副邊磁心均為緊密貼合方式，因此對于原邊磁心、線圈和緊固件，以及副邊磁心和線圈接觸面間的穩(wěn)態(tài)熱傳導(dǎo)過程可根據(jù)傅里葉定律進(jìn)行描述為

電磁耦合器的熱量來源于磁心以及線圈產(chǎn)生的電磁損耗，溫度升高致使其周圍空氣介質(zhì)的密度減小。根據(jù)阿基米德定律可知，電磁耦合器中絕大部分熱量通過周圍空氣循環(huán)散失。考慮到三維熱流場模型中的NS方程具有高度非線性且相互作用非常復(fù)雜，故在大范圍內(nèi)對空氣進(jìn)行自然對流散熱計算時，動量守恒方程可利用Boussinesq近似簡化分析。在開放空間的流體場內(nèi)，對電磁耦合器自然對流傳熱過程建立控制方程為

流體場中質(zhì)量守恒、動量守恒及能量守恒方程分別為

式中，、、分別為流體速度在、、軸方向的分量；為流體的體積密度；為流體的動力粘度；為流體受到的壓力；為流體的熱膨脹系數(shù)；、、分別為體積力在、、軸方向的加速度；為比定壓熱容；為流體介質(zhì)導(dǎo)熱系數(shù)；c為流體的內(nèi)熱源，即在忽略黏性耗散、輻射時，損耗熱量流向空氣中的部分。

在電磁耦合器中，緊固件與副邊線圈絕緣、副邊磁心外側(cè)無線圈接觸部分之間存在輻射換熱過程。本文基于面-面輻射方式，利用Stefan-Boltzmann方程對該過程進(jìn)行描述為

其中，面等效發(fā)射率r為

式中，該范圍邊界溫度設(shè)為常數(shù)；該范圍流體邊界設(shè)置為無滑移條件，對照電磁耦合器應(yīng)用位置而言，氣流速度在、方向無加速度，僅存在沿軸負(fù)方向的重力加速度。

3 計算結(jié)果與分析

3.1 電磁損耗計算分析

本文采用Matlab Simulink搭建WPT系統(tǒng)中DC-DC變換電路，模擬移相全橋控制下電流波形?？紤]到旋轉(zhuǎn)式電磁耦合器電流波形為非正弦，因此選擇Ansys Maxwell瞬態(tài)場求解器求解電磁損耗。得到的電磁耦合器磁場分布情況如圖5所示。

圖5 電磁耦合器磁場分布

觀察圖5可以發(fā)現(xiàn)，電磁耦合器的磁場分布基本呈上下對稱，磁場強度沿磁心中心位置向端部逐漸衰減。副邊磁心中間位置磁力線分布最密，最大磁感應(yīng)強度位于貼近副邊線圈幾何中心的磁心表面處，約為5.32×10-2T。根據(jù)文獻(xiàn)[21]可得到旋轉(zhuǎn)式電磁耦合器最大磁感應(yīng)強度，在此本文不再過多贅述。最大磁感應(yīng)強度計算公式為

式中，為線圈匝數(shù)；為原邊磁心內(nèi)邊緣半徑；為副邊磁心外邊緣半徑；為線圈高度。代入相應(yīng)數(shù)據(jù)可計算出旋轉(zhuǎn)式電磁耦合器最大磁感應(yīng)強度max=5.18×10-2T，與瞬態(tài)電磁場所得最大磁感應(yīng)強度數(shù)值的差值控制在誤差范圍內(nèi)。

圖6為電磁耦合器在一周期內(nèi)損耗密度的分布情況。通過觀察圖5和圖6a可以發(fā)現(xiàn)，旋轉(zhuǎn)式電磁耦合器磁心的平均損耗密度分布與磁感應(yīng)強度分布具有相同特性。副邊磁心靠近中間位置在一周期內(nèi)損耗最大，原、副邊磁心平均損耗密度沿中心位置向端部衰減，副邊磁心單位面積的平均損耗密度高于原邊磁心。從圖6b可以看出，副邊線圈的平均損耗密度高于原邊線圈，線圈的平均損耗密度coil計算公式為

式中，為電磁耦合器工作的一周期；i(t)為線圈內(nèi)流經(jīng)的總電流；Rcoil可通過式（3）計算得到。根據(jù)式（19）計算得到20℃時原、副邊線圈在一周期內(nèi)的平均損耗分別為0.381 W、3.133 W，有限元仿真得到原、副邊線圈在一周期內(nèi)的平均損耗分別為0.392 W、3.241 W，誤差分別為2.9%、3.4%，驗證了本文電磁場三維模型的準(zhǔn)確性。

3.2 電磁耦合器熱流場計算模型

電磁耦合器中每個部件產(chǎn)生的電磁損耗將會以熱能的形式消散，并最終表現(xiàn)為溫度升高。相較于電磁場中三維模型，線圈、磁心表面間的絕緣材料以及用于固定線圈和磁心的緊固件對熱流場分析影響較大，因此應(yīng)考慮在內(nèi)。求解電磁耦合器的三維熱流場模型如圖7所示，其中忽略基座和緊固件中對散熱影響較小的部分。

圖7 熱流場模型

本文根據(jù)文獻(xiàn)[22]中的變壓器溫度場模型建立方法對單層多匝線圈電磁耦合器熱流場模型進(jìn)行簡化，如圖8所示。對Litz線單股外部絕緣及單根外側(cè)尼龍絲一起等效，根據(jù)等效導(dǎo)熱系數(shù)計算公式

式中，eq為等效導(dǎo)熱系數(shù)；S和分別為各絕緣材料截面積及相應(yīng)的導(dǎo)熱系數(shù)。電磁耦合器中磁心、銅線、環(huán)氧改性漆、復(fù)合聚氨酯漆、PLA塑膠和尼龍絲的導(dǎo)熱系數(shù)分別為27.2、387.6、0.21、0.23、0.226和0.25 W/(m·℃)。

圖8 線圈等效模型

3.3 穩(wěn)態(tài)熱流場計算結(jié)果分析

將電磁損耗計算結(jié)果間接耦合至熱流場中作為初始熱源，模擬自然對流情況下旋轉(zhuǎn)式電磁耦合器達(dá)到穩(wěn)態(tài)時的溫升變化，以實驗測試環(huán)境為參考值，設(shè)置初始環(huán)境溫度為20℃，空氣相對壓強為0 Pa。為減少熱流場迭代步數(shù)和增強計算收斂性，在重力方向的入口速度設(shè)置為0.005 m/s。

圖9給出了電磁耦合器穩(wěn)態(tài)溫升與流速分布結(jié)果，可以看出旋轉(zhuǎn)式電磁耦合器溫升隨縱向高度增加呈上升趨勢。外部氣體主要從軸方向進(jìn)入，由于受電磁耦合器嵌套式結(jié)構(gòu)影響，原、副邊間氣道內(nèi)空氣徑向流通路徑受阻，導(dǎo)致氣流流速較低，降低對流換熱速率，使內(nèi)側(cè)磁心及線圈溫度較高。從圖9a中可以看出，線圈溫度整體高于磁心溫度，這是由于線圈的發(fā)熱功率遠(yuǎn)高于磁心。每匝線圈作為獨立發(fā)熱源而言，由于外側(cè)環(huán)繞線圈的散熱較好，靠近中間位置的線圈散熱受到外側(cè)線圈熱源的影響使得溫度略高。線圈緊密環(huán)繞在磁心表面，并且在建模時已經(jīng)假設(shè)線圈與磁心接觸良好，結(jié)合圖9b和圖9c可以觀察到，越靠近線圈、周圍空氣流速越高的位置，溫度越高；與線圈距離較遠(yuǎn)、周圍空氣流速低的位置，溫度越低，磁心的溫升分布特點符合熱力學(xué)定律。同時，由于緊固件的存在使得原邊磁心和原邊線圈增大了與緊固件之間的傳導(dǎo)熱源，降低和消除了與空氣之間的對流熱源。緊固件作為非導(dǎo)磁性材料不產(chǎn)生能量損耗，因此其發(fā)熱功率為零，圖9a中所示緊固件的溫升主要來源于原邊線圈和原邊磁心的熱傳導(dǎo)。

圖9 電磁耦合器溫度與流速分布結(jié)果

旋轉(zhuǎn)式電磁耦合器副邊磁心溫度分布如圖10所示，其熱點溫度出現(xiàn)在副邊線圈頂部對應(yīng)區(qū)域，相比環(huán)境溫度的最高溫升大約為21℃；最低溫度位于磁心邊緣處，與副邊線圈軸向距離最小值為 1.3 mm，溫差大約為3℃，平均溫度為39.6℃。副邊磁心最高溫升產(chǎn)生的原因一部分來自自身的鐵心，另一部分來自于副邊高溫線圈的熱傳導(dǎo)。溫度不均勻分布對電磁耦合器參數(shù)計算增加不確定性，同時長期過熱容易引發(fā)磁心變形或損壞。因此當(dāng)旋轉(zhuǎn)式電磁耦合器應(yīng)用于大功率情況時，在研發(fā)設(shè)計初期應(yīng)考慮副邊磁心的溫升情況。

圖10 副邊磁心溫升分布

3.4 溫升測試結(jié)果

為驗證旋轉(zhuǎn)式電磁耦合器多物理場耦合仿真模型的合理性與仿真結(jié)果的準(zhǔn)確性，本文以環(huán)境溫度為20℃，諧振頻率為50 kHz，最大穩(wěn)定輸出功率為1 kW的無線電能傳輸系統(tǒng)電磁耦合器進(jìn)行溫度分布的實驗測試。搭建圖11所示的無線電能傳輸系統(tǒng)實驗平臺，利用LW-1001KD為系統(tǒng)提供工頻整流后的100 V直流電壓輸入，通過RIGOL DS7024示波器采集原、副邊線圈輸出電壓、電流波形，同時JT6343A作為直流負(fù)載，同步顯示系統(tǒng)輸出功率。采用紅外熱成像儀實時觀測電磁耦合器整體的溫升變化過程，配合使用K型熱電偶多通道測溫儀和紅外測溫槍可更加直觀地監(jiān)測在仿真中出現(xiàn)的局部熱點。

圖11 WPT系統(tǒng)實驗平臺

通過LCR測試儀測量旋轉(zhuǎn)式電磁耦合器參數(shù)見表2，其中原、副邊線圈自感分別為161.2mH、152.4mH，互感為133.1mH，電阻值分別為70.9 mW、64 mW，與電磁仿真中數(shù)值偏差均小于10%，驗證了建模方法的有效性。其中存在誤差的主要原因有：①在仿真建模中忽略Litz線內(nèi)部絞合結(jié)構(gòu)，采用簡化模型會使得仿真值與實際測量值產(chǎn)生偏差；②物理模型中材料屬性默認(rèn)為各向同性對參數(shù)的計算也會產(chǎn)生影響。

表2 旋轉(zhuǎn)式電磁耦合器仿真與測量參數(shù)比較

Tab.2 Comparison of simulation and measurement parameters of rotary electromagnetis coupler

圖12為功率分析儀采集到的數(shù)據(jù)，由圖12可知，輸入功率dc為645 W時，輸出功率L為525 W，且原邊電流dc有效值約為6.34 A，副邊電流L有效值約為 16.18 A，其有效值約為原邊電流有效值的 2.5倍。根據(jù)損耗計算公式可知，旋轉(zhuǎn)式電磁耦合器副邊產(chǎn)生的損耗將大于原邊，即旋轉(zhuǎn)式電磁耦合器副邊產(chǎn)生的發(fā)熱功率將大于原邊。

圖12 功率分析儀采集數(shù)據(jù)

根據(jù)3.3節(jié)熱流場模擬結(jié)果選擇各部件的溫度測量位置，采用多路測溫儀對選擇熱點進(jìn)行溫度測量。溫升測試的結(jié)果可以明確系統(tǒng)中發(fā)熱最為嚴(yán)重的部分，從而在系統(tǒng)優(yōu)化時多加考慮，因此選取仿真中各部件的最高溫度點。分別選擇原邊磁心外側(cè)頂端中心位置，副邊磁心內(nèi)側(cè)頂端中心位置，原、副邊線圈軸向中心位置，即從外向內(nèi)的第11、12匝線圈間作為溫度測量點，如圖13所示。

圖13 松耦合變壓器各部件測溫位置

使用FLUKE紅外熱成像儀拍攝電磁耦合器的穩(wěn)態(tài)溫升分布，如圖14所示。熱電偶多路測溫儀記錄了在諧振狀態(tài)下的系統(tǒng)運行20 min后的實驗值，測量結(jié)果與仿真值大致相同，見表3。造成溫度誤差的主要原因有：

（1）為使銅線圈的熱源體積密度相同，將其簡化為與銅材料體積相同的薄板式結(jié)構(gòu)，忽略每匝繞線間的尼龍纖維對線圈的散熱影響。

（2）三維熱流場仿真中模擬的自然對流屬于理想狀態(tài)，實際狀況下的實驗平臺無法完全滿足條件。仿真溫度與實測值偏差≤6%，滿足工程計算需要，驗證了建模方法的合理性。

圖14 紅外線穩(wěn)態(tài)溫升分布結(jié)果

表3 熱點溫度仿真與實驗結(jié)果

Tab.3 Hot spot temperature simulation and experimental results

4 結(jié)論

鑒于電磁耦合器的工作溫度是決定WPT系統(tǒng)穩(wěn)定運行的關(guān)鍵因素，本文建立了考慮絕緣引起匝間距的旋轉(zhuǎn)式電磁耦合器三維電磁場仿真模型，并對其磁場與損耗進(jìn)行計算和分析。采用磁-熱-流多物理場耦合方法計算自然狀態(tài)下表面?zhèn)鳠嵯禂?shù)，實現(xiàn)對電磁耦合器的溫升計算，最終通過所搭建的實驗平臺驗證理論分析的正確性，可得到以下結(jié)論：

1）旋轉(zhuǎn)式電磁耦合器磁心的平均損耗密度與磁感應(yīng)強度的分布呈正相關(guān)，靠近線圈軸向中心位置的磁感應(yīng)強度最高，損耗密度最大，端部最小，基本呈遞減規(guī)律分布。

2）旋轉(zhuǎn)式電磁耦合器高溫區(qū)域位于副邊線圈頂端，最高溫度仿真值與測量值分別為40.7℃和42.3℃，誤差為3.93%，副邊線圈自身最大溫差約6.6℃。由于磁心的存在使得電磁耦合器極大地減小了漏磁，在提高傳輸效率時也引入了損耗，受制于嵌套結(jié)構(gòu)的影響，在線圈高溫傳導(dǎo)下產(chǎn)生局部過熱，其最高溫度相比于環(huán)境溫度高22.3℃。

3）考慮材料屬性溫度效應(yīng)的多物理場耦合計算，克服傳統(tǒng)磁路法計算不準(zhǔn)確的問題，能夠更好地研究旋轉(zhuǎn)式電磁耦合器的溫升分布。通過對電磁耦合器的溫升測試，驗證了建模方法的合理性與結(jié)果的準(zhǔn)確性。

本文的分析計算可為無線充電系統(tǒng)旋轉(zhuǎn)式電磁耦合器的損耗計算及溫升分布研究提供理論依據(jù)。后續(xù)將開展對旋轉(zhuǎn)式電磁耦合器應(yīng)用于水下大功率下的溫升試驗，以及針對高溫部分展開散熱優(yōu)化研究?；诒疚哪Ｐ?，可以對AUV的傳輸特性隨溫度的變化趨勢展開深入分析，本文的多物理場耦合分析也為AUV的結(jié)構(gòu)設(shè)計提供思路。

[1] Zeng Yingqin. Lu Conghui, Liu Renzhe, et al. Liu, Wireless power and data transfer system using multidirectional magnetic coupler for swarm AUVs[J]. IEEE Transactions on Power Electronics, 2023, 38(2): 1440-1444.

[2] 閆爭超, 胡謙宇, 趙晨旭, 等. 水下航行器感應(yīng)式無線電能傳輸技術(shù)研究綜述[J/OL]. 中國電機工程學(xué)報, 2022: 1-16[2023-02-15]. https://kns.cnki.net/ kns8/defaultresult/index. Yan Zhengchao, Hu Qianyu, Zhao Chenxu, et al. Summary of research on induction radio energy transmission technology for underwater vehicles[J/OL].Proceedings of the CSEE, 2022: 1-16[2023-02-15] https:// kns.cnki.net/kns8/defaultresult/index

[3] Yan Zhengchao, Zhang Yiming, Zhang Kehan, et al. Underwater wireless power transfer system with a curly coil structure for AUVs[J]. IET Power Electronics, 2019, 12(10): 2559-2565.

[4] Qiao Kangheng, Sun Pan, Rong Enguo, et al. Anti-misalignment and lightweight magnetic coupler with H‐shaped receiver structure for AUV wireless power transfer[J]. IET Power Electronics, 2022, 15(16): 1843-1857.

[5] Gerasimov V A,Komlev A V, Filozhenko A Y. improving the quality of the underwater robot’s contactless battery charging system[C]// Journal of Physics: Conference Series, 2021, 2096(1): 012107-. DOI:10.1088/1742-6596/2096/1/012107.

[6] Chae D J, Kim Y S. Analysis for electrical and thermal hazard in kW-class wireless power transmission environment[C]// 2018 IEEE PELS Workshop on Emerging Technologies: Wireless Power Transfer (Wow), Montreal, QC, Canada, 2018. DOI:10.1109/ WoW.2018.8450930.

[7] Hwang K, Chung S, Yoon U, et al. Thermal analysis for temperature robust wireless power transfer systems[C]//2013 IEEE Wireless Power Transfer (WPT), Perugia, Italy, 2013.DOI:10.1109/WPT.2013. 6556879.

[8] Li Yang, Zhao Jingtai, Yang Qingxin, et al. A novel coil with high misalignment tolerance for wireless power transfer[J]. IEEE Transactions on Magnetics, 2019, 55(6): 1-4.

[9] Campi T, Cruciani S, De Santis V, et al. EMF safety and thermal aspects in a pacemaker equipped with a wireless power transfer system working at low frequency[J]. IEEE Transactions on Microwave Theory and Techniques, 2016, 64(2): 375-382.

[10] 陳慶彬, 范峰, 汪金帥, 等. 無線電能傳輸磁耦合系統(tǒng)Litz線圈交流電阻精確評估方法[J]. 電工技術(shù)學(xué)報, 2022, 37(24): 6294-6305. Chen Qingbin, Fan Feng, Wang Jinshuai, et al. Accurate evaluation method for AC resistance of Litz coil of radio energy transmission magnetic coupling system [J].Transactions of China Electrotechnical Society, 2022, 37(24): 6294-6305.

[11] 張獻(xiàn), 韓大穩(wěn), 沙琳, 等. 一種共享磁通多耦合模式的無線電能傳輸系統(tǒng)抗偏移方法[J]. 電工技術(shù)學(xué)報, 2022, 37(21): 5359-5368. Zhang Xian, Han Dawen, Sha Lin, et al. An anti-offset method for radio energy transmission system with shared magnetic flux multi-coupling mode[J] Transactions of China Electrotechnical Society, 2022, 37(21): 5359-5368.

[12] 趙蒙蒙, 李巖松, 豐江波, 等. 軸式松耦合變壓器的熱行為特性及影響因素[J]. 電測與儀表, 2019, 56(7): 22-29. Zhao Mengmeng, Li Yansong, Feng Jiangbo, et al. Thermal behavior characteristics and influence factors of shaft-type loosely coupled transformer[J]. Electrical Measurement and Instrumentation, 2019, 56(7): 22-29.

[13] Wang Shuo, Dorrell David G. Loss analysis of circular wireless EV charging coupler[J]. IEEE Transactions on Magnetics, 2014, 50(11): 1-4.

[14] Liang Ce,Yang Guang, Yuan Feng, et al. Modeling and analysis of thermal characteristics of magnetic coupler for wireless electric vehicle charging system[J]. IEEE Access, 2020, 8: 173177-173185.

[15] 毛行奎, 陳為. 開關(guān)電源高頻功率平面變壓器熱設(shè)計研究[J]. 電工電能新技術(shù), 2007, 26(4): 21-25. Mao Xingkui, Chen Wei. Research on thermal design of high frequency power plane transformer for switching power supply[J]. New Electrical Energy Technology, 2007, 26(4): 21-25.

[16] Zhu Tianxin, Feng Panshang, Li Xuekun, et al. The study of the effect of magnetic flux concentrator to the induction heating system using coupled electromagnetic-thermal simulation model[C]//International Conference on Mechanical and Automation Engineering, Jiujiang, China,2013: 123-127.

[17] Kalra G R, Pearce M G S, Kim S, et al. A power loss measurement technique for inductive power transfer magnetic couplers[J]. IEEE Journal of Emerging and Selected Topics in Industrial Electronics, 2020, 1(2): 113-122.

[18] Tiemann M, Saifo M, Clemens M, et al. Magnetic and thermal coupled field analysis of wireless charging systems for electric vehicles[J]. IEEE Transactions on Magnetics, 2019, 55(6): 1-4.

[19] Kazimierczuk M. High-frequency magnetic components[M]. 2nd ed. Chichester, UK: John Wiley & Sons, 2014.

[20] Liu Jiangtao, Deng Qijun, Czarkowski D, et al. Frequency optimization for inductive power transfer based on AC resistance evaluation in Litz-wire coil[J]. IEEE Transactions on Power Electronics, 2019, 34(3): 2355-2363.

[21] 程志遠(yuǎn), 陳坤, 李東東, 等. 旋轉(zhuǎn)式無線充電系統(tǒng)偏移特性研究[J]. 電工技術(shù)學(xué)報, 2021, 36(22): 4648-4657. Cheng Zhiyuan, Chen Kun, Li Dongdong, et al. Research on the offset characteristics of rotary wireless charging system[J]. Transactions of China Electrotechnical Society, 2021, 36(22): 4648-4657.

[22] 李永建, 閆鑫笑, 張長庚, 等. 基于磁-熱-流耦合模型的變壓器損耗計算和熱點預(yù)測[J]. 電工技術(shù)學(xué)報, 2020, 35(21): 4483-4491. Li Yongjian, Yan Xinxiao, Zhang Changgeng, et al. Transformer loss calculation and hot spot prediction based on magnetic-thermal-current coupling model[J]. Transactions of China Electrotechnical Society, 2020, 35(21): 4483-4491.

Loss Calculation and Hot Spot Temperature Research of Rotary Electromagnetic Coupler in Wireless Charging System

Cheng Zhiyuan Song Xiaoyi WuXiaoting Gui Zihang Li Dongdong

（School of Electric Power Engineering Shanghai University of Electric Power Shanghai 200090 China）

Electromagnetic couplers are the core components of wireless power transfer (WPT), and their operating temperature determines the stability and service life of the system operation. Therefore, accurately estimating the coil winding and magnetic core losses of high-frequency electromagnetic couplers and studying the temperature rise distribution of electromagnetic couplers under normal modes are crucial for the heat dissipation design of underwater autonomous vehicles and improving the reliability of WPT systems. To address these issues, this article takes the rotary electromagnetic coupler with strong resistance to radial offset as the research object and establishes a three-dimensional electromagnetic thermal field simulation model that considers the influence of insulation induced turn spacing and temperature on material properties, achieving accurate calculation of hot spot temperature.

Firstly, theoretical analysis is conducted on the temperature rise of the electromagnetic coupler, which involves calculating the loss of the electromagnetic coupler, and then loading the loss as a heat source into the heat flow field to obtain the temperature distribution of the electromagnetic coupler; Secondly, based on the temperature rise principle of electromagnetic couplers and strictly following the temperature rise calculation process, the finite element Ansys Maxwell simulation software is used to solve the electromagnetic loss and hot spot distribution. In the simulation software, the first step is to calculate coil losses and magnetic core losses.Then, using the transient electromagnetic field heat flow field bidirectional indirect coupling method, a three-dimensional electromagnetic field simulation model is established. Finally, the electromagnetic loss calculation results are coupled as loads into the heat flux field and corresponding relationships are established to obtain the hot spot distribution of the rotating electromagnetic coupler.

The simulation results show that the magnetic field distribution of the electromagnetic coupler is basically symmetrical from top to bottom, and the magnetic field intensity gradually decreases along the center position of the magnetic core towards the end. The maximum magnetic density is located on the surface of the magnetic core close to the geometric center of the secondary coil. The average loss density distribution of the electromagnetic coupler has the same characteristics as the magnetic induction intensity distribution. At 20 ℃, the average loss errors of the primary and secondary coils within one cycle were obtained through theoretical analysis and finite element simulation, which were 2.9% and 3.4%. The experimental results show that the deviation between the parameters measured by the LCRand the numerical values in simulation is less than 10%. The thermocouple multi-channel thermometer records the temperature of the system when it reaches steady state at 20 ℃. The high-temperature area is located at the top of the secondary coil, and the deviation between the simulated temperature of the hot spot and the measured value is less than 6%, verifying the rationality of the modeling method.

The following conclusions can be drawn through experimental analysis:(1) The average loss density of the magnetic core of the rotating electromagnetic coupler is positively correlated with the distribution of magnetic induction intensity; (2) The relative position of the coil, flow velocity, and convective heat transfer coefficient all have an impact on the operating temperature of the system. At the same time, due to the influence of the nested structure, local overheating will occur under the high temperature conduction of the coil, resulting in the severe temperature rise area of the rotating electromagnetic coupling concentrated in the middle of the LCT; (3) In order to simulate the temperature rise distribution of rotary electromagnetic couplers more accurately, in addition to considering the turn spacing caused by insulation, the temperature effect of material properties should also be considered.

Wireless power transfer (WPT), rotary electromagnetic coupler, multi-physical field coupling, temperature rise characteristic, electromagnetic loss

國家自然科學(xué)基金資助項目（51877129）。

2023-02-17

2023-04-01

10.19595/j.cnki.1000-6753.tces.230205

TM724

程志遠(yuǎn) 男，1977年生，博士，高級工程師，研究方向為無線電能傳輸、高壓電纜無線取電。E-mail：chengzhiyuan@126.com

李東東 男，1976年生，博士，教授，研究方向為風(fēng)力發(fā)電與電力系統(tǒng)穩(wěn)定控制、智能用電等。0E-mail：powerldd@163.com（通信作者）

（編輯 郭麗軍）