基于搜索空间约束策略的MUSIC损伤定位方法技术

本发明专利技术公开一种基于搜索空间约束策略的MUSIC损伤定位方法，本发明专利技术通过单个阵元激励获得的损伤散射信号进行特征分解，直接获取信号噪声子空间，初定位结构损伤。根据结构复杂形式的先验信息，定义搜索邻域约束因子；依据初定位损伤位置及邻域约束因子确定算法进一步搜索约束空间；最后在约束空间内聚焦增强损伤散射信号，搜索精确损伤位置。传统MUSIC方法在结构健康监测中用于损伤定位时的搜索过程需要对监测结构进行全域细分、信号聚焦、特征值分解获取噪声子空间、空间谱峰比较，导致计算次数多且效率低。本发明专利技术方法可大幅减少MUSIC方法聚焦信号、噪声子空间及空间谱峰计算次数，提高MUSIC方法损伤定位效率。提高MUSIC方法损伤定位效率。提高MUSIC方法损伤定位效率。

【技术实现步骤摘要】
基于搜索空间约束策略的MUSIC损伤定位方法


[0001]本专利技术属于工程结构健康监测
，具体涉及了一种基于搜索空间约束策略的MUSIC损伤定位方法。

技术介绍

[0002]Boeing公司开发了飞机健康管理系统，通过实时故障监测和数据分析等诊断工具，提高运营效率，在飞行中也能够进行故障排除并作针对性维护。Airbus公司推出了飞行健康监测系统，可以实时收集飞机诊断信息，整理、跟踪并分析警报，从而能够快速有效地进行决策和准备最佳解决方案。因此，目前高效算法将是一个重要的需求导向。在实际工程应用中，人们逐渐发现巨大的计算量成为了阻碍MUSIC损伤定位方法实际工程化的瓶颈之一，尤其是对实时性要求较高的工程中的应用变得异常困难。
[0003]由于导波在结构中传播距离远且对小损伤敏感等优点，因此基于导波的结构健康监测方法被认为是最有工程应用潜力的在线监测方法之一。多重信号分类(Multiple signal classification,MUSIC)损伤定位方法是近年来新引入导波结构健康监测的一种阵列信号处理方法。该方法的基本思想是将任意阵列输出数据的协方差矩阵进行特征值分解，从而得到与信号分量相对应的信号子空间和与信号分量相正交的噪声子空间，然后利用这两个子空间的正交性来估计信号源的方向。
[0004]然而，现有的结构健康监测方法中，应用于压电导波密集阵列的MUSIC损伤定位方法进行损伤定位时所采用的搜索过程需要对监测结构全域细分，进行信号聚焦，实施特征值分解获取噪声子空间，再进行空间谱峰比较，进而进行损伤定位，计算次数多、计算效率低，使得对结构状态的在线实时损伤定位变得十分困难。因此，提升MUSIC损伤定位方法的效率是非常必要的。

技术实现思路

[0005]针对于上述现有技术的不足，本专利技术的目的在于提供一种基于搜索空间约束策略的MUSIC损伤定位方法，以解决现有技术中MUSIC损伤定位方法运算量大的问题；本专利技术方法提高了MUSIC损伤定位方法的运算效率。
[0006]为达到上述目的，本专利技术采用的技术方案如下：
[0007]本专利技术的一种基于搜索空间约束策略的MUSIC损伤定位方法，步骤如下：
[0008](1)采集结构健康状态传感信号：在结构健康状态下，将激励阵列中的单个阵元向结构中激发激励信号，传感阵列接收结构的响应信号；
[0009](2)采集结构损伤状态传感信号：在损伤监测过程中，将激励阵列中的单个阵元向结构中激发激励信号，传感阵列接收结构的响应信号；
[0010](3)获取单激励源下的损伤散射阵列信号：将同一激励传感通道下结构损伤状态传感信号与结构健康状态传感信号作差，得到单激励源下的损伤散射阵列信号；
[0011](4)获取噪声子空间：利用单激励源下的损伤散射阵列信号，计算阵列信号的协方
差矩阵，再对协方差矩阵进行特征值分解，得到由大到小排列的特征值，将小特征值对应的特征向量形成噪声子空间，大特征值对应的特征向量形成信号子空间；
[0012](5)设定搜索起始位置和搜索步长：将搜索起始位置设定为(r
s
,θ
s
)，即搜索位置至传感阵列中参考阵元的距离为r
s
，搜索位置相对于传感阵列的方向为θ
s
，距离和角度的搜索步长分别为r
int
和θ
int
；其中r
s
、r
int
、θ
s
、θ
int
的取值范围在监测范围内；
[0013](6)获取阵列导向矢量：计算搜索位置到达传感阵列中各阵元的距离，然后得到搜索位置至各传感阵元的距离相对于传感阵列中参考阵元的距离差，再根据导波传播速度得到搜索位置至各传感阵元的时间相对于参考阵元的时间差，定义为阵列时延，根据获得的阵列时延计算出阵列导向矢量；
[0014](7)计算基于单激励源下的空间谱，对空间谱成像后得到损伤初定位位置；
[0015](8)根据结构复杂形式的先验信息定义搜索邻域约束因子；
[0016](9)确定搜索约束空间：当损伤初定位位置与复杂结构形式的间距在R以内，则判定损伤初定位位置在复杂结构形式上；否则判定损伤初定位位置在非复杂结构形式上，进一步结合具体结构形式自适应设定搜索约束空间；
[0017](10)在约束空间内对多个单激励源下的损伤散射阵列信号进行聚焦增强；计算各激励阵元到达搜索位置处的时延，按照时延将各阵元激励时的散射信号前移或后移，再将施加时延后的损伤散射信号进行叠加，得到叠加后的损伤散射信号；
[0018](11)计算搜索约束空间谱并成像，成像结果峰值点为损伤位置。
[0019]进一步地，在步骤(1)中包含激励阵列E和传感阵列T，每个阵列均由2K+1个阵元组成，阵元等间距排列，激励阵列E中的每个阵元用E
g
表示，传感阵列T中的每个阵元用T
g
表示，g表示阵列的阵元编号，g＝
‑
K,
‑
(K
‑
1),
…
,0,
…
,K
‑
1,K。
[0020]进一步地，所述步骤(3)中基于单激励源E0得到的损伤散射阵列信号为：
[0021][0022]式中，H(t)表示结构健康状态下基于单激励源E0时的传感阵列响应信号，D(t)表示结构损伤状态下基于单激励源E0时的传感阵列响应信号；表示传感阵列T基于单激励源E0得到的损伤散射阵列信号。
[0023]进一步地，所述步骤(4)中，损伤散射阵列信号的协方差矩阵C：
[0024][0025]式中，为的Hermitian转置，Y为信号采样长度；
[0026]对协方差矩阵C进行特征分解：
[0027][0028]式中，Z
S
、Z
η
分别为信号子空间与噪声子空间，∑
S
、Σ
η
分别为信号子空间对应的大特征值与噪声子空间对应的小特征值。
[0029]进一步地，所述步骤(6)中阵列导向矢量A(r,θ)中的元素表示为
ω0为信号传播的中心频率，j为虚数单位，b为相邻阵元的间距，v为信号传播速度，r和θ分别为搜索位置相对于传感阵列T的距离和角度，τ
g
表示阵列T中各阵元相对参考阵元T0的信号到达时间延迟。
[0030]进一步地，所述步骤(7)中，利用步骤(4)中的噪声子空间Z
η
以及步骤(6)中的阵列导向矢量A(r,θ)，计算搜索位置(r,θ)对应的空间谱J
MUSIC
，空间谱表达式为：
[0031][0032]式中，A
H
(r,θ)为A(r,θ)的Hermitian转置，为Z
η
的Hermitian转置；
[0033]对空间谱进行成像，在空间谱图中存在一个本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于搜索空间约束策略的MUSIC损伤定位方法，其特征在于，步骤如下：(1)采集结构健康状态传感信号：在结构健康状态下，将激励阵列中的单个阵元向结构中激发激励信号，传感阵列接收结构的响应信号；(2)采集结构损伤状态传感信号：在损伤监测过程中，将激励阵列中的单个阵元向结构中激发激励信号，传感阵列接收结构的响应信号；(3)获取单激励源下的损伤散射阵列信号：将同一激励传感通道下结构损伤状态传感信号与结构健康状态传感信号作差，得到单激励源下的损伤散射阵列信号；(4)获取噪声子空间：利用单激励源下的损伤散射阵列信号，计算阵列信号的协方差矩阵，再对协方差矩阵进行特征值分解，得到由大到小排列的特征值，将小特征值对应的特征向量形成噪声子空间，大特征值对应的特征向量形成信号子空间；(5)设定搜索起始位置和搜索步长：将搜索起始位置设定为(r
s
,θ
s
)，即搜索位置至传感阵列中参考阵元的距离为r
s
，搜索位置相对于传感阵列的方向为θ
s
，距离和角度的搜索步长分别为r
int
和θ
int
；(6)获取阵列导向矢量：计算搜索位置到达传感阵列中各阵元的距离，然后得到搜索位置至各传感阵元的距离相对于传感阵列中参考阵元的距离差，再根据导波传播速度得到搜索位置至各传感阵元的时间相对于参考阵元的时间差，定义为阵列时延，根据获得的阵列时延计算出阵列导向矢量；(7)计算基于单激励源下的空间谱，对空间谱成像后得到损伤初定位位置；(8)根据结构复杂形式的先验信息定义搜索邻域约束因子；(9)确定搜索约束空间：当损伤初定位位置与复杂结构形式的间距在R以内，则判定损伤初定位位置在复杂结构形式上；否则判定损伤初定位位置在非复杂结构形式上，进一步结合具体结构形式自适应设定搜索约束空间；(10)在约束空间内对多个单激励源下的损伤散射阵列信号进行聚焦增强；计算各激励阵元到达搜索位置处的时延，按照时延将各阵元激励时的散射信号前移或后移，再将施加时延后的损伤散射信号进行叠加，得到叠加后的损伤散射信号；(11)计算搜索约束空间谱并成像，成像结果峰值点为损伤位置。2.根据权利要求1所述的基于搜索空间约束策略的MUSIC损伤定位方法，其特征在于，在步骤(1)中包含激励阵列E和传感阵列T，每个阵列均由2K+1个阵元组成，阵元等间距排列，激励阵列E中的每个阵元用E
g
表示，传感阵列T中的每个阵元用T
g
表示，g表示阵列的阵元编号，g＝
‑
K,
‑
(K
‑
1),
…
,0,
…
,K
‑
1,K。3...

【专利技术属性】
技术研发人员：郑飞，袁慎芳，徐秋慧，任元强，陈健，邱雷，徐钰森，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：