发明名称 风力发电机组齿轮箱的故障诊断方法
摘要 本发明涉及高故障设备的故障诊断技术领域,公开了一种风力发电机组齿轮箱的故障诊断方法,包括第一步:选择齿轮箱的若干个有效特征建立多维目标空间,模拟若干个故障类别,并对所述若干个故障类别的故障数据进行采集;第二步:对所述故障数据进行时域和幅域分析,提取频域参数和幅域参数成为训练数据集;第三步:将所述训练数据集看成一个种群,所述训练数据集中的每一个数据看成一个粒子,通过迭代算法优化数据子集;第四步:通过适应度函数f(i),更新粒子的速度和位置;第五步:得到最优子集。本发明解决了风电机组齿轮箱故障诊断中的数据不均衡问题,PSOEE算法提高了分类器对不均衡数据集的预报精度。
申请公布号 CN103278326A 申请公布日期 2013.09.04
申请号 CN201310234926.2 申请日期 2013.06.14
申请人 上海电机学院 发明人 刘天羽;邢飞
分类号 G01M13/02(2006.01)I 主分类号 G01M13/02(2006.01)I
代理机构 上海翼胜专利商标事务所(普通合伙) 31218 代理人 翟羽;黄燕石
主权项 1.一种风力发电机组齿轮箱的故障诊断方法,其特征在于:包括如下方法:第一步:选择齿轮箱的若干个有效特征建立多维目标空间,模拟若干个故障类别,并对所述若干个故障类别的故障数据进行采集;第二步:对所述故障数据进行时域和幅域分析,提取频域参数和幅域参数成为训练数据集;第三步:将所述训练数据集看成一个种群,所述训练数据集中的每一个数据看成一个粒子,通过下述迭代算法优化数据子集:v<sub>id</sub>(k+1)=v<sub>id</sub>(k)+c<sub>1</sub>r<sub>1</sub>(p<sub>id</sub>(k)-x<sub>id</sub>(k))+c<sub>2</sub>r<sub>2</sub>(p<sub>gd</sub>(k)-x<sub>id</sub>(k))x<sub>id</sub>(k+1)=x<sub>id</sub>(k)+v<sub>id</sub>(k+1)(i=1,2,…,m;d=1,2,…,D)其中第i个粒子表示成一个D维的向量<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>X</mi><mo>&RightArrow;</mo></mover><mi>i</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mi>iD</mi></msub><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>m</mi><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>i</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>v</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>v</mi><mi>iD</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>是第i个粒子的飞行速度,<img file="FDA00003343884300013.GIF" wi="439" he="90" />是第i个粒子的最优位置pbest,<img file="FDA00003343884300014.GIF" wi="478" he="92" />是整个种群历史搜索到的最优位置gbest,其中D为所述目标空间的维数,m粒子的总数,在式中,k是迭代代数;学习因子c<sub>1</sub>和c<sub>2</sub>为非负常数;r<sub>1</sub>和r<sub>2</sub>是均匀分布于[0,1]之间的两个随机任意数;v<sub>id</sub>∈[-v<sub>max</sub>,v<sub>max</sub>],v<sub>max</sub>是预先设定的常数;v<sub>id</sub>(k)为第k次迭代时粒子i飞行速度矢量的第d维的分量;x<sub>id</sub>(k)为第k次迭代时粒子i位置矢量的第d维的分量;p<sub>id</sub>(k)是粒子i个体最好位置的第d维分量;p<sub>gd</sub>(k)为种群最好位置的第d维分量;迭代终止条件设定为最大迭代数或(和)粒子群迄今为止搜索到的最优位置满足最小适应度值。第四步:通过适应度函数f(i),更新粒子的速度和位置,适应度函数如下:<maths num="0003"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mi>a</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mi>p</mi><mfrac><msub><mi>m</mi><mi>ic</mi></msub><msub><mi>m</mi><mi>all</mi></msub></mfrac></mrow></math>]]></maths>式中a(i)为粒子i所选的特征子集的分类正确率估计;p为调整参数,是分类准确率和所选特征数量的平衡系数;M<sub>ic</sub>为粒子i所选的特征子集中特征的数量;M<sub>all</sub>为所有待选特征的总数。第五步:得到最优子集,将所述最优子集作为第三步中的种群,重复第三步,直至得到集成模型。
地址 200240 上海市闵行区江川路690号