CN109580218B - 一种基于似然学习机的风机齿轮箱状态识别方法 - Google Patents
一种基于似然学习机的风机齿轮箱状态识别方法 Download PDFInfo
- Publication number
- CN109580218B CN109580218B CN201811498568.5A CN201811498568A CN109580218B CN 109580218 B CN109580218 B CN 109580218B CN 201811498568 A CN201811498568 A CN 201811498568A CN 109580218 B CN109580218 B CN 109580218B
- Authority
- CN
- China
- Prior art keywords
- gaussian
- model
- state
- gear box
- kurtosis
- 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
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/025—Test-benches with rotational drive means and loading means; Load or drive simulation
-
- 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)从风机齿轮箱振动信号中提取出峭度作为风机齿轮箱的故障特征量;2)利用似然学习机对峭度进行学习生成高斯模型作为状态监测模型,即高斯监测模型;3)利用生成的高斯监测模型对齿轮箱状态进行判断。与传统的信号处理进行故障诊断的方法相比,传统方法进行故障识别需要依赖运行人员的经验,本发明不依赖运行人员的经验,具有智能识别的特点;与现有智能算法相比,现有智能算法易出现过拟合的问题,本发明泛化性好,具有很好的可移植性。
Description
技术领域
本发明涉及电力系统故障监测与诊断领域,尤其是涉及一种基于似然学习机的风机齿轮箱状态识别方法。
背景技术
随着风电装机容量不断增大,风电场设备的维护使用定期维护、事后维修的方式,其维护费用占风电场总运行费用高达30%,行星齿轮箱作为风机传动系统的关键部件,长期运行在交变载荷下,其故障率超过风机机组总故障的60%,如何对风机齿轮箱进行运行状态监测与故障诊断以降低运维成本、提高风电机组安全性与可靠性是电力系统故障监测与诊断领域的热点。
风机传动系统中存在的不确定因素与非线性因素造成振动信号分析困难。正常的齿轮运行时,齿轮之间啮合完好,齿轮滑动接触表面之间存在润滑油膜,振动信号平稳。当齿轮发生局部故障,带有损伤的轮齿与其他轮齿啮合,导致润滑油膜破裂,产生振动冲击现象,冲击信号将随着齿轮的旋转按照一定的时间间隔重复出现。周期性或准周期性的振动冲击是齿轮箱故障的一个关键特征。
齿轮箱状态监测与故障诊断方法主要可以分为信号处理和智能算法两类。信号处理方法有:1)齿轮在工作过程中受力会出现方向不确定性的问题,利用全矢频带能量分析的方法可用于判定齿轮箱的早期微弱故障;2)针对EMD出现模态混叠,在采集信号中加入白噪声,根据信号本身特性自适应设定白噪声标准差改进分解,达到最优分解效果提高了故障诊断的效率;3)应用自适应最大相关峭度解卷积的方法解决早起信号微弱难以提取的问题。智能算法有:1)应用模糊C均值聚类对信号进行聚类可以提高故障诊断的准确度;2)隐马尔可夫模型算法应用于齿轮箱诊断系统能够提高平均识别度;3)卷积自动编码器进行齿轮箱故障特征是一种无监督的学习方式,具有效率高的特点。
目前,信号分析需要结合运行人员的经验进行故障判断与识别因而效率低,智能算法存在过拟合泛化性能低的问题。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于似然学习机的风机齿轮箱状态识别方法。
本发明的目的可以通过以下技术方案来实现:
一种基于似然学习机的风机齿轮箱状态识别方法,包括以下步骤:
1)从风机齿轮箱振动信号中提取出峭度作为风机齿轮箱的故障特征量;
2)利用似然学习机对峭度进行学习生成高斯模型作为状态监测模型,即高斯监测模型;
3)利用生成的高斯监测模型对齿轮箱状态进行判断。
所述的步骤1)中,峭度K的表达式为:
其中,s为信号值,μ为信号均值,σ为信号标注差,E(·)为期望。
所述的步骤2)具体为:
采用似然学习机对峭度数据进行概率密度估计并选择高斯模型的概率分布,将其作为先验概率,然后对新的峭度数据进行多次学习,每次学习将上一次的后验概率作为下一次的先验概率,最终得到学习好的状态监测模型。
所述的步骤2)中,高斯监测模型包括单一高斯模型和高斯混合模型,当风机齿轮箱正常时,高斯监测模型为单一高斯模型,当风机齿轮箱为故障时,高斯监测模型为混合系数为2的高斯混合模型。
所述故障的类型包括磨损、断齿和齿根裂缝。
当风机齿轮箱为正常状态时,其对应单一高斯模型的特征量峭度值为3,当风机齿轮箱为故障时,磨损、断齿和齿根裂缝对应高斯混合模型的特征量峭度值均大于3,且依次增大。
与现有技术相比,本发明具有以下优点:
本发明针对信号分析和人工智能算法存在的弊端,从概率论与数理统计角度进行机器学习研究,提出了似然学习机(likelihood learning machine,LLM)的学习算法能够自主学习风机齿轮箱的运行模式,降低了对运行人员的依赖,且其具有良好的泛化性能得以应对风机运行的各种复杂工况,达到对风机齿轮箱监测与故障诊断的目的,降低运维成本,提高经济性和风机运行的安全性和可靠性。
附图说明
图1为行星齿轮箱行星轮振动波形图,其中,图(1a)为正常状态的振动波形图,图(1b)为磨损状态的振动波形图,图(1c)为齿根断裂状态的振动波形图,图(1d)为断齿状态的振动波形图。
图2为转速1800r/min状态监测图。
图3为转速2400r/min状态监测图。
图4为转速3000r/min状态监测图。
图5为本发明的方法流程图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
本发明旨对复杂工况下的风机齿轮箱监测与故障诊断,首先从风机齿轮箱振动信号提取出峭度作为特征量,然后利用似然学习机对采集信号进行学习生成高斯模型,最后利用所生成的高斯模型可以对齿轮箱状态进行评估。
具体采取的技术方案如下:
(1)齿轮箱故障特征量提取。风机的实际运行过程中,由于存在振动信号传输路径复杂多变、运行环境噪声干扰严重、振动源与噪声互相耦合等诸多因素的影响,齿轮箱故障检测困难。峭度是传统的时域统计指标中使用最普遍,最可靠的指标之一,峭度(Kurtosis)是反映振动信号分布特性的数值统计量及描述波形尖峰度的无量纲参数,是归一化的四阶中心矩,数学描述为:
(2)似然学习机理论(likelihood learning machine)
采集一组振动信号数据为S=(s1,s2,...sn)T,用ω为控制参数的条件概率p(S|ω)来表示,根据贝叶斯定理可以表示为:
通过后验概率p(ω|S),对信号S进行观测可对ω进行量化的不确定性估计。变量p(S|ω)作为控制变量ω的函数,表达在不同ω下信号S出现可能性,称为似然函数(likelihood function)。利用似然函数,贝叶斯定理可以表述为:
posterior∝likelihood×prior (3)
posterior为后验概率,likelihood为似然函数,prior为先验概率。
将先验概率与后验概率作为表征两种系统状态的状态模型,将上式进行进一步改写:
posterior state∝likelihood learning machine*prior state (4)
posterior state称为后验状态,表示学习之后的状态;likelihood learningmachine被称为似然学习机;prior state称为先验状态,表示学习之前的状态。
从上式可以看出,最大似然函数将先验状态与后验状态进行链接,前后两种状态进行更新是似然函数学习的结果,将似然函数看作学习器,称之为似然学习机(likelihoodlearning machine,LLM),似然学习机能够对采集数据进行学习估计,以样本数据出现最大概率进行参数学习。
似然学习机算法首先对观测数据进行概率密度估计并选择合适的概率分布,将其作为先验概率即学习前的状态;然后应用似然学习机对新观察到的数据进行学习,先验概率转化为后验概率即得到学习后状态;每经过一次学习,状态模型便得到更新,原后验状态就作为下一次学习的先验状态。如此反复进行学习达到对运行状态不断监测的目的,依据不断学习的状态模型对风机维护做出最优决策。
式(4)中先验模型与后验模型的选择需要根据研究对象分布特性进行确定。本文研究对象为风机齿轮箱振动信号,根据上文特征提取研究,其振动信号的高斯分布特性,故选择高斯模型作为状态监测模型。
(3)高斯混合模型
高斯模型具有良好的分布统计特性,其在状态估计与故障预测方面有巨大优势,既可提供预测点的估计值,又可给出预测结果的概率密度分布即不确定性估计。通过将单一高斯模型进行线性组合,称为混合高斯模型,调节均值和方差以及线性组合的系数,可以对连续的概率密度给出任意精度近似,形式为:
N(s|μk,Σk)称为混合高斯模型成分,πk称为混合系数。对上式两侧进行积分,p(s)和各个高斯成分都是归一化的,得:
给定N(s|μk,Σk),满足p(s)≥0的充分必要条件是对于所有的k都有πk≥0,将此条件与条件式(13)结合,得:
0≤πk≤1 (7)
混合高斯模型似然学习机因引入混合系数的奇异性造成其参数无法直接求解,需使用EM算法进行求解。
(4)似然学习机的EM算法
似然学习机的EM算法求解混合高斯模型的步骤为:
(a)初始化高斯监测模型参数均值μk、协方差Σk、混合系数πk;
(b)E步骤,使用当前参数值计算监测隐变量“责任”
(c)M步骤,使用当前责任更新监测模型参数,
(d)计算似然监测模型的对数形式
检查其收敛性,若不收敛返回步骤(b)。
(5)利用似然学习机进行状态监测的具体步骤为:
(a)初始化状态监测模型;
(b)对采集信号进行峭度特征量的提取;
(c)利用似然学习机对峭度进行学习;
(d)对状态监测模型进行更新;
(e)根据状态监测模型对风机齿轮箱状态进行评估。
实施例:
为验证基于本发明的有效性,采用SpectraQuest Inc公司的DrivetrainDynamics Simulator进行实验验证,实验平台由6个部分组成,分别为:驱动电机,2级制行星齿轮箱,平行齿轮箱,模拟负载,传感器,转速控制器。实验所用的采样率为12000Hz。振动传感器安装在齿轮箱轴向方向。行星轮的故障的常见类型包括磨损、齿根有裂缝和断齿故障。
当电机转速为1800r/min时,行星齿轮在不同状态下的原始振动信号如图(1)所示。正常状态下,齿轮啮合完好,振动信号平稳,幅值小,没有明显的冲击信号;磨损状态下,磨损导致齿轮表面的油膜破裂,啮合度降低,振动信号幅值增大,出现准周期性的冲击信号;齿根裂缝是在行星轮齿根部出现了裂缝,裂缝同样导致油膜破裂,啮合度降低,相对正常状态下,其振动幅值增大,出现了周期的冲击信号;断齿状态下,与磨损与齿根断裂相比较,冲击信号出现的频次更低,且周期性特征明显。
正常状态下,采集信号的峭度值作为特征量进行特征提取;将特征量输入到似然学习机内,利用EM算法进行状态模型求解;得到正常状态下的状态监测图如图2中所示,正常状态下得到的监测模型为单一高斯模型,其众数值落在3附近,表明风机齿轮箱运行平稳。高斯模型表明峭度分布在3附近的占比高即风机处在正常状态,扰动造成的峭度值偏离3的数值占比少,扰动存在的概率小。
将磨损状态下的数据作为输入进行实验,得到的状态监测图如图2中所示。在磨损状态下,因为振动冲击的存在,特征量峭度值增大,EM算法求解得到了混合系数K为2的混合高斯模型。与正常状态下的单一峰值相比,混合模型出现两个峰值,两者存在明显差异,可用于进行状态监测。大峭度值表示冲击成分增多,且大峭度值在整体分布中的占比增多,表明出现故障概率增大。监测结果表明,大峭度值的占比增大时,风机齿轮箱出现故障。
将断齿状态下的数据作为输入进行实验,得到的状态监测图如图2中所示。与磨损状态相同,出现了K为2的混合高斯模型,但相比磨损状态,其出现的高峭度峰值大于磨损状态,表明断齿状态下,振动冲击出现的占比高于磨损状态,与原始状态图中结果分析一致。
将齿根裂缝状态下的数据作为输入进行实验,得到的状态监测图如图2中所示。与磨损、断齿两种状态相同,出现了K为2的混合高斯模型,但相比前两种状态,其出现的高峭度峰值更大,表明其冲击振动更剧烈,与原始状态图中结果分析一致。
进行电机转速为2400r/min的实验,状态监测图如图3所示。输入正常状态下的采集数据,得到K=1的监测模型,表明风机齿轮箱处于正常状态。分别输入3种故障状态下的采集数据,分别得到K=2的监测模型,可做出风机齿轮箱故障的报警。故障状态类别可根据分布峰值所处大小进行判断。
进一步将电机转速为3000r/min的试验数据输入,得到状态监测图如图4所示。正常与故障数据分别得到K=1与K=2的状态模型,能够给出故障警报。
根据上述实施实例可以看出,本发明能够在不同的转速下辨识齿轮的状态,当出现故障时,不仅能够有效地做出预警,而且能够判断故障类型,可以有效地进行风机齿轮箱运行状态监测,具有良好的实用价值。
Claims (1)
1.一种基于似然学习机的风机齿轮箱状态识别方法,其特征在于,包括以下步骤:
1)从风机齿轮箱振动信号中提取出峭度作为风机齿轮箱的故障特征量,峭度K的表达式为:
其中,s为信号值,μ为信号均值,σ为信号标注差,E(·)为期望;
2)利用似然学习机对峭度进行学习生成高斯模型作为状态监测模型,即高斯监测模型,高斯监测模型包括单一高斯模型和高斯混合模型,当风机齿轮箱正常时,高斯监测模型为单一高斯模型,当风机齿轮箱为故障时,高斯监测模型为混合系数为2的高斯混合模型,具体为:
采用似然学习机对峭度数据进行概率密度估计并选择高斯模型的概率分布,将其作为先验概率,然后对新的峭度数据进行多次学习,每次学习将上一次的后验概率作为下一次的先验概率,最终得到学习好的状态监测模型;
3)利用生成的高斯监测模型对齿轮箱状态进行判断,所述故障的类型包括磨损、断齿和齿根裂缝,当风机齿轮箱为正常状态时,其对应单一高斯模型的特征量峭度值为3,当风机齿轮箱为故障时,磨损、断齿和齿根裂缝对应高斯混合模型的特征量峭度值均大于3,且依次增大。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811498568.5A CN109580218B (zh) | 2018-12-08 | 2018-12-08 | 一种基于似然学习机的风机齿轮箱状态识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811498568.5A CN109580218B (zh) | 2018-12-08 | 2018-12-08 | 一种基于似然学习机的风机齿轮箱状态识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109580218A CN109580218A (zh) | 2019-04-05 |
CN109580218B true CN109580218B (zh) | 2020-12-22 |
Family
ID=65927919
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811498568.5A Active CN109580218B (zh) | 2018-12-08 | 2018-12-08 | 一种基于似然学习机的风机齿轮箱状态识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109580218B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110794683A (zh) * | 2019-11-27 | 2020-02-14 | 中国大唐集团科学技术研究院有限公司华中电力试验研究院 | 基于深度神经网络和峭度特征的风电齿轮箱状态评估方法 |
CN111721528B (zh) * | 2020-05-18 | 2022-04-05 | 浙江工业大学 | 基于cms系统大数据的风力发电机组齿轮箱故障预警方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103234753A (zh) * | 2013-04-11 | 2013-08-07 | 华北电力大学 | 基于高斯过程建模的风电机组轴系状态监测方法 |
CN108536971A (zh) * | 2018-04-13 | 2018-09-14 | 广州市建筑科学研究院有限公司 | 一种基于贝叶斯模型的结构损伤识别方法 |
CN108709745A (zh) * | 2018-05-02 | 2018-10-26 | 东北林业大学 | 一种基于增强型lpp算法和极限学习机快速轴承故障识别方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102434387A (zh) * | 2011-11-16 | 2012-05-02 | 三一电气有限责任公司 | 风机检测诊断系统 |
KR101300037B1 (ko) * | 2012-02-01 | 2013-08-29 | 동아대학교 산학협력단 | 초광대역 다중 사용자 시스템에서 일반화된 가우시언 라플라시안 분포 모델에 기반한 향상된 uwb 수신 방법 |
CN103808509A (zh) * | 2014-02-19 | 2014-05-21 | 华北电力大学(保定) | 一种基于人工智能算法的风机齿轮箱故障诊断方法 |
CN105136454A (zh) * | 2015-10-15 | 2015-12-09 | 上海电机学院 | 一种风电机组齿轮箱故障识别方法 |
CN105699080B (zh) * | 2015-12-18 | 2019-09-03 | 华北电力大学(保定) | 一种基于振动数据的风电机组轴承故障特征提取方法 |
CN105719023A (zh) * | 2016-01-24 | 2016-06-29 | 东北电力大学 | 一种基于混合高斯分布的风电功率实时预测误差分析方法 |
CN107563251B (zh) * | 2016-07-01 | 2021-11-09 | 华北电力大学(保定) | 基于极限学习机的风机故障诊断方法 |
CN107907327B (zh) * | 2017-11-14 | 2019-06-11 | 上海电力学院 | 一种风电机组行星齿轮箱故障诊断方法 |
-
2018
- 2018-12-08 CN CN201811498568.5A patent/CN109580218B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103234753A (zh) * | 2013-04-11 | 2013-08-07 | 华北电力大学 | 基于高斯过程建模的风电机组轴系状态监测方法 |
CN108536971A (zh) * | 2018-04-13 | 2018-09-14 | 广州市建筑科学研究院有限公司 | 一种基于贝叶斯模型的结构损伤识别方法 |
CN108709745A (zh) * | 2018-05-02 | 2018-10-26 | 东北林业大学 | 一种基于增强型lpp算法和极限学习机快速轴承故障识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109580218A (zh) | 2019-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jardine et al. | A review on machinery diagnostics and prognostics implementing condition-based maintenance | |
Zhang et al. | Fault diagnosis and prognosis using wavelet packet decomposition, Fourier transform and artificial neural network | |
Soualhi et al. | Prognosis of bearing failures using hidden Markov models and the adaptive neuro-fuzzy inference system | |
Feng et al. | A novel vibration-based prognostic scheme for gear health management in surface wear progression of the intelligent manufacturing system | |
Wang et al. | Intelligent diagnosis method for rolling element bearing faults using possibility theory and neural network | |
CN113469060A (zh) | 多传感器融合卷积神经网络航空发动机轴承故障诊断方法 | |
CN113614359A (zh) | 用于预测风力涡轮机齿轮箱组件中的可观察损坏的风险的方法和系统 | |
CN109214355B (zh) | 一种基于核估计lof的机械监测数据异常段检测方法 | |
Elforjani | Diagnosis and prognosis of real world wind turbine gears | |
CN111103137A (zh) | 基于深度神经网络的风力涡轮机齿轮箱故障诊断方法 | |
CN111459144A (zh) | 一种基于深度循环神经网络的飞机飞控系统故障预测方法 | |
CN109580218B (zh) | 一种基于似然学习机的风机齿轮箱状态识别方法 | |
Yang et al. | Gearbox fault diagnosis based on artificial neural network and genetic algorithms | |
CN112461543B (zh) | 一种基于多分类支持向量数据描述的旋转机械故障诊断方法 | |
Hu et al. | Intelligent condition assessment of industry machinery using multiple type of signal from monitoring system | |
CN114942139A (zh) | 齿轮箱中考虑轴承退化影响的齿轮剩余寿命预测方法 | |
Lahrache et al. | Anomaly detection in a cutting tool by k-means clustering and support vector machines | |
CN116296365B (zh) | 一种风力发电机组齿轮箱磨损自动预警方法及系统 | |
CN112380782A (zh) | 一种基于混合指标和神经网络的旋转设备故障预测方法 | |
CN117034001A (zh) | 一种风电机组故障预测方法、装置及电子设备 | |
CN114486262B (zh) | 一种基于cnn-at-lstm的轴承剩余使用寿命预测方法 | |
Wang et al. | Wavelet packet transform-assisted least squares support vector machine for gear wear degree diagnosis | |
CN114118720A (zh) | 一种核电厂循环水泵运行状态综合评价方法、装置及电子设备和存储介质 | |
Latuny | A sensitivity comparison of Neuro-fuzzy feature extraction methods from bearing failure signals | |
Ebrahimi | A method based on support vector machine and vibration analysis for fault detection in bevel gears (Case study: Differential) |
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 |