用MATLAB实现转子的故障诊断,单盘转子故障诊断报告(附matlab程序).doc

单盘转子故障诊断报告(附matlab程序).doc

单盘转子故障诊断报告1数据TEST1TXT，数据长度1556486，2组振动数据第一、二列为测点1涡流传感器X、Y振动数据；第三、四列为测点2涡流传感器X、Y振动数据；第五列为转速脉冲信号(每转1个)；第六列为转速值(采集仪输出转速，可能不准确)；采样频率2000256HZ实验转速3300RPM2故障分类转子试验台常见故障类型有不平衡、不对中、松动、裂纹、碰摩等。常见故障类型故障特征如下1、转子质量不平衡主要特征1一倍频率振动幅值大。2出现较小的高次谐波，整个频谱呈所谓的“枞树形”。2、转子不对中故障主要特征1从振动的时域波形上可以看出旋转基本频率的高次成分。2从振动信号的频谱图上可以发现工频的高次分量，如2FR和3FR振动，尤其是2FR振动非常明显。3当不对中比较轻微时，轴心轨迹呈椭圆形；当不对中故障达到中等程度时，轴心轨迹呈香蕉形；当不对中故障较严重时，轴心轨迹呈外“８”字形。3、转子部件脱落主要特征有1转子部件脱落后，转子的振动频谱与质量不平衡时的振动频谱类似。2转子部件脱落的前后，振动的幅值和相位突然发生变化。3部件脱落一段时间后，振动的幅值和相位趋于稳定。4轴心运动轨迹与质量不平衡时的轴心运动轨迹类似。4、动静碰磨故障特征有1振动的时域波形特征当转子发生碰磨故障时，振动的时域波形发生畸变，出现削波现象。另外，在振动信号中有奇异信号。2振动的频谱特征由动、静部分碰磨而产生的振动，具有丰富的频谱特征。振动有时还会随着时间发生缓慢的变化。3轴心运动轨迹特征①若发生的是整周碰磨故障，则轴心运动轨迹为圆形或椭圆形，且轴心轨迹比较紊乱。②若发生的是单点局部碰磨故障，则轴心运动轨迹呈内“８”字形。③若发生的是多点局部碰磨故障，则轴心运动轨迹呈花瓣形。4当转轴与静子发生碰磨时，会使转子产生振幅时大时小、振动相位也时大时小的旋转振动。5、转子裂纹故障主要特征1由于裂纹的存在改变了转子的刚度，从而使转子的各阶临界转速较正常值要小，裂纹越严重，各阶临界转速减小得越多。2由于裂纹造成刚度变化且不对称，从而使转子的共振转速扩展为一个区域。3裂纹转子在做强迫响应时，一次分量的分散度较无裂纹时大。4在恒定转速下，各阶谐波幅值及其相位不稳定，尤以二倍频最为突出。辅助特征5裂纹转子引起的刚度不对称，使转子动平衡发生困难，往往多次试加平衡质量也达不到所要求的平衡精度。3分析步骤31转速分析对第6列数据进行分析可得转速的时域和频域图，如图1所示图1转速时域和频域图由图1可看出单盘转子转速较为稳定，可认为其是恒定转速。32测点1处X和Y方向振动信号分析绘制测点1处X和Y方向振动信号时域图如图2图2测点1处X和Y方向振动信号图由图2看出X和Y方向振动幅值在9S左右明显变化，Y方向的振动位移明显比X方向大，由此可判断转子在此时发生了故障。继续对X方向时域图进行局部放大如图3所示图3X方向时域局部放大图由图可看出X方向时域波形产生了削波现象，由前文所列故障现象可初步判断该故障为碰摩。33轴心轨迹分析(故障前和故障后)前文得出9S左右发生故障，因此首先对9S前转子正常转动的轴心轨迹进行分析，轨迹图如图4，大致为一椭圆形。图4,故障前轴心轨迹此后继续对故障发生后的轴心轨迹进行分析，如图5，轨迹呈内八字形，进一步确定产生的故障为碰磨，且为单点局部碰磨。结论通过对单盘转子转速、X和Y方向振动信号的时域及频域的分析，得知其时域波形出现削波现象；此后对转子故障前后轴心轨迹进行分析，其轴心轨迹呈内8字形态，最终判断出故障类型为恒定转速下的单点碰磨。附MATLAB程序程序1(对单盘转子的转速分析)FID WKTXT , R AFSCANFFID, F ,6,INFFCLOSEFIDN1155648WK1DATAA1,1N1取第6列元素WK1DATAWK1DATAMEANWK1DATA去掉直流分量WK1DATAWK1DATA11N1每个十个点去一个数NN1/10WK1WK1DATA1N取得数新组成的矩阵N0N1FS5120TN/FSSUBPLOT2,1,1,PLOTT,WK1TITLE 转速时域图 XLABEL 时间 YLABEL 幅值 YFFTWK1DATA,NFFT变换MAGABSY2/N振幅FNFS/NSUBPLOT2,1,2,PLOTF1N/2,MAG1N/2TITLE 转速频域图 XLABEL 频率/HZ YLABEL 幅值 GRIDON程序2(对测点1X和Y方向的分析)FID WKTXT , R AFSCANFFID, F ,6,INFFCLOSEFIDN155648WK1DATAA1,1N取第1列元素WK1DATAWK1DATAMEANWK1DATA去掉直流分量WK2DATAA2,1N取第2列元素WK2DATAWK2DATAMEANWK2DATA去掉直流分量WK1WK1DATA1N取得数新组成的矩阵WK2WK2DATA1NN0N1FS5120TN/FSSUBPLOT2,2,1,PLOTT,WK1TITLE X方向时域图 XLABEL T YLABEL 幅值 SUBPLOT2,2,2,PLOTT,WK2TITLE Y方向时域图 XLABEL T YLABEL 幅值 YFFTWK1,NMAGABSY2/NFNFS/NSUBPLOT2,2,3,PLOTF1N/2,MAG1N/2TITLE X方向频域图 XLABEL 频率/HZ YLABEL 幅值 YFFTWK2,NMAGABSY2/NFNFS/NSUBPLOT2,2,4,PLOTF1N/2,MAG1N/2TITLE Y方向频域图 XLABEL 频率/HZ YLABEL 幅值 程序3(故障前轴心轨迹)FID WKTXT , R AFSCANFFID, F ,6,INFFCLOSEFIDN150000WK1DATAA1,1N取第一列的元素WK1DATAWK1DATAMEANWK1DATA去掉直流分量WK1DATAWK1DATA11N/4取故障发生前的数据WK2DATAA2,1N取第2列的元素WK2DATAWK2DATAMEANWK2DATA去掉直流分量WK2DATAWK2DATA11N/4NN/4WK1WK1DATA1N取得数新组成的矩阵WK2WK2DATA1N取得数新组成的矩阵N0N1FS5120TN/FSPLOTWK1,WK2TITLE 故障前轴心轨迹 程序4(故障后轴心轨迹)FID WKTXT , R AFSCANFFID, F ,6,INFFCLOSEFIDN150000WK1DATAA1,1N取第一列的元素WK1DATAWK1DATAMEANWK1DATA去掉直流分量WK1DATAWK1DATA1000001N取故障发生前的数据WK2DATAA2,1N取第2列的元素WK2DATAWK2DATAMEANWK2DATA去掉直流分量WK2DATAWK2DATA100001NN50000WK1WK1DATA1N取得数新组成的矩阵WK2WK2DATA1N取得数新组成的矩阵N0N1FS5120TN/FSPLOTWK1,WK2TITLE 故障轴心轨迹