首页/文章/ 详情

基于谐响应的大型异步电机电磁振动分析

13天前浏览572

研究思路:首先根据电机参数建立完整的电机三维模型,包含定子、转子、机壳等部分。通过模态分析求解出该电机的固有频率及振型。建立电机定转子的二维模型,采用有限元法进行电磁计算,通过麦克斯韦张量法求出电机定子受到的径向电磁力。将计算的径向电磁力加载至电机三维模型的定子齿部,通过谐响应分析得到电机振动速率的频域响应。为了验证所采用计算分析方法的有效性,对电机进行振动实验研究。采用振动传感器测量电机不同运行状态时的振动速率,得到电机的固有频率。将实验结果与模态分析结果及谐响应分析结果进行对比,通过模态分析和谐响应分析能够较为准确计算出电机固有频率和振动速率的频谱。

1 异步电机径向电磁力计算

电磁力的计算方法主要有安培力定律、虚位移法、麦克斯韦应力法等。其中,基于有限元的麦克斯韦应力法能够计算电磁力分布,可以得到电机定子所受径向电磁力。(这里忽略,需要的可以自行查询基础知识)

1.2 径向电磁力计算

采用有限元瞬态场对一台3.5MW的三相鼠笼异步电机进行电磁分析,基于电机二维模型进行计算,如图1所示。相较于二维模型,虽然三维模型可以更真实地反映电磁力的轴向分布特性,但本文所分析的电机不存在定子斜槽等设计,且主要关注电机径向电磁力,因此适合使用二维模型计算。

根据实际电机材料属性,对电机各部分材料进行设置,将电机定子外壳设置为零矢量边界。根据模型各部分的尺寸进行网格剖分,为提高气隙磁场计算的精度,对气隙附近网格进行加密。

图1 电机的二维模型

为求解电机空载情况下定子受到的径向电磁力,给电机施加10kV的工频50Hz电压源激励进行瞬态场分析。为保证电机达到稳态,计算总时长为1.5s,计算步长Δt=0.000 2s,对应采样频率为fs=1/Δt=5000Hz。用fF表示电磁力的频率,根据采样定律fs>2fF,可以得到fF最大为2500Hz,能够满足振动分析的需要。计算得到电机空载稳定运行时的磁场分布,如图2所示。

图2 电机稳态磁感线分布

为了计算电机定子所受径向电磁力,需要提取气隙磁感应强度,图3所示为电机稳定后某时刻气隙中磁感应强度沿圆周方向的分布。

根据旋转磁场理论,气隙旋转磁场为行波,即某一位置处磁感应强度在一个电周期内随时间交变。在电机稳态后选取一个周期,分析气隙中任意一点的磁感应强度,将其分解为径向分量和切向分量,如图4所示。可见,径向磁感应强度占主要部分,将导致径向电磁力较大,这与理论相符合。

图3 气隙磁感应强度的空间分布

图4 磁感应强度径向分量和切向分量的时间分布

根据电机磁场分布,采用麦克斯韦应力法计算电磁力,得到单位面积径向电磁力在一个周期内的变化,如图5所示。

对图5中的径向电磁力进行傅里叶变换,得到径向电磁力的主要阶次,如图6所示。各阶径向电磁力中0Hz和100Hz占比最大,其中0Hz的电磁力不随时间变化,作用于电机上表现为恒定的应变,不会引起电机的振动。100Hz的电磁力由基波磁场产生,是引起电机振动的主要原因。在电机稳定运行后,取一个完整电周期的径向电磁力进行谐响应分析,加载于定子齿部的表面,作为振动分析的载荷。

2 异步电机的模态分析

2.1电机模型及模态分析

模态分析是动力学分析的基本方法之一,它可以求解电机的各阶固有频率及对应振型。在模态分析中,采用Hamilton原理建立动力学方程。

基于上述理论,采用有限元法对电机定子铁心进行模态分析,可以得到电机各阶的固有频率和振型,图7所示为轴向0阶的部分振型,其中,径向2阶振型的固有频率为271.61Hz,径向3阶振型的固有频率为685.22Hz,径向4阶振型的固有频率为1138.10Hz

图7 定子铁心各阶模态及对应频率

从图7可以看出,若仅对定子铁心进行模态分析,计算的定子铁心固有频率高于电机正常工作时主要激励源的频率,无法反映出整个电机的低阶模态。因此有必要对电机整体进行模态分析。

完整的电机模型包括定子铁心、端盖、转子以及加强筋等结构。图8所示为完整的电机模型主要组成部分的俯视图。图8直角坐标系中,Z方向为电机的轴向,Y方向电机水平方向,X方向与YZ平面垂直。

在电机底座的螺栓孔处添加Bushing约束及弹簧单元,模拟电机在地面的固定状态。不考虑定子硅钢片、支撑结构间的相对滑动,在阻尼比为零的条件下进行模态分析。图9展示了前四阶固有频率对应的电机振型。

图8 完整的电机模型

图9 完整电机的低阶模态

与单独定子铁心模型相比,完整电机模型的振型更为复杂。在任一固有频率下,电机在xyz三个方向的振动位移和运动质量各不相同,可以据此总结出振型的特点。本文中使用参与系数γ描述电机在每个方向振动的程度,其定义为:

其中:ϕ为式动力学方程在频率下的归一化的特征向量;d为方向矢量。参与系数绝对值越大表示更为主要的振型形式。在xyz各分量下,将所有频率中绝对值最大的值设为1,其余值与之相比,就能清晰判断出振型形式。表1中记录了图9中4个模态对应的参与系数比值,其中第一阶、第三阶Y方向参与系数比值最大,振型主要以Y方向的运动为主;第二阶、第四阶X方向参与系数比值最大,以X方向的运动为主。由上述结果可得,相较于定子铁心的模型,完整电机模型的固有频率更低,更能反映出电机的实际振动特性。

表1 固有频率及其振型参数

固有频率HzX方向参与
系数比值
Y方向参与
系数比值
Z方向参与
系数比值
35.710.270.40
43.850.370.140
66.17
0.04
0.150.03
98.42
.0
0
0.04

2.2 电机模态实验结果与仿真结果对比

对于中小型电机,通过锤击法进行模态测试并取得了很好的效果,但对于大型电机而言很难通过锤击法激发出振动模态,而且复杂的机壳结构会导致激振点的选取很困难。案例根据电机实际应用的工况,通过电机升速和超速实验直接测量电机在多个转速下的振动响应,从而得到固有频率。

实验中,使电机在恒定磁通下从0加速至额定转速3000 r/min, 完成升速实验。然后在恒定电压10kV下电机从3000 r/min加速至3600 r/min, 完成超速实验。通过在电机端盖安装的速度型振动传感器测量X、Y两个方向的振动速率。实验所用电机的转子动平衡度为G0.7,因此可以忽略转子动平衡对实验结果的影响。根据升速、超速实验的测量数据可以得到各转速下电机的振动频谱,从而分析出电机的固有频率。实验得到的固有频率结果如表2所示。

表2 实验测量的固有频率

阶次固有频率Hz
主要振动方向
1
33.33

Y

2
38.91
X
3
60.40
Y
4
100.00
X


从实验结果中可以看出,第一阶、第三阶固有频率在Y方向的振动幅度大于其他方向;第二阶、第四阶固有频率在X方向的振动幅度大于其他方向。

将表1中的仿真结果与表2中的实验结果进行对比,如图10所示,固有频率比较接近,对应的振型方向也一致,仿真结果与实验结果具有较高的一致性。结果表明,对完整电机的模态分析能够较为准确地计算出电机的固有频率,电机每个结构都会对各阶模态的频率及振型产生一定的影响,因此建立完整的电机仿真模型可以得到更接近实际情况的分析结果。

图10 固有频率的实验与仿真结果对比

3 异步电机的谐响应分析

3.1谐响应分析

当电机稳定工作时,电机定子受到的径向电磁力可分解为一系列不同频率的正弦分量,是周期性简谐载荷,可以作为谐响应分析的激励。同时,正常工作中电机的形变很小,可忽略非线性特征。因此可以通过谐响应分析研究径向电磁力作用下电机的振动响应。

将磁场分析中计算得到的电磁力加载于电机三维模型的定子齿端,计算0至500 Hz范围内电机端盖处的频率响应,结果如图11所示。图11中纵轴为振动速率,采用对数刻度,单位为m/s;横轴为振动频率,单位为Hz。在X、Y方向的振动速率的频谱中,100Hz对应的分量占比最大,其他分量可以忽略,这与径向电磁力的分布相对应。

图11 端盖处X、Y方向振动速率仿真计算结果

3.2谐响应实验结果与仿真结果对比

在10kV额定电压激励下,测量电机空载状态时的振动响应。在电机端盖安装振动传感器进行测量,电机到达稳态后,连续10分钟记录X、Y方向的振动速率并计算其频谱。图12为10kV电压下X、Y方向振动速率的频谱。

从图12中可以看出,测量结果中50Hz和100Hz对应的振动分量占比较大,50Hz的振动主要由转子旋转的机械振动引起的,100Hz的振动主要由径向电磁力引起的。对比100Hz下实验和仿真的振动速率,如表3所示,可以看出X、Y方向上振动速率的仿真结果与实验结果都比较一致,证明了所采用的振动分析方法的有效性。

表3 100 Hz振动速率的实验与仿真结果对比

振动方向X(10^-4m/s)
Y(10^-4m/s)
实验3.8
6.0
仿真
3.75
6.5

图12 电机振动速率频谱


4 结 论

1)通过对比大型异步电机的定子三维模型和整机模型的模态分析结果,发现定子模型的固有频率较高;整机模型的固有频率更丰富,包含许多低阶固有频率,且电机的振型主要表现为径向的平移。因此在对大型异步电机进行模态分析时需要考虑支撑结构及转子带来的影响。

2)根据谐响应分析能够较为准确地计算出电机的电磁振动响应。从电机振动的频谱可以看出,大型异步电机的电磁振动主要表现为2倍电频率的振动,这与电磁力频谱中幅值最大的2倍频分量相对应;高阶电磁力引起的振动占比很小。

3)电磁振动的仿真分析结果与实际测量值存在误差。大型异步电机的振动实验测试结果受到测量位置、电机固定强度、部件连接刚度等因素影响,而有限元仿真中,虽然可以通过设置电机材料、固定方式、阻尼等参数对实际情况进行模拟,但无法完全与实际工况保持一致。

4)仿真和实验结果证明了所采用振动分析方法的有效性。振动分析能够为电机的振动抑制提供依据,评估电机的振动指标是否合格,从而为电机设计或控制算法的优化提供参考。


网络整理,禁止商用,仅供内部分享

公 众号:机电君



来源:机电君
振动非线性电磁力理论电机材料控制螺栓
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-09-01
最近编辑:13天前
ErNan.Chen🍃
硕士 | CAE工程师 即物而穷其理
获赞 28粉丝 35文章 49课程 0
点赞
收藏
作者推荐

摩擦型高强螺栓节点微动疲劳性能试验及有限元分析

引言:在现实工程结构中,摩擦型高强螺栓搭接板连接面不同的表面处理方式导致了表面摩擦系数的不同,在施工过程中,由于工人操作以及现场不确定因素较多,摩擦型高强螺栓的预紧力与设计值易出现偏差,这些情况都会影响构件的微动疲劳寿命。总结:基于Q235钢材、10. 9级摩擦型高强螺栓连接件微动疲劳寿命试验和断口扫描电镜分析结果,通过有限元方法建立Q235高强螺栓单面连接微动疲劳寿命数值预测模型。在此基础上开展螺栓预紧力和接触面摩擦系数对高强螺栓单面连接微动疲劳寿命的影响分析,考察螺栓预紧力的不足和连接板接触面摩擦系数的改变对接触面正应力、切应力和相对位移的影响以及数值的变化。分析结果表明: 螺栓预紧力和接触面摩擦系数对微动疲劳寿命都有一定的影响,其中螺栓预紧力较JGJ 82—2001《钢结构高强度螺栓连接技术规程》要求低10%时,微动疲劳寿命降低约12% ; 连接板接触面摩擦系数每减小0. 5,微动疲劳寿命降低约10% 。1 微动疲劳试验1 试件设计试验设计制作了9个试件,采用两块大小相同( 厚度为5mm,宽度为80 mm) 的Q235钢板组成连接节点,Q235钢材力学性能及化学成分见表1和表2。板材通过两个M20的10. 9级摩擦型高强螺栓搭接连接,试件详细尺寸见图1a。钢板表面经喷丸处理,并按照 GB 50017—2003《钢结构设计规范》及GB50205—2001《钢结构工程施工质量验收规范》要求,拧紧过程分为初拧和终拧,螺栓终拧预紧力为155kN。表 1 实测 Q235 钢材力学性能表 2 实测 Q235 钢材化学成分图 1 试件尺寸及试验装置2 试验装置及加载制度疲劳试验机采用MTS-50电液伺服万能试验机,试验采用的应力比为0. 1,频率为30Hz,试验环境为常温实验室内大气环境。不同工况对应的最大应力值及对应的微动疲劳寿命见表3。当连接板最终断裂时视为失效,当寿命达到10的7次即认为材料疲劳寿命无限大。表 3 试验荷载基本情况3 试验主要现象及分析试验现象为在疲劳荷载达到一定循环次数之后,连接件作动器一侧螺栓孔边缘处在先前没有预兆的情况下突然出现断口( 图 2a) ,疲劳试验机由于自身的位移保护机制而出现停机,摩擦型高强螺栓并未出现破坏,说明了摩擦型高强螺栓连接件的疲劳破坏属于脆性断裂。断口位置为靠近作动器一侧钢板在螺栓孔外5mm左右处( 图 2a) ,而不是发生在螺栓孔直径横截面积最小处。断口周围有暗红色粉末,可以断定此处两接触面之间发生了微动并产生了微动磨损,并对裂纹的产生和扩展有一定的影响。由扫描电子显微镜( SEM) 观察断裂表面发现,试件破坏属于微动疲劳断裂,有明显的微动疲劳源和裂纹扩展区。裂纹从接触面开始形成疲劳源,几个疲劳源相互连接并向外部扩展( 图 2b) ,最终导致试件断裂失效。图 2 破坏现象及断口扫描电镜试验结果2 微动疲劳寿命有限元分析2. 1 有限元模型采用ABAQUS有限元软件对试验进行模拟。由于试件在长度方向是轴对称试件,故取一半模型进行分析。模型中,垫片和螺栓的材料相同,可将其视为一体进行简化。模型中Q235连接板的材料本构通过试验获得,屈服强度为321. 5MPa,抗拉强度为482. 2MPa,弹性模量为1. 87 × 10^5 MPa,泊松比为0. 3。对于10. 9级摩擦型高强螺栓,根据JGJ 82—2011《钢结构高强度螺栓连接技术规程》可知: 螺栓预紧力为155kN,螺栓杆有效截面积为245mm^2,截面应力为632. 6 MPa;10. 9级摩擦型高强螺栓屈服强度根据GB /T1231—2006《钢结构用大六角螺栓》可知: 屈服强度为940MPa,泊松比为0. 3,故螺栓使用阶段处于弹性段。网格划分时,由于此模型中要考虑接触面之间的微动效应以及摩擦力,线 性减缩积分单元( C3D8R) 和非协调单元( C3D8I) 都可以使用。同样情况下,相比非协调单元,线性减缩单元可以大大缩短计算时间,所以在模型中全部使用线性减缩积分单元( C3D8R) 。考虑应力梯度的变化规律,应力集中较为明显的螺栓孔周围网格划分较密,而边缘处网格相对稀疏,这样既能保证运算精度同时又可提高计算效率。为了方便计算特征点之间的相对位移,网格划分时,接触面上两块板的结点是一一对应的,并对图3 中靠近模型加载端的接触面上0°,90°和180°方向的结点进行了标号( 图 4) 。模型建成后,需要施加相应的螺栓预紧力,设置不同的接触面摩擦系数,之后对试件施加与试验相等的拉应力。图 3 试件有限元模型应力云图图 4 有限元模型结点编号2. 2 试验数据与有限元分析结果对比通过式( 1a) 计算tFFD,并通过用文献中提出的预测微动疲劳寿命曲线 tFFD-N( 即式( 1b) ) 计算微动疲劳寿命。表 4 给出了试验寿命与微动疲劳的预测寿命,寿命的预测模型相对比较准确。图 5 给出了传统S-N与 tFFD-N 曲线与试验疲劳寿命之间的关系。文献:孙先磊. 摩擦型高强螺栓搭接连接微动疲劳试验研究[D]. 西安: 西安建筑科技大学,2015通过图5可以看出,tFFD-N曲线与试验疲劳寿命在总体趋势上表现大体一致。从图6可以看出:tFFD-N 寿命预测模型预测的疲劳寿命与试验寿命误差较小,基本呈现出线性关系,可以说明此模型对摩擦型高强螺栓微动疲劳寿命有一定的指导意义。表 4 试验寿命与模型预测寿命对比图 5 tFFD -N,S-N 关系试验结果图 6 试验寿命与拟合寿命2. 3 参数分析2. 3. 1 接触面摩擦系数材料表面处理方法有很多种,对应的表面摩擦系数亦不同。表面处理方式与摩擦系数的关系见表5。在施工过程中,也易出现高强螺栓扭矩不足的情况。在ABAQUS中,设置不同的接触面摩擦系数代表相应的表面处理方式; 不同的螺栓预紧力代表螺栓欠拧的程度。取试验时2种最大拉应力幅值( 225,237. 5MPa) ,通过改变螺栓预紧力和接触面的摩擦系数,分析每种情况下接触面特征点的接触正应力、接触切应力和相对位移改变量。表 5 Q235 钢材连接件不同摩擦面处理方法及摩擦系数在ABAQUS模型中,接触面摩擦系数取0. 45,0. 40,0. 35,0. 30,0. 25 共5种,以此代表构件表面的不同处理方法。此时螺栓预紧力保持不变,取规范值155kN。经计算,在螺栓孔的3个方向中0°方向是tFFD最大值出现的方向,也是最危险方向。因此重点分析模型 0°方向特征点的接触正应力( 图7) ,接触切应力( 图 8) 和相对位移( 图 9) 。由图7可知: 距螺栓孔越远,接触正应力越小,曲线变化规律相似。接触正应力来源于螺栓的预紧力,当预紧力不变,只改变接触面的摩擦系数时,对接触正应力影响不大。a—S = 237. 5 MPa; b—S = 225 MPa。—■—μ = 0. 45; —●—μ = 0. 40; —▲—μ = 0. 35; ——μ = 0. 30。图 7 不同摩擦系数下的接触正应力由图8可知: 接触切应力是由1号结点逐渐增大到10号结点处(到达最大值),之后逐渐减小。不同摩擦系数的切应力曲线在2 ~ 8号结点之间,接触切应力随着摩擦系数的增大而增大,在10号结点处,切应力虽然也是随着摩擦系数增大而增大,但最大切应力变化幅度很小。因此当接触面的摩擦系数改变时,对接触切应力最大值影响不大。a—S = 237. 5 MPa; b—S = 225 MPa。—■—μ = 0. 45; —●—μ = 0. 40; —▲—μ = 0. 35; ——μ = 0. 30。图 8 不同摩擦系数下的接触正切应力由图9可知: 同一疲劳荷载下,在1号结点处,摩擦系数的改变不影响接触面的相对位移,但在3号结点之后,接触面相对位移发生有规律的变化,即随摩擦系数的降低,接触面间的相对位移逐渐增加。a—S = 237. 5 MPa; b—S = 225 MPa。—■—μ = 0. 45; —●—μ = 0. 40; —▲—μ = 0. 35; ——μ = 0. 30。图 9 不同摩擦系数下的结点相对位移由图10可知: tFFD曲线变化明显,尤其是在tFFD值最大的10号结点,随着摩擦系数的减小,tFFD值增大。a—S = 237. 5 MPa; b—S = 225 MPa。—■—μ = 0. 45; —●—μ = 0. 40; —▲—μ = 0. 35; ——μ = 0. 30。图 10 不同摩擦系数下的 tFFD综上,其他外界条件不变的情况下,随着接触面的摩擦系数减小,导致接触面之间相对滑动量增加,进而导致了tFFD值的增大。根据tFFD-N拟合的试验寿命预测曲线,在螺栓预紧力和拉应力不变的情况下,对接触面不同摩擦系数的寿命进行预测。从图11可以看出: 当应力幅比较大,且低周疲劳时,摩擦系数与预测寿命呈线性关系,当表面摩擦系数μ 减小0. 05时,疲劳寿命减小约10% 。—■—S = 237. 5 MPa; —●—S = 225 MPa图11 摩擦系数对微动疲劳寿命的影响2. 3. 2 螺栓预紧力在ABAQUS有限元分析中,螺栓预紧力取JGJ82—2011要求预紧力的 100% ,90% ,85% ,80% ,70% 共5种情况,来代表不同欠拧情况下的螺栓预紧力。最大拉应力取237. 5,225MPa。此时表面摩擦系数取0.45不变。不同螺栓预紧力下0°方向特征点的接触正应力、接触切应力和相对位移由ABAQU 得出。案例只给出S=237. 5MPa时各应力与相对位移变化见图12、图13。由图12a可以看出: 接触正应力由螺栓孔边缘向外逐渐减小,在1~8号结点,预紧力对正应力影响较为明显,正应力随着螺栓预紧力减小而减小。从9号结点之后,出现了螺栓预紧力越大,接触面的正应力反而随着预紧力的增大而减小的现象,不同螺栓预紧力下各结点的接触正应力变化较小。由图12b可以明显看出: 当摩擦系数不变,在不同螺栓预紧力下,接触切应力的变化趋势是相同的。对于不同的螺栓预紧力,在9号结点之前,预紧力较大的接触面接触切应力较大,9 号结点之后,预紧力越小,接触切应力越大。螺栓预紧力改变对接触切应力的影响很小,不同预紧力的切应力曲线几乎重合。a—正应力; b—切应力。—■—100% ; —●—90% ; —▲—85% ; ——80% ; —◆—70% 。图 12 螺栓预紧力对接触应力的影响( 237. 5 MPa)由图13可以看出: 1~13号结点的相对位移逐渐变大。在1号结点,出现了负向的相对位移。螺栓预紧力对节点相对位移的影响在4~9号结点之间较为明显。螺栓预紧力为108. 5 kN,最大拉应力237. 5MPa工况下,拉应力已经接近了接触表面滑动摩擦力的值,在此情况下出现较大的相对位移值,而当其他情况下,接触表面滑动摩擦力极限值远大于拉应力,因此螺栓预紧力减小,接触面上相对滑移量有少量增大。—■—100% ; —●—90% ; —▲—85% ; ——80% ; —◆—70%图13 螺栓预紧力对接触面相对位移的影响由图14可以看出: 当拉应力值较大时,螺栓预紧力减小,0°方向的最大 tFFD值有所增大,tFFD值在一定程度上反映了磨损的程度。a—S = 237. 5 MPa; b—S = 225 MPa。—■—100% ; —●—90% ; —▲—85% ; ——80% ; —◆—70% 。图 14 螺栓预紧力对 tFFD的影响对不同螺栓预紧力进行寿命预测,由图15可以看出,随着螺栓预紧力的降低,微动疲劳寿命降低。当拉应力比较大,螺栓预紧力与寿命呈近似线性关系,且螺栓预紧力每降低10% ,寿命约降低12% 。—■—S = 237. 5 MPa; —●—S = 225 MPa。 图 15 螺栓预紧力对疲劳寿命的影响结束1) 试验表明: Q235钢材与10. 9级摩擦型高强螺栓构成的单面连接件,当疲劳荷载的应力比一定时,随着最大拉应力的增大,钢材的微动疲劳寿命降低。2) ABAQUS有限元分析表明: 随着螺栓预紧力减小,微动疲劳寿命降低; 接触面摩擦系数越小,微动疲劳寿命越低。3) 构件表面摩擦系数主要影响接触点之间的相对滑移量,从而影响tFFD的值,即影响微动疲劳寿命,当摩擦系数减小时,接触面上相对滑移量增加,微动疲劳寿命减小。当表面摩擦系数μ从0. 45每减小0. 05时,微动疲劳寿命约减小10% 。4) 当螺栓欠拧时,对接触面的接触正应力有一定影响,对接触切应力影响不明显,对相对位移略有影响,但几种影响相互累加,导致螺栓预紧力相对规范规定每降低10% ,寿命降低约12% 。网络整理,仅限内部分享,禁止商用来源:机电君

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈