[发明专利]一种风力发电机组振动故障诊断方法在审

专利信息
申请号: 201611071559.9 申请日: 2016-11-29
公开(公告)号: CN106596116A 公开(公告)日: 2017-04-26
发明(设计)人: 贾嵘;李涛涛;武桦;张惠智;尹浩霖 申请(专利权)人: 西安理工大学
主分类号: G01M15/12 分类号: G01M15/12;G06K9/62;G06N3/12
代理公司: 西安弘理专利事务所61214 代理人: 李娜
地址: 710048*** 国省代码: 陕西;61
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 风力 发电 机组 振动 故障诊断 方法
【权利要求书】:

1.一种风力发电机组振动故障诊断方法,其特征在于,包括以下步骤:

步骤1,采用振动传感器采集风力发电机组振动故障信号;

步骤2,利用变分模态分解方法对振动信号进行分解,得到包含不同频段信息的IMF分量;

步骤3,计算振动信号经变分模态分解得到的各个IMF分量的能量,并经归一化处理,得到该信号的故障特征向量;

步骤4,利用遗传算法优化支持向量机,构建基于遗传算法优化的支持向量机故障诊断模型;

步骤5,将提取的故障特征向量输入到步骤4建立的诊断模型中进行故障的识别与诊断。

2.根据权利要求1所述的风力发电机组振动故障诊断方法,其特征在于,所述步骤2具体包括以下子步骤:

步骤2.1,估计IMF分量频率带宽的目标,

(1)对每个模态函数uk(t)进行Hibert变换,得到每个模态函数uk(t)的解析信号,如式(1)所示:

<mrow><mo>&lsqb;</mo><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mi>j</mi><mrow><mi>&pi;</mi><mi>t</mi></mrow></mfrac><mo>&rsqb;</mo><mo>*</mo><msub><mi>u</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>

(2)利用指数修正,将每个模态函数的频谱调制到各自估算的中心频率,如式(2)所示:

<mrow><mo>&lsqb;</mo><mo>(</mo><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mi>j</mi><mi>n</mi></mfrac><mo>)</mo><mo>*</mo><msub><mi>u</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>jw</mi><mi>k</mi></msub><mi>t</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>

(3)计算以上解调信号的梯度的平方L2范数,估算出各模态函数的带宽,得到其对应的约束变分问题:

<mrow><munder><mi>min</mi><mrow><mrow><mo>{</mo><msub><mi>u</mi><mi>k</mi></msub><mo>}</mo></mrow><mo>,</mo><mrow><mo>{</mo><msub><mi>w</mi><mi>k</mi></msub><mo>}</mo></mrow></mrow></munder><mo>{</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mo>|</mo><mo>|</mo><msub><mo>&part;</mo><mi>t</mi></msub><mo>{</mo><mo>(</mo><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mi>j</mi><mrow><mi>&pi;</mi><mi>t</mi></mrow></mfrac><mo>)</mo><mo>*</mo><msub><mi>u</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>}</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>jw</mi><mi>k</mi></msub><mi>t</mi></mrow></msup><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>}</mo></mrow>

<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>u</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi></mrow></mtd></mtr></mtable></mfenced>

式中,uk={u1,u2,…uK}为各模态函数集;wk={w1,w2,…wK}为各模态中心频率;为对函数求时间t的偏导数;δt为单位脉冲函数;j为虚数单位;*表示卷积;

步骤2.2,将上述约束变分问题转化为非约束变分问题,引入增广拉格朗日函数L,如下式所示:

<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>L</mi><mrow><mo>(</mo><mo>{</mo><msub><mi>u</mi><mi>k</mi></msub><mo>}</mo><mo>,</mo><mo>{</mo><msub><mi>w</mi><mi>k</mi></msub><mo>}</mo><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>=</mo><mi>&alpha;</mi><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mo>|</mo><mo>|</mo><msub><mo>&part;</mo><mi>t</mi></msub><mo>{</mo><mo>(</mo><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mi>j</mi><mrow><mi>&pi;</mi><mi>t</mi></mrow></mfrac><mo>)</mo><mo>*</mo><msub><mi>u</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>}</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>jw</mi><mi>k</mi></msub><mi>t</mi></mrow></msup><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mo>|</mo><mo>|</mo><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>u</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><mo>&lt;</mo><mi>&lambda;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>u</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&gt;</mo></mrow></mtd></mtr></mtable></mfenced>

式中,α为带宽参数;λ(t)为拉格朗日乘子;

步骤2.3,采用交替方向乘子算法求取式(2)扩展的拉格朗日函数,得到K个IMF分量。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西安理工大学,未经西安理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201611071559.9/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top