CN107907327A - 一种风电机组行星齿轮箱故障诊断方法 - Google Patents
一种风电机组行星齿轮箱故障诊断方法 Download PDFInfo
- Publication number
- CN107907327A CN107907327A CN201711122369.XA CN201711122369A CN107907327A CN 107907327 A CN107907327 A CN 107907327A CN 201711122369 A CN201711122369 A CN 201711122369A CN 107907327 A CN107907327 A CN 107907327A
- Authority
- CN
- China
- Prior art keywords
- fault
- signal
- curve
- fractal dimension
- wind turbine
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/021—Gearings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
Abstract
本发明涉及一种风电机组行星齿轮箱故障诊断方法,包括以下步骤:1、利用奇异值差分谱法对采集的振动信号中的高频冲击信号和噪声信号的复合分量进行提取;S2、对采集的振动信号添加白噪声进行N次EEMD分解,得到一系列本征模态函数,其中,白噪声的添加准则为:3、计算步骤2得到的本征模态函数的分形维数,得到对应的分形维数曲线;4、将分形维数曲线与正常信号和不同实际故障曲线对比,根据曲线相似度判定振动信号所属的故障位置和故障类型。与现有技术相比,本发明使EEMD方法分解过程可以自适应设定白噪声大小和分解次数,解决了传统EEMD模态混叠和参数设定的盲目性问题,具有分解效果好、分解效率高等优点。
Description
技术领域
本发明涉及风电机组故障检测领域,尤其是涉及一种风电机组行星齿轮箱故障诊断方法。
背景技术
行星齿轮箱是风电机组中重要部件。由于其本身结构复杂,且经常处在复杂多变的运行工况下。导致风电机组的行星齿轮箱是风电机组中故障率最高的部件之一,且其故障引起停机时间最长。所以对于该部件实现准确、迅速的故障定位具有极大研究意义。
在众多故障检测方法中,EMD方法作为一种全新的、创造性的信号时频处理方法,对于非线性、非平稳信号的分析的效果显著,所以在风电故障检测中使用率极高,并且取得了较好的实际应用效果。但当信号的极值点分布不均匀时,EMD分解易导致“过冲”、“欠冲”现象,进而造成模态混叠现象。
为了解决模态混叠问题针对上述缺陷,研究人员提出了EEMD方法,从解决极值点分布问题出发,通过在初始信号中添加高斯白噪声来使极值点分布均匀化,以达到减小模态混叠影响提高分解效果的目地。该方法在一定程度上改善了模态混叠现象,但添加噪声所需两个重要参数即白噪声标准差和总体平均次数均需人为经验选择,由于掺杂了人为经验干预,大大影响了分解效果以及方法的自适应性。另一方面,分形维数作为衡量信号分形特征的手段,与信号的频率和振幅关系密切,该方法在传统的使用中均用于故障严重程度的分析中,在故障定位中鲜有应用。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种风电机组行星齿轮箱故障诊断方法。
本发明的目的可以通过以下技术方案来实现:
一种风电机组行星齿轮箱故障诊断方法,包括以下步骤:
S1、利用奇异值差分谱法对采集的振动信号中的高频冲击信号和噪声信号的复合分量进行提取;
S2、对采集的振动信号添加白噪声进行N次EEMD分解,得到一系列本征模态函数,其中,白噪声的添加准则为:
e表示添加白噪声的标准差,表示经过奇异值差分谱法提取出的高频冲击分量和噪声分量信号复合分量的标准差;
S3、计算步骤S2得到的本征模态函数的分形维数,得到对应的分形维数曲线;
S4、将分形维数曲线与正常信号和不同故障实际曲线对比,根据曲线相似度判定振动信号所属的故障位置和故障类型。
优选的,所述白噪声e的添加准则为:
优选的,所述EEMD分解的次数N为:
其中,ε表示允许分解误差。
优选的,所述允许分解误差ε为2%。
优选的,所述分形维数的计算方法为通过最小二乘法计算信号的网格维数。
优选的,所述步骤S4具体包括:
S41、将第一层分形维数曲线与正常信号和表示不同故障位置的实际曲线对比,若与第一层分形维数曲线最相似的实际曲线为故障曲线,则将其对应的故障位置作为振动信号所属的故障位置,进行步骤S42,否则判定振动信号为正常信号,结束判定;
S42、将第二层分形维数曲线与已确定故障位置的不同故障类型的实际曲线对比,将最相似的实际曲线对应的故障类型作为振动信号所属的故障类型。
优选的,所述故障位置包括:行星轮和太阳轮。
优选的,所述故障类型包括:断齿、点蚀、磨损和缺齿。
与现有技术相比,本发明具有以下优点:
1、利用奇异值差分谱法对信号进行分解重构,增加了预先滤波来确定分解参数的预处理步骤,使EEMD方法分解过程可以自适应设定白噪声大小和分解次数,解决了传统EEMD模态混叠和参数设定的盲目性问题,改善了分解效果,提高了分解效率。
2、使用分层的分形维数实现先判定故障位置、再判定故障类型,提高了故障诊断的精确度。
附图说明
图1为本发明方法流程图;
图2为本发明中分层分形维数判定的流程示意图;
图3为实施例中仿真信号示意图;
图4为实施例中采用传统EEMD分解结果图;
图5为实施例中采用本发明提出的EEMD分解结果图;
图6为实施例中实际信号与一层分形维数曲线对比图;
图7为实施例中实际信号与二层分形维数曲线对比图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例
在EMD分解过程中高频冲击信号是造成模态混叠现象的最主要因素,噪声信号则是次要因素。传统EEMD添加白噪声辅助分解方法的目的在于均匀信号极值点分布,但在高频冲击信号和噪声信号多变的实际工况下,传统EEMD方法无法自适应地调整分解参数值,这将直接影响分解效果,进而影响对于运行工况的诊断和分析。
本申请提出的自适应EEMD分解中添加白噪声标准差的准则如下:
其中,表示经过奇异值差分谱法提取出的高频冲击分量和噪声分量信号复合分量的标准差,为原始信号的标准差,e表示添加白噪声的标准差。
因为添加的是白噪声,所以e>0,且应同时保证e的上限为不改变高频冲击信号极值点分布的白噪声n(t)的标准差上限。由置信区间概念可知,对于参数θ,如果统计量分别为满足给定的T∈(0,1)则可以得到则分别为θ的置信下限和置信上限,1-T为θ的置信水平
因此,根据正态分布函数的概率分布P(-3e<x<+3e)=99.73%。如果取则添加的白噪声n(t)幅值分布在-3e<n(t)<3e区间内的概率为99.73%,即离散白噪声的幅值的绝对值小于高频冲击信号标准差的概率为99.73%,因此可以保证添加的白噪声不会影响高频冲击信号的极值点分布。由上述过程可知添加白噪声的准则为:
当上式取等号时可以取得良好的分解效果。
由上式可得EEMD分解的次数N为:
其中,ε表示允许分解误差,一般设定为2%。
如图1所示,一种风电机组行星齿轮箱故障诊断方法,包括以下步骤:
S1、利用奇异值差分谱法对采集的振动信号中的高频冲击信号和噪声信号的复合分量进行提取;
S2、对采集的振动信号添加白噪声进行N次EEMD分解,得到一系列本征模态函数(IMF),其中,白噪声的添加准则为:
EEMD分解的次数N为:
S3、计算步骤S2得到的本征模态函数的分形维数,得到对应的分形维数曲线;
S4、将分形维数曲线与正常信号和不同故障实际曲线对比,根据曲线相似度判定振动信号所属的故障位置和故障类型。
本实施例中分形维数的计算方法为通过最小二乘法计算信号的网格维数。因为分形维数对信号识别和分类具有良好效果,可以根据不同的分量求出其特有的分形维数,以达到故障诊断的效果。两个信号越相似,其分形维数值越接近;反之当两信号差别较大时,其分形维数将有很大差距。但经过分解的实际信号往往不会像仿真信号一样为单一频率,更多情况下分解后的信号仍然是复合信号。此时分形维数不仅受频率变化的影响,同时受到不同频率信号振幅变化的影响。如果要用一层分形维数区分多种类型故障的不同严重程度将十分困难,而且容易出现误判断。所以本申请提出了一种分层的分形维数来解决故障的准确定位问题,为传动系统保养和设备更换提供合理依据,其判定逻辑结构如图2所示。根据第一层分形维数判定故障位置:行星轮故障、太阳轮故障以及正常情况。本实施例中,第一层判定中行星轮故障和太阳轮故障具体选择点蚀故障。由第二层分形维数在故障位置基础上判定具体的故障严重程度,即故障类型,包括:断齿、点蚀、磨损和缺齿。
所以步骤S4具体包括:
S41、将第一层分形维数曲线与正常信号和表示不同故障位置的实际曲线对比,若与第一层分形维数曲线最相似的实际曲线为故障曲线,则将其对应的故障位置作为振动信号所属的故障位置,进行步骤S42,否则判定振动信号为正常信号,结束判定;
S42、将第二层分形维数曲线与已确定故障位置的不同故障类型的实际曲线对比,将最相似的实际曲线对应的故障类型作为振动信号所属的故障类型。
下面通过仿真实验对本方法进行验证。在风电机组运行中,行星齿轮箱正常工作时一般周期性地运转,所以可以使用正弦信号代表正常的旋转振动信号;行星轮本身运转具有通过效应,将对啮合信号产生调制作用形成调制信号;当齿轮箱出现故障时,通常表现为高频冲击信号,呈现的特性为冲击和调制;实际运行中通常伴随大量噪声的存在,所以噪声信号必不可少。因此,本实施例中,仿真一组由高频冲击信号、调制信号、周期信号以及噪声信号组成的模拟故障信号,如图3的分量1~4所示。
经过传统EEMD和本方法提出的EEMD分解的结果分别如图4和图5所示。可以看出本方法提出的EEMD分解结果中的imf1、imf2、imf3分别与仿真信号中的分量1、2、3相似,而传统EEMD分解结果中是imf2、imf4、imf6分别与仿真信号中分量1、2、3相似,说明本方法所提出的自适应EEMD方法不仅很好的解决了EEMD本身的模态混叠问题,同时提高了分解效果并缩短了分解所用时间,很大程度上提高EEMD的综合性能,为后续的信号识别打下基础。
如图6所示,本实施例将得到的第一层分形维数曲线与正常信号、太阳轮点蚀和行星轮点蚀曲线对比,可以清楚的对待诊断信号进行定位,所以根据曲线变化趋势和数值的相似性将故障定位为太阳轮故障。
如图7所示,将第二层分形维数曲线与太阳轮点蚀信号、太阳轮均匀磨损信号、太阳轮断齿信号、太阳轮缺齿信号曲线对比。由于四种故障类型的故障特征接近,所以四条曲线的形状相似度很高,无法通过曲线变化趋势进行识别,而选择通过比较曲线两两间的平均距离来对待诊断信号进行识别,从局部放大图中可以看出待诊断信号与太阳轮断齿信号分形维数曲线接近重合,而与太阳轮其他类型故障则差别较大,所以可以确定故障类型为太阳轮断齿。
使用两层分形维数方法,成功对仿真实验中待诊断信号进行了故障识别与诊断,证明了本方法能有效定位故障位置、确定故障类型。
Claims (8)
1.一种风电机组行星齿轮箱故障诊断方法,其特征在于,包括以下步骤:
S1、利用奇异值差分谱法对采集的振动信号中的高频冲击信号和噪声信号的复合分量进行提取;
S2、对采集的振动信号添加白噪声进行N次EEMD分解,得到一系列本征模态函数,其中,白噪声的添加准则为:
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>e</mi>
<mo>&le;</mo>
<mfrac>
<msub>
<mo>&part;</mo>
<mi>c</mi>
</msub>
<mn>3</mn>
</mfrac>
</mrow>
e表示添加白噪声的标准差,表示经过奇异值差分谱法提取出的高频冲击分量和噪声分量信号复合分量的标准差;
S3、计算步骤S2得到的本征模态函数的分形维数,得到对应的分形维数曲线;
S4、将分形维数曲线与正常信号和不同故障实际曲线对比,根据曲线相似度判定振动信号所属的故障位置和故障类型。
2.根据权利要求1所述的一种风电机组行星齿轮箱故障诊断方法,其特征在于,所述白噪声e的添加准则为:
<mrow>
<mi>e</mi>
<mo>=</mo>
<mfrac>
<msub>
<mo>&part;</mo>
<mi>c</mi>
</msub>
<mn>3</mn>
</mfrac>
<mo>.</mo>
</mrow>
3.根据权利要求1所述的一种风电机组行星齿轮箱故障诊断方法,其特征在于,所述EEMD分解的次数N为:
<mrow>
<mi>N</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>e</mi>
<mi>&epsiv;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
其中,ε表示允许分解误差。
4.根据权利要求3所述的一种风电机组行星齿轮箱故障诊断方法,其特征在于,所述允许分解误差ε为2%。
5.根据权利要求1所述的一种风电机组行星齿轮箱故障诊断方法,其特征在于,所述分形维数的计算方法为通过最小二乘法计算信号的网格维数。
6.根据权利要求1所述的一种风电机组行星齿轮箱故障诊断方法,其特征在于,所述步骤S4具体包括:
S41、将第一层分形维数曲线与正常信号和表示不同故障位置的实际曲线对比,若与第一层分形维数曲线最相似的实际曲线为故障曲线,则将其对应的故障位置作为振动信号所属的故障位置,进行步骤S42,否则判定振动信号为正常信号,结束判定;
S42、将第二层分形维数曲线与已确定故障位置的不同故障类型的实际曲线对比,将最相似的实际曲线对应的故障类型作为振动信号所属的故障类型。
7.根据权利要求1所述的一种风电机组行星齿轮箱故障诊断方法,其特征在于,所述故障位置包括:行星轮和太阳轮。
8.根据权利要求1所述的一种风电机组行星齿轮箱故障诊断方法,其特征在于,所述故障类型包括:断齿、点蚀、磨损和缺齿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711122369.XA CN107907327B (zh) | 2017-11-14 | 2017-11-14 | 一种风电机组行星齿轮箱故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711122369.XA CN107907327B (zh) | 2017-11-14 | 2017-11-14 | 一种风电机组行星齿轮箱故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107907327A true CN107907327A (zh) | 2018-04-13 |
CN107907327B CN107907327B (zh) | 2019-06-11 |
Family
ID=61845412
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711122369.XA Active CN107907327B (zh) | 2017-11-14 | 2017-11-14 | 一种风电机组行星齿轮箱故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107907327B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109297705A (zh) * | 2018-08-16 | 2019-02-01 | 东南大学 | 基于med和模糊熵的行星齿轮箱振动信号故障诊断方法 |
CN109443768A (zh) * | 2018-11-14 | 2019-03-08 | 中国直升机设计研究所 | 一种直升机主减速器行星齿轮振动信号分离方法 |
CN109580218A (zh) * | 2018-12-08 | 2019-04-05 | 上海电力学院 | 一种基于似然学习机的风机齿轮箱状态识别方法 |
CN109685253A (zh) * | 2018-12-07 | 2019-04-26 | 上海睿深电子科技有限公司 | 适用于行星齿轮的品质管理系统、方法、终端及存储介质 |
CN110727908A (zh) * | 2019-09-27 | 2020-01-24 | 宁夏凯晨电气集团有限公司 | 解决复杂电气故障的模态解析方法 |
CN112232321A (zh) * | 2020-12-14 | 2021-01-15 | 西南交通大学 | 一种振动数据干扰降噪方法、装置、设备及可读存储介质 |
US11344987B2 (en) * | 2019-09-04 | 2022-05-31 | Tsinghua Shenzhen International Graduate School | Method for monitoring chatter in machining process |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954443A (zh) * | 2014-04-29 | 2014-07-30 | 华电电力科学研究院 | 自适应总体平均经验模式分解eemd协助噪声大小确定方法 |
CN104237731A (zh) * | 2014-09-25 | 2014-12-24 | 福州大学 | 基于eemd与能量法的谐振接地配电网单相接地故障选线方法 |
CN104392444A (zh) * | 2014-11-18 | 2015-03-04 | 陕西师范大学 | 基于二维集合经验模式分解的医学mr图像特征提取方法 |
CN104748961A (zh) * | 2015-03-30 | 2015-07-01 | 中国矿业大学 | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 |
-
2017
- 2017-11-14 CN CN201711122369.XA patent/CN107907327B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954443A (zh) * | 2014-04-29 | 2014-07-30 | 华电电力科学研究院 | 自适应总体平均经验模式分解eemd协助噪声大小确定方法 |
CN104237731A (zh) * | 2014-09-25 | 2014-12-24 | 福州大学 | 基于eemd与能量法的谐振接地配电网单相接地故障选线方法 |
CN104392444A (zh) * | 2014-11-18 | 2015-03-04 | 陕西师范大学 | 基于二维集合经验模式分解的医学mr图像特征提取方法 |
CN104748961A (zh) * | 2015-03-30 | 2015-07-01 | 中国矿业大学 | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109297705A (zh) * | 2018-08-16 | 2019-02-01 | 东南大学 | 基于med和模糊熵的行星齿轮箱振动信号故障诊断方法 |
CN109443768A (zh) * | 2018-11-14 | 2019-03-08 | 中国直升机设计研究所 | 一种直升机主减速器行星齿轮振动信号分离方法 |
CN109685253A (zh) * | 2018-12-07 | 2019-04-26 | 上海睿深电子科技有限公司 | 适用于行星齿轮的品质管理系统、方法、终端及存储介质 |
CN109580218A (zh) * | 2018-12-08 | 2019-04-05 | 上海电力学院 | 一种基于似然学习机的风机齿轮箱状态识别方法 |
US11344987B2 (en) * | 2019-09-04 | 2022-05-31 | Tsinghua Shenzhen International Graduate School | Method for monitoring chatter in machining process |
CN110727908A (zh) * | 2019-09-27 | 2020-01-24 | 宁夏凯晨电气集团有限公司 | 解决复杂电气故障的模态解析方法 |
CN112232321A (zh) * | 2020-12-14 | 2021-01-15 | 西南交通大学 | 一种振动数据干扰降噪方法、装置、设备及可读存储介质 |
CN112232321B (zh) * | 2020-12-14 | 2021-03-19 | 西南交通大学 | 一种振动数据干扰降噪方法、装置、设备及可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107907327B (zh) | 2019-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107907327A (zh) | 一种风电机组行星齿轮箱故障诊断方法 | |
CN107016235B (zh) | 基于多特征自适应融合的设备运行状态健康度评估方法 | |
CN102866027B (zh) | 基于lmd和局域时频熵的旋转机械故障特征提取方法 | |
CN110259646B (zh) | 一种基于历史数据的风力发电机组部件状态预警方法 | |
CN110674842A (zh) | 一种风电机组主轴轴承故障预测方法 | |
CN104131950B (zh) | 一种风电机组温度特征量的阈值分区确定方法 | |
CN108072524B (zh) | 一种风电机组齿轮箱轴承故障预警方法 | |
CN102768115A (zh) | 一种风电机组齿轮箱健康状态实时动态监控方法 | |
CN109685371B (zh) | 基于贝叶斯网的动态概权风电机组运行状态综合评估方法 | |
CN110046562A (zh) | 一种风电系统健康监测方法及装置 | |
CN108520269A (zh) | 一种风速预测方法及风速预测系统 | |
CN101718634A (zh) | 基于多元概率模型的设备状态综合动态报警方法 | |
CN104748962A (zh) | 一种基于堆叠自动编码机的行星齿轮箱智能诊断方法 | |
Li et al. | A review on the signal processing methods of rotating machinery fault diagnosis | |
CN103926075B (zh) | 一种提高风电机组齿轮箱故障分析精度的方法 | |
CN107291991A (zh) | 一种基于动态网络标志的风电机组早期缺陷预警方法 | |
CN111077386A (zh) | 一种电气设备早期故障信号降噪方法 | |
CN113806893A (zh) | 基于工业互联网的风机状态监测与故障诊断方法及系统 | |
CN105065212A (zh) | 一种风电场风电机组校验方法及系统 | |
CN111878322B (zh) | 风力发电机装置 | |
CN115453356B (zh) | 一种动力设备运行状态监测分析方法、系统、终端及介质 | |
CN113250911B (zh) | 一种基于vmd分解算法的风机叶片故障诊断方法 | |
CN114386452A (zh) | 核电循环水泵太阳轮故障检测方法 | |
CN108242411B (zh) | 用于管控线上缺陷的方法和系统 | |
CN112766200A (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 |