基于谱图的桅杆结构杆件损伤识别.pdfVIP

基于谱图的桅杆结构杆件损伤识别 彭琦 武汉理工大学土木工程与建筑学院，武汉（430070 ） E-mail ：sinbaren@163.com 摘 要：桅杆结构是一种高柔结构，具有非常强的非线性，运用频域法对其损伤进行识别非 常困难。本文运用谱图的时频分布一阶频率矩的方法，对有杆件损伤的桅杆结构进行了分析， 并找到损伤弦杆的位置。 关键词：桅杆杆件，损伤识别，时频分析，谱图 1. 引言 桅杆结构的破坏事故时有发生，如何对桅杆进行损伤诊断近年来引起了工程界和学术界 的重视。总结过去桅杆结构的破坏事故，发现其破坏部位主要在杆身与纤绳连接的拉耳、弦 杆、腹杆与弦杆连接点[1] 。对损伤杆件的识别是对桅杆结构进行损伤识别的重要部分。由于 桅杆结构的动力特性相当复杂，如：纤绳力和纤绳位移的非线性、纤绳与杆身的耦合作用、 阻尼的非线性等，因此，传统的频域方法无法对桅杆进行损伤识别。时频分析的方法是时间 和频率的联合能量分布，能够反映非平稳信号的在时域和频域上的特征，为非线性结构的损 伤诊断提供了一条新的途径。 [2] 2. 短时傅立叶变换 与谱图一阶频率矩 桅杆结构产生的时程响应信号是非平稳信号，而传统的傅立叶分析只适用于平稳信号或 确定性信号。为克服传统的傅立叶变换在时域无任何分辨率的缺点，Gabor 于 1946 年提出 了短时傅立叶变换（Short-term Fourier Transform）的时频分析方法,其基本思想是在傅立叶 变换的基础上为了实现时域的局部化，首先强调围绕时刻t 的信号而忽略远离时刻t 的信号， 将信号与一个中心在时刻t 的宽度有限的窗函数h(t ) 相乘： z (τ) z (τ)h(τ−t) (1) t 在足够窄的窗口宽度中，信号可近似为平稳信号，则可以对加窗后的信号进行傅立叶变 换，得到只反映围绕时刻t 的一段信号的频率分布： ∞ X (t ,ω) ∫−∞z (τ)h (τ=−t )e −j ωτdτ (2) 上式中，t 表示有限时窗函数窗口的位置，ω表示局部化的频率，所以，X (t,ω) 是窗 口位置和频率的二元函数。其能量密度谱为： ∞ P (t ,ω) | X (t ,ω) |2 | ∫−∞z (τ)h(τ−t )e−j ωτdτ |2 (3) [3] P (t ,ω) 在时间和频率上的分布即为谱图分布 。 对于时频分析来说，矩是信号的一个重要信息，一个信号时频分布的一阶频率矩为： ∞ ∞ T (ω) ∫−∞t =×P(t,ω)dt / ∫−∞ P(t,ω)dt (4) 一阶频率矩对应了信号的群延迟。与时域信号 z (t ) 的瞬时频率相对应，频域信号的 Z (ω) 的群延迟g (ω) 也是一个重要的瞬时参数，它表示频谱Z (ω) 中频率为ω的各个分量 的延迟。 - 1 - 3. 数值分析 3.1 建立数值模型 考虑纤绳刚度硬化和大变形影响，用有限元分析 软件ANSYS[4]建立一个高度为150m的桅杆模型，有两 层三方纤绳。如图2 ，在A点施加沿Y轴作用的激励， 进行瞬态分析，得到该桅杆模型的响应，图2为位移 扩大10倍的响应图。