基于MATLAB语言的多自由度振动系统的固有频率和主振型计算分析.pdfVIP

78 2007年1月中国制造业信息化第36卷第1期 基于MATLAB语言的多自由度振动系统的 固有频率及主振型计算分析 文涛，胡青春 (华南理工大学机械工程学院，广东广州 510640) 摘要：多自由度振动系统固有频率及主振型计算分析是研究其振动特性的基础，矩阵迭代法是计 算固有频率及主振型的基本方法之一。根据矩阵迭代的方法，利用MATLAB编程并验证程序的 正确性。通过程序的运行，能快速获得多自由度振动系统的固有频率以及主振型，为设计人员提 供了防止系统共振的理论依据，也为初步分析各构件的振动情况以及解耦分析系统响应奠定了 基础。 关键词：凇TL蟠；多自由度；振动系统；固有频率；主振型 中图分类号：THll3 文献标识码：A 在工程振动中，确定系统固有频率与主振型是 采用矩阵迭代求解系统固有频率与主振型时，大部 非常重要的。固有频率是决定系统振动特性的重 分都是用Visiual 要物理量，它既是防止系统共振的依据，又是多自 序旧J。限于Visiual 由度系统解耦分析(模态分析)的前提，因此研究某 语法的局限性，用这种高级语言编写的程序涉及到 系统振动时，首先要求出系统的固有频率。主振型 选择合适的算法和编写冗长的语言代码以及键人 则为初步分析各构件的振动情况以及解耦分析奠 和调试等一系列问题。即使有现成的标准子程序 定了基础。 可供调用，要在一些较复杂的、科研问题中编写一 对于多自由度振动系统，计算系统固有频率与 个完整的程序仍然是一个复杂的、技巧性很强的工 主振型主要有2种方法…：(1)利用特征矩阵方程作。因此，用高级语言编写的程序一般代码段较 式与特征方程式求解；(2)矩阵迭代法求解。2种长，需要调用的子程序较多，整个程序的通读性较 方法各有各的特色。对于低自由度的振动系统，方 差。相反，MATLAB则有简洁、可读性强等优点。 法一容易、快捷。但是在实际工程中，大多数振动 系统都是自由度较多，用特征矩阵方程式与特征方 2 程式求解系统固有频率与主振型这种传统的计算 系统固有频率以及主振型分析 方法虽然从原则上可行，但当自由度增加时，惯性、 2．1 矩阵迭代求解系统固有频率以及主振 刚度阵的阶数增高，计算量也急剧加大，这显然很 型分析[4,5] 不方便。但采用矩阵迭代法，即使是自由度很大的 在多自由度正定系统的自由振动中，其固有频 振动系统，计算量也只不过是多进行矩阵迭代而 率及主振型可由迭代式 已，而且假设的初始矩阵愈接近实际状况，迭代的 次数愈少，相应的计算量也愈少。 ÷{A(，’}=[6][M][A(J’](1) ∞”．J 或 1 MATLAB语言的优点 (2) ∞i，，A‘J’=[M]一1[K][A‘J’] MATLAB作为一个以矩阵和数组为核心计算 求得。实践得知，用式(1)进行迭代将首先获得最 的软件，对矩阵迭代法中的矩阵迭代计算尤其适 低阶的固有频率及主振型，