基于信息融合技术的故障诊断方法举例课件

1、7笼型异步电机故障的信息融合诊断方法 7.1 笼型异步电机故障 v定子铁芯故障铁芯短路、振动v绕组绝缘故障老化、磨损、过热、受潮、污染、电晕v转子故障转子气隙偏心、笼条断裂（断条）、端环开裂$常见故障转子偏心、断条、端环开裂和绕组过热（匝间短路） 7.1.1 转子断条和端环开裂故障 v电机起动时，转子绕组内短时间流过很大电流，定子绕组电流可达（57）In转子导条承受冲击电动力、热应力、离心应力 v起动、运行、停转过程中的循环热应力和变形作用 v由于负载变化和电压波动时，笼条又要受到交变负荷的作用，容易产生疲劳v转子绕组铸造质量、导条与端环的材质和焊接质量存在问题v故障现象起动时间延长、滑差加大

2、，力矩减小；振动和噪声增加，电流出现摆动等转子断条故障的特征频率 v定子电流中出现边频分量：式中，转差率s（n1-n）/n1100电机磁场旋转速度n160f1/p转子转速n（r/min）p为电机磁场极对数f1为电源基频k为正整数v其它检测方法振动、轴向磁通、转速波动等方法 1)21 (fksfbb7.1.2 气隙偏心故障 v气隙偏心:气隙磁场的不对称分布v气隙偏心的类型静态偏心偏心的位置在空间是固定的，与转子位置无关定子铁芯内径的椭圆度或装配不正确 动态偏心偏心位置在空间是变化的，通常与转子的位置和旋转频率有关转轴弯曲、轴颈椭圆、临界转速时的机械共振、轴承磨损 气隙偏心故障的特征频率 v在定子

3、电流频谱中 Z2为转子槽数，n rt为任意整数，n d为任意整数 （静态偏心：n d0；动态偏心：n d1、2、3、） n w s为奇整数v简化的气隙偏心故障特征频率判据 转子旋转频率fr（1s）f1/p，m为正整数 v气隙偏心故障还会造成电机的振动异常，并出现振动特征频率 12)/)1)(fnpsnZnfwsdrtagragmfff17.1.3 绕组过热与匝间短路故障 v绕组过热是电机运行过程中经常出现的一种现象，通常是短时和整体性的v比较危险的是绕组的局部过热绝缘的局部损伤匝间短路v温度监测是诊断电机绕组过热与匝间短路故障的常用方法当绕组过热导致绝缘局部老化，出现非金属性轻微匝间短路时，电

4、机的三相电流会出现不同的变化，一般导致有匝间短路的相电流增加7.1.4 单相电流谱分析方法的局限性 v是诊断鼠笼异步电机转子断条、端环开裂和气隙偏心等故障的常用方法v局限性故障特征频率易被基频的“旁瓣”所湮没故障严重程度比较轻微负载频繁波动容易受到电源谐波干扰（变频调速电源供电）未能充分发掘三相电流中所携带的故障特征信息新方法：信息融合 7.2 转子断条故障的诊断 v应用扩展派克矢量（应用扩展派克矢量（Park Vector）方法）方法v基于基于Park矢量的融合电流小波分析方法矢量的融合电流小波分析方法 7.2.1 应用扩展派克矢量方法 v只有转子断条或端环开裂故障时 )4cos(3)2co

5、s(3)2cos(3)(2/3(2222rldrdlrdrfldlfdrdlftsiitsiitsiiiiijii同时存在断条和气隙偏心故障时 )/ )1 (1cos()/ )1 (1cos( )21cos()21cos()cos(111111ecpmecnmecpmecpmmkbnkbnkbpkbpkapatpsmItpsmItksItksItIi3/2)/ )1 (1cos(3/2)/ )1 (1cos( 3/2)21cos(3/2)21cos()3/2cos(111111ecpmecnmecpmecpmmkbnkbnkbpkbpkbpbtpsmItpsmItksItksItIi3/2)/

6、 )1 (1cos( 3/2)/ )1 (1cos( 3/2)21cos( 3/2)21cos() 3/2cos(111111ecpmecnmecpmecpmmkbnkbnkbpkbpkcpctpsmItpsmItksItksItIi两种故障同时发生时派克矢量的模k=1、m=1 / )1 (2cos3)2/ )1cos(3)2/ )1cos(3)2/ )1cos(3)2/ )1cos(3)4cos(3)/ )1cos(3)/ )1cos(3)2cos(3)2cos(3)(2/3111111111111111111111111111111111111111111111121212121212ec

7、necpecnecpecpbnecpbnecnbpecnbpecnbnecnbnecpbpecpbpbnbpbnbpecnecnpecpecppbnbnpbpbppecnecpbnbpptpsIItspsIItspsIItspsIItspsIItsIItpsIItpsIItsIItsIIIIIIIjii派克矢量的模中包含的频率成分 序号频率成分产生原因1直流三相基波分量、各故障特征分量自身作用产生22ksf1基波分量与断条故障特征分量相互作用产生34ksf1两种断条故障特征分量之间相互作用产生4mfr基波分量与气隙偏心故障特征分量相互作用产生52 mfr气隙偏心故障特征分量相互作用产生6mfr

8、+2ksf1气隙偏心故障特征分量与断条故障特征分量相互作用产生7mfr-2ksf11 故障特征矢量的提取01020304050-3-2-10123f/HzdB01020304050-1-0.500.51f/Hz无转子断条故障的特征矢量无转子断条故障的特征矢量 01020304050-1-0.500.51f/Hz01020304050-2-10123f/HzdB转子断条故障的特征矢量转子断条故障的特征矢量 2 运用BP网络实现断条故障的自动识别 v输出层有输出层有4个神经元，分别对应电机的个神经元，分别对应电机的4种故障模式种故障模式1 1、良好、良好2 2、断条、断条3 3、气隙偏心、气隙偏心

9、4 4、匝间短路、匝间短路BP网络训练误差曲线 0 5 10 15 10 -4 10 -2 10 0 10 2 训练步数误差曲线 训练目标BP网络对断条故障的诊断结果 12345678910正常0.99910.90060.97780.99850.97290.11210.00070.00050.00120.0071断条0.00170.39840.13960.01470.18270.93360.99800.9989099800.9900气隙偏心000.00010.0001000.00020.000100.0001匝间短路0.00010.00090.00060.00060.010200000实际故障

10、正常正常正常正常正常断条断条断条断条断条7.2.2基于Park矢量的融合电流小波分析方法 v小波变换不同分解尺度对应的频带宽度设信号采样频率为fs，则信号中的最高频率为fs/2采用DWT对信号进行分解时，设分解尺度为j（j1，2，3，）j1时，信号被分解为0，fs/4和fs/4，fs/2两个频带j2时，信号被分解为0，fs/8、fs/8，fs/4和fs/4，fs/2三个频带依次类推，jk时，分解信号的最低频带为0，fs/2k+1 小波变换不同分解尺度对应的最低频带宽度（fs5000Hz） 尺度j123456780,fs/2j+10,12500,6250,312.50,156.30,78.10,

11、39.10,19.50,9.8良好电机空载运行时DB5小波分解结果 01002003004005001.61.822.22.42.62.83x 10-3ti*i020406080100-0.0100.010.020.030.04（a）Park矢量融合电流信号 (b)尺度j8时的低频带信号 良好电机满载运行时DB5小波分解结果 （a）Park矢量融合电流信号矢量融合电流信号 (b)尺度尺度j8时的低频带信号时的低频带信号 01002003004005000.0060.0070.0080.0090.010.0110.0120.0130.014020406080100-0.0500.050.10.1

12、50.22根断条故障电机满载运行时DB5小波分解结果 （a）Park矢量融合电流信号矢量融合电流信号 (b)尺度尺度j8时的低频带信号时的低频带信号 020406080100-0.0500.050.10.150.201002003004005000.0060.0070.0080.0090.010.0110.0120.0130.0145根断条故障电机空载运行时DB5小波分解结果 （a）Park矢量融合电流信号矢量融合电流信号 (b)尺度尺度j8时的低频带信号时的低频带信号 020406080100-0.0500.050.10.150.2010020030040050056789101112x 1

13、0-35根断条故障电机满载运行时DB5小波分解结果 （a）Park矢量融合电流信号矢量融合电流信号 (b)尺度尺度j8时的低频带信号时的低频带信号 020406080100-0.0500.050.10.150.2010020030040050056789101112x 10-3020406080100-0.0100.010.020.030.04020406080100-0.0500.050.10.150.2020406080100-0.0500.050.10.150.2020406080100-0.0200.020.040.060.080.1020406080100-0.0500.050.10

14、.150.2良好空载良好满载5根断条空载5根断条满载2根断条满载7.3绕组匝间短路故障的诊断 （a）良好状态 （b）断条 （c）气隙偏心 （d）匝间短路 （e）断相 （f）断条匝间短路 7.3.1 故障特征矢量的提取 00111111111100000000000000101000000000001111000000011100000001111111111110000001100000000000011000011110000000011000001100000000000001100001100000000000110001110000000000000011001100000000000

15、00100011000000000000000110011000000000000011001100000000000000011011000000000000000111100000000000000000111100000000000000011110000000000000000011100000000000000001110000000000000000001110000000000000000011000v数字滤波数字滤波, ,去除信号中的奇异点去除信号中的奇异点v将极坐标系中将极坐标系中360极角均匀分为极角均匀分为4040个等角扇形，个等角扇形，对每个扇形，按照等面积原则，分为对

16、每个扇形，按照等面积原则，分为1010等分等分v根据根据ParkPark矢量轨迹的覆盖情况，即可获得对应的矢量轨迹的覆盖情况，即可获得对应的黑白象素表示黑白象素表示（1 1和和0 0为元素的矩阵）为元素的矩阵）运用RBF网络实现故障的自动识别 RBF网络对匝间短路故障的识别结果 12345678910正常1.0001.0001.0001.0001.0000.0000.0001.0000.0000.000断条0.0000.0000.0000.0000.0000.0000.0000.0000.0000.000气隙偏心0.0000.0000.0000.0000.0000.0000.0000.0000

17、.0000.000匝间短路0.0000.0000.0000.0000.0001.0001.0000.0001.0001.000诊断结果正确正确正确正确正确正确正确错误正确正确匝间短路故障相别判定方法 良好转子顺时针旋转、空载运行时的Park矢量平均轨迹 （三角形接法） Park矢量轨迹倾斜角绕组匝间短路相别判据 （一） 当电机三相绕组为三角形接法，不论转子旋转方向为顺时针还是逆时针，也不论电源频率是否为工频，当Park矢量轨迹长轴指向150度时，从电源端看匝间短路相为A相；当Park矢量轨迹长轴指向90度时，匝间短路故障相为B相；当Park矢量轨迹长轴指向30度时，匝间短路故障相为C相。 绕组

18、匝间短路相别判据 （二） 当电机三相绕组为星形接法，Park矢量轨迹长轴呈水平方向时，从电源端看匝间短路相为A相；当Park矢量轨迹长轴指向120度时，匝间短路故障相为B相；当Park矢量轨迹长轴指向60度时，匝间短路故障相为C相。7.4 气隙偏心故障的诊断 气隙偏心故障特征频率f1mfr （fr（1s）f1/p ） 三种诊断方法的对比： 单相电流频谱分析 扩展Park矢量方法 基于Park变换的融合电流频谱分析 单相电流频谱分析方法，简单易行并便于和基频进行比较同一气隙偏心故障电流 050100150200250300350-40-200204060I/dBf/Hz0501001502002

19、50300350-60-40-200204060I2/dBf/Hz050100150200250300350-40-200204060I/dBf/Hz a扩展Park矢量方法 b基于Park变换的方法 c单相电流分析方法 7.5 变频调速电机的故障诊断 绕组绝缘故障尤其是匝间短路等故障发生的几率大大增加 转子断条故障发生几率会明显下降 变频调速电机故障诊断的侧重点应该是绕组绝缘故障变频电流与工频电流波形比较 050100150200-0.03-0.02-0.0100.010.020.03t/msI/A 050100150200-0.4-0.3-0.2-0.100.10.20.30.4t/msI

20、/A 变频器输出电流变频器输出电流工频电源输出电流工频电源输出电流变频调速电机故障诊断方法评述v一定条件下，扩展Park矢量方法仍然能有效诊断转子断条故障vPark矢量方法依然是诊断变频调速电机绕组匝间短路故障的有效方法v常规方法受到严重影响，甚至无法进行有效诊断同一转子断条故障电机，由变频电源、工频电源供电时的频谱01020304050-20-1001020304050I2/dBf/Hz01020304050-20-100102030405060I/dBf/Hz图6-16 单相电流频谱（35Hz） 图6-17 扩展Park矢量频谱（35Hz） 01020304050-60-40-200204

21、0I2/dBf/Hz01020304050-40-200204060I2/dBf/Hz图图6-18扩展扩展Park矢量频谱（工频，空载）图矢量频谱（工频，空载）图6-19扩展扩展Park矢量频谱（工频，满载）矢量频谱（工频，满载） 基于小波变换的转子断条故障的诊断 v基于Park矢量的融合电流小波分析方法对于电机负载频繁变化的情况，也具有一定的适应性 v基于启动电流的小波诊断方法 基于启动电流的小波诊断方法 利用小波变换的时频局部化特性，捕捉电机启动时，断条特征频率（12ks）f1由0逐渐逼近工频f1的特性 根据不同尺度对应的低频带分解信号的变化，就可以判断出有没有频率逐渐变化的信号成分，进而

22、诊断出有无断条故障 采用连续复值Morlet小波bcfttfibeeft/221)(连续复值连续复值Morlet小波函数表达式：小波函数表达式：基于启动电流的小波诊断方法基于启动电流的小波诊断方法断条故障电机启动电流信号的小波分解结果 02004006008001000-1.5-1-0.500.51ti0200400600800100000.511.522.5bC(a,b)启动电流信号启动电流信号 a=40(频窗中心频窗中心50Hz) 0200400600800100000.10.20.30.40.50.60.7bC(a,b)0200400600800100000.050.10.150.20.

23、250.30.35bC(a,b) a=80(频窗中心频窗中心25Hz) a=200(频窗中心频窗中心10Hz) 良好电机启动电流信号的小波分解结果010002000300040005000-0.8-0.6-0.4-0.200.20.40.60.8ti01000200030004000500000.511.522.53bC(a,b) 启动电流信号启动电流信号 a=100(频窗中心频窗中心50Hz)01000200030004000500000.050.10.150.20.250.30.350.4bC(a,b)01000200030004000500000.050.10.150.20.250.30

24、.35bC(a,b)a=200(频窗中心频窗中心25Hz) a=250(频窗中心频窗中心20Hz)基于启动电流的小波诊断方法的优点 具有很高的诊断灵敏性 从根本上消除了负载变化带来的不利影响7.6 电机故障的信息融合诊断模型 C相电流B相电流A相电流三相电压绕组温度Park矢量融合方法Park矢量融合电流Park矢量轨迹图象电压Park矢量轨迹Park矢量电压融合方法单相电流FFT数数据据层层融融合合特特征征层层融融合合人工神经网络（RBF）电压校正Park矢量轨迹人工神经网络（BP）转子断条故障绕组匝间短路故障气隙偏心故障一级诊断一级诊断提取特征矢量极坐标系象素法提取特征矢量温度曲线绕组温度

25、异常小波分析启动电流匝间短路相别轨迹倾斜角二级诊断二级诊断证据理论证据理论决决策策层层融融合合实例一 表6.2中第2组故障模式矢量的BP网络输出为0.9006，0.3984，0，0.0009，初步诊断结果为电机正常。增加融合电流信号的DWT分解结果，和启动电流信号的连续复值小波分解结果，进行决策层融合诊断 分别用F1、F2、F3、F4表示电机的4种状态：正常、转子断条、气隙偏心和绕组匝间短路，构成故障诊断识别框架UF1，F2，F3，F4，用F表示框架U的不确定性。 BP网络输出的判据结果用基本概率赋值函数m1表示，融合电流信号的DWT分解和启动电流信号的连续复值小波分解的判决结果，分别用函数m

26、2和m3表示。 对BP网络的输出进行归一化，结果为0.6928，0.3065，0，0.0007。得到电机正常F1和转子断条F2的基本概率赋值（因BP网络不具有识别F3、F4的能力，对应F3、F4的输出归入F） 设Bel1和Bel2是同一识别框架U上的两个信任函数，m1、m2分别是其对应的基本概率赋值，焦元分别为A1，A2，Ak和B1，B2，Br，又设1)()(211jiBAjiBmAmKCUCKBmAmCCmCBAjiji1211)()(0)(实例一实例一 F1F2Fm10.69280.30650.0007m20.80000.10000.1000m30.90000.05000.0500 m2(

27、F1)m2(F2)m2(F)m1(F1)F1(0.5542)(0.0693F1(0.0693)m1(F2)(0.2452)F2(0.0307)F2(0.0307)m1(F)F1(0.0006)F2(0.0001)F(0.0001)K1=(0.2452+0.0693)=0.3145 m12（F1）（0.55420.00060.0693）/(1-K1)=0.9103m12（F2）（0.03070.00010.0307）/(1-K1)=0.0896m12（F）0.0001/(1-K1)=0.0001m123（F1）（0.81930.00010.0455）/(1-K2)=0.9897m123（F2）0

28、.0103，m123（F）0 m2(F1)m2(F2)m2(F)m1(F1)F1(0.5542)(0.0693F1(0.0693)m1(F2)(0.2452)F2(0.0307)F2(0.0307)m1(F)F1(0.0006)F2(0.0001)F(0.0001)实例二 表6.4中第8组故障诊断实例中，因样本取自启动电流信号，样本本身已经失去了匝间短路的主要特征，从而导致RBF网络的输出结果（m1）与实际不符。如果有电压Park矢量轨迹的比较结果（m2）和绕组温度异常的结果（m3），通过证据组合，可纠正RBF网络的错误诊断结果。 m1F1，F4，F0.9000，0.0800，0.0200m2

29、F1，F4，F0，0.99，0.01m3F1，F4，F0.05，0.9，0.05m123F1，F4，F0.0095，0.9904，0.0001v纠正了先前的错误诊断结论！ 实例二实例二 实例三 表6.2中第6组故障模式矢量的BP网络输出为0.1121，0.9336，0，0，初步诊断结果为转子断条。增加融合电流信号的DWT分解结果，和启动电流信号的连续复值小波分解结果，进行决策层融合诊断。 m1F1，F2，F0.1072，0.8928，0m2F1，F2，F0.1000，0.8000，0.1000m3F1，F2，F0.0010，0.9000，0.0990根据证据理论，进行证据组合，最终结果为： m123F1，F2，F0.0001，0.9878，0.0121 实例四 选取转子有2根断条的故障电机空载运行时的三相电流，提取特征矢量，送入训练好的BP网络。网络输出为0.9999，0.0015，0，0.