摘要
为了降低监测系统成本和提高桥梁损伤识别精度,本文利用少量传感器获得高分辨率模态振型. 在考虑路面粗糙度的情况下,基于柔度变化率进行损伤识别. 在不同路面粗糙度水平下,对移动车辆荷载作用下桥梁上所测位移响应进行特征正交分解,再将得到的主成分矩阵中高频动力分量信息过滤,获得高分辨率的模态振型. 利用少量传感器方法得到的高分辨率模态振型计算柔度变化率,将其作为桥梁损伤定位的指标. 为了验证该方法的有效性,对一简支梁进行模拟和试验. 结果显示,该方法在考虑不同粗糙度水平情况下,利用少量传感器得到高分辨率模态振型的柔度变化率能准确定位损伤.
温度变化、环境腐蚀和材料性能退化等因素,导致桥梁结构性能降低甚至出现损伤,严重影响交通正常运营和人民生命财产安全. 因此,桥梁结构的安全监测和损伤识别尤为重要,已成为当前研究的热
由于结构动力参数变化要比静力参数包含更多结构状态信息,所以多数研究基于结构的动力参数开展损伤识别. 而基于模态参数的方法原理简单、计算方便,在损伤识别研究中得到广泛关注,其中以“动力指纹”类为例,模态振
结构在发生损伤后,会导致刚度降低,柔度增加. 柔度具有对损伤更加敏感的特点,Raghavendrachar
无论是曲率模态还是模态柔度等动力指纹方法,均需要布置大量传感器才能获得较高分辨率的指标效果,高昂的成本导致其在实际工程中难以被推广应用. 而现有方法只能提供稀疏的、低分辨率的模态振型,难以满足基于模态振型的损伤定位. 因此利用少量传感器获取高分辨率模态振型具有重要意义. 一些学者对高分辨模态振型做了研究. Yang
针对上述问题,在考虑不同路面粗糙度的前提下,本文提出在移动车辆荷载作用下,通过对桥梁上少量传感器测得的位移信号进行特征正交分解(Proper Orthogonal Decomposition,POD)的方法来提取高分辨率的模态振型,然后利用高分辨率柔度变化率作为损伤指标进行损伤定位.
1 基本理论
1.1 桥梁高分辨率模态振型识别的少量传感器方法
一个移动荷载匀速通过长度为l的均质简支梁,如

图1 移动荷载作用下的简支梁模型
Fig.1 Simple bridge model under a moving load
(1) |
式中:分别表示简支梁的竖向位移、密度、阻尼系数、弹性模量、惯性矩;分别表示竖向移动荷载和移动速度.
采用振型叠加法,振型叠加的公式表示为:
(2) |
式中:表示第阶振型坐标;表示第阶振型.
简支梁在移动荷载作用下的竖向位移响
(3) |
式中:为第阶结构固有频率;,为无量纲参数;为第阶阻尼比;为移动荷载的第阶激励频率,;为阻尼振动频率.
对竖向位移响应的解析解进行POD分解,由于大多数桥梁由钢、混凝土组成,阻尼系数较小且材料密度较大,
(4) |
式中:为位移信号响应矩阵;为位移响应矩阵的协方差矩阵的特征向量.
由
Nie
(5) |
式中:为桥梁的第个频率;为采样频率;K是由计算的滤波长度.
利用低通滤波器把高频动力分量信息过滤可得到第阶模态振型.
该方法通过少量传感器即可获得高分辨率模态振型,如布置三个传感器即可获得前3阶高分辨率模态振型.
1.2 柔度变化率
Pandey
(6) |
式中:为第阶固有频率;为第阶振型.
Ko
(7) |
式中:与分别为结构损伤前后柔度矩阵,为取矩阵对角值.
2 数值模拟
2.1 仿真模型及损伤工况
有限元模型中,使用ANSYS中plane42单元建立简支梁模型. 一个质量为1 000 kg的移动质量块从简支梁的一端移动到另外一端,速度设置为1 m/s. 采用矩形截面钢材,矩形截面尺寸为0.1 m×0.2 m,密度为7 850 kg/

图2 简支梁模型
Fig.2 Simple beam model
路面粗糙度是影响系统分析的关键因素,较差的路面会显著增加车辆荷载和桥梁响应的振动幅度. 基于国际标
(8) |
(9) |
式中:为每单位长度的空间频率;为0.1 cycle/m;为位移功能密度函数值;为0~之间的随机相位角.
在梁的0.4和0.7位置处设置两个单损,通过沿梁高方向从桥底到中性面依次去除相应的单元来设置不同程度的损伤,截面损失率分别为10%、20%和30%,以及在梁的0.4l和0.7l位置处设立两个损伤,工况设置如
工况 | 损伤位置 | 损伤程度/% |
---|---|---|
1 | 0 | 0 |
2 | 0.4l | 10、20、30 |
3 | 0.7l | 10、20、30 |
4 | 0.4l、0.7l | 20 |
5 | 0.4l、0.7l | 40 |

图3 A粗糙度水平下工况1的位移响应
Fig.3 Displacement response at Scenario 1 of roughness A
2.2 数值模拟结果
在路面粗糙度的影响下,通过对桥梁上布置的7个传感器获取的位移响应进行POD处理,再通过低通滤波器将动力分量信息过滤,即可获得简支梁的前7阶模态振型.
特征值表示其对应分量在振动响应中所占据的能量,所以利用主成分分析中第n个分量的贡献率来估计能量意义上对应的贡献率(Contribution Ratio,CR
(10) |
式中:对应协方差矩阵的第个分量的特征值.
利用
阶数粗糙度 | 1st | 2nd | 3rd | 4th | 5th | 6th | 7th |
---|---|---|---|---|---|---|---|
A | 99.593 03 | 0.389 62 | 0.015 47 | 0.001 53 | 0.000 25 | 0.000 06 | 0.000 01 |
B | 99.576 88 | 0.405 87 | 0.015 34 | 0.001 56 | 0.000 25 | 0.000 06 | 0.000 01 |
模态阶数 | 1st | 2nd | 3rd |
---|---|---|---|
频率/Hz | 1.14 | 4.57 | 10.26 |
传统方法利用3个传感器只能获得3个振型数据,基于3个数据对工况3的损伤因子柔度变化率进行计算,如

图4 模态分析损伤识别结果
Fig.4 Damage identification results of modal analysis
以工况1为例,在车辆速度为1 m/s时,采样点个数为4 000,对桥梁位移信号进行POD处理,获得模态振型个数也为4 000,因此得到的桥梁模态振型具有较高的空间分辨率. 在A、B两个粗糙度水平下,3个传感器获得的前3阶高分辨率模态振型归一化后如

(a) A粗糙度前3阶振型

(b) B粗糙度前3阶振型
图5 A、B粗糙度识别的高分辨率模态振型
Fig.5 The identified high-resolution modal shapes of roughness A and B
利用A粗糙度水平的高分辨率模态振型计算柔度变化率,图中红色三角形表示损伤位置. 损伤识别结果如

(a) A粗糙度水平0.4l处10%、20%和30%损伤

(b) A粗糙度水平0.7l处10%、20%和30%损伤

(c) A粗糙度水平0.4l和0.7l处20%和40%损伤
图6 A粗糙度水平下损伤识别结果
Fig.6 Damage identification results of roughness A
增大路面粗糙度水平到B后,分别计算各个工况的柔度变化率损伤指标值. 如

(a) B粗糙度水平0.4l处10%、20%和30%损伤

(b) B粗糙度水平0.7l处10%、20%和30%损伤

(c) B粗糙度水平0.4l和0.7l处20%和40%损伤
图7 B粗糙度水平下损伤识别结果
Fig.7 Damage identification results of roughness B
以上结果均表明,在分别考虑了两个不同路面粗糙度水平影响的情况下,基于高分辨率模态振型计算得到的柔度变化率损伤因子仍然能够对损伤进行精确定位,而且随着损伤程度的增大,各个工况的柔度变化率也在增大,可以定性比较各个工况损伤程度的大小. 同时也可以看到,由于边界效应的影响,各个工况在梁两端的曲线会有比较大的波动.
2.3 噪声的鲁棒性
在土木工程结构健康监测中,不可避免地受到环境中各种噪声的影响. 为了验证该方法的噪声鲁棒性,把数值模拟中得到位移响应添加70 dB和 90 dB的高斯白噪音来模拟环境因素的影响. 通过信噪比衡量噪声的水平,信噪比公式如下:
(11) |
式中:表示原始信号;表示噪声;是信号的均方根;是信号的均方根;信噪比的单位为分贝.
以A粗糙度水平单损和双损工况为例,损伤识别结果分别如

(a) 90 dB噪声下A粗糙度水平0.4l处10%、20%和30%损伤

(b) 70 dB噪声下A粗糙度水平0.4l处10%、20%和30%损伤
图8 噪声影响下A粗糙度水平下单损识别结果
Fig.8 Single damage identification results of roughness A with noise

(a) 90 dB噪声下A粗糙度水平双损

(b) 70 dB噪声下A粗糙度水平双损
图9 噪声影响下A粗糙度水平下双损识别结果
Fig.9 Multiple damage identification results of roughness A with noise
3 试验研究
为了进一步验证本文提出的损伤识别方法的可行性,利用简支梁试验模型对其进行验证. 简支梁为一空心矩形钢梁,钢梁长度L=6 m,弹性模量E= 210 GPa,截面参数为200 mm×100 mm×3 mm(长×宽×厚),材料密度为7 850 kg/

图10 试验装置简图
Fig.10 Simple model of experimental device
由于试验梁表面并不是完全理想化的光滑,会存在加工以及环境因素造成梁表面的不平度,所以以此来等同于粗糙度,研究其对该方法的影响. 试验装置图如

图11 试验装置图
Fig.11 Experimental device
在梁0.68l位置处设立一处损伤,在梁0.38l、0.68l位置处设置两处损伤. 利用切割机对简支梁底部进行切割来制造损伤,在0.68l处设置一处损伤,损伤深度为7 mm,在0.38l处设置一处损伤,损伤深度分别为4 mm和7 mm,试验损伤工况如
工况 | 损伤位置 | 损伤深度/mm | 小车速度/(m· |
---|---|---|---|
1 | 0 | 0 | 0.5、0.75 |
2 | 0.68l | 7 | 0.5、0.75 |
3 | 0.38l、0.68l | 4、7 | 0.5、0.75 |
4 | 0.38l、0.68l | 7、7 | 0.5、0.75 |
以工况1中0.75 m/s小车速度为例,位移传感器采集到的位移响应如

(a) 位移响应

(b) 傅里叶频谱
图12 0.75 m/s位移响应及傅里叶频谱
Fig.12 Displacement responses and Fourier
spectrums of 0.75 m/s
以工况1中0.75 m/s小车速度为例,5个传感器采集位移信号并进行POD分析,可以得到前5阶高分辨率模态振型,如

图13 0.75 m/s工况1前5阶振型
Fig.13 First 5 mode shapes of scenarios 1 with 0.75 m/s
计算各个工况下的损伤因子柔度变化率,红色三角形标记点表示损伤位置,在小车速度为0.5 m/s时,从

(a) 单损工况

(b) 双损工况
图14 0.5 m/s的损伤识别结果
Fig.14 Damage identification results of 0.5 m/s
为了检验小车速度对该方法的影响,将小车速度增加到0.75 m/s时,损伤识别结果如

(a) 单损工况

(b) 双损工况
图15 0.75 m/s的损伤识别结果
Fig.15 Damage identification results of 0.75 m/s
为了进一步检验传感器数量减少时该方法的有效性,将上述试验中5个传感器数量减少到3个(分别取

(a) 单损工况

(b) 双损工况
图16 0.5 m/s时3个传感器的损伤识别结果
Fig.16 Damage identification results of 0.5 m/s with three sensors

(a) 单损工况

(b) 双损工况
图17 0.75 m/s时3个传感器的损伤识别结果
Fig.17 Damage identification results of 0.75 m/s with three sensors
4 结 论
本文在考虑了桥面粗糙度影响的情况下,获得移动荷载作用下桥梁少量传感器的位移响应,通过对位移响应进行特征正交分解(POD),得到高分辨率模态振型,并以此计算损伤因子柔度变化率,对桥梁进行损伤定位. 本文主要结论如下:
1)基于少量传感器的柔度变化率损伤指标能够精准识别损伤位置,而且在减少传感器的数量时,仍然具有很好的损伤识别效果.
2)该方法克服了传统模态方法只能获得稀疏振型的缺点,并且考虑了不同路面粗糙度水平和噪声对该方法的影响,更加适合实际工程.
3)利用少量传感器就能实现损伤识别,大大降低了成本.
参考文献
马宏伟, 聂振华.桥梁安全监测最新研究进展与思考[J].力学与实践,2015,37(2):161-170. [百度学术]
MA H W, NIE Z H.Recent advances and review of bridge safety monitoring[J].Mechanics in Engineering,2015,37(2):161-170.(in Chinese) [百度学术]
SHIH H W,THAMBIRATNAM D P,CHAN T H T.Damage detection in slab-on-girder bridges using vibration characteristics[J].Structural Control and Health Monitoring,2013,20(10):1271-1290. [百度学术]
WICKRAMASINGHE W R,THAMBIRATNAM D P,CHAN T H T,et al.Vibration characteristics and damage detection in a suspension bridge[J].Journal of Sound and Vibration,2016,375:254-274. [百度学术]
PANDEY A K,BISWAS M,SAMMAN M M.Damage detection from changes in curvature mode shapes[J].Journal of Sound and Vibration,1991,145(2):321-332. [百度学术]
PANDEY A K,BISWAS M.Experimental verification of flexibility difference method for locating damage in structures[J].Journal of Sound Vibration,1995,184(2):311-328. [百度学术]
HUNT D L.Applicrion of an enhanced coordinate modal assurance criteria[C]//10th International Modal Analysis Conference,San Diego,Califomia. CA: Society for Experiment Mechanics, Inc, 1992. [百度学术]
聂振华,程良彦,马宏伟.基于结构动力特性的损伤检测可视化方法[J].振动与冲击,2011,30(12):7-13. [百度学术]
NIE Z H,CHENG L Y,MA H W.Visualization method for structural damage detection based on its dynamic characteristics[J].Journal of Vibration and Shock,2011,30(12):7-13.(in Chinese) [百度学术]
RAGHAVENDRACHAR M,AKTAN A E.Flexibility by multireference impact testing for bridge diagnostics[J].Journal of Structural Engineering,1992,118(8):2186-2203. [百度学术]
YANG Y B,LI Y C, CHANG K C.Constructing the mode shapes of a bridge from a passing vehicle: a theoretical study[J]. Smart Structures and Systems,2014,13(5):797-819. [百度学术]
ESHKEVARI S S,MATARAZZO T J,PAKZAD S N.Bridge modal identification using acceleration measurements within moving vehicles[J].Mechanical Systems and Signal Processing,2020,141:106733. [百度学术]
NIE Z H,SHEN Z F,LI J,et al.Output-only complete mode shape identification of bridges using a limited number of sensors[J].Mechanical Systems and Signal Processing,2022,178:109246. [百度学术]
NIE Z H,NGO T,MA H W.Reconstructed phase space-based damage detection using a single sensor for beam-like structure subjected to a moving mass[J].Shock and Vibration,2017,2017:1-20. [百度学术]
PANDEY A K,BISWAS M.Damage detection in structures using changes in flexibility[J].Journal of Sound and Vibration,1994,169(1):3-17. [百度学术]
KO J M,SUN Z G,Ni Y Q.Multi-stage identification scheme for detecting damage in cable-stayed Kap Shui Mun Bridge[J].Engineering Structures,2002,24(7):857-868. [百度学术]
Mechanical vibration-Road surface profiles-Reporting of measured data:ISO 8608:2016[S]. London: BSI Standards Limited, 2016: 1-48. [百度学术]
王红岩,王钦龙,芮强,等.车辆行驶路面的数字化建模方法研究[J].兵工学报,2016,37(7):1153-1160. [百度学术]
WANG H Y,WANG Q L,RUI Q,et al.Research on digitized modeling method of riding road of vehicle[J]. Acta Armamentarii,2016,37(7):1153-1160.(in Chinese) [百度学术]
LIANG Y C,LEE H P,LIM S P,et al.Proper orthogonal decomposition and its applications:part Ⅰ:theory[J].Journal of Sound and Vibration, 2002, 252(3): 527-544. [百度学术]
WU C G,LIANG Y C,LIN W Z,et al.A note on equivalence of proper orthogonal decomposition methods[J].Journal of Sound and Vibration,2003,265(5):1103-1110. [百度学术]