CN100464328C - 改进希-黄变换的结构响应分析方法 - Google Patents
改进希-黄变换的结构响应分析方法 Download PDFInfo
- Publication number
- CN100464328C CN100464328C CNB2006100468977A CN200610046897A CN100464328C CN 100464328 C CN100464328 C CN 100464328C CN B2006100468977 A CNB2006100468977 A CN B2006100468977A CN 200610046897 A CN200610046897 A CN 200610046897A CN 100464328 C CN100464328 C CN 100464328C
- Authority
- CN
- China
- Prior art keywords
- intrinsic mode
- mode functions
- original signal
- response
- signal
- 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
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
改进希-黄变换的结构响应分析方法，涉及一种非稳态结构响应信号处理方法。该方法先将响应信号进行带通滤波，接着进行经验模态分解；根据本征模函数和原信号的相关系数m判定真正本征模函数；然后采用随机减量技术将每个本征模函数转换为自由振动响应，最后运用希尔伯特变换和最小二乘拟合方法提取结构响应的频率、阻尼比模态参数。该方法能够较好地处理多自由度、非自由振动、非线性、非稳态的响应信号处理技术与模态参数识别，可用于土木工程、航空航天、自动控制、机械工程、水利工程等领域的信号处理及故障诊断，具有简单、对高阶模态识别更加准确的特点。
Description
技术领域
本发明涉及一种非稳态结构响应信号处理方法，特别涉及一种改进希-黄变换的结构响应分析方法。
背景技术
随着地震、飓风等自然灾害的频繁发生与构筑物使用性能的下降，各种工程事故如房屋倒塌、桥梁突然出现裂缝乃至断裂等频繁发生，对国民经济及人民生命财产造成了极大地损失。如果能够对这些构筑物提前进行检测、诊断，以使相关管理部门提前做出维修、加固决策，从而将损失减少到最低。可见，重大工程的安全性与可靠性判定显得日益重要。
结构损伤检测较早采用的方法是结构静态检测，通过直接测量结构及构件的尺寸、强度和弹性模量(回弹法、钻芯法等)，进行结构分析以确定结构的工作性能与可靠性水平的方法。该方法的优点是测量结果直接且比较可靠，缺点是工作量巨大、工作效率低，受到应用条件限制。其后，基于振动的结构动力无损检测方法得到了重视。它利用结构的动力响应进行结构性能状态识别，它不受结构规模和构件隐蔽的限制，高效率模块化、数字化的结构动力响应量测技术为结构动力检测方法提供了有效的技术支持。动力检测方法的特点是应用条件限制少，效率高，既可以对复杂结构的整体也可以对结构的局部进行检测，对构筑物无破坏和无环境污染的优点而受到国内外重视，是国际社会研究的热点和前沿课题。
结构动力检测方法的基础和重点是信号分析，当前结构响应信号处理技术集中在线性结构和线性响应上面，实际上复杂结构及其响应通常是非线性、非稳态的，结构的冗余性和实测响应的噪声与误差是结构损伤检测与信号处理另一个技术难点。常用的时-频域信号处理技术包括(短时)快速傅立叶变换(FFT)、小波分析及近几年出现的希-黄变换(HHT)。其中，(短时)FFT实际上是加窗的FFT，其分析精度由滑移窗的时窗和谱窗的宽度联合确定，窗宽越窄，则相应的分辨率越高。小波(包)变换不但继承和发展了短时FFT的局部思想，而且克服了窗口大小不随频率变化、缺乏离散正交基的缺点，但小波(包)分析需要经验来确定小波基函数类型，进而进行小波(包)分解与处理。HHT是近年新提出的一种信号处理方法，它由经验模态分解(EMD)方法及希尔伯特变换(HT)两部分组成，其核心是EMD分解，为非线性和非稳态信号的分析提供了一种十分有效的方法HHT。虽然HHT对于分析广泛存在的非线性与非平稳信号十分有效，但也存在一些问题。第一，EMD方法在低频域会产生一些我们并不希望的低幅值的IMFs并且会产生一些不应该有的频率组分。第二，在高频域第一个产生的IMF可能会覆盖一个很广泛的频率范围，所以不能很好的满足单一组分的定义。第三，EMD过程经常不能够从被分析信号中区分一些低能量组分，所以这些组分在时频平面内就可能不会出现。此外，它只能处理自由振动响应。
传统HHT在EMD过程中存在高频域第一个产生的IMF可能会覆盖很广泛的频率范围而不能较好地满足单一组分的定义、在低频域会产生一些不希望的低幅值的IMFs且会产生不该有的频率组分、不能从信号中区分一些低能量组分以及只能处理自由振动响应等问题。
发明内容
本发明的目的在于提出一种改进希-黄变换的结构响应分析方法，该方法能够较好地处理多自由度、非自由振动、非线性、非稳态的响应信号处理技术与模态参数识别，具有简单、对高阶模态识别更加准确的特点。
本发明的目的是通过以下技术方案实现的：
改进希-黄变换的结构响应分析方法，先将响应信号进行带通滤波，接着进行经验模态分解；根据本征模函数和原信号的相关系数m判定真正本征模函数；然后采用随机减量技术将每个本征模函数转换为自由振动响应，最后运用希尔伯特变换和最小二乘拟合方法提取结构响应的频率、阻尼比。
一种改进希-黄变换的结构响应分析方法，该方法包括以下步骤：
(1)先将原信号根据快速傅立叶变换得出每一个固有频率的范围ωjL<ωj<ωjH(j＝1，2，...，n)，使原信号通过每一个以ωjL<ωj<ωjH作为频率带的带通滤波器；
(2)对通过第j个带通滤波器过滤得到的时间信号进行经验模态分解，获得本征模函数，对j＝1，2，...，n重复进行同样的过程，就可以获得n个所要求的本征模函数，它很容易地就分离出本征模函数并且可以将频率ωjL<ωj<ωjH以外的所有噪声全部去除；
(3)计算本征模函数和原信号的相关系数，以判定真正低频本征模函数组分，首先将所有的本征模函数和原信号归一化，计算所有本征模函数与原信号的相关系数covef(i)，将相关系数与初始值λ相比较，如果covef(i)≥λ，则作为本征模函数，否则将该本征模函数删除；
这里的初始值λ取
λ＝max(covef(i))/β (i＝1，2，...，n)
其中，β通常取1～20；
(4)确定本征模函数后，采用随机减量技术将每个本征模函数转换为自由振动响应；
(5)将每个自由振动响应进行希尔伯特变换，得到瞬时相位角和衰减幅值；
(6)使用最小二乘拟合方法对瞬时相位角和衰减幅值进行估计，得到固有频率和阻尼比。
本发明的优点与效果是：
1.本发明技术合理简单，适合处理非线性、非稳态的结构响应；
2.本发明可以较好地解决传统HHT在低频域错分IMF和频率，以及在高频域难以满足单一组分的要求；
3.本发明既可以用于自由振动结构响应信号，也可以用于受迫振动等外部或环境荷载作用下的响应信号处理。
4.本发明可以直接得到结构的频率和阻尼比，不需要分别独立计算。
附图说明
图1为本发明的原信号图；
图2为本发明的信号分解得到的IMFs示意图；
图3为本发明加速度信号时程记录示意图；
图4为本发明经过滤波和相关系数计算后的前5～7阶IMFs示意图；
图5为本发明平稳响应IMFs分量的自由衰减响应曲线图。
具体实施方式
下面参照附图对本发明进行详细说明。
实施例1
图1所示为一个多组分自由振动响应信号，它包含两个频率组分。图2给出了这个信号通过传统HHT的EMD后得到的4个IMFs。很明显，只有前两个IMFs是信号真正的组分，其余的两个是假组分，它们是低频的并在Hilbert谱中绘成低频组分来误导我们的分析。
经过本发明计算，在图2中只得到前2个IMFs，后2个本征模函数(IMF)被删除掉只是作为残余分量。
实施例2
某7层、2×1跨的钢结构框架，经过受迫振动后实测的加速度响应信号(图3)，采用本发明方法进行信号分析。
首先应用FFT对此加速度响应信号做快速傅立叶变换，来进行频率识别，然后，使用改进的HHT方法对原时程信号进行EMD，经过滤波和相关系数计算，得到前5～7阶IMF(本征模函数)如图4所示。
图5示出了由平稳响应分量IMFs经随机减量技术(RDT)法得到的衰减曲线，其截取阀值取为原始信号标准差的1.3倍，衰减长度取为0.512s，平均次数为197次。经最小二乘拟合得到频率值为58.1249Hz和阻尼比ξ＝0.0822％。
Claims (1)
1.一种改进希-黄变换的结构响应分析方法，其特征在于，该方法包括以下步骤：
(1)先将原信号根据快速傅立叶变换得出每一个固有频率的范围ωjL<ωj<ωjH(j＝1，2，...，n)，使原信号通过每一个以ωjL<ωj<ωjH作为频率带的带通滤波器；
(2)对通过第j个带通滤波器过滤得到的时间信号进行经验模态分解，获得本征模函数，对j＝1，2，...，n重复进行同样的过程，就可以获得n个所要求的本征模函数，它很容易地就分离出本征模函数并且可以将频率ωjL<ωj<ωjH以外的所有噪声全部去除；
(3)计算本征模函数和原信号的相关系数，以判定真正低频本征模函数组分，首先将所有的本征模函数和原信号归一化，计算所有本征模函数与原信号的相关系数covef(i)，将相关系数与初始值λ相比较，如果covef(i)≥λ，则作为本征模函数，否则将该本征模函数删除；
这里的初始值λ取
λ＝max(covef(i))/β (i＝1，2，...，n)
其中，β通常取1～20；
(4)确定本征模函数后，采用随机减量技术将每个本征模函数转换为自由振动响应；
(5)将每个自由振动响应进行希尔伯特变换，得到瞬时相位角和衰减幅值；
(6)使用最小二乘拟合方法对瞬时相位角和衰减幅值进行估计，得到固有频率和阻尼比。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2006100468977A CN100464328C (zh) | 2006-06-15 | 2006-06-15 | 改进希-黄变换的结构响应分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2006100468977A CN100464328C (zh) | 2006-06-15 | 2006-06-15 | 改进希-黄变换的结构响应分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1869972A CN1869972A (zh) | 2006-11-29 |
| CN100464328C true CN100464328C (zh) | 2009-02-25 |
Family
ID=37443639
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2006100468977A Expired - Fee Related CN100464328C (zh) | 2006-06-15 | 2006-06-15 | 改进希-黄变换的结构响应分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100464328C (zh) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101901209B (zh) * | 2010-02-04 | 2013-04-03 | 福州大学 | 基于改进emd和arma模型的结构响应分析方法 |
| BE1019445A3 (fr) * | 2010-08-11 | 2012-07-03 | Reza Yves | Procede d'extraction d'information audio. |
| CN101900708B (zh) * | 2010-08-18 | 2013-07-31 | 哈尔滨工业大学 | 一种基于振动和声频信号的高速列车轨道伤损探测方法 |
| CN102467654A (zh) * | 2010-11-03 | 2012-05-23 | 西安电子科技大学 | 一种结构模态参数识别方法 |
| CN102346064A (zh) * | 2011-09-08 | 2012-02-08 | 浙江工业大学 | 装载机称重装置 |
| CN102736561A (zh) * | 2012-07-10 | 2012-10-17 | 北京信息科技大学 | 一种面向机电装备的远程动态自适应规则获取方法 |
| CN102855408A (zh) * | 2012-09-18 | 2013-01-02 | 福州大学 | 基于ica的改进emd过程中imf判定方法 |
| CN103344448A (zh) * | 2013-06-26 | 2013-10-09 | 中国路桥工程有限责任公司 | 一种桥梁结构损伤识别方法和*** |
| CN104048677A (zh) * | 2014-07-03 | 2014-09-17 | 哈尔滨工业大学 | 基于k-s分布性检验和hht的陀螺仪故障诊断方法 |
| CN105587475B (zh) * | 2015-12-16 | 2018-12-21 | 北京金风科创风电设备有限公司 | 风力发电机组及其塔架***状态的检测方法和装置 |
| CN106092879B (zh) * | 2016-06-07 | 2019-07-12 | 西安向阳航天材料股份有限公司 | 基于振动响应信息的***复合管结合状态检测方法 |
| CN108956070B (zh) * | 2018-07-03 | 2020-01-31 | 武汉理工大学 | 一种基于磁流变弹性体变频阻尼器减振试验装置及其方法 |
| CN109214092A (zh) * | 2018-09-11 | 2019-01-15 | 吉林大学 | 基于希尔伯特黄变换的四分量钻孔应变数据异常提取方法 |
| CN113158785B (zh) * | 2021-03-11 | 2022-11-15 | 复旦大学 | 一种振荡信号模态参数的识别方法 |
| CN113009576B (zh) * | 2021-03-26 | 2022-04-12 | 中国石油大学(华东) | 一种基于本征模态滤波的储层固有频率原位非接触式检测方法 |
| CN114993451A (zh) * | 2022-06-23 | 2022-09-02 | 无锡中鼎集成技术有限公司 | 低频振动测试***及测量方法 |
-
2006
- 2006-06-15 CN CNB2006100468977A patent/CN100464328C/zh not_active Expired - Fee Related
Non-Patent Citations (8)
| Title |
|---|
| 基于主组分分析的概率神经网络损伤定位研究. 姜绍飞等.地震工程与工程振动，第24卷第2期. 2004 |
| 基于主组分分析的概率神经网络损伤定位研究. 姜绍飞等.地震工程与工程振动，第24卷第2期. 2004 * |
| 基于振动的结构健康监测技术. 姜绍飞等.沈阳建筑工程学院学报(自然科学版)，第19卷第4期. 2003 |
| 基于振动的结构健康监测技术. 姜绍飞等.沈阳建筑工程学院学报(自然科学版)，第19卷第4期. 2003 * |
| 小波函数的选择对结构损伤识别的影响. 杨晓楠等.沈阳建筑工程学院学报(自然科学版)，第21卷第6期. 2005 |
| 小波函数的选择对结构损伤识别的影响. 杨晓楠等.沈阳建筑工程学院学报(自然科学版)，第21卷第6期. 2005 * |
| 考虑不确定性因素的结构损伤检测方法. 姜绍飞等.沈阳建筑工程学院学报(自然科学版)，第18卷第2期. 2002 |
| 考虑不确定性因素的结构损伤检测方法. 姜绍飞等.沈阳建筑工程学院学报(自然科学版)，第18卷第2期. 2002 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1869972A (zh) | 2006-11-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100464328C (zh) | 改进希-黄变换的结构响应分析方法 | |
| McCormick et al. | Cyclostationarity in rotating machine vibrations | |
| Kijewski et al. | Wavelet transforms for system identification in civil engineering | |
| Mechefske et al. | Fault detection and diagnosis in low speed rolling element bearings Part I: The use of parametric spectra | |
| Yan et al. | A comparative study of modal parameter identification based on wavelet and Hilbert–Huang transforms | |
| CN110160765B (zh) | 一种基于声音或振动信号的冲击特征识别方法及*** | |
| CN104165925B (zh) | 随机共振的离心式压缩机半开式叶轮裂纹故障检测方法 | |
| CN106546892A (zh) | 基于深度学习的局部放电超声音频识别方法及*** | |
| CN106096242B (zh) | 一种基于改进emd分解的尾水管压力脉动综合评价方法 | |
| Wang et al. | An energy kurtosis demodulation technique for signal denoising and bearing fault detection | |
| Zhan et al. | Adaptive autoregressive modeling of non-stationary vibration signals under distinct gear states. Part 1: modeling | |
| Zhang et al. | Variational mode decomposition based modal parameter identification in civil engineering | |
| Ni et al. | Rolling element bearings fault diagnosis based on a novel optimal frequency band selection scheme | |
| CN101769787A (zh) | 电力变压器铁心松动故障监测方法 | |
| Dong et al. | A frequency-shifted bispectrum for rolling element bearing diagnosis | |
| Xiao et al. | Research on fault feature extraction method of rolling bearing based on NMD and wavelet threshold denoising | |
| Dias et al. | A study of spectra, structure and correlation functions and their implications for the stationarity of surface-layer turbulence | |
| Shi et al. | The VMD-scale space based hoyergram and its application in rolling bearing fault diagnosis | |
| Yang et al. | Distinguish coherence resonance and stochastic resonance in bearing fault evaluation | |
| CN101251445A (zh) | 旋转机械碰摩声发射信号的分形特征分析方法 | |
| Jiang et al. | Differential spectral amplitude modulation and its applications in rolling bearing fault diagnosis | |
| Sun et al. | Modal identification from non-stationary responses of high-rise buildings by variational mode decomposition and direct interpolation techniques | |
| Hua et al. | The methodology of modified frequency band envelope kurtosis for bearing fault diagnosis | |
| Chan et al. | A novel, fast, reliable data transmission algorithm for wireless machine health monitoring | |
| Zheng et al. | Zero-phase filter-based adaptive Fourier decomposition and its application to fault diagnosis of rolling bearing |
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: 20090225 Termination date: 20150615 |
|
| EXPY | Termination of patent right or utility model |