基于Matlab的发动机悬置系统的固有频率和主振型计算76662 .pdf

仅供个人参考 基于 Matlab 的发动机悬置系统的固有频率和主振型计算（二） 3 运用 MATLAB对动力总成悬置系统固有特性的计算 3.1 理论计算 动力总成系统固有特性的计算 , 即计算系统的固有频率和振型。动力总成悬置系统无阻 尼的自由振动微分方程 : 式中 : M ——对称正定惯性矩阵 ; K ——对称正定刚度矩阵。 求多自由度振动系统的固有频率 , 从数学上讲就是求特征值的问题 : 设式（ 13） 的解为 : X=Xsin （ωt+a ） 代入式（ 13） 化简后得 : KX= ω2MX 左乘 M- 1 得 : M- 1KX= ω2X （14） 令 M- 1K=A, 则 : AX= ω2X （15） ω2 即为 A 阵的特征值 , X 为其特征向量。 由于 M 对称正定 , K 也是对称阵 , 因而式（ 13） 是广义特征值问题。可用广义特征值 的方法求得特征值及特征向量 , 所求特征值即为系统的固有频率。 3.2 MATLAB 计算过程 Matlab 是 Matrix Laboratory （矩阵实验室） 的缩写 , 它是由美国 Mathwork 公司于 1967 年推出的软件包 , 已发展为一种功能强大的计算机语言 , 特别适合于科学与工程计算。 （1） 将动力总成系统质量参数代入式（ 6 ） 可得惯性矩阵 M。 （2 ） 将各悬置点的位置参数及悬置块的主刚度代入 , 可得 EiBiDi 。再根据式（ 12）求 得总体的刚度矩阵 K。 （3 ） 编制 Matlab 程序 , 由上述（ 1）、（2 ） 得到矩阵 M, K, 由式（ 14）、（15）即可求 得 A。 不得用于商业用途 仅供个人参考 （4 ） 由式（ 15） , 通过 Matlab 命令 eig （A）, 即可求出矩阵 A 的特征值 ω2。利用 公式 ω2=2 πf, 即可得到悬置系统的各个振动固有频率 f 。 4 振动占优方向的判定 在系统定坐标系中 , 根据系统的质量矩阵 [M] 及振型矩阵 , 可以求出系统在做各阶主振 动时的能量分布 , 将它写成矩阵形式 , 定义为能量分布矩阵 [EG]j 。当系统以第 j 阶固有频 率振动时 , 此矩阵的（ k, j ） 元素为 : 式中 : [M]kl ——质量矩阵的（ k, j ） 元素 ; {u （ j ）}k ——第 j 阶振型列阵的第 k 个元素 ; {u （ j ）}1 ——第 j 阶振型列阵的第 l 个元素 ; ωj ——为第 j 阶固有频率。 根据此矩阵各行元素值总和的大小便可以判别出在系统以第 j 阶固有频率振动时的占 优方向。 在确定了系统的固有频率和固有振型以及相应的振动占优方向后 , 便可以从避开共振频 率这一理论来初步评价系统的隔振性能。 目录 文尾口口口 m 工程地质计算机应用 2004 年 1 期 总 33 期 应用 MATLAB计算结构自振频率和振型的一种方法 关文阁 杨黎萌 魏翠玲（河北工程学院土木工程系 河北邯郸 056038 ） 【摘要】 本文从质量归一化原理出发，采用 MATLAB工具箱中 eig 函数，推导出一个新的求 解多质点弹性体系自振频率和振型的方法。算例表明该方法比传统的雅可比方法更简捷更易 于应用。 【关键词】自振频率 振型 雅可比法 特征值法 Matlab 不得用于商业用途 仅供个人参考 当采用地震反应谱方法计算水平地