斜风作用下大跨度桥梁抖振响应分析时域方法

摘要

斜风作用下大跨度桥梁抖振响应时域分析方法以现有基于平均风分解法的斜风作用下大跨度桥梁抖振响应分析方法为基础，结合SHMS实测风特性数据的特点对其进行了改进，建立了适于SHMS采集数据的风分解法。同时为了实现气动自激力的时域化，创建了专用于模拟气动刚度和气动阻尼的Aero-dyn18单元。在上述基础上引入现有正风下大跨度桥梁抖振响应的时域分析方法，发展了一套能够全面考虑各种复杂因素的斜风作用下大跨度桥梁抖振响应的时域分析实用方法，编制了全部相关程序，从而实现了直接由SHMS实测风环境数据得到结构的抖振响应，便于广大桥梁工程技术人员推广应用。

主权项

一种斜风作用下大跨度桥梁抖振响应时域分析方法，其特征在于该方法包括以下步骤：第一步：基于结构设计图纸和现场动静力响应实测数据，采用可编程参数化设计语言建立大跨度桥梁有限元计算模型，同时将主梁断面的颤振导数以可调用数组方式进行存储，第二步：调用实测风环境数据之后，采用斜风分解法对其进行处理，并基于分解后得到的风速数据进行大跨度桥梁桥址区三维脉动风场模拟，第三步：根据风速数据和颤振导数图确定式“4a、4b、5a、5b”中主梁断面的颤振导数： K ae e = K ae 1 e 0 0 K ae 1 e - - - 4 a 其中， K ae 1 e = ρ U 2 K 2 L e 2 0 0 0 0 0 0 0 P 6 * P 4 * B P 3 * 0 0 0 H 6 * H 4 * B H 3 * 0 0 0 B A 6 * B A 4 * B 2 A 3 * 0 0 0 0 0 0 0 0 0 0 0 0 0 0 - - - 4 b C ae e = C ae 1 e 0 0 C ae 1 e - - - 5 a 其中， C ae 1 e = ρUBK L e 2 0 0 0 0 0 0 0 P 5 * P 1 * BP 2 * 0 0 0 H 5 * H 1 * BH 2 * 0 0 0 BA 5 * BA 1 * B 2 A 2 * 0 0 0 0 0 0 0 0 0 0 0 0 0 0 - - - 5 b 式“4a、4b、5a、5b”中，Kaee和Caee分别为当采用集中气动力矩阵模拟桥面受到的气动力时，与桥面单元e所受气动自激力等价的气动刚度矩阵和气动阻尼矩阵；ρ为空气密度；U为平均风速；B为桥面宽度；K＝Bω/U为无量纲频率；Le为单元e的长度；Hi*为自激升力相关气动导数、Pi*为自激阻力相关气动导数、Ai*为自激升力矩相关气动导数，其中i＝1，2，…，6；第四步：采用专用于气动自激力模拟的Aero-dyn18单元模拟Kaee和Caee，并在有限元计算模型的基础上添加Aero-dyn18单元，得到大跨度桥梁用于抖振分析的有限元模型，第五步：由分解后得到的风速数据，采用结构静力三分力系数计算出静力风荷载，包括阻力、升力和升力矩，再利用第二步三维脉动风场模拟所得风速时程按照式“6a、6b、6c”计算出抖振力时程： L b ( t ) = 1 2 ρ U 2 B [ 2 C L ( α 0 ) u ( t ) U + ( C L ′ ( α 0 ) + C D ( α 0 ) ) w ( t ) U ] - - - 6 a D b ( t ) = 1 2 ρ U 2 B [ 2 C D ( α 0 ) u ( t ) U + C D ′ ( α 0 ) w ( t ) U ] - - - 6 b M b ( t ) = 1 2 ρ U 2 B 2 [ 2 C M ( α 0 ) u ( t ) U + C M ′ ( α 0 ) w ( t ) U ] - - - 6 c 式“6a、6b、6c”中，Lb(t)、Db(t)和Mb(t)分别表示升力时程、阻力时程和升力矩时程，ρ为空气密度，U为平均风速，B为桥面宽度，α0为平均风攻角，CL、CD、CM分别为升力、阻力和升力矩系数，CL′、CD′、CM′分别为升力、阻力和升力矩系数曲线斜率，u(t)、w(t)分别代表顺风向和竖向脉动风速；第六步：将上述所得静力风荷载以及抖振力时程同时施加到有限元模型上，采用纽马科数值积分法进行大跨度桥梁抖振响应非线性时程分析求解；第七步：分析处理大跨度桥梁的加速度、位移时程响应计算结果，求出加速度均方根响应值，并与现场实测结果进行对比分析；第二步所述的斜风分解法按照式“3a、3b”进行风速分解：UN＝Vcosβ                        3aUP＝Vsinβ                        3b式“3a、3b”中，UN和UP分别表示分解之后横桥向和顺桥向的瞬时风速；V和β分别表示风速仪实测瞬时风速大小和方位角；第四步所述的Aero-dyn18单元中具有两个节点，每个节点有6个自由度；Aero-dyn18单元的单元矩阵为12×12矩阵，其中包括C1，C2，…，C9总计9个不同的系数： 0 0 0 0 0 0 0 0 0 0 0 0 0 C 1 C 2 C 3 0 0 0 0 0 0 0 0 0 C 4 C 5 C 6 0 0 0 0 0 0 0 0 0 C 7 C 8 C 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 C 1 C 2 C 3 0 0 0 0 0 0 0 0 0 C 4 C 5 C 6 0 0 0 0 0 0 0 0 0 C 7 C 8 C 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 该单元没有固定的几何形状，它通过实常数的方式输入气动刚度矩阵和气动阻尼矩阵以模拟结构系统的气动刚度矩阵和气动阻尼矩阵。