Nonlinear Data Fitting for Reflective Continuous Heat Excited Thermography Testing
-
摘要: 为了在反射式持续热激励红外热像技术中对缺陷深度进行定量检测，基于实际情况建立了持续热激励下被检物内的热传导模型，并求解了被检物热激励面上的温差-时间关系。通过对求解出的温差-时间关系的分析，发现实际检测时取热激励面温差-时间数据依照所求解出的关系进行非线性曲线拟合，可测量缺陷的深度。为了检验方法的可行性，利用风电叶片玻璃钢制作了平底洞试件，并对平底洞试件进行了检测实验。实验表明，本文方法对缺陷深度具有较高的测量精度。Abstract: To quantitatively detect defects using reflective continuous-heat-excited thermography, a heat conduction model of an object under continuous heat excitation was established, and the temperature increment-time relationship on the thermal excitation surface of the object was derived. Based on an analysis of the temperature increment-time relationship on the thermal excitation surface, the depth of the defects could be measured by nonlinear fitting of the temperature increment-time data. To test the feasibility of this method, a GFRP flat-bottomed hole specimen was fabricated and analyzed using reflective continuous-heat-excited thermography. The results show that this method is highly accurate in measuring the depth of defects.
-
Keywords:
- continuous-heat excited /
- nonlinear data fitting /
- reflective /
- heat conduction
-
0. 引言
随着红外制导技术的不断发展，装有红外传感器的导弹是飞机的最主要威胁，据统计，1973-2001年有49%的飞机损失来自红外导弹，其中1984-2001年90%的飞机损失来自红外导弹[1-2]。为保护飞机免受红外制导导弹的攻击，红外干扰弹应运而生，现今发展的红外干扰弹包括多点源红外干扰弹和面源红外干扰弹。多点源红外干扰弹依靠烟火药柱产生一个与飞机红外辐射特性类似，但能量大于飞机辐射能量2~10倍的热源，以质心干扰的方式达到欺骗来袭红外制导导弹的目的[3]；面源红外干扰弹依靠在空间形成大面积红外辐射云团，使飞机“淹没”在红外辐射云团中，以掩盖飞机自身的红外辐射特征，导致红外导引头无法分辨飞机红外辐射源和面源红外干扰弹辐射源[4-5]，飞机采用机动规避的方式逃逸出导引头视场，从而免遭红外制导导弹的攻击；多点源红外干扰弹在空中形成强红外辐射源诱骗红外导引头，面红外辐射源在空中形成大面积红外辐射源掩盖飞机自身的红外辐射源，两者结合，可有效迷惑红外制导导引头，提高载机的自卫对抗能力。
飞机的红外辐射特性是红外制导导弹识别、捕获、跟踪、摧毁飞机的信号来源，飞机的红外辐射信号主要由机身蒙皮、尾喷管和尾焰组成[4]，其中尾焰和尾喷管是飞机红外辐射的主要来源，红外制导导弹从不同方向攻击飞机时，由于方向不同，导致导引头探测目标的辐射强度和辐射面积差异很大，飞机投放干扰弹的时机和种类也不同。因此，探测飞机在不同方向的辐射特征对干扰弹弹种选择具有重要意义。
本文采用光谱辐射计、中波红外热像仪、长波红外热像仪对飞机在不同方向上的红外辐射性能进行测试，研究飞机发动机工作时在不同方向上的红外辐射性能，指导飞机应对不同方向来袭的红外制导导弹时红外干扰弹投放使用策略。
1. 测试
1.1 测试场景
选择飞机0°、正侧向90°、尾后180°三个测试角度，选择合适测试点，飞机发动机处于正常工作状态，测试场景示意图如图 1所示，测试距离300 m，飞机包络尺寸为不大于50 m×50 m×10 m，飞机发动机直径不大于1 m，因此，在测试场景下飞机可当作点目标。
1.2 测试方法
1.2.1 光谱辐射计标定
试验前应对光谱辐射计进行标定，在同一测试条件下，对黑体和被测试目标进行测试，标定对应温度下的黑体温度。
$$ I（\lambda ）=\frac{({V}_{\rm{s}}-{V}_{\rm{b}}){R}^{2}}{({V}_{\rm{BB}}-{V}_{\rm{Bb}}){L}^{2}}\times M({T}_{\rm{BB}})\times \frac{\rm{π}{D}^{4}}{4}$$ (1) 式中：Vs、Vb分别为飞机的信号、飞机附近的背景信号；VBB、VBb分别为黑体的信号、背景信号；R为测试距离；L为黑体标定距离；TBB为黑体温度；M(TBB)为黑体在TBB的光谱辐射出射度；D为光阑口径。
1.2.2 测试方法
选定不同方位（迎头、正侧向、尾后），架好仪器设备，发动机处于正常工作状态下，测试飞机的中/长波红外辐射强度、辐射面积。
1.3 数据处理
1.3.1 光谱辐射计测试数据处理
辐射计测试时，首先采用标准黑体源对光谱辐射计进行标定，得到光谱辐射计的光谱响应曲线；随后，采用标定后的光谱辐射计对待测对象进行测试，经过光谱强度积分和大气修正得到测试飞机的辐射强度曲线[6]。
按式(2)计算光谱响应系数：
$$K\left( \lambda \right) = \frac{{\Delta {V_{\rm{c}}}\left( \lambda \right) \cdot {{\left( {{L_{\rm{c}}}} \right)}^2}}}{{{A_{{\rm{bb}}}} \cdot P\left( {\lambda ,{T_{\rm{c}}}} \right)}}$$ (2) 式中：K(λ)为标定系数，V/(W/(μm·m2))；ΔVc(λ)为标定电压，V；Lc为标定距离，m；Abb为黑体腔口面积，m2；P(λ, Tc)为黑体的光谱辐亮度，W/(μm·sr·m2)；Tc为黑体温度，K。
按式(3)计算红外辐射强度：
$$I = \frac{1}{{{\tau _{\rm{a}}}}}\int_{{\lambda _{\min }}}^{{\lambda _{\max }}} {\frac{{\Delta {V_{\rm{m}}}\left( \lambda \right) \cdot {{\left( {{L_{{\rm{ob}}}}} \right)}^2}}}{{K\left( \lambda \right)}}} {\rm{d}}\lambda $$ (3) 式中：I为辐射强度，W/sr；ΔVm(λ)为测试电压，V；Lob为测试距离，m；λmin、λmax分别为波段的波长下限和上限，μm；τa为大气透过率（采用Lowtran大气模型根据实时气象条件计算）。
1.3.2 中长波红外热像仪测试数据处理
热像仪测试时，首先采用标准黑体源对热像仪进行标定，得到热像仪的辐射响应定标文件（定期进行）；采用热像仪对待测对象进行测试，得到飞机在不同方向上的温度分布图像，根据测试距离、视场角分辨率及图像像元数计算特征辐射面积S：
$$ \mathit{S}{\rm{ = }}{\mathit{n}_{\rm{T}}}{\rm{\cdot}}\mathit{\gamma }{\rm{\cdot}}\mathit{d} $$ (4) 式中：nT为视场内温度为某特定值或特定范围内的像元数目；γ为仪器角分辨率，rad/像元；d为测试距离，m。可根据需要，提取热像中目标的温度分布特征和目标图像特征。
1.4 测试结果
1.4.1 测试数据
1）光谱辐射计测试结果（迎头）
按照上述测试方法和数据处理方法，对光谱辐射计测试数据进行处理，得到飞机在中长波红外的辐射强度随时间的关系如图 2所示。
![]()
图 2 迎头测试中长波红外辐射强度随时间的关系图Figure 2. Plot of middle-long wave infrared radiation intensity over time in head-on test由图 2可知，测试时，3~5 μm的光谱辐射强度和8~14 μm的光谱辐射强度值均较低，这是由于测试时，机身遮挡，导致尾焰的辐射强度大部分被遮挡。
2）中红外热像仪测试结果（迎头）
按照1.3.2测试方法和数据处理方法，对中红外热像仪测试数据进行处理，得到飞机在中红外的特征辐射面积时间变化关系如图 3所示。
![]()
图 3 中波红外特征辐射面积随时间的变化关系图Figure 3. Middle wave infrared characteristic radiation area over time由图 3可知，中红外热像仪测试得到辐射面积最大为3.68 m2，且辐射面积随时间变化。这是由两种原因造成的，第一是飞行员操作时，不可能保证飞机发动机一直处于同一种工作状态，而是在某一工作状态附近波动；第二是由于气候原因，由于风速、风向等时刻变化，导致发动机尾焰在与大气环境传热传质过程中一直处于变化状态，从而导致特征辐射面积发生变化。
3）长波红外热像仪测试结果（迎头）
按照1.3.2测试方法和数据处理方法，对长波红外热像仪测试数据进行处理，得到飞机在长波红外的辐射面积随时间变化关系如图 4所示。
![]()
图 4 长波红外辐射面积随时间的变化关系图Figure 4. Long wave infrared characteristic radiation area over time由图 4可知，长波红外热像仪测试得飞机的最大红外特征辐射面积约为5.4 m2，比中红外辐射面积高，这也是由于飞机发动机尾焰的大部分红外辐射被飞机机身遮挡，而机身蒙皮在长波红外的辐射特征比中红外的辐射特征明显。
1.4.2 数据分析
用同样的方式测试飞机在同样工作状态下，正测向、尾向的红外辐射特征。将测试数据用同样的方式处理，得到测试结果如表 1所示。
表 1 飞机红外辐射特征测试结果汇总Table 1. Summary of aircraft radiation characteristics test resultsTest direction and angle Wave length Radiation intensity/(W/Sr) Radiation area/m2 Head-on 0° 3 to 5 μm 50.00 3.68 8 to 14 μm 50.00 5.40 Lateral 90° 3 to 5 μm 68.00 25.90 8 to 14 μm 150.00 12.53 Stern 180° 3 to 5 μm 140.00 8.85 8 to 14 μm 80.00 1.05 测试数据可得如下结论：①飞机正侧向的中长波红外辐射面积比迎头和尾后都高，这是由于正侧向时，飞机最大截面与红外热像仪视场垂直，导致红外辐射面积最大；②飞机尾后的中长波红外辐射强度远远高于正侧向和迎头，尤其是中红外辐射强度达到最高，说明发动机工作时，尾焰和尾喷管能显著提高飞机的红外辐射特征，但是辐射面积较小；③由于飞机机身的遮挡作用，飞机的红外辐射强度和面积均最小。
2. 红外干扰弹使用策略
2.1 红外干扰机理
2.1.1 面源红外干扰弹干扰红外导弹机理分析
面源红外干扰弹投放后在载机附近迅速扩散开来形成红外干扰云团，其与被保护载机的红外图像相似，或改变载机的红外图像特征，欺骗红外成像制导导弹，继而诱使红外成像制导导弹偏离被保护的载机。面源红外干扰弹对抗红外成像制导导弹，各个阶段可以采用不同的方式进行干扰[7]。
目标搜索阶段，导弹距离目标通常很远，目标机形成的红外图像很小，处于点目标阶段，当面源红外干扰弹连续投放后，在真目标附近形成大面积红外辐射云团，掩盖目标的红外辐射特征，使红外成像制导导弹导引头无法锁定真目标，面源红外干扰弹起到“隐真示假”的作用。
目标跟踪阶段初期，弹目距离逐渐缩短，目标形成的红外图像逐渐变大，且逐渐显示目标的外形特征，此时，红外成像制导导弹使用形心跟踪或是质心跟踪方式跟踪目标。当连续投放面源红外干扰弹后，在目标附近形成大面积红外辐射云团，其红外特征与目标相似，并与目标融合在一起形成共同的目标信息，甚至完全掩盖目标的红外辐射特征，从而改变导弹视场内的目标红外辐射特征，常用的灰度、面积、长宽比、圆度等识别特征均发生较大变化，制导系统难以识别出真目标，此时，实际目标的特征或运动参数将以不稳定或“非目标”的方式变化（如机动逃逸），而面源红外干扰弹的特征保持稳定（持续燃烧或连续投放），此时，成像导弹的识别系统将会锁定假目标，而达到干扰目的。
目标跟踪阶段后期，弹目距离很近，目标在导引头视场内外形特征清楚，导引头可根据目标的外形特征（面积、长宽比、圆度等）对目标进行跟踪。此时，连续投放面源红外干扰弹后，面源红外干扰弹形成大面积红外辐射云团，与目标红外辐射特征融合，增大了整个红外辐射云团的辐射面积，改变目标的外形特征，增大导引头识别算法的运算量和识别难度，使其难以识别真目标，导致成像导引头无法通过目标的面目标特征识别出目标，从而有效降低其命中目标概率。
2.1.2 点源红外干扰弹干扰红外导引头机理分析
点源红外干扰弹发射后形成单个热点目标，连续投放后形成多个点目标，通过质心干扰的方式诱骗红外制导导引头哦。由于点源红外干扰弹能量高，远大于飞机的红外辐射强度，在对抗早期第一代、第二代红外制导导引头时有效；随着具有辐射强度阈值鉴别能力的第三代、第四代红外制导导引头的装备，点源红外干扰弹由于能量太高，且干扰源数量少，干扰效能不高的缺点越来越明显。因此，点源红外干扰弹逐步被多点源红外干扰弹取代。
2.1.3 多点源红外干扰弹干扰红外导弹机理分析
多点源红外干扰弹发射后形成多个分散的点源红外干扰弹，连续投放后形成多个点状发热体。采用多发齐射或多方位齐射时，可迅速在一定空域形成红外高辐射区，在导引头瞬时视场内形成持续的多个干扰源，将目标信号淹没，且多点源红外干扰弹的辐射能量比点源红外干扰弹的小，略高于飞机，导引头无法通过强度阈值鉴别滤除干扰，导引头必须处理多组脉冲信号，降低了其检测目标的概率，红外导引头即使启动了抗干扰措施，但因探测器的噪声几何级数增大，而难以提取有效的制导信号，从而起到保护载机的作用[8]。
2.2 干扰策略分析
根据飞机发动机工作时，飞机不同方向上的红外辐射特征数据，以及多点源红外干扰弹、面源红外干扰弹的特点，可分析出飞机在面临不同方向成像红外导弹导引头时红外干扰弹的使用策略[9]。
1）当来袭红外制导导弹从尾后攻击飞机时，开始导弹距飞机较远，飞机在红外导引头中显示为一个热点，此时采用多点源红外干扰弹，以质心干扰的方式诱骗红外导引头；随着来袭导弹与飞机距离缩短，飞机在导引头中面积逐渐增大，采用面源红外干扰弹，其红外辐射云团与飞机融合，改变目标的红外辐射特征和在导引头视角下的面目标特征，使导引头分辨不出目标，从而形成有效干扰。多点源红外干扰弹和面源红外干扰弹组合使用，两者结合可有效干扰红外导引头。
2）当来袭红外制导导弹侧向攻击飞机时，由于飞机在红外导引头中显示的面积较大，且能量较高，此时可采用面源红外干扰弹有效掩盖目标的红外辐射特征，干扰来袭导弹导引头。
3）当来袭红外制导导弹迎头攻击飞机时，由于此时飞机在红外导引头中显示的面积和能量均很低，此时采用多点源红外干扰弹和面源红外干扰弹组合使用，通过多点源红外干扰弹的质心干扰，形成假目标，结合面源红外干扰弹的大面积特征，掩盖目标的红外辐射特征，隐藏真实目标特征，点和面结合，达到“隐真示假”的效果，可有效干扰红外制导导弹导引头。
3. 结论
本文以飞机为研究对象，采用光谱辐射计、中长波红外热像仪测试了飞机发动机工作时，飞机迎头、正侧向、尾向3个不同方向上的红外辐射特征，根据红外辐射特征，结合多点源红外干扰弹、面源红外干扰弹的特点，分析了飞机在面临不同方向来袭红外制导导弹时干扰弹的投放策略。当来袭导弹从尾后攻击时，远距时采用多点源红外干扰弹、中近距时采用面源红外干扰弹；侧向攻击时，采用面源红外干扰弹；迎头攻击时，采用点源、面源组合使用，可形成有效干扰。
-
![]()
图 1 反射式持续热激励红外热像技术实验装置
Figure 1. Schematic diagram of the reflective continuous-heat excited thermography
![]()
图 2 玻璃钢平底洞试件：(a) 试件正面；(b) 试件背面
Figure 2. The GFRP specimen with flat-bottom holes: (a) The front of GFRP specimen; (b)The back of GFRP specimen
![]()
图 3 试件热激励面热图：(a) 热激励0 s热图；(b) 热激励60 s热图；(c) 热激励120 s热图；(d) 热激励180 s热图；(e) 热激励240 s热图；(f) 热激励300 s热图
Figure 3. The thermal image of GFRP specimen with flat-bottom holes: (a) The thermal image at 0 seconds of thermal excitation; (b) The thermal image at 60 seconds of thermal excitation; (c) The thermal image at 120 seconds of thermal excitation; (d) The thermal image at 180 seconds of thermal excitation; (e) The thermal image at 240 seconds of thermal excitation; (f) The thermal image at 300 seconds of thermal excitation
![]()
图 4 温差-时间最小二乘法多项式拟合曲线与温差-时间一阶微分实验曲线：(a) 温差-时间最小二乘法多项式拟合曲线；(b) 温差-时间一阶微分实验曲线
Figure 4. The least squares polynomial fitting curves of ΔT−t and the first order differential experimental curves of ΔT−t: (a) The least squares polynomial fitting curves of ΔT−t; (b) The first order differential experimental curves of ΔT−t
![]()
图 6 温差-时间的拟合曲线和实验曲线
Figure 6. The fitting curves and experimental curves of temperature increment-time
表 1 风电叶片玻璃钢热性能参数
Table 1 Thermal parameters of wind turbine blade GFRP
ρ/(kg·m−3) C/(J·(kg·K) −1) λ/(W·(m·K) −1) α/(m2·s−1) 1830 950 0.38 $ 2.19 \times {10^{ - 7}} $ 
下载: 导出CSV
表 2 各平底洞中心对应热激励面的特征时间
Table 2 The characteristic time of the thermal excitation surface corresponds to the center of each flat-bottom hole
Holes A B C D E F Characteristic time/s 150 200 146 165 205 149
下载: 导出CSV
表 3 各平底洞内厚度测量值及相对误差
Table 3 The measured thickness and relative error of each flat-bottom hole
Holes Thickness
/mmMeasured
thickness/mmRelative errors A 14 13.8 1.4% B 19 17.6 7.4% C 10 9.3 7.0% D 16 15.3 4.4% E 20 18.1 9.5% F 11 10.0 9.1%
下载: 导出CSV
-
[1] 李晓丽. 红外热波定量测量技术的研究及其应用[D]. 北京: 北京理工大学, 2018: 1-10. LI Xiaoli. Studies for Quantitatively Measuring Technique by Using Thermal Wave Imaging and Their Applications[D]. Beijing: Beijing Institute of Technology, 2018: 1-10.
[2] Worzewski T, Krankenhagen R, Doroshtnasir M, et al. Thermographic inspection of a wind turbine rotor blade segment utilizing natural conditions as excitation source, Part Ⅰ: Solar excitation for detecting deep structures in GFRP[J]. Infrared Phys. Technol. , 2016, 76: 756-766. DOI: 10.1016/j.infrared.2016.04.011
[3] Kalyanavalli V, Mithun P M, Sastikumar D. Analysis of long-pulse thermography methods for defect depth prediction in transmission mode[J]. Meas. Sci. Technol. , 2019, 31: 014002.
[4] Shepard S M, Lhota J R, Rubadeux B A, et al. Reconstruction and enhancement of active thermographic image sequences[J]. Opt. Eng. , 2003, 42(5): 1337-1342. DOI: 10.1117/1.1566969
[5] SUN J G. Pulsed thermal imaging measurement of thermal properties for thermal barrier coatings based on a multilayer heat transfer mode[J]. J. Heat Trans. , 2014, 136: 081601. DOI: 10.1115/1.4027551
[6] SUN J G. Thermal conductivity measurement for thermal barrier coatings based on one- and two-sided thermal imaging methods[J]. Review of Quantitative Nondestructive Evaluation, 2010, 29: 458-465.
[7] SUN J G. Quantitative three-dimensional imaging by thermal tomography method[J]. Review of Quantitative Nondestructive Evaluation, 2011, 30: 430-437.
[8] LI Xiaoli, TAO Ning, SUN J G, et al. Thickness measurement by two-sided step-heating thermal imaging[J]. Rev. Sci. Instrum., 2018, 89: 104902. DOI: 10.1063/1.5047049
[9] LI Xiaoli, SUN J G, TAO Ning, et al. An effective method to inspect adhesive quality of wind turbine blades using transmission thermography[J]. J. Nondstrct. Eval. , 2018, 37(19): 1-11
[10] LI Xiaoli, TAO Ning, SUN Jiangang, et al. Evaluation of an ancient cast-iron Buddha head by step-heating infrared thermography[J]. Infrared Phys. Technol. , 2019, 98: 223-229. DOI: 10.1016/j.infrared.2019.03.017
[11] 金学元. 持续热激励红外热像无损检测技术的研究及应用[D]. 北京: 首都师范大学, 2013: 12-18. JIN Xueyuan. Study on the Application of Non-destructive Testing Technology of Infrared Thermal Image with Continuous Thermal Excitation[D]. Beijing: Capital Normal University, 2013: 12-18.
[12] 李晓丽, 孙建刚, 陶宁, 等. 非线性拟合方法用于透射式脉冲红外技术测试碳/碳复合材料的热扩散系数[J]. 物理学报, 2017, 66(18): 188702. https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201718031.htm LI Xiaoli, SUN Jiangang, TAO Ning, et al. Application of nonlinear data fitting method to thermal diffusivity of carbon-carbon composite measured by transmission pulsed thermography[J]. Acta Physica Sinica, 2017, 66(18): 188702. https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201718031.htm
[13] 曾智, 陶宁, 冯立春, 等. 缺陷尺寸对红外热波技术缺陷深度测量的影响研究[J]. 红外与激光工程, 2011, 41(7): 1910-1915. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201207041.htm ZENG Z, TAO N, FENG L C, et al. Effect of defect size on the measurement of defect depth using thermal wave imaging[J]. Infrared and Laser Engineering, 2011, 41(7): 1910-1915. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201207041.htm
-
期刊类型引用(2)
1. 高程，唐超，童安炀，王文剑. 基于CNN和LSTM混合模型的红外人体行为识别. 合肥学院学报(综合版). 2023(05): 77-85 . 
百度学术
2. 赵普，武一. 面向社区医疗的跌倒检测算法. 中国医学物理学杂志. 2023(12): 1486-1493 . 百度学术
其他类型引用(12)
计量
- 文章访问数: 75
- HTML全文浏览量: 26
- PDF下载量: 1
- 被引次数: 14
