非均匀风下流线型箱梁悬索桥静风稳定性分析

单秋伟 张亮亮 言志超 杨建新

摘 要：大跨度悬索桥结构轻柔，风致响应明显。借助有限元方法，考虑悬索桥的几何非线性和位移荷载非线性，对非均风攻角来流沿主梁对称分布和非对称分布以及非均匀风速来流沿主梁对称分布和非对称分布时悬索桥的非线性静风稳定性展开研究。结果表明：正攻角来流会降低桥梁的静风失稳临界风速，负攻角来流有利于桥梁抵抗静风失稳，且负攻角来流对桥梁的静风稳定性影响程度比正攻角来流影响程度大;非均匀风攻角来流非对称分布时，可由小风攻角来流确定的静风失稳临界风速作为稳定性判断依据，非均匀风攻角来流对称分布时，可以选择平均攻角来流衡量桥梁的静风稳定性。非均匀风速来流对桥梁的静风稳定性有不利影响，且非均匀风速来流对称分布时的影响比非均匀风速来流非对称分布时的影响更大。在不同初始攻角下，这类来流对桥梁静风稳定性有着相似的影响，桥梁的静风失稳临界风速均随着风速非均匀程度增大而减小。

关键词：悬索桥;流线型箱梁;非线性;靜风稳定性;非均匀风

中图分类号：U448.25   文献标志码：A   文章编号：2096-6717（2022）01-0117-09

收稿日期：2020-07-08

基金项目：重庆市教委科学技术研究项目（KJZD-k201802501）;重庆市技术创新与应用示范（社会民生类）一般项目（cstc2018jscx-msybX0277）

作者简介：单秋伟（1994- ），男，主要从事大跨度桥梁抗风研究，E-mail：shanqiuwei@foxmail.com。

张亮亮（通信作者），男，教授，博士生导师，E-mail：zll200510@126.com。

Abstract： The long-span suspension bridges， which are notable flexibility， are sensitive to wind load. Using the finite element method， considering the geometric nonlinearity and displacement-load nonlinearity， the aerostatic stability of suspension bridges is studied when the non-uniform wind attack angle flow and non-uniform wind speed flow distributes symmetrically and asymmetrically along the main girder. The results show that the flow of positive angle of attack will reduce the critical stability wind speed of the bridges， the flow of negative angle of attack is conducive to the bridges to resist the instability of static wind. The impact of negative angle of attack on the static wind stability of the bridge is greater than that of positive angle of attack. When the flow is asymmetrically distributed， the critical stability wind speed can be determined by small wind angle of attack. When the flow is symmetrically distributed， the mean angle of attack can be selected to measure the aerostatic stability of the bridge. Incoming flow with non-uniform wind speed has an adverse effect on the aerostatic stability of the bridge. Moreover， the influence of symmetric distribution of non-uniform wind speed is greater than that of asymmetric distribution of non-uniform wind speed.Under different initial angles of attack， this kind of incoming flow has similar effects on the static wind stability of the bridge. That's to say， the critical wind speed of aerostatic stability of the bridge decreases with the increase of the unevenness of the wind speed.

Keywords： suspension bridge; streamlined box girder; nonlinearity; aerostatic stability; non-uniform wind

随着桥梁建设能力的提升，桥梁跨径越来越大，但跨径的增大往往会使桥梁趋于轻柔，进而降低桥梁的稳定性。桥梁在建造和使用过程中一般需要考虑风灾、撞击和震害等影响，其中，风灾对桥梁结构的影响最为剧烈，发生频率也最高，由此造成的社会经济损失最为严重[1]。因此，在修建大跨度桥梁时，必须要考虑到桥梁结构对风作用的响应，并采取相应的措施降低其对结构的影响。

许多学者对大跨度桥梁的静风稳定性开展了研究。Boonyapinyo等[2]采用有限元方法，考虑桥梁几何非线性以及位移相关风荷载非线性，对大跨度斜拉桥在风荷载作用下的弯扭屈曲失稳临界风速进行了计算。程进等[3]采用非线性方法开展了江阴长江大桥的参数分析和比较，发现悬索桥的静风稳定性随着初始攻角的增大而降低。Cheng等[4-6]提出了大跨度斜拉桥和悬索桥静风稳定性的非线性方法，并开发了相应的分析程序。Boonyapinyo等[7]提出了一种考虑位移非线性、几何非线性和材料非线性的大跨度悬索桥静风稳定性的非线性分析方法，并分析了影响大跨度悬索桥静风稳定性的因素。李永乐等[8]采用风荷载增量和双重迭代相结合的方法对大跨度悬索桥的静风失稳形态以及机理进行了研究。Xu等[9]研究了斜拉桥的非线性静风稳定性，发现初始风攻角对桥梁静风稳定性影响极大。管青海等[10-11]分析了大跨度悬索桥的静风失稳发展过程以及影响因素。张文明等[12-13]研究了平均风速空间分布下悬索桥的静风稳定性，结果表明，风速空间分布对悬索桥静风稳定性的影响不可忽略。Zhou等[14]研究了双主跨悬索桥的静风失稳模式和破坏机理，发现此类悬索桥存在两种失稳模式。张玉琢等[15]在进行静风稳定性分析时发现非线性效应对静风稳定性的影响不可忽视。Dong等[16]基于逆可靠度评估法分析了大跨度斜拉桥的静风稳定性安全系数。胡朋等[17-18]对山区峡谷非均匀风场下大跨度斜拉桥的静风稳定性进行了分析，结果表明，非均匀风攻角下大桥的静风失稳临界风速远远低于均匀风攻角下的静风失稳临界风速。

上述研究多集中在空间均匀风下的静风稳定性，对非均匀风下的静风稳定性研究较少，因此，有必要对大跨度悬索桥在非均匀风下的静风稳定性开展研究，从而保证大跨度悬索桥在建设和运营状态下的安全。笔者利用有限元分析软件ANSYS的APDL语言，采用非线性分析方法，对非均匀风攻角和非均匀风速下大跨度悬索桥的静风稳定性进行分析。

1 静风稳定性分析理论

1.1 风的静力作用

平均风作用会使处在风场中的结构产生一定的变形，风的作用相当于一个静荷载[1]。将此荷载沿桥梁主梁断面分解为阻力、升力和升力矩，称为风荷载的三分力，如图1所示。

借助ANSYS对大跨度悬索桥进行非线性静风稳定性分析的具体过程：

1）给出初始风攻角α0和初始风速U0以及风速步长ΔU;

2）根据此时的风攻角和风速，计算主梁受到风荷载的三分力;

3）在ANSYS的静态分析中，采用完全Newton-Rapson法求解式（5），得出模型主梁的位移，并从中提取出主梁各单元的扭转角，重新计算主梁受到风荷载的三分力;

4）计算三分力系数的欧几里得范数，当范数大于给定值时，重复执行过程3），使主梁达到该风速下的平衡状态，若范数小于给定值，则增加一级风速，重复执行过程2）和3）;

5）若在计算过程中，出現了迭代不收敛，则将风速调整到上一次迭代状态，减小风速步长，再进行计算，相邻两次迭代风速差值满足要求时便可终止计算，相对应的最后一级风速即是失稳临界风速。

2 大跨度悬索桥有限元模型

某大跨度悬索桥，全长为1 600 m，由主桥和引桥构成，主桥跨径为880 m。南引桥共两联，总长分别为230、240 m，其中，第1联为连续梁结构，第2联为连续刚构结构;北引桥长度为240 m，采用连续刚构结构。中跨主缆的矢跨比为1/8.8，桥梁的立面整体布置如图3所示。

在ANSYS中，采用单主梁模型模拟主梁，单元类型选用beam4;主塔因其为变截面，故采用3D渐变梁单元beam44模拟;采用link10杆单元模拟主缆和吊索;mass21质量单元模拟主梁的质量惯矩和二期恒载。

大桥的边界条件设置为：主塔底部固结，主缆锚固处固结，桥塔下横梁与主梁采用耦合自由度方式进行连接，主塔下横梁处设置combin14弹簧单元模拟纵向阻尼器，以限制主梁纵向位移。

全桥有限元模型共计705个节点，887个单元，如图5所示。

3 非均匀风攻角下的静风稳定性

自然界中的风场都是非均匀的，而前人的研究多以均匀流场为基础，这样的研究结果与真实情况可能存在偏差，因此，有必要开展非均匀风场下悬索桥的静风稳定性分析。

为了研究非均匀风场中的非均匀风攻角风场对流线型箱梁悬索桥静风稳定性的影响，根据风攻角分布区间的不同，设置了如图6所示的两类风攻角分布工况。

工况1考虑风攻角沿着桥跨方向非对称分布的情况，主梁右半跨来流风攻角为α1，左半跨来流风攻角为α2。工况2则考虑风攻角沿着桥跨方向对称分布的情况，主梁跨中1/2跨度范围来流风攻角为α1，主梁两端各1/4跨度范围来流风攻角为α2。

由图12（a）可知，上述4种衡量方式所对应的偏离幅度a、b、c、d中，衡量方式Ⅱ的偏离幅度b最接近0，方式Ⅰ的偏离程度a最大，说明此衡量方式偏离实际情况的程度最严重，方式Ⅲ的偏离幅度c和方式Ⅳ的偏离幅度d也较大。说明对非对称风攻角工况采用较小风攻角计算的静风失稳临界风速最接近实际情况，因此，可以采用小风攻角对非均匀风攻角来流进行静风稳定性分析。但需要注意的是，当半跨来流风攻角接近于0°，另半跨来流负风攻角较大时，方式Ⅱ的偏离幅度较大，说明偏离程度较远，此时，需要按照方式Ⅲ或方式Ⅳ来衡量静风稳定性。

由图12（b）可知，工况2下衡量方式Ⅳ的偏离幅度d在整体范围内最接近0，表明对于风攻角对称分布工况，可以采用大小风攻角的平均值进行非均匀风攻角来流下桥梁静风稳定性的分析。但当主梁端部攻角接近于0°，跨中附近负攻角较大时，方式Ⅳ的偏离幅度较大，此时，应该采用小风攻角进行静风稳定性分析。

4 非均匀风速下静风稳定性

作用在大跨度悬索桥上的风，除了风攻角的分布存在不均匀的情况外，风速大小也可能沿着桥跨方向出现不均匀的情况。当风速分布不均匀时，为了便于分析，设置了两类风速非均匀分布工况，如图13所示。

参考明渠流流层流速分布规律，采用指数规律来描述主梁上的风速变化。工况1考虑风速沿主梁非对称分布情况：主梁一端风速为U1，另一端风速为U2，主梁上的风速由U1到U2呈指数律变化。工况2考虑风速沿主梁对称分布情况：主梁端部风速为U1，跨中风速为U2，端部到跨中的风速由U1到U2呈指数律变化。

为了分析风速的非均匀程度对大跨度悬索桥静风稳定性的影响，设置了如表2所示的风速比值关系。

4.1 非均匀风速下静风失稳临界风速

非均匀分布风速下大桥静风稳定性的比较须在同一基准下进行。以主梁全长范围内非均匀分布风速下单位时间内总流量与均匀分布风速下单位时间内总流量相等为原则，在主梁全长范围内对风速积分后取平均值作为平均静风失稳临界风速Ucr，即式

由图14可以发现，在两类风速分布工况下，大桥的平均静风失稳临界风速随着风速比值U1/U2的降低而降低，表明风速非均匀分布不利于流线型箱梁悬索桥抵抗静风失稳，会降低其静风稳定性。此外，风速对称分布工况下，大桥的平均静风失稳临界风速降低幅度比非对称分布工况下的降低幅度大，说明风速对称分布工况更加不利于该类型悬索桥抵抗静风失稳。

4.2 不同风攻角下非均匀风速来流对静风稳定性的影响

上述分析仅考虑了给定攻角下的非均匀风速，为了说明不同初始风攻角下非均匀风速对静风稳定性的影响，以工况2为例，分析了-2°、0°、+2°、+4°、+8°攻角下的非均匀风速来流对桥梁静风稳定性的影响。其中较小的攻角主要是考虑大桥所在桥址以及周围环境，较大的攻角则是考虑一些极端气候因素以及局地强风影响。不同风攻角下不同风速比值的平均静风失稳临界风速如图15所示。

由图15可知，不同初始攻角下，非均匀风速来流对桥梁静风稳定性有着相似的影响，在-2°、0°、+2°、+4°、+8°攻角下，悬索桥的平均静风失稳临界风速均随着风速比值减小而降低，这与前述分析相吻合，也表明了风速非均匀分布不利于流线型箱梁悬索桥抵抗静风失稳。

5 结论

以某一流线型箱梁悬索桥为背景，采用有限元方法分析了悬索桥在非均匀风作用下的非线性静风稳定性，得到以下结论：

1）当非均匀风攻角来流作用时，随着正攻角的增大，桥梁的静风失稳临界风速逐渐降低，正攻角会降低桥梁的静风稳定性。随着负攻角的增大，桥梁的静风失稳临界风速逐步提高，负攻角有利于桥梁的静风稳定性。通过对比正负攻角的影响程度，发现负攻角对桥梁静风稳定性的影响程度比正攻角的影响程度大。

2）当非均匀风攻角来流非对称分布时，由小风攻角确定的桥梁静风失稳临界风速可以在一定程度上代表桥梁实际来流情况下的静风失稳临界風速;当非均匀风攻角来流对称分布时，可以选择平均攻角确定的静风失稳临界风速来衡量桥梁的静风稳定性。

3）非均匀风速来流会降低桥梁的静风失稳临界风速，对桥梁的静风稳定性有不利影响，并且来流风速对称分布时对该类型桥梁静风稳定性的影响比来流风速非对称分布时的影响更大。不同初始攻角下，非均匀风速来流对桥梁静风稳定性有着相似的影响，悬索桥的平均静风失稳临界风速均随着风速非均匀程度增大而降低。参考文献：

[1] 陈政清. 桥梁风工程[M]. 北京： 人民交通出版社， 2005.

CHEN Z Q. Bridge wind engineering [M]. Beijing： China Communications Press， 2005. （in Chinese）

[2]  BOONYAPINYO V， YAMADA H， MIYATA T. Wind-induced nonlinear lateral-torsional buckling of cable-stayed bridges [J]. Journal of Structural Engineering， 1994， 120（2）： 486-506.

[3] 程进， 肖汝诚， 项海帆. 大跨径悬索桥静风稳定性的参数研究[J]. 公路交通科技， 2001， 18（2）： 29-32.

CHENG J， XIAO R C， XIANG H F. Parametric studies on aerostatics stability for suspension bridges [J]. Journal of Highway and Transportation Reseach and Development， 2001， 18（2）： 29-32. （in Chinese）

[4]  CHENG J， JIANG J J， XIAO R C， et al. Advanced aerostatic stability analysis of cable-stayed bridges using finite-element method [J]. Computers & Structures， 2002， 80（13）： 1145-1158.

[5]  CHENG J， JIANG J J， XIAO R C. Aerostatic stability analysis of suspension bridges under parametric uncertainty [J]. Engineering Structures， 2003， 25（13）： 1675-1684.

[6]  CHENG J， XIAO R C， XIANG H F， et al. NASAB： A finite element software for the nonlinear aerostatic stability analysis of cable-supported bridges [J]. Advances in Engineering Software， 2003， 34（5）： 287-296.

[7]  BOONYAPINYO V， LAUHATANON Y， LUKKUNAPRASIT P. Nonlinear aerostatic stability analysis of suspension bridges [J]. Engineering Structures， 2006， 28（5）： 793-803.

[8] 李永樂， 欧阳韦， 郝超， 等. 大跨度悬索桥静风失稳形态及机理研究[J]. 空气动力学学报， 2009， 27（6）： 701-706.

LI Y L， OUYANG W， HAO C， et al. Study on shape and mechanism of aerostatic stability for long span suspension bridges [J]. Acta Aerodynamica Sinica， 2009， 27（6）： 701-706. （in Chinese）

[9]  XU M， GUO W W， XIA H， et al. Nonlinear aerostatic stability analysis of Hutong cable-stayed rail-cum-road bridge [J]. Wind and Structures， 2016， 23（6）： 485-503.

[10] 管青海， 张凯， 李加武， 等. 大跨人行悬索桥非线性静风失稳发展过程分析[J]. 桥梁建设， 2018， 48（1）： 76-81.

GUAN Q H， ZHANG K， LI J W， et al. Analysis of development process of nonlinear aerostatic instability of a long span pedestrian suspension bridge [J]. Bridge Construction， 2018， 48（1）： 76-81. （in Chinese）

[11] 管青海， 周燕， 李加武， 等. 主跨420 m人行悬索桥非线性静风稳定影响参数分析[J]. 振动与冲击， 2018， 37（9）： 155-160.

GUAN Q H， ZHOU Y， LI J W， et al. Effects of parameters on nonlinear aerostatic stability of a pedestrian suspension bridge with main span of 420 m [J]. Journal of Vibration and Shock， 2018， 37（9）： 155-160. （in Chinese）

[12] ZHANG W M， QIAN K R， WANG L， et al. Aerostatic instability mode analysis of three-tower suspension bridges via strain energy and dynamic characteristics [J]. Wind and Structures， 2019， 3（29）： 163-175.

[13] 张文明， 葛耀君. 考虑风速空间分布的三塔悬索桥静风稳定分析[J]. 东南大学学报（自然科学版）， 2018， 48（1）： 152-156.

ZHANG W M， GE Y J. Aerostatic stability analysis of suspension bridge with three towers considering spatial distribution of wind speed [J]. Journal of Southeast University （Natural Science Edition）， 2018， 48（1）： 152-156. （in Chinese）

[14]  ZHOU Q， LIAO H L， WANG T. Numerical study on aerostatic instability modes of the double-main-span suspension bridge [J]. Shock and Vibration， 2018， 2018： 7458521-7458529.

[15] 张玉琢， 刘海卿， 马凯. 大跨径悬索桥静风稳定性研究[J]. 辽宁工程技术大学学报（自然科学版）， 2019， 38（3）： 258-262.

ZHANG Y Z， LIU H Q， MA K. Research on aerostatic stability of long span bridge [J]. Journal of Liaoning Technical University （Natural Science）， 2019， 38（3）： 258-262. （in Chinese）

[16]  DONG F H， CHENG J. A new method for estimation of aerostatic stability safety factors of cable-stayed bridges [J]. Proceedings of the Institution of Civil Engineers-Structures and Buildings， 2019， 172（1）： 17-29.

[17] 胡朋， 颜鸿仁， 韩艳， 等. 山区峡谷非均匀风场下大跨度斜拉桥静风稳定性分析[J]. 中国公路学报， 2019， 32（10）： 158-168.

HU P， YAN H R， HAN Y， et al. Aerostatic stability of long-span cable-stayed bridge under inhomogeneous wind fields induced by mountain-gorge terrain [J]. China Journal of Highway and Transport， 2019， 32（10）： 158-168. （in Chinese）

[18]  HU P， HAN Y， XU G J， et al. Effects of inhomogeneous wind fields on the aerostatic stability of a long-span cable-stayed bridge located in a mountain-gorge terrain [J]. Journal of Aerospace Engineering， 2020， 33（3）： 04020006.

[19] 李文杰. 大跨径悬索桥静风稳定性分析[D]. 西安： 长安大学， 2009.

LI W J. Study on aerostatic stability of long-span suspension bridges [D]. Xian： Changan University， 2009. （in Chinese）

[20] 张志田， 葛耀君. 基于正交异性壳单元的悬索桥非线性静风稳定性分析[J]. 中国公路学报， 2004， 17（4）： 64-69.

ZHANG Z T， GE Y J. Nonlinear aerostatic instability analysis of suspension bridges based on orthotropic shell element models [J]. China Journal of Highway and Transport， 2004， 17（4）： 64-69. （in Chinese）

[21] 戴礼勇. 大跨度公铁两用悬索桥静风稳定性研究[D]. 成都： 西南交通大学， 2013.

DAI L Y. Aerostatic stability of long-span combined highway and railway suspension bridge [D]. Chengdu： Southwest Jiaotong University， 2013. （in Chinese）

（編辑 胡玲）