电机定子振动模态频率分裂特性分析.pdf
《电机定子振动模态频率分裂特性分析.pdf》由会员分享,可在线阅读,更多相关《电机定子振动模态频率分裂特性分析.pdf(10页珍藏版)》请在文库网上搜索。
1、第 36 卷第 4 期2023 年 8 月振 动 工 程 学 报Journal of Vibration EngineeringVol.36 No.4Aug.2023电机定子振动模态频率分裂特性分析徐港辉,祝长生(浙江大学电气工程学院,浙江 杭州 310027)摘要:有效抑制由电机径向电磁力激发的电机定子振动是实现电机减振降噪的一个重要途径,而对电机定子模态频率及模态振型的准确分析是抑制电机定子径向振动的基础。采用圆环的弹性力学解析模型作为电机定子振动的分析模型,对无约束状态下电机定子的模态进行分析,得到了电机定子径向振动模态频率和模态振型的解析解。以齿槽和底脚为典型附加结构,采用摄动法对电机
2、定子模态频率的分裂现象进行分析,总结了频率分裂与否以及分裂阶次的判定准则。通过 ANSYS 有限元软件验证了理论方法和计算的有效性。结果表明,所建立的二维圆环模型可以准确、高效地应用于电机定子模态特性的分析;附加结构的分布形式对定子频率分裂特性具有重要的影响。关键词:电机定子;模态频率;模态振型;圆环模型;频率分裂中图分类号:TM303.3;TB535 文献标志码:A 文章编号:1004-4523(2023)04-0953-10 DOI:10.16385/ki.issn.1004-4523.2023.04.008引言电机的振动和噪声是电机的一项重要技术指标,同时也是制约电机向高性能方向发展的一
3、个关键因素。由电机气隙磁场产生的径向电磁力所激发的电机定子的电磁振动是电机振动噪声的重要来源,在电机减振降噪领域引起了广泛关注1。为了减小电机的电磁振动,除了需要控制电磁力的幅值外,还需要避免电机定子的模态频率及其模态阶次与电磁力的频率及其阶次相等或者接近导致的定子结构的共振2。因此,电机定子的模态特性分析对电机的减振降噪具有重要意义,并且在电机的设计阶段就需要进行校核和规划。准确获得定子的模态频率和模态振型,一直是电机定子模态特性分析的重要内容。国外的研究在20世纪 30年代就已开始,国内学者对该问题的研究开始于 20世纪 80年代3,目前还有学者在继续开展相关的研究。电机定子的整体框架近似
4、为圆柱壳体结构46,但作为三维连续结构,圆柱壳体的振动特性分析相对复杂58。为此,在以往的诸多研究中,经常采用二维圆环模型作为电机定子振动的分析模型。邱家俊等9、于慎波等10分别通过解析计算和有限元仿真说明了两端无约束状态下圆柱壳体与圆环的径向模态频率十分接近,验证了采用圆环作为电机定子分析模型的合理性。基于圆环模型的机电类比法4,11具有公式简单、求解方便等优点,在电机定子模态频率计算中得到了广泛应用。李晓华等3考虑了绕组和浸漆结构对定子模态频率的影响,基于圆环模型计算了含绕组定子铁心的模态频率。王宏华等12考虑了定子凸极对磁轭质量和刚度的影响,采用等效圆环模型计算了开关磁阻电机的定子模态频
5、率。邢泽智等13考虑了定子齿槽结构对定子模态频率的影响,基于数据拟合法提出了一种定子铁心等效圆环尺寸计算公式,有效提升了定子径向模态频率的计算精度。然而以上研究主要围绕模态频率的准确计算开展,对模态振型的理论研究较少。并且机电类比法属于公式法,只能用于计算定子的模态频率,不能揭示定子的模态振型特性,方法的局限性较大。因此,为了深入分析定子的模态特性,需要建立相应的理论模型。在电机定子的模态试验中,文献 1416 均提到了“双共振峰”的现象。Girgis 等7指出定子结构中的不对称是导致这一现象产生的原因。陈永校等14提出了“广义齿槽”的概念来描述定子在尺寸、材料、约束等方面的不均匀性,从理论上
6、对“双共振峰”现象进行了解释。然而相关文献均没有给出“双共振峰”出现与否的判断标准。实际上,“双共振峰”是定子同一阶模态对应了两个接近但不相等的模态频率的一种现象(又称为“频率分裂”现象)。理想的对称结构往往具有重特征值(模态频率相等)的特性,此收稿日期:2021-11-25;修订日期:2022-01-30基金项目:国家自然科学基金资助项目(51975516);重点基础研究项目(2020-ZD-23200)。振 动 工 程 学 报第 36 卷时的系统又称为退化系统。当退化系统发生小参数变化后,原本的多重模态频率可能会分裂成多个数值不等的模态频率17,即所谓的频率分裂现象。在早期的研究中,胡海昌
7、17和陈塑寰18分别采用变分原理和摄动法实现了小参数变化下退化系统的模态特性分析。实际上,结构参数或者约束条件的变化均可能引起频率分裂现象。Wu等19采用摄动法和Galerkin 离散法对不同形式弹性支承下圆环的动力学特性进行了研究,得到了附加刚度对圆环模态频率分裂的影响规律。Bisegna 等20采用线性弹性理论和摄动法分析了结构缺陷引起的圆环频率分裂现象,得到了缺陷圆环模态频率和振型的闭环解,并通过 RayleighRitz 法 验 证 了 理 论 计 算 的 有 效 性。Wang 等21采用摄动法对含有附加质量和附加刚度的圆环的模态特性进行了分析,推导了圆筒形超声电机中模态频率分裂的一般
8、规律;并研究了分组对称形式的附加单元对旋转周期结构频率分裂特性的影响规律22。以上文献在分析结构的频率分裂特性时,广泛采用了摄动法的思想。在以上文献的基础上,为了进一步完善电机定子模态特性分析的理论,并探究电机定子振动中“双共振峰”现象的发生规律,考虑到机电类比法的局限性,本文首先采用圆环弹性力学解析模型对电机定子的模态特性进行分析,得到了各阶模态频率和模态振型的解析表达式。然后以齿槽和底脚为典型的附加结构,基于圆环解析模型,采用摄动法分析了这两类附加结构对定子模态频率的影响,总结了定子频率分裂与否的规律,为“双共振峰”现象提供了理论解释。最后,通过文献对比和有限元仿真对理论结果进行了验证。1
9、基于圆环解析模型的定子模态特性分析图 1为电机定子的圆环模型。其中 oz为建立在圆环中性线上的极坐标系,和 z分别为圆环的切向坐标和径向坐标,r 为圆环中性线的半径,h 为圆环的厚度,v 及 w 分别为圆环中性线上某点 P 在切向及径向方向上的位移。本文只考虑圆环的面内振动,圆环在面内的自由振动方程为23:h2vt2=()D+Kr2r42v2-Dr43w3+Kr2wh2wt2=Dr43v3-Kr2v-Dr44w4-Kr2w(1)式中 为圆环的密度,D=Eh312,K=Eh,E 为圆环的等效弹性模量。式(1)所示的圆环模型是关于时间和位置的偏微分方程组,直接求其解析解往往是困难的,通常借助变量分
10、离法进行求解。圆环的振动位移 v及 w 既是时间 t的函数,又是角度 的函数,根据变量分离法可设:v(),t=V()ejtw(),t=W()ejt(2)式中V()及W()为圆环的振型函数,与时间无关,表征了圆环不同角位置的相对振幅;为圆环面内振动的模态角频率。根据圆环在圆周方向上的周期性特征,假设圆环的各阶振型为如下两种三角函数形式:径向为余弦振型:V()=Vnsin(n)W()=Wncos(n)(3)径向为正弦振型:V()=-Vncos(n)W()=-Wnsin(n)(4)式中Vn及-Vn分别为切向 n 阶正弦及余弦振型的振幅;Wn及-Wn分别为径向 n 阶余弦及正弦振型的振幅;n为圆环模态
11、振型的周向阶次,表征了振型沿周向的波动特征。将式(2)(4)代入式(1),经过变量分离,将偏微分方程组简化为齐次常微分方程组。消去齐次方程组中相同的时间项和正余弦项,可以得到关于各阶振型振幅的方程组为:径向为余弦振型:P11P12P21P22VnWn=0(5)径向为正弦振型:Q11Q12Q21Q22 -Vn-Wn=0(6)式中P11=Q11=h2-(Dn2+Kr2n2)r4,P12=P21=-Q12=-Q21=-Dn3r4-Kn r2,P22=Q22=h2-Dn4r4-K r2。可以证明,式(5)与式(6)中系数矩阵的行列式是相等的。为了使各阶振型振幅有非零解,式(5)和(6)中系数矩阵的行列
12、式必等于零,从而得到圆环振图 1 电机定子的圆环模型示意图Fig.1 Schematic diagram of ring model of motor stator954第 4 期徐港辉,等:电机定子振动模态频率分裂特性分析动的模态频率方程为:4-12+2=0(7)式中1=(n2+1)(Dn2+Kr2)(hr4),2=n2(n2-1)2DK (2h2r6)。由式(7)可知,标准圆环径向的 n 阶正弦、余弦振型对应于相等的模态角频率。进一步求取式(7)的解析解(记为 n),有:2n=112-422(8)根据式(8)可知,每个 n值对应于两个大小不等的模态频率,记其中较小者为 n,s,较大者为 n
13、,b。当 n=0 或 1 时,由于2正比于 n2(n2-1)2,均有2=0,此时较小者 0,s=1,s=0,分别对应于圆环的两 个 刚 体 模 态;较 大 者 0,b=E/(r2)及 1,b2E/(r2)均不为 0。实际上,0,b对应于圆环的呼吸模态,即周向波数为 0的伸缩模态,此时圆环的振动形式在任意位置 都相同,因此又称为轴对称模态;而1,b对应于圆环周向波数为 1的伸缩模态23。当n 2时,结 合h2 r,对 式(8)进 行 化 简可得:2n,sn2(n2-1)2n2+1Eh212r4(9)2n,b(n2+1)Er2(10)可以验证,n,s对应于圆环 2阶及以上的弯曲模态;n,b对应于圆
14、环 2阶及以上的伸缩模态23。进一步地,根据式(9)和(10)可得:n,b/n,s12r2/(nh)2,由此可以判断,伸缩模态对应的频率通常要比弯曲模态的频率高很多。在求得定子的各阶模态频率之后,根据式(5)和(6)可以确定各阶振型的相对振幅。由式(5)和(6)的系数行列式为 0 可知,取其中任意一组方程均可确定各阶振幅之间的相对大小,有:VnWn=-Vn-Wn=-P12P11=Dn3+Kr2nh2nr4-()D+Kr2n2(11)联立式(11)与式(3),(4),即可确定圆环的各阶模态振型。根据式(11)可知,圆环同一阶次的正弦振型与余弦振型对应的振幅也是相等的。但是根据式(3)和(4)可知
15、两种振型之间存在/(2n)的相位差,并不完全重合。2电机定子模态的频率分裂现象频率分裂现象指的是同一阶模态(相同波动特征)对应于两个或两个以上不同频率的现象。以本文圆环模型为例,由式(7)可知,标准圆环径向的 n阶正弦、余弦模态对应的模态频率是相等的,此时两种模态对应的频率是重合的。然而,当圆环上存在附加结构时,由于圆环的对称性被改变,同一阶模态对应的频率可能会出现不相等的分裂现象。本节将在文献 21 的理论基础上,对电机定子模态频率的分裂现象做进一步的理论分析。2.1旋转对称结构引起的定子模态频率分裂现象在电机定子中,存在周向均布的齿槽、散热筋等旋转对称结构,这里以齿槽结构作为旋转对称结构的
16、典型代表,分析旋转对称结构引发模态频率分裂现象的规律。图 2 为含齿槽结构的定子圆环模型。基于式(1)所示的标准圆环面内振动方程,不考虑圆环的伸缩模态,有w=-v,代入式(1)并进行无量纲化处理,可以得到标准圆环面内振动的特征方程为:K()0v-()0nM()0v=0(12)式 中M(0)及K(0)为 微 分 算 子,且 有K()0=-(66+244+22),M()0=1-22;()0n=2n,n为标准圆环模态角频率;v 为无量纲化后的位移,无量纲化方法见文献 21。在式(12)的基础上,考虑引入齿槽结构带来的影响,根据摄动法可设微分算子的变化为:M()0 M()0+M()1+()2K()0
17、K()0+K()1+()2(13)根据摄动理论,对应的特征值和特征向量的变化为:()0n()0n+()1n+()2v()0n v()0n+v()1n+()2(14)式中为小参数,=0对应于无齿槽结构的标准圆环模型;v(0)n为标准圆环在切向的 n 阶振型函数,参见式(3)和(4)。根据齿槽结构沿圆周的分布特征,定义以下微分算子的摄动量进行表征:图 2 含齿槽结构的定子圆环模型示意图Fig.2 Schematic diagram of ring model of motor stator with cogging structure955振 动 工 程 学 报第 36 卷 M()1=m t/si
18、=1N()-iK()1=kt/si=1N()-i(15)式中m t/s和kt/s为无量纲化的齿质量;为 Dirac 函数;i为第 i 个齿或槽的方位角,i=2()i-1N,N 为总齿数或槽数;齿槽结构可以看作是标准圆环加齿(符号取+)或开槽(符号取-)形成的。将式(13)(15)代入式(12)中,忽略二阶及以上小量,比较的同次幂系数,可得:0:K()0v()0n=()0nM()0v()0n(16)1:K()0v()1n+K()1v()0n=()0nM()0v()1n+()0nM()1v()0n+()1nM()0v()0n(17)将式(17)两边与v()0n进行内积运算(x,y=02xy d),
19、根据 Dirac函数的积分性质可得:|()1n|=|kt/s-()0nm t/s2()n2+1|Ncos()4n()i-1N(18)从式(18)中可以看出,引入齿槽结构后,圆环模态频率的摄动量可能存在两个值,从而可能导致模态频率分裂现象。记此时模态频率分裂的判别式为:=i=1Ncos 4n()i-1N(19)当且仅当=0时,齿槽引起的模态频率摄动量为单值,否则摄动后模态频率存在两个值,发生频率分裂现象。根据三角级数的性质可知:0 频率分裂 2n N 为整数(20)因此,当且仅当阶次 n与总齿槽数 N 满足 2n/N为整数时,含齿槽的圆环发生频率分裂。对于任意一个参数既定的电机而言,总齿槽数N
20、是一个常数,因此总会存在 n 使得 2n/N 为整数(如 n=N),所以齿槽结构必然引起电机定子的模态频率分裂,只不过不同齿槽数下电机定子发生模态分裂的频率阶次有所不同。根据式(20)可知,电机定子的齿槽数目越多,发生模态分裂的频率阶次 n将越大;而对于齿槽数目较少的电机定子而言,发生模态分裂的频率阶次 n 较低,频率分裂对定子振动特性的影响相对较大。2.2轴对称结构引起的定子模态频率分裂现象在电机定子结构上,不仅有类似齿槽的旋转对称结构,还存在诸如底脚、加强肋和拉环等的轴对称结构。如图 3 所示,以对称分布的底脚为例,在 2.1节的理论基础上,研究附加的轴对称结构对电机定子模态频率分裂的影响
21、规律。底脚沿周向的分布方式与齿槽有所不同。设两个底脚与对称轴之间的夹角为((0,)),引入Dirac 函数,可以将底脚引起的摄动量M(1)及K(1)表示为:M()1=m t/s()-+()+K()1=kt/s()-+()+(21)不失一般性,式中采用表示摄动量。类似地,结合式(21),将式(17)两边与v()0n进行内积运算,根据 Dirac函数的积分性质,得到:|()1n|=C sin2(n)或C cos2(n)(22)式中C=|2()kt/s-()0nm t/s()n2+1。从式(22)中可以看出,底脚结构也可能引起定子的频率分裂现象。记此时频率分裂的判别式为:=cos2(n)-sin2(
22、n)=cos(2n)(23)则不发生模态频率分裂的条件为=0,进一步有:=0 频率不分裂 =()2l+1 (4n)(24)式中l为自然数。根据式(24)可知,底脚结构引起的模态频率分裂主要与底脚的位置角相关。对于任意一个既定的,仅当阶次n=(2l+1)(4)时,模态频率才不发生分裂。与含齿槽定子结构的模态分裂条件相比,底脚更容易引起定子结构的模态频率分裂。此外,当=2时,可以发现对于任意 n阶模态频率,均有 0,所以此时定子模态频率一定会发生分裂,而此时相当于旋转对称结构中 N=2 的情况,对应的 2n/N=n 恒为整数,从而根据式(20)也可知此时定子各阶模态频率均发生分裂,两种理论分析的结
23、论是一致的。图 3 含底脚结构的定子圆环模型示意图Fig.3 Schematic diagram of ring model of motor stator with foot structure956第 4 期徐港辉,等:电机定子振动模态频率分裂特性分析3理论模型与分析方法的验证针对本文提出的分析模型及分析方法,本节将通过对比文献和有限元仿真的方式验证其有效性和适用范围。3.1基于圆环模型的定子模态频率对比圆环模型因其简便、高效的优点在求解定子径向模态频率时得到了广泛应用。文献 13 采用数据拟合的方法确定了含齿槽定子铁心的等效圆环内外径,采用机电类比法计算了圆环模型的模态频率,并通过有限元
24、法和模态试验法验证了等效模型的有效性。其所得等效圆环的内直径为 227.5 mm,外直径为 260 mm,密 度 为 7600 kg/m3,弹 性 模 量 为 195 GPa。将以上参数代入本文所建模型中进行计算,并与文献 13 中的结果进行对比,如表 1 和图 4 所示。表 1中的误差是以模态试验结果为基准计算得到的,图 4中的误差归一化处理是以表 1中最大误差(18.03%)为基准进行的。表 1 和图 4 中等效前后分别表示不考虑齿槽和考虑齿槽时的情况,可以看到齿槽等效前的计算误差可达 18.03%,而齿槽等效后的计算误差都在 5%以内,因此齿槽结构对定子模态频率的影响不可忽略。根据等效后
25、的圆环参数,本文圆环模型的计算结果在 2阶频率处取得了极小的误差,在 3阶和 4阶时误差也均不超过 7.1%。进一步分析本文计算误差的形成原因可以发现,本文圆环模型是基于小厚径比的假设建立的,而根据文献 11 可知,机电类比法的计算公式中引入了一个影响系数来表征定子轭厚对定子模态频率的影响。该影响系数的取值通常小于 1,且厚径比越小,该影响系数越接近于 1。而当该影响系数取 1时,可以验证机电类比法所得结果与本文模型结果基本一致,这一方面验证了本文模型的有效性,另一方面也证明了厚度对定子或圆环的模态频率有重要影响。当定子及其等效圆环的厚度相对较大时,本文模型的计算误差也会有所增大。为了验证以上
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电机 定子 振动 频率 分裂 特性 分析