《1. 引言》

1. 引言

阻尼比与大跨度桥梁中涡激振动的发生有很强的相关性[1,2]。为了估算大跨度桥梁的阻尼比,许多研究人员已经使用了唯输出运行模态分析(OMA),而并非使用需要暂时关闭使用中的桥梁的重型激发器。唯输出OMA的重要假设之一是结构系统应处于平稳的环境振动之下。然而,非平稳负荷如地震、极端风力和交通是民用基础设施激励信号的主要来源,并且这些负荷会导致非平稳响应[3–5]。因此,对于经典OMA的基本假设,违反激励条件可能是阻尼比估测不准确的原因之一。

交通荷载是桥梁的主要荷载源[6,7],并且可以表示为平稳的随机过程,这是因为道路的粗糙度被模拟为零均值平稳高斯随机过程[8]。但是,当交通量较低时,传感器位置处的环境振动信号在车辆接近并远离时被视为信号包络。在特定位置观察到的交通导致振动(TIV)可以被看作是一个非平稳过程,其表示为静止白噪声和包络函数[9]的乘积。因此,受到交通负荷影响的桥梁荷载加速度会在特定位置含局部化特征[8]

交通荷载通常会激发与车辆频率相对应的结构模态。《安大略公路桥梁设计规范》[10]建议当结构的主要固有频率为2~5Hz时,增加放大系数,美国国家高速公路和交通运输协会(AASHTO)规定一般卡车的频率应该是2.5Hz[11]。Bartos[11]也提出,当结构频率为1.5~5Hz时,动态放大效应可能会超过限制。车桥相互作用也会受驱动频率的影响,这取决于车辆穿越单个纵梁的持续时间[12,13]。Brewick[14]发现这些模态信息的失真增加了基于OMA方案的阻尼估计的不确定性。

本研究的重点是交通荷载引起的非平稳性及其对基于OMA的阻尼估计的影响。交通荷载的影响通过悬索桥的现场测量数据进行检查。通过引入调幅(AM)函数,我们提出了一种信号平稳化算法,并应用于悬索桥的运行监测数据来进行阻尼估计。

《2.悬索桥的TIV研究》

2.悬索桥的TIV研究

《2.1.索罗克大桥》

2.1.索罗克大桥

索罗克大桥(图1)是韩国连接索罗克岛和大陆的自锚式单索悬索桥。桥梁全长470m,由250m主跨和两个110m的对称边跨组成。其钢箱梁总宽度为15.7m,提供两条车道。该桥于2009年3月通车。从那时起,研究者详细检测,并进行了一系列的动态测试,同时出于维护目的进行了模型更新。表1总结了桥梁的模态频率,由OMA使用从环境振动测试(AVT)获得的数据频域分解(FDD)来识别。采用基于手动调节和基于参数化灵敏度的模型更新方法,表1还包括从已更新的有限元模型获得的计算固有频率[15,16]。测量频率与计算频率的一致性见表1。

《图1》

图1. 索罗克大桥。

《表1》

表1 索罗克大桥的模态频率

《2.2.从AVT中获取数据》

2.2.从AVT中获取数据

利用主跨中心部署的两个加速度计来测量采样频率为100Hz的桥面的垂直加速度。使用激光位移传感器同时测量桥面的相应位移,该传感器配备为用于监测桥梁运营行为的内置传感器。通过安装在桥面上的超声风速计记录风向和风速。将现有数据以10min为间隔来划分以进行深入调查。

尽管现场测试基本上是为了确保AVT数据的安全,但重型卡车也准备好创造高幅激励的运行条件。我们已经分析了悬架桥的一系列加速度,以确定TIV的性能。特别是考虑了两种类型的测试——卡车荷载测试和使用普通车辆的AVT,以便根据不同类型的车辆调查TIV的特性。卡车装载测试使用总重为25t、速度为25km·h–1和40km·h–1的三轴卡车进行。这座桥没有关闭正常交通,但与卡车相比,普通车辆的重量相对较轻,并且在调查卡车引起的振动方面没有很大的价值。

《2.3.AVT中的非平稳特征》

2.3.AVT中的非平稳特征

为了检验交通荷载引起的非平稳效应,采用符合经典OMA假设的稳定响应作为参考。为了从一组操作监测数据中选择平稳的响应,我们对一系列振动数据评估了峰度值。峰度值数据低于10,并且不包含峰值响应的可归类为平稳[17]。具体而言,选择最高风速来表示稳定数据以确保较高的信噪比。

表2总结了TIV数据统计。在表2中,由卡车荷载引起的TIV数据被称为“重型卡车”,而在正常操作期间测量的振动被分类为“普通车辆”。

《表2》

表2 TIV数据统计和相应的风速

图2所示的平稳应激数据因过往的车辆而没有极端的峰值,并且峰度值被识别为4.17,这表示平稳状态。功率谱密度(PSD)表明主要结构模态第一个(0.415Hz)和第三个竖向模态(0.842Hz)占据优势。

《图2》

图2.(a)测量的加速度;(b)稳定应激的PSD。

相反,如图3所示,当重型卡车通过桥时,2Hz附近的结构模态被放大。当一辆卡车进入主跨时,振幅逐渐增加,并引起非平稳响应。峰度值为12.47,其已经超过了平稳模态的标准。

《图3》

图3.(a)测量的加速度;(b)重型卡车的PSD。

如图4所示,与平稳模态相比,普通车辆也放大了较高的结构模态。4~5Hz的结构模态由较轻的车辆激发,其激发相对较高的频率。峰度值为49.69,超过了值为10的平稳标准,并指出了非平稳模态。

 《图4》

图4.(a)测量的加速度;(b)普通车辆的PSD。

图5表示使用短时傅里叶变换计算的时频分布。Hamming窗口用于计算谱段长度为213和50%重叠的谱图。图5清楚地表明,当一辆汽车穿过桥时,2~5Hz的较高结构模态变得极其活跃。

《图5》

图5.(a)时间频率分布;(b)测量的相对加速度(2012年7月14日, 早6 :00)。

《3.使用AM功能的信号平稳化》

3.使用AM功能的信号平稳化

《3.1.AM功能》

3.1.AM功能

Chiang和Lin[18]提出了一种从非平稳环境振动下收集到的结构响应数据来识别模态参数的方法。如图6所示,非平稳测量加速度a(t)可以模拟为平稳加速度ã(t)和AM函数Γ(t)的乘积,如下所示:

《图6》

图6. 利用AM函数的非平稳过程建模。(a)测量的加速度a(t);(b)AM函数C(t);(c)静止加速度ã(t)。

在这个模型中,平稳过程可以通过AM函数Γ(t)的调制从非平稳响应估算得出。AM函数可以从非平稳数据的时间均方根(RMS)[19]估算如下:

式中,C为遍历过程部分的平方根的期望值;Tw为用于计算时间RMS的片段长度。因此,可以获得近似的稳态加速度:

《3.2.AM函数在操作监控数据中的应用》

3.2.AM函数在操作监控数据中的应用

拟议的平稳化过程适用于从索罗克大桥获得的运行监测数据。AM函数是通过计算测得的加速度的时间均方根函数来评估,然后通过将测量的加速度除以AM函数来提取近似平稳加速度。

图7(a)、(b)分别显示了信号平稳化前后的交通引起的加速度。在信号平稳化过程之前,时间历史中测得的加速度具局部特征。峰度值超过了10的标准。然而,在稳定过程之后,由于交通负载而引起的峰值被从时间历史中移除。峰度值也从12.83降至2.36。

《图7》

图7.利用AM函数测量加速度。(a)信号稳定之前;(b)信号稳定之后。

图8比较了信号稳定之前和之后的归一化PSD。在信号平稳化期间,通过应用AM函数,车辆频率2~5Hz周围的频率分量明显减少。因此,平稳化数据的PSD显示出与平稳情况的PSD类似的频率分量。

《图8》

图8. AM函数信号平稳化之前和之后的测量加速度的归一化PSD。

《4.基于OMA的信号平稳化阻尼估计》

4.基于OMA的信号平稳化阻尼估计

《4.1.基于OMA阻尼估计的参数选择》

4.1.基于OMA阻尼估计的参数选择

将自然激励技术(NExT)-特征系统实现算法(ERA)应用于记录了4天(自2016年6月13日至6月16日)的1h加速度数据。为计算NExT中测量加速度的自相关函数程序,我们应用了一个长度为215、重叠率为50%的Hanning窗口[2]。通过设定此参数,平均数为21的频率分辨率为0.0031Hz。图9显示了结构的模拟脉冲响应函数(IRF)的计算相关函数。

如图9所示,基于计算出的IRF的形状选择汉克尔矩阵的大小。例如,从IRF正/负包络函数中确定计算的IRF降低至50%水平所需的时间(t1=17.77s,t2=18.08s),然后选择时间较长的一组(18.08s)。

《图9》

图9.基于包络函数确定IRF的使用时间长度。

基于OMA的阻尼估计对ERA的系统阶数非常敏感。为了获得收敛的阻尼比,应该在稳定图的辅助下选择一个系统的阶数,如图10所示。我们应该从1到50执行一系列的模态分析来检查所有的系统阶数,并且识别和记录满足频率范围、正阻尼和扩展模态振幅相干性(EMAC)等标准的第一个竖向模态(0.406Hz)。接下来,我们应该确定识别的阻尼比的中值。

《图10》

图10.选择系统顺序和识别阻尼比的过程。

为了选择合适的信号平稳化段长度,我们对从索罗克大桥获得的一天运行监测数据进行了参数研究。图11显示了估算阻尼比的平均值和方差系数(COV),根据分段比率定义,该分段比率定义为分段长度和目标固有周期(Tn)之间的比率。如图11所示,估测的阻尼比的平均值从0.5的分段比收敛到1%,而COV值从分段比为1.5开始变得稳定。因此,对于索罗克大桥使用移动RMS的信号分段过程,建议采用比目标自然周期长1.5倍的分段长度。

《图11》

图11.根据分隔比率统计索罗克大桥的阻尼估计值。(a)估测阻尼比的平均值;(b)估测阻尼比的COV。

《4.2.信号平稳化的估测阻尼比》

4.2.信号平稳化的估测阻尼比

在估算索罗克大桥第一模态的阻尼比时采用了三种比较方法。第一种情况没有考虑平稳化,只是将NExT-ERA用于测量数据;在下文中,这种情况被称为“未平稳化”。具有信号平稳化的识别阻尼比在下文中被称为“平稳化”,并且“组合”表示信号平稳化是根据测得的加速度的峰度值选择性地应用。正如Guo[17]所建议的那样,峰度值10被用作区分平稳和非平稳测量数据的标准。如果测得的加速度的峰度值大于10,则应用信号平稳化程序。从96个数据集中,有80个测得的加速度数据被确定为非平稳数据。如图11所示,对目标固有周期采用5倍的分段长度进行信号平稳化。

图12显示了索罗克大桥在运行过程中的估测阻尼比。估测的阻尼比根据从早到晚的时间段而波动,阻尼比看起来取决于交通量。在应用信号平稳化过程时,7月14日和15日所记录的严重分散现象减少。

《图12》

图12. 基于信号平稳化的索罗克桥的模态阻尼比估测。

在信号平稳后,第一竖向模态的阻尼比的平均值和标准偏差减小。信号平稳化后估计的阻尼比的COV为34%,与以前的AVT实践值相比是可接受的[20–24]。

《5.结语》

5.结语

本研究主要集中在韩国索罗克桥上的交通荷载造成的非平稳性以及它对基于OMA的阻尼估计的影响。文献调查和现场监测数据证实,TIV是局部的,导致了非平稳响应。车辆频率从2Hz到5Hz的结构模态被交通荷载放大,但较低的结构模态由稳定激励主导。我们提出了一种采用AM函数的信号平稳化过程。AM函数是从运行车辆引起的仅包络输出数据中获得的,并用于调制非平稳成分。通过信号平稳化程序成功地减少了由车辆引起的非平稳效应和频率失真值。信号平稳化后,峰度值显著降低,PSD函数与平稳数据近似。结果,TIV引起的过高/低估阻尼比被恢复,并且通过应用信号平稳化过程来减小阻尼比估计的不确定性。尽管本研究仅考虑了平稳化的AM函数,但在进一步的研究中可以对Hilbert变换等其他方法进行广泛的研究。

《Acknowledgements》

Acknowledgements

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIP) (2017R1A2B4008973) and also supported by grants (17SCIPB119960-02) from the Ministry of Land, Infrastructure and Transport of the Korean Government. This work was also partially supported by the BK21 PLUS research program of the National Research Foundation of Korea and the New Technology Research Center for Bridges through the Institute of Engineering Research at Seoul National University. The authors are also grateful to the owner of the Sorok Bridge, TM E&C, and to the Korea Infrastructure Safety and Technology Corporation (KISTEC) for sharing operational monitoring data.

《Compliance with ethics guidelines》

Compliance with ethics guidelines

Sunjoong Kim and Ho-Kyung Kim declare that they have no conflict of interest or financial conflicts to disclose.