CN101614533B - 一种精确测量超薄工件厚度的方法及仪器 - Google Patents

一种精确测量超薄工件厚度的方法及仪器 Download PDF

Info

Publication number
CN101614533B
CN101614533B CN2008100120296A CN200810012029A CN101614533B CN 101614533 B CN101614533 B CN 101614533B CN 2008100120296 A CN2008100120296 A CN 2008100120296A CN 200810012029 A CN200810012029 A CN 200810012029A CN 101614533 B CN101614533 B CN 101614533B
Authority
CN
China
Prior art keywords
thickness
omega
testee
frequency
ultrasonic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN2008100120296A
Other languages
English (en)
Other versions
CN101614533A (zh
Inventor
蔡桂喜
韩晓华
刘畅
徐华
董瑞琪
贾中青
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Metal Research of CAS
Original Assignee
Institute of Metal Research of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Metal Research of CAS filed Critical Institute of Metal Research of CAS
Priority to CN2008100120296A priority Critical patent/CN101614533B/zh
Publication of CN101614533A publication Critical patent/CN101614533A/zh
Application granted granted Critical
Publication of CN101614533B publication Critical patent/CN101614533B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

一种能精确测量超薄工件厚度方法,采用信号处理的方法来提取表征工件厚度的超声波信息,既适合于测量超薄大曲率工件厚度也适合于测量普通工件厚度的超声测厚技术,它的核心是对接收到的回波信号进行快速傅里叶变换从而获得测量结果;一种能精确测量超薄工件厚度方法的仪器,由探头和主机两部分组成,主机部分包括DSP控制模块、超声波发射电路、超声波信号接收电路,信号放大电路、闸门电路、通讯接口,液晶显示器和键盘。本发明的优点:既适合于测量超薄大曲率工件厚度也适合于测量普通工件厚度,可实时采集超声回波信号,实时的进行快速傅里叶变换,实时计算和存储管壁厚度,可以实现机械化自动测量。

Description

一种精确测量超薄工件厚度的方法及仪器
技术领域
本发明涉及无损检测领域,特别提供了一种能精确测量超薄工件厚度方法及仪器。 
背景技术
超声波测厚技术是无损检测技术领域的重要部分,广泛应用于石油、化工、冶金、造船、航空、航天等各个领域。国内外现有的超声测厚仪都是根据超声波脉冲反射原理来进行厚度测量,超声探头发射的超声脉冲波经耦合剂进入被测物体,在被测物体内传播至底面时发生反射,反射回来的超声波(回波)被超声探头接收,超声脉冲在被测物体两表面之间多次往返形成多次回波,测得相邻两次回波之间的时间间隔t,根据超声波在材料中的传播声速c,即可计算出被测物体的厚度d=ct/2。国内外按回波原理设计的超声测厚仪一般只能测量0.75mm以上的材料厚度,在实际测量过程中,厚度越薄,其测量的准确性越低。而且显示分辨率一般为0.1mm,常常不能满足实际应用中测量精确度的要求。当被测物体较厚时,相邻两次回波之间的距离宽;当被测物体较薄时,相邻两次回波之间的距离窄。当被测物体的厚度薄到一定程度时,相邻两次回波发生重叠,二者的时间间隔无法测量,也就无法得到被测物体的厚度。大曲率工件由于曲率大,超声反射回波信号弱,不能形成回波系列,工件的厚度无法测量。因此对于超薄工件尤其是大曲率的超薄工件,现有的超声回波仪器无能为力。而这 样的工件,例如小口径薄壁管,在现代飞机和核技术中不仅是不可缺少的,而且需要一定的批量。 
发明内容
本发明的目的是为了测量超薄工件尤其是大曲率的超薄工件的厚度,提供一种能精确测量超薄工件厚度方法及仪器。 
本发明提供了一种能精确测量超薄工件厚度方法,其特征在于:能精确测量超薄工件厚度方法为采用信号处理的方法来提取表征工件厚度的超声波信息,既适合于测量超薄大曲率工件厚度也适合于测量普通工件厚度的超声测厚技术,它的核心是对接收到的回波信号进行快速傅里叶变换(FFT)从而获得测量结果;超声脉冲有一定的频谱范围,见图3,当被测物体较薄时,相邻的回波重叠在一起,必然产生干涉现象,某些频率的超声波振动相互加强,而另一些频率的超声波振动相互减弱或完全抵消;通过快速傅里叶变换将接收到的时域信号转换为频域信号,也叫做频谱曲线,此时频域信号就会显示出周期性的峰值和谷值;在频率为基频f0=v/4x的奇数倍时出现周期性极小值,偶数倍时出现周期性极大值,其中x为被测物体的厚度,v为超声波在被测物体中的传播速度,因此测出频域信号达到极大值的频率fn可以计算出被测物体的厚度x, x = nv 2 f n n=1,2,3,Λ式中,fn 是频谱曲线上第n个极大值所对应的频率。 
所述的能精确测量超薄工件厚度方法,超声波经被测物体底面反射的 一次回波用f1(t)表示,二次回波用f2(t)表示, f 2 ( t ) = α f 1 ( t - 2 x v ) , α<1; 
F 1 ( ω ) = ∫ - ∞ + ∞ f 1 ( t ) e - iωt dt ,
F 2 ( ω ) = ∫ - ∞ + ∞ f 2 ( t ) e - iωt dt = ∫ - ∞ + ∞ α f 1 ( t - 2 x v ) e - iωt dt = α e - iω 2 x c ∫ - ∞ + ∞ f 1 ( t ) e - iωt dt = α e - iω 2 x v F 1 ( ω )
其中,x为被测物体的厚度,v为超声波在被测物体中的传播速度;根据傅里叶变换的线性性质,函数一次回波与二次回波的叠加信号f1(t)+f2(t)的傅里叶变换为F1(ω)+F2(ω),即 
F ( ω ) = F 1 ( ω ) + F 2 ( ω ) = F 1 ( ω ) + α e - iω 2 x v F 1 ( ω ) = F 1 ( ω ) · ( 1 + α e - iω 2 x v ) ,
ω 2 x v = β = 2 πf v 2 x
F(ω)=F1(ω)·(1+αe-iβ)=F1(ω)·(1+αcosβ-iαsinβ) 
当β=2nπ,即f=2nv/4x(n=0,1,2Λ)时,F(ω)=(1+α)F1(ω)有最大值; 
当β=(2n+1)π,即f=(2n+1)v/4x (n=0,1,2Λ)时,F(ω)=1-α有最小值; 
即干涉信号的幅度谱在频率为基频f0=v/4x的偶数倍时出现周期性极大值点,奇数倍时出现周期性极小值点;当超声波在被测材料中的声速v已知,根据超声波的干涉规律对回波信号进行幅度谱分析,若f0为幅度谱的极大值时,计算出被测物体的厚度为x=v/4f0的偶数倍;若f0为幅度谱的极小值时,计算出被测物体的厚度为x=v/4f0的奇数倍,由此计算出被测物体的厚度x。 
一种按权利要求1所述的能精确测量超薄工件厚度方法的仪器,其特征在于:所述的能精确测量超薄工件厚度方法的仪器由探头和主机两部分组成,主机部分包括DSP控制模块、超声波发射电路、超声波信号接收电路,信号放大电路、闸门电路、通讯接口,液晶显示器和键盘。 
探头的选择上若采用接触工件的方式测厚始波和界面波重合,同时由 于测量薄工件的超声反射回波也淹没在始波中,因此需要将始波与界面波分离开,可选用水浸探头实现信号延迟,或选用直探头前加有机玻璃块将界面波延迟;由于探头的频率特性,不同频率的探头有不同的工作频带,通过对上述测厚方法的分析可知,被测物体的厚度越薄,其峰值频率fn越高,因此要选用频率高的探头,根据被测物体的厚度范围来选择合适的探头频率。 
本发明的软件功能包括超声仪器的控制、信号的采集和处理以及厚度的计算。本发明的测厚流程图如图9所示,厚度信息提取子程序的流程图如图2所示。本发明的测厚方法首先要对与被测材料同种材质的标准试块的厚度进行测量,并得出线性拟合曲线,然后通过信号处理的方法测量出被测材料的厚度值,根据标定的拟合曲线得出被测材料的厚度校正值。 
本发明的优点: 
本发明的特点是既适合于测量超薄大曲率工件厚度也适合于测量普通工件厚度,性能稳定,能满足实际应用中测量精确度的要求。该仪器在检测时,可实时采集超声回波信号,实时的进行快速傅里叶变换,实时计算和存储管壁厚度,可以实现机械化自动测量。 
附图说明
下面结合附图及实施方式对本发明作进一步详细的说明: 
图1为本发明所述的超声波测厚仪的示意框图; 
图2为DSP信号处理软件系统的示意框图; 
图3为5M探头频谱图; 
图4为塞尺1.0mm采集波形; 
图5为实施例1的FFT波形; 
图6为10M探头的频谱图; 
图7为Ф6*0.7薄壁管采集波形; 
图8为实施例3的FFT波形; 
图9为软件流程图。 
具体实施方式
实施例1: 
本发明提供了一种能精确测量超薄工件厚度方法,其特征在于:能精确测量超薄工件厚度方法为采用信号处理的方法来提取表征工件厚度的超声波信息,既适合于测量超薄大曲率工件厚度也适合于测量普通工件厚度的超声测厚技术,它的核心是对接收到的回波信号进行快速傅里叶变换(FFT)从而获得测量结果;超声脉冲有一定的频谱范围,见图3,当被测物体较薄时,相邻的回波重叠在一起,必然产生干涉现象,某些频率的超声波振动相互加强,而另一些频率的超声波振动相互减弱或完全抵消;通过快速傅里叶变换将接收到的时域信号转换为频域信号,也叫做频谱曲线,此时频域信号就会显示出周期性的峰值和谷值;在频率为基频f0=v/4x的奇数倍时出现周期性极小值,偶数倍时出现周期性极大值,其中x为被测物体的厚度,v为超声波在被测物体中的传播速度,因此测出频域信号达到极大值的频率fn可以计算出被测物体的厚度x, x = nv 2 f n n=1,2,3,Λ式中,fn 是频谱曲线上第n个极大值所对应的频率。 
所述的能精确测量超薄工件厚度方法,超声波经被测物体底面反射的 一次回波用f1(t)表示,二次回波用f2(t)表示, f 2 ( t ) = α f 1 ( t - 2 x v ) , α<1; 
F 1 ( ω ) = ∫ - ∞ + ∞ f 1 ( t ) e - iωt dt ,
F 2 ( ω ) = ∫ - ∞ + ∞ f 2 ( t ) e - iωt dt = ∫ - ∞ + ∞ α f 1 ( t - 2 x v ) e - iωt dt = α e - iω 2 x c ∫ - ∞ + ∞ f 1 ( t ) e - iωt dt = α e - iω 2 x v F 1 ( ω )
其中,x为被测物体的厚度,v为超声波在被测物体中的传播速度;根据傅里叶变换的线性性质,函数一次回波与二次回波的叠加信号f1(t)+f2(t)的傅里叶变换为F1(ω)+f2(ω),即 
F ( ω ) = F 1 ( ω ) + F 2 ( ω ) = F 1 ( ω ) + α e - iω 2 x v F 1 ( ω ) = F 1 ( ω ) · ( 1 + α e - iω 2 x v ) ,
ω 2 x v = β = 2 πf v 2 x
F(ω)=F1(ω)·(1+αe-iβ)=F1(ω)·(1+αcosβ-iαsinβ) 
当β=2nπ,即f=2nv/4x(n=0,1,2Λ)时,F(ω)=(1+α)F1(ω)有最大值; 
当β=(2n+1)π,即f=(2n+1)v/4x(n=0,1,2Λ)时,F(ω)=1-α有最小值; 
即干涉信号的幅度谱在频率为基频f0=v/4x的偶数倍时出现周期性极大值点,奇数倍时出现周期性极小值点;当超声波在被测材料中的声速v已知,根据超声波的干涉规律对回波信号进行幅度谱分析,  若f0为幅度谱的极大值时,计算出被测物体的厚度为x=v/4f0的偶数倍;若f0为幅度谱的极小值时,计算出被测物体的厚度为x=v/4f0的奇数倍,由此计算出被测物体的厚度x。 
所述的能精确测量超薄工件厚度方法首先要对与被测材料同种材质的标准试块的厚度进行测量,并得出线性拟合曲线,然后通过信号处理的方法测量出被测材料的厚度值,根据标定的拟合曲线得出被测材料的厚度校正值。 
一种按权利要求1所述的能精确测量超薄工件厚度方法的仪器,其特 征在于:所述的能精确测量超薄工件厚度方法的仪器由探头和主机两部分组成,主机部分包括DSP控制模块、超声波发射电路、超声波信号接收电路,信号放大电路、闸门电路、通讯接口,液晶显示器和键盘。 
探头的选择上采用MATEC公司5MHz频率水浸聚焦探头,测量塞尺的厚度。 
5MHz频率水浸聚焦探头的频谱图如图3所示。从图中可以看出其工作频率有一定的带宽范围。 
测量1.0mm塞尺的厚度,通过调节闸门选择超声回波信号需要采集的数据,其波形如图4所示,对采集数据进行快速傅里叶变换如图5所示。从图5可以看出在频率f1=2.83203MHz时出现峰值,取超声波在材料中的声速v=5900m/s,根据公式计算出厚度值x=1.0417mm。 
根据上述方法分别对厚度为0.75mm、0.5mm、0.4mm、0.3mm塞尺进行测量,测量值如表1所示。 
表一5M探头塞尺厚度测量值对比表 
实际值(mm) 峰值频率(MHz) 测量值(mm) 拟合值(mm) 偏差(mm)
 1.0  2.83203  1.0417  0.9943 0.0057
 0.75  3.71093  0.7949  0.7533 -0.0033
 0.5  5.37109  0.5492  0.5133 -0.0133
 0.4  6.83593  0.4315  0.3983  0.0017
 0.3  9.17968  0.3214  0.2908  0.0092
对厚度测量值进行线性拟合处理后得到被测物体厚度的拟合值,与实际值之间的偏差在几十μm数量级。 
实施例2: 
本发明提供了一种能精确测量超薄工件厚度方法,其特征在于:能精 确测量超薄工件厚度方法为采用信号处理的方法来提取表征工件厚度的超声波信息,既适合于测量超薄大曲率工件厚度也适合于测量普通工件厚度的超声测厚技术,它的核心是对接收到的回波信号进行快速傅里叶变换(FFT)从而获得测量结果;超声脉冲有一定的频谱范围,见图1,当被测物体较薄时,相邻的回波重叠在一起,必然产生干涉现象,某些频率的超声波振动相互加强,而另一些频率的超声波振动相互减弱或完全抵消;通过快速傅里叶变换将接收到的时域信号转换为频域信号,也叫做频谱曲线,此时频域信号就会显示出周期性的峰值和谷值;在频率为基频f0=v/4x的奇数倍时出现周期性极小值,偶数倍时出现周期性极大值,其中x为被测物体的厚度,v为超声波在被测物体中的传播速度,因此测出频域信号达到极大值的频率fn可以计算出被测物体的厚度x, x = nv 2 f n n=1,2,3,Λ式中,fn 是频谱曲线上第n个极大值所对应的频率。 
所述的能精确测量超薄工件厚度方法,超声波经被测物体底面反射的一次回波用f1(t)表示,二次回波用f2(t)表示, f 2 ( t ) = α f 1 ( t - 2 x v ) , α<1; 
F 1 ( ω ) = ∫ - ∞ + ∞ f 1 ( t ) e - iωt dt ,
F 2 ( ω ) = ∫ - ∞ + ∞ f 2 ( t ) e - iωt dt = ∫ - ∞ + ∞ α f 1 ( t - 2 x v ) e - iωt dt = α e - iω 2 x c ∫ - ∞ + ∞ f 1 ( t ) e - iωt dt = α e - iω 2 x v F 1 ( ω )
其中,x为被测物体的厚度,v为超声波在被测物体中的传播速度;根据傅里叶变换的线性性质,函数一次回波与二次回波的叠加信号f1(t)+f2(t)的傅里叶变换为F1(ω)+F2(ω),即 
F ( ω ) = F 1 ( ω ) + F 2 ( ω ) = F 1 ( ω ) + α e - iω 2 x v F 1 ( ω ) = F 1 ( ω ) · ( 1 + α e - iω 2 x v ) ,
ω 2 x v = β = 2 πf v 2 x
F(ω)=F1(ω)·(1+αe-iβ)=F1(ω)·(1+αcosβ-iαsinβ) 
当β=2nπ,即f=2nv/4x (n=0,1,2Λ)时,F(ω)=(1+α)F1(ω)有最大值; 
当β=(2n+1)π,即f=(2n+1)v/4x (n=0,1,2Λ)时,F(ω)=1-α有最小值; 
即干涉信号的幅度谱在频率为基频f0=v/4x的偶数倍时出现周期性极大值点,奇数倍时出现周期性极小值点;当超声波在被测材料中的声速v已知,根据超声波的干涉规律对回波信号进行幅度谱分析,  若f0为幅度谱的极大值时,计算出被测物体的厚度为x=v/4f0的偶数倍;若f0为幅度谱的极小值时,计算出被测物体的厚度为x=v/4f0的奇数倍,由此计算出被测物体的厚度x。 
一种按权利要求1所述的能精确测量超薄工件厚度方法的仪器,其特征在于:所述的能精确测量超薄工件厚度方法的仪器由探头和主机两部分组成,主机部分包括DSP控制模块、超声波发射电路、超声波信号接收电路,信号放大电路、闸门电路、通讯接口,液晶显示器和键盘。 
探头的选择上若采用接触工件的方式测厚始波和界面波重合,同时由于测量薄工件的超声反射回波也淹没在始波中,因此需要将始波与界面波分离开,可选用水浸探头实现信号延迟,或选用直探头前加有机玻璃块将界面波延迟;由于探头的频率特性,不同频率的探头有不同的工作频带,通过对上述测厚方法的分析可知,被测物体的厚度越薄,其峰值频率fn越高,因此要选用频率高的探头,根据被测物体的厚度范围来选择合适的探头频率。 
所述的能精确测量超薄工件厚度方法首先要对与被测材料同种材质的标准试块的厚度进行测量,并得出线性拟合曲线,然后通过信号处理的方 法测量出被测材料的厚度值,根据标定的拟合曲线得出被测材料的厚度校正值。 
根据实施例1和测厚原理可以得出,厚度越薄,其波形的FFT变换的峰值频率越高。若想要测出小于0.3mm的厚度,需要采用10MHz以上的探头。 
采用MATEC公司10MHz频率水浸聚焦探头,测量塞尺的厚度。探头的频谱图如图6所示。 
采用与实施例1同样的测厚方法得到0.75mm、0.5mm、0.4mm、0.3mm、0.25mm、0.2mm塞尺的厚度。 
表二10M探头塞尺测厚 
实际值(mm) 峰值频率(MHz) 测量值(mm) 拟合值(mm) 偏差(mm)
 0.75  3.71093  0.7949  0.7451  0.0049
 0.5  5.46875  0.5394  0.5059 -0.0059
 0.4  6.83593  0.4315  0.4049 -0.0049
 0.3  9.17969  0.3214  0.3019 -0.0019
 0.25  11.32813  0.2604  0.2448 0.0052
 0.2  14.0625  0.2098  0.1974 0.0026
对厚度测量值进行线性拟合处理后得到被测物体厚度的拟合值,与实际值之间的偏差在几个μm数量级。 
实施例3: 
本发明提供了一种能精确测量超薄工件厚度方法,其特征在于:能精确测量超薄工件厚度方法为采用信号处理的方法来提取表征工件厚度的超声波信息,既适合于测量超薄大曲率工件厚度也适合于测量普通工件厚度的超声测厚技术,它的核心是对接收到的回波信号进行快速傅里叶变换(FFT)从而获得测量结果;超声脉冲有一定的频谱范围,见图1,当被测 物体较薄时,相邻的回波重叠在一起,必然产生干涉现象,某些频率的超声波振动相互加强,而另一些频率的超声波振动相互减弱或完全抵消;通过快速傅里叶变换将接收到的时域信号转换为频域信号,也叫做频谱曲线,此时频域信号就会显示出周期性的峰值和谷值;在频率为基频f0=v/4x的奇数倍时出现周期性极小值,偶数倍时出现周期性极大值,其中x为被测物体的厚度,v为超声波在被测物体中的传播速度,因此测出频域信号达到极大值的频率fn可以计算出被测物体的厚度x, n = nv 2 f n n=1,2,3,Λ式中,fn 是频谱曲线上第n个极大值所对应的频率。 
所述的能精确测量超薄工件厚度方法,超声波经被测物体底面反射的一次回波用f1(t)表示,二次回波用f2(t)表示, f 2 ( t ) = α f 1 ( t - 2 x v ) , α<1; 
F 1 ( ω ) = ∫ - ∞ + ∞ f 1 ( t ) e - iωt dt ,
F 2 ( ω ) = ∫ - ∞ + ∞ f 2 ( t ) e - iωt dt = ∫ - ∞ + ∞ α f 1 ( t - 2 x v ) e - iωt dt = α e - iω 2 x c ∫ - ∞ + ∞ f 1 ( t ) e - iωt dt = α e - iω 2 x v F 1 ( ω )
其中,x为被测物体的厚度,v为超声波在被测物体中的传播速度;根据傅里叶变换的线性性质,函数一次回波与二次回波的叠加信号f1(t)+f2(t)的傅里叶变换为F1(ω)+F2(ω),即 
F ( ω ) = F 1 ( ω ) + F 2 ( ω ) = F 1 ( ω ) + α e - iω 2 x v F 1 ( ω ) = F 1 ( ω ) · ( 1 + α e - iω 2 x v ) ,
ω 2 x v = β = 2 πf v 2 x
F(ω)=F1(ω)·(1+αe-iβ)=F1(ω)·(1+αcosβ-iαsinβ) 
当β=2nπ,即f=2nv/4x(n=0,1,2Λ)时,F(ω)=(1+α)F1(ω)有最大值; 
当β=(2n+1)π,即f=(2n+1)v/4x(n=0,1,2Λ)时,F(ω)=1-α有最小值; 
即干涉信号的幅度谱在频率为基频f0=v/4x的偶数倍时出现周期性极 大值点,奇数倍时出现周期性极小值点;当超声波在被测材料中的声速v已知,根据超声波的干涉规律对回波信号进行幅度谱分析,  若f0为幅度谱的极大值时,计算出被测物体的厚度为x=v/4f0的偶数倍;若f0为幅度谱的极小值时,计算出被测物体的厚度为x=v/4f0的奇数倍,由此计算出被测物体的厚度x。 
所述的能精确测量超薄工件厚度方法首先要对与被测材料同种材质的标准试块的厚度进行测量,并得出线性拟合曲线,然后通过信号处理的方法测量出被测材料的厚度值,根据标定的拟合曲线得出被测材料的厚度校正值。 
一种按权利要求1所述的能精确测量超薄工件厚度方法的仪器,其特征在于:所述的能精确测量超薄工件厚度方法的仪器由探头和主机两部分组成,主机部分包括DSP控制模块、超声波发射电路、超声波信号接收电路,信号放大电路、闸门电路、通讯接口,液晶显示器和键盘。 
探头的选择上若采用接触工件的方式测厚始波和界面波重合,同时由于测量薄工件的超声反射回波也淹没在始波中,因此需要将始波与界面波分离开,可选用水浸探头实现信号延迟,或选用直探头前加有机玻璃块将界面波延迟;由于探头的频率特性,不同频率的探头有不同的工作频带,通过对上述测厚方法的分析可知,被测物体的厚度越薄,其峰值频率fn越高,因此要选用频率高的探头,根据被测物体的厚度范围来选择合适的探头频率。 
国产超声回波测厚仪对于管子的测量有一定的测量下限值Ф20mm*2.0mm,绝大多数航空或核能用小口径薄壁管的管径和壁厚远小于上 述下限值,国外名牌产品测厚的下限值为0.5mm,管径的下限值未见报道。用本专利发明的测厚方法能够精确测量目前航空或核能用的所有的小口径薄壁管的壁厚。 
测量Ф6mm*0.7mm的薄壁管,采集的数据波形如图7所示,通过对采集数据进行快速傅里叶变换如图8所示。从图8可以看出在频率f=3.90625MHz时出现峰值,取超声波在材料中的声速v=5900m/s,根据公式计算出厚度值x=0.7552mm。按前述的实施例得出的拟合方程:Y=A+B*X,A=-0.001073984,B=1.068268293,可计算出实际测量的壁厚值X=(0.7552-0.0011)/1.068268293=0.7059mm,可见测量的精确度很好。 
该仪器在检测时,可实时采集超声回波信号,实时的进行快速傅里叶变换,实时计算和存储管壁厚度,可以实现机械化自动测量。 

Claims (4)

1.一种精确测量超薄工件厚度的方法,其特征在于:该方法为采用信号处理的方法来提取表征工件厚度的超声波信息,是一种既适合于测量超薄大曲率工件厚度也适合于测量普通工件厚度的超声测厚技术,它的核心是对接收到的回波信号进行快速傅里叶变换从而获得测量结果;超声脉冲有一定的频谱范围,当被测物体较薄时,相邻的回波重叠在一起,必然产生干涉现象,某些频率的超声波振动相互加强,而另一些频率的超声波振动相互减弱或完全抵消;通过快速傅里叶变换将接收到的时域信号转换为频域信号,也叫做频谱曲线,此时频域信号就会显示出周期性的峰值和谷值;在频率为基频f0=v/4x的奇数倍时出现周期性极小值,偶数倍时出现周期性极大值,其中x为被测物体的厚度,v为超声波在被测物体中的传播速度,因此测出频域信号达到极大值的频率fn可以计算出被测物体的厚度x,
Figure FSB00000531830800011
式中,fn是频谱曲线上第n个极大值所对应的频率。
2.按照权利要求1所述的精确测量超薄工件厚度的方法,其特征在于:超声波经被测物体底面反射的一次回波用f1(t)表示,二次回波用f2(t)表示,
f 2 ( t ) = &alpha;f 1 ( t - 2 x v ) , &alpha; < 1 ;
F 1 ( &omega; ) = &Integral; - &infin; + &infin; f 1 ( t ) e - i&omega;t dt ,
F 2 ( &omega; ) = &Integral; - &infin; + &infin; f 2 ( t ) e - i&omega;t dt = &Integral; - &infin; + &infin; &alpha; f 1 ( t - 2 x v ) e - i&omega;t dt = &alpha; e - i&omega; 2 x c &Integral; - &infin; + &infin; f 1 ( t ) e - i&omega;t dt = &alpha; e - i&omega; 2 x v F 1 ( &omega; )
其中,x为被测物体的厚度,v为超声波在被测物体中的传播速度;根据傅里叶变换的线性性质,函数一次回波与二次回波的叠加信号f1(t)+f2(t)的傅里叶变换为F1(ω)+F2(ω),即
F ( &omega; ) = F 1 ( &omega; ) + F 2 ( &omega; ) = F 1 ( &omega; ) + &alpha; e - i&omega; 2 x v F 1 ( &omega; ) = F 1 ( &omega; ) &CenterDot; ( 1 + &alpha;e - i&omega; 2 x v ) ,
&omega; 2 x v = &beta; = 2 &pi;f v 2 x
F(ω)=F1(ω)·(1+αe-iβ)=F1(ω)·(1+αcosβ-iαsinβ)
当β=2nπ,即f=2nv/4x(n=0,1,2…)时,F(ω)=(1+α)F1(ω)有最大值;
当β=(2n+1)π,即f=(2n+1)v/4x(n=0,1,2…)时,F(ω)=1-α有最小值;
即干涉信号的幅度谱在频率为基频f0=v/4x的偶数倍时出现周期性极大值点,奇数倍时出现周期性极小值点;当超声波在被测材料中的声速v已知,根据超声波的干涉规律对回波信号进行幅度谱分析,若f0为幅度谱的极大值时,计算出被测物体的厚度为x=v/4f0的偶数倍;若f0为幅度谱的极小值时,计算出被测物体的厚度为x=v/4f0的奇数倍,由此计算出被测物体的厚度x。
3.按照权利要求1所述的精确测量超薄工件厚度的方法,其特征在于:首先要对与被测材料同种材质的标准试块的厚度进行测量,并得出线性拟合曲线,然后通过信号处理的方法测量出被测材料的厚度值,根据标定的拟合曲线得出被测材料的厚度校正值。
4.一种使用权利要求1的所述方法精确测量超薄工件厚度的仪器,其特征在于:所述仪器由探头和主机两部分组成,主机部分包括DSP控制模块、超声波发射电路、超声波信号接收电路、信号放大电路、闸门电路、通讯接口、液晶显示器和键盘。
CN2008100120296A 2008-06-26 2008-06-26 一种精确测量超薄工件厚度的方法及仪器 Expired - Fee Related CN101614533B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008100120296A CN101614533B (zh) 2008-06-26 2008-06-26 一种精确测量超薄工件厚度的方法及仪器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008100120296A CN101614533B (zh) 2008-06-26 2008-06-26 一种精确测量超薄工件厚度的方法及仪器

Publications (2)

Publication Number Publication Date
CN101614533A CN101614533A (zh) 2009-12-30
CN101614533B true CN101614533B (zh) 2012-01-11

Family

ID=41494313

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008100120296A Expired - Fee Related CN101614533B (zh) 2008-06-26 2008-06-26 一种精确测量超薄工件厚度的方法及仪器

Country Status (1)

Country Link
CN (1) CN101614533B (zh)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009060106A1 (de) * 2009-12-17 2011-06-22 Salzgitter Mannesmann Line Pipe GmbH, 57074 Verfahren zur Prüfung von Verbindungen metallischer Werkstücke mit Kunststoffmassen auf Hohlräume mittels Ultraschall
CN102721747B (zh) * 2012-06-08 2014-05-21 北京航空航天大学 一种非共线非线性超声无损检测方法
CN102865839B (zh) * 2012-09-21 2015-07-01 华南理工大学 一种基于宽带调频及接收补偿的超声波测厚方法及装置
CN103308012B (zh) * 2013-06-24 2016-04-27 广东惠利普路桥信息工程有限公司 混凝土路面厚度检测系统及检测方法
CN103486987B (zh) * 2013-10-14 2016-03-30 上海电力学院 一种利用曲线拟合提高超声测厚精度的方法
JP6388545B2 (ja) * 2015-01-16 2018-09-12 株式会社ディスコ 被加工物の研削方法
CN105180853B (zh) * 2015-09-14 2017-05-03 哈尔滨工业大学 一种电磁超声金属材料测厚方法
CN108279417A (zh) * 2018-01-08 2018-07-13 上海应用技术大学 超声波测量薄壁复合钢管厚度的方法
CN108759756A (zh) * 2018-04-04 2018-11-06 成都众宜天成科技有限公司 基于镍氢电池持续测量的涂料涂层厚度精准测量仪
CN108426519A (zh) * 2018-04-04 2018-08-21 成都众宜天成科技有限公司 一种基于温度补偿电路的涂料涂层厚度精准测量仪
CN108759757A (zh) * 2018-04-04 2018-11-06 成都众宜天成科技有限公司 一种涂料涂层厚度测量仪
CN109099864A (zh) * 2018-09-17 2018-12-28 中南大学 一种硅烷浸渍深度检测装置及方法
CN109781041B (zh) * 2019-02-28 2021-04-20 华中科技大学 一种基于频域分析的电磁超声测厚方法
CN110500974B (zh) * 2019-08-06 2021-01-05 天津大学 基于改进峰值识别的工件厚度检测方法
CN110702042B (zh) * 2019-10-15 2021-07-02 河海大学常州校区 一种超声脉冲回波测厚中信号混叠条件下的厚度估算方法
CN111486804B (zh) * 2020-06-15 2021-10-08 东莞职业技术学院 用于精密部件厚度测量的信号处理方法及测量方法
CN112815886A (zh) * 2020-12-31 2021-05-18 南京晨光集团有限责任公司 一种薄壁回转体工件壁厚测量系统
CN113188489B (zh) * 2021-04-29 2022-10-21 深圳市麒博精工科技有限公司 检测薄平板材料厚度一致性的超声反射谱方法
CN114111659B (zh) * 2021-11-12 2024-10-01 西安热工研究院有限公司 一种可提高超声测厚精度的超声检测方法及装置
CN114353656B (zh) * 2022-01-11 2024-01-30 中北大学 一种渐变曲率工件涂层厚度的涡流测量曲面校正方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5974886A (en) * 1997-09-15 1999-11-02 General Electric Company Method and apparatus for thickness determination in multilayer articles
CN1496475A (zh) * 2001-03-28 2004-05-12 ����ͨAnp���޹�˾ 通过迭加多个回波周期测量弱反射的子层厚度的超声方法
CN1825056A (zh) * 2005-02-22 2006-08-30 哈茨有限公司 对冶金炉等进行不间断和无损检查的系统、方法和装置
CN101109628A (zh) * 2006-07-20 2008-01-23 国际商业机器公司 用于测量固态器件中的薄层的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5974886A (en) * 1997-09-15 1999-11-02 General Electric Company Method and apparatus for thickness determination in multilayer articles
CN1496475A (zh) * 2001-03-28 2004-05-12 ����ͨAnp���޹�˾ 通过迭加多个回波周期测量弱反射的子层厚度的超声方法
CN1825056A (zh) * 2005-02-22 2006-08-30 哈茨有限公司 对冶金炉等进行不间断和无损检查的系统、方法和装置
CN101109628A (zh) * 2006-07-20 2008-01-23 国际商业机器公司 用于测量固态器件中的薄层的方法

Also Published As

Publication number Publication date
CN101614533A (zh) 2009-12-30

Similar Documents

Publication Publication Date Title
CN101614533B (zh) 一种精确测量超薄工件厚度的方法及仪器
EP2799820B1 (en) Liquid surface level measurement device, method, and program
CN103245311B (zh) 用超声检测多层吸波涂层测厚装置的测厚方法
CN110082428B (zh) 基于裂纹尖端效应结合激光超声测量材料弹性常数的方法
NO343972B1 (no) Metode og system for fluidanalyse
CN103615996A (zh) 一种超声信号频谱滤波技术无损测量涂层厚度的方法
CN103543208B (zh) 基于频谱分析原理减小tofd检测近表面盲区的方法
EP4130655A1 (en) Method and system for determining a thickness of an elongate or extended structure
CN102608212A (zh) 基于声压反射系数功率谱测量薄层声阻抗和声衰减的方法
JP2013130572A (ja) 超音波厚さ測定方法及び装置
CN104457635A (zh) 基于Welch法谱估计的超薄涂层厚度均匀性无损检测方法
CN104698089A (zh) 一种适用于倾斜裂纹定量和成像的超声相对时间传播技术
CN103075981B (zh) 一种超声波测厚方法
CN109696482A (zh) 一种弯头腐蚀的相控阵超声柔性探头检测方法
CN104897777A (zh) 基于Burg算法的自回归谱外推技术提高TOFD检测纵向分辨率的方法
EP3516329B1 (en) Determining a wall thickness of an object
CN107024535A (zh) 一种基于表面波的垂直缺陷的多系数深度检测方法
CN103616436B (zh) 一种接触刚度的高精度超声检测方法
RU108627U1 (ru) Система ультразвуковой дефектоскопии трубопровода
Kumar et al. Experimentation for sag and dimension measurement of thin-walled tubes and pipes using multi-channel ultrasonic imaging system
CN105044215A (zh) 一种非破坏性的材料声速现场测量方法
CN111189513B (zh) 一种超声波液位测量方法、装置及超声波液位计
CN202947698U (zh) 锅炉受热面管管壁金属层和氧化层厚度的测量装置
Silk Changes in ultrasonic defect location and sizing
CN202304777U (zh) 一种工件厚度测量装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120111

Termination date: 20170626