大跨度斜拉桥下击暴流风致振动响应实测

刘志文，李书琼，刘勇，许映梅，陈政清

摘   要：以蘇通长江公路大桥为工程背景，针对该桥风致振动响应监测系统实测的一次下击暴流风与桥梁结构振动加速度响应实测数据，对该桥在一次雷暴天气下风速、风向及主梁振动响应进行研究. 首先，对桥位处下击暴流实测风速、风向数据进行分析，获得了该桥主梁跨中、桥塔塔顶处下击暴流风的时变平均风与脉动风特性;然后，对下击暴流作用下主梁风致振动加速度响应数据进行分析. 结果表明：在下击暴流作用下，该桥主梁与塔顶高度处风速发生了明显突变，持续时间约为10～24 min;主跨跨中主梁外侧边缘处下游、上游侧最大瞬时风速分别为32.4 m/s和27.3 m/s，南、北桥塔塔顶高度处最大瞬时风速分别达60.5 m/s和62.9 m/s. 主梁高度处30 s时距湍流度约0.048～0.32，10 min时距湍流度约0.43～0.51;主梁下游与北塔处折减脉动风速符合高斯特性，其功率谱与Burlando等学者的实测结果吻合较好. 主梁跨中附近（即NJ26D、NJ32D拉索锚固处）发生了较为明显的短时竖向与横桥向振动，相应加速度响应幅值分别为0.25 m/s2和0.10 m/s2，对应位移幅值分别为0.12 m与0.03 m;主梁竖向振动响应明显大于横桥向振动响应，主梁竖向振动主频为0.183 Hz，与主梁全桥一阶正对称竖弯振型频率0.174 Hz接近;横桥向振动主频为0.117 Hz，与主梁全桥一阶正对称侧弯振型频率0.097 5 Hz接近.

关键词：大跨度斜拉桥;下击暴流;风特性;主梁风致振动;现场实测

中图分类号：U446.2                                 文献标志码：A

Field Measurement of Wind-induced Vibration Response

of Long-span Cable-stayed Bridge under Downburst

LIU Zhiwen1，2？覮，LI Shuqiong1，4，LIU Yong3，XU Yingmei3，CHEN Zhengqing1，2

（1. Hunan Provincial Key Laboratory for Wind Engineering & Bridge Engineering （Hunan University），Changsha 410082，China;

2. College of Civil Engineering，Hunan University，Changsha 410082，China;

3. Jiangsu Sutong Bridge Co，Ltd，Nantong 226001，China;

4. China Railway Siyuan Survey and Design Group Co，Ltd，Wuhan 430063，China）

Abstract：Taking Sutong Yangtze River Highway Bridge （STB） as an engineering background， a downburst wind and the vibration acceleration responses of the bridge structure were monitored based on the wind-induced vibration response monitoring system of the bridge. Furthermore， the wind speed and direction， vibration responses of the main girder under the thunderstorm were investigated. Firstly， the measured wind speed and wind direction data of the downburst wind at the bridge site were analyzed， the time-varying mean wind and fluctuating wind characteristics of the downburst at the mid-span of the main girder and the pylon tops were obtained. Then， the wind-induced vibration acceleration responses of the main girder under the downburst were analyzed. The research results show that the wind speeds at the main girder level and the pylon tops of the STB changes abruptly under downburst， which lasts approximately 10～24 minutes. The maximum instantaneous wind speeds at the leeward and windward side of the main girder in the middle of the main span of STB are 32.4 m/s and 27.3 m/s， respectively. Meanwhile， the maximum instantaneous wind speeds at the south and north pylon tops are 60.5 m/s and 62.9 m/s， respectively. The turbulence intensity for a time interval of 30 s at the main girder level is about 0.048～0.32， and the turbulence intensity for a time interval of 10 minutes is about 0.43～0.51. The reduced fluctuating wind speed at the downstream of the main girder and the north pylon conforms to Gaussian characteristics， and its power spectrum is in good agreement with the field measurement results by Burlando and other scholars. Obvious short-term vertical and horizontal vibration occurred near the middle of the main girder （that is， the anchorage of NJ26D and NJ32D cables）， and the corresponding acceleration response amplitudes were 0.25 m/s2 and 0.10 m/s2， corresponding to the displacement amplitudes of 0.12 m and 0.03 m， respectively. The vertical vibration response of the main girder is obviously larger than the transverse vibration response. The predominant frequency of vertical vibration of the main girder is 0.183 Hz， which is close to the 1st symmetric vertical bending mode frequency 0.174 Hz of the main girder. And the predominant frequency of horizontal vibration is 0.117Hz， which is close to the 1st symmetry lateral bending mode frequency 0.097 5 Hz of the main girder.

Key words：large-span cable-stayed bridge;downburst;wind characteristics;wind-induced vibration of main deck;field measurement

下击暴流是一种雷暴云中局部强下沉气流在到达地面后产生的直线型大风，在接近地面处风速达到最大，具有突发性、局部性和随机性等特点[1]. 我国是下击暴流多发国家之一，在全国较大范围内均有发生的可能性[2]. 下击暴流对工程结构影响较大，可能引起建筑结构、输电线塔、桥梁结构附属设施等破坏[3-5]. 因此，开展下击暴流风特性及其对工程结构影响的研究具有十分重要的意义.

国内外许多学者针对下击暴流风特性及其对结构的影响开展了大量的研究工作，主要工作有现场实测[6-9]、数值模拟[10-12]、风洞试验[13-18]和理论计算[19-22]等. 在实测研究方面，Choi[6]在一座高为150 m的塔上设置了5个观测层，对50多次雷暴的风速剖面进行了实测研究. Burlando等[7]对2012年10月在意大利利沃诺科斯特观测到的一次下击暴流进行了风场特性分析. Solari等[8]对地中海北部港口6年间发生的277个下击暴流风速记录进行了系统分析，获得了下击暴流风特性. Stengel等[9]对德国北部的一条输电线路进行实测，观测到悬索塔导线在一次下击暴流下的实测响应，并与有限元模型的时域模拟进行了比较. 在数值模拟研究方面，Wood等[10]进行了基于湍流模型的下击暴流计算流体力学数值模拟，其结果与下击暴流喷射装置试验结果吻合较好. Chay等[11]采用CFD数值模拟方法模拟了下击暴流平均风，重点考察了风速随下击暴流成熟和衰减强度的变化规律. 刘志文等[12]采用二维数值模拟方法在边界层风洞中设置倾斜平板进行了数值模拟研究，结果表明边界层风洞中设置倾斜平板可有效模拟下击暴流水平风速剖面. 在试验研究方面，曹曙阳等[13]在日本宫崎大学大型多风扇主动控制风洞中实现了雷暴冲击风模拟. Jesson等[14]为研究建筑结构在下击暴流瞬态气流作用下的压力分布，研发了下击暴流瞬态风场模拟试验装置，其试验结果表明采用该装置模拟的瞬态风速时程与实测下击暴流瞬态风速时程吻合较好. Aboutabikh等[15]设计并制造了两层带叶片的百叶窗下击暴流模拟试验装置，在风洞中模拟了下击暴流风场. 辛亚兵等[16]基于传统大气边界层风洞开发了下击暴流模拟装置，并对下击暴流风作用下大跨连续刚构桥最大双悬臂状态风致振动响应进行了试验研究. Elawady等[17]采用WindEEE多功能风洞模拟了下击暴流风场，并进行了多跨输电线路下击暴流风致振动响应气弹模型风洞试验研究. Junayed等[18]采用WindEEE多功能风洞模拟了缩尺比较大的下击暴流风场特性，并将试验模拟的下击暴流平均风场和脉动风场特性与实测下击暴流风特性进行了比较，两者吻合相对较好. 在理论计算方面，Chen等[19]提出了一种混合随机模型模拟下击暴流风速时程，即下击暴流平均风由Wood风剖面和Holmes时间函数得到，脉动风速由随时间变化的幅值调幅函数和服从标准正态分布的高斯随机过程得到. Hao等[20]采用基于冲击射流模型的CFD数值仿真技术模拟了下击暴流稳态风场和瞬态风场，分析了桥梁在模拟下击暴流作用下抖振时域响应，结果表明，下击暴流脉动风速对桥梁响应影响相对较小. 辛亚兵等[21]以赤石大桥桥址处实测下击暴流风时程数据为基础，采用谐波叠加法模拟桥址区脉动风速并加以调制，实现了桥址区下击暴流风时程模拟，对下击暴流作用下大跨度斜拉桥施工状态静风响应和非线性时域抖振响应进行了计算.

综上所述，目前国内外学者对下击暴流风场特性的观测与试验研究相对较多，而对下击暴流作用下的结构响应现场实测研究则相对较少，因此进行下击暴流风特性及其风致振动响应实测研究具有十分重要的价值和意义. 本文依托苏通大桥风致振动监测系统，对大跨度斜拉桥风致振动响应特性进行为期2年的现场实测研究，监测到2019年4月19日桥位处发生了一次突发大风，监测系统完整记录了本次突发大风风速、风向及主梁振动加速度响应数据，为大跨度斜拉桥下击暴流风效应研究积累了十分宝贵的实测数据. 本文重点对本次下击暴流大风的风特性与主梁风致振动响应特性进行分析.

1   苏通大桥风致振动监测系统

苏通长江公路大桥位于江苏省南通市和苏州市之间，是国家重点干线公路沈海高速（G1）跨越长江的重要通道. 据设计资料可知苏通大桥主桥为主跨1 088 m的双塔双索面斜拉桥，其跨径布置为100 m+100 m+300 m+1 088 m+300 m+100 m+100 m=2 088 m. 主梁采用闭口流线型钢箱梁，梁宽41.0 m，梁高4.0 m;斜拉索最大长度为577 m;采用倒Y形桥塔，塔高300.4 m. 考虑到大桥规模与运营期内的维护需要，该桥建成后安装了结构健康监测系统（Structural Health Monitoring System，SHMS）[23]. 该系统可对风、温度、车辆荷载及腐蚀作用等进行监测，并对桥梁结构的位移、索塔和桥墩倾斜度、支座位移、拉索索力和结构应变等进行监测.

为确保大桥在风荷载作用下安全运营，在苏通大桥结构健康监测系统基础上，建立了苏通大桥风致振动监测系统，对大桥部分拉索锚固处主梁竖向、横向振动加速度响应进行监测. 风速仪布置于主桥南、北塔塔顶，主跨跨中桥面上、下游距离桥面2.28 m高度处. 风向角0°对应北风，90°对应东风，采样频率为1 Hz（实际测量中开启二维模式，只采集水平风速、风向），由于苏通大桥桥轴线与正北方向有10.6°偏角，故需对风向角进行修正以得到以桥轴线为参考的风偏角. 在主梁上NJ26D与NJ32D号拉索锚固处分别布置了竖向和横桥向加速度传感器（型号为941B），以監测主梁竖向与横桥向振动加速度响应，采样频率设置为100 Hz. 图1所示为苏通大桥风致振动响应监测系统传感器布置示意图，图2所示为苏通大桥主梁横断面及风速仪位置示意图.

2   桥址处风特性分析

2.1   风速风向时程

2019年4月19日，苏通大桥桥位处出现了一次显著的大风天气. 图3所示为苏通大桥风致振动监测系统所记录的主跨跨中桥面上、下游及桥塔塔顶处风速风向时程. 表1所示为各风速监测点位置风速突变时刻风速、风向参数统计汇总. 结合图3及表1结果可知，在4月19日中午12：19时，南塔塔顶风速首先达到最大瞬时风速60.5 m/s;在12：55北塔塔顶风速出现最大瞬时风速62.9 m/s;12：57—12：58时主跨跨中桥面下游、上游风速分别达到最大值32.4 m/s和27.3 m/s;各监测点处突变大风持续时间约为10～24 min，与文献[24]所给出的单个微下击暴流平均持续时间约为13 min较为接近;在下击暴流发生前较长时间内，主跨跨中主梁上、下游侧及南北塔塔顶的风向变化趋势基本一致，即由南风变为西风，并在较长时间内风向保持不变，在下击暴流发生时段内各监测点风向发生了较为明显的突变. 综合以上现象可看出，本次大风天气具备下击暴流的气象特征，初步判断应为一次下击暴流.

由表1可知，大桥不同监测点处风速开始突变时刻以及峰值时刻存在较大差异，南塔出现下击暴流时刻较其他测点要提前30 min，其原因可能是此次下击暴流的尺度较小，影响范围有限，且其中心可能处于移动中. 由图3中风向与表1中时间参数可以推断出，4月9日出现的下击暴流位于苏通大桥上游侧，且其中心自南向北缓慢移动.

依据图3中风速可知，主梁上游风速在13：13以及13：45均存在特异性大风，此时风速可达68 m/s，远大于下击暴流峰值风速. 结合上游风向时程可看到，上述2个时刻风向突然转变，由-100°（垂直于桥轴向）突变至-168.6°与179.1°（基本沿着桥轴向）. 可以判断出，这2个时刻存在着沿桥轴线方向的强局部气流对上游风速造成干扰，并且其流向与桥轴线平行，故未对下游风速产生影响. 此外，也不排除13：10之后仪器出现短暂故障，具体原因有待进一步研究. 本文重点研究12：49—13：09发生的下击暴流风特性，故此次特异数据基本没有影响. 同时为了保证数据的真实可靠以及对比不同高度处下击暴流风特性，后文主要对主梁下游以及北塔塔顶风速进行分析.

2.2   时变平均风与脉动风速

由图3可知，在下击暴流时段主跨跨中主梁处、桥塔塔顶等风速为非平稳风速时程，故参考文献[25-26]中瞬态风信号的经典分解规则. 将瞬时风速U（t）分解为时变平均风速U（t）与非平稳脉动风速u（t）：

U（t） = U（t） + u（t）               （1）

式中：U（t）为t时刻的时变平均风速，m/s;u（t）为t时刻的脉动风速，m/s.

脉动风速u（t）为一非平稳随机过程，可表示为

u（t） = σu（t）·u′（t）                （2）

式中：u′（t）为折减脉动风速，即为平稳高斯随机过程;σu（t）为脉动风速u（t）的缓变根方差，即

σu（t） = ηU（t）              （3）

即

η = ■ = Iu（t）           （4）

式中：Iu（t）為时变湍流度.

为进一步研究下击暴流作用下主跨跨中桥面高度处以及桥塔顶部高度处的风速特性，需要对其进行风速分解. 其中时变平均风速的提取采用文献[25]中的滑动平均法，滑动平均风速定义如下：

Urm（j） = ■■U（t）dt           （5）

式中：Urm（j）为j时刻的瞬时风速;Trm为滑动平均的时间间隔，结合文献[7]中Burlando的建议取值，本文中Trm = 30 s.

图4所示分别为下击暴流作用下苏通大桥主跨跨中桥面高度处下游与北塔塔顶的时变平均风速和脉动风速分析结果.

由图4（a）（b）对比可知，在下击暴流时段，北塔塔顶的时变平均风速比主跨跨中桥面下游侧时变平均风速变化更为剧烈，且塔顶风速在下击暴流出现后平均风速比发生前增大数倍. 两者的脉动风速变化规律也显著不同，主梁跨中下游侧在下击暴流出现前脉动风速低，当其过境时脉动风速峰值为13.7 m/s，过境后脉动风速仍保持较大值;北塔塔顶脉动风速变化却与之相反，下击暴流出现前脉动风速已达10 m/s以上，下击暴流发生时达到16.2 m/s，过境后反而较低. 可以看到下击暴流对苏通大桥不同高度处风特性的影响存在较大差异.

除下击暴流风速突变特征外，下击暴流作用下主梁下游与北塔塔顶的风向变化规律也值得进一步分析. 故采用与式（5）一致的滑动平均法求取30 s时变平均风向，图5所示为下击暴流作用下苏通大桥主梁下游与北塔塔顶的瞬时风向和时变平均风向. 由图5可知，下击暴流发生前主梁下游主导风向约165°，北塔风向波动剧烈;下游风向突变后稳定至-96°～-110°，北塔风向为-100°，基本垂直于桥轴线;下击暴流衰退后，下游瞬时风向与时变平均风向差异显著，北塔风向保持稳定. 此次下击暴流未发生时风向变化剧烈，发生时风向稳定，发生后不同高度处风向变化规律不同.

2.3   湍流度

湍流度是描述脉动风速变化大小的一个重要的参数，为脉动风速根方差与平均风速之比. 考虑到下击暴流时空尺度小，具有突发性，且持续时间短暂，参考文献[7]中处理方法，选取时距30 s采用式（4）计算湍流度. 为便于比较，针对下击暴流时段风速分别取平均时距为T1 = 30 s和T2 = 10 min进行顺风向湍流度计算，图6所示为主跨跨中主梁高度处下游与北塔塔顶处湍流度随时间变化曲线.

由图6可知，当取平均时距为T1 = 30 s时，在下击暴流时段主跨跨中桥面高度处下游侧顺风向湍流度约为Iu（t）=0.048～0.32，北塔塔顶湍流度约为Iu（t） = 0.01 ～ 0.014. 当取平均时距为T2 = 10 min时，在下击暴流时段主跨跨中桥面高度处下游侧顺风向湍流度约为Iu（t） = 0.43 ～ 0.51，北塔塔顶湍流度约为Iu（t） = 0.31 ～ 0.48. 而文献[27]中台风作用下苏通大桥主梁高度处湍流度为0.10 ～ 0.30，可见下击暴流作用下主梁高度处湍流度总体上略大于台风作用下的湍流度. 主要原因是下击暴流风速突变较台风而言持续时间短且突变风速较大.

2.4   折减脉动风特性

由上文可知下击暴流风场中的脉动分量与常规台风气象差距较大，现着重对其脉动分量进行研究. 仍以30 s为基本时距，利用式（2）计算下击暴流过境时主跨跨中主梁桥面高度處下游与北塔塔顶处时变风速中的折减脉动风速成分，绘制折减脉动风速随时间变化曲线如图7（a）（d）所示. 由图7可知，折减脉动风速呈现出经典随机平稳高斯特性. 通过图7（b）（e）的频率直方图可见，虽然偏斜度不完全为0，峰度不为3，但是图形与参考高斯拟合曲线有良好的一致性，证实了这两个测点的折减脉动风的高斯特性. 按式（6）计算折减脉动风速的概率密度p（u′）.

p（u′） =  ■e■          （6）

式中：u′为脉动风速;σu′表示脉动风速标准差. 为研究折减脉动风速的频率成分特征，将实测顺风向脉动风谱与Von-Karman谱、Davenport谱、Simiu谱进行对比.

Von-Karman谱是1948年美国著名空气动力学专家Von-Karman提出的自由大气水平脉动风谱：

■ = ■        （7）

式中： f = nLxu/U，n为风的脉动频率，Lxu是湍流积分尺度，U是平均风速;β是摩擦因数;u*为摩阻速度.

Davenport谱是1961年加拿大著名风工程专家Davenport提出的自由大气水平脉动风谱：

■ = ■        （8）

式中： f = 1 200 n/U（10），U（10）是z = 10 m高度处的平均风速.

Simiu谱为1974年美国学者Simiu提出的与高度有关的水平脉动风速谱：

■ = ■        （9）

式中：f = nz/u（z），其中z为测点高度.

由图7（c）（f）可见，跨中桥面高度处下游与北塔塔顶的折减脉动风速功率谱趋势相同. 对比实测谱线与经验谱线可知，3种经验谱线与实测谱总体趋势一致，但Davenport谱与实测谱线在低频段拟合较好，高频差异较大;Simiu谱相反，低频差异大，高频吻合度高;Von-Karman谱曲线变化介于两者之间. 出现这种现象的原因主要是下击暴流风自身的非平稳特性、数据处理时采用非平稳风速分解模型以及各经验谱的特有参数与适用性的差异.

本文中实测谱线与文献[7]中Burlando等在2012年观测到的意大利利沃诺海岸遭遇的下击暴流风场下的折减脉动风速功率谱密度函数曲线相比，两者功率谱在曲线的变化趋势以及峰值大小上具有较好的相似性. 且两者低频段均较低，下降段符合经典的天气型风惯性子区间的n-5/3斜率下降.

3   主梁振动响应

3.1   主梁加速度响应时程

为研究苏通大桥主梁在风荷载作用下的振动特性，以保证大桥在运营期间安全运行，分别在NJ32D与NJ26D号拉索与主梁锚固处安装加速度传感器，监测主梁竖向、横桥向加速度响应. 图8所示为2019年4月19日12：00—14：00主梁在NJ32D、NJ26D拉索锚固处竖向、横桥向加速度响应时程曲线.

由图8可知，NJ32D和NJ26D号拉索锚固处主梁在该日中午12：49—13：09附近均发生了一次加速度响应较大的短时振动现象，主梁在NJ32D拉索锚固处的加速度响应与主梁在NJ26D拉索锚固处的加速度响应总体较为接近;主梁在这两处的竖向和横桥向最大加速度响应幅值分别约为0.25 m/s2和0.10 m/s2，主梁竖向加速度响应幅值约为横桥向加速度响应幅值的2.5倍.

为进一步了解下击暴流发生时主梁结构振动响应特征，对12：49—13：09时段主梁在NJ32D和NJ26D拉索锚固处的振动加速度响应进行时频分析. 考虑到下击暴流作用下的结构振动响应具有明显非平稳性，故本文采用可考虑信号时变特性的连续小波变换分析方法对主梁振动响应进行分析.

在小波变换中，一组形状相似的小波基函数是由母小波经过伸缩与移动进行转化后得到的. 任意信号Ψ（t）的小波变换就是把信号Ψ（t）以小波函数为基底将其展开，再把信号以小波函数的线性组合表示出来. 对于任意信号Ψ（t）∈L2（R）（其中L2（R）为能量有限的信号空间），Ψ（t）的傅里叶变换为■（ω），当■（ω）满足条件：

CΨ = ■■dω < ∞        （10）

时，则称Ψ（t）为一个基本小波或母小波，同时将Ψ（t）经过平移与伸缩后可以得到：

■a，b（t）=■Ψ■，a，b∈R;a≠0  （11）

称其为一个小波系列，称a为伸缩因子，b为平移因子.

任意信号f（t）（f（t）∈I2（R））的连续小波变换（Continue Wavelet Transform，CWT）表达式为：

CWTf （a，b） = 〈 f（t），Ψa，b（t）〉 =

■■f（t）Ψ*■dt           （12）

采用MATLAB中小波工具箱对主梁振动加速度数据做连续小波变换并绘制时频图. 选取CMOR小波为小波基函数，其中参数定义：带宽fb = 1.5 Hz，中心频率fc = 3 Hz. 图9、图10分别所示为NJ32D、NJ26D锚固处下击暴流作用时段主梁竖向、横桥向加速度响应不同时刻频谱特征.

同时为进一步分析下击暴流作用下苏通大桥的振动模态，使用有限元计算软件ANSYS建立全桥有限元模型进行动力特性分析，表2所示为苏通大桥主桥结构成桥状态动力特性分析结果汇总. 由表2及图9可知，下击暴流发生时段主梁在NJ32D号拉索锚固处竖向振动加速度响应主频为0.183 Hz，与苏通大桥主桥结构成桥状态一阶正对称竖弯频率f3 =0.174 4 Hz较为接近;由图10可知，下击暴流发生时段主梁NJ26D号拉索锚固处竖向振动加速度响应主频为0.183 Hz和0.217 Hz，与苏通大桥主桥结构成桥状态主梁一阶正对称竖弯频率f3=0.177 4 Hz和主梁结构一阶反对称竖弯频率f4=0.214 6 Hz接近. 综合表2及图9～图10可知，主梁在NJ32D和NJ26D处横桥向振动加速度响应主频均为0.116 7 Hz，与主梁一阶正对称侧弯频率f2=0.097 5 Hz较为接近. 综上可知，在下击暴流风作用下，主梁竖向振动响应明显大于横桥向振动响应，且均以低阶竖弯振型为主.

3.2   主梁位移响应

考虑到主梁竖向、横向加速度响应以低频为主，故对主梁加速度響应时程进行积分，以得到主梁竖向、横桥向位移响应时程. 由逆傅里叶变换公式，可将任意频率下的加速度信号的傅里叶分量表示为：

a（t） = Ae jωt         （13）

式中：a（t）为加速度信号在频率ω的傅里叶分量;A为对应a（t）的系数;j为虚数，即■.

当初速度与初位移分量均为0时，对加速度信号的傅里叶分量两次积分可得出位移分量：

x（t） =■■a（λ）dλdτ =■Ve jωtdτ = Xe jωt

（14）

式中：x（t）为速度信号在频率ω的傅里叶分量;X为对应x（t）的系数.

计算加速度信号的二次积分的数值公式为：

y（r）=■-■H（k）X（k）e j2πkr/N     （15）

H（k）=1， fd ≤ kΔf ≤ fu0，其他       （16）

式中：fd和fu分别为下限截止频率与上限截至频率，Hz，此处取fd = 0.1 Hz、 fu = 50 Hz;X（k）为x（r）的傅里叶变换;Δf 为频率分辨率，Hz.

图11所示为采用该方法所计算得到的主梁在NJ32D拉索、NJ26D拉索锚固处的竖向和横桥向位移响应时程曲线. 由图11可知，在下击暴流时段，主梁在NJ32D拉索锚固处竖向和横桥向位移响应最大幅值分别为0.111 m和0.027 m;主梁在NJ26D拉索锚固处竖向和横桥向位移响应最大幅值分别为0.116 m和0.020 m. 可见在下击暴流作用下，主桥结构主梁振动响应明显增大.

4   结   论

依托苏通大桥结构健康监测系统和风致振动响应监测系统，对2019年4月19日的一次突发大风风速、风向和主梁风致振动加速度响应数据进行分析，得到如下主要结论：

1）2019年4月19日12：19，苏通大桥桥位处出现了瞬时突发大风，南塔、北塔塔顶最大瞬时风速分别为60.5 m/s和62.9 m/s，主跨跨中桥面高度处下游、上游侧最大瞬时风速分别为32.4 m/s和27.3 m/s，突变大风持续时间约为10～24 min;各监测点风向存在较为明显的突变现象;不同监测点处风速最大值达到时刻存在一定差别，初步判断该突发大风应为一次下击暴流现象.

2）在下击暴流时段（约在13：00左右），当取平均时距为T1=30 s时，主跨跨中桥面高度处下游侧顺风向湍流度约为Iu（t）=0.048～0.32，北塔塔顶顺风向湍流度约为Iu（t）=0.01～0.014;当取平均时距为T2=10 min时，下击暴流时段内主跨跨中桥面高度处下游侧顺风向湍流度约为Iu（t）=0.43～0.51，北塔塔顶顺风向湍流度约为Iu（t）=0.31～0.48. 主梁下游与北塔塔顶处折减脉动风速符合高斯特性，其功率谱与Burlando等学者[7]的实测结果吻合较好.

3）在下击暴流作用下，苏通大桥主梁主跨跨中竖向振动响应明显大于横桥向振动响应，且均以低阶振动响应为主. 主梁竖向位移约0.12 m，横桥向位移0.03 m，下击暴流作用下大跨度斜拉桥主梁振动响应明显增大.

参考文献

[1]    FUJITA T T，BYERS H R. Spearhead echo and downburst in the crash of an airliner[J]. Monthly Weather Review，1977，105（2）：129—146.

[2]    李宏海，欧进萍. 我国下击暴流的时空分布特性[J]. 自然灾害学报，2015，24（6）：9—18.

LI H H，OU J P. Spatiotemporal distribution characteristics of downburst in China[J]. Journal of Natural Disasters，2015，24（6）：9—18. （In Chinese）

[3]    HAWES H. Review of recent Australian transmission line failures due to high intensity winds[R]. Buenos Aires：[s.n.]，1993.

[4]    谢强，张勇，李杰. 华东电网500 kV任上5237线飑线风致倒塔事故调查分析[J]. 电网技术，2006，30（10）：59—63.

XIE Q，ZHANG Y，LI J. Investigation on tower collapses of 500 kV Renshang 5237 transmission line caused by downburst[J]. Power System Technology，2006，30（10）：59—63. （In Chinese）

[5]    刘志文，丁冬，贾亚光，等. 山区大跨斜拉桥悬臂施工期风振响应实测[J]. 湖南大学学报（自然科学版），2017，44（1）：16—22.

LIU Z W，DING D，JIA Y G，et al. Field measurement of wind-induced vibration responses of a large cable-stayed bridge during construction stages in mountainous terrain[J]. Journal of Hunan University （Natural Sciences），2017，44（1）：16—22. （In Chinese）

[6]    CHOI E C C. Field measurement and experimental study of wind speed profile during thunderstorms[J]. Journal of Wind Engineering and Industrial Aerodynamics，2004，92（3/4）：275—290.

[7]    BURLANDO M，ROMANI■ D，SOLARI G，et al. Field data analysis and weather scenario of a downburst event in Livorno，Italy，on 1 October 2012[J]. Monthly Weather Review，2017，145（9）：3507—3527.

[8]    SOLARI G，BURLANDO M，DE GAETANO P，et al. Characteristics of thunderstorms relevant to the wind loading of structures[J]. Wind and Structures，2015，20（6）：763—791.

[9]    STENGEL D，THIELE K. Measurements of downburst wind loading acting on an overhead transmission line in Northern Germany[J]. Procedia Engineering，2017，199：3152—3157.

[10]  WOOD G S，KWOK K C S，MOTTERAM N A，et al. Physical and numerical modelling of thunderstorm downbursts[J]. Journal of Wind Engineering and Industrial Aerodynamics，2001，89（6）：535—552.

[11]  CHAY M T，ALBERMANI F，WILSON R. Numerical and analytical simulation of downburst wind loads[J]. Engineering Structures，2006，28（2）：240—254.

[12]  劉志文，陈以荣，辛亚兵，等. 基于边界层风洞的下击暴流稳态风场特性数值模拟[J]. 湖南大学学报（自然科学版），2020，47（7）：10—20.

LIU Z W，CHEN Y R，XIN Y B，et al. Numerical simulation on steady wind field characteristics of downburst based on atmosphere boundary layer wind tunnel[J]. Journal of Hunan University （Natural Sciences），2020，47（7）：10—20. （In Chinese）

[13]  CAO S Y，NISHI A，KIKUGAWA H，et al. Reproduction of wind velocity history in a multiple fan wind tunnel[J]. Journal of Wind Engineering and Industrial Aerodynamics，2002，90（12/13/14/15）：1719—1729.

[14]  JESSON M，STERLING M. A simple vortex model of a thunderstorm downburst—a parametric evaluation[J]. Journal of Wind Engineering and Industrial Aerodynamics，2018，174：1—9.

[15]  ABOUTABIKH M，GHAZAL T，CHEN J X，et al. Designing a blade-system to generate downburst outflows at boundary layer wind tunnel[J]. Journal of Wind Engineering and Industrial Aerodynamics，2019，186：169—191.

[16]  辛亞兵，刘志文，邵旭东，等. 大跨连续刚构桥下击暴流作用效应试验研究[J]. 中国公路学报，2019，32（10）：279—290.

XIN Y B，LIU Z W，SHAO X D，et al. Effects of downburst on long-span continuous rigid frame bridges[J]. China Journal of Highway and Transport，2019，32（10）：279—290. （In Chinese）

[17]  ELAWADY A，ABOSHOSHA H，EL DAMATTY A，et al. Aero-elastic testing of multi-spanned transmission line subjected to downbursts[J]. Journal of Wind Engineering and Industrial Aerodynamics，2017，169：194—216.

[18]  JUNAYED C，JUBAYER C，PARVU D，et al. Flow field dynamics of large-scale experimentally produced downburst flows[J]. Journal of Wind Engineering and Industrial Aerodynamics，2019，188：61—79.

[19]  CHEN L Z，LETCHFORD C W. A deterministic-stochastic hybrid model of downbursts and its impact on a cantilevered structure[J]. Engineering Structures，2004，26（5）：619—629.

[20]  HAO J M，WU T. Nonsynoptic wind-induced transient effects on linear bridge aerodynamics[J]. Journal of Engineering Mechanics，2017，143（9）：04017092.

[21]  辛亚兵，刘志文，邵旭东. 实测下击暴流作用下大跨斜拉桥抖振响应分析[J]. 空气动力学学报，2020，38（4）：705—714.

XIN Y B，LIU Z W，SHAO X D. Buffeting response analysis of long-span cable-stayed bridge under the measured downburst[J]. Acta Aerodynamica Sinica，2020，38（4）：705—714. （In Chinese）

[22]  牛华伟，洪飞，欧阳克俭，等. 拉线门塔在下击暴流作用下的响应分析[J]. 湖南大学学报（自然科学版），2018，45（3）：72—81.

NIU H W，HONG F，OUYANG K J，et al. Analysis on downburst induced response of guyed portal tower[J]. Journal of Hunan University （Natural Sciences），2018，45（3）：72—81. （In Chinese）

[23]  余波，邱洪兴，王浩，等. 苏通大桥结构健康监测系统设计[J]. 地震工程与工程振动，2009，29（4）：170—177.

YU B，QIU H X，WANG H，et al. Health monitoring system for Sutong Yangtze River bridge[J]. Journal of Earthquake Engineering and Engineering Vibration，2009，29（4）：170—177. （In Chinese）

[24]  HJELMFELT M R. Structure and life cycle of microburst outflows observed in Colorado[J]. Journal of Applied Meteorology，1988，27（8）：900—927.

[25]  CHOI E C，HIDAYAT F A. Dynamic response of structures to thunderstorm winds[J]. Progress in Structural Engineering and Materials，2002，4（4）：408—416.

[26]  TAO T Y，WANG H，WU T. Comparative study of the wind characteristics of a strong wind event based on stationary and nonstationary models[J]. Journal of Structural Engineering，2017，143（5）：04016230.

[27]  王浩，李爱群，谢静，等. 台风作用下超大跨度斜拉桥抖振响应现场实测研究[J]. 土木工程学报，2010，43（7）：71—78.

WANG H，LI A Q，XIE J，et al. Field measurement of the buffeting response of a super-long-span cable-stayed bridge under typhoon[J]. China Civil Engineering Journal，2010，43（7）：71—78. （In Chinese）