桥梁抗风设计的目的

Download Report

Transcript 桥梁抗风设计的目的

第八章 桥梁抗风设计
本章主要介绍桥梁动力特性、作用于桥梁上的
风荷载的计算及桥梁动力失稳判断等内容。此处的
桥梁动力特性主要涉及桥梁的自振周期及频率,本
章介绍了如何用结构动力学方法和一些经验公式进
行计算。风荷载计算在基准风压基础上考虑了重现
期、结构体型、地形、地理条件等因素的影响。桥
梁动力失稳包括颤振失稳和驰振失稳,本章介绍失
稳机理及如何用运动方程和经验公式来判断桥梁是
否可能发生动力失稳。
8.1 相关的基本概念
桥梁抗风设计有很多地方不同于建筑结构,在此先将要涉及的
一些术语介绍如下:
基本风速:桥梁所在地区中的开阔平坦地面以上10m高度处100
年重现期的10min平均年最大风速
设计基准风速:在桥梁所在地区基本风速的基础上,考虑桥位
局部地表粗糙度影响的桥面高度处100年重现期的10min平均年
最大风速。
设计风荷载:进行静力抗风设计所采用的风荷载。跨度较小、
刚性较大的桥梁可只考虑阵风荷载作用下的强度问题,较大跨
度的柔性桥梁应考虑风致振动引起的动力风荷载作用。
风的攻角:由于地形的影响,近地风的方向可能对水平面
产生一定的倾斜度,称为风的攻角。具有攻角的风可能对
桥梁的风致振动,如颤振,产生不利的影响。一般认为高
风速时的平均攻角约在±3°之间。
阵风系数:瞬时风速与10min平均风速的比值。计算阵风荷
载时应采用时距为1~3s的瞬时(阵风)风速,即由阵风系
数乘以设计基准风速求得。
静力扭转发散:在空气静力扭转力矩作用下,当风速超过
某一临界值时,悬吊桥梁主梁扭转变形的附加攻角所产生
的空气力矩增量超过了结构抵抗力矩的增量,使主梁出现
一种不稳定的扭转发散现象。
静力横向屈曲:作用于悬吊桥梁主梁上的横向静风载超过
主梁侧向屈曲的临界荷载时出现的一种静力失稳现象。
颤振:是一种危险性自激发散振动,当其达到临界风速时,
振动的桥梁通过气流的反馈作用不断吸取能量从而使振幅逐
步增大直至最后使结构破坏。
驰振:对于非圆形的边长比在一定范围内的类似矩形断面的
钝体结构及构件,由于升力曲线的负斜率效应,微幅振动的
(8-2)
结构能够从风流中不断吸取能量,当达到临界风速时,结构
吸收的能量将克服结构阻尼所消耗的能量,形成一种发散的
横风向单自由度弯曲自激振动。
涡激共振:风流经各种断面形状(圆形、矩形、多边形等)
的钝体结构时都有可能发生旋涡的脱落,出现两侧交替变化
的涡激力。当旋涡脱落频率接近或等于结构的自振频率时,
将由此激发出结构的共振。
抖振:大气中的紊流成分所激起的强迫振动,也称为紊流风
响应。抖振是一种限幅振动,由于它发生频度高,可能会引
起结构的疲劳。过大的抖振振幅会引起人感不适,甚至危及
桥上高速行车的安全。
静力三分试验:采用主梁或桥塔的刚性节段模型,在风
洞中测定平均风绕流的静作用力的三个分量,即阻力、升
力和扭转力矩。无量纲的三分力系数和攻角的关系曲线反
映出断面的基本气动性能,是分析桥梁各种风致振动和静
力稳定的重要参数。
节段模型试验:将主梁的代表性做成刚性模型,用弹簧
悬挂在支架上形成一个有竖向平动、转动(及侧向)自由
度的振动模型,在风洞中测定风的动力作用。满足相似条
件的节段模型试验可直接测定二维颤振的临界风速,也可
识别出用气动导数表示的非定常动力,是桥梁最重要的风
洞试验之一。
全桥气动弹性模型试验:将全桥按一定几何缩尺制成并
满足各种必要的空气动力学相似条件的三维弹性模型,在
大型边界层风洞中观测其在均匀流及紊流风场中的各种风
致振动现象,用于考察桥梁从施工期各阶段到成桥的抗风
性能。是研究桥梁风致振动最精确的试验方法。
8.2 概述
桥梁一般是交通运输的咽喉,在国民经济中占有极重
要的地位,应能经受各种自然灾害而不轻易使交通中断。
风作为一种主要自然灾害,每年都给人民的生命财产带来
巨大损失,作为重要交通设施的桥梁也经常受到风的威胁
甚至危害。人们正是在各种桥梁风毁事故发生后开始重视
桥梁结构抗风设计,并逐渐形成桥梁抗风设计理论。
桥梁的风毁事故最早可以追溯到1818年,苏格兰的
Dryburgh Abbey桥首先因风的作用而遭到毁坏。之后,英
国的Tay桥因未考虑风的静力作用垮掉,造成75人死亡的
惨剧。一系列桥梁的风毁事故,使人们开始重视风的作用,
最初人们只认识到考虑静风载的必要性,直到1940年美国
Tacoma悬索桥的风毁事故(图8-1),才使工程界注意到
桥梁风致振动的重要性。
图8-1
Tacoma悬索桥的风毁事故资料照片
一、 风对桥梁作用的现象及作用机制
风对桥梁的作用是一个十分复杂的现象,它受到风的自
然特性、结构动力性能以及风与结构的相互作用三方面的制
约。由于地表的起伏和各种建筑物的影响,使得近地风的风
速和风向及其空间分布都是非定常的(即随时间变化的)和
随机的。当这种带有脉动成份的风绕过非流线形截面的桥梁
结构时,就会产生旋涡和流动分离,形成复杂的空气作用力。
这种作用力可能引起桥梁的振动,而桥梁结构的振动又将引
起流场的改变,这种相互作用的机制使得问题更加复杂。
从工程抗风设计角度,可以把自然风分解成不随时间变
化的平均风和随时间变化的脉动风两部分的叠加,分别考虑
它们对桥梁的作用,即静力作用和动力作用两种作用的现象
和机制见表1。
表1 风对桥梁作用的现象及作用机制
分类
现象
作用机制
静风载引起的内力和变形
平均风的静风压产生的
阻力、升力和力矩作用
扭转发散
静(扭转)力矩作用
横向屈曲
静阻力作用
静力
作用
静力不稳定
抖振(紊流风响应)
紊流风作用
限幅振动
涡振
动力
作用
驰振
自激振动
扭转
颤振
古典耦
合颤振
单自
由度
二自
由度
发
散
振
动
旋涡脱落引起的涡激
力作用
自激力的气动负阻尼
效应——阻尼驱动
自激力的气动刚度驱
动
二、 桥梁抗风设计目的和基本过程
(1)桥梁抗风设计的目的
桥梁抗风设计的目的首先在于保证结构在施工阶
段和建成后的营运阶段能够安全承受可能发生的最大
风荷载的静力作用和由于风致振动引起的动力作用。
因此,首先应掌握架桥地点的风特性,决定桥梁的设
计风速,并据此推算风对桥梁的作用,校核抗风安全
性,如果有可能出现有害的振动或变形,就应考虑适
当的防止对策或进行设计变更。
(2)桥梁抗风设计中的重要因素
抗风设计中的重要因素有:
(1)风特性参数 应通过调查和收集气象资料掌握桥
址处的风特性,并采用正确的方法确定合理的参数供抗
风设计使用。特别要注意桥址处特殊的地形、地貌和风
向条件,以便对常规的取值进行必要的修正。
(2)桥梁的动力特性 需采用合理的力学模型,并注
意边界支承条件的正确处理。对计算结果要通过与相似
桥梁的比较检验其合理性和可靠性,其中特别是对于主
梁前二阶对称和反对称的竖向弯曲、侧向弯曲和扭转振
型要作出正确的判断。
(3)桥梁风荷载、颤振临界风速、抖振响应 抖振响
应的正确预测主要取决于桥梁的动力特性、主梁断面的
气动特性和紊流风特性。
(3)桥梁抗风设计的基本过程
桥梁抗风设计的过程见图8-2。
对于一般的大桥,初步设计阶段的抗风分析可采用近
似的公式对各方案的静风载内力和气动稳定性进行估算,
待方案确定后再通过节段模型的风洞试验测定各种参数,
进行抗风验算和风振分析。对于重要桥梁,宜在初步设计
阶段通过风洞试验进行气动选型,为确定主梁断面提供依
据。在技术设计阶段再对选定的断面方案进行详细的抗风
验算和风振分析,还应通过全桥模型的风洞试验对分析结
果予以确认。
重现期
结构型式
气象资料
动力特性
基本风速
截面选择
假定Th值
设计风速
假定阻力系数
颤振风速估计
不安全
阵风系数
很安全
稳定性验算
及格
设计风载
很安全
不安全
静力抗风设计算
及格
节段模型风洞试验
阶段I
三分力试验
各类风振分析
满意
静力抗风验算
否
是否要进行
全模型验算
不够满意
抗风措施
阶段II
是
全桥气弹模型试验
满意
阶段III
检验
不够满意
抗风措施
图8-2 桥梁抗风设计过程
8.3 风对桥梁的静力作用
一、作用在桥梁结构上的平均风荷载
1. 主梁静力三分力及静力三分力系数
FV
FL
y
FD
MZ
α
a
风
o
体轴坐标系和风轴坐标系
x
FH
在体轴坐标系下,静力三分力表达为:
阻力
FH
1
 V 2 HLC H
2
1
V 2 BLC V
2
升力
FV 
力矩
1
M  V 2 B 2 LC M
2
1
2

V
式中:为空气密度,H为梁高,B为梁宽,L为长度,2
为气流的动压。CH、CV、CM分别为主梁的阻力系数、升力
系数、力矩系数,它们由节段模型试验提供。
在风轴坐标系下,静力三分力表达为:
阻力
升力
力矩
1
FD  V 2 HLCD
2
1
FL  V 2 BLCL
2
1
M  V 2 B 2 LCM
2
1
2
式中:为空气密度,H为梁高,B为梁宽,L为长度,2 V
为气流的动压。CD、CL、CM分别为主梁的阻力系数、升力
系数、力矩系数。
2. 桥塔、主缆及拉索上平均风荷载
1
FD  V 2DLC D
2
式中:D为桥墩、塔柱宽度或拉索外径,其余参数意义
同上。计算桥塔和拉索承受的风荷载时,按风剖面变
化考虑不同高度的风速。由于桥墩、塔柱、拉索截面
较为规则其阻力系数CD可按《公路桥梁抗风设计指南》
取值或通过模型实测
a
二、扭转发散
Ka
a
V
弹性轴
扭转发散问题的几何位置与参数
令扭转弹簧刚度为Ka ,其含义为梁段发生单位转角所需的气动
力矩。扭转角为a,平均风速为V,桥面宽为B,则单位长度的气
动力矩为 :
Ma 
1
V 2 B 2 C M a
2
式中:C M a  为绕扭转轴转动的气动力参数。
代回风速计算式:
三、横向屈曲
对于单跨简支的悬索桥,可采用以下公式估算横向
屈曲临界风速
Vlb  K lb f t B
B
r
 ( ) ( )
H
b
K lb 
C L' Bc
1.88C D 4.54 
CD H
3
m
B r 1

;b  ; 
2
b
2 b b
Im
ft
; 
m
fb
四、桥梁空气静力稳定性的非线性分析
考虑结构的几何非线性及静力三分力随攻角的变化,
采用非线性有限元方法进行分析。该方法可以将横向屈
曲和静力扭转发散一并考虑,是研究桥梁空气静力稳定
性的较为完善的方法。
8.4 桥梁动力特性及其计算分析
桥梁结构动力特性是桥梁抗风设计中的重要
资料,进行桥梁风致振动计算分析和桥梁模型风
洞试验,都须以结构动力特性为依据。这里首先
介绍如何采用有限元方法进行悬索桥和斜拉桥的
结构动力特性分析,然后介绍 采用近似公式计算
桥梁基频。
一、 频率、振型及结构阻尼
频率——单位时间内系统简谐振动的次数,常记为f,单位为Hz(次/
秒)。简谐振动的频率等于周期的倒数。圆频率w2f,单位为(周/
秒)。
振型——结构以某一频率做箭谐振动时,结构各点相对位移的关系。
阻尼——结构在做有阻尼自由振动时振幅衰减的程度
对数衰减律
1 x0
  ln
n xn
阻尼比与对数衰减率的关系

2
1 
2
,  2
结构的频率和振型可以通过结构动力特性分析获得,
结构阻尼与材料、结构形式等多种因素有关,无法通过
计算取得。桥梁抗风设计中结构的阻尼比可以取以下经
验值:
桥梁种类
阻尼比
阻尼比的统计范围
钢桥
0.005
0.5%~1.0%
结合梁桥
0.01
1.0%~1.5%
混凝土桥
0.02
2.0%~3.0%
二、采用有限元方法计算桥梁结构动力特性
空间杆单元
空间梁单元
进行结构动力特性分析常用的商业有限元软件
1. ANSYS
2. ALGOR
3. SAP2000
4. ADINA
5. NASTRAN
6. ABAQUS
7. DIANA
悬索桥结构动力特性示例
阶次
频率(Hz)
振型特点
1
0.0693
纵漂
2
0.1419
L-S-1
3
0.2163
V-S-1
4
0.2690
V-A-1
5
0.3954
V-S-2
6
0.4106
L-A-1
7
0.4806
V-A-2
8
0.5020
主塔横摆
9
0.5129
主塔横摆
10
0.5614
T-S-1
11
0.6437
V-S-3
12
0.6620
V-A-3
13
0.7039
V-S-4
14
0.7045
L-S-2
15
0.7480
边跨竖向
16
0.7634
T-A-1
斜拉桥结构动力特性计算示例
1
2
3
4
第一对称竖弯
1
2
3
4
第一反对称竖弯
1
2
3
4
第一对称扭转
1
2
3
4
第一对反称扭转
桥梁颤振分析中重要的模态
•
•
•
•
第一对称竖弯
第一对称扭转
第一反对称竖弯
第一反对称扭转
桥梁抖振分析中重要的模态
• 结构前20-30阶模态
三、 桥梁基频的近似计算
单跨简支悬索桥一阶反对称竖弯
 2 

L


wba  2f ba  
2
EJ  2 H g (
L 2
)
2
m
其中:L——跨度(m);
EJ——加劲梁的竖弯刚度(N.m2);
Hg——恒载单根主缆水平拉力(N);
m——桥面和主缆的单位长度总质量,m=md+2mc;
md—— 桥面单位长度质量(kg/m);
mc——
单根主缆单位长度质量(kg/m)。
(md  2mc ) gL2
Hg 
16 f
对于跨度超过500米的悬索桥
f ba 
1.11
f
f-——主缆矢高。
单跨简支悬索桥一阶对称竖弯
0.3 Ec Ac
f 
L
W
s
b
其中:L——跨度(m);
Ec——主缆弹性模量(N.m2);
Ac——单根主缆截面积(N);
W——桥梁单位长度重量;
单跨简支悬索桥一阶反对称扭转
fta 
1
L
H g Bc2
2 2
EJw ( )  (GJ d 
)
L
2
Bc2
2
md r 
mc
2
其中:EJw,EJd——加劲梁的约束扭转和自由扭转刚度(N.m2和N.m4),
对闭口箱梁,约束扭转刚度可以忽略;
r——加劲梁惯性径(m);
Bc——两边索的中心距(m);
W——桥梁单位长度重量;
单跨简支悬索桥一阶对称扭转
2
GJ

0
.
05256
E
A
B
1
d
c c c
fts 
2
2L
B
md r 2  c mc
2
8.5 桥梁空气动力失稳
一、非定常气动力
当结构发生振动时,由于周围的绕流受结构变位的干扰而发生变化,
导致作用在结构上的空气力也随时间发生变化。这种作用力由于是伴随
结构振动产生的,称为自激气动力,它是非定常气动力的主要形式。
(1)Theodorson平板空气力公式
1935年,Th. Theodorson首先
从理论上研究了薄平版的非定
常气动力。他根据流体力学势
流理论求得了作用于振动平板
上的非定常气动力的解析表达
式。对于图示二维理想平板,
在均匀水平流场中作微小振动
时所受到的非定常空气升力和
力矩可表达为:


h 
b a 
L  2bv C (k ) a    [1  C (k )]

v
2
v






h 
b a 
2 2
M  2b v C (k ) a    [1  C (k )]

v
2 v 


2
式中
 ——空气密度;
b——平板的半宽度;
v——空气流速;
h和a分别为截面的竖向和扭转位移;
w b
k
为折算频率,w为圆频率;
v
C (k )为Theodorson函数, 当用Bessel函数表示时可以写成
C k   F (k )  iG(k )
(2) Scanlan桥道气动力公式
由于桥梁断面属于非流线型的钝体,其绕流形态及空气力比
平板复杂得多,无法通过解析方法获得桥道的非定常气动力。
为此,Scanlan提出了用颤振导数(flutter derivatives, 也称为
气动导数)表达的线化的桥梁断面非定常气动力模型,通过桥
梁节段模型风洞试验测取表达式中的颤振导数。
h  2 hw h h  w h2 h  H 1 h  H 2a  H 3a
a  2 w a  w 2a  A h  A a  A a
a
a
a
1
2
3
式中:
H1  H
*
1
b 2w
m
3
* b w
H2  H2
m
3
2
* b w
H3  H3
m
A1  A
*
1
b 3w
I
4
* b w
A2  A2
I
4
2
* b w
A3  A3
I
*
*
无量纲系数 Hi , Ai (i  1,2,3) 称为颤振导数(气动导数)。
颤振导数的物理意义
1.
H1*,H 2*,A1*,A2*均与速度项有关,具有气动阻尼的意义;
H 3*,A3*均与变位有关,具有气动刚度的意义。
2.
H1*表示竖向运动对升力的贡献,H 2*,H 3*表示扭转运动对升力的贡献;
A1*表示竖向运动对力矩的贡献,A2*,A3*表示扭转运动对力矩的贡献。
典型断面颤振导数曲线
典型断面颤振导数曲线
二、桥梁颤振理论
经典颤振(弯扭耦合颤振)
——竖弯模态和扭转模态相互耦合的弯扭
耦合颤振,常发于扁平流线型桥梁断面。
分离流扭转颤振(单自由度扭转颤振,失速颤振)
——以扭转模态为主的颤振,常发生于钝
体桥梁断面,如槽型、工字型断面。
颤振计算分析——确定颤振临界风速Vcr
(1)桥梁颤振计算理论的发展
1948年Bleich 第一次用Theodorson 的平板空气公式来解决悬索桥的颤振分析。
他认为在悬索桥中常用的桁架加劲梁的上承桥面接近于一块平板,此时悬索桥的
二维颤振微分方程可以写成:
上述颤振计算方法只能考虑桥梁两个模态的贡献,因而是近似的计算方法。
随着计算机和有限元方法的发展,建立了能够考虑桥梁全部模态贡献的多种频域
分析方法和时域分析方法,如状态空间法、复模态分析方法、PK法等。这些方法属
于“精确”的分析方法。
(2)影响桥梁颤振的主要因素
1. 重要因素
• 扭弯频率比

ft
fh
越大,颤振临界风速越高。
• 质量及质量惯矩
m、Im越大,颤振临界风速越高。
• 结构阻尼
对于分离流扭转颤振,扭转阻尼比越大,颤振临界风速越高;
对于弯扭耦合颤振,结构阻尼对颤振的影响较小。
• 主梁气动外形
2. 振型贡献率
(3) 抗风设计中的经验公式
对于非平板的实际桥梁断面,一般应通过风洞试验确定相对于平板临界风速的截
面形状折减系数hs和攻角折减系数ha,在初步设计阶段可按Van der Put公式或Selberg
公式计算平板颤振临界风速,再 照下表给出的折减系数进行修正
关于振型的选取
在二维颤振计算中,正确的选取竖弯模态和扭转模态十分重要。通常桥梁的
弯扭耦合颤振总是以第一对称竖弯模态与第一对称扭转模态的组合、或者第一反
对称竖弯模态和第一反对称模态的组合为控制模态。
用于纯扭转颤振的经验公式
对于非流线型的桥梁断面,由于A2*曲线的负阻尼倾向,容易发生扭转形式的
颤振,颤振的频率接近扭转频率,由结构阻尼和气动负阻尼总和为零的条件为:
时的折算速度即为临界值Kc,此时相应的颤振临界风速为
Bfa
Kc
无量纲的Kc值称为Theodorson值(Th),于是:
vcr 
1
vc  Th Bfa
(4) 抗风设计中颤振稳定性的检验
1. 颤振稳定性检验准则: Vcr  [Vcr ]
式中颤振检验风速:
[Vcr ]  K f Vd
K为安全系数,取1.2;
f为考虑风速脉动影响及水平相关特性的无量纲修正系数;
Vd为设计基准风速;
2. 颤振稳定性分级
为了快速判断大跨度桥梁的抗风稳定性,可采用近似公式估算桥梁
的第一扭转频率,根据桥址处的气象资料估计出颤振的检验风速,然后
按照下式计算主梁断面所需要的最小Th-1值:
1
minTh 
[Vcr ]
ft B
然后根据minTh-1的大小,参照下表可以判断桥梁的抗风稳定性。