CN114154367A - 一种电机齿部磁密比值和最优裂比的确定方法 - Google Patents

一种电机齿部磁密比值和最优裂比的确定方法 Download PDF

Info

Publication number
CN114154367A
CN114154367A CN202111360354.3A CN202111360354A CN114154367A CN 114154367 A CN114154367 A CN 114154367A CN 202111360354 A CN202111360354 A CN 202111360354A CN 114154367 A CN114154367 A CN 114154367A
Authority
CN
China
Prior art keywords
motor
ratio
stator
slot
obtaining
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
CN202111360354.3A
Other languages
English (en)
Other versions
CN114154367B (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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202111360354.3A priority Critical patent/CN114154367B/zh
Publication of CN114154367A publication Critical patent/CN114154367A/zh
Priority to US18/695,934 priority patent/US20240265164A1/en
Priority to PCT/CN2022/080819 priority patent/WO2023087581A1/zh
Application granted granted Critical
Publication of CN114154367B publication Critical patent/CN114154367B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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/22Rotating parts of the magnetic circuit
    • H02K1/24Rotor cores with salient poles ; Variable reluctance rotors
    • H02K1/246Variable reluctance rotors
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • 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/14Stator cores with salient poles
    • H02K1/146Stator cores with salient poles consisting of a generally annular yoke with salient poles
    • 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
    • H02K1/165Shape, form or location of the slots
    • 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/22Rotating parts of the magnetic circuit
    • H02K1/26Rotor cores with slots for windings
    • H02K1/265Shape, form or location of the slots
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K21/00Synchronous motors having permanent magnets; Synchronous generators having permanent magnets
    • H02K21/02Details
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K2213/00Specific aspects, not otherwise provided for and not covered by codes H02K2201/00 - H02K2211/00
    • H02K2213/03Machines characterised by numerical values, ranges, mathematical expressions or similar information

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Power Engineering (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Iron Core Of Rotating Electric Machines (AREA)
  • Mathematical Physics (AREA)

Abstract

本发明公开了一种电机齿部磁密比值和最优裂比的确定方法,属于电机设计领域。包括:在给定不同每槽串联匝数的条件下,通过解析计算,得出多条以电机裂比为自变量的定子槽面积及电流密度曲线;利用磁密比值,得到在不同裂比情况下电机的齿部及轭部厚度,获得电机的几何模型,然后对电机模型进行有限元仿真,得出多条以电机裂比为自变量的平均转矩曲线;获取每条平均转矩曲线最大值对应的电机裂比,代入对应电流密度曲线,得到不同每槽串联匝数下平均转矩最大时电流密度;获取不超过电流密度最大值且平均转矩最大时的电机裂比,作为电机的最优裂比。本发明通过有限元仿真结合解析计算,实现电机最优裂比的快速准确确定。

Description

一种电机齿部磁密比值和最优裂比的确定方法
技术领域
本发明属于电机设计领域,更具体地,涉及一种电机齿部磁密比值和最优裂比的确定方法。
背景技术
电机的高效节能已成为全社会关注的焦点。据统计,电机的耗电量约占全社会总用电量的64%左右,随着新能源汽车等新兴产业的飞速发展,电机在日常生活、工业等领域的覆盖量仍在稳步攀升。巨大市场的需求量和节能减排的政策要求,既给电机市场带来了新的朝气,也对其高效运行提出了新的挑战。因此,发展高效电机对于推动节能减排、实现我国经济绿色发展有着重要的意义。
传统的异步电机,无论是鼠笼型还是绕线型,转子电流引起的焦耳损耗使得其效率难以得到进一步提升,因此,无法满足高效节能的需求。同步磁阻电机转子上无鼠笼结构,显著降低了转子的损耗,提高了电机的效率。永磁电机在转子中加入永磁体,转矩密度得到显著提高,同时也改善了同步磁阻电机仅有定子绕组单边励磁导致功率因数低的不足。因此,在许多工业领域,传统的异步电机逐渐向同步磁阻电机或永磁电机过渡。
一般来说,在进行同基座号产品替换时,多会采用传统的异步电机定子。但由于不同类型电机的结构、工作原理、磁路饱和程度均有所不同,因此需要对电机的定转子尺寸重新设计,以最大程度发挥新电机拓扑的优势。这其中,电机裂比(即定子内外径之比)的选择就显得尤为重要,对于转矩的影响非常明显。传统的电机裂比计算方法大多依赖解析公式,建立转矩和电机裂比之间的表达式。但对于同步磁阻电机或永磁电机来说,其定转子局部饱和严重,影响了转矩计算结果的准确性,解析计算得到的最优裂比通常并非实际最优值,常常无法满足设计需求。且定子轭部厚度和齿宽均随电机裂比的变化而改变,电机的槽面积和定子铜耗因而发生变化,无法保证绕组温升在相同冷却条件仍然满足设计要求,需要经过验证从而导致整个设计流程十分繁杂。
发明内容
针对现有技术的缺陷和改进需求,本发明提供了一种电机齿部磁密比值和最优裂比的确定方法,其目的在于通过有限元仿真结合解析计算,实现电机最优裂比的快速准确确定。
为实现上述目的,按照本发明的第一方面,提供了一种电机齿部磁密比值的确定方法,所述电机为定子有槽径向磁通交流电机,该方法包括:
S1.在定子铜耗固定的前提下,获取以电机裂比为自变量的定子槽面积表达式;
S2.获取以定子轭部厚度和齿宽为自变量的定子槽面积表达式,获取以电机裂比和齿部磁密比值为自变量的定子轭部厚度表达式,获取以电机裂比和齿部磁密比值为自变量的齿宽表达式;
S3.将步骤S2得到定子轭部厚度表达式和齿宽表达式代入步骤S2得到定子槽面积表达式,并按齿部磁密比值进行降幂整理,得到以电机裂比和齿部磁密比值为自变量的定子槽面积表达式;
S4.当电机裂比给定时,利用固定定子铜耗条件下步骤S1得到定子槽面积与电机裂比之间的表达式,计算出相应的定子槽面积,再代入以步骤S3得到电机裂比和齿部磁密比值为自变量的定子槽面积表达式,反演出齿部磁密比值。
优选地,步骤S1中,定子槽面积的计算公式如下:
Figure BDA0003358998820000031
其中,
Figure BDA0003358998820000032
表示定子相电流幅值,ncs表示每槽串联匝数,Sslot表示定子槽面积,ρCu表示铜的电阻率,Lstk表示叠片长度,k表示端部长度与极距的比值,对于分布绕组,k取1.5-1.7,对于分数槽集中绕组,k取
Figure BDA0003358998820000033
De表示定子外径,χ表示电机裂比,p表示电机的极对数,Qs表示定子槽数,PCu表示定子铜耗,kfill表示槽满率。
有益效果:针对现有电机裂比的优化方案常需考虑平均转矩与效率的权衡问题,且忽略了在不同电机裂比下对应的定子电流密度不同而引起不同的绕组温升。本发明通过将电机的定子铜耗视为常数,并根据这一条件推导出定子槽面积关于电机裂比的表达式,进而得到电流密度随电机裂比的变化曲线。由于定子铜耗一定,电机的总损耗基本保持不变,因而在实现中仅需考虑电机的平均转矩最优,而无需进一步考虑电机的效率。
优选地,齿部磁密比值β的计算公式如下:
Figure BDA0003358998820000034
Figure BDA0003358998820000035
其中,fa表示二次幂系数,fb表示一次幂系数,fc表示常数项。
有益效果:针对现有在设计定子槽形尺寸时常需参考以往的设计经验,初步选取一个定子齿部磁密的值,在改变电机的裂比时,往往需要对初始值进行调整,计算流程较为繁杂。本发明通过将齿部磁密比值表示为关于电机裂比和定子槽面积的函数,由于电机裂比和定子槽面积的一一对应关系已在固定定子铜耗为限定下得到,因而实现了齿部磁密比值的快速确定。
为实现上述目的,按照本发明的第二方面,提供了一种电机最优裂比的确定方法,所述电机为定子有槽径向磁通交流电机,该方法包括:
S1.在给定不同每槽串联匝数的条件下,通过解析计算,得出多条以电机裂比为自变量的电流密度曲线,每条曲线和一个每槽串联匝数对应;通过对电机模型进行有限元仿真,得出多条以电机裂比为自变量的平均转矩曲线,每条曲线和一个每槽串联匝数对应;
S2.获取每条平均转矩曲线最大值对应的电机裂比,代入对应电流密度曲线,得到不同每槽串联匝数下平均转矩最大时电流密度;
S3.获取不超过电流密度最大值且平均转矩最大时的电机裂比,作为电机的最优裂比。
优选地,电流密度J的计算公式如下:
Figure BDA0003358998820000041
其中,
Figure BDA0003358998820000042
表示定子相电流幅值,ncs表示每槽串联匝数,Sslot表示定子槽面积,kfill表示槽满率。
优选地,所述以电机裂比为自变量的平均转矩曲线通过以下方式获得:
(1)给定每槽串联匝数,采用如第一方面所述的方法,得到该每槽串联匝数下齿部磁密比值;
(2)将该齿部磁密比值和给定电机裂比代入定子轭部厚度表达式和齿宽表达式,得到该每槽串联匝数下电机裂比对应定子轭部厚度和齿宽;
(3)根据该每槽串联匝数下该电机裂比对应定子轭部厚度和齿宽建立对应电机模型,进行有限元仿真,后处理得到该电机模型对应的平均转矩值;
(4)改变电机裂比,重复步骤(2)-(3),得到该每槽串联匝数下以电机裂比为自变量的平均转矩曲线。
有益效果:本发明通过建立电机的参数化模型,由于电机的定子轭部厚度、齿宽随电机裂比的变化而自动调节,实现电机建模及后处理的快速化。
优选地,所述给定不同每槽串联匝数通过以下方式确定:
(1)根据电机相电压U,确定每槽串联匝数的基数:
Figure BDA0003358998820000051
其中,m表示电机相数,f表示频率,Qs表示定子槽数,kN1表示绕组系数,Φ1表示每极磁通量;
(2)确定每槽串联匝数的变化步长:对于单层绕组,变化步长取值1/a;对于双层绕组,变化步长取值2/a,其中,a表示绕组并联支路数;
(3)对每槽串联匝数的基数,按确定的变化步长进行上下波动,从而确定每槽串联匝数的取值范围。
有益效果:现有每槽串联匝数变化范围未知的问题,导致有限元仿真的计算量过大。本发明通过将每槽串联匝数与相电压相联系,初步确定每槽串联匝数,实现了每槽串联匝数范围及步长的确定,减小了有限元仿真的计算时间。
为实现上述目的,按照本发明的第三方面,提供了一种计算机可读存储介质,包括存储的计算机程序;所述计算机程序被处理器执行时,控制所述计算机可读存储介质所在设备执行第一方面所述的电机齿部磁密比值的确定方法,或者,执行第二方面所述的电机最优裂比的确定方法。
总体而言,通过本发明所构思的以上技术方案,能够取得以下有益效果:
(1)现有设计定子槽形尺寸前需根据以往的设计经验,初步选取电机的定子轭部及齿部磁密,根据估计值计算定子轭部厚度及齿宽,并进行磁路计算及性能求解。当电机裂比的选取改变时,电机内部气隙磁密随之变化,根据估算值得到的定子轭部厚度和齿宽往往难以获得之前选取的定子轭部及齿部磁密,还需对初选值进行调整,重新计算性能,计算流程较为繁杂。本发明通过将定子轭部厚度和齿宽表示为电机裂比及齿部磁密比值的函数关系,而齿部磁密比值建立了定子轭部及齿部磁密与气隙磁密之间的关系,有效地避免了上述问题。然后利用定子轭部厚度和齿宽得到电机定子槽面积关于电机裂比和齿部磁密比值的表达式。通过相同定子铜耗的前提条件,得到不同电机裂比对应的定子槽面积,实现了齿部磁密比值及相应定子轭部厚度、齿宽的快速确定。
(2)现有电机裂比的计算方法常依赖于解析公式,建立转矩和裂比的函数关系。然而,电机内部通常存在定转子局部饱和问题,铁芯上的磁压降不可忽略,因而削弱了气隙磁场。而解析法难以考虑局部饱和问题,造成得到的转矩结果偏高,难以满足较高的计算精度需求。此外,定转子开槽引起的气隙磁场谐波含量较高,而解析公式中的计算量仅考虑了基波成分,忽略了其他次谐波含量,利用解析公式难以得到准确的结果。本发明通过有限元仿真与解析计算相结合的技术手段,考虑了铁芯饱和及谐波的影响,实现最大平均转矩下最优电机裂比的快速精确确定。
附图说明
图1是本发明提供的一种电机最优裂比的确定方法流程图。
图2是本发明提供的每极同步磁阻电机的截面图。
图3是本发明提供的梯形定子槽面积计算的示意图。
图4是本发明提供的48槽8极同步磁阻电机中,每槽串联匝数分别取3.5、4、4.5时定子槽面积随电机裂比的变化关系。
图5是本发明提供的48槽8极同步磁阻电机中,每槽串联匝数分别取3.5、4、4.5时平均转矩和电流密度随电机裂比的变化关系。
图6是本发明提供的斜肩圆底槽定子槽面积计算的示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明提供了一种电机齿部磁密比值的确定方法,所述电机为定子有槽径向磁通交流电机,该方法包括:
步骤S1.在定子铜耗固定的前提下(即电机总损耗近似保持不变),获取以电机裂比为自变量的定子槽面积表达式。
优选地,步骤S1中,定子槽面积的计算公式如下:
Figure BDA0003358998820000071
其中,
Figure BDA0003358998820000072
表示定子相电流幅值,ncs表示每槽串联匝数,Sslot表示定子槽面积,ρCu表示铜的电阻率,Lstk表示叠片长度,k表示端部长度与极距的比值,对于分布绕组,k取1.5-1.7,对于分数槽集中绕组,k取
Figure BDA0003358998820000073
De表示定子外径,χ表示电机裂比,p表示电机的极对数,Qs表示定子槽数,PCu表示定子铜耗,kfill表示槽满率。
步骤S2.获取以定子轭部厚度和齿宽为自变量的定子槽面积表达式,获取以电机裂比和齿部磁密比值为自变量的定子轭部厚度表达式,获取以电机裂比和齿部磁密比值为自变量的齿宽表达式。
步骤S3.将步骤S2得到定子轭部厚度表达式和齿宽表达式代入步骤S2得到定子槽面积表达式,并按齿部磁密比值进行降幂整理,得到以电机裂比和齿部磁密比值为自变量的定子槽面积表达式。
步骤S4.当电机裂比给定时,利用固定定子铜耗条件下步骤S1得到定子槽面积与电机裂比之间的表达式,计算出相应的定子槽面积,再代入以步骤S3得到电机裂比和齿部磁密比值为自变量的定子槽面积表达式,反演出齿部磁密比值。
优选地,齿部磁密比值β的计算公式如下:
Figure BDA0003358998820000081
Figure BDA0003358998820000082
其中,fa表示二次幂系数,fb表示一次幂系数,fc表示常数项。
如图1所示,本发明提供了一种电机最优裂比的确定方法,所述电机为定子有槽径向磁通交流电机,该方法包括:
步骤S1.在给定不同每槽串联匝数的条件下,通过解析计算,得出多条以电机裂比为自变量的定子槽面积及电流密度曲线,每条曲线和一个每槽串联匝数对应;通过对电机模型进行有限元仿真,得出多条以电机裂比为自变量的平均转矩曲线,每条曲线和一个每槽串联匝数对应。
优选地,电流密度J的计算公式如下:
Figure BDA0003358998820000083
其中,
Figure BDA0003358998820000084
表示定子相电流幅值,ncs表示每槽串联匝数,Sslot表示定子槽面积,kfill表示槽满率。
优选地,所述以电机裂比为自变量的平均转矩曲线通过以下方式获得:
(1)给定每槽串联匝数,采用如上述方法,得到该每槽串联匝数下齿部磁密比值;
(2)将该齿部磁密比值和给定电机裂比代入定子轭部厚度表达式和齿宽表达式,得到该每槽串联匝数下该电机裂比对应定子轭部厚度和齿宽;
(3)根据该每槽串联匝数下该电机裂比对应定子轭部厚度和齿宽建立对应电机模型,进行有限元仿真,后处理得到该电机模型对应的平均转矩值;
(4)改变电机裂比,重复步骤(2)-(3),得到该每槽串联匝数下以电机裂比为自变量的平均转矩曲线。
优选地,所述给定不同每槽串联匝数通过以下方式确定:
(1)根据电机相电压U,确定每槽串联匝数的基数:
Figure BDA0003358998820000091
其中,m表示电机相数,f表示频率,Qs表示定子槽数,kN1表示绕组系数,Φ1表示每极磁通量。
Ns取决于电机的相电压U,计算公式如下:
Figure BDA0003358998820000092
(2)确定每槽串联匝数的变化步长:对于单层绕组,变化步长取值1/a,其中,a是绕组并联支路数;对于双层绕组,变化步长取值2/a。
每槽串联匝数的计算公式如下:
Figure BDA0003358998820000093
其中,a为每相并联支路数,对于单层绕组,其最大值等于极对数。对于双层绕组,其最大值等于极数。因此,ncs的变化步长由每相并联支路数决定。
(3)对每槽串联匝数的基数,按确定的变化步长进行上下波动,从而确定每槽串联匝数的取值范围。
步骤S2.获取每条平均转矩曲线最大值对应的电机裂比,代入对应电流密度曲线,得到不同每槽串联匝数下平均转矩最大时电流密度。
步骤S3.获取不超过电流密度最大值且平均转矩最大时的电机裂比,作为电机的最优裂比。
以图2所示的梯形槽定子48槽8极同步磁阻电机为例,将定子内径Di、外径De、齿宽wt、轭厚hbi、定子槽面积Sslot标示在图中。
电机裂比定义为定子内径与外径之比:
χ=Di/De (8)
电机的转矩可以表示为:
Figure BDA0003358998820000101
式中,Lstk为电机叠片长度,Bg为基波气隙磁密幅值,Ks为线电流密度,可以进一步表示为:
Figure BDA0003358998820000102
式中,kw为基波绕组系数,Qs为定子槽数,ncs为每槽串联匝数,
Figure BDA0003358998820000103
为相电流幅值。
进一步地,将式(10)代入式(9)中,得到转矩方程:
Figure BDA0003358998820000104
可以发现,当电机的定子外径固定时,输出转矩与裂比大小有关,因此,要实现最大转矩的优化,需要找到合适的电机裂比。
气隙磁通:
Φg=BgπDiLstk (12)
定子齿部磁通:
Φt=BtQswtLstk (13)
忽略槽漏磁,则气隙磁通全部穿过齿部,即Φg=Φt,并将磁密比值β定义为气隙磁密Bg与定子齿磁密Bt的比值,此时,齿宽可以表示为:
Figure BDA0003358998820000105
通过定子轭部的磁通为气隙主磁通的一半:
Figure BDA0003358998820000111
式中,p为电机的极对数。
将式(14)和式(15)联立,得到定子轭部厚度的表达式为:
Figure BDA0003358998820000112
式中,α定义为齿部磁密与轭部磁密的比值,考虑到铁芯叠片的膝点,α为常数且约等于1.16。
根据图3,将梯形槽简化为平行槽,忽略槽楔尺寸,得到定子槽面积的表达式为:
Figure BDA0003358998820000113
将式(14)、式(16)代入式(17)中,并按β的降幂整理,得:
Figure BDA0003358998820000114
式中:
Figure BDA0003358998820000115
给定槽极配合,将槽数与极数固定,且规定定子外径为常数,则定子槽面积仅与β和χ有关。
进一步地,引入定子铜耗为常数的限制条件,并将铜耗表示为:
Figure BDA0003358998820000116
式中,ρcu为铜的电阻率,kfill为槽满率,m为相数,Lcond为铜导体的长度,包括叠片长度与端部长度。其中,在工艺上,分布绕组的端部长度可近似为极距的1.5-1.7倍,这里取1.6。因此,Lcond也可以表示为裂比的函数:
Figure BDA0003358998820000121
在式(20)中,除Sslot、ncs和χ外,其余参数均为常数,因此,当每槽串联匝数给定时,可以得到定子槽面积随电机裂比的变化关系:
Figure BDA0003358998820000122
48槽8极同步磁阻电机Sslot、ncs和χ的关系如图4所示。当电机裂比或每槽串联匝数增加时,定子槽面积也随之增大。
将式(18)与式(19)联立,解得:
Figure BDA0003358998820000123
式中:
Figure BDA0003358998820000124
进一步地,定子齿宽与轭部厚度可以分别由式(14)和式(16)表示为仅以电机裂比为变量的函数。当电机裂比固定时,诸如定子内径、定子齿宽与轭部厚度等几何参数可以全部确定。对于不同的电机裂比,建立电机模型,然后通过有限元仿真,得到其对应的平均转矩,如图5上半部所示,对应的电流密度如图5下半部所示。显然,对于不同的每槽串联匝数ncs取值,以最大平均转矩为目标,可以找到不同的电机裂比最优值。另外,考虑到定子槽内温升的限制,导体的电流密度存在最大值限制,在实例所示情况下,最大电流密度取22.8A/mm2。当ncs增加时,电流密度趋于增大并有超过这个限制的风险。在图5中,最优的电机裂比与每槽串联匝数分别为0.65和4,对应的最大平均转矩为170.9N.m。
上述实例采用的是梯形槽,为验证方法的普适性,下面采用圆底槽进行说明,如图6所示。此时,式(17)改写为:
Figure BDA0003358998820000131
式中,r为槽底半径,在电机设计中,在槽数确定的情况下,r通常为一常数。
同样地,将式(14)、式(16)代入式(25)中,并按β的降幂整理,得:
Figure BDA0003358998820000132
式中:
Figure BDA0003358998820000133
给定槽极配合,将槽数与极数固定,且规定定子外径为常数,定子槽面积仍然仅与β和χ有关,与梯形槽得到的结论一致。因此,本发明对其他槽型同样适用。
此外,本发明还提供了一种计算机可读存储介质,包括存储的计算机程序;所述计算机程序被处理器执行时,控制所述计算机可读存储介质所在设备执行上述电机齿部磁密比值的确定方法,或者,执行上述电机最优裂比的确定方法。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种电机齿部磁密比值的确定方法,所述电机为定子有槽径向磁通交流电机,其特征在于,该方法包括:
S1.在定子铜耗固定的前提下,获取以电机裂比为自变量的定子槽面积表达式;
S2.获取以定子轭部厚度和齿宽为自变量的定子槽面积表达式,获取以电机裂比和齿部磁密比值为自变量的定子轭部厚度表达式,获取以电机裂比和齿部磁密比值为自变量的齿宽表达式;
S3.将步骤S2得到定子轭部厚度表达式和齿宽表达式代入步骤S2得到定子槽面积表达式,并按齿部磁密比值进行降幂整理,得到以电机裂比和齿部磁密比值为自变量的定子槽面积表达式;
S4.当电机裂比给定时,利用固定定子铜耗条件下步骤S1得到定子槽面积与电机裂比之间的表达式,计算出相应的定子槽面积,再代入以步骤S3得到电机裂比和齿部磁密比值为自变量的定子槽面积表达式,反演出齿部磁密比值。
2.如权利要求1所述的方法,其特征在于,步骤S1中,定子槽面积的计算公式如下:
Figure FDA0003358998810000011
其中,
Figure FDA0003358998810000012
表示定子相电流幅值,ncs表示每槽串联匝数,Sslot表示定子槽面积,ρCu表示铜的电阻率,Lstk表示叠片长度,k表示绕组端部长度与极距的比值,对于分布绕组,k取1.5-1.7,对于分数槽集中绕组,k取
Figure FDA0003358998810000013
De表示定子外径,χ表示电机裂比,p表示电机的极对数,Qs表示定子槽数,PCu表示定子铜耗,kfill表示槽满率。
3.如权利要求1所述的方法,其特征在于,齿部磁密比值β的计算公式如下:
Figure FDA0003358998810000021
Figure FDA0003358998810000022
其中,fa表示二次幂系数,fb表示一次幂系数,fc表示常数项。
4.一种电机最优裂比的确定方法,所述电机为定子有槽径向磁通交流电机,其特征在于,该方法包括:
S1.在给定不同每槽串联匝数的条件下,通过解析计算,得出多条以电机裂比为自变量的电流密度曲线,每条曲线和一个每槽串联匝数对应;通过对电机模型进行有限元仿真,得出多条以电机裂比为自变量的平均转矩曲线,每条曲线和一个每槽串联匝数对应;
S2.获取每条平均转矩曲线最大值对应的电机裂比,代入对应电流密度曲线,得到不同每槽串联匝数下平均转矩最大时电流密度;
S3.获取不超过电流密度最大值且平均转矩最大时的电机裂比,作为电机的最优裂比。
5.如权利要求4所述的方法,其特征在于,电流密度J的计算公式如下:
Figure FDA0003358998810000023
其中,
Figure FDA0003358998810000024
表示定子相电流幅值,ncs表示每槽串联匝数,Sslot表示定子槽面积,kfill表示槽满率。
6.如权利要求4所述的方法,其特征在于,所述以电机裂比为自变量的平均转矩曲线通过以下方式获得:
(1)给定每槽串联匝数,采用如权利要求1至3任一项所述的方法,得到该每槽串联匝数下齿部磁密比值;
(2)将该齿部磁密比值和给定电机裂比代入定子轭部厚度表达式和齿宽表达式,得到该每槽串联匝数下该电机裂比对应定子轭部厚度和齿宽;
(3)根据该每槽串联匝数下该电机裂比对应定子轭部厚度和齿宽建立对应电机模型,进行有限元仿真,后处理得到该电机模型对应的平均转矩值;
(4)改变电机裂比,重复步骤(2)-(3),得到该每槽串联匝数下以电机裂比为自变量的平均转矩曲线。
7.如权利要求4至6任一项所述的方法,其特征在于,所述给定不同每槽串联匝数通过以下方式确定:
(1)根据电机相电压U,确定每槽串联匝数的基数:
Figure FDA0003358998810000031
其中,m表示电机相数,f表示频率,Qs表示定子槽数,kN1表示绕组系数,Φ1表示每极磁通量;
(2)确定每槽串联匝数的变化步长:对于单层绕组,变化步长取值1/a;对于双层绕组,变化步长取值2/a,其中,a表示绕组并联支路数;
(3)对每槽串联匝数的基数,按确定的变化步长进行上下波动,从而确定每槽串联匝数的取值范围。
8.一种计算机可读存储介质,其特征在于,包括存储的计算机程序;所述计算机程序被处理器执行时,控制所述计算机可读存储介质所在设备执行权利要求1至3任一项所述的电机齿部磁密比值的确定方法,或者,执行权利要求4至7任一项所述的一种电机最优裂比的确定方法。
CN202111360354.3A 2021-11-17 2021-11-17 一种电机齿部磁密比值和最优裂比的确定方法 Active CN114154367B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202111360354.3A CN114154367B (zh) 2021-11-17 2021-11-17 一种电机齿部磁密比值和最优裂比的确定方法
US18/695,934 US20240265164A1 (en) 2021-11-17 2022-03-15 Method for determining tooth magnetic flux density ratio and optimal split ratio of motor
PCT/CN2022/080819 WO2023087581A1 (zh) 2021-11-17 2022-03-15 一种电机齿部磁密比值和最优裂比的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111360354.3A CN114154367B (zh) 2021-11-17 2021-11-17 一种电机齿部磁密比值和最优裂比的确定方法

Publications (2)

Publication Number Publication Date
CN114154367A true CN114154367A (zh) 2022-03-08
CN114154367B CN114154367B (zh) 2024-07-16

Family

ID=80456317

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111360354.3A Active CN114154367B (zh) 2021-11-17 2021-11-17 一种电机齿部磁密比值和最优裂比的确定方法

Country Status (3)

Country Link
US (1) US20240265164A1 (zh)
CN (1) CN114154367B (zh)
WO (1) WO2023087581A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023087581A1 (zh) * 2021-11-17 2023-05-25 华中科技大学 一种电机齿部磁密比值和最优裂比的确定方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117195737B (zh) * 2023-09-21 2024-03-22 江苏中车电机有限公司 永磁风力发电机绕组结构及其智能优化方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110851949A (zh) * 2019-09-01 2020-02-28 天津工业大学 分析多层磁障永磁辅助式同步磁阻电机电磁性能的方法
WO2021189594A1 (zh) * 2020-03-24 2021-09-30 江苏大学 一种磁场调制混合励磁电机及其多工作波设计方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5479978B2 (ja) * 2010-03-30 2014-04-23 アイシン・エィ・ダブリュ株式会社 回転電機
US8829742B2 (en) * 2013-01-04 2014-09-09 Xinzhang Wu High efficiency permanent magnet machine
JP5647307B1 (ja) * 2013-06-04 2014-12-24 成田 憲治 直流励磁界磁型同期電動機
US10622875B2 (en) * 2017-06-07 2020-04-14 GM Global Technology Operations LLC Interior permanent magnet electric machine
US10700632B1 (en) * 2019-01-11 2020-06-30 GM Global Technology Operations LLC Method for motor and inverter temperature control
CN110492696B (zh) * 2019-09-05 2021-01-19 合肥工业大学 优化Halbach阵列永磁电机裂比和极弧比的方法
CN114154367B (zh) * 2021-11-17 2024-07-16 华中科技大学 一种电机齿部磁密比值和最优裂比的确定方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110851949A (zh) * 2019-09-01 2020-02-28 天津工业大学 分析多层磁障永磁辅助式同步磁阻电机电磁性能的方法
WO2021189594A1 (zh) * 2020-03-24 2021-09-30 江苏大学 一种磁场调制混合励磁电机及其多工作波设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
姜仁华;刘闯;宁银行;王承宇;: "分数槽集中绕组永磁同步电机设计与分析", 南京航空航天大学学报, no. 03, 15 June 2018 (2018-06-15) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023087581A1 (zh) * 2021-11-17 2023-05-25 华中科技大学 一种电机齿部磁密比值和最优裂比的确定方法

Also Published As

Publication number Publication date
WO2023087581A1 (zh) 2023-05-25
US20240265164A1 (en) 2024-08-08
CN114154367B (zh) 2024-07-16

Similar Documents

Publication Publication Date Title
Wang et al. Cogging torque minimization and torque ripple suppression in surface-mounted permanent magnet synchronous machines using different magnet widths
Labak et al. Designing and prototyping a novel five-phase pancake-shaped axial-flux SRM for electric vehicle application through dynamic FEA incorporating flux-tube modeling
CN114154367B (zh) 一种电机齿部磁密比值和最优裂比的确定方法
Kahourzade et al. Optimal design of Axial-Flux Induction Motors based on an improved analytical model
CN110968969A (zh) 异步电机铁心损耗分析方法
CN111666651A (zh) 永磁电机的退磁风险评估方法
Marfoli et al. Squirrel cage induction motor: A design-based comparison between aluminium and copper cages
Ashwin et al. Design optimization of 15 kW, 2-pole induction motor to achieve IE4 efficiency level with copper die-casting
Di Nardo et al. Rotor slot design of squirrel cage induction motors with improved rated efficiency and starting capability
Gercekcioglu et al. Optimal rotor design of novel axial flux synchronous reluctance motor and validation
Feyzi et al. A comprehensive review on the performance improvement in switched reluctance motor design
Yun et al. Influence of aluminum die-cast rotor porosity on the efficiency of induction machines
Van Khang et al. Loss minimization for form-wound stator winding of a high-speed induction motor
Resketi et al. Assessing the impact of three‐phase rewinding of a failed single‐phase motor on efficiency improvement and energy saving
Alidousti et al. A new rotor prototype for single phase line start permanent magnet synchronous motor based on amendments to a small industrial shaded pole induction motor
Di Nardo et al. Open and closed rotor slots design of single and double cages induction motor
Zhang et al. Optimization and experimental validation of amorphous alloy high-speed asynchronous motor for simultaneous reduction on core and copper losses
Cabuk Simulation of the effect of segmented axial direction magnets on the efficiency of in-wheel permanent magnet brushless DC motors used in light electric vehicles based on finite element method
Özsoy et al. FEM-based analysis of rotor cage material and slot geometry on double air gap axial flux induction motors
Piscini et al. Contribution on AC bar windings losses reduction for a high frequency and high performance machine for aeronautical application
Sprangers et al. Design and optimization tools for high-efficiency three-phase induction motors
CN110162892B (zh) 一种大型同步调相机定子端部漏感计算方法
Ugale et al. A hybrid interior rotor high-performance line start permanent magnet synchronous motor
El-Kharahi et al. Toward including the effect of manufacturing processes in the pre-estimated losses of the switched reluctance motor
Năvrăpescu et al. Computation methods for space harmonic effects on single-phase induction motor performance

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