CN106326595B - 一种改进的开关磁阻电机铁耗计算模型建立方法 - Google Patents
一种改进的开关磁阻电机铁耗计算模型建立方法 Download PDFInfo
- Publication number
- CN106326595B CN106326595B CN201610804714.7A CN201610804714A CN106326595B CN 106326595 B CN106326595 B CN 106326595B CN 201610804714 A CN201610804714 A CN 201610804714A CN 106326595 B CN106326595 B CN 106326595B
- Authority
- CN
- China
- Prior art keywords
- magnetic induction
- loss
- frequency
- factor
- iron
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开一种用于开关磁阻电机铁耗计算的改进模型,属于开关磁阻电机铁耗计算及分析技术领域。该模型是通过对不同频率和磁场强度下铁磁材料实测的损耗值进行数据分析,得到磁滞耗系数、涡流损耗系数等几个关键参数随磁感应强度和频率变化的规律,从而得到正弦交变磁化时铁磁材料铁耗计算的变系数两项式计算模型。针对开关磁阻电机中磁感应强度波形非正弦的特殊情况,运用有限元计算得到电机每一部分的切向和径向磁感应强度随时间变化波形,分别对其进行傅里叶分解将非正弦磁感应强度波形分解成一列正弦波形,对每次谐波运用变系数两项式计算模型计算得到每次谐波的铁损,累加求和即可得到计及n次谐波单位质量铁损模型。
Description
技术领域
本发明涉及一种改进的开关磁阻电机铁损计算模型建立方法,尤其适用于各种相数及其结构的开关磁阻电机。
背景技术
开关磁阻电机具有结构简单,坚固耐用,运行可靠,能量转换效率高等优点,但是由于开关磁阻直线电机的双凸极结构和磁饱和特性使得开关磁阻电机磁感应强度的非正弦特性,造成电机的铁耗计算困难,弄清开关磁阻电机铁耗产生的机理能够为提高电机的工作效率找到理论依据和改进方向;在建立电机控制系统模型时可以考虑把铁耗包含在内,从而得到更加准确的系统控制模型;对于选择合适的散热方式也有很大帮助,延长电机的适用寿命。因此研究精确的适用于开关磁阻电机非正弦特性的损耗模型是亟待解决的难点问题。目前开关磁阻电机的损耗计算模型主要有两种,一是根据经典的Steinmetz公式建立的两项式计算模型,结合电机所用的铁磁材料铁损数据,根据Steinmetz公式拟合得到磁滞损耗系数、涡流损耗系数和Steinmetz因子,然后建立电机的铁耗计算模型;二是根据Bertotti提出的考虑异常损耗建立的三项式计算模型,拟合得到磁滞损耗系数、涡流损耗系数和异常损耗系数,最后建立电机的铁耗计算模型;上述模型都是在固定频率情况下所确定的系数,且上述所确定的系数为常数,拟合误差较大,对于开关磁阻电机磁感应强度非正弦变化、频率不固定的情况计算误差较大,并不太适用。
发明内容
针对上述技术中存在问题,提供一种在不同频率下计算精确,并能够适用于开关磁阻电机磁感应强度非正弦特性的铁损计算模型。
为实现上述技术目的,本发明改进的开关磁阻电机铁损计算模型是在通过对不同频率和磁场强度下铁磁材料实测的损耗值进行数据分析,并结合Steinmetz公式,得到Steinmetz公式中磁滞耗系数kh、涡流损耗系数ke和Steinmetz因子α这几个关键参数随磁感应强度和频率变化的规律,从而得到了正弦交变磁化时铁磁材料铁耗计算的变系数两项式计算模型。针对开关磁阻电机中磁感应强度波形非正弦的特殊情况,运用有限元计算得到电机每一部分的切向和径向磁感应强度随时间变化波形,分别对其进行傅里叶分解将非正弦磁感应强度波形分解成一列正弦波形,对每次谐波运用变系数两项式计算模型计算得到每次谐波的铁损,累加求和即可得到计及n次谐波单位质量铁损模型,建立电机的铁损模型。具体模型建立过程如下。
在磁场为正弦交变磁场,磁感应强度峰值接近饱和以及不计高频下集肤效应对涡流的影响时,Steinmetz铁耗计算模型表达式为:
Pfe=khfBα+kef2B2 (1)
式中,B为磁感应强度,f为频率,khBαf为磁滞损耗,keB2f2为涡流损耗,kh、ke和α分别为磁滞损耗系数、涡流损耗系数和Steinmetz因子,它们的取值直接影响到计算模型的准确性。
令khBα=C,keB2=D,并使左右两边同除以f则可得
通过双频法即可求得上式中的未知数C和D,因此可求得不同频率下的磁滞损耗系数、不同磁感应强度下的磁滞损耗系数、不同频率下的涡流损耗系数、不同磁感应强度下的涡流损耗系数、不同频率下的Steinmetz因子值和不同磁感应强度下的Steinmetz因子值,分析磁滞损耗系数和涡流损耗系数随频率和磁感应强度的变化趋势,得到磁滞损耗系数kh的取值受频率和磁感应强度变化的影响很小,可近似认为kh的取值为常数。
涡流损耗系数ke随频率的增大基本呈直线下降趋势,涡流损耗系数ke随磁感应强度的变化趋势不大,可忽略不计,因此涡流损耗系数可用以下表达式拟合:
Ke=af+b (3)
随着频率的变化Steinmetz因子值基本保持不变,可以认为Steinmetz因子的取值与频率无关,Steinmetz因子值随磁感应强度的增大而增大,其变化规律可用以下多项式拟合:
α=k1B3+k2B2+k3B+k4 (4)
通过以上分析可得到铁磁材料的两项式变系数计算模型如下:
以上所述两项式变系数计算模型只适用于磁感应强度波形为正弦波的情况下铁耗的计算,而在开关磁阻电机铁心中磁感应强度恰恰就不是正弦波形,因此不能直接用上述模型计算开关磁阻电机中的铁耗。因此假设磁感应强度波形的基波频率为f,经过傅立叶分解后,第m次谐波的幅值为Bm,频率为mf,则由第m次谐波产生的铁耗为
通过有限元计算分析得到的电机铁心某处的磁感应强度是以切向分量和径向分量的形式给出的,假设径向分量用BN表示,切向分量用BT表示,则第m次分量的可表示为:
则考虑切向和径向分量的第m次谐波产生的铁耗为
因此计及前n次谐波铁心中某一部分铁损计算模型为
对电机每一部分铁损运用上述计及前n次谐波铁损计算模型计算得到电机铁损,求和即可得到电机的铁损模型。
有益效果:本发明对各种相数的开关磁阻电机适用,也适用于各种相数的开关磁阻直线电机。通过对不同频率和磁场强度下铁磁材料实测的损耗值进行数据分析,得到正弦交变磁化时铁磁材料铁耗计算的变系数两项式计算模型。针对开关磁阻电机中磁感应强度波形非正弦的特殊情况,运用有限元计算得到电机每一部分的切向和径向磁感应强度随时间变化波形,分别对其进行傅里叶分解将非正弦磁感应强度波形分解成一列正弦波形,对每次谐波运用变系数两项式计算模型计算得到每次谐波的铁损,累加求和即可得到计及n次谐波单位质量铁损模型,建立电机铁损的模型,达到本发明的目的。该种改进的开关磁阻电机铁损计算模型能准确的对磁滞损耗和涡流损耗进行分离,精度度高,求解方法简单、快速,具有良好的工程应用价值。
附图说明
图1是本发明的改进前的50W470铁耗曲线拟合与实验对比图。
图2是本发明的磁滞损耗系数随频率变化曲线图。
图3是本发明的磁滞损耗系数随磁感应强度变化曲线图。
图4是本发明的涡流损耗系数随频率变化曲线图。
图5是本发明的涡流损耗系数随磁感应强度变化曲线图。
图6是本发明的Steinmetz因子随频率变化曲线图。
图7是本发明的Steinmetz因子随磁感应强度变化曲线图。
图8是本发明的改进后的50W470铁耗拟合与实验结果对比图。
图9是本发明的开关磁阻电机区域分块图。
具体实施方式
下面结合附图对本发明的一个实施例作进一步的描述。
在磁场为正弦交变磁场,磁感应强度峰值接近饱和以及不计高频下集肤效应对涡流的影响时,Steinmetz铁耗计算模型表达式为:
Pfe=khfBα+kef2B2 (1)
式中,B为磁感应强度,f为频率,khBαf为磁滞损耗,keB2f2为涡流损耗,kh、ke和α分别为磁滞损耗系数、涡流损耗系数和Steinmetz因子,它们的取值直接影响到计算模型的准确性。
令khBα=C,keB2=D,并使左右两边同除以f则可得
通过双频法即可求得上式中的未知数C和D,进而求得涡流损耗系数、磁滞损耗系数和Steinmetz因子的值,如图1为铁磁材料50W470通过以上所述的系数拟合的铁耗曲线和实验对比图,由实测曲线与拟合曲线对比可以发现在频率较低且磁感应强度较小时损耗的拟合曲线与实验曲线吻合的较好,这说明在这个范围内拟合准确度较好。当频率低于200Hz时在每个频率下随磁感应强度的增大拟合值与实验值的偏差越来越大,这说明涡流损耗系数、磁滞损耗系数和Steinmetz因子三个参数中至少有一个参数不是固定常数,会跟着磁感应强度的变化而变化;另外在相同的磁感应强度范围内,当频率升高时吻合度变得越来越差,这也说明涡流损耗系数、磁滞损耗系数和Steinmetz因子三个参数中至少有一个会随频率的变化而变化。
根据铁磁材料50W470不同频率下的损耗数据和上述双频法拟合得到这三个参数更多的值,磁滞损耗系数kh的取值受频率和磁感应强度变化规律如图2和图3所示;涡流损耗系数ke的取值受频率和磁感应强度变化规律如图4和图5所示;Steinmetz因子α的取值受频率和磁感应强度变化规律如图6和图7所示;由图2和图3的拟合结果可以看出,磁滞损耗系数kh的取值受频率和磁感应强度变化的影响很小,因此可以近似认为kh的取值为常数,其值大小约为0.0287。
由图4的拟合结果可以发现,涡流损耗系数ke随频率的增大基本呈直线下降趋势,频率每升高50Hz ke的取值大约减小0.000002;在图5中当磁感应强度小于1.5T时ke的值基本保持为2.7E-4不变,当磁感应强度大于1.5T时ke的值上升为3.6E-4,因此ke可以分段取值。如果考虑到多数情况下铁心的磁感应强度都在1.5T,也可以认为ke的值不随磁感应强度变化。因此的变化规律可以用以下表达式表示
由图6可知,随着频率的变化Steinmetz因子基本保持不变,可以认为Steinmetz因子的取值与频率无关。在图7中Steinmetz因子随磁感应强度的增大而增大,其变化规律可以用下面的多项式来表示
α=0.38B3-0.02B2-0.33B+1.44 (4)
通过以上分析可得到铁磁材料的两项式变系数计算模型如下
使用变系数两项模型拟合得到的50W470铁耗数据与实验数据对比图如图8所示,由图可知变系数两项式模型相对与常系数两项式模型在拟合精度上有了很大程度的提高。
以上所述两项式变系数计算模型只适用于磁感应强度波形为正弦波的情况下铁耗的计算,而在开关磁阻电机铁心中磁感应强度恰恰就不是正弦波形,因此不能直接用上述模型计算开关磁阻电机中的铁耗。因此假设磁感应强度波形的基波频率为f,经过傅立叶分解后,第m次谐波的幅值为Bm,频率为mf,则由第m次谐波产生的铁耗为
通过有限元计算分析得到的电机铁心某处的磁感应强度是以切向分量和径向分量的形式给出的,假设径向分量用BN表示,切向分量用BT表示,则第m次分量的可表示为
则考虑切向和径向分量的第m次谐波产生的铁耗为
因此计及前10次谐波铁心中某一部分铁损计算模型为
对电机每一部分铁损运用上述计及前10次谐波铁损计算模型计算得到电机铁损,求和即可得到电机的铁损模型。以一台四相8/6结构开关磁阻电机为例根据磁感应强度的变化情况和大小对电机进行分块图如图9所示,其中ST表示定子齿部,SYAB、SYBC、SYCD、SYDA分别表示定子轭的四部分,RY表示转子轭部分,RT表示转子齿部分。
根据上述分块区域及其铁心模型对电机铁耗进行计算,当转矩保持为恒定值0.4N·m不变时,四相8/6结构开关磁阻电机在不同转速下的铁耗计算结果如表1所示。
表1不同转速下四相8/6结构开关磁阻电机铁耗
当转矩速保持为恒定值500r/min不变时,四相8/6结构开关磁阻电机在不同转矩下的铁耗计算结果如表2所示。
表2不同转矩下四相8/6结构SRM铁耗
Claims (1)
1.一种改进的开关磁阻电机铁耗计算模型建立方法,其特征在于,该方法是在通过对不同频率和磁场强度下铁磁材料实测的损耗值进行数据分析,并结合Steinmetz公式,得到Steinmetz公式中磁滞耗系数kh、涡流损耗系数ke和Steinmetz因子α这几个关键参数随磁感应强度和频率变化的规律,从而得到了正弦交变磁化时铁磁材料铁耗计算的变系数两项式计算模型;针对开关磁阻电机中磁感应强度波形非正弦的特殊情况,运用有限元计算得到电机每一部分的切向和径向磁感应强度随时间变化波形,分别对其进行傅里叶分解将非正弦磁感应强度波形分解成一列正弦波形,对每次谐波运用变系数两项式计算模型计算得到每次谐波的铁损,累加求和即可得到计及n次谐波单位质量铁损模型,建立电机的铁损模型;具体模型建立过程如下:
在磁场为正弦交变磁场,磁感应强度峰值接近饱和以及不计高频下集肤效应对涡流的影响时,Steinmetz铁耗计算模型表达式为:
Pfe=khfBα+kef2B2 (1)
式中,B为磁感应强度,f为频率,khBαf为磁滞损耗,keB2f2为涡流损耗,kh、ke和α分别为磁滞损耗系数、涡流损耗系数和Steinmetz因子,它们的取值直接影响到计算模型的准确性;令khBα=C,keB2=D,并使左右两边同除以f则可得
根据上式求得不同频率下的磁滞损耗系数、不同磁感应强度下的磁滞损耗系数、不同频率下的涡流损耗系数、不同磁感应强度下的涡流损耗系数、不同频率下的Steinmetz因子和不同磁感应强度下的Steinmetz因子,分析磁滞损耗系数和涡流损耗系数随频率和磁感应强度的变化趋势,得到磁滞损耗系数kh的取值受频率和磁感应强度变化的影响很小,可近似认为kh的取值为常数;涡流损耗系数ke随频率的增大基本呈直线下降趋势,涡流损耗系数ke随磁感应强度的变化趋势不大,可忽略不计,因此涡流损耗系数可用以下表达式拟合:
Ke=af+b (3)
随着频率的变化Steinmetz因子基本保持不变,可以认为Steinmetz因子的取值与频率无关,Steinmetz因子随磁感应强度的增大而增大,其变化规律可用以下多项式拟合:
α=k1B3+k2B2+k3B+k4 (4)
通过以上分析可得到铁磁材料的两项式变系数计算模型如下:
以上所述两项式变系数计算模型只适用于磁感应强度波形为正弦波的情况下铁耗的计算,而在开关磁阻电机铁心中磁感应强度恰恰就不是正弦波形,因此不能直接用上述模型计算开关磁阻电机中的铁耗;因此假设磁感应强度波形的基波频率为f,经过傅立叶分解后,第m次谐波的幅值为Bm,频率为mf,则由第m次谐波产生的铁耗为
通过有限元计算分析得到的电机铁心某处的磁感应强度是以切向分量和径向分量的形式给出的,假设径向分量用BN表示,切向分量用BT表示,则第m次分量的可表示为:
则考虑切向和径向分量的第m次谐波产生的铁耗为
因此计及前n次谐波铁心中某一部分铁损计算模型为
对电机每一部分铁损运用上述计及前n次谐波铁损计算模型计算得到电机铁损,求和即可得到电机的铁损模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610804714.7A CN106326595B (zh) | 2016-09-06 | 2016-09-06 | 一种改进的开关磁阻电机铁耗计算模型建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610804714.7A CN106326595B (zh) | 2016-09-06 | 2016-09-06 | 一种改进的开关磁阻电机铁耗计算模型建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106326595A CN106326595A (zh) | 2017-01-11 |
CN106326595B true CN106326595B (zh) | 2019-05-31 |
Family
ID=57787877
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610804714.7A Active CN106326595B (zh) | 2016-09-06 | 2016-09-06 | 一种改进的开关磁阻电机铁耗计算模型建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106326595B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112560301A (zh) * | 2020-11-26 | 2021-03-26 | 东南大学 | 一种基于磁感原理的磁性材料涡流损耗计算方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107862167B (zh) * | 2017-12-14 | 2021-04-30 | 中国矿业大学 | 一种考虑瞬时铁耗的开关磁阻电机建模方法 |
CN108595902A (zh) * | 2018-07-16 | 2018-09-28 | 山东大学 | 计及谐波影响的直线电机铁耗计算方法及模型 |
CN111044956B (zh) * | 2019-11-28 | 2021-05-11 | 浙江大学 | 一种磁滞损耗的估算方法 |
CN111967170B (zh) * | 2020-06-15 | 2023-07-04 | 中国矿业大学 | 一种开关磁阻电机铁芯损耗计算模型建立方法 |
CN112540330B (zh) * | 2020-11-26 | 2021-10-08 | 东南大学 | 一种基于磁感原理的磁性材料b-h曲线测量方法 |
CN112966371B (zh) * | 2021-02-08 | 2021-11-02 | 华北电力大学(保定) | 交直流混合激励下的铁磁材料的异常损耗计算方法 |
CN113049998B (zh) * | 2021-02-08 | 2021-11-02 | 华北电力大学(保定) | 多谐波激励作用下铁磁材料损耗预测方法 |
CN113111556B (zh) * | 2021-04-15 | 2022-06-17 | 浙江大学 | 一种超重力环境下考虑磁场与应力影响的感应电机铁耗分析方法 |
CN113692078A (zh) * | 2021-07-08 | 2021-11-23 | 华南理工大学 | 一种磁性复合材料的固化系统及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103197267A (zh) * | 2013-02-28 | 2013-07-10 | 南京邮电大学 | 一种磁芯损耗分离计算方法 |
CN105205234A (zh) * | 2015-09-09 | 2015-12-30 | 华北电力大学 | 一种用于交流电机损耗精细化分析的分段变系数铁耗模型 |
-
2016
- 2016-09-06 CN CN201610804714.7A patent/CN106326595B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103197267A (zh) * | 2013-02-28 | 2013-07-10 | 南京邮电大学 | 一种磁芯损耗分离计算方法 |
CN105205234A (zh) * | 2015-09-09 | 2015-12-30 | 华北电力大学 | 一种用于交流电机损耗精细化分析的分段变系数铁耗模型 |
Non-Patent Citations (3)
Title |
---|
偏转磁芯中磁滞损耗的估测;张晓兵 等;《真空科学与技术》;20010331;第21卷(第2期);全文 * |
开关磁阻电机磁通波形与铁秏计算;丁文 等;《电机与控制应用》;20061231;第33卷(第6期);全文 * |
转子结构对高速永磁电机转子损耗影响分析;王晓远 等;《微电机》;20160630;第49卷(第6期);全文 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112560301A (zh) * | 2020-11-26 | 2021-03-26 | 东南大学 | 一种基于磁感原理的磁性材料涡流损耗计算方法 |
CN112560301B (zh) * | 2020-11-26 | 2021-10-08 | 东南大学 | 一种基于磁感原理的磁性材料涡流损耗计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106326595A (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106326595B (zh) | 一种改进的开关磁阻电机铁耗计算模型建立方法 | |
Zhu et al. | Improved analytical modelling of rotor eddy current loss in brushless machines equipped with surface-mounted permanent magnets | |
Doering et al. | Torque, power, losses, and heat calculation of a transverse flux reluctance machine with soft magnetic composite materials and disk-shaped rotor | |
Reigosa et al. | Magnet temperature estimation in permanent magnet synchronous machines using the high frequency inductance | |
Bramerdorfer et al. | Accurate and easy-to-obtain iron loss model for electric machine design | |
Virtic et al. | Torque analysis of an axial flux permanent magnet synchronous machine by using analytical magnetic field calculation | |
Marfoli et al. | Rotor design optimization of squirrel cage induction motor-part i: Problem statement | |
Cham et al. | Brushless dc motor electromagnetic torque estimation with single-phase current sensing | |
CN108875168A (zh) | 一种考虑饱和的开关磁阻电机磁场解析计算方法 | |
Kumar et al. | Electromagnetic field analysis of 3-Phase induction motor drive under broken rotor bar fault condition using FEM | |
Chung et al. | Modeling of the armature slotting effect in the magnetic field distribution of a linear permanent magnet motor | |
CN106339563B (zh) | 一种开关磁阻电机铁耗快速计算方法 | |
Caruso et al. | Interior permanent magnet synchronous motors: Impact of the variability of the parameters on their efficiency | |
Cox et al. | End turn leakage reactance of concentrated modular winding stators | |
Verma et al. | Analytical prediction of static flux-linkage characteristics of switched reluctance machine | |
Somesan et al. | Two simple analytical models, direct and inverse, for switched reluctance motors | |
Zheng et al. | Thrust performance improvement of PMSLM based on lasso regression with embedded analytical model | |
Mekuria | Development of a high speed solid rotor asynchronous drive fed by a frequency converter system | |
Qiu et al. | The influence of inter-turn short circuit fault considering loop current on the electromagnetic field of permanent magnet synchronous motor | |
Ojeda et al. | Comparison of 3-phase and 5-phase flux switching machine using a stochastic approach | |
Taghizadeh Kakhki | Modeling of losses in a permanent magnet machine fed by a PWM supply | |
Ji et al. | Multi-objective optimization of interior permanent magnet machine for heavy-duty vehicle direct-drive applications | |
Yanli et al. | Analytical calculation for predicting the core loss of surface-mounted permanent magnet machine | |
Yetgin et al. | Effects of rotor slot area on squirrel cage induction motor performance | |
CN105676125B (zh) | 变频电机功率测量方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |