CN104215453B - 一种一级行星二级平行轴齿轮箱故障检测方法 - Google Patents

一种一级行星二级平行轴齿轮箱故障检测方法 Download PDF

Info

Publication number
CN104215453B
CN104215453B CN201410468968.7A CN201410468968A CN104215453B CN 104215453 B CN104215453 B CN 104215453B CN 201410468968 A CN201410468968 A CN 201410468968A CN 104215453 B CN104215453 B CN 104215453B
Authority
CN
China
Prior art keywords
frequency
gear
vibration
signal
modulation
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.)
Expired - Fee Related
Application number
CN201410468968.7A
Other languages
English (en)
Other versions
CN104215453A (zh
Inventor
万书亭
豆龙江
胡媛媛
绳晓玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wenzhou Huadian Industrial Co ltd
Original Assignee
North China Electric Power University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN201410468968.7A priority Critical patent/CN104215453B/zh
Publication of CN104215453A publication Critical patent/CN104215453A/zh
Application granted granted Critical
Publication of CN104215453B publication Critical patent/CN104215453B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

一种一级行星二级平行轴齿轮箱故障检测方法,所述方法首先利用振动加速度传感器采集齿轮箱输出轴端垂向振动加速度信号和水平振动加速度信号,并对采集的信号进行功率谱分析,然后将功率谱中的不同频率分量与齿轮箱正常振动数据样本中的对应频率分量相比较,从而判断风力发电机行星齿轮箱内部行星轮系的故障。本发明只需要在行星齿轮箱的输出轴端安装振动加速度传感器,而且只需对采集的振动信号进行功率谱分析就可以快速、准确地诊断出故障所在的位置。具有信号采集装置简单,操作方便,诊断结果准确等优点。

Description

一种一级行星二级平行轴齿轮箱故障检测方法
技术领域
本发明涉及一种仅利用齿轮箱输出轴端振动信号诊断一级行星二级平行轴齿轮箱故障的检测方法,属于风电技术领域。
背景技术
行星齿轮箱是风力发电机组中的重要传动部件之一,主要由太阳轮7、行星轮8、内齿圈6、行星架13和与太阳轮相接的齿轮增速机构组成,一级行星二级平行轴齿轮箱的齿轮增速机构包括与太阳轮轴2同轴连接的太阳轮轴输出齿轮9、与太阳轮轴输出齿轮9啮合的第一级平行轴输入齿轮10、通过第一级平行轴3与第一级平行轴输入齿轮10同轴固接的第一级平行轴输出齿轮11、固定在输出轴4上且与第一级平行轴输出齿轮11啮合的输出轴齿轮12(参看图1)。通常情况下,内齿圈固定不动,太阳轮绕自身的中心轴线旋转,行星轮不仅自转,还围绕太阳轮公转。行星轮既和太阳轮啮合,又和内齿圈啮合。行星齿轮箱中多个齿轮的复杂运动造成了振动信号成分的复杂性。其中包含太阳轮、行星轮的旋转频率;太阳轮和行星轮、行星轮和内齿圈之间的啮合频率;上述频率的倍频成分;成组出现在啮合频率及其倍频附近的边带成分。在进行振动监测时,传感器通常被安装在内齿圈及与之相连的箱体上,太阳轮和行星轮以及行星轮和内齿圈啮合副的啮合点相对于传感器的位置随行星架旋转而发生变化,使得啮合点至传感器之间的振动传递路径发生变化,这种时变的传递路径对振动信号产生幅值调制效应,进一步增加了信号的复杂性。
在风力发电机组运行过程中,行星齿轮箱不仅承受动态重载负荷,而且运行工况变化频繁,是整个动力传动链中的薄弱环节。其中太阳轮、行星轮、内齿圈等关键零部件容易出现损伤故障。一旦出现故障,不仅维修困难,而且维修费用也很高。因此,对其开展监测诊断对于保障风力发电机组安全高效稳定运行至关重要。但目前提出的一些风力发电机组故障诊断方法主要针对发电机、风力机叶片以及传动系统中的轴承和定轴轮系等,对于行星齿轮箱的研究则比较少见。而且齿轮箱信号采集时要使用大量的传感器,给采集工作带来了很大的困难。因此,如何简化采集装置以及快速判断齿轮箱行星轮系故障位置是有关技术人员研究的重点。
发明内容
本发明的目的在于针对现有技术之弊端,提供一种简单易行的一级行星二级平行轴齿轮箱故障检测方法,以正确诊断一级行星二级平行轴齿轮箱行星轮系的故障。
本发明所述问题是以下述技术方案实现的:
一种一级行星二级平行轴齿轮箱故障检测方法,所述方法首先利用振动加速度传感器采集齿轮箱输出轴端垂向振动加速度信号和水平振动加速度信号,并对采集的信号进行功率谱分析,然后将功率谱中的不同频率分量与齿轮箱正常振动数据样本中的对应频率分量相比较,从而判断风力发电机行星齿轮箱内部行星轮系的故障,所述方法的具体步骤如下:
①采用振动加速度传感器测取一级行星二级平行轴齿轮箱输出轴端振动信号;
②对振动信号进行功率谱分析;
③将功率谱信号中的不同频率分量与一级行星二级平行轴齿轮箱正常振动数据样本中的对应频率分量相比较,判断一级行星二级平行轴齿轮箱行星轮系的故障:
a.功率谱中若含有KZ1f1±mfr±f2±f3±f4的频率成分,则表明行星轮系的内齿圈发生了故障;
b.功率谱中若含有KZ1f1±mfp±f2±f3±f4的频率成分,则表明行星轮系的行星轮发生了故障;
c.功率谱中若含有KZ1f1±mfs±f2±f3±f4的频率成分,则表明行星轮系的太阳轮发生了故障;
其中,Z1为内齿圈的齿数,f1为输入轴转动频率实时监测值,f2为太阳轮轴转动频率实时监测值,f3为第一级平行轴转动频率实时监测值,f4为输出轴转动频率实时监测值,fr为内齿圈故障特征频率,fp为行星轮故障特征频率,fs为太阳轮故障特征频率,K、m为在(0,∞)范围内变化的整数。
上述一级行星二级平行轴齿轮箱故障检测方法,所述故障特征频率按如下方式获得:
内齿圈故障特征频率由下式获得:
行星轮故障特征频率由下式获得:
太阳轮故障特征频率由下式获得:
其中,N为行星轮的数量,Z2为行星轮的齿数,Z3为太阳轮的齿数。
本发明只需要在行星齿轮箱的输出轴端安装振动加速度传感器,而且只需对采集的振动信号进行功率谱分析就可以快速、准确地诊断出故障所在的位置。具有信号采集装置简单,操作方便,诊断结果准确等优点。通过对比试验和测定,其判断准确率达到93%以上。
附图说明
下面结合附图对本发明作进一步说明。
图1是一级行星二级平行轴齿轮箱传动结构图;
图2是实施方案图;
图3是故障判别流程图。
图中各标号清单为:1、输入轴;2、太阳轮轴;3、第一级平行轴;4、输出轴;5、振动加速度传感器;6、内齿圈;7、太阳轮;8、行星轮;9、太阳轮轴输出齿轮;10、第一级平行轴输入齿轮;11、第一级平行轴输出齿轮;12、输出轴齿轮;13、行星架。
文中各符号清单为:Z1、内齿圈的齿数;Z2、行星轮齿数;Z3、太阳轮齿数;f1、输入轴转动频率实时监测值;f2、太阳轮轴转动频率实时监测值;f3、第一级平行轴转动频率实时监测值;f4、输出轴转动频率实时监测值;fr、内齿圈故障特征频率;fp、行星轮故障特征频率;fs、太阳轮故障特征频率;K、m、在(0,∞)范围内变化的整数;N、行星轮数量。
具体实施方式
本发明是针对现有故障信息采集方法较为复杂,风力发电机组一级行星二级平行轴齿轮箱行星轮系故障诊断较为复杂这两个问题设计的。
本发明需要解决两个问题:一是如何从振动信号中分离出有效的频率成分;二是如何判断行星轮系故障的位置。
参看图1,针对本发明提出的问题一,测量一级行星二级平行轴齿轮箱输出轴端的振动信号。通过以下方法实现:在一级行星二级平行轴齿轮箱输出轴端安装振动加速度传感器5,振动加速度在输出轴端的安装分水平和垂直两个方向,通过不同方向采集齿轮箱输出轴端的振动信号,后续处理中不同方向的振动信号分类处理。
根据振动调制机理,那么太阳轴输出轴振动啮合频率调制模型为:
其中,A为信号幅值,f1为输入轴转动频率、f2为太阳轴的转动频率,β为调频的调制系数,B为调频的调制系数,
振动以上述调制形式作用于太阳轮2,经与太阳轮2连接的第一级平行轴3传递到太阳轮轴输出齿轮9,并间接作用在第一级平行轴输入齿轮10上,由于第一级平行轴3的转动频率为f3,在三种振动的共同作用下,此处忽略传递过程中振幅的变化仅考虑频率特征,将会出现以f3为调制信号,振动啮合频率调制模型为:
其中,f3为第一级平行轴3的转动频率,振动以上述调制形式经第一级平行轴3传递给第一级平行轴输出齿轮11,并间接作用在输出轴齿轮12上,由于输出轴4的转动频率为f4,在上述各种调制信号振动的共同作用下,此处忽略传递过程中振幅的变化仅考虑频率特征,将会出现以f4为调制信号,此时振动啮合频率调制模型为:
其中,f4为输出轴4的转动频率,上述振动模型经傅里叶变换得到:
其中,Jm(β)为变量β的第一类Bessel函数。
由上式可知,信号分别含有Z1f1、Z1f1±f2、Z1f1±f3、Z1f1±f4、Z1f1±f2±f3、Z1f1±f2±f4、Z1f1±f3±f4、Z1f1±f2±f3±f4的一组频率成分,对于齿轮的固有频率调制现象以及齿轮固箱有频率调制现象也有类似的结论。不过此处的Z1f1要替换为齿轮固有频率fg和齿轮箱固有频率fx
由上述啮合频率调制的推导类比出齿轮固有频率调制后,信号含有fg、fg±f2、fg±f3、fg±f4、fg±f2±f3、fg±f2±f4、fg±f3±f4、fg±f2±f3±f4的一组频率成分;箱体固有频率的调制后,信号含有fx、fx±f2、fx±f3、fx±f4、fx±f2±f3、fx±f2±f4、fx±f3±f4、fx±f2±f3±f4的一组频率成分。
考虑到太阳轮、内齿圈和行星轮同时存在的情况,以fs、fr和fp分别代表太阳轮、内齿圈和行星轮的故障特征频率,信号经调制后会出现Z1f1±K1fs±K2fr±K3fp、fg±K1fs±K2fr±K3fp和fx±K1fs±K2fr±K3fp的振动分量,此时可以诊断出故障的具体位置,其中K1、K2、K3为在(-∞,∞)范围内变化的整数。
综上,若振动信号功率谱中含有KZ1f1±mfr±f2±f3±f4的频率成分,则表明是行星轮系的内齿圈发生了故障;含有KZ1f1±mfp±f2±f3±f4的频率成分,则表明是行星轮系的行星轮发生了故障;含有KZ1f1±mfs±f2±f3±f4的频率成分,则表明是行星轮系的太阳轮发生了故障。

Claims (2)

1.一种一级行星二级平行轴齿轮箱故障检测方法,其特征是,所述方法首先利用振动加速度传感器采集齿轮箱输出轴端垂向振动加速度信号和水平振动加速度信号,并对采集的信号进行功率谱分析,然后将功率谱中的不同频率分量与齿轮箱正常振动数据样本中的对应频率分量相比较,从而判断风力发电机行星齿轮箱内部行星轮系的故障;
具体步骤如下:
①采用振动加速度传感器测取一级行星二级平行轴齿轮箱输出轴端振动信号;
②对振动信号进行功率谱分析;
③将功率谱信号中的不同频率分量与一级行星二级平行轴齿轮箱正常振动数据样本中的对应频率分量相比较,判断一级行星二级平行轴齿轮箱行星轮系的故障:
a.功率谱中若含有KZ1f1±mfr±f2±f3±f4的频率成分,则表明行星轮系的内齿圈发生了故障;
b.功率谱中若含有KZ1f1±mfp±f2±f3±f4的频率成分,则表明行星轮系的行星轮发生了故障;
c、功率谱中若含有KZ1f1±mfs±f2±f3±f4的频率成分,则表明行星轮系的太阳轮发生了故障;
其中,Z1为内齿圈的齿数,f1为输入轴转动频率实时监测值,f2为太阳轮轴转动频率实时监测值,f3为第一级平行轴转动频率实时监测值,f4为输出轴转动频率实时监测值,fr为内齿圈故障特征频率,fp为行星轮故障特征频率,fs为太阳轮故障特征频率,K、m为在(0,∞)范围内变化的整数;
所述步骤①中,测取振动信号的具体方法为:在一级行星二级平行轴齿轮箱输出轴端安装振动加速度传感器(5),振动加速度在输出轴端的安装分水平和垂直两个方向,通过不同方向采集齿轮箱输出轴端的振动信号,后续处理中不同方向的振动信号分类处理;
所述步骤②中,对振动信号进行功率谱分析的具体方法为:
根据振动调制机理,建立太阳轴输出轴振动啮合频率调制模型为:
其中,A为信号幅值,f1为输入轴转动频率、f2为太阳轴的转动频率,β为调频的调制系数,B为调频的调制系数;
振动以上述调制形式作用于太阳轮2,经与太阳轮2连接的第一级平行轴3传递到太阳轮轴输出齿轮9,并间接作用在第一级平行轴输入齿轮10上,由于第一级平行轴3的转动频率为f3,在三种振动的共同作用下,此处忽略传递过程中振幅的变化仅考虑频率特征,将会出现以f3为调制信号,振动啮合频率调制模型为:
其中,f3为第一级平行轴3的转动频率,振动以上述调制形式经第一级平行轴3传递给第一级平行轴输出齿轮11,并间接作用在输出轴齿轮12上,由于输出轴4的转动频率为f4,在上述各种调制信号振动的共同作用下,此处忽略传递过程中振幅的变化仅考虑频率特征,将会出现以f4为调制信号,此时振动啮合频率调制模型为:
其中,f4为输出轴4的转动频率,上述振动模型经傅里叶变换得到:
其中,Jm(β)为变量β的第一类Bessel函数;
由上式可知,信号分别含有Z1f1、Z1f1±f2、Z1f1±f3、Z1f1±f4、Z1f1±f2±f3、Z1f1±f2±f4、Z1f1±f3±f4、Z1f1±f2±f3±f4的一组频率成分,对于齿轮的固有频率调制现象以及齿轮固箱有频率调制现象也有类似的结论,不过此处的Z1f1要替换为齿轮固有频率fg和齿轮箱固有频率fx
由上述啮合频率调制的推导类比出齿轮固有频率调制后,信号含有fg、fg±f2、fg±f3、fg±f4、fg±f2±f3、fg±f2±f4、fg±f3±f4、fg±f2±f3±f4的一组频率成分;箱体固有频率的调制后,信号含有fx、fx±f2、fx±f3、fx±f4、fx±f2±f3、fx±f2±f4、fx±f3±f4、fx±f2±f3±f4的一组频率成分;
考虑到太阳轮、内齿圈和行星轮同时存在的情况,以fs、fr和fp分别代表太阳轮、内齿圈和行星轮的故障特征频率,信号经调制后会出现Z1f1±K1fs±K2fr±K3fp、fg±K1fs±K2fr±K3fp和fx±K1fs±K2fr±K3fp的振动分量,此时可以诊断出故障的具体位置,其中K1、K2、K3为在(-∞,∞)范围内变化的整数。
2.根据权利要求1所述的一种一级行星二级平行轴齿轮箱故障检测方法,其特征是,所述故障特征频率按如下方式获得:
内齿圈故障特征频率由下式获得:
行星轮故障特征频率由下式获得:
太阳轮故障特征频率由下式获得:
其中,N为行星轮的数量,Z2为行星轮的齿数,Z3为太阳轮的齿数。
CN201410468968.7A 2014-09-16 2014-09-16 一种一级行星二级平行轴齿轮箱故障检测方法 Expired - Fee Related CN104215453B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410468968.7A CN104215453B (zh) 2014-09-16 2014-09-16 一种一级行星二级平行轴齿轮箱故障检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410468968.7A CN104215453B (zh) 2014-09-16 2014-09-16 一种一级行星二级平行轴齿轮箱故障检测方法

Publications (2)

Publication Number Publication Date
CN104215453A CN104215453A (zh) 2014-12-17
CN104215453B true CN104215453B (zh) 2016-07-06

Family

ID=52097176

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410468968.7A Expired - Fee Related CN104215453B (zh) 2014-09-16 2014-09-16 一种一级行星二级平行轴齿轮箱故障检测方法

Country Status (1)

Country Link
CN (1) CN104215453B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110455528A (zh) * 2018-05-03 2019-11-15 国家能源投资集团有限责任公司 用于行星齿轮损伤诊断的方法和系统

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109073505A (zh) * 2016-03-30 2018-12-21 Ntn株式会社 齿轮装置的状态监视系统以及状态监视方法
CN107607321B (zh) * 2017-09-06 2019-11-05 成都大汇物联科技有限公司 一种设备故障精确定位方法
CN111259737B (zh) * 2020-01-08 2023-07-25 科大讯飞股份有限公司 车辆方向盘故障的预测方法、装置、电子设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201016849Y (zh) * 2007-03-14 2008-02-06 中国南方工业汽车股份有限公司四川建安车桥分公司 齿轮啮合振动测试装置
CN102288286A (zh) * 2011-06-16 2011-12-21 中国科学院沈阳自动化研究所 一种振动加速度传感器齿轮箱测点精度分析与评价方法
CN102539151A (zh) * 2012-01-18 2012-07-04 北京工业大学 一种汽车变速器智能化在线质量检测方法
CN103234750A (zh) * 2013-04-16 2013-08-07 北京工业大学 一种基于改进倒频谱法的等高齿锥齿轮故障诊断方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI447371B (zh) * 2010-11-18 2014-08-01 Univ Nat Taiwan Science Tech 即時檢測系統和方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201016849Y (zh) * 2007-03-14 2008-02-06 中国南方工业汽车股份有限公司四川建安车桥分公司 齿轮啮合振动测试装置
CN102288286A (zh) * 2011-06-16 2011-12-21 中国科学院沈阳自动化研究所 一种振动加速度传感器齿轮箱测点精度分析与评价方法
CN102539151A (zh) * 2012-01-18 2012-07-04 北京工业大学 一种汽车变速器智能化在线质量检测方法
CN103234750A (zh) * 2013-04-16 2013-08-07 北京工业大学 一种基于改进倒频谱法的等高齿锥齿轮故障诊断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于多传感器信息融合的行星齿轮箱故障诊断;雷亚国;《振动与冲击》;20101231;第29卷(第S2010期);第216-218页 *
机械系统故障信号提取技术研究;余红英;《中国优秀博士学位论文全文数据库,工程科技II辑》;20060831(第8期);第C029-14页,第2章 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110455528A (zh) * 2018-05-03 2019-11-15 国家能源投资集团有限责任公司 用于行星齿轮损伤诊断的方法和系统

Also Published As

Publication number Publication date
CN104215453A (zh) 2014-12-17

Similar Documents

Publication Publication Date Title
Teng et al. Multi-fault detection and failure analysis of wind turbine gearbox using complex wavelet transform
Shahriar et al. Electrical signature analysis-based detection of external bearing faults in electromechanical drivetrains
Bajric et al. Feature extraction using discrete wavelet transform for gear fault diagnosis of wind turbine gearbox
Assaad et al. Vibration based condition monitoring of a multistage epicyclic gearbox in lifting cranes
Yoon et al. On the use of a single piezoelectric strain sensor for wind turbine planetary gearbox fault diagnosis
CN105806613A (zh) 一种基于阶比复杂度的行星齿轮箱故障诊断方法
CN104215453B (zh) 一种一级行星二级平行轴齿轮箱故障检测方法
CN110044610B (zh) 齿轮故障诊断方法
Wang et al. Fault diagnosis for wind turbine planetary ring gear via a meshing resonance based filtering algorithm
CN105510023A (zh) 基于散度指标的变工况风电行星齿轮箱故障诊断方法
CN103857995A (zh) 用于变速箱的状态监控的方法和系统
Lu et al. Condition monitoring and fault diagnostics of wind turbines
CN106055734A (zh) 一种风电齿轮箱故障诊断方法
CN104568435A (zh) 一种基于时频能量谱提取风机状态信号调制成分的方法
CN109253244A (zh) 一种多齿轮系统大型机械齿轮箱故障检测方法
CN105738102A (zh) 一种风电齿轮箱故障诊断方法
Lin et al. A review and strategy for the diagnosis of speed-varying machinery
Hanna et al. Detection of wind turbine gear tooth defects using sideband energy ratio™
CN105548595A (zh) 一种提取风电齿轮箱各级轴转速检测方法
Yoon et al. Planetary gearbox fault diagnosis using a single piezoelectric strain sensor
CN103940607B (zh) 不依赖于时域平均的行星齿轮箱信号分离和诊断方法
Hong et al. Application of spectral kurtosis and Vold-Kalman filter based order tracking in wind turbine gearbox fault diagnosis
Wang et al. Tacholess order-tracking approach for wind turbine gearbox fault detection
Dziedziech et al. Epicyclic gearbox fault detection by instantaneous circular pitch cycle map
Gryllias et al. Advanced cyclostationary-based analysis for condition monitoring of complex systems

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20220506

Address after: 325200 No. 266, YuMeng Road, Dongshan Economic and Technological Development Zone, Rui'an City, Wenzhou City, Zhejiang Province

Patentee after: Wenzhou Huadian Industrial Co.,Ltd.

Address before: 071003 Hebei province Baoding Yonghua No. 619 North Street

Patentee before: NORTH CHINA ELECTRIC POWER University (BAODING)

TR01 Transfer of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160706

CF01 Termination of patent right due to non-payment of annual fee