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

王长林，任 利，钟永腾

〈无损检测〉

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

王长林1，任 利2，钟永腾2

（1. 江苏开放大学 信息工程学院，江苏 南京 210003；
2. 温州大学 机电工程学院，浙江 温州 325035）

基于Lamb波的损伤监测方法已在复合材料结构健康监测中得到了广泛的应用。然而，复合材料结构的服役环境复杂多变，超声导波信号极易被外界因素干扰。本文研究了玻璃纤维层合板在变温场下的超声导波特性及损伤成像修正方法。首先，利用红外热成像仪研究持续加热下超声导波信号幅值和相位的变化规律，并构建因变温引起的幅相误差矩阵。其次，将误差矩阵代入多重信号分类（2D-MUSIC）算法中修正导向矢量，并建立代价函数迭代出准确损伤位置。玻璃纤维复合材料层合板的模拟损伤实验表明该方法有效提高了2D-MUSIC算法在变温场下的损伤定位分辨率和精度。

变温场；
超声导波；
2D-MUSIC算法；
损伤成像

复合材料结构健康监测（structural heath monitoring, SHM）技术是利用集成在结构中的传感、驱动元件网络，实时地获取相关响应信号，结合先进的信号处理方法以及材料结构力学建模方法，提取特征参数，识别结构的健康状态，实现结构的安全评估。Lamb波作为超声导波传播距离很远，可以实现对大面积结构的状态监测。

基于压电传感器阵列信号处理技术作为一种高分辨率的损伤成像方法，在SHM中具有广阔的应用前景。Wang等人[1]结合准静态压痕试验方法提出了一种用相控阵超声方法，实现对复合材料板低速冲击损伤定位及损伤程度预测。王志凌等人[2]提出了一种基于超声相控阵的主动结构多损伤监测方法，实现了复合材料结构中多损伤定位。Yang等人[3]结合多重信号分类算法（MUltiple SIgnal Classification, MUSIC）和弹性动力学理论提出了一种基于阵列信号处理的多冲击源定位，利用MUSIC算法估计冲击源到达方向，结合Lamb波到达时间延迟估计信号源距离。苏永振等人[4]基于MUSIC算法和小波变换技术，提出了提高复合材料结构冲击定位方法，有效地提高了方法的精度和实时性。Zhong等人[5]结合优化的集成经验模态分解（ensemble empirical mode decomposition，EEMD）和二维多重信号分类（2D-MUSIC）算法，提出了一种复合材料结构冲击定位改进方法，该方法不需分析计算具体中心频率，大大节省了算法估计时间。Zuo等人[6]根据导波传播模型建立了实验接收到的残差信号和损伤散射模型接收到的散射信号的互相关函数，并提出了一种基于模型的板状结构损伤识别的2D-MUSIC算法。上述方法在实验室环境中体现了很好的优势，然而在航空航天、机械和土木工程领域使用的板类结构在其服役期间通常经历着不同的环境条件，这些环境变化将严重影响超声导波SHM系统的可靠性。

不同研究者均证实环境和操作条件变化，特别是温度会掩盖由缺陷引起的信号变化从而限制SHM系统的性能。Radecki等人[7]研究了温度对Lamb波振幅和到达时间的影响，并提出了聚类算法判断在温度变化的情况下损伤状态。Sikdar等人[8]研究了变环境温度条件下脱层对复合材料中Lamb波传播的影响。Marzania等人[9]提出了一种基于半解析有限元模型的数值方法来预测低轮廓压电陶瓷（PZT）传感器产生和接收导波的温度效应，并用于预测－40℃～＋60℃温度范围内铝板的响应谱。Bao等人[10]针对结构各向异性对波传播的影响问题，提出了一种各向异性补偿MUSIC方法来补偿复合材料结构冲击定位中的传感器相位误差，将结构各向异性和传感器位置误差引起的传感器相位误差代入补偿算法导向矢量，并通过实验进行全方位测量。Zhang等人[11]研究了温度变化对波传播信号的影响，利用正常温度与变化的温度的对比信号，提出了一种相位校正改进的MUSIC方法。然而，上述方法均需通过大量的实验修正传感器相位误差来提高算法在温度变化下的定位精度。利用传感器本身接收的信号对阵列进行修正或微调的自适应阵列技术可以有效地解决阵列受不确定性因素的影响。而目前传感器阵列的自校正技术很少应用在结构健康监测领域。

本文针对变温场的不确定因素对超声导波阵列信号的影响问题，结合传感器阵列的自校正技术提出了复合材料板在变温场环境下的超声导波损伤成像方法。首先，借助InfraTec红外热成像仪研究了玻璃纤维层合板在变温场环境下的超声导波传播特性，并提取直达波信号幅值和相位作为误差特征建立修正模型。基于阵列超声导波信号修正模型，提出了基于修正多重信号分类（2D-MUSIC）算法的复合材料损伤成像自校正方法，有效提高了原算法在变温场下的损伤定位分辨率和精度。

1.1 损伤信号的阵列超声导波建模

如图1所示，对称布置一间距为、传感器数目为的均匀线性阵列。当中心频率为0损伤散射波到达传感器阵列时，从第个阵元输出的观测数据可以表示为：

式中：和是由于变温因素引起的与第个阵元传感器相关的幅值和相位延迟。

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

定义a(,)为近场冲击信号的阵列导向向量，可以表示为：

在实际工程应用中存在时变不确定性因素，比如运行环境中的变温场，导致该压电陶瓷传感器阵列阵元信号与理想上的模型存在误差，此时需要对原有的理论模型进行修正。如图1所示。

从第个阵元输出的观测数据修正为：

假设作为第个阵元的幅相误差，可以表示为：

得到阵列传感器误差矩阵：

通过引入误差矩阵，可以将近场冲击信号的阵列导向矢量为：

对于均匀线性阵列，观察到的输入阵列信号矢量为：

其中：

1.2 局部加热下的损伤成像方法

变温场下传感器阵列观测信号矢量的协方差矩阵为：

式中：S表示信号子空间协方差矩阵且2是白噪声，H为共轭转置。对进行特征值分解：

式中：S表示对应于最大特征值S的信号子空间；
N表示对应于最小特征值N的噪声子空间。

为了描述正交特性，可以通过以下方法计算变温场下传感器阵列的空间谱：

构建代价函数为：

给定一个预设阈值，通过迭代计算损伤位置，变温场下补偿2D-MUSIC算法的复合材料冲击损伤成像定位方法如图2所示。

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

实验系统如图3所示。变温环境实验包括集成主控结构健康监测系统、玻璃纤维环氧树脂板、红外线热成像仪及加热装置。其中红外热像仪（型号：VC HD head）的温度范围为－40℃～1200℃，可扩展至2000℃，测量精度为1%，图像规格为680×480，热灵敏度为30mK@30℃。玻璃纤维环氧树脂板的尺寸为50cm×50cm×0.2cm，板子的铺层顺序为[02/904/02]S，各层厚度为0.125mm。板子被固定在铝框架上，在传感器阵列上方25cm处布置有一个传感器模拟损伤散射源，均匀线性阵列包含有7个压电陶瓷（PZT）传感器，相邻传感器中间的间距为10mm。为研究变温度环境对传感器信号的影响，采用电炉作为局部热源对环氧树脂板进行加热，电炉的功率为300W，采用铝皮和铝箔包成一个上方开孔直径5cm圆形的圆锥体，模拟结构局部受热。

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

2.1 局部加热下超声导波的传播特性

在研究变温场下玻璃纤维环氧树脂板的实验中，从室温21.87℃将其加热到80℃左右，持续加热660s后停止加热。在加热过程中，利用热成像仪对环氧树脂板上的温度场进行实时监测并记录下局部最高温度，同时集成主控结构健康监测系统不间断地采集由压电陶瓷传感器激发的超声导波信号数据。图4所示为单个阵元损伤散射响应信号在局部加热300s前后的幅值相位变化示意，由图可知超声导波信号的直达波部分幅值略有增大，且相位发生漂移。在无损伤的状态下，超声导波信号发生了幅值和相位误差。为了研究不同变温场下的传播特性，拟合出局部持续加热状态下超声导波信号幅值衰减与相位延迟曲线，如图5所示。随着环氧树脂板上局部温度的升高，各个传感器信号的直达波相位朝一个方向漂移，相位变化和温度成近似线性的关系。

2.2 局部加热下超声导波损伤成像

实验采用压电阵列进行采集，采样频率设置为2MHz，数据采样长度为5000，预采集点数为2000，传感器阵列接收以1#阵元为触发通道，触发电压为1V。以局部温度70.77℃时的响应作为典型信号，如图6所示。在0.2～0.3ms之间出现的为损伤散射波包，为了更好地满足算法要求，通过窗函数截取此段信号作为式(9)的输入。

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

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

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

针对基于2D-MUSIC算法在变温场的成像效果比较分析，图7给出了基于标准2D-MUSIC算法及其修正算法的损伤定位方法得到的空间谱估计图。图中像素值表示空间谱的大小，其横坐标对应冲击源的角度，纵坐标对应冲击源的距离。由标准2D-MUSIC算法得到的空间谱中谱峰不明显，较难找出峰值。而修正后的算法出现了“尖锐”谱峰，修正后算法明显提高了分辨率。表1给出了修正算法在不同变温度环境下的定位结果及其误差统计，7个变温场的估计位置与实际位置较符合，误差较小，其中最大的角度误差为3°，而最大的距离误差为1.2 cm。

1）本文首先研究了玻璃纤维复合材料板在变温场环境下的超声导波传播特性，建立以幅值-相位为传播误差特征，构建了变温场下的修正模型。

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

2）本文提出了基于修正2D-MUSIC算法的损伤成像方法，该方法有效提高了原算法在变温场下的损伤定位分辨率和精度。

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

[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]., 2020, 110: 106264.

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

WANG Z L, YUAN S F, QIU L, et al.Multi-damage monitoring method based on piezoelectric ultrasonic phased array[J].,, 2014(5): 796-801.

[3] YANG H, LEE Y J, Lee S K. Impact source localization in plate utilizing multiple signal classification[C]//-, 2013, 227: 703-713.

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

SU Y Z, YUAN S F, WANG Y. Impact localization in composite using multiple signal classification method[J]., 2010(3): 105-110.

[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]., 2018, 8(9): 1447.

[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]., 2018, 203: 164-171.

[7] Radecki R, Staszewski W J, Uhl T. Impact of changing temperature on lamb wave propagation for damage detection[J]., 2014, 588: 140-148.

[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]., 2019, 28: 015021.

[9] Francesco Lanza di Scalea, Salvatore Salamone. Temperature effects in ultrasonic lamb wave structural health monitoring systems[J]., 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]., 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]., 2020, 152: 107374.

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

WANGChanglin1，RENLi2，ZHONGYongteng2

(1.,,210036,;2.,,325035,)

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 aninfrared 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.

variable temperature field, ultrasonic guided wave, 2D-MUSIC algorithm, damage imaging

V219

A

1001-8891(2022)12-1338-06

2022-01-21；

2022-07-12.

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

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

国家自然科学基金青年科学基金项目（51505339），浙江省自然科学基金青年基金项目（LQ16E050005）。

猜你喜欢 导波变温幅值 室温下7050铝合金循环变形研究装备维修技术(2022年7期)2022-07-01多尺度串联非线性能量阱的减振效能及阻尼连接方式研究振动与冲击(2022年10期)2022-05-30卷簧缺陷检测的超声导波传感器研制电子制作(2019年9期)2019-05-30氯乙烯生产中变温吸附脱水工艺的使用及改进中国盐业(2018年16期)2018-12-23基于S变换的交流电网幅值检测系统计算机仿真研究电子制作(2017年7期)2017-06-05Prevention of aspiration of gastric contents during attempt in tracheal intubation in the semi-lateral and lateral positionsWorld journal of emergency medicine(2016年4期)2016-11-23冻融处理对甘薯变温压差膨化干燥动力学的影响中国粮油学报(2016年5期)2016-01-23非共面四频激光陀螺变温零偏周期性波动中国惯性技术学报(2015年1期)2015-12-19Mn掺杂ZnSe量子点变温发光性质研究中国光学(2015年5期)2015-12-09复合管道轴对称导波改进半解析有限元建模浙江大学学报（工学版）(2015年1期)2015-03-01