基于变分模态分解-传递熵的轴承退化状态识别预测方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN109580222A(43)申请公布日2019.04.05(21)申请号201811473007.X(22)申请日2018.12.04(71)申请人河北科技大学地址050018河北省石家庄市裕翔街26号(72)发明人武哲张嘉钰崔彦平常宏杰张付祥张新聚牛虎利(74)专利代理机构石家庄冀科专利商标事务所有限公司13108代理人曹淑敏(51)Int.Cl.G01M13/04(2019.01)权利要求书3页说明书5页附图2页(54)发明名称基于变分模态分解-传递熵的轴承退化状态识别预测方法(57)摘要本发明涉及一种基于变分模态分解-传递熵的轴承退化状态识别预测方法，属于齿轮故障分析技术领域。技术方案是：通过系统非线性和复杂程度的参数有效地反映旋转机械故障的严重程度；将变分模态分解与基于非线性动力学参数的信号复杂性的传递熵理论相结合，实现滚动轴承退化状态识别；建立基于变分模态分解-传递熵的滚动轴承状态评价指标，应用模型对滚动轴承全寿命试验数据进行预测。本发明为旋转机械的故障诊断、性能退化状态识别和趋势预测提供新的有效手段；建立变分模态分解-传递熵和SVR的滚动轴承故障演化趋势预测模型，应用该模型对滚动轴承全寿命试验数据进行预测，提高准确性和有效性。CN109580222ACN109580222A权利要求书1/3页1.基于变分模态分解-传递熵的轴承退化状态识别预测方法，其特征在于包含如下步骤：当机械设备发生故障时，系统的动力学行为表现出强非线性和非稳定性，通过系统非线性和复杂程度的参数有效地反映旋转机械故障的严重程度；针对旋转机械性能退化识别指标的建立问题，将变分模态分解与基于非线性动力学参数的信号复杂性的传递熵理论相结合，实现滚动轴承退化状态识别；定量描述轴承时间序列信号频段间非线性的耦合及信息传递特征，得到轴承故障振动信号的传递路径，建立基于变分模态分解-传递熵的滚动轴承状态评价指标，进行旋转机械的故障诊断、性能退化状态识别和趋势预测；针对滚动轴承退化趋势预测问题，建立变分模态分解-传递熵和SVR的滚动轴承故障演化趋势预测模型，应用该模型对滚动轴承全寿命试验数据进行预测。2.根据权利要求1所述的基于变分模态分解-传递熵的轴承退化状态识别预测方法，其特征在于具体步骤如下：步骤一：数据采集；(1)读取轴承振动信号x的时序序列x(t)，其中，t＝1,2,3…,N,N为轴承信号的采样点个数；(2)对轴承振动信号x(t)进行采用VMD分解为一系列变分模态分量，使各个模态的估计带宽之和最小，约束条件为各模态之和等于输入信号x(t)，约束变分模型描述为式(1)和(2)；式中：{uk}＝{u1,…,uk}表分解得到的K干IMF分量；{ωk}＝{ω1,…,ωk}表示各分量的中心频率；表示所有模态分量的求和，*表示卷积；(3)为求取上述约束变分问题，引入二次惩罚因子α和Lagrange乘法算子λ(t)，其中二次惩罚因子可在高斯噪声存在的情况下保证信号的重构精度，Lagrange约束条件保持严格性，扩展的Lagrange表达式如下：(4)利用交替方向乘子算法求取上述增广Lagrange函数的最优解，从而将原始信号分解为K个窄带IMF分量；VMD大致过程如下：1)将λ和n初始化；2)n＝n+1，执行整个循环；3)执行内层第一个循环，根据更新uk；2CN109580222A权利要求书2/3页4)k＝k+1重复步骤(3)，直至k＝K，结束内层第一个循环；5)执行内层第二个循环，根据6)k＝k+1重复步骤(5)，直至k＝K，结束内层第二个循环；7)根据更新λ；8)重复上述步骤(2)～(7)，直至满足迭代停止条件结束整个循环，最终得到K个窄带IMF分量；步骤二：使用量子粒子群优化算法确定此单通道轴承振动信号的变分模态分解算法的分量个数和惩罚因子两个参数；(1)基于PSO(粒子群算法)建立变分模态分解算法的分量个数和惩罚因子两个参数的优化模型，其模型如下：vij(t+1)＝wvij(t)+c1(t)(pij(t)-xij(t))+c2r2j(t)(pgj(t)-xij(t))xij(t+1)＝xij(t)+vij(t+1)式中：i＝1，2，…，m，m是该群体中粒子的数目j表示粒子的维数，即算法中所优化的第j个参数；t表示进化代数；vij(t)表示t时刻粒子；i在j维的空间速度；w为惯性权重，它使微粒保持运动的惯性，使其有能力探索新的区域；rij为[0，1]范围内变化的随机数；pij为粒子i的历史最好解(个体最优)的j维值；pgj为所有粒子的历史最好解(群体最优)的j维值；c1和c2为速度常数，它们使每个粒子向pij和pgj位置加速运动；xij(t)为t时刻粒子i在j维的空间位置，即第j个优化