無線充電系統旋轉式電磁耦合器損耗計算及熱點溫度研究

程志遠 宋曉逸 吳曉婷 貴子航 李東東

無線充電系統旋轉式電磁耦合器損耗計算及熱點溫度研究

程志遠 宋曉逸 吳曉婷 貴子航 李東東

（上海電力大學電氣工程學院 上海 200090）

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

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

0 引言

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

圖1 磁耦合諧振式WTP系統示意圖

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

圖2 AUV對接系統示意圖

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

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

1 電磁場計算模型

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

圖3 旋轉式電磁耦合器示意圖

表1 實驗樣機具體參數

Tab.1 Specific parameters of experimental prototype

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

2 多物理場耦合計算原理

圖4為求解旋轉式電磁耦合器三維磁-熱-流耦合計算流程。首先對三維電磁模型設置與溫度相關的材料屬性，在初始溫度下計算電磁耦合器的損耗密度；其次采用間接耦合方式將電磁損耗作為熱流場激勵，利用網格-節點單元對應關系對流熱場模型的流體、傳熱及能量方程進行耦合分析，計算出三維模型的表面傳熱系數；然后，根據溫度結果反饋至電磁場中得到對應的材料系數后，再反復迭代計算；最后，當迭代步數內監測點的相鄰兩步間溫度差值小于1 K時，終止迭代并得到電磁耦合器穩態溫升分布。

圖4 多場耦合計算流程

2.1 瞬態電磁場計算

2.2 電磁損耗計算

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

2.2.1 線圈損耗計算

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

假設除了該待求線圈外，系統中還存在個其他線圈，同時認為Litz線中每股細線導通電流及物理參數均相等，忽略線圈彎曲部分的內外半徑延展導致的輕微變化，根據文獻[20]式（3），skin、in_pro、turn_pro和coil_pro分別為

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

2.2.2 磁心損耗計算

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

2.3 熱流場計算方程

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

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

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

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

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

其中，面等效發射率r為

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

3 計算結果與分析

3.1 電磁損耗計算分析

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

圖5 電磁耦合器磁場分布

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

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

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

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

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

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

圖7 熱流場模型

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

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

圖8 線圈等效模型

3.3 穩態熱流場計算結果分析

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

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

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

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

圖10 副邊磁心溫升分布

3.4 溫升測試結果

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

圖11 WPT系統實驗平臺

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

表2 旋轉式電磁耦合器仿真與測量參數比較

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

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

圖12 功率分析儀采集數據

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

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

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

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

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

圖14 紅外線穩態溫升分布結果

表3 熱點溫度仿真與實驗結果

Tab.3 Hot spot temperature simulation and experimental results

4 結論

鑒于電磁耦合器的工作溫度是決定WPT系統穩定運行的關鍵因素，本文建立了考慮絕緣引起匝間距的旋轉式電磁耦合器三維電磁場仿真模型，并對其磁場與損耗進行計算和分析。采用磁-熱-流多物理場耦合方法計算自然狀態下表面傳熱系數，實現對電磁耦合器的溫升計算，最終通過所搭建的實驗平臺驗證理論分析的正確性，可得到以下結論：

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

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

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

本文的分析計算可為無線充電系統旋轉式電磁耦合器的損耗計算及溫升分布研究提供理論依據。后續將開展對旋轉式電磁耦合器應用于水下大功率下的溫升試驗，以及針對高溫部分展開散熱優化研究。基于本文模型，可以對AUV的傳輸特性隨溫度的變化趨勢展開深入分析，本文的多物理場耦合分析也為AUV的結構設計提供思路。

[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] 閆爭超, 胡謙宇, 趙晨旭, 等. 水下航行器感應式無線電能傳輸技術研究綜述[J/OL]. 中國電機工程學報, 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] 陳慶彬, 范峰, 汪金帥, 等. 無線電能傳輸磁耦合系統Litz線圈交流電阻精確評估方法[J]. 電工技術學報, 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] 張獻, 韓大穩, 沙琳, 等. 一種共享磁通多耦合模式的無線電能傳輸系統抗偏移方法[J]. 電工技術學報, 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] 毛行奎, 陳為. 開關電源高頻功率平面變壓器熱設計研究[J]. 電工電能新技術, 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] 程志遠, 陳坤, 李東東, 等. 旋轉式無線充電系統偏移特性研究[J]. 電工技術學報, 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] 李永建, 閆鑫笑, 張長庚, 等. 基于磁-熱-流耦合模型的變壓器損耗計算和熱點預測[J]. 電工技術學報, 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

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

2023-02-17

2023-04-01

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

TM724

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

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

（編輯 郭麗軍）