振动桥梁CFD数值模拟的网格运动算法.pdfVIP

第 28卷 第4期 计 算 物 理 Vo1．28．No．4 201i年 7月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Ju1．．2O11 文章编号 ：1001—246X(2011)04-0547-07 振动桥梁 CFD数值模拟的网格运动算法 索奇峰 ， 吴 晴 ， 钟易成 ， 李明水 (1．西南交通大学土木工程学院，四川 成都 610031；2 南京航空航天大学能源与动力学院，江苏 南京 210016) 摘 要 ：针对非流线型桥梁断面模型，采用 Delaunay背景网格方法、棱边弹簧及线性弹簧 网格方法 ，研究非结 构网格与结构网格条件下，计算域 网格运动的特点及算法优劣．研究表明，在正常的桥梁断面振动振幅 内，三 种 网格运动方法均能保持网格变形时不会失效，其中Delaunay背景网格方法耗时最少 ，线性弹簧 网格方法获得 网格质量最高；在大振幅状态下，结构化网格容易失效，线性弹簧网格方法可最大限度保证网格有效． 关键词 ：网格运动 ；非流线型断面 ；强迫振动 ；计算流体力学 中图分类号 ：TB126 文献标识码 ：A 0 引言 计算流体力学(CFD)方法 已广泛应用于桥梁颤振及涡激振研究中¨I3]，即通过数值模拟桥梁断面周 围 的绕流流场，获得桥梁上的非定常气动力变化 ，进而对桥梁振动特性进行分析．考虑到桥梁大展弦比的特点 ， 一 般的桥梁绕流流场模拟均采用二维桥梁断面模型．基于非定常气动力 因素的桥梁振动特性研究主要包括 两部分内容 ：①模拟强迫振动条件下非定常流场，获得桥梁表面气动力随频率及振幅的变化规律 ，用于识别 颤振导数等 ；②与振动模型相耦合的受迫振动 ，通过引入结构控制方程 ，模拟一定结构条件下桥梁模型受气 动力影响产生 的限幅振动 ，从而分析桥梁的涡激振动特性．不论是强迫振动还是受迫振动，在流场数值模拟 的过程中，桥梁模型的位置和姿态均是随时间变化 的，因此在实现桥梁模型非定常流场数值模拟时，还须处 理随时间变化的网格化计算域 ，即实现运动网格算法． 基于有限体积方法的非定常流场数值模拟，须对计算域进行网格划分，并保证网格单元一一相邻，且不 发生交错 ．模拟桥梁振动时，网格壁面边界变化，若不相应改变内部 网格点位置 ，则变形后的计算域网格 极易出现交错现象，造成网格失效 ，导致桥梁非定常流动模拟失败．不仅如此 ，在计算域 网格变形时，保持较 为稳定的网格质量 ，对于计算的稳定性和准确性也极为重要．特别是对于粘性 网格区域 ，需要尽可能保持其 原有的正交性和长宽 比例．除此之外 ，随边界运动的网格变形 ，还需考虑计算耗时的长短 ，非定常问题计算步 数多 ，耗时长 ，应尽量缩短网格变形时的时间消耗． 本文采用基于弹簧原理的运动网格算法 和基于 Delaunay背景图映射 的运动网格算法 ，对典型的 非流线型桥梁断面非结构化计算网格和粘性结构化计算 网格 ，在桥梁模型振动过程中，对网格运动算法及变 形后的网格质量进行研究．对 比分析了两种运动网格算法下，不同类型网格变形时的质量优劣及时耗长短． 1 弹簧网格模型 1．1 棱边 (Segment)弹簧模型 弹簧模型的本质就是依据弹簧受力平衡构建合理的相邻 网格点位移关系 ，即通过在网格边上虚拟一个 弹簧 ，通过各弹簧受力平衡构建线性方程组 Aq=B， (1) 其中q=[qt q … qr，为各节点位移量 (边界节点位移量为已知条件，由主动位移条件获得，为网格 点个数)，B=[曰t … ] 为非齐次项，平衡态下的弹簧网格模型B=0，A矩阵(n×n)为位移系数 收稿 日期 ：2010—05—12；修回 日期：2010—10—12 作者简介：索奇峰(1974一)，男，河南洛阳，博士生，从事桥梁计算风工程研究．E-mail：qfsuo@126．com 548 计 算 物 理