CN112906145A - 电机电磁力与电磁振动噪声分析方法及系统 - Google Patents
电机电磁力与电磁振动噪声分析方法及系统 Download PDFInfo
- Publication number
- CN112906145A CN112906145A CN202011426440.5A CN202011426440A CN112906145A CN 112906145 A CN112906145 A CN 112906145A CN 202011426440 A CN202011426440 A CN 202011426440A CN 112906145 A CN112906145 A CN 112906145A
- Authority
- CN
- China
- Prior art keywords
- motor
- stator
- vibration
- air gap
- slot
- 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
- 238000004458 analytical method Methods 0.000 title claims description 19
- 230000006870 function Effects 0.000 claims abstract description 79
- 238000004364 calculation method Methods 0.000 claims abstract description 68
- 230000004907 flux Effects 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 39
- 238000012546 transfer Methods 0.000 claims abstract description 34
- 230000001133 acceleration Effects 0.000 claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 18
- 230000005284 excitation Effects 0.000 claims abstract description 9
- 230000003068 static effect Effects 0.000 claims abstract description 9
- 238000009826 distribution Methods 0.000 claims description 20
- 238000004088 simulation Methods 0.000 claims description 19
- 238000004804 winding Methods 0.000 claims description 18
- 230000005855 radiation Effects 0.000 claims description 10
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 5
- 239000000126 substance Substances 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 239000004020 conductor Substances 0.000 description 6
- 230000008901 benefit Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Iron Core Of Rotating Electric Machines (AREA)
Abstract
本发明公开了一种电机电磁力及电磁振动噪声分析方法及系统,所述方法包括S100:确定电机的运行工况,获得电机任意时刻的相电流,确定电机的槽磁动势函数和槽矩阵及电机任意时刻的电枢磁动势;S200:建立电机静态有限元模型,并分析该电机永磁体磁动势以及该电机定子和转子的复数气隙磁导函数;S300:获得该电机的气隙磁密,获得所述电机每个定子齿上所受的电磁力;S400:建立电机的三维结构有限元模型,并分析电机的振动传递函数;S500:计算电机在任意电磁力激励下的振动速度和振动加速度,基于振动速度和振动加速度获得电机辐射的声压级和声功率。本发明可对电机各种工况下的电磁力及电磁振动噪声进行快速评估与分析,大幅提升电机振动噪声计算与优化效率。
Description
技术领域
本发明属于电机振动噪声技术领域,更具体地,涉及一种电机电磁力与电磁振动噪声分析方法及系统。
背景技术
电机在社会生产生活的各个领域得到了广泛的应用。近年来,电机的振动噪声逐渐成为衡量电机性能的重要指标。以新能源汽车为代表的交通装备领域和以全电舰船、水下潜航器为代表的国防领域都对电机的振动噪声提出了极高的要求。对新能源汽车而言,驱动电机的振动噪声不仅降低了车辆的乘坐舒适性,还会影响到用户对品牌的认可度和市场竞争力。在全电舰船、水下潜航器等领域,低振动噪声电机系统可提高舰艇的声隐身性和生命力,形成威慑力。在家用电器、医疗器械等领域,出于保护人体健康、降低环境噪声污染的要求,也要求将电机的噪声控制在一定的范围内。因此,降低电机振动和噪声对军民领域都有重要意义。
电机振动噪声与电机内部电磁力分布、电机结构特性以及运行工况密切相关,为了准确对电机振动噪声进行评估,给低振动电机的设计提供有力工具,急需一种能准确快速计算电机电磁力及电磁振动噪声的方法。现有电机振动计算方法主要有解析法、有限元法和半解析法。在解析法中,电机电磁力、固有频率和振动传递函数、电机振动变形及振动加速度都通过解析公式计算得到。解析法效率较高,但电机结构复杂导致解析法计算精度较低,不能准确评估电机的真实振动情况。因此,有限元法成为振动计算的主流,在利用有限元法计算电机振动时,首先需要建立电磁有限元模型,仿真得到电机的电磁力;然后建立电机的三维结构有限元模型,确定边界条件,导入电磁有限元计算得到的电磁力;最后利用有限元软件进行结构动力学仿真得到振动响应,并通过声学边界元仿真得到电机的噪声分布情况。有限元方法计算精度较高且适用于任何复杂结构的电机,但有限元方法涉及多场耦合计算,建模复杂导致计算效率较低,不适合在电机设计初期对电机振动性能进行快速评估,也不能胜任电机多转速工况振动计算的需求。而半解析法则结合了解析法和有限元的优势,在保证计算精度的同时也能提高振动计算效率。但现有振动半解析计算方法中只是将振动响应计算部分用模态叠加法代替,而最为耗时的电磁力计算和传递函数计算部分仍然依赖有限元计算数据,并未从根本上解决振动计算复杂的问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供一种电机电磁力与电磁振动噪声分析方法及系统,其目的在于在保证电机振动计算精度的同时能减小计算量,实现电机振动噪声的快速预测和评估,为电机振动噪声优化与多转速工况振动计算提供有力工具。
为了实现上述目的,按照本发明的一个方面,提供一种电机电磁力及电磁振动噪声分析方法,包括如下步骤:
S100:确定电机的运行工况,根据其运行工况获得电机任意时刻的相电流,根据电机槽极配合、绕组分布确定电机的槽磁动势函数和槽矩阵,并结合所述相电流获得电机任意时刻的电枢磁动势;
S200:建立所述电机的静态有限元模型,并分析获得该电机永磁体磁动势以及该电机定子和转子的复数气隙磁导函数;
S300:根据所述电枢磁动势、永磁体磁动势及定子和转子的复数气隙磁导函数获得该电机的气隙磁密,并基于该气隙磁密获得所述电机每个定子齿上所受的电磁力;
S400:建立所述电机的三维结构有限元模型,并分析获得电机的振动传递函数;
S500:根据所述电机定子齿电磁力和振动传递函数,计算获得所述电机在任意电磁力激励下的振动速度和振动加速度,基于所述振动速度和振动加速度获得电机辐射的声压级和声功率。
进一步地,步骤S100具体包括:
S101:确定电机的槽极配合与绕组分布;
S102:根据所述槽极配合与绕组分布,建立电机槽矩阵;
S103:建立所述电机的槽磁动势函数;
S104:根据电机的运行工况,确定电机的相电流瞬时值;
S105:结合所述相电流瞬时值、电机槽矩阵以及槽磁动势函数获得电机在任意时刻与电流下的电枢磁动势。
进一步地,所述槽磁动势函数为:
其中,Nq(θs)为槽磁动势函数,Δθ为槽宽,θq为第q个槽中心位置,其计算公式为:θq=(q-1)·2π/Qs,θs为任意时刻转子位置,Qs为电机槽数。
进一步地,步骤S200具体包括:
S201~S203:建立无槽定子和无槽转子有限元仿真模型,并施加幅值大小为1的恒定虚拟磁动势,仿真计算定子及转子均无槽时的气隙磁密;
S204~S205:建立有槽定子-无槽转子有限元模型,并施加幅值大小为1的恒定虚拟磁动势,仿真计算定子有槽、转子无槽时的气隙磁密;
S207~S208:建立无槽定子-有槽转子有限元模型,并施加幅值大小为1的恒定虚拟磁动势,仿真计算定子无槽、转子有槽时的气隙磁密;
S206、S209:根据所述定子及转子均无槽时的气隙磁密、定子有槽、转子无槽时的气隙磁密以及定子无槽、转子有槽时的气隙磁密,得到定子和转子的复数气隙磁导函数。
进一步地,所述复数气隙磁导函数为:
Br,slotless和Bt,slotless分别为定子及转子均无槽时的气隙磁密的径向和切向分量;
进一步地,步骤S400具体包括:
S401:确定电机定子铁芯及绕组的等效材料参数;
S402:建立所述电机的三维结构有限元模型;
S403:在电机定子齿上施加幅值为1,沿定子圆周呈周期分布的径向力和切向力;
S404:通过有限元仿真计算得到电机不同齿在单位激振力作用下的振动传递函数。
进一步地,步骤S500具体包括:
S502:将定子齿电磁力进行傅里叶分解,得到定子齿电磁力的频谱分布;
S503:将计算得到的各个定子齿上的电磁力谐波和对应齿的振动传递函数求积,将所有齿的计算结果叠加即可得到电机的总振动变形;
S505:基于位移、速度和加速度之间的关系获得电机表面的振动速度和振动加速度。
按照本发明的第二个方面,提供一种电机电磁力与电磁振动噪声分析系统,包括:
电枢磁动势分析模块:用于确定电机的运行工况,根据其运行工况获得电机任意时刻的相电流,根据电机槽极配合、绕组分布确定电机的槽磁动势函数和槽矩阵,并结合所述相电流获得电机任意时刻的电枢磁动势;
复数气隙磁导函数计算模块:用于建立所述电机的静态有限元模型,并分析获得该电机永磁体磁动势以及该电机定子和转子的复数气隙磁导函数;
电磁力计算模块:用于根据所述电枢磁动势、永磁体磁动势及定子和转子的复数气隙磁导函数获得该电机的气隙磁密,并基于该气隙磁密获得所述电机每个定子齿上所受的电磁力;
振动传递函数分析模块:用于建立所述电机的三维结构有限元模型,并分析获得电机的振动传递函数;
电机振动计算模块:用于根据所述电机定子齿电磁力和振动传递函数,计算获得所述电机在任意电磁力激励下的振动速度和振动加速度;
电机辐射的声压级和声功率计算模块:用于根据所述振动速度和振动加速度获得电机辐射的声压级和声功率。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
1.本发明提出了一种电机电磁力及电磁振动噪声的分析方法,基于此方法建立了电机电磁振动噪声快速计算系统,可以对电机各种工况下的电磁力及电磁振动噪声进行快速评估与分析,其计算精度可以达到传统有限元计算方法相当的水平,但计算时间大幅缩短,为电机振动优化和多转速工况振动快速评估提供了有力工具。
2.本发明提供一种基于静态有限元的电机气隙磁导计算方法,通过构建虚拟磁动势有限元仿真模型,可以快速计算任意结构电机的气隙磁导,解决了电磁气隙磁导计算复杂的难题。
3.本发明创新性的提出了槽磁动势和槽矩阵,用于计算电机在任意电流激励下的电枢磁动势。通过该方法,将考虑高频电流谐波的电枢磁动势的计算时间大幅缩短,同时提高了计算的精度。
4.本发明提出的电磁力及电磁振动计算方法同时考虑了转子磁导、磁场饱和、切向电磁力以及PWM电流谐波,在大幅加快计算速度的同时进一步提高计算精度。
附图说明
图1为本发明实施例电机电磁力及电磁振动噪声计算流程示意图;
图2为本发明实施例中电机电枢磁动势计算流程示意图;
图3为本发明实施例六相36槽6极整数槽双层短距绕组分布示意图;
图4为本发明实施例复数气隙磁导计算流程示意图;
图5为本发明实施例振动传递函数计算流程示意图;
图6为本发明实施例中电机振动计算流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,本发明实施例提供一种电机电磁力及电磁振动噪声计算方法,其包括如下步骤:
(1)整个计算流程开始前需要进行计算平台的初始化,包括设定电机的转速、转矩、功率等,以确定电机的运行工况。根据电机运行工况可以计算电机任意时刻下的相电流Iph,t。
(2)根据电机槽极配合、绕组分布确定电机的槽磁动势函数和槽矩阵,并结合步骤(1)中计算得到的相电流计算电机任意时刻下的电枢磁动势,具体流程如图2所示:
首先根据电机的槽极配合和绕组分布建立电机槽矩阵Mslot,用于描述任意电机槽内导体分布情况,该矩阵的维度为m×Qs,其中m为电机的相数,Qs为电机槽数。其建立的具体规则为:当电机某个槽内的导体均属于同一相时,槽矩阵中相应的行和列的元素记为1;当电机某个槽内的导体不属于同一相时,槽矩阵中相应的行和列的元素分别记为0.5。以六相36槽6极整数槽双层短距绕组为例,其一个极距内的绕组分布如图3所示:
其槽矩阵可表示为:
接下来,建立电机的槽磁动势函数,其定义为当槽内导体匝数为一匝,通入大小为1A的电流时,导体产生的磁动势沿定子圆周的分布,可以用公式表示为:
其中,Nq(θs)为槽磁动势函数,Δθ为槽宽,θq分别为第q个槽中心位置,其计算公式为:θq=(q-1)·2π/Qs,θs为任意时刻转子位置,Qs为电机槽数。
结合相电流、槽矩阵和槽磁动势函数,根据公式(3)则可以计算出电机在任意时刻与电流下的电枢磁动势。
farm=Iph,t·Mslot·Nq(θs) (3)
(4)通过静态有限元仿真计算电机的复数气隙磁导函数,具体计算流程如图4所示:
第一步,建立无槽定子、无槽转子有限元仿真模型,并施加幅值大小为1的恒定虚拟磁动势,仿真计算定、转子均无槽时的气隙磁密,其径向和切向分量分别为Br,slotless和Bt,slotless;第二步,建立有槽定子-无槽转子有限元模型,并施加幅值大小为1的恒定虚拟磁动势,仿真计算定子有槽、转子无槽时的气隙磁密,其径向和切向分量分别为和第三步,建立无槽定子-有槽转子有限元模型,并施加幅值大小为1的恒定虚拟磁动势,仿真计算定子无槽、转子有槽时的气隙磁密,其径向和切向分量分别为和将以上计算得到的三组数据代入公式(4)和(5),可以计算得到定子和转子的复数气隙磁导。
(5)将电枢磁动势、永磁体磁动势、定转子气隙磁导代入公式(6)中,可以计算出电机的气隙磁密。
(5)在气隙磁密的基础上,基于公式(7)所示麦克斯韦应力张量法,可以计算得到电机的气隙电磁力密度;由于电机振动和噪声是由于定子齿上所受电磁力引起,因此可以将气隙电磁力密度在一个齿距内进行积分来计算每个定子齿上所受的电磁力,具体计算方法如公式(8)所示:
其中,σr,σt分别表示径向、切向的气隙电磁力密度,Br,Bt分别表示气隙磁密的径向、切向分量,μ0为真空磁导率,θz表示第z个齿中心位置,Fr和Ft分别为定子齿上径向和切向电磁力,Lstk为定子铁芯长度,Rsi为定子内径。
(6)为了准确描述电机的机械结构特性,本发明通过建立三维结构有限元模型计算出电机的振动传递函数,具体计算流程如图5所示:
由于电机定子的叠压结构和导体在槽内分布的不均匀性,导致电机铁芯和绕组的材料参数呈现明显的正交异性。因此,本发明首先基于材料力学理论计算定子铁芯及绕组的等效材料参数;然后建立电机三维结构有限元模型,在电机定子齿上施加幅值为1,沿定子圆周呈周期分布的径向力和切向力,即Funit,z(f)=1,其中,Funit,z(f)表示施加的单位激振力,f为力的频率。然后通过有限元仿真计算得到电机不同齿在单位激振力作用下的振动传递函数其中,为传递函数的幅值,φH,z为传递函数相位。
(7)根据上述步骤获得的电机定子齿电磁力和振动传递函数,即可计算出电机在任意电磁力激励下的振动速度和振动加速度。具体计算流程如图6所示。首先,将定子齿电磁力进行傅里叶分解,得到定子齿电磁力的频谱分布,如公式(9)所示。其中,为电磁力谐波的幅值,φF,z为电磁力谐波的相位。
将计算得到的各个定子齿上的电磁力谐波和对应齿的振动传递函数相乘,并将所有齿的计算结果叠加即可得到电机的总振动变形,最后基于位移、速度和加速度之间的关系计算出电机表面的振动速度和振动加速度,如公式(10)-(11)所示。
其中,y,v,a分别为振动位移、振动速度和振动加速度,和为振动位移的一阶和二阶导数,为振动速度的一阶导数,电机不同齿在单位激振力作用下的振动传递函数,为传递函数的幅值,φH,z为传递函数相位;Fr(z,f)和Ft(z,f)分别为定子齿上径向和切向电磁力,分别表示径向、切向电磁力的幅值,φr,z,φt,z分别表示径向、切向电磁力相位,f为频率,t为时间。
(8)最后基于电机表面振动速度计算电机辐射的声压级和声功率,首先基于公式(12)计算不同频率的激振力下的电机表面声辐射效率。
其中,ω为声辐射的频率,c0为空气中声音的传播速度,k0=ω/c0为空气中声音的波长,d为定子机壳外径,m为轴向模态阶次,σm为m阶模态的声辐射效率,Jm,Jm+1分别为m阶和m+1阶第一类Bessel函数,Ym,Ym+1分别为m阶和m+1阶第二类Bessel函数。
然后将振动速度和声辐射效率代入公式(13)中,可以计算得到电机的声辐射功率,其中,Wm(ω)表示m阶模态辐射声功率,vmω为空间阶次为m、频率为ω的振动速度幅值,ρ0为空气密度,Sc为辐射面积。
则声功率级为:
其中,W0为基准声功率,W0=10-2W。
基于上述实施例的内容,作为一种可选的实施例,本发明实施例中提供的电机电磁力及电磁振动噪声分析系统包括电枢磁动势分析模块:用于确定电机的运行工况,根据其运行工况获得电机任意时刻的相电流,根据电机槽极配合、绕组分布确定电机的槽磁动势函数和槽矩阵,并结合所述相电流获得电机任意时刻的电枢磁动势。
基于上述实施例的内容,作为一种可选的实施例,本发明实施例中提供的电机电磁力及电磁振动噪声分析系统包括复数气隙磁导函数计算模块:用于建立所述电机的静态有限元模型,并分析获得该电机永磁体磁动势以及该电机定子和转子的复数气隙磁导函数。
基于上述实施例的内容,作为一种可选的实施例,本发明实施例中提供的电机电磁力及电磁振动噪声分析系统包括电磁力计算模块:用于根据所述电枢磁动势、永磁体磁动势及定子和转子的复数气隙磁导函数获得该电机的气隙磁密,并基于该气隙磁密获得所述电机每个定子齿上所受的电磁力。基于上述实施例的内容,作为一种可选的实施例,本发明实施例中提供的电机电磁力及电磁振动噪声分析系统包括振动传递函数分析模块:用于建立所述电机的三维结构有限元模型,并分析获得电机的振动传递函数。
基于上述实施例的内容,作为一种可选的实施例,本发明实施例中提供的电机电磁力及电磁振动噪声分析系统包括电机振动计算模块:用于根据所述电机定子齿电磁力和振动传递函数,计算获得所述电机在任意电磁力激励下的振动速度和振动加速度。
基于上述装置实施例的内容,作为一种可选的实施例,本发明实施例中提供的电机电磁振动噪声分析系统包括电机辐射的声压级和声功率计算模块:用于根据所述振动速度和振动加速度获得电机辐射的声压级和声功率。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来。
附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机软件产品的可能实现的体系架构、功能和操作。基于这种认识,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,有时也可以按相反的顺序执行,这依所涉及的功能而定。
在本专利中,术语"包括"、"包含"或者其任何其它变体意在涵盖非排它性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其它要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句"包括……"限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (8)
1.一种电机电磁力及电磁振动噪声分析方法,其特征在于,包括如下步骤:
S100:确定电机的运行工况,根据其运行工况获得电机任意时刻的相电流,根据电机槽极配合、绕组分布确定电机的槽磁动势函数和槽矩阵,并结合所述相电流获得电机任意时刻的电枢磁动势;
S200:建立所述电机的静态有限元模型,并分析获得该电机永磁体磁动势以及该电机定子和转子的复数气隙磁导函数;
S300:根据所述电枢磁动势、永磁体磁动势及定子和转子的复数气隙磁导函数获得该电机的气隙磁密,并基于该气隙磁密获得所述电机每个定子齿上所受的电磁力;
S400:建立所述电机的三维结构有限元模型,并分析获得电机的振动传递函数;
S500:根据所述电机定子齿电磁力和振动传递函数,计算获得所述电机在任意电磁力激励下的振动速度和振动加速度,基于所述振动速度和振动加速度获得电机辐射的声压级和声功率。
2.根据权利要求1所述的一种电机电磁力及电磁振动噪声分析方法,其特征在于,步骤S100具体包括:
S101:确定电机的槽极配合与绕组分布;
S102:根据所述槽极配合与绕组分布,建立电机槽矩阵;
S103:建立所述电机的槽磁动势函数;
S104:根据电机的运行工况,确定电机的相电流瞬时值;
S105:结合所述相电流瞬时值、电机槽矩阵以及槽磁动势函数获得电机在任意时刻与电流下的电枢磁动势。
4.根据权利要求1所述的一种电机电磁力及电磁振动噪声分析方法,其特征在于,步骤S200具体包括:
S201~S203:建立无槽定子和无槽转子有限元仿真模型,并施加幅值大小为1的恒定虚拟磁动势,仿真计算定子及转子均无槽时的气隙磁密;
S204~S205:建立有槽定子-无槽转子有限元模型,并施加幅值大小为1的恒定虚拟磁动势,仿真计算定子有槽、转子无槽时的气隙磁密;
S207~S208:建立无槽定子-有槽转子有限元模型,并施加幅值大小为1的恒定虚拟磁动势,仿真计算定子无槽、转子有槽时的气隙磁密;
S206、S209:根据所述定子及转子均无槽时的气隙磁密、定子有槽、转子无槽时的气隙磁密以及定子无槽、转子有槽时的气隙磁密,得到定子和转子的复数气隙磁导函数。
6.根据权利要求1所述的一种电机电磁力及电磁振动噪声分析方法,其特征在于,步骤S400具体包括:
S401:确定电机定子铁芯及绕组的等效材料参数;
S402:建立所述电机的三维结构有限元模型;
S403:在电机定子齿上施加幅值为1,沿定子圆周呈周期分布的径向力和切向力;
S404:通过有限元仿真计算得到电机不同齿在单位激振力作用下的振动传递函数。
7.根据权利要求1所述的一种电机电磁力及电磁振动噪声分析方法,其特征在于,步骤S500具体包括:
S502:将定子齿电磁力进行傅里叶分解,得到定子齿电磁力的频谱分布;
S503:将计算得到的各个定子齿上的电磁力谐波和对应齿的振动传递函数求积,将所有齿的计算结果叠加即可得到电机的总振动变形;
S505:基于位移、速度和加速度之间的关系获得电机表面的振动速度和振动加速度。
8.一种电机电磁力及电磁振动噪声分析系统,其特征在于,包括:
电枢磁动势分析模块:用于确定电机的运行工况,根据其运行工况获得电机任意时刻的相电流,根据电机槽极配合、绕组分布确定电机的槽磁动势函数和槽矩阵,并结合所述相电流获得电机任意时刻的电枢磁动势;
复数气隙磁导函数计算模块:用于建立所述电机的静态有限元模型,并分析获得该电机永磁体磁动势以及该电机定子和转子的复数气隙磁导函数;
电磁力计算模块:用于根据所述电枢磁动势、永磁体磁动势及定子和转子的复数气隙磁导函数获得该电机的气隙磁密,并基于该气隙磁密获得所述电机每个定子齿上所受的电磁力;
振动传递函数分析模块:用于建立所述电机的三维结构有限元模型,并分析获得电机的振动传递函数;
电机振动计算模块:用于根据所述电机定子齿电磁力和振动传递函数,计算获得所述电机在任意电磁力激励下的振动速度和振动加速度;
电机辐射的声压级和声功率计算模块:用于根据所述振动速度和振动加速度获得电机辐射的声压级和声功率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011426440.5A CN112906145B (zh) | 2020-12-09 | 2020-12-09 | 电机电磁力与电磁振动噪声分析方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011426440.5A CN112906145B (zh) | 2020-12-09 | 2020-12-09 | 电机电磁力与电磁振动噪声分析方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112906145A true CN112906145A (zh) | 2021-06-04 |
CN112906145B CN112906145B (zh) | 2024-05-14 |
Family
ID=76111530
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011426440.5A Active CN112906145B (zh) | 2020-12-09 | 2020-12-09 | 电机电磁力与电磁振动噪声分析方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112906145B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113364172A (zh) * | 2021-06-17 | 2021-09-07 | 江苏大学 | 一种减小内置式永磁电机振动噪声的转子辅助槽优化设计方法 |
CN113505495A (zh) * | 2021-07-29 | 2021-10-15 | 中国第一汽车股份有限公司 | 一种永磁同步电机定子电磁力分析方法 |
CN113765324A (zh) * | 2021-08-30 | 2021-12-07 | 江苏大学 | 针对新型不对称模块化永磁辅助同步磁阻电机及其振噪优化方法 |
CN114172435A (zh) * | 2021-11-30 | 2022-03-11 | 重庆长安汽车股份有限公司 | 异步电动机运行噪声与振动数据的后处理方法、系统及计算机可读存储介质 |
CN116205091A (zh) * | 2023-05-06 | 2023-06-02 | 广汽埃安新能源汽车股份有限公司 | 一种电机噪声的优化方法、装置、电子设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109039215A (zh) * | 2018-09-05 | 2018-12-18 | 上海电力学院 | 逆变器谐波对车用永磁同步电机振动噪声影响分析方法 |
CN109214125A (zh) * | 2018-10-22 | 2019-01-15 | 华中科技大学 | 一种电机电磁振动噪声的计算方法 |
EP3660522A1 (en) * | 2018-11-30 | 2020-06-03 | Zhejiang University | Method for evaluating electromagnetic performace of permanent-magnet machines |
CN111697892A (zh) * | 2020-06-24 | 2020-09-22 | 华中科技大学 | 一种短初级直线电机绕组磁动势分析方法 |
-
2020
- 2020-12-09 CN CN202011426440.5A patent/CN112906145B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109039215A (zh) * | 2018-09-05 | 2018-12-18 | 上海电力学院 | 逆变器谐波对车用永磁同步电机振动噪声影响分析方法 |
CN109214125A (zh) * | 2018-10-22 | 2019-01-15 | 华中科技大学 | 一种电机电磁振动噪声的计算方法 |
EP3660522A1 (en) * | 2018-11-30 | 2020-06-03 | Zhejiang University | Method for evaluating electromagnetic performace of permanent-magnet machines |
CN111697892A (zh) * | 2020-06-24 | 2020-09-22 | 华中科技大学 | 一种短初级直线电机绕组磁动势分析方法 |
Non-Patent Citations (4)
Title |
---|
YANG LU 等: "Electromagnetic Force and Vibration Analysis of Permanent-Magnet-Assisted Synchronous Reluctance Machines", IEEE TRANSACTIONS ON INDUSTRY APPLICATIONS, vol. 54, no. 5, 31 October 2018 (2018-10-31) * |
YANG LU 等: "Electromagnetic Force and Vibration Study on Axial Flux Permanent Magnet Synchronous Machines With Dual Three-Phase Windings", IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS, vol. 67, no. 1, 31 January 2020 (2020-01-31), XP011743681, DOI: 10.1109/TIE.2018.2890494 * |
林福 等: "考虑电流谐波的永磁同步电机电磁振动和噪声半解析模型", 电工技术学报, no. 09, 10 May 2017 (2017-05-10) * |
高鹏 等: "电动汽车用永磁同步电机电磁振动噪声分析及优化", 微电机, vol. 52, no. 12, 31 December 2019 (2019-12-31) * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113364172A (zh) * | 2021-06-17 | 2021-09-07 | 江苏大学 | 一种减小内置式永磁电机振动噪声的转子辅助槽优化设计方法 |
CN113364172B (zh) * | 2021-06-17 | 2023-02-17 | 江苏大学 | 一种减小内置式永磁电机噪声的转子辅助槽优化设计方法 |
CN113505495A (zh) * | 2021-07-29 | 2021-10-15 | 中国第一汽车股份有限公司 | 一种永磁同步电机定子电磁力分析方法 |
CN113765324A (zh) * | 2021-08-30 | 2021-12-07 | 江苏大学 | 针对新型不对称模块化永磁辅助同步磁阻电机及其振噪优化方法 |
CN114172435A (zh) * | 2021-11-30 | 2022-03-11 | 重庆长安汽车股份有限公司 | 异步电动机运行噪声与振动数据的后处理方法、系统及计算机可读存储介质 |
CN114172435B (zh) * | 2021-11-30 | 2023-06-02 | 重庆长安汽车股份有限公司 | 异步电动机运行噪声与振动数据的后处理方法、系统及计算机可读存储介质 |
CN116205091A (zh) * | 2023-05-06 | 2023-06-02 | 广汽埃安新能源汽车股份有限公司 | 一种电机噪声的优化方法、装置、电子设备及存储介质 |
CN116205091B (zh) * | 2023-05-06 | 2023-07-28 | 广汽埃安新能源汽车股份有限公司 | 一种电机噪声的优化方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112906145B (zh) | 2024-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112906145A (zh) | 电机电磁力与电磁振动噪声分析方法及系统 | |
Lee et al. | Particle swarm optimization algorithm with intelligent particle number control for optimal design of electric machines | |
CN109214125B (zh) | 一种电机电磁振动噪声的计算方法 | |
Lin et al. | Noise prediction and sound quality analysis of variable-speed permanent magnet synchronous motor | |
Le Besnerais et al. | Optimal slot numbers for magnetic noise reduction in variable-speed induction motors | |
CN113258696B (zh) | 一种用于降低分数槽集中绕组永磁电机电磁振动的方法 | |
Qian et al. | Interior noise and vibration prediction of permanent magnet synchronous motor | |
Frias et al. | Rotor and stator shape optimization of a synchronous machine to reduce iron losses and acoustic noise | |
CN106934162A (zh) | 一种基于磁路法与有限元法的电机噪声优化方法及装置 | |
CN113765251B (zh) | 一种确定削弱整数槽永磁电机振动的转子分段方式的方法 | |
CN104063556A (zh) | 一种轴向永磁电机等效磁路模型的建模方法 | |
CN108319768B (zh) | 一种基于计算机的永磁电机电枢反应磁场预测方法 | |
CN111199124A (zh) | 表贴式永磁同步电机空载电磁激振力波的计算方法 | |
CN111241735A (zh) | 内置式永磁同步电机负载电磁激振力波的计算方法 | |
Dupont et al. | Multiphysics modelling to simulate the noise of an automotive electric motor | |
Lee et al. | High-speed, high-power motor design for a four-legged robot actuator optimized using the weighted sum and response surface methods | |
CN112152400B (zh) | 永磁体涡流损耗的优化方法 | |
CN113452209B (zh) | 一种多相笼型感应电机磁场计算方法及系统 | |
CN105302999A (zh) | 一种发现汽轮发电机内负序分量对转矩与谐波影响规律的方法 | |
Martis et al. | Vibroacoustic behaviour analysis of a permanent magnet synchronous machine for automotive applications | |
Sutthiphornsombat et al. | Mitigation of acoustic noise and vibration in permanent magnet synchronous machines drive using field reconstruction method | |
Huang et al. | Tangential Electromagnetic Force Array on the Vibration and Noise of Electric Axle for New Energy Vehicle | |
Wu et al. | Numerical prediction and analysis of electromagnetic vibration and noise of claw pole alternator | |
Wu et al. | A review of the research progress of motor vibration and noise | |
Varghese | Electromagnetic noise from permanent magnet motors |
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 |