复合材料板在变温场下的超声导波传播特性及损伤成像方法

王长林, 任利, 钟永腾

王长林, 任利, 钟永腾. 复合材料板在变温场下的超声导波传播特性及损伤成像方法[J]. 红外技术, 2022, 44(12): 1338-1343.
引用本文: 王长林, 任利, 钟永腾. 复合材料板在变温场下的超声导波传播特性及损伤成像方法[J]. 红外技术, 2022, 44(12): 1338-1343.
WANG Changlin, REN Li, ZHONG Yongteng. Study on Ultrasonic Guided Wave Propagation Characteristics and Damage Imaging for Composite Structures Under Variable Temperature Field[J]. Infrared Technology , 2022, 44(12): 1338-1343.
Citation: WANG Changlin, REN Li, ZHONG Yongteng. Study on Ultrasonic Guided Wave Propagation Characteristics and Damage Imaging for Composite Structures Under Variable Temperature Field[J]. Infrared Technology , 2022, 44(12): 1338-1343.

复合材料板在变温场下的超声导波传播特性及损伤成像方法

基金项目: 

国家自然科学基金青年科学基金项目 51505339

浙江省自然科学基金青年基金项目 LQ16E050005

详细信息
    作者简介:

    王长林(1983-),男,江苏盐城人,硕士,讲师,主要从事目标检测方面的研究工作。E-mail:wangcl@jsou.edu.cn

    通讯作者:

    钟永腾(1984-),男,福建龙岩人,博士,副教授,主要从事结构健康监测技术方面的研究工作。E-mail:zhongyongteng@wzu.edu.cn

  • 中图分类号: V219

Study on Ultrasonic Guided Wave Propagation Characteristics and Damage Imaging for Composite Structures Under Variable Temperature Field

  • 摘要: 基于Lamb波的损伤监测方法已在复合材料结构健康监测中得到了广泛的应用。然而,复合材料结构的服役环境复杂多变,超声导波信号极易被外界因素干扰。本文研究了玻璃纤维层合板在变温场下的超声导波特性及损伤成像修正方法。首先,利用红外热成像仪研究持续加热下超声导波信号幅值和相位的变化规律,并构建因变温引起的幅相误差矩阵。其次,将误差矩阵代入多重信号分类(2D-MUSIC)算法中修正导向矢量,并建立代价函数迭代出准确损伤位置。玻璃纤维复合材料层合板的模拟损伤实验表明该方法有效提高了2D-MUSIC算法在变温场下的损伤定位分辨率和精度。
    Abstract: The Lamb-wave-based damage location method has been widely used for health monitoring of composite structures. However, it is easily interfered by external factors because its service environment is complex and changeable. To determine the influence of the temperature field on the propagation of Lamb waves in composite structures, in this study we first investigated the propagation characteristics of ultrasonic guided waves on glass fiber laminates under a temperature field using an infrared thermal imager. Subsequently, a corrected model was established using the extracted amplitude attenuation and phase delay errors. Consequently, a modified multiple signal classification(2D-MUSIC) algorithm-based damage imaging method is proposed for composite structures. The experimental results on glass fiber composite laminates show that the proposed method can effectively improve the resolution and accuracy of the original algorithm under a variable temperature field.
  • 红外成像具有非接触直接测量输变电设备状态的特征,无需停电,在电力领域具有广泛的应用[1-3]。输变电设备在野外严峻环境中长期运行,由于材料老化、污秽闪络、机械受损等因素,常伴随有局部放电、温度增高等现象。绝缘子作为输电线路重要部件,起到电气隔离和机械支撑的作用,当红外图像显示的温度较高,表明其有异常缺陷,这是因为场强不均匀发生局部放电,严重时可能导致线路故障,甚至造成停电故障。文献[4-5]开展了变压器高压套管的红外诊断研究,通过热像特征谱图分析了故障原因。红外图像可应用于交流滤波器的故障分析中,提取其跳闸的典型故障特征[6-8]

    输变电设备红外图谱特征提取主要为图像处理方法,诸如纹理、色彩、边框等图像特征[9-10]。随着无人机航拍技术的发展,输变电设备红外图谱呈现爆发式增长,深度学习提供了一种良好的手段[11-13],采用卷积神经网络(Convolutional Neural Networks,CNN)的方法,对海量图片进行训练学习,提取特征进行测试和验证。CNN通过卷积层计算,输入较深层的特征图,对于小目标绝缘子、套管,权重值相对较少,无法实现小目标的有效提取[14-15]。针对这一缺点,本文对Faster R-CNN方法进行改进,提高绝缘子红外图谱诊断的精度。

    区别于普通的CNN,Faster R-CNN增加了一个区域提取网络(Region Proposal Network,RPN),即图 1中的候选区域,摒弃传统的滑动窗口,可在GPU直接运行计算,极大地加快了计算速度。RPN判断每个像素点对应的多个不同尺度和宽高比的锚框是否为前景目标的二分类,形成候选区域。

    图  1  Faster R-CNN的算法流程
    Figure  1.  Algorithm flow of Faster R-CNN

    Faster R-CNN一般采用随机梯度下降法(Stochastic Gradient Descent,SGD)训练神经网络,见式(1):

    $$ h(x) = \sum\limits_{i = 0}^n {{w_i}} {x_i} = {\mathit{\boldsymbol{W}}^{\rm{T}}}\mathit{\boldsymbol{X}} $$ (1)

    式中:X为输入;W为权重;wixi分别表示第i个权重和输入;h(x)为对应的输出。

    损失函数S(W)基于平方误差实现,见式(2):

    $$ S(\mathit{\boldsymbol{W}}) = \frac{1}{2}\sum\limits_{i = 0}^m {{{({h_W}({x_i}) - {y_i})}^2}} $$ (2)

    式中:yi为真实输出。

    W的更新函数见式(3):

    $$ {W_j} = {W_j} - \alpha \frac{\partial }{{\partial {W_j}}}S(\mathit{\boldsymbol{W}}) $$ (3)

    式中:α为学习率,可设置步长。W通过梯度下降法进行求解,首先正向计算样本输出值,接着根据反向传递的误差迭代计算,常用在CNN训练中。

    RPN的选择本质是通过平移或者尺度变化的方法将合适的区域提取网络R变成$\hat C$,从而接近实际的候选框C

    $$ f({R_x}, {R_y}, {R_w}, {R_h}) = \left( {{{\hat C}_x}, {{\hat C}_y}, {{\hat C}_w}, {{\hat C}_h}} \right) \approx \left( {{C_x}, {C_y}, {C_w}, {C_h}} \right) $$ (4)

    式中:(x, y)、(w, h)分别为矩形区域的中心坐标和宽、高。

    令${t_*}$为矩形区域的平移和缩放量,则有:

    $$ \left\{ {\begin{array}{*{20}{l}} {{t_x} = ({C_x} - {R_x})/{R_{\rm{w}}}}\\ {{t_y} = ({C_y} - {R_y})/{R_{\rm{h}}}}\\ {{t_{\rm{w}}} = \lg ({C_{\rm{w}}}/{R_{\rm{w}}})}\\ {{t_{\rm{h}}} = \lg ({C_{\rm{h}}}/{R_{\rm{h}}})} \end{array}} \right. $$ (5)

    式中:txty为矩形区域的中心坐标平移量;twth分别为矩形区域的宽、高的缩放量。

    预测值计算过程为:

    $$ {d_*}(R) = w_*^{\rm{T}}\phi (R) $$ (6)

    式中:ϕ是最后一次卷积计算。

    损失函数的目标值计算见式(7),通过调整平移和缩放的尺度,确定最终的候选区域[16]

    $$ {L_{{\rm{oss}}}} = \sum\limits_{i = 1}^N {{{(t_*^i - w_*^{\rm{T}}\phi ({C^i}))}^2}} $$ (7)

    为了增强Faster R-CNN的小目标特征提取能力,引入压缩激励结构,即压缩和激励两大操作,设图像的特征参数设置为(H, W, K),分别表示为长、宽和通道数。

    压缩操作Fsq(·)基于各个通道实现特征图空间信息的压缩,见式(8):

    $$ {h_c} = {F_{{\rm{sq}}}}\left( {{k_c}} \right) = \frac{1}{{H \times W}}\sum\limits_{i = 1}^H {\sum\limits_{j = 1}^W {{k_c}(i, j)} } $$ (8)

    式中:kc表示第c个通道;hc表示经过压缩后输出向量h的第c个元素。

    激励操作分为激励Fex(·)和校准Fscale(·)两个过程,分别见式(9)和式(10):

    $$ s = {F_{{\rm{ex}}}}(h, w) = \sigma (g(z, w)) = \sigma ({w_2}\delta ({w_1}z)) $$ (9)

    式中:σ是sigmoid激活函数;w1为$\frac{c}{r} \times \vec C$的实数矩阵,表示通道的缩减,r为缩减因子;δ是ReLU激活函数,w2为$\vec C \times \frac{c}{r}$的实数矩阵,表示通道的恢复。

    $$ {\tilde h_c} = {F_{{\rm{scale}}}}({h_c}, {s_c}) = {s_c}.{h_c} $$ (10)

    式中:sc表示激活向量s的第c个元素;${\tilde h_c}$表示校准后的对应元素。

    本文基于改进的Faster R-CNN方法,对平台的环境搭建要求较高,具体配置见表 1。操作系统为开源Linux,数据库为MySQL;硬件配置较高,CPU采用Intel高端系列,内存和硬盘容量均较大,保证大量数据的高效运算。框架采用2018年初公开的目标检测平台Detectron,包含最具代表性的目标检测、图像分割、关键点检测算法。

    表  1  软硬件配置
    Table  1.  Hardware and software configuration
    Name Model
    Operating system Ubuntu 16.04.1
    Database mysql 5.5.20
    CPU Intel Xeon Silver 4114T 12C
    GPU NVIDIA GTX1080Ti
    Memory 32 G
    Hard disk 1 T
    Frame Detectron
    下载: 导出CSV 
    | 显示表格

    图像数据来源于多条输电线路无人机拍摄的大量绝缘子照片。在神经网络的训练过程中,对正负样本的判定见图 2,主要基于锚框映射图与真实目标框的交并比(Intersection over Union,IoU)来进行计算分类。首先对RPN形成的锚框进行排序筛选形成锚框序列,接着利用边框回归参数向量修正锚框的位置形成候选区域集合,然后计算所有感兴趣区域(Region of Interest,RoI)与真实目标框的IoU,求最大值,并判断其是否大于0.5,若满足,则为正样本,否则为负样本。

    图  2  正负样本判定
    Figure  2.  Positive and negative sample decision

    完成正负样本的判定后,为了使样本的采样尽量均衡,保证双方的训练集和验证集数量一致,同时采用迁移学习的方法,经过相关修正与补偿,扩充样本总量至2375,样本信息配置见表 2

    表  2  样本配置信息
    Table  2.  Information of sample configuration
    sample type training set verification set test set total
    positive 500 250 750 1500
    negative 500 250 125 875
    total 1000 500 875 2375
    下载: 导出CSV 
    | 显示表格

    普通的CNN方法,原始图像经过卷积层和池化层后,全链接层输出结果,本文方法的结构如图 3所示,引入压缩激励的过程,压缩特征图的空间信息,并通过激励操作学习通道间的依赖关系,可自适应分配每个通道的权重值,提取有利于任务的重要特征通道,最终能进一步增强网络模型的特征提取能力,采用SE-DenseNet-169框架的Faster R-CNN模型。

    图  3  改进的Faster R-CNN结构
    Figure  3.  Structure of improved Faster R-CNN

    改进模型主要完成绝缘子异常状态的精准识别,首先对红外原始图像进行相关修正与补偿实现样本扩充,然后采用本文方法进行训练,收敛后,获得最终的改进Faster R-CNN模型。

    CNN学习中,精确度的衡量一般会采用准确率(Precision)和召回率(Recall),其计算过程分别见式(11)和式(12):

    $$ {P_{{\rm{re}}}} = \frac{{{\rm{TP}}}}{{{\rm{TP}} + {\rm{FP}}}} $$ (11)
    $$ {P_{{\rm{ca}}}} = \frac{{{\rm{TP}}}}{{{\rm{TP + FN}}}} $$ (12)

    式中:TP表示是实际值和预测值均是异常绝缘子的个数;FP表示预测值是异常绝缘子,实际值却不是的个数;FN表示是实际值是异常绝缘子,预测值却不是的个数。

    为进一步衡量改进模型的优劣,这里采用平均检测精度(mean Average Precision,mAP),其中AP等价于召回率和准确率形成曲线与横轴包围的几何图形的面积,对所有类别的AP求平均值即可得到mAP。

    基于样本数据,开展BP、Faster R-CNN以及本文方法的精确度和效率对比研究,不同方法的实验结果统计见表 3。Faster R-CNN和本文方法在Precision方面均明显优于BP方法,本文方法的Recall最高,mAP也最高,相对于BP提高了近10%,这说明经过改进的Faster R-CNN对于小目标的特征提取具有明显的优势。本文方法通过压缩激励结果,减少了数据量的计算,所以相对于其他方法,有更高的效率。

    表  3  不同方法的实验结果统计
    Table  3.  Statistics of experimental results by different methods
    Name Precision Recall mAP Time/s
    BP 93.5% 90.4% 80.3% 2.3
    Faster R-CNN 98.7% 95.3% 88.7% 1.2
    BFEM 99.2% 97.6% 90.2% 0.9
    下载: 导出CSV 
    | 显示表格

    绘制其准确率-召回率关系曲线,如图 4所示,可更加形象直观地反映出本文方法对绝缘子异常特征的提取优势,因为另外两种方法的曲线均被完整的覆盖,说明本文方法改进效果明显。

    图  4  准确率-召回率关系曲线
    Figure  4.  Relation curves of precision and recall

    常见绝缘子排列有单Ⅰ型、双Ⅰ型以及Ⅴ型。本文开展这3类绝缘子的红外图像研究,比较不同排列方式的诊断准确率,如图 5所示。根据电力标准DL/T 664-2008[19],图(b)和图(c)绝缘子端部明显发热,属于异常情况。

    图  5  不同类型绝缘子的红外图像
    Figure  5.  Infrared image of different types of insulators

    不同类型绝缘子的异常诊断准确率见表 4,准确率均较高,均在90%以上;Ⅰ型和Ⅴ型绝缘子的准确率明显优于双Ⅰ型绝缘子,这是因为双Ⅰ型绝缘子会出现两排绝缘子重叠的情况,对红外图像的研究造成一定的影响,为此无人机对于该种类型绝缘子的线路需开展多种角度的拍摄。

    表  4  绝缘子异常诊断的准确率
    Table  4.  Accuracy of insulator anomaly diagnosis
    Insulator type Abnormal total Detected number Accuracy
    Single Ⅰ 62 61 98.4%
    Double Ⅰ 47 44 93.6%
    31 31 100.0%
    下载: 导出CSV 
    | 显示表格

    本文提出一种改进的Faster R-CNN方法,引入激励压缩环节,搭建训练模型,完成绝缘子红外图像的异常诊断,并成功应用于电力现场运维。本文方法可高效并精准地识别出绝缘子的异常缺陷,mAP达到90.2%。研究结果可为输电线路绝缘子缺陷识别研究提供一定的参考。

  • 图  1   变温下阵列信号传播模型示意图

    Figure  1.   Array signal model under temperature fields

    图  2   变温环境下补偿2D-MUSIC算法流程

    Figure  2.   2D-MUSIC based localization procedure

    图  3   变温场下玻璃纤维环氧树脂板的实验测试系统

    Figure  3.   Experimental setup of inconsistent temperature field

    图  4   传感器PZT S0超声导波响应信号在局部加热300 s前后的幅值相位变化

    Figure  4.   The direct wave of PZT S0 gain-phase contrast after heating 300 s

    图  5   局部持续加热状态下传感器PZT S0超声导波信号幅值衰减与相位延迟曲线

    Figure  5.   Gain-phase errors contrast of PZT S0 under variable temperature field

    图  6   局部温度70.77℃时的响应时域信号

    Figure  6.   Array signal at local temperature 70.77℃

    图  7   变温场下算法修正前后的定位效果比较

    Figure  7.   Impact location comparison under variable temperature

    表  1   不同变温度环境下损伤定位结果及误差统计

    Table  1   Damage localization results and errors under variable temperature

    No. Temperature/℃ ${\hat r_1}$/cm ${\hat \theta _1}$/° E1r E1θ
    1 21.87 25.7 93 0.7 3
    2 29.54 24.3 93 0.7 1
    3 39.80 24.4 93 0.6 3
    4 49.86 25.9 93 0.9 3
    5 60.65 24.6 91 0.4 1
    6 70.77 23.8 93 1.2 3
    7 78.67 25.7 93 0.7 3
    下载: 导出CSV
  • [1]

    WANG X Y, HE J J, GAO W H, et al. Three-dimensional damage quantification of low velocity impact damage in thin composite plates using phased-array ultrasound[J]. Ultrasonics, 2020, 110: 106264.

    [2] 王志凌, 袁慎芳, 邱雷, 等. 基于压电超声相控阵方法的结构多损伤监测[J]. 振动测试与诊断, 2014(5): 796-801. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCS201405002.htm

    WANG Z L, YUAN S F, QIU L, et al. Multi-damage monitoring method based on piezoelectric ultrasonic phased array[J]. Journal of Vibration, Measurement & Diagnosis, 2014(5): 796-801. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCS201405002.htm

    [3]

    YANG H, LEE Y J, Lee S K. Impact source localization in plate utilizing multiple signal classification[C]//Proceedings of the Institution of Mechanical Engineers Part C-Journal of Mechanical Engineering Science, 2013, 227: 703-713.

    [4] 苏永振, 袁慎芳, 王瑜. 基于多重信号分类算法的复合材料冲击定位[J]. 复合材料学报, 2010(3): 105-110. https://www.cnki.com.cn/Article/CJFDTOTAL-FUHE201003016.htm

    SU Y Z, YUAN S F, WANG Y. Impact localization in composite using multiple signal classification method[J]. Acta Materiae Compositae Sinica, 2010(3): 105-110. https://www.cnki.com.cn/Article/CJFDTOTAL-FUHE201003016.htm

    [5]

    ZHONG Y T, XIANG J W, CHEN X. et al. Multiple signal classification-based impact localization in composite structures using optimized ensemble empirical mode decomposition[J]. Appl. Sci. , 2018, 8(9): 1447. DOI: 10.3390/app8091447

    [6]

    ZUO H, YANG Z, Xu C, et al. Damage identification for plate-like structures using ultrasonic guided wave based on improved MUSIC method[J]. Composite Structure, 2018, 203: 164-171. DOI: 10.1016/j.compstruct.2018.06.100

    [7]

    Radecki R, Staszewski W J, Uhl T. Impact of changing temperature on lamb wave propagation for damage detection[J]. Key Engineering Materials, 2014, 588: 140-148. https://www.scientific.net/AMM.105-107.621

    [8]

    Sikdar S, Fiborek P, Kudela P, et al. Effects of debonding on Lamb wave propagation in a bonded composite structure under variable temperature conditions[J]. Smart Materials & Structures, 2019, 28: 015021.

    [9]

    Francesco Lanza di Scalea, Salvatore Salamone. Temperature effects in ultrasonic lamb wave structural health monitoring systems[J]. Journal of the Acoustical Society of America, 2008, 161: https://doi.org/10.1121/1.2932071.

    [10]

    BAO Q, YUAN S F, WANG Y W, et al. Anisotropy compensated MUSIC algorithm based composite structure damage imaging method[J]. Composite Structure, 2019, 214: 293-303.

    [11]

    ZHANG Z H, ZHONG Y T, XIANG J W, et al. Phase correction improved multiple signal classification for impact source localization under varying temperature condition[J]. Measurement, 2020, 152: 107374.

  • 期刊类型引用(8)

    1. 李明超,闫宽,张聪,胡记伟,欧锴,陈绪兵. 用于激光软钎焊温度测量的高精度红外辐射测温装置. 红外技术. 2025(01): 108-114 . 本站查看
    2. 秦沛强,聂传杰,吝曼卿,卢永雄,张岸勤,何家冰. 磷矿巷道岩爆风险的可视化及特征增强研究. 矿业研究与开发. 2025(02): 123-131 . 百度学术
    3. 李贞,魏勇. 基于BP神经网络的红外测温补偿算法研究. 机械制造与自动化. 2023(01): 170-172+176 . 百度学术
    4. 曾飞,胡文祥,高彦鑫,宋杰杰. 基于激光扫描的输送带横向跑偏检测系统. 制造业自动化. 2023(05): 21-24 . 百度学术
    5. 杨帆,徐军,吴振生,孙明晓,金添. 基于Web端多节点红外热成像传感系统设计. 激光杂志. 2022(02): 154-157 . 百度学术
    6. 何翔. 非接触式检测装置综合实验设计. 电子技术与软件工程. 2022(07): 148-151 . 百度学术
    7. 谢彬棽. 基于双CAN总线的露天矿带式输送机速度自动化控制方法. 煤矿机械. 2022(12): 214-217 . 百度学术
    8. 武存喜. 回转窑焚烧设备退火温度模糊自适应控制技术. 工业加热. 2022(12): 27-31 . 百度学术

    其他类型引用(2)

图(7)  /  表(1)
计量
  • 文章访问数:  118
  • HTML全文浏览量:  30
  • PDF下载量:  21
  • 被引次数: 10
出版历程
  • 收稿日期:  2022-01-20
  • 修回日期:  2022-07-11
  • 刊出日期:  2022-12-19

目录

/

返回文章
返回
x 关闭 永久关闭

尊敬的专家、作者、读者:

端午节期间因系统维护,《红外技术》网站(hwjs.nvir.cn)将于2024年6月7日20:00-6月10日关闭。关闭期间,您将暂时无法访问《红外技术》网站和登录投审稿系统,给您带来不便敬请谅解!

预计6月11日正常恢复《红外技术》网站及投审稿系统的服务。您如有任何问题,可发送邮件至编辑部邮箱(irtek@china.com)与我们联系。

感谢您对本刊的支持!

《红外技术》编辑部

2024年6月6日