《基于小波包能量变化率的激光超声处理》由会员分享，可在线阅读，更多相关《基于小波包能量变化率的激光超声处理（8页珍藏版）》请在金锄头文库上搜索。
1、基于小波包能量变化率的激光超声处理 宋潮 郑宾 郭华玲 刘辉 侯静 中北大学计算机与控制工程学院 中北大学电子测试技术重点实验室 摘 要： 为提供解决航空发动机原位裂纹的快速检测方法, 得到缺陷深度与特征量之间的关系, 使用激光器在一系列不同缺陷深度的航空板上激发声表面波, 对收集到的信号进行预处理, 并提出新的小波包变换能量变化率指标, 用于进行缺陷损伤识别。结果表明:反射波信号在 S80、S 81频带上的小波包能量变化率随缺陷深度的增加呈现增加趋势, 增长幅度分别为 96.8%、86.3%;而反射回波在S80、S 83频带上的小波包能量变化率较为显著, 增长幅度分别为23.8%、80.2%
2、。该分析方法为激光超声表征表面缺陷提供新的思路, 为今后从能量变化率指标分析裂纹深度奠定基础。关键词： 激光超声; 小波包分解; 反射波; 能量变化率; 缺陷深度; 作者简介：宋潮 (1991-) , 女, 山西晋城市人, 硕士研究生, 专业方向为激光超声无损检测技术。收稿日期：2017-06-20Laser ultrasonic treatment based on wavelet packet energy change rateSONG Chao ZHENG Bin GUO Hualing LIU Hui HOU Jing Computer Science and Control Eng
3、ineering, North University of China; Abstract： In order to provide a quick detection method for situ crack of aeroengine so as to get the relationship between detect depth and feature quantity, lasers are used on aviation plate with different defect depths to motivate surface acoustic wave, and carr
4、y out pre-treatment for collected signals, and provide a new wavelet packet transformation energy change rate index, which is used for defect damage identification.The results show that the change rate of wavelet packet energy of reflected wave signal in the S80 and S81 bands increases along with th
5、e defect depth with increasing rate respectively being 96.8% and 86.3%; while the change rate of wavelet packet energy of reflected wave in S80 and S83 bands is significant, with increasing rate respectively being 23.8% and 80.2%.The method provides a new idea for laser ultrasonic surface characteri
6、stic defects, which lays a solid foundation for analyzing the crack depth based on energy change rate index in the future.Keyword： laser ultrasonic; wave packet decomposition; reflected wave; energy rate change; defect depth; Received： 2017-06-200 引言航空发动机叶片内部结构复杂, 在成型过程中采用无余量精铸成型1-2。在工艺成型过程中, 经常会导致成
7、型叶片内部存在夹杂、裂缝、表面微裂纹等缺陷, 且叶片的叶身段型面复杂, 在后续的热处理过程中可能会产生各种各样的缺陷, 如气孔、应力集中等;在服役期间, 由于承受高温、高压以及巨大的离心力作用, 极易产生疲劳性裂纹3, 而这些裂纹通常产生在叶片的焊接部位以及根部的近表面位置。微小裂纹随着交变载荷周次的增加以及环境侵蚀时间的延长而逐渐扩展, 随着裂纹尺度的增大, 叶片关键部位的强度将逐渐减小, 最终将会导致空中停车, 带来巨大的安全隐患4-5。因此, 寻找有效的检测手段对于提高叶片质量、保障飞行安全有着非常重要的意义。常规的检测方法由于外场检测时的空间狭小不适合放置换能器, 且不易产生检测表面缺
8、陷的声表面波, 因此, 需要一种新的方法对其进行检测。激光超声检测方法相对于常规超声是一种新型检测方法, 该方法通过脉冲激光入射到固体或液体表面, 介质由于吸收激光能量产生大量热能, 这些热能短时间来不及扩散, 在介质表面附近形成一个很大的温度梯度, 进而引起材料热膨胀, 导致膨胀区域周围介质的约束将产生一个应力分布, 由于应力平衡的需要, 这一应力分布将以瞬态脉冲超声的形式在介质中传播出去, 进而形成超声波。相比于常规超声, 激光超声检测具有非接触、高空间分辨率、探测缺陷范围广以及可对复杂结构进行在线检测的优势6-7。因此, 采用激光超声检测技术对发动机的疲劳裂纹进行检测是可行的。近些年,
9、国内外的专家学者对激光超声应用于微缺陷检测进行了大量卓有成效的研究。Li 等8研究了激光超声产生的声表面波相速度以及频率的变化与软材料机械性能的关系, 总结了杨氏模量与声表面波频率的关系;Kartashev 等9研究了频率衰减和转换器的频率特征对超声探伤中的测量误差的影响;胡海峰等10采用振动声调制技术, 根据测量信号中是否具有振动-超声调制现象对材料损伤进行检测。但是, 现有的研究仍然没有分析裂纹深度与能量变化率之间的关系。针对这一问题, 本文对小波包能量变化率进行适当改进, 采用改进后的能量变换率对缺陷深度信息进行表征, 为这一技术的应用提供参考。1 小波包分析原理超声信号的小波包分析是通
10、过小波包的分解与重构, 提取隐含在微震信号中的特征分量, 并映射到不同频带上, 由于不同频带的能量对于信号特征的变化是十分敏感的, 因此, 可以用频带能量描述不同的超声信号, 从而达到提取信号特征的目的11-13。1.1 小波包分解能量理论小波包元素是由位置、尺度和频率 3 个参数来确定波形的, 根据给定的正交小波函数, 它会组成一组小波包基, 能够将信号的能量保留起来, 并且根据信号自身的特征进行准确的小波包重构14。准 (t) 为正交尺度函数, (t) 为小波函数, 两者之间的关系为其中 h0k、h 1k是在多分辨分析中的滤波器系数。将式 (1) 、式 (2) 进行推广之后, 可得到小波包
11、分解算法公式:小波包重构算法:1.2 小波能量变化率指标对激光超声信号 S (t) 进行小波包分解就是将 S (t) 投影到小波包基上, 得到一系列的小波包系数, 通过分析这些小波包系数, 得到反映缺陷信号的特征。S (t) 的表达式为式中 fi, j (tj) 为激光超声信号小波包分解到节点 (i, j) 上的重构信号, j=0, 1, 2, , 2-1, i=1, 2, , 8。由于小波基函数具有正交性, 因此, 可将小波包变换看作能量守恒, 信号经过小波包变换后能量不变, 由巴什瓦定理可得, 第 i 层信号分量的能量定义为式中:E i, j (tj) 激光超声信号小波包分解到第 i 层第
12、 j 个节点的频带能量;m激光超声信号的采样点数;xj, k重构信号 fi, j (tj) 离散采样点的幅值 (j=0, 1, 2, , 2-1, k=1, 2, , m) 。频带能量 Ei, j (tj) 是由小波函数所确定的频带内的信号能量, 信号的总能为对应于不同频带内的小波包组能量之和。由于频带能量对信号特征的变化十分敏感, 可用于揭示信号的固有特征15, 因此, 本文提出小波包能量变化率指标, 用于对不同缺陷深度的信号进行特征分析。设无损信号的频带能量为 (E i, j) a, 有损信号的频带能量为 (E i, j) b, 定义小波包能量变化率为本文用该指标对缺陷信号进行表征, 其变
13、化特征对缺陷深度变化敏感, 容易进行实际实施。2 实验系统实验装置主要由激励部分和接收部分构成。定义表面微缺陷顶端为 B 点, 缺陷深度为 d, 缺陷宽度为 w, 实验中反射波检测点超声传感器 A 距缺陷顶端 B 的距离为 (501) mm, 线光源 S 距离缺陷顶端 B 的距离为 (181) mm, 其中1 mm 的偏差是由于更换不同缺陷的铝板形成的。选择 Nd-YAG 激光器作为激发装置, 波长为 1 064 nm, 持续时间 7 ns 的激光脉冲经光学调制后为 10 mm30m 的线光源打在铝板上激发声表面波, 超声探头中心频率为 2.5 MHz, 带宽 5 MHz, 使用泰克科技 DP
14、O3034 数字荧光示波器进行数据采集, 每次试验重复 5 次, 采样频率 500 MHz, 采样点数 10 000 个。为避免其他不规则裂纹的干扰, 采用对航空铝件进行人工加工的方式, 以确保研究的对象只有裂纹深度, 选用的航空铝板总体尺寸为 200 mm50 mm8 mm, 6 块铝板的裂纹尺寸分别是无损、0.1 mm0.1 mm、0.1 mm0.2 mm、0.1 mm0.3 mm、0.1 mm0.4 mm、0.1 mm0.5 mm, 图 1 即为实验侧面示意图。图 1 实验过程侧面示意图 下载原图3 激光超声信号的小波包分析本文采用基于小波包分析的小波包能量变化率指标对裂纹深度进行表征,
15、 其处理步骤为:1) 对原始信号进行预处理, 减小噪声的影响;2) 分别对处理后的反射波与反射回波信号进行小波包分解, 得到其频带能量分布;3) 分别计算第 8层信号频带能量的小波包能量变化率指标;4) 分别绘制小波包能量变化率指标柱状图, 进行分析。3.1 信号预处理在超声波无损检测技术中, 信号噪声很大部分是由于结构噪声引起的, 噪声的存在必然会对数据的分析结果产生不良影响, 因此, 在进行数据分析之前, 对数据进行预处理是十分有必要的。本文对数据进行小波去噪处理, 图 2 为激光超声反射波处理前后的对比图, 可以看出处理后的图像更平滑, 波形特征更明显, 更能反映缺陷信息与能量的分布情况
16、, 椭圆中的信号为反射回波信号。图 2 预处理前后对比图 下载原图3.2 反射波信号的小波包分析小波函数的选择是用小波包方法对信号进行分析时必须要考虑的问题, 因为不同的小波基分析相同的信号会产生不同的结果, 会对信号的完整分析产生影响。在对激光超声信号进行小波包分析时, 小波包的选择一般要遵循具有紧支撑性、对称性和光滑性的原则。因此, 本文选择 db8 小波基, 对激光超声反射波信号进行 8 层分解, 通过统计分析得出缺陷深度为 0.10.5 mm 的信号在 S80S87的频带能量分别占总能量的 94.08%、92.07%、93.65%、93.66%、93.37%, 无损信号的频带能量为 96.62%, 即反射波的能量大部分都集中于 S80S87频带上, 因此本文着重对这 8 个频带上的能量进行分析。计算这 8 个频带上的小波包能量变化率, 并绘制柱状图进行对比分析, 见图3。由图可知:1) 反射波信号的小波包
