基于分数阶傅里叶变换的无刷直流电机逆变器故障诊断技术研究.docx

? ? 基于分数阶傅里叶变换的无刷直流电机逆变器故障诊断技术研究 ? ? 吴鑫，崔江，陈则王 (南京航空航天大学 自动化学院，江苏 南京 210016) 0 引言 无刷直流电机(BLDCM)因其结构简单、运行效率高、调速性能好，已经越来越广泛地应用在各个领域中。无刷直流电机驱动系统中逆变器的功率半导体管器件是最易发生故障的薄弱环节，其故障率约占整个系统故障的82.5%[1]。因此，准确定位逆变器故障，据此进行容错控制，是提高无刷直流电机驱动系统运行安全性和可靠性的根本。 逆变器功率管故障可分为短路和开路。功率管短路故障已有成熟的解决方案，即通过硬件电路检测其D-S压降，可准确定位故障管[2]。而功率管开路时，电机往往还能继续运行，不易发现，从而导致二次故障，引发更大的事故。文献[3]对现有的逆变器功率管开路故障诊断方法进行了归类和总结，可分为基于电流量的诊断方法和基于电压量的诊断方法。文献[4]提出了一种基于傅里叶变换和神经网络的逆变器故障检测与诊断方法，利用了加窗傅里叶变换提取故障特征。文献[5]对电机三相电流进行小波变换提取故障特征，并通过BP神经网络识别故障类型。 本文采用分数阶傅里叶变换(FRFT)对输出相电流进行故障特征提取，并与时域、傅里叶变换及小波变换提取的特征进行比较、分析，寻求最优的故障特征,最后运用支持向量机进行故障诊断。 1 分数阶傅里叶变换定义及特点 传统的傅里叶变换是分析和处理平稳信号的一种标准和有力的工具，而对于分析和处理时变的非平稳信号则显得乏力，这是由于傅里叶变换采用的全局性的基函数所决定的。FRFT方法已经在信号检测、滤波器、图像处理等领域得到广泛的重视和应用[6-9]。FRFT作为傅里叶变换的一种广义形式，可以解释为信号在时频平面内坐标轴绕原点逆时针旋转任意角度后构成的分数阶傅里叶域上的表示方法。如图1所示，其中α=pπ/2,p表示FRFT的分数阶次。 图1 分数阶傅里叶域示意图 信号的p阶FRFT可表示为： (1) 其中kα(t,u)为分数阶傅里叶变换的核函数： kα(t,u)= (2) 常用的核函数有线性核函数、多项式核函数、Sigmoid核函数和径向基核函数。p阶分数阶傅里叶域是在(t,w)平面上按逆时针方向旋转的α角度所产生的坐标空间。由于分数阶域的对称性和周期性,p取值一般在0～1之内,即α取0～π/2。当p=0时，零旋转对应于信号本身；当p=1时，信号的FRFT等于经典的傅里叶变换。当p从0变到1时，FRFT能给出信号逐渐从时域变化到频域的所有特性。 分数阶傅里叶变换同小波变换一样，具有良好的时频局部化特性。可由公式推导得，将式(2)带入式(1)得： (3) (4) 分数阶傅里叶变化的p值直接影响信号在新的分数阶空间的特征。文献[10]通过穷举法和最后的诊断结果确定最优p值。文献[11]提出一种基于二阶矩极值点的方法确定p值，具有最佳谱聚集性。文献[12]提出了通过遗传算法和类内类间距离确定最优p值,获得了较好的诊断结果。 2 实验研究 2.1 仿真实验 采用MATLAB搭建无刷直流电机逆变器模型，其原理图如图2所示。 图2 无刷直流电机逆变器原理图 T1～T6为功率管，D1～D6为续流二极管，G1～G6为功率管的驱动信号，由BLDCM的霍尔位置信号判断决定。 仿真参数设置为：电机定子相绕组电阻R=0.6Ω，定子相绕组自感与互感之差L-M=0.93mH，转动惯量J=0.00009kg·m2，额定转速n=3000r/min，极对数P=4，48V直流电源供电。无刷直流电机逆变器工作在两两导通方式，每隔60°电角度换相一次，每个功率管导通120°。采样频率为5MHz，逆变器各个功率管分别设置在0.03s时发生故障，本文只研究单个功率管开路情况。设定电机定子相绕组电阻容差10%，定子相绕组电感容差5%，采集A相电流值，对单个功率管开路及正常工作共7种故障模式各进行50次蒙特卡洛分析，共采样350组数组。可得到正常工作及单个功率管开路时的A相电流波形如图3所示。 图3 正常工作及单个功率管开路时的A相电流波形 由图3可以看出，T1开路时，A相不再与电源正极相接，A相电流不再为正值；T2开路时，A相电流不再为负值；其它管子开路时，由于反电动势的影响，A相电流都有明显变化。 2.2 故障特征提取 对于时域分析，通过采集电流波形信号一个周期的点，再等间隔采样32个点作为故障特征。对时域采集的特征进行FFT，取其幅值，作为频域特征。对于小波分析，根据电流信号波形和小波函数的相似度，并且考虑小波的消失矩、正则性、支撑长度等参数，经比较分析后，采用了db8小波对电流信号进行6层分解，将获得的频带能量值作为故障特征[13]。对于分数阶傅里叶域，以类内类间距离为判据，对信号进行分数阶傅里叶变换，寻