2023年4月19日发(作者：幼儿园墙绘)基于声压反射系数幅度谱匹配分析的薄层厚度和超声纵波声速
双参数反演
张伟;马志远;赫丽华;高剑英;罗文;林莉;雷明凯
【摘 要】For the overlapped signal and unknow n longitudinal velocity
leading to difficulty in measur‐ing the thickness of thin layers ,a method
bad on ultrasonic reflection coefficient amplitude spectrum (URCAS)
matching analysis technique was propod to characterize thickness and
ultrasonic longitu‐dinal velocity of thin layered structure
simultaneously .An inver algorithm bad on the correlation coefficient
method was introduced .The thickness and velocity were then determined
by matching the theoretical and experimental URCAS within the effective
frequency band of ultrasonic probe and inver‐sing the ultrasonic testing
parameters w hich corresponding to the maximum correlation coefficient .
Experimental validation was conducted on the radar absorbing coatings
(RAC) on the aluminum alloy substrates .The results show that this method
can realize extraction of the ultrasonic parameters from the overlapped
signal effectively ,the relative errors between ultrasonic and micrometer
measurement are in the range of 2 .53%‐3 .72% for the thickness ,and
2 .51%‐3 .75% for the longitudinal velocity .%针对薄层材料超声测厚过程中回
波信号混叠、超声纵波声速未知导致薄层厚度无法测量的问题，提出一种基于声压
反射系数幅度谱（Ultrasonic Reflection Coefficient Amplitude Spectrum ，
URCAS）匹配分析技术同时测量薄层厚度和超声纵波声速的方法。采用相关系数
法对薄层试样实测声压反射系数幅度谱和理论声压反射系数幅度谱在超声检测有效
频带范围内逐一进行匹配分析，通过反演计算得到相关系数最大值点对应的超声检
测参数，最终实现薄层厚度和超声纵波声速的同时表征。利用该方法对铝合金基体
上的雷达吸波涂层（Radar Absorbing Coatings ，RAC）进行实验测试及信号分
析。结果表明：该方法可以有效实现混叠信号中超声特征参量的提取，反演得到吸
波涂层厚度与千分尺测量厚度间相对误差为2．53％～3．72％、纵波声速测量相
对误差为2．51％～3．75％。
【期刊名称】《材料工程》
【年(卷),期】2016(044)010
【总页数】6页(P74-79)
【关键词】薄层材料;声压反射系数幅度谱;相关系数法;厚度;纵波声速;反演
【作 者】张伟;马志远;赫丽华;高剑英;罗文;林莉;雷明凯
【作者单位】大连理工大学材料科学与工程学院，辽宁大连116024;大连理工大
学材料科学与工程学院，辽宁大连116024;北京航空材料研究院，北京100095;
大连理工大学材料科学与工程学院，辽宁大连116024;北京航空材料研究院，北
京100095;大连理工大学材料科学与工程学院，辽宁大连116024;大连理工大学
材料科学与工程学院，辽宁大连116024
【正文语种】中 文
【中图分类】TB52
薄层材料已广泛应用于航空航天等诸多领域，如航空发动机热障涂层、可磨耗封严
涂层、先进战机表面涂装的雷达吸波涂层等。在表征薄层材料性能的各项指标中，神秘莫测
厚度不仅是表征薄层材料自身尺寸的几何参量，而且对薄层材料的性能和使用寿命
也有重要影响[1-6]。例如，在吸波涂层设计和涂装过程中，如果涂层厚度达不到
或超过设计值都无法实现预定的吸波效果，而且厚度过大或不均匀还会降低其结合
强度和使用寿命，引起涂层塑性变差甚至导致局部剥落。因此，采用无损检测技术
对厚度进行准确测量已成为保证薄层材料性能的重要手段。
现有的无损检测技术中，超声脉冲回波技术由于具有原理明确、实现简单、适用范
围广等优点，已广泛应用于层状介质材料的厚度测量之中[7,8]。但对于薄层材料，
超声波在薄层中的传播时间小于其脉冲宽度，薄层表面波和界面波将会发生混叠。
同时，由于受到工艺参数、薄层非均质性等因素的影响，薄层的纵波声速往往为未
知且在不同位置存在一定的波动[9]。这些问题都为薄层厚度的超声检测带来了困
难。Haines等[10]提出利用垂直入射条件下的声压反射系数幅度谱分析方法对铝
基体表面涂覆的环氧树脂层厚度进行测量，实现了混叠信号中超声特征参量的提取，
但该方法将薄层纵波声速设为常数，未考虑非均质性引起的声速波动，难以用于非
均质薄层厚度的测量。赵扬等[11]提出采用高斯-牛顿法对非均质热障涂层厚度和
超声纵波声速进行反演，该方法具有局部收敛速度快、匹配精度高等优点，但对初
始值的选取要求较高，如果初始值选取不当就可能导致只收敛到局部最优解甚至方
程不收敛，无法得到反演结果。
本工作在已有研究的基础上，首先利用声压反射系数幅度谱分析方法实现混叠信号
中超声特征参量的提取，有效解决了薄层中超声信号混叠导致超声特征参量无法提
取的问题。随后，引入统计学中的相关系数匹配法，对薄层试样实测声压反国际爱鸟日
射系数
幅度谱和理论声压反射系数幅度谱矩阵在超声检测有效频带范围内逐一进行匹配分
析，通过计算相关系数最大值点对应的反演参数获得薄层厚度和纵波声速的最佳匹
配结果，实现了超声纵波声速未知及非均质薄层声速波动条件下薄层厚度的测量，
避免了现有高斯-牛顿法由于反演初始值选取不当导致迭代发散、无法获得反演结
果的问题。
声压反射系数幅度谱(URCAS)[11-15]分析方法是一种适用于薄层介质材料超声检
测的信号后处理技术，其核心思想是利用回波信号的声压反射系数幅度谱信息提取
出与被检材料特性相关的超声特征参量，实现材料特性的超声表征。公式(1)为超
声波垂直入射条件下，薄层的声压反射系数幅度谱表达式[11]：
式中：下角1，2，3分别代表探头延迟块、薄层和基底材料；r12，r23为异质材
料界面的声压反射系数；d为薄层厚度；c2和(f)分别为薄层的声速和声衰减系
数。
探头接收到的回波信号实际上是声波经过多次反射和透射后综合叠加的结果。超声
波在薄层结构中多次反射过程中，频率为fn的信号将发生干涉，导致声压反射系
数幅度谱上出现多个极大值或极小值，对(1)式频率f求导，当导数为零时得到
URCAS的一系列极值点
式中fn和n分别为谐振频率和谐振频率阶数。
因此，当薄层纵波声速已知时，利用幅度谱上极值点的位置坐标，通过(2)式即可
测得薄层厚度。但由于受到工艺参数、薄层非均质性等因素的影响，薄层声速通常
为未知量且在不同位置存在一定的波动。如果探头延迟块及基底材料的声速、密度，
薄层的密度和声衰减系数(f)已知，公式(1)就变成了以薄层厚度d和纵波声速c2
为自变量的二元非线性方程，采用合适的方法对此方程进行求解，即可实现d和
c2的双参数同时表征。
本工作从参数反演的角度出发提出采用统计学中的相关系数匹配法对涂层厚度和纵
波声速进行表征。相关系数匹配法是以相关系数为匹配测度判断两组数据相似性的
统计学方法，相关系数越大说明二者的相似度蜜枣的功效
越高。常用的皮尔逊积矩相关系数定
义如下[16]：
式中：X,Y为总体相关系数；X，Y为变量X,Y的标准差；X，Y为X,Y的期
望。基于样本对协方差和标准差进行估计，可以得到样本相关系数
利用式(4)对实测得到的反射系数幅度谱|r(f)|和理论计算得到的反射系数幅度谱
|r(f)|*在-6dB有效频带内进行匹配分析：
式中：N表示时域信号经FFT变换后有效频带范围内的数据点数；下角标i表示第
i个频率值；与*分别为有效频带内实测与理论的声压反射系数幅度谱；与*分别为
有效频带内实测与理论的声压反射系数幅度谱的算术平均值。
通过赋予吸波涂层厚度d和纵波声速c2一系列连续变化值，得到理论声压反射系
数幅度谱矩阵，令其作为母体，将实测的声压反射系数幅度谱与其逐一进行匹配分
析，相关系数矩阵中最大位置对应的d和c2即为被测试样涂层厚度与纵波声速的
最佳反演结果。
2.1 实验样品
实验选用的样品为空气喷涂法制备的雷达吸波涂层试样，吸收剂类型为羰基铁粉，
黏结剂主要成分为环氧树脂，基底材料为铝合金平板。吸波涂层厚度约为0.2mm，
吸波涂层试样的材料属性如表1所示。其中，铝合金基底材料的声速由超声渡越
时间法测得，涂层和铝合金的密度采用阿基米德排水法测定，涂层声衰减系数由实
测数据采用最小二乘法拟合得到。
2.2 实验系统
超声脉冲回波检测系统连接示意图如图1所示。实验中采用延迟块探头作为超声
信号收发装置，探头晶片直径6mm，中心频率5MHz。延迟块材料为环氧树脂，
纵波声速2316.6m/s，密度1045.5kg/m3。使用USIP40超声波探伤仪激发探头，
DPO4032数字示波器以及计算机完成超声信号的观察、数据采集和反演分析工作。
实验过程中首先将延迟块探头连接到探伤仪进行系统与仪器校正。通过DPO4032
数字示波器采集的声源时域波形，其中心频率为5MHz，幅度谱-6dB有效频带宽
为2.95～8.35MHz，如图2所示。将探头耦合到吸波涂层试样上，采集到的涂层
试样时域回波信号和幅度谱如图3所示。
对超声检测实验中采集到的吸波涂层时域回波信号进行频谱分析处理，获得吸波涂
层的实测声压反射系数幅度谱，利用相关系数匹配法对实验得到的反射系数幅度谱
|r(f)|(见图4)和理论计算得到的反射系数幅度谱|r(f)|*在-6dB有效频带内进行匹配
分析，反演结果如图5所示，通过读取图中相关系数rp(d,c2)的最大值位置的坐标，
即可确定出涂层厚度和纵波声速的最佳反演结果。
同理，将延迟块探头耦合到试样的其他位置重复上述步骤，反演结果及测量误差如
表2所示，涂层厚度
反演结果与千分尺测量得到的真实段子手
值间绝对误差为5～8m、相对误差为
2.53%～3.72%，纵波声速与真实值间绝对误差为4.3～6.4m/s、相对误差为
2.51%～3.75%，满足工程检测要求。图6为涂层试样实测与最佳匹配位置理论声
压反射系数幅度谱的对比曲线，可以发现二者在低频段范围内符合较好，但随着频
率的增加二者的偏离程度逐渐增大。分析认为，这是由于涂层理论声压发射系数幅
度谱中密度和衰减系数采用的是涂层试样整体平均测量结果，并未充分考虑吸波涂
层非均质性引起的密度和声衰减系数在不同位置存在波动的问题，导致理论谱线发
生偏离，带来一定的反演误差。
(1)从参数反演的角度，提出采用相关系数匹配法对薄层试样实测和理论声压反射
系数幅度谱逐一进行匹配分析，无需已知薄层超声纵波声速即可实现薄层厚度和纵
波声速的同时表征，解决了薄层厚度和纵波声速均为未知情况下薄层厚度的测量问
题。
(2)利用本方法反演得到的吸波涂层厚度与实际测量值间最大绝对误差为8m、最
大相对误差为3.72%，纵波声速与实测值间最大绝对误差为6.4m/s、最大相对误
差为3.75%。
(3)本方法不仅适用于雷达吸波涂层的厚度测量，还可应用于航空发动机热障涂层、
可磨耗封严涂层等其他非均质薄层材料，具有良好的推广和应用前景。
参考文献
[1] QU C C, LI J, BAI L L, et al. Effects of the thickness of the pre-placed
layer on microstructural evolution and mechanical properties of the lar-
clad coatings[J]. Journal of Alloys and Compounds, 2015, 644: 450-463.
[2] CARVALHO P, SAMPAIO P, AZEVEDO S, et al. Influence of thickness and
coatings morphology in the antimicrobial performance of zinc oxide
coatings[J]. Applied Surface Science, 2014, 307(10): 548-557.
[3] PLOYPECH S, JEARANAISILAWONG P, BOONYONGMANEERAT Y.
Influence of thickness of intermetallic layers on fracture resistance of
galvanized coatings[J]. Surface and Coatings Technology, 2013, 223(6): 1-5.
[4] 赫丽华，刘平桂，王晓红. 雷达吸波涂层的失效行为[J]. 失效分析与预防, 2009,
4(3): 182-187.
HE L H, LIU P G, WANG X H. Degeneration behavior of radar absorbing
coatings[J]. Failure Analysis and Prevention小学音乐教学计划
, 2009, 4(3): 182-187.
[5] 韩萌，黄继华，陈树海. 热障涂层应力与失效机理若干关键问题的研究进展与
评述[J]. 航空材料学报, 2013, 33(5): 83-98.
HAN M, HUANG J H, CHEN S H. Rearch progress and review on key
problems of stress and failure mechanism of thermal barrier coating[J].
Journal of Aeronautical Materials, 2013, 33(5): 83-98.
[6] MORIDI A, AZADI M, FARRAHI G H. Thermo-mechanical stress analysis
of thermal barrier coating system considering thickness and roughness
effects[J]. Surface and Coatings Technology, 2014, 243(4): 91-99.
[7] MAEV R G, SHAO H, MAEVA E Y. Thickness measurement of a curved
multilayered polymer system by using an ultrasonic pul-echo method[J].
Materials Characterization, 1998, 41(2-3): 97-105.
[8] MLLER H, SCHALLER N, EGER T. Ultrasonic determination of thickness
of masticatory mucosa: a methodologic study[J]. Oral Surgery, Oral
Medicine, Oral Pathology, Oral Radiology, and Endodontology, 1999, 88(2):
248-253.
[9] 马志远，罗忠兵，林莉. 基于RVM表征热障涂层孔隙率与孔隙形貌对超声纵波
声速的影响[J]. 材料工程, 2014, (5): 86-90.
MA Z Y, LUO Z B, LIN L. Characterization of ultrasonic longitudinal velocity
evolution in TBCs with different porosity and void morphology bad on
random void model[J]. Journal of Materials Engineering, 2014, (5): 86-90.
[10] HAINES N F, BELL J C, MCINTYRE P J. The application of broadband
ultrasonic spectroscopy to the study of layered media [J]. The Journal of
the Acoustical Society of America, 1978, 64(6): 1645-1651.
[11] ZHAO Y, LIN L, LI X M, et al. Simultaneous determination of the
coating thickness and its longitudinal velocity by ultrasonic nondestructive
method[J]. NDT & E International, 2010, 43(7): 579-585.
[12] LI J C, LIN L, LI X M, et al. Ultrasonic characterization of modified
Cr2O3 coatings by reflection coefficient spectroscopy[J]. Transactions of
Nonferrous Metals Society of China, 2010, 20(3): 418-424.
[13] CONG S, GANG T. Ultrasonic thickness measurement for aluminum
alloy irregular surface parts bad on spectral analysis[J]. Transactions of
Nonferrous Metals Society of China, 2012, 22 (2): 323-328.
[14] LIN L, ZHAO Y, LI X M, et al. Ultrasonic characterization of EB-PVD
thermal barrier coatings irradiated by HIPIB[J]. Key Engineering Materials,
2008, 373-374: 358-362.
[15] MA Z Y, ZHAO Y, LUO Z B, et al. Ultrasonic characterization of
thermally grown oxide in thermal barrier coating by reflection coefficient
amplitude spectrum[J]. Ultrasonics, 2014, 54(4): 1005-国子监是什么机构
1009.
[16] 朱永松，国澄明. 基于相关系数的相关匹配算法的研究[J]. 信号处理,亚洲偶像盛典
2003,
19(6): 531-534.
ZHU Y S, GUO C M. The rearch of correlation matching algorithm bad
on correlation coefficient[J]. Journal of Signal Processing, 2003, 19(6): 531-
534.
