CN103279607A - 一种计算永磁电动机气隙磁通密度的方法 - Google Patents

一种计算永磁电动机气隙磁通密度的方法 Download PDF

Info

Publication number
CN103279607A
CN103279607A CN2013102070502A CN201310207050A CN103279607A CN 103279607 A CN103279607 A CN 103279607A CN 2013102070502 A CN2013102070502 A CN 2013102070502A CN 201310207050 A CN201310207050 A CN 201310207050A CN 103279607 A CN103279607 A CN 103279607A
Authority
CN
China
Prior art keywords
air gap
magnetic
flux density
permanent magnet
radially
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
CN2013102070502A
Other languages
English (en)
Other versions
CN103279607B (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.)
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Original Assignee
Changchun Institute of Optics Fine Mechanics and Physics of CAS
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 Changchun Institute of Optics Fine Mechanics and Physics of CAS filed Critical Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority to CN201310207050.2A priority Critical patent/CN103279607B/zh
Publication of CN103279607A publication Critical patent/CN103279607A/zh
Application granted granted Critical
Publication of CN103279607B publication Critical patent/CN103279607B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/80Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
    • Y02T10/82Elements for improving aerodynamics

Landscapes

  • Permanent Field Magnets Of Synchronous Machinery (AREA)

Abstract

本发明涉及一种计算永磁电动机气隙磁通密度的方法,属于电动机技术领域,包括以下步骤:根据永磁电动机物理模型,分别建立周向磁路计算模型、径向磁路计算模型,并对永磁电动机周向磁路气隙磁场和径向磁路气隙磁场分别进行电磁场求解,得到周向和径向磁路气隙磁通密度分布曲线,引入周向权数和径向权数,并结合周向权数和径向权数计算电动机气隙磁通密度。本发明的方法不受电动机结构形式、尺寸限制,而且计算准确,适合于永磁电动机的设计与分析。

Description

一种计算永磁电动机气隙磁通密度的方法
技术领域
本发明属于电动机技术领域,具体涉及一种计算永磁电动机气隙磁通密度的方法。
背景技术
永磁电动机广泛应用于工农业各领域,设计出满足性能要求的电动机是设计者的基本要求。气隙磁通密度是永磁电动机电磁设计最重要的参数之一,准确的计算气隙磁通密度不但可以提高永磁电动机设计效率、减少试制次数,而且能够减少设计和制造成本。
目前,气隙磁通密度的计算主要方法有等效磁路法,算法简单,但许多设计系数需要依靠经验取值,精度不高;解析法,关键在于建立正确的数学模型。一般仅适用于边界形状规则、媒质线性的场合。尤其是对于扁平式电动机,采用传统方法计算气隙磁通密度精度很差。近年来,随着有限元方法不断完善,采用有限元方法可以较准确的分析永磁电动机的电磁场。因此,借助有限元方法对永磁电动机气隙磁通密度进行计算,是一种行之有效的方法。利用有限元方法对电动机进行电磁场分析计算气隙磁通密度,有文献提出的利用径向极弧系数和周向极弧系数相加求平均值或相乘计算气隙磁通密度,对于某些特定结构尺寸的永磁电动机具有较好的精度,但具有一定局限性。
发明内容
本发明要解决现有永磁电动机气隙磁通密度计算不准确和通用性不强的问题,提出一种计算永磁电动机气隙磁通密度的方法。
为了解决上述技术问题,本发明的技术方案具体如下:
一种计算永磁电动机气隙磁通密度的方法,包括以下步骤:根据永磁电动机物理模型,分别建立周向磁路计算模型、径向磁路计算模型,并对永磁电动机周向磁路气隙磁场和径向磁路气隙磁场分别进行电磁场求解,得到周向和径向磁路气隙磁通密度分布曲线,引入周向权数和径向权数,并结合周向权数和径向权数计算电动机气隙磁通密度。
上述技术方案中,本发明的计算永磁电动机气隙磁通密度的方法的具体步骤为:
第I步:对永磁电动机周向磁路进行建模,进行电磁场求解,得到周向磁路气隙磁通密度分布曲线;
第II步:根据周向磁路气隙磁通密度分布曲线,计算周向磁路气隙磁通密度平均值Bcirc_av,最大值Bcirc_max,周向磁路极弧系数αcirc,其中
Figure BDA00003270903300021
第III步:对永磁电动机径向磁路进行建模,进行电磁场求解,得到径向气隙磁通密度分布曲线;
第IV步:根据径向气隙磁通密度分布曲线,计算径向磁路气隙磁通密度平均值Brad_av,最大值Brad_max和径向磁路极弧系数αrad;其中
Figure BDA00003270903300022
第V步:计算加权气隙磁通密度最大值Bweight_max,其中 B weight _ max = ω circ B circ _ max + ω rad B rad _ max ω circ + ω rad
所述周向权数ωcirc=1;所述径向权数ωrad等于电动机电枢直径与电枢长度的比值;
第VI步:计算加权极弧系数αweight,其中 α weight = ω circ α circ + ω rad α rad ω circ + ω rad ;
第VII步:计算电动机气隙磁通密度平均值Bav,其中Bav=Bweight_maxαweight
上述技术方案中,所述气隙磁通密度的计算,是通过对周向磁路与径向磁路分别进行二维电磁场分析综合计算得到;或者是通过周向磁路气隙磁场与径向磁路气隙磁场进行二维电磁场分析加权计算得到。
上述技术方案中,所述二维电磁场分析采用有限元方法进行求解。
上述技术方案中,所述有限元方法通过有限元分析软件实现。
上述技术方案中,所述加权极弧系数αweight是通过周向磁路极弧系数αcirc与径向磁路极弧系数αrad加权计算得到。
上述技术方案中,在建立周向磁路计算模型和径向磁路计算模型之前,还包括步骤:将气隙平均分成N等分,N为2、3、4、5或6;周向、径向气隙磁通密度为N个位置上的磁通密度的平均值。
本发明具有以下的有益效果:
本发明的计算永磁电动机气隙磁通密度的方法,基于有限元方法进行电磁场分析,考虑电动机结构尺寸变化对气隙磁通密度的影响,引入周向权数和径向权数来表征电枢尺寸变化的影响来计算气隙磁通密度。本发明的方法不受电动机结构形式、尺寸限制,而且计算准确,适合于永磁电动机的设计与分析。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1是本发明的计算永磁电动机气隙磁通密度的方法的流程图;
图2是径向磁场永磁电动机周向磁路磁场计算模型;
图3是径向磁场永磁电动机径向磁路磁场计算模型;
图4是轴向磁场永磁电动机径向磁路磁场计算模型;
图5是轴向磁场永磁电动机周向磁路磁场计算模型;
图中:1、磁轭,2、永磁体,3、气隙,4、铁芯,5、端盖。
具体实施方式
本发明的发明思想为,参照图1,根据电动机物理模型,分别建立周向磁路、径向磁路电动机二维计算模型,采用有限元方法进行电磁场分析,得到周向磁路气隙磁通密度分布曲线和径向磁路气隙磁通密度分布曲线,进而计算周向磁路气隙磁通密度最大值和周向磁路极弧系数,径向磁路气隙磁通密度最大值和径向磁路极弧系数。引入周向权数ωaxis和径向权数ωrad,结合周向权数和径向权数计算电动机气隙磁通密度;具体步骤为:
第I步:对永磁电动机周向磁路进行建模,进行电磁场求解,得到周向磁路气隙磁通密度分布曲线;
第II步:根据周向磁路气隙磁通密度分布曲线,计算周向磁路气隙磁通密度平均值Bcirc_av,最大值Bcirc_max,周向磁路极弧系数αcirc,其中
Figure BDA00003270903300041
第III步:对永磁电动机径向磁路进行建模,进行电磁场求解,得到径向气隙磁通密度分布曲线;
第IV步:根据径向气隙磁通密度分布曲线,计算径向磁路气隙磁通密度平均值Brad_av,最大值Brad_max和径向磁路极弧系数αrad;其中
Figure BDA00003270903300042
第V步:计算加权气隙磁通密度最大值Bweight_max,其中 B weight _ max = ω circ B circ _ max + ω rad B rad _ max ω circ + ω rad
所述周向权数ωcirc=1;所述径向权数ωrad等于电动机电枢直径与电枢长度的比值;
第VI步:计算加权极弧系数αweight,其中 α weight = ω circ α circ + ω rad α rad ω circ + ω rad ;
第VII步:计算电动机气隙磁通密度平均值Bav,其中Bav=Bweight_maxαweight
所述气隙磁通密度的计算是通过对周向磁路与径向磁路分别进行二维电磁场分析综合计算得到。
所述气隙磁通密度的计算是通过周向磁路气隙磁场与径向磁路气隙磁场进行二维电磁场分析加权计算得到。
所述电磁场分析采用有限元方法进行求解。
所述有限元方法通过有限元分析软件实现。
所述加权极弧系数是通过周向磁路气隙磁场极弧系数与径向磁路气隙磁场极弧系数加权计算得到。
对于具有大的气隙的永磁电动机,将大的气隙平均分成N等分,N为2、3、4、5或6;周向、径向气隙磁通密度为N个位置上的磁通密度的平均值。
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明,应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
以下结合附图对本发明的具体实施方式做详细说明。
具体实施例1
如图2,对于径向磁场永磁电动机,建立周向磁路物理模型,为了减小计算量,仅建立一对磁极模型,设置边界条件、材料属性和充磁方向,利用有限元软件,如ANSYS、MAXWELL、COMSOL、FEKO等进行电磁场求解,将气隙平均分成2等分并取平均值,得到周向磁路气隙磁通密度分布曲线,同理,建立径向磁路物理模型,如图3,得到径向磁路气隙磁通密度分布曲线。
根据周向磁路气隙磁通密度分布曲线,计算周向磁路气隙磁通密度平均值Bcirc_av,最大值Bcirc_max,周向磁路极弧系数αcirc,其中
Figure BDA00003270903300051
根据径向磁路气隙磁通密度分布曲线,计算径向磁路气隙磁通密度平均值Brad_av,最大值Brad_max和径向磁路极弧系数αrad;其中
Figure BDA00003270903300061
计算加权气隙磁通密度最大值Bweight_max,其中 B weight _ max = ω circ B circ _ max + ω rad B rad _ max ω circ + ω rad
所述周向权数ωcirc=1;所述径向权数ωrad等于电动机电枢直径与电枢长度的比值;
计算加权极弧系数αweight,其中 α weight = ω circ α circ + ω rad α rad ω circ + ω rad ;
计算电动机气隙磁通密度平均值Bav,其中Bav=Bweight_maxαweight
具体实施例2
如图4,对于轴向磁场永磁电动机,建立周向磁路物理模型,为了减小计算量,仅建立一对磁极模型,设置边界条件、材料属性和充磁方向,利用有限元软件,如ANSYS、MAXWELL、COMSOL、FEKO等进行电磁场求解,将气隙平均分成2等分并取平均值,得到周向磁路气隙磁通密度分布曲线,同理,如图5,建立径向磁路物理模型,得到径向磁路气隙磁通密度分布曲线。
根据周向磁路气隙磁通密度分布曲线,计算周向磁路气隙磁通密度平均值Bcirc_av,最大值Bcirc_max,周向磁路极弧系数αcirc,其中
Figure BDA00003270903300064
根据径向磁路气隙磁通密度分布曲线,计算径向磁路气隙磁通密度平均值Brad_av,最大值Brad_max和径向磁路极弧系数αrad;其中
Figure BDA00003270903300065
计算加权气隙磁通密度最大值Bweight_max,其中 B weight _ max = ω circ B circ _ max + ω rad B rad _ max ω circ + ω rad
所述周向权数ωcirc=1;所述径向权数ωrad等于电动机电枢直径与电枢长度的比值;
计算加权极弧系数αweight,其中 α weight = ω circ α circ + ω rad α rad ω circ + ω rad ;
计算电动机气隙磁通密度平均值Bav,其中Bav=Bweight_maxαweight
具体实施例3
以径向磁场永磁电动机为例,建立电动机三维实体模型,为了减小计算量,仅建立一对磁极模型,设置边界条件、材料属性和充磁方向,利用有限元软件,如ANSYS、MAXWELL、COMSOL、FEKO等进行电磁场求解,将三维磁场计算结果沿圆周方向和轴向分别进行处理,其中气隙平均分成4等分并取平均值,从而得到周向磁路气隙磁通密度分布曲线和径向磁路气隙磁通密度分布曲线。
根据周向磁路气隙磁通密度分布曲线,计算周向磁路气隙磁通密度平均值Bcirc_av,最大值Bcirc_max,周向磁路极弧系数αcirc,其中
Figure BDA00003270903300072
根据径向磁路气隙磁通密度分布曲线,计算径向磁路气隙磁通密度平均值Brad_av,最大值Brad_max和径向磁路极弧系数αrad;其中
Figure BDA00003270903300073
计算加权气隙磁通密度最大值Bweight_max,其中 B weight _ max = ω circ B circ _ max + ω rad B rad _ max ω circ + ω rad
所述周向权数ωcirc=1;所述径向权数ωrad等于电动机电枢直径与电枢长度的比值;
计算加权极弧系数αweight,其中 α weight = ω circ α circ + ω rad α rad ω circ + ω rad ;
计算电动机气隙磁通密度平均值Bav,其中Bav=Bweight_maxαweight
显然,上述实施例仅仅是为清楚地说明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (7)

1.一种计算永磁电动机气隙磁通密度的方法,其特征在于,包括以下步骤:根据永磁电动机物理模型,分别建立周向磁路计算模型、径向磁路计算模型,并对永磁电动机周向磁路气隙磁场和径向磁路气隙磁场分别进行电磁场求解,得到周向和径向磁路气隙磁通密度分布曲线,引入周向权数和径向权数,并结合周向权数和径向权数计算电动机气隙磁通密度。
2.根据权利要求1所述的计算永磁电动机气隙磁通密度的方法,其特征在于,具体步骤为:
第I步:对永磁电动机周向磁路进行建模,进行电磁场求解,得到周向磁路气隙磁通密度分布曲线;
第II步:根据周向磁路气隙磁通密度分布曲线,计算周向磁路气隙磁通密度平均值Bcirc_av,最大值Bcirc_max,周向磁路极弧系数αcirc,其中
Figure FDA00003270903200011
第III步:对永磁电动机径向磁路进行建模,进行电磁场求解,得到径向气隙磁通密度分布曲线;
第IV步:根据径向气隙磁通密度分布曲线,计算径向磁路气隙磁通密度平均值Brad_av,最大值Brad_max和径向磁路极弧系数αrad;其中
Figure FDA00003270903200012
第V步:计算加权气隙磁通密度最大值Bweight_max,其中 B weight _ max = ω circ B circ _ max + ω rad B rad _ max ω circ + ω rad
所述周向权数ωcirc=1;所述径向权数ωrad等于电动机电枢直径与电枢长度的比值;
第VI步:计算加权极弧系数αweight,其中 α weight = ω circ α circ + ω rad α rad ω circ + ω rad ;
第VII步:计算电动机气隙磁通密度平均值Bav,其中Bav=Bweight_maxαweight
3.根据权利要求1或2所述的计算永磁电动机气隙磁通密度的方法,其特征在于,所述气隙磁通密度的计算,是通过对周向磁路与径向磁路分别进行二维或三维电磁场分析综合计算得到;或者是通过周向磁路气隙磁场与径向磁路气隙磁场进行二维或三维电磁场分析加权计算得到。
4.根据权利要求3所述的计算永磁电动机气隙磁通密度的方法,其特征在于,所述二维或三维电磁场分析采用有限元方法进行求解。
5.根据权利要求4所述的计算永磁电动机气隙磁通密度的方法,其特征在于,所述有限元方法通过有限元分析软件实现。
6.根据权利要求2所述的计算永磁电动机气隙磁通密度的方法,其特征在于,所述加权极弧系数αweight是通过周向磁路极弧系数αcirc与径向磁路极弧系数αrad加权计算得到。
7.根据权利要求1或2所述的计算永磁电动机气隙磁通密度的方法,其特征在于:在建立周向磁路计算模型和径向磁路计算模型之前,还包括步骤:将气隙平均分成N等分,N为2、3、4、5或6;周向、径向气隙磁通密度为N个位置上的磁通密度的平均值。
CN201310207050.2A 2013-05-30 2013-05-30 一种计算永磁电动机气隙磁通密度的方法 Expired - Fee Related CN103279607B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310207050.2A CN103279607B (zh) 2013-05-30 2013-05-30 一种计算永磁电动机气隙磁通密度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310207050.2A CN103279607B (zh) 2013-05-30 2013-05-30 一种计算永磁电动机气隙磁通密度的方法

Publications (2)

Publication Number Publication Date
CN103279607A true CN103279607A (zh) 2013-09-04
CN103279607B CN103279607B (zh) 2015-12-02

Family

ID=49062124

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310207050.2A Expired - Fee Related CN103279607B (zh) 2013-05-30 2013-05-30 一种计算永磁电动机气隙磁通密度的方法

Country Status (1)

Country Link
CN (1) CN103279607B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104063556A (zh) * 2014-07-07 2014-09-24 东南大学 一种轴向永磁电机等效磁路模型的建模方法
CN105425174A (zh) * 2015-11-24 2016-03-23 广东电网有限责任公司电力科学研究院 一种获取心式变压器漏磁场的方法及装置
CN105720874A (zh) * 2014-12-04 2016-06-29 同济大学 基于分布参数的电机气隙磁场建模方法及其应用
CN105740590A (zh) * 2016-04-13 2016-07-06 上海信耀电子有限公司 一种获取伺服电机空载工作点的系统及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2138836A2 (en) * 2008-06-26 2009-12-30 Pratt & Whitney Canada Corp. Permanent magnet rotor crack detection
CN101957884A (zh) * 2010-03-01 2011-01-26 南京航空航天大学 混合励磁同步电机发电系统建模方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2138836A2 (en) * 2008-06-26 2009-12-30 Pratt & Whitney Canada Corp. Permanent magnet rotor crack detection
CN101957884A (zh) * 2010-03-01 2011-01-26 南京航空航天大学 混合励磁同步电机发电系统建模方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周兰: "基于ANSOFT电动汽车用永磁同步电机电磁设计及性能研究", 《中国优秀硕士学位论文全文数据库》 *
闫杰: "盘式无铁心永磁同步电机磁场分析及磁钢结构优化", 《中国优秀硕士学位论文全文数据库》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104063556A (zh) * 2014-07-07 2014-09-24 东南大学 一种轴向永磁电机等效磁路模型的建模方法
CN105720874A (zh) * 2014-12-04 2016-06-29 同济大学 基于分布参数的电机气隙磁场建模方法及其应用
CN105720874B (zh) * 2014-12-04 2018-07-06 同济大学 基于分布参数的电机气隙磁场建模方法及其应用
CN105425174A (zh) * 2015-11-24 2016-03-23 广东电网有限责任公司电力科学研究院 一种获取心式变压器漏磁场的方法及装置
CN105740590A (zh) * 2016-04-13 2016-07-06 上海信耀电子有限公司 一种获取伺服电机空载工作点的系统及方法

Also Published As

Publication number Publication date
CN103279607B (zh) 2015-12-02

Similar Documents

Publication Publication Date Title
CN109600006B (zh) 一种用于表贴式永磁电机电磁设计的求解方法
CN108563912B (zh) 一种表贴式永磁同步电机气隙磁场的解析方法
CN104063556A (zh) 一种轴向永磁电机等效磁路模型的建模方法
CN103279607B (zh) 一种计算永磁电动机气隙磁通密度的方法
CN103984864B (zh) 一种含有不等厚磁极的表贴式永磁电机磁场优化方法
CN104462624A (zh) 一种基于多场耦合的电机温度场数据处理方法
CN107612256A (zh) 一种磁极分段型表贴式永磁同步电机的优化设计方法
CN104182580A (zh) 一种永磁风力发电机中涡流损耗分析方法
CN107196433A (zh) 一种减少永磁电机齿槽转矩的装置、永磁电机及方法
CN112182844B (zh) 一种永磁磁力耦合器有效磁势计算方法
CN103823926B (zh) 解析算法对永磁电机护套磁导率优化设计
Yang et al. Nonlinear modeling of transmission performance for permanent magnet eddy current coupler
CN105871175B (zh) 一种简化的轴向磁通永磁涡流联轴器转矩计算方法
Guo et al. A performance prediction model for permanent magnet eddy-current couplings based on the air-gap magnetic field distribution
CN105787158A (zh) 基于等效磁路法的永磁同步直线电机建模与特性分析方法
CN113507189B (zh) 内置式永磁同步电机设计参数优化方法
CN106021863A (zh) 轴向磁通永磁涡流联轴器电磁转矩解析算法
Xu et al. Optimum design of interior permanent magnet synchronous motor using Taguchi method
CN104993626A (zh) 一种能够有效减小电机铁耗的内置式永磁电机
CN104362917A (zh) 车用交流发电机漏磁问题的优化设计方法
CN204361868U (zh) 永磁转子正弦充磁装置
CN202737606U (zh) 内置式转子磁钢凸极装置
CN105356804A (zh) 提高同步磁阻电机凸极比的方法
CN206894472U (zh) 非均匀分块式永磁直线同步电机励磁拓扑结构
Lim et al. Magnetising fixture design for optimal magnetisation orientation of ring‐type magnet in surface‐mounted permanent magnet motor

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151202

Termination date: 20170530