几种气动导数识别方法比较.pdfVIP

第十三届全国结构风工程学术会议论文集 几种气动导数识别方法的比较 姜保宋 周志勇 (同济大学土木工程防灾国家重点实验室，上海200092) 摘要在由sc．,a删lan建立起来的基于非变形结构和线性气动力模型的桥梁近似风振理论中，气动导数扮演着一 个较为重要的角色．本文对几种气动导数的时城识别方法从理论上和实例上进行了分析比较．对结果进行进一步的 分析并讨论指出算法的优势，不足． 关键词 气动导数、模态参数识别，桥梁结构 ，，宝l 暑 现代斜拉桥朝着超大跨度、极度柔性方向发展，桥梁跨径的加长使得结构对风的敏感程度大为 增强，避免结构发生发散性的颤振运动也因此成为大跨径桥梁设计过程中必须关注的问题。在桥梁 结构颤振和抖振分析中，气动导数扮演着非常重要的角色。 目前，对于桥梁气动导数识别的方法可以分为时域方法和频域方法，时域方法主要有最小二乘 法m娩删、随机子空间法慨sI)啪和特征系统实现方法锄呦协m1，本文主要对几种时域的识别方法从 理论上和实例上作一些分析比较。对结果进行进一步的分析讨论并进一步指出算法的优势和不足。 凳，算藩舞埂■端 雾，I，幡纂籍榜唇募曩朔昀蕾季万基 根据Scan]an的颤振分析理论，两自由度桥梁节段模型在均匀流场中的运动可以用如下的微分方 程描述： 朋语+旄蛾|I；+砰∞=L (1) “戊+2己吼d+畦叻=帆 (2) 其中k，虬分别为升力和力矩，可以写为如下形式： kt丢pU2(2口)【腰：僻吾+艇(幻可Bd+膏2巧(聊+置1E㈣吾 纲为一的系数，称为气动导数，它与桥梁断面的具体形状、结构自身的物理参数以及流体自身特性 有关，可通过风洞试验获得。令： 日。：坐Ⅳ：(量)日：；旦塑哦(七)马：坐E(七)矾：丛Z(的 ^=竽彳∞妒坐I咖)4=华4∞4=孚㈣ 则方程(1)，(2)变为： (3) (4) 菇+2{—％矗+霹矗=日．I；+月犯+脚+月r4^ 成+2蓦。nk戊+《口=4矗+彳2应+坞口+^^ 与方程(3)，(4)对应的状态方程为： 714 第十三届全国结构风工程学术会议论文集 矿=AY (5) 其中’』*三]c·=产Xq磕-一H2】肛瞄“-Ha】 几种时域识别方法就是通过求出不同风速下的特征矩阵A来实现气动导数的识别的。 雾，雾，■小昙墨港疗鞠 求解特征矩阵A的特征值和特征向量，对于欠阻尼系统，特征值和特征向量成共轭对出现，特 征值记为^，如，．．．九，^。，九‘，．．．九‘，对应的特征向量记为％，y：，．．．％，y。’。IP2*,．--V．‘，^和％分别称为系统 的复频率和复振型。体系的自由振动响应可以表示为： 撕)=窆胁‘+,fv／，·、) (6) 式中：c．，c|‘为与初始条件相关的常数。 如果我们能够得到矩阵A的特征值与特征向量，那么矩阵A便可以求得，从而可以得到结构的气 动导数。最d--乘法就是基于这一思想，事先假定初值，通过迭代的方法最终求得矩阵A的特征值 与特征向量，从而实现气动导数的识别。但是最小二乘法对初值假定的要求较高，初值不合理可能 造成迭代的不收敛，因此，一般先用ITD或者MITD的方法算得的结果作为初值。 雾．鲁，抟毽萋箍塞飘港难B^i纾辅 对前面得到的状态方程(5)，我们可以进一步转化为离散的线性时不变系统，可用状态方程表示 为： x(七+1)=ax(k)+肌(七) (7) 灭七)=c《妁 k=l，2…