CN109828033B - 基于振动响应相似度分析的损伤识别方法和*** - Google Patents

Abstract

本发明提供了一种基于振动响应相似度分析的损伤识别方法和***，包括：通过测试分别获得结构损伤前、后的振动响应时域数据，并通过拉式变换得到其频率响应数据；根据损伤前、后结构的频率响应数据，进行全频段相似度分析，得到结构上所有部位的全频段相似度和对应的阈值，全频段相似度小于自身阈值则为响应异常，对应结构的响应异常部位；对响应异常部位的频率响应数据进行分频段相似度分析，得到分频段相似度和对应的阈值，分频段相似度小于自身阈值的频段为响应异常频段；结合响应异常部位和异常频段评估结构损伤程度。本发明能有效的对不同形式和不同材料的结构进行单损伤、多损伤工况的识别。

Description

基于振动响应相似度分析的损伤识别方法和***

技术领域

本发明涉及测试测量技术领域，具体地，涉及基于振动响应相似度分析的损伤识别方法和***。

背景技术

损伤识别技术被广泛应用在土木工程、航空、航天领域，特别是桥梁结构、飞机机身、卫星主承力结构等。所涉及的结构形式复杂，有杆、梁、板、壳或其组合形式；涉及的材料多样化，有金属、非金属、复合材料，特别是复合材料因其优异的性能被大量应用在航空航天器上。而复合材料的损伤机理更加复杂，可能由于工艺、疲劳等产生脱层、脱粘、纤维断裂等损伤。因此损伤识别技术的研究意义很大。

损伤识别主要包含4个层次的内容：1)检测——识别到结构中的损伤是否发生；2)定位——识别结构的损伤位置；3)评估——识别结构的损伤程度；4)预测——对结构是否能继续服役的预测。但到目前为止，仅损伤检测与定位技术较为成熟，对损伤评估和结构寿命预测还有待进一步研究。

国内外学者提出了许多方法，包括静力识别法、动力指纹法、模型修正法、动静结合等，但现有方法仍存在诸多不足，如现有的损伤识别技术适用性差，大多仅适用形式简单的简单，如梁、板；仅适用于单一材料如金属，对复合材料损伤则无法识别。

发明内容

针对现有技术中的缺陷，本发明的目的是提供一种基于振动响应相似度分析的损伤识别方法和***。

根据本发明提供的一种基于振动响应相似度分析的损伤识别方法，包括：

数据获取步骤：通过测试分别获得结构损伤前、后的振动响应时域数据，并通过拉式变换得到其频率响应数据；

全频段相似度分析步骤：根据损伤前、后结构的频率响应数据，进行全频段相似度分析，得到结构上所有部位的全频段相似度和对应的阈值，全频段相似度小于自身阈值则为响应异常，对应结构的响应异常部位；

分频段相似度分析步骤：对响应异常部位的频率响应数据进行分频段相似度分析，得到分频段相似度和对应的阈值，分频段相似度小于自身阈值的频段为响应异常频段；

评估步骤：结合响应异常部位和异常频段评估结构损伤程度。

较佳的，获取结构损伤前、后的振动响应时域数据时，测点的数目和位置能代表结构的空间几何形状特征，至少能表征出结构前3阶模态振型。

较佳的，所述数据获取步骤中，通过测试分别获得结构损伤前、后的振动响应时域数据时的两次测试条件相同。

较佳的，所述数据获取步骤中，结构损伤前、后的振动响应时域数据的采样频率满足采样定理。

较佳的，所述全频段相似度分析步骤中，阈值的计算是在对结构的损伤评估基础上，基于选定的相似度方法进行导出或约定，依据如下公式：

式中，s_Eucld(X,Y)、s_Jaccard(X,Y)、s_Esim(X,Y)分别表示Eucld相似度、Jaccard相似度、Esim相似度，X＝(x₁x₂…x_n)、Y＝(y₁y₂…y_n)分别表示结构上损伤前、后的测点的频率响应数据，C_L表示相似度阈值；α为结构允许损伤程度下的响应幅值变化率；ω_i为比重因子。

较佳的，所述分频段相似度分析步骤中，对响应异常部位的频率响应数据进行分频段相似度分析，采用移动平均的方式，求各分频段频率响应数据的分频段相似度和阈值，依据所述全频段相似度分析步骤中公式，以频率为横轴，相似度为纵轴，绘制成分频段相似度和阈值曲线。

较佳的，所述评估步骤中，将计算得到的全频段相似度和对应的阈值，以及分频段相似度和对应的阈值转化得到损伤程度指标D用来衡量损伤程度，依据公式如下：

其中，s_i和C_iL分别表示测点i的相似度和对应的阈值，n为结构上所有的测点数。

较佳的，所述分频段相似度分析将全频段分成若干个子频段，子频段长度为全频段长度的0.01倍。

根据本发明提供的一种基于振动响应相似度分析的损伤识别***，包括：

数据获取模块：通过测试分别获得结构损伤前、后的振动响应时域数据，并通过拉式变换得到其频率响应数据；

全频段相似度分析模块：根据损伤前、后结构的频率响应数据，进行全频段相似度分析，得到结构上所有部位的全频段相似度和对应的阈值，全频段相似度小于自身阈值则为响应异常，对应结构的响应异常部位；

分频段相似度分析模块：对响应异常部位的频率响应数据进行分频段相似度分析，得到分频段相似度和对应的阈值，分频段相似度小于自身阈值的频段为响应异常频段；

评估模块：结合响应异常部位和异常频段评估结构损伤程度。

较佳的，获取结构损伤前、后的振动响应时域数据时，测点的数目和位置能代表结构的空间几何形状特征，至少能表征出结构前3阶模态振型。

与现有技术相比，本发明具有如下的有益效果：

1、适用于不同类型结构。不仅适用于常见的简单结构如梁、板等，也适用于复杂的结构如壳、桁架、板杆结构等。

2、适用于不同材料结构。同时适用于金属、非金属、复合材料等的损伤识别。

3、分频段相似度分析，从频率变化识别出结构动态特性改变，分析出损伤对结构***响应和安全性的影响。

4、可实现多损伤工况下的损伤定位和损伤程度评估。

附图说明

通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：

图1所示是本发明方法的流程框图；

图2所示是分频段相似度分析中“移动平均”方式的示意图；

图3所示是实施例1的卫星承力筒模型示意图；

图4所示是Eucld全频段相似度分析结果；

图5所示是Jaccard全频段相似度分析结果；

图6所示是Esim全频段相似度分析结果；

图7所示是响应异常部位(测点25)损伤前、后的频率响应曲线；

图8所示是响应异常部位(测点25)的Eucld分频段相似度曲线；

图9所示是响应异常部位(测点25)的Jaccard分频段相似度曲线；

图10所示是响应异常部位(测点25)的Esim分频段相似度曲线；

图11所示是5种损伤工况下的相似度损伤程度指标CL。

具体实施方式

下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。

如图1所示，本发明提供的一种基于振动响应相似度分析的损伤识别方法，包括：

数据获取步骤(以下简称：步骤(1))：通过测试分别获得结构损伤前、后的振动响应时域数据，并通过拉式变换得到其频率响应数据；

全频段相似度分析步骤(以下简称：步骤(2))：根据损伤前、后结构的频率响应数据，选取适当的相似度分析方法，进行全频段相似度分析，得到结构上所有部位的全频段相似度和对应的阈值，全频段相似度小于自身阈值则为响应异常，对应结构的响应异常部位；

分频段相似度分析步骤(以下简称：步骤(3))：选取适当的相似度分析方法，对响应异常部位的频率响应数据进行分频段相似度分析，得到分频段相似度和对应的阈值，分频段相似度小于自身阈值的频段为响应异常频段；

评估步骤(以下简称：步骤(4))：结合响应异常部位和异常频段评估结构损伤程度。

步骤(1)中，结构损伤前、后的振动响应数据采集，要求测点的数目和位置能代表结构的空间几何形状特征，至少能表征出结构前3阶模态振型。

步骤(1)中，结构损伤前、后的振动响应数据采集，要求前后两次试验条件完全相同，包括力学环境、采集参数设置等。

步骤(1)中，结构损伤前、后的振动响应数据可以为各种类型，包括位移、加速度、应变等。

步骤(1)中，结构损伤前、后的振动响应数据采样频率满足采样定理，使得振动响应数据包含足够的频域信息。

步骤(2)中，相似度分析方法可选Eucld相似度、Jaccard相似度、Esim相似度三种方法之一，但不限于该三种方法。若基于该三种方法，则结构损伤前、后的频率响应数据全频段相似度分析可依据如下公式：

式中，s_Eucld(X,Y)、s_Jaccard(X,Y)、s_Esim(X,Y)分别表示Eucld相似度、Jaccard相似度、Esim相似度，X＝(x₁x₂…x_n)、Y＝(y₁y₂…y_n)分别表示结构上损伤前、后的某测点的频率响应数据。

具体的，步骤(2)结构损伤前、后的频率响应数据全频段相似度阈值计算是在对结构的损伤评估基础上，基于选定的相似度方法进行合理导出或约定。此处可采用Eucld相似度、Jaccard相似度、Esim相似度三种方法之一，依据如下公式(4)，但不限于该三种方法,

其中，C_L表示相似度阈值；α为结构允许损伤程度下的响应幅值变化率；ω_i为比重因子。

步骤(3)中，对响应异常部位的频率响应数据进行分频段相似度分析需要对全频段分成若干个子频段，子频段长度大小要合适，一般取为全频段长度的0.01。

步骤(3)对响应异常部位的频率响应数据进行分频段相似度分析，应采用“移动平均”的方式，依据公式(1)～(4)，求各分频段频率响应数据的分频段相似度和阈值。以频率为横轴，相似度为纵轴，可以绘制成分频段相似度和阈值曲线。

其中，步骤(3)对响应异常部位的频率响应数据进行分频段相似度分析，“移动平均”的具体实现方式：全频段区间[f0,f1]内，从起始频率点f₀开始，将长度为w(w＜＜(f₁-f₀))的子频段编号为1，对其进行分频段相似度分析，然后以Δf(0＜Δf＜w)为步长向右移动该子频段，编号为2，继续进行相似度分析，以此类推，采用子频段“滑动”的方式进行分频段相似度分析，便可以得到不同子频段的相似度。如2所示表示子频段i和i+1的相对关系。

具体的，步骤(4)综合响应异常部位和异常频段评估结构损伤，步骤(2)计算得到的响应异常部位可以用来损伤定位；步骤(2)、(3)计算得到的相似度和其阈值，转化得到损伤程度指标D用来衡量损伤程度，依据公式如下：

其中，s_i和C_iL分别表示测点i的相似度和阈值，n为结构上所有的测点数。

下面结合具体工程实例对本发明进行说明。

实施例一：如图3所示为卫星承力筒，筒体柱段外径Φ1114mm，内径Φ1090mm，高3260mm，主体部分采用碳纤维/环氧-铝蜂窝夹层复合材料。承力筒上设置测点40个，从底部到顶部均分5周，每周布置8个测点。损伤区域位于测点25下方附近，损伤形式为蜂窝夹层结构的蒙皮芯子分层损伤。激励为承力筒底部正弦扫频激励，频带范围20～2000Hz，幅值大小10m/s²。

具体实施步骤如下：

步骤1：通过有限元模型仿真分析获得承力筒损伤前、后，所有测点处的振动加速度响应时域数据，然后通过拉式变换得到相应的响应频域数据；

步骤2：根据损伤前、后结构的振动频率响应数据，利用公式(1)～(4)进行全频段相似度分析，暂定α＝0.15，得到结构上所有部位的全频段相似度和对应的阈值。以测点编号为横轴，以相似度为纵轴，绘制成全频段相似度曲线，结果如图4～6所示。从图中可以看出，测点25的全频段相似度低于阈值，因此响应异常部位为测点25附近，准确识别出了损伤位置。

步骤3：对响应异常部位，即测点25的频率响应数据(见图7)进行分频段相似度分析，暂定α＝0.15，利用公式(1)～(4)，得到分频段相似度和阈值，以频率为横轴，以相似度为纵轴，得到分频段相似度曲线，如图8～10所示。从图中可以看到，阈值虚线下方的分频段相似度曲线对应的频段为响应异常频段，该信息对分析结构损伤有一定参考价值。

步骤4：设置如下表1的损伤工况，1～5分别表示损伤区域逐渐增大，利用公式(5)计算其损伤程度，得到结果图11。可以看出，随损伤区域增大，损伤程度值增大。该指标可以准确衡量损伤程度。

表1损伤工况

在上述一种基于振动响应相似度分析的损伤识别方法的基础上，本发明还提供一种基于振动响应相似度分析的损伤识别***，包括：

数据获取模块：通过测试分别获得结构损伤前、后的振动响应时域数据，并通过拉式变换得到其频率响应数据；

全频段相似度分析模块：根据损伤前、后结构的频率响应数据，进行全频段相似度分析，得到结构上所有部位的全频段相似度和对应的阈值，全频段相似度小于自身阈值则为响应异常，对应结构的响应异常部位；

分频段相似度分析模块：对响应异常部位的频率响应数据进行分频段相似度分析，得到分频段相似度和对应的阈值，分频段相似度小于自身阈值的频段为响应异常频段；

评估模块：结合响应异常部位和异常频段评估结构损伤程度。

本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的***及其各个装置、模块、单元以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的***及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的***及其各项装置、模块、单元可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

在本申请的描述中，需要理解的是，术语“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本申请和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本申请的限制。

以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。

Claims (7)

1.一种基于振动响应相似度分析的损伤识别方法，其特征在于，包括：

数据获取步骤：通过测试分别获得结构损伤前、后的振动响应时域数据，并通过拉式变换得到其频率响应数据；

全频段相似度分析步骤：根据损伤前、后结构的频率响应数据，进行全频段相似度分析，得到结构上所有部位的全频段相似度和对应的阈值，全频段相似度小于自身阈值则为响应异常，对应结构的响应异常部位；

分频段相似度分析步骤：对响应异常部位的频率响应数据进行分频段相似度分析，得到分频段相似度和对应的阈值，分频段相似度小于自身阈值的频段为响应异常频段；

评估步骤：结合响应异常部位和异常频段评估结构损伤程度；

所述全频段相似度分析步骤中，阈值的计算是在对结构的损伤评估基础上，基于选定的相似度方法进行导出或约定，依据如下公式：

式中，s_Eucld(X,Y)、s_Jaccard(X,Y)、s_Esim(X,Y)分别表示Eucld相似度、Jaccard相似度、Esim相似度，X＝(x₁ x₂…x_n)、Y＝(y₁ y₂…y_n)分别表示结构上损伤前、后的测点的频率响应数据，C_L表示相似度阈值；α为结构允许损伤程度下的响应幅值变化率；ω_i为比重因子；

所述分频段相似度分析步骤中，对响应异常部位的频率响应数据进行分频段相似度分析，采用移动平均的方式，求各分频段频率响应数据的分频段相似度和阈值，依据所述全频段相似度分析步骤中公式，以频率为横轴，相似度为纵轴，绘制成分频段相似度和阈值曲线；

所述评估步骤中，将计算得到的全频段相似度和对应的阈值，以及分频段相似度和对应的阈值转化得到损伤程度指标D用来衡量损伤程度，依据公式如下：

其中，s_i和C_iL分别表示测点i的相似度和对应的阈值，n为结构上所有的测点数。

2.根据权利要求1所述的基于振动响应相似度分析的损伤识别方法，其特征在于，获取结构损伤前、后的振动响应时域数据时，测点的数目和位置能代表结构的空间几何形状特征，至少能表征出结构前3阶模态振型。

3.根据权利要求1所述的基于振动响应相似度分析的损伤识别方法，其特征在于，所述数据获取步骤中，通过测试分别获得结构损伤前、后的振动响应时域数据时的两次测试条件相同。

4.根据权利要求1所述的基于振动响应相似度分析的损伤识别方法，其特征在于，所述数据获取步骤中，结构损伤前、后的振动响应时域数据的采样频率满足采样定理。

5.根据权利要求1所述的基于振动响应相似度分析的损伤识别方法，其特征在于，所述分频段相似度分析将全频段分成若干个子频段，子频段长度为全频段长度的0.01倍。

6.一种基于振动响应相似度分析的损伤识别***，其特征在于，包括：

数据获取模块：通过测试分别获得结构损伤前、后的振动响应时域数据，并通过拉式变换得到其频率响应数据；

全频段相似度分析模块：根据损伤前、后结构的频率响应数据，进行全频段相似度分析，得到结构上所有部位的全频段相似度和对应的阈值，全频段相似度小于自身阈值则为响应异常，对应结构的响应异常部位；

分频段相似度分析模块：对响应异常部位的频率响应数据进行分频段相似度分析，得到分频段相似度和对应的阈值，分频段相似度小于自身阈值的频段为响应异常频段；

评估模块：结合响应异常部位和异常频段评估结构损伤程度；

所述全频段相似度分析模块中，阈值的计算是在对结构的损伤评估基础上，基于选定的相似度方法进行导出或约定，依据如下公式：

式中，s_Eucld(X,Y)、s_Jaccard(X,Y)、s_Esim(X,Y)分别表示Eucld相似度、Jaccard相似度、Esim相似度，X＝(x₁ x₂…x_n)、Y＝(y₁ y₂…y_n)分别表示结构上损伤前、后的测点的频率响应数据，C_L表示相似度阈值；α为结构允许损伤程度下的响应幅值变化率；ω_i为比重因子；

所述分频段相似度分析步骤中，对响应异常部位的频率响应数据进行分频段相似度分析，采用移动平均的方式，求各分频段频率响应数据的分频段相似度和阈值，依据所述全频段相似度分析步骤中公式，以频率为横轴，相似度为纵轴，绘制成分频段相似度和阈值曲线；

所述评估步骤中，将计算得到的全频段相似度和对应的阈值，以及分频段相似度和对应的阈值转化得到损伤程度指标D用来衡量损伤程度，依据公式如下：

其中，s_i和C_iL分别表示测点i的相似度和对应的阈值，n为结构上所有的测点数。

7.根据权利要求6所述的基于振动响应相似度分析的损伤识别***，其特征在于，获取结构损伤前、后的振动响应时域数据时，测点的数目和位置能代表结构的空间几何形状特征，至少能表征出结构前3阶模态振型。

Images