一种风电齿轮箱故障诊断方法
【专利摘要】本发明提供了一种风电齿轮箱故障诊断方法。齿轮箱存在故障时其振动信号具有两个特点：1、故障信号比较微弱，容易被其他的信号所淹没。2、故障信号会因为故障冲击产生调制现象。本发明分别应用二维经验模态分解(BEMD)和希尔伯特平方解调技术解决以上两个问题。首先应用BEMD分解将原始复值振动信号分解为一系列的复值本征模态函数IMF，能准确的检测出复值信号的相位信息；然后筛选出峭度值最大的IMF分量IMFc，通常对应的就是故障特征分量；而后利用希尔伯特平方解调对IMFc的进行解调，对解调后的信号的实部与虚部分别做频谱分析；最后将频谱分析结果与齿轮箱各回转部件的回转频率对比，确定故障位置。
【专利说明】
一种风电齿轮箱故障诊断方法
技术领域
[0001] 本发明涉及一种风电齿轮箱故障诊断方法，尤其涉及一种风力发电机齿轮箱故障 诊断方法。
【背景技术】
[0002] 面对严峻的能源安全问题，以风力发电为代表的可再生能源得到飞速的发展。风 力发电机组大多安装在戈壁滩、海岛、高山等风口处，自然条件恶劣，承受的是无规律变向 风力作用及强阵风冲击的变载荷，所以风力发电机组是一个复杂、多变量、非线性的不确定 系统，机组故障时有发生。故障统计结果显示，风电系统的故障12%来自齿轮箱的失效，另 外，风电机组的安装位置交通不便，而齿轮箱多安装在距离地面几十米高的机舱内，空间狭 小，如果没有及时检测到故障，很容易导致风电齿轮箱的损坏，维修十分困难，往往需要下 塔处理，吊装维修费用高，并且会严重影响风场的经济效益，因此，风电齿轮箱的状态监测 和故障诊断在这种情况下具有重要的意义。
[0003] 随着我国大批风电机组的安装运行，风电齿轮箱的监测和诊断技术也得到了空前 的发展，但存在的主要问题是:监测对象单一，多数为一种状态量的监测；分析诊断功能不 够完善，缺乏有效的故障诊断能力;不适用于批量生产和现场安装。
[0004]现代风电齿轮箱传递功率大、结构复杂，故障信号为非平稳瞬态信号，并且故障 特征量容易被齿轮的啮合信号及其他的噪音所淹没。此外，当齿轮箱存在故障时，会产出冲 击振动，引起信号的调幅调制现象，对振动信号直接进行频谱分析时，频谱图上一般都会出 现调制载波频率而没有相应的故障频率，因此无法提取到故障信息。
【发明内容】
[0005] 为了克服上述现有技术的不足，本发明提供了一种风电齿轮箱故障诊断方法，其 目的在于快速、准确实现对齿轮箱故障地监测。
[0006] 为了实现上述目的，本发明提供的一种风电齿轮箱故障诊断方法，所述方法包括 以下步骤：
[0007] SI:通过振动信号采集系统对风电齿轮箱的振动信号进行采集，将采集到的振动 信号构造成复值信号s(t);
[0008] S2:利用BEMD对步骤SI中构造的复值信号s(t)进行分解，得到一系列的复值固有 模态函数頂Fs，这些頂Fs能准确的反应出复值信号的相位信息；
[0009] S3:计算上述各頂F的峭度值，并选出峭度值最大的頂F分量頂Fe;
[0010] S4:利用希尔伯特平方解调对頂F。进行解调；
[0011] S5:对解调后的信号的实部与虚部分别做频谱分析，比较验证频率成分是否相同。 将频谱分析结果与齿轮箱各回转部件的回转频率对比，确定故障位置。
[0012] 进一步地，所述步骤Sl中的振动信号采集系统包括三相压电加速度传感器、信号 采集卡和控制器机箱，信号采集卡与传感器采用导线相连，信号采集卡通过插拔接口插接 在控制器机箱上。三相压电式加速度传感器布置在箱体外壳、输出轴的前侧轴承、输入轴的 前侧轴承和后侧轴承上。将压电式加速传感器安装在这些位置或其它重要位置上，方便检 测出重要部位是否存在故障。
[0013] 进一步地，所述步骤S2中的所述的BEMD方法包括以下步骤：
[0014] S21:选取信号投影方向的数目N，并计算投影方向：
[0015]
(1)
[0016] S22:计算步骤Sl中得到的复制信号s(t)投影到方向和上的值：
[0017]
(2)
[0018] S23:找出(〇的局部极大值和极小值点，记录下位置和数值:{(f/", Pf ) h
[0019] S24:通过三次样条差值拟合方法对步骤S23中得到的{(if，i^) }进行插值得到包 络曲线~,.(〇;
[0020] S25:重复步骤S22-S24,直到获得所有投影方向上的包络曲线；
[0021 ] S26:计算步骤S25中得到的所有包络线的均值：
[0022]
(3)
[0023] S27:从原始信号s(t)中减去包络线均值瓦得到分量81(七）：
[0024]
(4)
[0025] S28:检测gl(t)是否满足固有模态函数的筛选停止条件，如果不满足，将gl(t)替换 成新的原始信号s(t)，并返回步骤S22;否则，将分量 81(1)作为第一个固有模态函数分量(^ (t)，用原始信号s(t)减去第一个固有模态函数分量(^(〇获得残余量n(t):
[0026] ci(t) = g(t)，ri(t) = s(t)_ci(t) (5)
[0027] S29:判断残余量n(t)是否为单调函数，若是，则进入步骤S210;若否，则将n(t)作 为新的原始信号并的返回步骤S22;
[0028] S210:重复步骤S22-S29直至获得所有的固有模态函数和一个最终的参与量^ (t);最终参与量rn(t)为振动信号s(t)的平均趋势;振动信号s(t)可以表示为：
[0029]
(6)
[0030]进一步地，所述步骤S4中利用希尔伯特平方解调对调制信号进行解调处理具体过 程为：
[0031] 调制信号具有式(7)所示的形式：
[0032] c(t)=f(t)s(t) (7)
[0033] 其中，f(t)为故障产生的低频冲击信号，s(t)为运行中正常的高频振动信号；
[0034] 构造函数X(t)和Y(t)如式(8)、（9)所示：
[0035] X(t)=Re2(c(t))+H2[Re(c(t))] (8)
[0036] Y(t)=Im2(c(t))+H2[Im(c(t))] (9)
[0037] 其中Re表示取复值信号的实部，Im表示取复制信号的虚部;H[ ·]是对信号进行希 尔伯特变换；
[0038] 根据 Bedrosian 定理：
[0039] H(Re(c(t)))=H[Re(f(t)s(t))]=Re(f(t))H[Re(s(t))] (10)
[0040] H[Im(c(t))]=H[Im(f(t)s(t))] = Im(f(t))H[Im(s(t))] (11)
[0041] 将式(10)、（11)分别带入式(8)、（9)可得：
[0042] X(t)=Re2(f(t)){Re2(s(t))+H2[Re(s(t))]} (12)
[0043] Y(t) = Im2(f(t)){Im2(s(t))+H2[Im(s(t))]} (13)
[0044] 令r(t)=s2(t)+H2[(t)]，高频振动信号s(t)的结构如式（14)、（15)所示：
[0045] Re(s(t)) =aicos(2jrfit) (14)
[0046] Im(s(t)) =a2C〇s(2jrf2t) (15)
[0047] 则s(t)的希尔伯特变换为：
[0048] H[Re(s(t)) ] =aisin(2jrfit) (16)
[0049] H[ Im(s(t)) ] =a2sin(2jrf2t) (17)
[0050] 则可得r(t):
[0051 ] /i (t) = a'( cos-(2^i/) + Uf sin ?-{2τ?[\?) = a'f (18)
[0052] r2(t) = a'j cos2(2/i/2〇 + 〇l Sxn1(ITrfit) ~ aj (19)
[0053] 将式（18)、（19)分别带入式（12)、（13)可得：
[0054] X(t)=ai2Re2(f(t)) (20)
[0055] Y(t)=a22Im2(f(t)) (21)
[0056] 对仅有低频部分f(t)的式(20)、（21)进行频谱分析;上述式中，t表示时间，心、5表 示频率。
[0057] 齿轮箱存在故障时其振动信号具有两个特点：1、故障信号比较微弱，容易被其他 的信号所淹没。2、故障信号会因为故障冲击产生调制现象。
[0058]针对故障信号的上述两个特点，本发明分别应用BEMD分解和希尔伯特平方解调技 术解决以上两个问题。首先应用BEMD分解将原始复值振动信号分解为一系列的复值本征模 态函数IMF，能准确的检测出复值信号的相位信息;然后筛选出峭度值最大的IMF分量IMF。， 通常对应的就是故障特征分量;而后利用希尔伯特平方解调对IMF。的进行解调，对解调后 的信号的实部与虚部分别做频谱分析，比较验证频率成分是否相同；最后将频谱分析结果 与齿轮箱各回转部件的回转频率对比，确定故障位置。通过本发明方法可以及时且方便地 检测出故障信息，及时对齿轮箱进行维修。
[0059]本发明中采用的BEMD方法是一种数据驱动式的自适应信号处理方法。其分解过程 是利用信号局部特征时间尺度信息，逐次将复值信号分解为若干个复值的固有模态函数， 每个固有模态函数反应了信号中内嵌的旋转模式，且其性质完全由信号本身性质决定，于 其他因素无关。该方法能够处理以傅里叶变换为基础的那些方法所不能处理的非线性非平 稳的数据，且能够比EMD更高效的来处理二维的信号，并反应出信号的相位偏移信息。。与 小波变换相比，因为其基函数是由信号本身分解得到的，并不需要事先选定小波基和分解 尺度，所以它是直观的、后验的和自适应的分解方法。
【附图说明】
[0060]下面结合附图对本发明作进一步说明：
[0061] 图1为四个复值信号的时域波形图；
[0062] 图2为四个复值信号直接加和后的时域波形图；
[0063] 图3为对构造的复值信号进行傅里叶变换得到的频谱结果图；
[0064] 图4为BEMD分解后的頂Fs结果图；
[0065] 图5为对解调后的信号的实部和虚部分别做频谱分析的结果图；
[0066] 图6为本发明提供的一种齿轮箱故障诊断方法的流程图；
【具体实施方式】
[0067]该齿轮箱故障诊断方法是基于齿轮箱运行过程中所产生的振动信号进行的。齿轮 箱在运行过程中无论是否存在故障都会有振动信号的产生，但是无故障和有故障时的振动 信号有不同的区别。齿轮箱没有故障时，振动信号的频率成分主要是齿轮的啮合频率和轴 承的自身的固有频率。当齿轮箱中的齿轮和轴承存在故障时，其运行过程中会因为故障产 生冲击效应，从而使振动信号发生改变，此时振动信号的频率成分不仅有齿轮的啮合频率 和轴承自身的固有频率，还会增加故障部件所在部件的回转频率。
[0068] 通过三相加速度传感器测得的振动加速度信号是时域信号。如图1所示；图1中的 复值信号是下面四个实部信号x(t)和四个虚部信号y(t)构造后的时域波形图，黑色实线表 示的是实部信号，黑色虚线表示的是虚部信号。
[0069]
[0070]
[0071]
[0072]
[0073]
[0074] 其中，X1(tWPyi(t)为模拟轴承外圈故障信号和系统结构谐振信号的调制;X2(t)、 x3(tWPy2(t)、y3(t)为周期性谐波信号成分;X4(t)和y 4(t)为模拟噪声信号;衰减频率(与系 统阻尼特性相关)α = 400;结构谐振频率5 = 800;高速级齿轮的啮合频率f2 = 100;低速级齿 轮的啮合频率f 3 = 50，故障特征频率fB = 33。
[0075]图2的信号是以上四个实部信号x(t)和四个虚部信号y(t)分别直接加和构造的复 值信号的结果。黑色实线表示的是实部信号，黑色虚线表示的是虚部信号。
[0076]诊断故障需要得到的是振动信号所对应的频率值。对于简单的信号，如上面的X2 (t)、X3(tWPy2(t)、y3(t)，直接进行傅里叶变换就可以得到振动信号的频率值，但是对振动 信号 X1(tWPyi(t)这种较为复杂的调制信号，直接进行傅里叶变换得到频率值是心=800， 但是实际中想要得到的是频率f B = 33。对振动信号X1(tWPyi(t)经过解调后，再进行傅里叶 变换后就能得到故障频率f B = 33。如图3所示，是对构造的信号直接进行傅里叶变换之后的 频谱图，从实部与虚部的频谱结果中都可以得到50Hz ,100Hz和800Hz的频率成分，但显示不 出33Hz的频率成分。
[0077]啮合频率和回转频率:假如一对啮合的齿轮对，主动轮的转速是20r/s，齿数是40， 从动轮的齿数是80,那么从动轮的转速可以根据主动轮计算得出，即10r/s。啮合频率是用 齿轮对中任意一个齿轮的回转频率乘以齿数，即主动轮啮合频率为20 X 40 = 800Hz，从动轮 啮合频率为10 X 80 = 800Hz。回转频率是指单个齿轮的回转频率，如主动轮20Hz，从动轮 10Hz。这样如果故障诊断所得结果是20Hz，说明是主动轮故障；如果是10Hz，说明是从动轮 故障。
[0078] 故障诊断的过程就是当齿轮箱中存在故障时，确定故障频率，进而确定出故障位 置的过程。
[0079] 齿轮箱存在故障时其振动信号具有两个特点：
[0080] 1、故障信号比较微弱，容易被其他的信号所淹没。
[0081] 2、故障信号会因为故障的原因产生调制现象。
[0082] 针对故障信号的上述两个特点，故障诊断方法分别应用BEMD分解和希尔伯特平方 解调技术解决以上两个问题。首先应用BEMD分解将原始复值振动信号分解为一系列的复值 本征模态函数頂F，在这些MF中，峭度值最大的那个頂F通常对应的就是故障特征分量。例 如从图2中的合信号很难直接找到故障特征信号，而对其进行BEMD分解后就能得到图4所示 的结果。其中IMF2就是故障特征信号，他的峭度值会比其他的IMF分量的大。从而就找到了 故障特征分量。然后寻找故障频率，进而确定出故障位置。用到的就是希尔伯特平方解调技 术。对信号进行解调后再进行傅里叶变换就能得到故障频率，进而确定出故障位置。如图5 所示，即对解调信号的实部和虚部分别进行傅里叶变换之后得到的频谱结果图，可以发现 大约在33Hz左右，实部与虚部的频率幅值都达到最大值。
[0083]参阅图6右侧，本发明中的方法主要包括:信号采集、信号处理和结果显示。
[0084] 参阅图6左侧，一种风电齿轮箱故障诊断方法，具体包括以下步骤：
[0085] SI:通过振动信号采集系统对风电齿轮箱的振动信号进行采集，将采集到的振动 信号构造成复值信号s(t);
[0086] S2:利用BEMD对步骤SI中构造的复值信号s(t)进行分解，得到一系列的复值固有 模态函数頂Fs，这些頂Fs能准确的反应出复值信号的相位偏移信息；
[0087] S3:计算上述各頂F的峭度值，并选出峭度值最大的頂F分量頂Fc;
[0088] S4:利用希尔伯特平方解调对IMF。进行解调；齿轮箱中存在故障时出现信号的调 制信号具有式(7)所示的形式：
[0089] c(t)=f(t)s(t) (7)
[0090] 其中，f(t)为故障产生的低频冲击信号，s(t)为运行中正常的高频振动信号；
[0091] 采用希尔伯特平方解调对调制信号进行解调处理，具体过程如下：
[0092] 构造函数X(t)和Y(t)如式(8)、（9)所示：
[0093] X(t)=Re2(c(t))+H2[Re(c(t))] (8)
[0094] Y(t) = Im2(c(t))+H2[Im(c(t))] (9)
[0095] 其中Re表示取复值信号的实部，Im表示取复制信号的虚部;H[ ·]是对信号进行希 尔伯特变换；
[0096] 根据 Bedrosian 定理：
[0097] H(Re(c(t)))=H[Re(f(t)s(t))]=Re(f(t))H[Re(s(t))] (10)
[0098] H[Im(c(t))]=H[Im(f(t)s(t))] = Im(f(t))H[Im(s(t))] (11)
[0099] 将式（10)、（11)分别带入式(8)、（9)可得：
[0100] X(t)=Re2(f(t)){Re2(s(t))+H 2[Re(s(t))]} (12)
[0101] Y(t) = Im2(f(t)){Im2(s(t))+H2[Im(s(t))]} (13)
[0102] 令r(t)=s2(t)+H2[(t)]，高频振动信号s(t)的结构如式（14)、（15)所示：
[0103] Re(s(t)) =aicos(2jrfit) (14)
[0104] Im(s(t)) =a2C〇s(2jrf2t) (15)
[0105] 则s(t)的希尔伯特变换为：
[0106] H[Re(s(t)) ] =aisin(2jrfit) (16)
[0107] H[ Im(s(t)) ] =a2sin(2jrf2t) (17)
[0108] 则可得 r(t):
[0109] (18)
[oho] .(19)
[0111] 将式(18)、（19)分别带入式(12)、（13)可得：
[0112] X(t)=ai2Re2(f(t)) (20)
[0113] Y(t)=a22Im2(f(t)) (21)
[0114] 对仅有低频部分f(t)的式(20)、（21)进行频谱分析;上述式中，t表示时间，心、5表 示频率。
[0115] S5:对解调后的信号的实部与虚部分别做频谱分析，比较验证频率成分是否相同。 将频谱分析结果与齿轮箱各回转部件的回转频率对比，确定故障位置。
[0116] 采集模块包括三相压电式加速度传感器、信号采集卡和控制器机箱，信号采集卡 与传感器采用导线相连，信号采集卡通过插拔接口插接在控制器机箱上。三相压电式加速 度传感器布置在箱体外壳、输出轴的前侧轴承、输入轴的前侧轴承和后侧轴承等重要位置 上。
[0117] c(t)=f(t)s(t)是对希尔伯特平方解调技术进行理论上解释，这里的c(t)是一个 泛指，可以指任意一个类似于图l* X1(t)Syi(t)的调制信号。而頂Fc是实际应用中，齿轮箱 故障时所产生的调制信号。
[0118] 以上所述仅为本发明的具体实施例，但本发明的技术特征并不局限于此，任何本 领域的技术人员在本发明的领域内，所作的变化或修饰皆涵盖在本发明的专利范围之中。
【主权项】
1. 一种风电齿轮箱故障诊断方法，其特征在于，所述方法包括以下步骤： S1:通过振动信号采集系统对风电齿轮箱的振动信号进行采集，将采集到的振动信号 构造成复值信号S(t); S2:利用BEMD对步骤S1中构造的复值信号s(t)进行分解，得到一系列的复值固有模态 函数頂Fs，这些頂Fs能准确的反应出复值信号的相位偏移信息； S3:计算上述各頂F的峭度值，并选出峭度值最大的頂F分量頂Fc; S4:利用希尔伯特平方解调对IMF。进行解调； S5:对解调后的信号的实部与虚部分别做频谱分析，比较验证频率成分是否相同。将频 谱分析结果与齿轮箱各回转部件的回转频率对比，确定故障位置。2. 根据权利要求1所述的风电齿轮箱故障诊断方法，其特征在于，所述步骤S1中的振动 信号采集系统包括三相压电加速度传感器、信号采集卡和控制器机箱，信号采集卡与传感 器采用导线相连，信号采集卡通过插拔接口插接在控制器机箱上。三相压电式加速度传感 器布置在箱体外壳、输出轴的前侧轴承、输入轴的前侧轴承和后侧轴承上。3. 根据权利要求2所述的风电齿轮箱故障诊断方法，其特征在于，所述步骤S2中的所述 的BEMD方法包括以下步骤： S21:选取信号投影方向的数目N，并计算投影方向： - Imn/N, I <m< N (1) S22:计算步骤S1中得到的复制信号s(t)投影到方向％上的值： (t) = Re[e-，v'- .y(r)j (2) S23:找出(?)的局部极大值和极小值点，记录下位置和数值:,尸r)}; S24:通过三次样条差值拟合方法对步骤S23中得到的{(?/%/?/"）}进行插值得到包络曲 线 ％" (f); S25:重复步骤S22-S24,直到获得所有投影方向上的包络曲线； S26:计算步骤S25中得到的所有包络线的均值：S27:从原始信号s(t)中减去包络线均值网]")得到分量幻(〇: g, (〇 = s(〇 - m,(f) (4) S28:检测gl(t)是否满足固有模态函数的筛选停止条件，如果不满足，将gl(t)替换成新 的原始信号s(t)，并返回步骤S22;否则，将分量81(1)作为第一个固有模态函数分量(^(〇， 用原始信号s(t)减去第一个固有模态函数分量(^(〇获得残余量n(t): ci(t)=g(t)，ri(t) = s(t)-ci(t) (5) S29:判断残余量n(t)是否为单调函数，若是，则进入步骤S210;若否，则将n(t)作为新 的原始信号并的返回步骤S22; S210:重复步骤S22-S29直至获得所有的固有模态函数和一个最终的参与量化(〇;最终 参与量rn(t)为振动信号s(t)的平均趋势;振动信号s(t)可以表示为：4.根据权利要求3所述的风电齿轮箱故障诊断方法，其特征在于，所述步骤S4中齿轮箱 存在故障时会出现信号的调制现象，调制信号具有式(7)所示的形式： c(t)=f (t)s(t) (7) 其中，f(t)为故障产生的低频冲击信号，s(t)为运行中正常的高频振动信号； 采用希尔伯特平方解调对调制信号进行解调处理，具体过程如下： 构造函数X(t)和Y(t)如式(8)、（9)所示： X(t)=Re2(c(t))+H2[Re(c(t))] (8) Y(t) = Im2(c(t))+H2[Im(c(t))] (9) 其中Re表示取复值信号的实部，Im表示取复制信号的虚部;H[ ?]是对信号进行希尔伯 特变换； 根据Bedrosian定理： H(Re(c(t)))=H[Re(f(t)s(t))]=Re(f(t))H[Re(s(t))] (10) H[Im(c(t))]=H[Im(f(t)s(t))] = Im(f(t))H[Im(s(t))] (11) 将式(1〇)、（11)分别带入式(8)、（9)可得： X(t)=Re2(f(t)){Re2(S(t))+H2[R e(S(t))]} (12) Y(t) = Im2(f(t)){Im2(s(t))+H2[Im(s(t))]} (13) 令r(t) = s2(t)+H2[(t)]，高频振动信号s(t)的结构如式(14)、（15)所示： Re(s(t)) = aicos(23ifit) (14) Im(s(t)) = a2C〇s(23if2t) (15) 则s(t)的希尔伯特变换为： H[Re(s(t)) ] =aisin(23ifit) (16) H[Im(s(t)) ] =a2sin(23if2t) (17) 则可得r(t):将式(18)、（19)分别带入式(12)、（13)可得： X(t)=ai2Re2(f(t)) (20) Y(t)=a22Im2(f (t)) (21) 对仅有低频部分f(t)的式(20)、（21)进行频谱分析;上述式中，t表示时间，表示频 率。
【文档编号】G01M13/02GK106055734SQ201610322826
【公开日】2016年10月26日
【申请日】2016年5月13日
【发明人】沈建洋, 陈换过, 陈文华, 陈培, 易永余, 钱佳诚
【申请人】浙江理工大学