基于激励波束成型和加权图像融合的MUSIC腐蚀监测方法技术

本发明专利技术公开了基于激励波束成型和加权图像融合的MUSIC腐蚀监测方法。该方法针对损伤散射信号弱并影响MUSIC精度，首先将激励波束成型和双阵列引入到MUSIC方法中，增强腐蚀损伤的散射信号和提高其信噪比。其次将双阵列轮流作为激励或传感阵列，通过加权融合双阵列的腐蚀损伤成像，实现了一维线阵盲区内的腐蚀损伤监测。最后提出了基于阵列信号协方差矩阵特征值的腐蚀损伤因子，用以判别腐蚀损伤深度。本发明专利技术提高了MUSIC腐蚀监测的定位精度，扩大了传统一维线阵的监测区域，实现了腐蚀损伤深度的评估，在真实航空结构的腐蚀损伤监测方法具有广泛的应用前景。 1

【技术实现步骤摘要】
基于激励波束成型和加权图像融合的MUSIC腐蚀监测方法
本专利技术属于工程结构健康监测
，特别涉及了一种MUSIC腐蚀损伤监测方法。
技术介绍
目前全球有超过14000架老龄化飞机，而老龄化飞机通常需要更加频繁的检测和维护。作为一种频繁发生的损伤，很多老龄化飞机正遭受着结构腐蚀。根据美国联邦航空管理局统计，老龄化飞机有关腐蚀的维护花费占总维护费用的30％。因此，腐蚀损伤监测成为当前结构健康监测中一个很重要的课题。由于Lamb波具有长距离的传播能力且对小损伤敏感，基于Lamb波的结构健康监测方法具有很好的应用前景。由于传感器阵列易于在结构上布置且具有方向扫描的功能，阵列信号处理方法逐渐被引入结构健康监测区域。近些年，MUSIC(多重信号分类)算法是其中一种代表性的算法。MUSIC算法属于子空间类算法，其算法的基本思想是将任意阵列输出数据的协方差矩阵进行特征值分解，从而得到与信号分量相对应的信号子空间和与信号分量相正交的噪声子空间，然后利用这两个子空间的正交性来估计信号的参数。在MUSIC腐蚀损伤监测中，腐蚀损伤并不能产生Lamb波，因而需要采用主动Lamb波方法并利用腐蚀损伤的散射信号来实现定位和程度评估。然而腐蚀损伤的散射信号很微弱，信噪比低，从而导致MUSIC方法的定位精度低。其次，通常的MUSIC损伤定位方法中采用一维均匀线阵来接受Lamb波，但一维线阵在靠近0°和180°角度范围内存在监测盲区问题。若腐蚀损伤发生在盲区内，即[0°30°]和[150°180°]，则对其不能定位。此外，目前MUSIC算法还不能评估腐蚀损伤的程度。
技术实现思路
为了解决上述
技术介绍
提出的技术问题，本专利技术旨在提供基于激励波束成型和加权图像融合的MUSIC腐蚀监测方法，增强腐蚀损伤的散射信号并提高其信噪比，消除传统一维均匀线阵的监测盲区，并实现腐蚀损伤的定位和深度评估为了实现上述技术目的，本专利技术的技术方案为：基于激励波束成型和加权图像融合的MUSIC腐蚀监测方法，包括以下步骤：(1)在结构处于健康状态时，轮流驱动阵列A和阵列S激发Lamb波，并采集对应传感器阵列的响应信号，记为基准信号；(2)在腐蚀损伤监测过程中，轮流驱动阵列A和阵列S激发Lamb波，并采集对应传感器阵列的响应信号，记为监测信号；将监测信号减去基准信号，记为腐蚀损伤的散射信号；(3)基于MUSIC算法和腐蚀损伤的散射信号，得到腐蚀成像结果，搜索腐蚀成像结果的峰值点，即为腐蚀损伤的初估位置；(4)根据此初估位置，计算激励阵元相对腐蚀损伤位置的时延，并根据该时延前移或后移腐蚀损伤的散射信号，叠加这些腐蚀损伤的散射信号获得增强的腐蚀损伤的散射信号；(5)将增强的腐蚀损伤的散射信号再次代入MUSIC算法，得到对应的腐蚀成像结果和散射信号协方差矩阵的大特征值；(6)对阵列A和阵列S分别作为激励源阵列时对应的腐蚀成像结果设置相应的权值，根据权值融合两次腐蚀成像结果，搜索融合成像的峰值点，即为腐蚀损伤位置；(7)根据步骤(6)得到的权值叠加散射信号协方差矩阵的大特征值，进而计算腐蚀因子，利用腐蚀因子评估腐蚀损伤深度。进一步地，所述阵列A和阵列S均为一维均匀线阵，两个阵列均包含2N+1个阵元，阵列A中的各个阵元用Ai表示，阵列S中的各个阵元用Si表示，i＝-N,-(N-1),…,0,…,N-1,N。进一步地，在步骤(2)中，驱动阵列A激发Lamb波，即依次将阵列A中每个阵元作为激励，得到对应的腐蚀损伤的散射信号：上式中，为阵列A中阵元Ai作为激励时，阵列S中所有阵元得到的散射信号组成的向量；N(t)为阵列S中各阵元采集通道的背景噪声组成的向量；xa(t)为阵列S中参考阵元获取的腐蚀损伤的散射信号；A(r,θ)为导向矢量，A(r,θ)中的元素aq(r,θ)＝exp(-jω0τq)，q＝N,-(N-1),...,0,...,N-1,N，ω0为传播信号的中心频率，j为虚数单位，d为阵列中相邻阵元的间距，c为信号传播速度，r和θ分别为损伤相对于阵列S的距离和角度；t表示信号的时刻；驱动阵列S激发Lamb波的过程与阵列A相同。进一步地，在步骤(3)中，首先计算的协方差矩阵上式中，上标H表示Hermitian转置，K为采集信号长度；然后对进行特征值分解：上式中，US、UN分别为信号子空间、噪声子空间，∑S、∑N分别为信号子空间对应的大特征值、噪声子空间对应的小特征值；基于MUSIC算法的空间谱估计：在空间谱图中存在一个明显的波峰，即为腐蚀损伤的初估位置。进一步地，在步骤(4)中，计算激励源阵列中各阵元与腐蚀损伤的初估位置的距离：p＝N,-(N-1),...,0,...,N-1,N上式中，l为激励源阵列与传感器阵列之间的间距，r1为腐蚀损伤初估位置的距离信息，θ1为腐蚀损伤初估位置的方向信息；计算激励源阵列中各阵元至腐蚀损伤的时延：上式中，为激励源阵列中参考阵元与腐蚀损伤的初估位置的距离；叠加各腐蚀损伤的散射信号获得增强的腐蚀损伤的散射信号：进一步地，在步骤(6)中，若其中一个腐蚀成像结果定位的腐蚀损伤位于其传感器阵列的监测盲区，则设置该腐蚀成像结果的权值为0，另一腐蚀成像结果的权值为1；若两个腐蚀成像结果定位的腐蚀损伤均处于各自传感器阵列的可监测区域内，则设置两个腐蚀成像结果的权值均为0.5。进一步地，在步骤(7)中，按照下式计算腐蚀因子：上式中，λmax为根据权值叠加的散射信号协方差矩阵的大特征值，uA和uS分别是激励信号和腐蚀散射信号的峰值。进一步地，在步骤(7)中，根据下式计算腐蚀损伤深度：D＝1.29×CI+0.115(cm)。采用上述技术方案带来的有益效果：本专利技术有效提高了基于MUSIC算法的腐蚀损伤定位精度；能够监测一维线阵监测盲区内的腐蚀损伤，扩大了监测区域；能够评估腐蚀损伤的深度，可有效地应用于航空结构的真实腐蚀损伤监测中。附图说明图1为本专利技术的方法流程图；图2为实施例中铝板及双阵列布置示意图；图3为实施例中单激励源下MUSIC腐蚀损伤定位结果图，其中(a)为单激励源下腐蚀损伤散射阵列信号图，(b)为空间谱图；图4为实施例中基于激励波束成型的MUSIC腐蚀损伤定位结果，其中(a)为聚焦的腐蚀损伤散射阵列信号图；(b)为空间谱图；图5为实施例中加权融合的空间谱图，其中(a)为阵列S作为激励源阵列时的空间谱图；(b)为加权融合空间谱；图6为实施例中腐蚀因子与腐蚀深度的关系图。具体实施方式以下将结合附图，对本专利技术的技术方案进行详细说明。本专利技术提出了基于激励波束成型和加权图像融合的MUSIC腐蚀监测方法，如图1所示，过程如下。步骤1：在结构处于健康状态时，轮流驱动阵列A和阵列S激发Lamb波，并采集对应传感器阵列的响应信号，记为基准信号；步骤2：在腐蚀损伤监测过程中，轮流驱动阵列A和阵列S激发Lamb波，并采集对应传感器阵列的响应信号，记为监测信号；将监测信号减去基准信号，记为腐蚀损伤的散射信号；步骤3：基于MUSIC算法和腐蚀损伤的散射信号，得到腐蚀成像结果，搜索腐蚀成像结果的峰值点，即为腐蚀损伤的初估位置；步骤4：根据此初估位置，计算激励阵元相对腐蚀损伤位置的时延，并根据该时延前移或后移腐蚀损伤的散射信号，叠加这些腐蚀损伤的散射信号获得增强的腐蚀损伤的散射本文档来自技高网...

【技术保护点】
1.基于激励波束成型和加权图像融合的MUSIC腐蚀监测方法，其特征在于，包括以下步

【技术特征摘要】
1.基于激励波束成型和加权图像融合的MUSIC腐蚀监测方法，其特征在于，包括以下步骤：(1)在结构处于健康状态时，轮流驱动阵列A和阵列S激发Lamb波，并采集对应传感器阵列的响应信号，记为基准信号；(2)在腐蚀损伤监测过程中，轮流驱动阵列A和阵列S激发Lamb波，并采集对应传感器阵列的响应信号，记为监测信号；将监测信号减去基准信号，记为腐蚀损伤的散射信号；(3)基于MUSIC算法和腐蚀损伤的散射信号，得到腐蚀成像结果，搜索腐蚀成像结果的峰值点，即为腐蚀损伤的初估位置；(4)根据此初估位置，计算激励阵元相对腐蚀损伤位置的时延，并根据该时延前移或后移腐蚀损伤的散射信号，叠加这些腐蚀损伤的散射信号获得增强的腐蚀损伤的散射信号；(5)将增强的腐蚀损伤的散射信号再次代入MUSIC算法，得到对应的腐蚀成像结果和散射信号协方差矩阵的大特征值；(6)对阵列A和阵列S分别作为激励源阵列时对应的腐蚀成像结果设置相应的权值，根据权值融合两次腐蚀成像结果，搜索融合成像的峰值点，即为腐蚀损伤位置；(7)根据步骤(6)得到的权值叠加散射信号协方差矩阵的大特征值，进而计算腐蚀因子，利用腐蚀因子评估腐蚀损伤深度。2.根据权利要求1所述MUSIC腐蚀监测方法，其特征在于，所述阵列A和阵列S均为一维均匀线阵，两个阵列均包含2N+1个阵元，阵列A中的各个阵元用Ai表示，阵列S中的各个阵元用Si表示，i＝-N,-(N-1),…,0,…,N-1,N。3.根据权利要求2所述MUSIC腐蚀监测方法，其特征在于，在步骤(2)中，驱动阵列A激发Lamb波，即依次将阵列A中每个阵元作为激励，得到对应的腐蚀损伤的散射信号：上式中，为阵列A中阵元Ai作为激励时，阵列S中所有阵元得到的散射信号组成的向量；N(t)为阵列S中各阵元采集通道的背景噪声组成的向量；xa(t)为阵列S中参考阵元获取的腐蚀损伤的散射信号；A(r,θ)为导向矢量，A(r,θ)中的元素aq(r,θ)＝exp(-jω0τq)，q＝N,-(...

【专利技术属性】
技术研发人员：袁慎芳，鲍峤，邱雷，郭方宇，任元强，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏,32