专利名称：基于多重信号分类算法的冲击载荷定位方法
技术领域：
本发明基于阵列信号处理技术中的MUSIC算法的结构冲击载荷位置识别方法，属 工程结构健康监测方法。
背景技术：
复合材料以其比强度高、比刚度大、抗疲劳性能好及材料性能可设计等一系列优 点，在航空、航天、汽车等工程领域得到了日益广泛的应用，尤其是在军用、民用飞机上已开 始越来越多地使用先进复合材料结构。然而层合板复合材料在服役过程中不可避免的要承 受各种能量物体的冲击，比如，冰雹的撞击、飞行器翼面与空中飞鸟的碰撞、飞行器受到枪 击及工具经常在维护过程中掉落在飞行器表面等。这些冲击极易造成复合材料结构的内 部分层、基体开裂和纤维断裂等损伤。这些内部损伤将使层合结构的力学性能严重退化， 强度可削弱35% 40%，导致承载能力大大降低，对结构的整体破坏和失效形成潜在的威 胁。而且这些损伤多发生在材料内部不易从表面发现，留下严重的隐患，使得具有损伤的复 合材料结构具有突发性和灾难性失效的潜在能力，并且它在损伤、失效等方面的表现却是 机理复杂，现象多样，判别困难。因此很有必要对复合材料结构进行全寿命的监测，以确保 结构的稳定性和安全性。目前已有许多传统的已被广泛应用的无损检测技术，例如敲击、超 声、X射线、电涡流射线、电位测量以及热应力场等方法。但是这些检测方法一般设备复杂、 耗时耗力，需要对损伤的位置有初步的了解，使用不方便，局限性大，不易做到服役环境下 的实时在线监测，不适合未来大型航空、航天飞行器结构的健康监测与诊断。对于在服役过 程中的飞行器，重要的是如何来实时监测冲击试件以保护飞行器的安全运行。而利用基于 压电智能结构的健康监测技术可以实现低速冲击的实时、在线监测，评估结构的损伤程度， 从而保证结构的稳定性和安全性以及降低结构的维护成本。目前冲击载荷识别的研究主要集中在实验室针对小试件进行，多是在模拟冲击载 荷的情况下进行的实验研究，而且实验对象也多以结构简单的梁或板结构为主，很少在复 合材料结构上进行真实的冲击实验。尽管已有各种冲击载荷识别方法，但是目前还没有形 成一种能够高精度、实时在线识别冲击载荷的方法。对于复合材料结构的冲击定位问题，原 理上可以利用声发射技术进行定位，但由于复合材料的各向异性和应力波在传播过程中的 反射、折射、散射及噪声的影响，很难准确确定波速及波达时间，因此单纯采用常规的声发 射技术很难对冲击源进行精确定位。而基于优化技术的识别方法在反演过程中需要不断地 进行结果迭代，求解时间过长。因此，对于范围稍大的监测区域来说，搜寻最优位置所需的 时间更长，单纯运用最优化方法进行实时在线监测已不大可能。并且基于优化技术的识别 方法通常直接将时域内许多离散点的冲击力时间历程本身当作未知参数进行求解。由于要 求解的未知参数过多，计算任务重，也很难满足在线监测的实时性要求。基于人工智能的荷 载识别方法可以解决荷载识别问题，但用于荷载识别的模型需要以实际加载方式或解析解 求出响应样本的方法来训练，只能用于已知结构而且对结构的训练工作量较大，训练后的 结构在实际荷载识别时亦存在不稳定现象，因此，应用于冲击监测这一实际问题仍有一定的难度。阵列信号处理技术作为信号处理领域一个重要的分支，近几十年来已在雷达、声 纳、通信、地震勘探、射电天文及医学诊断等多种国民经济和军事领域广泛应用。其与传统 的单个或少量定向传感器相比，阵列信号处理具有灵活的波束控制、高的信号增益、极强的 抗干扰能力及高的空间超分辨能力等优点。MUSIC是Schmidt和Bienvenu及Kopp于1979 年提出，由Schmdit在1986年重新发表。MUSIC算法的提出促进了特征结构分类算法的兴 起和发展，该算法已成为空间谱估计理论体系中的标志性算法。由于它能对多个空间信号 进行识别，故这种方法称为多重信号分类法。由于MUSIC方法发展成熟，分辨率高，对阵列 结构要求低等特点，以及超高速超大运算能力信号处理芯片的出现使其在阵列信号处理中 的应用非常广泛。因此，将阵列信号处理技术中的MUSIC算法引入到结构健康监测领域中， 研究冲击源的定位问题。
发明内容
技术问题本发明的技术问题是提供一种利用阵列信号处理技术中的MUSIC算法 和Lamb波频散特性的结构冲击载荷位置识别方法。本方法利用MUSIC方法实现冲击源波达 方向的估计，根据Lamb波的频散特性实现冲击源距离的估计，可以对冲击实时、在线监测， 准确、快速地识别出冲击载荷位置。技术问题本发明为实现上述目的，采用如下技术方案本发明基于多重信号分类算法的结构冲击载荷定位方法，其特征在于包括如下步 骤步骤一冲击响应信号的采集(1)在所述结构上布置由多个相同型号的压电元件组成的压电阵列；(2)将2块4通道的数采卡0通道并联在一起，保证两块数采卡同步触发采集；(3)设定采集参数，利用电荷放大器和数据采集卡采集压电阵列输出的冲击响应
信号；步骤二 基于MUSIC算法进行冲击源的波达方向估计首先用快速傅里叶变换对所述冲击响应信号进行频谱分析，然后根据选定的方向 估计中心频率利用小波变换从所述冲击响应信号中提取出中心频率下的窄带信号，最后取 窄带信号的幅值最大的信号运用多重信号分类(MUSIC)算法的进行D0A估计；步骤三基于Lamb波传播特性进行冲击源的距离估计首先用小波变换从任一压电元件接收到的冲击响应信号中提取出距离估计中心 频率下的Lamb波信号，获得同一频率的对称&模式的Lamb波和反对称&模式的Lamb波， 用最大峰值法求出对称&和反对称A0这两种模式的Lamb波到所述达压电元件的时间差， 最后由测得的这两种模式的Lamb波的传播速度求出冲击源到达压电元件的距离，将各压 电元件估计出的距离值转化为相对于压电阵列参考点的距离，取所有压电元件测得的距离 的平均值作为距离估计值。优选地，步骤一中，所述压电阵列中压电元件间的距离小于冲击源信号波长。有益效果本发明与现有技术相比，本方法简便可靠，可对冲击进行实时在线监 测，实时性强、定位精度高，能实际应用于工程结构中冲击载荷的监测及定位，确保工程结构的运行安全，防止事故的发生。
图1是一均勻线阵在窄带远场条件下是我信号模型的示意图。它是由M个阵元等 间距排列成一条直线的均勻线阵，以左边第一个阵元为参考点，以垂直于阵列的方向为法 线方向。N个远场的窄带信号以e角入射到此均勻线阵上。图2是铝板上的一组Lamb波频散曲线。纵坐标为群速度，横坐标为板厚与频率的 乘积(频厚积)。图3是压电阵列测得的复合材料结构上的冲击响应信号。图4是截取出的0. 4-0. 6ms时间段的小波信号波形。图5是用MUSIC方法估计出的方位谱。图6是用小波变换提取出的Lamb波，小波变换的中心频率为25kHz。图7是本发明方法流程图。
具体实施例方式下面结合附图对发明的技术方案进行详细说明如图7所示，本发明的方法包括如下步骤(1)冲击响应信号的采集包括如下步骤(1. 1)在结构上布置一个由7个相同型号的压电元件组成的均勻线阵；(1. 2)将2块4通道的数采卡0通道并联在一起，保证两块数采卡同步触发采集。(1. 3)设定采集参数，如采样频率、触发预置、预采点数、采集点数等。利用电荷放 大器和数据采集卡采集冲击响应信号。步骤二 基于MUSIC算法进行冲击源的波达方向估计。MUSIC方法在窄带远场条件下可实现信号源的D0A估计。但冲击响应信号频率成 分十分复杂，直接测得的信号不符合窄带信号模型要求，因此本发明实现D0A估计的基本 思路就是首先利用小波变换对由压电传感器阵列单元直接测得的冲击响应信号进行预处 理，提取出利用MUSIC算法进行D0A估计所需要的窄带信号。具体做法就是首先用快速傅 里叶变换对由压电传感器阵列直接测得的冲击响应信号进行频谱分析，然后选定某一中心 频率，利用小波变换从压电传感器阵列直接测得的冲击响应信号中提取出此中心频率下的 窄带信号，最后取此提取出的响应波信号的某一时间段内的幅值较大的信号运用多重信号 分类(MUSIC)算法的进行D0A估计。下面结合MUSIC方法的基本原理对本发明实现D0A估 计的方法进行说明。MUSIC算法的基本思想就是则是将任意阵列输出数据的协方差矩阵进行特征分 解，从而得到与信号分量相对应的信号子空间和与信号分量相正交的噪声子空间，然后利 用这两个子空间的正交性来估计信号的参数(入射方向、极化信息及信号强度等)。在窄带远场条件下，N个信号源以入射角e入射到由M个压电阵元组成的均勻线 阵上，如图1所示，D0A估计的数学模型为x(t) = Bs(t)+Noise(t)(1)
式中，x(t)为阵列的MX1维数据矢量，NoiSe(t)为阵列的MX1维噪声数据矢量， s(t)为空间信号的NX 1维矢量，M和N分别为阵元数目和信号源数目，B为空间阵列的MXN 维的流型矩阵(导向矢量阵)，且BztMejMe^HMeN)]⑵其中，导向矢量b( 0 i) = [1 exp(-j 0T2i)...exp(-j 0TMi)]T，i = 1,2, ...，N (3)其中，《0 = 2Jif，f为入射信号频率。下面根据窄带远场信号的数学模型对MUSIC算法实现D0A估计的过程进行分析。 由窄带远场信号的数学模型式(1)可知阵列输出数据的协方差矩阵为R = E[xxh] = BE[ssH]BH+BE[sNoiseH]+E[NoisesH]B+E[NoiseNoiseH] (4)由于假设信号与噪声是不相关的，且噪声为平稳的加性高斯白噪声，因此上式中 的二、三项为零，上式可化简为R = E[xxh] = be[ssh]bH+E[NoiseNoiseH] = brsbh+rn (5)式中，Rs，&分别为信号协方差矩阵和噪声协方差矩阵，设噪声功率为o2，则有下 式R = BRsBh+Rn = BRsBh+ o 2I(6)由于信号与噪声相互独立，数据协方差矩阵可分解为与信号、噪声相关的两部分， 其中Rs是信号的协方差矩阵，BRsBh是信号部分。对R进行特征值分解有 式中，特征值满足关系入2彡…入XN+1= =An= o2。信号协方差矩阵R特征值分解得到的特征值即为信号功率，在信号功率大于噪声 功率的条件下大特征值就对应信号子空间，小特征值对应噪声子空间。因此Us是由大特征 值对应的特征矢量张成的信号子空间，而UN是由小特征值对应的特征矢量张成的子空间即 噪声子空间。M个特征值对应的特征向量分别为V2，…，VN，VN+1，…，VM，其中前N个对应大 特征值，后M-N个对应小特征值。对于M-N个小特征值对应的特征向量，有(R-A jDVj = 0 i = N+l, ...，M即 因为B满秩，Rs非奇异，则有BHVi = 0 i = N+l, ...，M上式表明，矩阵B的各个列向量与M-N个最小特征值对应的特征向量正交，也即是 矩阵B的各个列向量与噪声子空间正交，故有(11)式。由于Rs非奇异，则阵列方向矩阵B， 与阵列输出向量的协方差矩阵的N个大特征值对应的特征向量所张成的信号子空间Us相 同，也就是信号子空间Us和噪声子空间UN正交。VHNb{9) = 0(11)MUSIC算法正是基于上述性质提出的，但在实际应用中，理想的协方差矩阵是未知 的，并且无法得到，所谓的协方差矩阵R是利用有限个数据向量X⑴，X⑵，…，X(L)估计 得到的，数据协方差矩阵的最大似然函数估计为 式中，X= [X⑴，X⑵，…，X(L)]，X为数据向量矩阵，L为采样数据的段数。对i进行特征分解可以计算得到噪声子空间Uy。在实际情况下，存在的噪声与信 号源并不是完全不相关，信号方向与噪声方向并不完全正交，即b( e)与并不能完全正
交，也就是式(11)并不严格成立。因此，实际上求doa是以最小优化搜索实现的，即 式中，arg min表示函数值最小时所对应的角度。于是可以定义MUSIC算法的空
间谱函数 根据正交性，上式会在信号D0A值附近出现峰值，与“谱峰”对应的9，即为D0A的 估计。由上述分析可知，利用MUSIC算法实现冲击源的D0A估计的步骤可归纳如下(1)根据公式(12)得到阵列接收数据的协方差矩阵估计R ；(2)对益进行特征值分解；(3)确定信号子空间仏和噪声子空间(4)根据信号参数范围由式(14)进行谱峰搜索；(5)找出极大值点对应的角度就是信号源的D0A。下面对利用MUSIC算法实现冲击源D0A估计的所需的运算量进行分析，主要包括 三部分首先需要估计阵列协方差矩阵，其运算量为0(M2L)；其次对阵列协方差矩阵作特征 值分解，其运算量为0(M3)；最后就是“谱峰”搜索需要的运算量。对于本发明的单冲击源的 情况，阵元数目为7，由运算量分析可知，由于阵元数目较少，MUSIC算法估计D0A所需的运 算量也就很小，因此该方法具有较高的实时性。步骤三基于Lamb波传播特性进行冲击源的距离估计。在无限均勻各向同性弹性介质中，只存在两种波——横波和纵波，二者分别以各 自的特征速度传播而无波形耦合，而在板中则不然，在板的某一点上激励声波，由于声波传 播到板的上、下界面时会发生波形转换，经过在板内一段时间的传播之后，因叠加而产生波 包，即所谓的板中Lamb波模态。Lamb波在板中传播时，存在不同的模式，各种模式的叠加效 果即为Lamb波。本质上，Lamb波是在具有两个平行表面的结构中由横波和纵波相互耦合 而成的一种弹性波。Lamb技术是目前研究较多的一种板类结构损伤监测技术。板波理论认为，对于工程上大量使用的板状结构，机械波在板中主要激励起扩展波(最低阶对称波)和 弯曲波(最低阶反对称波)两种模式。薄板结构在承受法向冲击产生的板波声发射信号从 本质上来说是一种Lamb波。Lamb波在板中传播时，板中质点产生振动，其振动方式十分复杂，随各种参数(频 率、板厚等)的变化而改变。根据薄板两表面质点的振动相位关系，可以把Lamb波的模式 分为对称波和反对称波两种，对称波又分为S0, S”…，3 等多个模式，反对称波又分为k0’ 、，…，Am等多个模式。Lamb波传播特性极为复杂，主要体现在其频散特性和多模现象。所谓的频散特性 就是Lamb波的速度随着频率的变化而发生改变。而多模现象就是同一个频率值对应着不 止一种Lamb波模式。根据板波理论，对于薄板结构，冲击产生的Lamb波信号，主要是&和 A0两种模式的Lamb波，而且在较低的频厚积范围内同一频率的&模式Lamb波比k0模式 Lamb波传播速度快，一组铝板中传播的Lamb波群速度频散曲线如图2所示。因此，可根据 同一频率的这两种模式的Lamb波的传播速度不同求出这两种模式的Lamb波到达同一阵列 单元的时间差，再根据测得的这两种模式的Lamb波的传播速度即可实现距离估计。薄板结构因受到冲击而产生的Lamb波信号极为复杂，含有较多的模式和丰富的 频率成分。这就使直接利用冲击响应信号实现距离估计变得非常困难。因此在估计冲击源 的距离时，也同样用小波变换提取出第1(1 = 0，…，M-1)个压电阵元的某一中心频率下的 Lamb波响应信号，可以获得同一频率的对称&模式的Lamb波和反对称&模式的Lamb波， 然后根据同一频率的对称&模式和反对称A0模式这两种模式的Lamb波的传播速度不同， 可用最大峰值法求出&和&这两种模式的Lamb波到达第此阵元的时间差At，最后由测 得的这两种模式的Lamb波的传播速度cs和cA就求出冲击源到达此阵元的距离屯，即 为了提高距离估计的精度，将各阵元估计出的距离值屯转化为相对于阵列参考点
的距离巧，取所有阵元测得的距离的平均值作为距离估计值R。 实施例试件为T-30碳纤维复合材料层合板，尺寸为610 X 610 X 2mm，四边通过螺钉紧固 在夹具上。碳纤维铺层的材料为T300/QY8911，铺层顺序为[45/0/-45/90/0/45/0/-45/0] so采用由7个相同的压电片组成的均勻线阵，相邻两个压电片中心间的距离为1.1cm。阵 元的标号均是从左到右依次为0 6，均以0号压电片的中心为坐标原点，以压电阵列的轴 线为实轴，以垂直于阵列的方向为法线方向，如图3所示。传感元件为PZT-5型压电片，直 径为8mm，厚度为0. 48mm。冲击工具为弹簧机械冲击杆，共有两档冲击能量，冲击头顶部为 半球形，直径为10mm，如图3所示。由于碳纤维复合材料层合板为各向异性，各方向的波速 不同，所以在进行D0A估计时波速取为0°、45°、90°及-45°方向上的测量值的均值；在 进行距离估计时，&和&两种模式的波速分别取为测得的各阵元接收到的同一中心频率的S0模式Lamb波的波速均值和k0模式Lamb波的波速均值。本实施例包括如下步骤步骤一冲击响应信号。当冲击锤头作用在试件上的(29.0°，47. 6cm)位置时，每个压电元件上的冲击响 应信号经过各自的电荷放大器，输入到数据采集卡上，当信号幅值超过设定的触发预置时， 数据采集卡开始对冲击响应信号进行采集，并显示、存储在计算机中，0 6号压电片采集 到的传感信号如图3所示，从上到下依次为0 6号压电片的信号波形。步骤二 首先利用快速傅里叶变换对步骤一中获得的冲击响应信号进行频谱分 析，然后用小波变换对由压电传感阵列直接测得的冲击响应信号做预处理，小波变换的中 心频率为35kHz。然后从提取出的小波信号中截取信号幅值较大的0. 4 0. 6ms时间段内 的数据进行D0A估计，截取获得的小波信号如图4所示。最后对截取得到的信号按照MUSIC 算法的步骤进行D0A估计(1)根据公式(12)中得到阵列接收数据的协方差矩阵估计R ；(2)对叙进行特征值分解；(3)确定信号子空间亡s和噪声子空间(4)根据信号参数范围由式(14)进行谱峰搜索；(5)找出极大值点对应的角度就是信号源的D0A。按照上述步骤得到冲击源的方位谱如图5所示，最大谱峰对应的角度为28. 7°， 角度估计误差为0. 3°，由误差可以看出D0A估计的精度非常高。步骤三冲击源距离的估计。按照技术路线步骤三中所述的冲击源距离估计方法 进行距离估计。以对0号压电片接收到的信号分析为例，首先用小波变换提取出中心频率 为25KHz的响应信号如图6所示。从图6可以看出提取出的中心频率为25KHz的响应有3 个幅值较大的波包。在较低的频厚积范围内，由于&模式的Lamb波传播速度较快，所以先 到达传感器，所以第一个波包为对称&模式波包，第二个波包为反对称&模式波包，第三个 较大的波包为由反射造成的两种模式的混叠。然后用最大峰值法可求出&模式和&模式 两个波包到达0号压电片的时间差，再根据测得的两种模式的速度，即可用(15)式估计出 冲击源到0号压电片的距离屯为46. 4cm。按照相同的方法可估计出冲击源到达其它压电 片的距离。然后运用公式(16)计算出的距离R为48. 8cm，即为冲击源到达坐标参考点的距 离，从而实现了冲击定位。
权利要求
一种基于多重信号分类算法的结构冲击载荷定位方法，其特征在于包括如下步骤步骤一冲击响应信号的采集(1)在所述结构上布置由多个相同型号的压电元件组成的压电阵列；(2)将2块4通道的数采卡0通道并联在一起，保证两块数采卡同步触发采集；(3)设定采集参数，利用电荷放大器和数据采集卡采集压电阵列输出的冲击响应信号；步骤二基于MUSIC算法进行冲击源的波达方向估计首先用快速傅里叶变换对所述冲击响应信号进行频谱分析，然后根据选定的方向估计中心频率利用小波变换从所述冲击响应信号中提取出中心频率下的窄带信号，最后取窄带信号的幅值最大的信号运用多重信号分类(MUSIC)算法的进行DOA估计；步骤三基于Lamb波传播特性进行冲击源的距离估计首先用小波变换从任一压电元件接收到的冲击响应信号中提取出距离估计中心频率下的Lamb波信号，获得同一频率的对称S0模式的Lamb波和反对称A0模式的Lamb波，用最大峰值法求出对称S0和反对称A0这两种模式的Lamb波到所述达压电元件的时间差，最后由测得的这两种模式的Lamb波的传播速度求出冲击源到达压电元件的距离，将各压电元件估计出的距离值转化为相对于压电阵列参考点的距离，取所有压电元件测得的距离的平均值作为距离估计值。
2.根据权利要求1所述的基于阵列信号处理技术的冲击载荷定位监测方法，其特征在 于步骤一中，所述压电阵列中压电元件间的距离小于冲击源信号波长。
全文摘要
本发明公布了一种基于多重信号分类算法的结构冲击载荷定位方法，包括如下步骤步骤一冲击响应信号的采集；步骤二利用多重信号分类(MUSIC)算法进行冲击源的波达方向(DOA)估计；步骤三根据Lamb波传播特性进行冲击源距离的估计。本发明方法简便可靠，可对冲击进行实时、在线监测，冲击定位精度高、实时性强。
文档编号G01B7/02GK101865758SQ20101019927
公开日2010年10月20日 申请日期2010年6月12日 优先权日2010年6月12日
发明者张炳良, 梁栋, 苏永振, 袁慎芳, 邱雷 申请人:南京航空航天大学