专利名称：一种改进功率谱峰值法的空间网格结构频率识别方法
技术领域：
本发明涉及一种结构的频率识别方法，特别涉及一种改进功率谱峰值法的空间网格结构频率识别方法。
背景技术：
空间网格结构体系迅速发展于20世纪80年代，截至目前其结构形式已经向多样化、新颖化和超大跨度化方向发展。2008年北京奥运会、2009年i^一届全运会、2010年上海世博会及2011年的深圳大运会后，大量新颖独特的空间网格结构奥运场馆即将向公众开发并长期投入使用，如鸟巢、水立方、深圳大运会场馆、各大城市的奥体中心等。这些新建或已建的大跨度空间网格结构工程，连同上个世纪建造的大跨空间结构无疑构成了我国社会和经济发展的主要载体。大型空间网格结构工程设计使用年限长达上百年，在其服役过程中由于环境荷载作用、疲劳效应、腐蚀效应、材料老化和突变效应等灾害因素影响，结构不可避免地会产生损伤的累积和抗力的衰减，从而导致抵抗自然灾害能力的下降，极端情况下会引发灾难性的突发事故。这样的重大工程结构一旦失效，将会造成巨大的经济损失和人员伤亡，由此可见为了保障结构的安全性和可靠性，已建成使用的许多大跨、超大跨空间结构的损伤诊断与健康评估显得日益重要。结构损伤检测较早采用的方法是结构静态检测，即通过直接测量结构及构件的尺寸、强度和内部缺陷(如超声波法、发射光谱法、射线法等)，进行计算分析以确定结构的工作性能与可靠性水平。该方法的优点是测量结果比较可靠，缺点是工作量大、效率低，受到应用条件的限制，且仅能测量结构的局部。随着科技与经济的发展，建筑结构不断向高层、超高层或大跨超大跨方向发展，结构局部的检测已逐渐不能满足结构健康判定的要求，因此基于振动的结构动力无损检测方法受到了重视。它利用动力响应进行结构性能的识别，不受结构规模和构件隐蔽的限制，另外高效率模块化、数字化的动力响应量测技术为结构动力检测方法提供了有效的技术支持。动力检测方法的特点是应用条件限制少，效率高，既可以对复杂结构的整体也可以对结构的局部进行检测，包括人工激振法和环境随机激振法。而环境激励下的结构检测(监测)由于具有对构筑物无破坏、无环境污染、不打断结构的正常运营等优点而受到国内外的重视，成为国际社会研究的热点和前沿课题。模态识别是结构进行长期健康监测的基础，环境激励下模态参数的识别是仅基于输出的响应的系统识别。国内外研究者提出了多种方法，主要有功率谱峰值法(PP)、自然激励技术(NEXT)、频域分解法(FDD)和随机子空间法(SSI)等。环境随机激振下的结构动力响应测试数据，具有幅值小、随机性强和数据量巨大的特点，给空间网格结构的模态参数识别带来很大的难度和挑战。传统的功率谱峰值法在识别模态参数时具有速度快，直观性强，容易操作的特点，因此在建筑工程中具有较多的应用，但也存在一些问题。例如(1)功率谱峰值法适用的前提是要求频率必须稀疏， 而空间网格结构具有自振频率密集、空间耦合振动明显的特点，同时自然环境振动下的结构动力响应测试数据幅值小、随机性强，因此功率谱峰值法用于空间网格结构的模态参数识别时会遇到遗漏模态和重叠频率无法筛选的问题；(2)传统的方法在频率识别时单纯采用取极大值的数学方法，峰值选取过于主观，没有与结构的理论模型联系起来。
发明内容
本发明针对上述问题，提出了一种改进功率谱峰值法的空间网格结构频率识别方法，其在传统功率谱峰值法的基础上进行了改进，根据空间网格结构模态密集的特点提出了若干判断准则，是一种以结构的理论振型为基础的辅助正则化功率谱算法，使其适用于空间网格结构的频率识别，同时根据此改进算法编制了功率谱分析程序。改进后的算法具有速度快、避免遗漏以及自动识别密集频率的特点。本发明是通过以下技术方案实现的一种改进功率谱峰值法的空间网格结构频率识别方法，其特征在于先将结构测点的加速度响应信号进行预处理，接着进行测点的自功率谱计算，然后形成平均正则化功率谱曲线和辅助功率谱曲线，分别拾取各曲线的峰值点对应的频率值，经过幅角和理论频率的筛选及采用相位准则和振型准则剔除虚假频率，然后取二者的并集，最终得到结构各个阶次的实测频率。上述方法具体包括以下几个步骤(I)将实测的加速度响应信号进行预处理，计算各测点的自功率谱曲线，以所有测点的自功率谱曲线为目标进 行计算，形成全部测点的平均正则化功率谱曲线，然后选择各个阶次目标模态中均振动较大的测点形成辅助功率谱曲线；(2)提取正则化功率谱曲线和辅助功率谱曲线的峰值点分别对应的频率值α与β，首先进行幅角判断，然后进行理论频率筛选，得到筛选后的频率值α ’与β ’ ；(3)若筛选后的频率值α ’与β ’结果不存在重叠频率，则此结果为最终的结构实测频率，若存在重叠频率则采用相位准则进行进一步筛选得到频率值α ’’与β ’’ ；(4)若进一步筛选后的频率值α ’’与β ’’不存在重叠频率，则此结果为最终的结构实测频率，若仍存在重叠频率则采用振型准则进行筛选，从而得到基于正则化功率谱曲线、辅助功率谱曲线最终的结构实测频率值α ’’’与β ’’’ ；(5)取曲线正则化功率谱曲线、辅助功率谱曲线最终的结构实测频率值α ’’’与β ’’’的并集，得到最终的结构实测频率。所述步骤⑴中加速度响应信号的预处理包括平滑、滤波。所述步骤(2)中幅角判断方法为正则化功率谱曲线和辅助功率谱曲线峰值点对应的传递函数幅角值(即选取振动较大的测点为参考点，计算它与其余测点的互功率谱，可得出互功率谱曲线(与自谱曲线相似，同样会在频率处出现峰值)和幅角曲线，此时可得到每个峰值下对应的幅角值，也称相位角)在0° ±30°或者180° ±30°时为结构频率值，否则为虚假值。所述步骤(2)中理论频率筛选方法为将权利要求4中幅角判断中筛选后的结构频率值与理论频率值进行比较，上述结构频率值在理论频率值的5%范围内时为结构频率值，否则为虚假值。所述步骤(3)中的相位准则筛选方法为选择与重叠频率对应的理论振型的两个测点，计算二者的相位角，并分析相位关系，与理论振型中两测点相位关系不一致的为虚假频率。所述步骤(4)中的振型准则筛选方法为根据重叠频率对应的理论振型，选择振型分量较大处对应的测点，计算该测点的自谱，所得自谱中该振型具有明显峰值的为结构频率值，否则为虚假值。有益效果本发明的优点及效果是1.克服了传统功率谱峰值法中频率取值过于主观的缺点，本发明的方法可快速、高效的得出结构目标模态中每阶频率的值。2.该发明在获得所有测点平均正则化功率谱的基础上，结合结构的理论振型特点计算辅助正则化功率谱，有效的解决了密集模态的遗漏问题。3.该发明针对空间网格结构模态密集的特点，结合结构的理论模态采用相位准则、振型准则有效的剔除了虚假频率，解决了模态重叠的问题。4.针对此改进算法的编制了对应的程序，可方便有效的应用到实际工程的频率测试中。
图1为本发明的流程框图；图2为某测点的 加速度时程曲线；图3为经平滑、滤波预处理后的曲线；图4为实施例1各测点的平均正则化功率谱曲线和部分响应较大测点的辅助功率谱曲线；图5为实施例2各测点的平均正则化功率谱曲线和部分响应较大测点的辅助功率谱曲线。
具体实施例方式本发明的具体实施方式
如下实施例1某正交正放四角锥网架结构，经白噪声激励后提取了测点的加速度时程曲线，为了更好的模拟环境激励，同时将10%的噪声加入结构响应，以下为结构前6阶频率识别步骤(I)图2为某测点的加速度时程曲线，图3为经平滑、滤波预处理后的曲线。(2)经计算，形成测点的平均正则化功率谱曲线和部分响应较大测点的辅助功率谱曲线，如图4所示。经幅角判断和理论频率筛选后α’与β’的结果分别如表1、2所示。β’的结果避免了第6阶模态的遗漏。表I筛选后的α ’结果(Hz)
权利要求
1.一种改进功率谱峰值法的空间网格结构频率识别方法，其特征在于先将结构测点加速度响应信号进行预处理，接着进行测点的自功率谱计算，然后形成平均正则化功率谱曲线和辅助功率谱曲线，分别拾取各曲线的峰值点对应的频率值，经过幅角和理论频率的筛选及采用相位准则和振型准则剔除虚假频率，然后取二者的并集，最终得到结构各个阶次的实测频率。
2.如权利要求1所述的改进功率谱峰值法的空间网格结构频率识别方法，其特征在于所述方法具体包括以下几个步骤(1)将实测的加速度响应信号进行预处理，计算各测点的自功率谱曲线，以所有测点的自功率谱曲线为目标进行计算，形成全部测点的平均正则化功率谱曲线，然后选择各个阶次目标模态中均振动较大的测点形成辅助功率谱曲线；(2)提取正则化功率谱曲线和辅助功率谱曲线的峰值点分别对应的频率值α与β，首先进行幅角判断，然后进行理论频率筛选，得到筛选后的频率值α ’与β ’ ；(3)若筛选后的频率值α’与β ’结果不存在重叠频率，则此结果为最终的结构实测频率，若存在重叠频率则采用相位准则进行进一步筛选得到频率值α ’’与β ’’ ；(4)若进一步筛选后的频率值α’’与β ’’不存在重叠频率，则此结果为最终的结构实测频率，若仍存在重叠频率则采用振型准则进行筛选，从而得到基于正则化功率谱曲线、辅助功率谱曲线最终的结构实测频率值α ’’’与β ’’’ ；(5)取曲线正则化功率谱曲线、辅助功率谱曲线最终的结构实测频率值α’’’与β ’’’的并集，得到最终的结构实测频率。
3.如权利要求1或2所述的改进功率谱峰值法的空间网格结构频率识别方法，其特征在于所述加速度响应信号的预处理包括平滑、滤波。
4.如权利要求3所述的改进功率谱峰值法的空间网格结构频率识别方法，其特征在于所述步骤(2)中幅角判断方法为正则化功率谱曲线和辅助功率谱曲线峰值点对应的传递函数幅角值在0° ±30°或者180° ±30°时为结构频率值，否则为虚假值。
5.如权利要求4所述的改进功率谱峰值法的空间网格结构频率识别方法，其特征在于所述步骤(2)中理论频率筛选方法为将权利要求4中幅角判断中筛选后的结构频率值与理论频率值进行比较，上述结构频率值在理论频率值的5%范围内时为结构频率值，否则为虚假值。
6.如权利要求5所述的改进功率谱峰值法的空间网格结构频率识别方法，其特征在于所述步骤(3)中的相位准则筛选方法为选择与重叠频率对应的理论振型的两个测点，计算二者的相位角，并分析相位关系，与理论振型中两测点相位关系不一致的为虚假频率。
7.如权利要求6所述的改进功率谱峰值法的空间网格结构频率识别方法，其特征在于所述步骤(4)中的振型准则筛选方法为根据重叠频率对应的理论振型，选择振型分量较大处对应的测点，计算该测点的自谱，所得自谱中该振型具有明显峰值的为结构频率值，否则为虚假值。
全文摘要
本发明涉及一种改进功率谱峰值法的空间网格结构频率识别方法。先将加速度响应信号进行预处理(平滑、滤波)，接着进行各个测点的自功率谱计算，然后形成平均正则化功率谱曲线和辅助功率谱曲线，分别拾取两曲线的峰值点对应的频率值，然后经过幅角和理论频率的筛选以及采用相位准则和振型准则剔除其虚假频率，最后取两条曲线对应频率值的并集，最终得到结构各个阶次的实测频率。同时根据此改进算法编制了功率谱分析程序。适用于大型桥梁、输电塔、高层建筑等结构形式的频率识别，尤其适用于具有密集模态的大跨空间网格结构的频率识别，改进后的算法具有速度快、避免频率遗漏及自动识别密集频率的特点。
文档编号G01M7/02GK103033329SQ20121056824
公开日2013年4月10日 申请日期2012年12月24日 优先权日2012年12月24日
发明者吴金志, 刘才玮, 张毅刚 申请人:北京工业大学