欢迎您访问长沙鹏翔电子科技有限公司官方网站

技术与应用

PCIE高速声发射仪/千兆网络声发射仪

技术与应用

基于EEMD和IMF能量分布的刀具破损识别术

发布日期:2016-01-19 13:44    浏览次数:

 杨明伦,邵华

(上海交通大学机械与动力工程学院,上海200240)
摘要:针对铣削过程中声发射信号非平稳的特点,提出了一种基于噪声辅助经验模态分解(EEMD)和本征模函数(IMF)能量分布的刀具破损识别方法。首先对经过滤波后的原信号进行EEMD分解,抽取本征模函数组(IMF),后计算每一阶模函数能量及总体能量分布,最后提取特征向量,通过特征向量的变化识别刀具破损。利用该方法,在立式铣削加工中心上对稳定切削中刀具破损和变参数铣削加工进行了系统的分析,结果表明此方法能够剔除切削参数变化的影响,准确的识别刀具破损,具有很高的稳定性和准确性。
0 引言
刀具状态的实时监控对提高生产效率,降低生产成本,生产自动化以及产品质量保证起到重要的作用。切削刀具破损作为一种主要的刀具失效形式,一直是各国学者研究的热点。目前在刀具破损监控上研究主要围绕着两方面,一是传感器技术,二是信号处理技术。多年来,随着传感器技术的不断发展,不同类型的传感器被应用在刀具状态监控上面‘1|,如力传感器‘川,功率变送器[3],声发射H1传感器等。在众多传感器中,声发射传感器灵敏性高,抗干扰性强以及经济性等特点”。,逐渐被各国研究人员所采用。多种信号处理与特征提取技术被开发并应用到监控系统中,从早期的时域特征分析∞1,到中期的基于傅里叶变换的时.频域特征分析"1,再到近年发展起来的小波分析"1,神经网络等方法,信号处理技术的发展对刀具破损监控起着推动作用。但是,针对铣削过程中刀具破损的监控一直没有很好的解决方法,原因主要在于两方面:首先,在铣削加工中,切削参数变化会对信号处理技术的稳定性产生影响;其次,铣削过程中的AE信号表现出很强的非线性和非平稳性,这种非平稳性会导致AE信号频域特征随着时间的推移而改变,无法用固定的频域特征描述铣削AE信号。经验模态分解(empirical mode decomposition,EMD)是一种基于瞬时频率的自适应的处理非线性系统和非平稳信号的处理方法"1。原信号被EMD分解成一组分量信号,这些分量信号相互正交,被定义为本征模函数(IMF)。这些模函数能够内在的描述原信号的特征,能够对原信号起到识别的作用。
而对于铣削AE信号而言,传统的经验模态分解最大的弊端就是模态混叠问题。针对模态混叠问题,噪声辅助经验模态分解(EEMD)被提出’1…,基本消除了模态混叠的影响。
本文首先利用EEMD对铣削AE信号和刀具破损信号进行分解,得到本征模函数(IMF),计算IMF的能量分布,建立信号类型与特征向量的对应关系。进一步通过系统的试验,验证了基于EEMD和IMF能量分布的破损识别方法对切削参数变化有效性。
1 基于EEMD与IMF能量分布的识别方法
1.1噪声辅助经验模态分解(EEMD)
经验模态分解(EMD)是一个自适应的筛选过程。经过END分解之后,原信号被分解成一些类的分量信号,定义为本征模函数(IMF),满足两条标准:①在整个数据长度内,极值点(包括极大值和极小值)数目与过零点的数目相等或者最多相差1;②在信号长度内的任意点,上包络线与下包络线的均值必须等于零。EMD分解出INFs的具体过程可以简要的描述为:
(1)对于原信号茗(t),信号长度为Ⅳ,用样条曲线分别链接局部最大值和局部最小值,形成上包络线M(t)和下包络线z(t)。上下包络线把原信号包含其中,可以得到局部平均值m.(t)和与原信号的差值hl(t):
(2)如果h.(t)满足本征模函数的两条定义,那么h.(t)就是第一阶IMF。否则,用h,(t)替代原信号菇(t),重新进行上述筛选过程,计算h.(t)的上下包络线(“。(t),f。(t))并计算局部平均值与h。(t)的差值:
(3)当相邻的两个筛选结果的标准差(SD)降低到设定好的阈值之下时,可以获得一个IMF,通用的标准差阈值为0.2~0.3。通常情况下,为了获得一阶IMF需要重新进行||2次筛选过程:
(4)在第一阶IMF被分解出之后,把原信号石(£)与c:(t)的差值作为新的原信号,按照相同的标准重复筛选过程直至得到全部的IMF。当第Ⅳ阶IgFc。(t)或者残余分量r。(t)下降超过设定好的阈值时,或者是残余分量函数r。(#)变成单调函数时,筛选过程停止。
尽管END分解是一种处理非平稳信号有效的工具,但是传统的算法仍然存在一定的缺陷。而对于铣削AE信号而言,模态混叠产生的影响最为严重。模态混叠的定义为:一阶INF包含多个时间尺度,或者一个时间尺度分布在不同的INF中。为了克服这个缺陷,噪声辅助经验模态分解(EEMO)被建立以消除模态混叠现象。
EEMD原理简单:①把白噪声混入原信号(混入白噪声的标准差与原信号标准差的比例为0.25);②对混合后信号应用EMD分解;③应用不同的自噪声序列不断的进行步骤①②;④对每一阶IMF,对所有分解过程分解出的结果取平均值,作为最终结果。
1.2 IMF能量分布
EEMD等价于一个自适应的带通滤波器,原信号戈(t)分解成一组分量信号,
每~个IMF imf,(t)是原信号在特定频率段上的分量,而残余分量r。(t)的能量很小,所以在计算信号总能量时可忽略不计。原信号的总能量可以表示为:
所有的IMFs都是相互近似正交的:
计算原信号在每一个自适应频段上的能量分布,进而定义信号的频率.能量分布状态向量:
当切削条件发生改变时,比如刀具突然破损、切削参数改变或者是其他因素的影响,IMF能量分布和状态向量囝将会同时产生一些改变。为了定量的描述这种切削状态,△g(t)and艿…(t)被引入来表示切削过程的变化:
2 识别方法的应用
2.1 试验系统与试验设计
试验系统如图2所示,铣削试验在DMG公司生产的DMU70五轴立式加工中上进行,传感器采用富士陶瓷公司生产的FUJICERA 1045s宽带式声发射传感器,前置放大器放大倍率为40DB,工件为45#钢,尺寸为50×50×400mm。传感器涂抹耦合剂之后,通过磁吸附装置固定在工件末端,如图2所示。传感器输出信号经过前置放大器之后,被送进数据采集卡(ADLINK PCI9816),采样频率5NHZ,后由基于NI—LABVIEW平台的数据采集保存程序处理,保存成二进制格式。被保存下来的原信号,经过20kHz高通滤波后,送人基于MATLAB平台的数据处理程序进行后处理。
为了研究铣削参数改变对识别方法的影响,在图2所示的试验系统上进行铣削试验,试验条件如表1。为了获取稳定切削过程中刀具突然破损的信号,将直径为0.5mm的硬质合金钻头镶嵌人45#钢基体内部,形成硬质点,当切削刀具接触硬质点时发生破损,如图2所示,试验条件如表2。
2.2 基于IMF能量分布的破损信号特征识别
经验模态分解相当于一个自适应的滤波器组,原信号经过EEMD分解之后,形成多阶本征模函数,每一阶模函数是原信号在特定频段范围内的分量。由于各分量信号之间相互正交,且每一阶IMF的能量都来自原信号,所以所有分量信号的能量总和等于原信号的能量¨…,各本征模函数的能量分布状态能够自适应的反映原信号的频率一能量状态。铣削AE信号具有很强的非平稳性,利用IMF能量分布能够有效地揭示原信号内在的特征,提取信号的特征向量。
图1为稳定切削中刀具突然破损的信号时域波形图。从波形图中,可以看出稳定切削过程中刀具破损信号分为三部分:稳定切削,刀具破损以及破损刀具切削。新刀具切削工件产生的AE信号平稳性高,突发型的信号少,信号幅值较低;刀具发生破损时,刀具几何形状改变的断裂能瞬间释放,产生能量很大,幅值很高的信号,信号平稳性差,衰减非常快;破损后的刀具切削工件产生的AE信号稳定性较新刀具差,幅值较高并伴随着随即突发型信号。这是由刀刃几何形状发生改变后,表面几何形貌不规则,刀具与工件接触面积增大所造成的。
 
分另0取0.008—0.018s,0.0244—0.0264s,0.035—0.045s三个时间片段进行分析。图4表示的是稳定切削信号,刀具破损信号以及破损刀具切削信号的功率谱。从功率谱中可以看出,稳定切削和破损刀具切削AE信号的主要能量都集中在30~
50kHz区域,而且相对于稳定切削,破损刀具切削信号的能量大,高频成分含量高。刀具破损信号的主要能量集中在100~400kHz区域内,能量远高于铣削信号。但是,仅仅从能量大小上来识别刀具破损信号是不准确的,因为切削参数的改变会剧烈的影响AE信号的能量,所以需要对图4所示的特征进行数学建模描述。图5为对三段信号分别进行EEMD分解得到的IMFs分量信号,从图中可以看出,刀具破损信号的主要能量集中在第1、2、3阶IMF上,而稳定切削信号主要能量集中在第3、4阶IMF。对于破损刀具切削,第1、2阶IMF主要分离出的是随机的突发型信号,但是主要能量仍然集中在第4阶IMF上。应用上一节提出的IMF能量分布算法,计算各IMF能量,得到三段信号的IMF能量百分比,数值见表3。在此基础上绘制出IMF能量分布曲线,如图6。

如图6所示,对于稳定切削和破损刀具切削信号,IMF能量分布曲线峰值在第4阶IMF处,信号的主要能量都集中在第四阶IMF上面,而刀具破损信号的曲线峰值出现在第3阶IMF处。主要能量集中在第三阶IMF上。从识别的角度,当刀具发生破损时,IMF分布曲线的峰值从第4阶移动到第3阶。通过提取IMF能量分布曲线的峰值位置,可以准确地从铣削信号中识别出刀具破损信号,实现刀具破损的监控。
2.3变切削参数的影响
在铣削过程中切削参数经常发生变化,而切削参数的变化会对AE信号产生明显的影响,所以识别方法必须对切削参数的变化不敏感,才能在实际应用中得到推广。针对切削速度,切削深度和进给率三个参数,分别进行单因素试验,对每一种切削参数下采集的信号提取IMF能量分布,绘制曲线,如图7、8、9。从图7中可以得到结论,切削速度不会改变IMF能量分布状态,也不会改变分布曲线的峰值位置,从图8、9中,也可以得到相同的结论,综上所述,切削参数的变化不会使IMF能量分布曲线峰值位置产生移动,只有刀具破损才能使其从第4阶移动到第3阶。可以通过监测曲线峰值位置来识别刀具破损信号,大量实验表明,识别方法能够剔除切削参数的干扰,准确的识别刀具破损。

 
3 结论
根据铣削过程切削参数改变以及铣削AE信号非平稳的特点,本文提出一种新的基于EEMD和IMF能量分布的刀具破损识别方法,该方法能够自适应的把原信号分解成本征模函数,并通过计算]MF的能量,刻画原信号的频率-能量状态。对提取的IMF能量进行归一化处理,绘制分布曲线,通过判断峰值位置来识别刀具破损的发生。进一步地,通过大量的实验,验证了该识别方法对变切削参数切削过程的有效性。综上所述,本文的主要结论为:(1)铣削信号的主要能量集中在30~50kHz范围内,刀具破损信号主要能量集中在100~400kHz范围内。
(2)稳定切削信号幅值低,突发型信号少,高频成分比重低;破损刀具切削信号幅值较高,并伴随有低幅值突发型信号,信号稳定性差,且高频成分比重较稳定切削信号有所升高;刀具破损信号能量最大,高频成分比重远高于铣削信号。
(3)铣削信号IMF分布曲线峰值位于第4阶,刀具破损信号曲线峰值位于第3阶段。
(4)切削速度,切削深度和进给率不会使分布曲线峰值移动,只有刀具破损能够对IMF能量分布状态产生影响。
(5)在铣削过程中,通过监测AE信号IMF能量分布曲线峰值的移动,可以有效地识别刀具破损。