基于模态应变能诊断结构破损的修正方法-东南大学学报.PDF

第３０卷第３期 东南大学学报（自然科 学版 ） Ｖｏｌ３０ Ｎｏ３ 　 ２０００年５月 ＪＯＵＲＮＡＬＯＦＳＯＵＴＨＥＡＳＴＵＮＩＶＥＲＳＩＴＹ（ＮａｔｕｒａｌＳｃｉｅｎｃｅＥｄｉｔｉｏｎ） 　 Ｍａｙ２０００ 　　　　　　 基于模态应变能诊断结构破损的修正方法 １ １ ２ 史治宇　　张令弥　　吕令毅 １ （南京航空航天大学飞行器系，南京２１００１６） ２ （东南大学土木工程学院，南京２１００９６） 摘　要　针对单元模态应变能法诊断结构破损时由于模态截断引起的计算误差等问 题，改进了求解系统特征向量导数的典型模态法，提出了仅用部分低阶模态来确定结 构破损位置和大小的修正方法，从而减小了模态截断对诊断结果的影响． 关键词　结构振动；模态分析；故障诊断；模态应变能 分类号　Ｏ３２７；ＴＵ３１１．３ 近年来，基于模态应变能来诊断结构破损的研究不断深入［１～７］，文［７］基于单元模态应变 能的变化导出了结构破损大小的确定方法，理论上需要完备的分析模态振型，而在实际计算中 由于仅用部分模态振型常常存在截断误差，更重要的是通常有限元模型能较好地吻合实际结 构的低阶模态，而高阶模态则有较大误差．基于以上分析，本文改进了文［７］的方法，导出了仅 用部分低阶模态来确定结构破损大小的方法，减小了模态截断所带来的计算误差． １　方法回顾［７］ 定义结构破损前后第ｊ个单元关于第ｉ阶模态的单元模态应变能变化（Ｃ）为 ｉｊ ｄ Ｔ ｄ Ｔ 　　Ｃ ＝（ ）Ｋ －（ ）Ｋ （１） φ φ φ φ ｉｊ ｉ ｊ ｉ ｉ ｊ ｉ 其中，Ｋ和 分别是结构第ｊ个单元的刚度矩阵和第ｉ阶模态振型；上标ｄ表示结构破损后的 φ ｊ ｉ ｄ 情况．由 ＝ ＋ ，展开式（１）并略去高阶项，则 φ φ Δφ ｉ ｉ ｉ Ｔ 　　Ｃ ＝２（ ）Ｋ （２） φ Δφ ｉｊ ｉ ｊ ｉ 式中， 是结构破损前后第ｉ阶模态振型的变化．由模态法，并考虑到结构破损通常是结构局 Δφ ｉ 部刚度的损失，而与质量无关， 可由下式求得： Δφ ｉ ｎ Ｔ －（ ） Ｋ ／（ － ）　　 ｋ ｉ φ Δφ λ λ ≠ 　　 ＝ ｄ ，　　ｄ ＝ ｋ ｉ ｋ ｉ （３） Δφ ∑ φ ｉ ｋ＝１ ｉｋ ｋ ｉｋ {０ ｋ＝ｉ 式中，ｎ和 分别是结构自由度总数和第ｉ阶的特征值；Ｋ是结构破损引起的刚度损失．通过 λ Δ ｉ