大跨桥梁非线性静风稳定优化迭代计算方法论文.pdfVIP

第十四届全国结构风工程学术会议论文集 大跨桥梁非线性静风稳定优化迭代计算方法 张文明 葛耀君 (同济大学土木工程防灾国家重点实验室，上海200092) 摘要综合考虑静风荷载非线性和几何非线性，引用大跨度桥梁非线性静风稳定分析理论，通过引入外层迭 代次数上限对增量与内外两重迭代相结合的方法进行改进，提出了优化迭代分析方法，编制相应计算程序对世界最 大跨径钢箱梁悬索桥西堠门大桥的静风稳定性进行全过程分析，研究了该桥的失稳形态和机理． 关键词桥梁静风稳定非线性优化迭代 随着交通事业的迅猛发展，桥梁结构日益大跨化、轻柔化，对风致响应变得更加敏感。近年来 风洞试验和研究结果表明，大跨缆索承重桥梁存在静风失稳的可能性【1~3】。所谓静风失稳是指结在给 定风速作用下，主梁发生弯曲和扭转变形，一方面改变了结构刚度，另一方面因结构姿态的变化改 变了风荷载的大小，并反过来增大结构的变形，最终导致结构失稳的现象‘41。当结构变形引起的抗 力增量小于外荷载增量时，就会发生静风失稳【51。静风失稳是静风荷载与结构变形耦合作用的体现。 荷载增量与内外两霞迭代相结合的方法【6】是分析静风稳定性的有效方法，所谓荷载增量就是指 将风速按一定比例增加，内层迭代主要是进行结构的非线性计算，而外层迭代则是为了寻找结构在 某一风速下的平衡位置。在该方法中，如果风速增量很大，在某一较高风速下外层迭代要进行多次 才能收敛，耗费大量计算时间，而且可能出现在高于临界风速的某一风速下内外层迭代均收敛的情 况，因此不能准确地搜索到临界风速；如果风速增量很小，虽能保证准确地搜索到临界风速，但仍 会耗费大量计算时间。通过引入外层迭代次数上限，本文提出了优化迭代分析方法克服了上述不足， 随着风速的增加，外层迭代次数越来越多，如果外层迭代次数达到了上限，则减小风速增量，这样 在较低风速时可采用较大的风速增量，在较高风速时自动采用较小的风速增量，能准确地搜索到临 界风速，而且节省了计算时间。本文编制了相应计算程序对西堠门大桥的静风稳定性进行了全过程 分析，研究了该桥的静风失稳形态和机理，并验证了优化迭代分析方法的有效性和先进性。 1．静风失稳优化迭代分析方法及其实现 作用在主梁单位长度的静风荷载可分解为横向风荷载尸小竖向风荷载尸矿和扭转力矩Rf，具 体表达式【6】如下： PH=pViCH(tz)h／2 (1) 昂：pV2Cv(oOb／2 匕=pV2％(功62／2 式中：cH(a)、cv(a)、c从a)分别表示在有效攻角下主梁的阻力、升力、升力矩系数，所谓有效攻角 是指静风初始攻角与静风作用引起的主梁扭转角之和；h、b分别表示主梁的侧向投影高度和宽度； P为空气密度；矿为风速。 按照杆系结构空间计算理论，问题可归结为求解如下形式的非线性方程： [E(“)+砰御(甜)Ⅳ=，【乃(叻，ev(a)，％(功】 (2) 式中：恐和乓分别为结构的线弹性和几何刚度矩阵；a为有效攻角；％、Pv、砌分别为风阻力、 升力和升力矩；上标G和形分别代表重力和风力；U为结构位移。 708 从式(2)可知，不仅结构的刚度是结构变形的函数，而且右端项所表示的静风荷载也是结构变形 的函数，为了求解该非线性方程，就必须采用迭代法。本文采用Newton．Rapson法进行求解，因为 它有较快的收敛速度。而为了跟踪结构变形的全过程，又必须采用增量法。为此，本文采用增量与 内外两重迭代相结合的方法f6】，并引入外层迭代次数上限，具体实施步骤如下： (1)在自重作用下非线性求解(荷载步1)。 (2)提取主梁扭转角(为左右两节点扭转位移的平均值)，计算该状态下的三分力系数：此时，主 梁的有效攻角均等于初始攻角嘞。 (3)设定初始风速％和风速步长AV,当前风速％=v0；设定外层迭代次数上限卜『m越。 解(内层迭代)。 (5)提取单元扭转角(为左右两节点扭转位移的平均值)，计算该状