一种基于端口电压包络线的永磁电机偏心故障诊断方法

1.本发明涉及电机故障诊断技术领域，具体是一种基于端口电压包络线的永磁电机偏心故障诊断方法。


背景技术：

2.永磁电机具有高转速、高转矩、高可靠性而被广泛应用于航空航天、国防军工、汽车工业、风力发电等领域。作为机械能与电能之间进行转换装置，保证其稳定、高效运行至关重要。由于加工、装配以及长时间的高负荷运行，永磁电机可能会出现各种故障，偏心故障即为常见的电机故障的一种，若没有及时检测出电机偏心故障，长时间的偏心状态运行在电机内部产生不平衡磁拉力，使电机轴承磨损，故障程度进一步加剧，甚至引起转子损坏。研究表明，永磁电机偏心故障后，电机内部气隙分布不均匀，导致电机内部磁场分布发生变化。
3.因此，在电机内部放置探测线圈，根据法拉第电磁感应定律，不均匀分布的磁场在不同位置探测线圈中产生的感应电压不同，监测探测线圈端口电压变化，可以得知电机内部气隙变化情况，实现电机偏心故障诊断。为此，本发明提出一种基于端口电压包络线的永磁电机偏心故障诊断方法。


技术实现要素：

4.本发明旨在至少解决现有技术中存在的技术问题之一。为此，本发明提出一种基于端口电压包络线的永磁电机偏心故障诊断方法，实现了永磁电机偏心故障诊断与类型判别，能够有效的识别永磁电机的静态偏心、动态偏心及混合偏心故障；该方法易于实现，成本低。
5.为实现上述目的，根据本发明的第一方面的实施例提出一种基于端口电压包络线的永磁电机偏心故障诊断方法，包括如下步骤：
6.步骤一：将三个探测线圈采用漆包线，绕制在待诊断永磁电定子齿上；其中，每个探测线圈相隔整数对极；
7.步骤二：电机平稳运行时，通过电压采样电路采集每个探测线圈端口电压信号ui(x)；
8.步骤三：根据电压信号计算得到静态偏心故障指标se和动态偏心故障指标de；
9.步骤四：将静态偏心指标se和动态偏心指标de与设置好的阈值t
se
和t
de
进行比较；判断电机静态偏心和动态偏心故障。
10.进一步地，其中，静态偏心故障指标se和动态偏心故障指标de的具体计算步骤如下：
11.s31：首先对得到的三组探测线圈电压信号取绝对值；即：
12.u(i)(x)＝|u
sc(i)
(x)|
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
13.s32：对得到的绝对值信号取上包络线，即：
14.u
up(i)
(x)＝up(u(i)(x))
ꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
15.s33：对u
up(i)
(x)取倒数，即：
[0016][0017]
s34：利用信号均值对三组信号进行归一化，即
[0018][0019][0020]
s35：对归一化的信号三组信号w
n(i)
(x)分别取上包络线和下包络线，即：
[0021]vup(i)
(x)＝up(w
n(i)
(x))
ꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0022]vlow(i)
(x)＝low(w
n(i)
(x))
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0023]
s36：计算每组上、下包络线的均值，即：
[0024][0025]
s37：取每组信号上、下包络线均值的方差作为静态偏心存在性指标，即：
[0026]
se＝s(v(i))
ꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0027]
s38：取三组信号上、下包络线距离的均值作为动态偏心存在性指标，即：
[0028][0029]
进一步地，其中，步骤四中判断电机静态偏心和动态偏心故障，具体为：
[0030]
若se≥t
se
，则存在静态偏心；若de≥t
de
，则存在动态偏心；若se≥t
se
和de≥t
de
同时成立，则电机存在混合偏心。
[0031]
与现有技术相比，本发明的有益效果是：
[0032]
本发明在电机内部放置探测线圈，根据法拉第电磁感应定律，不均匀分布的磁场在不同位置探测线圈中产生的感应电压不同，监测探测线圈端口电压变化，可以得知电机内部气隙变化情况，实现电机偏心故障诊断，能够有效的识别永磁电机的静态偏心、动态偏心及混合偏心故障；该方法易于实现，成本低。
附图说明
[0033]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0034]
图1为本发明中探测线圈的摆放示意图。
[0035]
图2为本发明的原理框图。
具体实施方式
[0036]
下面将结合实施例对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
[0037]
其中图1为具体的36槽8极永磁电机模型与探测线圈在槽内的具体摆放位置，可以看到每个探测线圈相隔整数对极；其中探测线圈放置在定子槽内，其需要在电机出厂时直接埋入。在电机正常下线之后，在槽口嵌入探测线圈之后，再插入槽楔封住槽口；
[0038]
如图2所示，一种基于端口电压包络线的永磁电机偏心故障诊断方法，包括如下步骤：
[0039]
步骤一：将三个探测线圈采用漆包线，绕制在待诊断永磁电定子齿上；其中，每个探测线圈相隔整数对极；
[0040]
步骤二：电机平稳运行时，通过电压采样电路采集每个探测线圈端口电压信号ui(x)；
[0041]
步骤三：根据电压信号计算得到静态偏心故障指标se和动态偏心故障指标de；具体计算步骤如下：
[0042]
s31：首先对得到的三组探测线圈电压信号取绝对值；即：
[0043]
u(i)(x)＝|u
sc(i)
(x)|
ꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0044]
s32：对得到的绝对值信号取上包络线，即：
[0045]uup(i)
(x)＝up(u(i)(x))
ꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0046]
s33：对u
up(i)
(x)取倒数，即：
[0047][0048]
s34：利用信号均值对三组信号进行归一化，即
[0049][0050][0051]
s35：对归一化的信号三组信号w
n(i)
(x)分别取上包络线和下包络线，即：
[0052]vup(i)
(x)＝up(w
n(i)
(x))
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0053]vlow(i)
(x)＝low(w
n(i)
(x))
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0054]
s36：计算每组上、下包络线的均值，即：
[0055][0056]
s37：取每组信号上、下包络线均值的方差作为静态偏心存在性指标：
[0057]
se＝s(v(i))
ꢀꢀꢀꢀꢀꢀꢀ
(9)
[0058]
s38：取三组信号上、下包络线距离的均值作为动态偏心存在性指标：
[0059][0060]
步骤四：将静态偏心指标se和动态偏心指标de与设置好的阈值t
se
和t
de
进行比较；判断电机静态偏心和动态偏心故障，具体为：
[0061]
若se≥t
se
，则存在静态偏心；若de≥t
de
，则存在动态偏心；若se≥t
se
和de≥t
de
同时成立，则电机存在混合偏心。
[0062]
在本实施例中，探测线圈原理如下：
[0063]
偏心故障下任意角度气隙长度为
[0064][0065]
根据等效磁路欧姆定律，第j条等效磁路磁通为
[0066][0067]
根据以上分析可知，当电机发生偏心故障时，等效磁路中的磁通发生变化，为了检测磁通变化，本文提出采用在电机定子齿上紧密地绕制三个探测线圈的方式，进行偏心检测，每个探测线圈相隔整数极。通过探测线圈的磁通可以近似看作为该齿所在等效磁路的磁通，由电磁感应定律，探测线圈两端电压为：
[0068][0069]
取这些时间点可近似的看作为第j个定子齿上探测线圈取得最大电动势的时间点，那么一个机械周期t内将会有p个这样的时间点
[0070][0071]
[0072][0073]
根据单元电机理论，定子槽数为z，极对数为p的电机可以视为由p个相同的单元电机构成，当电机发生偏心故障时，可以视为p个相同的单元电机仅气隙分布不同，绕组磁动势与磁阻rc分布完全相同，因此相隔整数对极的定子齿上安装探测线圈，即相当于每个探测线圈安装在不同单元电机的相同位置上；即：
[0074]rc(j(1))
＝r
c(j(2))
＝r
c(j(3))
＝...＝r
c(j(p))
[0075][0076]
为了更方便地计算故障程度，对u
max
(j,k)和u
min
(j,k)取倒数，并将其定义为故障指示量
[0077]
当电机无故障时
[0078][0079][0080]
定义得到
[0081]
结论：根据分析可得，当转速恒定时，w
max
(i,k)和w
min
(i,k)为直流量，且绝对值相等；
[0082]
当电机只存在静态偏心时
[0083][0084][0085]
当电机只存在动态偏心时
[0086]
[0087][0088]
结论：动态偏心时，三个探测线圈故障指示量w
max
(i,k)和w
min
(i,k)可视为直流量与余弦量的叠加，直流分量幅值相等，余弦分量的周期为机械周期t，时间相差2π
·
i/p，幅值相等，且幅值可反映动态偏心程度；
[0089]
当电机混合偏心状态时
[0090][0091][0092]
结论：混合偏心状态下故障指示量w
max
(i,k)和w
min
(i,k)相当于静态偏心和动态偏心的线性叠加，静态偏心导致三个故障指示量直流分量幅值不同，动态偏心产生余弦分量；
[0093]
其中，不同偏心故障电机诊断结果如表1所示：
[0094]
静态偏心动态偏心sede诊断结果007.55*10-4
1.07*10-2
健康18％02.56*10-2
2.20*10-3
静态偏心27％04.11*10-2
2.30*10-3
静态偏心36％06.30*10-2
2.40*10-3
静态偏心45％08.86*10-2
2.30*10-3
静态偏心54％01.06*10-1
6.20*10-3
静态偏心63％01.07*10-1
8.80*10-3
静态偏心036％4.90*10-3
3.03*10-1
动态偏心045％6.10*10-3
3.67*10-1
动态偏心054％8.30*10-3
4.37*10-1
动态偏心063％1.00*10-2
4.90*10-1
动态偏心18％18％3.19*10-2
1.58*10-1
混合偏心27％27％4.86*10-2
2.37*10-1
混合偏心36％18％3.65*10-2
3.01*10-1
混合偏心36％36％6.55*10-2
3.07*10-1
混合偏心
[0095]
表1
[0096]
上述公式均是去除量纲取其数值计算，公式是由采集大量数据进行软件模拟得到最接近真实情况的一个公式，公式中的预设参数和预设阈值由本领域的技术人员根据实际情况设定或者大量数据模拟获得。
[0097]
在本说明书的描述中，参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施
例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
[0098]
以上公开的本发明优选实施例只是用于帮助阐述本发明。优选实施例并没有详尽叙述所有的细节，也不限制该发明仅为的具体实施方式。显然，根据本说明书的内容，可作很多的修改和变化。本说明书选取并具体描述这些实施例，是为了更好地解释本发明的原理和实际应用，从而使所属技术领域技术人员能很好地理解和利用本发明。本发明仅受权利要求书及其全部范围和等效物的限制。