一种复合材料冲击损伤红外热成像检测特征提取方法

1.本发明涉及热成像检测技术领域，具体涉及一种复合材料冲击损伤红外热成像检测特征提取方法。


背景技术：

2.随着碳纤维增强环氧树脂基复合材料制造工艺的提升，碳纤维复合材料的比强度与比模量逐渐增加，在飞机结构中的应用也由非承力结构、次承力结构逐步过渡到承力结构。在20世纪60年代，碳纤维增强环氧树脂基复合材料最先应用在整流罩、舱门等非承力结构部位，到了20世纪末，碳纤维复合材料在垂尾、平尾、进气道等次承力结构部位得到广泛应用，再到现阶段，某些飞机的机身、中央翼盒等主承力结构部位也采用了碳纤维复合材料代替原来的金属材料。
3.碳纤维复合材料由于基体采用环氧类热固性树脂，韧性低、脆性大，使得碳纤维复合材料在制造与服役过程中容易产生冲击损伤，造成碳纤维复合材料的剩余强度急剧下降，严重威胁飞机的正常飞行。目前碳纤维复合材料冲击损伤缺陷一般采用超声检测技术和红外热成像技术对其实施无损检测，超声检测需要耦合剂，检测工艺复杂，检测效率低，而红外热成像检测由于检测工艺简单，一次成像面积大，非接触无污染，广泛应用在碳纤维复合材料无损检测项目中。
4.飞机碳纤维复合材料红外热成像无损检测技术采用主动式检测，由热激励源给被检对象主动加热，利用红外热像仪检测工件表面的温度变化来识别缺陷信息，而在实际检测过程中，由于受加热不均影响，碳纤维复合材料冲击损伤红外热图常出现信噪比不高、缺陷轮廓难识别，缺陷定量不准等问题。
5.如在中国申请号为cn202110357421.x；公布日为2021.06.18的专利文件公开了一种碳纤维层压板红外热成像缺陷尺寸检测方法，其包括以下步骤：s1、制作模拟缺陷试块；s2、对检测工件进行清洁；s3、调整和设置红外热成像检测系统；s4、选择缺陷最优图像进行保存；s5、确定缺陷的轮廓像素个数；s6、得到修正系数与深度的关系曲线；s7、确定待测件缺陷的轮廓像素个数；s8、基于对数二阶微分峰值法得到待测件的缺陷深度；s9、根据拟合的关系曲线读取修正系数；s10、计算得到待测件的缺陷尺寸。该发明采用与实际产品材料相同的基础材料制作缺陷模拟试块，保证了缺陷尺寸定量过程的一致性；本发明通过对不同深度缺陷的基准标定，降低了由于横向热扩散对尺寸定量的影响，极大地提高了缺陷定量的精度。
6.该文献的检测方法，其只是对检测工件进行红外热成像检测后，根据特征图像选择缺陷最优图像来进行缺陷尺寸的检测，而最优图像选择的原理是，通过区分缺陷区域与完好区域的温度差异来识别缺陷区域，选择缺陷最优显示时刻时采集到的图像，在现有针对复合材料宫外热图信噪比不高，缺陷特征提取不准确的问题，目前主要采用对比度增强与滤波处理、阈值分割、区域增长、主成分分析等单帧热图的图像处理方法，所有这些研究都是通过主观肉眼分辨出复合材料冲击损伤热图序列中的最优热图，并没有提出如何在整
个热图序列中选取一帧量化性能指标的最优热图，这样，容易降低对复合材料冲击损伤缺陷信息特征提取准确率。


技术实现要素：

7.本发明提供一种复合材料冲击损伤红外热成像检测特征提取方法；通过本发明的方法，可以解决在检测过程中，复合材料由于受加热不均影响而出现缺陷轮廓难识别，缺陷定量不准等问题，能够有效提高复合材料冲击损伤缺陷信息特征提取的准确率。
8.为达到上述目的，本发明的技术方案是：一种复合材料冲击损伤红外热成像检测特征提取方法，具体步骤包括：（1）通过热成像检测装置对试块进行检测，预设采集频率，按照预设采集频率采集试块的热成像信息图，并组成热成像信息图序列。
9.（2）选取具有冲击损伤区域信息量的热成像信息图，选取冲击损伤区域从开始出现到逐渐消失过程的两个以上热成像信息图形成最大热成像信息图序列，然后获取最大热成像信息图序列温差并形成温差时间历程曲线，根据温差时间历程曲线生成温差时间历程曲线最大梯度值，根据最大梯度值对应的帧数的热图为最大的热成像信息图。
10.（3）对最大的热成像信息图进行预处理，提升图像对比度。
11.（4）对最大的热成像信息图进行后处理。
12.（5）得出试块缺陷区域信息，即试块缺陷区域的冲击损伤面积、损伤长度、损伤宽度、最大损伤直径特征信息。
13.上述方法，通过热成像检测装置对复合材料试块冲击损伤情况进行检测，并选取两个以上热成像信息图形成最大热成像信息图序列，根据试块上有冲击损伤区域的温差-时间历程曲线，并根据温差时间历程曲线确定其中最大梯度值对应的帧数从而确定为最大热成像信息图，然后对最大热成像信息图进行图像处理后，提取对应的损伤参数，由于其是根据形成的温差时间历程曲线形成最大梯度值确定的最大热成像信息图，从而使得最大热成像信息图是实际受到损坏的区域，然后通过对图像进行处理然后提取对应的参数，从而能确保如何提取复合材料冲击损伤特征的单针最优热图提出了量化指标，可以有效提高复合材料冲击损伤缺陷信息特征提取的准确性。
14.进一步的，步骤（1）具体包括通过卤素灯发射热激励源照射试块3s-5s，采集频率为5-15hz，共采集900-1100帧热成像信息图。
15.进一步的，步骤（2）具体包括选择能够完整显示冲击损伤区域从出现到逐渐消失过程两个以上的热成像信息图帧数。以上方法能使得选择的热成像信息图帧数能够完整显示冲击损伤区域的整个过程，使得提取方法更加准确。
16.进一步的，步骤（3）具体包括：（3.1）对最大的热成像信息图进行灰度变换。
17.（3.2）采用直方图均衡化处理方式对最大的热成像信息图进行图像预处理，提升图像对比度。
18.（3.3）对直方图均衡化后的图像采用高通滤波处理。
19.以上方法，通过进行灰度变化之后进行滤波处理，使得热成像信息图能较小噪声影响，确保提取特征值的准确性。
20.进一步的，步骤（4.2）具体包括：通过钝化掩蔽对缺陷区域边缘进行处理。
21.所述钝化掩蔽具体步骤包括：（4.21）模糊缺陷区域图像得到模糊后的区域图像。
22.（4.22）从原区域图像中减去模糊后的区域图像。
23.（4.23）取原区域图像与模糊后的区域图像相加为模板。
24.（4.24）将模板与原区域图像相加；得出钝化掩蔽后的区域图像。
25.以上方法，通过钝化掩蔽处理对缺陷区域边缘进行处理，进一步提升图像的准确性。
26.进一步的，步骤（1）具体包括通过卤素灯发射热激励源照射试块4s，采集频率为10hz，共采集1000帧热成像信息图。
27.以上方法，通过采用该参数，能使得形成的热成像信息图的帧数最多。
28.进一步的，步骤（2）具体包括若在预设的采集频率内没有冲击损伤区域信息量的热成像信息图时，则增大采集频率。
29.以上方法，在没有冲击损伤区域信息量时，增大采集频率进一步扩大采集帧数，确保采集的准备性。
30.进一步的，步骤（2）具体还包括若在增大采集频率内仍然没有冲击损伤区域信息量的热成像信息图时，则判断该试块没有受到冲击损伤；若出现冲击损伤区域，则选取冲击损伤区域从开始出现到逐渐消失过程的两个以上热成像信息图并形成最大热成像信息图序列。
附图说明
31.图1为本发明实施例1的能够完整显示冲击损伤区域从出现到逐渐消失过程的热图。
32.图2为本发明实施例1的感兴趣区域的示意图。
33.图3为本发明的试块温差时间（帧数）历程曲线示意图。
34.图4为本发明实施例1的试块温差时间历程曲线上升沿的梯度值示意图。
35.图5为本发明实施例1的最大热成像信息图。
36.图6为本发明实施例1中经过直方图均衡化处理后的最佳热图。
37.图7为本发明实施例1中经过高通滤波处理后的最佳热图。
38.图8为本发明实施例1中经过钝化掩蔽方法处理后的最佳热图。
39.图9为本发明实施例2的能够完整显示冲击损伤区域从出现到逐渐消失过程的热图。
40.图10为本发明实施例2的感兴趣区域的示意图。
41.图11为本发明实施例2的试块温差时间历程曲线上升沿的梯度值示意图。
42.图12为本发明实施例2的最佳热图。
43.图13为本发明实施例2中经过直方图均衡化处理后的最佳热图。
44.图14为本发明实施例2中经过高通滤波处理后的最佳热图。
45.图15为本发明实施例2中经过钝化掩蔽方法处理后的最佳热图。具体实施方式下面结合附图和具体实施方式对本发明做进一步详细说明。
46.实施例1。
47.如图1-8所示，一种复合材料冲击损伤红外热成像检测特征提取方法，具体步骤包括：（1）通过热成像检测装置对试块进行检测，预设采集时间以及采集频率，按照预设采集频率在照射时间内采集试块的热成像信息图，并组成热成像信息图序列。具体地，预设采集照射时间为3s-5s，预设采集频率是5-15hz。
48.（2）选取具有冲击损伤区域信息量的热成像信息图，选取冲击损伤区域从开始出现到逐渐消失过程的两个以上热成像信息图形成最大热成像信息图序列，然后获取最大热成像信息图序列温差并形成温差时间历程曲线，根据温差时间历程曲线生成温差时间历程曲线最大梯度值，根据最大梯度值对应的帧数的热图为最大的热成像信息图。
49.（3）对最大的热成像信息图进行预处理，提升图像对比度。
50.（4）对最大的热成像信息图进行后处理。
51.（5）得出试块缺陷区域信息，即试块缺陷区域的冲击损伤面积、损伤长度、损伤宽度、最大损伤直径特征信息。
52.通过热成像检测装置对复合材料试块冲击损伤情况进行检测，并选取两个以上热成像信息图形成最大热成像信息图序列，根据试块上有冲击损伤区域的温差-时间历程曲线，并根据温差时间历程曲线确定其中最大梯度值对应的帧数从而确定为最大热成像信息图，然后对最大热成像信息图进行图像处理后，提取对应的损伤参数，由于其是根据形成的温差时间历程曲线形成最大梯度值确定的最大热成像信息图，从而使得最大热成像信息图是实际受到损坏的区域，然后通过对图像进行处理然后提取对应的参数，从而能确保如何提取复合材料冲击损伤特征的单针最优热图提出了量化指标，可以有效提高复合材料冲击损伤缺陷信息特征提取的准确性。
53.步骤（1）具体包括通过卤素灯发射热激励源照射试块3s-5s，采集频率为5-15hz，共采集900-1100帧热成像信息图。在本实施例中，所述试块a的大小为150*100mm，厚度为3.89mm，对试块造成的冲击能力为72.19j，所述热成像检测装置采用四个功率为500w的卤素灯作为热激励源对试块进行照射，卤素灯发射热激励源照射试块4s，采集频率为10hz，共采集1000帧热成像信息图（以下简称为热图）。
54.步骤（2）具体包括选择能够完整显示冲击损伤区域从出现到逐渐消失过程的热成像信息图帧数。在本实施例中，操作者通过观察热图的状态变化，从1000帧热成像信息图中选取第39帧至186帧能够完整显示冲击损伤区域从出现到逐渐消失的过程的热图，如图1所示，39th为第39帧损伤区域的热图，其损伤区域较小，而到第186帧时，其区域增大到较大，且在第39帧-186帧之间，损伤区域处于增大状态。
55.若在预设的采集频率内没有冲击损伤区域信息量的热成像信息图时，则增大采集频率。具体地，增大采集频率为10hz。
56.若在增大采集频率内仍然没有冲击损伤区域信息量的热成像信息图时，则判断该试块没有受到冲击损伤。
57.若出现冲击损伤区域，则选取冲击损伤区域从开始出现到逐渐消失过程的两个以上热成像信息图并形成最大热成像信息图序列。
58.在能够完整显示冲击损伤区域从出现到逐渐消失的过程的每一帧热图中具有冲
击损伤区域和无冲击损伤区域上拾取多个感兴趣区域，拾取的过程是操作者随机选取预设数量的感兴趣区域。在本实施例中，感兴趣区域即冲击损伤区域和无冲击损伤区域内多个区域，操作者在冲击损伤区域拾取8个感兴趣和在无冲击损伤区域捡取4个感兴趣区域，其中冲击损坏区域内8个感兴趣区域沿着类似圆周方向排列，如图2所示。
59.以冲击损伤区域和无冲击损伤区域的感兴趣区域的平均温度求得温差作为纵坐标，即冲击损区域的平均温度和无冲击损伤区域的平均温度之间的差值作为纵坐标，以帧数作为横坐标绘制试块的温差时间（帧数）历程曲线，如图3所示。在本实施例中，所述试块为图3中试块a表达的线段所示。
60.从图3中可以看出，试块在第110帧时达到最大温差，而通过试块的部分热成像信息图如图1可知，按照图1的趋势可知试块a在110帧时显示的冲击损伤面积处于变小状态，不如第110帧前面热图显示的冲击损伤面积大，说明最能展示试块冲击损伤面的热图，即热图序列中的最优热图在最大温差前，因此，能够提取试块冲击损伤特征的最大热成像信息图，应该在温差时间历程曲线的上升沿中。然后对处于上升沿的试块如试块的第39-110帧能够明显看出冲击损伤特征的热成像信息图对应的温差时间历程曲线进行微分，从而求取温差时间历程曲线上升沿的最大梯度值，微分求解公式为：gradu=
∂
u/
∂
x式中，grad为温差时间历程曲线上升沿的梯度值，
∂
u代表相邻帧数之间温度对应的温差值，
∂
x代表热成像信息图的间隔的帧数，比如对于第3帧的梯度值为第3帧的温度与第2帧的温度之间的差值除以第3帧减去第2帧数之间的差值，预设第0帧对应的温度值为0。如图4所示。
61.试块在第60帧时达到梯度最大值，即试块在第60帧时的热成像信息图为最大热成像信息图；如图5所示。
62.在一实施例中，步骤（3）具体包括：步骤（3.1）具体包括先将最大热成像信息图进行灰度变换，通过灰度变换改善最大热成像信息图的画质，使图像的显示效果更加清晰。
63.步骤（3.2）具体包括采用直方图均衡化处理方法对试块的最大热成像信息图进行图像预处理，提升对比度。如图6所示。
64.在本实施例中，所述直方图均衡化处理方法是图像处理领域中常见的利用图像直方图对对比度进行调整的方法，其具体的调整方法为现有技术，在此不再累述。
65.步骤（3.3）具体包括对直方图均衡化后的图像采用高通滤波处理，消除细小的噪音信号，增强缺陷区域与非缺陷区域的对比度。
66.在本实施例中，所述高通滤波处理方法是图像处理领域中常见的用于图像处理的减噪方法，具体为现有技术，在此不再累述。通过高通滤波处理后的图像如图7所示。
67.在一实施例中，步骤（4）具体包括：步骤（4.2）具体包括：通过钝化掩蔽对缺陷区域边缘进行处理。
68.所述钝化掩蔽具体步骤包括：（4.21）模糊缺陷区域图像得到模糊后的区域图像。
69.（4.22）从原区域图像中减去模糊后的区域图像。
70.（4.23）取原区域图像与模糊后的区域图像相加为模板。
71.令f1（x,y）表示模糊后的区域图像，公式为：g
mask
（x,y）=f（x,y）-f1（x,y）式中，f（x,y）代表原区域图像在坐标（x，y）所在位置的灰度值，g
mask
（x,y）为在坐标（x，y）所在位置的模板图像函数。
72.（4.24）将模板与原区域图像相加，得出钝化掩蔽后的区域图像如图8所示。
73.g（x,y）=f（x,y）+ g
mask
（x,y）式中，g（x,y）代表在坐标（x，y）所在位置钝化掩蔽后的区域图像。
74.在一实施例中，步骤（5）具体包括从钝化掩蔽处理后的区域图像上测量出缺陷区域的冲击损伤面积a1、损伤长度a2、损伤宽度a3、最大损伤直径a4特征信息。
75.为了进一步验证本方法的效果，通过测量进行，测量方式可以通过超声波c扫描技术或红外热成像（本实施例采用）的方式实现，通过借助超声波c扫描技术检测试块的冲击损伤面积、损伤长度、损伤宽度、最大损伤直径特征信息，再与本实施例的方法进行比较，计算误差，结果如下表1所示：表1从表中可以看出，整体误差较小，可控，因此，本实施例的方法是可行的。
76.实施例2。
77.如图3、图9-15所示，一种复合材料冲击损伤红外热成像检测特征提取方法，具体步骤包括：（1）通过热成像检测装置对试块进行检测，预设采集时间以及采集频率，按照预设采集频率在照射时间内采集试块的热成像信息图，并组成热成像信息图序列。具体地，预设采集照射时间为3s-5s，预设采集频率是5-15hz。
78.（2）选取具有冲击损伤区域信息量的热成像信息图，选取冲击损伤区域从开始出现到逐渐消失过程的两个以上热成像信息图形成最大热成像信息图序列，，然后获取最大热成像信息图序列温差并形成温差时间历程曲线，根据温差时间历程曲线生成温差时间历程曲线最大梯度值，根据最大梯度值对应的帧数的热图为最大的热成像信息图。
79.（3）对最大的热成像信息图进行预处理，提升图像对比度。
80.（4）对最大的热成像信息图进行后处理。
81.（5）得出试块缺陷区域信息，即试块缺陷区域的冲击损伤面积、损伤长度、损伤宽度、最大损伤直径特征信息。
82.上述方法，通过热成像检测装置对复合材料试块冲击损伤情况进行检测，并选取两个以上热成像信息图形成最大热成像信息图序列，根据试块上有冲击损伤区域的温差-时间历程曲线，并根据温差时间历程曲线确定其中最大梯度值对应的帧数从而确定为最大
热成像信息图，然后对最大热成像信息图进行图像处理后，提取对应的损伤参数，由于其是根据形成的温差时间历程曲线形成最大梯度值确定的最大热成像信息图，从而使得最大热成像信息图是实际受到损坏的区域，然后通过对图像进行处理然后提取对应的参数，从而能确保如何提取复合材料冲击损伤特征的单针最优热图提出了量化指标，可以有效提高复合材料冲击损伤缺陷信息特征提取的准确性。
83.步骤（1）具体包括通过卤素灯发射热激励源照射试块3s-5s，采集频率为5-15hz，共采集900-1100帧热成像信息图。在本实施例中，本实施例与实施例1采用的试块不同，所述试块b的大小为150*100mm，厚度为3.11mm，对试块造成的冲击能力为18.29j，所述热成像检测装置采用四个功率为500w的卤素灯作为热激励源对试块进行照射，卤素灯发射热激励源照射试块4s，采集频率为10hz，共采集1000帧热成像信息图（以下简称为热图）。
84.步骤（2）具体包括选择能够完整显示冲击损伤区域从出现到逐渐消失过程的热成像信息图帧数。在本实施例中，操作者通过观察热图的状态变化，从1000帧热成像信息图中选取第22帧至108帧能够完整显示冲击损伤区域从出现到逐渐消失的过程的热成像信息图，如图9所示，22th为第22帧损伤区域的热成像信息图，其损伤区域较小，而到第108帧时，其区域增大到较大，且在第22帧-108帧之间，损伤区域处于增大状态。
85.若在预设的采集频率内没有冲击损伤区域信息量的热成像信息图时，则增大采集频率。具体地，增大采集频率为10hz。
86.若在增大采集频率内仍然没有冲击损伤区域信息量的热成像信息图时，则判断该试块没有受到冲击损伤。
87.若出现冲击损伤区域，则选取冲击损伤区域从开始出现到逐渐消失过程的两个以上热成像信息图并形成最大热成像信息图序列。
88.在能够完整显示冲击损伤区域从出现到逐渐消失的过程的每一帧热图中具有冲击损伤区域和无冲击损伤区域上拾取多个感兴趣区域，拾取的过程是操作者随机选取预设数量的感兴趣区域。在本实施例中，操作者在冲击损伤区域捡取8个感兴趣和在无冲击损伤区域捡取4个感兴趣区域，如图10所示。
89.以冲击损伤区域和无冲击损伤区域的感兴趣区域的平均温度求得温差作为纵坐标，即冲击损区域的平均温度和无冲击损伤区域的平均温度之间的差值作为纵坐标，以帧数作为横坐标绘制试块的温差-时间（帧数）历程曲线，如图3所示。在本实施例中，所述试块为图3中试块b表达的线段所示。
90.从图3中可以看出，试块b在第110帧时达到最大温差，而通过试块的部分热图如图9可知，按照图9的趋势可知试块b在110帧时显示的冲击损伤面积处于变小状态，不如第110帧前面热图显示的冲击损伤面积大，说明最能展示试块冲击损伤面的热图，即热图序列中的最大热成像信息图在最大温差前，因此，能够提取试块冲击损伤特征的最优热图，应该在温差时间历程曲线的上升沿中。对试块b的第22帧-108帧能够明显看出冲击损伤特征的热成像信息图对应的温差时间历程曲线进行微分，从而求取温差时间历程曲线上升沿的梯度值，微分求解公式为：gradu=
∂
u/
∂
x式中，grad为温差时间历程曲线上升沿的梯度值，
∂
u代表相邻帧数之间温度对应的温差值，
∂
x代表热成像信息图的间隔的帧数，比如对于第3帧的梯度值为第3帧的温度与
第2帧的温度之间的差值除以第3帧减去第2帧数之间的差值，预设第0帧对应的温度值为0。如图11所示。
91.在另一实施例中，在取得最大温差之后结合图9所示的热成像信息图，可知在第80帧时热成像信息图的损失区域已开始变小，从而选择对试块b的第22-80帧能够明显看出冲击损伤特征的热图对应的温差时间历程曲线微分，从而求取温差时间历程曲线上升沿的梯度值，这样通过进一步结合热成像信息图使得计算结构更加准确，微分求解公式为：gradu=
∂
u/
∂
x式中，grad为温差时间历程曲线上升沿的梯度值，
∂
u代表相邻帧数之间温度对应的温差值，
∂
x代表热成像信息图的间隔的帧数，比如对于第3帧的梯度值为第3帧的温度与第2帧的温度之间的差值除以第3帧减去第2帧数之间的差值，预设第0帧对应的温度值为0。如图11所示。
92.试块在第52帧时达到梯度最大值，即试块在第52帧时的热图为最佳热图；如图12所示。
93.在一实施例中，步骤（3）具体包括：步骤（3.1）具体包括先将最最大热成像信息图进行灰度变换，通过灰度变换改善最大热成像信息图的画质，使图像的显示效果更加清晰。
94.步骤（3.2）具体包括采用直方图均衡化处理方法对试块的最大热成像信息图进行图像预处理，提升对比度。如图13所示。
95.在本实施例中，所述直方图均衡化处理方法是图像处理领域中常见的利用图像直方图对对比度进行调整的方法，其具体的调整方法为现有技术，在此不再累述。
96.步骤（3.3）具体包括对直方图均衡化后的图像采用高通滤波处理，消除细小的噪音信号，增强缺陷区域与非缺陷区域的对比度。
97.在本实施例中，所述高通滤波处理方法是图像处理领域中常见的用于图像处理的减噪方法，具体为现有技术，在此不再累述。通过高通滤波处理后的图像如图14所示。
98.在一实施例中，步骤（4）具体包括：步骤（4.2）具体包括：通过钝化掩蔽对缺陷区域边缘进行处理。
99.所述钝化掩蔽具体步骤包括：（4.21）模糊缺陷区域图像得到模糊后的区域图像。
100.（4.22）从原区域图像中减去模糊后的区域图像。
101.（4.23）取原区域图像与模糊后的区域图像相加为模板。
102.令f1（x,y）表示模糊后的图像，公式为：g
mask
（x,y）=f（x,y）-f1（x,y）式中，f（x,y）代表原区域图像在坐标（x，y）所在位置的灰度值，g
mask
（x,y）为在坐标（x，y）所在位置的模板图像函数。
103.（4.24）将模板与原区域图像相加，得出钝化掩蔽后的区域图像如图15所示。
104.g（x,y）=f（x,y）+ g
mask
（x,y）式中，g（x,y）代表在坐标（x，y）所在位置钝化掩蔽后的区域图像。
105.在一实施例中，步骤（5）具体包括从钝化掩蔽处理后的图像上测量出缺陷区域的冲击损伤面积b1、损伤长度b2、损伤宽度b3、最大损伤直径b4特征信息。
106.为了进一步验证本方法的效果，通过测量进行，测量方式可以通过超声波c扫描技术或红外热成像（本实施例采用）的方式实现，通过借助超声波c扫描技术检测试块的冲击损伤面积、损伤长度、损伤宽度、最大损伤直径特征信息，再与本实施例的方法进行比较，计算误差，结果如下表2所示：表2从表中可以看出，整体误差较小，可控，因此，本实施例的方法是可行的。