CN110555254B - 一种磁极分组降低永磁电机面外振动不稳定的方法 - Google Patents

一种磁极分组降低永磁电机面外振动不稳定的方法 Download PDF

Info

Publication number
CN110555254B
CN110555254B CN201910796164.2A CN201910796164A CN110555254B CN 110555254 B CN110555254 B CN 110555254B CN 201910796164 A CN201910796164 A CN 201910796164A CN 110555254 B CN110555254 B CN 110555254B
Authority
CN
China
Prior art keywords
vibration
grouping
group
instability
equation
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.)
Active
Application number
CN201910796164.2A
Other languages
English (en)
Other versions
CN110555254A (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201910796164.2A priority Critical patent/CN110555254B/zh
Publication of CN110555254A publication Critical patent/CN110555254A/zh
Application granted granted Critical
Publication of CN110555254B publication Critical patent/CN110555254B/zh
Active 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
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Reciprocating, Oscillating Or Vibrating Motors (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

本发明公开一种磁极分组降低永磁电机面外振动不稳定的方法,包括以下步骤:在随动坐标系下建立动力学模型;借助三角函数的运算性质,判断振动波数与分组数之间的关系,分类求解特征方程;根据基本振动理论,通过调节组内夹角,降低面外振动的不稳定。本发明解决了现有转子面外振动存在的不稳定现象,通过分组的形式,调节组内夹角的关系,从而降低不稳定现象,使得计算结果更好地满足工程需求。

Description

一种磁极分组降低永磁电机面外振动不稳定的方法
技术领域
本发明涉及抑制振动的领域,尤其针对一种磁极分组降低永磁电机面外振动不稳定的技术。
背景技术
电机转子是一种典型的环状周期结构,永磁电机又以其体积小、结构紧凑等特点,广泛应用于各工程领域,例如:航空航天、数控机床、汽车电气、家用电器、计算机外围设备等。在工程实际中,该类结构将出现显著的振动、噪声以及稳定性问题。因此特别需要提出一种降低系统不稳定的分析技术。
文献(Jamali Arand S,Ardebili M.Cogging torque reduction in axial-fluxpermanent magnet wind generators with yokeless and segmented armature byradially segmented and peripherally shifted magnet pieces.Renewable Energy,2016,99:95-106)通过三维有限元方法,提出了一种通过使永磁体偏斜,降低齿槽转矩。但是,由于永磁体偏斜的角度在加工安装中,对精度要求过高,所以该文献提出的方法较难实现。
文献(Jianfeng H,Shanming W,Yuguang S.An effective method with copperring for vibration reduction in permanent magnet brush DC motors.IEEETransactions on Magnetics,2018,54:8108105)通过在永磁体上增加铜圆环的方式,降低了电机振动的影响。但是,这种通过附加质量的方法,势必会对电机工作的效率造成影响。
此外,现有技术还通常采用数值方法进行预测和有限元进行模拟降低动力稳定性的方法,该方法的计算效率较低,且不能揭示普适规律。
发明内容
本发明的目的是为了克服现有技术中的不足,提供一种磁极分组降低永磁电机面外振动不稳定的方法,本发明解决了现有降低永磁电机转子振动不稳定技术的不足,从而提出一种便于实现的降低不稳定的技术,更好地满足工程需求。
本发明的目的是通过以下技术方案实现的:
一种磁极分组降低永磁电机面外振动不稳定的方法,包括以下步骤:
在随动坐标系下建立动力学模型;
借助三角函数的运算性质,判断振动波数与分组数之间的关系,分类求解动力学方程的特征方程;
根据基本振动理论,通过调节组内夹角,降低面外振动的不稳定,根据所获得的特征方程,求得特征值的虚实部,即可根据基本的振动理论预测转子面外振动的不稳定现象,然后通过调节组内夹角降低系统的不稳定现象;
考虑磁极分组降低永磁电机面外振动不稳定的方法的具体步骤如下:
借助随动坐标系,建立系统的动力学方程;
借助三角函数的运算性质,判断振动波数与分组数之间的关系,分类求解特征方程;
动力学方程的响应为公式(1),在不同波数与分组数下分别为公式(2)和公式(3),将式
Figure GDA0004168620220000021
分别代入式(2)和式(3)中,求解特征值虚实部,实现降低不稳定的方法;WRe和WIm分别表示实部和虚部的振幅,λ为特征值,β为相位;
Figure GDA0004168620220000022
式中,w为面外振动位移,
Figure GDA0004168620220000023
为位置角,Ω为转子转速,kt为离心刚度算子,krp和krs分别表示动支撑刚度和静支撑刚度算子,kp表示磁刚度算子;
当2n/N1=int时,特征方程为
Figure GDA0004168620220000024
当2n/N1≠int时,特征方程为
Figure GDA0004168620220000025
式中,n为振动波数,N1为分组数,int为整数,M为质量矩阵,G为陀螺矩阵,Kc和Ku分别为不受组合关系影响和受组合关系影响的刚度矩阵。
与现有技术相比,本发明的技术方案所带来的有益效果是:
1、本发明首先在随动坐标系下建立动力学模型;借助三角函数的运算性质,判断振动波数与分组数之间的关系,分类求解特征方程;根据基本振动理论,通过调节组内夹角,降低面外振动的不稳定;
2、本发明采用解析法给出了转子面外振动的特征值,根据特征值判断系统的动力稳定性,从而提出降低不稳定的方法;
3、与现有技术相比,本发明具有高效、准确和普适特征,根据该技术可揭示组内夹角对于降低不稳定的影响,从而指导旋转对称机械的动力学设计,提高运行稳定性和可靠性。
附图说明
图1为本发明提供的转子面外振动磁极分组及坐标系的示意图;
图2(a)为根据本发明提供的方法获得的组内磁极个数为2,波数为2时,转速与组内夹角对于不稳定的影响;
图2(b)为根据本发明提供的方法获得的组内磁极个数为2,波数为3时,转速与组内夹角对于不稳定的影响;
图3(a)为根据本发明提供的方法获得的组内磁极个数为2,波数为2时,不同组内夹角在影响较大转速作用下,特征值实部的值;
图3(b)为根据本发明提供的方法获得的组内磁极个数为2,波数为3时,不同组内夹角在影响较大转速作用下,特征值实部的值;
图4(a)为根据本发明提供的方法获得的组内磁极个数为2,波数为2时,不同组内夹角在影响较大转速作用下,特征值虚部的值;
图4(b)为根据本发明提供的方法获得的组内磁极个数为2,波数为3时,不同组内夹角在影响较大转速作用下,特征值虚部的值;
图5(a)为根据本发明提供的方法获得的组内磁极个数为3,波数为2时,转速与组内夹角对于不稳定的影响;
图5(b)为根据本发明提供的方法获得的组内磁极个数为3,波数为3时,转速与组内夹角对于不稳定的影响;
图6(a)为根据本发明提供的方法获得的组内磁极个数为3,波数为2时,不同组内夹角在影响较大转速作用下,特征值实部的值;
图6(b)为根据本发明提供的方法获得的组内磁极个数为3,波数为3时,不同组内夹角在影响较大转速作用下,特征值实部的值;
图7(a)为根据本发明提供的方法获得的组内磁极个数为3,波数为2时,不同组内夹角在影响较大转速作用下,特征值实部的值;
图7(b)为根据本发明提供的方法获得的组内磁极个数为3,波数为3时,不同组内夹角在影响较大转速作用下,特征值实部的值;
具体实施方式
以下结合附图和具体实施例对本发明作进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明实施例提出了一种磁极分组降低永磁电机面外振动不稳定的方法,本发明实施例首先借助随动坐标系建立环状周期结构的动力学模型,然后借助三角函数的运算性质,判断振动波数与分组数之间的关系,分类求解特征方程;根据基本振动理论,通过调节组内夹角,降低面外振动的不稳定。
本发明实施例可适用于航空航天、数控机床、汽车电气、家用电器、计算机外围设备等需要应用永磁电机转子的领域。
本发明实施例的方法方案如下:一种磁极分组降低永磁电机面外振动不稳定的方法,考虑了转子转动、磁极分组及轴向(w)弹性振动产生的动力稳定性降低方法。
永磁电机转子结构由薄圆环、辐板支撑和分组分布的磁极组成;降低不稳定的方法的基本特征在于:采用随动坐标系,调整组内夹角,降低转子的不稳定,具体步骤为:
(A1)图1为永磁电机转子面外振动示意图,
Figure GDA0004168620220000041
为随动坐标系。在外环上均布N1组磁极,每组由N2个磁极组成,
Figure GDA0004168620220000042
为第i组第j个磁极的初始位置角,其中
Figure GDA0004168620220000043
假设组内均布,α为组内相邻两个磁极上边缘之间的夹角。设第一组第一个磁极的下边缘位于极轴上,且磁极夹角为γ,则第i组第j个磁极的下边缘位置可以写作
Figure GDA0004168620220000044
上边缘为
Figure GDA0004168620220000045
用H(.)函数,即阶跃函数去描述其位置。
借助随动坐标系,根据Hamilton原理建立所述环状周期结构的动力学模型:
Figure GDA0004168620220000046
式中,w为面外振动位移,
Figure GDA0004168620220000047
为位置角,Ω为转子转速,kt为离心刚度算子,krp和krs分别表示动支撑刚度和静支撑刚度算子,kc表示空间旋转算子,kp表示磁刚度算子。
(A2)进一步地,根借助三角函数的运算性质,判断振动波数与分组数之间的组合关系,然后分类计算面外振动的特征值,特征方程分别为:
当2n/N1=int时,特征方程为
Figure GDA0004168620220000048
当2n/N1≠int时,特征方程为
Figure GDA0004168620220000049
式中,n为振动波数,N1为分组数,int为整数,M为质量矩阵,G为陀螺矩阵,Kc和Ku分别为不受组合关系影响和受组合关系影响的刚度矩阵。式中,
Figure GDA0004168620220000051
式中,轴向振动w为复函数,其实部用wRe表示,虚部用wIm表示。设刚度矩阵的形式为
Figure GDA0004168620220000052
其中,
Kc11=Kc22=[krp+n2(kt-1)]Ω2+n4+krs-kapγ (6)
Kc12=Kc21=0 (7)
Figure GDA0004168620220000053
Figure GDA0004168620220000054
式中,
Figure GDA0004168620220000055
式中,
Figure GDA0004168620220000056
Figure GDA0004168620220000057
Figure GDA0004168620220000058
Figure GDA0004168620220000059
Figure GDA00041686202200000510
Figure GDA00041686202200000511
式中,τ(0<τ<1)为辐板切向应力调节系数,ρ0为转子密度,E、I、R、h和b分别为外环的弹性模量、截面惯性矩、中性圆半径、厚度和宽度,hp、Rb和Ra分别为辐板的厚度、外径和内径。式中,
Figure GDA0004168620220000061
Figure GDA0004168620220000062
Figure GDA0004168620220000063
Figure GDA0004168620220000064
式中,Br和hm分别为永磁体的剩磁和充磁厚度,μ0为真空磁导率,δ为定转子之间的气隙长度。
(A3)求解永磁电机转子面外振动的特征值,设式(2)和(3)特征解的形式为
Figure GDA0004168620220000065
式中,WRe和WIm分别表示实部和虚部的振幅,λ为特征值,β为相位。
将式(14)式分别代入式(2)和(3)中,可以得到特征值求解方程
Figure GDA0004168620220000066
式中,
Figure GDA0004168620220000067
为了进一步解析分析,将特征值写作
λ=λRe+iλIm (17)
将式(16)代入式(15)中,即可得到不同组合下面外振动特征值实部和虚部。根据基本振动理论,即可预测系统的稳定性。
针对上述振动方程的特征,本发明实施例提出了一种磁极分组降低永磁电机面外振动不稳定的方法,并据此降低系统的动力稳定性,具体过程如下:
(B1)在随动坐标系下建立动力学模型;
(B2)借助三角函数的运算性质,判断振动波数与分组数之间的关系,分类求解特征方程;
(B3)根据基本振动理论,通过调节组内夹角,降低面外振动的不稳定。
考虑磁极分组降低永磁电机面外振动不稳定的方法的具体步骤如下:
(C1)借助随动坐标系,建立系统的动力学方程;
(C2)借助三角函数的运算性质,判断振动波数与分组数之间的关系,分类求解特征方程。
假设步骤(C1)中的动力学方程的响应为上述公式(1),在不同波数与分组数下分别为公式(2)和公式(3),将式(15)分别代入式(2)和式(3)中,求解特征值虚实部,提出降低不稳定的方法。
(C3)以表1中的环状周期结构参数为例,计算特征值。
表1永磁电机转子基本参数
Figure GDA0004168620220000071
(C4)根据步骤(C3)所求的特征值,可得组内磁极个数为2,波数为2时,转速与组内夹角对于不稳定的影响,具体结果如附图2(a)所示。
(C5)根据步骤(C3)所求的特征值,可得组内磁极个数为2,波数为3时,转速与组内夹角对于不稳定的影响,具体结果如附图2(b)所示。
(C6)根据步骤(C3)所求的特征值,可得组内磁极个数为2,波数为2时,不同组内夹角在影响较大转速作用下,特征值实部的值,具体结果如附图3(a)所示。
(C7)根据步骤(C3)所求的特征值,可得组内磁极个数为2,波数为3时,不同组内夹角在影响较大转速作用下,特征值实部的值,具体结果如附图3(b)所示。
(C8)根据步骤(C3)所求的特征值,可得组内磁极个数为2,波数为2时,不同组内夹角在影响较大转速作用下,特征值虚部的值,具体结果如附图4(a)所示。
(C9)根据步骤(C3)所求的特征值,可得组内磁极个数为2,波数为3时,不同组内夹角在影响较大转速作用下,特征值虚部的值,具体结果如附图4(b)所示。
(C10)根据步骤(C3)所求的特征值,可得组内磁极个数为3,波数为2时,转速与组内夹角对于不稳定的影响,具体结果如附图5(a)所示。
(C11)根据步骤(C3)所求的特征值,可得组内磁极个数为3,波数为3时,转速与组内夹角对于不稳定的影响,具体结果如附图5(b)所示。
(C12)根据步骤(C3)所求的特征值,可得组内磁极个数为3,波数为2时,不同组内夹角在影响较大转速作用下,特征值实部的值,具体结果如附图6(a)所示。
(C13)根据步骤(C3)所求的特征值,可得组内磁极个数为3,波数为3时,不同组内夹角在影响较大转速作用下,特征值实部的值,具体结果如附图6(b)所示。
(C14)根据步骤(C3)所求的特征值,可得组内磁极个数为3,波数为2时,不同组内夹角在影响较大转速作用下,特征值虚部的值,具体结果如附图7(a)所示。
(C15)根据步骤(C3)所求的特征值,可得组内磁极个数为3,波数为3时,不同组内夹角在影响较大转速作用下,特征值虚部的值,具体结果如附图7(b)所示。
综上所述,本发明提供了一种磁极分组降低永磁电机面外振动不稳定的方法。该方法运用随动坐标系,采用解析方法得到系统的特征值,提高了准确性、计算效率及普适性,更好地满足了工程实际的需求。
本领域方法人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
本发明并不限于上文描述的实施方式。以上对具体实施方式的描述旨在描述和说明本发明的方法方案,上述的具体实施方式仅仅是示意性的,并不是限制性的。在不脱离本发明宗旨和权利要求所保护的范围情况下,本领域的普通方法人员在本发明的启示下还可做出很多形式的具体变换,这些均属于本发明的保护范围之内。

Claims (1)

1.一种磁极分组降低永磁电机面外振动不稳定的方法,其特征在于,包括以下步骤:
在随动坐标系下建立动力学模型;
借助三角函数的运算性质,判断振动波数与分组数之间的关系,分类求解动力学方程的特征方程;
根据基本振动理论,通过调节组内夹角,降低面外振动的不稳定,根据所获得的特征方程,求得特征值的虚实部,即可根据基本的振动理论预测转子面外振动的不稳定现象,然后通过调节组内夹角降低系统的不稳定现象;
考虑磁极分组降低永磁电机面外振动不稳定的方法的具体步骤如下:
借助随动坐标系,建立系统的动力学方程;
借助三角函数的运算性质,判断振动波数与分组数之间的关系,分类求解特征方程;
动力学方程的响应为公式(1),在不同波数与分组数下分别为公式(2)和公式(3),将式
Figure FDA0004141728200000011
分别代入式(2)和式(3)中,求解特征值虚实部,实现降低不稳定的方法;WRe和WIm分别表示实部和虚部的振幅,λ为特征值,β为相位;
Figure FDA0004141728200000012
式中,w为面外振动位移,
Figure FDA0004141728200000013
为位置角,Ω为转子转速,kt为离心刚度算子,krp和krs分别表示动支撑刚度和静支撑刚度算子,kp表示磁刚度算子;
当2n/N1=int时,特征方程为
Figure FDA0004141728200000014
当2n/N1≠int时,特征方程为
Figure FDA0004141728200000015
式中,n为振动波数,N1为分组数,int为整数,M为质量矩阵,G为陀螺矩阵,Kc和Ku分别为不受组合关系影响和受组合关系影响的刚度矩阵。
CN201910796164.2A 2019-08-22 2019-08-22 一种磁极分组降低永磁电机面外振动不稳定的方法 Active CN110555254B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910796164.2A CN110555254B (zh) 2019-08-22 2019-08-22 一种磁极分组降低永磁电机面外振动不稳定的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910796164.2A CN110555254B (zh) 2019-08-22 2019-08-22 一种磁极分组降低永磁电机面外振动不稳定的方法

Publications (2)

Publication Number Publication Date
CN110555254A CN110555254A (zh) 2019-12-10
CN110555254B true CN110555254B (zh) 2023-06-06

Family

ID=68738051

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910796164.2A Active CN110555254B (zh) 2019-08-22 2019-08-22 一种磁极分组降低永磁电机面外振动不稳定的方法

Country Status (1)

Country Link
CN (1) CN110555254B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111563322B (zh) * 2020-04-24 2023-06-30 天津大学 一种消除永磁电机定子振动不稳定的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110046418A (zh) * 2019-04-09 2019-07-23 天津大学 一种永磁电机周期定子的振动特性分析方法
CN110119532A (zh) * 2019-04-09 2019-08-13 天津大学 一种旋转环状周期结构的弹性振动计算方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105181600B (zh) * 2015-09-02 2017-11-17 哈尔滨工业大学 一种大口径曲面光学元件表面微缺陷的检测与激光修复装置
CN106940735B (zh) * 2016-10-19 2021-05-14 天津大学 一种感应电机弹性振动稳定性预测方法
CN106547957B (zh) * 2016-10-19 2019-07-26 天津大学 一种旋转环状周期结构参激弹性振动分析方法
CN106737659B (zh) * 2016-11-15 2017-11-21 哈尔滨工程大学 一种水下无人航行器和机械手系统的手艇协调控制方法
CN208858869U (zh) * 2018-07-06 2019-05-14 华南理工大学 一种少齿差行星线齿轮减速器
CN109284569B (zh) * 2018-10-15 2023-04-07 天津大学 一种永磁电机周期定子的组合不稳定分析及预测方法
CN110011585B (zh) * 2019-05-20 2020-11-24 中国矿业大学 一种外部激励引起的永磁半直驱式传动轴系扭振控制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110046418A (zh) * 2019-04-09 2019-07-23 天津大学 一种永磁电机周期定子的振动特性分析方法
CN110119532A (zh) * 2019-04-09 2019-08-13 天津大学 一种旋转环状周期结构的弹性振动计算方法

Also Published As

Publication number Publication date
CN110555254A (zh) 2019-12-10

Similar Documents

Publication Publication Date Title
CN110046418B (zh) 一种永磁电机周期定子的振动特性分析方法
Zheng et al. Rotor balancing for magnetically levitated TMPs integrated with vibration self-sensing of magnetic bearings
Yang et al. Vibration reduction by magnetic slot wedge for rotating armature permanent magnet motors
CN110555254B (zh) 一种磁极分组降低永磁电机面外振动不稳定的方法
CN110580383A (zh) 一种分组拓扑径向受载圆环应力叠加的方法
CN109551026A (zh) 旋风铣床用电机直驱式磁悬浮铣刀盘
Li et al. Analysis, modeling, and verification of limited angle torque motors with irregular slot numbers for performance improvement
Cheng et al. Study on the design method of high speed permanent magnet synchronous machine
Sugimoto et al. A vibration reduction method of one-axis actively position regulated single-drive bearingless motor with repulsive passive magnetic bearings
CN110569574B (zh) 一种永磁电机正弦形磁极提高转子面外振动稳定性方法
Eshan et al. Rotodynamic study of a high speed switched reluctance generator
Werner et al. Rotor dynamic analysis of asynchronous machines including the finite-element-method for engineering low vibration motors
CN111563322B (zh) 一种消除永磁电机定子振动不稳定的方法
CN110569575B (zh) 一种永磁电机转子面外振动稳定性预测方法
Xing et al. Critical speed calculation and modal analysis of rotor for high-speed machine
Li et al. Modal analysis of electromagnetic resonance for multi-degree-of-freedom spherical motor
Sun et al. Electromagnetic Force Waves of Permanent Magnet Synchronous Motor Considering Eccentricity
Shi et al. Influence of static rotor eccentricity on suspension performance of bearingless doubly salient electromagnetic machine
CN110579706B (zh) 绕空间轴旋转永磁电机转子的面外振动稳定性预测方法
He et al. Rotor UMP & mechanical response in HSPMSM in typical running conditions
CN112436632B (zh) 压缩机及空调器
GB2574279A (en) Permanent magnet conical motors, stators and rotors thereof
CN114337016B (zh) 一种电机转子、电机及移动平台
Sugimoto et al. Low speed test in two-axis actively positioned bearingless machines with non-collocated structure for wind power applications
Filip et al. Multiphysics modeling of an permanent magnet synchronous machine

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