CN108195326B - 一种基于经验模态分解的撞击定位方法 - Google Patents
一种基于经验模态分解的撞击定位方法 Download PDFInfo
- Publication number
- CN108195326B CN108195326B CN201711340879.4A CN201711340879A CN108195326B CN 108195326 B CN108195326 B CN 108195326B CN 201711340879 A CN201711340879 A CN 201711340879A CN 108195326 B CN108195326 B CN 108195326B
- Authority
- CN
- China
- Prior art keywords
- signal
- shock
- mode decomposition
- localization method
- empirical mode
- 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
- 230000035939 shock Effects 0.000 title claims abstract description 36
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000004807 localization Effects 0.000 title claims abstract description 18
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- 239000011324 bead Substances 0.000 claims abstract description 6
- 238000012544 monitoring process Methods 0.000 claims description 19
- 230000006698 induction Effects 0.000 claims description 6
- 238000009792 diffusion process Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 description 9
- 239000002131 composite material Substances 0.000 description 4
- 239000003365 glass fiber Substances 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000008188 pellet Substances 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/0052—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes measuring forces due to impact
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于经验模态分解的撞击定位方法,通过布设在待测结构上设置的传感器阵列以及小球的自由落体运动，对应力波接收端采集到的撞击信号进行经验模态分解，得到有限多个本征模函数，再结合能量重心定位方法，确定撞击源的位置坐标。本发明公开的撞击定位方法是通过经验模态分解法与能量重心定位法的结合来确定撞击位置即非损伤定位。本发明公开的撞击定位方法与传统的傅里叶分析和小波分析相比精度更高。
Description
技术领域
本发明涉及一种基于经验模态分解的撞击定位方法，尤其涉及一种经验模态分解和能量重心定位联合作用下的撞击定位方法，属于被动式结构健康监测领域。
背景技术
撞击定位是结构健康监测技术中至关重要的一步。不管是各向同性还是各向异性结构，都不可避免会受到外来物体的撞击，严重时可能导致结构内部出现损伤，最终致使结构失效。遭到撞击之后，结构内部会产生应力波，应力波会以撞击源为起点沿结构向外扩散。如果在结构表面不同的位置粘贴一些传感器阵列，撞击产生的应力波便会被记录在各个传感器中，只要对传感器接收到的信号进行适当分析和处理，就能成功实现对撞击源的定位。
能够实现精确定位的前提是提取到真实可靠的特征信号。在以前的研究中，国内外都将重点放在了定位的方法上。比较典型的基于波达时间的定位方法，往往由于各向异性结构的频散问题难以提取到准确的波达时间；基于神经网络的定位方法需要大量的样本训练，费时费力。对于撞击信号的处理一般多用傅里叶变换和小波变换进行滤波以及提取信号特征。然而，这两种方法都需要利用特定的函数形式(傅里叶的三角函数，小波变换的小波基函数)对数据进行分解，且傅里叶变换不具备时域特征，小波变换容易产生虚假频率。因此目前需要一种不需要利用特定的函数形式对数据进行分解的具有自适应能力的先进的时频域信号处理方法以获得有效的撞击信号。
发明内容
根据现有问题，本发明公开了一种基于经验模态分解的撞击定位方法：该方法采用经验模态分解和能量重心法对撞击源进行定位，具体包括以下步骤：
S1：在待测结构上建立直角坐标系，设置由四个压电传感器组成的撞击监测区域，且区域中心位置坐标为(0,0)；
S2:将上述传感器分别记为是s1，s2，s3，s4，并将上述监测区域由内到外划分为49个小区域(0-48)；
S3：利用小球固定高度的自由落体运动使待测结构中产生应力波，应力波在监测区域中的某一点产生并迅速向外扩散传到所述压电传感器s1，s2，s3，s4，多通道信息采集系统将记录所述压电传感器接收到的撞击信号；
S4：按照S3的方式完成剩余监测区域的信号采集。每次撞击点均为S2所述小区域的中心位置；
S5：将采集到信号进行有效信号的提取，然后将有效信号分为四组、分别进行不同的分析：第一组信号不再作处理，第二组信号作经验模态分解，第三组作傅里叶分析，第四组作小波分析；
S6：通过S5中所述的处理后的信号，分别计算对应信号的能量p1，p2，p3，p4；
S7:利用所述四个压电传感器s1，s2，s3，s4的位置坐标以及S6所求的对应信号的能量p1，p2，p3，p4求出其能量重心位置坐标，即撞击源位置坐标。
进一步的，步骤S1中所述的待测结构为复合材料薄板结构。
进一步的，所述步骤S3中的小球自由落体运动由固定长度的塑料管完成且撞击事件不会致使结构出现损伤。
进一步的，步骤S5中所述的对信号进行经验模态分解得到有限多个本征模函数，取前五阶本征模函数进行定位且将误差最小的作为撞击源位置坐标。
本发明公开的一种基于经验模态分解的撞击定位方法，通过撞击定位方法能有效实现撞击源位置的识别，通过经验模态分解提高了传统信号处理后的定位精度，并且基于经验模态分解的撞击定位方法有效遏制了由能量重心定位法带来的额外误差。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明中布置玻璃纤维复合材料板及撞击监测区域的平面图；
图2为本发明中撞击监测区域内部划分的平面图；
图3为本发明实施例中的一条撞击信号波形图；
图4为本发明实施例中的一条撞击信号进行经验模态分解后的波形图。
具体实施方式
为使本发明的技术方案和优点更加清楚，下面结合本发明实施例中的附图，
对本发明实施例中的技术方案进行清楚完整的描述：
本发明公开了一种基于经验模态分解的撞击定位方法,采用经验模态分解和能量重心法相结合对撞击源进行定位，具体包括以下步骤：
S1：在待测结构上建立直角坐标系，设置由四个压电传感器组成的撞击监测区域，且区域中心位置坐标为(0,0)；
S2:将上述传感器分别记为是s1，s2，s3，s4，并将上述监测区域由内到外划分为49个小区域(0-48)；
S3：利用小球固定高度的自由落体运动使待测结构中产生应力波，应力波在监测区域中的某一点产生并迅速向外扩散传到所述压电传感器s1，s2，s3，s4，多通道信息采集系统将记录所述压电传感器接收到的撞击信号；
S4：按照S3的方式完成剩余监测区域的信号采集。每次撞击点均为S2所述小区域的中心位置；
S5：将采集到信号进行有效信号的提取，然后将有效信号分为四组、分别进行不同的分析：第一组信号不再作处理，第二组信号作经验模态分解，第三组作傅里叶分析，第四组作小波分析；
S6：得到S5中所述处理过的信号后，分别计算对应信号的能量p1，p2，p3，p4；
S7:利用所述四个压电传感器s1，s2，s3，s4的位置坐标以及S6所求的对应信号的能量p1，p2，p3，p4求出其能量重心位置坐标，即撞击源位置坐标。
进一步的，步骤S1中所述的待测结构为复合材料薄板结构。
进一步的，所述步骤S3中的小球自由落体运动由固定长度的塑料管完成且撞击事件不会致使结构出现损伤。
进一步的，步骤S5中所述的对信号进行经验模态分解得到有限多个本征模函数，取前五阶本征模函数进行定位且将误差最小的作为撞击源位置坐标。
实施例：
步骤1：取一块如图1所示复合材料薄板，本实例中选取厚度为0.3cm的玻璃纤维复合材料板，板1尺寸为123cm×103cm。板1上粘接有四个规格相同的压电传感器s1，s2，s3，s4，直径6.5mm，厚度0.25mm。四个传感器围成20cm×20cm的撞击监测区域2。所述撞击监测区域的中心与板1中心重合，为坐标原点，如图1所示。
步骤2：为方便记录数据，将所述撞击监测区域2由内到外有序划分为49个小的撞击区域，如图2所示；
步骤3：为了保证结构在试验中的稳定性，将步骤1中所述复合材料板1的四条边固定。利用小球的自由落体运动对板1中的监测区域2进行撞击。假定撞击源为点源，即每一次小球的落点均为各个小区域的中心；
步骤4：选取监测区域2中的某一撞击区域，利用多通道信息采集系统对步骤3中触发的撞击事件进行信号采集。一般地，信号超过预设阈值即定义为发生撞击事件，系统会自动采集信号。本实例中阈值设为0.25V，采样频率为40.54kHZ；
步骤5：按照步骤4的方式完成剩余小监测区域测的信号采集，典型的撞击信号如图3所示；
步骤6：选取某一个小监测区域采集到的撞击信号，将其作基本处理后分为四组，分别进行不同的分析：第一组信号不再作处理，第二组信号作经验模态分解，第三组作傅里叶分析，第四组作小波分析；
其中，信号作经验模态分解后可以用下式表示：
其中，s(t)为作基本处理后的撞击信号，c(t)为经验模态分解后的本征模函数，r(t)为分解后剩余的余项，n为筛选终止后获得的本征模函数的个数。典型的经验模态分解后的信号如图4所示(IMFi代表第i阶本征模函数，residue代表余项)；
本实例中的傅里叶变换采用快速傅里叶变换，小波变换采用db5小波基进行离散小波变换；
步骤7：按照步骤5和6，将剩余其他撞击区域的信号都分为四组进行处理；
步骤8：根据下式计算各个信号的能量：
其中，s(t)为处理后的撞击信号，t0，tn分别代表撞击信号的起始时间和终止时间。起始时间指信号值第一次超过预设的信号阈值。求出的信号能量分别为p1，p2，p3，p4；
步骤9：能量重心定位法确定撞击源位置坐标。利用四个传感器位置坐标和步骤8中所求信号能量p1，p2，p3，p4，按照下式估算出撞击源位置坐标：
其中，xi和yi分别代表四个传感器的横纵坐标，xc和yc代表估算出的当前撞击事件的位置坐标。
步骤10：按照步骤7,8,9将其他所有小撞击区域的位置坐标计算出来，与真实位置坐标进行比较，求出误差。本实例中的误差属于绝对误差，即估计撞击位置与真实撞击位置之间的直线距离。
下面是本实施例得到的定位结果(选取了10个代表性位置)：
表1原始信号定位结果
表2经验模态分解后定位结果
表3快速傅里叶变换后定位结果
表4离散小波变换后定位结果
根据表1和表2，经验模态分解可以应用于撞击定位并具有较好的定位精度，且有效遏制了由能量重心定位带来的额外误差；对比表2和表3，表4可知，信号经过经验模态分解后的定位结果要比传统的傅里叶分析和小波变换精度更高。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
Claims (2)
1.一种基于经验模态分解的撞击定位方法，其特征在于：该方法采用基于经验模态分解的信号处理方法和能量重心定位方法对撞击位置进行定位，具体包括以下步骤：
S1：在待测结构上建立直角坐标系，设置由四个压电传感器组成的撞击监测区域，设置该撞击监测区域的中心位置坐标为(0,0)；
S2:将上述传感器分别记为s1，s2，s3，s4，并将上述监测区域由内到外划分为49个小区域；
S3：利用小球固定高度的自由落体运动使待测结构中产生应力波，应力波在监测区域中的某一点产生并向外扩散传到所述压电传感器s1，s2，s3，s4，采用多通道信息采集系统记录所述压电传感器接收到的撞击信号；
S4：按照S3的方式完成剩余监测区域的信号采集，每次撞击点均为S2所述小区域的中心位置；
S5：将S4中采集到信号进行有效信号的提取，然后将有效信号分为四组、分别进行不同的分析：第一组信号不再作处理，第二组信号作经验模态分解，第三组作傅里叶分析，第四组作小波分析；
S6：对S5中分析后的信号分别计算第一组至第四组对应信号的能量p1，p2，p3，p4；
S7:利用所述四个压电传感器s1，s2，s3，s4的位置坐标以及S6所求的对应信号的能量p1，p2，p3，p4求出其能量重心位置坐标，即撞击源位置坐标。
2.根据权利要求1所述的一种基于经验模态分解的撞击定位方法，其特征还在于：步骤S5中所述信号经过经验模态分解后，统一取前五阶本征模函数进行定位,误差最小的作为撞击源位置坐标。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711340879.4A CN108195326B (zh) | 2017-12-14 | 2017-12-14 | 一种基于经验模态分解的撞击定位方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711340879.4A CN108195326B (zh) | 2017-12-14 | 2017-12-14 | 一种基于经验模态分解的撞击定位方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108195326A CN108195326A (zh) | 2018-06-22 |
| CN108195326B true CN108195326B (zh) | 2019-08-20 |
Family
ID=62574149
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711340879.4A Active CN108195326B (zh) | 2017-12-14 | 2017-12-14 | 一种基于经验模态分解的撞击定位方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108195326B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113616217B (zh) * | 2021-10-12 | 2021-12-21 | 深圳市倍轻松科技股份有限公司 | 基线漂移曲线的生成方法和装置 |
| CN114065487B (zh) * | 2021-11-03 | 2024-10-29 | 大连君晟科技有限责任公司 | 一种基于误差函数的结构冲击定位方法 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5779241A (en) * | 1995-06-02 | 1998-07-14 | D'costa; Joseph F. | Apparatus and method for measurement of position and velocity |
| CN103019446A (zh) * | 2012-12-20 | 2013-04-03 | 杜朝亮 | 基于波传播时间及能量函数的触击位置和能量测量方法 |
| CN103076914A (zh) * | 2012-12-20 | 2013-05-01 | 杜朝亮 | 一种基于能量分布向量比的触击位置和能量测量方法 |
| CN103438794A (zh) * | 2013-08-13 | 2013-12-11 | 南京航空航天大学 | 一种基于压电传感网络的工程结构冲击位置确定方法 |
| CN104931806A (zh) * | 2014-09-13 | 2015-09-23 | 国家电网公司 | 一种基于经验模态理论的扰动源自动化定位方法 |
| CN105067712A (zh) * | 2015-07-23 | 2015-11-18 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 复合材料结构的损伤监测方法、装置和系统 |
| CN106645406A (zh) * | 2016-12-02 | 2017-05-10 | 北京空间飞行器总体设计部 | 一种航天器遭受空间碎片撞击的定位系统及定位方法 |
| CN106643484A (zh) * | 2016-10-17 | 2017-05-10 | 南京航空航天大学 | 基于瞬态分形计算和类二维正态分布拟合的冲击定位方法 |
| CN107025343A (zh) * | 2017-04-01 | 2017-08-08 | 司靓 | 一种面向复合材料结构的冲撞击监控与能评估技术 |
| CN107228797A (zh) * | 2017-06-19 | 2017-10-03 | 南京航空航天大学 | 冲击定位方法及装置 |
-
2017
- 2017-12-14 CN CN201711340879.4A patent/CN108195326B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5779241A (en) * | 1995-06-02 | 1998-07-14 | D'costa; Joseph F. | Apparatus and method for measurement of position and velocity |
| CN103019446A (zh) * | 2012-12-20 | 2013-04-03 | 杜朝亮 | 基于波传播时间及能量函数的触击位置和能量测量方法 |
| CN103076914A (zh) * | 2012-12-20 | 2013-05-01 | 杜朝亮 | 一种基于能量分布向量比的触击位置和能量测量方法 |
| CN103438794A (zh) * | 2013-08-13 | 2013-12-11 | 南京航空航天大学 | 一种基于压电传感网络的工程结构冲击位置确定方法 |
| CN104931806A (zh) * | 2014-09-13 | 2015-09-23 | 国家电网公司 | 一种基于经验模态理论的扰动源自动化定位方法 |
| CN105067712A (zh) * | 2015-07-23 | 2015-11-18 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 复合材料结构的损伤监测方法、装置和系统 |
| CN106643484A (zh) * | 2016-10-17 | 2017-05-10 | 南京航空航天大学 | 基于瞬态分形计算和类二维正态分布拟合的冲击定位方法 |
| CN106645406A (zh) * | 2016-12-02 | 2017-05-10 | 北京空间飞行器总体设计部 | 一种航天器遭受空间碎片撞击的定位系统及定位方法 |
| CN107025343A (zh) * | 2017-04-01 | 2017-08-08 | 司靓 | 一种面向复合材料结构的冲撞击监控与能评估技术 |
| CN107228797A (zh) * | 2017-06-19 | 2017-10-03 | 南京航空航天大学 | 冲击定位方法及装置 |
Non-Patent Citations (7)
| Title |
|---|
| Detecting position dependent tremor with the Empirical mode decomposition;André Lee等;《Journal of Clinical Movement Disorders》;20151231;第1-6页 |
| 一种基于时间反转的冲击定位方法;吴早凤等;《武汉大学学报（工学版）》;20140831;第47卷(第4期);第520-526页 |
| 压电智能结构冲击定位方法的实验研究;周晚林等;《力学与实践》;20031231;第25卷;第37-40页 |
| 基于光纤传感网络的变截面复合材料低速冲击定位;陆观等;《传感技术学报》;20141231;第27卷(第12期);第1632-1636页 |
| 基于冲击能量大小的冲击力位置识别方法;王兴君等;《电子测量技术》;20151130;第38卷(第11期);第10-13页 |
| 基于小波包能量特征向量的光纤布拉格光栅低速冲击定位;郭飞等;《振动与冲击》;20170428;第36卷(第8期);第184-189页 |
| 经验模态分解方法在结构冲击信号分析中的应用;姚熊亮等;《中国舰船研究》;20060831;第1卷(第4期);第11-15页、第47页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108195326A (zh) | 2018-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pomponi et al. | Wavelet based approach to signal activity detection and phase picking: Application to acoustic emission | |
| CN102426195B (zh) | 一种结构损伤监测系统 | |
| CN105488795B (zh) | 一种复合材料损伤识别方法 | |
| CN108195326B (zh) | 一种基于经验模态分解的撞击定位方法 | |
| EP2846175A3 (en) | Seismic survey analysis | |
| CN103076914A (zh) | 一种基于能量分布向量比的触击位置和能量测量方法 | |
| CN105640545A (zh) | 一种胎儿心电信号提取方法及装置 | |
| CN105445697B (zh) | 一种低成本低功耗的声源定向方法 | |
| WO2019023941A1 (zh) | 触控位置的确定方法、电容触控装置以及电容触控终端 | |
| CN105091732A (zh) | 一种变压器绕组变形的检测方法及系统 | |
| CN105403622A (zh) | 一种板材损伤识别定位方法 | |
| CN108345033A (zh) | 一种微地震信号时频域初至检测方法 | |
| CN106646034A (zh) | 一种避雷器结构健康在线监测分析系统 | |
| CN103020269B (zh) | 一种数据验证的方法以及装置 | |
| CN104749249B (zh) | 一种基于超声波技术的检测种子净度的方法 | |
| CN105572232B (zh) | 低频电磁超声导波接收信号的放大方法和装置 | |
| CN104698080A (zh) | 一种利用兰姆波对结构损伤进行状态监测的方法 | |
| CN110057918B (zh) | 强噪声背景下的复合材料损伤定量识别方法及系统 | |
| CN110049423A (zh) | 一种利用广义互相关和能量谱检测麦克风的方法和系统 | |
| RU2011119089A (ru) | Прецизионное измерение формы колебаний | |
| CN103674177B (zh) | 一种信号处理方法和装置 | |
| CN102520075A (zh) | 基于谐波小波包的直升机部件声发射信号传播特性的测试方法 | |
| CN105527180B (zh) | 一种金属材料弯曲疲劳试验方法 | |
| CN104545890B (zh) | 一种基于心电信号的身份识别方法和装置 | |
| CN116881781A (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 |