发明名称 面向数控机床故障诊断与故障预报的知识库构建方法
摘要 本发明涉及一种面向数控机床故障诊断与故障预报的知识库构建方法,其步骤如下:步骤一、通过远程监测设备对高档车削加工中心实时在线监测,获得代表不同故障类型的多组振动数据X<sub>j</sub>(t),j为采集到的振动数据组数,n为正整数;步骤二、对实时在线监测的多组振动数据X<sub>j</sub>(t)依次进行时态粗糙小波包分析处理,得到能量特征向量T′作为条件属性,以故障类型为决策属性,构建故障知识原始决策表;步骤三、对故障知识原始决策表进行基于差别矩阵的故障特征属性约简,生成规则,构成知识库;步骤四、采用规则的可信度作为评价指标对最终规则进行度量和评价。本发明能为故障诊断和故障预报提供有效的保障,本发明可以广泛应用于高档车削加工中心。
申请公布号 CN102736562B 申请公布日期 2014.06.04
申请号 CN201210240271.5 申请日期 2012.07.10
申请人 北京信息科技大学 发明人 徐小力;吴国新;王少红;任彬
分类号 G05B19/406(2006.01)I 主分类号 G05B19/406(2006.01)I
代理机构 北京远大卓悦知识产权代理事务所(普通合伙) 11369 代理人 贺持缓
主权项 1.一种面向数控机床故障诊断与故障预报的知识库构建方法,其步骤如下:步骤一、通过远程监测设备对高档车削加工中心实时在线监测,获得代表不同故障类型的多组振动数据X<sub>j</sub>(t),j=1,2,…,n,j为采集到的振动数据组数,n为正整数;步骤二、对实时在线监测的多组振动数据X<sub>j</sub>(t)依次进行时态粗糙小波包分析处理,得到能量特征向量T′作为条件属性,以故障类型为决策属性,构建故障知识原始决策表;所述时态粗糙小波包分析方法如下:(1)选择正交小波包分解公式<maths num="0001"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&mu;</mi><mrow><mn>2</mn><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>k</mi></munder><msub><mi>h</mi><mrow><mi>k</mi><mo>-</mo><mn>2</mn><mi>x</mi></mrow></msub><msub><mi>&mu;</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>&mu;</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>k</mi></munder><msub><mi>g</mi><mrow><mi>k</mi><mo>-</mo><mn>2</mn><mi>x</mi></mrow></msub><msub><mi>&mu;</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>将采集到的多组振动信号X<sub>j</sub>(t)依次进行小波包分解,得到分解后的所有分解层的低频带系数u<sub>2n</sub>(x)数和高频系数u<sub>2n+1</sub>(x);其中h<sub>k</sub>、g<sub>k</sub>是一对共轭正交实系数滤波器,且h<sub>k</sub>=(-1)<sup>k</sup>g<sub>k</sub>(1-k);k为小波包分解层数,k=1,2,…;n=0、1、2…,其为小波包分解频率带序列号;(2)将小波包分解的最后一层利用粗糙集理论中的等价关系上、下近似的概念,划分成上、下两部分,分别表示粗糙集的上、下近似,此分解层的频率上边界为<img file="FDA0000401010460000014.GIF" wi="162" he="72" />称为上边界时态粗糙小波包低频系数;下边界为<u>u</u><sub>2n+1</sub>(x),称为下边界时态粗糙小波包高频系数;(3)对上边界时态粗糙小波包低频系数<img file="FDA0000401010460000015.GIF" wi="128" he="72" />和下边界时态粗糙小波包高频系数<u>u</u><sub>2n+1</sub>(x)进行重构,则得到重构信号S<sub>n</sub>为:<img file="FDA0000401010460000016.GIF" wi="471" he="72" />(4)分别对各重构信号S<sub>n</sub>求解相应的能量E<sub>nj</sub>=∫|S<sub>nj</sub>(x)|<sup>2</sup>dx,得到能量向量T=[E<sub>n0</sub>,E<sub>n1</sub>,E<sub>n2</sub>,…,E<sub>nj</sub>],j为采集到的振动数据组数;(5)将能量向量T进行归一化处理,得到能量特征向量T′,以能量特征向量T′为条件属性,故障类型为决策属性,构建故障知识原始决策表;步骤三、对故障知识原始决策表进行基于差别矩阵的故障特征属性约简,生成规则,构成知识库;采用所述基于差别矩阵的方法进行决策表属性约简,包括如下步骤:(1)基于所构建的故障知识原始决策表DT=(U,A=T′∪D,V,f),采用动态层次聚类算法对故障知识原始决策表DT进行离散化处理,其中,U={1,2,…,n}是对象集合,即论域的非空有限集合;D是决策属性集合,T′=[E'<sub>n0</sub>,E'<sub>n1</sub>,E'<sub>n2</sub>,…]是条件属性集合,A为包括条件属性和决策属性的全集,V表示信息值域;f={f<sub>α</sub>|f<sub>α</sub>:U→V<sub>α</sub>}表示决策表的信息函数,α=E′<sub>nj</sub>,f<sub>α</sub>为属性α的信息函数,α为条件属性中的一个简单属性;(2)判断离散化后的原始决策表DT的相容性,删除不一致或者重复的数据信息;(3)构建原始决策表DT的差别矩阵:给定T′={E<sub>n</sub>′|E<sub>n</sub>′∈T′}称为条件属性,每个E′<sub>nj</sub>∈T′称为T′的一个简单特征属性,差别矩阵为:<img file="FDA0000401010460000021.GIF" wi="1754" he="307" />式中E<sub>ij</sub>={α|(α∈A)∧(f<sub>α</sub>(x<sub>i</sub>)≠f<sub>α</sub>(x<sub>j</sub>))},<img file="FDA0000401010460000022.GIF" wi="328" he="66" />差别矩阵的元素E<sub>ij</sub>是区分对象E<sub>i</sub>和E<sub>j</sub>的所有简单属性组成的集合,当对象E<sub>i</sub>=E<sub>j</sub>时,<img file="FDA0000401010460000023.GIF" wi="175" he="72" />(4)根据差别矩阵建立相应的析取范式:<img file="FDA0000401010460000024.GIF" wi="669" he="103" />对P<sub>∧(∨)</sub>进行布尔运算,得到一个合取范式:<img file="FDA0000401010460000025.GIF" wi="337" he="103" />其中P<sub>k</sub>=α(x<sub>i</sub>,x<sub>j</sub>);(5)输出属性约简<img file="FDA0000401010460000026.GIF" wi="567" he="76" />即每个合取项P<sub>k</sub>为原始决策表的一个属性约简RED(C);(6)将得到的所有属性约简RED(C)进行属性值约简,删除属性约简中的冗余决策属性,即对每个合取项P<sub>k</sub>中的α计算相对决策属性D的重要度sig(α,P<sub>k</sub>,D):<maths num="0002"><![CDATA[<math><mrow><mi>sig</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><msub><mi>P</mi><mi>k</mi></msub><mo>,</mo><mi>D</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>|</mo><msub><mi>pos</mi><mrow><msub><mi>P</mi><mi>k</mi></msub><mo>&cup;</mo><mi>&alpha;</mi></mrow></msub><mrow><mo>(</mo><mi>D</mi><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>|</mo><msub><mi>pos</mi><msub><mi>P</mi><mi>k</mi></msub></msub><mrow><mo>(</mo><mi>D</mi><mo>)</mo></mrow><mo>|</mo></mrow><mrow><mo>|</mo><mi>U</mi><mo>|</mo></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>取sig(α,P<sub>k</sub>,D)=1的所有属性α构成最终规则,生成故障诊断知识库;步骤四、采用规则的可信度α作为评价指标对最终规则进行度量和评价,并将可信度α≥80%的最终规则形成知识库,并传输至故障诊断与预报知识库内;样本规则库和故障案例库存储的历史数据也传输至故障诊断与预报知识库内;故障诊断与预报知识库将接收到的各个数据综合处理后,输出故障信息,实现为高档车削加工中心进行故障诊断及故障预测;同时将故障信息作为历史数据存储至样本规则库和故障案例库。
地址 100092 北京市海淀区清河小营东路12号166信箱