绝缘材料

参数	数值
额定功率/kW	2 800
额定电压/kV	10
额定电流/A	211
额定频率/Hz	50
定子铁心外径/mm	850
定子铁心内径/mm	445
导体泊松比	0.326
导体密度/(kg/m³)	8 933
导体剪切模量/GPa	44.87
绝缘层厚度/mm	2
绝缘层弹性模量/GPa	1.1
绝缘层泊松比	0.41
绝缘层密度/(kg/m³)	2 500
绝缘层剪切模量/GPa	0.397

参数	数值
额定功率/kW	2 800
额定电压/kV	10
额定电流/A	211
额定频率/Hz	50
定子铁心外径/mm	850
定子铁心内径/mm	445
导体泊松比	0.326
导体密度/(kg/m³)	8 933
导体剪切模量/GPa	44.87
绝缘层厚度/mm	2
绝缘层弹性模量/GPa	1.1
绝缘层泊松比	0.41
绝缘层密度/(kg/m³)	2 500
绝缘层剪切模量/GPa	0.397

复杂电磁力作用下10 kV交流电机端部绕组绝缘损伤应力强度分析

PDF下载

张海军 ¹^,² , 张明杰 ¹ , 陈小勇 ² , 万少华 ¹ , 郭季 ¹

绝缘材料 | 绝缘技术 2023,56(1): 45-50

收起

绝缘材料 | 绝缘技术 2023, 56(1): 45-50

复杂电磁力作用下10 kV交流电机端部绕组绝缘损伤应力强度分析

全屏

张海军¹^,², 张明杰¹, 陈小勇², 万少华¹, 郭季¹

作者信息

¹湖北文理学院纯电动汽车动力系统设计与测试湖北省重点实验室,湖北襄阳 441053

²新昇电气有限公司,湖北襄阳 441899

张海军（1981-），男（汉族），河北保定人，副教授，博士，研究方向为电机多物理场数值计算与故障诊断。

Stress intensity analysis on insulation damage of 10 kV AC motor end winding under complex electromagnetic force

Haijun ZHANG¹^,², Mingjie ZHANG¹, Xiaoyong CHEN², Shaohua WAN¹, Ji GUO¹

Affiliations

¹Hubei Key Laboratory of Power System Design and Test for Electrical Vehicle, Hubei University of Arts and Science, Xiangyang 441053, China

²NRE Electric Co., Ltd., Xiangyang 441899, China

出版时间: 2023-01-20 doi: 10.16790/j.cnki.1009-9239.im.2023.01.008

文章导航

摘要

收起

针对10 kV大功率牵引电机端部绕组绝缘层的早期故障问题，基于断裂力学理论分析了绝缘层裂纹损伤在复杂电磁力、机械力作用下的破坏机理。利用有限元方法建立三维分析模型，并通过分段加载电磁力对不同位置的绝缘层裂纹应力应变状态进行数值计算，同时引入应力强度因子对裂纹损伤扩展的程度进行评估。结果表明：与鼻端位置绝缘层裂纹损伤相比，线棒关节处裂纹扩展的速度更快、断裂程度更大；裂纹的初始损伤状态对断裂及扩展程度也有重要影响，横向裂纹更容易扩展。

关键词

电机端部绕组 / 断裂损伤 / 有限元分析 / 裂纹扩张 / 故障机理

Abstract

收起

Aimed at the early fault of 10 kV high power traction motor end winding insulation layer, the failure mechanism of insulation layer crack damage under complex electromagnetic and mechanical forces was analyzed on the basis of fracture mechanics theory. A three-dimensional analysis model was established by using the finite element method, and the stress-strain states of the insulating layer crack at different positions were calculated numerically by piecewise loading electromagnetic force. At the same time, the stress intensity factor was introduced to evaluate the degree of crack damage expansion. The results show that compared with the insulating layer crack damage at the nose end, the crack expansion speed at the rod joint is faster and the fracture degree is greater. The initial damage state of the crack also has an important influence on the fracture and expansion degree, and the transverse crack is easier to expand.

Key words

motor end winding / fracture damage / finite element analysis / crack expansion / fault mechanism

引用本文

张海军, 张明杰, 陈小勇, 万少华, 郭季. 复杂电磁力作用下10 kV交流电机端部绕组绝缘损伤应力强度分析. 绝缘材料, 2023 , 56 (1) : 45 -50 . DOI: 10.16790/j.cnki.1009-9239.im.2023.01.008

Haijun ZHANG, Mingjie ZHANG, Xiaoyong CHEN, Shaohua WAN, Ji GUO. Stress intensity analysis on insulation damage of 10 kV AC motor end winding under complex electromagnetic force[J]. Insulating Materials, 2023 , 56 (1) : 45 -50 . DOI: 10.16790/j.cnki.1009-9239.im.2023.01.008

正文

收起

0 引言

收起

牵引电机作为交流电力机车的核心设备之一，其健康状态关系到整个列车的安全与稳定运行。一旦发生事故，除了电机本体检修周期比较长，还可能造成巨大经济损失，甚至带来严重的社会影响。统计表明，端部绕组故障是直接或间接导致牵引电机事故的主要原因之一^[1-2]。随着机车牵引速度与功率的增大，受复杂、恶劣工作环境因素（高温、负载波动大、高频振动）影响，牵引电机绕组将长期工作在以强电磁力为主要力，辅以机械应力和热应力等综合力的作用下，导致电机端部绕组故障概率显著提高^[3-4]。端部绕组绝缘层破坏是典型的电机故障之一，如图1所示。

目前关于端部绕组的研究主要集中在绕组电磁力计算、电磁力引起的动态变形以及振动分析方面，文献[5]基于有限元法分析了发电机定子端部绕组电磁力及端部绕组单匝的应力和应变分布。文献[6]采用有限元体积法研究了在额定载荷条件下定子端部绕组表面的传热系数分布、流体流动分布和端部绕组的温度分布。文献[7]提出了一种同时考虑端部绕组边缘效应和漏磁通的新方法，可用于计算端部绕组的漏电感。文献[8]采用有限元法研究了过渡状态和稳态的电磁场，用数值方法研究了电机端部绕组电磁力引起的振动。文献[9]通过三维有限元法得出了电机稳态运行时整个端部绕组的电磁力和应力分布，并和铜的屈服强度进行了比较。

近年来，国内外学者通过在电机端部绕组电磁力计算宏观受力及其振动规律方面的深入研究，获得了足够精度的端部绕组电磁力分布，然而绕组受电磁力作用引起导体变形以及绝缘层破坏的机理还有待进一步研究。统计表明^[10-11]，端部绕组早期故障形式主要表现在绝缘层的破坏，例如出现不同程度的裂纹损伤。随着牵引电机频繁启动和加减速的工况变化，端部绕组绝缘层裂纹在综合应力作用下将出现不同程度的扩展，导致绝缘加剧破坏，严重时甚至会引起绝缘击穿和短路故障。

本研究以某型号大功率牵引电机为例，建立其端部绕组的三维有限元模型，并利用有限元方法计算裂纹在电磁力作用下的应力应变状态。基于断裂力学理论，探究端部绕组绝缘层裂纹扩展和破坏机理，并通过应力强度因子对裂纹扩展程度进行评估。

1 绕组绝缘损伤断裂理论

收起

1.1 电磁与机械应力应变分析

电机端部绕组的电磁力分布复杂且具有交变特性，会引起绝缘层内部不规则的应力应变状态变化。因此，确定绕组导体和外层绝缘的内应力应变状态是进行端部绕组绝缘损伤断裂机理研究的前提^[12]。

根据洛伦兹力定律，电磁力可表示为式（1）。

F = J × B

（1）

式（1）中：F为电磁力；J为电流密度；B为磁通密度。

根据弹性形变分析理论，应力应变控制方程如式（2）所示^[13]。

∂ σ r ∂ r + ∂ τ z r ∂ z + σ r - σ θ r + f r = 0 ∂ σ z ∂ z + ∂ τ r z ∂ r + τ r z r + f z = 0

（2）

式（2）中：σ_z为z方向应力分量；σ_r为r方向应力分量；σ_θ为θ方向应力分量；τ_rz为剪切应力分量。

表征绕组铜导体材料内部应力应变状态关系的几何方程组如式（3）所示。

ε r = ∂ u r ∂ r ε θ = u r r ε z = ∂ u z ∂ z γ r z = ∂ u r ∂ z + ∂ u z ∂ r

（3）

式（3）中：ε_z为z方向应变分量；ε_r为r方向应变分量；ε_θ为θ方向应变分量；γ_rz为剪切应变分量；u_z为z方向位移；u_r为r方向位移。

一般应力状态下的各向同性材料本构方程组如式（4）所示。

ε r = 1 E [σ r - υ (σ θ + σ z)] ε z = 1 E [σ z - υ (σ r + σ θ)] ε θ = 1 E [σ θ - υ (σ z + σ r)] γ z r = 2 (1 + υ) E τ z r

（4）

式（4）中：E为弹性模量；υ为泊松比。

1.2 绝缘层损伤断裂的应力强度因子

电机端部绕组不可避免地会受到加工、装配过程中的机械碰撞以及质量缺陷的影响而产生初始损伤，主要表现为埋藏裂纹和表面裂纹。此外，受材料老化因素影响，绕组长期工作在复杂电磁-机械应力综合作用下，绝缘层裂纹将进一步扩展演化^[14]。在电机绝缘中天然材料云母占主要地位，为了量化评估和表征绕组绝缘层受应力作用后引起的裂纹断裂和扩展程度，本研究以云母材料的表面裂纹为例^[15]，引入应力强度因子对裂纹断裂程度进行量化描述，如式（5）所示。

K Ι = 2 π 2 G 1 + k V r K Ι Ι = 2 π 2 G 1 + k U r K Ι Ι Ι = 2 π 2 G W r

（5）

式（5）中：K_I、K_II、K_III分别对应张开型、滑开型、撕开型裂纹尖端的应力强度因子；U、V、W分别为局部笛卡尔坐标系下的位移分量；r为局部圆柱坐标，如图2所示；G为剪切量；k如式（6）所示。

k = 3 - 4 υ 平面 应变 3 - υ 1 + υ 平面 应力

（6）

在断裂理论力学中，应力强度因子是表征裂纹尖端应力应变场的重要参数，其值越大，裂纹越容易扩展。用应力强度因子表征裂纹尖端节点位移，如式（7）所示。

U = K Ι 4 G r 2 π 2 k - 1 c o s θ 2 - c o s θ 3 - K Ι Ι 4 G r 2 π 2 k + 3 s i n θ 2 + s i n 3 θ 2 V = K Ι 4 G r 2 π 2 k - 1 s i n θ 2 - s i n 3 θ 2 - K Ι Ι 4 G r 2 π 2 k + 3 c o s θ 2 + c o s 3 θ 2 W = 2 K Ι Ι Ι G r 2 π s i n θ 2

（7）

2 电机端部绕组三维有限元建模

收起

2.1 电机及端部绕组参数

本研究以一台三相交流牵引电机为例进行端部绕组绝缘损伤断裂及破坏强度分析，三维模型如图3所示，模型主要参数如表1所示。

2.2 端部线圈建模与裂纹初始状态模拟

在牵引电机实际运行中，其端部绕组长期处于不均衡、交变的电磁力作用下，端部绕组受力非常复杂且呈现不对称性，特别是绑扎的端部绕组线圈容易受振动影响出现松动，应力应变状态较为复杂^[16]。本研究重点针对端部绕组松动后，线圈绝缘层损伤断裂的规律进行分析，忽略线圈之间的接触约束影响，同时为了减少有限单元数和计算量，使用端部绕组单根绕组线圈的简化模型，如图4所示。为了提高计算精度，并考虑端部线圈结构和受力的不对称性，采用三维有限元法对其进行数值分析。

为了更准确地分析端部绕组在电磁力作用下的绝缘层损伤演化，本研究中线圈模型考虑了实际绕组的构成。图4中线圈模型包括铜导体和外部绝缘层，黄色部分为端部绕组内部的线圈导体，红色部分为线圈绝缘。建模中，将端部线圈外层绝缘设置为一个整体，而将导体分为16段，其中线棒部分各分为7段，线圈鼻端分为2段，以便更真实、精确地对导体加载电磁力。加载后的电磁力如图5所示。

2.3 绝缘层初始裂纹状态模拟

电机端部绕组短路早期故障主要表现为绝缘层裂纹损伤，本研究对绝缘层进行初始裂纹建模。通过在绝缘层表面设置半椭圆裂纹来描述绕组绝缘的初始损伤状态，如图6所示。图6中l表示初始裂纹的主半径，b表示初始裂纹的最大轮廓半径，h表示初始裂纹的次半径。考虑到裂纹与绝缘层有限单元结构的影响，绕组导体、绝缘层分别采用六面体和四面体单元进行剖分，网格剖分单元尺寸为10 mm，总单元数为34 853，总节点数为105 125。图7为初始裂纹尖端网格划分放大情况。

3 算例及结果分析

收起

利用有限元软件对电机端部绕组进行瞬态磁场仿真计算，求解得出绕组云图，其中整体定子端部绕组电磁力密度如图8所示，单根定子端部绕组电磁力密度如图9所示。由图8～9可知，整体电机端部绕组电磁力密度呈现对称且规律性变化，且关节部位电磁力密度大于鼻端部位。

3.1 端部线圈集中应力分析

在电机运行中端部绕组会受到较为复杂的电磁力作用，通常裂纹在应力集中位置更容易扩展。分析有限元数值计算结果，发现端部绕组受电磁力作用后应力集中，应变峰值主要出现在铁心槽口附近，如图10所示。图10中，在铁心槽口处的线圈绝缘层上设置不同的采样点进行应力应变状态比较和定量分析。

图11和图12分别为不同采样点对应的应力应变值，其应力应变值的结果变化与理论上公式（2）、（3）中应力应变的关系相同。从图11～12可以发现，各采样点应力应变状态变化规律基本一致，而采样点3、4、5的应力明显高于其他采样点的应力，且最大应力应变值约为最小应力应变值的2.3倍。

图13为绝缘层应力集中位置存在裂纹初始损伤下的应变状态分布。与图10比较可以发现，应变峰值大小和位置都发生变化：应变最大位置转移到绝缘层裂纹的尖端附近，应变最大值显著增大。

3.2 绝缘层裂纹强度分析

为了研究和定量表征裂纹损伤扩展的影响因素，并揭示裂纹断裂损伤扩展规律，分别对不同位置、不同方向以及不同h值的张开型裂纹进行数值计算和结果分析。

3.2.1 不同位置裂纹的强度分析

电机端部绕组线圈的振动变形与悬臂梁结构力学问题类似，一方面铁心槽口位置的应力应变集中比较明显，另一方面绕组鼻端位置的应力应变状态也值得关注^[17]。针对上述应力应变比较集中的位置，分别在铁心槽口（端部关节）和鼻端的绕组绝缘层设置相同方向、相同尺寸的裂纹，以分析裂纹位置对裂纹扩展程度的影响，线棒关节处裂纹（裂纹1）和鼻端裂纹（裂纹2）的位移分布计算结果如图14所示。

基于图14中不同位置裂纹尖端（裂纹1和裂纹2）的应变状态值分别计算得到其应力强度因子，如图15所示。由15中曲线数据可知，线棒关节处裂纹的应力强度因子最大值为4.0 MPa·mm²，鼻端处应力强度因子最大值为0.5 MPa·mm²。可以看出，端部关节处应力强度因子是鼻端裂纹应力强度因子的8倍，这表明端部绕组关节处的裂纹损伤扩展程度更大，随着电机启动产生的电磁力频繁冲击，端部绕组关节处裂纹损伤也更容易出现进一步断裂。

3.2.2 不同方向裂纹强度分析

实际牵引电机的端部绕组绝缘层初始裂纹裂开方向也具有一定的随机性，图16为不同方向的裂纹分布示意图。为了研究裂纹裂开方向对应力应变状态变化和裂纹扩展程度的影响，本研究针对应力最集中的铁心槽口位置、相同尺寸但不同方向（分别为0°、30°、60°、90°，其中定义0°裂纹为纵向裂纹，90°裂纹为横向裂纹）的绝缘层裂纹分别进行了有限元数值计算。

图17为不同方向的裂纹尖端应力强度因子。从图17可以发现，随着角度从0°逐渐变化到90°，应力强度因子逐渐变大，由60°到90°这一过程中变化幅度最大，且90°对应的横向裂纹应力强度因子最大。这表明90°方向的裂纹更易扩展，使绝缘层损伤程度更大。

3.2.3 不同h值时的裂纹强度分析

为了进一步分析初始裂纹h值对断裂程度及扩展趋势的影响，设定了不同h值的初始裂纹，并设定绕组绝缘总厚度为2 mm，在应力应变集中的同一位置设定大小相同，h值分别为0.4、0.7、1.0、1.3、1.6 mm的裂纹且未设定贯穿裂纹进行计算，分析得出应力强度因子随裂纹h值的变化曲线如图18所示。

由图18可以得出，在应力应变最集中的位置处设定大小相同h值不同的裂纹，当裂纹的h值越大时，其应力强度因子越大，在裂纹h值达到一定程度后其应力强度因子可能达到相对最大值（2.0 MPa·mm²），并在一定区域内上下浮动，与公式（5）、（7）中理论结果相对应。表明当初始裂纹的大小相同时，h值越大，绝缘层的损伤程度越大。

4 结论

收起

本研究分析了电机端部绕组绝缘层的裂纹损伤规律，主要得到以下结论：

（1）端部绕组线棒关节处的裂纹与鼻端位置裂纹的应力强度因子计算结果表明，端部绕组线棒关节处存在的初始裂纹比鼻端位置更容易扩展。

（2）不同方向的裂纹（横向裂纹或纵向裂纹）应力强度因子变化明显，当裂纹呈现为横向裂纹时，裂纹的应力强度因子和裂纹扩展程度更大。

（3）不同h值的初始裂纹应力强度因子也不相同，整体表现为初始裂纹h值越大，应力强度因子也越大，当h值达到一定程度后，其应力强度因子值可能达到最大值。

基金

收起

湖北省自然科学基金资助项目(2019CFC893)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

李仕云,邱忠才,叶洋,等.牵引电机定子绕组匝间故障特征仿真分析[J].电机与控制应用,2014,41(9):69-73.

[2]

郝亮亮,吴俊勇,陈占锋,等.转子匝间短路故障对大型汽轮发电机振动的影响机理[J].电力系统自动化,2014,38(4):25-31,50.

[3]

乔长帅,唐赢武,钟博.牵引电机定子端部绕组电磁力计算分析[J].电机技术,2020(5):15-18,22.

[4]

LIN

, LAIHO

A N

, HAAVISTO

, et al. End-winding vibrations caused by steady-state magnetic forces in an induction machine[J]. IEEE Transactions on Magnetics,2010,46(7):2665-2674.

[5]

张宇娇,吴刚梁,翁杰,等.发电电动机定子端部绕组电磁-结构耦合场计算[J].水电能源科学,2015,33(10):141-144.

[6]

HAN

, LI

, WANG

, et al. Calculation and analysis of the surface heat-transfer coefficient and temperature fields on the three-dimensional complex end windings of a large turbogenerator[J]. IEEE Transactions on Industrial Electronics,2014,61(10):5222-5231.

[7]

LIN

, ZHOU

, HU

, et al. Analytical computation of end-winding leakage inductance for multi-phase AC machines[C]//2017 IEEE International Electric Machines and Drives Conference. Maimi, USA:IEEE,2017:1-6.

[8]

SHAN

, CHENG

X W

, FANG

, et al. Analysis of vibration caused by electromagnetic force on stator end-winding of large dry submersible motor[J]. Applied Mechanics & Materials,2013,416-417:428-432.

[9]

FANG

, LÜ

, CHENG

, et al. Analysis of stress distribution on end winding of large water filling submersible motor during steady state operation[C]//2013 5th International Conference on Power Electronics Systems and Applications. Hong Kong, China:IEEE,2013,1-6.

[10]

杜协和,陈求索,王健,等.10kV防爆电机定子绕组匝间绝缘结构研究[J].绝缘材料,2020,53(10):51-54.

[11]

陈阳,孙玉田,孙永鑫.大型电机定子线圈及绕组绝缘介质损耗因数试验标准分析综述[J].绝缘材料,2020,53(1):9-15.

[12]

陆宝琦.交流变频电机的绝缘[J].绝缘材料,2001,34(3):29-34.

[13]

张海军,张华,马强,等.基于有限元的电力变压器绕组弹塑性变形分析[J].高压电器,2019,55(6):170-176.

[14]

胡宇达,邱家俊.大型发电机定子端部绕组的电磁特性[J].机械强度,2006,28(1):20-24.

[15]

方建国,吴梦艳,朱珊珊,等.电机绝缘技术发展历程与方向[J].绝缘材料,2011,44(2):43-48.

[16]

万书亭,姚肖方,豆龙江.发电机定子绕组端部电磁力特性与鼻端扭矩计算[J].振动、测试与诊断,2014,34(5):920-925,980.

[17]

LIN

, LAIHO

A N

, HAAVISTO

, et al. End-winding vibrations caused by steady-state magnetic forces in an induction machine[J]. IEEE Transactions on Magnetics,2010,46(7):2665-2674.

2023年第56卷第1期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.16790/j.cnki.1009-9239.im.2023.01.008

接收时间：2022-01-11
首发时间：2025-11-21
出版时间：2023-01-20

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2022-01-11
修回日期：2022-03-23

基金

湖北省自然科学基金资助项目(2019CFC893)

作者信息

¹湖北文理学院纯电动汽车动力系统设计与测试湖北省重点实验室,湖北襄阳 441053

²新昇电气有限公司,湖北襄阳 441899

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/jycl/CN/10.16790/j.cnki.1009-9239.im.2023.01.008

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT