两自由度体系颤振机理研究.pdfVIP

第五届全国风工程及二r=业空气动力学学术会议论文集 两自由度体系颤振机理研究 项海吼张若雪顾明 同济大学土木工程防灾国家重点实验室上海200092 摘要；本文针对具有竖弯和扭转运动的两自由度体系的颤振问题，提出一种通过迭代求颤振方程 复模态解的方法，利用模态阻尼比小于零的判定条件，可以得到从风速为零到颤振发生整个过程中系统 的振动形态．通过对5F板和钝体的计算与分析，揭示了颤振发生的机理及各种断面颤振形态的内 在统一性．本文还考查了各气动导数对颤振风速的影响，证明了气动耦台项对经典耦合颤振发生的必 要性． 关键词：颤振气动导数复模态 引言 桥梁的颤振特性随着其断面几何形状及动力特性而改变，Scanlarl将颤振分为类似于 机翼的经典颤振和单自由度颤振两种形式…。单自由度颤振是由扭转方向的气动阻尼随风速 增大由正变负，最终导致结构一气流系统的总阻尼为负值而引起的分离流扭转颤振；而对于 经典耦合颤振，代表竖向和扭转两个方向气动阻尼的气动导数确+和42’一直保持为负，被 认为两个自由度都各自受到正的气动阻尼作用，因而用系统总阻尼的变化不能解释这～现 象。研究表明，这种颤振是由于弯扭振型之间的耦合而产生，Scanlan又将其称之为“刚 度驱动”而非“阻尼驱动”的颤振。如果桥梁断面的钝体性质比较显著，则可能发生单自 由度颤振，著名的Tacoma桥风毁就是因为主粱的扭转颤振引起。现代大跨桥梁一般都具有 较好的气动外形，多发生因气动耦合项参与而产生的经典耦台颤振。 单自由度颤振临界风速可通过系统机械阻尼与气动阻尼之和等于零来判断。而经典耦 合颤振一般采用Scanlan的半逆解法求解”。：假设系统的两个自由度方向均以单一的颤振频 率进行阻尼为零的稳态振动，根据气动导数和结构参数，对系统运动方程通过逐次迫近的 方法或者采用图解法求解得出颤振临界风速。这种半逆解法只能得到临界点的振动频率与 风速，而不能完撼地描述结构在从稳定到发散的过程中其振动特性的变化，从而不能揭示 颤振发生的机理。 本章提出一种迭代方法，利用体系的复模态解来描述从风速为零到颤振发生的整个过 程中体系的振动形态，这种方法针对耦合颤振问题和一般认为的单自由度颤振均可以求解。 通过对这两种情况下系统模态特性的分析，揭示颤振发生的机理及各种断面颤振形态的内 I一7 在统一性。 Scanl8n用气动导数表示的自激力公式在颤振分析中得到广泛的应用，本章的方法也是 以此为基础并研究了各气动导数对颤振临界风速的影响． 1计算模型 U 小M 图1 二维结构断面运动示意图 考虑在平滑来流作用下的一个翼型或桥面的断面。如图1所示，假设它具有两个自由 度：竖向位移^和扭转a。结构单位展长质量为m，质量惯性矩为，，质量中心与弹性中心 重合，则其运动方程可写为： Ⅲ(^+％^+簖^)=三，。 … I蝤+2‘。a+∞：a、=Mse 式中靠、靠分别为竖向和扭转运动的阻尼比，％和‰为自振频率，三一和肘一是物体运 动与气流相互作用产生的升力和力矩。Scanlan用六个颤振导数米描述自激升力和力矩： k=pv2口晒台+蹦警槲嘲 ㈤ 帆。=∥282蚴专+黝：鲁十尼2彳；a】 式中K=Beo／U是折算频率，气动导数H：和A：是K的无量纲函数，Ⅳ：，爿；，』j表示竖向 和扭转运动分别对升力和力矩的贡献，称为正导数或直接导数，Ⅳ；，^仁一}称为副导数或 交叉导数，表示一个自由度上的运动在另一个自由度方向上产生的力。理想薄平板的气动 导数理论解己由Theodorson根据势流理论导出．钝体的气动导数与断面具体形状有关。需 通过风洞试验获得。 2颤振体系的复模态解 将(2)代入(i)式并将方程的右端项都移到左边．可以写成如下的矩阵形式： 【^f】{i)+【C：1{i)+【K。】{x)={0}