一种基于裕度指标的旋转机械故障诊断方法技术

本发明专利技术提供一种基于裕度指标的旋转机械故障诊断方法，用信号分离的思想构建裕度指标，克服现有裕度指标诊断技术的缺点，提高诊断的精确度和可靠性，通过标准振动信号将实时采集振动信号分为无故障振动信号和混合信号，用混合信号和无故障振动信号构建裕度指标，具有对故障诊断更加灵敏，更加能够通过信号的变化检测微弱故障特征信号。旋转机械设备正常运行和发生故障时，裕度指标的取值范围重叠少，不同状态下裕度指标变化明显，对偏心轴加弯轴类复合故障比较敏感，能够很好地将此故障与其他故障区分开来，旋转机械故障诊断的精确度和可靠性高。

【技术实现步骤摘要】

本专利技术涉及一种旋转机械故障诊断方法，特别涉及一种基于裕度指标的旋转机械故障诊断方法，属于故障诊断与信号处理分析

技术介绍
大型旋转机械设备(如汽轮机、旋转轴承、压缩机等)是石油、化工、机械制造、航空航天等重要工程领域的关键设备，旋转机械设备正朝着大型化、自动化、精密化的方向不断发展，其组成和结构也变得越来越复杂，发生故障的概率也越来越大，因此对大型旋转机械设备的安全性和可靠性的要求也越来越高。但是，大型旋转机械设备发生故障时，振动监测信号往往存在大量的非线性、随机、不可遍历的信息，给故障信号的分析和旋转机械故障的诊断带来很大的困难。现有技术的旋转机械故障诊断主要基于振动信号进行分析，而一般采用时域分析法，通过对机械振动信号的概率密度函数分析，推导出了幅值域中的有量纲指标和无量纲指标，有量纲指标如均值、均方根值等；无量纲指标如裕度指标、裕度指标、裕度指标等。在实际应用中，现有技术的有量纲指标对故障特征敏感，其数值会随着故障的发展而上升，同时因工作条件，如负载、转速等的变化而变化，并极易受环境干扰的影响，表现不够稳定。无量纲指标中的裕度指标对于振动监测信号中的扰动不敏感，能对中、高频故障的机械设备做出预判。特别地，裕度指标对信号的幅值和频率的变化不敏感，与旋转机械的工作条件关系不大，只依赖于概率密度函数的形状。因此，裕度指标在旋转机械故障诊断中得到了广泛的应用。但是，现有技术利用裕度指标来诊断旋转机械的故障主要有以下几点问题：一是没有采用信号分离的思想构建裕度指标，没有对标准振动信号分离，无法更好地体现旋转机械的故障，对故障的诊断不敏感；二是在利用裕度指标进行分析时，裕度指标随着故障的逐渐发展，敏感度反而会下降，稳定性不好，难以诊断和预测故障，设备正常运行和发生故障时的取值范围会部分重叠，易产生误判，旋转机械故障诊断的精确度和可靠性较差；三是采用磨合期的振动信号作为标准振动信号，没有考虑旋转机械设备在运行初期振动信号中噪声严重，不适宜作为标准振动信号的问题，诊断故障不准确。
技术实现思路
针对现有技术的不足，本专利技术提供一种基于裕度指标的旋转机械故障诊断方法，克服传统裕度指标诊断技术的缺点，提高诊断的精确度和可靠性，用信号分离的思想构建裕度指标，由于对标准振动信号的分离，剩余的混合信号能够更好地体现旋转机械的故障，对故障的诊断更加地敏感，所得出的诊断结果能够准确地诊断旋转机械的故障。为达到以上技术效果，本专利技术所采用的技术方案如下：一种基于裕度指标的旋转机械故障诊断方法，包括以下步骤：(1)通过振动加速度传感器采集旋转机械在运行磨合期后的正常运行参数，所述正常运行参数为按抽样频率fs采集的无故障振动信号s0(t)(t＝0,1,…,T-1)；(2)对无故障振动信号s0(t)进行归一化处理后得到标准振动信号s(t)；(3)对标准振动信号s(t)做快速傅里叶变换(FFT)得到标准振动频域信号S(k)(k＝0,1,…,K-1)；(4)通过振动加速度传感器采集旋转机械需要故障诊断时的运行参数，所述故障诊断时的运行参数为按抽样频率fs实时采集旋转机械运行中的实时采集振动信号z(t)(t＝0,1,…,T-1)；(5)对实时采集振动信号z(t)做快速傅里叶变换(FFT)得到实时振动频域信号Z(k)(k＝0,1,…,K-1)；(6)对标准振动频域信号S(k)取复共轭为S(k)*，将Z(k)与S(k)*相乘得到Y(k)(k＝0,1,…,K-1)，再对Y(k)(k＝0,1,…,K-1)做快速傅里叶逆变换(IFFT)得到两个信号Z(k)与S(k)*的相关函数I(t)(t＝0,1,…,T-1)；(7)在t＝0,1,…,T-1中取相关函数I(t)的模|I(t)|，|I(t)|的最大值所对应的时间点为实时采集振动信号和标准振动信号的延迟时间τ，τ＝argmax|I(t)|；(8)求实时采集振动信号z(t)与延迟后的标准振动信号s(t-τ)的相关系数c，c＝E[z(t)s(t-τ)]；(9)计算混合信号y(t),y(t)为故障特征提取信号和噪声的混合信号，y(t)＝z(t)-cs(t-τ)；(10)计算出裕度指标一种基于裕度指标的旋转机械故障诊断方法，进一步的，根据所述根据裕度指标CLys的值，判断旋转机械是否发生故障及故障类别，以下区间值均包括边界值：裕度指标CLys的值在4.8008到6.8580之间时，旋转机械运行正常；裕度指标CLys的值在7.6329到10.1778之间时，考虑旋转机械发生了裂轴故障；裕度指标CLys的值在7.1511到11.4458之间时，考虑旋转机械发生了弯轴故障；裕度指标CLys的值在7.0263到14.6501之间时，考虑旋转机械发生了偏心轴故障；裕度指标CLys的值在11.6733到13.9574之间时，考虑旋转机械发生了裂轴+弯轴故障；裕度指标CLys的值在6.8843到11.3631之间时，考虑旋转机械发生了裂轴+偏心轴故障；裕度指标CLys的值在7.3568到18.8755之间时，考虑旋转机械发生了弯轴+偏心轴故障；裕度指标CLys的值在9.8834到16.4384之间时，考虑旋转机械发生了裂轴+弯轴+偏心轴故障。一种基于裕度指标的旋转机械故障诊断方法，进一步的，快速傅里叶变换(FFT)采用有限序列离散傅里叶变换(DFT)的快速算法，所述快速傅里叶变换(FFT)采用频率抽取算法，在频域内把序列按照奇偶分组并利用周期性和对称性进行计算。一种基于裕度指标的旋转机械故障诊断方法，进一步的，快速傅里叶逆变换(IFFT)通过Y(k)频域中的频谱，将每个频率分量变换成时域正弦波，再全部叠加得到相关函数I(t)。一种基于裕度指标的旋转机械故障诊断方法，进一步的，步骤(1)和步骤(4)按1024个点为一组进行采样，采样频率为1000Hz。与现有技术相比，本专利技术的优点在于：1.本专利技术提供的一种基于裕度指标的旋转机械故障诊断方法，通过标准振动信号将实时采集的旋转机械振动信号分为无故障振动信号和混合信号，其中混合信号包含了故障特征信号和高斯噪声。用混合信号和无故障振动信号构建裕度指标，使用了信号分离的思想构建裕度指标，由于对标准振动信号的分离，剩余的混合信号能够更好地体现旋转机械的故障，对故障的诊断更加地敏感，与现有裕度指标相比，更加能够通过信号的变化检测微弱故障特征信号。2.本专利技术提供的一种基于裕度指标的旋转机械故障诊断方法，通过一系列的数学运算，使得旋转机械设备正常运行和发生故障时，得到的裕度指标的取值范围相对重叠少，不同状态下裕度指标变化明显，不容易产生误判，对偏心轴加弯轴类复合故障敏感，能够很好地将此故障与其他故障区分开来，旋转机械故障诊断的精确度和可靠性较高。3.本专利技术提供的一种基于裕度指标的旋转机械故障诊断方法，采用磨合期结束后的振动信号归一化后作为标准振动信号，考虑了旋转机械设备在运行初期振动信号中噪声严重，不适宜作为标准振动信号的问题，故障信号反应灵敏，诊断故障更加准确有效。附图说明图1是本专利技术提供的一种基于裕度指标的旋转机械故障诊断方法的流程图。具体实施方式下面结合附图，对本专利技术提供的一种基于裕度指标的旋转机械故障诊断方法的技术方案进行进一步的描述，使本领域的技术人员可以更好的理解本文档来自技高网...

【技术保护点】
一种基于裕度指标的旋转机械故障诊断方法，其特征在于，包括以下步骤：(1)通过振动加速度传感器采集旋转机械在运行磨合期后的正常运行参数，所述正常运行参数为按抽样频率fs采集的无故障振动信号s0(t)(t＝0,1,…,T‑1)；(2)对无故障振动信号s0(t)进行归一化处理后得到标准振动信号s(t)；(3)对标准振动信号s(t)做快速傅里叶变换得到标准振动频域信号S(k)(k＝0,1,…,K‑1)；(4)通过振动加速度传感器采集旋转机械需要故障诊断时的运行参数，所述故障诊断时的运行参数为按抽样频率fs实时采集旋转机械运行中的实时采集振动信号z(t)(t＝0,1,…,T‑1)；(5)对实时采集振动信号z(t)做快速傅里叶变换得到实时振动频域信号Z(k)(k＝0,1,…,K‑1)；(6)对标准振动频域信号S(k)取复共轭为S(k)*，将Z(k)与S(k)*相乘得到Y(k)(k＝0,1,…,K‑1)，再对Y(k)(k＝0,1,…,K‑1)做快速傅里叶逆变换得到两个信号Z(k)与S(k)*的相关函数I(t)(t＝0,1,…,T‑1)；(7)在t＝0,1,…,T‑1中取相关函数I(t)的模|I(t)|，|I(t)|的最大值所对应的时间点为实时采集振动信号和标准振动信号的延迟时间τ，τ＝arg max|I(t)|；(8)求实时采集振动信号z(t)与延迟后的标准振动信号s(t‑τ)的相关系数c，c＝E[z(t)s(t‑τ)]；(9)计算混合信号y(t),y(t)为故障特征提取信号和噪声的混合信号，y(t)＝z(t)‑cs(t‑τ)；(10)计算出裕度指标...

【技术特征摘要】
1.一种基于裕度指标的旋转机械故障诊断方法，其特征在于，包括以下步骤：(1)通过振动加速度传感器采集旋转机械在运行磨合期后的正常运行参数，所述正常运行参数为按抽样频率fs采集的无故障振动信号s0(t)(t＝0,1,…,T-1)；(2)对无故障振动信号s0(t)进行归一化处理后得到标准振动信号s(t)；(3)对标准振动信号s(t)做快速傅里叶变换得到标准振动频域信号S(k)(k＝0,1,…,K-1)；(4)通过振动加速度传感器采集旋转机械需要故障诊断时的运行参数，所述故障诊断时的运行参数为按抽样频率fs实时采集旋转机械运行中的实时采集振动信号z(t)(t＝0,1,…,T-1)；(5)对实时采集振动信号z(t)做快速傅里叶变换得到实时振动频域信号Z(k)(k＝0,1,…,K-1)；(6)对标准振动频域信号S(k)取复共轭为S(k)*，将Z(k)与S(k)*相乘得到Y(k)(k＝0,1,…,K-1)，再对Y(k)(k＝0,1,…,K-1)做快速傅里叶逆变换得到两个信号Z(k)与S(k)*的相关函数I(t)(t＝0,1,…,T-1)；(7)在t＝0,1,…,T-1中取相关函数I(t)的模|I(t)|，|I(t)|的最大值所对应的时间点为实时采集振动信号和标准振动信号的延迟时间τ，τ＝arg max|I(t)|；(8)求实时采集振动信号z(t)与延迟后的标准振动信号s(t-τ)的相关系数c，c＝E[z(t)s(t-τ)]；(9)计算混合信号y(t),y(t)为故障特征提取信号和噪声的混合信号，y(t)＝z(t)-cs(t-τ)；(10)计算出裕度指标2.根据权利要求1所述的一种基于裕度指标的旋转机械故障诊断方法，其...

【专利技术属性】
技术研发人员：肖明，谢侃，雷宝，王子为，谭巨兴，张清华，
申请(专利权)人：广东工业大学，
类型：发明
国别省市：广东;44