CN112444562B - 采用对称结构响应信号分析的无基准导波损伤监测方法 - Google Patents
采用对称结构响应信号分析的无基准导波损伤监测方法 Download PDFInfo
- Publication number
- CN112444562B CN112444562B CN202011339159.8A CN202011339159A CN112444562B CN 112444562 B CN112444562 B CN 112444562B CN 202011339159 A CN202011339159 A CN 202011339159A CN 112444562 B CN112444562 B CN 112444562B
- Authority
- CN
- China
- Prior art keywords
- signal
- damage
- response signal
- signals
- wave
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/50—Processing the detected response signal, e.g. electronic circuits specially adapted therefor using auto-correlation techniques or cross-correlation techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/10—Number of transducers
- G01N2291/106—Number of transducers one or more transducer arrays
Landscapes
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种采用对称结构响应信号分析的无基准导波损伤监测方法,在板结构件上对称布置PZT压电阵列,阵列中的任意两个PZT压电元件分别作为激励器和传感器而组成一条监测路径;获得滤波后的结构响应信号,对滤波后的信号进行包络提取,定位响应信号中的各波包中心点位置;计算当前响应信号中第一个相似度低于阈值的波包中心点时刻,记录为当前响应信号的损伤散射信号到达时间;交换另一组响应信号为计算基准,记录该响应信号的损伤散射信号到达时间;获得所有对称路径下的损伤散射信号到达时刻,并计算出损伤位置。本发明不需事先采集健康状态下的结构响应信号作为基准,消除了对基准信号的依赖,降低了时变条件引起的监测误差,具有实用意义。
Description
技术领域
本发明涉及一种采用工程结构的无基准在线损伤监测方法,具体涉及一种采用对称结构响应信号分析的无基准导波损伤监测方法。
背景技术
超声导波结构损伤监测方法一直是结构健康监测研究和应用的热点,其优势在于传播/监测距离远、对小损伤敏感。近些年来,这一技术方法的研究取得了长足的发展和进步,在损伤定位、成像、评估等方面均有很多富有成效的研究成果。另一方面,导波作为一种特殊的超声波,其传播过程又存在一定的特殊性,即多模特性和频散特性。这些特殊的传播特性一定程度上限制了导波损伤监测研究和应用的发展。前者在于响应信号中会同时存在多个模态的信号,造成目标特征提取的困难,后者在于不同模态和不同频率的信号传播速度会不同,尤其在宽带信号激励下,信号的理解十分困难,给信号分析和处理带来新的问题。为了缓解上述难题,在大多数的研究中,一般采用窄带激励以降低频散效应的影响,同时,采用基准信号的方式,通过分析对比损伤前后结构响应信号的变化,来实现损伤的判别和定位等。然而,由于监测环境、工况、载荷等条件的差异性,此类采用基准信号的监测方法在应用过程中面临一些新的问题,例如温度差异会导致结构响应信号的变化,从而会对监测结果的准确性带来影响,甚至造成误判。
现有的无基准导波结构损伤监测主要从传感和信号处理两类方式入手,例如自传感技术和双元激励传感,响应信号中仅存在反射或散射信号,由于损伤也是散射源,因此可以进而获得损伤信息;基于时间反转聚焦的损伤因子方法是一种典型的无基准损伤监测方法,其思想在于损伤引起的监测信号传播非线性效应,导致时间反转聚焦重构差异,进而获取监测通道及其附近的损伤发生情况。上述方法一定程度上实现了导波无基准损伤监测,但也存在一定的局限性,例如自传感技术中,损伤作为散射源之一,其引起的信号散射和反射往往较为微弱,难易提取;基于时间反转聚焦损伤因子方法仅对监测通道敏感,因此应用中也会面临一定的限制。
本发明在现有基于压电阵列的导波损伤监测方法基础上,从板结构导波传播的过程出发,根据损伤对信号传播的作用机理,利用损伤引起的对称响应信号差异性实现无基准结构损伤监测,为导波结构健康监测发展和应用提供新的思路。
发明内容
发明目的:本发明提供一种采用对称结构响应信号分析的无基准导波损伤监测方法,发现并定位损伤引起的散射信号波包及其中心时刻点位置,并基于定位出的损伤散射信号波包中心点位置,实现对损伤的定位,从而实现无基准损伤监测和定位的方法,消除对基准响应信号的依赖性。
技术方案:本发明所述的一种采用对称结构响应信号分析的无基准导波损伤监测方法,包括以下步骤:
(1)在板结构件上按照等传播距离对称布置PZT压电阵列,阵列中的任意两个PZT压电元件分别作为激励器和传感器而组成一条监测路径;
(2)选定同一激励器激励下具有相同传播路径和条件的任意两条监测路径M1、M2,在板结构件中激励以单一模式为主的窄带导波结构响应信号,传感采集两条监测路径下的结构响应,并对采集到的结构响应信号进行滤波处理,消除直流偏置以及高频噪声信号,获得滤波后的结构响应信号R1、R2;
(3)对滤波后的信号进行包络提取,计算定位两组结构响应信号R1、R2中的各波包中心点位置P;
(4)对R1、R2两组信号,以R2响应信号为计算基准,在定位出的波包中心点位置,以设定脉宽为窗口,逐个对比两组信号的波包相似度,计算得到R1响应信号中第一个相似度低于阈值的波包中心点时刻,记录为R1响应信号中的损伤散射信号到达时间t1;交换以R1响应信号为计算基准,重复信号比对过程,获得R2响应信号中第一个相似度低于阈值的波包中心点时刻,记录为R2响应信号中的损伤散射信号到达时间t2;
(5)重复步骤(2)至(4),获得所有对称路径下的损伤散射信号到达时刻ti,i=1,2,3,…,根据根据ti,计算得到损伤的位置。
进一步地,步骤(2)所述的激励信号均为窄带信号。
进一步地,所述步骤(3)实现过程如下:
进一步地,步骤(4)所述的波包相似度是采用皮尔逊相关系数法,通过信号间差异性比较方法计算获得。
进一步地,所述步骤(5)的实现过程如下:
根据信号沿激励器-损伤-传感器的传播时间,确立信号传播时间和损伤位置对应方程,对于坐标为(xi,yi)的传感器以及其感知得到的损伤散射信号到达时间ti,可以获得关于损伤位置(x,y)的方程多个方程连列求解,即可计算得到损伤位置。
有益效果:与现有技术相比,本发明的有益效果:1、本发明不依赖于结构健康状态下或前期测得的结构响应信号为基准,消除了时变条件等对基于基准响应信号的损伤监测带来的干扰;2、本发明方法在实现过程中无需更改或增加设备和参数,利用现有硬件系统就可以实现;3、本发明实现方法简单,与现有方法相比,只需要测得当前状态下所有具有对称监测通道的结构响应信号,通过两两比对就可得出损伤信息。
附图说明
图1是本发明中实验试件及对称布置压电激励器/传感器阵列布局图;
图2是本发明中窄带激励信号的波形时域图;
图3是本发明中选取的两条对称路径下的两组结构响应信号时域波形图;
图4是本发明中选取的对称路径下的两组结构响应信号包络时域波形及其波包中心点位置图;
图5是本发明中选取的对称路径下两组响应信号中信号fd12相对于fd14的波包相似度计算结果图。
图6是本发明中选取的对称路径下两组响应信号中信号fd14相对于fd12的波包相似度计算结果。
具体实施方式
下面结合附图对本发明作进一步详细说明。
本发明提供一种采用对称结构响应信号分析的无基准导波损伤监测方法,包括如下步骤:
步骤1,在板结构件上按照等传播距离对称布置PZT压电阵列,阵列中的任意两个PZT压电元件分别作为激励器和传感器而组成一条监测路径。
在本实施例中,采用一块T6061铝板作为实验结构,其尺寸为600mm*600mm*6mm。为了简化条件,设计了4压电片组成正方形阵列居中布置,每个压电片距离最近的边界均为200mm,以板结构中心为坐标原点,各压电片坐标分别为:1#(100mm,100mm)、2#(-100mm,100mm)、3#(-100mm,-100mm)、4#(100mm,-100mm),具体如图1所示。实施例中采用粘贴质量块的方式来模拟损伤的发生情况,中心坐标为(50mm,10mm)。
步骤2,选定同一激励器激励下具有相同传播路径和条件的任意两条监测路径M1、M2,在板结构件中激励以单一模式为主的窄带导波结构响应信号,传感采集两条监测路径下的结构响应,并对采集到的结构响应信号进行滤波处理,消除直流偏置以及高频噪声信号,获得滤波后的结构响应信号R1、R2。
选定同一激励器1#压电片激励下2#、4#压电片传感,具有相同传播路径和条件的两条监测路径M12、M14,在结构中激励以单一模式为主的窄带导波结构响应信号,传感采集两条监测路径下的结构响应,并对采集到的结构响应信号进行滤波处理,消除直流偏置以及高频噪声信号,获得滤波后的结构响应信号fd12、fd14。
采用如图2所示的余弦调制五周期正弦信号为激励信号,激励信号经功率放大器输出放大后在压电激励器1上加载,同时采集经电荷放大器放大输出的压电传感器2和4的输出信号,并按照激励信号的带宽范围进行带通滤波,消除直流偏置和高频噪声,最终获得结构响应信号fd12、fd14,如图3所示。为了激发单一模态为主的导波响应信号,经过实验测定,当激励信号中心频率为225kHz时,激发出的导波结构响应主要是以S0模态为主,简化了信号的复杂度。
步骤3,对滤波后的信号进行包络提取,计算定位两组结构响应信号R1、R2中的各波包中心点位置P。
步骤4,对R1、R2两组信号,以R2响应信号为计算基准,在定位出的波包中心点位置,以设定脉宽为窗口,逐个对比两组信号的波包相似度,计算得到R1响应信号中第一个相似度低于阈值的波包中心点时刻,记录为R1响应信号中的损伤散射信号到达时间t1;交换以R1响应信号为计算基准,重复信号比对过程,获得R2响应信号中第一个相似度低于阈值的波包中心点时刻,记录为R2响应信号中的损伤散射信号到达时间t2。
对fd12、fd14两组信号,首先以fd14响应信号为计算基准,在定位出的波包中心点位置,以设定脉宽为窗口,逐个对比两组信号的波包相似度,如图5所示,计算得到fd12响应信号中第一个相似度低于阈值,的波包为响应信号的第3个波包,本发明实施例中设定为0.7,将该波包中心点时刻0.0997ms记录为fd12响应信号中的损伤散射信号到达时间t1;交换以fd12响应信号为计算基准,重复上述信号比对过程,如图6所示,获得fd14响应信号中第一个相似度低于阈值的波包为第2个波包,将该中心点时刻0.0958ms记录为fd14响应信号中的损伤散射信号到达时间t2。
其中,在计算信号相似度中,本发明实施例中采用了皮尔逊相关系数(PearsonCorrelation Coefficient)法进行相似度运算,该算法被广泛认为是一种行之有效的信号间差异性比较的方法,对于滤除直流量之后的信号x1(t)及x2(t),两组信号的相似度r计算基本原理如下所示:
由于存在传感器布置位置的轻微差异,对称响应信号中的波包并不是严格对齐,严格按照两组响应信号中的一组作为基准信号,根据其提取出的各波包中心点进行波包相似度计算,容易带来一定的偏差,因此,实际计算中,固定基准信号提取出的波包不变,另外一组信号的波包做适当延迟和超前(本发明实施例中选取±5个采样点),以计算出的相似度最大值为准。
步骤5,重复步骤2至4,获得所有对称路径下的损伤散射信号到达时刻ti,i=1,2,3,…,根据根据ti,计算得到损伤的位置.
根据信号沿激励器-损伤-传感器的传播时间,确立信号传播时间和损伤位置对应方程,对于坐标为(xi,yi)的传感器以及其感知得到的损伤散射信号到达时间ti,可以获得关于损伤位置(x,y)的方程:
其中,ν为导波波速,多个方程连列求解,即可计算得到损伤位。
根据fd12、fd14两组信号得到的损伤波包到达时间,在已知波速和激励信号开始时刻的前提下,可以直接列出与损伤坐标有关的信号传播时间与距离方程,通过求解方程可以计算得到损伤监测结果为(31mm,22mm),与真实损伤中心点的误差为(-19mm,12mm)。重复上述过程,3#压电片激励,2#、4#压电片传感得到的结构响应信号经过同样的数据处理过程,可以判别出这两组响应信号中的损伤波包中心点出现时刻分别为0.1056ms和0.1075ms,据此,计算得到的损伤监测结果为(29mm,26mm),误差为(-21mm,16mm)。在误差允许范围内,实现了损伤的无基准监测和定位。
Claims (3)
1.一种采用对称结构响应信号分析的无基准导波损伤监测方法,其特征在于,包括如下步骤:
(1)在板结构件上按照等传播距离对称布置PZT压电阵列,阵列中的任意两个PZT压电元件分别作为激励器和传感器而组成一条监测路径;
(2)选定同一激励器激励下具有相同传播路径和条件的任意两条监测路径M1、M2,在板结构件中激励以单一模式为主的窄带导波结构响应信号,传感采集两条监测路径下的结构响应,并对采集到的结构响应信号进行滤波处理,消除直流偏置以及高频噪声信号,获得滤波后的结构响应信号R1、R2;
(3)对滤波后的信号进行包络提取,计算定位两组结构响应信号R1、R2中的各波包中心点位置P;
(4)对R1、R2两组信号,以R2响应信号为计算基准,在定位出的波包中心点位置,以设定脉宽为窗口,逐个对比两组信号的波包相似度,计算得到R1响应信号中第一个相似度低于阈值的波包中心点时刻,记录为R1响应信号中的损伤散射信号到达时间t1;交换以R1响应信号为计算基准,重复信号比对过程,获得R2响应信号中第一个相似度低于阈值的波包中心点时刻,记录为R2响应信号中的损伤散射信号到达时间t2;
(5)重复步骤(2)至(4),获得所有对称路径下的损伤散射信号到达时刻ti,i=1,2,3,…,根据ti,计算得到损伤的位置;所述损伤的位置计算过程如下:
3.根据权利要求1所述的采用对称结构响应信号分析的无基准导波损伤监测方法,其特征在于,步骤(4)所述的波包相似度是采用皮尔逊相关系数法,通过信号间差异性比较方法计算获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011339159.8A CN112444562B (zh) | 2020-11-25 | 2020-11-25 | 采用对称结构响应信号分析的无基准导波损伤监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011339159.8A CN112444562B (zh) | 2020-11-25 | 2020-11-25 | 采用对称结构响应信号分析的无基准导波损伤监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112444562A CN112444562A (zh) | 2021-03-05 |
CN112444562B true CN112444562B (zh) | 2022-09-23 |
Family
ID=74738084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011339159.8A Active CN112444562B (zh) | 2020-11-25 | 2020-11-25 | 采用对称结构响应信号分析的无基准导波损伤监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112444562B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113155971A (zh) * | 2021-03-22 | 2021-07-23 | 南京信息职业技术学院 | 一种导波双点传感管道结构损伤检测方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101451977B (zh) * | 2008-12-30 | 2011-01-19 | 南京航空航天大学 | 一种基于双元压电传感阵列和时间窗函数的无基准Lamb波损伤监测方法 |
CN102445496B (zh) * | 2011-10-20 | 2013-05-08 | 南京航空航天大学 | 一种基于Lamb波的板状结构无基准快速损伤检测方法 |
CN102998370A (zh) * | 2012-12-14 | 2013-03-27 | 南京邮电大学 | 一种工程结构无基准Lamb波时间反转损伤监测方法 |
CN107045017B (zh) * | 2017-04-25 | 2019-07-02 | 南京信息工程大学 | 基于超声兰姆波和时间反转理论的薄板裂纹深度分析方法 |
CN113720917B (zh) * | 2021-07-12 | 2022-08-12 | 北京航空航天大学 | 一种准静态未知时变载荷下结构损伤的自适应定位方法 |
-
2020
- 2020-11-25 CN CN202011339159.8A patent/CN112444562B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112444562A (zh) | 2021-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110389170B (zh) | 基于Lamb波成像的列车部件裂纹损伤检测方法和系统 | |
Poddar et al. | Time reversibility of a Lamb wave for damage detection in a metallic plate | |
CN110687208A (zh) | 一种基于双曲线定位的无基准Lamb波损伤监测方法 | |
CN111208206B (zh) | 一种温度变化环境下的复合材料超声导波结构健康监测方法 | |
CN111380957B (zh) | 一种基于频率-波数域滤波的超声导波检测方法 | |
CN109696480B (zh) | 一种基于改进时间反转算法的玻璃纤维复合材料声发射源定位成像方法 | |
CN110068613B (zh) | 一种结构导波响应群速度频散测试方法 | |
CN113358743B (zh) | 一种基于时频分布相似度分析的兰姆波模态分离方法 | |
CN112444562B (zh) | 采用对称结构响应信号分析的无基准导波损伤监测方法 | |
CN110208383A (zh) | 一种基于反转路径差信号的板结构兰姆波无参考成像方法 | |
CN110243945A (zh) | 基于合成孔径聚焦与模式转换波的超声tofd盲区抑制方法 | |
CN104749249B (zh) | 一种基于超声波技术的检测种子净度的方法 | |
CN116698986B (zh) | 一种基于调谐带通滤波器的吸音材料质控系统和方法 | |
CN111076806B (zh) | 一种基于聚偏氟乙烯(pvdf)压电薄膜的结构健康监测装置及方法 | |
CN102520075A (zh) | 基于谐波小波包的直升机部件声发射信号传播特性的测试方法 | |
Kabeya III | Structural health monitoring using multiple piezoelectric sensors and actuators | |
KR20100090912A (ko) | 유도초음파를 이용한 구조건전성 모니터링 방법 | |
Janarthan et al. | Damage detection in stiffened composite panels using Lamb wave | |
CN108646248A (zh) | 一种针对低速运动声源的被动声学测速测距方法 | |
CN113155971A (zh) | 一种导波双点传感管道结构损伤检测方法 | |
CN110441391B (zh) | 一种高频Lamb波传播特性补偿方法 | |
CN110057911B (zh) | 一种声表面波无损检测系统 | |
CN115236192A (zh) | 一种声发射信号检测方法及系统 | |
CN109884193A (zh) | 一种基于有效信号段自提取算法和分段截取相关系数算法的椭圆定位算法 | |
CN110672717B (zh) | 检查结构的方法以及检查系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |