CN110702786B - 基于多尺度奇异吸引子预测误差的梁结构损伤识别方法 - Google Patents
基于多尺度奇异吸引子预测误差的梁结构损伤识别方法 Download PDFInfo
- Publication number
- CN110702786B CN110702786B CN201910938918.3A CN201910938918A CN110702786B CN 110702786 B CN110702786 B CN 110702786B CN 201910938918 A CN201910938918 A CN 201910938918A CN 110702786 B CN110702786 B CN 110702786B
- Authority
- CN
- China
- Prior art keywords
- scale
- beam structure
- attractor
- prediction error
- mran
- 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/045—Analysing solids by imparting shocks to the workpiece and detecting the vibrations or the acoustic waves caused by the shocks
-
- 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/4472—Mathematical theories or simulation
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Algebra (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Acoustics & Sound (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
本发明公开了一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法,该方法具体包括以下步骤:通过传感器系统采集梁结构的加速度响应数据;采用软阈值法对加速度数据进行去噪;对去噪信号实施平稳离散小波分解,得到无降采样的多尺度子信号;应用相空间重构理论对多尺度子信号进行相空间重构,得到多尺度重构吸引子MRA,并对其进行归一化处理,得到MRAN;计算多尺度奇异吸引子预测误差PE,据其判断梁结构的健康状态,并识别其损伤程度。本发明建立的多尺度奇异吸引子预测误差法,突破了激励形式的限制,同时具有损伤敏感性强、识别准确度高的特点,仅采用单点加速度数据即可准确识别梁结构的损伤程度。
Description
技术领域
本发明涉及一种梁结构损伤识别方法,具体涉及一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法,属于结构损伤诊断领域。
背景技术
梁结构在长期服役中,受自身材料退化和复杂环境影响,不可避免发生损伤,早期局部损伤的演化累积会威胁整体结构的可靠运行和安全,甚至导致整体破坏。基于此,近二十年来结构无损检测技术取得迅速进步,在及时发现结构早期损伤方面发挥了不可替代的作用。奇异吸引子预测误差法是随着混沌理论的快速发展而新兴的一种结构无损检测方法,目前正逐步被应用于机械、土木领域的结构损伤诊断。然而,利用小波变换改进传统的奇异吸引子,得到新的多尺度奇异吸引子预测误差,并据其进行梁结构损伤识别的研究还未见报道。
传统的基于奇异吸引子预测误差的梁结构损伤检测方法对激励形式有严格要求,一般需根据结构特征选择具有特定李雅普诺夫指数的混沌激励;而常规激励,如锤击激励、简谐激励等,无法满足以上激励选择的要求,会降低奇异吸引子预测误差法在损伤识别中的准确性和稳定性。
发明内容
本发明所要解决的计算问题是:针对现有技术不足,改进了常规的基于奇异吸引子预测误差的结构损伤检测方法,对振动信号进行平稳小波变换和相空间重构处理,得到新的多尺度奇异吸引子,进而通过多尺度奇异吸引子预测误差判断梁结构健康状态并识别损伤程度。相比于传统方法,本发明建立的多尺度奇异吸引子预测误差法,突破了激励形式的限制,同时具有损伤敏感性强、识别准确度高的特点,仅采用单点加速度数据即可准确识别梁结构的损伤程度,提供了一种简单、快速、准确的梁结构损伤识别方法。
本发明所采用的技术方案为:
一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法,具体步骤包括:
步骤一、布置安装传感器量测系统,采集梁结构的加速度响应数据,记为x={x(tn)},n=1,2,…,N,其中,tn为第n个采样点所在的时刻,N为数据长度;
步骤二、采用小波软阈值法对加速度响应数据进行去噪处理,得到去噪后的加速度响应数据,记为y={y(tn)};
步骤三、对y实施平稳离散小波分解,得到无降采样的多尺度子信号,分别为近似系数aj={aj(tn)}和细节系数dj={dj(tn)},其中j是小波分解尺度;
步骤四、应用相空间重构理论对多尺度子信号进行相空间重构,得到多尺度重构奇异吸引子MRA,并对其进行归一化处理,得到MRAN;
步骤五、计算多尺度奇异吸引子MRAN预测误差,记为PE,据其判断梁结构的健康状态,并识别梁结构的损伤程度。
作为本发明的进一步优选方案,步骤二小波软阈值法具体计算方式如下:
(1)对x进行离散小波变换。
(2)计算信号噪声的方差σ:
其中,d1是第一层离散小波细节系数。
(3)计算通用软阈值h:
作为本发明的进一步优选方案,所述离散小波变换采用小波基函数为rbio4.4,小波分解层数为3。
作为本发明的进一步优选方案,所述步骤三对y实施平稳离散小波分解,选取小波基函数为db4,小波分解层数为9。
作为本发明的进一步优选方案,所述步骤四具体计算方式如下:
(1)选取平稳小波变换的细节系数d9作为重构信号;
(2)通过互信息量法和Cao法分别确定最佳延迟时间τ和最佳嵌入维数m,重构得到奇异吸引子MRA,吸引子相点坐标表示为:
其中,ti为第i个采样点所在的时刻,M=N-(m-1)τ+1,表示重构相空间中的相点数量;
(3)采用线性归一化方法,将MRA的所有维度都规范到区间[-1,1],得到MRAN。
作为本发明的进一步优选方案,所述最佳延迟时间τ和最佳嵌入维数m分别取τ=8,m=4。
作为本发明的进一步优选方案,所述步骤五中多尺度奇异吸引子MRAN预测误差PE的具体计算方式如下:
(7)计算MRANB基于MRANA的2-范数预测误差,记为PEBA,F:
(8)多次重复过程(2)~(7),取其平均值作为多尺度奇异吸引子预测误差PE。
作为本发明的进一步优选方案,所述邻近点数量p=3,时间演化步长Δt=2。
本发明具有以下有益效果:
本发明提供一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法,该方法改进了常规的基于奇异吸引子预测误差的结构损伤检测方法,对振动信号进行平稳小波变换和相空间重构处理,得到新的多尺度奇异吸引子,进而通过多尺度奇异吸引子预测误差判断梁结构健康状态并识别损伤程度。相比于传统方法,本发明建立的多尺度奇异吸引子预测误差法,突破了激励形式的限制,同时具有损伤敏感性强、识别准确度高的特点,仅采用单点加速度数据即可准确识别梁结构的损伤程度,提供了一种简单、快速、准确的梁结构损伤识别方法。
附图说明
图1是本发明方法流程示意图。
图2是原始量测梁结构加速度响应信号x。
图3是软阈值去噪后梁结构加速度响应信号y。
图4是无降采样的多尺度子信号d9。
图5是现有技术中的相空间重构奇异吸引子RAN。
图6是本发明所公开的多尺度奇异吸引子MRAN。
图7基于多尺度奇异吸引子预测误差的梁结构损伤识别结果。
具体实施方式
下面结合附图对本发明作进一步说明:
本发明在奇异吸引子预测误差法的基础上,发展了一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法。相比于常规奇异吸引子预测误差法,本发明方法不仅突破了激励形式的限制,还能加强损伤信息的特征,简单、快速、准确地识别梁结构损伤。请参照图1,本发明提出一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法,具体步骤包括:
1.布置安装传感器量测系统,采集梁结构的加速度响应数据,记为x={x(tn)},n=1,2,…,N,N为数据长度;
2.对进行软阈值去噪处理,具体为:
(1)选取小波基函数为rbio4.4,小波分解层数为3,对x进行离散小波变换;
(2)计算信号噪声的方差σ:
其中,d1是第一层离散小波系数;
(3)计算通用软阈值h:
3.选取小波基函数为db4,小波分解层数为9,对y实施平稳离散小波分解,得到无降采样的多尺度子信号,分别为近似系数aj={aj(tn)}和细节系数dj={dj(tn)},其中j=1,2,…,9是小波分解尺度;
4.应用相空间重构理论对多尺度子信号进行相空间重构,具体为:
(1)选取平稳小波变换的细节系数d9作为重构信号;
(2)通过互信息量法和Cao法分别确定最佳延迟时间τ=8和最佳嵌入维数m=4,重构得到奇异吸引子MRA,吸引子相点坐标表示为:
其中,ti为第i个采样点所在的时刻,M=N-(m-1)τ+1,表示重构相空间中的相点数量;
(3)采用线性归一化方法,将MRA的所有维度都规范到区间[-1,1],得到MRAN。
5.计算多尺度奇异吸引子预测误差,据其进行梁结构的损伤识别,具体为:
(7)计算MRANB基于MRANA的2-范数预测误差,记为PEBA,F:
(8)多次重复过程(2)~(7),取其平均值作为多尺度奇异吸引子预测误差PE;
(9)根据PE的数值大小判断梁结构的健康状态,并识别其损伤程度。
实施例:
为验证本发明一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法的有效性,开展物理模型实验,提取梁结构的加速度响应进行分析。
参照图1,在本实施例中,梁结构的几何尺寸为:长度为L=0.5m,横截面尺寸H×B=0.019m×0.012m;物理参数为:弹性模量E=1.84×1011Pa,泊松比ν=0.3,材料密度ρ=7750kg/m3;边界条件为一端固支,一端自由;荷载采用锤击激励,激励位置距自由端45mm;损伤位置距固定端133.3mm;加速度响应通过Utel数据采集系统收集得到,传感器距自由端90mm,采样频率为1.28kHz,分析频率5000Hz,采样时间5.12s。采集到的加速度响应数据记为x,如图2所示。
对x进行小波软阈值去噪处理,应用公式(1)~(3),得到去噪后的加速度响应数据y,如图3所示。选取小波基函数为“db4”,小波分解层数为9,对y实施平稳离散小波分解,得到无降采样的多尺度子信号,分别为近似系数a9和细节系数d9,d9如图4所示。应用相空间重构理论对多尺度子信号d9进行相空间重构,通过互信息量法和Cao法获得最佳延迟时间τ=8和最佳嵌入维数m=4,应用公式(4),得到多尺度奇异吸引子MRA;采用线性归一化方法,将MRA的所有维度都规范到区间[-1,1],得到MRAN。常规重构奇异吸引子RAN如图5所示,本发明公开的多尺度奇异吸引子MARN如图6所示,其中下标数字表示坐标维度。应用公式(5)~(7),计算多尺度奇异吸引子MRAN预测误差PE,据其判断梁结构的健康状态,并识别梁结构的损伤程度,识别结果如图7所示。从图7中可以看出,随着梁结构损伤程度的增加,PERAN没有显著的变化规律,无法据其判断结构损伤状态;相反地,PEMRAN随着损伤程度的增加呈显著的单调递增规律,因此PEMRAN的数值大小反映了梁结构当前的健康情况,仅采用单点加速度数据即可准确识别梁结构的损伤程度,提供了一种简单、快速、准确的梁结构损伤识别方法。
以上所述仅是本发明专利的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明专利原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明专利的保护范围。
Claims (7)
1.一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法,其特征在于,具体步骤包括:
步骤一、布置安装传感器量测系统,采集梁结构的加速度响应数据,记为x={x(tn)},其中,tn为第n个采样点所在的时刻,n=1,2,…,N,N为数据长度;
步骤二、采用小波软阈值法对加速度响应数据进行去噪处理,得到去噪后的加速度响应数据,记为y={y(tn)};
步骤三、对y实施平稳离散小波分解,得到无降采样的多尺度子信号,分别为近似系数aj={aj(tn)}和细节系数dj={dj(tn)},其中j是小波分解尺度;
步骤四、应用相空间重构理论对多尺度子信号进行相空间重构,得到多尺度重构奇异吸引子MRA,并对其进行归一化处理,得到MRAN;
步骤五、计算多尺度奇异吸引子MRAN预测误差,记为PE,据其判断梁结构的健康状态,并识别梁结构的损伤程度;其中,多尺度奇异吸引子MRAN预测误差PE的具体计算方式如下:
(7)计算MRANB基于MRANA的2-范数预测误差,记为PEBA,F:
(8)多次重复过程(2)~(7),取其平均值作为多尺度奇异吸引子预测误差PE。
3.根据权利要求2所述的一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法,其特征在于,所述离散小波变换采用小波基函数为rbio4.4,小波分解层数为3。
4.根据权利要求1所述的一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法,其特征在于,所述步骤三对y实施平稳离散小波分解,选取小波基函数为db4,小波分解层数为9。
6.根据权利要求5所述的一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法,其特征在于,所述最佳延迟时间τ和最佳嵌入维数m分别取τ=8,m=4。
7.根据权利要求1所述的一种基于多尺度奇异吸引子预测误差的梁结构损伤识别方法,其特征在于,所述步骤五中邻近点数量p=3,时间演化步长Δt=2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910938918.3A CN110702786B (zh) | 2019-09-30 | 2019-09-30 | 基于多尺度奇异吸引子预测误差的梁结构损伤识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910938918.3A CN110702786B (zh) | 2019-09-30 | 2019-09-30 | 基于多尺度奇异吸引子预测误差的梁结构损伤识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110702786A CN110702786A (zh) | 2020-01-17 |
CN110702786B true CN110702786B (zh) | 2022-06-17 |
Family
ID=69197929
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910938918.3A Active CN110702786B (zh) | 2019-09-30 | 2019-09-30 | 基于多尺度奇异吸引子预测误差的梁结构损伤识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110702786B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111239549A (zh) * | 2020-02-18 | 2020-06-05 | 国网信通亿力科技有限责任公司 | 一种基于离散小波变换的配电故障快速定位方法 |
CN112131718B (zh) * | 2020-09-07 | 2022-06-10 | 中国科学院力学研究所 | 一种基于结构时域响应的点阵夹层板损伤识别方法 |
CN118277719B (zh) * | 2024-06-04 | 2024-08-06 | 山东神力索具有限公司 | 一种用于计数型矿用索具计算的数据智能处理方法 |
CN118603617A (zh) * | 2024-08-05 | 2024-09-06 | 北京航空航天大学 | 一种基于吸引子的定常激励下机电设备的状态检测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000007327A1 (en) * | 1998-07-29 | 2000-02-10 | Jong Uk Choi | Method of evaluating encryption algorithms using chaos analysis |
CN103940905A (zh) * | 2014-04-16 | 2014-07-23 | 山东农业大学 | 基于平稳小波变换和分形分析的梁结构损伤检测方法 |
CN104849348A (zh) * | 2015-05-13 | 2015-08-19 | 河海大学 | 一种基于奇异值分解-小波变换的梁结构损伤检测方法 |
CN104881567A (zh) * | 2015-04-24 | 2015-09-02 | 浙江工业大学 | 一种基于统计模型的桥梁健康监测数据小波降噪方法 |
WO2017144757A1 (es) * | 2016-02-25 | 2017-08-31 | Universitat Politècnica De València | Método de ensayo del estado de materiales |
JP2019039787A (ja) * | 2017-08-24 | 2019-03-14 | 国立大学法人福井大学 | 構造物の打音検査装置及び方法 |
CN110231403A (zh) * | 2019-06-10 | 2019-09-13 | 中设设计集团股份有限公司 | 移动荷载响应小波包分析支持的长跨桥梁在线实时损伤识别方法 |
-
2019
- 2019-09-30 CN CN201910938918.3A patent/CN110702786B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000007327A1 (en) * | 1998-07-29 | 2000-02-10 | Jong Uk Choi | Method of evaluating encryption algorithms using chaos analysis |
CN103940905A (zh) * | 2014-04-16 | 2014-07-23 | 山东农业大学 | 基于平稳小波变换和分形分析的梁结构损伤检测方法 |
CN104881567A (zh) * | 2015-04-24 | 2015-09-02 | 浙江工业大学 | 一种基于统计模型的桥梁健康监测数据小波降噪方法 |
CN104849348A (zh) * | 2015-05-13 | 2015-08-19 | 河海大学 | 一种基于奇异值分解-小波变换的梁结构损伤检测方法 |
WO2017144757A1 (es) * | 2016-02-25 | 2017-08-31 | Universitat Politècnica De València | Método de ensayo del estado de materiales |
JP2019039787A (ja) * | 2017-08-24 | 2019-03-14 | 国立大学法人福井大学 | 構造物の打音検査装置及び方法 |
CN110231403A (zh) * | 2019-06-10 | 2019-09-13 | 中设设计集团股份有限公司 | 移动荷载响应小波包分析支持的长跨桥梁在线实时损伤识别方法 |
Non-Patent Citations (3)
Title |
---|
Damage detection using modal frequency curve and squared residual wavelet coefficients-based damage indicator;R. Sampaio et al;《Journal of Sound and Vibration》;19990419;第226卷(第5期);第1029-1042页 * |
ONTHEUSEOFATTRACTOR DIMENSIONASA FEATURE IN STRUCTURALHEALTHMONITORING;J. M. Nichols et al;《Mechanical Systems and Signal Processing》;20020805;第17卷(第6期);第1305-1320页 * |
基于结构响应的损伤诊断方法及其应用;刘毅 等;《东南大学学报(自然科学版)》;20100730;第40卷(第4期);第810-815页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110702786A (zh) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110702786B (zh) | 基于多尺度奇异吸引子预测误差的梁结构损伤识别方法 | |
CN108959778B (zh) | 一种基于退化模式一致性的航空发动机剩余寿命预测方法 | |
CN109932179B (zh) | 一种基于ds自适应谱重构的滚动轴承故障检测方法 | |
CN116108008A (zh) | 一种装饰材料甲醛检测数据处理方法 | |
Ince et al. | A machine learning approach for locating acoustic emission | |
CN102805620A (zh) | 一种房颤检测方法和装置 | |
CN107748209B (zh) | 一种结构损伤检测方法 | |
CN108645620A (zh) | 一种基于信息熵和多尺度形态学的滚动轴承早期故障诊断方法 | |
CN117454096A (zh) | 一种电机生产质量检测方法及系统 | |
CN111964909A (zh) | 滚动轴承运行状态检测方法、故障诊断方法及系统 | |
CN116718218B (zh) | 一种光电集成式传感芯片测试数据处理方法 | |
CN117665935B (zh) | 一种破碎岩体支护施工过程的监测数据处理方法 | |
CN109117450B (zh) | 振动测试数据最佳分析长度的确定方法 | |
CN116977708B (zh) | 一种基于自适应聚合可视图的轴承智能诊断方法及系统 | |
CN108629770B (zh) | 基于支持向量机的超声图像分割方法 | |
CN116304549A (zh) | 一种隧道健康监测数据小波阈值去噪方法 | |
CN114755010A (zh) | 一种旋转机械振动故障诊断方法及其系统 | |
CN115169401A (zh) | 一种基于多尺度DenseNet-ResNet-GRU模型的刀具磨损量预测方法 | |
CN116625678B (zh) | 基于平均峭度反卷积网络的故障诊断方法及系统 | |
CN112880726A (zh) | 基于变分模态分解样本熵的传感器故障诊断方法和装置 | |
CN116881781A (zh) | 一种运行模态阻尼识别方法、损伤检测方法、系统及设备 | |
CN104278990B (zh) | 测井数据质量恢复方法及装置 | |
CN114021275A (zh) | 一种基于深度卷积模糊系统的滚动轴承故障诊断方法 | |
JP7472012B2 (ja) | 異常信号抽出装置 | |
Deuschle et al. | Robust sensor spike detection method based on dynamic time warping |
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 |