CN110690770A - 定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法 - Google Patents

定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法 Download PDF

Info

Publication number
CN110690770A
CN110690770A CN201911066869.5A CN201911066869A CN110690770A CN 110690770 A CN110690770 A CN 110690770A CN 201911066869 A CN201911066869 A CN 201911066869A CN 110690770 A CN110690770 A CN 110690770A
Authority
CN
China
Prior art keywords
air gap
magnetic
permanent magnet
flux density
magnet motor
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
Application number
CN201911066869.5A
Other languages
English (en)
Other versions
CN110690770B (zh
Inventor
王明杰
武洁
郭磊磊
李从善
张志艳
杨存祥
邱洪波
申永鹏
刘小梅
李彦彦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhengzhou University of Light Industry
Original Assignee
Zhengzhou University of Light Industry
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhengzhou University of Light Industry filed Critical Zhengzhou University of Light Industry
Priority to CN201911066869.5A priority Critical patent/CN110690770B/zh
Publication of CN110690770A publication Critical patent/CN110690770A/zh
Application granted granted Critical
Publication of CN110690770B publication Critical patent/CN110690770B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K1/00Details of the magnetic circuit
    • H02K1/06Details of the magnetic circuit characterised by the shape, form or construction
    • H02K1/12Stationary parts of the magnetic circuit
    • H02K1/16Stator cores with slots for windings
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K21/00Synchronous motors having permanent magnets; Synchronous generators having permanent magnets
    • H02K21/02Details
    • H02K21/021Means for mechanical adjustment of the excitation flux
    • H02K21/022Means for mechanical adjustment of the excitation flux by modifying the relative position between field and armature, e.g. between rotor and stator
    • H02K21/023Means for mechanical adjustment of the excitation flux by modifying the relative position between field and armature, e.g. between rotor and stator by varying the amount of superposition, i.e. the overlap, of field and armature
    • H02K21/024Radial air gap machines
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K21/00Synchronous motors having permanent magnets; Synchronous generators having permanent magnets
    • H02K21/12Synchronous motors having permanent magnets; Synchronous generators having permanent magnets with stationary armatures and rotating magnets
    • H02K21/14Synchronous motors having permanent magnets; Synchronous generators having permanent magnets with stationary armatures and rotating magnets with magnets rotating within the armatures
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K2201/00Specific aspects not provided for in the other groups of this subclass relating to the magnetic circuits
    • H02K2201/03Machines characterised by aspects of the air-gap between rotor and stator

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Permanent Field Magnets Of Synchronous Machinery (AREA)
  • Iron Core Of Rotating Electric Machines (AREA)

Abstract

本发明提出了一种定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法,首先考虑定子开槽效应,建立复数气隙相对磁导与气隙磁密之间的关系;然后,以电机一个齿距为求解模型,在极坐标系下建立气隙和磁极区域的拉普拉斯方程,得到气隙和磁极区域的磁场通解;再利用槽口表面与铁心齿尖边角区附近磁场发散的数学规律,在槽口区域设置边界条件,求解出复数气隙相对磁导的径向分量和切向分量;最后,计算定子开槽时内转子、外转子永磁电机的有槽气隙磁密,有限元验证模型的准确性。本发明解决了齿槽效应对气隙磁场径向和切向方向上的影响,准确计算不同极槽配合下齿槽效应对永磁电机的径向和切向气隙磁场分布,具有较好的简捷性和通用性。

Description

定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法
技术领域
本发明涉及永磁电机技术领域,特别是指一种用于定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法。
背景技术
永磁电机具有高功率密度、高效率、高可靠性等优点,在电动汽车、风力发电、伺服驱动等领域应用广泛,引起了业界学者极大的研究兴趣。磁场分析是电机性能分析的基础,为了提高磁场解析的准确度必须考虑齿槽效应对气隙磁场的影响,考虑电机开槽对气隙磁场的影响一直是研究的热点问题,其结果的准确性直接影响电机的电磁性能。
考虑齿槽效应时,目前常采用的解析方法有:气隙相对磁导模型、复数气隙相对磁导模型、精确子域模型、磁网络模型。气隙相对磁导模型因求解方法简便,是目前应用最广的模型,但该模型忽略了齿槽效应对切向磁密造成的影响,未能准确计算切向磁密。目前的复数气隙相对磁导模型基于保角变换法,能准确考虑齿槽开槽影响的径向磁密和切向磁密,但计算方法类似数值解法,求解方法比较复杂且解析式不直观。精确子域模型在所有解析法中精度最高,但边界条件矩阵方程较多,求解方法复杂且与结构参数关系不直接,计算工作量较大。磁网络模型能考虑磁路局部饱和的影响,但求解方法依靠磁力线走向画磁网络,电机结构变化时模型缺乏通用性。
以上分析模型中各有各的优缺点和适用范围,其中能以结构参数直观表示,并以精简的求解方法反应齿槽效应的,仅有气隙相对磁导模型,但该模型存在缺点时无法计及切向方向齿槽效应,而其余模型存在的共同问题是计算相对复杂,无统一的结构参数关系来反应齿槽效应的影响。如以简单的通用结构参数给出并准确的反映齿槽效应,优化解析方法,创新解析模型,则能方便解决复杂气隙磁场的分析。
发明内容
针对上述背景技术中存在的不足,本发明提出了一种用于定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法,解决了现有齿槽效应对气隙磁场径向和切向方向上的影响以及计算复杂度高的技术问题。
本发明的技术方案是这样实现的:
一种定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法,其步骤如下:
S1、根据定子开槽效应,建立复数气隙相对磁导与气隙磁密之间的关系;
S2、以电机一个齿距为求解模型,在极坐标系下建立气隙和磁极区域的拉普拉斯方程,得到气隙和磁极区域的磁场通解;
S3、利用槽口表面与铁心齿尖边角区附近磁场发散的数学规律,在槽口区域设置边界条件,并结合卡氏系数确定出磁场系数,求解出复数气隙相对磁导的径向分量和切向分量;
S4、根据求解出的复数气隙相对磁导,计算出不同极槽配合及径向和平行充磁方式下,定子开槽时内转子永磁电机的有槽气隙磁密和外转子永磁电机的有槽气隙磁密。
所述气隙磁密包括有槽气隙磁密和无槽气隙磁密,且有槽气隙磁密等于无槽气隙磁密与复数气隙相对磁导的共轭复数之积,即:
Bs(r,θ)=Bnos(r,θ)·λ*(r,θ) (1),
λ(r,θ)=λa(r,θ)+jλb(r,θ) (2),
其中,Bs(r,θ)为有槽气隙磁密,Bnos(r,θ)为无槽气隙磁密,λ(r,θ)为复数气隙相对磁导,λ*(r,θ)为复数气隙相对磁导的共轭复数,j为虚数单位,λa(r,θ)为复数气隙相对磁导的实部,λb(r,θ)为复数气隙相对磁导的虚部,r为气隙半径,θ为机械角;
所述有槽气隙磁密的径向分量Br_s(R0,θ)和切向分量Bθ_s(R0,θ)分别为:
Br_s(r,θ)=Br_nos(r,θ)·λa(r,θ)+Bθ_nos(r,θ)·λb(r,θ) (3),
Bθ_s(r,θ)=Bθ_nos(r,θ)·λa(r,θ)-Br_nos(r,θ)·λb(r,θ) (4),
其中,Br_nos(r,θ)为无槽气隙磁密的径向分量,Bθ_nos(r,θ)为无槽气隙磁密的切向分量;
所述无槽气隙磁密的切向分量小,令切向分量Bθ_nos(r,θ)=0,则:
Br_s(r,θ)≈Br_nos(r,θ)·λa(r,θ) (5),
Bθ_s(r,θ)≈-Br_nos(r,θ)·λb(r,θ) (6);
令气隙半径r=R0为气隙处位置,Rr为转子铁心半径,Rs为定子铁心半径,则复数气隙相对磁导的实部λa(r,θ)和虚部λb(r,θ)分别为:
Figure BDA0002259647930000022
Figure BDA0002259647930000023
所述气隙和磁极区域的拉普拉斯方程为:
Figure BDA0002259647930000024
其中,Rr为转子铁心半径,Rs为定子铁心半径,θs为齿距角,A为气隙和磁极区域的矢量磁位通解,则
Figure BDA0002259647930000031
其中,A0、k为待定系数,An、Bn为磁矢位的n次谐波系数,z为槽数;
则极坐标系下气隙和磁极区域的磁密通解分别为:
Figure BDA0002259647930000033
其中,Br(r,θ)为径向磁密通解,Bθ(r,θ)为切向磁密通解。
所述根据槽口表面与铁心齿尖边角区附近磁场发散的数学规律,一个齿槽表面的切向磁密Bcc'(θ)为:
Figure BDA0002259647930000034
其中,c0为待定常数,sign为符号函数,αs0为槽开口角度;
对切向磁密Bcc'(θ)进行傅里叶级数展开:
Figure BDA0002259647930000035
其中,
Figure BDA0002259647930000036
根据齿槽与气隙区域边界之间的关系,在一个齿距区域内建立槽口区域的边界条件方程:
S31.1、在齿中心线处,磁力线垂直穿过气隙和磁极区域:
Figure BDA0002259647930000037
S31.2、在转子表面:
Figure BDA0002259647930000038
S31.3、在定子齿槽表面:
Bcc'(θ)=Bθ(Rs,θ) (17);
根据步骤S31.1至步骤S31.3的边界条件得到谐波磁密系数,进而得到气隙和磁极区域的磁密通解为:
Figure BDA0002259647930000041
Figure BDA0002259647930000042
电机开槽时一个齿距角θs内的平均气隙磁密为:
Figure BDA0002259647930000043
开槽时气隙磁密在齿中心线处最大,最大气隙磁密为:
Figure BDA0002259647930000044
根据卡氏系数得:
其中,
Figure BDA0002259647930000047
bs为齿宽,ts为齿距。
所述复数气隙相对磁导包括内转子永磁电机的复数气隙相对磁导和外转子永磁电机的复数气隙相对磁导,所述内转子永磁电机的复数气隙相对磁导的求解方法为:
S32.1、无槽气隙磁密分布比较均匀,齿中心线处气隙磁密最大,电机无槽时的平均气隙磁密为:
Figure BDA0002259647930000048
S32.2、内转子永磁电机的复数气隙相对磁导的径向分量和切向分量分别为:
Figure BDA0002259647930000049
Figure BDA00022596479300000410
其中,
Figure BDA0002259647930000051
g1=δ+hm为定转子之间的距离,δ为气隙长度,hm为永磁体磁化方向高度;
S32.3、根据步骤S32.1至步骤S32.2求得外转子永磁电机的复数气隙相对磁导为:
Figure BDA0002259647930000052
Figure BDA0002259647930000053
所述定子开槽时内、外转子结构永磁电机的径向和切向气隙磁密的求解方法为:
S41、由无槽磁场和复数气隙相对磁导,求出齿槽效应时的空载磁场分布,又永磁电机的无槽气隙磁密为:
当np≠1时:
Figure BDA0002259647930000054
其中,
Figure BDA0002259647930000055
KB(n)·fBr(r)为无槽径向磁密谐波幅值,KB(n)·f(r)为无槽切向磁密谐波幅值,
Figure BDA0002259647930000056
Mn=Mrn+npMθn;径向充磁方式下:
Figure BDA0002259647930000057
A3n=np;平行充磁方式下:
Figure BDA0002259647930000061
其中,p为极对数,αp为永磁体极弧系数,μr为永磁体相对磁导率,Br为永磁体剩磁;
S42、考虑齿槽效应时,根据步骤S32.2中内转子永磁电机的复数气隙相对磁导,则内转子永磁电机的有槽气隙磁密为:
Figure BDA0002259647930000062
S43、考虑齿槽效应时,根据步骤S32.3中外转子永磁电机的复数气隙相对磁导,则外转子永磁电机的有槽气隙磁密为:
Figure BDA0002259647930000063
本技术方案能产生的有益效果:本发明基于求解区域的控制方程和合理的边界条件,给出复数气隙相对磁导函数的径向分量和切向分量,能准确反应齿槽效应对气隙相对磁导的径向和切向分布影响,准确计算定子开槽永磁电机的径向和切向气隙磁场,简化了繁琐的计算过程,提高计算效率,对进一步分析电机空载反电动势、齿槽转矩、不平衡磁拉力等电机性能及电机优化设计具有重要意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的有槽永磁电机结构横截面示意图;
图2为本发明的无槽永磁电机结构横截面示意图;
图3为本发明的单一齿距区域的解析模型;
图4为本发明的单一齿槽表面的切向磁密分布;
图5为本发明的9槽结构半闭口槽的复数气隙相对磁导分布;
图6为本发明的9槽结构矩形槽的复数气隙相对磁导分布;
图7为本发明的24槽结构半闭口槽的复数气隙相对磁导分布;
图8为本发明的24槽结构矩形槽的复数气隙相对磁导分布;
图9为本发明的9槽8极半闭口槽气隙磁密分布;
图10为本发明的9槽8极矩形槽气隙磁密分布;
图11为本发明的24槽20极半闭口槽气隙磁密分布;
图12为本发明的24槽20极矩形槽气隙磁密分布。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明以内转子永磁电机为实施例提供了一种用于定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法,有槽和无槽的永磁电机结构分别如图1和图2所示,电机结构参数如下:气隙长度为δ,永磁体磁化方向高度为hm,定转子之间距离g1=δ+hm,Rs为定子铁心半径,Rr为转子铁心半径;内转子的定子铁心半径Rs=Rr+g1,外转子的定子铁心半径Rs=Rr-g1,齿距角为θs,槽开口角度为αs0
Figure BDA0002259647930000071
为气隙处位置。
以电机一个齿距为求解区域,对解析模型作如下基本假设:
1)定转子铁心磁导率无穷大;
2)忽略磁路饱和和端部效应的影响;
3)忽略临近齿槽效应,定子槽为无限深。
S1、考虑定子开槽效应时,建立复数气隙相对磁导与气隙磁密之间的关系;所述气隙磁密包括有槽气隙磁密和无槽气隙磁密,且有槽气隙磁密等于无槽气隙磁密与复数气隙相对磁导的共轭复数之积,即:
Bs(r,θ)=Bnos(r,θ)·λ*(r,θ) (1),
λ(r,θ)=λa(r,θ)+jλb(r,θ) (2),
其中,Bs(r,θ)为有槽气隙磁密,Bnos(r,θ)为无槽气隙磁密,λ(r,θ)为复数气隙相对磁导,λ*(r,θ)为复数气隙相对磁导的共轭复数,j为虚数单位,λa(r,θ)为复数气隙相对磁导的实部,用来考虑径向齿槽效应的影响,λb(r,θ)为复数气隙相对磁导的虚部,用来考虑切向齿槽效应的影响,r为气隙半径,θ为机械角。
根据公式(1)和公式(2)可知,有槽气隙磁密的径向分量Br_s(R0,θ)和切向分量Bθ_s(R0,θ)分别为:
Br_s(r,θ)=Br_nos(r,θ)·λa(r,θ)+Bθ_nos(r,θ)·λb(r,θ) (3),
Bθ_s(r,θ)=Bθ_nos(r,θ)·λa(r,θ)-Br_nos(r,θ)·λb(r,θ) (4),
其中,Br_nos(r,θ)为无槽气隙磁密的径向分量,Bθ_nos(r,θ)为无槽气隙磁密的切向分量。
所述无槽气隙磁密的切向分量小,主要为径向分量,为方便计算复数气隙相对磁导,令切向分量Bθ_nos(r,θ)=0,则:
Br_s(r,θ)≈Br_nos(r,θ)·λa(r,θ) (5),
Bθ_s(r,θ)≈-Br_nos(r,θ)·λb(r,θ) (6)。
令气隙半径r=R0
Figure BDA0002259647930000081
为气隙处位置,Rr为转子铁心半径,Rs为定子铁心半径,则复数气隙相对磁导的实部λa(r,θ)和虚部λb(r,θ)分别为:
Figure BDA0002259647930000082
Figure BDA0002259647930000083
由此知,需要先求出有槽和无槽时的气隙磁密,才能求出复数气隙相对磁导的分布。
S2、以电机一个齿距为求解模型,在极坐标系下建立气隙和磁极区域的拉普拉斯方程,得到气隙和磁极区域的磁场通解。为方便计算复数气隙相对磁导,只考虑电枢绕组单独作用,且电枢磁势幅值位于槽中心处,单一齿距内的求解区域如图3所示。
采用矢量磁位建立求解区域的控制方程,电枢绕组单独作用,在气隙和磁极区域为拉普拉斯方程:
Figure BDA0002259647930000084
其中,Rr为转子铁心半径,Rs为定子铁心半径,θs为齿距角,A为气隙和磁极区域的矢量磁位通解,则
Figure BDA0002259647930000085
其中,A0、k为待定系数,An、Bn为磁矢位的n次谐波系数,z为槽数;
则极坐标系下气隙和磁极区域的磁位通解分别为:
Figure BDA0002259647930000091
其中,Br(r,θ)为径向磁密通解,Bθ(r,θ)为切向磁密通解。
S3、利用槽口表面与铁心齿尖边角区附近磁场发散的数学规律,在槽口区域设置合适的边界条件,结合卡氏系数确定出磁场系数,求解出复数气隙相对磁导的径向分量和切向分量。
如图3所示,在齿槽表面BB'处,齿表面BC和B'C'切向磁密为零,槽口表面CC'与铁心齿尖边角区附近磁场发散,根据槽口表面与铁心齿尖边角区附近磁场发散的数学规律,如公式(13)所示的一个齿槽表面的切向磁密Bcc'(θ),图4为单一齿槽表面BB'的切向磁密分布Bcc'(θ):
其中,c0为未知常数,sign为符号函数,αs0为槽开口角度;
对切向磁密Bcc'(θ)进行傅里叶级数展开:
Figure BDA0002259647930000093
其中,
Figure BDA0002259647930000094
根据齿槽与气隙区域边界之间的关系,在一个齿距区域内建立边界条件系数方程:
S31.1、在齿中心线AB和A'B'处,磁力线垂直穿过气隙和磁极区域,切向磁密为零:
Figure BDA0002259647930000095
S31.2、在转子表面AA',切向磁密为零:
Figure BDA0002259647930000096
S31.3、在定子齿槽表面BB',切向磁密分布为Bcc'(θ):
Bcc'(θ)=Bθ(Rs,θ) (17);
根据步骤S31.1至步骤S31.3的边界条件得到谐波磁密系数,进而得到气隙和磁极区域的磁密通解为:
Figure BDA0002259647930000101
Figure BDA0002259647930000102
电机开槽时一个齿距角θs内的平均气隙磁密为:
Figure BDA0002259647930000103
开槽时气隙磁密在齿中心线处最大,最大气隙磁密为:
Figure BDA0002259647930000104
根据卡氏系数
Figure BDA0002259647930000105
得:
Figure BDA0002259647930000106
其中,
Figure BDA0002259647930000107
bs为齿宽,ts为齿距。
所述复数气隙相对磁导包括内转子永磁电机的复数气隙相对磁导和外转子永磁电机的复数气隙相对磁导,所述内转子永磁电机的复数气隙相对磁导的求解方法为:
S32.1、无槽气隙磁密分布比较均匀,齿中心线处气隙磁密最大,电机无槽时的平均气隙磁密为:
Figure BDA0002259647930000108
S32.2、根据式(7)、(8)、(23)、步骤S31.1至步骤S31.3,得到内转子永磁电机的复数气隙相对磁导的径向分量和切向分量分别为:
Figure BDA0002259647930000109
其中,
Figure BDA0002259647930000111
S32.3、根据步骤S32.1至步骤S32.2求得外转子永磁电机的复数气隙相对磁导为:
Figure BDA0002259647930000112
Figure BDA0002259647930000113
以上所有公式中,气隙处位置R0可取定转子之间的任意半径位置,得出复数气隙相对磁导与槽宽、齿距、及卡氏系数等结构参数有关,计算量小、求解方便。
S4、根据求解出的复数气隙相对磁导,计算出不同极槽配合下,径向和平行充磁方式下,定子开槽时内转子永磁电机的有槽气隙磁密和外转子永磁电机的有槽气隙磁密,再利用有限元验证模型的准确性。
为了验证提出方法的正确性,与传统的气隙相对磁导方法进行对比,需要说明的是传统的气隙相对磁导仅能考虑径向气隙相对磁导,传统的气隙相对磁导计算方法如下:
Figure BDA0002259647930000114
其中,
Figure BDA0002259647930000115
g'为等效气隙长度。
所述定子开槽时内、外转子结构永磁电机的径向和切向气隙磁密的求解方法为:
S41、由无槽磁场和复数气隙相对磁导,求出齿槽效应时的空载磁场分布,又永磁电机的无槽气隙磁密为:
当np≠1时:
Figure BDA0002259647930000117
其中,
Figure BDA0002259647930000121
KB(n)·fBr(r)为无槽径向磁密谐波幅值,KB(n)·f(r)为无槽切向磁密谐波幅值,
Figure BDA0002259647930000122
Mn=Mrn+npMθn;径向充磁方式下:
Figure BDA0002259647930000123
A3n=np;平行充磁方式下:其中,p为极对数,αp为永磁体极弧系数,μr为永磁体相对磁导率,Br为永磁体剩磁。
S42、考虑齿槽效应时,一般情况下无槽气隙切向分量Bθ_nos(r,θ)和复数气隙相对磁导的虚部λb(r,θ)小,对计算有槽气隙的径向分量Br_s(r,θ)没有影响,但是复数气隙相对磁导的实部λa(r,θ)大,对有槽气隙的切向分量Bθ_s(r,θ)影响大,因此Bθ_nos(r,θ)·λa(r,θ)不可忽略,根据步骤S32.2中内转子永磁电机的复数气隙相对磁导,则内转子永磁电机的有槽气隙磁密为:
Figure BDA0002259647930000125
S43、考虑齿槽效应时,根据步骤S32.3中外转子永磁电机的复数气隙相对磁导,则外转子永磁电机的有槽气隙磁密为:
Figure BDA0002259647930000126
复数气隙相对磁导及有槽气隙磁密的有限元验证:
为验证所提出方法的正确性,分别对一台内转子径向充磁9槽8极、外转子平行充磁24槽20极表贴式永磁电机进行计算,并且每台电机给出两种槽口宽度,半闭口槽的槽口宽小和矩形槽的槽口宽大,其主要参数如表1所示,并用有限元法进行验证。
表1永磁电机结构参数
当磁极中心线对齐槽中心线,R0可以在气隙和磁极区域处任意位置选取,此处选取
Figure BDA0002259647930000132
图5、图6、图7和图8分别给出了本发明方法和有限元法得到的9槽、24槽不同槽宽时一个齿距下的复数气隙相对磁导,同时与传统气隙磁导方法进行对比。有限元法验证时为方便计算,分别在电机有槽绕组区域和无槽槽口表面通入相同的电流密度,使绕组产生的电枢磁势幅值位于槽中心处,分别提取有槽结构和无槽结构气隙磁密,即可得到复数气隙相对磁导分布。可以看到本发明提出的方法拟合较好,证明了本发明提出的求解方法的正确性。
以磁极轴线对齐槽中心线处为坐标原点,利用表1的结构参数,内转子选取R0=Rr+hm+δ/2为气隙中心,外转子选取R0=Rr-hm-δ/2为气隙中心,得到磁密分布如图9、图10、图11和图12所示,可以看出,本发明方法得到的气隙磁密与有限元法结果拟合较好。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法,其特征在于,其步骤如下:
S1、根据定子开槽效应,建立复数气隙相对磁导与气隙磁密之间的关系;
S2、以电机一个齿距为求解模型,在极坐标系下建立气隙和磁极区域的拉普拉斯方程,得到气隙和磁极区域的磁场通解;
S3、利用槽口表面与铁心齿尖边角区附近磁场发散的数学规律,在槽口区域设置边界条件,并结合卡氏系数确定出磁场系数,求解出复数气隙相对磁导的径向分量和切向分量;
S4、根据求解出的复数气隙相对磁导,计算出不同极槽配合及径向和平行充磁方式下,定子开槽时内转子永磁电机的有槽气隙磁密和外转子永磁电机的有槽气隙磁密。
2.根据权利要求1所述的定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法,其特征在于,所述气隙磁密包括有槽气隙磁密和无槽气隙磁密,且有槽气隙磁密等于无槽气隙磁密与复数气隙相对磁导的共轭复数之积,即:
Bs(r,θ)=Bnos(r,θ)·λ*(r,θ)(1),
λ(r,θ)=λa(r,θ)+jλb(r,θ)(2),
其中,Bs(r,θ)为有槽气隙磁密,Bnos(r,θ)为无槽气隙磁密,λ(r,θ)为复数气隙相对磁导,λ*(r,θ)为复数气隙相对磁导的共轭复数,j为虚数单位,λa(r,θ)为复数气隙相对磁导的实部,λb(r,θ)为复数气隙相对磁导的虚部,r为气隙半径,θ为机械角;
所述有槽气隙磁密的径向分量Br_s(R0,θ)和切向分量Bθ_s(R0,θ)分别为:
Br_s(r,θ)=Br_nos(r,θ)·λa(r,θ)+Bθ_nos(r,θ)·λb(r,θ)(3),
Bθ_s(r,θ)=Bθ_nos(r,θ)·λa(r,θ)-Br_nos(r,θ)·λb(r,θ)(4),
其中,Br_nos(r,θ)为无槽气隙磁密的径向分量,Bθ_nos(r,θ)为无槽气隙磁密的切向分量;
所述无槽气隙磁密的切向分量小,令切向分量Bθ_nos(r,θ)=0,则:
Br_s(r,θ)≈Br_nos(r,θ)·λa(r,θ)(5),
Bθ_s(r,θ)≈-Br_nos(r,θ)·λb(r,θ)(6);
令气隙半径r=R0为气隙处位置,Rr为转子铁心半径,Rs为定子铁心半径,则复数气隙相对磁导的实部λa(r,θ)和虚部λb(r,θ)分别为:
Figure FDA0002259647920000012
Figure FDA0002259647920000013
3.根据权利要求1或2所述的定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法,其特征在于,所述气隙和磁极区域的拉普拉斯方程为:
Figure FDA0002259647920000021
其中,Rr为转子铁心半径,Rs为定子铁心半径,θs为齿距角,A为气隙和磁极区域的矢量磁位通解,则
Figure FDA0002259647920000022
其中,A0、k为待定系数,An、Bn为磁矢位的n次谐波系数,z为槽数;
则极坐标系下气隙和磁极区域的磁密通解分别为:
Figure FDA0002259647920000024
其中,Br(r,θ)为径向磁密通解,Bθ(r,θ)为切向磁密通解。
4.根据权利要求3所述的定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法,其特征在于,所述根据槽口表面与铁心齿尖边角区附近磁场发散的数学规律,一个齿槽表面的切向磁密Bcc'(θ)为:
Figure FDA0002259647920000025
其中,c0为待定常数,sign为符号函数,αs0为槽开口角度;
对切向磁密Bcc'(θ)进行傅里叶级数展开:
Figure FDA0002259647920000026
其中,
Figure FDA0002259647920000027
根据齿槽与气隙区域边界之间的关系,在一个齿距区域内建立槽口区域的边界条件方程:
S31.1、在齿中心线处,磁力线垂直穿过气隙和磁极区域:
Figure FDA0002259647920000031
S31.2、在转子表面:
Figure FDA0002259647920000032
S31.3、在定子齿槽表面:
Bcc'(θ)=Bθ(Rs,θ)(17);
根据步骤S31.1至步骤S31.3的边界条件得到谐波磁密系数,进而得到气隙和磁极区域的磁密通解为:
Figure FDA0002259647920000034
电机开槽时一个齿距角θs内的平均气隙磁密为:
Figure FDA0002259647920000035
开槽时气隙磁密在齿中心线处最大,最大气隙磁密为:
Figure FDA0002259647920000036
根据卡氏系数
Figure FDA0002259647920000037
得:
Figure FDA0002259647920000038
其中,
Figure FDA0002259647920000039
bs为齿宽,ts为齿距。
5.根据权利要求1或2或4所述的定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法,其特征在于,所述复数气隙相对磁导包括内转子永磁电机的复数气隙相对磁导和外转子永磁电机的复数气隙相对磁导,所述内转子永磁电机的复数气隙相对磁导的求解方法为:
S32.1、无槽气隙磁密分布比较均匀,齿中心线处气隙磁密最大,电机无槽时的平均气隙磁密为:
Figure FDA0002259647920000041
S32.2、内转子永磁电机的复数气隙相对磁导的径向分量和切向分量分别为:
Figure FDA0002259647920000042
Figure FDA0002259647920000043
其中,
Figure FDA0002259647920000044
g1=δ+hm为定转子之间的距离,δ为气隙长度,hm为永磁体磁化方向高度;
S32.3、根据步骤S32.1至步骤S32.2求得外转子永磁电机的复数气隙相对磁导为:
Figure FDA0002259647920000045
Figure FDA0002259647920000046
6.根据权利要求5所述的定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法,其特征在于,所述定子开槽时内、外转子结构永磁电机的径向和切向气隙磁密的求解方法为:
S41、由无槽磁场和复数气隙相对磁导,求出齿槽效应时的空载磁场分布,又永磁电机的无槽气隙磁密为:
当np≠1时:
Figure FDA0002259647920000047
其中,
Figure FDA0002259647920000051
KB(n)·fBr(r)为无槽径向磁密谐波幅值,KB(n)·f(r)为无槽切向磁密谐波幅值,
Figure FDA0002259647920000052
Mn=Mrn+npMθn;径向充磁方式下:
Figure FDA0002259647920000053
A3n=np;平行充磁方式下:
Figure FDA0002259647920000054
其中,p为极对数,αp为永磁体极弧系数,μr为永磁体相对磁导率,Br为永磁体剩磁;
S42、考虑齿槽效应时,根据步骤S32.2中内转子永磁电机的复数气隙相对磁导,则内转子永磁电机的有槽气隙磁密为:
Figure FDA0002259647920000055
S43、考虑齿槽效应时,根据步骤S32.3中外转子永磁电机的复数气隙相对磁导,则外转子永磁电机的有槽气隙磁密为:
Figure FDA0002259647920000056
CN201911066869.5A 2019-11-04 2019-11-04 定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法 Active CN110690770B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911066869.5A CN110690770B (zh) 2019-11-04 2019-11-04 定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911066869.5A CN110690770B (zh) 2019-11-04 2019-11-04 定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法

Publications (2)

Publication Number Publication Date
CN110690770A true CN110690770A (zh) 2020-01-14
CN110690770B CN110690770B (zh) 2021-04-20

Family

ID=69115694

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911066869.5A Active CN110690770B (zh) 2019-11-04 2019-11-04 定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法

Country Status (1)

Country Link
CN (1) CN110690770B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112016172A (zh) * 2020-09-15 2020-12-01 国网湖南省电力有限公司 一种实心转子感应电机的气隙磁场预测方法及系统
CN112257231A (zh) * 2020-09-22 2021-01-22 重庆交通大学 考虑硅钢片非线性特性的永磁电机解析方法、系统及应用
CN112329293A (zh) * 2020-10-28 2021-02-05 郑州轻工业大学 一种用于永磁直线同步电机空载反电势和推力的计算方法
CN112421919A (zh) * 2020-11-09 2021-02-26 华中科技大学 基于工作磁场谐波定向的永磁电机拓扑构造方法及电机
CN113076665A (zh) * 2021-05-14 2021-07-06 天津大学 一种永磁电机空载气隙磁场分析方法
CN113315424A (zh) * 2021-05-31 2021-08-27 沈阳工业大学 一种v型内置式永磁同步电机齿槽转矩解析方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011130522A (ja) * 2009-12-15 2011-06-30 Minebea Co Ltd 微小回転電気機械の磁気回路
CN103984864A (zh) * 2014-05-16 2014-08-13 天津大学 一种含有不等厚磁极的表贴式永磁电机磁场计算方法
CN105005692A (zh) * 2015-07-08 2015-10-28 三峡大学 一种基于解析法的永磁电机磁场分析与转矩计算方法
CN106682375A (zh) * 2017-03-17 2017-05-17 福州大学 一种考虑转子偏心的无轴承磁通切换电机磁路分析方法
CN110165961A (zh) * 2019-05-21 2019-08-23 东南大学 基于傅里叶拟合的记忆电机调磁电流预测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011130522A (ja) * 2009-12-15 2011-06-30 Minebea Co Ltd 微小回転電気機械の磁気回路
CN103984864A (zh) * 2014-05-16 2014-08-13 天津大学 一种含有不等厚磁极的表贴式永磁电机磁场计算方法
CN105005692A (zh) * 2015-07-08 2015-10-28 三峡大学 一种基于解析法的永磁电机磁场分析与转矩计算方法
CN106682375A (zh) * 2017-03-17 2017-05-17 福州大学 一种考虑转子偏心的无轴承磁通切换电机磁路分析方法
CN110165961A (zh) * 2019-05-21 2019-08-23 东南大学 基于傅里叶拟合的记忆电机调磁电流预测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ZARKO D等: ""Analytical calculation of magnetic field distribution in the slotted air gap of a surface permanent-magnet motor using complex relative air-gap permeance[J]"", 《IEEE TRANSACTIONS ON MAGNETICS》 *
于吉坤: ""定子开槽永磁同步电机气隙比磁导解析计算"", 《电工技术学报》 *
狄冲等: ""基于动态卡式系数的大型潜水电机偏心特征量的计算"", 《电力科学与技术学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112016172A (zh) * 2020-09-15 2020-12-01 国网湖南省电力有限公司 一种实心转子感应电机的气隙磁场预测方法及系统
CN112257231A (zh) * 2020-09-22 2021-01-22 重庆交通大学 考虑硅钢片非线性特性的永磁电机解析方法、系统及应用
CN112257231B (zh) * 2020-09-22 2023-03-21 重庆交通大学 考虑硅钢片非线性特性的永磁电机解析方法、系统及应用
CN112329293A (zh) * 2020-10-28 2021-02-05 郑州轻工业大学 一种用于永磁直线同步电机空载反电势和推力的计算方法
CN112329293B (zh) * 2020-10-28 2024-02-02 郑州轻工业大学 一种用于永磁直线同步电机空载反电势和推力的计算方法
CN112421919A (zh) * 2020-11-09 2021-02-26 华中科技大学 基于工作磁场谐波定向的永磁电机拓扑构造方法及电机
CN113076665A (zh) * 2021-05-14 2021-07-06 天津大学 一种永磁电机空载气隙磁场分析方法
CN113315424A (zh) * 2021-05-31 2021-08-27 沈阳工业大学 一种v型内置式永磁同步电机齿槽转矩解析方法
CN113315424B (zh) * 2021-05-31 2023-05-23 沈阳工业大学 一种v型内置式永磁同步电机齿槽转矩解析方法

Also Published As

Publication number Publication date
CN110690770B (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
CN110690770B (zh) 定子开槽永磁电机复数气隙相对磁导及其磁场的求解方法
Bianchi et al. Design techniques for reducing the cogging torque in surface-mounted PM motors
Aydin et al. Reduction of cogging torque in double-rotor axial-flux permanent-magnet disk motors: A review of cost-effective magnet-skewing techniques with experimental verification
Bianchi et al. An analytical approach to design the PM in PMAREL motors robust toward the demagnetization
Wang et al. Integrated optimization of two design techniques for cogging torque reduction combined with analytical method by a simple gradient descent method
Dutta et al. Winding inductances of an interior permanent magnet (IPM) machine with fractional slot concentrated winding
Xia et al. Cogging torque modeling and analyzing for surface-mounted permanent magnet machines with auxiliary slots
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
CN113507189B (zh) 内置式永磁同步电机设计参数优化方法
CN108206660B (zh) 一种新型五相表嵌式永磁同步电机的转矩解析方法
CN108319768B (zh) 一种基于计算机的永磁电机电枢反应磁场预测方法
Kwon et al. Steady-state modeling and analysis of a double-sided interior permanent-magnet flat linear brushless motor with slot-phase shift and alternate teeth windings
Marwani et al. Design of a coiled rotor synchronous motor dedicated to electric traction
Pina et al. Modeling of synchronous reluctance motors aided by permanent magnets with asymmetric rotor poles
CN117195633B (zh) 一种定转子双永磁容错电机的磁场建模与分析方法
CN110890827B (zh) 永磁直线电机改进的复数气隙相对磁导及磁场的求解方法
Hu et al. Analytical modeling of magnetic field considering the saturation in switched reluctance motor
CN111967170A (zh) 一种开关磁阻电机铁芯损耗计算模型建立方法
Credo et al. Model validation of synchronous motors with a new standstill measurement technique
Jannot et al. An analytical model for interior permanent-magnet synchronous machine with circumferential magnetization design
Liu et al. Design and optimization of permanent magnet synchronous motor based on finite element analysis
Mirahki et al. Analytical prediction of cogging torque for interior permanent magnet synchronous machines
CN112329293A (zh) 一种用于永磁直线同步电机空载反电势和推力的计算方法
CN106897507A (zh) 一种基于电感解析模型的内置式永磁同步电机设计方法
Chen et al. Cogging torque reduction in surface-mounted permanent magnet machines with nonuniform slot distributions

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