一种金属拉深件微小裂纹冲击性信号的识别方法技术

本发明专利技术公开一种金属拉深件微小裂纹冲击性信号的识别方法，先采用声发射传感器对金属拉深件微小裂纹声发射冲击性信号进行采集，对采集的信号进行前置放大、滤波和Ａ／Ｄ转换预处理，将预处理的信号输入计算机中，再进行小波包分析，对小波包分解后的不同频段的数据进行重构，采用时间序列的方法对去噪后的声发射信号进行时序分析，建立时序模型，最后在计算机中结合模糊综合评判方法与最大隶属度原则实现对金属拉深件的状态辨别。本发明专利技术能提高特征参数反映实际工况的准确度，在保证获得足够信息的前提下，可以大大减少采样数目，频率定位准确且清晰，适合于在线监测的高要求场合。

【技术实现步骤摘要】

本专利技术涉及一种基于金属件声发射技术的状态检测及故障诊断方法，特别涉及用 于在拉深状态下对金属制件的早期微小裂纹的识别方法。
技术介绍
拉深是利用专用模具将平板毛坯支撑开口空心零件的一种冲压工艺方法。用拉深 方法可以制成筒形、阶梯形、锥形、球形和其它不规则形状的薄壁零件，如果和其它冲压成 型工艺配合，还可以制造形状极为复杂的零件。用拉深方法来制造薄壁空心件，生产效率 高，省材料，零件的强度和刚度好，精度较高，拉深可加工范围非常广泛，直径从几毫米的小 零件直至2 3m的大型零件，因此，拉深在汽车、航空航天、国防、电器和电子等工业部门以 及日用品生产中占据相当重要的地位。在拉深件的加工检测中，应用最为普遍的是目测。但利用目测只能判断拉深件产 品显而易见的质量问题，如圆角拉裂、凸缘起皱等变形严重的问题。而金属材料在拉深成形 过程中，其不同区域的变形程度、受力情况以及摩擦与润滑等情况千差万别，拉深件在成形 过程中产生裂纹在所难免，由于这些裂纹是在材料内部产生，所以肉眼难以察觉。一旦制件 内部产生裂纹，如果不及时察觉，将导致大量次品的出现。由于早期裂纹比较细小，在产品 使用初期可能并不会影响到产品的性能，但在使用一段时间后，会因为受力、振动、环境等 因素使裂纹进一步扩展，当裂纹扩展到一定程度会造成较大的安全事故和经济损失。金属件拉深过程裂纹声发射的频率较高且在时域上是突发性的非平稳时序信号， 信号频率在100 300KHZ之间，属于高频率声发射信号。此外，金属微观裂纹扩展成为宏 观裂纹之前，需要经过裂纹的慢扩展阶段，理论计算表明，裂纹扩展所需要的能量比裂纹形 成需要的能量约大100 1000倍，因此对金属裂纹识别最好的办法是采用能量分析法。对金属拉深件声发射裂纹检测的传统方法是采用FFT (傅里叶)频谱分析，选取某 些特征频率的幅值来进行诊断，这种方法的缺陷是只考虑了正弦振动的能量，却没有考虑 到频带里信号的所有能量，包括非线性、非平稳的振动能量，而且对于采集的数据量要求较 大，难以满足现场在线监测的高要求场合。
技术实现思路
本专利技术的目的是为了克服现有技术中金属拉深件裂纹特征参数提取的不准确性 及低效率性等不足，提出一种准确性高、效率高、计算速度快的基于小波包分解与自回归谱 分析的金属拉深件早期微小裂纹冲击性声发射信号的识别方法。本专利技术采用如下步骤实现先采用声发射传感器对金属拉深件微小裂纹声发射冲 击性信号进行采集，对采集的信号进行前置放大、滤波和A/D转换预处理，将预处理的信号 输入计算机中，再进行小波包分析，对小波包分解后的不同频段的数据进行重构，然后采用 时间序列的方法对去噪后的声发射信号进行时序分析，建立时序模型，最后在计算机中结 合模糊综合评判方法与最大隶属度原则实现对金属拉深件的状态辨别。本专利技术的有益效果是1、本专利技术针对声发射采集信号的影响，运用时序分析的方法对瞬时冲击性声发射 信号进行特征量的提取，从金属件拉伸过程产生的裂纹声发射信号中提取微弱的缺陷信 息，即早期微小裂纹信号，主要识别过程可在微型计算机上进行，提高了特征参数反映实际 工况的准确度，易操作和实施。2、本专利技术采用小波包技术对金属拉深件早期微小裂纹冲击性声发射的非平稳信 号进行全面的精细的频率分解，可将包括正弦信号在内的任意信号分解到相应的频带内。3、本专利技术采用自回归谱分析建立动态时序模型，对动态数据具有外延特性，对周 期性较强的序列不要求严格按周期采样，在保证获得足够信息的前提下，可以大大减少采 样数目，频率定位准确且清晰，适合于在线监测的高要求场合。附图说明下面结合附图和具体实施方式对本专利技术作进一步详细说明。图1是本专利技术流程示意图。具体实施例方式如图1所示，本专利技术对采集到的金属件裂纹声发射信号进行放大、滤波、A/D转换 及归一化处理前期数据预处理之后，采用小波包技术分解信号至不同频带，根据金属裂纹 及环境噪声的频率特性，剔除噪声及摩擦磨损信号，即对信号进行小波消噪处理；对去噪后 的信号进行重构，而后进行自回归谱分析，选取能量参数，形成辨别缺陷所需的特征向量， 最后采用模糊综合评判实现对裂纹的精确识别。具体步骤如下先采用声发射传感器对金属早期微小裂纹声发射原始纹冲击性信号进行采集，对 采集的信号进行必要的预处理，预处理包括对信号进行前置放大、滤波和A/D转换，预处理 后将信号输入计算机中。再在计算机中对数据信号先进行归一化，而后再进行小波包分析。应用北京鹏翔 科技有限公司PXAES型号声发射系统软件，采用基于矩阵实验室环境下的信号分析软件进 行数据的小波包分解，选择恰当的小波基函数、分解次数。由于状态检测过程数据采集量 大、高频声发射信号具有脉冲衰减性、具有环境噪声等特点，故选用多贝西(Daubechies) 系列小波基，多贝西(Daubechies)系列小波基是典型的紧支撑性正交小波基，这样可以有 效地提取金属拉深件声发射信号的特征。根据金属拉深件声发射信号频率范围特点，将小 波包分解次数定为3次。对上述小波包分解后的不同频段的数据进行重构，重构的依据是根据分解所得不 同频带的信号，将环境噪声信号、摩擦磨损信号对应频带的权数置零，而后合成其它所有频 带的信号，如此完成对声发射信号去噪，提高了待辨别信号的真实性。然后在计算机中采用时间序列的方法对去噪后的声发射信号进行时序分析，建立 时序模型，建立时序模型的内容如下1、采用如下式(1)所示的自回归(AR)模型研究声发射信号的时域、频域及能量特 征，估计自回归模型的参数；Xk= Φ1Χ1 -1+Φ2Χ1 -2+··· + ΦηΧ1 -η+Β1 ,αΑ ^nid (0, CTa )(!)式(1)中χ^Χη，…，乂^表示信号幅值，Φ2，…，小 表示自回归模型参数， ak表示白噪声，σ^表示ak的方差。采用尤耳-沃克(Yule-Walker)方程估计法估计出钟，夠，鄉和σα2这η+1个参数。2、确定模型阶数，根据最小最终预报误差准则(FPE准则)，由样本对模型定阶。根 据模型的预报误差来判明自回归模型的阶数是否合适，最终预报误差定义为FPE(n) = σ · ^^Ν-η(2)式(2)中N是观测数据长度，η为自回归模型的阶数，FPE (η)表不模型阶数的函 数。当η增大时，模型残差方差下降，但(Ν+η)/(Ν-η)值增大，因此，取使FPE(η)最小的η 值作为模型的最佳阶数。FPE准则只适用于判断AR模型的阶次。3、根据下式(3)的系统输入、输出的自功率谱与传递函数的关系，得出时间序列 IxJ的自回归谱图。Sout(w) = Sin(w) |Hx(w) I2(3)这里Sin(W)为输入白噪声的功率谱密度，Sout(W)为输出白噪声的功率谱密度，Ts 为采样时间，Δ为采样间隔，Hx(W)为传递函数，也称为频率响应特性函数，其中Sin(W) = σ a2 · Ts(4)Hx(w) = ^-1 - Σj=l(5)4、根据时间序列{xk}的自回归(AR)谱图，提取声发射信号时序模型的特征参数， 在AR谱中要提取的特征参数有总能量、高频能量、高频能量比、谱峰峰值和中心频率。自回归谱图也是功率谱密度，从中提取能量由下式(6)实现E^fjPSD(I)Af'=I(6)式中N为AR谱图上特定分析频带本文档来自技高网...

【技术保护点】
一种金属拉深件微小裂纹冲击性信号的识别方法，采用声发射传感器对金属拉深件微小裂纹声发射冲击性信号进行采集，对采集的信号进行前置放大、滤波和Ａ／Ｄ转换预处理，将预处理的信号输入计算机中，其特征是采用如下步骤：（Ａ）将计算机中的预处理信号进行小波包分析，对小波包分解后的不同频段的数据进行重构，根据分解所得不同频带的信号，将环境噪声信号、摩擦磨损信号对应频带的权数置零，合成其它所有频带的信号，完成对声发射信号去噪；（Ｂ）在计算机中采用时间序列的方法对去噪后的声发射信号进行时序分析，［１］，ｘ↓［２］，…，ｘ↓［ｎ］｝；２）根据金属拉深件各种资料和实际的经验统计出各种故障状态构成的集合为Ω＝｛ｗ↓［１］，ｗ↓［２］，…，ｗ↓［ｍ］｝；３）首先对Ω中的一个故障ｗ↓［ｉ］（ｉ＝１，２，．．．，ｍ）作单故障模糊评判，确定被评判对象对征兆集ｘ↓［ｊ］（ｊ＝１，２，…，ｎ）的隶属度ｒ↓［ｉｊ］，得出第ｉ个故障ｗ↓［ｉ］的单故障模糊集：ｒ↓［ｉ］＝（ｒ↓［ｉ１］，ｒ↓［ｉ２］，…，ｒ↓［ｉｎ］）由ｍ个故障的评价集构造出一个总的评价模糊矩阵Ｒ＝［＊＊＊］Ｒ即是故障论域Ｕ到征兆域Ｘ的一个模糊关系，μ↓［Ｒ］（ω↓［ｉ］，ｘ↓［ｊ］）＝ｒ↓［ｉｊ］表示ω↓［ｉ］和ｘ↓［ｊ］之间隶属关系的程度；４）用ａ↓［ｉ］（ｉ＝１，２，．．．，ｍ）表示第ｉ个故障ｗ↓［ｉ］在总评价中重要程度的权数，各权数组成的集合为故障论域Ｕ上的模糊子集Ａ，记作：Ａ＝ａ１／ｗ１＋ａ２／ｗ２＋…＋ａｍ／ｗｍ或者Ａ＝（ａ↓［１］，ａ↓［２］，…，ａ↓［ｍ］）式中：ａ↓［ｉ］（０≤ａ↓［ｉ］≤１）为ｗ↓［ｉ］对Ａ的隶属度；５）作模糊变换来进行模糊综合评判：Ｂ＝Ａ.Ｒ＝（ａ↓［１］，ａ↓［２］，…，ａ↓［ｍ］）.［＊＊＊］＝（ｂ↓［１］，ｂ↓［２］，…，ｂ↓［ｍ］）Ｂ为征兆集Ｘ上的模糊子集，ｂ↓［ｊ］（ｊ＝１，２…，ｎ）为征兆ｘ↓［ｊ］综合评判所得模糊子集Ｂ的隶属度；按照最大隶属度原则选择最大的ｂ↓［ｊ］所对应的征兆ｘ↓［ｊ］作为综合评判的结果。建立时序模型，建模内容如下：１）采用式（１）的自回归模型研究声发射信号的时域、频域及能量特征，估计自回归模型的参数；ｘ↓［ｋ］＝φ↓［１］ｘ↓［ｋ－１］＋φ↓［２］ｘ↓［ｋ－２］＋…＋φ↓［ｎ］ｘ↓［ｋ－ｎ］＋ａ↓［ｋ］，ａ↓［ｋ］～ＮＩＤ（０，σ↓［ａ］↑［２］）（１）式中：ｘ↓［ｋ］，ｘ↓［ｋ－１］，…，ｘ↓［ｋ－ｎ］表示信号幅值，...

【技术特征摘要】
一种金属拉深件微小裂纹冲击性信号的识别方法，采用声发射传感器对金属拉深件微小裂纹声发射冲击性信号进行采集，对采集的信号进行前置放大、滤波和A/D转换预处理，将预处理的信号输入计算机中，其特征是采用如下步骤(A)将计算机中的预处理信号进行小波包分析，对小波包分解后的不同频段的数据进行重构，根据分解所得不同频带的信号，将环境噪声信号、摩擦磨损信号对应频带的权数置零，合成其它所有频带的信号，完成对声发射信号去噪；(B)在计算机中采用时间序列的方法对去噪后的声发射信号进行时序分析，建立时序模型，建模内容如下1)采用式(1)的自回归模型研究声发射信号的时域、频域及能量特征，估计自回归模型的参数；xk＝φ1xk 1+φ2xk 2+…+φnxk n+ak，式中xk，xk 1，…，xk n表示信号幅值，φ1，φ2，…，φn表示自回归模型参数，ak表示白噪声，表示ak的方差；采用尤耳 沃克方程估计法估计出和这n+1个参数；2)根据最小最终预报误差准则由样本对模型定阶，根据模型的预报误差来判明自回归模型的阶数是否合适，最终预报误差定义为 <mrow><mi>FPE</mi><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mo>=</mo><msubsup> <mi>&sigma;</mi> <mi>a</mi> <mn>2</mn></msubsup><mo>&CenterDot;</mo><mfrac> <mrow><mi>N</mi><mo>+</mo><mi>n</mi> </mrow> <mrow><mi>N</mi><mo>-</mo><mi>n</mi> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo></mrow> </mrow>式中N是观测数据长度，n为自回归模型的阶数，FPE(n)表示模型阶数的函数，当n增大时，模型残差方差下降，但(N+n)/(N n)值增大，取使FPE(n)最小的n值作为模型的最佳阶数；3)根据式(3)的系统输入、输出的自功率谱与传递函数的关系，得出时间序列{xk}的自回归谱图Sout(w)＝Sin(w)|Hx(w)|2 (3)式中Sin(w)为输入白噪声的功率谱密度，Sout(w)为输出白噪声的功率谱密度，Ts为采样时间，Δ为采样间隔，Hx(w)为传递函数，也称为频率响应特性函数，其中Sin(w)＝σa2·Ts (4)4)根据时间序列{xk}的自回归谱图，提取声发射信号时序模型的特征参数，在自回归谱图中要提取的特征参数有总能量、高频能量、高频能量比、谱峰峰值和中心频率；从自回归谱图提取能量由下式(6)实现 <mrow><mi>E</mi><mo>=</mo><munderover> <mi>&Sigma;</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>N</mi></munderover><mi>PSD</mi><mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo></mrow><mi>&Delta;f</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo></mrow> </mrow>式中N为自回归谱图上特定分析频带内谱线的数目；PSD(i)为对应于第i条谱线的功率谱密度函数值；Δf为频率分辨率；(C)在计算机中结合模糊综合评判方法与最大隶属度原则实现对金属拉深件的状态辨别，步骤如下1)根据金属拉深件各种资料统计出故障发生时所有可能表现出的各种征兆，构成征兆集X＝{x1，x2，…，xn}；2)根据金属拉深件各种资料和实际的经验统计出各种故障状态构成的集合为Ω＝{w1，w2，…，wm}；3)首先对Ω中的一个故障wi(i＝1，2，...，m)作单故障模糊评判，确定被评判对象对征兆集xj(j＝1，2，…，n)的隶属度rij，得出第i个故障wi的单故障模糊集ri＝(ri1，ri2，…，rin)由m个故障的评价集构造出一个总的评价模糊矩阵 <mrow><mi>R</mi><mo>=</mo><mfenced open='[' close=']'> <mtable><mtr> <mtd><msub> <mi>r</mi> <mn>11</mn></msub> </mtd> <mtd><msub> <mi>r</mi> <mn>12</mn></msub> </mtd> <mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo> </mtd> <mtd><msub> <mi>r</mi> <mrow><mn>1</mn><mi>n</mi> </mrow></msub> </mtd></mtr><mtr> <mtd><msub> <mi>r</mi> <mn>21</mn></msub> </mtd> <mtd><msub> <mi>r</mi> <mn>22</mn></msub> </mtd> <mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo> </mtd> <mtd><msub> <mi>r</mi> <mrow><mn>2</mn><mi>n</mi> </mrow></msub> </mtd></mtr><mtr> <mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo> </mtd> <mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo> </mtd> <mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo> </mtd> <mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo> </mtd></mtr><mtr> <mtd><msub> <mi>r</mi> <mrow><mi>m</mi><mn>1</mn> </mrow></msub> </mtd> <mtd><msub> <mi>r</mi> <mrow><mi>m</mi><mn>2</mn> </mrow></msub> </mtd> <mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo> </mtd> <mtd><msub> <mi>r</mi> <mi>mn</mi></msub> </mtd></mtr> </mtable></mfenced> </mrow>R即是故障论域U到征兆域X的一个模糊关系，μR(ωi，xj)＝rij表示ωi和xj之间隶属关系的程度；4)用ai(i＝1，2，.....

【专利技术属性】
技术研发人员：骆志高，胥爱成，何鑫，陈强，
申请(专利权)人：江苏大学，
类型：发明
国别省市：32[中国|江苏]