一种面向复合材料薄壁结构的被动式冲击损伤识别方法及系统技术方案

本发明专利技术提供了一种面向复合材料薄壁结构的被动式冲击损伤识别方法和系统，采集冲击产生的Lamb波信号；利用变分模态分解(VMD)将Lamb波信号分解成多模态信号分量，运用连续小波变换(CWT)获取中心频率最高的模态信号分量的时频谱；利用DBSCAN聚类算法判别冲击损伤有无；采用归一化方差序列测量Lamb波信号的波达时间；使用迭代圆弧轨迹成像算法确定冲击损伤位置在结构监测区域内的概率分布。本发明专利技术利用归一化方差序列提取Lamb波的波达时间特征过程，无需解释信号的成分，受到环境噪声干扰较小，采用迭代圆弧轨迹成像算法在减少了计算时间成本的情况下同时实现了冲击损伤的精确定位。位。位。

【技术实现步骤摘要】
一种面向复合材料薄壁结构的被动式冲击损伤识别方法及系统


[0001]本专利技术涉及复合材料结构健康监测领域，尤其涉及一种面向复合材料薄壁结构的被动式冲击损伤识别方法及系统。

技术介绍

[0002]复合材料薄壁结构在航空航天、能源动力和轨道交通等领域扮演着重要的角色。然而复合材料薄壁结构在运输和服役过程中不可避免的会遭受冲击载荷的作用，从而对结构造成损伤，直接威胁结构安全性和完整性。人工定期对结构进行无损检测维护成本较高，且需要消耗大量时间。因此，对复合材料薄壁结构进行实时在线的冲击损伤监测，对维护结构的健康安全具有重要意义。
[0003]基于Lamb波的薄壁结构损伤识别技术因其成本低、传播距离远、覆盖面积大等优势，在结构健康监测领域中得到了广泛应用。基于Lamb波的损伤识别技术主要分为主动式和被动式两种。其中，基于Lamb波的主动式损伤识别技术是通过在薄壁结构表面安装驱动器主动激励Lamb波信号，然后使用接收传感器采集经过监测区域的Lamb波信号，接着应用信号处理方法提取出Lamb波与结构损伤相互作用的损伤指数，最后利用损伤识别算法对监测结构中的损伤进行识别。而基于Lamb波的被动式损伤识别技术只需要设置接收传感器采集冲击在薄壁结构中产生的Lamb波信号，应用信号处理方法提取冲击产生的Lamb波信号特征结合特定算法实现冲击损伤识别。基于Lamb波的被动式损伤识别技术相比主动式损伤识别技术，只需接收器而无需驱动器，更加简单，快速和高效，在冲击损伤识别方面更加具有优势，而目前国内外关于被动式损伤识别技术研究较少。所以对于复合材料薄壁结构，利用被动式损伤识别技术判别冲击损伤的有无，确定冲击损伤位置在结构监测区域内的概率分布无疑是一个挑战。

技术实现思路

[0004]根据现有技术存在的问题，本专利技术公开了一种面向复合材料薄壁结构的被动式冲击损伤识别方法及系统，技术方案如下：
[0005]本专利技术在第一方面提供了一种面向复合材料薄壁结构的被动式冲击损伤识别方法，其包括步骤：采集冲击在所述复合材料薄壁结构的结构检测区域内产生的Lamb波信号；利用变分模态分解将所述Lamb波信号分解成多个模态分量，运用连续小波变换得到每个所述多个模态分量的时频谱；基于其中最高频率模态分的时频谱利用DBSCAN聚类算法判别所述冲击损伤有无；若有则采用归一化方差序列测量所述Lamb波信号的波达时间；以及使用迭代圆弧轨迹成像算法基于所述波达时间确定冲击损伤位置在所述结构监测区域内的概率分布，将概率最大的位置作为预测的冲击损伤位置。
[0006]在一些实施例中假设所述Lamb波信号f(t)被分解为k个模态分量，保证分解序列为具有中心频率的有限带宽的模态分量，同时满足各模态的估计带宽之和最小，约束条件
为所有模态之和与原始信号相等，则相应约束变分表达式为：
[0007][0008]其中，K为分解的总模态数，{u
k
}和{w
k
}分别为第k个模态分量和对应的中心频率，k越小中心频率越高，δ(t)是Dirac函数，t是时间采样点，j是复数，e是自然常数；求解上述约束变分表达式即可获得所述Lamb波信号的多个模态分量。
[0009]在一些实施例中，采用小波变换获取中心频率最高的模态分量u1(t)的时频谱；对于空间下信号u1(t)的连续小波变换定义如下：
[0010][0011]其中，CWT()为连续小波变换函数，a为尺度因子，b为平移因子，“<>”表示内积，
“‑”
表示共轭，ψ(t)为母小波函数。
[0012]在一些实施例中，于其中最高频率模态分的时频谱利用DBSCAN聚类算法判别所述冲击损伤有无具体包括：提取中心频率最高的模态分量u1(t)的时频谱像素最大部分对应的样本点；初始化DBSCAN聚类算法的参数R和MinPoints；根据上述样本点的自动分类，判定结构是否为无损。
[0013]在一些实施例中，采用归一化方差序列测量所述Lamb波信号的波达时间包括计算所述Lamb波信号的方差序列；按所述方差序列的最大值进行归一化得到经归一化的信号；以及选择一个数值较低的阈值，将所述经归一化的信号到达该阈值的时间作为测量的波达时间。
[0014]在一些实施例中，所述使用迭代圆弧轨迹成像算法基于所述波达时间确定冲击损伤位置在所述结构监测区域内的概率分布，将概率最大的位置作为预测的冲击损伤位置包括：将所述测量的波达时间进行排序，记录波达时间最早的多个波达时间；对结构监测区域进行网格划分，建立笛卡尔坐标系，标记所述最早的多个波达时间的对应的多个传感器的位置坐标，假定所述网格的某一点为冲击损伤点，计算其至所述多个传感器的距离；设置合理迭代区间，以设定的迭代初始值和迭代步长确定多个冲击发生时刻；基于所述距离迭代计算所述多个冲击发生时刻对应的圆弧轨迹概率成像结果；以及将成像结果图像中像素最大值的网格点位置作为概率最大的冲击损伤的位置。
[0015]在一些实施例中，所述多个传感器位于阵列设置的传感器网络中，所述传感器网络用于采集所述结构监测区域的信号。
[0016]本申请的另一些实施例还提供了对应的面向复合材料薄壁结构的被动式冲击损伤识别系统，其包括被配置为采集冲击在所述复合材料薄壁结构的结构检测区域内产生的Lamb波信号的信号采集模块；被配置为利用变分模态分解将所述Lamb波信号分解成多个模态分量，运用连续小波变换得到每个所述多个模态分量的时频谱的信号处理模块；被配置为基于其中最高频率模态分的时频谱利用DBSCAN聚类算法判别所述冲击损伤有无的损伤判别模块；被配置为采用归一化方差序列测量Lamb波信号的波达时间的特征提取模块；以
及被配置为使用迭代圆弧轨迹成像算法基于所述波达时间确定冲击损伤位置在所述结构监测区域内的概率分布，将概率最大的位置作为预测的冲击损伤位置的冲击损伤概率分布获取模块。
[0017]由于采用了上述技术方案，本专利技术提供的一种面向复合材料薄壁结构的被动式冲击损伤识别方法及系统，与现有技术相比，本专利技术的有益效果是：
[0018]一方面，本专利技术所述的方法、系统针对复合材料薄壁结构提出了一种被动式冲击损伤识别方法，无需驱动传感器，只需要接收传感器，方法简单、快速和高效。
[0019]另一方面，本专利技术所述的方法、系统既适用于复合材料平板，也适用于加筋的复合材料板。
[0020]又一方面，本专利技术所述的方法、系统结合变分模态分解与连续小波变换，提供了一种基于DBSCAN的聚类算法判别冲击损伤有无，为无监督学习，无需对数据进行标签。
[0021]再一方面，本专利技术所述的方法、系统在特征提取模块采用归一化方差序列测量Lamb波信号的波达时间，无需解释信号的成分，受环境噪声干扰小。
[0022]再一方面，本专利技术所述的方法、系统提供了一种迭代圆弧轨迹成像算法，可以在减小计算时间成本的前提下确定冲击损伤位置在监测区域内的概率分布。
附图说明
[0023]为了更清楚地本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种面向复合材料薄壁结构的被动式冲击损伤识别方法，其特征在于：包括：步骤S1，采集冲击在所述复合材料薄壁结构的结构检测区域内产生的Lamb波信号；步骤S2，利用变分模态分解将所述Lamb波信号分解成多个模态分量，运用连续小波变换得到每个所述多个模态分量的时频谱；步骤S3，基于其中最高频率模态分的时频谱利用DBSCAN聚类算法判别所述冲击损伤有无，若有损伤则执行步骤S4和步骤S5；步骤S4，采用归一化方差序列测量所述Lamb波信号的波达时间；以及步骤S5，使用迭代圆弧轨迹成像算法基于所述波达时间确定冲击损伤位置在所述结构监测区域内的概率分布，将概率最大的位置作为预测的冲击损伤位置。2.根据权利要求1所述的面向复合材料薄壁结构的被动式冲击损伤识别方法，其特征在于：步骤S2包括：假设所述Lamb波信号f(t)被分解为k个模态分量，保证分解序列为具有中心频率的有限带宽的模态分量，同时满足各模态的估计带宽之和最小，约束条件为所有模态之和与原始信号相等，则相应约束变分表达式为：其中，K为分解的总模态数，{u
k
}和{w
k
}分别为第k个模态分量和对应的中心频率，k越小中心频率越高，δ(t)是Dirac函数，t是时间采样点，j是复数，e是自然常数；求解上述约束变分表达式即可获得所述Lamb波信号的多个模态分量。3.根据权利要求2所述的面向复合材料薄壁结构的被动式冲击损伤识别方法，其特征在于：步骤S2包括采用小波变换获取中心频率最高的模态分量u1(t)的时频谱；对于空间下信号u1(t)的连续小波变换定义如下：其中，CWT()为连续小波变换函数，a为尺度因子，b为平移因子，“< >”表示内积，
“‑”
表示共轭，ψ(t)为母小波函数。4.根据权利要求1所述的面向复合材料薄壁结构的被动式冲击损伤识别方法，其特征在于：步骤S3具体包括：步骤S301，提取中心频率最高的模态分量u1(t)的时频谱像素最大部分对应的样本点；步骤S302，初始化DBSCAN聚类算法的参数R和MinPoints；步骤S303，根据上述样本点的自动分类，判定结构是否为无损，若是则步骤S3终止，若否则进行步骤S4和步骤S5。5.根据权利要求1所述的面向复合材料薄壁结构的被动式冲击损伤识别方法，其特征在于：所述步骤S4具体包括步骤S401，计算所述Lamb波信号的方差序列；步骤S402，按所述方差序列的最大值进行归一化得到经归一化的信号；以及步骤S403，选择一个数值较低的阈值，将所述经归一化的信号到达...

【专利技术属性】
技术研发人员：杨雷，邓德双，杨正岩，曾旭，武湛君，马书义，
申请(专利权)人：大连君晟科技有限责任公司，
类型：发明
国别省市：