一种基于未知脉冲激励响应的结构阻尼比识别方法与流程

1.本发明属于工程结构检测数据分析技术领域，涉及基于未知脉冲激励响应的结构阻尼比识别方法。


背景技术：

2.工程结构检测是定期维护结构，保障结构安全的重要手段，而模态参数反应结构的动力特性，可用于结构性能的评估，因此，利用结构检测数据来识别结构模态参数至关重要。
3.结构的模态参数包括频率、振型、阻尼比，阻尼可以使振动系统的能量随时间的变化而消散，在实际工程中，阻尼不仅有助于减小结构共振的振幅，避免结构因动力反应达到极限应力遭到破坏，而且有助于在收到动力作用时快速恢复到稳定的状态。因此，在工程检测中准确的识别出阻尼比是非常重要的。
4.工程中常用阻尼识别方法有多种，曾储惠于2003年在对数衰减法的基础上推倒出了使用能量估计阻尼比的方法，能有效降低噪声的干扰；ivan wang于2011年分析了传统半功率带宽求解公式误差与结构自身阻尼比大小的关系，提出了基于三阶修正的阻尼比识别公式，提高了大阻尼情况下的识别精度；t.kijewski于2001年利用基小波变换和fft变换计算出阻尼比；brincker,r.于2006年提出了efdd,将频域信号转换的时域用对数衰减法求得阻尼比。fdd方法是基于频域的一种模态参数识别方法，可以在未知激励的情况下识别出固有频率、振型等模态参数，并且具有很好的抗噪性能，但是由于截断误差的原因，fdd方法对于阻尼比的识别误差比较大。


技术实现要素：

5.本发明旨在提供一种在未知脉冲激励下结构的阻尼比识别方法，解决频域分解法(fdd)因为截断误差而导致阻尼比识别不准确的问题。
6.本发明的技术方案：
7.推倒一种阻尼比识别方法，首先获得未知脉冲激励下的结构的响应信号，对响应信号的自相关函数和互相关函数进行傅里叶变换获得结构的输出功率谱密度矩阵，通过频域分解法(fdd)对输出功率功率谱密度函数进行奇异值分解获得有相位误差的振型以及固有频率。使用得到的有误差的振型、固有频率以及阻尼修正参数x计算出修正项，对在固有频率处的输出功率谱进行修正，通过变化阻尼修正参数x的值，直到由修正后的输出功率谱奇异值分解得到的振型的相位得到修正，将这时的x值带入到奇异值的公式当中计算出阻尼比，再将计算出来的阻尼比带入到修正项当中对上述过程进行多次迭代，即可获得阻尼比。
8.一种基于未知脉冲激励响应的结构阻尼比识别方法，步骤如下：
9.(1)获得结构的位移响应y＝[y1(t),y2(t),y3(t)
…yn
(t)]，使用位移响应求得其
相关函数对相关函数进行傅里叶变换获得功率谱密度函数,进而获得功率谱密度函数矩阵其对角线上为各自由度的自功率谱密度函数，非对角线上不同自由度间的互谱密度函数；
[0010]
(2)对输出功率谱密度函数进行奇异值分解，保留最大奇异值s和其所对应的u阵的第一列，然后绘出最大奇异值关于频率ω的曲线，其最大奇异值曲线峰值点所对应的频率，即为结构的固有频率，该峰值点处最大奇异值对应的酉矩阵即为该固有频率下结构的有相位误差的振型，各阶阻尼比初始假定为0.001；
[0011]
(3)阻尼比修正参数x1～xn的取值范围为0～1，令x1～xn的初始值均为δ＝0.01；令初始求解阶次p＝1；
[0012]
(4)利用第(2)步获得的结构固有频率和带有相位误差的振型构建第p阶输出功率谱修正项：
[0013][0014][0015]
其中：分别为(2)中获得的有相位误差的第p、q阶振型；ω
p
、ωq即为(2)中获得的第p、q阶固有频率；ξ
p
为第p阶阻尼比；x
p
、xq为阻尼比修正参数，其取值范围为0～1。
[0016]
(5)使用第p阶输出功率谱修正项对在固有频率ω
p
的输出功率谱进行修正修正后的输出功率谱进行奇异值分解得到修正后的振型
[0017]
(6)令p＝p+1，重复(4)～(5)步，直至p＝n，获得各阶的修正后的振型
[0018]
(7)如果每个振型中各项满足在该振型内虚部与实部的比值前四位相等，则记录下x1～xn的值，如果不相等，则保持x1～x
n-1
不变，xn＝xn+δ，令p＝1，重复(4)～(6)步；
[0019]
(8)若xn变化到1过程中，振型中各项实虚比前四位若仍不相等，则令x
n-1
＝x
n-1
+δ，令p＝1，重复(4)～(7)步；
[0020]
(9)当x
n-1
变化到1过程中，若仍不相等，令x
n-2
＝x
n-2
+δ，令p＝1，重复(4)～(8)步；依此类推，直至x1变化到1，若以上过程中振型中各项实虚比前四位相等，则终止循环，记录x1～xn的值，如果不等，令重新计算(4)～(9)步，直至求得满足相等条件时的x1～xn值；
[0021]
(10)将(9)中记录下的x1～xn的值带入阻尼比求解公式中：其中：s
p
为第p阶最大奇异值；为第p阶修正振型；为第p阶阻尼比；依次获得各阶阻尼比
[0022]
(11)令带入到(4)当中，重复(4)-(10)步，直到两次迭代获得的阻尼比误差小于0.01停止，即为结构真实的阻尼比。
[0023]
本发明的有益效果：可以在仅有输出的情况下计算出结构各阶的阻尼比，并且与传统fdd方法求解阻尼比不同，不需要经过傅里叶逆变换转换由频域变换到时域，而是直接在频域通过迭代求解，消除了截断误差对阻尼比识别的影响。不仅有较好的抗噪性能，而且提高了阻尼比的识别精度。
具体实施方式
[0024]
以下结合技术方案，进一步阐明本发明的实施方式。
[0025]
以一个3自由度结构为例，设其个自由度的质量分别为1kg、1.5kg、2kg，其刚度矩阵和阻尼矩阵如下：
[0026][0027][0028]
激励形式为单位脉冲响应，响应信号为结构的每层位移。
[0029]
(1)获得结构的位移响应y＝[y1(t),y2(t),y3(t)
…yn
(t)]，使用位移响应求得其相关函数对相关函数进行傅里叶变换获得功率谱密度函数,进而获得功率谱密度函数矩阵其对角线上为各自由度的自功率谱密度函数，非对角线上不同自由度间的互谱密度函数。
[0030]
(2)对输出功率谱密度函数进行奇异值分解，保留最大奇异值s和其所对应的u阵的第一列，然后绘出最大奇异值关于频率ω的曲线，其最大奇异值曲线峰值点所对应的频率，即为结构的固有频率，该峰值点处最大奇异值对应的酉矩阵即为该固有频率下结构的有相位误差的振型，各阶阻尼比初始假定为0.001；
[0031]
由此即可获得结构的固有频率：
[0032]
ω＝[14.52，31.04,46.10]
[0033]
固有频率处的第一奇异值s：
[0034]
s＝[1.2630，0.0355,0.0028]
[0035]
有相位误差的振型
[0036]
[0037]
各阶段固有频率为：
[0038]
ξ＝[0.0010.001,0.001]
[0039]
(3)阻尼比修正参数x1～xn的取值范围为0～1，令x1～xn的初始值均为δ＝0.01；令初始求解阶次p＝1；
[0040]
(4)利用第(2)步获得的结构固有频率和带有相位误差的振型构建第p阶输出功率谱修正项：
[0041][0042][0043]
其中：分别为(2)中获得的有相位误差的第p、q阶振型；ω
p
、ωq即为(2)中获得的第p、q阶固有频率；ξ
p
为第p阶阻尼比；x
p
、xq为阻尼比修正参数，其取值范围为0～1。
[0044]
(5)使用第p阶输出功率谱修正项对在固有频率ω
p
的输出功率谱进行修正修正后的输出功率谱进行奇异值分解得到修正后的振型
[0045]
(6)令p＝p+1，重复(4)～(5)步，直至p＝n，获得各阶的修正后的振型
[0046]
(7)如果每个振型中各项满足在该振型内虚部与实部的比值前四位相等，则记录下x1～xn的值，如果不相等，则保持x1～x
n-1
不变，xn＝xn+δ，令p＝1，重复(4)～(6)步；
[0047]
(8)若xn变化到1过程中，振型中各项实虚比前四位若仍不相等，则令x
n-1
＝x
n-1
+δ，令p＝1，重复(4)～(7)步；
[0048]
(9)当x
n-1
变化到1过程中，若仍不相等，令x
n-2
＝x
n-2
+δ，令p＝1，重复(4)～(8)步；依此类推，直至x1变化到1，若以上过程中振型中各项实虚比前四位相等，则终止循环，记录x1～xn的值，如果不等，令重新计算(4)～(9)步，直至求得满足相等条件时的x1～xn值；
[0049]
得到修正后的振型为：
[0050][0051]
(10)将(9)中记录下的x1～xn的值带入阻尼比求解公式中：其中：s
p
为第p阶最大奇异值；为第p阶修正振型；为第p阶阻尼比；依次获得各阶阻尼比
[0052]
各阶阻尼比为：
[0053][0054]
(11)令带入到(4)当中，重复(4)-(10)步，直到两次迭代获得的阻尼比误差小
于0.01停止，即为结构真实的阻尼比。
[0055]