一种基于有效信号段自提取算法和分段截取相关系数算法的椭圆定位算法的制作方法

本发明属于机械损伤诊断技术领域，涉及一种基于有效信号段自提取算法和分段截取相关系数算法的椭圆定位算法。

背景技术：

lamb波是一种超声导波，其长距离传播，多模态，频散等特性使得lamb波广泛应用于无损检测领域。

现有损伤定位原理多采用rapid方法进行损伤定位，对其定位精度起决定作用的是传感器数量，即传感器数量越多，定位精度越精确。由于其成像算法要求，传感器数目一般应大于八块。rapid算法主要是采用的信号特征是sdc(signaldifferencecoefficient，信号差异系数)，采取对参考信号和对损伤信号在统计量上的比较，进而得出损伤部分的位置。

上述算法存在着许多不足：

1)对数据的整理方法不足，无法有效地利用数据。

2)所用传感器数目过多，成本较高，难度过大。

3)概率定位不准确，无法精确定位。

4)虽然原理简单，但在具体实施过程中操作困难，难以实现。

技术实现要素：

本发明通过运用有效信号段自提取算法以及分段截取相关系数算法来对所提取的信号进行操作处理，然后通过椭圆定位算法，根据多个椭圆的公共交点来确定损伤位置。在保证精度的同时使用更少的传感器，可以减少成本，降低实验难度。

本发明采用的技术方案包括如下步骤：

步骤1：利用激励传感器激励出lamb波，同时，接收传感器开始接收信号。记录接收传感器在采集时间内采集到的所有信号。

步骤2：lamb波在传播时中需要一段在结构中传播的时间，这段时间理论上表现在接收信号中就是一段幅值为0的信号段，由于实际环境中噪声、信号干扰的影响，实际这段信号的幅值通常不为0。但信号噪声通常频率相对较低，lamb波为高频激励信号，可以通过滤波去除噪声。这个信号段在计算信号相关系数时因为不携带任何与结构有关的信息，可以被认为是无关信号。同理，lamb波传播到边界或者其他传感器时，其模态会发生转变，甚至出现信号混叠等现象，较为复杂，难以处理。因此这段信号对于计算信号相关系数是不利的，也是无关的。所以，实际的有效信号段就是在检测范围内，从接收传感器接收到lamb波到接收传感器接收完lamb波的这段信号，也就是接收传感器在采集时间内采集到的信号去掉无关信号后，截取出来的信号。这段有效信号包含了lamb波在无损结构中传播的信号和遇到结构损伤时的散射信号。当以较短的信号段参与损伤指标的计算时，距离激励传感器较远的接收传感器接收到的信号中可能会缺失一部分有效信号；若以较长的信号段参与损伤指标的计算，距离激励传感器较近的接收传感器接收到的信号中可能会包含较多的无关信号。因此，信号段长度过长或过短都不利于基于信号相关系数的损伤指标的计算。信号中携带结构损伤信息的有效信号越多，得出的信号相关系数就越有效。

接收传感器在接收到lamb波之前采集到的信号和从边界反射回来的信号都属于无关信号，这种无关信号越少，对损伤指标的不利影响也就越小，采用等长度信号段来得到损伤指标时，改变信号段的长度，就会改变参与计算的信号段中的有效信号段和无关信号段的长度，使得基于信号相关系数的损伤指标发生改变，损伤识别结果中将有可能出现损伤识别虚假和识别不明显的现象。针对等长度信号段参与计算的不足，提出了一种有效信号段自提取的算法。有效信号段自提取算法提取的关键参数是接收传感器接收到lamb波的时刻和完全接收完不包括边界和其他传感器反射的lamb波的时刻。因此在传感器布置时，应考虑传感器之间的距离和传感器与边界的距离，即要求速度最快的s0波从激励传感器出发，经边界或相邻传感器反射到接收传感器所需的时间尽可能地多于速度最慢的a0波从出发到接收完全所需要的时间。根据lamb波在结构的传播速度，结合压电传感器之间的距离，计算速度最快的s0波从出发到刚接收所需要的时间，转化为相对应的信号采集点数m1，计算速度最慢的a0波从出发到刚接收所需要的时间，转化为相对应的信号采集点数m2，并补偿lamb波的波宽长度对应的信号采集点数n，为减小误差，乘一个放大倍数q对采集到的信号进行截断处理，提取出有效信号段的过程如下：从信号段第m1个点数作为起点，从信号段第m2+q×n个点数作为终点，起点至终点的这段信号段就是有效信号段，有效信号段长度m2+q×n-m1+1，如图1所示。

步骤3：对提取的信号利用分段截取相关系数算法计算出不同部分的信号相关系数。由于lamb波在经过厚度突变或分层损伤出会发生模态转换效应，通过对健康信号和损伤信号进行比对，可识别出因模态转换而新激发的模态的出现时间，进而判断结构中损伤的位置。

利用矩形窗对传感信号进行分段截取，得到分段后健康与损伤信号的相关系数，度量其线性相关性，以相关系数作为判断依据识别损伤出现时间，进而确定损伤位置。对传感信号分段处理后可极大增强损伤处信号的差异性，损伤指标的确定也将更加直观。

由于lamb波信号经过结构损伤时发生散射和能量衰减，信号幅值会发生明显的改变，通过比较结构在发生损伤前后的lamb波信号，就能确定幅值的变化量，进而判断出结构中是否存在损伤。另一种处理方法则是对lamb波信号进行时频分析，把时域内的lamb波信号转换为频域或者时频域内的特征进行分析，提取出的幅值变化量就作为损伤因子。利用损伤发生前后的信号幅值的相关系数作为损伤的判断依据来判断损伤和损伤成像。某条路径采集到的损伤前后的信号段经截断处理后分别用a和b表示，ρ表示信号a和b的相关系数，对于离散信号ρ可表示为：

式中：分别为信号a、b的均值；aj、bj分别为a、b的第j个值；num为有效信号段的长度。损伤指标di则定义为

di＝1-|ρ|

因此di的数值越大，则该路径损伤前后的信号相关性越小，损伤位置范围就越大。而后通过对di值的合理选择，可以直观的判断损伤位置。如图3所示。

步骤4：确定存在损伤之后，利用椭圆定位算法来确定损伤位置。在lamb波激励器s0处激发应力波，该应力波可从s0直接传播到信号接收传感器s1处，也可从s0经损伤处f反射后再到达s1。由此，结构损伤的反射信号ri(t)为

ri(t)＝xi(t)-xbj(t)(1.1)

(1.1)式中：xi(t)为结构中存在损伤时传感器接收到的信号；xbi(t)为结构中无损伤时传感器接收到的信号；i为传感器编号。进一步，反射信号的传播时间ti为

(1.2)式中：是lamb波从激励传感器通过结构损伤反射传播到接收器所经历的距离；c0为lamb波传播的群速度。lamb波在各向同性介质中的传播速度c0由材料的物理特性及lamb的频率决定，与传感器的布置位置无关。因此可事先进行标定。当c0确定后，ti可通过传感器测得，则‖s0f‖+‖sif‖为定值。由椭圆几何定理可知，损伤f在以s0和si为焦点的椭圆上。若s2为另一传感器，则损伤点f也在以s0和s2为焦点的椭圆上。两个椭圆的焦点即为损伤点f的所在位置。如图4所示。

当传感器为4块时，如图5所示，可构造6个椭圆，这些椭圆在4块传感器内部相交于较小的一个范围内，也可以确定损伤的位置。相较于使用3块传感器，使用4块传感器的时候，仅多使用了一个传感器，就可以将探测的范围扩大一倍。个数/面积的值增大，效率更高。

附图说明

图1为有效信号段的提取

图2为信号提取示意图

图3为自提取之后波形与分段相关系数图

图4为三块传感器椭圆定位法示意图

图5为四块传感器椭圆定位法示意图

具体实施方式

下面对本发明进行详细的说明，所述是对本发明的解释而不是限定。

一种基于有效信号段自提取算法和分段截取相关系数算法的椭圆定位算法是建立在原有的椭圆定位算法的基础上，充分利用有效信号段自提取算法和分段截取相关系数算法来更加简便的确定损伤的位置，包括以下步骤：

1.在待测机械结构上按矩形的方式安装四个压电传感器，并设置边界；

2.在机械结构件的初始无损伤状态下，依次以压电传感器中的一个为激励端进行信号的激励，以其余的压电传感器为接收端接收信号，每两个压电传感器的传感路径上传播的信号都被采集，所采集到的信号作为参考信号；

3.当机械结构件产生损伤后，分别依次以压电传感器中的一个为激励端进行信号的激励，以其余的压电传感器为接收端接收信号，每两个压电传感器的传感路径上传播的信号都被采集；

4.调整安装压电传感器之间的距离，应使速度最快的波s0从激励传感器出发，经损伤或边界到达接收传感器的时间晚于波a0从激励传感器直达接收传感器的时间的平方；

5.对信号数据进行拦截：以激励端与接收端的直线距离l除以lamb波在有效频段内的最大传播速度vmax得到的时间点，作为有效数据起始点t，以激励端与接收端的直线距离l乘rapid算法的尺度参数p得到计算区域内的最大路径长度j3l，用该长度除以lamb波在有效频段内的最小传播速度νmin加上激励信号的半个周期t，得到的时间点，作为有效数据终止点t1。如图2；

6.由理论分析与实验可得lamb波s0、a0模态波速，结合传感器之间的距离，可得s0由激励端-接收端直达波的到达时间t1；a0直达时间t2，将lamb波波宽长度记为δt，为减少误差，乘一个优化系数n，并提取数据；

7.对健康信号和损伤信号进行比对，利用矩形窗对传感信号进行分段截取，得到分段后健康与损伤信号的相关系数，度量其线性相关性，以相关系数作为判断依据识别损伤出现时间，进而确定损伤位置；

某路径采集到的损伤前后信号段经分段截取进行损伤指标di定义，所截取健康损伤信号段的差异性越大，相关性越小。即di值越大，损伤范围越大。通过对di值的合理选择，可以直观的判断是否存在损伤，如图3；

确定存在损伤之后，对传感器由两个传感器形成一个椭圆，更为精确的使用四个传感器来规划椭圆，多个椭圆的交点，即为损伤点。

具体地，在本实验方案中，实验研究对象分别采用3块和4块压电传感器，对本文所研究的基于有效信号段自提取算法和分段截取相关系数算法的椭圆定位算法监测方法进行实验研究。实验试件采用铝板。采用p2压电陶瓷片组成激励-传感阵列，压电片直径10mm，厚度1mm，在板的中心位置对压电陶瓷进行布置。损伤采用粘质量块/磁铁/打孔，产生局部结构异化，对lamb波传播而言，该损伤界面将发生与真实损伤相类似的信号散射现象。以铝板中心点为零点构建坐标系，压电阵列坐标为，模拟损伤位置坐标为，压电阵列坐标和模拟损伤位置坐标的位置及编号示意图如图所示。实验设备由函数发生器afg1062、数字示波器rtb2004，多通道阵列扫查系统，计算机和传感器阵列组成。

具体步骤如下：

1.在待测机械结构上安装传感器，以矩形进行布置，组成设备网络；

2.调整安装压电传感器之间的距离，应使速度最快的s0波从激励传感器出发，经损伤或边界到达接收传感器的时间晚于a0波从激励传感器直达接收传感器的时间的平方；

3.由理论分析与实验可得lamb波s0、a0模态波速，结合传感器之间的距离，可得s0波由激励端-接收端直达波的到达时间t1；a0波直达时间t2；将lamb波波宽长度记为δt，为减少误差，乘一个优化系数h，公式为

4.对接收信号进行提取处理，提取过程如下：将t1设为有效信号段起点，t2+δt×n设为终点，则有效信号段为t2+δt×n-t1，如图2所示。

5.使用分段截取相关系数的方法，利用矩形窗对传感信号进行分段截取，得到分段后健康与损伤信号的相关系数，度量其线性相关性，以相关系数作为判断依据识别损伤出现时间，通过对健康信号和损伤信号进行比对，可识别出因模态转换而新激发的模态的出现时间，进而判断结构中损伤的位置。

6.在没有出现异常比较差异时，对前面每段所截取信号的相关系数进行综合数据整理，即为下一段信号数据出现时的指导值，根据指导值和后一段的信号数据的相关系数比较，再根据相关系数比较的结果值进行综合分析，若所截取健康损伤信号段的差异性越大，相关性越小即di值越大，通过对di值的合理选择，可以直观的判断大体损伤位置，如图3所示；

7.根据所获得的数据，进一步的确定损伤位置，如图4所示，当传感器为3块时，在lamb波激励器s0处激发应力波，该应力波可从s0直接传播到信号接收传感器s1处，也可从s0经损伤处f反射后再到达s1。

8.由椭圆定位算法可知，损伤f在以s0和s1为焦点的椭圆上。若s2为另一传感器，则损伤点f也在以s0和s2为焦点的椭圆上，而且损伤点f也在以s1、s2为焦点的椭圆上。三个椭圆的交点即为损伤点f的所在位置。

如图5所示当传感器数量为4块时，损伤处也是在以两个传感器为焦点的多个椭圆的交点处。同时，可以明显的发现，使用四块传感器所检测的范围相较于使用三块整整扩大了一倍，个数/面积的值增大，效率显著提高。