CN113155973A - 一种基于自适应奇异值分解的梁损伤识别方法 - Google Patents
一种基于自适应奇异值分解的梁损伤识别方法 Download PDFInfo
- Publication number
- CN113155973A CN113155973A CN202110487297.9A CN202110487297A CN113155973A CN 113155973 A CN113155973 A CN 113155973A CN 202110487297 A CN202110487297 A CN 202110487297A CN 113155973 A CN113155973 A CN 113155973A
- Authority
- CN
- China
- Prior art keywords
- damage
- dimension
- singular value
- value decomposition
- matrix
- 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
- 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/12—Analysing solids by measuring frequency or resonance 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
- G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
- G01N29/4418—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with a model, e.g. best-fit, regression analysis
-
- 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/023—Solids
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Acoustics & Sound (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于自适应奇异值分解的梁损伤识别方法,包括:S1:通过激光位移传感器获取损伤梁的振动位移,利用经验模态分解提取损伤梁的模态振型和模态频率特征信号;S2:利用最优粒子逐维算法自动优化奇异值分解矩阵维度参数,然后将梁模态振型信号生成迹矩阵;S3:通过将提取的频率特征信号画在损伤结构的空间全域上,定义适度函数确定峰值点的位置定位损伤的位置;S4:通过基于识别的位置建立损伤结构BSWI模型计算损伤程度不同损伤程度的固有频率,建立固有频率与损伤程度对应关系数据库。该方法具有抗噪声干扰能力强和自适应优化奇异值分解矩阵维度参数的特点,能够对损伤机械结构进行快速精确定位损伤位置和深度。
Description
技术领域
本发明涉及梁损伤识别方法,具体是指一种基于自适应奇异值分解的梁损伤识别方法。
背景技术
梁是一种具有较好的抗弯能力构件,在工程上作为机械结构的受力承载体取得了广泛的应用。机械结构的突发损伤由人为引起或自然灾害导致,难以预料。将基于小波有限元损伤模型求解出的模态振型分解到不同频带中去,实现不同频段信号的分离,但小波分解得到的细节信号对结构损伤灵敏度不高,仍然无法定位损伤;然而,现有的分形维数具有度量信号复杂程度特点,利用适合工程应用的波形容量维数对小波分解的细节信号进行维数估计,根据所得到的波形容量维数的峰值点来定位损伤位置,过程复杂定位精度不高。
发明内容
本发明的目的是为了克服现有技术存在的缺点和不足,而提供一种基于自适应奇异值分解的梁损伤识别方法,该方法具有抗噪声干扰能力强和自适应优化奇异值分解矩阵维度参数的特点,能够对损伤机械结构进行快速精确定位损伤位置和深度。
为实现上述目的,本发明的技术方案是包括以下步骤:
S1:通过激光位移传感器获取损伤梁的振动位移,利用经验模态分解提取损伤梁的模态振型和模态频率特征信号;
S2:利用最优粒子逐维算法自动优化奇异值分解矩阵维度参数,然后将梁模态振型信号生成迹矩阵;
S3:通过将提取的频率特征信号画在损伤结构的空间全域上,定义适度函数确定峰值点的位置定位损伤的位置;
S4:通过基于识别的位置建立损伤结构BSWI模型计算损伤程度不同损伤程度的固有频率,建立固有频率与损伤程度对应关系数据库;
S5:将测量频率作为模糊系统的输入,通过智能搜索损伤程度诊断据库寻找最近匹配的损伤程度,实现梁损伤深度的定量识别。
进一步设置是所述的步骤S1中,通过激光位移传感器获取损伤梁的振动位移,利用经验模态分解提取损伤梁的模态振型和频率特征信号。
进一步设置是所述的步骤S2中,将前几阶振型信号生成迹矩阵,利用最优粒子逐维变异的粒子群算法优化最佳自适应奇异值分解矩阵维度参数是:
a.利用最优粒子逐维变异的粒子群算法优化最佳自适应奇异值分解矩阵维度参数;
b.计算任意一阶模态信息生成列向量X=[x(1)x(2)L x(i)L x(N)]i=1,2,L,N,通过最优矩阵维度参数将梁模态振型信号生成迹矩阵S:
进一步设置是所述的步骤S3中对最优粒子逐维算法采用重心反向学习的变异方法,降低维间干扰和提高反向学习能扩大搜索空间,收敛精度。
具体算法如下:
S3.1按照随机概率分布初始化多维的粒子形成初始种群;
S3.2计算当前种群中每个粒子的适应度来确定个体最优和全体最优个体;
S3.3.根据式更新种群中每个粒子的速度向量及位置向量;
S3.4计算当前种群中每个粒子的适应度函数值。
进一步设置是所述的步骤S4中,建立损伤结构BSWI模型中采用任意尺度j的m阶的尺度函数以及其相应小波函数在0、1边界上有m-1个边界尺度函数与小波函数,2j-m+1个内部尺度函数,以及2j-2m+2个内部小波函数。
根据以下公式求出:
其中j0为满足2j≥2m-1的尺度,对于任意j>j0,令式(2)和(3)中的l=0,可以求出任意尺度j的尺度函数和小波函数由上述公式可知,在0、1边界上有m-1个边界尺度函数与小波函数,2j-m+1个内部尺度函数,以及2j-2m+2个内部小波函数。
进一步设置是所述的步骤S5中模糊系统将损伤结构BSWI模型计算结果进行二维空间插值,然后将测量频率作为输入智能搜索损伤程度诊断据库寻找最近匹配的损伤程度,实现梁损伤深度的定量识别。
本发明的有益效果是:
本发明能够克服单一损伤特征提取方法对损伤敏感性不足的问题,本发明具有抗噪声干扰能力强和自适应优化奇异值分解矩阵维度参数的特点,能够对损伤机械结构进行快速精确定位损伤位置和深度。本发明通过自适应奇异值分解的损伤识别方法实现获取损伤结构的模态信息、自适应优化迹矩阵维度得到奇异值分析损伤位置、综合利用小波分解的多分辨分析特性建立损伤数据库计算损伤深度。该自适应损伤识别方法对损伤的敏感性得到了相应的增强,提高了损伤识别能力,可以有效地对结构损伤进行定位识别。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,根据这些附图获得其他的附图仍属于本发明的范畴。
图1本发明的放大的步骤运行图;
图2是本发明的流程PN图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
如图1至图2所示,为本发明实施例中,包括以下步骤:
步骤S1中,通过激光位移传感器沿着梁的轴向方向遍历测量获取梁的一维振动位移,利用经验模态分解提取损伤梁的振型和频率特征信号。可以得到无损伤结构的质量矩阵M与刚度矩阵K和损伤结构的质量矩阵与刚度矩阵接下来先对无损伤结构进行模态分析,在此基础上采用摄动法对损伤结构进行模态分析。构造方程的解x如下:x=Φ·aejωt,其中Φ和a是常量。当Φ不为零时,这转化为一个广义特征值问题,ω2是特征值,Φ是特征向量。则系统特征方程
|K-ω2M|=0
解此特征方程可以得到n个特征值以及与每个特征值相对应的特征向量Φi(i=1,2,L,n),在结构模态分析中,与Φi分别为结构的第i阶自然频率和与其相对应的模态振型。将n个特征值按从小到大的次序排成一个对角矩阵,同时将Φi也根据的次序按列排成一个矩阵:
式中:Φ称为系统特征向量矩阵,在模态分析又称为模态矩阵或模态振型;Λ称为谱矩阵。可以证明以下等式成立:
ΦTMΦ=diag[mi],ΦTKΦ=diag[ki]
式中:mi与ki分别为结构的模态质量与模态质量;上标T表示矩阵的转置;diag表示对角矩阵。可以得到:
模态质量mi与模态刚度ki均与Φi大小有关,而Φi中的各个元素大小不定,而比例是固定的。
步骤S2中,将前几阶振型信号生成迹矩阵,利用最优粒子逐维变异的粒子群算法优化最佳自适应奇异值分解矩阵维度参数是:
a.利用最优粒子逐维变异的粒子群算法优化最佳自适应奇异值分解矩阵维度参数;
b.计算任意一阶模态信息生成列向量X=[x(1)x(2)L x(i)L x(N)]i=1,2,L,N,通过最优矩阵维度参数将梁模态振型信号生成迹矩阵S。
步骤S3中:对最优粒子逐维算法采用重心反向学习的变异方法,降低维间干扰和提高反向学习能扩大搜索空间,具体算法如下:
a.按照随机概率分布初始化多维的粒子形成初始种群;
b.计算当前种群中每个粒子的适应度来确定个体最优和全体最优个体;
c.根据式更新种群中每个粒子的速度向量及位置向量;
d.计算当前种群中每个粒子的适应度函数值。
其中j0为满足2j≥2m-1的尺度,对于任意j>j0,令式(2)和(3)中的l=0,可以求出任意尺度j的尺度函数和小波函数由上述公式可知,在0、1边界上有m-1个边界尺度函数与小波函数,2j-m+1个内部尺度函数,以及2j-2m+2个内部小波函数。
步骤S5中模糊系统将损伤结构BSWI模型计算结果进行二维空间插值,然后将测量频率作为输入智能搜索损伤程度诊断据库寻找最近匹配的损伤程度,实现梁损伤深度的定量识别。
以上所揭露的仅为本发明较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。
Claims (6)
1.一种基于自适应奇异值分解的梁损伤识别方法,其特征在于包括以下步骤:
S1:通过激光位移传感器获取损伤梁的振动位移,利用经验模态分解提取损伤梁的模态振型和模态频率特征信号;
S2:利用最优粒子逐维算法自动优化奇异值分解矩阵维度参数,然后将梁模态振型信号生成迹矩阵;
S3:通过将提取的频率特征信号画在损伤结构的空间全域上,定义适度函数确定峰值点的位置定位损伤的位置;
S4:通过基于识别的位置建立损伤结构BSWI模型计算损伤程度不同损伤程度的固有频率,建立固有频率与损伤程度对应关系数据库;
S5:将测量频率作为模糊系统的输入,通过智能搜索损伤程度诊断据库寻找最近匹配的损伤程度,实现梁损伤深度的定量识别。
2.根据权利要求1所述的一种基于自适应奇异值分解的梁损伤识别方法,其特征在于:所述的步骤S1中,通过激光位移传感器获取损伤梁的振动位移,利用经验模态分解提取损伤梁的模态振型和频率特征信号。
4.根据权利要求1所述的一种基于自适应奇异值分解的梁损伤识别方法,其特征在于:所述的步骤S3中对最优粒子逐维算法采用重心反向学习的变异方法,降低维间干扰和提高反向学习能扩大搜索空间,收敛精度。
6.根据权利要求1所述的一种基于自适应奇异值分解的梁损伤识别方法,其特征在于:所述的步骤S5中模糊系统将损伤结构BSWI模型计算结果进行二维空间插值,然后将测量频率作为输入智能搜索损伤程度诊断据库寻找最近匹配的损伤程度,实现梁损伤深度的定量识别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110487297.9A CN113155973B (zh) | 2021-05-05 | 2021-05-05 | 一种基于自适应奇异值分解的梁损伤识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110487297.9A CN113155973B (zh) | 2021-05-05 | 2021-05-05 | 一种基于自适应奇异值分解的梁损伤识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113155973A true CN113155973A (zh) | 2021-07-23 |
CN113155973B CN113155973B (zh) | 2023-05-16 |
Family
ID=76873258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110487297.9A Active CN113155973B (zh) | 2021-05-05 | 2021-05-05 | 一种基于自适应奇异值分解的梁损伤识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113155973B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113945633A (zh) * | 2021-09-30 | 2022-01-18 | 西安交通大学 | 一种基于高空间分辨率振型分型维数的结构微弱损伤识别方法 |
CN113947569A (zh) * | 2021-09-30 | 2022-01-18 | 西安交通大学 | 一种基于计算机视觉的梁型结构多尺度微弱损伤定位方法 |
CN115655887A (zh) * | 2022-11-01 | 2023-01-31 | 广东建设职业技术学院 | 一种混凝土强度的预测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6779404B1 (en) * | 1999-11-03 | 2004-08-24 | Rune Brincker | Method for vibration analysis |
CN104807534A (zh) * | 2015-05-21 | 2015-07-29 | 华北电力大学(保定) | 基于在线振动数据的设备固有振动模式自学习识别方法 |
CN104849348A (zh) * | 2015-05-13 | 2015-08-19 | 河海大学 | 一种基于奇异值分解-小波变换的梁结构损伤检测方法 |
CN107300451A (zh) * | 2017-01-09 | 2017-10-27 | 温州大学 | 一种基于损伤梁固有频率快速估算的检测方法 |
CN108802176A (zh) * | 2018-04-08 | 2018-11-13 | 大连理工大学 | 一种基于pvdf传感器与应变模态的结构损伤识别实验方法 |
CN109145464A (zh) * | 2018-08-28 | 2019-01-04 | 暨南大学 | 融合多目标蚁狮优化与迹稀疏正则化的结构损伤识别方法 |
WO2020156348A1 (zh) * | 2019-01-31 | 2020-08-06 | 青岛理工大学 | 基于总体经验模态分解与卷积神经网络的结构损伤识别法 |
-
2021
- 2021-05-05 CN CN202110487297.9A patent/CN113155973B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6779404B1 (en) * | 1999-11-03 | 2004-08-24 | Rune Brincker | Method for vibration analysis |
CN104849348A (zh) * | 2015-05-13 | 2015-08-19 | 河海大学 | 一种基于奇异值分解-小波变换的梁结构损伤检测方法 |
CN104807534A (zh) * | 2015-05-21 | 2015-07-29 | 华北电力大学(保定) | 基于在线振动数据的设备固有振动模式自学习识别方法 |
CN107300451A (zh) * | 2017-01-09 | 2017-10-27 | 温州大学 | 一种基于损伤梁固有频率快速估算的检测方法 |
CN108802176A (zh) * | 2018-04-08 | 2018-11-13 | 大连理工大学 | 一种基于pvdf传感器与应变模态的结构损伤识别实验方法 |
CN109145464A (zh) * | 2018-08-28 | 2019-01-04 | 暨南大学 | 融合多目标蚁狮优化与迹稀疏正则化的结构损伤识别方法 |
WO2020156348A1 (zh) * | 2019-01-31 | 2020-08-06 | 青岛理工大学 | 基于总体经验模态分解与卷积神经网络的结构损伤识别法 |
Non-Patent Citations (3)
Title |
---|
JIAWEI XIANG等: "A two-step approach to multi-damage detection for plate structures", 《ENGINEERING FRACTURE MECHANICS》 * |
YONGYING JIANG等: "A two-step damage quantitative identification method for beam structures", 《MEASUREMENT》 * |
邱飞力等: "基于模态柔度矩阵的结构损伤识别", 《噪声与振动控制》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113945633A (zh) * | 2021-09-30 | 2022-01-18 | 西安交通大学 | 一种基于高空间分辨率振型分型维数的结构微弱损伤识别方法 |
CN113947569A (zh) * | 2021-09-30 | 2022-01-18 | 西安交通大学 | 一种基于计算机视觉的梁型结构多尺度微弱损伤定位方法 |
CN113945633B (zh) * | 2021-09-30 | 2022-10-28 | 西安交通大学 | 一种基于高空间分辨率振型分型维数的结构微弱损伤识别方法 |
CN113947569B (zh) * | 2021-09-30 | 2023-10-27 | 西安交通大学 | 一种基于计算机视觉的梁型结构多尺度微弱损伤定位方法 |
CN115655887A (zh) * | 2022-11-01 | 2023-01-31 | 广东建设职业技术学院 | 一种混凝土强度的预测方法 |
CN115655887B (zh) * | 2022-11-01 | 2023-04-21 | 广东建设职业技术学院 | 一种混凝土强度的预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113155973B (zh) | 2023-05-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113155973A (zh) | 一种基于自适应奇异值分解的梁损伤识别方法 | |
CN109655799B (zh) | 基于iaa的协方差矩阵向量化的非均匀稀疏阵列测向方法 | |
WO2016091017A1 (zh) | 一种高光谱图像分类中光谱向量互相关特征的抽取方法 | |
CN101609550A (zh) | 基于流形子空间的图像配准方法 | |
CN103648106A (zh) | 一种基于类别匹配的半监督流形学习的WiFi室内定位方法 | |
CN107656239B (zh) | 一种基于极化敏感阵列的相干信源测向方法 | |
CN107907855A (zh) | 一种互素阵列转化为均匀线阵的doa估计方法及装置 | |
CN110009745B (zh) | 根据平面基元和模型驱动对点云提取平面的方法 | |
CN107743296A (zh) | 一种基于压缩感知的rssi区域分段式定位方法 | |
CN112629786A (zh) | 一种工作模态参数识别方法及设备故障诊断方法 | |
Yan et al. | Computationally efficient direction finding using polynomial rooting with reduced-order and real-valued computations | |
CN115235769A (zh) | 一种基于马氏距离补偿因子的风机轴承故障诊断方法 | |
CN102609944A (zh) | 基于距离几何理论的高光谱遥感图像混合像元分解方法 | |
CN117119377A (zh) | 基于滤波Transformer的室内指纹定位方法 | |
CN103294647A (zh) | 基于正交张量近邻保持嵌入的头相关传输函数降维方法 | |
CN104407319A (zh) | 阵列信号的目标源测向方法和系统 | |
CN103258134A (zh) | 一种高维的振动信号的降维处理方法 | |
CN107909099A (zh) | 一种基于热核的三维模型识别与检索方法 | |
CN112597796A (zh) | 一种鲁棒的基于图卷积的点云表征学习方法 | |
CN116401528A (zh) | 一种基于全局-局部散度的多元时间序列无监督降维方法 | |
CN112528849B (zh) | 一种基于内积矩阵及深度学习的结构健康监测方法 | |
CN116486185A (zh) | 油井工况识别方法、装置、设备和存储介质 | |
CN103605880A (zh) | 一种精确诊断密集模态阻尼比的方法 | |
CN114299079A (zh) | 一种面向稠密点云数据的发动机叶片截面线数据获取方法 | |
CN118468093B (zh) | 基于Morse的欠定时不变工作模态参数识别方法及装置 |
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 |