《武汉工程大学学报》 2012年4期
64-68
出版日期：2012-05-30
ISSN:1674-2869
CN:42-1779/TQ
再制造对象磁记忆信号的特征提取方法
0引言 一种新的金属诊断技术——磁记忆检测技术（Metal Magnetic Memory,MMM）由Doubov提出［1］.它可准确可靠地测出被测对象上以应力集中为特征的危险部位，为金属构件寿命评估提供依据.MMM在再制造寿命预测领域也得到了广泛运用［2］，但在对再制造对象进行寿命预测时，常常需要考虑再制造对象的历史健康信息，采集构件在役运行状态的磁记忆信号.因此，采集得到的磁记忆信号往往容易受到现场环境磁场和构件表面状态的影响，附带有大量的噪声，对磁记忆信号真实性造成严重干扰，从而影响信号对缺陷识别的指示作用.由于磁记忆信号具有非平稳性，能量主要集中在低频段，为此，将小波变换用于分析环境磁场干扰下的磁记忆信号，去除高频噪声，提高信噪比，研究磁记忆检测信号的特征提取方法，为准确预测铁磁性材料再制造对象的剩余寿命提供依据.1磁记忆检测技术的机理 MMM能用于判断缺陷或应力集中的位置的根本原因在于铁磁性金属具有的磁记忆效应.受地球磁场激励作用，在铁磁性金属应力集中区域会发生磁畴组织的定向或不可逆的重新取向，形成类似缺陷的漏磁场分布，外载消除后该状态会继续保留，还与最大作用应力有关，金属构件的这一特性即所谓的磁记忆效应［3］.正是由于磁记忆效应的存在，实际应用的金属构件在应力集中区域的漏磁场强度切向分量Hp(x)具有最大值，法向分量Hp(y)改变符号并具有零值,如图1所示［4］.根据此特征进行漏磁信号的采集便可确定作为再制造对象的金属零件上的危险区，继而可以实现再制造对象的强度评估和寿命预测.图1磁记忆检测原理图
Fig.1Principle diagram of MMM 表征磁记忆检测信号特征的参数，除了法向分量Hp(y)外，还有一个重要的参数是磁场法向分量Hp(y)梯度值K.K可用于度量应力集中线附近的应力集中程度［5］.可定义为：K=|ΔHp(y)|2LK（1）
式（1）中, K为应力集中区金属磁场强度变化表征的漏磁场梯度值, (A/m)/mm；ΔHp(y)为位于应力集中线两侧同等线段上两检测点之间磁场Hp(y)的差，A/m；LK为两测点与应力集中线之间的距离，mm.第4期曾寿金，等：再制造对象磁记忆信号的特征提取方法
武汉工程大学学报第34卷
2磁记忆信号处理和小波分析 在再制造对象的磁记忆检测中，通过探头采集的信号包含着各种干扰和噪声，对后续的数据分析会产生很大的影响，因而需要采用适当的信号处理技术，将有用信号有效地分离出来.磁记忆信号所含的噪声来源于许多方面，但主要是受检测现场的强磁环境干扰、检测探头抖动造成的提离干扰以及被检对象的表面缺陷干扰等三方面影响，因而信号频率成分比较复杂，具有极强的不平稳性.由于傅立叶分析方法不具有时-频局部化分析能力，只适合用于分析和处理平稳信号，不适合用来分析磁记忆信号, 因此,引入小波分析方法对磁记忆信号进行处理.2.1小波变换理论 小波变换继承和发展了窗口傅里叶变换的局部化思想，具有良好的时频特性, 适用于处理局部或暂态信号，因而能有效地从信号中提取信息.通过小波基函数的伸缩和平移等运算对观测信号进行多尺度细化分析，既可以观察信号的低频概貌，又可以观察信号的高频细节，解决了傅里叶变换不能解决的许多问题，具有“数学显微镜”作用.因此, 小波变换在信号滤波、去噪、压缩、传递等领域应用广泛. 设Ψ(t)为平方可积函数空间L2(R)中满足下述“容许性条件”的一个函数或者信号:0＜∫0+∞|Ψ(ω)|2ωdω＜+∞（2） 则Ψ(t)称为小波函数或小波母函数，其中Ψ(ω)为Ψ(t)的傅里叶变换;ω为频域自变量. 对于任意的函数或者信号f(t), 其小波变换定义为［6］：Wf(a,b)=＜f(t),Ψa,b(t)＞=1a∫-∞+∞f(t)(t-ba)dt（3）
式（3）中a为尺度参数；b为平移参数；为小波函数的共轭函数.可见，小波变换是一个二元函数，这是小波变换和傅里叶变换最重要的不同点. 在工程实际中,小波变换的尺度参数实质上描述了观测信号的范围，所以小波变换可以理解为信号的多分辨率分析.多分辨率分析的思想是把信号分解成低频逼近和高频细节两部分，每次用不同的分辨率逼近越来越集中的光滑函数. 多分辨率分析只对低频做进一步分解，而不考虑高频部分.设aj表示分解中低频部分，dj表示分解中高频部分，fi∈aj表示分辨率为2-j的函数f的逼近，而dj代表逼近的误差，则有： f=fN+∑Ni=1di （4）
式（4）表明，某一函数或信号f(t)∈L2(R)可以根据分辨率为2-N时f的低频部分和分辨率2-j(1≤j≤N)下的高频部分完全重构.2.2小波阈值去噪方法 一个含有噪声的信号模型可以表示成如下形式：s(ti)=f(ti)+σn(ti)（5）
式（5）中，f(ti)为真实信号；n(ti)为噪声；σ为噪声强度，i=0,1,...,N-1. 小波变换的目的就是通过抑制噪声n(ti)来获取真实信号f(ti).利用小波变换进行去噪的方法主要有3种，即基于小波变换极大值原理去噪，基于小波系数相关性去噪和阈值去噪.由于阈值去噪是一种实现简单，效果较好的小波去噪方法，所以采用阈值法对磁记忆信号进行去噪处理，主要分3步进行： （1）信号的小波分解.选择合适的小波基函数，确定合理的分解层次N，对含噪信号s(ti)进行N层小波分解. （2）高频系数的阈值处理.对经分解的第1至第N层的每层高频系数，选择一个合适的阈值和阈值方法进行量化处理. （3）信号的小波重构.根据第N层的低频系数和经过量化处理后的第1至第N层高频系数，完成信号的重构.2.3磁记忆信号的Hilbert变换 利用磁记忆检测仪和软件直接获得采集信号的梯度, 有时会因为噪声的干扰使得梯度极大值所对应的位置变得难以判断.需要将小波消噪后的信号先进行Hilbert变换再求取梯度, 才能进一步提取应力集中处的梯度极值. 信号的突变信息往往体现在信号的包络里, 最常用的信号包络提取方法是Hilbert变换法.一个时域信号x(t)的Hilbert变换定义为 ［7］：(t)=1πtx(t)=1π∫-∞-∞x(τ)t-τdτ（6） 得到x(t)的解析信号为：W(t)=x(t)+j(t) （7） 对磁记忆信号而言，应力集中区的漏磁信号是奇异的，故其振幅为：A(t)=x(t)2+［(t)］2 （8） A(t)是原始信号的包络, 包含着信号的低频成分, 反映了信号的整体变化趋势.通过小波去噪的信号经过Hilbert变换后再求梯度值，更易于判断信号梯度极大值所对应的位置.3数据的采集与处理3.1数据采集系统及实验方法 采用某一废旧套类零件作为检测试件，材料为Q235钢，试件表面未作任何处理，形状及尺寸如图2所示.磁记忆检测设备选用厦门爱德森公司生产的EMS－2003金属磁记忆检测仪，磁记忆检测原理及试验系统如图3和图4所示.从图3可以看出，磁记忆检测仪主要由探头、信号处理电路以及计算机处理系统三个部分组成.其中，探头由磁敏传感器、温度传感器、距离传感器等传感器组成，信号处理电路由滤波器、放大器和A/D转换器等电子元件组成.图2试件形状及尺寸
Fig.2Shape and size of specimen图3磁记忆检测仪原理框图
Fig.3Principle diagram of MMM instrument 在检测时，EMS-2003金属磁记忆检测仪既可以采用点测的方式，也可以采用扫描检测的方式.笔者在试件表面选取一检测线，采用沿线扫描检测的方式进行检测.为了避免探头与试件表面距离发生改变造成提离干扰而影响检测结果，在测量的过程中要求试件不能挪动.实验分成两种不同的状态进行检测，一种是周边没有环境干扰的原始检测状态，一种是周边存在电机、铁磁性材料及其他强磁干扰环境的状态，两种状态检测的参数完全相同.由于试验的目的是探讨环境磁场干扰下磁记忆检测信号的去噪方法与特征提取方法，为消除探头抖动带来的提离干扰，利用扫描架夹持探头进行定向运动.图4试验系统示意图
Fig.4Schematic diagram of test system3.2数据分析处理 图5所示为由所建试验系统采集到的两种不同状态下的试件磁记忆信号.图5(a)为无干扰状态Hp(y)曲线，图5(b)为存在环境干扰下的Hp(y)曲线.可以看出，在没有干扰状态下，磁记忆信号曲线光滑，能清晰看出信号的变化，但在存在环境干扰的状态下，信号信噪比不高,很难辨认出应力集中部位，失去了对应力集中区的指示作用.因此，需要对环境干扰信号下的检测数据做进一步处理来反映构件的应力集中状况.首先对环境干扰状态下采集的磁记忆信号直接求取梯度值, 结果如图6所示.梯度值是磁记忆信号的另外一个重要特征参数，如本文式（1）所示，在数值上等于两点的磁感应强度模数差与两点距离的比值.由图6可见, 由于噪声干扰,直接求取的梯度值依然含有很强的噪声，还是无法反映真正的低频磁记忆信号的梯度变化. 为进一步提取信号特征，对采集的环境干扰状态磁记忆信号进行小波去噪.主要步骤如下：①选择合适小波基函数.小波基函数在选择的时候一般需要考虑正交性、正则性、对称性、紧支性和消失矩等五个方面因素，但任何一个小波函数都无法同时满足上述所有特性，对于磁记忆信号去噪而言，主要考虑正则性和紧支性.通过去噪效果对比数字实验，最终选择Db4 函数.②确定分解层数.噪声一般具有一定的频带宽度，存在于一定范围的尺度参数的分解结果中.因此需要确定合理的分解层数，使得既能保证完全去除噪声，又能减少处理工作量.经过实验发现，磁记忆信号从4次分解开始信号波形就基本没有变化，故分解层数选为4 层.③选择阈值.阈值的处理方法主要包括硬阈值和软阈值两类，由于硬阈值容易产生不连续性，丢失信号的有用成分，因此采用软阈值方法对采集的磁记忆信号进行处理.通过对小波分解的第1层到第4层的每层高频系数采用Heursure软阈值处理，得到新的小波分解系数, 最后重构信号.图5不同状态磁记忆检测数据
Fig.5Test data of MMM in different states图6环境干扰状态磁记忆信号梯度值
Fig.6Gradient of MMM signal with interference图7是小波滤波去噪后的信号，由图7可见信号已经得到了平滑处理，能够较清晰地看出试件各位置对应的信号变化.对比图5(a)曲线，去噪后的信号和未受干扰影响的信号基本一致，但法向磁场强度未出现过零点，仍然未出现明显的信号特征.进一步对降噪后的信号求取梯度值, 如图8所示.图8中在25附近已经出现了梯度极值点，但不够清晰. 对小波降噪后的图7信号进行Hilbert变换，可以得到信号的包络图.再对信号包络求取梯度值， 其梯度曲线如图9所示.对比图8和图9可以看出,对包络后的信号求取梯度，比直接对信号求取梯度更清晰的显示信号梯度极值特征，能够更准确的表征应力集中区域.图7小波去噪结果
Fig.7Result after wavelet de\|noising 图8小波去噪后的梯度
Fig.8Gradient after wavelet de\|noising图9对信号包络求取梯度
Fig.9Gradient obtained from signal envelope4结语 a.对试件采集的磁记忆信号进行分析可以发现，环境磁场干扰对磁记忆检测结果有着重要的影响.若检测时环境磁场过大，则采集的观测信号会失去对应力集中区的指示作用.因此在对再制造对象进行磁记忆检测时，应尽量避免环境磁场的干扰，或者对含噪信号采用适当的去噪方法，提高应力评价的准确性. b.在磁记忆信号处理中，利用小波阈值去噪技术和Hilbert变换，可以有效地提取磁记忆信号特征，从而进一步判断铁磁性构件的应力集中状况，有望实现金属磁记忆技术对再制造对象剩余寿命的定量检测.