基于线圈子单元的永磁同步电机健康与定子绕组短路故障数学模型

高彩霞1 苗 壮1 陈 昊2 司纪凯3 吕 珂4

(1.河南理工大学电气工程与自动化学院 焦作 454003 2.河南理工大学应急管理学院 焦作 454003 3.郑州大学电气工程学院 郑州 450001 4.海军工程大学舰船综合电力技术国防科技重点实验室 武汉 430033)

摘 要 针对相绕组为基本单元的永磁同步电机定子绕组短路故障数学模型,存在故障空间位置影响无法计及、不同短路故障需要重新建模,且无法研究健康状态时线圈对电机性能的贡献等问题,该文提出一种基于线圈子单元的永磁同步电机健康与定子绕组短路故障数学模型(APM)。首先,将每个线圈分割为多个线圈子单元,分别建立各子单元的电压、阻抗和空载反电动势等物理量的矩阵方程,利用基尔霍夫电压定律建立所有子单元的电压方程,建立电机的电磁功率和转矩方程。其次,利用有限元法和拟合法建立考虑槽内空间位置关系的线匝电感矩阵,通过矩阵变换得到计及位置信息的子单元电感矩阵。然后,在Matlab/Simulink 中建立APM 的仿真模型和图形化界面,通过修改抽头编号及图形化界面中短路电阻模块的连接位置,计算电机健康和不同类型定子绕组短路故障下的电磁特性。最后,有限元仿真和实验结果验证了模型的正确性和准确性,表明APM 不仅可以精细地分析健康状态时线圈对电机性能的贡献,也可以高效、准确地分析线圈内部不同位置短路故障对电机性能的影响。该文研究可以为永磁同步电机的设计、故障诊断和容错控制提供依据。

关键词:永磁同步电机 数学模型 线圈子单元 定子绕组短路故障 矩阵方程

0 引言

近年来,永磁同步电机(Permanent Magnet Synchronous Motor,PMSM)因其高转矩/惯量比、高功率密度、高效率、响应快等优点,被广泛应用于电动汽车、轨道交通、航空航天等领域[1-3]。PMSM在长期运行过程中,受到热应力、电应力、机械应力等因素影响[4-5],定子绕组容易因绝缘损坏而引发匝间短路故障(Interturn Short-circuit Fault,ISF)[6-7]。如果早期ISF 没有被及时监测到并采取相应措施,会迅速发展成多线圈短路故障、相间短路故障等更严重的故障,甚至可能引发灾难性事故[8-9]。定子绕组短路故障(Winding Short-circuit Fault,WSF)的精准诊断是提高电机可靠性、安全性和经济性的重要手段。建立模型对电机故障前后的性能分析并遴选故障特征量,是故障诊断的基础。目前,大量文献通过建立电机WSF 数学模型,实现对电机健康和WSF 状态的分析[10-17],但它们不能有效考虑绕组结构与故障空间位置的影响,难以分析线圈对电机性能的贡献及线圈内部不同位置短路故障对电机性能的影响。因此,有必要建立一个考虑绕组结构与故障空间位置影响的PMSM 健康状态与WSF 的数学模型。

目前,已有许多学者对电机WSF 模型进行了广泛而深入的研究,并取得一些重要成果。文献[10-13]采用绕组分区法,把故障相绕组分为健康和故障部分,以相绕组为基本单元计算模型参数,并建立电机短路故障模型。对健康状态分析时,该类模型仅能分析相绕组的电压、电流等参数;对故障状态分析时,该类模型忽略了相内故障空间位置的影响,无法区分相内不同位置ISF 对电机性能的影响。文献[14-16]从单个线圈出发计算电感和空载反电动势参数,再通过叠加得到实际回路的电感和空载反电动势参数,并根据电机的实际回路建立多回路模型。文献[14-16]不仅分析了健康状态时线圈对电机性能的贡献,也分析了支路内部不同位置ISF 对电机性能的影响。文献[17]基于分区法将相绕组分成多个线圈,并以线圈为基本单元,在Matlab/Simulink 中建立和求解模型。对健康状态分析时,该模型仍把相绕组作为一个整体,计算线圈的空载反电动势,忽略了线圈位置对空载反电动势的影响,难以准确分析线圈对电机性能的贡献;对故障状态分析时,该模型能够方便、快速地分析不同类型的WSF,但仍忽略了线圈内部故障线匝位置的影响。文献[18-21]建立了电机 WSF 的有限元模型(Finite Element Model,FEM),可以精确分析健康的及不同类型的WSF 下的电磁特性。FEM 考虑了故障空间位置、绕组结构等因素的影响,不仅能分析健康状态时线圈对电机性能的贡献,还能区分线圈内部不同位置ISF 对电机性能的影响。但FEM 对不同短路故障分析需要重新建模,而且其求解过程耗时。

为了在较短的时间内计算电机健康的和不同类型定子绕组短路故障下的电磁特性,并考虑电机绕组结构与故障空间位置影响,本文将每个线圈分割为多个线圈子单元,以线圈子单元为基本单元,建立一个基于线圈子单元的PMSM 健康与WSF 数学模型(A PMSM Mathematical Model,APM)。

本文的主要贡献为:①在Matlab/Simulink 中建立APM 的仿真模型和图形化界面,通过修改抽头编号及图形化界面中短路电阻模块的连接位置,方便、快捷地计算电机健康和不同类型定子绕组短路故障下的电磁特性。②考虑了电机绕组结构的影响,APM 能够准确分析健康状态时线圈对电机性能的影响。③考虑了故障空间位置的影响,APM 能够准确分析线圈内部不同位置短路故障对电机性能的影响。

1 基于线圈子单元的PMSM 数学建模

1.1 基于线圈子单元的PMSM 数学模型

为了高效、准确、精细化分析电机在健康及不同类型WSF 状态下的电磁性能,本文建立一个基于线圈子单元的PMSM 健康与WSF 数学模型,其电路示意图如图1 所示。

图1 PMSM 电路示意图
Fig.1 The circuit diagram of PMSM

该APM 将每个线圈分为三个线圈子单元,并引出四个抽头intXkjz 如图1a 所示,其中intXkj1 和intXkj4 为固定抽头,intXkj2 和intXkj3 为可变抽头。子单元的匝数随可变抽头位置的变化而变化,可变抽头的位置取决于待分析的短路故障的位置和匝数。根据上述划分规则,基于线圈子单元的等效电路如图1b 所示,其中每个子单元具有各自的电阻、电感及空载反电动势,在电机健康状态下,由同一线圈分割而成的三个子单元的空载反电动势、电感及电阻之和分别等于未分割线圈的空载反电动势、电感及电阻。

为简化分析,忽略涡流、磁滞损耗以及趋肤效应;电机铁心磁导率为无穷大;电机磁路为线性的。基于以上假设,该APM 的表达式为

式中,VsRsIsLsE0 分别为电压矩阵、电阻矩阵、电流矩阵、电感矩阵和空载反电动势矩阵。vXkjwiXkjwrXkjweXkjwLXkjw 分别为子单元 Xkjw的瞬时电压、瞬时电流、电阻、瞬时空载反电动势、瞬时自感。MXkjwYmno 为子单元Xkjw 与子单元Ymno之间的瞬时互感。PeP1pCupFe 分别为瞬时电磁功率、输入功率、定子铜耗和铁耗(根据建模假设,铁耗忽略不计)。Pe.acTave、f、nΩ 分别为电磁功率平均值、电磁转矩、电源频率、转速和机械角速度。其中XY 代表相编号,km 代表支路编号,jn 代表线圈编号,wo 代表子单元编号,z 代表抽头序号(X,Y=A,B,C,…;km=1,2,3,…;j,n=1,2,3,…;w,o=1,2,3;z=1,2,3,4;XkjwYmno)。基于式(1)~式(12),在Matlab/Simulink环境下的图形化界面中建立APM 仿真模型,如图2所示。

图2 建立的APM 仿真模型
Fig.2 The established APM simulation model

图2 所示的APM 仿真模型由交流电源AC、子单元电阻、子单元电感、短路电阻Rf、电流传感器、电压传感器等模块构成,其中IXkIfVXIX 分别为测量Xk 支路电流、短路电流、X 相电压、X 相电流的传感器端口。在分析电机健康状态性能时,只需要将图2 的Rf 开路,通过电源模块输入交流电源的幅值及频率设置运行工况,通过传感器模块获取相绕组和各线圈的电压、电流等参数,进而分析不同位置线圈对电机性能的贡献。对WSF 分析的步骤如下:

(1)设置运行工况:通过交流电源的幅值及频率模拟电机的负载及转速情况。

(2)设置不同的WSF:将Rf 并联到图2 所示的相应的抽头之间,完成WSF 的短路电阻、短路匝数及短路位置的设置,具体设置方式如下:

1)模拟槽口处发生ISF(短路匝数Nf<线圈匝数Nc),将Rf 连接到抽头intXkj1 和intXkj2 之间。

2)模拟槽底处发生ISF(NfNc),将Rf 连接到抽头intXkj3 和intXkj4 之间。

3)模拟单线圈短路时(Nf=Nc),将Rf 连接到抽头intXkjz 和intXknzjk)之间。

4)模拟相间短路故障,将Rf 连接到抽头intXkjz和intYkjzXY)之间。

(3)设置intXkj2 和intXkj3 的编号;沿槽深位置依次对定子线圈的线匝进行编号,如图3 所示。

图3 线匝位置编号
Fig.3 The position number of coil turns

根据待分析的短路类型、位置及匝数输入可变抽头intXkj2 和intXkj3 的编号。子单元的匝数由intXkj2和intXkj3 的编号自动确定,其电阻rXkjz、瞬时空载反电动势ec 和瞬时电感分别由式(13)、式(15)、和式(20)自动计算。

式中,rc 为单个线圈的电阻。ec 为单个线圈的瞬时空载反电动势矩阵。ep 为编号为p 的线匝的瞬时空载反电动势。C 为位置矩阵,位置矩阵第1 行第1列至第intXkj2 列数值为1,其余列为0;第2 行第intXkj2+1 列至第intXkj3 列数值为1,其余列为0;第3 行第intXkj3+1 列至第intXkj4 列数值为1,其余列为0。使用时输入可变抽头intXkj2 和intXkj3 的位置编号,可自动生成位置矩阵。eXk1w-hαh 为子单元 Xk1w的瞬时空载反电动势第h 次谐波分量的幅值和初相角。β 为相邻线圈在空间上相隔的电角度。结构正常的永磁体(不考虑退磁故障)产生空载反电动势的空间周期为1 对极,空载反电动势只包含基波和3、5 等奇数次谐波[15],即h=1,3,5,…。

(4)通过Simulink 环境下的传感器模块获取短路电流、支路电流、子单元电压、线圈电压、相电压等参数,并通过式(7)~式(12)计算功率及转矩。

在分析不同类型定子绕组短路故障时,传统数学建模方法需要改变模型的拓扑结构并重新建立方程;所提出的APM 只需要根据短路类型、位置及匝数在仿真模型的图形化界面中修改抽头编号及短路电阻模块的连接位置,具有方便、快速的优点。

1.2 考虑空间位置的电感计算

由于相同槽中不同位置线匝交链的槽漏磁通不同,所以相同槽、相同匝数、不同位置线圈的电感和空载反电动势也不同。对于图3 中线圈匝数为Nc的电机,对应的线匝电感计算为

式中,Lp 为位置编号为p 的线匝的瞬时自感;M pp为同一线圈内编号为p 的线匝与编号为p'的线匝之间的瞬时互感;M pp为线圈Xkj 内编号为p 的线匝与线圈Ymn 内编号为p'的线匝之间的瞬时互感;ψ p为位置编号为p 的线匝和永磁体共同在编号为p的线匝上产生的瞬时磁链; ψ p为位置编号为p 的线匝和永磁体共同在编号为p'的线匝上产生的瞬时磁链; ψ p为位置编号为p 的线匝和永磁体共同在编号为p'的线匝上产生的瞬时磁链;ψPMpψPMpψ PMp为永磁体分别在位置编号为p p'和p'的线匝上产生的瞬时磁链;i 为通入位置编号为p 的线匝的直流电流。

在线圈电感的有限元计算模型中,将1 A 的直流电通入到编号为 p 的线匝来计算 LpM ppM pp。为了减小电感计算的工作量,首先利用FEM计算出奇数线匝瞬时电感,其次利用拟合法计算出剩余线匝瞬时电感,即

最后通过变换矩阵计算出所有子单元电感为

2 有限元仿真验证

2.1 PMSM 关键参数

为了验证APM 的正确性和准确性,本文以一台电流源供电的 66 极72 槽星形联结的三相面贴式PMSM 为研究对象。其定子绕组采用分数槽集中绕组隔齿绕的型式。PMSM 的关键参数见表1。

表1 PMSM 的关键参数
Tab.1 The key parameters of PMSM

2.2 健康状态的解析与有限元结果比较分析

为了验证所建立的APM 在健康状态时不同工况下对电机性能分析的正确性与准确性,对由两种转速情况(n=100,200 r/min)和三种负载情况(IX=14,28,42 A)组成的六种工况下的电机性能进行比较分析。进行对比的六种工况列于表2 中。

表2 进行对比的工况
Tab.2 The compared conditions

表3 为不同工况下APM 与FEM 计算出的平均转矩及其误差。

表3 平均转矩
Tab.3 The average torque

由表3 可知,APM 计算出的平均转矩略微大于FEM 计算出的平均转矩,最大误差约为4.85 %,该误差是由于建模时忽略铁耗和磁路饱和引起的。

表4 为不同工况下APM 与FEM 计算出的A 相电压VA 峰峰值及其误差。

表4 VA 的峰峰值
Tab.4 The peak-peak value of VA

由表4 可知,APM 计算出的VA 的峰峰值略微小于FEM 计算出的VA 的峰峰值,最大误差约为2.4 %,该误差主要是由于建模时忽略磁路饱和引起的。

表5 为不同工况下APM 与FEM 计算出的Xk1线圈、Xk2 线圈电压VXk1VXk2 峰峰值及其误差。

由表5 可知,APM 与FEM 计算出的线圈电压峰峰值的最大误差约为2.2 %。Ak2 线圈的电压高于Ak1 线圈的电压,说明电机运行时Xk2 线圈比Xk1线圈产生的功率多,这是由于两类线圈空载反电动势相位不同造成的。因此,Xkjj=2、4)和Xknn=1、3)线圈可以设计成不同的绝缘等级,提高使用寿命,降低制造成本。与工况Ⅳ的VAk1 相比,工况Ⅴ的VAk1和工况Ⅵ的VAk1 分别增加了1.9 %和7.6 %,与工况Ⅳ的VAk2 相比,工况Ⅴ的VAk2 和工况Ⅵ的VAk2 分别增加了10.4 %和22.3 %,Xk2 线圈比Xk1 线圈对负载变化更敏感,说明电机运行时Xk2 线圈比Xk1 线圈承担的负载多。在恒定的转矩负载下,Xk2 线圈故障时对电机性能的影响大于 Xk1 线圈故障时。当APM 仿真模型完成建模后,对健康状态分析时,单次仿真计算需要2.3 s,而相同条件下FEM 的单次仿真计算需要4 141 s,APM 可以快速地对电机健康状态进行分析。

表5 VXk1VXk2 的峰峰值
Tab.5 The peak-peak value of VXk1,VXk2

因此,有限元仿真验证了所提出的APM 不仅能精确、高效地分析健康状态时电机在不同工况下的相电压、线圈电压、平均转矩等电磁特性,还能精细地分析线圈对电机性能的贡献。利用APM 对不同位置的线圈进行分析,能够发现同一支路不同位置线圈的电压不同,功率不同,承担的负载也不同。在设计线圈绝缘时,为降低电机制造成本,同一支路不同位置线圈采用不同的绝缘等级。上述分析为PMSM 的本体设计提供依据。由于同一支路不同线圈的电压不同,在线圈绝缘等级相同的情况下,同一支路不同线圈的故障概率不同;线圈内部同一短路匝数不同短路位置故障对电机性能的影响也不相同。上述分析可为PMSM 的故障诊断提供依据。

2.3 WSF 状态的解析与有限元结果比较分析

为了验证APM 在不同WSF 下特性分析的正确性与准确性,对三类WSF 下的电机电磁特性进行比较分析。仿真分析的故障类型列于表6 中,电机运行在额定工况下(n=200 r/min,IX=28 A)。

表6 仿真分析的故障类型
Tab.6 Fault types of simulation analysis

A12_01~48 表示A12 线圈48 匝短路故障,为分析该故障,将Rf 并联在A12 线圈的抽头intA121和intA124 之间。A11_01~24 和A11_25~48 分别为A11 线圈槽口处和槽底处24 匝短路,槽口24 匝短路时将抽头intA112 的值设置为24,并将Rf 连接在抽头intA111 和intA112 之间;槽底24 匝短路时将抽头intA113 的值设置为25,并将Rf 连接在抽头intA113 和intA114 之间。A34.48_C34.48 表示A34 线圈与C34线圈间发生相间短路故障,为分析该故障,将 Rf并联在A34 线圈的抽头intA341 和C34 线圈的抽头intC341 之间。

图4 所示为 PMSM 分别在 A12_01~48、A11_01~24 及A11_25~48 故障状态下由APM 与FEM 计算出的VAIA1If 和平均转矩TaveRf 的变化曲线。

图4 不同故障下APM 结果与FEM 结果比较
Fig.4 Comparison between APM results and FEM results under different faults

由图4a 和图4b 可知,APM 与FEM 计算结果吻合度较好,A12_01~48 故障情况下,最大误差约为3.0 %。该误差是由于建模时忽略磁路饱和与铁耗引起的。有限元仿真证明了APM 能够分析单线圈ISF,并且具有较高的精度。

由图4c~图4f 可知,APM 与FEM 计算结果吻合度较好,A11_01~24 故障情况下,最大误差为5.7 %;A11_25~48 故障情况下,最大误差为5.5 %。该误差是由于建模时忽略磁路饱和与铁耗引起的。有限元仿真证明了APM 能够分析线圈内部ISF,并且具有较高的精度。

图5 为PMSM 在A34.48_C34.48 故障状态下由APM 与FEM 计算出的VA、C 相电压VC、A3 支路电流IA3、C3 支路电流IC3IfTaveRf 的变化曲线。

图5 A34.48_C34.48 下APM 结果与FEM 结果比较
Fig.5 Comparison between APM results and FEM results under A34.48_C34.48

由图5 可知,APM 与FEM 计算结果吻合度较好,A34.48_C34.48 故障情况下,最大误差约为4.0 %。该误差是由于建模时忽略磁路饱和与铁耗引起的。有限元仿真验证了APM 也适用于分析相间短路故障,并且具有较高的精度。

图6 为A11 线圈内部不同位置发生24 匝金属性短路时的If

图6 不同位置发生24 匝金属性短路时的If
Fig.6 If of 24 turns metallic short-circuit at different positions

由图6 可知,槽口24 匝金属性短路时的If 幅值远大于槽底24 匝金属性短路时的If 幅值。这是由于槽内相同匝数不同位置线圈的电感沿槽深方向依次增大造成的。对比分析图4c~图4f 和图6 可知,在NfRf 相同的前提下,槽口故障的If 大于槽底故障的If;槽口故障的VATave 小于槽底故障的VATave,说明槽口故障对电机性能的影响大于槽底故障。可以看出,该APM 能够精确分析线圈内部不同位置ISF 对电机性能的影响。当APM 仿真模型完成建模后,对WSF 状态分析时,单次仿真计算需要10.8 s,而相同条件下FEM 的单次仿真计算需要4 323 s,APM 可以快速分析不同类型WSF。

APM 与FEM 的对比分析表明,该APM 不仅可以高效、精确地分析电机在不同类型WSF 下的相电压、支路电流、短路电流、平均转矩等电磁特性,还可以精细地分析槽内线圈不同位置的ISF 对电机性能的影响。利用APM 对槽内线圈不同位置ISF进行分析,发现在NfRf 相同的前提下,槽内不同位置ISF 下的短路电流和平均转矩不同,揭示了NfRf 作为故障程度判据的局限性。因此,为保证电机在槽内不同位置ISF 状态下仍能输出满足要求的转矩,其抑制短路电流的补偿策略也不相同。上述分析为PMSM 故障诊断及容错控制提供了依据。

3 实验验证

为了进一步验证所建立APM 的正确性和准确性,本文制造了2.1 节所述的PMSM 故障实验样机,并搭建了实验平台,如图7 所示。为了进行线圈内部不同位置的ISF 实验,将样机的每个线圈沿槽深方向横向分为四个线圈子单元,并将所有子单元的抽头引出,每个子单元的匝数为12 匝,如图7a 所示。实验平台由样机、转矩测量仪、磁粉制动器、工控机、驱动器、数字示波器、数字电桥、直流电源和短路电阻等构成。

图7 实验装置
Fig.7 Experiment setup

首先对健康状态下的APM 进行实验验证。图8为额定工况下电机的线圈电压VA11VA12 及转矩的APM 结果与实验结果。

由图8 可知,VA11VA12 和转矩的APM 结果与实验结果的吻合程度较好,VA11 的幅值明显小于VA12 的幅值,最大误差分别为3.4 %、3.1 %和1.8 %。该误差是由于建模时忽略铁耗、磁路饱和、机械摩擦等引起的。实验结果验证了APM 不仅能精确计算电机健康状态下的电磁特性,还能分析不同位置的线圈对电机性能的贡献。

图8 额定工况下电机的线圈电压和转矩
Fig.8 Coil voltage and torque of motor under rated condition

在进行电机WSF 状态实验测试时,将Rf 连接至相应子单元的抽头之间。另外,为了保证实验设备的安全,Rf 设置为2 Ω。图9 分别为额定工况时A12_01~48、A11_01~24、A11_25~48 和A34.48_C34.48 故障情况下If 的APM 结果与实验结果。由图9 可知,If 的APM 结果与实验结果的吻合程度较好,四种故障情况下最大误差分别为5.6 %、8.9 %、8.4 %及4.9 %。

图9 不同故障类型下的If
Fig.9 If under different fault types

图10 分别为额定工况时A12_01~48、A11_01~24、A11_25~48 和A34.48_C34.48 故障情况下故障支路电流IA1IA3IC3 的APM 结果与实验结果。

图10 不同故障类型下的故障支路电流
Fig.10 Fault branch current under different fault types

由图10 可知,APM 结果与实验结果的吻合程度较好,A12_01~48、A11_01~24 和A11_25~48故障情况下 IA1 最大误差分别为 4.9 %、3.4 %及3.7 %。A34.48_C34.48 故障情况下IA3IC3 的最大误差分别为4.2 %、4.7 %。

图11 分别为25 %负载、100 %转速和25 %负载、50 %转速两种运行工况时A12_01~48 故障情况下IfIA1 的APM 结果与实验结果。

由图11 可知,APM 结果与实验结果的吻合程度较好,25 %负载100 %转速与25 %负载50 %转速两种运行状态下,If 的最大误差分别为6.2 %及4.6 %;IA1 的最大误差分别为4.9 %及4.8 %。上述误差是由于建模时忽略铁耗、磁路饱和、短路电阻温升和机械摩擦引起的。实验结果验证了APM 在分析电机WSF 状态下电磁特性的正确性和准确性。实验结果验证了APM 不仅能精确计算电机WSF 状态下的电磁特性,还能分析线圈内部不同位置短路故障对电机性能的影响。

图11 两种运行工况下的IfIA1
Fig.11 If and torque under two operating conditions

4 结论

本文提出了一种基于线圈子单元的PMSM 健康与WSF 数学模型,并在Matlab/Simulink 中建立了其仿真模型和图形化界面。该模型考虑了绕组结构与短路空间位置的影响,可以分析健康状态时线圈对电机性能的贡献及线圈内部不同位置短路故障对电机性能的影响。通过将模型计算结果、有限元仿真结果和实验测量结果进行对比分析,验证了模型的正确性和准确性,并得到以下结论:

1)APM 通过修改抽头编号及图形化界面中短路电阻模块的连接位置,可以方便、快捷地计算电机健康和不同类型定子绕组短路故障下的电压、电流、转矩等电磁特性。

2)APM 可以精细地分析PMSM 健康状态时不同工况下线圈对电机性能的贡献。研究发现,PMSM绕组内不同位置的线圈对电机性能的贡献不同。

3)APM 在不改变内部结构的前提下,可以高效、准确、精细地分析ISF、相间短路故障等不同类型的WSF 下的电机性能。

4)APM 可以分析同一槽相同短路匝数槽口处短路故障与槽底处短路故障对电机性能的影响。研究发现,当短路匝数和短路电阻都相同时,PMSM 槽内线圈不同位置的短路故障对电机性能影响不同。

所提出的数学模型可为PMSM 的设计、故障诊断和容错控制奠定重要基础。

参考文献

[1] 李晓华,赵容健,田晓彤,等.逆变器供电对电动汽车内置式永磁同步电机振动噪声特性影响研究[J].电工技术学报,2020,35(21): 4455-4464.

Li Xiaohua,Zhao Rongjian,Tian Xiaotong,et al.Study on vibration and noise characteristics of interior permanent magnet synchronous machine for electric vehicles by inverter[J].Transactions of China Electrotechnical Society,2020,35(21): 4455-4464.

[2] 张伟伟,肖飞,刘计龙,等.轨道交通车辆永磁同步牵引系统断电区穿越控制策略[J].电工技术学报,2021,36(16): 3483-3492.

Zhang Weiwei,Xiao Fei,Liu Jilong,et al.Power-off area traversing control strategy of permanent magnet synchronous motor traction system in rail transit vehicle[J].Transactions of China Electrotechnical Society,2021,36(16): 3483-3492.

[3] 崔淑梅,匡志,杜博超,等.基于自抗扰控制原理的全电飞机用永磁同步电机转速闭环控制[J].电工技术学报,2017,32(增刊1): 107-115.

Cui Shumei,Kuang Zhi,Du Bochao,et al.Speed closed-loop control of permanent magnet synchronous motor for all-electric aircraft applications based on active disturbance rejection controller[J].Transactions of China Electrotechnical Society,2017,32(S1): 107-115.

[4] 郑大勇,张品佳.交流电机定子绝缘故障诊断与在线监测技术综述[J].中国电机工程学报,2019,39(2): 395-406,637.

Zheng Dayong,Zhang Pinjia.A review of fault diagnosis and online condition monitoring of stator insulation in AC electrical machine[J].Proceedings of the CSEE,2019,39(2): 395-406,637.

[5] Zafarani M,Bostanci E,Qi Yuan,et al.Interturn short-circuit faults in permanent magnet synchronous machines: an extended review and comprehensive analysis[J].IEEE Journal of Emerging and Selected Topics in Power Electronics,2018,6(4): 2173-2191.

[6] Hang Jun,Ding Shichuan,Zhang Jianzhong,et al.Detection of interturn short-circuit fault for PMSM with simple fault indicator[J].IEEE Transactions on Energy Conversion,2016,31(4): 1697-1699.

[7] 谢颖,胡圣明,陈鹏,等.永磁同步电机匝间短路故障温度场分析[J].电工技术学报,2022,37(2):322-331.

Xie Ying,Hu Shengming,Chen Peng,et al.Thermal field analysis on inter-turn short circuit fault of permanent magnet synchronous motor[J].Transactions of China Electrotechnical Society,2022,37(2):322-331.

[8] Qi Yuan,Bostanci E,Gurusamy V,et al.A comprehensive analysis of short-circuit current behavior in PMSM interturn short-circuit faults[J].IEEE Transactions on Power Electronics,2018,33(12):10784-10793.

[9] Choi G,Jahns T M.PM synchronous machine drive response to asymmetrical short-circuit faults[J].IEEE Transactions on Industry Applications,2016,52(3):2176-2185.

[10] Jeong I,Hyon B J,Nam K.Dynamic modeling and control for SPMSMs with internal turn short fault[J].IEEE Transactions on Power Electronics,2013,28(7):3495-3508.

[11] Hang Jun,Zhang Jianzhong,Cheng Ming,et al.Online interturn fault diagnosis of permanent magnet synchronous machine using zero-sequence components[J].IEEE Transactions on Power Electronics,2015,30(12): 6731-6741.

[12] Jeong H,Moon S,Kim S W.An early stage interturn fault diagnosis of PMSMs by using negative-sequence components[J].IEEE Transactions on Industrial Electronics,2017,64(7): 5701-5708.

[13] 张业成,刘国海,陈前.基于电流波动特征的永磁同步电机匝间短路与局部退磁故障分类诊断研究[J].电工技术学报,2022,37(7): 1634-1643,1653.

Zhang Yecheng,Liu Guohai,Chen Qian.Discrimination of interturn short-circuit and local demagnetization in permanent magnet synchronous motor based on current fluctuation characteristics[J].Transactions of China Electrotechnical Society,2022,37(7):1634-1643,1653.

[14] 赵洪森,戈宝军,陶大军,等.定子绕组匝间短路对发电机电磁转矩特性的影响[J].电工技术学报,2016,31(5): 192-198.

Zhao Hongsen,Ge Baojun,Tao Dajun,et al.Influence of stator winding inter-turn short-circuit fault on generator electromagnetic torque characteristics[J].Transactions of China Electrotechnical Society,2016,31(5): 192-198.

[15] 牛化敏,桂林,孙宇光,等.基于多回路理论的交流励磁电机定子绕组内部故障仿真与实验研究[J].中国电机工程学报,2019,39(12): 3676-3684.

Niu Huamin,Gui Lin,Sun Yuguang,et al.Simulation and experimental research on stator winding internal faults of AC excitation machine based on multi-loop theory[J].Proceedings of the CSEE,2019,39(12):3676-3684.

[16] 田代宗,孙宇光,王善铭,等.多相整流永磁同步发电机绕组内部相间短路的故障分析[J].电工技术学报,2020,35(6): 1262-1271.

Tian Daizong,Sun Yuguang,Wang Shanming,et al.Analysis of stator internal phase-to-phase shortcircuit in the multiphase permanent magnet synchronous generator with rectifier load system[J].Transactions of China Electrotechnical Society,2020,35(6): 1262-1271.

[17] Ben Khader Bouzid M,Champenois G,Maalaoui A,et al.Efficient simplified physical faulty model of a permanent magnet synchronous generator dedicated to the stator fault diagnosis part I: faulty model conception[J].IEEE Transactions on Industry Applications,2017,53(3): 2752-2761.

[18] Vaseghi B.Modelling and study of PM machines with inter-turn fault dynamic model-FEM model[J].Electric Power Systems Research,2011,81(8): 1715-1722.

[19] Kumar N P,Isha T B.FEM based electromagnetic signature analysis of winding inter-turn short-circuit fault in inverter fed induction motor[J].CES Transactions on Electrical Machines and Systems,2019,3(3):309-315.

[20] 张艳辉,郑晓钦,吴新振,等.基于有限元场路耦合的十二相整流发电机系统定子匝间短路分析[J].电工技术学报,2019,34(9): 1842-1849.

Zhang Yanhui,Zheng Xiaoqin,Wu Xinzhen,et al.Analysis of stator winding inter-turn short circuit on a 12-phase rectifier generator system based on finite element field-circuit coupling[J].Transactions of China Electrotechnical Society,2019,34(9): 1842-1849.

[21] Zhao Jing,Guan Xiaoqing,Li Chenghai,et al.Comprehensive evaluation of inter-turn short circuit faults in PMSM used for electric vehicles[J].IEEE Transactions on Intelligent Transportation Systems,2021,22(1): 611-621.

A Mathematical Model Based on Coil Sub-Element for Permanent Magnet Synchronous Motor with Health and Stator Winding Short-Circuit Fault

Gao Caixia1 Miao Zhuang1 Chen Hao2 Si Jikai3 Lü Ke4
(1.School of Electrical Engineering and Automation Henan Polytechnic University Jiaozuo 454003 China 2.School of Emergency Management Henan Polytechnic University Jiaozuo 454003 China 3.School of Electrical Engineering Zhengzhou University Zhengzhou 450001 China 4.National Key Laboratory for Vessel Integrated Power System Technology Naval University of Engineering Wuhan 430033 China)

Abstract To solve the issues existing in the mathematical model of stator winding short-circuit fault of permanent magnet synchronous motor (PMSM) with phase winding as the basic element,such as the influence of fault space location cannot be taken into account,different short-circuit faults need to be remodeled,and the contribution of coil to motor performance under healthy condition cannot be studied.This paper proposes a mathematical model based on coil sub-element for PMSM (APM) with health and stator winding short-circuit fault.Firstly,each coil of PMSM is divided into multiple coil sub-elements,and the matrix equations of physical quantities such as voltage,impedance and no-load back EMF of each sub-element are established respectively.The voltage equations of all sub-elements are established by using Kirchhoff's voltage law,and the electromagnetic power and torque equations of the motor are established.Secondly,the finite element method is used to establish the no-load back EMF matrix of each wire turn considering the spatial position relationship in the slot,and the sub-element no-load back EMF matrix considering the position information is obtained by matrix transformation.Then,the finite element method and fitting method are used to establish the inductance matrix of each wire turn considering the spatial position relationship in the slot,and the sub-element inductance matrix considering the position information is also obtained by matrix transformation.Finally,the simulation model and graphical interface of APM are established in Matlab/Simulink.By modifying the tap number and the connection position of the short-circuit resistance module in the graphical interface,the electromagnetic characteristics of motor health and different types of stator winding short-circuit fault (WSF) are calculated.By comparing the APM calculation results,finite element simulation results and experimental results,the correctness and accuracy of the APM are verified,and the following conclusions are obtained: (1) By modifying the tap number and the connection position of the short-circuit resistance module in the graphical interface,APM can conveniently and quickly calculate the voltage,current,torque and other electromagnetic characteristics of the motor health and different types of WSF.(2) Considering the influence of the motor winding structure,APM can subtly analyze the contribution of coils to motor performance under different operating conditions when PMSM is in a healthy state.Using APM to analyze coils at different positions,it can be found that the coils at different positions of the same branch have different voltages,different powers,and different loads.Due to the different voltages of different coils in the same branch,the fault probability of different coils in the same branch is different under the same coil insulation level.In the design of coil insulation,coils at different positions of the same branch can adopt different insulation levels to reduce motor manufacturing costs.(3) Without changing the internal structure,APM can efficiently,accurately and subtly analyze the motor performance under different types of WSF such as interturn short-circuit fault and interphase short-circuit fault.(4) Considering the influence of fault space position,APM can accurately analyze the influence of short-circuit fault at different positions of the coil in the PMSM slot on motor performance.It is found that when the short-circuit turns (Rf) and short-circuit resistance (Nf) are the same,the short-circuit current and average torque are different under ISF at different positions in the slot,which reveals the limitation of Nf and Rf as the criterion of fault severity.

Keywords:Permanent magnet synchronous motor,mathematical model,coil sub-element,stator winding short-circuit fault,matrix equations

中图分类号:TM351

DOI: 10.19595/j.cnki.1000-6753.tces.211762

国家自然科学基金项目(52177039)、河南省科技攻关项目(222102220017)和河南理工大学博士基金项目(B2021-22.B2018-48)资助。

收稿日期 2021-11-01

改稿日期 2022-09-06

作者简介

高彩霞 女,1981 年生,教授,硕士生导师,研究方向为电机智能故障诊断与预测。

E-mail: gcx@hpu.edu.cn

陈 昊 男,1979 年生,副教授,硕士生导师,研究方向为机电一体化技术。

E-mail: chenhao@hpu.edu.cn(通信作者)

(编辑 郭丽军)