1、：钢 箱 梁 悬 索 桥 风 致 响 应 与 静 力 分 析收稿日期：作者简介：文海泉（），男，在读硕士，从事大跨度桥梁抗风分析研究；马小淇（），男，在读硕士，从事桥梁裂缝分析研究 通信作者：江阿兰，女，博士，教授，从事桥梁结构损伤诊断与旧桥加固技术研究文海泉，马小淇，江阿兰（大连交通大学土木工程学院，辽宁 大连）摘 要：风载和静载是桥梁在成桥阶段最常见的两种荷载，而钢箱梁悬索桥在风致作用下的响应尤为剧烈，容易引发安全性问题。为分析其风致响应情况，对某流线型钢箱梁悬索桥进行了数值风洞研究，计算了桥梁结构的三分力系数，并以此为基础分析了桥梁的风致响应，得出了最大的内力与位移，验证了悬索桥的安全性
2、。研究结果表明：流线型钢箱梁三分力系数在 攻角中稳定变化；桥梁在驰振、颤振与涡振中均表现出良好稳定性；全桥在不同荷载工况下的位移与内力均符合规范要求，验证了桥梁结构的安全性与稳定性。关键词：钢箱梁悬索桥；静力三分力系数；数值风洞中图分类号：文献标识码：文章编号：（）现代桥梁不断朝着大跨度和轻质量的方向快速发展，导致桥梁结构刚度和阻尼不断降低，从而大大增加了结构对风荷载的敏感性。钢箱梁悬索桥同时具备这两样优点，如今在跨江大桥中得到广泛应用，随之而来我们更需要考虑风荷载对桥梁的影响，风荷载已经成为悬索桥结构设计和施工的主要控制荷载。风经过桥梁结构时，会将风的动能转化为桥梁的能量，造成了桥梁振动，反
3、之，振动起来的桥梁又会影响空气的流动，改变空气作用力，形成了风与结构的相互作用。例如，年美国塔科马大桥设计不当发生颤振，年虎门大桥因连续设置水马，改变了钢箱梁的气动外形，引起了涡振现象。主梁的选型与尺寸，整体桥梁的跨度都对桥梁抗风性能有着至关重要的影响。以某座钢箱梁悬索桥为例，采用有限元软件对主梁进行数值风洞模拟，分析桥梁风致响应情况，对整桥进行静力分析，验证桥梁的安全性。桥梁模型桥梁模型为（）地锚式两塔连续钢箱梁悬索桥。加劲梁共有两对竖向支座，设于两个锚碇上。主缆中跨矢跨比为 ，主梁中心距为 ，加劲梁是闭口扁平流线型钢箱梁，高 ，宽 （含风嘴）。桥梁总体布置和加劲梁断面尺寸见图，图。图 1桥
4、梁总体布置主缆和吊索设置为节点有三个线自由度的索单元，用于模拟仅能受拉的杆件。主缆锚固处和索塔底为固图 2加劲梁截面几何尺寸单位：mm1 80018 0001 8009001 4653 70014 2003 7009001 465结；主缆在塔顶主鞍中心处按永不脱离点考虑；加劲梁梁端约束竖向、横桥向的位移和绕 轴（顺桥向）的转动；约束索塔下加劲梁横桥向的位移。静力三分力平均风对桥梁作用力称为静力三分力，按照体轴坐标，体轴系分为阻力系数、升力系数、扭矩系数 三个分量，按照风轴坐标，风轴系分为阻力系数、升力系数、扭矩系数 三个分量。风轴坐标下的静力三分力系数可表示为：阻力系数：（）升力系数：（）扭矩
5、系数：（）其中，为断面来流平均风速；为空气密度；，分别为桥梁断面的高度和宽度。风轴坐标系和体轴坐标系下的三分力系数存在如下换算关系，转换式如式（）所示，示意图如图 所示。（）第 卷 第 期 年 月 山西建筑 图 3体轴和风轴坐标下的三分力系数CLCVCDCMCHV 计算模型采用有限元软件建立二维流场，模型尺寸大小与实际相同，为了提高计算效率，对路灯、栏杆等进行忽略处理。计算区域为 ，流场几何尺寸见图。图 4二维流场几何尺寸15b36b8b8b 网格划分将流场分为 个区域进行网格划分，对桥梁截面周围进行网格加密，网格尺寸取 ，同时也为了平衡计算精度和计算速率，离截面较远的流场区域，受到截面的干扰
6、小，网格尺寸取 ，在上下两端的边缘区域尺寸取 。攻角下流场的网格划分如图 所示。图 5二维流场网格划分 边界条件流场的湍流模型采用大涡模拟湍流模型，桥梁截面壁面函数选取无滑移壁面函数，流场上下端选取滑移壁面函数，线入口风速设置 ，线出口压强为，对应计算模型的雷诺数为，湍流动能 ，湍流长度尺寸为，湍流黏度比为 ，其他参数均采用默认设置。流场材料选用空气，定义密度为 ，黏度为 （）。亚松弛因子等计算参数均采用默认值，时间增量 ，步骤数 。流场边界条件设置如图 所示。图 6流场边界条件设置无滑移壁面滑移壁面线入口V=12 m/s线出口P=0 Pa 静力三分力计算结果静力三分力系数能够反映桥梁在静风作
7、用下三个分量的受力与截面的关系，体现截面在风荷载作用下的气动性能。模拟主梁断面 的三分力系数汇总图见图（），图（）。图 7CFD 模拟主梁断面-12+12的三分力系数汇总图3.02.52.01.51.00.50.0-0.5-1.0-1.5-2.0三分力系数-12-9-6-3036912攻角/（）阻力系数 CD升力系数 CL扭转系数 CM3.02.52.01.51.00.50.0-0.5-1.0-1.5-2.0-2.5三分力系数-12-9-6-3036912攻角/（）阻力系数 CH升力系数 CV扭转系数 CM（a）风轴坐标系（b）体轴坐标系根据风轴和体轴的坐标系图，看到钢箱梁的静力三分力系数数值
8、变化较为平滑，没有出现数值骤变情况。主梁设计阻力系数在 攻角时最小，说明此时主梁受到的风荷载阻力最小，角度越大受到的阻力越大；升力系数随攻角增大逐步上升，风荷载对桥梁的升力作用由向下到向上；扭转系数增大幅度并不明显，可以推测主梁在风荷载作用下扭转可能性不大。桥梁风致响应分析 桥梁结构动力特性为了保证桥梁结构的稳定性与安全性，我们需要确定结构的动力反应分析，即结构频率、自振周期等，为分析结构性能提供依据。采用多重 向量法进行结构动力特性求解，表 为桥梁前五阶基频、周期、振型及振型特点。表 结构动力特性阶数频率 周期 振动特性 主梁纵飘 一阶反对称竖弯 一阶对称竖弯 一阶对称横弯 二阶对称竖弯 驰
9、振稳定性分析驰振是桥梁在运动中不断吸收能量，产生负阻尼使加劲梁振动逐渐增大，发生空气动力失稳的一种现象。根据 驰振理论，桥梁阻尼项的正负决定结构振动是否稳定，当气动力引起的负阻尼大于结构自身阻第 卷 第 期 年 月 山西建筑 尼才会导致结构的驰振发散。本文模型采用混凝土桥塔，阻尼较大，只需要对钢箱梁进行验证。公式为：（）（）（）其中，为结构断面驰振力系数，用于判断驰振稳定性。代入桥梁断面静力三分力系数，可知宽高比大于，驰振力系数大于，力的变化与速度方向相同，不会产生负阻尼，桥梁驰振稳定性良好，与桥梁结构动力特性扭转振型靠后情况相符。颤振稳定性分析颤振稳定性 反映了颤振稳定性能要求与桥梁颤振稳定
10、性的关系，可以通过该系数判别桥梁颤振综合性能，按下式计算：（）（）其中，因闭口箱梁取；为桥梁结构与空气密度比，经计算 ；为主梁梁扭转基频，取主梁的扭转基频；为主梁特征长度，为 。计算得，因为 ，根据规范桥梁的颤振临界风速 按下式计算：（）（）其中，为形状系数，取 ；为攻角效应系数，取；为与主梁相同宽度的理想平板颤振临界风速；为桥梁惯性半径，为 ；。主梁颤振检验风速为：（）其中，为考虑风洞试验误差及设计、施工中不确定因素的综合安全系数，取 ；为考虑风速的脉动影响及水平相关特性的无量纲修正系数，取 。，所以本桥发生颤振的可能性很小，满足颤振稳定性要求。涡振稳定性分析涡激振动发生时，结构的振动将对漩
11、涡的脱落产生反馈作用，使涡脱频率在相当长的风速范围内被结构自振频率所俘获，这种现象被称为“锁定”，这段风速范围被称为“锁定风速区间”，在锁定风速区间内，涡振振幅先增加后减小，存在一个极大值（见图）。图 8风速区间涡振振幅振幅0U3U2UU1先计算加劲梁断面的斯特罗哈尔数，由下式计算：（）其中，为涡流脱落频率，由升力系数卓越频率经傅里叶变换得到；为桥梁断面的高度；为来流平均风速。再计算涡激共振发生风速：（）其中，为弯曲或扭转振动频率，取动力特性竖弯振型频率。将升力系数卓越频率的傅里叶变换，在 风速下涡流脱落频率为，斯特劳哈尔数为（见图）。图 9涡流脱落频率0.350.300.250.200.15
12、0.100.050幅值102030405060708090100频率/Hz（3.52，0.346 2）采用 软件对桥梁截面进行涡振分析，在 攻角范围内出现涡旋脱落时，最大振幅为，对桥梁整体安全性影响不大，同时由于桥梁跨度不大且风作用对桥梁影响较小，可以认为该桥在自然风作用下不会发生较大的涡振振幅（见图）。图 100攻角涡振流场图根据静力三分力系数对钢箱梁进行驰振、颤振与涡振分析，判断出本桥主梁风致响应情况较好，不会出现明显的风致振动。桥梁结构静力分析 桥梁风环境 和 两个风作用水平对应的风速分别为重现期 与重现期 对应的风速值。风作用水平体现了频遇荷载的概念，当风荷载与汽车荷载及相关荷载作用组
13、合时，风荷载按 风作用水平确定；在 风作用水平下进行相关极限状态设计时，汽车荷载不 第 卷 第 期 年 月 文海泉等：钢箱梁悬索桥风致响应与静力分析参与荷载组合。在 上查得该地区地面以上 高，频率分别为 和 平均年最大风速，按最大风速服从极值型分布，可得 的 平均年最大风速期望值为 ，最大风速期望值 。桥梁位于入海口，地表粗糙度类型按类场地计算。主梁构件高度高出水面 ，所以桥梁设计基准风速可按下式计算：（）（）其中，为抗风风险系数，取值 ；为地形条件系数，取值 ；为地表类别转换及风速高度修正系数，取值 。根据抗风设计规范，当 风作用在主梁上的风速值大于 时，取 ；风作用风速取 。风作用下等效静
14、阵风荷载如今在工程实例中一般将风荷载等效成静力作用进行分析，本文将风对加劲梁和桥塔的作用以梁单元荷载施加到结构上，对主缆和吊索的作用以节点荷载施加到单元两端。风作用下主梁等效静阵风风速 可表示为：（）其中，为等效静阵风系数，按照规范取值。加劲梁上单位长度等效静阵风荷载按下式计算：（）（）（）其中，主梁的静气动力系数，选用 范围内的最不利值；为空气密度；，分别为桥梁断面的高度和宽度。作用桥塔、主缆、吊索的横桥向等效静阵风荷载只需要计算风荷载引起的阻力，按下列公式计算：（）（）其中，桥塔的风速按 塔高处风速计算，为等效静阵风系数；为各构件阻力系数；为构件单位上顺风向的投影面积，都按规范取值。风作用
15、下等效静阵风荷载、温度荷载、汽车与人群荷载 风作用下主梁等效静阵风风速 可表示为：（）作用在加劲梁、桥塔、主缆、吊索的风荷载都按式（）式（）计算。风灾是自然灾害中发生最频繁的灾害之一，但桥梁在正常使用阶段，不仅仅要考虑风荷载的影响，同时温度荷载与车道荷载也是常见荷载，需要考虑不同荷载组合的工况。温度荷载按体系升温（钢结构升温 ，混凝土索塔升温 ）；体系降温（钢结构降温 ，混凝土索塔降温 ）。汽车荷载采用公路级车道荷载，按照规范换算成均布荷载，横向 个车道，横向车道布载系数为 ，纵向折减系数为 ，普通车换算集度为 ，重车换算集度为 。当 风荷载与汽车荷载组合时，作用在汽车上的横向风荷载增加值取为
16、 。人群荷载集度为 （两侧）。计算结果分析为了分析不同的荷载工况对桥梁的影响，将荷载工况分为以下三种：工况：升温荷载；工况二：升温荷载 风荷载；工况三：升温荷载 风荷载 车道荷载。通过分析表 表 的数值模拟结果得出结论。表 加劲梁与塔顶位移荷载工况主缆横向位移梁横向位移梁竖向位移塔顶位移 表 主缆内力及应力表项目内力 应力 最大值 最小值 最大值 最小值 表 桥塔内力表（轴力受压为正）项目内力 弯矩 （）弯矩 （）最大值 最小值 最大值 最小值 表 加劲梁内力表（轴力受压为正）项目内力 扭矩（）最大值 最小值最大值 最小值 ）荷载工况中升温作用下加劲梁最大纵向位移为 ，最大竖向位移向下且为 ，如图（）所示，温度对钢箱梁悬索桥变形的影响需要重点考虑。）荷载工况中，风荷载作用桥梁产生最大横第 卷 第 期 年 月 山西建筑 向位移在跨中 ，如图（）所示，横向位移挠跨比 满足规范要求最大横向位移挠跨比。）荷载工况中，荷载组合引起的最大竖向位移为 ，如图（）所示，竖向位移挠跨比为，大于规范要求的。图 11三种荷载工产品销售下的位移图31.8928.9926.0923.1920.2917.3914
