基于matlab的电机故障诊断
基于 matlab 的电机故障诊断 班级:电气 09-6 班 组员: fudongshan……….……….……….一.引言三相鼠笼式异步电动机因其结构简单、使用方便、成本较低等优点在工业生产中得到 广泛的应用。但是由于工作环境恶劣,或者频繁起动等原因,异步电机有很大一部分故障属 于转子故障,而且主要是断条和端环开裂。如何防止电机严重故障的发生,改变现有的电机 维修检修制度,降低维护成本,已经成为工程上具有安全和经济两方面重大意义的事情。 三相异步鼠笼电动机在三相交流电源的驱动下,产生旋转的磁场,并产生旋转地磁势,磁 路的磁阻均匀,电机在额定负载下的三相电流为正弦波,且互差 120 度的相位角。但当电 机出现故障,会导致磁路的磁阻出现变化,会导致电机在额定负载下的三相电流含有一定 的谐波。通过对电机三相电流谐波的分析,可对电机故障进行诊断。以下是对一根断条和 三根断条的实验数据分析。实验数据为采样频率为 5000hz 的电流数值。 如图 1,部分实验电流数据的分布图形图 1 传统的基于傅立叶变换的FFT 频谱分析对平稳随机信号分析和处理很有效, 然而它只 对信号中的正弦成分进行统计, 实际的振动诊断信号中可能包含早期的微弱信号与大量的 非平稳信号, 比如摩擦、旋转失速、机械松动、电磁故障等。利用基于傅立叶变换的频谱分析显得无能为力。后来采用加窗FFT 也就是STFT (短时FFT )。它的基本思想是把信号分成 很多段, 每段近似为平稳信号。但是STFT 是单分辨率的分析方法, 适用于分析具有固定不 变带宽的突变信号, 无法对非平稳信号进行完全的分析。近年来兴起的小波分析技术具有 良好的 时频局部化特性, 不仅可以分析平稳的随机信号,还可以分析非平稳的随机信号。 在实际应用中通常有以下几种方法进行检测: (1)起动电流时变频谱分析法 通过采集笼型异步电机在启动过程中的定子电流信号,对其进行频谱分析,获取时变 频谱,据此判断转子有无断条。具体原理是:在启动过程中,频谱图中(1-2s)f1频率分 量可以远离f1频率分量,降低对频率分辨力的要求;而且在启动过程中,(1-2s)f1频率 分量与f1频率分量幅值之比远大于稳态运行时的幅值之比,因而使(1-2s)f1频率分量易 于检测。该方法具有很高的诊断灵敏性,从根本上消除了负载变化带来的不利影响。不过 不足之处是必须等电机启动才能用。 (2)Park 矢量法 基于定子电流 Park 矢量轨迹的转子断条故障检测方法其基本思路是:在理想情况下, 笼型异步电机定子三相电流对称并按正弦规律变化,且仅含 f1 频率分量,其 Park 矢量复平面轨迹是一个以坐标原点为圆心的圆形;而在转子断条故障情况下,定子三相电流中将 出现(1-2s)f1 与(1+2s)f1 频率分量,Park 矢量复平面轨迹成为圆环形,据此即可判断转子断 条与否。Park 矢量方法与傅立叶频谱分析相结合,即成为扩展 Park 矢量方法,应用该方 法可以进一步判断故障严重程度。 (3) clerk 矢量法Clerk 基于定子电流 Park 矢量轨迹 在旋转坐标系下交流电机三相对称的静止绕组 A 、B 、C ,通以三相平衡的正弦电流时, 所产生的合成磁动势是旋转磁动势 F,它在空间呈正弦分布,以同步转速 1 (即电流 的角频率)顺着 电流的相序旋转,两相静止坐标系到两相旋转坐标系 M、T 变换称作两 相—两相旋转变换,简称 2s/2r 变换,其中 s 表示静止,r 表示旋转。另一种很有用的坐 标系是两相同步旋转坐标系,其坐标轴仍用 d,q 表示,只是坐标轴的旋转速度 dqs 等 于定子频率的同步角转速 1 。而转子的转速为 ,因此 dq 轴相对于转子的角转速 dqr = 1 - = s ,即转差。代入式(6-105) ,即得同步旋转坐标系上的电压方程。在进 行极坐标绘图,在理论上理想电机图形应为一点。(3)小波分析方法 此方法应用于转子断条故障检测的基本思路是:首先通过小波分解计算定子电流小波 系数,之后对小波系数作频谱分析,从而突出了转子断条故障特征分量---(1-2s) f1、(1+2s)f1等边频分量。 二.基于小波分析的转子断条故障诊断 2.1连续小波基函数 小波(wavelet) ,即小区域的波,是一种特殊的长度有限、均值为零的波形。其函 数的确切定义为:设 为一平方可积函数,若其傅里叶变换满足条件: φ( 𝑡 ) (1) 则称 为一个基本小波或小波母波。将小波母函数 进行平移和伸缩就 φ( 𝑡 ) 可以得到函数 式中a为伸缩因子(也称尺度因子) , 为平移因子,我们称 为依赖于参 τ 数a, 小波积函数。 τ 为 2.2连续小波变换与离散小波变换的定义 将任意平方可积函数f(t)在小波基下展开,称这种展开为函数f(t)的连续小波变换,其表达式为: 由上定义知,小波变换和傅里叶变换都属于积分变换, 为小波变换系 数。但二者不同的是,小波基具有尺度和平移两个参数,所以函数经小波变换,就意味 着时间函数投影到二维的时间—尺度相平面上,这样有利于提取信号的某些特征。相应 的离散小波变换定义为: 如果去a0=2, 0=1离散化,且如果 满足 τ 我们称 为二进正交小波。上理想电机图形应为一点。 b c f t t f i b e e f t / 2 2 1 ) ( 程序: 傅里叶变换: fs = 5000; N=20056; temdat= xlsread( d:\22.xls );A = temdat(1:N,1); A=A-mean(A); AFFT = fft(A,N); magx = abs(AFFT); k=[0:N/2-1]/N*fs; subplot(2,2,1) plot(k,magx(1:N/2)) axis([0,100,0,1000]); grid on; magx=magx ; n=N/5000; y=max(magx(1,1:n)) for i=1:49m=max(magx(1,n*i:n*(i+1)));y=[y,m]; end i=0:49; subplot(2,2,2) bar(i,y) Park程序: x = xlsread( d:\11.xls ) t=(0:0.0002:4.0110); T=t ; Ia=x(1:20056,1); Ib=x(1:20056,2); Ic=x(1:20056,3); ia=sqrt(2/3).*(Ia-Ib./2-Ic./2); ib=sqrt(2/3).*(sqrt(3)/2.*Ib-sqrt(3)/2.*Ic) is=sqrt(ia.*ia+ib.*ib); O=ib./ia; F=atan2(ib,ia); G=F; polar(G,is) Clerk程序 x = xlsread( d:\11.xls ) t=(0:0.0002:4.0110); T=t ; Ia=x(1:20056,1); Ib=x(1:20056,2); Ic=x(1:20056,3); ia=sqrt(2/3).*(Ia-Ib./2-Ic./2)