[发明专利]一种基于特征频率的金属构件脉冲涡流测厚方法及设备有效
申请号: | 202010052606.5 | 申请日: | 2020-01-17 |
公开(公告)号: | CN111351422B | 公开(公告)日: | 2021-05-18 |
发明(设计)人: | 武新军;李威 | 申请(专利权)人: | 华中科技大学 |
主分类号: | G01B7/06 | 分类号: | G01B7/06 |
代理公司: | 华中科技大学专利中心 42201 | 代理人: | 尚威;李智 |
地址: | 430074 湖北*** | 国省代码: | 湖北;42 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 基于 特征 频率 金属构件 脉冲 涡流 方法 设备 | ||
1.一种基于特征频率的金属构件脉冲涡流测厚方法,其特征在于,包括如下步骤:
离线拟合阶段:
(1)从已知壁厚的金属构件的脉冲涡流检测信号中提取分析信号u;
(2)对u进行快速傅里叶变换,根据快速傅里叶变换后的频谱设定Laplace小波原子φ(f,ζ,τ,n)的频率f、阻尼比ζ、起始时间τ取值的集合分别为F、Z、T;n∈[1,Ws]且n为正整数,Ws是小波原子φ(f,ζ,τ,n)的序列长度;
(3)对于时刻τ0∈T,计算所有起始时刻为τ0的Laplace小波原子φ(f,ζ,τ0,n)与分析信号u的互相关系数,τ0时刻对应的互相关系数的最大值记为对应的Laplace小波原子的频率记为将作为τ0时刻u的频率;
(4)对于每一时刻τ∈T,重复步骤(3),得到互相关系数的最大值κτ以及频率随时间变化的离散序列κτ(n)、
(5)选择κτ(n)>κ0的时段,在这些时段中统计F={1,2,...,Fmax}中的各频率在序列中出现的次数,κ0是预设的互相关系数值,κ0越大则筛选出的互相关系数越高;
(6)按照步骤(5)的统计结果,顺序计算F={1,2,...,Fmax}中任意三个连续的频率的出现次数之和,选取这些和的最大值对应的三个频率,以这三个频率分别出现的次数为各自的权重,计算这三个频率的加权平均数fw,将fw作为量化试件壁厚的特征量,即特征频率;
(7)选择同一种材料、不同壁厚的金属构件的脉冲涡流检测信号,重复步骤(1)~(6),获得该种材料的金属构件的特征频率-壁厚的拟合曲线;
在线检测阶段:
(8)对待测厚的金属构件,获取其脉冲涡流检测信号并按照步骤(1)~(6)进行处理,得到其特征频率,然后根据其材料对应的特征频率-壁厚的拟合曲线获得其壁厚。
2.如权利要求1所述的基于特征频率的金属构件脉冲涡流测厚方法,其特征在于,所述步骤(2)中,根据快速傅里叶变换后的频谱中的最大频率Fmax设定Laplace小波原子φ(f,ζ,τ,n)的频率f取值的集合为F={1,2,...,Fmax},单位Hz;
设定Laplace小波原子φ(f,ζ,τ,n)的阻尼比ζ取值的集合Z,
设定Laplace小波原子的起始时间τ取值的集合T:
T={1/Fs,2/Fs,...,(N-Ws)/Fs}
其中,Fs是采样频率,N是u的序列长度。
3.如权利要求2所述的基于特征频率的金属构件脉冲涡流测厚方法,其特征在于,所述步骤(3)包括如下子步骤:
(3.1)取Laplace小波的虚部作为小波原子φ(f,ζ,τ,n),小波原子表达式为:
(3.2)小波原子φ(f,ζ,τ,n)与分析信号u的互相关系数κ的计算公式为:
其中,*,*表示内积,||*||2表示向量的二阶范数,向量内积除以向量的二阶范数从而将κ归一化;
(3.3)对于时刻τ0∈T,计算所有起始时间为τ0的小波原子φ(f,ζ,τ0,n)与u的互相关系数选取其中互相关系数的最大值作为τ0时刻的互相关系数的计算结果
其中,且分别是所有起始时间为τ0的小波原子φ(f,ζ,τ0,n)中互相关系数取得最大值的小波原子的频率和阻尼比;
(3.4)将作为τ0时刻u的频率。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于华中科技大学,未经华中科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202010052606.5/1.html,转载请声明来源钻瓜专利网。