大跨度桥梁非线性颤振实验的几点研究

一、脉动风速的模拟试验

要进行抖振时程分析就必须首先模拟作用在桥梁上的脉动风速时程。本文采用经改进的谐波合成系列中的一种方法，大大提高了模拟效率，为在后文进行颤振时程分析中不断变换风速计算节约了时间。作用在大跨度桥梁上脉动风速可视为一维多变量随机过程。众所周知，用谐波合成法模拟一维多变量随机过程需要计算互谱密度矩阵的Cholesky分解。该分解通常采用迭代法求借，计算最大，常常影响模拟的规模的效率。本文作者利用桥梁上各点的互谱密度近似相等的特点，导出了显式的分解公式，并且采用了FFT技术，从而极大地提高了模拟效率。

二、风荷载计算

引起桥梁风振的荷载可以分为静力风荷载，抖振力和自激力。其中静力荷载按常规静力三分力系数计算，抖振力常按Scanlan的准定常理论计算。

自激力的计算一直是研究得较多的课题之一。传统频域抖振和颤振分析方法中的自激力都采用Scanlan提出的气动导数的线性表达式。由于该表达式是频域和时域的混合表达式，不能在时域中求解。为了在时域中顺利计算耦合自激力，Lin提出了一种用单位脉冲响应函数表达的统一自激力表达式。本文按Lin的理论计算耦合自激力。Lin的理论基于二自由度耦台。然而，三自由度耦合对结构振动的影响最近也引起了一些学者的关注。虽然并非所有的自由度之间都具有耦合特性，但从理论和形式完备的角度出发，本文将Lin的理论从二自由度推广到三自由度，成功地实现了时域内三自由度耦合自激力的计算。

用脉冲响应函数表达的自激力适合于任意形式的振动，也适用于正余弦振动(颤振)。根据在正余弦振动形式下，脉冲响应函数表达的自激力与气动导数表达的自激力相等价的关系，Lin导出了用脉冲响应函数表达的自激力的具体表达形式。

三、统一的额报和抖报时域分析方法

在传统的步域分析方法中，抖振和颜振是通过完全不同的方法来分析的。其中，抖振分析用的是基于随机振动理论的响应谱方法，颤振分析用的是与特征值问题有关的半逆解法或复模态解法。风振时程分析的初衷是为了解决非线性情况下的抖振响应计算。但是颤振分析中所需要的计算自激力的公式在抖振时程分析中都要用到，所以从理论上讲，利用计算抖振时程分析的方法同样可以在时域中计算颤振。实际上，抖振和颤振并不是完垒独立的。在任何风速之下，桥梁都受到抖振力和自激力的作用。当风速较低时，自激力很小，不起控制作用，桥粱的振动就体现为抖振。当风速增加到一定程度时，自激力逐渐发散，并控制桥梁的运动，桥梁就发生了颤振。因此，只要正确地描述了抖振力和自激力，运用时程分析这一仿真的分析方法，就可以算出一定风速之下桥梁的真实运动状态。如果表现为随机振动，则说明是抖振，我们就可以得到响应时程统计指标。如果是发散振动，就说明桥梁发生了颤振。只要不断进行搜索计算，我们就能在时域中找到桥梁的颤振临界风速。

根据以上设想，本文设计并首次成功地实现了时域中统一的颤振和抖振分析算法。

流程中，耦合自激力的计算是个关键。过去的一些抖振时程分析方法中常只近似考虑非耦合的自激力。而大跨度桥梁的颤振发散大多是受耦合自激力控制的，因此，过去的抖振时程分析方法不能用于计算颤振的原因就在于此。颤振发散的判断依据也是关键之一。考虑到结构在接近颤振临界状态时，振动形式逐渐从随机振动过渡到谐波发散振动，其振幅将逐渐增大，相应振动的阻尼将逐渐减小。因此，本文先通过位移时程曲线观察振幅的变化规律，当结构的振动明显过渡为谐波振动时，则恨据计算结构的阻尼系统，当阻尼系统为负时，则认为结构进入颤振临界状态。

四、非线性颤振和抖振时程分析的程序设计

除了在时域中统一颤振和抖振分析方法以外，本文研究时程分析方法的目的还在于分析不同非线性因素对桥梁颤振和抖振响应的影响。与大跨度桥梁抖振和颤振有关的非线性现象主要有：

(1)几何非线性，包括平均风荷载引起的位移：由于大跨度桥梁相对细长，几何非线性现象不能忽视；

(2)有效攻角效应；由平均风荷载引起的位移使风对桥梁的攻角发生变化，从而使静力三分力系统和气动导数发生变化，因此附加攻角对桥梁的影响不能忽视。

根据以上分析流程并考虑这些非线性因素，借鉴一些通用有限元程序的理论和源代码。

大跨度桥梁在非线性情况下的颤振和抖振分析是目前桥梁风工程研究的热点之一。本文着重提出了时域中统一的颤振和抖振方法，同时解决了脉动风速的高效率模拟、结构几何非线性和气动非线性的处理方法。在此基础上，本文编制了计算程序Nbuffete-用该程序分析了江阴长江大桥非线陛颤振和抖振响应。结果表明本文提出的方法及所编制的程序在理论和实践上都是正确的。

在此基础上，我们就可以在时域中增加考虑各种非线性因素对结构进行分析从而寻找结构对这些因素的敏感性，我们也可以根据时程计算来进行非线性的振动控制。而这些研究工作在频域范围内是难以开展的。如果与CFD技术相结合，将可望实现从参数识别到结构宏观计算和控制的全过程分析。从而达到与风洞试验互为补充的目的。

应该说，尽管以上方法和程序是成功的，但是更重要的是要利用这种方法对所关心的桥梁进行备种用过去的方法所不能进行的全过程参数分析，从而得到更具有普遍规律性的结论。因此，大量的实例计算和总结是必要的。