一种运用不完全模态试验数据的无溢出模型修正方法.docVIP

一种运用不完全模态试验数据的无溢出模型修正方法 第27卷第7期 振动与冲击 JOURNALOFVIBRATIONANDSHOCK 一 种运用不完全模态试验数据的无溢出模型修正方法 袁永新,蒋家尚 (江苏科技大学数理系,江苏镇江212003) 摘要:模型修正即为利用结构现场实测的振动信息修正不精确的有限元结构动力模型.鉴于测量的模态数据 一 般是不完整的,因此在进行模型修正时,希望测量的低阶模态数据融于修正模型而不改变原模型的高阶模态数据.这 样的一种修正,如果可能,称为无溢出模型修正.在假定有限元模型的质量矩阵与刚度矩阵均为对称非负定矩阵,并且质 量矩阵是精确的情况下,提出了一种无溢出有限元模型修正方法.该方法可使得修正模型仍为对称非负定,测试模态与 测试频率融于修正模型,并且修正模型的剩余模态和频率与原模型一致.最后数值例子验证了该方法的有效性. 关键词:有限元模型;无溢出;模型修正;模态数据 中图分类号:0327文献标识码:A 很多实际问题,如横梁,桥梁,建筑物等的振动问 题都可以用二阶微分方程来描述.直接求解这类连续 型的微分方程一般是比较困难的.因此常用的做法是 运用有限元技术将其离散化,从而得到如下的二阶常 系数微分方程: Moq(t)+q(t)=f(t)(1) 方程(1)称为系统的有限元模型,q(t)是凡维状态向 量,f(t)是凡维外激励向量,Mo,KoR分别称为分 析质量矩阵与分析刚度矩阵.通常M0,Kogt;--0,此 时式(1)称为对称非负定系统.与系统(1)对应的特 征值问题为 (AMo—K.)X=0(2) 使得关系式(2)成立的实数A和非零向量分别称为 这个系统的特征值与特征向量.若系统(1)的所有特 征值为A一,A,相应的特征向量为X一,X.记 A=digcA,…,A=()1)R×n, A1=dig(A1,…,A)∈R(3) X:[X一,X]=[l,2]R, X1=[X一,X]=R(4) 则式(2)可用矩阵形式表示为 MoXA—X=0(5) 另外,特征向量可以规范化使得下面的正交性 条件成立: 1[X]=J, 1[X]=diaglAA}(6): 基金项目:江苏省高校自然科学研究计划性项目(2oo5sL001J) 收稿日期:2007—08—06修改稿收到日期:2007—11—30 第一作者袁永新男,博士,副教授,1966年7月生 有限元方法由于具有适应性广,分析速度快,与结 构动力试验相比费用较低等优点,在实际工程中得到 了广泛的应用.然而,在多数情况下由有限元方法得 到的结果与实验测试得到的结果并不能很好地吻合. 导致这一现象的原因是有限元离散化建立的模型与实 际系统比较往往存在一定的误差.当这些误差较大 时,将导致由有限元方法所得到的结构模型动力特性 与实际测量的结果之间有较大的出入,甚至会超出工 程实际所允许的误差.近二十年来,随着建模技术和 动力测试技术的日益成熟,结构动力模型修正受到了 工程技术人员的广泛关注¨』.修正的结构动力模型 能够更加精确地预测结构的动力响应,还可以结合实 测结果对结构进行损伤检测和剩余寿命评估. 结构动力模型修正即为利用结构现场实测的振动 信息修正结构有限元模型,其中典型的工作是Bar— uch[] , Wei和Berman[m一12],戴华n等人提出的误差 矩阵范数极小化方法.鉴于测量的模态数据一般是不 完整的,因此在进行模型修正时,希望测量的低阶模态 数据能够融于修正模型而不改变原模型的高阶模态数 据.这样的一种修正,如果可能,称为无溢出模型修 正.本文在假定有限元模型的质量矩阵与刚度矩阵均 为对称非负定矩阵,并且质量矩阵是精确的情况下,提 出了一种无溢出有限元模型修正新方法.该方法可使 得修正模型仍为对称非负定,测试模态与测试频率融 于修正模型,并且修正模型的剩余模态和频率与原模 型一致. 具体地说,本文考虑如下问题: 问题P设测量的特征值组成的矩阵为=diag {l…}R,且y0,i=1,2,…,m,相应的 实测模态矩阵为Y=[Y…Y]R,求K0 使得 . y1=KYl,MaX2A2=KX2(7) 其中:,/I由式(3)与式(4)给出. 2振动与冲击2008年第27卷 通常的数值方法难以精确地求得特征值问题(2) 的高阶特征对,但可以较精确地求得特征值问题(2)的 低阶特征对.本文将避免用高阶特征值和相应的 特征向量表示问题P的解,而是借助于特征值问题 (2)的低阶特征对给出问题P的解的显式表达式. 1问题P的解 为了求解问题P,我们需要下面的引理. 引理1n设特征向量[]满足条件(6), 则对称矩阵满足矩阵方程MoA=KX当且仅当 为 K=K.+MaX1q,XlMa(8) 其中为任意的m阶对称矩阵. 下面我们将选取适当的=..使得式(8)定义的 矩阵满足式(7