CN109883693B - 基于低秩-稀疏逼近的行星齿轮箱故障诊断方法 - Google Patents
基于低秩-稀疏逼近的行星齿轮箱故障诊断方法 Download PDFInfo
- Publication number
- CN109883693B CN109883693B CN201910270414.9A CN201910270414A CN109883693B CN 109883693 B CN109883693 B CN 109883693B CN 201910270414 A CN201910270414 A CN 201910270414A CN 109883693 B CN109883693 B CN 109883693B
- Authority
- CN
- China
- Prior art keywords
- signal
- matrix
- sun
- gear ring
- sparse
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Retarders (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
基于低秩‑稀疏逼近的行星齿轮箱故障诊断方法,先利用内置编码器采集输出轴角度位置信号,然后以均方误差为损失函数,对输出轴角度位置信号拟合并差分,得到瞬时角速度信号,基于瞬时角速度信号,分别对行星轮、太阳轮及齿圈时域同步平均,随后对各自的时域同步平均信号做低秩‑稀疏逼近,得到行星轮、太阳轮及齿圈的稀疏恢复信号,并计算各自的峭度,峭度大于峭度阈值的则认为存在故障;本发明利用低秩‑稀疏逼近方法的特征,对采集到的复合信号解耦,分开行星齿轮箱的正常啮合信息和故障信息,凸显故障特征,有利于实现行星齿轮箱早期微弱故障的诊断。
Description
技术领域
本发明涉及行星齿轮箱故障诊断和信号分析技术领域,特别涉及基于低秩-稀疏逼近的行星齿轮箱故障诊断方法。
背景技术
行星齿轮箱作为重要传动部件,因承载效率高与传动比大等诸多优点,被广泛用于风电机组、直升机等重大机械装备,对保障设备的健康高效运转具有重大意义。然而其低速重载的恶劣工作环境导致行星轮、太阳轮及齿圈等关键部件易发生磨损、齿根裂纹甚至断齿等严重故障,致使设备运转振动加大以致运行系统整体瘫痪,产生严重后果。因此,针对行星齿轮箱的故障诊断尤其是早期微弱故障的诊断,具有重要的安全和经济价值。
随着信号处理理论的不断进步,针对行星齿轮箱的故障诊断技术已取得长足进展。然而因行星齿轮箱内部结构复杂、服役工况多变、信号传递路径长等多种因素,导致采集到的信号耦合性严重、信噪比较低,为故障的定位带来了巨大挑战。
中国专利(申请号为201810932320.9)为解决信号耦合性严重的问题,采用MED和模糊熵结合的方式进行行星齿轮箱的故障诊断,MED的基本原理是通过解反褶积突出尖脉冲,以峭度最大值作为迭代的终止条件,但直接用峭度作为迭代终止条件易受到随机冲击的影响,导致分辨不出冲击是由故障引起还是来源为随机冲击。中国专利(申请号为201711122369.X)先对信号进行奇异谱差分提取分量,随后进行EEMD分解,并计算得到的每个本征模态函数的分形维数,得到分形维数曲线,与正常信号与多种不同实际故障分形曲线对比来判断故障位置与类型;但在奇异谱分析中,得到的奇异值本身按照能量排序,如果仅挑选前N个奇异值,可能会将故障信息遗漏,从而造成误诊,且该方法最终和多种情况下的分形维数曲线对比,需要大量先验知识,限制了其应用范围。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供了一种基于低秩-稀疏逼近的行星齿轮箱故障诊断方法,分离故障信息和正常啮合信息,凸显故障冲击,实现行星齿轮箱的早期微弱故障诊断。
为达到上述目的,本发明采取的技术方案为:
基于低秩-稀疏逼近的行星齿轮箱故障诊断方法,包括以下步骤:
步骤一:针对行星齿轮箱,利用内置编码器采集输出轴角度位置信号x(t);
步骤二:以均方误差为损失函数,对输出轴角度位置信号x(t)进行局部多项式拟合,随后进行一阶差分,得到瞬时角速度信号w(t);
步骤三:基于步骤二得到的瞬时角速度信号w(t),分别对行星轮、太阳轮及齿圈时域同步平均,得到行星轮时域同步平均信号p(t)、太阳轮时域同步平均信号s(t)与齿圈时域同步平均信号r(t);
步骤四:针对步骤三得到的行星轮时域同步平均信号p(t)、太阳轮时域同步平均信号s(t)与齿圈时域同步平均信号r(t),分别用窗长为n的滑窗,以步长m构建信号矩阵,得到行星轮信号矩阵P、太阳轮信号矩阵S与齿圈信号矩阵R;
步骤五:依次对步骤四得到的行星轮信号矩阵P、太阳轮信号矩阵S与齿圈信号矩阵R进行低秩-稀疏逼近,得到行星轮稀疏矩阵P2、太阳轮稀疏矩阵S2、齿圈稀疏矩阵R2;
步骤六:对步骤五得到的行星轮稀疏矩阵P2、太阳轮稀疏矩阵S2、齿圈稀疏矩阵R2分别进行信号重构,得到行星轮稀疏矩阵恢复信号e1(t)、太阳轮稀疏矩阵恢复信号e2(t)、齿圈稀疏矩阵恢复信号e3(t);
步骤七:利用峭度指标,分别计算步骤六中行星轮稀疏矩阵恢复信号e1(t)、太阳轮稀疏矩阵恢复信号e2(t)、齿圈稀疏矩阵恢复信号e3(t)对应的行星轮峭度k1、太阳轮峭度k2、齿圈峭度k3,并设定峭度阈值α,大于峭度阈值α的则认为对应的齿轮存在故障。
所述的峭度阈值α为7。
本发明相比于现有技术,具有以下有益效果:
a)本发明在对信号进行时域同步平均的基础上,引入了低秩-稀疏逼近方法,对复杂的行星齿轮箱信号解耦,分离出复合信号中齿轮的啮合信息与存在的故障冲击成分,凸显故障,实现早期微弱故障的诊断。
b)不同于传统方法从信号时域、频域及时频域等多个维度对信号进行分析,本发明仅需对时域信号分析处理,可有效简化分析流程,减少计算资源。
c)本发明是一种自适应故障诊断方法,由数据本身即可实现诊断过程,无需训练数据与大量先验知识的支撑,有利于实现故障诊断的自动化。
附图说明
图1为本发明实施例实验台结构示意图。
图2为本发明实施例太阳轮故障示意图。
图3为本发明方法的流程图。
图4为本发明实施例的输出轴角度位置信号图。
图5为本发明实施例的瞬时角速度信号图。
图6为本发明实施例行星轮时域同步平均信号及稀疏恢复信号图。
图7为本发明实施例太阳轮时域同步平均信号及稀疏恢复信号图。
图8为本发明实施例齿圈时域同步平均信号及稀疏恢复信号图。
具体实施方式
下面结合附图和实施例对本发明做详细描述。
以某行星齿轮箱实验台为例进行说明,如图1所示,所用实验台由伺服电机、编码器、行星齿轮箱及磁粉制动器组成,其中行星齿轮箱由三个均布的行星轮、太阳轮和行星架组成,太阳轮与输入轴连接,通过啮合运动将动力通过行星架输出;行星齿轮箱的输入输出轴分别安装有编码器,由伺服电机驱动整个实验台。
实验台具体参数如下:
1)伺服电机额定功率:2.2kW,额定转速:600rpm;
2)行星齿轮箱传动比:5.1:1,内齿圈齿数:82,模数:1,行星轮齿数:31,模数:1,太阳轮齿数:20,模数:1;
针对太阳轮,人为引入故障,在太阳轮某个齿上粘贴尺寸为4mm*2mm*0.3mm的箔片,引入的故障太阳轮如图2所示。
如图3所示,基于低秩-稀疏逼近的行星齿轮箱故障诊断方法,包括以下步骤:
步骤一:针对行星齿轮箱,利用内置编码器采集输出轴角度位置信号x(t),如图4所示;
步骤二:以均方误差为损失函数,对输出轴角度位置信号x(t)进行局部多项式拟合,随后进行一阶差分,得到瞬时角速度信号w(t),如图5所示;
步骤三:基于步骤二得到的瞬时角速度信号w(t),分别对行星轮、太阳轮及齿圈进行时域同步平均,得到行星轮时域同步平均信号p(t)、太阳轮时域同步平均信号s(t)与齿圈时域同步平均信号r(t),时域同步平均的计算公式如下:
其中:q(t)表示求得的时域同步平均信号,w(t)表示瞬时角速度信号,P表示被平均的周期数,M表示每个周期内的采样点数;
步骤四:针对步骤三得到的行星轮时域同步平均信号p(t)、太阳轮时域同步平均信号s(t)与齿圈时域同步平均信号r(t),分别用窗长为n的滑窗,以步长m构建信号矩阵,得到行星轮信号矩阵P、太阳轮信号矩阵S与齿圈信号矩阵R;以行星轮时域同步平均信号p(t)为例,其构造过程如公式(2)所示:
其中ai为信号矩阵P中任一元素,i为大于等于1的正整数,n表示滑窗的窗长,m为移动步长;
步骤五:依次对步骤四得到的行星轮信号矩阵P、太阳轮信号矩阵S与齿圈信号矩阵R进行低秩-稀疏逼近,得到行星轮低秩矩阵P1、行星轮稀疏矩阵P2;太阳轮低秩矩阵S1、太阳轮稀疏矩阵S2;齿圈低秩矩阵R1、齿圈稀疏矩阵R2;
低秩-稀疏逼近方法同时将矩阵分解为低秩矩阵与稀疏矩阵两部分,如公式(3)所示:
其中:A表示低秩矩阵,即对应上述的行星轮低秩矩阵P1,太阳轮低秩矩阵S1与齿圈低秩矩阵R1,||A||*表示矩阵A的核范数;E表示稀疏矩阵,即对应上述的行星轮稀疏矩阵P2,太阳轮稀疏矩阵S2与齿圈稀疏矩阵R2,||E||l表示矩阵E的l1范数;D表示信号矩阵,即对应上述的行星轮信号矩阵P、太阳轮信号矩阵S与齿圈信号矩阵R;A+E=D表示约束条件;λ为常数,表示低秩项与稀疏项的平衡因子;
利用增广拉格朗日乘子法,将约束条件带入,公式(3)被转化为凸优化问题,其转化后形式如公式(4)所示:
其中:Y表示拉格朗日乘子组成的矢量;μ表示惩罚因子,初值大于0;||D-A-E||F表示矩阵的Frobenius范数;
之后利用交叉迭代法解公式(4)的凸优化问题,即得到低秩矩阵A与稀疏矩阵E;
步骤六:对步骤五得到的行星轮稀疏矩阵P2、太阳轮稀疏矩阵S2、齿圈稀疏矩阵R2分别进行信号重构恢复,得到行星轮稀疏矩阵恢复信号e1(t)、太阳轮稀疏矩阵恢复信号e2(t)、齿圈稀疏矩阵恢复信号e3(t);行星轮、太阳轮及齿圈的时域同步平均信号及各自的稀疏恢复信号如图6、图7、图8所示,其中行星轮及齿圈经时域同步平均后均无任何形式的冲击,它们对应的稀疏矩阵恢复信号也杂乱无章,并无周期性冲击,太阳轮经时域同步平均后,虽然能看到一定程度的冲击,但非常微弱,难以确定是否存在故障,而太阳轮的稀疏矩阵恢复信号存在非常明显的周期性冲击;
步骤七:利用峭度指标,分别计算步骤六中行星轮稀疏矩阵恢复信号e1(t)、太阳轮稀疏矩阵恢复信号e2(t)、齿圈稀疏矩阵恢复信号e3(t)对应的行星轮峭度k1、太阳轮峭度k2、齿圈峭度k3,并设定峭度阈值α,大于峭度阈值α的则认为对应的齿轮存在故障。
理论而言,服从高斯分布的数据,其峭度会在3附近浮动;但实际情况中,正常啮合信息和噪声的分布难以严格服从高斯分布,且低秩-稀疏逼近会变相放大冲击,导致其峭度值一般大于3,因此可适当将峭度阈值上调至7;大于峭度阈值7的则认为对应的齿轮存在故障;经计算,行星轮峭度k1为6.0014,太阳轮峭度k2为68.6953,齿圈峭度k3为6.1828,可知,太阳轮峭度k2远大于峭度阈值,而行星轮峭度k1与齿圈峭度k3则小于峭度阈值,因此可确定太阳轮发生故障。
Claims (2)
1.基于低秩-稀疏逼近的行星齿轮箱故障诊断方法,其特征在于,包括以下步骤:
步骤一:针对行星齿轮箱,利用内置编码器采集输出轴角度位置信号x(t);
步骤二:以均方误差为损失函数,对输出轴角度位置信号x(t)进行局部多项式拟合,随后进行一阶差分,得到瞬时角速度信号w(t);
步骤三:基于步骤二得到的瞬时角速度信号w(t),分别对行星轮、太阳轮及齿圈进行时域同步平均,得到行星轮时域同步平均信号p(t)、太阳轮时域同步平均信号s(t)与齿圈时域同步平均信号r(t);
步骤四:针对步骤三得到的行星轮时域同步平均信号p(t)、太阳轮时域同步平均信号s(t)与齿圈时域同步平均信号r(t),分别用窗长为n的滑窗,以步长m构建信号矩阵,得到行星轮信号矩阵P、太阳轮信号矩阵S与齿圈信号矩阵R;
步骤五:依次对步骤四得到的行星轮信号矩阵P、太阳轮信号矩阵S与齿圈信号矩阵R进行低秩-稀疏逼近,得到行星轮稀疏矩阵P2、太阳轮稀疏矩阵S2、齿圈稀疏矩阵R2;
步骤六:对步骤五得到的行星轮稀疏矩阵P2、太阳轮稀疏矩阵S2、齿圈稀疏矩阵R2分别进行信号重构,得到行星轮稀疏矩阵恢复信号e1(t)、太阳轮稀疏矩阵恢复信号e2(t)、齿圈稀疏矩阵恢复信号e3(t);
步骤七:利用峭度指标,分别计算步骤六中行星轮稀疏矩阵恢复信号e1(t)、太阳轮稀疏矩阵恢复信号e2(t)、齿圈稀疏矩阵恢复信号e3(t)对应的行星轮峭度k1、太阳轮峭度k2、齿圈峭度k3,并设定峭度阈值α,大于峭度阈值α的则认为对应的齿轮存在故障。
2.根据权利要求1所述的基于低秩-稀疏逼近的行星齿轮箱故障诊断方法,其特征在于:所述的峭度阈值α为7。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910270414.9A CN109883693B (zh) | 2019-04-04 | 2019-04-04 | 基于低秩-稀疏逼近的行星齿轮箱故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910270414.9A CN109883693B (zh) | 2019-04-04 | 2019-04-04 | 基于低秩-稀疏逼近的行星齿轮箱故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109883693A CN109883693A (zh) | 2019-06-14 |
CN109883693B true CN109883693B (zh) | 2020-01-21 |
Family
ID=66936121
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910270414.9A Active CN109883693B (zh) | 2019-04-04 | 2019-04-04 | 基于低秩-稀疏逼近的行星齿轮箱故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109883693B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110285970B (zh) * | 2019-07-18 | 2020-05-12 | 东北大学 | 基于矩阵恢复的滚动轴承微弱故障增强方法 |
CN110580471B (zh) * | 2019-09-12 | 2021-11-02 | 北京航空航天大学 | 基于编码器信号瞬变特征的机械设备故障诊断方法 |
CN111535999B (zh) * | 2020-05-22 | 2021-08-24 | 三一重能有限公司 | 风机脱落物监测方法、装置、系统和存储介质 |
CN111855208B (zh) * | 2020-07-02 | 2021-09-21 | 南京航空航天大学 | 一种基于阶次分析和稀疏编码的旋转机械微弱故障信号提取方法 |
CN113465916B (zh) * | 2021-06-24 | 2022-06-07 | 西安交通大学 | 一种行星轮系轮齿状态评估方法、装置、设备和介质 |
CN114942133B (zh) * | 2022-05-20 | 2023-04-14 | 大连理工大学 | 基于最优秩非负矩阵分解的行星齿轮箱早期故障诊断方法 |
CN115575488B (zh) * | 2022-11-21 | 2023-06-02 | 中国石油大学(华东) | 基于acfm与视觉的集成探头及协同检测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102636347A (zh) * | 2012-04-24 | 2012-08-15 | 西安交通大学 | 一种用于变转速齿轮箱的振动信号时域同步平均方法 |
CN103780520A (zh) * | 2014-01-07 | 2014-05-07 | 清华大学 | 基于先验信息辅助的压缩感知窄带干扰估计方法及装置 |
CN105716857A (zh) * | 2016-01-20 | 2016-06-29 | 中国人民解放军军械工程学院 | 一种行星齿轮箱健康状态评估方法 |
EP3284671A1 (en) * | 2016-08-17 | 2018-02-21 | Bell Helicopter Textron Inc. | Diagnostic method, system and device for a rotorcraft drive system |
CN108362492A (zh) * | 2018-01-11 | 2018-08-03 | 中国人民解放军国防科技大学 | 一种适用于低转速下行星轮系故障诊断的振动分离方法 |
-
2019
- 2019-04-04 CN CN201910270414.9A patent/CN109883693B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102636347A (zh) * | 2012-04-24 | 2012-08-15 | 西安交通大学 | 一种用于变转速齿轮箱的振动信号时域同步平均方法 |
CN103780520A (zh) * | 2014-01-07 | 2014-05-07 | 清华大学 | 基于先验信息辅助的压缩感知窄带干扰估计方法及装置 |
CN105716857A (zh) * | 2016-01-20 | 2016-06-29 | 中国人民解放军军械工程学院 | 一种行星齿轮箱健康状态评估方法 |
EP3284671A1 (en) * | 2016-08-17 | 2018-02-21 | Bell Helicopter Textron Inc. | Diagnostic method, system and device for a rotorcraft drive system |
CN108362492A (zh) * | 2018-01-11 | 2018-08-03 | 中国人民解放军国防科技大学 | 一种适用于低转速下行星轮系故障诊断的振动分离方法 |
Non-Patent Citations (1)
Title |
---|
基于时域同步平均法的行星齿轮振动信号分离技术;胡贵锋等;《中国机械工程》;20130331(第6期);第787-791页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109883693A (zh) | 2019-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109883693B (zh) | 基于低秩-稀疏逼近的行星齿轮箱故障诊断方法 | |
CN109029977B (zh) | 一种基于vmd-amckd的行星齿轮箱早期故障诊断方法 | |
Wang et al. | Unknown fault feature extraction of rolling bearings under variable speed conditions based on statistical complexity measures | |
Pan et al. | Nonlinear sparse mode decomposition and its application in planetary gearbox fault diagnosis | |
CN102879196B (zh) | 利用矩阵小波变换的行星齿轮箱复合故障诊断方法 | |
CN108844733B (zh) | 一种基于kl散度与均方根值的齿轮状态监测指标提取方法 | |
CN110132578B (zh) | 一种齿轮系统复合故障特征提取方法及故障试验装置 | |
CN105510023A (zh) | 基于散度指标的变工况风电行星齿轮箱故障诊断方法 | |
CN109060347B (zh) | 基于堆叠消噪自动编码器和门控循环单元神经网络的行星齿轮故障识别方法 | |
CN109682597A (zh) | 一种齿轮箱振动信号处理和分析方法 | |
CN105445022A (zh) | 一种基于双树复小波变换-熵特征融合的行星齿轮故障诊断方法 | |
CN105699072A (zh) | 一种基于级联经验模态分解齿轮故障诊断方法 | |
CN112577736A (zh) | 基于sanc和1d-cnn-lstm的风电机组行星齿轮箱故障诊断方法 | |
CN115901248B (zh) | 一种齿轮箱故障特征提取方法 | |
CN111413095A (zh) | 基于瞬时角速度的行星轴承分布式故障诊断分析方法 | |
CN109883692B (zh) | 基于内置编码器信息的广义差分滤波方法 | |
CN115434872A (zh) | 一种基于avmd与改进rssd的风电机组齿轮箱复合故障诊断方法 | |
CN107966287B (zh) | 一种自适应机电装备微弱故障特征提取方法 | |
Wang et al. | Research on gearbox composite fault diagnosis based on improved local mean decomposition | |
Xiangjun et al. | Anomaly detection of wind turbine gearbox based on digital twin drive | |
CN112364551B (zh) | 一种行星传动齿轮箱疲劳分析方法 | |
CN117370781A (zh) | 一种rv减速器复合缺陷特征提取方法、装置 | |
Wang et al. | Fusion deconvolution for reliability analysis of a flywheel-battery hybrid energy storage system | |
Lu et al. | Research on Transmission Efficiency Prediction of Heavy-Duty Tractors HMCVT Based on VMD and PSO–BP | |
CN114840920B (zh) | 一种混合动力乘用车机电耦合装置载荷谱编制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |