异步电机铁心损耗计算方法的分析与研究

（国网浙江省电力有限公司岱山县供电公司，浙江 舟山 316200）

文献[1-2]介绍了几种分析异步电机铁心损耗的方法。文献[3-4]分析了输入电压中谐波分量引入的铁心损耗。文献[5]采用有限元数值分析法分析了由PWM变换器驱动的异步电机内定子铁心的损耗。文献[6]中假定异步电机的转子铁心损耗为零，进而分析了定子铁心内的损耗，但是负载运行状态下，异步电机的转子铁心损耗无法被忽略。文献[7-8]分析了定子铁心内一个点的磁感应强度，进而计算定子铁心损耗，然而异步电机内各组成部件的磁感应强度分布并不相同。文献[9-10]分析了异步电机铁心损耗的系数，然而该方法未考虑异步电机各组成部件内电磁场分布的不同。

异步电机定、转子铁心中磁感应强度分布受到定、转子铁心开槽与电源内的高次谐波的影响，进而影响铁心损耗。文献[11-14]介绍了异步电机定、转子铁心内磁感应强度非正弦分布时，铁心损耗的计算方法。文献[15]计算了定、转子铁心内磁感应强度的变化率，进而采用一种等效频率法计算铁心损耗。文献[16]采用了一种广义斯坦梅茨方程分析磁感应强度非正弦分布时异步电机的铁心损耗。文献[17]采用了一种局部回路法推导的改进广义斯坦梅茨方程计算异步电机的铁心损耗。

本文提出了一种计算不同工况下异步电机各部分铁心的铁心损耗的通用方法。该方法首先采用有限元数值分析软件Ansoft分析异步电机内非正弦分布的磁感应强度。通常，制造商提供的硅钢片铁心损耗数据是某些特定的频率、特定的磁感应强度下，单位质量或单位体积的铁心损耗。本文采用了一种曲线拟合方法，可获得任意频率下的铁心损耗曲线，进而根据有限元软件分析的磁感应强度计算异步电机的铁心损耗。

本文搭建了试验平台，分析了不同负载功率下、不同开关频率下异步电机的铁心损耗。本文一共采用了4种方法计算异步电机的铁心损耗，分别为本文提出的方法、二维有限元法、改进的广义斯坦梅茨方程法以及基于速率模型的方法。分析对比结果表明，本文提出的方法分析的铁心损耗具有很高的准确度。

1 铁心损耗计算方法

异步电机铁心损耗包含2个部分，涡流损耗与磁滞损耗。部分文献中，将附加损耗也认为是铁心损耗。涡流损耗与磁滞损耗均取决于硅钢片中磁感应强度的变化率。由于定、转子铁心开槽，定、转子铁心饱和与变频器供电电源内高次谐波的影响，定、转子铁心内磁感应强度为非正弦分布。

1.1 斯坦梅茨方程

斯坦梅茨提出了分析铁心损耗的一般斯坦梅茨方程，如下式（1）所示：

式中：P为单位质量的铁心损耗；Bmax为磁感应强度幅值；f为磁化频率；Cm，α，β为常系数，α的取值范围为[1，3]，β的取值范围为[2，3]。

此后，许多学者对一般斯坦梅茨方程进行了改进，提出了许多不同形式的斯坦梅茨方程。如将铁心损耗分为磁滞损耗与涡流损耗时，斯坦梅茨方程可表示为

式中：kh，ke分别为磁滞损耗系数与涡流损耗系数。

为了减小式（2）计算铁心损耗的误差与计算铁心的异常损耗，斯坦梅茨方程可表示为

式中：ka为异常损耗系数。

式（1）～式（3）所示的方程中，磁感应强度是正弦分布的。

为分析磁感应强度非正弦分布的铁心内的损耗，可采用下式所示的改进斯坦梅茨方程：

式中：feq为铁心内磁感应强度的等效频率；B为磁感应强度；T为磁化周期；ΔB为一个周期内磁感应强度的变化值。

式（4）也可改写为另一种形式，如下式所示：

式（4）与式（5）所示的改进斯坦梅茨方程也可应用于磁滞损耗与涡流损耗的计算中。此时，铁心损耗可表示为

考虑铁心的异常损耗时，改进斯坦梅茨方程可表示为

考虑铁心的异常损耗、变频器供电电源内的谐波含量与磁感应强度幅值的变化，铁心损耗可表示为

式中：γ，δ均为考虑铁心异常损耗的常系数。

1.2 曲线外推法

本文提出了一种外推法，该方法认为斯坦梅茨方程中的系数为磁化频率f与磁感应强度幅值Bmax的函数。若斯坦梅茨方程中的系数与磁化频率间的关系已知，则可计算任意频率下的铁心损耗。

本文采用该外推法分析了硅钢片50WW470的斯坦梅茨方程中的系数和磁化频率f、磁感应强度幅值Bmax之间的关系。式（2）所示的斯坦梅茨方程中包含4个未知系数，即可通过4组数据获得未知系数。表1为磁化频率为50 Hz，100 Hz，200 Hz以及400 Hz时、磁感应强度幅值为0.1 T，0.2 T，0.3 T以及0.4 T时，硅钢片50WW470的铁心损耗。

表1 硅钢片50WW470的铁心损耗Tab.1 The loss data of silicon steel sheet 50WW470

图1所示为硅钢片50WW470在不同频率下的铁心损耗曲线。

图1 硅钢片50WW470的损耗曲线Fig.1 The loss data of silicon steel sheet 50WW470

根据表1所示的数据，本文计算了磁化频率为50 Hz，100 Hz，200 Hz以及400 Hz时，式（2）所示的斯坦梅茨方程中的4个系数，如表2所示。

表2 硅钢片50WW470的铁心损耗系数Tab.2 The coefficient of core loss of silicon steel sheet 50WW470

本文采用式（12）所示的拟合公式，分析了不同频率下kh，ke，α，β，如下式所示：

式中：a，b，c为系数，且a＞0，b＜0，c＞0。

图2为系数kh随频率的变化曲线。由图2可知，kh为幂函数，且不同频率下的kh满足下式：

图2 不同磁化频率下kh的值Fig.2 The value of khunder different magmerionriom magnetization frequency

通过该方法可确定不同频率下kh的大致取值范围。其他系数ke，α，β的取值方法类似。

利用外推法获得系数后，可采用斯坦梅茨方程计算磁化频率为任意值时，硅钢片的铁心损耗。表3为磁化频率为200 Hz与1 000 Hz时，硅钢片50WW470的铁心损耗。

表3 硅钢片50WW470不同频率下的铁心损耗Tab.3 The core loss of silicon steel sheet 50WW470 at different magnetization frequencies

由表3可知，采用该外推法计算的铁心损耗与硅钢片生产制造公司提供的实验数据的误差较小，均在10%以下。

1.3 异步电机定转子铁心

异步电机铁心如图3所示。为简化分析一台异步电机的铁心损耗，将异步电机分为如图3a所示的定子齿距区域与转子齿距区域。定子齿距区域可分为3个组成部分，分别为：定子轭部区域、定子齿区域、定子槽口区域；转子齿距区域可分为3个组成部分，分别为：转子轭部区域、转子齿区域、转子槽口区域，如图3b以及图3c所示。图3b以及图3c中，定子齿距区域与转子齿距区域划分为许多矩形单元，且各矩形单元足够小到可认为该矩形单元内的磁感应强度为常数。因此矩形单元越小，分析结果越精确。然而矩形单元越小，矩形单元数越多，计算量越大，而且矩形单元数量达到一定程度后，继续增加矩形单元的数量，对计算结果的精确度影响较小。因此需权衡矩形单元数量与计算量的大小。有限元分析法将异步电机铁心划分为三角形单元，三角形单元相比于矩形单元更为优越，然而采用三角形单元，计算量大。

图3 异步电机铁心Fig.3 The core of induction motor

异步电机铁心分块如图4所示。

图4 异步电机铁心分块Fig.4 The segment of iron core in induction motor

本文采用了一台1.1 kW，50 Hz的异步电机有限元模型，阐述该损耗计算方法的具体流程，如图4a所示。由于定、转子铁心内的电磁场频率不一样，转子铁心内的电磁场频率较低，因此需设置2种求解步长。并将定子齿距区域剖分为28个单元，转子齿距区域剖分为36个单元。定子齿距区域的剖分方法如图4b所示，转子齿距区域的剖分方法如图4c所示，每个剖分单元内的各点的电磁场为常数。

在某个时刻，图4b中单元13内的线A与线B的磁感应强度分布与单元3内的线C与线D的磁感应强度分布分别如图5a、图5b、图5c以及图5d所示。

图5 某一时刻沿着线段的磁感应强度分布Fig.5 The magnetic induction intensity distribution versus the line segments at a certain time

由图5可知，磁感应强度在单元3与单元13内的变化很小。因此，在稳态运行时，可通过提取各单元内某点的磁感应强度计算铁心损耗。

2 铁心损耗计算过程

2.1 feq与 Bmax分析

采用本文提出的方法计算异步电机内的铁心损耗时，第1步为计算feq与Bmax。

利用二维有限元软件Ansoft可分析异步电机内任意一点的随时间变化的磁感应强度波形，并可通过下式获得feq：

式中：B1，…，Bk，…，Bn分别为时间点 t1，…，tk，…tn时的磁感应强度；Bmin为磁感应强度最小值。

图6为转速为2 880 r/min时，异步电机定子铁心单元5内某一点的磁感应强度波形。由式（15）可知，feq为296 Hz，Bmax为0.724 T。

图6 1个周期内某一点的磁感应强度随时间的分布Fig.6 The magnetic induction intensity versus the time within one period at a certain point

2.2 硅钢片50WW470的损耗曲线

采用本文提出的方法计算异步电机铁心损耗时，第2步为分析等效磁化频率为296 Hz时硅钢片的损耗曲线。

由前文分析的外推法可知，50WW470在频率为296 Hz时的损耗曲线如图7所示。

图7 feq为296 Hz时，50WW470的损耗曲线Fig.7 The loss curve of 50WW470 at feq=296 Hz

2.3 计算铁心损耗

采用本文提出的方法计算异步电机内的铁心损耗时，第3步为计算各单元的铁心损耗。

由有限元仿真软件计算可知单元5的面积为4.8 mm2，而该电机的轴向长度为0.135 m。因此单元5的体积为651 mm3，质量为0.005 013 kg。利用外推法与式（1）～式（4）、式（5）～式（12），可分析获得单元5的铁心损耗，如表4所示。由表4可知，式（1）～式（4）、式（5）～式（12）所示的斯坦梅茨方程计算的铁心损耗的结果非常接近。

表4 单元5的铁心损耗Tab.4 The core loss of segment 5

3 实验方法分析

为了验证所提出的计算铁心损耗的方法，设计了2个异步电机实验平台，实验平台示意图如图8所示。其中一台异步电机有正弦电流供电，如图8a所示，另外一台异步电机由PWM变频器供电，如图8b所示，实验平台实物图如图8c所示。

图8 实验平台示意图Fig.8 The schematic diagram of experimental platform

对于由正弦电流供电的异步电机而言，计算了供电电压不变情况下，不同负载时异步电机的损耗。对于由PWM变频器供电的异步电机而言，计算了不同开关频率下异步电机的损耗。下式所示为利用实验平台计算铁心损耗的方法。

式中：Pcore为异步电机内的铁心损耗；Pin为输入功率；Pout为输出功率；Pmech为包含风磨损耗与摩擦损耗的机械损耗；Pcopper为总铜耗。

3.1 输入与输出功率

输入功率Pin可通过三相功率分析仪获得，Pout可通过下式获得。

式中：ωm为转子速度；TL为负载转矩。

异步电机的转速可通过数字式转速表测得，负载转矩由涡流制动器提供。

3.2 机械损耗

机械损耗与转速密切相关，速度不变时异步电机的机械损耗不变。因此可通过直流电机测量在固定转速下的异步电机的机械损耗。首先测量直流电机在固定转速下的空载输入功率Pdc-without-IM，然后由直流电机驱动异步电机旋转在固定转速下，测量此时直流电机的输入功率Pdc-with-IM。则异步电机的机械损耗可表示为

3.3 铜损耗

异步电机工作时，定子绕组与转子鼠笼条中均会产生铜损耗。利用功率分析仪可获得定子电流，利用万用表可测量定子绕组的电阻，则定子绕组的铜损耗可表示为

式中：Pcu-stator为定子绕组铜损耗；Is为定子电流有效值；Rs为定子绕组电阻；m为相数。

由异步电机的T型等效电路可知，转子鼠笼条的铜损耗可表示为

式中：Pcu-rotor为转子鼠笼条的铜损耗；s为转差率；ωs为同步转速。

4 实验结果分析

4.1 由正弦电流供电的异步电机铁心损耗分析

图9为正弦电流供电时，本文提出方法以及其他方法计算的异步电机铁心损耗实验结果图。由前文分析可知，式（1）～式（4）、式（5）～式（12）分析计算铁心损耗的差距很小，可忽略不计。因此，本文选择式（1）所示的斯坦梅茨方程计算了铁心损耗。其他计算损耗的方法选用了IGSE法、二维有限元软件分析法以及基于速率模型[18]的算法。

图9 异步电机铁心损耗Fig.9 The core loss of induction motor

采用不同方法计算的铁心损耗与实验相比的误差率E，可通过下式计算获得：

式中：Pcore-expt为采用实验计算的铁心损耗；Pcore-method为采用不同方法计算的铁心损耗。

表5为由正弦电流供电时，异步电机铁心损耗的误差率。

表5 异步电机不同转速下铁心损耗误差率Tab.5 The error rate of core loss of induction motor at different speed

由图8与表5所示的结果可知，与其他方法相比，采用本文所提出方法计算的铁心损耗与实验分析获得铁心损耗较为接近。而且异步电机低速运行或重载运行时，铁心损耗增加。

4.2 由PWM变换器供电的异步电机损耗分析

图10a为开关频率为300 Hz与5 000 Hz时，由PWM变换器供电的异步电机定子绕组电流，图10b所示为定子电流的谐波含量。由图10b可知，开关频率低时定子电流的谐波含量高。

图10 异步电机定子电流Fig.10 The stator current of induction motor

图11为采用本文提出方法与实验分析获得的随开关频率变化的异步电机铁心损耗以及损耗的误差比。由图11可知，不同的开关频率下，两种方法获得铁心损耗较为接近，即该方法可准确分析异步电机的铁心损耗。

图11 不同开关频率下的异步电机的铁心损耗Fig.11 The core loss of induction motor at different switch frequencies

5 结论

本文提出了一种计算各种工况下异步电机各部分铁心损耗的方法。该方法可分析任意频率下硅钢片的损耗曲线、斯坦梅茨方程系数。该方法需通过有限元仿真软件获得异步电机内各部分铁心的磁感应强度分布，确定异步电机内各部分铁心的等效磁化频率。然后根据该磁化频率下的损耗曲线、斯坦梅茨方程系数分析异步电机各部分的铁心损耗，最终获得异步电机的总铁心损耗。该方法考虑了异步电机内谐波磁场对铁心损耗的影响。

本文采用提出的方法分析了正弦激励的异步电机在不同负载工况下的铁心损耗与PWM变换器激励的异步电机在不同开关频率下的铁心损耗。与实验分析结果的对比表明该方法可准确计算异步电机的铁心损耗。