CN110880894A - 一种永磁电机pwm谐波损耗的快速计算方法 - Google Patents
一种永磁电机pwm谐波损耗的快速计算方法 Download PDFInfo
- Publication number
- CN110880894A CN110880894A CN201911064620.0A CN201911064620A CN110880894A CN 110880894 A CN110880894 A CN 110880894A CN 201911064620 A CN201911064620 A CN 201911064620A CN 110880894 A CN110880894 A CN 110880894A
- Authority
- CN
- China
- Prior art keywords
- permanent magnet
- loss
- magnet motor
- harmonic
- rotor
- 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.)
- Granted
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Iron Core Of Rotating Electric Machines (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
本发明公开了一种永磁电机PWM谐波损耗的快速计算方法。首先采用静态有限元计算得到永磁电机硅钢片中的磁密分布;然后基于冻结增量张量磁阻率法构建永磁电机交流小信号模型,采用频域有限元法计算d、q轴高频谐波电压分别激励下,定、转子上的损耗因子,实现直接分别利用定、转子上的PWM电压频谱,快速计算定、转子上的PWM谐波铁耗;最后通过二次多项式拟合涡流损耗因子随基波电流的变化关系,实现永磁电机在整个工作范围内PWM谐波铁耗的快速计算。本发明克服了现有技术耗时多的缺点,且兼具准确性与快速性。
Description
技术领域
本发明属于永磁电机领域,特别涉及了一种永磁电机PWM谐波损耗的计算方法。
背景技术
为实时调节永磁电机的输出功率和转速,永磁电机通常都是与PWM(Pulse WidthModulated,脉宽调制)电压源型逆变器配合使用的。而另一方面PWM电压源型逆变器所产生的电压谐波会在硅钢片中感应出额外的PWM谐波铁耗,从而增加电机的损耗,引起效率下降和温升提高。因此设计人员在设计电机时需要对这一现象进行准确的建模计算,为优化设计奠定基础。
目前通过以PWM电压源为输入的场路耦合有限元计算、并结合解析铁耗模型的方法需要采用很小的步长来分辨高频PWM谐波电压,这使得以该方法在计算铁耗时非常耗时。
发明内容
为了解决上述背景技术提到的技术问题,本发明提出了一种永磁电机PWM谐波损耗的快速计算方法。
为了实现上述技术目的,本发明的技术方案为:
一种永磁电机PWM谐波损耗的快速计算方法,包括以下步骤:
(1)对于永磁电机,已知其在某个工况下所需的d、q轴基波电流为id0和iq0,以d、q轴基波电流为输入,在转子位置角θe处进行一步静态有限元计算;
(2)基于冻结张量增量磁导率法构建永磁电机线性交流小信号模型;
(6)改变id0和iq0,重复步骤(1)-(5),计算不同工况、不同基波电流下的四个损耗因子的平均值;
(7)采用二次多项式拟合四个损耗因子随基波电流的变化关系;
(8)对于永磁电机转矩-转速图上的任意一个工作点,首先确定在此工作点的id0和iq0,代入步骤(7)中的拟合公式计算该工况下的损耗因子;
(9)对于步骤(8)中的工作点,根据逆变器参数和调制策略进一步确定输入永磁电机的PWM谐波频谱,并通过坐标变换得到定子坐标系下的PWM电压频谱和转子坐标系下的电压频谱;
(10)分别利用定、转子坐标系下的PWM电压频谱,计算各次谐波在定、转子上产生的PWM谐波损耗;
(11)将除基波分量之外的所有谐波分量所产生的损耗叠加,分别得到定子和转子上总的谐波损耗;
(12)对于永磁电机转矩-转速图上的各个工作点,重复步骤(8)-(11)得到所有工作点的PWM谐波损耗。
上式中,ke为硅钢片涡流损耗系数,Br,r和Br,i为硅钢片中在两倍开关频率2fc、幅值任意的d轴谐波电压udh激励下的交变磁密径向分量的实部与虚部,Bθ,r和Bθ,i分别为在两倍开关频率2fc、幅值任意的d轴谐波电压udh激励下的硅钢片中交变磁密切向分量的实部与虚部,stator表示定子区域,rotor表示转子区域,dv表示体积微元,|udh|为udh的幅值。
上式中,ke为硅钢片涡流损耗系数,Br,r和Br,i为硅钢片中在两倍开关频率2fc、幅值任意的q轴谐波电压uqh激励下的交变磁密径向分量的实部与虚部,Bθ,r和Bθ,i分别为在两倍开关频率2fc、幅值任意的q轴谐波电压uqh激励下的硅钢片中交变磁密切向分量的实部与虚部,stator表示定子区域,rotor表示转子区域,dv表示体积微元,|uqh|为uqh的幅值。
进一步地,在步骤(6)中,对于采用id0=0控制的表贴式永磁电机,计算iq0为0、iq0为最大电流以及iq0为1/2最大电流这三种工况。
上式中,a1、a2和a3为待拟合系数,由步骤(6)中三种工况计算得到的损耗因子拟合得到。
进一步地,在步骤(6)中,对于内嵌式永磁同步电机电机,计算空载、最大转矩、1/2最大转矩、1/4最大转矩、最高转速空载和1/2最高转速满载这六种工况。
进一步地,在步骤(7)中,对于内嵌式永磁电机采用的拟合关系式f(Im,α)如下:
上式中,Im为电流幅值,α为电流相角,b1、b2、b3、b4、b5和b6为待拟合系数,由步骤(6)中六种工况计算得到的损耗因子拟合得到。
进一步地,在步骤(10)中,定、转子上产生的PWM谐波损耗的计算式如下:
上式中,Uαm和Uβm为定子坐标系下m次谐波电压的幅值,为对应频率,为该次谐波电压在定子上产生的谐波损耗,kh为硅钢片磁滞损耗系数,ke为硅钢片涡流损耗系数,Udm和Uqm为转子坐标系下m次谐波电压的幅值,为对应频率,为该次谐波电压在转子中产生的谐波损耗。
采用上述技术方案带来的有益效果:
本发明设计的计算方法无需在每个工作点都采用PWM电压供电下的瞬态有限元计算来计算硅钢片中的高频PWM谐波铁耗,而只需要在在每个工作点重复数十次线性频域有限元计算即可。通过频域有限元计算结果,即可计算PWM谐波电压与对应谐波损耗之间关系,进而实现分别利用定、转坐标系下的电压频谱快速计算定、转子上的PWM谐波损耗。此外,表征谐波电压与谐波损耗之间关系的损耗因子随基波电流之间的变化关系可用简单二次多项式进行拟合,这样在计算PWM谐波铁耗图时,可直接利用拟合的公式计算损耗因子,而无需再重复线性频域有限元计算,这大大提高了计算速度。对于单一工作点的PWM谐波损耗计算,计算时间缩短了数十倍;而对于PWM谐波铁耗图的计算,计算时间缩短了数百倍。
附图说明
图1是id0=-64.9A,iq0=76.04A时,频率为10kHz的Udh或Uqh单独激励下,定子中涡流铁耗随转子位置变化关系图;
图2是id0=-64.9A,iq0=76.04A时,频率为10kHz的Udh或Uqh单独激励下,转子中涡流铁耗随转子位置变化关系图;
图3是id0=-64.9A,iq0=76.04A,fc=4.8kHz,转子转速为2000r/min,直流母线电压为650V,采用空间矢量PWM调制策略时,定子坐标系下的PWM电压频谱图;
图4是id0=-64.9A,iq0=76.04A,fc=4.8kHz,转子转速为2000r/min,直流母线电压为650V,采用空间矢量PWM调制策略时,转子坐标系下的PWM电压频谱图;
图5是Prius2010中内嵌式永磁电机在ke=0.585(W/m3/Hz2/T2),kh=140(W/m3/Hz/T2)时计算得到的PWM谐波铁耗分布图。
具体实施方式
以下将结合附图,对本发明的技术方案进行详细说明。
本实施例以在Prius2010中应用的内嵌式永磁同步电机为例进行计算。
快速计算其硅钢片中PWM谐波铁耗的方法方法,包括如下具体步骤:
步骤1,以id0=-64.9A,iq0=76.04A为输入,进行转子电角度θe从0变化到360°的静态有限元计算,得到每个网格的磁密。
步骤2,根据上步计算结果,采用冻结张量增量磁导率法构建θe位于不同转子位置时的永磁电机线性交流小信号模型。
步骤3,利用构建的永磁电机线性交流小信号模型,采用线性频域有限元法计算硅钢片在频率10kHz、幅值为100V的udh和uqh分别激励下,硅钢片中交变磁密。由此计算得到在θe位于不同位置时,100V的udh和uqh分别激励下,定子上的涡流铁耗如图1所示,转子上的涡流铁耗如图2所示。
步骤4,分析图1和2,得出相同谐波电压激励下得到的涡流铁耗随转子电角度变化周期为60°,在一个周期内取15个点的涡流铁耗,将其除以(100V)2即可到不同转子位置的涡流损耗因子。
步骤6,改变id0和iq0,重复步骤1-5,计算空载、最大转矩、1/2最大转矩、1/4最大转矩、最高转速空载和1/2最高转速满载这六个工况时对应的平均损耗因子。
步骤7,根据上步得到的六个结果拟合四个损耗因子随着基波电流的变化规律,所采用的拟合关系式为:
上式中,Im为电流幅值,α为电流相角,b1、b2、b3、b4、b5和b6为待拟合系数,由步骤6中六种工况计算得到的损耗因子拟合得到。
步骤8,对于永磁电机转矩-转速图上的任意一个工作点,首先确定在此工作点的id0和iq0,带入步骤7中的拟合公式计算该工况下的损耗因子。
步骤9,对于步骤8中的工作点,根据逆变器参数及调制策略进一步确定输入永磁电机的PWM谐波频谱,并通过坐标变换得到,定子坐标系下的PWM电压频谱如图3所示,转子坐标系下的PWM电压频谱如图4所示。
步骤10,分别利用定、转子坐标系下的PWM电压频谱,计算各次谐波在定、转子上产生的PWM谐波损耗,具体计算公式为:
其中,Uαm和Uβm为定子坐标系下m次谐波电压的幅值,为对应频率,为该次谐波电压在定子上产生的谐波损耗,kh为硅钢片磁滞损耗系数,Udm和Uqm为转子坐标系下m次谐波电压的幅值,为对应频率,为该次谐波电压在转子中产生的谐波损耗;
步骤11,将除基波分量之外的所有谐波分量所产生的损耗叠加,即可分别得到定、转子上总的谐波损耗;
步骤12,对于永磁电机转矩-转速图上的各个工作点,重复步骤8-11即可得到,所有工作点的PWM谐波损耗,进而可以绘制出永磁电机PWM谐波铁耗图,如图5所示。
上述实施例是以内嵌式永磁同步电机为例进行说明的,本发明并不限于此种永磁电机,本发明还适用于表贴式永磁电机等其他永磁电机。
实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (8)
1.一种永磁电机PWM谐波损耗的快速计算方法,其特征在于,包括以下步骤:
(1)对于永磁电机,已知其在某个工况下所需的d、q轴基波电流为id0和iq0,以d、q轴基波电流为输入,在转子位置角θe处进行一步静态有限元计算;
(2)基于冻结张量增量磁导率法构建永磁电机线性交流小信号模型;
(6)改变id0和iq0,重复步骤(1)-(5),计算不同工况、不同基波电流下的四个损耗因子的平均值;
(7)采用二次多项式拟合四个损耗因子随基波电流的变化关系;
(8)对于永磁电机转矩-转速图上的任意一个工作点,首先确定在此工作点的id0和iq0,代入步骤(7)中的拟合公式计算该工况下的损耗因子;
(9)对于步骤(8)中的工作点,根据逆变器参数和调制策略进一步确定输入永磁电机的PWM谐波频谱,并通过坐标变换得到定子坐标系下的PWM电压频谱和转子坐标系下的电压频谱;
(10)分别利用定、转子坐标系下的PWM电压频谱,计算各次谐波在定、转子上产生的PWM谐波损耗;
(11)将除基波分量之外的所有谐波分量所产生的损耗叠加,分别得到定子和转子上总的谐波损耗;
(12)对于永磁电机转矩-转速图上的各个工作点,重复步骤(8)-(11)得到所有工作点的PWM谐波损耗。
4.根据权利要求1所述永磁电机PWM谐波损耗的快速计算方法,其特征在于,在步骤(6)中,对于采用id0=0控制的表贴式永磁电机,计算iq0为0、iq0为最大电流以及iq0为1/2最大电流这三种工况。
6.根据权利要求1所述永磁电机PWM谐波损耗的快速计算方法,其特征在于,在步骤(6)中,对于内嵌式永磁同步电机电机,计算空载、最大转矩、1/2最大转矩、1/4最大转矩、最高转速空载和1/2最高转速满载这六种工况。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911064620.0A CN110880894B (zh) | 2019-11-04 | 2019-11-04 | 一种永磁电机pwm谐波损耗的快速计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911064620.0A CN110880894B (zh) | 2019-11-04 | 2019-11-04 | 一种永磁电机pwm谐波损耗的快速计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110880894A true CN110880894A (zh) | 2020-03-13 |
CN110880894B CN110880894B (zh) | 2020-11-06 |
Family
ID=69728929
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911064620.0A Active CN110880894B (zh) | 2019-11-04 | 2019-11-04 | 一种永磁电机pwm谐波损耗的快速计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110880894B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113791345A (zh) * | 2021-08-11 | 2021-12-14 | 华中科技大学 | 一种交流铜耗计算方法 |
WO2022029911A1 (ja) * | 2020-08-05 | 2022-02-10 | 三菱電機株式会社 | モータ鉄損演算装置およびそれを備えたモータ制御装置 |
CN114157198A (zh) * | 2021-11-02 | 2022-03-08 | 河海大学 | Pwm逆变器供电下永磁同步电机定子电磁激振力的快速计算方法 |
EP4184785A4 (en) * | 2020-07-30 | 2023-08-16 | Huawei Digital Power Technologies Co., Ltd. | METHOD AND DEVICE FOR TEMPERATURE PREDICTION |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675467A (zh) * | 2013-12-26 | 2014-03-26 | 北京交通大学 | 永磁电机损耗测试方法 |
CN105205234A (zh) * | 2015-09-09 | 2015-12-30 | 华北电力大学 | 一种用于交流电机损耗精细化分析的分段变系数铁耗模型 |
JP5959772B2 (ja) * | 2014-02-17 | 2016-08-02 | 三菱電機株式会社 | 制御装置 |
JP2016158313A (ja) * | 2015-02-23 | 2016-09-01 | 三菱電機株式会社 | 制御装置、これを用いた回転電機、および当該制御装置と当該回転電機とを備えた駆動システム |
CN108595902A (zh) * | 2018-07-16 | 2018-09-28 | 山东大学 | 计及谐波影响的直线电机铁耗计算方法及模型 |
CN109687802A (zh) * | 2018-12-21 | 2019-04-26 | 河海大学 | 一种快速计算永磁同步电机中永磁体高频涡流损耗的方法 |
-
2019
- 2019-11-04 CN CN201911064620.0A patent/CN110880894B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675467A (zh) * | 2013-12-26 | 2014-03-26 | 北京交通大学 | 永磁电机损耗测试方法 |
JP5959772B2 (ja) * | 2014-02-17 | 2016-08-02 | 三菱電機株式会社 | 制御装置 |
JP2016158313A (ja) * | 2015-02-23 | 2016-09-01 | 三菱電機株式会社 | 制御装置、これを用いた回転電機、および当該制御装置と当該回転電機とを備えた駆動システム |
CN105205234A (zh) * | 2015-09-09 | 2015-12-30 | 华北电力大学 | 一种用于交流电机损耗精细化分析的分段变系数铁耗模型 |
CN108595902A (zh) * | 2018-07-16 | 2018-09-28 | 山东大学 | 计及谐波影响的直线电机铁耗计算方法及模型 |
CN109687802A (zh) * | 2018-12-21 | 2019-04-26 | 河海大学 | 一种快速计算永磁同步电机中永磁体高频涡流损耗的方法 |
Non-Patent Citations (2)
Title |
---|
YUKIHIRO YOSHIDA ET AL.: "Calculation of eddy current loss in permanent magnet motor caused by carrier harmonics based on reluctance network analysis", 《2013 15TH EUROPEAN CONFERENCE ON POWER ELECTRONICS AND APPLICATIONS (EPE)》 * |
李明宇 等: "交交变频同步电机阻尼绕组损耗的分析与计算", 《大电机技术》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4184785A4 (en) * | 2020-07-30 | 2023-08-16 | Huawei Digital Power Technologies Co., Ltd. | METHOD AND DEVICE FOR TEMPERATURE PREDICTION |
WO2022029911A1 (ja) * | 2020-08-05 | 2022-02-10 | 三菱電機株式会社 | モータ鉄損演算装置およびそれを備えたモータ制御装置 |
CN113791345A (zh) * | 2021-08-11 | 2021-12-14 | 华中科技大学 | 一种交流铜耗计算方法 |
CN113791345B (zh) * | 2021-08-11 | 2022-05-20 | 华中科技大学 | 一种交流铜耗计算方法 |
CN114157198A (zh) * | 2021-11-02 | 2022-03-08 | 河海大学 | Pwm逆变器供电下永磁同步电机定子电磁激振力的快速计算方法 |
CN114157198B (zh) * | 2021-11-02 | 2023-09-22 | 河海大学 | Pwm逆变器供电下永磁同步电机定子电磁激振力的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110880894B (zh) | 2020-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110880894B (zh) | 一种永磁电机pwm谐波损耗的快速计算方法 | |
Lee et al. | Loss distribution of three-phase induction motor fed by pulsewidth-modulated inverter | |
CN107005194A (zh) | 多绕组电动机驱动控制装置 | |
CN110362890B (zh) | 一种计及pwm谐波条件下的变频电机铁耗电阻的计算方法 | |
Zhu et al. | Fast calculation of carrier harmonic loss in permanent magnet of IPMSM under PWM VSI supply over entire working range | |
Caruso et al. | A general mathematical formulation for the determination of differential leakage factors in electrical machines with symmetrical and asymmetrical full or dead-coil multiphase windings | |
Di Tommaso et al. | Efficiency maximization of permanent magnet synchronous generators coupled coupled to wind turbines | |
CN109687802B (zh) | 一种快速计算永磁同步电机中永磁体高频涡流损耗的方法 | |
Pando-Acedo et al. | Hybrid FEA-simulink modelling of permanent magnet assisted synchronous reluctance motor with unbalanced magnet flux | |
Joksimović et al. | Stator current spectral content of inverter-fed cage rotor induction motor | |
Saied et al. | New approach to cogging torque simulation using numerical functions | |
CN111737893B (zh) | 一种基于可预测铁耗的永磁同步电机建模方法 | |
CN112016197B (zh) | 一种永磁电机谐波电流的预测方法 | |
Azer et al. | Model-based spatial harmonics vector compensation method for three-phase mutually coupled switched reluctance machine with sinusoidal current excitation | |
CN114157198B (zh) | Pwm逆变器供电下永磁同步电机定子电磁激振力的计算方法 | |
Gupta et al. | Study and Analysis of Field Oriented Control of Brushless DC Motor Drive using Hysteresis Current Control Technique | |
Bramerdorfer et al. | PMSM for high demands on low torque ripple using optimized stator phase currents controlled by an iterative learning control algorithm | |
Shisha et al. | Loss distribution on solid pole plates of wound-rotor synchronous motors fed from inverters using direct torque control | |
Di Tommaso et al. | Optimum performance of permanent magnet synchronous generators coupled to wind turbines | |
Sizov | Design synthesis and optimization of permanent magnet synchronous machines based on computationally-efficient finite element analysis | |
Chen et al. | Maximum torque per ampere control of interior permanent magnet synchronous motor via optimal current excitation | |
Wex et al. | Fully automatized pwm harmonics analysis and loss calculation in multiphase pmsm with floating starpoint | |
Shami | A novel technique to determine the required capacitance for a single phase induction motor by applying roundness algorithm | |
CN111313783A (zh) | 计算表贴式永磁电机永磁体中pwm谐波损耗的方法 | |
Dutta | A segmented interior permanent magnet synchronous machine with wide field-weakening range |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |