一种基于峰值指标的旋转机械故障诊断方法
【专利摘要】本发明提供一种基于峰值指标的旋转机械故障诊断方法，用信号分离的思想构建峰值指标，克服传统峰值指标诊断技术的缺点，提高诊断的精确度和可靠性，通过标准振动信号将实时采集振动信号分为无故障振动信号和混合信号，用混合信号和无故障振动信号构建峰值指标，具有对故障诊断更加灵敏，更加能够通过信号的变化检测微弱故障特征信号。旋转机械设备正常运行和发生故障时，峰值指标的取值范围重叠少，不同状态下峰值指标变化明显，对偏心轴故障敏感，能够很好地将此故障与其他故障区分开来，旋转机械故障诊断的精确度和可靠性高。
【专利说明】
一种基于峰值指标的旋转机械故障诊断方法
技术领域
[0001] 本发明涉及一种旋转机械故障诊断方法，特别涉及一种基于峰值指标的旋转机械 故障诊断方法，属于故障诊断与信号处理分析技术领域。
【背景技术】
[0002] 大型旋转机械设备(如汽轮机、旋转轴承、压缩机等)是石油、化工、机械制造、航空 航天等重要工程领域的关键设备，旋转机械设备正朝着大型化、自动化、精密化的方向不断 发展，其组成和结构也变得越来越复杂，发生故障的概率也越来越大，因此对大型旋转机械 设备的安全性和可靠性的要求也越来越高。
[0003] 但是，大型旋转机械设备发生故障时，振动监测信号往往存在大量的非线性、随 机、不可遍历的信息，给故障信号的分析和旋转机械故障的诊断带来很大的困难。
[0004] 现有技术的旋转机械故障诊断主要基于振动信号进行分析，而一般采用时域分析 法，通过对机械振动信号的概率密度函数分析，推导出了幅值域中的有量纲指标和无量纲 指标，有量纲指标如均值、均方根值等;无量纲指标如峰值指标、裕度指标、脉冲指标等。
[0005] 在实际应用中，现有技术的有量纲指标对故障特征敏感，其数值会随着故障的发 展而上升，同时因工作条件，如负载、转速等的变化而变化，并极易受环境干扰的影响，表现 不够稳定。
[0006] 无量纲指标中的峰值指标对于振动信号中的扰动不敏感，对缓变故障敏感。特别 地，峰值指标对振动信号的幅值和频率的变化不敏感，即与机器的工作条件关系不大，只依 赖于概率密度函数的形状，因此峰值指标在旋转机械故障诊断中得到了广泛的应用。
[0007] 但是，现有技术采用峰值指标诊断是一种将"峰值参数"用于设备故障诊断的技术 方法，利用峰值指标来诊断旋转机械的故障主要有以下几点问题:一是没有采用信号分离 的思想构建峰值指标，没有对标准振动信号分离，无法更好地体现旋转机械的故障，对故障 的诊断不敏感;二是峰值指标要在信号有明显的脉动性时才可靠，当旋转机械的故障不断 扩展，峰值逐步达到极限值后，均方根值开始增大，峰值指标逐步减小，直至恢复到无故障 时的大小，因此峰值指标存在对故障判断失效，特别是一些微弱的故障信号反应不灵敏，易 产生误判，旋转机械故障诊断的精确度和可靠性较差;三是采用磨合期的振动信号作为标 准振动信号，没有考虑旋转机械设备在运行初期振动信号中噪声严重，不适宜作为标准振 动信号的问题，诊断故障不准确。
【发明内容】
[0008] 针对现有技术的不足，本发明提供一种基于峰值指标的旋转机械故障诊断方法， 克服传统峰值指标诊断技术的缺点，提高诊断的精确度和可靠性，用信号分离的思想构建 峰值指标，由于对标准振动信号的分离，剩余的混合信号能够更好地体现旋转机械的故障， 对故障的诊断更加地敏感，所得出的诊断结果能够准确地诊断旋转机械的故障。
[0009] 为达到以上技术效果，本发明所采用的技术方案如下：
[0010] -种基于峰值指标的旋转机械故障诊断方法，包括以下步骤：
[0011] (1)通过振动加速度传感器采集旋转机械在运行磨合期后的正常运行参数，所述 正常运行参数为按抽样频率fs采集的无故障振动信号so(t) (t = 0，1，…，T-1);
[0012] (2)对无故障振动信号so(t)进行归一化处理后得到标准振动信号s(t);
[0013] (3)对标准振动信号s(t)做快速傅里叶变换(FFT)得到标准振动频域信号S(k)(k =〇，1，…，K-1);
[0014] (4)通过振动加速度传感器采集旋转机械需要故障诊断时的运行参数，所述故障 诊断时的运行参数为按抽样频率fs实时采集旋转机械运行中的实时采集振动信号z(t)(t =〇，1，…，T-1);
[0015] (5)对实时采集振动信号z(t)做快速傅里叶变换(FFT)得到实时振动频域信号Z (k)(k = 0，l，…，K-1);
[0016] (6)对标准振动频域信号S(k)取复共辄为S(k)%将Z(k)与S(kf相乘得到Y(k)(k = 0,1，…，K-1)，再对Y(k) (k = 0，1，…，K-1)做快速傅里叶逆变换(IFFT)得到两个信号Z(k)与 S(kf 的相关函数 l(t)(t = 0，l，…，T-1);
[0017] (7)在t = 0，1，…，T-1中取相关函数I (t)的模11 (t) I，I I (t) I的最大值所对应的时 间点为实时采集振动信号和标准振动信号的延迟时间T，T = argmax| I(t) I ;
[0018] (8)求实时采集振动信号z(t)与延迟后的标准振动信号s(t-T)的相关系数c，c = E
[z(t)s(t-x)];
[0019] (9)计算混合信号y(t)，y(t)为故障特征提取信号和噪声的混合信号，y(t) = z (t)-cs(t-x);
[0020] (10)计算出峰值指标
[0021 ] -种基于峰值指标的旋转机械故障诊断方法，进一步的，根据峰值指标Cys的值，判 断旋转机械是否发生故障及故障类别，以下区间值包括边界值：
[0022]峰值指标Cys的值在2.1986到3.7719之间时，旋转机械运行正常；
[0023]峰值指标Cys的值在3.8474到6.0233之间时，考虑旋转机械发生了裂轴故障；
[0024]峰值指标Cys的值在3.8830到6.1981之间时，考虑旋转机械发生了弯轴故障；
[0025]峰值指标Cys的值在3.9094到7.5683之间时，考虑旋转机械发生了偏心轴故障； [0026]峰值指标Cys的值在4.1663到7.3144之间时，考虑旋转机械发生了裂轴+弯轴故障； [0027]峰值指标C ys的值在3.9924到7.4869之间时，考虑旋转机械发生了裂轴+偏心轴故 障；
[0028]峰值指标Cys的值在4.5938到6.2757之间时，考虑旋转机械发生了弯轴+偏心轴故 障；
[0029]峰值指标Cys的值在4.0169到6.1841之间时，考虑旋转机械发生了裂轴+弯轴+偏心 轴故障。
[0030] -种基于峰值指标的旋转机械故障诊断方法，进一步的，快速傅里叶变换(FFT)采 用有限序列离散傅里叶变换(DFT)的快速算法，所述快速傅里叶变换(FFT)采用频率抽取算 法，在频域内把序列按照奇偶分组并利用周期性和对称性进行计算。
[0031] -种基于峰值指标的旋转机械故障诊断方法，进一步的，快速傅里叶逆变换 (IFFT)通过Y(k)频域中的频谱，将每个频率分量变换成时域正弦波，再全部叠加得到相关 函数I(t)。
[0032] -种基于峰值指标的旋转机械故障诊断方法，进一步的，步骤（1)和步骤（4)按 1024个点为一组进行采样，采样频率为1000Hz。
[0033] 与现有技术相比，本发明的优点在于：
[0034] 1.本发明提供的一种基于峰值指标的旋转机械故障诊断方法，通过标准振动信号 将实时采集的旋转机械振动信号分为无故障振动信号和混合信号，其中混合信号包含了故 障特征信号和高斯噪声。用混合信号和无故障振动信号构建峰值指标，使用了信号分离的 思想构建峰值指标，由于对标准振动信号的分离，剩余的混合信号能够更好地体现旋转机 械的故障，对故障的诊断更加地敏感，与现有峰值指标相比，更加能够通过信号的变化检测 微弱故障特征信号。
[0035] 2.本发明提供的一种基于峰值指标的旋转机械故障诊断方法，通过一系列的数学 运算，使得旋转机械设备正常运行和发生故障时，得到的峰值指标的取值范围相对重叠少， 不同状态下峰值指标变化明显，对旋转机械轴承的局部剥落，压痕等故障非常敏感，且不 受振动信号的绝对水平影响，不容易产生误判，特别对偏心轴故障比较敏感，能够很好地将 此故障与其他故障区分开来，反应快、敏感性好，旋转机械故障诊断的精确度和可靠性较 尚。
[0036] 3.本发明提供的一种基于峰值指标的旋转机械故障诊断方法，采用磨合期结束后 的振动信号归一化后作为标准振动信号，考虑了旋转机械设备在运行初期振动信号中噪声 严重，不适宜作为标准振动信号的问题，故障信号反应灵敏，诊断故障更加准确有效。
【附图说明】
[0037] 图1是本发明提供的一种基于峰值指标的旋转机械故障诊断方法的流程图。
[0038]
【具体实施方式】
[0039] 下面结合附图，对本发明提供的一种基于峰值指标的旋转机械故障诊断方法的技 术方案进行进一步的描述，使本领域的技术人员可以更好的理解本发明并能予以实施。
[0040] 参见图1和表1，本发明提供一种基于峰值指标的旋转机械故障诊断方法，包括以 下步骤：
[0041] (1)通过振动加速度传感器采集旋转机械在运行磨合期后的正常运行参数，正常 运行参数为按抽样频率fs采集的无故障振动信号SQ(t)(t = 0，l，…，T-1)，按1024个点为一 组进行采样，峰值指标取50组，取峰值指标50组的最小值与最大值作为该峰值指标的取值 范围；
[0042] (2)对无故障振动信号so(t)进行归一化处理后得到标准振动信号s(t)，采用磨合 期后的振动信号归一化后作为标准振动信号，考虑了旋转机械设备在运行初期振动信号中 噪声严重，不适宜作为标准振动信号的问题，故障信号反应灵敏，诊断故障更加准确有效；
[0043] (3)对标准振动信号s(t)做快速傅里叶变换(FFT)得到标准振动频域信号S(k)(k =〇，1，…，κ-l);
[0044] (4)通过振动加速度传感器采集旋转机械需要故障诊断时的运行参数，所述故障 诊断时的运行参数为按抽样频率f s实时采集旋转机械运行中的实时采集振动信号z(t)(t =〇，1，…，T-1);
[0045] (5)对实时采集振动信号z(t)做快速傅里叶变换(FFT)得到实时振动频域信号Z (k)(k = 0，l，…，K-1);
[0046] (6)对标准振动频域信号S(k)求共辄复数为S(k)%共辄复数的两个实部相等，虚 部互为相反数。
[0047] 当虚部不为零时，共辄复数就是实部相等，虚部相反，如果虚部为零，其共辄复数 就是自身。将z(k)与s(kf相乘得到Y(k)(k=o，i，…，κ-1)，再对γ(ι〇α=ο，ι，···，κ-ι)做快 速傅里叶逆变换(IFFT)得到两个信号Z(k)与S(kf的相关函数l(t)(t = 0，l，…，T-1);
[0048] (7)在t = 0，1，…，T-1中取相关函数I (t)的模11 (t) I，I I (t) I的最大值所对应的时 间点为实时采集振动信号和标准振动信号的延迟时间T，T = argmax| I(t) I ;
[0049] (8)求实时采集振动信号z(t)与延迟后的标准振动信号s(t-T)的相关系数c，cSz (t)s(t-T)的数学期望， C = E[z(t)s(t-T)];
[0050] (9)计算Z(t)-cs(t-T)得到7(〇, 7(〇为故障特征提取信号和噪声的混合信号； [0051 ]实时米集振动彳目号z(t) =cs(t_T)+x(t)+v(t)，时间τ为实时振动彳目号和标准彳目号 的延迟时间，V(t)为高斯噪声，x(t)为故障特征信号，c为相关系数，y(t)为故障特征提取信 号和噪声的混合信号，y(t)=x(t)+u(t)，实时采集振动信号z(t)= CS(t-T)+y(t);
[0052] 通过标准振动信号s(t)将实时采集振动信号z(t)分为无故障振动信号和混合信 号y(t)，混合信号包含了故障特征信号x(t)和高斯噪声v(t)，用混合信号y(t)和无故障振 动信号构建峰值指标，具有对故障诊断更加灵敏，与现有峰值指标相比，更加能够通过信号 的变化检测微弱故障特征信号；
[0053] (10)计算出峰值指标
，E为数学期望；
[0054] 作为一种优选方案，本发明提供的一种基于峰值指标的旋转机械故障诊断方法， 根据峰值指标Cys的值，判断旋转机械是否发生故障及故障类别，以下区间值包括边界值： [0055]峰值指标C ys的值在2.1986到3.7719之间时，旋转机械运行正常；
[0056]峰值指标Cys的值在3.8474到6.0233之间时，考虑旋转机械发生了裂轴故障；
[0057]峰值指标Cys的值在3.8830到6.1981之间时，考虑旋转机械发生了弯轴故障；
[0058]峰值指标Cys的值在3.9094到7.5683之间时，考虑旋转机械发生了偏心轴故障； [0059]峰值指标Cys的值在4.1663到7.3144之间时，考虑旋转机械发生了裂轴+弯轴故障； [0060]峰值指标C ys的值在3.9924到7.4869之间时，考虑旋转机械发生了裂轴+偏心轴故 障；
[0061 ]峰值指标Cys的值在4.5938到6.2757之间时，考虑旋转机械发生了弯轴+偏心轴故 障；
[0062]峰值指标Cys的值在4.0169到6.1841之间时，考虑旋转机械发生了裂轴+弯轴+偏心 轴故障。
[0063]表 1
[0065] 本发明提供的一种基于峰值指标的旋转机械故障诊断方法，旋转机械设备正常运 行和发生故障时，峰值指标变化明显，对旋转机械轴承的局部剥落，压痕等故障非常敏感， 且不受振动信号的绝对水平影响，不容易产生误判，特别对偏心轴故障比较敏感，能够很好 地将此故障与其他故障区分开来，反应快、敏感性好，旋转机械故障诊断的精确度和可靠性 较高。
[0066] 作为一种优选方案，本发明提供的一种基于峰值指标的旋转机械故障诊断方法， 快速傅里叶变换(FFT)采用有限序列离散傅里叶变换(DFT)的快速算法，快速傅里叶变换 (FFT)采用频率抽取算法，在频域内把序列按照奇偶分组并利用周期性和对称性进行计算。
[0067] 作为一种优选方案，本发明提供的一种基于峰值指标的旋转机械故障诊断方法， 快速傅里叶逆变换（IFFT)通过Y( k)频域中的频谱，将每个频率分量变换成时域正弦波，再 全部叠加得到相关函数I(t)。快速傅里叶变换大大提升了计算机的运算效率，减少了运算 次数。离散傅里叶变换以及逆变换分别如下：
[0070] 其中0<k<N-l，令W=eTj2VN，则N点序列的离散傅里叶变换为：
[0071] X(k) = x{n)W''" (0 < k <N-\) n~Q
[0072] 具有周期性:wkn=wn(k+N)=w k(n+N)
[0073] Wkn 具有对称性:Wkn=-Wkn+N2
[0074] 通过周期性和对称性简化了离散傅里叶变换。
[0075] 作为一种优选方案，本发明提供的一种基于峰值指标的旋转机械故障诊断方法， 步骤(1)和步骤(4)按1024个点为一组进行采样，采样频率为1000Hz。
[0076] 本发明提供的一种峰值指标在不同故障下的敏感程度不同，对偏心轴故障比较敏 感，能够很好地将此故障与其他故障区分开来，能很好的反映了不同故障特征信息，构建的 峰值指标是以比值构成，具有不受机器工况影响的特点，对旋转机械轴承的局部剥落，压痕 等故障非常敏感，且不受振动信号的绝对水平影响，不容易产生误判。
[0077]上述实施方式为本发明较佳的实施方式，但本发明的实施方式并不受上述实施方 式的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简 化，均应为等效的置换方式，都包含在本发明的保护范围之内。
【主权项】
1. 一种基于峰值指标的旋转机械故障诊断方法，其特征在于，包括以下步骤： (1) 通过振动加速度传感器采集旋转机械在运行磨合期后的正常运行参数，所述正常 运行参数为按抽样频率fs采集的无故障振动信号so(t) (t = O，1，…，T-I); (2) 对无故障振动信号so(t)进行归一化处理后得到标准振动信号s(t); (3) 对标准振动信号s(t)做快速傅里叶变换得到标准振动频域信号S(k)(k = 0，l，…， K-I )； (4) 通过振动加速度传感器采集旋转机械需要故障诊断时的运行参数，所述故障诊断 时的运行参数为按抽样频率fs实时采集旋转机械运行中的实时采集振动信号z(t)(t = 0， 1，…，T-1); (5) 对实时采集振动信号z(t)做快速傅里叶变换得到实时振动频域信号Z(k)(k = 0， 1，…，K-1); (6) 对标准振动频域信号S(k)取复共辄为S(kf，将Z(k)与S(k)$相乘得到Y(k)(k = 0， 1，…，K-1)，再对Y(k) (k = 0,1，…，κ-1)做快速傅里叶逆变换得到两个信号Z(k)与S(kf的 相关函数I(t)(t = 0，l,…，τ-1); (7) 在t=〇，1，…，T-I中取相关函数I (t)的模11 (t) I，11 (t) I的最大值所对应的时间点 为实时采集振动信号和标准振动信号的延迟时间T，T = argmax| I(t) I ; (8) 求实时采集振动信号z(t)与延迟后的标准振动信号s(t-T)的相关系数c，c = E[z (t)s(t-T)]； (9) 计算混合信号y(t)，y(t)为故障特征提取信号和噪声的混合信号，y(t) = z(t)-CS (t-τ)； (10) 计算出峰值指相2. 根据权利要求1所述的一种基于峰值指标的旋转机械故障诊断方法，其特征在于，根 据所述峰值指标Cys的值，判断旋转机械是否发生故障及故障类别，以下区间值均包括边界 值： 峰值指标Cys的值在2.1986到3.7719之间时，旋转机械运行正常； 峰值指标Cys的值在3.8474到6.0233之间时，考虑旋转机械发生了裂轴故障； 峰值指标Cys的值在3.8830到6.1981之间时，考虑旋转机械发生了弯轴故障； 峰值指标Cys的值在3.9094到7.5683之间时，考虑旋转机械发生了偏心轴故障； 峰值指标Cys的值在4.1663到7.3144之间时，考虑旋转机械发生了裂轴+弯轴故障； 峰值指标Cys的值在3.9924到7.4869之间时，考虑旋转机械发生了裂轴+偏心轴故障； 峰值指标Cys的值在4.5938到6.2757之间时，考虑旋转机械发生了弯轴+偏心轴故障； 峰值指标Cys的值在4.0169到6.1841之间时，考虑旋转机械发生了裂轴+弯轴+偏心轴故 障。3. 根据权利要求1所述的一种基于峰值指标的旋转机械故障诊断方法，其特征在于:所 述步骤(3)和步骤(5)中的快速傅里叶变换采用有限序列离散傅里叶变换的快速算法，所述 快速傅里叶变换采用频率抽取算法，在频域内把序列按照奇偶分组并利用周期性和对称性 进行计算。4. 根据权利要求1所述的一种基于峰值指标的旋转机械故障诊断方法，其特征在于:所 述步骤(7)中的快速傅里叶逆变换通过Y(k)频域中的频谱，将每个频率分量变换成时域正 弦波，再全部叠加得到相关函数I(t)。5. 根据权利要求1或4所述的一种基于峰值指标的旋转机械故障诊断方法，其特征在 于:所述步骤(1)和步骤(4)按1024个点为一组进行采样，采样频率为1000Hz。
【文档编号】G01M13/02GK106017921SQ201610560721
【公开日】2016年10月12日
【申请日】2016年7月13日
【发明人】张清华, 王子为, 何俊, 谭巨兴, 肖明
【申请人】广东石油化工学院