CN103698748A - 复合材料中基于信号波速与衰减补偿的2d-music冲击定位方法 - Google Patents
复合材料中基于信号波速与衰减补偿的2d-music冲击定位方法 Download PDFInfo
- Publication number
- CN103698748A CN103698748A CN201310685400.6A CN201310685400A CN103698748A CN 103698748 A CN103698748 A CN 103698748A CN 201310685400 A CN201310685400 A CN 201310685400A CN 103698748 A CN103698748 A CN 103698748A
- Authority
- CN
- China
- Prior art keywords
- signal
- piezoelectric sensor
- wave
- impact
- velocity
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
- G01S5/22—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Description
技术领域
本发明涉及一种复合材料中基于信号波速与衰减补偿的2D-MUSIC(二维多重信号分类)冲击定位方法,属于工程复合材料结构健康监测方法技术领域。
背景技术
复合材料相比于传统材料具有高比模量、比强度等独特优势,可以有效地减轻飞机的重量并提高飞机的安全性和存活力,在航空领域具有良好的应用前景。然而复合材料对低速冲击很敏感,很容易造成复合材料的凹陷、脱层、纤维断裂等微小损伤,导致复合材料性能下降。因此,冲击监测成为当前结构健康监测中一个很重要的课题。
由于Lamb波具有长距离的传播能力,基于Lamb波的结构健康监测方法具有很好的应用前景。现有基于Lamb波冲击监测技术的研究算法主要有基于时差定位技术方法,基于人工智能的优化方法以及基于系统建模的方法。
近几年,许多学者提出了阵列信号处理方法,阵列信号处理技术属于信号处理的一个重要分支,即通过多个传感器同时接收(或激励)信号,相对于单个传感器,有效地提高了信号的信噪比。在这些方法中,由Schmidt R.O.提出的多重信号分类(MUSIC)算法具有代表性算法之一,其算法的基本思想是将任意阵列输出数据的协方差矩阵进行特征值分解,从而得到与信号分量相对应的信号子空间和与信号分量相正交的噪声子空间,然后利用这两个子空间的正交性来估计信号的参数。正是由于MUSCI算法在特定条件下具有超高的分辨力、估计精度及稳定性,从而吸引了大量的学者对其进行深入的研究和分析,目前算法主要应用于声波、电磁波等领域。目前已有学者考虑将远场1D-MUSCI算法运用于结构健康监测中,只实现冲击源的角度估计。在远场1D-MUSCI算法中,通常假设信号源足够远,这时信号的波阵面可以近似假设为平面波,减小信号模型的复杂程度。然而当冲击发生在传感器阵列附近时,此时远场模型将不再适用,通常在远场情况下可以忽略信号源的波前曲率,即近似为平面波,而在近场区域,信号源波前的曲率不能忽略,即波前的需要用球面波来描述,此时信号源的位置信息须由方位角 、距离两个参数来确定,此时将算法扩展成二维参数搜索。由冲击产生的Lamb波向传感器阵列各个阵元传播时,路径和距离都不一致,即信号的波速和衰减将影响冲击的定位。在常规的方法中,为了简化计算,通常把Lamb波的传播速度设定为一个常数,把Lamb波信号幅值衰减简单地假设为随距离线性衰减的关系。
然而,复合材料具有各向异性的特点,Lamb波沿结构各个方向传播时具有不同的波速。而且在复合材料结构中的主要部位,常有加筋、变厚度等复杂的结构形式,这些都将加强了复合材料的各向异性,增加冲击精确定位的难度。此时,Lamb波在复合材料结构传播时信号幅值衰减不单单于传播距离有关。另外在沿着平面扫描时,导向矢量同样与Lamb波的传播速度有关。在运用MUSIC算法时,应该考虑波速变化、信号幅值衰减的影响。常规方法中波速与信号衰减假设将不再适用于复合材料近场冲击源定位。
因此,为了保证复合材料中冲击的定位精度,需要考虑对Lamb信号的波速和衰减进行补偿,建立更加符合复合材料结构的冲击定位方法。
发明内容
本发明针对常规MUSIC算法的不足,依据Lamb波在复合材料中的各向异性与信号衰减规律,提出了一种复合材料中基于信号波速与幅值衰减补偿的2D-MUSIC冲击定位方法。
本发明为解决其技术问题采用如下技术方案:
(1)在被监测结构的相应位置布置压电传感器网络,组成不同的激励-传感通道;
(3)当触发通道中信号幅值大于触发阈值时,各个压电传感器开始采集信号,将这些信号合成为矩阵形式的冲击阵列信号;
(5)测量各个激励-传感通道所对应角度的波速,建立信号波速曲线,得到第个压电传感器相对于第个压电传感器的时间延迟关系,其中表示冲击源信号到达第个阵元时相对于第个阵元的时间延迟,, 为压电传感器编号,为冲击源到达第个压电传感器的距离,为冲击源传播到第个压电传感器直线路径与轴的夹角, 为冲击源到达第个压电传感器的距离,为冲击源传播到第个压电传感器直线路径与轴的夹角;为冲击源传播到第个压电传感器的波速,为冲击源传播到第个压电传感器的波速;
(6)测量各个激励-传感通道所对应距离的信号幅值,并用指数函数建立信号幅值衰减曲线,其中为第个压电传感器的信号幅值,为第个压电传感器的信号幅值,为幅值衰减系数,冲击源到第个压电传感器与到第个压电传感器的距离差,第个压电传感器相对于第个压电传感器的相位变化;
(7)将步骤()中时间关系和步骤()中的幅值衰减关系代入阵列导向矢量中,其中为冲击源到达坐标原点的距离,为冲击源与轴的夹角,为冲击窄带信号的角频率,得到信号波速与衰减补偿的冲击定位算法的阵列导向矢量和补偿后的空间谱估计公式,其中为空间谱,为阵列导向矢量的转置,为为噪声子空间的转置;
A、采集不同激励—传感通道上的激励信号和传感信号,根据已知激励源到压电传感器距离和激励信号与传感信号的时间差,得到相应传播角度上的波速;
B、将步骤A中得到的波速值,通过曲线拟合得到波速曲线。
a、选择同一个压电传感器用相同激励方式进行激励,再选择不同压电传感器接收信号,采集这些激励—传感通道上的激励信号和传感信号,获取相应传播距离的信号幅值;
本发明的有益效果如下:
本方法不需要复杂的系统建模和提取信号波达时间,对信号幅值衰减和波速进行了补偿,减少了复杂复合材料的各向异性对冲击定位的影响,较好地提高了复合材料中冲击源的定位精度和实时性,能有效地应用于真实的航空复杂结构的冲击载荷定位中。
附图说明
图1 航空复合材料油箱结构冲击定位示意图,其中1、模拟冲击;2、航空油箱结构;3、传感器均匀线性阵列。
图2 冲击阵列信号波形图。
图3 PZT1冲击信号的频谱图。
图4 结构监测区域压电传感器布置情况图。
图5 Lamb波在航空油箱结构中的波速曲线。
图6 Lamb波在航空油箱结构中的信号幅值衰减曲线。
图7 冲击信号传播模型图。
图8 冲击源(131 mm, 108o)的定位结果图。
图9 复合材料中基于信号波速与衰减补偿的2D-MUSIC冲击定位方法流程图。
具体实施方式
下面结合附图对本发明创造做进一步详细说明。
如图所示,油箱的尺寸为, 油箱的上、下表面均为碳纤维复合材料结构,其余个侧面为铝合金板。以铆钉的方式固定,个侧面各有加筋板,以上表面碳纤维复合材料板作为冲击监测区域,其铺层情况为:每层厚度为,铺层顺序为, 中间最厚区域厚度为,两端最薄区域厚度为。结构监测区域表面布置个压电传感器(),传感器的编号如图所示。
以冲击源的冲击分析,冲击阵列信号如图所示,图给出了冲击响应的时域信号和频谱图,可以看出冲击响应较复杂,为宽带信号。为了满足算法的窄带前提条件,在运用该算法之前必须将宽带信号转换为窄带信号。小波变换在时域和频域都具有良好的局部分析特性,适合冲击信号的分析。因此,选用小波变换提取冲击响应中某一中心频率窄带信号。
.波速曲线测量
采集各个激励-传感通道的信号,得到传感信号的最大幅值,同理,条激励-传感通道,得到波在航空油箱结构中经过不同传播距离上的信号幅值。将个波速离散值,通过函数拟合得到幅值随距离的变化曲线,如图6所示,得到信号幅值衰减系数和。
另外各阵元接收信号的幅度差异也非常明显,由于幅度差中携带了关于信号源方位的重要信息,因此在冲击信号传播模型中应该考虑幅度差。这里将信号幅值衰减定义为指数函数衰减形式,即
写成矩阵的形式为
Claims (3)
(1)在被监测结构的相应位置布置压电传感器网络,组成不同的激励-传感通道;
(3)当触发通道中信号幅值大于触发阈值时,各个压电传感器开始采集信号,将这些信号合成为矩阵形式的冲击阵列信号;
(5)测量各个激励-传感通道所对应角度的波速,建立信号波速曲线,得到第个压电传感器相对于第个压电传感器的时间延迟关系,其中表示冲击源信号到达第个阵元时相对于第个阵元的时间延迟,, 为压电传感器编号,为冲击源到达第个压电传感器的距离,为冲击源传播到第个压电传感器直线路径与轴的夹角, 为冲击源到达第个压电传感器的距离,为冲击源传播到第个压电传感器直线路径与轴的夹角;为冲击源传播到第个压电传感器的波速,为冲击源传播到第个压电传感器的波速;
(6)测量各个激励-传感通道所对应距离的信号幅值,并用指数函数建立信号幅值衰减曲线,其中为第个压电传感器的信号幅值,为第个压电传感器的信号幅值,为幅值衰减系数,冲击源到第个压电传感器与到第个压电传感器的距离差,第个压电传感器相对于第个压电传感器的相位变化;
(7)将步骤()中时间延迟关系和步骤()中的幅值衰减曲线代入阵列导向矢量中,其中为冲击源到达坐标原点的距离,为冲击源与轴的夹角,为冲击窄带信号的角频率,得到信号波速与衰减补偿的冲击定位算法的阵列导向矢量和补偿后的空间谱估计公式,其中为空间谱,为阵列导向矢量的转置,为为噪声子空间的转置;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310685400.6A CN103698748B (zh) | 2013-12-16 | 2013-12-16 | 复合材料中基于信号波速与衰减补偿的2d‑music冲击定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310685400.6A CN103698748B (zh) | 2013-12-16 | 2013-12-16 | 复合材料中基于信号波速与衰减补偿的2d‑music冲击定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103698748A true CN103698748A (zh) | 2014-04-02 |
CN103698748B CN103698748B (zh) | 2017-06-27 |
Family
ID=50360342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310685400.6A Active CN103698748B (zh) | 2013-12-16 | 2013-12-16 | 复合材料中基于信号波速与衰减补偿的2d‑music冲击定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103698748B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016037400A1 (zh) * | 2014-09-09 | 2016-03-17 | 南京航空航天大学 | 基于能量加权因子的复合材料结构冲击区域定位方法 |
CN109115894A (zh) * | 2018-09-04 | 2019-01-01 | 温州大学激光与光电智能制造研究院 | 一种基于Toeplitz近似法的复合材料板损伤检测方法 |
CN109374740A (zh) * | 2018-09-21 | 2019-02-22 | 南京航空航天大学 | 基于阵列误差校正的合成孔径music损伤定位方法 |
CN109406627A (zh) * | 2018-09-11 | 2019-03-01 | 温州大学苍南研究院 | 一种基于旋转不变子空间的复合材料板冲击定位方法 |
CN109406084A (zh) * | 2018-09-13 | 2019-03-01 | 温州大学 | 一种对复合材料产品实现全频域快速冲击监测的方法 |
CN109696480A (zh) * | 2018-05-10 | 2019-04-30 | 南昌航空大学 | 一种基于改进时间反转算法的玻璃纤维复合材料声发射源定位成像方法 |
WO2019113994A1 (zh) * | 2017-12-14 | 2019-06-20 | 南京航空航天大学 | 基于激励波束成型和加权图像融合的music腐蚀监测方法 |
CN110702785A (zh) * | 2019-09-24 | 2020-01-17 | 清华大学 | 频散Lamb波多项式时频域模态分解和缺陷定位方法及装置 |
CN112285651A (zh) * | 2020-10-22 | 2021-01-29 | 温州大学 | 基于特征提取矩阵相似度识别复合材料冲击源位置的方法及系统 |
CN112924542A (zh) * | 2021-01-19 | 2021-06-08 | 中南大学 | 一种岩石、混凝土及充填体的强度质量测量方法、装置、服务器及可读存储介质 |
CN112946070A (zh) * | 2021-01-08 | 2021-06-11 | 温州大学 | 基于传递路径误差补偿来识别冲击源位置的方法及系统 |
CN113687143A (zh) * | 2021-08-12 | 2021-11-23 | 国网上海市电力公司 | 一种电磁波信号幅值衰减与传播距离关系曲线的拟合方法 |
CN114136556A (zh) * | 2021-11-30 | 2022-03-04 | 北京卫星环境工程研究所 | 一种基于波速修正的航天器复材结构泄漏定位方法 |
CN114184895A (zh) * | 2021-12-02 | 2022-03-15 | 广西大学 | 不受检测距离影响的多芯线缺陷定位与分类系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010212818A (ja) * | 2009-03-08 | 2010-09-24 | Univ Of Tokyo | 複数のマイクロフォンにより受信された多チャンネル信号の処理方法 |
CN101865758A (zh) * | 2010-06-12 | 2010-10-20 | 南京航空航天大学 | 基于多重信号分类算法的冲击载荷定位方法 |
-
2013
- 2013-12-16 CN CN201310685400.6A patent/CN103698748B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010212818A (ja) * | 2009-03-08 | 2010-09-24 | Univ Of Tokyo | 複数のマイクロフォンにより受信された多チャンネル信号の処理方法 |
CN101865758A (zh) * | 2010-06-12 | 2010-10-20 | 南京航空航天大学 | 基于多重信号分类算法的冲击载荷定位方法 |
Non-Patent Citations (1)
Title |
---|
钟永腾 等: "基于MUSIC算法的复合材料近场冲击源定位", 《复合材料学报》, vol. 30, no. 2, 30 April 2013 (2013-04-30), pages 189 - 193 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016037400A1 (zh) * | 2014-09-09 | 2016-03-17 | 南京航空航天大学 | 基于能量加权因子的复合材料结构冲击区域定位方法 |
US10697937B2 (en) | 2017-12-14 | 2020-06-30 | Nanjing University Of Aeronautics And Astronautics | MUSIC corrosion monitoring method via excitation beam forming and weighted image fusing |
WO2019113994A1 (zh) * | 2017-12-14 | 2019-06-20 | 南京航空航天大学 | 基于激励波束成型和加权图像融合的music腐蚀监测方法 |
CN109696480A (zh) * | 2018-05-10 | 2019-04-30 | 南昌航空大学 | 一种基于改进时间反转算法的玻璃纤维复合材料声发射源定位成像方法 |
CN109696480B (zh) * | 2018-05-10 | 2021-06-29 | 南昌航空大学 | 一种基于改进时间反转算法的玻璃纤维复合材料声发射源定位成像方法 |
CN109115894B (zh) * | 2018-09-04 | 2020-12-29 | 温州大学激光与光电智能制造研究院 | 一种基于Toeplitz近似法的复合材料板损伤检测方法 |
CN109115894A (zh) * | 2018-09-04 | 2019-01-01 | 温州大学激光与光电智能制造研究院 | 一种基于Toeplitz近似法的复合材料板损伤检测方法 |
CN109406627B (zh) * | 2018-09-11 | 2021-07-20 | 温州大学苍南研究院 | 一种基于旋转不变子空间的复合材料板冲击定位方法 |
CN109406627A (zh) * | 2018-09-11 | 2019-03-01 | 温州大学苍南研究院 | 一种基于旋转不变子空间的复合材料板冲击定位方法 |
CN109406084A (zh) * | 2018-09-13 | 2019-03-01 | 温州大学 | 一种对复合材料产品实现全频域快速冲击监测的方法 |
CN109374740A (zh) * | 2018-09-21 | 2019-02-22 | 南京航空航天大学 | 基于阵列误差校正的合成孔径music损伤定位方法 |
CN110702785B (zh) * | 2019-09-24 | 2020-10-16 | 清华大学 | 频散Lamb波多项式时频域模态分解和缺陷定位方法及装置 |
CN110702785A (zh) * | 2019-09-24 | 2020-01-17 | 清华大学 | 频散Lamb波多项式时频域模态分解和缺陷定位方法及装置 |
CN112285651A (zh) * | 2020-10-22 | 2021-01-29 | 温州大学 | 基于特征提取矩阵相似度识别复合材料冲击源位置的方法及系统 |
CN112285651B (zh) * | 2020-10-22 | 2023-10-20 | 温州大学 | 基于特征提取矩阵相似度识别复合材料冲击源位置的方法及系统 |
CN112946070A (zh) * | 2021-01-08 | 2021-06-11 | 温州大学 | 基于传递路径误差补偿来识别冲击源位置的方法及系统 |
CN112946070B (zh) * | 2021-01-08 | 2023-01-03 | 温州大学 | 基于传递路径误差补偿来识别冲击源位置的方法及系统 |
CN112924542A (zh) * | 2021-01-19 | 2021-06-08 | 中南大学 | 一种岩石、混凝土及充填体的强度质量测量方法、装置、服务器及可读存储介质 |
CN113687143A (zh) * | 2021-08-12 | 2021-11-23 | 国网上海市电力公司 | 一种电磁波信号幅值衰减与传播距离关系曲线的拟合方法 |
CN113687143B (zh) * | 2021-08-12 | 2024-04-05 | 国网上海市电力公司 | 一种电磁波信号幅值衰减与传播距离关系曲线的拟合方法 |
CN114136556A (zh) * | 2021-11-30 | 2022-03-04 | 北京卫星环境工程研究所 | 一种基于波速修正的航天器复材结构泄漏定位方法 |
CN114136556B (zh) * | 2021-11-30 | 2024-03-08 | 北京卫星环境工程研究所 | 一种基于波速修正的航天器复材结构泄漏定位方法 |
CN114184895A (zh) * | 2021-12-02 | 2022-03-15 | 广西大学 | 不受检测距离影响的多芯线缺陷定位与分类系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103698748B (zh) | 2017-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103698748A (zh) | 复合材料中基于信号波速与衰减补偿的2d-music冲击定位方法 | |
CN106960068B (zh) | 一种基于脉冲激励响应频谱的模态阻尼比快速计算方法 | |
CN101865758B (zh) | 基于多重信号分类算法的冲击载荷定位方法 | |
Qiu et al. | A time reversal focusing based impact imaging method and its evaluation on complex composite structures | |
CN105510444B (zh) | 基于超声导波的双向时间反演损伤成像方法 | |
CN105842345B (zh) | 一种光纤布拉格光栅传感网络结构多裂纹定位方法 | |
CN101586997A (zh) | 一种拉索振动基频的计算方法 | |
CN104776966A (zh) | 基于分形理论的板结构冲击监测方法 | |
CN103344709A (zh) | 基于线性映射频散补偿的压电超声相控阵多损伤成像方法 | |
CN113933392B (zh) | 一种基于超声导波的特征融合概率重建损伤定位成像方法 | |
CN106526678A (zh) | 一种反射声波测井的波场分离方法及装置 | |
CN113359183B (zh) | 一种针对极地冰层的震源定位方法 | |
CN104656123A (zh) | 一种测定区域岩体等效波速的方法 | |
Yuan et al. | A single frequency component-based re-estimated MUSIC algorithm for impact localization on complex composite structures | |
CN114235962B (zh) | 一种面向各向异性结构的超声导波成像方法及系统 | |
Zhu et al. | A two-step impact localization method for composite structures with a parameterized laminate model | |
CN104792865A (zh) | 基于分数维的管道超声导波小缺陷识别与定位方法 | |
CN103913722B (zh) | 复合材料格栅结构低速冲击定位方法 | |
CN104374532A (zh) | 航天器在轨泄漏定向方法 | |
CN105259571A (zh) | 一种地层倾角检测方法 | |
CN105717198A (zh) | 一种面向结构冲击定位的单频-重估计music方法 | |
Kundu et al. | A generic framework for application of machine learning in acoustic emission-based damage identification | |
Gorgin et al. | Probability based impact localization in plate structures using an error index | |
Mu et al. | Time reversal damage localization method of ocean platform based on particle swarm optimization algorithm | |
Siano et al. | Acoustic optimization of a high-speed train composite sandwich panel based on analytical and experimental Transmission Loss evaluation integrated by FE/Test correlation analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |