小波包分析在结构损伤识别中应用.docVIP

小波包分析在结构损伤识别中应用 　　摘要:提出了基于小波包节点能量相对变化量的结构损伤识别方法,采用节点能量平方比作为损伤定位依据,从而能对结构损伤进行无模型定位,最后通过数值算例验证了该方法的可行性。 　　关键词:小波包分析 损伤识别 小波包节点能量 节点能量平方比 　　 　　0 引言 　　小波分析是近几年来发展起来的新的数学分支,非常适合对非稳态信号进行处理,被誉为数学显微镜,在土木工程领域有着广阔的应用前景[1]。小波包分析比小波分析具有更高的频率分辨率,在结构损伤识别中已经得到广泛的应用。本文采用小波包分析技术,对简支梁的实时损伤报警进行了数值研究,基于小波包节点能量给出了小波包节点能量平方比结构实时损伤预警判据,通过数值算例验证了该方法的可行性。 　　1 小波包分析理论 　　小波包是由一系列线性组合的小波函数ψi(t)组成,即 　　式中,i、j、k分别为频率因子?p尺度因子和平移因子。波包与其相应的小波函数具有相同的性质,如正交性?p时频定位性等。 　　小波包的定义如下: 　　这样对于任何一个信号f(t)∈L2(R)都可以进行如下分解: 　　其中f ij(t)为第j层小波包分解的第i个频带的信号。小波包分析树形结构如下图所示: 　　可以看出,小波包分解每次得到的是两个序列,而这两个序列都要进行进一步的分解,即同时对低频部分和高频部分进行分解。 　　2 基于小波包分析的结构损伤识别方法 　　基于小波包的结构损伤识别主要是基于小波包节点能量变化来识别结构损伤,它是利用小波包对发生损伤结构的动力反应进行分解,从而得到信号在各个尺度上的能量分布状况,并以此作为判别结构损伤的依据[2~5]。本文针对简支梁结构的能量分布特点提出了一种针对简支结构的小波包节点能量平方比损伤识别方法,其主要分析过程如下: 　　2.1 对测得的第n个结构响应信号进行j 层小波包分解,则末层形成等带宽的2j个频带,提取从低频到高频各频率成分的信号分解系数Wij(fn(t)),i=0,1,2,…(2j-1)。 　　2.2 对小波包分解系数重构,提取各频带范围的信号。以f ijn(t)表示Wij(fn(x))的重构信号,以此类推。则总信号f (t)可表示为 　　2.3 求各频带信号的能量。令第i个频带的能量为E ij(n),Anik为重构信号f ijn(t)各离散点的幅值,m为离散点的个数,其中,k=1,2,…,m,则 　　2.4 构造能量向量E ij(n)可对向量进行归一化处理。令 　　其中E ij(n)为第n个结构测点信号进行j尺度小波包分解时第i个频段内的信号能量,即为求出在整个结构中的能量沿结构的分布情况,m为测点总数; 　　2.5 构造WPCE的平方比为基础的平方比向量α ij(n):α ij(n)= 　　;式中,点n处的结构无损和结构损伤时的WPCE响应定义为 , 。 　　3 算例 　　3.1 模型概述 选用如图2所示的18单元的简支梁桥模型,模型的基本参数如下:单元长度L=0.1m,截面面积A=2.5×10-3m2,截面惯性矩I=0.76×10-8m4,材料密度ρ=7500kg/m3,弹性模量E=200Gpa,结构阻尼比ε=0.03。 　　在跨中10#节点上作用f(t)=60sin10tN的竖直力,运用编制的损伤模拟程序STD计算在荷载f(t)作用下的结构动力时程反应。计算时间取为30s,在15s时模拟结构突然出现一处或者多处损伤,即相应单元的刚度降低。通过程序分析得到结构的位移时程或者加速度时程反应,且用率模拟传感器采集数据,再对采集的数据进行小波包变换,分析不同频段上的结构损伤特征。损伤工况,如表1所示 　　3.2 单一损伤识别 　　将各质点发生损伤稳定以后的10s时间的加速度时程反应采用Db6小波作为小波包基函数进行4层分解?p重构以后进行归一化处理,得到平方比向量。这里仅列出低频信号E04(n)的分析结果,其它频段内也有类似的结果。 　　图3中,平方比向量α04(n)的曲线图在7#?p8#节点处达到了极大,可以初步判定此时结构的这两个节点处的刚度发生损伤,通过结构刚度矩阵的形成特点可以判定此时7#单元的刚度矩阵k7发生了损伤,以此可以判定结构损伤的发生及损伤位置的确定。同时,从图中也可以看出来随着损伤程度的增加,所选取的平方比向量变化非常明显,能够较好区分各不同的损伤程度,可以作为构实时损伤预警判据。 　　3.3 多处损伤识别 当结构出现工况11～工况13等多处损伤时,也可按单一损伤时的损伤识别方法对结构的时程反应采用Db6小波作为小波包基函数对时频信号进行4层分解?p重构以后进行归一化处理,得到平方比向量。这里仅列出低频信号E04(n)的分析结果,其它频段内也有类似的结果。