CN104215453B - 一种一级行星二级平行轴齿轮箱故障检测方法 - Google Patents
一种一级行星二级平行轴齿轮箱故障检测方法 Download PDFInfo
- 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
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为太阳轮的齿数。
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110455528A (zh) * | 2018-05-03 | 2019-11-15 | 国家能源投资集团有限责任公司 | 用于行星齿轮损伤诊断的方法和系统 |
Families Citing this family (3)
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)
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI447371B (zh) * | 2010-11-18 | 2014-08-01 | Univ Nat Taiwan Science Tech | 即時檢測系統和方法 |
-
2014
- 2014-09-16 CN CN201410468968.7A patent/CN104215453B/zh not_active Expired - Fee Related
Patent Citations (4)
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)
Title |
---|
基于多传感器信息融合的行星齿轮箱故障诊断;雷亚国;《振动与冲击》;20101231;第29卷(第S2010期);第216-218页 * |
机械系统故障信号提取技术研究;余红英;《中国优秀博士学位论文全文数据库,工程科技II辑》;20060831(第8期);第C029-14页,第2章 * |
Cited By (1)
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 |