基于随机森林的电机异音故障诊断方法.pdf

Jun.2023Journal of Yancheng Institute of Technology(Natural Science Edition)2023年0 6 月Vol.36 No.2盐城工学院学自然科学版）第36 卷第2 期doi:10.16018/32-1650/n.202302007基于随机森林的电机异音故障诊断方法钟张豪1-2,丑永新,侯千红1.21.盐城工学院机械工程学院，江苏盐城224051;2.常熟理工学院电气与自动化工程学院，江苏苏州215000摘要：针对电机异音故障检测技术存在准确率低、模型复杂等问题，提出一种基于随机森林的电机异音故障诊断方法。通过自行研制的汽车智能座椅靠背电机振动测试平台，分析电机故障产生过程及异音的特征，并从时域中提取11个特征表征异音信号的变化；通过主成分分析法对所提取的特征进行降维，在训练基于随机森林和概率神经网络的电机故障智能识别方法基础上，通过自行研制的汽车智能座椅靠背电机振动测试平台采集数据，得到随机森林的平均识别准确率为9 5.11%2.17%，概率神经网络的平均识别准确率为9 3.9 0%2.16%。关键词：电机；振动信号；异音识别；随机森林；主成分分析；特征提取中图分类号：TP181文献标志码：A文章编号：16 7 1-532 2(2 0 2 3)0 2-0 0 37-0 7电机是现代生产和生活中十分重要的机电设备，其工作状态好坏直接影响整个系统的稳定性。由于电机转子不平衡、轴承磨损等产生的异常震动会在电机外部产生“滴滴答答”等的异音，研究表明8 0%以上的电机故障都会产生异音1-3,因此可将异音判别作为电机故障诊断的重要依据。目前，我国电机出厂故障检测主要依靠人工在静音房听音，这种方式耗时耗力且容易出现误判。因此，实现电机异音故障的智能检测十分重要。到目前为止，国内外已经有了电机异音诊断方法的研究，这些方法主要从信号的频谱分析、时频分析、机器学习建模3个方面展开。在频谱分析方面，通过傅里叶变换得到电机音频信号的频谱图，再根据各种故障下的故障特征频率信息对异音信号进行分辨4。由于电机振动信号属于非平稳信号，使用一般的频谱分析法如傅里叶变换等的效果欠佳，并且傅里叶变换只能进行整段信号的频谱分析，不能对异音发生的时刻进行定位。因此，时频分析方法如短时傅里叶变换、小波变换和经验模态分解等被引人电机异音的检测中。短时傅里叶变换虽然拥有一定的时间分辨率，但窗口一旦选定，则不能再变化，不具有自动调节的功能5；小波变换作为一种新的时频分析方法，继承和发展了短时傅里叶变换局部化的思想，同时克服了窗口大小不随频率变化等的缺点，很适合运用在非平稳信号中6-7 ；经验模态分解将振动信号进行自适应分解后，提取方差贡献率高的分量重构新的信号，可以有效去除噪声以及不包含故障信息的频带，在电机异音故障诊断方面效果很好8 ，但是传统的经验模态分解容易产生模态混叠和虚假分量，影响最终结果。随着机器学习分析技术的发展，电机异音故障诊断问题迎来新的解决方法。反向传播（backpropagation,BP)神经网络算法通过自身训练学习规则，对电机异音识别的准确率可以达到9 0%9 ,但是BP神经网络也拥有学习速度慢，容易陷人局部极小值等缺点。支持向量机是一类按监督学习方法对数据进行二元分类的广义线性分类器，使用支持向量机对电机异音信号进行识别得到了电磁异音9 1.7 6%、轴承异音10 0%、摩擦异音99.24%以及松动异音9 8.7 5%的识别结果【10 ，但收稿日期：2 0 2 2-11-11基金项目：国家自然科学基金青年基金项目（6 19 0 10 6 2）。作者简介：钟张豪（19 9 9 一），男，浙江温州人，硕士生，主要研究方向为振动信号处理、电机异音故障诊断等。38第36 卷盐城工学院学报（自然科学版）是支持向量机对线性分类问题具有较高识别率，对非线性分类问题的识别率较低。卷积神经网络（convolutional neural networks，CNN）是一类包含卷积计算、具有深度结构的前馈神经网络，是深度学习的代表算法之一，其代表特点是局部连接、权值共享以及时间或空间的下采样。将卷积神经网络与迁移学习结合后对电机轴承故障进行诊断，准确率达到9 9.8%1。卷积神经网络优点是不需要进行特征提取，缺点是需要大量样本进行训练，培训成本高、花费时间长，而实际生产中产出的次品数量很少，难以满足卷积神经网络需要的大量样本的要求。与BP、CNN等机器学习方法相比，随机森林（r a n d o m f o r e s t，R F)是基于规则的学习，具有较高的分类精度，能够清楚地了解训练模型的结构和特征，并且RF是基于小样本学习的一种机器学习方法，因此本文提出一种基于随机森林模型的电机异音故障诊断方法。该方法首先搭建一套汽车智能座椅靠背电机振动测试平台，用于采集实验数据；然后提取电机振动信号特征，并通过主成分分析（principal component analysis，PCA）法对故障特征进行降维；最后分别训练RF分类模型和概率神经网络（probabilistic neuralnet-work，PNN），对其识别性能进行评估。1实验数据实验数据由搭建的汽车智能座椅靠背电机振动测试平台采集，内部结构如图1所示，包括电柜与PC机、显示器、打印机、测试平台、按钮、阀岛、传感器分线盒和一个气泵组成。其中按钮用于控制测试的启动、停止、复位或者调整手自动模式；阀岛用于将气泵的气传输至各气缸；传感器分线盒是数字信号传感器的集中连接器，用于将多个数字信号传感器的信号集合并统一与PC进行信息交互；电柜用于集中的电气控制；打印机用于测试结束后将测试结果打印出来；测试平台用于将所测汽车靠背固定、测试并采集数据。数据采集时，将汽车靠背放在测试位置并将电机插头插上，按下启动按钮后左右两边的定位梢气缸伸出将汽车靠背固定住；左右两边的传感器夹紧气缸夹紧振动传感器，伸缩气缸伸出并将头部装配有磁铁的传感器紧贴在汽车靠背的转轴处；再次按下启动按钮，上位机程序控制电机的插头输出+13.5V或者-13.5V的电压，以控制显示器电柜与PC机测试平台打印机按钮阀岛传感器分线盒图1汽车靠背电机振动测试平台结构Fig.1Structure of vibration testing platform forautomobilebackrestmotor靠背电机的正反转，使靠背实现图2 所示的运动周期。图2 中当汽车靠背处在上料位置与下料位置时，靠背电机的角度是相同的，此时通过贴在汽车靠背转轴处的两个振动传感器采集振动信号；振动信号经过放大器放大后，被USB-6212的数据采集卡接收，然后通过上位机LabVIEW程序中的DAQMX将信号采集并以TDMS形式保存。信号采集与传输流程如图3所示。前转电机后转电机上料位置前极限后极限前转后转电机下料位置前极限图2运动周期流程图Fig.2Flowchart ofmotion period传感器接信号USB-6212收电机振放大器数据采集卡动信号上位机保存数据LabVIEW数据采集图3信号采集与传输流程图Fig.3Flow chart of signal acquisition andtransmission汽车靠背电机振动测试平台所使用的传感器为丹麦B&K公司生产的4534-B-002型传感39钟张豪，等：基于随机森林的电机异音故障诊断方法第2 期器，如图4所示。该振动传感器使用密封钛外壳和绝缘底座，采用10-32 UNF螺纹安装，不仅重量轻，还非常坚固，且拥有频率范围广、安装范围广、噪音低、对环境敏感度低的特性，适合多种环境条件下使用。4534-B-002型传感器参数如表1所示。图44535-B-002型振动传感器Fig.44535-B-002vibration sensor表14534-B-002型传感器主要参数Table14534-B-002 main parameters of the sensor参数数值电压灵敏度/(mV/ms-2)50(110%)振幅频率范围/Hz0.312.800相位频率范围/Hz21500安装谐振频率/kHz38最大横向灵敏度/%5测量范围/g14数据采集时，将采样频率设为2 0 kHz，每30 35s采集1组数据，共采集9 2 组数据，其中2 4组数据为不合格数据。某组不合格数据如图5所示。64一2K异音段0-2-4一异音段-60102030时间/s图5不合格电机振动信号图Fig.5Vibration signal diagram of unqualified motor2电机异音故障诊断电机异音故障诊断流程如图6 所示，包括数据采集、特征提取、PCA降维、划分训练集和测试集、建立RF和PNN异音识别分类器模型、评估分类性能几个部分。靠背电机特征提取PCA降维振动信号结果评估机器学习图6电机异音故障诊断方法流程图Fig.6Flow chart of motor abnormal sound faultdiagnosis method2.1特征提取特征提取时，采用窗口法，将整个信号以窗口为单位（窗宽设定为采样频率的1/5，即w=f/5），对每个窗口里的小信号分别提取峰峰值、平均值、均方值、标准差、有效值、峰值因子、脉冲因子、波形因子、裕度因子、偏度因子和峭度因子等11个特征指标，从而形成一个特征向量作为分类器的输人。峰峰值描述了电压信号变化范围的大小，是电压信号最大值和最小值之间的差值，计算公式如式(1)所示。Xpp=xmax-Xmin(1)式中:XpvXmaxvXmin分别为电压信号的峰峰值、最大值和最小值，V。平均值是一段时间内电机振动电压信号的平均值，计算公式如式（2）所示。1N=;(2)N=1式中：x为离散信号点的电压，V；为电压信号的平均值，V；N为信号的采样点个数，个。均方值是一段时间内电机振动电压信号平方的平均值，在工程上表示信号的平均功率，计算公式如式（3）所示。1x=(3)N=1式中：x。为电压信号的均方值，V?标准差是方差的算术平方根，反映电机振动信号的离散程度，计算公式如式（4)所示。(x;一x)(4)rmse式中：xmse为电压信号的标准差，V。有效值是电机振动电压信号的均方根，也称为均方根值，计算公式如式（5）所示40第3 6 卷盐城工学院学（自然科学版）(5)2Xrms式中：xm为电压信号的有效值，V。峰值因子是电机振动电压信号最大值与有效值的比值，表示电机振动电压信号峰值与整体波形相比的极端程度，计算公式如式（6)所示。Xmaxx;max(6)之rms2式中：C,为信号的峰值因子。脉冲因子是电压信号峰值与整流平均值的比值，表示信号的冲击特性，计算公式如下：Xmaxmaxx;(7)1NxN=1式中：1为信号的脉冲因子。波形因子是电压信号的均方根与绝对平均值的比值，也即电压信号脉冲因子与峰值因子的比值，表示振动信号波形与标准正弦波相比的失真程度，计算公式如式（8）所示。1N2N=1(8)Sms1N之N=1式中：S为信号的波形因子。裕度因子是电压信号峰值与方根振幅的比值，表示电机振动信号的饱满程度，计算公式如式(9)所示。XmaxL=Xmax(9)1Y式中：L,为信号的裕度因子；x为方根幅值，V。偏度因子也叫偏斜度、偏态，是电压信号三阶中心矩和标准差的三次方的比值，代表电机振动信号的峰值分布，计算公式如式（10）所示。13NSk=(10)32式中：S.为信号的偏度因子。峭度因子是电压信号四阶中心矩除以标准差的四次幂，表示信号波形的平缓程度，计算公式如式(11)所示。4NKv=(11).44XrmsN2式中：K为信号的峭度因子；为信号的四阶中心矩，V4。将所有信号以f/5为窗宽进行分割，然后进行特征提取，得到特征向量空间X，大小为12 x14238（行数表示特征，列数表示特征样本数）；再采用人工方法进行标定，其中正常段标为0，异音段标为1，从而形成标签向量C，大小为1142 38（142 38 为标签数量）。2.2PCA降维PCA降维用于降低所提取的11个电机异音信号特征的维数。在正确区分电机正常振动信号与异音信号的情况下，保留特征分类精度，通过降低网络输人的维数，减少数据穴余并提高模型的训练效率。采用MATLAB软件库中的函数c,s,l,t，e =pca(X)进行主成分分析，并实现降维。其中，c是特征向量空间X对应的协方差矩阵的所有特征向量组成的矩阵，即变换矩阵或投影矩阵，其每一列按特征值从大到小排序，分别代表特征值所对应的特征向量；s表示原数据在各主成分向量上的投影；1表示主成分方差，也就是各特征向量对应的特征值，按照从大到小排列；t.表示特征空间X中每个观察值的T平方统计量，用来衡量多变量之间的距离；e表示每个主成分的贡献率。通过搜索贡献率大于1的特征，舍弃其余特征，从而降低特征维数12 2.3基于随机森林算法的电机异音诊断随机森林是由多个决策树组成的分类器在每棵树的节点上，从所有原始特征中随机选择若干特征。特征数据在决策树中被分割产生叶节点（表示类的标签），从根节点到叶节点的路径表示一个确定的决策过程。当样本在随机森林中分类时，每棵树都会得到一个分类结果，即对类的投票，随机森林根据“少数服从多数的原则，将投票数量多的类作为最终的类别值输出。在MATLAB中，可以使用“RANDOMFORESTMATLAB工具箱实现RF分类器模型的训练与测试，其中RF分类器模型的训练可以使用函数 R.（p,u,n u r e）”（p 是特征矩阵中作为训练41钟张豪，等：基十机森杯的用机开音故障诊断方法第2 期集的部分，u是标签中作为训练集的部分，nure是决策树的数量，设为30)进行。训练完成后输出变量表示已经训练好的随机森林模型，随后使用函数 R.。（Y,mz）13-141（Y为特征矩阵中划分为测试集的部分）对训练好的模型进行测试，并对模型的准确性进行评估。2.4分类器分类性能评价指标本文选择特异性Sp、敏感度S准确性A。和卡帕系数K来评估分类器的分类性能。在汽车智能座椅靠背电机振动信号识别过程中，特异性S.表示原信号属于异音信号而被正确识别的比例，计算公式如式（12）所示。T.S.100%(12)F,+T,P式中：T表示原信号是异音信号的情况下识别为异音信号的数量，Fp表示原信号是正常信号的情况下识别为异音信号的数量。敏感度S.表示原信号属于正常信号而被正确识别的比例，计算公式如式（13）所示T,S.100%(13)T,+F.e式中：T，表示原信号是正常信号的情况下识别结果归类为正常信号的数量；F，表示原信号是异音信号的情况下识别为正常信号的数量。T,、T,、F,、F.分别代表样本检测的真阳性、真阴性、假阳性、假阴性，它们共同构成了混淆矩阵。在评价过程中，T，和T,的值越高，代表模型识别的精确度越高；F,和F.的值越高，代表模型识别的精确度越低。准确度A是指分类器预测结果正确的比例，计算公式如式（14）所示T,+T.A.=100%(14)T,+F,+F,+T.Kappa系数K用来评估模型的综合表现，取值在0,1 之间，且越接近1,分类器的分类精度越高，计算公式如式（15）所示。K=P.P.(15)1-Pe其中，P。=n=1Pem2m式中：p。是每一类正确分类的样本数量之和除以总样本数，也就是总体的分类精度；p。表示所有类别对应的实际数量和预测数量的乘积再和总样本数的平方相除；m为测试集的总样本数；r为分类的类别数，此处为2;q,为混淆矩阵中的主对角线元素；9.表示混淆矩阵中各行向量所有元素的和；9+表示混淆矩阵中各列向量所有元素的和。K,是每个分类结果所对应的Kappa系数,计算公式如式（16）所示。K=P.-P.P.P.-P.,P.(16)式中:Pu=Iu,P.,P,=9+一Ommm3结果与讨论将采集到的9 2 组数据以f/5为窗宽进行分割，得到大小为142 38 x11的特征空间；再采用人工方法进行标定，其中正常段标为0,异音段标为1,然后对提取的11个时域特征采用主成分分析法进行分析，得到各主成分的方差；根据方差得到每个主成分的贡献值，以返回主成分空间。最终得到各特征的贡献比例如表2 所示，表2 各特征贡献比例表Table2Table of contribution ratio of each feature特征编号贡献率%160.126 9220.7421311.495043.280051.844061.280670.801 280.309390.084 7100.028 4110.0078从表2 可以看出，贡献率大于1%的共有6 个特征，即这6 个特征对区分正常信号和异音信号的作用较大。将得到的6 个特征归一化后，从142 38 个样本中随机选取2 0 0 0 个样本作为测试集，剩余12238个样本作为训练集，进行RF分类器分类；再从训练集中随机抽取样本进行模型训练，然后从测试集中随机抽取样本对模型进行测试，总共进行了10 0 次重复实验；最后用特异性S。灵敏度S。、准确度A.和Kappa系数K的平均值及其标准42第3 6 卷盐城工学院学报（自然科学版）差，对基于随机森林的电机异音故障诊断模型和基于PNN的电机异音故障诊断模型的分类性能进行评价，结果如表3所示，某一组异音识别结果如图7 所示。表3RF和PNN分类结果Table3Classification results of RF and PNN参数RF分类PNN分类7,192571926742733322685665$,/%94.312.9191.003.72S./%99.860.1399.910.09A/%99.660.1899.590.16K/%95.112.1793.902.16K,/%94.103.0190.683.84K,/%96.223.0997.482.32图7 中，虚线为标记线，值为0 或1。若分类器对信号段识别为异音信号时，其值为1；若分类器对信号段识别为正常信号时，其值为0。由表3可知，基于RF模型的电机异音故障诊断分类器能够准确地对数据集中样本的正常状态和异音状态进行分类，且识别精度较高，总体分类性能也高于基于PNN的电机振动信号异音故障诊断分类器。15电机振动信号10标记505-10-150102030时间/s图7异音信号识别结果图Fig.7The recognition result of abnormal sound signal4结论结合产品生产线上的异音检测步骤，采用基于振动信号的检测方法，搭建了汽车智能座椅靠背电机振动测试平台。实验结果表明，本文提出的基于RF分类器的异音故障诊断能够在不需要大样本训练的情况下获得较好的效果，可以方便企业在实际生产中对模型进行训练；通过已训练好的模型对实际生产中的汽车靠背电机进行测试，识别精度较高，总体分类性能也高于基于PNN的电机振动信号异音故障诊断分类器。参考文献：1】付伟男。电机异音检测技术的研究D.杭州：浙江大学，2 0 18.【2 孟玉龙新能源汽车模块化永磁同步电机振动噪声分析与优化D.石家庄：河北科技大学，2 0 2 1.3】张决决，张沛全基于音频分析的电机轴承异音故障诊断研究J计量与测试技术，2 0 2 2，49(7：46-49，53.4蒋威.基于虚拟仪器的电机振动测试分析系统D.杭州：浙江大学，2 0 0 6.5】宋秀，尉宇.基于短时分数阶傅里叶变换的异步电机转子断条故障诊断J武汉科技大学学报，2 0 16，39（2）：145-149.6】张新，郑燕萍，AUGEIXA，等基于小波包分解与机器学习的汽车调光电机异音识别J.森林工程，2 0 19，35(1):59-63,112.7 AWADA E,AL-QAISI A,RADWAN E,et al.Motor fault detection using sound signature and wavelet transformJ.In-ternational Journal of Power Electronics and Drive Systems(IJPEDS),2022,13(1):247-255.8 陈宗祥，陈明星，焦民胜，等。基于改进EMD和双谱分析的电机轴承故障诊断实现J。电机与控制学报，2 0 18，2 2(5):78-83.9江毓，郑燕萍，张新，等，基于改进BP神经网络的电机异音诊断J重庆理工大学学报（自然科学），2 0 2 0，34(1):242-246,262.10】李明超.基于异音检测的电机故障诊断方法D.江门：五邑大学，2 0 14.11李俊卿，刘静。结合卷积神经网络和迁移学习的电机轴承故障诊断方法J华北电力大学学报（自然科学版），（责任编辑：李华云43钟张豪，等：基于随机森林的电机异音故障诊断方法第2 期2023,50(1):76-83,91.12 CHOU Y X,WANG P,FENG Y F.The optimal morphological model for arterial blood pressurewave related classifica-tion:comparison of two types of kernel function mixturesJ.IEEE Access,2019,8:4133-4148.13 WANG S Z,ZHANG Z F,GENG SS,et al.Research on optimization of random forest algorithm based on spark J.Computers,Materials&Continua,2022,71(2):3721-3731.14 MANTAS C J,CASTELLANO JG,MORAL-GARCIA S,et al.A comparison of random forest based algorithms:randomcredal random forest versus oblique random forestJJ.Soft Computing,2019,23(21):10739-10754.Fault Diagnosis Method of Motor Abnormal Sound Based on RandomForestZHONG Zhanghao-2,CHOU Yongxin,HOU Qianhong1，2(1.School of Mechanical Engineering,Yancheng Institute of Technology,Yancheng Jiangsu 224051,China;2.School of Electrical and Automation Engineering,Changshu Institute of Technology,Suzhou Jiangsu215000,ChinaAbstract:Aiming at the problems of low accuracy and complex model of motor abnormal sound fault detection technology,a mo-tor abnormal sound fault diagnosis method based on random forest is proposed.Based on the self-developed vibration test plat-form of automotive intelligent seat back motor,the process of motor failure and the characteristics of abnormal sound are ana-lyzed,and 11 features are extracted from the time domain to characterize the changes of abnormal sound signals.The dimensional-ity of the extracted features is reduced by the principal component analysis method.On the basis of training the motor fault intelli-gent identification method based on random forest and probabilistic neural network,data is collected through the self-developedvibration test platform of automotive intelligent seat back motor,and the average recognition accuracy of random forest is 95.11%2.17%.The average recognition accuracy of probabilistic neural network is 93.90%2.16%.Keywords:motor;vibration signal;abnormal sound recognition;random forest;principal component analysis;feature extrac-tion