发明名称 一种数控机床切削加工的刀具磨损状态识别方法
摘要 本发明公开了一种数控机床切削加工的刀具磨损状态识别方法,具体包括如下步骤:(1)数据采集(2)采集数据区间化,获得区间数据集<img file="DDA00001669935800011.GIF" wi="787" he="89" />(3)对采集信号进行小波包分解,获得各频段的能量百分比和相应各频段的上界小波包系数和下界小波包系数;(4)获取多观察序列(5)求取广义隐马尔科夫初始模型λ=(A,B,π);(6)模型训练获得最优模型库λ=(λ<sub>1</sub>,...,λ<sub>n</sub>);(7)把待识别磨损状态的刀具磨损信号数据作为多观察输入代入上述最优模型库,识别刀具磨损状态。本发明通过广义区间概率来解决数据机床加工中研究刀具磨损状态时,信息处理中出现的不确定性问题,识别准确率显著提高。
申请公布号 CN103105820A 申请公布日期 2013.05.15
申请号 CN201210160153.3 申请日期 2012.05.22
申请人 华中科技大学 发明人 谢锋云;胡友民;吴波;贾广飞;金超;王小岑;程瑶
分类号 G05B19/4065(2006.01)I 主分类号 G05B19/4065(2006.01)I
代理机构 华中科技大学专利中心 42201 代理人 李佑宏
主权项 1.一种数控机床切削加工的刀具磨损状态识别方法,具体包括如下步骤:(1)数据采集测量获取数控机床切削加工中刀具在各种磨损状态下的测量数据,形成刀具磨损信号数据集X={x<sub>1</sub>,x<sub>2</sub>,...,x<sub>k</sub>},其中k=1,2,...,m,m为测量数据的个数;(2)采集数据区间化将所述数据集X中的每个元素转换成区间形式,即<img file="FDA00001669935500011.GIF" wi="286" he="61" />则数据集X转换为区间数据集<maths num="0001"><![CDATA[<math><mrow><msup><mi>X</mi><mo>&prime;</mo></msup><mo>=</mo><mo>{</mo><mo>[</mo><msub><munder><mi>x</mi><mo>&OverBar;</mo></munder><mn>1</mn></msub><mo>,</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>]</mo><mo>,</mo><mo>[</mo><msub><munder><mi>x</mi><mo>&OverBar;</mo></munder><mn>2</mn></msub><mo>,</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>]</mo><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mo>[</mo><msub><munder><mi>x</mi><mo>&OverBar;</mo></munder><mi>k</mi></msub><mo>,</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub><mo>]</mo><mo>}</mo><mo>;</mo></mrow></math>]]></maths>(3)对采集信号进行小波包分解,具体为:对所述数据集X={x<sub>1</sub>,x<sub>2</sub>,...,x<sub>k</sub>}进行小波包分解,获得各频段的能量百分比;对区间数据集X′进行小波包分解,即分别对X′各元素的上下界分别进行小波包分解,得到相应各频段的上界小波包系数和下界小波包系数;(4)获取多观察序列(4.1对数据集X进行小波包分解后得到各频段能量百分比进行区间化,然后将区间化了的各频段的能量百分比的上下界分别进行量化和编码,得到编码后的能量百分比的上界序列<u>O</u><sub>1</sub>和下界序列<img file="FDA00001669935500013.GIF" wi="97" he="58" />(4.2)对各频段的上界小波包系数和下界小波包系数,求得相应的上界小波包系数的均方根值和下界小波包系数的均方根值,从而得到上界的均方根值序列<u>O</u><sub>2</sub>和下界的均方根值序列<img file="FDA00001669935500014.GIF" wi="102" he="58" />(4.3)组合步骤(4.1)和(4.2)的序列,得到多观察序列<maths num="0002"><![CDATA[<math><mrow><mi>O</mi><mo>=</mo><mo>[</mo><msub><munder><mi>O</mi><mo>&OverBar;</mo></munder><mn>1</mn></msub><mo>,</mo><msub><mover><mi>O</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>,</mo><msub><munder><mi>O</mi><mo>&OverBar;</mo></munder><mn>2</mn></msub><mo>,</mo><msub><mover><mi>O</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>]</mo><mo>;</mo></mrow></math>]]></maths>(5)模型训练将所述多观察序列<img file="FDA00001669935500022.GIF" wi="411" he="66" />作为输入,代入数控机床的广义隐马尔科夫初始模型中,逐步训练改进初始模型参数,直至多观察序列在模型中的区间概率p(O|λ)的上、下界都收敛为止,即可得到刀具在各个磨损状态相应的最优广义隐马尔科夫模型<img file="FDA00001669935500023.GIF" wi="392" he="75" />并构成最优模型库<img file="FDA00001669935500024.GIF" wi="383" he="75" />其中,i为最优广义隐马尔科夫模型序号,代表刀具磨损状态类型,且i=1,2,...,n,n为刀具磨损状态数;(6)刀具磨损状态识别把待识别磨损状态的刀具的磨损信号数据,经过步骤(2)~(4)的处理,得到多观察序列<img file="FDA00001669935500025.GIF" wi="495" he="79" />将其划分为上界观测序列<u>O</u>′=[<u>O</u>′<sub>1</sub>,<u>O</u>′<sub>2</sub>]和下界观测序列<img file="FDA00001669935500026.GIF" wi="333" he="74" />分别作为多观察输入,代入所述最优模型库<img file="FDA00001669935500027.GIF" wi="342" he="75" />中,分别计算得到区间概率<img file="FDA00001669935500028.GIF" wi="205" he="108" />和<img file="FDA00001669935500029.GIF" wi="252" he="118" />进一步得到区间概率的对数似然值<img file="FDA000016699355000210.GIF" wi="322" he="108" />其中<maths num="0003"><![CDATA[<math><mrow><mi>p</mi><mrow><mo>(</mo><msup><mi>O</mi><mo>&prime;</mo></msup><mo>|</mo><msub><mover><mi>&lambda;</mi><mo>~</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mo>:</mo><mo>=</mo><mi>p</mi><mrow><mo>(</mo><msup><munder><mi>O</mi><mo>&OverBar;</mo></munder><mo>&prime;</mo></msup><mo>|</mo><msub><mover><mi>&lambda;</mi><mo>~</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mo>&times;</mo><mi>p</mi><mrow><mo>(</mo><msup><mover><mi>O</mi><mo>&OverBar;</mo></mover><mo>&prime;</mo></msup><mo>|</mo><msub><mover><mi>&lambda;</mi><mo>~</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>由区间比较法求出最大的<img file="FDA000016699355000212.GIF" wi="314" he="108" />此时,i对应的状态即为待识别样本的刀具磨损状态。
地址 430074 湖北省武汉市洪山区珞瑜路1037号