CN111190128A - 一种磁阻电机铁磁材料bh特性检测算法 - Google Patents
一种磁阻电机铁磁材料bh特性检测算法 Download PDFInfo
- Publication number
- CN111190128A CN111190128A CN201811357853.5A CN201811357853A CN111190128A CN 111190128 A CN111190128 A CN 111190128A CN 201811357853 A CN201811357853 A CN 201811357853A CN 111190128 A CN111190128 A CN 111190128A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- characteristic
- stator
- magnetic field
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/12—Measuring magnetic properties of articles or specimens of solids or fluids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/12—Measuring magnetic properties of articles or specimens of solids or fluids
- G01R33/14—Measuring or plotting hysteresis curves
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Iron Core Of Rotating Electric Machines (AREA)
- Synchronous Machinery (AREA)
Abstract
本发明属于电机铁磁材料技术领域,特别涉及一种磁阻电机铁磁材料BH特性检测算法。将开关磁阻电机对齐位置磁链特性划分为两个区域,区域I通过磁路的基尔霍夫定律反推出线性区BH特性,区域II采用递推法与安培环路定律计算出非线性区BH特性。本方法计算速度快,方法简单,易于实现,仅需要对齐位置磁链特性;适用性好,三相12/8极SRM与四相8/6极SRM均得到了较为准确的结果;利用计算得到的BH特性曲线可以大幅度提高仿真精度。
Description
技术领域
本发明属于电机铁磁材料技术领域,特别涉及一种磁阻电机铁磁材料BH特性检测算法。
背景技术
电机的电磁特性受到电机尺寸与材料的影响。当前加工工艺精度较高,使得电机尺寸加工引入的误差较少。而大量的研究表明,电机的材料在加工成硅钢片时,加工过程会对材料特性产生影响。例如,机械切割产生的塑性变型和激光切割产生的热应力,这些都会使电机实际的BH特性与原本的材料特性不同。因此,在使用有限元软件对电机电磁特性进行仿真时,会出现仿真和实验获得的电磁特性不能很好的吻合的情况。为了能够更加精确地分析电机的性能,在利用有限元法对电机的性能进行分析之前,通过实验或是其他方法获得电机材料特性是十分必要的。
发明内容
针对目前存在的问题,提出一种磁阻电机铁磁材料BH特性检测算法。
为了实现这一目的,本发明采取的技术方案是:
一种磁阻电机铁磁材料BH特性检测算法,将开关磁阻电机对齐位置磁链特性划分为两个区域,区域I通过磁路的基尔霍夫定律反推出线性区BH特性,区域II采用递推法与安培环路定律计算出非线性区BH特性,步骤一:测量对齐位置导通相电压、电流值,利用式(1)计算得到对齐位置磁链数据ψa;
其中,ψ(0)是初始磁链,由于硅钢材料剩磁较小,通常将ψ(0)取为0;u、i和r分别为开关磁阻电机的相电压、相电流和相电阻;
步骤二:利用辅助线将开关磁阻电机对齐位置磁链特性与相电流关系划分为两个区域,分别命名为非饱和区与饱和区;
步骤三:在非饱和区时,磁场未饱和,磁路中所有铁芯部分的磁导率是相同的,利用式(2)计算;
式中,uu为线性区磁导率,ug为气隙磁导率,Asy为定子磁轭截面积,Asp为定子磁极截面积,Arp为转子磁极截面积,Ary为转子磁轭截面积,Ag为气隙截面积,lsy为磁路通过的定子磁轭长度,lsp为磁路通过的定子磁极的长度,lrp为磁路通过的转子磁极长度,lry为磁路通过的转子磁轭长度,lg为磁路通过的气隙长度,φ为磁通,N为电机绕组匝数,i为相电流;
步骤四:通过式(3),利用对齐位置磁链特性推导出Bsp,并利用步骤三计算得到的磁导率推到出线性区的磁场强度,计算公式见式(4);
式中,Bsp为定子磁密,Asp为定子截面积,N为绕组匝数,ψ为磁链,Hsp为磁场强度,uu为线性区磁导率;
步骤五:利用靠近线性区的第一个磁链值计算得到Bsp,随后推导出其他位置的磁密与磁场强度,见式(5)~(12);计算当前磁链值下的定子磁密,从而使得电机其他部分的磁密处于已知的BH特性中;
式中,Bsp为定子磁密,Asp为定子截面积,Bsy为定子磁轭磁密,Asy为定子磁轭截面积,Bry为转子磁轭磁密,Ary为转子磁轭截面积,Brp为转子磁密,Arp为转子截面积,Bg为气隙磁密,Ag为气隙截面积,Hsy为定子磁轭磁场强度,Hry为转子磁轭磁场强度,Hrp为转子磁场强度,Hg为气隙磁场强度,Ag为磁路穿过气隙的截面积,ug为气隙磁导率,Φ为磁通;
步骤六:利用式(13)计算得到定子极磁场强度,
式中,lsy为磁路通过的定子磁轭长度,lsp为磁路通过的定子磁极的长度,lrp为磁路通过的转子磁极长度,lry为磁路通过的转子磁轭长度,lg为磁路通过的气隙长度,N为电机绕组匝数,i为相电流;Hsy为定子磁轭磁场强度,Hry为转子磁轭磁场强度,Hrp为转子磁场强度,Hg为气隙磁场强度;
步骤七:不断重复步骤五与步骤六,并逐步选择下一个靠近计算得到的已知区域的BH特性,计算出电机的BH材料特性;处于非线性区的计算得到的已知BH特性,可通过线性插值得到磁场强度,用于式(13)的计算,非线性区磁场强度计算公式如(14)所示;式中,B处于区间[Bsp0,Bsp1]中,其他区间可通过类似公式进行推导;应注意式(14)中的Hsp0,Hsp1,Bsp0,Bsp1是根据B所处区间进行选择的;
一种磁阻电机铁磁材料BH特性检测算法,所属步骤五中,开关磁阻电机定子磁密通常为最大值,选取靠近线性区的磁链特性。
本发明的有益效果为:
1.方法计算速度快。
2.方法简单,易于实现。仅需要对齐位置磁链特性;
3.适用性好。三相12/8极SRM与四相8/6极SRM均得到了较为准确的结果;
4.利用计算得到的BH特性曲线可以大幅度提高仿真精度。
附图说明
图1为对齐位置磁链特性线性区与非线性区划分;
图2为图1中对应P0,P1与P2的计算结果。
具体实施方式
本发明具体实施例如下:
一种磁阻电机铁磁材料BH特性检测算法,将开关磁阻电机对齐位置磁链特性划分为两个区域,区域I为非饱和区域,通过磁路的基尔霍夫定律反推出线性区BH特性,区域II为饱和区域,采用递推法与安培环路定律计算出非线性区BH特性。具体步骤如下:
步骤一:测量对齐位置导通相电压、电流值,利用式(1)计算得到对齐位置磁链数据ψa。
其中,ψ(0)是初始磁链,由于硅钢材料剩磁较小,通常将ψ(0)取为0;u、i和r分别为开关磁阻电机的相电压、相电流和相电阻。
步骤二:利用辅助线将开关磁阻电机对齐位置磁链特性与相电流关系划分为两个区域,分别命名为非饱和区与饱和区。
步骤三:在非饱和区时,磁场未饱和,磁路中所有铁芯部分的磁导率是相同的,利用式(2)计算。
式中,uu为线性区磁导率,ug为气隙磁导率,Asy为定子磁轭截面积,Asp为定子磁极截面积,Arp为转子磁极截面积,Ary为转子磁轭截面积,Ag为气隙截面积,lsy为磁路通过的定子磁轭长度,lsp为磁路通过的定子磁极的长度,lrp为磁路通过的转子磁极长度,lry为磁路通过的转子磁轭长度,lg为磁路通过的气隙长度,φ为磁通,N为电机绕组匝数,i为相电流。
步骤四:通过式(3),利用对齐位置磁链特性推导出Bsp,并利用步骤三计算得到的磁导率推到出线性区的磁场强度,计算公式见式(4)。
式中,Bsp为定子磁密,Asp为定子截面积,N为绕组匝数,ψ为磁链,Hsp为磁场强度,uu为线性区磁导率。
步骤五:利用靠近线性区的第一个磁链值计算得到Bsp,随后推导出其他位置的磁密与磁场强度,见式(5)~(12)。应注意,开关磁阻电机定子磁密通常为最大值,要选取靠近线性区的磁链特性,计算当前磁链值下的定子磁密,从而使得电机其他部分的磁密处于已知的BH特性中。
式中,Bsp为定子磁密,Asp为定子截面积,Bsy为定子磁轭磁密,Asy为定子磁轭截面积,Bry为转子磁轭磁密,Ary为转子磁轭截面积,Brp为转子磁密,Arp为转子截面积,Bg为气隙磁密,Ag为气隙截面积,Hsy为定子磁轭磁场强度,Hry为转子磁轭磁场强度,Hrp为转子磁场强度,Hg为气隙磁场强度,Ag为磁路穿过气隙的截面积,ug为气隙磁导率,Φ为磁通。
步骤六:利用式(13)计算得到定子极磁场强度
式中,lsy为磁路通过的定子磁轭长度,lsp为磁路通过的定子磁极的长度,lrp为磁路通过的转子磁极长度,lry为磁路通过的转子磁轭长度,lg为磁路通过的气隙长度,N为电机绕组匝数,i为相电流。Hsy为定子磁轭磁场强度,Hry为转子磁轭磁场强度,Hrp为转子磁场强度,Hg为气隙磁场强度。
步骤七:不断重复步骤五与步骤六,并逐步选择下一个靠近计算得到的已知区域的BH特性,即可计算出电机的BH材料特性。处于非线性区的计算得到的已知BH特性,可通过线性插值得到磁场强度,用于式(13)的计算,非线性区磁场强度计算公式如(14)所示。式中,B处于区间[Bsp0,Bsp1]中,其他区间可通过类似公式进行推导。应注意式(14)中的Hsp0,Hsp1,Bsp0,Bsp1是根据B所处区间进行选择的。
Claims (9)
1.一种磁阻电机铁磁材料BH特性检测算法,将开关磁阻电机对齐位置磁链特性划分为两个区域,区域I为非饱和区域,通过磁路的基尔霍夫定律反推出线性区BH特性,区域II为饱和区域,采用递推法与安培环路定律计算出非线性区BH特性,其特征在于:本方法包括七个步骤,步骤一,测量电压电流、并计算对齐位置磁链数据,步骤二,利用辅助线划分区域,步骤三,非饱和区磁导率计算,步骤四,计算线性区的磁场强度,步骤五,计算当前磁链值下的定子磁密,步骤六,计算定子极磁场强度,步骤七,计算出电机的BH材料特性。
3.如权利要求1所述的一种磁阻电机铁磁材料BH特性检测算法,其特征在于:所述步骤二:利用辅助线将开关磁阻电机对齐位置磁链特性与相电流关系划分为两个区域,分别命名为非饱和区与饱和区。
6.如权利要求1所述的一种磁阻电机铁磁材料BH特性检测算法,其特征在于:所述步骤五:利用靠近线性区的第一个磁链值计算得到Bsp,随后推导出其他位置的磁密与磁场强度,见式(5)~(12);计算当前磁链值下的定子磁密,从而使得电机其他部分的磁密处于已知的BH特性中;
式中,Bsp为定子磁密,Asp为定子截面积,Bsy为定子磁轭磁密,Asy为定子磁轭截面积,Bry为转子磁轭磁密,Ary为转子磁轭截面积,Brp为转子磁密,Arp为转子截面积,Bg为气隙磁密,Ag为气隙截面积,Hsy为定子磁轭磁场强度,Hry为转子磁轭磁场强度,Hrp为转子磁场强度,Hg为气隙磁场强度,Ag为磁路穿过气隙的截面积,ug为气隙磁导率,Φ为磁通。
9.如权利要求6所述的一种磁阻电机铁磁材料BH特性检测算法,其特征在于:所属步骤五中,开关磁阻电机定子磁密通常为最大值,选取靠近线性区的磁链特性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811357853.5A CN111190128B (zh) | 2018-11-15 | 2018-11-15 | 一种磁阻电机铁磁材料bh特性检测算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811357853.5A CN111190128B (zh) | 2018-11-15 | 2018-11-15 | 一种磁阻电机铁磁材料bh特性检测算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111190128A true CN111190128A (zh) | 2020-05-22 |
CN111190128B CN111190128B (zh) | 2022-10-18 |
Family
ID=70705295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811357853.5A Active CN111190128B (zh) | 2018-11-15 | 2018-11-15 | 一种磁阻电机铁磁材料bh特性检测算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111190128B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112540330A (zh) * | 2020-11-26 | 2021-03-23 | 东南大学 | 一种基于磁感原理的磁性材料b-h曲线测量方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101769998A (zh) * | 2010-01-18 | 2010-07-07 | 南通大学 | 开关磁阻电机磁链特性测量方法 |
CN102253352A (zh) * | 2011-04-22 | 2011-11-23 | 襄樊五二五泵业有限公司 | 一种磁性联轴器漏磁系数检测方法及装置 |
US20130049751A1 (en) * | 2010-04-14 | 2013-02-28 | Peter Hamberger | Method and apparatus for detecting a magnetic characteristic variable in a core |
CN104063556A (zh) * | 2014-07-07 | 2014-09-24 | 东南大学 | 一种轴向永磁电机等效磁路模型的建模方法 |
CN104316876A (zh) * | 2014-10-08 | 2015-01-28 | 西北工业大学 | 一种考虑互感耦合的获取三相12/8极srm磁链特性的快速方法 |
CN104330755A (zh) * | 2014-08-26 | 2015-02-04 | 南通大学 | 一种开关磁阻电机磁链特性在线测量方法 |
CN104539212A (zh) * | 2015-01-20 | 2015-04-22 | 中国矿业大学 | 一种开关磁阻电机磁链特性在线测量方法 |
CN104849587A (zh) * | 2015-04-30 | 2015-08-19 | 国网四川省电力公司电力科学研究院 | 一种分析变压器在直流偏磁影响下的励磁特性变化方法 |
CN106872895A (zh) * | 2017-04-25 | 2017-06-20 | 天津工业大学 | 一种开关磁阻电机磁链特性检测装置及检测方法 |
CN107425781A (zh) * | 2016-01-06 | 2017-12-01 | 西北工业大学 | 一种基于线性磁链模型和线性回归分析的srm位置预估方法 |
CN108282125A (zh) * | 2018-01-30 | 2018-07-13 | 深圳市默贝克驱动技术有限公司 | 一种基于磁密分布特征的凸极电机磁链计算方法 |
CN108595772A (zh) * | 2018-03-26 | 2018-09-28 | 江苏大学 | 一种定转子双永磁游标电机的2d磁路剖分建模方法 |
-
2018
- 2018-11-15 CN CN201811357853.5A patent/CN111190128B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101769998A (zh) * | 2010-01-18 | 2010-07-07 | 南通大学 | 开关磁阻电机磁链特性测量方法 |
US20130049751A1 (en) * | 2010-04-14 | 2013-02-28 | Peter Hamberger | Method and apparatus for detecting a magnetic characteristic variable in a core |
CN102253352A (zh) * | 2011-04-22 | 2011-11-23 | 襄樊五二五泵业有限公司 | 一种磁性联轴器漏磁系数检测方法及装置 |
CN104063556A (zh) * | 2014-07-07 | 2014-09-24 | 东南大学 | 一种轴向永磁电机等效磁路模型的建模方法 |
CN104330755A (zh) * | 2014-08-26 | 2015-02-04 | 南通大学 | 一种开关磁阻电机磁链特性在线测量方法 |
CN104316876A (zh) * | 2014-10-08 | 2015-01-28 | 西北工业大学 | 一种考虑互感耦合的获取三相12/8极srm磁链特性的快速方法 |
CN104539212A (zh) * | 2015-01-20 | 2015-04-22 | 中国矿业大学 | 一种开关磁阻电机磁链特性在线测量方法 |
CN104849587A (zh) * | 2015-04-30 | 2015-08-19 | 国网四川省电力公司电力科学研究院 | 一种分析变压器在直流偏磁影响下的励磁特性变化方法 |
CN107425781A (zh) * | 2016-01-06 | 2017-12-01 | 西北工业大学 | 一种基于线性磁链模型和线性回归分析的srm位置预估方法 |
CN106872895A (zh) * | 2017-04-25 | 2017-06-20 | 天津工业大学 | 一种开关磁阻电机磁链特性检测装置及检测方法 |
CN108282125A (zh) * | 2018-01-30 | 2018-07-13 | 深圳市默贝克驱动技术有限公司 | 一种基于磁密分布特征的凸极电机磁链计算方法 |
CN108595772A (zh) * | 2018-03-26 | 2018-09-28 | 江苏大学 | 一种定转子双永磁游标电机的2d磁路剖分建模方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112540330A (zh) * | 2020-11-26 | 2021-03-23 | 东南大学 | 一种基于磁感原理的磁性材料b-h曲线测量方法 |
US11965942B2 (en) | 2020-11-26 | 2024-04-23 | Southeast University | Measurement method for B-H curve of magnetic material based on magnetic-inductance |
Also Published As
Publication number | Publication date |
---|---|
CN111190128B (zh) | 2022-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bianchi et al. | Experimental comparison of PM-assisted synchronous reluctance motors | |
Azzouzi et al. | Quasi-3-D analytical modeling of the magnetic field of an axial flux permanent-magnet synchronous machine | |
Bianchi et al. | Effect of stator and rotor saturation on sensorless rotor position detection | |
Walker et al. | Flux-linkage calculation in permanent-magnet motors using the frozen permeabilities method | |
Bianchi et al. | An analytical approach to design the PM in PMAREL motors robust toward the demagnetization | |
CN108563912B (zh) | 一种表贴式永磁同步电机气隙磁场的解析方法 | |
JP5228730B2 (ja) | 永久磁石式モータの磁石内渦電流損失解析方法 | |
Kim et al. | Inductance calculation in IPMSM considering magnetic saturation | |
CN107565782A (zh) | 一种混合转子永磁同步电机的等效磁网络分析方法 | |
Li et al. | Open-circuit field prediction of interior permanent-magnet motor using hybrid field model accounting for saturation | |
Alberti et al. | High-frequency $ d $–$ q $ model of synchronous machines for sensorless control | |
Alberti et al. | Finite-element analysis of electrical machines for sensorless drives with high-frequency signal injection | |
Lin et al. | A new nonlinear anisotropic model for soft magnetic materials | |
Kellner et al. | Impact of iron losses on parameter identification of permanent magnet synchronous machines | |
CN111190128B (zh) | 一种磁阻电机铁磁材料bh特性检测算法 | |
Zhou et al. | Determination of $ d {\hbox {-}} q $ Axis Parameters of Interior Permanent Magnet Machine | |
Yang et al. | Analysis of iron losses in induction motor with an improved iron-loss model | |
CN112906259B (zh) | 基于精确保角映射的永磁游标电机电磁性能分析方法 | |
Lee et al. | Computation of inductance and static thrust of a permanent-magnet-type transverse flux linear motor | |
Faradonbeh et al. | Analytical modeling of steel shielding cylinder in surface mounted permanent magnet machines | |
Si et al. | Development of the equivalent magnetic circuit model for a surface-interior permanent magnet synchronous motor | |
Lu et al. | Modelling a single phase surface mounted permanent magnet transverse flux machine based on Fourier series method | |
Prieto et al. | Analytical model for a saturated permanent magnet assisted synchronous reluctance motor | |
Schulte et al. | Co-simulation of an interior permanent magnet synchronous motor with segmented rotor structure | |
Lim et al. | Characteristic analysis of 5-phase hybrid stepping motor considering the saturation effect |
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 |