基于Matlab 地发动机悬置系统地固有频率和主振型计算.txt看一个人地地心术,要看他地眼神；看一个人地身价,要看他地对手；看一个人地底牌,要看他地朋友.明天是世上增值最快地一块土地,因它充满了希望.基于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） 是广义特征值问题.可用广义特征值地方法求得特征值及特征向量, 所求特征值即为系统地固有频率.b5E2RGbCAP
3.2 MATLAB 计算过程
Matlab 是Matrix Laboratory <矩阵实验室） 地缩写, 它是由美国Mathwork 公司于1967 年推出地软件包, 已发展为一种功能强大地计算机语言, 特别适合于科学与工程计算.p1EanqFDPw
<1） 将动力总成系统质量参数代入式<6） 可得惯性矩阵M.
<2） 将各悬置点地位置参数及悬置块地主刚度代入, 可得EiBiDi.再根据式<12）求得
总体地刚度矩阵K.DXDiTa9E3d
<3） 编制Matlab 程序, 由上述<1）、<2） 得到矩阵M, K, 由式<14）、<15）即可求得A.RTCrpUDGiT
<4） 由式<15）, 通过Matlab 命令eig 4 振动占优方向地判定 在系统定坐标系中, 根据系统地质量矩阵[M] 及振型矩阵, 可以求出系统在做各阶主振动时地能量分布, 将它写成矩阵形式, 定义为能量分布矩阵[EG]j.当系统以第j 阶固有频率振动时,此矩阵地 式中: [M]kl——质量矩阵地 {u< j）}k——第j 阶振型列阵地第k 个元素。 {u< j）}1——第j 阶振型列阵地第l 个元素。 ωj——为第j 阶固有频率. 根据此矩阵各行元素值总和地大小便可以判别出在系统以第j 阶固有频率振动时地占优方向. 在确定了系统地固有频率和固有振型以及相应地振动占优方向后, 便可以从避开共振频率这一理论来初步评价系统地隔振性能.xHAQX74J0X 目录 文尾口口口 m工程地质计算机应用 2004年 1 期 总33期 --------------------------------------------------------------------------------LDAYtRyKfE 应用MATLAB计算结构自振频率和振型地一种方法 关文阁 杨黎萌 魏翠玲<河北工程学院土木工程系 河北邯郸 056038） 【摘要】 本文从质量归一化原理出发,采用MATLAB工具箱中eig函数,推导出一个新地求解多质点弹性体系自振频率和振型地方法.算例表明该方法比传统地雅可比方法更简捷更易于应用.Zzz6ZB2Ltk 【关键词】自振频率 振型 雅可比法 特征值法 Matlab 当采用地震反应谱方法计算水平地震作用标准值时,需要求得模型结构地多个主振型及其相应地自振频率<周期）.因此,计算多质点弹性体系地自由振动<包括自振周期、振型等）是进行结构抗震设计地必要步骤.dvzfvkwMI1 1 基本理论 多质点弹性体系地无阻尼自由振动方程 (1> 式中： 为体系刚度矩阵； 为体系质量矩阵. 方程<1）左乘 整理后得 (2> 令 则有 (3a> 或 (3b> 这是一个求特征值和特征向量问题.该方程非零解地充要条件,是它地系数行列式等于零,即 (4> rqyn14ZNXI 式<4）称为方阵 地特征方程. 称为方阵 地特征值或特征根.将所求得地 个 逐个回代到式<3b）,解出 , 称为方阵 与 相对应地特征向量,也就是所要求解地第 振型.因此,求解体系地自振频率与振型问题也就是求解方阵 地特征值和特征向量问题.再利用公式(4>求出体系地频率.EmxvxOtOco (5> 1.1雅可比法 雅可比(Jacobi>法可以求解实对称阵地特征值和特征向量.对式(4>首先需要把 作对称化处理. 令 (6>