一种基于阵列超声表面波的构件损伤概率成像定位方法与流程

本发明专利涉及无损检测及结构健康监测领域，特别是一种基于阵列超声表面波的构件损伤概率成像定位方法。

背景技术：

航空航天结构件、压力容器、风电叶片等关键构件在服役运行过程中，常因工作载荷、冲击、振动、及环境变化等原因，产生疲劳裂纹、腐蚀、冲击损伤等近表面损伤。为避免出现重大故障乃至安全事故，对这些构件进行损伤定位和监测是非常必要的。损伤定位有助于获取构件的损伤分布情况，进而评估构件安全状态或通过局部检测手段进一步探知损伤特性，对于保障关键构件的运行安全具有重要意义。

目前通过常规超声检测手段进行逐点扫查可获得较准确的损伤位置，但检测费时，不适用于在役复杂构件。超声相控阵检测法在检测效率上有较大提高，但仍需要进行扫查检测，获取的大量数据所需处理时间长，且定位效率仍不高。声发射定位是一种高效的损伤定位方法，适用于大型复杂构件，但声发射检测需要对构件加载使其内部裂纹扩展产生声发射信号，易对被测构件造成破坏。兰姆波传播距离远，但仅能在薄板和薄管等结构中传播，且传播时存在频散现象，模态较为复杂。超声表面波相对于体波，能量集中于构件近表面，传播距离较远，对构件近表面损伤敏感，对构件形状不敏感，已广泛应用于构件的近表面损伤如裂纹、腐蚀、冲击损伤等的检测。

在构件的损伤定位算法方面，申请公开号cn102998369a，申请公开日2013年3月27日的专利文献公布了一种二维损伤定量化检测方法，采用渡越时间的定位方法对损伤进行定位，但由于定位曲线数量少，不能检测定位出形状复杂的损伤，且系统对测量噪声、渡越时间的测量误差及不确定因素的抗干扰能力不强，难以获得良好定位效果；申请公开号cn104343043a，申请公开日2016年2月24日的专利公布了一种基于abaqus的金属薄板微裂纹时间反转定位方法，采用基于时间反转的损伤定位方法，用时反特征信号对仿真构件重新激励，在时间和空间上的聚焦进行损伤定位，但需要理想的时反函数才能还原构件中的能量空间分布，且仿真构件模型与实际构件材料结构上的偏差，往往不能很好地聚焦定位。

技术实现要素：

本发明针对现有构件损伤定位方法存在的检测定位效率不高，受测量误差影响较大，且只能检测到单个简单损伤的问题，提出一种基于阵列超声表面波的构件损伤概率成像定位方法，实现构件损伤的高效检测和快速定位。

为了实现上述技术目的，本发明的技术方案是，一种基于阵列超声表面波的构件损伤概率成像定位方法，包括以下步骤：

步骤一、构建压电换能器检测阵列，在构件表面矩形被测范围四周均匀布置n个带楔块的压电换能器，记为1,2,...,n，以矩形中心为原点建立平面直角坐标系xoy，并分别获取发射换能器中心与损伤边界、损伤边界与接收换能器中心、发射换能器中心与接收换能器中心之间的距离

步骤二、激励和采集超声表面波信号，采用两种激励方式激励压电换能器ai(xi,yi)(i＝1,2,...,n)产生超声信号，在构件表面产生超声表面波信号，其中一种激励方式采用任意信号发生器产生汉宁窗调制的正弦信号，并经射频功率放大器后激励超声换能器，另一种激励方式用超声脉冲发生接收仪产生脉冲信号激励超声换能器，采集正弦信号激励方式的由其他换能器rj(xj,yj)(j＝1,2,...,n,j≠i)直接接收的超声表面波信号，记为si-j，采集脉冲信号激励方式的经损伤边界反射后由其他压电换能器rj(xj,yj)(j＝1,2,...,n,j≠i)接收的超声表面波信号，记为ei-j；

步骤三、计算损伤存在概率源值，将所述步骤二中的直接接收模式信号si-j和反射回波模式信号ei-j，分别按以下(1)、(2)进行基于信号能量的损伤概率成像算法和基于渡越时间的损伤概率成像算法，计算检测点(x,y)处的损伤存在概率源值：

(1)对所述步骤二中的直接接收模式信号si-j进行基于信号能量的损伤概率成像算法，按以下步骤进行：

步骤1、信号si-j经傅里叶变换后，获取构件无损伤时中心频率上阵列超声表面波的时间-能量波形分布曲线，作为基准信号e0，保存每一对ai-rj压电换能器的能量波形数据；

步骤2、重复上述步骤1，获取并保存构件有损伤时的每一对ai-rj压电换能器的时间-能量波形数据，作为检测信号e'，然后根据公式计算检测信号e'与基准信号e0的相异系数q；

步骤3、设发射压电换能器坐标为(xi,yi)，接收压电换能器坐标为(xj,yj)，待测点坐标为(x,y)，待测点到两压电换能器连线的距离为d,定义距离系数δ＝0.05，距离函数根据上述步骤2所得结果，再根据公式i′i-j(x,y)＝q×ψ(x,y)，计算得到ai-rj检测路线的检测范围内的损伤存在概率值i′i-j(x,y)，式中

步骤4、对其他每一对ai-rj压电换能器重复进行上述步骤1-3的处理；

步骤5、定义强化系数β为1.02，将上述步骤4得到的损伤存在概率值i′i-j(x,y)根据对应检测点(x,y)按公式进行累加，重建直接接收模式下所有ai-rj检测路线确定的损伤存在概率源值p′i-j(x,y)和损伤存在概率源图ⅰ；

(2)对所述步骤二中的反射回波模式信号ei-j基于渡越时间的损伤概率成像算法，按以下步骤进行：

步骤1、设反射回波模式信号ei-j的激励时刻为ti，接收时刻为tj，从激励时刻起，经损伤边界反射，到接收时刻止，表面波信号经历的渡越时间为t″i-j，路程为统计所有ai-rj检测路线的反射回波模式信号峰值时刻t的方差σ，将反射回波模式信号时刻tj的瞬时值按正态分布处理，取接收回波信号98％能量部分对应的时间区间(ti-j,t′i-j)，将接收时间区间减去激励时刻，得到渡越时间区间(ti-j,t′i-j)，由公式绘得椭圆环检测区；

步骤2、由上述步骤1确定的不同渡越时间t″i-j值所占整个渡越时间区间(ti-j,t′i-j)的权重为一条ai-rj检测路线确定的定位椭圆环检测区内检测点(x,y)处的损伤存在概率值

步骤3、对其他换能器对ai-rj重复进行上述步骤1-2的处理；

步骤4、将上述步骤3所得的所有ai-rj检测路线确定的定位椭圆环检测区包含的检测点(x,y)的损伤存在概率值i″i-j(x,y)，对应到被测工件表面坐标系上，按公式进行累加计算，得到反射回波模式所有ai-rj检测路线的损伤存在概率源值p″i-j(x,y)和损伤存在概率源图ⅱ；

步骤四、融合计算损伤概率值。将上述步骤四所得的损伤存在概率源值p′i-j(x,y)和p″i-j(x,y)按公式进行并集、求和、平均等计算，得到待检测点(x,y)处的损伤概率值p(x,y)；

步骤五、构件损伤概率成像定位。将所述步骤五所得p(x,y)值作为直角坐标系坐标点(x,y)的像素值。像素值越大的区域表示该处损伤存在的可能性越大，即为损伤的位置。

本发明的技术效果在于，采用超声表面波检测阵列，综合使用基于信号能量的损伤概率成像法和基于渡越时间的损伤概率成像法，只需少量的压电换能器，即可实现构件损伤的快速检测和高效定位。

附图说明

图1为本发明超声表面波损伤概率成像定位法流程图

图2为本发明带楔块的压电换能器产生表面波的示意图

图3为本发明压电换能器4和压电换能器7的直接接收模式信号s4-7和反射回波模式信号e4-7示意图

图4为本发明压电换能器4和压电换能器7直接接收模式损伤检测示意图

图5为本发明压电换能器4和压电换能器7反射回波模式损伤检测示意图

具体实施方式

下面结合附图对本发明作进一步说明。

本发明具体实施方式以铝板损伤的阵列超声表面波检测定位为例，铝板尺寸为1200mm×1000mm，待测成像区面积设为200mm×200mm的正方形,对待测区进行阵列超声表面波损伤概率成像的步骤包括：

步骤一、构建压电换能器检测阵列，如图3所示，在铝板表面正方形被测范围四周均匀布置16个带楔块的压电换能器，记为1,2,...,16，以正方形中心为原点建立平面直角坐标系xoy，并分别获取发射换能器中心与损伤边界、损伤边界与接收换能器中心、发射换能器中心与接收换能器中心之间的距离

步骤二、激励和采集超声表面波信号，采用两种激励方式激励压电换能器ai(xi,yi)(i＝1,2,...,n)产生超声信号，在构件表面产生超声表面波信号，其中一种激励方式采用tektronixafg2021任意信号发生器产生汉宁窗调制的正弦信号，并经rfpower2100l射频功率放大器放大后激励超声换能器，另一种激励方式采用olympus5072pr超声脉冲发生接收仪产生脉冲信号激励超声换能器，采集正弦信号激励方式的由其他换能器rj(xj,yj)(j＝1,2,...,n,j≠i)直接接收的超声表面波信号，记为si-j，采集脉冲信号激励方式的经损伤边界反射后由其他压电换能器rj(xj,yj)(j＝1,2,...,n,j≠i)接收的超声表面波信号，记为ei-j，如图4为压电换能器4作为发射器，压电换能器作为接收器的直接接收模式信号s4-7和反射回波模式信号e4-7；

步骤三、计算损伤存在概率源值，将所述步骤二中的直接接收模式信号si-j和反射回波模式信号ei-j，分别按以下(1)、(2)进行基于信号能量的损伤概率成像算法和基于渡越时间的损伤概率成像算法，计算检测点(x,y)处的损伤存在概率源值：

(1)对所述步骤二中的直接接收模式信号si-j进行基于信号能量的损伤概率成像算法，按以下步骤进行：

步骤1、如图4以压电换能器4和压电换能器7为例，信号s4-7经傅里叶变换后，获取构件无损伤时中心频率上阵列超声表面波的时间-能量波形分布曲线，作为基准信号e0，保存a4-s7压电换能器对的能量波形数据；

步骤2、重复上述步骤1，获取并保存构件有损伤时的a4-s7压电换能器的时间-能量波形数据，作为检测信号e'，然后根据公式计算检测信号e'与基准信号e0的相异系数q；

步骤3、设发射压电换能器坐标为(x4,y4)，接收压电换能器坐标为(x7,y7)，待测点坐标为(x,y)，待测点到两压电换能器连线的距离为d,定义距离系数δ＝0.05，距离函数根据上述步骤2所得结果，再根据公式i'4-7(x,y)＝q×ψ(x,y)，计算得到a4-s7检测路线的检测范围内的损伤存在概率值i'4-7(x,y)，式中

步骤4、对其他每一对ai-rj压电换能器重复进行上述步骤1-3的处理；

步骤5、定义强化系数β为1.02，将上述步骤4得到的损伤存在概率值i′i-j(x,y)根据对应检测点(x,y)按公式进行累加，重建直接接收模式下所有ai-rj检测路线确定的损伤存在概率源值p′i-j(x,y)和损伤存在概率源图ⅰ；

(2)对所述步骤二中的反射回波模式信号ei-j基于渡越时间的损伤概率成像算法，按以下步骤进行：

步骤1、如图5以压电换能器4和压电换能器7为例，设反射回波模式信号e4-7的激励时刻为t4，接收时刻为t7，从激励时刻起，经损伤边界反射，到接收时刻止，表面波信号经历的渡越时间为t″4-7，路程为统计所有ai-rj检测路线的反射回波模式信号峰值时刻t的方差σ，将反射回波模式信号时刻tj的瞬时值按正态分布处理，取接收回波信号98％能量部分对应的时间区间(t4-7,t'4-7)，将接收时间区间减去激励时刻，得到渡越时间区间(t4-7,t′4-7)，由公式绘得椭圆环检测区；

步骤2、由上述步骤1确定的不同渡越时间t″4-7值所占整个渡越时间区间(t4-7,t′4-7)的权重为a4-r7检测路线确定的定位椭圆环检测区内检测点(x,y)处的损伤存在概率值

步骤3、对其他换能器对ai-rj重复进行上述步骤1-2的处理；

步骤4、将上述步骤3所得的所有ai-rj检测路线确定的定位椭圆环检测区包含的检测点(x,y)的损伤存在概率值i″i-j(x,y)，对应到被测工件表面坐标系上，按公式进行累加计算，得到反射回波模式所有ai-rj检测路线的损伤存在概率源值p″i-j(x,y)和损伤存在概率源图ⅱ；

步骤四、融合计算损伤概率值。将上述步骤四所得的损伤存在概率源值p′i-j(x,y)和p″i-j(x,y)按公式进行并集、求和、平均等计算，得到待检测点(x,y)处的损伤概率值p(x,y)；

步骤五、构件损伤概率成像定位。将所述步骤五所得p(x,y)值作为直角坐标系坐标点(x,y)的像素值。像素值越大的区域表示该处损伤存在的可能性越大，即为损伤的位置。