分别为位置变化前与位置变化后的次级线圈。Tx与Rx同轴且不存在偏转，属于对称分布。图1b～图1f中Tx与均属于非对称分布，其中包括水平偏移、水平角度偏转、垂直角度偏转、水平角度偏转加偏移和垂直角度偏转加偏移五种情况。图1中线圈位置变化参数的符号及具体含义见表1。摘要 互感是无线电能传输系统的一个关键参数，在电动汽车线圈耦合机构中添加磁屏蔽材料，可以有效减少磁泄漏，增强线圈之间的互感。目前还没有方法计算带磁屏蔽任意位置矩形线圈之间的互感。为此该文基于双傅里叶变换和麦克斯韦方程组推导了带磁屏蔽条件下的磁通密度。并在此基础上，利用空间直角坐标变换法，得到了任意相对位置两个矩形线圈的互感计算式。与传统的近似计算方法不同，所提出的方法是计算互感的解析方法，可以对线圈间互感进行精确的数值求解。最后，以两个具有多种位置变化情况的矩形螺旋线圈为例进行了模型验证，互感的理论计算值和有限元仿真结果及实验结果吻合良好，验证了所提方法的有效性。
关键词：无线电能传输 矩形线圈 互感计算 磁屏蔽
无线电能传输（Wireless Power Transfer, WPT）系统由于其便利性、可靠性和安全性等优点[1-3]，被广泛应用于电动汽车[4-6]、轨道交通[7]、消费电子[8]、医疗器械[9]、水下供电[10-11]和自动化制造装备等领域[12-13]。然而，在实际应用过程中，初级线圈和次级线圈之间不可避免地会现位置偏移情况，导致线圈之间的互感发生变化[14]，从而导致效率降低和输出电压波动[15]。因此，研究互感随线圈位置变化的规律具有重要意义。
在WPT系统中，常用的传输线圈是圆形线圈或矩形线圈。其中关于两个圆形线圈间的互感计算方法已经有大量文献研究，这些工作分别基于近似公式[16-17]、兰布达函数[18]、贝塞尔和斯特鲁夫函数[19-20]以及Biot-Savart定律[21-22]。然而，从实际的角度来看，具有两个矩形线圈的WPT系统在某些情况下更为必要和实用，因为在发生偏移时，矩形线圈比圆形线圈的互感波动更小。尽管这样，关于矩形线圈的互感计算方法的研究仍然较少[23]。文献[24]通过诺依曼公式，计算了同轴的两个矩形线圈之间的互感。文献[25]则利用Biot-Savart定律，提供了一种考虑线圈水平偏移和角度偏转的互感计算方法。在此基础上，文献[26-27]针对矩形线圈常见的相对位置变化情况，给出了对应的互感表达式。文献[28]将导线抽象为截面无穷小的理想线段，以任意位置两根通电导线的互感计算式为基础，解决了任意位置两个矩形线圈的互感计算问题。文献[29]从聂以曼公式入手，推导出了任意位置多边形线圈互感的一般表达式。尽管这些研究已经分析了线圈间相对位置各种变化情况，但是并没有分析磁屏蔽层对线圈互感的影响，从已有研究可知，在WPT线圈结构中添加磁屏蔽材料后，磁通路径发生改变，使得非工作区的磁场强度减小，同时线圈间耦合得到增强[30]。文献[31]基于二阶矢量位，推导了带磁屏蔽回折线圈阻抗的通用解析表达式。文献[32]则通过双傅里叶变换建立了铁磁板上矩形线圈的三维磁场解析模型。在此基础上，文献[33]从傅里叶分析的扩散方程出发，计算平行导体板上四个初级线圈与一个次级线圈在不同位置时的互感。文献[34]利用傅里叶-贝塞尔变换和双傅里叶变换，建立了双边带磁屏蔽的两个平行矩形螺旋线圈的互感解析模型，但是没有考虑线圈间的相对位置变化情况。综上所述，对于带磁屏蔽任意位置矩形线圈的互感计算方法还有待进一步研究。
本文建立了带磁介质矩形平面螺旋线圈在任意相对位置的互感解析模型。首先，利用傅里叶积分变量表达矩形线圈之间的磁场分布；其次，利用空间直角坐标变换法和磁通密度法，求解磁介质上任意相对位置的两个单匝矩形线圈之间的互感；然后，将矩形螺线线圈近似为一组同心的单匝矩形线圈，利用叠加法推导了矩形螺线线圈的互感计算公式；最后，通过仿真与实验验证了所推导互感计算式的正确性。
初级线圈与次级线圈之间的互感与相对位置有关：当两线圈对称分布且传输距离较近时，互感较大；当两线圈垂直分布且传输距离较远时，互感较小。事实上，两个线圈的位置常常不是固定不变的。本文总结了矩形线圈常见的几种位置变化情况，并针对每种位置变化情况，通过仿真与实验对所得互感计算公式进行了验证。
矩形线圈常见的几种位置变化如图1所示。图中Tx为初级线圈，Rx与分别为位置变化前与位置变化后的次级线圈。Tx与Rx同轴且不存在偏转，属于对称分布。图1b～图1f中Tx与均属于非对称分布，其中包括水平偏移、水平角度偏转、垂直角度偏转、水平角度偏转加偏移和垂直角度偏转加偏移五种情况。图1中线圈位置变化参数的符号及具体含义见表1。
图1 矩形线圈常见位置变化示意图
Fig.1 Schematic diagram of general position changes of rectangular coils
表1 线圈位置变化参数
Tab.1 The parameters ofcoils position variation
符号含义 ∆xX方向的偏移量 ∆yY方向的偏移量 ∆z垂直偏移量 θ水平偏转角度 γ垂直偏转角度
带磁屏蔽的两个矩形线圈示意图如图2所示，为了方便分析，将整个空间分为4个区间，其中区间2为厚度有限的磁介质，代表线圈的屏蔽层，其余为空气介质。区间3的上表面与XY坐标平面重合，s1和s2为初级线圈和次级线圈几何中心的高度；t、μr和σ分别表示屏蔽层的厚度、相对磁导率和电导率；a1、a2为初级线圈的半长与半宽，b1、b2为次级线圈的半长与半宽。与空芯矩形线圈模型不同，该模型里既存在Bi（入射磁通密度）又存在Br（反射磁通密度）。其中Bi由流过初级线圈的电流产生，Br则由磁屏蔽层中的感应涡流产生。通常认为区间3和区间4中只存在Br，而区间1和区间2中既含有Br也含有Bi[34]。
图2 带磁屏蔽矩形线圈的结构
Fig.2 Rectangular coils with magnetic shielding
对于空气中电流密度为J的载流导体，其在任意点P（x, y, z）处产生的磁矢位为[35]
（1）
式中，v为导体的电流分布；R为任意点P（x, y, z）到点源（
）之间的距离，可表示为
（2）
为求解上述方程，使用双傅里叶变换公式如下：
（3）
（4）
式中，
、
为双傅里叶积分变量。
将式（1）代入式（3）可得
（5）
其中
（6）
根据式（1）中的磁矢位表达式，可以求出入射磁通密度Bi的表示为[36]
式（7）经傅里叶变换可得
（8）
式中，符号“^”表示双傅里叶变换运算。
把初级线圈等效为四段导线l1、l2、l3和l4，导线l1、l3与X轴平行，导线l2、l4与Y轴平行。因此，磁矢位在X轴方向上的分量ax仅由导线l1、l3产生，而与导线l2、l4无关。根据式（5）得到ax的表达式为
式中，ax1、ax3分别为导线l1、l3产生的磁矢位分量。
同理得到磁矢位在Y轴方向上的分量为
（10）
将式（9）、式（10）代入式（8）得到区间3和区间4的入射磁通密度可表示为式（11）～式（16）。
区间3（0≤z＜s1）有
（11）
（12）
（13）
区间4（z≥s1）：
（14）
（15）
（16）
2.2.1 区间2
区间2为线圈的磁屏蔽层，根据麦克斯韦方程组可以得到关系式
（17）
（18）
式中，
为拉普拉斯算子；w为电源角频率；μ0为真空磁导率；
为相对磁导率。
将式（17）、式（18）代入式（3）可得
（19）
（20）
根据式（19）可设区间2的反射磁通密度的垂直分量为
（21）
式中，
；C2rz和C'2rz为待定系数。
已知
，由区间中电流密度满足J2rz=0可得b2rx和b2ry的关系式为
（22）
联立式（21）、式（22）可以得到b2rx、b2ry与b2rz的关系式为
（23）
（24）
2.2.2 区间1、区间3和区间4
区间1、区间3和区间4为空气介质所在区域，根据麦克斯韦方程组可以得到关系式
（25）
（26）
式中，n=1,3,4。
同理，将式（25）和式（26）代入式（3）可得
（27）
（28）
根据线圈无穷远处的磁通密度为零、式（27）可以得到区间1、区间3和区间4中反射磁通密度的垂直分量为
（29）
（30）
（31）
式中，C1rz、C3rz和C4rz为待定系数，且C3rz=C4rz。
已知该区域内不存在涡流，由
可得bnrx和bnry的关系式为
（32）
联立式（28）和式（32），得到bnrx、bnry和bnrz之间的关系为
（33）
（34）
根据电磁场边界条件：若分界面处不存在自由面电流，则磁通密度的法向分量和磁场强度的切向分量在分界面的两侧保持连续，即Bz和Hx在z=0和z=-t2平面上连续，其中Hx=Bx/m，由此得到
（35）
联立式（11）～式（13）、式（21）、式（23）、式（29）、式（30）、式（33）和式（35）可得
（36）
（37）
（38）
式中，
；
。
将式（14）～式（16）和式（31）～式（34）代入式（5）进行傅里叶逆变换，得到区域4中的磁通密度为
（39）
（40）
（41）
空间内同一点在不同直角坐标系中的坐标之间的关系，称为空间直角坐标变换。对式（39）～式（41）运用空间直角坐标变换可以得到任意直角坐标系下的磁通密度，进而得到任意相对位置的两个矩形线圈之间的互感。将空间直角坐标变换分为移轴变换、转轴变换和一般坐标变换。
（1）移轴变换：若原坐标系沿x轴、y轴和z轴分别平移Dx、Dy和Dz后得到新坐标系，则P点的原坐标（x, y, z）与其新坐标（x', y', z'）满足关系式
（42）
（2）转轴变换：若新坐标系是通过原坐标系绕x轴逆时针旋转角度y后得到，则P点的原坐标（x, y, z）与其新坐标（x', y', z'）满足
（43）
同理可得
（44）
（45）
式中，j为坐标系绕y轴旋转的角度；θ为坐标系绕z轴旋转的角度。
（3）一般坐标变换：指既包含移轴变换又包含转轴变换的坐标变换。若原空间直角坐标系在移轴变换的基础上，先后绕x轴、y轴和z轴旋转角y、角j和角θ后得到新坐标系，则
（46）
通过磁通密度计算线圈之间互感的公式[37]为
（47）
式中，s为次级线圈所围的区域；I为初级线圈的电流。
联立式（39）～式（41）、式（46）和式（47）可得任意相对位置带磁心的单匝矩形线圈的互感计算式为
（48）
式中，C4ix、C4iy、C4iz、C4rx、C4ry和C4rz可由式（14）～ 式（16）和式（36）～式（38）得到。
多匝矩形螺旋线圈如图3所示，矩形螺旋线圈可以近似为一组同心的单匝矩形线圈。根据叠加原理，两组同心的单匝矩形线圈之间的互感可以表示为
（49）
式中，Np、Ns分别为初级线圈和次级线圈的匝数；Mmn为第m匝初级线圈和第n匝次级线圈间的互感。
图3 多匝矩形螺旋线圈
Fig.3 Multi-turn rectangular coils
仿真在Ansys Maxwell 15.0有限元工具上进行，仿真模型如图4所示。实验采用文献[19]中使用的测量方法，在非磁性的木桌和亚克力框架上进行，实验装置如图5所示，初级线圈端设有尺寸66cm×66cm、厚1.5cm的铁氧体块作为磁屏蔽板，阻抗分析仪用于测量电感，其电流频率为 85 kHz。
图4 仿真模型
Fig.4 Simulation model
图5 实验装置
Fig.5 Experimental device
互感测量实验原理如图6所示，图中LT和LR分别表示Tx和Rx的自感。图6a中，线圈同相连接，线路两端的电感为L1=LT+LR+2M，图6b中，线圈反相连接，线路两端的电感为L2=LT+LR-2M。因此，Tx和Rx之间的互感M=|L1-L2|/4。
图6 互感测量实验原理图
Fig.6 Schematic diagram of mutual inductance measurement experiment
仿真和实验过程中，初级线圈与磁屏蔽层保持不动，线圈间相对位置的改变由调整次级线圈来实现。图7为线圈位置变化的实验模型。线圈和磁屏蔽的参数见表2。
图7 线圈位置变化实验模型图
Fig.7 Experimental model of coils position variation
表2 线圈与磁屏蔽参数
Tab.2 Parameters of coil and magnetic shields
参数数值 线圈的半长ɑ1, b1/mm50 线圈的半宽ɑ2, b2/mm40 初级线圈的匝数Np12 次级线圈的匝数Ns10 初级线圈的高度s1/mm17 次级线圈的高度s2/mm80 铜线直径/mm2 相邻线圈半长与半宽的改变量/mm2.4 磁屏蔽厚度t/mm15
图1a为线圈发生垂直偏移的示意图。当线圈以步长10mm发生垂直偏移，线圈的传输距离D（次级线圈与初级线圈几何中心的垂直距离）从100mm增加到180mm时，线圈的互感见表3。表3中，ε1表示仿真误差，ε2表示实验误差，两者表达式为
（50）
（51）
式中，Mc为本文所得互感计算公式得到的互感计算值；Ms为互感仿真值；Me为实验测量值。
表3 垂直偏移时的互感值及误差
Tab.3 Mutual inductance for vertical misalignment
D/mmMc/μHMs/μHMe/μHε1(%)ε2(%) 1002.694 6 2.703 52.737 8 0.331.58 1102.263 6 2.270 62.300 5 0.311.60 1201.912 1 1.920 11.942 2 0.421.55 1301.628 5 1.634 41.657 5 0.361.75 1401.389 6 1.399 21.420 5 0.692.18 1501.198 5 1.205 91.223 7 0.612.06 1601.033 7 1.042 11.053 7 0.811.90 1700.899 8 0.908 60.918 8 0.972.07 1800.786 4 0.794 40.803 7 1.012.15
根据表3，可以看到，传输距离D从100mm变化到180mm的过程中，仿真误差ε1均小于2%，实验误差ε2均小于3%，计算结果与仿真及实验结果基本保持一致。图8所示为表3对应的互感变化曲线。图8中，互感随着传输距离的增大而逐渐减小，这是因为线圈间的耦合磁链随着传输距离的增大而减小，从而导致了线圈间的互感变小。过程中，仿真误差
图8 垂直偏移时互感变化曲线
Fig.8 Diagram of the mutual inductance for vertical misalignment
图1b为发生了水平偏移的矩形线圈。当线圈的传输距离为80mm，次级线圈沿X轴以步长10mm发生水平偏移，偏移量∆x从0mm增大到70 mm时，线圈的互感见表4。
表4 水平偏移时的互感值及误差
Tab.4 Mutual inductance for horizontal misalignment
Dx/mmMc/μHMs/μHMe/μHε1(%)ε2(%) 03.910 8 3.924 13.986 1 0.341.89 103.859 6 3.8733.934 9 0.351.91 203.709 7 3.725 23.786 1 0.422.02 303.469 2 3.489 73.544 4 0.592.12 403.164 6 3.184 43.229 3 0.622.00 502.808 6 2.830 22.881 9 0.762.54 602.425 8 2.444 62.489 8 0.772.57 702.024 62.049 72.077 0 1.232.52
根据表4，可以看到，水平偏移量∆x从0mm变化到70mm的过程中，仿真误差ε1均小于2%，实验误差ε2均小于3%，计算结果与仿真及实验结果基本保持一致。图9所示为水平偏移时对应的互感变化曲线。
图9 水平偏移时互感变化曲线
Fig.9 Diagram of the mutual inductance for horizontal misalignment
图9中，互感随着水平偏移量∆x的增大而迅速减小，这是因为线圈间重合的区域随着∆x增大而减小，从而导致线圈之间的耦合减少。
图1c为发生了水平角度偏转的矩形线圈。当传输距离为80mm，水平偏转角度θ从0°增大到90°时，线圈的互感见表5。根据表5，可以看到，水平偏转角度θ从0°变化到90°的过程中，仿真误差ε1均小于1%，实验误差ε2均小于2%，计算结果与仿真及实验结果吻合良好。图10所示为表5对应的互感变化曲线。图10中，θ从0°变化到90°而互感基本恒定，这是因为θ变化前后
几乎重合，耦合磁链变化很小，从而导致线圈间互感基本不变。
表5 水平角度偏转时的互感值及误差
Tab.5 Mutual inductance for horizontal deflection
θ/(°)Mc/μHMs/μHMe/μHε1(%)ε2(%) 03.910 8 3.924 13.948 50.340.96 153.900 6 3.918 8 3.946 3 0.461.16 303.887 2 3.903 6 3.929 5 0.421.08 453.872 8 3.889 3 3.919 6 0.421.19 603.865 9 3.882 3 3.909 9 0.421.13 753.861 7 3.881 0 3.900 4 0.500.99 903.860 4 3.880 8 3.890 1 0.530.76
图10 水平角度偏转时互感变化曲线
Fig.10 Diagram of the mutual inductance for horizontal deflection
图1d为线圈发生垂直角度偏转的示意图。实验分析了传输距离为80mm，偏转角度γ从0°增大到90°时互感的变化情况，互感计算值、仿真值与实验测量值的结果及误差分析见表6。
根据表6，可以看到，垂直偏转角度γ从0°变化到90°的过程中，仿真误差ε1均小于1%，实验误差ε2均小于4%，计算结果与仿真及实验结果保持一致。图11所示为水平角度偏转时对应的互感变化曲线。
表6 垂直角度偏转时的互感值及误差
Tab.6 Mutual inductance for vertical deflection
γ/(°)Mc/μHMs/μHMe/μHε1(%)ε2(%) 03.910 8 3.924 13.978 6 0.341.70 153.934 4 3.945 4 4.010 8 0.281.91 303.946 7 3.957 3 4.005 1 0.271.46 453.782 2 3.794 8 3.854 3 0.331.87 603.122 2 3.1330 3.201 60.352.48 751.760 3 1.767 6 1.815 60.413.04 900.000 0 0.001 9———
图11 垂直角度偏转时互感变化曲线
Fig.11 Diagram of the mutual inductance for vertical deflection
图11中，当γ＜45°时，随着γ的增大互感变化较小；而当γ＞45°时，磁链随着γ增大急剧减小，从而导致了互感的急剧减小。
图1e为线圈发生水平角度偏转加偏移时的示意图。传输距离为80 mm，水平偏转角θ分别为15°和45°的矩形线圈，在沿x轴方向发生偏移时的互感见表7。从表7可以看出，在偏移量∆x从-40mm变化到40 mm的过程中，仿真误差均小于1%，实验误差均小于2%，计算结果与仿真及实验结果吻合良好。图12所示为水平角度偏转加水平偏移情况对应的互感变化曲线。图12中，15°和45°时互感的变化曲线基本重合，说明水平角度偏转对线圈互感的影响不大，线圈互感主要受水平偏移影响，与图9、图10的结论一致。
图12 水平角度偏转加偏移时互感变化曲线
Fig.12 Diagram of the mutual inductance for horizontal deflection with horizontal misalignment
表7 水平角度偏转加偏移时的互感值及误差
Tab.7 Mutual inductance for horizontal deflection with horizontal misalignment
Dx/mm绕Z轴偏转15°绕Z轴偏转45° Mc/μHMs/μHMe/μHε1(%)ε2(%)Mc/μHMs/μHMe/μHε1(%)ε2(%) -403.169 23.185 53.207 80.511.203.171 53.194 03.207 80.701.13 -303.470 73.487 93.519 50.491.393.456 63.481 63.519 50.721.79 -203.701 33.719 83.748 30.501.253.680 13.700 13.748 30.541.82 -103.848 63.866 53.891 60.461.113.811 33.839 93.878 60.751.74 03.900 63.918 83.955 90.461.403.872 83.889 33.945 90.421.85 103.849 33.868 63.889 60.501.043.820 33.845 43.889 60.651.78 203.710 93.723 83.757 50.351.243.681 83.712 43.757 50.822.02 303.478 53.493 33.518 30.421.133.466 93.498 23.518 30.901.46 403.177 33.193 33.239 70.501.933.190 13.215 83.239 70.801.53
图1f为线圈垂直角度偏转加偏移的位置变化示意图。表8给出了传输距离为80mm，垂直偏转角γ分别为15°和45°的矩形线圈在发生偏移时的互感值。由表8可得，在偏移量∆x从-40mm变化到40mm的过程中，仿真误差均小于1%，实验误差均小于3%，计算结果与仿真及实验结果基本保持一致。图13所示为垂直角度偏转加偏移情况对应的互感变化曲线。图13中，15°曲线和45°曲线存在交叉，因为线圈在此处垂直偏转15°后的磁链与垂直偏转45°后的磁链相等，互感相等。当∆x≤-10mm时，45°的互感大于15°；当垂直偏移量在∆x＞-10mm时，15°的互感大于45°。15°曲线在∆x=0时取得最大值，而45°曲线在∆x=-20mm时取得最大值。
图13 垂直角度偏转加偏移时互感变化曲线
Fig.13 Diagram of the mutual inductance for vertical deflection with horizontal misalignment
表8 垂直角度偏转加偏移时的互感值及误差
Tab.8 Mutual inductance for vertical deflection with horizontal misalignment
Dx/mm绕Y轴偏转15°绕Y轴偏转45° Mc/μHMs/μHMe/μHε1(%)ε2(%)Mc/μHMs/μHMe/μHε1(%)ε2(%) -403.282 93.299 13.348 70.491.973.595 93.619 83.668 60.661.98 -303.570 33.584 63.651 20.402.223.857 13.879 53.938 80.582.07 -203.784 93.799 73.862 70.392.023.994 24.013 84.084 30.492.21 -103.911 63.924 23.998 70.322.183.972 13.988 34.069 30.412.39 03.934 43.945 94.023 20.292.213.782 23.794 83.870 90.332.29 103.839 93.858 63.932 00.492.343.442 63.458 23.525 80.452.36 203.644 53.665 33.720 90.572.053.000 33.011 03.072 90.362.36 303.367 93.382 63.440 70.442.122.480 52.488 92.539 30.342.32 403.011 03.024 43.074 60.442.071.909 91.918 51.950 20.452.07
本文推导了磁屏蔽环境中的任意位置矩形螺旋线圈之间互感的解析模型。该模型基于双傅里叶变换和麦克斯韦方程组，求解了次级线圈所在区间的入射磁通密度和反射磁通密度，然后利用空间直角坐标变换法，得到了任意位置带磁屏蔽矩形线圈的互感计算式。分析了两个矩形线圈在多种位置变化的情况，互感的计算值、仿真值和实验测量值具有良好的一致性，并且各种误差均小于5%，充分说明了本文所推导的互感计算公式的正确性。研究结果为带磁屏蔽的无线电能传输系统提供了线圈设计和参数优化的理论依据，对接下来双边带磁屏蔽层情况的研究也具有一定借鉴意义。
参考文献
[1] 贾金亮, 闫晓强. 磁耦合谐振式无线电能传输特性研究动态[J]. 电工技术学报, 2020, 35(20): 4217-4231.
Jia Jinliang, Yan Xiaoqiang. Research tends of magnetic coupling resonant wireless power transfer characteristics[J]. Transactions of China Electrotechnical Society, 2020, 35(20): 4217-4231.
[2] Choi S Y, Gu B W, Jeong S Y, et al. Advances in wireless power transfer systems for roadway-powered electric vehicles[J]. IEEE Journal of Emerging and Selected Topics in Power Electronics, 2015, 3(1): 18-36.
[3] Li Zhenjie, Zhu Chunbo, Jiang Jinhai, et al. A 3-kW wireless power transfer system for sightseeing car supercapacitor charge[J]. IEEE Transactions on Power Electronics, 2017, 32(5): 3301-3316.
[4] 周玮, 蓝嘉豪, 麦瑞坤, 等. 无线充电电动汽车V2G模式下光储直流微电网能量管理策略[J]. 电工技术学报, 2022, 37(1): 82-91.
Zhou Wei, Lan Jiahao, Mai Ruikun, et al. Research on power management strategy of DC microgrid with photovoltaic, energy storage and EV-wireless power transfer in V2G mode[J]. Transactions of China Electrotechnical Society, 2022, 37(1): 82-91.
[5] 赵靖英, 张振远, 张珂. 基于H∞非线性控制器的电动汽车无线充电系统的副边控制设计与参数优化[J]. 电工技术学报, 2022, 37(3): 566-577.
Zhao Jingying, Zhang Zhenyuan, Zhang Ke. Control design and parameter optimization on secondary side of electric vehicle wireless charging system based on H∞ nonlinear controller[J]. Transactions of China Electrotechnical Society, 2022, 37(3): 566-577.
[6] 吴丽君, 李冠西, 张朱浩伯, 等. 一种具有恒流恒压输出自切换特性的电动汽车无线电能传输系统拓扑[J]. 电工技术学报, 2020, 35(18): 3781-3790.
Wu Lijun, Li Guanxi, Zhang Zhuhaobo, et al. A wireless power transfer system topology with automatic switching characteristics of constant current and constant voltage output for electric vehicle charging[J]. Transactions of China Electrotechnical Society, 2020, 35(18): 3781-3790.
[7] 麦瑞坤, 李勇, 何正友, 等. 无线电能传输技术及其在轨道交通中研究进展[J]. 西南交通大学学报, 2016, 51(3): 446-461.
Mai Ruikun, Li Yong, He Zhengyou, et al. Wireless power transfer technology and its research progress in rail transportation[J]. Journal of Southwest Jiaotong University, 2016, 51(3): 446-461.
[8] Jang Y, Jovanovic M M. A contactless electrical energy transmission system for portable-telephone battery chargers[J]. INTELEC Twenty-Second International Telecommunications Energy Conference (Cat No 00CH37131), 2000: 726-732.
[9] Zargham M, Gulak P G. Maximum achievable efficiency in near-field coupled power-transfer systems[J]. IEEE Transactions on Biomedical Circuits and Systems, 2012, 6(3): 228-245.
[10] 吴旭升, 孙盼, 杨深钦, 等. 水下无线电能传输技术及应用研究综述[J]. 电工技术学报, 2019, 34(8): 1559-1568.
Wu Xusheng, Sun Pan, Yang Shenqin, et al. Review on underwater wireless power transfer technology and its application[J]. Transactions of China Electrotechnical Society, 2019, 34(8): 1559-1568.
[11] Yang Lei, Zhang Yuanqi, Li Xiaojie, et al. Analysis and design of four-plate capacitive wireless power transfer system for undersea applications[J]. CES Transactions on Electrical Machines and Systems, 2021, 5(3): 202-211.
[12] 卿晓东, 苏玉刚. 电场耦合无线电能传输技术综述[J]. 电工技术学报, 2021, 36(17): 3649-3663.
Qing Xiaodong, Su Yugang. An overview of electric-filed coupling wireless power transfer technology[J]. Transactions of China Electrotechnical Society, 2021, 36(17): 3649-3663.
[13] 薛明, 杨庆新, 章鹏程, 等. 无线电能传输技术应用研究现状与关键问题[J]. 电工技术学报, 2021, 36(8): 1547-1568.
Xue Ming, Yang Qingxin, Zhang Pengcheng, et al. Application status and key issues of wireless power transmission technology[J]. Transactions of China Electrotechnical Society, 2021, 36(8): 1547-1568.
[14] 李中启, 李上游, 李晶, 等. 动态无线电能传输系统多接收线圈正反串联结构的互感计算与优化[J]. 电工技术学报, 2021, 36(24): 5153-5164.
Li Zhongqi, Li Shangyou, Li Jing, et al. Mutual inductance calculation and optimization of multi-receiver positive and negative series coil structure in dynamic wireless power transfer systems[J]. Transactions of China Electrotechnical Society, 2021, 36(24): 5153-5164.
[15] Wang Junhua, Li Jiangui, Ho S L, et al. Lateral and angular misalignments analysis of a new PCB circular spiral resonant wireless charger[J]. IEEE Transactions on Magnetics, 2012, 48(11): 4522-4525.
[16] Raju S, Wu Rongxiang, Chan Mansun, et al. Modeling of mutual coupling between planar inductors in wireless power applications[J]. IEEE Transactions on Power Electronics, 2014, 29(1): 481-490.
[17] Zhang Dianhai, Koh C S. An efficient semianalytic computation method of magnetic field for a circular coil with rectangular cross section[J]. IEEE Transactions on Magnetics, 2012, 48(1): 62-68.
[18] Babic S, Sirois F, Akyel C, et al. New formulas for mutual inductance and axial magnetic force between a thin wall solenoid and a thick circular coil of rectangular cross-section[J]. IEEE Transactions on Magnetics, 2011, 47(8): 2034-2044.
[19] Acero J, Carretero C, Lope I, et al. Analysis of the mutual inductance of planar-lumped inductive power transfer systems[J]. IEEE Transactions on Industrial Electronics, 2013, 60(1): 410-420.
[20] Conway J T. Exact solutions for the mutual inductance of circular coils and elliptic coils[J]. IEEE Transactions on Magnetics, 2012, 48(1): 81-94.
[21] Ravaud R, Lemarquand G, Babic S, et al. Cylindrical magnets and coils: fields, forces, and inductances[J]. IEEE Transactions on Magnetics, 2010, 46(9): 3585-3590.
[22] Misakian M. Equations for the magnetic field produced by one or more rectangular loops of wire in the same plane[J]. Journal of Research of the National Institute of Standards and Technology, 2000, 105(4): 557-564.
[23] 吴德会, 何天府, 王晓红, 等. 感应电能传输中矩形螺线线圈互感耦合的解析建模与分析[J]. 电工技术学报, 2018, 33(3): 680-688.
Wu Dehui, He Tianfu, Wang Xiaohong, et al. Analytical modeling and analysis of mutual inductance coupling of rectangular spiral coils in inductive power transfer[J]. Transactions of China Electrotechnical Society, 2018, 33(3): 680-688.
[24] 邝向军, 贾连宝. 共轴矩形载流线圈间的互感系数和相互作用力计算[J]. 西华师范大学学报(自然科学版), 2017, 38(4): 426-429.
Kuang Xiangjun, Jia Lianbao. On the mutual induction coefficient and magnetic force between two coaxial rectangle coils carrying currents[J]. Journal of China West Normal University (Natural Sciences), 2017, 38(4): 426-429.
[25] Kim J, Son H C, Kim D H, et al. Efficiency of magnetic resonance WPT with two off-axis self-resonators[C]//2011 IEEE MTT-S International Microwave Workshop Series on Innovative Wireless Power Transmission: Technologies, Systems, and Applications, Kyoto, Japan, 2011: 127-130.
[26] Joy E R, Dalal A, Kumar P. Accurate computation of mutual inductance of two air core square coils with lateral and angular misalignments in a flat planar surface[J]. IEEE Transactions on Magnetics, 2014, 50(1): 1-9.
[27] 全倡辉. 无线电能传输系统中带磁屏蔽矩形线圈间互感的计算[D]. 株洲: 湖南工业大学, 2021.
[28] Cheng Yuhua, Shu Yaming. Mutual inductance calculation between arbitrarily positioned rectangular filaments[J]. International Journal of Applied Electromagnetics and Mechanics, 2014, 46(1): 287-298.
[29] 刘洋, 何圣民. 任意空间位置两多边形多匝线圈的互感计算[J]. 价值工程, 2013, 32(31): 235-237.
Liu Yang, He Shengmin. The calculation of mutual inductance of two polygons with multiturn coils at arbitrarily position[J]. Value Engineering, 2013, 32(31): 235-237.
[30] 许乔迪, 徐叶飞, 麦瑞坤. 具有目标面最优磁屏蔽效果的IPT谐振式无功屏蔽系统研究[J]. 中国电机工程学报, 2019, 39(18): 5490-5498, 5597.
Xu Qiaodi, Xu Yefei, Mai Ruikun. IPT resonant reactive shielding systems with the characteristics of optimal magnetic shielding effect on the target surface[J]. Proceedings of the CSEE, 2019, 39(18): 5490-5498, 5597.
[31] 郝宽胜, 黄松岭, 赵伟, 等. 基于二阶矢量位的矩形截面回折线圈阻抗和脉冲磁场的解析建模与计算[J]. 物理学报, 2011, 60(7): 784-793.
Hao Kuansheng, Huang Songling, Zhao Wei, et al. Analytical modelling and calculation of impedance and pulsed magnetic field for rectangular meander coil based on second order potential[J]. Acta Physica Sinica, 2011, 60(7): 784-793.
[32] Smeets J P C, Overboom T T, Jansen J W, et al. Mode-matching technique applied to three-dimensional magnetic field modeling[J]. IEEE Transactions on Magnetics, 2012, 48(11): 3383-3386.
[33] Smeets J P C, Overboom T T, Jansen J W, et al. Inductance calculation nearby conducting material[J]. IEEE Transactions on Magnetics, 2014, 50(11): 1-4.
[34] Luo Zhichao, Wei Xuezhe. Analysis of square and circular planar spiral coils in wireless power transfer system for electric vehicles[J]. IEEE Transactions on Industrial Electronics, 2018, 65(1): 331-341.
[35] 崔翔. 电流连续的细导体段模型的磁场及电感[J]. 物理学报, 2020, 69(3): 87-98.
Cui Xiang. Magnetic field and inductance of filament conductor segment model with current continuity[J]. Acta Physica Sinica, 2020, 69(3): 87-98.
[36] Wu Linlin, Lu Kaiyuan, Xia Yongming. Mutual inductance calculation of two coaxial solenoid coils with iron core[C]//2018 21st International Conference on Electrical Machines and Systems (ICEMS), Jeju, Korea (South), 2018: 1804-1808.
[37] Piri M, Jaros V, Frivaldsky M. Verification of a mutual inductance calculation between two helical coils[C]//2015 16th International Scientific Conference on Electric Power Engineering (EPE), Kouty nad Desnou, Czech Republic, 2015: 712-717.
Mutual Inductance Calculation of Arbitrarily Positioned Rectangular Coils with Magnetic Shielding in Wireless Power Transfer Systems
Abstract Mutual inductance is a key parameter of the wireless power transfer (WPT) system. Themagnetic leakage can be reducedand the mutual inductance between coupling coils can be enhancedby adding the magnetic shielding for the electric vehicle. However, there is currently no method to calculate the mutual inductance between arbitrarily positioned rectangular coils with magnetic shielding. In this paper, the magnetic flux density of coils with magnetic shielding is derived by the Dual Fourier transform and Maxwell's equations. On this basis, the mutual inductance of two arbitrarily positioned rectangular coils is obtained by applying the spatial rectangular coordinate transformation method. The obtained method is different from the traditional approximate calculation method. It is an analytical method that can accurately numerically solve the mutual inductance between coils. Finally, the mutual inductance calculation formula is verified by two rectangular spiral coils with various position changes. The calculated results are in good agreement with the simulated and experimental results, which verifies the effectiveness of the proposed method.
keywords：Wireless power transfer, rectangular coils, mutual inductance calculation, magnetic shielding
DOI:10.19595/j.cnki.1000-6753.tces.211293
中图分类号：TM724
国家自然科学基金（11901188）、湖南省教育厅项目（20B186, 18A272）和湖南省自然科学基金（2019JJ60066）资助。
收稿日期 2021-08-15
改稿日期 2021-12-16
李中启 男，1985年生，博士，助理教授，研究生导师，研究方向为无线电能传输技术。E-mail：my3eee@126.com
黄守道 男，1962年生，博士，博士生导师，研究方向为特种电机本体及控制、无线电能传输技术。E-mail：hsd1962@hnu.edu.cn（通信作者）
（编辑 郭丽军）