留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

脉冲红外热波用于半透明涂层厚度测量

何涛 沈中华 陈飞 陈力

何涛, 沈中华, 陈飞, 陈力. 脉冲红外热波用于半透明涂层厚度测量[J]. 红外与激光工程, 2022, 51(2): 20210890. doi: 10.3788/IRLA20210890
引用本文: 何涛, 沈中华, 陈飞, 陈力. 脉冲红外热波用于半透明涂层厚度测量[J]. 红外与激光工程, 2022, 51(2): 20210890. doi: 10.3788/IRLA20210890
He Tao, Shen Zhonghua, Chen Fei, Chen Li. Measurement of translucent coating thickness based on pulsed thermography[J]. Infrared and Laser Engineering, 2022, 51(2): 20210890. doi: 10.3788/IRLA20210890
Citation: He Tao, Shen Zhonghua, Chen Fei, Chen Li. Measurement of translucent coating thickness based on pulsed thermography[J]. Infrared and Laser Engineering, 2022, 51(2): 20210890. doi: 10.3788/IRLA20210890

脉冲红外热波用于半透明涂层厚度测量

doi: 10.3788/IRLA20210890
基金项目: 国家自然科学基金(61975080)
详细信息
    作者简介:

    何涛,男,硕士生,主要从事无损检测方面的研究

    沈中华,女,教授,博士生导师,主要从事激光超声技术及应用、激光与材料相互作用及应用等方面的研究

  • 中图分类号: TN219

Measurement of translucent coating thickness based on pulsed thermography

  • 摘要: 为了快速、准确地测量半透明涂层的厚度,提出了一种基于脉冲红外热波的测量方法。建立了半透明涂层半无限大脉冲热传导简化理论模型和半透明涂层的脉冲红外加热双层物理模型,理论分析和数值计算结果表明:半透明涂层的厚度与表面温度的峰值时刻在对数坐标上呈现线性关系,利用这种线性关系可以直接测量半透明涂层的厚度,而不再需要在样品表面喷涂黑漆以避免半透明性的影响。实验上建立脉冲红外热波系统并制作了厚度连续变化的半透明涂层试件,得到的厚度误差小于5%。结果显示该技术具有快速和非接触式测量半透明涂层厚度的潜力。
  • 图  1  理论分析模型示意图

    Figure  1.  Schematic diagram of theoretical analysis model

    图  2  双对数坐标下不同涂层厚度样品的表面温度变化解析结果

    Figure  2.  Analytical surface temperature of samples with different coating thickness vs time in double logarithm scale

    图  3  涂层厚度-解析峰值时刻对数图

    Figure  3.  Coating thickness vs peak time of surface temperature in logarithm scale

    图  4  数值计算模型示意图

    Figure  4.  Schematic diagram of numerical calculation model

    图  5  涂层厚度-透过率关系

    Figure  5.  Coating thickness vs transmittance

    图  6  双对数坐标下不同涂层厚度样品的温度变化仿真结果

    Figure  6.  Simulation surface temperature of samples with different coating thickness vs time in double logarithm scale

    图  7  涂层厚度-模拟峰值时刻对数关系

    Figure  7.  Coating thickness vs simulation peak time of surface temperature in logarithm scale

    图  8  涂层试件示意图

    Figure  8.  Schematic diagram of the coating specimen

    图  9  部分红外图像序列

    Figure  9.  Selected frames from the IR image sequence

    图  10  双对数坐标下不同涂层厚度样品的温度变化实验结果

    Figure  10.  Experimental surface temperature of samples with different coating thickness vs time in double logarithm scale

    图  11  涂层厚度-实验峰值时刻对数图

    Figure  11.  Thickness vs experimental peak time in double logarithm scale

    图  12  试件厚度检测结果分布图

    Figure  12.  Test result of thickness distribution of the specimen

    表  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}$)−1
    Thermal
    conductivity/
    $\mathrm{W}\cdot(\mathrm{m}\cdot \mathrm{K}$)−1
    Coating120010000.2
    Substrate785047580
    下载: 导出CSV

    表  2  涂层厚度检测结果

    Table  2.   Measurement result of coating thickness

    Thickness by
    micrometer/
    μm
    Thickness by
    thermography/
    μm
    Difference
    110399.73.3%
    2136139.22.3%
    3158160.61.6%
    4197189.44.0%
    下载: 导出CSV
  • [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)
  • [1] 殷博, 薛常喜, 李闯.  碳化钨合金磨削仿真与工艺优化方法 . 红外与激光工程, 2023, 52(7): 20220898-1-20220898-10. doi: 10.3788/IRLA20220898
    [2] 朱强, 周维虎, 陈晓梅, 李冠楠, 石俊凯.  复杂衬底下的涂层厚度定量检测方法研究 . 红外与激光工程, 2022, 51(12): 20220156-1-20220156-10. doi: 10.3788/IRLA20220156
    [3] 张鹏辉, 赵扬, 李鹏, 周志权, 白雪, 马健.  基于有限元法的激光声磁检测系统优化研究 . 红外与激光工程, 2022, 51(7): 20210533-1-20210533-9. doi: 10.3788/IRLA20210533
    [4] 孙权, 莫德锋, 刘大福, 龚海梅.  深低温大功率电阻阵列封装结构研究 . 红外与激光工程, 2022, 51(8): 20210721-1-20210721-9. doi: 10.3788/IRLA20210721
    [5] 李鹏, 赵扬, 周志权, 张鹏辉, 白雪, 马健.  一种跨介质的空中-水下激光致声探测技术研究 . 红外与激光工程, 2021, 50(5): 20200310-1-20200310-9. doi: 10.3788/IRLA20200310
    [6] 李晟, 范斌, 王伟刚, 李康.  深低温SiC空间反射镜背部与侧面支撑结构对比 . 红外与激光工程, 2020, 49(2): 0214003-0214003. doi: 10.3788/IRLA202049.0214003
    [7] 汪奎, 辛宏伟, 徐宏, 任天赐.  空间相机快速反射镜的结构轻量化设计 . 红外与激光工程, 2019, 48(4): 418001-0418001(7). doi: 10.3788/IRLA201948.0418001
    [8] 杨正伟, 谢星宇, 李胤, 张炜, 田干.  激光扫描热成像无损检测关键参数影响分析 . 红外与激光工程, 2019, 48(11): 1105008-1105008(11). doi: 10.3788/IRLA201948.1105008
    [9] 郭兴旺, 管和清, 刘颖韬, 唐佳.  半透明材料红外热像检测的光谱特性和光源选择 . 红外与激光工程, 2017, 46(1): 104001-0104001(9). doi: 10.3788/IRLA201746.0104001
    [10] 郭伟, 董丽虹, 王海斗, 徐雅薇, 徐滨士.  基于小波分解的热波相位特征提取及喷涂层厚度评价 . 红外与激光工程, 2017, 46(9): 904003-0904003(7). doi: 10.3788/IRLA201746.0904003
    [11] 李永君, 肖俊峰, 朱立春, 张炯, 高斯峰, 唐文书, 南晴.  热障涂层厚度激光透射法红外热波检测技术研究 . 红外与激光工程, 2017, 46(7): 704003-0704003(5). doi: 10.3788/IRLA201746.0704003
    [12] Sekou Singare, 陈盛贵, 钟欢欢.  激光透射焊接聚碳酸酯的有限元分析 . 红外与激光工程, 2016, 45(2): 206005-0206005(6). doi: 10.3788/IRLA201645.0206005
    [13] 张耀平, 樊峻棋, 龙国云.  变形镜在激光辐照下热畸变有限元模拟 . 红外与激光工程, 2016, 45(11): 1136002-1136002(5). doi: 10.3788/IRLA201645.1136002
    [14] 肖冬明, 何宽芳, 王迪.  基于多层有限元模型的激光选区熔化多层瞬态温度场演化规律研究 . 红外与激光工程, 2015, 44(9): 2672-2678.
    [15] 许兆美, 汪通悦, 裴旭, 蒋素琴, 李伯奎, 王庆安, 洪宗海.  Al2O3陶瓷激光多道铣削温度场有限元模拟 . 红外与激光工程, 2015, 44(2): 477-481.
    [16] 陈林, 杨立, 范春利, 石宏臣, 赵小龙.  基于相位的热障涂层厚度及其脱粘缺陷红外定量识别 . 红外与激光工程, 2015, 44(7): 2050-2056.
    [17] 赵宏超, 张景旭, 杨飞.  TMT 三镜系统tilt 轴轴承方案研究 . 红外与激光工程, 2015, 44(1): 122-126.
    [18] 唐庆菊, 王扬, 刘俊岩, 秦雷, 龚金龙, 齐立涛.  耐热合金涂层结构板材内部缺陷的脉冲红外热成像检测 . 红外与激光工程, 2013, 42(7): 1685-1690.
    [19] 朱敏, 陈宇, 杨春玲.  红外诱饵弹干扰特性有限元建模 . 红外与激光工程, 2013, 42(8): 1979-1986.
    [20] 周于鸣, 赵野, 王海超, 闻广泉.  大口径轻质反射镜光轴水平卸载支撑方法 . 红外与激光工程, 2013, 42(5): 1285-1290.
  • 加载中
图(12) / 表(2)
计量
  • 文章访问数:  335
  • HTML全文浏览量:  171
  • PDF下载量:  35
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-11-24
  • 修回日期:  2022-01-29
  • 网络出版日期:  2022-03-04
  • 刊出日期:  2022-02-28

脉冲红外热波用于半透明涂层厚度测量

doi: 10.3788/IRLA20210890
    作者简介:

    何涛,男,硕士生,主要从事无损检测方面的研究

    沈中华,女,教授,博士生导师,主要从事激光超声技术及应用、激光与材料相互作用及应用等方面的研究

基金项目:  国家自然科学基金(61975080)
  • 中图分类号: TN219

摘要: 为了快速、准确地测量半透明涂层的厚度,提出了一种基于脉冲红外热波的测量方法。建立了半透明涂层半无限大脉冲热传导简化理论模型和半透明涂层的脉冲红外加热双层物理模型,理论分析和数值计算结果表明:半透明涂层的厚度与表面温度的峰值时刻在对数坐标上呈现线性关系,利用这种线性关系可以直接测量半透明涂层的厚度,而不再需要在样品表面喷涂黑漆以避免半透明性的影响。实验上建立脉冲红外热波系统并制作了厚度连续变化的半透明涂层试件,得到的厚度误差小于5%。结果显示该技术具有快速和非接触式测量半透明涂层厚度的潜力。

English Abstract

    • 随着现代工业的快速发展,涂层和薄膜扮演着越来越重要的角色。涂层厚度直接影响涂层的光学性能、热学性能等,因此,现代工业对涂层厚度的测量提出了更高的要求。涂层厚度测量技术可以分为直接测量和间接测量。常见的直接测量技术有金相显微镜法、千分尺法和探针法等。这些测量技术大多具有破坏性,且检测精度有限。间接测量是通过测量材料的声、光、热、电、磁等物理性质的变化,将相关的物理量通过计算转化为涂层厚度。超声检测的优势在于穿透力强、方向性好,但需要在待测件表面涂抹耦合剂,且材料中的分层、气孔等缺陷也会导致测量误差,在工业检测中会受到一定限制[1-3];涡流检测的优势在于检测速度快,无需耦合剂,但要求待测件是导电金属材料[4-6],且并不具备大面积检测的能力,不适用于现代工业大面积的涂层厚度检测;红外检测具备检测速度快、非接触、检测面积大等优势,能够胜任现代工业对涂层厚度的检测工作。

      目前对于大多数红外检测技术,待测件大多是不透明的,激励能量被材料吸收引起表层温度迅速上升。然而,对于一些透明度较高的涂层材料,激励能量并不仅仅作用在材料表层,而是对涂层整体进行加热,因此,广泛使用的表面加热模型已不再适用。

      Ptaszek 进行了关于热障涂层缺陷检测的理论和实验分析,提出使用长波红外相机可以消除半透明性对温度采集的不利影响[7]。Shepard等人研究了基于脉冲热成像的半透明热障涂层厚度检测技术,认为传统的TSR方法在一定程度上仍然可以用于检测涂层厚度[8]。李艳红等人通过在样品表面喷涂黑色涂料来消除涂层的半透明性,并使用二阶导数的峰值时刻来确定涂层厚度[9]。赵世斌等人在TBC表面喷涂胶体石墨,通过脉冲热成像检测热障涂层的疲劳缺陷 [10]。刘颖韬等人建立半透明试件内部缺陷的脉冲红外物理加热模型,结果显示试件的半透明性会降低表面温度信号或改变信号变化趋势,这可能会影响检测能力 [11-12]。郭兴旺等人通过有限元仿真进行闪光灯色温对缺陷检测效果影响的研究,结果显示低色温和高色温闪光灯对半透明材料红外检测更有利[13]。唐庆菊等人建立了不透明涂层的脉冲红外加热双层热传导模型,结果显示提高脉冲能量有助于涂层厚度不均匀性检测,提高采样频率有利于采集最大对比度红外图像[14]。Sreedhar等人基于脉冲相位热成像和太赫兹时域光谱进行热障涂层厚度退化对比研究,结果显示脉冲法检测误差小于10%,有助于飞机发动机等部件的预防性维修措施[15]。Shrestha等人通过有限元仿真和实验进行脉冲相位法和锁相法的涂层测厚对比研究,结果显示脉冲相位法的检测速度和检测精度都更胜一筹[16-17]。Caruso等人基于脉冲红外热成像对半透明纸制品的表面温度场进行研究,有限元数值模拟和实验结果显示样品表面温度变化与半透明纸品厚度呈现相关性[18]。何棱云等人基于长脉冲红外加热进行环氧树脂涂层测厚研究,提出一种基于 k 范围温度变化斜率曲线的涂层厚度检测方法,实验结果显示该方法对于环氧树脂涂层的检测误差小于10%[19]

      上述大多数研究都是用不透明材料覆盖待测件表面以避免涂层透明性的干扰,或者忽略涂层的透明性影响。基于红外无损检测直接对半透明涂层进行测厚的研究较少,且检测精度不高。因此,文中建立了半透明涂层的半无限大脉冲红外加热理论模型和脉冲红外加热双层物理模型,以中等吸光度和较薄涂层作为重点,分析不同厚度下的表面温度场,制作样品并通过实验验证了理论模拟结果。

    • 对于半透明涂层的传热问题,由于激励方式是脉冲激励,只考虑热流在厚度方向的传热过程。假设材料各向同性且处处均匀,基底厚度无限大,这样将原来的三维双层结构热传导问题简化为一维双层结构热传导问题。 图1所示为脉冲红外加热的半透明涂层几何模型,其中l为涂层的厚度。

      图  1  理论分析模型示意图

      Figure 1.  Schematic diagram of theoretical analysis model

      当脉冲光激励在涂层表面时,一部分光被涂层表面反射,另一部分光透过涂层上表面,经过涂层吸收后剩余能量辐照在基底表面被基底吸收。为了分析计算方便,假定材料和基底的热力学参数相同,不考虑材料热力学参数随温度变化的情况,仅考虑脉冲光能量和材料的热传导。建立热传导方程:

      $$ 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显示了表面温度峰值时刻与涂层厚度在对数坐标的关系。结果显示涂层厚度与表面温度峰值时刻在对数坐标上呈现线性关系。材料热力学参数固定的情况下,热波在固体中一维热传导的距离和时间呈现指数关系。表面温度峰值出现的原因是交界面向涂层上表面的热传导作用,因此涂层厚度和表面温度在对数坐标上呈现线性关系。

      图  3  涂层厚度-解析峰值时刻对数图

      Figure 3.  Coating thickness vs peak time of surface temperature in logarithm scale

    • 图4为半透明涂层的红外加热模型。其中a为模型长度,bl分别为基底和涂层的厚度。

      图  4  数值计算模型示意图

      Figure 4.  Schematic diagram of numerical calculation model

    • 当脉冲光激励在涂层表面后,一部分光被涂层上表面反射,另一部分光透过涂层上表面,经过涂层吸收后剩余能量辐照在基底表面,被基底吸收。涂层和基底吸收能量后温度迅速升高。为了分析方便,不考虑材料热力学参数随温度变化的情况,仅考虑脉冲光能量和材料的热传导。因此,热传导方程表示为:

      $$ {{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)

      图  5  涂层厚度-透过率关系

      Figure 5.  Coating thickness vs transmittance

      故半透明涂层的消光系数$\;\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}$)−1
      Thermal
      conductivity/
      $\mathrm{W}\cdot(\mathrm{m}\cdot \mathrm{K}$)−1
      Coating120010000.2
      Substrate785047580
    • 在仿真中,涂层厚度范围为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显示了表面温度峰值时刻与涂层厚度在对数坐标的关系。图中的线性关系表明,该方法可应用于半透明涂层的厚度测量。

      图  7  涂层厚度-模拟峰值时刻对数关系

      Figure 7.  Coating thickness vs simulation peak time of surface temperature in logarithm scale

    • 将环氧树脂和碳粉混合作为半透明涂层的材料,将其浇注在深色钢板的表面上,钢板右侧向上倾斜0.5 mm。因此,涂层的厚度从左到右逐渐减少。值得注意的是,与常规红外无损检测不同,样品表面不经过喷黑漆处理,图8为厚度连续变化的涂层试件示意图。

      图  8  涂层试件示意图

      Figure 8.  Schematic diagram of the coating specimen

      基底长200 mm,宽100 mm,厚度5 mm。在浇注涂层前后,分别用千分尺多点测量试件的厚度,取其差值作为该点涂层厚度。

    • 为了提升脉冲能量的均匀性,使用两个U型闪光灯管作为激励热源,闪光灯能量为9.6 kJ,脉宽为3 ms。在U型闪光灯管前方放置磨砂毛玻璃片,进一步匀化脉冲激励能量。使用非制冷热像仪(Flir A655 sc)采集试件表面温度,其等效噪声温差(NETD)为30 mK,图像分辨率为640×120,采集帧频为200 Hz。计算机控制闪光灯电源及同步触发器,在红外热像仪单帧图像采集结束后,电源驱动闪光灯对样品表面进行加热。红外热像仪采集样品表面温度信息,并将其传输回计算机进行数据处理分析。选取部分热成像序列如图9所示。

      图  9  部分红外图像序列

      Figure 9.  Selected frames from the IR image sequence

      图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  涂层厚度-实验峰值时刻对数图

      Figure 11.  Thickness vs experimental peak time in double logarithm scale

      同一样本在另外四个区域取样。根据图11的校准曲线计算涂层厚度,结果如表2所示,厚度误差在5%范围内。

      表 2  涂层厚度检测结果

      Table 2.  Measurement result of coating thickness

      Thickness by
      micrometer/
      μm
      Thickness by
      thermography/
      μm
      Difference
      110399.73.3%
      2136139.22.3%
      3158160.61.6%
      4197189.44.0%

      图12为试件厚度分布的图像。厚度从右侧较薄的一侧逐渐过渡为左侧较厚的一侧。由于红外热像仪帧频和闪光灯能量强度的限制,检测区域仅包括试样的右半部分,涂层厚度在80~220 µm之间。

      图  12  试件厚度检测结果分布图

      Figure 12.  Test result of thickness distribution of the specimen

    • 文中建立了半透明涂层半无限大脉冲热传导简化理论模型,推导出半透明涂层表面温度场表达式。建立了半透明涂层的脉冲红外加热双层物理模型,分析了半透明涂层的脉冲红外加热表面温度场。证明了对于中等吸收系数和较薄的涂层,涂层厚度与表面温度的峰值时刻在对数坐标上存在线性关系,这为脉冲热成像技术测量半透明涂层厚度提供了理论依据。

      利用脉冲热成像技术对制备的厚度连续变化的半透明涂层试件进行了测试。检测结果与模拟结果吻合较好,实测厚度与实际厚度相差在5%以内,在一定程度上提高了直接进行半透明涂层厚度检测的精度。与常规脉冲红外热成像对于半透明涂层厚度检测相比,不再需要在样品表面喷涂黑漆或黑色石墨胶体,有利于在工业中直接对半透明涂层厚度进行检测。

参考文献 (19)

目录

    /

    返回文章
    返回