CN109990968A - 一种基于时间反转法原理的板结构冲击定位方法 - Google Patents

一种基于时间反转法原理的板结构冲击定位方法 Download PDF

Info

Publication number
CN109990968A
CN109990968A CN201910223725.XA CN201910223725A CN109990968A CN 109990968 A CN109990968 A CN 109990968A CN 201910223725 A CN201910223725 A CN 201910223725A CN 109990968 A CN109990968 A CN 109990968A
Authority
CN
China
Prior art keywords
impact
signal
sensor
calibration
reconstruct
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.)
Granted
Application number
CN201910223725.XA
Other languages
English (en)
Other versions
CN109990968B (zh
Inventor
陈春林
李名锐
周刚
马坤
初哲
王可慧
冯娜
赵南
王长利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwest Institute of Nuclear Technology
Original Assignee
Northwest Institute of Nuclear Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwest Institute of Nuclear Technology filed Critical Northwest Institute of Nuclear Technology
Priority to CN201910223725.XA priority Critical patent/CN109990968B/zh
Publication of CN109990968A publication Critical patent/CN109990968A/zh
Application granted granted Critical
Publication of CN109990968B publication Critical patent/CN109990968B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00—Vibration-testing of structures; Shock-testing of structures
    • G01M7/08—Shock-testing
    • 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/045—Analysing solids by imparting shocks to the workpiece and detecting the vibrations or the acoustic waves caused by the shocks

Abstract

本发明公开一种基于时间反转法原理的板结构冲击定位方法，其属于结构健康监测领域。包括以下步骤：步骤一：压电传感器布局；步骤二：获得板结构冲击响应传递函数；步骤三：对冲击事件的传感器的冲击响应信号进行时间反转数据处理；步骤四：提取重构冲击信号特征值参数；步骤五：利用特征值参数确定冲击位置。发明所述定位方法充分利用了传递函数包含的信息进行冲击定位，不需要已知结构中的波速、结构材料性质及力学边界条件，具有良好的稳定性，该定位方法中时间和空间上的信息可以互补，仅采用1‑2个传感器就能够实现板结构冲击定位，在航空航天结构的冲击监测方面具有广泛的应用前景。

Description

一种基于时间反转法原理的板结构冲击定位方法
技术领域
本发明属于结构健康监测的冲击监测技术领域，提出了一种基于时间反转法原理的板结构冲击定位方法。
背景技术
航空航天器结构在服役过程中不可避免地会受到外界冲击的威胁，如空间碎片、飞鸟、弹丸/片、冰雹和掉落的工具等。航空航天器结构多为板结构，冲击损伤会显著降低结构的力学性能并减少其使用寿命，尤其是复合材料结构。常规的无损检测技术，如射线检测、超声波检测、磁粉检测等，大多需在地面进行并对局部结构进行拆卸，需耗费较多的人力和时间。
通过传感器采集的冲击响应信号来识别冲击位置和载荷时间历程，结合结构已有的参数和实验数据，就能对结构的冲击损伤状态和剩余强度进行评估。结构健康监测系统能实时识别冲击结果，迅速为结构的安全评估和维修提供重要依据。
三角定位法是目前国内外最常用的冲击定位方法，该方法基于传感器几何关系和应力波传播时间差，其前提是波从激励源到传感器沿直线传播。三角定位法原理简单，但需要精确获得波传播到各传感器的时间差和波在结构中的传播速度。实际应用中测量信号总会受到噪音的干扰，使得波的传播时间和速度难以精确测定。复杂结构中的波会在不连续处(接头、加强筋等)发生反射和散射，这也阻碍了波传播时间的精确测量。此外，各向异性复合材料中波的传播速度还具有方向性，都使得普通的三角定位法难以成立。
目前的优化方法是先根据板结构的激励和响应之间的关系建立优化模型，然后把位置识别问题转化为优化求逆问题来定位。该方法需建立精确的模型，且在每个候选位置做多次运算，对于大尺度复杂结构，该方法可能由于计算过程过于复杂而不实用。
人工神经网络方法能用于解决任意结构上的冲击定位问题，但其训练样本数据通常较多。与神经网络类似的人工智能方法，如模糊逻辑方法、支持向量机、特征正交分法和免疫遗传算法等也会面临同样的问题。冲击过程中，冲击响应信号的幅值随传播距离的增大而降低，基于这一特性，可采用均匀分布传感器阵列采集信号的能量分布趋势进行冲击位置识别，此类方法原理简单，但大都需要数目较多的传感器。
综上，冲击识别存在的普遍问题是冲击定位和冲击载荷重构的独立进行需要较多的已知条件和复杂的计算过程。
发明内容
为了解决冲击定位和冲击载荷重构的独立进行需要较多已知条件和计算过程复杂的问题，本发明提供一种基于时间反转法原理的冲击定位方法，基于时间反转法原理，与其他识别方法相比，实验标定获得的传递函数能同时用于冲击定位和冲击载荷重构，冲击位置的识别不需要已知结构中的波速等信息，仅需要1-2个传感器就能进行冲击定位，抗噪性能好，方便可靠，实用性强，适用于不同材料不同形状的板结构冲击，如加筋板、复合材料板等复杂板结构冲击定位，可应用于航空航天器外结构冲击定位。
该方法采用压电传感器网络，获得板结构的冲击响应信号，实验标定获得的传递函数，冲击发生后，利用时间反转数据处理方法重构冲击响应信号，然后提取重构信号的特征值，最后利用重构信号特征值的空间分布特征实现冲击定位。
本发明的技术解决方案是提供一种基于时间反转法原理的板结构冲击定位方法，包括以下步骤：
步骤一、通过标定实验，获得待测板结构冲击响应传递函数Gi(xjk,ω)；
在待测板结构上布置n个传感器，通过标定实验获得待测板结构不同位置上的冲击响应传递函数Gi(xjk,ω)；其中，xjk为第j行第k列个标定点；
步骤二、根据Gi(xjk,ω)，对冲击事件的传感器的冲击响应信号si(t)进行时间反转数据处理，获得第i个传感器的冲击响应信号si(t)在标定点xjk处的重构冲击信号fi(xjk,t)；
2.1)、待测板结构发生冲击后，在频域中对冲击事件的传感器的冲击响应信号si(t)进行时间反转数据处理，获得频域中该传感器信号si(ω)在xjk处的重构冲击信号fi(xjk,ω)；
其中，si(ω)由si(t)进行傅里叶变换获得；上标(*)表示复数的共轭运算，对应于时域中的时间反转处理；
2.2)、将fi(xjk,ω)进行逆傅里叶变换，获得时域中，第i个传感器的冲击响应信号si(t)在标定点xjk处的重构冲击信号fi(xjk,t)；
步骤三、提取重构冲击信号fi(xjk，t)的特征值；
3.1)、单一传感器重构冲击信号特征值提取；单一传感器重构冲击信号特征值包括峰值Pi(xjk)和相对峰值
(1)重构冲击信号的峰值为：
T为时间反转法处理时间长度；
(2)重构冲击信号的相对峰值为：
T1为信号聚集时段，T2＝T-T1，为聚集后时段；
(3)重构冲击信号的相对峰值为：
Mean表示平均值计算；
(4)重构冲击信号的相对峰值为：
3.2)、传感器的冲击响应信号融合，包括以下三种融合方式：
(1)针对每个标定点，将步骤3.1)获得的所有传感器相对应的特征值相加，获得每个标定点处，融合后的重构冲击信号的峰值与相对峰值；
(2)针对每个标定点，先将步骤二获得的所有传感器的重构冲击信号相加，然后根据步骤3.1)的提取方法提取特征值，获得每个标定点处，融合后的重构冲击信号的峰值与相对峰值；
(3)针对每个标定点，先对每个传感器的重构冲击信号进行归一化处理，然后相加，最后根据步骤3.1)的提取方法提取出特征值，获得每个标定点处，融合后的重构冲击信号的峰值与相对峰值；
步骤四、根据步骤三获得的融合后的峰值与相对峰值，优选时反重构冲击信号特征值C(xjk)；
将步骤三获得的融合后的峰值与相对峰值相乘，确定用于冲击定位的时反重构信号特征值C(xjk)：
C(xjk)＝[Pl(xjk)]o[γhg(xjk)]q
o＝0,1；q＝0,1；l＝1,2,3；h＝1,2,3；g＝1,2,3；其中，o和q不能同时为0；
步骤五、利用步骤四获得的特征值确定冲击位置；
5.1、确定初始的冲击位置：
5.2、确定冲击区域：
冲击区域由以初始的冲击位置为中心的m个标定点组成，m为正整数；
5.3、确定最终的冲击位置
C(xr，yr)为冲击区域内标定点处时反重构冲击信号特征值C(xjk)；
进一步地，步骤一中，n个压电传感器布置在待测板结构正面或背面的非检测区域。
进一步地，所述标定实验包括：
1)、在待测板结构的监测区域内标记多个标定点xjk，多个标记点均匀分布；
2)、采用脉冲力锤对标定点xjk施加低速冲击载荷，xjk为第j行第k列个标定点，记录标定载荷F(xjk，t)和n个压电传感器的冲击响应信号Si(t)，其中i＝1，2，...，n为传感器编号，t为时间；
3)、通过数据处理获得不同标定点与不同传感器之间的冲击响应传递函数，频域中冲击响应传递函数定义为：
Gi(xjk,ω)＝Si(ω)/F(xjk,ω)
F(xjk,ω)为点xjk处的频域中的标定载荷，Si(ω)为该标定载荷下第i个传感器的冲击响应，F(xjk,ω)、Si(ω)分别由F(xjk，t)和Si(t)进行傅里叶变换获得。
进一步地，所述传感器为压电传感器。
进一步地，n等于1或2。
本发明的有益效果是：
(1)本发明提出采用时间反转法原理对冲击信号进行重构，然后提取重构冲击信号的特征值进行定位，该方法通过测量不同位置处冲击载荷引起的冲击响应，实现对冲击载荷位置的辨识，并不需要精确测量应变值。
(2)本发明利用了时间反转法的自聚焦特性，采用重构冲击信号的特征值进行冲击载荷定位，该方法抗噪性能好、简单可靠。
(3)本发明适用于不同材料(金属板、复合材料板)和不同形状(加筋板)板结构的冲击定位，不需要知道板结构的波速和边界条件。
(4)本发明提出的方法不仅仅采用重构冲击信号峰值进行位置识别，还提出了相对峰值的概念和传感器的冲击响应信号融合方法，用于位置识别的特征值体系丰富，拓展了时间反转法的应用。
(5)本发明的冲击定位方法与冲击载荷识别关系紧密，标定实验获得的传递函数不仅能够用于冲击位置识别，也能用于冲击载荷时间曲线重构。
(6)本发明采用压电传感器构成网络，不需要电源和信号调理，成本较低；
(7)本发明仅采用1-2个传感器就能够进行冲击载荷定位，对传感器的数目要求低。
(8)本发明采用压电传感器，但不限于此类传感器，也可以采用其它类型用于应变测量的传感器，如PVDF压电薄膜、电阻应变片、光纤应变计等传感器。
(9)本发明中用于冲击定位的时间上(信号处理时长，聚集时段，峰值和相对峰值)和空间上(标定间距，传感器数目等)的信息具有互补性，适用性强。
(10)本发明中不同传感器的冲击响应信号的识别结果可以互相印证，能够识别传感器的状态是否正常。
附图说明
图1为本发明方法流程图；
图2为板结构示意图；
图3为重构冲击信号的相对峰值与时间关系图；
图4为冲击位置识别结果示意图。
具体实施方式
以下结合附图及具体实施例对本发明做进一步地描述。
如图1所示，本发明通过以下步骤实现板结构冲击定位：
步骤一、压电传感器布局；
压电传感器用于测量板结构的应变响应，压电传感器布置在板结构正面或背面，理论上传感器可以布置在任意位置，为避免传感器被冲击载荷损坏，将传感器布置在监测区域四周，依次记作PZT1、PZT2、…、PZTn；
步骤二、通过标定实验获得板结构不同位置上的冲击响应传递函数；
2-1、如图2所示，在板结构监测区域内标记标定点，各标定点均匀分布，标定点的间距为D；
2-2、采用脉冲力锤对板结构监测区域内均匀分布的标定点xjk施加低速冲击载荷，xjk为第j行第k列个标定点。记录标定载荷F(xjk，t)和n个压电传感器的冲击响应信号Si(t)，其中i＝1，2，...，n为传感器编号；
2-3、通过数据处理获得不同标定点与不同压电传感器之间的冲击响应传递函数，频域中冲击响应传递函数定义为
Gi(xjk,ω)＝Si(ω)/F(xjk,ω) (1)
F(xjk,ω)为频域中点xjk处的标定载荷，Si(ω)为频域中该标定载荷下第i个传感器的冲击响应，分别由F(xjk，t)和Si(t)进行傅里叶变换获得；
步骤三、对冲击事件压电传感器的冲击响应信号进行时间反转数据处理；获得fi(xjk，t)。
板结构发生冲击后，第i个压电传感器记录的冲击响应信号为si(t)。在频域中时间反转数据处理过程表示为
fi(xjk，ω)为由应变信号si(ω)在xjk处重构的冲击信号，上标(*)表示复数的共轭运算，它对应于时域中的时间反转处理。将fi(xjk，ω)进行逆傅里叶变换，获得fi(xjk，t)。
步骤四、提取重构冲击信号特征值参数；
包括选取单一传感器重构冲击信号的特征值参数和选用传感器的冲击响应信号融合方式。单一传感器重构冲击信号的特征值包括峰值和相对峰值。
4-1、选取重构冲击信号特征值参数的类型，单一传感器重构冲击信号特征值包括峰值Pi(xjk)和相对峰值
(1)重构冲击信号的峰值为：
fi(xjk,t)为由第i个传感器采集的信号在标定点xjk处重构的冲击信号，T为时间反转法处理时间长度。
(2)重构冲击信号的相对峰值为：
T1为信号聚集时段，T2＝T-T1为聚集后时段，如图3所示。
(3)重构冲击信号的相对峰值为：
Mean表示平均值计算。
(4)重构冲击信号的相对峰值为：
4-2、传感器的冲击响应信号融合方式包括以下三种：
(1)传感器的冲击响应信号融合方式R1，先提取每个传感器重构冲击信号的特征值Pi(xjk)(或γi 1(xjk)，γi 2(xjk)，γi 3(xjk))，然后将特征值相加，融合后的特征值P1(xjk)(或γ11(xjk)，γ21(xjk)，和γ31(xjk))表示为：
(2)传感器的冲击响应信号融合方式R2，先将每个传感器重构的冲击信号相加然后按照单一传感器重构冲击信号特征值提取的过程提取特征值，获得每个标定点处，融合后的重构冲击信号峰值与相对峰值。
(3)传感器的冲击响应信号融合方式R3，先对每个传感器重构冲击信号进行归一化处理然后相加最后提取出特征值，融合后的特征值表示为：
传感器的数目为1时，上述三种融合方式的结果是一致的。相对峰值的形式和传感器冲击响应信号融合的方式都有3种，多传感器融合后的相对峰值有9种情况。多传感器融合后的峰值有3种情况。
4-3、优选时反重构冲击信号特征值
重构冲击信号的峰值和相对峰值都能体现时间反转法的聚焦特性，同时利用重构冲击信号的这两个特征能够进一步提高位置识别的结果，融合后的峰值和相对峰值的组合采用相乘的方式。
时反重构冲击信号特征值C(xjk)可以为：
C(xjk)＝[Pl(xjk)]o[γhg(xjk)]q (10)
o＝0，1；q＝0，1；l＝1，2，3；h＝1，2，3；g＝1，2，3，其中o和q不能同时选0。
在具体结构中，可以采用不同信号融合方式和特征值进行位置识别，然后优选出合适的多重特征值。
步骤五、利用特征值参数确定冲击位置
5-1、确定初始的冲击位置
重构冲击信号在离冲击位置最近的标定点上有最大特征值，由此确定初始的冲击位置表示为
5-2、确定冲击区域
冲击区域由以初步识别位置为中心的a×b(a＝2，3；b＝2，3)个标定点组成，如图2所示，图中I1，I2和I3表示不同的初步识别位置。
5-3、确定最终的冲击位置
最终识别的位置由冲击区域的9(或6，或4)个标定点处重构冲击信号特征值的重心确定，计算公式如下：
C(xr，yr)为标定点处重构冲击信号的特征值，m＝9(或6，4)。
冲击位置识别结果如图4所示。

Claims (5)

1.一种基于时间反转法原理的板结构冲击定位方法，其特征在于，包括以下步骤：
步骤一、通过标定实验，获得待测板结构冲击响应传递函数Gi(xjk,ω)；
在待测板结构上布置n个传感器，通过标定实验获得待测板结构不同位置上的冲击响应传递函数Gi(xjk,ω)；其中，xjk为第j行第k列个标定点；
步骤二、根据Gi(xjk,ω)，对冲击事件的传感器的冲击响应信号si(t)进行时间反转数据处理，获得第i个传感器的冲击响应信号si(t)在标定点xjk处的重构冲击信号fi(xjk,t)；
2.1)、待测板结构发生冲击后，在频域中对冲击事件的传感器的冲击响应信号si(t)进行时间反转数据处理，获得频域中该传感器信号si(ω)在xjk处的重构冲击信号fi(xjk,ω)；
其中，si(ω)由si(t)进行傅里叶变换获得；上标(*)表示复数的共轭运算，对应于时域中的时间反转处理；
2.2)、将fi(xjk,ω)进行逆傅里叶变换，获得时域中，第i个传感器的冲击响应信号si(t)在标定点xjk处的重构冲击信号fi(xjk,t)；
步骤三、提取重构冲击信号fi(xjk，t)的特征值；
3.1)、单一传感器重构冲击信号特征值提取；单一传感器重构冲击信号特征值包括峰值Pi(xjk)和相对峰值
(1)重构冲击信号的峰值为：
T为时间反转法处理时间长度；
(2)重构冲击信号的相对峰值为：
T1为信号聚集时段，T2＝T-T1，为聚集后时段；
(3)重构冲击信号的相对峰值为：
Mean表示平均值计算；
(4)重构冲击信号的相对峰值为：
3.2)、传感器的冲击响应信号融合，包括以下三种融合方式：
(1)针对每个标定点，将步骤3.1)获得的所有传感器相对应的特征值相加，获得每个标定点处，融合后的重构冲击信号的峰值与相对峰值；
(2)针对每个标定点，先将步骤二获得的所有传感器的重构冲击信号相加，然后根据步骤3.1)的提取方法提取特征值，获得每个标定点处，融合后的重构冲击信号的峰值与相对峰值；
(3)针对每个标定点，先对每个传感器的重构冲击信号进行归一化处理，然后相加，最后根据步骤3.1)的提取方法提取出特征值，获得每个标定点处，融合后的重构冲击信号的峰值与相对峰值；
步骤四、根据步骤三获得的融合后的峰值与相对峰值，优选时反重构冲击信号特征值C(xjk)；
将步骤三获得的融合后的峰值与相对峰值相乘，确定用于冲击定位的时反重构信号特征值C(xjk)：
C(xjk)＝[Pl(xjk)]o[γhg(xjk)]q
o＝0,1；q＝0,1；l＝1,2,3；h＝1,2,3；g＝1,2,3；其中，o和q不能同时为0；
步骤五、利用步骤四获得的特征值确定冲击位置；
5.1、确定初始的冲击位置：
5.2、确定冲击区域：
冲击区域由以初始的冲击位置为中心的m个标定点组成，m为正整数；
5.3、确定最终的冲击位置
C(xr，yr)为冲击区域内标定点处时反重构冲击信号特征值C(xjk)。
2.根据权利要求1所述的基于时间反转法原理的板结构冲击定位方法，其特征在于，步骤一中，n个压电传感器布置在待测板结构正面或背面的非检测区域。
3.根据权利要求1所述的基于时间反转法原理的板结构冲击定位方法，其特征在于，所述标定实验包括：
1)、在待测板结构的监测区域内标记多个标定点xjk，多个标记点均匀分布；
2)、采用脉冲力锤对标定点xjk施加低速冲击载荷，xjk为第j行第k列个标定点，记录标定载荷F(xjk，t)和n个压电传感器的冲击响应信号Si(t)，其中i＝1，2，...，n为传感器编号，t为时间；
3)、通过数据处理获得不同标定点与不同传感器之间的冲击响应传递函数，频域中冲击响应传递函数定义为：
Gi(xjk,ω)＝Si(ω)/F(xjk,ω)
F(xjk,ω)为点xjk处的频域中的标定载荷，Si(ω)为该标定载荷下第i个传感器的冲击响应，F(xjk,ω)、Si(ω)分别由F(xjk，t)和Si(t)进行傅里叶变换获得。
4.根据权利要求1-3任一所述的基于时间反转法原理的板结构冲击定位方法，其特征在于：所述传感器为压电传感器。
5.根据权利要求4所述的基于时间反转法原理的板结构冲击定位方法，其特征在于：n等于1或2。
CN201910223725.XA 2019-03-22 2019-03-22 一种基于时间反转法原理的板结构冲击定位方法 Active CN109990968B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910223725.XA CN109990968B (zh) 2019-03-22 2019-03-22 一种基于时间反转法原理的板结构冲击定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910223725.XA CN109990968B (zh) 2019-03-22 2019-03-22 一种基于时间反转法原理的板结构冲击定位方法

Publications (2)

Publication Number Publication Date
CN109990968A true CN109990968A (zh) 2019-07-09
CN109990968B CN109990968B (zh) 2021-01-05

Family

ID=67130874

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910223725.XA Active CN109990968B (zh) 2019-03-22 2019-03-22 一种基于时间反转法原理的板结构冲击定位方法

Country Status (1)

Country Link
CN (1) CN109990968B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112729742A (zh) * 2020-12-25 2021-04-30 西安理工大学 一种基于概率成像的冲击定位方法
CN113866023A (zh) * 2021-08-27 2021-12-31 北京工业大学 一种预测岩石杆中应力波大小的方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL107159A (en) * 1992-10-02 1996-05-14 Univ Paris Vii Method and device for acoustic examination
US6161434A (en) * 1996-06-13 2000-12-19 Fink; Mathias Method and device for detecting and locating a reflecting sound source
CN101986756A (zh) * 2010-10-25 2011-03-16 电子科技大学 一种基于时间反转信号的无线定位方案
CN102128881A (zh) * 2010-12-22 2011-07-20 南京邮电大学 利用信号分解的Lamb波工程结构损伤监测方法
CN102998370A (zh) * 2012-12-14 2013-03-27 南京邮电大学 一种工程结构无基准Lamb波时间反转损伤监测方法
CN103389339A (zh) * 2013-07-22 2013-11-13 北京工业大学 一种基于时间反转周向Lamb波的大口径厚壁管道缺陷定位方法
CN103941232A (zh) * 2014-04-26 2014-07-23 南昌航空大学 一种基于时间反转聚焦技术的声发射源定位方法
CN104181235A (zh) * 2014-04-14 2014-12-03 中国商用飞机有限责任公司北京民用飞机技术研究中心 一种基于虚拟时间反转的声发射被动成像方法和装置
CN104698086A (zh) * 2015-03-13 2015-06-10 哈尔滨工程大学 一种基于时间反转及多径效应的一维构件应力波无损探伤方法
CN105954358A (zh) * 2016-05-04 2016-09-21 中国计量大学 一种TR与Duffing系统相结合的超声导波小缺陷定位检测方法
CN109374748A (zh) * 2018-10-30 2019-02-22 上海市特种设备监督检验技术研究院 声发射源时反定位成像方法
EP3477336A1 (en) * 2017-10-31 2019-05-01 BAE SYSTEMS plc Improvements in and relating to sonar

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL107159A (en) * 1992-10-02 1996-05-14 Univ Paris Vii Method and device for acoustic examination
US6161434A (en) * 1996-06-13 2000-12-19 Fink; Mathias Method and device for detecting and locating a reflecting sound source
CN101986756A (zh) * 2010-10-25 2011-03-16 电子科技大学 一种基于时间反转信号的无线定位方案
CN102128881A (zh) * 2010-12-22 2011-07-20 南京邮电大学 利用信号分解的Lamb波工程结构损伤监测方法
CN102998370A (zh) * 2012-12-14 2013-03-27 南京邮电大学 一种工程结构无基准Lamb波时间反转损伤监测方法
CN103389339A (zh) * 2013-07-22 2013-11-13 北京工业大学 一种基于时间反转周向Lamb波的大口径厚壁管道缺陷定位方法
CN104181235A (zh) * 2014-04-14 2014-12-03 中国商用飞机有限责任公司北京民用飞机技术研究中心 一种基于虚拟时间反转的声发射被动成像方法和装置
CN103941232A (zh) * 2014-04-26 2014-07-23 南昌航空大学 一种基于时间反转聚焦技术的声发射源定位方法
CN104698086A (zh) * 2015-03-13 2015-06-10 哈尔滨工程大学 一种基于时间反转及多径效应的一维构件应力波无损探伤方法
CN105954358A (zh) * 2016-05-04 2016-09-21 中国计量大学 一种TR与Duffing系统相结合的超声导波小缺陷定位检测方法
EP3477336A1 (en) * 2017-10-31 2019-05-01 BAE SYSTEMS plc Improvements in and relating to sonar
CN109374748A (zh) * 2018-10-30 2019-02-22 上海市特种设备监督检验技术研究院 声发射源时反定位成像方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112729742A (zh) * 2020-12-25 2021-04-30 西安理工大学 一种基于概率成像的冲击定位方法
CN112729742B (zh) * 2020-12-25 2023-08-22 西安理工大学 一种基于概率成像的冲击定位方法
CN113866023A (zh) * 2021-08-27 2021-12-31 北京工业大学 一种预测岩石杆中应力波大小的方法
CN113866023B (zh) * 2021-08-27 2023-11-10 北京工业大学 一种预测岩石杆中应力波大小的方法

Also Published As

Publication number Publication date
CN109990968B (zh) 2021-01-05

Similar Documents

Publication Publication Date Title
CN101865758B (zh) 基于多重信号分类算法的冲击载荷定位方法
Shrestha et al. Impact localization on composite structure using FBG sensors and novel impact localization technique based on error outliers
CN104142195B (zh) 基于超声波法的钢结构构件内部初始应力检测装置和方法
Shrestha et al. Low velocity impact localization on composite wing structure using error outlier based algorithm and FBG sensors
CN102422154B (zh) 一种结构损伤检测系统、设备以及结构损伤检测方法
Wang et al. Structural damage detection using cross correlation functions of vibration response
CN105510444B (zh) 基于超声导波的双向时间反演损伤成像方法
CN107014668A (zh) 一种基于压电和智能涂层传感器的疲劳裂纹综合监测方法
Saravanan et al. Damage detection in structural element through propagating waves using radially weighted and factored RMS
Amezquita-Sanchez et al. Nonlinear measurements for feature extraction in structural health monitoring
CN103438794B (zh) 一种基于压电传感网络的工程结构冲击位置确定方法
Al-Jumaili et al. Characterisation of fatigue damage in composites using an Acoustic Emission Parameter Correction Technique
CN109990968A (zh) 一种基于时间反转法原理的板结构冲击定位方法
CN108051126B (zh) 一种变厚度复合材料层合板低速冲击能量识别系统及其工作方法
Jang et al. Detection of impact damage in composite structures using high speed FBG interrogator
CN104316167B (zh) 一种低信噪比下结构响应异常检测方法
CN112285651B (zh) 基于特征提取矩阵相似度识别复合材料冲击源位置的方法及系统
CN107628268A (zh) 基于低频逼近系数变化幅度的单边固支机翼冲击定位方法
CN105445724A (zh) 单水听器自由场被动测距方法
Kessler et al. Hybrid passive/active impact detection & localization for aerospace structures
CN106482638A (zh) 基于全频段信号幅值能量和反函数求解的冲击判位方法
CN116136441A (zh) 基于非线性Lamb波的螺栓早期松动预紧状态监测方法、系统及装置
Kessler et al. Hybrid coherent/incoherent beam forming diagnostic approach to naval assets
CN109283071A (zh) 一种cfrp低速冲击损伤样本低试验代价获取方法
Maseras-Gutierrez et al. Detection of impacts in composite materials using piezoceramic sensors and neural networks

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