第44卷第7期 VO1.44 NO.7 红外与激光工程 Infrared and Laser En ̄in 2015年7月 Jul_2015 基于相位的热障涂层厚度及其脱粘缺陷红外定量识别 陈林,杨立,范春利,石宏臣,赵小龙 (海军工程大学动力工程学院,湖北武汉430033) 摘 要:为了促进热障涂层红外无损检测的定量检测研究,建立了轴对称圆柱坐标下的热障涂层脉 冲相位检测模型,针对研究模型,采用有限体积法求解出脉冲热激励下的温度场,将温度进行FFT变 化得到相位分布,分析了不同因素对检测表面相位差分布的影响。在此基础上,采用LM算法研究了 轴对称圆柱坐标下对热障涂层厚度的大小和脱粘缺陷的位置进行定量化检测的方法,分析了不同因 素对检测结果的影响。研究结果表明:当不存在测温误差时,不同的初始假设、采样窗口时间下都能 得到很高的识别精度,其对定量识别的影响不大,当测温仅存在均匀误差时,涂层厚度和脱粘缺陷位 置识别精度都很高,均匀误差对识别无影响,识别结果的精度会随测温随机误差的增大而降低,但在 较大的随机误差下仍有较高的识别精度。 关键词:热障涂层; 无损检测; 脱粘缺陷; 反演识别; 脉冲相位 中图分类号:TK38 文献标志码:A 文章编号:1007—2276(2015)07—2050—07 Quantitative identiifcation of coating thickness and debonding defects of TBC by pulse phase technology Chen Lin,Yang Li,Fan Chunli,Shi Hongchen,Zhao Xiaolong (College of Power Engineering,Naval University of Engineering,Wuhan 430033,China) Abstract:In order to promote quantitative detection of thermal barrier coating by pulse phase,the pulse phase nondestructive testing model of thermal barrier coating in axisymmetric cylindrical coordinates was established.Temperature field was solved by using a finite volume method in the heat pulse excitation. And the phase of the model was accepted by FFT.The impacts of different factors on the phase were analysised.Based on this,quantitative detection of coating thickness and posiion of tthe bonding defect of hertmal barrier coating were researched in axisymmetric cylindrical coordinates by using the method of LM.The impacts of different factors on quantitative detection were analysised.The research results show hatt the recognition accuracy of results are high in different initil assumption and diafferent modulation time when there is no measurement error,the effect of initial assumption and modulation time on quantitative recognition is little,the recognition accuracy of results are high in different uniform temperature measurement error,uniform temperature error would not change the result of recognition,the accuracy of the identiicatfion results would decrease with the increase of the random temperature error, 收稿日期:2014—11—11: 修订日期:2014—12—17 基金项目:国家自然科学基金(50906099,51479203) 作者简介:陈林(1986一),男,博士生,主要从事定量红外缺陷检测方面的研究。Email:chenlinhonghu@163.com 导师简介:杨 ̄(1962一),男,教授,博士生导师,博士,主要从事传热传质热流体学应用及量测方面的研究。Email:lyang39@hotmail.tom 第7期 陈林等:基于相位的热障涂层厚度及其脱粘缺陷红外定量识别 2051 but there is still a high recognition accuracy under large random temperature error,the effectiveness of he quatntitative identification method was proved. Key words:thermal barrier coating;non-destructive testing;bonding defect;inversion identification; pulse phase 外无损检测模型。采用有限体积法对所建立的模型进 O引言 行数值模拟,得到研究模型表厩温度场,将得到的温 度进行FFT变换。从而得到相位,分析表面热流强度、 热障涂层(Thermal Barrier Coating,TBC)已广泛 应用于各种燃气轮机和柴油机上【”,其常见结构由陶 瓷层、粘结层和耐高温基体构成,其中陶瓷层是机械 或部件正常工作的保证.而脱粘缺陷是热障涂层常 见的缺陷之一,因此对涂层厚度及其脱粘缺陷的检 测具有重要的现实价值,对于此类问题,红外无损检 测相对于微波检测【 、超声检测【s】等方法具有不接 触、适应大面积检测等优点。参考文献[4】研究了涂层 缺陷的红外检测方法,参考文献[5】对涂层厚度的检 测方法进行了研究,参考文献【6—7】通过数值模拟的 方法对热障涂层红外热无损检测的建模和热障涂层 厚度不均的热无损检测进行了研究,参考文献【4—7】 的研究都只是基于表面温度,且研究还都处于定性 阶段。远不能满足工程应用的需要,而脉冲相位检测 作为一种新的红外检测方法,相比于表面温度法,具 有对表面发射率变化和加热不均不敏感、抗反射干 扰强等优点【。_9】,然而现有文献大都是在一维导热下 对脉冲相位检测进行研究的,而一维导热的简化处理往 往会带来较大的检测误差,所以参考文献[10—11]在确 定真实的缺陷时都要通过曲线拟合的方法对结果进 行修正,而这种拟合修正得到的结果属于一种半经 验的方法,有时误差很大,并且不同的检测材料需要 分别进行拟合。现有的脉冲相位检测法在缺陷定量 识别时需要准确确定峰值频率和“盲频”,但在脉冲 相位检测中。频率是一些离散的值,有时可能因频率 的细化程度不足而无法准确确定峰值频率和“盲频”, 从而难以准确定量识别缺陷,因而参考文献【l2一l3】 提出了对频率进行细化的不同方法,然而这种细化 后的频率值仍然是一些有限的离散值。因此,研究基 于相位的热障涂层定量检测对于提高检测精度及脉 冲相位检测的应用范围具有重要意义。 文中主要研究热障涂层厚度及其脱粘缺陷的定 量识别。针对所研究的问题建立了一个热障涂层红 涂层厚度、缺陷宽度以及调制窗口时间对相位差的影 响,在此基础上,将相位作为定量检测的依据,通过 LM反演算法实现对热障涂层厚度及其脱粘缺陷位置 的定量识别,分析了不同因素对检测结果的影响。 1红外检测正问题模型 如图l所示.文中研究问题为轴对称圆柱坐标下 的热障涂层检测。模拟试件的半径总长为Lr=20 ITlln, 陶瓷层厚度为 ,参数 决定缺陷的深度,粘结层 的厚度为Ld=0.1 min。耐高温金属合金基体的厚度 为Ls=2nlin,在粘结层中有圆弧状的脱粘缺陷(其间 充满空气),缺陷厚度与粘结层厚度相同,为Ld=0.1 nlin, 缺陷左侧距圆柱体中心距离为Lrl,缺陷右侧距圆柱 体边界距离为Lr2。参数Lrl和Lr2决定缺陷的位置 大小。模型的加热面与周围空气的换热系数为 ,为 简化分析,模型侧面Lr=20inln处为绝热,非加热面 考虑为绝热化处理,也视为绝热,假设物体是均匀 的.各向同性,则轴对称圆柱坐标下瞬态导热偏微分 方程为: 图1 TBC脱粘缺陷红外检测模型 Fig.1 Nondestructive testing model of TBC by infrared technology 10(rk )+专 ㈣ 2052 红外与激光工程 44卷 初始条件: Tl ̄o=20℃ (2) (P)<8 ( 7) 式中:s为一个给定的极小的正数。 边界条件: 一 (厶 )=otqqs+ot(Tr- ̄ (3) 为求得目标函数 (P)的最小值,对公式(6)进行 偏微分求导: 一 昙 l 。:0 (4) (5) VS(P)=2 f-- 的迭代公式: AP=pk“一 =【( ) + ][y- 0 (8) 将 通过泰勒级数展开并去掉高阶项,再将衰 一I :0 o’ 减因数 添加到计算结果中去,从而得到LM算法 式中: 为试件内各离散点的温度,℃;k为试件的 导热系数,W/(m・℃);P为密度,kg/m3;c为比热容, J/ ・℃(各材料的物性参数如表1所示);乃为环境 温度取为20 oC;%=0.9为试件表面对热量的吸收系 ( ) [y_ ( )】 (9) 式中: =f旦 1 ; =diag[( 】;k为迭代次 数;衰减因数 为调节系数。迭始时,初始假 数;a=10W/(m ̄・℃)为表面换热系数;在t=0 ̄0.1 s 时,热流强度q,=106W/m2,在t>O.1 s时, =0。如无 设结果与真实值相差较大,为提高迭代收敛速度。 取较大值,随着迭代计算的进行,其结果越来越接 近真实值, 将逐渐减小直至得到精确的辨识结 果。 特殊说明,文中参数值都如上所述。 表l不同材料的物性参数 Tab.1 Thermal properties of different materials 2.2迭代求解过程 (1)给定一个脱粘缺陷位置和陶瓷层厚度的初 始假设Po,将P0作为输入条件带入迭代计算。 (2)在识别参数向量 下求解方程(1)一(5)得到 试件表面温度,将表面温度进行FFT变换,得到各个 离散点的相位 ( )( =0,l,2,…)。 (3)求解目标函数 ( ),并判断其是否满足收 敛条件(公式(7)),如果满足收敛条件则停止迭代,将 其作为定量辨识结果,否则依据公式(9)对参数 进 行修正得到下次的识别参数P“,并返回步骤(2)。 2基于脉冲相位的反演识别法 2.1 Levenberg-Marquardt反演算法 3检测表面相位差变化规律及分析 3.1不同热流强度对相位差的影响 Levenberg-Marquardt反演优化算法通常是对目 标函数 )进行最小化求解fH-巧】: ," 保持缺陷左侧距圆柱体中心距离为Lrl=12 mill, 缺陷右侧距圆柱体边界距离为Lr2=6mm,陶瓷层的 厚度 =0.3 mm不变,热流强度分别取为10 、2× (P)=[y一 (P)] ・[1 (P)]= [ 一 】 f=l (6) 106、3xl06W/m2,采样频率为10Hz(下同),采样窗口 时间(采样的时间总长)为50 s(下同),计算得到缺陷 式中:P= ,P ,…,pn)为待识别的参数向量,对于文 中的热障涂层厚度及其脱粘缺陷位置的定量识别, P=(Lrl,Lr2,Ly); 为实际测温值FFT变换后得到 的离散相位值; 为理论计算温度值FFT变换后得 到的相位值;m为离散数值点。 LM算法的收敛条件为目标函数S(P)达到最小 化.即 表面(取脱粘缺陷中心对应的检测表面)与非缺陷表 面(取圆柱体中心)的相位差变化如图2所示。由图2 可知,随着热流强度的增加,缺陷表面和非缺陷对应 表面的相位差没有变化,可见热流强度不影响缺陷 表面和非缺陷对应表面间的相位差.较小的热流强 度就能准确识别缺陷的存在。 4 红外与激光工程 值的真实值,由方程(1)一(5)计算得到表面的真实温 度,将测量误差添加到真实值上,如公式(1O)所示, 以此值作为实际检测时试件表面测得的温度。对得 到的温度进行FFT变化,得到表面的相位(反演识别 时的相位值为角度值),以此作为定量识别的依据, 然后假设一组热障涂层厚度及其脱粘缺陷位置值的 初始值,通过迭代计算得到最终定量检测结果。 = .tm,+e+O.)or(10) 式中: 为在第f个离散点的测得温度; . 为计算 得到的真实温度值;e为均匀误差; 为偏差程度;‘£I 为一1.0一1.0内的随机数。 为了判断识别结果的准确性。定义如下形式的 相对误差: l l :l1 生二生Il ×100% (11) 式中:Er/为第f个参数的识别相对误差;P:为第f个 参数的识别结果;P 为第f个参数的真实值。 在不考虑测温误差时,选择不同的窗口时间为 l0、30、50 S,脱粘缺陷位置和陶瓷层厚度的真实值为 Lrl=12mill、Lr2=6nllll、Ly=0.3mill,脱粘缺陷位置和 陶瓷层厚度的初始假设为Lrl=8 mln、Lr2=10 m/n、 =0.5mln,其识别结果如表2所示。从表中可以看 出,当不存在测温误差时,不同的窗I=I时间下都能准 确地识别涂层厚度和脱粘缺陷位置,可见窗口时间的 选取对定量识别的影响不大,而在传统脉冲相位检测 中,在确定的采样频率下采样窗口大小决定着频率的 细化程度,采样窗口大小对检测结果影响较大。 表2不同窗口时间下的检测结果 Tab.2 Identification results in different modulation time 在不考虑测温误差时,选取不同初始假设值,采 样窗口时间为50 S.脱粘缺陷位置和陶瓷层厚度的 真实值为Lrl=12i/lln、Lr2=6toni、 =0.3mln,其识别 结果如表3所示。从表中可以看出,当不存在测温误 差时,不同的初始假设下识别精度都很高,即不存在 测温误差时。不同初始假设下都能准确地识别涂层 厚度和脱粘缺陷位置,可见初始假设对检测结果的 影响较小。 表3不同初始假设下的检测结果 Tab.3 Identification results in diferent initial guess 在不考虑测温误差时.脱粘缺陷位置和陶瓷层 厚度在不同值下的识别结果如表4所示。识别时采 样窗口时间为50 S,脱粘缺陷位置和陶瓷层厚度的 初始假设为Lrl=8mln、Lr2=10I1Mn、Ly=0.5HIln。从表 中可以看出,当不存在测温误差时.脱粘缺陷位置和 陶瓷层厚度在不同值下的识别结果精度都很高.可 见脱粘缺陷位置和陶瓷层厚度在不同值下都能准确 识别。 表4不同缺陷位置下的检测结果 Tab.4 Identiifcation results in different position 实际检测中会存在测温误差,依据公式(10)的 误差模拟方法研究不同类型测温误差对检测情况。 采样窗口时间为50 S,脱粘缺陷位置和陶瓷层厚度的 初始假设分别为Lrl=8mln、Lr2=10mm、Ly=0.5mill。 当只有存在均匀误差而无随机误差时.识别结果如 表5所示。从表中可以看出,当仅存在均匀误差时, 涂层厚度和脱粘缺陷位置识别精度都很高.不同的 均匀误差不改变识别结果。可见均匀误差对基于相 位的检测无影响。 第7期 陈林等:基于相位的热障涂层厚度及其脱粘缺陷红外定量识剐 2055 表5不同均匀误差下的检测结果 Tab.5 IdentiifcatiOn results in distinct uniform temperature measurement error 表6为只考虑随机误差而无均匀误差时的检测 情况,采样窗口时间为50 S,脱粘缺陷位置和陶瓷层 厚度的初始假设分别为Lrl=8min、Lr2=10mill、Ly= 0.5IIllTI。从表中可以看出,当不存在测温误差时,涂 层厚度和脱粘缺陷位置识别精度都很高。随着测温 随机误差的增大。识别结果的精度降低,即测温误差 越小,越能准确识别涂层厚度和脱粘缺陷位置。 表6不同随机误差下的检测结果 Tab.6 Identification results in distinct random temperature measurement error 4结论 文中利用有限体积法建立了轴对称圆柱坐标下 的热障涂层脉冲相位检测模型,分析不同因素对表 面相位差分布的影响.在此基础上采用LM算法对 热障涂层厚度的大小和脱粘缺陷的位置进行了基于 相位的定量化检测,分析了不同因素对检测结果的 影响,得到的结论主要有: (1)相位差随着热流强度等增大而保持不变,热 流强度的改变对相位差的分布无影响.相位差随着 陶瓷层厚度的增加而减小。随着脱粘缺陷宽度的增 加而增大。 (2)当不存在测温误差时,在不同的初始假设、 采样窗口时间下涂层厚度和脱粘缺陷位置参数都有 很高的识别精度,初始假设和采样窗口时间对定量 识别的影响不大;当不存在测温误差时,脱粘缺陷位 置和陶瓷层厚度在不同值下都能准确识别。 (3)当仅存在均匀误差时,涂层厚度和脱粘缺陷 位置识别精度都很高,不同的均匀误差不改变识别 结果,可见均匀误差对基于相位的检测无影响。 (4)当测温存在随机误差时,涂层厚度和脱粘缺 陷位置参数的识别结果精度会随着测温随机误差的 增大而降低,但在较大的测温误差下仍有较高的识 别精度.表明了定量识别方法的有效性。 参考文献: 【1】Li Yong,Chen Zhenmao,Mao Ying,et a1.Quantitative evaluation of thermal barirer coating based on eddy current technique[J].NDT&E International,2012,50:29-35. [2】 Sayar M,Seo D,Ogawa K.Non—destructive microwave detection of layer thickness in degraded thermal barrier coatingsusingK—and W-band frequency range【J1.NDT&E International,2009,42:398—403. 【3】 Roge B,Fahr A,Giguere J S R,et a1.Nondestructive measurement of porosity in thermal barrier coatings【J】. Journal of Thermal Spray Technology,2003,12(4):530- 535. [4】Liu Huanan.Michim Sakamoto,Kazushi Kishi,et a1. Detection of defects in thermal barirer coatings by htermography analyses【J】.Materials Transactions,2003,44 (9):1845-1850. 【5】 Franke B,Sohn Y H,Chen X.Monitoring damage evolution in thermal barrier coatings wiht thermal wave imaging【J】. Surface and Coatings Technology,2005,200(5):1292—1297. 【6】 Guo Xingwang,Ding Mengmeng.Modeling and finite element analysis of thermal barrier coatings in IR NDT【J】.Journal of Beijing University ofAeronautics and Astronautics,2009, 35(2):174-178.(in Chinese) 【7】Guo Xingwang,Ding Mengmeng.Simulation of thermal NDT of thickness and its unevenness of thermal barrier coatings 【J】.Acta Aeronautica JEt Astronautica Sinica,2010,31(1): 198-203.(in Chinese) 【8]Vitali L,Fustinon D,Gramazio P,et a1.Estimation of trial parameters for pulse phase thermography with low power heat sources[J].,ournal of Physics:Conference Sedes, 2014,501:012008. 【9】 Stotter B,Gresslehner K H,Mayr G,et a1.Estimation of 2056 红外与激光工程 material parameters from pulse phase thermography ̄ata[C]// AIP Conference Proceedings,2014,1581:1126—1133. 【1O】 Li Yanhong,Zhao Yuejin,Feng Lichun,et a1.Measurement of defect depth by infrared thermal wave nondestructive evaluation based on pulsed phase[J】.Opacs and Precision Engineering,2008,16(1):56—59.(in Chinese) [11】 Yu Jiajie,Wu Naiming,Zeng Zhi,et a1.Fl冲depth measurement based on pulsed phase thermography【J】. Infrared and Laser Engineering,2012,41(7):1894—1896. (in Chinese) 【l2】 Ma Shuohan,Ma Qishuang.Pulse phase thermography chirp Z transform[J]..tOurnal of Belling University ofAeronautics ㈣ and Astronautics,2011,37(11):1436-1438.(in Chinese) Ma Shuohan,Ma Qishuang.Algorithm of infrared pulse phase thermography mulitple modulation zoom spectrum inspection[J].Infrared and Laser Engineering,2012,4l(8): 2223-2226.(in Chinese) Liu F B.A hybrid method for the inverse heat transfer of estimating fluid thermal conductivity and heat capacity【J】. InternationalJournalofThermalSciences,2011,50:718-724. Fan Chunli,Sun Fengmi,Yang Li.A quantitative identiifcation technique for a two—dimensional subsurface defect based on surface temperature measurement【J】.Heat Transfer-Asian Research,2009,38(4):223-233.