CN109768685A - 基于单槽模型的表贴式永磁电机电枢磁场重构方法 - Google Patents
基于单槽模型的表贴式永磁电机电枢磁场重构方法 Download PDFInfo
- Publication number
- CN109768685A CN109768685A CN201910082678.1A CN201910082678A CN109768685A CN 109768685 A CN109768685 A CN 109768685A CN 201910082678 A CN201910082678 A CN 201910082678A CN 109768685 A CN109768685 A CN 109768685A
- Authority
- CN
- China
- Prior art keywords
- slot
- magnetic field
- current
- magnetic
- field
- 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
Landscapes
- Iron Core Of Rotating Electric Machines (AREA)
- Permanent Magnet Type Synchronous Machine (AREA)
- Windings For Motors And Generators (AREA)
Abstract
本发明涉及一种基于单槽模型的表贴式永磁电机电枢磁场重构方法,包括下列步骤:第一步.建立永磁电机单槽模型;第二步.单层绕组磁场重构;第三步.双层绕组磁场重构。本发明将永磁电机总的电枢反应磁场表示为单槽通电时的气隙磁场,对单个槽磁场进行子域法求解,降低矩阵的维数,降低永磁电机磁场分析的复杂性并能够重构出与有限元计算相吻合的电枢反应磁场。
Description
技术领域
本发明属于永磁电机气隙磁场领域,涉及到基于单槽模型获取永磁体气隙磁场分布重构问题。
背景技术
永磁电机使用周期长,工作效率高,维护费用低,应用十分广泛。其中表贴式永磁同步电机具有体积 小、效率高、转矩密度大等优点,被广泛应用到众多工程领域中,是一类极具发展潜力的节能环保型电机。 获得永磁电机气隙磁场的分布,是计算转矩特性、涡流损耗、噪声和振动等的关键。在气隙磁场的求解方 法中,有限元法适用于多种结构复杂的电机,但速度较慢,在电机的优化设计中不便使用;解析法速度快, 使用价值高,因此,进一步优化解析方法,创新解析模型,减小计算量,对电机进一步优化计算具有重要 意义。
电机磁场解析法是利用数学理论给出特定边值问题的精确解,包括分离变量法、保角变换法和镜像法 等。通常是将电机二维平面划分为有限个子域,如:永磁体、气隙、定子槽。在每个子域内结合相关磁场 边界条件,通过求解偏微分方程,得到傅里叶级数形式的矢量磁位解析解。已有文献多是对电机整体区域 建立数学模型,在每个槽内建立子域,方程维数与电机槽数成正比,各次谐波系数与每个槽内激励有关, 计算复杂。
发明内容
本发明的目的是提供一种重构表贴式永磁电机电枢磁场的新方法。本发明将永磁电机总的电枢反应磁 场表示为单槽通电时的气隙磁场,对单个槽磁场进行子域法求解,降低矩阵的维数,降低永磁电机磁场分 析的复杂性并能够重构出与有限元计算相吻合的电枢反应磁场。为达到上述目的,本发明所采用的技术方 案如下:
一种基于单槽模型的表贴式永磁电机电枢磁场重构方法,基于如下设定:
a)分析区域在二维平面,不计端部效应。
b)定子和转子铁芯磁导率为无穷大,铁芯磁阻产生的影响忽略不计。
c)永磁体区域的相对磁导率与空气区域相同。
d)定子槽为径向开口槽,槽内电流密度均匀分布。
其特征在于,包括下列步骤:
第一步.建立永磁电机单槽模型
在永磁电机横截面上建立平面极坐标系,极坐标系中的动点坐标为(r,θ),其中r为该点到极点的距离, θ为从极轴到该点与极点连线所成的角度,将单槽模型求解模型划分为定子槽、气隙2个区域。电机的定 子槽数为Ns,极轴θ=0处位于一个槽的中心,将其编号为1号槽,逆时针一次编号为1号槽,2号槽,......, Ns号槽,编号为i的槽对应的位置为θi,且θi=2π(i-1)/Ns;选定1号槽内通有电枢电流,2~Ns号槽不通电, 此模型为单槽模型。
第二步.单层绕组磁场重构
1号槽的定子槽内存在电流密度激励,该区域中的泊松方程:
式中,A11(r,θ)为矢量磁位,J1为槽内电流密度,μ0为空气磁导率;
气隙子域内的拉普拉斯方程:
式中,A2(r,θ)为矢量磁位;
单槽模型的气隙磁场分布关于θ=0对称,电枢总磁场是各定子槽内电流产生的单个磁场之和;因此, 气隙的磁通密度的径向分量Br和切向分量Bt为:
式中,Ii为第i号槽内的电流值,B1r *、B1t *为1号槽通单位电流时的磁感应强度表达式,n为谐波阶次, En和Fn为第n次谐波系数;
第三步.双层绕组磁场重构
(1)分布绕组
在分布绕组中,槽内两相电流分别处于槽顶部和底部,将槽中的两部分电流取代数和的形式,等效电 流仍处在槽中心位置;将第i个槽内的两个绕组电流取为等效净电流Ii,即:
Ii=Ii1+Ii2
式中Ii1、Ii2分别为槽内上层电流和下层电流;
(2)集中绕组
将电流Ii平均分解为Ii1和Ii2,即
Ii1=Ii2=0.5Ii
此时两个电流单独产生的磁场将关于槽中线对称;因此,在集中绕组槽内电流Ii1、Ii2为任意值时,可 等效为Ii1、Ii2两部分独立工作;当θi=0时,满足
B1r=B1r1+B1r2
B1r1(r,θ)=-B1r2(r,-θ)
将各个槽内电流产生的的磁场累加即可获得总的磁场分布;径向磁感应强度为:
同理,切向磁感应强度表达式为:
其中,γ表示电流Ii1、Ii2关于槽中线逆时针或顺时针偏移的角度。
本发明联立各个子域内磁场表达式及边界条件方程,减少运算矩阵维数和计算复杂度,传统电机磁 场解析求解过程中,联立的方程维数和电机槽数成正比,本发明采用的基于单槽模型的磁场重构法所列写 的矩阵的维数只与单个槽内的谐波阶数相关,降低了计算复杂度。此外,经过仿真可知,采用本发明的方 法进行的电枢磁场重构与有限元方法具有很好的吻合度。
附图说明
图1:表贴式永磁电机横截面结构示意图
图2:单槽通电下的磁力线分布
图3:单槽通电下的径向磁通分布
图4(a)为单槽电机模型示意图,图4(b)揭示了定子槽、气隙、相邻槽之间的分界面边界条件设置
图5(a)和(b)分别为分布式双层绕组电流示意图及等效净电流示意图
图6:双层绕组磁场分布
图7:集中绕组电流示意图
图8:表贴式永磁电机10极12槽单层绕组分布示意图
图9:表贴式永磁电机10极12槽双层绕组分布示意图
图10:单层绕组电机径向磁感应强度
图11:单层绕组电机切向磁感应强度
图12:单层绕组电机径向磁感应强度谐波分量
图13:单层绕组电机切向磁感应强度谐波分量
图14:双层绕组电机径向磁感应强度
图15:双层绕组电机切向磁感应强度
图16:双层绕组电机径向磁感应强度谐波分量
图17:双层绕组电机切向磁感应强度谐波分量
具体实施方式
本发明提供了一种求解永磁电机电枢磁场的方法,下面结合附图和实施例对本发明进行详细的描述。 本发明具体实施步骤如下:
1.单层绕组
(1)电机横截面如图1所示,在电机横截面上建立平面极坐标系。极轴(图中θ=0处)位于一个槽的中 心。定子槽数为Ns,第i个槽对应的位置为
仅1号槽通电时,激励电流为1A,磁场分布如图2、3所示。不考虑饱和效应和端部效应,第i个槽 电流Ii产生的磁场分布将旋转i个槽距角,且磁感应强度与绕组电流成正比。根据图4(b),单槽通电时的 磁场分布关于θ=0对称,考虑到电枢总磁场是各定子槽内电流产生的单个磁场之和,即,定子绕组磁通密 度的法向分量和切向分量可表示为:
由静磁场的唯一性原理(uniqueness theory),并考虑到电机结构的对称性,电机模型可简化为如图4(a) 所示的单槽模型,其定子槽、气隙、相邻槽之间的分界面边界条件设置如图4(b)所示。
(2)将单槽求解模型划分为定子槽、气隙2个区域。设1号槽内电流密度为J1,且该区域中矢量磁位为 A11(r,θ)。由于定子槽内存在电流密度激励,在定子槽区域中列写泊松方程:
其边界条件可表示为:
由分离变量法可得满足边界条件的定子槽子域通解:
式中:
K0、Am为定子槽子域的谐波系数;m为定子槽子域磁场谐波阶数。
(3)气隙子域的拉普拉斯方程如下:
气隙子域的边界条件为:
由于气隙磁场以槽中心为轴对称分布,因此附加边界条件为:
由以上边界条件得出,气隙子域的通解如下:
式中:Cn为气隙子域的谐波系数;n为气隙子域磁场谐波阶数。
(4)在两类子域通解已知的情况下,根据定子槽与气隙子域交界面边界条件得出谐波系数。即r=R2处, 根据法向磁密相等,可得:
A11(R2,θ)=A2(R2,θ)
结合式定子槽通解和气隙子域通解,在-β/2≤θ≤β/2范围内,可得定子槽子域中谐波系数如下:
为便于谐波系数的求解,定义函数:
则各次谐波可表示为:
在气隙与定子槽交界面处有:
即:
将定子槽子域、气隙子域通解中的傅里叶级数取为有限次,联立上式并改写为矩阵形式,并进行线性 方程组求解,即:
其中,B1为1×1的矩阵,B2为1×n矩阵,B3为m×m矩阵,B4为n×m矩阵,B5为m×n矩阵,B6为n×n 矩阵,K0为1×1的矩阵,Am为m×1矩阵,Cn为n×1矩阵。L0为常数,Ln为n×1矩阵。通过求解上述方 程得到每个子域未知的各阶谐波系数K0、Am、Cn。
2.双层绕组
(1)分布式绕组
在分布绕组中,槽内两相电流分别处于槽顶部和底部,如图5(a)所示,将槽中的两部分电流取代数和 的形式,等效电流仍处在槽中心位置,如图5(b)所示。
将第i个槽内的两个绕组电流取为等效净电流Ii,即:
Ii=Ii1+Ii2
式中Ii1、Ii2分别为槽内上层电流和下层电流。将等效电流代入单层绕组磁场表达式中,即可得到所有 槽通电情况下的气隙磁感应强度。
(2)集中绕组
将电流Ii平均分解为Ii1和Ii2,即
Ii1=Ii2=0.5Ii
此时两个电流单独产生的磁场将关于槽中线对称。如图6所示,当θi=0时,满足:
B1r=B1r1+B1r2
B1r1(r,θ)=-B1r2(r,-θ)
考虑到B1r为n阶正弦函数之和,则B1r1和B1r2表示为:
其中,γ表示电流Ii1、Ii2关于槽中线逆时针或顺时针偏移的角度,Gn为谐波系数,可表示为:
在集中绕组槽内电流Ii1、Ii2为任意值时,可等效为Ii1、Ii2两部分独立工作,如图7所示。则原槽内电 流Ii产生的磁场可由Ii1和Ii2各自产生的磁场叠加而成。进一步,将各个槽内电流产生的的磁场累加即可 获得总的磁场分布。
径向磁感应强度为:
同理,切向磁感应强度表达式为:
其中:
本发明在单层绕组电机且单槽通电的解析模型中施加附加边界条件,利用子域法及叠加原理,获取 电机的电枢反应磁场分布。对于双层绕组结构,同一槽内两绕组通以相同幅值的电流,推导出每个绕组产 生的磁场,进而重构出电机电枢反应的总磁场,得到磁场表达式。已有文献多是对电机整体区域建立数学 模型,在每个槽内建立子域,若与本发明划分子域相同,则其方程维数至少为(Ns+m×Ns+2×n),本发明采 用的方法可将求解矩阵维数可减少至(1+m+n)。
3.仿真实验
分别以一台10极12槽单层绕组和一台10极12槽双层绕组表贴式永磁同步电机为例,其绕组结构分 布如图8、9所示,在双层绕组电机中,γ=β/4。分析其电枢反应磁场,设气隙磁感应强度、槽磁感应强度 的谐波阶数均为18。三相电流瞬态情况分别为:Ia=5A,Ib=Ic=-2.5A。
图10-图17为采用本发明得到的气隙内半径为29mm处电枢反应磁感应强度计算结果及分相应谐波分 量法所得的结果与有限元法所得结果的对比。由图可以看出,本发明的结果与有限元计算的结果能够很好 地吻合。
Claims (1)
1.一种基于单槽模型的表贴式永磁电机电枢磁场重构方法,基于如下设定:
a)分析区域在二维平面,不计端部效应;
b)定子和转子铁芯磁导率为无穷大,铁芯磁阻产生的影响忽略不计;
c)永磁体区域的相对磁导率与空气区域相同;
d)定子槽为径向开口槽,槽内电流密度均匀分布。
其特征在于,包括下列步骤:
第一步.建立永磁电机单槽模型
在永磁电机横截面上建立平面极坐标系,极坐标系中的动点坐标为(r,θ),其中r为该点到极点的距离,θ为从极轴到该点与极点连线所成的角度,将单槽模型求解模型划分为定子槽、气隙2个区域。电机的定子槽数为Ns,极轴θ=0处位于一个槽的中心,将其编号为1号槽,逆时针一次编号为1号槽,2号槽,......,Ns号槽,编号为i的槽对应的位置为θi,且θi=2π(i-1)/Ns;选定1号槽内通有电枢电流,2~Ns号槽不通电,此模型为单槽模型;
第二步.单层绕组磁场重构
1号槽的定子槽内存在电流密度激励,该区域中的泊松方程:
式中,A11(r,θ)为矢量磁位,J1为槽内电流密度,μ0为空气磁导率;
气隙子域内的拉普拉斯方程:
式中,A2(r,θ)为矢量磁位;
单槽模型的气隙磁场分布关于θ=0对称,电枢总磁场是各定子槽内电流产生的单个磁场之和;因此,气隙的磁通密度的径向分量Br和切向分量Bt为:
式中,Ii为第i号槽内的电流值,B1r *、B1t *为1号槽通单位电流时的磁感应强度表达式,n为谐波阶次,En和Fn为第n次谐波系数;
第三步.双层绕组磁场重构
(1)分布绕组
在分布绕组中,槽内两相电流分别处于槽顶部和底部,将槽中的两部分电流取代数和的形式,等效电流仍处在槽中心位置;将第i个槽内的两个绕组电流取为等效净电流Ii,即:
Ii=Ii1+Ii2
式中Ii1、Ii2分别为槽内上层电流和下层电流;
(2)集中绕组
将电流Ii平均分解为Ii1和Ii2,即
Ii1=Ii2=0.5Ii
此时两个电流单独产生的磁场将关于槽中线对称;因此,在集中绕组槽内电流Ii1、Ii2为任意值时,可等效为Ii1、Ii2两部分独立工作;当θi=0时,满足
B1r=B1r1+B1r2
B1r1(r,θ)=-B1r2(r,-θ)
将各个槽内电流产生的的磁场累加即可获得总的磁场分布;径向磁感应强度为:
同理,切向磁感应强度表达式为:
其中,γ表示电流Ii1、Ii2关于槽中线逆时针或顺时针偏移的角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910082678.1A CN109768685B (zh) | 2019-01-28 | 2019-01-28 | 基于单槽模型的表贴式永磁电机电枢磁场重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910082678.1A CN109768685B (zh) | 2019-01-28 | 2019-01-28 | 基于单槽模型的表贴式永磁电机电枢磁场重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109768685A true CN109768685A (zh) | 2019-05-17 |
CN109768685B CN109768685B (zh) | 2021-03-02 |
Family
ID=66455480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910082678.1A Expired - Fee Related CN109768685B (zh) | 2019-01-28 | 2019-01-28 | 基于单槽模型的表贴式永磁电机电枢磁场重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109768685B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110661468A (zh) * | 2019-09-29 | 2020-01-07 | 福州大学 | 一种无轴承磁通切换电机的电磁场性能获取方法 |
CN112100814A (zh) * | 2020-08-13 | 2020-12-18 | 河海大学 | 一种分数槽四层绕组永磁电机永磁体涡流损耗计算方法 |
CN113076665A (zh) * | 2021-05-14 | 2021-07-06 | 天津大学 | 一种永磁电机空载气隙磁场分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005692A (zh) * | 2015-07-08 | 2015-10-28 | 三峡大学 | 一种基于解析法的永磁电机磁场分析与转矩计算方法 |
CN108319768A (zh) * | 2018-01-23 | 2018-07-24 | 湖北西浦电机科技有限责任公司 | 一种基于计算机的永磁电机电枢反应磁场预测方法 |
-
2019
- 2019-01-28 CN CN201910082678.1A patent/CN109768685B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005692A (zh) * | 2015-07-08 | 2015-10-28 | 三峡大学 | 一种基于解析法的永磁电机磁场分析与转矩计算方法 |
CN108319768A (zh) * | 2018-01-23 | 2018-07-24 | 湖北西浦电机科技有限责任公司 | 一种基于计算机的永磁电机电枢反应磁场预测方法 |
Non-Patent Citations (1)
Title |
---|
董川: ""表贴式永磁电机子域模型及涡流损耗分析"", 《中国优秀硕士学位论文全文数据库 工程科技2辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110661468A (zh) * | 2019-09-29 | 2020-01-07 | 福州大学 | 一种无轴承磁通切换电机的电磁场性能获取方法 |
CN110661468B (zh) * | 2019-09-29 | 2021-04-27 | 福州大学 | 一种无轴承磁通切换电机的电磁场性能获取方法 |
CN112100814A (zh) * | 2020-08-13 | 2020-12-18 | 河海大学 | 一种分数槽四层绕组永磁电机永磁体涡流损耗计算方法 |
CN112100814B (zh) * | 2020-08-13 | 2022-09-16 | 河海大学 | 一种分数槽四层绕组永磁电机永磁体涡流损耗计算方法 |
CN113076665A (zh) * | 2021-05-14 | 2021-07-06 | 天津大学 | 一种永磁电机空载气隙磁场分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109768685B (zh) | 2021-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tiegna et al. | Overview of analytical models of permanent magnet electrical machines for analysis and design purposes | |
Roberts | A study of brushless doubly-fed (induction) machines | |
Wang et al. | Rotor eddy-current loss in permanent-magnet brushless AC machines | |
Zarko et al. | Analytical calculation of magnetic field distribution in the slotted air gap of a surface permanent-magnet motor using complex relative air-gap permeance | |
Dubas et al. | Analytical solution of the magnetic field in permanent-magnet motors taking into account slotting effect: No-load vector potential and flux density calculation | |
Xu et al. | Transient model of a doubly excited reluctance motor | |
El-Refaie et al. | Analysis of surface permanent magnet machines with fractional-slot concentrated windings | |
Holm et al. | Analytical modeling of a permanent-magnet synchronous machine in a flywheel | |
Devillers et al. | A review of subdomain modeling techniques in electrical machines: Performances and applications | |
Ionel et al. | Assessment of torque components in brushless permanent-magnet machines through numerical analysis of the electromagnetic field | |
Brahim et al. | Cogging torque minimization of surface-mounted permanent magnet synchronous machines using hybrid magnet shapes | |
CN109768685A (zh) | 基于单槽模型的表贴式永磁电机电枢磁场重构方法 | |
de La Barriere et al. | An analytical model for the computation of no-load eddy-current losses in the rotor of a permanent magnet synchronous machine | |
Zhao et al. | A novel quasi-3D analytical model for axial flux motors considering magnetic saturation | |
Dalal et al. | Analytical model for permanent magnet motor with slotting effect, armature reaction, and ferromagnetic material property | |
Parthasarathy | Linear PM generator for wave energy conversion | |
Guo et al. | Cogging force optimization of double-sided tubular linear machine with tooth-cutting | |
Echle et al. | Design and comparison of radial flux and axial flux brushless dc motors for power tool applications | |
Røkke et al. | Analytical calculation of yoke flux patterns in fractional-slot permanent magnet machines | |
Hu et al. | Analytical modeling of magnetic field considering the saturation in switched reluctance motor | |
Ali et al. | Cogging torque Mitigation for PMSM using stator slots design and Magnets skewing | |
Liu | Critical wavelength of eddy currents and its influence on harmonic losses in solid poles | |
Deng et al. | Comprehensive salient-pole synchronous machine parametric design analysis using time-step finite element-state space modeling techniques | |
Tikellaline et al. | Magnetic field analysis of double excited synchronous motor using numerical conformal mapping | |
CN104362917A (zh) | 车用交流发电机漏磁问题的优化设计方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210302 Termination date: 20220128 |
|
CF01 | Termination of patent right due to non-payment of annual fee |