鲁棒性幅值指数自适应谱幅值调制法_林云.pdf

年第 卷 月第 期机 械 科 学 与 技 术 ：收稿日期：基金项目：国家自然科学基金项目（）作者简介：林云（），硕士研究生，研究方向为旋转机械故障特征提取，通信作者：郭瑜，教授，博士生导师，林云，郭瑜，刘珍鲁棒性幅值指数自适应谱幅值调制法机械科学与技术，（）：鲁棒性幅值指数自适应谱幅值调制法林云，郭瑜，刘珍（昆明理工大学 机电工程学院，昆明）摘要：近期提出的谱幅值调制（）法可通过调节幅值指数适应不同故障信号的特征提取，具有较好的应用前景。但该方法目前幅值指数的取值需人工判断，尚无法用其实现对故障特征的自动优化提取，且当故障特征受到复杂干扰时，通过人工选择较难选取较佳的幅值指数。为此，本文研究提出了一种鲁棒性幅值指数自适应谱幅值调制法。该方法首先利用角域重采样将信号转换到角域，再通过多点最优最小熵解卷积（）对故障弱冲击特征进行增强，最后利用 阶循环平稳指标（）自适应选取 中的倒谱幅值指数，以该优化幅值指数计算倒谱信号，实现故障特征的自动提取。在行星轴承内圈故障特征提取上进行了验证研究，实验结果表明，本文所提方法能够实现复杂干扰下行星轴承内圈故障特征的自适应提取。关键词：谱幅值调制；鲁棒性；自适应；阶循环平稳指标；行星轴承中图分类号：.文献标志码：文章编号：（），（，）：（），（），（），：；齿轮箱故障引起的振动信号通常为多源耦合信号，为提取故障特征，常采用信号处理技术将信号中包含的特征分量有效分离出来。目前常用的信号分离方法主要有同步平均、自参考自适应消噪技术（）、离散随机分离（）和倒谱预白化（）等，其中 不需要针对不同的故障情况提前设置参数即可有效降低齿轮信号等离散分量的影响，在应用中取得了较好的效果。但值得指出的是，该技术将倒谱幅值置零来提升故障冲机 械 科 学 与 技 术第 卷：击成分的同时，会无差别地提升噪声成分的幅值。等近期提出的谱幅值调制（，）通过可变幅值指数弥补了 的不足，可以较灵活地适应不同的故障情况。尽管 技术在 基础上作了很好的改进，但目前 技术在选取幅值指数值上仍需人工判断，不适合实现对故障特征的自动提取，且当故障特征受到复杂干扰时，人工选择较难选取出较佳的幅值指数，易造成误诊。为实现在干扰环境下 对故障特征的自适应提取，本文提出一种鲁棒性幅值指数自适应 法。该方法先利用角域重采样将信号转换到角域，通过（，）对故障弱冲击特征进行增强，再利用 阶循环平稳指标（）自适应选取 中的倒谱幅值指数，以该优化幅值指数计算倒谱信号，实现故障特征的自动提取。研究中通过行星轴承内圈故障特征提取实验验证了本文方法的有效性。谱幅值调制法简介谱幅值调制法（）利用原始信号实倒谱和相位谱来重构信号，通过设置不同的实倒谱幅值的幂次来增强信号中的特定周期成分。该方法本质上相当于一个非线性滤波器，具有较好的鲁棒性，不需要选择优化解调频带也可以有效提取轴承故障特征。该方法是在倒谱预白化的基础上提出的。不同的是，该方法中实倒谱幅值的幂次是可变的。首先利用实倒谱和信号相位，得到修正信号，其表达式为（，）（）（）（）式中 为幅值幂次，取值范围为.。（，）（，）（，）（）（，）（，）（）式中：表示修正信号对应的解析信号；表示修正信号的希尔伯特变换；表示取绝对值。通过式（）求修正信号的平方包络谱。此方法包含了两个主要步骤：）通过设置不同范围的 值对信号产生不同的滤波效果。当.时，原始信号中振幅较大的频率成分被放大。随着 的增加，幅振幅较小的频率成分会被包含更高幅值水平频率的部分逐渐掩盖；当.时，情况则相反；当 时，振幅较大的频率成分被一定程度抑制。而振幅较小的频率成分被一定程度增强，成为占据主导地位的信号成分。）通过归一化平方包络谱降噪，提取具有较大振幅信号中的冲击特征。简介假设传感器采集到的振动信号为 （）式中：为传递函数；为故障冲击序列；为噪声。解卷积的目的是获得合适的滤波器对信号 进行滤波，最大程度恢复故障冲击相关成分。其原理简述如下：首先，定义多点 范数，即（，）（）式中目标向量 是用来定义解卷积目标冲击位置和权重的连续向量，该向量的大小与解卷积信号大小相同，冲击的位置和权重可分别通过其元素位置和大小来反映。该算法以多点 范数最大化为目标，即：（，）（）通过对滤波器组 求 阶导数的方式求式（）的极值，同时假设（）存在，最终可求得优化滤波器组 和输出结果分别为：（）（）（）算法的详细介绍可参考文献，本文不再赘述。鲁棒性幅值指数自适应谱幅值调制法 的优势在于通过调节幅值指数 可以非常灵活适应不同故障程度及类型的振动信号。但目前幅值指数 的取值需要人工根据谱幅值调制图综合各谐波的能量水平来判断，无法实现其对滚动轴承故障特征的自动优化提取。另一方面，在旋转机械故障诊断中，发展故障特征的自适应提取是实现故障诊断自动化的重要保障。为实现 中幅值指数 的自适应选取，需要利用某一指标有效地反映信号中包含故障冲击成分的情况。常规的峭度指标易受信号中随机冲击噪声等因素的干扰；文献提出的相关峭度较好克服了峭度易受随机冲击影响的缺点，但轴承故障产生的冲击在时域并非严格的周期信号而具有较强的 阶循环平稳性。因此，近年来 阶循环平稳指标（）被提出并在滚动轴承故障特征提取中取得较好效果。据此，本文选用 作为幅值指数自适应的评价指标，其定义为：第 期 林云，等：鲁棒性幅值指数自适应谱幅值调制法：，（）（）（）（）式中：（）为高斯平稳信号；为循环频率（预设为故障频率）；为信号 的均方功率；为选取的故障特征谐波数。由于轴承存在滑移，因此该式中计算 时需要设置一个容错区间。研究表明，将该指标应用于式（）计算所得的信号自适应获得幅值指数 可取得较好效果。利用 技术可有效抑制无关周期信号的干扰，并提升故障冲击。但研究中发现，对于受时变传递路径影响的行星齿轮箱多源耦合信号，并不直接适用（见本文图），需要对信号进行前处理来增强故障冲击特征。是一种较为理想的信号前处理方法，其优势主要包括两个方面：）以多点 范数最大为目标进行求解可以有效降低随机冲击带来的影响；）采用非迭代方式求解滤波器，可以避免陷入局部最优解。总体而言，该方法通过获取周期冲击的全局最优解，可以有效增强故障冲击特征。基于上述分析，本文在对 进行改进的基础上，结合 技术，提出鲁棒性幅值指数自适应谱幅值调制法。该方法的技术路线如图 所示。图 鲁棒性幅值指数自适应谱幅值调制法技术路线基本步骤为：）采集获得原始时域信号，通过等角度重采样得到角域信号，消除转速波动的影响。）选取合适的滤波器长度 和周期点数范围，利用 识别并提升角域信号中行星轴承内圈故障对应的微弱冲击特征。）经过 反卷积滤波后的信号作为 的输入信号，设定 值取值范围为.，.，研究中取步长为.，分别求修正信号。计算每个 值对应的修正信号的 值，取故障特征谐波数。输出其中最大的 对应的 值，为优化 值。）该优化 值自动输入 程序，得到其对应的归一化平方包络阶比谱，通过包络谱提取故障特征。行星轴承内圈故障特征提取实验 实验简介研究中用如图 所示的 型行星齿轮箱综合试验台验证本文方法的有效性。实验中使用 型加速度传感器（灵敏度：.）及配套 电荷放大器拾取振动信号；使用 型电涡流传感器（灵敏度：.）获取转速信号。使用 采集卡采集信号。图 行星齿轮箱综合试验台传感器的具体安装位置如图 所示，其中电涡流传感器安装位置为图 中 处，加速度传感器共 个，安装位置分别为图 中、处。因位置 处的传感器采集的振动信号传递路径最短，本研究中使用该处采集的振动信号。图 传感器安装位置行星轴承故障为电火花加工的人造故障，加工宽度约为 ，深度约为.，用以模拟轴承局部故障，如图 所示。机 械 科 学 与 技 术第 卷：图 行星轴承内圈人造故障行星齿轮箱 内所用行星轴承、齿轮参数分别如表、表 所示。在输入端（对应太阳轮）转速为 时，以行星架为参考，利用表 和表 相关参数，由行星轮系的传动理论，计算得到的行星轴承内圈故障相关阶次如表 所示。表 行星轴承参数滚子数 滚子直径 节圆直径 接触角 （）表 齿轮参数齿轮 太阳轮行星轮（个）齿圈齿数 模数.变位系数.表 行星轴承内圈故障相关阶次参数数值太阳轮绝对旋转阶次.行星架旋转阶次.行星轮绝对旋转阶次.行星轮轴承内圈故障特征阶次.行星轮轮齿啮合阶次 故障特征提取行星轴承内圈发生故障时，采集到的振动信号及实际转速曲线如图 所示，从图 中可以看出在该工况下，转速存在微小波动。为消除转速波动的影响，以行星架转轴为参考轴，对原始振动信号进行等角度重采样，将信号转换到角域。本文设置角域重采样率 点，为固定齿圈齿数的 倍。图 内圈故障行星轴承振动等角度重采样后得到的原始角域信号如图 所示，其包络阶比谱如图 所示。图 角域信号图 内圈故障行星轴承振动的包络阶比谱从图 中可见，太阳轮绝对旋转阶次（.）、行星轮绝对旋转阶次（.）和行星架旋转阶次（.）相对突出，而行星轴承内圈故障特征阶次（.）完全被淹没。直接使用 识别并增强原始角域信号中行星轴承内圈故障冲击。计算得行星轴承故障理论周期点数为 ，其计算公式为 （）式中 为行星轴承内圈故障特征阶次。为了 减 小 与 故 障 无 关 的 信 号 干 扰，应 用 技术时，周期点数范围应尽可能取小，但考虑到滚动轴承在工作过程中存在 的滑移，周期点数范围过小会导致故障信息利用不足。综合考虑，选取周期点数范围为（，），而第 期 林云，等：鲁棒性幅值指数自适应谱幅值调制法：滤波器长度 通常取在 到 范围内较为合适，本文取。图 所示为行星轴承内圈故障下振动信号经过 冲击增强后的信号包络阶比谱。从图中可见，低阶次频段仍受到噪声干扰，难以有效获得故障特征信息，需要进一步对信号进行处理。图 内圈故障行星轴承振动的包络阶比谱经过 处理后的信号，复杂调制对故障冲击造成的影响难以被有效降低。后续经过 处理后的谱幅值调制图如图 所示。图 内圈故障行星轴承振动的谱幅值调制图 中在.时可观察到故障特征谱线。需要指出的是，当 的取值在此区间内时，各阶谐波能量最大时对应的 值并不相同，需要考虑各阶谐波情况进行综合判断，通过人工观察无法实现特征提取的自动化且较难选取优化的 值。采用本文提出的方法，用 对 值进行自适应选取，行星轴承内圈故障信号对应的 值随 值的变化情况如图 所示，可知 .最合适。结合图 所示结果，说明了该指标能够较好地反映故障特征阶次及各阶谐波综合能量水平随 值的变化情况。图 内圈故障行星轴承振动的 曲线利用 优化 处理后得到的行星轴承内圈故障振动对应的信号平方包络阶比谱如图 所示。需要说明的是，因轴承滑移现象的存在，实际故障特征阶次与理论值存在约.的偏差，但不影响故障特征提取。从图 中可见较为突出的行星轴承内圈故障特征阶次及其谐波，而无关周期信号被抑制，表明本文所提方法可以有效提取行星轴承故障信号的冲击特征。图 内圈故障行星轴承振动的平方包络阶比谱（.）对比实验为说明利用 对信号进行前处理的必要性，直接使用 对同一原始信号进行处理，结果如图 和图 所示。从图 中可见，当 在.，.范围变化时，齿轮啮合阶次及谐波始终占据主导地位，而行星轴承内圈故障特征不易观察到。进一步地，求特定 值时的平方包络谱。根据文献中理论，在齿轮啮合信号干扰时，取 值在（.，）范围较为合适，本文取.得到如图 所示的平方包络阶比谱，图中可见轴承故障特征被噪声淹没，难以读取有用的故障信息。结合图 和图 可知，直接通过 不易有效提取行星轴承故障特征。图 内圈故障行星轴承振动的谱幅值调制图 内圈故障行星轴承振动的平方包络阶比谱（.）机 械 科 学 与 技 术第 卷：结论本文针对 需要人工判断幅值指数 值，尚无法用其实现对故障特征的自动优化提取的不足，提出了鲁棒性幅值指数自适应谱幅值调制法，使用自适应地选取 值。结合 可以有效增强故障冲击特征的优势，该方法先将 技术用于故障振动信号的前处理，以增强故障冲击特征，再利用参数自适应 技术减小无关周期信号的干扰，提取故障特征。通过行星轴承内圈故障的实验，验证了本文方法的有效性。参考文献 李杰，赵建民基于时域同步平均与分离技术的齿轮箱振