CN102402647B - 一种在窄频范围内分析动力信号预测结构变形的方法 - Google Patents
一种在窄频范围内分析动力信号预测结构变形的方法 Download PDFInfo
- Publication number
- CN102402647B CN102402647B CN201110306294.7A CN201110306294A CN102402647B CN 102402647 B CN102402647 B CN 102402647B CN 201110306294 A CN201110306294 A CN 201110306294A CN 102402647 B CN102402647 B CN 102402647B
- Authority
- CN
- China
- Prior art keywords
- frequency
- frequency domain
- modal
- transfer function
- parameter
- 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.)
- Expired - Fee Related
Links
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种在窄频范围内分析动力信号预测结构变形的方法,通过土木结构冲击振动测试观测结构反应动力信号,经傅里叶变换和谱分析后得到结构的频域传递函数;将所述频域传递函数的整个频率空间划分为2-5个连续的狭窄子空间;对各子空间进行模态参数识别;根据模态叠加理论,对各子空间的识别结果进行叠加,计算出结构模态参数并进一步计算出整体的柔度矩阵。本发明首先将整个频域分割为数个窄频区域,然后在每个窄频区域分析实验数据,最后综合各窄频区域的分析结果识别结构的柔度矩阵,从而可预测结构在任何静力荷载下的结构变形。具有计算时间少、预测结果精度高的特点,适用于对各种工程结构,尤其是土木结构的健康诊断与安全评估。
Description
技术领域
本发明涉及一种在窄频范围内分析结构振动测试数据的方法,可实现土木工程结构的变形预测和健康诊断。
背景技术
振动测试为现有土木结构健康诊断与安全评估的主要方法。通过布置加速度计和其它类型传感器到土木结构,进行振动测试然后分析测试数据可识别土木结构主要参数从而可进行结构的健康诊断。研究及工程实践表明,由于传感器灵敏度和实验环境影响等原因,所测振动数据包含各种噪音,如何从信噪比差的测试信号中精确分析结构特征为结构健康诊断工作中的一项重要内容。另外,现有基于加速度信号的健康诊断方法大多仅输出结构的模态参数(频率、阻尼、阵型),而不能直接进行结构的变形预测。这些原因限制了结构健康诊断技术在土木工程中的应用,因此,开发一种稳定可靠的动力信号分析方法,从振动测试所得的动力信号中识别结构详细结构特征(如柔度矩阵),可以预测结构的变形,将促进健康诊断技术在工程实际中的应用。
发明内容
本发明的目的是提供一种用于土木工程结构健康诊断的动力信号分析方法,通过将测试动力信号的整个频域范围分割为数个窄频范围,实现结构详细特征(如柔度矩阵)的快速和准确识别,从而实现利用振动测试和信号分析预测结构变形。该方法计算快,精度高,抗噪音能力强,稳定性和重复性良好。
本发明的技术解决方案如下:
一种用于土木工程结构健康诊断的动力信号分析方法:第一步、通过土木结构冲击振动测试观测结构反应动力信号,经傅里叶变换和谱分析后得到结构的频域传递函数;第二步、将所述频域传递函数的整个频率空间划分为2-5个连续的狭窄子空间;第三步、对各子空间进行参数识别;第四步、根据模态叠加理论,对各子空间的识别结果进行叠加,计算出结构模态参数并进一步计算出整体的柔度矩阵。
所述的模态参数识别方法是:
Jsθs=0 (1)
在公式(1)中, ns为一用户定义的参数.其中,θs包含未知参数, 为由测量数据得到的结构频域--函数,ωka表示频率,Js由已知的频域传递函数数据组成,未知参数θs可由最小二平方方法计算得出,θs得出后进行特征值分析可得到结构的频率阻尼系数,振型参与系数和振型,γ,L和 各窄频范围内的模态参数得出后进行叠加可得出整体结构的模态参数矩阵,柔度矩阵方程为:
公式2中f为柔度矩阵, 为残余项,符号*表示相应变量的共轭复数。与现有技术相比,本发明具有的有益效果是:计算速度快(只利用一个窄频范围内的数据),识别结果精度高(各窄频范围范围内的计算不互相影响)。特别的是,本发明通过在窄频范围内动力信号分析,可识别出结构的柔度矩阵,从而可进行结构的变形预测,这是绝大多数现有基于振动的结构健康诊断方法做不到的。
附图说明
图1是利用本发明在窄频范围识别结构参数的结果图。
图2是利用本发明结果预测的结构变形和实际结果比较图。
具体实施方式
以下结合附图对本发明作详细说明,但不应以此限制本发明的保护范围。
通过土木结构冲击振动测试观测结构反应动力信号,经傅里叶变换和谱分析后得到结构的频域传递函数,该频域传感函数为一系列分布在宽幅频域范围(如0.1-100赫兹)的数据。常规处理方法直接利用整个宽频范围内的频域传递函数数据进行信号分析以识别结 构模态参数(如PolyMAX方法),导致运算时间长,识别的模态参数不精确(特别是低幅值频域传递函数对应的高频区域)。本发明方法将整个频率空间划分为数个连续的狭窄频率子空间(如图1所示),每个子空间(一个窄频范围如0.1-5赫兹)仅包含数量有限的频域传递函数数据。然后在各个窄频范围内进行如下所示的参数识别:
Jsθs=0 (1)
在公式(1)中, ns为一用户定义的参数.其中,θs包含未知参数,Js由已知的频域传递函数数据组成。未知参数θs可由最小二平方方法计算得出。值得指出的是公式(1)是建立在分析窄频范围内动力信号的基础上,因此矩阵Js仅包含一个窄频范围内的数量很少的数据,所需的ns参数又小,从而保证了矩阵尺寸小,计算速度快。另一方面,各个窄频范围内未知参数的计算是独立的,从而保证了在计算低幅值频域传递函数对应的结构模态时,其不受其它频带内高幅值频域传递函数数据的影响,从而可得出精确的识别结果。各频域范围内的未知参数通过公式(1)求出后,根据模态叠加理论进行综合,可计算出结构模态参数(频率、阻尼和阵型)并进一步计算出柔度矩阵:
公式2中f为柔度矩阵, 为残余项,γ,L和 为由前面所述的数据分析方法得出的频率阻尼系数、阵型参与系数和阵型。根据公式2识别出的结构柔度矩阵可计算结构在任何静力荷载作用下的结构变形。工程实践中已广泛接受利用结构变形来进行结构健 康诊断的方法,所以本发明通过在窄频范围内进行动力信号分析和结构变形预测,有望在实际土木工程的健康诊断和评估中得到广泛应用和推广。
在土木工程结构的健康监测与诊断中,首先通过布置各类型传感器(如加速度计)到所测试桥梁上进行环境振动或冲击振动测试。经过数据预处理、傅里叶变换和谱分析后,可得到整个宽频范围内的频域传递函数数据,如图1所示(该例中整个频域范围为4-32赫兹)。利用本发明方法,将该宽频范围划分为两个窄频范围(4-18赫兹和18-32赫兹),随后在每个窄频范围内利用公式(1)识别未知参数。例如图1显示了在两个窄频范围内真实的频率传递函数和由识别的未知参数计算复原得出的频率传递函数。可以看出他们相当吻合。随后,根据模态理论融合在两个窄频范围内识别的未知参数,识别结构的模态参数(频率、阻尼和阵型)尤其是结构的阵型缩放系数和柔度矩阵。得到结构的柔度矩阵,并进而可预测结构在任何静力荷载下的结构变形。例如图2显示了在一静力荷载下利用本发明方法预测的结构变形以及和真实值的比较,可以发现当利用足够多的模态数后(该例中为前六阶模态),所预测的结构变形和真实值基本一致。
由图1可以看出,在划分的两个窄频范围内,结构的频率传递函数峰值差别很大。当利用传统方法利用整个频域范围内的数据进行结构识别是,低幅值频域传递函数对应的频率区域(如本例中的18-32赫兹区域)中的模态参数就不易精确识别。而本发明方法在各个窄频范围内识别结构参数,各个窄频范围内的计算互不影响,所以低幅值频域传递函数对应的频率区域内的模态参数也能精确识别。另外,本发明方法可精确预测结构反应,可有效用来进行土木结构的健康诊断和安全评估。为本发明的独特之处。
Claims (1)
1.一种在窄频范围内分析动力信号预测结构变形的方法,其特征在于:方法如下:第一步、通过土木结构冲击振动测试观测结构反应动力信号,经傅里叶变换和谱分析后得到结构的频域传递函数;第二步、将所述频域传递函数的整个频率空间划分为2-5个连续的狭窄子空间;第三步、对各子空间进行模态参数识别;第四步、根据模态叠加理论,对各子空间的识别结果进行叠加,计算出结构模态参数并进一步计算出整体的柔度矩阵;
所述的模态参数识别方法是:
Jsθs=0 (1)
在公式(1)中, n2为一用户定义的参数,其中,θ3包含未知参数,为由测量数据得到的结构频域函数,ωka表示频率,Js由已知的频域传递函数数据组成,未知参数θs可由最小二平方方法计算得出,θs得出后进行特征值分析可得到结构的频率阻尼系数,振型参与系数和振型,γ,L和各窄频范围内的模态参数得出后进行叠加可得出整体结构的模态参数矩阵,柔度矩阵方程为:
其中f为柔度矩阵,为残余项,符号*表示相应变量的共轭复数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110306294.7A CN102402647B (zh) | 2011-10-11 | 2011-10-11 | 一种在窄频范围内分析动力信号预测结构变形的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110306294.7A CN102402647B (zh) | 2011-10-11 | 2011-10-11 | 一种在窄频范围内分析动力信号预测结构变形的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102402647A CN102402647A (zh) | 2012-04-04 |
CN102402647B true CN102402647B (zh) | 2014-10-08 |
Family
ID=45884846
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110306294.7A Expired - Fee Related CN102402647B (zh) | 2011-10-11 | 2011-10-11 | 一种在窄频范围内分析动力信号预测结构变形的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102402647B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106950047B (zh) * | 2017-03-09 | 2018-12-18 | 西安交通大学 | 振动加速度信号频谱的可视化分析方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101266190A (zh) * | 2008-04-25 | 2008-09-17 | 天津市市政工程设计研究院 | 正常交通流下斜拉桥模态参数测定装置及方法 |
CN101487763B (zh) * | 2009-02-23 | 2010-12-08 | 西北工业大学 | 一种测量大噪声环境下振动结构频率响应函数的方法 |
-
2011
- 2011-10-11 CN CN201110306294.7A patent/CN102402647B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102402647A (zh) | 2012-04-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Amezquita-Sanchez et al. | Signal processing techniques for vibration-based health monitoring of smart structures | |
CN102353509B (zh) | 一种基于分块冲击振动测试的中小型桥梁快速检测方法 | |
Zhang et al. | Estimation of mode shapes of beam-like structures by a moving lumped mass | |
CN102928514B (zh) | 一种基于频率特征的木材应力波无损检测方法 | |
Dong et al. | Structural damage detection using empirical-mode decomposition and vector autoregressive moving average model | |
US20200089730A1 (en) | Modal identification method for non-proportionally damped structures based on extended sparse component analysis | |
Garcia-Pozuelo et al. | Development and experimental validation of a real-time analytical model for different intelligent tyre concepts | |
Nie et al. | Bridge condition monitoring under moving loads using two sensor measurements | |
Gillich et al. | A method to enhance frequency readability for early damage detection | |
CN104112072A (zh) | 基于小波阈值去噪的主成分分析的工作模态参数识别方法 | |
CN104406757B (zh) | 适用于中小桥梁快速安全诊断的一体化装置 | |
Mao et al. | The construction and comparison of damage detection index based on the nonlinear output frequency response function and experimental analysis | |
US20190376874A1 (en) | A method of estimating the number of modes for the sparse component analysis based modal identification | |
CN109839441A (zh) | 一种桥梁模态参数识别方法 | |
CN104165742A (zh) | 一种基于互谱函数的运行模态分析实验方法及装置 | |
CN105862935A (zh) | 一种用于挡土墙结构系统的损伤识别方法 | |
Li et al. | Covariance of dynamic strain responses for structural damage detection | |
Li et al. | Structural damage identification based on short-time temporal coherence using free-vibration response signals | |
CN110057918B (zh) | 强噪声背景下的复合材料损伤定量识别方法及系统 | |
CN103543026A (zh) | 基于振动传递率函数和支持向量机的结构损伤识别方法 | |
CN102402647B (zh) | 一种在窄频范围内分析动力信号预测结构变形的方法 | |
CN107145620B (zh) | 一种基于随机减量技术的结构动力特性识别方法 | |
CN106679911A (zh) | 一种基于多尺度数据融合理论的梁型结构损伤识别方法 | |
CN116522085A (zh) | 一种全自动拉索频率提取、定阶和索力识别方法与应用 | |
Abdullah et al. | Fatigue features extraction of road load time data using the S-transform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141008 Termination date: 20181011 |
|
CF01 | Termination of patent right due to non-payment of annual fee |