静风荷载引起的超大跨度桥梁关键问题研究

V ol118　N o16公　路　交　通　科　技2001年12月JOURNA L OF HIGHWAY AND TRANSPORT ATION RESEARCH AND DEVE LOPMENT 文章编号:1002Ο0268(2001)06Ο0034Ο05
收稿日期:2000Ο11Ο27静风荷载引起的超大跨度桥梁关键问题研究
肖汝诚,贾丽君,程　进,宋　馨,孙　斌,项海帆
(同济大学土木工程学院,上海　200092)
摘要:在超大跨径桥梁中,静风荷载不仅会引起结构的动力特性改变,还将导致结构强度破坏和失稳。

本文以超大跨径桥梁为研究对象,计入几何、材料以及静风荷载的非线性的三重影响,对它们在静风荷载作用下的关键问题进行了研究,揭示了结构的静风荷载响应与三分力系数曲线关系,给出了研究方法和一系列研究成果。

关键词:超大跨径桥梁;静风荷载;非线性
中图分类号:U44112 文献标识码:A
Re search on the K ey Points Caused by Aero static Load in SuperΟlong Bridge s
XIAO RuΟcheng,JIA LiΟjun,CHENG Jin,SONG Xin,SUN Bin,XIANG HaiΟfan
(Institute of Civil Engineering,T ongji Univ1,Shanghai　200092,China)
Abstract:Aerostatic load will not only cause alteration of dynamic character of super-long-span bridge structures,but als o will lead the structure to collapse or buckling1T aking the multiple in fluence of the geometrical nonlinear and the material nonlinear and the nonlin2 ear caused by aerostatic load into consideration,the key points of the structure under aerotatic load are researched in this paper1The re2 lationship between the response of aerostatic load and three components of static aerostatic coefficients are als o discussed1A series of re2 search methods and results are presented1
K ey words:Super-long-span bridge;Aerostatic load;N onlinear
0　引言
国内外跨海工程的建设离不开建造超大跨径桥梁,但是,桥梁跨径增大,势必会带来一系列新问题,静风荷载问题就是其例。

在动力特性方面,现有理论计算大跨径桥梁的动力特性,一般忽略随时间变化的动力荷载非线性影响,但是,风速变化时大跨度桥梁的几何变形与内力状态都将发生变化,结构的几何刚度和质量矩阵也随之变化,从而可能影响到结构的动力特性。

在强度方面,过去人们普遍认为大跨径桥梁的强度主要是受恒活载或地震荷载控制的,但我们在对香港青龙大桥(主跨1418m悬索桥)进行设计复核时发现,主塔构件的强度是由静风荷载控制的。

稳定方面,人们传统认为大跨径桥梁颤振临界风速一般都低于其静风失稳临界风速,但是,1967年日本东京大学Hirai教授在悬索桥的全桥模型风洞试验中观察到了静力扭转发散的现象,同济大学风洞实验室在对汕头海湾二桥的风洞试验中,发现了斜拉桥由静风引起的弯扭失稳现象[1]。

因此出现了静风荷载引起的超大跨度桥梁动力特性、强度与稳定的新问题,这些都是超大跨径桥梁在静风荷载作用下的关键问题。

本文以超大跨径桥梁为研究对象,计入几何、材料以及静风荷载的三重非线性影响,对悬索桥和斜拉
桥在静风荷载作用下的关键问题进行了研究。

给出了研究方法和研究成果。

1　静风作用下的结构计算理论
111　静风荷载的描述
静风荷载对大跨径桥梁的作用一般简化为风对结构的阻力、升力和升力矩的三分力的共同作用。

作用在主梁上的三分力(图1)表达式为
P H=1/2・ρV d2・C H0・D
P V=1/2・ρV d2・C V0・B(1)
P M=1/2・ρV d2・C M0・B2
式中,ρ为空气密度;D、B为主梁截面的高度和宽度;C H0、C V0、C M0为初始攻角时主梁沿体轴坐标系各方向的三分力系数。

风洞试验结果表明,三分力系数是风的有效攻角(图2,3)的函数。

大跨径桥梁是柔性结构,在静风作用下,结构的姿态将发生改变,导致静风与主梁截面间的有效攻角变化,其三分力也随有效攻角而改变。

这样,不仅风速自身的增长会引起静风荷载呈非线性变化,三分力系数的变化也会导致静风荷载的非线性变化。

因此,将式(1)用于超大跨径桥梁的静风响应分析,将无法获得结构的准确静风平衡点。

合理的分析方法应考虑三分力系数随有效攻角改变的影响。

有效攻角α为静风初始攻角θ0与静风作用引起的主梁扭转角θ之和。

静风荷载可表示为
P H=1/2・ρV d2・C H(α)・D
P V=1/2・ρV d2・C V(α)・B(2)
P M=1/2・ρV d2・C M(α)・B2
式中,C H(α)、C V(α)、C M(α)为随攻角变化的三分力函数,可通过节段模型风洞试验实测得;D、B为主梁截面的高度与宽度。

图1　
静力三分力示意图图2　
广东虎门桥三分力实测曲线图3　有效攻角示意图
112　静风荷载作用下的平衡方程
考虑静风荷载受有效攻角的影响,静风三分力引
起的等效节点力可以写成结构变形的函数。

大跨径桥
梁在静风荷载作用下的非线性有限元分析可归结为求
解以下的非线性平衡方程
[K(δ)]{δ}={F(α,v)}(3)
式中,[K(δ)]为大跨径桥梁的总体切线刚度矩阵;
{F(α,v)}为风速v和有效攻角α时的风载等效节
点力向量。

对式(3)采用U L增量法求解,相应非线性增
量平衡方程组如下
([K0]+[Kσ
j-1
(δj-1)])・{Δδj}={F j(αj,v i)}
-{F j-1(αj-1,v i)}(4)
式中,[K0]为大跨度桥梁的线弹性刚度矩阵;[Kσ
j-1
(δj-1)]为第j-1步状态时,单元的几何刚度矩阵;
{F j(αj,v i)}为i级风载j步有效攻角αj的风载等效节
点力向量;{F j-1(αj-1,v i)}为i级风载j-1步有效攻
角αj-1的风载等效节点力向量。

113　静风作用下的振动方程
为了求解桥梁结构的动力特性(包括频率和振
型),首先要建立结构的振动方程,在此方程中,结构的
几何刚度矩阵和质量矩阵随结构姿态和内力状态的变
化而改变,表现为结构位移{δ}的函数。

但是,只要风
速v给定,就可以根据方程(4)计算出结构在此风速下
的平衡位置,其结构刚度矩阵和质量矩阵也随之而定。

此时大跨度桥梁自由振动微分方程组为
[M(δv)]・{¨
δ}+([K0]+[Kσ(δv)])・{δ}=0(5)
式中,δv———风速为v时结构的节点静位移向量;
δ———结构的节点振动位移向量。

容易得到相应的频率方程为
[K(δv)]{x}-ω2v[M(δv)]{x}=0(6)
2　气静力稳定性
大跨径桥梁在静风荷载作用下,主梁发生弯曲和
扭转,一方面改变了结构刚度,另一方面改变了风荷
载的大小,并反过来增大结构的变形,最终导致结构
静风荷载引起的超大跨度桥梁关键问题研究　肖汝诚等
失稳的现象称为空气静力失稳简称静风失稳。

211　静风稳定性计算方法
考察式(3)可知,结构刚度和静风荷载都是结构变形的函数,为了求解该非线性方程,本文在综合考虑结构几何、材料非线性和静风荷载非线性的基础上提出了采用增量与内外两重迭代相结合的方法。

风速按一定比例增加的过程中,内层迭代完成结构的非线性计算,外层迭代寻找结构在某一风速下的平衡位置。

该方法的具体实施步骤如下:
(1)假定初始风速V 0、荷载参数λ及荷载参数增量Δ
λ。

(2)计算在风速V =V 0×λ下结构所受的静风荷载。

(3)采用Newton ΟRaps on 法求解(3)式,得到结构位移δ。

如果结构中单元出现塑性铰就进行总刚重组。

(4)从结构位移δ中提取单元扭转角(为左右两节点扭转位移之和的平均值),重新计算结构的静风荷载。

(5)检查三分力系数的欧几里得范数是否小于允许值,如下式所示
∑
N α
j =1
[C k (αj )-C k (αj -1)]
2
∑
N α
j =1
[C k
(αj -
1)]2
1/2
Φεk
(k =X ,Y ,Z )(7)
式中,N α为受到静风荷载作用的节点总数;C k 为阻
力、升力和升力矩系数;εk 为阻力、升力和升力矩系数的允许误差。

(6)若小于允许值,判断结构中是否有单元出现塑性铰,如有,记录出现塑性铰的单元号及相应状态,调整荷载参数λ,重新计算至该塑性铰处的弯矩值等于该断面处的极限弯矩值,重复步骤(2)～(5)步;如未出现塑性铰,令λ=λ+Δλ,重复步骤(2)～(5)进行计算。

(7)若大于允许值,则重复步骤(3)～(5)。

(8)若在某一级风速V 下出现迭代不收敛,恢
复到上一级风速状态,缩短步长,重新计算,直至相
邻两次风速之差小于预定值为止。

212　算例分析
根据现有静风稳定计算方法,分别计算了两座大跨度桥梁的静风稳定性。

(1)斜拉桥:本文以日本学者T 1Miyata 设计的
1000m 跨径的斜拉桥为结构模型(具体数据可参考
文献[5]),主梁截面取南京二桥形式进行静风稳定
计算。

计算结果如表1所示。

图4为各种非线性下主梁跨中处扭转变形随风速变化的比较。

各种非线性下斜拉桥的静风临界风速比较表1
几何非线性
几何非线性+静风荷载非线性
所有非线性影响
临界风速(m/s )
177
129
114
图4　各种非线性下主梁跨中处扭转变形随风速变化的比较
(2)悬索桥:以虎门桥为例,采用两种不同的分
析方法对其空气静力稳定性进行了计算。

虎门大桥是中跨88810m 的钢箱梁悬索桥,主梁梁高31012m ,桥宽3516m ,桥塔为门式框架结构,塔高为150m ,主缆间距为33m ,吊杆共2×72对,间距为1210m ,其主要构件截面材料及几何特性见表2。

计算结果列于表3。

图5为主梁跨中点处扭转变形随风速的变化历程。

截面材料及几何特性
表2
类型J d (m 4)I 2(m 4)I 3(m 4)m (t/m )E (MPa )
γ
钢箱梁5109551241392119786181336211E05013主　缆---21397210E05
-吊　杆
-
-
-
01172
116E05
-
不同方法的临界风速比较
表3
计算方法方法1
方法2
临界风速(m/s )
136
119
注:方法1是线性静力扭转发散公式计算出的结果;方法2为本文方法计算出的结果。

图5　主梁扭转变形变化过程(方法4)
计算表明:
(1)一般情况下,采用线性方法计算出的大跨度
公路交通科技　2001年　第6期
桥梁失稳临界风速比非线性结果高。

(2)作用在结构上的静风荷载是非线性的,所以结构的变形随风速的变化呈明显的非线性。

(3)大跨径桥梁的空气静力失稳表现为空间弯扭耦合失稳。

(4)计入材料非线性计算静风临界风速较不计入的结果小,但失稳时结构并不变成机构,这是因为在一般情况下,材料非线性降低了结构切线刚度,但不
是引起静风失稳的主要原因。

213　参数研究与特殊现象分析
21311　参数研究
为了避免静风失稳,必须提高其临界风速。

因此有必要考察设计参数变化对结构静风失稳风速的影响,以获得改善大跨径桥梁空气静力稳定性的方法。

作者分别就各种参数对斜拉桥和悬索桥空气静力稳定性的影响进行了研究,限于篇幅,本文仅给出研究结果。

a1结构宽跨比增大,其空气静力稳定性提高,宽跨比增大25%,临界风速提高1712%。

b1增加桥面均布荷载可以提高桥梁的空气静力稳定性。

c1采用不同的主梁断面将明显改变斜拉桥的静风稳定性。

d1初始攻角增大,斜拉桥的静风稳定性有所下降。

e1增加斜拉桥主塔高度,结构的静风稳定性降低。

f1改变斜拉桥边跨跨径对结构的静风稳定性影响不大。

g1考虑斜拉索上的静风荷载将降低结构的静风稳定性。

h1斜拉索的垂度效应会明显降低结构的抗静风能力,仅采用Ernst公式计入拉索垂度效应是不够的,只有采用悬链线索单元考虑斜拉索垂度效应才能比较真实地反映其空气静力稳定性。

i1悬索桥主缆垂度效应对结构的静风稳定性影响不大,计算时可以不计缆索的垂度效应。

21312　特殊现象分析
作者在研究中观察到两个特殊现象。

一是在对主跨518m的汕头海湾二桥进行分析时,发现该桥在0
°攻角下的空气静力失稳风速(129m/s
)低于颤振临界风速(140m/s),即空气静力失稳先于动力失稳,该现象已在同济大学风洞试验室里得到证实。

但一般而言,发生这种现象的斜拉桥主跨跨径应在800m以上。

图6中曲线1为汕头海湾二桥主梁断面升力矩系数实测值,曲线2为常规主梁断面升力矩系数曲线。

用上述两种升力矩系数分别对汕头海湾二桥进行分析,结果如表4所示。

图7为主梁跨中点处扭转角随风速变化过程。

汕头海湾二桥空气静力失稳风速比较表4计算结果升力矩系数1升力矩系数2临界风速(m/s)129140
图6　升力矩系数曲线
图7　主梁跨中点处扭转角
计算结果表明,产生这种现象的主要原因是该桥主梁断面的升力矩曲线在攻角大于3°后的形状与常规断面不同。

另一个现象是江阴大桥静风临界风速线性结果(97m/s)比非线性结果(113m/s)低,考察江阴桥截面的升力矩系数曲线可知(图8所示),产生这种现象的主要原因是线性公式近似采用0°攻角下的升力矩系数曲线斜率作为临界风速时的斜率,此时该曲线的斜率最大。

而非线性分析方法考虑了斜率随攻角的变化,临界风速作用下升力矩系数曲线斜率比0°攻角下的斜率小,从而导致了上述结果的发生。

为了验证这一解释,分别采用江阴桥和虎门桥的三分力系数计算曲线,对江阴大桥的静风稳定性进行计算。

计算结果列于表5。

计算结果证明了我们的判断。

综合以上两种现象可以看出,大跨度桥梁的静风临界风速与结构主梁断面升力矩系数曲线的形状密切静风荷载引起的超大跨度桥梁关键问题研究　肖汝诚等
江阴长江大桥临界风速比较(单位:m/s)表5不同方法线性方法非线性方法
临界风速(江阴桥
三分力系数)
97113
临界风速(虎门桥
三分力系数)
10614393
相关,改善升力矩系数曲线形状可以有效改善大跨径桥梁空气静力稳定性。

3　静风荷载对动力特性的影响
根据式(6),按如下方法容易求出风速v时的n 个自振频率和振型{x}。

(1)根据施工方法,求出结构的恒载内力和构形以确定零风速下成桥初始状态,以此状态下的[Kσ]G形成结构总刚([K0]+[Kσ]G)。

(2)根据给定风速,增加一级风速v1=v0+Δv,计算三分力及其等效节点力{F(δ,v1)},通过NewtonΟRaphs on法与增量法计算方程(2),获得在此风速下的结构状态(位移,内力等),进而求得几何刚度矩阵[Kσ(δv
1
)]。

(3)重复第(2)步,直至v1=v,求出相应的[Kσ(δv)]和[M(δv)]。

(4)将计算得到的几何刚度矩阵[Kσ(δv)]和质量矩阵[M(δv)]代入式(6),得到在风速v下桥梁的自振频率[ωv]和自振的振型。

(5)输出结果。

为了研究静风荷载对大跨度桥梁动力特性的影响,作者以虎门大桥为例进行了静风荷载作用下动力特性的分析。

图9、图10给出了虎门大桥的动力特性随风速的变化曲线,当风速较小时,各阶频率随风速变化不大。

但在接近80m/s风速时,频率有增大的趋势。

这是因为此时风的升力作用方向向下,与升力矩共同作用,使得主缆总体索力增加,从而使频率增加。

而当风速超过80m/s后,升力反向向上作用,使结构丧失了一部分重力刚度,频率开始下降。

当风速临近120m/s时,频率急剧下降。

图8　
虎门大桥和江阴长江大桥升力矩系数图9　
各阶频率随风速变化过程图10　弯曲频率随风速变化过程
通过以上以及其它桥梁的分析比较,我们得出如
下结论:
(1)大跨度桥梁的动力特性与静风荷载有关,斜
拉桥与悬索桥相比,静风荷载对悬索桥的动力特性影
响更大。

(2)常风速风载对结构动力特性影响较小,可以
忽略不计。

(3)在接近静风失稳阶段,结构的弯、扭频率急
速下降,计算其动力特性时必须计入静风效应。

(4)在频域内分析大跨度颤振临界风速时,如果
其临界风速与静风临界风速接近时,必须考虑静风对
动力特性的影响。

4　小结
本文以超大跨径桥梁为研究对象,计入几何、材
料以及静风荷载的非线性的三重影响,对它们在静风
荷载作用下的关键问题进行了研究,揭示了结构的静
风响应、动力特性与结构设计参数和三分力系数之间
的内在联系,得出了以下结论:
(1)结构形式、主梁断面形式、风的初始攻角等
因素对大跨度桥梁的静风响应都有不同程度的影
响。

(2)大跨径桥梁静风失稳时的构形表现为空间弯
扭耦合失稳,扭转变形对结构静风响应的影响是明显
的。

(3)计入材料非线性,静风临界风速较不计入的
结果小,但失稳时结构并不变成机构,即材料非线
性降低了结构切线刚度,但非引起失稳的主要原
因。

(4)大跨度桥梁主梁断面的升力矩曲线斜率与其
静风临界风速关系密切,升力矩曲线梯度小,结构的
空气静力稳定性就好,改善主梁断面的升力矩曲线,
可以改善大跨径桥梁空气静力稳定性。

(5)常风速风载对结构动力特性(下转第54页)
控制方向皆为y 方向;控制点挠度δ=-010776-010775
,
δ0=
01080
01080
,则调索后理论上控制点挠度应为Δ=
δ+δ0=
010024
10025。

其调整结果见表2,调索过程的施工简图及调索后的挠度变化图见图2、图3。

图2　
图3　由表2可见:Δ′=
-010056-010130
,而Δ=
010024
010025。

比较Δ′与Δ,010155m 调索结果(单位:m )表2
调整次数　
索单元号控制点Ο控制方向
的挠度
河跨
岸跨
18Οy
20Οy 12322-010725-010********-010604-010********-010440-010********-010260-0103205
27
18
-010056
-010130
4　结语
11根据连续刚构桥和斜拉桥各自的结构特点、
受力形式、施工过程特点,运用工程控制论的思想,建立了适合两种桥型的数学模型、目标函数,对象明确,力学概念清晰,实现程序化计算方便。

21从数学的角度,采用统计的方法,对两种桥
型的状态变量、控制变量进行了较为准确的计算,为有效地实施施工控制提供了保证。

31在建立模型的基础上,编制了相应的计算程序,使计算机技术与最优控制方法有效地结合,实现了对连续刚构桥和斜拉桥施工的实时跟踪分析。

41通过对两种桥型施工控制方法的比较以及在两座大桥上的实际应用,证明了随机最优控制系统的实用性,同时也证明了随机最优控制系统理论在桥梁施工控制中值得进一步推广,具有广阔的应用前景。

参考文献:
[1]　[美]小沃尔特波多尔尼,[法]J M 米勒尔1预应力混凝土桥
梁分段施工和设计1人民交通出版社,19861
[2]　林元培1斜拉桥1人民交通出版社,19941
[3]　钱学森,宋健1工程控制论(修订本)(上、下册)1科学出版
社,19831
[4]　[瑞典]K J 奥斯特隆姆1随机控制理论导论1科学出版社,
19801
(上接第38页)影响较小,可以忽略不计,在接近静
风失稳阶段,结构的弯、扭频率急速下降,计算其动力特性时必须计入静风效应。

参考文献:
[1]　M Nagai.P ossibility and Limitations of Long Οspan Cable Οstayed Bridges
Based on S tatic and Dynamic Instability Analyses 1Japan ,IABSE Sym 2posium K obe 1998,IABSE Reports ,(V ol79,K obe )1998.
[2]　方明山,项海帆,肖汝诚1超大跨径悬索桥空气静力非线性行
为研究1重庆交通学院学报,1999,18(2).
[3]　李国豪1桥梁结构稳定与振动1中国铁道出版社,1996.[4]　肖汝诚,薛惠钰1桥梁结构线性非线性动力分析程序系统1华
东公路,1995(3).
[5]　程进,肖汝诚,项海帆1大跨径斜拉桥非线性静风稳定性全过
程分析1中国公路学报,2000(2).
[6]　程进,肖汝诚,项海帆1大跨径悬索桥静风扭转发散的级数解
法1同济大学学报,1999(2).。