发明名称 一种智能刀具故障诊断方法
摘要 本发明涉及一种智能刀具故障诊断方法,其步骤:在机床的CNC端与机床端之间设置PLC;由PLC获取机床加工过程中产生的若干刀具振动信号x(t),并对采集到的各刀具振动信号x(t)进行处理并提取刀具特征;根据获得的刀具特征信息建立刀具振动信号和刀具磨损故障的对应关系,以完成刀具磨损故障的智能识别:对每个刀具振动信号进行训练得到隐马尔可夫模型λ;利用训练好的隐马尔可夫模型λ进行似然率的计算,将刀具振动信号的序列O输入到以训练好的各个隐马尔可夫模型中,得到各刀具振动信号在相应隐马尔可夫模型下的似然率,进而实现刀具状态识别。本发明结合机床本身特性和刀具加工过程振动信号特征,实现刀具磨损故障的自动检测,并能自动识别刀具磨损状态。
申请公布号 CN106002483A 申请公布日期 2016.10.12
申请号 CN201610289610.7 申请日期 2016.05.04
申请人 北京信息科技大学 发明人 黄民;高延;吴国新;孙巍伟;马超
分类号 B23Q17/09(2006.01)I 主分类号 B23Q17/09(2006.01)I
代理机构 北京远创理想知识产权代理事务所(普通合伙) 11513 代理人 张素妍
主权项 一种智能刀具故障诊断方法,其特征在于,该方法包括以下步骤:1)在机床的CNC端与机床端之间设置PLC,由PLC实现对机床侧和CNC侧的输入输出信号处理;2)由PLC获取机床加工过程中产生的若干刀具振动信号x(t),并对采集到的各刀具振动信号x(t)进行处理并提取刀具特征:2.1)对采集到的刀具振动信号x(t)进行经验模态分解,将刀具振动信号x(t)中不同尺度的波动或趋势逐级分解出来,产生具有不同特征尺度的数据序列,并将这些序列定义为本征模态函数IMF;2.2)设置IMF分量的两个判定条件为:(1)IMF的极值点和过零点的个数差值不能大于一;(2)IMF的上包络线和下包络相对于时间轴是局部对称的;若同时满足上述两个条件,则为IMF分量;2.3)对IMF分量进行希尔伯特‑黄变换,利用IMF分量c<sub>i</sub>(t)构造时间复信号z<sub>i</sub>(t):<maths num="0001"><math><![CDATA[<mrow><msub><mi>z</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>c</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>j</mi><msub><mover><mi>c</mi><mo>&OverBar;</mo></mover><mi>l</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><msub><mi>j&phi;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></msup></mrow>]]></math><img file="FDA0000981386650000011.GIF" wi="662" he="79" /></maths>式中,φ<sub>i</sub>(t)为相角;a<sub>i</sub>(t)为幅值;i=1、2、…、n;使时间复信号z<sub>i</sub>(t)的实部c<sub>i</sub>(t)和虚部c<sub>l</sub>(t)存在如下关系:<maths num="0002"><math><![CDATA[<mrow><msub><mover><mi>c</mi><mo>&OverBar;</mo></mover><mi>l</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>&pi;</mi></mfrac><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mi>&infin;</mi></mrow><mrow><mo>+</mo><mi>&infin;</mi></mrow></msubsup><mfrac><mrow><msub><mi>c</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>t</mi><mo>-</mo><mi>&tau;</mi></mrow></mfrac><mi>d</mi><mi>&tau;</mi><mo>;</mo></mrow>]]></math><img file="FDA0000981386650000012.GIF" wi="461" he="127" /></maths>2.4)对时间复信号的幅值进行积分运算求取希尔伯特谱H(ω,t):<maths num="0003"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>Re</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><mo>&Integral;</mo><msub><mi>&phi;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></msup><mo>;</mo></mrow>]]></math><img file="FDA0000981386650000013.GIF" wi="598" he="127" /></maths>2.5)通过上述步骤,提取得到刀具振动信号x(t)的有效特征量;3)根据获得的刀具特征信息建立刀具振动信号和刀具磨损故障的对应关系,以完成刀具磨损故障的智能识别:3.1)对每个刀具振动信号进行训练得到隐马尔可夫模型λ,λ=(π,A,B),π为初始状态概率分布的矢量值;A为状态转移概率矩阵,各元素为a<sub>ij</sub>;B为观察值概率矩阵,各元素为b<sub>ij</sub>:3.2)利用训练好的隐马尔可夫模型λ进行似然率的计算,将刀具振动信号的序列O输入到以训练好的各个隐马尔可夫模型中,得到各刀具振动信号在相应隐马尔可夫模型下的似然率,进而实现刀具状态识别。
地址 100192 北京市海淀区清河小营东路12号166信箱