《热障涂层厚度及其涂敷下材料缺陷的红外定量识别.DOC》由会员分享，可在线阅读，更多相关《热障涂层厚度及其涂敷下材料缺陷的红外定量识别.DOC（11页珍藏版）》请在三一文库上搜索。
1、热障涂层厚度及其涂敷下材料内部缺陷的红外定量识别 收稿日期： 基金项目：国家自然科学基金资助项目(50906099、51479203)作者简介：陈林（1986），男，湖北洪湖人，博士研究生，Email：chenlinhonghu杨立（1962），男，教授，博士生导师，博士， Email：lyang39陈林，杨立，范春利，石宏臣，赵小龙（海军工程大学 动力工程学院 湖北 武汉 430033）摘 要：针对热障涂层结构材料红外定量检测存在的不足，提出了基于脉冲相位的LM（Levenberg-Marquardt）识别算法。针对研究问题，建立了轴对称圆柱坐标下热障涂层结构材料的瞬态导热模型，利用有限体积
2、法计算出检测表面的温度场，经 FFT变换得到检测表面的相位分布，分析了待检测参数对检测表面相位差的影响。以相位为识别条件，采用LM反演算法对热障涂层厚度及其涂敷下材料内部缺陷位置大小进行定量识别，并采用数值实验的方法对基于相位的LM识别方法进行了验证，数值实验结果验证了识别方法的有效性。关键词：热障涂层；红外检测；定量识别；脉冲相位中图分类号：TK38 文献标志码：A 文章编号：Quantitative identification of coating thickness and inner defects of the material by infrared testing techno
3、logyChen Lin, Yang Li, Fan Chunli, Shi Hongchen Zhao Xiaolong (College of Power Engineering, Naval University of Engineering, Wuhan 430033, China)Abstract：Aiming at the problems of infrared quantitative detection of thermal barrier coating structures, the LM (Levenberg-Marquardt) recognition algorit
4、hm based on pulse phase was put forward. The transient thermalmodel of thermal barrier coating structures by an axisymmetric cylindrical coordinates was established in view of the research problem. The surface temperature field of the research model was calculated by the finite volume method. And th
5、en the phase of the inspection surface was obtained through FFT. The impacts of parameters to be detected on the phase of the inspection surface were analyzed. Based on the phase for the input parameter, quantitative identification method of coating thickness and inner defects of the metal was resea
6、rched by using LM recognition algorithm. The quantitative identification method based on phase detection was verified by using numerical experiment. The effectiveness of the quantitative identification method was evidenced on the basis of the numerical experiment results.Key words：thermal barrier co
7、ating；infrared testing； Quantitative identification；pulse phase热障涂层（TBC，Thermal barrier coating）由于其良好的隔热和防腐蚀性能在航空航天等领域得到广泛的应用。常见的热障涂层结构材料由陶瓷层、粘结层和合金基体构成，为了保证热障涂层结构材料的耐高温性和抗腐蚀性，热障涂层结构材料都有其合理的涂层厚度，而材料的完好是机械或部件正常工作的基本要求，因此对热障涂层厚度及其涂敷下材料内部缺陷的检测具有重要的意义。与微波检测1、超声检测2等传统检测方法相比，红外无损检测法对热障涂层结构材料的检测具有检测面积大、速度快
8、等优点。红外热像检测法主要有脉冲检测法3、锁相检测法4、脉冲相位检测法5等，其中脉冲相位检测法因兼具有脉冲检测法和锁相检测法的优点而在红外检测中得到广泛的应用。文献6采用脉冲红外热像法对涂层厚度检测进行了研究，文献7将脉冲红外热像法应用于涂层厚度及其失效过程的监测，文献8-9 对涂层结构板材内部缺陷的脉冲相位红外热成像法进行了研究。然而现有文献对相位检测的定量化研究都是通过热传导机理来进行，由于受复杂边界条件的影响，在实际传热模型的求解时，多数情况下被检测参数对检测表面相位分布的影响都不能得到精确的解析解，因而为了得到被检测参数的定量化结果，大都是在一维简化模型下进行研究的，这样往往带来检测结
9、果与实际情况存在较大的偏差。而依据导热反问题理论来求解被检测参数的方法能很好地满足定量化检测的要求，因此将反问题算法与红外检测技术相结合对于红外检测技术从定性到定量的转变具有重要的意义。目前导热反问题优化算法主要有共轭梯度法10、数值搜索算法11、遗传算法12等方法，其中，LM（Levenberg-Marquardt）因求解方法简单、收敛速度快等优点在反问题求解方面得到较广的应用13-14。本文主要通过将反问题算法与脉冲相位检测技术相结合来研究热障涂层厚度及其涂敷下材料内部缺陷的定量识别。 1 红外检测正问题模型在轴对称柱坐标系下，三维导热问题可以简化为二维导热问题来研究，这样将大大节省计算时
10、间，因而本文研究问题为轴对称圆柱坐标下的热障涂层结构材料的红外无损检测，如图1所示。研究模型的半径，Ly3为陶瓷层厚度，为粘结层的厚度，在一般的热障涂层结构材料中粘结层的典型厚度为15，基体的厚度，在合金基体中有圆弧状的内部缺陷，缺陷内充满空气，Lr1为缺陷距圆柱体中心的距离，Lr2为缺陷距圆柱体外边界的距离，Ly1为缺陷距非加热面的距离，Ly2为缺陷距加热面的距离，参数Lr1、Lr2、Ly1和Ly2决定缺陷的位置和大小。为简化分析，模型外侧面处视为绝热，非加热面考虑为绝热化处理，也视为绝热，假设材料是均质的，各向同性，则可建立如下的传热方程：图1热障涂层结构材料红外无损检测模型Fig.1 T
11、he nondestructive testing model of thermal barrier coating structures by infrared technology (1)初始条件： (2)边界条件： (3) (4) (5)其中：T为温度（）；为导热系数（）；为密度（）；c为比热容（）（各材料的物性参数如表1所示）； 为环境温度； 为检测表面的热量吸收系数；为检测表面与环境的换热系数；在，热流强度，在时，。表1 不同材料的物性参数Tab.1 Thermal properties of different materials材料密度比热容热导率陶瓷层60375001.5粘结层
12、73204504.3合金基体83004409.5空气1.210050.02572反演识别方法 热障涂层厚度及其涂敷下材料内部缺陷的反演识别可描述为：依据检测表面的相位来识别热障涂层厚度及材料内部缺陷的位置大小。在反问题求解时，公式（1）（5）中除了参数Lr1，Lr2，Ly1，Ly2，Ly3为未知量外，其他参数均为已知。检测表面的温度通过红外热像仪测得，检测表面的相位通过对温度进行FFT变换得到。2.1 Levenberg-Marquardt反演识别算法反演优化算法通常是对目标函数S(P) 进行最小化求解： （6）其中，是待识别的参数向量，对于本文热障涂层厚度及材料内部缺陷位置大小的定量识别，Y
13、i为实际检测时得到的相位值，为在参数向量下通过公式（1）（5）计算得到的相位值； m为检测表面的离散点数。由于在“峰值”频率（检测表面相位差最大时对应的频率）下缺陷的可检测性最好，因而在反演识别时，只选取“峰值”频率下检测表面的离散相位值作为识别输入参数。 从偏差原理出发，当目标函数满足如下收敛准则时，可以停止迭代求解并认为得到最终的识别结果： （7）其中，为一个给定的极小的正数。 为了使目标函数达到最小化，将目标函数对每个识别参数进行偏微分求导，并令每一个偏导数等于零，即： （8） 将通过泰勒级数展开并略去高阶项，再将一个调节系数添加到计算结果中去，从而得到LM算法的迭代公式： （9）其中：
14、， ， ， 为迭代次数。2.2 迭代求解过程（1） 给定一个热障涂层厚度及材料内部缺陷位置大小的初始假设，将作为开始迭代计算的输入条件。（2） 在参数向量下求解方程得到检测表面各测温点的温度，对检测表面测温点温度进行FFT变换，得到测量点处相位（k=0,1,2,）。（3） 在参数向量下求解目标函数S()，并依据收敛标准（7）判断其是否满足收敛条件，如果满足收敛条件则停止迭代，将其作为最终定量识别结果，否则依据公式（9）对参数进行修正得到下次的识别参数，并返回步骤（2）。3 检测表面的相位相位差分布规律为了获得检测表面的相位分布，利用有限体积法对建立的模型进行脉冲激励下温度场的模拟计算，对模拟计
15、算得到的温度进行FFT变换，得到缺陷对应检测表面和非缺陷对应检测表面温度变化的相位信息，采样窗口时间为10s，采样频率为10Hz。材料内部缺陷位置大小参数Lr1=12mm、Lr2=6mm、Ly1=15mm、Ly2=0.6mm，陶瓷层厚度参数Ly3=0.3mm，下文的分析都是在变化其中某一参数值而保持其他参数值不变下进行的。3.1参数Ly1对检测表面相位差的影响材料内部缺陷位置大小参数 Ly1分别取为1.3mm、1.5mm、1.6mm，缺陷（缺陷中心）对应检测表面与非缺陷（圆柱中心）对应检测表面的相位差分布如图2所示，由图可知，随着参数Ly1值的增大，检测表面的相位差最大值（指绝对值）减小，但是
16、检测表面相位差的改变量较小，可见参数Ly1对检测表面相位分布的影响较小。图2 参数Ly1对检测表面相位差的影响Fig.2 Impact of parameter Ly1 on the phase difference3.2参数Ly2对检测表面相位差的影响材料内部缺陷位置大小参数 Ly2分别取为0.5mm、0.6mm、0.7mm，缺陷（缺陷中心）对应检测表面与非缺陷（圆柱中心）对应检测表面的相位差分布如图3所示，由图可知，与参数Ly1相同，检测表面的相位差最大值（指绝对值）随着参数Ly2值的增大而减小，并且检测表面相位差的改变量相对较大，可见参数Ly2对检测表面相位分布的影响较大。图3 参数Ly
17、2对检测表面相位差的影响Fig.3 Impact of parameter Ly2 on the phase difference 3.3参数Ly3对检测表面相位差的影响陶瓷层厚度参数Ly3分别取为0.2mm、0.3mm、0.4mm，缺陷（缺陷中心）对应检测表面与非缺陷（圆柱中心）对应检测表面的相位差分布如图4所示，由图可知，随着参数Ly3值的增大，检测表面的相位差最大值（指绝对值）增大，检测表面相位差的改变量相对较大，其对检测表面相位分布的影响较大。图4 参数Ly3对检测表面相位差的影响Fig.4 Impact of parameter Ly3 on the phase difference
18、 3.4参数Lr1对检测表面相位差的影响参数Lr1分别取为12.0mm、12.2mm、12.4mm，缺陷（与上一小节相同）对应检测表面与非缺陷（圆柱中心）对应检测表面的相位差分布如图5所示，由图可知，随着参数Lr1值的增大，检测表面的相位差最大值（指绝对值）减小。 图5 参数Lr1对检测表面相位差的影响Fig.5 Impact of parameter Lr1 on the phase difference e3.5参数Lr2对检测表面相位差的影响参数Lr2分别取为 6.0mm、6.2mm、6.4mm，缺陷（与上一小节相同）对应检测表面与非缺陷（圆柱中心）对应检测表面的相位差分布如图6所示，由
19、图可知，与参数Lr1相同，检测表面的相位差最大值（指绝对值）随着参数Lr2值的增大而减小。图6 参数Lr2对检测表面相位差的影响Fig.6 Impact of parameter Lr2 on the phase difference e4 定量检测算例及分析在进行热障涂层厚度和材料内部缺陷位置大小的定量识别时，首先假定一组热障涂层厚度和材料内部缺陷位置大小的真实值，由导热方程(1)（5）计算出检测表面的真实温度，如式（10）所示，将测温误差添加到真实值上，以此值作为实际测得的温度，对温度进行FFT变换，得到检测表面的相位（反演识别时的相位值为角度值），以此作为定量识别的依据，然后再假设一组热
20、障涂层厚度和材料内部缺陷位置大小的初始值，通过迭代计算得到最终定量检测结果。 （10）式中：为在第i个离散点的测得温度，为第i个离散点的真实温度值，e为均匀误差，为偏差程度，为-1.0到1.0内的随机数。为了判断识别结果的准确性，定义如下形式的相对误差： （11）为第i个参数的识别相对误差，为第i个参数的识别结果，为第i个参数的真实值。在实际红外无损检测时，检测表面可能存在加热不均，为了研究检测表面加热不均对识别结果的影响，将一个随机偏差添加到加载的热流上，以此来模拟实际红外无损检测时检测表面加载的热流。 （12） 为第i个离散点的实际热流强度，为加载的热流强度，为相对偏差程度，为-1.0到1
21、.0内的随机数。当测温不存在误差时，在不同初始假设值下，热障涂层厚度和材料内部缺陷位置大小的识别结果如表2所示，从表中可以看出，当不存在测温误差时，不同初始假设下都能准确地识别热障涂层厚度和材料内部缺陷位置大小，初始假设对识别结果的影响不大。表2 无测温误差时初始假设对识别结果的影响Tab.2 Impact of initial guess on identification results in the presence of no measurement error 初始假设/mm 识别结果/mm 相对误差(%): Er15 , 3, 1.8, 0.3, 0.1 12.00, 6.00,
22、1.50, 0.60, 0.30 ,8, 10, 1.0, 0.8, 0.5 12.00, 6.00, 1.50, 0.60, 0.30 ,5, 12, 1.2, 1.0, 0.6 12.00, 6.00, 1.50, 0.60, 0.30 ,在实际红外无损检测时，测温会存在测量误差，依据公式（10）研究不同形式的测温误差对识别结果精度的影响, 热障涂层厚度和材料内部缺陷位置大小的初始假设分别为Lr1=8mm、Lr2=10mm、Ly1=1.0mm、Ly2=0.8mm 、Ly3=0.5mm。当只存在均匀误差而无随机误差时，识别结果如表3所示，从表中可以看出，当仅存在均匀误差时，在不同的测温误差下
23、，热障涂层厚度和材料内部缺陷位置大小都能准确识别，均匀误差不改变识别结果，可见均匀误差对识别结果无影响。表3 均匀误差对识别结果的影响Tab.3 Impact of uniform temperature measurement error on identification results均匀误差/0C 识别结果/mm 相对误差(%): Er0 12.00, 6.00, 1.50, 0.60, 0.30 , 0.2 12.00, 6.00, 1.50, 0.60, 0.30 , 0.5 12.00, 6.00, 1.50, 0.60, 0.30 ,1 12.00, 6.00, 1.50, 0
24、.60, 0.30 ,当测温存在随机误差而无均匀误差时，热障涂层厚度和材料内部缺陷位置大小的识别结果如表4所示，从表中可以看出，测温误差越小，热障涂层厚度和材料内部缺陷位置大小的识别结果精度越高，随测温误差的增大，识别结果的精度将降低，并且测温误差对不同参数的影响不同，参数Ly1的识别结果受测温误差的影响最大，这种差异性是由于不同参数的扰动对检测表面相位分布的影响程度不同造成的。 表4 随机误差对识别结果的影响Tab.4 Impact of random temperature measurement error on identification results 识别结果/mm 相对误差(%
25、): Er 0 12.00, 6.00, 1.50, 0.60, 0.30 , 0.2 12.035, 6.031, 1.532, 0.592, 0.298 0.29, 0.52, 2.13, 1.33, 0.670.5 12.410, 5,737, 1.691, 0.578, 0.303 3.42, 4.38, 12.73, 3.67, 1.01 11.297, 6.386, 1.274, 0.569, 0.305 5.86, 6.43, 15.07, 5.17, 1.67 表5为相同的测温误差（）下，不同初始假设对识别结果的影响，从表中可以看出在相同的测温误差下，不同初始假设下的识别结果精
26、度相差不大，可见初始假设对识别结果影响较小。表5当存在测温误差时不同初始假设对识别结果的影响Tab.5 Impact of initial guess on identification results in the presence of measurement error初始假设/mm 识别结果/mm 相对误差(%): Er15 , 3, 1.8, 0.3, 0.1 11.295, 6.561, 1.307, 0.638, 0.306 5.88, 7.32, 12.87, 6.33, 2.08, 10, 1.0, 0.8, 0.5 11.297, 6.386, 1.274, 0.569,
27、0.305 5.86, 6.43, 15.07, 5.17, 1.675, 12, 1.2, 1.0, 0.6 11.497, 6.456, 1.340, 0.628, 0.303 4.19, 7.60, 10.67, 4.67, 1.0在实际红外无损检测时，加载的热流可能存在加热不均，依据公式（12）的模拟方法研究在相同的测温误差（）下不同程度的加热不均对识别结果的影响，识别结果如表6所示。热障涂层厚度和材料内部缺陷位置大小的初始假设分别为Lr1=8mm、Lr2=10mm、Ly1=1.0mm、Ly2=0.8mm 、Ly3=0.5mm。从表中可以看出，在相同的测温误差下，随表面加热不均的增大，
28、热障涂层厚度和材料内部缺陷位置大小的识别结果精度未出现较明显的降低，可见识别结果精度对加热不均不敏感。表6 检测表面加热不均对识别结果的影响Tab.6 Impact of heating uneven on identification results 相对偏差(%) 识别结果/mm 相对误差(%): Er0 11.297, 6.386, 1.274, 0.569, 0.305 5.86, 6.43, 15.07, 5.17, 1.675 12.464, 5.506, 1.747, 0.571, 0.297 3.87, 8.23, 16.47, 4.83, 1.010 12.803, 5.55
29、4, 1.659, 0.565, 0.305 6.69, 7.43, 10.0, 5.83, 1.675 结论本文建立了热障涂层厚度及其涂敷下材料内部缺陷的红外检测模型，分析了待检测参数对检测表面相位差的影响，采用LM算法实现了对热障涂层厚度及其涂敷下材料内部缺陷的定量化识别，分析了初始假设、不同测温误差和加热不均对识别结果的影响，得到的结论主要有：（1）检测表面相位差最大值（指绝对值）随材料内部缺陷位置大小参数Lr1、Lr2、Ly1和Ly2的增大而减小，而随着热障涂层厚度参数Ly3的增大，检测表面相位差最大值（指绝对值）将增大，但不同参数的扰动对检测表面相位差分布的影响不同，参数Ly1的扰动
30、对相位差分布的影响最弱。（2）当不存在测温误差时，热障涂层厚度和材料内部缺陷位置大小在不同的初始假设下都有很高的识别精度；当存在测温误差时，在相同的测温误差下，热障涂层厚度和材料内部缺陷位置大小在不同初始假设下的识别结果精度相差不大，可见初始假设对定量识别的影响不大。（3）当测温仅存在均匀误差时，热障涂层厚度和材料内部缺陷位置大小的识别结果精度都很高，识别结果不随不同测温误差的增大而改变，可见均匀误差对识别结果无影响。（4）当测温仅存在随机误差时，测温误差较小时热障涂层厚度和材料内部缺陷位置大小的识别精度很高，随测温误差的增大，识别结果的精度将降低，但在较大的测温下仍有较高的识别精度，表明了识
31、别方法的有效性，并且测温误差对不同参数识别精度的影响不同，参数Ly1的识别精度最易受测温误差的影响。（5）在相同的测温误差下，随表面加热不均的增大，识别结果的精度未出现较明显的降低，可见识别结果精度对加热不均不敏感。参考文献1 Sayar M, Seo D, Ogawa K. Non-destructive microwave detection of layer thickness in degraded thermal barrier coatingsusingK-and W-band frequency rangeJ.NDT&E International, 2009,42:39
32、8-403.2 Roge B, Fahr A, Giguere J S R, et al. Nondestructive Measurement of Porosity in Thermal Barrier CoatingsJ.Journal of Thermal Spray Technology, 2003, 12(4): 530-535.3 肖劲松, 严天鹏. 风力机叶片的红外热成像无损检测的数值研究J.北京工业大学学报, 2006, 32(1): 48-52.Xiao Jingsong, Yan Tianpeng. Numerical Analysis of NDT of Wind Tu
33、rbine Blades Using Infrared ImagingJ. Journal of Beijing University Technology ,2006, 32(1): 48-52. (in Chinese)4 Bai W, Wong B S. Evaluation of defects in composite plate under connective environments using lock-in thermographyJ.Measurement Science and Technology, 2001(12): 142-150.5 马说邯, 马齐爽. 红外脉冲
34、相位复调制细化检测算法J.红外与激光工程, 2012, 41(8):2223-2225.Ma Shuohan, Ma Qishuang. Algorithm of infrared pulse phase thermography multiple modulation zoom spectrum inspectionJ. Infrared and Laser Engineering, 2012, 41(8):2223-2225. (in Chinese)6 Shepard S M, Hou Y L, James R, et al. Thermographic measurement of t
35、hermal barrier coating thicknessC/ Thermosense XXVII, Proceeding of the SPIE, 2005:407-410.7 Franke B, Sohn Y H, Chen X. Monitoring damage evolution in thermal barrier coatings with thermal wave imagingJ. Surface and Coatings Technology, 2005, 200(5):1292-1297.8 Tang Q J, Liu J Y, Wang Y, et al. Exp
36、erimental study of inspection on SiC coated high-temperature alloy plates with defects using pulsed thermographic techniqueJ. Infrared Physics & Technology, 2013,57(3):21-27.9 Tang Q J, Liu J Y, Wang Y, et al. Inspection on SiC coated carboncarbon composite with subsurface defects using pulsed t
37、hermographyJ. Infrared Physics & Technology, 2013,60:183-189.10 Chen W L, Yang Y C. Inverse estimation for unknown fouling-layer profiles with arbitrary geometries on the inner wall of a forced-convection ductJ.International Journal of Thermal Sciences,2010,49:86-98.11 Hsieh C K, Su K U. A metho
38、dology of predicting cavity geometry based on the scanned surface temperature data-prescribed heat flux at the cavity sideJ. ASME Journal of Heat Transfer, 1981,103(1):42-46.12 梅林, 吴立德, 王裕文. 一种基于遗传算法的红外无损检测缺陷的新方法J. 光学学报, 2002, 22(12):1452-1456.Mei Lin, Wu Lide, Wang Yuwen. A novel defect evaluation
39、method in infrared nondestructive testing based on genetic algorithmJ. Acta Optica Sinica, 2002, 22(12): 1452-1456.(in Chinese)13 Liu F G. A hybrid method for the inverse heat transfer of estimating fluid thermal conductivity and heat capacityJ.International Journal of Thermal Sciences, 2011,50:718-
40、724.14 Fan C L, Sun F R, Yang Li. A quantitative identification technique for a two-dimensional subsurface defect based on surface temperature measurementJ.Heat Transfer -Asian Research, 2009,38(4):223-233.15 丁彰雄. 热障涂层的研究动态及应用J.中国表面工程, 1999, 2:31-37.Ding Zhangxiong. Research status of thermal barrier coatings and applicationJ.China Surface Engineering , 1999, 2:31-37.(in Chinese)11 / 11文档可自由编辑打印