发明名称 一种多传感器融合的风电机组齿轮箱故障定位方法
摘要 一种多传感器融合的风电机组齿轮箱故障定位方法,它有七大步骤:一、依据齿轮箱内部结构以及振动传感器的位置和数量来确定故障定位模型应用条件;二、依据现有的故障诊断技术对采集的振动信号进行分析,提取相应的故障特征值;三、借鉴节约覆盖集理论建立故障部件与特征值集合的关系;四、根据齿轮箱内窥镜的结果确定齿轮箱故障的部位、类型以及故障的程度;五、分析不同故障部件对多个测点振动信号的影响,得到振动的衰减系数;六、利用故障特征值的差异以及BP神经网络算法对得到的数据进行故障诊断,将故障定位到某类型的部件;七、求取和分析振动信号的衰减系数β'<sub>k</sub>,采用对比的方法,通过求误差δ=β′<sub>k</sub>-β<sub>k</sub>最小的方法来确定故障部件的位置。
申请公布号 CN104048825B 申请公布日期 2016.05.25
申请号 CN201410230862.3 申请日期 2014.05.28
申请人 北京航空航天大学 发明人 钱政;田双蜀;张颖;王现伟
分类号 G01M13/02(2006.01)I 主分类号 G01M13/02(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种多传感器融合的风电机组齿轮箱故障定位方法,其特征在于:该方法的具体实现步骤如下:步骤一:依据齿轮箱的内部结构以及振动传感器的安装位置和数量来确定故障定位模型应用的条件;风电机组齿轮箱均由齿轮、轴承、轴及相应附件组成,由于其传动功率以及生产厂家的不同,齿轮箱内部器件的参数不同,因此故障定位的部件数量及振动特征存在着差别,不同的振动传感器安装位置和安装数量会导致故障特征值的大小及同种部件定位算法的不同;步骤二:依据现有的故障诊断技术包括时域分析、频域分析、小波分析现有的故障诊断方法对采集的振动信号进行分析,提取相应的故障特征值{a<sub>1</sub>,a<sub>2</sub>…a<sub>n</sub>};步骤三:借鉴节约覆盖集理论建立故障部件与特征值集合的关系;不同部件的故障对应不同的特征值集合,借鉴节约覆盖集理论得到最佳的故障与特征值的对应关系;利用节约覆盖集理论分析过程如下:设F代表齿轮箱可能发生的故障部件集合,|F|代表可能的故障数目,则故障发生共有2<sup>|F|</sup>种可能性;每种可能性看作F的一个子集F<sub>I</sub>,若F<sub>I</sub>={f<sub>1</sub>,f<sub>2</sub>,…,f<sub>n</sub>}表示故障部件{f<sub>1</sub>,f<sub>2</sub>,…,f<sub>n</sub>}发生故障,而其他f<img file="FDA0000512329340000011.GIF" wi="140" he="71" />的故障则没有发生;解决齿轮箱中多故障诊断问题就是在所有2<sup>|F|</sup>个可能性中寻找最可能的组合,利用符号推理—节约覆盖理论来寻找最佳的组合;节约覆盖理论是一种利用诱导推理来阐明诊断专家系统理论,它将诊断问题描述为一个四元组P=&lt;F,A,R,A<sup>+</sup>&gt;,其中F={f<sub>1</sub>,f<sub>2</sub>,…,f<sub>n</sub>}表示故障部件的有限非空集;A={a<sub>1</sub>,a<sub>2</sub>,…,a<sub>n</sub>}表示故障特征值的有限非空集;<img file="FDA0000512329340000021.GIF" wi="214" he="60" />表示定义在F×A上的有序关系子集;<img file="FDA0000512329340000022.GIF" wi="166" he="62" />表示已知的特征值集合;符号R代表故障部件和故障特征值之间的直接因果关系,&lt;f<sub>i</sub>,a<sub>j</sub>&gt;∈R表示f<sub>i</sub>能够引起a<sub>j</sub>,它并不表示当f<sub>i</sub>存在,a<sub>j</sub>总是发生,而仅仅是可能发生;A<sup>+</sup>是A的一个特殊子集,代表了在一个特定问题下已知存在的特征,不在A<sup>+</sup>的特征可认为是不存在的;此外,可定义两个函数:对于所有a<sub>j</sub>∈A,parts(a<sub>j</sub>)={f<sub>i</sub>|&lt;f<sub>i</sub>,a<sub>j</sub>&gt;∈R}代表特征值a<sub>j</sub>所有可能引起的故障部件,features(f<sub>i</sub>)={a<sub>j</sub>|&lt;f<sub>i</sub>,a<sub>j</sub>&gt;∈R}代表可能被f<sub>i</sub>引起的所有特征,同时定义<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>parts</mi><mrow><mo>(</mo><msub><mi>A</mi><mi>J</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mo>&cup;</mo><mrow><msub><mi>a</mi><mi>j</mi></msub><mo>&Element;</mo><msub><mi>A</mi><mi>J</mi></msub></mrow></msub><mi>parts</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000512329340000023.GIF" wi="577" he="87" /></maths>和<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>features</mi><mrow><mo>(</mo><msub><mi>F</mi><mi>I</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mo>&cup;</mo><mrow><msub><mi>f</mi><mi>i</mi></msub><mo>&Element;</mo><msub><mi>F</mi><mi>I</mi></msub></mrow></msub><mi>features</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000512329340000024.GIF" wi="725" he="78" /></maths>对于诊断问题,parts(a<sub>1</sub>)={f<sub>1</sub>,f<sub>2</sub>,f<sub>3</sub>},features(f<sub>1</sub>)={a<sub>1</sub>,a<sub>3</sub>},parts({a<sub>3</sub>,a<sub>4</sub>})={f<sub>1</sub>,f<sub>3</sub>,f<sub>4</sub>},如果<img file="FDA0000512329340000025.GIF" wi="408" he="71" />则称故障部件集合F<sub>I</sub>是特征值<img file="FDA0000512329340000026.GIF" wi="157" he="66" />的一个覆盖;利用最小准则确定齿轮箱中每个部件故障对应最少特征值覆盖;步骤四:根据现场齿轮箱内窥镜的结果确定齿轮箱故障的部位、类型以及故障的程度;利用与内窥镜检查相近时间内齿轮箱的振动数据来分析单一故障的不同程度对故障特征值的影响以及不同部位故障对不同振动测点特征值的影响;步骤五:分析不同故障部件对多个测点振动信号的影响,得到振动的衰减系数;故障部件的振动信号传播到不同振动测点的衰减率不同,进而导致同一故障部件振动信号在不同测点的特征值不同;以其中的某个测点作为基准,其它测点的值表示为<img file="FDA0000512329340000027.GIF" wi="342" he="79" />其中i为第i次测量的结果,k为第k个故障部件,m为第m个振动测点;对多次测量的数据进行结果平均得到第k故障部件振动信号传输到测点的衰减系数{β<sub>k,1</sub>,β<sub>k,2</sub>,β<sub>k,3</sub>,…β<sub>k,m</sub>},其中<img file="FDA0000512329340000031.GIF" wi="348" he="135" />其中l为重复测量的次数,则多个部件的衰减系数表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&beta;</mi><mo>=</mo><msub><mfenced open='{' close='}'><mtable><mtr><mtd><msub><mi>&beta;</mi><mn>1,1</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>&beta;</mi><mrow><mn>1</mn><mo>,</mo><mi>M</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>&beta;</mi><mrow><mi>K</mi><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>&beta;</mi><mrow><mi>K</mi><mo>,</mo><mi>M</mi></mrow></msub></mtd></mtr></mtable></mfenced><mrow><mi>K</mi><mo>&times;</mo><mi>M</mi></mrow></msub></mrow>]]></math><img file="FDA0000512329340000032.GIF" wi="542" he="245" /></maths>其中,K为齿轮箱中部件的个数,M为振动测点数,β<sub>k</sub>={β<sub>k,1</sub>…β<sub>kM</sub>};步骤六:不同类型的部件具有相异的故障特征值,利用特征值的差异以及BP神经网络算法对得到的数据进行故障诊断,将故障定位到某类型的部件;步骤七:对于同种类型的部件,求取振动信号的衰减系数β'<sub>k</sub>,分析衰减系数,采用对比的方法,通过求误差δ=β′<sub>k</sub>‑β<sub>k</sub>最小的方法来确定故障部件的位置。
地址 100191 北京市海淀区学院路37号