-
对于半透明涂层的传热问题,由于激励方式是脉冲激励,只考虑热流在厚度方向的传热过程。假设材料各向同性且处处均匀,基底厚度无限大,这样将原来的三维双层结构热传导问题简化为一维双层结构热传导问题。 图1所示为脉冲红外加热的半透明涂层几何模型,其中l为涂层的厚度。
当脉冲光激励在涂层表面时,一部分光被涂层表面反射,另一部分光透过涂层上表面,经过涂层吸收后剩余能量辐照在基底表面被基底吸收。为了分析计算方便,假定材料和基底的热力学参数相同,不考虑材料热力学参数随温度变化的情况,仅考虑脉冲光能量和材料的热传导。建立热传导方程:
$$ k \frac{\partial^{2} T}{\partial y^{2}}+q=\rho c \frac{\partial T}{\partial t} $$ (1) 式中:$ T(y, t) $为材料的温度;$ k $为材料的导热系数;$\rho $为材料的密度;$c$为材料的比热容;$ q $为材料的内热源,代表单位时间内脉冲光作用在材料单位体积上产生的能量。表示为:
$$ q(y, t)=p_{0} \cdot f(y) \cdot g(t) $$ (2) 式中:$ p_{0} $为脉冲光的功率密度;$ f(y) $为脉冲光在空间上的分布函数;$ g(t) $为脉冲光在时间上的分布函数。脉冲光在$ y $方向上呈现指数规律衰减,由于半透明涂层的穿透深度较大,不能忽视脉冲光直接对基底表面的热作用,因此脉冲光的时间和空间分布函数分别表示为:
$$ f(y)=\mu {\rm{e}}^{-\mu {y}} \cdot \eta(l-y)+{\rm{e}}^{-\mu l} \cdot \delta(y-l) $$ (3) $$ g(t)=\delta\left(t-t_{0}\right) $$ (4) 式中:$\;\mu $为半透明材料的消光系数;$ \eta(y) $为阶跃函数;$ \delta(y) $为脉冲函数;$ t_{0} $为脉冲光激励时刻。
初始条件:
$$ \left.T(y, t)\right|_{t=0}=0 $$ (5) 边界条件:
$$ \left.T(y, t)\right|_{y=\infty}=0 $$ (6) $$ \left.\frac{\partial T(y, t)}{\partial y}\right|_{y=0}=0 $$ (7) 由于$ t_{0} \rightarrow 0 $,通过Laplace变换:
$y \rightarrow \omega$,$t \rightarrow s$,$ T \rightarrow U $
将公式(1)改写为:
$$ U(\omega, s)=\sqrt{\frac{2}{\pi}} \frac{p_{0}}{\rho c} \frac{{\rm{e}}^{-s t_{0}}}{s+k \omega^{2} / \rho c} \frac{\mu+i \omega {\rm{e}}^{-(\mu+i \omega) l}}{\mu+i \omega} $$ (8) 首先对s做Laplace逆变换:
$$\begin{split} U(\omega, t)=&p_{0} \sqrt{\frac{2}{\pi}} \left[{\rm{exp}}{-\frac{k}{\rho c} \omega^{2}{\left(t-t_{0}\right)}}\right]{ \frac{\mu+i \omega {\rm{e}}^{-(\mu+i \omega) l}}{\mu+i \omega}}=\\ &p_{0} \sqrt{\frac{2}{\pi}} F(\omega) G(\omega) \end{split} $$ (9) 其中
$$ F(\omega)=\exp \left[\frac{-k}{\rho c} \omega^{2}\left(t-t_{0}\right)\right] $$ (10) $$ G(\omega)=\frac{\mu+i \omega {\rm{e}}^{-(\mu+\omega) i}}{\mu+i \omega} $$ (11) $ F(\omega) $的逆变换为:
$$f(y)=\sqrt{2} \lambda {\rm{e}}^{-\lambda^{2} y^{2}},\lambda=\sqrt{\frac{\rho c}{4 k t}} $$ (12) $ G(\omega) $的逆变换为:
$$ g(y)=\sqrt{2 \pi} {\rm{e}}^{-\mu y}[\delta(y-l)-\mu \eta(y-l)+\mu \eta(y)] $$ (13) 将公式(12)和公式(13)卷积后代入公式(9)得到半透明涂层表面温度场:
$$ \begin{split} T(y=0, t)=& p_{0} \sqrt{\frac{2}{\pi}} {\rm{e}}^{-\lambda^{2}l^{2}-\mu l} \cdot\left\{\pi \mu {\rm{e}}^{\left(\frac{2 \lambda^{2}l+\mu}{2 \lambda}\right)^{2}} \cdot\left[{{\rm{erf}}}\left(\lambda l+\frac{\mu}{2 \lambda}\right)-\right.\right. \\ &\left.\left.{{\rm{erf}}}\left(\frac{\mu}{2 \lambda}\right)\right]+2 \sqrt{\pi} \lambda\right\} \end{split} $$ (14) 假定闪光灯功率500 W,半透明材料为环氧树脂,其消光系数为0.5 mm−1,得到图2所示的不同厚度的半透明涂层表面温度时间变化曲线。
图 2 双对数坐标下不同涂层厚度样品的表面温度变化解析结果
Figure 2. Analytical surface temperature of samples with different coating thickness vs time in double logarithm scale
从图2中可以看出,因为脉冲光对不同涂层厚度的样品表面激励能量相同,不同涂层厚度的样品表面温度在初始时刻温度一致。由于脉冲光同时激励交界面,交界面具有更高的温度,热流同时向涂层上表面和基底下表面传导,样品表面温度持续升高,直到热流到达涂层表面,此时表面温度达到峰值。由于基底半无限大,热流持续向下传导,样品表面温度持续下降,最终不同厚度的样品表面温度以同样的速率降温,降温曲线斜率近似为−0.5,说明此时模型和半无限大脉冲激励一维热传导模型的降温过程近似。
图3显示了表面温度峰值时刻与涂层厚度在对数坐标的关系。结果显示涂层厚度与表面温度峰值时刻在对数坐标上呈现线性关系。材料热力学参数固定的情况下,热波在固体中一维热传导的距离和时间呈现指数关系。表面温度峰值出现的原因是交界面向涂层上表面的热传导作用,因此涂层厚度和表面温度在对数坐标上呈现线性关系。
-
图4为半透明涂层的红外加热模型。其中a为模型长度,b和l分别为基底和涂层的厚度。
-
当脉冲光激励在涂层表面后,一部分光被涂层上表面反射,另一部分光透过涂层上表面,经过涂层吸收后剩余能量辐照在基底表面,被基底吸收。涂层和基底吸收能量后温度迅速升高。为了分析方便,不考虑材料热力学参数随温度变化的情况,仅考虑脉冲光能量和材料的热传导。因此,热传导方程表示为:
$$ {{k}}_{{i}}\frac{{\partial }^{2}{{T}}_{{i}}({y},{t})}{\partial {{y}}^{2}}+{{q}}_{{i}}={{\rho }}_{{i}}{{c}}_{{i}}\frac{\partial {{T}}_{{i}}\left({y},{t}\right)}{\partial {t}}({i}={1,2}) $$ (15) 式中:$ T_{1}(y, t) $和$ {{T}}_{2}({y,t}) $分别表示t时刻时涂层和基底的温度场; $ k_{1} $ 和$ k_{2} $分别表示涂层和基底的导热系数;$ \;\rho_{1} $和$\; \rho_{2} $分别表示涂层和基底的密度;$ c_{1} $和$ c_{2} $分别表示涂层和基底的比热容;$ q_{1} $和$ q_{2} $为涂层和基底的内热源。
在闪光灯辐照下,半透明涂层吸收闪光灯能量转化成热能引起涂层温度迅速升高,剩余能量被基底金属材料吸收转化成热能引起表层温度迅速升高。金属材料的表层厚度一般在纳米级,因此可以将闪光灯对金属基底的热作用看作边界条件,内热源可以表示为:
$$ q_{1}=p_{0} \cdot f(y) \cdot g(t) $$ (16) $$ {{q}}_{2}=0 $$ (17) 式中:$ p_{0} $为脉冲光功率密度;$ f(y) $为闪光灯在涂层内部的空间分布函数;$ g(t) $为闪光灯的时间分布函数。分别表示为:
$$ f(y)=\mu \cdot {\rm{e}}^{-\mu {y}} $$ (18) $$ g(t)=\varepsilon(t) \cdot \varepsilon\left(t_{0}-t\right) $$ (19) 式中:$\;\mu $为材料的消光系数;$\varepsilon (t)$为关于时间的阶跃函数;${t_0}$为闪光灯脉冲持续时间。在涂层与基底的交界面满足连续条件:
$$ {\left.-{{k}}_{1}\frac{\partial {{T}}_{1}\left({y},{t}\right)}{\partial {y}}\right|}_{{y}={l}}+{{p}}_{0}\cdot {{{\rm{e}}}}^{-{\mu }{l}}\cdot {g}\left({t}\right)={\left.-{{k}}_{2}\frac{\partial {{T}}_{2}\left({y},{t}\right)}{\partial {y}}\right|}_{{y}={l}} $$ (20) $$ T_{1}(l, t)=T_{2}(l, t) $$ (21) 模型中基底下表面的边界条件按照绝热条件处理:
$$ \left.\frac{\partial T(y, t)}{\partial y}\right|_{y=l+b}=0 $$ (22) 温度场的初始条件:
$$ \left.T(y, t)\right|_{t=0}=293.15 \;{\rm{K}} $$ (23) -
由于环氧树脂具备高透明性,将其与碳粉混合作为半透明涂层的材料,结构钢作为基底。
为了确定涂层材料的消光系数$\;\mu $,制作了不同厚度的半透明试片,通过测量不同厚度试片的透射光强,对不同厚度的透射光强进行拟合,效果如图5所示。得到半透明涂层厚度与光强的关系如下:
$$ I=I_{0} {\rm{e}}^{-0.496\;5l} $$ (24) 故半透明涂层的消光系数$\;\mu = 0.496\;5$ mm−1,材料的其他热力学参数如表1所示。
表 1 模型使用的热力学参数
Table 1. Parameters used in the model
Density/
$\mathrm{k}\mathrm{g}\cdot{\mathrm{m} }^{-3}$Heat
capacity/
$\mathrm{J}\cdot(\mathrm{k}\mathrm{g}\cdot \mathrm{K}$)−1Thermal
conductivity/
$\mathrm{W}\cdot(\mathrm{m}\cdot \mathrm{K}$)−1Coating 1200 1000 0.2 Substrate 7850 475 80 -
在仿真中,涂层厚度范围为50~200 µm,厚度步长25 µm, 基底厚度5 mm,脉冲峰值功率104 W,闪光灯脉冲持续时间为3 ms。样品表面温度T在对数坐标的时间变化如图6所示。
图 6 双对数坐标下不同涂层厚度样品的温度变化仿真结果
Figure 6. Simulation surface temperature of samples with different coating thickness vs time in double logarithm scale
结果表明,虽然涂层厚度不同,但由于涂层表面吸收的能量相等,涂层的表面温度在一开始基本相同。由于热传导作用,金属基底表层吸收的能量传递到涂层表面,导致涂层表面温度升高达到峰值温度。随着时间推移,涂层表面温度逐渐降低,直到试件整体达到热平衡状态。
分析图6中温度信号变化规律可知,涂层表面温度到达峰值的时刻与涂层的厚度有关。图7显示了表面温度峰值时刻与涂层厚度在对数坐标的关系。图中的线性关系表明,该方法可应用于半透明涂层的厚度测量。
-
将环氧树脂和碳粉混合作为半透明涂层的材料,将其浇注在深色钢板的表面上,钢板右侧向上倾斜0.5 mm。因此,涂层的厚度从左到右逐渐减少。值得注意的是,与常规红外无损检测不同,样品表面不经过喷黑漆处理,图8为厚度连续变化的涂层试件示意图。
基底长200 mm,宽100 mm,厚度5 mm。在浇注涂层前后,分别用千分尺多点测量试件的厚度,取其差值作为该点涂层厚度。
-
为了提升脉冲能量的均匀性,使用两个U型闪光灯管作为激励热源,闪光灯能量为9.6 kJ,脉宽为3 ms。在U型闪光灯管前方放置磨砂毛玻璃片,进一步匀化脉冲激励能量。使用非制冷热像仪(Flir A655 sc)采集试件表面温度,其等效噪声温差(NETD)为30 mK,图像分辨率为640×120,采集帧频为200 Hz。计算机控制闪光灯电源及同步触发器,在红外热像仪单帧图像采集结束后,电源驱动闪光灯对样品表面进行加热。红外热像仪采集样品表面温度信息,并将其传输回计算机进行数据处理分析。选取部分热成像序列如图9所示。
图10为在对数坐标中绘制了不同涂层厚度的表面温度与时间的关系及拟合曲线。结果表明,不同涂层厚度下的表面温度变化趋势与模拟结果基本一致。表面温度的对数峰值时刻与涂层厚度之间存在单调关系。
图 10 双对数坐标下不同涂层厚度样品的温度变化实验结果
Figure 10. Experimental surface temperature of samples with different coating thickness vs time in double logarithm scale
实验结果与理论结果的数据差异在于:实验结果中的表面温度峰值时刻与涂层厚度在对数坐标上的线性关系存在误差。造成该现象的主要原因有以下方面:(1)热像仪帧频导致温度数据丢失,对于较薄涂层,表面温度的峰值出现的较早,温度上升较快。热像仪采集到的表面温度峰值时刻存在误差,曲线拟合在一定程度上降低了表面温度峰值时刻的误差。(2)数值计算中将半透明涂层样品简化成一维热传导模型,其材料均匀各向同性,但实验上的样品厚度连续变化,不如理论样品理想。(3)理论计算中涂层表面均匀加热,实验时闪光灯能量的分布并不如理论理想。
在对数坐标中,表面温度的峰值时刻与涂层厚度之间的关系曲线如图11所示。可以看出,它们之间存在线性关系,与图7所示的模拟结果一致,这些数据点与拟合直线的偏差在5%以内。
同一样本在另外四个区域取样。根据图11的校准曲线计算涂层厚度,结果如表2所示,厚度误差在5%范围内。
表 2 涂层厚度检测结果
Table 2. Measurement result of coating thickness
Thickness by
micrometer/
μmThickness by
thermography/
μmDifference 1 103 99.7 3.3% 2 136 139.2 2.3% 3 158 160.6 1.6% 4 197 189.4 4.0% 图12为试件厚度分布的图像。厚度从右侧较薄的一侧逐渐过渡为左侧较厚的一侧。由于红外热像仪帧频和闪光灯能量强度的限制,检测区域仅包括试样的右半部分,涂层厚度在80~220 µm之间。
Measurement of translucent coating thickness based on pulsed thermography
-
摘要: 为了快速、准确地测量半透明涂层的厚度,提出了一种基于脉冲红外热波的测量方法。建立了半透明涂层半无限大脉冲热传导简化理论模型和半透明涂层的脉冲红外加热双层物理模型,理论分析和数值计算结果表明:半透明涂层的厚度与表面温度的峰值时刻在对数坐标上呈现线性关系,利用这种线性关系可以直接测量半透明涂层的厚度,而不再需要在样品表面喷涂黑漆以避免半透明性的影响。实验上建立脉冲红外热波系统并制作了厚度连续变化的半透明涂层试件,得到的厚度误差小于5%。结果显示该技术具有快速和非接触式测量半透明涂层厚度的潜力。Abstract: To measure the thickness of translucent coating quickly and accurately, the method based on pulse thermography was proposed. A simplified theoretical semi-infinite model of pulsed infrared heating conduction for translucent coating and a two-layer physical model of pulsed infrared heating for translucent coatings were established. The theoretical analysis and numerical results show that the thickness of the translucent coating is linearly proportional to the peak time of the surface temperature in logarithmic coordinates, through which the thickness of the translucent coating can be derived without spraying black ink on the surface of the sample to avoid the effect of translucency. A translucent coating specimen with continuously varying thickness was fabricated and tested with a flash thermographic system. The result of the measured thicknesses is within 5% of the actual values. It demonstrates that the technique has the potential for fast and non-contact measurement of the thickness for translucent coatings.
-
表 1 模型使用的热力学参数
Table 1. Parameters used in the model
Density/
$\mathrm{k}\mathrm{g}\cdot{\mathrm{m} }^{-3}$Heat
capacity/
$\mathrm{J}\cdot(\mathrm{k}\mathrm{g}\cdot \mathrm{K}$)−1Thermal
conductivity/
$\mathrm{W}\cdot(\mathrm{m}\cdot \mathrm{K}$)−1Coating 1200 1000 0.2 Substrate 7850 475 80 表 2 涂层厚度检测结果
Table 2. Measurement result of coating thickness
Thickness by
micrometer/
μmThickness by
thermography/
μmDifference 1 103 99.7 3.3% 2 136 139.2 2.3% 3 158 160.6 1.6% 4 197 189.4 4.0% -
[1] 索会迎. 超声波无损检测技术应用研究[D]. 南京邮电大学, 2012. Suo Huiying. The application and research of ultrasonic non-destructive testing technology [D]. Nanjing: Najing University of Posts and Telecommunications, 2012. (in Chinese) [2] 姚志伟, 魏永庚, 石勇. 基于超声波测量油膜厚度方法的研究[J]. 黑龙江大学工程学报, 2021, 12(01): 54-61. Yao Zhiwei, Wei Yonggeng, Shi Yong. Research on the method of measuring oil film thickness based on ultrasonic measurement [J]. Journal of Engineering of Heilongjiang University, 2021, 12(1): 54-61. (in Chinese) [3] 潘翔, 张永康, 常国强. 超声波测厚仪校准结果超差问题的分析[J]. 计量与测试技术, 2020, 47(10): 89-91. Pan Xiang, Zhang Yongkang, Chang Guoqiang. Analysis of the out-of-tolerance of calibration results of ultrasonic thickness gauge [J]. Metrology & Measurement Technique, 2020, 47(10): 89-91. (in Chinese) [4] 孙吉伟, 冯辅周, 闵庆旭, 徐超, 朱俊臻. 涡流脉冲热像技术中检测条件的粒子群优化[J]. 红外与激光工程, 2018, 47(08): 240-246. doi: 10.3788/IRLA201847.0818005 Sun Jiwei, Feng Fuzhou, Min Qingxu, et al. Optimization of eddy current pulsed thermography detection condition using particle swarm optimization [J]. Infrared and Laser Engineering, 2018, 47(8): 0818005. (in Chinese) doi: 10.3788/IRLA201847.0818005 [5] 王安泉, 刘禹, 郝宪锋, 戴永寿, 闫泰松, 姜朝. 基于功率谱密度峰值的大提离下脉冲涡流测厚方法[J]. 无损检测, 2021, 43(10): 29-33. doi: 10.11973/wsjc202110006 Wang Anquan, Liu Yu, Hao Xianfeng, et al. Pulsed eddy current thickness measurement method under large lift based on power spectral density peak [J]. NDT, 2021, 43(10): 29-33. (in Chinese) doi: 10.11973/wsjc202110006 [6] 谢建红, 石剑, 习小文, 肖力伟, 宋凯. 基于BP算法的热障涂层厚度反演研究[J]. 兵器装备工程学报, 2021, 42(4): 228-232+243. Xie Jianhong, Shi Jian, Xi Xiaowen, et al. Research on inversion of thermal barrier coating thickness based on BP neural network [J]. Journal of Ordnance Equipment Engineering, 2021, 42(4): 228-232, 243. (in Chinese) [7] Ptaszek G, Cawley P, Almond D, et al. Transient thermography testing of unpainted thermal barrier coating (TBC) systems [J]. NDTE International, 2013, 59: 48-56. doi: 10.1016/j.ndteint.2013.05.001 [8] Shepard S M, Hou Y L, Lhota J R, et al. Thermographic measurement of thermal barrier coating thickness [C]//Proceedings of SPIE, 2005, 5782: 407-410. [9] 李艳红, 赵跃进, 张存林. 半透明材料涂层的红外检测[J]. 光学精密工程, 2009, 17(07): 1502-1506. Li Yanhong, Zhao Yuejin, Zhang Cunlin. Infrared testing for translucent coatings [J]. Optics and Precision Engineering, 2009, 17(7): 1502-1506. (in Chinese) [10] Zhao S , Wang H , Wu N , et al. Nondestructive testing of the fatigue properties of air plasma sprayed thermal barrier coatings by pulsed thermography[J]. Russian Journal of Nondestructive Testing, 2015, 51: 445-456. [11] 刘颖韬, 唐佳, 郭兴旺, 杨党纲, 岳翔. 半透明性对闪光灯激励红外热像检测的影响[J]. 激光与红外, 2017, 47(10): 1264-1270. doi: 10.3969/j.issn.1001-5078.2017.10.014 Liu Yingtao, Tang Jia, Guo Xingwang, et al. Effects of composites translucence on flash in frared thermal detection [J]. Laser & Infrared, 2017, 47(10): 1264-1270. (in Chinese) doi: 10.3969/j.issn.1001-5078.2017.10.014 [12] 陈栋, 郭兴旺, 刘颖韬. 半透明复合材料脉冲热像检测的有限元仿真分析[J]. 无损检测, 2016, 38(10): 42-47+52. Chen Dong, Guo Xingwang, Liu Yingtao. Modeling of infrared pulsed thermography on translucent composites[J] Nondestructive Testing, 2016, 38(10): 42-47, 52. (in Chinese) [13] 郭兴旺, 管和清, 刘颖韬, 唐佳. 半透明材料红外热像检测的光谱特性和光源选择[J]. 红外与激光工程, 2017, 46(01): 59-67. Guo Xingwang, Guan Heqing, Liu Yingtao, et al. Spectrum characteristics and light source selection for infrared thermal imaging testing of semitransparent materials [J]. Infrared and Laser Engineering, 2017, 46(1): 0104001. (in Chinese) [14] Tang Q, Liu J, Dai J, et al. Theoretical and experimental study on thermal barrier coating (TBC) uneven thickness detection using pulsed infrared thermography technology [J]. Applied Thermal Engineering, 2017, 114: 770-775. doi: 10.1016/j.applthermaleng.2016.12.032 [15] Unnikrishnakurup S, Dash J, Ray S, et al. Nondestructive evaluation of thermal barrier coating thickness degradation using pulsed IR thermography and THz-TDS measurements: A comparative study [J]. NDT & E International, 2020, 116: 102367. [16] Shrestha R, Kim W. Evaluation of coating thickness by thermal wave imaging: A comparative study of pulsed and lock-in infrared thermography - Part I: Simulation [J]. Infrared Physics & Technology, 2017, 83: 124-131. [17] Shrestha R, Kim W. Evaluation of coating thickness by thermal wave imaging: A comparative study of pulsed and lock-in infrared thermography – Part II: Experimental investigation [J]. Infrared Physics & Technology, 2018, 92: 24-29. [18] Caruso G, Paoloni S, Orazi N, et al. Quantitative evaluations by infrared thermography in optically semi-transparent paper-based artefacts [J]. Measurement, 2019, 143: 258-266. [19] He Lengyun. Research on coating thickness measurement based on infrared nondestructive testing technology[D]. Chengdu: University of Electronic Science and Technology of China, 2018. (in Chinese)