D01:10.13374j.isml00103x.2007.2.B9 第29卷第12期 北京科技大学学报 Vol.29 No.12 2007年12月 Journal of University of Science and Technology Beijing Dec.2007 小波分析在钢筋混凝土梁损伤识别中的应用 薛刚)蔡美峰) 1)北京科技大学土木与环境工程学院.北京1000832)内蒙古科技大学建筑与土木工程学院,包头014010 摘要为识别钢筋混凝土简支梁的损伤状况,对简支梁进行了逐级加载实验.每级加荷后卸荷.观测梁的裂缝.并测定梁的 动力反应.将钢筋混凝土简支梁作为无限自由度体系.并对该体系的动力方程进行小波变换.得到多尺度下的结构动力响应 表达式.信号经多尺度分解后,包含了信号中更多的结构损伤信息.基于此,用DASP信号处理系统对钢筋混凝土简支梁各损 伤阶段的动态信号进行二进制小波分解,通过分析各频段的波形.确定了梁的损伤. 关键词钢筋混凝土梁:损伤识别:小波分析:多尺度分析 分类号TU3923 传统的傅里叶分析是一种纯频率的分析方法, 1小波变换的多尺度辨析基本原理 它在频域的定位性是完全准确的(即频域分辨率最 高),而在时域无任何定位性(或分辨能力).短时傅 将钢筋混凝土梁分解成n个单元的无限自由 里叶变换以固定大小的窗函数对时域特性有一定的 度体系.每个单元的质量为m,刚度为k,阻尼为 分析能力,而小波变换通过小波函数的伸缩和平移, c,如图1所示. 克服了短时傅里叶变换窗口大小不随频率变化、缺 (0 f0 f.( 乏离散正交基的缺点,是一种非常适合对非稳态信 号进行处理的数学工具一列,被誉为数学显微镜. 小波函数定义如下. 设Ψ(t)为一平方可积函数,即其傅里叶变换 F(四)满足条件: FC)12 daK oo R W (1) 图1无限自由度体系的动力系统模型 称平(1)为一个基本小波或小波母函数,并称式(1) Fig 1 Dynamic mode of multi-degree-of freedom systems 为小波函数的可允许性条件. 无限自由度体系的动力方程为: 损伤结构的动力特性与结构的损伤有着直接的 M(t)+G(t)+Ky(t)=F(t)(2) 关系,所以结构的动力响应信号是结构损伤识别的 原始资料4.近年来小波分析技术逐渐被应用于 其中M、C和K分别为整体质量、阻尼及刚度矩 结构损伤识别领域:韩建刚等人提出了用于损 阵,均为nXn对角矩阵;y(1)、y(t)、y(t)分别为 加速度向量、速度向量及位移向量;F()为荷载向 伤识别的小波包能量变化率指标,并根据此指标提 量,F(t)=(f1(1),f2(t,f3t,:fm(t). 出了用其进行损伤识别的实施步骤1o:Kim等人采 由式(2)得: 用Gabor小波识别简支梁的裂缝不仅能识别损伤 位置而且也能估计损伤程度山.本文对三根钢筋 y(t)+M Cy(t)+MKy(1)=MF(t). 用正交小波基中t)对上式两边进行小波变换1四, 混凝土简支梁在逐级加荷一卸荷后的动力特性进行 得到: 了测试,研究了在各级荷载作用下的小波原始信号 经过多尺度分解后的波形图,分析了基于小波变换 WT(a,b+品wT,(a,b+点wT,(a,b= 多尺度带通滤波器特性在结构损伤识别中的应用. T,(a,b创,i=l.2n (3) 收稿日期:200609-26修回日期:2007-0-05 基金项目:内蒙古自治区教有厅资助项目(No.NJ06075) 根据小波变换的多分辨率分析,式(3)中各式的 作者简介:薛刚(1968一),男,副教授.博士研究生:蔡美峰 每一项都可以分解为不同尺度上的细节部分和概貌 (1943一).男,教授,博士生导师
小波分析在钢筋混凝土梁损伤识别中的应用 薛 刚1, 2) 蔡美峰1) 1) 北京科技大学土木与环境工程学院, 北京 100083 2) 内蒙古科技大学建筑与土木工程学院, 包头 014010 摘 要 为识别钢筋混凝土简支梁的损伤状况, 对简支梁进行了逐级加载实验, 每级加荷后卸荷, 观测梁的裂缝, 并测定梁的 动力反应.将钢筋混凝土简支梁作为无限自由度体系, 并对该体系的动力方程进行小波变换, 得到多尺度下的结构动力响应 表达式.信号经多尺度分解后, 包含了信号中更多的结构损伤信息.基于此, 用 DASP 信号处理系统对钢筋混凝土简支梁各损 伤阶段的动态信号进行二进制小波分解, 通过分析各频段的波形, 确定了梁的损伤. 关键词 钢筋混凝土梁;损伤识别;小波分析;多尺度分析 分类号 TU 392.3 收稿日期:2006-09-26 修回日期:2007-01-05 基金项目:内蒙古自治区教育厅资助项目( No .NJ06075) 作者简介:薛 刚( 1968—) , 男, 副教授, 博士 研究生;蔡美峰 ( 1943—) , 男, 教授, 博士生导师 传统的傅里叶分析是一种纯频率的分析方法, 它在频域的定位性是完全准确的(即频域分辨率最 高), 而在时域无任何定位性( 或分辨能力) .短时傅 里叶变换以固定大小的窗函数对时域特性有一定的 分析能力, 而小波变换通过小波函数的伸缩和平移, 克服了短时傅里叶变换窗口大小不随频率变化 、缺 乏离散正交基的缺点, 是一种非常适合对非稳态信 号进行处理的数学工具[ 1-3] , 被誉为数学显微镜 . 小波函数定义如下. 设 Χ( t )为一平方可积函数, 即其傅里叶变换 F ( ω)满足条件: ∫R F( ω) 2 ω d ω<∞ ( 1) 称 Χ( t )为一个基本小波或小波母函数, 并称式( 1) 为小波函数的可允许性条件. 损伤结构的动力特性与结构的损伤有着直接的 关系, 所以结构的动力响应信号是结构损伤识别的 原始资料 [ 4] .近年来, 小波分析技术逐渐被应用于 结构损伤识别领域[ 1, 5-9] :韩建刚等人提出了用于损 伤识别的小波包能量变化率指标, 并根据此指标提 出了用其进行损伤识别的实施步骤[ 10] ;Kim 等人采 用Gabo r 小波识别简支梁的裂缝, 不仅能识别损伤 位置, 而且也能估计损伤程度[ 11] .本文对三根钢筋 混凝土简支梁在逐级加荷-卸荷后的动力特性进行 了测试, 研究了在各级荷载作用下的小波原始信号 经过多尺度分解后的波形图, 分析了基于小波变换 多尺度带通滤波器特性在结构损伤识别中的应用 . 1 小波变换的多尺度辨析基本原理 将钢筋混凝土梁分解成 n 个单元的无限自由 度体系.每个单元的质量为 mi , 刚度为 ki , 阻尼为 ci , 如图 1 所示. 图 1 无限自由度体系的动力系统模型 Fig.1 Dynamic mode of multi-degree-of-freedom systems 无限自由度体系的动力方程为: My ·· ( t) +Cy · ( t) +Ky( t ) =F( t) ( 2) 其中 M 、C 和 K 分别为整体质量、阻尼及刚度矩 阵, 均为 n ×n 对角矩阵;y ( t ) 、y · ( t) 、y ·· ( t) 分别为 加速度向量、速度向量及位移向量 ;F( t) 为荷载向 量, F( t ) =( f 1( t) , f 2( t), f 3( t), …, f n ( t) ) . 由式( 2)得: y ·· ( t) +M -1 Cy · ( t ) +M -1 Ky( t) =M -1 F( t) . 用正交小波基 ψ( t) 对上式两边进行小波变换[ 12] , 得到: W Ty ·· i ( a, b) + ci mi WT y · i ( a, b) + k i mi WTy i ( a, b) = 1 mi WT f i ( a, b) , i =1, 2, …, n ( 3) 根据小波变换的多分辨率分析, 式( 3)中各式的 每一项都可以分解为不同尺度上的细节部分和概貌 第 29 卷 第 12 期 2007 年 12 月 北 京 科 技 大 学 学 报 Journal of University of Science and Technology Beijing Vol.29 No.12 Dec.2007 DOI :10.13374/j .issn1001 -053x.2007.12.039
。1192 北京科技大学学报 第29卷 部分,即为: (dd,d,(d,d,d)及( (d+c+(d,+c)+点d.+= 可,,可,可以对不同损伤率(,52,)造 或+c,=l2yn 成的结构的刚度损伤进行识别. (4) 由式(5)可知:(1)结构不同位置的加速度信号 利用正交小波变换不同尺度间的正交性,得到 y(t)对一定的损伤有不同的敏感性:(2)同一位置 式(4)中各式在各尺度上的关系如下: 的加速度信号在不同尺度上的分解信号,对损伤的 (d,+,++c+(d,+d,+d,+ 敏感性不同:(3)各加速度信号在不同尺度上的能量 分布与激励荷载F(t)的频带能量分布有关:(4)不 c+鼎d++汁= 同位置上的加速度信号和不同尺度上的分解信号其 信息含量是不同的,结构不同位置的传感器信号决 m(dh,+di,+dj,-ei). 定了所能提供的系统最大信息,而信号的多尺度分 解能够尽可能地挖掘出原始信号中所包含的结构损 伤信息. (d,+d,++d,+)+ 2基于小波分析的钢筋混凝土梁损伤 (,+明++4,++(d,+ 识别 d+叶d+g)=d+近,++ 2.1实验概况 制作了三根钢筋混凝土简支梁,梁的几何参数 或+,i=l,2,;n (5) 为150mm×250mm×2200mm.配筋为:纵向受拉 式中J表示尺度.当结构刚度参数(k1,k2,;km) 钢筋为二根16mm的钢筋,受压钢筋为二根 由未损伤状态(k,k3,:k9)变化到损伤状态(k, 中12mm的钢筋,箍筋6mm,双肢箍,间距为 k经,…时,可以用一个损伤率向量(,2,, 100mm.钢筋级别为HPB235.混凝土强度等级为 n)来表示其中0≤≤1.可以得到: C25.采用集中力逐级加载的方法,每级荷载为 k 「1-51 5kN,在每一级荷载下实验步骤为:(1)每一级荷载 0 0 下持荷,用DH3816静态应变测试系统测定控制截 0 1-52 0 (6) 面受拉钢筋及受压边缘混凝土的应变,并观察裂缝 : 位置、高度和宽度;(2)卸荷至零,观察裂缝位置、高 k 0 0 度和宽度:(3)每级荷载卸除后,用北京东方振动研 应用多尺度分析,通过观察结构系统输出的加 究所的DASP系统进行动态信号采集.实验在内蒙 速度信号(少1(t以,少2(t),少n(t)》在不同小波空 古科技大学结构实验室完成图2为实验装置示 间和尺度空间上的分解即(ddd), 意图. 力锤激振位置 8 9 10只支座压梁 支座压梁 、钢轴 钢轴 加速度传感器位置 地锚螺栓 地锚螺栓 试件 钢轴 名厚钢板 厚钢板 IA 图2实验装置 Fig 2 Experimental setup 2.2小波变换结果及分析 从图3和图4比较可以看到:随着外荷载的逐 对上述采集信号进行二进制小波分解,可得到 渐增加,频带N5(96~192Hz、N6(192~384Hz、 分解的波形图.图3和图4为实验梁的各级荷载作 N7(384~769Hz)和N8(769~1538Hz)的幅值越来 用下,原采集信号的二进制小波分解波形图. 越小,包含的能量降低所以波形衰减越来越快:对
部分, 即为 : ( d 1 y ·· i +c 1 y ·· i ) + ci mi ( d 1 y · i +c 1 y · i ) + ki mi ( d 1 y i +c 1 y i ) = 1 mi ( d 1 f i +c 1 f i ), i =1, 2, …, n ( 4) 利用正交小波变换不同尺度间的正交性, 得到 式( 4)中各式在各尺度上的关系如下: ( d 1 y ·· i +d 2 y ·· i +d 3 y ·· i +c 3 y ·· i ) + ci mi ( d 1 y · i +d 2 y · i +d 3 y · i + c 3 y · i ) + ki mi ( d 1 y i +d 2 y i +d 3 y i +c 3 y i ) = 1 mi ( d 1 f i +d 2 f i +d 3 f i +c 3 f i ), …, ( d 1 y ·· i +d 2 y ·· i +…+d J y ·· i +c J y ·· i ) + ci mi ( d 1 y · i +d 2 y · i +…+d J y · i +c J y · i ) + ki mi ( d 1 y i + d 2 y i +…+d J y i +c J y i ) = 1 mi ( d 1 f i +d 2 f i +…+ d J f i +c J f i ), i =1, 2, …, n ( 5) 式中 J 表示尺度 .当结构刚度参数( k 1, k 2, …, kn ) 由未损伤状态( k 0 1, k 0 2, …, k 0 n )变化到损伤状态( k d 1, k d 2, …, k d n) 时, 可以用一个损伤率向量( ξ1, ξ2, …, ξn )来表示, 其中 0 ≤ξi ≤1 .可以得到 : k d 1 k d 2 k d n = 1 -ξ1 0 … 0 0 1 -ξ2 … 0 0 … 0 1 -ξn k 0 1 k 0 2 k 0 n ( 6) 应用多尺度分析, 通过观察结构系统输出的加 速度信号( y ·· 1( t), y ·· 2( t ), …, y · n ( t)) 在不同小波空 间和尺度空间上的分解, 即( d 1 y ·· 1 , d 2 y ·· 1 , …, d J y ·· 1 ), ( d 1 y ·· 2 , d 2 y ·· 2 , …, d J y ·· 2 ), …, ( d 1 y ·· n , d 2 y ·· n , …, d J y ·· n )及( c J y ·· 1 , c J y ·· 2 , …, c J y ·· n ), 可以对不同损伤率( ξ1, ξ2, …, ξn ) 造 成的结构的刚度损伤进行识别. 由式( 5)可知:( 1) 结构不同位置的加速度信号 y ·· ( t)对一定的损伤有不同的敏感性;( 2) 同一位置 的加速度信号在不同尺度上的分解信号, 对损伤的 敏感性不同 ;( 3)各加速度信号在不同尺度上的能量 分布与激励荷载 F( t )的频带能量分布有关 ;( 4)不 同位置上的加速度信号和不同尺度上的分解信号其 信息含量是不同的, 结构不同位置的传感器信号决 定了所能提供的系统最大信息, 而信号的多尺度分 解能够尽可能地挖掘出原始信号中所包含的结构损 伤信息 . 2 基于小波分析的钢筋混凝土梁损伤 识别 2.1 实验概况 制作了三根钢筋混凝土简支梁, 梁的几何参数 为 150 mm ×250 mm ×2 200 mm .配筋为 :纵向受拉 钢筋 为二 根 16 mm 的 钢筋, 受 压钢 筋 为二 根 12 mm 的 钢 筋, 箍 筋 6 mm, 双 肢 箍, 间距 为 100 mm .钢筋级别为 HPB235 .混凝土强度等级为 C25 .采用集中力逐级加载的方法, 每级荷载为 5 kN, 在每一级荷载下实验步骤为 :( 1) 每一级荷载 下持荷, 用 DH3816 静态应变测试系统测定控制截 面受拉钢筋及受压边缘混凝土的应变, 并观察裂缝 位置 、高度和宽度;( 2) 卸荷至零, 观察裂缝位置、高 度和宽度 ;( 3)每级荷载卸除后, 用北京东方振动研 究所的 DASP 系统进行动态信号采集 .实验在内蒙 古科技大学结构实验室完成, 图 2 为实验装置示 意图. 图 2 实验装置 Fig.2 Experimental setup 2.2 小波变换结果及分析 对上述采集信号进行二进制小波分解, 可得到 分解的波形图.图 3 和图 4 为实验梁的各级荷载作 用下, 原采集信号的二进制小波分解波形图. 从图 3 和图 4 比较可以看到 :随着外荷载的逐 渐增加, 频带 N5( 96 ~ 192 Hz) 、N6( 192 ~ 384 Hz) 、 N7( 384 ~ 769 Hz)和 N8( 769 ~ 1 538 Hz) 的幅值越来 越小, 包含的能量降低, 所以波形衰减越来越快;对 · 1192 · 北 京 科 技 大 学 学 报 第 29 卷
第12期 薛刚等:小波分析在钢筋混凝土梁损伤识别中的应用 。1193。 于N4(48~96Hz)及N3(24~48Hz)频带,幅值则越 显大于在第2激振点的幅值.在实验过程中观察 来越大,相应的能量增加,波形衰减较慢.钢筋混凝 到,实验梁在第6激振点产生了裂缝.可见在相同 土梁在不同损伤状态下,由二进制小波分解后的波 的条件下,经过小波分解的波形图幅值较大处即为 形图有较明显的区别,可以判断损伤的存在 裂缝位置. (a) 时域波形 可 28 i 时域波形 -28E 0.050.100.150.200.250.30 0.050.100.150.200.250.30 拟小波分频和重构波形 拟小波分频和重构波形 NI NI N2 N2 N3 3 N4 N5 AAA N6 N6 N7 N7 N8 N8 Cn Cn 0.05 0.100.150.200.250.30 y.ww 0 0.050.100.150.200.250.30 时间s 时间s (b) 时域波形 b 时域波形 20F -20 0 0.050.100.150.200.250.30 0 0.050.100.15 0.200.250.30 拟小波分频和重构波形 拟小波分频和重构波形 NI N2 N3 N4 N3 N4 N5 N6 N7 N7 N8 N8 Cn Cn 0 0.050.100.150.200.250.30 n等 0.050.100.150.200.250.30 时间s 时间s (c) (c) 时域波形 时域波形 gn 0 0.050.100.150.200.250.30 0 0.050.100.150.200.250.30 拟小波分频和重构波形 拟小波分频和重构波形 NI NI N2 N2 N3 N3 N4 N4 VAAAAAAA N5 N5 A N6 N6 N7 N7 N8 N8 Cn Cn 0.050.100.150200.250.30 W 0 0.05 0.10 0.150.20 0.25 0.30 时间s 时间s 图3测点2处小波多尺度分解波形图.()加载前:(b)外荷载 图4测点6处小波多尺度分解波形图.(加载前:()外荷载 30kN(c外荷载50kN 30kN(c外荷载50kN Fig.3 Wavelet multi-scale decomposition signal of Point 2:(a) Fig.4 Wavelet multi-scale decomposition signal of Point 6:(a) before loading:(b)loading for 30kN;(c)loading for 50kN before loading (b)loading for 30 kN;(c)loading for 50 kN 从实验梁的第2激振点和第6激振点可以看 到:在各阶荷载下,两个激振点处经过小波分解的波 3结论 形图是不同的,在第6激振点处的波形图的幅值明 (1)本文分析了结构发生损伤时其动力参数的
于N4( 48 ~ 96Hz)及 N3( 24 ~ 48 Hz)频带, 幅值则越 来越大, 相应的能量增加, 波形衰减较慢 .钢筋混凝 土梁在不同损伤状态下, 由二进制小波分解后的波 形图有较明显的区别, 可以判断损伤的存在. 图3 测点 2 处小波多尺度分解波形图.( a) 加载前;( b) 外荷载 30 kN;( c) 外荷载 50 kN Fig.3 Wavel et multi-scale decomposition signal of Point 2:( a) before loading;( b) loading for 30 kN ;( c) loading for 50 kN 从实验梁的第 2 激振点和第 6 激振点可以看 到:在各阶荷载下, 两个激振点处经过小波分解的波 形图是不同的, 在第 6 激振点处的波形图的幅值明 显大于在第 2 激振点的幅值 .在实验过程中观察 到, 实验梁在第 6 激振点产生了裂缝.可见, 在相同 的条件下, 经过小波分解的波形图幅值较大处即为 裂缝位置. 图 4 测点6 处小波多尺度分解波形图.( a) 加载前;( b) 外荷载 30 kN;( c) 外荷载 50 kN Fig.4 Wavelet multi-scale decomposition signal of Point 6:( a) before loading;( b) loading for 30 kN;( c) loading for 50 kN 3 结论 ( 1) 本文分析了结构发生损伤时其动力参数的 第 12 期 薛 刚等:小波分析在钢筋混凝土梁损伤识别中的应用 · 1193 ·
。1194 北京科技大学学报 第29卷 变化特征.应用二进制小波变换,把结构动力系统 社,200436 的数学方程在不同尺度进行了分解,得到不同尺度 [4 Ismail Z,Razak HA,Rahman A G A.Detemination of damage location in RC beams using mode shape derivatives.En Struct, 上的动力参数的描述.通过对比不同尺度下的频率 2006.28:1566 波形,证明了应用小波分析方法可以得到更多的结 【习郭健。基于小波分析的结构损伤识别方法研剂学位论刘.杭 构损伤信息. 州:浙江大学,20046 (2)对不同损伤状态下钢筋混凝土简支梁的动 【(张和生,薛松海,陈铉等.基于小波多分辨滤波特性的结构损 力信号进行小波分解分解后的波形图可较清晰的 伤识别.力学季刊,2004,25(1):152 【7孙增寿,韩建刚,任伟新,等.基于小波分析的结构损伤检测研 显示,随着裂缝出现固有频率降低,低频控制的带 究进展.地震工程与工程振动,2005,25(2):93 宽部分波形振幅越来越大,衰减越来越慢,而高频带 【习唐和生,薛松涛,陈铭等.基于小波变换在线结构损伤检测分 宽波形衰减较快. 析研究.固体力学学报,2004,25(1):47 (3)不同位置的加速度传感器在不同分解尺度 [9 Khadem S E.Rezace M.Development of vibration signature anah 上所反映的损伤信息量是不同的.在离损伤位置近 ysis using multiw avelet systems.J Sound Vb,2003.261:613 的地方,频带上的损伤指标更为明显. 10 Han JG.Sun Z S,Ren W X.Wavelet based damage identifica tion of beams//Proceedings of the Eighth Intemational Sympo 参考文献 sium on Structural Engineering for Young Experts.Beijing:Sci- ence Press,2004:356 [刂冉启文.小波变换与分数傅里叶变换理论及应用哈尔滨:哈尔 Il]Kim H,Melhem H.Damage detection of stnuctures by wavelet 滨工业大学出版社,2001:1 analy sis.Eng Struct,2004.26:347 【程正兴.小波分析算法与应用.西安:西安交通大学出版社, 【1☑孙延奎.小波分析及其应用.北京:机械工业出版社,2005: 2001:3 171 【习徐长发,李国宽.实用小波分析方法.武汉:华中科技大学出版 Application of wavelet analysis to damage detection of reinforced concrete beams XUE Gang 2),CAI Meifeng 1)Civil and Environmental Engineering School,University of Science and Techmology Beiing,Beijing 100083.Chim 2)Architecture and Civil Engineering School,Inner Mongolia University of Scienceand Technology.Baotou 014010.China ABSTRACT In order to detect the damage of simple-supported reinforced concrete (RC)beams,experiment set-up for RC beams dynamic test was designed.In the test,the beam was subjected to an increasing static load in the middle to introduce cracks,and the last load step corresponded to failure.After each load step and unload- ing,the beam crack width was investigated and an experimental dynamic monitoring was performed.The multi- ple degree-of-freedom dynamic equation of the testing beam was wavelet transfo med,and the structural dy nam- ic response expression on multi-scale was acquired.It is show n that multi-scale decomposition of signal comprises more structural damage information.By means of DASP signal processing system,binary wavelet transformation is applied to the dynamic signal of reinforced concrete beams on different damage states.Through analyzing the wave-figure of all frequency scales,the damage of reinforced co ncrete beams is determined. KEY WORDS reinforced concrete beams;damage identification;wavelet analy sis;multi-scale analysis
变化特征.应用二进制小波变换, 把结构动力系统 的数学方程在不同尺度进行了分解, 得到不同尺度 上的动力参数的描述 .通过对比不同尺度下的频率 波形, 证明了应用小波分析方法可以得到更多的结 构损伤信息. ( 2) 对不同损伤状态下钢筋混凝土简支梁的动 力信号进行小波分解, 分解后的波形图可较清晰的 显示, 随着裂缝出现, 固有频率降低, 低频控制的带 宽部分波形振幅越来越大, 衰减越来越慢, 而高频带 宽波形衰减较快 . ( 3) 不同位置的加速度传感器在不同分解尺度 上所反映的损伤信息量是不同的 .在离损伤位置近 的地方, 频带上的损伤指标更为明显. 参 考 文 献 [ 1] 冉启文.小波变换与分数傅里叶变换理论及应用.哈尔滨:哈尔 滨工业大学出版社, 2001:1 [ 2] 程正兴.小波分析算法与应用.西安:西安交通大学出版社, 2001:3 [ 3] 徐长发, 李国宽.实用小波分析方法.武汉:华中科技大学出版 社, 2004:36 [ 4] Ismail Z, Razak H A, Rahman A G A.Determination of damage location in RC beams using mode shape derivatives.Eng Struct, 2006, 28:1566 [ 5] 郭健.基于小波分析的结构损伤识别方法研究[ 学位论文] .杭 州:浙江大学, 2004:6 [ 6] 张和生, 薛松涛, 陈 , 等.基于小波多分辨滤波特性的结构损 伤识别.力学季刊, 2004, 25( 1) :152 [ 7] 孙增寿, 韩建刚, 任伟新, 等.基于小波分析的结构损伤检测研 究进展.地震工程与工程振动, 2005, 25( 2) :93 [ 8] 唐和生, 薛松涛, 陈 , 等.基于小波变换在线结构损伤检测分 析研究.固体力学学报, 2004, 25( 1) :47 [ 9] Khadem S E, Rezaee M .Development of vibration signature analysis using multiw avelet systems.J Sound Vib, 2003, 261:613 [ 10] Han J G, S un Z S, Ren W X.Wavelet based damage identification of beam s∥Proceedings of the Eighth In ternational Symposium on S tructural Engineering for Young Experts.Beijing :Science Press, 2004:356 [ 11] Kim H, Melhem H .Damage detection of structures by wavelet analysis.Eng Struct, 2004, 26:347 [ 12] 孙延奎.小波分析及其应用.北京:机械工业出版社, 2005: 171 Application of wavelet analysis to damage detection of reinforced concrete beams X UE Gang 1, 2) , CAI Meifeng 1) 1) Civil and Environment al Engineering S chool, Uni versit y of Science and Tech nology Beijing, Beijing 100083, C hina 2) Architecture and Civil Engineering S chool, Inner Mongolia Universit y of Science and Technology, Baot ou 014010, China ABSTRACT In o rder to detect the damage of simple-supported reinfo rced concrete ( RC) beams, experiment set-up fo r RC beams dynamic test w as desig ned .In the test, the beam w as subjected to an increasing static load in the middle to introduce cracks, and the last load step corresponded to failure.After each load step and unloading , the beam crack w idth w as investigated and an experimental dynamic monitoring was perfo rmed .The multiple deg ree-of-freedom dynamic equation of the testing beam was wavelet transfo rmed, and the structural dy namic response expression on multi-scale w as acquired .It is show n that multi-scale decomposition of signal comprises more structural damage information .By means of DASP signal processing system, binary w avelet transformation is applied to the dynamic sig nal of reinforced concrete beams on different damage states .Through analyzing the w ave-figure of all frequency scales, the damage of reinforced co ncrete beams is determined . KEY WORDS reinforced concrete beams ;damage identification ;wavelet analy sis ;multi-scale analysis · 1194 · 北 京 科 技 大 学 学 报 第 29 卷