桥梁气动阶跃函数及颤振时域算法中的数值问题研究.pdfVIP

第十四届全国结构风工程学术会议论文集 桥梁气动阶跃函数及颤振时域算法中的数值问题 张志田 胡海波 陈政清 (湖南大学风工程试验研究中心，410082) 摘要从航空学中机翼断面的阶跃升力函数着手，阐述了其用于桥梁断面气动力表述的演变过程及存在的问 题，选取了适合桥梁断面气动力时域表达的升力以及升力矩阶跃函数形式并讨论了其合理性．以正确模拟基于试验 的桥梁断面自激力线性表达式为目标，提出阶跃函数各控制参数非线性拟合的线性搜索方法并给出合理的参数控制 区间，有效地避免了过大的阶跃函数瞬态值与长时间衰减过程造成的自激力时域模拟失真．通过某大桥的三维颧振 时域分析，对数值算法中颤振临界风速的关键影响因素进行了讨论． 关键词桥梁、自激力、阶跃函数、时域、颤振 1．引 言 桥梁颤振稳定性分析有时域与频域两大类算法，在频域算法中，通常的做法是通过振型分解法 建立多自由度矩阵运动方程，然后将以运动状态表达的桥梁断面自激力表达式合并至结构运动方程 的左端，形成考虑气动刚度与气动阻尼矩阵的结构有效刚度矩阵与阻尼矩阵，最后通过参数搜索的 方法得出结构颤振临界风速。在时域算法中，关键的问题是获取桥梁断面的自激力时域表达式，目 前，自激力时域表达有两种方法，其一是采用脉冲响应函数(Impulse 函数的表达方法Ⅲ；其二是采用阶跃响应函数(Indicialresponse 表达的方法捌。对于钝体桥梁断面，由于分离流的存在，必须采用风洞试验的方法来确定白激力时 域表达式中的相关参数，目前有直接与间接两种方法进行参数识别，直接方法针对自激力的时域表 达式在风洞试验中直接测定，文献[3]与文献[4]分别尝试了上述两种自激力时域表达式中的参数测 激力表达式与时域表达式功率谱等效的原则进行时域参数识别晗一。目前使用较多的是间接方法。本 文采用阶跃函数进行时域自激力的表达并通过Scanlan自激力表达式进行参数识别，针对时域自激 力表达式的识别以及使用该法进行桥梁三维颤振分析过程中的数值问题进行深入地分析与讨论。 2．桥梁断面自激力 在航空学关于机翼断面气动力的早期研究中，将机翼与流体相互作用分成几种基本类型，其中 比较典型的有：1)机翼断面有效风攻角的阶跃变化引起的气动力问题；2)结构围绕其平衡位置小 幅振动引起的气动力问题；3)结构以某一速度穿越具有某种垂直风速特性的半无限空间引起的气动 力问题。前两种基本情况都与结构的运动状态有关，围绕其建立起来的理论为后来的飞机机翼颤振、 钝体桥梁的颤振分析奠定了基础；而第三种情况则是后来的采用Sears函数进行气动导纳研究【6】、结 构抖振分析的理论基础。 以升力为例，在均匀流场中，假设结构是固定不动的，则单位长度的流线型(无特征紊流产生) 结构断面所受的定常升力为： pu2Bq (1) 荆=三1 式中p为空气密度，u为来流风速，B为机翼弦长(或某种表征结构尺寸的参考宽度)，q为 637 第十四届全国结构风工程学术会议论文集 升力系数，为风攻角a的函数。 假设断面在某一时刻相对于来流的攻角姿态突然产生一个阶跃位移伽，则由此形成的气动升力 时程为： ， 1 ． 三(s)=去户u2Bq‰认s) (2) Z 纲时间变量。在机翼的瞬态升力理论中，阶跃升力函数常采用以下形式： 烈s)=1一ae一缸-ce一出 (3) 后由Jones在1940年的NACA报告中给出其解析表达式如下【8】： 研s)=1一O．165e-o·045如一0．355e-o。35